MIRELLA: a mathematical model explains the effect of microRNA-mediated synthetic genes regulation on intracellular resource allocation.

Competition for intracellular resources, also known as gene expression burden, induces coupling between independently co-expressed genes, a detrimental effect on predictability and reliability of gene circuits in mammalian cells. We recently showed that microRNA (miRNA)-mediated target downregulation correlates with the upregulation of a co-expressed gene, and by exploiting miRNAs-based incoherent-feed-forward loops (iFFLs) we stabilise a gene of interest against burden. Considering these findings, we speculate that miRNA-mediated gene downregulation causes cellular resource redistribution. Despite the extensive use of miRNA in synthetic circuits regulation, this indirect effect was never reported before. Here we developed a synthetic genetic system that embeds miRNA regulation, and a mathematical model, MIRELLA, to unravel the miRNA (MI) RolE on intracellular resource aLLocAtion. We report that the link between miRNA-gene downregulation and independent genes upregulation is a result of the concerted action of ribosome redistribution and 'queueing-effect' on the RNA degradation pathway. Taken together, our results provide for the first time insights into the hidden regulatory interaction of miRNA-based synthetic networks, potentially relevant also in endogenous gene regulation. Our observations allow to define rules for complexity- and context-aware design of genetic circuits, in which transgenes co-expression can be modulated by tuning resource availability via number and location of miRNA target sites.

Distinct retrograde microtubule motor sets drive early and late endosome transport.

Although subcellular positioning of endosomes significantly impacts on their functions, the molecular mechanisms governing the different steady-state distribution of early endosomes (EEs) and late endosomes (LEs)/lysosomes (LYs) in peripheral and perinuclear eukaryotic cell areas, respectively, are still unsolved. We unveil that such differences arise because, while LE retrograde transport depends on the dynein microtubule (MT) motor only, the one of EEs requires the cooperative antagonism of dynein and kinesin-14 KIFC1, a MT minus end-directed motor involved in cancer progression. Mechanistically, the Ser-x-Ile-Pro (SxIP) motif-mediated interaction of the endoplasmic reticulum transmembrane protein stromal interaction molecule 1 (STIM1) with the MT plus end-binding protein 1 (EB1) promotes its association with the p150Glued subunit of the dynein activator complex dynactin and the distinct location of EEs and LEs/LYs. The peripheral distribution of EEs requires their p150Glued-mediated simultaneous engagement with dynein and SxIP motif-containing KIFC1, via HOOK1 and HOOK3 adaptors, respectively. In sum, we provide evidence that distinct minus end-directed MT motor systems drive the differential transport and subcellular distribution of EEs and LEs in mammalian cells.

MicroRNAs organize intrinsic variation into stem cell states.

Pluripotent embryonic stem cells (ESCs) contain the potential to form a diverse array of cells with distinct gene expression states, namely the cells of the adult vertebrate. Classically, diversity has been attributed to cells sensing their position with respect to external morphogen gradients. However, an alternative is that diversity arises in part from cooption of fluctuations in the gene regulatory network. Here we find ESCs exhibit intrinsic heterogeneity in the absence of external gradients by forming interconverting cell states. States vary in developmental gene expression programs and display distinct activity of microRNAs (miRNAs). Notably, miRNAs act on neighborhoods of pluripotency genes to increase variation of target genes and cell states. Loss of miRNAs that vary across states reduces target variation and delays state transitions, suggesting variable miRNAs organize and propagate variation to promote state transitions. Together these findings provide insight into how a gene regulatory network can coopt variation intrinsic to cell systems to form robust gene expression states. Interactions between intrinsic heterogeneity and environmental signals may help achieve developmental outcomes.

Rab11 activity and PtdIns(3)P turnover removes recycling cargo from endosomes.

Directional transport of recycling cargo from early endosomes (EE) to the endocytic recycling compartment (ERC) relies on phosphatidylinositol 3-phosphate (PtdIns(3)P) hydrolysis and activation of the small GTPase Rab11. However, how these events are coordinated is yet unclear. By using a novel genetically-encoded FRET biosensor for Rab11, we report that generation of endosomal PtdIns(3)P by the clathrin-binding phosphoinositide 3-kinase class 2 alpha (PI3K-C2α) controls the activation of Rab11. Active Rab11, in turn, prompts the recruitment of the phosphatidylinositol 3-phosphatase myotubularin 1 (MTM1), eventually enabling the release of recycling cargo from the EE and its delivery toward the ERC. Our findings thus define that delivery of recycling cargo toward the ERC requires spatial and sequential coupling of Rab11 activity with PtdIns(3)P turnover.

On the role of extrinsic noise in microRNA-mediated bimodal gene expression.

Several studies highlighted the relevance of extrinsic noise in shaping cell decision making and differentiation in molecular networks. Bimodal distributions of gene expression levels provide experimental evidence of phenotypic differentiation, where the modes of the distribution often correspond to different physiological states of the system. We theoretically address the presence of bimodal phenotypes in the context of microRNA (miRNA)-mediated regulation. MiRNAs are small noncoding RNA molecules that downregulate the expression of their target mRNAs. The nature of this interaction is titrative and induces a threshold effect: below a given target transcription rate almost no mRNAs are free and available for translation. We investigate the effect of extrinsic noise on the system by introducing a fluctuating miRNA-transcription rate. We find that the presence of extrinsic noise favours the presence of bimodal target distributions which can be observed for a wider range of parameters compared to the case with intrinsic noise only and for lower miRNA-target interaction strength. Our results suggest that combining threshold-inducing interactions with extrinsic noise provides a simple and robust mechanism for obtaining bimodal populations without requiring fine tuning. Furthermore, we characterise the protein distribution's dependence on protein half-life.

Integrated transcriptional and competitive endogenous RNA networks are cross-regulated in permissive molecular environments.

Competitive endogenous (ce)RNAs cross-regulate each other through sequestration of shared microRNAs and form complex regulatory networks based on their microRNA signature. However, the molecular requirements for ceRNA cross-regulation and the extent of ceRNA networks remain unknown. Here, we present a mathematical mass-action model to determine the optimal conditions for ceRNA activity in silico. This model was validated using phosphatase and tensin homolog (PTEN) and its ceRNA VAMP (vesicle-associated membrane protein)-associated protein A (VAPA) as paradigmatic examples. A computational assessment of the complexity of ceRNA networks revealed that transcription factor and ceRNA networks are intimately intertwined. Notably, we found that ceRNA networks are responsive to transcription factor up-regulation or their aberrant expression in cancer. Thus, given optimal molecular conditions, alterations of one ceRNA can have striking effects on integrated ceRNA and transcriptional networks.

A combination of transcriptional and microRNA regulation improves the stability of the relative concentrations of target genes.

It is well known that, under suitable conditions, microRNAs are able to fine tune the relative concentration of their targets to any desired value. We show that this function is particularly effective when one of the targets is a Transcription Factor (TF) which regulates the other targets. This combination defines a new class of feed-forward loops (FFLs) in which the microRNA plays the role of master regulator. Using both deterministic and stochastic equations, we show that these FFLs are indeed able not only to fine-tune the TF/target ratio to any desired value as a function of the miRNA concentration but also, thanks to the peculiar topology of the circuit, to ensure the stability of this ratio against stochastic fluctuations. These two effects are due to the interplay between the direct transcriptional regulation and the indirect TF/Target interaction due to competition of TF and target for miRNA binding (the so called "sponge effect"). We then perform a genome wide search of these FFLs in the human regulatory network and show that they are characterized by a very peculiar enrichment pattern. In particular, they are strongly enriched in all the situations in which the TF and its target have to be precisely kept at the same concentration notwithstanding the environmental noise. As an example we discuss the FFL involving E2F1 as Transcription Factor, RB1 as target and miR-17 family as master regulator. These FFLs ensure a tight control of the E2F/RB ratio which in turns ensures the stability of the transition from the G0/G1 to the S phase in quiescent cells.

Sodium acetate increases the productivity of HEK293 cells expressing the ECD-Her1 protein in batch cultures: experimental results and metabolic flux analysis.

Human Embryonic Kidney cells (HEK293) are a popular host for recombinant protein expression and production in the biotechnological industry. This has driven within both, the scientific and the engineering communities, the search for strategies to increase their protein productivity. The present work is inserted into this search exploring the impact of adding sodium acetate (NaAc) into a batch culture of HEK293 cells. We monitored, as a function of time, the cell density, many external metabolites, and the supernatant concentration of the heterologous extra-cellular domain ECD-Her1 protein, a protein used to produce a candidate prostate cancer vaccine. We observed that by adding different concentrations of NaAc (0, 4, 6 and 8 mM), the production of ECD-Her1 protein increases consistently with increasing concentration, whereas the carrying capacity of the medium decreases. To understand these results we exploited a combination of experimental and computational techniques. Metabolic Flux Analysis (MFA) was used to infer intracellular metabolic fluxes from the concentration of external metabolites. Moreover, we measured independently the extracellular acidification rate and oxygen consumption rate of the cells. Both approaches support the idea that the addition of NaAc to the culture has a significant impact on the metabolism of the HEK293 cells and that, if properly tuned, enhances the productivity of the heterologous ECD-Her1 protein.

Mammalian cell growth characterisation by a non-invasive plate reader assay.

Automated and non-invasive mammalian cell analysis is currently lagging behind due to a lack of methods suitable for a variety of cell lines and applications. Here, we report the development of a high throughput non-invasive method for tracking mammalian cell growth and performance based on plate reader measurements. We show the method to be suitable for both suspension and adhesion cell lines, and we demonstrate it can be adopted when cells are grown under different environmental conditions. We establish that the method is suitable to inform on effective drug treatments to be used depending on the cell line considered, and that it can support characterisation of engineered mammalian cells over time. This work provides the scientific community with an innovative approach to mammalian cell screening, also contributing to the current efforts towards high throughput and automated mammalian cell engineering.

The role of incoherent microRNA-mediated feedforward loops in noise buffering.

MicroRNAs are endogenous non-coding RNAs which negatively regulate the expression of protein-coding genes in plants and animals. They are known to play an important role in several biological processes and, together with transcription factors, form a complex and highly interconnected regulatory network. Looking at the structure of this network, it is possible to recognize a few overrepresented motifs which are expected to perform important elementary regulatory functions. Among them, a special role is played by the microRNA-mediated feedforward loop in which a master transcription factor regulates a microRNA and, together with it, a set of target genes. In this paper we show analytically and through simulations that the incoherent version of this motif can couple the fine-tuning of a target protein level with an efficient noise control, thus conferring precision and stability to the overall gene expression program, especially in the presence of fluctuations in upstream regulators. Among the other results, a nontrivial prediction of our model is that the optimal attenuation of fluctuations coincides with a modest repression of the target expression. This feature is coherent with the expected fine-tuning function and in agreement with experimental observations of the actual impact of a wide class of microRNAs on the protein output of their targets. Finally, we describe the impact on noise-buffering efficiency of the cross-talk between microRNA targets that can naturally arise if the microRNA-mediated circuit is not considered as isolated, but embedded in a larger network of regulations.

microRNA-Mediated Encoding and Decoding of Time-Dependent Signals in Tumorigenesis.

microRNAs, pivotal post-transcriptional regulators of gene expression, in the past decades have caught the attention of researchers for their involvement in different biological processes, ranging from cell development to cancer. Although lots of effort has been devoted to elucidate the topological features and the equilibrium properties of microRNA-mediated motifs, little is known about how the information encoded in frequency, amplitude, duration, and other features of their regulatory signals can affect the resulting gene expression patterns. Here, we review the current knowledge about microRNA-mediated gene regulatory networks characterized by time-dependent input signals, such as pulses, transient inputs, and oscillations. First, we identify the general characteristic of the main motifs underlying temporal patterns. Then, we analyze their impact on two commonly studied oncogenic networks, showing how their dysfunction can lead to tumorigenesis.

RAB11-Mediated Trafficking and Human Cancers: An Updated Review.

Many disorders block and subvert basic cellular processes in order to boost their progression. One protein family that is prone to be altered in human cancers is the small GTPase RAB11 family, the master regulator of vesicular trafficking. RAB11 isoforms function as membrane organizers connecting the transport of cargoes towards the plasma membrane with the assembly of autophagic precursors and the generation of cellular protrusions. These processes dramatically impact normal cell physiology and their alteration significantly affects the survival, progression and metastatization as well as the accumulation of toxic materials of cancer cells. In this review, we discuss biological mechanisms ensuring cargo recognition and sorting through a RAB11-dependent pathway, a prerequisite to understand the effect of RAB11 alterations in human cancers.

Initial cell density encodes proliferative potential in cancer cell populations.

Individual cells exhibit specific proliferative responses to changes in microenvironmental conditions. Whether such potential is constrained by the cell density throughout the growth process is however unclear. Here, we identify a theoretical framework that captures how the information encoded in the initial density of cancer cell populations impacts their growth profile. By following the growth of hundreds of populations of cancer cells, we found that the time they need to adapt to the environment decreases as the initial cell density increases. Moreover, the population growth rate shows a maximum at intermediate initial densities. With the support of a mathematical model, we show that the observed interdependence of adaptation time and growth rate is significantly at odds both with standard logistic growth models and with the Monod-like function that governs the dependence of the growth rate on nutrient levels. Our results (i) uncover and quantify a previously unnoticed heterogeneity in the growth dynamics of cancer cell populations; (ii) unveil how population growth may be affected by single-cell adaptation times; (iii) contribute to our understanding of the clinically-observed dependence of the primary and metastatic tumor take rates on the initial density of implanted cancer cells.

microRNA-mediated noise processing in cells: A fight or a game?

In the past decades, microRNAs (miRNA) have much attracted the attention of researchers at the interface between life and theoretical sciences for their involvement in post-transcriptional regulation and related diseases. Thanks to the always more sophisticated experimental techniques, the role of miRNAs as "noise processing units" has been further elucidated and two main ways of miRNA noise-control have emerged by combinations of theoretical and experimental studies. While on one side miRNAs were thought to buffer gene expression noise, it has recently been suggested that miRNAs could also increase the cell-to-cell variability of their targets. In this Mini Review, we focus on the role of miRNAs in molecular noise processing and on the advantages as well as current limitations of theoretical modelling.

From Endogenous to Synthetic microRNA-Mediated Regulatory Circuits: An Overview.

MicroRNAs are short non-coding RNAs that are evolutionarily conserved and are pivotal post-transcriptional mediators of gene regulation. Together with transcription factors and epigenetic regulators, they form a highly interconnected network whose building blocks can be classified depending on the number of molecular species involved and the type of interactions amongst them. Depending on their topology, these molecular circuits may carry out specific functions that years of studies have related to the processing of gene expression noise. In this review, we first present the different over-represented network motifs involving microRNAs and their specific role in implementing relevant biological functions, reviewing both theoretical and experimental studies. We then illustrate the recent advances in synthetic biology, such as the construction of artificially synthesised circuits, which provide a controlled tool to test experimentally the possible microRNA regulatory tasks and constitute a starting point for clinical applications.

Kinetic Modelling of Competition and Depletion of Shared miRNAs by Competing Endogenous RNAs.

Non-coding RNAs play a key role in the post-transcriptional regulation of mRNA translation and turnover in eukaryotes. miRNAs, in particular, interact with their target RNAs through protein-mediated, sequence-specific binding, giving rise to extended and highly heterogeneous miRNA-RNA interaction networks. Within such networks, competition to bind miRNAs can generate an effective positive coupling between their targets. Competing endogenous RNAs (ceRNAs) can in turn regulate each other through miRNA-mediated crosstalk. Albeit potentially weak, ceRNA interactions can occur both dynamically, affecting, e.g., the regulatory clock, and at stationarity, in which case ceRNA networks as a whole can be implicated in the composition of the cell's proteome. Many features of ceRNA interactions, including the conditions under which they become significant, can be unraveled by mathematical and in silico models. We review the understanding of the ceRNA effect obtained within such frameworks, focusing on the methods employed to quantify it, its role in the processing of gene expression noise, and how network topology can determine its reach.

Stochastic sequestration dynamics: a minimal model with extrinsic noise for bimodal distributions and competitors correlation.

Many biological processes are known to be based on molecular sequestration. This kind of dynamics involves two types of molecular species, namely targets and sequestrants, that bind to form a complex. In the simple framework of mass-action law, key features of these systems appear to be threshold-like profiles of the amounts of free molecules as a function of the parameters determining their possible maximum abundance. However, biochemical processes are probabilistic and take place in stochastically fluctuating environments. How these different sources of noise affect the final outcome of the network is not completely characterised yet. In this paper we specifically investigate the effects induced by a source of extrinsic noise onto a minimal stochastic model of molecular sequestration. We analytically show how bimodal distributions of the targets can appear and characterise them as a result of noise filtering mediated by the threshold response. We then address the correlations between target species induced by the sequestrant and discuss how extrinsic noise can turn the negative correlation caused by competition into a positive one. Finally, we consider the more complex scenario of competitive inhibition for enzymatic kinetics and discuss the relevance of our findings with respect to applications.

Somatic mutagenesis in satellite cells associates with human skeletal muscle aging.

Human aging is associated with a decline in skeletal muscle (SkM) function and a reduction in the number and activity of satellite cells (SCs), the resident stem cells. To study the connection between SC aging and muscle impairment, we analyze the whole genome of single SC clones of the leg muscle vastus lateralis from healthy individuals of different ages (21-78 years). We find an accumulation rate of 13 somatic mutations per genome per year, consistent with proliferation of SCs in the healthy adult muscle. SkM-expressed genes are protected from mutations, but aging results in an increase in mutations in exons and promoters, targeting genes involved in SC activity and muscle function. In agreement with SC mutations affecting the whole tissue, we detect a missense mutation in a SC propagating to the muscle. Our results suggest somatic mutagenesis in SCs as a driving force in the age-related decline of SkM function.

RNAs competing for microRNAs mutually influence their fluctuations in a highly non-linear microRNA-dependent manner in single cells.

BACKGROUND: Distinct RNA species may compete for binding to microRNAs (miRNAs). This competition creates an indirect interaction between miRNA targets, which behave as miRNA sponges and eventually influence each other's expression levels. Theoretical predictions suggest that not only the mean expression levels of targets but also the fluctuations around the means are coupled through miRNAs. This may result in striking effects on a broad range of cellular processes, such as cell differentiation and proliferation. Although several studies have reported the functional relevance of this mechanism of interaction, detailed experiments are lacking that study this phenomenon in controlled conditions by mimicking a physiological range. RESULTS: We used an experimental design based on two bidirectional plasmids and flow cytometry measurements of cotransfected mammalian cells. We validated a stochastic gene interaction model that describes how mRNAs can influence each other's fluctuations in a miRNA-dependent manner in single cells. We show that miRNA-target correlations eventually lead to either bimodal cell population distributions with high and low target expression states, or correlated fluctuations across targets when the pool of unbound targets and miRNAs are in near-equimolar concentration. We found that there is an optimal range of conditions for the onset of cross-regulation, which is compatible with 10-1000 copies of targets per cell. CONCLUSIONS: Our results are summarized in a phase diagram for miRNA-mediated cross-regulation that links experimentally measured quantities and effective model parameters. This phase diagram can be applied to in vivo studies of RNAs that are in competition for miRNA binding.

Modelling Competing Endogenous RNA Networks.

MicroRNAs (miRNAs) are small RNA molecules, about 22 nucleotide long, which post-transcriptionally regulate their target messenger RNAs (mRNAs). They accomplish key roles in gene regulatory networks, ranging from signaling pathways to tissue morphogenesis, and their aberrant behavior is often associated with the development of various diseases. Recently it has been experimentally shown that the way miRNAs interact with their targets can be described in terms of a titration mechanism. From a theoretical point of view titration mechanisms are characterized by threshold effect at near-equimolarity of the different chemical species, hypersensitivity of the system around the threshold, and cross-talk among targets. The latter characteristic has been lately identified as competing endogenous RNA (ceRNA) effect to mark those indirect interactions among targets of a common pool of miRNAs they are in competition for. Here we propose a stochastic model to analyze the equilibrium and out-of-equilibrium properties of a network of [Formula: see text] miRNAs interacting with [Formula: see text] mRNA targets. In particular we are able to describe in detail the peculiar equilibrium and non-equilibrium phenomena that the system displays in proximity to the threshold: (i) maximal cross-talk and correlation between targets, (ii) robustness of ceRNA effect with respect to the model's parameters and in particular to the catalyticity of the miRNA-mRNA interaction, and (iii) anomalous response-time to external perturbations.