The structure of the HIV-1 Vpu ion channel: modelling and simulation studies.

Vpu is an 81 amino acid auxiliary protein in HIV-1 which exhibits channel activity. We used two homo-pentameric bundles with the helical transmembrane segments derived from FTIR spectroscopy in combination with a global molecular dynamics search protocol: (i) tryptophans (W) pointing into the pore, and (ii) W facing the lipids. Two equivalent bundles have been generated using a simulated annealing via a restrained molecular dynamics simulations (SA/MD) protocol. A fifth model was generated via SA/MD with all serines facing the pore. The latter model adopts a very stable structure during the 2 ns of simulation. The stability of the models with W facing the pore depends on the starting structure. A possible gating mechanism is outlined.

Transmembrane peptide NB of influenza B: a simulation, structure, and conductance study.

The putative transmembrane segment of the ion channel forming peptide NB from influenza B was synthesized by standard solid-phase peptide synthesis. Insertion into the planar lipid bilayer revealed ion channel activity with conductance levels of 20, 61, 107, and 142 pS in a 0.5 M KCl buffer solution. In addition, levels at -100 mV show conductances of 251 and 413 pS. A linear current-voltage relation reveals a voltage-independent channel formation. In methanol and in vesicles the peptide appears to adopt an alpha-helical-like structure. Computational models of alpha-helix bundles using N = 4, 5, and 6 NB peptides per bundle revealed water-filled pores after 1 ns of MD simulation in a solvated lipid bilayer. Calculated conductance values [using HOLE (Smart et al. (1997) Biophys. J. 72, 1109-1126)] of ca. 20, 60, and 90 pS, respectively, suggested that the multiple conductance levels seen experimentally must correspond to different degrees of oligomerization of the peptide to form channels.

Molecular contacts in the transmembrane c-subunit oligomer of F-ATPases identified by tryptophan substitution mutagenesis.

When isolated in its monomeric form, subunit c of the proton transporting ATP synthase of Escherichia coli was shown to fold in a hairpin-like structure consisting of two hydrophobic membrane spanning helices and a short connecting hydrophilic loop. In the plasma membrane of Escherichia coli, however, about 9-12 c-subunit monomers form an oligomeric complex that functions in transmembrane proton conduction and in energy transduction to the catalytic F1 domain. The arrangement of the monomers and the molecular architecture of the complex were studied by tryptophan scanning mutagenesis and restrained MD simulations. Residues 12-24 of the N-terminal transmembrane segment of subunit c were individually substituted by the large and moderately hydrophobic tryptophan side chain. Effects on the activity of the mutant proteins were studied in selective growth experiments and various ATP synthase specific activity assays. The results identify potential intersubunit contacts and structurally non-distorted, accessible residues in the c-oligomer and add constraints to the arrangement of monomers in the oligomeric complex. Results from our mutagenesis experiments were interpreted in structural models of the c-oligomer that have been obtained by restrained MD simulations. Different stoichiometries and monomer orientations were applied in these calculations. A cylindrical complex consisting of 10 monomers that are arranged in two concentric rings with the N-terminal helices of the monomers located at the periphery shows the best match with the experimental data.

Homology modeling and molecular dynamics simulation studies of an inward rectifier potassium channel.

A homology model has been generated for the pore-forming domain of Kir6.2, a component of an ATP-sensitive K channel, based on the x-ray structure of the bacterial channel KcsA. Analysis of the lipid-exposed and pore-lining surfaces of the model reveals them to be compatible with the known features of membrane proteins and Kir channels, respectively. The Kir6.2 homology model was used as the starting point for nanosecond-duration molecular dynamics simulations in a solvated phospholipid bilayer. The overall drift from the model structure was comparable to that seen for KcsA in previous similar simulations. Preliminary analysis of the interactions of the Kir6.2 channel model with K(+) ions and water molecules during these simulations suggests that concerted single-file motion of K(+) ions and water through the selectivity filter occurs. This is similar to such motion observed in simulations of KcsA. This suggests that a single-filing mechanism is conserved between different K channel structures and may be robust to changes in simulation details. Comparison of Kir6.2 and KcsA suggests some degree of flexibility in the filter, thus complicating models of ion selectivity based upon a rigid filter.

Membrane simulations: bigger and better?

Molecular dynamics simulations of biological membranes have come of age. Simulations of pure lipid bilayers are extending our understanding of both optimal simulation procedures and the detailed structural dynamics of lipids in these systems. Simulation methods established using simple bilayer-embedded peptides are being extended to a wide range of membrane proteins and membrane protein models, and are beginning to reveal some of the complexities of membrane protein structural dynamics and their relationship to biological function.

Transmembrane domains of viral ion channel proteins: a molecular dynamics simulation study.

Nanosecond molecular dynamics simulations in a fully solvated phospholipid bilayer have been performed on single transmembrane alpha-helices from three putative ion channel proteins encoded by viruses: NB (from influenza B), CM2 (from influenza C), and Vpu (from HIV-1). alpha-Helix stability is maintained within a core region of ca. 28 residues for each protein. Helix perturbations are due either to unfavorable interactions of hydrophobic residues with the lipid headgroups or to the need of the termini of short helices to extend into the surrounding interfacial environment in order to form H-bonds. The requirement of both ends of a helix to form favorable interactions with lipid headgroups and/or water may also lead to tilting and/or kinking of a transmembrane alpha-helix. Residues that are generally viewed as poor helix formers in aqueous solution (e.g., Gly, Ile, Val) do not destabilize helices, if located within a helix that spans a lipid bilayer. However, helix/bilayer mismatch such that a helix ends abruptly within the bilayer core destabilizes the end of the helix, especially in the presence of Gly and Ala residues. Hydrogen bonding of polar side-chains with the peptide backbone and with one another occurs when such residues are present within the bilayer core, thus minimizing the energetic cost of burying such side-chains.

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.

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

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.

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.

Defining the transmembrane helix of M2 protein from influenza A by molecular dynamics simulations in a lipid bilayer.

Integral membrane proteins containing at least one transmembrane (TM) alpha-helix are believed to account for between 20% and 30% of most genomes. There are several algorithms that accurately predict the number and position of TM helices within a membrane protein sequence. However, these methods tend to disagree over the beginning and end residues of TM helices, posing problems for subsequent modeling and simulation studies. Molecular dynamics (MD) simulations in an explicit lipid and water environment are used to help define the TM helix of the M2 protein from influenza A virus. Based on a comparison of the results of five different secondary structure prediction algorithms, three different helix lengths (an 18mer, a 26mer, and a 34mer) were simulated. Each simulation system contained 127 POPC molecules plus approximately 3500-4700 waters, giving a total of approximately 18,000-21,000 atoms. Two simulations, each of 2 ns duration, were run for the 18mer and 26mer, and five separate simulations were run for the 34mer, using different starting models generated by restrained in vacuo MD simulations. The total simulation time amounted to 11 ns. Analysis of the time-dependent secondary structure of the TM segments was used to define the regions that adopted a stable alpha-helical conformation throughout the simulation. This analysis indicates a core TM region of approximately 20 residues (from residue 22 to residue 43) that remained in an alpha-helical conformation. Analysis of atomic density profiles suggested that the 18mer helix revealed a local perturbation of the lipid bilayer. Polar side chains on either side of this region form relatively long-lived H-bonds to lipid headgroups and water molecules.

Lipid properties and the orientation of aromatic residues in OmpF, influenza M2, and alamethicin systems: molecular dynamics simulations.

Molecular dynamics simulations allow a direct study of the structure and dynamics of membrane proteins and lipids. We describe the behavior of aromatic residues and lipid properties in POPE and POPC bilayer models with the Escherichia coli OmpF trimer, single alamethicin and Influenza M2 helices, 4-helix M2 bundles, and two alamethicin 6-helix channel models. The total simulation time is over 24 ns, of systems containing solvent, protein, and between 104 and 318 lipids. Various types of adjustment between lipids and proteins occur, depending on the size of the protein and the degree of hydrophobic mismatch between lipid and protein. Single helices cause little measurable effect on nearby lipids whereas the 4-helix bundles, 6-helix channel models, and OmpF cause a significant lowering of order parameters in nearby lipid chains, an increased difference between odd and even chain dihedrals in the magnitude of the trans dihedral fractions and dihedral transition rates, and in most cases a decreased gauche population and a decrease in bilayer thickness. An increased tilt of the lipid chains near the proteins can account for most of the observed decrease in order parameters. The orientation of tryptophans and tyrosines on the outside of the proteins is determined by packing at the protein exterior and non-specific hydrogen bonding with lipids and solvent. The tyrosines in the broad bands that delimit the hydrophobic exterior of OmpF show little change in orientation over one nanosecond. Their rings are oriented predominantly perpendicular to the bilayer plane, with the hydroxyl group pointing toward the lipid-water interface. Phenylalanines in OmpF, alamethicin, and Influenza M2 are more mobile and assume a variety of orientations.

Viral ion channels: molecular modeling and simulation.

In a number of membrane-bound viruses, ion channels are formed by integral membrane proteins. These channel proteins include M2 from influenza A, NB from influenza B, and, possibly, Vpu from HIV-1. M2 is important in facilitating uncoating of the influenza A viral genome and is the target of amantadine, an anti-influenza drug. The biological roles of NB and Vpu are less certain. In all cases, the protein contains a single transmembrane alpha-helix close to its N-terminus. Channels can be formed by homo-oligomerization of these proteins, yielding bundles of transmembrane helices that span the membrane and surround a central ion-permeable pore. Molecular modeling may be used to integrate and interpret available experimental data concerning the structure of such transmembrane pores. This has proved successful for the M2 channel domain, where two independently derived models are in agreement with one another, and with solid-state nuclear magnetic resonance (NMR) data. Simulations based on channel models may yield insights into possible ion conduction and selectivity mechanisms.

Two models of the influenza A M2 channel domain: verification by comparison.

BACKGROUND: The influenza M2 protein is a simple membrane protein, containing a single transmembrane helix. It is representative of a very large family of single-transmembrane helix proteins. The functional protein is a tetramer, with the four transmembrane helices forming a proton-permeable channel across the bilayer. Two independently derived models of the M2 channel domain are compared, in order to assess the success of applying molecular modelling approaches to simple membrane proteins. RESULTS: The Calpha RSMD between the two models is 1.7 A. Both models are composed of a left-handed bundle of helices, with the helices tilted roughly 15 degrees relative to the (presumed) bilayer normal. The two models have similar pore radius profiles, with a pore cavity lined by the Ser31 and Gly34 residues and a pore constriction formed by the ring of His37 residues. CONCLUSIONS: Independent studies of M2 have converged on the same structural model for the channel domain. This model is in agreement with solid state NMR data. In particular, both model and NMR data indicate that the M2 helices are tilted relative to the bilayer normal and form a left-handed bundle. Such convergence suggests that, at least for simple membrane proteins, restraints-directed modelling might yield plausible models worthy of further computational and experimental investigation.