Source author record

Pieter Rein ten Wolde

Pieter Rein ten Wolde appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

35works
10topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

35 published item(s)

preprint2023arXiv

Trade-offs between cost and information in cellular prediction

Living cells can leverage correlations in environmental fluctuations to predict the future environment and mount a response ahead of time. To this end, cells need to encode the past signal into the output of the intracellular network from which the future input is predicted. Yet, storing information is costly while not all features of the past signal are equally informative on the future input signal. Here, we show, for two classes of input signals, that cellular networks can reach the fundamental bound on the predictive information as set by the information extracted from the past signal: push-pull networks can reach this information bound for Markovian signals, while networks that take a temporal derivative can reach the bound for predicting the future derivative of non-Markovian signals. However, the bits of past information that are most informative about the future signal are also prohibitively costly. As a result, the optimal system that maximizes the predictive information for a given resource cost is, in general, not at the information bound. Applying our theory to the chemotaxis network of Escherichia coli reveals that its adaptive kernel is optimal for predicting future concentration changes over a broad range of background concentrations, and that the system has been tailored to predicting these changes in shallow gradients.

preprint2021arXiv

Energetic constraints on filament mediated cell polarization

Cell polarization underlies many cellular processes, such as differentiation, migration, and budding. Many living cells, such as budding yeast and fission yeast, use cytoskeletal structures to actively transport proteins to one location on the membrane and create a high density spot of membrane-bound proteins. Yet, the thermodynamic constraints on filament-based cell polarization remain unknown. We show by mathematical modeling that cell polarization requires detailed balance to be broken, and we quantify the free-energy cost of maintaining a polarized state of the cell. Our study reveals that detailed balance cannot only be broken via the active transport of proteins along filaments, but also via a chemical modification cycle, allowing detailed balance to be broken by the shuttling of proteins between the filament, membrane, and cytosol. Our model thus shows that cell polarization can be established via two distinct driving mechanisms, one based on active transport and one based on non-equilibrium binding. Furthermore, the model predicts that the driven binding process dissipates orders of magnitude less free-energy than the transport-based process to create the same membrane spot. Active transport along filaments may be sufficient to create a polarized distribution of membrane-bound proteins, but an additional chemical modification cycle of the proteins themselves is more efficient and less sensitive to the physical exclusion of proteins on the transporting filaments, providing insight in the design principles of the Pom1/Tea1/Tea4 system in fission yeast and the Cdc42 system in budding yeast.

preprint2019arXiv

Diffusible crosslinkers cause superexponential friction forces

The mitotic spindle lies at the heart of the spatio-temporal control over cellular components during cell division. The spindle consists of microtubules, which are not only crosslinked by motor proteins but also by passive binding proteins. These passive crosslinkers stabilize the highly dynamic mitotic spindle by generating friction forces between sliding filaments. However, it remains unclear how the friction coefficient depends on the number of crosslinkers and the size of the overlap between the microtubules. Here, we use theory and computer simulations to study the friction between two filaments that are crosslinked by passive proteins, which can hop between neighboring binding sites while physically excluding each other. The simulations reveal that the movement of one microtubule relative to the other is limited by free-energy barrier crossings, causing rare and discrete jumps of the microtubule that span the distance between adjacent crosslinker binding sites. We derive an exact analytical expression for the free-energy landscape and identify the reaction coordinate that governs the relative movement, which allows us to determine the effective barrier height for the microtubule jumps. Both through simulations and reaction rate theory, we make the experimentally testable prediction that the friction between the microtubules increases superexponentially with the density of crosslinkers.

preprint2019arXiv

The behaviour of information flow near criticality

Recent experiments have indicated that many biological systems self-organise near their critical point, which hints at a common design principle. While it has been suggested that information transmission is optimized near the critical point, it remains unclear how information transmission depends on the dynamics of the input signal, the distance over which the information needs to be transmitted, and the distance to the critical point. Here we employ stochastic simulations of a driven 2D Ising system and study the instantaneous mutual information and the information transmission rate between a driven input spin and an output spin. The instantaneous mutual information varies non-monotonically with the temperature, but increases monotonically with the correlation time of the input signal. In contrast, the information transmission rate exhibits a maximum as a function of the input correlation time. Moreover, there exists an optimal temperature that maximizes this maximum information transmission rate. It arises from a tradeoff between the necessity to respond fast to changes in the input so that more information per unit amount of time can be transmitted, and the need to respond to reliably. The optimal temperature lies above the critical point, but moves towards it as the distance between the input and output spin is increased.

preprint2018arXiv

Non-equilibrium correlations in minimal dynamical models of polymer copying

Living systems produce "persistent" copies of information-carrying polymers, in which template and copy sequences remain correlated after physically decoupling. We identify a general measure of the thermodynamic efficiency with which these non-equilibrium states are created, and analyze the accuracy and efficiency of a family of dynamical models that produce persistent copies. For the weakest chemical driving, when polymer growth occurs in equilibrium, both the copy accuracy and, more surprisingly, the efficiency vanish. At higher driving strengths, accuracy and efficiency both increase, with efficiency showing one or more peaks at moderate driving. Correlations generated within the copy sequence, as well as between template and copy, store additional free energy in the copied polymer and limit the single-site accuracy for a given chemical work input. Our results provide insight in the design of natural self-replicating systems and can aid the design of synthetic replicators.

preprint2016arXiv

The accuracy of telling time via oscillatory signals

Circadian clocks are the central timekeepers of life, allowing cells to anticipate changes between day and night. Experiments in recent years have revealed that circadian clocks can be highly stable, raising the question how reliably they can be read out. Here, we combine mathematical modeling with information theory to address the question how accurately a cell can infer the time from an ensemble of protein oscillations, which are driven by a circadian clock. We show that the precision increases with the number of oscillations and their amplitude relative to their noise. Our analysis also reveals that their exists an optimal phase relation that minimizes the error in the estimate of time, which depends on the relative noise levels of the protein oscillations. Lastly, our work shows that cross-correlations in the noise of the protein oscillations can enhance the mutual information, which suggests that cross-regulatory interactions between the proteins that read out the clock can be beneficial for temporal information transmission.

preprint2015arXiv

Discrete gene replication events drive coupling between the cell cycle and circadian clocks

Many organisms possess both a cell cycle to control DNA replication and a circadian clock to anticipate changes between day and night. In some cases, these two rhythmic systems are known to be coupled by specific, cross-regulatory interactions. Here, we use mathematical modeling to show that, additionally, the cell cycle generically influences circadian clocks in a non-specific fashion: The regular, discrete jumps in gene-copy number arising from DNA replication during the cell cycle cause a periodic driving of the circadian clock, which can dramatically alter its behavior and impair its function. A clock built on negative transcriptional feedback either phase locks to the cell cycle, so that the clock period tracks the cell division time, or exhibits erratic behavior. We argue that the cyanobacterium Synechococcus elongatus has evolved two features that protect its clock from such disturbances, both of which are needed to fully insulate it from the cell cycle and give it its observed robustness: a phosphorylation-based protein modification oscillator, together with its accompanying push-pull read-out circuit that responds primarily to the ratios of the different phosphoforms, makes the clock less susceptible to perturbations in protein synthesis; and the presence of multiple, asynchronously replicating copies of the same chromosome diminishes the effect of replicating any single copy of a gene.

preprint2015arXiv

Fundamental Limits to Cellular Sensing

In recent years experiments have demonstrated that living cells can measure low chemical concentrations with high precision, and much progress has been made in understanding what sets the fundamental limit to the precision of chemical sensing. Chemical concentration measurements start with the binding of ligand molecules to receptor proteins, which is an inherently noisy process, especially at low concentrations. The signaling networks that transmit the information on the ligand concentration from the receptors into the cell have to filter this noise extrinsic to the cell as much as possible. These networks, however, are also stochastic in nature, which means that they will also add noise to the transmitted signal. In this review, we will first discuss how the diffusive transport and binding of ligand to the receptor sets the receptor correlation time, and then how downstream signaling pathways integrate the noise in the receptor state; we will discuss how the number of receptors, the receptor correlation time, and the effective integration time together set a fundamental limit on the precision of sensing. We then discuss how cells can remove the receptor noise while simultaneously suppressing the intrinsic noise in the signaling network. We describe why this mechanism of time integration requires three classes of resources---receptors and their integration time, readout molecules, energy---and how each resource class sets a fundamental sensing limit. We also briefly discuss the scheme of maximum-likelihood estimation, the role of receptor cooperativity, and how cellular copy protocols differ from canonical copy protocols typically considered in the computational literature, explaining why cellular sensing systems can never reach the Landauer limit on the optimal trade-off between accuracy and energetic cost.

preprint2015arXiv

Optimal Prediction by Cellular Signaling Networks

Living cells can enhance their fitness by anticipating environmental change. We study how accurately linear signaling networks in cells can predict future signals. We find that maximal predictive power results from a combination of input-noise suppression, linear extrapolation, and selective readout of correlated past signal values. Single-layer networks generate exponential response kernels, which suffice to predict Markovian signals optimally. Multilayer networks allow oscillatory kernels that can optimally predict non-Markovian signals. At low noise, these kernels exploit the signal derivative for extrapolation, while at high noise, they capitalize on signal values in the past that are strongly correlated with the future signal. We show how the common motifs of negative feedback and incoherent feed-forward can implement these optimal response functions. Simulations reveal that E. coli can reliably predict concentration changes for chemotaxis, and that the integration time of its response kernel arises from a trade-off between rapid response and noise suppression.

preprint2014arXiv

How biochemical resources determine fundamental limits in cellular sensing

Living cells deploy many resources to sense their environments, including receptors, downstream signaling molecules, time and fuel. However, it is not known which resources fundamentally limit the precision of sensing, like weak links in a chain, and which can compensate each other, leading to trade-offs between them. We show by modeling that in equilibrium systems the precision is limited by the number of receptors; the downstream network can never increase precision. This limit arises from a trade-off between the removal of extrinsic noise in the receptor and intrinsic noise in the downstream network. Non-equilibrium systems can lift this trade-off by storing the receptor state over time in chemical modification states of downstream molecules. As we quantify for a push-pull network, this requires i) time and receptors; ii) downstream molecules; iii) energy (fuel turnover) to drive modification. These three resource classes cannot compensate each other, and it is the limiting class which sets the fundamental sensing limit. Within each class, trade-offs are possible. Energy allows a power-speed trade-off, while time can be traded against receptors.

preprint2014arXiv

The lower bound on the precision of transcriptional regulation

The diffusive arrival of transcription factors at the promoter sites on the DNA sets a lower bound on how accurately a cell can regulate its protein levels. Using results from the literature on diffusion-influenced reactions, we derive an analytical expression for the lower bound on the precision of transcriptional regulation. In our theory, transcription factors can perform multiple rounds of 1D diffusion along the DNA and 3D diffusion in the cytoplasm before binding to the promoter. Comparing our expression for the lower bound on the precision against results from Green's Function Reaction Dynamics simulations shows that the theory is highly accurate under biologically relevant conditions. Our results demonstrate that, to an excellent approximation, the promoter switches between the transcription-factor bound and unbound state in a Markovian fashion. This remains true even in the presence of sliding, i.e. with 1D diffusion along the DNA. This has two important implications: (1) minimizing the noise in the promoter state is equivalent to minimizing the search time of transcription factors for their promoters; (2) the complicated dynamics of 3D diffusion in the cytoplasm and 1D diffusion along the DNA can be captured in a well-stirred model by renormalizing the promoter association and dissociation rates, making it possible to efficiently simulate the promoter dynamics using Gillespie simulations. Based on the recent experimental observation that sliding can speed up the promoter search by a factor of 4, our theory predicts that sliding can enhance the precision of transcriptional regulation by a factor of 2.

preprint2014arXiv

The robustness of proofreading to crowding-induced pseudo-processivity in the MAPK pathway

Double phosphorylation of protein kinases is a common feature of signalling cascades. This motif may reduce cross-talk between signalling pathways, as the second phosphorylation site allows for proofreading, especially when phosphorylation is distributive rather than processive. Recent studies suggest that phosphorylation can be `pseudo-processive' in the crowded cellular environment, as rebinding after the first phosphorylation is enhanced by slow diffusion. Here, we use a simple model with unsaturated reactants to show that specificity for one substrate over another drops as rebinding increases and pseudo-processive behavior becomes possible. However, this loss of specificity with increased rebinding is typically also observed if two distinct enzyme species are required for phosphorylation, i.e. when the system is necessarily distributive. Thus the loss of specificity is due to an intrinsic reduction in selectivity with increased rebinding, which benefits inefficient reactions, rather than pseudo-processivity itself. We also show that proofreading can remain effective when the intended signalling pathway exhibits high levels of rebinding-induced pseudo-processivity, unlike other proposed advantages of the dual phosphorylation motif.

preprint2013arXiv

Multiplexing oscillatory biochemical signals

In recent years it is increasingly being recognized that biochemical signals are not necessarily constant in time and that the temporal dynamics of a signal can be the information carrier. Moreover, it is now well established that components are often shared between signaling pathways. Here we show by mathematical modeling that living cells can multiplex a constant and an oscillatory signal: they can transmit these two signals through the same signaling pathway simultaneously, and yet respond to them specifically and reliably. We find that information transmission is reduced not only by noise arising from the intrinsic stochasticity of biochemical reactions, but also by crosstalk between the different channels. Yet, under biologically relevant conditions more than 2 bits of information can be transmitted per channel, even when the two signals are transmitted simultaneously. These observations suggest that oscillatory signals are ideal for multiplexing signals.

preprint2013arXiv

Prediction and Dissipation in Biochemical Sensing

Cells sense and predict their environment via energy-dissipating pathways. However, it is unclear whether dissipation helps or harms prediction. Here we study dissipation and prediction for a minimal sensory module of receptors that reversibly bind ligand. We find that the module performs short-term prediction optimally when operating in an adiabatic regime where dissipation vanishes. In contrast, beyond a critical forecast interval, prediction becomes most precise in a regime of maximal dissipation, suggesting that dissipative sensing in biological systems can serve to enhance prediction performance.

preprint2013arXiv

Spatial partitioning improves the reliability of biochemical signaling

Spatial heterogeneity is a hallmark of living systems, even at the molecular scale in individual cells. A key example is the partitioning of membrane-bound proteins via lipid domain formation or cytoskeleton-induced corralling. Yet the impact of this spatial heterogeneity on biochemical signaling processes is poorly understood. Here we demonstrate that partitioning improves the reliability of biochemical signaling. We exactly solve a stochastic model describing a ubiquitous motif in membrane signaling. The solution reveals that partitioning improves signaling reliability via two effects: it moderates the non-linearity of the switching response, and it reduces noise in the response by suppressing correlations between molecules. An optimal partition size arises from a trade-off between minimizing the number of proteins per partition to improve signaling reliability and ensuring sufficient proteins per partition to maintain signal propagation. The predicted optimal partition size agrees quantitatively with experimentally observed systems. These results persist in spatial simulations with explicit diffusion barriers. Our findings suggest that molecular partitioning is not merely a consequence of the complexity of cellular substructures, but also plays an important functional role in cell signaling.

preprint2013arXiv

The importance of crowding in signaling, genetic, and metabolic networks

It is now well established that the cell is a highly crowded environment. Yet, the effects of crowding on the dynamics of signaling pathways, gene regulation networks and metabolic networks are still largely unknown. Crowding can alter both molecular diffusion and the equilibria of biomolecular reactions. In this review, we first discuss how diffusion can affect biochemical networks. Diffusion of transcription factors can increase noise in gene expression, while diffusion of proteins between intracellular compartments or between cells can reduce concentration fluctuations. In push-pull networks diffusion can impede information transmission, while in multi-site protein modification networks diffusion can qualitatively change the macroscopic response of the system, such as the loss or emergence of bistability. Moreover, diffusion can directly change the metabolic flux. We describe how crowding affects diffusion, and thus how all these phenomena are influenced by crowding. Yet, a potentially more important effect of crowding on biochemical networks is mediated via the shift in the equilibria of bimolecular reactions, and we provide computational evidence that supports this idea. Finally, we discuss how the effects of crowding can be incorporated in models of biochemical networks.

preprint2012arXiv

Fundamental Limits on Sensing Chemical Concentrations with Linear Biochemical Networks

Living cells often need to extract information from biochemical signals that are noisy. We study how accurately cells can measure chemical concentrations with signaling networks that are linear. For stationary signals of long duration, they can reach, but not beat, the Berg-Purcell limit, which relies on uniformly averaging in time the fluctuations in the input signal. For short times or nonstationary signals, however, they can beat the Berg-Purcell limit, by non-uniformly time-averaging the input. We derive the optimal weighting function for time averaging and use it to provide the fundamental limit of measuring chemical concentrations with linear signaling networks.

preprint2012arXiv

Mutual Repression enhances the Steepness and Precision of Gene Expression Boundaries

Embryonic development is driven by spatial patterns of gene expression that determine the fate of each cell in the embryo. While gene expression is often highly erratic, embryonic development is usually exceedingly precise. In particular, gene expression boundaries are robust not only against intrinsic noise from gene expression and protein diffusion, but also against embryo-to-embryo variations in the morphogen gradients, which provide positional information to the differentiating cells. How development is robust against intra- and inter-embryonic variations is not understood. A common motif in the gene regulation networks that control embryonic development is mutual repression between pairs of genes. To assess the role of mutual repression in the robust formation of gene expression patterns, we have performed large-scale stochastic simulations of a minimal model of two mutually repressing gap genes in Drosophila, hunchback (hb) and knirps (kni). Our model includes not only mutual repression between hb and kni, but also the stochastic and cooperative activation of hb by the anterior morphogen Bicoid (Bcd) and of kni by the posterior morphogen Caudal (Cad), as well as the diffusion of Hb and Kni. Our analysis reveals that mutual repression can markedly increase the steepness and precision of the gap gene expression boundaries. In contrast to spatial averaging and cooperative gene activation, mutual repression thus allows for gene-expression boundaries that are both steep and precise. Moreover, mutual repression dramatically enhances their robustness against embryo-to-embryo variations in the morphogen levels. Finally, our simulations reveal that gap protein diffusion plays a critical role not only in reducing the width of gap gene expression boundaries via spatial averaging, but also in repairing patterning errors that could arise due to the bistability induced by mutual repression.

preprint2012arXiv

Non-Stationary Forward Flux Sampling

We present a new method, Non-Stationary Forward Flux Sampling, that allows efficient simulation of rare events in both stationary and non-stationary stochastic systems. The method uses stochastic branching and pruning to achieve uniform sampling of trajectories in phase space and time, leading to accurate estimates for time-dependent switching propensities and time-dependent phase space probability densities. The method is suitable for equilibrium or non-equilibrium systems, in or out of stationary state, including non-Markovian or externally driven systems. We demonstrate the validity of the technique by applying it to a one-dimensional barrier crossing problem that can be solved exactly, and show its usefulness by applying it to the time-dependent switching of a genetic toggle switch.

preprint2012arXiv

Protein logic: a statistical mechanical study of signal integration at the single-molecule level

Information processing and decision making is based upon logic operations, which in cellular networks has been well characterized at the level of transcription. In recent years however, both experimentalists and theorists have begun to appreciate that cellular decision making can also be performed at the level of a single protein, giving rise to the notion of protein logic. Here we systematically explore protein logic using a well known statistical mechanical model. As an example system, we focus on receptors which bind either one or two ligands, and their associated dimers. Notably, we find that a single heterodimer can realize any of the 16 possible logic gates, including the XOR gate, by variation of biochemical parameters. We then introduce the novel idea that a set of receptors with fixed parameters can encode functionally unique logic gates simply by forming different dimeric combinations. An exhaustive search reveals that the simplest set of receptors (two single-ligand receptors and one double-ligand receptor) can realize several different groups of three unique gates, a result for which the parametric analysis of single receptors and dimers provides a clear interpretation. Both results underscore the surprising functional freedom readily available to cells at the single-protein level.

preprint2012arXiv

Rare switching events in non-stationary systems

Physical systems with many degrees of freedom can often be understood in terms of transitions between a small number of metastable states. For time-homogeneous systems with short-term memory these transitions are fully characterized by a set of rate constants. We consider the question how to extend such a coarse-grained description to non-stationary systems and to systems with finite memory. We identify the physical regimes in which time-dependent rates are meaningful, and state microscopic expressions that can be used to measure both externally time-dependent and history-dependent rates in microscopic simulations.

preprint2012arXiv

Reliability of frequency- and amplitude-decoding in gene regulation

In biochemical signaling, information is often encoded in oscillatory signals. However, the advantages of such a coding strategy over an amplitude encoding scheme of constant signals remain unclear. Here we study the dynamics of a simple model gene promoter in response to oscillating and constant transcription factor signals. We find that in biologically-relevant parameter regimes an oscillating input can produce a more constant protein level than a constant input. Our results suggest that oscillating signals may be used to minimize noise in gene regulation.

preprint2012arXiv

Signalling noise enhances chemotactic drift of E. coli

Noise in transduction of chemotactic stimuli to the flagellar motor of E. coli will affect the random run-and-tumble motion of the cell and the ability to perform chemotaxis. Here we use numerical simulations to show that an intermediate level of noise in the slow methylation dynamics enhances drift while not compromising localisation near concentration peaks. A minimal model shows how such an optimal noise level arises from the interplay of noise and the dependence of the motor response on the network output. Our results suggest that cells can exploit noise to improve chemotactic performance.

preprint2012arXiv

The macroscopic effects of microscopic heterogeneity

Over the past decade, advances in super-resolution microscopy and particle-based modeling have driven an intense interest in investigating spatial heterogeneity at the level of single molecules in cells. Remarkably, it is becoming clear that spatiotemporal correlations between just a few molecules can have profound effects on the signaling behavior of the entire cell. While such correlations are often explicitly imposed by molecular structures such as rafts, clusters, or scaffolds, they also arise intrinsically, due strictly to the small numbers of molecules involved, the finite speed of diffusion, and the effects of macromolecular crowding. In this chapter we review examples of both explicitly imposed and intrinsic correlations, focusing on the mechanisms by which microscopic heterogeneity is amplified to macroscopic effect.

preprint2011arXiv

Membrane clustering and the role of rebinding in biochemical signaling

In many cellular signaling pathways, key components form clusters at the cell membrane. Although much work has focused on the mechanisms behind such cluster formation, the implications for downstream signaling remain poorly understood. Here, motivated by recent experiments, we study via particle-based simulation a covalent modification network in which the activating component is either clustered or randomly distributed on the membrane. We find that while clustering reduces the response of a single-modification network, clustering can enhance the response of a double-modification network. The reduction is a bulk effect: a cluster presents a smaller effective target to a substrate molecule in the bulk. The enhancement, on the other hand, is a local effect: a cluster promotes the rapid rebinding and second activation of singly active substrate molecules. As such, the enhancement relies upon frequent collisions on a short timescale, which leads to a diffusion coefficient at which the enhancement is optimal. We complement simulation with analytic results at both the mean-field and first-passage distribution levels. Our results emphasize the importance of spatially resolved models, showing that significant effects of spatial correlations persist even in spatially averaged quantities such as response curves.

preprint2010arXiv

Multiplexing Biochemical Signals

In this paper we show that living cells can multiplex biochemical signals, i.e. transmit multiple signals through the same signaling pathway simultaneously, and yet respond to them very specifically. We demonstrate how two binary input signals can be encoded in the concentration of a common signaling protein, which is then decoded such that each of the two output signals provides reliable information about one corresponding input. Under biologically relevant conditions the network can reach the maximum amount of information that can be transmitted, which is 2 bits.

preprint2010arXiv

Mutual information in time-varying biochemical systems

Cells must continuously sense and respond to time-varying environmental stimuli. These signals are transmitted and processed by biochemical signalling networks. However, the biochemical reactions making up these networks are intrinsically noisy, which limits the reliability of intracellular signalling. Here we use information theory to characterise the reliability of transmission of time-varying signals through elementary biochemical reactions in the presence of noise. We calculate the mutual information for both instantaneous measurements and trajectories of biochemical systems for a Gaussian model. Our results indicate that the same network can have radically different characteristics for the transmission of instantaneous signals and trajectories. For trajectories, the ability of a network to respond to changes in the input signal is determined by the timing of reaction events, and is independent of the correlation time of the output of the network. We also study how reliably signals on different time-scales can be transmitted by considering the frequency-dependent coherence and gain-to-noise ratio. We find that a detector that does not consume the ligand molecule upon detection can more reliably transmit slowly varying signals, while an absorbing detector can more reliably transmit rapidly varying signals. Furthermore, we find that while one reaction may more reliably transmit information than another when considered in isolation, when placed within a signalling cascade the relative performance of the two reactions can be reversed. This means that optimising signal transmission at a single level of a signalling cascade can reduce signalling performance for the cascade as a whole.

preprint2010arXiv

Robust circadian clocks from coupled protein modification and transcription-translation cycles

The cyanobacterium Synechococcus elongatus uses both a protein phosphorylation cycle and a transcription-translation cycle to generate circadian rhythms that are highly robust against biochemical noise. We use stochastic simulations to analyze how these cycles interact to generate stable rhythms in growing, dividing cells. We find that a protein phosphorylation cycle by itself is robust when protein turnover is low. For high decay or dilution rates (and co mpensating synthesis rate), however, the phosphorylation-based oscillator loses its integrity. Circadian rhythms thus cannot be generated with a phosphorylation cycle alone when the growth rate, and consequently the rate of protein dilution, is high enough; in practice, a purely post-translational clock ceases to function well when the cell doubling time drops below the 24 hour clock period. At higher growth rates, a transcription-translation cycle becomes essential for generating robust circadian rhythms. Interestingly, while a transcription-translation cycle is necessary to sustain a phosphorylation cycle at high growth rates, a phosphorylation cycle can dramatically enhance the robustness of a transcription-translation cycle at lower protein decay or dilution rates. Our analysis thus predicts that both cycles are required to generate robust circadian rhythms over the full range of growth conditions.

preprint2010arXiv

The effect of feedback on the fidelity of information transmission of time-varying signals

Living cells are continually exposed to environmental signals that vary in time. These signals are detected and processed by biochemical networks, which are often highly stochastic. To understand how cells cope with a fluctuating environment, we therefore have to understand how reliably biochemical networks can transmit time-varying signals. To this end, we must understand both the noise characteristics and the amplification properties of networks. In this manuscript, we use information theory to study how reliably signalling cascades employing autoregulation and feedback can transmit time-varying signals. We calculate the frequency-dependence of the gain-to-noise ratio, which reflects how reliably a network transmits signals at different frequencies. We find that the gain-to-noise ratio may differ qualitatively from the power spectrum of the output, showing that the latter does not directly reflect signaling performance. Moreover, we find that auto-activation and auto-repression increase and decrease the gain-to-noise ratio for all of frequencies, respectively. Positive feedback specifically enhances information transmission at low frequencies, while negative feedback increases signal fidelity at high frequencies. Our analysis not only elucidates the role of autoregulation and feedback in naturally-occurring biological networks, but also reveals design principles that can be used for the reliable transmission of time-varying signals in synthetic gene circuits.

preprint2009arXiv

DNA looping provides stability and robustness to the bacteriophage lambda switch

The bistable gene regulatory switch controlling the transition from lysogeny to lysis in bacteriophage lambda presents a unique challenge to quantitative modeling. Despite extensive characterization of this regulatory network, the origin of the extreme stability of the lysogenic state remains unclear. We have constructed a stochastic model for this switch. Using Forward Flux Sampling simulations, we show that this model predicts an extremely low rate of spontaneous prophage induction in a recA mutant, in agreement with experimental observations. In our model, the DNA loop formed by octamerization of CI bound to the O_L and O_R operator regions is crucial for stability, allowing the lysogenic state to remain stable even when a large fraction of the total CI is depleted by nonspecific binding to genomic DNA. DNA looping also ensures that the switch is robust to mutations in the order of the O_R binding sites. Our results suggest that DNA looping can provide a mechanism to maintain a stable lysogenic state in the face of a range of challenges including noisy gene expression, nonspecific DNA binding and operator site mutations.

preprint2009arXiv

Forward Flux Sampling for rare event simulations

Rare events are ubiquitous in many different fields, yet they are notoriously difficult to simulate because few, if any, events are observed in a conventiona l simulation run. Over the past several decades, specialised simulation methods have been developed to overcome this problem. We review one recently-developed class of such methods, known as Forward Flux Sampling. Forward Flux Sampling uses a series of interfaces between the initial and final states to calculate rate constants and generate transition paths, for rare events in equilibrium or nonequilibrium systems with stochastic dynamics. This review draws together a number of recent advances, summarizes several applications of the method and highlights challenges that remain to be overcome.

preprint2009arXiv

Mutual information between in- and output trajectories of biochemical networks

Biochemical networks can respond to temporal characteristics of time-varying signals. To understand how reliably biochemical networks can transmit information we must consider how an input signal as a function of time--the input trajectory--can be mapped onto an output trajectory. Here we estimate the mutual information between in- and output trajectories using a Gaussian model. We study how reliably the chemotaxis network of E. coli can transmit information on the ligand concentration to the flagellar motor, and find the input power spectrum that maximizes the information transmission rate.

preprint2009arXiv

Spatio-temporal correlations can drastically change the response of a MAPK pathway

Multisite covalent modification of proteins is omnipresent in eukaryotic cells. A well-known example is the mitogen-activated protein kinase (MAPK) cascade, where in each layer of the cascade a protein is phosphorylated at two sites. It has long been known that the response of a MAPK pathway strongly depends on whether the enzymes that modify the protein act processively or distributively: distributive mechanism, in which the enzyme molecules have to release the substrate molecules in between the modification of the two sites, can generate an ultrasensitive response and lead to hysteresis and bistability. We study by Green's Function Reaction Dynamics, a stochastic scheme that makes it possible to simulate biochemical networks at the particle level and in time and space, a dual phosphorylation cycle in which the enzymes act according to a distributive mechanism. We find that the response of this network can differ dramatically from that predicted by a mean-field analysis based on the chemical rate equations. In particular, rapid rebindings of the enzyme molecules to the substrate molecules after modification of the first site can markedly speed up the response, and lead to loss of ultrasensitivity and bistability. In essence, rapid enzyme-substrate rebindings can turn a distributive mechanism into a processive mechanism. We argue that slow ADP release by the enzymes can protect the system against these rapid rebindings, thus enabling ultrasensitivity and bistability.

preprint2009arXiv

The role of spatial averaging in the precision of gene expression patterns

During embryonic development, differentiating cells respond via gene expression to positional cues from morphogen gradients. While gene expression is often highly erratic, embryonic development is precise. We show by theory and simulations that diffusion of the expressed protein can enhance the precision of its expression domain. While diffusion lessens the sharpness of the expression boundary, it also reduces super-Poissonian noise by washing out bursts of gene expression. Balancing these effects yields an optimal diffusion constant maximizing the precision of the expression domain.

preprint2008arXiv

Regulatory control and the costs and benefits of biochemical noise

Experiments in recent years have vividly demonstrated that gene expression can be highly stochastic. How protein concentration fluctuations affect the growth rate of a population of cells, is, however, a wide open question. We present a mathematical model that makes it possible to quantify the effect of protein concentration fluctuations on the growth rate of a population of genetically identical cells. The model predicts that the population's growth rate depends on how the growth rate of a single cell varies with protein concentration, the variance of the protein concentration fluctuations, and the correlation time of these fluctuations. The model also predicts that when the average concentration of a protein is close to the value that maximizes the growth rate, fluctuations in its concentration always reduce the growth rate. However, when the average protein concentration deviates sufficiently from the optimal level, fluctuations can enhance the growth rate of the population, even when the growth rate of a cell depends linearly on the protein concentration. The model also shows that the ensemble or population average of a quantity, such as the average protein expression level or its variance, is in general not equal to its time average as obtained from tracing a single cell and its descendants. We apply our model to perform a cost-benefit analysis of gene regulatory control. Our analysis predicts that the optimal expression level of a gene regulatory protein is determined by the trade-off between the cost of synthesizing the regulatory protein and the benefit of minimizing the fluctuations in the expression of its target gene. We discuss possible experiments that could test our predictions.