Critical slowing down as an indicator of transitions in two-species models

Transitions in ecological systems often occur without apparent warning, and may represent shifts between alternative persistent states. Decreasing ecological resilience (the size of the basin of attraction around a stable state) can signal an impending transition, but this effect is difficult to measure in practice. Recent research has suggested that a decreasing rate of recovery from small perturbations (critical slowing down) is a good indicator of ecological resilience. Here we use analytical techniques to draw general conclusions about the conditions under which critical slowing down provides an early indicator of transitions in two-species predator-prey and competition models. The models exhibit three types of transition: the predator-prey model has a Hopf bifurcation and a transcritical bifurcation, and the competition model has two saddle-node bifurcations (in which case the system exhibits hysteresis) or two transcritical bifurcations, depending on the parameterisation. We find that critical slowing down is an earlier indicator of the Hopf bifurcation in predator-prey models in which prey are regulated by predation rather than by intrinsic density-dependent effects and an earlier indicator of transitions in competition models in which the dynamics of the rare species operate on slower timescales than the dynamics of the common species. These results lead directly to predictions for more complex multi-species systems, which can be tested using simulation models or real ecosystems.     

Stability switches and double Hopf bifurcation in a two-neural network system with multiple delays

Time delay is an inevitable factor in neural networks due to the finite propagation velocity and switching speed. Neural system may lose its stability even for very small delay. In this paper, a two-neural network system with the different types of delays involved in self- and neighbor- connection has been investigated. The local asymptotic stability of the equilibrium point is studied by analyzing the corresponding characteristic equation. It is found that the multiple delays can lead the system dynamic behavior to exhibit stability switches. The delay-dependent stability regions are illustrated in the delay-parameter plane, followed which the double Hopf bifurcation points can be obtained from the intersection points of the first and second Hopf bifurcation, i.e., the corresponding characteristic equation has two pairs of imaginary eigenvalues. Taking the delays as the bifurcation parameters, the classification and bifurcation sets are obtained in terms of the central manifold reduction and normal form method. The dynamical behavior of system may exhibit the quasi-periodic solutions due to the Neimark- Sacker bifurcation. Finally, numerical simulations are made to verify the theoretical results.     

由三态跃迁机制产生的两种神经放电节律

在实验性神经起步点自发放电中,发现了两种三态跃迁节律,其特征为静息、周期n及周期n+1(n=1,2)簇放电随机交替出现。应用随机Chay模型数值仿真,分别得到了与实验模型中相似的两种三态跃迁节律,这两种节律都是在两个紧邻的分岔点附近,由噪声驱动而产生的。理论分析提示,当神经元系统接近从静息经分岔到放电的临界状态,且从静息到周期n的分岔点,与从周期n到周期n+1的分岔点非常接近时,在噪声的作用下,系统运动会在静息、周期n和周期n+1三种状态之间随机跃迁,从而形成了这种三态跃迁节律。基于这种三态跃迁放电的随机共振,还有待进一步深入研究。     

Hysteresis,oscillations, and pattern formation in realistic immobilized enzyme systems

Summary Hysteresis, oscillations, and pattern formation in realistic biochemical systems governed by P.D.E.s are considered from both numerical and mathematical points of view. Analysis of multiple steady states in the case of hysteresis, and bifurcation theory in the cases of oscillations and pattern formation, account for the observed numerical results. The possibility to realize these systems experimentally is their main interest, thus bringing further arguments in favor of theories explaining basic biological phenomena by diffusion and reaction.     

Computational multiple steady states for enzymatic esterification of ethanol and oleic acid in an isothermal CSTR

The capacity of complex biochemical reaction networks (consisting of 11 coupled non-linear ordinary differential equations) to show multiple steady states, was investigated. The system involved esterification of ethanol and oleic acid by lipase in an isothermal continuous stirred tank reactor (CSTR). The Deficiency One Algorithm and the Subnetwork Analysis were applied to determine the steady state multiplicity. A set of rate constants and two corresponding steady states are computed. The phenomena of bistability, hysteresis and bifurcation are discussed. Moreover, the capacity of steady state multiplicity is extended to the family of the studied reaction networks.     

The influence of depolarization block on seizure-like activity in networks of excitatory and inhibitory neurons

The inhibitory restraint necessary to suppress aberrant activity can fail when inhibitory neurons cease to generate action potentials as they enter depolarization block. We investigate possible bifurcation structures that arise at the onset of seizure-like activity resulting from depolarization block in inhibitory neurons. Networks of conductance-based excitatory and inhibitory neurons are simulated to characterize different types of transitions to the seizure state, and a mean field model is developed to verify the generality of the observed phenomena of excitatory-inhibitory dynamics. Specifically, the inhibitory population’s activation function in the Wilson-Cowan model is modified to be non-monotonic to reflect that inhibitory neurons enter depolarization block given strong input. We find that a physiological state and a seizure state can coexist, where the seizure state is characterized by high excitatory and low inhibitory firing rate. Bifurcation analysis of the mean field model reveals that a transition to the seizure state may occur via a saddle-node bifurcation or a homoclinic bifurcation. We explain the hysteresis observed in network simulations using these two bifurcation types. We also demonstrate that extracellular potassium concentration affects the depolarization block threshold; the consequent changes in bifurcation structure enable the network to produce the tonic to clonic phase transition observed in biological epileptic networks.     

Computational Multiple Steady States of the Penicillin G Acylase‐Catalyzed Hydrolysis in an Isothermal CFSTR

A family of an enzymatically catalyzed reaction network was studied, which involves the hydrolysis of penicillin G by penicillin G acylase in an isothermal continuous flow stirred tank reactor (CFSTR). This system consisted of 10 coupled non‐linear equations and was found to be capable of exhibiting computational multiple steady states. A set of kinetic parameters determined from the existing experimental data were used to compute a set of rate constants and two corresponding steady states. This suggested that multiple steady states may occur in the system studied. The phenomena of bistability, hysteresis and bifurcation were discussed. Moreover, the capacity of steady state multiplicity was extended to its family of reaction networks.     

Catastrophic shifts in vertical distributions of phytoplankton The existence of a bifurcation set

A model of phytoplankton dynamics within a water column was analyzed with special consideration on the existence of a bifurcation set in the parameter space. We considered two resources, light and a limiting nutrient, for phytoplankton growth and assumed that the water column is separated into two layers by thermal and/or density stratification. It was shown that there exists a bifurcation set in the parameter space when the growth function meets several conditions that are general for growth functions of two essential resources. Specifically, these conditions include that a less abundant of the two resources limits the growth while the effect of the other is sufficiently small. Folded structure with two stable states separated by one unstable state appears in the catastrophe manifold when parameters move to a certain direction with a certain curvature from a point in the bifurcation set. These results suggest that occurrence of discontinuous transition between two alternative vertical patterns is possible nature of phytoplankton dynamics within a stratified water column.     

心肌细胞团搏动整数倍节律的非线性动力学机制

利用心肌细胞耦合模型研究心肌整数倍节律的动力学机理。确定性模型仿真揭示了心肌细胞团同步搏动加周期分岔的节律变化规律;随机模型仿真发现在加周期分岔序列中分岔点附近会出现整数倍节律,其中,0-1整数倍节律产生于从静息到周期1的Hopf分岔点附近,1-2整数倍节律产生于周期1和周期2极限环间的加周期分岔点附近;对系统相空间轨道的分析进一步揭示出整数倍节律是由系统运动在相邻的两个轨道之间随机跃迁形成的。上述分析结果不仅阐明了心肌整数倍节律的机理,并且揭示了各种整数倍节律与加周期分岔序列中相邻节律的内在联系,为重新认识心律变化的规律开辟了新的途径。     

From bistability to oscillations in a model for the isocitrate dehydrogenase reaction

Considered is a bienzymatic system consisting of isocitrate dehydrogenase (IDH, EC 1.1.1.42), which transforms NADP(+) into NADPH, and of diaphorase (DIA, EC 1.8.1.4), which catalyzes the reverse reaction. Experimental evidence as well as a theoretical model show the possibility of a coexistence between two stable steady states in this reaction system. The phenomenon originates from the regulatory properties of IDH. We extend the analysis of a theoretical model proposed for the IDH-DIA bienzymatic system and investigate the occurrence of different modes of bistability, with or without hysteresis, i.e. in the presence of two or only one limit point bounding the domain of multiple steady states. The analysis indicates that the two types of bistability may sometimes be observed sequentially as a given control parameter is progressively increased. We further obtain conditions in which sustained oscillations develop in the model. These results establish the isocitrate dehydrogenase reaction coupled to diaphorase as a suitable candidate for further experimental and theoretical studies of bistability and oscillations in biochemical systems.     

Critical Slowing Down Governs the Transition to Neuron Spiking

Many complex systems have been found to exhibit critical transitions, or so-called tipping points, which are sudden changes to a qualitatively different system state. These changes can profoundly impact the functioning of a system ranging from controlled state switching to a catastrophic break-down; signals that predict critical transitions are therefore highly desirable. To this end, research efforts have focused on utilizing qualitative changes in markers related to a system’s tendency to recover more slowly from a perturbation the closer it gets to the transition—a phenomenon called critical slowing down. The recently studied scaling of critical slowing down offers a refined path to understand critical transitions: to identify the transition mechanism and improve transition prediction using scaling laws.Here, we outline and apply this strategy for the first time in a real-world system by studying the transition to spiking in neurons of the mammalian cortex. The dynamical system approach has identified two robust mechanisms for the transition from subthreshold activity to spiking, saddle-node and Hopf bifurcation. Although theory provides precise predictions on signatures of critical slowing down near the bifurcation to spiking, quantitative experimental evidence has been lacking. Using whole-cell patch-clamp recordings from pyramidal neurons and fast-spiking interneurons, we show that 1) the transition to spiking dynamically corresponds to a critical transition exhibiting slowing down, 2) the scaling laws suggest a saddle-node bifurcation governing slowing down, and 3) these precise scaling laws can be used to predict the bifurcation point from a limited window of observation. To our knowledge this is the first report of scaling laws of critical slowing down in an experiment. They present a missing link for a broad class of neuroscience modeling and suggest improved estimation of tipping points by incorporating scaling laws of critical slowing down as a strategy applicable to other complex systems.     

The influence of spermine on the structural dynamics of yeast tRNAPhe

A tRNAPhe derivative carrying ethidium at position 37 in the anticodon loop has been used to study the effect of spermine on conformational transitions of the tRNA. As previously reported (Ehrenberg, M., Rigler, R. and Wintermeyer, W. (1979) Biochemistry 18, 4588-4599) in the tRNA derivative the ethidium is present in three states (T1-T3) characterized by different fluorescence decay rates. T-jump experiments show two transitions between the states, a fast one (relaxation time 10-100 ms) between T1 and T2, and a slow one (100-1000 ms) between T2 and T3. In the presence of spermine the fast transition shows a negative temperature coefficient indicating the existence of a preequilibrium with a negative reaction enthalpy. Spermine shifts the distribution of states towards T3, as does Mg2+, but the final ratio [T2]/[T1] obtained with spermine is higher than with Mg2+, which we tentatively interpret to mean that spermine stabilizes one particular conformation of the anticodon loop.     

Two stable steady states in the Hodgkin-Huxley axons.

Two stable steady states were found in the numerical solution of the Hodgkin-Huxley equations for the intact squid axon bathed in potassium-rich sea water with an externally applied inward current. Under the conditions the two stable steady-states exist, the Hodgkin-Huxley equations have a complex bifurcation structure including, in addition to the two stable steady-states, a stable limit cycle, two unstable equilibrium points, and one asymptotically stable equilibrium point. It was also concluded that two stable steady states can appear in the Hodgkin-Huxley axons when the leak current is comparable to the currents through the Na and K channels.     

Method of determining the relative stability of different conformational states of biological macromolecules

A general approach to the determination of relative stability of any pair of con-formational states of biological macromolecules or their complexes (in particular, to the determination of relative stability of native and disordered states of the macromolecule) has been suggested. For determining the free energy difference of the two states under the conditions when one of them is considerably more advantageous than the other, it is necessary for the macromolecule to be influenced by the transforming agent which levels free energies of both the conformational states, and to determine the external parameter derivative of the free energy difference in the region of the conformational transition induced by the change in this parameter. If the character of the dependence of this derivative on the external parameter (temperature, solvent composition, etc.) is known, then this allows the determination of the free energy difference of the two states under the conditions considered, even including conditions far from the transition region. The value of the derivative of the free energy difference in the transition region in many cases can be measured directly (in particular, when using calorimetry), while in cases when a direct measurement of the derivative is impossible, it can often he estimated experimentally from the steepness of the conformational transition. The methods of this estimation and also a possible character of the change of the considered derivative during variation of the external parameter are considered for the case when the transforming agent is one of the components of the solvent and, consequently, the derivative of the free energy difference is equal to the difference of number of molecules of this component hound with the macromolecule in two conformational states.     

Dynamics and bifurcations of a discrete-time prey-predator model with Allee effect on the prey population

This paper studies the dynamic behavior of a discrete-time prey-predator model. It is shown that this model undergoes codimension one and codimension two bifurcations such as transcritical, flip (period-doubling), Neimark-Sacker and strong resonances 1:2, 1:3 and 1:4. The bifurcation analysis is based on the numerical normal form method and the bifurcation scenario around the bifurcation point is determined by their critical normal form coefficients. The advantage of this method is that there is no need to calculate the center manifold and to convert the linear part of the map to a Jordan form. The bifurcation curves of fixed points under variation of one and two parameters are obtained, and the codimensions one and the two bifurcations on the corresponding curves are computed.     

Equilibrium analysis and phase synchronization of two coupled HR neurons with gap junction

The properties of equilibria and phase synchronization involving burst synchronization and spike synchronization of two electrically coupled HR neurons are studied in this paper. The findings reveal that in the non-delayed system the existence of equilibria can be turned into intersection of two odd functions, and two types of equilibria with symmetry and non-symmetry can be found. With the stability and bifurcation analysis, the bifurcations of equilibria are investigated. For the delayed system, the equilibria remain unchanged. However, the Hopf bifurcation point is drastically affected by time delay. For the phase synchronization, we focus on the synchronization transition from burst synchronization to spike synchronization in the non-delayed system and the effect of coupling strength and time delay on spike synchronization in delayed system. In addition, corresponding firing rhythms and spike synchronized regions are obtained in the two parameters plane. The results allow us to better understand the properties of equilibria, multi-time-scale properties of synchronization and temporal encoding scheme in neuronal systems.     

Electronic structure and decomposition reaction mechanism of cyclopropenone, phenylcylopropenone and their sulfur analogues: a theoretical study

The electronic structure, the origin of the extraordinary stability and the reaction mechanisms of the decomposition reaction of the three-membered ring cyclopropenone (IO), its phenyl derivative (IIO) and its sulfur analogues (IS and IIS) have been investigated at the B3LYP/6-311?+?G** level of theory. All critical points on the reaction surface, reactants, transition states and intermediates were determined. Reaction rate constants and half-lives have been computed. Natural bond orbital (NBO) analysis has been used to investigate the type and extent of interaction in the studied species. Results indicate that the decomposition reaction occurs via a stepwise mechanism, with the formation of a short-lived intermediate. The characters of the intermediates for the decomposition of IIO and IIS are different. In case of IIO decomposition, the intermediate structure is of prevailing zwitterionic character, whereas that for the decomposition of IIS is of prevailing carbene character. Solvent effects are computed, analyzed and discussed.     

Quantification of the non-planarity of the human carotid bifurcation

The carotid bifurcation has been a region of particular interest due to its predilection for clinically significant atherosclerosis. It has been shown that the vessel geometry is a major determinant of the local haemodynamic properties which are believed to be associated with the location of atherosclerotic lesions. Current knowledge of the geometry of the carotid bifurcation is insufficient and restricted to basic geometric parameters. To provide some means of quantifying the degree of complexity of the 3D shape of the bifurcation, we made an initial attempt by evaluating the non-planarity of an arterial bifurcation based upon the singular value decomposition theorem.In this paper we present our results obtained on the right carotid bifurcations of six normal subjects, each of whom was scanned twice using the 2D time-of-flight MR sequence. The acquired 2D cross sectional images were processed by using our in-house software which comprises 2D segmentation, 3D reconstruction and smoothing. The centroids of each transverse slices were determined and used as input data for the non-planarity analysis. Our results using the singular value decomposition method have demonstrated discernible differences in non-planarity among individuals. Comparisons with the planarity definition proposed by other investigators suggest that the singular value decomposition method offers more information about the linearity and planarity of the bifurcation. However, it is also realised that a single measure of non-planarity can never fully characterise a bifurcation owing to the great variety of geometries.     

Effect of hunting cooperation and fear in a predator-prey model

Dynamics of predator-prey systems under the influence of cooperative hunting among predators and the fear thus imposed on the prey population is of great importance from ecological point of view. The role of hunting cooperation and the fear effect in the predator-prey system is gaining considerable attention by the researchers recently. But the study on combined effect of hunting cooperation and fear in the predator-prey system is not yet studied. In the present paper, we investigate the impact of hunting cooperation among predators and predator induced fear in prey population by using the classical predator-prey model. We consider that predator populations cooperate during hunting. We also consider that hunting cooperation induces fear among prey, which has far richer and complex dynamics. We observe that without hunting cooperation, the unique coexistence equilibrium point is globally asymptotically stable. However, an increase in the hunting cooperation induced fear may destabilize the system and produce periodic solution via Hopf-bifurcation. The stability of the Hopf-bifurcating periodic solution is obtained by computing the Lyapunov coefficient. The limit cycles thus obtained may be supercritical or subcritical. We also observe that the system undergoes the Bogdanov-Takens bifurcation in two-parameter space. Further, we observe that the system exhibits backward bifurcation between predator-free equilibrium and coexisting equilibrium. The system also exhibits two different types of bi-stabilities due to subcritical Hopf-bifurcation (between interior equilibrium and stable limit cycle) and backward bifurcation (between predator-free and interior equilibrium points). Further, we observe strong demographic Allee phenomenon in the system. To visualize the dynamical behavior of the system, extensive numerical experiments are performed by using MATLAB and MATCONT softwares.     

Metastability in a temperature-dependent model system for predator-prey mite outbreak interactions on fruit trees

The nonlinear behavior of a particular Kolmogorov-type exploitation differential equation system assembled by May (1973,Stability and Complexity in Model Ecosystems, Princeton University Press) from predator and prey components developed by Leslie (1948,Biometrica 35, 213–245) and Holling (1973,Mem. Entomol. Soc. Can. 45, 1–60), respectively, is re-examined by means of the numerical bifurcation code AUTO 86 with model parameters chosen appropriately for a temperature dependent mite interaction on fruit trees. The most significant result of this analysis is that, in addition to the temperature ranges over which the single community equilibrium point of the system iseither globally stableor gives rise to a globally stable limit cycle, there can also exist a range wherein multiple stable states occur. These stable states consist of a focus (spiral point) and a limit cycle, separated from each other in the phase plane by an unstable limit cycle. The ecological implications of such metastability, hysteresis and threshold behavior for the occurrence of outbreaks, the persistence of oscillations, the resiliency of the system and the biological control of mite populations are discussed. It is further suggested that a model of this sort which possesses a single community equilibrium point may be more useful for representing outbreak phenomena, especially in the presence of oscillations, than the non-Kolmogorov predator-prey systems possessing three community equilibrium points, two of which are stable and the other a saddle point, traditionally employed for this purpose.