Dietrich Flockerzi

Subnetwork analysis reveals dynamic features of complex (bio)chemical networks

In analyzing and mathematical modeling of complex (bio)chemical reaction networks, formal methods that connect network structure and dynamic behavior are needed because often, quantitative knowledge of the networks is very limited. This applies to many important processes in cell biology. Chemical reaction network theory allows for the classification of the potential network behavior—for instance, with respect to the existence of multiple steady states—but is computationally limited to small systems. Here, we show that by analyzing subnetworks termed elementary flux modes, the applicability of the theory can be extended to more complex networks. For an example network inspired by cell cycle control in budding yeast, the approach allows for model discrimination, identification of key mechanisms for multistationarity, and robustness analysis. The presented methods will be helpful in modeling and analyzing other complex reaction networks.

Continuous Influenza Virus Production in Cell Culture Shows a Periodic Accumulation of Defective Interfering Particles

Influenza viruses are a major public health burden during seasonal epidemics and a continuous threat due to their potential to cause pandemics. Annual vaccination provides the best protection against the contagious respiratory illness caused by influenza viruses. However, the current production capacities for influenza vaccines are insufficient to meet the increasing demands. We explored the possibility to establish a continuous production process for influenza viruses using the duck-derived suspension cell line AGE1.CR. A two-stage bioreactor setup was designed in which cells were cultivated in a first stirred tank reactor where an almost constant cell concentration was maintained. Cells were then constantly fed to a second bioreactor where virus infection and replication took place. Using this two-stage reactor system, it was possible to continuously produce influenza viruses. Surprisingly, virus titers showed a periodic increase and decrease during the run-time of 17 days. These titer fluctuations were caused by the presence of defective interfering particles (DIPs), which we detected by PCR. Mathematical modeling confirmed this observation showing that constant virus titers can only emerge in the absence of DIPs. Even with very low amounts of DIPs in the seed virus and very low rates for de novo DIP generation, defective viruses rapidly accumulate and, therefore, represent a serious challenge for continuous vaccine production. Yet, the continuous replication of influenza virus using a two-stage bioreactor setup is a novel tool to study aspects of viral evolution and the impact of DIPs.

Multistationarity in mass action networks with applications to ERK activation

Ordinary Differential Equations (ODEs) are an important tool in many areas of Quantitative Biology. For many ODE systems multistationarity (i.e. the existence of at least two positive steady states) is a desired feature. In general establishing multistationarity is a difficult task as realistic biological models are large in terms of states and (unknown) parameters and in most cases poorly parameterized (because of noisy measurement data of few components, a very small number of data points and only a limited number of repetitions). For mass action networks establishing multistationarity hence is equivalent to establishing the existence of at least two positive solutions of a large polynomial system with unknown coefficients. For mass action networks with certain structural properties, expressed in terms of the stoichiometric matrix and the reaction rate-exponent matrix, we present necessary and sufficient conditions for multistationarity that take the form of linear inequality systems. Solutions of these inequality systems define pairs of steady states and parameter values. We also present a sufficient condition to identify networks where the aforementioned conditions hold. To show the applicability of our results we analyse an ODE system that is defined by the mass action network describing the extracellular signal-regulated kinase (ERK) cascade (i.e. ERK-activation).

Equilibrium and rate-based approaches to liquid–liquid phase splitting calculations

Abstract Several approaches are presented that can be used to compute liquid–liquid phase splitting. These approaches are based on homotopy continuation methods or on non-equilibrium thermodynamics. It is shown that they exhibit quite different computational demands with respect to computation times and a suggestion is made as to which algorithm to use. For chemical systems with only a few components, the computation times encountered on a normal PC in most cases are only a few milliseconds, in some cases even less—suggesting the general use of such algorithms to assure that phase splitting is always detected when present.

Multistationarity in Sequential Distributed Multisite Phosphorylation Networks

Multisite phosphorylation networks are encountered in many intracellular processes like signal transduction, cell-cycle control, or nuclear signal integration. In this contribution, networks describing the phosphorylation and dephosphorylation of a protein at n sites in a sequential distributive mechanism are considered. Multistationarity (i.e., the existence of at least two positive steady state solutions of the associated polynomial dynamical system) has been analyzed and established in several contributions. It is, for example, known that there exist values for the rate constants where multistationarity occurs. However, nothing else is known about these rate constants.Here, we present a sign condition that is necessary and sufficient for multistationarity in n-site sequential, distributive phosphorylation. We express this sign condition in terms of linear systems, and show that solutions of these systems define rate constants where multistationarity is possible. We then present, for n≥2, a collection of feasible linear systems, and hence give a new and independent proof that multistationarity is possible for n≥2. Moreover, our results allow to explicitly obtain values for the rate constants where multistationarity is possible. Hence, we believe that, for the first time, a systematic exploration of the region in parameter space where multistationarity occurs has become possible. One consequence of our work is that, for any pair of steady states, the ratio of the steady state concentrations of kinase-substrate complexes equals that of phosphatase-substrate complexes.

A thermodynamically consistent explicit competitive adsorption isotherm model based on second order single component behaviour

A competitive adsorption isotherm model is derived for binary mixtures of components characterized by single component isotherms which are second-order truncations of higher order equilibrium models suggested by multi-layer theory and statistical thermodynamics. The competitive isotherms are determined using the ideal adsorbed solution (IAS) theory which, in case of complex single component isotherms, does not generate explicit expressions to calculated equilibrium loadings and causes time consuming iterations in simulations of adsorption processes. The explicit model derived in this work is based on an analysis of the roots of a cubic polynomial resulting from the set of IAS equations. The suggested thermodynamically consistent and widely applicable competitive isotherm model can be recommended as a flexible tool for efficient simulations of fixed-bed adsorber dynamics.

Switching in Mass Action Networks Based on Linear Inequalities

Many biochemical processes can successfully be described by dynamical systems allowing some form of switching when, depending on their initial conditions, solutions of the dynamical system end up in different regions of state space (associated with different biochemical functions). Switching is often realized by a bistable system (i.e., a dynamical system allowing two stable steady state solutions) and, in the majority of cases, bistability is established numerically. In our view, this approach is too restrictive. On the one hand, due to predominant parameter uncertainty, numerical methods are generally difficult to apply to realistic models originating in systems biology. On the other hand, switching already arises with the occurrence of a saddle-type steady state (characterized by a Jacobian where exactly one eigenvalue is positive and the remaining eigenvalues have negative real part). Consequently we derive conditions based on linear inequalities that allow the analytic computation of states and param...

Size-exclusion chromatography as a linear transfer system: Purification of human influenza virus as an example

Preparative size-exclusion chromatography suffers from low selectivity and productivity. Empirical optimization of operating conditions constitutes a laborious task due to many parameters. Here, a modeling framework based on linear systems theory is presented for predicting the influence of volume overloading. Impulse-responses characterizing system behavior are derived from experimental data by maximum entropy deconvolution. Theoretical derivations are validated experimentally by study of a model system and chromatography of human influenza virus. By application of the theory it is demonstrated how group separation operations can be optimized with respect to yield, purity, productivity and dilution of the product.

Nested autoinhibitory feedbacks alter the resistance of homeostatic adaptive biochemical networks

Negative feedback control is a ubiquitous feature of biochemical systems, as is time delay between a signal and its response. Negative feedback in conjunction with time delay can lead to oscillations. In a cellular context, it might be beneficial to mitigate oscillatory behaviour to avoid recurring stress situations. This can be achieved by increasing the distance between the parameters of the system and certain thresholds, beyond which oscillations occur. This distance has been termed resistance. Here, we prove that in a generic three-dimensional negative feedback system the resistance of the system is modified by nested autoinhibitory feedbacks. Our system features negative feedbacks through both input-inhibition as well as output-activation, a signalling component with mass conservation and perfect adaptation. We show that these features render the system applicable to biological data, exemplified by the high osmolarity glycerol system in yeast and the mammalian p53 system. Output-activation is better supported by data than input-inhibition and also shows distinguished properties with respect to the systems stimulus. Our general approach might be useful in designing synthetic systems in which oscillations can be tuned by synthetic autoinhibitory feedbacks.

Modeling the Light- and Redox-Dependent Interaction of PpsR/AppA in Rhodobacter sphaeroides

Facultative photosynthetic bacteria switch their energy generation mechanism from respiration to photosynthesis depending on oxygen tension and light. Part of this transition is mediated by the aerobic transcriptional repressor PpsR. In Rhodobacter sphaeroides, the repressive action of PpsR is antagonized by the redox- and blue-light-sensitive flavoprotein AppA which results in a unique phenotype: the repression of photosynthesis genes at intermediate oxygen levels and high light intensity, which is believed to reduce the risk of photooxidative stress. To analyze the underlying mechanism we developed a simple mathematical model based on the AppA-dependent reduction of a disulfide bond in PpsR and the light-sensitive complex formation between the reduced forms of AppA and PpsR. A steady-state analysis shows that high light repression can indeed occur at intermediate oxygen levels if PpsR is reduced on a faster timescale than AppA and if the electron transfer from AppA to PpsR is effectively irreversible. The model further predicts that if AppA copy numbers exceed those of PpsR by at least a factor of two, the transition from aerobic to anaerobic growth mode can occur via a bistable regime. We provide necessary conditions for the emergence of bistability and discuss possible experimental verifications.