Eol Han and
Hwankyu Lee*
Department of Chemical Engineering, Dankook University, Yongin, 448-701, South Korea. E-mail: leeh@dankook.ac.kr
First published on 1st December 2014
We performed coarse-grained molecular dynamics simulations of antimicrobial peptides PGLa and magainin 2 in lipid bilayers. PGLa peptides or mixtures of PGLa and magainin 2 were initially widely spaced or clustered above the bilayer surface with different heterodimeric orientations (parallel or antiparallel). Simulations show that the presence of magainin 2 promotes more tilting and insertion of PGLa into the bilayer, indicating the synergistic effect. Magainin 2 interact with lipid headgroups and thus stay horizontally on the bilayer surface, while PGLa insert into the bilayer, leading to more tilted conformation, in agreement with recent NMR experiments. In particular, for the systems with the initially antiparallel-oriented heterodimers or with the neutrally mutated magainin 2, much fewer parallel heterodimers form, and PGLa peptides are less tilted and inserted, indicating that the formation of parallel heterodimers is important for the PGLa insertion, as suggested in experiments. Peptides aggregate in the mixture of PGLa and magainin 2, but not in the system without magainin 2, indicating that magainin 2 induce the peptide aggregation, which is required for the pore formation. These simulation findings agree with the experimental observations of the heterodimer formation as well as different positions of PGLa and magainin 2 in the bilayer, which seem to conflict. These conflicting results can be explained by the synergistic mechanism that magainin 2 form parallel heterodimers with PGLa and induce the aggregation of heterodimers, leading to the formation of pores, where magainin 2 tend to interact with the bilayer surface, while PGLa are tilted and inserted into the hydrophobic region of the bilayer.
Williams et al. found that a mixture of magainin 2 and PGLa induces much stronger antimicrobial activity than does either magainin 2 or PGLa, indicating their synergistic effect.23 Also, Vaz Gomes et al. and Westerhoff et al. found that magainin 2 and PGLa form the complex, yielding the greatest synergistic effect at a molar ratio of 1:
1.24,25 To understand the mechanism for this synergistic effect, the Matsuzaki group characterized the conformation of PGLa–magainin complexes and their interactions with lipid bilayers, which showed that parallel-oriented heterodimers significantly increase the pore formation and stability, suggesting that heterodimers of PGLa and magainin 2 induce pore formation in lipid bilayers.26–28 To investigate the position and orientation of the PGLa–magainin 2 heterodimers in the bilayer, Ulrich and coworkers performed solid-state NMR experiments and showed that the presence of magainin 2 increases tilt angles of PGLa, indicating the deeper insertion of the C-terminal group of PGLa.29,30 In particular, they observed that magainin 2 peptides mostly stay on the bilayer surface, while PGLa peptides insert into the hydrophobic region of the bilayer,31 implying that the inserted peptides may not retain the side-by-side heterodimeric structure, which was also supported by Salnikov and Bechinger's experiments.32 They also showed that these synergistic effects on the PGLa insertion depend on membrane curvature (positive vs. negative) rather than membrane thickness, which conflicts with Salnikov and Bechinger's experiments that suggested that magainin 2 may reduce membrane thickness and thus promote the deeper PGLa insertion and pore formation. These experiments clearly showed that magainin 2 and PGLa form the parallel heterodimer and promote the deeper insertion of PGLa into the bilayer, but also showed different positions of PGLa and magainin 2, suggesting that magainin 2 and PGLa may not retain the side-by-side heterodimeric structure. Thus, these conflicting experimental observations and the mechanism for the synergistic effect are still not well understood. To resolve this, the interaction between those AMPs and lipid bilayers needs to be studied at nearly the atomic scale.
Molecular dynamics (MD) simulations have been able to explore the interactions of magainin 2 and PGLa with lipid bilayers. Simulations have shown that magainin peptides bind to the bilayer surface and insert into the lipid-tail region, which disorder lipids and thus induce pore formation, depending on the structure, sequence, and concentration of peptides,33–41 ion concentration,42 and polymer conjugation.43 In particular, coarse-grained (CG) MD simulations captured the experimentally proposed toroidal pores.36–39 For PGLa, Ulrich and coworker's simulations recently showed the formation of PGLa monomers and dimers, respectively, at low and high concentrations.44 They also found that dimers have the higher tilt angles than do monomers, indicating the effect of the dimeric structure on the insertion into the bilayer. These MD simulations have revealed atomic-scale insights into the structure and dynamics of either PGLa or magainin 2 and its interaction with the lipid bilayer, but the mixtures of PGLa and magainin 2 have not yet been simulated to understand their synergistic effect on the insertion and pore formation.
In this study, we therefore perform MD simulations of the mixture of magainin 2 and PGLa in lipid bilayers to understand their synergistic effects on the insertion into the bilayer. The PGLa–magainin 2 dimers were initially aligned in parallel or antiparallel, and the systems with mutated magainin 2 or without magainin 2 are also simulated. The orientation and position of peptides in lipid bilayers are analyzed by calculating the tilt angles and the distances from the bilayer center, which compare favorably with experiments. In particular, we find the dependence of the tilted conformation, aggregation, and insertion of PGLa on the presence of magainin 2 and the heterodimer orientation (parallel or antiparallel), which is rationalized by analyzing the interactions of PGLa, magainin 2, and lipids. We will show that these results help explain the conflicting experimental results regarding the dimer formation and different positions of PGLa and magainin 2.
![]() | ||
Fig. 1 Structures and helical wheel diagrams for PGLa (left) and magainin 2 (right). Gray, red and blue colors represent hydrophobic, anionic and cationic residues, respectively. For the wheel diagram, amino acid sequences are plotted clockwise. The structure images were created using Visual Molecular Dynamics.58 |
Twelve peptides (12 PGLa or a mixture of 6 PGLa and 6 magainin 2) were placed above the equilibrated bilayer surface with a distance of ∼2.4 nm between the peptide and bilayer centers, leading to the peptide/lipid molar ratio of 0.023, close to the recent NMR experimental condition.31 To check the effect of the initial configuration, 12 peptides were aligned in parallel or antiparallel with different initial distributions (widely spaced or clustered), shown in Fig. 2 and Table 1. These 12 peptides were added to both leaflets of the bilayer to obtain more sampling. The final simulation system consists of 24 peptides (12 peptides/leaflet), 512 DLPC lipids (256 lipids/leaflet), ∼10300 CG water beads (representing ∼41
200 real waters), 96–120 counterions (Cl−) in a periodic box of size 12.5 × 12.5 × 13 nm3. A temperature of 308 K and a pressure of 1 bar were maintained by applying the velocity-rescale thermostat and Berendsen barostat in the NPxyPzT ensemble (semi-isotropic pressure coupling).55,56 A real space cutoff of 12 Å was used for Lennard-Jones (LJ) and electrostatic forces. The LJ and Coulomb potentials were smoothly shifted to zero between 9 and 12 Å, and between 0 and 12 Å, respectively. All simulations were performed for 7 μs with a time step of 20 fs on computational facilities supported by the National Institute of Supercomputing and Networking/Korea Institute of Science and Technology Information with supercomputing resources including technical support (KSC-2013-C2-46).
No. of molecules | Initial orientation of peptides | Initial distribution of peptides | Simulation time (μs) | |||
---|---|---|---|---|---|---|
PGLa | Magainin 2 | DLPC | ||||
PGLa-as | 24 | — | 512 | Antiparallel | Separated | 7 |
PGLa-ac | 24 | — | 512 | Antiparallel | Clustered | 7 |
PGLa–MAG-ps | 12 | 12 | 512 | Parallel | Separated | 7 |
PGLa–MAG-pc | 12 | 12 | 512 | Parallel | Clustered | 7 |
PGLa–MAG-as | 12 | 12 | 512 | Antiparallel | Separated | 7 |
PGLa–mMAG-ps | 12 | 12 | 512 | Parallel | Separated | 7 |
The insertion of peptides into the bilayer was further quantified by calculating the distance between the bilayer center and the center of mass (COM) of each peptide in the bilayer normal direction (z-direction). In Fig. 4, all the distances reach steady-state values within 3 μs, again indicating that simulations are well equilibrated, although distances fluctuate greatly. For the systems without magainin 2, most PGLa peptides are closely located in the lipid headgroup region, while for those with magainin 2 many PGLa peptides are much closer to the bilayer center, again indicating the synergistic effect of magainin 2 on the insertion of PGLa into the bilayer. Also, distances between magainin 2 and bilayer centers are very close to the distances between DLPC phosphate and bilayer centers, indicating that magainin 2 molecules interact with the bilayer surface and do no insert, consistent with Fig. 2. These configurations of PGLa and magainin 2 in the bilayer are also confirmed by calculating mass densities. In Fig. 5, PGLa-as and PGLa-ac (systems without magainin 2) show that most PGLa peptides are positioned around the bilayer surface, while other systems with a mixture of PGLa and magainin 2 show that PGLa peptides are in the hydrophobic region of lipid bilayers, consistent with Fig. 2 and 4. Magainin 2 peptides are observed mostly in the DLPC–phosphate region, indicating that their strong electrostatic interactions with DLPC phosphates make magainin 2 stay on the bilayer surface. Fig. 5 also shows less insertion of PGLa in PGLa–MAG-as and PGLa–mMAG-ps than in PGLa–MAG-ps and PGLa–MAG-pc, indicating the heterodimer formation and orientation of PGLa and magainin 2 influence the insertion of PGLa, which will be discussed in detail later. Here, our simulations show that the presence of magainin 2 promotes the deeper insertion of PGLa into the bilayer, indicating the synergistic effect, as also observed in experiments.26,29,31,32 In particular, extents of the insertion of PGLa and magainin 2 differ, which agrees well with recent NMR experiments that showed that PGLa insert into the bilayer, while magainin 2 stay on the bilayer surface.31,32
Note that although the position and insertion extent of PGLa and magainin 2 in their mixtures agree with experiments,32 our simulations with only PGLa peptides show that PGLa mostly stay on the surface of DLPC bilayers, which is different from experiments that showed the insertion of PGLa into the hydrophobic region of the DLPC bilayer.32 This discrepancy may occur, since four methylene units are lumped into each CG tail-bead, and hence lipids with a difference of methylene units can be only semiquantitatively distinguished. To test this, we also simulated PGLa in the thinner bilayer, which has two CG beads per tail (recall that DLPC has three CG beads per tail). In Fig. 6, density profiles and configurations clearly show the much deeper insertion of PGLa into the thinner bilayer, leading to the cross-membrane configuration, as observed in experiments.32 These indicate that although the effect of the tail length cannot be accurately captured, simulations can reasonably reproduce the experimentally observed conformation of the inserted PGLa, which also supports the above-mentioned results of the insertion extent of PGLa and magainin 2 in their mixtures.
![]() | ||
Fig. 6 Mass density profiles of peptides, DLPC lipids, and Cl− ions for the systems composed of 24 PGLa peptides in the DLPC bilayer (PGLa-as; top) and those in the thinner bilayer (bottom). |
Fig. 7 shows tilt angles of the 12 individual peptides for each system as a function of time. Magainin 2 peptides show lower tilt angles than do PGLa peptides, indicating the less tilted conformation of magainin 2, apparently because magainin 2 peptides interact with lipid headgroups and thus stay horizontally on the bilayer surface, as observed in Fig. 4 and 5. All 12 PGLa peptides for each system (24 peptides for PGLa-as and PGLa-ac) were categorized according to their individual tilt angles in three different ranges of <105°, 106–120°, and >121°, which are tabulated in Table 2. PGLa–MAG-ps and PGLa–MAG-pc show that more than half of the peptides are in the tilt-angle range of >121°, while other systems show that most peptides have the lower tilt angles, indicating that the parallel orientation of dimers and the charge interaction between PGLa and magainin 2 are important to increase the tilting of peptides.
No. of PGLa in terms of the tilt angle | |||
---|---|---|---|
<105° | 106–120° | 121°< | |
PGLa-as | 9 | 15 | — |
PGLa-ac | 3 | 21 | — |
PGLa–MAG-ps | — | 5 | 7 |
PGLa–MAG-pc | 2 | 3 | 7 |
PGLa–MAG-as | 3 | 8 | 1 |
PGLa–mMAG-ps | 4 | 8 | — |
These results indicate that the presence of magainin 2 significantly increase the tilting of PGLa, in qualitative agreement with Tremouilhac et al.'s NMR experiments,29 although simulation and experimental values do not exactly match. Since the larger tilt angle of PGLa indicates deeper insertion into the bilayer, these results, combined with the observations in Fig. 4 and 5, indicate that magainin 2 interact strongly with lipid headgroups and thus stay horizontally on the bilayer surface, while PGLa insert into the bilayer, leading to the tilted orientation. In particular, PGLa is less tilted in the system with the initially antiparallel heterodimer orientation (PGLa–MAG-as) or with mutated magainin 2 (PGLa–mMAG-ps), indicating that the parallel orientation and the stronger electrostatic interaction between magainin 2 and PGLa induce deeper insertion of PGLa, which agrees with the Matsuzaki group's experiments that showed the greater extent of antimicrobial activity for parallel heterodimers.27,28
Fig. 8 shows the number of “clusters” of peptides as a function of time, where a cluster is either a complex of any size or a free peptide. Here, if the distance between beads of different peptides is less than 0.5 nm, then those are considered to be a cluster. Other criteria with the distance from 0.6 to 0.8 nm produce similar qualitative trends. Fig. 8 also shows the number of peptides in the largest cluster as a function of time. These both reach steady-state values at ∼3 μs, indicating the formation of the equilibrated complexes in the bilayer. Peptides are much more complexed in the systems with magainin 2 than in those without magainin 2 (PGLa-as and PGLa-ac). In particular, more than 7–8 peptides are clustered for each leaflet, leading to the largest complex with 14–16 peptides, similar to the pore formation, as visualized in Fig. 8. Fig. 9 shows radial distribution functions (RDFs) between anionic phosphates of DLPC and cationic residues of PGLa. The systems without magainin 2 (PGLa-as and PGLa-ac) or with mutated magainin 2 (PGLa–mMAG-ps) show the higher peaks than do others, indicating the stronger charge interactions between PGLa and lipid headgroups, which can inhibit the PGLa insertion. These imply that magainin 2 peptides induce the clustering of peptides, which may help the deeper insertion of PGLa.
![]() | ||
Fig. 8 Number of “clusters”, where a “cluster” is an aggregate of any size, including single peptides (top), and number of aggregated peptides in the largest cluster (bottom) as a function of time. |
To understand the interaction of heterodimers in the peptide complex visualized in Fig. 8, we calculated RDFs of five cationic residues of the clustered PGLa with respect to two anionic residues (glutamic acid (19Glu) and C-terminal (23Ser)) of magainin 2 in the complex. Recall from Fig. 1 that magainin 2 consists of two anionic and five cationic residues, while PGLa consists of only five cationic residues. Experiments suggested that parallel heterodimers are induced by the electrostatic interactions between anionic residues of magainin 2 and cationic residues of PGLa.29 Fig. 10 shows that the 19Glu residues of magainin 2 mostly interact with 12th, 15th, and 19th lysines (Lys) of PGLa, and that 23Ser residues of magainin 2 interact with 15th and 19th Lys of PGLa, indicating the formation of parallel heterodimers in the peptide complex, as shown in experiments.27 These results are observed only in the simulations PGLa–MAG-ps and PGLa–MAG-pc, which show the deepest insertion and the most tilted orientation of PGLa. Note that when peptides were initially aligned in antiparallel (PGLa–MAG-as) or mutated with neutral residues (PGLa–mMAG-ps), much fewer parallel dimers form, leading to less insertion and tilted orientation of PGLa. These indicate that magainin 2 peptides form parallel heterodimers with PGLa and induce the aggregation of heterodimers, where magainin 2 stay on the bilayer surface, but PGLa insert into the bilayer.
These simulation results, combined with observations in Fig. 4–7, support the experimental observations and suggestions. Parallel heterodimers of PGLa and magainin 2 promote more tilting and insertion of PGLa into the bilayer than do antiparallel heterodimers or only PGLa peptides, indicating that the parallel orientation of heterodimers is important for the insertion of PGLa, in agreement with experiments that showed the formation of magainin 2–PGLa heterodimers and their synergistic effects on insertion.27,28 Also, magainin 2 peptides tend to stay on the bilayer surface, while PGLa peptides insert into the bilayer, which also qualitatively agrees with solid-state MNR experiments.31,32 In particular, magainin 2 induces the aggregation of parallel heterodimers, which may be important for stable pore formation, as suggested from experiments.26,28,29 These findings explain the conflicting experimental observations of the heterodimer formation as well as different positions of PGLa and magainin 2 in the bilayer. PGLa and magainin 2 form parallel heterodimers and aggregate in the bilayer, leading to the formation of pores, where magainin 2 peptides tend to interact with the bilayer surface, while PGLa peptides insert into the hydrophobic region of the bilayer, which allows PGLa to be more tilted. To more carefully check the stability of parallel heterodimers in the peptide complex, all-atom simulation studies of the interaction between magainin 2 and PGLa ought to be performed, which we hope to report on elsewhere.
Although simulations capture the experimentally observed conformation and insertion extent of PGLa and magainin 2 in lipid bilayers, the pore formation was not explicitly observed over our simulation time scale. Despite this, simulations show that PGLa and magainin 2 form parallel heterodimers, and that magainin 2 and PGLa are positioned respectively on the bilayer surface and in the hydrophobic region of the bilayer, leading to the tilted conformation of PGLa. In particular, the extents of the PGLa insertion and tilting increase only in the presence of parallel heterodimers, showing the dependence of the PGLa insertion on the formation of parallel heterodimers, as observed in experiments. These indicate that leakage or even pore formation could well eventually occur in the systems with the mixture of PGLa and magainin 2, since experiments have shown that magainin 2–PGLa heterodimers and their tilting conformations induce more stable pore formation than do PGLa without magainin 2.26,28,29 Also, note that the heterodimer formation and different positions of PGLa and magainin 2 seem to conflict. Our simulations show that heterodimers aggregate, but PGLa peptides (without magainin 2) do not aggregate, indicating that magainin 2 induce the peptide aggregation, which may play an important role in pore formation. Thus, heterodimers form, aggregate, and induce pore formation, where magainin 2 are positioned on the bilayer surface, while PGLa are more tilted and inserted into the bilayer, which helps explain the conflicting experimental observations.
This journal is © The Royal Society of Chemistry 2015 |