27,259 research outputs found
Sapo: Reachability Computation and Parameter Synthesis of Polynomial Dynamical Systems
Sapo is a C++ tool for the formal analysis of polynomial dynamical systems.
Its main features are: 1) Reachability computation, i.e., the calculation of
the set of states reachable from a set of initial conditions, and 2) Parameter
synthesis, i.e., the refinement of a set of parameters so that the system
satisfies a given specification. Sapo can represent reachable sets as unions of
boxes, parallelotopes, or parallelotope bundles (symbolic representation of
polytopes). Sets of parameters are represented with polytopes while
specifications are formalized as Signal Temporal Logic (STL) formulas
Identifying dynamical systems with bifurcations from noisy partial observation
Dynamical systems are used to model a variety of phenomena in which the
bifurcation structure is a fundamental characteristic. Here we propose a
statistical machine-learning approach to derive lowdimensional models that
automatically integrate information in noisy time-series data from partial
observations. The method is tested using artificial data generated from two
cell-cycle control system models that exhibit different bifurcations, and the
learned systems are shown to robustly inherit the bifurcation structure.Comment: 16 pages, 6 figure
Analysis of parametric biological models with non-linear dynamics
In this paper we present recent results on parametric analysis of biological
models. The underlying method is based on the algorithms for computing
trajectory sets of hybrid systems with polynomial dynamics. The method is then
applied to two case studies of biological systems: one is a cardiac cell model
for studying the conditions for cardiac abnormalities, and the second is a
model of insect nest-site choice.Comment: In Proceedings HSB 2012, arXiv:1208.315
Numerical algebraic geometry for model selection and its application to the life sciences
Researchers working with mathematical models are often confronted by the
related problems of parameter estimation, model validation, and model
selection. These are all optimization problems, well-known to be challenging
due to non-linearity, non-convexity and multiple local optima. Furthermore, the
challenges are compounded when only partial data is available. Here, we
consider polynomial models (e.g., mass-action chemical reaction networks at
steady state) and describe a framework for their analysis based on optimization
using numerical algebraic geometry. Specifically, we use probability-one
polynomial homotopy continuation methods to compute all critical points of the
objective function, then filter to recover the global optima. Our approach
exploits the geometric structures relating models and data, and we demonstrate
its utility on examples from cell signaling, synthetic biology, and
epidemiology.Comment: References added, additional clarification
Algebraic Systems Biology: A Case Study for the Wnt Pathway
Steady state analysis of dynamical systems for biological networks give rise
to algebraic varieties in high-dimensional spaces whose study is of interest in
their own right. We demonstrate this for the shuttle model of the Wnt signaling
pathway. Here the variety is described by a polynomial system in 19 unknowns
and 36 parameters. Current methods from computational algebraic geometry and
combinatorics are applied to analyze this model.Comment: 24 pages, 2 figure
Maximizing Protein Translation Rate in the Ribosome Flow Model: the Homogeneous Case
Gene translation is the process in which intracellular macro-molecules,
called ribosomes, decode genetic information in the mRNA chain into the
corresponding proteins. Gene translation includes several steps. During the
elongation step, ribosomes move along the mRNA in a sequential manner and link
amino-acids together in the corresponding order to produce the proteins.
The homogeneous ribosome flow model(HRFM) is a deterministic computational
model for translation-elongation under the assumption of constant elongation
rates along the mRNA chain. The HRFM is described by a set of n first-order
nonlinear ordinary differential equations, where n represents the number of
sites along the mRNA chain. The HRFM also includes two positive parameters:
ribosomal initiation rate and the (constant) elongation rate. In this paper, we
show that the steady-state translation rate in the HRFM is a concave function
of its parameters. This means that the problem of determining the parameter
values that maximize the translation rate is relatively simple. Our results may
contribute to a better understanding of the mechanisms and evolution of
translation-elongation. We demonstrate this by using the theoretical results to
estimate the initiation rate in M. musculus embryonic stem cell. The underlying
assumption is that evolution optimized the translation mechanism.
For the infinite-dimensional HRFM, we derive a closed-form solution to the
problem of determining the initiation and transition rates that maximize the
protein translation rate. We show that these expressions provide good
approximations for the optimal values in the n-dimensional HRFM already for
relatively small values of n. These results may have applications for synthetic
biology where an important problem is to re-engineer genomic systems in order
to maximize the protein production rate
A Computational Algebra Approach to the Reverse Engineering of Gene Regulatory Networks
This paper proposes a new method to reverse engineer gene regulatory networks
from experimental data. The modeling framework used is time-discrete
deterministic dynamical systems, with a finite set of states for each of the
variables. The simplest examples of such models are Boolean networks, in which
variables have only two possible states. The use of a larger number of possible
states allows a finer discretization of experimental data and more than one
possible mode of action for the variables, depending on threshold values.
Furthermore, with a suitable choice of state set, one can employ powerful tools
from computational algebra, that underlie the reverse-engineering algorithm,
avoiding costly enumeration strategies. To perform well, the algorithm requires
wildtype together with perturbation time courses. This makes it suitable for
small to meso-scale networks rather than networks on a genome-wide scale. The
complexity of the algorithm is quadratic in the number of variables and cubic
in the number of time points. The algorithm is validated on a recently
published Boolean network model of segment polarity development in Drosophila
melanogaster.Comment: 28 pages, 5 EPS figures, uses elsart.cl
- …