A theoretical study of calcium microdomains in turtle hair cells.

Confocal imaging has revealed microdomains of intracellular free Ca2+ in turtle hair cells evoked by depolarizing pulses and has delineated factors affecting the growth and dissipation of such domains. However, imaging experiments have limited spatial and temporal resolution. To extend the range of the results we have developed a three-dimensional model of Ca2+ diffusion in a cylindrical hair cell, allowing part of the Ca2+ influx to occur over a small circular region (radius 0.125-1.0 micron) representing a high-density array of voltage-dependent channels. The model incorporated experimental information about the number of channels, the fixed and mobile Ca2+ buffers, and the Ca2+ extrusion mechanism. A feature of the calculations was the use of a variable grid size depending on the proximity to the Ca2+ channel cluster. The results agreed qualitatively with experimental data on the localization of the Ca2+ transients, although the experimental responses were smaller and slower, which is most likely due to temporal and spatial averaging in the imaging. The model made predictions about 1) the optimal Ca2+ channel number and density within a cluster, 2) the conditions to ensure independence of neighboring clusters, and 3) the influence of the Ca2+ buffers on the kinetics and localization of the microdomains. We suggest that an increase in the mobile Ca2+ buffer concentration in high-frequency hair cells (which possess a larger number of release sites) would allow lower amplitude and faster Ca2+ responses and promote functional independence of the sites.     

Analytical steady-state solution to the rapid buffering approximation near an open Ca2+ channel.

We derive an analytical steady-state solution for the Ca2+ profile near an open Ca2+ channel based on a transport equation which describes the buffered diffusion of Ca2+ in the presence of rapid stationary and mobile Ca2+ buffers (Wagner and Keizer, 1994). This steady-state rapid buffering approximation gives an upper bound on local Ca2+ elevations such as Ca2+ puffs or sparks when conditions for the validity of the rapid buffering approximation are met and is an alternative to approximations that assume that mobile buffers are unsaturable. This result also provides an analytical estimate of the cytosolic Ca2+ domain concentration ([Ca2+]d) near a channel pore and shows the dependence of [Ca2+]d on moderate concentrations of endogenous mobile buffer, Ca2+ indicator dye, and bulk cytosolic Ca2+. Assuming a simple relationship between [Ca2+]d and the lumenal depletion domain of an intracellular Ca2+ channel, lumenal and cytosolic Ca2+ profiles are matched to give an implicit analytical expression for the effect of bulk lumenal Ca2+ on [Ca2+]d.     

A simple numerical model of calcium spark formation and detection in cardiac myocytes.

The elementary events of excitation-contraction coupling in heart muscle are Ca2+ sparks, which arise from one or more ryanodine receptors in the sarcoplasmic reticulum (SR). Here a simple numerical model is constructed to explore Ca2+ spark formation, detection, and interpretation in cardiac myocytes. This model includes Ca2+ release, cytosolic diffusion, resequestration by SR Ca2+-ATPases, and the association and dissociation of Ca2+ with endogenous Ca2+-binding sites and a diffusible indicator dye (fluo-3). Simulations in a homogeneous, isotropic cytosol reproduce the brightness and the time course of a typical cardiac Ca2+ spark, but underestimate its spatial size (approximately 1.1 micron vs. approximately 2.0 micron). Back-calculating [Ca2+]i by assuming equilibrium with indicator fails to provide a good estimate of the free Ca2+ concentration even when using blur-free fluorescence data. A parameter sensitivity study reveals that the mobility, kinetics, and concentration of the indicator are essential determinants of the shape of Ca2+ sparks, whereas the stationary buffers and pumps are less influential. Using a geometrically more complex version of the model, we show that the asymmetric shape of Ca2+ sparks is better explained by anisotropic diffusion of Ca2+ ions and indicator dye rather than by subsarcomeric inhomogeneities of the Ca2+ buffer and transport system. In addition, we examine the contribution of off-center confocal sampling to the variance of spark statistics.     

Effects of rapid buffers on Ca2+ diffusion and Ca2+ oscillations.

Based on realistic mechanisms of Ca2+ buffering that include both stationary and mobile buffers, we derive and investigate models of Ca2+ diffusion in the presence of rapid buffers. We obtain a single transport equation for Ca2+ that contains the effects caused by both stationary and mobile buffers. For stationary buffers alone, we obtain an expression for the effective diffusion constant of Ca2+ that depends on local Ca2+ concentrations. Mobile buffers, such as fura-2, BAPTA, or small endogenous proteins, give rise to a transport equation that is no longer strictly diffusive. Calculations are presented to show that these effects can modify greatly the manner and rate at which Ca2+ diffuses in cells, and we compare these results with recent measurements by Allbritton et al. (1992). As a prelude to work on Ca2+ waves, we use a simplified version of our model of the activation and inhibition of the IP3 receptor Ca2+ channel in the ER membrane to illustrate the way in which Ca2+ buffering can affect both the amplitude and existence of Ca2+ oscillations.     

Modeling study of the effects of overlapping Ca2+ microdomains on neurotransmitter release.

Although single-channel Ca2+ microdomains are capable of gating neurotransmitter release in some instances, it is likely that in many cases the microdomains from several open channels overlap to activate vesicle fusion. We describe a mathematical model in which transmitter release is gated by single or overlapping Ca2+ microdomains produced by the opening of nearby Ca2+ channels. This model accounts for the presence of a mobile Ca2+ buffer, provided either that the buffer is unsaturable or that it is saturated near an open channel with Ca2+ binding kinetics that are rapid relative to Ca2+ diffusion. We show that the release time course is unaffected by the location of the channels (at least for distances up to 50 nm), but paired-pulse facilitation is greater when the channels are farther from the release sites. We then develop formulas relating the fractional release following selective or random channel blockage to the cooperative relationship between release and the presynaptic Ca2+ current. These formulas are used with the transmitter release model to study the dependence of this form of cooperativity, which we call Ca2+ current cooperativity, on mobile buffers and on the local geometry of Ca2+ channels. We find that Ca2+ current cooperativity increases with the number of channels per release site, but is considerably less than the number of channels, the theoretical upper bound. In the presence of a saturating mobile buffer the Ca2+ current cooperativity is greater, and it increases more rapidly with the number of channels. Finally, Ca2+ current cooperativity is an increasing function of channel distance, particularly in the presence of saturating mobile buffer.     

Phycobilisome diffusion is required for light-state transitions in cyanobacteria

Phycobilisomes are the major accessory light-harvesting complexes of cyanobacteria and red algae. Studies using fluorescence recovery after photobleaching on cyanobacteria in vivo have shown that the phycobilisomes are mobile complexes that rapidly diffuse on the thylakoid membrane surface. By contrast, the PSII core complexes are completely immobile. This indicates that the association of phycobilisomes with reaction centers must be transient and unstable. Here, we show that when cells of the cyanobacterium Synechococcus sp. PCC7942 are immersed in buffers of high osmotic strength, the diffusion coefficient for the phycobilisomes is greatly decreased. This suggests that the interaction between phycobilisomes and reaction centers becomes much less transient under these conditions. We discuss the possible reasons for this. State transitions are a rapid physiological adaptation mechanism that regulates the way in which absorbed light energy is distributed between PSI and PSII. Immersing cells in high osmotic strength buffers inhibits state transitions by locking cells into whichever state they were in prior to addition of the buffer. The effect on state transitions is induced at the same buffer concentrations as the effect on phycobilisome diffusion. This implies that phycobilisome diffusion is required for state transitions. The main physiological role for phycobilisome mobility may be to allow such flexibility in light harvesting.     

Endogenous calcium buffers regulate fast exocytosis in the synaptic terminal of retinal bipolar cells.

Calcium-triggered exocytosis at the synapse is suppressed by addition of calcium chelators, but the effects of endogenous Ca(2+) buffers have not been tested. We find that 80% of Ca(2+) binding sites in the synaptic terminal of retinal bipolar cells were associated with mobile molecules that suppressed activation of Ca(2+)-sensitive K(+) channels with an efficiency equivalent to approximately 1.2 mM BAPTA. Removing these buffers caused a 30-fold increase in the number of vesicles released by Ca(2+) tail currents lasting approximately 0.5 ms and a 2-fold increase in the rapidly releasable pool of vesicles (RRP). The effects of BAPTA and EGTA indicate that vesicles comprising the RRP were docked at variable distances from Ca(2+) channels. We propose that endogenous Ca(2+) buffers regulate the size of the RRP by suppressing the release of vesicles toward the periphery of the active zone.     

Dynamics of intracellular calcium in hair cells isolated from the semicircular canal of the frog.

Changes in cytosolic free Ca(2+) concentration ([Ca(2+)]i) were monitored optically in hair cells mechanically isolated from frog semicircular canals using the membrane-impermeant form of the Ca(2+)-selective dye Oregon Green 488 BAPTA-1 (OG, 100 microM). Cells stimulated by depolarization under whole-cell voltage clamp conditions revealed Ca(2+) entry at selected sites (hotspots) located mostly in the lower (synaptic) half of the cell body. [Ca(2+)]i at individual hotspots rose with a time constant tau1 approximately 70 ms and decayed with a bi-exponential time-course (tau2 approximately 160, tau3 approximately 2500 ms) following a 160 ms depolarization to -20 mV. With repeated stimulation [Ca(2+)]i underwent independent amplitude changes at distinct hotspots, suggesting that the underlying Ca(2+) channel clusters can be regulated differentially by intracellular signalling pathways. Block by nifedipine indicated that the L-type Ca(2+)channels are distributed at different densities in distinct hotspots. No diffusion barrier other than the nuclear region was found in the cytosol, so that, during a prolonged depolarization (lasting up to 1s), Ca(2+) was able to reach the cell apical ciliated pole. The effective Ca(2+) diffusion constant, measured from the progression of Ca(2+) wavefronts in the cytosol, was approximately 57 microm(2)/s. Our results indicate that in these hair cells, buffered diffusion of Ca(2+) proceeds evenly from the source point to the cell interior and is dominated by the diffusion constant of the endogenous mobile buffers.     

Calcium dependence of bleb formation and cell death in hepatocytes

Calcium dependence of bleb formation and cell death was evaluated in rat hepatocytes following ATP depletion by metabolic inhibition with KCN and iodoacetate ('chemical hypoxia'). Cytosolic free Ca2+ was measured in single cells by ratio imaging of Fura-2 fluorescence using multiparameter digitized video microscopy. Cells formed surface blebs within 10 to 20 minutes after chemical hypoxia and most cells lost viability within an hour. An increase of cytosolic free Ca2+ was not required for bleb formation to occur. One to a few minutes prior to the onset of cell death, free Ca2+ increased rapidly in high Ca2+ buffer (1.2 mM) but not in low Ca2+ buffer (less than 1 microM). In either buffer, the rate of cell killing was the same. As the onset of cell death was approached in both high and low Ca2+ buffers, Fura-2 began to leak from the cells at an accelerating rate indicating rapidly increasing plasma membrane permeability. In high Ca2+ buffer, cytosolic free Ca2+ increased in parallel with dye leakage. No regional changes in cytosolic free Ca2+ were observed during this metastable period of increased membrane permeability. In many experiments, actual rupture of cell surface blebs could be observed which led to micron-size discontinuities of the cell surface and cell death. We conclude that a metastable period characterized by increasing plasma membrane permeability marked the onset of cell death in cultured hepatocytes which culminated in rupture of a cell surface bleb. An increase of cytosolic free Ca2+ was not required for the metastable state to develop or cell death to occur.     

Mitogen-induced oscillations of cytosolic Ca2+ and transmembrane Ca2+ current in human leukemic T cells.

A rapid rise in the level of cytosolic free calcium ([Ca2+]i) is believed to be one of several early triggering signals in the activation of T lymphocytes by antigen. Although Ca2+ release from intracellular stores and its contribution to Ca2+ signaling in many cell types is well documented, relatively little is known regarding the role and mechanism of Ca2+ entry across the plasma membrane. We have investigated mitogen-triggered Ca2+ signaling in individual cells of the human T-leukemia-derived line, Jurkat, using fura-2 imaging and patch-clamp recording techniques. Phytohemagglutinin (PHA), a mitogenic lectin, induces repetitive [Ca2+]i oscillations in these cells peaking at micromolar levels with a period of 90-120 s. The oscillations depend critically upon Ca2+ influx across the plasma membrane, as they are rapidly terminated by removal of extracellular Ca2+, addition of Ca(2+)-channel blockers such as Ni2+ or Cd2+, or membrane depolarization. Whole-cell and perforated-patch recording methods were combined with fura-2 measurements to identify the mitogen-activated Ca2+ conductance involved in this response. A small, highly selective Ca2+ conductance becomes activated spontaneously in whole-cell recordings and in response to PHA in perforated-patch experiments. This conductance has properties consistent with a role in T-cell activation, including activation by PHA, lack of voltage-dependent gating, inhibition by Ni2+ or Cd2+, and regulation by intracellular Ca2+. Moreover, a tight temporal correlation between oscillations of Ca2+ conductance and [Ca2+]i suggests a role for the membrane Ca2+ conductance in generating [Ca2+]i oscillations in activated T cells.     

Facilitation through buffer saturation: constraints on endogenous buffering properties

Synaptic facilitation (SF) is a ubiquitous form of short-term plasticity, regulating synaptic dynamics on fast timescales. Although SF is known to depend on the presynaptic accumulation of Ca(2+), its precise mechanism is still under debate. Recently it has been shown that at certain central synapses SF results at least in part from the progressive saturation of an endogenous Ca(2+) buffer (Blatow et al., 2003), as proposed by Klingauf and Neher (1997). Using computer simulations, we study the magnitude of SF that can be achieved by a buffer saturation mechanism (BSM), and explore its dependence on the endogenous buffering properties. We find that a high SF magnitude can be obtained either by a global saturation of a highly mobile buffer in the entire presynaptic terminal, or a local saturation of a completely immobilized buffer. A characteristic feature of BSM in both cases is that SF magnitude depends nonmonotonically on the buffer concentration. In agreement with results of Blatow et al. (2003), we find that SF grows with increasing distance from the Ca(2+) channel cluster, and increases with increasing external Ca(2+), [Ca(2+)](ext), for small levels of [Ca(2+)](ext). We compare our modeling results with the experimental properties of SF at the crayfish neuromuscular junction, and find that the saturation of an endogenous mobile buffer can explain the observed SF magnitude and its supralinear accumulation time course. However, we show that the BSM predicts slowing of the SF decay rate in the presence of exogenous Ca(2+) buffers, contrary to experimental observations at the crayfish neuromuscular junction. Further modeling and data are required to resolve this aspect of the BSM.     

Simulation of calcium sparks in cut skeletal muscle fibers of the frog

Spark mass, the volume integral of Delta F/F, was investigated theoretically and with simulations. These studies show that the amount of Ca2+ bound to fluo-3 is proportional to mass times the total concentration of fluo-3 ([fluo-3T]); the proportionality constant depends on resting Ca2+ concentration ([Ca2+]R). In the simulation of a Ca2+ spark in an intact frog fiber with [fluo-3T] = 100 microM, fluo-3 captures approximately one-fourth of the Ca2+ released from the sarcoplasmic reticulum (SR). Since mass in cut fibers is several times that in intact fibers, both with similar values of [fluo-3T] and [Ca2+]R, it seems likely that SR Ca2+ release is larger in cut fiber sparks or that fluo-3 is able to capture a larger fraction of the released Ca2+ in cut fibers, perhaps because of reduced intrinsic Ca2+ buffering. Computer simulations were used to identify these and other factors that may underlie the differences in mass and other properties of sparks in intact and cut fibers. Our spark model, which successfully simulates calcium sparks in intact fibers, was modified to reflect the conditions of cut fiber measurements. The results show that, if the protein Ca2+-buffering power of myoplasm is the same as that in intact fibers, the Ca2+ source flux underlying a spark in cut fibers is 5-10 times that in intact fibers. Smaller source fluxes are required for less buffer. In the extreme case in which Ca2+ binding to troponin is zero, the source flux needs to be 3-5 times that in intact fibers. An increased Ca2+ source flux could arise from an increase in Ca2+ flux through one ryanodine receptor (RYR) or an increase in the number of active RYRs per spark, or both. These results indicate that the gating of RYRs, or their apparent single channel Ca2+ flux, is different in frog cut fibers--and, perhaps, in other disrupted preparations--than in intact fibers.     

Functional expression of the human cardiac Na+/Ca2+ exchanger in Sf9 cells: rapid and specific Ni2+ transport

Although inhibition of the Na+/Ca2+ exchanger normally increases [Ca2+]i in neonatal cardiac myocytes, application of the inhibitor Ni2+ appears to reduce [Ca2+] measured by fluo-3. To investigate how the apparent reduction in [Ca2+]i occurs we examined Ca2+ transport by the human Na+/Ca2+ exchanger expressed in Sf9 cells. Transport of Ca2+ by the Na+/Ca2+ exchanger was examined using a laser-scanning confocal microscope and the fluorescent Ca2+ indicator fluo-3, and the electrogenic function was determined by measuring the Na+/Ca2+ exchange current (INaCa) using patch clamp methods. INaCa was elicited with voltage-clamp steps or flash photolysis of caged Ca2+. We show significant expression of Na+/Ca2+ exchanger function in Sf9 cells infected with a recombinant Baculovirus carrying the Na+/Ca2+ exchanger. In addition to measurements of INaCa, characterization includes Ca2+ transport via the Na+/Ca2+ exchanger and the voltage dependence of Ca2+ transport. Application of Ni2+ blocked INaCa but, contrary to expectation, decreased fluo-3 fluorescence. Experiments with infected Sf9 cells suggested that Ni2+ was transported via the Na+/Ca2+ exchanger at a rate comparable to the Ca2+ transport. Once inside the cells, Ni2+ reduced fluorescence, presumably by quenching fluo-3. We conclude that Ni2+ does indeed block INaCa, but is also rapidly translocated across the cell membrane by the Na+/Ca2+ exchanger itself, most likely via an electroneutral partial reaction of the exchange cycle.     

Development and dissipation of Ca(2+) gradients in adrenal chromaffin cells

We used pulsed laser imaging to measure the development and dissipation of Ca(2+) gradients evoked by the activation of voltage-sensitive Ca(2+) channels in adrenal chromaffin cells. Ca(2+) gradients appeared rapidly (<5 ms) upon membrane depolarization and dissipated over several hundred milliseconds after membrane repolarization. Dissipation occurred with an initial fast phase, as the steep gradient near the membrane collapsed, and a slower phase as the remaining shallow gradient dispersed. Inhibition of active Ca(2+) uptake by the endoplasmic reticulum (thapsigargin) and mitochondria (carbonylcyanide p-trifluoro-methoxyphenylhydrazone/oligomycin) had no effect on the size of Ca(2+) changes or the rate of gradient dissipation, suggesting that passive endogenous Ca(2+) buffers are responsible for the slow Ca(2+) redistribution. We used a radial diffusion model incorporating Ca(2+) diffusion and binding to intracellular Ca(2+) buffers to simulate Ca(2+) gradients. We included a 3D optical sectioning model, simulating the effects of out-of-focus light, to allow comparison with the measured gradients. Introduction of a high-capacity immobile Ca(2+) buffer, with a buffer capacity on the order of 1000 and appropriate affinity and kinetics, approximated the size of the Ca(2+) increases and rate of dissipation of the measured gradients. Finally, simulations without exogenous buffer suggest that the Ca(2+) signal due to Ca(2+) channel activation is restricted by the endogenous buffer to a space less than 1 microm from the cell membrane.     

A comparative assessment of fluo Ca2+ indicators in rat ventricular myocytes

The fluo family of indicators is frequently used in studying Ca(2+) physiology; however, choosing which fluo indicator to use is not obvious. Indicator properties are typically determined in well-defined aqueous solutions. Inside cells, however, the properties can change markedly. We have characterized each of three fluo variants (fluo-2MA, fluo-3 and fluo-4) in two forms-the acetoxymethyl (AM) ester and the K(+) salt. We loaded indicators into rat ventricular myocytes and used confocal microscopy to monitor depolarization-induced fluorescence changes and fractional shortening. Myocytes loaded with the indicator AM esters showed significantly different Ca(2+) transients and fractional shortening kinetics. Loading the K(+) salts via whole-cell patch-pipette eliminated differences between fluo-3 and fluo-4, but not fluo-2MA. Cells loaded with different indicator AM esters showed different staining patterns-suggesting differential loading into organelles. Ca(2+) dissociation constants (K(d,Ca)), measured in protein-rich buffers mimicking the cytosol were significantly higher than values determined in simple buffers. This increase in K(d,Ca) (decrease in Ca(2+) affinity) was greatest for fluo-3 and fluo-4, and least for fluo-2MA. We conclude that the structurally-similar fluo variants differ with respect to cellular loading, subcellular compartmentalization, and intracellular Ca(2+) affinity. Therefore, judicious choice of fluo indicator and loading procedure is advisable when designing experiments.     

Calcium microdomains at presynaptic active zones of vertebrate hair cells unmasked by stochastic deconvolution

Signal transduction by auditory and vestibular hair cells involves an impressive ensemble of finely tuned control mechanisms, strictly dependent on the local intracellular Ca(2+) concentration ([Ca(2+)](i)). The study of Ca(2+) dynamics in hair cells typically combines Ca(2+)-sensitive fluorescent indicators (dyes), patch clamp and optical microscopy to produce images of the patterns of fluorescence of a Ca(2+) indicator following various stimulation protocols. Here we describe a novel method that combines electrophysiological recordings, fluorescence imaging and numerical simulations to effectively deconvolve Ca(2+) signals within cytoplasmic microdomains that would otherwise remain inaccessible to direct observation. The method relies on the comparison of experimental data with virtual signals derived from a Monte Carlo reaction-diffusion model based on a realistic reconstruction of the relevant cell boundaries in three dimensions. The model comprises Ca(2+) entry at individual presynaptic active zones followed by diffusion, buffering, extrusion and release of Ca(2+). Our results indicate that changes of the hair cell [Ca(2+)](i) during synaptic transmission are primarily controlled by the Ca(2+) endogenous buffers both at short (<1mu) and at long (tens of microns) distances from the active zones. We provide quantitative estimates of concentration and kinetics of the hair cell endogenous Ca(2+) buffers and Ca(2+)-ATPases. We finally show that experimental fluorescence data collected during Ca(2+) influx are not interpreted correctly if the [Ca(2+)](i) is estimated by assuming that Ca(2+) equilibrates instantly with its reactants. In our opinion, this approach is of potentially general interest as it can be easily adapted to the study of Ca(2+) dynamics in diverse biological systems.     

The spatial relationship between Ca2+ channels and Ca2+-activated channels and the function of Ca2+-buffering in avian sensory neurons

In order to learn about the endogenous Ca2+-buffering in the cytoplasm of chick dorsal root ganglion (DRG) neurons and the distance separating the ryanodine receptor Ca2+ release channels (RyRs) from the plasma membrane, we monitored the amplitude and time course of Ca2+-activated Cl- currents (I(ClCa)) in protocols that manipulated Ca2+-buffering. I(ClCa)was activated by Ca2+ influx via voltage-gated Ca2+ channels or by Ca2+ release via RyRs activated by 10 mM caffeine. I(ClCa)was measured in neurons at 20 degrees C and 35 degrees C using the amphotericin perforated patch technique that preserves endogenous Ca2+-buffering, or at 20 degrees C in neurons dialyzed with pipette solutions designed to replace the endogenous Ca2+ buffers. The amplitude of I(ClCa)activated by Ca2+ influx or Ca2+ at 20 degrees C was similar in the amphotericin neurons and neurons dialyzed with an 'unbuffered' pipette solution containing 10 mM citrate and 3 mM ATP as the only Ca2+ binding molecules. Thus, endogenous mobile Ca2+ buffers are relatively unimportant in chick DRG neurons. Warming the neurons from 20 degrees C to 35 degrees C increased the amplitude and the rate of deactivation of I(ClCa)consistent with an increased rate of Ca2+ buffering by fixed endogenous Ca2+-buffers. Dialysis with 2 mM EGTA/0.1 microM free Ca2+ reduced the amplitude and increased the rate of deactivation of I(ClCa)activated by Ca2+ influx and abolished I(ClCa)activated by Ca2+ release. Dialysis with 2 mM BAPTA/0.1 microM free Ca2+ abolished I(ClCa)activated by Ca2+ influx or release. Dialysis with 42 mM HEEDTA/0.5 microM free Ca2+ caused the persistent activation of I(ClCa). Calculations using a Ca2+-diffusion model suggest that the voltage-gated Ca2+ channels and the Ca2+-activated Cl- channels are separated by 50-400 nm and that the RyRs are more than 600 nm from the plasma membrane.     

On the roles of Ca2+ diffusion, Ca2+ buffers, and the endoplasmic reticulum in IP3-induced Ca2+ waves.

We have investigated the effects of Ca2+ diffusion, mobile and stationary Ca2+ buffers in the cytosol, and Ca2+ handling by the endoplasmic reticulum on inositol 1,4,5-trisphosphate-induced Ca2+ wave propagation. Rapid equilibration of free and bound Ca2+ is used to describe Ca2+ sequestration by buffers in both the cytosol and endoplasmic reticulum (ER) lumen. Cytosolic Ca2+ regulation is based on a kinetic model of the inositol 1,4,5-trisphosphate (IP3) receptor of De Young and Keizer that includes activation and inhibition of the IP3 receptor Ca2+ channel in the ER membrane and SERCA Ca2+ pumps in the ER. Diffusion of Ca2+ in the cytosol and the ER and the breakdown and diffusion of IP3 are also included in our calculations. Although Ca2+ diffusion is severely limited because of buffering, when conditions are chosen just below the threshold for Ca2+ oscillations, a pulse of IP3 or Ca2+ results in a solitary trigger wave that requires diffusion of Ca2+ for its propagation. In the oscillatory regime repetitive wave trains are observed, but for this type of wave neither the wave shape nor the speed is strongly dependent on the diffusion of Ca2+. Local phase differences lead to waves that are predominately kinematic in nature, so that the wave speed (c) is related to the wavelength (lambda) and the period of the oscillations (tau) approximately by the formula c = lambda/tau. The period is determined by features that control the oscillations, including [IP3] and pump activity, which are related to recent experiments. Both solitary waves and wave trains are accompanied by a Ca2+ depletion wave in the ER lumen, similar to that observed in cortical preparations from sea urchin eggs. We explore the effect of endogenous and exogenous Ca2+ buffers on wave speed and wave shape, which can be explained in terms of three distinct effects of buffering, and show that exogenous buffers or Ca2+ dyes can have considerable influence on the amplitude and width of the waves.     

Intracellular proton mobility and buffering power in cardiac ventricular myocytes from rat,rabbit, and guinea pig

Intracellular pH (pHi) is an important modulator of cardiac function. The spatial regulation of pH within the cytoplasm depends, in part, on intracellular H+ (Hi+) mobility. The apparent diffusion coefficient for Hi+, DHapp, was estimated in single ventricular myocytes isolated from the rat, guinea pig, and rabbit. DHapp was derived by best-fitting predictions of a two-dimensional model of H+ diffusion to the local rise of intracellular [H+], recorded confocally (ratiometric seminaphthorhodafluor fluorescence) downstream from an acid-filled, whole cell patch pipette. Under CO2/HCO3--free conditions, DHapp was similar in all three species (mean values: 8-12.5 x 10-7 cm2/s) and was over 200-fold lower than that for H+ in water. In guinea pig myocytes, DHapp was increased 2.5-fold in the presence of CO2/HCO3- buffer, in agreement with previous observations in rabbit myocytes. Hi+ mobility is therefore low in cardiac cells, a feature that may predispose them to the generation of pHi gradients in response to sarcolemmal acid/base transport or local cytoplasmic acid production. Low Hi+ mobility most likely results from H+ shuttling among cytoplasmic mobile and fixed buffers. This hypothesis was explored by comparing the pHi dependence of intrinsic, intracellular buffering capacity, measured for all three species, and subdividing buffering into mobile and fixed fractions. The proportion of buffer that is mobile will be the main determinant of DHapp. At a given pHi, this proportion appeared to be similar in all three species, consistent with a common value for DHapp. Over the pHi range of 6.0-8.0, the proportion is expected to change, predicting that DHapp may display some pHi sensitivity.     

Endogenous and exogenous Ca2+ buffers differentially modulate Ca2+-dependent inactivation of Ca(v)2.1 Ca2+ channels

Voltage-gated Ca2+ channels undergo a negative feedback regulation by Ca2+ ions, Ca2+-dependent inactivation, which is important for restricting Ca2+ signals in nerve and muscle. Although the molecular details underlying Ca2+-dependent inactivation have been characterized, little is known about how this process might be modulated in excitable cells. Based on previous findings that Ca2+-dependent inactivation of Ca(v)2.1 (P/Q-type) Ca2+ channels is suppressed by strong cytoplasmic Ca2+ buffering, we investigated how factors that regulate cellular Ca2+ levels affect inactivation of Ca(v)2.1 Ca2+ currents in transfected 293T cells. We found that inactivation of Ca(v)2.1 Ca2+ currents increased exponentially with current amplitude with low intracellular concentrations of the slow buffer EGTA (0.5 mm), but not with high concentrations of the fast Ca2+ buffer BAPTA (10 mm). However, when the concentration of BAPTA was reduced to 0.5 mm, inactivation of Ca2+ currents was significantly greater than with an equivalent concentration of EGTA, indicating the importance of buffer kinetics in modulating Ca2+-dependent inactivation of Ca(v)2.1. Cotransfection of Ca(v)2.1 with the EF-hand Ca2+-binding proteins, parvalbumin and calbindin, significantly altered the relationship between Ca2+ current amplitude and inactivation in ways that were unexpected from behavior as passive Ca2+ buffers. We conclude that Ca2+-dependent inactivation of Ca(v)2.1 depends on a subplasmalemmal Ca2+ microdomain that is affected by the amplitude of the Ca2+ current and differentially modulated by distinct Ca2+ buffers.