Source author record

Steffen Waldherr

Steffen Waldherr appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

9works
9topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

9 published item(s)

preprint2020arXiv

Numerical Gaussian process Kalman filtering

In this manuscript we introduce numerical Gaussian process Kalman filtering (GPKF). Numerical Gaussian processes have recently been developed to simulate spatiotemporal models. The contribution of this paper is to embed numerical Gaussian processes into the recursive Kalman filter equations. This embedding enables us to do Kalman filtering on infinite-dimensional systems using Gaussian processes. This is possible because i) we are obtaining a linear model from numerical Gaussian processes, and ii) the states of this model are by definition Gaussian distributed random variables. Convenient properties of the numerical GPKF are that no spatial discretization of the model is necessary, and manual setting up of the Kalman filter, that is fine-tuning the process and measurement noise levels by hand is not required, as they are learned online from the data stream. We showcase the capability of the numerical GPKF in a simulation study of the advection equation.

preprint2015arXiv

Ensemble Observability of Linear Systems

We address the observability problem for ensembles that are described by probability distributions. The problem is to reconstruct a probability distribution of the initial state from the time-evolution of the probability distribution of the output under a classical finite-dimensional linear system. We present two solutions to this problem, one based on formulating the problem as an inverse problem and the other one based on reconstructing all the moments of the distribution. The first approach leads us to a connection between the reconstruction problem and mathematical tomography problems. In the second approach we use the framework of tensor systems to describe the dynamics of the moments which leads to a more systems theoretic treatment of the reconstruction problem. Furthermore we show that both frameworks are inherently related. The appeal of having two dual view points, the first being more geometric and the second one being more systems theoretic, is illuminated in several examples of theoretical or practical importance.

preprint2014arXiv

Dynamic optimization of metabolic networks coupled with gene expression

The regulation of metabolic activity by tuning enzyme expression levels is crucial to sustain cellular growth in changing environments. Metabolic networks are often studied at steady state using constraint-based models and optimization techniques. However, metabolic adaptations driven by changes in gene expression cannot be analyzed by steady state models, as these do not account for temporal changes in biomass composition. Here we present a dynamic optimization framework that integrates the metabolic network with the dynamics of biomass production and composition, explicitly taking into account enzyme production costs and enzymatic capacity. In contrast to the established dynamic flux balance analysis, our approach allows predicting dynamic changes in both the metabolic fluxes and the biomass composition during metabolic adaptations. We applied our algorithm in two case studies: a minimal nutrient uptake network, and an abstraction of core metabolic processes in bacteria. In the minimal model, we show that the optimized uptake rates reproduce the empirical Monod growth for bacterial cultures. For the network of core metabolic processes, the dynamic optimization algorithm predicted commonly observed metabolic adaptations, such as a diauxic switch with a preference ranking for different nutrients, re-utilization of waste products after depletion of the original substrate, and metabolic adaptation to an impending nutrient depletion. These examples illustrate how dynamic adaptations of enzyme expression can be predicted solely from an optimization principle.

preprint2014arXiv

Generalization of the construction method for multistability-equivalent gene regulatory networks to systems with multi-input multi-output loopbreaking

The problem of equivalence in terms of multistability properties between gene regulatory network models of different dimensionality has been recently addressed by Schittler et al. (2013). The authors in that work proposed construction rules for a high-dimensional dynamical system, when given a low-dimensional dynamical system and the high-dimensional network structure. However, the proof therein was restricted to the class of systems for which all internal feedback loops can be broken by a loopbreaking approach yielding a single-input single-output (SISO) system. In this report, we present the generalization of the proof to systems with any number of internal feedback loops, which will be broken by a generalized loopbreaking approach resulting in a multi-input multi-output (MIMO) system. This generalization of the method renders the construction method applicable to a broad class of gene regulatory network models, thus promoting the transfer of results from core motif models to more realistic, high-dimensional models of gene regulation. We demonstrate the potential and value of our method by applying it to an example of a gene regulatory network in mesenchymal stem cell differentiation.

preprint2012arXiv

Efficient parametric analysis of the chemical master equation through model order reduction

Background: Stochastic biochemical reaction networks are commonly modelled by the chemical master equation, and can be simulated as first order linear differential equations through a finite state projection. Due to the very high state space dimension of these equations, numerical simulations are computationally expensive. This is a particular problem for analysis tasks requiring repeated simulations for different parameter values. Such tasks are computationally expensive to the point of infeasibility with the chemical master equation. Results: In this article, we apply parametric model order reduction techniques in order to construct accurate low-dimensional parametric models of the chemical master equation. These surrogate models can be used in various parametric analysis task such as identifiability analysis, parameter estimation, or sensitivity analysis. As biological examples, we consider two models for gene regulation networks, a bistable switch and a network displaying stochastic oscillations. Conclusions: The results show that the parametric model reduction yields efficient models of stochastic biochemical reaction networks, and that these models can be useful for systems biology applications involving parametric analysis problems such as parameter exploration, optimization, estimation or sensitivity analysis.

preprint2012arXiv

Network-level dynamics of diffusively coupled cells

We study molecular dynamics within populations of diffusively coupled cells under the assumption of fast diffusive exchange. As a technical tool, we propose conditions on boundedness and ultimate boundedness for systems with a singular perturbation, which extend the classical asymptotic stability results for singularly perturbed systems. Based on these results, we show that with common models of intracellular dynamics, the cell population is coordinated in the sense that all cells converge close to a common equilibrium point. We then study a more specific example of coupled cells which behave as bistable switches, where the intracellular dynamics are such that cells may be in one of two equilibrium points. Here, we find that the whole population is bistable in the sense that it converges to a population state where either all cells are close to the one equilibrium point, or all cells are close to the other equilibrium point. Finally, we discuss applications of these results for the robustness of cellular decision making in coupled populations.

preprint2011arXiv

Structural requirements and discrimination of cell differentiation networks

Mathematical models of stem cell differentiation are commonly based upon the concept of subsequent cell fate decisions, each controlled by a gene regulatory network. These networks exhibit a multistable behavior and cause the system to switch between qualitatively distinct stable steady states. However, the network structure of such a switching module is often uncertain, and there is lack of knowledge about the exact reaction kinetics. In this paper, we therefore perform an elementary study of small networks consisting of three interacting transcriptional regulators responsible for cell differentiation: We investigate which network structures can reproduce a certain multistable behavior, and how robustly this behavior is realized by each network. In order to approach these questions, we use a modeling framework which only uses qualitative information about the network, yet allows model discrimination as well as to evaluate the robustness of the desired multistability properties. We reveal structural network properties which are necessary and sufficient to realize distinct steady state patterns required for cell differentiation. Our results also show that structural and robustness properties of the networks are related to each other.

preprint2010arXiv

Bridging Time Scales in Cellular Decision Making with a Stochastic Bistable Switch

Cellular transformations which involve a significant phenotypical change of the cell's state use bistable biochemical switches as underlying decision systems. In this work, we aim at linking cellular decisions taking place on a time scale of years to decades with the biochemical dynamics in signal transduction and gene regulation, occuring on a time scale of minutes to hours. We show that a stochastic bistable switch forms a viable biochemical mechanism to implement decision processes on long time scales. As a case study, the mechanism is applied to model the initiation of follicle growth in mammalian ovaries, where the physiological time scale of follicle pool depletion is on the order of the organism's lifespan. We construct a simple mathematical model for this process based on experimental evidence for the involved genetic mechanisms. Despite the underlying stochasticity, the proposed mechanism turns out to yield reliable behavior in large populations of cells subject to the considered decision process. Our model explains how the physiological time constant may emerge from the intrinsic stochasticity of the underlying gene regulatory network. Apart from ovarian follicles, the proposed mechanism may also be of relevance for other physiological systems where cells take binary decisions over a long time scale.

preprint2010arXiv

Searching bifurcations in high-dimensional parameter space via a feedback loop breaking approach

Bifurcations leading to complex dynamical behaviour of non-linear systems are often encountered when the characteristics of feedback circuits in the system are varied. In systems with many unknown or varying parameters, it is an interesting, but difficult problem to find parameter values for which specific bifurcations occur. In this paper, we develop a loop breaking approach to evaluate the influence of parameter values on feedback circuit characteristics. This approach allows a theoretical classification of feedback circuit characteristics related to possible bifurcations in the system. Based on the theoretical results, a numerical algorithm for bifurcation search in a possibly high-dimensional parameter space is developed. The application of the proposed algorithm is illustrated by searching for a Hopf bifurcation in a model of the mitogen activated protein kinase (MAPK) cascade, which is a classical example for biochemical signal transduction.