Search

Recently developed DNA-based analogues of membrane proteins have advanced synthetic biology. A fundamental question is how hydrophilic nanostructures reside in the hydrophobic environment of the membrane. Here, we use multiscale molecular dynamics (MD) simulations to explore the structure, stability and dynamics of an archetypical DNA nanotube inserted via a ring of membrane anchors into a phospholipid bilayer. Coarse-grained MD reveals that the lipids reorganize locally to interact closely with the membrane-spanning section of the DNA tube. Steered simulations along the bilayer normal establish the metastable nature of the inserted pore, yielding a force profile with barriers for membrane exit due to the membrane anchors. Atomistic, equilibrium simulations at two salt concentrations confirm the close packing of lipid around of the stably inserted DNA pore and its cation selectivity, while revealing localized structural fluctuations. The wide-ranging and detailed insight informs the design of next-generation DNA pores for synthetic biology or biomedicine.

Neutron and X-ray reflectivity of model membranes is increasingly used as a tool for the study of membrane structures and dynamics. As the systems under study become more complex, and as long, all-atom molecular-dynamics (MD) simulations of membranes become more available, there is increasing interest in the use of MD simulations in the analysis of reflectometry data from membranes. In order to perform this, it is necessary to produce a model of the complete interface, including not only the MD-derived structure of the membrane, but also the supporting substrate and any other interfacial layers that may be present. Here, it is shown that this is best performed by first producing a model of the occupied volume across the entire interface, and then converting this into a scattering length density (SLD) profile, rather than by splicing together the separate SLD profiles from the substrate layers and the membrane, since the latter approach can lead to discontinuities in the SLD profile and subsequent artefacts in the reflectivity calculation. It is also shown how the MD-derived membrane structure should be corrected to account for lower than optimal coverage and out-of-plane membrane fluctuations. Finally, the method of including the entire membrane structure in the reflectivity calculation is compared with an alternative approach in which the membrane components are approximated by functional forms, with only the component volumes being extracted from the simulation. It is shown that using only the fragment volumes is insufficient for a typical neutron data set of a single deuteration measured at several water contrasts, and that either weighting the model by including more structural information from the fit, or a larger data set involving a range of deuterations, are required to satisfactorily define the problem.

Characterizing the nanoscale dynamic organization within lipid bilayer membranes is central to our understanding of cell membranes at a molecular level. We investigate phase separation and communication across leaflets in ternary lipid bilayers, including saturated lipids with between 12 and 20 carbons per tail. Coarse-grained molecular dynamics simulations reveal a novel two-step kinetics due to hydrophobic mismatch, in which the initial response of the apposed leaflets upon quenching is to increase local asymmetry (antiregistration), followed by dominance of symmetry (registration) as the bilayer equilibrates. Antiregistration can become thermodynamically preferred if domain size is restricted below ∼20 nm, with implications for the symmetry of rafts and nanoclusters in cell membranes, which have similar reported sizes. We relate our findings to theory derived from a semimicroscopic model in which the leaflets experience a "direct" area-dependent coupling, and an "indirect" coupling that arises from hydrophobic mismatch and is most important at domain boundaries. Registered phases differ in composition from antiregistered phases, consistent with a direct coupling between the leaflets. Increased hydrophobic mismatch purifies the phases, suggesting that it contributes to the molecule-level lipid immiscibility. Our results demonstrate an interplay of competing interleaflet couplings that affect phase compositions and kinetics, and lead to a length scale that can influence lateral and transverse bilayer organization within cells.

Lipid molecules can bind to specific sites on integral membrane proteins, modulating their structure and function. We have undertaken coarse-grained simulations to calculate free energy profiles for glycolipids and phospholipids interacting with modulatory sites on the transmembrane helix dimer of the EGF receptor within a lipid bilayer environment. We identify lipid interaction sites at each end of the transmembrane domain and compute interaction free energy profiles for lipids with these sites. Interaction free energies ranged from ca. -40 to -4 kJ/mol for different lipid species. Those lipids (glycolipid GM3 and phosphoinositide PIP2) known to modulate EGFR function exhibit the strongest binding to interaction sites on the EGFR, and we are able to reproduce the preference for interaction with GM3 over other glycolipids suggested by experiment. Mutation of amino acid residues essential for EGFR function reduce the binding free energy of these key lipid species. The residues interacting with the lipids in the simulations are in agreement with those suggested by experimental (mutational) studies. This approach provides a generalizable tool for characterizing the interactions of lipids that bind to specific sites on integral membrane proteins.

Viruses typically pack their genetic material within a protein capsid. Enveloped viruses also have an outer membrane made up of a lipid bilayer and membrane-spanning glycoproteins. X-ray diffraction and cryoelectron microscopy provide high resolution static views of viral structure. Molecular dynamics (MD) simulations may be used to provide dynamic insights into the structures of viruses and their components. There have been a number of simulations of viral capsids and (in some cases) of the inner core of RNA or DNA packaged within them. These simulations have generally focussed on the structural integrity and stability of the capsid and/or on the influence of the nucleic acid core on capsid stability. More recently there have been a number of simulation studies of enveloped viruses, including HIV-1, influenza A, and dengue virus. These have addressed the dynamic behaviour of the capsid, the matrix, and/or of the outer envelope. Analysis of the dynamics of the lipid bilayer components of the envelopes of influenza A and of dengue virus reveals a degree of biophysical robustness, which may contribute to the stability of virus particles in different environments. Significant computational challenges need to be addressed to aid simulation of complex viruses and their membranes, including the need to integrate structural data from a range of sources to enable us to move towards simulations of intact virions. This article is part of a Special Issue entitled: Membrane Proteins edited by J.C. Gumbart and Sergei Noskov., (Copyright © 2016 The Authors. Published by Elsevier B.V. All rights reserved.)

Understanding the energetics of peripheral protein-membrane interactions is important to many areas of biophysical chemistry and cell biology. Estimating free-energy landscapes by molecular dynamics (MD) simulation is challenging for such systems, especially when membrane recognition involves complex lipids, e.g., phosphatidylinositol phosphates (PIPs). We combined coarse-grained MD simulations with umbrella sampling to quantify the binding of the well-explored GRP1 pleckstrin homology (PH) domain to model membranes containing PIP molecules. The experimentally observed preference of GRP1-PH for PIP3 over PIP2 was reproduced. Mutation of a key residue (K273A) within the canonical PIP-binding site significantly reduced the free energy of PIP binding. The presence of a noncanonical PIP-interaction site, observed experimentally in other PH domains but not previously in GRP1-PH, was also revealed. These studies demonstrate how combining coarse-grained simulations and umbrella sampling can unmask the molecular basis of the energetics of interactions between peripheral membrane proteins and complex cellular membranes.

A macromolecular nanopore inserted into a membrane may perturb the dynamic organization of the surrounding lipid bilayer. To better understand the nature of such perturbations, we have undertaken a systematic molecular dynamics simulation study of lipid bilayer structure and dynamics around three different classes of nanopore: a carbon nanotube, three related cyclic peptide nanotubes differing in the nature of their external surfaces, and a model of a β-barrel nanopore protein. Periodic spatial distributions of several lipid properties as a function of distance from the nanopore were observed. This was especially clear for the carbon nanotube system, for which the density of lipids, the bilayer thickness, the projection of lipid head-to-tail vectors onto the membrane plane, and lipid lateral diffusion coefficients exhibited undulatory behavior as a function of the distance from the surface of the channel. Overall, the differences in lipid behavior as a function of the nanopore structure reveal local adaptation of the bilayer structure and dynamics to different embedded nanopore structures. Both the local structure and dynamic behavior of lipids around membrane-embedded nanopores are sensitive to the geometry and nature of the outer surface of the macromolecule/molecular assembly forming the pore.

Ephs are transmembrane receptors that mediate cell-cell signaling. The N-terminal ectodomain binds ligands and enables receptor clustering, which activates the intracellular kinase. Relatively little is known about the function of the membrane-proximal fibronectin domain 2 (FN2) of the ectodomain. Multiscale molecular dynamics simulations reveal that FN2 interacts with lipid bilayers via a site comprising K441, R443, R465, Q462, S464, S491, W467, F490, and P459-461. FN2 preferentially binds anionic lipids, a preference that is reduced in the mutant K441E + R443E. We confirm these results by measuring the binding of wild-type and mutant FN2 domains to lipid vesicles. In simulations of the complete EphA2 ectodomain plus the transmembrane region, we show that FN2 anchors the otherwise flexible ectodomain at the surface of the bilayer. Altogether, our data suggest that FN2 serves a dual function of interacting with anionic lipids and constraining the structure of the EphA2 ectodomain to adopt membrane-proximal configurations., (Copyright © 2016 The Authors. Published by Elsevier Ltd.. All rights reserved.)

Glycolipids are key components of mammalian cell membranes, influencing a diverse range of cellular functions. For example, a number of receptor tyrosine kinases, including the epidermal growth factor receptor (EGFR), are allosterically regulated by the glycolipid monosialodihexosylganglioside (GM3). Recent advances in molecular dynamics methods, especially the development of coarse-grained models, have enabled simulations of increasingly complex models of cell membranes. We demonstrate these methodological developments via a case study of a coarse-grained model for the ganglioside GM3. This glycolipid is included in simulations of a mixed lipid bilayer model reflecting the compositional complexity of a mammalian cell membrane. The resultant membrane model is used to simulate the interactions of GM3 with the transmembrane domain of the EGFR., (Copyright © 2015 The Authors. Published by Elsevier B.V. All rights reserved.)

The resemblance of lipid membrane models to physiological membranes determines how well molecular dynamics (MD) simulations imitate the dynamic behavior of cell membranes and membrane proteins. Physiological lipid membranes are composed of multiple types of phospholipids, and the leaflet compositions are generally asymmetric. Here we describe an approach for self-assembly of a Coarse-Grained (CG) membrane model with physiological composition and leaflet asymmetry using the MARTINI force field. An initial set-up of two boxes with different types of lipids according to the leaflet asymmetry of mammalian cell membranes stacked with 0.5 nm overlap, reliably resulted in the self-assembly of bilayer membranes with leaflet asymmetry resembling that of physiological mammalian cell membranes. Self-assembly in the presence of a fragment of the plasma membrane protein syntaxin 1A led to spontaneous specific positioning of phosphatidylionositol(4,5)bisphosphate at a positively charged stretch of syntaxin consistent with experimental data. An analogous approach choosing an initial set-up with two concentric shells filled with different lipid types results in successful assembly of a spherical vesicle with asymmetric leaflet composition. Self-assembly of the vesicle in the presence of the synaptic vesicle protein synaptobrevin 2 revealed the correct position of the synaptobrevin transmembrane domain. This is the first CG MD method to form a membrane with physiological lipid composition as well as leaflet asymmetry by self-assembly and will enable unbiased studies of the incorporation and dynamics of membrane proteins in more realistic CG membrane models.

There has been exponential growth in the number of membrane protein structures determined. Nevertheless, these structures are usually resolved in the absence of their lipid environment. Coarse-grained molecular dynamics (CGMD) simulations enable insertion of membrane proteins into explicit models of lipid bilayers. We have automated the CGMD methodology, enabling membrane protein structures to be identified upon their release into the PDB and embedded into a membrane. The simulations are analyzed for protein-lipid interactions, identifying lipid binding sites, and revealing local bilayer deformations plus molecular access pathways within the membrane. The coarse-grained models of membrane protein/bilayer complexes are transformed to atomistic resolution for further analysis and simulation. Using this automated simulation pipeline, we have analyzed a number of recently determined membrane protein structures to predict their locations within a membrane, their lipid/protein interactions, and the functional implications of an enhanced understanding of the local membrane environment of each protein., (Copyright © 2015 The Authors. Published by Elsevier Ltd.. All rights reserved.)

The Escherichia coli UraA H+-uracil symporter is a member of the nucleobase/ascorbate transporter (NAT) family of proteins, and is responsible for the proton-driven uptake of uracil. Multiscale molecular dynamics simulations of the UraA symporter in phospholipid bilayers consisting of: 1) 1-palmitoyl 2-oleoyl-phosphatidylcholine (POPC); 2) 1-palmitoyl 2-oleoyl-phosphatidylethanolamine (POPE); and 3) a mixture of 75% POPE, 20% 1-palmitoyl 2-oleoyl-phosphatidylglycerol (POPG); and 5% 1-palmitoyl 2-oleoyl-diphosphatidylglycerol/cardiolipin (CL) to mimic the lipid composition of the bacterial inner membrane, were performed using the MARTINI coarse-grained force field to self-assemble lipids around the crystal structure of this membrane transport protein, followed by atomistic simulations. The overall fold of the protein in lipid bilayers remained similar to the crystal structure in detergent on the timescale of our simulations. Simulations were performed in the absence of uracil, and resulted in a closed state of the transporter, due to relative movement of the gate and core domains. Anionic lipids, including POPG and especially CL, were found to associate with UraA, involving interactions between specific basic residues in loop regions and phosphate oxygens of the CL head group. In particular, three CL binding sites were identified on UraA: two in the inner leaflet and a single site in the outer leaflet. Mutation of basic residues in the binding sites resulted in the loss of CL binding in the simulations. CL may play a role as a "proton trap" that channels protons to and from this transporter within CL-enriched areas of the inner bacterial membrane.

Structural studies of membrane proteins have highlighted the likely influence of membrane mimetic environments (i.e., lipid bilayers versus detergent micelles) on the conformation and dynamics of small α-helical membrane proteins. We have used molecular dynamics simulations to compare the conformational dynamics of BM2 (a small α-helical protein from the membrane of influenza B) in a model phospholipid bilayer environment with its behavior in protein-detergent complexes with either the zwitterionic detergent dihexanoylphosphatidylcholine (DHPC) or the nonionic detergent dodecylmaltoside (DDM). We find that DDM more closely resembles the lipid bilayer in terms of its interaction with the protein, while the short-tailed DHPC molecule forms "nonphysiological" interactions with the protein termini. We find that the intrinsic micelle properties of each detergent are conserved upon formation of the protein-detergent complex. This implies that simulations of detergent micelles may be used to help select optimal conditions for experimental studies of membrane proteins.

Cell membranes are complex multicomponent systems, which are highly heterogeneous in the lipid distribution and composition. To date, most molecular simulations have focussed on relatively simple lipid compositions, helping to inform our understanding of in vitro experimental studies. Here we describe on simulations of complex asymmetric plasma membrane model, which contains seven different lipids species including the glycolipid GM3 in the outer leaflet and the anionic lipid, phosphatidylinositol 4,5-bisphophate (PIP2), in the inner leaflet. Plasma membrane models consisting of 1500 lipids and resembling the in vivo composition were constructed and simulations were run for 5 µs. In these simulations the most striking feature was the formation of nano-clusters of GM3 within the outer leaflet. In simulations of protein interactions within a plasma membrane model, GM3, PIP2, and cholesterol all formed favorable interactions with the model α-helical protein. A larger scale simulation of a model plasma membrane containing 6000 lipid molecules revealed correlations between curvature of the bilayer surface and clustering of lipid molecules. In particular, the concave (when viewed from the extracellular side) regions of the bilayer surface were locally enriched in GM3. In summary, these simulations explore the nanoscale dynamics of model bilayers which mimic the in vivo lipid composition of mammalian plasma membranes, revealing emergent nanoscale membrane organization which may be coupled both to fluctuations in local membrane geometry and to interactions with proteins.

Many cellular signalling and related events are triggered by the association of peripheral proteins with anionic lipids in the cell membrane (e.g. phosphatidylinositol phosphates or PIPs). This association frequently occurs via lipid-binding modules, e.g. pleckstrin homology (PH), C2 and four-point-one, ezrin, radixin, moesin (FERM) domains, present in peripheral and cytosolic proteins. Multiscale simulation approaches that combine coarse-grained and atomistic MD simulations may now be applied with confidence to investigate the molecular mechanisms of the association of peripheral proteins with model bilayers. Comparisons with experimental data indicate that such simulations can predict specific peripheral protein-lipid interactions. We discuss the application of multiscale MD simulation and related approaches to investigate the association of peripheral proteins which contain PH, C2 or FERM-binding modules with lipid bilayers of differing phospholipid composition, including bilayers containing multiple PIP molecules.

Membrane proteins are frequently present in crowded environments, which favour lateral association and, on occasions, two-dimensional crystallization. To better understand the non-specific lateral association of a membrane protein we have characterized the free energy landscape for the dimerization of a bacterial outer membrane protein, NanC, in a phospholipid bilayer membrane. NanC is a member of the KdgM-family of bacterial outer membrane proteins and is responsible for sialic acid transport in E. coli. Umbrella sampling and coarse-grained molecular dynamics were employed to calculate the potentials of mean force (PMF) for a variety of restrained relative orientations of two NanC proteins as the separation of their centres of mass was varied. We found the free energy of dimerization for NanC to be in the range of -66 kJ mol(-1) to -45 kJ mol(-1). Differences in the depths of the PMFs for the various orientations are related to the shape of the proteins. This was quantified by calculating the lipid-inaccessible buried surface area of the proteins in the region around the minimum of each PMF. The depth of the potential well of the PMF was shown to depend approximately linearly on the buried surface area. We were able to resolve local minima in the restrained PMFs that would not be revealed using conventional umbrella sampling. In particular, these features reflected the local organization of the intervening lipids between the two interacting proteins. Through a comparison with the distribution of lipids around a single freely-diffusing NanC, we were able to predict the location of these restrained local minima for the orientational configuration in which they were most pronounced. Our ability to make this prediction highlights the important role that lipid organization plays in the association of two NanCs in a bilayer.

The Ras family of small membrane-associated GTP-ases are important components in many different cell signalling cascades. They are thought to cluster on the cell membrane through association with cholesterol-rich nanodomains. This process remains poorly understood. Here we test the effect of adding multiple copies of NRas, one of the canonical Ras proteins, to a three-component lipid bilayer that rapidly undergoes spinodal decomposition (i.e. unmixing), thereby creating ordered and disordered phases. Coarse-grained molecular dynamics simulations of a large bilayer containing 6000 lipids, with and without protein, are compared. NRas preferentially localises to the interface between the domains and slows the rate at which the domains grow. We infer that this doubly-lipidated cell signalling protein is reducing the line tension between the ordered and disordered regions. This analysis is facilitated by our use of techniques borrowed from image-processing. The conclusions above are contingent upon several assumptions, including the use of a model lipid with doubly unsaturated tails and the limited structural data available for the C-terminus of NRas, which is where the lipid anchors are found.

The twin-arginine translocase (Tat) carries out the remarkable process of translocating fully folded proteins across the cytoplasmic membrane of prokaryotes and the thylakoid membrane of plant chloroplasts. Tat is required for bacterial pathogenesis and for photosynthesis in plants. TatA, the protein-translocating element of the Tat system, is a small transmembrane protein that assembles into ring-like oligomers of variable size. We have determined a structural model of the Escherichia coli TatA complex in detergent solution by NMR. TatA assembly is mediated entirely by the transmembrane helix. The amphipathic helix extends outwards from the ring of transmembrane helices, permitting assembly of complexes with variable subunit numbers. Transmembrane residue Gln8 points inward, resulting in a short hydrophobic pore in the center of the complex. Simulations of the TatA complex in lipid bilayers indicate that the short transmembrane domain distorts the membrane. This finding suggests that TatA facilitates protein transport by sensitizing the membrane to transient rupture.

Phosphatase and tensin-homolog deleted on chromosome 10 (PTEN) is a tumor-suppressor protein that regulates phosphatidylinositol 3-kinase (PI3-K) signaling by binding to the plasma membrane and hydrolyzing the 3' phosphate from phosphatidylinositol (3,4,5)-trisphosphate (PI(3,4,5)P3) to form phosphatidylinositol (4,5)-bisphosphate (PI(4,5)P2). Several loss-of-function mutations in PTEN that impair lipid phosphatase activity and membrane binding are oncogenic, leading to the development of a variety of cancers, but information about the membrane-associated state of PTEN remains sparse. We have modeled a membrane-associated state of the truncated PTEN structure bound to PI(3,4,5)P3 via multiscale molecular dynamics simulations. We show that the location of the membrane-binding surface agrees with experimental observations and is robust to changes in lipid composition. The level of membrane interaction is substantially reduced in the phosphatase domain for the triple mutant R161E/K163E/K164E, in line with experimental results. We observe clustering of anionic lipids around the C2 domain in preference to the phosphatase domain, suggesting that the C2 domain is involved in nonspecific interactions with negatively charged lipid headgroups. Finally, our simulations suggest that the oncogenicity of the R335L mutation may be due to a reduction in the interaction of the mutant PTEN with anionic lipids., (Copyright © 2013 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

The antimicrobial peptide maculatin 1.1 (M1.1) is an amphipathic α-helix that permeabilizes lipid bilayers. In coarse-grained molecular dynamics (CG MD) simulations, M1.1 has previously been shown to form membrane-spanning aggregates in DPPC bilayers. In this study, a simple multiscale methodology has been applied to allow sampling of important regions of the free energy surface at higher resolution. Thus, by back-converting the CG configurations to atomistic representations, it is shown that water is able to permeate through the M1.1 aggregates. Investigation of aggregate stoichiometry shows that at least six peptides are required for water permeation. The aggregates are dynamically disordered structures, and water flux occurs through irregular, fluctuating channels. The results are discussed in relation to experimental data and other simulations of antimicrobial peptides.

Interactions between protein domains and lipid molecules play key roles in controlling cell membrane signalling and trafficking. The pleckstrin homology (PH) domain is one of the most widespread, binding specifically to phosphatidylinositol phosphates (PIPs) in cell membranes. PH domains must locate specific PIPs in the presence of a background of approximately 20% anionic lipids within the cytoplasmic leaflet of the plasma membrane. We investigate the mechanism of such recognition via a multiscale procedure combining Brownian dynamics (BD) and molecular dynamics (MD) simulations of the GRP1 PH domain interacting with phosphatidylinositol (3,4,5)-trisphosphate (PI(3,4,5)P₃). The interaction of GRP1-PH with PI(3,4,5)P₃ in a zwitterionic bilayer is compared with the interaction in bilayers containing different levels of anionic 'decoy' lipids. BD simulations reveal both translational and orientational electrostatic steering of the PH domain towards the PI(3,4,5)P₃-containing anionic bilayer surface. There is a payoff between non-PIP anionic lipids attracting the PH domain to the bilayer surface in a favourable orientation and their role as 'decoys', disrupting the interaction of GRP1-PH with the PI(3,4,5)P₃ molecule. Significantly, approximately 20% anionic lipid in the cytoplasmic leaflet of the bilayer is nearly optimal to both enhance orientational steering and to localise GRP1-PH proximal to the surface of the membrane without sacrificing its ability to locate PI(3,4,5)P₃ within the bilayer plane. Subsequent MD simulations reveal binding to PI(3,4,5)P₃, forming protein-phosphate contacts comparable to those in X-ray structures. These studies demonstrate a computational framework which addresses lipid recognition within a cell membrane environment, offering a link between structural and cell biological characterisation.

Dimerization of transmembrane (TM) α helices of membrane receptors plays a key role in signaling. We show that molecular dynamics simulations yield models of integrin TM helix heterodimers, which agree well with available NMR structures. We use a multiscale simulation approach, combining coarse-grained and subsequent atomistic simulation, to model the dimerization of wild-type (WT) and mutated sequences of the αIIb and β3 integrin TM helices. The WT helices formed a stable, right-handed dimer with the same helix-helix interface as in the published NMR structure (PDB: 2K9J). In contrast, the presence of disruptive mutations perturbed the interface between the helices, altering the conformational stability of the dimer. The αIIb/β3 interface was more flexible than that of, e.g., glycophorin A. This is suggestive of a role for alternative packing modes of the TM helices in transbilayer signaling., (Copyright © 2011 Elsevier Ltd. All rights reserved.)

Voltage sensors (VS) domains couple the activation of ion channels/enzymes to changes in membrane voltage. We used molecular dynamics simulations to examine interactions with lipids of several VS homologs. VSs in intact channels in the activated state are exposed to phospholipids, leading to a characteristic local distortion of the lipid bilayer which decreases its thickness by ∼10 Å. This effect is mediated by a conserved hydrophilic stretch in the S4-S5 segment linking the VS and the pore domains, and may favor gating charges crossing the membrane. In cationic lipid bilayers lacking phosphate groups, VSs form fewer contacts with lipid headgroups. The S3-S4 paddle motifs show persistent interactions of individual lipid molecules, influenced by the hairpin loop. In conclusion, our results suggest common interactions with phospholipids for various VS homologs, providing insights into the molecular basis of their stabilization in the membrane and how they are altered by lipid modification., (Copyright © 2011 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

Most membrane proteins contain a transmembrane (TM) domain made up of a bundle of lipid-bilayer-spanning α-helices. TM α-helices are generally composed of a core of largely hydrophobic amino acids, with basic and aromatic amino acids at each end of the helix forming interactions with the lipid headgroups and water. In contrast, the S4 helix of ion channel voltage sensor (VS) domains contains four or five basic (largely arginine) side chains along its length and yet adopts a TM orientation as part of an independently stable VS domain. Multiscale molecular dynamics simulations are used to explore how a charged TM S4 α-helix may be stabilized in a lipid bilayer, which is of relevance in the context of mechanisms of translocon-mediated insertion of S4. Free-energy profiles for insertion of the S4 helix into a phospholipid bilayer suggest that it is thermodynamically favorable for S4 to insert from water to the center of the membrane, where the helix adopts a TM orientation. This is consistent with crystal structures of Kv channels, biophysical studies of isolated VS domains in lipid bilayers, and studies of translocon-mediated S4 helix insertion. Decomposition of the free-energy profiles reveals the underlying physical basis for TM stability, whereby the preference of the hydrophobic residues of S4 to enter the bilayer dominates over the free-energy penalty for inserting charged residues, accompanied by local distortion of the bilayer and penetration of waters. We show that the unique combination of charged and hydrophobic residues in S4 allows it to insert stably into the membrane., (Copyright © 2011 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

Free energy profiles for insertion of a hydrophobic transmembrane protein α-helix (M2 from CFTR) into a lipid bilayer have been calculated using coarse-grained molecular dynamics simulations and umbrella sampling to yield potentials of mean force along a reaction path corresponding to translation of a helix across a lipid bilayer. The calculated free energy of insertion is smaller when a bilayer with a thinner hydrophobic region is used. The free energies of insertion from the potentials of mean force are compared with those derived from a number of hydrophobicity scales and with those derived from translocon-mediated insertion. This comparison supports recent models of translocon-mediated insertion and in particular suggests that: 1), helices in an about-to-be-inserted state may be located in a hydrophobic region somewhat thinner than the core of a lipid bilayer; and/or 2), helices in a not-to-be-inserted state may experience an environment more akin (e.g., in polarity/hydrophobicity) to the bilayer/water interface than to bulk water., (Copyright © 2010 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

Integrins are cell surface receptors crucial for cell migration and adhesion. They are activated by interactions of the talin head domain with the membrane surface and the integrin β cytoplasmic tail. Here, we use coarse-grained molecular dynamic simulations and nuclear magnetic resonance spectroscopy to elucidate the membrane-binding surfaces of the talin head (F2-F3) domain. In particular, we show that mutations in the four basic residues (K258E, K274E, R276E, and K280E) in the F2 binding surface reduce the affinity of the F2-F3 for the membrane and modify its orientation relative to the bilayer. Our results highlight the key role of anionic lipids in talin/membrane interactions. Simulation of the F2-F3 in complex with the α/β transmembrane dimer reveals information for its orientation relative to the membrane. Our studies suggest that the perturbed orientation of talin relative to the membrane in the F2 mutant would be expected to in turn perturb talin/integrin interactions., (Copyright © 2010 Elsevier Ltd. All rights reserved.)

A set of 49 protein nanopore-lipid bilayer systems was explored by means of coarse-grained molecular-dynamics simulations to study the interactions between nanopores and the lipid bilayers in which they are embedded. The seven nanopore species investigated represent the two main structural classes of membrane proteins (alpha-helical and beta-barrel), and the seven different bilayer systems range in thickness from approximately 28 to approximately 43 A. The study focuses on the local effects of hydrophobic mismatch between the nanopore and the lipid bilayer. The effects of nanopore insertion on lipid bilayer thickness, the dependence between hydrophobic thickness and the observed nanopore tilt angle, and the local distribution of lipid types around a nanopore in mixed-lipid bilayers are all analyzed. Different behavior for nanopores of similar hydrophobic length but different geometry is observed. The local lipid bilayer perturbation caused by the inserted nanopores suggests possible mechanisms for both lipid bilayer-induced protein sorting and protein-induced lipid sorting. A correlation between smaller lipid bilayer thickness (larger hydrophobic mismatch) and larger nanopore tilt angle is observed and, in the case of larger hydrophobic mismatches, the simulated tilt angle distribution seems to broaden. Furthermore, both nanopore size and key residue types (e.g., tryptophan) seem to influence the level of protein tilt, emphasizing the reciprocal nature of nanopore-lipid bilayer interactions.

Coarse-grained molecular dynamics simulations have been used to explore the interactions of C(60) and its derivatives with lipid bilayers. Pristine C(60) partitions into the bilayer core, whilst C(60)(OH)(20) experiences a central energetic barrier to permeation across the bilayer. For intermediate levels of derivatization, e.g. C(60)(OH)(10), this central barrier is smaller and there is an energetic well at the bilayer/water interface, thus promoting entry into cells via bilayer permeation whilst maintaining solubility in water.

Monotopic membrane proteins bind tightly to cell membranes but do not generally span the lipid bilayer. Their interactions with lipid bilayers may be studied via coarse-grained molecular dynamics (CG-MD) simulations. Understanding such interactions is important as monotopic enzymes frequently act on hydrophobic substrates, while X-ray structures rarely provide direct information about their interactions with membranes. CG-MD self-assembly simulations enable prediction of the orientation and depth of insertion into a lipid bilayer of a monotopic protein, and also of the interactions of individual protein residues with lipid molecules. The CG-MD method has been evaluated via comparison with extended (>30 ns) atomistic simulations of monoamine oxidase, revealing good agreement between the results of coarse-grained and atomistic simulations. CG-MD simulations have been applied to a set of 11 monotopic proteins for which three-dimensional structures are available. These proteins may be divided into two groups on the basis of the results of the simulations. One group consists of those proteins which are inserted into the lipid bilayer to a limited extent, interacting mainly at the phospholipid-water interface. The second group consists of those which are inserted more deeply into the bilayer. Those monotopic proteins which are inserted more deeply cause significant local perturbation of bilayer properties such as bilayer thickness. Deeper insertion seems to correlate with a greater number of basic residues in the "foot" whereby a monotopic protein interacts with the membrane.

DNA-lipid complexes are of biomedical importance as delivery vectors for gene therapy. To gain insight into the interactions of DNA with zwitterionic and cationic (dimyristoyltrimethylammonium propane (DMTAP)) lipids, we have used coarse-grained molecular dynamics simulations to study the self-assembly of DPPC and DPPC/DMTAP lipid bilayers in the presence of a DNA dodecamer. We observed the spontaneous formation of lipid bilayers from initial systems containing randomly placed lipids, water-counterions and DNA. In both the DPPC and DPPC/DMTAP simulations, the DNA molecule is located at the water-lipid headgroup interface, lying approximately parallel to the plane of the bilayer. We have also calculated the potential of mean force for transferring a DNA dodecamer through a DPPC/DMTAP bilayer. A high energetic barrier to DNA insertion into the hydrophobic core of the bilayer is observed. The DNA adopts a transmembrane orientation only in this region. Local bilayer deformation in the vicinity of the DNA molecule is observed, largely as a result of the DNA-DMTAP headgroup attraction.

Membrane protein function and stability has been shown to be dependent on the lipid environment. Recently, we developed a high-throughput computational approach for the prediction of membrane protein/lipid interactions. In the current study, we enhanced this approach with the addition of a new measure of the distortion caused by membrane proteins on a lipid bilayer. This is illustrated by considering the effect of lipid tail length and headgroup charge on the distortion caused by the integral membrane proteins MscS and FLAP, and by the voltage sensing domain from the channel KvAP. Changing the chain length of lipids alters the extent but not the pattern of distortion caused by MscS and FLAP; lipid headgroups distort in order to interact with very similar but not identical regions in these proteins for all bilayer widths investigated. Introducing anionic lipids into a DPPC bilayer containing the KvAP voltage sensor does not affect the extent of bilayer distortion.

Experimental and computational studies have indicated that hydrophobicity plays a key role in driving the insertion of transmembrane alpha-helices into lipid bilayers. Molecular dynamics simulations allow exploration of the nature of the interactions of transmembrane alpha-helices with their lipid bilayer environment. In particular, coarse-grained simulations have considerable potential for studying many aspects of membrane proteins, ranging from their self-assembly to the relation between their structure and function. However, there is a need to evaluate the accuracy of coarse-grained estimates of the energetics of transmembrane helix insertion. Here, three levels of complexity of model system have been explored to enable such an evaluation. First, calculated free energies of partitioning of amino acid side chains between water and alkane yielded an excellent correlation with experiment. Second, free energy profiles for transfer of amino acid side chains along the normal to a phosphatidylcholine bilayer were in good agreement with experimental and atomistic simulation studies. Third, estimation of the free energy profile for transfer of an arginine residue, embedded within a hydrophobic alpha-helix, to the center of a lipid bilayer gave a barrier of approximately 15 kT. Hence, there is a substantial barrier to membrane insertion for charged amino acids, but the coarse-grained model still underestimates the corresponding free energy estimate (approximately 29 kT) from atomistic simulations (Dorairaj, S., and Allen, T. W. (2007) Proc. Natl. Acad. Sci. U.S.A. 104, 4943-4948). Coarse-grained simulations were then used to predict the free energy profile for transfer of a simple model transmembrane alpha-helix (WALP23) across a lipid bilayer. The results indicated that a transmembrane orientation was favored by about -70 kT.

A number of membrane proteins act via binding at the water/lipid bilayer interface. An important example of such proteins is provided by the gating-modifier toxins that act on voltage-gated potassium (Kv) channels. They are thought to partition to the headgroup region of lipid bilayers, and so provide a good system for probing the nature of interactions of a protein with the water/bilayer interface. We used coarse-grained molecular dynamics simulations to compute the one-dimensional potential of mean force (i.e., free energy) profile that governs the interaction between a Kv channel gating-modifier toxin (VSTx1) and model phospholipid bilayers. The reaction coordinate sampled corresponds to the position of the toxin along the bilayer normal. The course-grained representation of the protein and lipids enabled us to explore extended time periods, revealing aspects of toxin/bilayer dynamics and energetics that would be difficult to observe on the timescales currently afforded by atomistic molecular dynamics simulations. In particular, we show for this model system that the bilayer deforms as it interacts with the toxin, and that such deformations perturb the free energy profile. Bilayer deformation therefore adds an additional layer of complexity to be addressed in investigations of membrane/protein systems. In particular, one should allow for local deformations that may arise due to the spatial array of charged and hydrophobic elements of an interfacially located membrane protein.

Maculatin 1.1 (M1.1) is a membrane-active antimicrobial peptide (AMP) from an Australian tree frog that forms a kinked amphipathic alpha-helix in the presence of a lipid bilayer or bilayer-mimetic environment. To help elucidate its mechanism of membrane-lytic activity, we performed a total of approximately 8 micros of coarse-grained molecular dynamics (CG-MD) simulations of M1.1 in the presence of zwitterionic phospholipid membranes. Several systems were simulated in which the peptide/lipid ratio was varied. At a low peptide/lipid ratio, M1.1 adopted a kinked, membrane-interfacial location, consistent with experiment. At higher peptide/lipid ratios, we observed spontaneous, cooperative membrane insertion of M1.1 peptide aggregates. The minimum size for formation of a transmembrane (TM) aggregate was just four peptides. The absence of a simple and well-defined central channel, along with the exclusion of lipid headgroups from the aggregates, suggests that a pore-like model is an unlikely explanation for the mechanism of membrane lysis by M1.1. We also performed an extended 1.25 micros simulation of the permeabilization of a complete liposome by multiple peptides. Consistent with the simpler bilayer simulations, formation of monomeric interfacial peptides and TM peptide clusters was observed. In contrast, major structural changes were observed in the vesicle membrane, implicating induced membrane curvature in the mechanism of active antimicrobial peptide lysis. This contrasted with the behavior of the nonpore-forming model peptide WALP23, which inserted into the vesicle to form extended clusters of TM alpha-helices with relatively little perturbation of bilayer properties.

Carbon nanotubes (CNTs) are possible nanoinjectors for the introduction of therapeutic agents into cells. To explore their interactions with a lipid bilayer membrane and to model the nanoinjection process, we used coarse-grained molecular dynamics to simulate the penetration of dipalmitoylphosphatidylcholine (DPPC) bilayers by single-walled CNTs. Lipids are extracted from a bilayer during CNT penetration and reside on both the inner and the outer tube surfaces. Lipids that interact with the CNT interior wall spread out and hence can "block" the tube. However, the degree of lipid lining of the inner surface is strongly dependent upon the tube penetration velocity, with fewer lipids extracted from the bilayer at higher rates. There is no apparent effect on bilayer integrity after CNT penetration, with the bilayer able to self-seal. Our findings reveal some of the complexities of the interactions of lipids with CNT nanoinjectors and suggest a need to further characterize the influence of, for example, CNT functionalization and cargo on lipid blocking of CNTs.

A number of membrane-active enzymes act in a complex environment formed by the interface between a lipid bilayer and bulk water. Although x-ray diffraction studies yield structures of isolated enzyme molecules, a detailed characterization of their interactions with the interface requires a measure of how deeply such a membrane-associated protein penetrates into a lipid bilayer. Here, we apply coarse-grained (CG) molecular dynamics (MD) simulations to probe the interaction of porcine pancreatic phospholipase A2 (PLA2) with a lipid bilayer containing palmitoyl-oleoyl-phosphatidyl choline and palmitoyl-oleoyl-phosphatidyl glycerol molecules. We also used a configuration from a CG-MD trajectory to initiate two atomistic (AT) MD simulations. The results of the CG and AT simulations are evaluated by comparison with available experimental data. The membrane-binding surface of PLA2 consists of a patch of hydrophobic residues surrounded by polar and basic residues. We show this proposed footprint interacts preferentially with the anionic headgroups of the palmitoyl-oleoyl-phosphatidyl glycerol molecules. Thus, both electrostatic and hydrophobic interactions determine the location of PLA2 relative to the bilayer. From a general perspective, this study demonstrates that CG-MD simulations may be used to reveal the orientation and location of a membrane-surface-bound protein relative to a lipid bilayer, which may subsequently be refined by AT-MD simulations to probe more detailed interactions.

Complete determination of a membrane protein structure requires knowledge of the protein position within the lipid bilayer. As the number of determined structures of membrane proteins increases so does the need for computational methods which predict their position in the lipid bilayer. Here we present a coarse-grained molecular dynamics approach to lipid bilayer self-assembly around membrane proteins. We demonstrate that this method can be used to predict accurately the protein position in the bilayer for membrane proteins with a range of different sizes and architectures.

The structures of three bacterial outer membrane proteins (OmpA, OmpX and PagP) have been determined by both X-ray diffraction and NMR. We have used multiple (7 x 15 ns) MD simulations to compare the conformational dynamics resulting from the X-ray versus the NMR structures, each protein being simulated in a lipid (DMPC) bilayer. Conformational drift was assessed via calculation of the root mean square deviation as a function of time. On this basis the 'quality' of the starting structure seems mainly to influence the simulation stability of the transmembrane beta-barrel domain. Root mean square fluctuations were used to compare simulation mobility as a function of residue number. The resultant residue mobility profiles were qualitatively similar for the corresponding X-ray and NMR structure-based simulations. However, all three proteins were generally more mobile in the NMR-based than in the X-ray simulations. Principal components analysis was used to identify the dominant motions within each simulation. The first two eigenvectors (which account for >50% of the protein motion) reveal that such motions are concentrated in the extracellular loops and, in the case of PagP, in the N-terminal alpha-helix. Residue profiles of the magnitude of motions corresponding to the first two eigenvectors are similar for the corresponding X-ray and NMR simulations, but the directions of these motions correlate poorly reflecting incomplete sampling on a approximately 10 ns timescale.

Molecular dynamics simulations of lactose permease (LacY) in a phospholipid bilayer reveal the conformational dynamics of the protein. In inhibitor-bound simulations (i.e., those closest to the X-ray structure) the protein was stable, showing little conformational change over a 50 ns timescale. Movement of the bound inhibitor, TDG, to an alternative binding mode was observed, so that it interacted predominantly with the N-terminal domain and with residue E269 from the C-terminal domain. In multiple ligand-free simulations, a degree of domain closure occurred. This switched LacY to a state with a central cavity closed at both the intracellular and periplasmic ends. This may resemble a possible intermediate in the transport mechanism. Domain closure occurs by a combination of rigid-body movements of domains and of intradomain motions of helices, especially TM4, TM5, TM10, and TM11. A degree of intrahelix flexibility appears to be important in the conformational change.

An implicit-membrane representation based on the generalized Born theory of solvation has been developed. The method was parameterized against the water-to-cyclohexane insertion free energies of hydrophobic side-chain analogs. Subsequently, the membrane was compared with experimental data from translocon inserted polypeptides and validated by comparison with an independent dataset of six membrane-associated peptides and eight integral membrane proteins of known structure and orientation. Comparison of the insertion energy of alpha-helical model peptides with the experimental values from the biological hydrophobicity scale of Hessa et al. gave a correlation of 93% with a mean unsigned error of 0.64 kcal/mol, when charged residues were ignored. The membrane insertion energy was found to be dependent on residue position. This effect is particularly pronounced for charged and polar residues, which strongly prefer interfacial locations. All integral membrane proteins investigated orient and insert correctly into the implicit-membrane model. Remarkably, the membrane model correctly predicts a partially inserted configuration for the monotopic membrane protein cyclooxygenase, matching experimental and theoretical predictions. To test the applicability and usefulness of the implicit-membrane method, molecular simulations of influenza A M2 as well as the glycophorin A dimer were performed. Both systems remain structurally stable and integrated into the membrane.

Monotopic proteins make up a class of membrane proteins that bind tightly to, but do not span, cell membranes. We examine and compare how two monotopic proteins, monoamine oxidase B (MAO-B) and cyclooxygenase-2 (COX-2), interact with a phospholipid bilayer using molecular dynamics simulations. Both enzymes form between three and seven hydrogen bonds with the bilayer in our simulations with basic side chains accounting for the majority of these interactions. By analyzing lipid order parameters, we show that, to a first approximation, COX-2 disrupts only the upper leaflet of the bilayer. In contrast, the top and bottom halves of the lipid tails surrounding MAO-B are more and less ordered, respectively, than in the absence of the protein. Finally, we identify which residues are important in binding individual phospholipids by counting the number and type of lipid atoms that come close to each amino acid residue. The existing models that explain how these proteins bind to bilayers were proposed following inspection of the X-ray crystallographic structures. Our results support these models and suggest that basic residues contribute significantly to the binding of these monotopic proteins to bilayers through the formation of hydrogen bonds with phospholipids.

Coarse-grained molecular dynamics simulations are used to explore the interaction with a phospholipid bilayer of the voltage sensor (VS) domain and the S4 helix from the archaebacterial voltage-gated potassium (Kv) channel KvAP. Multiple 2-mus self-assembly simulations reveal that the isolated S4 helix may adopt either interfacial or transmembrane (TM) locations with approximately equal probability. In the TM state, the insertion of the voltage-sensing region of S4 is facilitated via local bilayer deformation that, combined with side chain "snorkeling," enables its Arg side chains to interact with lipid headgroups and water. Multiple 0.2-mus self-assembly simulations of the VS domain are also performed, along with simulations of MscL and KcsA, to permit comparison with more "canonical" integral membrane protein structures. All three stably adopt a TM orientation within a bilayer. For MscL and KcsA, there is no significant bilayer deformation. In contrast, for the VS, there is considerable local deformation, which is again primarily due to the lipid-exposed S4. It is shown that for both the VS and isolated S4 helix, the positively charged side chains of S4 are accommodated within the membrane through a combination of stabilizing interactions with lipid glycerol and headgroup regions, water, and anionic side chains. Our results support the possibility that bilayer deformation around key gating charge residues in Kv channels may result in "focusing" of the electrostatic field, and indicate that, when considering competing models of voltage-sensing, it is essential to consider the dynamics and structure of not only the protein but also of the local lipid environment.

The nature of voltage sensing by voltage-activated ion channels is a key problem in membrane protein structural biology. The way in which the voltage-sensor (VS) domain interacts with its membrane environment remains unclear. In particular, the known structures of Kv channels do not readily explain how a positively charged S4 helix is able to stably span a lipid bilayer. Extended (2 x 50 ns) molecular dynamics simulations of the high-resolution structure of the isolated VS domain from the archaebacterial potassium channel KvAP, embedded in zwitterionic and in anionic lipid bilayers, have been used to explore VS/lipid interactions at atomic resolution. The simulations reveal penetration of water into the center of the VS and bilayer. Furthermore, there is significant local deformation of the lipid bilayer by interactions between lipid phosphate groups and arginine side chains of S4. As a consequence of this, the electrostatic field is "focused" across the center of the bilayer.

VSTx1 is a tarantula venom toxin which binds to the archaebacterial voltage-gated potassium channel KvAP. VSTx1 is thought to access the voltage sensor domain of the channel via the lipid bilayer phase. In order to understand its mode of action and implications for the mechanism of channel activation, it is important to characterize the interactions of VSTx1 with lipid bilayers. Molecular dynamics (MD) simulations (for a total simulation time in excess of 0.2 micros) have been used to explore VSTx1 localization and interactions with zwitterionic (POPC) and with anionic (POPE/POPG) lipid bilayers. In particular, three series of MD simulations have been used to explore the net drift of VSTx1 relative to the center of a bilayer, starting from different locations of the toxin. The preferred location of the toxin is at the membrane/water interface. Although there are differences between POPC and POPE/POPG bilayers, in both cases the toxin forms favorable interactions at the interface, maximizing H-bonding to lipid headgroups and to water molecules while retaining interactions with the hydrophobic core of the bilayer. A 30 ns unrestrained simulation reveals dynamic partitioning of VSTx1 into the interface of a POPC bilayer. The preferential location of VSTx1 at the interface is discussed in the context of Kv channel gating models and provides support for a mode of action in which the toxin interacts with the Kv voltage sensor "paddle" formed by the S3 and S4 helices.

Ion channel-forming peptides enable us to study the conformational dynamics of a transmembrane helix as a function of sequence and environment. Molecular dynamics simulations are used to study the conformation and dynamics of three 22-residue peptides derived from the second transmembrane domain of the glycine receptor (NK4-M2GlyR-p22). Simulations are performed on the peptide in four different environments: trifluoroethanol/water; SDS micelles; DPC micelles; and a DMPC bilayer. A hierarchy of alpha-helix stabilization between the different environments is observed such that TFE/water < micelles < bilayers. Local clustering of trifluoroethanol molecules around the peptide appears to help stabilize an alpha-helical conformation. Single (S22W) and double (S22W,T19R) substitutions at the C-terminus of NK4-M2GlyR-p22 help to stabilize a helical conformation in the micelle and bilayer environments. This correlates with the ability of the W22 and R19 side chains to form H-bonds with the headgroups of lipid or detergent molecules. This study provides a first atomic resolution comparison of the structure and dynamics of NK4-M2GlyR-p22 peptides in membrane and membrane-mimetic environments, paralleling NMR and functional studies of these peptides.

Interactions of lipids are central to the folding and stability of membrane proteins. Coarse-grained molecular dynamics simulations have been used to reveal the mechanisms of self-assembly of protein/membrane and protein/detergent complexes for representatives of two classes of membrane protein, namely, glycophorin (a simple alpha-helical bundle) and OmpA (a beta-barrel). The accuracy of the coarse-grained simulations is established via comparison with the equivalent atomistic simulations of self-assembly of protein/detergent micelles. The simulation of OmpA/bilayer self-assembly reveals how a folded outer membrane protein can be inserted in a bilayer. The glycophorin/bilayer simulation supports the two-state model of membrane folding, in which transmembrane helix insertion precedes dimer self-assembly within a bilayer. The simulations also suggest that a dynamic equilibrium exists between the glycophorin helix monomer and dimer within a bilayer. The simulated glycophorin helix dimer is remarkably close in structure to that revealed by NMR. Thus, coarse-grained methods may help to define mechanisms of membrane protein (re)folding and will prove suitable for simulation of larger scale dynamic rearrangements of biological membranes.

A computational method to calculate the orientation of membrane-associated alpha-helices with respect to a lipid bilayer has been developed. It is based on a previously derived implicit membrane representation, which was parameterized using the structures of 46 alpha-helical membrane proteins. The method is validated by comparison with an independent data set of six transmembrane and nine antimicrobial peptides of known structure and orientation. The minimum energy orientations of the transmembrane helices were found to be in good agreement with tilt and rotation angles known from solid-state NMR experiments. Analysis of the free-energy landscape found two types of minima for transmembrane peptides: i), Surface-bound configurations with the helix long axis parallel to the membrane, and ii), inserted configurations with the helix spanning the membrane in a perpendicular orientation. In all cases the inserted configuration also contained the global energy minimum. Repeating the calculations with a set of solution NMR structures showed that the membrane model correctly distinguishes native transmembrane from nonnative conformers. All antimicrobial peptides investigated were found to orient parallel and bind to the membrane surface, in agreement with experimental data. In all cases insertion into the membrane entailed a significant free-energy penalty. An analysis of the contributions of the individual residue types confirmed that hydrophobic residues are the main driving force behind membrane protein insertion, whereas polar, charged, and aromatic residues were found to be important for the correct orientation of the helix inside the membrane.

Octyl glucoside (OG) is a detergent widely employed in structural and functional studies of membrane proteins. To better understand the nature of protein-OG interactions, molecular dynamics simulations (duration 10 ns) have been used to explore an alpha-helical membrane protein, GlpF, in OG micelles and in DMPC bilayers. Greater conformational drift of the extramembraneous protein loops, from the initial X-ray structure, is seen for the GlpF-OG simulations than for the GlpF-DMPC simulation. The mobility of the transmembrane alpha-helices is approximately 1.3x higher in the GlpF-OG than the GlpF-DMPC simulations. The detergent is seen to form an irregular torus around the protein. The presence of the protein leads to a small perturbation in the behavior of the alkyl chains in the OG micelle, namely an approximately 15% increase in the trans-gauche(-)-gauche(+) transition time. Aromatic side chains (Trp, Tyr) and basic side chains (Arg, Lys) play an important role in both protein-detergent (OG) and protein-lipid (DMPC) interactions.