Mutual interaction in network motifs robustly sharpens gene expression in developmental processes

The gene regulatory network of a developmental process contains many mutually repressive interactions between two genes. They are often regulated by or regulate an additional factor, which constitute prominent network motifs, called regulated and regulating mutual loops. Our database analysis on the gene regulatory network for Drosophila melanogaster indicates that those with mutual repression are working specifically for the segmentation process. To clarify their biological roles, we mathematically study the response of the regulated mutual loop with mutual repression to input stimuli. We show that the mutual repression increases the response sensitivity without affecting the threshold input level to activate the target gene expression, as long as the network output is unique for a given input level. This high sensitivity of the motif can contribute to sharpening the spatial domain pattern without changing its position, assuring a robust developmental process. We also study transient dynamics that shows shift of domain boundary, agreeing with experimental observations. Importance of mutual repression is addressed by comparing with other types of regulations.     

Optimal Feedback Strength for Noise Suppression in Autoregulatory Gene Networks

Autoregulatory feedback loops, where the protein expressed from a gene inhibits or activates its own expression are common gene network motifs within cells. In these networks, stochastic fluctuations in protein levels are attributed to two factors: intrinsic noise (i.e., the randomness associated with mRNA/protein expression and degradation) and extrinsic noise (i.e., the noise caused by fluctuations in cellular components such as enzyme levels and gene-copy numbers). We present results that predict the level of both intrinsic and extrinsic noise in protein numbers as a function of quantities that can be experimentally determined and/or manipulated, such as the response time of the protein and the level of feedback strength. In particular, we show that for a fixed average number of protein molecules, decreasing response times leads to attenuation of both protein intrinsic and extrinsic noise, with the extrinsic noise being more sensitive to changes in the response time. We further show that for autoregulatory networks with negative feedback, the protein noise levels can be minimal at an optimal level of feedback strength. For such cases, we provide an analytical expression for the highest level of noise suppression and the amount of feedback that achieves this minimal noise. These theoretical results are shown to be consistent and explain recent experimental observations. Finally, we illustrate how measuring changes in the protein noise levels as the feedback strength is manipulated can be used to determine the level of extrinsic noise in these gene networks.     

Analytical and Simulation Results for Stochastic Fitzhugh-Nagumo Neurons and Neural Networks

An analytical approach is presented for determining the response of a neuron or of the activity in a network of connected neurons, represented by systems of nonlinear ordinary stochastic differential equations—the Fitzhugh-Nagumo system with Gaussian white noise current. For a single neuron, five equations hold for the first- and second-order central moments of the voltage and recovery variables. From this system we obtain, under certain assumptions, five differential equations for the means, variances, and covariance of the two components. One may use these quantities to estimate the probability that a neuron is emitting an action potential at any given time. The differential equations are solved by numerical methods. We also perform simulations on the stochastic Fitzugh-Nagumo system and compare the results with those obtained from the differential equations for both sustained and intermittent deterministic current inputs withsuperimposed noise. For intermittent currents, which mimic synaptic input, the agreement between the analytical and simulation results for the moments is excellent. For sustained input, the analytical approximations perform well for small noise as there is excellent agreement for the moments. In addition, the probability that a neuron is spiking as obtained from the empirical distribution of the potential in the simulations gives a result almost identical to that obtained using the analytical approach. However, when there is sustained large-amplitude noise, the analytical method is only accurate for short time intervals. Using the simulation method, we study the distribution of the interspike interval directly from simulated sample paths. We confirm that noise extends the range of input currents over which (nonperiodic) spike trains may exist and investigate the dependence of such firing on the magnitude of the mean input current and the noise amplitude. For networks we find the differential equations for the means, variances, and covariances of the voltage and recovery variables and show how solving them leads to an expression for the probability that a given neuron, or given set of neurons, is firing at time t. Using such expressions one may implement dynamical rules for changing synaptic strengths directly without sampling. The present analytical method applies equally well to temporally nonhomogeneous input currents and is expected to be useful for computational studies of information processing in various nervous system centers.     

Exponential Signaling Gain at the Receptor Level Enhances Signal-to-Noise Ratio in Bacterial Chemotaxis

Cellular signaling systems show astonishing precision in their response to external stimuli despite strong fluctuations in the molecular components that determine pathway activity. To control the effects of noise on signaling most efficiently, living cells employ compensatory mechanisms that reach from simple negative feedback loops to robustly designed signaling architectures. Here, we report on a novel control mechanism that allows living cells to keep precision in their signaling characteristics – stationary pathway output, response amplitude, and relaxation time – in the presence of strong intracellular perturbations. The concept relies on the surprising fact that for systems showing perfect adaptation an exponential signal amplification at the receptor level suffices to eliminate slowly varying multiplicative noise. To show this mechanism at work in living systems, we quantified the response dynamics of the E. coli chemotaxis network after genetically perturbing the information flux between upstream and downstream signaling components. We give strong evidence that this signaling system results in dynamic invariance of the activated response regulator against multiplicative intracellular noise. We further demonstrate that for environmental conditions, for which precision in chemosensing is crucial, the invariant response behavior results in highest chemotactic efficiency. Our results resolve several puzzling features of the chemotaxis pathway that are widely conserved across prokaryotes but so far could not be attributed any functional role.     

Multivariate analysis of noise in genetic regulatory networks

Stochasticity is an intrinsic property of genetic regulatory networks due to the low copy numbers of the major molecular species, such as, DNA, mRNA, and regulatory proteins. Therefore, investigation of the mechanisms that reduce the stochastic noise is essential in understanding the reproducible behaviors of real organisms and is also a key to design synthetic genetic regulatory networks that can reliably work. We use an analytical and systematic method, the linear noise approximation of the chemical master equation along with the decoupling of a stoichiometric matrix. In the analysis of fluctuations of multiple molecular species, the covariance is an important measure of noise. However, usually the representation of a covariance matrix in the natural coordinate system, i.e. the copy numbers of the molecular species, is intractably complicated because reactions change copy numbers of more than one molecular species simultaneously. Decoupling of a stoichiometric matrix, which is a transformation of variables, significantly simplifies the representation of a covariance matrix and elucidates the mechanisms behind the observed fluctuations in the copy numbers. We apply our method to three types of fundamental genetic regulatory networks, that is, a single-gene autoregulatory network, a two-gene autoregulatory network, and a mutually repressive network. We have found that there are multiple noise components differently originating. Each noise component produces fluctuation in the characteristic direction. The resulting fluctuations in the copy numbers of the molecular species are the sum of these fluctuations. In the examples, the limitation of the negative feedback in noise reduction and the trade-off of fluctuations in multiple molecular species are clearly explained. The analytical representations show the full parameter dependence. Additionally, the validity of our method is tested by stochastic simulations.     

Noise propagation and signaling sensitivity in biological networks: a role for positive feedback

Interactions between genes and proteins are crucial for efficient processing of internal or external signals, but this connectivity also amplifies stochastic fluctuations by propagating noise between components. Linear (unbranched) cascades were shown to exhibit an interplay between the sensitivity to changes in input signals and the ability to buffer noise. We searched for biological circuits that can maintain signaling sensitivity while minimizing noise propagation, focusing on cases where the noise is characterized by rapid fluctuations. Negative feedback can buffer this type of noise, but this buffering comes at the expense of an even greater reduction in signaling sensitivity. By systematically analyzing three-component circuits, we identify positive feedback as a central motif allowing for the buffering of propagated noise while maintaining sensitivity to long-term changes in input signals. We show analytically that noise reduction in the presence of positive feedback results from improved averaging of rapid fluctuations over time, and discuss in detail a particular implementation in the control of nutrient homeostasis in yeast. As the design of biological networks optimizes for multiple constraints, positive feedback can be used to improve sensitivity without a compromise in the ability to buffer propagated noise.     

Design Principles of a Genetic Alarm Clock

Turning genes on and off is a mechanism by which cells and tissues make phenotypic decisions. Gene network motifs capable of supporting two or more steady states and thereby providing cells with a plurality of possible phenotypes are referred to as genetic switches. Modeled on the bases of naturally occurring genetic networks, synthetic biologists have successfully constructed artificial switches, thus opening a door to new possibilities for improvement of the known, but also the design of new synthetic genetic circuits. One of many obstacles to overcome in such efforts is to understand and hence control intrinsic noise which is inherent in all biological systems. For some motifs the noise is negligible; for others, fluctuations in the particle number can be comparable to its average. Due to their slowed dynamics, motifs with positive autoregulation tend to be highly sensitive to fluctuations of their chemical environment and are in general very noisy, especially during transition (switching). In this article we use stochastic simulations (Gillespie algorithm) to model such a system, in particular a simple bistable motif consisting of a single gene with positive autoregulation. Due to cooperativety, the dynamical behavior of this kind of motif is reminiscent of an alarm clock – the gene is (nearly) silent for some time after it is turned on and becomes active very suddenly. We investigate how these sudden transitions are affected by noise and show that under certain conditions accurate timing can be achieved. We also examine how promoter complexity influences the accuracy of this timing mechanism.     

Brownian dynamics simulation of helix-capping motifs

Helix-capping motifs are believed to play an important role in stabilizing alpha-helices and defining helix start and stop signals. We performed microsecond scale Brownian dynamics simulations to study ten XAAD sequences, with X = (A,E,I,L,N,Q,S,T,V,Y), to examine their propensity to form helix capping motifs and correlate these results with those obtained from analyzing a structural database of proteins. For the widely studied capping box motif S**D, where the asterisk can be any amino acid residue, the simulations suggested that one of the two hydrogen bonds proposed earlier as a stabilizing factor might not be as important. On the other hand, side-chain interactions between the capping residue and the third residue downstream on the polypeptide chain might also play a role in stabilizing this motif. These results are consistent with explicit-solvent molecular dynamics simulations of two capping box motifs found in the proteins BPTI and alpha-dendrotoxin. Principal component analysis of the SAAD trajectory showed that the first three principal components, after those corresponding to translational-rotational motion were removed, accounted for more than half of the conformational fluctuations. The first component separated the conformational space into two parts with the all-helical conformation and the capping box motif lying largely in one part. The second component, on the other hand, could be used to describe conformational transitions between the all-helical form and the capping box motif.     

Compensation for nucleotide bias in a genome by representation as a discrete channel with noise

MOTIVATION: Calculation of the information content of motifs in genomes highly biased in nucleotide composition is likely to lead to overestimates of the amount of useful information in the motif. Calculating relative information can compensate for biases, however the resulting information content is the amount seen by an observer and not by a macromolecule binding to the motif. The latter is needed to calculate the discriminatory power of the motif and to compare motifs between species. RESULTS: By treating a biased genome as a discrete channel with noise, in accordance with Shannon Information Theory, we were able to remove both 'Distortion' and 'Noise' from the motif and recover a more instructive biological 'signal.' A Java application, LogoPaint, was developed to remove nucleotide bias distortion and triplet frequency noise from motifs, calculate information content and present the motif as a logo. We demonstrate how this technique can 'unmask' motifs in the translation initiation regions of bacteria that are obscured by strong sequence biases. AVAILABILITY: LogoPaint is available to all users from the authors as an executable JAR file. Source code is available by arrangement.     

Listening to the noise: random fluctuations reveal gene network parameters

The cellular environment is abuzz with noise originating from the inherent random motion of reacting molecules in the living cell. In this noisy environment, clonal cell populations show cell‐to‐cell variability that can manifest significant phenotypic differences. Noise‐induced stochastic fluctuations in cellular constituents can be measured and their statistics quantified. We show that these random fluctuations carry within them valuable information about the underlying genetic network. Far from being a nuisance, the ever‐present cellular noise acts as a rich source of excitation that, when processed through a gene network, carries its distinctive fingerprint that encodes a wealth of information about that network. We show that in some cases the analysis of these random fluctuations enables the full identification of network parameters, including those that may otherwise be difficult to measure. This establishes a potentially powerful approach for the identification of gene networks and offers a new window into the workings of these networks.     

Assessing the Exceptionality of Coloured Motifs in Networks

Various methods have been recently employed to characterise the structure of biological networks. In particular, the concept of network motif and the related one of coloured motif have proven useful to model the notion of a functional/evolutionary building block. However, algorithms that enumerate all the motifs of a network may produce a very large output, and methods to decide which motifs should be selected for downstream analysis are needed. A widely used method is to assess if the motif is exceptional, that is, over- or under-represented with respect to a null hypothesis. Much effort has been put in the last thirty years to derive -values for the frequencies of topological motifs, that is, fixed subgraphs. They rely either on (compound) Poisson and Gaussian approximations for the motif count distribution in Erdös-Rényi random graphs or on simulations in other models. We focus on a different definition of graph motifs that corresponds to coloured motifs. A coloured motif is a connected subgraph with fixed vertex colours but unspecified topology. Our work is the first analytical attempt to assess the exceptionality of coloured motifs in networks without any simulation. We first establish analytical formulae for the mean and the variance of the count of a coloured motif in an Erdös-Rényi random graph model. Using simulations under this model, we further show that a Pólya-Aeppli distribution better approximates the distribution of the motif count compared to Gaussian or Poisson distributions. The Pólya-Aeppli distribution, and more generally the compound Poisson distributions, are indeed well designed to model counts of clumping events. Altogether, these results enable to derive a -value for a coloured motif, without spending time on simulations.     

Intrinsic noise, gene regulation and steady-state statistics in a two-gene network

The intrinsic noise in a two-gene network model is analysed. The technique of the Fokker-Planck approximation is used to investigate the statistics of noise when the system state is near a stable equilibrium. This is called also the steady-state statistics. The relative size of noise is measured by the Fano factor that is defined as the ratio of the variance to the mean. Our main result shows that in general, the noise control in a two-gene network might be a very complicated process, but for the repressor-repressor system that is a very important case in investigating the genetic switch, the relative size of noise, i.e. the Fano factor, must be bigger than one for both the repressor proteins.