Drug block of the hERG potassium channel: insight from modeling

Many commonly used, structurally diverse, drugs block the human ether-a-go-go-related gene (hERG) K(+) channel to cause acquired long QT syndrome, which can lead to sudden death via lethal cardiac arrhythmias. This undesirable side effect is a major hurdle in the development of safe drugs. To gain insight about the structure of hERG and the nature of drug block we have produced structural models of the channel pore domain, into each of which we have docked a set of 20 hERG blockers. In the absence of an experimentally determined three-dimensional structure of hERG, each of the models was validated against site-directed mutagenesis data. First, hERG models were produced of the open and closed channel states, based on homology with the prokaryotic K(+) channel crystal structures. The modeled complexes were in partial agreement with the mutagenesis data. To improve agreement with mutagenesis data, a KcsA-based model was refined by rotating the four copies of the S6 transmembrane helix half a residue position toward the C-terminus, so as to place all residues known to be involved in drug binding in positions lining the central cavity. This model produces complexes that are consistent with mutagenesis data for smaller, but not larger, ligands. Larger ligands could be accommodated following refinement of this model by enlarging the cavity using the inherent flexibility about the glycine hinge (Gly648) in S6, to produce results consistent with the experimental data for the majority of ligands tested.     

Probing the outer mouth structure of the HERG channel with peptide toxin footprinting and molecular modeling

Previous studies have shown that the unusually long S5-P linker lining human ether a-go-go related gene's (hERG's) outer vestibule is critical for its channel function: point mutations at high-impact positions here can interfere with the inactivation process and, in many cases, also reduce the pore's K+ selectivity. Because no data are available on the equivalent region in the available K channel crystal structures to allow for homology modeling, we used alternative approaches to model its three-dimensional structure. The first part of this article describes mutant cycle analysis used to identify residues on hERG's outer vestibule that interact with specific residues on the interaction surface of BeKm-1, a peptide toxin with known NMR structure and a high binding affinity to hERG. The second part describes molecular modeling of hERG's pore domain. The transmembrane region was modeled after the crystal structure of KvAP pore domain. The S5-P linker was docked to the transmembrane region based on data from previous NMR and mutagenesis experiments, as well as a set of modeling criteria. The models were further restrained by contact points between hERG's outer vestibule and the bound BeKm-1 toxin molecule deduced from the mutant cycle analysis. Based on these analyses, we propose a working model for the open conformation of the outer vestibule of the hERG channel, in which the S5-P linkers interact with the pore loops to influence ion flux through the pore.     

The free energy landscape of small molecule unbinding

The spontaneous dissociation of six small ligands from the active site of FKBP (the FK506 binding protein) is investigated by explicit water molecular dynamics simulations and network analysis. The ligands have between four (dimethylsulphoxide) and eleven (5-diethylamino-2-pentanone) non-hydrogen atoms, and an affinity for FKBP ranging from 20 to 0.2 mM. The conformations of the FKBP/ligand complex saved along multiple trajectories (50 runs at 310 K for each ligand) are grouped according to a set of intermolecular distances into nodes of a network, and the direct transitions between them are the links. The network analysis reveals that the bound state consists of several subbasins, i.e., binding modes characterized by distinct intermolecular hydrogen bonds and hydrophobic contacts. The dissociation kinetics show a simple (i.e., single-exponential) time dependence because the unbinding barrier is much higher than the barriers between subbasins in the bound state. The unbinding transition state is made up of heterogeneous positions and orientations of the ligand in the FKBP active site, which correspond to multiple pathways of dissociation. For the six small ligands of FKBP, the weaker the binding affinity the closer to the bound state (along the intermolecular distance) are the transition state structures, which is a new manifestation of Hammond behavior. Experimental approaches to the study of fragment binding to proteins have limitations in temporal and spatial resolution. Our network analysis of the unbinding simulations of small inhibitors from an enzyme paints a clear picture of the free energy landscape (both thermodynamics and kinetics) of ligand unbinding.     

Saxitoxin is a gating modifier of HERG K+ channels

Potassium (K+) channels mediate numerous electrical events in excitable cells, including cellular membrane potential repolarization. The hERG K+ channel plays an important role in myocardial repolarization, and inhibition of these K+ channels is associated with long QT syndromes that can cause fatal cardiac arrhythmias. In this study, we identify saxitoxin (STX) as a hERG channel modifier and investigate the mechanism using heterologous expression of the recombinant channel in HEK293 cells. In the presence of STX, channels opened slower during strong depolarizations, and they closed much faster upon repolarization, suggesting that toxin-bound channels can still open but are modified, and that STX does not simply block the ion conduction pore. STX decreased hERG K+ currents by stabilizing closed channel states visualized as shifts in the voltage dependence of channel opening to more depolarized membrane potentials. The concentration dependence for steady-state modification as well as the kinetics of onset and recovery indicate that multiple STX molecules bind to the channel. Rapid application of STX revealed an apparent "agonist-like" effect in which K+ currents were transiently increased. The mechanism of this effect was found to be an effect on the channel voltage-inactivation relationship. Because the kinetics of inactivation are rapid relative to activation for this channel, the increase in K+ current appeared quickly and could be subverted by a decrease in K+ currents due to the shift in the voltage-activation relationship at some membrane potentials. The results are consistent with a simple model in which STX binds to the hERG K+ channel at multiple sites and alters the energetics of channel gating by shifting both the voltage-inactivation and voltage-activation processes. The results suggest a novel extracellular mechanism for pharmacological manipulation of this channel through allosteric coupling to channel gating.     

Na+ permeation and block of hERG potassium channels

The inactivation gating of hERG channels is important for the channel function and drug-channel interaction. Whereas hERG channels are highly selective for K+, we have found that inactivated hERG channels allow Na+ to permeate in the absence of K+. This provides a new way to directly monitor and investigate hERG inactivation. By using whole cell patch clamp method with an internal solution containing 135 mM Na+ and an external solution containing 135 mM NMG+, we recorded a robust Na+ current through hERG channels expressed in HEK 293 cells. Kinetic analyses of the hERG Na+ and K+ currents indicate that the channel experiences at least two states during the inactivation process, an initial fast, less stable state followed by a slow, more stable state. The Na+ current reflects Na+ ions permeating through the fast inactivated state but not through the slow inactivated state or open state. Thus the hERG Na+ current displayed a slow inactivation as the channels travel from the less stable, fast inactivated state into the more stable, slow inactivated state. Removal of fast inactivation by the S631A mutation abolished the Na+ current. Moreover, acceleration of fast inactivation by mutations T623A, F627Y, and S641A did not affect the hERG Na+ current, but greatly diminished the hERG K+ current. We also found that external Na+ potently blocked the hERG outward Na+ current with an IC50 of 3.5 mM. Mutations in the channel pore and S6 regions, such as S624A, F627Y, and S641A, abolished the inhibitory effects of external Na+ on the hERG Na+ current. Na+ permeation and blockade of hERG channels provide novel ways to extend our understanding of the hERG gating mechanisms.     

Structural Properties of PAS Domains from the KCNH Potassium Channels

KCNH channels form an important family of voltage gated potassium channels. These channels include a N-terminal Per-Arnt-Sim (PAS) domain with unknown function. In other proteins PAS domains are implicated in cellular responses to environmental queues through small molecule binding or involvement in signaling cascades. To better understand their role we characterized the structural properties of several channel PAS domains. We determined high resolution structures of PAS domains from the mouse EAG (mEAG), drosophila ELK (dELK) and human ERG (hERG) channels and also of the hERG domain without the first nine amino acids. We analyzed these structures for features connected to ligand binding and signaling in other PAS domains. In particular, we have found cavities in the hERG and mEAG structures that share similarities with the ligand binding sites from other PAS domains. These cavities are lined by polar and apolar chemical groups and display potential flexibility in their volume. We have also found that the hydrophobic patch on the domain β-sheet is a conserved feature and appears to drive the formation of protein-protein contacts. In addition, the structures of the dELK domain and of the truncated hERG domain revealed the presence of N-terminal helices. These helices are equivalent to the helix described in the hERG NMR structures and are known to be important for channel function. Overall, these channel domains retain many of the PAS domain characteristics known to be important for cell signaling.     

hERG K+ channel blockade by the antipsychotic drug thioridazine: An obligatory role for the S6 helix residue F656

The phenothiazine antipsychotic agent thioridazine has been linked with prolongation of the QT interval on the electrocardiogram, ventricular arrhythmias, and sudden death. Although thioridazine is known to inhibit cardiac hERG K(+) channels there is little mechanistic information on this action. We have investigated in detail hERG K(+) channel current (I(hERG)) blockade by thioridazine and identified a key molecular determinant of blockade. Whole-cell I(hERG) measurements were made at 37 degrees C from human embryonic kidney (HEK-293) cells expressing wild-type and mutant hERG channels. Thioridazine inhibited I(hERG) tails at -40mV following a 2s depolarization to +20mV with an IC(50) value of 80nM. Comparable levels of I(hERG) inhibition were seen with physiological command waveforms (ventricular and Purkinje fibre action potentials). Thioridazine block of I(hERG) was only weakly voltage-dependent, though the time dependence of I(hERG) inhibition indicated contingency of blockade upon channel gating. The S6 helix point mutation F656A almost completely abolished, and the Y652A mutation partially attenuated, I(hERG) inhibition by thioridazine. In summary, thioridazine is one of the most potent hERG K(+) channel blockers amongst antipsychotics, exhibiting characteristics of a preferential open/activated channel blocker and binding at a high affinity site in the hERG channel pore.     

Mechanism of block of the hERG K+ channel by the scorpion toxin CnErg1

The scorpion toxin CnErg1 binds to human ether-a-go-go related gene (hERG) K(+) channels with a 1:1 stoichiometry and high affinity. However, in contrast to other scorpion toxin-ion channel interactions, the inhibition of macroscopic hERG currents by high concentrations of CnErg1 is incomplete. In this study, we have probed the molecular basis for this incomplete inhibition. High concentrations of CnErg1 had only modest effects on hERG gating that could not account for the incomplete block. Furthermore, the residual current in the presence of 1 microM CnErg1 had normal single channel conductance. Analysis of the kinetics of CnErg1 interaction with hERG indicated that CnErg1 binding is not diffusion-limited. A bimolecular binding scheme that incorporates an initial encounter complex and permits normal ion conduction was able to completely reproduce both the kinetics and steady-state level of CnErg1-hERG binding. This scheme provides a simple kinetic explanation for incomplete block; that is, relatively fast backward compared to forward rate constants for the interconversion of the toxin-channel encounter complex and the blocked toxin-channel complex. We have also examined the temperature-dependence of CnErg1 binding to hERG. The dissociation constant, K(d), for CnErg1 increases from 7.3 nM at 22 degrees C to 64 nM at 37 degrees C (i.e., the affinity decreases as temperature increases) and the proportion of binding events that lead to channel blockade decreases from 70% to 40% over the same temperature range. These temperature-dependent effects on CnErg1 binding correlate with a temperature-dependent decrease in the stability of the putative CnErg1 binding site, the amphipathic alpha-helix in the outer pore domain of hERG, assayed using circular dichroism spectropolarimetry. Collectively, our data provides a plausible kinetic explanation for incomplete blockade of hERG by CnErg1 that is consistent with the proposed highly dynamic conformation of the outer pore domain of hERG.     

The thermodynamic landscape of testosterone binding to cytochrome P450 3A4: ligand binding and spin state equilibria

Human cytochrome P450 (CYP) 3A4 catalyzes the oxygen-dependent metabolism of greater than 60% of known drugs. CYP3A4 binds multiple ligands simultaneously, and this contributes to complex allosteric kinetic behavior. Substrates that bind to this enzyme change the ferric spin state equilibrium of the heme, which can be observed by optical absorbance and electron paramagnetic resonance (EPR) spectroscopy. The ligand-dependent spin state equilibrium has not been quantitatively understood for any ligands that exhibit multiple binding. The CYP3A4 substrate testosterone (TST) has been shown previously by absorbance spectroscopy to induce spin state changes that are characteristic of a low spin to high spin conversion. Here, EPR was used to examine the equilibrium binding of TST to CYP3A4 at [CYP3A4] > K(D), which allows for characterization of the singly occupied state (i.e., CYP3A4.TST). We also have used absorbance spectroscopy to examine equilibrium binding, where [CYP3A4] < K(D), which allows for determination of K(D)'s. The combination of absorbance and EPR spectroscopy at different CYP3A4 concentrations relative to K(D) and curve fitting of the resultant equilibrium binding titration curves to the Adair-Pauling equations, and modifications of it, reveals that the first equivalent of TST binds with higher affinity than the second equivalent of TST and its binding is positively cooperative with respect to ligand-dependent spin state conversion. Careful analysis of the EPR and absorbance spectral results suggests that the binding of the second TST induces a shift to the high spin state and thus that the second TST binding causes displacement of the bound water. A model involving six thermodynamic states is presented and this model is related to the turnover of the enzyme.     

The effect of multiple binding modes on empirical modeling of ligand docking to proteins.

A popular approach to the computational modeling of ligand/receptor interactions is to use an empirical free energy like model with adjustable parameters. Parameters are learned from one set of complexes, then used to predict another set. To improve these empirical methods requires an independent way to study their inherent errors. We introduce a toy model of ligand/receptor binding as a workbench for testing such errors. We study the errors incurred from the two state binding assumption--the assumption that a ligand is either bound in one orientation, or unbound. We find that the two state assumption can cause large errors in free energy predictions, but it does not affect rank order predictions significantly. We show that fitting parameters using data from high affinity ligands can reduce two state errors; so can using more physical models that do not use the two state assumption. We also find that when using two state models to predict free energies, errors are more severe on high affinity ligands than low affinity ligands. And we show that two state errors can be diagnosed by systematically adding new binding modes when predicting free energies: if predictions worsen as the modes are added, then the two state assumption in the fitting step may be at fault.     

Opening mechanism of a cyclic nucleotide-gated channel based on analysis of single channels locked in each liganded state.

Cyclic nucleotide-gated channels contain four subunits, each with a binding site for cGMP or cAMP in the cytoplasmic COOH-terminal domain. Previous studies of the kinetic mechanism of activation have been hampered by the complication that ligands are continuously binding and unbinding at each of these sites. Thus, even at the single channel level, it has been difficult to distinguish changes in behavior that arise from a channel with a fixed number of ligands bound from those that occur upon the binding and unbinding of ligands. For example, it is often assumed that complex behaviors like multiple conductance levels and bursting occur only as a consequence of changes in the number of bound ligands. We have overcome these ambiguities by covalently tethering one ligand at a time to single rod cyclic nucleotide-gated channels (Ruiz, ML., and J.W. Karpen. 1997. Nature. 389:389-392). We find that with a fixed number of ligands locked in place the channel freely moves between three conductance states and undergoes bursting behavior. Furthermore, a thorough kinetic analysis of channels locked in doubly, triply, and fully liganded states reveals more than one kinetically distinguishable state at each conductance level. Thus, even when the channel contains a fixed number of bound ligands, it can assume at least nine distinct states. Such complex behavior is inconsistent with simple concerted or sequential allosteric models. The data at each level of liganding can be successfully described by the same connected state model (with different rate constants), suggesting that the channel undergoes the same set of conformational changes regardless of the number of bound ligands. A general allosteric model, which postulates one conformational change per subunit in both the absence and presence of ligand, comes close to providing enough kinetically distinct states. We propose an extension of this model, in which more than one conformational change per subunit can occur during the process of channel activation.     

Activation gating of hERG potassium channels: S6 glycines are not required as gating hinges

The opening of ion channels is proposed to arise from bending of the pore inner helices that enables them to pivot away from the central axis creating a cytosolic opening for ion diffusion. The flexibility of the inner helices is suggested to occur either at a conserved glycine located adjacent to the selectivity filter (glycine gating hinge) and/or at a second site occupied by glycine or proline containing motifs. Sequence alignment with other K+ channels shows that hERG possesses glycine residues (Gly648 and Gly657) at each of these putative hinge sites. In apparent contrast to the hinge hypotheses, substitution of both glycine residues for alanine causes little effect on either the voltage-dependence or kinetics of channel activation, and open state block by intracellular blockers. Substitution of the glycines with larger hydrophobic residues causes a greater propensity for the channel to open. We propose that in contrast to Shaker the pore of hERG is intrinsically more stable in the open than the closed conformation and that substitution at Gly648 or Gly657 further shifts the gating equilibrium to favor the open state. Molecular dynamics simulations indicate the S6 helices of hERG are inherently flexible, even in the absence of the glycine residues. Thus hERG activation gating exhibits important differences to other Kv channels. Our findings indicate that the hERG inner helix glycine residues are required for the tight packing of the channel helices and that the flexibility afforded by glycine or proline residues is not universally required for activation gating.     

Interaction simulation of hERG K+ channel with its specific BeKm-1 peptide: insights into the selectivity of molecular recognition

Potassium channels show a huge variability in the affinity when recognizing enormous bioactive peptides, and the elucidation of their recognition mechanism remains a great challenge due to an undetermined peptide-channel complex structure. Here, we employed combined computation methods to study the specific binding of BeKm-1 peptide to the hERG potassium channel, which is an essential determinant of the long-QT syndrome. By the use of a segment-assembly homology modeling method, the closed-state hERG structure containing unusual longer S5P linker was successfully constructed. It has a "petunia" shape, while four "petals" of symmetrically distributed S5P segments always decentralize. Starting from the hERG and BeKm-1 structures, a considerably reasonable BeKm-1-hERG complex structure was then screened out and identified by protein-protein docking, molecular dynamics (MD) simulations, and calculation of relative binding free energies. The validity of this predicted complex was further assessed by computational alanine-scanning, with the results correlating reasonably well with experimental data. In the novel complex structure, four considerably flexible S5P linkers are far from the BeKm-1 peptide. The BeKm-1 mainly uses its helical region to associate the channel outer vestibule, except for the S5P linker region; however, structural analysis further implies this neutral pore region with wiggling S5P linker is highly beneficial to the binding of BeKm-1 with lower positive charges. The most critical Lys18 of BeKm-1 plugs its side chain into the channel selectivity filter, while the secondarily important Arg20 forms three hydrogen bonds with spatially neighboring residues in the hERG channel. Different from the classical peptide-K+ channel interaction mainly induced by electrostatic interaction, a synergetic effect of the electrostatic and van der Waals interactions was found to mediate the molecular recognition between BeKm-1 and the hERG channel. And this specific binding process is revealed to be a dynamic change of reduction of binding free energy and conformational rearrangement mainly in the interface of both BeKm-1 and the hERG channel. All these structural and energy features yield deep insights on the high selective binding mechanism of hERG-specific peptides, present a diversity of peptide-K+ channel interactions, and also provide important clues to further study structure-function relationships of the hERG channel.     

Rb+ flux through hERG channels affects the potency of channel blocking drugs: correlation with data obtained using a high-throughput Rb+ efflux assay

The nonradioactive Rb+ efflux assay has become a reliable and efficient high-throughput hERG screening method, but it is limited by its low sensitivity for potent hERG blockers. Using the patch clamp technique, the authors found that the low sensitivity is due in part to the use of Rb+ as the permeating cation in the assay. The affinities of the drugs measured by patch clamp technique in the presence of Rb+ were 3- to 10-fold lower than when measured by the same method in the presence of K+ ions. The apparent affinity of the drugs decreased even further when monitored by the Rb+ efflux assay. It was also observed that Rb+ had minimal effects on the activation properties of channels while there was a significant change in the half-inactivation potential. This voltage shift reduces hERG channel inactivation at efflux assay potentials, and will reduce the affinity of hERG-blocking drugs that bind to inactivated states of the channel. In combination with the effects of elevated extracellular ion concentrations, it is likely that Rb+ modulation of hERG channel inactivation is largely responsible for the reduced drug potencies observed in the Rb+ efflux assay.     

Unique interaction of scorpion toxins with the hERG channel

ERG potassium channels specify one component of the delayed rectifier in the heart and are likely to play an important functional role in other excitable cells. Compared to other K+ channels, the human ERG (hERG) channel possesses an unusually long S5-P linker that presumably forms an alpha-helix important for channel function. hERG-specific toxins bind to the outer mouth of the hERG channel. Channel residues in the middle of the S5-P linker and at the pore entrance are critical for toxin binding. One of these scorpion toxins is BeKm-1. Residues critical for BeKm-1 binding to the hERG channel are located in the alpha-helix and the following loop, whereas the "traditional" interaction surface of other short scorpion toxins is formed by residues on the beta-sheet. This unique localization of BeKm-1's interaction surface and its specific action on the hERG channel suggest a unique outer mouth structure of the hERG channel. We used the mutant cycle analysis approach to define contacts in the toxin-channel complex. This information provides critical constraints and is important for molecular modeling of the hERG pore structure.     

Structural refinement of the hERG1 pore and voltage‐sensing domains with ROSETTA‐membrane and molecular dynamics simulations

The hERG1 gene (Kv11.1) encodes a voltage‐gated potassium channel. Mutations in this gene lead to one form of the Long QT Syndrome (LQTS) in humans. Promiscuous binding of drugs to hERG1 is known to alter the structure/function of the channel leading to an acquired form of the LQTS. Expectably, creation and validation of reliable 3D model of the channel have been a key target in molecular cardiology and pharmacology for the last decade. Although many models were built, they all were limited to pore domain. In this work, a full model of the hERG1 channel is developed which includes all transmembrane segments. We tested a template‐driven de‐novo design with ROSETTA‐membrane modeling using side‐chain placements optimized by subsequent molecular dynamics (MD) simulations. Although backbone templates for the homology modeled parts of the pore and voltage sensors were based on the available structures of KvAP, Kv1.2 and Kv1.2‐Kv2.1 chimera channels, the missing parts are modeled de‐novo. The impact of several alignments on the structure of the S4 helix in the voltage‐sensing domain was also tested. Herein, final models are evaluated for consistency to the reported structural elements discovered mainly on the basis of mutagenesis and electrophysiology. These structural elements include salt bridges and close contacts in the voltage‐sensor domain; and the topology of the extracellular S5‐pore linker compared with that established by toxin foot‐printing and nuclear magnetic resonance studies. Implications of the refined hERG1 model to binding of blockers and channels activators (potent new ligands for channel activations) are discussed. Proteins 2010. © 2010 Wiley‐Liss, Inc.     

In silico design of novel hERG-neutral sildenafil-like PDE5 inhibitors

Cyclic nucleotide phosphodiesterase enzymes (PDEs) have functions in regulating the levels of intracellular second messengers, 3′, 5′-cyclic adenosine monophosphate (cAMP) and 3′, 5′-cyclic guanosine monophosphate (cGMP), via hydrolysis and decomposing mechanisms in cells. They take essential roles in modulating various cellular activities such as memory and smooth muscle functions. PDE type 5 (PDE5) inhibitors enhance the vasodilatory effects of cGMP in the corpus cavernosum and they are used to treat erectile dysfunction. Patch clamp experiments showed that the IC₅₀ values of the human ether-à-go-go-related gene (hERG1) potassium (K) ion channel blocking affinity of PDE5 inhibitors sildenafil, vardenafil, and tadalafil as 33, 12, and 100 μM, respectively. hERG1 channel is responsible for the regulation of the action potential of human ventricular myocyte by contributing the rapid component of delayed rectifier K⁺ current (I_Kr) component of the cardiac action potential. In this work, interaction patterns and binding affinity predictions of selected PDE5 inhibitors against the hERG1 channel are studied. It is attempted to develop PDE5 inhibitor analogs with lower binding affinity to hERG1 ion channel while keeping their pharmacological activity against their principal target PDE5 using in silico methods. Based on detailed analyses of docking poses and predicted interaction energies, novel analogs of PDE5 inhibitors with lower predicted binding affinity to hERG1 channels without loosing their principal target activity were proposed. Moreover, molecular dynamics (MD) simulations and post-processing MD analyses (i.e. Molecular Mechanics/Generalized Born Surface Area calculations) were performed. Detailed analysis of molecular simulations helped us to better understand the PDE5 inhibitor–target binding interactions in the atomic level. Results of this study can be useful for designing of novel and safe PDE5 inhibitors with enhanced activity and other tailored properties.     

Modulation of high affinity phenylalkylamine binding sites on cultured human embryonal vascular smooth muscle cells.

Two phenylalkylamine Ca2+ channel ligands, (+/-)-[3H]verapamil ((+/-)-[3H]V) (-)-[3H]desmethoxyverapamil ((-)-[3H]DV), were employed in whole cell binding assays to characterize the specific high affinity binding sites on Ca2+ channels, their cooperativity and modulations induced on cultured human embryonal vascular smooth muscle preparation (VSM) by: 1) Beta-adrenergic stimulation of the cell, 2) exposure to high K+ concentration, 3) exposure to high concentration of Mg2+ ions, 4) the presence of a benzothiazepine Ca2+ channel antagonist and modulator d-cis-diltiazem, and 5) guanylylimidodiphosphate. The total amounts of specific (+/-)-[3H]V and (-)-[3H]DV binding sites present on VSM cells increased significantly after beta-adrenergic receptor activation, following cell membrane depolarization induced by high concentrations of K+, in the presence of Ca2+ chelator Na3EDTA, and after incubation of VSM cells with a benzothiazine-type Ca2+ channel blocker d-cis-diltiazem. A marked reduction of (-)-[3H]DV binding was observed after permanent G-protein activation by a nonhydrolyzable analog of guanylylimidodiphosphate, after incubation of the cells with norepinephrine, and after incubation of VSM cells with millimolar concentration of Mg2+. The results suggest the existence of multiple modulations of specific (-)-[3H]DV binding sites on Ca2+ channel corresponding to the way of activation of the cell and also to the immediate "state" of the membrane bound Ca2+ channels present on VSM cells, the positive heterotropic interaction after beta-adrenergic stimulation, the homotropic positive allosteric interaction induced by d-cis-diltiazem and pure noncompetitive inhibition induced by guanylylimidodiphosphate. The presence of high concentrations of Mg2+ inhibited whereas the presence of Ca2+ chelator, of ethylenediamine-tetraacetic acid sodium salt, significantly increased the total number of specific high affinity (-)-[3H]DV binding sites on VSM cells.     

Molecular determinants of loperamide and N-desmethyl loperamide binding in the hERG cardiac K+ channel

Abuse of the common anti-diarrheal loperamide is associated with QT interval prolongation as well as development of the potentially fatal arrhythmia torsades de pointes. The mechanism underlying this cardiotoxicity is high affinity inhibition of the human ether-a-go-go-related gene (hERG) cardiac K+ channel. N-Desmethyl loperamide is the major metabolite of loperamide and is a close structural relative of the parent molecule. To date no information is available regarding the affinity of N-desmethyl loperamide for human cardiac ion channels. The effects of N-desmethyl loperamide on various cloned human cardiac ion channels including hERG, KvLQT1/mink and Na_v1.5 were studied and compared to that of the parent. N-Desmethyl loperamide was a much weaker (7.5-fold) inhibitor of hERG compared to loperamide. However, given the higher plasma levels of the metabolite relative to the parent, it is likely that N-desmethyl loperamide can contribute, at least secondarily, to the cardiotoxicity observed with loperamide abuse. We used the recently solved cryo-EM structure of the hERG channel together with previously published inhibitors, to understand the basis of the interactions as well as the difference that a single methyl plays in the hERG channel blocking affinities of these two compounds.     

Structures of the MthK RCK domain and the effect of Ca2+ on gating ring stability

MthK is a Ca2+-gated K+ channel from Methanobacterium autotrophicum. The crystal structure of the MthK channel in a Ca2+-bound open state was previously determined at 3.3 A and revealed an octameric gating ring composed of eight intracellular ligand-binding RCK (regulate the conductance of K+) domains. It was suggested that Ca2+ binding regulates the gating ring conformation, which in turn leads to the opening and closing of the channel. However, at 3.3 AA resolution, the molecular details of the structure are not well defined, and many of the conclusions drawn from that structure were hypothetical. Here we have presented high resolution structures of the MthK RCK domain with and without Ca2+ bound from three different crystals. These structures revealed a dimeric architecture of the RCK domain and allowed us to visualize the Ca2+ binding and protein-protein contacts at atomic detail. The dimerization of RCK domains is also conserved in other RCK-regulated K+ channels and transporters, suggesting that the RCK dimer serves as a basic unit in the gating ring assembly. A comparison of these dimer structures confirmed that the dimer interface is indeed flexible as suggested previously. However, the conformational change at the flexible interface is of an extent smaller than the previously hypothesized gating ring movement, and a reconstruction of these dimers into octamers by applying protein-protein contacts at the fixed interface did not generate enclosed gating rings. This indicated that there is a high probability that the previously defined fixed interface may not be fixed during channel gating. In addition to the structural studies, we have also carried out biochemical analyses and have shown that near physiological pH, isolated RCK domains form a stable octamer in solution, supporting the notion that the formation of octameric gating ring is a functionally relevant event in MthK gating. Additionally, our stability studies indicated that Ca2+ binding stabilizes the RCK domains in this octameric state.