Refine
Year of publication
Document Type
- Article (31)
- ZIB-Report (15)
- In Proceedings (7)
- Poster (2)
- Book (1)
- Book chapter (1)
- In Collection (1)
Is part of the Bibliography
- no (58)
Keywords
- systems biology (4)
- differential equations (3)
- Bayesian inference (2)
- DS-MLE (2)
- EM algorithm (2)
- Jeffreys prior (2)
- MPLE (2)
- NPMLE (2)
- cow (2)
- hormone patterns (2)
Institute
- Computational Systems Biology (58) (remove)
Background.
The chemical master equation is the fundamental equation of stochastic chemical kinetics. This differential-difference equation describes temporal evolution of the probability density function for states of a chemical system. A state of the system, usually encoded as a vector, represents the number of entities or copy numbers of interacting species, which are changing according to a list of possible reactions. It is often the case, especially when the state vector is high-dimensional, that the number of possible states the system may occupy is too large to be handled computationally. One way to get around this problem is to consider only those states that are associated with probabilities that are greater than a certain threshold level.
Results.
We introduce an algorithm that significantly reduces computational resources and is especially powerful when dealing with multi-modal distributions. The algorithm is built according to two key principles. Firstly, when performing time integration, the algorithm keeps track of the subset of states with significant probabilities (essential support). Secondly, the probability distribution that solves the equation is parametrised with a small number of coefficients using collocation on Gaussian radial basis functions. The system of basis functions is chosen in such a way that the solution is approximated only on the essential support instead of the whole state space.
Discussion.
In order to demonstrate the effectiveness of the method, we consider four application examples: a) the self-regulating gene model, b) the 2-dimensional bistable toggle switch, c) a generalisation of the bistable switch to a 3-dimensional tristable problem, and d) a 3-dimensional cell differentiation model that, depending on parameter values, may operate in bistable or tristable modes. In all multidimensional examples the manifold containing the system states with significant probabilities undergoes drastic transformations over time. This fact makes the examples especially challenging for numerical methods.
Conclusions.
The proposed method is a new numerical approach permitting to approximately solve a wide range of problems that have been hard to tackle until now. A full representation of multi-dimensional distributions is recovered. The method is especially attractive when dealing with models that yield solutions of a complex structure, for instance, featuring multi-stability.
Electronic version: http://www.biomedcentral.com/1752-0509/9/67
One of the main goals of mathematical modelling in systems biology related to medical applications is to obtain patient-specific parameterisations and model predictions.
In clinical practice, however, the number of available measurements for single patients is usually limited due to time and cost restrictions. This hampers the process of making patient-specific predictions about the outcome of a treatment. On the other hand, data are often available for many patients, in particular if extensive clinical studies have been performed. Using these population data, we propose an iterative algorithm for contructing an informative prior distribution, which then serves as the basis for computing patient-specific posteriors and obtaining individual predictions. We demonsrate the performance of our method by applying it to a low-dimensional parameter estimation problem in a toy model as well as to a high-dimensional ODE model of the human menstrual cycle, which represents a typical example from systems biology modelling.
A deterministic model of tuberculosis in Cameroon is designed and analyzed with respect to its transmission dynamics. The model includes lack of access to treatment and weak diagnosis capacity as well as both frequency- and density-dependent transmissions. It is shown that the model is mathematically well-posed and epidemiologically reasonable. Solutions are non-negative and bounded whenever the initial values are non-negative. A sensitivity analysis of model parameters is performed and the most sensitive ones are identified by means of a state-of-the-art Gauss-Newton method. In particular, parameters representing the proportion of individuals having access to medical facilities are seen to have a large impact on the dynamics of the disease. The model predicts that a gradual increase of these parameters could significantly reduce the disease burden on the population within the next 15 years.
A deterministic model of tuberculosis in sub-Saharan Africa in general and Cameroon in particular including lack of access to the treatment and weak diagnose capacity is designed and analyzed with respect to its transmission dynamics.
The model includes both frequency- and density-dependent transmissions. It is shown that the model is mathematically well-posed and epidemiologically reasonable. Solutions are non-negative and bounded whenever the initial values are non-negative.
A sensitivity analysis of model parameters is performed and most sensitive parameters of the model are identified using a state-of-the-art Gauss-Newton
Method. In particular, parameters representing the proportion of individuals having access to medical facilities have a large impact on the dynamics of the disease. It has been shown that an increase of these parameter values over the
time can significantly reduce the disease burden in the population within the next 15 years.
When dealing with Bayesian inference the choice of the prior often remains a debatable question. Empirical Bayes methods offer a data-driven solution to this problem by estimating the prior itself from an ensemble of data. In the nonparametric case, the maximum likelihood estimate is known to overfit the data, an issue that is commonly tackled by regularization. However, the majority of regularizations are ad hoc choices which lack invariance under reparametrization of the model and result in inconsistent estimates for equivalent models. We introduce a nonparametric, transformation-invariant estimator for the prior distribution. Being defined in terms of the missing information similar to the reference prior, it can be seen as an extension of the latter to the data-driven setting. This implies a natural interpretation as a trade-off between choosing the least informative prior and incorporating the information provided by the data, a symbiosis between the objective and empirical Bayes methodologies.
Potassium is fundamental for cell functioning including signal transduction, acid-base- and water-metabolism.
Since diet of dairy cows is generally rich in potassium, hypokalemia was not in the focus of research for long time. Furthermore, hypokalemia was not frequently diagnosed because blood potassium content is difficult to measure. In recent years, measurement methods have been improved. Nowadays hypokalemia is increasingly diagnosed in cows with disorders such as abomasal displacement, ketosis or down cow syndrome, calling for intensified research on this topic.
In this report we describe the development of a basic mechanistic, dynamic model of potassium balance based on ordinary differential and algebraic equations. Parameter values are obtained from data of a clinical trial in which potassium balance and the influence of therapeutic intervention in glucose and electrolyte metabolism on potassium balance in non-lactating dairy cows were studied. The model is formulated at a high abstraction level and includes information and hypotheses from literature. This work represents a first step towards the understanding and design of effective prophylactic feed additives and treatment strategies.