Search

The history of Earth’s carbon cycle reflects trends in atmospheric composition convolved with the evolution of photosynthesis. Fortunately, key parts of the carbon cycle have been recorded in the carbon isotope ratios of sedimentary rocks. The dominant model used to interpret this record as a proxy for ancient atmospheric CO 2 is based on carbon isotope fractionations of modern photoautotrophs, and longstanding questions remain about how their evolution might have impacted the record. Therefore, we measured both biomass (ε p ) and enzymatic (ε Rubisco ) carbon isotope fractionations of a cyanobacterial strain ( Synechococcus elongatus PCC 7942) solely expressing a putative ancestral Form 1B rubisco dating to ≫1 Ga. This strain, nicknamed ANC, grows in ambient pCO 2 and displays larger ε p values than WT, despite having a much smaller ε Rubisco (17.23 ± 0.61‰ vs. 25.18 ± 0.31‰, respectively). Surprisingly, ANC ε p exceeded ANC ε Rubisco in all conditions tested, contradicting prevailing models of cyanobacterial carbon isotope fractionation. Such models can be rectified by introducing additional isotopic fractionation associated with powered inorganic carbon uptake mechanisms present in Cyanobacteria, but this amendment hinders the ability to accurately estimate historical pCO 2 from geological data. Understanding the evolution of rubisco and the CO 2 concentrating mechanism is therefore critical for interpreting the carbon isotope record, and fluctuations in the record may reflect the evolving efficiency of carbon fixing metabolisms in addition to changes in atmospheric CO 2 .

Cyanobacteria rely on CO 2 -concentrating mechanisms (CCMs) to grow in today’s atmosphere (0.04% CO 2 ). These complex physiological adaptations require ≈15 genes to produce two types of protein complexes: inorganic carbon (Ci) transporters and 100+ nm carboxysome compartments that encapsulate rubisco with a carbonic anhydrase (CA) enzyme. Mutations disrupting any of these genes prohibit growth in ambient air. If any plausible ancestral form—i.e., lacking a single gene—cannot grow, how did the CCM evolve? Here, we test the hypothesis that evolution of the bacterial CCM was “catalyzed” by historically high CO 2 levels that decreased over geologic time. Using an E. coli reconstitution of a bacterial CCM, we constructed strains lacking one or more CCM components and evaluated their growth across CO 2 concentrations. We expected these experiments to demonstrate the importance of the carboxysome. Instead, we found that partial CCMs expressing CA or Ci uptake genes grew better than controls in intermediate CO 2 levels (≈1%) and observed similar phenotypes in two autotrophic bacteria, Halothiobacillus neapolitanus and Cupriavidus necator . To understand how CA and Ci uptake improve growth, we model autotrophy as colimited by CO 2 and HCO 3 − , as both are required to produce biomass. Our experiments and model delineated a viable trajectory for CCM evolution where decreasing atmospheric CO 2 induces an HCO 3 − deficiency that is alleviated by acquisition of CA or Ci uptake, thereby enabling the emergence of a modern CCM. This work underscores the importance of considering physiology and environmental context when studying the evolution of biological complexity.

From a metabolic perspective, molecular oxygen (O

Enzymes are well known for their catalytic abilities, some even reaching "catalytic perfection" in the sense that the reaction they catalyze has reached the physical bound of the diffusion rate. However, our growing understanding of enzyme superfamilies has revealed that only some share a catalytic chemistry while others share a substrate-handle binding motif, for example, for a particular phosphate group. This suggests that some families emerged through a "substrate-handle-binding-first" mechanism ("binding-first" for brevity) instead of "chemistry-first" and we are, therefore, left to wonder what the role of non-catalytic binders might have been during enzyme evolution. In the last of their eight seminal, back-to-back articles from 1976, John Albery and Jeremy Knowles addressed the question of enzyme evolution by arguing that the simplest mode of enzyme evolution is what they defined as "uniform binding" (parallel stabilization of all enzyme-bound states to the same degree). Indeed, we show that a uniform-binding proto-catalyst can accelerate a reaction, but only when catalysis is already present, that is, when the transition state is already stabilized to some degree. Thus, we sought an alternative explanation for the cases where substrate-handle-binding preceded any involvement of a catalyst. We find that evolutionary starting points that exhibit negative catalysis can redirect the reaction's course to a preferred product without need for rate acceleration or product release; that is, if they do not stabilize, or even destabilize, the transition state corresponding to an undesired product. Such a mechanism might explain the emergence of "binding-first" enzyme families like the aldolase superfamily.

Significance The Calvin cycle is the most important carbon fixation pathway in the biosphere. However, its carboxylating enzyme Rubisco also accepts oxygen, thus producing 2-phosphoglycolate. Phosphoglycolate salvage pathways were extensively studied in photoautotrophs but remain uncharacterized in chemolithoautotrophs using the Calvin cycle. Here, we study phosphoglycolate salvage in the chemolithoautotrophic model bacterium Cupriavidus necator H16. We demonstrate that this bacterium mainly reassimilates 2-phosphoglycolate via the glycerate pathway. Upon disruption of this pathway, a secondary route, which we term the malate cycle, supports photorespiration by completely oxidizing 2-phosphoglycolate to CO2. While the malate cycle was not previously known to metabolize 2-phosphoglycolate in nature, a bioinformatic analysis suggests that it may support phosphoglycolate salvage in diverse chemoautotrophic bacteria., Carbon fixation via the Calvin cycle is constrained by the side activity of Rubisco with dioxygen, generating 2-phosphoglycolate. The metabolic recycling of phosphoglycolate was extensively studied in photoautotrophic organisms, including plants, algae, and cyanobacteria, where it is referred to as photorespiration. While receiving little attention so far, aerobic chemolithoautotrophic bacteria that operate the Calvin cycle independent of light must also recycle phosphoglycolate. As the term photorespiration is inappropriate for describing phosphoglycolate recycling in these nonphotosynthetic autotrophs, we suggest the more general term “phosphoglycolate salvage.” Here, we study phosphoglycolate salvage in the model chemolithoautotroph Cupriavidus necator H16 (Ralstonia eutropha H16) by characterizing the proxy process of glycolate metabolism, performing comparative transcriptomics of autotrophic growth under low and high CO2 concentrations, and testing autotrophic growth phenotypes of gene deletion strains at ambient CO2. We find that the canonical plant-like C2 cycle does not operate in this bacterium, and instead, the bacterial-like glycerate pathway is the main route for phosphoglycolate salvage. Upon disruption of the glycerate pathway, we find that an oxidative pathway, which we term the malate cycle, supports phosphoglycolate salvage. In this cycle, glyoxylate is condensed with acetyl coenzyme A (acetyl-CoA) to give malate, which undergoes two oxidative decarboxylation steps to regenerate acetyl-CoA. When both pathways are disrupted, autotrophic growth is abolished at ambient CO2. We present bioinformatic data suggesting that the malate cycle may support phosphoglycolate salvage in diverse chemolithoautotrophic bacteria. This study thus demonstrates a so far unknown phosphoglycolate salvage pathway, highlighting important diversity in microbial carbon fixation metabolism.

Flamholz and Shih explain how photosynthetic organisms on earth have evolved carbon dioxide concentrating mechanisms to contend with an increased abundance of oxygen in the atmosphere. This critical innovation has occurred numerous times and comes in two basic flavors.

Making sense of the metabolism of microbial communities is a daunting task. Using denitrification as a model metabolism, a new paper shows that the rate of denitrification can often be predicted from genome contents, and dynamical models can be composed to predict denitrification rates of communities of two to five species.

SummaryThe Human Impacts Database (www.anthroponumbers.org) is a curated searchable resource housing quantitative data relating to the diverse anthropogenic impacts on our planet, with topics ranging from sea level rise, to livestock populations, greenhouse gas emissions, fertilizer use, and beyond. Each entry in the database relates a quantitative value (or a time-series of values) along with clear referencing of the primary source, the method of measurement or estimation, an assessment of uncertainty, links to the underlying data, as well as a permanent identifier called an Human Impacts ID (“HuID”). While there are other databases that house some of these values, they are typically focused on a single topic area like energy usage or greenhouse gas emissions. The Human Impacts Database provides centralized access to quantitative information about the myriad ways in which humans impact the Earth, giving links to more specialized databases for interested readers. Here, we outline the structure of the database and describe our curation procedures. Finally, we use this database to generate a graphical summary of the current state of human impacts on the Earth, illustrating both their numerical values and their dense interconnections.The Bigger PictureOver the last 10,000 years, human activities have transformed the Earth through farming, forestry, mining and industry. The complex results of these activities are now observed and quantified as “human impacts” on Earth’s atmosphere, oceans, biosphere and geochemistry. While myriad studies have explored facets of human impacts on the planet, they are necessarily technical and often tightly-focused. Thus, finding reliable quantitative information requires a significant investment of time to assess each quantity, its methods of determination, and associated uncertainty. We present the Human Impacts Database (www.anthroponumbers.org), which houses a diverse array of such quantities. We review a subset of these values and how they help build intuition for understanding the Earth-human system. While collation alone does not tell us how to best ameliorate human impacts, we contend that any future plans should be made in light of a quantitative understanding of the interconnected ways in which humans impact the planet.

eQuilibrator (equilibrator.weizmann.ac.il) is a database of biochemical equilibrium constants and Gibbs free energies, originally designed as a web-based interface. While the website now counts around 1,000 distinct monthly users, its design could not accommodate larger compound databases and it lacked a scalable Application Programming Interface (API) for integration into other tools developed by the systems biology community. Here, we report on the recent updates to the database as well as the addition of a new Python-based interface to eQuilibrator that adds many new features such as a 100-fold larger compound database, the ability to add novel compounds, improvements in speed and memory use, and correction for Mg2+ ion concentrations. Moreover, the new interface can compute the covariance matrix of the uncertainty between estimates, for which we show the advantages and describe the application in metabolic modelling. We foresee that these improvements will make thermodynamic modelling more accessible and facilitate the integration of eQuilibrator into other software platforms., Nucleic Acids Research, 50 (D1), ISSN:1362-4962, ISSN:0301-5610

Wastewater surveillance for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) RNA can be integrated with COVID-19 case data to inform timely pandemic response. However, more research is needed to apply and develop systematic methods to interpret the true SARS-CoV-2 signal from noise introduced in wastewater samples (e.g., from sewer conditions, sampling and extraction methods, etc.). In this study, raw wastewater was collected weekly from five sewersheds and one residential facility. The concentrations of SARS-CoV-2 in wastewater samples were compared to geocoded COVID-19 clinical testing data. SARS-CoV-2 was reliably detected (95% positivity) in frozen wastewater samples when reported daily new COVID-19 cases were 2.4 or more per 100,000 people. To adjust for variation in sample fecal content, four normalization biomarkers were evaluated: crAssphage, pepper mild mottle virus, Bacteroides ribosomal RNA (rRNA), and human 18S rRNA. Of these, crAssphage displayed the least spatial and temporal variability. Both unnormalized SARS-CoV-2 RNA signal and signal normalized to crAssphage had positive and significant correlation with clinical testing data (Kendall's Tau-b (τ)=0.43 and 0.38, respectively), but no normalization biomarker strengthened the correlation with clinical testing data. Locational dependencies and the date associated with testing data impacted the lead time of wastewater for clinical trends, and no lead time was observed when the sample collection date (versus the result date) was used for both wastewater and clinical testing data. This study supports that trends in wastewater surveillance data reflect trends in COVID-19 disease occurrence and presents tools that could be applied to make wastewater signal more interpretable and comparable across studies., Graphical abstract Image, graphical abstract

Coenzymes distribute a variety of chemical moieties throughout cellular metabolism, participating in group (e.g., phosphate, acyl) and electron transfer. For a variety of reactions requiring acceptors or donors of specific resources, there often exist degenerate sets of molecules (e.g., NAD(H) and NADP(H)) that carry out similar functions. Although the physiological roles of various coenzyme systems are well established, it is unclear what selective pressures may have driven the emergence of coenzyme redundancy. Here we use genome-wide metabolic modeling approaches to decompose the selective pressures driving enzymatic specificity for either NAD(H) or NADP(H) in the metabolic network ofEscherichia coli. We found that few enzymes are thermodynamically constrained to using a single coenzyme, and in principle, a metabolic network relying on only NAD(H) is feasible. However, structural and sequence analyses revealed widespread conservation of residues that retain selectivity for either NAD(H) or NADP(H), suggesting that additional forces may shape specificity. Using a model accounting for the cost of oxidoreductase enzyme expression, we found that coenzyme redundancy universally reduces the minimal amount of protein required to catalyze coenzyme-coupled reactions, inducing individual reactions to strongly prefer one coenzyme over another when reactions are near thermodynamic equilibrium. We propose that protein minimization generically promotes coenzyme redundancy, and that coenzymes typically thought to exist in a single pool (e.g., CoA) may exist in more than one form (e.g., dephospho-CoA).Significance statementMetabolism relies on a small class of molecules (coenzymes) that serve as universal donors and acceptors of key chemical groups and electrons. Although metabolic networks crucially depend on structurally redundant coenzymes (e.g., NAD(H) and NADP(H)) associated with different enzymes, the criteria that led to the emergence of this redundancy remain poorly understood. Our combination of modeling, and structural and sequence analysis indicates that coenzyme redundancy is not essential for metabolism, but rather an evolved strategy promoting efficient usage of enzymes when biochemical reactions are near equilibrium. Our work suggests that early metabolism may have operated with fewer coenzymes, and that adaptation for metabolic efficiency may have driven the rise of coenzyme diversity in living systems.

Wastewater surveillance for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) RNA can be integrated with COVID-19 case data to inform timely pandemic response. However, more research is needed to apply and develop systematic methods to interpret the true SARS-CoV-2 signal from noise introduced in wastewater samples (e.g., from sewer conditions, sampling and extraction methods, etc.). In this study, raw wastewater was collected weekly from five sewersheds and one residential facility, and wastewater SARS-CoV-2 concentrations were compared to geocoded COVID-19 clinical testing data. SARS-CoV-2 was reliably detected (95% positivity) in frozen wastewater samples when reported daily new COVID-19 cases were 2.4 or more per 100,000 people. To adjust for variation in sample fecal content, crAssphage, pepper mild mottle virus, Bacteroides ribosomal RNA (rRNA), and human 18S rRNA were evaluated as normalization biomarkers, and crAssphage displayed the least spatial and temporal variability. Both unnormalized SARS-CoV-2 RNA signal and signal normalized to crAssphage had positive and significant correlation with clinical testing data (Kendall’s Tau-b (τ)=0.43 and 0.38, respectively). Locational dependencies and the date associated with testing data impacted the lead time of wastewater for clinical trends, and no lead time was observed when the sample collection date (versus the result date) was used for both wastewater and clinical testing data. This study supports that trends in wastewater surveillance data reflect trends in COVID-19 disease occurrence and presents approaches that could be applied to make wastewater signal more interpretable and comparable across studies.

Viral genome sequencing has guided our understanding of the spread and extent of genetic diversity of SARS-CoV-2 during the COVID-19 pandemic. SARS-CoV-2 viral genomes are usually sequenced from nasopharyngeal swabs of individual patients to track viral spread., Viral genome sequencing has guided our understanding of the spread and extent of genetic diversity of SARS-CoV-2 during the COVID-19 pandemic. SARS-CoV-2 viral genomes are usually sequenced from nasopharyngeal swabs of individual patients to track viral spread. Recently, RT-qPCR of municipal wastewater has been used to quantify the abundance of SARS-CoV-2 in several regions globally. However, metatranscriptomic sequencing of wastewater can be used to profile the viral genetic diversity across infected communities. Here, we sequenced RNA directly from sewage collected by municipal utility districts in the San Francisco Bay Area to generate complete and nearly complete SARS-CoV-2 genomes. The major consensus SARS-CoV-2 genotypes detected in the sewage were identical to clinical genomes from the region. Using a pipeline for single nucleotide variant calling in a metagenomic context, we characterized minor SARS-CoV-2 alleles in the wastewater and detected viral genotypes which were also found within clinical genomes throughout California. Observed wastewater variants were more similar to local California patient-derived genotypes than they were to those from other regions within the United States or globally. Additional variants detected in wastewater have only been identified in genomes from patients sampled outside California, indicating that wastewater sequencing can provide evidence for recent introductions of viral lineages before they are detected by local clinical sequencing. These results demonstrate that epidemiological surveillance through wastewater sequencing can aid in tracking exact viral strains in an epidemic context.

Significance Knowing the absolute numbers of virions in an infection promotes better understanding of disease dynamics and response of the immune system. Here we use current knowledge on the concentrations of virions in infected individuals to estimate the total number and mass of SARS-CoV-2 virions in an infected person. Although each infected person carries an estimated 1 billion to 100 billion virions during peak infection, their total mass is no more than 0.1 mg. This curiously implies that all SARS-CoV-2 virions currently in all human hosts have a mass of between 100 g and 10 kg. Combining the known mutation rate and our estimate of the number of infectious virions, we quantify the formation rate of genetic variants., Quantitatively describing the time course of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection within an infected individual is important for understanding the current global pandemic and possible ways to combat it. Here we integrate the best current knowledge about the typical viral load of SARS-CoV-2 in bodily fluids and host tissues to estimate the total number and mass of SARS-CoV-2 virions in an infected person. We estimate that each infected person carries 109 to 1011 virions during peak infection, with a total mass in the range of 1 μg to 100 μg, which curiously implies that all SARS-CoV-2 virions currently circulating within human hosts have a collective mass of only 0.1 kg to 10 kg. We combine our estimates with the available literature on host immune response and viral mutation rates to demonstrate how antibodies markedly outnumber the spike proteins, and the genetic diversity of virions in an infected host covers all possible single nucleotide substitutions.

Many photosynthetic organisms employ a CO2concentrating mechanism (CCM) to increase the rate of CO2fixation via the Calvin cycle. CCMs catalyze ≈50% of global photosynthesis, yet it remains unclear which genes and proteins are required to produce this complex adaptation. We describe the construction of a functional CCM in a non-native host, achieved by expressing genes from an autotrophic bacterium in anEscherichia colistrain engineered to depend on rubisco carboxylation for growth. Expression of 20 CCM genes enabledE. colito grow by fixing CO2from ambient air into biomass, with growth in ambient air depending on the components of the CCM. Bacterial CCMs are therefore genetically compact and readily transplanted, rationalizing their presence in diverse bacteria. Reconstitution enabled genetic experiments refining our understanding of the CCM, thereby laying the groundwork for deeper study and engineering of the cell biology supporting CO2assimilation in diverse organisms.

Many photosynthetic organisms employ a CO2concentrating mechanism (CCM) to increase the rate of CO2fixation via the Calvin cycle. CCMs catalyze ≈50% of global photosynthesis, yet it remains unclear which genes and proteins are required to produce this complex adaptation. We describe the construction of a functional CCM in a non-native host, achieved by expressing genes from an autotrophic bacterium in an engineeredE. colistrain. Expression of 20 CCM genes enabledE. colito grow by fixing CO2from ambient air into biomass, with growth depending on CCM components. Bacterial CCMs are therefore genetically compact and readily transplanted, rationalizing their presence in diverse bacteria. Reconstitution enabled genetic experiments refining our understanding of the CCM, thereby laying the groundwork for deeper study and engineering of the cell biology supporting CO2assimilation in diverse organisms.One Sentence SummaryA bacterial CO2concentrating mechanism enablesE. colito fix CO2from ambient air.

Carbon fixation via the Calvin cycle is constrained by the side activity of Rubisco with dioxygen, generating 2-phosphoglycolate. The metabolic recycling of 2-phosphoglycolate, an essential process termed photorespiration, was extensively studied in photoautotrophic organisms, including plants, algae, and cyanobacteria, but remains uncharacterized in chemolithoautotrophic bacteria. Here, we study photorespiration in the model chemolithoautotroph Cupriavidus necator (Ralstonia eutropha) by characterizing the proxy-process of glycolate metabolism, performing comparative transcriptomics of autotrophic growth under low and high CO2 concentrations, and testing autotrophic growth phenotypes of gene deletion strains at ambient CO2. We find that the canonical plant-like C2 cycle does not operate in this bacterium and instead the bacterial-like glycerate pathway is the main photorespiratory pathway. Upon disruption of the glycerate pathway, we find that an oxidative pathway, which we term the malate cycle, supports photorespiration. In this cycle, glyoxylate is condensed with acetyl-CoA to give malate, which undergoes two oxidative decarboxylation steps to regenerate acetyl-CoA. When both pathways are disrupted, autotrophic growth is abolished at ambient CO2. We present bioinformatic data suggesting that the malate cycle may support photorespiration in diverse chemolithoautotrophic bacteria. This study thus demonstrates a so-far unknown photorespiration pathway, highlighting important diversity in microbial carbon fixation metabolism.

Bacterial autotrophs often rely on CO2 concentrating mechanisms (CCMs) to assimilate carbon. Although many CCM proteins have been identified, a systematic screen of the components of CCMs is lacking. Here, we performed a genome-wide barcoded transposon screen to identify essential and CCM-related genes in the γ-proteobacterium Halothiobacillus neapolitanus. Screening revealed that the CCM comprises at least 17 and probably no more than 25 genes, most of which are encoded in 3 operons. Two of these operons (DAB1 and DAB2) contain a two-gene locus that encodes a domain of unknown function (Pfam: PF10070) and a putative cation transporter (Pfam: PF00361). Physiological and biochemical assays demonstrated that these proteins-which we name DabA and DabB, for DABs accumulate bicarbonate-assemble into a heterodimeric complex, which contains a putative β-carbonic anhydrase-like active site and functions as an energy-coupled inorganic carbon (Ci) pump. Interestingly, DAB operons are found in a diverse range of bacteria and archaea. We demonstrate that functional DABs are present in the human pathogens Bacillus anthracis and Vibrio cholerae. On the basis of these results, we propose that DABs constitute a class of energized Ci pumps and play a critical role in the metabolism of Ci throughout prokaryotic phyla.

Significance The k cat values of enzymes are important for the study of metabolic systems. However, the current use of k cat presents major difficulties, as values for most enzymes have not been experimentally measured, and experimentally available values are often measured under nonphysiological conditions, thereby casting doubt on the relevance of k cat under in vivo conditions. We present an approach that utilizes omics data to quantitatively analyze the relationship between in vitro k cat values and the maximal catalytic rate of enzymes in vivo. Our approach offers a high-throughput method to obtain enzyme kinetic constants, which reflect in vivo conditions, and are useful for more accurate and complete cellular metabolic models.

Rubisco is the primary carboxylase of the Calvin cycle, the most abundant enzyme in the biosphere, and one of the best-characterized enzymes. On the basis of correlations between Rubisco kinetic parameters, it is widely posited that constraints embedded in the catalytic mechanism enforce trade-offs between CO2 specificity, SC/O, and maximum carboxylation rate, kcat,C. However, the reasoning that established this view was based on data from ≈20 organisms. Here, we re-examine models of trade-offs in Rubisco catalysis using a data set from ≈300 organisms. Correlations between kinetic parameters are substantially attenuated in this larger data set, with the inverse relationship between kcat,C and SC/O being a key example. Nonetheless, measured kinetic parameters display extremely limited variation, consistent with a view of Rubisco as a highly constrained enzyme. More than 95% of kcat,C values are between 1 and 10 s-1, and no measured kcat,C exceeds 15 s-1. Similarly, SC/O varies by only 30% among Form I Rubiscos and

Many bacterial autotrophs rely on CO2concentrating mechanisms (CCMs) to assimilate carbon. Although many CCM proteins have been identified, including a 200+ MDa protein organelle called the carboxysome, a systematic screen of CCM components has not been carried out. Here, we performed a genome-wide barcoded transposon screen to identify essential and CCM-related genes in the ɣ-proteobacteriumH. neapolitanus. Our screen revealed an operon critical for CCM function which encodes a domain of unknown function (PFAM:PF10070) and putative cation transporter subunit (PFAM:PF00361). These two proteins, which we name DabA and DabB for “DABs accumulate bicarbonate,” function as a heterodimeric, energy-coupled inorganic carbon pump inE. coli. Furthermore, DabA binds zinc and has a an active site homologous to a β-carbonic anhydrase. Based on these results, we propose that DABs function as vectorial CAs coupled to cation gradients and serve as inorganic carbon pumps throughout prokaryotic phyla.

Rubisco is the most abundant enzyme in the biosphere and one of the best-characterized enzymes. Based on correlations between Rubisco kinetic parameters, it is widely posited that tradeoffs embedded in the catalytic mechanism constrain its specificity and maximum catalytic rate. However, the reasoning that established this view was based on data from ≈20 organisms. We re-examine these tradeoff models using a dataset from ≈300 organisms. Most correlations are substantially attenuated, with the inverse relationship between carboxylation kcat and specificity SC/O being a key example. Only one tradeoff model survives in our dataset. In this model, increasing catalytic efficiency (kcat/KM) for carboxylation requires increased catalytic efficiency for the competing oxygenation reaction, evidenced by strong power-law correlation between catalytic efficiencies. Our results imply that Rubisco evolution is constrained primarily by the physicochemical limits of O2/CO2 discrimination, which should reframe efforts to engineer this very central enzyme.

Thermodynamics dictates the structure and function of metabolism. Redox reactions drive cellular energy and material flow. Hence, accurately quantifying the thermodynamics of redox reactions should reveal key design principles that shape cellular metabolism. However, only a limited number of redox potentials have been measured experimentally, and mostly with inconsistent, poorly-reported experimental setups. Here, we develop a quantum chemistry approach for the calculation of redox potentials of biochemical reactions. We demonstrate that our method predicts experimentally measured potentials with unparalleled accuracy. We calculate the reduction potentials of all redox pairs that can be generated from biochemically relevant compounds and highlight fundamental thermodynamic trends that define cellular redox biochemistry. We further use the calculated potentials to address the question of why NAD/NADP are used as the primary cellular electron carriers, demonstrating how their physiological redox range specifically fits the reactions of central metabolism and minimizes the concentration of reactive carbonyls. The use of quantum chemistry tools, as demonstrated in this study, can revolutionize our understanding of key biochemical phenomena by enabling fast and accurate calculation of large datasets of thermodynamic values.

The microscopic world of a cell can be as alien to our human-centered intuition as the confinement of quarks within protons or the event horizon of a black hole. We are prone to thinking by analogy—Golgi cisternae stack like pancakes, red blood cells look like donuts—but very little in our human experience is truly comparable to the immensely crowded, membrane-subdivided interior of a eukaryotic cell or the intricately layered structures of a mammalian tissue. So in our daily efforts to understand how cells work, we are faced with a challenge: how do we develop intuition that works at the microscopic scale?

Electrosynthesis is a promising approach that enables the biological production of commodities, like fuels and fine chemicals, using renewably produced electricity. Several techniques have been proposed to mediate the transfer of electrons from the cathode to living cells. Of these, the electroproduction of formate as a mediator seems especially promising: formate is readily soluble, of low toxicity and can be produced at relatively high efficiency and at reasonable current density. While organisms that are capable of formatotrophic growth, i.e. growth on formate, exist naturally, they are generally less suitable for bulk cultivation and industrial needs. Hence, it may be helpful to engineer a model organism of industrial relevance, such as E. coli, for growth on formate. There are numerous metabolic pathways that can potentially support formatotrophic growth. Here we analyze these diverse pathways according to various criteria including biomass yield, thermodynamic favorability, chemical motive force, kinetics and the practical challenges posed by their expression. We find that the reductive glycine pathway, composed of the tetrahydrofolate system, the glycine cleavage system, serine hydroxymethyltransferase and serine deaminase, is a promising candidate to support electrosynthesis in E. coli. The approach presented here exemplifies how combining different computational approaches into a systematic analysis methodology provides assistance in redesigning metabolism. This article is part of a Special Issue entitled: Metals in Bioenergetics and Biomimetics Systems.

Many carbon-fixing bacteria rely on a CO2 concentrating mechanism (CCM) to elevate the CO2 concentration around the carboxylating enzyme ribulose bisphosphate carboxylase/oxygenase (RuBisCO). The CCM is postulated to simultaneously enhance the rate of carboxylation and minimize oxygenation, a competitive reaction with O2 also catalyzed by RuBisCO. To achieve this effect, the CCM combines two features: active transport of inorganic carbon into the cell and colocalization of carbonic anhydrase and RuBisCO inside proteinaceous microcompartments called carboxysomes. Understanding the significance of the various CCM components requires reconciling biochemical intuition with a quantitative description of the system. To this end, we have developed a mathematical model of the CCM to analyze its energetic costs and the inherent intertwining of physiology and pH. We find that intracellular pH greatly affects the cost of inorganic carbon accumulation. At low pH the inorganic carbon pool contains more of the highly cell-permeable H2CO3, necessitating a substantial expenditure of energy on transport to maintain internal inorganic carbon levels. An intracellular pH ≈8 reduces leakage, making the CCM significantly more energetically efficient. This pH prediction coincides well with our measurement of intracellular pH in a model cyanobacterium. We also demonstrate that CO2 retention in the carboxysome is necessary, whereas selective uptake of HCO3 (-) into the carboxysome would not appreciably enhance energetic efficiency. Altogether, integration of pH produces a model that is quantitatively consistent with cyanobacterial physiology, emphasizing that pH cannot be neglected when describing biological systems interacting with inorganic carbon pools.

The cyanobacterium Synechococcus elongatus relies upon photosynthesis to drive metabolism and growth. During darkness, Synechococcus stops growing, derives energy from its glycogen stores, and greatly decreases rates of macromolecular synthesis via unknown mechanisms. Here, we show that the stringent response, a stress response pathway whose genes are conserved across bacteria and plant plastids, contributes to this dark adaptation. Levels of the stringent response alarmone guanosine 3'-diphosphate 5'-diphosphate (ppGpp) rise after a shift from light to dark, indicating that darkness triggers the same response in cyanobacteria as starvation in heterotrophic bacteria. High levels of ppGpp are sufficient to stop growth and dramatically alter many aspects of cellular physiology, including levels of photosynthetic pigments and polyphosphate, DNA content, and the rate of translation. Cells unable to synthesize ppGpp display pronounced growth defects after exposure to darkness. The stringent response regulates expression of a number of genes in Synechococcus, including ribosomal hibernation promoting factor (hpf), which causes ribosomes to dimerize in the dark and may contribute to decreased translation. Although the metabolism of Synechococcus differentiates it from other model bacterial systems, the logic of the stringent response remains remarkably conserved, while at the same time having adapted to the unique stresses of the photosynthetic lifestyle.

The clustered, regularly interspaced, short palindromic repeats (CRISPR)-associated protein Cas9 from Streptococcus pyogenes is an RNA-guided DNA endonuclease with widespread utility for genome modification. However, the structural constraints limiting the engineering of Cas9 have not been determined. Here we experimentally profile Cas9 using randomized insertional mutagenesis and delineate hotspots in the structure capable of tolerating insertions of a PDZ domain without disruption of the enzyme's binding and cleavage functions. Orthogonal domains or combinations of domains can be inserted into the identified sites with minimal functional consequence. To illustrate the utility of the identified sites, we construct an allosterically regulated Cas9 by insertion of the estrogen receptor-α ligand-binding domain. This protein showed robust, ligand-dependent activation in prokaryotic and eukaryotic cells, establishing a versatile one-component system for inducible and reversible Cas9 activation. Thus, domain insertion profiling facilitates the rapid generation of new Cas9 functionalities and provides useful data for future engineering of Cas9.

Thermodynamics impose a major constraint on the structure of metabolic pathways. Here, we use carbon fixation pathways to demonstrate how thermodynamics shape the structure of pathways and determine the cellular resources they consume. We analyze the energetic profile of prototypical reactions and show that each reaction type displays a characteristic change in Gibbs energy. Specifically, although carbon fixation pathways display a considerable structural variability, they are all energetically constrained by two types of reactions: carboxylation and carboxyl reduction. In fact, all adenosine triphosphate (ATP) molecules consumed by carbon fixation pathways – with a single exception – are used, directly or indirectly, to power one of these unfavorable reactions. When an indirect coupling is employed, the energy released by ATP hydrolysis is used to establish another chemical bond with high energy of hydrolysis, e.g. a thioester. This bond is cleaved by a downstream enzyme to energize an unfavorable reaction. Notably, many pathways exhibit reduced ATP requirement as they couple unfavorable carboxylation or carboxyl reduction reactions to exergonic reactions other than ATP hydrolysis. In the most extreme example, the reductive acetyl coenzyme A (acetyl-CoA) pathway bypasses almost all ATP-consuming reactions. On the other hand, the reductive pentose phosphate pathway appears to be the least ATP-efficient because it is the only carbon fixation pathway that invests ATP in metabolic aims other than carboxylation and carboxyl reduction. Altogether, our analysis indicates that basic thermodynamic considerations accurately predict the resource investment required to support a metabolic pathway and further identifies biochemical mechanisms that can decrease this requirement.

Motivation: The laws of thermodynamics describe a direct, quantitative relationship between metabolite concentrations and reaction directionality. Despite great efforts, thermodynamic data suffer from limited coverage, scattered accessibility and non-standard annotations. We present a framework for unifying thermodynamic data from multiple sources and demonstrate two new techniques for extrapolating the Gibbs energies of unmeasured reactions and conditions. Results: Both methods account for changes in cellular conditions (pH, ionic strength, etc.) by using linear regression over the ΔG○ of pseudoisomers and reactions. The Pseudoisomeric Reactant Contribution method systematically infers compound formation energies using measured K′ and pKa data. The Pseudoisomeric Group Contribution method extends the group contribution method and achieves a high coverage of unmeasured reactions. We define a continuous index that predicts the reversibility of a reaction under a given physiological concentration range. In the characteristic physiological range 3μM–3mM, we find that roughly half of the reactions in Escherichia coli's metabolism are reversible. These new tools can increase the accuracy of thermodynamic-based models, especially in non-standard pH and ionic strengths. The reversibility index can help modelers decide which reactions are reversible in physiological conditions. Availability: Freely available on the web at: http://equilibrator.weizmann.ac.il. Website implemented in Python, MySQL, Apache and Django, with all major browsers supported. The framework is open-source (code.google.com/p/milo-lab), implemented in pure Python and tested mainly on Linux. Contact: ron.milo@weizmann.ac.il Supplementary Information: Supplementary data are available at Bioinformatics online.

The laws of thermodynamics constrain the action of biochemical systems. However, thermodynamic data on biochemical compounds can be difficult to find and is cumbersome to perform calculations with manually. Even simple thermodynamic questions like 'how much Gibbs energy is released by ATP hydrolysis at pH 5?' are complicated excessively by the search for accurate data. To address this problem, eQuilibrator couples a comprehensive and accurate database of thermodynamic properties of biochemical compounds and reactions with a simple and powerful online search and calculation interface. The web interface to eQuilibrator (http://equilibrator.weizmann.ac.il) enables easy calculation of Gibbs energies of compounds and reactions given arbitrary pH, ionic strength and metabolite concentrations. The eQuilibrator code is open-source and all thermodynamic source data are freely downloadable in standard formats. Here we describe the database characteristics and implementation and demonstrate its use.

Thermodynamics dictates the structure and function of metabolism. Redox reactions drive cellular energy and material flow. Hence, accurately quantifying the thermodynamics of redox reactions should reveal design principles that shape cellular metabolism. However, only few redox potentials have been measured, and mostly with inconsistent experimental setups. Here, we develop a quantum chemistry approach to calculate redox potentials of biochemical reactions and demonstrate our method predicts experimentally measured potentials with unparalleled accuracy. We then calculate the potentials of all redox pairs that can be generated from biochemically relevant compounds and highlight fundamental trends in redox biochemistry. We further address the question of why NAD/NADP are used as primary electron carriers, demonstrating how their physiological potential range fits the reactions of central metabolism and minimizes the concentration of reactive carbonyls. The use of quantum chemistry can revolutionize our understanding of biochemical phenomena by enabling fast and accurate calculation of thermodynamic values., PLoS Computational Biology, 14 (10), ISSN:1553-734X, ISSN:1553-7358

Single-fluorescent protein biosensors (SFPBs) are an important class of probes that enable the single-cell quantification of analytes in vivo. Despite advantages over other detection technologies, their use has been limited by the inherent challenges of their construction. Specifically, the rational design of green fluorescent protein (GFP) insertion into a ligand-binding domain, generating the requisite allosteric coupling, remains a rate-limiting step. Here, we describe an unbiased approach, termed domain-insertion profiling with DNA sequencing (DIP-seq), that combines the rapid creation of diverse libraries of potential SFPBs and high-throughput activity assays to identify functional biosensors. As a proof of concept, we construct an SFPB for the important regulatory sugar trehalose. DIP-seq analysis of a trehalose-binding-protein reveals allosteric hotspots for GFP insertion and results in high-dynamic range biosensors that function robustly in vivo. Taken together, DIP-seq simultaneously accelerates metabolite biosensor construction and provides a novel tool for interrogating protein allostery., In the construction of single fluorescent protein biosensors, selection of the insertion point of a fluorescent protein into a ligand-binding domain is a rate-limiting step. Here, the authors develop an unbiased, high-throughput approach, called domain insertion profiling with DNA sequencing (DIP-seq), to generate a novel trehalose biosensor.

Proteomics techniques generate an avalanche of data and promise to satisfy biologists' long-held desire to measure absolute protein abundances on a genome-wide scale. However, can this knowledge be translated into a clearer picture of how cells invest their protein resources? This article aims to give a broad perspective on the composition of proteomes as gleaned from recent quantitative proteomics studies. We describe proteomaps, an approach for visualizing the composition of proteomes with a focus on protein abundances and functions. In proteomaps, each protein is shown as a polygon-shaped tile, with an area representing protein abundance. Functionally related proteins appear in adjacent regions. General trends in proteomes, such as the dominance of metabolism and protein production, become easily visible. We make interactive visualizations of published proteome datasets accessible at www.proteomaps.net. We suggest that evaluating the way protein resources are allocated by various organisms and cell types in different conditions will sharpen our understanding of how and why cells regulate the composition of their proteomes.

Michaelis and Menten’s mechanism for enzymatic catalysis is remarkable both in its simplicity and its wide applicability. The extension for reversible processes, as done by Haldane, makes it even more relevant as most enzymes catalyze reactions that are reversible in nature and carry in vivo flux in both directions. Here, we decompose the reversible Michaelis–Menten equation into three terms, each with a clear physical meaning: catalytic capacity, substrate saturation and thermodynamic driving force. This decomposition facilitates a better understanding of enzyme kinetics and highlights the relationship between thermodynamics and kinetics, a relationship which is often neglected. We further demonstrate how our separable rate law can be understood from different points of view, shedding light on factors shaping enzyme catalysis.

Contrary to the textbook portrayal of glycolysis as a single pathway conserved across all domains of life, not all sugar-consuming organisms use the canonical Embden–Meyerhoff–Parnass (EMP) glycolytic pathway. Prokaryotic glucose metabolism is particularly diverse, including several alternative glycolytic pathways, the most common of which is the Entner–Doudoroff (ED) pathway. The prevalence of the ED pathway is puzzling as it produces only one ATP per glucose—half as much as the EMP pathway. We argue that the diversity of prokaryotic glucose metabolism may reflect a tradeoff between a pathway’s energy (ATP) yield and the amount of enzymatic protein required to catalyze pathway flux. We introduce methods for analyzing pathways in terms of thermodynamics and kinetics and show that the ED pathway is expected to require several-fold less enzymatic protein to achieve the same glucose conversion rate as the EMP pathway. Through genomic analysis, we further show that prokaryotes use different glycolytic pathways depending on their energy supply. Specifically, energy-deprived anaerobes overwhelmingly rely upon the higher ATP yield of the EMP pathway, whereas the ED pathway is common among facultative anaerobes and even more common among aerobes. In addition to demonstrating how protein costs can explain the use of alternative metabolic strategies, this study illustrates a direct connection between an organism’s environment and the thermodynamic and biochemical properties of the metabolic pathways it employs.

What governs the concentrations of metabolites within living cells? Beyond specific metabolic and enzymatic considerations, are there global trends that affect their values? We hypothesize that the physico-chemical properties of metabolites considerably affect their in-vivo concentrations. The recently achieved experimental capability to measure the concentrations of many metabolites simultaneously has made the testing of this hypothesis possible. Here, we analyze such recently available data sets of metabolite concentrations within E. coli, S. cerevisiae, B. subtilis and human. Overall, these data sets encompass more than twenty conditions, each containing dozens (28-108) of simultaneously measured metabolites. We test for correlations with various physico-chemical properties and find that the number of charged atoms, non-polar surface area, lipophilicity and solubility consistently correlate with concentration. In most data sets, a change in one of these properties elicits a ∼100 fold increase in metabolite concentrations. We find that the non-polar surface area and number of charged atoms account for almost half of the variation in concentrations in the most reliable and comprehensive data set. Analyzing specific groups of metabolites, such as amino-acids or phosphorylated nucleotides, reveals even a higher dependence of concentration on hydrophobicity. We suggest that these findings can be explained by evolutionary constraints imposed on metabolite concentrations and discuss possible selective pressures that can account for them. These include the reduction of solute leakage through the lipid membrane, avoidance of deleterious aggregates and reduction of non-specific hydrophobic binding. By highlighting the global constraints imposed on metabolic pathways, future research could shed light onto aspects of biochemical evolution and the chemical constraints that bound metabolic engineering efforts., Author Summary What governs the identity and concentrations of metabolites within living cells? The first part of this question has received much attention. Organisms were found to qualitatively prefer hydrophilic and charged metabolites, a phenomenon that was explained to be a result of constraints imposed by contemporary as well as archaic metabolism. However, among the metabolites that are used, a quantitative preference has never been analyzed systematically. Here we use the most comprehensive data sets of metabolite concentrations available to explore such trends. We find that in various organisms and growth conditions, living cells minimize the concentrations of non-polar, un-charged metabolites. More specifically, metabolites' hydrophobicity alters concentrations by two orders of magnitudes on average and explains up to half of the variation of metabolite concentrations within cells. We suggest that this can be attributed to an evolutionary pressure to avoid an unspecific hydrophobic effect: the preference of hydrophobic surfaces in an aqueous environment to adhere to other hydrophobic surfaces. Our findings shed light on the evolution of the internal makeup of living cells and can assist in establishing metabolic models that support synthetic biology and metabolic engineering efforts.

Background The availability of microarrays measuring thousands of genes simultaneously across hundreds of biological conditions represents an opportunity to understand both individual biological pathways and the integrated workings of the cell. However, translating this amount of data into biological insight remains a daunting task. An important initial step in the analysis of microarray data is clustering of genes with similar behavior. A number of classical techniques are commonly used to perform this task, particularly hierarchical and K-means clustering, and many novel approaches have been suggested recently. While these approaches are useful, they are not without drawbacks; these methods can find clusters in purely random data, and even clusters enriched for biological functions can be skewed towards a small number of processes (e.g. ribosomes). Results We developed Nearest Neighbor Networks (NNN), a graph-based algorithm to generate clusters of genes with similar expression profiles. This method produces clusters based on overlapping cliques within an interaction network generated from mutual nearest neighborhoods. This focus on nearest neighbors rather than on absolute distance measures allows us to capture clusters with high connectivity even when they are spatially separated, and requiring mutual nearest neighbors allows genes with no sufficiently similar partners to remain unclustered. We compared the clusters generated by NNN with those generated by eight other clustering methods. NNN was particularly successful at generating functionally coherent clusters with high precision, and these clusters generally represented a much broader selection of biological processes than those recovered by other methods. Conclusion The Nearest Neighbor Networks algorithm is a valuable clustering method that effectively groups genes that are likely to be functionally related. It is particularly attractive due to its simplicity, its success in the analysis of large datasets, and its ability to span a wide range of biological functions with high precision.

In metabolism research, thermodynamics is usually used to determine the directionality of a reaction or the feasibility of a pathway. However, the relationship between thermodynamic potentials and fluxes is not limited to questions of directionality: thermodynamics also affects the kinetics of reactions through the flux-force relationship, which states that the logarithm of the ratio between the forward and reverse fluxes is directly proportional to the change in Gibbs energy due to a reaction (ΔrG′). Accordingly, if an enzyme catalyzes a reaction with a ΔrG′ of -5.7 kJ/mol then the forward flux will be roughly ten times the reverse flux. As ΔrG′ approaches equilibrium (ΔrG′ = 0 kJ/mol), exponentially more enzyme counterproductively catalyzes the reverse reaction, reducing the net rate at which the reaction proceeds. Thus, the enzyme level required to achieve a given flux increases dramatically near equilibrium. Here, we develop a framework for quantifying the degree to which pathways suffer these thermodynamic limitations on flux. For each pathway, we calculate a single thermodynamically-derived metric (the Max-min Driving Force, MDF), which enables objective ranking of pathways by the degree to which their flux is constrained by low thermodynamic driving force. Our framework accounts for the effect of pH, ionic strength and metabolite concentration ranges and allows us to quantify how alterations to the pathway structure affect the pathway's thermodynamics. Applying this methodology to pathways of central metabolism sheds light on some of their features, including metabolic bypasses (e.g., fermentation pathways bypassing substrate-level phosphorylation), substrate channeling (e.g., of oxaloacetate from malate dehydrogenase to citrate synthase), and use of alternative cofactors (e.g., quinone as an electron acceptor instead of NAD). The methods presented here place another arrow in metabolic engineers' quiver, providing a simple means of evaluating the thermodynamic and kinetic quality of different pathway chemistries that produce the same molecules., Author Summary Given data about enzyme kinetics and reaction thermodynamics, traditional metabolic control analysis (MCA) can pinpoint the enzymes whose expression will have the largest effect on steady-state flux through the pathway. These analyses can aid experimentalists in tuning enzyme expression levels along a metabolic pathway. In this work, we offer a framework that is complementary to MCA. Rather than focusing on the relationship between enzyme levels and pathway flux, we examine a pathway's stoichiometry and thermodynamics and ask whether it is likely to support high flux in cellular conditions. Our framework calculates a single thermodynamically-derived metric (the MDF) for each pathway, which is convenient for selecting the promising pathways from a large collection. This approach has several advantages. First, enzyme kinetic properties are laborious to measure and differ between organisms and isozymes, but no kinetic data is required to calculate the MDF. Second, as our framework accounts for pH, ionic strength and allowed concentration ranges, it is simple to model the effect of these parameters on the MDF. Finally, as it can be difficult to control the exact expression level of enzymes within cells, the MDF helps identify alternative pathways that are less sensitive to the levels of their constituent enzymes.

Protein levels are a dominant factor shaping natural and synthetic biological systems. Although proper functioning of metabolic pathways relies on precise control of enzyme levels, the experimental ability to balance the levels of many genes in parallel is a major outstanding challenge. Here, we introduce a rapid and modular method to span the expression space of several proteins in parallel. By combinatorially pairing genes with a compact set of ribosome-binding sites, we modulate protein abundance by several orders of magnitude. We demonstrate our strategy by using a synthetic operon containing fluorescent proteins to span a 3D color space. Using the same approach, we modulate a recombinant carotenoid biosynthesis pathway in Escherichia coli to reveal a diversity of phenotypes, each characterized by a distinct carotenoid accumulation profile. In a single combinatorial assembly, we achieve a yield of the industrially valuable compound astaxanthin 4-fold higher than previously reported. The methodology presented here provides an efficient tool for exploring a high-dimensional expression space to locate desirable phenotypes.

CO 2 is converted into biomass almost solely by the enzyme rubisco. The poor carboxylation properties of plant rubiscos have led to efforts that made it the most kinetically characterized enzyme, yet these studies focused on 100) were active in vitro, with the fastest having a turnover number of 22 ± 1 s−1—sixfold faster than the median plant rubisco and nearly twofold faster than the fastest measured rubisco to date. Unlike rubiscos from plants and cyanobacteria, the fastest variants discovered here are homodimers and exhibit a much simpler folding and activation kinetics. Our pipeline can be utilized to explore the kinetic space of other enzymes of interest, allowing us to get a better view of the biosynthetic potential of the biosphere., Metagenomic and biochemical analysis of previously uncharacterized naturally‐occurring form‐II and II/III rubiscos leads to identification of an enzyme with the fastest CO2 fixation rate described to date.