16,133 research outputs found
Reduction of dynamical biochemical reaction networks in computational biology
Biochemical networks are used in computational biology, to model the static
and dynamical details of systems involved in cell signaling, metabolism, and
regulation of gene expression. Parametric and structural uncertainty, as well
as combinatorial explosion are strong obstacles against analyzing the dynamics
of large models of this type. Multi-scaleness is another property of these
networks, that can be used to get past some of these obstacles. Networks with
many well separated time scales, can be reduced to simpler networks, in a way
that depends only on the orders of magnitude and not on the exact values of the
kinetic parameters. The main idea used for such robust simplifications of
networks is the concept of dominance among model elements, allowing
hierarchical organization of these elements according to their effects on the
network dynamics. This concept finds a natural formulation in tropical
geometry. We revisit, in the light of these new ideas, the main approaches to
model reduction of reaction networks, such as quasi-steady state and
quasi-equilibrium approximations, and provide practical recipes for model
reduction of linear and nonlinear networks. We also discuss the application of
model reduction to backward pruning machine learning techniques
A Taxonomy of Causality-Based Biological Properties
We formally characterize a set of causality-based properties of metabolic
networks. This set of properties aims at making precise several notions on the
production of metabolites, which are familiar in the biologists' terminology.
From a theoretical point of view, biochemical reactions are abstractly
represented as causal implications and the produced metabolites as causal
consequences of the implication representing the corresponding reaction. The
fact that a reactant is produced is represented by means of the chain of
reactions that have made it exist. Such representation abstracts away from
quantities, stoichiometric and thermodynamic parameters and constitutes the
basis for the characterization of our properties. Moreover, we propose an
effective method for verifying our properties based on an abstract model of
system dynamics. This consists of a new abstract semantics for the system seen
as a concurrent network and expressed using the Chemical Ground Form calculus.
We illustrate an application of this framework to a portion of a real
metabolic pathway
An Introduction to Rule-based Modeling of Immune Receptor Signaling
Cells process external and internal signals through chemical interactions.
Cells that constitute the immune system (e.g., antigen presenting cell, T-cell,
B-cell, mast cell) can have different functions (e.g., adaptive memory,
inflammatory response) depending on the type and number of receptor molecules
on the cell surface and the specific intracellular signaling pathways activated
by those receptors. Explicitly modeling and simulating kinetic interactions
between molecules allows us to pose questions about the dynamics of a signaling
network under various conditions. However, the application of chemical kinetics
to biochemical signaling systems has been limited by the complexity of the
systems under consideration. Rule-based modeling (BioNetGen, Kappa, Simmune,
PySB) is an approach to address this complexity. In this chapter, by
application to the FcRI receptor system, we will explore the
origins of complexity in macromolecular interactions, show how rule-based
modeling can be used to address complexity, and demonstrate how to build a
model in the BioNetGen framework. Open source BioNetGen software and
documentation are available at http://bionetgen.org.Comment: 5 figure
Model reduction for stochastic CaMKII reaction kinetics in synapses by graph-constrained correlation dynamics.
A stochastic reaction network model of Ca(2+) dynamics in synapses (Pepke et al PLoS Comput. Biol. 6 e1000675) is expressed and simulated using rule-based reaction modeling notation in dynamical grammars and in MCell. The model tracks the response of calmodulin and CaMKII to calcium influx in synapses. Data from numerically intensive simulations is used to train a reduced model that, out of sample, correctly predicts the evolution of interaction parameters characterizing the instantaneous probability distribution over molecular states in the much larger fine-scale models. The novel model reduction method, 'graph-constrained correlation dynamics', requires a graph of plausible state variables and interactions as input. It parametrically optimizes a set of constant coefficients appearing in differential equations governing the time-varying interaction parameters that determine all correlations between variables in the reduced model at any time slice
Dynamic Influence Networks for Rule-based Models
We introduce the Dynamic Influence Network (DIN), a novel visual analytics
technique for representing and analyzing rule-based models of protein-protein
interaction networks. Rule-based modeling has proved instrumental in developing
biological models that are concise, comprehensible, easily extensible, and that
mitigate the combinatorial complexity of multi-state and multi-component
biological molecules. Our technique visualizes the dynamics of these rules as
they evolve over time. Using the data produced by KaSim, an open source
stochastic simulator of rule-based models written in the Kappa language, DINs
provide a node-link diagram that represents the influence that each rule has on
the other rules. That is, rather than representing individual biological
components or types, we instead represent the rules about them (as nodes) and
the current influence of these rules (as links). Using our interactive DIN-Viz
software tool, researchers are able to query this dynamic network to find
meaningful patterns about biological processes, and to identify salient aspects
of complex rule-based models. To evaluate the effectiveness of our approach, we
investigate a simulation of a circadian clock model that illustrates the
oscillatory behavior of the KaiC protein phosphorylation cycle.Comment: Accepted to TVCG, in pres
- …