Publications

Autotrophy is the basis for complex life on Earth. Central to this process is rubiscothe enzyme that catalyzes almost all carbon fixation on the planet. Yet, with only a small fraction of rubisco diversity kinetically characterized so far, the underlying biological factors driving the evolution of fast rubiscos in nature remain unclear. We conducted a high-throughput kinetic characterization of over 100 bacterial form I rubiscos, the most ubiquitous group of rubisco sequences in nature, to uncover the determinants of rubiscos carboxylation velocity. We show that the presence of a carboxysome CO₂ concentrating mechanism correlates with faster rubiscos with a median fivefold higher rate. In contrast to prior studies, we find that rubiscos originating from α-cyanobacteria exhibit the highest carboxylation rates among form I enzymes (≈10 s⁻¹ median versus

Synthetic autotrophy is a promising avenue to sustainable bioproduction from CO₂. Here, we use iterative laboratory evolution to generate several distinct autotrophic strains. Util-ising this genetic diversity, we identify that just three mutations are sufficient for Escherichia coli to grow autotrophically, when introduced alongside non-native energy (formate dehydrogenase) and carbon-fixing (RuBisCO, phosphoribulokinase, carbonic anhydrase) modules. The mutated genes are involved in glycolysis (pgi), central-carbon regulation (crp), and RNA transcription (rpoB). The pgi mutation reduces the enzymes activity, thereby stabilising the carbon-fixing cycle by capping a major branching flux. For the other two mutations, we observe down-regulation of several metabolic pathways and increased expression of native genes associated with the carbon-fixing module (rpiB) and the energy module (fdoGH), as well as an increased ratio of NADH/NAD⁺-the cycles electron-donor. This study demonstrates the malleability of metabolism and its capacity to switch trophic modes using only a small number of genetic changes and could facilitate transforming other hetero-trophic organisms into autotrophs.

The immune system is a complex network of cells with critical functions in health and disease. However, a comprehensive census of the cells comprising the immune system is lacking. Here, we estimated the abundance of the primary immune cell types throughout all tissues in the human body. We conducted a literature survey and integrated data from multiplexed imaging and methylome-based deconvolution. We also considered cellular mass to determine the distribution of immune cells in terms of both number and total mass. Our results indicate that the immune system of a reference 73 kg man consists of 1.8 x 10¹² cells (95% CI 1.5-2.3 x 10¹²), weighing 1.2 kg (95% CI 0.8-1.9). Lymphocytes constitute 40% of the total number of immune cells and 15% of the mass and are mainly located in the lymph nodes and spleen. Neutrophils account for similar proportions of both the number and total mass of immune cells, with most neutrophils residing in the bone marrow. Macrophages, present in most tissues, account for 10% of immune cells but contribute nearly 50% of the total cellular mass due to their large size. The quantification of immune cells within the human body presented here can serve to understand the immune function better and facilitate quantitative modeling of this vital system.

Following evidence of waning immunity against both infection and severe disease after 2 doses of the BNT162b2 vaccine, Israel began administering a 3rd BNT162b2 dose (booster) in July 2021. Recent studies showed that the 3rd dose provides a much lower protection against infection with the Omicron variant compared to the Delta variant and that this protection wanes quickly. However, there is little evidence regarding the protection of the 3rd dose against Omicron (BA.1/BA.2) severe disease. In this study, we estimate the preservation of immunity from severe disease up to 7 months after receiving the booster dose. We calculate rates of severe SARS-CoV-2 disease between groups of individuals aged 60 and above, comparing those who received two doses at least 4 months previously to those who received the 3rd dose (stratified by the time from vaccination), and to those who received a 4th dose. The analysis shows that protection conferred by the 3rd dose against Omicron severe disease did not wane over a 7-month period. Moreover, a 4th dose further improved protection, with a severe disease rate approximately 3-fold lower than in the 3-dose cohorts.

Insects and other arthropods are central to terrestrial ecosystems. However, data are lacking regarding their global population abundance. We synthesized thousands of evaluations from around 500 sites worldwide, estimating the absolute biomass and abundance of terrestrial arthropods across different taxa and habitats. We found that there are approximate to 1 x 10(19) (twofold uncertainty range) soil arthropods on Earth, approximate to 95% of which are soil mites and springtails. The soil contains approximate to 200 (twofold uncertainty range) million metric tons (Mt) of dry biomass. Termites contribute approximate to 40% of the soil biomass, much more than ants at approximate to 10%. Our estimate for the global biomass of above-ground arthropods is more uncertain, highlighting a knowledge gap that future research should aim to close. We estimate the combined dry biomass of all terrestrial arthropods at approximate to 300 Mt (uncertainty range, 100 to 500), similar to the mass of humanity and its livestock. These estimates enhance the quantitative understanding of arthropods in terrestrial ecosystems and provide an initial holistic benchmark on their decline.

The 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 reports 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, and links to the underlying data, as well as a permanent identifier called a 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 facilitates access to carefully curated data, acting as a quantitative resource pertaining to the myriad ways in which humans have an impact on the Earth, for practicing scientists, the general public, and those involved in education for sustainable development alike. We outline the structure of the database, describe our curation procedures, and 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 intimate interconnections.

Quantifying the temporal dynamics of infectiousness of individuals infected with SARS-CoV-2 is crucial for understanding the spread of COVID-19 and for evaluating the effectiveness of mitigation strategies. Many studies have estimated the infectiousness profile using observed serial intervals. However, statistical and epidemiological biases could lead to underestimation of the duration of infectiousness. We correct for these biases by curating data from the initial outbreak of the pandemic in China (when mitigation was minimal), and find that the infectiousness profile of the original strain is longer than previously thought. Sensitivity analysis shows our results are robust to model structure, assumed growth rate and potential observational biases. Although unmitigated transmission data is lacking for variants of concern (VOCs), previous analyses suggest that the alpha and delta variants have faster within-host kinetics, which we extrapolate to crude estimates of variant-specific unmitigated generation intervals. Knowing the unmitigated infectiousness profile of infected individuals can inform estimates of the effectiveness of isolation and quarantine measures. The framework presented here can help design better quarantine policies in early stages of future epidemics.

Israel began administering a BNT162b2 booster dose to restore protection following the waning of the 2-dose vaccine. Biological studies have shown that a "fresh" booster dose leads to increased antibody levels compared to a fresh 2-dose vaccine, which may suggest increased effectiveness. To compare the real-world effectiveness of a fresh (up to 60 days) booster dose with that of a fresh 2-dose vaccine, we took advantage of a quasi-experimental study that compares populations that were eligible to receive the vaccine at different times due to age-dependent policies. Specifically, we compared the confirmed infection rates in adolescents aged 12-14 (215,653 individuals) who received the 2-dose vaccine and in adolescents aged 16-18 (103,454 individuals) who received the booster dose. Our analysis shows that the confirmed infection rate was lower by a factor of 3.7 (95% CI: 2.7 to 5.2) in the booster group.

BACKGROUND On July 30, 2021, the administration of a third (booster) dose of the BNT162b2 messenger RNA vaccine (Pfizer-BioNTech) was approved in Israel for persons who were 60 years of age or older and who had received a second dose of vaccine at least 5 months earlier. Data are needed regarding the effect of the booster dose on the rate of confirmed coronavirus 2019 disease (Covid-19) and the rate of severe illness. METHODS We extracted data for the period from July 30 through August 31, 2021, from the Israeli Ministry of Health database regarding 1,137,804 persons who were 60 years of age or older and had been fully vaccinated (i.e., had received two doses of BNT162b2) at least 5 months earlier. In the primary analysis, we compared the rate of confirmed Covid-19 and the rate of severe illness between those who had received a booster injection at least 12 days earlier (booster group) and those who had not received a booster injection (nonbooster group). In a secondary analysis, we evaluated the rate of infection 4 to 6 days after the booster dose as compared with the rate at least 12 days after the booster. In all the analyses, we used Poisson regression after adjusting for possible confounding factors. RESULTS At least 12 days after the booster dose, the rate of confirmed infection was lower in the booster group than in the nonbooster group by a factor of 11.3 (95% confidence interval [CI], 10.4 to 12.3); the rate of severe illness was lower by a factor of 19.5 (95% CI, 12.9 to 29.5). In a secondary analysis, the rate of confirmed infection at least 12 days after vaccination was lower than the rate after 4 to 6 days by a factor of 5.4 (95% CI, 4.8 to 6.1). CONCLUSIONS In this study involving participants who were 60 years of age or older and had received two doses of the BNT162b2 vaccine at least 5 months earlier, we found that the rates of confirmed Covid-19 and severe illness were substantially lower among those who received a booster (third) dose of the BNT162b2 vaccine.

Population growth and changes in dietary patterns place an ever-growing pressure on the environment. Feeding the world within sustainable boundaries therefore requires revolutionizing the way we harness natural resources. Microbial biomass can be cultivated to yield protein-rich feed and food supplements, collectively termed single-cell protein (SCP). Yet, we still lack a quantitative comparison between traditional agriculture and photovoltaic-driven SCP systems in terms of land use and energetic efficiency. Here, we analyze the energetic efficiency of harnessing solar energy to produce SCP from air and water. Our model includes photovoltaic electricity generation, direct air capture of carbon dioxide, electrosynthesis of an electron donor and/or carbon source for microbial growth (hydrogen, formate, or methanol), microbial cultivation, and the processing of biomass and proteins. We show that, per unit of land, SCP production can reach an over 10-fold higher protein yield and at least twice the caloric yield compared with any staple crop. Altogether, this quantitative analysis offers an assessment of the future potential of photovoltaic-driven microbial foods to supplement conventional agricultural production and support resource-efficient protein supply on a global scale.

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 10⁹ to 10¹¹ 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.

Humanity has become a dominant force in shaping the face of Earth(1-9). An emerging question is how the overall material output of human activities compares to the overall natural biomass. Here we quantify the human-made mass, referred to as 'anthropogenic mass', and compare it to the overall living biomass on Earth, which currently equals approximately 1.1 teratonnes(10,11). We find that Earth is exactly at the crossover point; in the year 2020 (+/- 6), the anthropogenic mass, which has recently doubled roughly every 20 years, will surpass all global living biomass. On average, for each person on the globe, anthropogenic mass equal to more than his or her bodyweight is produced every week. This quantification of the human enterprise gives a mass-based quantitative and symbolic characterization of the human-induced epoch of the Anthropocene.

Many photosynthetic organisms employ a CO2 concentrating mechanism (CCM) to increase the rate of CO2 fixation via the Calvin cycle. CCMs catalyze approximate to 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 Escherichia coli strain engineered to depend on rubisco carboxylation for growth. Expression of 20 CCM genes enabled E. coli to grow by fixing CO2 from 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 CO2 assimilation in diverse organisms.

The current SARS-CoV-2 pandemic is a harsh reminder of the fact that, whether in a single human host or a wave of infection across continents, viral dynamics is often a story about the numbers. In this snapshot, our aim is to provide a one-stop, curated graphical source for the key numbers that help us understand the virus driving our current global crisis. The discussion is framed around two broad themes: 1) the biology of the virus itself and 2) the characteristics of the infection of a single human host. Our one-page summary provides the key numbers pertaining to SARS-CoV-2, based mostly on peer-reviewed literature. The numbers reported in summary format are substantiated by the annotated references below. Readers are urged to remember that much uncertainty remains and knowledge of this pandemic and the virus driving it is rapidly evolving. In the paragraphs below we provide \u201cback of the envelope\u201d calculations that exemplify the insights that can be gained from knowing some key numbers and using quantitative logic. These calculations serve to improve our intuition through sanity checks, but do not replace detailed epidemiological analysis.

The living world is largely divided into autotrophs that convert CO2 into biomass and heterotrophs that consume organic compounds. In spite of widespread interest in renewable energy storage and more sustainable food production, the engineering of industrially relevant heterotrophic model organisms to use CO2 as their sole carbon source has so far remained an outstanding challenge. Here, we report the achievement of this transformation on laboratory timescales. We constructed and evolved Escherichia coli to produce all its biomass carbon from CO2. Reducing power and energy, but not carbon, are supplied via the one-carbon molecule formate, which can be produced electrochemically. Rubisco and phosphoribulokinase were co-expressed with formate dehydrogenase to enable CO2 fixation and reduction via the Calvin-Benson-Bassham cycle. Autotrophic growth was achieved following several months of continuous laboratory evolution in a chemostat under intensifying organic carbon limitation and confirmed via isotopic labeling.

Biofilms are a fundamental form of microbial life and occur in diverse environments, ranging from the mammalian gut to deep subsurface rocks. It is often claimed that most bacteria and archaea live in biofilms, but this claim awaits quantification. Recent updates on global microbial cell numbers prompt a revisiting of this question.

A census of the biomass on Earth is key for understanding the structure and dynamics of the biosphere. However, a global, quantitative view of how the biomass of different taxa compare with one another is still lacking. Here, we assemble the overall biomass composition of the biosphere, establishing a census of the ≈550 gigatons of carbon (Gt C) of biomass distributed among all of the kingdoms of life. We find that the kingdoms of life concentrate at different locations on the planet; plants (≈450 Gt C, the dominant kingdom) are primarily terrestrial, whereas animals (≈2 Gt C) are mainly marine, and bacteria (≈70 Gt C) and archaea (≈7 Gt C) are predominantly located in deep subsurface environments. We show that terrestrial biomass is about two orders of magnitude higher than marine biomass and estimate a total of ≈6 Gt C of marine biota, doubling the previous estimated quantity. Our analysis reveals that the global marine biomass pyramid contains more consumers than producers, thus increasing the scope of previous observations on inverse food pyramids. Finally, we highlight that the mass of humans is an order of magnitude higher than that of all wild mammals combined and report the historical impact of humanity on the global biomass of prominent taxa, including mammals, fish, and plants.

Carbon fixation is the gateway of inorganic carbon into the biosphere. Our ability to engineer carbon fixation pathways in living organisms is expected to play a crucial role in the quest towards agricultural and energetic sustainability. Recent successes to introduce non-native carbon fixation pathways into heterotrophic hosts offer novel platforms for manipulating these pathways in genetically malleable organisms. Here, we focus on past efforts and future directions for engineering the dominant carbon fixation pathway in the biosphere, the Calvin-Benson cycle, into the well-known model organism Escherichia coli. We describe how central carbon metabolism of this heterotrophic bacterium can be manipulated to allow directed evolution of carbon fixing enzymes. Finally, we highlight future directions towards synthetic autotrophy.

Viruses are incapable of autonomous energy production. Although many experimental studies make it clear that viruses are parasitic entities that hijack the molecular resources of the host, a detailed estimate for the energetic cost of viral synthesis is largely lacking. To quantify the energetic cost of viruses to their hosts, we enumerated the costs associated with two very distinct but representative DNA and RNA viruses, namely, T4 and influenza. We found that, for these viruses, translation of viral proteins is the most energetically expensive process. Interestingly, the costs of building a T4 phage and a single influenza virus are nearly the same. Due to influenza's higher burst size, however, the overall cost of a T4 phage infection is only 2.3% of the cost of an influenza infection. The costs of these infections relative to their host's estimated energy budget during the infection reveal that a T4 infection consumes about a third of its host's energy budget, whereas an influenza infection consumes only ≈ 1%. Building on our estimates for T4, we show how the energetic costs of double-stranded DNA phages scale with the capsid size, revealing that the dominant cost of building a virus can switch from translation to genome replication above a critical size. Last, using our predictions for the energetic cost of viruses, we provide estimates for the strengths of selection and genetic drift acting on newly incorporated genetic elements in viral genomes, under conditions of energy limitation.

Feeding a growing population while minimizing environmental degradation is a global challenge requiring thoroughly rethinking food production and consumption. Dietary choices control food availability and natural resource demands. In particular, reducing or avoiding consumption of low production efficiency animal-based products can spare resources that can then yield more food. In quantifying the potential food gains of specific dietary shifts, most earlier research focused on calories, with less attention to other important nutrients, notably protein. Moreover, despite the well-known environmental burdens of livestock, only a handful of national level feed-to-food conversion efficiency estimates of dairy, beef, poultry, pork, and eggs exist. Yet such high level estimates are essential for reducing diet related environmental impacts and identifying optimal food gain paths. Here we quantify caloric and protein conversion efficiencies for US livestock categories. We then use these efficiencies to calculate the food availability gains expected from replacing beef in the US diet with poultry, a more efficient meat, and a plant-based alternative. Averaged over all categories, caloric and protein efficiencies are 7%-8%. At 3% in both metrics, beef is by far the least efficient. We find that reallocating the agricultural land used for beef feed to poultry feed production can meet the caloric and protein demands of ≈120 and ≈140 million additional people consuming the mean American diet, respectively, roughly 40% of current US population.

Livestock farming incurs large and varied environmental burdens, dominated by beef. Replacing beef with resource efficient alternatives is thus potentially beneficial, but may conflict with nutritional considerations. Here we show that protein-equivalent plant based alternatives to the beef portion of the mean American diet are readily devisible, and offer mostly improved nutritional profile considering the full lipid profile, key vitamins, minerals, and micronutrients. We then show that replacement diets require on average only 10% of land, 4% of greenhouse gas (GHG) emissions, and 6% of reactive nitrogen (Nr) compared to what the replaced beef diet requires. Applied to 320 million Americans, the beef-to-plant shift can save 91 million cropland acres (and 770 million rangeland acres), 278 million metric ton CO_2e, and 3.7 million metric ton Nr annually. These nationwide savings are 27%, 4%, and 32% of the respective national environmental burdens.

Data of gene expression levels across individuals, cell types, and disease states is expanding, yet our understanding of how expression levels impact phenotype is limited. Here, we present a massively parallel system for assaying the effect of gene expression levels on fitness in Saccharomyces cerevisiae by systematically altering the expression level of 100 genes at 100 distinct levels spanning a 500-fold range at high resolution. We show that the relationship between expression levels and growth is gene and environment specific and provides information on the function, stoichiometry, and interactions of genes. Wild-type expression levels in some conditions are not optimal for growth, and genes whose fitness is greatly affected by small changes in expression level tend to exhibit lower cell-to-cell variability in expression. Our study addresses a fundamental gap in understanding the functional significance of gene expression regulation and offers a framework for evaluating the phenotypic effects of expression variation.

Pyruvate formate-lyase (PFL) is a ubiquitous enzyme that supports increased ATP yield during sugar fermentation. While the PFL reaction is known to be reversible in vitro, the ability of PFL to support microbial growth by condensing acetyl-CoA and formate in vivo has never been directly tested. Here, we employ Escherichia coli mutant strains that cannot assimilate acetate via the glyoxylate shunt and use carbon labeling experiments to unequivocally demonstrate PFL-dependent co-assimilation of acetate and formate. Moreover, PFL-dependent growth is faster than growth on acetate using the glyoxylate shunt. Hence, growth via the reverse activity of PFL could have substantial ecological and biotechnological significance.

Turnover numbers, also known as kcat values, are fundamental properties of enzymes. However, kcat data are scarce and measured in vitro, thus may not faithfully represent the in vivo situation. A basic question that awaits elucidation is: how representative are kcat values for the maximal catalytic rates of enzymes in vivo? Here, we harness omics data to calculate kvivomax, the observed maximal catalytic rate of an enzyme inside cells. Comparison with kcat values from Escherichia coli, yields a correlation of r2= 0.62 in log scale (p

The high environmental costs of raising livestock are now widely appreciated, yet consumption of animal-based food items continues and is expanding throughout the world. Consumers' ability to distinguish among, and rank, various interchangeable animal-based items is crucial to reducing environmental costs of diets. However, the individual environmental burdens exerted by the five dominant livestock categories - beef, dairy, poultry, pork and eggs - are not fully known. Quantifying those burdens requires splitting livestock's relatively well-known total environmental costs (e.g. land and fertilizer use for feed production) into partial categorical costs. Because such partitioning quantifies the relative environmental desirability of various animal-based food items, it is essential for environmental impact minimization efforts to be made. Yet to date, no such partitioning method exists. The present paper presents such a partitioning method for feed production-related environmental burdens. This approach treated each of the main feed classes individually - concentrates (grain, soy, by-products; supporting production of all livestock), processed roughage (mostly hay and silage) and pasture - which is key given these classes' widely disparate environmental costs. It was found that for the current US food system and national diet, concentrates are partitioned as follows: beef 0·21±0·112, poultry 0·27±0·046, dairy 0·24±0·041, pork 0·23±0·093 and eggs 0·04±0·018. Pasture and processed roughage, consumed only by cattle, are 0·92±0·034 and 0·87±0·031 due to beef, with the remainder due to dairy. In a follow-up paper, the devised methodology will be employed to partition total land, irrigated water, greenhouse gases and reactive nitrogen burdens incurred by feed production among the five edible livestock categories.

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.

Most genes change expression levels across conditions, but it is unclear which of these changes represents specific regulation and what determines their quantitative degree. Here, we accurately measured activities of ∼900 S. cerevisiae and ∼1800 E. coli promoters using fluorescent reporters. We show that in both organisms 60-90% of promoters change their expression between conditions by a constant global scaling factor that depends only on the conditions and not on the promoter's identity. Quantifying such global effects allows precise characterization of specific regulation - promoters deviating from the global scale line. These are organized into few functionally related groups that also adhere to scale lines and preserve their relative activities across conditions. Thus, only several scaling factors suffice to accurately describe genome-wide expression profiles across conditions. We present a parameter-free passive resource allocation model that quantitatively accounts for the global scaling factors. It suggests that many changes in expression across conditions result from global effects and not specific regulation, and provides means for quantitative interpretation of expression profiles.

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 over-whelmingly 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.

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.

Background: Constraint-based modeling is increasingly employed for metabolic network analysis. Its underlying assumption is that natural metabolic phenotypes can be predicted by adding physicochemical constraints to remove unrealistic metabolic flux solutions. The loopless-COBRA approach provides an additional constraint that eliminates thermodynamically infeasible internal cycles (or loops) from the space of solutions. This allows the prediction of flux solutions that are more consistent with experimental data. However, it is not clear if this approach over-constrains the models by removing non-loop solutions as well.Results: Here we apply Gordan's theorem from linear algebra to prove for the first time that the constraints added in loopless-COBRA do not over-constrain the problem beyond the elimination of the loops themselves.Conclusions: The loopless-COBRA constraints can be reliably applied. Furthermore, this proof may be adapted to evaluate the theoretical soundness for other methods in constraint-based modeling.

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 nonstandard 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 pK _a 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.

Metabolic engineering of plants can reduce the cost and environmental impact of agriculture while providing for the needs of a growing population. Although our understanding of plant metabolism continues to increase at a rapid pace, relatively few plant metabolic engineering projects with commercial potential have emerged, in part because of a lack of principles for the rational manipulation of plant phenotype. One underexplored approach to identifying such design principles derives from analysis of the dominant constraints on plant fitness, and the evolutionary innovations in response to those constraints, that gave rise to the enormous diversity of natural plant metabolic pathways.

While the reductive pentose phosphate cycle is responsible for the fixation of most of the carbon in the biosphere, it has several natural substitutes. In fact, due to the characterization of three new carbon fixation pathways in the last decade, the diversity of known metabolic solutions for autotrophic growth has doubled. In this review, the different pathways are analysed and compared according to various criteria, trying to connect each of the different metabolic alternatives to suitable environments or metabolic goals. The different roles of carbon fixation are discussed; in addition to sustaining autotrophic growth it can also be used for energy conservation and as an electron sink for the recycling of reduced electron carriers. Our main focus in this review is on thermodynamic and kinetic aspects, including thermodynamically challenging reactions, the ATP requirement of each pathway, energetic constraints on carbon fixation, and factors that are expected to limit the rate of the pathways. Finally, possible metabolic structures of yet unknown carbon fixation pathways are suggested and discussed.

Latency and ongoing replication have both been proposed to explain the drug-insensitive human immunodeficiency virus (HIV) reservoir maintained during antiretroviral therapy. Here we explore a novel mechanism for ongoing HIV replication in the face of antiretroviral drugs. We propose a model whereby multiple infections per cell lead to reduced sensitivity to drugs without requiring drug-resistant mutations, and experimentally validate the model using multiple infections per cell by cell-free HIV in the presence of the drug tenofovir. We then examine the drug sensitivity of cell-to-cell spread of HIV, a mode of HIV transmission that can lead to multiple infection events per target cell. Infections originating from cell-free virus decrease strongly in the presence of antiretrovirals tenofovir and efavirenz whereas infections involving cell-to-cell spread are markedly less sensitive to the drugs. The reduction in sensitivity is sufficient to keep multiple rounds of infection from terminating in the presence of drugs. We examine replication from cell-to-cell spread in the presence of clinical drug concentrations using a stochastic infection model and find that replication is intermittent, without substantial accumulation of mutations. If cell-to-cell spread has the same properties in vivo, it may have adverse consequences for the immune system, lead to therapy failure in individuals with risk factors, and potentially contribute to viral persistence and hence be a barrier to curing HIV infection.

Background: C ₄ plants such as corn and sugarcane assimilate atmospheric CO ₂ into biomass by means of the C ₄ carbon fixation pathway. We asked how PEP formation rate, a key step in the carbon fixation pathway, might work at a precise rate, regulated by light, despite fluctuations in substrate and enzyme levels constituting and regulating this process.Results: We present a putative mechanism for robustness in C ₄ carbon fixation, involving a key enzyme in the pathway, pyruvate orthophosphate dikinase (PPDK), which is regulated by a bifunctional enzyme, Regulatory Protein (RP). The robust mechanism is based on avidity of the bifunctional enzyme RP to its multimeric substrate PPDK, and on a product-inhibition feedback loop that couples the system output to the activity of the bifunctional regulator. The model provides an explanation for several unusual biochemical characteristics of the system and predicts that the system's output, phosphoenolpyruvate (PEP) formation rate, is insensitive to fluctuations in enzyme levels (PPDK and RP), substrate levels (ATP and pyruvate) and the catalytic rate of PPDK, while remaining sensitive to the system's input (light levels).Conclusions: The presented PPDK mechanism is a new way to achieve robustness using product inhibition as a feedback loop on a bifunctional regulatory enzyme. This mechanism exhibits robustness to protein and metabolite levels as well as to catalytic rate changes. At the same time, the output of the system remains tuned to input levels.

Although the quantitative description of biological systems has been going on for centuries, recent advances in the measurement of phenomena ranging from metabolism to gene expression to signal transduction have resulted in a new emphasis on biological numeracy. This article describes the confluence of two different approaches to biological numbers. First, an impressive array of quantitative measurements make it possible to develop intuition about biological numbers ranging from how many gigatons of atmospheric carbon are fixed every year in the process of photosynthesis to the number of membrane transporters needed to provide sugars to rapidly dividing Escherichia coli cells. As a result of the vast array of such quantitative data, the BioNumbers web site has recently been developed as a repository for biology by the numbers. Second, a complementary and powerful tradition of numerical estimates familiar from the physical sciences and canonized in the so-called "Fermi problems" calls for efforts to estimate key biological quantities on the basis of a few foundational facts and simple ideas from physics and chemistry. In this article, we describe these two approaches and illustrate their synergism in several particularly appealing case studies. These case studies reveal the impact that an emphasis on numbers can have on important biological questions.

The sun's spectrum harvested through photosynthesis is the primary source of energy for life on earth. Plants, green algae, and cyanobacteria-the major primary producers on earth-utilize reaction centers that operate at wavelengths of 680 and 700 nm. Why were these wavelengths "chosen" in evolution? This study analyzes the efficiency of light conversion into chemical energy as a function of hypothetical reaction center absorption wavelengths given the sun's spectrum and the overpotential cost associated with charge separation. Surprisingly, it is found here that when taking into account the empirical charge separation cost the range 680-720 nm maximizes the conversion efficiency. This suggests the possibility that the wavelengths of photosystem I and II were optimized at some point in their evolution for the maximal utilization of the sun's spectrum.

Myelination in the peripheral nervous system requires close contact between Schwann cells and the axon, but the underlying molecular basis remains largely unknown. Here we show that cell adhesion molecules (CAMs) of the nectin-like (Necl, also known as SynCAM or Cadm) family mediate Schwann cell-axon interaction during myelination. Necl4 is the main Necl expressed by myelinating Schwann cells and is located along the internodes in direct apposition to Necl1, which is localized on axons. Necl4 serves as the glial binding partner for axonal Necl1, and the interaction between these two CAMs mediates Schwann cell adhesion. The disruption of the interaction between Necl1 and Necl4 by their soluble extracellular domains, or the expression of a dominant-negative Necl4 in Schwann cells, inhibits myelination. These results suggest that Necl proteins are important for mediating axon-glia contact during myelination in peripheral nerves.

Physiological and evolutionary adaptations operate at very different time scales. Nevertheless, there are reasons to believe there should be a strong relationship between the two, as together they modify the phenotype. Physiological adaptations change phenotype by altering certain microscopic parameters; evolutionary adaptation can either alter genetically these same parameters or others to achieve distinct or similar ends. Although qualitative discussions of this relationship abound, there has been very little quantitative analysis. Here, we use the hemoglobin molecule as a model system to quantify the relationship between physiological and evolutionary adaptations. We compare measurements of oxygen saturation curves of 25 mammals with those of human hemoglobin under a wide range of physiological conditions. We fit the data sets to the Monod-Wyman-Changeux model to extract microscopic parameters. Our analysis demonstrates that physiological and evolutionary change act on different parameters. The main parameter that changes in the physiology of hemoglobin is relatively constant in evolution, whereas the main parameter that changes in the evolution of hemoglobin is relatively constant in physiology. This orthogonality suggests continued selection for physiological adaptability and hints at a role for this adaptability in evolutionary change.

Protein expression is a stochastic process that leads to phenotypic variation among cells. The cell-cell distribution of protein levels in microorganisms has been well characterized but little is known about such variability in human cells. Here, we studied the variability of protein levels in human cells, as well as the temporal dynamics of this variability, and addressed whether cells with higher than average protein levels eventually have lower than average levels, and if so, over what timescale does this mixing occur. We measured fluctuations over time in the levels of 20 endogenous proteins in living human cells, tagged by the gene for yellow fluorescent protein at their chromosomal loci. We found variability with a standard deviation that ranged, for different proteins, from about 15% to 30% of the mean. Mixing between high and low levels occurred for all proteins, but the mixing time was longer than two cell generations (more than 40 h) for many proteins. We also tagged pairs of proteins with two colours, and found that the levels of proteins in the same biological pathway were far more correlated than those of proteins in different pathways. The persistent memory for protein levels that we found might underlie individuality in cell behaviour and could set a timescale needed for signals to affect fully every member of a cell population.

Genes and proteins generate molecular circuitry that enables the cell to process information and respond to stimuli. A major challenge is to identify characteristic patterns in this network of interactions that may shed light on basic cellular mechanisms. Previous studies have analyzed aspects of this network, concentrating on either transcription-regulation or protein-protein interactions. Here we search for composite network motifs: characteristic network patterns consisting of both transcription-regulation and protein-protein interactions that recur significantly more often than in random networks. To this end we developed algorithms for detecting motifs in networks with two or more types of interactions and applied them to an integrated data set of protein-protein interactions and transcription regulation in Saccharomyces cerevisiae. We found a two-protein mixed-feedback loop motif, five types of three-protein motifs exhibiting coregulation and complex formation, and many motifs involving four proteins. Virtually all four-protein motifs consisted of combinations of smaller motifs. This study presents a basic framework for detecting the building blocks of networks with multiple types of interactions.

Biological and technological networks contain patterns, termed network motifs, which occur far more often than in randomized networks. Network motifs were suggested to be elementary building blocks that carry out key functions in the network. It is of interest to understand how network motifs combine to form larger structures. To address this, we present a systematic approach to define \u201cmotif generalizations\u201d: families of motifs of different sizes that share a common architectural theme. To define motif generalizations, we first define \u201croles\u201d in a subgraph according to structural equivalence. For example, the feedforward loop triada motif in transcription, neuronal, and some electronic networkshas three roles: an input node, an output node, and an internal node. The roles are used to define possible generalizations of the motif. The feedforward loop can have three simple generalizations, based on replicating each of the three roles and their connections. We present algorithms for efficiently detecting motif generalizations. We find that the transcription networks of bacteria and yeast display only one of the three generalizations, the multi-output feedforward generalization. In contrast, the neuronal network of C. elegans mainly displays the multi-input generalization. Forward-logic electronic circuits display a multi-input, multi-output hybrid. Thus, networks which share a common motif can have very different generalizations of that motif. Using mathematical modeling, we describe the information processing functions of the different motif generalizations in transcription, neuronal, and electronic networks.