Ribulose 1,5-bisphosphate (RuBP) carboxylase-oxygenase (Rubisco) plays a fundamental role in the carbon cycle by fixing the atmospheric CO2 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.

In the carboxylation reaction catalyzed by ribulose 1,5-bisphosphate (RuBP) carboxylase-oxygenase (Rubisco), which is fundamental to photosynthesis, scission of a C-C bond in the six-carbon gemdiolate intermediate forms a carbanion that must be protonated stereospecifically to form product. It is thought that a conserved lysine side chain (LYS175 in spinach Rubisco), in the immediate vicinity of the carbanion, provides the necessary proton. Here, we endeavor to determine from the electronic-structure calculations whether protonation via this route is energetically possible. The two-dimensional energy surface was mapped to determine the minimum energy path (MEP) using density functional theory (B3LYP) and incorporating basis set superposition and classical (London) dispersion corrections. The potential of mean force (free energy) was then calculated from ab initio molecular dynamics simulations with umbrella sampling in the vicinity of the MEP on the scission-protonation reaction coordinate. MEP calculations were also carried out to evaluate the possibility of an active-site water near the phosphate (P1) of RuBP, with an excess proton positioned at P1, as an alternative facilitator of stereospecific protonation via a classical Grotthuss mechanism. In both cases, the C-C bond scission in the six-carbon intermediate and proton transfer from the donor was found to be concerted and highly asynchronous, without a stable carbanion intermediate. However, the free energy change was unfavorable for direct protonation by the LYS175 side chain. In contrast, the Grotthuss mechanism yielded stable products and an activation energy in good agreement with experiment. It also provides a plausible mechanism for alternative product formed in enzyme mutations at the LYS175 position and is consistent with the observed deuterium isotope effects.

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

Ribulose 1,5-bisphosphate carboxylase-oxygenase (Rubisco) is the primary carbon-fixing enzyme in photosynthesis, fixing CO2 to a 5-carbon sugar, RuBP, in a series of five reactions. However, it also catalyzes an oxygenase reaction by O2 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 CO2 or O2 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 (3O2), its reaction with near-universal singlet-state molecules is formally spin forbidden. Although in oxygenase enzymes, 3O2 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 3O2 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 Mg2+ and its ligands extends the region for charge delocalization and stabilization of its negative charge to allow formation of a caged biradical enediolate-O2 complex. Thus, Rubisco is a unique type of oxygenase without precedent in the literature. Second, for the O2 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 3O2. This finding supports our recent proposal that deoxygenation, that is, reversal of gas binding, is possible. Third, as neither CO2 nor O2 binds to the enzyme, our findings support the proposal from our recent carboxylase study that the observed KC or KO (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 3O2 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.

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 2 . Hydration of the RuBP.CO 2 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 2 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 2 ; both hydration and CC bond scission require early protonation of CO 2 in the RuBP.CO 2 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.

Based on our critique of requirements for performing an efficient molecular dynamics simulation with the particle-mesh Ewald (PME) implementation in GROMACS 4.5, we present a computational tool to enable the discovery of parameters that produce a given accuracy in the PME approximation of the full electrostatics. Calculations on two parallel computers with different processor and communication structures showed that a given accuracy can be attained over a range of parameter space, and that the attributes of the hardware and simulation system control which parameter sets are optimal. This information can be used to find the fastest available PME parameter sets that achieve a given accuracy. We hope that this tool will stimulate future work to assess the impact of the quality of the PME approximation on simulation outcomes, particularly with regard to the trade-off between cost and scientific reliability in biomolecular applications. (C) 2011 Wiley Periodicals, Inc. J Comput Chem 32:2031-2040, 2011

Disulfide torsional energy, a good predictor of disulfide redox potential in proteins, may be estimated by interpolation on a potential energy surface (PES) describing the twisting of diethyl disulfide through its three central dihedral angles. Here we update PES calculations at the M05-2X level of theory with the 6-31G(d) basis set. Although the surface shows no qualitative differences from an earlier MP2(full) PES, energy differences greater than 1 kJ mol–1 were seen for conformations with χ2 between –60° and 30°, or with χ3 below 60° or above 130°. This is particularly significant for highly strained disulfides that are likely to be spontaneously reduced by mechanical means. In benchmarking against the high-level G3X method, M05-2X showed significantly reduced root mean squared deviation compared with MP2(full) (1.0 versus 2.0 kJ mol–1 respectively). Results are incorporated into a web application that calculates relative torsional energies from disulfide dihedral angles (http://www.sbinf.org/applications/pes.html).

Methyltetrahydrofolate:corrinoid/iron-sulfur protein methyltransferase (MeTr) catalyzes the transfer of the N5-methyl group from N5-methyltetrahydrofolate (CH(3)THF) to the cobalt center of a corrinoid/iron-sulfur protein, a reaction similar to that of cobalamin-dependent methionine synthase (MetH). For such a reaction to occur, CH(3)THF is expected to be activated by a stereospecific protonation at the N5 position. It has been shown experimentally that binding to MeTr is associated with a pK(a) increase and proton uptake. The enzyme could achieve this by binding the unprotonated form of CH(3)THF, followed by specific protonation at the correct orientation. Here we have used computational approaches to investigate the protonation state of the ligand and active-site residues in MeTr. First, quantum mechanical (QM) methods with the PCM solvation model were used to predict protonation positions and pK(a) values of pterin, folate, and their analogues in an aqueous environment. After a reliable calibration of computational and experimental results was obtained, the effect of the protein environment was then considered. Different protonation states of CH(3)THF and active-site aspartic residues (D75 and D160) were investigated using QM calculations of active-site fragment complexes and the perturbed quantum atom (PQA) approach within QM/MM simulations. The final free energy results indicate that the N5 position of the tetrahydropterin ring is the preferred protonation position of CH(3)THF when bound to the active site of MeTr, followed by Asp160. We also found that the active-site environment is likely to increase the pK(a) of N5 by about 3 units, leading to proton uptake upon CH(3)THF binding, as observed experimentally for MeTr. Some implications of the results are discussed for the MetH mechanism.

Picoeukaryotes are a taxonomically diverse group of organisms less than 2 micrometers in diameter. Photosynthetic marine picoeukaryotes in the genus Micromonas thrive in ecosystems ranging from tropical to polar and could serve as sentinel organisms for biogeochemical fluxes of modern oceans during climate change. These broadly distributed primary producers belong to an anciently diverged sister clade to land plants. Although Micromonas isolates have high 18S ribosomal RNA gene identity, we found that genomes from two isolates shared only 90% of their predicted genes. Their independent evolutionary paths were emphasized by distinct riboswitch arrangements as well as the discovery of intronic repeat elements in one isolate, and in metagenomic data, but not in other genomes. Divergence appears to have been facilitated by selection and acquisition processes that actively shape the repertoire of genes that are mutually exclusive between the two isolates differently than the core genes. Analyses of the Micromonas genomes offer valuable insights into ecological differentiation and the dynamic nature of early plant evolution.

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.

We address the question of constructing a protocol for replica-exchange molecular dynamics (REMD) simulations that make efficient use of the replica space, assess whether published applications are achieving such "mixing" efficiency, and provide a how-to guide to assist users to plan efficient REMD simulations. To address our first question, we introduce and discuss three metrics for assessing the number of replica-exchange attempts required to justify the use of a replica scheme and define a "transit number" as the lower bound for the length of an efficient simulation. Our literature survey of applications of REMD simulations of peptides in explicit solvent indicated that authors are not routinely reporting sufficient details of their simulation protocols to allow readers to make independent assessments of the impact of the method on their results, particularly whether mixing efficiency has been achieved. Necessary details include the expected or observed replica-exchange probability, together with the total number of exchange attempts, the exchange period, and estimates of the autocorrelation time of the potential energy. Our analysis of cases where the necessary information was reported suggests that in many of these simulations there are insufficient exchanges attempted or an insufficiently long period between them to provide confidence that the simulation length justifies the size of the replica scheme. We suggest guidelines for designing REMD simulation protocols to ensure mixing efficiency. Two key recommendations are that the exchange period should in general be larger than 1 ps and the number of exchange attempts should be chosen to significantly exceed the transit number for the replica scheme.

Using comparative genomics and in-silico analyses, we previously identified a new member of the prion-protein (PrP) family, the gene SPRN, encoding the protein Shadoo (Sho), and suggested its functions might overlap with those of PrP. Extended bioinformatics and conceptual biology studies to elucidate Sho's functions now reveal Sho has a conserved RGG-box motif, a well-known RNA-binding motif characterized in proteins such as FragileX Mental Retardation Protein. We report a systematic comparative analysis of RGG-box containing proteins which highlights the motif's functional versatility and supports the suggestion that Sho plays a dual role in cell signaling and RNA binding in brain. These findings provide a further link to PrP, which has well-characterized RNA-binding properties.

Crystallographic, mutagenesis, kinetic, and computational studies on Rubisco over three decades have revealed much about its catalytic mechanism and the role played by several active-site residues. However, key questions remain unanswered. Specific details of the carboxylase and oxygenase mechanisms, required to underpin the rational re-engineering of Rubisco, are still speculative. Here we address critical gaps in knowledge with a definitive comprehensive computational investigation of the mechanism of carboxylase activity at the Rubisco active site. Density functional theory calculations (B3LYP/6-31G(d,p)) were performed on active-site fragment models of a size up to 77 atoms, not previously possible computationally. All amino acid residues suspected to play roles in the acid-base chemistry in the multistep reaction, and interacting directly with the central Mg2+ atom and the reactive moiety of substrate and intermediates, were included. The results provide a firm basis for us to propose a novel mechanism for the entire sequence of reactions in the carboxylase catalysis and to define precise roles for the activesite residues, singly and in concert. In this mechanism, the carbamylated LYS201 plays a more limited role than previously proposed but is crucial for initiating the reaction by acting as a base in the enolization. We suggest a wider role for HIS294, with involvement in the carboxylation, hydration, and C2-C3 bond-scission steps, consistent with the suggestion of Harpel et al. (1998) but contrary to the consensus view of Cleland et al. (1998). In contrast to the common assumption that the water molecule for the hydration step comes from within the active site, we propose that the Mg-coordinated water is not dissociated at the start of the gas-addition reaction but rather remains coordinated and is used for the hydration of the C3 carbon atom. New roles are also proposed for LYS175, GLU204, and HIS294. The mechanism suggests roles in the gas-addition step for residues in three spatially distinct regions of the active site, HIS294 and LYS334 in the C-terminal domain of the large subunit (LSU), but also hitherto unsuspected roles for a cluster of three residues (ASN123, GLU60, and TYR20) in the N-terminal domain of the partner LSU of the dinner containing the active site. Our new mechanism is supported by existing experimental data, provides new convincing interpretations of previously puzzling data, and allows new insights into mutational strategies for improving Rubisco activity.