Adaptive stochastic-deterministic chemical kinetic simulations

MOTIVATION: Biochemical signaling pathways and genetic circuits often involve very small numbers of key signaling molecules. Computationally expensive stochastic methods are necessary to simulate such chemical situations. Single-molecule chemical events often co-exist with much larger numbers of signaling molecules where mass-action kinetics is a reasonable approximation. Here, we describe an adaptive stochastic method that dynamically chooses between deterministic and stochastic calculations depending on molecular count and propensity of forward reactions. The method is fixed timestep and has first order accuracy. We compare the efficiency of this method with exact stochastic methods. RESULTS: We have implemented an adaptive stochastic-deterministic approximate simulation method for chemical kinetics. With an error margin of 5%, the method solves typical biologically constrained reaction schemes more rapidly than exact stochastic methods for reaction volumes >1-10 micro m(3). We have developed a test suite of reaction cases to test the accuracy of mixed simulation methods. AVAILABILITY: Simulation software used in the paper is freely available from http://www.ncbs.res.in/kinetikit/download.html     

STOCHSIM: modelling of stochastic biomolecular processes

SUMMARY: STOCHSIM is a stochastic simulator for chemical reactions. Molecules are represented as individual software objects that react according to probabilities derived from concentrations and rate constants. Version 1.2 of STOCHSIM provides a novel cross-platform graphical interface written in Perl/Tk. A simple two-dimensional spatial structure has also been implemented, in which nearest-neighbour interactions of molecules in a 2-D lattice can be simulated.     

Cyto-Sim: a formal language model and stochastic simulator of membrane-enclosed biochemical processes

MOTIVATION: Compartments and membranes are the basis of cell topology and more than 30% of the human genome codes for membrane proteins. While it is possible to represent compartments and membrane proteins in a nominal way with many mathematical formalisms used in systems biology, few, if any, explicitly model the topology of the membranes themselves. Discrete stochastic simulation potentially offers the most accurate representation of cell dynamics. Since the details of every molecular interaction in a pathway are often not known, the relationship between chemical species in not necessarily best described at the lowest level, i.e. by mass action. Simulation is a form of computer-aided analysis, relying on human interpretation to derive meaning. To improve efficiency and gain meaning in an automatic way, it is necessary to have a formalism based on a model which has decidable properties. RESULTS: We present Cyto-Sim, a stochastic simulator of membrane-enclosed hierarchies of biochemical processes, where the membranes comprise an inner, outer and integral layer. The underlying model is based on formal language theory and has been shown to have decidable properties (Cavaliere and Sedwards, 2006), allowing formal analysis in addition to simulation. The simulator provides variable levels of abstraction via arbitrary chemical kinetics which link to ordinary differential equations. In addition to its compact native syntax, Cyto-Sim currently supports models described as Petri nets, can import all versions of SBML and can export SBML and MATLAB m-files. AVAILABILITY: Cyto-Sim is available free, either as an applet or a stand-alone Java program via the web page (http://www.cosbi.eu/Rpty_Soft_CytoSim.php). Other versions can be made available upon request.     

Differential Evolution of Advertisement Call Traits in Dart-Poison Frogs (Anura: Dendrobatidae)

The ability to recognise and discriminate between heterospecific and conspecific individuals plays an essential role in mate choice, reproductive isolation and thus species diversification. Many animals discriminate based on advertisement calls, whose evolution may be driven by a variety of forces such as natural selection, sexual selection or stochastic processes. The relative importance of stochastic processes acting on a given trait is usually correlated with its phylogenetic signal. Mate-recognition signals are complex traits composed of multiple features that could potentially respond independently to evolutionary forces. The advertisement call of anurans is used in species recognition and mate choice. In this study, we estimate the phylogenetic signal for body size and a suite of traits describing the male advertisement call from dart-poison frogs (Anura: Dendrobatidae). We found a surprisingly high phylogenetic signal for all call traits. In addition, call traits varied in their degree of phylogenetic signal, suggesting that evolutionary forces have been acting differently on different traits. Pulse duration showed the strongest phylogenetic signal. Peak frequency and body size were correlated and presented high phylogenetic signal indicating that the evolution of one trait may be driving or constraining the other. Since most variation in call traits can be explained by the phylogenetic history of the species, we cannot reject the hypothesis that stochastic processes account for significant evolutionary divergence in frog calls.     

StochPy: A Comprehensive,User-Friendly Tool for Simulating Stochastic Biological Processes

Single-cell and single-molecule measurements indicate the importance of stochastic phenomena in cell biology. Stochasticity creates spontaneous differences in the copy numbers of key macromolecules and the timing of reaction events between genetically-identical cells. Mathematical models are indispensable for the study of phenotypic stochasticity in cellular decision-making and cell survival. There is a demand for versatile, stochastic modeling environments with extensive, preprogrammed statistics functions and plotting capabilities that hide the mathematics from the novice users and offers low-level programming access to the experienced user. Here we present StochPy (Stochastic modeling in Python), which is a flexible software tool for stochastic simulation in cell biology. It provides various stochastic simulation algorithms, SBML support, analyses of the probability distributions of molecule copy numbers and event waiting times, analyses of stochastic time series, and a range of additional statistical functions and plotting facilities for stochastic simulations. We illustrate the functionality of StochPy with stochastic models of gene expression, cell division, and single-molecule enzyme kinetics. StochPy has been successfully tested against the SBML stochastic test suite, passing all tests. StochPy is a comprehensive software package for stochastic simulation of the molecular control networks of living cells. It allows novice and experienced users to study stochastic phenomena in cell biology. The integration with other Python software makes StochPy both a user-friendly and easily extendible simulation tool.     

A coalescent framework for comparing alternative models of population structure with genetic data: evolution of Celebes toads

Isolation of populations eventually leads to divergence by genetic drift, but if connectivity varies over time, its impact on diversification may be difficult to discern. Even when the habitat patches of multiple species overlap, differences in their demographic parameters, molecular evolution and stochastic events contribute to differences in the magnitude and distribution of their genetic variation. The Indonesian island of Sulawesi, for example, harbours a suite of endemic species whose intraspecific differentiation or interspecific divergence may have been catalysed by habitat fragmentation. To further test this hypothesis, we have performed phylogenetic and coalescent-based analyses on molecular variation in mitochondrial and nuclear DNA of the Celebes toad (Bufo celebensis). Results support a role for habitat fragmentation that led to a population structure in these toads that closely matches distributions of Sulawesi macaque monkeys. Habitat fragmentation, therefore, may also have affected other groups on this island.     

Stochastic resonance in a single neuron model: theory and analog simulation.

Here, we consider a noisy, bistable, single neuron model in the presence of periodic external modulation. The modulation induces a correlated switching between states driven by the noise. The information flow through the system, from the modulation, or signal, to the output switching events, leads to a succession of strong peaks in the power spectrum. The signal-to-noise ratio (SNR) obtained from this power spectrum is a measure of the information content in the neuron response. With increasing noise intensity, the SNR passes through a maximum: an effect which has been called stochastic resonance, and which was first advanced as a possible explanation of the observed periodicity in the recurrences of the Earth's ice ages. We treat the problem within the framework of a recently developed approximate theory, valid in the limits of weak noise intensity, weak periodic forcing and low forcing frequency, for both additive and multiplicative noise. Moreover, we have constructed an analog simulator of the neuron which demonstrates the stochastic resonance effect, and with which we have measured the SNRs for comparison with the theoretical results. Our model should be of interest in situations where a single inherently noisy neuron is the receptor of a periodic signal, which is itself noisy, either from the network or from an external source.     

Bayesian History Matching of Complex Infectious Disease Models Using Emulation: A Tutorial and a Case Study on HIV in Uganda

Advances in scientific computing have allowed the development of complex models that are being routinely applied to problems in disease epidemiology, public health and decision making. The utility of these models depends in part on how well they can reproduce empirical data. However, fitting such models to real world data is greatly hindered both by large numbers of input and output parameters, and by long run times, such that many modelling studies lack a formal calibration methodology. We present a novel method that has the potential to improve the calibration of complex infectious disease models (hereafter called simulators). We present this in the form of a tutorial and a case study where we history match a dynamic, event-driven, individual-based stochastic HIV simulator, using extensive demographic, behavioural and epidemiological data available from Uganda. The tutorial describes history matching and emulation. History matching is an iterative procedure that reduces the simulator''s input space by identifying and discarding areas that are unlikely to provide a good match to the empirical data. History matching relies on the computational efficiency of a Bayesian representation of the simulator, known as an emulator. Emulators mimic the simulator''s behaviour, but are often several orders of magnitude faster to evaluate. In the case study, we use a 22 input simulator, fitting its 18 outputs simultaneously. After 9 iterations of history matching, a non-implausible region of the simulator input space was identified that was times smaller than the original input space. Simulator evaluations made within this region were found to have a 65% probability of fitting all 18 outputs. History matching and emulation are useful additions to the toolbox of infectious disease modellers. Further research is required to explicitly address the stochastic nature of the simulator as well as to account for correlations between outputs.     

COPASI--a COmplex PAthway SImulator

MOTIVATION: Simulation and modeling is becoming a standard approach to understand complex biochemical processes. Therefore, there is a big need for software tools that allow access to diverse simulation and modeling methods as well as support for the usage of these methods. RESULTS: Here, we present COPASI, a platform-independent and user-friendly biochemical simulator that offers several unique features. We discuss numerical issues with these features; in particular, the criteria to switch between stochastic and deterministic simulation methods, hybrid deterministic-stochastic methods, and the importance of random number generator numerical resolution in stochastic simulation. AVAILABILITY: The complete software is available in binary (executable) for MS Windows, OS X, Linux (Intel) and Sun Solaris (SPARC), as well as the full source code under an open source license from http://www.copasi.org.     

Modeling biological systems using Dynetica--a simulator of dynamic networks

We present Dynetica, a user-friendly simulator of dynamic networks for constructing, visualizing, and analyzing kinetic models of biological systems. In addition to generic reaction networks, Dynetica facilitates construction of models of genetic networks, where many reactions are gene expression and interactions among gene products. Further, it integrates the capability of conducting both deterministic and stochastic simulations. AVAILABILITY AND SUPPLEMENTARY INFORMATION: Dynetica 1.0, example models, and the user's guide are available at http://www.its.caltech.edu/~you/Dynetica/Dynetica_page.htm     

Multi-dimensional, mesoscopic Monte Carlo simulations of inhomogeneous reaction-drift-diffusion systems on graphics-processing units

For many biological applications, a macroscopic (deterministic) treatment of reaction-drift-diffusion systems is insufficient. Instead, one has to properly handle the stochastic nature of the problem and generate true sample paths of the underlying probability distribution. Unfortunately, stochastic algorithms are computationally expensive and, in most cases, the large number of participating particles renders the relevant parameter regimes inaccessible. In an attempt to address this problem we present a genuine stochastic, multi-dimensional algorithm that solves the inhomogeneous, non-linear, drift-diffusion problem on a mesoscopic level. Our method improves on existing implementations in being multi-dimensional and handling inhomogeneous drift and diffusion. The algorithm is well suited for an implementation on data-parallel hardware architectures such as general-purpose graphics processing units (GPUs). We integrate the method into an operator-splitting approach that decouples chemical reactions from the spatial evolution. We demonstrate the validity and applicability of our algorithm with a comprehensive suite of standard test problems that also serve to quantify the numerical accuracy of the method. We provide a freely available, fully functional GPU implementation. Integration into Inchman, a user-friendly web service, that allows researchers to perform parallel simulations of reaction-drift-diffusion systems on GPU clusters is underway.     

Spatial stochastic modelling of the phosphoenolpyruvate-dependent phosphotransferase (PTS) pathway in Escherichia coli

MOTIVATION: Many biochemical networks involve reactions localized on the cell membrane. This can give rise to spatial gradients of the concentration of cytosolic species. Moreover, the number of membrane molecules can be small and stochastic effects can become relevant. Pathways usually consist of a complex interaction network and are characterized by a large set of parameters. The inclusion of spatial and stochastic effects is a major challenge in developing quantitative and dynamic models of pathways. RESULTS: We have developed a particle-based spatial stochastic method (GMP) to simulate biochemical networks in space, including fluctuations from the diffusion of particles and reactions. Gradients emerging from membrane reactions can be resolved. As case studies for the GMP method we used a simple gene expression system and the phosphoenolpyruvate:glucose phosphotransferase system pathway. AVAILABILITY: The source code for the GMP method is available at http://www.science.uva.nl/research/scs/CellMath/GMP.     

One-hit stochastic decline in a mechanochemical model of cytoskeleton-induced neuron death III: diffusion pulse death zones

This is the third of three papers in which we study a mathematical model of cytoskeleton-induced neuron death. In the first two papers of this suite [Lomasko, T., Clarke, G., Lumsden, C., 2007a. One-hit stochastic decline in a mechanochemical model of cytoskeleton-induced neuron death I: cell fate arrival times. J. Theor. Biol. 249, 1-17, doi:10.1016/j.jtbi.2007.05.031; Lomasko, T., Clarke, G., Lumsden, C., 2007b. One-hit stochastic decline in a mechanochemical model of cytoskeleton-induced neuron death II: transition state metastability. J. Theor. Biol. 249, 18-28, doi:10.1016/j.jtbi.2007.05.032], we established that the mean-field limit of our model relates the known patterns of neuron decline to specific scales of cytoskeleton reorganization and cell-cell interaction by diffusible death factors. In the mean-field limit, the spatially variable concentration of diffusing death factor is replaced by a constant average value. Recent empirical advances now permit the actual diffusion of such factors to be followed in intact neuropil. In this paper we therefore extend the model beyond the mean-field limit, to include the diffusion dynamics of death factor bursts released from dying neurons. A range of novel tissue degeneration patterns is observed, for which we confirm and extend the mean-field prediction that sigmoidal patterns of neuron population decay are a principal hallmark of cell death in the presence of death factor release.     

Comparative gene marker selection suite

MOTIVATION: An important step in analyzing expression profiles from microarray data is to identify genes that can discriminate between distinct classes of samples. Many statistical approaches for assigning significance values to genes have been developed. The Comparative Marker Selection suite consists of three modules that allow users to apply and compare different methods of computing significance for each marker gene, a viewer to assess the results, and a tool to create derivative datasets and marker lists based on user-defined significance criteria. AVAILABILITY: The Comparative Marker Selection application suite is freely available as a GenePattern module. The GenePattern analysis environment is freely available at http://www.broad.mit.edu/genepattern.     

A cooperative particle swarm optimizer with migration of heterogeneous probabilistic models

Particle Swarm Optimization (PSO) is a stochastic optimization approach that originated from simulations of bird flocking, and that has been successfully used in many applications as an optimization tool. Estimation of distribution algorithms (EDAs) are a class of evolutionary algorithms which perform a two-step process: building a probabilistic model from which good solutions may be generated and then using this model to generate new individuals. Two distinct research trends that emerged in the past few years are the hybridization of PSO and EDA algorithms and the parallelization of EDAs to exploit the idea of exchanging the probabilistic model information. In this work, we propose the use of a cooperative PSO/EDA algorithm based on the exchange of heterogeneous probabilistic models. The model is heterogeneous because the cooperating PSO/EDA algorithms use different methods to sample the search space. Three different exchange approaches are tested and compared in this work. In all these approaches, the amount of information exchanged is adapted based on the performance of the two cooperating swarms. The performance of the cooperative model is compared to the existing state-of-the-art PSO cooperative approaches using a suite of well-known benchmark optimization functions.     

A stochastic simulator of a blood product donation environment with demand spikes and supply shocks

The availability of an adequate blood supply is a critical public health need. An influenza epidemic or another crisis affecting population mobility could create a critical donor shortage, which could profoundly impact blood availability. We developed a simulation model for the blood supply environment in the United States to assess the likely impact on blood availability of factors such as an epidemic. We developed a simulator of a multi-state model with transitions among states. Weekly numbers of blood units donated and needed were generated by negative binomial stochastic processes. The simulator allows exploration of the blood system under certain conditions of supply and demand rates, and can be used for planning purposes to prepare for sudden changes in the public's health. The simulator incorporates three donor groups (first-time, sporadic, and regular), immigration and emigration, deferral period, and adjustment factors for recruitment. We illustrate possible uses of the simulator by specifying input values for an 8-week flu epidemic, resulting in a moderate supply shock and demand spike (for example, from postponed elective surgeries), and different recruitment strategies. The input values are based in part on data from a regional blood center of the American Red Cross during 1996-2005. Our results from these scenarios suggest that the key to alleviating deficit effects of a system shock may be appropriate timing and duration of recruitment efforts, in turn depending critically on anticipating shocks and rapidly implementing recruitment efforts.     

The Importance of Stimulus Noise Analysis for Self-Motion Studies

Motion simulators are widely employed in basic and applied research to study the neural mechanisms of perception and action during inertial stimulation. In these studies, uncontrolled simulator-introduced noise inevitably leads to a disparity between the reproduced motion and the trajectories meticulously designed by the experimenter, possibly resulting in undesired motion cues to the investigated system. Understanding actual simulator responses to different motion commands is therefore a crucial yet often underestimated step towards the interpretation of experimental results. In this work, we developed analysis methods based on signal processing techniques to quantify the noise in the actual motion, and its deterministic and stochastic components. Our methods allow comparisons between commanded and actual motion as well as between different actual motion profiles. A specific practical example from one of our studies is used to illustrate the methodologies and their relevance, but this does not detract from its general applicability. Analyses of the simulator’s inertial recordings show direction-dependent noise and nonlinearity related to the command amplitude. The Signal-to-Noise Ratio is one order of magnitude higher for the larger motion amplitudes we tested, compared to the smaller motion amplitudes. Simulator-introduced noise is found to be primarily of deterministic nature, particularly for the stronger motion intensities. The effect of simulator noise on quantification of animal/human motion sensitivity is discussed. We conclude that accurate recording and characterization of executed simulator motion are a crucial prerequisite for the investigation of uncertainty in self-motion perception.     

popRange: a highly flexible spatially and temporally explicit Wright-Fisher simulator

Background

Sequencing and genotyping technology advancements have led to massive, growing repositories of spatially explicit genetic data and increasing quantities of temporal data (i.e., ancient DNA). These data will allow more complex and fine-scale inferences about population history than ever before; however, new methods are needed to test complex hypotheses.

Results

This article presents popRange, a forward genetic simulator, which incorporates large-scale genetic data with stochastic spatially and temporally explicit demographic and selective models. Features such as spatially and temporally variable selection coefficients and demography are incorporated in a highly flexible manner. popRange is implemented as an R package and presented with an example simulation exploring a selected allele’s trajectory in multiple subpopulations.

Conclusions

popRange allows researchers to evaluate and test complex scenarios by simulating large-scale data with complicated demographic and selective features. popRange is available for download at http://cran.r-project.org/web/packages/popRange/index.html.