QM/MM modeling of the hydroxylation of the androstenedione substrate catalyzed by cytochrome P450 aromatase (CYP19A1)

CYP19A1 aromatase is a member of the Cytochrome P450 family of hemeproteins, and is the enzyme responsible for the final step of the androgens conversion into the corresponding estrogens, via a three‐step oxidative process. For this reason, the inhibition of this enzyme plays an important role in the treatment of hormone‐dependent breast cancer. The first catalytic subcycle, corresponding to the hydroxilation of androstenedione, has been proposed to occur through a first hydrogen abstraction and a subsequent oxygen rebound step. In present work, we have studied the mechanism of the first catalytic subcycle by means of hybrid quantum mechanics/molecular mechanics methods. The inclusion of the protein flexibility has been achieved by means of Free Energy Perturbation techniques, giving rise to a free energy of activation for the hydrogen abstraction step of 13.5 kcal/mol. The subsequent oxygen rebound step, characterized by a small free energy barrier (1.5 kcal/mol), leads to the hydroxylated products through a highly exergonic reaction. In addition, an analysis of the primary deuterium kinetic isotopic effects, calculated for the hydrogen abstraction step, reveals values (∼10) overpassing the semiclassical limit for the CH, indicating the presence of a substantial tunnel effect. Finally, a decomposition analysis of the interaction energy for the substrate and cofactor in the active site is also discussed. According to our results, the role of the enzymatic environment consists of a transition state stabilization by means of dispersive and polarization effects. © 2015 Wiley Periodicals, Inc.

[1]  W. L. Jorgensen,et al.  Development and Testing of the OPLS All-Atom Force Field on Conformational Energetics and Properties of Organic Liquids , 1996 .

[2]  Shaik,et al.  On the "Rebound" Mechanism of Alkane Hydroxylation by Cytochrome P450: Electronic Structure of the Intermediate and the Electron Transfer Character in the Rebound Step. , 1999, Angewandte Chemie.

[3]  Michael T. Green,et al.  X-ray absorption spectroscopy of chloroperoxidase compound I: Insight into the reactive intermediate of P450 chemistry. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[4]  R. Obach,et al.  Drug Metabolizing Enzymes : Cytochrome P450 and Other Enzymes in Drug Discovery and Development , 2003 .

[5]  Jan H. Jensen,et al.  Very fast empirical prediction and rationalization of protein pKa values , 2005, Proteins.

[6]  S. Sligar,et al.  Kinetic solvent isotope effect in human P450 CYP17A1-mediated androgen formation: evidence for a reactive peroxoanion intermediate. , 2013, Journal of the American Chemical Society.

[7]  J. Shimada,et al.  Unique structural and electronic features of perferryl-oxo oxidant in Cytochrome P450. , 2011, The journal of physical chemistry. B.

[8]  Michael T. Green,et al.  Cytochrome P450 Compound I: Capture, Characterization, and C-H Bond Activation Kinetics , 2010, Science.

[9]  I. Williams.  Force-constant computations in cartesian coordinates. Elimination of translational and rotational contributions , 1983 .

[10]  C. Che,et al.  Bis(sulfonylimide)ruthenium(VI) porphyrins: X-ray crystal structure and mechanism of C-H bond amination by density functional theory calculations. , 2013, Chemistry.

[11]  J Berendzen,et al.  The catalytic pathway of cytochrome p450cam at atomic resolution. , 2000, Science.

[12]  M. Calder,et al.  Mechanistic studies on C-19 demethylation in oestrogen biosynthesis. , 1982, The Biochemical journal.

[13]  K. Sen,et al.  Coupled electron transfer and proton hopping in the final step of CYP19-catalyzed androgen aromatization. , 2012, Biochemistry.

[14]  Multiple Low-Lying States for Compound I of P450cam and Chloroperoxidase Revealed from Multireference Ab Initio QM/MM Calculations. , 2010, Journal of chemical theory and computation.

[15]  W. L. Jorgensen,et al.  Comparison of simple potential functions for simulating liquid water , 1983 .

[16]  S. Sligar,et al.  Kinetic solvent isotope effect in steady‐state turnover by CYP19A1 suggests involvement of Compound 1 for both hydroxylation and aromatization steps , 2014, FEBS letters.

[17]  S. Skinner,et al.  The intermediary role of a 19-oxoandrogen in the biosynthesis of oestrogen. , 1968, The Biochemical journal.

[18]  John T. Groves,et al.  Aliphatic hydroxylation via oxygen rebound. Oxygen transfer catalyzed by iron , 1976 .

[19]  T. Arunachalam,et al.  Biosynthesis of estrogens: the steric mode of the initial C-19 hydroxylation of androgens by human placental aromatase , 1983 .

[20]  W. Pangborn,et al.  X-ray structure of human aromatase reveals an androgen-specific active site , 2010, The Journal of Steroid Biochemistry and Molecular Biology.

[21]  S. Sligar,et al.  The critical iron-oxygen intermediate in human aromatase. , 2009, Biochemical and biophysical research communications.

[22]  Ortiz de Montellano,et al.  Cytochrome P-450: Structure, Mechanism, and Biochemistry , 1986 .

[23]  Weitao Yang,et al.  Free energy calculation on enzyme reactions with an efficient iterative procedure to determine minimum energy paths on a combined ab initio QM/MM potential energy surface , 2000 .

[24]  Chun Wu,et al.  How accurate are the popular PCM/GB continuum solvation models for calculating the solvation energies of amino acid side-chain analogs? , 2014, Theoretical Chemistry Accounts.

[25]  S. Sligar,et al.  Thirty years of microbial P450 monooxygenase research: peroxo-heme intermediates--the central bus station in heme oxygenase catalysis. , 2005, Biochemical and biophysical research communications.

[26]  S. Shaik,et al.  Two-state reactivity mechanisms of hydroxylation and epoxidation by cytochrome P-450 revealed by theory. , 2002, Current opinion in chemical biology.

[27]  S. Sligar.  Glimpsing the Critical Intermediate in Cytochrome P450 Oxidations , 2010, Science.

[28]  S. Davis,et al.  Minireview: Aromatase and the Regulation of Estrogen Biosynthesis-Some New Perspectives. , 2001, Endocrinology.

[29]  Yong Wang,et al.  P450 enzymes: their structure, reactivity, and selectivity-modeled by QM/MM calculations. , 2010, Chemical reviews.

[30]  W. Boon,et al.  Aromatase--a brief overview. , 2002, Annual review of physiology.

[31]  D. Nelson.  The Cytochrome P450 Homepage , 2009, Human Genomics.

[32]  M. Newcomb,et al.  Spectra and kinetic studies of the compound I derivative of cytochrome P450 119. , 2008, Journal of the American Chemical Society.

[33]  S. Shaik,et al.  A two-state reactivity model explains unusual kinetic isotope effect patterns in C-H bond cleavage by nonheme oxoiron(IV) complexes. , 2009, Angewandte Chemie.

[34]  Ruth Nussinov,et al.  QM/MM study of the active species of the human cytochrome P450 3A4, and the influence thereof of the multiple substrate binding. , 2007, The journal of physical chemistry. B.

[35]  Jon Baker,et al.  The generation and use of delocalized internal coordinates in geometry optimization , 1996 .

[36]  S. Skinner,et al.  The stereospecific removal of a C-19 hydrogen atom in oestrogen biosynthesis. , 1969, The Biochemical journal.

[37]  D. Arigoni,et al.  Stereospecificity of oxidation at C-19 in oestrogen biosynthesis , 1975 .

[38]  P. Kollman,et al.  Atomic charges derived from semiempirical methods , 1990 .

[39]  Sergio Martí,et al.  Improving the QM/MM Description of Chemical Processes:  A Dual Level Strategy To Explore the Potential Energy Surface in Very Large Systems. , 2005, Journal of chemical theory and computation.

[40]  Constrained optimization in delocalized internal coordinates , 1997 .

[41]  T. Arunachalam,et al.  Biosynthesis of estrogens: aromatization of (19R)-, (19S)-, and (19RS)-[19-3H,2H,1H]-3.beta.-hydroxyandrost-5-en-17-ones by human placental aromatase , 1986 .

[42]  Sam P. de Visser,et al.  Propene activation by the oxo-iron active species of taurine/α- ketoglutarate dioxygenase (TauD) enzyme. How does the catalysis compare to heme-enzymes? , 2006 .

[43]  F. Neese,et al.  Multireference ab initio quantum mechanics/molecular mechanics study on intermediates in the catalytic cycle of cytochrome P450(cam). , 2008, The journal of physical chemistry. A.

[44]  J. Fishman,et al.  Mechanism of estrogen biosynthesis. Participation of multiple enzyme sites in placental aromatase hydroxylations. , 1981, The Journal of biological chemistry.

[45]  Marc Albe,et al.  The dynamo library for molecular simulations using hybrid quantum mechanical and molecular mechanical potentials , 2000 .

[46]  F. Guengerich,et al.  Common and uncommon cytochrome P450 reactions related to metabolism and chemical toxicity. , 2001, Chemical research in toxicology.

[47]  S. Shaik,et al.  Axial ligand tuning of a nonheme iron(IV)–oxo unit for hydrogen atom abstraction , 2007, Proceedings of the National Academy of Sciences.

[48]  Sason Shaik,et al.  Theoretical perspective on the structure and mechanism of cytochrome P450 enzymes. , 2005, Chemical reviews.

[49]  W. Pangborn,et al.  Structural basis for androgen specificity and oestrogen synthesis in human aromatase , 2009, Nature.

[50]  Jeffrey P. Jones,et al.  A New Mechanistic Probe for Cytochrome P450: An Application of Isotope Effect Profiles , 1997 .

[51]  Jorge Nocedal,et al.  A Limited Memory Algorithm for Bound Constrained Optimization , 1995, SIAM J. Sci. Comput..

[52]  W. Pryor.  Cytochrome P450: Structure, mechanism, and biochemistry , 1996 .

[53]  I. Schlichting,et al.  The status of high-valent metal oxo complexes in the P450 cytochromes. , 2006, Journal of inorganic biochemistry.

[54]  Qin Wang,et al.  Quantitative production of compound I from a cytochrome P450 enzyme at low temperatures. Kinetics, activation parameters, and kinetic isotope effects for oxidation of benzyl alcohol. , 2009, Journal of the American Chemical Society.

[55]  Laxmikant V. Kalé,et al.  Scalable molecular dynamics with NAMD , 2005, J. Comput. Chem..

[56]  M. J. Coon,et al.  Cytochrome P450-Catalyzed Hydroxylation of Mechanistic Probes that Distinguish between Radicals and Cations. Evidence for Cationic but Not for Radical Intermediates , 2000 .

[57]  W. R. Wadt,et al.  Ab initio effective core potentials for molecular calculations. Potentials for K to Au including the outermost core orbitals , 1985 .

[58]  S. Shaik,et al.  Alkane Hydroxylation by Cytochrome P450: Is Kinetic Isotope Effect a Reliable Probe of Transition State Structure? , 2000 .

[59]  P. Kollman,et al.  An approach to computing electrostatic charges for molecules , 1984 .

[60]  Direct Combined ab Initio/Transition State Theory Study of the Kinetics of the Abstraction Reactions of Halogenated Methanes with Hydrogen Atoms , 2004 .

[61]  Kenneth B. Wiberg,et al.  Application of the pople-santry-segal CNDO method to the cyclopropylcarbinyl and cyclobutyl cation and to bicyclobutane , 1968 .

[62]  Ilme Schlichting,et al.  Structure and chemistry of cytochrome P450. , 2005, Chemical reviews.

[63]  M. J. Coon,et al.  Two distinct electrophilic oxidants effects hydroxylation in cytochrome P-450-catalyzed reactions [20] , 1998 .

[64]  Sason Shaik,et al.  What factors affect the regioselectivity of oxidation by cytochrome p450? A DFT study of allylic hydroxylation and double bond epoxidation in a model reaction. , 2002, Journal of the American Chemical Society.

[65]  Hermann Stoll,et al.  Results obtained with the correlation energy density functionals of becke and Lee, Yang and Parr , 1989 .

[66]  R. Santen,et al.  History of aromatase: saga of an important biological mediator and therapeutic target. , 2009, Endocrine reviews.

[67]  J. Oláh,et al.  Enolization as an alternative proton delivery pathway in human aromatase (P450 19A1). , 2014, The journal of physical chemistry. B.

[68]  A. Becke.  Density-functional thermochemistry. III. The role of exact exchange , 1993 .

[69]  R. Brueggemeier,et al.  Aromatase inhibitors in the treatment of breast cancer. , 2005, Endocrine reviews.

[70]  F. Guengerich.  Cytochrome p450 and chemical toxicology. , 2008, Chemical research in toxicology.

[71]  P. R. Montellano.  Hydrocarbon hydroxylation by cytochrome P450 enzymes. , 2010 .

[72]  S. Sligar,et al.  Cytochrome p450 compound I. , 2006, Journal of the American Chemical Society.

[73]  M. Akhtar,et al.  A review of mechanistic studies on aromatase (CYP19) and 17α-hydroxylase-17,20-lyase (CYP17) , 2011, The Journal of Steroid Biochemistry and Molecular Biology.

[74]  M. J. Coon,et al.  Hydroxylation by the hydroperoxy-iron species in cytochrome P450 enzymes. , 2004, Journal of the American Chemical Society.

[75]  R. P. Bell,et al.  The tunnel effect in chemistry , 1959 .

[76]  M. J. Coon,et al.  Kinetic isotope effects implicate two electrophilic oxidants in cytochrome p450-catalyzed hydroxylations. , 2003, Journal of the American Chemical Society.