Theoretical prediction of early afterdepolarization-evoked triggered activity formation initiating ventricular reentrant arrhythmias.

BACKGROUND AND OBJECTIVE: Excessive prolongation of QT interval on ECGs in patients with congenital/acquired long QT syndrome and heart failure is a sign suggesting the development of early afterdepolarization (EAD), an abnormal repolarization in the action potential of ventricular cardiomyocytes. The development of EAD has been believed to be a trigger for fatal tachyarrhythmia, which can be a risk for sudden cardiac death. The role of EAD in triggering ventricular tachycardia (VT) remains unclear. The aim of this study was to elucidate the mechanism of EAD-induced triggered activity formation that leads to the VT such as Torsades de Pointes. METHODS: We investigated the relationship between EAD and tachyarrhythmia initiation by constructing homogeneous myocardial sheet models consisting of the mid-myocardial cell version of a human ventricular myocyte model and performing simulations of excitation propagation. RESULTS: A solitary island-like (clustering) occurrence of EADs in the homogeneous myocardial sheet could induce a focal excitation wave. However, reentrant excitation, an entity of tachyarrhythmia, was not able to be triggered regardless of the EAD cluster size when the focal excitation wave formed a repolarization potential difference boundary consisting of only a convex surface. The discontinuous distribution of multiple EAD clusters in the ventricular tissue formed a specific repolarization heterogeneity due to the repolarization potential difference, the shape of which depended on EAD cluster size and placed intervals. We found that the triggered activity was formed in such a manner that the repolarization potential difference boundary included a concave surface. CONCLUSIONS: The formation of triggered activity that led to tachyarrhythmia required not only the occurrence of EAD onset-mediated focal excitation wave but also a repolarization heterogeneity-based specific repolarization potential difference boundary shape formed within the tissue.

Bifurcations and Proarrhythmic Behaviors in Cardiac Electrical Excitations.

The heart is a hierarchical dynamic system consisting of molecules, cells, and tissues, and acts as a pump for blood circulation. The pumping function depends critically on the preceding electrical activity, and disturbances in the pattern of excitation propagation lead to cardiac arrhythmia and pump failure. Excitation phenomena in cardiomyocytes have been modeled as a nonlinear dynamical system. Because of the nonlinearity of excitation phenomena, the system dynamics could be complex, and various analyses have been performed to understand the complex dynamics. Understanding the mechanisms underlying proarrhythmic responses in the heart is crucial for developing new ways to prevent and control cardiac arrhythmias and resulting contractile dysfunction. When the heart changes to a pathological state over time, the action potential (AP) in cardiomyocytes may also change to a different state in shape and duration, often undergoing a qualitative change in behavior. Such a dynamic change is called bifurcation. In this review, we first summarize the contribution of ion channels and transporters to AP formation and our knowledge of ion-transport molecules, then briefly describe bifurcation theory for nonlinear dynamical systems, and finally detail its recent progress, focusing on the research that attempts to understand the developing mechanisms of abnormal excitations in cardiomyocytes from the perspective of bifurcation phenomena.

Specific decreasing of Na+ channel expression on the lateral membrane of cardiomyocytes causes fatal arrhythmias in Brugada syndrome.

Reduced cardiac sodium (Na+) channel current (INa) resulting from the loss-of-function of Na+ channel is a major cause of lethal arrhythmias in Brugada syndrome (BrS). Inspired by previous experimental studies which showed that in heart diseases INa was reduced along with expression changes in Na+ channel within myocytes, we hypothesized that the local decrease in INa caused by the alteration in Na+ channel expression in myocytes leads to the occurrence of phase-2 reentry, the major triggering mechanism of lethal arrhythmias in BrS. We constructed in silico human ventricular myocardial strand and ring models, and examined whether the Na+ channel expression changes in each myocyte cause the phase-2 reentry in BrS. Reducing Na+ channel expression in the lateral membrane of each myocyte caused not only the notch-and-dome but also loss-of-dome type action potentials and slowed conduction, both of which are typically observed in BrS patients. Furthermore, the selective reduction in Na+ channels on the lateral membrane of each myocyte together with spatial tissue heterogeneity of Na+ channel expression caused the phase-2 reentry and phase-2 reentry-mediated reentrant arrhythmias. Our data suggest that the BrS phenotype is strongly influenced by expression abnormalities as well as genetic abnormalities of Na+ channels.

Integrative and theoretical research on the architecture of a biological system and its disorder.

An organism stems from assemblies of a variety of cells and proteins. This complex system serves as a unit, and it exhibits highly sophisticated functions in response to exogenous stimuli that change over time. The complete sequencing of the entire human genome has allowed researchers to address the enigmas of life and disease at the gene- or molecular-based level. The consequence of such studies is the rapid accumulation of a multitude of data at multiple levels, ranging from molecules to the whole body, that has necessitated the development of entirely new concepts, tools, and methodologies to analyze and integrate these data. This necessity has given birth to systems biology, an advanced theoretical and practical research framework that has totally changed the directions of not only basic life science but also medicine. During the symposium of the 95th Annual Meeting of The Physiological Society of Japan 2018, five researchers reported on their respective studies on systems biology. The topics included reactions of drugs, ion-transport architecture in an epithelial system, multi-omics in renal disease, cardiac electrophysiological systems, and a software platform for computer simulation. In this review article these authors have summarized recent achievements in the field and discuss next-generation studies on health and disease.

Facilitation of I Kr current by some hERG channel blockers suppresses early afterdepolarizations.

Drug-induced block of the cardiac rapid delayed rectifying potassium current (I Kr), carried by the human ether-a-go-go-related gene (hERG) channel, is the most common cause of acquired long QT syndrome. Indeed, some, but not all, drugs that block hERG channels cause fatal cardiac arrhythmias. However, there is no clear method to distinguish between drugs that cause deadly arrhythmias and those that are clinically safe. Here we propose a mechanism that could explain why certain clinically used hERG blockers are less proarrhythmic than others. We demonstrate that several drugs that block hERG channels, but have favorable cardiac safety profiles, also evoke another effect; they facilitate the hERG current amplitude in response to low-voltage depolarization. To investigate how hERG facilitation impacts cardiac safety, we develop computational models of I Kr block with and without this facilitation. We constrain the models using data from voltage clamp recordings of hERG block and facilitation by nifekalant, a safe class III antiarrhythmic agent. Human ventricular action potential simulations demonstrate the ability of nifekalant to suppress ectopic excitations, with or without facilitation. Without facilitation, excessive I Kr block evokes early afterdepolarizations, which cause lethal arrhythmias. When facilitation is introduced, early afterdepolarizations are prevented at the same degree of block. Facilitation appears to prevent early afterdepolarizations by increasing I Kr during the repolarization phase of action potentials. We empirically test this prediction in isolated rabbit ventricular myocytes and find that action potential prolongation with nifekalant is less likely to induce early afterdepolarization than action potential prolongation with dofetilide, a hERG channel blocker that does not induce facilitation. Our data suggest that hERG channel blockers that induce facilitation increase the repolarization reserve of cardiac myocytes, rendering them less likely to trigger lethal ventricular arrhythmias.

Multiple Dynamical Mechanisms of Phase-2 Early Afterdepolarizations in a Human Ventricular Myocyte Model: Involvement of Spontaneous SR Ca 2+ Release.

Early afterdepolarization (EAD) is known to cause lethal ventricular arrhythmias in long QT syndrome (LQTS). In this study, dynamical mechanisms of EAD formation in human ventricular myocytes (HVMs) were investigated using the mathematical model developed by ten Tusscher and Panfilov (Am J Physiol Heart Circ Physiol 291, 2006). We explored how the rapid (IKr) and slow (IKs) components of delayed-rectifier K+ channel currents, L-type Ca2+ channel current (ICa L), Na+/Ca2+ exchanger current (INCX), and intracellular Ca2+ handling via the sarcoplasmic reticulum (SR) contribute to initiation, termination and modulation of phase-2 EADs during pacing in relation to bifurcation phenomena in non-paced model cells. Parameter-dependent dynamical behaviors of the non-paced model cell were determined by calculating stabilities of equilibrium points (EPs) and limit cycles, and bifurcation points to construct bifurcation diagrams. Action potentials (APs) and EADs during pacing were reproduced by numerical simulations for constructing phase diagrams of the paced model cell dynamics. Results are summarized as follows: (1) A modified version of the ten Tusscher-Panfilov model with accelerated ICaL inactivation could reproduce bradycardia-related EADs in LQTS type 2 and ß-adrenergic stimulation-induced EADs in LQTS type 1. (2) Two types of EADs with different initiation mechanisms, ICaL reactivation-dependent and spontaneous SR Ca2+ release-mediated EADs, were detected. (3) Termination of EADs (AP repolarization) during pacing depended on the slow activation of IKs. (4) Spontaneous SR Ca2+ releases occurred at higher Ca2+ uptake rates, attributable to the instability of steady-state intracellular Ca2+ concentrations. Dynamical mechanisms of EAD formation and termination in the paced model cell are closely related to stability changes (bifurcations) in dynamical behaviors of the non-paced model cell, but they are model-dependent. Nevertheless, the modified ten Tusscher-Panfilov model would be useful for systematically investigating possible dynamical mechanisms of EAD-related arrhythmias in LQTS.

Hysteretic Dynamics of Multi-Stable Early Afterdepolarisations with Repolarisation Reserve Attenuation: A Potential Dynamical Mechanism for Cardiac Arrhythmias.

Some cardiovascular and non-cardiovascular drugs frequently cause excessive prolongation of the cardiac action potential (AP) and lead to the development of early afterdepolarisations (EADs), which trigger lethal ventricular arrhythmias. Combining computer simulations in APs with numerical calculations based on dynamical system theory, we investigated stability changes of APs observed in a paced human ventricular myocyte model by decreasing and/or increasing the rapid (I Kr) and slow (I Ks) components of delayed rectifying K+ current. Upon reducing I Kr, the APs without EADs (no-EAD response) showed gradual prolongation of AP duration (APD), and were annihilated without AP configuration changes due to the occurrence of saddle-node bifurcations. This annihilation caused a transition to an AP with EADs as a new stable steady state. Furthermore, reducing repolarisation currents (repolarisation reserve attenuation) evoked multi-stable states consisting of APs with different APDs, and caused multiple hysteretic dynamics. Depending on initial ion circumstances within ventricular myocytes, these multi-stable AP states might increase the local/global heterogeneity of AP repolarisations in the ventricle. Thus, the EAD-induced arrhythmias with repolarisation reserve attenuation might be attributed to the APD variability caused by multi-stability in cardiac AP dynamics.

Integrative Approach with Electrophysiological and Theoretical Methods Reveals a New Role of S4 Positively Charged Residues in PKD2L1 Channel Voltage-Sensing.

Numerical model-based simulations provide important insights into ion channel gating when experimental limitations exist. Here, a novel strategy combining numerical simulations with patch clamp experiments was used to investigate the net positive charges in the putative transmembrane segment 4 (S4) of the atypical, positively-shifted voltage-dependence of polycystic kidney disease 2-like 1 (PKD2L1) channel. Charge-neutralising mutations (K452Q, K455Q and K461Q) in S4 reduced gating charges, positively shifted the Boltzmann-type activation curve [i.e., open probability (P open)-V curve] and altered the time-courses of activation/deactivation of PKD2L1, indicating that this region constitutes part of a voltage sensor. Numerical reconstruction of wild-type (WT) and mutant PKD2L1-mediated currents necessitated, besides their voltage-dependent gating parameters, a scaling factor that describes the voltage-dependence of maximal conductance, G max. Subsequent single-channel conductance (γ) measurements revealed that voltage-dependence of G max in WT can be explained by the inward-rectifying property of γ, which is greatly changed in PKD2L1 mutants. Homology modelling based on PKD2 and NaVAb structures suggest that such voltage dependence of P open and γ in PKD2L1 could both reflect the charged state of the S4 domain. The present conjunctive experimental and theoretical approaches provide a framework to explore the undetermined mechanism(s) regulating TRP channels that possess non-classical voltage-dependent properties.

HD Physiology Project-Japanese efforts to promote multilevel integrative systems biology and physiome research.

The HD Physiology Project is a Japanese research consortium that aimed to develop methods and a computational platform in which physiological and pathological information can be described in high-level definitions across multiple scales of time and size. During the 5 years of this project, an appropriate software platform for multilevel functional simulation was developed and a whole-heart model including pharmacokinetics for the assessment of the proarrhythmic risk of drugs was developed. In this article, we outline the description and scientific strategy of this project and present the achievements and influence on multilevel integrative systems biology and physiome research.

Dynamical mechanisms of phase-2 early afterdepolarizations in human ventricular myocytes: insights from bifurcation analyses of two mathematical models.

Early afterdepolarization (EAD) is known as a cause of ventricular arrhythmias in long QT syndromes. We theoretically investigated how the rapid (IKr) and slow (IKs) components of delayed-rectifier K+ channel currents, L-type Ca2+ channel current (ICaL), Na+/Ca2+ exchanger current (INCX), Na+-K+ pump current (INaK), intracellular Ca2+ (Cai) handling via sarcoplasmic reticulum (SR), and intracellular Na+ concentration (Nai) contribute to initiation, termination, and modulation of phase-2 EADs, using two human ventricular myocyte models. Bifurcation structures of dynamical behaviors in model cells were explored by calculating equilibrium points, limit cycles (LCs), and bifurcation points as functions of parameters. EADs were reproduced by numerical simulations. The results are summarized as follows: 1) decreasing IKs and/or IKr or increasing ICaL led to EAD generation, to which mid-myocardial cell models were especially susceptible; the parameter regions of EADs overlapped the regions of stable LCs. 2) Two types of EADs (termination mechanisms), IKs activation-dependent and ICaL inactivation-dependent EADs, were detected; IKs was not necessarily required for EAD formation. 3) Inhibiting INCX suppressed EADs via facilitating Ca2+-dependent ICaL inactivation. 4) Cai dynamics (SR Ca2+ handling) and Nai strongly affected bifurcations and EAD generation in model cells via modulating ICaL, INCX, and INaK Parameter regions of EADs, often overlapping those of stable LCs, shifted depending on Cai and Nai in stationary and dynamic states. 5) Bradycardia-related induction of EADs was mainly due to decreases in Nai at lower pacing rates. This study demonstrates that bifurcation analysis allows us to understand the dynamical mechanisms of EAD formation more profoundly. NEW & NOTEWORTHY: We investigated mechanisms of phase-2 early afterdepolarization (EAD) by bifurcation analyses of human ventricular myocyte (HVM) models. EAD formation in paced HVMs basically depended on bifurcation phenomena in non-paced HVMs, but was strongly affected by intracellular ion concentrations in stationary and dynamic states. EAD generation did not necessarily require IKs.

Ischemia-related subcellular redistribution of sodium channels enhances the proarrhythmic effect of class I antiarrhythmic drugs: a simulation study.

BACKGROUND: Cardiomyocytes located at the ischemic border zone of infarcted ventricle are accompanied by redistribution of gap junctions, which mediate electrical transmission between cardiomyocytes. This ischemic border zone provides an arrhythmogenic substrate. It was also shown that sodium (Na+) channels are redistributed within myocytes located in the ischemic border zone. However, the roles of the subcellular redistribution of Na+ channels in the arrhythmogenicity under ischemia remain unclear. METHODS: Computer simulations of excitation conduction were performed in a myofiber model incorporating both subcellular Na+ channel redistribution and the electric field mechanism, taking into account the intercellular cleft potentials. RESULTS: We found in the myofiber model that the subcellular redistribution of the Na+ channels under myocardial ischemia, decreasing in Na+ channel expression of the lateral cell membrane of each myocyte, decreased the tissue excitability, resulting in conduction slowing even without any ischemia-related electrophysiological change. The conventional model (i.e., without the electric field mechanism) did not reproduce the conduction slowing caused by the subcellular Na+ channel redistribution. Furthermore, Na+ channel blockade with the coexistence of a non-ischemic zone with an ischemic border zone expanded the vulnerable period for reentrant tachyarrhythmias compared to the model without the ischemic border zone. Na+ channel blockade tended to cause unidirectional conduction block at sites near the ischemic border zone. Thus, such a unidirectional conduction block induced by a premature stimulus at sites near the ischemic border zone is associated with the initiation of reentrant tachyarrhythmias. CONCLUSIONS: Proarrhythmia of Na+ channel blockade in patients with old myocardial infarction might be partly attributable to the ischemia-related subcellular Na+ channel redistribution.

Modeling light adaptation in circadian clock: prediction of the response that stabilizes entrainment.

Periods of biological clocks are close to but often different from the rotation period of the earth. Thus, the clocks of organisms must be adjusted to synchronize with day-night cycles. The primary signal that adjusts the clocks is light. In Neurospora, light transiently up-regulates the expression of specific clock genes. This molecular response to light is called light adaptation. Does light adaptation occur in other organisms? Using published experimental data, we first estimated the time course of the up-regulation rate of gene expression by light. Intriguingly, the estimated up-regulation rate was transient during light period in mice as well as Neurospora. Next, we constructed a computational model to consider how light adaptation had an effect on the entrainment of circadian oscillation to 24-h light-dark cycles. We found that cellular oscillations are more likely to be destabilized without light adaption especially when light intensity is very high. From the present results, we predict that the instability of circadian oscillations under 24-h light-dark cycles can be experimentally observed if light adaptation is altered. We conclude that the functional consequence of light adaptation is to increase the adjustability to 24-h light-dark cycles and then adapt to fluctuating environments in nature.

Roles of subcellular Na+ channel distributions in the mechanism of cardiac conduction.

The gap junction and voltage-gated Na(+) channel play an important role in the action potential propagation. The purpose of this study was to elucidate the roles of subcellular Na(+) channel distribution in action potential propagation. To achieve this, we constructed the myocardial strand model, which can calculate the current via intercellular cleft (electric-field mechanism) together with gap-junctional current (gap-junctional mechanism). We conducted simulations of action potential propagation in a myofiber model where cardiomyocytes were electrically coupled with gap junctions alone or with both the gap junctions and the electric field mechanism. Then we found that the action potential propagation was greatly affected by the subcellular distribution of Na(+) channels in the presence of the electric field mechanism. The presence of Na(+) channels in the lateral membrane was important to ensure the stability of propagation under conditions of reduced gap-junctional coupling. In the poorly coupled tissue with sufficient Na(+) channels in the lateral membrane, the slowing of action potential propagation resulted from the periodic and intermittent dysfunction of the electric field mechanism. The changes in the subcellular Na(+) channel distribution might be in part responsible for the homeostatic excitation propagation in the diseased heart.

Transmural dispersion of repolarization determines scroll wave behavior during ventricular tachyarrhythmias.

BACKGROUND: Ventricular tachyarrhythmia is the leading cause of sudden cardiac death, and scroll wave re-entry is known to underlie this condition. Class III antiarrhythmic drugs are commonly used worldwide to treat ventricular tachyarrhythmias; however, these drugs have a proarrhythmic adverse effect and can cause Torsade de Pointes or ventricular fibrillation. Transmural dispersion of repolarization (TDR) has been suggested to be a strong indicator of ventricular tachyarrhythmia induction. However, the role of TDR during sustained scroll wave re-entry is poorly understood. The purpose of the present study was to investigate how TDR affects scroll wave behavior and to provide a novel analysis of the mechanisms that sustain tachyarrhythmias, using computer simulations. METHODS AND RESULTS: Computer simulations were carried out to quantify the TDR and QT interval under a variety of I(Ks) and I(Kr) during transmural conduction. Simulated scroll wave re-entries were done under a variety of I(Ks) and I(Kr) in a ventricular wall slab model, and the scroll wave behavior and the filament dynamics (3-dimensional organizing center) were analyzed. A slight increase in TDR, but not in the QT interval, reflected antiarrhythmic properties resulting from the restraint of scroll wave breakup, whereas a marked increase in TDR was proarrhythmic, as a result of scroll wave breakup. CONCLUSIONS: The TDR determines the sustainment of ventricular tachyarrhythmias, through control of the scroll wave filament dynamics.

Synchronization of firing in cortical fast-spiking interneurons at gamma frequencies: a phase-resetting analysis.

Fast-spiking (FS) cells in the neocortex are interconnected both by inhibitory chemical synapses and by electrical synapses, or gap-junctions. Synchronized firing of FS neurons is important in the generation of gamma oscillations, at frequencies between 30 and 80 Hz. To understand how these synaptic interactions control synchronization, artificial synaptic conductances were injected in FS cells, and the synaptic phase-resetting function (SPRF), describing how the compound synaptic input perturbs the phase of gamma-frequency spiking as a function of the phase at which it is applied, was measured. GABAergic and gap junctional conductances made distinct contributions to the SPRF, which had a surprisingly simple piecewise linear form, with a sharp midcycle break between phase delay and advance. Analysis of the SPRF showed how the intrinsic biophysical properties of FS neurons and their interconnections allow entrainment of firing over a wide gamma frequency band, whose upper and lower frequency limits are controlled by electrical synapses and GABAergic inhibition respectively.

Nonlinear dynamics and bifurcations of a coupled oscillator model for calling behavior of Japanese tree frogs (Hyla japonica).

Synchronization has been observed in various systems, including living beings. In a previous study, we reported a new phenomenon with antisynchronization in calling behavior of two interacting Japanese tree frogs. In this paper, we theoretically analyse nonlinear dynamics in a system of three coupled oscillators, which models three interacting frogs, where the oscillators of each pair have the property of antisynchronization; in particular, we perform bifurcation analysis and Lyapunov function analysis.

Bidirectional modulation of neuronal responses by depolarizing GABAergic inputs.

The reversal potential of GABAA receptor channels is known to be less negative than the resting membrane potential under some cases. Recent electrophysiological experiments revealed that a GABAergic unitary conductance with such a depolarized reversal potential could not only prevent but also facilitate action potential generation depending on the timing of its application relative to the excitatory unitary conductance. Using a two-dimensional point neuron model, we simulate the experiments regarding the integration of unitary conductances, and execute bifurcation analysis. Then we extend our analysis to the case in which the neuron receives two kinds of periodic input trains-an excitatory one and a GABAergic one. We show that the periodic depolarizing GABAergic input train can modulate the output time-averaged firing rate bidirectionally, namely as an increase or a decrease, in a devil's-staircase-like manner depending on the phase difference with the excitatory input train. Bifurcation analysis reveals the existence of a wide variety of phase-locked solutions underlying such a graded response of the neuron. We examine how the input time-width and the value of the GABAA reversal potential affect the response. Moreover, considering a neuronal population, we show that depolarizing GABAergic inputs bidirectionally modulate the amplitude of the oscillatory population activity.

Bifurcations in a mathematical model for circadian oscillations of clock genes.

Circadian oscillations with a period of about 24h are observed in nearly all living organisms as conspicuous biological rhythms. In this paper, we investigate various kinds of bifurcation phenomena produced in a circadian oscillator model of Drosophila. In Drosophila, it is known that circadian oscillations in the levels of two proteins, PER and TIM, result from the negative feedback exerted by a PER-TIM complex on the expression of the per and tim genes that code for the two proteins. For studying circadian oscillations of proteins in Drosophila, a mathematical model has been proposed. The model cannot only account for regular circadian oscillations in environmental conditions such as constant darkness, but also give rise to more complex oscillatory phenomena including chaos and birhythmicity. By calculating bifurcations using Kawakami's method, we obtain detailed bifurcation diagrams related to stable and unstable invariant sets, and identify parameter regions in which the model generates complex oscillations as well as regular circadian oscillations. Moreover, we study bifurcations observed in the model incorporating the effect on a light-dark (LD) cycle and show that the waveform of the periodic variation in the light-induced parameter has a marked influence on the global bifurcation structure or the type of dynamic behavior resulting from the forcing term of the circadian oscillator by the LD cycles.

Possible effects of depolarizing GABAA conductance on the neuronal input-output relationship: a modeling study.

Recent in vitro experiments revealed that the GABAA reversal potential is about 10 mV higher than the resting potential in mature mammalian neocortical pyramidal cells; thus GABAergic inputs could have facilitatory, rather than inhibitory, effects on action potential generation under certain conditions. However, how the relationship between excitatory input conductances and the output firing rate is modulated by such depolarizing GABAergic inputs under in vivo circumstances has not yet been understood. We examine herewith the input-output relationship in a simple conductance-based model of cortical neurons with the depolarized GABAA reversal potential, and show that a tonic depolarizing GABAergic conductance up to a certain amount does not change the relationship between a tonic glutamatergic driving conductance and the output firing rate, whereas a higher GABAergic conductance prevents spike generation. When the tonic glutamatergic and GABAergic conductances are replaced by in vivo-like highly fluctuating inputs, on the other hand, the effect of depolarizing GABAergic inputs on the input-output relationship critically depends on the degree of coincidence between glutamatergic input events and GABAergic ones. Although a wide range of depolarizing GABAergic inputs hardly changes the firing rate of a neuron driven by noncoincident glutamatergic inputs, a certain range of these inputs considerably decreases the firing rate if a large number of driving glutamatergic inputs are coincident with them. These results raise the possibility that the depolarized GABAA reversal potential is not a paradoxical mystery, but is instead a sophisticated device for discriminative firing rate modulation.

Bifurcations of synchronized responses in synaptically coupled Bonhöffer-van der Pol neurons.

The Bonhöffer-van der Pol (BvdP) equation is considered as an important model for studying dynamics in a single neuron. In this paper, we investigate bifurcations of periodic solutions in model equations of four and five BvdP neurons coupled through the characteristics of synaptic transmissions with a time delay. The model can be considered as a dynamical system whose solution includes jumps depending on a condition related to the behavior of the trajectory. Although the solution is discontinuous, we can define the Poincaré map as a synthesis of successive submaps, and give its derivatives for obtaining periodic points and their bifurcations. Using our proposed numerical method, we clarify mechanisms of bifurcations among synchronized oscillations with phase-locking patterns by analyzing periodic solutions observed in the coupling system and its subsystems. Moreover, we show that a global behavior of chaotic itinerancy or a phenomenon of chaotic transitions among several quasiattracting states can be observed in higher-dimensional systems of the synaptically four and five coupled neurons.