A solvent model for simulations of peptides in bilayers. II. Membrane-spanning alpha-helices

We describe application of the implicit solvation model (see the first paper of this series), to Monte Carlo simulations of several peptides in bilayer- and water-mimetic environments, and in vacuum. The membrane-bound peptides chosen were transmembrane segments A and B of bacteriorhodopsin, the hydrophobic segment of surfactant lipoprotein, and magainin2. Their conformations in membrane-like media are known from the experiments. Also, molecular dynamics study of surfactant lipoprotein with different explicit solvents has been reported (Kovacs, H., A. E. Mark, J. Johansson, and W. F. van Gunsteren. 1995. J. Mol. Biol. 247:808-822). The principal goal of this work is to compare the results obtained in the framework of our solvation model with available experimental and computational data. The findings could be summarized as follows: 1) structural and energetic properties of studied molecules strongly depend on the solvent; membrane-mimetic media significantly promote formation of alpha-helices capable of traversing the bilayer, whereas a polar environment destabilizes alpha-helical conformation via reduction of solvent-exposed surface area and packing; 2) the structures calculated in a membrane-like environment agree with the experimental ones; 3) noticeable differences in conformation of surfactant lipoprotein assessed via Monte Carlo simulation with implicit solvent (this work) and molecular dynamics in explicit solvent were observed; 4) in vacuo simulations do not correctly reproduce protein-membrane interactions, and hence should be avoided in modeling membrane proteins.     

Conformation of membrane fusion-active 20-residue peptides with or without lipid bilayers. Implication of alpha-helix formation for membrane fusion

Fusion of small unilamellar vesicles of egg phosphatidylcholine can be triggered with synthetic 20-residue peptides. Taking the N-terminal amino acid sequence of HA-2 polypeptide of influenza virus as a guideline, we designed and synthesized several peptides having amphiphilic structures. Among the peptides so far studied, those active to induce membrane fusion took an alpha-helical conformation in the presence of phospholipid bilayers, while a peptide which was unable to induce membrane fusion was in a beta-structure. Mixing of a pair of positively and negatively charged peptides, which had a complementary arrangement of electric charges to each other, resulted in alpha-helix formation at neutral pH, the condition of forming a randomly coiled conformation for each peptide. We concluded that alpha-helix formation was one of the necessary conditions to trigger a process of membrane fusion, at least in the present set of peptides. Characteristic features of these amphiphilic peptides are also described.     

A model for the coupling of alpha-helix and tertiary contact formation

Peptides corresponding to excised alpha-helical segments of natural proteins can spontaneously form helices in solution. However, peptide helices are usually substantially less stable in solution than in the structural context of a folded protein, because of the additional interactions possible between helices in a protein. Such interactions can be thought of as coupling helix formation and tertiary contact formation. The relative energetic contributions of the two processes to the total energy of the folded state of a protein is a matter of current debate. To investigate this balance, an extended helix-coil model (XHC) that incorporates both effects has been constructed. The model treats helix formation with the Lifson-Roig formalism, which describes helix initiation and propagation through cooperative local interactions. The model postulates an additional parameter representing participation of a site in a tertiary contact. In the model, greater helix stability can be achieved through combinations of these short-range and long-range interactions. For instance, stronger tertiary contacts can compensate for helices with little intrinsic stability. By varying the strength of the nonlocal interactions, the model can exhibit behavior consistent with a variety of qualitative models describing the relative importance of secondary and tertiary structure. Moreover, the model is explicit in that it can be used to fit experimental data to individual peptide sequences, providing a means to quantify the two contributions on a common energetic basis.     

Free-energy determinants of alpha-helix insertion into lipid bilayers.

A detailed treatment is provided of the various free-energy terms that contribute to the transfer of a polyalanine alpha-helix from the aqueous phase into lipid bilayers. In agreement with previous work, the hydrophobic effect is found to provide the major driving force for helix insertion. However, an opposing effect of comparable magnitude is also identified and is attributed to the large free-energy penalty associated with the desolvation of peptide hydrogen bonds on transfer to the low dielectric environment of the bilayer. Lipid perturbation effects as well as the entropy loss associated with helix immobilization in the bilayer are also evaluated. Two configurations of a membrane-bound 25mer polyalanine helix were found to be lower in free energy than the isolated helix in the aqueous phase. The first corresponds to the case of vertical insertion, in which a helix terminus protrudes from each side of the bilayer. The second minimum is for the case of horizontal insertion, for which the helix is adsorbed upon the surface of the bilayer. The calculated free-energy minima are found to be in good agreement with recent measurements of related systems. Large free-energy barriers resulting from desolvation of unsatisfied hydrogen-bonding groups at the helix termini are obtained for both insertion processes. The barriers for insertion are significantly reduced if the helix termini are assumed to be "capped" through the formation of hydrogen bonds with polar sidechains. For uncapped helices, our results support recently proposed models in which helices are inserted by first adsorbing on the membrane surface and then having one terminus "swing around" so as to penetrate the bilayer.     

A theoretical model for the association of amphiphilic transmembrane peptides in lipid bilayers

A theoretical model is proposed for the association of trans-bilayer peptides in lipid bilayers. The model is based on a lattice model for the pure lipid bilayer, which accounts accurately for the most important conformational states of the lipids and their mutual interactions and statistics. Within the lattice formulation the bilayer is formed by two independent monolayers, each represented by a triangular lattice, on which sites the lipid chains are arrayed. The peptides are represented by regular objects, with no internal flexibility, and with a projected area on the bilayer plane corresponding to a hexagon with seven lattice sites. In addition, it is assumed that each peptide surface at the interface with the lipid chains is partially hydrophilic, and therefore interacts with the surrounding lipid matrix via selective anisotropic forces. The peptides would therefore assemble in order to shield their hydrophilic residues from the hydrophobic surroundings. The model describes the self-association of peptides in lipid bilayers via lateral and rotational diffusion, anisotropic lipid-peptide interactions, and peptide-peptide interactions involving the peptide hydrophilic regions. The intent of this model study is to analyse the conditions under which the association of trans-bilayer and partially hydrophilic peptides (or their dispersion in the lipid matrix) is lipid-mediated, and to what extent it is induced by direct interactions between the hydrophilic regions of the peptides. The model properties are calculated by a Monte Carlo computer simulation technique within the canonical ensemble. The results from the model study indicate that direct interactions between the hydrophilic regions of the peptides are necessary to induce peptide association in the lipid bilayer in the fluid phase. Furthermore, peptides within each aggregate are oriented in such a way as to shield their hydrophilic regions from the hydrophobic environment. The average number of peptides present in the aggregates formed depends on the degree of mismatch between the peptide hydrophobic length and the lipid bilayer hydrophobic thickness: The lower the degree of mismatch is the higher this number is. Received: 30 December 1996 / Accepted: 9 May 1997     

Molecular dynamics simulations of model trans-membrane peptides in lipid bilayers: a systematic investigation of hydrophobic mismatch

Hydrophobic mismatch, which is the difference between the hydrophobic length of trans-membrane segments of a protein and the hydrophobic width of the surrounding lipid bilayer, is known to play a role in membrane protein function. We have performed molecular dynamics simulations of trans-membrane KALP peptides (sequence: GKK(LA)nLKKA) in phospholipid bilayers to investigate hydrophobic mismatch alleviation mechanisms. By varying systematically the length of the peptide (KALP15, KALP19, KALP23, KALP27, and KALP31) and the lipid hydrophobic length (DLPC, DMPC, and DPPC), a wide range of mismatch conditions were studied. Simulations of durations of 50-200 ns show that under positive mismatch, the system alleviates the mismatch predominantly by tilting the peptide and to a smaller extent by increased lipid ordering in the immediate vicinity of the peptide. Under negative mismatch, alleviation takes place by a combination of local bilayer bending and the snorkeling of the lysine residues of the peptide. Simulations performed at a higher peptide/lipid molar ratio (1:25) reveal slower dynamics of both the peptide and lipid relative to those at a lower peptide/lipid ratio (1:128). The lysine residues have favorable interactions with specific oxygen atoms of the phospholipid headgroups, indicating the preferred localization of these residues at the lipid/water interface.     

Evaluation of a fast implicit solvent model for molecular dynamics simulations.

A solvation term based on the solvent accessible surface area (SASA) is combined with the CHARMM polar hydrogen force field for the efficient simulation of peptides and small proteins in aqueous solution. Only two atomic solvation parameters are used: one is negative for favoring the direct solvation of polar groups and the other positive for taking into account the hydrophobic effect on apolar groups. To approximate the water screening effects on the intrasolute electrostatic interactions, a distance-dependent dielectric function is used and ionic side chains are neutralized. The use of an analytical approximation of the SASA renders the model extremely efficient (i.e., only about 50% slower than in vacuo simulations). The limitations and range of applicability of the SASA model are assessed by simulations of proteins and structured peptides. For the latter, the present study and results reported elsewhere show that with the SASA model it is possible to sample a significant amount of folding/unfolding transitions, which permit the study of the thermodynamics and kinetics of folding at an atomic level of detail.     

Side-chain interactions determine amyloid formation by model polyglutamine peptides in molecular dynamics simulations

The pathological manifestation of nine hereditary neurodegenerative diseases is the presence within the brain of aggregates of disease-specific proteins that contain polyglutamine tracts longer than a critical length. To improve our understanding of the processes by which polyglutamine-containing proteins misfold and aggregate, we have conducted molecular dynamics simulations of the aggregation of model polyglutamine peptides. This work was accomplished by extending the PRIME model to polyglutamine. PRIME is an off-lattice, unbiased, intermediate-resolution protein model based on an amino acid representation of between three and seven united atoms, depending on the residue being modeled. The effects of hydrophobicity on the system are studied by varying the strength of the hydrophobic interaction from 12.5% to 5% of the hydrogen-bonding interaction strength. In our simulations, we observe the spontaneous formation of aggregates and annular structures that are made up of beta-sheets starting from random configurations of random coils. This result was interesting because tubular protofibrils were recently found in experiments on polyglutamine aggregation and because of Perutz's prediction that polyglutamine would form water-filled nanotubes.     

alpha-helix formation: discontinuous molecular dynamics on an intermediate-resolution protein model.

An intermediate-resolution model of small, homogeneous peptides is introduced, and discontinuous molecular dynamics simulation is applied to study secondary structure formation. Physically, each model residue consists of a detailed three-bead backbone and a simplified single-bead side-chain. Excluded volume and hydrogen bond interactions are constructed with discontinuous (i.e., hard-sphere and square-well) potentials. Simulation results show that the backbone motion of the model is limited to realistic regions of Phi-Psi conformational space. Model polyalanine chains undergo a locally cooperative transition to form alpha-helices that are stabilized by backbone hydrogen bonding, while model polyglycine chains tend to adopt nonhelical structures. When side-chain size is increased beyond a critical diameter, steric interactions prevent formation of long alpha-helices. These trends in helicity as a function of residue type have been well documented by experimental, theoretical, and simulation studies and demonstrate the ability of the intermediate-resolution model developed in this work to accurately mimic realistic peptide behavior. The efficient algorithm used permits observation of the complete helix-coil transition within 15 min on a single-processor workstation, suggesting that simulations of very long times are possible with this model.     

Phosphatidylcholesterol bilayers. A model for phospholipid-cholesterol interaction

Aqueous dispersions of monovalent and divalent cation salts of O-(1,2-dipalmitoyl-sn-glycero-3-phosphoryl) cholesterol form multilamellar vesicles as shown by freeze-fracture electron microscopy, by electron micrographs of the negatively stained liposomes, and by swelling curves of liposomes in hypoosmotic medium. Differential scanning calorimetry reveals that aqueous dispersions of divalent metal salts of O-(1,2-dipalmitoyl-sn-glycero-3-phosphoryl)-cholesterol undergo a characteristic thermotropic phase transition with a relatively large cooperative unit (n > 250 for the calcium salt). In contrast, monovalent cation salts of O-(1,2-dipalmitoyl-sn-glycerol-3-phosphoryl)cholesterol do not show a thermotropic phase transition under comparable conditions. The molecular area of O-(1,2-dipalmitoyl-sn-glycero-3-phosphoryl)cholesterol in a monolayer is the same in the presence and absence of Ca²⁺, and is virtually equal to the area of an equimolar mixture of dipalmitoyl phosphatidic acid and cholesterol. To account for the novel state induced by Ca²⁺ on aqueous dispersions of O-(1,2-dipalmitoyl-sn-glycero-3-phosphoryl)cholesterol (i.e., bilayer organization and highly cooperative phase transition), a linear array model is proposed in which Ca²⁺ bridges adjacent arrays of O-(1,2-dipalmitoyl-sn-glycero-3-phosphoryl)cholesterol molecules, thus freezing the acyl chains in their normal state. One of the main corollaries of the model is that the cooperative unit for a thermotropic phase transition is essentially one-dimensional, rather than a two-dimensional matrix. O-(1,2-Dipalmitoyl-sn-glycero-3-phosphoryl)cholesterol is proposed as an orientationally and conformationally restricted analog of glycerophospholipid plus cholesterol in bilayers.     

Molecular dynamics simulations of indolicidin association with model lipid bilayers

Identifying the mechanisms responsible for the interaction of peptides with cell membranes is critical to the design of new antimicrobial peptides and membrane transporters. We report here the results of a computational simulation of the interaction of the 13-residue peptide indolicidin with single-phase lipid bilayers of dioleoylphosphatidylcholine, distearoylphosphatidylcholine, dioleoylphosphatidylglycerol, and distearoylphosphatidylglycerol. Ensemble analysis of the membrane-bound peptide revealed that, in contrast to the extended, linear backbone structure reported for indolicidin in sodium dodecyl sulphate detergent micelles, the peptide adopts a boat-shaped conformation in both phosphatidylglycerol and phosphatidylcholine lipid bilayers, similar to that reported for dodecylphosphocholine micelles. In agreement with fluorescence and NMR experiments, simulations confirmed that the peptide localizes in the membrane interface, with the distance between phosphate headgroups of each leaflet being reduced in the presence of indolicidin. These data, along with a concomitant decrease in lipid order parameters for the upper-tail region, suggest that indolicidin binding results in membrane thinning, consistent with recent in situ atomic force microscopy studies.     

Binding and insertion of alpha-helical anti-microbial peptides in POPC bilayers studied by molecular dynamics simulations

We have performed molecular dynamics simulations of the interactions of two alpha-helical anti-microbial peptides, magainin2 and its synthetic analog of MSI-78, with palmitoyl-oleoyl-phosphatidylcholine (POPC) lipid bilayers. We used various initial positions and orientations of the peptide with respect to the lipid bilayer, including a surface-bound state parallel to the interface, a trans-membrane state, and a partially inserted state. Our 20 ns long simulations show that both magainin2 and MSI-78 are most stable in the lipid environment, with the peptide destabilized to different extents in both aqueous and lipid/water interfacial environments. We found that there are strong specific interactions between the lysine residues of the peptides and the lipid head-group regions. MSI-78, owing to its large number of lysines, shows better binding characteristics and overall stability when compared to magainin2. We also find that both peptides destabilize the bilayer environment, as observed by the increase in lipid tail disorder and the induction of local curvature on the lipid head-groups by the peptides. From all the simulations, we conclude that the hydrogen bonding interactions between the lysines of the peptides and the oxygens of the polar lipid head-groups are the strongest and determine the overall peptide binding characteristics to the lipids.     

A polyalanine-based peptide cannot form a stable transmembrane alpha-helix in fully hydrated phospholipid bilayers

The conformation and amide proton exchangeability of the peptide acetyl-K(2)-A(24)-K(2)-amide (A(24)) and its interaction with phosphatidylcholine bilayers were examined by a variety of physical techniques. When dissolved in or cast from methanol as a dried film, A(24) is predominantly alpha-helical. In aqueous media, however, A(24) exists primarily as a mixture of helical (though not necessarily alpha-helical) and random coiled structures, both of which allow rapid H-D exchange of all amide protons. When incorporated into phospholipids in the absence of water, A(24) also exists primarily as a transmembrane alpha-helix. However, upon hydration of that system, rapid exchange of all amide protons also occurs along with a marked change in the amide I absorption band of the peptide. Also, when dispersed with phosphatidylcholine in aqueous media, the conformation and thermal stability of A(24) are not significantly altered by the presence of the phospholipid or by its gel/liquid-crystalline phase transition. Differential scanning calorimetric and electron spin resonance spectroscopic studies indicate that A(24) has relatively minor effects on the thermodynamic properties of the lipid hydrocarbon chain-melting phase transition, that it does not abolish the lipid pretransition, and that its presence has no significant effect on the orientational order or rates of motion of the phospholipid hydrocarbon chains. We therefore conclude that A(24) has sufficient alpha-helical propensity, but insufficient hydrophobicity, to maintain a stable transmembrane association with phospholipid bilayers in the presence of water. Instead, it exists primarily as a dynamic mixture of helices and other conformers and resides mostly in the aqueous phase where it interacts weakly with the bilayer surface or with the polar/apolar interfacial region of phosphatidylcholine bilayers. Thus, polyalanine-based peptides are not good models for the transmembrane alpha-helical segments of natural membrane proteins.     

Interaction of small peptides with lipid bilayers.

Molecular dynamics simulations of the tripeptide Ala-Phe-Ala-O-tert-butyl interacting with dimyristoylphosphatidylcholine lipid bilayers have been carried out. The lipid and aqueous environments of the peptide, the alkyl chain order, and the lipid and peptide dynamics have been investigated with use of density profiles, radial distribution functions, alkyl chain order parameter profiles, and time correlation functions. It appears that the alkyl chain region accommodates the peptides in the bilayer with minimal perturbation to this region. The peptide dynamics in the bilayer bound form has been compared with that of the free peptide in water. The peptide structure does not vary on the simulation time scale (of the order of hundreds of picoseconds) compared with the solution structure in which a random structure is observed.     

Monte Carlo studies of a model for lipid-gramicidin A bilayers.

This paper presents results of Monte Carlo simulations of a full bilayer of 200 lipid chains and one gramicidin A dimer. Simulations are described for systems with lipid chains of 14, 16, and 18 carbons, respectively. Using accepted potential functions to calculate interactions between all non-hydrogen atoms a Monte Carlo configuration sampling is generated from which order parameter profiles are calculated and specific configurations are displayed. Results are compared with experimental data for lipid-gramicidin bilayers.     

Sequence-dependence of secondary structure formation: conformational studies of host-guest peptides in alpha-helix and beta-structure supporting media

A newly designed host–guest approach is introduced as a experimental tool to explore the relationship between the sequence of peptides and their secondary structure. From the CD spectra of the host–guest peptides studied, a tentative scale for the α-helix potential in 2,2,2-trifluorethanol of guest amino acids is delineated. The conformational preferences are also examined in β-structure supporting media (solid state, CH₂Cl₂, CH₃OH, H₂O) using ir-absorption and CD techniques. Scales for the β-forming tendency of guest amino acid residues in the different media are delineated. It is shown that the preferred conformation of the host–guest peptides is a function of the medium, the chain length, and the protecting groups. Given the fact that conformational effects are important in peptide synthesis, the tentative scales may serve as a guideline to predict secondary structures of side-chain-protected or -deprotected peptides in a given solvent, complementing the well-known empirical conformational prediction parameters.     

Interactions of hydrophobic peptides with lipid bilayers: Monte Carlo simulations with M2delta

We introduce here a novel Monte Carlo simulation method for studying the interactions of hydrophobic peptides with lipid membranes. Each of the peptide's amino acids is represented as two interaction sites: one corresponding to the backbone alpha-carbon and the other to the side chain, with the membrane represented as a hydrophobic profile. Peptide conformations and locations in the membrane and changes in the membrane width are sampled using the Metropolis criterion, taking into account the underlying energetics. Using this method we investigate the interactions between the hydrophobic peptide M2delta and a model membrane. The simulations show that starting from an extended conformation in the aqueous phase, the peptide first adsorbs onto the membrane surface, while acquiring an ordered helical structure. This is followed by formation of a helical-hairpin and insertion into the membrane. The observed path is in agreement with contemporary understanding of peptide insertion into biological membranes. Two stable orientations of membrane-associated M2delta were obtained: transmembrane (TM) and surface, and the value of the water-to-membrane transfer free energy of each of them is in agreement with calculations and measurements on similar cases. M2delta is most stable in the TM orientation, where it assumes a helical conformation with a tilt of 14 degrees between the helix principal axis and the membrane normal. The peptide conformation agrees well with the experimental data; average root-mean-square deviations of 2.1 A compared to nuclear magnetic resonance structures obtained in detergent micelles and supported lipid bilayers. The average orientation of the peptide in the membrane in the most stable configurations reported here, and in particular the value of the tilt angle, are in excellent agreement with the ones calculated using the continuum-solvent model and the ones observed in the nuclear magnetic resonance studies. This suggests that the method may be used to predict the three-dimensional structure of TM peptides.     

Domain formation in phosphatidylcholine bilayers containing transmembrane peptides: specific effects of flanking residues

Lateral segregation in biological membranes leads to the formation of domains. We have studied the lateral segregation in gel-state model membranes consisting of supported dipalmitoylphosphatidylcholine (DPPC) bilayers with various model peptides, using atomic force microscopy (AFM). The model peptides are derivatives of the Ac-GWWL(AL)(n)WWA-Etn peptides (the so-called WALP peptides) and have instead of tryptophans, other flanking residues. In a previous study, we found that WALP peptides induce the formation of extremely ordered, striated domains in supported DPPC bilayers. In this study, we show that WALP analogues with other uncharged residues (tyrosine, phenylalanine, or histidine at pH 9) can also induce the formation of striated domains, albeit in some cases with a slightly different pattern. The WALP analogues with positively charged residues (lysine or histidine at low pH) cannot induce striated domains and give rise to a completely different morphology: they induce irregularly shaped depressions in DPPC bilayers. The latter morphology is explained by the fact that the positively charged peptides repel each other and hence are not able to form striated domains in which they would have to be in close vicinity. They would reside in disordered, fluidlike lipid areas, appearing below the level of the ordered gel-state lipid domains, which would account for the irregularly shaped depressions.