Gene regulatory networks: a coarse-grained, equation-free approach to multiscale computation

We present computer-assisted methods for analyzing stochastic models of gene regulatory networks. The main idea that underlies this equation-free analysis is the design and execution of appropriately-initialized short bursts of stochastic simulations; the results of these are processed to estimate coarse-grained quantities of interest, such as mesoscopic transport coefficients. In particular, using a simple model of a genetic toggle switch, we illustrate the computation of an effective free energy and of a state-dependent effective diffusion coefficient that characterize an unavailable effective Fokker-Planck equation. Additionally we illustrate the linking of equation-free techniques with continuation methods for performing a form of stochastic "bifurcation analysis"; estimation of mean switching times in the case of a bistable switch is also implemented in this equation-free context. The accuracy of our methods is tested by direct comparison with long-time stochastic simulations. This type of equation-free analysis appears to be a promising approach to computing features of the long-time, coarse-grained behavior of certain classes of complex stochastic models of gene regulatory networks, circumventing the need for long Monte Carlo simulations.Comment: 33 pages, submitted to The Journal of Chemical Physic

Megawatt solar power systems for lunar surface operations

The work presented here shows that a solar power system can provide power on the order of one megawatt to a lunar base with a fairly high specific power. The main drawback to using solar power is still the high mass, and therefore, cost of supplying energy storage through the solar night. The use of cryogenic reactant storage in a fuel cell system, however, greatly reduces the total system mass over conventional energy storage schemes

Stronger computational modelling of signalling pathways using both continuous and discrete-state methods

Starting from a biochemical signalling pathway model expresses in a process algebra enriched with quantitative information, we automatically derive both continuous-space and discrete-space representations suitable for numerical evaluation. We compare results obtained using approximate stochastic simulation thereby exposing a flaw in the use of the differentiation procedure producing misleading results

Sampling rare switching events in biochemical networks

Bistable biochemical switches are ubiquitous in gene regulatory networks and signal transduction pathways. Their switching dynamics, however, are difficult to study directly in experiments or conventional computer simulations, because switching events are rapid, yet infrequent. We present a simulation technique that makes it possible to predict the rate and mechanism of flipping of biochemical switches. The method uses a series of interfaces in phase space between the two stable steady states of the switch to generate transition trajectories in a ratchet-like manner. We demonstrate its use by calculating the spontaneous flipping rate of a symmetric model of a genetic switch consisting of two mutually repressing genes. The rate constant can be obtained orders of magnitude more efficiently than using brute-force simulations. For this model switch, we show that the switching mechanism, and consequently the switching rate, depends crucially on whether the binding of one regulatory protein to the DNA excludes the binding of the other one. Our technique could also be used to study rare events and non-equilibrium processes in soft condensed matter systems.Comment: 9 pages, 6 figures, last page contains supplementary informatio

Individual Differences and Stability of Dynamics among Self-Concept Clarity, Impatience, and Negative Affect

Self-concept clarity (SCC) is associated with behavioral and emotion regulation, although the nature of this link is unclear. SCC may serve as a self-regulatory resource or it may be a product of well-regulated behaviors and emotions. In two studies using experience sampling among undergraduates (n = 46 and n = 36), we investigate whether models representing relationships among SCC, impatience, and negative affect (NA) states conform to these theories, are similar across individuals, and are stable across a one-month period. Results reveal substantial variation between persons in these dynamic relationships, suggesting that multiple SCC-relevant regulatory processes exist. These patterns were not stable from one month to the next, but changes in them related to changes in stress, suggesting higher-order regulation of these dynamics

Number of adaptive steps to a local fitness peak

We consider a population of genotype sequences evolving on a rugged fitness landscape with many local fitness peaks. The population walks uphill until it encounters a local fitness maximum. We find that the statistical properties of the walk length depend on whether the underlying fitness distribution has a finite mean. If the mean is finite, all the walk length cumulants grow with the sequence length but approach a constant otherwise. Experimental implications of our analytical results are also discussed

Analytical study of non Gaussian fluctuations in a stochastic scheme of autocatalytic reactions

A stochastic model of autocatalytic chemical reactions is studied both numerically and analytically. The van Kampen perturbative scheme is implemented, beyond the second order approximation, so to capture the non Gaussianity traits as displayed by the simulations. The method is targeted to the characterization of the third moments of the distribution of fluctuations, originating from a system of four populations in mutual interaction. The theory predictions agree well with the simulations, pointing to the validity of the van Kampen expansion beyond the conventional Gaussian solution.Comment: 15 pages, 8 figures, submitted to Phys. Rev.