A Steroid in a Lipid Bilayer: Localization, Orientation, and Energetics

Steroid hormones are known to freely partition into lipid bilayers. As a case study, we investigated the behavior of the steroid hormone cortisone in a model lipid bilayer. First, we looked at energy barriers involved in the partitioning of a single molecule into a bilayer using umbrella sampling molecular dynamics simulations. A rather wide well of −4.5 kcal/mol was observed in the interfacial region between the lipid headgroup and tailgroup. Next, using two unconstrained molecular dynamics simulations with cortisone initially positioned at distinct locations within a bilayer, we studied the preferred location and orientation of the molecule. Finally, we observed how cortisone molecules could spontaneously insert and localize in a bilayer from bulk solution. The three independent approaches produced a converged picture of how cortisone behaves in a model lipid bilayer.     

The effect of cholesterol on short- and long-chain monounsaturated lipid bilayers as determined by molecular dynamics simulations and X-ray scattering

We investigate the structure of cholesterol-containing membranes composed of either short-chain (diC14:1PC) or long-chain (diC22:1PC) monounsaturated phospholipids. Bilayer structural information is derived from all-atom molecular dynamics simulations, which are validated via direct comparison to x-ray scattering experiments. We show that the addition of 40 mol % cholesterol results in a nearly identical increase in the thickness of the two different bilayers. In both cases, the chain ordering dominates over the hydrophobic matching between the length of the cholesterol molecule and the hydrocarbon thickness of the bilayer, which one would expect to cause a thinning of the diC22:1PC bilayer. For both bilayers there is substantial headgroup rearrangement for lipids directly in contact with cholesterol, supporting the so-called umbrella model. Importantly, in diC14:1PC bilayers, a dynamic network of hydrogen bonds stabilizes long-lived reorientations of some cholesterol molecules, during which they are found to lie perpendicular to the bilayer normal, deep within the bilayer’s hydrophobic core. Additionally, the simulations show that the diC14:1PC bilayer is significantly more permeable to water. These differences may be correlated with faster cholesterol flip-flop between the leaflets of short-chain lipid bilayers, resulting in an asymmetric distribution of cholesterol molecules. This asymmetry was observed experimentally in a case of unilamellar vesicles (ULVs), and reproduced through a set of novel asymmetric simulations. In contrast to ULVs, experimental data for oriented multilamellar stacks does not show the asymmetry, suggesting that it results from the curvature of the ULV bilayers.     

Vstx1, a modifier of Kv channel gating, localizes to the interfacial region of lipid bilayers

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.     

Effect of cholesterol on the structure and dynamic properties of unsaturated phospholipid bilayers

Molecular dynamics (MD) computer simulations of five different hydrated unsaturated phosphatidylcholine lipid bilayers built up by 18:0/18:1(n-9)cis PC, 18:0/18:2(n-6)cis PC, 18:0/18:3(n-3)cis PC, 18:0/20:4(n-6)cis PC, and 18:0/22:6(n-3)cis PC molecules with 40 mol% cholesterol, and the same five pure phosphatidylcholine bilayers have been performed at 303 K. The simulation box of a lipid bilayer contained 96 phosphatidylcholines, 64 cholesterols, and 3840 water molecules (48 phosphatidylcholine molecules and 32 cholesterols per layer and 24 water molecules per phospholipid or cholesterol in each case). The lateral self-diffusion coefficients of the lipids in these systems and mass density profiles with respect to the bilayer normal have been analyzed. It has been found that the lateral diffusion coefficients of phosphatidylcholine molecules increase with increasing number of double bonds in one of the lipid chains, both in pure bilayers and in bilayers with cholesterol. It has been found as well that the lateral diffusion coefficient of phosphatidylcholine molecules of a lipid bilayer with 40 mol% cholesterol is smaller than that for the corresponding pure phosphatidylcholine bilayer.     

The influence of cholesterol on interactions and dynamics of ibuprofen in a lipid bilayer

In this work, molecular dynamics (MD) simulations with atomistic details were performed to examine the influence of the cholesterol on the interactions and the partitioning of the hydrophobic drug ibuprofen in a fully hydrated 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) bilayer. Analysis of MD simulations indicated that ibuprofen molecules prefer to be located in the hydrophobic acyl chain region of DMPC/cholesterol bilayers. This distribution decreases the lateral motion of lipid molecules. The presence of ibuprofen molecules in the bilayers with 0 and 25 mol% cholesterol increases the ordering of hydrocarbon tails of lipids whereas for the bilayers with 50 mol% cholesterol, ibuprofen molecules perturb the flexible chains of DMPC lipids which leads to the reduction of the acyl chain order parameter. The potential of the mean force (PMF) method was used to calculate the free energy profile for the transferring of an ibuprofen molecule from the bulk water into the DMPC/cholesterol membranes. The PMF studies indicated that the presence of 50 mol% cholesterol in the bilayers increases the free energy barrier and slows down the permeation of the ibuprofen drug across the DMPC bilayer. This can be due to the condensing and ordering effects of the cholesterol on the bilayer.     

Properties of unsaturated phospholipid bilayers: Effect of cholesterol

Properties of hydrated unsaturated phosphatidylcholine (PC) lipid bilayers containing 40 mol % cholesterol and of pure PC bilayers have been studied. Various methods were applied, including molecular dynamics simulations, self-consistent field calculations, and the pulsed field gradient nuclear magnetic resonance technique. Lipid bilayers were composed of 18:0/18:1(n-9)cis PC, 18:0/18:2(n-6)cis PC, 18:0/18:3(n-3)cis PC, 18:0/20:4(n-6)cis PC, and 18:0/22:6(n-3)cis PC molecules. Lateral self-diffusion coefficients of the lipids in all these bilayers, mass density distributions of atoms and atom groups with respect to the bilayer normal, the C-H and C-C bond order parameter profiles of each phospholipid hydrocarbon chain with respect to the bilayer normal were calculated. It was shown that the lateral self-diffusion coefficient of PC molecules of the lipid bilayer containing 40 mol % cholesterol is smaller than that for a corresponding pure PC bilayer; the diffusion coefficients increase with increasing the degree of unsaturation of one of the PC chains in bilayers of both types (i.e., in pure bilayers or in bilayers with cholesterol). The presence of cholesterol in a bilayer promoted the extension of saturated and polyunsaturated lipid chains. The condensing effect of cholesterol on the order parameters was more pronounced for the double C=C bonds of polyunsaturated chains than for single C-C bonds of saturated chains.     

Molecular dynamics simulations and 2H NMR study of the GalCer/DPPG lipid bilayer

Molecular dynamics simulations were performed on a two-component lipid bilayer system in the liquid crystalline phase at constant pressure and constant temperature. The lipid bilayers were composed of a mixture of neutral galactosylceramide (GalCer) and charged dipalmitoylphosphatidylglycerol (DPPG) lipid molecules. Two lipid bilayer systems were prepared with GalCer:DPPG ratio 9:1 (10%-DPPG system) and 3:1 (25%-DPPG system). The 10%-DPPG system represents a collapsed state lipid bilayer, with a narrow water space between the bilayers, and the 25%-DPPG system represents an expanded state with a fluid space of approximately 10 nm. The number of lipid molecules used in each simulation was 1024, and the length of the production run simulation was 10 ns. The simulations were validated by comparing the results with experimental data for several important aspects of the bilayer structure and dynamics. Deuterium order parameters obtained from (2)H NMR experiments for DPPG chains are in a very good agreement with those obtained from molecular dynamics simulations. The surface area per GalCer lipid molecule was estimated to be 0.608 +/- 0.011 nm(2). From the simulated electron density profiles, the bilayer thickness defined as the distance between the phosphorus peaks across the bilayer was calculated to be 4.21 nm. Both simulation systems revealed a tendency for cooperative bilayer undulations, as expected in the liquid crystalline phase. The interaction of water with the GalCer and DPPG oxygen atoms results in a strong water ordering in a spherical hydration shell and the formation of hydrogen bonds (H-bonds). Each GalCer lipid molecule makes 8.6 +/- 0.1 H-bonds with the surrounding water, whereas each DPPG lipid molecule makes 8.3 +/- 0.1 H-bonds. The number of water molecules per GalCer or DPPG in the hydration shell was estimated to be 10-11 from an analysis of the radial distribution functions. The formation of the intermolecular hydrogen bonds was observed between hydroxyl groups from the opposing GalCer sugar headgroups, giving an energy of adhesion in the range between -1.0 and -3.4 erg/cm(2). We suggest that this value is the contribution of the hydrogen-bond component to the net adhesion energy between GalCer bilayers in the liquid crystalline phase.     

Association of ethanol with lipid membranes containing cholesterol, sphingomyelin and ganglioside: a titration calorimetry study.

The association of ethanol at physiologically relevant concentrations with lipid bilayers of different lipid composition has been investigated by use of isothermal titration calorimetry (ITC). The liposomes examined were composed of combinations of lipids commonly found in neural cell membranes: dimyristoyl phosphatidylcholine (DMPC), ganglioside (GM(1)), sphingomyelin and cholesterol. The calorimetric results show that the interaction of ethanol with fluid lipid bilayers is endothermic and strongly dependent on the lipid composition of the liposomes. The data have been used to estimate partitioning coefficients for ethanol into the fluid lipid bilayer phase and the results are discussed in terms of the thermodynamics of partitioning. The presence of 10 mol% sphingomyelin or ganglioside in DMPC liposomes enhances the partitioning coefficient by a factor of 3. Correspondingly, cholesterol (30 mol%) reduces the partitioning coefficient by a factor of 3. This connection between lipid composition and partitioning coefficient correlates with in vivo observations. Comparison of the data with the molecular structure of the lipid molecules suggests that ethanol partitioning is highly sensitive to changes in the lipid backbone (glycerol or ceramide) while it appears much less sensitive to the nature of the head group.     

Arginine in Membranes: The Connection Between Molecular Dynamics Simulations and Translocon-Mediated Insertion Experiments

Several laboratories have carried out molecular dynamics (MD) simulations of arginine interactions with lipid bilayers and found that the energetic cost of placing arginine in lipid bilayers is an order of magnitude greater than observed in molecular biology experiments in which Arg-containing transmembrane helices are inserted across the endoplasmic reticulum membrane by the Sec61 translocon. We attempt here to reconcile the results of the two approaches. We first present MD simulations of guanidinium groups alone in lipid bilayers, and then, to mimic the molecular biology experiments, we present simulations of hydrophobic helices containing single Arg residues at different positions along the helix. We discuss the simulation results in the context of molecular biology results and show that the energetic discrepancy is reduced, but not eliminated, by considering free energy differences between Arg at the interface and at the center of the model helices. The reduction occurs because Arg snorkeling to the interface prevents Arg from residing in the bilayer center where the energetic cost of desolvation is highest. We then show that the problem with MD simulations is that they measure water-to-bilayer free energies, whereas the molecular biology experiments measure the energetics of partitioning from translocon to bilayer, which raises the fundamental question of the relationship between water-to-bilayer and water-to-translocon partitioning. We present two thermodynamic scenarios as a foundation for reconciliation of the simulation and molecular biology results. The simplest scenario is that translocon-to-bilayer partitioning is independent of water-to-bilayer partitioning; there is no thermodynamic cycle connecting the two paths.     

Molecular dynamics simulations of charged and neutral lipid bilayers: treatment of electrostatic interactions

Molecular dynamics (MD) simulations complement experimental methods in studies of the structure and dynamics of lipid bilayers. The choice of algorithms employed in this computational method represents a trade-off between the accuracy and real calculation time. The largest portion of the simulation time is devoted to calculation of long-range electrostatic interactions. To speed-up evaluation of these interactions, various approximations have been used. The most common ones are the truncation of long-range interactions with the use of cut-offs, and the particle-mesh Ewald (PME) method. In this study, several multi-nanosecond cut-off and PME simulations were performed to establish the influence of the simulation protocol on the bilayer properties. Two bilayers were used. One consisted of neutral phosphatidylcholine molecules. The other was a mixed lipid bilayer consisting of neutral phosphatidylethanolamine and negatively charged phosphatidylglycerol molecules. The study shows that the cut-off simulation of a bilayer containing charge molecules generates artefacts; in particular the mobility and order of the charged molecules are vastly different from those determined experimentally. In the PME simulation, the bilayer properties are in general agreement with experimental data. The cut-off simulation of bilayers containing only uncharged molecules does not generate artefacts, nevertheless, the PME simulation gives generally better agreement with experimental data.     

NMR study of the interaction of retinoids with phospholipid bilayers

The interaction of three vitamin A derivatives or retinoids: all-trans-retinoic acid, 13-cis-retinoic acid and retinol with multilamellar phospholipid bilayers was studied using a combination of 2H- and 31P-NMR measurements. The following model membrane systems were used: (1) dipalmitoylphosphatidylcholine (DPPC) bilayers; (2) bilayers composed of a mixture of DPPC and bovine heart phosphatidylcholine (PC); (3) mixed PC/phosphatidylethanolamine (PE) bilayers. Only a weak interaction was observed between 13-cis-retinoic acid and DPPC membranes. Addition of all-trans-retinoic acid at a molar ratio of 1:2 to the lipid causes a small decrease (5 C degrees) in the gel to liquid crystalline phase-transition temperature of DPPC, a small increase in the order parameters of the lipid side-chains of single component bilayers and no measurable effect in the other lipid systems studied. Considerably larger perturbation in the lipid bilayer structure is introduced by addition of retinol which, at a molar ratio of 1:2 to the lipid, lowered the gel to liquid crystalline phase-transition temperature of DPPC by 21 C degrees and caused a decrease of order parameters of the lipid side-chains in all three lipid bilayer systems. These effects are consistent with intercalation of retinol molecules into the bilayer interior. The results for the mixed PC/PE bilayers indicate that the presence of retinol caused lateral separation of PE- and retinol-enriched regions.     

Distribution of halothane in a dipalmitoylphosphatidylcholine bilayer from molecular dynamics calculations

We report a 2-ns constant pressure molecular dynamics simulation of halothane, at a mol fraction of 50%, in the hydrated liquid crystal bilayer phase of dipalmitoylphosphatidylcholine. Halothane molecules are found to preferentially segregate to the upper part of the lipid acyl chains, with a maximum probability near the C(5) methylene groups. However, a finite probability is also observed along the tail region and across the methyl trough. Over 95% of the halothane molecules are located below the lipid carbonyl carbons, in agreement with photolabeling experiments. Halothane induces lateral expansion and a concomitant contraction in the bilayer thickness. A decrease in the acyl chain segment order parameters, S(CD), for the tail portion, and a slight increase for the upper portion compared to neat bilayers, are in agreement with several NMR studies on related systems. The decrease in S(CD) is attributed to a larger accessible volume per lipid in the tail region. Significant changes in the electric properties of the lipid bilayer result from the structural changes, which include a shift and broadening of the choline headgroup dipole (P-N) orientation distribution. Our findings reconcile apparent controversial conclusions from experiments on diverse lipid systems.     

Molecular packing in 1-hexanol-DMPC bilayers studied by molecular dynamics simulation

The structure and molecular packing density of a "mismatched" solute, 1-hexanol, in lipid membranes of dimyristoyl phosphatidylcholine (DMPC) was studied by molecular dynamics simulations. We found that the average location and orientation of the hexanol molecules matched earlier experimental data on comparable systems. The local density or molecular packing in DMPC-hexanol was elucidated through the average Voronoi volumes of all heavy (non-hydrogen) atoms. Analogous analysis was conducted on trajectories from simulations of pure 1-hexanol and pure (hydrated) DMPC bilayers. The results suggested a positive volume change, DeltaV(m), of 4 cm(3) mol(-1) hexanol partitioned at 310 K in good accordance with experimental values. Analysis of the apparent volumes of each component in the pure and mixed states further showed that DeltaV(m) reflects a balance between a substantial increase in the packing density of the alcohol upon partitioning and an even stronger loosening in the packing of the lipid. Furthermore, analysis of Voronoi volumes along the membrane normal identifies a distinctive depth dependence of the changes in molecular packing. The outer (interfacial) part of the lipid acyl chains (up to C8) is stretched by about 4%. Concomitantly, the average lateral area per chain decreases and these two effects compensate so that the overall packing density in the outer region, where the hexanol molecules are located, remains practically constant. The core of the bilayer (C9-C13) is slightly thinned. The average lateral area per chain in this region expands, resulting in a looser packing density. The net effect in the core is a 2-3% decrease in density corresponding to a total volume increase of approximately 14 cm(3) mol(-1) hexanol partitioned.     

Partitioning of 2,6-Bis(1H-Benzimidazol-2-yl)pyridine fluorophore into a phospholipid bilayer: complementary use of fluorescence quenching studies and molecular dynamics simulations

Successful use of fluorescence sensing in elucidating the biophysical properties of lipid membranes requires knowledge of the distribution and location of an emitting molecule in the bilayer. We report here that 2,6-bis(1H-benzimidazol-2-yl)pyridine (BBP), which is almost non-fluorescent in aqueous solutions, reveals a strong emission enhancement in a hydrophobic environment of a phospholipid bilayer, making it interesting for fluorescence probing of water content in a lipid membrane. Comparing the fluorescence behavior of BBP in a wide variety of solvents with those in phospholipid vesicles, we suggest that the hydrogen bonding interactions between a BBP fluorophore and water molecules play a crucial role in the observed “light switch effect”. Therefore, the loss of water-induced fluorescence quenching inside a membrane are thought to be due to deep penetration of BBP into the hydrophobic, water-free region of a bilayer. Characterized by strong quenching by transition metal ions in solution, BBP also demonstrated significant shielding from the action of the quencher in the presence of phospholipid vesicles. We used the increase in fluorescence intensity, measured upon titration of probe molecules with lipid vesicles, to estimate the partition constant and the Gibbs free energy (ΔG) of transfer of BBP from aqueous buffer into a membrane. Partitioning BBP revealed strongly favorable ΔG, which depends only slightly on the lipid composition of a bilayer, varying in a range from − 6.5 to − 7.0 kcal/mol. To elucidate the binding interactions of the probe with a membrane on the molecular level, a distribution and favorable location of BBP in a POPC bilayer were modeled via atomistic molecular dynamics (MD) simulations using two different approaches: (i) free, diffusion-driven partitioning of the probe molecules into a bilayer and (ii) constrained umbrella sampling of a penetration profile of the dye molecule across a bilayer. Both of these MD approaches agreed with regard to the preferred location of a BBP fluorophore within the interfacial region of a bilayer, located between the hydrocarbon acyl tails and the initial portion of the lipid headgroups. MD simulations also revealed restricted permeability of water molecules into this region of a POPC bilayer, determining the strong fluorescence enhancement observed experimentally for the membrane-partitioned form of BBP.     

Membrane Partitioning of Anionic,Ligand-Coated Nanoparticles Is Accompanied by Ligand Snorkeling,Local Disordering,and Cholesterol Depletion

Intracellular uptake of nanoparticles (NPs) may induce phase transitions, restructuring, stretching, or even complete disruption of the cell membrane. Therefore, NP cytotoxicity assessment requires a thorough understanding of the mechanisms by which these engineered nanostructures interact with the cell membrane. In this study, extensive Coarse-Grained Molecular Dynamics (MD) simulations are performed to investigate the partitioning of an anionic, ligand-decorated NP in model membranes containing dipalmitoylphosphatidylcholine (DPPC) phospholipids and different concentrations of cholesterol. Spontaneous fusion and translocation of the anionic NP is not observed in any of the 10-µs unbiased MD simulations, indicating that longer timescales may be required for such phenomena to occur. This picture is supported by the free energy analysis, revealing a considerable free energy barrier for NP translocation across the lipid bilayer. 5-µs unbiased MD simulations with the NP inserted in the bilayer core reveal that the hydrophobic and hydrophilic ligands of the NP surface rearrange to form optimal contacts with the lipid bilayer, leading to the so-called snorkeling effect. Inside cholesterol-containing bilayers, the NP induces rearrangement of the structure of the lipid bilayer in its vicinity from the liquid-ordered to the liquid phase spanning a distance almost twice its core radius (8–10 nm). Based on the physical insights obtained in this study, we propose a mechanism of cellular anionic NPpartitioning, which requires structural rearrangements of both the NP and the bilayer, and conclude that the translocation of anionic NPs through cholesterol-rich membranes must be accompanied by formation of cholesterol-lean regions in the proximity of NPs.     

Triglyceride Blisters in Lipid Bilayers: Implications for Lipid Droplet Biogenesis and the Mobile Lipid Signal in Cancer Cell Membranes

Triglycerides have a limited solubility, around 3%, in phosphatidylcholine lipid bilayers. Using millisecond-scale course grained molecular dynamics simulations, we show that the model lipid bilayer can accommodate a higher concentration of triolein (TO) than earlier anticipated, by sequestering triolein molecules to the bilayer center in the form of a disordered, isotropic, mobile neutral lipid aggregate, at least 17 nm in diameter, which forms spontaneously, and remains stable on at least the microsecond time scale. The results give credence to the hotly debated existence of mobile neutral lipid aggregates of unknown function present in malignant cells, and to the early biogenesis of lipid droplets accommodated between the two leaflets of the endoplasmic reticulum membrane. The TO aggregates give the bilayer a blister-like appearance, and will hinder the formation of multi-lamellar phases in model, and possibly living membranes. The blisters will result in anomalous membrane probe partitioning, which should be accounted for in the interpretation of probe-related measurements.     

The effect of benzyl alcohol and cholesterol on the acyl chain order and alkane solubility of bimolecular phosphatidylcholine membranes

An investigation was made of the effects of cholesterol and benzyl alcohol on the partitioning of n-alkanes between lipid bilayer membranes and bulk lipid/alkane solutions (in the torus). Bilayers were formed from solutions containing alkanes of different chain lengths, together with phosphatidylcholine and cholesterol in varying proportions. The partitioning of the alkanes was determined from measurements of the very low frequency (1 Hz) capacitance of the membranes. Perturbation of the internal membrane structure by the inclusion of cholesterol and benzyl alcohol produced very significant changes in the n-alkane partition coefficient, cholesterol causing a decrease and benzyl alcohol an increase in the alkane partitioning into the bilayer. A correlation exists between the effects of these compounds on the alkane partitioning and their effect on the segmental chain order of the acyl chains in the bilayer and this correlation is consistent with a statistical-mechanical model of the lipid/alkane bilayers in the liquid crystalline state. The perturbation by cholesterol and benzyl alcohol of the internal structure of membranes bears on the conflicting reports of the effects of these substances on artificial lipid bilayers and could also be relevant to their known physiological effects.     

Interaction Between Amyloid-β (1-42) Peptide and Phospholipid Bilayers: A Molecular Dynamics Study

The amyloid-β (Aβ) peptide is a key aggregate species in Alzheimer's disease. Although important aspects of Aβ peptide aggregation are understood, the initial stage of aggregation from monomer to oligomer is still not clear. One potential mediator of this early aggregation process is interactions of Aβ with anionic cell membranes. We used unconstrained and umbrella sampling molecular dynamics simulations to investigate interactions between the 42-amino acid Aβ peptide and model bilayers of zwitterionic dipalmitoylphosphatidylcholine (DPPC) lipids and anionic dioleoylphosphatidylserine (DOPS) lipids. Using these methods, we determined that Aβ is attracted to the surface of DPPC and DOPS bilayers over the small length scales used in these simulations. We also found supporting evidence that the charge on both the bilayer surface and the peptide affects the free energy of binding of the peptide to the bilayer surface and the distribution of the peptide on the bilayer surface. Our work demonstrates that interactions between the Aβ peptide and lipid bilayer promotes a peptide distribution on the bilayer surface that is prone to peptide-peptide interactions, which can influence the propensity of Aβ to aggregate into higher-order structures.     

Implicit solvent model estimates of the stability of model structures of the alamethicin channel

Alamethicin is a hydrophobic helical peptide of 20 residues, which oligomerizes to form ion-conducting channels in membranes. The behavior of an intact alamethicin channel in POPC bilayers was recently studied, using 2 ns molecular dynamics (MD) simulations of a model hexameric channel. These simulations produced numerous conformations of the channel. In the present study, we used 11 of these channel conformations and carried out continuum-solvent model calculations, similar to those used for the monomers in our previous studies, to investigate the energetics of the channel inside the lipid bilayer. Our results suggest that, out of the 11 channel conformations produced by the MD simulations, only four are stable inside the lipid bilayer, with water-to-membrane free energies of transfer ranging from ~–6 to ~–10 kcal/mol. Analysis of the results suggests two causes for the apparent instability of the remainder of the structures inside the lipid bilayer, both resulting from the desolvation of channel polar groups (i.e. their transfer from the aqueous phase into the bilayer). The first is specific, uncompensated backbone hydrogen bonds, which exist in the region of the channel exposed to the hydrocarbon of the lipid bilayer. The second is exposure of intra-pore water molecules to the surrounding lipid. Thus, the association of these structures with the membrane involves a large electrostatic desolvation free-energy penalty. The apparent conflict between continuum-solvent and MD calculations, and its significance for the interpretation of membrane proteins simulations, are discussed.