Modeling the dynamic features of the electrogenic Na,K pump of cardiac cells

A Thermodynamic Model of Monovalent Cation Homeostasis in the Yeast Saccharomyces cerevisiae

Cationic and heavy metal toxicity is involved in a substantial number of diseases in mammals and crop plants. Therefore, the understanding of tightly regulated transporter activities, as well as conceiving the interplay of regulatory mechanisms, is of substantial interest. A generalized thermodynamic description is developed for the complex interplay of the plasma membrane ion transporters, membrane potential and the consumption of energy for maintaining and restoring specific intracellular cation concentrations. This concept is applied to the homeostasis of cation concentrations in the yeast cells of S. cerevisiae. The thermodynamic approach allows to model passive ion fluxes driven by the electrochemical potential differences, but also primary or secondary active transport processes driven by the inter- play of different ions (symport, antiport) or by ATP consumption (ATPases). The model-confronted with experimental data-reproduces the experimentally observed potassium and proton fluxes induced by the external stimuli KCl and glucose. The estimated phenomenological constants combine kinetic parameters and transport coefficients. These are in good agreement with the biological understanding of the transporters thus providing a better understanding of the control exerted by the coupled fluxes. The model predicts the flux of additional ion species, like e.g. chloride, as a potential candidate for counterbalancing positive charges. Furthermore, the effect of a second KCl stimulus is simulated, predicting a reduced cellular response for cells that were first exposed to a high KCl stimulus compared to cells pretreated with a mild KCl stimulus. By describing the generalized forces that are responsible for a given flow, the model provides information and suggestions for new experiments. Furthermore, it can be extended to other systems such as e.g. Candida albicans, or selected plant cells.

Coupled left-shift of Nav channels: modeling the Na⁺-loading and dysfunctional excitability of damaged axons

Injury to neural tissue renders voltage-gated Na⁺ (Nav) channels leaky. Even mild axonal trauma initiates Na⁺-loading, leading to secondary Ca²⁺-loading and white matter degeneration. The nodal isoform is Nav1.6 and for Nav1.6-expressing HEK-cells, traumatic whole cell stretch causes an immediate tetrodotoxin-sensitive Na⁺-leak. In stretch-damaged oocyte patches, Nav1.6 current undergoes damage-intensity dependent hyperpolarizing- (left-) shifts, but whether left-shift underlies injured-axon Nav-leak is uncertain. Nav1.6 inactivation (availability) is kinetically limited by (coupled to) Nav activation, yielding coupled left-shift (CLS) of the two processes: CLS should move the steady-state Nav1.6 "window conductance" closer to typical firing thresholds. Here we simulated excitability and ion homeostasis in free-running nodes of Ranvier to assess if hallmark injured-axon behaviors--Na⁺-loading, ectopic excitation, propagation block--would occur with Nav-CLS. Intact/traumatized axolemma ratios were varied, and for some simulations Na/K pumps were included, with varied in/outside volumes. We simulated saltatory propagation with one mid-axon node variously traumatized. While dissipating the [Na⁺] gradient and hyperactivating the Na/K pump, Nav-CLS generated neuropathic pain-like ectopic bursts. Depending on CLS magnitude, fraction of Nav channels affected, and pump intensity, tonic or burst firing or nodal inexcitability occurred, with [Na⁺] and [K⁺] fluctuating. Severe CLS-induced inexcitability did not preclude Na⁺-loading; in fact, the steady-state Na⁺-leaks elicited large pump currents. At a mid-axon node, mild CLS perturbed normal anterograde propagation, and severe CLS blocked saltatory propagation. These results suggest that in damaged excitable cells, Nav-CLS could initiate cellular deterioration with attendant hyper- or hypo-excitability. Healthy-cell versions of Nav-CLS, however, could contribute to physiological rhythmic firing.

A mathematical model of ion movements in grey matter during a stroke

The development of cytotoxic oedema during a stroke consists in cell swelling and shrinking of the extracellular space. This phenomenon is triggered by ion movements through voltage-gated channels, exchangers and pumps. During ischaemia, sodium, calcium and chloride enter the neurons whereas potassium and glutamate are expelled out of the cells. A mathematical model is proposed to represent the long-term dynamics of membrane potentials, cell volumes and ionic concentrations in intracellular and extracellular spaces during a stroke and to study the influence of each ionic current on cell swelling. The model relies on electrophysiological mechanisms and takes into account the behaviour of two types of cells: neurons and also astrocytes known to play a key role in the excitotoxic process in grey matter. The results obtained when a severe or a moderate ischaemia is simulated are consistent with those observed in the in vitro and in vivo experiments. As this model appears to be robust, it is used to perform illustrative simulations aimed at studying the effect of some channel blockers on cell swelling. This approach may help to explore new therapeutic strategies in order to reduce stroke damage.

Osmotic forces and gap junctions in spreading depression: a computational model

In a computational model of spreading depression (SD), ionic movement through a neuronal syncytium of cells connected by gap junctions is described electrodiffusively. Simulations predict that SD will not occur unless cells are allowed to expand in response to osmotic pressure gradients and K+ is allowed to move through gap junctions. SD waves of [K+]out approximately 25 to approximately 60 mM moving at approximately 2 to approximately 18 mm/min are predicted over the range of parametric values reported in gray matter, with extracellular space decreasing up to approximately 50%. Predicted waveform shape is qualitatively similar to laboratory reports. The delayed-rectifier, NMDA, BK, and Na+ currents are predicted to facilitate SD, while SK and A-type K+ currents and glial activity impede SD. These predictions are consonant with recent findings that gap junction poisons block SD and support the theories that cytosolic diffusion via gap junctions and osmotic forces are important mechanisms underlying SD.

Hysteresis and bistability in the direct transition from 1:1 to 2:1 rhythm in periodically driven single ventricular cells

The transmembrane potential of a single quiescent cell isolated from rabbit ventricular muscle was recorded using a suction electrode in whole-cell recording mode. The cell was then driven with a periodic train of current pulses injected into the cell through the same recording electrode. When the interpulse interval or basic cycle length (BCL) was sufficiently long, 1:1 rhythm resulted, with each stimulus pulse producing an action potential. Gradual decrease in BCL invariably resulted in loss of 1:1 synchronization at some point. When the pulse amplitude was set to a fixed low level and BCL gradually decreased, N+1:N rhythms (N>/=2) reminiscent of clinically observed Wenckebach rhythms were seen. Further decrease in BCL then yielded a 2:1 rhythm. In contrast, when the pulse amplitude was set to a fixed high level, a period-doubled 2:2 rhythm resembling alternans rhythm was seen before a 2:1 rhythm occurred. With the pulse amplitude set to an intermediate level (i.e., to a level between those at which Wenckebach and alternans rhythms were seen), there was a direct transition from 1:1 to 2:1 rhythm as the BCL was decreased: Wenckebach and alternans rhythms were not seen. When at that point the BCL was increased, the transition back to 1:1 rhythm occurred at a longer BCL than that at which the {1:1-->2:1} transition had initially occurred, demonstrating hysteresis. With the BCL set to a value within the hysteresis range, injection of a single well-timed extrastimulus converted 1:1 rhythm into 2:1 rhythm or vice versa, providing incontrovertible evidence of bistability (the coexistence of two different periodic rhythms at a fixed set of stimulation parameters). Hysteresis between 1:1 and 2:1 rhythms was also seen when the stimulus amplitude, rather than the BCL, was changed. Simulations using numerical integration of an ionic model of a single ventricular cell formulated as a nonlinear system of differential equations provided results that were very similar to those found in the experiments. The steady-state action potential duration restitution curve, which is a plot of the duration of the action potential during 1:1 rhythm as a function of the recovery time or diastolic interval immediately preceding that action potential, was determined. Iteration of a finite-difference equation derived using the restitution curve predicted the direct {1:1<-->2:1} transition, as well as bistability, in both the experimental and modeling work. However, prediction of the action potential duration during 2:1 rhythm was not as accurate in the experiments as in the model. Finally, we point out a few implications of our findings for cardiac arrhythmias (e.g., Mobitz type II block, ischemic alternans). (c) 1999 American Institute of Physics.

Modeling cell volume regulation in nonexcitable cells: the roles of the Na+ pump and of cotransport systems

The purpose of this study is to contribute to understanding the role of Na+-K+-ATPase and of ionic cotransporters in the regulation of cell volume, by employing a model that describes the rates of change of the intracellular concentrations of Na+, K+, and Cl-, of the cell volume, and of the membrane potential. In most previous models of dynamic cellular phenomena, Na+-K+-ATPase is incorporated via phenomenological formulations; the enzyme is incorporated here via an explicit kinetic scheme. Another feature of the present model is the capability to perform short-term cell volume regulation mediated by cotransporters of KCl and NaCl. The model is employed to perform numerical simulations for a "typical" nonpolarized animal cell. Basically, the results are consistent with the view that the Na+ pump mainly plays a long-term role in the maintenance of the electrochemical gradients of Na+ and K+ and that short-term cell volume regulation is achieved via passive transport, exemplified in this case by the cotransport of KCl and NaCl.

Pump and exchanger mechanisms in a model of smooth muscle

A novel approach to modelling pump and exchanger mechanisms is presented. In this approach, new thermodynamic expressions for the calcium pump, sodium-calcium exchanger and sodium-potassium pump are developed using statistical rate theory (SRT). This theory is well-defined and is not derived empirically. This is in contrast to previous thermodynamic pump expressions which used a simple linear relationship or relied on empirical data for their functional form. The functional form of these new expressions does not require assumptions of steady state or particular forms of voltage dependencies in specific steps. Also, the explicit reaction scheme is not required. Instead, assumptions of a rate-limiting step in the scheme and a near-equilibrium ratio of intermediate substrates are required. These expressions are incorporated into an overall model of gastric smooth muscle. This model presents a novel approach whereby thermodynamic representations for calcium pumps, sodium-calcium exchangers and sodium-potassium pumps have been included together in a model of ionic transport mechanisms for smooth muscle. Variations in basal metabolic concentrations are used to explain the observed amplitude variation in the transmembrane voltage of gastric smooth muscle. The interaction of the various mechanisms are used to illustrate the large depolarization obtained in smooth muscle with ouabain as well as the forward and reverse modes of the sodium-calcium exchanger.