# Probing the Limits to MicroRNA-Mediated Control of Gene Expression

Araks Martirosyan, Matteo Figliuzzi, Enzo Marinari, Andrea de Martino

► To cite this version:

Araks Martirosyan, Matteo Figliuzzi, Enzo Marinari, Andrea de Martino. Probing the Limits to MicroRNA-Mediated Control of Gene Expression. PLoS Computational Biology, 2016, 12 (1), pp.e1004715. 10.1371/journal.pcbi.1004715. hal-01527410

HAL Id: hal-01527410

https://hal.sorbonne-universite.fr/hal-01527410v1

Submitted on 24 May 2017

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons CC BY 4.0 - Attribution - International License
<figure><img src="image_2.png" /></figure>

CrossMark

click for updates

<figure><img src="image_3.png" /></figure>

OPEN ACCESS

**Citation:** Martirosyan A, Figliuzzi M, Marinari E, De Martino A (2016) Probing the Limits to MicroRNA-Mediated Control of Gene Expression. PLoS Comput Biol 12(1): e1004715. doi:10.1371/journal.pcbi.1004715

**Editor:** Edwin Wang, National Research Council of Canada, CANADA

**Received:** July 15, 2015

**Accepted:** December 20, 2015

**Published:** January 26, 2016

**Copyright:** © 2016 Martirosyan et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** All relevant data are within the paper and its Supporting Information files.

**Funding:** This work was supported by Marie Curie Action Initial Training Networks NETworks Across Disciplines, grant agreement n 290038 of the European Union Seventh Framework Programme 2007-2013. www.netadis.eu. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

**Competing Interests:** The authors have declared that no competing interests exist.

RESEARCH ARTICLE

# Probing the Limits to MicroRNA-Mediated Control of Gene Expression

Araks Martirosyan¹,²*, Matteo Figliuzzi³,⁴,⁵, Enzo Marinari¹,⁶*, Andrea De Martino¹,²,⁷*

1 Dipartimento di Fisica, Sapienza Università di Roma, Rome, Italy, 2 Soft and Living Matter Lab, Institute of Nanotechnology (CNR-NANOTEC), Consiglio Nazionale delle Ricerche, Rome, Italy, 3 Sorbonne Universités, UPMC, Institut de Calcul et de la Simulation, Paris, France, 4 Sorbonne Universités, UPMC, UMR 7238, Computational and Quantitative Biology, Paris, France, 5 CNRS, UMR 7238, Computational and Quantitative Biology, Paris, France, 6 INFN, Sezione di Roma 1, Rome, Italy, 7 Center for Life Nano Science @Sapienza, Istituto Italiano di Tecnologia, Rome, Italy

* These authors contributed equally to this work.

araks.martirosyan@cnr.it

## Abstract

According to the 'ceRNA hypothesis', microRNAs (miRNAs) may act as mediators of an effective positive interaction between long coding or non-coding RNA molecules, carrying significant potential implications for a variety of biological processes. Here, inspired by recent work providing a quantitative description of small regulatory elements as information-conveying channels, we characterize the effectiveness of miRNA-mediated regulation in terms of the optimal information flow achievable between modulator (transcription factors) and target nodes (long RNAs). Our findings show that, while a sufficiently large degree of target derepression is needed to activate miRNA-mediated transmission, (a) in case of differential mechanisms of complex processing and/or transcriptional capabilities, regulation by a post-transcriptional miRNA-channel can outperform that achieved through direct transcriptional control; moreover, (b) in the presence of large populations of weakly interacting miRNA molecules the extra noise coming from titration disappears, allowing the miRNA-channel to process information as effectively as the direct channel. These observations establish the limits of miRNA-mediated post-transcriptional cross-talk and suggest that, besides providing a degree of noise buffering, this type of control may be effectively employed in cells both as a failsafe mechanism and as a preferential fine tuner of gene expression, pointing to the specific situations in which each of these functionalities is maximized.

## Author Summary

The discovery of RNA interference has revolutionized the decades' old view of RNAs as mere intermediaries between DNA and proteins in the gene expression workflow. MicroRNAs (or miRNAs), in particular, have been shown to be able to both stabilize the protein output by buffering transcriptional noise and to create an effective positive interaction between the levels of their target RNAs through a simple competition mechanism known as 'ceRNA effect'. With miRNAs commonly targeting multiple species of RNAs, the
potential implication is that RNAs could regulate each other through extended miRNA-mediated interaction networks. Such cross-talk is certainly active in many specific cases (like cell differentiation), but it's unclear whether the degree of regulation of gene expression achievable through post-transcriptional miRNA-mediated coupling can effectively overcome the one obtained through other mechanisms, e.g. by direct transcriptional control via DNA-binding factors. This work quantifies the maximal post-transcriptional regulatory power achievable by miRNA-mediated cross-talk, characterizing the circumstances in which indirect control outperforms direct one. The emerging scenario suggests that, in addition to its widely recognized noise-buffering role, miRNA-mediated control may indeed act as a master regulator of gene expression.

## Introduction

The problem of tuning protein expression levels is central for eukaryotic cell functionality. A variety of molecular mechanisms are implemented to guarantee, on one hand, that protein copy numbers stay within a range that is optimal in the given conditions and, on the other, that shifts in expression levels can be achieved efficiently whenever necessary [1–3] (whereby ‘efficiency’ here encompasses both a dynamical characterization, in terms of the times required to shift, and a static one, in terms of moving as precisely as possible from one functional range to another). Quantifying and comparing their effectiveness in different conditions is an important step to both deepen our fundamental understanding of regulatory circuits and to get case-by-case functional insight about why a specific biochemical network has been selected over the others.

As the major direct regulators of gene expression, transcription factors (TFs) are most immediately identified as the key potential modulators of protein levels [4]. In a somewhat simplified picture, one may imagine that a change in amount of a TF can induce a change in the expression level of the corresponding gene, and that the ability to regulate the latter (the output node) via the former (the input node) can be assessed by how strongly the two levels correlate. The effectiveness of a regulatory element is however limited by the stochasticity of intracellular processes, from the TF-DNA binding dynamics to translation [5]. A convenient framework to analyze how noise constrains regulation is provided by information theory [6, 7]. In particular, the simplest situation in which a single TF modulates the expression of a single protein can be characterized analytically under the assumption that the noise affecting the input-output channel is sufficiently small. The mutual information between modulator and target –a convenient quantity through which regulatory effectiveness can be characterized– depends on the distribution of modulator levels and can be maximized over it. Remarkably, in at least one case this maximum has been found to be almost saturated by the actual information flow measured in a living system (for more details see [8, 9]). In other terms, for sufficiently small noise levels in the channel that links TFs to their functional products, one may quantify the optimal regulatory performance achievable in terms of the maximum number of bits of mutual information that can be exchanged between modulator and target.

Several control mechanisms however act at the post-transcriptional level [10–12]. Among these, regulation by small regulatory RNAs like eukaryotic microRNAs (miRNAs) has attracted considerable attention over the past few years [13–15]. In short, miRNAs are small non-coding RNA molecules encoded by nuclear DNA, that can inhibit translation or catalyze degradation of mRNAs when bound to them via protein-mediated base-pairing. miRNAs appear to be crucial in an increasing number of situations ranging from development to disease [16–18]. Their function however can differ significantly from case to case. For instance, they have been well
characterized as noise buffering agents in protein expression [19, 20] or as key signaling molecules in stress response [21]. Recently, though, investigations of the miRNA-mediated post-transcriptional regulatory (PTR) network have hinted at a possibly more subtle and complex role. It is indeed now clear [22–27] that the miRNA-RNA network describing the potential couplings stretches across a major fraction of the transcriptome, with a large heterogeneity both in the number of miRNA targets and in the number of miRNA regulators for a given mRNA. The competition effects that may emerge in such conditions suggest that miRNAs may act as channels through which perturbations in the levels of one RNA could be transmitted to other RNA species sharing the same miRNA regulator(s). Such a scenario has been termed the ‘ceRNA effect’, whereby ceRNA stands for ‘competing endogenous RNA’ [28]. In view of its considerable regulatory and therapeutic implications, the ceRNA effect has been extensively analyzed both theoretically and experimentally [29–44].

The apparent ubiquity of potentially cross-talking ceRNAs however raises a number of fundamental questions about the effectiveness of “regulation via competition” *per se*. Although hundreds of targets are predicted for a single miRNA, observations show that only few of them are sensitive to changes in miRNA expression levels. Most targets are likely to provide a global buffering mechanism through which miRNA levels are overall stabilized [28, 29]. Effective competition between miRNA targets requires that the ratio of miRNA molecules to the number of target sites lies in a specific range, so that the relative abundance of miRNA and RNA species must be tightly regulated for the ceRNA mechanism to operate [29–34]. On the other hand, the magnitude of the ceRNA effect is tunable by the miRNA binding and mRNA loss rates [33, 34, 42]. The performance of a regulatory element, however, does not only depend on kinetic parameters, but also on the range of variability (and possibly on the distribution itself) of modulator levels (e.g. TFs) [8, 9]. The maximal regulatory effectiveness of a given genetic circuit – quantifying how precisely the output level can be determined by the input level– can therefore generically be obtained by solving an optimization problem over the distribution of inputs. This type of approach provides an upper bound to the effectiveness of a regulatory mechanism as well as indications concerning which parameters, noise sources and/or interactions most hamper its performance.

It would be especially important to understand in which conditions the degree of control of the output variable (i.e. the ceRNA/protein level) that can be accomplished through post-transcriptional miRNA-mediated cross-talk may exceed that obtainable by different regulatory mechanisms. In this work we characterize the maximal regulatory power achievable by miRNA-mediated control and compare it with that of a direct, TF-based transcriptional unit [45]. In principle, since fluctuations can be reduced by increasing the number of molecules, an (almost) arbitrary amount of information can be transmitted through a biochemical network. However, cells have to face the burden of macromolecular synthesis [46–48]. Optimality is therefore the result of a trade-off between the benefits of reduced fluctuations and the drawbacks of the associated metabolic costs. For this reason, we start by fixing a maximal rate of transcription (or, alternatively, the maximal number of output molecules) so as to have a simple but reasonable framework to characterize and compare the capacities of the different regulatory channels. Next, we quantify how an input signal is processed by the transcriptional (TF-based) and post-transcriptional (miRNA-mediated) regulatory elements by characterizing the response in the output ceRNA’s expression levels. In such a setting, information flow is hampered by intrinsic noise if the target gene is weakly derepressed by the activation of its competitor. Otherwise, target derepression appears to have a strong impact on a regulatory element’s capacity. Upon varying the magnitude of derepression by tuning the kinetic parameters, we then show that in certain regimes miRNA-mediated regulation can indeed outperform direct control of gene expression. Finally, we argue that the presence of miRNA molecules in large
copy numbers notably reduces the level of intrinsic noise on weakly targeted transcripts. In this case, the mutual regulation of ceRNA molecules by miRNA-mediated channels may become a primary mechanism to finely tune gene expression.

Besides providing a quantitative characterization of the maximal regulatory power achievable through miRNA-based post-transcriptional control, these results provide important hints on the circumstances in which regulation by small RNAs may function as the main tuner of gene expression in cells.

# Results

## Mathematical model of ceRNA competition

We consider (see Fig 1) a system formed by two ceRNA species (whose levels are labeled $m_1$ and $m_2$, respectively) and one miRNA species (with level labeled $\mu$), whose transcriptions are activated by a single TF each (with levels labeled, respectively, $f_1, f_2$ and $f_\mu$). Both ceRNAs are in turn targeted by the miRNA. miRNA-ceRNA complexes (levels labeled $c_i$ with $i = 1,2$) assemble and disassemble at rates $k_i^+$, respectively, whereas complexes can be degraded both stoichio-metrically (i.e. without miRNA recycling) at rates $\sigma_i$ and catalytically (i.e. with miRNA recy-cling) at rates $\kappa_i$ [49]. In addition, ceRNA and miRNA molecules degrade (resp. synthesize) at rates $d_i$ (resp. $b_i$) and $\delta$ (resp. $\beta$), respectively. Steps leading to the formation of the RNA-induced silencing complex (RISC), allowing for the miRNA-ceRNA binding, are neglected for simplicity. TF levels are treated as externally controlled parameters.

Denoting the TF-DNA binding/unbinding rates by $k_{in}$ and $k_{out}$, respectively (for simplicity, these parameters are taken to be the same for all involved TFs), the TF binding sites' fractional

<figure><img src="image_4.png" /><figcaption>**Fig 1.** Schematic representation of the basic ceRNA network model. The model comprises three transcription factors (TF₁, TF₂ and TFµ) controlling the synthesis of two ceRNA species (ceRNA₁ and ceRNA₂) and one miRNA species, respectively. The fractional occupancies of TF binding sites are denoted by $n_1$, $n_2$ and $n_\mu$, respectively. Both ceRNAs are targets for the miRNA, with whom they form complexes denoted respectively as $C_1$ and $C_2$. The amount of molecules for each species is denoted by the variable next to the corresponding node. Reaction rates are reported next to the corresponding arrow.</figcaption></figure>

doi:10.1371/journal.pcbi.1004715.g001
occupancies $n_\ell$ ($0 \le n_\ell \le 1$, $\ell \in \{1, 2, \mu\}$) obey the dynamics

$$ \frac{dn_\ell}{dt} = k_{\text{in}}(1 - n_\ell)f_\ell^h - k_{\text{out}}n_\ell, \qquad (1) $$

according to which transcriptional activation requires the cooperative binding of h TF molecules for each RNA species involved. In general, the occupancies by different TFs equilibrate on different timescales [50]. However it is often assumed that the transcriptional on/off dynamics is much faster than transcription itself [51, 52]. As a consequence, each $n_\ell$ can be fixed at its ‘equilibrium’ value.

$$ \bar{n}_{\ell} = \frac{k_{\text{in}}f_{\ell}^{h}}{k_{\text{in}}f_{\ell}^{h} + k_{\text{out}}}. \qquad (2) $$

Every process in the above scheme contributes to the overall level of noise. We represent the mass-action kinetics of the system through the set of coupled Langevin processes ($i=1,2$) [41, 53]

$$ \begin{aligned} \frac{dm_i}{dt} &= -d_i m_i + b_i \bar{n}_i - k_i^+ \mu m_i + k_i^- c_i + \xi_i - \xi_i^+ + \xi_i^-, \\ \frac{dc_i}{dt} &= k_i^+ \mu m_i - (k_i^- + \kappa_i + \sigma_i) c_i + \xi_i^\sigma + \xi_i^+ - \xi_i^- - \xi_i^\kappa, \\ \frac{d\mu}{dt} &= -\delta\mu + \beta \bar{n}_\mu - \sum_i k_i^+ \mu m_i + \sum_i (k_i^- + \kappa_i) c_i + \xi_\mu - \sum_i \xi_i^+ + \sum_i \xi_i^- + \sum_i \xi_i^\kappa, \end{aligned} \qquad (3) $$

where the mutually independent random (Poisson) ‘forces’ $\xi_b$, $\xi_i^\pm$, $\xi_\mu$, $\xi_i^\kappa$ and $\xi_i^\sigma$ denote, respectively, the intrinsic noise in ceRNA levels (due to random synthesis and degradation events), in the association/dissociation processes of complexes, in the miRNA level, in the catalytic complex decay and in the stoichiometric complex decay. Each of the above noise terms has zero mean, while correlations are given by

$$ \langle \xi_i(t) \xi_i(t') \rangle = (d_i \bar{m}_i + b_i \bar{n}_i) \delta(t - t'), \qquad (4) $$

$$ \langle \xi_i^+(t) \xi_i^+(t') \rangle = k_i^+ \bar{m}_i \bar{\mu} \delta(t - t'), \qquad (5) $$

$$ \langle \xi_i^-(t) \xi_i^-(t') \rangle = k_i^- \bar{c}_i \delta(t - t'), \qquad (6) $$

$$ \langle \xi_\mu(t) \xi_\mu(t') \rangle = (\delta\bar{\mu} + \beta\bar{n}_\mu) \delta(t - t'), \qquad (7) $$

$$ \langle \xi_i^\kappa(t) \xi_i^\kappa(t') \rangle = \kappa_i \bar{c}_i \delta(t - t'), \qquad (8) $$

$$ \langle \xi_i^\sigma(t) \xi_i^\sigma(t') \rangle = \sigma_i \bar{c}_i \delta(t - t'), \qquad (9) $$

where we introduced the steady state molecule numbers

$$ \bar{m}_i = \frac{b_i \bar{n}_i + k_i^- \bar{c}_i}{d_i + k_i^+ \bar{\mu}}, \qquad (10) $$

$$ \bar{\mu} = \frac{\beta \bar{n}_{\mu} + \sum_{i} (k_{i}^{-} + \kappa_{i}) \bar{c}_{i}}{\delta + \sum_{i} k_{i}^{+} \bar{m}_{i}}, \qquad (11) $$

$$ \bar{c}_i = \frac{k_i^+ \bar{\mu} \bar{m}_i}{\sigma_i + k_i^- + \kappa_i}. \qquad (12) $$
<figure><img src="image_5.png" /><figcaption>**Fig 2. Cross-talk scenario in the ceRNA network.** Steady state values of (A) the ceRNA and miRNA levels, and (B) the corresponding Fano factors. Markers denote GA computations, curves correspond to approximate analytical solutions obtained by the linear noise approximation. Values of the kinetic parameters are reported in Table 1.</figcaption></figure>

doi:10.1371/journal.pcbi.1004715.g002

We shall be interested in the fluctuations of molecular levels around the steady state that arise due to intrinsic noise sources. The stochastic dynamics of the system can be simulated via the Gillespie algorithm (GA, see Materials and Methods). Fig 2A shows typical GA results for $\bar{m}_1$, $\bar{m}_2$ and $\bar{\mu}$, with the corresponding Fano factors (FFs) shown in Fig 2B, as functions of $f_1$. Numerical results are matched against analytical estimates obtained by the linear noise approximation (see Materials and Methods).

In particular, in Fig 2 one sees that an increase of $\bar{m}_1$ is accompanied by a concomitant increase of $\bar{m}_2$ and by a decrease of the average number of free miRNA molecules. This is an instance of miRNA-mediated ceRNA cross-talk. Indeed, upon up-regulating $m_1$ by injecting the corresponding TF (i.e. by increasing $f_1$), the level of free miRNAs will decrease as more and more molecules will be actively repressing ceRNA$_1$, causing in turn an up-regulation of ceRNA$_2$. $f_1$ will thus positively correlate with $m_2$. One easily sees that steady-state ceRNA levels depend on $\bar{\mu}$ through a sigmoidal function, namely [42]

$$ \bar{m}_i = \frac{b_i \bar{n}_i}{d_i} F_i[\bar{\mu}], \quad F_i[\bar{\mu}] = \frac{\mu_{0,i}}{\mu_{0,i} + \bar{\mu}}. \qquad (13) $$

The constant $\mu_{0,i} = \frac{d_i}{k_i^+} \left( 1 + \frac{k_i^-}{\sigma_{i+\kappa_i}} \right)$ acts as a 'soft' threshold for the miRNA level, allowing to distinguish three situations:

* if $\bar{\mu} \ll \mu_{0,i}$, ceRNA$_i$ is free or unrepressed: spontaneous ceRNA degradation dominates over miRNA-mediated decay channels, so that, effectively, the ceRNA level is weakly sensitive to small changes in $\bar{\mu}$;

* if $\bar{\mu} \gg \mu_{0,i}$, ceRNA$_i$ is bound or repressed: miRNA-mediated ceRNA decay dominates over spontaneous ceRNA degradation but, again, the ceRNA level is weakly sensitive to small changes in $\bar{\mu}$ as most ceRNAs are bound in complexes with the miRNA;

* if $\bar{\mu} \simeq \mu_{0,i}$, ceRNA$_i$ is susceptible to $\mu$: spontaneous and miRNA-mediated decay channels have comparable weights and the ceRNA level is very sensitive to small changes in $\bar{\mu}$.

The behaviour of the FFs (see Fig 2B) emphasizes how noise patterns change in the different regimes. In the displayed example, ceRNA$_1$ and ceRNA$_2$ become susceptible for $\ln(f_1) \simeq 2.5$ and $\ln(f_1) \simeq 2.1$, respectively. Indeed, one observes that the corresponding FFs peak close to these
<figure><img src="image_6.png" /><figcaption>Fig 3. Schematic representation of the channels under consideration. (A) miRNA-channel, (B) TF-channel.</figcaption></figure>

doi:10.1371/journal.pcbi.1004715.g003

values, in accordance with the observation that stochastic fluctuations are enhanced when the rates of substrate supply are adequately balanced in a stoichiometrically coupled system [34, 54, 55]. The FF for ceRNA₂ appears to approach one for large values of $f_1$, as expected for the pure Poisson birth/death process that characterizes the free regime [56]. On the other hand, for very small but nonzero mean fractional occupancy $\bar{n}_1$ (corresponding to small values of $f_1$), $m_1$ will with high probability only take on the values 0 or 1, as a transcribed molecule will quickly undergo degradation or sequestration in a complex. In such a situation, the mean and variance of $m_1$ will be calculated by summing up zeros and ones over time, leading to a FF equal to one.

## Measuring the performance of regulatory elements

**The amplitude of ceRNA derepression.** We shall focus on quantifying the regulatory capacity of the post-transcriptional and transcriptional channels that link, respectively, TF₁ to ceRNA₂ and TF₂ to ceRNA₂, as depicted in Fig 3.

In what follows, we shall take the number of transcription factors as the key input variables ($f_1$ for the miRNA-mediated channel and $f_2$ for the direct channel), and the number of target RNAs $m_2$ as the key output variable. It is reasonable to think that an effective regulator should have a sensible impact on the target molecule. As a first measure of regulatory performance we will therefore consider the *amplitude of variation* (AOV) of $m_i$ via $f_j$, defined as

$$ \Delta_{ij} = \bar{m}_i(f_j^{\max}) - \bar{m}_i(f_j^{\min}), \qquad (14) $$

corresponding to the difference between the maximum and the minimum steady state level of ceRNAᵢ obtained by varying $f_j$ from $f_j^{\min}$ to $f_j^{\max}$. Clearly, larger values of $\Delta_{ij}$ imply stronger degrees of derepression of ceRNAᵢ induced by the corresponding TF. For sake of simplicity, we shall adopt the following notation to distinguish the AOV induced by transcriptional and miRNA-mediated control channels:

$$ \Delta_{22} \equiv \Delta_{\text{TF}}, \qquad (15) $$

$$ \Delta_{21} \equiv \Delta_{\text{miRNA}}, \qquad (16) $$

Note that the AOV only involves mean levels. As such, it is unable to account for the intrinsic degree of stochasticity of the system.

**Maximal mutual information.** In order to evaluate the impact of noise, as a further quantifier of regulatory power we shall employ the mutual information (MI), a standard
information-theoretic measure of channel performance given by $(i, j = 1, 2)$ [52]

$$I(m_i, f_j) = \int p(m_i, f_j) \log_2 \frac{p(m_i, f_j)}{p(m_i)p(f_j)} dm_i df_j, \quad (17)$$

where $p(m_i, f_j)$ stands for the joint probability distribution of $m_i$ and $f_j$, while $p(m_i)$ and $p(f_j)$ denote its marginals. The MI is a non-negative quantity that vanishes when $m_i$ and $f_j$ are independent. Conversely, when $I(m_i, f_j) = \mathcal{I}$ one can reliably distinguish (roughly) $2^\mathcal{I}$ different values of the output variable upon tuning the input (see [9, 45, 52] for further details about the usefulness of this scheme for the analysis of transcriptional regulatory or biochemical modules). The maximal achievable information flow for a fixed channel (i.e. for given kinetic parameters) is called the *channel capacity* and can be computed by maximising the MI over the input variable (in this case, over the distribution of the TF level). In the simplest case, one can assume that the channel is Gaussian, meaning that for every given input $f_j$ the output $m_2$ is described by the conditional probability density

$$p(m_2|f_j) = \frac{1}{\sqrt{2\pi\sigma_{m_2}^2}} \exp\left[-\frac{(m_2 - \bar{m}_2(f_j))^2}{2\sigma_{m_2}^2(f_j)}\right], \quad (18)$$

where both the average $\bar{m}_2(f_j)$ and variance $\sigma_{m_2}^2(f_j) = \langle (m_2 - \bar{m}_2(f_j))^2 \rangle$ depend on the input variable $f_j$. In the limit of small variance, the distribution $p_{\text{opt}}(f_j)$ of TF levels that maximises the MI between $m_2$ and $f_j$ reads [45]

$$p_{\text{opt}}(f_j) = \frac{1}{Z} \left( \frac{1}{2\pi e \sigma_{m_2}^2(f_j)} \left[ \frac{d\bar{m}_2(f_j)}{df_j} \right]^2 \right)^{1/2}, \quad (19)$$

where

$$Z = \int_{f_j^{\text{min}}}^{f_j^{\text{max}}} df_j \left( \frac{1}{2\pi e \sigma_{m_2}^2(f_j)} \left[ \frac{d\bar{m}_2(f_j)}{df_j} \right]^2 \right)^{1/2} \quad (20)$$

is a normalizing factor. For $p = p_{\text{opt}}$ one obtains the channel capacity

$$I_{\text{opt}}(m_2, f_j) = \log_2 Z. \quad (21)$$

S1 Fig provides a qualitative yet intuitive representation of how the shape of the input/output relation and the size of fluctuations affect information flow. Intuitively, the normalizing factor $Z$ ‘counts’ the number of distinguishable output levels comparing the local slope of the input/output curve with the noise strength: if the slope is too low or the fluctuations are too high no information can be transmitted through the channel. We shall use the following notation to distinguish the capacities of the direct and indirect channels:

$$I_{\text{opt}}(m_2, f_2) \equiv I_{\text{TF}}, \quad (22)$$

$$I_{\text{opt}}(m_2, f_1) \equiv I_{\text{miRNA}}. \quad (23)$$

We shall furthermore denote their difference by $\Delta I = I_{\text{TF}} - I_{\text{miRNA}}$.

We are henceforth going to adopt the following protocol. After generating the output data $(\bar{m}_2, \sigma_{m_2})$ by simulating the circuit’s dynamics via the GA, the optimal input distribution will be constructed according to Eq (19). We shall then let the circuit process an input signal sampled
<figure><img src="image_7.png" /><figcaption>**Fig 4.** The flowchart of the method. A C++ implementation is available at https://github.com/araksm/ceRNA.</figcaption></figure>

doi:10.1371/journal.pcbi.1004715.g004

from that distribution and calculate the MI for both the transcriptional and post-transcriptional channels (see the flowchart in Fig 4). Channel capacities will then be compared, with the goal of clarifying the conditions under which the former can be more effective than the latter.

Comparing indirect miRNA-mediated regulation with direct transcriptional control

Target prediction algorithms suggest that miRNA binding affinities vary significantly among their targets [57]. Such heterogeneities indeed are mapped to the miRNA binding kinetics and are shown to influence the susceptibility of the targets to the miRNA molecules [36]. Moreover, the level of complementarity between the regulator and the target seems to be decisive for the selection of a decay channel (catalytic or stoichiometric) for the miRNA-ceRNA complex [15, 27]. One may therefore expect that the effectiveness of miRNA-mediated post-transcriptional control depends strongly on the kinetic parameters characterizing the network.

Here in particular, we are going to investigate how the capacities of these regulatory ele-
ments are affected by changes in (i) miRNA-ceRNA binding kinetics, (ii) miRNA recycling
rates, and (iii) effective transcription rates of all RNA species involved.
<figure><img src="image_8.png" /><figcaption>**Fig 5.** Dependence of the transcriptional and post-transcriptional channel capacities on miRNA-ceRNA association rates. (A) *I*<sub>miRNA</sub>. (B) *I*<sub>TF</sub>. (C) *Δ*I/(*I*<sub>TF</sub> − *I*<sub>miRNA</sub>). Values of the kinetic parameters are reported in Table 1.</figcaption></figure>

doi:10.1371/journal.pcbi.1004715.g005

In order to contrast the performances of miRNA- and TF-channels we shall start by analyzing their respective responses to the same input signal. More precisely, the input variable *f*<sub>j</sub> will be varied from 0 to a value *f*<sub>j</sub><sup>max</sup> defined by the condition *n*<sub>j</sub>(*f*<sub>j</sub><sup>max</sup>) = 0.99 (i.e. from a situation in which the promoter is always free to one in which it is essentially always occupied).

**The capacity of the miRNA-channel is maximal in a specific range of miRNA-ceRNA binding rates.** Available estimates of RNA binding energetics indicate that the target sites of an individual small non coding RNA can span a wide range of affinity both in mammalian cells [57] and in bacteria [58]. To explore the functional consequences of such heterogeneities we analyse the behaviour of the regulatory channels changing the values of the complex association rates *k*<sub>1</sub><sup>+</sup> and *k*<sub>2</sub><sup>+</sup> while keeping fixed all others parameters. As displayed in Fig 5, the maximal amount of information transmitted through the post-transcriptional channel (Fig 5A) increases monotonically with *k*<sub>1</sub><sup>+</sup> (for fixed *k*<sub>2</sub><sup>+</sup>) and displays a maximum versus *k*<sub>2</sub><sup>+</sup> (for fixed *k*<sub>1</sub><sup>+</sup>). Indeed, the higher *k*<sub>1</sub><sup>+</sup>, the higher *c*<sub>1</sub>, and the more information about changes in *f*<sub>1</sub> can be transmitted through the miRNA-channel. On the other hand, ceRNA<sub>2</sub> is unrepressed by the miRNA for small values of *k*<sub>2</sub><sup>+</sup> while it is strongly repressed for high values of *k*<sub>2</sub><sup>+</sup>, for any value of the input. In both regimes, ceRNA<sub>2</sub> is only weakly sensitive to changes in the level of its competitor. In such conditions, no information can be transmitted through the miRNA-channel. The information flow is therefore maximal between these two extremes.

For the transcriptional channel, we expect to see maximum information transmission when the target, ceRNA₂, is unrepressed by the miRNA, in which case the channel is effectively decoupled from the rest of the network, while no information transmission is expected when the target is fully repressed by the miRNA. Indeed, in Fig 5B one sees that, for small values of *k*<sub>2</sub><sup>+</sup> (unrepressed target), *I*<sub>TF</sub> is the highest and it decreases to zero as *k*<sub>2</sub><sup>+</sup> increases together with the degree of target repression. Moreover, since a higher *k*<sub>1</sub><sup>+</sup> moves *m*<sub>2</sub> to larger values (i.e. towards the free regime), *I*<sub>TF</sub> increases with *k*<sub>1</sub><sup>+</sup> and decreases with *k*<sub>2</sub><sup>+</sup> in the susceptible regime. Notice that, when post-transcriptional capacity is highest, miRNA-mediated regulation is as effective as transcriptional control (see Fig 5C).

The intuition that the channel capacity is strongly coupled to the target's size of derepression (its AOV) is confirmed by the fact that the latter indeed displays a very similar behaviour (see S2 Fig).
<figure><img src="image_9.png" /><figcaption>**Fig 6.** Dependence of the transcriptional and post-transcriptional channel capacities on miRNA recycling rates. (A) *I*<sub>miRNA</sub>. (B) *I*<sub>TF</sub>. (C) *Δ*I* = *I*<sub>TF</sub> - *I*<sub>miRNA</sub>. Values of the kinetic parameters are reported in Table 1.</figcaption></figure>

doi:10.1371/journal.pcbi.1004715.g006

**miRNA-mediated regulation may represent the sole control mechanism in case of differential complex processing.** Direct experimental measurements of miRNA-recycling rates are challenging [59]. However it is well known that the degree of complementarity of a miRNA to its target lead to different repression pathways [60] and that differential complex processing strongly affects reciprocity of competition [33]. In order to better quantify the degree of regulatory control achievable in case of different processing mechanisms we changed the values of miRNA-recycling rates κ₁ and κ₂ keeping fixed all the other parameters, obtaining the results displayed in Fig 6. Both *I*<sub>TF</sub> and *I*<sub>miRNA</sub> are higher for small values of the recycling rates, and decrease as miRNA recycling increases. The observed dependence of *I*<sub>miRNA</sub> on κ₁ agrees with the results of [42], where the sensitivity of miRNA molecules to changes in the levels of their targets decreases as the miRNA recycling rate increases, eventually vanishing for κ₁ ≫ σ₁. Indeed, for large κ₁ miRNAs become insensitive to changes in *f*₁, and hence no information can be transmitted.

Remarkably (see Fig 6C), for large κ₂ and small κ₁ miRNA-mediated regulation may be the sole mechanism able to effectively control the target level. In this case, in fact, *I*<sub>TF</sub> becomes especially small as the target is completely repressed by the miRNA while *I*<sub>miRNA</sub> stays finite since the activation of ceRNA₁ can titrate away the miRNA and therefore derepress *m*₂.

As for the previous case, the behaviour of the AOV correlates strongly with this scenario, i.e. larger AOV corresponds to higher achievable information flow in the regulatory element (see S2 Fig). This suggests that, although targets that are degraded purely catalytically cannot compete for miRNA at steady state (as also shown for bacterial small RNAs [61]), they can be notably derepressed if one of their competitors decays stoichiometrically.

**Effective derepression of a repressed target may be conditional on the activation of its competitor and of miRNA.** Transcriptional and post-transcriptional channel capacities are expected to depend on the effective transcription rates of the ceRNAs and the miRNA, described respectively by the quantities *b*ᵢ*n*ᵢ and *β*nᵢ, since these bear a direct impact on the regime (repressed, susceptible or unrepressed) to which each ceRNA belongs. In particular, upon changing the fractional occupancies, situations where miRNA-mediated regulation outperforms TF-based control can occur. Such a possibility might indeed be realized for relatively large values of *k*₁⁺, where the miRNA-mediated channel works optimally. This regime is explored in Fig 7.
<figure><img src="image_10.png" /><figcaption>**Fig 7.** Dependence of Δ*I* on the fractional occupancy of the TF binding site in the optimality range for the post-transcriptional miRNA-mediated channel. (A) Case of targets with weak catalytic degradation (small κ₁ and κ₂). (B) Case of a strongly catalytically degraded target (large κ₂) with a weakly catalytically degraded competitor (small κ₁). Note that *n̄*ᵢ ≈ 1 for (A) while *n̄*ᵢ ≈ 0.5 for (B). Values of the kinetic parameters are reported in Table 1.</figcaption></figure>

doi:10.1371/journal.pcbi.1004715.g007

Remarkably, one finds that post-transcriptional control outperforms transcriptional one when

1. the effective transcription rate of ceRNA₂ exceeds that of ceRNA₁ ($b_2\bar{n}_2 > b_1\bar{n}_1$, see in Table 1) in presence of weak catalytic degradation (Fig 7A);

2. the effective transcription rate of ceRNA₁ is smaller than that of the miRNA ($b_1\bar{n}_1 < \beta\bar{n}_\mu$, see in Table 1) in presence of strongly catalytically degraded target. (see Fig 7B).

A more detailed analysis shows that, in these scenarios, whenever $I_{\text{miRNA}} > I_{\text{TF}}$ the target's derepression size is larger when the signal is processed post-transcriptionally rather than through the direct transcriptional channel, as shown in S3 Fig. Again, the AOV appears therefore to be a good proxy for the channel capacity.

Previous work had pointed to near-equimolarity of ceRNAs and miRNAs within a network as one of the central preconditions for significant ceRNA regulation [34, 42, 62]. On the other hand, experiments on liver cells have shown that such conditions may not be met in physiological conditions, where the number of miRNA target sites can vastly exceed that of miRNA molecules [29, 30]. Our results suggest that target derepression may be significant even if the competitor is in low copy numbers, provided a certain heterogeneity in kinetic parameters (e.g. for a catalytically degraded target and a stoichiometrically degraded competitor) is present, in line with recent experiments pointing to a key role of parameter diverseness [36].

**A sufficiently high degree of target derepression is required for information transmission.** We have seen that AOV displays a dependence on kinetic parameters that very closely matches the one found for the channel capacities. In particular, whenever $I_{\text{miRNA}} > I_{\text{TF}}$, we found $\Delta_{\text{miRNA}} > \Delta_{\text{TF}}$. One may be lead to think that the AOV by itself may fully describe the capacity of miRNA-mediated post-transcriptional information processing and that stochasticity could, to a large degree, be neglected. In some situations, however, the size of the output range is poorly informative about the channel performance. This becomes clear by exploring the relationship between the AOV and the channel capacity at a quantitative level. In Fig 8 the dependence of the optimal MI on the AOV is discussed for two cases: (i) direct transcriptional
<table><caption>**Table 1. Parameters used to obtain each figure.** When the parameter was varied over a range of values, the range is shown in brackets. The parameter $\omega$ appearing in the last column denotes the re-scaling factor defined in Eq (24), while the channel to which parameters refer is specified in brackets. Values of $f_1, f_2$ and $f_\mu$ are mapped to those of $\bar{n}_1, \bar{n}_2, \bar{n}_\mu$ via the steady state conditions, Eq (2).</caption><thead><tr><th>Parameter</th><th>Fig 2A and 2B</th><th>Fig 5A-5C</th><th>Fig 6A-6C</th><th>Fig 7A</th><th>Fig 7B</th><th>Fig 8</th><th>Fig 9</th></tr></thead><tbody><tr><td>b<sub>1</sub> [molecule min<sup>-1</sup>]</td><td>120</td><td>100</td><td>100</td><td>100</td><td>100</td><td>[1, 147] (miRNA); 0 (TF)</td><td>90 (miRNA); 0 (TF)</td></tr><tr><td>b<sub>2</sub> [molecule min<sup>-1</sup>]</td><td>100</td><td>100</td><td>100</td><td>100</td><td>100</td><td>110.2 (miRNA); [10, 98] (TF)</td><td>22.3 (miRNA); 22(TF)</td></tr><tr><td>β [molecule min<sup>-1</sup>]</td><td>100</td><td>100</td><td>100</td><td>100</td><td>100</td><td>100 (miRNA); 0 (TF)</td><td>80 (miRNA); 0(TF)</td></tr><tr><td>d<sub>1</sub> [min<sup>-1</sup>]</td><td>0.1</td><td>0.1</td><td>0.1</td><td>0.1</td><td>0.1</td><td>0.1</td><td>0.1</td></tr><tr><td>d<sub>2</sub> [min<sup>-1</sup>]</td><td>0.1</td><td>0.1</td><td>0.1</td><td>0.1</td><td>0.1</td><td>0.1</td><td>0.1</td></tr><tr><td>δ [min<sup>-1</sup>]</td><td>0.1</td><td>0.1</td><td>0.1</td><td>0.1</td><td>0.1</td><td>0.1</td><td>0.1 ω</td></tr><tr><td>σ<sub>1</sub> [min<sup>-1</sup>]</td><td>0.5</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td></tr><tr><td>σ<sub>2</sub> [min<sup>-1</sup>]</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td></tr><tr><td>k<sub>1</sub><sup>+</sup> [molecule<sup>-1</sup> min<sup>-1</sup>]</td><td>0.002</td><td>[e<sup>-10</sup>, e<sup>2</sup>]</td><td>e<sup>2</sup></td><td>e<sup>2</sup></td><td>e<sup>2</sup></td><td>e<sup>3</sup> (miRNA); 0 (TF)</td><td>e<sup>3.39</sup> ω (miRNA); 0 (TF)</td></tr><tr><td>k<sub>2</sub><sup>+</sup> [molecule<sup>-1</sup> min<sup>-1</sup>]</td><td>0.001</td><td>[e<sup>-10</sup>, e<sup>2</sup>]</td><td>e<sup>-3</sup></td><td>e<sup>-3</sup></td><td>e<sup>-3</sup></td><td>e<sup>0</sup>(miRNA); 0 (TF)</td><td>e<sup>-5.77</sup> ω (miRNA); 0 (TF)</td></tr><tr><td>k<sub>1</sub><sup>-</sup> [min<sup>-1</sup>]</td><td>0.0005</td><td>0.001</td><td>0.001</td><td>0.001</td><td>0.001</td><td>0.001</td><td>0.001</td></tr><tr><td>k<sub>2</sub><sup>-</sup> [min<sup>-1</sup>]</td><td>0.001</td><td>0.001</td><td>0.001</td><td>0.001</td><td>0.001</td><td>0.001</td><td>0.001</td></tr><tr><td>κ<sub>1</sub> [min<sup>-1</sup>]</td><td>0.0003</td><td>0.002</td><td>[e<sup>-5</sup>, e<sup>5</sup>]</td><td>0.002</td><td>e<sup>-8</sup></td><td>0.002</td><td>0.002</td></tr><tr><td>κ<sub>2</sub> [min<sup>-1</sup>]</td><td>0.0004</td><td>0.002</td><td>[e<sup>-5</sup>, e<sup>5</sup>]</td><td>0.002</td><td>e<sup>5</sup></td><td>0.002</td><td>0.002</td></tr><tr><td>n̄<sub>1</sub> [adim.]</td><td>(0,0.99)</td><td>(0,0.99) (miRNA); 0.99(TF)</td><td>(0, 0.99) (miRNA); 0.21(TF)</td><td>(0.1,0.9)</td><td>(0.1,0.9)</td><td>(0, 0.99) (miRNA)</td><td>(0, 0.99) (miRNA)</td></tr><tr><td>n̄<sub>2</sub> [adim.]</td><td>0.71</td><td>0.99 (miRNA); (0,0.99) (TF)</td><td>0.21 (miRNA); (0, 0.99) (TF)</td><td>(0.1,0.9)</td><td>(0.1,0.9)</td><td>0.99 (miRNA); (0, 0.99) (TF)</td><td>0.99 (miRNA); (0, 0.99) (TF)</td></tr><tr><td>n̄<sub>μ</sub> [adim.]</td><td>0.9</td><td>0.99</td><td>0.99</td><td>0.99</td><td>0.5</td><td>0.99</td><td>0.99</td></tr><tr><td>h [adim.]</td><td>5</td><td>5</td><td>5</td><td>5</td><td>5</td><td>5</td><td>5</td></tr><tr><td>k<sub>out</sub>/k<sub>in</sub> [molecules<sup>h</sup>]</td><td>63300</td><td>63300</td><td>63300</td><td>63300</td><td>63300</td><td>63300</td><td>63300</td></tr><tr><td>f<sub>max</sub> [molecules]</td><td>30</td><td>30</td><td>30</td><td>30</td><td>30</td><td>30</td><td>30</td></tr></tbody></table>

doi:10.1371/journal.pcbi.1004715.t001

regulation by TF in absence of miRNA molecules (black curve), (ii) miRNA-mediated post-transcriptional regulation (blue curve).

In Fig 8 we indeed see that for large enough AOV, where fluctuations on each level are small compared to $\Delta_{TF}$, the TF-channel capacity approaches the theoretical limit of a Poissonian channel worked out in S1 Text. On the other hand, for small values of $\Delta_{TF}$, the capacity of TF-channel becomes vanishingly small.

Similarly in Fig 8 one can see that $I_{miRNA}$ increases (roughly) logarithmically with $\Delta_{miRNA}$ if the latter exceeds a certain threshold, whereas it is essentially zero below this value. In other words, the miRNA-mediated channel is unable to convey any information unless the degree of target derepression is sufficiently high. Hence $I_{miRNA}$ can be close to zero even though $\Delta_{miRNA}$ is finite. Notice that the derepression threshold for miRNA-mediated information transmission is larger than the corresponding threshold for direct transcriptional information flow.

A close analysis reveals that such a feature is due to the random fluctuations on the lowest expression level of the target, $m_2^{\min}$. Below threshold, in fact, the degree of derepression of the target is smaller than the fluctuations on $m_2^{\min}$ caused by intrinsic noise (reported in Fig 8).
<figure><img src="image_11.png" /><figcaption>Fig 8. Channel capacities as a function of the target's degree of derepression (AOV). The black curve corresponds to direct transcriptional control in absence of miRNAs, while the blue curve describes the behavior of the miRNA-mediated post-transcriptional channel. For both channels In $m_2^{min} = 4.62$. The predicted maximal MI in the Poissonian limit given in S1 Text is shown as a dashed line. Values of the kinetic parameters are reported in Table 1.</figcaption></figure>

doi:10.1371/journal.pcbi.1004715.g008

Therefore the response observed at the level of the output is indistinguishable from the noise. In the miRNA-mediated post-transcriptional channel the threshold AOV is larger since noise on $m_2^{min}$ is higher. As soon as the AOV overcomes random fluctuations on $m_2^{min}$, the channel capacities starts to increase with $\Delta_{miRNA}$.

**Noise reduction makes the capacity of the post-transcriptional channel comparable to that of the transcriptional channel in the limit of large populations of weakly interacting miRNAs.** We will now focus specifically on the influence of the miRNA-ceRNA binding component of the intrinsic noise, which will be studied upon varying miRNA-ceRNA binding rates and miRNA population size. We shall consider two different scenarios, namely (a) direct, TF-based transcriptional regulation in absence of miRNAs (corresponding to $k_i^+ = 0$ for $i = 1,2$), and (b) indirect post-transcriptional miRNA-mediated regulation (corresponding to non-zero complex binding rates). Kinetic parameters will however be re-scaled as

$$ 
\begin{align} 
\delta &\rightarrow \delta^* &= \omega\delta, \nonumber \\ 
k_i^+ &\rightarrow (k_i^+)^* &= \omega k_i^+, \tag{24}
\end{align} 
$$

where $\omega > 0$ is the re-scaling factor. Note that the above choice does not change $\bar{m}_2$ as long as $k_i^- \bar{c}_i \ll b_i \bar{n}_i$ and $(k_i^- + \kappa_i) \bar{c}_i \ll \beta \bar{n}_\mu$ (both of which are true in the case we consider). By changing $\omega$ one may therefore characterize the dependence of the channel capacities on the kinetic parameters $\delta$ and $k_i^+$ while keeping the output range approximately fixed, such that $\Delta_{miRNA} \approx \Delta_{TF}$ for every $\omega$.

Fig 9A shows how the channel capacities change with $\omega$. One sees that in general $I_{miRNA} < I_{TF}$, so that at fixed AOV the miRNA-mediated channel typically has smaller capacity than the direct one. However, in the limit of small $\omega$, i.e. when miRNA populations are sufficiently large and miRNA-ceRNA couplings are weak, the miRNA channel can perform as effectively as a pure transcriptional channel. This result can be understood by examining the noise levels (see Fig 9B). For small $\omega$ (or high $\mu_{0,2}$), in particular, the FF approaches the Poissonian limit since (i) for $\mu \ge \mu_{02}$ relative fluctuations in miRNA levels are small, extra noise coming from
<figure><img src="image_12.png" /><figcaption>**Fig 9.** Comparison of transcriptional and post-transcriptional regulation for a fixed output variation range. (A) Dependence of $I_{TF}$ and $I_{miRNA}$ on $\omega$. (B) FF of $m_2$ at the corresponding steady state. AOV and $m_2^{min}$ are the same for both channels, namely $\Delta_{miRNA} = \Delta_{TF} = 209 \pm 1$ and $m_2^{min} = 11 \pm 1$. Values of the kinetic parameters are reported in Table 1.</figcaption></figure>

doi:10.1371/journal.pcbi.1004715.g009

molecular titration disappears and the target degradation constant is rescaled as $d_2 \rightarrow d_2 + k_2^+ \bar{\mu}(f_1)$; (ii) $\mu \ll \mu_{0,2}$ target ceRNA is seen in free regime, where it is not sensitive to miRNA molecules. In such conditions, miRNA-mediated noise buffering [20] cannot be observed, nevertheless transcriptional and miRNA-mediated regulation are effectively equivalent.

Observed heterogeneities in miRNA-target interaction strength may
enable regulation by miRNA-mediated cross-talk

The quality of miRNA-target interaction influences the binding kinetics and may be decisive for the activation of the target decay channel [15, 27, 63, 64]. Estimations of the miRNA-ceRNA binding affinities are experimentally challenging. However, computational methods predict a considerable degree of heterogeneity across different miRNA-ceRNA pairs [57]. Remarkably, for the majority of cases reported in the literature, the predicted miRNA-binding energies of the RNAs on the ‘input side’ of the channel are lower than those of the RNAs on the ‘output side’, in line with the optimal conditions highlighted by our model. For example, binding affinities between the long non-coding RNA linc-MD1 and its regulatory miRNAs (playing a central role in skeletal muscle cell differentiation) are significantly lower than those characterizing miRNA interactions with linc-MD1’s competitors, namely the MAML1 and MEF2C mRNAs, as predicted by miRanda algorithm [39, 64]. Likewise, the circular RNA CDR1 has been found to contain around 70 binding sites with high complementarity for miR-7, corresponding to a strong effective coupling through which it can regulate the expression of miR-7’s target genes [37, 65, 66]. Finally, the high sequence homology of pseudogenes (long non-coding RNA genes developed from protein-coding genes but unable to produce proteins) with their parental gene allows them to compete for a large number of shared miRNAs [28, 38, 67].

Our study also points to the potential relevance of miRNA-ceRNA complex decay channels. It is known that, in case of sufficient complementarity, miRNA molecules can function as siRNAs and cleave their targets after binding them [15, 27]. Such targets decay catalytically and are therefore effectively degraded by the miRNAs. Our model predicts that miRNA-mediated control may be the preferred regulatory mechanism in presence of kinetic heterogeneities at the level of miRNA recycling rates. In particular, targets that undergo catalytic degradation
may be efficiently derepressed by their competitors. This type of scenario has been observed in experiments concerning the bacterial small RNA Qrr [61]. Qrr represses its targets by distinct mechanisms. For instance, luxR is repressed catalytically, luxM stoichiometrically, while luxO is silenced through translational repression. luxM and luxO are however able to derepress LuxR in the presence of Qrr. In the light of our results, these observations may therefore point to a higher than expected role for the ceRNA effect in vivo, especially in cases in which heterogeneities in kinetic parameters are thought to be strong [34, 59, 62, 68, 69].

## Discussion

Non-coding RNA molecules, and miRNAs specifically, are increasingly associated to regulatory functions. Besides making it mandatory to characterize the specific role of ncRNAs on a case-by-case basis, especially for situations like disease or differentiation, this fact also raises the question of what ingredients can make miRNAs a preferred tool to regulate the level of a target RNA over, for instance, the target's TF. A possible answer lies in the noise-buffering role that miRNAs can play, which is especially evident in genetic circuitries like incoherent feed-forward loops [19, 20]. By reducing relative fluctuations in the output level, miRNAs can confer robustness to gene expression profiles. However the so-called 'ceRNA hypothesis' opens the way to the possibility that their regulatory functions are carried out at a broader, though more subtle, level. In short, according to the ceRNA scenario miRNAs can mediate an effective positive interaction between their target RNAs driven by the targets' competition to bind them. In this sense, miRNAs can be seen as a sort of 'channel of communication' between RNAs through which RNA levels can be altered and noise can be processed (both buffered and amplified). Previous work [42] has shown that the ceRNA effect may generate both highly plastic and highly selective ceRNA-ceRNA couplings, thereby representing a potentially powerful mechanism to implement gene regulation at the post-transcriptional level.

Although predicted theoretically, the extent and relevance of ceRNA effect *in vivo* is poorly understood. On one hand, considerable evidence points to the ceRNA effect playing a major role in certain dis-regulated or transient cellular states. For instance, it has been shown that the expression of the tumor-suppressor gene PTEN can be regulated by its miRNA-mediated competitors VAPA, CNOT6L, SERINC1 or ZNF460 [70]. Furthermore, many pseudogenes have been found to compete with their parental genes for a shared pool of common microRNAs, thus regulating their expression as competitive endogenous RNA [28, 38, 67, 71]. Such mechanisms seem to be of particular relevance in cancer. For instance, murine models engineered to overexpress the pseudogenes of the proto-oncogene BRAF develop an aggressive malignancy resembling human B cell lymphoma since, by functioning as ceRNAs, they elevate BRAF expression both *in vitro* and *in vivo* [38]. Likewise, the long noncoding RNA linc-MD1 has been shown to regulate the skeletal muscle cell differentiation clock by sponging miRNAs from its competitors, thereby enacting a ceRNA mechanism. In particular, MAML1 and MEF2C (coding for transcription factors that activate muscle-specific gene expression) compete with linc-MD1 for miR-133 and miR-135 respectively [39]. Taken together, the available evidence indicates that miRNA activity depends on the miRNA:target ratio, on miRNA target site abundance and on miRNA binding affinities. Further analyses of high throughput datasets confirm this observation [30, 36]. One may therefore question how said factors may influence miRNA-mediated post-transcriptional control.

The problem however arises of quantifying the degree of control that can be exerted through miRNAs. Taking the 'channel' analogy more strictly (as done before for simpler regulatory elements [9, 45, 52]), one may resort to information theoretic concepts and tools to characterize
precisely how well miRNAs can process fluctuations coming from the modulator nodes and transfer them to the target nodes. This issue goes beyond noise buffering, specifically including the ability to respond to large changes in the mean levels as well as to changes in the structure of fluctuations. As the properties of a channel are conveniently encoded in the mutual information between the input and output nodes, asking how well a channel can function amounts to asking what is the channel's capacity, i.e. the maximum value of the input-output (or modulator-target) mutual information achievable through that channel. This work aimed precisely at quantifying the effectiveness of microRNA-mediated post-transcriptional control of gene expression by computing the capacity of the corresponding regulatory channel and comparing it to that of direct, TF-driven transcriptional regulation. Evidently, multiple factors can influence the flow of information across nodes in a biochemical network, starting from the intrinsic noisiness of each reaction step.

Our basic challenge was therefore understanding in which circumstances miRNA-mediated control can outperform the TF-based one, thereby obtaining insight on why the ceRNA effect appears to be so often employed by cells in situations where accurate tuning and/or shifts of expression levels are required. We have therefore considered, along the lines of [33, 34, 41–44], a mathematical model of the ceRNA effect and characterized its steady state in terms of both mean molecular levels and regulatory capacities of the miRNA-mediated and TF-based channels via stochastic simulations.

We have first considered how the two channels process inputs (the TF levels $f_1$ and $f_2$) that vary in the same range. We have shown that, while the capacity of the TF-channel depends monotonously on each miRNA-ceRNA binding rate and is largest when the target is unrepressed by miRNAs (as might have been expected), the capacity of the post-transcriptional channel is maximal in a specific range of values of the miRNA-ceRNA binding rates. In agreement with [33], we found that miRNA-channel's efficiency is tunable to optimality by the binding kinetics. Furthermore, our model suggests that both capacities decrease as the miRNA recycling rates increase, confirming previous indications obtained by different analytical techniques [42]. Consistently with the scenario observed experimentally for the bacterial small RNA Qrr [61], our model finally suggests that catalytically regulated targets are weakly capable of competing for miRNAs but might be significantly derepressed by their competitors.

In addition, post-transcriptional miRNA-mediated information processing was shown to be characterized by a threshold behaviour as a function of the AOV. In other terms, no information can be transmitted across the channel unless the target's degree of derepression is sufficiently large. This implies that the regulatory effectiveness of the channel is well encoded by the degree of target derepression when the latter is sufficiently high, in which case it is possible to identify regimes in which post-transcriptional regulation is more accurate than transcriptional control.

To get a deeper insight on the origin of the observed threshold behaviour one must however go beyond the AOV and consider more carefully how the miRNA-ceRNA binding noise affects the overall picture. After showing that miRNA-ceRNA binding noise is indeed at the origin of the threshold behaviour that limits the miRNA-channel capacity, we have uncovered the rather remarkable property that in presence of large but weakly interacting miRNA populations the ceRNA effect can regulate gene expression as effectively as the target's modulator node itself.

The present work has focused on a small genetic circuit made up of a single miRNA species and two target RNA species at steady state. Previous work has however shown that cross-talk is possible even during transients [44]. Going beyond stationarity is therefore likely to bring to light new scenarios where miRNA-mediated regulation plays possibly a yet more prominent role.

On the other hand in a typical eukaryotic cell there are thousands ceRNAs, hundreds miRNAs and a rich structure of conserved targeting patterns [72]. Moreover, cells might be
interested in tightly controlling not only each output individually but also particular combina-
tions of output levels (which might be required e.g. for the efficient operation of metabolic
pathways). In such a scenario, miRNA-mediated control could represent a powerful mecha-
nism to increase robustness and flexibility in specific directions of the output space. In view of
this, it would be important to consider a more general multi-source network coding problem
in which a large number of transcription processes are seen as mutually independent informa-
tion sources, and each of the information sources is multicast to sets of output nodes through
the effective network of miRNA-mediated cross-talk interactions. The information-theoretic
scheme employed in this work is easily generalized to deal with more complex networked situa-
tions. Novel insight might finally shed light on the partly controversial picture unveiled by
recent experiments addressing the relevance of the ceRNA effect in vivo [30, 35–39].

**Materials and Methods**

Gillespie algorithm

Numerical simulations have been carried out using the Gillespie algorithm (GA), a standard stochastic method to analyze the time evolution of a system of chemical reactions which is exact for spatially homogeneous systems [73]. In short, based on the reaction rates, GA calculates when the next reaction will occur and what reaction it will be, and then modifies the amount of each molecular species in the system according to the process that took place. If we denote the probability of reaction *r* to occur in the time interval (*τ*, *τ* + *dτ*) by *P(r, τ)*, the algorithm proceeds through the following steps:

1. Initiate the number of reactants in the system and the termination time;

2. Generate a random pair (r, τ) according to P(r, τ);

3. Using the pair (r, τ) just generated, advance time by τ and change number of species involved in reaction r accordingly;

4. Read out the molecular population values. If the termination time is reached, stop the simulations, otherwise return to Step 2.

After a long run, independently on the initial setup, the system of chemical reactants will
come to the equilibrium state.

Linear noise approximation

The mathematical model of ceRNA competition can be solved numerically in the so-called lin-
ear noise approximation. Letting **x** = (*m*<sub>1</sub>, *m*<sub>2</sub>, *µ*, *c*<sub>1</sub>, *c*<sub>2</sub>) stand for the vector of molecular levels,
the kinetic mass action Eq (3) can be re-cast in compact form as

$$
\frac{d\mathbf{x}}{dt} = \mathbf{g}(\mathbf{x}) + \eta, \qquad (25)
$$

where the vector **g** encodes for the deterministic part of the dynamics, while the vector **η** repre-
sents the aggregate noise terms. Each element of **η** has zero mean, and we shall denote its corre-
lations by ⟨**η**<sub>a</sub>(t) **η**<sub>b</sub>(t')⟩ = Γ<sub>ab</sub> δ(t − t'). Denoting by x̄ the steady state, small deviations from it
(i.e. δ**x** = **x**(t) − x̄) relax, in the linear regime, according to

$$
\frac{d}{dt} \delta \mathbf{x} = \mathbf{A} \delta \mathbf{x} + \eta, \qquad (26)
$$

where $\mathbf{A} = \frac{dg}{dx}|_{x=x}$. In this approximation, the correlation matrix $C_{ab} = \langle \delta x_a \delta x_b \rangle$ is given by
[74]

$$C_{ab} = -\sum_{p,q,r,s} B_{ap}B_{br} \frac{\Gamma_{qs}}{\lambda_{p} + \lambda_{r}} (B^{-1})_{pq} (B^{-1})_{rs}, \qquad (27)$$

where $\lambda$'s and $B$'s are, respectively, eigenvalues and eigenvectors of the matrix A.

## Model parameters

Table 1 reports the values of the kinetic parameters (or of their range of variability) employed in the different figures. As $k_1^+$, $k_2^+$, $\kappa_1$, $\kappa_2$, $b_1$, $b_2$, $b_\mu$ and $d_\mu$ are varied in wide ranges in order to test how channel capacities depend on them, our choice was mainly guided by the need to focus the analysis on regimes where ceRNA cross-talk is established (so that the miRNA-mediated channel can actually convey information).

## Supporting Information

**S1 Fig. Information transmission and noise.** Qualitative depiction of how noise affects information flow. Three different situations for the input-output curve $m_2$ (target) vs $f_j$ (modulator) are shown, namely, long (A) target expression level independent of the modulator; (B) derepressed target, with large fluctuations; (C) derepressed target, with small fluctuations. In the first two cases, little or no information can be transmitted from modulator to target (hence the regulatory effectiveness is severely limited), as either the target is insensitive to changes in modulator levels, or its response is strongly hindered by noise. On the other hand, in case (C) information will be transmitted, since modulating the input one can clearly distinguish different output levels. The number of distinguishable levels is linked to the local slope of the input/output curve, as shown mathematically in Eqs (20) and (21) and is mainly limited by the noise strength. In the limit of vanishing noise, when the input-output relationship becomes deterministic, the mutual information between $m_2$ and $f_j$ diverges. (TIF)

**S2 Fig. Dependence of the AOV on kinetic parameters.** (A) $\ln \Delta_{\text{miRNA}}$, (B) $\ln \Delta_{\text{TF}}$, (C) $\ln \Delta_{\text{TF}} - \ln \Delta_{\text{miRNA}}$. Values of the kinetic parameters are as in Fig 5 for the panels in the top row and as in Fig 6 for the panels in the bottom row. (TIF)

**S3 Fig. Dependence of $\ln \Delta_{\text{TF}} - \ln \Delta_{\text{miRNA}}$ on the fractional occupancy of the TF binding site.** Values of the kinetic parameters are as in Fig 7A for panel (A), and as in Fig 7B for panel (B). (TIF)

**S1 Text. Capacity of a Poissonian channel.** (PDF)

**S1 Dataset. Dataset used in the figures.** (ODS)

## Author Contributions

Conceived and designed the experiments: MF EM ADM. Performed the experiments: AM. Analyzed the data: AM MF EM ADM. Contributed reagents/materials/analysis tools: AM MF EM ADM. Wrote the paper: AM MF EM ADM.
## References

1. López-Maury L, Marguerat S, Bähler J. Tuning gene expression to changing environments: from rapid responses to evolutionary adaptation. Nature Reviews Genetics 2008; 9(8): 583–593. doi: 10.1038/ nrg2398 PMID: 18591982

2. Guil S, Esteller M. RNA-RNA interactions in gene regulation: the coding and noncoding players. Trends in Biochemical Sciences 2015; 40(5): 248–256. doi: 10.1016/j.tibs.2015.03.001 PMID: 25818326

3. Cech TR, Steitz JA. The noncoding RNA revolution-trashing old rules to forge new ones. Cell 2014; 157(1): 77–94. doi: 10.1016/j.cell.2014.03.008 PMID: 24679528

4. Latchman DS. Transcription factors: an overview. The International Journal of Biochemistry and Cell Biology 1997; 29(12): 1305–1312. doi: 10.1016/S1357-2725(97)00085-X PMID: 9570129

5. Raser JM, Oheja EK. Noise in gene expression: origins, consequences, and control. Science 2005; 309(5743): 2010–2013. doi: 10.1126/science.1105891 PMID: 16179466

6. Shannon CE. A Mathematical Theory of Communication. The Bell System Technical Journal 1948; 27 (3): 379–423. doi: 10.1002/j.1538-7305.1948.tb01338.x

7. Shannon CE. A Mathematical Theory of Communication. The Bell System Technical Journal 1948; 27 (4): 623–656. doi: 10.1002/j.1538-7305.1948.tb00917.x

8. Tkačik G, Callan CG, Bialek W. Information flow and optimization in transcriptional regulation. PNAS 2008; 105(34): 12265–12270. doi: 10.1073/pnas.0806077105 PMID: 18719112

9. Bialek W. Biophysics: Searching for Principles. Princeton University Press; 2012.

10. Fire A, Xu S, Montgomery MK, Kostas SA, Driver SE, Mello CC. Potent and specific genetic interference by double-stranded RNA in Caenorhabditis elegans. Nature 1998; 391(6669): 806–811. doi: 10.1038/35888 PMID: 9486653

11. Mello CC, Darryl C Jr. Revealing the world of RNA interference. Nature 2004; 431(7006): 338–342. doi: 10.1038/nature02872 PMID: 15372040

12. Baker KE, Parker R. Nonsense-mediated mRNA decay: Terminating erroneous gene expression. Current Opinion in Cell Biology 2004; 16(3): 293–299. doi: 10.1016/j.ceb.2004.03.003 PMID: 15145354

13. Bartel DP. MicroRNAs: genomics, biogenesis, mechanism, and function. Cell 2004; 116(2): 281–297. doi: 10.1016/S0092-8674(04)00045-5 PMID: 14744438

14. Chekulaeva M, Filipowicz W. Mechanisms of miRNA-mediated post-transcriptional regulation in animal cells. Current Opinion in Cell Biology 2009; 21(3): 452–460. doi: 10.1016/j.ceb.2009.04.009 PMID: 19450959

15. Valencia-Sanchez MA, Liu J, Hannon GJ, Parker R. Control of translation and mRNA degradation by miRNAs and siRNAs. Genes and Development 2006; 20(5): 515–524. doi: 10.1101/gad.1399806 PMID: 16510870

16. Hwang H-W, Mendell JT. MicroRNAs in cell proliferation, cell death, and tumorigenesis. British Journal of Cancer 2006; 94(6): 776–780. doi: 10.1038/sj.bjc.6603023 PMID: 16495913

17. Lynn FC. Meta-regulation: microRNA regulation of glucose and lipid metabolism. Trends in Endocrinology and Metabolism 2009; 20(9): 452–459. doi: 10.1016/j.tem.2009.05.007 PMID: 19800254

18. Shi Y, Zhao X, Hsieh J, Wichterle H, Impey S, Banerjee S, et al. MicroRNA Regulation of Neural Stem Cells and Neurogenesis. The Journal of Neuroscience 2010; 30(45): 14931–14936. doi: 10.1523/ JNEUROSCI.4280-10.2010 PMID: 21068294

19. Osella M, Bosia C, Corá D, Caselle M. The role of incoherent microRNA-mediated feedforward loops in noise buffering. PLoS Computational Biology 2011; 7(3): e1001101. doi: 10.1371/journal.pcbi. 1001101 PMID: 21423718

20. Siciliano V, Garzilli I, Fracassi C, Criscuolo S, Ventre S, di Bernardo D. miRNAs confer phenotypic robustness to gene networks by suppressing biological noise. Nature Communications 2013; 4(9): 2364. doi: 10.1038/ncomms3364 PMID: 24077216

21. Leung AK, Sharp PA. MicroRNA functions in stress responses. Molecular Cell 2010; 40(2): 205–215. doi: 10.1016/j.molcel.2010.09.027 PMID: 20965416

22. Lee RC, Feinbaum RL, Ambros V. The C. elegans heterochronic gene lin-4 encodes small RNAs with antisense complementarity to lin-14. Cell 1993; 75(5): 843–854. doi: 10.1016/0092-8674(93)90529-Y PMID: 8252621

23. Lewis BP, Shih I, Jones-Rhoades MW, Bartel DP, Burge CB. Prediction of mammalian microRNA targets. Cell 2003; 115(7): 787–798. doi: 10.1016/S0092-8674(03)01018-3 PMID: 14697198

24. Rhoades MW, Reinhart BJ, Lim LP, Burge CB, Bartel B, Bartel DP. Prediction of plant microRNA targets. Cell 2002; 110(4): 513–520. doi: 10.1016/S0092-8674(02)00863-2 PMID: 12202040
25. Selbach M, Schwanhäusser B, Thierfelder N, Fang Z, Khanin R, Rajewsky N. Widespread changes in protein synthesis induced by microRNAs. Nature 2008; 455(7209): 58–63. doi: 10.1038/nature07228 PMID: 18668040

26. Baek D, Villén J, Shin C, Camargo FD, Gygi SP, Bartel DP. The impact of microRNAs on protein output. Nature 2008; 455(7209): 64–71. doi: 10.1038/nature07242 PMID: 18668037

27. Bartel DP. MicroRNAs: target recognition and regulatory functions. Cell 2009; 136(2): 215–233. doi: 10.1016/j.cell.2009.01.002 PMID: 19167326

28. Salmena L, Poliseno L, Tay Y, Kats L, Pandolfi PP. A ceRNA hypothesis: the Rosetta Stone of a hidden, RNA language? Cell 2011; 146(3): 353–358. doi: 10.1016/j.cell.2011.07.014 PMID: 21802130

29. Jens M, Rajewsky N. Competition between target sites of regulators shapes post-transcriptional gene regulation. Nature Reviews Genetics 2015; 16(2): 113–126. doi: 10.1038/nrg3853 PMID: 25488579

30. Denzler R, Agarwal V, Stefano J, Bartel DP. Stoffel M. Assessing the ceRNA Hypothesis with Quantitative Measurements of miRNA and Target Abundance. Molecular Cell 2014; 54(5): 766–776. doi: 10.1016/j.molcel.2014.03.045 PMID: 24793693

31. Chiu Y-C, Hsiao T-H, Chen Y, Chuang EY. Parameter optimization for constructing competing endogenous RNA regulatory network in glioblastoma multiforme and other cancers. BMC Genomics 2015; 16(4): S1. doi: 10.1186/1471-2164-16-S4-S1 PMID: 25917195

32. Tay Y, Rinn J, Pandolfi PP. The multilayered complexity of ceRNA crosstalk and competition. Nature 2014; 505(7483): 344–352. doi: 10.1038/nature12986 PMID: 24429633

33. Yuan Y, Liu B, Xie P, Zhang MQ, Li Y, Xie Z, et al. Model-guided quantitative analysis of microRNA-mediated regulation on competing endogenous RNAs using a synthetic gene circuit. PNAS 2015; 112(10): 3158–3163. doi: 10.1073/pnas.1413896112 PMID: 25713348

34. Bosia C, Pagnani A, Zecchina R. Modelling Competing Endogenous RNA Networks. PLoS ONE 2013; 8(6): e66609. doi: 10.1371/journal.pone.0066609 PMID: 23840508

35. Ebert MS, Neilson JR, Sharp PA. MicroRNA sponges: Competitive inhibitors of small RNAs in mammalian cells. Nature Methods 2007; 4(9): 721–726. doi: 10.1038/nmeth1079 PMID: 17694064

36. Bosson AD, Zamudio JR, Sharp PA. Endogenous miRNA and target concentrations determine susceptibility to potential ceRNA competition. Molecular Cell 2015; 56(3): 347–359. doi: 10.1016/j.molcel.2014.09.018

37. Memczak S, Jens M, Elefsinioti A, Torti F, Krueger J, Rybak A, et al. Circular RNAs are a large class of animal RNAs with regulatory potency. Nature 2013; 495(7441): 333–338. doi: 10.1038/nature11928 PMID: 23446348

38. Karreth FA, Tay Y, Perna D, Ala U, Tan SM, Rust AG, et al. In vivo identification of tumor-suppressive PTEN ceRNAs in an oncogenic BRAF-induced mouse model of melanoma. Cell 2011; 147(2): 382–395. doi: 10.1016/j.cell.2011.09.032 PMID: 22000016

39. Cesana M, Cacchiarelli D, Legnini I, Santini T, Sthandier O, Chinappi M, et al. A long noncoding RNA controls muscle differentiation by functioning as a competing endogenous RNA. Cell 2011; 147(2): 358–369. doi: 10.1016/j.cell.2011.09.028 PMID: 22000014

40. Levine E, Zhang Z, Kuhlman T, Hwa T. Quantitative characteristics of gene regulation by small RNA. PLoS Biology 2007; 5(9): e229. doi: 10.1371/journal.pbio.0050229 PMID: 17713988

41. Mehta P, Goyal S, Wingreen NS. A quantitative comparison of sRNA-based and protein-based gene regulation. Molecular Systems Biology 2008; 4(1): 221. doi: 10.1038/msb.2008.58 PMID: 18854820

42. Figliuzzi M, Marinari E, De Martino A. MicroRNAs as a selective channel of communication between competing RNAs: a steady-state theory. Biophysical Journal 2013; 104(5): 1203–1213. doi: 10.1016/j.bpj.2013.01.012 PMID: 23473503

43. Bosia C, Sgrò F, Conti L, Baldassi C, Cavallo F, Di Cunto F et al. Quantitative study of crossregulation, noise and synchronization between microRNA targets in single cells. 2015. Preprint. Available: arXiv:1503.06696. Accessed 23 March 2015.

44. Figliuzzi M, De Martino A, Marinari E. RNA-based regulation: dynamics and response to perturbations of competing RNAs. Biophysical journal 2014; 107(4): 1011–1022. doi: 10.1016/j.bpj.2014.06.035 PMID: 25140437

45. Tkačik G, Walczak AM, Bialek W. Optimizing information flow in small genetic networks. Physical Review E 2009; 80(3): 031920. doi: 10.1103/PhysRevE.80.031920

46. Tkačik G, Callan CG Jr, Bialek W. Information Capacity of genetic regulatory networks. Physical Review E 2008; 78(1): 011910. doi: 10.1103/PhysRevE.78.011910

47. Dekel E, Alon U. Optimality and evolutionary tuning of expression level of a protein. Nature 2005; 436(7050): 588–592. doi: 10.1038/nature03842 PMID: 16049495
48. T nase-Nicola S, ten Wolde PR. Regulatory control and the costs and benefit of biochemical noise. PLoS Computational Biology 2008; 4(8): e1000125. doi: 10.1371/journal.pcbi.1000125 PMID: 18716677

49. Muers M. Small RNAs: Recycling for silencing. Nature Reviews Genetics 2011; 12(4): 227. doi: 10.1038/nrg2977 PMID: 21403677

50. Tkačik G, Gregor T, Bialek W. The role of input noise in transcriptional regulation. PLoS One 2007; 3 (7): e2774.

51. Alon U. An Introduction to Systems Biology: Design Principles of Biological Circuits. CRC Press; 2006.

52. Tkačik G, Walczak AM. Information transmission in genetic regulatory networks: a review. Journal of Physics: Condensed Matter 2011; 23(15): 153102. PMID: 21460423

53. van Kampen NG. Stochastic Processes in Physics and Chemistry. North Holland; 1981.

54. Elf J, Paulsson J, Berg OG, Ehrenberg M. Near-critical phenomena in intracellular metabolite pools. Biophysical Journal 2003; 84(1): 154–170. doi: 10.1016/S0006-3495(03)74839-5 PMID: 12524272

55. Noorbakhsh J, Lang AH, Mehta P. Intrinsic Noise of microRNA-Regulated Genes and the ceRNA Hypothesis. PLoS ONE 2013; 8(8): e72676. doi: 10.1371/journal.pone.0072676 PMID: 23991139

56. Rash AM, Winkel BJ. Birth and Death Process Modeling Leads to the Poisson Distribution: A Journey Worth Taking. PRIMUS 2009; 19(1): 57–73. doi: 10.1080/10511970802475157

57. Breda J, Rzepiela AJ, Gumienny R, van Nimwegen E, Zavolan M. Quantifying the strength of miRNA-target interactions. Methods 2015; 85(1): 90–99. doi: 10.1016/j.ymeth.2015.04.012 PMID: 25892562

58. Hao Y, Zhang ZJ, Erickson DW, Huang M, Huang Y, Li J, et al. Quantifying the sequence-function relation in gene silencing by bacterial small RNAs. PNAS 2011; 108(30): 12473–12478. doi: 10.1073/pnas.1100432108 PMID: 21742981

59. Baccarini A, Chauhan H, Gardner TJ, Jayaprakash AD, Sachidanandam R, Brown BD. Kinetic analysis reveals the fate of a microRNA following target regulation in mammalian cells. Current Biology 2011; 21(5): 369–376. doi: 10.1016/j.cub.2011.01.067 PMID: 21353554

60. Hutvágner G, Zamore PD. A microRNA in a multiple-turnover RNAi enzyme complex. Science 2002; 297(5589): 2056–2060. doi: 10.1126/science.1073827 PMID: 12154197

61. Feng L, Rutherford ST, Papenfort K, Bagert JD, van Kessel JC, Tirrell DA, et al. A Qrr-non-coding RNA deploys four different regulatory mechanisms to optimize quorum-sensing dynamics. Cell 2015; 160(1-2): 228–240. doi: 10.1016/j.cell.2014.11.051 PMID: 25579683

62. Ala U, Karreth FA, Bosia C, Pagnani A, Taulli R, Léopold V, et al. Integrated transcriptional and competitive endogenous RNA networks are cross-regulated in permissive molecular environments. PNAS 2013; 110(18): 7154–7159. doi: 10.1073/pnas.1222509110 PMID: 23536298

63. Morozova N, Zinovyev A, Nonne N, Pritchard L-L, Gorban AN, Harel-Bellan A. Kinetic signatures of microRNA modes of action. RNA 2012; 18(9): 1635–1655. doi: 10.1261/rna.032284.112 PMID: 22850425

64. Enright AJ, John B, Gaul U, Tuschl T, Sander C and Marks DS. MicroRNA targets in Drosophila. Genome Biology 2003; 5:R1 doi: 10.1186/gb-2003-5-1-r1 PMID: 14709173

65. Taulli R, Loretelli C, Pandolfi PP, From pseudo-ceRNAs to circ-ceRNAs: a tale of cross-talk and competition. Nature Structural and Molecular Biology 2013; 20(5): 541–543. doi: 10.1038/nsmb.2580 PMID: 23649362

66. Hansen TB, Jensen TI, Clausen BH, Bramsen JB, Finsen B, Damgaard CK, Kjems J. Natural RNA circles function as efficient microRNA sponges. Nature 2013; 495(7441): 384–388. doi: 10.1038/nature11993 PMID: 23446346

67. Poliseno L, Salmena L, Zhang J, Carver B, Haveman WJ, Pandolfi PP. A coding-independent function of gene and pseudogene mRNAs regulates tumour biology. Nature 2010; 465(7301): 1033–1038. doi: 10.1038/nature09144 PMID: 20577206

68. Wee LM, Flores-Jasso CF, Salomon WE, Zamore PD. Argonaute divides its RNA guide into domains with distinct functions and RNA-binding properties. Cell 2012; 151(5): 1055–1067. doi: 10.1016/j.cell.2012.10.036 PMID: 23178124

69. Haley B, Zamore PD. Kinetic analysis of the RNAi enzyme complex. Nature Structural and Molecular Biology 2004; 11(7): 599–606. doi: 10.1038/nsmb780 PMID: 15170178

70. Tay Y, Kats L, Salmena L, Weiss D, Tan SM, Ala U, et al. Coding-Independent Regulation of the Tumor Suppressor PTEN by Competing Endogenous mRNAs. Cell 2011; 147(2), 344–357. doi: 10.1016/j.cell.2011.09.029 PMID: 22000013

71. Kalyana-Sundaram S, Kumar-Sinha C, Shankar S, Robinson DR, Wu Y-M, Cao X. Expressed Pseudo-genes in the Transcriptional Landscape of Human Cancers. Cell 2012; 149(7): 1622–1634. doi: 10.1016/j.cell.2012.04.041 PMID: 22726445
72. Obermayer B, Levine E. Exploring the miRNA Regulatory Network Using Evolutionary Correlations. PLoS Computational Biology 2014; 10(10): e1003860. doi: 10.1371/journal.pcbi.1003860 PMID: 25299225

73. Gillespie DT. Exact stochastic simulation of coupled chemical reactions. Journal of Physical Chemistry 1977; 81(25): 2340–2361. doi: 10.1021/j100540a008

74. Swain PS. Efficient Attenuation of Stochasticity in Gene Expression Through Post-transcriptional Control. Journal of Molecular Biology 2004; 344(4): 965–976. doi: 10.1016/j.jmb.2004.09.073 PMID: 15544806