Revised mechanism of carboxylation of ribulose‐1,5‐biphosphate by rubisco from large scale quantum chemical calculations

Here, we describe a computational approach for studying enzymes that catalyze complex multi‐step reactions and apply it to Ribulose 1,5‐bisphosphate carboxylase–oxygenase (Rubisco), the enzyme that fixes atmospheric carbon dioxide within photosynthesis. In the 5‐step carboxylase reaction, the substrate Ribulose‐1,5‐bisphosphate (RuBP) first binds Rubisco and undergoes enolization before binding the second substrate, CO2. Hydration of the RuBP.CO2 complex is followed by CC bond scission and stereospecific protonation. However, details of the roles and protonation states of active‐site residues, and sources of protons and water, remain highly speculative. Large‐scale computations on active‐site models provide a means to better understand this complex chemical mechanism. The computational protocol comprises a combination of hybrid semi‐empirical quantum mechanics and molecular mechanics within constrained molecular dynamics simulations, together with constrained gradient minimization calculations using density functional theory. Alternative pathways for hydration of the RuBP.CO2 complex and associated active‐site protonation networks and proton and water sources were investigated. The main findings from analysis of the resulting energetics advocate major revision to existing mechanisms such that: hydration takes place anti to the CO2; both hydration and CC bond scission require early protonation of CO2 in the RuBP.CO2 complex; CC bond scission and stereospecific protonation reactions are concerted and, effectively, there is only one stable intermediate, the C3‐gemdiolate complex. Our main conclusions for interpreting enzyme kinetic results are that the gemdiolate may represent the elusive Michaelis–Menten‐like complex corresponding to the empirical Km (=Kc) with turnover to product via bond scission concerted with stereospecific protonation consistent with the observed catalytic rate. © 2018 Wiley Periodicals, Inc.

[1]  Andrews Tj,et al.  Quantum chemical analysis of the enolization of ribulose bisphosphate: the first hurdle in the fixation of CO2 by Rubisco. , 1998 .

[2]  G. Groenhof.  Introduction to QM/MM simulations. , 2013, Methods in molecular biology.

[3]  I. Andersson,et al.  Large structures at high resolution: the 1.6 A crystal structure of spinach ribulose-1,5-bisphosphate carboxylase/oxygenase complexed with 2-carboxyarabinitol bisphosphate. , 1996, Journal of molecular biology.

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

[5]  Alfredo Mayall Simas,et al.  RM1: A reparameterization of AM1 for H, C, N, O, P, S, F, Cl, Br, and I , 2006, J. Comput. Chem..

[6]  D. Hanson.  Breaking the rules of Rubisco catalysis , 2016, Journal of experimental botany.

[7]  S. Wildman.  Along the trail from Fraction I protein to Rubisco (ribulose bisphosphate carboxylase-oxygenase) , 2002 .

[8]  P. Krüger,et al.  Free energy as the potential of mean constraint force , 1996 .

[9]  I. Andersson,et al.  Structure of a product complex of spinach ribulose-1,5-bisphosphate carboxylase/oxygenase. , 1997, Biochemistry.

[10]  Luca Monticelli,et al.  Biomolecular Simulations: Methods and Protocols , 2019, Methods in Molecular Biology.

[11]  D. Jordan,et al.  Species variation in the specificity of ribulose biphosphate carboxylase/oxygenase , 1981, Nature.

[12]  Babu Kannappan,et al.  Redefinition of rubisco carboxylase reaction reveals origin of water for hydration and new roles for active-site residues. , 2008, Journal of the American Chemical Society.

[13]  P. J. Andralojc,et al.  Manipulation of Rubisco: the amount, activity, function and regulation. , 2003, Journal of experimental botany.

[14]  Michiel Sprik,et al.  Free energy from constrained molecular dynamics , 1998 .

[15]  G. Farquhar,et al.  Despite slow catalysis and confused substrate specificity, all ribulose bisphosphate carboxylases may be nearly perfectly optimized. , 2006, Proceedings of the National Academy of Sciences of the United States of America.

[16]  W. Briels,et al.  THE CALCULATION OF FREE-ENERGY DIFFERENCES BY CONSTRAINED MOLECULAR-DYNAMICS SIMULATIONS , 1998 .

[17]  P. Kollman,et al.  A Second Generation Force Field for the Simulation of Proteins, Nucleic Acids, and Organic Molecules , 1995 .

[18]  C. Orengo,et al.  Stability-activity tradeoffs constrain the adaptive evolution of RubisCO , 2014, Proceedings of the National Academy of Sciences.

[19]  F. C. Hartman,et al.  Mechanism of Rubisco: The Carbamate as General Base. , 1998, Chemical reviews.

[20]  G. Farquhar,et al.  D2O solvent isotope effects suggest uniform energy barriers in ribulose-1,5-bisphosphate carboxylase/oxygenase catalysis. , 2013, Biochemistry.

[21]  P. Saalfrank,et al.  Quantum chemical modeling of the kinetic isotope effect of the carboxylation step in RuBisCO , 2012, Journal of Molecular Modeling.

[22]  Eric F Darve,et al.  Calculating free energies using average force , 2001 .

[23]  J. Gready,et al.  CO(2) fixation by Rubisco: computational dissection of the key steps of carboxylation, hydration, and C-C bond cleavage. , 2001, Journal of the American Chemical Society.

[24]  MNDO method for calculations of magnesium clusters , 1994 .

[25]  Tjerk P. Straatsma,et al.  NWChem: A comprehensive and scalable open-source solution for large scale molecular simulations , 2010, Comput. Phys. Commun..

[26]  Natércia F. Brás,et al.  Application of quantum mechanics/molecular mechanics methods in the study of enzymatic reaction mechanisms , 2017 .

[27]  Tsvi Tlusty,et al.  Cross-species analysis traces adaptation of Rubisco toward optimality in a low-dimensional landscape , 2010, Proceedings of the National Academy of Sciences.

[28]  U. Rothlisberger,et al.  Mixed Quantum Mechanical/Molecular Mechanical Molecular Dynamics Simulations of Biological Systems in Ground and Electronically Excited States. , 2015, Chemical reviews.

[29]  J. Gready,et al.  Coupled semiempirical molecular orbital and molecular mechanics model (QM/MM) for organic molecules in aqueous solution , 1997, J. Comput. Chem..

[30]  J. Gready,et al.  Directions for Optimization of Photosynthetic Carbon Fixation: RuBisCO's Efficiency May Not Be So Constrained After All , 2018, Front. Plant Sci..

[31]  P. C. Hariharan,et al.  The influence of polarization functions on molecular orbital hydrogenation energies , 1973 .

[32]  M Mezei,et al.  Free Energy Simulations a , 1986, Annals of the New York Academy of Sciences.

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

[34]  W. D. Otter,et al.  Thermodynamic integration of the free energy along a reaction coordinate in Cartesian coordinates , 2000 .

[35]  G. Tcherkez.  Modelling the reaction mechanism of ribulose-1,5-bisphosphate carboxylase/oxygenase and consequences for kinetic parameters. , 2013, Plant, cell & environment.

[36]  D. Jordan,et al.  Species variation in kinetic properties of ribulose 1,5-bisphosphate carboxylase/oxygenase. , 1983, Archives of biochemistry and biophysics.

[37]  G. Lorimer,et al.  Reaction intermediate partitioning by ribulose-bisphosphate carboxylases with differing substrate specificities. , 1986, The Journal of biological chemistry.

[38]  G. Farquhar,et al.  Determining RuBisCO activation kinetics and other rate and equilibrium constants by simultaneous multiple non-linear regression of a kinetic model. , 2006, Journal of experimental botany.

[39]  J. Gready,et al.  Coupled semiempirical quantum mechanics and molecular mechanics (QM/MM) calculations on the aqueous solvation free energies of ionized molecules , 1999, J. Comput. Chem..

[40]  G. Ciccotti,et al.  Numerical Integration of the Cartesian Equations of Motion of a System with Constraints: Molecular Dynamics of n-Alkanes , 1977 .