Molecular dynamics simulations of water within models of ion channels.

The transbilayer pores formed by ion channel proteins contain extended columns of water molecules. The dynamic properties of such waters have been suggested to differ from those of water in its bulk state. Molecular dynamics simulations of ion channel models solvated within and at the mouths of their pores are used to investigate the dynamics and structure of intra-pore water. Three classes of channel model are investigated: a) parallel bundles of hydrophobic (Ala20) alpha-helices; b) eight-stranded hydrophobic (Ala10) antiparallel beta-barrels; and c) parallel bundles of amphipathic alpha-helices (namely, delta-toxin, alamethicin, and nicotinic acetylcholine receptor M2 helix). The self-diffusion coefficients of water molecules within the pores are reduced significantly relative to bulk water in all of the models. Water rotational reorientation rates are also reduced within the pores, particularly in those pores formed by alpha-helix bundles. In the narrowest pore (that of the Ala20 pentameric helix bundle) self-diffusion coefficients and reorientation rates of intra-pore waters are reduced by approximately an order of magnitude relative to bulk solvent. In Ala20 helix bundles the water dipoles orient antiparallel to the helix dipoles. Such dipole/dipole interaction between water and pore may explain how water-filled ion channels may be formed by hydrophobic helices. In the bundles of amphipathic helices the orientation of water dipoles is modulated by the presence of charged side chains. No preferential orientation of water dipoles relative to the pore axis is observed in the hydrophobic beta-barrel models.     

Electrostatics and the ion selectivity of ligand-gated channels.

The nicotinic acetylcholine receptor (nAChR) is a cation-selective ion channel that opens in response to acetylcholine binding. The related glycine receptor (GlyR) is anion selective. The pore-lining domain of each protein may be modeled as a bundle of five parallel M2 helices. Models of the pore-lining domains of homopentameric nAChR and GlyR have been used in continuum electrostatics calculations to probe the origins of ion selectivity. Calculated pKA values suggest that "rings" of acidic or basic side chains at the mouths of the nAChR or GlyR M2 helix bundles, respectively, may not be fully ionized. In particular, for the nAChR the ring of glutamate side chains at the extracellular mouth of the pore is predicted to be largely protonated at neutral pH, whereas those glutamate side chains in the intracellular and intermediate rings (at the opposite mouth of the pore) are predicted to be fully ionized. Inclusion of the other domains of each protein represented as an irregular cylindrical tube in which the M2 bundles are embedded suggests that both the M2 helices and the extramembrane domains play significant roles in determining ion selectivity.     

The nicotinic acetylcholine receptor: from molecular model to single-channel conductance

The nicotinic acetylcholine receptor (nAChR) is the archetypal ligand-gated ion channel. A model of the α7 homopentameric nAChR is described in which the pore-lining M2 helix bundle is treated atomistically and the remainder of the molecule is treated as a “low resolution” cylinder. The surface charge on the cylinder is derived from the distribution of charged amino acids in the amino acid sequence (excluding the M2 segments). This model is explored in terms of its predicted single-channel properties. Based on electrostatic potential profiles derived from the model, the one-dimensional Poisson-Nernst-Planck equation is used to calculate single-channel current/voltage curves. The predicted single-channel conductance is three times higher (ca. 150 pS) than that measured experimentally, and the predicted ion selectivity agrees with the observed cation selectivity of nAChR. Molecular dynamics (MD) simulations are used to estimate the self-diffusion coefficients (D) of water molecules within the channel. In the narrowest region of the pore, D is reduced ca. threefold relative to that of bulk water. Assuming that the diffusion of ions scales with that of water, this yields a revised prediction of the single-channel conductance (ca. 50 pS) in good agreement with the experimental value. We conclude that combining atomistic (MD) and continuum electrostatics calculations is a promising approach to bridging the gap between structure and physiology of ion channels. Received: 2 August 1999 / Revised version: 5 November 1999 / Accepted: 9 November 1999     

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.     

A hydrophobic gate in an ion channel: the closed state of the nicotinic acetylcholine receptor

The nicotinic acetylcholine receptor (nAChR) is the prototypic member of the 'Cys-loop' superfamily of ligand-gated ion channels which mediate synaptic neurotransmission, and whose other members include receptors for glycine, gamma-aminobutyric acid and serotonin. Cryo-electron microscopy has yielded a three-dimensional structure of the nAChR in its closed state. However, the exact nature and location of the channel gate remains uncertain. Although the transmembrane pore is constricted close to its center, it is not completely occluded. Rather, the pore has a central hydrophobic zone of radius about 3 A. Model calculations suggest that such a constriction may form a hydrophobic gate, preventing movement of ions through a channel. We present a detailed and quantitative simulation study of the hydrophobic gating model of the nicotinic receptor, in order to fully evaluate this hypothesis. We demonstrate that the hydrophobic constriction of the nAChR pore indeed forms a closed gate. Potential of mean force (PMF) calculations reveal that the constriction presents a barrier of height about 10 kT to the permeation of sodium ions, placing an upper bound on the closed channel conductance of 0.3 pS. Thus, a 3 A radius hydrophobic pore can form a functional barrier to the permeation of a 1 A radius Na+ ion. Using a united-atom force field for the protein instead of an all-atom one retains the qualitative features but results in differing conductances, showing that the PMF is sensitive to the detailed molecular interactions.     

Molecular dynamics simulations of a bacterial autotransporter: NalP from Neisseria meningitidis

NalP is an autotransporter secretory protein found in the outer membrane of Neisseria meningitidis. The crystal structure of the NalP translocator domain revealed a transmembrane beta-barrel containing a central alpha-helix. The role of this alpha-helix, and of the conformational dynamics of the beta-barrel pore have been studied via atomistic molecular dynamics simulations. Three simulations, each of 10 ns duration, of NalP embedded within a solvated DMPC bilayer were performed. The helix was removed from the barrel interior in one simulation. The conformational stability of the protein is similar to that of other outer membrane proteins, e.g., OmpA, in comparable simulations. The transmembrane beta-barrel is stable even in the absence of the alpha-helix. Removal of the helix results in an influx of water into the pore region, suggesting the helix acts as a 'plug'. Water molecules entering the resultant pore form hydrogen bonds with the barrel lining that compensate for the loss of helix-barrel hydrogen bonds. The dimensions of the pore fluctuate over the course of the simulation revealing it to be flexible, but only wide enough to allow transport of the passenger domain in an unfolded or extended conformation. The simulations help us to understand the role of the central helix in plugging the pore and in maintaining the width of the barrel, and show that the NalP monomer is sufficient for the transport of the passenger domain in an unfolded or extended conformation.     

Probing the structure of the affinity-purified and lipid-reconstituted torpedo nicotinic acetylcholine receptor

The Torpedo nicotinic acetylcholine receptor (nAChR) is the only member of the Cys-loop superfamily of ligand-gated ion channels (LGICs) that is available in high abundance in a native membrane preparation. To study the structure of the other LGICs using biochemical and biophysical techniques, detergent solubilization, purification, and lipid reconstitution are usually required. To assess the effects of purification on receptor structure, we used the hydrophobic photoreactive probe 3-trifluoromethyl-3-(m-[(125)I]iodophenyl)diazirine ([(125)I]TID) to compare the state-dependent photolabeling of the Torpedo nAChR before and after purification and reincorporation into lipid. For the purified nAChR, the agonist-sensitive photolabeling within the M2 ion channel domain of positions M2-6, M2-9, and M2-13, the agonist-enhanced labeling of deltaThr274 (deltaM2-18) within the delta subunit helix bundle, and the labeling at the lipid-protein interface (alphaMu4) were the same as for the nAChR in native membranes. However, addition of agonist did not enhance [(125)I]TID photolabeling of deltaIle288 within the deltaM2-M3 loop. These results indicate that after purification and reconstitution of the Torpedo nAChR, the difference in structure between the resting and desensitized states within the M2 ion channel domain was preserved, but not the agonist-dependent change of structure of the deltaM2-M3 loop. To further characterize the pharmacology of [(125)I]TID binding sites in the nAChR in the desensitized state, we examined the effect of phencyclidine (PCP) on [(125)I]TID photolabeling. PCP inhibited [(125)I]TID labeling of amino acids at the cytoplasmic end of the ion channel (M2-2 and M2-6) while potentiating labeling at M2-9 and M2-13 and allosterically modulating the labeling of amino acids within the delta subunit helix bundle.     

Molecular dynamics simulations of a bacterial autotransporter: NalP from Neisseria meningitidis

NalP is an autotransporter secretory protein found in the outer membrane of Neisseria meningitidis. The crystal structure of the NalP translocator domain revealed a transmembrane β-barrel containing a central α-helix. The role of this α-helix, and of the conformational dynamics of the β-barrel pore have been studied via atomistic molecular dynamics simulations. Three simulations, each of 10 ns duration, of NalP embedded within a solvated DMPC bilayer were performed. The helix was removed from the barrel interior in one simulation. The conformational stability of the protein is similar to that of other outer membrane proteins, e.g., OmpA, in comparable simulations. The transmembrane β-barrel is stable even in the absence of the α-helix. Removal of the helix results in an influx of water into the pore region, suggesting the helix acts as a ‘plug’. Water molecules entering the resultant pore form hydrogen bonds with the barrel lining that compensate for the loss of helix-barrel hydrogen bonds. The dimensions of the pore fluctuate over the course of the simulation revealing it to be flexible, but only wide enough to allow transport of the passenger domain in an unfolded or extended conformation. The simulations help us to understand the role of the central helix in plugging the pore and in maintaining the width of the barrel, and show that the NalP monomer is sufficient for the transport of the passenger domain in an unfolded or extended conformation.     

Parallel helix bundles and ion channels: molecular modeling via simulated annealing and restrained molecular dynamics.

A parallel bundle of transmembrane (TM) alpha-helices surrounding a central pore is present in several classes of ion channel, including the nicotinic acetylcholine receptor (nAChR). We have modeled bundles of hydrophobic and of amphipathic helices using simulated annealing via restrained molecular dynamics. Bundles of Ala20 helices, with N = 4, 5, or 6 helices/bundle were generated. For all three N values the helices formed left-handed coiled coils, with pitches ranging from 160 A (N = 4) to 240 A (N = 6). Pore radius profiles revealed constrictions at residues 3, 6, 10, 13, and 17. A left-handed coiled coil and a similar pattern of pore constrictions were observed for N = 5 bundles of Leu20. In contrast, N = 5 bundles of Ile20 formed right-handed coiled coils, reflecting loosened packing of helices containing beta-branched side chains. Bundles formed by each of two classes of amphipathic helices were examined: (a) M2a, M2b, and M2c derived from sequences of M2 helices of nAChR; and (b) (LSSLLSL)3, a synthetic channel-forming peptide. Both classes of amphipathic helix formed left-handed coiled coils. For (LSSLLSL)3 the pitch of the coil increased as N increased from 4 to 6. The M2c N = 5 helix bundle is discussed in the context of possible models of the pore domain of nAChR.     

Pores formed by the nicotinic receptor m2delta Peptide: a molecular dynamics simulation study

The M2delta peptide self-assembles to form a pentameric bundle of transmembrane alpha-helices that is a model of the pore-lining region of the nicotinic acetylcholine receptor. Long (>15 ns) molecular dynamics simulations of a model of the M2delta(5) bundle in a POPC bilayer have been used to explore the conformational dynamics of the channel assembly. On the timescale of the simulation, the bundle remains relatively stable, with the polar pore-lining side chains remaining exposed to the lumen of the channel. Fluctuations at the helix termini, and in the helix curvature, result in closing/opening transitions at both mouths of the channel, on a timescale of approximately 10 ns. On average, water within the pore lumen diffuses approximately 4x more slowly than water outside the channel. Examination of pore water trajectories reveals both single-file and path-crossing regimes to occur at different times within the simulation.     

The dielectric properties of water within model transbilayer pores.

Ion channels contain extended columns of water molecules within their transbilayer pores. The dynamic properties of such intrapore water have been shown to differ from those of water in its bulk state. In previous molecular dynamics simulations of two classes of model pore (parallel bundles of Ala20 alpha-helices and antiparallel barrels of Ala10 beta-strands), a substantially reduced translational and rotational mobility of waters was observed within the pore relative to bulk water. Molecular dynamics simulations in the presence of a transpore electrostatic field (i.e., a voltage drop along the pore axis) have been used to estimate the resultant polarization (due to reorientation) of the intrapore water, and hence to determine the local dielectric behavior within the pore. It is shown that the local dielectric constant of water within a pore is reduced for models formed by parallel alpha-helix bundles, but not by those formed by beta-barrels. This result is discussed in the context of electrostatics calculations of ion permeation through channels, and the effect of the local dielectric of water within a helix bundle pore is illustrated with a simple Poisson-Boltzmann calculation.     

A novel mechanism of pore formation: membrane penetration by the N-terminal amphipathic region of equinatoxin

Equinatoxin II is a representative of actinoporins, eukaryotic pore-forming toxins from sea anemones. It creates pores in natural and artificial lipid membranes by an association of three or four monomers. Cysteine-scanning mutagenesis was used to study the structure of the N terminus, which is proposed to be crucial in transmembrane pore formation. We provide data for two steps of pore formation: a lipid-bound monomeric intermediate state and a final oligomeric pore. Results show that residues 10-28 are organized as an alpha-helix in both steps. In the first step, the whole region is transferred to a lipid-water interface, laying flat on the membrane. In the pore-forming state, the hydrophilic side of the amphipathic helix lines the pore lumen. The pore has a restriction around Asp-10, according to the permeabilization ratio of ions flowing through pores formed by chemically modified mutants. A general model was introduced to derive the tilt angle of the helix from the ion current data. This study reveals that actinoporins use a unique single helix insertion mechanism for pore formation.     

Modeling noncompetitive antagonism of a nicotinic acetylcholine receptor

Models of closed and open channel pores of a muscle-type nicotinic acetylcholine receptor (nAChR) channel comprising M1 and M2 segments are presented. A model of the closed channel is proposed in which hydrophobic residues of the Equatorial Leucine ring screen the oxygen domain formed by the Serine ring, thereby preventing ion flux without completely occluding the pore. This model demonstrates a high similarity with the structure derived from a recent electron microscopy study. We propose that hydrophobic residues of the Equatorial Leucine ring are retracted when the pore is open. Our models provide a possible resolution of the nAChR gate controversy. We have also obtained explanations for the complex mechanisms underlying inhibition of nAChR by philanthotoxins (PhTXs). PhTX-343, containing a spermine moiety with a charge of +3, binds deep in the pore near the Serine ring where classical open channel blockers of nAChR bind. In contrast, PhTX-(12), which has a single charged amino group is unable to reach deeply located rings because of steric restrictions. Both philanthotoxins may bind to a hydrophobic site located close to the external entrance of the pore in a region that includes residues associated with the regulation of desensitization.     

Genetic disorders caused by mutated acetylcholine receptors

The nicotinic acetylcholine receptors (nAChRs) are members of the large family of ligand-gated ion channels and are constituted by the assembly of five subunits arranged pseudosymmetrically around the central axis that forms a cation-selective ion pore. They are widely distributed in both the nervous system and non-neuronal tissues, and can be activated by endogenous agonists such as acetylcholine or exogenous ligands such as nicotine. Mutations in neuronal nAChRs are found in a rare form of familial nocturnal frontal lobe epilepsy (ADNFLE), while mutations in the neuromuscular subtype of the nAChR are responsible for either congenital myasthenia syndromes (adult subtype of neuromuscular nAChR) or a form of arthrogryposis multiplex congenita type Escobar (fetal subtype of neuromuscular nAChR).     

Mechanics of channel gating of the nicotinic acetylcholine receptor

The nicotinic acetylcholine receptor (nAChR) is a key molecule involved in the propagation of signals in the central nervous system and peripheral synapses. Although numerous computational and experimental studies have been performed on this receptor, the structural dynamics of the receptor underlying the gating mechanism is still unclear. To address the mechanical fundamentals of nAChR gating, both conventional molecular dynamics (CMD) and steered rotation molecular dynamics (SRMD) simulations have been conducted on the cryo-electron microscopy (cryo-EM) structure of nAChR embedded in a dipalmitoylphosphatidylcholine (DPPC) bilayer and water molecules. A 30-ns CMD simulation revealed a collective motion amongst C-loops, M1, and M2 helices. The inward movement of C-loops accompanying the shrinking of acetylcholine (ACh) binding pockets induced an inward and upward motion of the outer β-sheet composed of β9 and β10 strands, which in turn causes M1 and M2 to undergo anticlockwise motions around the pore axis. Rotational motion of the entire receptor around the pore axis and twisting motions among extracellular (EC), transmembrane (TM), and intracellular MA domains were also detected by the CMD simulation. Moreover, M2 helices undergo a local twisting motion synthesized by their bending vibration and rotation. The hinge of either twisting motion or bending vibration is located at the middle of M2, possibly the gate of the receptor. A complementary twisting-to-open motion throughout the receptor was detected by a normal mode analysis (NMA). To mimic the pulsive action of ACh binding, nonequilibrium MD simulations were performed by using the SRMD method developed in one of our laboratories. The result confirmed all the motions derived from the CMD simulation and NMA. In addition, the SRMD simulation indicated that the channel may undergo an open-close (O ↔ C) motion. The present MD simulations explore the structural dynamics of the receptor under its gating process and provide a new insight into the gating mechanism of nAChR at the atomic level.     

Molecular dynamics studies of the archaeal translocon

The translocon is a protein-conducting channel conserved over all domains of life that serves to translocate proteins across or into membranes. Although this channel has been well studied for many years, the recent discovery of a high-resolution crystal structure opens up new avenues of exploration. Taking advantage of this, we performed molecular dynamics simulations of the translocon in a fully solvated lipid bilayer, examining the translocation abilities of monomeric SecYEbeta by forcing two helices comprised of different amino acid sequences to cross the channel. The simulations revealed that the so-called plug of SecYEbeta swings open during translocation, closing thereafter. Likewise, it was established that the so-called pore ring region of SecYEbeta forms an elastic, yet tight, seal around the translocating oligopeptides. The closed state of the channel was found to block permeation of all ions and water molecules; in the open state, ions were blocked. Our results suggest that the SecYEbeta monomer is capable of forming an active channel.     

Multiple transmembrane binding sites for p-trifluoromethyldiazirinyl-etomidate, a photoreactive Torpedo nicotinic acetylcholine receptor allosteric inhibitor

Photoreactive derivatives of the general anesthetic etomidate have been developed to identify their binding sites in γ-aminobutyric acid, type A and nicotinic acetylcholine receptors. One such drug, [(3)H]TDBzl-etomidate (4-[3-(trifluoromethyl)-3H-diazirin-3-yl]benzyl-[(3)H]1-(1-phenylethyl)-1H-imidazole-5-carboxylate), acts as a positive allosteric potentiator of Torpedo nACh receptor (nAChR) and binds to a novel site in the transmembrane domain at the γ-α subunit interface. To extend our understanding of the locations of allosteric modulator binding sites in the nAChR, we now characterize the interactions of a second aryl diazirine etomidate derivative, TFD-etomidate (ethyl-1-(1-(4-(3-trifluoromethyl)-3H-diazirin-3-yl)phenylethyl)-1H-imidazole-5-carboxylate). TFD-etomidate inhibited acetylcholine-induced currents with an IC(50) = 4 μM, whereas it inhibited the binding of [(3)H]phencyclidine to the Torpedo nAChR ion channel in the resting and desensitized states with IC(50) values of 2.5 and 0.7 mm, respectively. Similar to [(3)H]TDBzl-etomidate, [(3)H]TFD-etomidate bound to a site at the γ-α subunit interface, photolabeling αM2-10 (αSer-252) and γMet-295 and γMet-299 within γM3, and to a site in the ion channel, photolabeling amino acids within each subunit M2 helix that line the lumen of the ion channel. In addition, [(3)H]TFD-etomidate photolabeled in an agonist-dependent manner amino acids within the δ subunit M2-M3 loop (δIle-288) and the δ subunit transmembrane helix bundle (δPhe-232 and δCys-236 within δM1). The fact that TFD-etomidate does not compete with ion channel blockers at concentrations that inhibit acetylcholine responses indicates that binding to sites at the γ-α subunit interface and/or within δ subunit helix bundle mediates the TFD-etomidate inhibitory effect. These results also suggest that the γ-α subunit interface is a binding site for Torpedo nAChR negative allosteric modulators (TFD-etomidate) and for positive modulators (TDBzl-etomidate).     

Exploring models of the influenza A M2 channel: MD simulations in a phospholipid bilayer

The M2 protein of influenza A virus forms homotetrameric helix bundles, which function as proton-selective channels. The native form of the protein is 97 residues long, although peptides representing the transmembrane section display ion channel activity, which (like the native channel) is blocked by the antiviral drug amantadine. As a small ion channel, M2 may provide useful insights into more complex channel systems. Models of tetrameric bundles of helices containing either 18 or 22 residues have been simulated while embedded in a fully hydrated 1-palmitoyl-2-oleoyl-sn-glycerol-3-phosphatidylcholine bilayer. Several different starting models have been used. These suggest that the simulation results, at least on a nanosecond time scale, are sensitive to the exact starting structure. Electrostatics calculations carried out on a ring of four ionizable aspartate residues at the N-terminal mouth of the channel suggest that at any one time, only one will be in a charged state. Helix bundle models were mostly stable over the duration of the simulation, and their helices remained tilted relative to the bilayer normal. The M2 helix bundles form closed channels that undergo breathing motions, alternating between a tetramer and a dimer-of-dimers structure. Under these conditions either the channel forms a pocket of trapped waters or it contains a column of waters broken predominantly at the C-terminal mouth of the pore. These waters exhibit restricted motion in the pore and are effectively "frozen" in a way similar to those seen in previous simulations of a proton channel formed by a four-helix bundle of a synthetic leucine-serine peptide (, Biophys. J. 77:2400-2410).     

Molecular dynamics simulations and KcsA channel gating

The gating mechanism of a bacterial potassium channel, KcsA, has been investigated via multi-nanosecond molecular dynamic simulations of the channel molecules embedded in a fully solvated palmitoyloleoylphosphatidylcholine bilayer. Four events are seen in which a cation (K(+) or, in one case, Na(+)) initially present in the central cavity exits through the intracellular mouth (the presumed gate) of the channel. Whilst in the cavity a cation interacts with the sidechain T107 O gamma atom of one of the subunits prior to its exit from the channel. Secondary structure analysis as a function of time reveals a break in the helicity of one of the M2 helices. This break is expected to lend flexibility to the helices, enabling them to "open" (minimum pore radius >0.13 nm) and "close" (minimum pore radius <0.13 nm) the channel. Fluctuations in the pore radius at the intracellular gate region are of the order of 0.05 nm, with an average radius in the region of the gate of ca. 0.1 nm. However, around the time of exit of a cation, the pore widens to about 0.15 nm. The distances between the C alpha atoms of the inner helices M2 reveal a coupled increase and decrease between the opposite pair of helices at about the time of exit of the ion. This suggests a breathing motion of the M2 helices that may form the basis for a gating mechanism.