Different roles for aspartates and glutamates for cation permeation in bacterial sodium channels

A key driving force for ion channel selectivity is represented by the negative charge of the Selectivity Filter carried by aspartate (D) and glutamate (E) residues. However, the structural effects and specific properties of D and E residues have not been extensively studied. In order to investigate this issue we studied the mutants of NaChBac channel with all possible combinations of D and E in the charged rings in position 191 and 192. Electrophysiological measurements showed significant Ca²⁺ currents only when position 191 was occupied by E. Equilibrium Molecular Dynamics simulations revealed the existence of two binding sites, corresponding to the charged rings and another one, more internal, at the level of L190. The simulations showed that the ion in the innermost site can interact with the residue in position 191 only when this is glutamate. Based on the MD simulations, we suggest that a D in position 191 leads to a high affinity Ca²⁺ block site resulting from a significant drop in the free energy of binding for an ion moving between the binding sites; in contrast, the free energy change is more gradual when an E residue occupies position 191, resulting in Ca²⁺ permeability. This scenario is consistent with the model of ion channel selectivity through stepwise changes in binding affinity proposed by Dang and McCleskey. Our study also highlights the importance of the structure of the selectivity filter which should contribute to the development of more detailed physical models for ion channel selectivity.     

Simulation of the effect of an external GHz electric field on the potential energy profile of Ca2+ ions in the selectivity filter of the CaVAb channel

Ca_V channels are transmembrane proteins that mediate and regulate ion fluxes across cell membranes, and they are activated in response to action potentials to allow Ca²⁺ influx. Since ion channels are composed of charge or polar groups, an external alternating electric field may affect the ion‐selective membrane transport and the performance of the channel. In this article, we have investigated the effect of an external GHz electric field on the dynamics of calcium ions in the selectivity filter of the Ca_VAb channel. Molecular dynamics (MD) simulations and the potential of mean force (PMF) calculations were carried out, via the umbrella sampling method, to determine the free energy profile of Ca²⁺ ions in the Ca_VAb channels in presence and absence of an external field. Exposing Ca_VAb channel to 1, 2, 3, 4, and 5 GHz electric fields increases the depth of the potential energy well and this may result in an increase in the affinity and strength of Ca²⁺ ions to binding sites in the selectivity filter the channel. This increase of strength of Ca²⁺ ions binding in the selectivity filter may interrupt the mechanism of Ca²⁺ ion conduction, and leads to a reduction of Ca²⁺ ion permeation through the Ca_VAb channel.     

Structural implications of weak Ca2+ block in Drosophila cyclic nucleotide–gated channels

Calcium permeability and the concomitant calcium block of monovalent ion current (“Ca²⁺ block”) are properties of cyclic nucleotide–gated (CNG) channel fundamental to visual and olfactory signal transduction. Although most CNG channels bear a conserved glutamate residue crucial for Ca²⁺ block, the degree of block displayed by different CNG channels varies greatly. For instance, the Drosophila melanogaster CNG channel shows only weak Ca²⁺ block despite the presence of this glutamate. We previously constructed a series of chimeric channels in which we replaced the selectivity filter of the bacterial nonselective cation channel NaK with a set of CNG channel filter sequences and determined that the resulting NaK2CNG chimeras displayed the ion selectivity and Ca²⁺ block properties of the parent CNG channels. Here, we used the same strategy to determine the structural basis of the weak Ca²⁺ block observed in the Drosophila CNG channel. The selectivity filter of the Drosophila CNG channel is similar to that of most other CNG channels except that it has a threonine at residue 318 instead of a proline. We constructed a NaK chimera, which we called NaK2CNG-Dm, which contained the Drosophila selectivity filter sequence. The high resolution structure of NaK2CNG-Dm revealed a filter structure different from those of NaK and all other previously investigated NaK2CNG chimeric channels. Consistent with this structural difference, functional studies of the NaK2CNG-Dm chimeric channel demonstrated a loss of Ca²⁺ block compared with other NaK2CNG chimeras. Moreover, mutating the corresponding threonine (T318) to proline in Drosophila CNG channels increased Ca²⁺ block by 16 times. These results imply that a simple replacement of a threonine for a proline in Drosophila CNG channels has likely given rise to a distinct selectivity filter conformation that results in weak Ca²⁺ block.     

Ryanodine receptor allosteric coupling and the dynamics of calcium sparks

Puffs and sparks are localized intracellular Ca²⁺ elevations that arise from the cooperative activity of Ca²⁺-regulated inositol 1,4,5-trisphosphate receptors and ryanodine receptors clustered at Ca²⁺ release sites on the surface of the endoplasmic reticulum or the sarcoplasmic reticulum. While the synchronous gating of Ca²⁺-regulated Ca²⁺ channels can be mediated entirely though the buffered diffusion of intracellular Ca²⁺, interprotein allosteric interactions also contribute to the dynamics of ryanodine receptor (RyR) gating and Ca²⁺ sparks. In this article, Markov chain models of Ca²⁺ release sites are used to investigate how the statistics of Ca²⁺ spark generation and termination are related to the coupling of RyRs via local [Ca²⁺] changes and allosteric interactions. Allosteric interactions are included in a manner that promotes the synchronous gating of channels by stabilizing neighboring closed-closed and/or open-open channel pairs. When the strength of Ca²⁺-mediated channel coupling is systematically varied (e.g., by changing the Ca²⁺ buffer concentration), simulations that include synchronizing allosteric interactions often exhibit more robust Ca²⁺ sparks; however, for some Ca²⁺ coupling strengths the sparks are less robust. We find no evidence that the distribution of spark durations can be used to distinguish between allosteric interactions that stabilize closed channel pairs, open channel pairs, or both in a balanced fashion. On the other hand, the changes in spark duration, interspark interval, and frequency observed when allosteric interactions that stabilize closed channel pairs are gradually removed from simulations are qualitatively different than the changes observed when open or both closed and open channel pairs are stabilized. Thus, our simulations clarify how changes in spark statistics due to pharmacological washout of the accessory proteins mediating allosteric coupling may indicate the type of synchronizing allosteric interactions exhibited by physically coupled RyRs. We also investigate the validity of a mean-field reduction applicable to the dynamics of a ryanodine receptor cluster coupled via local [Ca²⁺] and allosteric interactions. In addition to facilitating parameter studies of the effect of allosteric coupling on spark statistics, the derivation of the mean-field model establishes the correct functional form for cooperativity factors representing the coupled gating of RyRs. This mean-field formulation is well suited for use in computationally efficient whole cell simulations of excitation-contraction coupling.     

Ionic selectivity in L-type calcium channels by electrostatics and hard-core repulsion

A physical model of selective “ion binding” in the L-type calcium channel is constructed, and consequences of the model are compared with experimental data. This reduced model treats only ions and the carboxylate oxygens of the EEEE locus explicitly and restricts interactions to hard-core repulsion and ion–ion and ion–dielectric electrostatic forces. The structural atoms provide a flexible environment for passing cations, thus resulting in a self-organized induced-fit model of the selectivity filter. Experimental conditions involving binary mixtures of alkali and/or alkaline earth metal ions are computed using equilibrium Monte Carlo simulations in the grand canonical ensemble. The model pore rejects alkali metal ions in the presence of biological concentrations of Ca²⁺ and predicts the blockade of alkali metal ion currents by micromolar Ca²⁺. Conductance patterns observed in varied mixtures containing Na⁺ and Li⁺, or Ba²⁺ and Ca²⁺, are predicted. Ca²⁺ is substantially more potent in blocking Na⁺ current than Ba²⁺. In apparent contrast to experiments using buffered Ca²⁺ solutions, the predicted potency of Ca²⁺ in blocking alkali metal ion currents depends on the species and concentration of the alkali metal ion, as is expected if these ions compete with Ca²⁺ for the pore. These experiments depend on the problematic estimation of Ca²⁺ activity in solutions buffered for Ca²⁺ and pH in a varying background of bulk salt. Simulations of Ca²⁺ distribution with the model pore bathed in solutions containing a varied amount of Li⁺ reveal a “barrier and well” pattern. The entry/exit barrier for Ca²⁺ is strongly modulated by the Li⁺ concentration of the bath, suggesting a physical explanation for observed kinetic phenomena. Our simulations show that the selectivity of L-type calcium channels can arise from an interplay of electrostatic and hard-core repulsion forces among ions and a few crucial channel atoms. The reduced system selects for the cation that delivers the largest charge in the smallest ion volume.     

Energetics of divalent selectivity in a calcium channel: the ryanodine receptor case study

A model of the ryanodine receptor (RyR) calcium channel is used to study the energetics of binding selectivity of Ca²⁺ versus monovalent cations. RyR is a calcium-selective channel with a DDDD locus in the selectivity filter, similar to the EEEE locus of the L-type calcium channel. While the affinity of RyR for Ca²⁺ is in the millimolar range (as opposed to the micromolar range of the L-type channel), the ease of single-channel measurements compared to L-type and its similar selectivity filter make RyR an excellent candidate for studying calcium selectivity. A Poisson-Nernst-Planck/density functional theory model of RyR is used to calculate the energetics of selectivity. Ca²⁺ versus monovalent selectivity is driven by the charge/space competition mechanism in which selectivity arises from a balance of electrostatics and the excluded volume of ions in the crowded selectivity filter. While electrostatic terms dominate the selectivity, the much smaller excluded-volume term also plays a substantial role. In the D4899N and D4938N mutations of RyR that are analyzed, substantial changes in specific components of the chemical potential profiles are found far from the mutation site. These changes result in the significant reduction of Ca²⁺ selectivity found in both theory and experiments.     

Structural effects of divalent calcium cations on the α7 nicotinic acetylcholine receptor: A molecular dynamics simulation study

The α7 subtype of neuronal nicotinic acetylcholine receptor (nAChR) is a ligand-gated ion channel protein that is vital to various neurological functions, including modulation of neurotransmitter release. A relatively high concentration of extracellular Ca²⁺ in the neuronal environment is likely to exert substantial structural and functional influence on nAChRs, which may affect their interactions with agonists and antagonists. In this work, we employed atomistic molecular dynamics (MD) simulations to examine the effects of elevated Ca²⁺ on the structure and dynamics of α7 nAChR embedded in a model phospholipid bilayer. Our results suggest that the presence of Ca²⁺ in the α7 nAChR environment results in closure of loop C-in the extracellular ligand-binding domain, a motion normally associated with agonist binding and receptor activation. Elevated Ca²⁺ also alters the conformation of key regions of the receptor, including the inter-helical loops, pore-lining helices and the “gate” residues, and causes partial channel opening in the absence of an agonist, leading to an attendant reduction in the free energy of Ca²⁺ permeation through the pore as elucidated by umbrella sampling simulations. Overall, the structural and permeability changes in α7 nAChR suggest that elevated Ca²⁺ induces a partially activated receptor state that is distinct from both the resting and the agonist-activated states. These results are consistent with the notion that divalent ions can serve as a potentiator of nAChRs, resulting in a higher rate of receptor activation (and subsequent desensitization) in the presence of agonists, with possible implications for diseases involving calcium dysregulation.     

Probing ion binding in the selectivity filter of the Cav1.1 channel with molecular dynamics

Ca_v1.1 is the voltage-gated calcium channel essential for the contraction of skeletal muscles upon membrane potential changes. Structural determination of the Ca_v1.1 channel opens the avenue toward understanding of the structure-function relationship of voltage-gated calcium channels. Here, we show that there exist two Ca²⁺-binding sites, termed S1 and S2, within the selectivity filter of Ca_v1.1 through extensive molecular dynamics simulations on various initial ion arrangement configurations. The formation of both binding sites is associated with the four Glu residues (Glu292/614/1014/1323) that constitute the so-called EEEE locus. At the S1 site near the extracellular side, the Ca²⁺ ion is coordinated with the negatively charged carboxylic groups of these Glu residues and of the Asp615 residue either in a direct way or via an intermediate water molecule. At the S2 site, Ca²⁺ binding shows two distinct states: an upper state involving two out of the four Glu residues in the EEEE locus and a lower state involving only one Glu residue. In addition, there exist two recruitment sites for Ca²⁺ above the entrance of the filter. These findings promote the understanding of mechanism for ion permeation and selectivity in calcium channels.     

Close agreement between deterministic versus stochastic modeling of first-passage time to vesicle fusion

Ca²⁺-dependent cell processes, such as neurotransmitter or endocrine vesicle fusion, are inherently stochastic due to large fluctuations in Ca²⁺ channel gating, Ca²⁺ diffusion, and Ca²⁺ binding to buffers and target sensors. However, previous studies revealed closer-than-expected agreement between deterministic and stochastic simulations of Ca²⁺ diffusion, buffering, and sensing if Ca²⁺ channel gating is not Ca²⁺ dependent. To understand this result more fully, we present a comparative study complementing previous work, focusing on Ca²⁺ dynamics downstream of Ca²⁺ channel gating. Specifically, we compare deterministic (mean-field/mass-action) and stochastic simulations of vesicle exocytosis latency, quantified by the probability density of the first-passage time (FPT) to the Ca²⁺-bound state of a vesicle fusion sensor, following a brief Ca²⁺ current pulse. We show that under physiological constraints, the discrepancy between FPT densities obtained using the two approaches remains small even if as few as ～50 Ca²⁺ ions enter per single channel-vesicle release unit. Using a reduced two-compartment model for ease of analysis, we illustrate how this close agreement arises from the smallness of correlations between fluctuations of the reactant molecule numbers, despite the large magnitude of fluctuation amplitudes. This holds if all relevant reactions are heteroreaction between molecules of different species, as is the case for bimolecular Ca²⁺ binding to buffers and downstream sensor targets. In this case, diffusion and buffering effectively decorrelate the state of the Ca²⁺ sensor from local Ca²⁺ fluctuations. Thus, fluctuations in the Ca²⁺ sensor’s state underlying the FPT distribution are only weakly affected by the fluctuations in the local Ca²⁺ concentration around its average, deterministically computable value.     

Mechanisms of Beat-to-Beat Regulation of Cardiac Pacemaker Cell Function by Ca Cycling Dynamics

Whether intracellular Ca²⁺ cycling dynamics regulate cardiac pacemaker cell function on a beat-to-beat basis remains unknown. Here we show that under physiological conditions, application of low concentrations of caffeine (2–4 mM) to isolated single rabbit sinoatrial node cells acutely reduces their spontaneous action potential cycle length (CL) and increases Ca²⁺ transient amplitude for several cycles. Numerical simulations, using a modified Maltsev-Lakatta coupled-clock model, faithfully reproduced these effects, and also the effects of CL prolongation and dysrhythmic spontaneous beating (produced by cytosolic Ca²⁺ buffering) and an acute CL reduction (produced by flash-induced Ca²⁺ release from a caged Ca²⁺ buffer), which we had reported previously. Three contemporary numerical models (including the original Maltsev-Lakatta model) failed to reproduce the experimental results. In our proposed new model, Ca²⁺ releases acutely change the CL via activation of the Na⁺/Ca²⁺ exchanger current. Time-dependent CL reductions after flash-induced Ca²⁺ releases (the memory effect) are linked to changes in Ca²⁺ available for pumping into sarcoplasmic reticulum which, in turn, changes the sarcoplasmic reticulum Ca²⁺ load, diastolic Ca²⁺ releases, and Na⁺/Ca²⁺ exchanger current. These results support the idea that Ca²⁺ regulates CL in cardiac pacemaker cells on a beat-to-beat basis, and suggest a more realistic numerical mechanism of this regulation.     

Cation Permeability and Selectivity of a Root Plasma Membrane Calcium Channel

Calcium channels in the plasma membrane of root cells fulfill both nutritional and signaling roles. The permeability of these channels to different cations determines the magnitude of their cation conductances, their effects on cell membrane potential and their contribution to cation toxicities. The selectivity of the rca channel, a Ca²⁺-permeable channel from the plasma membrane of wheat (Triticum aestivum L.) roots, was studied following its incorporation into planar lipid bilayers. The permeation of K⁺, Na⁺, Ca²⁺ and Mg²⁺ through the pore of the rca channel was modeled. It was assumed that cations permeated in single file through a pore with three energy barriers and two ion-binding sites. Differences in permeation between divalent and monovalent cations were attributed largely to the affinity of the ion binding sites. The model suggested that significant negative surface charge was present in the vestibules to the pore and that the pore could accommodate two cations simultaneously, which repelled each other strongly. The pore structure of the rca channel appeared to differ from that of L-type calcium channels from animal cell membranes since its ion binding sites had a lower affinity for divalent cations. The model adequately accounted for the diverse permeation phenomena observed for the rca channel. It described the apparent submillimolar K _m for the relationship between unitary conductance and Ca²⁺ activity, the differences in selectivity sequences obtained from measurements of conductance and permeability ratios, the changes in relative cation permeabilities with solution ionic composition, and the complex effects of Ca²⁺ on K⁺ and Na⁺ currents through the channel. Having established the adequacy of the model, it was used to predict the unitary currents that would be observed under the ionic conditions employed in patch-clamp experiments and to demonstrate the high selectivity of the rca channel for Ca²⁺ influx under physiological conditions. Received: 23 August 1999/Revised: 12 November 1999     

Use-dependent block of the voltage-gated Na+ channel by tetrodotoxin and saxitoxin: Effect of pore mutations that change ionic selectivity

Voltage-gated Na⁺ channels (NaV channels) are specifically blocked by guanidinium toxins such as tetrodotoxin (TTX) and saxitoxin (STX) with nanomolar to micromolar affinity depending on key amino acid substitutions in the outer vestibule of the channel that vary with NaV gene isoforms. All NaV channels that have been studied exhibit a use-dependent enhancement of TTX/STX affinity when the channel is stimulated with brief repetitive voltage depolarizations from a hyperpolarized starting voltage. Two models have been proposed to explain the mechanism of TTX/STX use dependence: a conformational mechanism and a trapped ion mechanism. In this study, we used selectivity filter mutations (K1237R, K1237A, and K1237H) of the rat muscle NaV1.4 channel that are known to alter ionic selectivity and Ca²⁺ permeability to test the trapped ion mechanism, which attributes use-dependent enhancement of toxin affinity to electrostatic repulsion between the bound toxin and Ca²⁺ or Na⁺ ions trapped inside the channel vestibule in the closed state. Our results indicate that TTX/STX use dependence is not relieved by mutations that enhance Ca²⁺ permeability, suggesting that ion–toxin repulsion is not the primary factor that determines use dependence. Evidence now favors the idea that TTX/STX use dependence arises from conformational coupling of the voltage sensor domain or domains with residues in the toxin-binding site that are also involved in slow inactivation.     

Divergence of Ca2+ selectivity and equilibrium Ca2+ blockade in a Ca2+ release-activated Ca2+ channel

Prevailing models postulate that high Ca²⁺ selectivity of Ca²⁺ release-activated Ca²⁺ (CRAC) channels arises from tight Ca²⁺ binding to a high affinity site within the pore, thereby blocking monovalent ion flux. Here, we examined the contribution of high affinity Ca²⁺ binding for Ca²⁺ selectivity in recombinant Orai3 channels, which function as highly Ca²⁺-selective channels when gated by the endoplasmic reticulum Ca²⁺ sensor STIM1 or as poorly Ca²⁺-selective channels when activated by the small molecule 2-aminoethoxydiphenyl borate (2-APB). Extracellular Ca²⁺ blocked Na⁺ currents in both gating modes with a similar inhibition constant (K_i; ∼25 µM). Thus, equilibrium binding as set by the K_i of Ca²⁺ blockade cannot explain the differing Ca²⁺ selectivity of the two gating modes. Unlike STIM1-gated channels, Ca²⁺ blockade in 2-APB–gated channels depended on the extracellular Na⁺ concentration and exhibited an anomalously steep voltage dependence, consistent with enhanced Na⁺ pore occupancy. Moreover, the second-order rate constants of Ca²⁺ blockade were eightfold faster in 2-APB–gated channels than in STIM1-gated channels. A four-barrier, three–binding site Eyring model indicated that lowering the entry and exit energy barriers for Ca²⁺ and Na⁺ to simulate the faster rate constants of 2-APB–gated channels qualitatively reproduces their low Ca²⁺ selectivity, suggesting that ion entry and exit rates strongly affect Ca²⁺ selectivity. Noise analysis indicated that the unitary Na⁺ conductance of 2-APB–gated channels is fourfold larger than that of STIM1-gated channels, but both modes of gating show a high open probability (P_o; ∼0.7). The increase in current noise during channel activation was consistent with stepwise recruitment of closed channels to a high P_o state in both cases, suggesting that the underlying gating mechanisms are operationally similar in the two gating modes. These results suggest that both high affinity Ca²⁺ binding and kinetic factors contribute to high Ca²⁺ selectivity in CRAC channels.     

A Structural Model of the Pore-Forming Region of the Skeletal Muscle Ryanodine Receptor (RyR1)

Ryanodine receptors (RyRs) are ion channels that regulate muscle contraction by releasing calcium ions from intracellular stores into the cytoplasm. Mutations in skeletal muscle RyR (RyR1) give rise to congenital diseases such as central core disease. The absence of high-resolution structures of RyR1 has limited our understanding of channel function and disease mechanisms at the molecular level. Here, we report a structural model of the pore-forming region of RyR1. Molecular dynamics simulations show high ion binding to putative pore residues D4899, E4900, D4938, and D4945, which are experimentally known to be critical for channel conductance and selectivity. We also observe preferential localization of Ca²⁺ over K⁺ in the selectivity filter of RyR1. Simulations of RyR1-D4899Q mutant show a loss of preference to Ca²⁺ in the selectivity filter as seen experimentally. Electrophysiological experiments on a central core disease mutant, RyR1-G4898R, show constitutively open channels that conduct K⁺ but not Ca²⁺. Our simulations with G4898R likewise show a decrease in the preference of Ca²⁺ over K⁺ in the selectivity filter. Together, the computational and experimental results shed light on ion conductance and selectivity of RyR1 at an atomistic level.     

Cations Residing at the Selectivity Filter of the Voltage-Gated Ca2+-Channel Modify Fusion-Pore Kinetics

Strontium (Sr²⁺), Barium (Ba²⁺), and Lanthanum (La³⁺) can substitute for Ca²⁺ in driving synaptic transmission during membrane depolarization. Ion recognition at the polyglutamate motif (EEEE), comprising the channel selectivity-filter, during voltage-driven transitions, controls the kinetics of the voltage-gated calcium channel (VGCC) and its interactions with the synaptic proteins. We tested the effect of different charge carriers on evoked-release, as a means of exploring the involvement of VGCC in the fusion pore configuration. Employing amperometry recordings in single bovine chromaffin cells we show that the size of the fusion pore, designated by the 'foot'-amplitude, was increased when Ba²⁺ substituted for Ca²⁺ and decreased, with La³⁺. The fusion pore stability, indicated by 'foot'-width, was decreased in La³⁺. Also, the mean open time of the fusion pore (tfp) was significantly lower in Sr²⁺ and La³⁺ compared to Ba²⁺ and Ca²⁺. These cations when occupying the selectivity filter reduced the spike frequency in the order of Ca²⁺ > Sr²⁺ > Ba²⁺ > La³⁺, which is parallel to the reduction in total catecholamine release. The correlation between ion binding at the selectivity filter and fusion pore properties supports a model in which the Ca²⁺ channel regulates secretion through a site at the selectivity filter, upstream to cation entry into the cell.     

Kinetic Analysis of Ca2+/K+ Selectivity of an Ion Channel by Single-Binding-Site Models

Current-voltage relationships of a cation channel in the tonoplast of Beta vulgaris, as recorded in solutions with different activities of Ca²⁺ and K⁺ (from Johannes & Sanders 1995, J. Membrane Biol. 146:211–224), have been reevaluated for Ca²⁺/K⁺ selectivity. Since conversion of reversal voltages to permeability ratios by constant field equations is expected to fail because different ions do not move independently through a channel, the data have been analyzed with kinetic channel models instead. Since recent structural information on K⁺ channels show one short and predominant constriction, selectivity models with only one binding site are assumed here to reflect this region kinetically. The rigid-pore model with a main binding site between two energy barriers (nine free parameters) had intrinsic problems to describe the observed current-saturation at large (negative) voltages. The alternative, dynamic-pore model uses a selectivity filter in which the binding site alternates its orientation (empty, or occupied by either Ca²⁺ or K⁺) between the cytoplasmic side and the luminal side within a fraction of the electrical distance and in a rate-limiting fashion. Fits with this model describe the data well. The fits yield about a 10% electrical distance of the selectivity filter, located about 5% more cytoplasmic than the electrical center. For K⁺ translocation, reorientation of the unoccupied binding site (with a preference of about 6:5 to face the lumenal side) is rate limiting. For Ca²⁺, the results show high affinity to the binding site and low translocation rates (<1% of the K⁺ translocation rate). With the fitted model Ca²⁺ entry through the open channel has been calculated for physiological conditions. The model predicts a unitary open channel current of about 100 fA which is insensitive to cytoplasmic Ca²⁺ concentrations (between 0.1 and 1 μm) and which shows little sensitivity to the voltage across the tonoplast. Received: 19 February 1997/Revised: 19 May 1997     

Mechanisms of Beat-to-Beat Regulation of Cardiac Pacemaker Cell Function by Ca2+ Cycling Dynamics

Whether intracellular Ca²⁺ cycling dynamics regulate cardiac pacemaker cell function on a beat-to-beat basis remains unknown. Here we show that under physiological conditions, application of low concentrations of caffeine (2–4 mM) to isolated single rabbit sinoatrial node cells acutely reduces their spontaneous action potential cycle length (CL) and increases Ca²⁺ transient amplitude for several cycles. Numerical simulations, using a modified Maltsev-Lakatta coupled-clock model, faithfully reproduced these effects, and also the effects of CL prolongation and dysrhythmic spontaneous beating (produced by cytosolic Ca²⁺ buffering) and an acute CL reduction (produced by flash-induced Ca²⁺ release from a caged Ca²⁺ buffer), which we had reported previously. Three contemporary numerical models (including the original Maltsev-Lakatta model) failed to reproduce the experimental results. In our proposed new model, Ca²⁺ releases acutely change the CL via activation of the Na⁺/Ca²⁺ exchanger current. Time-dependent CL reductions after flash-induced Ca²⁺ releases (the memory effect) are linked to changes in Ca²⁺ available for pumping into sarcoplasmic reticulum which, in turn, changes the sarcoplasmic reticulum Ca²⁺ load, diastolic Ca²⁺ releases, and Na⁺/Ca²⁺ exchanger current. These results support the idea that Ca²⁺ regulates CL in cardiac pacemaker cells on a beat-to-beat basis, and suggest a more realistic numerical mechanism of this regulation.     

Malignant Hyperthermia: An Inherited Disorder of Skeletal Muscle Ca2+ Regulation

Malignant hyperthermia (MH) is a pharmacogenetic disorder of skeletal muscle characterized by muscle contracture and life-threatening hypermetabolic crisis following exposure to halogenated anesthetics and depolarizing muscle relaxants during surgery. Susceptibility to MH results from mutations in Ca²⁺ channel proteins that mediate excitation–contraction (EC) coupling, with the ryanodine receptor Ca²⁺ release channel (RyR1) representing the major locus. Here we review recent studies characterizing the effects of MH mutations on the sensitivity of the RyR1 to drugs and endogenous channel effectors including Ca²⁺ and calmodulin. In addition, we present a working model that incorporates these effects of MH mutations on the isolated RyR1 with their effects on the physiologic mechanism that activates Ca²⁺ release during EC coupling in intact muscle.     

Calcium-induced potassium pathway in sided erythrocyte membrane vesicles

We have characterized the asymmetric effect of Ca²⁺ on passive K⁺ permeability in erythrocyte membranes, using inside out and right-side out vesicles. Ca²⁺, but not Mg²⁺, can induce an increase in K⁺ uptake in inside out vesicles. The half-maximal concentration of Ca²⁺ required to induce the K⁺ uptake is 0.2 mM, and the permeability increase is not specific for K⁺. Thus, the Ca²⁺-induced permeation process in inside out vesicles is changed from that in the energy-depleted intact cell which requires only micromolar concentrations of Ca²⁺ and is specific for K⁺. Removal of spectrin had no effect on the vesicle permeability increase due to Ca²⁺. Studies with $\text{N-}\text{ethylmaleimide}$ show that the vesicle channel opening is mediated by a protein and passage is controlled by sulfhydryl groups; furthermore, the Ca²⁺-induced vesicle pathway is distinct from the normal channel for passive K⁺ leak in the absence of Ca²⁺. The protein is sensitive to its phospholipid environment since removal of easily accessible phospholipid head groups on the cytoplasmic face of the vesicles inhibits the Ca²⁺-stimulated channel opening.