Cytochrome P450-The Wonderful Nanomachine Revealed through Dynamic Simulations of the Catalytic Cycle

Cytochrome P450 electrochemical biosensors transforming in vitro metabolism testing - Opportunities and challenges

The ability of the living world to flourish in the face of constant exposure to dangerous chemicals depends on the management ability of a widespread group of enzymes known as heme-thiolate monooxygenases or cytochrome P450 superfamily. About three-quarters of all reactions determining the metabolism of endogenous compounds, of those carried in foods, of taken drugs, or even of synthetic chemicals discarded into the environment depend on their catalytic performance. The chromatographic and (photo)luminometric methods routinely used as predictive and analytical tools in laboratories have significant drawbacks ranging from limited shelf-life of reagents, use of synthetic substrates, laborious and tedious procedures for highly sensitive detection. In this review, alternative electrochemical biosensors using the cytochrome P450 enzymes as bio-element are emphasized in their main aspects as well regarding their implementation and usefulness. Despite the various schemes proposed for the implementation, reports on real applications are scant for several reasons, including low reaction rates, broad substrate specificity, uncoupling reactions occurrence, and the need for expensive electron transfer partners to promote electron transfer. Finally, the prospect for future developments is introduced, focusing on integrating miniaturized systems with electrochemical techniques, alongside optimizing enzyme immobilization methods and electrode modifications to improve enzymatic stability and enhance sensor reliability. This progress represents a crucial step towards the creation of portable biosensors that mimic human physiological responses, supporting the precision medicine approach.

Debate of Nucleophilic versus Electrophilic Oxidative Aldehyde Deformylation by Mononuclear Nonheme Iron(III)-Peroxo and Iron(IV)-Oxo Complexes

High-valent iron(IV)-oxo species are fleeting intermediates that perform vital reactions in enzymatic catalysis. In contrast, heme and nonheme iron(III)-peroxo intermediates usually act as nucleophiles and are converted to high-valent iron-oxo intermediates for electrophilic oxidation reactions. Herein, we report a study on aldehyde deformylation reactions of 2-phenylpropionaldehyde (2-PPA) and its derivatives by iron(III)-peroxo complexes bearing tetramethylated cyclam (TMC) analogues, including [FeIII(O2)(12-TMC)]+ (1), [FeIII(O2)(13-TMC)]+ (2), and [FeIII(O2)(14-TMC)]+ (3). Reactivity studies by employing deuterated substrates, such as α-[D1]-2-phenylpropionaldehyde and aldehyde-[D]-2-phenylpropionaldehyde, demonstrate that deformylation of 2-PPA by the nonheme iron(III)-peroxo complexes occurs via abstraction of the stronger aldehyde C-H atom, rather than the expected nucleophilic attack or weaker α-C-H atom abstraction reactions. Interestingly, the preference for aldehyde C-H atom abstraction is retained during the deformylation of 2-PPA by iron(IV)-oxo complexes, i.e., [FeIV(O)(13-TMC)]2+ (4) and [FeIV(O)(N4Py)]2+ (5). DFT calculations reproduce the experimental trends in reactivity and reveal that the peroxide O-O bond is cleaved to form an iron(III)-dioxyl species that conducts aldehyde C-H bond abstraction; this chemoselectivity is achieved through stabilizing noncovalent interactions between the oxidants and the aromatic ring of the substrate that positions the aldehyde in close proximity to the FeIII-O2/FeIV═O cores. These new experimental and theoretical findings together with the previous demonstrations of the ability of 1-3 in hydrogen atom transfer, oxygen atom transfer, and cis-dihydroxylation reactions highlight that iron(III)-peroxo cores are not inherently nucleophiles and can have more important functions in chemical and biological oxidation reactions, rather than acting as transient species en route to high-valent metal-oxo intermediates.

Fusing Engineered CYP109E1 and a New Reductase Domain for 25(OH)VD3 Biosynthesis

25-hydroxyvitamin D3 (25(OH)VD3) is an important daily nutritional supplement, and direct enzymatic C25 hydroxylation of vitamin D3 (VD3) to 25(OH)VD3 is a green and sustainable method. However, the low catalytic activity and electron transfer efficiency of redox partners of cytochrome P450 limited its production. Here, structure-guided semirational design of CYP109E1 led to the mutant CYP109E1M2, which showed a 2-fold increase in 25(OH)VD3 production compared to the wild-type. Furthermore, the novel reductase domain BmRD was first employed to construct a fusion protein with CYP109E1M2. Notably, truncating only two amino acids at the N-terminus of BmRD generated a fusion protein Chimera-2, resulting in a 38.5% increase in 25(OH)VD3 production. Under optimized conditions, the engineered strain Escherichia coli 03-2 produced 491.3 mg/L 25(OH)VD3 with a conversion of 49.0% and a space-time yield of 61.4 mg/(L·h). This work demonstrates the modification and optimization potential of P450 and lays a theoretical foundation for the industrialization of the 25(OH)VD3 biosynthesis.

Brain Cytochrome P450: Navigating Neurological Health and Metabolic Regulation

Human cytochrome P450 (CYP) enzymes in the brain represent a crucial frontier in neuroscience, with far-reaching implications for drug detoxification, cellular metabolism, and the progression of neurodegenerative diseases. The brain's complex architecture, composed of interconnected cell types and receptors, drives unique neuronal signaling pathways, modulates enzyme functions, and leads to distinct CYP gene expression and regulation patterns compared to the liver. Despite their relatively low levels of expression, brain CYPs exert significant influence on drug responses, neurotoxin susceptibility, behavior, and neurological disease risk. These enzymes are essential for maintaining brain homeostasis, mediating cholesterol turnover, and synthesizing and metabolizing neurochemicals, neurosteroids, and neurotransmitters. Moreover, they are key participants in oxidative stress responses, neuroprotection, and the regulation of inflammation. In addition to their roles in metabolizing psychotropic drugs, substances of abuse, and endogenous compounds, brain CYPs impact drug efficacy, safety, and resistance, underscoring their importance beyond traditional drug metabolism. Their involvement in critical physiological processes also links them to neuroprotection, with significant implications for the onset and progression of neurodegenerative diseases. Understanding the roles of cerebral CYP enzymes is vital for advancing neuroprotective strategies, personalizing treatments for brain disorders, and developing CNS-targeting therapeutics. This review explores the emerging roles of CYP enzymes, particularly those within the CYP1-3 and CYP46 families, highlighting their functional diversity and the pathological consequences of their dysregulation on neurological health. It also examines the potential of cerebral CYP-based biomarkers to improve the diagnosis and treatment of neurodegenerative disorders, offering new avenues for therapeutic innovation.

Mechanism of Regio- and Enantioselective Hydroxylation of Arachidonic Acid Catalyzed by Human CYP2E1: A Combined Molecular Dynamics and Quantum Mechanics/Molecular Mechanics Study

Regio- and enantioselective hydroxylation of free fatty acids by human cytochrome P450 2E1 (CYP2E1) plays an important role in metabolic regulation and has significant pathological implications. Despite extensive research, the detailed hydroxylation mechanism of CYP2E1 remains incompletely understood. To clarify the origins of regioselectivity and enantioselectivity observed for CYP2E1-mediated fatty acid hydroxylation, molecular dynamics (MD) simulations and quantum mechanics/molecular mechanics (QM/MM) calculations were performed. MD simulations provided key insights into the proximity of arachidonic acid's carbon atoms to the reactive iron(IV)-oxo moiety in compound I (Cpd I), with the ω-1 position being closest, indicating higher reactivity at this site. QM/MM calculations identified hydrogen abstraction as the rate-determining step, with the ω-1S transition state exhibiting the lowest energy barrier, consistent with experimentally observed enantioselectivity. Energy decomposition analysis revealed that variations in quantum mechanical energy (ΔEQM) significantly influence reaction barriers, with the most efficient hydrogen abstraction occurring at the ω-1S and ω-2R positions. These findings underscore the importance of substrate positioning within the active site in determining product selectivity. Comparisons with two related P450s, P450BM3 and P450SPα, further highlighted the critical role of active site architecture and substrate positioning in modulating selectivity. While surrounding residues do not directly dictate product selectivity, they shape the active site environment and influence substrate positioning. Furthermore, our analysis revealed a previously unrecognized catalytic role of Ala299. These findings provide a deeper mechanistic understanding of human CYP2E1 and offer valuable insights for its precise engineering in targeted C-H functionalization.

Mechanistic Perspective on C-N and C-S Bond Construction Catalyzed by Cytochrome P450 Enzymes

Cytochrome P450 enzymes catalyze a large number of oxidative transformations that are responsible for natural product synthesis. Recent studies have revealed their unique ability to catalyze the formation of C-N and C-S bonds, broadening their biosynthetic applications. However, the enzymatic mechanisms behind these reactions are still unclear. This review focuses on theoretical insights into the mechanisms of P450-catalyzed C-N and C-S bond formation. The key roles of the conformational dynamics of substrate radicals, influenced by the enzyme environment, in modulating selectivity and reactivity are highlighted. Understanding these reaction mechanisms offers valuable guidance for P450 enzyme engineering and the design of biosynthetic applications.

Revealing the Monooxygenase Mechanism for Selective Ring Cleavage of Anthraquinone by BTG13 through Multiscale Simulations

BTG13, a non-heme iron-dependent enzyme with a distinctive coordination environment of four histidines and a carboxylated lysine, has been found to catalyze the cleavage of the C4a-C10 bond in anthraquinone. Contrary to typical dioxygenase mechanisms, our quantum mechanical/molecular mechanical (QM/MM) calculations reveal that BTG13 functions more like a monooxygenase. It selectively inserts an oxygen atom into the C10-C4a bond, creating a lactone species that subsequently undergoes hydrolysis, leading to the formation of a ring-opened product. This discovery highlights the unique catalytic properties of BTG13 and expands our understanding of non-heme iron enzyme mechanisms.

Methylcobalamin protects against liver failure via engaging gasdermin E

Gasdermin E (GSDME) is a pyroptotic cell death effector and a promising target for pyroptotic tissue injury. Here we perform high-throughput screening and demonstrate that methylcobalamin (MeCbl), an endogenous coenzyme form of vitamin B12, is a specific GSDME inhibitor and highly effective against cholestatic liver failure. MeCbl specifically blocks GSDME cleavage by directly binding with GSDME. In cholestasis-, cisplatin- or concanavalin A (Con A)-induced male mouse models, MeCbl significantly suppresses liver transaminase activities and inflammation, alleviates hepatocyte death, and reduces mortality of mice by blocking GSDME cleavage. The conserved Cys180 residue in GSDME is essential for caspase-3/GzmB recognition. MeCbl in base-off conformation coordinates to Cys180 to prevent caspase-3/GzmB-GSDME interactions and thereby GSDME-mediated pyroptosis. In summary, our study discovers MeCbl as a specific GSDME inhibitor that is promisingly to be developed as an effective drug against cholestatic liver failure, and other GSDME triggered sterile inflammation and/or organ failure.

Development of a Ru-porphyrin-based supramolecular framework catalyst for styrene epoxidation

A new microporous supramolecular-framework Ru(II)-porphyrin catalyst containing non-covalent interactions between pyrenylphenyl moieties at the meso-position of the porphyrin ring is synthesised and structurally characterised. This recyclable catalyst expedites styrene epoxidation more efficiently than homogeneous Ru-porphyrin catalytic systems.

Tutorial Review on the Set-Up and Running of Quantum Mechanical Cluster Models for Enzymatic Reaction Mechanisms

Enzymes turnover substrates into products with amazing efficiency and selectivity and as such have great potential for use in biotechnology and pharmaceutical applications. However, details of their catalytic cycles and the origins surrounding the regio- and chemoselectivity of enzymatic reaction processes remain unknown, which makes the engineering of enzymes and their use in biotechnology challenging. Computational modelling can assist experimental work in the field and establish the factors that influence the reaction rates and the product distributions. A popular approach in modelling is the use of quantum mechanical cluster models of enzymes that take the first- and second coordination sphere of the enzyme active site into consideration. These QM cluster models are widely applied but often the results obtained are dependent on model choice and model selection. Herein, we show that QM cluster models can give highly accurate results that reproduce experimental product distributions and free energies of activation within several kcal mol-1, regarded that large cluster models with >300 atoms are used that include key hydrogen bonding interactions and charged residues. In this tutorial review, we give general guidelines on the set-up and applications of the QM cluster method and discuss its accuracy and reproducibility. Finally, several representative QM cluster model examples on metal-containing enzymes are presented, which highlight the strength of the approach.

Computational Study on Flavin-Catalyzed Aerobic Dioxygenation of Alkenyl Thioesters: Decomposition of Anionic Peroxides

Flavin-dependent catalysts are widely applied to aerobic monooxygenation/oxidation reactions. In contrast, flavin-catalyzed aerobic dioxygenation reactions exhibit higher atomic economy but are less reported, not to mention the relevant mechanistic studies. Herein, a density functional theory study on flavin-catalyzed aerobic epoxidation-oxygenolysis of alkenyl thioesters was performed for the first time. Different from the previous mechanistic proposal, a pathway featuring two catalytic stages, monoanionic flavin-C(4a)-peroxide/oxide intermediates, and a reverse reaction sequence (epoxidation goes prior to oxygenolysis) was revealed. In comparison, the pathways involving dianionic flavin catalysts, monoanionic flavin-N(5)-(hydro)peroxide/C(10a)-peroxide, or neutral flavin-C(4a)-hydroperoxide/hydroxide/N(5)-oxide, and the pathways where oxygenolysis goes prior to epoxidation are less favored. Epoxidation goes through intramolecular substitution of the O-O bond of anionic flavin-C(4a)-peroxide by β-carbon, while the resulting flavin-C(4a)-oxide accomplishes the oxygenolysis. Furthermore, two other reaction modes, i.e., concerted O-O cleavage/1,2-shift of α-substituents and dyotropic rearrangement were discovered for the decomposition of other anionic peroxides, and preliminary rules were summarized for understanding the chemoselectivity for this process. This study sheds light on the different reaction features of numerous flavin-dioxygen derivatives, providing deeper insights into flavin-catalyzed dioxygenation reactions, and is expected to inspire experimental design based on unconventional anionic peroxides.

Amplified and distinctive genotoxicity of titanium dioxide nanoparticles in transformed yeast reporters with human cytochrome P450 (CYP) genes

Titanium dioxide nanoparticles (nTiO2) have been considered a possible carcinogen to humans, but most existing studies have overlooked the role of human enzymes in assessing the genotoxicity of nTiO2. Here, a toxicogenomics-based in vitro genotoxicity assay using a GFP-fused yeast reporter library was employed to elucidate the genotoxic potential and mechanisms of nTiO2. Moreover, two new GFP-fused yeast reporter libraries containing either human CYP1A1 or CYP1A2 genes were constructed by transformation to investigate the potential modulation of nTiO2 genotoxicity in the presence of human CYP enzymes. This study found a lack of appreciable nTiO2 genotoxicity as indicated by the yeast reporter library in the absence of CYP expression but a significantly elevated indication of genotoxicity in either CYP1A1- or CYP1A2-expressing yeast. The intracellular reactive oxygen species (ROS) measurement indicated significantly higher ROS in yeast expressing either enzyme. The detected mitochondrial DNA damage suggested mitochondria as one of the target sites for oxidative damage by nTiO2 in the presence of either one of the CYP enzymes. The results thus indicated that the genotoxicity of nTiO2 was enhanced by human CYP1A1 or CYP1A2 enzyme and was associated with elevated oxidative stress, which suggested that the similar mechanisms could occur in human cells.

Molecular Basis for Chemoselectivity Control in Oxidations of Internal Aryl-Alkenes Catalyzed by Laboratory Evolved P450s

P450 enzymes naturally perform selective hydroxylations and epoxidations of unfunctionalized hydrocarbon substrates, among other reactions. The adaptation of P450 enzymes to a particular oxidative reaction involving alkenes is of great interest for the design of new synthetically useful biocatalysts. However, the mechanism that these enzymes utilize to precisely modulate the chemoselectivity and distinguishing between competing alkene double bond epoxidations and allylic C-H hydroxylations is sometimes not clear, which hampers the rational design of specific biocatalysts. In a previous work, a P450 from Labrenzia aggregata (P450LA1) was engineered in the laboratory using directed evolution to catalyze the direct oxidation of trans-β-methylstyrene to phenylacetone. The final variant, KS, was able to overcome the intrinsic preference for alkene epoxidation to directly generate a ketone product via the formation of a highly reactive carbocation intermediate. Here, additional library screening along this evolutionary lineage permitted to serendipitously detect a mutation that overcomes epoxidation and carbonyl formation by exhibiting a large selectivity of 94 % towards allylic C-H hydroxylation. A multiscalar computational methodology was applied to reveal the molecular basis towards this hydroxylation preference. Enzyme modelling suggests that introduction of a bulky substitution dramatically changes the accessible conformations of the substrate in the active site, thus modifying the enzymatic selectivity towards terminal hydroxylation and avoiding the competing epoxidation pathway, which is sterically hindered.

Enhanced Reactivities of Iron(IV)-Oxo Porphyrin Species in Oxidation Reactions Promoted by Intramolecular Hydrogen-Bonding

High-valent iron-oxo species are one of the common intermediates in both biological and biomimetic catalytic oxidation reactions. Recently, hydrogen-bonding (H-bonding) has been proved to be critical in determining the selectivity and reactivity. However, few examples have been established for mechanistic insights into the H-bonding effect. Moreover, intramolecular H-bonding effect on both C-H activation and oxygen atom transfer (OAT) reactions in synthetic porphyrin model system has not been investigated yet. In this study, a series of heme-containing iron(IV)-oxo porphyrin species with or without intramolecular H-bonding are synthesized and characterized. Kinetic studies revealed that intramolecular H-bonding can significantly enhance the reactivity of iron(IV)-oxo species in OAT, C-H activation, and electron-transfer reactions. This unprecedented unified H-bonding effect is elucidated by theoretical calculations, which showed that intramolecular H-bonding interactions lower the energy of the anti-bonding orbital of iron(IV)-oxo porphyrin species, resulting in the enhanced reactivities in oxidation reactions irrespective of the reaction type. To the best of the knowledge, this is the first extensive investigation on the intramolecular H-bonding effect in heme system. The results show that H-bonding interactions have a unified effect with iron(IV)-oxo porphyrin species in all three investigated reactions.

Substrate Conformational Switch Enables the Stereoselective Dimerization in P450 NascB: Insights from Molecular Dynamics Simulations and Quantum Mechanical/Molecular Mechanical Calculations

P450 NascB catalyzes the coupling of cyclo-(l-tryptophan-l-proline) (1) to generate (-)-naseseazine C (2) through intramolecular C-N bond formation and intermolecular C-C coupling. A thorough understanding of its catalytic mechanism is crucial for the engineering or design of P450-catalyzed C-N dimerization reactions. By employing MD simulations, QM/MM calculations, and enhanced sampling, we assessed various mechanisms from recent works. Our study demonstrates that the most favorable pathway entails the transfer of a hydrogen atom from N7-H to Cpd I. Subsequently, there is a conformational change in the substrate radical, shifting it from the Re-face to the Si-face of N7 in Substrate 1. The Si-face conformation of Substrate 1 is stabilized by the protein environment and the π-π stacking interaction between the indole ring and heme porphyrin. The subsequent intermolecular C3-C6' bond formation between Substrate 1 radical and Substrate 2 occurs via a radical attack mechanism. The conformational switch of the Substrate 1 radical not only lowers the barrier of the intermolecular C3-C6' bond formation but also yields the correct stereoselectivity observed in experiments. In addition, we evaluated the reactivity of the ferric-superoxide species, showing it is not reactive enough to initiate the hydrogen atom abstraction from the indole NH group of the substrate. Our simulation provides a comprehensive mechanistic insight into how the P450 enzyme precisely controls both the intramolecular C-N cyclization and intermolecular C-C coupling. The current findings align with the available experimental data, emphasizing the pivotal role of substrate dynamics in governing P450 catalysis.

Enhanced Reactivity through Equatorial Sulfur Coordination in Nonheme Iron(IV)-Oxo Complexes: Insights from Experiment and Theory

Sulfur ligation in metalloenzymes often gives the active site unique properties, whether it is the axial cysteinate ligand in the cytochrome P450s or the equatorial sulfur/thiol ligation in nonheme iron enzymes. To understand sulfur ligation to iron complexes and how it affects the structural, spectroscopic, and intrinsic properties of the active species and the catalysis of substrates, we pursued a systematic study and compared sulfur with amine-ligated iron(IV)-oxo complexes. We synthesized and characterized a biomimetic N4S-ligated iron(IV)-oxo complex and compared the obtained results with an analogous N5-ligated iron(IV)-oxo complex. Our work shows that the amine for sulfur replacement in the equatorial ligand framework leads to a rate enhancement for oxygen atom and hydrogen atom transfer reactions. Moreover, the sulfur-ligated iron(IV)-oxo complex reacts through a different reaction mechanism as compared to the N5-ligated iron(IV)-oxo complex, where the former reacts through hydride transfer with the latter reacting via radical pathways. We show that the reactivity differences are caused by a dramatic change in redox potential between the two complexes. Our studies highlight the importance of implementing a sulfur ligand into the equatorial ligand framework of nonheme iron(IV)-oxo complexes and how it affects the physicochemical properties of the oxidant and its reactivity.

On the engineering of reductase-based-monooxygenase activity in CYP450 peroxygenases

Recent bioengineering of CYP450OleT shows that peroxide-based CYP450OleT can be converted to a reductase-based self-sufficient enzyme, which is capable of showing efficient hydroxylation and decarboxylation activity for a wide range of substrates. The so-generated enzyme creates several mechanistic puzzles: (A) as CYP450 peroxygenases lack the conventional acid-alcohol pair, what is the source of two protons that are required to create the ultimate oxidant Cpd I? (B) Why is it only CYP450OleT that shows the reductase-based activity but no other CYP members? The present study provides a mechanistic solution to these puzzles using comprehensive MD simulations and hybrid QM/MM calculations. We show that the fusion of the reductase domain to the heme-binding domain triggers significant conformational rearrangement, which is gated by the propionate side chain, which constitutes a new water aqueduct via the carboxylate end of the substrate that ultimately participates in Cpd I formation. Importantly, such well-synchronized choreographies are controlled by remotely located Tyr359, which senses the fusion of reductase and communicates to the heme domain via non-covalent interactions. These findings provide crucial insights and a broader perspective which enables us to make a verifiable prediction: thus, the catalytic activity is not only limited to the first or second catalytic shell of an enzyme. Furthermore, it is predicted that reinstatement of tyrosine at a similar position in other members of CYP450 peroxygenases can convert these enzymes to reductase-based monooxygenases.

A frontier-orbital view of the initial steps of lytic polysaccharide monooxygenase reactions

Lytic polysaccharide monooxygenases (LPMOs) are copper enzymes that oxidatively cleave the strong C-H bonds in recalcitrant polysaccharide substrates, thereby playing a crucial role in biomass degradation. Recently, LPMOs have also been shown to be important for several pathogens. It is well established that the Cu(II) resting state of LPMOs is inactive, and the electronic structure of the active site needs to be altered to transform the enzyme into an active form. Whether this transformation occurs due to substrate binding or due to a unique priming reduction has remained speculative. Starting from four different crystal structures of the LPMO LsAA9A with well-defined oxidation states, we use a frontier molecular orbital approach to elucidate the initial steps of the LPMO reaction. We give an explanation for the requirement of the unique priming reduction and analyse electronic structure changes upon substrate binding. We further investigate how the presence of the substrate could facilitate an electron transfer from the copper active site to an H2O2 co-substrate. Our findings could help to control experimental LPMO reactions.

Unusual catalytic strategy by non-heme Fe(ii)/2-oxoglutarate-dependent aspartyl hydroxylase AspH

Biocatalytic C-H oxidation reactions are of important synthetic utility, provide a sustainable route for selective synthesis of important organic molecules, and are an integral part of fundamental cell processes. The multidomain non-heme Fe(ii)/2-oxoglutarate (2OG) dependent oxygenase AspH catalyzes stereoselective (3R)-hydroxylation of aspartyl- and asparaginyl-residues. Unusually, compared to other 2OG hydroxylases, crystallography has shown that AspH lacks the carboxylate residue of the characteristic two-His-one-Asp/Glu Fe-binding triad. Instead, AspH has a water molecule that coordinates Fe(ii) in the coordination position usually occupied by the Asp/Glu carboxylate. Molecular dynamics (MD) and quantum mechanics/molecular mechanics (QM/MM) studies reveal that the iron coordinating water is stabilized by hydrogen bonding with a second coordination sphere (SCS) carboxylate residue Asp721, an arrangement that helps maintain the six coordinated Fe(ii) distorted octahedral coordination geometry and enable catalysis. AspH catalysis follows a dioxygen activation-hydrogen atom transfer (HAT)-rebound hydroxylation mechanism, unusually exhibiting higher activation energy for rebound hydroxylation than for HAT, indicating that the rebound step may be rate-limiting. The HAT step, along with substrate positioning modulated by the non-covalent interactions with SCS residues (Arg688, Arg686, Lys666, Asp721, and Gln664), are essential in determining stereoselectivity, which likely proceeds with retention of configuration. The tetratricopeptide repeat (TPR) domain of AspH influences substrate binding and manifests dynamic motions during catalysis, an observation of interest with respect to other 2OG oxygenases with TPR domains. The results provide unique insights into how non-heme Fe(ii) oxygenases can effectively catalyze stereoselective hydroxylation using only two enzyme-derived Fe-ligating residues, potentially guiding enzyme engineering for stereoselective biocatalysis, thus advancing the development of non-heme Fe(ii) based biomimetic C-H oxidation catalysts, and supporting the proposal that the 2OG oxygenase superfamily may be larger than once perceived.

Computational Characterization of the Reactivity of Compound I in Unspecific Peroxygenases

Unspecific peroxygenases (UPOs) are emerging as promising biocatalysts for selective oxyfunctionalization of unactivated C-H bonds. However, their potential in large-scale synthesis is currently constrained by suboptimal chemical selectivity. Improving the selectivity of UPOs requires a deep understanding of the molecular basis of their catalysis. Recent molecular simulations have sought to unravel UPO's selectivity and inform their design principles. However, most of these studies focused on substrate-binding poses. Few researchers have investigated how the reactivity of CpdI, the principal oxidizing intermediate in the catalytic cycle, influences selectivity in a realistic protein environment. Moreover, the influence of protein electrostatics on the reaction kinetics of CpdI has also been largely overlooked. To bridge this gap, we used multiscale simulations to interpret the regio- and enantioselective hydroxylation of the n-heptane substrate catalyzed by Agrocybe aegerita UPO (AaeUPO). We comprehensively characterized the energetics and kinetics of the hydrogen atom-transfer (HAT) step, initiated by CpdI, and the subsequent oxygen rebound step forming the product. Notably, our approach involved both free energy and potential energy evaluations in a quantum mechanics/molecular mechanics (QM/MM) setting, mitigating the dependence of results on the choice of initial conditions. These calculations illuminate the thermodynamics and kinetics of the HAT and oxygen rebound steps. Our findings highlight that both the conformational selection and the distinct chemical reactivity of different substrate hydrogen atoms together dictate the regio- and enantio-selectivity. Building on our previous study of CpdI's formation in AaeUPO, our results indicate that the HAT step is the rate-limiting step in the overall catalytic cycle. The subsequent oxygen rebound step is swift and retains the selectivity determined by the HAT step. We also pinpointed several polar and charged amino acid residues whose electrostatic potentials considerably influence the reaction barrier of the HAT step. Notably, the Glu196 residue is pivotal for both the CpdI's formation and participation in the HAT step. Our research offers in-depth insights into the catalytic cycle of AaeUPO, which will be instrumental in the rational design of UPOs with enhanced properties.

Structural Characterization and Molecular Dynamics Study of the REPI Fusion Protein from Papaver somniferum L

REPI is a pivotal point enzyme in plant benzylisoquinoline alkaloid metabolism as it promotes the evolution of the biosynthetic branch of morphinan alkaloids. Experimental studies of its activity led to the identification of two modules (DRS and DRR) that catalyze two sequential steps of the epimerization of (S)- to (R)-reticuline. Recently, special attention has been paid to its genetic characterization and evolutionary history, but no structural analyses of the REPI protein have been conducted to date. We present here a computational structural characterization of REPI with heme and NADP cofactors in the apo state and in three complexes with substrate (S)-reticuline in DRS and intermediate 1,2-dehydroreticuline in DRS and in DRR. Since no experimental structure exists for REPI, we used its AlphaFold model as a scaffold to build up these four systems, which were submitted to all-atom molecular dynamics (MD) simulations. A comparison of MD results for the four systems revealed key dynamic changes associated with cofactor and ligand binding and provided a dynamic picture of the evolution of their structures and interactions. We also explored the possible dynamic occurrence of tunnels and electrostatic highways potentially involved in alternative mechanisms for channeling the intermediate from DRS to DRR.

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.

Understanding the Multifaceted Mechanism of Compound I Formation in Unspecific Peroxygenases through Multiscale Simulations

Unspecific peroxygenases (UPOs) can selectively oxyfunctionalize unactivated hydrocarbons by using peroxides under mild conditions. They circumvent the oxygen dilemma faced by cytochrome P450s and exhibit greater stability than the latter. As such, they hold great potential for industrial applications. A thorough understanding of their catalysis is needed to improve their catalytic performance. However, it remains elusive how UPOs effectively convert peroxide to Compound I (CpdI), the principal oxidizing intermediate in the catalytic cycle. Previous computational studies of this process primarily focused on heme peroxidases and P450s, which have significant differences in the active site from UPOs. Additionally, the roles of peroxide unbinding in the kinetics of CpdI formation, which is essential for interpreting existing experiments, have been understudied. Moreover, there has been a lack of free energy characterizations with explicit sampling of protein and hydration dynamics, which is critical for understanding the thermodynamics of the proton transport (PT) events involved in CpdI formation. To bridge these gaps, we employed multiscale simulations to comprehensively characterize the CpdI formation in wild-type UPO from Agrocybe aegerita (AaeUPO). Extensive free energy and potential energy calculations were performed in a quantum mechanics/molecular mechanics setting. Our results indicate that substrate-binding dehydrates the active site, impeding the PT from H2O2 to a nearby catalytic base (Glu196). Furthermore, the PT is coupled with considerable hydrogen bond network rearrangements near the active site, facilitating subsequent O-O bond cleavage. Finally, large unbinding free energy barriers kinetically stabilize H2O2 at the active site. These findings reveal a delicate balance among PT, hydration dynamics, hydrogen bond rearrangement, and cosubstrate unbinding, which collectively enable efficient CpdI formation. Our simulation results are consistent with kinetic measurements and offer new insights into the CpdI formation mechanism at atomic-level details, which can potentially aid the design of next-generation biocatalysts for sustainable chemical transformations of feedstocks.

Cytochromes P450 Associated with the Biosyntheses of Ribosomally Synthesized and Post-translationally Modified Peptides

Ribosomally synthesized and post-translationally modified peptides (RiPPs) are a class of exponentially increased natural products with characteristic chemical structures, topologies, and biosynthetic mechanisms as well as exceptional bioactivities including antibacteria, antitumors, and antiviruses. The biosynthesis of RiPP proceeds via a ribosomally assembled precursor peptide that undergoes varied post-translational modifications to generate a mature peptide. Cytochrome P450 (CYP or P450) monooxygenases are a superfamily of heme-containing enzymes that span a wide range of secondary metabolite biosynthetic pathways due to their broad substrate scopes and excellent catalytic versatility. In contrast to the enormous quantities of RiPPs and P450s, the P450 associated RiPP biosynthesis is comparatively limited, with most of their functions and timings remaining mysterious. Herein, this Review aims to provide an overview on the striking roles of P450s in RiPP biosyntheses uncovered to date and to illustrate their remarkable functions, mechanisms, as well as remaining challenges. This will shed light on novel P450 discovery and characterizations in RiPP biosyntheses.

Insights into Cytochrome P450 Enzyme Catalyzed Defluorination of Aromatic Fluorides

Density functional calculations establish a novel mechanism of aromatic defluorination by P450 Compound I. This is achieved via either an initial epoxide intermediate or through a 1,2-fluorine shift in an electrophilic intermediate, which highlights that the P450s can defluorinate fluoroarenes. However, in the absence of a proton donor a strong Fe-F bond can be obtained as shown from the calculations.

Melatonin Activation by Human Cytochrome P450 Enzymes: A Comparison between Different Isozymes

Cytochrome P450 enzymes in the human body play a pivotal role in both the biosynthesis and the degradation of the hormone melatonin. Melatonin plays a key role in circadian rhythms in the body, but its concentration is also linked to mood fluctuations as well as emotional well-being. In the present study, we present a computational analysis of the binding and activation of melatonin by various P450 isozymes that are known to yield different products and product distributions. In particular, the P450 isozymes 1A1, 1A2, and 1B1 generally react with melatonin to provide dominant aromatic hydroxylation at the C6-position, whereas the P450 2C19 isozyme mostly provides O-demethylation products. To gain insight into the origin of these product distributions of the P450 isozymes, we performed a comprehensive computational study of P450 2C19 isozymes and compared our work with previous studies on alternative isozymes. The work covers molecular mechanics, molecular dynamics and quantum mechanics approaches. Our work highlights major differences in the size and shape of the substrate binding pocket amongst the different P450 isozymes. Consequently, substrate binding and positioning in the active site varies substantially within the P450 isozymes. Thus, in P450 2C19, the substrate is oriented with its methoxy group pointing towards the heme, and therefore reacts favorably through hydrogen atom abstraction, leading to the production of O-demethylation products. On the other hand, the substrate-binding pockets in P450 1A1, 1A2, and 1B1 are tighter, direct the methoxy group away from the heme, and consequently activate an alternative site and lead to aromatic hydroxylation instead.

Fluorescence-Based High-Throughput Assays for Investigating Cytochrome P450 Enzyme-Mediated Drug-Drug Interactions

The cytochrome P450 enzymes (CYPs), a group of heme-containing enzymes, catalyze oxidative metabolism of a wide range of drugs and xenobiotics, as well as different endogenous molecules. Strong inhibition of human CYPs is the most common cause of clinically associated pharmacokinetic drug-drug/herb-drug interactions (DDIs/HDIs), which may result in serious adverse drug reactions, even toxicity. Accurate and rapid assessing of the inhibition potentials on CYP activities for therapeutic agents is crucial for the prediction of clinically relevant DDIs/HDIs. Over the past few decades, significant efforts have been invested into developing optical substrates for the human CYPs, generating a variety of powerful tools for high-throughput assays to detect CYP activities in biologic specimens and for screening of CYP inhibitors. This minireview focuses on recent advances in optical substrates developments for human CYPs, as well as their applications in screening CYP inhibitors and DDIs/HDIs studies. The examples for rational design and optimization of highly specific optical substrates for the target CYP enzyme, as well as applications in investigating CYP-mediated DDIs, are illustrated. Finally, the challenges and future perspectives in this field are proposed. Collectively, this review summarizes the reported optical-based biochemical assays for highly efficient CYP activities detection, which strongly facilitated the discovery of CYP inhibitors and the investigations on CYP-mediated DDIs. SIGNIFICANCE STATEMENT: Optical substrates for cytochrome P450 enzymes (CYPs) have emerged as powerful tools for the construction of high-throughput assays for screening of CYP inhibitors. This mini-review covers the advances and challenges in the development of highly specific optical substrates for sensing human CYP isoenzymes, as well as their applications in constructing fluorescence-based high-throughput assays for investigating CYP-mediated drug-drug interactions.

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.

Rieske Oxygenases and Other Ferredoxin-Dependent Enzymes: Electron Transfer Principles and Catalytic Capabilities

Enzymes that depend on sophisticated electron transfer via ferredoxins (Fds) exhibit outstanding catalytic capabilities, but despite decades of research, many of them are still not well understood or exploited for synthetic applications. This review aims to provide a general overview of the most important Fd-dependent enzymes and the electron transfer processes involved. While several examples are discussed, we focus in particular on the family of Rieske non-heme iron-dependent oxygenases (ROs). In addition to illustrating their electron transfer principles and catalytic potential, the current state of knowledge on structure-function relationships and the mode of interaction between the redox partner proteins is reviewed. Moreover, we highlight several key catalyzed transformations, but also take a deeper dive into their engineerability for biocatalytic applications. The overall findings from these case studies highlight the catalytic capabilities of these biocatalysts and could stimulate future interest in developing additional Fd-dependent enzyme classes for synthetic applications.

Nanodisc-embedded cytochrome P450 P3A4 binds diverse ligands by distributing conformational dynamics to its flexible elements

Cytochrome P450 3A4 (CYP3A4) metabolizes a wide range of drugs and toxins. Interactions of CYP3A4 with ligands are difficult to predict due to promiscuity and conformational flexibility. To better understand CYP3A4 conformational responses to ligands we use hydrogen deuterium exchange mass spectrometry (HDX-MS) to investigate the effect of ligands on nanodisc-embedded CYP3A4. For a subset of CYP3A4-ligand complexes, differences in the low-frequency modes derived by principal component analyses of molecular dynamics trajectories mirrored the HDX-MS results. The effects of ligands are distributed to flexible elements of CYP3A4 between stretches of secondary structure. The largest effects occur in the F- and G-helices, where most ligands increase the flexibility of the F-helix and connecting loops and decrease the flexibility of the C-term of the G-helix. Most ligands affect the E-F-G, CD and HI regions of the protein. Ligand-dependent differences are observed in the A"-A' loop, BC region, E-helix, K-β1 region, proximal loop, and C-term loop. Correlated HDX responses were observed in the CD region and the C-term of the G-helix that were most pronounced for Type II ligands. Collectively, the HDX and molecular dynamics results suggest that CYP3A4 accommodates diverse binding partners by propagating local backbone fluctuations from the binding site onto the flexible regions of the enzyme via long-range interactions that are differentially modulated by ligands. In contrast to the paradigm wherein ligands decrease protein dynamics at their binding site, a wide range of ligands modestly increase CYP3A4 dynamics throughout the protein including effects remote from the active site.

Bioinspired Active Site with a Coordination-Adaptive Organosulfonate Ligand for Catalytic Water Oxidation at Neutral pH

Many enzymes use adaptive frameworks to preorganize substrates, accommodate various structural and electronic demands of intermediates, and accelerate related catalysis. Inspired by biological systems, a Ru-based molecular water oxidation catalyst containing a configurationally labile ligand [2,2':6',2″-terpyridine]-6,6″-disulfonate was designed to mimic enzymatic framework, in which the sulfonate coordination is highly flexible and functions as both an electron donor to stabilize high-valent Ru and a proton acceptor to accelerate water dissociation, thus boosting the catalytic water oxidation performance thermodynamically and kinetically. The combination of single-crystal X-ray analysis, various temperature NMR, electrochemical techniques, and DFT calculations was utilized to investigate the fundamental role of the self-adaptive ligand, demonstrating that the on-demand configurational changes give rise to fast catalytic kinetics with a turnover frequency (TOF) over 2000 s^-1, which is compared to oxygen-evolving complex (OEC) in natural photosynthesis.

Conversion and synthesis of chemicals catalyzed by fungal cytochrome P450 monooxygenases: A review

Cytochrome P450s (also called CYPs or P450s) are a superfamily of heme-containing monooxygenases. They are distributed in all biological kingdoms. Most fungi have at least two P450-encoding genes, CYP51 and CYP61, which are housekeeping genes that play important roles in the synthesis of sterols. However, the kingdom fungi is an interesting source of numerous P450s. Here, we review reports on fungal P450s and their applications in the bioconversion and biosynthesis of chemicals. We highlight their history, availability, and versatility. We describe their involvement in hydroxylation, dealkylation, oxygenation, C═C epoxidation, C-C cleavage, C-C ring formation and expansion, C-C ring contraction, and uncommon reactions in bioconversion and/or biosynthesis pathways. The ability of P450s to catalyze these reactions makes them promising enzymes for many applications. Thus, we also discuss future prospects in this field. We hope that this review will stimulate further study and exploitation of fungal P450s for specific reactions and applications.

β-Cryptoxanthin Production in Escherichia coli by Optimization of the Cytochrome P450 CYP97H1 Activity

Cytochromes P450, forming a superfamily of monooxygenases containing heme as a cofactor, show great versatility in substrate specificity. Metabolic engineering can take advantage of this feature to unlock novel metabolic pathways. However, the cytochromes P450 often show difficulty being expressed in a heterologous chassis. As a case study in the prokaryotic host Escherichia coli, the heterologous synthesis of β-cryptoxanthin was addressed. This carotenoid intermediate is difficult to produce, as its synthesis requires a monoterminal hydroxylation of β-carotene whereas most of the classic carotene hydroxylases are dihydroxylases. This study was focused on the optimization of the in vivo activity of CYP97H1, an original P450 β-carotene monohydroxylase. Engineering the N-terminal part of CYP97H1, identifying the matching redox partners, defining the optimal cellular background and adjusting the culture and induction conditions improved the production by 400 times compared to that of the initial strain, representing 2.7 mg/L β-cryptoxanthin and 20% of the total carotenoids produced.

Melatonin Activation by Cytochrome P450 Isozymes: How Does CYP1A2 Compare to CYP1A1?

Cytochrome P450 enzymes are versatile enzymes found in most biosystems that catalyze mono-oxygenation reactions as a means of biosynthesis and biodegradation steps. In the liver, they metabolize xenobiotics, but there are a range of isozymes with differences in three-dimensional structure and protein chain. Consequently, the various P450 isozymes react with substrates differently and give varying product distributions. To understand how melatonin is activated by the P450s in the liver, we did a thorough molecular dynamics and quantum mechanics study on cytochrome P450 1A2 activation of melatonin forming 6-hydroxymelatonin and N-acetylserotonin products through aromatic hydroxylation and O-demethylation pathways, respectively. We started from crystal structure coordinates and docked substrate into the model, and obtained ten strong binding conformations with the substrate in the active site. Subsequently, for each of the ten substrate orientations, long (up to 1 μs) molecular dynamics simulations were run. We then analyzed the orientations of the substrate with respect to the heme for all snapshots. Interestingly, the shortest distance does not correspond to the group that is expected to be activated. However, the substrate positioning gives insight into the protein residues it interacts with. Thereafter, quantum chemical cluster models were created and the substrate hydroxylation pathways calculated with density functional theory. These relative barrier heights confirm the experimental product distributions and highlight why certain products are obtained. We make a detailed comparison with previous results on CYP1A1 and identify their reactivity differences with melatonin.

Human Cytochrome P450 1, 2, 3 Families as Pharmacogenes with Emphases on Their Antimalarial and Antituberculosis Drugs and Prevalent African Alleles

Precision medicine gives individuals tailored medical treatment, with the genotype determining the therapeutic strategy, the appropriate dosage, and the likelihood of benefit or toxicity. Cytochrome P450 (CYP) enzyme families 1, 2, and 3 play a pivotal role in eliminating most drugs. Factors that affect CYP function and expression have a major impact on treatment outcomes. Therefore, polymorphisms of these enzymes result in alleles with diverse enzymatic activity and drug metabolism phenotypes. Africa has the highest CYP genetic diversity and also the highest burden of malaria and tuberculosis, and this review presents current general information on CYP enzymes together with variation data concerning antimalarial and antituberculosis drugs, while focusing on the first three CYP families. Afrocentric alleles such as CYP2A6*17, CYP2A6*23, CYP2A6*25, CYP2A6*28, CYP2B6*6, CYP2B6*18, CYP2C8*2, CYP2C9*5, CYP2C9*8, CYP2C9*9, CYP2C19*9, CYP2C19*13, CYP2C19*15, CYP2D6*2, CYP2D6*17, CYP2D6*29, and CYP3A4*15 are implicated in diverse metabolic phenotypes of different antimalarials such as artesunate, mefloquine, quinine, primaquine, and chloroquine. Moreover, CYP3A4, CYP1A1, CYP2C8, CYP2C18, CYP2C19, CYP2J2, and CYP1B1 are implicated in the metabolism of some second-line antituberculosis drugs such as bedaquiline and linezolid. Drug-drug interactions, induction/inhibition, and enzyme polymorphisms that influence the metabolism of antituberculosis, antimalarial, and other drugs, are explored. Moreover, a mapping of Afrocentric missense mutations to CYP structures and a documentation of their known effects provided structural insights, as understanding the mechanism of action of these enzymes and how the different alleles influence enzyme function is invaluable to the advancement of precision medicine.

Efficient Production of Epoxy-Norbornane from Norbornene by an Engineered P450 Peroxygenase

Epoxy-norbornane (EPO-NBE) is a crucial building block for the synthesis of various biologically active heterocyclic systems. To develop an efficient protocol for producing EPO-NBE using norbornene (NBE) as a substrate, cytochrome P450 enzyme from Pseudomonas putida (CYP238A1) was examined and its crystal structure (PDB code: 7X53) was resolved. Molecular mechanism analysis showed a high energy barrier related to iron-alkoxy radical complex formation. Therefore, a protein engineering strategy was developed and an optimal CYP238A1^NPV variant containing a local hydrophobic "fence" at the active site was obtained, which increased the H₂ O₂ -dependent epoxidation activity by 7.5-fold compared with that of CYP238A1^WT . Among the "fence", Glu255 participates in an efficient proton transfer system. Whole-cell transformation using CYP238A1^NPV achieved an EPO-NBE yield of 77.6 g ⋅ L^-1 in a 30-L reactor with 66.3 % conversion. These results demonstrate the potential of this system for industrial production of EPO-NBE and provides a new biocatalytic platform for epoxidation chemistry.

Biotransformation of Bisphenol by Human Cytochrome P450 2C9 Enzymes: A Density Functional Theory Study

Bisphenol A (BPA, 2,2-bis-(4-hydroxyphenyl)propane) is used as a precursor in the synthesis of polycarbonate and epoxy plastics; however, its availability in the environment is causing toxicity as an endocrine-disrupting chemical. Metabolism of BPA and their analogues (substitutes) is generally performed by liver cytochrome P450 enzymes and often leads to a mixture of products, and some of those are toxic. To understand the product distributions of P450 activation of BPA, we have performed a computational study into the mechanisms and reactivities using large model structures of a human P450 isozyme (P450 2C9) with BPA bound. Density functional theory (DFT) calculations on mechanisms of BPA activation by a P450 compound I model were investigated, leading to a number of possible products. The substrate-binding pocket is tight, and as a consequence, aliphatic hydroxylation is not feasible as the methyl substituents of BPA cannot reach compound I well due to constraints of the substrate-binding pocket. Instead, we find low-energy pathways that are initiated with phenol hydrogen atom abstraction followed by OH rebound to the phenolic ortho- or para-position. The barriers of para-rebound are well lower in energy than those for ortho-rebound, and consequently, our P450 2C9 model predicts dominant hydroxycumyl alcohol products. The reactions proceed through two-state reactivity on competing doublet and quartet spin state surfaces. The calculations show fast and efficient substrate activation on a doublet spin state surface with a rate-determining electrophilic addition step, while the quartet spin state surface has multiple high-energy barriers that can also lead to various side products including C⁴-aromatic hydroxylation. This work shows that product formation is more feasible on the low spin state, while the physicochemical properties of the substrate govern barrier heights of the rate-determining step of the reaction. Finally, the importance of the second-coordination sphere is highlighted that determines the product distributions and guides the bifurcation pathways.

The Performance of Different Water Models on the Structure and Function of Cytochrome P450 Enzymes

Modeling approaches and modern simulations to investigate the biomolecular structure and function rely on various methods. Since water molecules play a crucial role in all sorts of chemistry, the accurate modeling of water molecules is vital for such simulations. In cytochrome P450 (CYP450), in particular, water molecules play a key role in forming active oxidant that ultimately performs oxidation and metabolism. In the present study, we have highlighted the behavior of the three most widely used water models─TIP3P, SPC/E, and OPC─for three different CYP450 enzymes─CYP450_BM3, CYP450_OleT, and CYP450_BSβ─during MD simulations and QM/MM calculations. We studied the various properties, such as RMSD, RMSF, H-bond, water occupancy, and hydrogen atom transfer (HAT), using QM/MM calculations and compared them for all three water models. Our study shows that the stabilities of the enzyme complexes are well maintained in all three water models. However, the OPC water model performs well for the polar active sites, that is, in CYP450_OleT and CYP450_BSβ, while the TIP3P water model is superior for the hydrophobic site, such as CYP450_BM3.

Local Electric Fields: From Enzyme Catalysis to Synthetic Catalyst Design

This Mini-Review Article outlines recent advances in the study of local electric field (LEF) governed enzyme catalysis and the application of the LEF principle in synthetic catalyst design. We start by discussing the electrostatics principles that drive enzyme catalysis, and its experimental verifications through vibrational Stark spectroscopy. Subsequently, we describe aspects of LEFs other than catalysis, i.e., induction of mechanistic crossovers, among others. Here, we focus on the early work done using computational tools, along with some recent contributions. Following an in-depth discussion of the role of LEFs in enzyme catalysis, we then highlight some recent works on designed local electric fields (D-LEF) and their applications in organic synthesis. Subsequently, we turn to D-LEFs in synthetic enzymes and supramolecular systems (cf. the work by the Head-Gordon group). We end by discussing some of the software packages that have been developed to analyze local electric fields computationally. Overall, the present Mini-Review Article paints an insightful picture of the current state of the art using LEF in enzyme catalysis and its application for further bioengineering and synthetic organic frameworks in a broad perspective.

Computations reveal a crucial role of an aromatic dyad in the catalytic function of plant cytochrome P450 mint superfamily

Enzymes are highly specific for their native functions, however with advances in bioengineering tools such as directed evolution, several enzymes are being repurposed for the secondary function of contemporary significance(Khersonsky and Tawfik, 2010 [1]⁾. Due to the functional versatility, the Cytochrome P450 (CYP450) superfamily has become the ideal scaffold for such bioengineering. In the current study, using MD (molecular dynamics) simulations and hybrid QM/MM (Quantum mechanics/molecular mechanics) calculations, we have studied the mechanism of spontaneous emergence of a secondary function due to a single site mutation in two plant CYP450 enzymes from the mint family. The MD simulations of WT (wild type) CYP71D18 and CYP71D13 enzymes and their variants show a crucial gating mechanism by aromatic dyad formed by Phe121 and Phe363 which regulates the substrate recognition. The QM/MM calculations reveal that the hydroxylation reactions at C3 and C6 positions in WT CYP71D18 and CYP71D13 enzymes as well as their variants follow a hydrogen atom transfer (HAT) followed by a single electron transfer (SET) mechanism, which is different from the typical rebound mechanism shown by most of the CYP450 enzymes.

Mechanism of Melatonin Metabolism by CYP1A1: What Determines the Bifurcation Pathways of Hydroxylation versus Deformylation?

Melatonin, a widely applied cosmetic active ingredient, has a variety of uses as a skin protector through antioxidant and anti-inflammatory functions as well as giving the body UV-induced defenses and immune system support. In the body, melatonin is synthesized from a tryptophan amino acid in a cascade of reactions, but as melatonin is toxic at high concentrations, it is metabolized in the human skin by the cytochrome P450 enzymes. The P450s are diverse heme-based mono-oxygenases that catalyze oxygen atom-transfer processes that trigger metabolism and detoxification reactions in the body. In the catalytic cycle of the P450s, a short-lived high-valent iron(IV)-oxo heme cation radical is formed that has been proposed to be the active oxidant. How and why it activates melatonin in the human body and what the origin of the product distributions is, are unknown. This encouraged us to do a detailed computational study on a typical human P450 isozyme, namely CYP1A1. We initially did a series of molecular dynamics simulations with substrate docked into several orientations. These simulations reveal a number of stable substrate-bound positions in the active site, which may lead to differences in substrate activation channels. Using tunneling analysis on the full protein structures, we show that two of the four binding conformations lead to open substrate-binding pockets. As a result, in these open pockets, the substrate is not tightly bound and can escape back into the solution. In the closed conformations, in contrast, the substrate is mainly oriented with the methoxy group pointing toward the heme, although under a different angle. We then created large quantum cluster models of the enzyme and focused on the chemical reaction mechanisms for melatonin activation, leading to competitive O-demethylation and C⁶-aromatic hydroxylation pathways. The calculations show that active site positioning determines the product distributions, but the bond that is activated is not necessarily closest to the heme in the enzyme-substrate complex. As such, the docking and molecular dynamics positioning of the substrate versus oxidant can give misleading predictions on product distributions. In particular, in quantum mechanics cluster model I, we observe that through a tight hydrogen bonding network, a preferential 6-hydroxylation of melatonin is obtained. However, O-demethylation becomes possible in alternative substrate-binding orientations that have the C⁶-aromatic ring position shielded. Finally, we investigated enzymatic and non-enzymatic O-demethylation processes and show that the hydrogen bonding network in the substrate-binding pocket can assist and perform this step prior to product release from the enzyme.

Cytochrome P450 monooxygenase-mediated tailoring of triterpenoids and steroids in plants

The cytochrome P450 monooxygenase (CYP) superfamily comprises hemethiolate enzymes that perform remarkable regio- and stereospecific oxidative chemistry. As such, CYPs are key agents for the structural and functional tailoring of triterpenoids, one of the largest classes of plant natural products with widespread applications in pharmaceuticals, food, cosmetics, and agricultural industries. In this review, we provide a full overview of 149 functionally characterised CYPs involved in the biosynthesis of triterpenoids and steroids in primary as well as in specialised metabolism. We describe the phylogenetic distribution of triterpenoid- and steroid-modifying CYPs across the plant CYPome, present a structure-based summary of their reactions, and highlight recent examples of particular interest to the field. Our review therefore provides a comprehensive up-to-date picture of CYPs involved in the biosynthesis of triterpenoids and steroids in plants as a starting point for future research.

Bioimaging agents based on redox-active transition metal complexes

Detecting the fluctuation and distribution of various bioactive species in biological systems is of great importance in determining diseases at their early stages. Metal complex-based probes have attracted considerable attention in bioimaging applications owing to their unique advantages, such as high luminescence, good photostability, large Stokes shifts, low toxicity, and good biocompatibility. In this review, we summarized the development of redox-active transition metal complex-based probes in recent five years with the metal ions of iron, manganese, and copper, which play essential roles in life and can avoid the introduction of exogenous metals into biological systems. The designing principles that afford these complexes with optical or magnetic resonance (MR) imaging properties are elucidated. The applications of the complexes for bioimaging applications of different bioactive species are demonstrated. The current challenges and potential future directions of these probes for applications in biological systems are also discussed.

Oriented External Electric Fields Regurating the Reaction Mechanism of CH4 Oxidation Catalyzed by Fe(IV)-Oxo-Corrolazine: Insight from Density Functional Calculations

Methane is the simplest alkane and can be used as an alternative energy source for oil and coal, but the greenhouse effect caused by its leakage into the air is not negligible, and its conversion into liquid methanol not only facilitates transportation, but also contributes to carbon neutrality. In order to find an efficient method for converting methane to methanol, CH4 oxidation catalyzed by Fe(IV)-Oxo-corrolazine (Fe(IV)-Oxo-Cz) and its reaction mechanism regulation by oriented external electric fields (OEEFs) are systematically studied by density functional calculations. The calculations show that Fe(IV)-Oxo-Cz can abstract one H atom from CH4 to form the intermediate with OH group connecting on the corrolazine ring, with the energy barrier of 25.44 kcal mol-1. And then the product methanol is formed through the following rebound reaction. Moreover, the energy barrier can be reduced to 20.72 kcal mol-1 through a two-state reaction pathway. Furthermore, the effect of OEEFs on the reaction is investigated. We found that OEEFs can effectively regulate the reaction by adjusting the stability of the reactant and the transition state through the interaction of electric field-molecular dipole moment. When the electric field is negative, the energy barrier of the reaction decreases with the increase of electric intensity. Moreover, the OEEF aligned along the intrinsic Fe‒O reaction axis can effectively regulate the ability of forming the OH on the corrolazine ring by adjusting the charges of O and H atoms. When the electric field intensity is -0.010 a.u., the OH can be directly rebounded to the CH3· before it is connecting on the corrolazine ring, thus forming the product directly from the transition state without passing through the intermediate with only an energy barrier of 17.34 kcal mol-1, which greatly improves the selectivity of the reaction.

Single-Site Mutation Induces Water-Mediated Promiscuity in Lignin Breaking Cytochrome P450GcoA

Cytochrome P450_GcoA is an enzyme that catalyzes the guaiacol unit of lignin during the lignin breakdown via an aryl-O-demethylation reaction. This reaction is intriguing and is of commercial importance for its potential applications in the production of biofuel and plastic from biomass feedstock. Recently, the F169A mutation in P450_GcoA elicits a promiscuous activity for syringol while maintaining the native activity for guaiacol. Using comprehensive MD simulations and hybrid QM/MM calculations, we address, herein, the origin of promiscuity in P450_GcoA and its relevance to the specific activity toward lignin-derived substrates. Our study shows a crucial role of an aromatic dyad of F169 and F395 by regulating the water access to the catalytic center. The F169A mutation opens a water aqueduct and hence increases the native activity for G-lignin. We show that syringol binds very tightly to the WT enzyme, which blocks the conformational rearrangement needed for the second step of O-demethylation. The F169A creates an extra room favoring the conformational rearrangement in the 3-methoxycatechol (3MC) and second dose of the dioxygen insertion. Therefore, using MD simulations and complemented by thorough QM/MM calculations, our study shows how a single-site mutation rearchitects active site engineering for promiscuous syringol activity.

Combined MD and QM/MM Calculations Reveal Allostery-Driven Promiscuity in Dipeptide Epimerases of Enolase Family

The adaptability of the active site to amplify the secondary function is supposed to be the fundamental cause of the promiscuity and the evolution of new functions in enzymes. In most cases, mutations occur close to the active site and/or in the catalytic site to change the active site plasticity to accommodate the non-native substrate. In the present study, using MD simulations and hybrid QM/MM calculations, we have shown a way to enhance the promiscuity, i. e., the allostery-driven promiscuity. Using a case study of the AEE enzyme where the capping loop recognizes the substrate, herein, we show that a single site mutation (D321G) far from the capping loop can induce a large conformational change in the capping loop to recognize different substrates for different functions. The QM/MM calculations for the WT and mutated enzyme provide a first validation of the mechanism of 1,1-proton transfer and dehydration by the AEE enzyme. Since AEE epimerase possesses a highly conserved TIM-barrel fold, we believe that our study provides a crucial lead to understanding the mechanism of emergence of secondary function which can be useful to repurpose ancient enzymes for modern usage.

Heme compound II models in chemoselectivity and disproportionation reactions

Heme compound II models bearing electron-deficient and -rich porphyrins, [FeIV(O)(TPFPP)(Cl)]- (1a) and [FeIV(O)(TMP)(Cl)]- (2a), respectively, are synthesized, spectroscopically characterized, and investigated in chemoselectivity and disproportionation reactions using cyclohexene as a mechanistic probe. Interestingly, cyclohexene oxidation by 1a occurs at the allylic C-H bonds with a high kinetic isotope effect (KIE) of 41, yielding 2-cyclohexen-1-ol product; this chemoselectivity is the same as that of nonheme iron(iv)-oxo intermediates. In contrast, as observed in heme compound I models, 2a yields cyclohexene oxide product with a KIE of 1, demonstrating a preference for C[double bond, length as m-dash]C epoxidation. The latter result is interpreted as 2a disproportionating to form [FeIV(O)(TMP+˙)]+ (2b) and FeIII(OH)(TMP), and 2b becoming the active oxidant to conduct the cyclohexene epoxidation. In contrast to 2a, 1a does not disproportionate under the present reaction conditions. DFT calculations confirm that compound II models prefer C-H bond hydroxylation and that disproportionation of compound II models is controlled thermodynamically by the porphyrin ligands. Other aspects, such as acid and base effects on the disproportionation of compound II models, have been discussed as well.