The genetic basis for the adaptation of E. coli to sugar synthesis from CO2.

Understanding the evolution of a new metabolic capability in full mechanistic detail is challenging, as causative mutations may be masked by non-essential "hitchhiking" mutations accumulated during the evolutionary trajectory. We have previously used adaptive laboratory evolution of a rationally engineered ancestor to generate an Escherichia coli strain able to utilize CO2 fixation for sugar synthesis. Here, we reveal the genetic basis underlying this metabolic transition. Five mutations are sufficient to enable robust growth when a non-native Calvin-Benson-Bassham cycle provides all the sugar-derived metabolic building blocks. These mutations are found either in enzymes that affect the efflux of intermediates from the autocatalytic CO2 fixation cycle toward biomass (prs, serA, and pgi), or in key regulators of carbon metabolism (crp and ppsR). Using suppressor analysis, we show that a decrease in catalytic capacity is a common feature of all mutations found in enzymes. These findings highlight the enzymatic constraints that are essential to the metabolic stability of autocatalytic cycles and are relevant to future efforts in constructing non-native carbon fixation pathways.

Design principles of autocatalytic cycles constrain enzyme kinetics and force low substrate saturation at flux branch points.

A set of chemical reactions that require a metabolite to synthesize more of that metabolite is an autocatalytic cycle. Here, we show that most of the reactions in the core of central carbon metabolism are part of compact autocatalytic cycles. Such metabolic designs must meet specific conditions to support stable fluxes, hence avoiding depletion of intermediate metabolites. As such, they are subjected to constraints that may seem counter-intuitive: the enzymes of branch reactions out of the cycle must be overexpressed and the affinity of these enzymes to their substrates must be relatively weak. We use recent quantitative proteomics and fluxomics measurements to show that the above conditions hold for functioning cycles in central carbon metabolism of E. coli. This work demonstrates that the topology of a metabolic network can shape kinetic parameters of enzymes and lead to seemingly wasteful enzyme usage.

Sugar Synthesis from CO2 in Escherichia coli.

Can a heterotrophic organism be evolved to synthesize biomass from CO2 directly? So far, non-native carbon fixation in which biomass precursors are synthesized solely from CO2 has remained an elusive grand challenge. Here, we demonstrate how a combination of rational metabolic rewiring, recombinant expression, and laboratory evolution has led to the biosynthesis of sugars and other major biomass constituents by a fully functional Calvin-Benson-Bassham (CBB) cycle in E. coli. In the evolved bacteria, carbon fixation is performed via a non-native CBB cycle, while reducing power and energy are obtained by oxidizing a supplied organic compound (e.g., pyruvate). Genome sequencing reveals that mutations in flux branchpoints, connecting the non-native CBB cycle to biosynthetic pathways, are essential for this phenotype. The successful evolution of a non-native carbon fixation pathway, though not yet resulting in net carbon gain, strikingly demonstrates the capacity for rapid trophic-mode evolution of metabolism applicable to biotechnology. PAPERCLIP.

A Minimalistic Resource Allocation Model to Explain Ubiquitous Increase in Protein Expression with Growth Rate.

Most proteins show changes in level across growth conditions. Many of these changes seem to be coordinated with the specific growth rate rather than the growth environment or the protein function. Although cellular growth rates, gene expression levels and gene regulation have been at the center of biological research for decades, there are only a few models giving a base line prediction of the dependence of the proteome fraction occupied by a gene with the specific growth rate. We present a simple model that predicts a widely coordinated increase in the fraction of many proteins out of the proteome, proportionally with the growth rate. The model reveals how passive redistribution of resources, due to active regulation of only a few proteins, can have proteome wide effects that are quantitatively predictable. Our model provides a potential explanation for why and how such a coordinated response of a large fraction of the proteome to the specific growth rate arises under different environmental conditions. The simplicity of our model can also be useful by serving as a baseline null hypothesis in the search for active regulation. We exemplify the usage of the model by analyzing the relationship between growth rate and proteome composition for the model microorganism E.coli as reflected in recent proteomics data sets spanning various growth conditions. We find that the fraction out of the proteome of a large number of proteins, and from different cellular processes, increases proportionally with the growth rate. Notably, ribosomal proteins, which have been previously reported to increase in fraction with growth rate, are only a small part of this group of proteins. We suggest that, although the fractions of many proteins change with the growth rate, such changes may be partially driven by a global effect, not necessarily requiring specific cellular control mechanisms.

Global characterization of in vivo enzyme catalytic rates and their correspondence to in vitro kcat measurements.

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 kmax(vivo), the observed maximal catalytic rate of an enzyme inside cells. Comparison with kcat values from Escherichia coli, yields a correlation ofr(2)= 0.62 in log scale (p < 10(-10)), with a root mean square difference of 0.54 (3.5-fold in linear scale), indicating that in vivo and in vitro maximal rates generally concur. By accounting for the degree of saturation of enzymes and the backward flux dictated by thermodynamics, we further refine the correspondence between kmax(vivo) and kcat values. The approach we present here characterizes the quantitative relationship between enzymatic catalysis in vitro and in vivo and offers a high-throughput method for extracting enzyme kinetic constants from omics data.

Promoters maintain their relative activity levels under different growth conditions.

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.

Quantifying translational coupling in E. coli synthetic operons using RBS modulation and fluorescent reporters.

Translational coupling is the interdependence of translation efficiency of neighboring genes encoded within an operon. The degree of coupling may be quantified by measuring how the translation rate of a gene is modulated by the translation rate of its upstream gene. Translational coupling was observed in prokaryotic operons several decades ago, but the quantitative range of modulation translational coupling leads to and the factors governing this modulation were only partially characterized. In this study, we systematically quantify and characterize translational coupling in E. coli synthetic operons using a library of plasmids carrying fluorescent reporter genes that are controlled by a set of different ribosome binding site (RBS) sequences. The downstream gene expression level is found to be enhanced by the upstream gene expression via translational coupling with the enhancement level varying from almost no coupling to over 10-fold depending on the upstream gene's sequence. Additionally, we find that the level of translational coupling in our system is similar between the second and third locations in the operon. The coupling depends on the distance between the stop codon of the upstream gene and the start codon of the downstream gene. This study is the first to systematically and quantitatively characterize translational coupling in a synthetic E. coli operon. Our analysis will be useful in accurate manipulation of gene expression in synthetic biology and serves as a step toward understanding the mechanisms involved in translational expression modulation.

Spanning high-dimensional expression space using ribosome-binding site combinatorics.

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.