Conformational flexibility of voltage gated dihydropyridine sensitive calcium channel in hydrated DMPC bilayer.

We have built a model for Ca2+ channel using amino acid sequence from S3 helix of the fourth internal repeat of alpha 1 subunit of dihydropyridine sensitive calcium channel from rabbit skeletal muscle, on the basis of X-ray crystallographic data on four helix bundle. The assembling of the geometry of the pore was achieved using a sixteen residues peptide fragment from short SSI/II loop (residues 1010-1025) which had F1013 and E1014 residues, considered to be important for the drug induced activity of the channel. This had hairpin bend between F1013 to W1016. The drug 2,6-dimethyl 3,5-dicarbomethoxy-4 (2-nitrobenzyl) 1,4 dihydropyridine (DHP) (nifedipine), which is a calcium channel inhibitor used in the treatment of cardiovascular diseases, was introduced, interacting with these two residues via Ca2+ ion. Two more Ca2+ ions were introduced in the pore. The model was incorporated in the bilayer of 36 dimyristoyl phosphatidyl choline (DMPC) molecules with 1201 water molecules and simulated for 200 picoseconds (ps) after equilibration for 120 ps. We also simulated the channel model in vacuum and in aqueous environment for comparison. The latter was unstable after 120 ps. The geometric parameters of the pore are analysed by MOLMOL, PCURVE 3.1 and a special program ANHELIX developed by us. Stability of the pore dimensions during simulations is discussed in this paper.     

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.     

Molecular dynamics investigation of membrane-bound bundles of the channel-forming transmembrane domain of viral protein U from the human immunodeficiency virus HIV-1

Molecular dynamics (MD) simulations have been carried out on bundles of the channel-forming transmembrane (TM) domain of the viral protein U (VPU(1-27) and VPU(6-27)) from the human immunodeficiency virus (HIV-1). Simulations of hexameric and pentameric bundles of VPU(6-27) in an octane/water membrane mimetic system suggested that the pentamer is the preferred oligomer. Accordingly, an unconstrained pentameric helix bundle of VPU(1-27) was then placed in a hydrated palmitoyl-oleyl-3-n-glycero-phosphatidylethanolamine (POPE) lipid bilayer and its structural properties calculated from a 3-ns MD run. Some water molecules, initially inside the channel lumen, were expelled halfway through the simulation and the bundle adopted a conical structure reminiscent of previous MD results obtained for VPU(6-27) in an octane/water system. The pore constriction generated may correspond to a closed state of the channel and underlies the relocation of the W residue toward the pore lumen. The relative positions of the helices with respect to the bilayer and their interactions with the lipids are discussed. The observed structure is stabilized via specific interactions between the VPU helices and the carbonyl oxygen atoms of the lipid molecules, particularly at the Q and S residues.     

Molecular dynamics of synthetic leucine-serine ion channels in a phospholipid membrane

Molecular dynamics calculations were carried out on models of two synthetic leucine-serine ion channels: a tetrameric bundle with sequence (LSLLLSL)(3)NH(2) and a hexameric bundle with sequence (LSSLLSL)(3)NH(2). Each protein bundle is inserted in a palmitoyloleoylphosphatidylcholine bilayer membrane and solvated by simple point charge water molecules inside the pore and at both mouths. Both systems appear to be stable in the absence of an electric field during the 4 ns of molecular dynamics simulation. The water motion in the narrow pore of the four-helix bundle is highly restricted and may provide suitable conditions for proton transfer via a water wire mechanism. In the wider hexameric pore, the water diffuses much more slowly than in bulk but is still mobile. This, along with the dimensions of the pore, supports the observation that this peptide is selective for monovalent cations. Reasonable agreement of predicted conductances with experimentally determined values lends support to the validity of the simulations.     

An alamethicin channel in a lipid bilayer: molecular dynamics simulations

We present the results of 2-ns molecular dynamics (MD) simulations of a hexameric bundle of Alm helices in a 1-palmitoyl-2-oleoylphosphatidylcholine bilayer. These simulations explore the dynamic properties of a model of a helix bundle channel in a complete phospholipid bilayer in an aqueous environment. We explore the stability and conformational dynamics of the bundle in a phospholipid bilayer. We also investigate the effect on bundle stability of the ionization state of the ring of Glu18 side chains. If all of the Glu18 side chains are ionised, the bundle is unstable; if none of the Glu18 side chains are ionized, the bundle is stable. pKA calculations suggest that either zero or one ionized Glu18 is present at neutral pH, correlating with the stable form of the helix bundle. The structural and dynamic properties of water in this model channel were examined. As in earlier in vacuo simulations (Breed et al., 1996 .Biophys. J. 70:1643-1661), the dipole moments of water molecules within the pore were aligned antiparallel to the helix dipoles. This contributes to the stability of the helix bundle.     

Free energy of a potassium ion in a model of the channel formed by an amphipathic leucine-serine peptide

We use molecular dynamics simulations to investigate the position-dependent free energy of a potassium ion in a model of an ion channel formed by the synthetic amphipathic leucine-serine peptide, LS3. The channel model is a parallel bundle of six LS3 helices around which are packed 146 methane-like spheres in order to mimic a membrane. At either end of and within the channel are 1051 water molecules, plus four ions (two potassium and two chloride). The free energy of a potassium ion in the channel was estimated using the weighted histogram analysis (WHAM) method. This is the first time to our knowledge that such a calculation has been carried out as a function of the position of an ion in three dimensions within a channel. The results indicate that for this channel, which is lined by hydrophilic serine sidechains, there is a relatively weak dependence of the free energy on the axial/off-axial position of the ion. There are some off-axis local minima, especially in the C-terminal half of the channel. Using the free energy results, a single channel current-voltage curve was estimated using a one-dimensional Nernst-Planck equation. Although reasonable agreement with experiment is achieved for K(+) ions flowing from the N-terminal to the C-terminal mouth, in the opposite direction the current is underestimated. This underestimation may be a consequence of under-sampling of the conformational dynamics of the channel. We suggest that our simulations may have captured, for example, a sub-conductance level (i.e. an incompletely open state) of the LS3 channel.     

Molecular Dynamics Simulations of Homo-oligomeric Bundles Embedded Within a Lipid Bilayer

Using molecular dynamics simulations, we studied the structure, interhelix interactions, and dynamics of transmembrane proteins. Specifically, we investigated homooligomeric helical bundle systems consisting of synthetic α-helices with either the sequence Ac-(LSLLLSL)₃-NH₂ (LS2) or Ac-(LSSLLSL)₃-NH₂ (LS3). The LS2 and LS3 helical peptides are designed to have amphipathic characteristics that form ion channels in membrane. We simulated bundles containing one to six peptides that were embedded in palmitoyl-oleoyl-phosphatidylcholine (POPC) lipid bilayer and placed between two lamellae of water. We aim to provide a fundamental understanding of how amphipathic helical peptides interact with each other and their dynamical behaviors in different homooligomeric states. To understand structural properties, we examined the helix lengths, tilt angles of individual helices and the entire bundle, interhelix distances, interhelix cross-angles, helix hydrophobic-to-hydrophilic vector projections, and the average number of interhelix hydrophilic (serine–serine) contacts lining the pore of the transmembrane channel. To analyze dynamical properties, we calculated the rotational autocorrelation function of each helix and the cross-correlation of the rotational velocity between adjacent helices. The observed structural and dynamical characteristics show that higher order bundles containing four to six peptides are composed of multiple lower order bundles of one to three peptides. For example, the LS2 channel was found to be stable in a tetrameric bundle composed of a “dimer of dimers.” In addition, we observed that there is a minimum of two strong hydrophilic contacts between a pair of adjacent helices in the dimer to tetramer systems and only one strong hydrophilic interhelix contact in helix pairs of the pentamer and hexamer systems. We believe these results are general and can be applied to more complex ion channels, providing insight into ion channel stability and assembly.     

Molecular Dynamics Simulations of Homo-oligomeric Bundles Embedded Within a Lipid Bilayer

Using molecular dynamics simulations, we studied the structure, interhelix interactions, and dynamics of transmembrane proteins. Specifically, we investigated homooligomeric helical bundle systems consisting of synthetic α-helices with either the sequence Ac-(LSLLLSL)₃-NH₂ (LS2) or Ac-(LSSLLSL)₃-NH₂ (LS3). The LS2 and LS3 helical peptides are designed to have amphipathic characteristics that form ion channels in membrane. We simulated bundles containing one to six peptides that were embedded in palmitoyl-oleoyl-phosphatidylcholine (POPC) lipid bilayer and placed between two lamellae of water. We aim to provide a fundamental understanding of how amphipathic helical peptides interact with each other and their dynamical behaviors in different homooligomeric states. To understand structural properties, we examined the helix lengths, tilt angles of individual helices and the entire bundle, interhelix distances, interhelix cross-angles, helix hydrophobic-to-hydrophilic vector projections, and the average number of interhelix hydrophilic (serine–serine) contacts lining the pore of the transmembrane channel. To analyze dynamical properties, we calculated the rotational autocorrelation function of each helix and the cross-correlation of the rotational velocity between adjacent helices. The observed structural and dynamical characteristics show that higher order bundles containing four to six peptides are composed of multiple lower order bundles of one to three peptides. For example, the LS2 channel was found to be stable in a tetrameric bundle composed of a “dimer of dimers.” In addition, we observed that there is a minimum of two strong hydrophilic contacts between a pair of adjacent helices in the dimer to tetramer systems and only one strong hydrophilic interhelix contact in helix pairs of the pentamer and hexamer systems. We believe these results are general and can be applied to more complex ion channels, providing insight into ion channel stability and assembly.     

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.     

A de novo peptide hexamer with a mutable channel

The design of new proteins that expand the repertoire of natural protein structures represents a formidable challenge. Success in this area would increase understanding of protein structure and present new scaffolds that could be exploited in biotechnology and synthetic biology. Here we describe the design, characterization and X-ray crystal structure of a new coiled-coil protein. The de novo sequence forms a stand-alone, parallel, six-helix bundle with a channel running through it. Although lined exclusively by hydrophobic leucine and isoleucine side chains, the 6-? channel is permeable to water. One layer of leucine residues within the channel is mutable, accepting polar aspartic acid and histidine side chains, which leads to subdivision and organization of solvent within the lumen. Moreover, these mutants can be combined to form a stable and unique (Asp-His)(3) heterohexamer. These new structures provide a basis for engineering de novo proteins with new functions.     

The structure of the M2 channel-lining segment from the nicotinic acetylcholine receptor

The structures of functional peptides corresponding to the predicted channel-lining M2 segment of the nicotinic acetylcholine (AChR) were determined using solution NMR experiments on micelle samples, and solid-state NMR experiments on bilayer samples. The AChR M2 peptide forms a straight transmembrane alpha-helix, with no kinks. M2 inserts in the lipid bilayer at an angle of 12 degrees relative to the bilayer normal, with a rotation about the helix long axis such that the polar residues face the N-terminus of the peptide, which is assigned to be intracellular. A molecular model of the AChR channel pore, constructed from the solid-state NMR 3-D structure of the AChR M2 helix in the membrane assuming a pentameric organization, results in a funnel-like architecture for the channel with the wide opening on the N-terminal intracellular side. A central narrow pore has a diameter ranging from about 3.0 A at its narrowest, to 8.6 A at its widest. Nonpolar residues are predominantly on the exterior of the bundle, while polar residues line the pore. This arrangement is in fair agreement with evidence collected from permeation, mutagenesis, affinity labeling and cysteine accessibility measurements. A pentameric M2 helical bundle may, therefore, represent the structural blueprint for the inner bundle that lines the channel of the nicotinic AChR.     

Intracellular coiled-coil domain engaged in subunit interaction and assembly of melastatin-related transient receptor potential channel 2

TRPM2 channels, activated by adenosine diphosphoribose and related molecules, are assembled as oligomers and most likely tetramers. However, the molecular determinants driving the subunit interaction and assembly of the TRPM2 channels are not well defined. Here we examined, using site-directed mutagenesis in conjunction with co-immunoprecipitation and patch clamp recording, the role of a coiled-coil domain in the intracellular C terminus of TRPM2 subunit in subunit interaction and channel assembly. Deletion of the coiled-coil domain resulted in severe disruption of the subunit interaction and substantial loss of the adenosine diphosphoribose-evoked channel currents. Individual or combined mutations to glutamine of the hydrophobic residues at positions a and d of the abcdef heptad repeat, key residues for protein-protein interaction, significantly reduced the subunit interaction and channel currents; the mutational effects on the subunit interaction and channel currents were clearly correlated. Furthermore, deletion of the coiled-coil domain in a pore mutant subunit abolished its dominant negative phenotypic functional suppression. These results provide strong evidence that the coiled-coil domain is critically engaged in the TRPM2 subunit interaction and such interaction is required for assembly of functional TRPM2 channel. The coiled-coil domain, which is highly conserved within the TRPM subfamily, may serve as a general structural element governing the assembly of TRPM channels.     

Molecular resolution visualization of a pore formed by trichogin,an antimicrobial peptide,in a phospholipid matrix

Electrochemical scanning tunneling microscopy (EC-STM) was employed to study the aggregation of trichogin OMe (TCG), an antimicrobial peptide, incorporated into a lipid monolayer. High-resolution EC-STM images show that trichogin molecules aggregate to form channels in the lipid monolayer. Two types of aggregates were observed in the images. The first consisted of a bundle of six TCG molecules surrounding a central pore. The structure and dimensions of this channel are similar to aggregates that in bilayers are described by the barrel-stave model. The EC-STM images also reveal that channels aggregate further to form a hexagonal lattice of a two dimensional (2D) nanocrystal. The model of 2D lattice was built from trimers of TCG molecules that alternatingly are oriented with either hydrophilic or hydrophobic faces to each other. In this way each TCG molecule is oriented partially with its hydrophilic face towards the hexameric pore allowing the formation of the column of water inside this pore.     

Dynamic properties of individual water molecules in a hydrophobic pore lined with acyl chains: a molecular dynamics study

Recently, a certain class of synthetic molecules has been shown to form ion channels, the pore of which is lined with hydrophobic acyl chains [M. Sokabe, in: F. Oosawa, H. Hayashi, T. Yoshioka (Eds.), Transmembrane Signaling and Sensation, JSSP/VNU Science Press BV, Tokyo, 1984, p. 119; F. Hayashi, M. Sokabe, M. Takagi, K. Hayashi, U. Kishimoto, Biochim. Biophys. Acta, 510 (1978) 305; M.J. Pregel, L. Jullien, J. Canceill, L. Lacombe, J.M. Lehn, J. Chem. Soc. Perkin Trans., 2 (1995) 417; Y. Tanaka, Y. Kobuke, M. Sokabe, Angew. Chem. Int. Ed. Engl., 34 (1995) 693; M. Sokabe, Z. Qi, K. Donowaki, H. Ishida, K. Okubo, Biophys. J., 70 (1996) A201; H. Ishida, K. Donowaki, Y. Inoue, Z. Qi, M. Sokabe, Chem. Lett. (1997) p. 953]. As an initial step towards understanding the physical mechanisms of ion permeation across such a hydrophobic pore, systematic molecular dynamics simulations were performed to investigate dynamic and energetic properties of water molecules inside the pore using a dimer of alanine-N'-acylated cyclic peptide as a channel model. Dynamic energy profiles for water molecules indicated that the energy barrier at the middle region of the pore is approximately 2-3 kcal/mol higher than that in the cap water region which was defined as a vicinity region of the channel entrance. Energetics analyses demonstrated that the mutual interactions among intrapore water molecules are the major factor to give favorable interaction (negative energy contribution) for themselves. The pore, despite being lined with acyl chains, has a favorable van der Waals interaction with intrapore water molecules. These results may help to explain why water-filled channels can be formed by the hydrophobic helices in natural channels.     

The structure of the M2 channel-lining segment from the nicotinic acetylcholine receptor

The structures of functional peptides corresponding to the predicted channel-lining M2 segment of the nicotinic acetylcholine (AChR) were determined using solution NMR experiments on micelle samples, and solid-state NMR experiments on bilayer samples. The AChR M2 peptide forms a straight transmembrane α-helix, with no kinks. M2 inserts in the lipid bilayer at an angle of 12° relative to the bilayer normal, with a rotation about the helix long axis such that the polar residues face the N-terminus of the peptide, which is assigned to be intracellular. A molecular model of the AChR channel pore, constructed from the solid-state NMR 3-D structure of the AChR M2 helix in the membrane assuming a pentameric organization, results in a funnel-like architecture for the channel with the wide opening on the N-terminal intracellular side. A central narrow pore has a diameter ranging from about 3.0 Å at its narrowest, to 8.6 Å at its widest. Nonpolar residues are predominantly on the exterior of the bundle, while polar residues line the pore. This arrangement is in fair agreement with evidence collected from permeation, mutagenesis, affinity labeling and cysteine accessibility measurements. A pentameric M2 helical bundle may, therefore, represent the structural blueprint for the inner bundle that lines the channel of the nicotinic AChR.     

Homotetrameric structure of the SNAP-23 N-terminal coiled-coil domain

SNARE proteins mediate intracellular membrane fusion by forming a coiled-coil complex to merge opposing membranes. A "fusion-active" neuronal SNARE complex is a parallel four-helix bundle containing two coiled-coil domains from SNAP-25 and one coiled-coil domain each from syntaxin-1a and VAMP-2. "Prefusion" assembly intermediate complexes can also form from these SNAREs. We studied the N-terminal coiled-coil domain of SNAP-23 (SNAP-23N), a non-neuronal homologue of SNAP-25, and its interaction with other coiled-coil domains. SNAP-23N can assemble spontaneously with the coiled-coil domains from SNAP-23C, syntaxin-4, and VAMP-3 to form a heterotetrameric complex. Unexpectedly, pure SNAP-23N crystallizes as a coiled-coil homotetrameric complex. The four helices have a parallel orientation and are symmetrical about the long axis. The complex is stabilized through the interaction of conserved hydrophobic residues comprising the a and d positions of the coiled-coil heptad repeats. In addition, a central, highly conserved glutamine residue (Gln-48) is buried within the interface by hydrogen bonding between glutamine side chains derived from adjacent subunits and to solvent molecules. A comparison of the SNAP-23N structure to other SNARE complex structures reveals how a simple coiled-coil motif can form diverse SNARE complexes.     

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).     

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.     

Two possible conducting states of the influenza A virus M2 ion channel

Molecular dynamics simulations have been performed on protonated four-helix bundles based on the 25-residue Duff-Ashley transmembrane sequence of the M2 channel of the influenza A virus. Well-equilibrated tetrameric channels, with one, two and four of the H37 residues protonated, were investigated. The protonated peptide bundles were immersed in the octane portion of a phase-separated water/octane system, which provided a membrane-mimetic environment. The simulations suggest that there could be two conducting states of the M2 channel corresponding to tetramers containing one or two protonated histidines. The more open structure of the doubly protonated state suggests it would have the higher conductance.     

Simulation of molecular dynamics of water movement in ion channels

The motion of water molecules in a gramicidin-like channel was studied by the molecular dynamics method. Water molecules are presented in the ST2 model. The structure of the channel was presented in the form of channel's helix frame possessing mobile dipole groups. The interaction of all mobile particles with the membrane channel's walls was taken into account. The calculation consisted of 50,000 integration steps of delta t = 5 x 10(-16) s which corresponded to a total elapsed time of 25 ps. It was shown that water molecules in the channel did not possess rigid spatial structure but exhibited a structure oriented along the channel axis. The motion of water molecules in the channel occurred smoothly, i.e. all water molecules did not have any deep, stable potential wells in the channel. The distribution of water molecules along the radial coordinate of the channel was estimated. Water density was shown to be maximal near the channel axis.