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.     

Slow Passage Through a Hopf Bifurcation in Excitable Nerve Cables: Spatial Delays and Spatial Memory Effects

It is well established that in problems featuring slow passage through a Hopf bifurcation (dynamic Hopf bifurcation) the transition to large-amplitude oscillations may not occur until the slowly changing parameter considerably exceeds the value predicted from the static Hopf bifurcation analysis (temporal delay effect), with the length of the delay depending upon the initial value of the slowly changing parameter (temporal memory effect). In this paper we introduce new delay and memory effect phenomena using both analytic (WKB method) and numerical methods. We present a reaction–diffusion system for which slowly ramping a stimulus parameter (injected current) through a Hopf bifurcation elicits large-amplitude oscillations confined to a location a significant distance from the injection site (spatial delay effect). Furthermore, if the initial current value changes, this location may change (spatial memory effect). Our reaction–diffusion system is Baer and Rinzel’s continuum model of a spiny dendritic cable; this system consists of a passive dendritic cable weakly coupled to excitable dendritic spines. We compare results for this system with those for nerve cable models in which there is stronger coupling between the reactive and diffusive portions of the system. Finally, we show mathematically that Hodgkin and Huxley were correct in their assertion that for a sufficiently slow current ramp and a sufficiently large cable length, no value of injected current would cause their model of an excitable cable to fire; we call this phenomenon “complete accommodation.”     

Bifurcations in a white-blood-cell production model

We study the dynamics of a model of white-blood-cell (WBC) production. The model consists of two compartmental differential equations with two discrete delays. We show that from normal to pathological parameter values, the system undergoes supercritical Hopf bifurcations and saddle-node bifurcations of limit cycles. We characterize the steady states of the system and perform a bifurcation analysis. Our results indicate that an increase in apoptosis rate of either hematopoietic stem cells or WBC precursors induces a Hopf bifurcation and an oscillatory regime takes place. These oscillations are seen in some hematological diseases.     

Diffusively coupled bursters: Effects of cell heterogeneity

The interaction of a pair of weakly coupled biological bursters is examined. Bursting refers to oscillations in which an observable slowly alternates between phases of relative quiescence and rapid oscillatory behavior. The motivation for this work is to understand the role of electrical coupling in promoting the synchronization of bursting electrical activity (BEA) observed in the β-cells of the islet of Langerhans, which secrete insulin in response to glucose. By studying the coupled fast subsystem of a model of BEA, we focus on the interaction that occurs during the rapid oscillatory phase. Coupling is weak, diffusive and non-scalar. In addition, non-identical oscillators are permitted. Using perturbation methods with the assumption that the uncoupled oscillators are near a Hopf bifurcation, a reduced system of equations is obtained. A detailed bifurcation study of this reduced system reveals a variety of patterns but suggests that asymmetrically phase-locked solutions are the most typical. Finally, the results are applied to the unreduced full bursting system and used to predict the burst pattern for a pair of cells with a given coupling strength and degree of heterogeneity. An erratum to this article is available at .     

Resetting Behavior in a Model of Bursting in Secretory Pituitary Cells: Distinguishing Plateaus from Pseudo-Plateaus

We study a recently discovered class of models for plateau bursting, inspired by models for endocrine pituitary cells. In contrast to classical models for fold-homoclinic (square-wave) bursting, the spikes of the active phase are not supported by limit cycles of the frozen fast subsystem, but are transient oscillations generated by unstable limit cycles emanating from a subcritical Hopf bifurcation around a stable steady state. Experimental time courses are suggestive of such fold-subHopf models because the spikes tend to be small and variable in amplitude; we call this pseudo-plateau bursting. We show here that distinct properties of the response to attempted resets from the silent phase to the active phase provide a clearer, qualitative criterion for choosing between the two classes of models. The fold-homoclinic class is characterized by induced active phases that increase towards the duration of the unperturbed active phase as resets are delivered later in the silent phase. For the fold-subHopf class of pseudo-plateau bursting, resetting is difficult and succeeds only in limited windows of the silent phase but, paradoxically, can dramatically exceed the native active phase duration. J.V. Stern and H.M. Osinga contributed equally.     

Irreversible transitions in the 6-phosphofructokinase/fructose 1,6-bisphosphatase cycle.

The dynamics of the fructose 6-phosphate fructose-1,6-bisphosphate cycle operating in an open and homogeneous system reconstituted from purified enzymes was extensively studied. In addition to 6-phosphofructokinase and fructose-1,6-bisphosphatase, pyruvate kinase, adenylate kinae and glucose-6-phosphate isomerase were involved. In that multi-enzyme system, the main source of non-linearity is the reciprocal effect of AMP on the activities of 6-phosphofructokinase and fructose-1,6-bisphosphatase. Depending upon the experimental parameter values, stable attractors, various types of multiple states and sustained oscillations were shown to occur. In the present report we show that irreversible transitions are also likely to occur for realistic operating conditions. Two parameters of the system, that is the adenylate energy charge of the influx and the fructose-1,6-bisphosphatase maximal activity, are potential candidates to provoke such irreversible transitions from one steady state to the other: (a) when varying the maximal activity of fructose-1,6-bisphosphatase, the system can jump irreversibly from a low to a high stable steady state, and (b) when the adenylate energy charge of the influx is the changing parameter, irreversible transitions occur from a high stable steady state to a stable oscillatory state (limit cycle motion). This behavior can be predicted by constructing the loci of limit points and Hopf bifurcation points.     

Modelling the spatio-temporal dynamics of multi-species host-parasitoid interactions: heterogeneous patterns and ecological implications

A mathematical model of the spatio-temporal dynamics of a two host, two parasitoid system is presented. There is a coupling of the four species through parasitism of both hosts by one of the parasitoids. The model comprises a system of four reaction-diffusion equations. The underlying system of ordinary differential equations, modelling the host-parasitoid population dynamics, has a unique positive steady state and is shown to be capable of undergoing Hopf bifurcations, leading to limit cycle kinetics which give rise to oscillatory temporal dynamics. The stability of the positive steady state has a fundamental impact on the spatio-temporal dynamics: stable travelling waves of parasitoid invasion exhibit increasingly irregular periodic travelling wave behaviour when key parameter values are increased beyond their Hopf bifurcation point. These irregular periodic travelling waves give rise to heterogeneous spatio-temporal patterns of host and parasitoid abundance. The generation of heterogeneous patterns has ecological implications and the concepts of temporary host refuge and niche formation are considered.     

Delay-induced model for tumor-immune interaction and control of malignant tumor growth

The paper deals with the qualitative analysis of the solutions of a system of delay differential equations describing the interaction between tumor and immune cells. The asymptotic stability of the possible steady states is showed and the occurrence of limit cycle of the system around the interior equilibrium is proved by the application of Hopf bifurcation theorem by using the delay as a bifurcation parameter. The length of the delay parameter for preserving stability of the system is also estimated, which gives the idea about the mode of action for controlling oscillations in malignant tumor cell growth. The theoretical and numerical outcomes have been supported through experimental results from literatures. This approach gives new insight of modeling tumor-immune interactions and provides significant control strategies to overcome the large oscillations in tumor cells.     

Eye movement instabilities and nystagmus can be predicted by a nonlinear dynamics model of the saccadic system

The study of eye movements and oculomotor disorders has, for four decades, greatly benefitted from the application of control theoretic concepts. This paper is an example of a complementary approach based on the theory of nonlinear dynamical systems. Recently, a nonlinear dynamics model of the saccadic system was developed, comprising a symmetric piecewise-smooth system of six first-order autonomous ordinary differential equations. A preliminary numerical investigation of the model revealed that in addition to generating normal saccades, it could also simulate inaccurate saccades, and the oscillatory instability known as congenital nystagmus (CN). By varying the parameters of the model, several types of CN oscillations were produced, including jerk, bidirectional jerk and pendular nystagmus. The aim of this study was to investigate the bifurcations and attractors of the model, in order to obtain a classification of the simulated oculomotor behaviours. The application of standard stability analysis techniques, together with numerical work, revealed that the equations have a rich bifurcation structure. In addition to Hopf, homoclinic and saddlenode bifurcations organised by a Takens-Bogdanov point, the equations can undergo nonsmooth pitchfork bifurcations and nonsmooth gluing bifurcations. Evidence was also found for the existence of Hopf-initiated canards. The simulated jerk CN waveforms were found to correspond to a pair of post-canard symmetry-related limit cycles, which exist in regions of parameter space where the equations are a slow-fast system. The slow and fast phases of the simulated oscillations were attributed to the geometry of the corresponding slow manifold. The simulated bidirectional jerk and pendular waveforms were attributed to a symmetry invariant limit cycle produced by the gluing of the asymmetric cycles. In contrast to control models of the oculomotor system, the bifurcation analysis places clear restrictions on which kinds of behaviour are likely to be associated with each other in parameter space, enabling predictions to be made regarding the possible changes in the oscillation type that may be observed upon changing the model parameters. The analysis suggests that CN is one of a range of oculomotor disorders associated with a pathological saccadic braking signal, and that jerk and pendular nystagmus are the most probable oscillatory instabilities. Additionally, the transition from jerk CN to bidirectional jerk and pendular nystagmus observed experimentally when the gaze angle or attention level is changed is attributed to a gluing bifurcation. This suggests the possibility of manipulating the waveforms of subjects with jerk CN experimentally to produce waveforms with an extended foveation period, thereby improving visual resolution.     

Dissection of a model for neuronal parabolic bursting

We have obtained new insight into the mechanisms for bursting in a class of theoretical models. We study Plant's model [24] for Aplysia R-15 to illustrate our view of these so-called parabolic bursters, which are characterized by low spike frequency at the beginning and end of a burst. By identifying and analyzing the fast and slow processes we show how they interact mutually to generate spike activity and the slow wave which underlies the burst pattern. Our treatment is essentially the first step of a singular perturbation approach presented from a geometrical viewpoint and carried out numerically with AUTO [12]. We determine the solution sets (steady state and oscillatory) of the fast subsystem with the slow variables treated as parameters. These solutions form the slow manifold over which the slow dynamics then define a burst trajectory. During the silent phase of a burst, the solution trajectory lies approximately on the steady state branch of the slow manifold and during the active phase of spiking, the trajectory sweeps through the oscillation branch. The parabolic nature of bursting arises from the (degenerate) homoclinic transition between the oscillatory branch and the steady state branch. We show that, for some parameter values, the trajectory remains strictly on the steady state branch (to produce a resting steady state or a pure slow wave without spike activity) or strictly in the oscillatory branch (continuous spike activity without silent phases). Plant's model has two slow variables: a calcium conductance and the intracellular free calcium concentration, which activates a potassium conductance. We also show how bursting arises from an alternative mechanism in which calcium inactivates the calcium current and the potassium conductance is insensitive to calcium. These and other biophysical interpretations are discussed.     

Zooplankton mortality and the dynamical behaviour of plankton population models

We investigate the dynamical behaviour of a simple plankton population model, which explicitly simulates the concentrations of nutrient, phytoplankton and zooplankton in the oceanic mixed layer. The model consists of three coupled ordinary differential equations. We use analytical and numerical techniques, focusing on the existence and nature of steady states and unforced oscillations (limit cycles) of the system. The oscillations arise from Hopf bifurcations, which are traced as each parameter in the model is varied across a realistic range. The resulting bifurcation diagrams are compared with those from our previouswork, where zooplankton mortality was simulated by a quadratic function—here we use a linear function, to represent alternative ecological assumptions. Oscillations occur across broader ranges of parameters for the linear mortality function than for the quadratic one, although the two sets of bifurcation diagrams show similar qualitative characteristics. The choice of zooplankton mortality function, or closure term, is an area of current interest in the modelling community, and we relate our results to simulations of other models.     

Asymmetry in the Presence of Migration Stabilizes Multistrain Disease Outbreaks

We study the effect of migration between coupled populations, or patches, on the stability properties of multistrain disease dynamics. The epidemic model used in this work displays a Hopf bifurcation to oscillations in a single, well-mixed population. It is shown numerically that migration between two non-identical patches stabilizes the endemic steady state, delaying the onset of large amplitude outbreaks and reducing the total number of infections. This result is motivated by analyzing generic Hopf bifurcations with different frequencies and with diffusive coupling between them. Stabilization of the steady state is again seen, indicating that our observation in the full multistrain model is based on qualitative characteristics of the dynamics rather than on details of the disease model.     

Impact of fear effect on the growth of prey in a predator-prey interaction model

Several field data and experiments on a terrestrial vertebrates exhibited that the fear of predators would cause a substantial variability of prey demography. Fear for predator population enhances the survival probability of prey population, and it can greatly reduce the reproduction of prey population. Based on the experimental evidence, we proposed and analyzed a prey-predator system introducing the cost of fear into prey reproduction with Holling type-II functional response. We investigate all the biologically feasible equilibrium points, and their stability is analyzed in terms of the model parameters. Our mathematical analysis exhibits that for strong anti-predator responses can stabilize the prey-predator interactions by ignoring the existence of periodic behaviors. Our model system undergoes Hopf bifurcation by considering the birth rate r₀ as a bifurcation parameter. For larger prey birth rate, we investigate the transition to a stable coexisting equilibrium state, with oscillatory approach to this equilibrium state, indicating that the greatest characteristic eigenvalues are actually a pair of imaginary eigenvalues with real part negative, which is increasing for r₀. We obtained the conditions for the occurrence of Hopf bifurcation and conditions governing the direction of Hopf bifurcation, which imply that the prey birth rate will not only influence the occurrence of Hopf bifurcation but also alter the direction of Hopf bifurcation. We identify the parameter regions associated with the extinct equilibria, predator-free equilibria and coexisting equilibria with respect to prey birth rate, predator mortality rates. Fear can stabilize the predator-prey system at an interior steady state, where all the species can exists together, or it can create the oscillatory coexistence of all the populations. We performed some numerical simulations to investigate the relationship between the effects of fear and other biologically related parameters (including growth/decay rate of prey/predator), which exhibit the impact that fear can have in prey-predator system. Our numerical illustrations also demonstrate that the prey become less sensitive to perceive the risk of predation with increasing prey growth rate or increasing predators decay rate.     

Generation of Very Slow Neuronal Rhythms and Chaos Near the Hopf Bifurcation in Single Neuron Models

We have presented a new generation mechanism of slow spiking or repetitive discharges with extraordinarily long inter-spike intervals using the modified Hodgkin-Huxley equations (Doi and Kumagai, 2001). This generation process of slow firing is completely different from that of the well-known potassium A-current in that the steady-state current-voltage relation of the neuronal model is monotonic rather than the N-shaped one of the A-current. In this paper, we extend the previous results and show that the very slow spiking generically appears in both the three-dimensional Hodgkin-Huxley equations and the three dimensional Bonhoeffer-van der Pol (or FitzHugh-Nagumo) equations. The generation of repetitive discharges or the destabilization of the unique equilibrium point (resting potential) is a simple Hopf bifurcation. We also show that the generation of slow spiking does not depend on the stability of the Hopf bifurcation: supercritical or subcritical. The dynamics of slow spiking is investigated in detail and we demonstrate that the phenomenology of slow spiking can be categorized into two types according to the type of the corresponding bifurcation of a fast subsystem: Hopf or saddle-node bifurcation.     

Dynamical analysis of periodic bursting in piece-wise linear planar neuron model

A piece-wise linear planar neuron model, namely, two-dimensional McKean model with periodic drive is investigated in this paper. Periodical bursting phenomenon can be observed in the numerical simulations. By assuming the formal solutions associated with different intervals of this non-autonomous system and introducing the generalized Jacobian matrix at the non-smooth boundaries, the bifurcation mechanism for the bursting solution induced by the slowly varying periodic drive is presented. It is shown that, the discontinuous Hopf bifurcation occurring at the non-smooth boundaries, i.e., the bifurcation taking place at the thresholds of the stimulation, leads the alternation between the rest state and spiking state. That is, different oscillation modes of this non-autonomous system convert periodically due to the non-smoothness of the vector field and the slow variation of the periodic drive as well.     

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.     

Analysis of bifurcations in continuous fermentation using recombinant microorganisms with delayed responses

When the feed rate to a fermenter is varied periodically in order to favor the growth of plasmid-containing cells, a transition may occur from the starting stationary state to another state. The resulting state may be constant or oscillatory. A generalised model based on the adaption times of plasmid-free and plasmid-harboring cells has been used. Analytical conditions have been derived for bifurcation from one nonoscillatory state to another or to an oscillatory state (Hopf bifurcation). The frequency of oscillation is shown to have an upper bound, which can be controlled by manipulating certain process parameters. The production of tryptophan synthetase by the plasmid pPLc23trpAl in E. coli is used as an example to determine the nature of the Hopf bifurcations.     

Method for constructing the boundary of the bursting oscillations region in the neuron model

We examine the problem of constructing the boundary of bursting oscillations on a parameter plane for the system of equations describing the electrical behaviour of the membrane neuron arising from the interaction of fast oscillations of the cytoplasma membrane potential and slow oscillations of the intracellular calcium concentration. As the boundary point on the parameter plane we consider the values at which the limit cycle of the slow subsystem is tangent to the Hopf bifurcation curve of the fast subsystem. The method suggested for determining the boundary is based on the dissection of the system variables into slow and fast. The strong point of the method is that it requires the integration of the slow subsystem only. An example of the application of the method for the stomatogastric neuron model [Guckenheimer J, Gueron S, Harris-Warrick RM (1993) Philos Trans R Soc Lond B 341: 345–359] is given. Received: 31 May 1999 / Accepted in revised form: 19 November 1999     

Synchronization properties of networks of electrically coupled neurons in the presence of noise and heterogeneities

We investigate how synchrony can be generated or induced in networks of electrically coupled integrate-and-fire neurons subject to noisy and heterogeneous inputs. Using analytical tools, we find that in a network under constant external inputs, synchrony can appear via a Hopf bifurcation from the asynchronous state to an oscillatory state. In a homogeneous net work, in the oscillatory state all neurons fire in synchrony, while in a heterogeneous network synchrony is looser, many neurons skipping cycles of the oscillation. If the transmission of action potentials via the electrical synapses is effectively excitatory, the Hopf bifurcation is supercritical, while effectively inhibitory transmission due to pronounced hyperpolarization leads to a subcritical bifurcation. In the latter case, the network exhibits bistability between an asynchronous state and an oscillatory state where all the neurons fire in synchrony. Finally we show that for time-varying external inputs, electrical coupling enhances the synchronization in an asynchronous network via a resonance at the firing-rate frequency.

Fluctuation-driven rhythmogenesis in an excitatory neuronal network with slow adaptation

We study an excitatory all-to-all coupled network of N spiking neurons with synaptically filtered background noise and slow activity-dependent hyperpolarization currents. Such a system exhibits noise-induced burst oscillations over a range of values of the noise strength (variance) and level of cell excitability. Since both of these quantities depend on the rate of background synaptic inputs, we show how noise can provide a mechanism for increasing the robustness of rhythmic bursting and the range of burst frequencies. By exploiting a separation of time scales we also show how the system dynamics can be reduced to low-dimensional mean field equations in the limit N → ∞. Analysis of the bifurcation structure of the mean field equations provides insights into the dynamical mechanisms for initiating and terminating the bursts.