Interaction of anesthetics with open and closed conformations of a potassium channel studied via molecular dynamics and normal mode analysis

A variety of experiments suggest that membrane proteins are important targets of anesthetic molecules, and that ion channels interact differently with anesthetics in their open and closed conformations. The availability of an open and a closed structural model for the KirBac1.1 potassium channel has made it possible to perform a comparative analysis of the interactions of anesthetics with the same channel in its open and closed states. To this end, all-atom molecular dynamics simulations supplemented by normal mode analysis have been employed to probe the interactions of the inhalational anesthetic halothane with both an open and closed conformer of KirBac1.1 embedded in a lipid bilayer. Normal mode analysis on the closed and open channel, in the presence and absence of halothane, reveals that the anesthetic modulates the global as well as the local dynamics of both conformations differently. In the case of the open channel, the observed reduction of flexibility of residues in the inner helices suggests a functional modification action of anesthetics on ion channels. In this context, preferential quenching of the aromatic residue motion and modulation of global dynamics by halothane may be seen as steps toward potentiating or favoring open state conformations. These molecular dynamics simulations provide the first insights into possible specific interactions between anesthetic molecules and ion channels in different conformations.     

Two different conformational states of the KirBac3.1 potassium channel revealed by electron crystallography

Potassium channels allow the selective flow of K(+) ions across membranes. In response to external gating signals, the potassium channel can move reversibly through a series of structural conformations from a closed to an open state. 2D crystals of the inwardly rectifying K(+) channel KirBac3.1 from Magnetospirillum magnetotacticum have been captured in two distinct conformations, providing "snap shots" of the gating process. Analysis by electron cryomicroscopy of these KirBac3.1 crystals has resulted in reconstructed images in projection at 9 A resolution. Kir channels are tetramers of four subunits arranged as dimers of dimers. Each subunit has two transmembrane helices (inner and outer). In one crystal form, the pore is blocked; in the other crystal form, the pore appears open. Modeling based on the KirBac1.1 (closed) crystal structure shows that opening of the ion conduction pathway could be achieved by bending of the inner helices and significant movements of the outer helices.     

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.     

Structural locus of the pH gate in the Kir1.1 inward rectifier channel

The closed-state crystal structure of prokaryotic inward rectifier, KirBac1.1, has implicated four inner helical phenylalanines near the cytoplasmic side as a possible locus of the channel gate. In the present study, we investigate whether this structural feature corresponds to the physiological pH gate of the renal inward rectifier, Kir1.1 (ROMK, KCNJ1). Kir1.1 is endogenous to the mammalian renal collecting duct and the thick ascending limb of Henle and is strongly gated by internal pH in the physiological range. It has four leucines (L160-Kir1.1b), homologous to the phenylalanines of KirBac1.1, which could function as steric gates near the convergence of the inner (M2) helices. Replacing these Leu-160 residues of Kir1.1b by smaller glycines abolished pH gating; however, replacement with alanines, whose side chains are intermediate in size between leucine and glycine, did not eliminate normal pH gating. Furthermore, a double mutant, constructed by adding the I163M-Kir1.1b mutation to the L160G mutation, also lacked normal pH gating, although the I163M mutation by itself enhanced the pH sensitivity of the channel. In addition to size, side-chain hydrophobicity at 160-Kir1.1b was also important for normal pH gating. Mutants with polar side chains (L160S, L160T) did not gate normally and were as insensitive to internal pH as the L160G mutant. Hence, either small or highly polar side chains at 160-Kir1.1b stabilize the open state of the channel. A homology model of the Kir1.1 closed state, based on the crystal structure of KirBac1.1, was consistent with our electrophysiological data and implies that closure of the Kir1.1 pH gate results from steric occlusion of the permeation path by the convergence of four leucines at the cytoplasmic apex of the inner transmembrane helices. In the open state, K crosses the pH gate together with its hydration shell.     

Role of aromatic localization in the gating process of a potassium channel

Position of the transmembrane aromatic residues of the KirBac1.1 potassium channel shifts from an even distribution in the closed state toward the membrane/solute interface in the open state model. This is the first example of an integral membrane protein making use of the observed preference for transmembrane aromatic residues to reside at the interfaces. The process of aromatic localization is proposed as a means of directing and stabilizing structural changes during conformational transitions within the transmembrane region of integral membrane proteins. All-atom molecular dynamics simulations of the open and closed conformers in a membrane environment have been carried out to take account of the interactions between the aromatic residues and the lipids, which may be involved in the conformational change, e.g., the gating of the channel.     

Role of Aromatic Localization in the Gating Process of a Potassium Channel

Position of the transmembrane aromatic residues of the KirBac1.1 potassium channel shifts from an even distribution in the closed state toward the membrane/solute interface in the open state model. This is the first example of an integral membrane protein making use of the observed preference for transmembrane aromatic residues to reside at the interfaces. The process of aromatic localization is proposed as a means of directing and stabilizing structural changes during conformational transitions within the transmembrane region of integral membrane proteins. All-atom molecular dynamics simulations of the open and closed conformers in a membrane environment have been carried out to take account of the interactions between the aromatic residues and the lipids, which may be involved in the conformational change, e.g., the gating of the channel.     

A homology model of the pore region of HCN channels

HCN channels are activated by membrane hyperpolarization and regulated by cyclic nucleotides, such as cyclic adenosine-mono-phosphate (cAMP). Here we present structural models of the pore region of these channels obtained by using homology modeling and validated against spatial constraints derived from electrophysiological experiments. For the construction of the models we make two major assumptions, justified by electrophysiological observations: i), in the closed state, the topology of the inner pore of HCN channels is similar to that of K(+) channels. In particular, the orientation of the S5 and S6 helices of HCN channels is very similar to that of the corresponding helices of the K(+) KcsA and K(+) KirBac1.1 channels. Thus, we use as templates the x-ray structure of these K(+) channels. ii), In the open state, the S6 helix is bent further than it is in the closed state, as suggested (but not proven) by experimental data. For this reason, the template of the open conformation is the x-ray structure of the MthK channel. The structural models of the closed state turn out to be consistent with all the available electrophysiological data. The model of the open state turned out to be consistent with all the available electrophysiological data in the filter region, including additional experimental data performed in this work. However, it required the introduction of an appropriate, experimentally derived constraint for the S6 helix. Our modeling provides a structural framework for understanding several functional properties of HCN channels: i), the cysteine ring at the inner mouth of the pore may act as a sensor of the intracellular oxidizing/reducing conditions; ii), the bending amplitude of the S6 helix upon gating appears to be significantly smaller than that found in MthK channels; iii), the reduced ionic selectivity of HCN channels, relative to that of K(+) channels, may be caused, at least in part, by the larger flexibility of the inner pore of HCN channels.     

Filter flexibility and distortion in a bacterial inward rectifier K+ channel: simulation studies of KirBac1.1

The bacterial channel KirBac1.1 provides a structural homolog of mammalian inward rectifier potassium (Kir) channels. The conformational dynamics of the selectivity filter of Kir channels are of some interest in the context of possible permeation and gating mechanisms for this channel. Molecular dynamics simulations of KirBac have been performed on a 10-ns timescale, i.e., comparable to that of ion permeation. The results of five simulations (total simulation time 50 ns) based on three different initial ion configurations and two different model membranes are reported. These simulation data provide evidence for limited (<0.1 nm) filter flexibility during the concerted motion of ions and water molecules within the filter, such local changes in conformation occurring on an approximately 1-ns timescale. In the absence of K(+) ions, the KirBac selectivity filter undergoes more substantial distortions. These resemble those seen in comparable simulations of other channels (e.g., KcsA and KcsA-based homology models) and are likely to lead to functional closure of the channel. This suggests filter distortions may provide a mechanism of K-channel gating in addition to changes in the hydrophobic gate formed at the intracellular crossing point of the M2 helices. The simulation data also provide evidence for interactions of the "slide" (pre-M1) helix of KirBac with phospholipid headgroups.     

Molecular dynamics simulations of inwardly rectifying (Kir) potassium channels: a comparative study

Inward rectifier potassium (Kir) channels regulate cell excitability and transport K+ ions across membranes. Homotetrameric models of three mammalian Kir channels (Kir1.1, Kir3.1, and Kir6.2) have been generated, using the KirBac3.1 transmembrane and rat Kir3.1 intracellular domain structures as templates. All three models have been explored by 10 ns molecular dynamics simulations in phospholipid bilayers. Analysis of the initial structures revealed conservation of potential lipid interaction residues (Trp/Tyr and Arg/Lys side chains near the lipid headgroup-water interfaces). Examination of the intracellular domains revealed key structural differences between Kir1.1 and Kir6.2 which may explain the difference in channel inhibition by ATP. The behavior of all three models in the MD simulations revealed that they have conformational stability similar to that seen for comparable simulations of, for example, structures derived from cryoelectron microscopy data. Local distortions of the selectivity filter were seen during the simulations, as observed in previous simulations of KirBac and in simulations and structures of KcsA. These may be related to filter gating of the channel. The intracellular hydrophobic gate does not undergo any substantial changes during the simulations and thus remains functionally closed. Analysis of lipid-protein interactions of the Kir models emphasizes the key role of the M0 (or "slide") helix which lies approximately parallel to the bilayer-water interface and forms a link between the transmembrane and intracellular domains of the channel.     

Functional Complementation and Genetic Deletion Studies of KirBac Channels: ACTIVATORY MUTATIONS HIGHLIGHT GATING-SENSITIVE DOMAINS*

The superfamily of prokaryotic inwardly rectifying (KirBac) potassium channels is homologous to mammalian Kir channels. However, relatively little is known about their regulation or about their physiological role in vivo. In this study, we have used random mutagenesis and genetic complementation in K⁺-auxotrophic Escherichia coli and Saccharomyces cerevisiae to identify activatory mutations in a range of different KirBac channels. We also show that the KirBac6.1 gene (slr5078) is necessary for normal growth of the cyanobacterium Synechocystis PCC6803. Functional analysis and molecular dynamics simulations of selected activatory mutations identified regions within the slide helix, transmembrane helices, and C terminus that function as important regulators of KirBac channel activity, as well as a region close to the selectivity filter of KirBac3.1 that may have an effect on gating. In particular, the mutations identified in TM2 favor a model of KirBac channel gating in which opening of the pore at the helix-bundle crossing plays a far more important role than has recently been proposed.     

An NMDA Receptor Gating Mechanism Developed from MD Simulations Reveals Molecular Details Underlying Subunit-Specific Contributions

N-methyl-D-aspartate (NMDA) receptors are obligate heterotetrameric ligand-gated ion channels that play critical roles in learning and memory. Here, using targeted molecular dynamics simulations, we developed an atomistic model for the gating of the GluN1/GluN2A NMDA receptor. Upon agonist binding, lobe closure of the ligand-binding domain produced outward pulling of the M3-D2 linkers, leading to outward movements of the C-termini of the pore-lining M3 helices and opening of the channel. The GluN2A subunits, similar to the distal (B/D) subunits in the homotetrameric GluA2 α-amino-3-hydroxy-5-methyl-4-isoxazoleproprionate receptor, had greater M3 outward movements and thus contributed more to channel gating than the GluN1 subunits. Our gating model is validated by functional studies, including cysteine modification data indicating wider accessibility to the GluN1 M3 helices than to the GluN2A M3 helices from the lumen of the open channel, and reveals why the Lurcher mutation in GluN1 has a stronger ability in maintaining channel opening than the counterpart in GluN2A. The resulting structural model for the open state provides an explanation for the Ca²⁺ permeability of NMDA receptors, and the structural differences between the closed and open states form the basis for drug design.     

New insights of membrane environment effects on MscL channel mechanics from theoretical approaches

The prokaryotic mechanosensitive channel of large conductance (MscL) is a remarkable integral membrane protein. During hypo-osmotic shock, it responses to membrane tension through large conformational changes, that lead to an open state of the pore. The structure of the channel from Mycobacterium tuberculosis has been resolved in the closed state. Numerous experiments have attempted to trap the channel in its open state but they did not succeed in obtaining a structure. A gating mechanism has been proposed based on different experimental data but there is no experimental technique available to follow this process in atomic details. In addition, it has been shown that a decrease of the lipid bilayer thickness lowered MscL activation energy and stabilized a structurally distinct closed channel intermediate. Here, we use atomistic molecular dynamics simulations to investigate the effect of the lipid bilayer thinning on our model of the structure of the Escherichia coli. We thoroughly analyze simulations of the channel embedded in two pre-equilibrated membranes differing by their hydrophobic tail length (DMPE and POPE). The MscL structure remains stable in POPE, whereas a distinct structural state is obtained in DMPE in response to hydrophobic mismatch. This latter is obtained by tilts and kinks of the transmembrane helices, leading to a widening and a diminution of the channel height. Part of these motions is guided by a competition between solvent and lipids for the interaction with the periplasmic loops. We finally conduct a principal component analysis of the simulation and compare anharmonic motions with harmonic ones, previously obtained from a coarse-grained normal mode analysis performed on the same structural model. Significant similarities exist between low-frequency harmonic motions and those observed with essential dynamics in DMPE. In summary, change in membrane thickness permits to accelerate the conformational changes involved in the mechanics of the E. coli channel, providing a closed structural intermediate en route to the open state. These results give clues for better understanding why the channel activation energy is lowered in a thinner membrane.     

Closing in on the resting state of the Shaker K(+) channel

Membrane depolarization causes voltage-gated ion channels to transition from a resting/closed conformation to an activated/open conformation. We used voltage-clamp fluorometry to measure protein motion at specific regions of the Shaker Kv channel. This enabled us to construct new structural models of the resting/closed and activated/open states based on the Kv1.2 crystal structure using the Rosetta-Membrane method and molecular dynamics simulations. Our models account for the measured gating charge displacement and suggest a molecular mechanism of activation in which the primary voltage sensors, S4s, rotate by approximately 180 degrees as they move "outward" by 6-8 A. A subsequent tilting motion of the S4s and the pore domain helices, S5s, of all four subunits induces a concerted movement of the channel's S4-S5 linkers and S6 helices, allowing ion conduction. Our models are compatible with a wide body of data and resolve apparent contradictions that previously led to several distinct models of voltage sensing.     

The gating mechanism of the bacterial mechanosensitive channel MscL revealed by molecular dynamics simulations: From tension sensing to channel opening

One of the ultimate goals of the study on mechanosensitive (MS) channels is to understand the biophysical mechanisms of how the MS channel protein senses forces and how the sensed force induces channel gating. The bacterial MS channel MscL is an ideal subject to reach this goal owing to its resolved 3D protein structure in the closed state on the atomic scale and large amounts of electrophysiological data on its gating kinetics. However, the structural basis of the dynamic process from the closed to open states in MscL is not fully understood. In this study, we performed molecular dynamics (MD) simulations on the initial process of MscL opening in response to a tension increase in the lipid bilayer. To identify the tension-sensing site(s) in the channel protein, we calculated interaction energy between membrane lipids and candidate amino acids (AAs) facing the lipids. We found that Phe78 has a conspicuous interaction with the lipids, suggesting that Phe78 is the primary tension sensor of MscL. Increased membrane tension by membrane stretch dragged radially the inner (TM1) and outer (TM2) helices of MscL at Phe78, and the force was transmitted to the pentagon-shaped gate that is formed by the crossing of the neighboring TM1 helices in the inner leaflet of the bilayer. The radial dragging force induced radial sliding of the crossing portions, leading to a gate expansion. Calculated energy for this expansion is comparable to an experimentally estimated energy difference between the closed and the first subconductance state, suggesting that our model simulates the initial step toward the full opening of MscL. The model also successfully mimicked the behaviors of a gain of function mutant (G22N) and a loss of function mutant (F78N), strongly supporting that our MD model did simulate some essential biophysical aspects of the mechano-gating in MscL.     

The gating mechanism of the bacterial mechanosensitive channel MscL revealed by molecular dynamics simulations

One of the ultimate goals of the study on mechanosensitive (MS) channels is to understand the biophysical mechanisms of how the MS channel protein senses forces and how the sensed force induces channel gating. The bacterial MS channel MscL is an ideal subject to reach this goal owing to its resolved 3D protein structure in the closed state on the atomic scale and large amounts of electrophysiological data on its gating kinetics. However, the structural basis of the dynamic process from the closed to open states in MscL is not fully understood. In this study, we performed molecular dynamics (MD) simulations on the initial process of MscL opening in response to a tension increase in the lipid bilayer. To identify the tension-sensing site(s) in the channel protein, we calculated interaction energy between membrane lipids and candidate amino acids (AAs) facing the lipids. We found that Phe78 has a conspicuous interaction with the lipids, suggesting that Phe78 is the primary tension sensor of MscL. Increased membrane tension by membrane stretch dragged radially the inner (TM1) and outer (TM2) helices of MscL at Phe78, and the force was transmitted to the pentagon-shaped gate that is formed by the crossing of the neighboring TM1 helices in the inner leaflet of the bilayer. The radial dragging force induced radial sliding of the crossing portions, leading to a gate expansion. Calculated energy for this expansion is comparable to an experimentally estimated energy difference between the closed and the first subconductance state, suggesting that our model simulates the initial step toward the full opening of MscL. The model also successfully mimicked the behaviors of a gain of function mutant (G22N) and a loss of function mutant (F78N), strongly supporting that our MD model did simulate some essential biophysical aspects of the mechano-gating in MscL.     

Conformational dynamics of the ligand-binding domain of inward rectifier K channels as revealed by molecular dynamics simulations: toward an understanding of Kir channel gating

Inward rectifier (Kir) potassium channels are characterized by two transmembrane helices per subunit, plus an intracellular C-terminal domain that controls channel gating in response to changes in concentration of various ligands. Based on the crystal structure of the tetrameric C-terminal domain of Kir3.1, it is possible to build a homology model of the ATP-binding C-terminal domain of Kir6.2. Molecular dynamics simulations have been used to probe the dynamics of Kir C-terminal domains and to explore the relationship between their dynamics and possible mechanisms of channel gating. Multiple simulations, each of 10 ns duration, have been performed for Kir3.1 (crystal structure) and Kir6.2 (homology model), in both their monomeric and tetrameric forms. The Kir6.2 simulations were performed with and without bound ATP. The results of the simulations reveal comparable conformational stability for the crystal structure and the homology model. There is some decrease in conformational flexibility when comparing the monomers with the tetramers, corresponding mainly to the subunit interfaces in the tetramer. The beta-phosphate of ATP interacts with the side chain of K185 in the Kir6.2 model and simulations. The flexibility of the Kir6.2 tetramer is not changed greatly by the presence of bound ATP, other than in two loop regions. Principal components analysis of the simulated dynamics suggests loss of symmetry in both the Kir3.1 and Kir6.2 tetramers, consistent with "dimer-of-dimers" motion of subunits in C-terminal domains of the corresponding Kir channels. This is suggestive of a gating model in which a transition between exact tetrameric symmetry and dimer-of-dimers symmetry is associated with a change in transmembrane helix packing coupled to gating of the channel. Dimer-of-dimers motion of the C-terminal domain tetramer is also supported by coarse-grained (anisotropic network model) calculations. It is of interest that loss of exact rotational symmetry has also been suggested to play a role in gating in the bacterial Kir homolog, KirBac1.1, and in the nicotinic acetylcholine receptor channel.     

Moving the pH Gate of the Kir1.1 Inward Rectifier Channel

Both structural and functional studies suggest that pH gating of the inward rectifier potassium (K) channel, Kir1.1 (ROMK), is mediated by the convergence of 4 hydrophobic leucines (one from each subunit) near the cytoplasmic bundle-crossing of the inner transmembrane helices. We tested this hypothesis by moving the putative leucine gate from the L160-Kir1.1b to other positions along the inner transmembrane helix, and measuring inward current and conductance as functions of internal pH, using the Xenopus oocyte heterologous expression system. Results of these studies indicated that it was possible to replace the putative inward rectifier pH gate at L160-Kir1.1b by either a leucine or methionine at 157-Kir1.1b (G157L-L160G or G157M-L160G). Although both leucine and methionine gated the channel at 157-Kir1.1b, residues of similar hydrophobicity (tyrosine and valine) did not. Hence, hydrophobicity was a necessary but not a sufficient condition for steric gating at 157. This was in contrast to the 160-Kir1.1b locus, where side-chain hydrophobicity was both a necessary and sufficient property for steric gating. Homology models were constructed for all mutants that expressed significant whole-cell currents, using the closed-state coordinates of the prokaryotic inward rectifier, KirBac1.1. Models of mutants that retained pH gating were too narrow at the bundle crossing to permit hydrated K ion permeation in the closed-state. On the other hand, mutants that lost pH gating had ample space at the bundle crossing for hydrated K permeation in the closed-state. These results support our hypothesis that hydrophobic leucines at the cytoplasmic end of the inner transmembrane helices comprise the principal pH gate of Kir1.1, a gate that can be relocated from 160-Kir1.1b to 157-Kir1.1b.     

Voltage-dependent conformational changes in connexin channels

Channels formed by connexins display two distinct types of voltage-dependent gating, termed V(j)- or fast-gating and loop- or slow-gating. Recent studies, using metal bridge formation and chemical cross-linking have identified a region within the channel pore that contributes to the formation of the loop-gate permeability barrier. The conformational changes are remarkably large, reducing the channel pore diameter from 15 to 20? to less than 4?. Surprisingly, the largest conformational change occurs in the most stable region of the channel pore, the 3(10) or parahelix formed by amino acids in the 42-51 segment. The data provide a set of positional constraints that can be used to model the structure of the loop-gate closed state. Less is known about the conformation of the V(j)-gate closed state. There appear to be two different mechanisms; one in which conformational changes in channel structure are linked to a voltage sensor contained in the N-terminus of Cx26 and Cx32 and a second in which the C-terminus of Cx43 and Cx40 may act either as a gating particle to block the channel pore or alternatively to stabilize the closed state. The later mechanism utilizes the same domains as implicated in effecting pH gating of Cx43 channels. It is unclear if the two V(j)-gating mechanisms are related or if they represent different gating mechanisms that operate separately in different subsets of connexin channels. A model of the V(j)-closed state of Cx26 hemichannel that is based on the X-ray structure of Cx26 and electron crystallographic structures of a Cx26 mutation suggests that the permeability barrier for V(j)-gating is formed exclusively by the N-terminus, but recent information suggests that this conformation may not represent a voltage-closed state. Closed state models are considered from a thermodynamic perspective based on information from the 3.5? Cx26 crystal structure and molecular dynamics (MD) simulations. The applications of computational and experimental methods to define the path of allosteric molecular transitions that link the open and closed states are discussed. This article is part of a Special Issue entitled: The Communicating junctions, composition, structure and characteristics.     

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.     

Gating of MscL studied by steered molecular dynamics

Steered molecular dynamics simulations of the mechanosensitive channel of large conductance, MscL, were used to investigate how forces arising from membrane tension induce gating of the channel. A homology model of the closed form of MscL from Escherichia coli was subjected to external forces of 35-70 pN applied to residues near the membrane-water interface. The magnitude and location of these forces corresponded to those determined from the lateral pressure profile computed from a lipid bilayer simulation. A fully expanded state was obtained on the 10-ns timescale that revealed the mechanism for transducing membrane forces into channel opening. The expanded state agrees well with proposed models of MscL gating, in that it entails an irislike expansion of the pore accompanied by tilting of the transmembrane helices. The channel was most easily opened when force was applied predominantly on the cytoplasmic side of MscL. Comparison of simulations in which gating progressed to varying degrees identified residues that pose steric hindrance to channel opening.