QM/MM study of the taxadiene synthase mechanism

Combined quantum mechanics/molecular mechanics (QM/MM) calculations were used to investigate the reaction mechanism of taxadiene synthase (TXS). TXS catalyzes the cyclization of geranylgeranyl diphosphate (GGPP) to taxadiene (T) and four minor cyclic products. All these products originate from the deprotonation of carbocation intermediates. The reaction profiles for the conversion of GGPP to T as well as to minor products were calculated for different configurations of relevant TXS carbocation complexes. The QM region was treated at the M06‐2X/TZVP level, while the CHARMM27 force field was used to describe the MM region. The QM/MM calculations suggest a reaction pathway for the conversion of GGPP to T, which slightly differs from previous proposals regarding the number of reaction steps and the conformation of the carbocations. The QM/MM results also indicate that the formation of minor products via water‐assisted deprotonation of the carbocations is highly exothermic, by about −7 to −23 kcal/mol. Curiously, however, the computed barriers and reaction energies indicate that the formation of some of the minor products is more facile than the formation of T. Thus, the present QM/MM calculations provide detailed insights into possible reaction pathways and into the origin of the promiscuity of TXS, but they do not reproduce the product distribution observed experimentally. © 2019 Wiley Periodicals, Inc.

[1]  S. C. Rogers,et al.  QUASI: A general purpose implementation of the QM/MM approach and its application to problems in catalysis , 2003 .

[2]  Meimei Xu,et al.  Following evolution's lead to a single residue switch for diterpene synthase product outcome , 2007, Proceedings of the National Academy of Sciences.

[3]  D. Major,et al.  Catalytic Control in the Facile Proton Transfer in Taxadiene Synthase , 2017 .

[4]  Alexander D. MacKerell,et al.  Extending the treatment of backbone energetics in protein force fields: Limitations of gas‐phase quantum mechanics in reproducing protein conformational distributions in molecular dynamics simulations , 2004, J. Comput. Chem..

[5]  Frank Neese,et al.  The ORCA program system , 2012 .

[6]  Jorge Nocedal,et al.  On the limited memory BFGS method for large scale optimization , 1989, Math. Program..

[7]  Stefan Grimme,et al.  Semiempirical GGA‐type density functional constructed with a long‐range dispersion correction , 2006, J. Comput. Chem..

[8]  W. Thiel,et al.  Amine oxidation mediated by lysine-specific demethylase 1: quantum mechanics/molecular mechanics insights into mechanism and role of lysine 661. , 2013, Journal of the American Chemical Society.

[9]  R. Croteau,et al.  Stereochemistry of the macrocyclization and elimination steps in taxadiene biosynthesis through deuterium labeling. , 2005, Journal of Organic Chemistry.

[10]  Adrian J Mulholland,et al.  Conformational change and ligand binding in the aristolochene synthase catalytic cycle. , 2013, Biochemistry.

[11]  J. Nocedal.  Updating Quasi-Newton Matrices With Limited Storage , 1980 .

[12]  C. Brooks,et al.  Multiple intermediates, diverse conformations, and cooperative conformational changes underlie the catalytic hydride transfer reaction of dihydrofolate reductase. , 2013, Topics in current chemistry.

[13]  D. Tantillo,et al.  Quantum chemical dissection of the classic terpinyl/pinyl/bornyl/camphyl cation conundrum-the role of pyrophosphate in manipulating pathways to monoterpenes. , 2010, Organic & biomolecular chemistry.

[14]  J. Pople,et al.  Self—Consistent Molecular Orbital Methods. XII. Further Extensions of Gaussian—Type Basis Sets for Use in Molecular Orbital Studies of Organic Molecules , 1972 .

[15]  M. Mirjalili,et al.  Production of the anticancer drug taxol in Taxus baccata suspension cultures: A review , 2011 .

[16]  S. Benkovic,et al.  Flexibility, diversity, and cooperativity: pillars of enzyme catalysis. , 2011, Biochemistry.

[17]  Walter Thiel,et al.  Quantum mechanical/molecular mechanical study on the mechanism of the enzymatic Baeyer-Villiger reaction. , 2012, Journal of the American Chemical Society.

[18]  Dean J Tantillo,et al.  The taxadiene-forming carbocation cascade. , 2011, Journal of the American Chemical Society.

[19]  Danna Zhou,et al.  d. , 1934, Microbial pathogenesis.

[20]  W. Eisenreich,et al.  Identification of amino acid networks governing catalysis in the closed complex of class I terpene synthases , 2016, Proceedings of the National Academy of Sciences.

[21]  M. Head‐Gordon,et al.  Long-range corrected hybrid density functionals with damped atom-atom dispersion corrections. , 2008, Physical chemistry chemical physics : PCCP.

[22]  Walter Thiel,et al.  Quantum Mechanics/Molecular Mechanics Insights into the Enantioselectivity of the O-Acetylation of (R,S)-Propranolol Catalyzed by Candida antarctica Lipase B , 2017 .

[23]  Walter Thiel,et al.  A combined QM/MM study on the reductive half-reaction of xanthine oxidase: substrate orientation and mechanism. , 2009, Journal of the American Chemical Society.

[24]  G. Loake,et al.  Paclitaxel: biosynthesis, production and future prospects. , 2014, New biotechnology.

[25]  D. Tantillo.  Biosynthesis via carbocations: theoretical studies on terpene formation. , 2011, Natural product reports.

[26]  D. Cane,et al.  Exploring the Influence of Domain Architecture on the Catalytic Function of Diterpene Synthases. , 2017, Biochemistry.

[27]  J. Baker.  An algorithm for the location of transition states , 1986 .

[28]  D. Major,et al.  Electrostatically guided dynamics--the root of fidelity in a promiscuous terpene synthase? , 2012, Journal of the American Chemical Society.

[29]  J. Noel,et al.  Pre-steady-state study of recombinant sesquiterpene cyclases. , 1997, Biochemistry.

[30]  Electrostatic Control of Chemistry in Terpene Cyclases , 2017 .

[31]  Ajit Banerjee,et al.  Search for stationary points on surfaces , 1985 .

[32]  H. Williams,et al.  Studies on taxadiene synthase: interception of the cyclization cascade at the isocembrene stage with GGPP analogues. , 2005, The Journal of organic chemistry.

[33]  Jiali Gao,et al.  Chemical Control in the Battle against Fidelity in Promiscuous Natural Product Biosynthesis: The Case of Trichodiene Synthase. , 2017, ACS catalysis.

[34]  David W. Christianson,et al.  Structural and Chemical Biology of Terpenoid Cyclases , 2017, Chemical reviews.

[35]  D. Truhlar,et al.  The M06 suite of density functionals for main group thermochemistry, thermochemical kinetics, noncovalent interactions, excited states, and transition elements: two new functionals and systematic testing of four M06-class functionals and 12 other functionals , 2008 .

[36]  R. Croteau,et al.  Taxadiene synthase-catalyzed cyclization of 6-fluorogeranylgeranyl diphosphate to 7-fluoroverticillenes. , 2005, Journal of the American Chemical Society.

[37]  D. Cane,et al.  Pre-steady-state kinetic analysis of the trichodiene synthase reaction pathway. , 1997, Biochemistry.

[38]  Tsuyoshi Murata,et al.  {m , 1934, ACML.

[39]  David W. Christianson,et al.  Taxadiene Synthase Structure and Evolution of Modular Architecture in Terpene Biosynthesis , 2010, Nature.

[40]  R. Peters,et al.  Electrostatic effects on (di)terpene synthase product outcome. , 2011, Chemical communications.

[41]  D. Whittington,et al.  Bornyl diphosphate synthase: Structure and strategy for carbocation manipulation by a terpenoid cyclase , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[42]  R. Croteau,et al.  Intramolecular proton transfer in the cyclization of geranylgeranyl diphosphate to the taxadiene precursor of taxol catalyzed by recombinant taxadiene synthase. , 2000, Chemistry & biology.

[43]  Luigi Di Costanzo,et al.  X-ray crystal structure of aristolochene synthase from Aspergillus terreus and evolution of templates for the cyclization of farnesyl diphosphate. , 2007, Biochemistry.

[44]  Walter Thiel,et al.  Molecular dynamics study of taxadiene synthase catalysis , 2018, J. Comput. Chem..

[45]  W Smith,et al.  DL_POLY_2.0: a general-purpose parallel molecular dynamics simulation package. , 1996, Journal of molecular graphics.

[46]  R. Croteau,et al.  Heterologous expression and characterization of a "Pseudomature" form of taxadiene synthase involved in paclitaxel (Taxol) biosynthesis and evaluation of a potential intermediate and inhibitors of the multistep diterpene cyclization reaction. , 2000, Archives of biochemistry and biophysics.

[47]  Dan Thomas Major,et al.  Catalytic control in terpenoid cyclases: multiscale modeling of thermodynamic, kinetic, and dynamic effects. , 2014, Current opinion in chemical biology.

[48]  R. Croteau,et al.  Structure of limonene synthase, a simple model for terpenoid cyclase catalysis , 2007, Proceedings of the National Academy of Sciences.

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

[50]  Alexander D. MacKerell,et al.  All-atom empirical potential for molecular modeling and dynamics studies of proteins. , 1998, The journal of physical chemistry. B.

[51]  D. Major,et al.  Slow-Starter Enzymes: Role of Active-Site Architecture in the Catalytic Control of the Biosynthesis of Taxadiene by Taxadiene Synthase. , 2018, Biochemistry.

[52]  F. Weigend,et al.  Balanced basis sets of split valence, triple zeta valence and quadruple zeta valence quality for H to Rn: Design and assessment of accuracy. , 2005, Physical chemistry chemical physics : PCCP.

[53]  R. Patel,et al.  Tour de paclitaxel: biocatalysis for semisynthesis. , 1998, Annual review of microbiology.

[54]  J. Noel,et al.  Structural basis for cyclic terpene biosynthesis by tobacco 5-epi-aristolochene synthase. , 1997, Science.

[55]  A. Schäfer,et al.  Fully optimized contracted Gaussian basis sets of triple zeta valence quality for atoms Li to Kr , 1994 .

[56]  Jianpeng Ma,et al.  CHARMM: The biomolecular simulation program , 2009, J. Comput. Chem..

[57]  Susan C. Roberts,et al.  Production and engineering of terpenoids in plant cell culture. , 2007, Nature chemical biology.

[58]  D. Christianson.  Unearthing the roots of the terpenome. , 2008, Current opinion in chemical biology.

[59]  Gordon M. Cragg.  Paclitaxel (Taxol®): A success story with valuable lessons for natural product drug discovery and development , 1998, Medicinal research reviews.

[60]  Frank Neese,et al.  Sparse maps--A systematic infrastructure for reduced-scaling electronic structure methods. II. Linear scaling domain based pair natural orbital coupled cluster theory. , 2016, The Journal of chemical physics.

[61]  Jiali Gao,et al.  Comment on "Substrate Folding Modes in Trichodiene Synthase: A Determinant of Chemo- and Stereoselectivity". , 2018, ACS catalysis.

[62]  W. Thiel,et al.  Hybrid Models for Combined Quantum Mechanical and Molecular Mechanical Approaches , 1996 .

[63]  N. Lewis,et al.  Cyclization of Geranylgeranyl Diphosphate to Taxa-4 (5), 11 (12) -diene Is the Committed Step of Taxol Biosynthesis in Pacific Yew (*) , 1995, The Journal of Biological Chemistry.