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

SimService: a lightweight library for building simulation services in Python

SUMMARY

Integrative biological modeling requires software infrastructure to launch, interconnect, and execute simulation software components without loss of functionality. SimService is a software library that enables deploying simulations in integrated applications as memory-isolated services with interactive proxy objects in the Python programming language. SimService supports customizing the interface of proxies so that simulation developers and users alike can tailor generated simulation instances according to model, method, and integrated application.

AVAILABILITY AND IMPLEMENTATION

SimService is written in Python, is freely available on GitHub under the MIT license at https://github.com/tjsego/simservice, and is available for download via the Python Package Index (package name "simservice") and conda (package name "simservice" on the conda-forge channel).

Morphogenetic systems: Models and experiments

M systems are mathematical models of morphogenesis developed to gain insights into its relations to phenomena such as self-assembly, self-controlled growth, homeostasis, self-healing and self-reproduction, in both natural and artificial systems. M systems rely on basic principles of membrane computing and self-assembly, as well as explicit emphasis on geometrical structures (location and shape) in 2D, 3D or higher dimensional Euclidean spaces. They can be used for principled studies of these phenomena, both theoretically and experimentally, at a computational level abstracted from their detailed implementation. In particular, they afford 2D and 3D models to explore biological morphogenetic processes. Theoretical studies have shown that M systems are powerful tools (e.g., computational universal, i.e. can become as complex as any computer program) and their parallelism allows for trading space for time in solving efficiently problems considered infeasible on conventional computers (NP-hard problems). In addition, they can also exhibit properties such as robustness to injuries and degrees of self-healing. This paper focuses on the experimental side of M systems. To this end, we have developed a high-level morphogenetic simulator, Cytos, to implement and visualize M systems in silico in order to verify theoretical results and facilitate research in M systems. We summarize the software package and make a brief comparison with some other simulators of membrane systems. The core of the article is a description of a range of experiments inspired by aspects of morphogenesis in both prokaryotic and eukaryotic cells. The experiments explore the regulatory role of the septum and of the cytoskeleton in cell fission, the robustness of cell models against injuries, and, finally, the impact of changing nutrient concentration on population growth.

Inflammatory Bowel Disease Meets Systems Biology: A Multi-Omics Challenge and Frontier

The inflammatory bowel disease (IBD) is a systemic disease that is characterized by the inflammation of the gastrointestinal tract. It includes ulcerative colitis and the Crohn's disease. Presently, IBD is one of the most investigated common complex human disorders, although its causes remain unclear. Multi-omics mechanisms involving genomic, transcriptomic, proteomic, and epigenomic variations, not to forget the miRNome, together with environmental contributions, result in an impairment of the immune system in persons with IBD. Such interactions at multiple levels of biology and in concert with the environment constitute the actual engine of this complex disease, demanding a multifactorial and multi-omics perspective to better understand the root causes of IBD. This expert analysis reviews and examines the latest literature and underscores, from the perspective of systems biology, the value of multi-omics technologies as opportunities to unravel the "IBD integrome." We anticipate that multi-omics research will accelerate the new discoveries and insights on IBD in the near future. It shall also pave the way for early diagnosis and help clinicians and families with IBD to forecast and make informed decisions about the prognosis and, possibly, personalized therapeutics in the future.

Fasting cycles potentiate the efficacy of gemcitabine treatment in in vitro and in vivo pancreatic cancer models

Background/aims

Pancreatic cancer (PC) is ranked as the fourth leading cause of cancer-related deaths worldwide. Despite recent advances in treatment options, a modest impact on the outcome of the disease is observed so far. Short-term fasting cycles have been shown to potentiate the efficacy of chemotherapy against glioma. The aim of this study was to assess the effect of fasting cycles on the efficacy of gemcitabine, a standard treatment for PC patients, in vitro and in an in vivo pancreatic cancer mouse xenograft model.

Materials and Methods

BxPC-3, MiaPaca-2 and Panc-1 cells were cultured in standard and fasting mimicking culturing condition to evaluate the effects of gemcitabine. Pancreatic cancer xenograft mice were subjected to 24h starvation prior to gemcitabine injection to assess the tumor volume and weight as compared to mice fed ad libitum.

Results

Fasted pancreatic cancer cells showed increased levels of equilibrative nucleoside transporter (hENT1), the transporter of gemcitabine across the cell membrane, and decreased ribonucleotide reductase M1 (RRM1) levels as compared to those cultured in standard medium. Gemcitabine was more effective in inducing cell death on fasted cells as compared to controls. Consistently, xenograft pancreatic cancer mice subjected to fasting cycles prior to gemcitabine injection displayed a decrease of more than 40% in tumor growth.

Conclusion

Fasting cycles enhance gemcitabine effect in vitro and in the in vivo PC xenograft mouse model. These results suggest that restrictive dietary interventions could enhance the efficacy of existing cancer treatments in pancreatic cancer patients.

Stochastic and deterministic multiscale models for systems biology: an auxin-transport case study

BACKGROUND

Stochastic and asymptotic methods are powerful tools in developing multiscale systems biology models; however, little has been done in this context to compare the efficacy of these methods. The majority of current systems biology modelling research, including that of auxin transport, uses numerical simulations to study the behaviour of large systems of deterministic ordinary differential equations, with little consideration of alternative modelling frameworks.

RESULTS

In this case study, we solve an auxin-transport model using analytical methods, deterministic numerical simulations and stochastic numerical simulations. Although the three approaches in general predict the same behaviour, the approaches provide different information that we use to gain distinct insights into the modelled biological system. We show in particular that the analytical approach readily provides straightforward mathematical expressions for the concentrations and transport speeds, while the stochastic simulations naturally provide information on the variability of the system.

CONCLUSIONS

Our study provides a constructive comparison which highlights the advantages and disadvantages of each of the considered modelling approaches. This will prove helpful to researchers when weighing up which modelling approach to select. In addition, the paper goes some way to bridging the gap between these approaches, which in the future we hope will lead to integrative hybrid models.

Snazer: the simulations and networks analyzer

BACKGROUND

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.From a pure dynamical perspective, simulation represents a relevant way-out. 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.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.

RESULTS

We designed and implemented Snazer, 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.

CONCLUSIONS

Snazer is a solid prototype that integrates biological network and simulation time-course data analysis techniques.

Rate and extent of protein localization is controlled by peptide-binding domain association kinetics and morphology

Protein localization is an important regulatory mechanism in many cell signaling pathways such as cytoskeletal organization and genetic regulation. The specific mechanism of protein localization determines the kinetics and morphological constraints of protein translocation, and thus affects the rate and extent of localization. To investigate the affect of localization kinetics and morphology on protein localization, we designed a protein localization system based on Ca(2+)-calmodulin and Src homology 3 domain binding peptides that can translocate between specific localizations in response to a Ca(2+) signal. We used a stochastic biomolecular simulator to predict that such a protein localization system will exhibit slower and less complete translocations when the association kinetics of a binding domain and peptide are reduced. As well, we predicted that increasing the diffusion resistance by manipulating the morphology of the system would similarly impair translocation speed and completeness. We then constructed a network of synthetic fusion proteins and showed that these predictions could be qualitatively confirmed in vitro. This work provides a basis for explaining the different characteristics (rate and extent) of protein transport and localization in cells as a consequence of the kinetics and morphology of the transport mechanism.

Cascading signaling pathways improve the fidelity of a stochastically and deterministically simulated molecular RS latch

Background

While biological systems have often been compared with digital systems, they differ by the strong effect of crosstalk between signals due to diffusivity in the medium, reaction kinetics and geometry. Memory elements have allowed the creation of autonomous digital systems and although biological systems have similar properties of autonomy, equivalent memory mechanisms remain elusive. Any such equivalent memory system, however, must silence the effect of crosstalk to maintain memory fidelity.

Results

Here, we present a system of enzymatic reactions that behaves like an RS latch (a simple memory element in digital systems). Using both a stochastic molecular simulator and ordinary differential equation simulator, we showed that crosstalk between two latches operating in the same spatial localization disrupts the memory fidelity of both latches. Crosstalk was reduced or silenced when simple reaction loops were replaced with multiple step or cascading reactions, showing that cascading signaling pathways are less susceptible to crosstalk.

Conclusion

Thus, the common biological theme of cascading signaling pathways is advantageous for maintaining the fidelity of a memory latch in the presence of crosstalk. The experimental implementation of such a latch system will lead to novel approaches to cell control using synthetic proteins and will contribute to our understanding of why cells behave differently even when given the same stimulus.

ALC: automated reduction of rule-based models

Background

Combinatorial complexity is a challenging problem for the modeling of cellular signal transduction since the association of a few proteins can give rise to an enormous amount of feasible protein complexes. The layer-based approach is an approximative, but accurate method for the mathematical modeling of signaling systems with inherent combinatorial complexity. The number of variables in the simulation equations is highly reduced and the resulting dynamic models show a pronounced modularity. Layer-based modeling allows for the modeling of systems not accessible previously.

Results

ALC (Automated Layer Construction) is a computer program that highly simplifies the building of reduced modular models, according to the layer-based approach. The model is defined using a simple but powerful rule-based syntax that supports the concepts of modularity and macrostates. ALC performs consistency checks on the model definition and provides the model output in different formats (C MEX, MATLAB, Mathematica and SBML) as ready-to-run simulation files. ALC also provides additional documentation files that simplify the publication or presentation of the models. The tool can be used offline or via a form on the ALC website.

Conclusion

ALC allows for a simple rule-based generation of layer-based reduced models. The model files are given in different formats as ready-to-run simulation files.

A structured approach for the engineering of biochemical network models, illustrated for signalling pathways

Quantitative models of biochemical networks (signal transduction cascades, metabolic pathways, gene regulatory circuits) are a central component of modern systems biology. Building and managing these complex models is a major challenge that can benefit from the application of formal methods adopted from theoretical computing science. Here we provide a general introduction to the field of formal modelling, which emphasizes the intuitive biochemical basis of the modelling process, but is also accessible for an audience with a background in computing science and/or model engineering. We show how signal transduction cascades can be modelled in a modular fashion, using both a qualitative approach--qualitative Petri nets, and quantitative approaches--continuous Petri nets and ordinary differential equations (ODEs). We review the major elementary building blocks of a cellular signalling model, discuss which critical design decisions have to be made during model building, and present a number of novel computational tools that can help to explore alternative modular models in an easy and intuitive manner. These tools, which are based on Petri net theory, offer convenient ways of composing hierarchical ODE models, and permit a qualitative analysis of their behaviour. We illustrate the central concepts using signal transduction as our main example. The ultimate aim is to introduce a general approach that provides the foundations for a structured formal engineering of large-scale models of biochemical networks.