MolPharm

Home Help [Feedback] [For Subscribers] [Archive] [Search] [Contents]
 QUICK SEARCH:   [advanced]


     


Molecular Pharmacology Fast Forward
First published on March 7, 2008; DOI: 10.1124/mol.107.042085


0026-895X/08/7306-1643-1651$20.00
Mol Pharmacol 73:1643-1651, 2008

This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Data Supplement
Right arrow All Versions of this Article:
mol.107.042085v1
73/6/1643    most recent
Right arrow Submit a response
Right arrow Alert me when this article is cited
Right arrow Alert me when eLetters are posted
Right arrow Alert me if a correction is posted
Services
Right arrow Related articles in MolPharm
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Myokai, T.
Right arrow Articles by Oiki, S.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Myokai, T.
Right arrow Articles by Oiki, S.

Topological Mapping of the Asymmetric Drug Binding to the Human Ether-à-go-go-Related Gene Product (HERG) Potassium Channel by Use of Tandem DimersFormula

Toshihiko Myokai, Sunghi Ryu, Hirofumi Shimizu, and Shigetoshi Oiki

Department of Molecular Physiology and Biophysics, Faculty of Medical Sciences, University of Fukui, Fukui, Japan (T.M., H.S., S.O.); and Pharmaceutical Research Division, Takeda Pharmaceutical Company Limited, Osaka, Japan (T.M., S.R.)

Received for publication September 25, 2007.

Accepted for publication March 7, 2008.


    Abstract
 Top
 Abstract
 Materials and Methods
 Results
 Discussion
 Conclusion
 References
 
The human ether-à-go-go related gene product (HERG) channel is essential for electrical activity of heart cells, and block of this channel by many drugs leads to lethal arrhythmias. Tyr652 and Phe656 of the sixth transmembrane helix are candidates for the drug binding site. In the tetrameric HERG channel, a drug with asymmetric structure should interact unevenly with multiple residues from different subunits. To elucidate the topology of the drug-binding site, we constructed tandem dimers of HERG channels and the aromatic Tyr652 and Phe656 residues were replaced by alanine singly or doubly. Eight types of HERG channels, including homotetrameric mutants, having different numbers and arrangements of aromatic residues at the blocking site, were studied. Effects of cisapride on channels expressed in Xenopus laevis oocytes were examined electrophysiologically. The inhibition constants (Ki) were increased significantly as the diagonal Tyr652 were deleted, whereas those for the diagonal Phe656-deleted mutant were not changed. These results suggest that Tyr652 residues from adjacent subunits contributed to the binding. Two types of double mutants of tandem dimers showed significantly distinct affinities, suggesting that the coexistence of Tyr652 and Phe656 on a subunit in diagonal position is crucial to having a high affinity. Thermodynamic double-mutant cycle analyses revealed interactions between Tyr652 and Phe656 upon binding. The kinetics and voltage-dependence of blocking suggested transitions of the binding site from low to high affinity. These approaches using a set of mutant HERG channels gave a dynamic picture of the spatial arrangements of residues that contribute to the drug-channel interaction.


Malfunction of the human ether-à-go-go related gene product (HERG) potassium channels by either hereditary disorders or drug block leads to the long QT syndrome, which increases the risk of lethal arrhythmias (Ashcroft, 2000Go; Haverkamp et al., 2000Go; Sanguinetti and Tristani-Firouzi, 2006Go). A unique property of the HERG channel is that it is blocked by a variety of drugs (Haverkamp et al., 2000Go). Understanding the mechanism of the susceptibility to the wide spectra of molecular species is of primary importance in developing drugs without side effects to HERG channels (Mitcheson et al., 2000Go). Experimental data have accumulated about the residues that contribute to the drug binding, and the most important residues are Tyr652 and Phe656 in the sixth transmembrane helix (Mitcheson et al., 2000Go, 2005Go). The HERG channel is composed of four identical subunits that are assembled into a 4-fold symmetric structure (Morais Cabral et al., 1998Go); however, the blocking drugs are asymmetric in their structure and large enough to span multiple subunits. Thus, aromatic residues (Tyr652 and Phe656) from different subunits could contribute unequally to the binding during drug-channel interaction. To identify the contributing residues and their positions in a tetrameric channel, a logical and powerful strategy is to reduce the symmetry of the HERG channel from 4-fold into 2-fold. We constructed, for the first time, a tandem dimer of the HERG channel and introduced aromatic ring-deleting mutations (F656A and Y652A). Here, we focused on a typical HERG channel blocker, cisapride (Mohammad et al., 1997Go; Walker et al., 1999Go; Mitcheson et al., 2000Go), as a probe for the drug-binding sites (Mitcheson et al., 2000Go). In tandem dimers, double mutations generated different positional arrangements of mutated sites in a channel. A set of mutants, including homotetrameric mutants, permitted a systematic approach to elucidate the topological map of the cisapride binding in HERG channel.


    Materials and Methods
 Top
 Abstract
 Materials and Methods
 Results
 Discussion
 Conclusion
 References
 
Construction of Tandem Dimers. The wild-type HERG cDNA was cloned at Takeda Pharmaceutical Company Limited (Osaka, Japan), and subcloned into the NheI and EcoRI sites in pCDNA3.1(+) plasmid expression vector (Invitrogen, Carlsbad, CA). Single or double mutations (Y652A and/or F656A) were generated by the overlap extension polymerase chain reaction using Platinum Taq DNA Polymerase High Fidelity (Invitrogen, Carlsbad, CA) and confirmed by sequencing the mutated region and by restriction enzyme analysis. The homotetrameric mutant harboring Y652A or F656A was made by introducing the mutated fragment into the HERG channel gene in the pcDNA3.1(+) vector. The HERG tandem dimers were constructed by linking two monomers that were prepared separately (Christie et al., 1990Go; Isacoff et al., 1990Go; Ruppersberg et al., 1990Go; Liman, 1992Go). For the N-terminal side monomer, the wild-type or mutant (Y652A) HERG channel gene in pcDNA3.1(+) was modified by deleting its stop codon. For the C-terminal side monomer, the HERG channel gene (wild type, harboring single or double mutations) was subcloned downstream of the in-frame EcoRI site (note that the original NheI site in front of the first methionine of the channel gene was modified to a unique EcoRI site when subcloned) in the pKF3 vector (Takara Bio, Kyoto, Japan). To generate tandem dimers, the entire channel coding region in the pKF3 vector was digested with EcoRI and ligated into the pcDNA3.1(+) vector harboring the stop codon-deleted channel gene. The two monomers were linked by the nucleotide sequence of GAATTC. For expression in Xenopus laevis oocytes, a stretch of 30 polyA (dA:dT) residues was inserted behind the stop codon.

Western Blot Analysis. HEK293 cells stably expressing HERG channels were homogenized at 4°C in lysis buffer (50 mM Tris-HCl, pH 7.5, 150 mM NaCl, and 1 mM EDTA), containing a protease inhibitor mix (Roche Applied Science, Indianapolis, IN) and spun at 500g for 10 min. Pellets of the membrane fractions were produced from the low-speed supernatants by centrifugation at 150,000 rpm for 30 min. For Western blots, membrane proteins were separated on 7% SDS-polyacrylamide gels and transferred onto polyvinylidene difluoride membranes. Membranes were probed with polyclonal rabbit anti-HERG antibody (Millipore Bioscience Research Reagents, Temecula, CA) and the antibody was detected with an ECL detection kit (GE Healthcare, Chalfont St. Giles, Buckinghamshire, UK).

Electrophysiological Recordings. Two-electrode voltage-clamp experiments were performed. Procedures for the isolation of X. laevis oocytes, injection of cRNA and electrophysiological experiments have been described in detail previously (Shimizu et al., 2003Go) and were approved by the Animal Research Committee of University of Fukui. The bath solution contained 4 mM KCl, 96 mM N-methyl-D-glucamine Cl, 1.8 mM CaCl2, 1 mM MgCl2, and 5 mM HEPES, pH 7.40. Cisapride (Accurate Chemicals, Westbury, NY) was dissolved into dimethyl sulfoxide to make 10 mM stock solutions and kept at -20°C until use. Voltage clamp experiments were performed on the oocytes 2 to 6 days after injection using a Dagan CA1 amplifier (Dagan Corporation, Minneapolis, IL) at 22-24°C. For the recordings from oocytes expressing the F656A mutant, which has been known as poor expression of channels (Mitcheson et al., 2000Go; Milnes et al., 2003Go), electrophysiological experiments were performed within 3 days after the injection until the oocyte became damaged. The oocytes were perfused continually with the perfusates at the rate of 1 ml/min using an infusion pump (TE-331; TERUMO, Tokyo, Japan). Oocytes exhibiting endogenous currents were excluded from analyses.

Activation and Inactivation Gate. To characterize the activation gate for mutant channels, the peak amplitudes for the tail currents at -50 mV were plotted as a function of the preceding activation voltages. This tail I-V curve was fitted with the Boltzmann function.

Formula(1)
where V is the membrane potential, V1/2 is the half activation voltage, dx is the slope factor, and A1 and A2 are the amplitude factors. For the inactivation gate, the current amplitudes at the end of the 3-s depolarization pulses (isochronal I-V curves) were drawn, which were superimposed on the tail I-V curves for each channel species (Fig. S1). The isochronal I-V curves show a typical bell shape, indicating the characteristic inactivation of HERG channel. The double-Boltzmann function was used for the fitting.

Formula(2)
where Vinac and dxinac values were the half inactivation voltage and the slope factor. A1 and A2 are the amplitude factors. In this fitting, the fitted parameters from the activation data (V1/2 and dx) were used.

Blocking Inhibition Constants. Block of the HERG currents at each cisapride concentration was evaluated after the steady state of blocking was reached. This was performed by eliciting repeatedly 3.5-s depolarization pulses to +20 mV followed by repolarization at -50 mV an interval of 20 s (the holding potential was -80 mV) and checking the current amplitudes reached at the steady state. At low concentrations, it took 10 to 20 minutes. Inhibition constants (Ki) were obtained from the peak amplitudes of the tail currents at -50 mV preceded by the depolarizing pulse to +20 mV for 3 s. From concentration dependence of current amplitude, Ki values and the Hill coefficients (n) were obtained from the concentration dependence of current amplitude by fits to the equation:

Formula(3)

Voltage-dependence of the binding affinity was obtained by fitting the Ki values with the function.

Formula(4)
where K0 represents the inhibition constant at 0 mV, z is the gating charge, and e, k, and T are the elementary charge, the Boltzmann constant, and the absolute temperature, respectively.

Thermodynamic Double-Mutant Cycle. To estimate the interaction free energy upon binding, thermodynamic double-mutant cycle analysis was applied (Schreiber and Fersht, 1995Go; Ranganathan et al., 1996Go; Fersht, 1999Go; Yan et al., 2006Go). First, the binding free energy (G) was evaluated from the inhibition constant.

Formula(5)

The differences in the affinities between two channel species were obtained from the differences in G. For example, the change in the binding free energy by a single mutation relative to the wild type (WT) is

Formula(6)

A double mutant can be regarded as a combination of two single mutations. Therefore, {Delta}G for the double mutant compared with WT is:

Formula(7)
The {Delta}G for mutation of Phe656 site introduced into either WT or Y652A is {Delta}GsingleM-WT or {Delta}GdoubleM-singleM. The difference in these energies defines the coupling energy, {Delta}{Delta}Gint, which represents interaction between two sites.

Formula(8)
The values of {Delta}{Delta}Gint indicate cooperativity between residues.

Blocking Kinetics. The apparent blocking rate (equal to the reciprocal of the time constant) for a simple three-state model can be expressed as the following equation, if the binding and unbinding kinetics are fast relative to the subsequent transition represented by the rate constants of k3 and koff.

Formula(9)
Kf (= k2/k1) is the dissociation constant for the initial binding site (see Supplemental Data).


Figure 1
View larger version (37K):
[in this window]
[in a new window]

 
Fig. 1. Mutants of HERG channels. A, a schematic diagram of the tandem dimer construct of the HERG gene. B, Western blot analyses of homotetrameric and tandem dimeric HERG gene products.

 
Data Analysis. Experimental data were analyzed with Origin software (OriginLab Corp., Northampton, MA) and Mathematica (Wolfram Research Inc., Champaign, IL). Data are presented as mean ± S.E.M.


    Results
 Top
 Abstract
 Materials and Methods
 Results
 Discussion
 Conclusion
 References
 
Expression of Tandem Dimers. We constructed seven mutants (Table 1), which were expressed in X. laevis oocytes and human embryonic kidney 293 cells (Ando et al., 2005Go). Western blot was performed on human embryonic kidney 293 cells to verify the expression of the mutants (Fig. 1B). Tandem dimers showed protein bands at the expected locations. For monomeric channels, double protein bands were observed at around 135/155 kDa, as has been reported (Zhou et al., 1998Go; Ficker et al., 2004Go). For tandem dimers, double bands were seen at 270/310 kDa.


View this table:
[in this window]
[in a new window]

 
TABLE 1 Mutant constructs of HERG channels

In this study, the mutants were named in the conventional manner for homotetrameric channels (such as F656A). For tandem dimers, single or double mutations introduced into the N- or C-terminal subunits separated by a colon for the linker.

 

Electrophysiological Characterization of Mutant Channels. The electrophysiological properties of mutant and WT HERG channels expressed in X. laevis oocytes were examined. Gating properties for WT and all the mutants, including the tandem WT (td[wt:wt]), were evaluated from the current traces elicited by depolarizing pulses (Fig. 2A). The activation and inactivation gating were measured by plotting the tail I-V curves and the isochronal I-V curves at the end of the depolarization pulses (Fig. 2A, inset, and Supplemental Fig. S1). The voltage dependence of activation and inactivation gating for WT and td[wt:wt] were basically similar, suggesting that concatenation did not affect the steady-state gating. In all mutants, the activation curves were shifted slightly in the negative direction (Table 2), although the slope factor was not changed significantly. Deactivation for td[wt:wt] was faster compared with the WT. This tendency seemed to hold for all the tandem dimers. A contribution of the N-terminal domain to the deactivation gating of HERG channel has been reported (Morais Cabral et al., 1998Go; Wang et al., 2000Go). Then a cytoplasmic domain between subunits in the tandem dimers might be restricted by the short linker, which may accelerate the deactivation rate.


Figure 2
View larger version (28K):
[in this window]
[in a new window]

 
Fig. 2. Representative current traces for WT and mutant HERG channels. Currents were elicited by 3-s step pulses from -100 to +40 mV with 10-mV increments followed by 3-s pulses to -50 mV (the holding potential was -80 mV). In this figure, the current amplitudes were normalized to emphasize the differences among the various channels. A, current traces for HERG channels in the absence of cisapride. Cartoons represent the arrangements of Tyr652 and Phe656 residues around the symmetrical axis of the channel. The partial overlap of the symbols for Tyr652 and Phe656 show that they are on the same subunit. Small circles represent mutations to alanine. Schemes with linkers represent tandem dimers. For tandem dimers, N- and C-terminal subunits were arranged counter-clockwise. Inset, the tail I-V curves (open symbols) and the isochronal I-V curves (closed symbols) at the end of the depolarizing pulses for WT (black) and td[wt:wt] (red). B, current traces in the presence of cisapride. Cisapride concentrations (micromolar) were 0.63 for WT, 0.25 for td[wt:wt], 0.63 for td[wt:F656A], 0.63 for td[wt:Y652AF656A], 1.6 for F656A, 0.63 for td[wt:Y652A], 1.6 for td[Y652A:F656A], and 10 for Y652A. Inset, the voltage command.

 

View this table:
[in this window]
[in a new window]

 
TABLE 2 Gating properties of HERG channels

The half activation and inactivation voltages (V1/2) and their slopes are shown. N represents the number of observations.

 

The Steady-State Blocking of Cisapride. To examine the concentration dependence of current blocking to cisapride, an oocyte was perfused continually throughout the experiment, and HERG currents were evaluated after the block reached the steady-level. The current traces showed (Fig. 2B) a slow decay during depolarizing pulses for WT and some of the other mutants and mostly reached the steady state at the end of the depolarization pulses. This decay was not observed for Y652A, F656A, and td[Y652A:F656A]. Ki values were obtained from the concentration dependence of the peak amplitudes of the tail currents at -50 mV (Table 3). The Hill coefficients were nearly one, suggesting the binding stoichiometry of the drug-HERG channel as one to one.


View this table:
[in this window]
[in a new window]

 
TABLE 3 Blocking parameters of cisapride to HERG channels at +20 mV

These Ki values represent drug binding at depolarized potential of +20 mV, since unblocking was slow and could not be detected at the time of peak-current measurements for most of the channel species except for F656A. The Ki value for F656A at +20 mV was obtained by correcting for the unblocking component at the tail potential (see Supplemental Data). N represents the number of observations

 

The free energy of the binding of cisapride to HERG channels, which was calculated from the Ki values as G = kT log Ki, is shown in Fig. 3. It is clearly seen that not only the number of the aromatic residues but also the spatial location of each residue is crucial for the binding affinity. The affinities were reduced dramatically by complete deletion of either of the aromatic residues (F656A and Y652A), confirming previous reports (Mitcheson et al., 2000Go). The levels of the binding energy for WT and td[wt:wt] were similar, indicating that the concatenation did not affect the binding of the blocker. The affinity for td[wt:F656A] was similar to that of WT, which indicates that, although Phe656 contributed to the binding significantly, at most two Phe656 residues from diagonal subunits were required for binding. On the other hand, the affinity of td[wt:Y652A] was significantly reduced, suggesting that Tyr652 residues from adjacent subunits contributed to the binding. This is the first experimental observation that cisapride may interact with multiple Tyr652 residues.


Figure 3
View larger version (28K):
[in this window]
[in a new window]

 
Fig. 3. The binding energetics. The vertical axis represents the free energy of binding. The binding energy for each mutant is shown as a level bar. Each column represents channel species having the same numbers of aromatic residues.

 
The sequence of the binding energies among mutants with the same number of aromatic residues was most informative (see the rightmost column of Fig. 3). First, to have both Phe656 and Tyr652 aromatic residues on a channel is more important for high-affinity binding than to have single-species aromatic rings. Second, the 4-fold higher affinity (the difference of the binding free energy: {Delta}G = 3.5 kJ/mol) of td[wt:Y652AF656A] than td[Y652A:F656A] indicates that having both Tyr652 and Phe656 on the same subunit is crucial for a channel to show a high affinity. Energetic considerations provided further clues on the contributions of residues to the binding and interactions between residues upon binding, which will be discussed below.

Voltage-Dependent Blocking. Voltage-dependent blocking was evaluated from current traces elicited by depolarizing pulses (Fig. 2B). In WT and some mutants, the slow decay of the currents seen at depolarized potentials was accelerated as the membrane potential was depolarized more. In the presence of cisapride, the tail I-V curves (the peak amplitudes of the tail currents as a function of the preceding depolarizing voltages) showed a maximum at around 0 mV for WT (Fig. 4A). The current amplitudes were depressed slightly as the membrane potential was depolarized further. This pattern is an indication of the presence of voltage-dependent block (Walker et al., 1999Go; Sánchez-Chapula et al., 2003Go). At each voltage, the Ki values were calculated from the concentration-dependence curves (Fig. 4B).


Figure 4
View larger version (22K):
[in this window]
[in a new window]

 
Fig. 4. Voltage-dependent block. A, the tail I-V curves of WT in the absence and in the presence of increased concentrations of cisapride. B, voltage-dependence of Ki values for WT and mutant HERG channels. The scale of the vertical axis is logarithmic. Symbols indicate: {circ}, WT; {blacktriangledown}, td[wt:F656A]; bullet, td[wt:wt]; {blacksquare}, td[wt:Y652AF656A]; {blacktriangleup}, td[wt:Y652A]; {diamondsuit}, td[Y652A:F656A]; {triangledown}, F656A; {triangleup}, Y652A. Voltage dependence of the block was expressed by the gating charge (see Supplemental Data): 0.40 e for WT, 0.48 e for td[wt:wt], 0.37 e for td[wt:F656A], 0.03 e for td[wt:Y652A], 0.01 e for td[Y652A:F656A], 0.21 e for td[wt:Y652AF656A], 0.33 e for F656A, and -0.10 e for Y652A (e is the elementary charge). The value for F656A was corrected for the unblocking (see Supplemental Data).

 
Likewise, the Ki values for mutants at different membrane potentials were obtained (see Supplemental Data) and are also shown in Fig. 4B. Ki values were voltage-dependent for WT, td[wt:wt], td[wt:F656A], td[wt:Y652AF656A], and F656A. td[Y652A:F656A], td[wt:Y652A], and Y652A showed no voltage dependence. It should be noted that voltage-dependence was almost abolished even when only two diagonal Tyr652 residues, not four of them, were deleted.

Kinetics of Blocking. Current traces of some of the mutants as well as WT exhibited slow decays at depolarized potentials, whereas others did not. To examine the mechanisms underlying the different blocking kinetics, we focused on F656A. F656A showed voltage-dependent block, but the currents did not show a decay at depolarized potentials. We found that the shapes of the tail currents were changed in the presence of cisapride (Fig. 5). Expansion of the tail current traces showed that the currents reached their peak with a slower rising phase in the presence of cisapride (Fig. 5, right), which was not observed in the absence of cisapride (Fig. 5, left). The slowly rising phases became more prominent as the preceding voltages were more depolarized. The tail currents of F656A were fitted with a double-exponential function in the absence of cisapride (the recovery from the inactivation and deactivation components) and with a triple-exponential function in the presence of cisapride. The novel kinetic component having an intermediate time constant (~100 ms) between those of recovery from the inactivation and deactivation seems to be related to blocking. In fact, the time constant of the intermediate component was reduced as the cisapride concentration was increased (see Supplemental Data). This is in contrast to the time constant for recovery from the inactivation being unaltered. These results suggest that the intermediate component represented unblocking kinetics. The amplitudes of the unblocking component increased as the prepotential was more depolarized (see Supplemental Data). This is because more channels had been blocked at more depolarized potentials, and those blocked channels released more drug upon repolarization. These are typical voltage-dependent unblocking kinetics (Hille, 2001Go). The plot of the apparent rate of block (equal to the reciprocal of the unblocking time constant) as a function of cisapride concentration showed a nearly linear relationship (Fig. 7A, open diamonds). This kinetic feature will be discussed with respect to the blocking mechanism below.


Figure 5
View larger version (23K):
[in this window]
[in a new window]

 
Fig. 5. Unblocking of cisapride from F656A. Left, current traces at different membrane potentials (top) along with those of tail currents on an expanded time scale in the absence of cisapride (bottom). The right panel shows current traces at 1.6 µM cisapride. A slow rising phase appeared in the presence of cisapride. Kinetic analyses were performed on the tail current traces (see detail in Supplemental Data).

 

Figure 7
View larger version (10K):
[in this window]
[in a new window]

 
Fig. 7. Blocking kinetics. A, the apparent rate of blocking as a function of cisapride concentration. Curves represent fits with eq. 9 (see Materials and Methods). The intercepts of the y-axis represent the values of koff. The koff values (per second) for the mutants were as follows: 0.10 ± 0.13 for WT, 0.01 ± 0.15 for td[wt:wt], 0.01 ± 0.08 for td[wt:F656A], 0.03 ± 0.00 for td[wt:Y652AF656A], 0.32 ± 0.16 for td[wt:Y652A], and 3.1 ± 0.4 for F656A. The k3 values (per second) were 0.015 for WT, 0.022 for td[wt:wt], 0.034 for td[wt:F656A], 0.017 for td[wt:Y652AF656A], 0.017 for td[wt:Y652A], and 0.012 for F656A. The symbols for channel species are shown in A. B, Kf values (micromolar) obtained from fitting the data shown in A. The fit parameters roughly reproduced the Ki values obtained from the concentration dependence of the steady-state block. Nearly linear relationships suggest that experiments were performed in the concentration range below the Kf value for each channel species.

 

For other types of channels, the kinetics of tail current were examined, but kinetic components relevant to the unblocking of cisapride were not found. However, the blocking kinetics was evaluated on a longer time scale from a pulse protocol shown in Fig. 6. Depolarizing pulses to +20 mV were prolonged progressively from 50 to 7000 ms, and the tail currents at -50 mV were recorded. The peaks of the tail currents were plotted as a function of the duration of the depolarization (the envelope current; Fig. 6B, inset). In the absence of cisapride, the current envelopes reached saturation as the depolarization pulses were prolonged (inset, filled squares). These represent the time courses of the activation gating. In the presence of cisapride the envelope currents were depressed as the duration of depolarization was prolonged. This depression was seen for WT, td[wt:wt], td[wt: F656A], td[wt:Y652A], and td[wt:Y652AF656A] (these channels are named here as the high-affinity channels) but not for td[Y652A:F656A], F656A, and Y652A. The envelope currents were fitted with a double-exponential function (Fig. 6B, inset). The fast and slow components represent the time course of activation and block. The envelope experiments were performed at different cisapride concentrations. The time constants of activation for all the channel species were similar and did not change significantly with cisapride concentration. On the other hand, the time constants for the blocking changed significantly. For channel species exhibiting depression of the current envelopes, plots of the apparent rates of block (equal to reciprocal of the blocking time constants) as a function of cisapride concentration showed nearly linear relationships (Fig. 7A).


Figure 6
View larger version (34K):
[in this window]
[in a new window]

 
Fig. 6. Envelope analyses of HERG channels without (A) and with cisapride (B). A, representative traces of envelope current for WT and mutants. Inset, the command potentials. The longest depolarization pulse was 7 s. B, representative traces of envelope current for WT and mutants at the following cisapride concentrations (micromolar): 0.25 for WT, 0.25 for td[wt:wt], 0.63 for td[wt:F656A], 0.25 for td[wt:Y652AF656A], 4.0 for F656A, 0.63 for td[wt:Y652A], 1.6 for td[Y652A:F656A], and 10.0 for Y652A. It should be noted that current traces early in depolarization pulses for pulses of different durations overlapped completely. This signifies that the currents were recorded in the steady state of blocking. Inset, the time course of block for WT. The peaks of the tail currents were plotted as a function of the duration of depolarizing pulses. The time courses of the currents were fitted with a double-exponential function (curved lines). Current data with a fitted line were shown for different cisapride concentrations. As the cisapride concentrations were increased, current amplitudes were depressed and the decay of current (blocking) was accelerated.

 


    Discussion
 Top
 Abstract
 Materials and Methods
 Results
 Discussion
 Conclusion
 References
 
In this study, tandem dimers of HERG channels were constructed for the first time to investigate the mechanisms underlying the block of HERG channels. Introducing single or double mutations into tandem dimers, including single-mutation homotetramers, produced channels with variable arrangements of aromatic residues. The tandem dimers exhibited channel activity with the fast deactivation. This might be a signature that both subunits of the concatenated dimers were incorporated into tetrameric channels: a cytoplasmic domain between concatenated subunits may be restricted in the closely packed subunits which affected the deactivation gating (Morais Cabral et al., 1998Go; Tu and Deutsch, 1999Go; Wang et al., 2000Go). Interactions between cisapride and these mutant HERG channels were examined electrophysiologically, and the steady-state binding, its voltage dependence, and the kinetics of block were analyzed.

Binding Energetics and Topology. The set of mutants revealed the importance of the spatial arrangement of distinct residues to the binding. In this study, we found that cisapride binding was not confined to a single subunit but spans multiple subunits in a tetrameric channel: Tyr652 residues of adjacent subunits contribute to the binding. In contrast, at most two Phe656 residues in diagonal subunits, rather than those in adjacent subunits, were enough for binding, although complete loss of Phe656 (F656A) destabilized the binding dramatically. For the high-affinity binding, retaining intact subunits seems to be important.

For the two double mutants (td[wt:Y652AF656A] and td[Y652A:F656A]), the binding energy differed significantly, which indicates the importance of the arrangements of the aromatic residues in a channel. The contribution of residues to the binding and the interactions between residues upon binding can be evaluated quantitatively by examining the coupling energies ({Delta}{Delta}Gint) in the thermodynamic double-mutant cycle analysis (see Materials and Methods; Schreiber and Fersht, 1995Go; Ranganathan et al., 1996Go; Fersht, 1999Go; Yan et al., 2006Go). For each of two double mutants, a cyclic diagram was drawn (Fig. 8A). In these diagrams, a double mutant can be regarded as the consequences of two successive single mutations reached through two different routes. Inspection of a set of cartoons for a cycle gives ideas how the residues interact upon drug binding.


Figure 8
View larger version (20K):
[in this window]
[in a new window]

 
Fig. 8. Energetics and orientation of cisapride binding. A, the thermodynamic double-mutant cycle analysis. Cartoons for mutant channels are presented without the linker to show differences in special arrangements of residues between mutant channels. Orange arrows represent mutations from Tyr652 to Ala and blue ones for Phe656 -> Ala. There are two pathways to get double mutants: first, Tyr652 -> Ala (the upper orange arrow), then Phe656 -> A (the right blue arrow). Another route is Phe656 -> A (the left blue arrow) first and Tyr652 -> A (the lower orange arrow) later. The upper panel shows the cycle to td[wt:Y652AF656A] and the lower to td[Y652A:F656A]. Numerical values indicate {Delta}Gs (see Materials and Methods). The parallelogram indicates the level of the binding free energy for td[wt:wt], from which up and down arrows were drawn for {Delta}Gs. The changes in the binding energy from td[wt:wt] to td[wt:Y652AF656A] were evaluated by the sums of {Delta}Gs through either of the pathways ({Delta}GsingleM1-WT +{Delta}GdoubleM-singleM1 = {Delta}GsingleM2-WT +{Delta}GdoubleM-singleM2). If two residues contribute to the binding independently, {Delta}G of a mutation at one site should not depend on whether another site is intact or mutated ({Delta}GsingleM1-WT ={Delta}GdoubleM-singleM2). This is seen in the mutations from td[wt:wt] to td[wt:Y652A] ({Delta}G = 3.0 kJ/mol) and from td[wt:F656A] to td[wt:Y652AF656A] ({Delta}G = 2.7 kJ/mol) and the {Delta}{Delta}Gint value was 0.3 kJ/mol. For the double-mutant cycle in the bottom panel, the {Delta}{Delta}Gint value was -3.1 kJ/mol. As seen from the double mutants for both cycles, the coupling energy represents interactions between Tyr652 and Phe656 in the same subunit and in the adjacent subunits. Significant interaction suggests that Tyr652 and Phe656 from different subunits located in proximity. B, chemical structure of cisapride and a hypothetical binding topology. Red lines indicate positive interaction, whereas the blue dotted line indicates negative contributions. In this scheme, bindings are restricted to the adjacent two subunits.

 
For the upper cycle (toward td[wt:Y652AF656A]), the cartoons for four channel species exhibit a common feature; i.e., a diagonal pair of subunits is intact. This means that the double mutant was generated by introducing two mutations successively into subunits adjacent to the intact subunits. From a thermodynamic relationship (eq. 8), {Delta}{Delta}Gint of 0.3 kJ/mol was obtained. This value of {Delta}{Delta}Gint did not deviate statistically from zero, indicating that no interaction took place between Tyr652 and Phe656 upon binding. The results indicate that, as far as intact subunits were retained in a channel, Tyr652 and Phe656 residues on the subunits adjacent to the intact subunits contributed to the binding independently of each other. This additivity is important because it suggests that deleting bulky aromatic rings in two sites did not significantly change the structures of the binding site (Fersht, 1999Go). This is an implicit assumption for studies of mutational analysis in general but has not been supported for HERG channels until this study.

Next, look at another cycle with td[Y652A: F656A] as a double mutant (Fig. 8A, lower panel). The {Delta}{Delta}Gint value of -3.1 kJ/mol indicates significant interactions. The negative sign of {Delta}{Delta}Gint represents that Phe656 and Tyr652 contribute to the binding with negative interactions. What is the origin of this cooperativity? In the double-mutant cycle analysis, significant coupling is explained as being attributable to the proximity of mutated residues (Schreiber and Fersht, 1995Go). The set of cartoons demonstrates that two mutational sites are proximally positioned on the two adjacent subunits. Then the strong interaction can be interpreted that Phe656 and Tyr652 residues on two adjacent subunits are closely apposed.

In the thermodynamic double-mutant cycle analysis, two different types of interactions were elucidated. Tyr652 and Phe656 on the same subunits interact additively (upper cycle) and Tyr652 and Phe656 on adjacent subunits interact cooperatively (lower cycle). One may imagine that the Tyr652 and Phe656 residues from adjacent subunits are placed closer than those placed on the same subunits.

Hypothetical Binding Pose of Cisapride in the High-Affinity Site of HERG Channels. Cisapride contains two aromatic moieties and one basic tertiary nitrogen (Fig. 8B) that can interact with the aromatic amino acids by cation-{pi} and {pi}-{pi} interactions or by hydrophobic interactions (Mitcheson et al., 2000Go; Fernandez et al., 2004Go). Here we propose a possible binding mode of cisapride. Two Tyr652 residues from adjacent subunits and one Phe656 form a favorable binding site, whereas Phe656 in the adjacent subunit hinders the binding (Fig. 8B).

This mode of interaction was compared with the computational docking model built on a homology model of HERG channel (Aronov and Goldman, 2004Go; Osterberg and Aqvist, 2005Go; Farid et al., 2006Go). It was predicted that multiple simultaneous aromatic ring stacking and/or hydrophobic interactions between Tyr652 and Phe656 side chains and aromatic/hydrophobic blocker groups are involved in the docking poses (Farid et al., 2006Go). It was also predicted that Tyr652 and Phe656 from multiple subunits might be involved in cisapride binding (Farid et al., 2006Go). This docking pose is in general agreement with our experimental data.

Voltage-Dependent Block Suggested a Sequential Binding Mechanism. The voltage-dependence of cisapride block has been reported (Walker et al., 1999Go), and in this study, it was shown that the Ki values for WT and some mutants were voltage-dependent (Fig. 4). Not only Y652A but also td[Y652A:F656A] and td[wt:Y652A] exhibited voltage independence. The presence or absence of the voltage dependence for different mutants indicates that voltage-dependent block occurs through voltage-dependent conformational changes of the binding site rather than the binding of a charged drug to the site under the influence of the membrane electric field. To introduce voltage-dependence into the blocking mechanism, a sequential binding model was applied:Formula

in which O represents the open state, B1 represents the first blocked state, and B2 is the subsequent blocked state. [B] indicates the concentration of blocker. In this model, drug-bound channels undergo conformational changes between the B1 and B2 states in a voltage-dependent manner, leading to more stable binding of the drug in the B2 state. A similar model has been proposed for blocking by chloroquine (Sánchez-Chapula et al., 2002Go).

The behavior of channels exhibiting voltage-independent blocking (Y652A, td[wt:Y652A] and td[Y652A:F656A]) can be explained by the absence of transitions from B1 to B2. In these mutants, some of Tyr652 residues were deleted, suggesting that they contribute to the voltage-dependence through changes in their orientation upon depolarization (Sánchez-Chapula et al., 2003Go). In this study, we found that the voltage dependence was almost abolished even for td[wt: Y652A]. On the other hand, td[wt:Y652AF656A] retained voltage dependence. It is suggested that rearrangements of Tyr652 residues seem to be affected by surrounding Phe656 residues.

Kinetic Data Suggest a Low-Affinity Binding Site. Blocking kinetics were evaluated by two different methods: unblocking kinetics after repolarization and slow blocking kinetics by the envelope analyses. In either case, the apparent rate constants of block were found to be a nearly linear function of the cisapride concentrations (Fig. 7A). If the blocking kinetics is described by a simple two-state model, the rate constant of block should be linear. Sublinearity of the blocking rate constant can be accounted for by the presence of an intermediate state of block. Therefore at least three states are necessary to describe the kinetic data, which is compatible with the linear 3-state model proposed for the voltage-dependent block. Sublinearity of the kinetic data demands that the model has a fast low-affinity binding process followed by a slow transition.Formula

where k1 through k4 represent the rate constants, and the white arrow indicates fast transitions.


Figure 9
View larger version (23K):
[in this window]
[in a new window]

 
Fig. 9. A hypothetical binding model for cisapride. A scheme for the longitudinal section of the HERG pore from the cytoplasmic entryway to the central cavity is shown. B1 and B2 represent the low- and high-affinity binding states. Orange arrows indicate voltage-dependent rearrangements of Tyr652 residues.

 
The eq. 9 (see Materials and Methods and Supplemental data) for the apparent blocking rate was used to fit the data and values of Kf (the dissociation constant for the initial binding site), k3 and koff (= k4) were obtained (Fig. 7B, also see Supplemental Data). For the high-affinity channels (WT, td[wt:wt], td[wt:F656A], td[wt:Y652A], and td[wt:Y652AF656A]), the koff values were very slow (< 0.32/s; Fig. 7A legend). This slow unblocking accounts for the apparent absence of the unblocking at the time of peak for tail current measurements: the high affinity blocked state (B2) did not return to the open state through B1 during the early course of the tail currents. Ki values evaluated from the tail currents thus reflect blocking at depolarized potentials. k3 was so slow that the high-affinity binding site scarcely appeared during short depolarizing pulses. Slow appearance of the high-affinity site and slow release from it led to slow accumulation of blocking during repeated depolarizations.

The relationships between the inactivated state and the high-affinity binding state have been discussed. The inactivation state is a state for a nonconducting selectivity filter, whereas the high-affinity binding state represents arrangements of aromatic residues in the central cavity. From a spatial perspective, they are closely located. However, the kinetic analyses suggest that the high-affinity site appears later (Fig. 7A, legends) than the fast inactivation. Then it is likely that the inactivation may induce the transition to the high-affinity state. This is in agreement with the previous papers that inactivation controlled the high-affinity binding (Numaguchi et al., 2000Go).

For F656A, the Kf and koff (Fig. 7, B and C) values were significantly greater, which led to the low affinity of this mutant. Voltage dependence of this channel indicates that the transition between B1 and B2 was retained. Then, the absence of current decay in the envelope experiment is related to the fast transition from the B1 to the B2 state. This feature suggests that Phe656 residues may retard the rearrangements of Tyr652 residues.

An overview of the kinetic features of channel species revealed that only in F656A did the initial binding site show a reduced affinity and the other channels did not. All of the latter possess at least two diagonal Phe656 residues in common. Therefore, it is suggested that the diagonal Phe656 residues may form the initial binding site for cisapride.


    Conclusion
 Top
 Abstract
 Materials and Methods
 Results
 Discussion
 Conclusion
 References
 
Combining the steady-state block data, its voltage dependence and the kinetic data suggest the hypothetical binding process shown in Fig. 9. Cisapride binds to the low-affinity site constituted by diagonal Phe656 residues. Meanwhile, voltage-dependent conformational changes reorient the Tyr652 residues with Phe656 and form the high-affinity site produced by two adjacent Tyr652 residues and Phe656 but hindered by adjacent Phe656.

In this study, a topological picture of the interaction between HERG and cisapride was obtained. Systematic approaches using a set of mutants and thermodynamic and kinetic analyses mapped the binding site and the blocking process. This approach will be applied to understand better the mechanism of action between HERG channels and the broad spectrum of blocking substances.


    Acknowledgements
 
We thank R. Horn (Jefferson Medical College, Philadelphia, PA), Y. Imai (Takeda Pharmaceutical Company Limited, Osaka, Japan), S. Morishima, M. Iwamoto, and C. Edwards for discussion and T. Goto for secretarial assistance.


    Footnotes
 
This work is supported by Grant-in-Aid for Scientific Research from the Ministry of Education, Science, Sports, and Culture, Japan.

ABBREVIATIONS: HERG, human ether-à-go-go related gene product; V1/2, the half activation voltage; Ki, the blocking inhibition constant.

Formula The online version of this article (available at http://molpharm.aspetjournals.org) contains supplemental material. Back

Address correspondence to: Shigetoshi Oiki, Department of Molecular Physiology and Biophysics, Faculty of Medical Sciences, University of Fukui, 23-3 Matsuokashimoaizuki, Eiheiji-cho, Yoshida-gun, Fukui 910-1193, Japan. E-mail: oiki-fki{at}umin.ac.jp


    References
 Top
 Abstract
 Materials and Methods
 Results
 Discussion
 Conclusion
 References
 
Ando H, Kuno M, Shimizu H, Muramatsu I, and Oiki S (2005) Coupled K+-water flux through the HERG potassium channel measured by an osmotic pulse method. J Gen Physiol 126: 529-538.[Abstract/Free Full Text]

Aronov AM and Goldman BB (2004) A model for identifying HERG K+ channel blockers. Bioorg Med Chem 12: 2307-2315.[CrossRef][Medline]

Ashcroft FM (2000) Ion Channels and Disease: Channelopathies. Academic Press, London.

Christie MJ, North RA, Osborne PB, Douglass J, and Adelman JP (1990) Heteropolymeric potassium channels expressed in Xenopus oocytes from cloned subunits. Neuron 4: 405-411.[CrossRef][Medline]

Farid R, Day T, Friesner RA, and Pearlstein RA (2006) New insights about HERG blockade obtained from protein modeling, potential energy mapping, and docking studies. Bioorg Med Chem 14: 3160-3173.[CrossRef][Medline]

Fernandez D, Ghanta A, Kauffman GW, and Sanguinetti MC (2004) Physicochemical features of the hERG channel drug binding site. J Biol Chem 279: 10120-10127.[Abstract/Free Full Text]

Fersht AR (1999) Structure and Mechanism in Protein Science: A Guide to Enzyme Catalysis and Protein Folding. Freeman and Company, New York.

Ficker E, Kuryshev YA, Dennis AT, Obejero-Paz C, Wang L, Hawryluk P, Wible BA, and Brown AM (2004) Mechanisms of arsenic-induced prolongation of cardiac repolarization. Mol Pharmacol 66: 33-44.[Abstract/Free Full Text]

Haverkamp W, Breithardt G, Camm AJ, Janse MJ, Rosen MR, Antzelevitch C, Escande D, Franz M, Malik M, Moss A, et al. (2000) The potential for QT prolongation and pro-arrhythmia by non-anti-arrhythmic drugs: clinical and regulatory implications. Report on a Policy Conference of the European Society of Cardiology. Cardiovasc Res 47: 219-233.[Free Full Text]

Hille B (2001) Ion Channels of Excitable Membranes, 3rd Ed. Sinauer Associates Inc, Sunderland, MA.

Isacoff EY, Jan YN, and Jan LY (1990) Evidence for the formation of heteromultimeric potassium channels in Xenopus oocytes. Nature 345: 530-534.[CrossRef][Medline]

Liman E, Tytgat J, and Hess P (1992) Subunit stoichiometry of a mammalian K+ channel determined by construction of multimeric cDNAs. Neuron 9: 861-871.[CrossRef][Medline]

Milnes JT, Crociani O, Arcangeli A, Hancox JC, and Witchel HJ (2003) Blockade of HERG potassium currents by fluvoxamine: incomplete attenuation by S6 mutations at F656 or Y652. Br J Pharmacol 139: 887-898.[CrossRef][Medline]

Mitcheson JS, Chen J, Lin M, Culberson C, and Sanguinetti MC (2000) A structural basis for drug-induced long QT syndrome. Proc Natl Acad Sci U S A 97: 12329-12333.[Abstract/Free Full Text]

Mitcheson J, Perry M, Stansfeld P, Sanguinetti MC, Witchel H, and Hancox J (2005) Structural determinants for high-affinity block of hERG potassium channels. Novartis Found Symp 266: 136-150.[Medline]

Mohammad S, Zhou Z, Gong Q, and January CT (1997) Blockage of the HERG human cardiac K+ channel by the gastrointestinal prokinetic agent cisapride. Am J Physiol 273: H2534-H2538.[Medline]

Morais Cabral JH, Lee A, Cohen SL, Chait BT, Li M, and Mackinnon R (1998) Crystal structure and functional analysis of the HERG potassium channel N terminus: a eukaryotic PAS domain. Cell 95: 649-655.[CrossRef][Medline]

Numaguchi H, Mullins FM, Johnson JP Jr., Johns DC, Po SS, Yang IC-H, Tomaselli GF, and Balser JR (2000). Circ Res 87: 1012-1018.[Abstract/Free Full Text]

Osterberg F and Aqvist J. (2005) Exploring blocker binding to a homology model of the open hERG K+ channel using docking and molecular dynamics methods. FEBS Lett 579: 2939-2944.[CrossRef][Medline]

Ranganathan R, Lewis JH, and MacKinnon R (1996) Spatial localization of the K+ channel selectivity filter by mutant cycle-based structure analysis. Neuron 16: 131-139.[CrossRef][Medline]

Ruppersberg JP, Schröter KH, Sakmann B, Stocker M, Sewing S, and Pongs O (1990) Heteromultimeric channels formed by rat brain potassium-channel proteins. Nature 345: 535-537.[CrossRef][Medline]

Sánchez-Chapula JA, Navarro-Polanco RA, Culberson C, Chen J, and Sanguinetti MC (2002) Molecular determinants of voltage-dependent human ether-a-go-go related gene (HERG) K+ channel block. J Biol Chem 277: 23587-23595.[Abstract/Free Full Text]

Sánchez-Chapula JA, Ferrer T, Navarro-Polanco RA, and Sanguinetti MC (2003) Voltage-dependent profile of human ether-a-go-go-related gene channel block is influenced by a single residue in the S6 transmembrane domain. Mol Pharmacol 63: 1051-1058.[Abstract/Free Full Text]

Sanguinetti MC and Tristani-Firouzi M (2006) hERG potassium channels and cardiac arrhythmia. Nature 440: 463-469.[CrossRef][Medline]

Schreiber G and Fersht AR (1995) Energetics of protein-protein interactions: analysis of the barnase-barstar interface by single mutaitons and double mutant cycles. J Mol Biol 248: 478-486.[Medline]

Shimizu H, Toyoshima C, and Oiki S (2003) Interaction between tetraethylammonium and permeant cations at the inactivation gate of the HERG potassium channel. Jpn J Physiol 53: 25-34.[CrossRef][Medline]

Tu L and Deutsch C (1999) Evidence for dimerization of dimers in K+ channel assembly Biophys J 76: 2004-2017.[Medline]

Walker BD, Singleton CB, Bursill JA, Wyse KR, Valenzuela SM, Qiu MR, Breit SN, and Campbell TJ (1999) Inhibition of the human ether-a-go-go-related gene (HERG) potassium channel by cisapride: affinity for open and inactivated states. Br J Pharmacol 128: 444-450.[CrossRef][Medline]

Wang J, Myers CD, and Robertson GA (2000) Dynamic control of deactivation gating by a soluble amino-terminal domain in HERG K+ channels. J Gen Physiol 115: 749-758.[Abstract/Free Full Text]

Yan D, Meyer JK, and White MM (2006) Mapping residues in the ligand-binding domain of the 5-HT3 receptor onto d-tubocurarine structure. Mol Pharmacol 70: 571-578.[Abstract/Free Full Text]

Zhou Z, Gong Q, Ye B, Fan Z, Makielski JC, Robertson GA, and January CT (1998) Properties of HERG channels stably expressed in HEK 293 cells studied at physiological temperature. Biophys J 74: 230-241.[Medline]


Related articles in MolPharm:

Refining Insights into High-Affinity Drug Binding to the Human Ether-à-go-go-Related Gene Potassium Channel
Jules C. Hancox and Andrew F. James
MolPharm 2008 73: 1592-1595. [Abstract] [Full Text]  



This article has been cited by other articles:


Home page
Mol. Pharmacol.Home page
M. J. Perrin, P. W. Kuchel, T. J. Campbell, and J. I. Vandenberg
Drug Binding to the Inactivated State Is Necessary but Not Sufficient for High-Affinity Binding to Human Ether-a-go-go-Related Gene Channels
Mol. Pharmacol., November 1, 2008; 74(5): 1443 - 1452.
[Abstract] [Full Text] [PDF]


Home page
Mol. Pharmacol.Home page
J. C. Hancox and A. F. James
Refining Insights into High-Affinity Drug Binding to the Human Ether-a-go-go-Related Gene Potassium Channel
Mol. Pharmacol., June 1, 2008; 73(6): 1592 - 1595.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Data Supplement
Right arrow All Versions of this Article:
mol.107.042085v1
73/6/1643    most recent
Right arrow Submit a response
Right arrow Alert me when this article is cited
Right arrow Alert me when eLetters are posted
Right arrow Alert me if a correction is posted
Services
Right arrow Related articles in MolPharm
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Myokai, T.
Right arrow Articles by Oiki, S.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Myokai, T.
Right arrow Articles by Oiki, S.


Home Help [Feedback] [For Subscribers] [Archive] [Search] [Contents]
All ASPET Journals Molecular Pharmacology Pharmacological Reviews
 Molecular Interventions Drug Metabolism and Disposition

Copyright © 2008 by the American Society for Pharmacology and Experimental Therapeutics