Home | Help | Feedback | Subscriptions | Archive | Search | Table of Contents

This Article Abstract PDF (Full Text) PPT slides of all figures Supplemental Material Index Alert me when this article is cited Citation Map Services Email this article Similar articles in this journal Similar articles in PubMed Alert me to new content in the JGP Download to citation manager Citing Articles Citing Articles via HighWire Citing Articles via CrossRef Citing Articles via Google Scholar Google Scholar Articles by Yamashita, M. Articles by Prakriya, M. Search for Related Content PubMed PubMed Citation Articles by Yamashita, M. Articles by Prakriya, M. Pubmed/NCBI databases Gene GEO Profiles HomoloGene UniGene Compound via MeSH Substance via MeSH Hazardous Substances DB CALCIUM COMPOUNDS CALCIUM, ELEMENTAL Social Bookmarking                      What's this?

Department of Molecular Pharmacology and Biological Chemistry, Northwestern University Feinberg School of Medicine, Chicago, IL 60611

Correspondence to Murali Prakriya: m-prakriya{at}northwestern.edu

TOP ABSTRACT INTRODUCTION MATERIALS AND METHODS RESULTS DISCUSSION REFERENCES

 
Ca²⁺ entry through store-operated Ca²⁺ release-activated Ca²⁺ (CRAC) channels is an essential trigger for lymphocyte activation and proliferation. The recent identification of Orai1 as a key CRAC channel pore subunit paves the way for understanding the molecular basis of Ca²⁺ selectivity, ion permeation, and regulation of CRAC channels. Previous Orai1 mutagenesis studies have indicated that a set of conserved acidic amino acids in trans membrane domains I and III and in the I–II loop (E106, E190, D110, D112, D114) are essential for the CRAC channel's high Ca²⁺ selectivity. To further dissect the contribution of Orai1 domains important for ion permeation and channel gating, we examined the role of these conserved acidic residues on pore geometry, properties of Ca²⁺ block, and channel regulation by Ca²⁺. We find that alteration of the acidic residues lowers Ca²⁺ selectivity and results in striking increases in Cs⁺ permeation. This is likely the result of enlargement of the unusually narrow pore of the CRAC channel, thus relieving steric hindrance for Cs⁺ permeation. Ca²⁺ binding to the selectivity filter appears to be primarily affected by changes in the apparent on-rate, consistent with a rate-limiting barrier for Ca²⁺ binding. Unexpectedly, the mutations diminish Ca²⁺-mediated fast inactivation, a key mode of CRAC channel regulation. The decrease in fast inactivation in the mutant channels correlates with the decrease in Ca²⁺ selectivity, increase in Cs⁺ permeability, and enlargement of the pore. We propose that the structural elements involved in ion permeation overlap with those involved in the gating of CRAC channels.

	INTRODUCTION

TOP ABSTRACT INTRODUCTION MATERIALS AND METHODS RESULTS DISCUSSION REFERENCES

 
Depletion of Ca²⁺ from the ER following stimulation of cell surface receptors produces Ca²⁺ influx through store-operated channels (SOCs) in a wide variety of cells. One well-characterized SOC is the calcium release-activated Ca²⁺ (CRAC) channel, which is expressed by mast cells, T cells, and related hematopoetic cells. CRAC channels mediate important biological functions including gene expression during T cell activation, release of histamine and serotonin from mast cells, and the exocytosis of lytic granules from cytotoxic T cells during target cell killing (Lewis, 2001; Parekh and Putney, 2005). In light of the critical roles of CRAC channels for cellular Ca²⁺ signaling and the devastating immunodeficiencies that result from the loss of their activity in human patients (Partiseti et al., 1994; Feske et al., 2005, 2006), there is considerable interest in understanding the phenotypic properties, modes of regulation, and the structural bases of CRAC channel function.

Electrophysiological studies have provided a wealth of information on the CRAC channel's pore properties and point to a unique biophysical fingerprint (for reviews see Lewis, 1999; Prakriya and Lewis, 2003; Parekh and Putney, 2005). The key characteristics of this fingerprint include (1) an extraordinarily high selectivity for Ca²⁺ over monovalent ions (P_Ca/P_Na > 1,000). As in voltage-gated Ca²⁺ (Ca_v) channels, this appears to arise from a high affinity Ca²⁺ binding site within the pore whose properties are finely honed to block permeation of monovalent ions, (2) a relatively narrow pore size of 3.9 Å (Prakriya and Lewis, 2006), which is in striking contrast to the much larger pores of Ca_v channels (Cataldi et al., 2002), (3) a very low permeability to Cs⁺ (P_Cs/P_Na0.1) that is also in contrast to the higher Cs⁺ permeability of Ca_v channels (Hess et al., 1986) and Ca²⁺-selective TRP channels (Voets et al., 2001), (4) an extremely low unitary conductance for both divalent and monovalent ions, and (5) multiple modes of modulation by Ca²⁺. The recent identification of the Orai family of proteins, thought to be the pore-forming subunits of store-operated channels (Feske et al., 2006; Vig et al., 2006b; Zhang et al., 2006), and the expression of functional CRAC currents from cloned Orai1 provide the necessary tools to elucidate the molecular underpinnings of these unique characteristics.

The Orai family consists of three closely conserved cell surface proteins (Orai1–3), each of which contains four predicted transmembrane domains with intracellular N and C termini (Feske et al., 2006; Prakriya et al., 2006). There is considerable evidence that Orai1 is a key subunit of the CRAC channel. First, a mutation in Orai1 (R91W) that causes a severe immunodeficiency in human patients eliminates I_CRAC in T cells (Feske et al., 2006). Second, overexpression of Orai1 together with STIM1 in HEK293 cells generates a large Ca²⁺ current similar to native I_CRAC in its biophysical and pharmacological profile (Mercer et al., 2006; Peinelt et al., 2006). Third, Orai1 contains a set of acidic residues that when altered greatly diminishes the Ca²⁺ selectivity of CRAC channels (Prakriya et al., 2006; Vig et al., 2006a; Yeromin et al., 2006). This is reminiscent of Ca_v channels, whose selectivity filters contain glutamate residues that are essential for high Ca²⁺ selectivity.

The acidic residues examined in the recent mutational studies include glutamates at positions 106 and 190 in TM1 and TM3 and aspartates at positions 110, 112, and 114 in the linker region between TM1 and TM2 (Prakriya et al., 2006; Vig et al., 2006a; Yeromin et al., 2006). By analogy to Ca_v channels, it is speculated that these residues form the key elements of the CRAC channel selectivity filter with the carboxylate-bearing side chains projecting into the pore to form one or more Ca²⁺-binding sites (Prakriya et al., 2006; Vig et al., 2006a; Yeromin et al., 2006). However, in the absence of knowledge of the number of subunits that make up the pore and lack of any structural information on the channel, a number of fundamental uncertainties persist about the characteristics of the selectivity filter. These uncertainties include the locations of the acidic residues within the pore and their role in shaping the overall pore architecture, the relationship between the pore diameter and ion permeation, and the affinity and kinetic properties of the Ca²⁺ binding site(s).

Much less is known about the molecular features of CRAC channel gating. In addition to its dependence on the luminal ER [Ca²⁺], the activity of CRAC channels is regulated at multiple levels. One mode of regulation, termed fast inactivation, arises from feedback inhibition of I_CRAC during brief hyperpolarizing steps and is influenced by the local [Ca²⁺]_i around individual CRAC channels (Hoth and Penner, 1993; Zweifach and Lewis, 1995a; Fierro and Parekh, 1999). Although the inactivation Ca²⁺-binding site has been mapped to the cytoplasmic face of the channel in close proximity to the pore (Zweifach and Lewis, 1995a), important details of the inactivation mechanism such as the molecular identities of the inactivation gate and the Ca²⁺ binding site(s), and the mechanisms that transduce Ca²⁺ binding to closure of the inactivation gate remain unknown.

In this study, we show that the same Orai1 acidic residues implicated in the control of Ca²⁺ selectivity are also important for regulating the CRAC channel pore geometry, Ca²⁺ block, and channel gating. We find that the loss of Ca²⁺ selectivity that results from alteration of the acidic residues is accompanied by striking increases in the CRAC channel pore diameter and Cs⁺ permeability, suggesting that the architecture of the pore as shaped by these residues strongly influences ion permeation. The mutations do not affect the CRAC channel's requirement for store depletion for activation, nor do they affect the spatial redistribution of Orai1 localization that occurs following store depletion, suggesting that the overall Orai1 structure is intact. Unexpectedly, alteration of the acidic residues also produces striking loss of Ca²⁺-dependent fast inactivation. These results suggest that the structural determinants of CRAC channel gating overlap with those involved in ion permeation, and raise the possibility that the opening and closing of the channel is regulated by the selectivity filter itself.

	MATERIALS AND METHODS

TOP ABSTRACT INTRODUCTION MATERIALS AND METHODS RESULTS DISCUSSION REFERENCES

 
Cells
HEK293 cells were grown in a medium consisting of 44% DMEM (Mediatech) and 44% Ham's F12 (Mediatech), supplemented with 10% FCS (HyClone), 1% 200 mM glutamine, 1% 5,000 U/ml penicillin, and 5,000 µg/ml streptomycin. The cells were maintained in log-phase growth at 37°C in 5% CO₂.

In some experiments, HEK293 cells grown in a medium consisting of CD293 (Invitrogen) supplemented with 2% 200 mM GlutaMAX (Invitrogen) were used. These cells were grown in suspension at 37° in 5% CO₂, plated onto poly-L-lysined coverslips at the time of passage, and grown in FBS-supplemented media (described above) until the time of transfection (24–48 h later).

Plasmids and Transfections
The IRES-eGFP-Orai1 plasmid was provided by S. Feske (Harvard University, Boston, MA). Site-directed mutagenesis to generate the Orai1 mutants was performed using the QuikChange Site-Directed Mutagenesis Kit (Stratagene) according to manufacturer's instructions. The following mutations were generated: E106D, E190Q, and the D110/112/114A triple mutant. HEK293 cells were transfected with Lipofectamine 2000 (Invitrogen), with 100 ng of Orai1 and 500 ng STIM1 per 12-mm coverslip. Cells were used for electrophysiology 24–48 h after transfection.

Solutions and Chemicals
The standard extracellular Ringer's solution contained (in mM) 135 NaCl, 4.5 KCl, 20 CaCl₂, 1 MgCl₂, 10 D-glucose, and 5 Na-HEPES (pH 7.4). The standard divalent-free (DVF) Ringer's solution contained (in mM) 150 NaCl, 10 HEDTA, 1 EDTA, and 10 HEPES (pH 7.4 with N-methyl-D-glucamine [NMDG] hydroxide). The 110 mM Ca²⁺ solution contained 110 CaCl₂, 10 D-glucose, and 5 HEPES (pH 7.4). Where indicated, BaCl₂ was substituted for CaCl₂ in this solution. The 12 mM Ba²⁺ solution contained (in mM) 12 mM BaCl₂, 118 NMDG chloride, 10 D-glucose, and 5 HEPES (pH 7.4). For experiments examining block of Na⁺-I_CRAC by Ca²⁺ of the WT Orai1 channels, CaCl₂ was added to the standard DVF solution at the appropriate amount calculated from the MaxChelator software (WEBMAXC 2.10, available at http://www.stanford.edu/cpatton/webmaxc2.htm). For the experiments examining block of Na⁺-I_CRAC by Ca²⁺ in E106D-expressing HEK cells, MgCl₂ was omitted from the standard extracellular solution and CaCl₂ was added to the indicated concentration (i.e., this solution did not contain HEDTA or EDTA). Where indicated, the following organic compounds (purchased from Sigma-Aldrich) were substituted for NaCl in the standard DVF solution: ammonium chloride (NH₄Cl), methylamine HCl (CH₃NH₂-HCl), dimethylamine HCl ((CH₃)₂NH-HCl), trimethylamine HCl ((CH₃)₃N-HCl), and tetramethylammonium chloride ((CH₃)₄NCl), hydroxylamine HCl (NH₂OH-HCl), and hydrazine HCl (NH₂NH₂-HCl). pH was adjusted to 7.4 with NMDG except in the case of hydrazine HCl (pH 6.4) and hydroxylamine HCl (pH 6.2), which were studied at acidic pH to increase the ionized concentration of the test ion. 10 mM TEA-Cl was included in all extracellular solutions to prevent contamination from K⁺_v channels. The standard internal solution contained (in mM) 135 Cs aspartate, 8 mM MgCl₂, 8 BAPTA, and 10 Cs-HEPES (pH 7.2). Where indicated, 10 mM EGTA was substituted for BAPTA in the internal solution.

Stock solutions of thapsigargin (Sigma-Aldrich) and 2-aminoethyldiphenyl borate (2-APB) (Sigma, 20 µM) were prepared in DMSO at concentrations of 1 mM and 20 mM. All solutions were applied using a multibarrel local perfusion pipette with a common delivery port. Reversal potential measurements with 150 mM extracellular KCl applications indicated that the solution exchange time was <1 s.

Patch-Clamp Measurements
Patch-clamp recordings were performed using an Axopatch 200B amplifier (Axon Instruments) interfaced to an ITC-18 input/output board (Instrutech) and an iMac G5 computer. Currents were filtered at 1 kHz with a 4-pole Bessel filter and sampled at 5 kHz. Recording electrodes were pulled from 100-µl pipettes, coated with Sylgard, and fire polished to a final resistance of 2–5 M. Stimulation, data acquisition, and analysis were performed using in-house routines developed on the Igor Pro platform (Wavemetrics). Data are corrected for the liquid junction potential of the pipette solution relative to Ringer's in the bath (–10 mV). The holding potential was +30 mV unless otherwise indicated. Two types of stimuli were usually employed as indicated in the figure legends: (1) a 100-ms step to –100 mV followed by a 100-ms ramp from –100 to +100 mV usually applied every 1 s, or (2) a 300-ms step to –100 mV applied every 1 s. For variance/mean analysis, 200-ms sweeps were acquired at the rate of 4 Hz at a constant holding potential of –100 mV, digitized at 5 kHz, low-pass filtered using a 2 kHz Bessel filter, and recorded directly to hard disk. The mean current and variance were calculated from each sweep.

Data Analysis
The analysis of currents transfected with E190Q Orai1 presented special problems. A large fraction of HEK293 cells transfected with E190Q Orai1 cDNA (typically 30%) exhibited a small I_CRAC-like current in 20 mM Ca²⁺_o (current density 2 pA/pF) and with I-V characteristics similar to native CRAC channels (reversal potential >+70 mV in 20 mM Ca²⁺_o and >+40 mV in DVF solutions). The remaining cells exhibited large currents (typically 10 pA/pF) with a reversal potential of +25 mV in 20 mM Ca²⁺_o and 0 mV in DVF solutions. Because previous reports have demonstrated significant loss of Ca²⁺ selectivity and gain of Cs⁺ permeation by the E190Q Orai1 mutation (Prakriya et al., 2006; Vig et al., 2006a), we assumed that the cells with the small CRAC-like currents had failed to express E190Q Orai1, and manifested predominantly the endogenous CRAC current. Therefore, these cells were excluded from the analysis. For all mutants, only those cells exhibiting current densities 10 pA/pF were included for the analysis to exclude potential problems arising from the contamination of the ectopically expressed Orai1 channels by native CRAC channels of HEK293 cells.

Unless noted otherwise, all data were corrected for leak currents collected in 20 mM Ca²⁺ + 10–100 µM La³⁺. In cells expressing the D110/112/114A Orai1 triple mutant, 100 µM of La³⁺ was employed to account for diminished lanthanide sensitivity of these currents (Yeromin et al., 2006). Averaged results are presented as the mean value ± SEM. All curve fitting was done by least-squares methods using built-in functions in Igor Pro 5.0.

Relative permeabilities were calculated from changes in the reversal potential using the Goldman-Hodgkin-Katz (GHK) voltage equation:

(1)

where R, T, and F have their usual meanings, and P_x and P_Na are the permeabilities of the test ion and Na⁺, respectively, [X] and [Na] are the ionic concentrations, and E_rev is the shift in reversal potential when the test cation is exchanged for Na⁺.

To analyze the voltage dependence of block of Na⁺-I_CRAC by Ca²⁺_o, the model described by Guo and Lu (2000) was employed. This model does not specify the location of the binding site within the field, but instead expresses voltage dependence in terms of an apparent valence, an empirical factor that encompasses the effects of blocker valence and the coupled movements of conducting and blocking ions within the field.

Ca²⁺ binds to a site within the pore according to the reaction:

(1)

where Ch is the channel, Ca²⁺_o and Ca²⁺_i are extracellular and intracellular Ca²⁺, k₁ and k_–1 are the binding and unbinding rates from the extracellular side, k₂ is the unbinding rate from the intracellular side, and each z_i represents the apparent valence for the corresponding transition. We assumed that block from the intracellular compartment is negligible due to nM intracellular Ca²⁺ concentrations. With these assumptions, the fraction of unblocked current is given by Guo and Lu (2000):

(2)

where K₁ = k_–1/k₁ is the equilibrium dissociation constant at 0 applied voltage, and Z₁ and z_i are the apparent valences. Z₁ (= z₁ + z_–1) provides a measure of the overall voltage dependence of Ca²⁺ block and arises from the movement of the charged blocker (Ca²⁺) within the field as well as the possible displacement of permeant ions (Na⁺) within the pore. k₂/k_–1 is the ratio of the rates of Ca²⁺ escaping into the cytoplasm vs. returning to the extracellular solution from the pore, and thus provides a measure of Ca²⁺ permeation. The quantities k₂/k_–1, and z_–1 + z₂ were treated as single adjustable parameters for fitting the data.

Online Supplemental Material
The online supplemental material (Figs. S1–S3, available at http://www.jgp.org/cgi/content/full/jgp.200709872/DC1) contains additional information about the store dependence of activation of the wild-type and mutant CRAC channels, and about the kinetics of fast inactivation in currents arising from WT Orai1 and E190Q Orai1.

	RESULTS

TOP ABSTRACT INTRODUCTION MATERIALS AND METHODS RESULTS DISCUSSION REFERENCES

 
Mutations in the Putative Selectivity Filter Increase Cs⁺ Permeability of CRAC Channels
Recent studies have implicated a set of acidic residues, glutamates at positions 106 and 190 in TM1 and TM3, respectively, and aspartates (D110, D112, and D114) in the TM1–TM2 linker region of human Orai1, as important determinants of the Ca²⁺ selectivity of CRAC channels (Prakriya et al., 2006; Vig et al., 2006a; Yeromin et al., 2006). To uncover additional roles of these residues on ion selectivity, pore geometry, and channel gating, we constructed the E106D, E190Q, and D110/112/114A mutants. The mutations did not change the requirement for Ca²⁺ store depletion to occur before CRAC channel activation (see Fig. S1, available at http://www.jgp.org/cgi/content/full/jgp.200709872/DC1), nor did they affect the spatial redistribution of Orai1 into discrete puncta that follows store depletion (not depicted), suggesting that the alterations do not perturb the overall Orai1 structure. However, as previously shown for some residues, an important unexpected consequence of the mutations is a dramatic increase in Cs⁺ permeability (Prakriya et al., 2006; Vig et al., 2006a; Yeromin et al., 2006). This is illustrated in Fig. 1. HEK293 cells overexpressing either wild-type (WT) or mutant Orai1 together with STIM1 were pretreated with thapsigargin (TG, 1 µM) for 5–10 min before seal formation to deplete Ca²⁺ stores and activate I_CRAC. As previously shown, overexpression of WT Orai1 and STIM1 in HEK293 cells reconstituted a Ca²⁺-selective inward current with key properties similar to native CRAC channels (Mercer et al., 2006; Peinelt et al., 2006). As in native CRAC channels, application of a Na-based divalent free (DVF) solution elicited a prominent Na⁺ current (Fig. 1 A), indicating that Orai1-encoded channels readily conduct Na⁺ in the absence of extracellular divalent ions. However, application of a Cs⁺-based DVF solution elicited only a small monovalent current (Fig. 1 A), consistent with the low Cs⁺ permeability of native CRAC channels (Lepple-Wienhues and Cahalan, 1996; Kozak et al., 2002; Prakriya and Lewis, 2002).

View larger version (17K): [in this window] [in a new window] [Download PPT slide]   Figure 1. Mutation of key Orai1 acidic residues enhances Cs+ permeation through CRAC channels. (A–D) Na+ and Cs+ currents were measured from HEK293 cells overexpressing WT Orai1 or mutant Orai1. Cells were pretreated with 1 µM TG to deplete stores and activate CRAC channels. Leak-corrected currents at –100 mV are plotted against time as the extracellular solution was periodically switched between 20 mM Ca2+ and Na+- or Cs+-based DVF solutions. Application of a Cs+-based DVF solution to cells expressing WT Orai1 fails to elicit significant monovalent CRAC current, whereas cells expressing the mutant Orai1 channels manifest significant Cs+ currents. (E) Plots of the I-V relationship of monovalent currents through CRAC channels in the presence of Na+-based DVF solution. The internal solution was Cs+-aspartate in all experiments.

Orai1 Mutations that Affect Cs⁺ Selectivity Alter CRAC Channel Pore Geometry
What structural change in the pore is responsible for the increased Cs⁺ permeability in the mutants? The limiting dimension of the narrowest region of the pore of the native CRAC channel (3.9 Å) is very close to the atomic diameter of Cs⁺ (3.8 Å), suggesting that steric hindrance could explain the CRAC channel's unusually low permeability to Cs⁺ (Prakriya et al., 2006). We hypothesized that the increased Cs⁺ permeability conferred by the above mutations arises from an increase in the apparent pore diameter of the CRAC channel.

To test this hypothesis, we estimated the narrowest region of the pore of CRAC channels arising from overexpression of wild-type and mutant Orai1. In each case, the minimal dimension of the pore was estimated by examining permeation of a series of organic monovalent cations of increasing size (Prakriya and Lewis, 2006). Ammonium and its methylated derivatives are commonly used to estimate pore size because the progressive addition of methyl groups results in a gradual increase in ion size without introducing gross modifications in the overall ion structure (Dwyer et al., 1980; Burnashev et al., 1996). The diameters of these ions estimated from Corey-Pauling-Koltun space-filling models are as follows (Liu and Adams, 2001): NH₄, 3.2 Å; hydroxylammonium, 3.30 Å; methylammonium, 3.78 Å; dimethylammonium, 4.6 Å; trimethylammonium, 5.34 Å; and tetramethylammonium (TMA), 5.6 Å. Experiments were performed in buffered Ca²⁺-free solutions to avoid the potent blocking effects of Ca²⁺ ions on monovalent CRAC channel currents.

Like native CRAC channels (Prakriya and Lewis, 2006), channels arising from overexpression of WT Orai1 in HEK293 cells failed to appreciably conduct methylated ammonium derivatives (Fig. 2 A). By contrast, large stable currents carried by methylammonium were easily detectable in channels arising from overexpression of E106D, E190Q, and the D110/112/114A triple mutants of Orai1 (Fig. 2, B–D). N-methyl-D-glucamine⁺ was not measurably permeant through any of the Orai1 mutants. The relative permeabilities for the various cations were calculated from changes in reversal potential using the GHK equation (Eq. 1). If we assume that these permeabilities are influenced primarily by steric hindrance rather than ion interactions with the pore, the permeabilities should follow the hydrodynamic relationship (Dwyer et al., 1980; Burnashev et al., 1996)

(3)

A fit of this relationship to the plot of permeability against ion diameter gives an estimated d_pore of 3.8 Å for WT Orai1, close to the value reported previously for native CRAC channels in Jurkat T cells (3.9 Å; Prakriya and Lewis, 2006). By contrast, estimates of d_pore for the Orai1 mutants D110/112/114A Orai1, E106D Orai1, and E190Q Orai1 were 4.4, 5.3, and 7.0 Å, respectively (Fig. 2 E). Thus, Orai1 mutations that diminish the Ca²⁺ selectivity of CRAC channels significantly alter the geometry of the pore. Importantly, the enhancement of the relative Cs⁺ permeability of mutant channels was well correlated with the increase in d_pore (Fig. 2 F), suggesting that the augmentation of Cs⁺ permeation in the mutant channels arises from relief of steric hindrance to Cs⁺ mobility.

View larger version (23K): [in this window] [in a new window] [Download PPT slide]   Figure 2. Orai1 mutations increase the pore diameter of CRAC channels. (A–D) TG-treated HEK293 cells expressing either WT Orai1 or mutant Orai1 were exposed to DVF solutions containing either Na+ or the indicated organic monovalent cation. Methylammonium and larger derivatives of ammonium fail to carry significant current through WT Orai1 channels, but these cations produce significant currents through the mutant channels. (E) The permeabilities of NH4+ and its methylated derivatives relative to Na+ are plotted against the size of each cation. The solid lines are least-squares fit to Eq. 3. The values of dpore estimated from the fits are 3.8 Å (WT Orai1), 4.4 Å, (D110/112/114A Orai1), 5.3 Å (E106D Orai1), and 7.0 Å (E190Q Orai1). (F) Plot of the relative Cs+ permeability (PCs/PNa) against the estimated dpore for WT and mutant Orai1 channels. The straight line is a linear regression with a Pearson's correlation coefficient (R) of 0.91.

View larger version (6K): [in this window] [in a new window] [Download PPT slide]   Figure 3. A rapid phase of current decay in E106D Orai1 channels arising from Ca2+ block of Na+ permeation. (A) HEK293 cells overexpressing E106D Orai1 were treated with thapsigargin to activate ICRAC. Traces show the sequential recordings of current during 300-ms hyperpolarizations to –100 mV from one cell in the presence of 20 mM Ca2+o (with 130 mM Na+) then shortly after Ca2+ concentration was increased to 110 mM Ca2+o (with 0 Na+). A rapid phase of current decay seen in the 20 mM Ca2+o (+130 mM Na+o) solution (indicated by the arrow) is not seen in 110 mM Ca2+ (+0 mM Na+o), suggesting that this transient current arises from Ca2+ block of Na+ permeation. (B) [Ca2+]o dependence of the time constant of Na+-ICRAC decay. The fast decay of currents in 20 mM Ca2+o was fitted by a single exponential function to obtain the time constant of block.

The appearance of time-dependent block in 20 mM Ca²⁺_o and diminished Ca²⁺ selectivity (Prakriya et al., 2006) indicates that the E106D Orai1 substitution alters the energetics of Ca²⁺ binding to the CRAC channel pore. This was examined further by two approaches. First, we measured the voltage dependence of Ca²⁺ blockade of Na-I_CRAC. Voltage-dependent blockade of ion channels generally arises from the combined interaction of blocking molecules with specific pore sites and conducting ions. Thus, alteration of the voltage dependence would provide evidence that E106 regulates Ca²⁺ binding to the pore. Second, we quantified the rates of Ca²⁺ entry and exit from its binding site to determine the role of E106 in controlling the kinetic properties of Ca²⁺ block.

Ca²⁺_o blocked Na⁺-I_CRAC through WT Orai1 channels in a dose-dependent manner with a K_i of 23 µM and a Hill coefficient of 1.0 (n = 4 cells, Fig. 4 A). As shown previously, the E106D substitution increased the K_i of block, in this case, to 490 µM with a Hill coefficient of 1.3 (n = 4 cells, Fig. 4 A). Hyperpolarizing steps revealed that the block was strongly voltage dependent, resulting in a rapid, time-dependent decrease in Na⁺-I_CRAC amplitude (Fig. 4 C). The magnitude of block was estimated as (1-I_ss/I_pk), where I_ss and I_pk are the steady-state and extrapolated peak currents during each voltage step. The results, plotted for five cells in Fig. 4 D, show voltage-dependent blockade that increases with hyperpolarization. In WT Orai1-mediated currents, the block peaked at –100 mV and declined at more negative potentials, a phenomenon also seen in native CRAC channels and likely reflecting escape of Ca²⁺ into the cytoplasm (Prakriya et al., 2006). However, in E106D Orai1 currents, relief of block at negative voltages was not seen.

View larger version (19K): [in this window] [in a new window] [Download PPT slide]   Figure 4. The E106D Orai1 substitution alters the properties of Ca2+ block of Na+-ICRAC. (A) Inhibition of E106D Orai1 Na+ currents by Ca2+o. A TG-treated cell was exposed to various concentrations of [Ca2+]o as indicated by the bars. 20 mM Ca2+o was present in between the bars. The leak-corrected steady-state current measured at the end of 300-ms hyperpolarizing pulses to –100 mV is plotted against time. (B) Concentration dependence of Na+-ICRAC block by Ca2+o in WT and E106D Orai1 currents. The solid line is a least-square fit of the Hill equation, block = 1/[1 + (Ki/[Ca])n] with Ki = 23 µM and n = 1.06 (WT Orai1), and Ki = 490 µM and n = 1.3 (E106D Orai1). Each point shows mean ± SEM from four to five cells from recordings similar to that shown in A. (C) Effects of 300 and 600 µM Ca2+o on Na+ currents during voltage steps. Hyperpolarizing voltage steps from –120 to 0 mV lasting 250 ms were applied from a holding potential of +50 mV in each solution. Note the time-dependent reduction of Na+ currents in the presence of Ca2+o within the first 20 ms following the hyperpolarizing step (indicated by the arrow). (D) Voltage dependence of Ca2+o block. Block was quantified, as [1 – Iss/Ipeak], where Ipeak and Iss are the currents at the start and end of the voltage pulse. Ipeak was determined from a fit of the current decay with a single exponential function and extrapolating back to the start of the voltage pulse. The solid lines are a least-squares fit to Eq. 2. Fit parameters were as follows: WT Orai1 (20 µM): Z1 = 0.78, K1 = 440 µM, z–1 + z2 = 1.7, and k2/k–1 = 0.0007. E106D Orai1([Ca] = 600 µM): Z1 = 1.4, K1 = 820 µM, z–1 + z2 = 1.2, and k2/k–1 = 0.09. E106D Orai1([Ca] = 300 µM): Z1 = 1.27, K1 = 1200 µM, z–1 + z2 = 1.14, and k2/k–1 = 0.09. (E) The time constant of block plotted against [Ca2+]o. Tau was obtained by fitting the current decay in the presence of Ca2+o with a single exponential function. In D and E, each point shows mean ± SEM from four to five cells.

Further, the kinetics of block revealed from experiments in Fig. 4 (D and E) were used to obtain information about the on- and off-rates of Ca²⁺ binding to the block site. At 600 µM [Ca²⁺] the observed time constant of block is 7.2 ms at –100 mV (Fig. 4 E). Assuming that Ca²⁺ accesses a single binding site from the extracellular side at a rate k_on and exits in both directions at a combined rate k_off then it follows that

(4)

The fraction of channels that are blocked is given by

(5)

Substituting the measured values of and the fractional block in Eqs. 4 and 5, we get k_on = 9.8 x 10⁴ M^–1s^–1 and k_off = 94 s^–1 for E106D Orai1-mediated currents. By contrast, k_on and k_off values for currents arising from WT Orai1 were 4.0 x 10⁶ M^–1s^–1 and 116 s^–1. Thus, these results indicate that the E106D Orai1 substitution results in a striking reduction in the k_on for Ca²⁺ binding. Taken together, these results suggest that the E106 regulates high affinity Ca²⁺ binding within the CRAC channel pore, probably by directly interacting with the conducting ions.

Alteration of Orai1 Acidic Residues Diminishes Ca²⁺-dependent Fast Inactivation
Fast inactivation is a prominent hallmark of CRAC channels and occurs from feedback inhibition of channel activity by the local [Ca²⁺]_i around CRAC channels (Hoth and Penner, 1993; Zweifach and Lewis, 1995a; Fierro and Parekh, 1999). The lack of detectable inactivation during hyperpolarizing steps in the above experiments with 110 mM Ca²⁺_o (Fig. 3 A) hinted at an alteration of this process in currents arising from overexpression of E106D Orai1 and prompted us to investigate fast inactivation in greater detail in all the mutants. Fast inactivation was assessed by measuring the extent of current decay during 300-ms hyperpolarizing pulses to –100 mV in 20 and 110 mM Ca²⁺_o solutions. The 110 mM Ca²⁺_o solution confers the advantage that because Ca²⁺ carries all of the current, any confounding effects of Na⁺ permeation that might occur in the 20 mM Ca²⁺_o solution are eliminated. Additionally, the driving force for Ca²⁺ entry, and therefore the local [Ca²⁺] around the cytoplasmic face of the channel, is enhanced in isotonic Ca²⁺_o, thereby accentuating Ca²⁺-mediated fast inactivation.

We have previously shown that currents arising from the overexpression of WT Orai1 in human T cells derived from SCID patients exhibit fast inactivation with properties resembling the inactivation of CRAC currents in Jurkat T cells (Feske et al., 2006). I_CRAC arising from overexpression of WT Orai1 + STIM1 in HEK293 cells also exhibited fast inactivation that was visible both during short (300 ms) and long (2 s) hyperpolarizing voltage steps (Fig. 5 A and Fig. S2). In the D110/112/114A Orai1 triple mutant, there was a small but significant reduction in the extent of fast inactivation in 20 and 110 mM [Ca²⁺]_o solutions. In currents arising from E106D Orai1, fast inactivation could only be assessed in 110 mM Ca²⁺_o as currents in 20 mM Ca²⁺_o + 130 mM Na⁺ resulted in rapid block of Na⁺ permeation (Fig. 3). In these mutant channels, fast inactivation in 110 mM Ca²⁺_o was virtually absent. The most striking effect was in currents arising from the E190Q Orai1 substitution. Instead of inactivation, hyperpolarizing steps in 20 mM Ca²⁺_o triggered a gradual, time-dependent enhancement of the current (Fig. 5 D) in these channels. Examination of currents in 110 mM Ca²⁺_o revealed an additional degree of complexity. In this condition, an initial decay of the current lasting 40 ms was followed by robust enhancement resulting in a significantly larger current at the end of the hyperpolarizing step (Fig. 5 D). Switching the external solution to a DVF solution eliminated the time-dependent enhancement of the current (Fig. 5 D), indicating that the facilitation during hyperpolarization is Ca²⁺ dependent. Additionally, increasing [Ca²⁺]_o from 2 to 110 mM progressively increased the degree of enhancement (Fig. S3). Thus, it appears that the E190Q substitution alters Ca²⁺ regulation of the channel such that feedback inactivation during hyperpolarizing pulses is overcome by strong Ca²⁺-mediated facilitation.

View larger version (20K): [in this window] [in a new window] [Download PPT slide]   Figure 5. Alteration of Orai1 acidic residues diminishes fast inactivation. (A–D) Representative traces showing fast inactivation in WT and mutant Orai1 channels during hyperpolarizing steps to –100 mV applied every second. Note that the current measured in 20 mM Ca2+ from E106D channels shows Ca2+ block of Na+-ICRAC (Fig. 3). Currents measured through E190Q channels exhibit complex kinetics with a small initial decay followed by strong enhancement during the voltage step. (E and F) Alteration of fast inactivation in the mutants in the presence of 20 and 110 mM [Ca2+]o. Inactivation of ICRAC was quantified here as fraction of the remaining current, Iss/Ipeak, where Ipeak and Iss are the currents at the start and end of the 300-ms voltage pulses. Inactivation of current arising from E106D Orai1 channels could not be assessed in 20 mM Ca2+o due to rapid block of Na+ permeation (Fig. 3). * indicates significant (P < 0.05) difference between mutant and WT currents. Each point shows mean ± SEM from at least five cells.

View larger version (7K): [in this window] [in a new window] [Download PPT slide]   Figure 6. Alteration of inactivation is correlated with changes in ion selectivity and pore diameter. (A and B) Inactivation of ICRAC measured in 110 Ca2+ is plotted against the reversal potential of currents measured in either 20 mM Ca2+ (A) or in DVF solution (B) for WT and mutant Orai1 channels. Inactivation of ICRAC was quantified as Iss/Ipeak, where Ipeak and Iss are the currents at the start and end of a 300-ms voltage pulse to –100 mV, respectively. The Pearson's correlation coefficients (R) for plots in A and B are 0.81 and 0.92, respectively. (C) Inactivation ICRAC measured in 110 mM Ca2+o is plotted against the measured pore diameter determined from Fig. 2 E. The Pearson's correlation coefficient (R) for this plot is 0.99.

Diffusion models predict that the local [Ca²⁺]_i a short distance away from the Ca²⁺ channel pore is directly related to the unitary Ca²⁺ current (Neher, 1986). Thus, any reduction in the unitary conductance would strongly diminish the local [Ca²⁺] in the vicinity of the inactivation binding site and reduce fast inactivation. To test for alteration of unitary Ca²⁺ current, we used nonstationary noise analysis to compare the unitary conductance of CRAC channels arising from WT Orai1 and E106D Orai1. The mean (I) and variance (²) of the macroscopic currents were calculated for a series of 200-ms epochs in TG-treated cells in 110 mM Ca²⁺. As shown previously for endogenous CRAC channels in Jurkat T cells (Zweifach and Lewis, 1993), a switch of the extracellular solution from 20 to 110 mM Ca²⁺ causes I_CRAC to increase rapidly and then decay slowly (Fig. 7 A), possibly from slow Ca²⁺-dependent inactivation (Zweifach and Lewis, 1995b). Plots of the current variance against the mean current amplitude obtained during the slow decay of the current were well fitted by straight lines (Fig. 7). The average slopes (²/I) in WT Orai1 and E106D Orai1 currents were –6.6 ± 1.6 (n = 5) and –6.5 ± 1.6 fA (n = 4), respectively. Furthermore, application of a low concentration of 2-APB (5 µM), which is known to potentiate native CRAC channels (Prakriya and Lewis, 2001), enhanced both WT Orai1 and E106D Orai1 currents with similar ²/I slopes (Fig. 7). If the open probability of CRAC channels in 110 mM Ca²⁺ during the current decay is low (P_o << 0.5), the slope (²/I) provides a direct estimate of the unitary CRAC channel current (Prakriya and Lewis, 2006). If the P_o were higher, as might be expected for the E106D Orai1 channels due to lack of inactivation, the true unitary conductance would, in fact, be larger by a factor of 1/(1 – P_o) (Prakriya and Lewis, 2006). Thus, the close similarity in the unitary current argues that, at least under these conditions, the local [Ca²⁺] in the vicinity of the inactivation binding site should be comparable in E106D Orai1 and WT Orai1 channels.

View larger version (20K): [in this window] [in a new window] [Download PPT slide]   Figure 7. Noise analysis indicates that the unitary Ca2+ conductance of E106D Orai1-mediated CRAC channels is not diminished. (A) HEK293 cells overexpressing either WT Orai1 or E106D Orai1 were pretreated with TG and held at a constant potential of –100 mV to record ICRAC. Application of 110 mM [Ca2+]o causes a slow decline of ICRAC, likely due to slow inactivation (Zweifach and Lewis, 1995b). After the current reached a steady amplitude in 110 mM Ca2+, a low concentration (5 µM) of 2-APB was applied to enhance ICRAC. The average and the variance of the Ca2+ current was measured from 200-ms sweeps collected every 0.25 s. (B and C) Mean-variance analysis of ICRAC measured during the slow decline in 110 mM Ca2+o and enhancement of ICRAC by 2-APB from the experiment shown in A. The data are well fit by straight lines with slopes of 7 fA (B) and 5 fA (C). (D) As with WT Orai1 currents, E106D Orai1 currents exhibit slow inactivation in the presence of 110 mM Ca2+o. Application of a low concentration (5 µM) of 2-APB enhances E106D Orai1 currents. (E and F) Mean-variance analysis of ICRAC measured during the slow decline in 110 mM Ca2+o and 2-APB potentiation from the experiment shown in D. Plots were fit with straight lines with slopes of 4 fA (E) and 7 fA (F).

View larger version (15K): [in this window] [in a new window] [Download PPT slide]   Figure 8. Fast inactivation of Ba2+ currents is also diminished by the E106D Orai1 mutation. (A) Ba2+ currents through WT Orai1 channels exhibit fast inactivation. A TG-pretreated cell was exposed to 110 mM Ca2+ and 110 mM Ba2+ solutions as indicated by the bars. Exposure to 110 mM Ba2+ results in a transient Ba2+ current that declines within 20 s, presumably due to CRAC channel depotentiation. The lower traces show selected current responses to hyperpolarizing voltage steps (to –100 mV, 300 ms) at the time points indicated by arrowheads in the upper plots. (B) Ba2+ currents through E106D Orai1 channels. Exposure to 110 mM Ba2+ results in a large Ba2+ current that gradually declines over tens of seconds. The lower traces show inactivation of selected currents at the time points indicated by the arrowheads in the upper plots. Note the presence of inactivation in the 110 mM Ba2+ trace. (C) In a different cell, when the external Ba2+ concentration was reduced to 12 mM to produce comparable Ca2+ and Ba2+ current amplitudes, no difference in fast inactivation was observed between the Ca2+ and Ba2+ currents. To prevent monovalent permeation, NMDG+ was used instead of Na+ in the 12 mM Ba2+ solution (see Materials and methods). (D) The extent of inactivation (1 – Iss/Ipeak) in currents arising from WT Orai1 and E106D Orai1 channels in Ca2+ and Ba2+.

Intracellular dialysis with EGTA resulted in only a minor reappearance of fast inactivation in E106D Orai1 currents (Fig. 9 A). Hyperpolarization, which would be expected to increase the size of the unitary Ca²⁺ current, sped up the kinetics of current decay and increased the extent of inactivation to a small degree, consistent with the known Ca²⁺ dependence of the process (Fig. 9 C). At a membrane potential of –100 mV, the degree of inactivation (1– I_ss/I_peak) with intracellular EGTA in 110 mM Ca²⁺_o in E106D Orai1 currents was only 34 ± 3% compared with 72 ± 5% in WT Orai1-mediated currents. These results indicate that elevating the local [Ca²⁺] in the vicinity of the CRAC channels fails to restore inactivation. Together with the lack of change of the unitary Ca²⁺ conductance noted in the preceding sections, these results suggest that the loss of inactivation elicited by the Orai1 mutations arises from alteration of the inactivation mechanism itself.

View larger version (10K): [in this window] [in a new window] [Download PPT slide]   Figure 9. Intracellular dialysis of EGTA does not restore fast inactivation in E106D Orai1-expressing cells. (A) ICRAC measured from cells expressing E106D Orai1 during hyperpolarizing pulses to –100 mV in two different cells dialyzed with either BAPTA (8 mM) or EGTA (10 mM). Replacing intracellular BAPTA with EGTA restores fast inactivation only to a small degree. (B) In cells expressing WT Orai1, replacing intracellular BAPTA with EGTA strongly enhances fast inactivation. In A and B, cells with comparable peak currents were chosen to compare inactivation rates of cells dialyzed with EGTA or BAPTA. (C) Responses of one cell to a graded series of hyperpolarizing steps to the indicated potential from a holding potential of +30 mV with 10 mM EGTA in pipette.

	DISCUSSION

TOP ABSTRACT INTRODUCTION MATERIALS AND METHODS RESULTS DISCUSSION REFERENCES

Relationship between Cs⁺ Permeability and Pore Geometry
A key hallmark of CRAC channels is their unusually low permeability to Cs⁺ (P_Cs/P_Na 0.1; Lepple-Wienhues and Cahalan, 1996; Kozak et al., 2002; Prakriya and Lewis, 2002). This is in contrast to Ca²⁺-selective L-type Ca_v channels, which readily pass Cs⁺ under DVF conditions (P_Cs/P_Na 0.6; Hess et al., 1986). Similarly, the Ca²⁺-selective TRPV6 channel is also highly permeable to Cs⁺ (P_Cs/P_Na of 0.5; Voets et al., 2001). A recent study provides a clue for the basis for the CRAC channel's unusually low Cs⁺ permeability. Whereas the L-type Ca_v and TRPV6 channels have relatively large pore diameters of 6.2 Å (Cataldi et al., 2002) and 5.4 Å (Voets et al., 2004), respectively, the narrowest region of the CRAC channel pore is only 3.9 Å (Prakriya and Lewis, 2006). Given that the atomic diameter of a naked Cs⁺ ion is 3.8 Å, this raised the possibility that the low permeability to Cs⁺ reflects steric hindrance to its permeation.

In the present study, we show that mutations of acidic residues that diminish Ca²⁺ selectivity and enhance Cs⁺ permeation cause significant widening of the pore. In particular the E106D and E190Q substitutions, which strongly elevate Cs⁺ permeation, increase the apparent pore diameter from 3.8 Å (WT Orai1) to 5.3 and 7.0 Å, respectively. It is important to note that these estimates reflect the dimensions of the pore when conducting monovalent ions and may differ under conditions of Ca²⁺ permeation. This is because as Ca²⁺ diffuses into the pore, its binding to the flexible carboxylate side chains of the Glu and Asp residues may provide the necessary countercharges to stabilize and alter the conformation of the selectivity filter in a manner analogous to the change in the conformation of the KcsA channel selectivity filter by the entry of K⁺ (Zhou and MacKinnon, 2003). Nevertheless, the profound alterations of Cs⁺ permeability and pore diameter by the mutations suggests that regardless of the precise side chain orientation of the carboxylates, the acidic residues examined here regulate ion selectivity by influencing the structural and architectural features of the CRAC channel pore.

At present, we do not know the structural basis of the substantial alterations in pore geometry caused by alterations of the acidic residues. If we assume that these residues are all arranged at a single locus within the selectivity filter, analogous to the EEEE locus of L- type Ca_v channels, then one possibility is that the location of the minimal pore diameter is at the position of the acidic residues themselves. In this paradigm, the GluAsp substitution widens the selectivity filter due to shortening of the side chains. Similarly, elimination of the Asp side chains could underlie the widening seen with the D110/112/114A substitutions. However, this argument does not readily explain the widening seen with the E190Q substitution, which does not alter the chain length.

A second possibility is that the acidic residues are not arranged at a single locus, but instead mediate multiple binding sites along the permeation pathway and contribute to differing extents to the apparent pore diameter. It has been proposed that the CRAC channel permeation pathway might contain multiple binding sites with the TM glutamate residues (E106 and E190) contributing sites within the pore and the D110/112 residues forming a site at the extracellular pore mouth (Vig et al., 2006a; Yeromin et al., 2006). If true, the large alteration in the apparent pore diameter seen with the E190Q substitution might suggest that this residue contributes most significantly to the apparent CRAC channel pore diameter.

A further explanation is that the narrowest part of the pore occurs not at the locus of the acidic residues themselves, but rather at a neighboring site in the ion conduction pathway, with the mutations triggering changes in geometry at this distal location. Although simplistic and speculative, this idea embodies the notion that there may be a barrier for ion permeation in series with the selectivity filter responsible for sculpting the unique hallmarks of the CRAC channel pore such as its low Cs⁺ permeability. Regardless of the exact mechanism, the increase in pore diameter and accompanying changes in key Ca²⁺-binding properties such as the k_on and the voltage dependence of Ca²⁺ block collectively indicate that the mutations trigger complex changes in the tertiary structure of the selectivity filter with profound consequences for ion permeation, selectivity, and channel gating.

Ca²⁺ Block of Na⁺-I_CRAC
Characterization of the kinetics and extent of Ca²⁺ block of Na⁺ currents indicates that the association rate of Ca²⁺ block of Na⁺ conduction through WT Orai1 channels is 4 x 10⁶ M^–1s^–1. This is similar to the value found in the endogenous CRAC channels of Jurkat T cells (6 x 10⁶ M^–1s^–1; Prakriya and Lewis, 2006). The measured pore diameter of 3.8 Å for channels arising from WT Orai1 is also very similar to that determined for the endogenous CRAC channels of Jurkat cells (d_pore of 3.9 Å; Prakriya and Lewis, 2006). Together with similarities in the affinity of Ca²⁺ block of Na⁺-I_CRAC (20 µM; Prakriya and Lewis, 2006; Prakriya et al., 2006) and of the voltage dependence of Ca²⁺ block (apparent valence of 0.78 in WT Orai1 channels vs. 0.71 in Jurkat T cell CRAC channels; Prakriya and Lewis, 2006), these findings indicate that the Ca²⁺ binding properties of the native CRAC