-
Predicting concentration changes via discrete sampling
Authors:
Age J. Tjalma,
Pieter Rein ten Wolde
Abstract:
To successfully navigate chemical gradients, microorganisms need to predict how the ligand concentration changes in space. Due to their limited size, they do not take a spatial derivative over their body length but rather a temporal derivative, comparing the current signal with that in the recent past, over the so-called adaptation time. This strategy is pervasive in biology, but it remains unclea…
▽ More
To successfully navigate chemical gradients, microorganisms need to predict how the ligand concentration changes in space. Due to their limited size, they do not take a spatial derivative over their body length but rather a temporal derivative, comparing the current signal with that in the recent past, over the so-called adaptation time. This strategy is pervasive in biology, but it remains unclear what determines the accuracy of such measurements. Using a generalized version of the previously established sampling framework, we investigate how resource limitations and the statistics of the input signal set the optimal design of a well-characterized network that measures temporal concentration changes: the Escherichia coli chemotaxis network. Our results show how an optimal adaptation time arises from the trade-off between the sampling error, caused by the stochastic nature of the network, and the dynamical error, caused by uninformative fluctuations in the input. A larger resource availability reduces the sampling error, which allows for a smaller adaptation time, thereby simultaneously decreasing the dynamical error. Similarly, we find that the optimal adaptation time scales inversely with the gradient steepness, because steeper gradients lift the signal above the noise and reduce the sampling error. These findings shed light on the principles that govern the optimal design of the E. coli chemotaxis network specifically, and any system measuring temporal changes more broadly.
△ Less
Submitted 8 February, 2024;
originally announced February 2024.
-
Predicting a noisy signal: the costs and benefits of time averaging as a noise mitigation strategy
Authors:
Jenny Poulton,
Age Tjalma,
Lotte Slim,
Pieter Rein ten Wolde
Abstract:
One major challenge for living cells is the measurement and prediction of signals corrupted by noise. In general, cells need to make decisions based on their compressed representation of noisy, time-varying signals. Strategies for signal noise mitigation are often tackled using Wiener filtering theory, but this theory cannot account for systems that have limited resources and hence must compress t…
▽ More
One major challenge for living cells is the measurement and prediction of signals corrupted by noise. In general, cells need to make decisions based on their compressed representation of noisy, time-varying signals. Strategies for signal noise mitigation are often tackled using Wiener filtering theory, but this theory cannot account for systems that have limited resources and hence must compress the signal. To study how accurately linear systems can predict noisy, time-varying signals in the presence of a compression constraint, we extend the information bottleneck method. We show that the optimal integration kernel reduces to the Wiener filter in the absence of a compression constraint. This kernel combines a delta function at short times and an exponential function that decays on a timescale that sets the integration time. Moreover, there exists an optimal integration time, which arises from a trade-off between time averaging signal noise and dynamical error. As the level of compression is increased, time averaging becomes harder, and as a result the optimal integration time decreases and the delta peak increases. We compare the behaviour of the optimal system with that of a canonical motif in cell signalling, the push-pull network, finding that the system reacts to signal noise and compression in a similar way.
△ Less
Submitted 6 July, 2023;
originally announced July 2023.
-
Stable developmental patterns of gene expression without morphogen gradients
Authors:
Maciej Majka,
Nils B. Becker,
Pieter Rein ten Wolde,
Marcin Zagorski,
Thomas R. Sokolowski
Abstract:
Gene expression patterns (GEPs) are established by cross-regulating target genes that interpret morphogen gradients. However, as development progresses, morphogen activity is reduced, leaving the emergent GEP without stabilizing positional cues. The GEP then can be deteriorated by the intrinsically noisy biochemical processes acting at the cellular level. However, the established GEPs remain spati…
▽ More
Gene expression patterns (GEPs) are established by cross-regulating target genes that interpret morphogen gradients. However, as development progresses, morphogen activity is reduced, leaving the emergent GEP without stabilizing positional cues. The GEP then can be deteriorated by the intrinsically noisy biochemical processes acting at the cellular level. However, the established GEPs remain spatio-temporally stable in many biological systems. Here we combine spatial-stochastic simulations with an enhanced sampling method (Non-Stationary Forward Flux Sampling) and a recently developed stability theory to address how spatiotemporal integrity of a GEP is maintained without morphogen gradients. Using a minimal embryo model consisting of spatially coupled biochemical reactor volumes, we study a stripe pattern in which weak cross-repression between nearest neighbor domians alternates with strong repression between next-nearest neighbor domains, inspired by the gap gene system in the Drosophila embryo. We find that fine-tuning of the weak repressive interactions to an optimal level increases temporal stability of GEPs by orders of magnitude, providing stability over developmentally relevant times, without morphogen gradients. The numerically determined optimal parameter regime closely agrees with the predictions of the stability theory. By analizing the dynamics of GEP asymmetry factors, we trace back the pattern stability enhancement to the emergence of a metastable basin and restoring forces that counteract pattern perturbations. The origin of these forces is further explained by the effective model, describing the emergent deterministic dynamics of the system. Altogether, we show that metastable attractors can emerge as a property of stochastic GEPs even without system-wide positional cues, provided that the gene regulatory interactions shaping the pattern are optimally tuned.
△ Less
Submitted 5 April, 2024; v1 submitted 1 June, 2023;
originally announced June 2023.
-
Trade-offs between cost and information in cellular prediction
Authors:
Age J. Tjalma,
Vahe Galstyan,
Jeroen Goedhart,
Lotte Slim,
Nils B. Becker,
Pieter Rein ten Wolde
Abstract:
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 sig…
▽ More
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.
△ Less
Submitted 10 January, 2023;
originally announced January 2023.
-
Constriction of actin rings by passive crosslinkers
Authors:
Alexander Cumberworth,
Pieter Rein ten Wolde
Abstract:
In many organisms, cell division is driven by the constriction of a cytokinetic ring, which consists of actin filaments and crosslinking proteins. While it has long been believed that the constriction is driven by motor proteins, it has recently been discovered that passive crosslinkers that do not turn over fuel are able to generate enough force to constrict actin filament rings. To study the rin…
▽ More
In many organisms, cell division is driven by the constriction of a cytokinetic ring, which consists of actin filaments and crosslinking proteins. While it has long been believed that the constriction is driven by motor proteins, it has recently been discovered that passive crosslinkers that do not turn over fuel are able to generate enough force to constrict actin filament rings. To study the ring constriction dynamics, we develop a model that includes the driving force of crosslinker condensation and the opposing forces of friction and filament bending. We analyze the constriction force as a function of ring topology and crosslinker concentration, and predict forces that are sufficient to constrict an unadorned plasma membrane. Our model also predicts that actin-filament sliding arises from an interplay between filament rotation and crosslinker hopping, producing frictional forces that are low compared to those of crosslinker-mediated microtubule sliding.
△ Less
Submitted 17 March, 2023; v1 submitted 8 March, 2022;
originally announced March 2022.
-
Path Weight Sampling: Exact Monte Carlo Computation of the Mutual Information between Stochastic Trajectories
Authors:
Manuel Reinhardt,
Gašper Tkačik,
Pieter Rein ten Wolde
Abstract:
Most natural and engineered information-processing systems transmit information via signals that vary in time. Computing the information transmission rate or the information encoded in the temporal characteristics of these signals, requires the mutual information between the input and output signals as a function of time, i.e. between the input and output trajectories. Yet, this is notoriously dif…
▽ More
Most natural and engineered information-processing systems transmit information via signals that vary in time. Computing the information transmission rate or the information encoded in the temporal characteristics of these signals, requires the mutual information between the input and output signals as a function of time, i.e. between the input and output trajectories. Yet, this is notoriously difficult because of the high-dimensional nature of the trajectory space, and all existing techniques require approximations. We present an exact Monte Carlo technique called Path Weight Sampling (PWS) that, for the first time, makes it possible to compute the mutual information between input and output trajectories for any stochastic system that is described by a master equation. The principal idea is to use the master equation to evaluate the exact conditional probability of an individual output trajectory for a given input trajectory, and average this via Monte Carlo sampling in trajectory space to obtain the mutual information. We present three variants of PWS, which all generate the trajectories using the standard stochastic simulation algorithm. While Direct PWS is a brute-force method, Rosenbluth-Rosenbluth PWS exploits the analogy between signal trajectory sampling and polymer sampling, and Thermodynamic Integration PWS is based on a reversible work calculation in trajectory space. PWS also makes it possible to compute the mutual information between input and output trajectories for systems with hidden internal states as well as systems with feedback from output to input. Applying PWS to the bacterial chemotaxis system, consisting of 182 coupled chemical reactions, demonstrates not only that the scheme is highly efficient, but also that the number of receptor clusters is much smaller than hitherto believed, while their size is much larger.
△ Less
Submitted 4 October, 2023; v1 submitted 7 March, 2022;
originally announced March 2022.
-
Energetic constraints on filament mediated cell polarization
Authors:
Harmen Wierenga,
Pieter Rein ten Wolde
Abstract:
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. W…
▽ More
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.
△ Less
Submitted 30 April, 2021;
originally announced April 2021.
-
Diffusible crosslinkers cause superexponential friction forces
Authors:
Harmen Wierenga,
Pieter Rein ten Wolde
Abstract:
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 ho…
▽ More
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.
△ Less
Submitted 6 June, 2019;
originally announced June 2019.
-
The behaviour of information flow near criticality
Authors:
Matthijs Meijers,
Sosuke Ito,
Pieter Rein ten Wolde
Abstract:
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,…
▽ More
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.
△ Less
Submitted 3 June, 2019;
originally announced June 2019.
-
Theory for the optimal detection of time-varying signals in cellular sensing systems
Authors:
G. Malaguti,
P. R. ten Wolde
Abstract:
Living cells often need to measure chemical concentrations that vary in time. To this end, they deploy many resources, e.g. receptors, downstream signaling molecules, time and energy. Here, we present a theory for the optimal design of a large class of sensing systems that need to detect time-varying signals, a receptor driving a push-pull network. The theory is based on the concept of the dynamic…
▽ More
Living cells often need to measure chemical concentrations that vary in time. To this end, they deploy many resources, e.g. receptors, downstream signaling molecules, time and energy. Here, we present a theory for the optimal design of a large class of sensing systems that need to detect time-varying signals, a receptor driving a push-pull network. The theory is based on the concept of the dynamic input-output relation, which describes the mapping between the current ligand concentration and the average receptor occupancy over the past integration time. This concept is used to develop the idea that the cell employs its push-pull network to estimate the receptor occupancy and then uses this estimate to infer the current ligand concentration by inverting the dynamic input-output relation. The theory reveals that the sensing error can be decomposed into two terms: the sampling error in the estimate of the receptor occupancy and the dynamical error that arises because the average ligand concentration over the past integration time may not reflect the current ligand concentration. The theory generalizes the design principle of optimal resource allocation previously identified for static signals, which states that in an optimally designed sensing system the three fundamental resource classes of receptors and their integration time, readout molecules, and energy are equally limiting so that no resource is wasted. However, in contrast to static signals, receptors and power cannot be traded freely against time to reach a desired sensing precision: there exists an optimal integration time that maximizes the sensing precision, which depends on the number of receptors, the receptor correlation time, and the correlation time and variance of the input signal. Applying our theory to the chemotaxis system of Escherichia coli indicates that this bacterium has evolved to optimally sense shallow gradients.
△ Less
Submitted 25 February, 2019;
originally announced February 2019.
-
Theory of circadian metabolism
Authors:
Michele Monti,
David K. Lubensky,
Pieter Rein ten Wolde
Abstract:
Many organisms repartition their proteome in a circadian fashion in response to the daily nutrient changes in their environment. A striking example is provided by cyanobacteria, which perform photosynthesis during the day to fix carbon. These organisms not only face the challenge of rewiring their proteome every 12 hours, but also the necessity of storing the fixed carbon in the form of glycogen t…
▽ More
Many organisms repartition their proteome in a circadian fashion in response to the daily nutrient changes in their environment. A striking example is provided by cyanobacteria, which perform photosynthesis during the day to fix carbon. These organisms not only face the challenge of rewiring their proteome every 12 hours, but also the necessity of storing the fixed carbon in the form of glycogen to fuel processes during the night. In this manuscript, we extend the framework developed by Hwa and coworkers (Scott et al., Science 330, 1099 (2010)) for quantifying the relatinship between growth and proteome composition to circadian metabolism. We then apply this framework to investigate the circadian metabolism of the cyanobacterium Cyanothece, which not only fixes carbon during the day, but also nitrogen during the night, storing it in the polymer cyanophycin. Our analysis reveals that the need to store carbon and nitrogen tends to generate an extreme growth strategy, in which the cells predominantly grow during the day, as observed experimentally. This strategy maximizes the growth rate over 24 hours, and can be quantitatively understood by the bacterial growth laws. Our analysis also shows that the slow relaxation of the proteome, arising from the slow growth rate, puts a severe constraint on implementing this optimal strategy. Yet, the capacity to estimate the time of the day, enabled by the circadian clock, makes it possible to anticipate the daily changes in the environment and mount a response ahead of time. This significantly enhances the growth rate by counteracting the detrimental effects of the slow proteome relaxation.
△ Less
Submitted 11 May, 2018;
originally announced May 2018.
-
Quantifying fluctuations in reversible enzymatic cycles and clocks
Authors:
Harmen Wierenga,
Pieter Rein ten Wolde,
Nils B. Becker
Abstract:
Biochemical reactions are fundamentally noisy at a molecular scale. This limits the precision of reaction networks, but also allows fluctuation measurements which may reveal the structure and dynamics of the underlying biochemical network. Here, we study non-equilibrium reaction cycles, such as the mechanochemical cycle of molecular motors, the phosphorylation cycle of circadian clock proteins, or…
▽ More
Biochemical reactions are fundamentally noisy at a molecular scale. This limits the precision of reaction networks, but also allows fluctuation measurements which may reveal the structure and dynamics of the underlying biochemical network. Here, we study non-equilibrium reaction cycles, such as the mechanochemical cycle of molecular motors, the phosphorylation cycle of circadian clock proteins, or the transition state cycle of enzymes. Fluctuations in such cycles may be measured using either of two classical definitions of the randomness parameter, which we show to be equivalent in general microscopically reversible cycles. We define a stochastic period for reversible cycles and present analytical solutions for its moments. Furthermore, we associate the two forms of the randomness parameter with the thermodynamic uncertainty relation, which sets limits on the timing precision of the cycle in terms of thermodynamic quantities. Our results should prove useful also for the study of temporal fluctuations in more general networks.
△ Less
Submitted 9 March, 2018; v1 submitted 16 January, 2018;
originally announced January 2018.
-
Robustness of clocks to input noise
Authors:
Michele Monti,
David K Lubensky,
Pieter Rein ten Wolde
Abstract:
To estimate the time, many organisms, ranging from cyanobacteria to animals, employ a circadian clock which is based on a limit-cycle oscillator that can tick autonomously with a nearly 24h period. Yet, a limit-cycle oscillator is not essential for knowing the time, as exemplified by bacteria that possess an 'hourglass': a system that when forced by an oscillatory light input exhibits robust oscil…
▽ More
To estimate the time, many organisms, ranging from cyanobacteria to animals, employ a circadian clock which is based on a limit-cycle oscillator that can tick autonomously with a nearly 24h period. Yet, a limit-cycle oscillator is not essential for knowing the time, as exemplified by bacteria that possess an 'hourglass': a system that when forced by an oscillatory light input exhibits robust oscillations from which the organism can infer the time, but that in the absence of driving relaxes to a stable fixed point. Here, using models of the Kai system of cyanobacteria, we compare a limit- cycle oscillator with two hourglass models, one that without driving relaxes exponentially and one that does so in an oscillatory fashion. In the limit of low input-noise, all three systems are equally informative on time, yet in the regime of high input-noise the limit-cycle oscillator is far superior. The same behavior is found in the Stuart-Landau model, indicating that our result is universal.
△ Less
Submitted 11 May, 2018; v1 submitted 27 September, 2017;
originally announced October 2017.
-
eGFRD in all dimensions
Authors:
Thomas R. Sokolowski,
Joris Paijmans,
Laurens Bossen,
Martijn Wehrens,
Thomas Miedema,
Nils B. Becker,
Kazunari Kaizu,
Koichi Takahashi,
Marlieen Dogterom,
Pieter Rein ten Wolde
Abstract:
Biochemical reactions typically occur at low copy numbers, but at once in crowded and diverse environments. Space and stochasticity therefore play an essential role in biochemical networks. Spatial-stochastic simulations have become a prominent tool for understanding how stochasticity at the microscopic level influences the macroscopic behavior of such systems. However, while particle-based models…
▽ More
Biochemical reactions typically occur at low copy numbers, but at once in crowded and diverse environments. Space and stochasticity therefore play an essential role in biochemical networks. Spatial-stochastic simulations have become a prominent tool for understanding how stochasticity at the microscopic level influences the macroscopic behavior of such systems. However, while particle-based models guarantee the level of detail necessary to accurately describe the microscopic dynamics at very low copy numbers, the algorithms used to simulate them oftentimes imply trade-offs between computational efficiency and accuracy. eGFRD (enhanced Green's Function Reaction Dynamics) is an exact algorithm that evades such trade-offs by partitioning the N-particle system into M<N analytically tractable one- and two-particle systems; the analytical solutions (Green's functions) then are used to implement an event-driven particle-based scheme that allows particles to make large jumps in time and space while retaining access to their state variables at any moment. Here we present "eGFRD2", a new eGFRD version that implements the principle of eGFRD in all dimensions, enabling efficient simulation of biochemical reaction-diffusion processes in the 3D cytoplasm, on 2D planes representing membranes, and on 1D elongated cylinders representative of, e.g., cytoskeletal tracks or DNA; in 1D, it also incorporates convective motion used to model active transport. We find that, for low particle densities, eGFRD2 is up to 3 orders of magnitude faster than optimized Brownian Dynamics. We exemplify the capabilities of eGFRD2 by simulating an idealized model of Pom1 gradient formation, which involves 3D diffusion, active transport on microtubules, and autophosphorylation on the membrane, confirming recent results on this system and demonstrating that it can efficiently operate under genuinely stochastic conditions.
△ Less
Submitted 30 August, 2017;
originally announced August 2017.
-
Optimal entrainment of circadian clocks in the presence of noise
Authors:
Michele Monti,
David K. Lubensky,
Pieter Rein ten Wolde
Abstract:
Circadian clocks are biochemical oscillators that allow organisms to estimate the time of the day. These oscillators are inherently noisy due to the discrete nature of the reactants and the stochastic character of their interactions. To keep these oscillators in sync with the daily day-night rhythm in the presence of noise, circadian clocks must be coupled to the dark-light cycle. In this manuscri…
▽ More
Circadian clocks are biochemical oscillators that allow organisms to estimate the time of the day. These oscillators are inherently noisy due to the discrete nature of the reactants and the stochastic character of their interactions. To keep these oscillators in sync with the daily day-night rhythm in the presence of noise, circadian clocks must be coupled to the dark-light cycle. In this manuscript, we study the entrainment of phase oscillators as a function of the intrinsic noise in the system. Using stochastic simulations, we compute the optimal coupling strength, intrinsic frequency and shape of the phase-response curve, that maximize the mutual information between the phase of the clock and time. We show that the optimal coupling strength and intrinsic frequency increase with the noise, but that the shape of the phase-response curve varies non-monotonically with the noise: in the low-noise regime, it features a deadzone that increases in width as the noise increases, while in the high-noise regime, the width decreases with the noise. These results arise from a trade-off between maximizing stability -noise suppression - and maximizing linearity of the input-output, i.e. time- phase, relation. We also show that three analytic approximations -the linear-noise approximation, the phase-averaging method, and linear-response theory- accurately describe different regimes of the coupling strength and the noise.
△ Less
Submitted 7 June, 2017;
originally announced June 2017.
-
Spatial-Stochastic Simulation of Reaction-Diffusion Systems
Authors:
Thomas R. Sokolowski,
Pieter Rein ten Wolde
Abstract:
Biochemical networks play a crucial role in biological systems, implementing a broad range of vital functions. They normally operate at low copy numbers and in spatial settings, but this is often ignored and well-stirred conditions are assumed. Yet, it is increasingly becoming clear that even microscopic spatial inhomogeneities oftentimes can induce significant differences on the macroscopic level…
▽ More
Biochemical networks play a crucial role in biological systems, implementing a broad range of vital functions. They normally operate at low copy numbers and in spatial settings, but this is often ignored and well-stirred conditions are assumed. Yet, it is increasingly becoming clear that even microscopic spatial inhomogeneities oftentimes can induce significant differences on the macroscopic level. Since experimental observation of single-molecule behavior is extremely challenging, theoretical modeling of biochemical reactions on the single-particle level is an important tool for understanding spatial effects in biochemical systems. While purely analytical models quickly become intractable here, spatial-stochastic simulations can capture a wide range of biochemical processes with the necessary levels of detail. Here we discuss different techniques for spatial-stochastic simulation of reaction-diffusion systems, and explain important precautions required to make them biochemically accurate and efficient. We illustrate non-negligible accuracy issues arising even in the most simple approaches to biochemical simulation, and present methods to deal with them. We first explain how Brownian Dynamics, a widely used particle-based diffusion simulation technique with fixed propagation time, can be adapted to incorporate chemical reactions, and portray a range of schemes that elaborate on this idea. We then introduce event-driven spatial-stochastic simulation methods, in which system updates are performed asynchronously with situation-dependent, varying time steps; here we particularly focus on eGFRD, a computationally efficient particle-based algorithm that makes use of analytical functions to accurately sample interparticle reactions and diffusive motion with large jumps in time and space. We end by briefly presenting recent developments in the field of spatial-stochastic biochemical simulation.
△ Less
Submitted 24 May, 2017;
originally announced May 2017.
-
Period Robustness and Entrainability of the Kai System to Changing Nucleotide Concentrations
Authors:
Joris Paijmans,
David K Lubensky,
Pieter Rein ten Wolde
Abstract:
Circadian clocks must be able to entrain to time-varying signals to keep their oscillations in phase with the day-night rhythm. On the other hand, they must also exhibit input compensation: their period must remain about one day in different constant environments. The post-translational oscillator of the Kai system can be entrained by transient or oscillatory changes in the ATP fraction, yet is in…
▽ More
Circadian clocks must be able to entrain to time-varying signals to keep their oscillations in phase with the day-night rhythm. On the other hand, they must also exhibit input compensation: their period must remain about one day in different constant environments. The post-translational oscillator of the Kai system can be entrained by transient or oscillatory changes in the ATP fraction, yet is insensitive to constant changes in this fraction. We study in three different models of this system how these two seemingly conflicting criteria are met: the Van Zon model (Van Zon et al., PNAS, 2007), the Rust model (Phong et al., PNAS, 2013), and our new model (Paijmans et al., arXiv:1612.02715). We find that the new model exhibits the best trade-off between input compensation and entrainability: on the footing of equal phase-response curves, it exhibits the strongest input compensation. Performing stochastic simulations at the level of individual hexamers allows us to identify a new mechanism, which is employed by the new model to achieve input compensation: At lower ATP fraction, the individual hexamers make a shorter cycle in the phosphorylation state space, which compensates for the slower pace at which they traverse the cycle.
△ Less
Submitted 22 June, 2017; v1 submitted 25 December, 2016;
originally announced December 2016.
-
Robustness of synthetic oscillators in growing and dividing cells
Authors:
Joris Paijmans,
David K Lubensky,
Pieter Rein ten Wolde
Abstract:
Synthetic biology sets out to implement new functions in cells, and to develop a deeper understanding of biological design principles. In 2000, Elowitz and Leibler showed that by rational design of the reaction network, and using existing biological components, they could create a network that exhibits periodic gene expression, dubbed the repressilator (Elowitz and Leibler, Nature, 2000). More rec…
▽ More
Synthetic biology sets out to implement new functions in cells, and to develop a deeper understanding of biological design principles. In 2000, Elowitz and Leibler showed that by rational design of the reaction network, and using existing biological components, they could create a network that exhibits periodic gene expression, dubbed the repressilator (Elowitz and Leibler, Nature, 2000). More recently, Stricker et al. presented another synthetic oscillator, called the dual-feedback oscillator (Stricker et al., 2008), which is more stable. How the stability of these oscillators is affected by the intrinsic noise of the interactions between the components and the stochastic expression of their genes, has been studied in considerable detail. However, as all biological oscillators reside in growing and dividing cells, an important question is how these oscillators are perturbed by the cell cycle. In previous work we showed that the periodic doubling of the gene copy numbers due to DNA replication can couple not only natural, circadian oscillators to the cell cycle (Paijmans et al., PNAS, \textbf{113}, 4063, (2016)), but also these synthetic oscillators. Here we expand this study. We find that the strength of the locking between oscillators depends not only on the positions of the genes on the chromosome, but also on the noise in the timing of gene replication: noise tends to weaken the coupling. Yet, even in the limit of high levels of noise in the replication times of the genes, both synthetic oscillators show clear signatures of locking to the cell cycle. This work enhances our understanding of the design of robust biological oscillators inside growing and diving cells.
△ Less
Submitted 23 December, 2016;
originally announced December 2016.
-
A Thermodynamically consistent model of the post-translational Kai circadian clock
Authors:
Joris Paijmans,
David K Lubensky,
Pieter Rein ten Wolde
Abstract:
The principal pacemaker of the circadian clock of the cyanobacterium S. elongatus is a protein phosphorylation cycle consisting of three proteins, KaiA, KaiB and KaiC. KaiC forms a homohexamer, with each monomer consisting of two domains, CI and CII. Both domains can bind and hydrolyze ATP, but only the CII domain can be phosphorylated, at two residues, in a well-defined sequence. While this syste…
▽ More
The principal pacemaker of the circadian clock of the cyanobacterium S. elongatus is a protein phosphorylation cycle consisting of three proteins, KaiA, KaiB and KaiC. KaiC forms a homohexamer, with each monomer consisting of two domains, CI and CII. Both domains can bind and hydrolyze ATP, but only the CII domain can be phosphorylated, at two residues, in a well-defined sequence. While this system has been studied extensively, how the clock is driven thermodynamically has remained elusive. Inspired by recent experimental observations and building on ideas from previous mathematical models, we present a new, thermodynamically consistent, statistical-mechanical model of the clock. At its heart are two main ideas: i) ATP hydrolysis in the CI domain provides the thermodynamic driving force for the clock, switching KaiC between an active conformational state in which its phosphorylation level tends to rise and an inactive one in which it tends to fall; ii) phosphorylation of the CII domain provides the timer for the hydrolysis in the CI domain. The model also naturally explains how KaiA, by acting as a nucleotide exchange factor, can stimulate phosphorylation of KaiC, and how the differential affinity of KaiA for the different KaiC phosphoforms generates the characteristic temporal order of KaiC phosphorylation. As the phosphorylation level in the CII domain rises, the release of ADP from CI slows down, making the inactive conformational state of KaiC more stable. In the inactive state, KaiC binds KaiB, which not only stabilizes this state further, but also leads to the sequestration of KaiA, and hence to KaiC dephosphorylation. Using a dedicated kinetic Monte Carlo algorithm, we show that the model can describe a wealth of experimental data.
△ Less
Submitted 8 December, 2016;
originally announced December 2016.
-
Multiscale simulations of anisotropic particles combining Brownian Dynamics and Green's Function Reaction Dynamics
Authors:
Adithya Vijaykumar,
Thomas E. Ouldridge,
Pieter Rein ten Wolde,
Peter G. Bolhuis
Abstract:
The modeling of complex reaction-diffusion processes in, for instance, cellular biochemical networks or self-assembling soft matter can be tremendously sped up by employing a multiscale algorithm which combines the mesoscopic Green's Function Reaction Dynamics (GFRD) method with explicit stochastic Brownian, Langevin, or deterministic Molecular Dynamics to treat reactants at the microscopic scale…
▽ More
The modeling of complex reaction-diffusion processes in, for instance, cellular biochemical networks or self-assembling soft matter can be tremendously sped up by employing a multiscale algorithm which combines the mesoscopic Green's Function Reaction Dynamics (GFRD) method with explicit stochastic Brownian, Langevin, or deterministic Molecular Dynamics to treat reactants at the microscopic scale [A. Vijaykumar, P.G. Bolhuis and P.R. ten Wolde, J. Chem. Phys. {\bf 43}, 21: 214102 (2015)]. Here we extend this multiscale BD-GFRD approach to include the orientational dynamics that is crucial to describe the anisotropic interactions often prevalent in biomolecular systems. We illustrate the novel algorithm using a simple patchy particle model. After validation of the algorithm we discuss its performance. The rotational BD-GFRD multiscale method will open up the possibility for large scale simulations of e.g. protein signalling networks.
△ Less
Submitted 28 November, 2016;
originally announced November 2016.
-
Fundamental costs in the production and destruction of persistent polymer copies
Authors:
Thomas E. Ouldridge,
Pieter Rein ten Wolde
Abstract:
Living cells use readout molecules to record the state of receptor proteins, similar to measurements or copies in typical computational devices. But is this analogy rigorous? Can cells be optimally efficient, and if not, why? We show that, as in computation, a canonical biochemical readout network generates correlations; extracting no work from these correlations sets a lower bound on dissipation.…
▽ More
Living cells use readout molecules to record the state of receptor proteins, similar to measurements or copies in typical computational devices. But is this analogy rigorous? Can cells be optimally efficient, and if not, why? We show that, as in computation, a canonical biochemical readout network generates correlations; extracting no work from these correlations sets a lower bound on dissipation. For general input, the biochemical network cannot reach this bound, even with arbitrarily slow reactions or weak thermodynamic driving. It faces an accuracy-dissipation trade-off that is qualitatively distinct from and worse than implied by the bound, and more complex steady-state copy processes cannot perform better. Nonetheless, the cost remains close to the thermodynamic bound unless accuracy is extremely high. Additionally, we show that biomolecular reactions could be used in thermodynamically optimal devices under exogenous manipulation of chemical fuels, suggesting an experimental system for testing computational thermodynamics.
△ Less
Submitted 11 March, 2017; v1 submitted 18 September, 2016;
originally announced September 2016.
-
Biochemical machines for the interconversion of mutual information and work
Authors:
Thomas McGrath,
Nick S. Jones,
Pieter Rein ten Wolde,
Thomas E. Ouldridge
Abstract:
We propose a physically-realisable biochemical device that is coupled to a biochemical reservoir of mutual information, fuel molecules and a chemical bath. Mutual information allows work to be done on the bath even when the fuel molecules appear to be in equilibrium; alternatively, mutual information can be created by driving from the fuel or the bath. The system exhibits diverse behaviour, includ…
▽ More
We propose a physically-realisable biochemical device that is coupled to a biochemical reservoir of mutual information, fuel molecules and a chemical bath. Mutual information allows work to be done on the bath even when the fuel molecules appear to be in equilibrium; alternatively, mutual information can be created by driving from the fuel or the bath. The system exhibits diverse behaviour, including a regime in which the information, despite increasing during the reaction, enhances the extracted work. We further demonstrate that a modified device can function without the need for external manipulation, eliminating the need for a complex and potentially costly control.
△ Less
Submitted 26 November, 2016; v1 submitted 19 April, 2016;
originally announced April 2016.
-
The accuracy of telling time via oscillatory signals
Authors:
Michele Monti,
Pieter Rein ten Wolde
Abstract:
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…
▽ More
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.
△ Less
Submitted 14 January, 2016;
originally announced January 2016.
-
Optimal Prediction by Cellular Signaling Networks
Authors:
Nils B. Becker,
Andrew Mugler,
Pieter Rein ten Wolde
Abstract:
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, whic…
▽ More
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.
△ Less
Submitted 7 December, 2015;
originally announced December 2015.
-
Molecular clustering digitizes signaling and increases fidelity
Authors:
Edward Roob III,
Nicola Trendel,
Pieter Rein ten Wolde,
Andrew Mugler
Abstract:
Many membrane-bound molecules in cells form small clusters. It has been hypothesized that these clusters convert an analog extracellular signal into a digital intracellular signal and that this conversion increases signaling fidelity. However, the mechanism by which clusters digitize a signal and the subsequent effects on fidelity remain poorly understood. Here we demonstrate using a stochastic mo…
▽ More
Many membrane-bound molecules in cells form small clusters. It has been hypothesized that these clusters convert an analog extracellular signal into a digital intracellular signal and that this conversion increases signaling fidelity. However, the mechanism by which clusters digitize a signal and the subsequent effects on fidelity remain poorly understood. Here we demonstrate using a stochastic model of cooperative cluster formation that sufficient cooperation leads to digital signaling. We show that despite reducing the number of output states, which decreases fidelity, digitization also reduces noise in the system, which increases fidelity. The tradeoff between these effects leads to an optimal cluster size that agrees with experimental measurements.
△ Less
Submitted 2 November, 2015;
originally announced November 2015.
-
Discrete gene replication events drive coupling between the cell cycle and circadian clocks
Authors:
Joris Paijmans,
Mark Bosman,
Pieter Rein ten Wolde,
David K. Lubensky
Abstract:
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 regula…
▽ More
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.
△ Less
Submitted 20 October, 2015;
originally announced October 2015.
-
Fundamental Limits to Cellular Sensing
Authors:
Pieter Rein ten Wolde,
Nils B. Becker,
Thomas E. Ouldridge,
A. Mugler
Abstract:
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…
▽ More
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.
△ Less
Submitted 25 May, 2015;
originally announced May 2015.
-
The thermodynamics of computational copying in biochemical systems
Authors:
Thomas E. Ouldridge,
Christopher C. Govern,
Pieter Rein ten Wolde
Abstract:
Living cells use readout molecules to record the state of receptor proteins, similar to measurements or copies in typical computational devices. But is this analogy rigorous? Can cells be optimally efficient, and if not, why? We show that, as in computation, a canonical biochemical readout network generates correlations; extracting no work from these correlations sets a lower bound on dissipation.…
▽ More
Living cells use readout molecules to record the state of receptor proteins, similar to measurements or copies in typical computational devices. But is this analogy rigorous? Can cells be optimally efficient, and if not, why? We show that, as in computation, a canonical biochemical readout network generates correlations; extracting no work from these correlations sets a lower bound on dissipation. For general input, the biochemical network cannot reach this bound, even with arbitrarily slow reactions or weak thermodynamic driving. It faces an accuracy-dissipation trade-off that is qualitatively distinct from and worse than implied by the bound, and more complex steady-state copy processes cannot perform better. Nonetheless, the cost remains close to the thermodynamic bound unless accuracy is extremely high. Additionally, we show that biomolecular reactions could be used in thermodynamically optimal devices under exogenous manipulation of chemical fuels, suggesting an experimental system for testing computational thermodynamics.
△ Less
Submitted 9 March, 2017; v1 submitted 3 March, 2015;
originally announced March 2015.
-
Positive feedback can lead to dynamic nanometer-scale clustering on cell membranes
Authors:
Martijn Wehrens,
Pieter Rein ten Wolde,
Andrew Mugler
Abstract:
Clustering of molecules on biological membranes is a widely observed phenomenon. In some cases, such as the clustering of Ras proteins on the membranes of mammalian cells, proper cell signaling is critically dependent on the maintenance of these clusters. Yet, the mechanism by which clusters form and are maintained in these systems remains unclear. Recently, it has been discovered that activated R…
▽ More
Clustering of molecules on biological membranes is a widely observed phenomenon. In some cases, such as the clustering of Ras proteins on the membranes of mammalian cells, proper cell signaling is critically dependent on the maintenance of these clusters. Yet, the mechanism by which clusters form and are maintained in these systems remains unclear. Recently, it has been discovered that activated Ras promotes further Ras activation. Here we show using particle-based simulation that this positive feedback is sufficient to produce persistent clusters of active Ras molecules at the nanometer scale via a dynamic nucleation mechanism. Furthermore, we find that our cluster statistics are consistent with experimental observations of the Ras system. Interestingly, we show that our model does not support a Turing regime of macroscopic reaction-diffusion patterning, and therefore that the clustering we observe is a purely stochastic effect, arising from the coupling of positive feedback with the discrete nature of individual molecules. These results underscore the importance of stochastic and dynamic properties of reaction diffusion systems for biological behavior.
△ Less
Submitted 26 August, 2014;
originally announced August 2014.
-
The robustness of proofreading to crowding-induced pseudo-processivity in the MAPK pathway
Authors:
Thomas E. Ouldridge,
Pieter Rein ten Wolde
Abstract:
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…
▽ More
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.
△ Less
Submitted 7 November, 2014; v1 submitted 6 July, 2014;
originally announced July 2014.
-
The lower bound on the precision of transcriptional regulation
Authors:
Joris Paijmans,
Pieter Rein ten Wolde
Abstract:
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 o…
▽ More
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.
△ Less
Submitted 3 July, 2014;
originally announced July 2014.
-
Prediction and Dissipation in Biochemical Sensing
Authors:
Nils B. Becker,
Andrew Mugler,
Pieter Rein ten Wolde
Abstract:
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,…
▽ More
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.
△ Less
Submitted 19 December, 2013;
originally announced December 2013.
-
The importance of crowding in signaling, genetic, and metabolic networks
Authors:
Pieter Rein ten Wolde,
Andrew Mugler
Abstract:
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…
▽ More
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.
△ Less
Submitted 18 November, 2013;
originally announced November 2013.
-
Multiplexing oscillatory biochemical signals
Authors:
Wiet de Ronde,
Pieter Rein ten Wolde
Abstract:
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…
▽ More
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.
△ Less
Submitted 15 September, 2013;
originally announced September 2013.
-
How biochemical resources determine fundamental limits in cellular sensing
Authors:
Christopher C. Govern,
Pieter Rein ten Wolde
Abstract:
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…
▽ More
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.
△ Less
Submitted 30 January, 2014; v1 submitted 6 August, 2013;
originally announced August 2013.
-
Spatial partitioning improves the reliability of biochemical signaling
Authors:
Andrew Mugler,
Filipe Tostevin,
Pieter Rein ten Wolde
Abstract:
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 bioc…
▽ More
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.
△ Less
Submitted 26 February, 2013;
originally announced February 2013.
-
Fundamental Limits on Sensing Chemical Concentrations with Linear Biochemical Networks
Authors:
Christopher C. Govern,
Pieter Rein ten Wolde
Abstract:
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 nonstatio…
▽ More
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.
△ Less
Submitted 18 October, 2012;
originally announced October 2012.
-
The macroscopic effects of microscopic heterogeneity
Authors:
Andrew Mugler,
Pieter Rein ten Wolde
Abstract:
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 o…
▽ More
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.
△ Less
Submitted 3 October, 2012; v1 submitted 29 September, 2012;
originally announced October 2012.
-
Protein logic: a statistical mechanical study of signal integration at the single-molecule level
Authors:
Wiet de Ronde,
Pieter Rein ten Wolde,
Andrew Mugler
Abstract:
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…
▽ More
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.
△ Less
Submitted 10 September, 2012;
originally announced September 2012.
-
Signalling noise enhances chemotactic drift of E. coli
Authors:
Marlo Flores,
Thomas S. Shimuzu,
Pieter Rein ten Wolde,
Filipe Tostevin
Abstract:
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 opt…
▽ More
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.
△ Less
Submitted 4 October, 2012; v1 submitted 27 August, 2012;
originally announced August 2012.
-
Mutual Repression enhances the Steepness and Precision of Gene Expression Boundaries
Authors:
Thomas R. Sokolowski,
Thorsten Erdmann,
Pieter Rein ten Wolde
Abstract:
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 i…
▽ More
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.
△ Less
Submitted 13 August, 2012; v1 submitted 1 August, 2012;
originally announced August 2012.
-
Reliability of frequency- and amplitude-decoding in gene regulation
Authors:
Filipe Tostevin,
Wiet de Ronde,
Pieter Rein ten Wolde
Abstract:
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 i…
▽ More
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.
△ Less
Submitted 7 February, 2012;
originally announced February 2012.
-
Rare switching events in non-stationary systems
Authors:
Nils B. Becker,
Pieter Rein ten Wolde
Abstract:
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 id…
▽ More
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.
△ Less
Submitted 18 January, 2012;
originally announced January 2012.
-
Non-Stationary Forward Flux Sampling
Authors:
Nils B. Becker,
Rosalind J. Allen,
Pieter Rein ten Wolde
Abstract:
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 probabilit…
▽ More
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.
△ Less
Submitted 18 January, 2012;
originally announced January 2012.
-
Membrane clustering and the role of rebinding in biochemical signaling
Authors:
Andrew Mugler,
Aimee Gotway Bailey,
Koichi Takahashi,
Pieter Rein ten Wolde
Abstract:
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…
▽ More
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.
△ Less
Submitted 25 October, 2011;
originally announced October 2011.
-
Multiplexing Biochemical Signals
Authors:
Wiet de Ronde,
Filipe Tostevin,
Pieter Rein ten Wolde
Abstract:
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 informat…
▽ More
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.
△ Less
Submitted 3 August, 2010;
originally announced August 2010.
-
Robust circadian clocks from coupled protein modification and transcription-translation cycles
Authors:
David Zwicker,
David K. Lubensky,
Pieter Rein ten Wolde
Abstract:
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 prote…
▽ More
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.
△ Less
Submitted 16 April, 2010;
originally announced April 2010.
-
Mutual information in time-varying biochemical systems
Authors:
Filipe Tostevin,
Pieter Rein ten Wolde
Abstract:
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 s…
▽ More
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.
△ Less
Submitted 16 June, 2010; v1 submitted 23 February, 2010;
originally announced February 2010.
-
The effect of feedback on the fidelity of information transmission of time-varying signals
Authors:
Wiet de Ronde,
Filipe Tostevin,
Pieter Rein ten Wolde
Abstract:
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 charact…
▽ More
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.
△ Less
Submitted 12 February, 2010;
originally announced February 2010.
-
The role of spatial averaging in the precision of gene expression patterns
Authors:
Thorsten Erdmann,
Martin Howard,
Pieter Rein ten Wolde
Abstract:
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,…
▽ More
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.
△ Less
Submitted 24 July, 2009;
originally announced July 2009.