Role of a conserved glutamine in the function of voltage-gated Ca2+ channels revealed by a mutation in human CACNA1D

Voltage-gated Cav Ca2+ channels play crucial roles in regulating gene transcription, neuronal excitability, and synaptic transmission. Natural or pathological variations in Cav channels have yielded rich insights into the molecular determinants controlling channel function. Here, we report the consequences of a natural, putatively disease-associated mutation in the CACNA1D gene encoding the pore-forming Cav1.3 α1 subunit. The mutation causes a substitution of a glutamine residue that is highly conserved in the extracellular S1–S2 loop of domain II in all Cav channels with a histidine and was identified by whole-exome sequencing of an individual with moderate hearing impairment, developmental delay, and epilepsy. When introduced into the rat Cav1.3 cDNA, Q558H significantly decreased the density of Ca2+ currents in transfected HEK293T cells. Gating current analyses and cell-surface biotinylation experiments suggested that the smaller current amplitudes caused by Q558H were because of decreased numbers of functional Cav1.3 channels at the cell surface. The substitution also produced more sustained Ca2+ currents by weakening voltage-dependent inactivation. When inserted into the corresponding locus of Cav2.1, the substitution had similar effects as in Cav1.3. However, the substitution introduced in Cav3.1 reduced current density, but had no effects on voltage-dependent inactivation. Our results reveal a critical extracellular determinant of current density for all Cav family members and of voltage-dependent inactivation of Cav1.3 and Cav2.1 channels.

Within the human genome, hundreds of genes encode ion channels, mutations which cause diseases affecting a variety of organ systems (1). In many cases, electrophysiological analysis of the mutations has uncovered mechanistic insights into channel function. Mutations causing long QT syndrome have revealed molecular determinants controlling voltage-depen-dent activation (2) and inactivation (3). In addition, the functional characterization of a mutation causing congenital stationary night blindness in the CACNA1F gene encoding voltage-gated Ca v 1.4 Ca 2ϩ channels uncovered a role for a C-terminal automodulatory domain that suppresses Ca 2ϩ -dependent inactivation of these channels (4). Thus, the study of channelopathy-causing mutations can advance understanding of the structure-function relationships of ion channel proteins as well as inform about the physiological roles of these channels in humans.
In addition to congenital stationary night blindness, mutations affecting voltage-gated Ca 2ϩ channels cause migraine, ataxia, malignant hypothermia, and cardiac arrhythmia (5,6). Within the Ca v channel 2 superfamily, Ca v 1.3 L-type channels have emerged as targets of heterogeneous disorders including deafness and autism (7,8). In the ear, Ca v 1.3 channels are localized at inner hair cell synapses where they mediate Ca 2ϩ influx required for transmission of auditory information into the brain via the auditory nerve (9,10). A mutation involving a glycine insertion in the domain IS6 transmembrane helix prevents voltage-dependent gating of Ca v 1.3 channels and causes deafness in homozygous individuals carrying this mutation (11). Ca v 1.3 channels are also expressed in the brain, where they regulate mood-related behaviors (12) and fear consolidation (13). De novo mutations causing a gain of function in Ca v 1.3 are associated with intellectual disability and autism spectrum disorders (7,8).
Here, we report the functional effects of a mutation in the CACNA1D gene encoding Ca v 1.3 found by exome sequencing of an individual with intellectual disability, developmental delay, and mild hearing impairment. We show that this mutation, which causes a histidine substitution for a glutamine residue in the extracellular loop connecting the S1 and S2 helices in domain II (IIS1-S2L), causes a significant reduction in Ca v 1.3 current density and yet suppresses voltage-dependent inactivation of these channels in transfected cells. The glutamine residue is conserved in all Ca v family members and its mutation to histidine has similar effects when introduced into Ca v 2.1 and Ca v 3.1 channels. Our findings implicate IIS1-S2L as a key determinant for regulating the cell-surface density and inactivation gating of Ca v channels.

Q567H mutation in an individual with moderate hearing impairment and intellectual disability
The male proband was the first-born child of healthy parents who were first cousins of Arabic descent. The father had a history of febrile seizures since childhood. The proband was born after an uneventful pregnancy and delivered at 38 weeks of gestation by cesarean section because of fetal distress. He failed the hearing test at birth and had seizures beginning at 4 months of age. At the age of 4 1 ⁄ 2 years, he was reassessed at the genetic clinic for the purpose of family planning. He did not speak and his comprehension was limited. He exhibited autistic behaviors including very limited eye contact and hyperactivity. He wore hearing aids from the age of 11 months. At 5 years of age, he underwent hearing tests which revealed moderate hearing loss (40 -50 db SPL across all frequencies from 500 Hz to 40 kHz).
Whole-exome sequencing and filtering for homozygous variants revealed two potential candidate genes: A homozygous mutation in the otogelin gene (OTOG): Chr11(GRCh37): g.17631250CϾT; NM_001277269.1: c.4439CϾT (p.Thr1480Ile); and a homozygous missense variant in CACNA1D: Chr3 (GRCh37): g.53756476GϾC; NM_000720.3: c.1701GϾC (p.Gln567His). The OTOG variant was present in a heterozygous state in the gnomAD browser (0.00003348), whereas the CACNA1D was found neither in gnomAD nor in an in-house exome database. Both parents were heterozygous for the OTOG and CACNA1D variants. The high degree of consanguinity within the family was the basis for assuming a recessive model of inheritance. A full pedigree molecular analysis was not possible because of difficulties in obtaining consent from other family members.
OTOG encodes otogelin, an N-glycosylated protein that is enriched in cellular membranes of the inner ear (14). Mutations in OTOG cause autosomal-recessive moderate hearing impairment (15). CACNA1D encodes the pore-forming subunit of Ca v 1.3 Ca 2ϩ channels (Fig. 1A), which mediate the presynaptic release of glutamate from inner hair cells in the cochlea (9,16). A mutation in CACNA1D that causes a complete loss of Ca v 1.3 function causes congenital deafness and sinoatrial node dysfunction (11). CACNA1D variants have also been associated with autism spectrum disorders (8) and epilepsy (8,17,18). Because of the coincident involvement of CACNA1D in the proband's symptoms (i.e. hearing loss, epilepsy, autism), we analyzed the functional consequences of the CACNA1D in a heterologous expression system.

The Gln to His mutation decreases current density and cell-surface levels of Ca v 1.3
The mutation in the affected individual causes substitution of a histidine for a glutamine residue in the proximal third of the extracellular loop connecting transmembrane helices S1 and S2 in domain II (Q567H) (Fig. 1B). To elucidate the functional consequences of Q567H, we introduced this mutation in the corresponding region of the rat Ca v 1.3 channel (Q558H). We used whole-cell patch clamp electrophysiology to compare the properties of the wild-type (WT) and mutant channels coexpressed with the auxiliary ␤ 2a and ␣ 2 ␦ 1 subunits in human embryonic kidney cells transformed with T-antigen (HEK293T). In current-voltage (I-V) relationships, the most prominent difference between cells transfected with WT and Q558H was a significant reduction in the peak current density of the mutant channel as compared with the WT channel ( Fig. 2A; Table 1). To determine whether this was because of alterations in the number of functional Q558H channels and/or their open probability, we analyzed gating currents evoked by depolarization to the reversal potential (I gating ) and tail currents (I tail ) upon repolarization to the holding voltage (Fig. 2B). As was described previously (19,20), the time integral of I gating (Q max ) is a measure of the number of channels in the cell membrane, and the ratio of I tail and Q max describes the relative channel open probability. Q558H caused a significant reduction in Q max (56.5 Ϯ 6.4 fC for WT versus 16 Ϯ 2.1 fC for Q558H; p Ͻ 0.01 by t test), but no change in the slope of the I tail versus Q max relationship (9.8 Ϯ 0.8 for WT versus 12.4 Ϯ 1.2 for Q558H; p ϭ 0.62 by F-test) (Fig. 2C). Thus, Q558H causes a decrease in the number of functional channels, but not their open probability.
Q558H could cause greater degradation of the channel because of improper folding. Alternatively, the mutation could alter the trafficking of the channel to or from the plasma membrane. To distinguish between these possibilities, we compared the amount of plasma membrane-associated WT and mutant channel protein in cell-surface biotinylation assays. Cell-surface proteins were subject to biotinylation and pulldown by streptavidin beads, and the levels of WT and mutant channels were compared by Western blotting with Ca v 1.3 antibodies.

A conserved glutamine in voltage-gated Ca 2؉ channels
Blotting with antibodies of Na ϩ /K ϩ ATPase was used as a positive control for the cell-surface biotinylation reaction. Expression of the WT or mutant channels had no effect on the amount of cell-surface Na ϩ /K ϩ ATPase as the corresponding signals measured by densitometry were not significantly different in untransfected control cells (0.28 Ϯ 0.1 arbitrary units (a.u.)) or those transfected with WT (0.33 Ϯ 0.11 a.u.) or Q558H (0.24 Ϯ 0.07 a.u.; p ϭ 0.83 by analysis of variance (ANOVA)).
Although there was no difference in the total levels of WT and mutant channel protein in whole-cell lysates, there was a significant reduction in the levels of Q558H channel protein at the cell surface compared with WT channels (Fig. 3, A and B). To confirm these results, we also compared the cell-surface localization of WT and Q558H mutant channels containing an extracellular FLAG epitope between S5 and S6 of domain II. Consistent with the lack of effect of epitopes in this region on channel function (21), these FLAG-tagged channels exhibited similar electrophysiological properties as the untagged channels (data not shown). Cells transfected with FLAG-tagged WT and mutant channels were labeled first with FLAG antibodies under nonpermeabilizing conditions followed by permeabilization and labeling with Ca v 1.3 antibodies to mark the distribution of cell-surface and total channels, respectively. Because of  Table 1. Results are from four independent experiments. B, representative traces showing I gating evoked by a 20-ms step depolarization from Ϫ80 mV to ϩ70 mV and I tail upon repolarization to Ϫ80 mV. Insets show I gating and I tail on expanded time scale. C, left, average maximal gating charge (Q max ) for WT and mutant; p value was determined by t test. Right, I tail for each cell is plotted against Q max . The fit by linear regression and coefficient of determination (r 2 ) are shown for WT (n ϭ 19 cells) and Q558H (n ϭ 27 cells). Results are from four independent experiments. Error bars represent S.E.

Table 1 Parameters of I-V relationships
Values for V1 ⁄ 2 and k were obtained from Boltzmann fits of data in Fig. 2A. Results were obtained from six independent experiments, and p values were determined by t test. pA/pF, picoamperes/picofarads.

Channel
WT n Mutant n pvalue Ϫ14.5 Ϯ 2.8 0.02 Ϫ30.1 Ϯ 7.8 0.03 Ϫ67.9 Ϯ 13 Ͻ0.01 our transient transfection strategy, channel expression levels varied considerably between cells as detected by Ca v 1.3 immunofluorescence, and this impacted the intensity of cell-surface FLAG-labeled channels. However, in cells selected for a similar level of total Ca v 1.3 labeling, FLAG labeling was less apparent at the plasma membrane when transfected with Q558H than with the WT channel (Fig. 3C). Taken together, our results show that Q558H causes a loss of function because of improper cell-surface trafficking of Ca v 1.3.

Q558H alters voltage-dependent inactivation of Ca v 1.3
Compared with WT currents, Ca 2ϩ currents (I Ca ) mediated by Q558H channels were more prolonged ( Fig. 2A), suggesting defects in inactivation kinetics. To characterize this more thoroughly, we performed exponential curve fitting of I Ca evoked by sustained (3-s) pulses to ϩ10 mV from a holding voltage of Ϫ80 mV. Traces corresponding to WT and mutant I Ca were best fit with a double-exponential function, indicating the presence of a fast and slow phase of inactivation. Compared with WT channels, the time constant ( fast ) and fractional amplitude of the fast phase of inactivation (F fast ) were significantly reduced in mutant channels (Fig. 4, A-C). Conversely, the time constant ( slow ) and fractional amplitude of slow inactivation (F slow ) for the mutant channels were significantly greater than for WT channels (Fig. 4, A-C). These results indicate that Q558H prolongs Ca v 1.3 currents by stabilizing slow inactivation.
To investigate effects of Q558H on the voltage-dependence of inactivation, we measured the impact of a 5-s conditioning prepulse to various voltages on inactivation of a test current evoked by a 40-ms pulse to ϩ10 mV. The amplitude of the test current was normalized to that with a Ϫ60 mV prepulse (Normalized I Ca ) and plotted against prepulse voltage. Note that although we used Ca 2ϩ as the charge carrier, this voltage protocol measures voltage-dependent inactivation (VDI) rather than Ca 2ϩ -dependent inactivation given the relatively long duration of the conditioning prepulse (22). Consistent with a voltage-dependent process, Boltzmann fits of the data showed a sigmoidal dependence on prepulse voltage. Parameters from these fits indicated no difference in the voltage of half-maximal inactivation (V1 ⁄ 2 ) but a significantly greater slope (k) for Q558H compared with WT channels (Fig. 5A; Table 2). To determine whether this effect of Q558H might be secondary to effects on voltage-dependent activation, we analyzed the amplitude of tail currents measured by 2-ms repolarizations to Ϫ60 mV from various test voltages. The amplitudes of tail currents were normalized to that using a ϩ40 mV test pulse and plotted against test voltage. Boltzmann fits of the data yielded similar parameters for WT and Q558H channels in terms of V h and k ( Fig. 5B; Table 3), which argued against an effect of Q558H on voltagedependent activation. Overlay of the inactivation and activation curves revealed a sizeable increase in window currents generated by Q558H (Fig. 5C). Thus, despite the decrease in peak current density caused by Q558H, the mutation also has a gain of function effect in allowing Ca 2ϩ entry at voltages at which WT channels are inactivated.

Alterations in Ca v 2.1 and Ca v 3.1 properties by similar Gln to His mutations
The strong conservation of Gln-567 (Fig. 1B) suggests that it may play a key role in regulating current density and steadystate inactivation of all Ca v channel family members. To test this, we introduced the Gln to His mutation in the corresponding sites of Ca v 2.1 (Q511H) and the more distantly related Ca v 3.1 (Q767H). Similar to the effects of Q558H in Ca v 1.3 (Fig.  2), Q511H and Q767H caused significant reductions in peak current densities but without a change in Boltzmann parameters for the I-V relationship (Fig. 6, A and B; Table 1). We were unable to test whether the reduced peak current density for Q511H and Q767H was due to decreased cell-surface channel protein, as was done for Ca v 1.3 and Q567H (Fig. 3, A and B) because of low levels of cell-surface biotinylation of the mutant channels. However, gating current analyses revealed a signifi-   (Fig. 7). Thus, the Gln to His mutations in all three classes of Ca v channel decrease the number of functional channels, but not their open probability, which could be explained by decreased levels of channels in the plasma membrane.
Steady-state inactivation was also suppressed by Q511H as compared with WT Ca v 2.1, but was more severe than Q558H in Ca v 1.3 in that there was a strong positive shift in V1 ⁄ 2 because of Q511H ( Fig. 8A; Table 2). Although voltage-dependent activation of Q511H was not affected (Fig. 8B; Table 3), the window current generated by the mutant channels was greater than for WT channels because of lack of inactivation at positive voltages (Fig. 8C). In contrast to the robust effects of Q558H and Q511H on inactivation in Ca v 1.3 and Ca v 2.1, respectively, VDI of Ca v 3.1 was not affected by Q767H (Fig. 9A; Table 2). Although A, left, representative traces of currents evoked by 5-s prepulses to various voltages prior to a 40-ms test pulse to ϩ10 mV from a holding voltage of Ϫ80 mV. Right, test current amplitudes were normalized with the test current amplitude obtained to Ϫ60 mV and plotted against prepulse voltage for WT (n ϭ 7 cells) and Q558H mutant channels (n ϭ 7 cells). B, left, representative traces of currents evoked by 10-ms steps to various voltages prior to 2-ms repolarization to Ϫ60 mV to better resolve tail currents. The amplitudes of tail currents were normalized to those obtained with ϩ40 mV test pulse and plotted against test voltage for WT (n ϭ 9 cells) and Q558H mutant channels (n ϭ 10 cells). C, left, overlay of activation and inactivation curves for WT and mutant channels from data in A and B. Right, data from C are shown on expanded vertical scale. Window current of Q558H (gray shading) extends to depolarized voltages at which WT window current (dark shading) is absent. Results are from six independent experiments. Error bars represent S.E.  (n ϭ 14 cells, B). Voltage protocol was the same as in Fig. 2A. Results are from four independent experiments. Error bars represent S.E.

Table 2 Parameters for voltage-dependent inactivation
Values for V1 ⁄ 2 and k were obtained from Boltzmann fits of data in Fig. 5A. Results were obtained from nine independent experiments, and p values were determined by t test.

A conserved glutamine in voltage-gated Ca 2؉ channels
Q767H modestly suppressed voltage-dependent activation of Ca v 3.1 ( Fig. 9B; Table 3), there was limited impact of this mutation on the window current (Fig. 9C). We conclude that Gln-558 and Gln-511 are determinants regulating VDI of high-voltage activated (i.e. Ca v 1 and Ca v 2) Ca 2ϩ channels, and of the cell-surface density of all Ca v channel family members.

Discussion
VDI is a crucial feature of most Ca v channels, limiting Ca 2ϩ influx that could otherwise cause aberrant cellular excitability (23,24). A variety of factors regulate VDI of Ca v channels, including the auxiliary Ca v ␤ subunit (reviewed in Ref. 25). Within the main ␣ 1 subunit, multiple elements contributing to VDI of Ca v 1 and Ca v 2 channels have been identified in cytoplasmic regions (26 -28) or transmembrane segments (29,30). Considering that Gln-558 and Gln-511 are located in the extracellular S1-S2 loop of domain II (IIS1-S2L) (Fig. 1), how these residues might promote VDI (Figs. 5 and 7) is unclear. Together with S3 and S4, S1 and S2 form the voltage sensor domain (VSD), a conserved module in voltage-gated ion channels that senses changes in the cell membrane potential in ways that alter the opening, closing, and inactivation of the channel (31). Within Ca v 1 and Ca v 2 channels, VSDs in domains I through IV may be functionally heterogeneous, with domains II and III playing the most prominent roles in voltage-dependent activation (32) and VDI (33). Thus, Gln-558 and Gln-511 may modulate the contribution of the VSD in domain II (VSDII) to gating of VDI. In the structure of the skeletal muscle Ca v 1.1 complex by cryo-EM, S1-S2 loops are closely associated with the porelining S5-S6 loop (34). The substitution of a charged histidine residue at this site might produce electrostatic alterations that destabilize VDI by maintaining the pore in an open rather than inactivated conformation (Figs. 5 and 7). Determinants for VDI

A conserved glutamine in voltage-gated Ca 2؉ channels
in Ca v 3 channels differ from those in Ca v 1 and Ca v 2 channels (35), which may explain the lack of effect of the Gln to His mutation on Ca v 3.1 VDI (Fig. 9).
At first glance, the inhibitory effects of the Gln to His mutations on current density in all three classes of Ca v channels (Figs. 2 and 6) might simply result from a decrease in the stability of the channel protein. However, in the context of Ca v 1.3, Q558H diminished the amount of Ca v 1.3 protein on the cell surface but had no effect on total protein Ca v 1.3 levels (Fig. 3). Moreover, gating current analyses indicate a decrease in the number of functional Ca v 1.3 channels with the Q558H mutation (Fig. 2, B and C). Taken together, these results suggest that the decrease in current density caused by the Gln to His mutation results from an impairment in the trafficking and/or maintenance of channels in the plasma membrane. Although auxiliary ␣ 2 ␦ subunit promotes the cell-surface trafficking of these channels through interactions with extracellular sites in the ␣ 1 subunit (reviewed in Ref. 36), the Gln to His mutations are unlikely to disrupt the interaction of Ca v 1.3 and Ca v 2.1 with ␣ 2 ␦ given that domain II S1-S2 is distant from sites in Ca v ␣ 1 that are found in contact with ␣ 2 ␦ based on the cryo-EM structure of Ca v 1.1 (34). Moreover, Q767H also significantly reduces the current density of Ca v 3.1 expressed without ␣ 2 ␦ (Fig. 6B).
Within K v channels, a conserved threonine in the S1-S2 linker (S1-S2L) is required for efficient trafficking of the channel to the plasma membrane (37,38). This residue is thought to interact with a glutamate in the S5-S6L (39), which may stabilize VSDs and gating domains in ways that support K v channel trafficking out of the endoplasmic reticulum and into the plasma membrane (37). In the cryo-EM structure of Ca v 1.1, the S1-S2L in all the VSDs is well resolved, indicative of rigidity that may be important for its interaction with LS5-S6 (34). Similar to the role of S1-S2L in K v channels (37,38), efficient trafficking of Ca v channels may depend on the contact between S1-S2L and S5-S6L, and this is disrupted by the Gln to His mutation.
It is important to note that both the proband and his parents possess a variant in OTOG, a gene that is associated with autosomal-recessive moderate hearing impairment (15). Thus we cannot conclude whether it is the OTOG and/or CACNA1D variant that contributes to the hearing impairment of the affected individual. However, there are 12 reported mutations in OTOG of which 11 are lack of function, and only 1 reported missense mutation, making the mutation in this gene less of a candidate for the individual's symptoms. Considering the additional phenotypes of the proband, it is noteworthy that gain of function mutations in CACNA1D are associated with autism spectrum disorder (ASD) and epilepsy (7,8). Like Q567H, these mutations cause an increase in the window current mediated by Ca v 1.3. However, the mechanism involves a slowing of VDI and/or enhancement in voltage-dependent activation (7,8), in contrast to the selective effect of Q558H on suppressing VDI (Figs. 4 and 5). The gain of function in Ca v 1.3 would be expected to cause abnormalities in the morphology of dendritic spines (40), which have been linked to ASD (41) and epilepsy (42). Ca v 1.3 regulates neuronal excitability and synaptic plasticity in the striatum (43)(44)(45)(46), a brain region implicated in the neuropathology of ASD (47). However, we are cautious in noting that the impact of Q558H on Ca v 1.3 properties may be influenced by additional factors present in neurons that are absent from HEK293T cells. As has been demonstrated for ASD-linked mutations in Ca v 1.2 (48 -50), an understanding of how these gain of function mutations in Ca v 1.3 contribute neurodevelopmental and neurological disorders will be facilitated by the development of appropriate animal models.

Genetic analysis
This study abides by the Declaration of Helsinki principles and was approved by the Institutional Review Board at Hadassah-Hebrew University Hospital. Affected patients' guardians gave informed consent for exome sequencing and to publication of this study.

A conserved glutamine in voltage-gated Ca 2؉ channels
Whole exome sequencing was performed as described previously (51). Briefly, exonic sequences were enriched using Sure-Select Human All Exon 50 megabase kit (Agilent Technologies, Santa Clara, CA). Sequences were determined by HiSeq 2500 (Illumina, San Diego, CA) as 100-bp paired-end runs. Data analysis including read alignment and variant calling was performed by DNAnexus software (Palo Alto, CA) using the default parameters with the human genome assembly hg19 (GRCh37) as reference. Exome analysis of the analyzed individuals yielded 44.1 million mapped reads with a coverage of 86.0ϫ. Following alignment to the reference genome (hg19) and variant calling, variants were filtered out if they were offtarget (Ͼ8 bp from splice junction), synonymous (Ͼ3 bp from splice junction), or had minor allele frequency Ͼ0.01 in the ExAC database (Exome Aggregation Consortium, Cambridge, MA; URL: http://exac.broadinstitute.org) 3 or in our in-house exome database comprising ϳ2500 exomes. All potentially causative variants were confirmed using Sanger sequencing.

Molecular biology
The following cDNAs were used: Ca V 1.3 (GenBank TM no. AF370009), Ca V 2.1 (GenBank TM no. NM_001127221), Ca v 3.1 (AF190860), ␤ 2a (GenBank TM no. NC013684), and ␣ 2 ␦ 1 (Gen-Bank TM no. M76559.1). The Q558H mutation was inserted into the corresponding regions of the domain II S1-S2 linker (IIS1-S2L) in Ca V 1.3, Ca V 2.1, and Ca V 3.1 using the NEBuilder HiFi DNA Assembly cloning system (New England Biolabs). Channel fragments were amplified by PCR with appropriate primers and ligated into the parent plasmid (Ca V 1.3/pcDNA6, Ca V 2.1/pcDNA3.1, Ca V 3.1/pDsRed Express-N1). For the FLAG-tagged Ca V 1.3 WT and Q567H constructs, a FLAG epitope with spacer sequences (TRH-DYKDDDDK-VTFDEMQT) was added to the extracellular loop just C-terminal to residue Gly-693 by PCR using NEBuilder HiFi DNA Assembly cloning system (New England Biolabs). The vector was linearized by digestion with BsiWI and BamHI enzymes (New England Biolabs). The entire open reading frames of all the plasmids were verified by sequencing to confirm the presence of the mutation.

Cell culture and transfection
Human embryonic kidney cells transformed with SV40 T-antigen (ATCC, cat. no. CRL-3216, RRID: CVCL_0063) were maintained in Dulbecco's modified Eagle medium (Gibco, 11965) supplemented with 10% fetal bovine serum (Atlanta Biologicals, S11150), at 37°C in 5% CO 2 -95% air humidified atmosphere. For electrophysiological recordings, the cells were grown at ϳ70% confluent in 35-mm Petri dish prior to transient transfection with Ca v channel cDNAs using FuGene 6 transfection reagent (Promega, E2691). The following amounts of cDNAs were used: WT or mutant Ca V 1.3 or Ca V 2.1 (2.5 g) and 1 g each of the auxiliary subunits ␤ 2a and ␣ 2 ␦ 1 . Ca v 3.1 (2.5 g) was transfected without auxiliary subunit cDNAs. In all cases, cDNA encoding the GFP (pEGFP) was cotransfected to aid fluorescent detection of transfected cells. cDNA was incubated with 6 l of FuGene transfection reagent in 100 l of FBS-free DMEM in a sterile microfuge tube (1.5 ml). The transfection solution was added to the dish and cells were maintained at 37°C for 48 h. Two h before its electrophysiological recording the cells were dissociated using Versene reagent (Thermo Fisher Scientific, 15040066) in 35-mm Petri dishes with fresh medium at low density ϳ20%.
For cell-surface biotinylation assays, HEK293T cells were transfected using Lipofectamine 3000 reagent (Life Technologies, L3000-008) following the manufacturer's protocol. Plasmid DNA (WT or mutant Ca v 1.3 (1.2 g); ␤ 2a and ␣ 2 ␦ 1 (0.4 g each)) was diluted in Opti-MEM (50 l, Life Technologies) and 4 l of P3000 Reagent, and added to a mixture of Opti-MEM (50 l) and Lipofectamine 3000 Reagent (3 l) and allowed to incubate for 10 min at room temperature. The solution was then added to the cells. The cell culture media were changed after 24 h and cell-surface biotinylation assays were performed 48 h after transfection.

Electrophysiology
Transfected HEK293T cells were recorded with the wholecell patch clamp technique at room temperature (21-24°C) using an EPC-8 patch clamp amplifier and Patchmaster software (HEKA Elektronik). Electrophysiological data were sampled at 20 kHz and filtered at 5 kHz. Leak current was subtracted using an online P/4 protocol. The external solution contained the following (in mM): CaCl 2 10, TEA-Cl 125, HEPES 10, and glucose 10, pH 7.3, with TEA-OH. The internal solution contained (in mM): CsCl 110, MgCl 2 5, EGTA 10, HEPES 10, Na-ATP 4, GTP 0.1, pH 7.3, with CsOH. The recording pipettes had a resistance in the bath solution in the range of 4 -7 megohms. All reagents were purchased from Sigma-Aldrich.
Current-voltage relationships were obtained from the peak amplitude of currents evoked by 200-ms depolarizations from Ϫ80 mV to various voltages in 5-mV increments. I-V data were fit with the Boltzmann equation: I ϭ G⅐(V Ϫ E)/(1 ϩ exp(Ϫ (V Ϫ V1 ⁄ 2 )/k)), where G is conductance, E is the reversal potential, V1 ⁄ 2 is the voltage of half-maximal activation, and k is the slope factor for activation. The current density was obtained by normalizing the current amplitude with the membrane capacitance.
The kinetics of inactivation were measured by fitting traces of currents with a double exponential function: I ϭ A fast exp(Ϫt/ fast ) ϩ F slow exp(Ϫt/ slow ), where t is time, and A fast and A slow are the amplitudes of the current decaying with corresponding time constants (). The fraction of the total inactivating current (A fast ϩ A slow ) corresponding to fast inactivation was determined as F fast ϭ A fast /(A fast ϩ A slow ) and that for slow inactivation (F slow ) was determined as A slow /(A fast ϩ A slow ).
For measuring steady-state inactivation, a conditioning prepulse to various voltages preceded a test depolarization to ϩ10 mV for Ca V 1.3 and Ca V 2.1 and to Ϫ30 mV for Ca V 3.1. Inactivation curves were fitted with the Boltzmann function: I ϭ I max / (1 ϩ exp (V Ϫ V1 ⁄ 2 )/k), where I represents test current amplitude, V is test voltage, V1 ⁄ 2 is the voltage of half-maximal inactivation, and k is the slope factor. Voltage dependence of activation was measured by plotting the normalized tail currents evoked by 2-ms repolarizations to Ϫ60 mV from 10-ms pulses to vari-3 Please note that the JBC is not responsible for the long-term archiving and maintenance of this site or any other third party hosted site.
A conserved glutamine in voltage-gated Ca 2؉ channels ous voltages, followed by fitting with the above Boltzmann function. Data were analyzed using custom routines in Igor Pro software (WaveMetrics) and SigmaPlot (Systat Software Inc., San Jose, CA). Averaged data are shown as the mean Ϯ S.E. For statistical analysis, data were first analyzed for normality by Shapiro-Wilk test followed by Student's t test with differences considered statically significant when p Ͻ 0.05. For gating current analysis, the data were fit by linear regression and the slopes compared by F-test.

Double-label immunofluorescence
HEK293T cells were transfected as described for cell-surface biotinylation experiments except that FLAG-tagged Ca v 1.3 and Q567H cDNAs were used, and cells were plated on glass cover-slips. Sequential double-labeling for cell-surface and total channels was performed as described previously (53). Twenty-four h post transfection, the cells were rinsed in PBS and incubated with mouse anti-FLAG M2 antibody (Sigma-Aldrich) (1:50 diluted in DMEM media with 10% FBS) at room temperature for 30 min. The cells were then rinsed in PBS and fixed with 4% paraformaldehyde in PBS. After rinsing three times with PBS the cells were then blocked using 5% normal goat serum and 0.3% Triton X-100 in PBS for 30 min at room temperature and incubated with rabbit anti Ca v 1.3 antibodies to label the total intracellular and cell-surface channels for 1 h at room temperature. After three times with PBS, cells were incubated in secondary antibodies (RhoRed-X 568 anti-rabbit and Alexa Fluor 488 anti-mouse, Jackson ImmunoResearch Laboratories) (1:1000 each diluted in 5% normal goat serum and 0.1% Triton X-100 in PBS) for 1 h at room temperature. Coverslips were rinsed again three times in PBS and one time with water before mounting in Fluoromount G (Electron Microscopy Sciences, 17984 -25). Labeled cells were visualized with a 60ϫ oil-immersion objective on a confocal microscope (Fluoview 1000, Olympus). Images were acquired at 510 for 488 channel and 500 for 568 channel laser intensity and processed in ImageJ.