|
|
||||||||
VASCULAR BIOLOGY
Department of Biomedical Engineering, Florida International University, Miami, Florida
Submitted 22 October 2006 ; accepted in final form 20 March 2007
| ABSTRACT |
|---|
|
|
|---|
microcirculation; vascular tone regulation; calcium influx pathway(s), plasma membrane Ca2+-ATPase
Mounting experimental evidence shows that EC Vm hyperpolarization itself, independently of its impact on EC Ca2+, is an important signal in resistance vessels. This electrical signal is generated by Ca2+-activated K+(KCa) currents and most likely comprises a key component of the EDHF response, which is transmitted to adjacent SMCs via myoendothelial gap junctions to cause vessel relaxation (14, 22, 52, 53). Another candidate for EDHF is simply K+ expelled by ECs into the vascular interstitial gap affecting K+-sensitive channels and transporters in both SMCs and ECs. The EC layer is most likely the main electrotonic current conduction route (for signals generated by either SMC or EC stimulation) along the vessel wall during conducted vasomotor responses (24, 31, 32, 65, 73). In addition, EC KCa channel inhibition can greatly affect arterial vasomotion (53), suggesting an important role of EC electrophysiology in this phenomenon.
Isolated ECs can be classified into two subtypes according to their resting Vm, namely, K-type and Cl-type (56). K-type EC resting Vm levels fall between 70 and 60 mV, which is close to the Nernst potential of K+ (EK) and thus indicates a dominant K+ membrane conductance, mainly due to the inward rectifier K+ (Kir) channel. On the other hand, Cl-type EC potential at rest is usually between 40 and 10 mV, which is close to the Nernst potential for Cl (ECl) and suggests a Cl conductance dominance under resting conditions. Experiments on isolated vessels suggest that ECs under hypoosmotic stress (which activates volume-sensitive Cl conductance and brings Vm toward ECl) will not be able to hyperpolarize in response to raised extracellular K+ concentration ([K+]o), via Kir activation, and thus cannot relax the artery (25). Experimental data on isolated ECs, however, have shown that Cl-type and not K-type cells hyperpolarize in response to K+ challenge (55). To understand the complex Vm role in EC behavior and vascular tone regulation, the theoretical understanding of EC electrophysiology must be enhanced.
Previous EC mathematical modeling efforts have made important contributions toward describing both Ca2+ signaling and plasmalemmal electrical activity. Two models by Wiesner et al. (77, 78) simulated the response of human umbilical vein ECs (HUVECs) to thrombin stimulation and fluid shear stress. Wong and Klassen (80) developed a model describing both [Ca2+]i and electrical responses of vascular ECs to shear stress. Korngreen et al. (44) successfully simulated Ca2+ transients of electrically nonexcitable cells following agonist washout and intracellular Ca2+ store depletion. Schuster et al. (67) modeled the changes in Vm and [Ca2+]i following bradykinin stimulation of coronary artery ECs.
However, ECs regulate Ca2+ entry and Vm by expressing an abundant and diverse collection of plasmalemmal ion channels (55), which are for the most part absent in previous EC models. In addition, considering the balance of the other major intracellular ionic species (i.e., Na+, K+, and Cl) is essential to modeling both single cell Vm behavior and cell-to-cell electrochemical coupling as Nernst potentials are concentration gradient dependent and gap junctions are nonselective, allowing the passage of small cytoplasmic ions. Consequently, there is a need to build on previous EC modeling work to provide a more biophysically detailed model of EC plasma membrane electrophysiology and further analyze the extent to which, if any, Vm affects EC intracellular Ca2+ dynamics.
In the last four decades, mathematical modeling of biological systems has contributed to the basic understanding of physiological behavior. For some organs (i.e., the heart), multiscale models have been developed that describe function at the macroscale level while incorporating mechanisms and events at the subcellular and molecular levels. A similar theoretical progress has not been paralleled in the vasculature. This mathematical model presents a first step toward this direction. The aims of this study were 1) to deliver a mathematical model that captures experimentally observed behavior of vascular ECs (and particularly ECs from rat mesenteric arterioles) and 2) to analyze how these cell responses emerge from the nonlinear interactions of individual cellular components.
| METHODS |
|---|
|
|
|---|
|
Model mathematical expressions are given in the APPENDIX, while model standard parameters and initial conditions are shown in Tables 1 and 2, respectively. Model components and parameter values were chosen to best describe rat mesenteric artery ECs. All abbreviations and symbol definitions can be found in the text and in Tables 1 and 2.
|
|
The EC plasma membrane electrical activity is modeled according to the classical Hodgkin-Huxley model as it is applied to a variety of single cell models, including ECs (37, 47, 67, 80, 83). The EC plasma membrane is thus regarded as a capacitor with its capacitance (Cm) shunted by several ionic currents that are described in detail in this section. Kirchhoff's current law describes dynamic Vm changes by:
![]() | (1) |
Kir channel current. Kir channels are considered key for resting cell Vm regulation, but their expression is rather heterogeneous among vascular EC types and some do not express it at all (55). IKir is presently modeled following Hodgkin-Huxley formalism for the instantaneous voltage-dependent gating of the channel's conductance (Eqs. A1A3 in the APPENDIX) (50). Maximal Kir conductance is described by an increasing function of [K+]o and fits well Kir conductance data from guinea pig coronary ECs (data not shown) (75). Equations A1A3 were fitted to rat mesenteric EC Kir current-voltage (I-V) data in the 100 to +40 mV range (Fig. 2A, solid line) (17). Kir conductance and Kir velocity fitted values are in agreement with data from Ref. 75. The Kir I-V model predictions under physiological intracellular K+ concentration ([K+]i) with control and elevated [K+]o displayed the reported "crossover" effect (Fig. 2A, dashed and dotted curves) (34, 55).
|
Ca2+-activated Cl channel current.
The physiological significance of CaCC in vascular ECs remains controversial, but CaCCs can function in feedback control loops for Ca2+ entry into the cell via membrane depolarization, antagonizing the action of KCa channels (33, 62). The model ICaCC (Eqs. A6A9) was based on and fitted to rescaled HUVEC data (Fig. 3A) (27). Calcium-dependent activation of CaCC was captured by a Hill function (Eq. A6). High positive Vm enhanced Ca2+-dependent activation, producing outward rectification of ICaCC (Fig. 3A) (27, 33, 62). This Vm-dependent activation was modeled using a Boltzmann function with a large positive half-activation Vm (Eq. A8). The time constant (
) of the voltage activation was fitted with a Gaussian function (Eq. A9) (data not shown).
|
|
NSC current. NSCs permeable to Ca2+ (39, 40) constitute Ca2+ entry pathways other than store operated in ECs (55). A Ca2+-permeable NSC is present in rat mesenteric ECs and has been characterized at the single channel level. It is permeable to major intracellular cations (i.e., Ca2+, Na+, and K+), displays PKG sensitivity, and plays a key function in the EDHF response of small mesenteric arterioles to KT-5823 (a potent PKG inhibitor) (23). However, since the study (23) did not measure whole cell INSC and its gating mechanisms were not thoroughly analyzed, these data were not used in the present INSC formulation. The model INSC contributes to the balance of Na+, K+, and Ca2+ and provides an alternative pathway for Ca2+ entry into the cytosol, as observed experimentally (55). The formulation developed here to model INSC (Eqs. A16A20) is based on rescaled I-V data from a constitutively open NSC found in rabbit aorta ECs (63). This study (63) found a basally active INSC having a permeability ratio of 1:0.40:0.18 for K+:Na+:Ca2+, respectively, and this current was the major resting Vm determinant in these cells, which lack Kir channels. Like in SOC, external Ca2+also had a blocking effect on INSC Na+ permeability (PNSC,Na), whose description fits a Hill relationship similar to the one used to model PSOC,Na (Eq. A20). INSC was described by the sum of three GHK equations (Eq. A19), one for each permeant species, namely, Na+ (INSC,Na), K+ (INSC,K), and Ca2+ (INSC,Ca). GHK has been used before to model rat brain endothelial NSC (18). Additionally, INSC is not gated by store depletion or [Ca2+]i, as evidenced in calf pulmonary artery ECs (CPAECs) (55). Model INSC,Na is the major inward NSC current (data not shown), as found in EC experiments (63).
NCX current.
The NCX is an electrogenic transporter of Na+ and Ca2+ with a stoichiometry of 3:1, respectively, as found in different cell types including ECs (6, 20, 30, 77, 83). Although the presence of NCX in ECs has been confirmed, its involvement in the modulation of Ca2+ signaling is controversial (55). In CPAECs, NCX contributes to cytosolic Ca2+ removal as a low-affinity system due to the high [Ca2+]i needed for significant activation, and, in effect, this exchanger counteracts large and rapid rises in EC [Ca2+]i during early stages of stimulation (68). In cultured bovine pulmonary artery ECs, NCX was found to be voltage dependent but not responsible for resting [Ca2+]i maintenance (64). The model's INCX formulation (Eqs. A21A23) is based on previously reported equations (30, 47, 76). The INCX expression given in Ref. 47 has fewer parameters from the one given in Ref. 76, and it was adapted in this study. A Hill-type instantaneous [Ca2+]i-dependent activation term is also included in this INCX description (Eq. A21), as used in modeling NCX in cardiac myocytes and pancreatic
-cells (30, 76). INCX equations were fitted to rescaled data from mouse pancreatic
-cells (30).
NaK current.
The NaK pump has been found in various EC types under a multitude of experimental protocols, and it represents the main transmembrane ionic pump present in virtually all mammalian cell membranes (1, 70). This pump has a coupling ratio of 3Na+:2K+, and its inhibition causes a 5 to 10 mV depolarization in cultured ECs (1, 70). NaK activity maintains intracellular Na+ concentration ([Na+]i) and [K+]i low and high, respectively, thus keeping physiological Na+ and K+ gradients across the plasmalemma (70). NaK activity has been described mathematically in detail, taking into account the metabolic state of the cell (ATP levels), which is suitable for modeling studies concerned with the effects of ischemia on the functionality of cardiac myocytes (70). However, a simpler INaK formulation is applied as the effects of low ATP in EC behavior were not pursued at present. The INaK formula used here (Eq. A24) was based on previous models (47, 83) and fitted to rescaled NaK I-V data from human embryonic kidney cells expressing the rat brain NaK
1-subunit (84).
NaKCl cotransport flux. The NaKCl mechanism is electroneutral, pivotally involved in the regulation of resting cell volume and tissue permeability, and believed to constitute the major influx pathway for K+ in vascular ECs, even surpassing the NaK contribution (1, 41, 58). It also largely contributes to EC Cl transport (21). In cultured bovine corneal ECs, the NaKCl function needs extracellular Na+, K+ and Cl, and it is modulated by the cytosolic Cl concentration ([Cl]i), extracellular bicarbonate, protein kinases, and the cytoskeleton state (21). The mathematical formulation of INaKCl (Eqs. A25 and A26) is based on the work by Strieter et al. (71), which described cotransport fluxes using nonequilibrium thermodynamics. This NaKCl description requires the specification of only one parameter, namely, the cotransport coefficient (L) (71), which was adjusted to the present model (Model Parameter Estimation). Model INaKCl follows an inherent 1:1:2 stoichiometry for Na+:K+:Cl, respectively, as shown experimentally.
PMCA current.
The function of the PMCA is to extrude Ca2+ to the extracellular medium (77). In both ECs and rabbit atrial cells, the PMCA is considered a low-capacity high-affinity transport system (47, 66). In CPAECs, the PMCA function was experimentally observed to be high affinity and [Ca2+]i, time, and CaM regulated as well as linked to NCX activity via [Ca2+]i, comprising an intricate dynamic compensatory mechanism (68). The PMCA is responsible for determining resting [Ca2+]i levels and counteracting moderate perturbations in [Ca2+]i, having a higher affinity for Ca2+than NCX (68). The IPMCA formulation (Eq. A27) is based on previous descriptions (47, 83). The average Hill coefficient and [Ca2+]i for half activation of PMCA current values were calculated from cultured CPAE cell data (68).
PMCA was estimated by rescaling the value reported for a rabbit atrial cell model (47).
Fluid Compartment Model
In this EC model component, the balance of intracellular chemical species takes place (Eqs. A33A41). The model assumes that Na+, K+, and Cl can occupy the total intracellular volume while Ca2+ is restricted to discrete cellular compartmental volumes, namely, the cytosolic volume available to free Ca2+ and the IP3-sensitive store volume. Another major model assumption is that chemical species are homogeneously distributed throughout the intracellular environment (67).
Intracellular Ca2+ store. The model contains an intracellular Ca2+store mainly composed by the ER (74). It is sensitive to IP3 and considered the principal Ca2+-release site in ECs (55) via IP3R activation, releasing Ca2+ from the store into the cytosol once IP3 binds IP3R (74). Additionally, EC experimental data have suggested that ryanodine-sensitive stores coexist with IP3-sensitive ones and that the former plays an essential role in shaping EC Ca2+ signals (36, 60). However, the ryanodine-sensitive stores are not implemented in the model as they lack a quantitative description level in ECs.
IP3 generation and IP3R current.
Upon agonist binding to cell surface receptors, a G protein activates PLC, which in turn hydrolyzes phosphatidylinositol (4,5)-bisphosphate to give IP3 and diacylglycerol (not shown) (Fig. 1B). IP3 dynamics are described in the model (Eqs. A40 and A41) like in Ref. 43. Thus, agonist stimulation in the present model is simulated by a step increase in the steady-state IP3 generation rate constant, which causes an exponential increase in IP3 generation (Eq. A41) with time constant
IP3, and the slower IP3 breakdown occurs via a first-order reaction with rate constant kDIP3 (Eq. A40). This approach allows the model to be easily adapted to simulate EC responses to different types of agonists (e.g., thrombin or ACh). The IIP3R formulation selected here (Eqs. A29 and A30) is based on two receptor models found in Refs. 28 and 80. These IP3R models were combined and modified based on experimental data from porcine aortic ECs (11). The data suggest that IP3-mediated receptor activation follows a Hill function with a coefficient of 3.8, suggesting approximately four binding sites for IP3, which reflects the channel's tetrameric nature (74). The data gathered from porcine ECs also underscore an inhibitory effect of [Ca2+]i on IIP3R. Ca2+ inhibits IP3R according to a modified version of the inactivation function (Pi,IP3R) from the IP3R model in Ref. 28. The Pi,IP3R function was fitted to porcine EC data (11) (Fig. 4). Although cytosolic Ca2+ can also activate IP3R (5) and this effect is included in other IP3R models (28, 80), such activation is absent in vascular ECs (11) and thus not considered here.
|
Ca2+ buffering. Both cytosolic and ER Ca2+ buffering are included in the present model, unlike in previous EC modeling efforts. The rate of [Ca2+]i buffering by CaM (Eq. A35) was based on the HUVEC model (77). The buffering mechanism inside the IP3-sensitive store is provided by CSQN kinetics, which is described mathematically by the rapid buffering approximation (Eq. A36) as applied in Ref. 79. Mitochondrial buffering was excluded from the model since this organelle does not modulate Ca2+ signals in cultured CPAECs (42).
Sensitivity Analysis
Sensitivity analysis of the EC model was performed utilizing the Latin hypercube sampling (LHS) method and multiple regression techniques as previously described (8, 19). Parameters with significant uncertainty and/or contribution to model's behavior were chosen for the present analysis (Table 3). A variation range of 20% was assigned to parameters obtained from rat mesenteric ECs. Uncertainties of 40% and 60% were assigned to parameters obtained from EC types other than rat mesentery and cell types different from ECs, respectively. The highest uncertainty range (80%) was applied to parameters whose experimental method of determination did not allow an absolute measurement of the parameter. Although extracellular ionic concentrations can be accurately determined, uncertainty ranges of 20% were assigned to them to simulate variation in experimental media preparations in different studies. Fifty simulations were performed by selecting parameter values randomly and without replacement. Sensitivity of each output to selected model parameters is quantified by the sensitivity index (
i,j) as calculated from linear least-squares multiple regression of Eq. 2:
![]() | (2) |
)% confidence region for all the sensitivity indexes to determine whether they were statistically different from zero using the F-statistic (26). Five key simulation outputs were selected for analysis based on their impact on EC physiology: 1) resting Vm (Vm,rest), 2) resting [Ca2+]i levels ([Ca2+]i,rest), 3) maximum Vm hyperpolarization (Vm,max), 4) peak [Ca2+]i level achieved ([Ca2+]i,peak), and 5) [Ca2+]i level at the end of stimulation ([Ca2+]i,end).
|
All the parameters used in model simulations along with their reference sources are defined and shown in Table 1. A Cm of 14 pF was assumed for the model as this value falls within the range of capacitances (1025 pF) reported in the literature for cultured vascular ECs (1). A capacitance to surface area ratio of 1 µF/cm2 was assumed for the estimation of EC plasma membrane surface area. Although parameter optimization was not performed to fit integrated cellular responses, a few parameters had to be modified from their initial (literature derived) values as explained below. Adjustments in parameters such as channel conductances and permeabilities are justified when values are based on experimental data from cells other than rat mesenteric ECs, given that channel expression can significantly vary among different cell types or even in the same cell type at different cell cycle stages (56).The parameters describing [Ca2+]i-dependent activation of KCa channels were selected from Refs. 2 and 81 to render them dormant at rest, as reported in rat mesenteric ECs (16, 54). Moreover, both channel conductances were increased, aiming to achieve physiologically comparable agonist-induced rat mesenteric EC hyperpolarization (54). The CaCC conductance was decreased 10-fold to allow significant agonist-induced hyperpolarization, as normally observed in rat mesenteric arteries. The [Ca2+]IS necessary for about half-maximal SOC activation was assumed to be
50% of the total content of the model IP3 store, as reported experimentally (69). ISOC parameters (Eq. A15) were modified from their initial values to achieve a better agreement between simulations and experimentally recorded [Ca2+]i transients in rat mesenteric ECs (Figs. 2C and 6C in Ref. 54). PNSC,K,
NaK, and LNaKCl were also adjusted to achieve physiological intracellular ionic concentrations and Vm at rest.
IP3 was reduced because the ACh-mediated [Ca2+]i transient in rat mesenteric ECs reached peak
40 s earlier than the reference HUVEC responses, which were caused by thrombin (77). The maximum current via IP3R was changed due to modifications in the original IP3R mathematical formulations (28, 80).
Numerical Methods
The model contains a total of 10 nonlinear differential equations, which were coded in FORTRAN 90 and solved numerically using Gear's backward differentiation formula method for stiff differential equation systems (IMSL Numerical Library routine). The maximum time step was 4 ms, and tolerance for convergence was 0.0005.
| RESULTS |
|---|
|
|
|---|
EC model initial conditions (Table 2) were obtained by letting the system reach equilibrium using the parameters shown in Table 1 according to adjustments discussed in Model Parameter Estimation. The simulated EC had a resting Vm of approximately 20 mV, which is close to the value reported (19.2 mV) for rat mesenteric ECs in Ref. 54. Moreover, the membrane resistance, calculated by injecting a known transmembrane current (Istim = ±1 pA) and observing the magnitude of Vm change, was close to 2 G
, which is in the range (15 G
) reported for ECs (56). Blockade of NaK under resting conditions caused a
4 mV depolarization in Vm, which is close to the range (510 mV) recorded from ECs after [K+]o-free solution or ouabain-induced NaK inhibition protocols (1). Resting model [Ca2+]IS was physiological (
2.25 mM) (74).
Resting Vm Modulation
Figure 5 presents model responses for different inhibition levels of VRAC. Blockade of VRAC was simulated by multiplying the GVRAC parameter by 0.0, 0.25, and 0.50 for 100%, 75%, and 50% levels of inhibition, respectively. Inhibition of VRAC was performed for 120 s starting at three different time points, namely, time = 100, 320, and 680 s to mimic the experimental protocol in Fig. 11D from Ref. 55. As Fig. 5 shows, the model was able to replicate the experimentally observed EC Vm hyperpolarization response when a compound such as NPPB or mibefradil blocked GVRAC (Fig. 11D in Ref. 55). Model outputs rapidly reached a stable Vm condition, and the magnitude of the hyperpolarization for 100% inhibition of VRAC was similar to the experiments (Fig. 11D in Ref. 55).
|
100 s) were qualitatively similar to experimental results in cultured ECs (Fig. 8C from Ref. 55). Vm hyperpolarization (solid line) was not as large as recorded in this particular experiment, and the new poststimulation resting Vm value was established more rapidly than in the experimental tracing. Nonetheless, Vm hyperpolarization magnitudes achieved in the model (Fig. 6A, solid and dotted lines) were within the range reported for nine ECs under similar [K+]o challenge conditions (Fig. 8D in Ref. 55). Doubling the Kir conductance (dotted line) helped achieve a larger Vm hyperpolarization level. However, when [K+]o was restored to the control value, a short negative transient in Vm was observed (due to faster EK increase than inactivation of IKir), and resting Vm became more negative than its original level. Blockade of VRAC drove resting Vm to hyperpolarized potentials (Fig. 6A, dashed-dotted line) compared with control. Raising [K+]o under VRAC blockade caused EC Vm depolarization instead of the hyperpolarization seen previously. Depolarization upon [K+]o stimulation also occurred when Kir was blocked (dashed line), but its magnitude was lower than for VRAC inhibition.
|
|
8 mM. Blockade of Kir channels (dashed line) completely abolished the hyperpolarizing response, and only depolarization was obtained by increasing [K+]o. Inhibition of VRAC (dashed-dotted line) greatly increased the maximum Vm hyperpolarization (64 mV) achieved by raising [K+]o to
3 mM, which was lower than the control case (8 mM). These model results may relate to the relaxation profiles for [K+]o stimulation obtained from rat mesenteric artery data (Figs. 2 and 3A in Refs. 45 and 22, respectively). The optimum [K+]o value for maximum EC Vm hyperpolarization under control conditions (8 mM; Fig. 6B) was close to the range of 10 mM (Fig. 2 in Ref. 45) to 14 mM [K+]o (Fig. 3A in Ref. 22) necessary for maximum rat mesenteric artery relaxation for the control case. Figure 6B also shows that EC Vm responses to [K+]o changes, when Kir is functional, depend on how K+ levels are varied and on GVRAC. For instance, changing [K+]o from 0.5 to 5 mM depolarized Vm under control conditions but hyperpolarized it under VRAC blockade. In contrast, adjusting [K+]o from 5 to 10 mM caused Vm hyperpolarization for control but depolarization under VRAC block. Model Responses to Agonist Stimulation
Figure 7 depicts the overall dynamics of the main components of the EC model, namely, Vm, ionic currents, and intracellular species concentrations, upon short-term agonist stimulation. The agonist-induced EC response was simulated by changing the steady-state IP3 generation rate constant from 0 to 5.5 x 108 mM/ms at time = 20 s and maintaining for 100 s. Thus, as IP3 slowly approached a constant value (Fig. 7D), [Ca2+]i levels increased (Fig. 7D) prior to the activation of KCa channels (Fig. 7H) and before Vm began to hyperpolarize (Fig. 7A). As [Ca2+]i approached its peak, INCX and IPMCA (Fig. 7C) and ISERCA (Fig. 7F) achieved maximum activity, which decayed as cytosolic Ca2+ decreased, thus reflecting their Ca2+ dependence. The depolarizing ISOC (Fig. 7B) was enhanced during the response as Vm hyperpolarized and to a small extent by IP3-sensitive store depletion, with ISOC,Ca magnitude being larger than ISOC,Na at rest and during stimulation. INSCs (Fig. 7E) also exerted a depolarizing effect on Vm as INSC,K decreased, and both INSC,Ca and INSC,Na became more negative, following changes in electrochemical gradients. The largest hyperpolarizing currents came from the activity of KCa and Kir channels (Fig. 7H). ICaCC and IVRAC (Fig. 7I), on the other hand, were mostly responsible for Vm depolarization (more so than NSC and SOC). Given the activity of ionic channels, [K+]i and [Cl]i levels decreased while [Na+]i increased during agonist stimulation (Fig. 7G). NaK activity (Fig. 7C) was minimal as Vm reached maximum hyperpolarization, which hampered normal NaK function by hindering Na+ efflux, but increased following the rise in [Na+]i and Vm repolarization. The major intracellular Ca2+ current came from IP3R release of Ca2+ (Fig. 7F), which peaked around the same time as [Ca2+]i. The Cl component of the cotransporter (INaKCl,Cl) slightly increased during stimulation (Fig. 7I). Moreover, as [Ca2+]i reached peak, [Ca2+]IS continued to decrease, suggesting a dynamic uncoupling of intracellular Ca2+ in this phase of the transient (Fig. 7, D and G).
Figure 8, A and C, illustrates long-term model Ca2+ and Vm dynamics, respectively, under agonist stimulation (at time = 20 s) as performed in Fig. 7. Figure 8 also shows how simulated transmembrane Ca2+ influx (via NSC and SOC; Fig. 8B) and its driving force (Fig. 8D) change during long-term stimulation. The model predicted that [Ca2+]i changes are followed by variations in Vm (Fig. 8, A and C), and both responses displayed a similar transient pattern, as reported in experiments (55, 72). These simulations replicate characteristic Ca2+ and Vm responses experimentally recorded from rat mesenteric ECs under ACh stimulation (Figs. 2D and 6C in Ref. 54) at least for the first 100 s of stimulation reported in the study. Blockade of Vm hyperpolarization, by either direct KCa inhibition (i.e., SKCa conductance and IKCa conductance set to 0) or raised [K+]o from 5 to 35 mM, caused a small reduction in [Ca2+]i during the transient but did not affect the shape and dynamics of the agonist-induced [Ca2+]i signal (Fig. 8A), as found experimentally (Fig. 6C in Ref. 54). Once KCa channels were inhibited, model Vm was unable to hyperpolarize, remaining essentially constant (Fig. 8C), and thus [Ca2+]i alterations could no longer affect Vm significantly (Fig. 8, A and C). High [K+]o was not as effective in reducing the [Ca2+]i transient as direct KCa blockade, which agrees with experimental data (Fig. 6C in Ref. 54). Vm hyperpolarization increased [Ca2+]i. The relative effect of Vm, however, was not homogeneous throughout the transient, becoming more significant as the later part of the transient approached (i.e., [Ca2+]i increased
24% and
64% at peak and plateau, respectively, from hyperpolarization block to control conditions). Control transmembrane Ca2+ currents comprised
45% of the total Ca2+ flux (mainly NSC, SOC, and net store release) at 100 s after the initiation of the stimulation but comprised
87% of the total at 1,000 s. Agonist-stimulated Ca2+ influx currents increased more than twofold from hyperpolarization block to control conditions (Fig. 8B), although the maximal Ca2+ entry driving force change was only
35 mV (Fig. 8D), which was just 20% of the total resting driving force. Activation of SOC by store depletion became significant in the latter part of the Ca2+ transient (Fig. 8B).
To better examine the impact of Vm on [Ca2+]i, voltage-clamp simulations were performed under resting and agonist stimulation conditions. Figure 9 depicts [Ca2+]i as a function of clamped Vm (i.e., the Vm was set to a given value and [Ca2+]i was estimated at rest and following stimulation). The resting, peak, and plateau [Ca2+]i in Fig. 9 were normalized with respect to their corresponding values after being clamped at 20 mV. (Note that 20 mV is the Vm under control conditions or during agonist stimulation with KCa channels blocked.) Plateau and resting Ca2+ levels were more sensitive to Vm (and consequently to transmembrane Ca2+ influx) than peak Ca2+ levels, with resting levels a little more sensitive than plateau levels.
|
The outputs of all 50 simulations as well as the model control response are shown in Fig. 10. Only one of the simulations, indicated by the arrow, exhibited extraneous behavior, and it was discarded from the further analysis as an outlier. The remaining 49 simulations were used for the application of the LHS method. Table 4 presents semirelative sensitivity coefficients calculated according to Eq. 2 for the five selected outputs. Vm,rest and Vm,max were significantly sensitive to 4 and 6 parameters, respectively, whereas [Ca2+]i,rest, [Ca2+]i,peak, and [Ca2+]i,end were significantly affected by 2, 10, and 6 of the parameters, respectively.
PMCA was the parameter that had the most significant impact on model predictions, and [Ca2+]i,peak was the output with the largest sensitivity on model parameters.
|
|
| DISCUSSION |
|---|
|
|
|---|
Model Verification
The model's performance was assessed by its ability to reproduce established features of rat mesenteric EC behavior while utilizing physiological parameter values. The model's resting values for Vm, for the four intracellular ionic concentrations, and for store Ca2+ agree with experimental data. Whole cell conductance was also within the physiological range, and this provides an overall validation for the descriptions of transmembrane currents. Vm responses to NaK blockade, to VRAC blockade (Fig. 5), extracellular K+ challenge (Fig. 6), and agonist stimulation (Fig. 7A) also validate plasma membrane electrophysiology. Intracellular Ca2+-handling machinery was tested against experimental data with agonist-induced stimulation (Fig. 8A). Most importantly, proper model behavior in these scenarios was accomplished by adjusting (within physiological ranges) a limited number of unknown parameters (i.e., SKCa conductance, IKCa conductance, CaCC conductance, KSOC,CaIS,
NaK, PNSC,K, LNaKCl,
IP3, and
IP3R).
Resting Vm Modulation
Experimental data from bovine pulmonary artery ECs (Fig. 11D in Ref. 55) have shown that VRAC is crucial for establishing resting EC Vm. The three levels of simulated VRAC blockade (Fig. 5) can be interpreted as different levels or types of inhibition (i.e., distinct blocker concentrations or degrees of hyperosmotic stress) or dissimilar levels of VRAC expression in ECs. Model simulations concur with VRAC inhibition experiments (25, 55). The model also supports the notion that different levels of VRAC and Kir channel expression or activation can lead to EC switching from K type (resting Vm close to EK) to Cl type (resting Vm near ECl), or vice versa, as proposed to occur throughout the cell cycle (56) or using channel inhibition protocols (Fig. 7b in Ref. 25). New experimental studies are required to verify model predictions (Fig. 5) and confirm basal VRAC activity in rat mesenteric ECs as found in other ECs (55, 72).
The EC model was used to investigate two related hypotheses: one regarding the ability of ECs to hyperpolarize in response to raised [K+]o under certain conditions (25) and another suggesting that K+ is the EDHF (22, 45) (Fig. 6). The former hypothesis was formulated based on rat mesenteric artery experiments showing that these vessels could relax to increased [K+]o (from 5.88 to 10.58 mM) only if the VRAC in ECs was blocked. To explain these results, it was postulated that ECs in intact vessels can only hyperpolarize to raised [K+]o if they are switched from Cl- to K-type cells by VRAC inhibition, using either channel antagonists or hyperosmotic media. However, isolated EC experiments (Fig. 8, C and D, in Ref. 55) have demonstrated that when EC Vm is close to ECl (Cl type), the cell responds to raised [K+]o (from 6 to 12 mM) by hyperpolarizing. This Vm behavior was captured by the model (Fig. 6A) and is attributed to the crossover effect of the Kir I-V relationship. Figure 6A also shows that a K-type EC (i.e., after VRAC block) depolarizes upon [K+]o challenge (from 5 to 12 mM). The model Vm (Fig. 6A, solid and dotted lines) achieved a new resting value (during [K+]o stimulation) faster than what was observed experimentally (Fig. 8C in Ref. 55). This may attributed to a slower increase of [K+]o in the experiments relative to the simulations. Despite these minor differences, there is an overall agreement between simulated EC Vm behavior and experiments in isolated cells.
The model suggests that hyperpolarization and depolarization can occur in response to K+ increase, under both control and VRAC block conditions, depending on the initial and final K+ concentration (Fig. 6B). Although the model predicts depolarization after VRAC blockade for the particular K+ concentration range examined in Ref. 25, simulations do not directly translate to EC behavior in intact vessels as coupling to SMC may alter EC electrophysiology. Interestingly, Fig. 6B shows that a large Vm hyperpolarization is possible also under VRAC block, when the initial value of [K+]o is low (36 mV, from 1 to
3 mM [K+]o). Under control or Kir block conditions (i.e., solid and dotted lines in Fig. 6B), resting Vm at low starting [K+]o (e.g., 0.5 mM) is rather hyperpolarized (i.e., about 40 mV). This is attributed to a very negative EK. As the [K+]o increases, the increase in EK drives Vm to more depolarized values. As [K+]o increases between 5 and 8 mM, Vm hyperpolarizes, when Kir are active, due to the crossover effect. When VRAC is blocked (i.e., dashed-dotted line in Fig. 6B) and at very low [K+]o, the model predicts a significant depolarization as a result of a decrease in INaKCl that leads to the establishment of a larger transmembrane Na+ gradient and an increase in INSC,Na.
Assuming that a similar behavior as predicted in Fig. 6B occurs in intact vessels, the model suggests that by inhibiting VRAC, the EC Vm hyperpolarization response to [K+]o gets larger (for a particular [K+]o range), which may explain why rat mesenteric artery relaxations are only observed after this channel is blocked (25). Additionally, experimental setup differences could create different levels of VRAC activation and thus dramatically affect the predicted profiles (Fig. 6B). This may explain in part the variability observed in the response of rat small mesenteric arteries to [K+]o (25). This issue merits further experimental and theoretical (via coupled EC-SMC models) investigations. Larger hyperpolarizing responses in the model might be also possible with higher values for resting [K+]i.
Regarding the K+ as EDHF hypothesis, model simulations (Fig. 6) support the idea that Kir channels on ECs may play an important role in the relaxation of rat mesenteric arteries upon [K+]o stimulation and thus in part of the EDHF response (23, 45). As SMCs and ECs in rat mesenteric arteries are electrically coupled via myoendothelial gap junctions, EC Vm hyperpolarization caused by [K+]o-induced Kir activation is then able to be transmitted to the adjacent SMCs and thus relax the vessel (23, 25, 45, 65). However, simulations are based on an isolated EC model and may not correspond to the in vivo EC-SMC coupling situation, which might shift model response profiles (Fig. 6B) and may explain why maximum artery relaxations (23, 45) are observed at [K+]o levels higher than predicted for maximum EC Vm hyperpolarization. Additionally, elevated [K+]o may activate NaK in SMCs, causing further SMC Vm hyperpolarization and thus affect relaxation profiles (22). Therefore, theoretical studies using EC-SMC integrated models are necessary to further analyze the role of EC Kir channels and K+ in the EDHF response.
Model Responses to Agonist Stimulation
Following the examination of EC plasma membrane electrophysiological behavior, the next task was to investigate if the model was able to replicate EC responses to agonist stimulation. The overall dynamic agonist-induced responses of the main model components are displayed in Fig. 7. All EC model components performed their expected roles in the agonist-induced response, in keeping with literature (1, 55, 61). Model behavior can be explained in the following way. IP3 activates IP3R, releasing store Ca2+ into the cytosol, which in turn activates KCa channels and thus causes Vm hyperpolarization, increasing the driving force for Ca2+ entry. A significant increase in IKir during hyperpolarization (a result of the negative slope of the Kir I-V relationship; Fig. 2A, dashed line) acts as a positive feedback to hyperpolarization. NCX, PMCA, and SERCA work to counteract this sudden elevation in [Ca2+]i, shaping the Ca2+ transient peak. As [Ca2+]i reaches a peak, [Ca2+]IS continues to decrease during agonist stimulation, and such intracellular Ca2+ dynamic uncoupling has been observed in CHO cells stimulated with the IP3-generating agonist ATP (7). As IP3-sensitive stores are emptied and Vm hyperpolarizes, ISOCs are activated and contribute to the sustained [Ca2+]i rise. INSCs, on the other hand, tend to offset the hyperpolarizing effect of KCa activation by passively following their electrochemical gradients. The overall depolarizing effect on Vm exerted by NSC during agonist stimulation, as INSC,K decreases and both INSC,Ca and INSC,Na become more negative (Fig. 7E), has been observed experimentally (56). Moreover, the [Ca2+]i increase and Vm hyperpolarization activate CaCCs, while IVRACs increase by Vm hyperpolarization alone (Fig. 7I). These two Cl conductances are the main players in the Vm repolarization process, which decreases the Ca2+ influx driving force. As a result of the agonist-induced processes described above, the levels of [K+]i, [Cl]i, and [Na+]i are altered (Fig. 7G). During stimulation, NCX, NaK, and NaKCl activities tend to compensate intracellular ionic species changes and thus restore cellular homeostasis, thus reflecting the importance of these model components.
Figure 8 illustrates that the model reproduces dynamic [Ca2+]i and Vm signals recorded from ACh-stimulated rat mesenteric ECs under both control and KCa blockade (54)