Interpreting Frequency Responses to Dose-Conserved Pulsatile Input Signals in Simple Cell Signaling Motifs

Many hormones are released in pulsatile patterns. This pattern can be modified, for instance by changing pulse frequency, to encode relevant physiological information. Often other properties of the pulse pattern will also change with frequency. How do signaling pathways of cells targeted by these hormones respond to different input patterns? In this study, we examine how a given dose of hormone can induce different outputs from the target system, depending on how this dose is distributed in time. We use simple mathematical models of feedforward signaling motifs to understand how the properties of the target system give rise to preferences in input pulse pattern. We frame these problems in terms of frequency responses to pulsatile inputs, where the amplitude or duration of the pulses is varied along with frequency to conserve input dose. We find that the form of the nonlinearity in the steady state input-output function of the system predicts the optimal input pattern. It does so by selecting an optimal input signal amplitude. Our results predict the behavior of common signaling motifs such as receptor binding with dimerization, and protein phosphorylation. The findings have implications for experiments aimed at studying the frequency response to pulsatile inputs, as well as for understanding how pulsatile patterns drive biological responses via feedforward signaling pathways.     

Mathematical analysis of a metapopulation model with space-limited recruitment

We investigate a mathematical aspect of a multi-species' sessile metapopulation model with space-limited recruitment proposed by Iwasa et al. in 1986. We define some basic reproduction numbers to show the threshold condition for the stability of trivial steady state and the existence of coexistent steady state. We show the existence of steady state where all species exist when some reproduction numbers are greater than one by the fixed point theorem. And we construct the Lyapunov function to show the global stability of trivial steady state when some basic reproduction numbers are not greater than one.     

Numerical bifurcation analysis of delay differential equations arising from physiological modeling

This paper has a dual purpose. First, we describe numerical methods for continuation and bifurcation analysis of steady state solutions and periodic solutions of systems of delay differential equations with an arbitrary number of fixed, discrete delays. Second, we demonstrate how these methods can be used to obtain insight into complex biological regulatory systems in which interactions occur with time delays: for this, we consider a system of two equations for the plasma glucose and insulin concentrations in a diabetic patient subject to a system of external assistance. The model has two delays: the technological delay of the external system, and the physiological delay of the patient's liver. We compute stability of the steady state solution as a function of two parameters, compare with analytical results and compute several branches of periodic solutions and their stability. These numerical results allow to infer two categories of diabetic patients for which the external system has different efficiency.     

Oscillatory viral dynamics in a delayed HIV pathogenesis model

We consider an HIV pathogenesis model incorporating antiretroviral therapy and HIV replication time. We investigate the existence and stability of equilibria, as well as Hopf bifurcations to sustained oscillations when drug efficacy is less than 100%. We derive sufficient conditions for the global asymptotic stability of the uninfected steady state. We show that time delay has no effect on the local asymptotic stability of the uninfected steady state, but can destabilize the infected steady state, leading to a Hopf bifurcation to periodic solutions in the realistic parameter ranges.     

Separable models in age-dependent population dynamics

A class of population models is considered in which the parameters such as fecundity, mortality and interaction coefficients are assumed to be age-dependent. Conditions for the existence, stability and global attractivity of steady-state and periodic solutions are derived. The dependence of these solutions on the maturation periods is analyzed. These results are applied to specific single and multiple population models. It is shown that periodic solutions cannot occur in a general class of single population age-dependent models. Conditions are derived that determine whether increasing the maturation period has a stabilizing effect. In specific cases, it is shown that any number of switches in stability can occur as the maturation period is increased. An example is given of predator-prey model where each one of these stability switches corresponds to a stable steady state losing its stability via a Hopf bifurcation to a periodic solution and regaining its stability upon further increase of the maturation period.     

Spatial and spatio-temporal patterns in a cell-haptotaxis model

We investigate a cell-haptotaxis model for the generation of spatial and spatio-temporal patterns in one dimension. We analyse the steady state problem for specific boundary conditions and show the existence of spatially hetero-geneous steady states. A linear analysis shows that stability is lost through a Hopf bifurcation. We carry out a nonlinear multi-time scale perturbation procedure to study the evolution of the resulting spatio-temporal patterns. We also analyse the model in a parameter domain wherein it exhibits a singular dispersion relation.     

Mathematical analysis of an age-structured population model with space-limited recruitment

In this paper, we investigate structured population model of marine invertebrate whose life stage is composed of sessile adults and pelagic larvae, such as barnacles contained in a local habitat. First we formulate the basic model as an Cauchy problem on a Banach space to discuss the existence and uniqueness of non-negative solution. Next we define the basic reproduction number R0 to formulate the invasion condition under which the larvae can successfully settle down in the completely vacant habitat. Subsequently we examine existence and stability of steady states. We show that the trivial steady state is globally asymptotically stable if R0 < or = 1, whereas it is unstable if R0 > 1. Furthermore, we show that a positive (non-trivial) steady state uniquely exists if R0 > 1 and it is locally asymptotically stable as far as absolute value of R0 - 1 is small enough.     

Non-Linear Bifurcation Analysis of Reaction-Diffusion Activator-Inhibator System

The paper first deals with the linear stability analysis of an activator-inhibitor reaction diffusion system to determine the nature of the bifurcation point of the system. The non-linear bifurcation analysis determining the steady state solution beyond the critical point enables us to determine characteristic features of the spatial inhomogeneous pattern formation arising out of the bifurcation of the state of the system.     

Explosive route to chaos through a fractal torus in a generalized Lotka-Volterra model

The behavior of a model that generalizes the Lotka-Volterra problem into three dimensions is presented. The results show the analytic derivation of stability diagrams that describe the system's qualitative features. In particular, we show that for a certain value of the bifurcation parameter the system instantly jumps out of a steady state solution into a chaotic solution that portrays a fractal torus in the three-dimensional phase space. This scenario, is referred to as the explosive route to chaos and is attributed to the non-transversal saddle connection type bifurcation. The stability diagrams also present a region in which the Hopf type bifurcation leads to periodic and chaotic solutions. In addition, the bifurcation diagrams reveal a qualitative similarity to the data obtained in the Texas and Bordeaux experiments on the Belousov-Zhabotinskii chemical reaction. The paper is concluded by showing that the model can be useful for representing dynamics associated with biological and chemical phenomena.     

Existence and stability of equilibria in age-structured population dynamics

The existence of positive equilibrium solutions of the McKendrick equations for the dynamics of an age-structured population is studied as a bifurcation phenomenon using the inherent net reproductive rate n as a bifurcation parameter. The local existence and uniqueness of a branch of positive equilibria which bifurcates from the trivial (identically zero) solution at the critical value n=1 are proved by implicit function techniques under very mild smoothness conditions on the death and fertility rates as functional of age and population density. This first requires the development of a suitable linear theory. The lowest order terms in the Liapunov-Schmidt expansions are also calculated. This local analysis supplements earlier global bifurcation results of the author. The stability of both the trivial and the positive branch equilibria is studied by means of the principle of linearized stability. It is shown that in general the trivial solution losses stability as n increases through one while the stability of the branch solution is stable if and only if the bifurcation is supercritical. Thus the McKendrick equations exhibit, in the latter case, a standard exchange of stability with regard to equilibrium states as they depend on the inherent net reproductive rate. The derived lower order terms in the Liapunov-Schmidt expansions yield formulas which explicitly relate the direction of bifurcation to properties of the age-specific death and fertility rates as functionals of population density. Analytical and numerical results for some examples are given which illustrate these results.     

Growth of nonnecrotic tumors in the presence and absence of inhibitors

In this article a model for the evolution of a spherically symmetric, nonnecrotic tumor is presented. The effects of nutrients and inhibitors on the existence and stability of time-independent solutions are studied. With a single nutrient and no inhibitors present, the trivial solution, which corresponds to a state in which no tumor is present, persists for all parameter values, whereas the nontrivial solution, which corresponds to a tumor of finite size, exists for only a prescribed range of parameters, which corresponds to a balance between cell proliferation and cell death. Stability analysis, based on a two-timing method, suggests that, where it exists, the nontrivial solution is stable and the trivial solution unstable. Otherwise, the trivial solution is stable. Modification to these characteristic states brought about by the presence of different types of inhibitors are also investigated and shown to have significant effect. Implications of the model for the treatment of cancer are also discussed.     

Selecting a common direction

We present results and analysis of models for contact-induced turning responses and alignment in populations of interacting individuals. Such models describe distributions of orientation, and how these evolve under different assumptions about the turning behaviour of individuals. One of these models was first used to describe interactions between mammalian cells called fibroblasts in Edelstein-Keshet and Ermentrout (1990) J. Math. Biol. 29: 33–58 (henceforth abbreviated EKE 1990). Here the original model is generalized to encompass motion in both 2 and 3 dimensions. Two modifications of this model are introduced: in one, the turning is described by a gradual direction shift (rather than abrupt transition). In another variant, the interactions between individuals changes as the density of the population increases to include the effects of crowding. Using linear stability analysis and synergetics analysis of interacting modes we describe the nature and stability properties of the steady state solutions. We investigate how nonhomogeneous pattern evolves close to the bifurcation point. We find that individuals tend to cluster together in one direction of alignment.     

Phase Transitions in a Logistic Metapopulation Model with Nonlocal Interactions

The presence of one or more species at some spatial locations but not others is a central matter in ecology. This phenomenon is related to ecological pattern formation. Nonlocal interactions can be considered as one of the mechanisms causing such a phenomenon. We propose a single-species, continuous time metapopulation model taking nonlocal interactions into account. Discrete probability kernels are used to model these interactions in a patchy environment. A linear stability analysis of the model shows that solutions to this equation exhibit pattern formation if the dispersal rate of the species is sufficiently small and the discrete interaction kernel satisfies certain conditions. We numerically observe that traveling and stationary wave-type patterns arise near critical dispersal rate. We use weakly nonlinear analysis to better understand the behavior of formed patterns. We show that observed patterns arise through both supercritical and subcritical bifurcations from spatially homogeneous steady state. Moreover, we observe that as the dispersal rate decreases, amplitude of the patterns increases. For discontinuous transitions to instability, we also show that there exists a threshold for the amplitude of the initial condition, above which pattern formation is observed.     

A condition for setting off ectopic waves in computational models of excitable cells

The purpose of this paper is to study the stability of steady state solutions of the Monodomain model equipped with Luo-Rudy I kinetics. It is well established that re-entrant arrhythmias can be created in computational models of excitable cells. Such arrhythmias can be initiated by applying an external stimulus that interacts with a partially refractory region, and spawn breaking waves that can eventually generate extremely complex wave patterns commonly referred to as fibrillation. An ectopic wave is one possible stimulus that may initiate fibrillation. Physiologically, it is well known that ectopic waves exist, but the mechanism for initiating ectopic waves in a large collection of cells is poorly understood. In the present paper we consider computational models of collections of excitable cells in one and two spatial dimensions. The cells are modeled by Luo-Rudy I kinetics, and we assume that the spatial dynamics is governed by the Monodomain model. The mathematical analysis is carried out for a reduced model that is known to provide good approximations of the initial phase of solutions of the Luo-Rudy I model. A further simplification is also introduced to motivate and explain the results for the more complicated models. In the analysis the cells are divided into two regions; one region (N) consists of normal cells as model by the standard Luo-Rudy I model, and another region (A) where the cells are automatic in the sense that they would act as pacemaker cells if they where isolated from their surroundings. We let delta denote the spatial diffusion and a denote a characteristic length of the automatic region. It has previously been shown that reducing diffusion or increasing the automatic region enhances ectopic activity. Here we derive a condition for the transition from stable resting state to ectopic wave spread. Under suitable assumptions on the model we provide mathematical and computational arguments indicating that there is a constant eta such that a steady state solution of this system is stable whenever delta approximately > etaa(2), and unstable whenever delta approximately < etaa(2).     

Phase resetting and bifurcation in the ventricular myocardium.

With the dynamic differential equations of Beeler, G. W., and H. Reuter (1977, J. Physiol. [Lond.]. 268:177-210), we have studied the oscillatory behavior of the ventricular muscle fiber stimulated by a depolarizing applied current I app. The dynamic solutions of BR equations revealed that as I app increases, a periodic repetitive spiking mode appears above the subthreshold I app, which transforms to a periodic spiking-bursting mode of oscillations, and finally to chaos near the suprathreshold I app (i.e., near the termination of the periodic state). Phase resetting and annihilation of repetitive firing in the ventricular myocardium were demonstrated by a brief current pulse of the proper magnitude applied at the proper phase. These phenomena were further examined by a bifurcation analysis. A bifurcation diagram constructed as a function of I app revealed the existence of a stable periodic solution for a certain range of current values. Two Hopf bifurcation points exist in the solution, one just above the lower periodic limit point and the other substantially below the upper periodic limit point. Between each periodic limit point and the Hopf bifurcation, the cell exhibited the coexistence of two different stable modes of operation; the oscillatory repetitive firing state and the time-independent steady state. As in the Hodgkin-Huxley case, there was a low amplitude unstable periodic state, which separates the domain of the stable periodic state from the stable steady state. Thus, in support of the dynamic perturbation methods, the bifurcation diagram of the BR equation predicts the region where instantaneous perturbations, such as brief current pulses, can send the stable repetitive rhythmic state into the stable steady state.     

Transient spatio-temporal dynamics of a diffusive plant–herbivore system with Neumann boundary conditions

In many existing predator–prey or plant–herbivore models, the numerical response is assumed to be proportional to the functional response. In this paper, without such an assumption, we consider a diffusive plant–herbivore system with Neumann boundary conditions. Besides stability of spatially homogeneous steady states, we also derive conditions for the occurrence of Hopf bifurcation and steady-state bifurcation and provide geometrical methods to locate the bifurcation values. We numerically explore the complex transient spatio-temporal behaviours induced by these bifurcations. A large variety of different types of transient behaviours including oscillations in one or both of space and time are observed.     

Periodic oscillations in leukopoiesis models with two delays

The term leukopoiesis describes processes leading to the production and regulation of white blood cells. It is based on stem cells differentiation and may exhibit abnormalities resulting in severe diseases, such as cyclical neutropenia and leukemias. We consider a nonlinear system of two equations, describing the evolution of a stem cell population and the resulting white blood cell population. Two delays appear in this model to describe the cell cycle duration of the stem cell population and the time required to produce white blood cells. We establish sufficient conditions for the asymptotic stability of the unique nontrivial positive steady state of the model by analysing roots of a second degree exponential polynomial characteristic equation with delay-dependent coefficients. We also prove the existence of a Hopf bifurcation which leads to periodic solutions. Numerical simulations of the model with parameter values reported in the literature demonstrate that periodic oscillations (with short and long periods) agree with observations of cyclical neutropenia in patients.     

Solution of a multinephron,multisolute model of the mammalian kidney by Newton and continuation methods

A method of numerically solving the differential equations specifying solute and water flow in a multinephron, multisolute model of the mammalian kidney by a combination of Newton and continuation techniques is described. This method is used to generate a connected component of the steady state solution manifold of the model. A three dimensional section of this manifold is shown to be convoluted, with upper and lower sheets of stable solutions connected by an unstable middle sheet. Two dimensional sections of this surface are followed from a trivial constant profile of concentrations in the nonconcentrating kidney to the profiles of the maximally concentrating kidney. Study of these sections shows that for a given choice of model parameters there may exist no solution, there may be a unique solution, or there may be multiple solutions. A study of the time dependent solutions shows that the dynamic transition from the lower to the upper state and return may be via a hysteresis loop.     

Analysis of operating principles with S-system models

Operating principles address general questions regarding the response dynamics of biological systems as we observe or hypothesize them, in comparison to a priori equally valid alternatives. In analogy to design principles, the question arises: Why are some operating strategies encountered more frequently than others and in what sense might they be superior? It is at this point impossible to study operation principles in complete generality, but the work here discusses the important situation where a biological system must shift operation from its normal steady state to a new steady state. This situation is quite common and includes many stress responses. We present two distinct methods for determining different solutions to this task of achieving a new target steady state. Both methods utilize the property of S-system models within Biochemical Systems Theory (BST) that steady states can be explicitly represented as systems of linear algebraic equations. The first method uses matrix inversion, a pseudo-inverse, or regression to characterize the entire admissible solution space. Operations on the basis of the solution space permit modest alterations of the transients toward the target steady state. The second method uses standard or mixed integer linear programming to determine admissible solutions that satisfy criteria of functional effectiveness, which are specified beforehand. As an illustration, we use both methods to characterize alternative response patterns of yeast subjected to heat stress, and compare them with observations from the literature.     

Classification and stability of global inhomogeneous solutions of a macroscopic model of cell motion

Many micro-organisms use chemotaxis for aggregation, resulting in stable patterns. In this paper, the amoeba Dictyostelium discoideum serves as a model organism for understanding the conditions for aggregation and classification of resulting patterns. To accomplish this, a 1D nonlinear diffusion equation with chemotaxis that models amoeba behavior is analyzed. A classification of the steady state solutions is presented, and a Lyapunov functional is used to determine conditions for stability of inhomogenous solutions. Changing the chemical sensitivity, production rate of the chemical attractant, or domain length can cause the system to transition from having an asymptotic steady state, to having asymptotically stable single-step solution and multi-stepped stable plateau solutions.