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

Hydrogen-Bond-Assisted Catalysis: Hydroxylation of Paclitaxel by Human CYP2C8

Paclitaxel (PTX, or Taxol), a chemotherapeutic agent widely employed in the treatment of various cancers, undergoes metabolic transformations through the cytochrome P450 enzymes CYP3A4 and CYP2C8. CYP3A4 catalyzes the aromatic hydroxylation reaction of PTX, whereas CYP2C8 demonstrates a distinct reactivity pattern, producing 6α-hydroxypaclitaxel via alkane hydroxylation. Despite the significant impact of PTX metabolism on its anticancer efficacy, the detailed mechanisms underlying these transformations have remained largely unclear. In this study, we employed hybrid quantum mechanics and molecular mechanics (QM/MM) calculations to elucidate the mechanism of PTX metabolism by human CYP2C8. Our QM/MM results reveal that the hydroxylation of PTX by CYP2C8 follows an atypical rebound mechanism. Either of the two hydrogen atoms at the C6 position of PTX can be abstracted, leading to a common radical intermediate. Although the subsequent rebound barrier is unusually high, stereochemical scrambling is unlikely, as the rebound barrier for the formation of the 6α-hydroxylated PTX─the actual product─is significantly lower than that for the 6β-hydroxylated metabolite. Thus, product selectivity is determined by the non-rate-determining rebound step. Furthermore, the hydroxyl group at the C7 position of PTX plays a catalytic role by facilitating the hydrogen abstraction and rebound steps. Our study also confirms a pronounced stability of the transition state in the high-spin sextet spin state, enabled by the enzyme's specific substrate positioning.

Exploring the selectivity of cytochrome P450 for enhanced novel anticancer agent synthesis

Cytochrome P450 monooxygenases are an extensive and unique class of enzymes, which can regio- and stereo-selectively functionalise hydrocarbons by way of oxidation reactions. These enzymes are naturally occurring but have also been extensively applied in a synthesis context, where they are used as efficient biocatalysts. Recently, a biosynthetic pathway where a cytochrome P450 monooxygenase catalyses a critical step of the pathway was uncovered, leading to the production of a number of products that display high antitumour potency. In this work, we use computational techniques to gain insight into the factors that determine the relative yields of the different products. We use conformational search algorithms to understand the substrate stereochemistry. On a machine-learned 3D protein structure, we use molecular docking to obtain a library of favourable poses for substrate-protein interaction. With molecular dynamics, we investigate the most favourable poses for reactivity on a molecular level, allowing us to investigate which protein-substrate interactions favour a given product and thus gain insight into the product selectivity.

Identification of volatile biomarkers for lung cancer from different histological sources: A comprehensive study

The identification of noninvasive volatile biomarkers for lung cancer is a significant clinical challenge. Through in vitro studies, the recognition of altered metabolism in cell volatile organic compound (VOC) emitting profile, along with the occurrence of oncogenesis, provides insight into the biochemical pathways involved in the production and metabolism of lung cancer volatile biomarkers. In this research, for the first time, a comprehensive comparative analysis of the volatile metabolites in NSCLS cells (A549), SCLC cells (H446), lung normal cells (BEAS-2B), as well as metabolites in both the oxidative stress (OS) group and control group. Specifically, the combination of eleven VOCs, including n-dodecane, acetaldehyde, isopropylbenzene, p-ethyltoluene and cis-1,3-dichloropropene, exhibited potential as volatile biomarkers for lung cancer originating from two different histological sources. Furthermore, the screening process in A549 cell lines resulted in the identification of three exclusive biomarkers, isopropylbenzene, formaldehyde and bromoform. Similarly, the exclusive biomarkers 1,2,4-trimethylbenzene, p-ethyltoluene, and cis-1,3-dichloropropene were present in the H446 cell line. Additionally, significant changes in trans-2-pentene, acetaldehyde, 1,2,4-trimethylbenzene, and bromoform were observed, indicating a strong association with OS. These findings highlight the potential of volatile biomarkers profiling as a means of noninvasive identification for lung cancer diagnosis.

Tyrosyl Radical-Mediated Sequential Oxidative Decarboxylation of Coproporphyrinogen III through PCET: Theoretical Insights into the Mechanism of Coproheme Decarboxylase ChdC

The peroxide-dependent coproheme decarboxylase ChdC from Geobacillus stearothermophilus catalyzes two key steps in the synthesis of heme b, i.e., two sequential oxidative decarboxylations of coproporphyrinogen III (coproheme III) at propionate groups P2 and P4. In the binding site of coproheme III, P2 and P4 are anchored by different residues (Tyr144, Arg217, and Ser222 for P2 and Tyr113, Lys148, and Trp156 for P4); however, strong experimental evidence supports that the generated Tyr144 radical acts as an unique intermediary for hydrogen atom transfer (HAT) from both reactive propionates. So far, the reaction details are still unclear. Herein, we carried out quantum mechanics/molecular mechanics calculations to explore the decarboxylation mechanism of coproheme III. In our calculations, the coproheme Cpd I, Fe(IV) = O coupled to a porphyrin radical cation (por^•+) with four propionate groups, was used as a reactant model. Our calculations reveal that Tyr144 is directly involved in the decarboxylation of propionate group P2. First, the proton-coupled electron transfer (PCET) occurs from Tyr144 to P2, generating a Tyr144 radical, which then abstracts a hydrogen atom from the Cβ of P2. The β-H extraction was calculated to be the rate-limiting step of decarboxylation. It is the porphyrin radical cation (por^•+) that makes the PCET from Tyr144 to P2 to be quite easy to initiate the decarboxylation. Finally, the electron transfers from the Cβ^• through the porphyrin to the iron center, leading to the decarboxylation of P2. Importantly, the decarboxylation of P4 mediated by Lys148 was calculated to be very difficult, which suggests that after the P2 decarboxylation, the generated harderoheme III intermediate should rebind or rotate in the active site so that the propionate P4 occupies the binding site of P2, and Tyr144 again mediates the decarboxylation of P4. Thus, our calculations support the fact that Tyr144 is responsible for the decarboxylation of both P2 and P4.

Reduction Potentials of P450 Compounds I and II: Insight into the Thermodynamics of C-H Bond Activation

We present a mixed experimental/theoretical determination of the bond strengths and redox potentials that define the ground-state thermodynamics for C-H bond activation in cytochrome P450 catalysis. Using redox titrations with [Ir(IV)Cl₆]^2-, we have determined the compound II/ferric (or Fe(IV)OH/Fe(III)OH₂) couple and its associated D(O-H)_Ferric bond strength in CYP158. Knowledge of this potential as well as the compound II/ferric (or Fe(IV)O/Fe(III)OH) reduction potential in horseradish peroxidase and the two-electron compound I/ferric (or Fe(IV)O(Por^•)/Fe(III)OH₂(Por)) reduction potential in aromatic peroxidase has allowed us to gauge the accuracy of theoretically determined bond strengths. Using the restricted open shell (ROS) method as proposed by Wright and co-workers, we have obtained O-H bond strengths and associated redox potentials for charge-neutral H-atom reductions of these iron(IV)-hydroxo and -oxo porphyrin species that are within 1 kcal/mol of experimentally determined values, suggesting that the ROS method may provide accurate values for the P450-II O-H bond strength and P450-I reduction potential. The efforts detailed here indicate that the ground-state thermodynamics of C-H bond activation in P450 are best described as follows: E^0'_Comp-I = 1.22 V (at pH 7, vs NHE) with D(O-H)_Comp-II = 95 kcal/mol and E^0'_Comp-II = 0.99 V (at pH 7, vs NHE) with D(O-H)_Ferric = 90 kcal/mol.

Biochemical Characterization of the Cytochrome P450 CYP107CB2 from Bacillus lehensis G1

The bioconversion of vitamin D3 catalyzed by cytochrome P450 (CYP) requires 25-hydroxylation and subsequent 1α-hydroxylation to produce the hormonal activated 1α,25-dihydroxyvitamin D3. Vitamin D3 25-hydroxylase catalyses the first step in the vitamin D3 biosynthetic pathway, essential in the de novo activation of vitamin D3. A CYP known as CYP107CB2 has been identified as a novel vitamin D hydroxylase in Bacillus lehensis G1. In order to deepen the understanding of this bacterial origin CYP107CB2, its detailed biological functions as well as biochemical characteristics were defined. CYP107CB2 was characterized through the absorption spectral analysis and accordingly, the enzyme was assayed for vitamin D3 hydroxylation activity. CYP-ligand characterization and catalysis optimization were conducted to increase the turnover of hydroxylated products in an NADPH-regenerating system. Results revealed that the over-expressed CYP107CB2 protein was dominantly cytosolic and the purified fraction showed a protein band at approximately 62 kDa on SDS-PAGE, indicative of CYP107CB2. Spectral analysis indicated that CYP107CB2 protein was properly folded and it was in the active form to catalyze vitamin D3 reaction at C25. HPLC and MS analysis from a reconstituted enzymatic reaction confirmed the hydroxylated products were 25-hydroxyitamin D3 and 1α,25-dihydroxyvitamin D3 when the substrates vitamin D3 and 1α-hydroxyvitamin D3 were used. Biochemical characterization shows that CYP107CB2 performed hydroxylation activity at 25 °C in pH 8 and successfully increased the production of 1α,25-dihydroxyvitamin D3 up to four fold. These findings show that CYP107CB2 has a biologically relevant vitamin D3 25-hydroxylase activity and further suggest the contribution of CYP family to the metabolism of vitamin D3.

Hydrogen Abstraction of Camphor Catalyzed by Cytochrome P450cam: A QM/MM Study

A combined quantum mechanics and molecular mechanics (QM/MM, QM = UB3LYP-D3, MM = AMBER) method is used to study the hydrogen abstraction reaction in P450_cam catalyzed hydroxylation of camphor in the quartet state. Compared to QM/MM calculations in the literature, this study uses larger basis sets for the most important atoms at the active site and QM/MM Hessian harmonic frequency calculations to determine the standard Gibbs free energy of activation and kinetic isotope effect. The QM/MM covalent boundary is treated with a capping hydrogen atom method, which is simple and robust. An energy barrier of 21.3 kcal/mol and a standard free energy of activation of 16.8 kcal/mol are obtained for this hydrogen abstraction reaction. These values are similar to those reported in the literature, suggesting that when a general protocol is followed, QM/MM results are reproducible. It is found that using a sufficiently large basis set is important to minimize basis set errors.

Molecular characterization, modeling and docking of CYP107CB2 from Bacillus lehensis G1, an alkaliphile

Cytochrome P450s are a superfamily of heme monooxygenases which catalyze a wide range of biochemical reactions. The reactions involve the introduction of an oxygen atom into an inactivated carbon of a compound which is essential to produce an intermediate of a hydroxylated product. The diversity of chemical reactions catalyzed by cytochrome P450s has led to their increased demand in numerous industrial and biotechnology applications. A recent study showed that a gene sequence encoding a CYP was found in the genome of Bacillus lehensis G1, and this gene shared structural similarity with the bacterial vitamin D hydroxylase (Vdh) from Pseudonocardia autotrophica. The objectives of present study was to mine, for a novel CYP from a new isolate B. lehensis G1 alkaliphile and determine the biological properties and functionalities of CYP in this bacterium. Our study employed the usage of computational methods to search for the novel CYP from CYP structural databases to identify the conserved pattern, functional domain and sequence properties of the uncharacterized CYP from B. lehensis G1. A computational homology model of the protein's structure was generated and a docking analysis was performed to provide useful structural knowledge on the enzyme's possible substrate and their interaction. Sequence analysis indicated that the newly identified CYP, termed CYP107CB2, contained the fingerprint heme binding sequence motif FxxGxxxCxG at position 336-345 as well as other highly conserved motifs characteristic of cytochrome P450 proteins. Using docking studies, we identified Ser-79, Leu-81, Val-231, Val-279, Val-383, Ala-232, Thr-236 and Thr-283 as important active site residues capable of stabilizing interactions with several potential substrates, including vitamin D3, 25-hydroxyvitamin D3 and 1α-hydroxyvitamin D3, in which all substrates docked proximally to the enzyme's heme center. Biochemical analysis indicated that CYP107CB2 is a biologically active protein to produce 1α,25-dihydroxyvitamin D3 from 1α-hydroxyvitamin D3. Based on these results, we conclude that the novel CYP107CB2 identified from B. lehensis G1 is a putative vitamin D hydroxylase which is possibly capable of catalyzing the bioconversion of parental vitamin D3 to calcitriol, or related metabolic products.

A multiscale approach to modelling drug metabolism by membrane-bound cytochrome P450 enzymes

Cytochrome P450 enzymes are found in all life forms. P450s play an important role in drug metabolism, and have potential uses as biocatalysts. Human P450s are membrane-bound proteins. However, the interactions between P450s and their membrane environment are not well-understood. To date, all P450 crystal structures have been obtained from engineered proteins, from which the transmembrane helix was absent. A significant number of computational studies have been performed on P450s, but the majority of these have been performed on the solubilised forms of P450s. Here we present a multiscale approach for modelling P450s, spanning from coarse-grained and atomistic molecular dynamics simulations to reaction modelling using hybrid quantum mechanics/molecular mechanics (QM/MM) methods. To our knowledge, this is the first application of such an integrated multiscale approach to modelling of a membrane-bound enzyme. We have applied this protocol to a key human P450 involved in drug metabolism: CYP3A4. A biologically realistic model of CYP3A4, complete with its transmembrane helix and a membrane, has been constructed and characterised. The dynamics of this complex have been studied, and the oxidation of the anticoagulant R-warfarin has been modelled in the active site. Calculations have also been performed on the soluble form of the enzyme in aqueous solution. Important differences are observed between the membrane and solution systems, most notably for the gating residues and channels that control access to the active site. The protocol that we describe here is applicable to other membrane-bound enzymes.

Use of density functional theory in drug metabolism studies

INTRODUCTION

The cytochrome P450 enzymes (CYPs) metabolize many drug compounds. They catalyze a wide variety of reactions, and potentially, a large number of different metabolites can be generated. Density functional theory (DFT) has, over the past decade, been shown to be a powerful tool to rationalize and predict the possible metabolites generated by the CYPs as well as other drug-metabolizing enzymes.

AREAS COVERED

We review applications of DFT on reactions performed by the CYPs and other drug-metabolizing enzymes able to perform oxidation reactions, with an emphasis on predicting which metabolites are produced. We also cover calculations of binding energies for complexes in which the ligands interact directly with the heme iron atom.

EXPERT OPINION

DFT is a useful tool for prediction of the site of metabolism. The use of small models of the enzymes work surprisingly well for most CYP isoforms. This is probably due to the fact that the binding of the substrates is not the major determinant. When binding of the substrate plays a significant role, the well-known issue of determining the free energy of binding is the challenge. How approaches taking the protein environment into account, like docking, MD and QM/MM, can be used are discussed.

Quantum mechanics/molecular mechanics modeling of regioselectivity of drug metabolism in cytochrome P450 2C9

Cytochrome P450 enzymes (P450s) are important in drug metabolism and have been linked to adverse drug reactions. P450s display broad substrate reactivity, and prediction of metabolites is complex. QM/MM studies of P450 reactivity have provided insight into important details of the reaction mechanisms and have the potential to make predictions of metabolite formation. Here we present a comprehensive study of the oxidation of three widely used pharmaceutical compounds (S-ibuprofen, diclofenac, and S-warfarin) by one of the major drug-metabolizing P450 isoforms, CYP2C9. The reaction barriers to substrate oxidation by the iron-oxo species (Compound I) have been calculated at the B3LYP-D/CHARMM27 level for different possible metabolism sites for each drug, on multiple pathways. In the cases of ibuprofen and warfarin, the process with the lowest activation energy is consistent with the experimentally preferred metabolite. For diclofenac, the pathway leading to the experimentally observed metabolite is not the one with the lowest activation energy. This apparent inconsistency with experiment might be explained by the two very different binding modes involved in oxidation at the two competing positions. The carboxylate of diclofenac interacts strongly with the CYP2C9 Arg108 side chain in the transition state for formation of the observed metabolite—but not in that for the competing pathway. We compare reaction barriers calculated both in the presence and in the absence of the protein and observe a marked improvement in selectivity prediction ability upon inclusion of the protein for all of the substrates studied. The barriers calculated with the protein are generally higher than those calculated in the gas phase. This suggests that active-site residues surrounding the substrate play an important role in controlling selectivity in CYP2C9. The results show that inclusion of sampling (particularly) and dispersion effects is important in making accurate predictions of drug metabolism selectivity of P450s using QM/MM methods.

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.

Testing high-level QM/MM methods for modeling enzyme reactions: acetyl-CoA deprotonation in citrate synthase

Combined quantum mechanics/molecular mechanics (QM/MM) calculations with high levels of correlated ab initio theory can now provide benchmarks for enzyme-catalyzed reactions. Here, we use such methods to test various QM/MM methods and the sensitivity of the results to details of the models for an important enzyme reaction, proton abstraction from acetyl-coenzyme A in citrate synthase. We calculate multiple QM/MM potential energy surfaces up to the local coupled cluster theory (LCCSD(T0)) level, with structures optimized at hybrid density functional theory and Hartree-Fock levels. The influence of QM methods, basis sets, and QM region size is shown to be significant. Correlated ab initio QM/MM calculations give barriers in agreement with experiment for formation of the acetyl-CoA enolate intermediate. In contrast, B3LYP fails to identify the enolate as an intermediate, whereas BH&HLYP does. The results indicate that QM/MM methods and setup should be tested, ideally using high-level calculations, to draw reliable mechanistic conclusions.

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.

A rational basis for the axial ligand effect in C-H oxidation by [MnO(porphyrin)(X)]+ (X = H2O, OH-, O2-) from a DFT study

Oxyl radical character in the MnO group of the title system is shown from a density functional theory study to be essential for efficient C-H cleavage, which is a key step in C-H oxidation. Since oxyl species have elongated Mn-O bonds relative to the more usual oxo species of type MnO, the normal expectation would be that high trans-influence ligands X should facilitate oxyl character by elongating the Mn-O bond and thus enhance both oxyl character and reactivity. Contrary to this expectation, but in line with the experimental data (Jin, N.; Ibrahim, M.; Spiro, T. G.; Groves, J. T. J. Am. Chem. Soc. 2007, 129, 12416), we find that reactivity increases along the series X = O(2-) < OH(-) < H2O for the following reasons. The ground-state singlet (S) is unreactive for all X, and only the higher-energy triplet (T) and quintet (Q) states have the oxyl character needed for reactivity, but the higher trans-influence X ligands are also shown to increase the S/T and S/Q gaps, thus making attainment of the needed T and Q states harder. The latter effect is dominant, and high trans-influence X ligands thus disfavor reaction. The higher reactivity in the presence of acid noted by Groves and co-workers is thus rationalized by the preference for having X = H2O over OH(-) or O(2-).

Structures of prostacyclin synthase and its complexes with substrate analog and inhibitor reveal a ligand-specific heme conformation change

Prostacyclin synthase (PGIS) is a cytochrome P450 (P450) enzyme that catalyzes production of prostacyclin from prostaglandin H(2). PGIS is unusual in that it catalyzes an isomerization rather than a monooxygenation, which is typical of P450 enzymes. To understand the structural basis for prostacyclin biosynthesis in greater detail, we have determined the crystal structures of ligand-free, inhibitor (minoxidil)-bound and substrate analog U51605-bound PGIS. These structures demonstrate a stereo-specific substrate binding and suggest features of the enzyme that facilitate isomerization. Unlike most microsomal P450s, where large substrate-induced conformational changes take place at the distal side of the heme, conformational changes in PGIS are observed at the proximal side and in the heme itself. The conserved and extensive heme propionate-protein interactions seen in all other P450s, which are largely absent in the ligand-free PGIS, are recovered upon U51605 binding accompanied by water exclusion from the active site. In contrast, when minoxidil binds, the propionate-protein interactions are not recovered and water molecules are largely retained. These findings suggest that PGIS represents a divergent evolution of the P450 family, in which a heme barrier has evolved to ensure strict binding specificity for prostaglandin H(2), leading to a radical-mediated isomerization with high product fidelity. The U51605-bound structure also provides a view of the substrate entrance and product exit channels.

Computational enzymology: modelling the mechanisms of biological catalysts

Simulations and modelling [e.g. with combined QM/MM (quantum mechanics/molecular mechanics) methods] are increasingly important in investigations of enzyme-catalysed reaction mechanisms. Calculations offer the potential of uniquely detailed, atomic-level insight into the fundamental processes of biological catalysis. Highly accurate methods promise quantitative comparison with experiments, and reliable predictions of mechanisms, revolutionizing enzymology.

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.

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.