Redefinition of rubisco carboxylase reaction reveals origin of water for hydration and new roles for active-site residues

Red Rubiscos and opportunities for engineering green plants

Nature's vital, but notoriously inefficient, CO2-fixing enzyme Rubisco often limits the growth of photosynthetic organisms including crop species. Form I Rubiscos comprise eight catalytic large subunits and eight auxiliary small subunits and can be classified into two distinct lineages-'red' and 'green'. While red-type Rubiscos (Form IC and ID) are found in rhodophytes, their secondary symbionts, and certain proteobacteria, green-type Rubiscos (Form IA and IB) exist in terrestrial plants, chlorophytes, cyanobacteria, and other proteobacteria. Eukaryotic red-type Rubiscos exhibit desirable kinetic properties, namely high specificity and high catalytic efficiency, with certain isoforms outperforming green-type Rubiscos. However, it is not yet possible to functionally express a high-performing red-type Rubisco in chloroplasts to boost photosynthetic carbon assimilation in green plants. Understanding the molecular and evolutionary basis for divergence between red- and green-type Rubiscos could help us to harness the superior CO2-fixing power of red-type Rubiscos. Here we review our current understanding about red-type Rubisco distribution, biogenesis, and sequence-structure, and present opportunities and challenges for utilizing red-type Rubisco kinetics towards crop improvements.

Theoretical Study on the Kinetics of the Rubisco Carboxylase Reaction by a Model Based on Quantum Chemistry and Absolute Reaction Rate Theory

The rate of the Rubisco carboxylase reaction is evaluated by statistical mechanics and hybrid density functional theory (DFT). The Rubisco molecular model given by Kannappan et al. was modified and used in the present calculation. The activation energies of CO₂ addition reaction, H₂O addition reaction, C2-C3 bond scission, and C2 protonation are estimated. We calculated the turnover number (TON) for each of the four reaction steps based on a revised absolute reaction rate theory, which became applicable to soft matter reactions. The molecular parameters used in TON calculations were obtained by DFT calculations. The TON of the total Rubisco reaction was finally evaluated using rate equations. The calculation in a vacuum gave the total TON to be around 5 × 10^-5, which was much lower than the experimental value. The DFT calculation in water solvent gave the total TON to be around 0.1, which agreed reasonably well with experimentally reported values (∼2.71). The rate-limiting process was the scission reaction. The present calculation showed that both the phosphate groups in the substrate accelerate each reaction step. The present calculation showed that a more comprehensive molecular model including enolization and quantum chemical methods is necessary to make a more precise reaction model including the irreversibility of some reactions.

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

The Mechanism of Rubisco Catalyzed Carboxylation Reaction: Chemical Aspects Involving Acid-Base Chemistry and Functioning of the Molecular Machine

In recent years, a great deal of attention has been paid by the scientific community to improving the efficiency of photosynthetic carbon assimilation, plant growth and biomass production in order to achieve a higher crop productivity. Therefore, the primary carboxylase enzyme of the photosynthetic process Rubisco has received considerable attention focused on many aspects of the enzyme function including protein structure, protein engineering and assembly, enzyme activation and kinetics. Based on its fundamental role in carbon assimilation Rubisco is also targeted by the CO2-fertilization effect, which is the increased rate of photosynthesis due to increasing atmospheric CO2-concentration. The aim of this review is to provide a framework, as complete as possible, of the mechanism of the RuBP carboxylation/hydration reaction including description of chemical events occurring at the enzyme “activating” and “catalytic” sites (which involve Broensted acid-base reactions) and the functioning of the complex molecular machine. Important research results achieved over the last few years providing substantial advancement in understanding the enzyme functioning will be discussed.

Molecular Evolution of rbcL in Orthotrichales (Bryophyta): Site Variation, Adaptive Evolution, and Coevolutionary Patterns of Amino Acid Replacements

Molecular evolution of the large subunit of the RuBisCO enzyme is understudied in early diverging land plants. These groups show morphological and eco-physiological adaptations to the uneven and intermittent distribution of water in the terrestrial environment. This might have prompted a continuous fine-tuning of RuBisCO under a selective pressure modifying the species-specific optima for photosynthesis in contrasting microdistributions and environmental niches. To gain a better insight into the molecular evolution of RuBisCO large subunits, the aim of this study was to assess the pattern of evolutionary change in the amino acid residues in a monophyletic group of Bryophyta (Orthotrichaceae). Tests for positive, neutral, or purifying selection at the amino acid level were assessed by comparing rates (ω) of non-synonymous (d_N) and synonymous (d_S) nucleotide substitutions along a Maximum Likelihood phylogenetic tree. Molecular adaptation tests using likelihood ratio tests, reconstruction of ancestral amino acid sites, and intra-protein coevolution analyses were performed. Variable amino acid sites (39) were unevenly distributed across the LSU. The residues are located on rbcL sites that are highly variable in higher plants and close to key regions implying dimer-dimer (L₂L₂), RuBisCO-activase interactions, and conformational functions during catalysis. Ten rbcL sites (32, 33, 91, 230, 247, 251, 255, 424, 449 and 475) have been identified by the Bayesian Empirical Bayes inference to be under positive selection and under adaptive evolution under the M8 model. The pattern of amino acid variation suggests that it is not lineage specific, but rather representative of a case of convergent evolution, suggesting recurrent changes that potentially favor the same amino acid substitutions that are likely optimized the RuBisCO activity.

The Coevolution of RuBisCO, Photorespiration, and Carbon Concentrating Mechanisms in Higher Plants

Ribulose-1,5-bisphosphate (RuBP) carboxylase/oxygenase (RuBisCO) is the carbon-fixing enzyme present in most photosynthetic organisms, converting CO₂ into organic matter. Globally, photosynthetic efficiency in terrestrial plants has become increasingly challenged in recent decades due to a rapid increase in atmospheric CO₂ and associated changes toward warmer and dryer environments. Well adapted for these new climatic conditions, the C₄ photosynthetic pathway utilizes carbon concentrating mechanisms to increase CO₂ concentrations surrounding RuBisCO, suppressing photorespiration from the oxygenase catalyzed reaction with O₂. The energy efficiency of C₃ photosynthesis, from which the C₄ pathway evolved, is thought to rely critically on an uninterrupted supply of chloroplast CO₂. Part of the homeostatic mechanism that maintains this constancy of supply involves the CO₂ produced as a byproduct of photorespiration in a negative feedback loop. Analyzing the database of RuBisCO kinetic parameters, we suggest that in genera (Flaveria and Panicum) for which both C₃ and C₄ examples are available, the C₄ pathway evolved only from C₃ ancestors possessing much lower than the average carboxylase specificity relative to that of the oxygenase reaction (S _C/O=S _C/S _O), and hence, the higher CO₂ levels required for development of the photorespiratory CO₂ pump (C₂ photosynthesis) essential in the initial stages of C₄ evolution, while in the later stage (final optimization phase in the Flaveria model) increased CO₂ turnover may have occurred, which would have been supported by the higher CO₂ levels. Otherwise, C₄ RuBisCO kinetic traits remain little changed from the ancestral C₃ species. At the opposite end of the spectrum, C₃ plants (from Limonium) with higher than average S _C/O, which may be associated with the ability of increased CO₂, relative to O₂, affinity to offset reduced photorespiration and chloroplast CO₂ levels, can tolerate high stress environments. It is suggested that, instead of inherently constrained by its kinetic mechanism, RuBisCO possesses the extensive kinetic plasticity necessary for adaptation to changes in photorespiration that occur in the homeostatic regulation of CO₂ supply under a broad range of abiotic environmental conditions.

Ribulose 1,5-bisphosphate carboxylase/oxygenase activates O2 by electron transfer

Ribulose 1,5-bisphosphate carboxylase/oxygenase (Rubisco) is the cornerstone of atmospheric CO2 fixation by the biosphere. It catalyzes the addition of CO2 onto enolized ribulose 1,5-bisphosphate (RuBP), producing 3-phosphoglycerate which is then converted to sugars. The major problem of this reaction is competitive O2 addition, which forms a phosphorylated product (2-phosphoglycolate) that must be recycled by a series of biochemical reactions (photorespiratory metabolism). However, the way the enzyme activates O2 is still unknown. Here, we used isotope effects (with 2H, 25Mg, and 18O) to monitor O2 activation and assess the influence of outer sphere atoms, in two Rubisco forms of contrasted O2/CO2 selectivity. Neither the Rubisco form nor the use of solvent D2O and deuterated RuBP changed the 16O/18O isotope effect of O2 addition, in clear contrast with the 12C/13C isotope effect of CO2 addition. Furthermore, substitution of light magnesium (24Mg) by heavy, nuclear magnetic 25Mg had no effect on O2 addition. Therefore, outer sphere protons have no influence on the reaction and direct radical chemistry (intersystem crossing with triplet O2) does not seem to be involved in O2 activation. Computations indicate that the reduction potential of enolized RuBP (near 0.49 V) is compatible with superoxide (O2•-) production, must be insensitive to deuteration, and yields a predicted 16O/18O isotope effect and energy barrier close to observed values. Overall, O2 undergoes single electron transfer to form short-lived superoxide, which then recombines to form a peroxide intermediate.

Electronic structure benchmark calculations of CO2 fixing elementary chemical steps in RuBisCO using the projector-based embedding approach

Ribulose 1,5-bisphosphate carboxylase-oxygenase (RuBisCO) is the main enzyme involved in atmospheric carbon dioxide (CO₂ ) fixation in the biosphere. This enzyme catalyzes a set of five chemical steps that take place in the same active-site within magnesium (II) coordination sphere. Here, a set of electronic structure benchmark calculations have been carried out on a reaction path proposed by Gready et al. by means of the projector-based embedding approach. Activation and reaction energies for all main steps catalyzed by RuBisCO have been calculated at the MP2, SCS-MP2, CCSD, and CCSD(T)/aug-cc-pVDZ and cc-pVDZ levels of theory. The treatment of the magnesium cation with post-HF methods is explored to determine the nature of its involvement in the mechanism. With the high-level ab initio values as a reference, we tested the performance of a set of density functional theory (DFT) exchange-correlation (xc) functionals in reproducing the reaction energetics of RuBisCO carboxylase activity on a set of model fragments. Different DFT xc-functionals show large variation in activation and reaction energies. Activation and reaction energies computed at the B3LYP level are close to the reference SCS-MP2 results for carboxylation, hydration and protonation reactions. However, for the carbon-carbon bond dissociation reaction, B3LYP and other functionals give results that differ significantly from the ab initio reference values. The results show the applicability of the projector-based embedding approach to metalloenzymes. This technique removes the uncertainty associated with the selection of different DFT xc-functionals and so can overcome some of inherent limitations of DFT calculations, complementing, and potentially adding to modeling of enzyme reaction mechanisms with DFT methods.

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.

Kohn-Sham Density Functional Calculations Reveal Proton Wires in the Enolization and Carboxylase Reactions Catalyzed by Rubisco

Ribulose 1,5-bisphosphate (RuBP) carboxylase-oxygenase (Rubisco) plays a fundamental role in the carbon cycle by fixing the atmospheric CO₂ used in photosynthesis. Rubisco is all the more remarkable because it must catalyze some difficult multistep reaction chemistry involving proton transfers within the one active site. In the present study, we have used Kohn-Sham density functional theory at the B3LYP/6-31G* level with basis set superposition error and dispersion corrections (B3LYP-gCP-D3) to examine the possibility that the proton transfers can take place through molecular wires (including active-site water molecules) via the classical Grotthuss proton-shuttle mechanism. The results support an essential role for water molecules found in the crystal structures of Rubisco complexes as facilitators of proton transport in all the rate-limiting (catalytic) reaction steps through a network of short proton wires within the Rubisco active site. We suggest that completion of the initial product turnover (cycle) requires two excess protons produced in the initial carbamylation that is required for Rubisco activation. By use of proton wires, a large number of reaction steps may be accommodated within a single active site without necessitating the input of excessive conformational strain energy arising from the movement of residue side chains into positions where direct protonation of substrates can occur. The involvement of the identified types of proton wires in the kinetic mechanism is capable of providing a unique explanation for various experimental observations, including deuterium isotope effects and the results of site-directed mutagenesis experiments, and may thus provide a realistic solution to the problem of Rubisco's challenging chemistry.

Electronic structure benchmark calculations of inorganic and biochemical carboxylation reactions

Carboxylation reactions represent a very special class of chemical reactions that is characterized by the presence of a carbon dioxide (CO₂ ) molecule as reactive species within its global chemical equation. These reactions work as fundamental gear to accomplish the CO₂ fixation and thus to build up more complex molecules through different technological and biochemical processes. In this context, a correct description of the CO₂ electronic structure turns out to be crucial to study the chemical and electronic properties associated with this kind of reactions. Here, a systematic study of CO₂ electronic structure and its contribution to different carboxylation reaction electronic energies has been carried out by means of several high-level ab initio post-Hartree Fock (post-HF) and density functional theory (DFT) calculations for a set of biochemistry and inorganic systems. We have found that for a correct description of the CO₂ electronic correlation energy it is necessary to include post-CCSD(T) contributions (beyond the gold standard). These high-order excitations are required to properly describe the interactions of the four π-electrons associated with the two degenerated π-molecular orbitals of the CO₂ molecule. Likewise, our results show that in some reactions it is possible to obtain accurate reaction electronic energy values with computationally less demanding methods when the error in the electronic correlation energy compensates between reactants and products. Furthermore, the provided post-HF reference values allowed to validating different DFT exchange-correlation functionals combined with different basis sets for chemical reactions that are relevant in biochemical CO₂ fixing enzymes. © 2019 Wiley Periodicals, Inc.

Mechanism of Oxygenase-Pathway Reactions Catalyzed by Rubisco from Large-Scale Kohn-Sham Density Functional Calculations

Ribulose 1,5-bisphosphate carboxylase-oxygenase (Rubisco) is the primary carbon-fixing enzyme in photosynthesis, fixing CO₂ to a 5-carbon sugar, RuBP, in a series of five reactions. However, it also catalyzes an oxygenase reaction by O₂ addition to the same enolized RuBP substrate in an analogous reaction series in the same active site, producing a waste product and loss of photosynthetic efficiency. Starting from RuBP, the reactions are enolization to the enediolate form, addition of CO₂ or O₂ to form the carboxy or peroxo adduct, hydration to form a gemdiolate, scission of the C2-C3 bond of the original RuBP, and stereospecific or nonstereospecific protonation to form two molecules of the 3-carbon PGA product, or one molecule of PGA, one of 2-carbon PG (waste product), and one water molecule. Reducing the loss of efficiency from the oxygenase reaction is an attractive means to increase crop productivity. However, lack of understanding of key aspects of the catalytic mechanisms for both the carboxylase and oxygenase reactions, particularly those involving proton exchanges and roles of water molecules, has stymied efforts at re-engineering Rubisco to reduce losses from the oxygenation reaction. As the stable form of molecular oxygen is the triplet biradical state (³O₂), its reaction with near-universal singlet-state molecules is formally spin forbidden. Although in oxygenase enzymes, ³O₂ activation is usually achieved by one-electron transfers using transition-metal ions or organic cofactors, recently, cofactor-less oxygenases in which the substrate itself is the source of the electron for ³O₂ activation have been identified, but in all such cases an aromatic ring stabilizes the substrate's negative charge. Here we present the first large-scale Kohn-Sham density functional theory study of the reaction mechanism of the Rubisco oxygenase pathway. First, we show that the enediolate substrate complexed to Mg²⁺ and its ligands extends the region for charge delocalization and stabilization of its negative charge to allow formation of a caged biradical enediolate-O₂ complex. Thus, Rubisco is a unique type of oxygenase without precedent in the literature. Second, for the O₂ addition to proceed to the singlet peroxo-adduct intermediate, the system must undergo an intersystem crossing. We found that the presence of protonated LYS334 is required to stabilize this intermediate and that both factors (strongly stabilized anion and protonated LYS334) facilitate a barrier-less activation of ³O₂. This finding supports our recent proposal that deoxygenation, that is, reversal of gas binding, is possible. Third, as neither CO₂ nor O₂ binds to the enzyme, our findings support the proposal from our recent carboxylase study that the observed K_C or K_O (Michaelis-Menten constants) in the steady-state kinetics reflect the respective adducts, carboxy or peroxo. Fourth, after computing hydration pathways with water addition both syn and anti to C3, we found, in contrast to the results of our carboxylation study indicating anti addition, that in the oxygenation reaction only syn-hydration is capable of producing a stable gemdiolate that facilitates the rate-limiting C2-C3 bond scission to final products. Fifth, we propose that an excess proton we previously found was required in the carboxylation reaction for activating the C2-C3 bond scission is utilized in the oxygenation reaction for the required elimination of a water molecule. In summary, despite its oxygenase handicap, Rubisco's success in directing 75% of its substrate through the carboxylation pathway can be considered impressively effective. Although native C3 Rubiscos are in a fix with unwanted activity of ³O₂ hampering its primary carboxylase function, mechanistic differences presented here with findings in our recent carboxylase study for both the gas-addition and subsequent reactions provide some clues as to how creative Rubisco re-engineering may offer a solution to reducing the oxygenase activity.

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.

Directed -in vitro- evolution of Precambrian and extant Rubiscos

Rubisco is an ancient, catalytically conserved yet slow enzyme, which plays a central role in the biosphere's carbon cycle. The design of Rubiscos to increase agricultural productivity has hitherto relied on the use of in vivo selection systems, precluding the exploration of biochemical traits that are not wired to cell survival. We present a directed -in vitro- evolution platform that extracts the enzyme from its biological context to provide a new avenue for Rubisco engineering. Precambrian and extant form II Rubiscos were subjected to an ensemble of directed evolution strategies aimed at improving thermostability. The most recent ancestor of proteobacteria -dating back 2.4 billion years- was uniquely tolerant to mutagenic loading. Adaptive evolution, focused evolution and genetic drift revealed a panel of thermostable mutants, some deviating from the characteristic trade-offs in CO2-fixing speed and specificity. Our findings provide a novel approach for identifying Rubisco variants with improved catalytic evolution potential.

Directions for Optimization of Photosynthetic Carbon Fixation: RuBisCO's Efficiency May Not Be So Constrained After All

The ubiquitous enzyme Ribulose 1,5-bisphosphate carboxylase-oxygenase (RuBisCO) fixes atmospheric carbon dioxide within the Calvin-Benson cycle that is utilized by most photosynthetic organisms. Despite this central role, RuBisCO's efficiency surprisingly struggles, with both a very slow turnover rate to products and also impaired substrate specificity, features that have long been an enigma as it would be assumed that its efficiency was under strong evolutionary pressure. RuBisCO's substrate specificity is compromised as it catalyzes a side-fixation reaction with atmospheric oxygen; empirical kinetic results show a trend to tradeoff between relative specificity and low catalytic turnover rate. Although the dominant hypothesis has been that the active-site chemistry constrains the enzyme's evolution, a more recent study on RuBisCO stability and adaptability has implicated competing selection pressures. Elucidating these constraints is crucial for directing future research on improving photosynthesis, as the current literature casts doubt on the potential effectiveness of site-directed mutagenesis to improve RuBisCO's efficiency. Here we use regression analysis to quantify the relationships between kinetic parameters obtained from empirical data sets spanning a wide evolutionary range of RuBisCOs. Most significantly we found that the rate constant for dissociation of CO₂ from the enzyme complex was much higher than previous estimates and comparable with the corresponding catalytic rate constant. Observed trends between relative specificity and turnover rate can be expressed as the product of negative and positive correlation factors. This provides an explanation in simple kinetic terms of both the natural variation of relative specificity as well as that obtained by reported site-directed mutagenesis results. We demonstrate that the kinetic behaviour shows a lesser rather than more constrained RuBisCO, consistent with growing empirical evidence of higher variability in relative specificity. In summary our analysis supports an explanation for the origin of the tradeoff between specificity and turnover as due to competition between protein stability and activity, rather than constraints between rate constants imposed by the underlying chemistry. Our analysis suggests that simultaneous improvement in both specificity and turnover rate of RuBisCO is possible.

Structure of Rubisco from Arabidopsis thaliana in complex with 2-carboxyarabinitol-1,5-bisphosphate

The crystal structure of ribulose-1,5-bisphosphate carboxylase/oxygenase (Rubisco) from Arabidopsis thaliana is reported at 1.5 Å resolution. In light of the importance of A. thaliana as a model organism for understanding higher plant biology, and the pivotal role of Rubisco in photosynthetic carbon assimilation, there has been a notable absence of an A. thaliana Rubisco crystal structure. A. thaliana Rubisco is an L8S8 hexadecamer comprising eight plastome-encoded catalytic large (L) subunits and eight nuclear-encoded small (S) subunits. A. thaliana produces four distinct small-subunit isoforms (RbcS1A, RbcS1B, RbcS2B and RbcS3B), and this crystal structure provides a snapshot of A. thaliana Rubisco containing the low-abundance RbcS3B small-subunit isoform. Crystals were obtained in the presence of the transition-state analogue 2-carboxy-D-arabinitol-1,5-bisphosphate. A. thaliana Rubisco shares the overall fold characteristic of higher plant Rubiscos, but exhibits an interesting disparity between sequence and structural relatedness to other Rubisco isoforms. These results provide the structural framework to understand A. thaliana Rubisco and the potential catalytic differences that could be conferred by alternative A. thaliana Rubisco small-subunit isoforms.

RubisCO Early Oxygenase Activity: A Kinetic and Evolutionary Perspective

RubisCO (D-ribulose 1,5-bisphosphate carboxylase/oxygenase) is Earth's main enzyme responsible for CO₂ fixation via carboxylation of ribulose-1,5-bisphosphate (RuBP) into organic matter. Besides the carboxylation reaction, RubisCO also catalyzes the oxygenation of RuBP by O₂ , which is probably as old as its carboxylation properties. Based on molecular phylogeny, the occurrence of the reactive oxygen species (ROS)-removing system and kinetic properties of different RubisCO forms, we postulated that RubisCO oxygenase activity appeared in local microoxic areas, yet before the appearance of oxygenic photosynthesis. Here, in reviewing the literature, we present a novel hypothesis: the RubisCO early oxygenase activity hypothesis. This hypothesis may be compared with the exaptation hypothesis, according to which latent RubisCO oxygenase properties emerged later during the oxygenation of the Earth's atmosphere. The reconstruction of ancestral RubisCO forms using ancestral sequence reconstruction (ASR) techniques, as a promising way for testing of RubisCO early oxygenase activity hypothesis, is presented.

The Diverse AAA+ Machines that Repair Inhibited Rubisco Active Sites

Gaseous carbon dioxide enters the biosphere almost exclusively via the active site of the enzyme ribulose 1,5-bisphosphate carboxylase/oxygenase (Rubisco). This highly conserved catalyst has an almost universal propensity to non-productively interact with its substrate ribulose 1,5-bisphosphate, leading to the formation of dead-end inhibited complexes. In diverse autotrophic organisms this tendency has been counteracted by the recruitment of dedicated AAA+ (ATPases associated with various cellular activities) proteins that all use the energy of ATP hydrolysis to remodel inhibited Rubisco active sites leading to release of the inhibitor. Three evolutionarily distinct classes of these Rubisco activases (Rcas) have been discovered so far. Green and red-type Rca are mostly found in photosynthetic eukaryotes of the green and red plastid lineage respectively, whereas CbbQO is associated with chemoautotrophic bacteria. Ongoing mechanistic studies are elucidating how the various motors are utilizing both similar and contrasting strategies to ultimately perform their common function of cracking the inhibited Rubisco active site. The best studied mechanism utilized by red-type Rca appears to involve transient threading of the Rubisco large subunit C-terminal peptide, reminiscent of the action performed by Clp proteases. As well as providing a fascinating example of convergent molecular evolution, Rca proteins can be considered promising crop-improvement targets. Approaches aiming to replace Rubisco in plants with improved enzymes will need to ensure the presence of a compatible Rca protein. The thermolability of the Rca protein found in crop plants provides an opportunity to fortify photosynthesis against high temperature stress. Photosynthesis also appears to be limited by Rca when light conditions are fluctuating. Synthetic biology strategies aiming to enhance the autotrophic CO₂ fixation machinery will need to take into consideration the requirement for Rubisco activases as well as their properties.

Magnesium-based systems for carbon dioxide capture, storage and recycling: from leaves to synthetic nanostructured materials

Magnesium is used as leitmotif in this review in order to explore the systems involved in natural and artificial CO₂ cycles.

Highly sensitive detection of carbon dioxide by a pyrimido[1,2-a]benzimidazole derivative: combining experimental and theoretical studies

In the present paper, a “light-up” chemsensor with a high specificity for carbon dioxide detection using a pyrimido[1,2-a]benzimidazole derivative (P1H) in liquid media has been developed.

Recent advances in riboflavin biosynthesis

Riboflavin is biosynthesized from GTP and ribulose 5-phosphate. Whereas the early reactions conducing to 5-amino-6-ribitylamino-2,4(1H,3H)-pyrimidinedione 5'-phosphate show significant taxonomic variation, the subsequent reaction steps are universal in all taxonomic kingdoms. With the exception of a hitherto elusive phosphatase, all enzymes of the pathway have been characterized in some detail at the structural and mechanistic level. Some of the pathway enzymes (GTP cycloyhdrolase II, 3,4-dihydroxy-2-butanone 4-phosphate synthase, riboflavin synthase) have exceptionally complex reaction mechanisms. The commercial production of the vitamin is now entirely based on highly productive fermentation processes. Due to their absence in animals, the pathway enzymes are potential targets for the development of novel anti-infective drugs.

Modelling the reaction mechanism of ribulose-1,5-bisphosphate carboxylase/oxygenase and consequences for kinetic parameters

Although ribulose-1,5-bisphosphate carboxylase/oxygenase (Rubisco) was discovered nearly 60 years ago, the associated chemical mechanism of the reaction is still incompletely understood. The catalytic cycle consists of four major steps: ribulose-1,5-bisphosphate binding, enolization, CO₂ or O₂ addition and hydration, and cleavage of the intermediate. The use of individual rate constants for these elemental steps yields mathematical expressions for usual kinetic constants (k(cat), K(m)), CO₂ versus O₂ specificity (S(c/o)) as well as other chemical parameters such as the ¹²C/¹³C isotope effect. That said, most of them are not simple and thus the interpretation of experimental and observed values of kcat , Km and Sc/o may be more complicated than expected. That is, Rubisco effective catalysis depends on several kinetic parameters that are influenced by both the biological origin and the cellular medium (which, in turn, can vary with environmental conditions). In this brief review, we present the basic model of Rubisco kinetics and describe how subtle biochemical changes (which may have occurred along Evolution) can easily modify Rubisco catalysis.

CO2 on a tightrope: stabilization, room-temperature decarboxylation, and sodium-induced carboxylate migration

A sterically shielded 3-substituted zwitterionic N,N-dimethylisotryptammonium carboxylate has been synthesized by consecutive chemoselective double alkylation of indole. The carboxylate undergoes a quantitative and unusually facile decarboxylation in dimethyl sulfoxide (DMSO) or dimethyl formamide (DMF) at room temperature. The breaking of a nearly equidistant hydrogen bond by solvent molecules initiates heterolytic C-C cleavage. The decarboxylation rate decreases with increasing CO(2) partial pressure, proving the competitiveness of protonation and re-carboxylation of the carbanionic intermediate. Corresponding spiro compounds containing silylene and stannylene moieties show high thermal stability. Addition of an excess of methyllithium to the sodium salt triggers a reaction sequence comprising a deprotonation, carboxylate transfer, and nucleophilic trapping of the rearranged carboxylate by another equivalent of methyllithium. Hydrolytic work-up of the geminal diolate leads to an acetyl product. The role of the sodium counterion and the mechanism of the rearrangement have been unraveled by deuteration experiments.

Evolutionary context for understanding and manipulating plant responses to past, present and future atmospheric [CO2]

Variation in atmospheric [CO(2)] is a prominent feature of the environmental history over which vascular plants have evolved. Periods of falling and low [CO(2)] in the palaeo-record appear to have created selective pressure for important adaptations in modern plants. Today, rising [CO(2)] is a key component of anthropogenic global environmental change that will impact plants and the ecosystem goods and services they deliver. Currently, there is limited evidence that natural plant populations have evolved in response to contemporary increases in [CO(2)] in ways that increase plant productivity or fitness, and no evidence for incidental breeding of crop varieties to achieve greater yield enhancement from rising [CO(2)]. Evolutionary responses to elevated [CO(2)] have been studied by applying selection in controlled environments, quantitative genetics and trait-based approaches. Findings to date suggest that adaptive changes in plant traits in response to future [CO(2)] will not be consistently observed across species or environments and will not be large in magnitude compared with physiological and ecological responses to future [CO(2)]. This lack of evidence for strong evolutionary effects of elevated [CO(2)] is surprising, given the large effects of elevated [CO(2)] on plant phenotypes. New studies under more stressful, complex environmental conditions associated with climate change may revise this view. Efforts are underway to engineer plants to: (i) overcome the limitations to photosynthesis from today's [CO(2)] and (ii) benefit maximally from future, greater [CO(2)]. Targets range in scale from manipulating the function of a single enzyme (e.g. Rubisco) to adding metabolic pathways from bacteria as well as engineering the structural and functional components necessary for C(4) photosynthesis into C(3) leaves. Successfully improving plant performance will depend on combining the knowledge of the evolutionary context, cellular basis and physiological integration of plant responses to varying [CO(2)].

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

Ribulose-1,5-bisphosphate carboxylase/oxygenase (RuBisCO), the most important enzyme for the assimilation of carbon into biomass, features a well-known isotope effect with regards to the CO(2) carbon atom. This kinetic isotope effect α = k(12)/k(13) for the carboxylation step of the RuBisCO reaction sequence, and its microscopic origin, was investigated with the help of cluster models and quantum chemical methods [B3LYP/6-31G(d,p)]. We use a recently proposed model for the RuBisCO active site, in which a water molecule remains close to the reaction center during carboxylation of ribulose-1,5-bisphosphate [B. Kannappan, J.E. Gready, J. Am. Chem. Soc. 130 (2008), 15063]. Alternative active-site models and/or computational approaches were also tested. An isotope effect alpha for carboxylation is found, which is reasonably close to the one measured for the overall reaction, and which originates from a simple frequency shift of the bending vibration of (12)CO(2) compared to (13)CO(2). The latter is the dominant mode for the product formation at the transition state.

Slow deactivation of ribulose 1,5-bisphosphate carboxylase/oxygenase elucidated by mathematical models

Ribulose 1,5-bisphosphate carboxylase/oxygenase (RuBisCO) is the key enzyme of the Calvin cycle, catalyzing the fixation of inorganic carbon dioxide to organic sugars. Unlike most enzymes, RuBisCO is extremely slow, substrate unspecific, and catalyzes undesired side-reactions, which are considered to be responsible for the slow deactivation observed in vitro, a phenomenon known as fallover. Despite the fact that amino acid sequences and the 3D structures of RuBisCO from a variety of species are known, the precise molecular mechanisms for the various side reactions are still unclear. In the present study, we investigate the kinetic properties of RuBisCO using mathematical models. Initially, we formulate a minimal model that quantitatively reflects the kinetic behavior of RuBisCOs from different organisms. By relating rate parameters for single molecular steps to experimentally determined K(m) and V(max) values, we can examine mechanistic differences among species. The minimal model further demonstrates that two inhibitor producing side reactions are sufficient to describe experimentally determined fallover kinetics. To explain the observed kinetics of the limited capacity of RuBisCO to accept xylulose 1,5-bisphosphate as substrate, the inclusion of other side reactions is necessary. Our model results suggest a yet undescribed alternative enolization mechanism that is supported by the molecular structure. Taken together, the presented models serve as a theoretical framework to explain a wide range of observed kinetic properties of RuBisCOs derived from a variety of species. Thus, we can support hypotheses about molecular mechanisms and can systematically compare enzymes from different origins.

Substrate-induced assembly of Methanococcoides burtonii D-ribulose-1,5-bisphosphate carboxylase/oxygenase dimers into decamers

Like many enzymes, the biogenesis of the multi-subunit CO(2)-fixing enzyme ribulose-1,5-bisphosphate (RuBP) carboxylase/oxygenase (Rubisco) in different organisms requires molecular chaperones. When expressed in Escherichia coli, the large (L) subunits of the Rubisco from the archaeabacterium Methanococcoides burtonii assemble into functional dimers (L(2)). However, further assembly into pentamers of L(2) (L(10)) occurs when expressed in tobacco chloroplasts or E. coli producing RuBP. In vitro analyses indicate that the sequential assembly of L(2) into L(10) (via detectable L(4) and L(6) intermediates) occurs without chaperone involvement and is stimulated by protein rearrangements associated with either the binding of substrate RuBP, the tight binding transition state analog carboxyarabinitol-1,5-bisphosphate, or inhibitory divalent metal ions within the active site. The catalytic properties of L(2) and L(10) M. burtonii Rubisco (MbR) were indistinguishable. At 25 degrees C they both shared a low specificity for CO(2) over O(2) (1.1 mol x mol(-1)) and RuBP carboxylation rates that were distinctively enhanced at low pH (approximately 4 s(-1) at pH 6, relative to 0.8 s(-1) at pH 8) with a temperature optimum of 55 degrees C. Like other archaeal Rubiscos, MbR also has a high O(2) affinity (K(m)(O(2)) = approximately 2.5 microM). The catalytic and structural similarities of MbR to other archaeal Rubiscos contrast with its closer sequence homology to bacterial L(2) Rubisco, complicating its classification within the Rubisco superfamily.

Putative functional role for the invariant aspartate 263 residue of Rhodospirillum rubrum Rubisco

Although aspartate residue D263 of Rhodospirillum rubrum Rubisco is close to the active site and invariant in all reported Rubiscos, its possible functional and structural roles in Rubisco activity have not been investigated. We have mutagenised D263 to several selected amino acids (asparagine, alanine, serine, glutamate, and glutamine) to probe possible roles in facilitating proton movements within the active site and maintaining structural positioning of key active-site groups. The mutants have been characterized by kinetic methods and by differential scanning calorimetry (DSC) to examine the effects of the substitutions on the stability of the folded state. We show that D263 is essential for maintaining effective levels of catalysis with the mutations reducing carboxylation variously by up to 100-fold but having less than 10% effect on the carboxylase/oxygenase specificity of the catalytic reaction. Removing the charge of the residue 263 side chain significantly strengthens binding of the activating (carbamylating) CO(2) molecule. In contrast, a charge on the 263 site has only a small influence on binding of the positively charged Mg(2+) ion, suggesting that the local protein structure provides different shielding of the formal charges on the Mg(2+) ion and the epsilon-lysine group of K191. Interestingly, introduction of an internal cavity (D263S and D263A) and insertion of an extra -CH(2)- group (D263E and D263Q) have opposite effects on catalysis, the former relatively small and the latter much larger, suggesting that the extra side-chain group induces a specific structural distortion that inhibits formation of the transition state. As the DSC results show that the mutations only slightly increase the kinetic stability of the folded state, we conclude that the rate-limiting (activated) step of unfolding involves substantial unfolding of the structure but not in the region of site 263. In summary, interaction of D263 with H287 of a largely electrostatic nature appears critical for maintaining correct positioning of catalytic groups in the active site. The conservation of D263 can thus be accounted for by its contribution to the maintenance of a finely tuned structure in this region abutting the active site.