In the yeast
Saccharomyces cerevisiae, the Ras/cAMP/PKA pathway plays a major role in the regulation of metabolism, stress resistance and cell cycle progression. We extend here a mechanistic model of ...the Ras/cAMP/PKA pathway that we previously defined by describing the molecular interactions and post-translational modifications of proteins, and perform a computational analysis to investigate the dynamical behaviors of the components of this pathway, regulated by different control mechanisms. We carry out stochastic simulations to consider, in particular, the effect of the negative feedback loops on the activity of both Ira2 (a Ras-GAP) and Cdc25 (a Ras-GEF) proteins. Our results show that stable oscillatory regimes for the dynamics of cAMP can be obtained only through the activation of these feedback mechanisms, and when the amount of Cdc25 is within a specific range. In addition, we highlight that the levels of guanine nucleotides pools are able to regulate the pathway, by influencing the transition between stable steady states and oscillatory regimes.
Reaction systems represent a theoretical framework based on the regulation mechanisms of facilitation and inhibition of biochemical reactions. The dynamic process defined by a reaction system is ...typically derived by hand, starting from the set of reactions and a given context sequence. However, this procedure may be error-prone and time-consuming, especially when the size of the reaction system increases. Here we present HERESY, a simulator of reaction systems accelerated on Graphics Processing Units (GPUs). HERESY is based on a fine-grained parallelization strategy, whereby all reactions are simultaneously executed on the GPU, therefore reducing the overall running time of the simulation. HERESY is particularly advantageous for the simulation of large-scale reaction systems, consisting of hundreds or thousands of reactions. By considering as test case some reaction systems with an increasing number of reactions and entities, as well as an increasing number of entities per reaction, we show that HERESY allows up to 29× speed-up with respect to a CPU-based simulator of reaction systems. Finally, we provide some directions for the optimization of HERESY, considering minimal reaction systems in normal form.
To investigate the behavior of biochemical systems, many runs of Gillespie’s Stochastic Simulation Algorithm (SSA) are generally needed, causing excessive computational costs on Central Processing ...Units (CPUs). Since all SSA runs are independent, the Intel Xeon Phi coprocessors based on the Many Integrated Core (MIC) architecture can be exploited to distribute the workload. We considered two execution modalities on MIC: one consisted in running exactly the same CPU code of SSA, while the other exploited MIC’s vector instructions to reuse the CPU code with only few modifications. MIC performance was compared with Graphics Processing Units (GPUs), specifically implemented in CUDA to optimize the use of memory hierarchy. Our results show that GPU largely outperforms MIC and CPU, but required a complete redesign of SSA. MIC allows a relevant speedup, especially when vector instructions are used, with the additional advantage of requiring minimal modifications to CPU code.
Metapopulations, or multi-patch systems, are models describing the interactions and the behavior of populations living in fragmented habitats. Dispersal, persistence and extinction are some of the ...characteristics of interest in ecological studies of metapopulations. In this paper, we propose a novel method to analyze metapopulations, which is based on a discrete and stochastic modelling framework in the area of Membrane Computing. New structural features of membrane systems, necessary to appropriately describe a multi-patch system, are introduced, such as the reduction of the maximal parallel consumption of objects, the spatial arrangement of membranes and the stochastic creation of objects. The role of the additional features, their meaning for a metapopulation model and the emergence of relevant behaviors are then investigated by means of stochastic simulations. Conclusive remarks and ideas for future research are finally presented.
The introduction of general-purpose Graphics Processing Units (GPUs) is boosting scientific applications in Bioinformatics, Systems Biology, and Computational Biology. In these fields, the use of ...high-performance computing solutions is motivated by the need of performing large numbers of in silico analysis to study the behavior of biological systems in different conditions, which necessitate a computing power that usually overtakes the capability of standard desktop computers. In this work we present coagSODA, a CUDA-powered computational tool that was purposely developed for the analysis of a large mechanistic model of the blood coagulation cascade (BCC), defined according to both mass-action kinetics and Hill functions. coagSODA allows the execution of parallel simulations of the dynamics of the BCC by automatically deriving the system of ordinary differential equations and then exploiting the numerical integration algorithm LSODA. We present the biological results achieved with a massive exploration of perturbed conditions of the BCC, carried out with one-dimensional and bi-dimensional parameter sweep analysis, and show that GPU-accelerated parallel simulations of this model can increase the computational performances up to a 181× speedup compared to the corresponding sequential simulations.
In the yeast
Saccharomyces cerevisiae
, the Ras/cAMP/PKA pathway is involved in the regulation of cell growth and proliferation in response to nutritional sensing and stress conditions. The pathway ...is tightly regulated by multiple feedback loops, exerted by the protein kinase A (PKA) on a few pivotal components of the pathway. In this article, we investigate the dynamics of the second messenger cAMP by performing stochastic simulations and parameter sweep analysis of a mechanistic model of the Ras/cAMP/PKA pathway, to determine the effects that the modulation of these feedback mechanisms has on the establishment of stable oscillatory regimes. In particular, we start by studying the role of phosphodiesterases, the enzymes that catalyze the degradation of cAMP, which represent the major negative feedback in this pathway. Then, we show the results on cAMP oscillations when perturbing the amount of protein Cdc25 coupled with the alteration of the intracellular ratio of the guanine nucleotides (GTP/GDP), which are known to regulate the switch of the GTPase Ras protein. This multi-level regulation of the amplitude and frequency of oscillations in the Ras/cAMP/PKA pathway might act as a fine tuning mechanism for the downstream targets of PKA, as also recently evidenced by some experimental investigations on the nucleocytoplasmic shuttling of the transcription factor Msn2 in yeast cells.
The modeling of biochemical reaction networks is a fundamental but complex task in Systems Biology, which is traditionally performed exploiting human expertise and the available experimental data. ...Because of the general lack of knowledge on the molecular mechanisms occurring in living cells, an intense research activity focused on the development of reverse engineering methodologies is currently underway. This problem is further complicated by the fact that a proper parameterization needs to be associated to the reaction network, in order to investigate its dynamical behavior. In this work we propose a novel computational methodology for the reverse engineering of fully parameterized kinetic networks, based on the combined use of two evolutionary programming techniques: Cartesian Genetic Programming (CGP) and Particle Swarm Optimization (PSO). In particular, CGP is used to infer the network topology, while PSO performs the parameter estimation task. To the purpose of applying our methodology in routine laboratory environments, we designed it to exploit a small set of experimental time series as target. We show that our methodology is able to reconstruct kinetic networks that perfectly fit with the target data.