Quantum Mechanical/Molecular Mechanical Investigation of the Mechanism of C−H Hydroxylation of Camphor by Cytochrome P450cam: Theory Supports a Two-State Rebound Mechanism

The rate-limiting step in P450 hydroxylation of hydrocarbons a direct comparison of the "somersault" versus the "consensus" mechanism involving compound I

Model theoretical quantum mechanical (QM) calculations are described for the P-450 hydroxylation of methane, isobutane, and camphor that compare the concerted somersault H-abstraction mechanism with the oxidation step involving Cpd I. Special emphasis has been placed on maintaining a balanced basis set in the oxidation step. QM calculations, employing the 6-311+G(d,p) basis set on the Fe atom and all of the key surrounding atoms involved in the C-H abstraction step, reaffirm a mechanism involving rearrangement of the iron hydroperoxide group (FeO-OH --> FeO...HO(*)) in concert with hydrogen abstraction from the C-H bond of the substrate by the incipient bound hydroxyl radical HO(*). The barrier for the somersault rearrangement of model Cpd 0 (FeO-OH) is calculated to be 21.4 kcal/mol in the absence of substrate. The overall activation energy for the oxidation of camphor involving the somersault motion of the FeO-OH group of P450 model porphyrin iron(III) hydroperoxide [Por(SH)Fe(III)-OOH(-)] --> [Por(SH)Fe(III)-O....HO(-)] in concert with hydrogen abstraction is DeltaE(++) = 12.4 kcal/mol. The corresponding abstraction of the hydrogen atom from the C-H bond of camphor by Cpd I has an activation barrier of 17.6 kcal/mol. Arguments are presented that the somersault rearrangement is induced by steric compression at the active site. Kinetic isotope effect data are discussed that provides compelling evidence for a rate-limiting step involving C-H bond cleavage.

What factors influence the rate constant of substrate epoxidation by compound I of cytochrome P450 and analogous iron(IV)-oxo oxidants?

The cytochromes P450 are a versatile range of mono-oxygenase enzymes that catalyze a variety of different chemical reactions, of which the key reactions include aliphatic hydroxylation and C=C double bond epoxidation. To establish the fundamental factors that govern substrate epoxidation by these enzymes we have done a systematic density functional theory study on substrate epoxidation by the active species of P450 enzymes, namely the iron(IV)-oxo porphyrin cation radical oxidant or Compound I. We show here, for the first time, that the rate constant of substrate epoxidation, and hence the activation energy, correlates with the ionization potential of the substrate as well as with intrinsic electronic properties of the active oxidant such as the polarizability volume. To explain these findings we present an electron-transfer model for the reaction mechanism that explains the factors that determine the barrier heights and developed a valence bond (VB) curve crossing mechanism to rationalize the observed trends. In addition, we have found a correlation for substrate epoxidation reactions catalyzed by a range of heme and nonheme iron(IV)-oxo oxidants with the strength of the O-H bond in the iron-hydroxo complex, i.e. BDE(OH), which is supported by the VB model. Finally, the fundamental factors that determine the regioselectivity change between substrate hydroxylation and epoxidation are discussed. It is shown that the regioselectivity of aliphatic hydroxylation versus double bond epoxidation is not influenced by the choice of the oxidant but is purely substrate dependent.

Manganese catalysts for C-H activation: an experimental/theoretical study identifies the stereoelectronic factor that controls the switch between hydroxylation and desaturation pathways

We describe competitive C-H bond activation chemistry of two types, desaturation and hydroxylation, using synthetic manganese catalysts with several substrates. 9,10-Dihydrophenanthrene (DHP) gives the highest desaturation activity, the final products being phenanthrene (P1) and phenanthrene 9,10-oxide (P3), the latter being thought to arise from epoxidation of some of the phenanthrene. The hydroxylase pathway also occurs as suggested by the presence of the dione product, phenanthrene-9,10-dione (P2), thought to arise from further oxidation of hydroxylation intermediate 9-hydroxy-9,10-dihydrophenanthrene. The experimental work together with the density functional theory (DFT) calculations shows that the postulated Mn oxo active species, [Mn(O)(tpp)(Cl)] (tpp = tetraphenylporphyrin), can promote the oxidation of dihydrophenanthrene by either desaturation or hydroxylation pathways. The calculations show that these two competing reactions have a common initial step, radical H abstraction from one of the DHP sp(3) C-H bonds. The resulting Mn hydroxo intermediate is capable of promoting not only OH rebound (hydroxylation) but also a second H abstraction adjacent to the first (desaturation). Like the active Mn(V)=O species, this Mn(IV)-OH species also has radical character on oxygen and can thus give H abstraction. Both steps have very low and therefore very similar energy barriers, leading to a product mixture. Since the radical character of the catalyst is located on the oxygen p orbital perpendicular to the Mn(IV)-OH plane, the orientation of the organic radical with respect to this plane determines which reaction, desaturation or hydroxylation, will occur. Stereoelectronic factors such as the rotational orientation of the OH group in the enzyme active site are thus likely to constitute the switch between hydroxylase and desaturase behavior.

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.

Quantum tunneling in testosterone 6beta-hydroxylation by cytochrome P450: reaction dynamics calculations employing multiconfiguration molecular-mechanical potential energy surfaces

Testosterone hydroxylation is a prototypical reaction of human cytochrome P450 3A4, which metabolizes about 50% of oral drugs on the market. Reaction dynamics calculations were carried out for the testosterone 6beta-hydrogen abstraction and the 6beta-d(1)-testosterone 6beta-duterium abstraction employing a model that consists of the substrate and the active oxidant compound I. The calculations were performed at the level of canonical variational transition state theory with multidimensional tunneling and were based on a semiglobal full-dimensional potential energy surface generated by the multiconfiguration molecular mechanics technique. The tunneling coefficients were found to be around 3, indicating substantial contributions by quantum tunneling. However, the tunneling made only modest contributions to the kinetic isotope effects. The kinetic isotope effects were computed to be about 2 in the doublet spin state and about 5 in the quartet spin state.

Insights into the (superoxo)Fe(III)Fe(III) intermediate and reaction mechanism of myo-inositol oxygenase: DFT and ONIOM(DFT:MM) study

The (superoxo)Fe(III)Fe(III) reactive species and the catalytic reaction mechanism of a diiron enzyme, myo-inositol oxygenase (MIOX), were theoretically investigated by means of density functional theory (DFT) and ONIOM quantum mechanical/molecular mechanical (QM/MM) approaches. The ground state of the (superoxo)Fe(III)Fe(III) intermediate was shown to have a side-on coordination geometry and an S = 1/2 spin state, wherein the two iron sites are antiferromagnetically coupled while the superoxide site and the nearest iron are ferromagnetically coupled. A full reaction pathway leading to a D-glucuronate product from myo-inositol was proposed based on ONIOM computational results. Two major roles of the enzyme surrounding during the catalytic reaction were identified. One is to facilitate the initial H-abstraction step, and the other is to restrict the movement of the substrate via H-bonding interactions in order to avoid unwanted side reactions. In our proposed mechanism, O-O bond cleavage has the highest barrier, thus constituting the rate-limiting step of the reaction. The unique role of the bridging hydroxide ligand as a catalytic base was also identified.

Finite-temperature effects in enzymatic reactions — Insights from QM/MM free-energy simulations

We report potential-energy and free-energy data for three enzymatic reactions: carbon–halogen bond formation in fluorinase, hydrogen abstraction from camphor in cytochrome P450cam, and chorismate-to-prephenate Claisen rearrangement in chorismate mutase. The results were obtained by combined quantum mechanics/molecular mechanics (QM/MM) optimizations and two types of QM/MM free-energy simulations (free-energy perturbation and umbrella sampling) using semi-empirical or density-functional QM methods. Based on these results and our previously published free-energy data on electrophilic substitution in para-hydroxybenzoate hydroxylase, we discuss the importance of finite-temperature effects in the chemical step of enzyme reactions. We find that the entropic contribution to the activation barrier is generally rather small, usually of the order of 5 kJ mol–1 or less, consistent with the notion that enzymes bind and pre-organize the reactants in the active site. A somewhat larger entropic contribution is encountered in the case of chorismate mutase where the pericyclic transition state is intrinsically more rigid than the chorismate reactant (also in the enzyme). The present results suggest that barriers from QM/MM geometry optimization may often be close to free-energy barriers for the chemical step in enzymatic reactions.

QM/MM Simulation on P450 BM3 Enzyme Catalysis Mechanism

Using a structure generated by induced fit modeling of the protein-ligand complex, the reaction path for hydrogen atom abstraction in P450 BM3 is studied by means of mixed QM/MM methods to determine the structures and energetics along the reaction path. The IFD structure is suitable for hydrogen atom abstraction at the ω-1 position. The electronic structures obtained are similar to those observed in P450 cam. We show that the barrier for the hydrogen abstraction step from QM/MM modeling is 13.3 kcal/mol in quartet and 15.6 kcal/mol in doublet. Although there is some strain energy present in the ligand, the activation barrier is not dramatically affected. A crystal water molecule, HOH502, plays a role as catalyst and decreases the activation barrier by about 2 kcal/mol and reaction energy by about 3-4 kcal/mol. In order to achieve reactive chemistry at the remaining experimentally observed positions in the hydrocarbon tail of the ligand, other structures would have to be utilized as a starting point for the reaction. Finally, the present results still leave open the question of whether DFT methods provide an accurate computation of the barrier height in the P450 hydrogen atom abstraction reaction.

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.

Progress in ab initio QM/MM free-energy simulations of electrostatic energies in proteins: accelerated QM/MM studies of pKa, redox reactions and solvation free energies

Hybrid quantum mechanical/molecular mechanical (QM/MM) approaches have been used to provide a general scheme for chemical reactions in proteins. However, such approaches still present a major challenge to computational chemists, not only because of the need for very large computer time in order to evaluate the QM energy but also because of the need for proper computational sampling. This review focuses on the sampling issue in QM/MM evaluations of electrostatic energies in proteins. We chose this example since electrostatic energies play a major role in controlling the function of proteins and are key to the structure-function correlation of biological molecules. Thus, the correct treatment of electrostatics is essential for the accurate simulation of biological systems. Although we will be presenting different types of QM/MM calculations of electrostatic energies (and related properties) here, our focus will be on pKa calculations. This reflects the fact that pKa's of ionizable groups in proteins provide one of the most direct benchmarks for the accuracy of electrostatic models of macromolecules. While pKa calculations by semimacroscopic models have given reasonable results in many cases, existing attempts to perform pKa calculations using QM/MM-FEP have led to discrepancies between calculated and experimental values. In this work, we accelerate our QM/MM calculations using an updated mean charge distribution and a classical reference potential. We examine both a surface residue (Asp3) of the bovine pancreatic trypsin inhibitor and a residue buried in a hydrophobic pocket (Lys102) of the T4-lysozyme mutant. We demonstrate that, by using this approach, we are able to reproduce the relevant side chain pKa's with an accuracy of 3 kcal/mol. This is well within the 7 kcal/mol energy difference observed in studies of enzymatic catalysis, and is thus sufficient accuracy to determine the main contributions to the catalytic energies of enzymes. We also provide an overall perspective of the potential of QM/MM calculations in general evaluations of electrostatic free energies, pointing out that our approach should provide a very powerful and accurate tool to predict the electrostatics of not only solution but also enzymatic reactions, as well as the solvation free energies of even larger systems, such as nucleic acid bases incorporated into DNA.

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.

Density functional studies on thromboxane biosynthesis: mechanism and role of the heme-thiolate system

Reaction mechanisms for the isomerization of prostaglandin H(2) to thromboxane A(2), and degradation to 12-L-hydroxy-5,8,10-heptadecatrienoic acid (HHT) and malondialdehyde (MDA), catalyzed by thromboxane synthase, were investigated using the unrestricted Becke-three-parameter plus Lee-Yang-Parr (UB3LYP) density functional level theory. In addition to the reaction pathway through Fe(IV)-porphyrin intermediates, a new reaction pathway through Fe(III)-porphyrin pi-cation radical intermediates was found. Both reactions proceed with the homolytic cleavage of endoperoxide O-O to give an alkoxy radical. This intermediate converts into an allyl radical intermediate by a C-C homolytic cleavage, followed by the formation of thromboxane A(2) having a 6-membered ring through a one electron transfer, or the degradation into HHT and MDA. The proposed mechanism shows that an iron(III)-containing system having electron acceptor ability is essential for the 6-membered ring formation leading to thromboxane A(2). Our results suggest that the step of the endoperoxide O-O homolytic bond cleavage has the highest activation energy following the binding of prostaglandin H(2) to thromboxane synthase.

Human cytochromes P450 in the metabolism of drugs: new molecular models of enzyme-substrate interactions

The overall predictive ability of molecular modelling, as applied to the cytochrome P450 (CYP) system, is analysed in the light of current developments in a variety of techniques, including X-ray crystallography, molecular biology, enzyme kinetics, molecular mechanics and dynamics, in relation to its relevance to drug metabolism in humans. This review demonstrates that it is possible to generate realistic models for the major human CYPs, which metabolise xenobiotics that compare favourably with crystal structures, and thus may be used to derive substrate binding energies that agree closely with experimental K(m) values obtained from enzyme kinetics.

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-).