Quantum chemical modeling of the kinetic isotope effect of the carboxylation step in RuBisCO

Possible Steps of the Carboxylation of Ribulose-1,5-biphosphate from Intermediates: 2,3-Enediol versus 1,2-Enol

Ribulose 1,5-bisphosphate (RuBP) undergoes enolization to initiate fixation of atmospheric carbon dioxide in the plant carbon cycle. The known model assumes the binding of RuBP to the Rubisco active site with the subsequent formation of 2,3-enediol (2,3,4-trihydroxypent-2-ene-1,5-diyl diphosphate). In the present study, it is assumed that 1,2-enol (2,3,4-trihydroxypent-1-ene-1,5-diyl diphosphate) can be formed in the enolization step to initiate the carboxylation reaction. We have used Kohn–Sham density functional theory on WB97X-D3/Def2-TZVP levels to compare the reaction barriers in the two ways. We considered the pathways of carboxylation of 1/2-ene (mono/di)ol via the C1 and C2 carbons without taking into account the binding of RuBP to the magnesium ion. Calculations of Gibbs free energies confirm the equal probability of the formation of 2,3-enediol and 1,2-enol. Quantum–chemical modeling of enolization and carboxylation reactions supports the important role of the bridging water molecule and diphosphate groups, which provide proton transfer and lower reaction barriers. The results show that carbon dioxide fixation can occur without a magnesium ion, and binding with C1 can have a lower barrier (~12 kcal/mol) than with C2 (~23 kcal/mol).

Theoretical analysis of the kinetic isotope effect on carboxylation in RubisCO

Ribulose-1,5-bisphosphate carboxylase/oxygenase (RubisCO) fixes atmospheric carbon dioxide into bioavailable sugar molecules. It is also well known that a kinetic isotope effect (KIE; CO₂ carbon atoms) accompanies the carboxylation process. To describe the reaction and the KIE α, two different types of molecular dynamics (MD) simulations (ab initio MD and classical MD) have been performed with an Own N-layered Integrated molecular Orbitals and molecular Mechanics (ONIOM)-hybrid model. A channel structure for CO₂ transport has been observed during the MD simulation in RubisCO, and assuming the reaction path from the inlet to the product through the coordinate complex with Mg²⁺ , simulations have been performed on several molecular configuration models fixing several distances between CO₂ and ribulose-1,5-bisphosphate along the channel. Free energy analysis and diffusion coefficient analysis have been evaluated for different phases of the process. It is confirmed that the isotopic fractionation effect for CO₂ containing either ¹³ C or ¹² C would appear through the transiting path in the channel structure identified in RubisCO. The estimated isotope fractionation constant was quite close to the experimental value.

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, CO₂ . Hydration of the RuBP.CO₂ 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.CO₂ 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 CO₂ ; both hydration and CC bond scission require early protonation of CO₂ in the RuBP.CO₂ 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 K_m (=K_c ) with turnover to product via bond scission concerted with stereospecific protonation consistent with the observed catalytic rate. © 2018 Wiley Periodicals, Inc.