help button home button Biophys. J.
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
Originally published as Biophys J. BioFAST on August 11, 2006.
doi:10.1529/biophysj.106.088633
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow All Versions of this Article:
biophysj.106.088633v1
91/9/3206    most recent
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Feenstra, K. A.
Right arrow Articles by Vermeulen, N. P. E.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Feenstra, K. A.
Right arrow Articles by Vermeulen, N. P. E.
Biophysical Journal 91:3206-3216 (2006)
© 2006 The Biophysical Society

Enantioselective Substrate Binding in a Monooxygenase Protein Model by Molecular Dynamics and Docking

K. Anton Feenstra *, Karin Hofstetter {ddagger}, Rolien Bosch *, Andreas Schmid {dagger} {ddagger}, Jan N. M. Commandeur * and Nico P. E. Vermeulen *

* Division of Molecular Toxicology, Department of Pharmacochemistry, Leiden/Amsterdam Center for Drug Research (LACDR), Vrije Universiteit, Amsterdam; {dagger} ISAS-Institute for Analytical Sciences, Dortmund, Germany; and {ddagger} Department of Biochemical and Chemical Engineering, University of Dortmund, Dortmund, Germany

Correspondence: Address reprint requests to N. P. E. Vermeulen, Division of Molecular Toxicology, Dept. of Pharmacochemistry, Leiden/Amsterdam Center for Drug Research (LACDR), Vrije Universiteit, De Boelelaan 1083, 1081HV Amsterdam, The Netherlands. Tel.: 31-20-5987590; Fax: -31-20-5987610; E-mail: NPE.Vermeulen{at}few.vu.nl.


   ABSTRACT
TOP
 ABSTRACT
INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 CONCLUSIONS
 ACKNOWLEDGEMENTS
 REFERENCES
 
The two-component flavoenzyme styrene monooxygenase (SMO) is an efficient alternative to several chemical epoxidation catalysts on a preparative scale. A first homology model of the catalytic domain (StyA) of SMO was constructed (Protein Data Bank ID 2HD8) based on the structure of para-hydroxybenzoate hydroxylase. The StyA protein structure was optimized by restrained molecular dynamics to reproduce specific pre-S binding orientations of styrene. Effects of all 10 point mutations examined were explained by the distance of the site to the styrene and FAD binding sites. Thirteen of 20 ligands could be accommodated in a catalytically active binding orientation, and predicted affinities correlated well with experimental turnover and inhibition. The binding cavity is almost completely hydrophobic except for a hydrogen-bonded network formed by three water molecules, the backbone of residues 300–302, and the flavin ribityl, similar to P293, and three crystal waters in para-hydroxybenzoate hydroxylase suggest that P302, T47, and the waters in StyA are a vital component of the catalytic mechanism. The current optimized and validated StyA model provides a good starting point for elucidation of the structural basis of StyA ligand binding and catalysis. Novel insights in the binding of ligands to SMO/StyA, provided by the current protein model, will aid the rational design of mutants with specific, altered enantioselective properties.


   INTRODUCTION
TOP
 ABSTRACT
 INTRODUCTION
METHODS
 RESULTS
 DISCUSSION
 CONCLUSIONS
 ACKNOWLEDGEMENTS
 REFERENCES
 
Recently, interest in the class of two-component flavoenzyme hydroxylases has emerged (1Go,2Go). Two-component flavoenzyme hydroxylases consist of a smaller reductase component that uses NAD(P)H to reduce flavin and a larger hydroxylase component that uses the reduced flavin to perform the hydroxylation reaction. Styrene monooxygenase (SMO) from Pseudomonas sp. VLB120 is such a two-component flavoenzyme, consisting of a larger oxygenase domain StyA and a smaller reductase domain StyB, each with an independent catalytic center (2Go–4Go). StyA is not dependent on its native reductase; also, different reductases, chemical mediators that generate reduced flavin, and even direct electrochemical reduction on an electrode can be used as electron donor (3Go,5Go). The SMO enzyme is highly enantioselective, producing >99% enantiopure S-styrene oxide. Of further interest is the development of methods to rationalize and design mutants with altered substrate and/or enantioselectivity (6Go), which can provide a basis to explore possibly interesting applications in the chiral synthesis of epoxide building blocks for pharmaceuticals or agrochemicals (1Go). For this purpose, a structural model of the protein based on crystallographic data or from homology model building is indispensable. The class of flavoenzymes is well studied (2Go,7Go,8Go). For 50 FAD-containing oxidoreductase/hydroxylases there are x-ray structures in the Protein Data Bank (PDB) (9Go), but for SMO no x-ray structures have yet been published.

The current article describes a newly constructed, optimized, and validated homology structure of the StyA oxygenase component of SMO, deposited in the PDB under ID 2HD8 (9Go). Some general features of the binding cavity were described before, based on a preliminary StyA model (5Go), but that model was not capable of accommodating substrates in a catalytically active position, nor could it be used to describe or explain substrate binding in detail. For further optimization and for validation of a StyA protein model, experimental data on binding affinity and enzyme turnover for a range of substrates and for selected site-directed mutants were used. The construction and analysis of the current StyA structure comprise a vital step in the rationalization of activity and (enantio)selectivity of catalysis by StyA wild-type (WT) and mutants.


   METHODS
TOP
 ABSTRACT
 INTRODUCTION
 METHODS
RESULTS
 DISCUSSION
 CONCLUSIONS
 ACKNOWLEDGEMENTS
 REFERENCES
 
Sequence alignment
The primary sequence of StyA (accession code O50214) was obtained from the Sequence Retrieval System accessible from the EMBL (www.embl-heidelberg.de). The search for homologous sequences with resolved crystal structure was carried out using BLAST (10Go), yielding para-hydroxybenzoate hydroxylase (pHBH), PDB identifier (9Go) 1PBE (11Go), and 2-phenol-hydroxylase (PhHy), PDB identifier 1FOH (12Go). A multiple sequence alignment (MSA) of StyA with pHBH and PhHy was manually made based on conserved sequence motifs and substrate recognition sites (SRSs) in flavoproteins (see Table 1) (13Go), with the homology module of InsightII (Biosym, San Diego, CA). Further alignment was done on consensus using LOOPP (14Go). The MSA was refined using data on amino acids relevant for the binding of FAD (Table 2) (11Go,15Go–17Go). Ten different 3D strucures of StyA were generated with the restraint-based comparative modeling program Modeller 4.0 (18Go). From the structures generated by the Modeller program, the one with the best loop conformations and stereochemical parameters as determined by ProCheck (19Go) was selected for further modeling.


View this table:
[in this window]
[in a new window]
  		TABLE 1  Conserved sequence motifs in flavoproteins used for aligning StyA to pHBH

 

View this table:
[in this window]
[in a new window]
  		TABLE 2  FAD interacting residues in pHBH with correspondingly aligned residues in StyA, and secondary structure elements

 
Protein structure optimization
An optimization scheme employing steepest descents energy minimization (EM) and molecular dynamics (MD) simulations with positional restraints applied to different sets of atoms was used, which could be described as a "controlled release" method. The Gromacs MD package, version 3.1.4 (20Go,21Go), and the GROMOS 43a1 forcefield (22Go,23Go) were used with a 0.8/1.2-nm cutoff and a neighbor list update frequency of five steps. The LINCS algorithm (24Go) was used to constrain the length of all covalent bonds, and a time step of 2 fs was used. The protein structure, including the substrate styrene and cofactor FAD, was energy minimized in vacuum, first with harmonic position restraints on all nonhydrogen atoms with a force constant of 1000 kJ/mol/nm and then without position restraints. Subsequently, preequilibrated simple point-charge (SPC) water was added in a periodic dodecahedral box with a minimum distance of 1.2 nm between the protein and box edges, and the system was minimized again. Starting velocities were randomly generated from a 300 K Maxwell distribution. Temperature was maintained at 300 K and pressure at 1 bar by weak coupling to a bath using relaxation times of 1 ps and 0.1 ps, respectively.

During a 1-ps MD run the water was allowed to relax while the whole protein was position restrained with a force constant of 1000 kJ/mol/nm. In a series of short MD runs, subsequent parts of the protein were allowed to relax by removing the position restraints for the appropriate atoms. First, the side chains only of the residues not involved in either FAD or styrene binding were released for 1 ps. Next, also the backbone of the binding residues was released for 10 ps. During this whole procedure the positions of both FAD and styrene were also restrained. Finally, only the backbone atoms of the binding residues were restrained for 100 ps, while the rest of the protein as well as FAD and styrene were "free". The set of binding residues was defined as all residues that had at least one atom closer than 6 Å to one atom of either FAD or styrene in the nonoptimized structure. These were Ala13, Glu38, Tyr39, Arg43, Thr47, Val48, Glu113, Ala209, Val211, Leu220, Val222, Ala224, Arg233, Phe235, Asp295, Pro302, Gly305, Ala308, Asn309, and Phe409 and include all known FAD-interacting residues, as listed in Table 2, with the exception of Asp33. Several 1-ns-long MD simulations without position restraints were performed with styrene and FAD bound in the StyA protein: six started from the nonoptimized structure and six from the optimized structure, each with a different starting velocity randomly generated from a 300 K Maxwell distribution. Dynamic stability of the structure before and after the procedure was assessed by determining the atomic position root mean-square deviation (RMSD) with respect to the starting structure, averaged over the final 0.5 ns of each simulation, and differences in dynamic behavior are assessed from ED projections of the trajectories on the first two eigenvectors of all trajectories combined (25Go,26Go).

Candidate model structures were extracted from the 1-ns MD simulations started from the optimized structure by taking average structures over the first 100 ps and over the last 100 ps of the trajectories and energy minimizing to remove unphysical geometries that might result from the averaging, and the substrate styrene was removed. This resulted in a total of 12 candidate structures. Subsequently the substrate styrene was docked into each of the candidate structures. The optimal candidate was selected based on the possibility of reactivity of the docked conformations as judged from the distance to the reactive carbon atom FC4A in the FAD cofactor, the correct prochirality of the docked orientations of styrene as judged from the angle between the styrene plane and that of the isoalloxazine ring of FAD, and the specificity of the docked conformations as judged from the RMSD between different docked conformations within one structure.

StyA mutant structures
Structural models of 14 single (point) mutants of StyA were generated: V274Y/C/H, P275H/A, P302A/H, A298H, F409S, T12P/H, Q341R, L45C/H, and F173H. From the WT StyA structure, the selected amino-acid side chains were replaced by the mutant forms using the maximum overlap principle (retaining all overlapping atoms between WT and mutant residue) (27Go). The resulting structures were energy minimized to remove unfavorable interactions between the modified side chains and the surrounding protein.

Substrate binding conformation prediction
Twenty substances (styrene, 1,2-dihydronaphthaline, 3/4-bromo/chloro/methyl/nitro-stryrene, {alpha}/ß-methyl-styrene, indene, 1-phenyl-1-cyclohexen, 2/4-vinylpyridine, allylbenzol, 1,3-cyclo-hepta-/octadiene, naphthalene, vinylcyclopentane) were docked into the StyA protein structure using the automated docking program Gold (28Go). Of each compound, 50 docked conformations were generated and analyzed in terms of "active" binding, i.e., at the right distance and orientation with respect to the reactive carbon atom FC4A of FAD. A maximum distance of 6.5 Å was taken, based on FC4A-Op-Od bond lengths and estimates of distances in the reaction coordinate (15Go). The prediction of active binding was correlated with measured activities of these substrates to validate details of the active site geometry and chemical composition as well as the bound orientation of the FAD cofactor (29Go).

Substrate binding affinity prediction
The binding free energy ({Delta}G) of styrene and several of the other substrates for the StyA WT and mutant proteins were calculated using the linear interaction energy (LIE) approximation (30Go,31Go) as follows,

Formula 1(1)
where EVdW,bound and ECoul,bound are the van der Waals and the Coulomb interaction energies, respectively, between the substrate and its surroundings in the bound state, i.e., protein, FAD, and water, and EVdW,free and ECoul,free between substrate and surroundings in the free state, i.e., only water. These energies are obtained from MD simulations of the respective states, and error margins of the calculated {Delta}G are obtained from the respective energy fluctuations during the simulations. Using the solvent-accessible surface area of styrene calculated using the program MSMS (32Go) and the method proposed by Wang et al. (33Go) based on the correlation between weighted nonpolar desolvation ratio (WNDR) and {alpha}, an {alpha} of 0.98 was obtained. For the calculation of binding affinities for uncharged ligands without hydroxyl groups, the suggested value of 0.43 for ß is used (34Go). From the free energy of binding ({Delta}GCalc), the binding affinity is calculated as follows,

Formula 2(2)
where kB is Boltzmann's constant and T = 300 K.

Selected docked conformations according to the criteria mentioned above of styrene, 3-/4-/{alpha}-/ß-methylstyrene, 2-/4-vinylpyridine, vinylcyclopentane, 1,2-dihidronaphthaline, indene, and naphthalene, for the StyA-substrate complexes of WT and for styrene in the WT and the F235A/S, V274C, P275A/H, P302A/H, and F409S mutants were energy minimized, solvated, and equilibrated as described under Protein structure optimization. Four hundred picoseconds of production MD simulations were performed, and average interaction energies were collected over the final 300 ps. For interaction energies of the free substrate in solvent, the same procedure was followed, with 100 ps of production MD simulation and energies averaged over the last 50 ps.


   RESULTS
TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
DISCUSSION
 CONCLUSIONS
 ACKNOWLEDGEMENTS
 REFERENCES
 
StyA sequence alignment
A homology model of StyA (SMO catalytic domain) was constructed based on the x-ray crystallographic structure of pHBH, the closest homolog to StyA with 23% pairwise sequence identity. An overview of this protein model is presented in Fig. 1 A. The higher local similarity with PhHy (pairwise identity 22%) was also used for the MSA. Conserved motifs are summarized in Table 1. The conserved sequence motif GxGxxG (StyA residues 9–14) is thought to play an important role in the interaction with the ribose of the adenosyl moiety of FAD in pHBH (35Go) and in FAD-binding proteins in general (29Go,36Go). It is localized in the structurally conserved ß-{alpha}-ß unit at the C-terminal end of hydroxylases and monooxygenases (13Go). The short DG motif is located in loop ßA4-{alpha}H7 of the FAD-binding domain of pHBH and is situated near the cleft leading toward the active site (13Go); in StyA only the Gly171 is conserved. The GDx6P motif (StyA residues 294–302) is conserved in flavoprotein hydroxylases with a putative dual function in FAD and/or NAD(P)H binding (13Go) and is located in strand ßD3 (note that StyA binds FAD, not NAD(P)H). The conserved GxNx8L motif (13Go) corresponds to StyA residues 307–318 and is located in helix {alpha}H10. Further alignment was done on consensus using the LOOPP program output (14Go).


Figure 1
View larger version (69K):
[in this window]
[in a new window]
  		FIGURE 1  StyA, the substrate styrene (STY) and the cofactor FAD. (A) Overview of the whole StyA structure. Secondary structure elements are labeled according to the template pHBH (11Go). (B) Mutated residues around the StyA active site described in this study. Orange lines indicate the closest distances between the mutated residues and styrene or FAD. (C) Styrene binding in the StyA active site. The styrene binding cavity surface is shown, and the frontmost part is clipped to allow an inside view. Bound waters are highlighted with transparent blue spheres. The catalytic FC4A in the isoalloxazine ring of FAD and the vinyl group of Sty are highlighted in orange; the line in between shows their relative orientation. Hydrogen bonds among FAD, binding residues, and water are indicated with dashed yellow lines.

 
The MSA was refined by comparing the FAD binding residues in pHBH (11Go,13Go,15Go,16Go) (see Table 2) with residues in StyA. Ala13, Asp33, Arg43, Glu113, Arg233, Asp295, Pro302, Ala308, and Asn309 in StyA were found to correspond with Ser13, Glu32, Arg44, Gln102, Arg220, Asp286, Pro293, Leu299, and Asn300 in pHBH, respectively. The complete MSA of pHBH, PhHy, and StyA is presented in Table 3. In this final MSA between StyA and pHBH, 173 (38%) are strongly conserved (conservation score ≥8, e.g., Val-Ala, see Table 3).


View this table:
[in this window]
[in a new window]
  		TABLE 3  Multiple sequence alignment (MSA) among StyA (SMO_LV12), pHBH (1PBE), and PhHy (1FOH) as produced by CLUSTAL W (1.82)

 
In Table 4, the effects of the single amino-acid substitutions V274Y/C/H, P275H/A, P302A/H, A298H, F409S, T12P/H, Q341R, L45C/H, F173H in the StyA protein on the activity of the enzyme are listed as well as the closest separation between any of the atoms of the corresponding amino acids in the model structure and styrene and FAD. These respective residues are shown in Fig. 1 B in relation to the location of styrene and FAD in the model. For all residues closer than 5 Å to either styrene or FAD (T12, L45, V274, P302 and F409), a large effect of at least one of the substitutions (<50% activity versus wild-type) was observed. Conversely, mutations in residues at more than 5.5 Å from styrene and FAD (F173, P275, A298, and Q341) showed small effects of all substitutions (at least 50% activity versus wild-type). In addition, for mutations in six positions, the stabilization of bound reduced FAD was determined (see Table 4), and it was changed by at least one mutation in five of them (T12, L45, F173, F235, and P302), all closer than 5.5 Å to FAD in the model. With the effects of all 10 mutated residue positions in StyA explained by distances in the StyA model structure, the alignment of the StyA sequence to the template structures pHBH and PhHy is properly validated, especially in the most important areas of the enzyme: the binding sites of styrene and FAD.


View this table:
[in this window]
[in a new window]
  		TABLE 4  Mutants of StyA, shortest distances between atoms of the mutated residues and the styrene (Sty) and FAD binding pockets in the protein structure, measured relative activities, and binding affinities of reduced FAD, average distances between FC4A of FAD and the vinyl carbons (CH) of styrene during LIE-MD simulations, and calculated free energies ({Delta}GCalc) for styrene in selected mutants

 
StyA structure and optimization
The initial protein model structure was found to be moderately stable in EM and MD simulations (20Go–23Go), with C{alpha} atoms RMSD with respect to the starting conformation rising during six 1-ns free MD simulations to plateau values of 4.2 ± 0.5 Å. The C{alpha} atom RMSD between the protein structure before and after the controlled-release optimization was 2.8 Å for the whole protein (3.4 Å for all atoms) and 0.9 Å for the C{alpha}s of the active site only (1.7 Å for all atoms). The optimized protein structure was more stable with corresponding RMSD plateau values of 3.6 ± 0.3 Å during six 1-ns free MD simulations, which indicates a good dynamic stability for a protein of this size and type on a timescale of nanoseconds. The C{alpha} RMSD plateau of the binding residues was reduced from 2.0 ± 0.3 Å in the simulations before optimization to 1.5 ± 0.3 Å in the simulations after optimization. This indicates a good dynamic stability of the active site region. Six candidate model structures were extracted from six simulations of the optimized structure by averaging the conformations of the first 100 ps, and another six structures from the final 100 ps. ProCheck analysis (19Go) of the distribution of backbone dihedral angles of these candidate structures indicated the structures after the optimization simulations still possessed good stereochemical quality (19Go): in all cases 94 ± 1% of residues were in core or allowed regions. Table 5 shows a summary of Ramachandran distributions at different stages in the optimization. For the initial model, 78% of main-chain bond lengths were within limits, and two were off the graph; for all other models 100% were within limits. As a general trend, it was observed that, during the five stages of the controlled-release optimization, the quality decreased somewhat; i.e., the Ramachandran plot shifted from 74% core and 20% allowed to 63% and 32%, respectively. However, much of this decrease was made up again by taking the minimized average structure over 100 ps of simulation in the final stage.


View this table:
[in this window]
[in a new window]
  		TABLE 5  ProCheck Ramachandran distributions and other quality factors for the StyA structures in various stages of optimization and for the final model

 
In Fig. 2, the essential dynamics (ED) projection (25Go,26Go) is shown for the changes in the protein and active site structure on optimization and during the MD simulations. Distances between two points in this plot are closely related to the familiar RMSD measure, but in addition, differences between the trajectories visualize motion in different directions. Motions along the two axes in the plot correspond to the two largest collective modes of motion of the enzyme, also known as "essential modes" (25Go). This illustrates the difference in behavior between the simulations started from the nonoptimized structure and those started from the optimized structure. The stretched appearance of the traces of the nonoptimized simulations is caused by a ballistic type of motion indicative of high internal strain in the protein structure and corresponds to the high RMSD plateau level of 4.2 ± 0.5 Å for the whole protein and 2.0 ± 0.3 Å for the active site. In the projection of the whole protein, the difference in the nonoptimized and optimized structures of 2.9 Å can be seen as a shift of the starting point of the corresponding simulations, whereas in contrast, in projection of the active site, the smaller structural shift of 0.9 Å is hardly visible. The optimized structure, in contrast, shows a much more diffusive nature of the projections of the simulations, which is indicative of a protein structure with native-like properties (37Go) and corresponds to the lower RMSDs of 3.6 ± 0.3 Å and 1.5 ± 0.3 Å for protein and active site, respectively.


Figure 2
View larger version (32K):
[in this window]
[in a new window]
  		FIGURE 2  MD simulations for optimization and evaluation of optimized structures projected onto the two essential modes (the two largest principal eigenvectors) of (a) C{alpha} positions of all 415 residues of the protein; and (b) C{alpha} positions of the binding residues only. The dotted lines indicate differences in conformation before and after the controlled-release optimization process. These are prominent for the whole protein but hardly visible for the active site. Note the difference in scale between the plots.

 
Of all optimized StyA protein structures, the predictions from automated docking showed one structure to bind styrene in a single tight cluster with a maximum RMSD between members of the cluster of 0.5 Å, in a catalytically active position with respect to the FAD cofactor reactive center (<6.5 Å between the ethylene of styrene and FC4A, see Fig. 1 C), and an orientation prochiral to the formation of S-styrene oxide exclusively. In the other structures, styrene was bound in a noncatalytically relevant position, in widely varying positions, or in a mixed prochiral orientation. No protein structure was found that would support exclusive pro(R) orientation. The one structure showing the correct, tight, and pro(S) binding characteristics for styrene was chosen for subsequent validation. The final StyA structure is shown in Fig. 1 A, and details of the styrene binding pocket and the bound orientation of styrene are shown in Fig. 1 C.

StyA structure validation
Binding free energies ({Delta}G) for styrene in the wild-type (WT) StyA and the mutants F235A/S, V274C, P275A/H, P302A/H, and F409S were calculated from MD simulations using the LIE method, using Eq. 1 (30Go,31Go). The docked orientation of styrene in the StyA WT structure was also used as the starting orientation for the MD simulations in the mutant structures (see Fig. 1, A and C). Calculated free energies ({Delta}G) of binding are listed in Table 4. The average distances between the catalytically active carbon FC4A of FAD and the vinyl group of styrene (see Fig. 1 C) are listed in Table 4. A distance of <6.5 Å can be considered reactive (15Go).

In Table 6, experimental activities classified as high (>20%), low (1–20%), or no activity (<1%) relative to styrene for the 20 selected styrene-homologous substrates are listed. Sixteen of these could be docked at a distance to the catalytically active carbon FC4A and in an orientation with respect to the isoalloxazine ring of FAD, corresponding to an active conformation. However, 3-bromo-, 4-chloro-, and 3- and 4-nitrostyrene were predicted by the docking algorithm to bind in an orientation relative to FAD in which the aromatic rings were stacked.


View this table:
[in this window]
[in a new window]
  		TABLE 6  Substrates and other compounds with known activity in StyA*

 
For styrene, 3-, 4-, {alpha}- and ß-methylstyrene, 1,2-dihydronaphthaline, naphthalene, indene, 2- and 4- vinylpyridine, vinylcyclopentane, allylbenzol, and 1,3-cycloheptadiene binding free energies ({Delta}G) were determined using the LIE method and are listed in Table 6. In Fig. 3, it can be seen that in general the substrates with low activity have low predicted binding affinities (≤–8.7 kcal/mol), whereas most with higher activity ("+" or "o") have higher affinity (≥–8.7 kcal/mol). For the individual substrates, however, the error margins are such that the differences are not significant. Finally, the relatively high inhibition of styrene catalysis by 4- and ß-methylstyrene and the relatively low inhibition by {alpha}-methylstyrene determined experimentally, are reproduced by the calculated (relative) binding affinities of these substrates (see Table 6). The final StyA structure has been deposited in the PDB under ID 2HD8.


Figure 3
View larger version (8K):
[in this window]
[in a new window]
  		FIGURE 3  Calculated binding free energies ({Delta}Gcalc) for 13 of the substrates, as listed in Table 6, following the ordering and classification of experimental activities as in Table 6.

 
Structural features of styrene, FAD, and water binding in StyA
Styrene binds in the StyA binding pocket stacked with Phe235 at the bottom and with Val303 at the top, as is shown in Fig. 1 C. Leu220 and the di-methyl-phenyl ring of FAD make up the left and right sides of the binding pocket, respectively. Finally, Ile198 makes up the front, and Val222 and the hydroxyl group of Thr47 the rear side. Of these six residues, five are at ~3.5 Å from styrene, whereas Val303 and FAD are at ~4 Å. Thr47 is the only polar side chain in the styrene binding pocket. It is in a position close to the catalytic FC4A of FAD and the ethylene group of styrene, where it may be involved in stabilization of an intermediate in the reaction. Ten additional residues can be found that have at least one atom within 6.0 Å of styrene, shown as well in Fig. 1 C, but none of these has any hydrophilic parts within 6 Å of styrene. Furthermore, the backbone carbonyls of Ile198, Leu220, Pro302, and Val303 are located inside the pocket. These active site residues and secondary structure elements are summarized in Table 7. In the strongly hydrophobic StyA styrene binding cavity, no water is present at all. A putative substrate access channel can be identified going through the FAD binding site, which is blocked when FAD is bound.


View this table:
[in this window]
[in a new window]
  		TABLE 7  Active site residues with at least one atom within 6 Å from at least one atom of the bound styrene*

 
The interactions of FAD with StyA are shown in Fig. 4. The isoalloxazine ring is H-bonded to Thr47, Val48, and Asn309, and additional H-bonds are formed with Ala13, Asn46, and Ala308 (not shown). Hydrophobic contacts are made with Gln306 and Gly307, and additionally with Leu45 and Gly305 (not shown). The flavin ribityl hydroxyl groups are H-bonded to the pyrophosphate and to Glu113, Asp295, and some water molecules (not all waters are shown). The pyrophosphate is H-bonded to Thr12, via water to Arg34, and in addition to Arg43 (not shown). The adenosine ribose is H-bonded to the backbone of Arg34 and additionally to Asp33 and the backbone of Ala10 and Gly11. Finally, the adenine double ring is H-bonded to Leu139 and Lys174 (not shown) and makes hydrophobic contacts with Arg34. The front (as in Fig. 4) of the FAD binding pocket of StyA is open and accessible to the solvent; ~16 water molecules hydrate the bound FAD, several of which are worth mentioning. A string of three waters runs from the active site along the backbone carbonyls of Pro302, Asp301, Glu306, and Gly307 to the ribityl of FAD. Two waters meditate H-bonds between Arg34 and Asp295 and the ribityl and pyrophosphate, and three more (not shown) extend this H-bonded network. An additional three waters (not shown) mediate H-bonds between the adenine and Leu138, Leu139, and Glu178.


Figure 4
View larger version (18K):
[in this window]
[in a new window]
  		FIGURE 4  FAD and water binding in the StyA active site. Residues in the read are labeled with a smaller font. The substrate styrene is shown as well. Hydrogen bonds among FAD, selected binding residues, and selected water molecules are indicated with dashed yellow lines. Waters are highlighted with transparent light-blue spheres.

 

   DISCUSSION
TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
CONCLUSIONS
 ACKNOWLEDGEMENTS
 REFERENCES
 
StyA structure building and optimization
We have examined the substrate binding of the oxygenase domain (StyA) of the styrene monooxygenase enzyme from Pseudomonas sp. VLB120, in its wild-type form and several site-directed mutants, using a newly constructed homology model. For this purpose, ligand binding orientations and affinities were predicted for the wild-type StyA and 17 mutants and for 20 other substrates for the wild-type StyA. The predictions were used to rationalize the observed relative activities of styrene catalysis determined experimentally.

The StyA protein model was based on the crystallized para-hydroxybenzoate hydroxylase (pHBH) and 2-phenol hydroxylase (PhHy). The pairwise identity between pHBH and StyA was high enough (23%) to make homology modeling feasible and was further facilitated by the strong sequence similarity between pHBH and StyA of ~38% (conservation score ≥8 in the MSA of Table 3). Strongly conserved structural and functional features within the oxidoreductase enzyme family can be traced down to several key patterns of conservation (Table 1) (13Go,29Go), certainly for the most important active site region of the enzyme (Table 2) (11Go,15Go–17Go,38Go), lending additional reliability to the alignment. All known pHBH active site residues were very strongly conserved throughout the MSA (conservation score ≥8 in Table 3), and their counterparts in the StyA sequence (Table 2) could be incorporated in the active site region of the StyA protein structure (Fig. 4). For several additional residues that are known to influence the StyA catalytic function (5Go), their close proximity to the substrate and/or cofactor is shown (Table 4).

This first StyA homology model has been demonstrated to be stable in MD simulations (Fig. 2). For the whole StyA protein, the change in structure from the nonoptimized to the optimized structure is rather large (2.9 Å RMSD on C{alpha}), whereas for the active site only, the concomitant structural shift is relatively small (1.7 Å RMSD on C{alpha}). In addition, from the ED projections (Fig. 2), it can be seen that the structural changes during free MD simulation of the nonoptimized structure are ballistic and indicative of high levels of internal strain in the protein structure, in contrast to the simulations of the optimized structure, which are diffusive, which is a characteristic of stable, native-state protein structures (37Go). These are clear effects of the controlled-release optimization method employed that minimizes structural changes in the active site region while allowing the rest of the protein structure to relax significantly.

StyA structure validation with experimental data
For the 13 mutants of StyA studied (Table 4 and Fig. 1 B), a clear correlation was observed between the distance of a mutation site to the closest of the styrene and FAD binding sites and the maximum effect of the substitutions on the measured enzyme activity, indicating a good overall alignment of StyA with the template structure. The P302 substitutions are particularly worth mentioning. The P302H mutant has a low measured activity for styrene (<10% versus WT), and styrene is predicted to bind very tightly ({Delta}G < –11 kcal/mol, Kd < 1 nM) but too far away for catalysis (~10 Å) (see Table 4). It appears that the histidine at position 302, which is located just above the isoalloxazine ring of FAD, may block access to the catalytic site in StyA (see Fig. 1 C). Alternately, the P302H mutation could affect the structure or stability of the ßD3-{alpha}H10 loop. The P302A mutant has a low measured activity (10%), but the predicted binding affinity of styrene is comparable to WT ({Delta}G {approx} –9 kcal/mol, Kd {approx} 100 nM), and it binds close enough for catalysis (~6.5 Å). Possibly, changing the stiff proline for the more flexible alanine adversely affects the protein's stability, as is indeed also seen from the increased appearance of inclusion bodies (5Go). In the V274C mutant, styrene is predicted to bind strongly ({Delta}G = –10.1 kcal/mol, Kd = 24 nM), but at too large distance (>7 Å) to be actively metabolized (Table 4).

All 21 compounds studied (Table 6) could be docked into the active site cavity. For 13 compounds (Table 6), affinities were predicted using the LIE method. As a trend, the active compounds showed a higher predicted affinity, although for individual cases the error margins are such that significance is minimal. In addition, it must be borne in mind that in other cases too, the actual binding affinity could well be the limiting factor for efficient turnover. Only 3-bromo-, 4-chloro-, and 3- and 4-nitrostyrene were predicted to bind in an orientation unsuited for metabolism, whereas they were nevertheless metabolized in our experiment. This might be attributed to the chemical nature of the relatively large, electron-rich, and polarizable bromo-, chloro- and nitro- substituents, whose properties are difficult to describe appropriately in a force field or scoring function. The two substrates with lowest predicted affinity (>250 nM; cycloheptadiene and vinylcyclopentane) lack the characteristic aromatic six-membered ring and are also experimentally inactive. The other substrate that lacks the aromatic ring is cyclooctadiene, which, although predicted to bind in the right place and at good affinity, is not catalyzed (Table 6). This leads to the tentative conclusion that the aromatic six-membered ring might be important for binding and essential for catalysis. The final StyA structure has been deposited in the PDB under ID 2HD8.

The StyA active site and styrene and FAD binding regions
A close look at the styrene binding cavity and active site revealed that the hydroxyl group of Tyr47 and carbonyl group of Pro302 are at ~4 Å from the catalytic FC4A of FAD and the ethylene group of styrene in an otherwise completely hydrophobic binding cavity (at the rear and front, respectively, in Fig. 1 C). The corresponding Pro293 in pHBH (Table 2) is implicated in stabilization of intermediate state(s) of the reaction (15Go), which suggests that the Pro302 and Thr47 are important for catalytic activity in StyA. Most discussion in literature on the importance of water molecules in the pHBH catalytic function, appears to have been focused on a hydrogen-bonded water network in the active site cavity bridging the substrate and the putative catalytic Pro293 with the bulk water through the FAD binding pocket (39Go,40Go). Likewise, the strongly hydrophobic StyA styrene binding cavity of the StyA model contains no water molecules (Fig. 4), but three water molecules were hydrogen-bonded to the FAD side of the P302 carbonyl and form a hydrogen-bonded network with the backbone of residues 301, 302, 306, and 307 and the flavin ribityl of FAD (Fig. 4). In the crystal structures of pHBH an arrangement of water molecules is present, involving the active-site P293 carbonyl, which is very similar to that found in the current model. It should be noted that in the StyA model the crystal waters of the pHBH template structure were excluded, and the StyA protein structure was resolvated and equilibrated in water for the MD simulations. These water molecules have been suggested to be involved in the catalytic mechanism of pHBH for proton donation as well as stabilization of intermediate state(s) (40Go), which leads to the compelling suggestion that also in StyA, these waters are critical for catalysis.


   CONCLUSIONS
TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 CONCLUSIONS
ACKNOWLEDGEMENTS
 REFERENCES
 
A first homology model of StyA has been constructed and optimized using molecular dynamics simulations. For all key residues for ligand and cofactor binding and for catalytic activity, known from the template pHBH protein, corresponding counterparts in the StyA sequence and structure could be identified. Known experimental data on enantioselectivity and catalytic activity of styrene in StyA wild types and mutants and of different substrates in the wild-type StyA corresponded well with predicted substrate binding orientations and affinities based on this StyA homology model. The StyA active site residues Thr47 and Pro302 are critical in StyA catalysis and are possibly involved in the stabilization of reaction intermediate(s). The remainder of the styrene binding pocket is exclusively hydrophobic and consists of Phe235, Val303, Leu220, the di-methyl-phenyl ring of FAD, Ile198, and Val222, where styrene binds stacked between Phe235 and Val303. The model structure of StyA is deposited in the PDB under ID 2HD8. The current StyA protein model thus provides a good starting point for further validation of the role of the active site residues and FAD cofactor binding residues in ligand binding and catalysis. Novel insights in the active site topology and in the binding of ligands in SMO/StyA may facilitate the rational design of mutants with specific, altered enantioselective properties.


   ACKNOWLEDGEMENTS
TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 CONCLUSIONS
 ACKNOWLEDGEMENTS
REFERENCES
 
We thank Dr. Chris Oostenbrink and Dr. Katja Otto for helpful discussions and Simon Allioth for assisting with protein purifications.

Financial support from the European Community grant "Industrial biocatalysis with new oxygenases in a novel electro-enzyme reactor" in the "Quality of Life and Management of Living Resources" program is acknowledged.

Submitted on May 9, 2006; accepted for publication July 12, 2006.


   REFERENCES
TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 CONCLUSIONS
 ACKNOWLEDGEMENTS
 REFERENCES
 
1. Panke, S., M. G. Wubbolts, A. Schmid, and B. Witholt. 2000. Production of enantiopure styrene oxide by recombinant Escherichia coli synthesizing a two-component styrene monooxygenase. Biotechnol. Bioeng. 69:91–100.[CrossRef][Medline]

2. Kantz, A., F. Chin, N. Nallamothu, T. Nguyen, and G. T. Gassner. 2005. Mechanism of flavin transfer and oxygen activation by the two-component flavoenzyme styrene monooxygenase. Arch. Biochem. Biophys. 442:102–116.[CrossRef][Medline]

3. Hollmann, F., P.-C. Lin, B. Witholt, and A. Schmid. 2003. Stereospecific biocatalytic epoxidation: the first example of direct regeneration of a FAD-dependent monooxygenase for catalysis. J. Am. Chem. Soc. 125:8209–8217.[CrossRef][Medline]

4. Otto, K., K. Hofstetter, M. Röthlisberger, B. Witholt, and A. Schmid. 2004. Biochemical characterization of StyAB from Pseudomonas sp. VLB120 as a two-component flavin-diffusible monooxygenase. J. Bacteriol. 186:5292–5302.[Abstract/Free Full Text]

5. Hofstetter, K. 2006. Biocatalytic asymmetric epoxidation of styrene and derivatives using isolated styrene monooxygenase. PhD thesis. ETH, Zürich. http://e-collection.ethbib.ethz.ch/show?type=diss&nr=16400.

6. Schmid, A., K. Hofstetter, H.-J. Feiten, F. Hollmann, and B. Witholt. 2001. Integrated biocatalytic synthesis on gram scale: the highly enantioselective preparation of chrial oxirans with styrene monooxygenase. Adv. Synth. Catal. 343:732–737.[CrossRef]

7. Entsch, B., L. J. Cole, and D. P. Ballou. 2005. Protein dynamics and electrostatics in the function of p-hydroxybenzoate hydroxylase. Arch. Biochem. Biophys. 433:297–311.[CrossRef][Medline]

8. Entsch, B., and W. J. H. Van Berkel. 1995. Flavoprotein structure and mechanism. 1. structure and mechanism of para-hydroxybenzoate hydroxylase. FASEB J. 9:476–483.[Abstract]

9. Berman, H. M., J. Westbrook, Z. Feng, G. Gilliland, T. N. Bhat, H. Weissig, I. N. Shindyalov, and P. E. Bourne. 2000. The Protein Data Bank. Nucleic Acids Res. 28:235–242.[Abstract/Free Full Text]

10. Altschul, S. F., W. Gish, W. Miller, E. W. Meyers, and D. J. Lipman. 1990. Basic local alignment search tool. J. Mol. Biol. 215:403–410.[CrossRef][Medline]

11. Schreuder, H. A., P. A. Prick, R. K. Wierenga, G. Vriend, K. S. Wilson, W. G. Hol, and J. Drenth. 1989. Crystal structure of the p-hydroxybenzoate hydroxylase-substrate complex refined at 1.9 A resolution. Analysis of the enzyme-substrate and enzyme-product complexes. J. Mol. Biol. 208:679–696.[CrossRef][Medline]

12. Enroth, C., H. Neujahr, G. Schneider, and Y. Lindqvist. 1998. The crystal structure of phenol hydroxylase in complex with FAD and phenol provides evidence for a concerted conformational change in the enzyme and its cofactor during catalysis. Structure. 6:605–617.[Medline]

13. Eppink, M. H., H. A. Schreuder, and W. J. Van Berkel. 1997. Identification of a novel conserved sequence motif in flavoprotein hydroxylases with a putative dual function in FAD/NAD(P)H binding. Protein Sci. 6:2454–2458.[Abstract]

14. Meller, J., and R. Elber. 2001. Linear programming optimization and a double statistical filter for protein threading protocols. Proteins. 45:241–261.[CrossRef][Medline]

15. Ridder, L., A. J. Mulholland, I. Rietjens, and J. Vervoort. 2000. A quantum mechanical/molecular mechanical study of the hydroxylation of phenol and halogenated derivatives by phenol hydroxylase. J. Am. Chem. Soc. 122:8728–8738.[CrossRef]

16. Nishiya, Y., and T. Imanaka. 1996. Analysis of the interaction between the Arthrobacter sarcosine oxidase and the coenzyme flavin adenine dinucleotide by site-directed mutagenisis. Appl. Environ. Microbiol. 62:2405–2410.[Abstract]

17. Gatti, D., B. Pafley, M. Lah, B. Entsch, V. Maasey, D. Ballou, and M. Ludwig. 1994. The mobile flavin of 4-OH benzoate hydroxylase. Science. 266:110–114.[Abstract/Free Full Text]

18. Sali, A., and T. L. Blundell. 1993. Comparative protein modelling by satisfaction of spatial restraints. J. Mol. Biol. 234:779–815.[CrossRef][Medline]

19. Laskowski, R. A., M. W. MacArthur, D. S. Moss, and J. M. Thornton. 1993. PROCHECK: a program to check the stereochemical quality of protein structures. J. Appl. Crystallogr. 26:283–291.[CrossRef]

20. Berendsen, H. J. C., D. van der Spoel, and R. van Drunen. 1995. Gromacs: A message-passing parallel molecular dynamics implementation. Comput. Phys. Commun. 91:43–56.[CrossRef]

21. Lindahl, E., B. Hess, and D. van der Spoel. 2001. GROMACS 3.0: a package for molecular simulation and trajectory analysis. J. Mol. Model. (Online). 7:306–317.

22. Daura, X., A. E. Mark, and W. F. v. Gunsteren. 1998. Parametrisation of aliphatic CHn united atoms of GROMOS96 force field. J. Comput. Chem. 19:535–547.[CrossRef]

23. Van Gunsteren, W. F., S. R. Billeter, A. A. Eising, P. H. Hunenberger, P. Kruger, A. E. Mark, W. R. P. Scott, and I. G. Tironi. 1996. Biomolecular Simulation: Gromos96 Manual and User Guide. Biomos b.v., Zurich, Groningen.

24. Hess, B., H. Bekker, H. J. C. Berendsen, and J. G. E. M. Fraaije. 1997. LINCS: A linear constraint solver for molecular simulations. J. Comput. Chem. 18:1463–1472.[CrossRef]

25. Amadei, A., A. B. Linssen, and H. J. Berendsen. 1993. Essential dynamics of proteins. Proteins. 17:412–425.[CrossRef][Medline]

26. Hünenberger, P. H., A. E. Mark, and W. F. Van Gunsteren. 1995. Fluctuation and cross-correlation analysis of protein motions observed in nanosecond molecular dynamics simulations. J. Mol. Biol. 252:492–503.[CrossRef][Medline]

27. Summers, N. L., and M. Karplus. 1991. Modeling of side chains, loops, and insertions in proteins. Methods Enzymol. 202:156–204.[Medline]

28. Jones, G., P. Willett, R. C. Glen, A. R. Leach, and R. Taylor. 1997. Development and validation of a genetic algorithm for flexible docking. J. Mol. Biol. 267:727–748.[CrossRef][Medline]

29. Babor, M., V. Sobolev, and M. Edelman. 2002. Conserved positions for ribose recognition: Importance of water bridging interactions among ATP, ADP and FAD-protein complexes. J. Mol. Biol. 323:523–532.[CrossRef][Medline]

30. Åqvist, J., C. Medina, and J. E. Samuelsson. 1994. A new method for predicting binding affinity in computer-aided drug design. Protein Eng. 7:385–391.[Abstract/Free Full Text]

31. Hansson, T., J. Marelius, and J. Åqvist. 1998. Ligand binding affinity prediction by linear interaction energy methods. J. Comput. Aided Mol. Des. 12:27–35.[CrossRef][Medline]

32. Sanner, M. F., A. J. Olson, and J. C. Spehner. 1996. Reduced surface: an efficient way to compute molecular surfaces. Biopolymers. 38:305–320.[CrossRef][Medline]

33. Wang, W., J. Wang, and P. A. Kollman. 1999. What determines the van der Waals coefficient beta in the LIE (linear interaction energy) method to estimate binding free energies using molecular dynamics simulations? Proteins. 34:395–402.[CrossRef][Medline]

34. Åqvist, J., and T. Hansson. 1996. On the validity of electrostatic linear response in polar solvents. J. Phys. Chem. 100:9512–9521.[CrossRef]

35. Wierenga, R. K., P. Terpstra, and W. G. J. Hol. 1986. Prediction of the occurrence of the ADP-binding beta-alpha-beta-fold in proteins, using an amino-acid-sequence fingerprint. J. Mol. Biol. 187:101–107.[CrossRef][Medline]

36. Kleiger, G., and D. Eisenberg. 2002. GXXXG and GXXXA motifs stabilize FAD and NAD(P)-binding Rossmann folds through C-alpha-H ··· O hydrogen bonds and van der Waals interactions. J. Mol. Biol. 323:69–76.[CrossRef][Medline]

37. Hess, B. 2002. Convergence of sampling in protein simulations. Phys. Rev. E. 65:31911.[CrossRef]

38. Palfey, B. A., D. P. Ballou, and V. Massey. 1997. Flavin conformational changes in the catalytic cycle of p-hydroxybenzoate hydroxylase substituted with 6-azido- and 6-aminoflavin adenine dinucleotide. Biochemistry. 36:15713–15723.[CrossRef][Medline]

39. Ortiz-Maldonado, M., B. Entsch, and D. P. Ballou. 2004. Oxygen reactions in p-hydroxybenzoate hydroxylase utilize the H-bond network during catalysis. Biochemistry. 43:15246–15257.[CrossRef][Medline]

40. Schreuder, H. A., W. G. Hol, and J. Drenth. 1990. Analysis of the active site of the flavoprotein p-hydroxybenzoate hydroxylase and some ideas with respect to its reaction mechanism. Biochemistry. 29:3101–3108.[CrossRef][Medline]

41. Clamp, M., J. Cuff, S. M. Searle, and G. J. Barton. 2004. The Jalview Java Alignment Editor. Bioinformatics. 12:426–427.





This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow All Versions of this Article:
biophysj.106.088633v1
91/9/3206    most recent
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Feenstra, K. A.
Right arrow Articles by Vermeulen, N. P. E.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Feenstra, K. A.
Right arrow Articles by Vermeulen, N. P. E.

HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
Copyright © 2006 by the Biophysical Society.