Analysis of an excitable dendritic spine with an activity-dependent stem conductance

 Dendritic spines are the major target for excitatory synaptic inputs in the vertebrate brain. They are tiny evaginations of the dendritic surface consisting of a bulbous head and a tenuous stem. Spines are considered to be an important locus for plastic changes underlying memory and learning processes. The findings that synaptic morphology may be activity-dependent and that spine head membrane may be endowed with voltage-dependent (excitable) channels is the motivation for this study. We first explore the dynamics, when an excitable, yet morphologically fixed spine receives a constant current input. Two parameter Andronov–Hopf bifurcation diagrams are constructed showing stability boundaries between oscillations and steady-states. We show how these boundaries can change as a function of both the spine stem conductance and the conductance load of the attached dendrite. Building on this reference case an idealized model for an activity-dependent spine is formulated and analyzed. Specifically we examine the possibility that the spine stem resistance, the tunable “synaptic weight” parameter identified by Rall and Rinzel, is activity-dependent. In the model the spine stem conductance depends (slowly) on the local electrical interactions between the spine head and the dendritic cable; parameter regimes are found for bursting, steady states, continuous spiking, and more complex oscillatory behavior. We find that conductance load of the dendrite strongly influences the burst pattern as well as other dynamics. When the spine head membrane potential exhibits relaxation oscillations a simple model is formulated that captures the dynamical features of the full model. Received: 10 January 1997/Revised version: 25 March 1997     

Threshold for repetitive activity for a slow stimulus ramp: a memory effect and its dependence on fluctuations

We have obtained new insights into the behavior of a class of excitable systems when a stimulus, or parameter, is slowly tuned through a threshold value. Such systems do not accommodate no matter how slowly a stimulus ramp is applied, and the stimulus value at onset of repetitive activity shows a curious, nonmonotonic dependence on ramp speed. (Jakobsson, E. and R. Guttman. Biophys. J. 1980. 31:293-298.) demonstrated this for squid axon and for the Hodgkin-Huxley (HH) model. Furthermore, they showed theoretically that for moderately slow ramps the threshold increases as the ramp speed decreases, but for much slower ramp speeds threshold decreases as the ramp speed decreases. This latter feature was found surprising and it was suggested that the HH model, and squid axon in low calcium, exhibits reverse accommodation. We have found that reverse accommodation reflects the influence of persistent random fluctuations, and is a feature of all such excitable systems. We have derived an analytic condition which yields an approximation for threshold in the case of a slow ramp when the effect of fluctuations are negligible. This condition predicts, and numerical calculations confirm, that the onset of oscillations occurs beyond the critical stimulus value which is predicted by treating the stimulus intensity as a static parameter, i.e., the dynamic aspect of a ramp leads to a delay in the onset. The condition further demonstrates a memory effect, i.e., firing threshold is dependent on the initial state of the system. For very slow ramps then, fluctuations diminish both the delay and memory effects.(ABSTRACT TRUNCATED AT 250 WORDS)     

Hopf bifurcation and bursting synchronization in an excitable systems with chemical delayed coupling

In this paper we consider the Hopf bifurcation and synchronization in the two coupled Hindmarsh–Rose excitable systems with chemical coupling and time-delay. We surveyed the conditions for Hopf bifurcations by means of dynamical bifurcation analysis and numerical simulation. The results show that the coupled excitable systems with no delay have supercritical Hopf bifurcation, while the delayed system undergoes Hopf bifurcations at critical time delays when coupling strength lies in a particular region. We also investigated the effect of the delay on the transition of bursting synchronization in the coupled system. The results are helpful for us to better understand the dynamical properties of excitable systems and the biological mechanism of information encoding and cognitive activity.     

Spontaneous secondary spiking in excitable cells

Kepler & Marder (1993, Biol. Cybern.68, 209-214) proposed a model describing the electrical activity of a crab neuron in which a train of directly induced action potentials is sometimes followed by one or more spontaneous action potentials, referred to as spontaneous secondary spikes. We reduce their five-dimensional model to three dimensions in two different ways in order to gain insight into the mechanism underlying the spontaneous spikes. We then treat a slowly varying current as a parameter in order to give a qualitative explanation of the phenomenon using phase-plane and bifurcation analysis. We demonstrate that a three-dimensional model, consisting of a two-dimensional excitable system plus a slow inward current, is sufficient to produce the behaviour observed in the original model. The exact dynamics of the excitable system are not important, but the relative time constant and amplitude of the slow inward current are crucial. Using the numerical bifurcation analysis package AUTO (Doedel & Kernevez, 1986, AUTO: Software for Continuation and Bifurcation Problems in Ordinary Differential Equations. California Institute of Technology), we compute bifurcation diagrams using the maximum amplitude of the slow inward current as the bifurcation parameter. The full and reduced models have a stable resting potential for all values of the bifurcation parameter. At a critical value of the bifurcation parameter, a stable tonic firing mode arises via a saddle-node of periodics bifurcation. Whether or not the models can exhibit transient or continuous spontaneous spiking depends on their position in parameter space relative to this saddle-node of periodics.     

Understanding bursting oscillations as periodic slow passages through bifurcation and limit points

We consider a biochemical system consisting of two allosteric enzyme reactions coupled in series. The system has been modeled by Decroly and Goldbeter (J. Theor. Biol. 124, 219 (1987)) and is described by three coupled, first-order, nonlinear, differential equations. Bursting oscillations correspond to a succession of alternating active and silent phases. The active phase is characterized by rapid oscillations while the silent phase is a period of quiescence. We propose an asymptotic analysis of the differential equations which is based on the limit of large allosteric constants. This analysis allows us to construct a time-periodic bursting solution. This solution is jumping periodically between a slowly varying steady state and a slowly varying oscillatory state. Each jump follows a slow passage through a bifurcation or limit point which we analyze in detail. Of particular interest is the slow passage through a supercritical Hopf bifurcation. The transition is from an oscillatory solution to a steady state solution. We show that the transition is delayed considerably and characterize this delay by estimating the amplitude of the oscillations at the Hopf bifurcation point.     

A delay reaction-diffusion model of the dynamics of botulinum in fish

A model of interaction between fish and a bacterium (Clostridium botulinum) responsible for avian botulism is introduced, considering diffusion of both fish and bacterium in water. The fish population moves randomly in water. Death fish disintegrate in water, at different locations, causing bacteria to diffuse through water and infect other fish. Existence of uniform steady states is investigated and the linearized stability of the positive uniform steady state is analyzed. A Hopf bifurcation is proved to occur from the uniform steady state when the bifurcation parameter, here the time delay, passes through a critical value and diffusion coefficients satisfy some conditions, that induces time oscillations of the populations. Comments on diffusion-driven instability are provided, and numerical simulations are carried out to illustrate the results.     

Ionic channel density of excitable membranes can act as a bifurcation parameter

As the maximal K⁺-conductance (or K⁺-channel density) of the Hodgkin-Huxley equations is reduced, the stable resting membrane potential bifurcates at a subcritical Hopf bifurcation into small amplitude unstable oscillations. These small amplitude solutions jump to large amplitude periodic solutions that correspond to a repetitive discharge of action potentials. Thus the specific channel density can act as a bifurcation parameter, and can control the excitability and autorhythmicity of excitable membranes.     

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.     

Bursting Regimes in a Reaction-Diffusion System with Action Potential-Dependent Equilibrium

The equilibrium Nernst potential plays a critical role in neural cell dynamics. A common approximation used in studying electrical dynamics of excitable cells is that the ionic concentrations inside and outside the cell membranes act as charge reservoirs and remain effectively constant during excitation events. Research into brain electrical activity suggests that relaxing this assumption may provide a better understanding of normal and pathophysiological functioning of the brain. In this paper we explore time-dependent ionic concentrations by allowing the ion-specific Nernst potentials to vary with developing transmembrane potential. As a specific implementation, we incorporate the potential-dependent Nernst shift into a one-dimensional Morris-Lecar reaction-diffusion model. Our main findings result from a region in parameter space where self-sustaining oscillations occur without external forcing. Studying the system close to the bifurcation boundary, we explore the vulnerability of the system with respect to external stimulations which disrupt these oscillations and send the system to a stable equilibrium. We also present results for an extended, one-dimensional cable of excitable tissue tuned to this parameter regime and stimulated, giving rise to complex spatiotemporal pattern formation. Potential applications to the emergence of neuronal bursting in similar two-variable systems and to pathophysiological seizure-like activity are discussed.     

Roles of L-type Ca2+ and delayed-rectifier K+ currents in sinoatrial node pacemaking: insights from stability and bifurcation analyses of a mathematical model

To elucidate the dynamical mechanisms of the sinoatrial (SA) node pacemaker activity, we investigated the roles of L-type Ca2+ (ICa,L) and delayed-rectifier K+ (IKr) currents in pacemaking by stability and bifurcation analyses of our rabbit SA node model (Kurata Y, Hisatome I, Imanishi S, and Shibamoto T. Am J Physiol Heart Circ Physiol 283: H2074-H2101, 2002). Equilibrium points (EPs), periodic orbits, stability of EPs, and Hopf bifurcation points were calculated as functions of conductance or gating time constants of the currents for constructing bifurcation diagrams. Structural stability (robustness) of the system was also evaluated by computing stability and dynamics during applications of constant bias currents (Ibias). Blocking ICa,L or IKr caused stabilization of an EP and cessation of pacemaking via a Hopf bifurcation. The unstable zero-current potential region determined with Ibias applications, where spontaneous oscillations appear, shrunk and finally disappeared as ICa,L diminished, but shrunk little when IKr was eliminated. The reduced system, including no time-dependent current except ICa,L, exhibited pacemaker activity. These results suggest that ICa,L is responsible for EP instability and pacemaker generation, whereas IKr is not necessarily required for constructing a pacemaker cell system. We further explored the effects of various K+ currents with different kinetics on stability and dynamics of the model cell. The original IKr of delayed activation and inward rectification appeared to be most favorable for generating large-amplitude oscillations with stable frequency, suggesting that IKr acts as an oscillation amplifier and frequency stabilizer. IKr may also play an important role in preventing bifurcation to quiescence of the system.     

Global Hopf bifurcation analysis of an susceptible-infective-removed epidemic model incorporating media coverage with time delay

An susceptible-infective-removed epidemic model incorporating media coverage with time delay is proposed. The stability of the disease-free equilibrium and endemic equilibrium is studied. And then, the conditions which guarantee the existence of local Hopf bifurcation are given. Furthermore, we show that the local Hopf bifurcation implies the global Hopf bifurcation after the second critical value of delay. The obtained results show that the time delay in media coverage can not affect the stability of the disease-free equilibrium when the basic reproduction number is less than unity. However, the time delay affects the stability of the endemic equilibrium and produces limit cycle oscillations while the basic reproduction number is greater than unity. Finally, some examples for numerical simulations are included to support the theoretical prediction.     

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.     

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     

Hopf bifurcation in three-species food chain models with group defense.

Three-species food-chain models, in which the prey population exhibits group defense, are considered. Using the carrying capacity of the environment as the bifurcation parameter, it is shown that the model without delay undergoes a sequence of Hopf bifurcations. In the model with delay it is shown that using a delay as a bifurcation parameter, a Hopf bifurcation can also occur in this case. These occurrences may be interpreted as showing that a region of local stability (survival) may exist even though the positive steady states are unstable. A computer code BIFDD is used to determine the stability of the bifurcation solutions of a delay model.     

Stability of periodic solutions for a model of genetic repression with delays

A technique is discussed for locating the Hopf bifurcation of an n-dimensional system of delay differential equations which arises from a model for control of protein biosynthesis. Certain parameter values are shown to allow a Hopf bifurcation to periodic orbits. At the Hopf bifurcation the periodic orbits are shown to be stable either analytically or numerically depending on the parameter values.On leave from North Carolina State University.Supported in part by N.S.F. Grant # MCS 81-02828     

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.     

Deconstructing the role of myosin contractility in force fluctuations within focal adhesions

Force fluctuations exhibited in focal adhesions that connect a cell to its extracellular environment point to the complex role of the underlying machinery that controls cell migration. To elucidate the explicit role of myosin motors in the temporal traction force oscillations, we vary the contractility of these motors in a dynamical model based on the molecular clutch hypothesis. As the contractility is lowered, effected both by changing the motor velocity and the rate of attachment/detachment, we show analytically in an experimentally relevant parameter space, that the system goes from decaying oscillations to stable limit cycle oscillations through a supercritical Hopf bifurcation. As a function of the motor activity and the number of clutches, the system exhibits a rich array of dynamical states. We corroborate our analytical results with stochastic simulations of the motor-clutch system. We obtain limit cycle oscillations in the parameter regime as predicted by our model. The frequency range of oscillations in the average clutch and motor deformation compares well with experimental results.     

External Noise and External Signal Induced Transition of Gene Switch and Coherence Resonance in the Genetic Regulatory System

The transition of gene switch induced by external noises (multiplicative external noise and additive external noise) and external signals is investigated in the genetic regulatory system. Results show that the state-to-state transition of gene switch as well as resonant behaviors, such as the explicit coherence resonance (ECR), implicit coherence resonance (ICR) and control parameter coherence biresonance (CPCBR), can appear when noises are injected into the genetic regulatory system. The ECR is increased with the increase of the control parameter value when starting from the supercritical Hopf bifurcation parameter point, and there exists a critical control parameter value for the occurrence of ECR. However, the ICR is decreased as the control parameter value is increased when starting from the subcritical Hopf bifurcation point. In particular, the coherence of ECR is higher and more sensitive to noise than that of ICR. When an external signal is introduced into the system, the enhancement or suppression of the CPCBR and the number of peaks strongly depend on the frequency and amplitude of the external signal. Furthermore, the gene regulation system can selectively enhance or decrease the noise-induced oscillation signals at preferred frequency and amplitude of an external signal.