Search

Colvars is an open-source C++ library that provides a modular toolkit for collective-variable-based molecular simulations. It allows practitioners to easily create and implement descriptors that best fit a process of interest and to apply a wide range of biasing algorithms in collective variable space. This paper reviews several features and improvements to Colvars that were added since its original introduction. Special attention is given to contributions that significantly expanded the capabilities of this software or its distribution with major MD simulation packages. Collective variables can now be optimized either manually or by machine-learning methods, and the space of descriptors can be explored interactively using the graphical interface included in VMD. Beyond the spatial coordinates of individual molecules, Colvars can now apply biasing forces to mesoscale structures and alchemical degrees of freedom and perform simulations guided by experimental data within ensemble averages or probability distributions. It also features advanced computational schemes to boost the accuracy, robustness, and general applicability of simulation methods, including extended-system and multiple-walker adaptive biasing force, boundary conditions for metadynamics, replica exchange with biasing potentials, and adiabatic bias molecular dynamics. The library is made available directly within the main distributions of the academic software GROMACS, LAMMPS, NAMD, Tinker-HP, and VMD. The robustness of the software and the reliability of the results are ensured through the use of continuous integration with a test suite within the source repository.

Membrane transporters mediate the passage of molecules across membranes and are essential for cellular function. While the transmembrane region of these proteins is responsible for substrate transport, often the cytoplasmic regions are required for modulating their activity. However, it can be difficult to obtain atomic-resolution descriptions of these autoregulatory domains by classical structural biology techniques, especially if they lack a single, defined structure. The betaine permease, BetP, a homotrimer, is a prominent and well-studied example of a membrane protein whose autoregulation depends on cytoplasmic N- and C-terminal segments. These domains sense and transduce changes in K+concentration and in lipid bilayer properties caused by osmotic stress. However, structural data for these terminal domains is incomplete, which hinders a clear description of the molecular mechanism of autoregulation. Here we used microsecond-scale molecular simulations of the BetP trimer to compare reported conformations of the 45-amino-acid long C-terminal tails. The simulations provide support for the idea that the conformation derived from electron microscopy (EM) data represents a more stable global orientation of the C-terminal segment under downregulating conditions while also providing a detailed molecular description of its dynamics and highlighting specific interactions with lipids, ions, and neighboring transporter subunits. A missing piece of the molecular puzzle is the N-terminal segment, whose dynamic nature has prevented structural characterization. Using Rosetta to generate ensembles of de novo conformations in the context of the EM-derived structure robustly identifies two features of the N-terminal tail, namely 1) short helical elements and 2) an orientation that would confine potential interactions to the protomer in the counterclockwise direction (viewed from the cytoplasm). Since each C-terminal tail only contacts the protomer in the clockwise direction, these results indicate an intricate interplay between the three protomers of BetP in the downregulated protein and a multidirectionality that may facilitate autoregulation of transport.

The cytoplasmic heme binding protein from Pseudomonas aeruginosa, PhuS, plays two essential roles in regulating heme uptake and iron homeostasis. First, PhuS shuttles exogenous heme to heme oxygenase (HemO) for degradation and iron release. Second, PhuS binds DNA and modulates the transcription of the prrF/H small RNAs (sRNAs) involved in the iron-sparing response. Heme binding to PhuS regulates this dual function, as the unliganded form binds DNA, whereas the heme-bound form binds HemO. Crystallographic studies revealed nearly identical structures for apo- and holo-PhuS, and yet numerous solution-based measurements indicate that heme binding is accompanied by large conformational rearrangements. In particular, hydrogen-deuterium exchange mass spectrometry (HDX-MS) of apo- versus holo-PhuS revealed large differences in deuterium uptake, notably in α-helices 6, 7, and 8 (α6,7,8), which contribute to the heme binding pocket. These helices were mostly labile in apo-PhuS but largely protected in holo-PhuS. In contrast, in silico-predicted deuterium uptake levels of α6,7,8 from molecular dynamics (MD) simulations of the apo- and holo-PhuS structures are highly similar, consistent only with the holo-PhuS HDX-MS data. To rationalize this discrepancy between crystal structures, simulations, and observed HDX-MS, we exploit a recently developed computational approach (HDXer) that fits the relative weights of conformational populations within an ensemble of structures to conform to a target set of HDX-MS data. Here, a combination of enhanced sampling MD, HDXer, and dimensionality reduction analysis reveals an apo-PhuS conformational landscape in which α6, 7, and 8 are significantly rearranged compared to the crystal structure, including a loss of secondary structure in α6 and the displacement of α7 toward the HemO binding interface. Circular dichroism analysis confirms the loss of secondary structure, and the extracted ensembles of apo-PhuS and of heme-transfer-impaired H212R mutant, are consistent with known heme binding and transfer properties. The proposed conformational landscape provides structural insights into the modulation by heme of the dual function of PhuS.

Hydrogen-deuterium exchange combined with mass spectrometry (HDX-MS) is a widely applied biophysical technique that probes the structure and dynamics of biomolecules without the need for site-directed modifications or bio-orthogonal labels. The mechanistic interpretation of HDX data, however, is often qualitative and subjective, owing to a lack of quantitative methods to rigorously translate observed deuteration levels into atomistic structural information. To help address this problem, we have developed a methodology to generate structural ensembles that faithfully reproduce HDX-MS measurements. In this approach, an ensemble of protein conformations is first generated, typically using molecular dynamics simulations. A maximum-entropy bias is then applied post hoc to the resulting ensemble such that averaged peptide-deuteration levels, as predicted by an empirical model, agree with target values within a given level of uncertainty. We evaluate this approach, referred to as HDX ensemble reweighting (HDXer), for artificial target data reflecting the two major conformational states of a binding protein. We demonstrate that the information provided by HDX-MS experiments and by the model of exchange are sufficient to recover correctly weighted structural ensembles from simulations, even when the relevant conformations are rarely observed. Degrading the information content of the target data—e.g., by reducing sequence coverage, by averaging exchange levels over longer peptide segments, or by incorporating different sources of uncertainty—reduces the structural accuracy of the reweighted ensemble but still allows for useful insights into the distinctive structural features reflected by the target data. Finally, we describe a quantitative metric to rank candidate structural ensembles according to their correspondence with target data and illustrate the use of HDXer to describe changes in the conformational ensemble of the membrane protein LeuT. In summary, HDXer is designed to facilitate objective structural interpretations of HDX-MS data and to inform experimental approaches and further developments of theoretical exchange models.

Voltage-dependence of Na+-coupled phosphate cotransporters of the SLC34 family arises from displacement of charges intrinsic to the protein and the binding/release of one Na+ion in response to changes in the transmembrane electric field. Candidate coordination residues for the cation at the Na1 site were previously predicted by structural modeling using the x-ray structure of dicarboxylate transporter VcINDY as template and confirmed by functional studies. Mutations at Na1 resulted in altered steady-state and presteady-state characteristics that should be mirrored in the conformational changes induced by membrane potential changes. To test this hypothesis by functional analysis, double mutants of the flounder SLC34A2 protein were constructed that contain one of the Na1-site perturbing mutations together with a substituted cysteine for fluorophore labeling, as expressed in Xenopusoocytes. The locations of the mutations were mapped onto a homology model of the flounder protein. The effects of the mutagenesis were characterized by steady-state, presteady-state, and fluorometric assays. Changes in fluorescence intensity (ΔF) in response to membrane potential steps were resolved at three previously identified positions. These fluorescence data corroborated the altered presteady-state kinetics upon perturbation of Na1, and furthermore indicated concomitant changes in the microenvironment of the respective fluorophores, as evidenced by changes in the voltage dependence and time course of ΔF. Moreover, iodide quenching experiments indicated that the aqueous nature of the fluorophore microenvironment depended on the membrane potential. These findings provide compelling evidence that membrane potential and cation interactions induce significant large-scale structural rearrangements of the protein.

Secondary transporters use alternating-access mechanisms to couple uphill substrate movement to downhill ion flux. Most known transporters use a 'rocking bundle' motion, wherein the protein moves around an immobile substrate-binding site. However, the glutamate-transporter homolog GltPhtranslocates its substrate-binding site vertically across the membrane, through an 'elevator' mechanism. Here, we used the 'repeat swap' approach to computationally predict the outward-facing state of the Na+/succinate transporter VcINDY, from Vibrio cholerae. Our model predicts a substantial elevator-like movement of VcINDY's substrate-binding site, with a vertical translation of ~15 Å and a rotation of ~43°. Our observation that multiple disulfide cross-links completely inhibit transport provides experimental confirmation of the model and demonstrates that such movement is essential. In contrast, cross-links across the VcINDY dimer interface preserve transport, thus revealing an absence of large-scale coupling between protomers.

Phosphate plays essential biological roles and its plasma level in humans requires tight control to avoid bone loss (insufficiency) or vascular calcification (excess). Intestinal absorption and renal reabsorption of phosphate are mediated by members of the SLC34 family of sodium-coupled transporters (NaPi-IIa,b,c) whose membrane expression is regulated by various hormones, circulating proteins, and phosphate itself. Consequently, NaPi-II proteins are also potentially important pharmaceutical targets for controlling phosphate levels. Their crucial role in Pihomeostasis is underscored by pathologies resulting from naturally occurring SLC34 mutations and SLC34 knockout animals. SLC34 isoforms have been extensively studied with respect to transport mechanism and structure-function relationships; however, the three-dimensional structure is unknown. All SLC34 transporters share a duplicated motif comprising a glutamine followed by a stretch of threonine or serine residues, suggesting the presence of structural repeats as found in other transporter families. Nevertheless, standard bioinformatic approaches fail to clearly identify a suitable template for molecular modeling. Here, we used hydrophobicity profiles and hidden Markov models to define a structural repeat common to all SLC34 isoforms. Similar approaches identify a relationship with the core regions in a crystal structure of Vibrio choleraeNa+-dicarboxylate transporter VcINDY, from which we generated a homology model of human NaPi-IIa. The aforementioned SLC34 motifs in each repeat localize to the center of the model, and were predicted to form Na+and Picoordination sites. Functional relevance of key amino acids was confirmed by biochemical and electrophysiological analysis of expressed, mutated transporters. Moreover, the validity of the predicted architecture is corroborated by extensive published structure-function studies. The model provides key information for elucidating the transport mechanism and predicts candidate substrate binding sites.

As new atomic structures of membrane proteins are resolved, they reveal increasingly complex transmembrane topologies and highly irregular surfaces with crevices and pores. In many cases, specific interactions formed with the lipid membrane are functionally crucial, as is the overall lipid composition. Compounded with increasing protein size, these characteristics pose a challenge for the construction of simulation models of membrane proteins in lipid environments; clearly, that these models are sufficiently realistic bears upon the reliability of simulation-based studies of these systems. Here, we introduce GRIFFIN (GRId-based Force Field INput), which uses a versatile framework to automate and improve a widely used membrane-embedding protocol. Initially, GRIFFIN carves out lipid and water molecules from a volume equivalent to that of the protein, to conserve the system density. In the subsequent optimization phase GRIFFIN adds an implicit grid-based protein force field to a molecular dynamics simulation of the precarved membrane. In this force field, atoms inside the implicit protein volume experience an outward force that will expel them from that volume, whereas those outside are subject to electrostatic and van der Waals interactions with the implicit protein. At each step of the simulation, these forces are updated by GRIFFIN and combined with the intermolecular forces of the explicit lipid−water system. This procedure enables the construction of realistic and reproducible starting configurations of the protein−membrane interface within a reasonable time frame and with minimal intervention. GRIFFIN is a stand-alone tool designed to work alongside any existing molecular dynamics package, such as NAMD or GROMACS.

In this study, we investigate the extent to which techniques for homology modeling that were developed for water-soluble proteins are appropriate for membrane proteins as well. To this end we present an assessment of current strategies for homology modeling of membrane proteins and introduce a benchmark data set of homologous membrane protein structures, called HOMEP. First, we use HOMEP to reveal the relationship between sequence identity and structural similarity in membrane proteins. This analysis indicates that homology modeling is at least as applicable to membrane proteins as it is to water-soluble proteins and that acceptable models (with Cα-RMSD values to the native of 2Å or less in the transmembrane regions) may be obtained for template sequence identities of 30% or higher if an accurate alignment of the sequences is used. Second, we show that secondary-structure prediction algorithms that were developed for water-soluble proteins perform approximately as well for membrane proteins. Third, we provide a comparison of a set of commonly used sequence alignment algorithms as applied to membrane proteins. We find that high-accuracy alignments of membrane protein sequences can be obtained using state-of-the-art profile-to-profile methods that were developed for water-soluble proteins. Improvements are observed when weights derived from the secondary structure of the query and the template are used in the scoring of the alignment, a result which relies on the accuracy of the secondary-structure prediction of the query sequence. The most accurate alignments were obtained using template profiles constructed with the aid of structural alignments. In contrast, a simple sequence-to-sequence alignment algorithm, using a membrane protein-specific substitution matrix, shows no improvement in alignment accuracy. We suggest that profile-to-profile alignment methods should be adopted to maximize the accuracy of homology models of membrane proteins.

This article presents the results of a detailed analysis of helix-helix interactions in membrane and soluble proteins. A data set of interacting pairs of helices in membrane proteins of known structure was constructed and a structure alignment algorithm was used to identify pairs of helices in soluble proteins that superimpose well with pairs of helices in the membrane-protein data set. Most helix pairs in membrane proteins are found to have a significant number of structural homologs in soluble proteins, although in some cases, primarily involving irregular helices, no close homologs exist. An analysis of geometric relationships between interacting helices in the two sets of proteins identifies some differences in the distributions of helix length, interfacial area, packing angle, and distance between the polypeptide backbones. However, a subset of soluble-protein helix pairs that are close structural homologs to membrane-protein helix pairs exhibits distributions that mirror those observed in membrane proteins. The larger average interface size and smaller distance of closest approach seen for helices in membrane proteins appears due in part to a relative enrichment of alanines and glycines, particularly as components of the AxxxA and GxxxG motifs. It is argued that membrane helices are not on average more tightly packed than helices in soluble proteins; they are simply able to approach each other more closely. This enables them to interact over longer distances, which may in turn facilitate their remaining in contact over much of the width of the lipid bilayer. The close structural similarity seen between some pairs of helices in membrane and soluble proteins suggests that packing patterns observed in soluble proteins may be useful in the modeling of membrane proteins. Moreover, there do not appear to be fundamental differences between the magnitude of the forces that drive helix packing in membrane and soluble proteins, suggesting that strategies to make membrane proteins more soluble by mutating surface residues are likely to encounter success, at least in some cases.

Nanosecond molecular dynamics simulations in a fully solvated phospholipid bilayer have been performed on single transmembrane α-helices from three putative ion channel proteins encoded by viruses: NB (from influenza B), CM2 (from influenza C), and Vpu (from HIV-1). α-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 α-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. © 2000 John Wiley & Sons, Inc. Biopoly 53: 529–538, 2000

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.

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°, 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 pKAs 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. © 2000 Wiley-Liss, Inc.

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°, 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 pKAs 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. © 2000 Wiley‐Liss, Inc.

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°, 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 pKAs 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. © 2000 Wiley-Liss, Inc.

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°, 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 pKAs 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. © 2000 Wiley-Liss, Inc.

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 ∼4000 water molecules, yielding ∼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 10ns. 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 (Kerr et al., 1996, 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 (<10ps) and long (>100ps) time scales.

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 (Randa et al., 1999, Biophys. J. 77:2400–2410).

Molecular dynamics calculations were carried out on models of two synthetic leucine-serine ion channels: a tetrameric bundle with sequence (LSLLLSL)3NH2 and a hexameric bundle with sequence (LSSLLSL)3NH2. 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 4ns 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.

Integral membrane proteins containing at least one transmembrane (TM) α-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 ∼3500–4700 waters, giving a total of ∼18,000–21,000 atoms. Two simulations, each of 2ns 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 11ns. Analysis of the time-dependent secondary structure of the TM segments was used to define the regions that adopted a stable α-helical conformation throughout the simulation. This analysis indicates a core TM region of ∼20 residues (from residue 22 to residue 43) that remained in an α-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.