Stability-activity tradeoffs constrain the adaptive evolution of RubisCO

Significance How enzymes acquire new functions is a key question in evolutionary biology. Here, we studied the evolution of some forms of ribulose-1,5-bisphosphate carboxylase, the enzyme responsible for CO2 fixation in photosynthesis, which has evolved enhanced activity in multiple groups of plants. We showed that the evolution of this enzyme was constrained by tradeoffs between activity and stability, two key properties of enzymes. The acquisition of enhanced activity was mediated by mutations destabilizing the structure. However, these were preceded and followed by periods in which stabilizing mutations were predominant, so that global stability was always maintained. This work shows that the natural evolution of enzymes is subject to strong biophysical constraints, and evolution follows perilous paths toward adaptation. A well-known case of evolutionary adaptation is that of ribulose-1,5-bisphosphate carboxylase (RubisCO), the enzyme responsible for fixation of CO2 during photosynthesis. Although the majority of plants use the ancestral C3 photosynthetic pathway, many flowering plants have evolved a derived pathway named C4 photosynthesis. The latter concentrates CO2, and C4 RubisCOs consequently have lower specificity for, and faster turnover of, CO2. The C4 forms result from convergent evolution in multiple clades, with substitutions at a small number of sites under positive selection. To understand the physical constraints on these evolutionary changes, we reconstructed in silico ancestral sequences and 3D structures of RubisCO from a large group of related C3 and C4 species. We were able to precisely track their past evolutionary trajectories, identify mutations on each branch of the phylogeny, and evaluate their stability effect. We show that RubisCO evolution has been constrained by stability-activity tradeoffs similar in character to those previously identified in laboratory-based experiments. The C4 properties require a subset of several ancestral destabilizing mutations, which from their location in the structure are inferred to mainly be involved in enhancing conformational flexibility of the open-closed transition in the catalytic cycle. These mutations are near, but not in, the active site or at intersubunit interfaces. The C3 to C4 transition is preceded by a sustained period in which stability of the enzyme is increased, creating the capacity to accept the functionally necessary destabilizing mutations, and is immediately followed by compensatory mutations that restore global stability.

[1]  G. Lorimer,et al.  Plant Photorespiration—An Inevitable Consequence of the Existence of Atmospheric Oxygen , 1973, Nature.

[2]  R. Chollet,et al.  Structural characterization and the determination of negative cooperativity in the tight binding of 2-carboxyarabinitol bisphosphate to higher plant ribulose bisphosphate carboxylase. , 1985, The Journal of biological chemistry.

[3]  M S Chapman,et al.  Tertiary structure of plant RuBisCO: domains and their contacts. , 1988, Science.

[4]  R. J. Spreitzer,et al.  Chloroplast intragenic suppression enhances the low CO2/O2 specificity of mutant ribulose-bisphosphate carboxylase/oxygenase. , 1989, The Journal of biological chemistry.

[5]  R. Jensen,et al.  Status of the substrate binding sites of ribulose bisphosphate carboxylase as determined with 2-C-carboxyarabinitol 1,5-bisphosphate. , 1990, Plant physiology.

[6]  M. Badger,et al.  Comparisons of rbcL genes for the large subunit of ribulose-bisphosphate carboxylase from closely related C3 and C4 plant species. , 1990, The Journal of biological chemistry.

[7]  T. Blundell,et al.  Comparative protein modelling by satisfaction of spatial restraints. , 1993, Journal of molecular biology.

[8]  R. J. Spreitzer,et al.  Directed Mutagenesis of Chloroplast Ribulose-1,5-bisphosphate Carboxylase/Oxygenase , 1996, The Journal of Biological Chemistry.

[9]  R. J. Spreitzer,et al.  Complementing Substitutions at the Bottom of the Barrel Influence Catalysis and Stability of Ribulose-bisphosphate Carboxylase/Oxygenase* , 1997, The Journal of Biological Chemistry.

[10]  S. Whelan,et al.  A general empirical model of protein evolution derived from multiple protein families using a maximum-likelihood approach. , 2001, Molecular biology and evolution.

[11]  D. M. Taverna,et al.  Why are proteins marginally stable? , 2002, Proteins.

[12]  Brian K Shoichet,et al.  Evolution of an antibiotic resistance enzyme constrained by stability and activity trade-offs. , 2002, Journal of molecular biology.

[13]  Beth M Beadle,et al.  Structural bases of stability-function tradeoffs in enzymes. , 2002, Journal of molecular biology.

[14]  R. Sage Variation in the k(cat) of Rubisco in C(3) and C(4) plants and some implications for photosynthetic performance at high and low temperature. , 2002, Journal of experimental botany.

[15]  R. Sage Variation in the kcat of Rubisco in C3 and C4 plants and some implications for photosynthetic performance at high and low temperature , 2002 .

[16]  R. J. Spreitzer,et al.  Assessment of Structural and Functional Divergence Far from the Large Subunit Active Site of Ribulose-1,5-bisphosphate Carboxylase/Oxygenase* , 2003, Journal of Biological Chemistry.

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

[18]  I. Andersson,et al.  Structural framework for catalysis and regulation in ribulose-1,5-bisphosphate carboxylase/oxygenase. , 2003, Archives of biochemistry and biophysics.

[19]  Conrad C. Huang,et al.  UCSF Chimera—A visualization system for exploratory research and analysis , 2004, J. Comput. Chem..

[20]  M. DePristo,et al.  Missense meanderings in sequence space: a biophysical view of protein evolution , 2005, Nature Reviews Genetics.

[21]  François Stricher,et al.  The FoldX web server: an online force field , 2005, Nucleic Acids Res..

[22]  Piet Herdewijn,et al.  Investigation of the DNA-dependent cyclohexenyl nucleic acid polymerization and the cyclohexenyl nucleic acid-dependent DNA polymerization , 2005, Nucleic acids research.

[23]  L. Serrano,et al.  Prediction of water and metal binding sites and their affinities by using the Fold-X force field. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[24]  Sriram Satagopan,et al.  Phylogenetic engineering at an interface between large and small subunits imparts land-plant kinetic properties to algal Rubisco. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[25]  F. Arnold,et al.  Protein stability promotes evolvability. , 2006, Proceedings of the National Academy of Sciences of the United States of America.

[26]  M. Kapralov,et al.  Widespread positive selection in the photosynthetic Rubisco enzyme , 2007, BMC Evolutionary Biology.

[27]  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.

[28]  Shoshana J. Wodak,et al.  Relating destabilizing regions to known functional sites in proteins , 2007, BMC Bioinformatics.

[29]  K. Henrick,et al.  Inference of macromolecular assemblies from crystalline state. , 2007, Journal of molecular biology.

[30]  Ziheng Yang PAML 4: phylogenetic analysis by maximum likelihood. , 2007, Molecular biology and evolution.

[31]  I. Andersson,et al.  Structure and function of Rubisco. , 2008, Plant physiology and biochemistry : PPB.

[32]  Flavien Russier,et al.  Evolutionary switch and genetic convergence on rbcL following the evolution of C4 photosynthesis. , 2008, Molecular biology and evolution.

[33]  François Stricher,et al.  How Protein Stability and New Functions Trade Off , 2008, PLoS Comput. Biol..

[34]  Geoffrey J. Barton,et al.  Jalview Version 2—a multiple sequence alignment editor and analysis workbench , 2009, Bioinform..

[35]  Yasuro Kadono,et al.  Molecular Adaptation of rbcL in the Heterophyllous Aquatic Plant Potamogeton , 2009, PloS one.

[36]  Richard A. Goldstein,et al.  Identifying Changes in Selective Constraints: Host Shifts in Influenza , 2009, PLoS Comput. Biol..

[37]  Philip A. Romero,et al.  Exploring protein fitness landscapes by directed evolution , 2009, Nature Reviews Molecular Cell Biology.

[38]  John C. Avise,et al.  In the Light of Directed Evolution: Pathways of Adaptive Protein Evolution , 2009 .

[39]  G. Schreiber,et al.  Assessing computational methods for predicting protein stability upon mutation: good on average but not in the details. , 2009, Protein engineering, design & selection : PEDS.

[40]  Bartek Wilczynski,et al.  Biopython: freely available Python tools for computational molecular biology and bioinformatics , 2009, Bioinform..

[41]  Dan S. Tawfik,et al.  Stability effects of mutations and protein evolvability. , 2009, Current opinion in structural biology.

[42]  R. Freckleton,et al.  Can phylogenetics identify C(4) origins and reversals? , 2010, Trends in ecology & evolution.

[43]  B. Shoichet,et al.  Structural bases for stability-function tradeoffs in antibiotic resistance. , 2010, Journal of molecular biology.

[44]  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.

[45]  Mauno Vihinen,et al.  Performance of protein stability predictors , 2010, Human mutation.

[46]  M. Harms,et al.  Analyzing protein structure and function using ancestral gene reconstruction. , 2010, Current opinion in structural biology.

[47]  Dan S. Tawfik,et al.  Mutational effects and the evolution of new protein functions , 2010, Nature Reviews Genetics.

[48]  D. van der Spoel,et al.  Subunit interface dynamics in hexadecameric rubisco. , 2011, Journal of molecular biology.

[49]  Bo Wang,et al.  Molecular evolution of rbcL in three gymnosperm families: identifying adaptive and coevolutionary patterns , 2011, Biology Direct.

[50]  Thomas A. Hopf,et al.  Protein 3D Structure Computed from Evolutionary Sequence Variation , 2011, PloS one.

[51]  I. Andersson,et al.  Changes in Rubisco kinetics during the evolution of C4 photosynthesis in Flaveria (Asteraceae) are associated with positive selection on genes encoding the enzyme. , 2011, Molecular biology and evolution.

[52]  R. Sage,et al.  The C(4) plant lineages of planet Earth. , 2011, Journal of experimental botany.

[53]  J. Galmés,et al.  Isoleucine 309 acts as a C4 catalytic switch that increases ribulose-1,5-bisphosphate carboxylase/oxygenase (rubisco) carboxylation rate in Flaveria , 2011, Proceedings of the National Academy of Sciences.

[54]  M. Anisimova,et al.  Coevolution of amino acid residues in the key photosynthetic enzyme Rubisco , 2011, BMC Evolutionary Biology.

[55]  Y. Kai,et al.  Crystal structure of rice Rubisco and implications for activation induced by positive effectors NADPH and 6-phosphogluconate. , 2012, Journal of molecular biology.

[56]  J. A. Smith,et al.  Rubisco Evolution in C4 Eudicots: An Analysis of Amaranthaceae Sensu Lato , 2012, PloS one.

[57]  T. Hodkinson,et al.  New grass phylogeny resolves deep evolutionary relationships and discovers C4 origins. , 2012, The New phytologist.

[58]  R. Rickaby,et al.  Adaptive signals in algal Rubisco reveal a history of ancient atmospheric carbon dioxide , 2012, Philosophical Transactions of the Royal Society B: Biological Sciences.

[59]  R. Sage,et al.  Photorespiration and the evolution of C4 photosynthesis. , 2012, Annual review of plant biology.

[60]  Jianbo Wang,et al.  Adaptive evolution of the rbcL gene in Brassicaceae , 2012 .

[61]  Songnian Hu,et al.  EvolView, an online tool for visualizing, annotating and managing phylogenetic trees , 2012, Nucleic Acids Res..

[62]  R Core Team,et al.  R: A language and environment for statistical computing. , 2014 .