Bistability: requirements on cell-volume, protein diffusion, and thermodynamics

Compartmental exchange regulates steady states and stochastic switching of a phosphorylation network

The phosphoregulation of proteins with multiple phosphorylation sites is governed by biochemical reaction networks that can exhibit multistable behavior. However, the behavior of such networks is typically studied in a single reaction volume, while cells are spatially organized into compartments that can exchange proteins. In this work, we use stochastic simulations to study the impact of compartmentalization on a two-site phosphorylation network. We characterize steady states and fluctuation-driven transitions between them as a function of the rate of protein exchange between two compartments. Surprisingly, the average time spent in a state before stochastically switching to another depends nonmonotonically on the protein exchange rate, with the most frequent switching occurring at intermediate exchange rates. At sufficiently small exchange rates, the state of the system and mean switching time are controlled largely by fluctuations in the balance of enzymes in each compartment. This leads to negatively correlated states in the compartments. For large exchange rates, the two compartments behave as a single effective compartment. However, when the compartmental volumes are unequal, the behavior differs from a single compartment with the same total volume. These results demonstrate that exchange of proteins between distinct compartments can regulate the emergent behavior of a common signaling motif.

Decorating the surface of Escherichia coli with bacterial lipoproteins: a comparative analysis of different display systems

Background

The display of recombinant proteins on cell surfaces has a plethora of applications including vaccine development, screening of peptide libraries, whole-cell biocatalysts and biosensor development for diagnostic, industrial or environmental purposes. In the last decades, a wide variety of surface display systems have been developed for the exposure of recombinant proteins on the surface of Escherichia coli, such as autotransporters and outer membrane proteins.

Results

In this study, we assess three approaches for the surface display of a panel of heterologous and homologous mature lipoproteins in E. coli: four from Neisseria meningitidis and four from the host strain that are known to be localised in the inner leaflet of the outer membrane. Constructs were made carrying the sequences coding for eight mature lipoproteins, each fused to the delivery portion of three different systems: the autotransporter adhesin involved in diffuse adherence-I (AIDA-I) from enteropathogenic E. coli, the Lpp’OmpA chimaera and a truncated form of the ice nucleation protein (INP), InaK-NC (N-terminal domain fused with C-terminal one) from Pseudomonas syringae. In contrast to what was observed for the INP constructs, when fused to the AIDA-I or Lpp’OmpA, most of the mature lipoproteins were displayed on the bacterial surface both at 37 and 25 °C as demonstrated by FACS analysis, confocal and transmission electron microscopy.

Conclusions

To our knowledge this is the first study that compares surface display systems using a number of passenger proteins. We have shown that the experimental conditions, including the choice of the carrier protein and the growth temperature, play an important role in the translocation of mature lipoproteins onto the bacterial surface. Despite all the optimization steps performed with the InaK-NC anchor motif, surface exposure of the passenger proteins used in this study was not achieved. For our experimental conditions, Lpp’OmpA chimaera has proved to be an efficient surface display system for the homologous passenger proteins although cell lysis and phenotype heterogeneity were observed. Finally, AIDA-I was found to be the best surface display system for mature lipoproteins (especially heterologous ones) in the E. coli host strain with no inhibition of growth and only limited phenotype heterogeneity.

Quantifying the roles of space and stochasticity in computer simulations for cell biology and cellular biochemistry

Most of the fascinating phenomena studied in cell biology emerge from interactions among highly organized multimolecular structures embedded into complex and frequently dynamic cellular morphologies. For the exploration of such systems, computer simulation has proved to be an invaluable tool, and many researchers in this field have developed sophisticated computational models for application to specific cell biological questions. However, it is often difficult to reconcile conflicting computational results that use different approaches to describe the same phenomenon. To address this issue systematically, we have defined a series of computational test cases ranging from very simple to moderately complex, varying key features of dimensionality, reaction type, reaction speed, crowding, and cell size. We then quantified how explicit spatial and/or stochastic implementations alter outcomes, even when all methods use the same reaction network, rates, and concentrations. For simple cases, we generally find minor differences in solutions of the same problem. However, we observe increasing discordance as the effects of localization, dimensionality reduction, and irreversible enzymatic reactions are combined. We discuss the strengths and limitations of commonly used computational approaches for exploring cell biological questions and provide a framework for decision making by researchers developing new models. As computational power and speed continue to increase at a remarkable rate, the dream of a fully comprehensive computational model of a living cell may be drawing closer to reality, but our analysis demonstrates that it will be crucial to evaluate the accuracy of such models critically and systematically.

Thermodynamics of switching in multistable non-equilibrium systems

Multistable non-equilibrium systems are abundant outcomes of nonlinear dynamics with feedback, but still relatively little is known about what determines the stability of the steady states and their switching rates in terms of entropy and entropy production. Here, we will link fluctuation theorems for the entropy production along trajectories with the action obtainable from the Freidlin-Wentzell theorem to elucidate the thermodynamics of switching between states in the large volume limit of multistable systems. We find that the entropy production at steady state plays no role, but the entropy production during switching is key. Steady-state entropy and diffusive noise strength can be neglected in this limit. The relevance to biological, ecological, and climate models is apparent.

Stochastic thermodynamics and entropy production of chemical reaction systems

We investigate the nonequilibrium stationary states of systems consisting of chemical reactions among molecules of several chemical species. To this end, we introduce and develop a stochastic formulation of nonequilibrium thermodynamics of chemical reaction systems based on a master equation defined on the space of microscopic chemical states and on appropriate definitions of entropy and entropy production. The system is in contact with a heat reservoir and is placed out of equilibrium by the contact with particle reservoirs. In our approach, the fluxes of various types, such as the heat and particle fluxes, play a fundamental role in characterizing the nonequilibrium chemical state. We show that the rate of entropy production in the stationary nonequilibrium state is a bilinear form in the affinities and the fluxes of reaction, which are expressed in terms of rate constants and transition rates, respectively. We also show how the description in terms of microscopic states can be reduced to a description in terms of the numbers of particles of each species, from which follows the chemical master equation. As an example, we calculate the rate of entropy production of the first and second Schlögl reaction models.

Entrance and escape dynamics for the typical set

According to the asymptotic equipartition property, sufficiently long sequences of random variables converge to a set that is typical. While the size and probability of this set are central to information theory and statistical mechanics, they can often only be estimated accurately in the asymptotic limit due to the exponential growth in possible sequences. Here we derive a time-inhomogeneous dynamics that constructs the properties of the typical set for all finite length sequences of independent and identically distributed random variables. These dynamics link the finite properties of the typical set to asymptotic results and allow the typical set to be applied to small and transient systems. The main result is a geometric mapping-the triangle map-relating sequences of random variables of length n to those of length n+1. We show that the number of points in this map needed to quantify the properties of the typical set grows linearly with sequence length, despite the exponential growth in the number of typical sequences. We illustrate the framework for the Bernoulli process and the Schlögl model for autocatalytic chemical reactions and demonstrate both the convergence to asymptotic limits and the ability to reproduce exact calculations.

Entropy production selects nonequilibrium states in multistable systems

Far-from-equilibrium thermodynamics underpins the emergence of life, but how has been a long-outstanding puzzle. Best candidate theories based on the maximum entropy production principle could not be unequivocally proven, in part due to complicated physics, unintuitive stochastic thermodynamics, and the existence of alternative theories such as the minimum entropy production principle. Here, we use a simple, analytically solvable, one-dimensional bistable chemical system to demonstrate the validity of the maximum entropy production principle. To generalize to multistable stochastic system, we use the stochastic least-action principle to derive the entropy production and its role in the stability of nonequilibrium steady states. This shows that in a multistable system, all else being equal, the steady state with the highest entropy production is favored, with a number of implications for the evolution of biological, physical, and geological systems.

A minimal model of burst-noise induced bistability

We investigate the influence of intrinsic noise on stable states of a one-dimensional dynamical system that shows in its deterministic version a saddle-node bifurcation between monostable and bistable behaviour. The system is a modified version of the Schlögl model, which is a chemical reaction system with only one type of molecule. The strength of the intrinsic noise is varied without changing the deterministic description by introducing bursts in the autocatalytic production step. We study the transitions between monostable and bistable behavior in this system by evaluating the number of maxima of the stationary probability distribution. We find that changing the size of bursts can destroy and even induce saddle-node bifurcations. This means that a bursty production of molecules can qualitatively change the dynamics of a chemical reaction system even when the deterministic description remains unchanged.

Effects of cell cycle noise on excitable gene circuits

We assess the impact of cell cycle noise on gene circuit dynamics. For bistable genetic switches and excitable circuits, we find that transitions between metastable states most likely occur just after cell division and that this concentration effect intensifies in the presence of transcriptional delay. We explain this concentration effect with a three-states stochastic model. For genetic oscillators, we quantify the temporal correlations between daughter cells induced by cell division. Temporal correlations must be captured properly in order to accurately quantify noise sources within gene networks.