RESUMO
The use of enzymes in preparative biocatalysis often requires tailoring enzyme selectivity by protein engineering. Herein we explore the use of computational library design and molecular dynamics simulations to create variants of limonene epoxide hydrolase that produce enantiomeric diols from meso-epoxides. Three substrates of different sizes were targeted: cis-2,3-butene oxide, cyclopentene oxide, and cis-stilbene oxide. Most of the 28 designs tested were active and showed the predicted enantioselectivity. Excellent enantioselectivities were obtained for the bulky substrate cis-stilbene oxide, and enantiocomplementary mutants produced (S,S)- and (R,R)-stilbene diol with >97 % enantiomeric excess. An (R,R)-selective mutant was used to prepare (R,R)-stilbene diol with high enantiopurity (98 % conversion into diol, >99 %â ee). Some variants displayed higher catalytic rates (kcat ) than the original enzyme, but in most cases KM values increased as well. The results demonstrate the feasibility of computational design and screening to engineer enantioselective epoxide hydrolase variants with very limited laboratory screening.
Assuntos
Álcoois/metabolismo , Epóxido Hidrolases/metabolismo , Álcoois/química , Sítios de Ligação , Biocatálise , Epóxido Hidrolases/genética , Cinética , Simulação de Dinâmica Molecular , Mutagênese , Proteínas Recombinantes/biossíntese , Proteínas Recombinantes/isolamento & purificação , Estereoisomerismo , Estilbenos/química , Estilbenos/metabolismo , Especificidade por SubstratoRESUMO
The P450 monooxygenase CYP109A2 from Bacillus megaterium DSM319 was previously found to convert vitamin D3 (VD3) to 25-hydroxyvitamin D3. Here, we show that this enzyme is also able to convert testosterone in a highly regio- and stereoselective manner to 16ß-hydroxytestosterone. To reveal the structural determinants governing the regio- and stereoselective steroid hydroxylation reactions catalyzed by CYP109A2, two crystal structures of CYP109A2 were solved in similar closed conformations, one revealing a bound testosterone in the active site pocket, albeit at a nonproductive site away from the heme-iron. To examine whether the closed crystal structures nevertheless correspond to a reactive conformation of CYP109A2, docking and molecular dynamics (MD) simulations were performed with testosterone and vitamin D3 (VD3) present in the active site. These MD simulations were analyzed for catalytically productive conformations, the relative occurrences of which were in agreement with the experimentally determined stereoselectivities if the predicted stability of each carbon-hydrogen bond was taken into account. Overall, the first-time determination and analysis of the catalytically relevant 3D conformation of CYP109A2 will allow for future small molecule ligand screening in silico, as well as enabling site-directed mutagenesis toward improved enzymatic properties of this enzyme.
Assuntos
Bacillus megaterium , Sistema Enzimático do Citocromo P-450 , Sistema Enzimático do Citocromo P-450/metabolismo , Bacillus megaterium/metabolismo , Hidroxilação , Cristalografia por Raios X , Esteroides/metabolismo , Simulação de Dinâmica Molecular , Colecalciferol/metabolismo , Testosterona/metabolismoRESUMO
CYP105AS1 is a cytochrome P450 from Amycolatopsis orientalis that catalyzes monooxygenation of compactin to 6-epi-pravastatin. For fermentative production of the cholesterol-lowering drug pravastatin, the stereoselectivity of the enzyme needs to be inverted, which has been partially achieved by error-prone PCR mutagenesis and screening. In the current study, we report further optimization of the stereoselectivity by a computationally aided approach. Using the CoupledMoves protocol of Rosetta, a virtual library of mutants was designed to bind compactin in a pro-pravastatin orientation. By examining the frequency of occurrence of beneficial substitutions and rational inspection of their interactions, a small set of eight mutants was predicted to show the desired selectivity and these variants were tested experimentally. The best CYP105AS1 variant gave >99% stereoselective hydroxylation of compactin to pravastatin, with complete elimination of the unwanted 6-epi-pravastatin diastereomer. The enzyme-substrate complexes were also examined by ultrashort molecular dynamics simulations of 50 × 100 ps and 5 × 22 ns, which revealed that the frequency of occurrence of near-attack conformations agreed with the experimentally observed stereoselectivity. These results show that a combination of computational methods and rational inspection could improve CYP105AS1 stereoselectivity beyond what was obtained by directed evolution. Moreover, the work lays out a general in silico framework for specificity engineering of enzymes of known structure.