From the * Departments of Medicine, Surgery, and Physiology, University of Connecticut, Farmington, Connecticut 06030; Department of Physiology and Biophysics, University of Washington, Seattle, Washington 98195; and § Department of Pharmacology and Cellular & Molecular Physiology, Yale University, New Haven, Connecticut 06520
![]() |
ABSTRACT |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Cytosolic calcium acts as both a coagonist and an inhibitor of the type 1 inositol 1,4,5-trisphosphate (InsP3)-gated Ca channel, resulting in a bell-shaped Ca dependence of channel activity (Bezprozvanny, I., J. Watras, and B.E. Ehrlich. 1991. Nature. 351:751-754; Finch, E.A., T.J. Turner, and S.M. Goldin. 1991. Science. 252: 443-446; Iino, M. 1990. J. Gen. Physiol. 95:1103-1122). The ability of Ca to inhibit channel activity, however, varies dramatically depending on InsP3 concentration (Combettes, L., Z. Hannaert-Merah, J.F. Coquil, C. Rousseau, M. Claret, S. Swillens, and P. Champeil. 1994. J. Biol. Chem. 269:17561-17571; Kaftan, E.J., B.E. Ehrlich, and J. Watras. 1997. J. Gen. Physiol. 110:529-538). In the present report, we have extended the characterization of the effect of cytosolic Ca on both InsP3 binding and InsP3-gated channel kinetics, and incorporated these data into a mathematical model capable of simulating channel kinetics. We found that cytosolic Ca increased the Kd of InsP3 binding ~3.5-fold, but did not influence the maximal number of binding sites. The ability of Ca to decrease InsP3 binding is consistent with the rightward shift in the bell-shaped Ca dependence of InsP3-gated Ca channel activity. High InsP3 concentrations are able to overcome the Ca-dependent inhibition of channel activity, apparently due to a low affinity InsP3 binding site (Kaftan, E.J., B.E. Ehrlich, and J. Watras. 1997. J. Gen. Physiol. 110:529-538). Constants from binding analyses and channel activity determinations were used to develop a mathematical model that fits the complex Ca-dependent regulation of the type 1 InsP3-gated Ca channel. This model accurately simulated both steady state data (channel open probability and InsP3 binding) and kinetic data (channel activity and open time distributions), and yielded testable predictions with regard to the regulation of this intracellular Ca channel. Information gained from these analyses, and our current molecular model of this Ca channel, will be important for understanding the basis and regulation of intracellular Ca waves and oscillations in intact cells.
Key words: cerebellum; inositol 1,4,5-trisphosphate receptor; intracellular calcium channel; calcium signaling; kinetic model ![]() |
INTRODUCTION |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
The inositol 1,4,5-trisphosphate (InsP3)1-gated Ca
channel plays a critical role in the regulation of intracellular Ca concentrations after hormonal activation of
the phosphoinositide cascade in a wide variety of cell
types (Abdel-Latif, 1986; Berridge and Irvine, 1989
;
Clapham, 1995
). The mechanisms underlying regulation of the InsP3-gated channel are poorly understood;
however, the involvement of both cytosolic Ca and
InsP3 in the regulation of this channel is well established (Iino, 1990
; Payne et al., 1990
; Bezprozvanny et
al., 1991
; Finch et al., 1991
; Bootman et al., 1995
; Kaftan et al., 1997
). That is, InsP3 is absolutely required for
the opening of the channel, though channel activity
can be modified by cytosolic Ca. Elevating the cytosolic
Ca concentration from 0.01 to 0.3 µM in the presence
of 2 µM InsP3, for example, resulted in a dramatic increase in channel activity (Bezprozvanny et al., 1991
).
Further elevation of cytosolic Ca, however, decreased
channel activity, with near complete inhibition of the
channel in the presence of 5 µM Ca (Bezprozvanny
et al., 1991
). Three basic models have been proposed
to explain this bell-shaped Ca dependence of the channel differing in terms of the effects of InsP3 on the bell-shaped Ca dependence. In one model (Bezprozvanny,
1994
), the peak of the bell-shaped Ca dependence is
unaffected by InsP3 concentration. The other two models predict either a shift to the left (Tang et al., 1996
) or
the right (De Young and Keizer, 1992
) in the peak of
the bell-shaped Ca dependence of channel activity as
InsP3 concentrations increase.
Recent studies indicate that the InsP3 concentration
influences the Ca dependence of InsP3-induced Ca release, with an increased sensitivity to Ca-dependent inhibition occurring at low InsP3 concentrations (Hannaert-Merah et al., 1995; Kaftan et al., 1997
). Thus, the
peak of the bell-shaped Ca dependence of InsP3-gated channel activity shifted to the right as InsP3 concentration was increased from 0.2 to 2 µM (Kaftan et al.,
1997
). Further elevation of InsP3 concentration to 180 µM overcame the Ca-dependent inhibition of channel
activity (Kaftan et al., 1997
), which may reflect the presence of a second low affinity InsP3 binding site on the
channel complex. Recently, we modeled this complex
steady state behavior of the InsP3-gated channel by extending a basic 8-state model proposed previously (De
Young and Keizer, 1992
) to 16-states, assuming that
there are two InsP3 (with dissociation constants of 0.3 and 10 µM) and two Ca binding sites per InsP3 receptor in the tetrameric channel complex (Kaftan et al.,
1997
). In the present report, we explored in further detail the Ca dependence of InsP3 binding, as well as the
influence of Ca and InsP3 on the mean open times of
the InsP3-gated channel, and incorporated these data
into a novel mathematical model of the channel. The proposed model differs conceptually from previous
models in that it focuses on the tetrameric channel
complex as a unit, but more importantly it provides a
better fit of the steady state channel and binding data,
and simulates the kinetics of channel activation.
![]() |
METHODS |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Isolation of Microsomes
Microsomes were isolated from canine cerebellum, as described
previously (Watras et al., 1991). To minimize proteolysis, all solutions used for the isolation/storage of the cerebellar microsomes included a protease inhibitor cocktail (5 µg/ml leupeptin, 2 µg/ml aprotinin, 1 µg/ml pepstatin A, 0.1 mM PMSF). In brief, cerebellum was homogenized in 6 vol of buffer A (20 mM HEPES, 0.1 mM
EGTA, 5 mM NaN3, pH 7.2, 4°C) using a Brinkman polytron. The
suspension was centrifuged at 4,000 g for 20 min, and then the
supernatant was centrifuged at 100,000 g for 30 min. The pellet
from the latter centrifugation was resuspended in 3 vol of buffer
B (20 mM HEPES, 0.6 mM KCl, 20 mM Na4P2O5, 5 mM NaN3; pH
7.2, 4°C), and then centrifuged at 4,000 g for 20 min. The resulting supernatant was centrifuged at 100,000 g for 30 min. The final pellet was resuspended in buffer C (10 mM MOPS, 10% sucrose, pH 7.0) to a protein concentration of ~6 mg/ml (as determined using a protein assay kit; Bio-Rad Laboratories; Bradford,
1976
), quick frozen in liquid nitrogen, and stored at
80°C.
InsP3 Binding
InsP3 binding was determined by a centrifugation technique (Benevolensky et al., 1994), using buffer identical to that used for
the InsP3-gated channel activity measurements (Bezprozvanny et al., 1991
). Binding media contained 250 mM HEPES, 110 mM
Tris, 1 mM ATP, 0.5 mM EGTA, pH 7.35, 0.6-1,200 nM 3H-InsP3,
and CaCl2 (to yield the specified free Ca concentrations). InsP3
binding was initiated by addition of cerebellar membranes to the
above buffer. After an 8-min incubation, InsP3 binding was terminated by centrifugation (40,000 g for 10 min). The supernatant was removed by aspiration. The pellets were dissolved in Soluene 350, and then assayed for radioactivity by liquid scintillation counting. Nonspecific binding was determined in the presence
of 40-160 µM nonradioactive InsP3. Scatchard analysis of InsP3
binding was performed using an iterative curve-fitting routine
for two binding sites (SigmaPlot; Jandel Scientific). All binding
studies were replicated three to five times, with different microsomal preparations.
Calculations of total Ca needed to obtain a specified free Ca
concentration were based on previously published association constants (Fabiato, 1988), adjusted to pH 7.35 (at the specified temperature; 0° or 22°C). The apparent association constants
(M
1) for Ca-EGTA and Ca-ATP used for these calculations were
1.055 × 107 and 7.595 × 103, respectively, at 0°C, and 1.244 × 107
and 6.312 × 103, respectively, at 22°C.
Measurement of Channel Activity
InsP3-gated channel activity was measured as described previously
(Bezprozvanny et al., 1991; Kaftan et al., 1997
), using 53 mM Ba in
the trans chamber as current carrier. The planar lipid bilayer
contained phosphatidylethanolamine and phosphatidylserine (3:1
wt:wt; Avanti Polar Lipids), dissolved in decane (20 mg/ml). Cytosolic bilayer solutions contained 500 µM Na-ATP, 500 µM
EGTA, 110 mM Tris, and 250 mM HEPES, pH 7.35. Luminal solutions contained 53 mM Ba(OH)2, 250 mM HEPES, pH 7.35. Calibrated CaCl2 was added to the cytoplasmic solution to obtain
the desired free Ca concentration (Fabiato, 1988
). The free Ca
concentration was checked by spectrofluorometric measurements using BTC (Molecular Probes, Inc.). InsP3-gated channel activity was initiated by addition of InsP3 (at the specified concentration). The number of channels in each experiment was estimated from the maximum number of channels observed simultaneously in the bilayer (Horn, 1991
). The InsP3 dependence of
the open probability, measured at a fixed Ca concentration, was
used to correct for variations in the maximum open probability
among individual channels. Transmembrane voltage was maintained at 0 mV and the single channel current was amplified
(Warner Instruments) and stored on VHS tape (Instrutech
Corp.). Data were filtered at 1 kHz and digitized at 5 kHz for
computer analysis using pClamp 6.0 (Axon Instruments).
Modeling
Modeling InsP3-gated channel activity was done by considering
the tetrameric channel complex as one functional unit having four medium affinity InsP3 binding sites (Kd = 0.3 µM at 22°C), four low affinity InsP3 binding sites (Kd = 10 µM at 22°C), and four Ca binding sites. The mathematical model assumes mass action kinetics with first order on (k+) and off (k) rates for the
binding of any ligand (L) to a specific site on the receptor monomer subunit (R), which was described by ordinary differential equations of the form:
![]() |
At equilibrium, the distribution of the receptor states is a function of the ligand dissociation constant (Kd = k/k+) and free
concentration of ligand. The four InsP3 receptor subunits comprising the channel complex were considered identical and positionally equivalent, and a state transition reaction scheme was
constructed for the 125 possible states of the channel tetramer
complex (T) (see RESULTS). The various Kd values fully define
the steady state distribution of states and allowed us to calculate
the relative abundance ti, c, j = [Ti, c, j]/[Ttotal] (of tetramers with i and j InsP3 molecules bound to the medium and low affinity site, respectively, and c Ca molecules bound) by solving a linear system of first degree equations that result from the equilibrium
condition for formation/degradation of each state. The no-flux
equilibrium condition, as well as several assumptions regarding
binding interactions, were used to constrain the various dissociation constants that were not directly inferred from binding experiments (see RESULTS). Channel states characterized by the
binding of one to two Ca and two to four InsP3 molecules were
assumed to be "active;" i.e., exhibit a high open probability. InsP3
binding to the low affinity site was assumed to create an active
state, independent of InsP3 and/or Ca binding to other sites
within the channel complex. Channel activity was described as a
transition between two substates, open and closed (Scheme I).
First order reaction rates were assumed for both transitions, and
the resulting open probability for any given state of the channel is:
![]() |
(SCHEME I)
Iterative curve-fitting of the InsP3-gated Ca channel open
probability data given the Ca dependence of InsP3 binding using a standard quasi-Newtonian algorithm with forward differencing and quadratic extrapolation was then performed to determine
the values of unknown parameters. The open time histograms
were used to further constrain the model and facilitate kinetic
simulations. Diffusion-limited on rates of binding were assumed
in all steps for both Ca and InsP3 binding and each off rate was
determined by multiplication with the respective dissociation
constant. Channel opening and closing rates were determined
from the open probabilities of the respective active states given
the constraints imposed on closing rates in order to match the
experimental distribution histogram of open states. Simulations
of the single channel recordings were performed by stochastic
modeling of state transition events of the channel complex, using a Markovian model of probability matrix decomposition
(Colquhoun and Hawkes, 1983) to predict state transitions over
a given time period. The behavior of the channel is calculated
from the matrix differential equation:
![]() |
where the elements pm,n(t) of the matrix P(t) are the probabilities of a channel to get from state m to n over the time period t and the elements of the matrix are the rates of transition from state m to n. The solution of this equation is:
![]() |
where I is a unit matrix. Relaxation from a random state (picked
according to steady state distribution probability) was simulated for the desired time length to obtain channel state data sets at
the chosen time step. Predicted traces were generated for the simultaneous presence of up to 10 channels in the membrane
from individual current functions: i(t) = A ·
(t) + S(t), where
A is the current amplitude through the open channel,
(t) is either 0 or 1, depending on whether the channel is in one of the
closed or one of the open states, respectively, and S(t) is a standard noise spectrum function chosen to match the experimentally observed noise in any given trace.
![]() |
RESULTS |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Ca Inhibition of InsP3 Binding
There have been conflicting reports on the influence
of Ca on InsP3 binding. Initial reports indicated that
micromolar Ca concentrations could completely inhibit InsP3 binding, with half-maximal inhibition occurring at 0.3 µM Ca (Worley et al., 1987; Supattapone et al.,
1988
). On the other hand, Ca-dependent inhibition of
cerebellar membranes can also occur as an artifact resulting from activation of an endogenous phospholipase C (Mignery et al., 1992
). We have re-investigated
the influence of Ca on InsP3 binding, with the goal of
relating Ca-dependent modulation of InsP3 binding to
the bell-shaped Ca dependence of channel activity. As
shown in Fig. 1 A, micromolar Ca decreased InsP3 binding when using media identical to that used for channel activity measurements. Ca concentrations as high as
100 µM resulted in only partial inhibition of InsP3 binding (Fig. 1 A, filled bars), in contrast to the near
complete inhibition previously reported (Worley et al.,
1987
). The inhibitory effects of Ca on InsP3 binding
were reversed by addition of EGTA (Fig. 1 A, hatched
bars), ruling out contributions of Ca-dependent activation of an endogenous phospholipase C (Mignery et al.,
1992
; Benevolensky et al., 1994
). Moreover, the presence of 100 µM free barium, an inhibitor of phospholipase C (Benevolensky et al., 1994
), did not alter this
Ca-dependent inhibition of InsP3 binding (data not
shown).
|
Scatchard analyses of InsP3 binding (Fig. 1 B) indicated that this decreased InsP3 binding was attributable
to a Ca-dependent increase in the Kd for InsP3 binding.
The Scatchard plots were nonlinear, suggesting the
presence of multiple InsP3 binding sites. Data (which
were obtained over an InsP3 concentration range of 1-600 nM) could be fit assuming two InsP3 binding
sites (with Kd values of 1 and 50 nM in the absence of
Ca). The concentration of InsP3 used was too low to observe the low affinity InsP3 binding site discussed below
(Kaftan et al., 1997). The 50-nM binding site was abundant, representing 99% of the total sites in the absence
of Ca. As the Ca concentration was increased from
0.0004 to 10 µM, the apparent Kd for this site increased
3.5-fold (Fig. 1, B and C). Half maximal inhibition
occurred at 0.3 µM Ca with a Hill coefficient of 1.5. The high affinity InsP3 binding site, by contrast, was in
very low abundance, representing <1% of the total sites
at all Ca concentrations examined, and its Kd (1 nM) appeared to be Ca independent (consistent with our
previous observations in aortic smooth muscle microsomes; Benevolensky et al., 1994
).
To more closely resemble the conditions used for
channel activity determinations, InsP3 binding assays
were repeated at room temperature again using media
identical to that employed in the single channel experiments. Increasing the temperature to 22°C resulted in a
67% reduction in InsP3 binding; an inhibition that was
reversed by reducing the temperature to 0°C (Fig. 2 A).
Scatchard analyses of InsP3 binding at 0° and 22°C (Fig.
2 B) show that the increase in temperature results in a
decrease in InsP3 binding affinity, with no apparent effect on the maximal number of InsP3 binding sites. The
calculated Q10 for this temperature-dependent change
in InsP3 binding affinity was 2.3, consistent with that calculated for InsP3 binding to rat basophilic leukemia
cells (Watras et al., 1994).
|
The Ca dependence of InsP3 binding was also reexamined at 22°C (Fig. 3), and yielded results similar to
those observed at 0°C (Fig. 1). As shown in Fig. 3 A, micromolar Ca concentrations reduced InsP3 binding
(measured in the presence of 170 nM InsP3) by 55%
(), and this inhibition was completely reversed by Ca chelation with EGTA (
, dotted line). Increasing the
Ca concentration to 100 µM did not increase the extent of the inhibition (data not shown), consistent with
our previous observations in aortic smooth muscle microsomes (Benevolensky et al., 1994
). The Ca-dependent inhibition at 22°C showed a high level of cooperativity (Hill coefficient, 2.8; half-maximal inhibition at
0.45 µM Ca) (Fig. 3 A), and was attributable to a Ca-
dependent increase in the Kd for InsP3 binding (Fig.
3 B). The Kd for InsP3 binding at 22°C increased from
0.30 (at 0.0005 µM Ca;
) to 1.78 µM in the presence of a high Ca concentration (
; 50 µM free Ca) (Fig. 3
B). There was no evidence for an effect of Ca on the
maximum number of InsP3 binding sites at either 0° or
22°C.
|
The Scatchard analyses in Figs. 2 and 3 show the presence of only one InsP3 binding site (Kd = 50 nM at 0°C
and 312 nM at 22°C) due to the range of InsP3 concentrations used (5 nM-1 µM), where neither the high (Kd = 1 nM) nor the low (Kd = 10 µM) affinity site can be resolved. These sites are observed only if one extends the
InsP3 concentration range below 1 nM (as in Fig. 1)
and above 5 µM (Kaftan et al., 1997). Three InsP3 binding sites were observed in cerebellar microsomes and
purified InsP3 receptor preparations, purified by either
heparin-agarose chromatography or immunoprecipitation (Kaftan et al., 1997
). In all assays, the 1-nM InsP3
binding site was a minor component, whereas the medium and low affinity InsP3 binding sites were abundant. At 0°C, the calculated dissociation constants for
InsP3 binding to these two predominant InsP3 binding sites were 50 nM (medium affinity site) and 10 µM (low
affinity site) (Kaftan et al., 1997
). At 22°C, the dissociation constant of the 50-nM site increased to ~300 nM,
whereas the dissociation constant for the 10-µM InsP3
binding site could not be accurately ascertained. It appeared to remain in the 10-µM range, suggesting that
the low affinity InsP3 binding site was largely temperature insensitive.
InsP3 Dependence of Channel Activity
We previously reported that the Hill coefficient for
InsP3-induced Ca release was 1.3, whereas the single
channel activity determinations suggested a Hill coefficient 1.0 (Watras et al., 1991). The previous experiments underestimated the Hill coefficient because the
lowest concentration of InsP3 used activated the channel to 10% of the maximum, a level of activity too high
for accurate measurements of cooperativity. Reinvestigation of the channel activity at lower InsP3 concentrations indicate that the Hill coefficient for channel activity measurements is at least 1.7 (Fig. 4). This value for
the Hill coefficient may reflect the requirement for
InsP3 to bind to at least two of the four InsP3 receptors
to activate the channel.
|
Open-Time Distribution of InsP3-gated Channels
Open-time histograms of InsP3-gated channels at two different cytoplasmic Ca concentrations were fit assuming an exponential with two time constants (Fig. 5). Increasing the number of time constants to three did not statistically improve the fit. In the presence of 0.01 µM Ca and 2 µM InsP3 (Fig. 5 A), the appropriate time constants were 1.8 and 9.3 ms. As the Ca concentration was increased to 0.1 µM (Fig. 5 B), the time constants decreased to 0.6 and 4.7 ms. Thus, the channel open times are Ca dependent. Similar analyses were done at a number of Ca and InsP3 concentrations, with the appropriate time constants incorporated into the mathematical model (described below).
|
Modeling Channel Activity Data and InsP3 Binding
Recent studies indicate that the cerebellar InsP3-gated
Ca channel exhibits a complex regulation by Ca and
InsP3 (Kaftan et al., 1997). Data points from this earlier
report are shown in Fig. 6; lines through the points represent the fit using the model described in the present
report (Fig. 6 A). For comparison, the fit obtained using the 16-state model (Kaftan et al., 1997
) is also
shown (Fig. 6 B), but fails to accurately predict channel activity at low InsP3 concentrations. The current model
provides an excellent fit of the channel data over a
broad range of InsP3 and Ca concentrations.
|
Simple mass-action kinetics was assumed for ligand binding to the sites on the receptor monomer subunits (see METHODS). Since the functional unit of the channel is the tetrameric complex (T), this analysis was extended by introducing apparent dissociation constants for sequential binding to the four monomer subunits. A total of 24 = 16 possible different states of the channel can be generated by binding of up to four ligand molecules to the same site on the four different monomer subunits of the tetramer. We considered the subunits identical and positionally equivalent, which reduces the 16 states to only 5 different states according to Scheme II.
|
The resulting apparent dissociation constants for sequential binding are:
![]() |
This general framework was applied to InsP3 and Ca binding to the channel complex, assuming that each monomer in the tetrameric channel complex contains the equivalent of one Ca regulatory site (C), one medium affinity InsP3 site (I), and one low affinity InsP3 site (J). This results in five different states of the tetrameric channel complex (zero, one, two, three, or four sites occupied) for each of the three classes of sites (C, I, and J). A total of 53 = 125 different states per channel complex are predicted (Ti,c,j, where i, c, and j represent the number (zero to four) of ligand molecules bound to the I, C, and J sites, respectively). Fig. 7 shows the state transition scheme for the 25 states in the absence of InsP3 binding to the low affinity site (Ti,c,0). The apparent dissociation constants of binding of InsP3 to the I site were labeled K'Idi, c (for a tetramer that has i other I sites and c of the C sites occupied) and of Ca binding to the C site were labeled K'Cdi, c (for a tetramer that has c other C sites and i of the I sites occupied). This reaction scheme was extended upwards by drawing four additional planes to include binding to the low affinity InsP3 binding site (J).
|
The "gating" or "regulatory" action of ligand binding
to the various sites on the channel is in effect determined by the differences in the opening/closing transition rates of the different states (see METHODS; i.e., a
state with
0 is an "inactive" state). Essentially, "active" states were hypothesized to be channels having either two or more medium affinity InsP3 sites and one or
two Ca sites occupied (Fig. 7, highlighted states; a total of eight active states in the absence of binding to the J
site). Additionally, active states were channels having
two or more low affinity InsP3 sites (J) and any number
of Ca sites occupied (in upper planes of the full 125-state model, not shown in Fig. 7). The activation of the
channel through J was thus assumed to overcome Ca-dependent inhibition of the channel, consistent with
experimental observations (Kaftan et al., 1997
).
At any given Ca and InsP3 concentration, the only parameters needed to calculate the steady state relative distribution of the different states and the overall channel open probability were the various apparent dissociation constants
![]() |
and the open probabilities of the various statesPopeni, c, j
. Some of these parameters could be inferred from measured dissociation constants based on several assumptions
and constraints. In accordance with the thermodynamic principle of detailed equilibrium (Hille, 1992), reaction
rates between any combination of states that form a complete cycle (e.g., the highlighted set of four transitions in
Fig. 7) were constrained to yield no net flux at equilibrium. Thus, the apparent dissociation constants (K'd)
of
binding to the C and I site were constrained as follows:
![]() |
With regard to interactions when binding to the same
site on different monomer subunits, we assumed that
Ca binding is cooperative (KCdc, i>KCdc+1, i), whereas InsP3
binding is independent KCdi, c=KCdi+1, c ; for simplicity,
we assumed that the cooperative effect is the same for
each step during sequential binding of Ca; i.e., the ratio KCdc, i>KCdc+1, i is the same for c = 0...3 (the value of
this ratio is a parameter we named cooperativity factor,
Ca). The dissociation constants for InsP3 were derived
directly from the binding data (Fig. 1 C, and Table I).
For the I site, we used K'Idi, 0 = 0.3 µM (no Ca bound)
and K'Idi, c
= 1.05 µM (for c > 0, reflecting the noncompetitive inhibition of InsP3 binding by Ca). Binding
of InsP3 to the low affinity site appears independent of
binding to any other sites with a dissociation constant
KJd
= 10 µM. Given the above constraints and InsP3
binding parameters, all the dissociation constants for
Ca binding could be derived from only two parameters,
KCd0, 0 and
Ca. The values of these as well as of Popeni, c, j were obtained through standard curve-fitting algorithms for the predicted channel open probability data
(Table I). The predicted Kd of Ca binding to the first
site is 0.7 µM, which approximates the reported Kd for
Ca binding to a cytoplasmic region of the type 1 InsP3
receptor (Sienaert et al., 1996
). Table I summarizes all
the parameters used to construct the model.
|
Simulation of Channel Kinetics
A sample recording of InsP3-gated channel activity obtained in the presence of 0.2 µM InsP3 and 0.1 µM Ca (Fig. 8 A) is shown for comparison with channel activity predicted by the model (B). Using constants derived from the InsP3 binding analyses, channel open time histograms, and the steady state open probability data, the mathematical model described above was used to generate a time course of channel activity in the presence of 0.2 µM InsP3 and 0.1 µM Ca. The predicted channel activity data (Fig. 8 B) closely resemble the actual channel recording (A); moreover, quantitative comparisons show that the dwell time histogram of a simulated trace (Fig. 8 C) could be fitted with a double exponential curve using identical time constants as those used to fit the experimental data obtained at the same Ca and InsP3 concentrations (Fig. 5 B).
|
The parameter list used for the steady state analysis
(and values) is not sufficient to calculate the values required to fill the transition rate matrix (see METHODS).
This requires explicit numbers for the on and off rates
(k+ and k values) and opening and closing rates (
and
values). The on rates for binding were based assuming diffusion limitation for both Ca and InsP3 (10 and 250 µM
1 s
1, respectively; this includes the effect
of buffering for Ca) binding. The channel closing rates
were based on the single channel open time distribution histograms. A minimum of two different closing
rates is demonstrated by the biexponential fits of the dwell times. To accurately predict the apparent time
constants and relative ratios, as well as their change
with increasing Ca concentration (Fig. 5), we used four
different closing rates ranging from 150 to 1,000 s
1
(Table I) whereby the faster rates were assigned to the
active states with two Ca bound and the slower rates to
the active states with one Ca bound. The corresponding binding off rates and channel opening rates were
then calculated using the parameter values from the
steady state simulations:
![]() |
To maintain the necessary matrix algebra at a manageable size, we limited ourselves to simulations of channel
activity under conditions with [InsP3] 2 µM, so that
the effect of the J site on channel activity is negligible
and only the 25 states depicted in Fig. 7 (and the corresponding substates of the active states) have to be included in the calculations.
To produce simulated traces of channel current recording, we generated predicted relaxation data. An
initial state was chosen for time 0, either arbitrarily or
the most probable state for the given Ca and InsP3 concentration (as indicated by the steady state distribution
values). We calculated the probability matrix P(t) for a
time value equal to the chosen sampling interval in single channel records ( = 40 µs), and generated data
sets of predicted channel states at the given sampling
rate for a desired length of simulation. Fig. 8 B shows a
representative segment of a simulated trace with 2.1 pA
open channel current amplitude and normally distributed baseline noise of 1 kHz median frequency.
The current mathematical model is thus capable of modeling both steady state open probability data, ligand binding data, and channel kinetics.
![]() |
DISCUSSION |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Recent studies in our laboratory indicate that Ca and
InsP3 interact to increase the dynamic range of the
InsP3-gated Ca channel (Kaftan et al., 1997). At a resting intracellular Ca concentration of 0.05-0.1 µM, agonist-induced production of InsP3 would be expected to
activate the InsP3-gated Ca channel, with a subsequent
rise of intracellular Ca (up to 0.3 µM) promoting further activation of the channel through a positive feedback mechanism (Iino, 1990
; Bezprozvanny et al., 1991
;
Finch et al., 1991
). Negative feedback of Ca on the
InsP3-gated channel, however, will limit the rise of intracellular Ca, with near complete inhibition of the channel at 5 µM Ca, even in the presence of 2 µM InsP3
(Bezprozvanny et al., 1991
). Increasing the InsP3 concentration results in a rightward shift in the bell-shaped
curve (Kaftan et al., 1997
). Thus, the channel is more
sensitive to the inhibitory effects of Ca when the InsP3
concentration is low. This may be important physiologically, as it would further limit the rise in intracellular Ca
in the presence of low InsP3 concentrations. That is, it
would promote a more gradual rise in intracellular
Ca concentration over the InsP3 concentration range
0.02-2 µM. Conversely, very high InsP3 concentrations
can largely overcome this Ca-dependent inhibition of
the channel (Kaftan et al., 1997
), preventing closure of
the channel in local/confined areas such as those near
the site of InsP3 production.
An important question concerns the mechanism(s)
responsible for this complex regulation of the InsP3-gated channel. Our working hypothesis has been that
the previously reported Ca-dependent decrease in
InsP3 binding is involved. Although micromolar Ca has
been reported to completely inhibit InsP3 binding
to cerebellar membranes (Worley et al., 1987; Supattapone et al., 1988
), our results reveal only a partial
(<50%) inhibition of InsP3 binding. As 95-99% of the
InsP3 receptors in the cerebellum have been shown to
be of the type 1 isoform (Sudhof et al., 1991
; Wojcikiewicz, 1995
), it is unlikely that this partial inhibition
is due to heterogeneity of InsP3 receptors. Instead, the
Ca-dependent decrease in InsP3 binding appears to reflect a change in the affinity for InsP3. This was first suggested in a report showing that the Kd for InsP3 binding
to rat cerebellar membranes increased from 0.04 to
0.12 µM (Joseph et al., 1989
). Such a shift in InsP3 affinity has also been observed in aortic smooth muscle
membranes (Benevolensky et al., 1994
). Using binding
conditions identical to those used for the channel activity determinations, the present study also shows a Ca-dependent shift in InsP3 affinity, though the calculated Kds are substantially higher than those previously reported. Thus, the binding affinities for InsP3 reported
in the present study are comparable to the calculated
InsP3 affinity for channel activation. The reason for the
relatively high Kd values in the present study involves in
large part the high temperature dependence of the binding constant (Q10 = 2.3 for the cerebellar membranes). That is, InsP3 binding studies are typically
done on ice (0°C), whereas channel activity determinations are done at room temperature (23°C). In the
present study, however, efforts were made to closely
match both the InsP3 binding and channel analyses. A
similar Q10 (2.0) has been reported for the Kd of InsP3
binding to rat basophilic leukemia cells (Watras et al.,
1994
), which was useful for modeling Ca release using
InsP3 binding data.
The eight-state model originally proposed for the
regulation of the InsP3-gated Ca channel (De Young
and Keizer, 1992) represents one possible explanation
for the bell-shaped Ca dependence of channel activity
at low InsP3 concentrations. Yet this eight-state model is
unable to explain the relief of Ca-dependent inhibition of InsP3-gated channel activity at high InsP3 concentrations. We hypothesized that a low affinity InsP3 binding
site may be involved, and provided evidence for its existence using ligand binding assays. To fit the data at
high InsP3 concentrations, the eight-state model was
therefore extended to include the 10 µM InsP3 binding
site, yielding a 16-state model (Kaftan et al., 1997
).
In an effort to both explain the steady state channel
data and extract kinetic information on InsP3-gated Ca
release, we used a novel approach that involves consideration of the channel as a unit. Previous models have
focused on an individual InsP3 receptor, predicting various states of the monomer, and then calculating probabilities of each combination of states in a tetrameric
complex. The latter approach was used when extending the previously reported eight-state model of the
InsP3 receptor (De Young and Keizer, 1992) to 16 states
(Kaftan et al., 1997
), and we obtained reasonable fits of
the channel data. A stochastic model has recently been proposed (Swillens et al., 1998
), though it assumed
only one InsP3 binding site per channel complex,
which is likely an oversimplification. The present
model provides a more realistic representation of the
channel with four medium-affinity InsP3 binding sites
per channel complex (Maeda et al., 1988
; Supattapone
et al., 1988
; Chadwick et al., 1990
), and yields superior
fits of the data over a broad range of InsP3 and Ca concentrations compared with our previous 16-state model
(Kaftan et al., 1997
; see Fig. 8 A). Moreover, rate constants (rather than dissociation constants) and mean
open times were incorporated into the current model,
which facilitates kinetic analyses of channel activity. By
focusing on the channel complex (rather than individual InsP3-receptor subunits), the number of possible
states of the channel could be greatly reduced, thereby
allowing for the kinetic simulations of channel activity.
The mathematical model in the present study requires just four Ca binding sites per tetrameric channel
complex. This is consistent with our initial observation
(Bezprozvanny et al., 1991) that in the presence of 2 µM
InsP3 the bell-shaped Ca dependence of channel activity could be fit assuming that two molecules of Ca were
needed to open the channel, whereas two molecules of
Ca appeared to close the channel. A similar stoichiometry of four divalent cation sites per channel complex
was obtained when Mn was substituted for Ca (Striggow
and Ehrlich, 1996
). Examination of Ca binding to the
cloned/expressed receptor fragments raises the possibility of up to eight Ca binding sites per InsP3 receptor
monomer (Mignery et al., 1992
; Sienaert et al., 1996
,
1997
). Seven of these Ca binding sites appear to reside
on the cytosolic region of the InsP3 receptor (Sienaert
et al., 1997
), though the Ca affinity of only one of these
sites (Kd 0.8 µM Ca, Hill coefficient 1.0) has been determined (Sienaert et al., 1996
). The predicted Kd for
Ca binding in the present study (0.7 µM, Table I) is
consistent with the latter report. It has been hypothesized, based on the cation selectivity of the stimulation
and inhibition of InsP3-induced Ca release (Marshall
and Taylor, 1994
; Sienaert et al., 1996
; Striggow and
Ehrlich, 1996
), that Ca binding to at least three distinct
cytosolic sites regulates InsP3-gated channel activity (Sienaert et al., 1997
). For simplicity of the mathematical model, it is assumed that there are four distinct Ca
binding sites on the channel, with Ca binding to one
site altering the Ca affinity of the other sites in the
channel complex. Inclusion of multiple Ca binding
sites per InsP3 receptor monomer in the model is possible, and would be expected to increase the freedom in
the model and provide even better fits of the channel
and binding data, but at the cost of a considerably
larger matrix calculation. As the identities of regulatory
Ca binding sites become available, the mathematical
model will be refined to more accurately represent the
molecular details of channel structure.
Although the proposed model provides an accurate
simulation of steady state InsP3-gated channel activity
and InsP3 binding data, it is important to ascertain the
physiological relevance of an InsP3 binding site with a
dissociation constant of ~10 µM, particularly if the medium affinity InsP3 binding site has an affinity 30×
lower (0.3 µM) at 22°C. It should be appreciated, however, that the medium affinity InsP3 binding site exhibited a strong temperature dependence (Q10 = 2.3), so
that at 37°C, the dissociation constant of the medium
affinity InsP3 binding site is predicted to approach 1 µM.
The temperature dependence of the low affinity InsP3
binding site was difficult to measure, but it appeared to
be relatively temperature independent. There are reports of ligand binding sites with low (or reversed) temperature dependence, where binding affinity is mainly
entropy driven; it is not possible to accurately predict a
priori if a given binding site will have a high or low temperature dependence of ligand binding (Testa et al.,
1987). The low level of specificity of the low affinity site
for 1,4,5-InsP3 versus 1,3,4,5-InsP4 (Kaftan et al., 1997
)
is in part consistent with the apparent low temperature dependence of this site. Thus, the difference in dissociation constants between the medium and low affinity
InsP3 binding sites may be as low as 10-fold at the
physiological temperature of 37°C (i.e., 1 vs. 10 µM).
At first glance, dissociation constants of 1-10 µM
seem to be unexpectedly high, although review of the
literature indicates that at physiological temperatures,
rather high InsP3 concentrations may be needed to release Ca (Putney, 1990). In the case of cerebellar Purkinje cells, flash photolysis studies indicated the
need for at least 9 µM intracellular InsP3 to induce Ca
release, with Ca release increasing progressively as
InsP3 concentration was elevated up to 80 µM (Khodakhah and Ogden, 1995
). Moreover, it has been reported that sustained elevation of intracellular Ca can
be observed using high concentrations (100 µM) of either InsP3 or the poorly hydrolyzed InsP3 analogue
InsPS3 (Wakui et al., 1989
; Petersen et al., 1991
). Furthermore, relatively high InsP3 concentrations have
been noted in some cell types under basal (0.1-3 µM)
and agonist-induced (1-20 µM) conditions (Putney,
1990
). Thus, InsP3 dissociation constants of 1-10 µM at
37°C may be of physiological significance.
In summary, the present study provides a simple mathematical model of the cerebellar InsP3-gated channel, which provides an accurate simulation of InsP3-gated channel activity over a broad range of InsP3 and Ca concentrations. The model incorporates the well known "medium affinity" InsP3 binding site, and takes into account its temperature and Ca dependences. The model also includes a low affinity InsP3 binding site that exhibits little Ca or temperature dependence. Individual rate constants were incorporated into the model to allow kinetic simulations. It is expected that this novel model will provide a foundation for subsequent studies aimed at elucidating the molecular mechanisms underlying the complex behavior of the InsP3-gated channel.
![]() |
FOOTNOTES |
---|
Original version received 17 June 1998 and accepted version received 19 March 1999.
The authors thank William Dyckman for his assistance in obtaining canine cerebella and Dr. Leslie Loew for valuable comments on the kinetic model.
This work was supported by National Institutes of Health grants HL-33026 (J. Watras and B.E. Ehrlich), GM-51480 (B.E. Ehrlich), RR13186 (J. Watras, I.I. Moraru), and a grant from the University of Connecticut Health Center (J. Watras).
![]() |
REFERENCES |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
1. | Abdel-Latif, A.. 1986. Calcium-mobilizing receptors, polyphosphoinositides, and the generation of second messengers. Pharmacol. Rev. 38: 227-272 [Medline]. |
2. | Benevolensky, D., I. Moraru, and J. Watras. 1994. Micromolar calcium reduces the affinity of the inositol 1,4,5-trisphosphate receptor in smooth muscle. Biochem. J. 299: 631-636 [Medline]. |
3. | Berridge, M.J., and R.F. Irvine. 1989. Inositol phosphates and cell signalling. Nature. 341: 197-205 [Medline]. |
4. | Bezprozvanny, I.. 1994. Theoretical analysis of calcium wave propagation based on inositol (1,4,5)-trisphosphate (InsP3) receptor functional properties. Cell Calc. 16: 151-166 [Medline]. |
5. | Bezprozvanny, I., J. Watras, and B.E. Ehrlich. 1991. Bell-shaped calcium-response curves of Ins(1,4,5)P3- and calcium-gated channels from endoplasmic reticulum of cerebellum. Nature. 351: 751-754 [Medline]. |
6. | Bootman, M.D., L. Missiaen, J.B. Parys, H. DeSmedt, and R. Casteels. 1995. Control of inositol 1,4,5-trisphosphate-induced Ca2+ release by cytosolic Ca2+. Biochem. J. 306: 445-451 [Medline]. |
7. | Bradford, M.M.. 1976. A rapid and sensitive method for the quantitation of microgram quantities of protein utilizing the principle of protein-dye binding. Anal. Biochem. 72: 248-254 [Medline]. |
8. | Chadwick, C.C., A. Saito, and S. Fleischer. 1990. Isolation and characterization of the inositol trisphosphate receptor from smooth muscle. Proc. Natl. Acad. Sci. USA. 87: 2132-2136 [Abstract]. |
9. | Clapham, D.E.. 1995. Calcium signaling. Cell. 80: 259-268 [Medline]. |
10. | Colquhoun, D., and A.G. Hawkes. 1983. The principles of the stochastic interpretation of ion-channel mechanisms. In Single Channel Recording. Plenum Publishing Corp., New York. 135-176. |
11. | De Young, G.W., and J. Keizer. 1992. A single-pool inositol 1,4,5-trisphosphate-receptor-based model for agonist-stimulated oscillations in Ca concentration. Proc. Natl. Acad. Sci. USA. 89: 9895-9899 [Abstract]. |
12. | Fabiato, A.. 1988. Computer programs for calculating total from specified free or free from specified total ionic concentrations in aqueous solutions containing multiple metals and ligands. Methods Enzymol. 157: 378-417 [Medline]. |
13. | Finch, E.A., T.J. Turner, and S.M. Goldin. 1991. Calcium as a coagonist of inositol 1,4,5-trisphosphate-induced calcium release. Science. 252: 443-446 [Medline]. |
14. | Hannaert-Merah, Z., L. Combettes, J.-F. Coquil, S. Swillens, J.-P. Mauger, M. Claret, and P. Champeil. 1995. Characterization of the co-agonist effects of strontium and calcium on myo-inositol trisphosphate-dependent ion fluxes in cerebellar microsomes. Cell Calc. 18: 390-399 [Medline]. |
15. | Hille, B. 1992. Gating mechanisms. In Ionic Channels of Excitable Membranes. Sinauer Associates, Inc., Sunderland, MA. 472-504. |
16. | Horn, R.. 1991. Estimating the number of channels in patch recordings. Biophys. J. 60: 433-439 . |
17. | Iino, M.. 1990. Biphasic Ca2+ dependence of inositol 1,4,5-trisphosphate-induced Ca release in smooth muscle cells of the guinea pig taenia caeci. J. Gen. Physiol. 95: 1103-1122 . |
18. | Joseph, S.K., H.L. Rice, and J.R. Williamson. 1989. The effect of external calcium and pH on inositol trisphosphate-mediated calcium release from cerebellum microsomal fractions. Biochem. J. 258: 261-265 [Medline]. |
19. |
Kaftan, E.J.,
B.E. Ehrlich, and
J. Watras.
1997.
Inositol 1,4,5-trisphosphate (InsP3) and calcium interact to increase the dynamic
range of InsP3 receptor-dependent calcium signaling.
J. Gen.
Physiol.
110:
529-538
|
20. | Khodakhah, K., and D. Ogden. 1995. Fast activation and inactivation of inositol trisphosphate-evoked Ca2+ release in rat cerebellar Purkinje neurones. J. Physiol. 487: 343-358 [Abstract]. |
21. | Maeda, N., M. Niinobe, K. Nakahira, and K. Mikoshiba. 1988. Purification and characterization of P400 protein, a glycoprotein characteristic of Purkinje cell, from mouse cerebellum. J. Neurochem. 51: 1724-1730 [Medline]. |
22. | Marshall, I., and C. Taylor. 1994. Two calcium binding sites mediate the interconversion of liver inositol 1,4,5-trisphosphate receptors between three conformations states. Biochem. J. 301: 591-598 [Medline]. |
23. |
Mignery, G.A.,
P.A. Johnston, and
T.C. Sudhof.
1992.
Mechanism of
Ca2+ inhibition of inositol 1,4,5-trisphosphate (InsP3) binding to
the cerebellar InsP3 receptor.
J. Biol. Chem.
267:
7450-7455
|
24. | Payne, R., T.M. Flores, and A. Fein. 1990. Feedback inhibition by calcium limits the release of calcium by inositol trisphosphate in Limulus ventral photoreceptors. Neuron. 4: 547-555 [Medline]. |
25. | Petersen, C.C.H., E.C. Toescu, B.V.L. Potter, and O.H. Petersen. 1991. Inositol triphosphate produces different patterns of cytoplasmic Ca2+ spiking depending on its concentration. FEBS Lett. 293: 179-182 [Medline]. |
26. | Putney, J.. 1990. The integration of receptor-regulated intracellular calcium release and calcium entry across the plasma membrane. Curr. Top. Cell Regul. 31: 111-127 [Medline]. |
27. |
Sienaert, I.,
H. De Smedt,
J.B. Parys,
L. Missiaen,
S. Vanlingen,
H. Sipma, and
R. Casteels.
1996.
Characterization of a cytosolic and
a luminal Ca2+ binding site in the type I inositol 1,4,5-trisphosphate receptor.
J. Biol. Chem.
271:
27005-27012
|
28. |
Sienaert, I.,
L. Missiaen,
H. De Smedt,
J.B. Parys,
H. Sipma, and
R. Casteels.
1997.
Molecular and functional evidence for multiple
Ca2+-binding domains in the type 1 inositol 1,4,5-trisphosphate
receptor.
J. Biol. Chem.
272:
25899-25906
|
29. | Striggow, F., and B.E. Ehrlich. 1996. The inositol 1,4,5-trisphosphate receptor of cerebellum. Mn2+ permeability and regulation by cytosolic Mn2+. J. Gen. Physiol. 108: 115-124 [Abstract]. |
30. | Sudhof, T.C., C.L. Newton, B.T. Archer III, Y.A. Ushkaryov, and G.A. Mignery. 1991. Structure of a novel InsP3 receptor. EMBO (Eur. Mol. Biol. Organ.) J. 10: 3199-3206 [Abstract]. |
31. |
Supattapone, S.,
P.F. Worley,
J.M. Baraban, and
S.H. Snyder.
1988.
Solubilization, purification, and characterization of an inositol
trisphosphate receptor.
J. Biol. Chem.
263:
1530-1534
|
32. | Swillens, S., P. Champeil, L. Combettes, and G. Dupont. 1998. Stochastic simulation of a single inositol 1,4,5-trisphosphate-sensitive Ca2+ channel reveals repetitive openings during `blip-like' Ca2+ transients. Cell Calc. 23: 291-302 [Medline]. |
33. | Tang, Y., J.L. Stephenson, and H.G. Othmer. 1996. Simplification and analysis of models of calcium dynamics based on IP3-sensitive calcium channel kinetics. Biophys. J. 70: 246-263 [Abstract]. |
34. | Testa, B., P. Jenner, G.J. Kilpatrick, N. el Tayar, H. Van de Waterbeemd, and C.D. Marsden. 1987. Do thermodynamic studies provide information on both the binding to and the activation of dopaminergic and other receptors? Biochem. Pharmacol. 36:4041- 4046. |
35. | Wakui, M., B.V.L. Potter, and O.H. Petersen. 1989. Pulsatile intracellular calcium release does not depend on fluctuations in inositol trisphosphate concentration. Nature. 339: 317-320 [Medline]. |
36. |
Watras, J.,
I. Bezprozvanny, and
B.E. Ehrlich.
1991.
Inositol 1,4,5-trisphosphate-gated channels in cerebellum![]() |
37. | Watras, J., I. Moraru, D.J. Costa, and L.A. Kindman. 1994. Two inositol 1,4,5-trisphosphate binding sites in rat basophilic leukemia cells: relationship between receptor occupancy and calcium release. Biochemistry. 33: 14359-14367 [Medline]. |
38. |
Wojcikiewicz, R.J..
1995.
Type I, II, and III inositol 1,4,5-trisphosphate receptors are unequally susceptible to down-regulation
and are expressed in markedly different proportions in different
cell types.
J. Biol. Chem.
270:
11678-11683
|
39. |
Worley, P.F.,
J.M. Baraban,
S. Supattapone,
V.S. Wilson, and
S.H. Snyder.
1987.
Characterization of inositol trisphosphate receptor
binding in brain. Regulation by pH and calcium.
J. Biol. Chem.
262:
12132-12136
|