© 2001 by European Society of Cardiology
Copyright © 2001, European Society of Cardiology
Cellular consequences of HERG mutations in the long QT syndrome: precursors to sudden cardiac death
Cardiac Bioelectricity Research and Training Center, Department of Physiology and Biophysics, Department of Biomedical Engineering, 509 Wickenden Building, Case Western Reserve University, Cleveland, OH 44106-7207, USA
* Corresponding author. Tel.: +1-216-368-4051; fax: +1-216-368-8672 yxr{at}po.cwru.edu
Received 1 September 2000; accepted 15 November 2000
| Abstract |
|---|
|
|
|---|
Background: A variety of mutations in HERG, the major subunit of the rapidly activating component of the cardiac delayed rectifier IKr, have been found to underlie the congenital Long-QT syndrome, LQT2. LQT2 may give rise to severe arrhythmogenic phenotypes leading to sudden cardiac death. Objective: We attempt to elucidate the mechanisms by which heterogeneous LQT2 genotypes can lead to prolongation of the action potential duration (APD) and consequently the QT interval on the ECG. Methods: We develop Markovian models of wild-type (WT) and mutant IKr channels and incorporate these models into a comprehensive model of the cardiac ventricular cell. Results: Using this virtual transgenic cell model, we describe the effects of HERG mutations on the cardiac ventricular action potential (AP) and provide insight into the mechanism by which each defect results in a net loss of repolarizing current and prolongation of APD. Conclusions: This study demonstrates which mutations can prolong APD sufficiently to generate early afterdepolarizations (EADs), which may trigger life-threatening arrhythmias. The severity of the phenotype is shown to depend on the specific kinetic changes and how they affect IKr during the time course of the action potential. Clarifying how defects in HERG can lead to impaired cellular electrophysiology can improve our understanding of the link between channel structure and cellular function.
KEYWORDS Arrhythmia (mechanisms); Computer modelling; K-channel; Long QT syndrome; Sudden death; Ventricular arrhythmias
| 1 Introduction |
|---|
|
|
|---|
Increasingly, seemingly idiopathic syndromes are being linked to genetic defects. A recent example is the congenital Long-QT (LQT) syndrome [1–3]. The syndrome is associated with prolongation of the Q-T interval on the ECG, which may precede syncope and sudden cardiac death. In 1995, patients from a family afflicted with a genetically linked form of the Long-QT syndrome (LQT2) were found to have abnormalities in the HERG gene [4]. This gene underlies a potassium channel that is highly expressed in the heart with properties similar to the E-4031 sensitive current [5,6]. Many experiments have since demonstrated that HERG encodes the major subunit of the rapidly activating component of the cardiac delayed rectifier current (IKr). Very recently the minK-related peptide 1 (MiRP1) was shown to stably assemble with HERG and give rise to a current strikingly similar to native IKr [7].
We develop a structurally based model of cardiac IKr in order to elucidate the functional consequences of gene defects in cardiac HERG that lead to abnormalities in repolarization. We have drawn from a variety of experiments presented in the peer-reviewed literature in order to build a robust model that accounts for the breadth of published data. More than 50 mutations in HERG have been linked to the congenital Long-QT syndrome [3]. Defects in the HERG pore have been shown to have heterogeneous cellular phenotypes. Pore mutations may result in a loss of function, sometimes due to a trafficking defect [8], and may or may not coassemble with WT HERG subunits to exert dominant negative effects [9]. Other defects in the channel pore give rise to altered channel kinetics leading to decreased repolarizing current [10,11]. Nearby mutations in the S4–S5 linker have been shown to variably affect activation [12].
Here, we use Markov models of wild-type (WT) and mutant IKr in a cardiac ventricular cell model to evaluate the cellular arrhythmogenic consequences of several HERG mutations associated with LQT2. The Markov models allow for the introduction of state-specific kinetic changes due to mutation induced structural defects. These changes affect adjacent channel states through coupling. We choose as examples two well-characterized defects that alter cell electrophysiology through their effects on channel kinetics. The first, the T474I point mutation, results in an altered voltage dependence of activation coupled with a reduction in macroscopic current density. The second, the R56Q mutation, is a point mutation in the Per—Arnt–Sim (PAS) domain located in the amino-terminus (N-terminus) region of HERG, which normally interacts with the channel and reduces the rate of deactivation [13]. Mutations in this region presumably prevent proper association of the N-terminus with the channel and act to increase the rate of channel deactivation. We choose these two defects as examples of mutations that exert their effects on cellular behavior by reducing the IKr repolarizing current through alterations in its channel kinetics.
We have also investigated the effects of the N629D mutation, a gain of function defect in HERG that results in loss of C-type inactivation coupled with loss of K+ selectivity [14]. HERG contains a pore selectivity sequence that is changed by the mutation from GFGN to GFGD, allowing for non-specific passage of monovalent cations. In this study, we incorporate these changes into our model of IKr and investigate the effects on the morphology of the action potential and action potential duration.
In addition, we elucidate the effects of a range of reduction in current on the action potential duration (APD), since many HERG defects exert their effects through a dominant negative mechanism, trafficking defects, or reduced mutant protein expression [3]. Finally we evaluate the electrophysiological consequences of a functional knockout mutation, where IKr current is absent, a situation described clinically in a patient born with such a defect [15]. In all, we provide several examples of heterogeneous genotypes and the mechanism by which each prolongs APD leading to a prolonged QT interval in the ECG, increased dispersion of repolarization, and possibly the development of early afterdepolarizations (EADs) and cardiac arrhythmias.
| 2 Methods |
|---|
|
|
|---|
The general approach to modeling the action potential (AP) is the same as that described for the dynamic Luo—Rudy (LRd) model of a ventricular cell [16] with modifications described in later publications [17,18]. The fast sodium current (INa) was recently reformulated using a Markovian approach and incorporated into the whole-cell model [19]. We use a similar approach to reformulate IKr. Macroscopic current density is given by:
|
| (1) |
PO,s is the sum of all channel open probabilities, Vm is the membrane potential, and Erev is the reversal potential. Gs is the maximum membrane conductance (channel density (
) times the unitary channel conductance (gs)). The channel state probabilities are described by first order differential equations using the methods described in [20]. The model parameters are fit using averaged measured activation, inactivation and deactivation data [5,6,21–23]. Parameters are then optimized to account for the proper behavior of macroscopic currents compared to data recorded at physiological temperatures and ion concentrations from guinea pig ventricular myocytes (see Figs. 3 and 4
).
|
|
All the simulations were encoded in C/C++. Simulations were implemented (double precision) on a Sun Workstation Ultra 1. A time step of 0.005 ms was used during the stimulus and AP upstroke. At all other times, a 0.01-ms time step was used.
The model for cardiac IKr is shown in Fig. 1. The model includes three closed states (C3, C2, C1), an open state (O), and an inactivation state (I). The transition between C2 and C1 is voltage independent [21]. Closed state inactivation (C1
I) is included from a single closed state, allowing for reduced computation time. The transition rate between C1 and I is the same as from C1 to O [24]. IKr exhibits strong inward rectification due to a rapid voltage dependent C-type inactivation at positive membrane potentials [11] (transition from O to I). The transitions between O and I are also dependent on extracellular potassium concentration [K+]out [25]. As [K+]out is increased the transition rates decrease. In addition, the maximum conductance of IKr (GKr) depends on [K+]out. These dependences on [K+]out are introduced into the model.
|
Action potential simulations are carried out in epicardial and midmyocardial (M) cells. Epicardial and M cells are simulated by varying the membrane conductance of IKs (GKs) as described previously [18]. Equations of the LRd model can be found in previous publications [16,17,26,27] and can be downloaded from the research section of http://www.cwru.edu/med/CBRTC
| 3 Results |
|---|
|
|
|---|
3.1 Simulated wild type (WT) behavior
The overall result of decreased rates of inactivation and recovery from inactivation at high [K+]out is an inactivation state that is more absorbing. This is demonstrated in the single channel gating of IKr in Fig. 2 (left panels) and channel mean open times (MOTs) (right panels). Panel A (left) shows experimentally recorded single channel gating of HERG channels at 21°C and 100 mM [K+]out [24]. Single channels are depolarized from a holding potential of –80 to 100 mV for 300 ms and subsequently hyperpolarized (arrow) to –120 mV for 100 ms (inset, bottom of Fig. 2). Experimental unitary HERG currents (panel A, left) are compared to simulated single channel gating of IKr (panel B, left) at 100 mM [K+]out [24]. The simulated single channel gating is at 37°C. The single channel gating behavior in the experiments and the simulations demonstrates channel inactivation from closed states and recovery from inactivation through the open state. Channels may fail to open during depolarization (single channel traces before arrow), due to inactivation from closed states, but conduct an inward current during subsequent hyperpolarization (traces after arrow) due to recovery from channel inactivation through the open state. The experimental traces show a longer latency to first opening compared to the simulations, an expected behavior at the lower temperature under which the experiments are conducted. The model predicts gating behavior at physiological [K+]out of 4.5 mM in panel C. Here, the inactivation state is less absorbing compared to panels A and B, an overall result of increased rates of inactivation and recovery from inactivation at lower [K+]out. Hence, more channel openings are observed during the depolarization step and the subsequent hyperpolarization (panel C, left). This behavior is consistent with macroscopic current recordings at physiological [K+]out, which have a persistent component of outward current during depolarization (see Fig. 4: experiment, top traces; simulation, bottom traces). This persistent current rapidly reaches a pseudo-steady state due to the rapid equilibrium between open and inactivation states. The MOTs of experimentally recorded HERG channels at [K+]out=100 mM [24] (panel A, right) and simulated IKr channels at [K+]out=100 mM and [K+]out=4.5 mM (right panels B and C, respectively) are shown for comparison. The simulated MOTs are computed as 1/(sum of transition rates leaving the open state, ββ+βi). At physiological [K+]out=4.5 mM, MOTs are shorter due to the increased inactivation rate. This comparison demonstrates that both the simulated and experimentally obtained MOTs exhibit weak voltage dependence over a wide range of voltage regardless of [K+]out [6,24].
|
The voltage dependence of activation and rectification (inactivation) of macroscopic IKr is shown in Fig. 3. The activation curve (panel A) is constructed by normalizing peak tail current at the indicated test potential to the maximum tail current during repolarization (–100 mV). The solid line is the simulation; filled boxes are experimental data obtained from the E-4031 sensitive current in guinea pig ventricular myocytes [5]. The rectification (R) curve (panel B) is computed as:
|
|
Computation of the rectification (inactivation) curve using this method more accurately describes channel inactivation than a standard inactivation protocol since IKr availability can be contaminated by deactivation at membrane potentials below 0 mV [7,11]. Rectification is removed in the simulation by making the O
I transition rate equal to zero. The simulated rectification curve is compared to experimental data from Sanguinetti and Jurkiewicz (filled circles) [5]. Rectification is obtained experimentally by extrapolating the linear portion of the peak tail current I–V curve into the positive potential range [5].
GKr is affected by variations in [K+]out. In Fig. 3C, the peak tail current—voltage relationship is shown for [K+]in=150 mM and [K+]out=5.4, 50 and 150 mM. Simulated I–V curves (top) are compared to experimentally recorded rabbit IKr (bottom) using the same protocol as in Shibasaki [6]. From a holding potential of –43 mV, the cell is pulsed to the test potential for 500 ms. A subsequent pulse to –43 mV for [K+]out=5.4, and –73 mV for [K+]out=50 and 150 mM elicits a tail current that is plotted as a function of the test potential. Notice that at [K+]out=5.4 mM, the tail current is outward at all membrane potentials (Erev=–89 mV) and exhibits strong inward rectification at positive potentials that compares well with the degree of rectification seen experimentally. At [K+]out=50 mM and [K+]out=150 mM, the tail current is inward in both the simulation and the experiment due to the equilibrium potentials of –29 and 0 mV at 50 and 150 mM [K+]out, respectively.
Fig. 4 (panel A) contains experimentally recorded IKr current traces from guinea pig ventricular myocytes [5]. Cells are depolarized for 250 ms from a holding potential of –40 mV and are subsequently repolarized to –40 mV for 750 ms (see protocol, inset). Simulated currents following the same protocol are shown in panel B. The tail current is plotted on a reduced time scale in order to be consistent with the experimental traces [5]. At a test potential of 10 mV or greater, hyperpolarization elicits an outward tail current hook experimentally and in the simulation (arrow). This hook is absent at test potentials 0 mV (arrow) or smaller. The hook reflects recovery from inactivation through the open state. Since recovery from inactivation (I
O) is very fast compared to deactivation (O
C1), a large outward tail current is conducted as channels slowly return to the closed resting states.
3.2 Simulated HERG mutations
The T474I mutation, a point mutation in the S2–S3 linker, alters gating properties and reduces macroscopic IKr [28,29]. We simulate these changes by altering the voltage dependence of the activation transition rates (
and 
, see Fig. 1 and Appendix) and reducing GKr by 35%. Fig. 5 compares simulated (panel B) activation and I–V curves for WT and T474I mutant IKr to experimentally measured (panel A) curves. The activation curves (protocol described for Fig. 3) are shown in the left panels. The experimentally measured V1/2 for WT and T474I is –15.9 and –43.2 mV, respectively (a shift of –27.3 mV for T474I relative to WT) [28]. The T474I simulated activation curve is shifted by –24 mV (V1/2=–47 mV) relative to the WT curve (V1/2=–23 mV). Experimental recordings were performed at 22–23°C, while the simulations were performed at physiological temperature (37°C). The leftward shift of the simulated activation curves relative to measured curves is an expected result of the increased temperature [22]. The fully activated I–V relationships of simulated IKr and experimentally recorded HERG are shown in the right panels of Fig. 5 (protocol described in Zhou et al. [28]). The mutant peak tail current is shifted by –20 mV relative to WT in the experimentally recorded and simulated curves. The maximum mutant current is 65% of WT levels.
|
Mutations in the PAS domain of the amino terminal of HERG act to increase the rate of deactivation [30]. We focus on the R56Q mutation since this mutant can increase the rate of deactivation most profoundly [30]. The mutation is simulated by increasing the O
C1 transition rate by a factor of 6.3 and the C2
C3 rate by a factor of 10.5 [30]. These increased rates of deactivation result in a rightward shift in the voltage dependence of activation and rectification (not shown), due to the strong coupling between the discrete states. The simulated kinetic changes result in shifts of +18 mV in the activation curve and +16 mV in the rectification curve. These shifts are also seen experimentally in HERG channels containing the R56Q mutation expressed in Xenopus oocytes [30] (not shown).
Fig. 6 compares action potentials from WT, T474I and R56Q simulated M-cells. The 1000th AP is shown at a cycle length (CL) of 750 ms. The left panel illustrates a WT M-cell action potential (AP). The AP duration at 90% repolarization (APD90) is 164 ms. The corresponding IKr current during the AP is shown in panel B. Channel state probabilities during the AP are shown in panels C through G (see model in Fig. 1). During the AP upstroke, channels move from the leftmost closed state (C3, panel C) and pass through C2 and C1 (panels D and E, respectively). From the C1 state, channels can pass directly to the open state (O, panel F) or to an unavailable inactivation state I (panel G). Once the channels open, they rapidly inactivate (transition from O
I). A pseudo-equilibrium is reached between the inactivation transition (O
I) and recovery from inactivation (I
O) that begins to favor recovery as the AP plateau repolarizes. As recovery from inactivation becomes more likely, IKr increases to a pronounced peak (arrow) and contributes to the faster repolarization late in the AP. At this time channels slowly deactivate (O
C1 transition).
|
In Fig. 6, the middle panel shows IKr and associated channel state probabilities during the AP for a T474I mutant M-cell (compare to WT). In addition to the kinetic changes, GKr is reduced by 35% in the simulation. Mutant APD is 16 ms longer than WT. Interestingly, the reduction in GKr, rather than altered channel kinetics, is the major factor in lengthening APD. The alterations in activation kinetics alone have only a minor effect on APD (not shown). The morphology of IKr early in the AP is altered due to the shift in the voltage dependence of activation. The leftward shift of the activation curve results in early channel openings at less depolarized membrane potentials (panel F, compare to WT, panel F) and an increased rate of IKr activation (panel B). Channels in the open state are available to inactivate earlier (panel G), resulting in the large notch following activation at the beginning of the current trace (panel B) which reflects fast inactivation. The channel open probability late in the AP (panel F), when IKr plays a major role in repolarization, displays a pronounced peak (arrow) and is similar to that of WT IKr.
The effect of the R56Q mutation on the cardiac AP is shown in an M-cell in the right panel of Fig. 6. While initial IKr elicited during the upstroke of the AP is comparable between the WT and R56Q mutant cells (panel B, compare to WT), the late current in the mutant is decreased due to increased probability of the O
C1 transition. This is seen in the open probability during the AP, shown in panel F. In contrast to WT and T474I mutant currents, the R56Q mutant lacks a peak in P(O) late in the AP (panel F, arrow). Deactivation occurs earlier, resulting in less repolarizing current late in the AP where IKr usually plays a major role in repolarization and determining APD.
Fig. 7 shows the rate dependence of APD in WT (panel A) and the two kinetic mutants T474I (panel B) and R56Q (panel C). There is only a slight prolongation at slow rate in mutant epicardial cells compared to WT (top). However, mutant M cells show significantly greater APD prolongation with decreasing rate (Fig. 7, bottom and table). The R56Q mutation gives rise to a prolonged APD which develops EADs (panel C, bottom inset) at CL=2000 ms (bradycardia).
|
Fig. 8 demonstrates the effects of the N629D gain of function mutation on the action potential in epicardial (panel A) and M-cells (panels C and E). All cells are paced at a CL=1000 ms for 1000 beats. The corresponding IKr is shown beneath each AP trace (panels B, D and F). The mutation is characterized by a loss of C-type inactivation and of ion selectivity, allowing for the passage of Na+ ions in addition to K+ ions. The relative selectivity of the channel for Na+ and K+ (PNa/PK) is 0.65 [14]. We simulate the mutation by eliminating C-type inactivation (O
I transition rate=0) and by altering selectivity to the ratio above. Using the Nernst equation, we compute the reversal potential of the N629D IKr to be –13 mV. This reversal potential for mutant IKr falls within the range of physiological plateau potentials. As a result, early IKr is an outward repolarizing current carried by K+. However, when the membrane potential repolarizes below –13 mV, the current is inward and carried by Na+. This is demonstrated in the IKr traces during epicardial APs shown in Fig. 8, panel B. During the initial phase of the action potential IKr is outward, but once the repolarizing membrane potential passes the reversal potential threshold, the current becomes inward (arrows). The inward current late in the plateau results in severe APD prolongation. Even more severe, are the effects of this mutation in M-cells (panel C and E). In panel C, we show an M-cell with GKs:GKr=7:1, and in panel E an M-cell with GKs:GKr=12:1, both of which fall within the physiological range for currents measured in M-cells. In both M-cells, the inward IKr carried by Na+ on the background of the smaller IKs of M-cells, acts to sufficiently prolong the AP plateau to allow for L-type Ca2+ channels recovery and reactivation resulting in EADs [26,31]. The morphology of the AP is profoundly disrupted as a result of this mutation.
|
Fig. 9 demonstrates the effects on the AP of HERG mutations resulting in decreased density of functional channels. We simulate the effect of the resulting current suppression on the morphology and duration of the AP in M-cells. In epicardial cells, these reductions act to prolong APD but do not give rise to the development of EADs (not shown). We consider 25, 50 and 75% reductions as corresponding to defects due to defective protein expression [32], impaired transport to the cell surface [28], or dominant negative suppression [9], respectively. A simulation showing the effects of reduced macroscopic current has been conducted previously in a Hodgkin—Huxely model of IKr [33]. Here we investigate the effects of a range of current reduction using the Markov IKr model in M-cells. APDs of the mutant M-cells are markedly affected over the range of reductions tested (Fig. 9). For a 75% reduction at slow rates (CL=1000 ms) the AP is marked by EADs (Fig. 9, arrow).
|
Fig. 10 (panels A and B) demonstrates the effects of changes in pacing rate in a cell with a 50% loss of functional IKr channels. In epicardial cells (panel A) the resulting reduction in current gives rise to a subtle prolongation of the APD that becomes more pronounced at slow rates. In panel B, the rate-dependent effect in an M-cell is more pronounced. At a slow CL=2000 ms, the M-cell APD is markedly prolonged and a characteristic shoulder-morphology AP is observed.
|
When IKr is absent as a consequence of a homozygous mutation (functional IKr knockout), the patient phenotype is marked by severe prolongation of the QT interval and polymorphic ventricular tachyarrhythmias [15]. In Fig. 10, panels C and D, simulation of this condition results in the most severe disruption of cellular repolarization. Epicardial cells (panel A) have prolonged APD (by 29 ms at CL of 750 ms; 35 ms at a CL of 2000 ms). The 999th and 1000th beats at the indicated cycle length are shown for an M-cell in panel D. The M-cell develops EADs even at a normal human physiological cycle length of 750 ms. During severe bradycardia (CL=2000 ms) multiple EADs develop.
| 4 Discussion |
|---|
|
|
|---|
In this study, functional defects in IKr due to several HERG mutations are introduced into a comprehensive model of the cardiac ventricular myocyte. The results demonstrate that distinct mutations can have variable effects on current morphology and lead to varying degrees of altered electrophysiology, depending on kinetic changes induced by the mutation. The use of a Markovian model to represent IKr deviates from the traditional Hodgkin—Huxely approach [34]. The Markovian scheme represents distinct channel states and coupling between these states, allowing us to relate state-specific kinetic properties of ion channels to the electrophysiological behavior of the whole cell.
Heterogeneous mutations in HERG that give rise to the same syndrome, congenital LQTS, may differ in severity and differentially affect the morphology of IKr during the AP. The T474I mutation predominantly affects the current at the onset of the AP through a leftward shift in the voltage dependence of activation. However, recovery from inactivation and deactivation are unaffected. Hence, the morphology of the late current is preserved in the T474I mutation, as evidenced by the preservation of the pronounced peak in channel open probability (P(O)) late in cellular repolarization seen in WT cells (Fig. 6, arrow). The result is a relatively small effect on APD.
In contrast, the R56Q mutation affects deactivation kinetics that decrease IKr late in the AP. The delicate equilibrium between recovery from inactivation and deactivation is altered as a result of this defect. The increased rate of deactivation acts to pull channels from the open state into closed states before they can conduct a substantial current. It is clear from examination of the current during the AP (Fig. 6, right), that late IKr is diminished due to the lack of the late maximum in channel P(O), which markedly slows AP repolarization and prolongs APD. Here, we focus on the most severe of the PAS mutations, but this behavior can be generalized to various defects that act to speed deactivation [30].
The N629D mutation gives rise to IKr that fails to inactivate. This gain of function would seemingly manifest as a larger IKr that would act to shorten APD. However, the disruption of the pore selectivity sequence that also results from the mutation allows for passage of Na+ ions into the cell. The loss of ion selectivity gives rise to a shift in the reversal potential of IKr to –13 mV. This elevated reversal potential is crossed during the plateau repolarization phase of the action potential. Once the membrane potential falls below –13 mV, the channel passes an inward current carried by Na+ (see Fig. 8, panel B, arrows). This inward current acts to depolarize the membrane late in the AP and prolongs APD. It results in a severe cellular phenotype similar to the effects of the
KPQ gain of function mutation described previously in the cardiac Na+ channel [19]. We also investigated the effects of loss of C-type inactivation and of loss of ion selectivity individually (not shown). Loss of inactivation caused a reduction in APD due to increased IKr while loss of ion selectivity was sufficient to cause APD prolongation.
To further explore the effects of HERG mutations, we have characterized the effects of a range of reductions in GKr on the morphology of the AP. A 25, 50, or 75% reduction in current could arise in a heterozygous mutation from reduced protein expression, processing or transport defect, or from dominant negative effects, respectively [15,28,32]. A 50% reduction in IKr acts to prolong the epicardial AP slightly, but has a major prolongation effect on the APD of M-cells especially at slow rates. A functional knockout (homozygous mutation) results in the most severe disruption of cellular repolarization, consistent with patient phenotypes observed clinically. The mutation acts to prolong APD at plateau potentials, providing sufficient time for reactivation of the L-type Ca2+ current, which depolarizes the cell and generates the EAD [26,31]. It should be emphasized that L-type reactivation during the prolonged plateau is common to all mutations investigated in this paper and provides the depolarizing charge during the EAD upstroke.
In general, structural and functional information is obtained from isolated membrane patches and/or cloned channels in expression systems (e.g. Xenopus oocytes) removed from the environment where they function physiologically. The natural cellular milieu is highly interactive and modulates single channel kinetic behavior. Cellular models are unique tools that allow for reinstitution of single channels into the cell in order to investigate integrated electrophysiological behavior. That the physiological effect of a given channel defect depends on the cellular environment is demonstrated here by the greater severity of APD prolongation and disruption of repolarization in M cells, suggesting that such mutations can increase APD differences and dispersion of repolarization in the myocardium [33]. A substrate with increased dispersion of repolarization is susceptible to the development of unidirectional conduction block and reentrant arrhythmias [2,35]. Such arrhythmias can be triggered by the EADs associated with the mutations which provide premature stimuli to a highly susceptible substrate [36].
In the context of arrhythmia, it is interesting to note that arrhythmias associated with HERG mutations and LQT2 are not rate dependent to the same extent as other forms of inherited LQTS [37]. Our simulations predict that HERG defects can lead to the development of cellular arrhythmogenic responses (APD prolongation and early afterdepolarizations) at slow rates in the absence of neural influences and may explain the development of arrhythmias in patients during sleep and relaxation. However, HERG defects have also been linked to arrhythmias in patients startled by loud noise, such as a telephone or alarm clock, suggesting a role for the autonomic nervous system in increasing susceptibility to arrhythmic events [38,39]. Previous modeling studies in our laboratory have incorporated the effects of β-adrenergic stimulation in a model of the WT cell [31]. Future efforts will expand the previous studies to investigate effects of the autonomic nervous system in the context of channelopathies including LQTS.
Time for primary review 27 days.
| Appendix. IKr formulation |
|---|
|
|
|---|
![]() |
Wild-type rate constants
![]() |
T474I rates
|
|
R56Q rates
|
|
N629D rates
|
|
N629D loss of selectivity
![]() |
| Acknowledgements |
|---|
This work was supported by the NIH National Heart, Lung and Blood Institute grants R01-HL49054 and R37-HL33343 (to YR) and a Whitaker Foundation Development Award.
| References |
|---|
|
|
|---|
- Keating M.T., Sanguinetti M.C. Molecular genetic insights into cardiovascular disease. Science (1996) 272:681–685.[Abstract]
- Schwartz P.J., Priori S., Napolitano C. Cardiac electrophysiology: from cell to bedside. Zipes D.P., Jalife J., eds. (2000) 3rd ed. Philidelphia: W.B. Saunders. 597–615.
- Roden D.M., Balser J.R. A plethora of mechanisms in the HERG-related long-QT syndrome: Genetics meets electrophysiology. Cardiovasc Res (1999) 44:242–246.
[Free Full Text] - Curran M.E., Splawski I., Timothy K.W., Vincent G.M., Green E.D., Keating M.T. A molecular basis for cardiac arrhythmia: HERG mutations cause long QT syndrome. Cell (1995) 80:795–803.[CrossRef][Web of Science][Medline]
- Sanguinetti M.C., Jurkiewicz N.K. Two components of cardiac delayed rectifier K+ current. Differential sensitivity to block by class III antiarrhythmic agents. J Gen Physiol (1990) 96:195–215.
[Abstract/Free Full Text] - Shibasaki T. Conductance and kinetics of delayed rectifier potassium channels in nodal cells of the rabbit heart. J Physiol (1987) 387:227–250.
[Abstract/Free Full Text] - Abbott G.W., Sesti F., Splawski I., et al. MiRP1 forms IKr potassium channels with HERG and is associated with cardiac arrhythmia. Cell (1999) 97:175–187.[CrossRef][Web of Science][Medline]
- Petrecca K., Atanasiu R., Akhavan A., Shrier A. N-linked glycosylation sites determine HERG channel surface membrane expression. J Physiol (1999) 515:41–48.
[Abstract/Free Full Text] - Sanguinetti M.C., Curran M.E., Spector P.S., Keating M.T. Spectrum of HERG K+ channel dysfunction in an inherited cardiac arrhythmia. Proc Natl Acad Sci USA (1996) 93:2208–2212.
[Abstract/Free Full Text] - Ficker E., Jarolimek W., Kiehn J., Baumann A., Brown A.M. Molecular determinants of dofetilide block of HERG K+ channels. Circ Res (1998) 82:386–395.
[Abstract/Free Full Text] - Smith P.L., Baukrowitz T., Yellen G. The inward rectification mechanism of the HERG cardiac potassium channel. Nature (1996) 379:833–836.[CrossRef][Web of Science][Medline]
- Sanguinetti M.C., Xu Q.P. Mutations of the S4-S linker alter activation properties of HERG potassium channels expressed in Xenopus oocytes. J Physiol (1999) 514:667–675.
[Abstract/Free Full Text] - Cabral J.H.M., Lee A., Cohen S.L., Chait B.T., Li M., Mackinnon R. Crystal structure and functional analysis of the HERG potassium channel N terminus: a eukaryotic PAS domain. Cell (1998) 95:649–656.[CrossRef][Web of Science][Medline]
- Lees-Miller J.P., Duan Y., Teng G.Q., Thorstad K., Duff H.J. Novel gain of function mechanism in K+ channel related Long-QT Syndrome: altered gating and selectivity in the HERG1 N629D mutant. Circ Res (2000) 86:507–513.
[Abstract/Free Full Text] - Hoorntje T., Alders M., van Tintelen P., et al. Homozygous premature truncation of the HERG protein: the human HERG knockout. Circulation (1999) 100:1264–1267.
[Abstract/Free Full Text] - Luo C.H., Rudy Y. A dynamic model of the cardiac ventricular action potential. I. Simulations of ionic currents and concentration changes. Circ Res (1994) 74:1071–1096.
[Abstract/Free Full Text] - Zeng J., Laurita K.R., Rosenbaum D.S., Rudy Y. Two components of the delayed rectifier K+ current in ventricular myocytes of the guinea pig type. Theoretical formulation and their role in repolarization. Circ Res (1995) 77:140–152.
[Abstract/Free Full Text] - Viswanathan P.C., Shaw R.M., Rudy Y. Effects of IKr and IKs heterogeneity on action potential duration and its rate dependence. Circulation (1999) 99:2466–2474.
[Abstract/Free Full Text] - Clancy C.E., Rudy Y. Linking a genetic defect to its cellular phenotype in a cardiac arrhythmia. Nature (1999) 400:566–569.[CrossRef][Medline]
- Balser J.R., Roden D.M., Bennett P.B. Global parameter optimization for cardiac potassium channel gating models. Biophys J (1990) 57:433–444.[Web of Science][Medline]
- Wang S., Liu S., Morales M.J., Strauss H.C., Rasmusson R.L. A quantitative analysis of the activation and inactivation kinetics of HERG expressed in Xenopus oocytes. J Physiol (1997) 502:45–60.
[Abstract/Free Full Text] - Zhou Z., Gong Q., Ye B., et al. Properties of HERG channels stably expressed in HEK 293 cells studied at physiological temperature. Biophys J (1998) 74:230–241.[Web of Science][Medline]
- Hancox J.C., Levi A.J., Witchel H.J. Time course and voltage dependence of expressed HERG current compared with the native rapid delayed rectifier K current during the cardiac ventricular action potential. Pflüger's Arch (1998) 436:843–853.[CrossRef][Web of Science][Medline]
- Kiehn J., Lacerda A.E., Brown A.M. Pathways of HERG inactivation. Am J Physiol (1999) 277:H199–H210.[Web of Science][Medline]
- Yang T., Snyders D.J., Roden D.M. Rapid inactivation determines the rectification and [K+]o dependence of the rapid component of the delayed rectifier K+ current in cardiac cells. Circ Res (1997) 80:782–789.
[Abstract/Free Full Text] - Viswanathan P.C., Rudy Y. Pause induced early afterdepolarizations in the long QT syndrome: a simulation study. Cardiovasc Res (1999) 42:530–542.
[Abstract/Free Full Text] - Faber GM, Rudy Y. Action potential and contractility changes in [Na+] overloaded cardiac myocytes: A simulation study. Biophys J 2000 (in press).
- Zhou Z., Gong Q., Epstein M.L., January C.T. HERG channel dysfunction in human long QT syndrome. Intracellular transport and functional defects. J Biol Chem (1998) 273:21061–21066.
[Abstract/Free Full Text] - Nakajima T., Furukawa T., Tanaka T., et al. Novel mechanism of HERG current suppression in LQT2: shift in voltage dependence of HERG inactivation. Circ Res (1998) 83:415–422.
[Abstract/Free Full Text] - Chen J., Zou A., Splawski I., Keating M.T., Sanguinetti M.C. Long QT syndrome-associated mutations in the Per—Arnt–Sim (PAS) domain of HERG potassium channels accelerate channel deactivation. J Biol Chem (1999) 274:10113–10118.
[Abstract/Free Full Text] - Zeng J., Rudy Y. Early afterdepolarizations in cardiac myocytes: mechanism and rate dependence. Biophys J (1995) 68:949–964.[Web of Science][Medline]
- Li X., Xu J., Li M. The human delta1261 mutation of the HERG potassium channel results in a truncated protein that contains a subunit interaction domain and decreases the channel expression. J Biol Chem (1997) 272:705–708.
[Abstract/Free Full Text] - Viswanathan P., Rudy Y. Cellular arrhythmogenic effects of congenital and acquired Long-QT syndrome in the heterogeneous Myocardium. Circulation (2000) 101:1192–1198.
[Abstract/Free Full Text] - Hodgkin A.L., Huxely A.F. A quantitative description of membrane current and its application to conduction and excitation in nerve. J Physiol (1952) 117:500–544.
[Free Full Text] - Antzelevitch C., Yan G., Shimuzu W., Burashnikov P. Cardiac electrophysiology: from cell to bedside. Zipes D.P., Jalife J., eds. (2000) Philidelphia: W.B. Saunders. 222–238.
- El-Sherif N., Chinushi M., Caref E.B., Restivo M. Electrophysiological mechanism of the characteristic electrocardiographic morphology of torsade de pointes tachyarrhythmias in the long-QT syndrome: detailed analysis of ventricular tridimensional activation patterns. Circulation (1997) 96:4392–4399.
[Abstract/Free Full Text] - Schwartz P.J., Priori S.G., Locati E.H., et al. Long QT syndrome patients with mutations of the SCN5A and HERG genes have differential responses to Na+ channel blockade and to increases in heart rate. Implications for gene-specific therapy. Circulation (1995) 92:3381–3386.
[Abstract/Free Full Text] - Wilde A.A.M., Jongbloed R.J.E., Doevendans P.A., et al. Auditory stimuli as a trigger for arrhythmic events differentiate HERG-related (LQTS2) patients from KVLQT1-related patients (LQTS1). J Am Coll Cardiol (1999) 33:327–332.
[Abstract/Free Full Text] - Moss A.J., Robinson J.L., Gessman L., et al. Comparison of clinical and genetic variables of cardiac events associated with loud noise versus swimming among subjects with the Long QT syndrome. Am J Cardiol (1999) 84:876–879.[CrossRef][Web of Science][Medline]
This article has been cited by other articles:
![]() |
H. Itoh, T. Sakaguchi, W.-G. Ding, E. Watanabe, I. Watanabe, Y. Nishio, T. Makiyama, S. Ohno, M. Akao, Y. Higashi, et al. Latent Genetic Backgrounds and Molecular Pathogenesis in Drug-Induced Long-QT Syndrome Circ Arrhythm Electrophysiol, October 1, 2009; 2(5): 511 - 523. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. S. Gustina and M. C. Trudeau A recombinant N-terminal domain fully restores deactivation gating in N-truncated and long QT syndrome mutant hERG potassium channels PNAS, August 4, 2009; 106(31): 13082 - 13087. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Fink and D. Noble Markov models for ion channels: versatility versus identifiability and speed Phil Trans R Soc A, June 13, 2009; 367(1896): 2161 - 2179. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. Sims, S. Reisenweber, P. C. Viswanathan, B.-R. Choi, W. H. Walker, and G. Salama Sex, Age, and Regional Differences in L-Type Calcium Current Are Important Determinants of Arrhythmia Phenotype in Rabbit Hearts With Drug-Induced Long QT Type 2 Circ. Res., May 9, 2008; 102(9): e86 - e100. [Abstract] [Full Text] [PDF] |
||||
![]() |
Y. Rudy Modelling the molecular basis of cardiac repolarization Europace, November 1, 2007; 9(suppl_6): vi17 - vi19. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. S. Hansen, S.-P. Olesen, and M. Grunnet Pharmacological Activation of Rapid Delayed Rectifier Potassium Current Suppresses Bradycardia-Induced Triggered Activity in the Isolated Guinea Pig Heart J. Pharmacol. Exp. Ther., June 1, 2007; 321(3): 996 - 1002. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. B. Huffaker, J. N. Weiss, and B. Kogan Effects of early afterdepolarizations on reentry in cardiac tissue: a simulation study Am J Physiol Heart Circ Physiol, June 1, 2007; 292(6): H3089 - H3102. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. Vecchietti, E. Grandi, S. Severi, I. Rivolta, C. Napolitano, S. G. Priori, and S. Cavalcanti In silico assessment of Y1795C and Y1795H SCN5A mutations: implication for inherited arrhythmogenic syndromes Am J Physiol Heart Circ Physiol, January 1, 2007; 292(1): H56 - H65. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. E. Clancy, Z. I. Zhu, and Y. Rudy Pharmacogenetics and anti-arrhythmic drug therapy: a theoretical investigation Am J Physiol Heart Circ Physiol, January 1, 2007; 292(1): H66 - H75. [Abstract] [Full Text] [PDF] |
||||
![]() |
C.-u. Choe, E. Schulze-Bahr, A. Neu, J. Xu, Z. I. Zhu, K. Sauter, R. Bahring, S. Priori, P. Guicheney, G. Monnig, et al. C-terminal HERG (LQT2) mutations disrupt IKr channel regulation through 14-3-3{epsilon} Hum. Mol. Genet., October 1, 2006; 15(19): 2888 - 2902. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Fink, W. R Giles, and D. Noble Contributions of inwardly rectifying K+ currents to repolarization assessed using mathematical models of human ventricular myocytes Phil Trans R Soc A, May 15, 2006; 364(1842): 1207 - 1222. [Abstract] [Full Text] [PDF] |
||||
![]() |
H. Sun, X. Liu, Q. Xiong, S. Shikano, and M. Li Chronic Inhibition of Cardiac Kir2.1 and hERG Potassium Channels by Celastrol with Dual Effects on Both Ion Conductivity and Protein Trafficking J. Biol. Chem., March 3, 2006; 281(9): 5877 - 5884. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. E. Clancy and R. S. Kass Inherited and Acquired Vulnerability to Ventricular Arrhythmias: Cardiac Na+ and K+ Channels Physiol Rev, January 1, 2005; 85(1): 33 - 47. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. Royer, S. Demolombe, A. El Harchi, K. Le Quang, J. Piron, G. Toumaniantz, D. Mazurais, C. Bellocq, G. Lande, C. Terrenoire, et al. Expression of human ERG K+ channels in the mouse heart exerts anti-arrhythmic activity Cardiovasc Res, January 1, 2005; 65(1): 128 - 137. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. J. Saucerman, S. N. Healy, M. E. Belik, J. L. Puglisi, and A. D. McCulloch Proarrhythmic Consequences of a KCNQ1 AKAP-Binding Domain Mutation: Computational Models of Whole Cells and Heterogeneous Tissue Circ. Res., December 10, 2004; 95(12): 1216 - 1224. [Abstract] [Full Text] [PDF] |
||||
![]() |
V. E. Bondarenko, G. P. Szigeti, G. C. L. Bett, S.-J. Kim, and R. L. Rasmusson Computer model of action potential of mouse ventricular myocytes Am J Physiol Heart Circ Physiol, September 1, 2004; 287(3): H1378 - H1403. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. G. KLEBER and Y. RUDY Basic Mechanisms of Cardiac Impulse Propagation and Associated Arrhythmias Physiol Rev, April 1, 2004; 84(2): 431 - 488. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Firouzi and W. A. Groenewegen Gene polymorphisms and cardiac arrhythmias Europace, January 1, 2003; 5(3): 235 - 242. [Full Text] [PDF] |
||||
![]() |
Y. Kurata, I. Hisatome, S. Imanishi, and T. Shibamoto Dynamical description of sinoatrial node pacemaking: improved mathematical model for primary pacemaker cell Am J Physiol Heart Circ Physiol, November 1, 2002; 283(5): H2074 - H2101. [Abstract] [Full Text] [PDF] |
||||
![]() |
E. Pietila, H. Fodstad, E. Niskasaari, P.a. J. Laitinen, H. Swan, M. Savolainen, Y. A. Kesaniemi, K. Kontula, and H. V. Huikuri association between HERG K897T polymorphism and QT interval in middle-aged finnish women J. Am. Coll. Cardiol., August 7, 2002; 40(3): 511 - 514. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. Gima and Y. Rudy Ionic Current Basis of Electrocardiographic Waveforms: A Model Study Circ. Res., May 3, 2002; 90(8): 889 - 896. [Abstract] [Full Text] [PDF] |
||||
![]() |
Y. Lu, M. P Mahaut-Smith, A. Varghese, C. L-H Huang, P. R Kemp, and J. I Vandenberg Effects of premature stimulation on HERG K+ channels J. Physiol., December 15, 2001; 537(3): 843 - 851. [Abstract] [Full Text] [PDF] |
||||
![]() |
Members of the Sicilian Gambit New Approaches to Antiarrhythmic Therapy, Part I: Emerging Therapeutic Applications of the Cell Biology of Cardiac Arrhythmias Circulation, December 4, 2001; 104(23): 2865 - 2873. [Abstract] [Full Text] [PDF] |
||||
![]() |
Members of the Sicilian Gambit New approaches to antiarrhythmic therapy; emerging therapeutic applications of the cell biology of cardiac arrhythmias Eur. Heart J., December 1, 2001; 22(23): 2148 - 2163. [Abstract] [PDF] |
||||
![]() |
Members of the Sicilian Gambit New approaches to antiarrhythmic therapy: emerging therapeutic applications of the cell biology of cardiac arrhythmias Cardiovasc Res, December 1, 2001; 52(3): 345 - 360. [Abstract] [Full Text] [PDF] |
||||
![]() |
Y. Rudy Multiple interactions determine cellular electrical processes in the multicellular tissue Cardiovasc Res, July 1, 2001; 51(1): 1 - 3. [Full Text] [PDF] |
||||
![]() |
P. M Spooner, S. G Priori, and R. J Myerburg Spotlight on sudden cardiac death Cardiovasc Res, May 1, 2001; 50(2): 173 - 176. [Full Text] [PDF] |
||||
![]() |
C. E. Clancy and Y. Rudy Na+ Channel Mutation That Causes Both Brugada and Long-QT Syndrome Phenotypes: A Simulation Study of Mechanism Circulation, March 12, 2002; 105(10): 1208 - 1213. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. Gima and Y. Rudy Ionic Current Basis of Electrocardiographic Waveforms: A Model Study Circ. Res., May 3, 2002; 90(8): 889 - 896. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||



























