Narrative-based computational modelling of the Gp130/JAK/STAT signalling pathway.

BACKGROUND: Appropriately formulated quantitative computational models can support researchers in understanding the dynamic behaviour of biological pathways and support hypothesis formulation and selection by "in silico" experimentation. An obstacle to widespread adoption of this approach is the requirement to formulate a biological pathway as machine executable computer code. We have recently proposed a novel, biologically intuitive, narrative-style modelling language for biologists to formulate the pathway which is then automatically translated into an executable format and is, thus, usable for analysis via existing simulation techniques. RESULTS: Here we use a high-level narrative language in designing a computational model of the gp130/JAK/STAT signalling pathway and show that the model reproduces the dynamic behaviour of the pathway derived by biological observation. We then "experiment" on the model by simulation and sensitivity analysis to define those parameters which dominate the dynamic behaviour of the pathway. The model predicts that nuclear compartmentalisation and phosphorylation status of STAT are key determinants of the pathway and that alternative mechanisms of signal attenuation exert their influence on different timescales. CONCLUSION: The described narrative model of the gp130/JAK/STAT pathway represents an interesting case study showing how, by using this approach, researchers can model biological systems without explicitly dealing with formal notations and mathematical expressions (typically used for biochemical modelling), nevertheless being able to obtain simulation and analysis results. We present the model and the sensitivity analysis results we have obtained, that allow us to identify the parameters which are most sensitive to perturbations. The results, which are shown to be in agreement with existing mathematical models of the gp130/JAK/STAT pathway, serve us as a form of validation of the model and of the approach itself

Snazer: the simulations and networks analyzer

<p>Abstract</p> <p>Background</p> <p>Networks are widely recognized as key determinants of structure and function in systems that span the biological, physical, and social sciences. They are static pictures of the interactions among the components of complex systems. Often, much effort is required to identify networks as part of particular patterns as well as to visualize and interpret them.</p> <p>From a pure dynamical perspective, simulation represents a relevant <it>way</it>-<it>out</it>. Many simulator tools capitalized on the "noisy" behavior of some systems and used formal models to represent cellular activities as temporal trajectories. Statistical methods have been applied to a fairly large number of replicated trajectories in order to infer knowledge.</p> <p>A tool which both graphically manipulates reactive models and deals with sets of simulation time-course data by aggregation, interpretation and statistical analysis is missing and could add value to simulators.</p> <p>Results</p> <p>We designed and implemented <it>Snazer</it>, the simulations and networks analyzer. Its goal is to aid the processes of visualizing and manipulating reactive models, as well as to share and interpret time-course data produced by stochastic simulators or by any other means.</p> <p>Conclusions</p> <p><it>Snazer </it>is a solid prototype that integrates biological network and simulation time-course data analysis techniques.</p

Fast network centrality analysis using GPUs

<p>Abstract</p> <p>Background</p> <p>With the exploding volume of data generated by continuously evolving high-throughput technologies, biological network analysis problems are growing larger in scale and craving for more computational power. General Purpose computation on Graphics Processing Units (GPGPU) provides a cost-effective technology for the study of large-scale biological networks. Designing algorithms that maximize data parallelism is the key in leveraging the power of GPUs.</p> <p>Results</p> <p>We proposed an efficient data parallel formulation of the All-Pairs Shortest Path problem, which is the key component for shortest path-based centrality computation. A betweenness centrality algorithm built upon this formulation was developed and benchmarked against the most recent GPU-based algorithm. Speedup between 11 to 19% was observed in various simulated scale-free networks. We further designed three algorithms based on this core component to compute closeness centrality, eccentricity centrality and stress centrality. To make all these algorithms available to the research community, we developed a software package <it>gpu</it>-<it>fan </it>(GPU-based Fast Analysis of Networks) for CUDA enabled GPUs. Speedup of 10-50× compared with CPU implementations was observed for simulated scale-free networks and real world biological networks.</p> <p>Conclusions</p> <p><it>gpu</it>-<it>fan </it>provides a significant performance improvement for centrality computation in large-scale networks. Source code is available under the GNU Public License (GPL) at <url>http://bioinfo.vanderbilt.edu/gpu-fan/</url>.</p

Simulating Food Web Dynamics along a Gradient: Quantifying Human Influence

Realistically parameterized and dynamically simulated food-webs are useful tool to explore the importance of the functional diversity of ecosystems, and in particular relations between the dynamics of species and the whole community. We present a stochastic dynamical food web simulation for the Kelian River (Borneo). The food web was constructed for six different locations, arrayed along a gradient of increasing human perturbation (mostly resulting from gold mining activities) along the river. Along the river, the relative importance of grazers, filterers and shredders decreases with increasing disturbance downstream, while predators become more dominant in governing eco-dynamics. Human activity led to increased turbidity and sedimentation which adversely impacts primary productivity. Since the main difference between the study sites was not the composition of the food webs (structure is quite similar) but the strengths of interactions and the abundance of the trophic groups, a dynamical simulation approach seemed to be useful to better explain human influence. In the pristine river (study site 1), when comparing a structural version of our model with the dynamical model we found that structurally central groups such as omnivores and carnivores were not the most important ones dynamically. Instead, primary consumers such as invertebrate grazers and shredders generated a greater dynamical response. Based on the dynamically most important groups, bottom-up control is replaced by the predominant top-down control regime as distance downstream and human disturbance increased. An important finding, potentially explaining the poor structure to dynamics relationship, is that indirect effects are at least as important as direct ones during the simulations. We suggest that our approach and this simulation framework could serve systems-based conservation efforts. Quantitative indicators on the relative importance of trophic groups and the mechanistic modeling of eco-dynamics could greatly contribute to understanding various aspects of functional diversity

Transcriptional Regulation Is a Major Controller of Cell Cycle Transition Dynamics

DNA replication, mitosis and mitotic exit are critical transitions of the cell cycle which normally occur only once per cycle. A universal control mechanism was proposed for the regulation of mitotic entry in which Cdk helps its own activation through two positive feedback loops. Recent discoveries in various organisms showed the importance of positive feedbacks in other transitions as well. Here we investigate if a universal control system with transcriptional regulation(s) and post-translational positive feedback(s) can be proposed for the regulation of all cell cycle transitions. Through computational modeling, we analyze the transition dynamics in all possible combinations of transcriptional and post-translational regulations. We find that some combinations lead to ‘sloppy’ transitions, while others give very precise control. The periodic transcriptional regulation through the activator or the inhibitor leads to radically different dynamics. Experimental evidence shows that in cell cycle transitions of organisms investigated for cell cycle dependent periodic transcription, only the inhibitor OR the activator is under cyclic control and never both of them. Based on these observations, we propose two transcriptional control modes of cell cycle regulation that either STOP or let the cycle GO in case of a transcriptional failure. We discuss the biological relevance of such differences

A framework for parameter estimation and model selection from experimental data in systems biology using approximate Bayesian computation.

As modeling becomes a more widespread practice in the life sciences and biomedical sciences, researchers need reliable tools to calibrate models against ever more complex and detailed data. Here we present an approximate Bayesian computation (ABC) framework and software environment, ABC-SysBio, which is a Python package that runs on Linux and Mac OS X systems and that enables parameter estimation and model selection in the Bayesian formalism by using sequential Monte Carlo (SMC) approaches. We outline the underlying rationale, discuss the computational and practical issues and provide detailed guidance as to how the important tasks of parameter inference and model selection can be performed in practice. Unlike other available packages, ABC-SysBio is highly suited for investigating, in particular, the challenging problem of fitting stochastic models to data. In order to demonstrate the use of ABC-SysBio, in this protocol we postulate the existence of an imaginary reaction network composed of seven interrelated biological reactions (involving a specific mRNA, the protein it encodes and a post-translationally modified version of the protein), a network that is defined by two files containing 'observed' data that we provide as supplementary information. In the first part of the PROCEDURE, ABC-SysBio is used to infer the parameters of this system, whereas in the second part we use ABC-SysBio's relevant functionality to discriminate between two different reaction network models, one of them being the 'true' one. Although computationally expensive, the additional insights gained in the Bayesian formalism more than make up for this cost, especially in complex problems

Subliminally Perceived Odours Modulate Female Intrasexual Competition: An Eye Movement Study

Background: Evidence suggests that subliminal odorants influence human perception and behavior. It has been hypothesized that the human sex-steroid derived compound 4,16-androstadien-3-one (androstadienone) functions as a human chemosignal. The most intensively studied steroid compound, androstadienone is known to be biologically relevant since it seems to convey information about male mate quality to women. It is unclear if the effects of androstadienone are menstrual cycle related. Methodology/Principal Findings: In the first experiment, heterosexual women were exposed to androstadienone or a control compound and asked to view stimuli such as female faces, male faces and familiar objects while their eye movements were recorded. In the second experiment the same women were asked to rate the level of stimuli attractiveness following exposure to the study or control compound. The results indicated that women at high conception risk spent more time viewing the female than the male faces regardless of the compound administered. Women at a low conception risk exhibited a preference for female faces only following exposure to androstadienone. Conclusions/Significance: We contend that a woman’s level of fertility influences her evaluation of potential competitors (e.g., faces of other women) during times critical for reproduction. Subliminally perceived odorants, such as androstadienone, might similarly enhance intrasexual competition strategies in women during fertility phases not critica

A High-Level Language for Rule-Based Modelling

Rule-based languages such as Kappa excel in their support for handling the combinatorial complexities prevalent in many biological systems, including signalling pathways. But Kappa provides little structure for organising rules, and large models can therefore be hard to read and maintain. This paper introduces a high-level, modular extension of Kappa called LBS-κ. We demonstrate the constructs of the language through examples and three case studies: a chemotaxis switch ring, a MAPK cascade, and an insulin signalling pathway. We then provide a formal definition of LBS-κ through an abstract syntax and a translation to plain Kappa. The translation is implemented in a compiler tool which is available as a web application. We finally demonstrate how to increase the expressivity of LBS-κ through embedded scripts in a general-purpose programming language, a technique which we view as generally applicable to other domain specific languages