| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
|
|
|
|||||||
|
|||||||||
1Department of Medicine, University of Wisconsin, Madison, Wisconsin 53711; and 2Divisions of Cardiovascular Diseases and Pediatric Cardiology, Departments of Medicine, Pediatrics, and Molecular Pharmacology and Experimental Therapeutics, Mayo Clinic, Rochester, Minnesota 55905
Submitted 15 September 2003 ; accepted in final form 13 February 2004
 |
ABSTRACT |
|---|
 TOP ABSTRACT MATERIALS AND METHODSRESULTSDISCUSSIONGRANTSREFERENCES |
|---|
-subunit have challenged this viewpoint. In this report, we test the hypothesis that common missense HERG polymorphisms alter channel physiology. Comprehensive mutational analysis of HERG was performed on genomic DNA derived from a population-based cohort of sudden infant death syndrome and two reference allele cohorts derived from 100 African American and 100 Caucasian individuals. Amino acid-encoding variants were considered common polymorphisms if they were present in at least two of the three study cohorts with an allelic frequency >0.5%. Four nSNPs were identified: K897T, P967L, R1047L, and Q1068R. Wild-type (WT) and polymorphic channels were heterologously expressed in human embryonic kidney cells, and biochemical and voltage-clamp techniques were used to characterize their functional properties. All channel types were processed similarly, but several electrophysiological differences were identified: 1) K897T current density was lower than the other polymorphic channels; 2) K897T channels activated at more negative potentials than WT and R1047L; 3) K897T and Q1068R channels inactivated and recovered from inactivation faster than WT, P967L, and R1047L channels; and 4) K897T channels showed subtle differences compared with WT channels when stimulated with an action potential waveform. In contrast to K897T and Q1068R channels, P967L and R1047L channels were electrophysiologically indistinguishable from WT channels. All HERG channels had similar sensitivity to block by cisapride. Therefore, some HERG polymorphic channels are electrophysiologically different from WT channels.
ion channels; genetics; long QT syndrome; arrhythmia; sudden death
- or -subunits or ankyrin-B, with K+ channel defects accounting for the majority of LQTS. LQT2 is caused by mutations in human ether-a-go-go-related gene (HERG; KCNH2), which encodes the -subunit of the channel mediating the rapidly activating component of the delayed rectifier K+ current (IKr) (5, 21, 29). Most LQT2-causing gene defects are missense mutations, resulting in pathogenic single amino acid residue changes. The functional consequence of LQT2-linked mutations is thought to be a net reduction in IKr leading to a diminished "repolarization reserve" (19).
In addition to rare mutations, common polymorphisms also exist. Common polymorphisms have been defined as nucleotide substitutions found in both control and patient populations, usually at a frequency of 
1% or greater (31). When viewed in the context of pathological mutations, the presence of common nonsynonymous single nucleotide polymorphisms (nSNPs) in apparently healthy populations suggests that they are well tolerated and likely to have wild-type (WT)-like physiology. However, the identification of common nSNPs in the cardiac SCN5A Na+ channel and the MiRP1 K+ channel -subunit that alter channel physiology and drug sensitivity has challenged this viewpoint (25, 28). Indeed, these particular nSNPs in critical channels of the heart have a functional phenotype in vitro and may mediate genetic susceptibility to fatal ventricular arrhythmias in the setting of myocardial infarction or exposure to QT-prolonging medications.
Two nSNPs have been found within HERG (10, 11, 31). While more prevalent than disease-linked mutations, polymorphisms may produce subtle functional changes in HERG channel physiology that may alter the QT interval. In fact, the most common nSNP identified to date, HERG K897T (lysine, K, to threonine, T, at amino acid position 897) has been associated with altered channel biophysics and QT interval, although results differ between investigative groups (3, 10, 17, 18, 23).
In the present work, we studied whether nSNPs alter HERG channel function. To test this, we functionally characterized four nSNPs found in genomic DNA derived from a population-based cohort of sudden infant death syndrome (SIDS) and subsequently in two anonymous reference allele cohorts of distinct ethnicity (100 African-American and 100 Caucasian individuals) (1). Two of these amino acid variants were the previously identified arginine to leucine at amino acid position 1047 (R1047L) and K897T polymorphisms. The other two variants, proline to leucine at amino acid position 967 (P967L) and glutamine to arginine at amino acid position 1068 (Q1068R), are novel nSNPs. Western blot and whole cell patch-clamp techniques were used to investigate the effects of the amino acid changes on HERG channels heterologously expressed in a human cell line. The data show that some polymorphisms alter HERG channel physiology.
|
MATERIALS AND METHODS |
|---|
|
TOPABSTRACTMATERIALS AND METHODSRESULTSDISCUSSIONGRANTSREFERENCES |
|---|
Site-directed mutagenesis and heterologous expression. The appropriate nucleotide changes (A2690C, C2900T, G3140T, and A3203G) resulting in K897T, P967L, R1047L, and Q1068R polymorphisms, respectively, were engineered into WT HERG cDNA in the pCDNA3 vector (Invitrogen; Carlsbad, CA), and the integrity of the construct was verified by DNA sequencing. WT and polymorphic HERG channels were stably and transiently expressed in human embryonic kidney (HEK)-293 cells as described previously (32). HEK-293 cells were cultured in modified MEM (GIBCO/Invitrogen; Carlsbad, CA) with G418 (Invitrogen) antibiotic and kept at 37°C in 5% CO2.
Western blotting, electrophysiology, and cisapride. HEK-293 cells were grown to similar confluency before transient transfection with equal amounts of cDNA. Forty-eight hours after transfection, cells expressing each channel variant were prepared in parallel for Western blot analysis as previously described (32). Equal volumes of cells or total protein were electrophoresed on a 7.5% SDS-polyacrylamide gel, transferred to nitrocellulose membranes, and probed with an antibody specific for the distal COOH terminal of HERG (32).
A tight-seal whole cell recording technique was used for biophysical analysis of WT and polymorphic HERG channels (33). Pipettes had resistances between 1.5 and 2.5 M
when filled with recording solution. Series resistance compensation was 
70% in all experiments. Cell capacitance was estimated from the area within the capacitance transient. During action potential voltage-clamp recordings, HERG currents were separated from endogenous HEK-293 currents through off-line subtraction of the current response after application of 5 µM E-4031 (a specific HERG channel blocker) from the response obtained before drug application (33). Data were acquired with an Axopatch-2C amplifier controlled by Clampex 8.0 (Axon Instruments; Union City, CA). Current traces fit with exponential functions were acquired at 20 kHz with all other data acquired at 2 kHz.
Cisapride (Research Diagnostics; Flanders, NJ) stock concentration was 10 µM in ethanol. Working concentrations were obtained through serial dilutions with recording saline and used within 48 h. E-4031 (Alomone Labs; Jerusalem, Israel) was diluted in recording saline to a working concentration of 5 µM.
Statistics.
Data are given as means ± SE. ANOVA analyses determined whether significant biophysical differences existed between the channel variants at P < 0.05. When the ANOVA yielded a significant F-value, differences between specific channel variants were ascertained through pair-wise comparisons with Fisher's protected least-significant difference (LSD) test, which takes into account the effect of multiple comparisons. For the test, group i was determined to be different from group j if |
i – j| LSD, where i and j are the means of groups i and j and LSD is calculated as (t/2, dferror) x MSerror(1/ni+1/nj), where t/2 is the t-statistic at P < 0.05/2, dferror is the degree of freedom (total number of measurements in the pairwise comparison – 2), MSerror is the intragroup error calculated from ANOVA, and ni and nj are the numbers of measurements in groups i and j, respectively. Differences in allelic frequency were analyzed using Fisher's exact test.
|
RESULTS |
|---|
|
TOPABSTRACTMATERIALS AND METHODSRESULTSDISCUSSIONGRANTSREFERENCES |
|---|
|
|
Channel biophysical properties. During voltage-clamp protocols, as well as cardiac action potentials, HERG channels undergo distinct voltage- and time-dependent gating transitions. When heterologously expressed in HEK-293 cells and subjected to voltage-clamp analysis, WT and polymorphic channels showed qualitatively similar current waveforms (Fig. 2A). Cells were depolarized from a holding potential of –80 mV to voltages between –70 and 50 mV in 10-mV increments for 3 s, followed by a step to –50 mV for 2.85 s to elicit Itail. All channel variants exhibited outward currents in response to steps positive to –50 mV that peaked around 0 mV, declined upon further depolarization, and showed a prominent Itail upon repolarization to –50 mV.
|
To measure channel inactivation, cells were depolarized from a holding potential of –80 to 50 mV for 1.5 s, hyperpolarized to –100 mV for 2.5 ms, and then stepped to test voltages between 60 and –40 mV in 10-mV decrements for 1.5 s. The test voltages generate large amplitude outward currents, the decline of which (arrow, Fig. 3A, inset) was fit as a single-exponential process to derive a time constant for channel inactivation (33). Inactivation time constants were plotted as a function of test voltage (Fig. 3A), and statistical analyses showed differences between channel variants (n = 13, 12, 12, 13, and 9 cells for WT, K897T, P967L, R1047L, and Q1068R channels, respectively). The most striking difference was that both K897T and Q1068R channels inactivated more rapidly than WT, P967L, or R1047L channels (P < 0.05) while being indistinguishable from each other. Specifically, K897T channel inactivation was faster than WT channels between –30 and 60 mV, faster than P967L channels between –40 and 60 mV, and faster than R1047L channels between –30 and 60 mV (P < 0.05). Q1068R channel inactivation was faster than WT channels at –20 mV and between 0 and 60 mV, faster than P967L channels between –30 and 60 mV, and faster than R1047L channels between –20 and 60 mV (P < 0.05). No significant differences were present between P967L, R1047L, or WT channel inactivation.
|
Deactivation was measured using the same voltage-clamp protocol as recovery from inactivation except that the decay of Itail (arrow, Fig. 3C, inset) was fit as a double-exponential process (33). The fast and slow time constants derived from the fits represent the fast and slow components of channel deactivation. The mean fast and slow deactivation time constants were plotted as a function of test voltage and showed no consistent difference between the five channel variants (Fig. 3C).
An action potential voltage-clamp protocol was used to further compare the biophysical properties of K897T and WT channels (Fig. 4). As previously described (33), depolarization by the action potential protocol first modestly increases HERG current through channel activation and subsequent inactivation. During repolarization, HERG current increases further to its maximal level through rapid recovery of channels from inactivation. After repolarization, HERG current gradually decreases through channel deactivation. For each cell, current was normalized to total charge (pC) calculated by integrating the area under the current response (from 0 to 700 ms; upward arrows, Fig. 4A). For comparisons between WT and K897T channels, normalized current responses were averaged for all cells expressing each channel type (n = 8 cells each). The mean current response for WT and K897T HERG channels is shown in Fig. 4A, with the voltage-clamp protocol superimposed in light gray. Recordings from single cells expressing WT (top) or K897T (bottom) channels are shown in the left insets. Whereas the normalized peak current values are virtually identical (4.54 ± 0.24 and 4.49 ± 0.21 pA/pC for WT and K897T, respectively), the two waveforms are not. Subtracting the normalized mean WT response from that of K897T (Diff trace, Fig. 4A, right inset) suggests that K897T currents are smaller during the upstroke of the current response with depolarization and larger during the decay of the current response after repolarization relative to WT currents. The statistical significance of this difference was analyzed by partitioning the current responses into sequential 100-ms windows (i.e., 0–100 ms, 100–200 ms, etc), determining the time point of maximal difference between the two current waveforms (downward arrows, Fig. 4A), and then assaying whether K897T was smaller during the upstroke and larger during the downstroke of the current response relative to the WT current level. During the upstroke of the current response, K897T current levels of 2.20 ± 0.15 and 3.43 ± 0.17 pA/pC at 256.3 and 329.2 ms, respectively, were significantly smaller than the corresponding WT current levels of 2.53 ± 0.06 and 3.80 ± 0.11 pA/pC (P < 0.05). During the decay of the current response after repolarization, the K897T current level of 1.66 ± 0.28 pA/pC at 418.3 ms was significantly larger than the WT current level of 1.12 ± 0.09 pA/pC (P < 0.05). All other assayed time points showed statistically similar current levels (Fig. 4B).
|
|
|
DISCUSSION |
|---|
|
TOPABSTRACTMATERIALS AND METHODSRESULTSDISCUSSIONGRANTSREFERENCES |
|---|
Functional analysis is a key step to understanding the impact of amino acid coding variants on protein function, and a strength of our study is that we compared functional characteristics of five channels (four polymorphisms and WT) under identical experimental conditions. Western blot analysis showed the 135- and 155-kDa bands associated with normal HERG channel processing. Thus unlike many LQT2-associated disease-causing mutations that fail to traffic normally (32), all of the polymorphic channels underwent biochemical processing similar to WT channels.
Current density measurements from polymorphic and WT channels showed significant differences. The most important finding was that HERG K897T channels generated less current than the other polymorphic channels, whereas they did not achieve significantly different levels with respect to WT channels.
Biophysical differences between WT and polymorphic channels were also present. The V1/2 of channel activation for K897T is 5–6 mV negative to that of WT and R1047L channels, and thus K897T channels open at more negative voltages. Furthermore, K897T and Q1068R channels inactivate and recover from inactivation faster than the other channel variants, showing 10- to 20-mV shifts in the voltage dependence of these time constants. Responses to the action potential voltage-clamp protocol suggest that the biophysical changes present in K897T channels may result in slight alterations (10–30%) of HERG current during the action potential. Compared with WT channels, the small decrease in K897T channel current during the action potential could cause a subtle increase in action potential duration, whereas after repolarization the small increase in current could alter terminal repolarization and the subsequent return of excitability. Channel deactivation was not consistently different between WT and polymorphic channels, nor were the IC50 values or Hill coefficients for block by cisapride. Thus in addition to varying from each other, some polymorphic channels (K897T and Q1068R), but not all (P967L and R1047L), have detectable differences when compared with WT channels. Our data also suggest that the allelic frequency of each polymorphism was not a reliable predictor of altered channel function.
The four HERG polymorphisms reside in the distal COOH terminal of the channel protein. Eleven LQT2-linked missense mutations have also been identified within the COOH terminal (2, 6, 13, 22, 27, 31). Three of these, R752W, S818L, and V822M, have been functionally characterized and shown to be trafficking defective (6, 16, 32). A portion of the COOH terminal region containing amino acids 1018–1122 has been shown to be necessary for HERG channel expression (9), whereas an overlapping region consisting of amino acids 1063–1159, as well as the point mutation V822M, disrupt binding to GM-130, a resident Golgi protein thought to be involved in protein trafficking (20). Thus noncontinuous regions of the COOH terminal are involved in regulating channel trafficking and the small differences in current densities observed in this study could arise from perturbations in positive and/or negative regulation of channel trafficking.
HERG channel inactivation is altered by missense mutations in the pore region as well as deletions of the NH2 terminal (24, 26, 30). The effects of K897T and Q1068R demonstrated here provide evidence that the COOH terminal may also affect HERG channel inactivation. Thus, similar to the cardiac SCN5A Na+ channel where the COOH terminal of the channel is host to multiple functional domains that regulate inactivation and current expression (4), the COOH terminal of HERG channels may also be involved in multiple aspects of channel physiology. This unexpected role of the HERG COOH terminal region stresses the importance of functional characterization of genetic variants.
Our data are the first to provide functional analysis regarding the biophysical properties of three of the HERG polymorphisms, P967L, R1047L, and Q1068R. Other studies have provided conflicting data regarding K897T channel biophysical properties. Scherer and colleagues (23) did not detect differences between K897T and WT current amplitude, biophysical properties, or sensitivity to drug block by the antihistamine terfenadine. Bezzina and colleagues (3) also found no differences in current amplitude or channel inactivation properties but did find K897T channels to have a hyperpolarized shift in channel activation as well as faster activation and deactivation kinetics at specific voltages. Paavonen and colleagues (17) did not detect changes in channel activation properties but did find decreased K897T protein levels on Western blot, slower channel deactivation and inactivation kinetics, and a hyperpolarized shift in steady-state inactivation. While our results provide new data, they do not resolve the reported differences, and in part the discrepancies may arise from differences in expression systems and/or experimental protocols.
The potential clinical impacts of P967L, R1047L, or Q1068R are unknown, whereas those associated with K897T status have been reported. Similar to discrepancies in K897T channel biophysics, the reports of K897T clinical phenotypes are not in complete agreement. Laitinen et al. (10) suggested a correlation between K897T and QT interval in female LQT1 patients. Pietilla et al. (18) found K897T to be associated with a longer maximal QTc interval and increased dispersion of ventricular repolarization in females but not males within a random middle-aged study group. Paavonen et al. (17) did not find gender-specific differences in resting QTc intervals for K897T carriers but did show an increased QTc interval during exercise for a subset of carriers with a prolonged baseline QTc. Bezzina et al. (3), on the other hand, found females homozygous for T897 to have shorter QTc intervals than either heterozygous or noncarriers. While our results provide potential mechanisms for altered IKr in K897T carriers (altered channel kinetics and slight changes in current levels), they do not resolve the apparent clinical or gender differences. The different biophysical and clinical phenotypes attributed to HERG K897T highlight the need for consistent experimental protocols when investigating the potentially subtle effects of polymorphisms.
Our data also provide insight in another important area, which is that the four HERG polymorphisms we studied do not directly convey heightened sensitivity to block by cisapride, a drug known to cause acquired LQTS by high-affinity HERG channel block (15). Thus unlike the T8A polymorphism in MiRP1 (a putative HERG channel -subunit) that conveys heightened channel sensitivity to drug block (25), these HERG -subunit polymorphisms do not do this. However, subtle modification of IKr by HERG polymorphisms such as K897T, whether through altered current density or channel kinetics, could indirectly contribute to proarrhythmic effects of HERG K+ channel blockers by reducing the repolarization reserve.
In summary, we studied WT and four nonsynonymous single nucleotide polymorphisms found in HERG K+ channels. Functional characterization of these polymorphisms demonstrated their WT-like biochemical processing with normal trafficking to the cell surface membrane. The expressed currents resemble WT current, although small differences in current density and biophysical properties were detectable within some of the polymorphisms. Sensitivity to block by cisapride was similar for WT and all polymorphic channels. These data suggest that some HERG polymorphisms are WT channel like in their physiological properties, whereas others have detectable differences. These findings may provide a basis for new insight into clinical phenotypes and arrhythmia risk prediction of polymorphism carriers; however, further correlative studies are clearly needed.
Limitations.
All polymorphisms were studied as homomeric and not hetermeric channels. Extrapolation of our results to heterozygous carriers must done cautiously and is beyond the scope of this work. DNA samples for these studies were obtained from either a SIDS-based tissue cohort or genomic samples obtained from apparently healthy individuals. Relevant clinical data do not exist for these cohorts to permit comparison of the HERG channel physiology with cardiac phenotypes.
|
GRANTS |
|---|
|
TOPABSTRACTMATERIALS AND METHODSRESULTSDISCUSSIONGRANTSREFERENCES |
|---|
|
ACKNOWLEDGMENTS |
|---|
|
FOOTNOTES |
|---|
The costs of publication of this article were defrayed in part by the payment of page charges. The article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
* B. D. Anson and M. J. Ackerman contributed equally to this work. 
|
REFERENCES |
|---|
|
TOPABSTRACTMATERIALS AND METHODSRESULTSDISCUSSIONGRANTSREFERENCES |
|---|
This article has been cited by other articles:
|
|
 |
|
  M. F. Sinner, A. Pfeufer, M. Akyol, B.-M. Beckmann, M. Hinterseer, A. Wacker, S. Perz, W. Sauter, T. Illig, M. Nabauer, et al. The non-synonymous coding IKr-channel variant KCNH2-K897T is associated with atrial fibrillation: results from a systematic candidate gene-based analysis of KCNH2 (HERG) Eur. Heart J., April 1, 2008; 29(7): 907 - 914. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
P. T. Ellinor and D. J. Milan Polymorphisms and atrial fibrillation: sorting the wheat from the chaff Eur. Heart J., April 1, 2008; 29(7): 843 - 845. [Full Text] [PDF]
|
|
|
|
 |
|
 C. Newton-Cheh, C.-Y. Guo, M. G. Larson, S. L. Musone, A. Surti, A. L. Camargo, J. A. Drake, E. J. Benjamin, D. Levy, R. B. D'Agostino Sr, et al. Common Genetic Variation in KCNH2 Is Associated With QT Interval Duration: The Framingham Heart Study Circulation, September 4, 2007; 116(10): 1128 - 1136. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
M. Arnestad, L. Crotti, T. O. Rognum, R. Insolia, M. Pedrazzini, C. Ferrandi, A. Vege, D. W. Wang, T. E. Rhodes, A. L. George Jr, et al. Prevalence of Long-QT Syndrome Gene Variants in Sudden Infant Death Syndrome Circulation, January 23, 2007; 115(3): 361 - 367. [Abstract] [Full Text] [PDF]
|
|
|
|
 |
|
 E. Schulze-Bahr Arrhythmia Predisposition: Between Rare Disease Paradigms and Common Ion Channel Gene Variants J. Am. Coll. Cardiol., October 27, 2006; 48(9_Suppl_A): A67 - A78. [Abstract] [Full Text] [PDF]
|
|
|
|
 |
|
 S. Rajamani, C. L. Anderson, C. R. Valdivia, L. L. Eckhardt, J. D. Foell, G. A. Robertson, T. J. Kamp, J. C. Makielski, B. D. Anson, and C. T. January Specific serine proteases selectively damage KCNH2 (hERG1) potassium channels and IKr Am J Physiol Heart Circ Physiol, March 1, 2006; 290(3): H1278 - H1288. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
C. L. Anderson, B. P. Delisle, B. D. Anson, J. A. Kilby, M. L. Will, D. J. Tester, Q. Gong, Z. Zhou, M. J. Ackerman, and C. T. January Most LQT2 Mutations Reduce Kv11.1 (hERG) Current by a Class 2 (Trafficking-Deficient) Mechanism Circulation, January 24, 2006; 113(3): 365 - 373. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
L. Crotti, A. L. Lundquist, R. Insolia, M. Pedrazzini, C. Ferrandi, G. M. De Ferrari, A. Vicentini, P. Yang, D. M. Roden, A. L. George Jr, et al. KCNH2-K897T Is a Genetic Modifier of Latent Congenital Long-QT Syndrome Circulation, August 30, 2005; 112(9): 1251 - 1258. [Abstract] [Full Text] [PDF]
|
|
|
|
 |
|
 D. J. Tester and M. J. Ackerman Sudden infant death syndrome: How significant are the cardiac channelopathies? Cardiovasc Res, August 15, 2005; 67(3): 388 - 396. [Abstract] [Full Text] [PDF]
|
|
|
|
 |
|
 A. Pfeufer, S. Jalilzadeh, S. Perz, J. C. Mueller, M. Hinterseer, T. Illig, M. Akyol, C. Huth, A. Schopfer-Wendels, B. Kuch, et al. Common Variants in Myocardial Ion Channel Genes Modify the QT Interval in the General Population: Results From the KORA Study Circ. Res., April 1, 2005; 96(6): 693 - 701. [Abstract] [Full Text] [PDF]
|
|
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH< |
) for channel inactivation are plotted as a function of test voltage. Inset shows a portion of the voltage-clamp protocol (top traces) and representative currents (bottom traces) highlighting inactivation (arrow). Scale is 1 nA by 10 ms. B: channel recovery from inactivation. Mean time constants for channel recovery from inactivation are plotted as a function of test voltage. Inset shows a portion of the voltage-clamp protocol (top traces) and representative currents (bottom traces) highlighting channel recovery (arrow). Symbols are as in A. Scale is 500 pA by 10 ms. C: channel deactivation. Mean time constants for channel deactivation are plotted as a function of test potential. Inset shows a portion of the voltage-clamp protocol (top traces) and representative currents (bottom traces) highlighting channel deactivation (arrow). Symbols are as in A with open and filled symbols representing the slow and fast time constants, respectively. Scale is 500 pA by 100 ms.
