Introduction

The development of drug resistance by many bacteria against traditional antibiotics poses an important challenge in curing infectious disease. Extensive research has been performed to develop antimicrobial peptides into powerful antibiotics to kill bacteria1,2,3,4,5,6,7. Because most antimicrobial peptides disrupt the cell membranes of bacteria, it is difficult for bacteria to develop drug resistance against antimicrobial peptides. However, the detailed interaction mechanisms between many antimicrobial peptides and bacterial cell membranes remain unclear.

LL-37, the only cathelicidin member in humans, plays an important role in human innate immunity system8,9. LL-37 exhibits a broad-spectrum antimicrobial activity and lipopolysaccharide-neutralizing effects. There is considerable therapeutic interest in utilizing LL-37 to overcome the bacterial resistance against traditional antibiotics and therefore there is significant interest in understanding its mechanism of action. Studies have reported the biological effects of LL-37 as well as the interactions of LL-37 with various types of lipid membranes. It was found that LL-37 readily disrupts the negative charged 1,2-dipalmitoyl-sn-glycero-3-phospho-(1′-rac-glycerol) (sodium salt) (DPPG) monolayer but exerts no effect on neutral charged 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) and 1,2-dipalmitoyl-sn-glycero-3-phosphoethanolamine (DPPE) monolayers by specular X-ray reflectivity10. NMR techniques have been used to determine 3D structures of LL-37 associated with neutral n-dodecylphosphocholine (DPC)11 and negative-charged sodium dodecyl sulfate (SDS) micelles12. Although the peptide structures are not completely the same in these two environments, they both feature a kink in the middle of the peptide. Solid-state NMR studies revealed the oligomerization13, membrane orientation and carpet mechanism action for the peptide14. While solid-state NMR and calorimetric studies have provided insights into the mechanism of action for LL-37, obtaining such information at very low, physiologically relevant concentrations have been a major limitation. Methods to overcome this limitation not only can be used to study other AMPs but also other membrane active peptides/proteins including cell penetrating peptides and amyloid proteins.

In this study, we demonstrate the power of the sum frequency generation (SFG) vibrational spectroscopic technique to study the interactions between LL-37 and a single lipid bilayer containing different ratios of negative-charged 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoglycerol (POPG), neutral-charged 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) and cholesterol for various peptide concentrations. SFG spectroscopy is an intrinsic surface-sensitive technique15,16,17,18,19,20,21,22 and has been widely used to study the structure and orientation of peptides and proteins at bio-interfaces23,24,25,26,27. Because of the excellent sensitivity, SFG can be used to investigate peptide-membrane interactions in the physiologically-relevant peptide concentration range, which cannot be done using most other biophysical techniques28,29,30,31.

Most of the previous orientation studies on α-helical peptides using SFG were focused on the linear α-helical structure32,33. For peptides that are not linear, data analysis on linear peptides was still used to determine orientation34. Here, we studied a non-linear α-helical structure using SFG and developed approaches to analyze the data for the first time in the literature. We considered two types of non-linear α-helical structure in this study: a bent structure and a disrupted structure. The bent structure changes in the helix axis direction with all the residues remaining helical, an example of which is LL-37 associated with SDS vesicles12. The disrupted structure also bears a change in the axial direction but with a loss of the helical character of the residues around the kink location, such as melittin, alamethicin and LL-37 associated with DPC vesicles11. For both “bend and disruption” models, we treat the helical part as two adjacent segments. The way we treat the first segment is the same as for an ideal linear helix previously reported32,33, but the calculation for the second segment is different for two models. For the bend model, since there is no random structure between the two segments, we ought to consider that the entire helical LL-37 structure is continuous. Especially, there is a phase term involved in the calculation of the E1 vibrational mode which should also be continuous for the entire LL-37 helical structure. On the contrary, the phases for the vibrations of the two segments in the disrupted structure are unrelated and we do not consider the phase continuation. More calculation details (i.e. the values for the dipole moments and the Raman tensors) are included in the supporting information. The peptide membrane orientation deduced using SFG is compared to results from NMR studies. Our results demonstrate that at the high peptide concentrations, SFG results are compatible to the NMR data. Different from NMR, SFG has the capability to investigate very low concentrations of complex α-helical structures and their membrane orientations.

Results

Structures of LL-37 in SDS12 and DPC11 micelles were reported based on solution NMR studies (Figure 1). The continuous helical region includes the residues 2–31 in SDS with a bend between residues 14–16. The bend is caused by the hydrophobic interaction between residues Ile13 and Phe17 and the membrane. The angle between the two helical segments connected by the bend is about 143° (Figure 1a). We developed the “bend model” to determine the peptide orientation in POPG lipid bilayers for our SFG study. The structure of LL-37 associated with zwitterionic DPC micelles was also reported based on a solution NMR study11. The well-structured region is from residues 4 to 33 and the helix-break-helix motif was highlighted with a break at residue Lys-12 (Figure 1b). Both structural studies suggest a tight cluster formed by residue Ile13 and Phe17. In contrast to the LL-37 structure in SDS vesicles, where the N-terminal region is rigid, the N terminus of LL-37 was found to be dynamic in DPC micelles. Since there is a disruption in LL-37 helicity when associated with zwitterionic lipid vesicles, we adopt the “disruption model” to address the disruption in the LL-37 helical structure for SFG data analysis. Since the N terminus is dynamic, we calculated the orientation dependence of the nonlinear optical susceptibility components for structures with different residue numbers in the N terminal α-helical region to ensure that our model is reliable (Figure 2). A schematic of the LL-37 molecule and the molecular axes used for SFG data analysis is shown in Figure 1c.

Figure 1
figure 1

Structures of LL-37 (a) associated with SDS vesicles12 (PDB: 2K6O) (b) associated with DPC vesicles11 (PDB: reported by Ramamoorthy group) (c) A cartoon representation of the LL-37 molecular structure and the rotation axis.

Figure 2
figure 2

(a) Dependence of the calculated SFG ratio of LL-37 associated with a POPC/POPC bilayer on the helix tilt angle (first segment) relative to the bilayer surface normal.The first segment has 21 amino acid residues and the second segment has 6, 7 and 8 amino acid residues respectively (the disruption model). (b) Dependence of the calculated SFG ratio of LL-37 associated with a POPG/POPG bilayer on the helix tilt angle (first segment) relative to the surface normal. 18–12, 17–13, 16–14 and 15–15 refer to the peptides with the first α-helical segment of 18, 17, 16 and 15 amino acid residues and the second α-helical segment of 12, 13, 14 and 15 amino acid residues respectively (the bend model).

Since LL-37 has a random coil structure for the N-terminal residues, we need to ensure that the detected SFG amide I signal is dominated by the contributions from the α-helical component. Contribution to the entire SFG signal from the unstructured region of the peptide was calculated by NLOPredict developed by the Simpson group35 for both the bent and disrupted structures of LL-37. We found that for both the membrane-surface and transmembrane orientations, the SFG signal contribution from random coil structured regions is less than 5% of the total SFG signal (supporting information, even under the assumption that the peak of the random coil overlaps with that of the α-helical components). Therefore, it is reasonable to assume that the signal attributed to the α-helical structured region of the peptide is the main contributor to the measured SFG signal in the amide I region.

SFG spectra collected with ssp and ppp polarization combinations from LL-37 associated with POPC/POPC bilayers in the amide I frequency region at 0.46 μM and 1.6 μM concentrations are shown in Figure 3a. Here, for the SFG signals in the amide I frequency range, we collected spectra using two polarization combinations, ssp and ppp, in order to deduce the LL-37 orientation from the spectral fitting result χeff,pppeff,ssp. The SFG spectra in Figure 3a exhibit a single peak centered around ~1647 cm−1, indicating an α-helical structure. The SFG spectral intensities for the two studied concentrations are very similar and the fitting results of χeff,ppp and χeff,ssp are displayed in Table 1a. Using the spectral fitting results and the “disruption model” (Figure 2a), the orientation analysis shows that the LL-37 peptide lies more or less parallel to the membrane surface with a tilt angle between 56° and 90° relative to the lipid bilayer normal at the peptide concentration of 0.46 μM and 68° to 90° at the peptide concentration of 1.6 μM. The similar orientation and SFG spectral intensities indicate that the adsorption amounts of LL-37 at 0.46 μM and 1.6 μM are similar.

Table 1 Fitting parameters of the SFG amide I spectra of LL-37 associated with (a) POPC/POPC, (b) POPG/POPG, (c) POPC:POPG = 7:3 and (d) POPC:POPG = 3:7 bilayers
Figure 3
figure 3

(a) SFG amide I spectra of LL-37 associated with POPC/POPC bilayer; (b) SFG O-H/N-H stretching signals collected from the POPC (top) and POPG (bottom) bilayers in contact with LL-37 solution (1.6 μM); SFG amide I spectra of LL-37 associated with (c) POPG/POPG bilayer, (d) 3:7 POPC:POPG and (e) 7:3 POPC:POPG lipid bilayers.

SFG signals detected in the frequency range of 3000–3600 cm−1 before and after the addition of LL-37 to the POPC/POPC bilayer subphase are shown in Figure 3b. Since the POPC headgroup is neutral, the POPC/POPC bilayer associated water molecules do not have a preferred orientation, therefore no water O-H stretching SFG signal was detected prior to the addition of the LL-37 stock solution into the subphase. After the addition of LL-37, a broad SFG signal in this frequency range was detected, which is mostly attributed to the water O-H stretching mode and the N-H stretching peak (centered at 3280 cm−1) of LL-37. Because LL-37 is positively charged, the ordered LL-37 molecules associated with the POPC/POPC bilayers can induce ordering of surrounding water molecules that generate the SFG O-H stretching signal. Since the POPC/POPC bilayer is neutral, the observed SFG signal in this range suggests that LL-37 potentially interacts with the POPC/POPC bilayer by hydrophobic interaction instead of electrostatic interaction.

The membrane surface orientation of LL-37 associated with the POPC bilayer deduced by the SFG amide I signal in this study is in excellent agreement with solid-state NMR results8,14.

The amide I signals observed for LL-37 associated with POPG/POPG lipid bilayers (Figure 3c) feature a single peak at ~1649 cm−1, similar to the POPC/POPC case. The SFG intensities for the high and low peptide concentrations are similar in both polarization combinations in a POPC/POPC system. However, there is a significant increase in the SFG intensity when the concentration is increased in the POPG/POPG system. The spectra were fitted with parameters displayed in Table 1b and the ratios of χeff,ppp and χeff,ssp were used to deduce the orientation angles. Since previous NMR studies indicated that LL-37 has a bend in the middle when associated with negatively charged vesicles, we adopt the “bend model” discussed above to examine LL-37's orientation associated with POPG/POPG lipid bilayers (Figure 2b). This model considers the phase continuity between the vibrational modes of the two segments in the data analysis; details can be found in the supporting information. Our analysis shows that LL-37 orients with a tilt angle of 44° to 58° relative to the POPG/POPG surface normal for the 0.46 μM peptide concentration case and 62° to 90° for the 1.6 μM concentration case. These results reveal that the membrane orientation of LL-37 in POPG/POPG bilayer is dependent on the peptide concentration and, at low concentrations, the peptides tilt to form a transmembrane orientation. This change in the peptide orientation could be attributed to the oligomerization of the peptide as observed from solid-state NMR experiments13.

SFG signals collected in the frequency range of 3000–3600 cm−1 before and after the addition of LL-37 (with the final concentration of 1.6 μM) to the POPG/POPG bilayer subphase are shown in Figure 3b. Before the addition of the LL-37, the spectra showed two broad peaks at ~3200 cm−1 and ~3500 cm−1 which can be explained as follows: the negatively charged lipid headgroups of the POPG/POPG bilayer facilitate the ordering of the associated water molecules, resulting in prominent water O-H stretching SFG signals. However, after LL-37 was added to the subphase, the two broad peaks diminished and a new peak centered at 3300 cm−1 appeared (Figure 3b). With LL-37 added, we expect the cationic peptide molecules associated with the negatively-charged POPG/POPG bilayer partially neutralized the bilayer charge. The originally ordered water molecules induced by the charged lipids would be less ordered, leading to a substantial decrease or even disappearance of the SFG water signal as shown in Figure 3b. The new 3300 cm−1 is likely from the N-H stretching mode of the peptide molecules due to the existence of ordered LL-37 molecules on the membrane surface.

In order to better simulate the real cell membrane which contains mixed lipids, we also investigated LL-37 interacting with mixed lipid bilayers with different negatively charged and zwitterionic lipid ratios. For a POPC:POPG = 3:7 bilayer, the interaction result with LL-37 has a similar trend as that with the pure POPG system (Figure 3d). After increasing the LL-37 peptide concentration from 0.46 μM to 1.6 μM, the peptide orientation changed from perpendicular to the membrane surface (0~15° vs. the surface normal) to parallel to the surface (70~90° vs. the surface normal) if we assume the peptide adopts the same structure as in negative charged SDS vesicles. This shows that when LL-37 interacts with the mixed bilayer, LL-37 molecules target the POPG component. The peptides inserted into the bilayer at low concentrations and previous research has shown that LL-37 has a tendency to oligomerize8. It is likely the inserted peptides can be pulled out from the bilayer by other peptide molecules at a higher peptide concentration via hydrophobic-hydrophobic interaction.

For a POPC:POPG = 7:3 lipid bilayer, the orientation information deduced is somewhat complicated (Figure 3e). When the peptide concentration is 0.46 μM, the peptide molecules show a parallel orientation to the surface (65~75° relative to the bilayer normal) if we assume that the peptide adopts the same structure as in the neutral charged DPC vesicles. Interestingly, at 1.6 μM peptide concentration, the peptide molecules have a transmembrane orientation (0~15° relative to the bilayer normal). After increasing the peptide concentration to 4.8 μM, 6.4 μM and 7.9 μM, the peptide molecules resumed the parallel orientation, nearly lying down on the bilayer surface.

We further investigated interactions between LL-37 and two lipid bilayers containing cholesterol (CHO). 1:1 POPC:CHO lipid bilayer was used and we monitored the SFG signal in the amide I and O-H stretching frequency regions in the same way as in the above studies without cholesterol. The experimental results showed that after adding LL-37 peptides to the subphase, the SFG water signal increases (Figure 4a), but it is not as substantial as that observed from the pure POPC lipid bilayer system (Figure 3b). The SFG water signal is induced by the positive charge of the adsorbed LL-37 molecules and therefore it is reasonable to believe it is related to the amount of adsorbed LL-37 molecules regardless of the ordering of these LL-37 molecules. With the presence of cholesterol in the lipid bilayer, fewer LL-37 molecules interact with the lipid bilayer. In addition to the much weaker water signal, there is a striking difference in the spectral feature since there is no 3300 cm−1 peak in the cholesterol-containing system. No amide I signal was detected in ssp or ppp spectrum (Figure 4b). Although we believe that the POPC/CHO system has fewer LL37 molecules adsorbed, we could not assess the ordering of the adsorbed LL-37 molecules due to the low SFG amide I signal we detected.

Figure 4
figure 4

(a) SFG O-H/N-H stretching signals collected from the 1:1 POPC:CHO (top) and 0.3:0.7:1 POPG:POPC:CHO (bottom) bilayers in contact with LL-37 solution (1.6 μM); (b) SFG amide I spectra of LL-37 associated with 1:1 POPC:CHO (top) and 0.3:0.7:1 POPG:POPC:CHO (bottom) bilayers in contact with LL-37 solution (1.6 μM).

We also investigated the incorporation of cholesterol into a lipid bilayer containing an anionic lipid: POPG:POPC:CHO = 0.3:0.7:1. After the injection of the LL-37 peptide solution into the subphase, neither discernible SFG N-H stretching signal (Figure 4a) nor SFG amide I ssp signal (Figure 4b) was detected. However, a small ppp SFG signal for the amide I band was observed. Compared to the 1:1 POPC:CHO lipid bilayer, this shows that the addition of negatively charged POPG to the lipid bilayer increases the surface coverage of associated LL-37 peptides. Still, the LL-37 surface coverage is significantly lower than that of LL-37 associated with the 3:7 POPG:POPC bilayer without CHO, suggesting that cholesterol suppresses the interaction of LL-37 with lipid bilayers. These results are in excellent agreement with solid-state NMR studies14,36.

Discussion

We compared the relative adsorption amounts of LL-37 molecules on POPC and POPG bilayers according to the observed SFG signal intensities, as summarized in Table 1. The number of adsorbed molecules is proportional to , where is deduced from the experimental SFG ssp spectrum (related to signal intensity) and is the microscopic hyperpolarizability component of the amide I signal of LL-37. The details of calculating the relative adsorption amounts can be found in the supporting information. The amount of LL-37 adsorption on the POPG bilayer is comparable to that on the POPC bilayer at 0.46 μM but it is ~2.6 times larger compared to that on the POPC bilayer at 1.6 μM (supporting information). The adsorption amount of LL-37 interacting with POPG/POPG increased as the peptide concentration was raised, which can be attributed to the electrostatic interaction between the anionic lipid and the cationic peptide. The different adsorption amounts of LL-37 on POPG and POPC lipid bilayers were confirmed by surface plasmon resonance (SPR) experiments (Figure 5). It was shown by SPR that at the concentration of 1.6 μM, the initial adsorption amount of LL-37 on POPG is ~3 times larger than that on POPC. Also, the peptides experience a certain degree of desorption on POPC but the adsorption amount remains stable on POPG. This indicates a weaker interaction between LL-37 and POPC compared to POPG. Since POPC is neutral, LL-37 only interacts with POPC via hydrophobic interactions. However, POPG is negatively charged and electrostatic interactions between the positively charged peptides and the membrane induce more peptide adsorption when the peptide concentration increases.

Figure 5
figure 5

Time-dependent SPR signals observed before and after LL-37 peptide solutions with 1.6 μM concentration in water were injected at 50 μl/min into the flow chambers to interact with the POPG (blue) and POPC (green) bilayers.

For the POPG bilayer, the adsorption curve will reach a plateau at around 220 s. While for the POPC bilayer, peptides start to desorb from the lipid bilayer after 150 s, indicating a weaker interaction between LL-37 and the POPC bilayer compared to that with the POPG lipid bilayer.

The differences in the membrane orientation of LL-37 on POPC and POPG lipid bilayers facilitate the understanding of its membrane interaction mode. It was reported that lipid headgroup perturbation induced by LL-37 is larger in bilayers containing a negatively charged lipid (DMPG:DMPC = 4:1) than in the zwitterionic DMPC bilayer, but the LL-37 molecules adopt the same parallel orientation when associated with the lipid bilayers deduced using NMR spectroscopy37. It was proposed that LL-37 exhibits a non-pore formation carpet mode on zwitterionic DPC vesicles11. Our results for LL-37 associated with the POPC/POPC and POPG/POPG lipid bilayers with a high peptide concentration of 1.6 μM supports the above conclusions: for both cases, the peptides adopt an approximately parallel orientation. As we discussed above, the unique advantage of SFG study is that SFG can be used to study the interactions between LL-37 and model cell membranes at much lower (physiologically relevant) peptide concentrations. At a lower peptide concentration of 0.46 μM, SFG results indicate that the LL-37 molecules associated with the POPG/POPG bilayer tilted towards the membrane normal. Previous LL-37 research revealed that at ~0.5 and 1.5 μM, LL-37 molecules have both monomer and dimer forms in aqueous solution38. Trimers were detected when peptide concentration is increased to 50 μM. At our experimental condition, i.e. 0.46 μM and 1.6 μM, which is similar to 0.5 and 1.5 μM, it is most likely that LL-37 dissolves in water as monomers or/and dimers. We believe that for the POPC/POPC bilayer, the peptide molecules saturated at a low concentration of 0.46 μM as monomers/dimers and remain so at a higher peptide concentration of 1.6 μM. However, for the POPG/POPG bilayer, the peptides penetrate into the membrane at a low concentration. While at a higher peptide concentration, the domination of the peptide-peptide interaction could drag the initially inserted peptide molecules out of the membrane, which not only changes the overall membrane orientation of the peptides, but also induces the increase of the peptide adsorption amount.

For the systems with mixed lipids, we believe that the concentration-dependent behavior can be explained as follows: Initially, at a low peptide concentration, most of the LL-37 molecules interact with the POPC lipid, which is the major component of the POPC:POPG = 7:3 bilayer. Therefore, the overall orientation of LL-37 is similar to that for the associated LL-37 with a pure POPC bilayer. At a slightly higher peptide concentration, more peptides can interact with the POPG lipids, which is the minor component in the bilayer. In this case, the interaction between LL-37 with the POPC:POPG = 7:3 bilayer is similar to the situation when a pure POPG bilayer interacts with a lower concentration of LL-37 in which LL-37 can insert into the bilayer. At higher peptide concentrations, additional peptides that associated with the bilayer pulled the inserted LL-37 out of the membrane, as in the pure POPG bilayer case which leads to a parallel orientation.

In summary, we have examined the molecular interaction of the LL-37 peptide with a variety of lipid bilayers using SFG (Figure 6) and have developed a SFG orientation analysis methodology for bent and disrupted α-helices. We have demonstrated that SFG is sensitive enough to study peptide-lipid molecular interaction at low-peptide concentrations, which is beyond other techniques such as NMR. LL-37 is shown to saturate the pure POPC lipid bilayer at a low concentration (0.46 μM) with an orientation parallel to the membrane surface. However, in pure POPG or POPC/POPG mixed lipid bilayers, LL-37 exhibits a reorientation upon changing the peptide concentration, suggesting the peptide aggregation process. In cholesterol-containing systems, SFG results demonstrate that cholesterol has a significant suppression effect on the peptide-membrane interaction. We strongly believe that the experimental and data analysis approaches developed in this study would be highly applicable in studying other membrane active systems including other AMPs, cell penetrating peptides, fusion peptides and amyloid proteins.

Figure 6
figure 6

Schematics showing interactions between LL-37 and different lipid bilayers.

(1) POPC bilayer at a low concentration (left, 1a) and a high concentration (right, 1b); (2) POPG bilayer at a low concentration (left, 2a) and a high concentration (right, 2b); (3) 3:7 POPC:POPG lipid bilayer at a low concentration (left, 3a) and a high concentration (right, 3b); (4) 7:3 POPC:POPG lipid bilayer at a low concentration (left, 4a) and a high concentration (right, 4b) and even higher concentrations (bottom, 4c); (5) 1:1 POPC:CHO lipid bilayer at a high concentration; (6) 0.3:0.7:1 POPG:POPC:CHO lipid bilayer at a high concentration.

Methods

Detailed calculation and experimental details can be found in the supporting information. For most lipid bilayer systems studied in this study, two final LL-37 peptide concentrations of 0.46 μM and 1.6 μM were used. The SFG spectra were collected with ssp and ppp polarization combinations after adding the LL-37 stock solution into the subphase of the model cell membrane bilayer for 1.5 hours. The SFG signal intensity remained unchanged for the next two hours, suggesting that the equilibration of the peptide-lipid bilayer interaction was reached. All the SFG spectra exhibited a single peak at ~1650 cm−1 (Figure 3). The spectral fitting results as well as the ratios of deduced from the ppp and ssp SFG spectra after deconvoluting the Fresnel factors are shown in Table 1.