AJP - Renal Watch the video to learn how APS reaches out to developing nations.
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


Am J Physiol Renal Physiol 293: F1518-F1532, 2007. First published August 1, 2007; doi:10.1152/ajprenal.00251.2007
0363-6127/07 $8.00
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow All Versions of this Article:
293/5/F1518    most recent
00251.2007v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via ISI Web of Science (2)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Edwards, A.
Right arrow Articles by Pallone, T. L.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Edwards, A.
Right arrow Articles by Pallone, T. L.

Ouabain modulation of cellular calcium stores and signaling

Aurélie Edwards1 and Thomas L. Pallone2

1Department of Chemical and Biological Engineering, Tufts University, Medford, Massachusetts; and 2Departments of Medicine and Physiology, University of Maryland School of Medicine, Baltimore, Maryland

Submitted 30 May 2007 ; accepted in final form 27 July 2007


    ABSTRACT
 TOP
 ABSTRACT
 METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
Ouabain-like factors modulate intracellular Ca2+ concentrations and Ca2+ stores. Recently, a role for Na+-K+-ATPase Na+ transport inhibition as a pivotal event in ouabain signaling was questioned (Kaunitz JD. Am J Physiol Renal Physiol 290: F995–F996, 2006). In the present study, we used a mathematical model of Ca2+ trafficking in cytoplasm and subplasmalemmal microdomains to simulate the pathways through which ouabain can affect Ca2+ signaling: inhibition of active transport by Na+-K+-ATPase {alpha}1- and {alpha}2-isoforms, activation of inositol trisphosphate (IP3) production, and increased IP3 receptor (IP3R) conductance. A fundamental prediction is that Na+-K+-ATPase inhibition favors sarcoplasmic reticulum Ca2+ store loading, whereas Src-mediated increases in IP3 production and IP3R sensitization favor store depletion. The model predicts that {alpha}2-isoform inhibition generates a peak-and-plateau pattern of cytosolic Ca2+ concentration ([Ca2+]cyt) elevation, whereas {alpha}1-isoform inhibition yields a monophasic rise. The effects of ouabain-mediated increases in IP3 production or IP3R conductance on [Ca2+]cyt depend on their relative distributions between cellular microdomains and the bulk cytoplasm. Simulations suggest that the intracellular localization of IP3 production is a pivotal determinant of the changes in compartmental Ca2+ concentrations that can be induced by ouabain. As a consequence of sequestration of the ouabain-sensitive {alpha}2-isoform into microdomains, inhibition of the {alpha}2-isoform in rodents is not predicted to significantly affect cytosolic Na+ concentration. Model simulations support the hypothesis that ouabain can enhance agonist-evoked [Ca2+]cyt transients when its predominant effect is to inhibit {alpha}2-isoform Na+ transport and, thereby, increase Ca2+ loading into sarcoplasmic reticulum stores.

sodium-potassium-adenosine triphosphatase; sodium-calcium exchange; inositol trisphosphate; pericyte; mathematical model


OUABAIN AND OTHER CARDIOTONIC steroids bind to the extracellular NH2-terminal loop of the {alpha}-subunit of Na+-K+-ATPase at a site that is highly conserved in evolution. That binding can mediate at least two separate effects: 1) ouabain inhibits the exchange of Na+ for K+, i.e., the pump function of Na+-K+-ATPase, and 2) ouabain activates Src kinase, leading to phosphorylation of phospholipase C-{gamma}1 (PLC-{gamma}1) and inositol trisphosphate (IP3) receptor (IP3R) type 2, which are tethered to Na+-K+-ATPase in caveolae to initiate several signaling cascades (36, 39, 42). Recently, the importance of Na+ pump function has been questioned (20). Blaustein and colleagues (3, 5) hypothesized that the inhibition of Na+ export from cells raises subplasmalemmal microdomain Na+ concentration to inhibit forward-mode Na+/Ca2+ exchange (NCX), elevate microdomain Ca2+ concentration ([Ca2+]md), and increase the mass of Ca2+ within the underlying sarcoplasmic reticulum (SR) stores. SR Ca2+ loading, as well as the role of NCX as a mediator, has been repeatedly confirmed in myocytes and endothelial cells (2, 6, 18, 31). In parallel, Tian et al. (36), Xie and Cai (39), and Yuan et al. (42) elegantly and convincingly demonstrated that ouabain-mediated activation of Src kinase leads to phosphorylation of PLC-{gamma}1 to stimulate generation of IP3, accompanied by IP3R phosphorylation, thereby raising cytosolic Ca2+ concentration ([Ca2+]cyt) by release of Ca2+ from stores. The relative balance of the effects of ouabain to inhibit pump function and/or activate Src/PLC-{gamma}1 signaling may be specific to the cell type and species under study.

Recently, we provided a mathematical simulation of Ca2+ trafficking between the extracellular space, subplasmalemmal microdomains, SR stores, and cytoplasm in vascular smooth muscle cells (VSMC) (13). The model accounted for major ion transport pathways, IP3 production, IP3R and ryanodine receptor (RyR) sensitization, intracellular Ca2+ buffering, and localization of ouabain-sensitive {alpha}2-isoform Na+ pumps into putative subplasmalemmal microdomains. The model successfully predicted that [Ca2+]cyt and the magnitude of SR store loading with Ca2+ can be modulated by variation of the activity of microdomain Na+ pumps. To achieve the latter effect, it proved critical to sequester microdomains, so as to limit diffusional exchange of Na+ with the "bulk" cytoplasm (4, 6, 19). Given the importance of ouabain-like factors (OLF) in the modulation of Ca2+ signaling and hypertension, we have investigated the predictions of the model with respect to the actions of ouabain to inhibit microdomain {alpha}2-isoform Na+ pumps (6) and/or activate IP3 production after Src kinase activation (42). The simulations support the contention that isolated inhibition of Na+ pumps by ouabain can elevate [Ca2+]cyt and enhance loading of Ca2+ into the SR. Interestingly, however, the model predicts that the SR Ca2+ content will be limited, or even depleted, if IP3 production and IP3R activation are the dominant effects of the interaction of ouabain and Na+-K+-ATPase.


    METHODS
 TOP
 ABSTRACT
 METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
Our model accounts for the characteristics of the channels, transporters, and pumps that exchange ions between the extracellular space, the bulk cytosol, the microdomains, and SR stores. These compartments are denoted by the superscripts or subscripts "out," "cyt," "md," and "sr," respectively.

As described previously (see Fig. 1 in Ref. 13), the following channels and pumps are taken to be uniformly distributed over the plasma membrane: inward rectifier K+ channels, delayed rectifier K+ channels, ATP-activated K+ channels, Ca2+-activated K+ channels, voltage-activated Na+ channels, nonselective cation channels, Ca2+ pumps, and L-type voltage-dependent Ca2+ channels. In contrast, we assume that the {alpha}2-isoforms of Na+-K+-ATPase pumps (INaK,{alpha}2) and the NCX are expressed exclusively in the region of the cell membrane directly above the microdomains, whereas the {alpha}1-isoforms (INaK,{alpha}1) are restricted to the region of the cell membrane directly above the bulk cytosol (4, 6, 19). The interfaces between the SR and cytosol and between the SR and microdomains are populated by sarco(endo)plasmic reticulum Ca2+-ATPase (SERCA) pumps, RyR, and IP3R.


Figure 1
View larger version (6K):
[in this window]
[in a new window]

 
Fig. 1. Intracellular Ca2+ and Na+ concentrations ([Ca2+] and [Na+]) as a function of time (in seconds). At 300 s, the {alpha}2-isoform of Na+-K+-ATPase is entirely inhibited. Inhibition of the {alpha}2-isoform elevates microdomain [Na+] ([Na+]md), thereby reducing current through Na+/Ca2+ exchangers (NCX) and raising microdomain [Ca2+] ([Ca2+]md). Ensuing increase in Ca2+ uptake by microdomain sarco(endo)plasmic reticulum Ca2+ (SERCA) pumps results in sarcoplasmic reticulum (SR) Ca2+ loading, which increases Ca2+ release by the inositol trisphosphate (IP3) receptor (IP3R) at the SR-cytosol interface to favor elevation of cytosolic [Ca2+] ([Ca2+]cyt).

 
As described above, a pivotal assumption is that diffusional exchange between microdomains and bulk cytosol is limited. Otherwise, as shown in our previous study (13), resting cytosol-to-microdomain Na+ gradients comparable to experimental observations (38) cannot be maintained. On the basis of experimental measurements in VSMC (21), the fraction of the cell membrane directly above the cytosol is taken as 85.8% and that above the microdomains as 14.2%.

Where possible, inputs for the model were taken from recent studies of descending vasa recta pericytes (7, 8, 28) and cerebrovascular arteries (40). When parameters were not available for VSMC, they were extrapolated from data for cardiac cells (24, 33). Equations describing ionic currents across Na+-K+-ATPase pumps, NCX, store-operated channels (SOC), IP3R, and SERCA pumps are summarized below. All other model equations and parameter values can be found in the previous study (13). The predicted resting values of all 44 variables of the model are summarized in Table 1. The adopted convention is that exit of positive charge from the cell is a positive current.


View this table:
[in this window]
[in a new window]

 
Table 1. Model variables and resting values

 
Na+-K+-ATPase pump current. As a baseline, we assume that the {alpha}1-isoform of the pump is expressed only in the plasma membrane above the bulk cytosol, whereas the {alpha}2-isoform is confined to the membrane region above the microdomains (4, 6, 19). The general form of the current across Na+-K+-ATPase pumps (INaK) is that proposed by Luo and Rudy (24), which we adapt to each isoform as follows

Formula 1(1)

Formula 2(2)

Formula 3(3)

Formula 4(4)
where [Na+] and [K+] represent Na+ and K+ concentration, respectively, Vm is transmembrane potential, Km,K is K+ half-saturation constant, Km,Na,{alpha}1 and Km,Na,{alpha}2 are Na+ half-saturation constants for Na+-K+-ATPase {alpha}1- and {alpha}2-isoform currents, {psi} is ionic potential, F is Faraday's constant, R is gas constant, T is temperature, and {sigma} is electrical conductivity. The maximum currents through the {alpha}1-isoform (INaK,{alpha}1,max) and the {alpha}2-isoform (INaK,{alpha}2,max) are estimated as 15.8 and 2.37 pA, respectively (13). Implicit in the calculations is the assumption that the ratio of the number of {alpha}2- to {alpha}1-isoforms is 20:80, as reported for astrocytes (14) and mesenteric arteries (44), and that the ratio of {alpha}2- to {alpha}1-isoform turnover is 0.6 (32). The values of Km,K, Km,Na,{alpha}1, and Km,Na,{alpha}2 are taken as 1.5, 12, and 22 mM, respectively (24, 43).

NCX current. We assume that all NCX are localized above the microdomains (4, 6, 19). We assume that NCX has asymmetric affinities on the internal and external sides of the membrane and that the NCX current (INaCa) is given by

Formula 5(5)

Formula 6(6)

Formula 7(7)

Formula 8(8)

Formula 9(9)
The maximum current (INaCa,max) is taken as 200 µA/µF (i.e., 2.42 x 103 pA), the external and internal half-saturation constants for Na+ and Ca2+ (Km,Nao, Km,Nai, Km,Cao, and Km,Cai) as 87.5 mM, 12.29 mM, 1.3 mM, 3.59 µM, respectively, the constant for NCX (Kd-Act) as 0.256 µM, the NCX saturation factor at negative potential (ksat) as 0.27, and the constant for voltage dependence of NCX ({gamma}) as 0.35 (33).

SOC current. Decreases in SR Ca2+ load stimulate SOC activity. The SOC Ca2+ current in compartment j [IFormula 9(j = cyt or md)] is calculated as

Formula 10(10)
where GFormula 10 is maximum conductance of the SOC Ca2+ current and the Nernst potential (E) of ion i (valence zi) is calculated on the basis of the concentration difference between the extracellular compartment and the intracellular compartment j (j = cyt or md)

Formula 11(11)
where [i] is ion concentration. We assume a Ca2+-to-Na+ permeability ratio (PFormula 11/PFormula 11) = 8, and we use the Goldman-Huxley-Katz current equation to relate the SOC Na+ and Ca2+ currents in compartment j

Formula 12(12)
The maximum Ca2+ conductances of cytosolic and microdomain SOC are taken as 20 and 3.5 pS, respectively, and the value of KSOC is taken as 100 µM (13).

IP3R release current. Following the approach of De Young and Keiser (9), the net generation of IP3 in compartment j is calculated as

Formula 13(13)
where {nu}4 is the maximum rate of IP3 production, Ir is the rate constant for IP3 consumption, [IPFormula 13] is a constant, and {alpha}4 determines the strength of the Ca2+ feedback on IP3 production. The parameters Ir and k4 are taken as 1 s–1 and 1.1 µM (9), respectively, and we assume that {alpha}4 = 0.5. The baseline value of {nu}4 is chosen as 1.84 s–1, so that the resting value of [IP3]cyt is equal to [IPFormula 13], i.e., 240 nM (9). IIP3R in compartment j [IFormula 13(j = cyt, md)] is then calculated as

Formula 14(14)
where {nu}1 is the Ca2+ conductivity of IP3R, x010 is the fraction of receptors bound by one activating Ca2+, d1 is a ratio of kinetic rate constants (0.13 µM), and volsr is the SR volume (0.07 pl). The baseline value of {nu}1 is taken as 6 s–1 (9), and x010 is calculated as described previously (13).

SERCA pump current. We assume that the SERCA pumps are located at the cytosol-SR interface and the microdomain-SR interface and that the uptake current depends on concentrations on the internal and external sides of the SR membrane

Formula 15(15)
where IFormula 15 is SERCA pump current (j = md or cyt) and f is fraction of membrane surface above microdomain or cytosol. The maximum current (ISERCA,max) is chosen as 100 pA, a midrange value (33, 40). The forward and reverse half-saturation constants (Kmf and Kmr) are taken as 400 nM and 1.7 mM, respectively, and the Ca2+-ATPase constant (H) as 1.787 (13).

Electrodiffusive flux from microdomain to cytosol. The electrodiffusive flux of ion i (Ji,diff) is calculated as follows, such that Ji,diff yields a positive current into the cytosol when cations flow down their electrochemical gradient from microdomain to cytosol

Formula 16(16)

Formula 17(17)

Formula 18(18)
where A is the cross-sectional area at the interface between the microdomains and the cytosol (taken as 7.8 x 10–9 cm2), Di is the diffusivity of ion i, L is the distance from the center of the microdomains to the center of the cytosol (estimated as 0.5 µm), {xi}i is normalized electrical potential difference, h is a hindrance factor (0 ≤ h ≤ 1) that lumps together steric and charge-related effects and Pi,diff is a permeability (in m/s). The whole cell diffusivity of Ca2+ is taken as 0.3 x 10–5 cm2/s. The Na+-to-Ca2+ and K+-to-Ca2+ diffusivity ratios are assumed to be 1.33:0.79 and 1.96:0.79, i.e., the ratios of the diffusivities in dilute solution. The baseline value of h is chosen as 2.5 x 10–3 (13) to reflect significant sequestration of the microdomains from the bulk cytosol. Note that the electrodiffusive current of ion i is given by Ii,diff = ziF·Ji,diff.

Numerical methods. As described previously (13), the values of the 44 model variables are obtained by solving a system of ordinary differential equations. The latter is programmed with MATLAB and solved numerically on a personal computer with an Intel-based processor.


    RESULTS
 TOP
 ABSTRACT
 METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
To investigate the relative importance of hypothetical mechanisms by which ouabain might alter Ca2+ signaling, we sequentially examined its predicted effects on Na+-K+ATPase activity, IP3 production, and IP3R conductance. For convenient reference, to facilitate the interpretation of the results, the progression of the simulations is summarized in Table 2.


View this table:
[in this window]
[in a new window]

 
Table 2. Summary of simulations

 
{alpha}2-Isoform inhibition. To simulate the isolated effects of nanomolar ouabain on ouabain-sensitive Na+-K+-ATPase, we first examined the effects of a complete, isolated {alpha}2-isoform inhibition (Fig. 1). A "peak-and-plateau" pattern of [Ca2+]cyt and [Ca2+]md elevation is predicted by the model. Similar, biphasic [Ca2+]cyt elevations have been observed in endothelia and LLC-PK1 cells (31, 42). As postulated by Blaustein and colleagues (3, 6), isolated {alpha}2-isoform inhibition raises [Na+]md, thereby inhibiting Ca2+ export from the microdomains by NCX and favoring loading of Ca2+ into the SR. The result is a rise in the mass of Ca2+ within the SR, as reflected by an increase in [Ca2+]sr. The rise in [Na+]md is not accompanied by an increase in [Na+]cyt that is large enough to be experimentally observable with cytosolic fluorescent Na+ probes (i.e., Na+-binding benzofuran isophthalate), a prediction that agrees with experimental observations in vascular smooth muscle (2). The overall effect of 100% {alpha}2-isoform inhibition on SR Ca2+ stores is an increase in the mass of Ca2+ from 0.674 to 0.830 fg (Table 3).


View this table:
[in this window]
[in a new window]

 
Table 3. Mass of Ca2+ in SR stores at equilibrium

 
The mechanisms underlying the [Ca2+]cyt profile are best understood through examination of the concomitant changes in ionic currents that accompany 100% {alpha}2-isoform inhibition (Fig. 2). At baseline, the model predicts that NCX exports Ca2+ and imports Na+ (i.e., "forward-mode" operation). Since the {alpha}2-isoform is exclusively expressed in the plasmalemma above the microdomains, its inhibition elevates [Na+]md (Fig. 1), thereby reducing NCX current and the associated export of Ca2+ from the microdomains. The result is an increase in [Ca2+]md. The rise in [Ca2+]md is augmented through stimulation of Ca2+-induced Ca2+ release via IP3R at the microdomain-SR interface, partially accounting for the transient peak in Fig. 1. The ensuing increase in Ca2+ uptake by microdomain SERCA pumps results in SR Ca2+ loading, which increases Ca2+ release by IP3R at the SR-cytosol interface to favor elevation of [Ca2+]cyt. The loading of Ca2+ into the SR, however, inhibits the SOC currents, which favor [Na+]md elevation, thereby providing a brake on the rise of [Ca2+]md, [Ca2+]sr, and [Ca2+]cyt, so that they settle at plateau values.


Figure 2
View larger version (13K):
[in this window]
[in a new window]

 
Fig. 2. Ionic currents (in pA) as a function of time (in seconds). At 300 s, the {alpha}2-isoform of Na+-K+-ATPase is entirely inhibited. Solid and dashed-dotted lines denote cytosolic and microdomain currents, respectively. INaK, current through Na+-K+ATPase pumps; INaCa, current through NCX; Idiff, electrodiffusive current between microdomains and cytosol; ISOC, current through store-operated channels (SOC); ISERCA, current through SERCA pumps; IIP3R, current through IP3 receptors (IP3R). Current through the ryanodine receptor (RyR) is negligible in comparison and not shown.

 
Because of the availability of fluorescent Ca2+ probes that load into the cytoplasm, [Ca2+]cyt is the most readily obtained experimental measurement. The predicted effects of graded {alpha}2-isoform inhibition on [Ca2+]cyt are summarized in Fig. 3. The peak and plateau elevations relative to the baseline are small, even at 100% inhibition (14% and 13%, respectively). In contrast, inhibition of the {alpha}1-isoform has a much greater effect to increase [Ca2+]cyt (Fig. 3). An accounting of the mechanisms that underlie the pronounced effect of {alpha}1-isoform inhibition is provided below.


Figure 3
View larger version (9K):
[in this window]
[in a new window]

 
Fig. 3. Increases in [Ca2+]cyt, relative to its resting value, after graded inhibition of the {alpha}1- or {alpha}2-isoform.

 
{alpha}1-Isoform inhibition. Millimolar concentrations of ouabain inhibit {alpha}1- and {alpha}2-isoforms of Na+-K+-ATPase. With the use of such high concentrations of ouabain, {alpha}1-isoform Na+ pumps cannot be inhibited without simultaneous induction of {alpha}2-isoform blockade. Nonetheless, we simulated isolated {alpha}1-isoform inhibition to obtain a better understanding of its independent effects.

Inhibition of the {alpha}1-isoform raises [Na+]cyt sufficiently (cf. Fig. 4 with Fig. 1) to reverse the direction of the intercompartmental electrodiffusive Na+ flux (Eqs. 1618), which then carries Na+ from the cytosol into the microdomains. The subsequent large rise in [Na+]md reverses the operating mode of NCX, which then exports Na+ and imports Ca2+ into the microdomains. The ensuing rise in [Ca2+]md enhances Ca2+ loading into the SR, thereby stimulating RyR- and IP3R-mediated Ca2+ release into the cytosol, markedly elevating [Ca2+]cyt. Thus, despite the presumed lack of NCX above the cytosol, {alpha}1-isoform inhibition is predicted to have a marked effect to elevate [Ca2+]cyt.


Figure 4
View larger version (7K):
[in this window]
[in a new window]

 
Fig. 4. Intracellular [Ca2+] and [Na+] as a function of time (in seconds). At 300 s, the {alpha}1-isoform of Na+-K+-ATPase is entirely inhibited. Time scale is longer than in Fig. 1. Inhibition of the {alpha}1-isoform raises [Na+]cyt sufficiently to reverse the direction of the intercompartmental electrodiffusive Na+ flux, thereby significantly raising [Na+]md. Subsequent NCX-mediated rise in [Ca2+]md enhances Ca2+ loading into the SR, thereby stimulating RyR- and IP3R-mediated Ca2+ release into the cytosol and markedly elevating [Ca2+]cyt.

 
In the case of 100% inhibition of the {alpha}1-isoform, [Ca2+]md, [Ca2+]sr, and [Ca2+]cyt increase to plateau levels that are approximately twice as high as resting (preinhibition) values (Fig. 4, Table 3). The rate of increase is much more gradual than that associated with inhibition of the {alpha}2-isoform (cf. Fig. 1 with Fig. 4). The rise is slower, because the volume of the cytoplasm is larger than that of the microdomains, thus requiring more transport of Ca2+ to achieve elevation. In agreement with experimental observations (2), [Na+]cyt is also predicted to increase more than threefold. As summarized in Fig. 3, a peak phase of [Ca2+]cyt elevation does not exist (cf. Fig. 4 with Fig. 1), and the plateau value of [Ca2+]cyt increases exponentially with the degree of {alpha}1-isoform inhibition.

IP3 production. In the above-described simulations, there is only limited stimulation of IP3 production, resulting solely from [Ca2+]md and [Ca2+]cyt elevations and the resultant feedback effect of Ca2+ that enhances PLC activity. In contrast, it has been shown that ouabain can independently phosphorylate PLC-{gamma}1 downstream of Src kinase activation (30, 42). In LLC-PK1 cells, this has been shown to increase the production of IP3 and is accompanied by tyrosine phosphorylation of the type 2 isoform of IP3R, presumably enhancing its sensitivity to IP3 (42). To distinguish the consequences of these two Src-stimulated effects of ouabain (henceforth denoted "Src-related effects") from ouabain-mediated inhibition of Na+ pump isoforms, we separately examined stimulation of IP3 production and enhancement of IP3R conductance as isolated events, without the accompanying {alpha}1- or {alpha}2-isoform inhibition (Table 2).

In LLC-PK1 cells, 100 nM ouabain increased phosphatidylinositol bisphosphate hydrolysis and IP3 production by ~50–100% (see Fig. 5 in Ref. 42). We simulated this increase in phosphatidylinositol bisphosphate hydrolysis and IP3 production by increasing the production rate of IP3 through an increase in the parameter {nu}4 in Eq. 13. Figure 5 shows the effect of a 100% increase in {nu}4 starting at ts = 300 s on intracellular IP3 and ion concentrations. Data describing the time dependence of changes in IP3 concentration after ouabain stimulation are not available. We therefore examined two cases: 1) a long-lasting, 100% increase

Formula 19A(19a)
and 2) a very rapid increase in {nu}4 followed by a slow exponential decrease toward a plateau level that remains 10% higher than the preinhibition level

Formula 19B(19b)
After the step change in {nu}4, IP3 concentrations in cytosol and microdomain ([IP3]cyt and [IP3]md) rise rapidly in both cases, from 240 and 251 nM to ~500 and 570 nM, respectively. Note that [IP3]md > [IP3]cyt, since [Ca2+]md > [Ca2+]cyt and, therefore, exerts a greater stimulation of IP3 production. As IP3R-mediated Ca2+ release subsequently increases, [Ca2+]cyt and [Ca2+]md rise rapidly and [Ca2+]sr drops sharply to ~45 µM. The subsequent increase in the rate of Ca2+ uptake by SERCA pumps halts the [Ca2+]cyt and [Ca2+]md increase, accounting for the first, brief peak. The first [Ca2+]cyt and [Ca2+]md transient is followed by a more significant peak-and-plateau effect, as unloading of SR Ca2+ stores increases transmembrane SOC currents, thereby greatly increasing Ca2+ import into the bulk cytosol. In the absence of cytosolic SOC channels, this secondary peak would be negligible (dashed-dotted line on the [Ca2+]cyt profile in Fig. 5). Furthermore, the increase of the SOC current significantly raises [Na+]md and [Na+]cyt, activating Na+-K+-ATPase. A pertinent conclusion is that, for nanomolar ouabain to induce a change in [Na+]cyt, activation of IP3 production would have to accompany the inhibition of the ouabain-sensitive Na+ pump {alpha}2-isoform (cf. Fig. 1 with Fig. 5).


Figure 5
View larger version (13K):
[in this window]
[in a new window]

 
Fig. 5. Intracellular [Ca2+], [Na+], and IP3 concentration ([IP3]) as a function of time (in seconds). At 300 s, IP3 production ({nu}4) is increased everywhere by a factor of 2 (A). Increase in {nu}4 is sustained in time (see Eq. 19) (B). Step increase in {nu}4 is followed by a slow exponential decrease (see Eq. 19b). Note varying time scales. Dashed-dotted line represents [Ca2+]cyt temporal variations predicted in the absence of cytosolic SOC channels. These simulations suggest that IP3 production reduces SR Ca2+ stores.

 
In case 1, where IP3 production is sustained, [Ca2+]cyt remains elevated long after the step change (Fig. 5A). Under these conditions, the unloading of SR Ca2+ stores is predicted to be very significant, as the mass of Ca2+ in the SR falls from 0.674 to 0.247 fg at equilibrium. In case 2, where IP3 production undergoes exponential decay, intracellular concentrations return to levels that are within 10% of their resting values (Fig. 5B). The mass of Ca2+ in the SR at equilibrium is then calculated as 0.598 fg. In either case, stimulation of IP3 production acts to reduce, rather than increase, SR Ca2+ stores (Table 3). Thus a fundamental prediction of this model is that the extent to which SR stores can be loaded with Ca2+ by ouabain (2, 31) depends, in a particular cell type, on the degree to which IP3 production is concomitantly stimulated (42).

The model also predicts that the intracellular localization of IP3 production is a pivotal determinant of the changes in compartmental Ca2+ concentrations that can be induced by ouabain. Figure 6 shows a simulation in which ouabain-stimulated IP3 production (i.e., 100% increase in {nu}4) occurs only in the cytosol, and not in the microdomains. In that case, the ouabain-induced increase in [Ca2+]cyt is predicted to be much greater. The peak rise is ~260 vs. 165 nM, stimulation in the cytosol vs. both compartments. The reason is that when IP3 production increases in both compartments, the depletion of SR stores is greater, thereby limiting the availability of Ca2+ that can be released to the cytosol. Specifically, when IP3 production is stimulated in both compartments, store Ca2+ falls from 0.674 fg to 0.247 or 0.598 fg (sustained vs. transient stimulation, Eq. 19, a and b, respectively). In contrast, when IP3 production is stimulated in the cytosol alone, store Ca2+ drops to only 0.381 or 0.639 fg.


Figure 6
View larger version (8K):
[in this window]
[in a new window]

 
Fig. 6. Effect of a transient 2-fold increase in {nu}4 (i.e., IP3 production) on [Ca2+]cyt when IP3 production increases in cytosol and microdomains, in cytosol only, or in microdomains only. In isolation, microdomain IP3 production reduces [Ca2+]cyt. When combined with cytosolic IP3 production, it blunts [Ca2+]cyt elevation. Intracellular location of IP3 production is predicted to be a pivotal determinant of ouabain-stimulated changes in compartmental [Ca2+].

 
Conversely, if the twofold increase in {nu}4 occurred only in the microdomains, [Ca2+]md would rise higher, achieving a peak of 880 nM (not shown) vs. 550 nM (Fig. 5), [Ca2+]cyt would decrease (Fig. 6), and SR Ca2+ content would decrease from the baseline of 0.674 fg to 0.407 or 0.631 fg (Eq. 19, a and b, respectively).

Increase in IP3R conductance. Studies suggest that ouabain also increases the sensitivity of the IP3R through its tyrosine phosphorylation by Src kinase (42). We modeled this effect by varying {nu}1, the Ca2+ conductivity of IP3R (Eq. 14). The dependence of {nu}1 on time was modeled according to Eq. 19, a and b (i.e., a sustained elevation or an abrupt increase followed by an exponential decay). Again, simulations were performed as an isolated effect, without simultaneous inhibition of transport by Na+-K+-ATPase (Table 2).

The effects of a twofold increase in {nu}1 on intracellular Ca2+ and Na+ concentrations are similar to those of a twofold increase in {nu}4, albeit smaller in magnitude. After an increase in {nu}1, Ca2+ release via IP3R rises, leading to elevation of [Ca2+]cyt and [Ca2+]md, which results in concomitant increases in [IP3]cyt and [IP3]md via Ca2+ stimulation of IP3 production. The model predicts that the {nu}1-mediated increases in IFormula 19B and IFormula 19B immediately following the step change are significantly smaller than the equivalent {nu}4-mediated increases (3.6 and 3.7 vs. 13.1 and 11.8 pA, respectively) and, hence, the lower [Ca2+]md peak (300 nM) and the smaller [Ca2+]cyt rise (115 nM) mediated by the former. Since the depletion of SR Ca2+ stores and the subsequent SOC stimulation are less significant (Table 3), the elevations of [Na+]cyt and [Na+]md are also lower.

As in the case of increasing IP3 production (see above) the [Ca2+]cyt increase would be higher if the enhancement of IP3R conductance occurred only in the cytosol, and not in the microdomains. Conversely, [Ca2+]cyt would decrease if the change in IP3R conductance occurred only in the microdomains (results not shown).

Concentration-dependent effects of ouabain on IP3 production and IP3R conductance. Graded increases in ouabain concentration, from the physiological (pM to nM) to the pharmacological (mM) range, are often used to test its experimental effects. To examine parallel modeling predictions, we performed simulations driven by the following assumptions.

The effect of ouabain on {nu}4 and {nu}1 was estimated on the basis of the experimental data of Yuan et al. (42) in porcine LLC-PK1 cells, where the Na+-K+-ATPase {alpha}1-isoform is ouabain sensitive. In their experiments with 100 nM ouabain, they found that 1) ouabain enhanced the interaction between the {alpha}1-isoform and PLC-{gamma}1 approximately twofold (see Fig. 3 in Ref. 42), which we simulated as a rise in {nu}4 of 13, and 2) ouabain enhanced tyrosine phosphorylation of IP3R and the interaction between the {alpha}1-isoform and IP3R approximately twofold, which we simulated as a rise in {nu}1 of Eq. 14.

Given that the effects of ouabain on the interactions between the {alpha}1-isoform and PLC-{gamma}1 and between the {alpha}1-isoform and IP3R are comparable in magnitude (cf. Figs. 3C and 6D in Ref. 42), we calculated their average and assumed for simplicity that variations in {nu}4 and {nu}1 are equal. We then extrapolated the results to correlate relative increases in cytosolic {nu}4 and {nu}1 (RI{nu},cyt) to the percentage of {alpha}1-isoform inhibition (PI{alpha}1) as described in Table 4. The concentration of ouabain that inhibits 50% of the pump activity (IC50) in LLC-PK1 cells was taken as 1 µM (15). The data were best fitted to a power law, and the following regression equation was obtained

Formula 20(20)


View this table:
[in this window]
[in a new window]

 
Table 4. Relationship between percent {alpha}1-isoform inhibition and {nu}4 and {nu}1 increases in the cytosol

 
In the simulations described below, IC50 was taken as 10 µM and 1.3 nM for {alpha}1- and {alpha}2-isoforms in rodents, respectively (27, 44). The percent inhibition of the {alpha}i-isoform (i = 1, 2) as a function of ouabain concentration (Couabain) was then calculated as [1 – 1/(1 + Couabain/IC50,{alpha}i)] x 100. In addition, we assumed that, in the cytosol, variations in {nu}4 and {nu}1 depend on percent inhibition of the {alpha}1-isoform, whereas, in the microdomains, they depend on percent inhibition of the {alpha}2-isoform. Our hypotheses regarding the graded effect of ouabain on {alpha}1- and {alpha}2-isoforms and {nu}1 and {nu}4 are summarized in Table 5.


View this table:
[in this window]
[in a new window]

 
Table 5. Assumed effect of ouabain concentration on Na+-K+-ATPase inhibition, IP3 production, and IP3R conductance

 
Src-related effects of ouabain on the {alpha}2-isoform. Effects of isolated inhibition of {alpha}1- and {alpha}2-isoforms or stimulation of IP3 generation or IP3R conductance have been simulated separately (Figs. 16). We next investigated the more physiologically meaningful predictions of the model when ouabain simultaneously inhibits Na+-K+-ATPase and activates Src signaling to generate IP3 and increase IP3R conductance (Table 2). Whether the stimulation of signaling cascades via binding of ouabain to the {alpha}2-isoform is comparable to that of the {alpha}1-isoform has not been specifically tested. Therefore, the effects of 10 nM ouabain on inhibition of Na+ pump activity while simultaneously stimulating IP3 production and IP3R conductance were simulated for three variations of the ability of ouabain to induce Src signaling via the {alpha}2-isoform (Fig. 7). In case A, we assumed that the IP3-related effects of ouabain binding to the {alpha}2-isoform are identical to those of the {alpha}1-isoform. We used the same correlation to relate the percentage of {alpha}2-isoform inhibition (PI{alpha}2) to the relative increases in microdomain {nu}4 and {nu}1 (RI{nu},md), i.e., Eq. 20, RI{nu},md = 0.2052 (PIFormula 20). In case B, we assumed that the IP3-related effects of ouabain binding to the {alpha}2-isoform are 10% of those that result from binding to the {alpha}1-isoform, i.e., RI{nu},md = (0.1) x 0.2052(PIFormula 20). In case C, we assumed that the binding of ouabain to the {alpha}2-isoform has no effect on microdomain IP3 production or IP3R conductance. In each case, ouabain stimulation was assumed to evoke a transient increase in cytosolic {nu}4 and {nu}1 at ts = 300 s with an exponential decay as in Eq. 19b.


Figure 7
View larger version (8K):
[in this window]
[in a new window]

 
Fig. 7. Predicted effects of stimulation with 10 nM ouabain (at 300 s) on intracellular [Ca2+] and [Na+], where ouabain inhibits Na+-K+-ATPase, stimulates IP3 production, and enhances IP3R conductance (Table 2). In case A, IP3-related effects of ouabain binding to the {alpha}2-isoform are identical to those of the {alpha}1-isoform. In case B, IP3-related effects of ouabain binding to the {alpha}2-isoform are 10% smaller than those of the {alpha}1-isoform. In case C, ouabain stimulates cytosolic ({alpha}1-isoform-mediated) but not microdomain ({alpha}2-isoform-mediated) IP3 production and IP3R conductance. Ouabain-mediated effects on {nu}4 and {nu}1 were taken to decrease with time (according to Eq. 19b) in these simulations and the remainder of this study. These simulations predict that microdomain IP3 production (case A) decreases [Ca2+]cyt, because significant diminution of SR Ca2+ stores accompanies [Ca2+]md elevation.

 
For cases A–C, 10 nM ouabain is predicted to increase cytosolic IP3 production and IP3R conductivity by 7.7%. In the microdomains, {nu}4 and {nu}1 are calculated to increase by 138, 13.8, or 0% for cases A, B, and C, respectively. Thus, in case A, the model predicts a substantial, rapid elevation of [Ca2+]md (and [Na+]md via NCX) and significant diminution of SR Ca2+ stores. Because of significant diminution of SR Ca2+ stores, [Ca2+]cyt is predicted to decrease; it drops sharply and then rises slightly as SR Ca2+ stores partly refill, finally remaining below pre-ouabain levels. Simulations with 100 nM ouabain yield similar results (not shown). In case B, where the {alpha}2-isoform effects on microdomain {nu}4 and {nu}1 are 90% smaller, the reduction in SR Ca2+ content immediately after the step change is attenuated. As a result, [Ca2+]cyt drops initially and then rises to yield a substantial net increase. In case C, microdomain {nu}4 and {nu}1 are unaffected by ouabain, so that no depletion of SR Ca2+ stores occurs and [Ca2+]cyt increases smoothly after ouabain application.

Clearly, model predictions based on the assumptions of case A are not compatible with those experimental measurements (31, 42) that have shown a peak-and-plateau pattern of [Ca2+]cyt elevation on ouabain application. The model predicts such a pattern at 10 and 100 nM only when microdomain IP3 production and IP3R sensitization are minimal, as in cases B and C.

We then examined the effects of varying other parameters: the affinity of the {alpha}1- or {alpha}2-isoform to ouabain, the time constant characterizing the exponential decay of {nu}4 and {nu}1 following the step increase (Eq. 19b), and the distribution of NCX between the microdomains and the cytosol. Since the initial [Ca2+]cyt drop is due to the transient depletion of SR Ca2+ stores, the subsequent variations in IP3R-mediated Ca2+ release into the cytosol, and Ca2+ uptake into the SR, we also varied the rate of SERCA uptake. In all those cases, the predicted [Ca2+]cyt pattern corresponded to experimental observations only in cases B and C (results not shown). We therefore subsequently focused on simulations where the binding of ouabain to the {alpha}2-isoform has a minimal effect to stimulate microdomain IP3 production or microdomain IP3R conductance.

Dose-dependent effects of ouabain on [Ca2+]cyt and [Na+]cyt. With baseline parameters, [Ca2+]cyt is predicted to peak at 121 nM after application of 10 nM ouabain (Fig. 7, case C), whereas [Na+]cyt variations are negligible. Indeed, the {alpha}1-isoform is inhibited by only 0.1% at that concentration. As expected, stimulation with 100 nM ouabain (i.e., 1% {alpha}1-isoform inhibition) is found to elicit a higher [Ca2+]cyt peak (145 nM) and equilibrium value (117 nM); [Na+]cyt is not predicted to vary significantly.

At >10 µM ouabain, {alpha}1-isoform inhibition becomes significant, and the [Ca2+]cyt profile is concomitantly altered (Fig. 8A). As an isolated effect, {alpha}1-isoform inhibition raises [Ca2+]cyt monotonically (Fig. 4), whereas a peak-and-plateau pattern emerges when it is combined with a step increase in cytosolic {nu}1 and/or {nu}4 (IP3R conductance and IP3 production, respectively; Fig. 5). Stimulation with <1 µM ouabain does not significantly inhibit the {alpha}1-isoform, hence, the predicted [Ca2+]cyt peak-and-plateau pattern. If ouabain concentration is >10 µM (the value of IC50 for the {alpha}1-isoform in rodents), then {alpha}1-isoform inhibition is significant, and the model predicts a secondary increase in [Ca2+]cyt after the initial peak.


Figure 8
View larger version (8K):
[in this window]
[in a new window]

 
Fig. 8. Predicted effects of ouabain stimulation on [Ca2+]cyt and [Na+]cyt as a function of ouabain concentration (1 nM–1 mM). To optimize display and interpretation, timing of ouabain stimulation is staggered, so that the resultant curves do not overlap. Ouabain concentration that inhibits 50% of Na+-K+-ATPase activity (IC50) is taken as 1.3 nM for {alpha}2-isoforms and 10 µM (A and B; baseline case) or 1.3 nM (C and D) for {alpha}1-isoforms.

 
As illustrated in Fig. 8B, simulations suggest that only at >1 µM does ouabain evoke a significant [Na+]cyt increase. These predictions are similar to the observations of Arnon et al. (2), who found that low concentrations (3–100 nM) of ouabain do not affect [Na+]cyt.

In the preceding simulations (Figs. 7, 8A, and 8B), we assumed that the affinity for ouabain of the {alpha}1- and {alpha}2-isoforms of the Na+-K+-ATPase pump differs by several orders of magnitude, as observed in rodents. In many other species, including primates, the affinity difference between the two isoforms is much smaller. Therefore, we also investigated the model predictions of the effects of ouabain on [Ca2+]cyt and [Na+]cyt when the affinities of the {alpha}1- and {alpha}2-isoforms are high and equal (IC50,{alpha}1 = IC50,{alpha}2 = 1.3 nM). As shown in Fig. 8, C and D, the general shape of the [Ca2+]cyt and [Na+]cyt profiles is the same as that for rodent isoforms ({alpha}1-isoform affinity << {alpha}2-isoform affinity), but the curves are shifted with respect to ouabain concentration. Thus, if it is assumed that other cellular components are similar in primates, our results suggest that picomolar ouabain concentrations may be sufficient to elicit significant increases in [Ca2+]cyt (Fig. 8C), whereas in rodents, ouabain concentration must be in the nanomolar range to exert comparable changes (Fig. 8A). Saturation is also predicted to occur at much lower ouabain concentrations in primates: the threshold ouabain concentration beyond which no further changes in [Ca2+]cyt are seen is ~100 nM in primates and ~100 µM in rodents (Fig. 8, A and C). Similarly, our model suggests that, to significantly raise [Na+]cyt, ouabain concentration must be >1 µM in rodents (Fig. 8B) but only 100 pM in primates (Fig. 8D).

Dose-dependent effects of ouabain on SR Ca2+ stores. As described above, in isolation, Na+-K+-ATPase inhibition increases Ca2+ loading in the SR, because it leads to inhibition of NCX, favoring an increase in [Ca2+]md. Conversely, increases in IP3 production and IP3R conductivity ({nu}1 and {nu}4, Eqs. 13 and 14) tend to deplete SR Ca2+ stores by augmenting Ca2+ release through IP3R. These trends are summarized in Fig. 9. If it is assumed that the increase in {nu}4 and {nu}1 is transient, the model predicts that, at equilibrium, the mass of Ca2+ in the SR will rise after ouabain exposure (Table 3). If the increase is sustained, however, the SR Ca2+ content of rodents ({alpha}1-isoform affinity << {alpha}2-isoform affinity) will achieve a net increase only at <1 µM ouabain.


Figure 9
View larger version (12K):
[in this window]
[in a new window]

 
Fig. 9. A and B: mass of Ca2+ in SR stores at equilibrium (ordinate, in 10–18 g) after exposure to increasing concentrations of ouabain, considering the effects of ouabain on Na+-K+-ATPase inhibition only, Src-mediated increases in cytosolic IP3 production and IP3R conductance only, or both. Value of IC50 for {alpha}1-isoforms is taken as 10 µM (A; baseline case) or 1.3 nM (B). C: mass of Ca2+ in SR stores after isolated increases in IP3 production ({nu}4) and IP3R conductivity ({nu}1) in the cytosol only, the microdomains only, or both. Abscissa denotes increase in {nu}4 and {nu}1 relative to baseline. Model predicts that Na+-K+-ATPase inhibition per se increases Ca2+ loading in the SR, whereas increases in IP3 production and IP3R sensitivity per se tend to deplete SR Ca2+ stores by augmenting Ca2+ release through IP3R.

 
Effect of agonist on [Ca2+]cyt with/without prior exposure to ouabain. Experiments performed in rodents have shown that ouabain can enhance cytosolic Ca2+ transients evoked by agonists such as bradykinin (31) and serotonin (2). The ouabain-mediated enhancement is thought to result from the increase in SR Ca2+ content achieved during prior ouabain exposure. To simulate those experiments, we compared the [Ca2+]cyt transient peak that follows an increase in the cytosolic IP3 production rate (i.e., via {nu}4, resulting from agonist stimulation) with and without prior exposure to 100 nM or 1 µM ouabain. In agreement with experimental observations, the model predicts that the elevation of [Ca2+]cyt, relative to its level immediately before agonist stimulation (i.e., the increase in {nu}4), is higher if the cell was previously exposed to ouabain (Table 6).


View this table:
[in this window]
[in a new window]

 
Table 6. Effect of prior exposure to ouabain on agonist-induced [Ca2+]cyt elevation

 
The mass of Ca2+ in the SR at equilibrium does not increase monotonically with ouabain concentration (Tables 3 and 6). That relationship is complex because of the opposing effects of increases in IP3R current (which favor Ca2+ store depletion) and Na+-K+-ATPase inhibition (which favors SR Ca2+ loading). Similarly, the nonlinear relationship between stimulation of IP3 production and IP3R conductance (i.e., increases in {nu}1 and {nu}4) and percent pump inhibition contributes to overall variability.


    DISCUSSION
 TOP
 ABSTRACT
 METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
Cardiotonic steroids such as digitalis and ouabain bind to the NH2-terminal extracellular loop of the Na+-K+-ATPase {alpha}-subunit at a site that is conserved in evolution. Endogenous OLF are synthesized by the adrenal gland and hypothalamus and circulate in picomolar-to-nanomolar concentrations. Taken together, these observations support the interpretation that ouabain-Na+ pump interactions serve a fundamental physiological role (16, 41). Blaustein and colleagues (3, 5) hypothesized that ouabain-induced inhibition of Na+ export from cells elevates localized subplasmalemmal Na+ concentrations, leading to secondary reduction of Ca2+ export, or enhanced import, via the NCX. Those events were proposed to favor an increase in cellular Ca2+ concentrations and the overall mass of Ca2+ in SR stores (3, 5). In rodents, the abundant {alpha}1-isoform of Na+-K+-ATPase, which ubiquitously maintains transcellular Na+ and K+ gradients, is ouabain insensitive. In contrast, other less abundant rodent isoforms ({alpha}2{alpha}4) retain ouabain sensitivity (6, 23) and possess an NH2-terminal sorting motif that tethers them to cellular microdomains (35). In turn, those microdomains provide an interface between the plasma membrane and SR stores. Colocalization of {alpha}2-isoform Na+ pumps with endoplasmic reticulum/SR protrusions that abut the plasma membrane has been verified (4, 6, 19).

Secondary blockade of NCX is not the only means by which ouabain influences Ca2+ signaling. Tian et al. (36) and Xie and Cai (39) elegantly demonstrated that binding of ouabain to Na+-K+-ATPase stimulates tyrosine phosphorylation through activation of Src kinase. In LLC-PK1 cells, stimulation of PLC-{gamma}1, leading to IP3 generation and [Ca2+]cyt elevation, was demonstrated (42). Those events appear to be dependent on spatial relationships and protein-protein interactions within caveolae and can involve surface-expressed Na+-K+-ATPase, which resides in a "nonpumping" pool (22). Recently, on the basis of those observations, the importance of Na+ pump inhibition in ouabain signaling has been questioned (20). The predictions of this model bear directly on this issue; the IP3 generation and IP3R sensitization that occur downstream of Src and PLC-{gamma}1 activation do not readily account for the enhancement of Ca2+ store loading observed in a number of studies (2, 6, 31). That is best ascribed to secondary inhibition of NCX related to Na+-K+-ATPase inhibition (Figs. 8 and 9). Although this model provides a perspective from which to explore implications of Src- and NCX-mediated ouabain signaling, in combination and separately, it seems likely that various cell types may function in vivo by emphasizing one or the other.

A great deal of evidence exists to demonstrate control of Ca2+ signaling by ouabain. For example, ouabain enhances agonist-induced Ca2+ release from cellular stores in smooth muscle and endothelium (2, 31) and increases resting cytoplasmic Ca2+ and myogenic tone in isolated mesenteric arterioles (44). The role of ouabain-sensitive Na+ pumps and NCX in the control of microvessel contractility and salt-dependent hypertension has been extensively studied. Chronic infusion of ouabain into rodents induces hypertension (25), and a substantial fraction of humans with essential hypertension have high circulating OLF (37). Recent observations in transgenic mice have provided key evidence that supports a causal role for OLF in the induction of salt-dependent hypertension (6, 17). Reduction of expression of {alpha}2-isoform Na+ pumps, mirroring ouabain inhibition of Na+-K+-ATPase in rodents, leads to salt-dependent hypertension. Expression of a mutant, ouabain-insensitive {alpha}2-isoform eliminates the ability of ACTH infusion to induce hypertension in mice (11). SEA-0400, a selective inhibitor of NCX, eliminates ouabain-induced hypertension and selectively treats salt-sensitive, but not salt-insensitive, hypertension in rat models (14, 17). Similarly, smooth muscle overexpression of an SEA-0400-insensitive NCX1.3 isoform eliminates the antihypertensive properties of that NCX inhibitor (10, 18). These topics have been reviewed recently (5, 6, 17, 25, 37, 39).

Motivated by the importance of ouabain and OLF in modulation of myocyte contractility and generation of hypertension, we simulated the pathways through which ouabain affects Ca2+ signaling in cells: 1) inhibition of active transport by Na+-K+-ATPase {alpha}1- and {alpha}2-isoforms; 2) activation of PLC-{gamma}1, leading to IP3 production; and 3) stimulation of tyrosine phosphorylation of IP3R, leading to an increase in its conductance.

Model simplifications. We recognize that our investigation into the effects of ouabain on Ca2+ signaling, IP3 generation, and IP3R activation is limited by the paucity of detailed data and the need to combine and extrapolate information that was obtained in studies of effects in different cell types. Nonetheless, it provides an opportunity to explore the implications of signaling pathways that may exist separately or, to some degree, in combination and dispute contentions that one or the other may be predominant in a given setting (6, 20). Stated another way, it may provide some insight into which actions of ouabain best account for the relative plethora of experimental observations concerning [Ca2+]cyt changes, store Ca2+ loading, and IP3-mediated effects.

To the best of our knowledge, the effects of ouabain to enhance IP3 production and IP3R conductance have not been fully quantified with respect to their spatial and temporal changes within cells or as a function of ouabain concentration. We therefore extrapolated the data of Yuan et al. (42) concerning the ouabain dependence of the interaction of the Na+-K+-ATPase {alpha}1-isoform with PLC-{gamma}1 and IP3R to yield a correlation between increases in IP3 production ({nu}4) and IP3R conductivity ({nu}1) with {alpha}1-isoform inhibition. Furthermore, our assumption that tyrosine phosphorylation of IP3R simply raises its conductance is most likely an oversimplification. It is probable that the spatiotemporal Ca2+ release properties of the channel are altered in more complex ways, as reported for serine and threonine phosphorylation events (29). The extent to which phosphorylation leads to modulation of IP3R sensitivity to Ca2+ or IP3 is similarly uncertain. In the absence of more specific data, we chose to simply retain the kinetics of the model of De Young and Keiser (9) previously used to simulate IP3R Ca2+ currents, while varying the Ca2+ conductivity of IP3R ({nu}1).

Comparisons with experimental measurements of [Ca2+]cyt. Our simulations indicate that the net effect of ouabain on intracellular concentrations depends on the relative changes in {alpha}1- and {alpha}2-isoform transport, IP3 production (i.e., {nu}4), and IP3R conductivity (i.e., {nu}1). The model predicts that {alpha}2-isoform inhibition generates a peak-and-plateau pattern for [Ca2+]cyt (Fig. 1), whereas {alpha}1-isoform inhibition results in a monophasic increase (Fig. 4). Perhaps most importantly, the effects of increases in {nu}4 or {nu}1 on [Ca2+]cyt depend on how they are distributed between the microdomains and the bulk cytoplasm. Specifically, when we assumed that the Src-related effects of ouabain, mediated through its binding to the {alpha}2-isoform (i.e., PLC-{gamma}1 activation and the subsequent increase in IP3 production and IP3R conductance), are identical to those of the {alpha}1-isoform, the model predicted that nanomolar ouabain would induce a significant, transient drop in [Ca2+]cyt (Fig. 7). The reason for this prediction is that {alpha}2-isoforms are exclusively expressed above the microdomains, and their sensitivity to ouabain in rodents is ≥103 times higher than that of the {alpha}1-isoform, so they extensively bind ouabain, even at nanomolar concentrations. Hence, if it is assumed that variation of {nu}4 and {nu}1 is regulated by that binding, nanomolar ouabain increases {nu}4 and {nu}1 to a much greater extent in the microdomains than in the cytosol. Therefore, the model predicts a substantial elevation of [Ca2+]md, a diminution of SR Ca2+ stores, and a reduction of [Ca2+]cyt. Given that such a reduction of [Ca2+]cyt has not been observed experimentally, we surmise that the binding of ouabain to the {alpha}2-isoform may have minimal effects on IP3 production and IP3R sensitivity in microdomains. That might be explained, for example, if microdomain protein-protein interactions are such that {alpha}2-isoform binding partners do not include Src, PLC-{gamma}1, or IP3R or if the {alpha}2-isoform resides outside signaling complexes in caveolae (39). No current data address this issue, and measurement of IP3 release and its effects in microdomains, separate from those in the bulk cytosol, might be prohibitively difficult.

When cells are acutely exposed to ouabain, investigators observed a variety of temporal patterns of [Ca2+]cyt changes. Monophasic increases have been described in some cells (12, 44), whereas peak elevation followed by a stable plateau has been observed in others (31, 42). Finally, [Ca2+]cyt oscillations have been described in cell lines of renal epithelial origin (1, 26). The predictions of this model illustrate a few paradigms involving cellular transporter distributions and signaling molecule generation through which cells might display such variable patterns in experimental settings.

SR Ca2+ stores. Our results suggest that Na+-K+-ATPase inhibition per se favors SR Ca2+ loading (Fig. 9, A and B), since it induces NCX to inhibit Ca2+ export or enhance its import, depending on the balance of membrane potential and transmembrane Ca2+ and Na+ concentrations. On the other hand, increases in IP3 production and IP3R conductance are predicted to augment Ca2+ release by the receptors and, therefore, partly deplete SR Ca2+ stores (Fig. 9C).

In rats, nanomolar concentrations of ouabain enhance cytosolic Ca2+ transients evoked by agonists such as bradykinin and serotonin (2, 31). A likely explanation for this enhancement is that the SR Ca2+ load increases after ouabain application, thereby increasing receptor-mediated Ca2+ release on exposure to agonists. Our predictions support this hypothesis (Table 6), with the caveat that it depends on minimization of the extent of stimulation of IP3 production and IP3R conductance by ouabain in the cell type under study (42).

[Na+]cyt variations. Our model also provides insight into which experimental maneuvers might be expected to affect [Na+]cyt. It predicts that isolated inhibition of the ouabain-sensitive {alpha}2-isoform in rodents, as might be accomplished with nanomolar ouabain or genetic manipulations, does not significantly affect [Na+]cyt (Fig. 1). This is largely a consequence of {alpha}2-isoform sequestration above the microdomains. In contrast, inhibition of the predominant {alpha}1-isoform, above the bulk cytosol, leads to a progressive and significant [Na+]cyt elevation. Selective increases in {nu}4 and/or {nu}1, whether in the bulk cytosol and/or the microdomains, also raise [Na+]cyt, because they trigger the partial depletion of SR Ca2+ stores, thereby stimulating SOC activity and increasing Na+ import into the cell.

Arnon et al. (2) observed that low concentrations (3–100 nM) of ouabain augment Ca2+ transients in rat arterial smooth muscle without raising cytosolic Na+. This is also predicted by our model, if we assume that the binding of ouabain to the {alpha}2-isoform does not elicit significant changes in microdomain IP3 production and IP3R conductance (Fig. 8B). If, however, the latter changes were significant, [Na+]cyt would be predicted to rise by >1 mM (Fig. 7). In other words, comparison between experimental and theoretical profiles of [Ca2+]cyt and [Na+]cyt evoked by nanomolar ouabain strongly suggests that binding of the {alpha}2-isoform by ouabain has small effects on IP3R-mediated Ca2+ release in the microdomains, at least in some rodent cell types.

An interesting and testable prediction of the model is that, in species where the ouabain affinity of the {alpha}1-isoform is high and comparable to that of the {alpha}2-isoform, ouabain concentrations as low as 100 pM–1 nM might be sufficient to raise cytosolic Na+ (Fig. 8D).

Conclusion. We simulated in the present study the pathways through which ouabain affects Ca2+ signaling in cells. The model predicts that {alpha}2-isoform inhibition generates a peak-and-plateau pattern for [Ca2+]cyt, whereas {alpha}1-isoform inhibition results in a monophasic increase. The effects of ouabain-mediated increases in IP3 production or IP3R conductance on [Ca2+]cyt depend on their relative distributions between the microdomains and the bulk cytoplasm. Comparison between experimental data and predictions suggests that the binding of ouabain to the {alpha}2-isoform may have minimal effects on IP3 production and IP3R conductance in microdomains. Our model also predicts that isolated inhibition of the ouabain-sensitive {alpha}2-isoform in rodents does not significantly affect [Na+]cyt, given sequestration of the {alpha}2-isoform above the microdomains. Finally, our results suggest that Na+-K+-ATPase inhibition favors SR Ca2+ store loading, whereas increases in IP3 production and IP3R conductance favor store depletion, pointing to possible explanations for experimental observations obtained in studies of various cell types. In a particular cell type, the balance of those two effects might govern the characteristics of responses to ouabain and OLF.


    FOOTNOTES
 

Address for reprint requests and other correspondence: A. Edwards, Dept. of Chemical and Biological Engineering, Tufts Univ., 4 Colby St., Medford, MA 02155 (e-mail: aurelie.edwards{at}tufts.edu)

The costs of publication of this article were defrayed in part by the payment of page charges. The article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.


    REFERENCES
 TOP
 ABSTRACT
 METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 

  1. Aizman O, Uhlén P, Lal M, Brismar H, Aperia A. Ouabain, a steroid hormone that signals with slow calcium oscillations. Proc Natl Acad Sci USA 98: 13420–13424, 2001.[Abstract/Free Full Text]
  2. Arnon A, Hamlyn JM, Blaustein MP. Ouabain augments Ca2+ transients in arterial smooth muscle without raising cytosolic Na+. Am J Physiol Heart Circ Physiol 279: H679–H691, 2000.[Abstract/Free Full Text]
  3. Blaustein MP. Sodium ions, calcium ions, blood pressure regulation, and hypertension: a reassessment and a hypothesis. Am J Physiol Cell Physiol 232: C165–C173, 1977.[Abstract/Free Full Text]
  4. Blaustein MP, Golovina VA. Structural complexity and functional diversity of endoplasmic reticulum Ca2+ stores. Trends Neurosci 24: 602–608, 2001.[CrossRef][Web of Science][Medline]
  5. Blaustein MP, Robinson SW, Gottlieb SS, Balke CW, Hamlyn JM. Sex, digitalis, and the sodium pump. Mol Interv 3: 68–72, 2003.[Abstract/Free Full Text]
  6. Blaustein MP, Zhang J, Chen L, Hamilton BP. How does salt retention raise blood pressure? Am J Physiol Regul Integr Comp Physiol 290: R514–R523, 2006.[Abstract/Free Full Text]
  7. Cao C, Lee-Kwon W, Silldorff EP, Pallone TL. KATP channel conductance of descending vasa recta pericytes. Am J Physiol Renal Physiol 289: F1235–F1245, 2005.[Abstract/Free Full Text]
  8. Cao C, Goo JH, Lee-Kwon W, Pallone TL. Vasa recta pericytes express a strong inward rectifier K+ conductance. Am J Physiol Regul Integr Comp Physiol 290: R1601–R1607, 2006.[Abstract/Free Full Text]
  9. De Young GW, Keiser J. A single-pool inositol 1,4,5-trisphosphate-receptor-based model for agonist-stimulated oscillations in Ca2+ concentrations. Proc Natl Acad Sci USA 89: 9895–9899, 1992.[Abstract/Free Full Text]
  10. Dostanic I, Paul RJ, Lorenz JN, Theriault S, Van Huysse JW, Lingrel JB. The {alpha}2-isoform of Na-K-ATPase mediates ouabain-induced hypertension in mice and increased vascular contractility in vitro. Am J Physiol Heart Circ Physiol 288: H477–H485, 2005.[Abstract/Free Full Text]
  11. Dostanic-Larson I, Van Huysse JW, Lorenz JN, Lingrel JB. The highly conserved cardiac glycoside binding site of Na,K-ATPase plays a role in blood pressure regulation. Proc Natl Acad Sci USA 102: 15845–15850, 2005.[Abstract/Free Full Text]
  12. Echevarria-Lima J, Giestal de Araújo E, de Meis L, Rumjanek VM. Ca2+ mobilization induced by ouabain in thymocytes involves intracellular and extracellular Ca2+ pools. Hypertension 41: 1386–1392, 2003.[Abstract/Free Full Text]
  13. Edwards A, Pallone TL. Modification of cytosolic calcium signaling by subplasmalemmal microdomains. Am J Physiol Renal Physiol 292: F1827–F1845, 2007.[Abstract/Free Full Text]
  14. Golovina VA, Song H, James PF, Lingrel JB, Blaustein MP. Na+ pump {alpha}2-subunit expression modulates Ca2+ signaling. Am J Physiol Cell Physiol 284: C475–C486, 2003.[Abstract/Free Full Text]
  15. Haas M, Wang H, Tian J, Xie Z. Src-mediated inter-receptor cross-talk between the Na+/K+-ATPase and the epidermal growth factor receptor relays the signal from ouabain to mitogen-activated protein kinases. J Biol Chem 277: 18694–18702, 2002.[Abstract/Free Full Text]
  16. Hamlyn JM, Hamilton BP, Manunta P. Endogenous ouabain, sodium balance and blood pressure: a review and a hypothesis. J Hypertens 14: 151–167, 1996.[CrossRef][Web of Science][Medline]
  17. Iwamoto T, Kita S. Hypertension, Na+/Ca2+ exchanger, Na-ATPase. Kidney Int 69: 2148–2154, 2006.[CrossRef][Web of Science][Medline]
  18. Iwamoto T, Kita S, Zhang J, Blaustein MP, Arai Y, Yoshida S, Wakimoto K, Komuro I, Katsuragi T. Salt-sensitive hypertension is triggered by Ca2+ entry via Na+/Ca2+ exchanger type-1 in vascular smooth muscle. Nat Med 10: 1193–1199, 2004.[CrossRef][Web of Science][Medline]
  19. Juhaszova M, Blaustein MP. Na+ pump low and high ouabain affinity {alpha}-subunit isoforms are differently distributed in cell. Proc Natl Acad Sci USA 94: 1800–1805, 1997.[Abstract/Free Full Text]
  20. Kaunitz JD. Membrane transport proteins: not just for transport anymore. Am J Physiol Renal Physiol 290: F995–F996, 2006.[Free Full Text]
  21. Lee CH, Kuo KH, Dai J, Leo JM, Seow CY, van Breemen C. Calyculin-A disrupts subplasmalemmal junction and recurring Ca2+ waves in vascular smooth muscle. Cell Calcium 37: 9–16, 2005.[CrossRef][Web of Science][Medline]
  22. Liang M, Tian J, Liu L, Pierre S, Liu J, Shapiro J, Xie ZJ. Identification of a pool of non-pumping Na/K-ATPase. J Biol Chem 282: 10585–10593, 2007.[Abstract/Free Full Text]
  23. Lingrel J, Moseley A, Dostanic I, Cougnon M, He S, James P, Woo A, O'Connor K, Neumann J. Functional roles of the {alpha}-isoforms of the Na,K-ATPase. Ann NY Acad Sci 986: 354–359, 2003.[CrossRef][Web of Science][Medline]
  24. Luo CH, Rudy Y. A dynamic model of the cardiac ventricular action potential. I. Simulations of ionic currents and concentration changes. Circ Res 74: 1071–1096, 1994.[Abstract/Free Full Text]
  25. Manunta P, Hamilton J, Rogowski AC, Hamilton BP, Hamlyn JM. Chronic hypertension induced by ouabain but not digoxin in the rat: antihypertensive effect of digoxin and digitoxin. Hypertens Res 23 Suppl: S77–S85, 2000.
  26. Miyakawa-Naito A, Uhlén P, Lal M, Aizman O, Mikoshiba K, Brismar H, Zelenin S, Aperia A. Cell signaling microdomain with Na,K-ATPase and inositol 1,4,5-trisphosphate receptor generates calcium oscillations. J Biol Chem 278: 50355–50361, 2003.[Abstract/Free Full Text]
  27. O'Brien WJ, Lingrel JB, Wallick ET. Ouabain binding kinetics of the rat {alpha}2 and {alpha}3 isoforms of the sodium-potassium adenosine triphosphate. Arch Biochem Biophys 310: 32–39, 1994.[CrossRef][Web of Science][Medline]
  28. Pallone TL, Cao C, Zhang Z. Inhibition of K+ conductance in descending vasa pericytes by ANG II. Am J Physiol Renal Physiol 287: F1213–F1222, 2004.[Abstract/Free Full Text]
  29. Patterson RL, Boehning D, Snyder SH. Inositol 1,4,5-triphosphate receptors as signal integrators. Annu Rev Biochem 73: 437–465, 2004.[CrossRef][Web of Science][Medline]
  30. Pierre SV, Yang C, Yuan Z, Seminerio J, Mouas C, Garlid KD, Dos-Santos P, Xie Z. Ouabain triggers preconditioning through activation of the Na+,K+-ATPase signaling cascade in rat hearts. Cardiovasc Res 73: 488–496, 2007.[Abstract/Free Full Text]
  31. Pittner J, Rhinehart K, Pallone TL. Ouabain modulation of endothelial calcium signaling in descending vasa recta. Am J Physiol Renal Physiol 291: F761–F769, 2006.[Abstract/Free Full Text]
  32. Segall L, Daly SE, Blostein R. Mechanistic basis for kinetic differences between the rat {alpha}1, {alpha}2, and {alpha}3 isoforms of the Na,K-ATPase. J Biol Chem 276: 31535–31541, 2001.[Abstract/Free Full Text]
  33. Shannon TR, Wang F, Puglisi J, Weber C, Bers DM. A mathematical treatment of integrated Ca2+ dynamics within the ventricular myocyte. Biophys J 87: 3351–3371, 2004.[CrossRef][Web of Science][Medline]
  34. Shelly DA, He S, Moseley A, Weber C, Stegemeyer M, Lynch RM, Lingrel JB, Paul RJ. Na+ pump {alpha}2-isoform specifically couples to contractility in vascular smooth muscle: evidence from gene-targeted neonatal mice. Am J Physiol Cell Physiol 286: C813–C820, 2004.[Abstract/Free Full Text]
  35. Song H, Lee MY, Kinsey SP, Weber DJ, Blaustein MP. An N-terminal sequence targets and tethers Na+ pump {alpha}2 subunits to specialized plasma membrane microdomains. J Biol Chem 281: 12929–12940, 2006.[Abstract/Free Full Text]
  36. Tian J, Cai T, Yuan Z, Wang H, Liu L, Haas M, Maksimova E, Huang XY, Xie ZJ. Binding of Src to Na+/K+-ATPase forms a functional signaling complex. Mol Biol Cell 17: 317–326, 2006.[Abstract/Free Full Text]
  37. Wang JG, Staessen JA, Messaggio E, Nawrot T, Fagard R, Hamlyn JM, Bianchi G, Manunta P. Salt, endogenous ouabain and blood pressure interactions in the general population. J Hypertens 21: 1475–1481, 2003.[CrossRef][Web of Science][Medline]
  38. Wendt-Gallitelli MF, Voigt T, Isenberg G. Microheterogeneity of subsarcolemmal sodium gradients. Electron probe microanalysis in guinea-pig ventricular myocytes. J Physiol 472: 33–44, 1993.[Web of Science][Medline]
  39. Xie Z, Cai T. Na+-K+-ATPase-mediated signal transduction: from protein interaction to cellular function. Mol Interv 3: 157–168, 2003.[Abstract/Free Full Text]
  40. Yang J, Clark JW Jr, Bryan RM, Robertson C. The myogenic response in isolated rat cerebrovascular arteries: smooth muscle cell model. Med Eng Phys 25: 691–709, 2003.[CrossRef][Web of Science][Medline]
  41. Yuan C, Manunta P, Chen S, Hamlyn JM, Haddy FJ, Pamnani MB. Role of ouabain-like factors in hypertension: effects of ouabain and certain endogenous ouabain-like factors in hypertension. J Cardiovasc Pharmacol 22 Suppl 2: S10–S12, 1993.
  42. Yuan Z, Cai T, Tian J, Ivanov AV, Giovannucci DR, Xie Z. Na/K-ATPase tethers phospholipase C and IP3 receptor into a calcium-regulatory complex. Mol Biol Cell 16: 4034–4045, 2005.[Abstract/Free Full Text]
  43. Zahler R, Zhang ZT, Manor M, Boron WF. Sodium kinetics of Na,K-ATPase {alpha} isoforms in intact transfected HeLa cells. J Gen Physiol 110: 201–213, 1997.[Abstract/Free Full Text]
  44. Zhang J, Lee MY, Cavalli M, Chen L, Berra-Romani R, Balke CW, Bianchi G, Ferrari P, Hamlyn JM, Iwamoto T, Lingrel JB, Matteson DR, Wier WG, Blaustein MP. Sodium pump {alpha}2 subunits control myogenic tone and blood pressure in mice. J Physiol 569: 243–256, 2005.[Abstract/Free Full Text]



This article has been cited by other articles:


Home page
Am. J. Physiol. Renal Physiol.Home page
C. Cao, K. Payne, W. Lee-Kwon, Z. Zhang, S. W. Lim, J. Hamlyn, M. P. Blaustein, H. M. Kwon, and T. L. Pallone
Chronic ouabain treatment induces vasa recta endothelial dysfunction in the rat
Am J Physiol Renal Physiol, January 1, 2009; 296(1): F98 - F106.
[Abstract] [Full Text] [PDF]


Home page
Am. J. Physiol. Gastrointest. Liver Physiol.Home page
M. W. Musch, A. Lucioni, and E. B. Chang
Aldosterone regulation of intestinal Na absorption involves SGK-mediated changes in NHE3 and Na+ pump activity
Am J Physiol Gastrointest Liver Physiol, November 1, 2008; 295(5): G909 - G919.
[Abstract] [Full Text] [PDF]


Home page
JCBHome page
T. Cai, H. Wang, Y. Chen, L. Liu, W. T Gunning, L. E. M. Quintas, and Z.-J. Xie
Regulation of caveolin-1 membrane trafficking by the Na/K-ATPase
J. Cell Biol., September 22, 2008; 182(6): 1153 - 1169.
[Abstract] [Full Text] [PDF]


Home page
PhysiologyHome page
J. Tian and Z.-j. Xie
The Na-K-ATPase and Calcium-Signaling Microdomains
Physiology, August 1, 2008; 23(4): 205 - 211.
[Abstract] [Full Text] [PDF]


Home page
Am. J. Physiol. Renal Physiol.Home page
A. Edwards and T. L. Pallone
Mechanisms underlying angiotensin II-induced calcium oscillations
Am J Physiol Renal Physiol, August 1, 2008; 295(2): F568 - F584.
[Abstract] [Full Text] [PDF]


Home page
Am. J. Physiol. Renal Physiol.Home page
Q. Zhang, C. Cao, Z. Zhang, W. G. Wier, A. Edwards, and T. L. Pallone
Membrane current oscillations in descending vasa recta pericytes
Am J Physiol Renal Physiol, March 1, 2008; 294(3): F656 - F666.
[Abstract] [Full Text] [PDF]


Home page
Am. J. Physiol. Heart Circ. Physiol.Home page
R. M. Lynch, C. S. Weber, K. D. Nullmeyer, E. D. W. Moore, and R. J. Paul
Clearance of store-released Ca2+ by the Na+-Ca2+ exchanger is diminished in aortic smooth muscle from Na+-K+-ATPase {alpha}2-isoform gene-ablated mice
Am J Physiol Heart Circ Physiol, March 1, 2008; 294(3): H1407 - H1416.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow All Versions of this Article:
293/5/F1518    most recent
00251.2007v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via ISI Web of Science (2)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Edwards, A.
Right arrow Articles by Pallone, T. L.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Edwards, A.
Right arrow Articles by Pallone, T. L.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
Visit Other APS Journals Online
Copyright © 2007 by the American Physiological Society.