New insights about HERG blockade obtained from protein modeling, potential energy mapping, and docking studies

We created a homology model of the homo-tetrameric pore domain of HERG using the crystal structure of the bacterial potassium channel, KvAP, as a template. We docked a set of known blockers with well-characterized effects on channel function into the lumen of the pore between the selectivity filter and extracellular entrance using a novel docking and refinement procedure incorporating Glide and Prime. Key aromatic groups of the blockers are predicted to form multiple simultaneous ring stacking and hydrophobic interactions among the eight aromatic residues lining the pore. Furthermore, each blocker can achieve these interactions via multiple docking configurations. To further interpret the docking results, we mapped hydrophobic and hydrophilic potentials within the lumen of each refined docked complex. Hydrophilic iso-potential contours define a 'propeller-shaped' volume at the selectivity filter entrance. Hydrophobic contours define a hollow 'crown-shaped' volume located above the 'propeller', whose hydrophobic 'rim' extends along the pore axis between Tyr652 and Phe656. Blockers adopt conformations/binding orientations that closely mimic the shapes and properties of these contours. Blocker basic groups are localized in the hydrophilic 'propeller', forming electrostatic interactions with Ser624 rather than a generally accepted pi-cation interaction with Tyr652. Terfenadine, cisapride, sertindole, ibutilide, and clofilium adopt similar docked poses, in which their N-substituents bridge radially across the hollow interior of the 'crown' (analogous to the hub and spokes of a wheel), and project aromatic/hydrophobic portions into the hydrophobic 'rim'. MK-499 docks with its longitudinal axis parallel to the axis of the pore and 'crown', and its hydrophobic groups buried within the hydrophobic 'rim'.     

High affinity HERG K(+) channel blockade by the antiarrhythmic agent dronedarone: resistance to mutations of the S6 residues Y652 and F656

Pharmacological inhibition of human-ether-a-go-go-related gene (HERG) K(+) channels by structurally and therapeutically diverse drugs is associated with the 'acquired' form of long QT syndrome and with potentially lethal cardiac arrhythmias. Two aromatic amino-acid residues (Y652 and F656) on the inner (S6) helices are considered to be key constituents of a high affinity drug binding site within the HERG channel pore cavity. Using wild-type (WT) and mutant HERG channels expressed in mammalian cell lines, we have investigated HERG channel current (I(HERG)) blockade at 37+/-1 degrees C by dronedarone (DRONED), a non-iodinated analogue of the Class III antiarrhythmic agent amiodarone (AMIOD). Under our conditions WT I(HERG) tails, measured at -40 mV following activating pulses to +30 mV, were blocked with IC(50) values of approximately 59 and 70 nM for DRONED and AMIOD, respectively. I(HERG) inhibition by DRONED was contingent upon channel gating, with block developing rapidly on membrane depolarization, but with no preference for activated over inactivated channels. High external [K(+)] (94 mM) reduced the potency of I(HERG) inhibition by both DRONED and AMIOD. Strikingly, mutagenesis to alanine of the S6 residue F656 (F656A) failed to eliminate blockade by both DRONED and AMIOD, whilst Y652A had comparatively little effect on DRONED but some effect on AMIOD. These findings demonstrate that high affinity drug blockade of I(HERG) can occur without a strong dependence on the Y652 and F656 aromatic amino-acid residues.     

Molecular determinants of voltage-dependent human ether-a-go-go related gene (HERG) K+ channel block

The structural determinants for the voltage-dependent block of ion channels are poorly understood. Here we investigate the voltage-dependent block of wild-type and mutant human ether-a-go-go related gene (HERG) K(+) channels by the antimalarial compound chloroquine. The block of wild-type HERG channels expressed in Xenopus oocytes was enhanced as the membrane potential was progressively depolarized. The IC(50) was 8.4 +/- 0.9 microm when assessed during 4-s voltage clamp pulses to 0 mV. Chloroquine also slowed the apparent rate of HERG deactivation, reflecting the inability of drug-bound channels to close. Mutation to alanine of aromatic residues (Tyr-652 or Phe-656) located in the S6 domain of HERG greatly reduced the potency of channel block by chloroquine (IC(50) > 1 mm at 0 mV). However, mutation of Tyr-652 also altered the voltage dependence of the block. In contrast to wild-type HERG, block of Y652A HERG channels was diminished by progressive membrane depolarization, and complete relief from block was observed at +40 mV. HERG channel block was voltage-independent when the hydroxyl group of Tyr-652 was removed by mutating the residue to Phe. Together these findings indicate a critical role for Tyr-652 in voltage-dependent block of HERG channels. Molecular modeling was used to define energy-minimized dockings of chloroquine to the central cavity of HERG. Our experimental findings and modeling suggest that chloroquine preferentially blocks open HERG channels by cation-pi and pi-stacking interactions with Tyr-652 and Phe-656 of multiple subunits.     

Inhibition of the HERG K+ channel by the antifungal drug ketoconazole depends on channel gating and involves the S6 residue F656

The mechanism of human ether-à-go-go-related gene (HERG) K+ channel blockade by the antifungal agent ketoconazole was investigated using patch-clamp recording from mammalian cell lines. Ketoconazole inhibited whole-cell HERG current (IHERG) with a clinically relevant half-maximal inhibitory drug concentration (IC50) value of 1.7 microM. The voltage- and time-dependent characteristics of IHERG blockade by ketoconazole indicated dependence of block on channel gating, ruling out a significant role for closed-state channel inhibition. The S6 HERG mutations Y652A and F656A produced approximately 4-fold and approximately 21-fold increases in IC50 for IHERG blockade, respectively. Thus, ketoconazole accesses the HERG channel pore-cavity on channel gating, and the S6 residue F656 is an important determinant of ketoconazole binding.     

Characterization of HERG potassium channel inhibition using CoMSiA 3D QSAR and homology modeling approaches

A data set consisting of twenty-two sertindole analogues and ten structurally diverse inhibitors, spanning a wide range in potency, was analyzed using CoMSiA. A homology model of HERG was constructed from the crystal structure of the open MthK potassium channel. A complementary relationship between our CoMSiA and homology models is apparent when the long inhibitor axis is oriented parallel to the longitudinal axis of the pore, with the tail region pointed toward the selectivity filter. The key elements of the pharmacophore, the CoMSiA and the homology model are: (1) The hydrophobic feature optimally consists of an aromatic group that is capable of engaging in pi-stacking with a Phe656 side chain. Optionally, a second aromatic or hydrophobic group present in some inhibitors may contact an additional Phe656 side chain. (2) The basic nitrogen appears to undergo a pi-cation interaction with Tyr652. (3) The pore diameter (12A+), and depth of the selectivity loop relative to the intracellular opening, act as constraints on the conformation-dependent inhibitor dimensions.     

Physicochemical features of the HERG channel drug binding site

Blockade of hERG K(+) channels in the heart is an unintentional side effect of many drugs and can induce cardiac arrhythmia and sudden death. It has become common practice in the past few years to screen compounds for hERG channel activity early during the drug discovery process. Understanding the molecular basis of drug binding to hERG is crucial for the rational design of medications devoid of this activity. We previously identified 2 aromatic residues, Tyr-652 and Phe-656, located in the S6 domain of hERG, as critical sites of interaction with structurally diverse drugs. Here, Tyr-652 and Phe-656 were systematically mutated to different residues to determine how the physicochemical properties of the amino acid side group affected channel block by cisapride, terfenadine, and MK-499. The potency for block by all three drugs was well correlated with measures of hydrophobicity, especially the two-dimensional approximation of the van der Waals hydrophobic surface area of the side chain of residue 656. For residue 652, an aromatic side group was essential for high affinity block, suggesting the importance of a cation-pi interaction between Tyr-652 and the basic tertiary nitrogen of these drugs. hERG also lacks a Pro-Val-Pro motif common to the S6 domain of most other voltage-gated K(+) channels. Introduction of Pro-Val-Pro into hERG reduced sensitivity to drugs but also altered channel gating. Together, these findings assign specific residues to receptor fields predicted by pharmacophore models of hERG channel blockers and provide a refined molecular understanding of the drug binding site.     

Erythromycin block of the HERG K+ channel: accessibility to F656 and Y652

The HERG potassium channel might have a non-canonical drug binding site, distinct from the channel's inner cavity, that could be responsible for elements of closed-state pharmacological inhibition of the channel. The macrolide antibiotic erythromycin is a drug that may block unconventionally because of its size. Here we used whole-cell patch-clamp recording at 37 degrees C from heterologously expressed HERG channels in a mammalian cell line to show that erythromycin either produces a rapid open-state-dependent HERG channel inhibition, or components of both open-state-dependent and closed-state-dependent inhibition. Alanine-substitution of HERG's canonical determinants of blockade revealed that Y652 was not important as a molecular determinant of blockade, and that mutation of F656 resulted in only weak attenuation of inhibition. In computer models of the channel, erythromycin could make several direct contacts with F656, but not with Y652, in the open-state model, and erythromycin was unable to fit into a closed-state channel model.     

Mutational analysis of block and facilitation of HERG current by a class III anti-arrhythmic agent, nifekalant

Chemicals and toxins are useful tools to elucidate the structure-function relationship of various proteins including ion channels. The HERG channel is blocked by many compounds and this may cause life-threatening cardiac arrhythmia. Besides block, some chemicals such as the class III anti-arrhythmic agent nifekalant stimulate HERG at low potentials by shifting its activation curve towards hyperpolarizing voltages. This is called "facilitation". Here, we report mutations and simulations analyzing the association between nifekalant and channel pore residues for block and facilitation. Alanine-scanning mutagenesis was performed in the pore region of HERG. The mutations at the base of the pore helix (T623A), the selectivity filter (V625A) and the S6 helix (G648A, Y652A and F656A) abolished and S624A attenuated both block and facilitation induced by the drug. On the other hand, the mutation of other residues caused either an increase or a decrease in nifekalant-induced facilitation without affecting block. An open-state homology model of the HERG pore suggested that T623, S624, Y652 and F656 faced the central cavity, and were positioned within geometrical range for the drug to be able to interact with all of them at the same time. Of these, S649 was the only polar residue located within possible interaction distance from the drug held in its blocking position. Further mutations and flexible-docking simulations suggest that the size, but not the polarity, of the side chain at S649 is critical for drug induced facilitation.     

Mutational Analysis of Block and Facilitation of HERG Current by A Class III Anti-Arrhythmic Agent,Nifekalant

Chemicals and toxins are useful tools to elucidate the structure-function relationship of various proteins including ion channels. The HERG channel is blocked by many compounds and this may cause life-threatening cardiac arrhythmia. Besides block, some chemicals such as the class III anti-arrhythmic agent nifekalant stimulate HERG at low potentials by shifting its activation curve towards hyperpolarizing voltages. This is called "facilitation". Here, we report mutations and simulations analyzing the association between nifekalant and channel pore residues for block and facilitation. Alanine-scanning mutagenesis was performed in the pore region of HERG. The mutations at the base of the pore helix (T623A), the selectivity filter (V625A) and the S6 helix (G648A, Y652A and F656A) abolished and S624A attenuated both block and facilitation induced by the drug. On the other hand, the mutation of other residues caused either an increase or a decrease in nifekalant-induced facilitation without affecting block. An open-state homology model of the HERG pore suggested that T623, S624, Y652 and F656 faced the central cavity, and were positioned within geometrical range for the drug to be able to interact with all of them at the same time. Of these, S649 was the only polar residue located within possible interaction distance from the drug held in its blocking position. Further mutations and flexible-docking simulations suggest that the size, but not the polarity, of the side chain at S649 is critical for drug induced facilitation.     

In silico analysis of conformational changes induced by mutation of aromatic binding residues: consequences for drug binding in the hERG K+ channel

Pharmacological inhibition of cardiac hERG K(+) channels is associated with increased risk of lethal arrhythmias. Many drugs reduce hERG current by directly binding to the channel, thereby blocking ion conduction. Mutation of two aromatic residues (F656 and Y652) substantially decreases the potency of numerous structurally diverse compounds. Nevertheless, some drugs are only weakly affected by mutation Y652A. In this study we utilize molecular dynamics simulations and docking studies to analyze the different effects of mutation Y652A on a selected number of hERG blockers. MD simulations reveal conformational changes in the binding site induced by mutation Y652A. Loss of π-π-stacking between the two aromatic residues induces a conformational change of the F656 side chain from a cavity facing to cavity lining orientation. Docking studies and MD simulations qualitatively reproduce the diverse experimentally observed modulatory effects of mutation Y652A and provide a new structural interpretation for the sensitivity differences.     

The synergic modeling for the binding of fluoroquinolone antibiotics to the hERG potassium channel

The fluoroquinolone antibiotic binding site in the hERG potassium channel was examined for the residues involved and their position in the tetrameric channel. The blocking effect of the two fluoroquinolones levofloxacin and sparfloxacin to tandem dimers of the hERG mutants were evaluated electrophysiologically. The results indicated that two Tyr652s in the neighboring subunits and one or two Phe656s in the diagonal subunits contributed to the blockade in the case of both compounds, and Ser624 was also involved. The docking studies suggested that the protonated carboxyl group in the compounds strongly interacts with Phe656 as a π acceptor.     

Exploring blocker binding to a homology model of the open hERG K+ channel using docking and molecular dynamics methods

Binding of blockers to the human voltage-gated hERG potassium channel is studied using a combination of homology modelling, automated docking calculations and molecular dynamics simulations, where binding affinities are evaluated using the linear interaction energy method. A homology model was constructed based on the available crystal structure of the bacterial KvAP channel and the affinities of a series of sertindole analogues predicted using this model. The calculations reproduce the relative binding affinities of these compounds very well and indicate that both polar interactions near the intracellular opening of the selectivity filter as well as hydrophobic complementarity in the region around F656 are important for blocker binding. These results are consistent with recent alanine scanning mutation experiments on the blocking of the hERG channel by other compounds.     

Mapping the binding site of a human ether-a-go-go-related gene-specific peptide toxin (ErgTx) to the channel's outer vestibule

The goals of this study are to investigate the mechanism and site of action whereby a human ether-a-go-go-related gene (HERG)-specific scorpion peptide toxin, ErgTx, suppresses HERG current. We apply cysteine-scanning mutagenesis to the S5-P and P-S6 linkers of HERG and examine the resulting changes in ErgTx potency. Data are compared with the characteristics of charybdotoxin (ChTx, or its analogs) binding to the Shaker channel. ErgTx binds to the outer vestibule of HERG but may not physically occlude the pore. In contrast to ChTx.Shaker interaction, elevating [K](o) (from 2 to 98 mm) does not affect ErgTx potency, and through-solution electrostatic forces only play a minor role in influencing ErgTx.HERG interaction. Cysteine mutations of three positions in S5-P linker (Trp-585, Gly-590, and Ile-593) and 1 position in P-S6 linker (Pro-632) induce profound changes in ErgTx binding (DeltaDeltaG > 2 kcal/mol). We propose that the long S5-P linker of the HERG channel forms an amphipathic alpha-helix that, together with the P-S6 linker, forms a hydrophobic ErgTx binding site. This study paves the way for future mutant cycle analysis of interacting residues in the ErgTx.HERG complex, which, in conjunction with NMR determination of the ErgTx solution structure, will yield information about the topology of HERG's outer vestibule.     

Stereospecific interactions of proline residues in protein structures and complexes

The constrained backbone torsion angle of a proline (Pro) residue has usually been invoked to explain its three-dimensional context in proteins. Here we show that specific interactions involving the pyrrolidine ring atoms also contribute to its location in a given secondary structure and its binding to another molecule. It is adept at participating in two rather non-conventional interactions, C-H...pi and C-H...O. The geometry of interaction between the pyrrolidine and aromatic rings, vis-à-vis the occurrence of the C-H...pi interactions has been elucidated. Some of the secondary structural elements stabilized by Pro-aromatic interactions are beta-turns, where a Pro can interact with an adjacent aromatic residue, and in antiparallel beta-sheet, where a Pro in an edge strand can interact with an aromatic residue in the adjacent strand at a non-hydrogen-bonded site. The C-H groups at the Calpha and Cdelta positions can form strong C-H...O interactions (as seen from the clustering of points) and such interactions involving a Pro residue at C' position relative to an alpha-helix can cap the hydrogen bond forming potentials of the free carbonyl groups at the helix C terminus. Functionally important Pro residues occurring at the binding site of a protein almost invariably engage aromatic residues (with one of them being held by C-H...pi interaction) from the partner molecule in the complex, and such aromatic residues are highly conserved during evolution.     

On the origin of the stronger binding of PIB over thioflavin T to protofibrils of the Alzheimer amyloid-β peptide: a molecular dynamics study

Pittsburgh compound B (PIB) is a neutral derivative of the fluorescent dye Thioflavin T (ThT), which displays enhanced hydrophobicity and binding affinity to amyloid fibrils. We present molecular dynamics simulations of binding of PIB and ThT to a common cross-β-subunit of the Alzheimer Amyloid-β peptide (Aβ). Our simulations of binding to Aβ(9-40) protofibrils show that PIB, like ThT, selectively binds to the hydrophobic or aromatic surface grooves on the β-sheet surface along the fibril axis. The lack of two methyl groups and charge in PIB not only improves its hydrophobicity but also leads to a deeper insertion of PIB compared to ThT into the surface grooves. This significantly increases the steric, aromatic, and hydrophobic interactions, and hence leads to stronger binding. Simulations on protofibrils consisting of the more-toxic Aβ(17-42) revealed an additional binding mode in which PIB and ThT insert into the channel that forms in the loop region of the protofibril, sandwiched between two sheet layers. Our simulations indicate that the rotation between the two ring parts of the dyes is significantly more restricted when the dyes are bound to the surface of the cross-β-subunits or to the channel inside the Aβ(17-42) cross-β-subunit, compared with free solution. The specific conformations of the dyes are influenced by small chemical modifications (ThT versus PIB) and by the environment in which the dye is placed.     

酮色林对人类ether-a-go-go相关基因钾通道的阻断作用

采用双电极电压钳技术,研究酮色林对表达在非洲爪蟾卵母细胞上的野生型和Y652突变型人类ether-a-go-go相关基因(human ether-a-go-go-related gene,HERG)钾通道的阻断效应,观测HERG通道的分子位点特性改变对其阻断效应的影响.结果显示,酮色林以电压依赖性和浓度依赖性的方式阻断野生型的HERG钾通道电流.尾电流包裹程序记录电流显示酮色林对HERG钾通道微小的张力性阻断.阻断特征符合对开放状态通道的阻断特征.酮色林也能调节失活状态的HERG钾通道.位于孔道S6区的氨基酸位点突变Y652A和Y652R可显著减弱酮色林对HERG通道的阻断作用.同野生犁HERG钾通道的阻断相比,Y652A突变使阻断的IC50提高72倍,而Y652R突变使阻断的IC50提高53倍.Y652A和Y652R的阴断效应之间没有明显的差别.以上结果提示,酮色林优先阻断开放状态的HERG钾通道,而Y652是酮色林与通道结合的关键位点之一.     

Spectroscopic and structural study of proton and halide ion cooperative binding to gfp

This study reports the influence of halogens on fluorescence properties of the Aequorea victoria Green Fluorescent Protein variant S65T/T203Y (E(2)GFP). Halide binding forms a specific nonfluorescent complex generating a substantial drop of the fluorescence via static quenching. Spectroscopic analysis under different solution conditions reveals high halogen affinity, which is strongly dependent on the pH. This evidences the presence in E(2)GFP of interacting binding sites for halide ions and for protons. Thermodynamic link and cooperative interaction are assessed demonstrating that binding of one halide ion is associated with the binding of one proton in a cooperative fashion with the formation, in the pH range 4.5-10, of a single fully protonated E(2)GFP.halogen complex. To resolve the structural determinants of E(2)GFP sensitivity to halogens, high-resolution crystallographic structures were obtained for the halide-free and I(-), Br(-), and Cl(-) bound E(2)GFP. Remarkably the first high-resolution (1.4 A) crystallographic structure of a chloride-bound GFP is reported. The chloride ion occupies a specific and unique binding pocket in direct contact (3.4 A) with the chromophore imidazolidinone aromatic ring. Unanticipated flexibility, strongly modulated by halide ion interactions, is observed in the region surrounding the chromophore. Furthermore molecular dynamics simulations identified E222 residue (along with the chromophore Y66 residue) being in the protonated state when E(2)GFP.halogen complex is formed. The impact of these results on high-sensitivity biosensor design will be discussed.     

Polyamines and the NMDA receptor: modifying intrinsic activities with aromatic substituents

Thirty-four spermidine (SPD) and spermine (SPM) derivatives with aromatic substituents were synthesized and tested as inhibitors of specific binding of the NMDA channel blocker [3H]MK-801 to membranes prepared from rat hippocampus and cerebral cortex. SPD and SPM derivatives with aromatic substituents at the primary amino groups were the most potent inhibitors (IC50 3.9-4.7 microM). These compounds most likely act directly at the NMDA ion channel, since 30 microM SPM had no pronounced influence on their inhibiting activities. SPD derivatives with aromatic substituents at the secondary amino group were either inactive or highly SPM-sensitive inhibitors (IC50 10-82 microM), depending on the size of the substituent. Our results support the hypothesis that an aromatic interaction site near the center of polyamine derivatives contributes to polyamine inverse agonism.     

Structural and functional analysis of natrin, a venom protein that targets various ion channels

Cysteine-rich secretory proteins (CRISPs) are secreted single-chain proteins found in different sources. Natrin is a member of the CRISP family purified from the snake venom of Naja naja atra, which has been reported as a BKca channel blocker. In our study, crystals of natrin were obtained in two different crystal forms and the structure of one of them was solved at a resolution of 1.68A. Our electrophysiological experiments indicated that natrin can block the ion channel currents of the voltage-gated potassium channel Kv1.3. Docking analyses of the interaction between natrin and Kv1.3 revealed a novel interaction pattern different from the two previously reported K(+) channel inhibition models termed "functional dyad" and "basic ring". These findings offered new insights into the function of natrin and how the specific interactions between CRISPs and different ion channels can be achieved.     

A cation-pi interaction between extracellular TEA and an aromatic residue in potassium channels

Open-channel blockers such as tetraethylammonium (TEA) have a long history as probes of the permeation pathway of ion channels. High affinity blockade by extracellular TEA requires the presence of an aromatic amino acid at a position that sits at the external entrance of the permeation pathway (residue 449 in the eukaryotic voltage-gated potassium channel Shaker). We investigated whether a cation-pi interaction between TEA and such an aromatic residue contributes to TEA block using the in vivo nonsense suppression method to incorporate a series of increasingly fluorinated Phe side chains at position 449. Fluorination, which is known to decrease the cation-pi binding ability of an aromatic ring, progressively increased the inhibitory constant K(i) for the TEA block of Shaker. A larger increase in K(i) was observed when the benzene ring of Phe449 was substituted by nonaromatic cyclohexane. These results support a strong cation-pi component to the TEA block. The data provide an empirical basis for choosing between Shaker models that are based on two classes of reported crystal structures for the bacterial channel KcsA, showing residue Tyr82 in orientations either compatible or incompatible with a cation-pi mechanism. We propose that the aromatic residue at this position in Shaker is favorably oriented for a cation-pi interaction with the permeation pathway. This choice is supported by high level ab initio calculations of the predicted effects of Phe modifications on TEA binding energy.