Qualitative modelling and analysis of regulations in multi-cellular systems using Petri nets and topological collections

In this paper, we aim at modelling and analyzing the regulation processes in multi-cellular biological systems, in particular tissues. The modelling framework is based on interconnected logical regulatory networks a la Rene Thomas equipped with information about their spatial relationships. The semantics of such models is expressed through colored Petri nets to implement regulation rules, combined with topological collections to implement the spatial information. Some constraints are put on the the representation of spatial information in order to preserve the possibility of an enumerative and exhaustive state space exploration. This paper presents the modelling framework, its semantics, as well as a prototype implementation that allowed preliminary experimentation on some applications.Comment: In Proceedings MeCBIC 2010, arXiv:1011.005

Bringing Dicynodonts Back to Life: Paleobiology and Anatomy of a New Emydopoid Genus from the Upper Permian of Mozambique

Dicynodontia represent the most diverse tetrapod group during the Late Permian. They survived the Permo-Triassic extinction and are central to understanding Permo-Triassic terrestrial ecosystems. Although extensively studied, several aspects of dicynodont paleobiology such as, neuroanatomy, inner ear morphology and internal cranial anatomy remain obscure. Here we describe a new dicynodont (Therapsida, Anomodontia) from northern Mozambique: Niassodon mfumukasi gen. et sp. nov. The holotype ML1620 was collected from the Late Permian K5 formation, Metangula Graben, Niassa Province northern Mozambique, an almost completely unexplored basin and country for vertebrate paleontology. Synchrotron radiation based micro-computed tomography (SRµCT), combined with a phylogenetic analysis, demonstrates a set of characters shared with Emydopoidea. All individual bones were digitally segmented allowing a 3D visualization of each element. In addition, we reconstructed the osseous labyrinth, endocast, cranial nerves and vasculature. The brain is narrow and the cerebellum is broader than the forebrain, resembling the conservative, "reptilian-grade" morphology of other non-mammalian therapsids, but the enlarged paraflocculi occupy the same relative volume as in birds. The orientation of the horizontal semicircular canals indicates a slightly more dorsally tilted head posture than previously assumed in other dicynodonts. In addition, synchrotron data shows a secondary center of ossification in the femur. Thus ML1620 represents, to our knowledge, the oldest fossil evidence of a secondary center of ossification, pushing back the evolutionary origins of this feature. The fact that the specimen represents a new species indicates that the Late Permian tetrapod fauna of east Africa is still incompletely known.Mozambique (Ministério dos Recursos Minerais), National Geographic Society, TAP airlines and other anonymous patrons, financial support from DESY through the I-20110184 EC project

Dynamic modeling and simulation of leukocyte integrin activation through an electronic design automation framework

Model development and analysis of biological systems is recognized as a key requirement for integrating in-vitro and in-vivo experimental data. In-silico simulations of a biochemical model allows one to test different experimental conditions, helping in the discovery of the dynamics that regulate the system. Several characteristics and issues of biological system modeling are common to the electronics system modeling, such as concurrency, reactivity, abstraction levels, as well as state space explosion during verification. This paper proposes a modeling and simulation framework for discrete event-based execution of biochemical systems based on SystemC. SystemC is the reference language in the electronic design automation (EDA) field for modeling and verifying complex systems at different abstraction levels. SystemC-based verification is the de-facto an alternative to model checking when such a formal verification technique cannot deal with the state space complexity of the model. The paper presents how the framework has been applied to model the intracellular signalling network controlling integrin activation mediating leukocyte recruitment from the blood into the tissues, by handling the solution space complexity through different levels of simulation accuracy

Characterization of Reachable Attractors Using Petri Net Unfoldings

International audienceAttractors of network dynamics represent the long-term behaviours of the modelled system. Their characterization is therefore crucial for understanding the response and differentiation capabilities of a dynamical system. In the scope of qualitative models of interaction networks, the computation of attractors reachable from a given state of the network faces combinatorial issues due to the state space explosion. In this paper, we present a new algorithm that exploits the concurrency between transitions of parallel acting components in order to reduce the search space. The algorithm relies on Petri net unfoldings that can be used to compute a compact representation of the dynamics. We illustrate the applicability of the algorithm with Petri net models of cell signalling and regulation networks, Boolean and multi-valued. The proposed approach aims at being complementary to existing methods for deriving the attractors of Boolean models, while being %so far more generic since it applies to any safe Petri net

Systems Biology Markup Language (SBML): Language Specification for Level 3 Version 2 Core Release 2

Computational models can help researchers to interpret data, understand biological functions, and make quantitative predictions. The Systems Biology Markup Language (SBML) is a file format for representing computational models in a declarative form that different software systems can exchange. SBML is oriented towards describing biological processes of the sort common in research on a number of topics, including metabolic pathways, cell signaling pathways, and many others. By supporting SBML as an input/output format, different tools can all operate on an identical representation of a model, removing opportunities for translation errors and assuring a common starting point for analyses and simulations. This document provides the specification for Release 2 of Version 2 of SBML Level 3 Core. The specification defines the data structures prescribed by SBML as well as their encoding in XML, the eXtensible Markup Language. Release 2 corrects some errors and clarifies some ambiguities discovered in Release 1. This specification also defines validation rules that determine the validity of an SBML document, and provides many examples of models in SBML form. Other materials and software are available from the SBML project website at http://sbml.org/

A critical review on modelling formalisms and simulation tools in computational biosystems

Integration of different kinds of biological processes is an ultimate goal for whole-cell modelling. We briefly review modelling formalisms that have been used in Systems Biology and identify the criteria that must be addressed by an integrating framework capable of modelling, analysing and simulating different biological networks. Aware that no formalism can fit all purposes we realize Petri nets as a suitable model for Metabolic Engineering and take a deeper perspective on the role of this formalism as an integrating framework for regulatory and metabolic networks.Research supported by PhD grant SFRH/BD/35215/2007 from the Fundacao para a Ciencia e a Tecnologia (FCT) and the MIT-Portugal program

Random sampling of elementary flux modes in large-scale metabolic networks

Motivation: The description of a metabolic network in terms of elementary (flux) modes (EMs) provides an important framework for metabolic pathway analysis. However, their application to large networks has been hampered by the combinatorial explosion in the number of modes. In this work, we develop a method for generating random samples of EMs without computing the whole set. Results: Our algorithm is an adaptation of the canonical basis approach, where we add an additional filtering step which, at each iteration, selects a random subset of the new combinations of modes. In order to obtain an unbiased sample, all candidates are assigned the same probability of getting selected. This approach avoids the exponential growth of the number of modes during computation, thus generating a random sample of the complete set of EMs within reasonable time. We generated samples of different sizes for a metabolic network of Escherichia coli, and observed that they preserve several properties of the full EM set. It is also shown that EM sampling can be used for rational strain design. A well distributed sample, that is representative of the complete set of EMs, should be suitable to most EM-based methods for analysis and optimization of metabolic networks

Systematic reconstruction of TRANSPATH data into Cell System Markup Language

<p>Abstract</p> <p>Background</p> <p>Many biological repositories store information based on experimental study of the biological processes within a cell, such as protein-protein interactions, metabolic pathways, signal transduction pathways, or regulations of transcription factors and miRNA. Unfortunately, it is difficult to directly use such information when generating simulation-based models. Thus, modeling rules for encoding biological knowledge into system-dynamics-oriented standardized formats would be very useful for fully understanding cellular dynamics at the system level.</p> <p>Results</p> <p>We selected the TRANSPATH database, a manually curated high-quality pathway database, which provides a plentiful source of cellular events in humans, mice, and rats, collected from over 31,500 publications. In this work, we have developed 16 modeling rules based on hybrid functional Petri net with extension (HFPNe), which is suitable for graphical representing and simulating biological processes. In the modeling rules, each Petri net element is incorporated with Cell System Ontology to enable semantic interoperability of models. As a formal ontology for biological pathway modeling with dynamics, CSO also defines biological terminology and corresponding icons. By combining HFPNe with the CSO features, it is possible to make TRANSPATH data to simulation-based and semantically valid models. The results are encoded into a biological pathway format, Cell System Markup Language (CSML), which eases the exchange and integration of biological data and models.</p> <p>Conclusion</p> <p>By using the 16 modeling rules, 97% of the reactions in TRANSPATH are converted into simulation-based models represented in CSML. This reconstruction demonstrates that it is possible to use our rules to generate quantitative models from static pathway descriptions.</p