Landscape and global stability of nonadiabatic and adiabatic oscillations in a gene network

We quantify the potential landscape to determine the global stability and coherence of biological oscillations. We explore a gene network motif in our experimental synthetic biology studies of two genes that mutually repress and activate each other with self-activation and self-repression. We find that in addition to intrinsic molecular number fluctuations, there is another type of fluctuation crucial for biological function: the fluctuation due to the slow binding/unbinding of protein regulators to gene promoters. We find that coherent limit cycle oscillations emerge in two regimes: an adiabatic regime with fast binding/unbinding and a nonadiabatic regime with slow binding/unbinding relative to protein synthesis/degradation. This leads to two mechanisms of producing the stable oscillations: the effective interactions from averaging the gene states in the adiabatic regime; and the time delays due to slow binding/unbinding to promoters in the nonadiabatic regime, which can be tested by forthcoming experiments. In both regimes, the landscape has a topological shape of the Mexican hat in protein concentrations that quantitatively determines the global stability of limit cycle dynamics. The oscillation coherence is shown to be correlated with the shape of the Mexican hat characterized by the height from the oscillation ring to the central top. The oscillation period can be tuned in a wide range by changing the binding/unbinding rate without changing the amplitude much, which is important for the functionality of a biological clock. A negative feedback loop with time delays due to slow binding/unbinding can also generate oscillations. Although positive feedback is not necessary for generating oscillations, it can make the oscillations more robust.     

Conformational switching upon phosphorylation: a predictive framework based on energy landscape principles

We investigate how post-translational phosphorylation modifies the global conformation of a protein by changing its free energy landscape using two test proteins, cystatin and NtrC. We first examine the changes in a free energy landscape caused by phosphorylation using a model containing information about both structural forms. For cystatin the free energy cost is fairly large indicating a low probability of sampling the phosphorylated conformation in a perfectly funneled landscape. The predicted barrier for NtrC conformational transition is several times larger than the barrier for cystatin, indicating that the switch protein NtrC most probably follows a partial unfolding mechanism to move from one basin to the other. Principal component analysis and linear response theory show how the naturally occurring conformational changes in unmodified proteins are captured and stabilized by the change of interaction potential. We also develop a partially guided structure prediction Hamiltonian which is capable of predicting the global structure of a phosphorylated protein using only knowledge of the structure of the unphosphorylated protein or vice versa. This algorithm makes use of a generic transferable long-range residue contact potential along with details of structure short range in sequence. By comparing the results obtained with this guided transferable potential to those from the native-only, perfectly funneled Hamiltonians, we show that the transferable Hamiltonian correctly captures the nature of the global conformational changes induced by phosphorylation and can sample substantially correct structures for the modified protein with high probability.     

Potential landscape and probabilistic flux of a predator prey network

Predator-prey system, as an essential element of ecological dynamics, has been recently studied experimentally with synthetic biology. We developed a global probabilistic landscape and flux framework to explore a synthetic predator-prey network constructed with two Escherichia coli populations. We developed a self consistent mean field method to solve multidimensional problem and uncovered the potential landscape with Mexican hat ring valley shape for predator-prey oscillations. The landscape attracts the system down to the closed oscillation ring. The probability flux drives the coherent oscillations on the ring. Both the landscape and flux are essential for the stable and coherent oscillations. The landscape topography characterized by the barrier height from the top of Mexican hat to the closed ring valley provides a quantitative measure of global stability of system. The entropy production rate for the energy dissipation is less for smaller environmental fluctuations or perturbations. The global sensitivity analysis based on the landscape topography gives specific predictions for the effects of parameters on the stability and function of the system. This may provide some clues for the global stability, robustness, function and synthetic network design.     

Noise-reduction through interaction in gene expression and biochemical reaction processes

We demonstrate that interaction in gene expression and biochemical reaction processes has a significant influence on reducing fluctuations. Especially, we have found that the interaction between synthesized proteins and background molecules can reduce the fluctuation level in gene expression, which is a counter example to the intuition that background factors disturb information processing in genetic networks by increasing the noise level. This fact also indicates that the macromolecular crowding observed in actual cells can contribute to reduce the noise level. In addition, the noise-reduction phenomenon is not limited to the interaction between the proteins and background molecules, but can be applied to other reactions such as a dimerization process and the coupling of reactions with large fluctuations by intrinsic noise. Finally, on the basis of these results, we propose a new and plausible method for reducing the fluctuations generated in synthesized genetic networks, and also discuss the applicability of this method to the stabilization of system dynamics by using a toggle switch model.     

A novel approach to decoy set generation: designing a physical energy function having local minima with native structure characteristics

We suggest a new approach to the generation of candidate structures (decoys) for ab initio prediction of protein structures. Our method is based on random sampling of conformation space and subsequent local energy minimization. At the core of this approach lies the design of a novel type of energy function. This energy function has local minima with native structure characteristics and wide basins of attraction. The current work presents our motivation for deriving such an energy function and also tests the derived energy function.Our approach is novel in that it takes advantage of the inherently rough energy landscape of proteins, which is generally considered a major obstacle for protein structure prediction. When local minima have wide basins of attraction, the protein's conformation space can be greatly reduced by the convergence of large regions of the space into single points, namely the local minima corresponding to these funnels. We have implemented this concept by an iterative process. The potential is first used to generate decoy sets and then we study these sets of decoys to guide further development of the potential. A key feature of our potential is the use of cooperative multi-body interactions that mimic the role of the entropic and solvent contributions to the free energy.The validity and value of our approach is demonstrated by applying it to 14 diverse, small proteins. We show that, for these proteins, the size of conformation space is considerably reduced by the new energy function. In fact, the reduction is so substantial as to allow efficient conformational sampling. As a result we are able to find a significant number of near-native conformations in random searches performed with limited computational resources.     

A dynamical systems hypothesis of schizophrenia

We propose a top-down approach to the symptoms of schizophrenia based on a statistical dynamical framework. We show that a reduced depth in the basins of attraction of cortical attractor states destabilizes the activity at the network level due to the constant statistical fluctuations caused by the stochastic spiking of neurons. In integrate-and-fire network simulations, a decrease in the NMDA receptor conductances, which reduces the depth of the attractor basins, decreases the stability of short-term memory states and increases distractibility. The cognitive symptoms of schizophrenia such as distractibility, working memory deficits, or poor attention could be caused by this instability of attractor states in prefrontal cortical networks. Lower firing rates are also produced, and in the orbitofrontal and anterior cingulate cortex could account for the negative symptoms, including a reduction of emotions. Decreasing the GABA as well as the NMDA conductances produces not only switches between the attractor states, but also jumps from spontaneous activity into one of the attractors. We relate this to the positive symptoms of schizophrenia, including delusions, paranoia, and hallucinations, which may arise because the basins of attraction are shallow and there is instability in temporal lobe semantic memory networks, leading thoughts to move too freely round the attractor energy landscape.     

Funneled landscape leads to robustness of cellular networks: MAPK signal transduction

We uncover the underlying potential energy landscape for a cellular network. We find that the potential energy landscape of the mitogen-activated protein-kinase signal transduction network is funneled toward the global minimum. The funneled landscape is quite robust against random perturbations. This naturally explains robustness from a physical point of view. The ratio of slope versus roughness of the landscape becomes a quantitative measure of robustness of the network. Funneled landscape is a realization of the Darwinian principle of natural selection at the cellular network level. It provides an optimal criterion for network connections and design. Our approach is general and can be applied to other cellular networks.     

Exploring the Free Energy Landscape: From Dynamics to Networks and Back

Knowledge of the Free Energy Landscape topology is the essential key to understanding many biochemical processes. The determination of the conformers of a protein and their basins of attraction takes a central role for studying molecular isomerization reactions. In this work, we present a novel framework to unveil the features of a Free Energy Landscape answering questions such as how many meta-stable conformers there are, what the hierarchical relationship among them is, or what the structure and kinetics of the transition paths are. Exploring the landscape by molecular dynamics simulations, the microscopic data of the trajectory are encoded into a Conformational Markov Network. The structure of this graph reveals the regions of the conformational space corresponding to the basins of attraction. In addition, handling the Conformational Markov Network, relevant kinetic magnitudes as dwell times and rate constants, or hierarchical relationships among basins, completes the global picture of the landscape. We show the power of the analysis studying a toy model of a funnel-like potential and computing efficiently the conformers of a short peptide, dialanine, paving the way to a systematic study of the Free Energy Landscape in large peptides.     

Landscape, flux, correlation, resonance, coherence, stability, and key network wirings of stochastic circadian oscillation

Circadian rhythms with a period of ∼24 h, are natural timing machines. They are broadly distributed in living organisms, such as Neurospora, Drosophila, and mammals. The underlying natures of the rhythmic behavior have been explored by experimental and theoretical approaches. However, the global and physical natures of the oscillation under fluctuations are still not very clear. We developed a landscape and flux framework to explore the global stability and robustness of a circadian oscillation system. The potential landscape of the network is uncovered and has a global Mexican-hat shape. The height of the Mexican-hat provides a quantitative measure to evaluate the robustness and coherence of the oscillation. We found that in nonequilibrium dynamic systems, not only the potential landscape but also the probability flux are important to the dynamics of the system under intrinsic noise. Landscape attracts the systems down to the oscillation ring while flux drives the coherent oscillation on the ring. We also investigated the phase coherence and the entropy production rate of the system at different fluctuations and found that dissipations are less and the coherence is higher for larger number of molecules. We also found that the power spectrum of autocorrelation functions show resonance peak at the frequency of coherent oscillations. The peak is less prominent for smaller number of molecules and less barrier height and therefore can be used as another measure of stability of oscillations. As a consequence of nonzero probability flux, we show that the three-point correlations from the time traces show irreversibility, providing a possible way to explore the flux from the observations. Furthermore, we explored the escape time from the oscillation ring to outside at different molecular number. We found that when barrier height is higher, escape time is longer and phase coherence of oscillation is higher. Finally, we performed the global sensitivity analysis of the underlying parameters to find the key network wirings responsible for the stability of the oscillation system.     

Quantifying Cell Fate Decisions for Differentiation and Reprogramming of a Human Stem Cell Network: Landscape and Biological Paths

Cellular reprogramming has been recently intensively studied experimentally. We developed a global potential landscape and kinetic path framework to explore a human stem cell developmental network composed of 52 genes. We uncovered the underlying landscape for the stem cell network with two basins of attractions representing stem and differentiated cell states, quantified and exhibited the high dimensional biological paths for the differentiation and reprogramming process, connecting the stem cell state and differentiated cell state. Both the landscape and non-equilibrium curl flux determine the dynamics of cell differentiation jointly. Flux leads the kinetic paths to be deviated from the steepest descent gradient path, and the corresponding differentiation and reprogramming paths are irreversible. Quantification of paths allows us to find out how the differentiation and reprogramming occur and which important states they go through. We show the developmental process proceeds as moving from the stem cell basin of attraction to the differentiation basin of attraction. The landscape topography characterized by the barrier heights and transition rates quantitatively determine the global stability and kinetic speed of cell fate decision process for development. Through the global sensitivity analysis, we provided some specific predictions for the effects of key genes and regulation connections on the cellular differentiation or reprogramming process. Key links from sensitivity analysis and biological paths can be used to guide the differentiation designs or reprogramming tactics.     

Robustness and Coherence of a Three-Protein Circadian Oscillator: Landscape and Flux Perspectives

Three-protein circadian oscillations in cyanobacteria sustain for weeks. To understand how cellular oscillations function robustly in stochastic fluctuating environments, we used a stochastic model to uncover two natures of circadian oscillation: the potential landscape related to steady-state probability distribution of protein concentrations; and the corresponding flux related to speed of concentration changes which drive the oscillations. The barrier height of escaping from the oscillation attractor on the landscape provides a quantitative measure of the robustness and coherence for oscillations against intrinsic and external fluctuations. The difference between the locations of the zero total driving force and the extremal of the potential provides a possible experimental probe and quantification of the force from curl flux. These results, correlated with experiments, can help in the design of robust oscillatory networks.     

Effects of intermediate bound states in dynamic force spectroscopy

We revisit some aspects of the interpretation of dynamic force spectroscopy experiments. The standard theory predicts that the typical unbinding force f* is linearly proportional to the logarithm of the loading rate r when a single energy barrier controls the unbinding process. For a more complex situation of N barriers, it predicts at most N linear segments for the f* vs. log(r) curve, each segment characterizing a different barrier. Here we extend this existing picture using a refined approximation, provide a more general analytical formula, and show that in principle up to N(N + 1) / 2 segments can show up experimentally. As a consequence, the determination of the positions and even the number of the energy barriers from the experimental data can be ambiguous. A further possible consequence of a multiple-barrier landscape is a bimodal or multimodal distribution of the unbinding force at certain loading rates, a feature recently observed experimentally.     

The influence of intrinsic folding mechanism of an unfolded protein on the coupled folding-binding process during target recognition

Intrinsically disordered proteins (IDPs) are extensively involved in dynamic signaling processes which require a high association rate and a high dissociation rate for rapid binding/unbinding events and at the same time a sufficient high affinity for specific recognition. Although the coupled folding-binding processes of IDPs have been extensively studied, it is still impossible to predict whether an unfolded protein is suitable for molecular signaling via coupled folding-binding. In this work, we studied the interplay between intrinsic folding mechanisms and coupled folding-binding process for unfolded proteins through molecular dynamics simulations. We first studied the folding process of three representative IDPs with different folded structures, that is, c-Myb, AF9, and E3 rRNase. We found the folding free energy landscapes of IDPs are downhill or show low barriers. To further study the influence of intrinsic folding mechanism on the binding process, we modulated the folding mechanism of barnase via circular permutation and simulated the coupled folding-binding process between unfolded barnase permutant and folded barstar. Although folding of barnase was coupled to target binding, the binding kinetics was significantly affected by the intrinsic folding free energy barrier, where reducing the folding free energy barrier enhances binding rate up to two orders of magnitude. This accelerating effect is different from previous results which reflect the effect of structure flexibility on binding kinetics. Our results suggest that coupling the folding of an unfolded protein with no/low folding free energy barrier with its target binding may provide a way to achieve high specificity and rapid binding/unbinding kinetics simultaneously.