Choose Your Own Adventure: A Comprehensive Database of Reactions Catalyzed by Cytochrome P450 BM3 Variants

Cytochrome P450 BM3 monooxygenase is the topic of extensive research as many researchers have evolved this enzyme to generate a variety of products. However, the abundance of information on increasingly diversified variants of P450 BM3 that catalyze a broad array of chemistry is not in a format that enables easy extraction and interpretation. We present a database that categorizes variants by their catalyzed reactions and includes details about substrates to provide reaction context. This database of >1500 P450 BM3 variants is downloadable and machine-readable and includes instructions to maximize ease of gathering information. The database allows rapid identification of commonly reported substitutions, aiding researchers who are unfamiliar with the enzyme in identifying starting points for enzyme engineering. For those actively engaged in engineering P450 BM3, the database, along with this review, provides a powerful and user-friendly platform to understand, predict, and identify the attributes of P450 BM3 variants, encouraging the further engineering of this enzyme.

Production of Mono-Hydroxylated Derivatives of Terpinen-4-ol by Bacterial CYP102A1 Enzymes

CYP102A1 from Bacillus megaterium is an important enzyme in biotechnology, because engineered CYP102A1 enzymes can react with diverse substrates and produce human cytochrome P450-like metabolites. Therefore, CYP102A1 can be applied to drug metabolite production. Terpinen-4-ol is a cyclic monoterpene and the primary component of essential tea tree oil. Terpinen-4-ol was known for therapeutic effects, including antibacterial, antifungal, antiviral, and anti-inflammatory. Because terpenes are natural compounds, examining novel terpenes and investigating the therapeutic effects of terpenes represent responses to social demands for eco-friendly compounds. In this study, we investigated the catalytic activity of engineered CYP102A1 on terpinen-4-ol. Among CYP102A1 mutants tested here, the R47L/F81I/F87V/E143G/L188Q/N213S/E267V mutant showed the highest activity to terpinen-4-ol. Two major metabolites of terpinen-4-ol were generated by engineered CYP102A1. Characterization of major metabolites was confirmed by liquid chromatography-mass spectrometry (LC-MS), gas chromatography-MS, and nuclear magnetic resonance spectroscopy (NMR). Based on the LC-MS results, the difference in mass-to-charge ratio of an ion (m/z) between terpinen-4-ol and its major metabolites was 16. One major metabolite was defined as 1,4-dihydroxy-p-menth-2-ene by NMR. Given these results, we speculate that another major metabolite is also a mono-hydroxylated product. Taken together, we suggest that CYP102A1 can be applied to make novel terpene derivatives.

Prokaryotic expression and characterization of artificial self-sufficient CYP120A monooxygenases

Cytochrome P450 monooxygenases CYP120As are the unique non-membrane P450s, which are extensively involved in retinoid biodegradation. As the O-functionalized 1,3,3-trimethylcyclohex-1-ene moiety exists in many bioactive compounds which could only be catalyzed by Class II P450s, exploration of the catalytic repertoire of CYP120As is therefore highly attractive. However, up to date, only one bacteriogenic candidate (CYP120A1) was demonstrated for the hydroxylation of C16 and C17 of retinoic acid, by utilizing the integral membrane protein cytochrome P450 reductase redox partner for the electron transfer. Herein, we provided an efficient prokaryotic functional expression system of CYP120As in E. coli by expression of the CYP120A1 coupled with several reductase partners. Fusion redox partners to the C-terminal of the heme-domain are also working on other CYP120A members. Among them, the fusion protein of CYP120A29 and FAD/FMN reductase from Bacillus megaterium P450BM3 (CYP101A2) showed the highest expression level. Based on the available translational fusion systems, the regioselectivity and the substrate scope of the CYP120As have also been explored. This work represents a good starting point for further expanding the catalytic potential of CYP120 family. KEY POINTS: • Characterization of CYP120As in E. coli is firstly achieved by constructing fusion proteins. • The feasibility of three P450 reductase domains to CYP120As was evaluated. • Hydroxylated products of retinoic acid by six CYP120As were sorted and analyzed.

MicroGelzymes: pH-Independent Immobilization of Cytochrome P450 BM3 in Microgels

Microgels are an emerging class of "ideal" enzyme carriers because of their chemical and process stability, biocompatibility, and high enzyme loading capability. In this work, we synthesized a new type of permanently positively charged poly(N-vinylcaprolactam) (PVCL) microgel with 1-vinyl-3-methylimidazolium (quaternization of nitrogen by methylation of N-vinylimidazole moieties) as a comonomer (PVCL/VimQ) through precipitation polymerization. The PVCL/VimQ microgels were characterized with respect to their size, charge, swelling degree, and temperature responsiveness in aqueous solutions. P450 monooxygenases are usually challenging to immobilize, and often, high activity losses occur after the immobilization (in the case of P450 BM3 from Bacillus megaterium up to 100% loss of activity). The electrostatic immobilization of P450 BM3 in permanently positively charged PVCL/VimQ microgels was achieved without the loss of catalytic activity at the pH optimum of P450 BM3 (pH 8; ∼9.4 nmol 7-hydroxy-3-carboxy coumarin ethyl ester/min for free and immobilized P450 BM3); the resulting P450-microgel systems were termed P450 MicroGelzymes (P450 μ-Gelzymes). In addition, P450 μ-Gelzymes offer the possibility of reversible ionic strength-triggered release and re-entrapment of the biocatalyst in processes (e.g., for catalyst reuse). Finally, a characterization of the potential of P450 μ-Gelzymes to provide resistance against cosolvents (acetonitrile, dimethyl sulfoxide, and 2-propanol) was performed to evaluate the biocatalytic application potential.

Metabolic Engineering Strategies for Sustainable Terpenoid Flavor and Fragrance Synthesis

Terpenoids derived from plant material are widely applied in the flavor and fragrance industry. Traditional extraction methods are unsustainable, but microbial synthesis offers a promising solution to attain efficient production of natural-identical terpenoids. Overproduction of terpenoids in microbes requires careful balancing of the synthesis pathway constituents within the constraints of host cell metabolism. Advances in metabolic engineering have greatly facilitated overcoming the challenges of achieving high titers, rates, and yields (TRYs). The review summarizes recent development in the molecular biology toolbox to achieve high TRYs for terpenoid biosynthesis, mainly in the two industrial platform microorganisms: Escherichia coli and Saccharomyces cerevisiae. The biosynthetic pathways, including alternative pathway designs, are briefly introduced, followed by recently developed methodologies used for pathway, genome, and strain optimization. Integrated applications of these tools are important to achieve high "TRYs" of terpenoid production and pave the way for translating laboratory research into successful commercial manufacturing.

A 96-multiplex capillary electrophoresis screening platform for product based evolution of P450 BM3

The main challenge that prevents a broader application of directed enzyme evolution is the lack of high-throughput screening systems with universal product analytics. Most directed evolution campaigns employ screening systems based on colorimetric or fluorogenic surrogate substrates or universal quantification methods such as nuclear magnetic resonance spectroscopy or mass spectrometry, which have not been advanced to achieve a high-throughput. Capillary electrophoresis with a universal UV-based product detection is a promising analytical tool to quantify product formation. Usage of a multiplex system allows the simultaneous measurement with 96 capillaries. A 96-multiplexed capillary electrophoresis (MP-CE) enables a throughput that is comparable to traditional direct evolution campaigns employing 96-well microtiter plates. Here, we report for the first time the usage of a MP-CE system for directed P450 BM3 evolution towards increased product formation (oxidation of alpha-isophorone to 4-hydroxy-isophorone; highest reached total turnover number after evolution campaign: 7120 mol_4-OH mol_P450⁻¹). The MP-CE platform was 3.5-fold more efficient in identification of beneficial variants than the standard cofactor (NADPH) screening system.

Exploring the substrate specificity of Cytochrome P450cin

Cytochromes P450 are enzymes that catalyse the oxidation of a wide variety of compounds that range from small volatile compounds, such as monoterpenes to larger compounds like steroids. These enzymes can be modified to selectively oxidise substrates of interest, thereby making them attractive for applications in the biotechnology industry. In this study, we screened a small library of terpenes and terpenoid compounds against P450_cin and two P450_cin mutants, N242A and N242T, that have previously been shown to affect selectivity. Initial screening indicated that P450_cin could catalyse the oxidation of most of the monoterpenes tested; however, sesquiterpenes were not substrates for this enzyme or the N242A mutant. Additionally, both P450_cin mutants were found to be able to oxidise other bicyclic monoterpenes. For example, the oxidation of (R)- and (S)-camphor by N242T favoured the production of 5-endo-hydroxycamphor (65-77% of the total products, dependent on the enantiomer), which was similar to that previously observed for (R)-camphor with N242A (73%). Selectivity was also observed for both (R)- and (S)-limonene where N242A predominantly produced the cis-limonene 1,2-epoxide (80% of the products following (R)-limonene oxidation) as compared to P450_cin (23% of the total products with (R)-limonene). Of the three enzymes screened, only P450_cin was observed to catalyse the oxidation of the aromatic terpene p-cymene. All six possible hydroxylation products were generated from an in vivo expression system catalysing the oxidation of p-cymene and were assigned based on ¹H NMR and GC-MS fragmentation patterns. Overall, these results have provided the foundation for pursuing new P450_cin mutants that can selectively oxidise various monoterpenes for biocatalytic applications.

Structural analysis of Cytochrome P450 BM3 mutant M11 in complex with dithiothreitol

The bacterial Cytochrome P450 (CYP) BM3 (CYP102A1) is one of the most active CYP isoforms. BM3 mutants can serve as a model for human drug-metabolizing CYPs and/or as biocatalyst for selective formation of drug metabolites. Hence, molecular and computational biologists have in the last two decades shown strong interest in the discovery and design of novel BM3 variants with optimized activity and selectivity for substrate conversion. This led e.g. to the discovery of mutant M11 that is able to metabolize a variety of drugs and drug-like compounds with relatively high activity. In order to further improve our understanding of CYP binding and reactions, we performed a co-crystallization study of mutant M11 and report here the three-dimensional structure M11 in complex with dithiothreitol (DTT) at a resolution of 2.16 Å. The structure shows that DTT can coordinate to the Fe atom in the heme group. UV/Vis spectroscopy and molecular dynamics simulation studies underline this finding and as first structure of the CYP BM3 mutant M11 in complex with a ligand, it offers a basis for structure-based design of novel mutants.

Characterization and structure-guided engineering of the novel versatile terpene monooxygenase CYP109Q5 from Chondromyces apiculatus DSM436

One of the major challenges in chemical synthesis is the selective oxyfunctionalization of non-activated C-H bonds, which can be enabled by biocatalysis using cytochrome P450 monooxygenases. In this study, we report on the characterization of the versatile CYP109Q5 from Chondromyces apiculatus DSM436, which is able to functionalize a wide range of substrates (terpenes, steroids and drugs), including the ring of β-ionone in non-allylic positions. The crystal structure of CYP109Q5 revealed flexibility within the active site pocket that permitted the accommodation of bulky substrates, and enabled a structure-guided approach to engineering the enzyme. Some variants of CYP109Q5 displayed a switch in selectivity towards the non-allylic positions of β-ionone, allowing the simultaneous production of 2- and 3-hydroxy-β-ionone, which are chemically challenging to synthesize and are important precursors for carotenoid synthesis. An efficient whole-cell system finally enabled the production of up to 0.5 g l-1 hydroxylated products of β-ionone; this system can be applied to product identification in further biotransformations. Overall, CYP109Q5 proved to be highly evolvable and active. The studies in this work demonstrate that, using rational mutagenesis, the highly versatile CYP109Q5 generalist can be progressively evolved to be an industrially valuable specialist for the synthesis of specific products.

Fungi-Mediated Biotransformation of the Isomeric Forms of the Apocarotenoids Ionone, Damascone and Theaspirane

In this work, we describe a study on the biotransformation of seven natural occurring apocarotenoids by means of eleven selected fungal species. The substrates, namely ionone (α-, β- and γ-isomers), 3,4-dehydroionone, damascone (α- and β-isomers) and theaspirane are relevant flavour and fragrances components. We found that most of the investigated biotransformation reactions afforded oxidized products such as hydroxy- keto- or epoxy-derivatives. On the contrary, the reduction of the keto groups or the reduction of the double bond functional groups were observed only for few substrates, where the reduced products are however formed in minor amount. When starting apocarotenoids are isomers of the same chemical compound (e.g., ionone isomers) their biotransformation can give products very different from each other, depending both on the starting substrate and on the fungal species used. Since the majority of the starting apocarotenoids are often available in natural form and the described products are natural compounds, identified in flavours or fragrances, our biotransformation procedures can be regarded as prospective processes for the preparation of high value olfactory active compounds.

Process development for oxidations of hydrophobic compounds applying cytochrome P450 monooxygenases in-vitro

Cytochrome P450 monooxygenases are a unique family of enzymes that are able to catalyze regio- and stereospecific oxidations for a broad substrate range. However, due to limited enzyme activities and stabilities, hydrophobicity of substrates, as well as the necessity of a continuous electron and oxygen supply the implementation of P450s for industrial processes remains challenging. Aim of this study was to point out key aspects for the development of an efficient synthesis concept for cytochrome P450 catalyzed oxidations. In order to regenerate the natural cofactor NADPH, a glucose dehydrogenase was applied. The low water soluble terpene α-ionone was used as substrate for the model reaction system. The studies reveal that an addition of surfactants in combination with low volumetric amounts of co-solvent can significantly increase substrate availability and reaction rates. Furthermore, these additives facilitated a reliable sampling procedure during the process. Another key factor for the process design was the oxygen supply. Based on various investigations, a bubble-aerated stirred tank reactor in batch mode represents a promising reactor concept for P450 oxidations. Main restriction of the investigated reaction system was the low process stability of the P450 monooxygenase, characterized by maximum total turnover numbers of ∼4100molα-ionone/molP450.

Insights into regioselective metabolism of mefenamic acid by cytochrome P450 BM3 mutants through crystallography, docking, molecular dynamics, and free energy calculations

Cytochrome P450 BM3 (CYP102A1) mutant M11 is able to metabolize a wide range of drugs and drug-like compounds. Among these, M11 was recently found to be able to catalyze formation of human metabolites of mefenamic acid and other nonsteroidal anti-inflammatory drugs (NSAIDs). Interestingly, single active-site mutations such as V87I were reported to invert regioselectivity in NSAID hydroxylation. In this work, we combine crystallography and molecular simulation to study the effect of single mutations on binding and regioselective metabolism of mefenamic acid by M11 mutants. The heme domain of the protein mutant M11 was expressed, purified, and crystallized, and its X-ray structure was used as template for modeling. A multistep approach was used that combines molecular docking, molecular dynamics (MD) simulation, and binding free-energy calculations to address protein flexibility. In this way, preferred binding modes that are consistent with oxidation at the experimentally observed sites of metabolism (SOMs) were identified. Whereas docking could not be used to retrospectively predict experimental trends in regioselectivity, we were able to rank binding modes in line with the preferred SOMs of mefenamic acid by M11 and its mutants by including protein flexibility and dynamics in free-energy computation. In addition, we could obtain structural insights into the change in regioselectivity of mefenamic acid hydroxylation due to single active-site mutations. Our findings confirm that use of MD and binding free-energy calculation is useful for studying biocatalysis in those cases in which enzyme binding is a critical event in determining the selective metabolism of a substrate.

Selective oxidation of carotenoid-derived aroma compounds by CYP260B1 and CYP267B1 from Sorangium cellulosum So ce56

Due to their bioactive properties as well as their application as precursors in chemical synthesis, hydroxylated isoprenoids and norisoprenoids are very valuable compounds. The efficient hydroxylation of such compounds remains a challenge in organic chemistry caused by the formation of a variety of side products and lack of overall regio- and stereoselectivity. In contrast, cytochromes P450 are known for their selective oxidation under mild conditions. Here, we demonstrate for the first time the ability of myxobacterial CYP260B1 and CYP267B1 from Sorangium cellulosum So ce56 to oxidize such carotenoid-derived aroma compounds. A focused library of 14 substrates such as ionones, damascones, as well as some of their isomers and derivatives was screened in vitro. Both P450s were capable of an efficient oxidation of all tested compounds. CYP260B1-dependent conversions mainly formed multiple products, whereas conversions by CYP267B1 resulted predominantly in a single product. To identify the main products by NMR spectroscopy, an Escherichia coli-based whole-cell system was used. CYP267B1 showed a hydroxylase activity towards the formation of allylic alcohols. Likewise, CYP260B1 performed the allylic hydroxylation of β-damascone [(E)-1-(2,6,6-trimethylcyclohex-1-enyl)but-2-en-1-one] and δ-damascone [(E)-1-(2,6,6-trimethylcyclohex-3-enyl)but-2-en-1-one]. Moreover, CYP260B1 showed an epoxidase activity towards β-ionone [(E)-4-(2,6,6-trimethylcyclohex-1-enyl)but-3-en-2-one] as well as the methyl-substituted α-ionone derivatives raldeine [(E)-1-(2,6,6-trimethylcyclohex-2-enyl)pent-1-en-3-one] and isoraldeine [(E)-3-methyl-4-(2,6,6-trimethylcyclohex-2-enyl)but-3-en-2-one]. In addition, to known products, also novel products such as 2-OH-δ-damascone [(E)-1-(5-hydroxy-2,6,6-trimethylcyclohex-3-enyl)but-2-en-1-one], 3-OH-allyl-α-ionone [(E)-1-(4-hydroxy-2,6,6-trimethylcyclohex-2-enyl)hepta-1,6-dien-3-one], and 4-OH-allyl-β-ionone [(E)-1-(3-hydroxy-2,6,6-trimethylcyclohex-1-enyl)hepta-1,6-dien-3-one] were identified during our studies.

Application of a cocktail approach to screen cytochrome P450 BM3 libraries for metabolic activity and diversity

In the present study, the validity of using a cocktail screening method in combination with a chemometrical data mining approach to evaluate metabolic activity and diversity of drug-metabolizing bacterial Cytochrome P450 (CYP) BM3 mutants was investigated. In addition, the concept of utilizing an in-house-developed library of CYP BM3 mutants as a unique biocatalytic synthetic tool to support medicinal chemistry was evaluated. Metabolic efficiency of the mutant library towards a selection of CYP model substrates, being amitriptyline (AMI), buspirone (BUS), coumarine (COU), dextromethorphan (DEX), diclofenac (DIC) and norethisterone (NET), was investigated. First, metabolic activity of a selection of CYP BM3 mutants was screened against AMI and BUS. Subsequently, for a single CYP BM3 mutant, the effect of co-administration of multiple drugs on the metabolic activity and diversity towards AMI and BUS was investigated. Finally, a cocktail of AMI, BUS, COU, DEX, DIC and NET was screened against the whole in-house CYP BM3 library. Different validated quantitative and qualitative (U)HPLC-MS/MS-based analytical methods were applied to screen for substrate depletion and targeted product formation, followed by a more in-depth screen for metabolic diversity. A chemometrical approach was used to mine all data to search for unique metabolic properties of the mutants and allow classification of the mutants. The latter would open the possibility of obtaining a more in-depth mechanistic understanding of the metabolites. The presented method is the first MS-based method to screen CYP BM3 mutant libraries for diversity in combination with a chemometrical approach to interpret results and visualize differences between the tested mutants.

The crystal structure of the versatile cytochrome P450 enzyme CYP109B1 from Bacillus subtilis

The crystal structure of the versatile CYP109B1 enzyme from Bacillus subtilis has been solved at 1.8 Å resolution. This is the first structure of an enzyme from this CYP family, whose members are prevalent across diverse species of bacteria. In the crystal structure the enzyme has an open conformation with an access channel leading from the heme to the surface. The substrate-free structure reveals the location of the key residues in the active site that are responsible for binding the substrate in the correct orientation for regioselective oxidation. Importantly, there are significant differences among these residues in members of the CYP109 and closely related CYP106 families and these likely account for the variations in substrate binding and oxidation profiles observed with these enzymes. A whole-cell oxidation biosystem was developed, which contains CYP109B1 and a phthalate family oxygenase reductase (PFOR), from Pseudomonas putida KT24440, as the electron transfer partner. This electron transfer system is able to support CYP109B1 activity resulting in the regioselective hydroxylation of both α- and β-ionone in vivo and in vitro. The PFOR is therefore a versatile electron transfer partner that is able to support the activity of CYP enzymes from other bacterium. The crystal structure of CYP109B1 has a positively charged proximal face and this explains why it can interact with PFOR and adrenodoxin which are predominantly negatively charged around their [2Fe-2S] clusters.

The efficient and selective biocatalytic oxidation of norisoprenoid and aromatic substrates by CYP101B1 from Novosphingobium aromaticivorans DSM12444

CYP101B1 fromNovosphingobium aromaticivoransoxidises ionone derivatives and phenylcyclohexane with high activity and regioselectivity.

Expanding the toolbox of organic chemists: directed evolution of P450 monooxygenases as catalysts in regio- and stereoselective oxidative hydroxylation

Cytochrome P450 enzymes (CYPs) have been used for more than six decades as catalysts for the CH-activating oxidative hydroxylation of organic compounds with formation of added-value products.

Selective allylic hydroxylation of acyclic terpenoids by CYP154E1 from Thermobifida fusca YX

Allylic alcohols are valuable precursors in the synthesis of pharmaceutical intermediates, agrochemicals and natural products. Regioselective oxidation of parental alkenes is a challenging task for chemical catalysts and requires several steps including protection and deprotection. Many cytochrome P450 enzymes are known to catalyse selective allylic hydroxylation under mild conditions. Here, we describe CYP154E1 from Thermobifida fusca YX that enables this type of oxidation. Several acyclic terpenoids were tested as possible substrates for CYP154E1, and the regio- and chemoselectivity of their oxidation was investigated. Using a previously established bioinformatics approach we identified position 286 in the active site of CYP154E1 which is putatively involved in substrate binding and thereby might have an effect on enzyme selectivity. To tune regio- and chemoselectivity of the enzyme three mutants at position 286 were constructed and used for substrate oxidation. All formed products were analysed with GC-MS and identified using chemically synthesised authentic samples and known compounds as references. Best regioselectivity towards geraniol and nerol was observed with the wild type enzyme mainly leading to 8-hydroxy derivatives (8-hydroxygeraniol or 8-hydroxynerol) with high selectivity (100% and 96% respectively). Highest selectivities during the oxidation of geranylacetone and nerylacetone were observed with the following variants: V286F led mainly to 7-hydroxygeranylacetone (60% of the total product) and V286A produced predominantly 12-hydroxynerylacetone (75% of total product). Thus, CYP154E1 and its mutants expand the tool-box for allylic hydroxylation in synthetic chemistry.

Selective aliphatic carbon–hydrogen bond activation of protected alcohol substrates by cytochrome P450 enzymes

Protected cyclohexanol and cyclohex-2-enol substrates were efficiently and selectively oxidised by different P450cam mutants providing a general methodology for generating substituted diols using biocatalysts.

Free energy calculations give insight into the stereoselective hydroxylation of α-ionones by engineered cytochrome P450 BM3 mutants

Previously, stereoselective hydroxylation of α-ionone by Cytochrome P450 BM3 mutants M01 A82W and M11 L437N was observed. While both mutants hydroxylate α-ionone in a regioselective manner at the C3 position, M01 A82W catalyzes formation of trans-3-OH-α-ionone products whereas M11 L437N exhibits opposite stereoselectivity, producing trans-(3S,6S)-OH-α-ionone and cis-(3S,6R)-OH-α-ionone. Here, we explore the stereoselective C3 hydroxylation of α-ionone by Cytochrome P450 BM3 mutants M01 A82W and M11 L437N using molecular dynamics-based free energy calculations to study the interaction between the enzyme and both the substrates and the products. The one-step perturbation approach is applied using an optimized reference state for substrates and products. While the free energy differences between the substrates free in solution amount to ~0 kJ mol(-1), the differences in mutant M01 A82W agree with the experimentally obtained dissociation constants K(d). Moreover, a correlation with experimentally observed trends in product formation is found in both mutants. The trans isomers show the most favorable relative binding free energy in the range of all four possible hydroxylated diastereomers for mutant M01 A82W, while the trans product from (6S)-α-ionone and the cis product from (6R)-α-ionone show highest affinity for mutant M11 L437N. Marcus theory is subsequently used to relate the thermodynamic stability to transition state energies and rates of formation.