Quantitative Control for Stoichiometric Protein Synthesis

Cross-evaluation of E. coli's operon structures via a whole-cell model suggests alternative cellular benefits for low- versus high-expressing operons

Many bacteria use operons to coregulate genes, but it remains unclear how operons benefit bacteria. We integrated E. coli's 788 polycistronic operons and 1,231 transcription units into an existing whole-cell model and found inconsistencies between the proposed operon structures and the RNA-seq read counts that the model was parameterized from. We resolved these inconsistencies through iterative, model-guided corrections to both datasets, including the correction of RNA-seq counts of short genes that were misreported as zero by existing alignment algorithms. The resulting model suggested two main modes by which operons benefit bacteria. For 86% of low-expression operons, adding operons increased the co-expression probabilities of their constituent proteins, whereas for 92% of high-expression operons, adding operons resulted in more stable expression ratios between the proteins. These simulations underscored the need for further experimental work on how operons reduce noise and synchronize both the expression timing and the quantity of constituent genes. A record of this paper's transparent peer review process is included in the supplemental information.

A continuous epistasis model for predicting growth rate given combinatorial variation in gene expression and environment

Quantifying and predicting growth rate phenotype given variation in gene expression and environment is complicated by epistatic interactions and the vast combinatorial space of possible perturbations. We developed an approach for mapping expression-growth rate landscapes that integrates sparsely sampled experimental measurements with an interpretable machine learning model. We used mismatch CRISPRi across pairs and triples of genes to create over 8,000 titrated changes in E. coli gene expression under varied environmental contexts, exploring epistasis in up to 22 distinct environments. Our results show that a pairwise model previously used to describe drug interactions well-described these data. The model yielded interpretable parameters related to pathway architecture and generalized to predict the combined effect of up to four perturbations when trained solely on pairwise perturbation data. We anticipate this approach will be broadly applicable in optimizing bacterial growth conditions, generating pharmacogenomic models, and understanding the fundamental constraints on bacterial gene expression. A record of this paper's transparent peer review process is included in the supplemental information.

Decoding stoichiometric protein synthesis in E. coli through translation rate parameters

E. coli is one of the most widely used organisms for understanding the principles of cellular and molecular genetics. However, we are yet to understand the origin of several experimental observations related to the regulation of gene expression in E. coli. One of the prominent examples in this context is the proportional synthesis in multiprotein complexes where all of their obligate subunits are produced in proportion to their stoichiometry. In this work, by combining the next-generation sequencing data with the stochastic simulations of protein synthesis, we explain the origin of proportional protein synthesis in multicomponent complexes. We find that the estimated initiation rates for the translation of all subunits in those complexes are proportional to their stoichiometry. This constraint on protein synthesis kinetics enforces proportional protein synthesis without requiring any feedback mechanism. We also find that the translation initiation rates in E. coli are influenced by the coding sequence length and the enrichment of A and C nucleotides near the start codon. Thus, this study rationalizes the role of conserved and nonrandom features of genes in regulating the translation kinetics and unravels a key principle of the regulation of protein synthesis.

Selective RNA Processing and Stabilization are Multi-Layer and Stoichiometric Regulators of Gene Expression in Escherichia coli

Selective RNA processing and stabilization (SRPS) facilitates the differential expression of multiple genes in polycistronic operons. However, how the coordinated actions of SRPS-related enzymes affect stoichiometric regulation remains unclear. In the present study, the first genome-wide targetome analysis is reported of these enzymes in Escherichia coli, at a single-nucleotide resolution. A strictly linear relationship is observed between the RNA pyrophosphohydrolase processing ratio and scores assigned to the first three nucleotides of the primary transcript. Stem-loops associated with PNPase targetomes exhibit a folding free energy that is negatively correlated with the termination ratio of PNPase at the 3' end. More than one-tenth of the RNase E processing sites in the 5'-untranslated regions（UTR） form different stem-loops that affect ribosome-binding and translation efficiency. The effectiveness of the SRPS elements is validated using a dual-fluorescence reporter system. The findings highlight a multi-layer and quantitative regulatory method for optimizing the stoichiometric expression of genes in bacteria and promoting the application of SRPS in synthetic biology.

Regulatory response to a hybrid ancestral nitrogenase in Azotobacter vinelandii

Biological nitrogen fixation, the microbial reduction of atmospheric nitrogen to bioavailable ammonia, represents both a major limitation on biological productivity and a highly desirable engineering target for synthetic biology. However, the engineering of nitrogen fixation requires an integrated understanding of how the gene regulatory dynamics of host diazotrophs respond across sequence-function space of its central catalytic metalloenzyme, nitrogenase. Here, we interrogate this relationship by analyzing the transcriptome of Azotobacter vinelandii engineered with a phylogenetically inferred ancestral nitrogenase protein variant. The engineered strain exhibits reduced cellular nitrogenase activity but recovers wild-type growth rates following an extended lag period. We find that expression of genes within the immediate nitrogen fixation network is resilient to the introduced nitrogenase sequence-level perturbations. Rather the sustained physiological compatibility with the ancestral nitrogenase variant is accompanied by reduced expression of genes that support trace metal and electron resource allocation to nitrogenase. Our results spotlight gene expression changes in cellular processes adjacent to nitrogen fixation as productive engineering considerations to improve compatibility between remodeled nitrogenase proteins and engineered host diazotrophs. IMPORTANCE Azotobacter vinelandii is a key model bacterium for the study of biological nitrogen fixation, an important metabolic process catalyzed by nitrogenase enzymes. Here, we demonstrate that compatibilities between engineered A. vinelandii strains and nitrogenase variants can be modulated at the regulatory level. The engineered strain studied here responds by adjusting the expression of proteins involved in cellular processes adjacent to nitrogen fixation, rather than that of nitrogenase proteins themselves. These insights can inform future strategies to transfer nitrogenase variants to non-native hosts.

Precise programming of multigene expression stoichiometry in mammalian cells by a modular and programmable transcriptional system

Context-dependency of mammalian transcriptional elements has hindered the quantitative investigation of multigene expression stoichiometry and its biological functions. Here, we describe a host- and local DNA context-independent transcription system to gradually fine-tune single and multiple gene expression with predictable stoichiometries. The mammalian transcription system is composed of a library of modular and programmable promoters from bacteriophage and its cognate RNA polymerase (RNAP) fused to a capping enzyme. The relative expression of single genes is quantitatively determined by the relative binding affinity of the RNAP to the promoters, while multigene expression stoichiometry is predicted by a simple biochemical model with resource competition. We use these programmable and modular promoters to predictably tune the expression of three components of an influenza A virus-like particle (VLP). Optimized stoichiometry leads to a 2-fold yield of intact VLP complexes. The host-independent orthogonal transcription system provides a platform for dose-dependent control of multiple protein expression which may be applied for advanced vaccine engineering, cell-fate programming and other therapeutic applications.

Toward the Complete Functional Characterization of a Minimal Bacterial Proteome

Recently, we presented a whole-cell kinetic model of the genetically minimal bacterium JCVI-syn3A that described the coupled metabolic and genetic information processes and predicted behaviors emerging from the interactions among these networks. JCVI-syn3A is a genetically reduced bacterial cell that has the fewest number and smallest fraction of genes of unclear function, with approximately 90 of its 452 protein-coding genes (that is less than 20%) unannotated. Further characterization of unclear JCVI-syn3A genes strengthens the robustness and predictive power of cell modeling efforts and can lead to a deeper understanding of biophysical processes and pathways at the cell scale. Here, we apply computational analyses to elucidate the functions of the products of several essential but previously uncharacterized genes involved in integral cellular processes, particularly those directly affecting cell growth, division, and morphology. We also suggest directed wet-lab experiments informed by our analyses to further understand these "missing puzzle pieces" that are an essential part of the mosaic of biological interactions present in JCVI-syn3A. Our workflow leverages evolutionary sequence analysis, protein structure prediction, interactomics, and genome architecture to determine upgraded annotations. Additionally, we apply the structure prediction analysis component of our work to all 452 protein coding genes in JCVI-syn3A to expedite future functional annotation studies as well as the inverse mapping of the cell state to more physical models requiring all-atom or coarse-grained representations for all JCVI-syn3A proteins.

Massively parallel characterization of engineered transcript isoforms using direct RNA sequencing

Transcriptional terminators signal where transcribing RNA polymerases (RNAPs) should halt and disassociate from DNA. However, because termination is stochastic, two different forms of transcript could be produced: one ending at the terminator and the other reading through. An ability to control the abundance of these transcript isoforms would offer bioengineers a mechanism to regulate multi-gene constructs at the level of transcription. Here, we explore this possibility by repurposing terminators as 'transcriptional valves' that can tune the proportion of RNAP read-through. Using one-pot combinatorial DNA assembly, we iteratively construct 1780 transcriptional valves for T7 RNAP and show how nanopore-based direct RNA sequencing (dRNA-seq) can be used to characterize entire libraries of valves simultaneously at a nucleotide resolution in vitro and unravel genetic design principles to tune and insulate termination. Finally, we engineer valves for multiplexed regulation of CRISPR guide RNAs. This work provides new avenues for controlling transcription and demonstrates the benefits of long-read sequencing for exploring complex sequence-function landscapes.

Regulation of gene expression via translational buffering

Translation of an mRNA represents a critical step during the expression of protein-coding genes. As mechanisms governing post-transcriptional regulation of gene expression are progressively unveiled, it is becoming apparent that transcriptional programs are not fully reflected in the proteome. Herein, we highlight a previously underappreciated post-transcriptional mode of regulation of gene expression termed translational buffering. In principle, translational buffering opposes the impact of alterations in mRNA levels on the proteome. We further describe three types of translational buffering: compensation, which maintains protein levels e.g. across species or individuals; equilibration, which retains pathway stoichiometry; and offsetting, which acts as a reversible mechanism that maintains the levels of selected subsets of proteins constant despite genetic alteration and/or stress-induced changes in corresponding mRNA levels. While mechanisms underlying compensation and equilibration have been reviewed elsewhere, the principal focus of this review is on the less-well understood mechanism of translational offsetting. Finally, we discuss potential roles of translational buffering in homeostasis and disease.