Article Information

• PDF (631 kb)
• Full Text with Large Figures
• Cited by in Scopus (4)

PubMed

• Articles by Santi Esteban-Martín
• Articles by Jesús Salgado

Related Articles

• Lipid Membrane Expansion and Micelle Formation by Polymer-Gr...

  Lipid Membrane Expansion and Micelle Formation by Polymer-Grafted Lipids: Scaling with Polymer Length Studied by Spin-Label Electron Spin Resonance Biophysical Journal, Volume 80, Issue 3, 1 March 2001, Pages 1372-1383 Giuseppina Montesano, Rosa Bartucci, Salvatore Belsito, Derek Marsh and Luigi Sportelli Abstract Spin-label electron spin resonance (ESR) spectroscopy and auxiliary optical density measurements are used to study lipid dispersions of -poly(ethylene glycol)-dipalmitoyl phosphatidylethanolamine (PEG:5000-DPPE) mixed with dipalmitoyl phosphatidylcholine (DPPC). PEG:5000-DPPE bears a large hydrophilic polymer headgroup (with ∼114 oxyethylene monomers) and is commonly used for steric stabilization of liposomes used in drug delivery. Comparison is made with results from mixtures of DPPC with polymer lipids bearing shorter headgroups (∼45 and 8 oxyethylene monomers). ESR spectra of phosphatidylcholine spin-labeled on the 5-C atom position of the -2 chain are shown to reflect the area expansion of the lipid membranes by the lateral pressure exerted in the polymer brush, in a way that is consistent with theory. The lipid chain packing density at the onset of micelle formation is the same for all three PEG-lipids, although the mole fraction at which this occurs differs greatly. The mole fraction at onset scales inversely with the size of the polymer headgroup, where the experimental exponent of 0.7 is close to theoretical predictions (viz. 0.55–0.6). The mole fraction of PEG-lipid at completion of micelle formation is more weakly dependent on polymer size, which conforms with theoretical predictions. At high mole fractions of PEG:5000-DPPE the dependence of lipid packing density on mole fraction is multiphasic, which differs qualitatively from the monotonic decrease in packing density found with the shorter polymer lipids. Lipid spin-label ESR is an experimental tool that complements theoretical analysis using polymer models combined with the lipid equation of state. Abstract | Full Text | PDF (242 kb)

• Origin of Laurdan Sensitivity to the Vesicle-to-Micelle Tran...

  Origin of Laurdan Sensitivity to the Vesicle-to-Micelle Transition of Phospholipid-Octylglucoside System: A Time-Resolved Fluorescence Study Biophysical Journal, Volume 80, Issue 1, 1 January 2001, Pages 347-359 Mathias Viard, Jacques Gallay, Michel Vincent and Maïté Paternostre Abstract The fluorescent probe laurdan has been shown to be sensitive to the vesicle-to-micelle transition of phosphatidylcholine/octylglucoside (M. Paternostre, O. Meyer, C. Grabielle-Madelmont, S. Lesieur, and M. Ollivon, 1995, 69:2476–2488). On the other hand, a study on the photophysics of laurdan in organic solvents has shown that the complex de-excitation pathway of the probe can be described by two successive processes, i.e., an intramolecular charge transfer followed by dielectric relaxation of the solvent if polar. These two excited-state reactions lead to three emitting states, i.e., a locally excited state, a charge transfer state, and a solvent relaxed state (M. Viard, J. Gallay, M. Vincent, B. Robert and M. Paternostre, 1997, 73:2221–2234). Experiments have been performed using time-resolved fluorescence on the probe inserted in amphiphile aggregates (mixed liposomes, mixed micelles) different in detergent-to-lipid ratios. The results have been compared with those obtained for laurdan inserted in dipalmitoyl phosphatidylcholine liposomes in the gel and in the fluid lamellar phase. Except for laurdan in dipalmitoyl phosphatidylcholine liposomes in the gel lamellar phase, the red part of the emission spectra originates from the de-excitation of the relaxed excited state of laurdan, indicating that indeed the dielectric relaxation process is an important phenomena in the ground-state return pathway of this probe. On the other hand, the maximization entropy method (MEM) analysis of the fluorescence decay recorded in the blue part of the emission spectra indicates that the dielectric relaxation is not the only reaction occurring to the excited state of laurdan. Moreover, the analysis of the fluorescence decays of laurdan inserted in gel lamellar dipalmitoylphosphatidylcholine (DPPC) liposomes indicates excited-state reactions, although dielectric relaxation is impossible. These results are in agreement with the de-excitation pathway determined from laurdan behavior in organic solvent even if, in most of the aggregates studied in this work, the major phenomenon is the dielectric relaxation of the solvent. All along the vesicle-to-micelle transition, we have observed that the lifetime of the relaxed excited state of laurdan continuously decreases probably due to a dynamic quenching process by water molecules. On the other hand, the time constant of the dielectric relaxation process remains almost unchanged in the lamellar part of the transition but abruptly decreases as soon as the first mixed micelle is formed. This decrease is continuous all over the rest of the transition even if it is more pronounced in the mixed liposomes’ and mixed micelles’ coexistence. The increase of the octylglucoside-to-lipid ratio of the mixed micelles via the change of the size and the shape of the aggregates may facilitate the penetration and the mobility of water molecules. Therefore, during the vesicle-to-micelle transition, laurdan probes the evolution of both the amphiphile packing in the aggregates and the increase of the interface polarity. This study finally shows that the detergent-to-lipid ratio of the mixed micelles is an important parameter to control to limit the penetration and the mobility of water within the amphiphile aggregates and that laurdan is a nice tool to monitor this phenomenon. Abstract | Full Text | PDF (317 kb)

• Liposomes, Disks, and Spherical Micelles: Aggregate Structur...

  Liposomes, Disks, and Spherical Micelles: Aggregate Structure in Mixtures of Gel Phase Phosphatidylcholines and Poly(Ethylene Glycol)-Phospholipids Biophysical Journal, Volume 85, Issue 6, 1 December 2003, Pages 3839-3847 Markus Johnsson and Katarina Edwards Abstract Poly(ethylene glycol) (PEG) decorated lipid bilayers are widely used in biomembrane and pharmaceutical research. The success of PEG-lipid stabilized liposomes in drug delivery is one of the key factors for the interest in these polymer/lipid systems. From a more fundamental point of view, it is essential to understand the effect of the surface grafted polymers on the physical-chemical properties of the lipid bilayer. Herein we have used cryo-transmission electron microscopy and dynamic light scattering to characterize the aggregate structure and phase behavior of mixtures of PEG-lipids and distearoylphosphatidylcholine or dipalmitoylphosphatidylcholine. The PEG-lipids contain PEG of molecular weight 2000 or 5000. We show that the transition from a dispersed lamellar phase (liposomes) to a micellar phase consisting of small spherical micelles occurs via the formation of small discoidal micelles. The onset of disk formation already takes place at low PEG-lipid concentrations (<5mol %) and the size of the disks decreases as more PEG-lipid is added to the lipid mixture. We show that the results from cryo-transmission electron microscopy correlate well with those obtained from dynamic light scattering and that the disks are well described by an ideal disk model. Increasing the temperature, from 25°C to above the gel-to-liquid crystalline phase transition temperature for the respective lipid mixtures, has a relatively small effect on the aggregate structure. Abstract | Full Text | PDF (436 kb)

• …more

Copyright © 2007 The Biophysical Society. All rights reserved.
Biophysical Journal, Volume 92, Issue 3, 903-912, 1 February 2007

doi:10.1529/biophysj.106.093013

Membranes

Self-Assembling of Peptide/Membrane Complexes by Atomistic Molecular Dynamics Simulations

Santi Esteban-Martín and Jesús Salgado^,  

Institute of Molecular Science, University of Valencia, 46980 Paterna, Valencia, Spain

Address reprint requests to Jesús Salgado, Instituto de Ciencia Molecular, Universitat de València, Polígono de La Coma s/n 46980 Paterna, Valencia, Spain. Tel.: 34-96-3543016.

Biological membranes are highly dynamic supramolecular complexes composed mainly of weakly interacting lipid and protein molecules. Because of the intrinsic disorder of their biologically relevant liquid-crystalline state, experimental methods encounter severe limitations for obtaining models of membranes at atomic detail 1. In contrast, the fast molecular motions of bilayer lipids make these complex structures particularly amenable to molecular dynamics (MD) simulations 2,3,4. Such computational methods allow the description of the spatial organization and temporal dynamics of the system, also providing a useful framework to interpret experimental results.

Among the biomembrane systems for which MD simulations can be applied, membrane/protein complexes are of special interest, due to the importance of protein-lipid interactions in understanding numerous biological functions. Particularly relevant are the fundamental processes of protein insertion and folding in the lipid membrane and the assembly of transmembrane protein segments. These have been studied by MD using simplified models of hydrophobic or amphipathic peptides and lipid bilayers, and binding of peptides to the membrane/water interface, their insertion into the bilayer and packing of peptide α-helices have been investigated (as recent examples see 5,6,7,8). Typical strategies involve the use of, often arbitrary, preformed membrane systems, where several lipids are removed from one or the two leaflets to accommodate the protein inclusion. However, such preassembled systems present important limitations: i), They are difficult to build without imparting unrealistic bilayer stress. ii), The initial position of the peptide with respect to the membrane, i.e., whether it is bound to the surface or inserted, its tilt angle or the depth of insertion is usually unknown and must be assumed for the starting configuration. iii), In the peptide/membrane complexes, the diffusion of the helix from its initial position is slow and limits the number of configurations and binding modes accessible during commonly used simulation times.

As an alternative to the preformed systems, the spontaneous aggregation of bilayer components provides a way to obtain unbiased lipid membrane systems 9. A self-aggregation strategy has been tried with success to generate micellar complexes, either of surfactant molecules alone 10 or in the presence of hydrophobic peptides or proteins 11,12,13. Similarly, the spontaneous formation of pure lipid bilayers by means of MD simulations has also been described 9,14. Guided by these studies, we have generated in this work self-assembled lipid-bilayer/peptide supramolecular complexes by allowing a number of randomly distributed dipalmitoyl-phosphatidylcholine (DPPC) lipids and an α-helical peptide to self-organize freely in an aqueous environment. We have chosen two model peptides flanked by Trp-anchoring residues and with different hydrophobicities, depending on the presence of a poly-Leu or a poly-Ala central sequence: CH₃-CO-W₂L₁₈W₂-NH₂ (WL22) or CH₃-CO-W₂A₁₈W₂-NH₂ (WA22), respectively, expected to be α-helical in a membrane environment. The molecules self-organize freely in the simulation box, giving characteristic peptide/membrane complexes. The assembly process occurs through a series of steps that are kinetically distinguishable but depends barely on the type of peptide or the initial or final peptide positions. In the final complexes, the two simulated peptides display characteristic behaviors with respect to their membrane-binding mode. Thus, although the very hydrophobic WL22 is always inserted across the lipid bilayer, the more polar WA22 is found in two alternative configurations: either transmembrane (TM) inserted or embedded close to the interface region, parallel to the membrane plane. These findings show that current MD simulations are able to predict the correct binding mode of characteristic model peptides depending on their sequence. Moreover, the strategy also reports the response of the lipid bilayer to the different ways of peptide binding.

All simulations were performed using the GROMACS suite of programs, version 3.2 15. The united atom lipid parameters were adapted from the work of Berger and co-workers 16 and the peptide used the GROMOS force field. For water, the single point charge (SPC) model 17 was used, which has been shown to behave well in lipid bilayer/water simulations 18. The simulations were carried out using periodic boundary conditions with constant pressure and temperature. A Berendsen thermostat 19, with a coupling constant of 0.1ps, was used. The reference temperature was set to 323K, well above the phase transition temperature of DPPC (315K). Pressure coupling was applied anisotropically, also using the Berendsen scheme 19, with a coupling constant of 1.0ps. The reference pressure was 1 bar in all directions. Lipids, solvent, and peptide were separately coupled to the temperature bath. Simulations were run with a 4-fs time step. Bond lengths were constrained using the LINCS algorithm 20. Short-range electrostatic and Lennard-Jones interactions were cut off at 1.0nm, and long-range electrostatic interactions were calculated by using the particle mesh Ewald algorithm 21.

Coordinates for the DPPC molecules were taken from http://moose.bio.ucalgary.ca/. DPPC has been extensively studied by MD methods using force fields and setups similar to ours 22,23,24,25. The coordinates for the peptides CH₃-CO-W₂L₁₈W₂-NH₂ (WL22) and CH₃-CO-W₂A₁₈W₂-NH₂ (WA22) in α-helical structures were generated using standard tools of the WHATIF software package 26.

The starting configurations were prepared through several steps: i), A peptide aligned with the z axis was placed in the center of a cubic box (8.0-nm edge length). ii), A total of 128 DPPC molecules were distributed in the box. For that, a lipid molecule was chosen randomly from a pool of different configurations and was randomly translated and rotated. The chosen lipid was placed in the box at any available space, using as many attempts as necessary. iii), The system was solvated with 6,000 SPC molecules. iv), Then, the system was energy minimized and an MD run with position restrained on the peptide backbone was performed. At this stage isotropic pressure coupling was used, shrinking the box to an edge length of 7.0nm in ∼50ps. v), Finally, the production run was performed by allowing the system to evolve freely under anisotropic pressure coupling. Similar protocols have been used before 27.

Trajectories were visualized with the help of the VMD program 28. Lipid clusters were defined using a general criterion based on the distance between the centers of mass of the lipids (or the lipids and the peptide). Two lipid molecules are defined to be in the same cluster if the distance between their centers of mass is smaller than a cutoff distance. The cutoff distance is chosen as that which detects two clusters (one per monolayer) along a trajectory of a simulated bilayer made of 128 DPPC lipids. We find acceptable cutoff values within the range 0.9nm<cutoff<1.2nm. A value of 1.1nm was used for the analysis, although the results do not depend significantly on the precise choice of the cutoff within the indicated range. A similar strategy has been applied to cluster dodecyl phosphatidylcholine (DPC) molecules while aggregating into a unique micelle 10.

The state of lipid organization throughout self-assembly at a given time is reflected by an orientational order parameter of the lipid acyl chains, S_l, defined as the average value of the instantaneous molecular order parameters, S_mol(n), of each n segment of the chain 29. In turn, S_mol(n) is calculated from S_mol(n)=1/2; 〈3cos²θ_n-1〉, where θ_n is the angle between the nth segmental vector linking carbon atoms C_n−1 and C_n+1 in the acyl chain and the normal of the membrane (eventually formed during the assembly process), and the brackets denote an ensemble average calculated for each frame.

The tilt of lipid acyl chains is the angle formed between the membrane normal and unit vectors pointing from the midpoint of C1 and C2 toward the midpoint of C15 and C16 30. The peptide tilt is the angle between the molecular axis of the peptide α-helix and the membrane normal.

The starting system consisted of a cubic box with an 8-nm edge, filled with 128 DPPC lipids, an α-helical peptide (WL22 or WA22), and 6,000 water molecules. The lipids were placed with random conformations of their acyl chains, random orientations, and random rotations about their long axis (see Methods for more details). Because initially all three main axes are equivalent, the peptide helix was placed aligned with the z axis. We found that this initial alignment of the peptide influenced neither the direction eventually chosen by the system as the director axis (the bilayer normal) nor the final position occupied by the peptide with respect to the lipid membrane. A total of eight simulations of peptide/membrane complexes were performed, plus two control simulations with only lipids and water, all of them listed in Table 1.

For the sake of clarity, we base our description mainly on simulations 1 and 7, as they correspond to paradigmatic cases of TM-inserted WL22 and interface-bound WA22, respectively. The most important facts corresponding to each simulation are summarized in Table 1. Additionally, comparisons are made between the different cases when worth commenting. In all cases studied the self-assembling process evolves at an irregular pace, similar to that described for the self-formation of pure DPPC bilayers 9 and DPC micelles 10. We can distinguish during this process a number of characteristic stages, as follows:

In the starting random configuration of the lipids/peptide/water mixture (Figure 1aae and Figure 2aae), there is a large total hydrophobic surface exposed to the polar solvent, which drives the rapid reorganization of the system toward aggregation of the peptide and lipid molecules and exclusion of water. Thus, at the beginning of the simulation the exposed area of lipids is quickly reduced (Figure 3a). However, the accessible surface of the peptide is initially low and does not suffer important variations until the bilayer is completely formed (Figure 3b, see also below). This is a consequence of the crowdedness in the starting random system, where the peptide is closely surrounded by lipids facilitating it acting as an aggregation nucleus. The hydrophobic effect gives rise to the rapid formation of a first generation of small, irregular lipid groupings, together with a large cluster where the peptide recruits 20–40 lipids (Figure 4c). The number of clusters decreases slowly during this stage, and it appears that the initial reduction of the lipid apolar surface is mainly due to the reorganization of the first set of clusters, rather than to their fusion into larger ones (Figure 4ab). It is interesting to note that despite the facilitated nucleation of a unique large cluster in the presence of the peptide, no important global kinetic variations were observed with respect to the aggregation of only lipids in control simulations (not shown).

Display large version of this figure

Figure 1

Snapshots of the spontaneous aggregation of a mixture of DPPC lipids, water, and the WL22 peptide (simulation 1 of Table 1). Headgroup atoms are depicted blue, atoms of the lipid tails are depicted light gray, and water molecules are drawn in green. The backbone of the peptide is shown in a simplified tube representation with dark gray color. The side chains of Trp residues are in red. For clarity, the figures do not exactly correspond to actual simulation boxes but show part of the system repeated in space. The initial random distribution of molecules (a) evolves into micelle-like clusters that concentrate in a distinct area, where water starts to be excluded ((b), 1.5ns). Extensive fusion and ordering yields a metastable bilayer structure where the two leaflets are fused at the level of a transbilayer lipid pore, here shown at time 35ns (c). The pore eventually closes (at 37ns), and the membrane is further equilibrated up to 50ns (d). During the aggregation process, the peptide accompanies the lipids through hydrophobic interactions, first as part of a big, micelle-like cluster (b). As the primordial bilayer forms, the peptide reorients (c) and keeps inserted across the membrane for the rest of the simulation (d).

Display large version of this figure

Figure 2

Density distribution of characteristic groups during the self-assembling of peptide/bilayer complexes, resolved along the direction normal to the membrane formed at the end of the process. Boxes a–d and e–h correspond to simulations 1 (complex with the WL22 peptide) and 7 (complex with the WA22 peptide), respectively. The solid line marks densities of lipid headgroup atoms, dotted lines are for lipid acyl tails, dashed lines for water molecules, and dotted-dashed lines for peptide atoms. Four different averages are shown, with labels corresponding to the same stages as in Fig. 1: (a and e), time interval 0–200ps; (b and f), time 1–2ns; (c and g), time 30–40ns; and (d and h), final equilibrated complexes (time 45–50ns and 135–140ns, respectively). The positional order characteristic of a lipid bilayer can be seen in c, d, g, and h, although in the first two cases the water density across the membrane indicates the presence of a pore. The density of peptide atoms in the final stages shows clearly its position across the membrane (WL22, (d)) or bound parallel to the membrane reaching both the hydrocarbon and interface regions (WA22, (h)).

Display large version of this figure

Figure 3

Water accessible surface for hydrophobic groups in lipid (a) and WL22 peptide (b) molecules in simulation 1. (a) Lipid acyl tail accessibility: the exposed carbon tail surface is quickly lowered within the first nanosecond (inset, logarithmic timescale) and subsequent reduction attenuates. A final small step toward lower accessibility occurs after ∼15ns, due to the ordering of lipids into a liquid crystal bilayer. (b) Hydrophobic peptide accessibility: From the beginning, the peptide is closely surrounded by nonordered lipids, which largely excludes water and keeps accessible area low and fairly constant during the first ∼30ns. After this time, accessibility decreases to a lower rung, coinciding with a fluctuation of the peptide orientation (Figure 6b) just before the pore closes (37ns).

Display large version of this figure

Figure 4

Evolution of the number of clusters throughout the self-assembly process for simulations 10 (no peptide, (a)) and 1 (with peptide WL22, (b and c)). The inset in b corresponds to the number of clusters of simulation 1 but calculated considering only lipids. Two lipids (or a peptide and a lipid) are defined to be in the same cluster if the distance between their center of mass is smaller than 1.1nm. (a and b) The number of clusters decreases at an irregular pace. At the end of the process two clusters are formed, corresponding to the two monolayers (a and inset in b), which reduce to one cluster if the peptide is considered in the peptide/membrane complex (b). (c) Size of the cluster including the peptide. The peptide cluster size does not increase significantly until the first nanosecond, when small and irregular lipid-only clusters coalesce massively into it.

After the first 200ps the fusion of clusters typically increases, reducing their total number (Figure 4ab) and giving rise to bigger aggregates with micelle-like structure (Figure 1b). In parallel, all lipid groupings tend to concentrate in a centered region from which water is gradually excluded (Figure 1bbf and Figure 2bbf). In this lipid-rich area, fusion events are more probable and most lipid domains coalesce quickly into the biggest, peptide-containing cluster (Figure 4bc), which in turn accelerates the aggregation process and accentuates the confinement of water and lipids to distinct regions of space. At ∼3ns only a few lipids remain outside a central big cluster to which the peptide is bound (Figure 4c). Although in this aggregate there appears to be some bilayer-like patches, the micellar organization still dominates, as is witnessed by the low acyl-chain lipid order parameter, S_l (Fig. 5; see Methods for a definition of S_l). A higher level ordering of the phospholipids is a much slower process, which we distinguish as a separate stage.

Display large version of this figure

Figure 5

Formation of the liquid-crystal bilayer seen as the time evolution of an order parameter of lipid acyl tails (S_l, see Methods). The increase of S_l proceeds in a step-wise manner, with fast and slow phases overlapping with other important events of self-assembly (see text).

In the big main cluster, the lipids tend to be confined to a centered layer of the box (Figure 1c), increasing intermolecular connectivities in the two directions of a plane and decreasing connectivities in the direction perpendicular to that plane. This corresponds to an improvement of the positional order of the lipids (Figure 2cg) that occurs simultaneously to an increase of their orientational order (Fig. 5). As a consequence, lipid packing improves and the dimensions of the box readjust into prism shapes, allowed by the anisotropic pressure-coupling scheme used during the simulation 9,14.

The time evolution of S_l indicates formation of a liquid-crystal bilayer at ∼30ns (Fig. 5, see a more detailed description below). This reorganization also affects the hydrophobic surface accessibility, which is further reduced during the same time window (15–25ns, Figure 3a). At the end of this stage, the aggregated structure consists of a deformed lipid bilayer traversed by a water-filled lipid pore (Figure 1ccg and Figure 2ccg). The peptide molecule orders together with the lipids and in most simulations at this stage is already aligned close to its final position with respect to the membrane normal, whether TM inserted (WL22 in simulations 1–3, and WA22 in simulations 4–6) or lying flat in the bilayer plane (WA22 in simulations 7 and 8). Such orientations involved in some cases a large realignment of the peptide from its initial position (see below).

A bilayer with a water-filled lipid pore has been observed before as a characteristic intermediate state of the self-organization of lipid bilayers in MD simulations 9. It constitutes a metastable structure with variable lifetime and shape and a radius of 1–2nm. For the pore to close, the lipids forming its wall must flip from their position and water must leave the volume corresponding to the center of the bilayer. This involves passing through a high energy state which makes the process of pore closure slow, being the limiting step for complete bilayer formation 9. Under these conditions, the pore may exist for a variable time, spanning from 40 to 100ns. We did not observe significant differences in the mechanism of pore closure or in the pore lifetimes that can be attributed to the presence of a peptide during the self-assembly process (see Table 1). On the other hand, the metastable pore facilitates the equilibration of the primordial bilayer, sometimes deformed by an asymmetric number of lipids per monolayer. Although under normal circumstances the transbilayer redistribution of lipids is very slow, here it is facilitated by lateral diffusion at the pore wall. Such an equilibration is accompanied by fluctuations of the pore shape and size, rearrangements of the box dimensions, and a better organization of the bilayer lipids. Nevertheless, in most cases the pore closes before complete equilibration of the number of lipids, which suggests some tolerance for intermonolayer asymmetry (Table 1).

A number of poly-Ala-, poly-Leu-, and Leu-Ala-based peptides have been shown to adopt predominantly α-helical conformations in membrane environments 31,32,33,34,35,36. Thus, for both peptides, WL22 and WA22, we use α-helical rods as the starting structures. Although these should be barely stable in aqueous media, the peptides are well surrounded by lipids from the very beginning of the simulation, as we have seen above (Figure 4c). This helps maintain an organized secondary structure throughout the simulations, with typically only 1–2 unfolded residues at the N- or C-terminal parts. However, in some cases we observed a partial loss of helical content. For instance, in simulation 2, corresponding to WL22, a kinked peptide appears as a result of early interactions with two large micelle-like clusters, but the complete helix structure is recovered after a few nanoseconds. On the other hand, during complex formation in simulation 8, WA22 curves as it interacts interfacially at the pore wall in the primordial bilayer. This deformation persists after the pore closes and the peptide shifts to the interface of one monolayer, giving rise to a kinked helix in the final configuration (not shown).

The assembly of peptide/membrane complexes from randomly oriented molecules involves increasing positional (in plane) and orientational order to form a smectic liquid crystal lamellar phase. As described above, in our simulated systems ordering occurs mainly after coalescence of most lipid clusters into a single central big cluster. This can be seen after the increase of the order parameter of lipids from an initial value of S_l ∼ 0, corresponding to a random distribution, to a value of S_l ∼ 0.35 when the membrane is formed (Fig. 5). During the complete ordering process there is typically a characteristic fast phase at ∼20ns coinciding with the collapse of micelle-like clusters to form the lipid bilayer.

One can also follow the evolution with time of the tilt of a molecular axis for lipids with respect to the director axis, taken as the normal of the bilayer that is formed at the end of the process. As a representative lipid molecular axis, we take that of the hydrocarbon tails, which makes an angle with the reference normal axis. Because this angle averages close to 0 throughout the simulation, we take the most populated value of the distribution of () as a more informative magnitude (see Figure 6a for a typical distribution of lipid tilt angles). Illustrative examples are shown in Figure 6bc, where we analyze the evolution of the lipid and peptide orientations in simulations 1 and 7, respectively. Starting from a random arrangement of lipids, evolves slowly from ∼90° (averaged over 256 lipid tails) to ∼31°. This latter value is characteristic of an equilibrated lipid bilayer in all our simulations and is in agreement with other reported simulation data 30,37.

Display large version of this figure

Figure 6

Tilt angle of lipid chains () and peptide molecules () during self-assembly. (a) Distribution of of sn1 and sn2 lipid acyl-chains in a 10-ns time window of the equilibrated WL22/membrane complex corresponding to simulation 1. Time evolution of (black) and the mean value of the distribution of (gray) for simulations 1 (b) and 7 (c).

It is also interesting to analyze the evolution of the orientation of the helical peptides with respect to the bilayer normal, , along the simulation. In the case shown in Figure 6b corresponding to WL22, the bilayer forms in the y-z plane, with the normal being the x axis of the system. The peptide starts tilted 90° and, as the simulation proceeds, it rotates to accommodate itself in a TM fashion with ±10° fluctuations from the equilibrium value. Comparable complexes, always with a TM peptide, are produced in two other simulations with WL22. However for WA22 we observe a variable behavior. Simulations 4–6 produce complexes with a TM-bound peptide in much the same way as for WL22. In contrast, simulations 7 and 8 give peptides bound parallel to the membrane plane in the final configurations. The time evolution of the peptide tilt for simulation 7 is shown in Figure 6c. In this case the bilayer normal corresponds to the y axis and the peptide is initially perpendicular to this direction. As soon as the lipids start to aggregate, the peptide rotates slightly, but it again positions perpendicular to y as the primordial lipid bilayer is formed (at time 15ns in this case). During the rest of the simulation the peptide keeps essentially parallel to the membrane with ±10° fluctuations. Therefore, it appears that the alignment of peptides in the membrane is guided by the orientation of lipids but depends also on the peptide sequence.

Atomistic simulations of self-assembly of surfactant molecules have proved successful for generating detailed and realistic models of micelles 10,38, protein/micelle complexes 11,12,13, and pure lipid bilayers 9,14. Here we apply an atomistic MD simulation approach to investigate the spontaneous formation of peptide/lipid bilayer complexes by allowing randomly distributed DPPC lipids and a peptide molecule to self-organize freely in an aqueous environment. Under the strong thermodynamic gradient imposed mainly by the hydrophobic effect, formation of peptide/membrane complexes proceeds fast and allows us to obtain equilibrium configurations in timescales of ∼100ns. The autoorganized systems so produced are not biased by a chosen starting configuration. Instead, they are expected to depend only on objective factors, like the quality of the force field or details of the simulation methodology, and should be valid for selecting the preferred binding mode of peptides according to their physiochemical properties, codified in their sequence. An added value of this strategy is the inherent mechanistic information that can be obtained from the time evolution of the system during the formation of the complex, which helps us understand the molecular interactions defining each final configuration.

The simulated spontaneous formation of detergent micelles and liquid-crystal lipid bilayers follows similar mechanisms during the initial steps of the aggregation processes. However, the final organization is more complex for the lipid phases, often showing diverse structures and metastable intermediate states 9,14. Pure DPC and mixtures corresponding to the human bile aggregate into small micelles in ∼3ns. At time ∼10ns the small clusters coalesce to form large micelles of ∼50 molecules, which then reorganize internally and become spherical during time 10–20ns 10,38. This has been extended to the self-assembly of surfactants together with hydrophobic proteins, with examples showing the formation of a sodium dodecyl sulfate micelle around dimeric glycophorin 11,12 or the β-barrel porins OmpF and OmpX 11,13. Similarly, in a series of landmark studies Marrink and co-workers have reported the spontaneous formation of fluid lipid bilayers. DPPC 9 or a mixture of dioleoylphosphatidylcholine (DOPC) and dioleoylphosphatidylethanolamine (DOPE) 27 form a large cluster within 2–3ns, which then reorders internally to form a bilayer crossed by a metastable water-filed pore of variable size and stability. Eventually, the pore shrinks and collapses, completing bilayer formation within ∼50ns from the starting random distribution of lipids in water 9.

For the work reported here we included a hydrophobic peptide in the initial random mixture of water and lipids. From the beginning, the peptide facilitates the formation of a unique large cluster (Figure 4c), which is not observed in the absence of peptides (control simulations 9 and 10, Table 1). Similar nucleation effects have been observed during the spontaneous formation of protein/micelle complexes 11,12. However, in our simulations of peptide/membrane complexes, the presence of the peptide does not alter significantly the global kinetics of self-assembling or the principal landmarks of this process, regardless of the type of peptide or the configuration finally achieved. Thus, as in the absence of peptides, after a rapid formation of the first generation of clusters, in ∼200ps (Figure 1abf and Figure 2abf), coalescence into a big singular cluster is an order of magnitude slower (∼3ns), and the increase of molecular order to form a primordial bilayer does not happen before ∼30ns (Figure 1c and Figure 5c). This indicates that the process is limited by a number of slow phases, like diffusion of clusters, coalescence to form bigger clusters, and internal reorganizations within each cluster, including the final ordering of the bilayer, which appears to be independent of the presence of the peptide.

Again, similar to that reported in the absence of peptides, the bilayer organizes first as a metastable structure characterized by the presence of a lipidic pore, and pore closure is the rate-limiting step toward formation of a defect-free membrane 9. Previous studies have shown that the stability of the pore can be increased under specific stress conditions, up to a lifetime of more than 150ns 39. Interestingly, we do not find significant variations of pore lifetime in the presence of WL22 or WA22 regardless of their binding mode. Both in the absence and in the presence of peptides, the pore seems to act as a mechanism for relaxing the stress that often builds up as the primordial bilayer organizes, coming from an asymmetric distribution of molecules between the two leaflets. However, in most cases the pore is closed before complete equilibration of the intermonolayer number of lipids, allowing maximum differences of three lipids in the absence of peptides, four lipids in complexes with a TM-bound peptide (WL22 or WA22), or six lipids for cases of a WA22 peptide bound parallel to the membrane (Table 1). For membranes with no peptide and with a TM peptide, this indicates a tolerance of up to 6% asymmetric area expansion. On the other hand, in the two cases where WA22 binds asymmetrically to only one monolayer (simulations 7 and 8), the defect of six lipids is found always at the side of the peptide. Interestingly, in these latter complexes, for an area per lipid of 0.65nm² calculated from the peptide-free monolayer, a defect of six lipids corresponds to 3.90nm², which fits well with an excess area of ∼3.96nm² contributed by one WA22 molecule at the opposite monolayer (estimated by assuming an ideal alanine-based α-helix of 22 residues).

One of our goals with this study was to test the atomistic self-assembling MD methods for the production of models of characteristic peptide/membrane complexes, sensitive to relevant physicochemical properties of the peptide, like hydrophobicity. To this aim, we chose two simple model membrane-binding α-helical peptides, WL22 and WA22. In both cases, the terminal Trp residues provide well-known interfacial anchoring of the peptides to the membrane 40,41,42,43,44. The core 18 Leu residues of WL22, compared to the 18 Ala residues of WA22, make the first peptide much more hydrophobic than the second 45, allowing us to test the influence of this property in the observed peptide-membrane-binding modes. To provide for a lipid environment, we chose DPPC for which the bilayer hydrophobic thickness 1 should match the hydrophobic length of our peptides.

Although leucine is clearly a very hydrophobic residue, hydrophobicity scales consider alanine from moderately apolar 46,47,48 to even slightly polar, if one accounts for the contribution of the hydrophilic peptide backbone 40,45. Such ambiguity of Ala is supported both experimentally 40,49,50,51,52,53,54 and from free energy calculations, which predict poly-Ala-based peptides in transmembrane- (TM-) and surface-bound configurations with similar probabilities 55. Thus, we take the poly-Leu-based WL22 as a paradigm for strongly interacting TM peptides and the poly-Ala-based WA22 as an ambivalent case with alternative binding modes. In agreement with such expectations, after simulated self-assembling, WL22 showed an absolute preference for TM binding. In contrast, WA22 was found in two distinct configurations: TM in three out of five simulations, and parallel to the membrane in two simulations (Table 1). Before discussing these facts in more detail, we elaborate briefly on their statistical significance. In principle it could be argued that chances are that all simulations, including cases of WL22 and WA22 (eight in total), are drawn from the same distribution. This would give a relatively high probability for a random sample of three measurements, as for WL22, to be exclusively TM (0.75³=0.4). However, we find this latter possibility very unlikely, as any additional self-assembly simulation run for peptides having in common a central 18-Leu stretch (up to a total of six simulations; S. Estaban-Martin and J. Salgado, unpublished) gave a TM configuration. Thus, the cases of poly-Leu- and poly-Ala-based peptides should follow different distributions, and their differentiated behavior through the simulated self-assembly process appears to reflect the underlying thermodynamic driving forces relevant for their membrane-binding state.

The ratio of alternative states found here for WA22 is in agreement with theoretical predictions from free energy calculations 55 and with results obtained for similar poly-Ala-based peptides by solid-state ¹⁵N-NMR 52. The latter study shows an equilibrium between in-plane and TM states for K₃A₁₈K₃ peptides, similar to our WA22 peptide, where the TM alignment can be stabilized by replacing a few Ala residues with Leu. Thus, it appears that the hydrophobicity of poly-Ala peptides is close to the threshold for TM insertion 52. Similarly, a recent investigation of the integration of GGPG-(L_nA_19−n)-GPGG peptides in endoplasmic reticulum membranes using the natural translocon machinery 56 yields a small, positive apparent free energy of insertion per Ala residue, that can be counteracted by a few Leu residues in the peptide sequence (). For the case of WA22 studied here, we should also consider the contribution of the flanking Trp residues to the free energy of insertion in the membrane. In agreement with the preference of these latter residues for interfaces, it has been shown that their contribution to a TM state is strongly dependent on their position with respect to the center of the helix, with maximum reduction of total ΔG_app for Trp-to-Trp separation ≥10 residues 56. Therefore, the flanking Trp residues in WA22 can effectively counteract the marginal polarity of the central Ala₁₈ stretch, and the peptide can be stable both in TM and parallel to the membrane-binding modes, as shown by the simulations. It is interesting that despite the change in orientation of the peptide between the two states, the central Ala-stretch and flanking Trp residues are found in similar membrane regions in the two alternative configurations. Thus, even in the parallel-binding state, the peptide is immersed and most of its volume occupies the hydrocarbon region (Figure 1h), whereas the Trp residues always reside close to the membrane interface.

In summary, self-assembly of peptide/membrane complexes by MD is able to distinguish the preferred binding mode of two model peptides depending on physicochemical properties codified by their sequence. Moreover, it correctly predicts available configurations for the WA22 peptide, which are almost equivalent in terms of their thermodynamic stability.

The spontaneous assembly, starting from random mixtures of DPPC lipids with a hydrophobic peptide in water by means of MD simulations, is an unbiased method for generating peptide/membrane complexes. The aggregation process proceeds through a number of distinct steps, characterized by the formation of lipid clusters of growing size and increasing positional and orientational order, until a lipid bilayer with a peptide bound to it is finally formed. The number of clusters change discontinuously in two main phases: a first one consisting of an increased packing of small micelle-like clusters, with an initial reduction of the accessible surface area due to the internal reorganization of existent clusters; and a second one characterized by the massive fusion of the clusters, accompanied by a further reduction of the hydrophobic lipid accessibility. After the main aggregate is formed, the reorganization of lipids to constitute a liquid-crystal bilayer and the closure of a lipid pore, which always accompanies bilayer formation, are slower processes and constitute rate-limiting steps.

The peptide binds strongly to the lipids from the beginning of the simulation, which facilitates aggregation of an initial big cluster. However, the mechanism and kinetics of formation of the lipid bilayer do not change significantly with respect to a system in the absence of the peptide. Likewise, the peptides do not seem to stabilize the intermediate metastable state characterized by a lipidic pore. The fact that the diffusion of lipids through the pore wall is faster than the lifetime of the pore itself allows for an efficient mechanism to compensate the tension induced by either an initial asymmetric random distribution of lipids, asymmetric inserted peptides, or peptides adsorbed at one monolayer interface. Particularly, when WA22 binds parallel to the membrane in only one monolayer, the surface occupied by the peptide is compensated by an equivalent reduction of the number of lipids.

During the self-assembly process the peptides are free to accommodate to their preferred configuration in the emerging bilayer, depending on their physiochemical properties, such as hydrophobicity, codified in their sequence. Thus, the very hydrophobic WL22 acquires a TM-inserted state, whereas the borderline apolar WA22 is found in TM-bound and parallel-bound states. Such a distinction of states with marginal differences of stability cannot be made starting from preformed bilayers, which makes the spontaneous assembly of peptides and lipids an unbiased strategy to produce reliable models of equilibrated peptide/lipid complexes. This opens the possibility of studying peptide/membrane complexes with peptides of unknown membrane-binding mode and topology, as well as systems where the equilibrium configuration depends on complex dynamic processes, like pore-forming peptides.