Modelling and measuring intracellular competition for finite resources during gene expression

Cooperative insulation of regulatory domains by CTCF-dependent physical insulation and promoter competition

The specificity of gene expression during development requires the insulation of regulatory domains to avoid inappropriate enhancer-gene interactions. In vertebrates, this insulator function is mostly attributed to clusters of CTCF sites located at topologically associating domain (TAD) boundaries. However, TAD boundaries allow some physical crosstalk across regulatory domains, which is at odds with the specific and precise expression of developmental genes. Here we show that developmental genes and nearby clusters of CTCF sites cooperatively foster the robust insulation of regulatory domains. By genetically dissecting a couple of representative loci in mouse embryonic stem cells, we show that CTCF sites prevent undesirable enhancer-gene contacts (i.e. physical insulation), while developmental genes preferentially contribute to regulatory insulation through non-structural mechanisms involving promoter competition rather than enhancer blocking. Overall, our work provides important insights into the insulation of regulatory domains, which in turn might help interpreting the pathological consequences of certain structural variants.

Optimal control of ribosome population for gene expression under periodic nutrient intake

Translation of proteins is a fundamental part of gene expression that is mediated by ribosomes. As ribosomes significantly contribute to both cellular mass and energy consumption, achieving efficient management of the ribosome population is also crucial to metabolism and growth. Inspired by biological evidence for nutrient-dependent mechanisms that control both ribosome-active degradation and genesis, we introduce a dynamical model of protein production, that includes the dynamics of resources and control over the ribosome population. Under the hypothesis that active degradation and biogenesis are optimal for maximizing and maintaining protein production, we aim to qualitatively reproduce empirical observations of the ribosome population dynamics. Upon formulating the associated optimization problem, we first analytically study the stability and global behaviour of solutions under constant resource input, and characterize the extent of oscillations and convergence rate to a global equilibrium. We further use these results to simplify and solve the problem under a quasi-static approximation. Using biophysical parameter values, we find that optimal control solutions lead to both control mechanisms and the ribosome population switching between periods of feeding and fasting, suggesting that the intense regulation of ribosome population observed in experiments allows to maximize and maintain protein production. Finally, we find some range for the control values over which such a regime can be observed, depending on the intensity of fasting.

The miRNA-target interactions: An underestimated intricacy

MicroRNAs (miRNAs) play indispensable roles in posttranscriptional gene regulation. Their cellular regulatory impact is determined not solely by their sheer number, which likely amounts to >2000 individual miRNAs in human, than by the regulatory effectiveness of single miRNAs. Although, one begins to develop an understanding of the complex mechanisms underlying miRNA-target interactions (MTIs), the overall knowledge of MTI functionality is still rather patchy. In this critical review, we summarize key features of mammalian MTIs. We especially highlight latest insights on (i) the dynamic make-up of miRNA binding sites including non-canonical binding sites, (ii) the cooperativity between miRNA binding sites, (iii) the adaptivity of MTIs through sequence modifications, (iv) the bearing of intra-cellular miRNA localization changes and (v) the role of cell type and cell status specific miRNA interaction partners. The MTI biology is discussed against the background of state-of-the-art approaches with particular emphasis on experimental strategies for evaluating miRNA functionality.

Emergence of division of labor in tissues through cell interactions and spatial cues

Most cell types in multicellular organisms can perform multiple functions. However, not all functions can be optimally performed simultaneously by the same cells. Functions incompatible at the level of individual cells can be performed at the cell population level, where cells divide labor and specialize in different functions. Division of labor can arise due to instruction by tissue environment or through self-organization. Here, we develop a computational framework to investigate the contribution of these mechanisms to division of labor within a cell-type population. By optimizing collective cellular task performance under trade-offs, we find that distinguishable expression patterns can emerge from cell-cell interactions versus instructive signals. We propose a method to construct ligand-receptor networks between specialist cells and use it to infer division-of-labor mechanisms from single-cell RNA sequencing (RNA-seq) and spatial transcriptomics data of stromal, epithelial, and immune cells. Our framework can be used to characterize the complexity of cell interactions within tissues.

Translation in the cell under fierce competition for shared resources: a mathematical model

During translation, mRNAs 'compete' for shared resources. Under stress conditions, during viral infection and also in high-throughput heterologous gene expression, these resources may become scarce, e.g. the pool of free ribosomes is starved, and then the competition may have a dramatic effect on the global dynamics of translation in the cell. We model this scenario using a network that includes m ribosome flow models (RFMs) interconnected via a pool of free ribosomes. Each RFM models ribosome flow along an mRNA molecule, and the pool models the shared resource. We assume that the number of mRNAs is large, so many ribosomes are attached to the mRNAs, and the pool is starved. Our analysis shows that adding an mRNA has an intricate effect on the total protein production. The new mRNA produces new proteins, but the other mRNAs produce less proteins, as the pool that feeds these mRNAs now has a smaller abundance of ribosomes. As the number of mRNAs increases, the marginal utility of adding another mRNA diminishes, and the total protein production rate saturates to a limiting value. We demonstrate our approach using an example of insulin protein production in a cell-free system.

Predicting Composition of Genetic Circuits with Resource Competition: Demand and Sensitivity

The design of genetic circuits typically relies on characterization of constituent modules in isolation to predict the behavior of modules' composition. However, it has been shown that the behavior of a genetic module changes when other modules are in the cell due to competition for shared resources. In order to engineer multimodule circuits that behave as intended, it is thus necessary to predict changes in the behavior of a genetic module when other modules load cellular resources. Here, we introduce two characteristics of circuit modules: the demand for cellular resources and the sensitivity to resource loading. When both are known for every genetic module in a circuit library, they can be used to predict any module's behavior upon addition of any other module to the cell. We develop an experimental approach to measure both characteristics for any circuit module using a resource sensor module. Using the measured resource demand and sensitivity for each module in a library, the outputs of the modules can be accurately predicted when they are inserted in the cell in arbitrary combinations. These resource competition characteristics may be used to inform the design of genetic circuits that perform as predicted despite resource competition.

Accelerating Whole-Cell Simulations of mRNA Translation Using a Dedicated Hardware

In recent years, intracellular biophysical simulations have been used with increasing frequency not only for answering basic scientific questions but also in the field of synthetic biology. However, since these models include networks of interaction between millions of components, they are extremely time-consuming and cannot run easily on parallel computers. In this study, we demonstrate for the first time a novel approach addressing this challenge by using a dedicated hardware designed specifically to simulate such processes. As a proof of concept, we specifically focus on mRNA translation, which is the process consuming most of the energy in the cell. We design a hardware that simulates translation in Escherichia coli and Saccharomyces cerevisiae for thousands of mRNAs and ribosomes, which is in orders of magnitude faster than a similar software solution. With the sharp increase in the amount of genomic data available today and the complexity of the corresponding models inferred from them, we believe that the strategy suggested here will become common and can be used among others for simulating entire cells with all gene expression steps.

RBS and Promoter Strengths Determine the Cell-Growth-Dependent Protein Mass Fractions and Their Optimal Synthesis Rates

Models of gene expression considering host–circuit interactions are relevant for understanding both the strategies and associated trade-offs that cell endogenous genes have evolved and for the efficient design of heterologous protein expression systems and synthetic genetic circuits. Here, we consider a small-size model of gene expression dynamics in bacterial cells accounting for host–circuit interactions due to limited cellular resources. We define the cellular resources recruitment strength as a key functional coefficient that explains the distribution of resources among the host and the genes of interest and the relationship between the usage of resources and cell growth. This functional coefficient explicitly takes into account lab-accessible gene expression characteristics, such as promoter and ribosome binding site (RBS) strengths, capturing their interplay with the growth-dependent flux of available free cell resources. Despite its simplicity, the model captures the differential role of promoter and RBS strengths in the distribution of protein mass fractions as a function of growth rate and the optimal protein synthesis rate with remarkable fit to the experimental data from the literature for Escherichia coli. This allows us to explain why endogenous genes have evolved different strategies in the expression space and also makes the model suitable for model-based design of exogenous synthetic gene expression systems with desired characteristics.

Gene Expression Space Shapes the Bioprocess Trade-Offs among Titer, Yield and Productivity

Optimal gene expression is central for the development of both bacterial expression systems for heterologous protein production, and microbial cell factories for industrial metabolite production. Our goal is to fulfill industry-level overproduction demands optimally, as measured by the following key performance metrics: titer, productivity rate, and yield (TRY). Here we use a multiscale model incorporating the dynamics of (i) the cell population in the bioreactor, (ii) the substrate uptake and (iii) the interaction between the cell host and expression of the protein of interest. Our model predicts cell growth rate and cell mass distribution between enzymes of interest and host enzymes as a function of substrate uptake and the following main lab-accessible gene expression-related characteristics: promoter strength, gene copy number and ribosome binding site strength. We evaluated the differential roles of gene transcription and translation in shaping TRY trade-offs for a wide range of expression levels and the sensitivity of the TRY space to variations in substrate availability. Our results show that, at low expression levels, gene transcription mainly defined TRY, and gene translation had a limited effect; whereas, at high expression levels, TRY depended on the product of both, in agreement with experiments in the literature.

Context-aware synthetic biology by controller design: Engineering the mammalian cell

The rise of systems biology has ushered a new paradigm: the view of the cell as a system that processes environmental inputs to drive phenotypic outputs. Synthetic biology provides a complementary approach, allowing us to program cell behavior through the addition of synthetic genetic devices into the cellular processor. These devices, and the complex genetic circuits they compose, are engineered using a design-prototype-test cycle, allowing for predictable device performance to be achieved in a context-dependent manner. Within mammalian cells, context effects impact device performance at multiple scales, including the genetic, cellular, and extracellular levels. In order for synthetic genetic devices to achieve predictable behaviors, approaches to overcome context dependence are necessary. Here, we describe control systems approaches for achieving context-aware devices that are robust to context effects. We then consider cell fate programing as a case study to explore the potential impact of context-aware devices for regenerative medicine applications.

Molecules to Microbes

How did life begin on Earth? And is there life elsewhere in the Cosmos? Challenging questions, indeed. The series of conferences established by NoR CEL in 2013 addresses these very questions. This paper comprises a summary report of oral presentations that were delivered by NoR CEL’s network members during the 2018 Athens conference and, as such, disseminates the latest research which they have put forward. More in depth material can be found by consulting the contributors referenced papers. Overall, the outcome of this conspectus on the conference demonstrates a case for the existence of “probable chemistry” during the prebiotic epoch.

An endoribonuclease-based feedforward controller for decoupling resource-limited genetic modules in mammalian cells

Synthetic biology has the potential to bring forth advanced genetic devices for applications in healthcare and biotechnology. However, accurately predicting the behavior of engineered genetic devices remains difficult due to lack of modularity, wherein a device's output does not depend only on its intended inputs but also on its context. One contributor to lack of modularity is loading of transcriptional and translational resources, which can induce coupling among otherwise independently-regulated genes. Here, we quantify the effects of resource loading in engineered mammalian genetic systems and develop an endoribonuclease-based feedforward controller that can adapt the expression level of a gene of interest to significant resource loading in mammalian cells. Near-perfect adaptation to resource loads is facilitated by high production and catalytic rates of the endoribonuclease. Our design is portable across cell lines and enables predictable tuning of controller function. Ultimately, our controller is a general-purpose device for predictable, robust, and context-independent control of gene expression.

Computational discovery and modeling of novel gene expression rules encoded in the mRNA

The transcript is populated with numerous overlapping codes that regulate all steps of gene expression. Deciphering these codes is very challenging due to the large number of variables involved, the non-modular nature of the codes, biases and limitations in current experimental approaches, our limited knowledge in gene expression regulation across the tree of life, and other factors. In recent years, it has been shown that computational modeling and algorithms can significantly accelerate the discovery of novel gene expression codes. Here, we briefly summarize the latest developments and different approaches in the field.

Whole cell biophysical modeling of codon-tRNA competition reveals novel insights related to translation dynamics

The importance of mRNA translation models has been demonstrated across many fields of science and biotechnology. However, a whole cell model with codon resolution and biophysical dynamics is still lacking. We describe a whole cell model of translation for E. coli. The model simulates all major translation components in the cell: ribosomes, mRNAs and tRNAs. It also includes, for the first time, fundamental aspects of translation, such as competition for ribosomes and tRNAs at a codon resolution while considering tRNAs wobble interactions and tRNA recycling. The model uses parameters that are tightly inferred from large scale measurements of translation. Furthermore, we demonstrate a robust modelling approach which relies on state-of-the-art practices of translation modelling and also provides a framework for easy generalizations. This novel approach allows simulation of thousands of mRNAs that undergo translation in the same cell with common resources such as ribosomes and tRNAs in feasible time. Based on this model, we demonstrate, for the first time, the direct importance of competition for resources on translation and its accurate modelling. An effective supply-demand ratio (ESDR) measure, which is related to translation factors such as tRNAs, has been devised and utilized to show superior predictive power in complex scenarios of heterologous gene expression. The devised model is not only more accurate than the existing models, but, more importantly, provides a framework for analyzing complex whole cell translation problems and variables that haven't been explored before, making it important in various biomedical fields.

mRNA translation is a fundamental process in all living organisms and the importance of its modeling has been demonstrated across many fields of science and biotechnology. Specifically, modeling a whole cell context with a high resolution has been a great challenge in the field, making many important problems un-addressable. In this study we devised a novel model, which allows, for the first time, simultaneous simulation of thousands of mRNAs, along with various bio-physical aspects that affect translation (such as codon-resolution dynamics and shared resources pool of both ribosomes and tRNAs). We demonstrated (using experimental data) that this model is more accurate than existing ones, and, more importantly, provides a framework for addressing complex translation problems (such as heterologous expression) at whole cell scale and in reasonable time. We demonstrated the model using E. coli data, but the model can be easily tailored to other organisms as well. Our model addresses an urgent unmet need for biophysically accurate whole cell translation model with resources coupling and has potential applications in many fields, including medicine and biotechnology.

Widespread non-modular overlapping codes in the coding regions

Messenger RNAs (mRNAs) consist of a coding region (open reading frame (ORF)) and two untranslated regions (UTRs), 5'UTR and 3'UTR. Ribosomes travel along the coding region, translating nucleotide triplets (called codons) to a chain of amino acids. The coding region was long believed to mainly encode the amino acid content of proteins, whereas regulatory signals reside in the UTRs and in other genomic regions. However, in recent years we have learned that the ORF is expansively populated with various regulatory signals, or codes, which are related to all gene expression steps and additional intracellular aspects. In this paper, we review the current knowledge related to overlapping codes inside the coding regions, such as the influence of synonymous codon usage on translation speed (and, in turn, the effect of translation speed on protein folding), ribosomal frameshifting, mRNA stability, methylation, splicing, transcription and more. All these codes come together and overlap in the ORF sequence, ensuring production of the right protein at the right time.

Molecules to Microbes

How did life begin on Earth? And is there life elsewhere in the Cosmos? Challenging questions, indeed. The series of conferences established by NoR CEL in 2013, addresses these very same questions. The basis for this paper is the summary report of oral presentations that were delivered by NoR CEL’s network members during the 2018 Athens conference and, as such, disseminates the latest research which they have put forward. More in depth material can be found by consulting the contributors referenced papers. Overall, the outcome of this conspectus on the conference demonstrates a case for the existence of “probable chemistry” during the prebiotic epoch.

Heterogeneity coordinates bacterial multi-gene expression in single cells

For a genetically identical microbial population, multi-gene expression in various environments requires effective allocation of limited resources and precise control of heterogeneity among individual cells. However, it is unclear how resource allocation and cell-to-cell variation jointly shape the overall performance. Here we demonstrate a Simpson’s paradox during overexpression of multiple genes: two competing proteins in single cells correlated positively for every induction condition, but the overall correlation was negative. Yet this phenomenon was not observed between two competing mRNAs in single cells. Our analytical framework shows that the phenomenon arises from competition for translational resource, with the correlation modulated by both mRNA and ribosome variability. Thus, heterogeneity plays a key role in single-cell multi-gene expression and provides the population with an evolutionary advantage, as demonstrated in this study.

Microbes perform multitasking for a wide range of purposes, including survival, adaptation, colonization, and evolution. Both modelling and experimental results at the ensemble level reveal trade-offs between different tasks due to resource competition, but it is unclear how single cells allocate limited intracellular resources to perform multitasking, and how does a population coordinate single cell performances during multitasking to maximize population efficiencies. In this study, we address this question by using bacterial multi-gene overexpression as the basic form of multitasking. We discovered and analyzed a statistical phenomenon called Simpson’s paradox, where competing proteins in single cells correlate positively at each constant condition, although the proteins correlate negatively when all conditions are combined. We demonstrate that the phenomenon arises from competition for translational resources, with the correlation modulated by heterogeneity of both mRNA and ribosomes. We further show that heterogeneity coordinates multiple functional modules, conferring an evolutionary advantage on the population. Our work discloses that heterogeneity in the form of Simpson’s paradox is an important phenomenon in coordinating multi-gene expression.

An equilibrium model for ribosome competition

The number of ribosomes in a cell is considered as limiting, and gene expression is thus largely determined by their cellular concentration. In this work we develop a toy model to study the trade-off between the ribosomal supply and the demand of the translation machinery, dictated by the composition of the transcript pool. Our equilibrium framework is useful to highlight qualitative behaviours and new means of gene expression regulation determined by the fine balance of this trade-off. We also speculate on the possible impact of these mechanisms on cellular physiology.

Contemporary Tools for Regulating Gene Expression in Bacteria

Insights from novel mechanistic paradigms in gene expression control have led to the development of new gene expression systems for bioproduction, control, and sensing applications. Coupled with a greater understanding of synthetic burden and modern creative biodesign approaches, contemporary bacterial gene expression tools and systems are emerging that permit fine-tuning of expression, enabling greater predictability and maximisation of specific productivity, while minimising deleterious effects upon cell viability. These advances have been achieved by using a plethora of regulatory tools, operating at all levels of the so-called 'central dogma' of molecular biology. In this review, we discuss these gene regulation tools in the context of their design, prototyping, integration into expression systems, and biotechnological application.

Molecules to Microbes

How did life begin on Earth? And is there life elsewhere in the Cosmos? Challenging questions, indeed. The series of conferences established by NoR CEL in 2013, addresses these very same questions. The basis for this paper is the summary report of oral presentations that were delivered by NoR CEL’s network members during the 2018 Athens conference and, as such, disseminates the latest research which they have put forward. More in depth material can be found by consulting the contributors referenced papers. Overall, the outcome of this conspectus on the conference demonstrates a case for the existence of “probable chemistry” during the prebiotic epoch.