Modeling Of Biological Regulatory Networks Research Articles

Revision of Boolean Models of Regulatory Networks Using Stable State Observations.

Models of biological regulatory networks are essential to understand cellular processes. However, the definition of such models is still mostly manually performed, and consequently prone to error. Moreover, as new experimental data are acquired, models need to be revised and updated. Here, we propose a model revision procedure and associated tool, capable of providing the set of minimal repairs to render a model consistent with a set of experimental observations. We consider four possible repair operations, using a lexicographic optimization criterion, giving preference to function repairs over topological ones. Also, we consider observations at stable state discarding the model dynamics. In this article, we extend our previous work to tackle the problem of repairing nodes with multiple reasons of inconsistency. We evaluate our tool on five publicly available logical models. We perform random changes considering several parameter configurations to assess the tool repairing capabilities. Whenever a model is repaired under the time limit, the tool successfully produces the optimal solutions to repair the model. Instances were generated without the previous limitation to validate this extended approach.

BioLQM: A Java Toolkit for the Manipulation and Conversion of Logical Qualitative Models of Biological Networks.

Here we introduce bioLQM, a new Java software toolkit for the conversion, modification, and analysis of Logical Qualitative Models of biological regulatory networks. BioLQM provides core modeling operations as building blocks for the development of integrated modeling software, or for the assembly of heterogeneous analysis workflows involving several complementary tools. Based on the definition of multi-valued logical models, bioLQM implements import and export facilities, notably for the recent SBML qual exchange format, as well as for formats used by several popular tools, facilitating the design of workflows combining these tools. Model modifications enable the definition of various perturbations, as well as model reduction, easing the analysis of large models. Another modification enables the study of multi-valued models with tools limited to the Boolean case. Finally, bioLQM provides a framework for the development of novel analysis tools. The current version implements various updating modes for model simulation (notably synchronous, asynchronous, and random asynchronous), as well as some static analysis features for the identification of attractors. The bioLQM software can be integrated into analysis workflows through command line and scripting interfaces. As a Java library, it further provides core data structures to the GINsim and EpiLog interactive tools, which supply graphical interfaces and additional analysis methods for cellular and multi-cellular qualitative models.

Identification of bifurcation transitions in biological regulatory networks using Answer-Set Programming

BackgroundNumerous cellular differentiation processes can be captured using discrete qualitative models of biological regulatory networks. These models describe the temporal evolution of the state of the network subject to different competing transitions, potentially leading the system to different attractors. This paper focusses on the formal identification of states and transitions that are crucial for preserving or pre-empting the reachability of a given behaviour.MethodsIn the context of non-deterministic automata networks, we propose a static identification of so-called bifurcations, i.e., transitions after which a given goal is no longer reachable. Such transitions are naturally good candidates for controlling the occurrence of the goal, notably by modulating their propensity. Our method combines Answer-Set Programming with static analysis of reachability properties to provide an under-approximation of all the existing bifurcations.ResultsWe illustrate our discrete bifurcation analysis on several models of biological systems, for which we identify transitions which impact the reachability of given long-term behaviour. In particular, we apply our implementation on a regulatory network among hundreds of biological species, supporting the scalability of our approach.ConclusionsOur method allows a formal and scalable identification of transitions which are responsible for the lost of capability to reach a given state. It can be applied to any asynchronous automata networks, which encompass Boolean and multi-valued models. An implementation is provided as part of the Pint software, available at http://loicpauleve.name/pint.

Learning Delays in Biological Regulatory Networks from Time Series Data

Models of Biological Regulatory Networks are generally based on prior knowledge, either derived from literature and/or the manual analysis of biological observations. With the development of high-throughput data, there is a growing need for methods that automatically generate admissible models. To have a better understanding of the dynamical phenomena at stake in the influences between biological components, it would be necessary to include delayed influences in the model. The main purpose of this work is to have a resulting network as consistent as possible with the observed datasets regarding the conflicts and the simultaneity between transitions. The originality of our work is threefold: (i) the identification the sign of the interactions, (ii) the direct integration of quantitative time delays in the learning approach and (iii) the identification of the qualitative discrete levels that lead to the systems dynamics.In this work the precision of our automatic approach is discussed by applying it on dynamical biological models coming from the DREAM4 Challenge datasets.

Modeling Delayed Dynamics in Biological Regulatory Networks from Time Series Data

Background: The modeling of Biological Regulatory Networks (BRNs) relies on background knowledge, deriving either from literature and/or the analysis of biological observations. However, with the development of high-throughput data, there is a growing need for methods that automatically generate admissible models. Methods: Our research aim is to provide a logical approach to infer BRNs based on given time series data and known influences among genes. Results: We propose a new methodology for models expressed through a timed extension of the automata networks (well suited for biological systems). The main purpose is to have a resulting network as consistent as possible with the observed datasets. Conclusion: The originality of our work is three-fold: (i) identifying the sign of the interaction; (ii) the direct integration of quantitative time delays in the learning approach; and (iii) the identification of the qualitative discrete levels that lead to the systems’ dynamics. We show the benefits of such an automatic approach on dynamical biological models, the DREAM4(in silico) and DREAM8 (breast cancer) datasets, popular reverse-engineering challenges, in order to discuss the precision and the computational performances of our modeling method.

A comparative study of qualitative and quantitative dynamic models of biological regulatory networks

Mathematical modeling of biological regulatory networks provides valuable insights into the structural and dynamical properties of the underlying systems. While dynamic models based on differential equations provide quantitative information on the biological systems, qualitative models that rely on the logical interactions among the components provide coarse-grained descriptions useful for systems whose mechanistic underpinnings remain incompletely understood. The middle ground class of piecewise affine differential equation models was proven informative for systems with partial knowledge of kinetic parameters. In this work we provide a comparison of the dynamic characteristics of these three approaches applied on several biological regulatory network motifs. Specifically, we compare the attractors and state transitions in asynchronous Boolean, piecewise affine and Hill-type continuous models. Our study shows that while the fixed points of asynchronous Boolean models are observed in continuous Hill-type and piecewise affine models, these models may exhibit different attractors under certain conditions. Overall, qualitative models are suitable for systems with limited knowledge of quantitative information. On the other hand, when practical, using quantitative models can provide detailed information about additional real-valued attractors not present in the qualitative models.

Parametric linear hybrid automata for complex environmental systems modeling

Environmental systems, whether they be weather patterns or predator-prey relationships, are dependent on a number of different variables, each directly or indirectly affecting the system at large. Since not all of these factors are known, these systems take on non-linear dynamics, making it difficult to accurately predict meaningful behavioral trends far into the future. However, such dynamics do not warrant complete ignorance of different efforts to understand and model close approximations of these systems. Towards this end, we have applied a logical modeling approach to model and analyze the behavioral trends and systematic trajectories that these systems exhibit without delving into their quantification. This approach, formalized by Rene Thomas for discrete logical modeling of Biological Regulatory Networks (BRNs) and further extended in our previous studies as parametric biological linear hybrid automata (Bio-LHA), has been previously employed for the analyses of different molecular regulatory interactions occurring across various cells and microbial species. As relationships between different interacting components of a system can be simplified as positive or negative influences, we can employ the Bio-LHA framework to represent different components of the environmental system as positive or negative feedbacks. In the present study, we highlight the benefits of hybrid (discrete/continuous) modeling which lead to refinements among the fore-casted behaviors in order to find out which ones are actually possible. We have taken two case studies: an interaction of three microbial species in a freshwater pond, and a more complex atmospheric system, to show the applications of the Bio-LHA methodology for the timed hybrid modeling of environmental systems. Results show that the approach using the Bio-LHA is a viable method for behavioral modeling of complex environmental systems by finding timing constraints while keeping the complexity of the model at a minimum.

Identification of biological regulatory networks from Process Hitting models

Qualitative formalisms offer a well-established alternative to the more traditionally used differential equation models of Biological Regulatory Networks (BRNs). These formalisms led to numerous theoretical works and practical tools to understand emerging behaviors. The analysis of the dynamics of very large models is however a rather hard problem, which led us to previously introduce the Process Hitting framework (PH), which is a particular class of nondeterministic asynchronous automata network (or safe Petri nets). Its major advantage lies in the efficiency of several static analyses recently designed to assess dynamical properties, making it possible to tackle very large models.In this paper, we address the formal identification of qualitative models of BRNs from PH models. First, the inference of the Interaction Graph from a PH model summarizes the signed influences between the components that are effective for the dynamics. Second, we provide the inference of all René Thomas models of BRNs that are compatible with a given PH. As the PH allows the specification of nondeterministic interactions between components, our inference emphasizes the ability of PH to deal with large BRNs with incomplete knowledge on interactions, where Thomas' approach fails because of the combinatorics of parameters.The inference of corresponding Thomas models is implemented using Answer Set Programming, which allows in particular an efficient enumeration of (possibly numerous) compatible parameterizations.

A REDUCTION METHOD FOR BOOLEAN NETWORK MODELS PROVEN TO CONSERVE ATTRACTORS.

Boolean models, wherein each component is characterized with a binary (ON or OFF) variable, have been widely employed for dynamic modeling of biological regulatory networks. However, the exponential dependencse of the size of the state space of these models on the number of nodes in the network can be a daunting prospect for attractor analysis of large-scale systems. We have previously proposed a network reduction technique for Boolean models and demonstrated its applicability on two biological systems, namely, the abscisic acid signal transduction network as well as the T-LGL leukemia survival signaling network. In this paper, we provide a rigorous mathematical proof that this method not only conserves the fixed points of a Boolean network, but also conserves the complex attractors of general asynchronous Boolean models wherein at each time step a randomly selected node is updated. This method thus allows one to infer the long-term dynamic properties of a large-scale system from those of the corresponding reduced model.

Boolean modeling of biological regulatory networks: A methodology tutorial

Given the complexity and interactive nature of biological systems, constructing informative and coherent network models of these systems and subsequently developing efficient approaches to analyze the assembled networks is of immense importance. The integration of network analysis and dynamic modeling enables one to investigate the behavior of the underlying system as a whole and to make experimentally testable predictions about less-understood aspects of the processes involved. In this paper, we present a tutorial on the fundamental steps of Boolean modeling of biological regulatory networks. We demonstrate how to infer a Boolean network model from the available experimental data, analyze the network using graph-theoretical measures, and convert it into a predictive dynamic model. For each step, the pitfalls one may encounter and possible ways to circumvent them are also discussed. We illustrate these steps on a toy network as well as in the context of the Drosophila melanogaster segment polarity gene network.

Analysis and Characterization of Asynchronous State Transition Graphs Using Extremal States

Logical modeling of biological regulatory networks gives rise to a representation of the system's dynamics as a so-called state transition graph. Analysis of such a graph in its entirety allows for a comprehensive understanding of the functionalities and behavior of the modeled system. However, the size of the vertex set of the graph is exponential in the number of the network components making analysis costly, motivating development of reduction methods. In this paper, we present results allowing for a complete description of an asynchronous state transition graph of a Thomas network solely based on the analysis of the subgraph induced by certain extremal states. Utilizing this notion, we compare the behavior of a simple multivalued network and a corresponding Boolean network and analyze the conservation of dynamical properties between them. Understanding the relation between such coarser and finer models is a necessary step toward meaningful network reduction as well as model refinement methods.