QM/MM modeling of compound I active species in cytochrome P450, cytochrome C peroxidase, and ascorbate peroxidase

QM/MM Study of the Reaction Mechanism of L-Tyrosine Hydroxylation Catalyzed by the Enzyme CYP76AD1

We have studied the hydroxylation mechanism of l-Tyr by the heme-dependent enzyme CYP76AD1 from the sugar beet (Beta vulgaris). This enzyme has a promising biotechnological application in modified yeast strains to produce medicinal alkaloids, an alternative to the traditional opium poppy harvest. A generative machine learning software based on AlphaFold was used to build the structure of CYP76AD1 since there are no structural data for this specific enzyme. After model validation, l-Tyr was docked in the active site of CYP76AD1 to assemble the reactive complex, whose catalytic distances remained stable throughout the 100 ns of MD simulation. Subsequent QM/MM calculations elucidated that l-Tyr hydroxylation occurs in two steps: hydrogen abstraction from l-Tyr by CpdI, forming an l-Tyr radical, and subsequent radical rebound, corresponding to a rate-limiting step of 16.0 kcal·mol-1. Our calculations suggest that the hydrogen abstraction step should occur in the doublet state, while the radical rebound should happen in the quartet state. The clarification of the reaction mechanism of CYP76AD1 provides insights into the rational optimization of the biosynthesis of alkaloids to eliminate the use of opium poppy.

Correlating Structure with Spectroscopy in Ascorbate Peroxidase Compound II

Structural and spectroscopic investigations of compound II in ascorbate peroxidase (APX) have yielded conflicting conclusions regarding the protonation state of the crucial Fe(IV) intermediate. Neutron diffraction and crystallographic data support an iron(IV)-hydroxo formulation, whereas Mössbauer, X-ray absorption (XAS), and nuclear resonance vibrational spectroscopy (NRVS) studies appear consistent with an iron(IV)-oxo species. Here we examine APX with spectroscopy-oriented QM/MM calculations and extensive exploration of the conformational space for both possible formulations of compound II. We establish that irrespective of variations in the orientation of a vicinal arginine residue and potential reorganization of proximal water molecules and hydrogen bonding, the Fe-O distances for the oxo and hydroxo forms consistently fall within distinct, narrow, and nonoverlapping ranges. The accuracy of geometric parameters is validated by coupled-cluster calculations with the domain-based local pair natural orbital approach, DLPNO-CCSD(T). QM/MM calculations of spectroscopic properties are conducted for all structural variants, encompassing Mössbauer, optical, X-ray absorption, and X-ray emission spectroscopies and NRVS. All spectroscopic observations can be assigned uniquely to an Fe(IV)═O form. A terminal hydroxy group cannot be reconciled with the spectroscopic data. Under no conditions can the Fe(IV)═O distance be sufficiently elongated to approach the crystallographically reported Fe-O distance. The latter is consistent only with a hydroxo species, either Fe(IV) or Fe(III). Our findings strongly support the Fe(IV)═O formulation of APX-II and highlight unresolved discrepancies in the nature of samples used across different experimental studies.

QM/MM Modeling Aided Enzyme Engineering in Natural Products Biosynthesis

Natural products and their derivatives are widely used across various industries, particularly pharmaceuticals. Modern engineered biosynthesis provides an alternative way of producing and meeting the growing need for diverse natural products. Natural enzymes, on the other hand, often exhibit unsatisfactory catalytic characteristics and necessitate further enzyme engineering modifications. QM/MM, as a powerful and extensively used computational tool in the field of enzyme catalysis, has been increasingly applied in rational enzyme engineering over the past decade. In this review, we summarize recent advances in QM/MM computational investigation on enzyme catalysis and enzyme engineering for natural product biosynthesis. The challenges and perspectives for future QM/MM applications aided enzyme engineering in natural product biosynthesis will also be discussed.

Potential Role of Lysine Acetylation and Autophagy in Brown Film Formation and Postripening of Lentinula edodes Mycelium

Lentinula edodes is one of the most widely cultivated edible mushrooms in the world. When cultivated in sawdust, the surface mycelium of L. edodes needs a long postripening stage wherein it forms a brown film (BF) by secreting and accumulating pigments. BF formation is critical for the high quality and yield of fruiting bodies. Protein lysine acetylation (KAC) is an important post-translational modification that regulates growth and development. Previous studies have shown that deacetylase levels are significantly increased during BF formation in the postripening stage of L. edodes. The aim of this study was to assess the role of protein acetylation during BF formation. To this end, we compared the acetylome of L. edodes mycelia before and after BF formation using anti-acetyl antibody-based label-free quantitative proteomics. We identified 5,613 acetylation sites in 1,991 proteins, and quantitative information was available for 4,848 of these sites in 1,815 proteins. Comparative acetylome analysis showed that the modification of 699 sites increased and that of 562 sites decreased during BF formation. Bioinformatics analysis of the differentially acetylated proteins showed significant enrichment in the tricarboxylic acid (TCA) cycle and proteasome pathways. Furthermore, functional assays showed that BF formation is associated with significant changes in the activities of proteasome, citrate synthase, and isocitrate dehydrogenase. Consistent with this hypothesis, the lysine deacetylase inhibitor trichostatin (TSA) delayed autophagy and BF formation in L. edodes. Taken together, KAC and autophagy play important roles in the mycelial BF formation and postripening stage of L. edodes. IMPORTANCE Mycelial BF formation and postripening of L. edodes affects the quality and quantity of its edible fruiting bodies. In this study, we explored the role of protein KAC in this biological process, with the aim of optimizing the cultivation and yield of L. edodes.

Catalysis and Electron Transfer in De Novo Designed Metalloproteins

One of the hallmark advances in our understanding of metalloprotein function is showcased in our ability to design new, non-native, catalytically active protein scaffolds. This review highlights progress and milestone achievements in the field of de novo metalloprotein design focused on reports from the past decade with special emphasis on de novo designs couched within common subfields of bioinorganic study: heme binding proteins, monometal- and dimetal-containing catalytic sites, and metal-containing electron transfer sites. Within each subfield, we highlight several of what we have identified as significant and important contributions to either our understanding of that subfield or de novo metalloprotein design as a discipline. These reports are placed in context both historically and scientifically. General suggestions for future directions that we feel will be important to advance our understanding or accelerate discovery are discussed.

The redox potential of a heme cofactor in Nitrosomonas europaea cytochrome c peroxidase: a polarizable QM/MM study

Redox reactions are crucial to biological processes that protect organisms against oxidative stress. Metalloenzymes, such as peroxidases which reduce excess reactive oxygen species into water, play a key role in detoxification mechanisms. Here we present the results of a polarizable QM/MM study of the reduction potential of the electron transfer heme in the cytochrome c peroxidase of Nitrosomonas europaea. We have found that environment polarization does not substantially affect the computed value of the redox potential. Particular attention has been given to analyzing the role of electrostatic interactions within the protein environment and the solvent on tuning the redox potential of the heme co-factor. We have found that the electrostatic interactions predominantly explain the fluctuations of the vertical ionization/attachment energies of the heme for the sampled configurations, and that the long range electrostatic interactions (up to 40 Å) contribute substantially to the absolute values of the vertical energy gaps.

Product Distributions of Cytochrome P450 OleTJE with Phenyl-Substituted Fatty Acids: A Computational Study

There are two types of cytochrome P450 enzymes in nature, namely, the monooxygenases and the peroxygenases. Both enzyme classes participate in substrate biodegradation or biosynthesis reactions in nature, but the P450 monooxygenases use dioxygen, while the peroxygenases take H2O2 in their catalytic cycle instead. By contrast to the P450 monooxygenases, the P450 peroxygenases do not require an external redox partner to deliver electrons during the catalytic cycle, and also no external proton source is needed. Therefore, they are fully self-sufficient, which affords them opportunities in biotechnological applications. One specific P450 peroxygenase, namely, P450 OleTJE, reacts with long-chain linear fatty acids through oxidative decarboxylation to form hydrocarbons and, as such, has been implicated as a suitable source for the biosynthesis of biofuels. Unfortunately, the reactions were shown to produce a considerable amount of side products originating from Cα and Cβ hydroxylation and desaturation. These product distributions were found to be strongly dependent on whether the substrate had substituents on the Cα and/or Cβ atoms. To understand the bifurcation pathways of substrate activation by P450 OleTJE leading to decarboxylation, Cα hydroxylation, Cβ hydroxylation and Cα-Cβ desaturation, we performed a computational study using 3-phenylpropionate and 2-phenylbutyrate as substrates. We set up large cluster models containing the heme, the substrate and the key features of the substrate binding pocket and calculated (using density functional theory) the pathways leading to the four possible products. This work predicts that the two substrates will react with different reaction rates due to accessibility differences of the substrates to the active oxidant, and, as a consequence, these two substrates will also generate different products. This work explains how the substrate binding pocket of P450 OleTJE guides a reaction to a chemoselectivity.

Local Electric Fields as a Natural Switch of Heme-Iron Protein Reactivity

Heme-iron oxidoreductases operating through the high-valent Fe^IVO intermediates perform crucial and complicated transformations, such as oxidations of unreactive saturated hydrocarbons. These enzymes share the same Fe coordination, only differing by the axial ligation, e.g., Cys in P450 oxygenases, Tyr in catalases, and His in peroxidases. By examining ~200 heme-iron proteins, we show that the protein hosts exert highly specific intramolecular electric fields on the active sites, and there is a strong correlation between the direction and magnitude of this field and the protein function. In all heme proteins, the field is preferentially aligned with the Fe-O bond ( F_z ). The Cys-ligated P450 oxygenases have the highest average F_z of 28.5 MV cm^-1, i.e., most enhancing the oxyl-radical character of the oxo group, and consistent with the ability of these proteins to activate strong C-H bonds. In contrast, in Tyr-ligated proteins, the average F_z is only 3.0 MV cm^-1, apparently suppressing single-electron off-pathway oxidations, and in His-ligated proteins, F_z is -8.7 MV cm^-1. The operational field range is given by the trade-off between the low reactivity of the Fe^IVO Compound I at the more negative F_z , and the low selectivity at the more positive F_z . Consequently, a heme-iron site placed in the field characteristic of another heme-iron protein class loses its canonical function, and gains an adverse one. Thus, electric fields produced by the protein scaffolds, together with the nature of the axial ligand, control all heme-iron chemistry.

Histidine versus Cysteine-Bearing Heme-Dependent Halogen Peroxidases: Parallels and Differences for Cl- Oxidation

The homodimeric myeloperoxidase (MPO) features a histidine as a proximal ligand and a sulfonium linkage covalently attaching the heme porphyrin ring to the protein. MPO is able to catalyze Cl^- oxidation with about the same efficiency as chloroperoxidase at pH 7.0. In this study, we seek to explore the parallels and differences between the histidine and cysteine heme-dependent halogen peroxidases. Transition states, reaction barriers, and relevant thermodynamic properties are calculated on protein models. Together with electronic structure calculations, it gives an overview of the reaction mechanisms and of the factors that determine the selectivity between one- and two-electron paths. Conclusions point to the innate oxidizing nature of MPO with the ester and sulfonium linkages hiking up the reactivity to enable chloride oxidation. The installation of a deprotonated imidazolate as a proximal ligand does not shift the equilibrium from one- to two-electron events without influencing the chemistry of the oxidation reaction.

How Does Replacement of the Axial Histidine Ligand in Cytochrome c Peroxidase by Nδ-Methyl Histidine Affect Its Properties and Functions? A Computational Study

Heme peroxidases have important functions in nature related to the detoxification of H₂O₂. They generally undergo a catalytic cycle where, in the first stage, the iron(III)-heme-H₂O₂ complex is converted into an iron(IV)-oxo-heme cation radical species called Compound I. Cytochrome c peroxidase Compound I has a unique electronic configuration among heme enzymes where a metal-based biradical is coupled to a protein radical on a nearby Trp residue. Recent work using the engineered N_δ-methyl histidine-ligated cytochrome c peroxidase highlighted changes in spectroscopic and catalytic properties upon axial ligand substitution. To understand the axial ligand effect on structure and reactivity of peroxidases and their axially N_δ-methyl histidine engineered forms, we did a computational study. We created active site cluster models of various sizes as mimics of horseradish peroxidase and cytochrome c peroxidase Compound I. Subsequently, we performed density functional theory studies on the structure and reactivity of these complexes with a model substrate (styrene). Thus, the work shows that the N_δ-methyl histidine group has little effect on the electronic configuration and structure of Compound I and little changes in bond lengths and the same orbital occupation is obtained. However, the N_δ-methyl histidine modification impacts electron transfer processes due to a change in the reduction potential and thereby influences reactivity patterns for oxygen atom transfer. As such, the substitution of the axial histidine by N_δ-methyl histidine in peroxidases slows down oxygen atom transfer to substrates and makes Compound I a weaker oxidant. These studies are in line with experimental work on N_δ-methyl histidine-ligated cytochrome c peroxidases and highlight how the hydrogen bonding network in the second coordination sphere has a major impact on the function and properties of the enzyme.

QM Calculations in ADMET Prediction

In recent years, there has been an increase in the application of quantum mechanics (QM) methods to describe properties related to the ADMET profile of small molecules. The application of these methods allows calculating useful descriptors and physiochemical properties contributing to ADMET prediction. Considering that QM methods are the only one that describe the electronic state of a molecules, such methods are particularly useful for studying the metabolism of drugs; furthermore, the introduction of mixed QM and molecular mechanics (QM/MM) is also increasing the understanding of drug interaction with cytochromes from a mechanistic point of view. Finally, combining the increase number of experimental data with machine learning algorithms and QM-derived descriptors allowed the creation of an end-user software capable of affecting the drug discovery process.

Molecular Dynamics Simulation Framework to Probe the Binding Hypothesis of CYP3A4 Inhibitors

The Cytochrome P450 family of heme-containing proteins plays a major role in catalyzing phase I metabolic reactions, and the CYP3A4 subtype is responsible for the metabolism of many currently marketed drugs. Additionally, CYP3A4 has an inherent affinity for a broad spectrum of structurally diverse chemical entities, often leading to drug-drug interactions mediated by the inhibition or induction of the metabolic enzyme. The current study explores the binding of selected highly efficient CYP3A4 inhibitors by docking and molecular dynamics (MD) simulation protocols and their binding free energy calculated using the WaterSwap method. The results indicate the importance of binding pocket residues including Phe57, Arg105, Arg106, Ser119, Arg212, Phe213, Thr309, Ser312, Ala370, Arg372, Glu374, Gly481 and Leu483 for interaction with CYP3A4 inhibitors. The residue-wise decomposition of the binding free energy from the WaterSwap method revealed the importance of binding site residues Arg106 and Arg372 in the stabilization of all the selected CYP3A4-inhibitor complexes. The WaterSwap binding energies were further complemented with the MM(GB/PB)SA results and it was observed that the binding energies calculated by both methods do not differ significantly. Overall, our results could guide towards the use of multiple computational approaches to achieve a better understanding of CYP3A4 inhibition, subsequently leading to the design of highly specific and efficient new chemical entities with suitable ADMETox properties and reduced side effects.

The Nature and Reactivity of Ferryl Heme in Compounds I and II

Aerobic organisms have evolved to activate oxygen from the atmosphere, which allows them to catalyze the oxidation of different kinds of substrates. This activation of oxygen is achieved by a metal center (usually iron or copper) buried within a metalloprotein. In the case of iron-containing heme enzymes, the activation of oxygen is achieved by formation of transient iron-oxo (ferryl) intermediates; these intermediates are called Compound I and Compound II. The Compound I and II intermediates were first discovered in the 1930s in horseradish peroxidase, and it is now known that these same species are used across the family of heme enzymes, which include all of the peroxidases, the heme catalases, the P450s, cytochrome c oxidase, and NO synthase. Many years have passed since the first observations, but establishing the chemical nature of these transient ferryl species remains a fundamental question that is relevant to the reactivity, and therefore the usefulness, of these species in biology. This Account summarizes experiments that were conceived and conducted at Leicester and presents our ideas on the chemical nature, stability, and reactivity of these ferryl heme species. We begin by briefly summarizing the early milestones in the field, from the 1940s and 1950s. We present comparisons between the nature and reactivity of the ferryl species in horseradish peroxidase, cytochrome c peroxidase, and ascorbate peroxidase; and we consider different modes of electron delivery to ferryl heme, from different substrates in different peroxidases. We address the question of whether the ferryl heme is best formulated as an (unprotonated) Fe^IV═O or as a (protonated) Fe^IV-OH species. A range of spectroscopic approaches (EXAFS, resonance Raman, Mossbauer, and EPR) have been used over many decades to examine this question, and in the last ten years, X-ray crystallography has also been employed. We describe how information from all of these studies has blended together to create an overall picture, and how the recent application of neutron crystallography has directly identified protonation states and has helped to clarify the precise nature of the ferryl heme in cytochrome c peroxidase and ascorbate peroxidase. We draw comparisons between the Compound I and Compound II species that we have observed in peroxidases with those found in other heme systems, notably the P450s, highlighting possible commonality across these heme ferryl systems. The identification of proton locations from neutron structures of these ferryl species opens the door for understanding the proton translocations that need to occur during O-O bond cleavage.

Cytochrome c peroxidase regulates intracellular reactive oxygen species and methylglyoxal via enzyme activities of erythroascorbate peroxidase and glutathione-related enzymes in Candida albicans

D-erythroascorbate peroxidase (EAPX1) deficiency causes glutathione deprivation, leading to the accumulation of methylglyoxal and reactive oxygen species (ROS), and especially, induction of cytochrome c peroxidase (Ccp1) in Candida albicans. Nevertheless, reciprocal effects between changes in Ccp1 activity and the antioxidative D-erythroascorbic acid- and glutathione-dependent redox status, which reflects methylglyoxal biosynthesis altering pathophysiology are unclear in eukaryotes. To elucidate the effect of CCP1 expression on EAPX1 and glutathione reductase (Glr1) activity-mediated D-erythroascorbic acid biosynthesis and redox homeostasis, the CCP1 gene was disrupted and overexpressed. First, we demonstrated both glutathione-independent and-dependent metabolite contents and their corresponding gene transcripts and enzyme activities (i.e., Ccp1, catalase-peroxidase [KatG], superoxide dismutase [Sod], Eapx1, and Glr1) in CCP1 mutants. Second, methylglyoxal-oxidizing alcohol dehydrogenase (Adh1) and methylglyoxal-reducing oxidoreductase activity on glycolytic methylglyoxal and pyruvate production and NAD(P)H content were determined in these mutants. Contrary to our expectation, CCP1 disruption (42.19±3.22nmolO₂h^-1mgwetcell^-1) failed to affect cell respiration compared to the wild-type strain (41.62±7.11nmolO₂h^-1mgwetcell^-1) under cyanide treatment, and in contrast to hydrogen peroxide (H₂O₂) treatment (21.74±1.03nmol O₂h^-1mgwetcell^-1). Additionally, Ccp1 predominantly detoxified H₂O₂ rather than negligible scavenging activities towards methylglyoxal and other oxidants. CCP1 deficiency stimulated Sod and Adh1 activity but downregulated Glr1, Eapx1, catalase, and peroxidase activity while enhancing KatG, EAPX1, and GLR1 transcription by decreasing glutathione and D-erythroascorbic acid and increasing pyruvate. Noticeably, the ROS-accumulating CCP1-deficient mutant maintained steady-state levels of methylglyoxal, which was revealed to be regulated by methylglyoxal-oxidizing and -reducing activity with drastic changes in NAD(P)H. We confirmed and clarified our results by showing that CCP1/EAPX1 double disruptants underwent severe growth defects due to the D-erythroascorbic acid and glutathione depletion because of pyruvate overaccumulation. These observations were made in both budding and hyphal-growing CCP1 mutants. The revealed metabolic network involving Ccp1 and other redox regulators affected ROS and methylglyoxal through D-erythroascorbic acid and glutathione-dependent metabolites, thereby influencing dimorphism. This is the first report of the Ccp1-mediated D-erythroascorbic acid and glutathione biosynthesis accompanying methylglyoxal scavengers for full fungal virulence.

Construction and in vivo assembly of a catalytically proficient and hyperthermostable de novo enzyme

Although catalytic mechanisms in natural enzymes are well understood, achieving the diverse palette of reaction chemistries in re-engineered native proteins has proved challenging. Wholesale modification of natural enzymes is potentially compromised by their intrinsic complexity, which often obscures the underlying principles governing biocatalytic efficiency. The maquette approach can circumvent this complexity by combining a robust de novo designed chassis with a design process that avoids atomistic mimicry of natural proteins. Here, we apply this method to the construction of a highly efficient, promiscuous, and thermostable artificial enzyme that catalyzes a diverse array of substrate oxidations coupled to the reduction of H₂O₂. The maquette exhibits kinetics that match and even surpass those of certain natural peroxidases, retains its activity at elevated temperature and in the presence of organic solvents, and provides a simple platform for interrogating catalytic intermediates common to natural heme-containing enzymes.Catalytic mechanisms of enzymes are well understood, but achieving diverse reaction chemistries in re-engineered proteins can be difficult. Here the authors show a highly efficient and thermostable artificial enzyme that catalyzes a diverse array of substrate oxidations coupled to the reduction of H₂O₂.

How do Enzymes Utilize Reactive OH Radicals? Lessons from Nonheme HppE and Fenton Systems

The iron(IV)-oxo (ferryl) intermediate has been amply established as the principal oxidant in nonheme enzymes and the key player in C-H bond activations and functionalizations. In contrast to this status, our present QM/MM calculations of the mechanism of fosfomycin biosynthesis (a broad range antibiotic) by the nonheme HppE enzyme rule out the iron(IV)-oxo as the reactive species in the hydrogen abstraction (H-abstraction) step of the pro-R hydrogen from the (S)-2-hydroxypropylphosphonic substrate. Moreover, the study reveals that the ferryl species is bypassed in HppE, while the actual oxidant is an HO(•) radical hydrogen-bonded to a ferric-hydroxo complex, resulting via the homolytic dissociation of the hydrogen peroxide complex, Fe(II)-H2O2. The computed energy barrier of this pathway is 12.0 kcal/mol, in fair agreement with the experimental datum of 9.8 kcal/mol. An alternative mechanism involves the iron-complexed hydroxyl radical (Fe(III)-(HO(•))) that is obtained by protonation of the iron(IV)-oxo group via the O-H group of the substrate. The barrier for this pathway, 23.0 kcal/mol, is higher than the one in the first mechanism. In both mechanisms, the HO(•) radical is highly selective; its H-abstraction leading to the final cis-fosfomycin product. It appears that HppE is prone to usage of HO(•) radicals for C-H bond activation, because the ferryl oxidant requires a specific H-abstraction trajectory (∠FeOH ∼ 180°) that cannot be met for intramolecular H-abstraction. Thus, this work broadens the landscape of nonheme iron enzymes and makes a connection to Fenton chemistry, with implications on new potential biocatalysts that may harness hydroxyl radicals for C-H bond functionalizations.

Applications of density functional theory to iron-containing molecules of bioinorganic interest

The past decades have seen an explosive growth in the application of density functional theory (DFT) methods to molecular systems that are of interest in a variety of scientific fields. Owing to its balanced accuracy and efficiency, DFT plays particularly useful roles in the theoretical investigation of large molecules. Even for biological molecules such as proteins, DFT finds application in the form of, e.g., hybrid quantum mechanics and molecular mechanics (QM/MM), in which DFT may be used as a QM method to describe a higher prioritized region in the system, while a MM force field may be used to describe remaining atoms. Iron-containing molecules are particularly important targets of DFT calculations. From the viewpoint of chemistry, this is mainly because iron is abundant on earth, iron plays powerful (and often enigmatic) roles in enzyme catalysis, and iron thus has the great potential for biomimetic catalysis of chemically difficult transformations. In this paper, we present a brief overview of several recent applications of DFT to iron-containing non-heme synthetic complexes, heme-type cytochrome P450 enzymes, and non-heme iron enzymes, all of which are of particular interest in the field of bioinorganic chemistry. Emphasis will be placed on our own work.

Effect of protein environment within cytochrome P450cam evaluated using a polarizable-embedding QM/MM method

Metalloenzymes accommodate cofactors and substrates in their active sites, thereby exerting powerful catalytic effects. Understanding the key elements of the mechanism via which such binding is accomplished using a number of atoms in a protein is a fundamental challenge. To address this issue computationally, here we used mechanical-embedding (ME), electronic-embedding (EE), and polarizable-embedding (PE) hybrid quantum mechanics and molecular mechanics (QM/MM) methods and performed an energy decomposition analysis (EDA) of the nonbonding protein environmental effect in the "compound I" intermediate state of cytochrome P450cam. The B3LYP and AMBER99/QP302 methods were used to deal with the QM and MM subsystems, respectively, and the nonbonding interaction energy between these subsystems was decomposed into electrostatic, van der Waals, and polarization contributions. The PE-QM/MM calculation was performed using polarizable force fields that were capable of describing induced dipoles within the MM subsystem, which arose in response to the electric field generated by QM electron density, QM nuclei, and MM point charges. The present QM/MM EDA revealed that the electrostatic term constituted the largest stabilizing interaction between the QM and MM subsystems. When proper adjustment was made for the point charges of the MM atoms located at the QM-MM boundary, EE-QM/MM and PE-QM/MM calculations yielded similar QM electron density distributions, indicating that the MM polarization effect does not have a significant influence on the extent of QM polarization in this particular enzyme system.

Combined quantum mechanics/molecular mechanics (QM/MM) methods in computational enzymology

Computational enzymology is a rapidly maturing field that is increasingly integral to understanding mechanisms of enzyme-catalyzed reactions and their practical applications. Combined quantum mechanics/molecular mechanics (QM/MM) methods are important in this field. By treating the reacting species with a quantum mechanical method (i.e., a method that calculates the electronic structure of the active site) and including the enzyme environment with simpler molecular mechanical methods, enzyme reactions can be modeled. Here, we review QM/MM methods and their application to enzyme-catalyzed reactions to investigate fundamental and practical problems in enzymology. A range of QM/MM methods is available, from cheaper and more approximate methods, which can be used for molecular dynamics simulations, to highly accurate electronic structure methods. We discuss how modeling of reactions using such methods can provide detailed insight into enzyme mechanisms and illustrate this by reviewing some recent applications. We outline some practical considerations for such simulations. Further, we highlight applications that show how QM/MM methods can contribute to the practical development and application of enzymology, e.g., in the interpretation and prediction of the effects of mutagenesis and in drug and catalyst design.

Effects of Dispersion in Density Functional Based Quantum Mechanical/Molecular Mechanical Calculations on Cytochrome P450 Catalyzed Reactions

Density functional theory (DFT) based quantum mechanical/molecular mechanical (QM/MM) calculations have provided valuable insight into the reactivity of the cytochrome P450 family of enzymes (P450s). A failure of commonly used DFT methods, such as B3LYP, is the neglect of dispersion interactions. An empirical dispersion correction has been shown to improve the accuracy of gas phase DFT calculations of P450s. The current work examines the effect of the dispersion correction in QM/MM calculations on P450s. The hydrogen abstraction from camphor, and hydrogen abstraction and C-O addition of cyclohexene and propene by P450cam have been modeled, along with the addition of benzene to Compound I in CYP2C9, at the B3LYP-D2/CHARMM27 level of theory. Single point energy calculations were also performed at the B3LYP-D3//B3LYP-D2/CHARMM27 level. The dispersion corrections lower activation energy barriers significantly (by ∼5 kcal/mol), as seen for gas phase calculations, but has a small effect on optimized geometries.These effects are likely to be important in modeling reactions catalyzed by other enzymes also. Given the low computational cost of including such dispersion corrections, we recommend doing so in all B3LYP based QM/MM calculations.

A practical guide to modelling enzyme-catalysed reactions

Molecular modelling and simulation methods are increasingly at the forefront of elucidating mechanisms of enzyme-catalysed reactions, and shedding light on the determinants of specificity and efficiency of catalysis. These methods have the potential to assist in drug discovery and the design of novel protein catalysts. This Tutorial Review highlights some of the most widely used modelling methods and some successful applications. Modelling protocols commonly applied in studying enzyme-catalysed reactions are outlined here, and some practical implications are considered, with cytochrome P450 enzymes used as a specific example.

Does compound I vary significantly between isoforms of cytochrome P450?

The cytochrome P450 (CYP) enzymes are important in many areas, including pharmaceutical development. Subtle changes in the electronic structure of the active species, Compound I, have been postulated previously to account partly for the experimentally observed differences in reactivity between isoforms. Current predictive models of CYP metabolism typically assume an identical Compound I in all isoforms. Here we present a method to calculate the electronic structure and to estimate the Fe-O bond enthalpy of Compound I, and apply it to several human and bacterial CYP isoforms. Conformational flexibility is accounted for by sampling large numbers of structures from molecular dynamics simulations, which are subsequently optimized with density functional theory (B3LYP) based quantum mechanics/molecular mechanics. The observed differences in Compound I between human isoforms are small: They are generally smaller than the spread of values obtained for the same isoform starting from different initial structures. Hence, it is unlikely that the variation in activity between human isoforms is due to differences in the electronic structure of Compound I. A larger difference in electronic structure is observed between the human isoforms and P450(cam) and may be explained by the slightly different hydrogen-bonding environment surrounding the cysteinyl sulfur. The presence of substrate in the active site of all isoforms studied appears to cause a slight decrease in the Fe-O bond enthalpy, apparently due to displacement of water out of the active site, suggesting that Compound I is less stable in the presence of substrate.

Multireference and multiconfiguration ab initio methods in heme-related systems: what have we learned so far?

This work reviews the recent applications of ab initio multireference/multiconfiguration (MR/MC) electronic structure methods to heme-related systems, involving tetra-, penta-, and hexa-coordinate species, as well as the high-valent iron-oxo species. The current accuracy of these methods in the various systems is discussed, with special attention to potential sources of systematic errors. Thus, the review summarizes and tries to rationalize the key elements of MR/MC calculations, namely, the choice of the employed active space, especially the so-called double-shell effect that has already been recognized to be important in transition-metal-containing systems, and the impact of these elements on the spin-state energetics of heme species, as well as on the bonding mechanism of small molecules to the heme. It is shown that expansion of the MC wave function into one based on localized orbitals provides a compact and insightful view on some otherwise complex electronic structures. The effects of protein environment on the MR/MC results are summarized for the few available quantum mechanical/molecular mechanical (QM/MM) studies. Comparisons with corresponding DFT results are also made wherever available. Potential future directions are proposed.

QM/MM methods: looking inside heme proteins biochemistry

Mixed quantum mechanics/molecular mechanics methods offer a valuable computational tool for understanding biochemical events. When combined with conformational sampling techniques, they allow for an exhaustive exploration of the enzymatic mechanism. Heme proteins are ubiquitous and essential for every organism. In this review we summarize our efforts towards the understanding of heme biochemistry. We present: 1) results on ligand migration on globins coupled to the ligand binding event, 2) results on the localization of the spin density in compound I of cytochromes and peroxidases, 3) novel methodologies for mapping the electron transfer pathways and 4) novel data on Tryptophan 2,3-dioxygenase. For this enzyme our results strongly indicate that the distal oxygen will end up on the C3 indole carbon, whereas the proximal oxygen will end up in the C2 position. Interestingly, the process involves the formation of an epoxide and a heme ferryl intermediate. The overall energy profile indicates an energy barrier of approximately 18 kcal/mol and an exothermic driving force of almost 80 kcal/mol.

New insights into the structural characteristics and functional relevance of the human cytochrome P450 2D6 enzyme

To date, the crystal structures of at least 12 human CYPs (1A2, 2A6, 2A13, 2C8, 2C9, 2D6, 2E1, 2R1, 3A4, 7A1, 8A1, and 46A1) have been determined. CYP2D6 accounts for only a small percentage of all hepatic CYPs (< 2%), but it metabolizes approximately 25% of clinically used drugs with significant polymorphisms. CYP2D6 also metabolizes procarcinogens and neurotoxins, such as 1-methyl-4-phenyl-1,2,3,6-tetrahydropyridine, 1,2,3,4-tetrahydroquinoline, and indolealkylamines. Moreover, the enzyme utilizes hydroxytryptamines and neurosteroids as endogenous substrates. Typical CYP2D6 substrates are usually lipophilic bases with an aromatic ring and a nitrogen atom, which can be protonated at physiological pH. Substrate binding is generally followed by oxidation (5-7 A) from the proposed nitrogen-Asp301 interaction. A number of homology models have been constructed to explore the structural features of CYP2D6, while antibody studies also provide useful structural information. Site-directed mutagenesis studies have demonstrated that Glu216, Asp301, Phe120, Phe481, and Phe483 play important roles in determining the binding of ligands to CYP2D6. The structure of human CYP2D6 has been recently determined and shows the characteristic CYP fold observed for other members of the CYP superfamily. The lengths and orientations of the individual secondary structural elements in the CYP2D6 structure are similar to those seen in other human CYP2 members, such as CYP2C9 and 2C8. The 2D6 structure has a well-defined active-site cavity located above the heme group with a volume of approximately 540 A(3), which is larger than equivalent cavities in CYP2A6 (260 A(3)), 1A2 (375 A(3)), and 2E1 (190 A(3)), but smaller than those in CYP3A4 (1385 A(3)) and 2C8 (1438 A(3)). Further studies are required to delineate the molecular mechanisms involved in CYP2D6 ligand interactions and their implications for drug development and clinical practice.

QM/MM methods for biomolecular systems

Combined quantum-mechanics/molecular-mechanics (QM/MM) approaches have become the method of choice for modeling reactions in biomolecular systems. Quantum-mechanical (QM) methods are required for describing chemical reactions and other electronic processes, such as charge transfer or electronic excitation. However, QM methods are restricted to systems of up to a few hundred atoms. However, the size and conformational complexity of biopolymers calls for methods capable of treating up to several 100,000 atoms and allowing for simulations over time scales of tens of nanoseconds. This is achieved by highly efficient, force-field-based molecular mechanics (MM) methods. Thus to model large biomolecules the logical approach is to combine the two techniques and to use a QM method for the chemically active region (e.g., substrates and co-factors in an enzymatic reaction) and an MM treatment for the surroundings (e.g., protein and solvent). The resulting schemes are commonly referred to as combined or hybrid QM/MM methods. They enable the modeling of reactive biomolecular systems at a reasonable computational effort while providing the necessary accuracy.

Mapping protein electron transfer pathways with QM/MM methods

Mixed quantum mechanics/molecular mechanics (QM/MM) methods offer a valuable computational tool for understanding the electron transfer pathway in protein-substrate interactions and protein-protein complexes. These hybrid methods are capable of solving the Schrödinger equation on a small subset of the protein, the quantum region, describing its electronic structure under the polarization effects of the remainder of the protein. By selectively turning on and off different residues in the quantum region, we are able to obtain the electron pathway for short- and large-range interactions. Here, we summarize recent studies involving the protein-substrate interaction in cytochrome P450 camphor, ascorbate peroxidase and cytochrome c peroxidase, and propose a novel approach for the long-range protein-protein electron transfer. The results on ascorbate peroxidase and cytochrome c peroxidase reveal the importance of the propionate groups in the electron transfer pathway. The long-range protein-protein electron transfer has been studied on the cytochrome c peroxidase-cytochrome c complex. The results indicate the importance of Phe82 and Cys81 on cytochrome c, and of Asn196, Ala194, Ala176 and His175 on cytochrome c peroxidase.

QM/MM study of the active species of the human cytochrome P450 3A4, and the influence thereof of the multiple substrate binding

Cytochrome P450 3A4 is involved in the metabolism of 50% of all swallowed drugs. The enzyme functions by means of a high-valent iron-oxo species, called compound I (Cpd I), which is formed after entrance of the substrate to the active site. We explored the features of Cpd I using hybrid quantum mechanical/molecular mechanical calculations on various models that are either substrate-free or containing one and two molecules of diazepam as a substrate. Mössbauer parameters of Cpd I were computed. Our major finding shows that without the substrate, Cpd I tends to elongate its Fe-S bond, localize the radical on the sulfur, and form hydrogen bonds with A305 and T309, which may hypothetically lead to Cpd I consumption by H-abstraction. However, the positioning of diazepam close to Cpd I, as enforced by the effector molecule, was found to strengthen the NH...S interactions of the conserved I443 and G444 residues with the proximal cysteinate ligand. These interactions are known to stabilize the Fe-S bond, and as such, the presence of the substrate leads to a shorter Fe-S bond and it prevents the localization of the radical on the sulfur. This diazepam-Cpd I stabilization was manifested in the 1W0E conformer. The effector substrate did not influence Cpd I directly but rather by positioning the active substrate close to Cpd I, thus displacing the hydrogen bonds with A305 and T309, and thereby giving preference to substrate oxidation. It is hypothesized that these effects on Cpd I, promoted by the restrained substrate, may be behind the special metabolic behavior observed in cases of multiple substrate binding (also called cooperative binding). This restraint constitutes a mechanism whereby substrates stabilize Cpd I sufficiently long to affect monooxygenation by P450s at the expense of Cpd I destruction by the protein residues.

What is the Active Species of Cytochrome P450 during Camphor Hydroxylation? QM/MM Studies of Different Electronic States of Compound I and of Reduced and Oxidized Iron−Oxo Intermediates

We have investigated C-H hydroxylation of camphor by Compound I (Cpd I) of cytochrome P450cam in different electronic states and by its one-electron reduced and oxidized forms, using QM/MM calculations in the native protein/solvent environment. Cpd I species with five unpaired electrons (pentaradicaloids) are ca. 12 kcal/mol higher in energy than the ground state Cpd I species with three unpaired electrons (triradicaloids). The H-abstraction transition states of pentaradicaloids lie ca. 21 (9) kcal/mol above the triradicaloid (pentaradicaloid) reactants. Hydroxylation via pentaradicaloids is thus facile provided that they can react before relaxing to the ground-state triradicaloids. Excited states of Cpd I with an Fe(V)-oxo moiety lie more than 20 kcal/mol above the triradicaloid ground state in single-point gas-phase calculations, but these electronic configurations are not stable upon including the point-charge protein environment which causes SCF convergence to the triradicaloid ground state. One-electron reduced species (Cpd II) show sluggish reactivity compared with Cpd I in agreement with experimental model studies. One-electron oxidized species are more reactive than Cpd I but seem too high in energy to be accessible. The barriers to hydrogen abstraction for the various forms of Cpd I are generally not affected much by the chosen protonation states of the Asp297 and His355 residues near the propionate side chains of the heme or by the appearance of radical character at Asp297, His355, or the propionates.

Reaction mechanisms of 15-hydroperoxyeicosatetraenoic acid catalyzed by human prostacyclin and thromboxane synthases

Prostacyclin synthase (PGIS) and thromboxane synthase (TXAS) are atypical cytochrome P450s. They do not require NADPH or dioxygen for isomerization of prostaglandin H(2) (PGH(2)) to produce prostacyclin (PGI(2)) and thromboxane A(2) (TXA(2)). PGI(2) and TXA(2) have opposing actions on platelet aggregation and blood vessel tone. In this report, we use a lipid hydroperoxide, 15-hydroperoxyeicosatetraenoic acid (15-HPETE), to explore the active site characteristics of PGIS and TXAS. The two enzymes transformed 15-HPETE not only into 13-hydroxy-14,15-epoxy-5,8,11-eicosatrienoic acid (13-OH-14,15-EET), like many microsomal P450s, but also to 15-ketoeicosatetraenoic acid (15-KETE) and 15-hydroxyeicosatetraenoic acid (15-HETE). 13-OH-14,15-EET and 15-KETE result from homolytic cleavage of the O-O bond, whereas 15-HETE results from heterolytic cleavage, a common peroxidase pathway. About 80% of 15-HPETE was homolytically cleaved by PGIS and 60% was homolytically cleaved by TXAS. The V(max) of homolytic cleavage is 3.5-fold faster than heterolytic cleavage for PGIS-catalyzed reactions (1100 min(-1)vs. 320 min(-1)) and 1.4-fold faster for TXAS (170 min(-1)vs. 120 min(-1)). Similar K(M) values for homolytic and heterolytic cleavages were found for PGIS ( approximately 60 microM 15-HPETE) and TXAS ( approximately 80 microM 15-HPETE), making PGIS a more efficient catalyst for the 15-HPETE reaction.

Ab initio QM/MM modelling of acetyl-CoA deprotonation in the enzyme citrate synthase

The first step of the reaction catalysed by the enzyme citrate synthase is studied here with high level combined quantum mechanical/molecular mechanical (QM/MM) methods (up to the MP2/6-31+G(d)//6-31G(d)/CHARMM level). In the first step of the reaction, acetyl-CoA is deprotonated by Asp375, producing an intermediate, which is the nucleophile for attack on the second substrate, oxaloacetate, prior to hydrolysis of the thioester bond of acetyl-CoA and release of the products. A central question has been whether the nucleophilic intermediate is the enolate of acetyl-CoA, the enol, or an 'enolic' intermediate stabilized by a 'low-barrier' hydrogen bond with His274 at the active site. The imidazole sidechain of His274 is neutral, and donates a hydrogen bond to the carbonyl oxygen of acetyl-CoA in substrate complexes. We have investigated the identity of the nucleophilic intermediate by QM/MM calculations on the substrate (keto), enolate, enol and enolic forms of acetyl-CoA at the active site of citrate synthase. The transition states for proton abstraction from acetyl-CoA by Asp375, and for transfer of the hydrogen bonded proton between His274 and acetyl-CoA have been modelled approximately. The effects of electron correlation are included by MP2/6-31G(d) and MP2/6-31+G(d) calculations on active site geometries produced by QM/MM energy minimization. The results do not support the hypothesis that a low-barrier hydrogen bond is involved in catalysis in citrate synthase, in agreement with earlier calculations. The acetyl-CoA enolate is identified as the only intermediate consistent with the experimental barrier for condensation, stabilized by conventional hydrogen bonds from His274 and a water molecule.

Common system setup for the entire catalytic cycle of cytochrome P450cam in quantum mechanical/molecular mechanical studies

We describe a system setup that is applicable to all species in the catalytic cycle of cytochrome P450(cam). The chosen procedure starts from the X-ray coordinates of the ferrous dioxygen complex and follows a protocol that includes the careful assignment of protonation states, comparison between different conceivable hydration schemes, and system preparation through a series of classical minimizations and molecular dynamics (MD) simulations. The resulting setup was validated by quantum mechanical/molecular mechanical (QM/MM) calculations on the resting state, the pentacoordinated ferric and ferrous complexes, Compound I, the transition state and hydroxo intermediate of the C--H hydroxylation reaction, and the product complex. The present QM/MM results are generally consistent with those obtained previously with individual setups. Concerning hydration, we find that saturating the protein interior with water is detrimental and leads to higher structural flexibility and catalytically inefficient active-site geometries. The MD simulations favor a low water density around Asp251 that facilitates side chain rotation of protonated Asp251 during the conversion of Compound 0 to Compound I. The QM/MM results for the two preferred hydration schemes (labeled SE-1 and SE-4) are similar, indicating that slight differences in the solvation close to the active site are not critical as long as camphor and the crystallographic water molecules preserve their positions in the experimental X-ray structures.

Mechanisms of reaction in cytochrome P450: Hydroxylation of camphor in P450cam

The fundamental nature of reactivity in cytochrome P450 enzymes is currently controversial. Modelling of bacterial P450cam has suggested an important role for the haem propionates in the catalysis, though this finding has been questioned. Understanding the mechanisms of this enzyme family is important both in terms of basic biochemistry and potentially in the prediction of drug metabolism. We have modelled the hydroxylation of camphor by P450cam, using combined quantum mechanics/molecular mechanics (QM/MM) methods. A set of reaction pathways in the enzyme was determined. We were able to pinpoint the source of the discrepancies in the previous results. We show that when a correct ionization state is assigned to Asp297, no spin density appears on the haem propionates and the protein structure in this region remains preserved. These results indicate that the haem propionates are not involved in catalysis.