A mathematical model of cardiocyte Ca(2+) dynamics with a novel representation of sarcoplasmic reticular Ca(2+) control

Facilitation by intracellular carbonic anhydrase of Na+ -HCO3- co-transport but not Na+ / H+ exchange activity in the mammalian ventricular myocyte

Carbonic anhydrase enzymes (CAs) catalyse the reversible hydration of CO2 to H+ and HCO3- ions. This catalysis is proposed to be harnessed by acid/base transporters, to facilitate their transmembrane flux activity, either through direct protein-protein binding (a 'transport metabolon') or local functional interaction. Flux facilitation has previously been investigated by heterologous co-expression of relevant proteins in host cell lines/oocytes. Here, we examine the influence of intrinsic CA activity on membrane HCO3- or H+ transport via the native acid-extruding proteins, Na+ -HCO3- cotransport (NBC) and Na+ / H+ exchange (NHE), expressed in enzymically isolated mammalian ventricular myocytes. Effects of intracellular and extracellular (exofacial) CA (CAi and CAe) are distinguished using membrane-permeant and -impermeant pharmacological CA inhibitors, while measuring transporter activity in the intact cell using pH and Na+ fluorophores. We find that NBC, but not NHE flux is enhanced by catalytic CA activity, with facilitation being confined to CAi activity alone. Results are quantitatively consistent with a model where CAi catalyses local H+ ion delivery to the NBC protein, assisting the subsequent (uncatalysed) protonation and removal of imported HCO3- ions. In well-superfused myocytes, exofacial CA activity is superfluous, most likely because extracellular CO2/HCO3- buffer is clamped at equilibrium. The CAi insensitivity of NHE flux suggests that, in the native cell, intrinsic mobile buffer-shuttles supply sufficient intracellular H+ ions to this transporter, while intrinsic buffer access to NBC proteins is restricted. Our results demonstrate a selective CA facilitation of acid/base transporters in the ventricular myocyte, implying a specific role for the intracellular enzyme in HCO3- transport, and hence pHi regulation in the heart.

Models of calcium dynamics in cerebellar granule cells

Intracellular calcium dynamics is critical for many functions of cerebellar granule cells (GrCs) including membrane excitability, synaptic plasticity, apoptosis, and regulation of gene transcription. Recent measurements of calcium responses in GrCs to depolarization and synaptic stimulation reveal spatial compartmentalization and heterogeneity within dendrites of these cells. However, the main determinants of local calcium dynamics in GrCs are still poorly understood. One reason is that there have been few published studies of calcium dynamics in intact GrCs in their native environment. In the absence of complete information, biophysically realistic models are useful for testing whether specific Ca(2+) handling mechanisms may account for existing experimental observations. Simulation results can be used to identify critical measurements that would discriminate between different models. In this review, we briefly describe experimental studies and phenomenological models of Ca(2+) signaling in GrC, and then discuss a particular biophysical model, with a special emphasis on an approach for obtaining information regarding the distribution of Ca(2+) handling systems under conditions of incomplete experimental data. Use of this approach suggests that Ca(2+) channels and fixed endogenous Ca(2+) buffers are highly heterogeneously distributed in GrCs. Research avenues for investigating calcium dynamics in GrCs by a combination of experimental and modeling studies are proposed.

Modeling calcium waves in cardiac myocytes: importance of calcium diffusion

Under certain conditions, cardiac myocytes engage in a mode of calcium signaling in which calcium release from the sarcoplasmic reticulum (SR) to myoplasm occurs in self-propagating succession along the length of the cell. This event is called a calcium wave and is fundamentally a diffusion-reaction phenomenon. We present a simple, continuum mathematical model that simulates calcium waves. The framework features calcium diffusion within the SR and myoplasm, and dual modulation of ryanodine receptor (RyR) release channels by myoplasmic and SR calcium. The model is used to illustrate the effect of varying RyR permeability, sarco/endoplasmic reticulum Ca2+ ATPase (SERCA) activity and calcium ion mobility in myoplasm and SR on wave velocity. The model successfully reproduces calcium waves using experimentally-derived variables. It also supports the proposal for wave propagation driven by the diffusive spread of myoplasmic calcium, and highlights the importance of SR calcium load on wave propagation.

Partial potentials of selected cardiac muscle regions and heart activity model based on single fibres

We present single fibre heart activity model (SFHAM) based on the current flow through the five bunches of fibres of the cardiac muscle (CM). The five effective fibres are identified and assigned to the appropriate segments of CM. Analytical functions describing ionic flows along the fibres are derived and proposed. The parameters determining the shapes and amplitudes of the functions proposed are obtained on the basis of standard 12-lead ECG measurements after numerical fitting procedures concentrating on the QRS-waves. As a consequence, five independent courses of partial, transient potentials are obtained representing: anterior, inferior, lateral, posterior walls, and interventricular septum activities, respectively. Moreover, to check our theoretical results we compare the potentials calculated with those from physical measurements performed on the patient's body surface. We expect that SFHAM will permit detection of pathological changes in particular fragments of CM.

Calsequestrin-mediated mechanism for cellular calcium transient alternans

Intracellular calcium transient alternans (CTA) has a recognized role in arrhythmogenesis, but its origin is not yet fully understood. Recent models of CTA are based on a steep relationship between calcium release from the sarcoplasmic reticulum (SR) and its calcium load before release. This mechanism alone, however, does not explain recent observations of CTA without diastolic SR calcium content alternations. In addition, nanoscopic imaging of calcium dynamics has revealed that the elementary calcium release units of the SR can become refractory independently of their local calcium content. Here we show using a new physiologically detailed mathematical model of calcium cycling that luminal gating of the calcium release channels (RyRs) mediated by the luminal buffer calsequestrin (CSQN) can cause CTA independently of the steepness of the release-load relationship. In this complementary mechanism, CTA is caused by a beat-to-beat alternation in the number of refractory RyR channels and can occur with or without diastolic SR calcium content alternans depending on pacing conditions and uptake dynamics. The model has unique features, in that it treats a realistic number of spatially distributed and diffusively coupled dyads, each one with a realistic number of RyR channels, and that luminal CSQN buffering and gating is incorporated based on experimental data that characterizes the effect of the conformational state of CSQN on its buffering properties. In addition to reproducing observed features of CTA, this multiscale model is able to describe recent experiments in which CSQN expression levels were genetically altered as well as to reproduce nanoscopic measurements of spark restitution properties. The ability to link microscopic properties of the calcium release units to whole cell behavior makes this model a powerful tool to investigate the arrhythmogenic role of abnormal calcium handling in many pathological settings.

Alternans of cardiac calcium cycling in a cluster of ryanodine receptors: a simulation study

Mechanical alternans in cardiac muscle is associated with intracellular Ca(2+) alternans. Mechanisms underlying intracellular Ca(2+) alternans are unclear. In previous experimental studies, we produced alternans of systolic Ca(2+) under voltage clamp, either by partially inhibiting the Ca(2+) release mechanism, or by applying small depolarizing pulses. In each case, alternans relied on propagating waves of Ca(2+) release. The aim of this study is to investigate by computer modeling how alternans of systolic Ca(2+) is produced. A mathematical model of a cardiac cell with 75 coupled elements is developed, with each element contains L-type Ca(2+) current, a subspace into which Ca release takes place, a cytoplasmic space, sarcoplasmic reticulum (SR) release channels [ryanodine receptor (RyR)], and uptake sites (SERCA). Interelement coupling is via Ca(2+) diffusion between neighboring subspaces via cytoplasmic spaces and network SR spaces. Small depolarizing pulses were simulated by step changes of cell membrane potential (20 mV) with random block of L-type channels. Partial inhibition of the release mechanism is mimicked by applying a reduction of RyR open probability in response to full stimulation by L-type channels. In both cases, systolic alternans follow, consistent with our experimental observations, being generated by propagating waves of Ca(2+) release and sustained through alternation of SR Ca(2+) content. This study provides novel and fundamental insights to understand mechanisms that may underlie intracellular Ca(2+) alternans without the need for refractoriness of L-type Ca or RyR channels under rapid pacing.

Moment closure for local control models of calcium-induced calcium release in cardiac myocytes

In prior work, we introduced a probability density approach to modeling local control of Ca2+-induced Ca2+ release in cardiac myocytes, where we derived coupled advection-reaction equations for the time-dependent bivariate probability density of subsarcolemmal subspace and junctional sarcoplasmic reticulum (SR) [Ca2+] conditioned on Ca2+ release unit (CaRU) state. When coupled to ordinary differential equations (ODEs) for the bulk myoplasmic and network SR [Ca2+], a realistic but minimal model of cardiac excitation-contraction coupling was produced that avoids the computationally demanding task of resolving spatial aspects of global Ca2+ signaling, while accurately representing heterogeneous local Ca2+ signals in a population of diadic subspaces and junctional SR depletion domains. Here we introduce a computationally efficient method for simulating such whole cell models when the dynamics of subspace [Ca2+] are much faster than those of junctional SR [Ca2+]. The method begins with the derivation of a system of ODEs describing the time-evolution of the moments of the univariate probability density functions for junctional SR [Ca2+] jointly distributed with CaRU state. This open system of ODEs is then closed using an algebraic relationship that expresses the third moment of junctional SR [Ca2+] in terms of the first and second moments. In simulated voltage-clamp protocols using 12-state CaRUs that respond to the dynamics of both subspace and junctional SR [Ca2+], this moment-closure approach to simulating local control of excitation-contraction coupling produces high-gain Ca2+ release that is graded with changes in membrane potential, a phenomenon not exhibited by common pool models. Benchmark simulations indicate that the moment-closure approach is nearly 10,000-times more computationally efficient than corresponding Monte Carlo simulations while leading to nearly identical results. We conclude by applying the moment-closure approach to study the restitution of Ca2+-induced Ca2+ release during simulated two-pulse voltage-clamp protocols.

A calcium-induced calcium release mechanism mediated by calsequestrin

Calcium (Ca(2+))-induced Ca(2+) release (CICR) is widely accepted as the principal mechanism linking electrical excitation and mechanical contraction in cardiac cells. The CICR mechanism has been understood mainly based on binding of cytosolic Ca(2+) with ryanodine receptors (RyRs) and inducing Ca(2+) release from the sarcoplasmic reticulum (SR). However, recent experiments suggest that SR lumenal Ca(2+) may also participate in regulating RyR gating through calsequestrin (CSQ), the SR lumenal Ca(2+) buffer. We investigate how SR Ca(2+) release via RyR is regulated by Ca(2+) and calsequestrin (CSQ). First, a mathematical model of RyR kinetics is derived based on experimental evidence. We assume that the RyR has three binding sites, two cytosolic sites for Ca(2+) activation and inactivation, and one SR lumenal site for CSQ binding. The open probability (P(o)) of the RyR is found by simulation under controlled cytosolic and SR lumenal Ca(2+). Both peak and steady-state P(o) effectively increase as SR lumenal Ca(2+) increases. Second, we incorporate the RyR model into a CICR model that has both a diadic space and the junctional SR (jSR). At low jSR Ca(2+) loads, CSQs are more likely to bind with the RyR and act to inhibit jSR Ca(2+) release, while at high SR loads CSQs are more likely to detach from the RyR, thereby increasing jSR Ca(2+) release. Furthermore, this CICR model produces a nonlinear relationship between fractional jSR Ca(2+) release and jSR load. These findings agree with experimental observations in lipid bilayers and cardiac myocytes.

Ca2+-mobility in the sarcoplasmic reticulum of ventricular myocytes is low

The sarcoplasmic reticulum (SR) in ventricular myocytes contains releasable Ca(2+) for activating cellular contraction. Recent measurements of intra-SR (luminal) Ca(2+) suggest a high diffusive Ca(2+)-mobility constant (D(CaSR)). This could help spatially to unify SR Ca(2+)-content ([Ca(2+)](SRT)) and standardize Ca(2+)-release throughout the cell. But measurements of localized depletions of luminal Ca(2+) (Ca(2+)-blinks), associated with local Ca(2+)-release (Ca(2+)-sparks), suggest D(CaSR) may actually be low. Here we describe a novel method for measuring D(CaSR). Using a cytoplasmic Ca(2+)-fluorophore, we estimate regional [Ca(2+)](SRT) from localized, caffeine-induced SR Ca(2+)-release. Caffeine microperfusion of one end of a guinea pig or rat myocyte diffusively empties the whole SR at a rate indicating D(CaSR) is 8-9 microm(2)/s, up to tenfold lower than previous estimates. Ignoring background SR Ca(2+)-leakage in our measurement protocol produces an artifactually high D(CaSR) (>40 microm(2)/s), which may also explain the previous high values. Diffusion-reaction modeling suggests that a low D(CaSR) would be sufficient to support local SR Ca(2+)-signaling within sarcomeres during excitation-contraction coupling. Low D(CaSR) also implies that [Ca(2+)](SRT) may readily become spatially nonuniform, particularly under pathological conditions of spatially nonuniform Ca(2+)-release. Local control of luminal Ca(2+), imposed by low D(CaSR), may complement the well-established local control of SR Ca(2+)-release by Ca(2+)-channel/ryanodine receptor couplons.

Measurement and modeling of Ca2+ waves in isolated rabbit ventricular cardiomyocytes

The time course and magnitude of the Ca(2+) fluxes underlying spontaneous Ca(2+) waves in single permeabilized ventricular cardiomyocytes were derived from confocal Fluo-5F fluorescence signals. Peak flux rates via the sarcoplasmic reticulum (SR) release channel (RyR2) and the SR Ca(2+) ATPase (SERCA) were not constant across a range of cellular [Ca(2+)] values. The Ca(2+) affinity (K(mf)) and maximum turnover rate (V(max)) of SERCA and the peak permeability of the RyR2-mediated Ca(2+) release pathway increased at higher cellular [Ca(2+)] loads. This information was used to create a computational model of the Ca(2+) wave, which predicted the time course and frequency dependence of Ca(2+) waves over a range of cellular Ca(2+) loads. Incubation of cardiomyocytes with the Ca(2+) calmodulin (CaM) kinase inhibitor autocamtide-2-related inhibitory peptide (300 nM, 30 mins) significantly reduced the frequency of the Ca(2+) waves at high Ca(2+) loads. Analysis of the Ca(2+) fluxes suggests that inhibition of CaM kinase prevented the increases in SERCA V(max) and peak RyR2 release flux observed at high cellular [Ca(2+)]. These data support the view that modification of activity of SERCA and RyR2 via a CaM kinase sensitive process occurs at higher cellular Ca(2+) loads to increase the maximum frequency of spontaneous Ca(2+) waves.

Modelling calcium microdomains using homogenisation

Microdomains of calcium (i.e., areas on the nanometer scale that have qualitatively different calcium concentrations from that in the bulk cytosol) are known to be important in many situations. In cardiac cells, for instance, a calcium microdomain between the L-type channels and the ryanodine receptors, the so-called diadic cleft, is where the majority of the control of calcium release occurs. In other cell types that exhibit calcium oscillations and waves, the importance of microdomains in the vicinity of clusters of inositol trisphosphate receptors, or between the endoplasmic reticulum (ER) and other internal organelles or the plasma membrane, is clear. Given the limits of computational power, it is not currently realistic to model an entire cellular cytoplasm by incorporating detailed structural information about the ER throughout the entire cytoplasm. Hence, most models use a homogenised approach, assuming that both cytoplasm and ER coexist at each point of the domain. Conversely, microdomain models can be constructed, in which detailed structural information can be incorporated, but, until now, methods have not been developed for linking such a microdomain model to a model at the level of the entire cell. Using the homogenisation approach we developed in an earlier paper [Goel, P., Friedman, A., Sneyd, J., 2006. Homogenization of the cell cytoplasm: the calcium bidomain equations. SIAM J. Multiscale Modeling Simulation, in press] we show how a multiscale model of a calcium microdomain can be constructed. In this model a detailed model of the microdomain (in which the ER and the cytoplasm are separate compartments) is coupled to a homogenised model of the entire cell in a rigorous way. Our method is illustrated by a simple model of the diadic cleft of a cardiac half-sarcomere.

A probability density approach to modeling local control of calcium-induced calcium release in cardiac myocytes

We present a probability density approach to modeling localized Ca2+ influx via L-type Ca2+ channels and Ca2+-induced Ca2+ release mediated by clusters of ryanodine receptors during excitation-contraction coupling in cardiac myocytes. Coupled advection-reaction equations are derived relating the time-dependent probability density of subsarcolemmal subspace and junctional sarcoplasmic reticulum [Ca2+] conditioned on "Ca2+ release unit" state. When these equations are solved numerically using a high-resolution finite difference scheme and the resulting probability densities are coupled to ordinary differential equations for the bulk myoplasmic and sarcoplasmic reticulum [Ca2+], a realistic but minimal model of cardiac excitation-contraction coupling is produced. Modeling Ca2+ release unit activity using this probability density approach avoids the computationally demanding task of resolving spatial aspects of global Ca2+ signaling, while accurately representing heterogeneous local Ca2+ signals in a population of diadic subspaces and junctional sarcoplasmic reticulum depletion domains. The probability density approach is validated for a physiologically realistic number of Ca2+ release units and benchmarked for computational efficiency by comparison to traditional Monte Carlo simulations. In simulated voltage-clamp protocols, both the probability density and Monte Carlo approaches to modeling local control of excitation-contraction coupling produce high-gain Ca2+ release that is graded with changes in membrane potential, a phenomenon not exhibited by so-called "common pool" models. However, a probability density calculation can be significantly faster than the corresponding Monte Carlo simulation, especially when cellular parameters are such that diadic subspace [Ca2+] is in quasistatic equilibrium with junctional sarcoplasmic reticulum [Ca2+] and, consequently, univariate rather than multivariate probability densities may be employed.

Contribution of the Na+/Ca2+ exchanger to rapid Ca2+ release in cardiomyocytes

Trigger Ca(2+) is considered to be the Ca(2+) current through the L-type Ca(2+) channel (LTCC) that causes release of Ca(2+) from the sarcoplasmic reticulum. However, cell contraction also occurs in the absence of the LTCC current (I(Ca)). In this article, we investigate the contribution of the Na(+)/Ca(2+) exchanger (NCX) to the trigger Ca(2+). Experimental data from rat cardiomyocytes using confocal microscopy indicating that inhibition of reverse mode Na(+)/Ca(2+) exchange delays the Ca(2+) transient by 3-4 ms served as a basis for the mathematical model. A detailed computational model of the dyadic cleft (fuzzy space) is presented where the diffusion of both Na(+) and Ca(2+) is taken into account. Ionic channels are included at discrete locations, making it possible to study the effect of channel position and colocalization. The simulations indicate that if a Na(+) channel is present in the fuzzy space, the NCX is able to bring enough Ca(2+) into the cell to affect the timing of release. However, this critically depends on channel placement and local diffusion properties. With fuzzy space diffusion in the order of four orders of magnitude lower than in water, triggering through LTCC alone was up to 5 ms slower than with the presence of a Na(+) channel and NCX.

Modulatory effect of calmodulin-dependent kinase II (CaMKII) on sarcoplasmic reticulum Ca2+ handling and interval-force relations: a modelling study

We hypothesize that slow inactivation of Ca2+/calmodulin-dependent kinase II (CaMKII) and its modulatory effect on sarcoplasmic reticulum (SR) Ca2+ handling are important for various interval-force (I-F) relations, in particular for the beat interval dependency in transient alternans during the decay of post-extrasystolic potentiation. We have developed a mathematical model of a single cardiomyocyte to integrate various I-F relations, including alternans, by incorporating a conceptual CaMKII kinetics model into the SR Ca2+ handling model. Our model integrates I-F relations, such as the beat interval-dependent twitch force duration, restitution and potentiation, positive staircase phenomenon and alternans. We found that CaMKII affects more or less all I-F relations, and it is a key factor for integration of the various I-F relations in our model. Alternans arises, in the model, out of a steep relation between SR Ca2+ load and release, owing to SR load-dependent changes in the releasability of Ca2+ via the ryanodine receptor. Beat interval-dependent CaMKII activity, owing to its kinetic properties and amplifying effect on SR Ca2+ load dependency of Ca2+ release, replicated the beat interval dependency of alternans, as observed experimentally. Additionally, our model enabled reproduction of the effects of various interventions on alternans, such as the slowing or accelerating of Ca2+ release and/or uptake. We conclude that a slow time-dependent factor, represented in the model by CaMKII, is important for the integration of I-F relations, including alternans, and that our model offers a useful tool for further analysis of the roles of integrative Ca2+ handling in myocardial I-F relations.

Model of intracellular calcium cycling in ventricular myocytes

We present a mathematical model of calcium cycling that takes into account the spatially localized nature of release events that correspond to experimentally observed calcium sparks. This model naturally incorporates graded release by making the rate at which calcium sparks are recruited proportional to the whole cell L-type calcium current, with the total release of calcium from the sarcoplasmic reticulum (SR) being just the sum of local releases. The dynamics of calcium cycling is studied by pacing the model with a clamped action potential waveform. Experimentally observed calcium alternans are obtained at high pacing rates. The results show that the underlying mechanism for this phenomenon is a steep nonlinear dependence of the calcium released from the SR on the diastolic SR calcium concentration (SR load) and/or the diastolic calcium level in the cytosol, where the dependence on diastolic calcium is due to calcium-induced inactivation of the L-type calcium current. In addition, the results reveal that the calcium dynamics can become chaotic even though the voltage pacing is periodic. We reduce the equations of the model to a two-dimensional discrete map that relates the SR and cytosolic concentrations at one beat and the previous beat. From this map, we obtain a condition for the onset of calcium alternans in terms of the slopes of the release-versus-SR load and release-versus-diastolic-calcium curves. From an analysis of this map, we also obtain an understanding of the origin of chaotic dynamics.

Analysing cardiac excitation–contraction coupling with mathematical models of local control

Cardiac excitation-contraction (E-C) coupling describes the process that links sarcolemmal Ca2+ influx via L-type Ca2+ channels to Ca2+ release from the sarcoplasmic reticulum via ryanodine receptors (RyRs). This process has proven difficult to study experimentally, and complete descriptions of how the cell couples surface membrane and intracellular signal transduction proteins to achieve both stable and sensitive intracellular calcium release are still lacking. Mathematical models provide a framework to test our understanding of how this is achieved. While no single model is yet capable of describing all features of cardiac E-C coupling, models of increasing complexity are revealing unexpected subtlety in the process. In particular, modelling has established a general failure of 'common-pool' models and has emphasized the requirement for 'local control' so that microscopic sub-cellular domains can separate local behaviour from the whole-cell average (common-pool) behaviour. The micro-architecture of the narrow diadic cleft in which the local control takes place is a key factor in determining local Ca2+ dynamics. There is still considerable uncertainty about the number of Ca2+ ions required to open RyRs within the cleft and various gating models have been proposed, many of which are in reasonable agreement with available experimental data. However, not all models exhibit a realistic voltage dependence of E-C coupling gain. Furthermore, it is unclear which model features are essential to producing reasonable gain properties. Thus, despite the success of local-control models in explaining many features of cardiac E-C coupling, more work will be needed to provide a sound theoretical basis of cardiac E-C coupling.

Dynamical description of sinoatrial node pacemaking: improved mathematical model for primary pacemaker cell

We developed an improved mathematical model for a single primary pacemaker cell of the rabbit sinoatrial node. Original features of our model include 1) incorporation of the sustained inward current (I(st)) recently identified in primary pacemaker cells, 2) reformulation of voltage- and Ca(2+)-dependent inactivation of the L-type Ca(2+) channel current (I(Ca,L)), 3) new expressions for activation kinetics of the rapidly activating delayed rectifier K(+) channel current (I(Kr)), and 4) incorporation of the subsarcolemmal space as a diffusion barrier for Ca(2+). We compared the simulated dynamics of our model with those of previous models, as well as with experimental data, and examined whether the models could accurately simulate the effects of modulating sarcolemmal ionic currents or intracellular Ca(2+) dynamics on pacemaker activity. Our model represents significant improvements over the previous models, because it can 1) simulate whole cell voltage-clamp data for I(Ca,L), I(Kr), and I(st); 2) reproduce the waveshapes of spontaneous action potentials and ionic currents during action potential clamp recordings; and 3) mimic the effects of channel blockers or Ca(2+) buffers on pacemaker activity more accurately than the previous models.

Cardiac nitric oxide synthase 1 regulates basal and beta-adrenergic contractility in murine ventricular myocytes

BACKGROUND

Evidence indicates that myocardial NO production can modulate contractility, but the source of NO remains uncertain. Here, we investigated the role of a type 1 NO synthase isoform (NOS1), which has been recently localized to the cardiac sarcoplasmic reticulum, in the regulation of basal and beta-adrenergic myocardial contraction.

METHODS AND RESULTS

Contraction was assessed in left ventricular myocytes isolated from mice with NOS1 gene disruption (NOS1(-/-) mice) and their littermate controls (NOS1(+/+) mice) at 3 stimulation frequencies (1, 3, and 6 Hz) in basal conditions and during beta-adrenergic stimulation with isoproterenol (2 nmol/L). In addition, we examined the effects of acute specific inhibition of NOS1 with vinyl-L-N-5-(1-imino-3-butenyl)-L-ornithine (L-VNIO, 500 micromol/L). NOS1((-/-)) myocytes exhibited greater contraction at all frequencies (percent cell shortening at 6 Hz, 10.7+/-0.92% in NOS1(-/-) myocytes versus 7.21+/-0.8% in NOS1(+/+) myocytes; P<0.05) with a flat frequency-contraction relationship. Time to 50% relaxation was increased in NOS1(-/-) myocytes at all frequencies (at 6 Hz, 26.53+/-1.4 ms in NOS1(-/-) myocytes versus 21.27+/-1.3 ms in NOS1(+/+) myocytes; P<0.05). L-VNIO prolonged time to 50% relaxation at all frequencies (at 6 Hz, 21.28+/-1.7 ms in NOS1(+/+) myocytes versus 26.45+/-1.4 ms in NOS1(+/+)+L-VNIO myocytes; P<0.05) but did not significantly increase basal contraction. However, both NOS1(-/-) myocytes and NOS1(+/+) myocytes treated with L-VNIO showed a greatly enhanced contraction in response to beta-adrenergic stimulation (percent increase in contraction at 6 Hz, 25.2+/-10.8 in NOS1(+/+) myocytes, 68.2+/-11.2 in NOS1(-/-) myocytes, and 65.1+/-13.2 in NOS1(+/+)+L-VNIO myocytes; P<0.05).

CONCLUSIONS

NOS1 disruption enhances basal contraction and the inotropic response to beta-adrenergic stimulation in murine ventricular myocytes. These findings indicate that cardiac NOS1-derived NO plays a significant role in the autocrine regulation of myocardial contractility.

Role of creatine kinase in cardiac excitation-contraction coupling: studies in creatine kinase-deficient mice

To understand the role of creatine kinase (CK) in cardiac excitation-contraction coupling, CK-deficient mice (CK-/-) were studied in vitro and in vivo. In skinned fibers, the kinetics of caffeine-induced release of Ca2+ was markedly slowed in CK-/- mice with a partial restoration when glycolytic substrates were added. These abnormalities were almost compensated for at the cellular level: the responses of Ca2+ transient and cell shortening to an increased pacing rate from 1 Hz to 4 Hz were normal with a normal post-rest potentiation of shortening. However, the post-rest potentiation of the Ca2+ transient was absent and the cellular contractile response to isoprenaline was decreased in CK-/- mice. In vivo, echocardiographically determined cardiac function was normal at rest but the response to isoprenaline was blunted in CK-/- mice. Previously described compensatory pathways (glycolytic pathway and closer sarcoplasmic reticulum-mitochondria interactions) allow a quasi-normal SR function in isolated cells and a normal basal in vivo ventricular function, but are not sufficient to cope with a large and rapid increase in energy demand produced by beta-adrenergic stimulation. This shows the specific role of CK in excitation-contraction coupling in cardiac muscle that cannot be compensated for by other pathways.

Role of the Na(+)-Ca(2+) exchanger as an alternative trigger of CICR in mammalian cardiac myocytes

Ca(2+) influx through the L-type Ca(2+) channels is the primary pathway for triggering the Ca(2+) release from the sarcoplasmic reticulum (SR). However, several observations have shown that Ca(2+) influx via the reverse mode of the Na(+)-Ca(2+) exchanger current (I(Na-Ca)) could also trigger the Ca(2+) release. The aim of the present study was to quantitate the role of this alternative pathway of Ca(2+) influx using a mathematical model. In our model 20% of the fast sodium channels and the Na(+)-Ca(2+) exchanger molecules are located in the restricted subspace between the sarcolemma and the SR where triggering of the calcium-induced calcium release (CICR) takes place. After determining the strengths of the alternative triggers with simulated voltage-clamps in varied membrane voltages and resting [Na](i) values, we studied the CICR in simulated action potentials, where fast sodium channel current contributes [Na](i) of the subspace. In low initial [Na](i) the Ca(2+) influx via the L-type Ca(2+) channels is the major trigger for Ca(2+) release from the SR, and the Ca(2+) influx via the reverse mode of the Na(+)-Ca(2+) exchanger cannot trigger the CICR. However, depending on the initial [Na](i), the contribution of the Ca(2+) entry via the exchanger may account for 25% (at [Na](i) = 10 mM) to nearly 100% ([Na](i) = 30 mM) of the trigger Ca(2+). The shift of the main trigger from L-type calcium channels to the exchanger reduced the delay between the action potential upstroke and the intracellular calcium transient. This may contribute to the function of the myocyte in physiological situations where [Na](i) is elevated. These main results remain the same when using different estimates for the most crucial parameters in the modeling or different models for the exchanger.

Dynamic regulation of sarcoplasmic reticulum Ca(2+) content and release by luminal Ca(2+)-sensitive leak in rat ventricular myocytes

In cardiac muscle, excitation-contraction (E-C) coupling is determined by the ability of the sarcoplasmic reticulum (SR) to store and release Ca(2+). It has been hypothesized that the Ca(2+) sequestration and release mechanisms might be functionally linked to optimize the E-C coupling process. To explore the relationships between the loading status of the SR and functional state of the Ca(2+) release mechanism, we examined the effects of changes in SR Ca(2+) content on spontaneous Ca(2+) sparks in saponin-permeabilized and patch-clamped rat ventricular myocytes. SR Ca(2+) content was manipulated by pharmacologically altering the capacities of either Ca(2+) uptake or leak. Ca(2+) sparks were recorded using a confocal microscope and Fluo-3 and were quantified considering missed events. SR Ca(2+) content was assessed by application of caffeine. Exposure of permeabilized cells to anti-phospholamban antibodies elevated the SR Ca(2+) content and increased the frequency of sparks. Suppression of the SR Ca(2+) pump by thapsigargin lowered [Ca(2+)](SR) and reduced the frequency of sparks. The ryanodine receptor (RyR) blockers tetracaine and Mg(2+) transiently suppressed the frequency of sparks. Upon washout of the drugs, sparking activity transiently overshot control levels. Low doses of caffeine transiently potentiated sparking activity upon application and transiently depressed the sparks upon removal. In patch-clamped cardiac myocytes, exposure to caffeine produced only a transient increase in the probability of sparks induced by depolarization. We interpret these results in terms of a novel dynamic control scheme for SR Ca(2+) cycling. A central element of this scheme is a luminal Ca(2+) sensor that links the functional activity of RyRs to the loading state of the SR, allowing cells to auto-regulate the size and functional state of their SR Ca(2+) pool. These results are important for understanding the regulation of intracellular Ca(2+) release and contractility in cardiac muscle.