Enantioselective Benzylic Hydroxylation Catalysed by P450 Monooxygenases: Characterisation of a P450cam Mutant Library and Molecular Modelling

Computationally guided bioengineering of the active site, substrate access pathway, and water channels of thermostable cytochrome P450, CYP175A1, for catalyzing the alkane hydroxylation reaction

Understanding structure-function relationships in proteins is pivotal in their development as industrial biocatalysts. In this regard, rational engineering of protein active site access pathways and various tunnels and channels plays a central role in designing competent enzymes with high stability and enhanced efficiency. Here, we report the rational evolution of a thermostable cytochrome P450, CYP175A1, to catalyze the C-H activation reaction of longer-chain alkanes. A strategy combining computational tools with experiments has shown that the substrate scope and enzymatic activity can be enhanced by rational engineering of certain important channels such as the substrate entry and water channels along with the active site of the enzyme. The evolved enzymes showed an improved catalytic rate for hexadecane hydroxylation with high regioselectivity. The Q67L/Y68F mutation showed binding of the substrate in the active site, water channel mutation L80F/V220T showed improved catalytic activity through the peroxide shunt pathway and substrate entry channel mutation W269F/I270A showed better substrate accessibility to the active pocket. All-atom MD simulations provided the rationale for the inactivity of the wild-type CYP175A1 for hexadecane hydroxylation and predicted the above hot-spot residues to enhance the activity. The reaction mechanism was studied by QM/MM calculations for enzyme-substrate complexes and reaction intermediates. Detailed thermal and thermodynamic stability of all the mutants were analyzed and the results showed that the evolved enzymes were thermally stable. The present strategy showed promising results, and insights gained from this work can be applied to the general enzymatic system to expand substrate scope and improve catalytic activity.

Regio- and stereoselective steroid hydroxylation by CYP109A2 from Bacillus megaterium explored by X-ray crystallography and computational modeling

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.

Mechanism-Guided Computational Design of ω-Transaminase by Reprograming of High-Energy-Barrier Steps

ω-Transaminases (ω-TAs) show considerable potential for the synthesis of chiral amines. However, their low catalytic efficiency towards bulky substrates limits their application, and complicated catalytic mechanisms prevent precise enzyme design. Herein, we address this challenge using a mechanism-guided computational enzyme design strategy by reprograming the transition and ground states in key reaction steps. The common features among the three high-energy-barrier steps responsible for the low catalytic efficiency were revealed using quantum mechanics (QM). Five key residues were simultaneously tailored to stabilize the rate-limiting transition state with the aid of the Rosetta design. The 14 top-ranked variants showed 16.9-143-fold improved catalytic activity. The catalytic efficiency of the best variant, M9 (Q25F/M60W/W64F/I266A), was significantly increased, with a 1660-fold increase in k_cat /K_m and a 1.5-26.8-fold increase in turnover number (TON) towards various indanone derivatives.

Computation-Aided Engineering of Cytochrome P450 for the Production of Pravastatin

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.

Directed evolution of cytochrome P450DA hydroxylase activity for stereoselective biohydroxylation

A colorimetric high throughput screening method was developed based on a dual-enzyme cascade and used for the directed evolution of cytochrome P450 hydroxylase activity.

Computational Prediction of ω-Transaminase Specificity by a Combination of Docking and Molecular Dynamics Simulations

ω-Transaminases (ω-TAs) catalyze the conversion of ketones to chiral amines, often with high enantioselectivity and specificity, which makes them attractive for industrial production of chiral amines. Tailoring ω-TAs to accept non-natural substrates is necessary because of their limited substrate range. We present a computational protocol for predicting the enantioselectivity and catalytic selectivity of an ω-TA from Vibrio fluvialis with different substrates and benchmark it against 62 compounds gathered from the literature. Rosetta-generated complexes containing an external aldimine intermediate of the transamination reaction are used as starting conformations for multiple short independent molecular dynamics (MD) simulations. The combination of molecular docking and MD simulations ensures sufficient and accurate sampling of the relevant conformational space. Based on the frequency of near-attack conformations observed during the MD trajectories, enantioselectivities can be quantitatively predicted. The predicted enantioselectivities are in agreement with a benchmark dataset of experimentally determined ee% values. The substrate-range predictions can be based on the docking score of the external aldimine intermediate. The low computational cost required to run the presented framework makes it feasible for use in enzyme design to screen thousands of enzyme variants.

Conformational Motion of Ferredoxin Enables Efficient Electron Transfer to Heme in the Full-Length P450TT

Cytochrome P450 monooxygenases (P450s) are versatile biocatalysts used in natural products biosynthesis, xenobiotic metabolisms, and biotechnologies. In P450s, the electrons required for O₂ activation are supplied by NAD(P)H through stepwise electron transfers (ETs) mediated by redox partners. While much is known about the machinery of the catalytic cycle of P450s, the mechanisms of long-range ET are largely unknown. Very recently, the first crystal structure of full-length P450_TT was solved. This enables us to decipher the interdomain ET mechanism between the [2Fe-2S]-containing ferredoxin and the heme, by use of molecular dynamics simulations. In contrast to the "distal" conformation characterized in the crystal structure where the [2Fe-2S] cluster is ∼28 Å away from heme-Fe, our simulations demonstrated a "proximal" conformation of [2Fe-2S] that is ∼17 Å [and 13.7 Å edge-to-edge] away from heme-Fe, which may enable the interdomain ET. Key residues involved in ET pathways and interdomain complexation were identified, some of which have already been verified by recent mutation studies. The conformational transit of ferredoxin between "distal" and "proximal" was found to be controlled mostly by the long-range electrostatic interactions between the ferredoxin domain and the other two domains. Furthermore, our simulations show that the full-length P450_TT utilizes a flexible ET pathway that resembles either P450_Scc or P450_cam. Thus, this study provides a uniform picture of the ET process between reductase domains and heme domain.

Structural insight into the electron transfer pathway of a self-sufficient P450 monooxygenase

Cytochrome P450 monooxygenases are versatile heme-thiolate enzymes that catalyze a wide range of reactions. Self-sufficient cytochrome P450 enzymes contain the redox partners in a single polypeptide chain. Here, we present the crystal structure of full-length CYP116B46, a self-sufficient P450. The continuous polypeptide chain comprises three functional domains, which align well with the direction of electrons traveling from FMN to the heme through the [2Fe-2S] cluster. FMN and the [2Fe-2S] cluster are positioned closely, which facilitates efficient electron shuttling. The edge-to-edge straight-line distance between the [2Fe-2S] cluster and heme is approx. 25.3 Å. The role of several residues located between the [2Fe-2S] cluster and heme in the catalytic reaction is probed in mutagenesis experiments. These findings not only provide insights into the intramolecular electron transfer of self-sufficient P450s, but are also of interest for biotechnological applications of self-sufficient P450s.

Self-sufficient cytochrome P450 monooxygenases, which contain all redox partners in a single polypeptide chain, are of interest for biotechnological applications. Here, the authors present the crystal structure of full-length Thermobispora bispora CYP116B46 and discuss the potential electron transfer pathway.

Characterization of a Self-Sufficient Cytochrome P450 Monooxygenase from Deinococcus apachensis for Enantioselective Benzylic Hydroxylation

A self-sufficient cytochrome P450 monooxygenase from Deinococcus apachensis (P450DA) was identified and successfully overexpressed in Escherichia coli BL21(DE3). P450DA would be a member of the CYP102D subfamily and assigned as CYP102D2 according to the phylogenetic tree and sequence alignment. Purification and characterization of the recombinant P450DA indicated both NADH and NADPH could be used by P450DA as a reducing cofactor. The recombinant E. coli (P450DA) strain was functionally active, showing excellent enantioselectivity for benzylic hydroxylation of methyl 2-phenylacetate. Further substrate scope studies revealed that P450DA is able to catalyze benzylic hydroxylation of a variety of compounds, affording the corresponding chiral benzylic alcohols in 86-99 % ee and 130-1020 total turnover numbers.

Chemoenzymatic o-Quinone Methide Formation

Generation of reactive intermediates and interception of these fleeting species under physiological conditions is a common strategy employed by Nature to build molecular complexity. However, selective formation of these species under mild conditions using classical synthetic techniques is an outstanding challenge. Here, we demonstrate the utility of biocatalysis in generating o-quinone methide intermediates with precise chemoselectivity under mild, aqueous conditions. Specifically, α-ketoglutarate-dependent non-heme iron enzymes, CitB and ClaD, are employed to selectively modify benzylic C-H bonds of o-cresol substrates. In this transformation, biocatalytic hydroxylation of a benzylic C-H bond affords a benzylic alcohol product which, under the aqueous reaction conditions, is in equilibrium with the corresponding o-quinone methide. o-Quinone methide interception by a nucleophile or a dienophile allows for one-pot conversion of benzylic C-H bonds into C-C, C-N, C-O, and C-S bonds in chemoenzymatic cascades on preparative scale. The chemoselectivity and mild nature of this platform is showcased here by the selective modification of peptides and chemoenzymatic synthesis of the chroman natural product (-)-xyloketal D.

Arylation Chemistry for Bioconjugation

Bioconjugation chemistry has been used to prepare modified biomolecules with functions beyond what nature intended. Central to these techniques is the development of highly efficient and selective bioconjugation reactions that operate under mild, biomolecule compatible conditions. Methods that form a nucleophile-sp2 carbon bond show promise for creating bioconjugates with new modifications, sometimes resulting in molecules with unparalleled functions. Here we outline and review sulfur, nitrogen, selenium, oxygen, and carbon arylative bioconjugation strategies and their applications to modify peptides, proteins, sugars, and nucleic acids.

A Structural Model of a P450-Ferredoxin Complex from Orientation-Selective Double Electron-Electron Resonance Spectroscopy

Cytochrome P450 (CYP) monooxygenases catalyze the oxidation of chemically inert carbon-hydrogen bonds in diverse endogenous and exogenous organic compounds by atmospheric oxygen. This C-H bond oxy-functionalization activity has huge potential in biotechnological applications. Class I CYPs receive the two electrons required for oxygen activation from NAD(P)H via a ferredoxin reductase and ferredoxin. The interaction of Class I CYPs with their cognate ferredoxin is specific. In order to reconstitute the activity of diverse CYPs, structural characterization of CYP-ferredoxin complexes is necessary, but little structural information is available. Here we report a structural model of such a complex (CYP199A2-HaPux) in frozen solution derived from distance and orientation restraints gathered by the EPR technique of orientation-selective double electron-electron resonance (os-DEER). The long-lived oscillations in the os-DEER spectra were well modeled by a single orientation of the CYP199A2-HaPux complex. The structure is different from the two known Class I CYP-Fdx structures: CYP11A1-Adx and CYP101A1-Pdx. At the protein interface, HaPux residues in the [Fe₂S₂] cluster-binding loop and the α3 helix and the C-terminus residue interact with CYP199A2 residues in the proximal loop and the C helix. These residue contacts are consistent with biochemical data on CYP199A2-ferredoxin binding and electron transfer. Electron-tunneling calculations indicate an efficient electron-transfer pathway from the [Fe₂S₂] cluster to the heme. This new structural model of a CYP-Fdx complex provides the basis for tailoring CYP enzymes for which the cognate ferredoxin is not known, to accept electrons from HaPux and display monooxygenase activity.

Stereoselective hydroxylation of isophorone by variants of the cytochromes P450 CYP102A1 and CYP101A1

The stereoselective oxidation of hydrocarbons is an area of research where enzyme biocatalysis can make a substantial impact. The cyclic ketone isophorone was stereoselectively hydroxylated (≥95%) by wild-type CYP102A1 to form (R)-4-hydroxyisophorone, an important chiral synthon and flavour and fragrance compound. CYP102A1 variants were also selective for 4-hydroxyisophorone formation and the product formation rate increased over the wild-type enzyme by up to 285-fold, with the best mutants being R47L/Y51F/I401P and A74G/F87V/L188Q. The latter variant, which contained mutations in the distal substrate binding pocket, was marginally less selective. Combining perfluorodecanoic acid decoy molecules with the rate accelerating variant R47L/Y51F/I401P engendered further improvement with the purified enzymes. However when the decoy molecules were used with A74G/F87V/L188Q the amount of product generated by the enzyme was reduced. Addition of decoy molecules to whole-cell turnovers did not improve the productivity of these CYP102A1 systems. WT CYP101A1 formed significant levels of 7-hydroxyisophorone as a minor product alongside 4-hydroxyisophorone. However the F87W/Y96F/L244A/V247L CYP101A1 mutant was ≥98% selective for (R)-4-hydroxyisophorone. A comparison of the two enzyme systems using whole-cell oxidation reactions showed that the best CYP101A1 variant was able to generate more product. We also characterised that the further oxidation metabolite 4-ketoisophorone was produced and then subsequently reduced to levodione by an endogenous Escherichia coli ene reductase.

Recent advances in enzymatic oxidation of alcohols

Enzymatic alcohol oxidation plays an important role in chemical synthesis. In the past two years, new alcohol oxidation enzymes were developed through genome-mining and protein engineering, such as new copper radical oxidases with broad substrate scope, alcohol dehydrogenases with altered cofactor preference and a flavin-dependent alcohol oxidase with enhanced oxygen coupling. New cofactor recycling methods were reported for alcohol dehydrogenase-catalyzed oxidation with photocatalyst and coupled glutaredoxin-glutathione reductase as promising examples. Different alcohol oxidation systems were used for the oxidation of primary and secondary alcohols, especially in the cascade conversion of alcohols to lactones, lactams, chiral amines, chiral alcohols and hydroxyketones. Among them, biocatalyst with low enantioselectivity demonstrated an interesting feature for complete conversion of racemic secondary alcohols through non-enantioselective oxidation.

P450 BM3-Catalyzed Regio- and Stereoselective Hydroxylation Aiming at the Synthesis of Phthalides and Isocoumarins

Cytochrome P450 BM3 monooxygenases are able to catalyze the regio- and stereoselective oxygenation of a broad range of substrates, with promising potential for synthetic applications. To study the suitability of P450 BM3 variants for stereoselective benzylic hydroxylation of 2-alkylated benzoic acid esters, the biotransformation of methyl 2-ethylbenzoate, resulting in both enantiomeric forms of 3-methylphthalide, was investigated. In the case of methyl 2-propylbenzoate as a substrate the regioselectivity of the reaction was shifted towards β-hydroxylation, resulting in the synthesis of enantioenriched R- and S-configured 3-methylisochroman-1-one. The potential of P450 BM3 variants for regio- and stereoselective synthesis of phthalides and isocoumarins offers a new route to a class of compounds that are valuable synthons for a variety of natural compounds.

Enzymatic site-selectivity enabled by structure-guided directed evolution

This review covers recent advances in the directed evolution of enzymes for controlling site-selectivity of hydroxylation, amination and chlorination.

Structural basis of steroid binding and oxidation by the cytochrome P450 CYP109E1 from Bacillus megaterium

Cytochrome P450 monooxygenases (P450s) are attractive enzymes for the pharmaceutical industry, in particular, for applications in steroidal drug synthesis. Here, we report a comprehensive functional and structural characterization of CYP109E1, a novel steroid‐converting cytochrome P450 enzyme identified from the genome of Bacillus megaterium DSM319. In vitro and whole‐cell in vivo turnover experiments, combined with binding assays, revealed that CYP109E1 is able to hydroxylate testosterone at position 16β. Related steroids with bulky substituents at carbon C17, like corticosterone, bind to the enzyme without being converted. High‐resolution X‐ray structures were solved of a steroid‐free form of CYP109E1 and of complexes with testosterone and corticosterone. The structural analysis revealed a highly dynamic active site at the distal side of the heme, which is wide open in the absence of steroids, can bind four ordered corticosterone molecules simultaneously, and undergoes substantial narrowing upon binding of single steroid molecules. In the crystal structures, the single bound steroids adopt unproductive binding modes coordinating the heme‐iron with their C3‐keto oxygen. Molecular dynamics (MD) simulations suggest that the steroids may also bind in ~180° reversed orientations with the C16 carbon and C17‐substituents pointing toward the heme, leading to productive binding of testosterone explaining the observed regio‐ and stereoselectivity. The X‐ray structures and MD simulations further identify several residues with important roles in steroid binding and conversion, which could be confirmed by site‐directed mutagenesis. Taken together, our results provide unique insights into the CYP109E1 activity, substrate specificity, and regio/stereoselectivity.

Database

The atomic coordinates and structure factors have been deposited in the Protein Data Bank with accession codes 5L90 (steroid‐free CYP109E1), 5L91 (CYP109E1‐COR4), 5L94 (CYP109E1‐TES), and 5L92 (CYP109E1‐COR).

Enzymes

Cytochrome P450 monooxygenase CYP109E1, EC 1.14.14.1, UniProt ID: D5DKI8, Adrenodoxin reductase EC 1.18.1.6.