Conformational dynamics of M2 helices in KirBac channels: helix flexibility in relation to gating via molecular dynamics simulations

KirBac1.1 and 3.1 are bacterial homologues of mammalian inward rectifier K channels. We have performed extended molecular dynamics simulations (five simulations, each of >20 ns duration) of the transmembrane domain of KirBac in two membrane environments, a palmitoyl oleoyl phosphatidylcholine bilayer and an octane slab. Analysis of these simulations has focused on the conformational dynamics of the pore-lining M2 helices, which form the cytoplasmic hydrophobic gate of the channel. Principal components analysis reveals bending of M2, with a molecular hinge at the conserved glycine (Gly134 in KirBac1.1, Gly120 in KirBac3.1). More detailed analysis reveals a dimer-of-dimers type motion. The first two eigenvectors describing the motions of M2 correspond to helix kink and swivel motions. The conformational flexibility of M2 seen in these simulations correlates with differences in M2 conformation between that seen in the X-ray structures of closed channels (KcsA and KirBac) in which the helix is undistorted, and in open channels (e.g. MthK) in which the M2 helix is kinked. Thus, the simulations, albeit on a time scale substantially shorter than that required for channel gating, suggest a gating model in which the intrinsic flexibility of M2 about a molecular hinge is coupled to conformational transitions of an intracellular 'gatekeeper' domain, the latter changing conformation in response to ligand binding.     

Structure and dynamics of K channel pore-lining helices: a comparative simulation study

Isolated pore-lining helices derived from three types of K-channel have been analyzed in terms of their structural and dynamic features in nanosecond molecular dynamics (MD) simulations while spanning a lipid bilayer. The helices were 1) M1 and M2 from the bacterial channel KcsA (Streptomyces lividans), 2) S5 and S6 from the voltage-gated (Kv) channel Shaker (Drosophila melanogaster), and 3) M1 and M2 from the inward rectifier channel Kir6.2 (human). In the case of the Kv and Kir channels, for which x-ray structures are not known, both short and long models of each helix were considered. Each helix was incorporated into a lipid bilayer containing 127 palmitoyloleoylphosphatidylcholine molecules, which was solvated with approximately 4000 water molecules, yielding approximately 20, 000 atoms in each system. Nanosecond MD simulations were used to aid the definition of optimal lengths for the helix models from Kv and Kir. Thus the study corresponds to a total simulation time of 10 ns. The inner pore-lining helices (M2 in KcsA and Kir, S6 in Shaker) appear to be slightly more flexible than the outer pore-lining helices. In particular, the Pro-Val-Pro motif of S6 results in flexibility about a molecular hinge, as was suggested by previous in vacuo simulations (, Biopolymers. 39:503-515). Such flexibility may be related to gating in the corresponding intact channel protein molecules. Analysis of H-bonds revealed interactions with both water and lipid molecules in the water/bilayer interfacial region. Such H-bonding interactions may lock the helices in place in the bilayer during the folding of the channel protein (as is implicit in the two-stage model of membrane protein folding). Aromatic residues at the extremities of the helices underwent complex motions on both short (<10 ps) and long (>100 ps) time scales.     

Conformational dynamics of helix S6 from Shaker potassium channel: simulation studies

Prolines in transmembrane (TM) alpha-helices are believed to play an important structural and/or functional role in membrane proteins. At a structural level a proline residue distorts alpha-helical structure due to the loss of at least one stabilizing backbone hydrogen bond, and introduces flexibility in the helix that may result in substantial kink and swivel motions about the effective "hinge." At a functional level, for example in Kv channels, it is believed that proline-induced molecular hinges may have a direct role in gating, i.e., the conformational change linked to opening/closing the channel to movement of ions. In this article we study the conformational dynamics of the S6 TM helix from of the Kv channel Shaker, which possesses the motif PVP--a motif that is conserved in Kv channels. We perform multiple molecular dynamics simulations of single S6 helices in a membrane-mimetic environment in order to effectively map the kink-swivel conformational space of the protein, exploiting the ability of multiple simulations to achieve greater sampling. We show that the presence of proline locally perturbs the helix, disrupting local dihedral angles and producing local twist and unwinding in the region of the hinge--an effect that is relaxed with distance from the PVP motif. We furthermore show that motions about the hinge are highly anisotropic, reflecting a preferred region of kink-swivel conformation space that may have implications for the gating process.     

A prokaryotic glutamate receptor: homology modelling and molecular dynamics simulations of GluR0

GluR0 is a prokaryotic homologue of mammalian glutamate receptors that forms glutamate-activated, potassium-selective ion channels. The topology of its transmembrane (TM) domain is similar to that of simple potassium channels such as KcsA. Two plausible alignments of the sequence of the TM domain of GluR0 with KcsA are possible, differing in the region of the P helix. We have constructed homology models based on both alignments and evaluated them using 6 ns duration molecular dynamics simulations in a membrane-mimetic environment. One model, in which an insertion in GluR0 relative to KcsA is located in the loop between the M1 and P helices, is preferred on the basis of lower structural drift and maintenance of the P helix conformation during simulation. This model also exhibits inter-subunit salt bridges that help to stabilise the TM domain tetramer. During the simulation, concerted K(+) ion-water movement along the selectivity filter is observed, as is the case in simulations of KcsA. K(+) ion exit from the central cavity is associated with opening of the hydrophobic gate formed by the C-termini of the M2 helices. In the intact receptor the opening of this gate will be controlled by interactions with the extramembranous ligand-binding domains.     

KcsA crystal structure as framework for a molecular model of the Na(+) channel pore

The crystal structure of the pore-forming part of the KcsA bacterial K(+)-selective channel suggests a possible motif for related voltage-gated channels. We examined the hypothesis that the spacial orientation of the KcsA M1 and M2 alpha-helices also predicts the backbone location of S5 and S6 helices of the voltage-gated Na(+) channel. That channel's P region structure is expected to be different because selectivity is determined by side-chain interactions rather than by main-chain carbonyls, and its outer vestibule accommodates relatively large toxin molecules, tetrodotoxin (TTX) and saxitoxin (STX), which interact with selectivity ring residues. The Na(+) channel P loop was well-modeled by the alpha-helix-turn-beta-strand motif, which preserves the relationships for toxin interaction with the Na(+) channel found experimentally. This outer vestibule was docked into the extracellular part of the inverted teepee structure formed by the S5 and S6 helices that were spacially located by coordinates of the KcsA M1 and M2 helix main chains [Doyle et al. (1998) Science 280, 69-74], but populated with side chains of the respective S5 and S6 structures. van der Waals contacts were optimized with minimal adjustment of the S5, S6, and P loop structures, forming a densely packed pore structure. Nonregular external S5-P and P-S6 segments were not modeled here, except the P-S6 segment of domain II. The resulting selectivity region structure is consistent with Na(+) channel permeation properties, offering suggestions for the molecular processes involved in selectivity. The ability to construct a Na(+) channel pore model consistent with most of the available biophysical and mutational information suggests that the KcsA structural framework may be conserved in voltage-gated channels.     

Proline-induced hinges in transmembrane helices: possible roles in ion channel gating

A number of ion channels contain transmembrane (TM) alpha-helices that contain proline-induced molecular hinges. These TM helices include the channel-forming peptide alamethicin (Alm), the S6 helix from voltage-gated potassium (Kv) channels, and the D5 helix from voltage-gated chloride (CLC) channels. For both Alm and KvS6, experimental data implicate hinge-bending motions of the helix in an aspect of channel gating. We have compared the hinge-bending motions of these TM helices in bilayer-like environments by multi-nanosecond MD simulations in an attempt to describe motions of these helices that may underlie possible modes of channel gating. Alm is an alpha-helical channel-forming peptide, which contains a central kink associated with a Gly-x-x-Pro motif in its sequence. Simulations of Alm in a TM orientation for 10 ns in an octane slab indicate that the Gly-x-x-Pro motif acts as a molecular hinge. The S6 helix from Shaker Kv channels contains a Pro-Val-Pro motif. Modeling studies and recent experimental data suggest that the KvS6 helix may be kinked in the vicinity of this motif. Simulations (10 ns) of an isolated KvS6 helix in an octane slab and in a POPC bilayer reveal hinge-bending motions. A pattern-matching approach was used to search for possible hinge-bending motifs in the TM helices of other ion channel proteins. This uncovered a conserved Gly-x-Pro motif in TM helix D5 of CLC channels. MD simulations of a model of hCLC1-D5 spanning an octane slab suggest that this channel also contains a TM helix that undergoes hinge-bending motion. In conclusion, our simulations suggest a model in which hinge-bending motions of TM helices may play a functional role in the gating mechanisms of several different families of ion channels.     

GSK1702934A and M085 directly activate TRPC6 via a mechanism of stimulating the extracellular cavity formed by the pore helix and transmembrane helix S6

Transient receptor potential canonical (TRPC) channels, as important membrane proteins regulating intracellular calcium (Ca²⁺_i) signaling, are involved in a variety of physiological and pathological processes. Activation and regulation of TRPC are more dependent on membrane or intracellular signals. However, how extracellular signals regulate TRPC6 function remains to be further investigated. Here, we suggest that two distinct small molecules, M085 and GSK1702934A, directly activate TRPC6, both through a mechanism of stimulation of extracellular sites formed by the pore helix (PH) and transmembrane (TM) helix S6. In silico docking scanning of TRPC6 identified three extracellular sites that can bind small molecules, of which only mutations on residues of PH and S6 helix significantly reduced the apparent affinity of M085 and GSK1702934A and attenuated the maximal response of TRPC6 to these two chemicals by altering channel gating of TRPC6. Combing metadynamics, molecular dynamics simulations, and mutagenesis, we revealed that W679, E671, E672, and K675 in the PH and N701 and Y704 in the S6 helix constitute an orthosteric site for the recognition of these two agonists. The importance of this site was further confirmed by covalent modification of amino acid residing at the interface of the PH and S6 helix. Given that three structurally distinct agonists M085, GSK1702934A, and AM-0883, act at this site, as well as the occupancy of lipid molecules at this position found in other TRP subfamilies, it is suggested that the cavity formed by the PH and S6 has an important role in the regulation of TRP channel function by extracellular signals.     

A refined atomic scale model of the Saccharomyces cerevisiae K+-translocation protein Trk1p combined with experimental evidence confirms the role of selectivity filter glycines and other key residues

Potassium ion (K⁺) uptake in yeast is mediated mainly by the Trk1/2 proteins that enable cells to survive on external K⁺ concentration as low as a few μM. Fungal Trks are related to prokaryotic TRK and Ktr and plant HKT K⁺ transport systems. Overall sequence similarity is very low, thus requiring experimental verification of homology models. Here a refined structural model of the Saccharomyces cerevisiae Trk1 is presented that was obtained by combining homology modeling, molecular dynamics simulation and experimental verification through functional analysis of mutants. Structural models and experimental results showed that glycines within the selectivity filter, conserved among the K-channel/transporter family, are not only important for protein function, but are also required for correct folding/membrane targeting.A conserved aspartic acid in the P_A helix (D79) and a lysine in the M2_D helix (K1147) were proposed earlier to interact. Our results suggested individual roles of these residues in folding, structural integrity and function. While mutations of D79 completely abolished protein folding, mutations at position 1147 were tolerated to some extent. Intriguingly, a secondary interaction of D79 with R76 could enhance folding/stability of Trk1 and enable a fraction of Trk1[K1147A] to fold.The part of the ion permeation path containing the selectivity filter is shaped similar to that of ion channels. However below the selectivity filter it is obstructed or regulated by a proline containing loop. The presented model could provide the structural basis for addressing the long standing question if Trk1 is a passive or active ion-translocation system.     

Properties and Molecular Determinants of the Natural Flavone Acacetin for Blocking hKv4.3 Channels

The natural flavone acacetin has been demonstrated to inhibit transient outward potassium current (I_to) in human atrial myocytes. However, the molecular determinants of acacetin for blocking I_to are unknown. The present study was designed to investigate the properties and molecular determinants of this compound for blocking hKv4.3 channels (coding I_to) stably expressed in HEK 293 cells using the approaches of whole-cell patch voltage-clamp technique and mutagenesis. It was found that acacetin inhibited hKv4.3 current by binding to both the closed and open channels, and decreased the recovery from inactivation. The blockade of hKv4.3 channels by acacetin was use- and frequency-dependent, and IC₅₀s of acacetin for inhibiting hKv4.3 were 7.9, 6.1, 3.9, and 3.2 µM, respectively, at 0.2, 0.5, 1, and 3.3 Hz. The mutagenesis study revealed that the hKv4.3 mutants T366A and T367A in the P-loop helix, and V392A, I395A and V399A in the S6-segment had a reduced channel blocking efficacy of acacetin (IC₅₀, 44.5 µM for T366A, 25.8 µM for T367A, 17.6 µM for V392A, 16.2 µM for I395A, and 19.1 µM for V399A). These results demonstrate the novel information that acacetin may inhibit the closed channels and block the open state of the channels by binding to their P-loop filter helix and S6 domain. The use- and rate-dependent blocking of hKv4.3 by acacetin is likely beneficial for managing atrial fibrillation.     

Fine-tuning of Voltage Sensitivity of the Kv1.2 Potassium Channel by Interhelix Loop Dynamics

Many proteins function by changing conformation in response to ligand binding or changes in other factors in their environment. Any change in the sequence of a protein, for example during evolution, which alters the relative free energies of the different functional conformations changes the conditions under which the protein will function. Voltage-gated ion channels are membrane proteins that open and close an ion-selective pore in response to changes in transmembrane voltage. The charged S4 transmembrane helix transduces changes in transmembrane voltage into a change in protein internal energy by interacting with the rest of the channel protein through a combination of non-covalent interactions between adjacent helices and covalent interactions along the peptide backbone. However, the structural basis for the wide variation in the V₅₀ value between different voltage-gated potassium channels is not well defined. To test the role of the loop linking the S3 helix and the S4 helix in voltage sensitivity, we have constructed a set of mutants of the rat K_v1.2 channel that vary solely in the length and composition of the extracellular loop that connects S4 to S3. We evaluated the effect of these different loop substitutions on the voltage sensitivity of the channel and compared these experimental results with molecular dynamics simulations of the loop structures. Here, we show that this loop has a significant role in setting the precise V₅₀ of activation in K_v1 family channels.     

Structure and dynamics of the pore-lining helix of the nicotinic receptor: MD simulations in water, lipid bilayers, and transbilayer bundles

Multiple nanosecond duration molecular dynamics simulations on the pore-lining M2 helix of the nicotinic acetylcholine receptor reveal how its structure and dynamics change as a function of environment. In water, the M2 helix partially unfolds to form a molecular hinge in the vicinity of a central Leu residue that has been implicated in the mechanism of ion channel gating. In a phospholipid bilayer, either as a single transmembrane helix, or as part of a pentameric helix bundle, the M2 helix shows less flexibility, but still exhibits a kink in the vicinity of the central Leu. The single M2 helix tilts relative to the bilayer normal by 12 degrees, in agreement with recent solid state NMR data (Opella et al., Nat Struct Biol 6:374-379, 1999). The pentameric helix bundle, a model for the pore domain of the nicotinic receptor and for channels formed by M2 peptides in a bilayer, is remarkably stable over a 2-ns MD simulation in a bilayer, provided one adjusts the pK(A)s of ionizable residues to their calculated values (when taking their environment into account) before starting the simulation. The resultant transbilayer pore shows fluctuations at either mouth which transiently close the channel. Proteins 2000;39:47-55.     

Decomposition of Slide Helix Contributions to ATP-dependent Inhibition of Kir6.2 Channels

Regulation of inwardly rectifying potassium channels by intracellular ligands couples cell membrane excitability to important signaling cascades and metabolic pathways. We investigated the molecular mechanisms that link ligand binding to the channel gate in ATP-sensitive Kir6.2 channels. In these channels, the “slide helix” forms an interface between the cytoplasmic (ligand-binding) domain and the transmembrane pore, and many slide helix mutations cause loss of function. Using a novel approach to rescue electrically silent channels, we decomposed the contribution of each interface residue to ATP-dependent gating. We demonstrate that effective inhibition by ATP relies on an essential aspartate at residue 58. Characterization of the functional importance of this conserved aspartate, relative to other residues in the slide helix, has been impossible because of loss-of-function of Asp-58 mutant channels. The Asp-58 position exhibits an extremely stringent requirement for aspartate because even a highly conservative mutation to glutamate is insufficient to restore normal channel function. These findings reveal unrecognized slide helix elements that are required for functional channel expression and control of Kir6.2 gating by intracellular ATP.     

Effective diffusion coefficients of K+ and Cl- ions in ion channel models.

We have used molecular dynamics simulations, corresponding to a total simulation time of 11 ns, to investigate the effective short-time local diffusion coefficient of potassium and chloride ions in a series of model ion channels. These models, which include channels formed by the fungal peptide alamethicin, by a synthetic leucine-serine peptide, and by the pore-lining M2 helix bundle of the nicotinic acetylcholine receptor, have a range of different secondary structures, diameters and hydrophobicities. We find that the diffusion coefficients of both ions are appreciably reduced in the narrower channels, the extent of the reduction being similar for both the anionic and cationic species. This suggests that a difference in mobility cannot be the source of the ion selectivity exhibited by some of the channels (for example, the leucine-serine peptide). We find no evidence for a reduction in mobility of either ion in the nAChR model. These results are broadly in line with a previous similar study of Na+ ions, and may be useful in Poisson-Nernst-Planck, Eyring rate theory or Brownian dynamics calculations of channel conductance.     

A Kv channel with an altered activation gate sequence displays both "fast" and "slow" activation kinetics

The Kv1–4 families of K⁺ channels contain a tandem proline motif (PXP) in the S6 helix that is crucial for channel gating. In human Kv1.5, replacing the first proline by an alanine resulted in a nonfunctional channel. This mutant was rescued by introducing another proline at a nearby position, changing the sequence into AVPP. This resulted in a channel that activated quickly (ms range) upon the first depolarization. However, thereafter, the channel became trapped in another gating mode that was characterized by slow activation kinetics (s range) with a shallow voltage dependence. The switch in gating mode was observed even with very short depolarization steps, but recovery to the initial "fast" mode was extremely slow. Computational modeling suggested that switching occurred during channel deactivation. To test the effect of the altered PXP sequence on the mobility of the S6 helix, we used molecular dynamics simulations of the isolated S6 domain of wild type (WT) and mutants starting from either a closed or open conformation. The WT S6 helix displayed movements around the PXP region with simulations starting from either state. However, the S6 with a AVPP sequence displayed flexibility only when started from the closed conformation and was rigid when started from the open state. These results indicate that the region around the PXP motif may serve as a "hinge" and that changing the sequence to AVPP results in channels that deactivate to a state with an alternate configuration that renders them "reluctant" to open subsequently. voltage-gated potassium channel     

Transmembrane structure of an inwardly rectifying potassium channel

Inwardly rectifying potassium channels (K(ir)), comprising four subunits each with two transmembrane domains, M1 and M2, regulate many important physiological processes. We employed a yeast genetic screen to identify functional channels from libraries of K(ir) 2.1 containing mutagenized M1 or M2 domains. Patterns in the allowed sequences indicate that M1 and M2 are helices. Protein-lipid and protein-water interaction surfaces identified by the patterns were verified by sequence minimization experiments. Second-site suppressor analyses of helix packing indicate that the M2 pore-lining inner helices are surrounded by the M1 lipid-facing outer helices, arranged such that the M1 helices participate in subunit-subunit interactions. This arrangement is distinctly different from the structure of a bacterial potassium channel with the same topology and identifies helix-packing residues as hallmark sequences common to all K(ir) superfamily members.     

Thermal‐ and urea‐induced unfolding processes of glutathione S‐transferase by molecular dynamics simulation

The Schistosoma juponicum 26 kDa glutathione S‐transferase (sj26GST) consists of the N‐terminal domain (N‐domain), containing three alpha‐helices (named H1‐H3) and four anti‐parallel beta‐strands (S1‐S4), and the C‐terminal domain (C‐domain), comprising five alpha‐helices (named H4‐H8). In present work, molecular dynamics simulations and fluorescence spectroscopic were used to gain insights into the unfolding process of sj26GST. The molecular dynamics simulations on sj26GST subunit both in water and in 8 M urea were carried out at 300 K, 400 K and 500 K, respectively. Spectroscopic measurements were employed to monitor structural changes. Molecular dynamics simulations of sj26GST subunit induced by urea and temperature showed that the initial unfolding step of sj26GST both in water and urea occurred on N‐domain, involving the disruption of helices H2, H3 and strands S3 and S4, whereas H6 was the last region exposed to solution and was the last helix to unfold. Moreover, simulations analyses combining with fluorescence and circular dichroism spectra indicated that N‐domain could not fold independent, suggesting that correct folding of N‐domain depended on its interactions with C‐domain. We further proposed that the folding of GSTs could begin with the hydrophobic collapse of C‐domain whose H4, H5, H6 and H7 could move close to each other and form a hydrophobic core, especially H6 wrapped in the hydrophobic center and beginning spontaneous formation of the helix. S3, S4, H3, and H2 could form in the wake of the interaction between C‐domain and N‐domain. The paper can offer insights into the molecular mechanism of GSTs unfolding. © 2014 Wiley Periodicals, Inc. Biopolymers 103: 247–259, 2015.     

Exploring coumarin egress channels in human cytochrome P450 2A6 by random acceleration and steered molecular dynamics simulations

The kinetic analysis of coumarin oxidation by CYP2A6 suggested that substrate binding and release occurred in the multiple steps and such events proceeded rapidly. However, the crystal structure of the CYP2A6-coumarin complex reveals that no obvious channel is open enough to allow coumarin to pass through. Thus, an intriguing and important question arises: how coumarin enters and exits the active site, which is deeply buried at the center of CYP2A6 fold. In this study, geometric analysis of the potential openings was first performed on all the available crystal structures of CYP2A6. And then, random acceleration molecular dynamics simulations were used to explore the possible substrate egress channels in CYP2A6. Two channels were most frequently observed. Afterwards, steered molecular dynamics simulations were performed and potentials of mean force were constructed to compare the preference of the two channels serving as the substrate egress channel. The results showed that channel 2c, which is located between helices I and G and the helix B'-C region, was the most likely channel for coumarin egress. The opening of channel 2c was characterized by a rotation of Phe111 together with a bending of helix B'. Our findings will not only be helpful for understanding the unbinding mechanism of coumarin and for identifying structural determinants related to the biological function of CYP2A6, but also provide further insight into the channel selectivity of P450s.     

Inactivation and Secondary Structure in the D4/S4-5 Region of the SkM1 Sodium Channel

The D4/S4-5 interhelical region plays a role in sodium channel fast inactivation. Examination of S4-5 primary structure in all domains suggests a possible amphipathic helical conformation in which a conserved group of small hydrophobic residues occupies one contiguous surface with a more variable complement of nonpolar and polar residues on the opposite face. We evaluated this potential structure by replacing each residue in D4/S4-5 of the rat SkM1 skeletal muscle sodium channel with substitutions having different side chain properties. Of the 63 mutations analyzed, 44 produced functional channels. P1473 was intolerant of substitutions. Nonpolar substitutions in the conserved hydrophobic region were functionally similar to wild type, while charged mutations in this region before P1473 were nonfunctional. Charged mutations at F1466, M1469, M1470, and A1474, located on the opposite surface of the predicted helix, produced functional channels with pronounced slowing of inactivation, shifted voltage dependence of steady-state inactivation, and increased rate of recovery from inactivation. The substituted-cysteine-accessibility method was used to probe accessibility at each position. Residues L1465, F1466, A1467, M1469, M1470, L1472, A1474, and F1476C were easily accessible for modification by sulfhydryl reagents; L1464, L1468, S1471, and L1475 were not accessible within the time frame of our measurements. Molecular dynamics simulations of residues A1458 to N1477 were then used to explore energetically favorable local structures. Based on mutagenesis, substituted-cysteine-accessibility method, and modeling results, we suggest a secondary structure for the D4/S4-5 region in which the peptide chain is α-helical proximal to P1473, bends at this residue, and may continue beyond this point as a random coil. In this configuration, the entire resultant loop is amphipathic; four residues on one surface could form part of the binding site for the inactivation particle.     

Molecular dynamics simulations of a K+ channel blocker: Tc1 toxin from Tityus cambridgei

Toxins that block voltage-gated potassium (Kv) channels provide a possible template for improved homology models of the Kv pore. In assessing the interactions of Kv channels and their toxins it is important to determine the dynamic flexibility of the toxins. Multiple 10 ns duration molecular dynamics simulations combined with essential dynamics analysis have been used to explore the flexibility of four different Kv channel-blocking toxins. Three toxins (Tc1, AgTx and ChTx) share a common fold. They also share a common pattern of conformational dynamics, as revealed by essential dynamics analysis of the simulation results. This suggests that some aspects of dynamic behaviour are conserved across a single protein fold class. In each of these three toxins, the residue exhibiting minimum flexibility corresponds to a conserved lysine residue that is suggested to interact with the filter domain of the channel. Thus, comparative simulations reveal functionally important conservation of molecular dynamics as well as protein fold across a family of related toxins.     

Tuning of EAG K+ channel inactivation: Molecular determinants of amplification by mutations and a small molecule

Ether-à-go-go (EAG) and EAG-related gene (ERG) K(+) channels are close homologues but differ markedly in their gating properties. ERG1 channels are characterized by rapid and extensive C-type inactivation, whereas mammalian EAG1 channels were previously considered noninactivating. Here, we show that human EAG1 channels exhibit an intrinsic voltage-dependent slow inactivation that is markedly enhanced in rate and extent by 1-10 μM 3-nitro-N-(4-phenoxyphenyl) benzamide, or ICA105574 (ICA). This compound was previously reported to have the opposite effect on ERG1 channels, causing an increase in current magnitude by inhibition of C-type inactivation. The voltage dependence of 2 μM ICA-induced inhibition of EAG1 current was half-maximal at -73 mV, 62 mV negative to the half-point for channel activation. This finding suggests that current inhibition by the drug is mediated by enhanced inactivation and not open-channel block, where the voltage half-points for current inhibition and channel activation are predicted to overlap, as we demonstrate for clofilium and astemizole. The mutation Y464A in the S6 segment also induced inactivation of EAG1, with a time course and voltage dependence similar to that caused by 2 μM ICA. Several Markov models were investigated to describe gating effects induced by multiple concentrations of the drug and the Y464A mutation. Models with the smallest fit error required both closed- and open-state inactivation. Unlike typical C-type inactivation, the rate of Y464A- and ICA-induced inactivation was not decreased by external tetraethylammonium or elevated [K(+)](e). EAG1 channel inactivation introduced by Y464A was prevented by additional mutation of a nearby residue located in the S5 segment (F359A) or pore helix (L434A), suggesting a tripartite molecular model where interactions between single residues in S5, S6, and the pore helix modulate inactivation of EAG1 channels.