Next Article in Journal
New Technological Approaches for Dental Caries Treatment: From Liquid Crystalline Systems to Nanocarriers
Next Article in Special Issue
Adjusting Heterodimeric Coiled-Coils (K/E Zipper) to Connect Autophagy-Inducing Peptide with Cell-Penetrating Peptide
Previous Article in Journal
3D-Printing of Silk Nanofibrils Reinforced Alginate for Soft Tissue Engineering
Previous Article in Special Issue
Mechanisms of a Mycobacterium tuberculosis Active Peptide
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Pharmaceutics
Volume 15
Issue 3
10.3390/pharmaceutics15030761
pharmaceutics-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Mechanism of Action of SAAP-148 Antimicrobial Peptide as Studied with NMR and Molecular Dynamics Simulations

by
Morgane Adélaïde
1,
Evgeniy Salnikov
2,
Francisco Ramos-Martín
1,*,
Christopher Aisenbrey
2,
Catherine Sarazin
1,
Burkhard Bechinger
2 and
Nicola D’Amelio
1,*
1
Unité de Génie Enzymatique et Cellulaire UMR 7025 CNRS, Université de Picardie Jules Verne, 80039 Amiens, France
2
Institut de Chimie, UMR7177, Université de Strasbourg/CNRS, 67000 Strasbourg, France
*
Authors to whom correspondence should be addressed.
Pharmaceutics 2023, 15(3), 761; https://doi.org/10.3390/pharmaceutics15030761
Submission received: 28 January 2023 / Revised: 18 February 2023 / Accepted: 21 February 2023 / Published: 24 February 2023
(This article belongs to the Special Issue State of the Art of Membrane Active Peptides)
Browse Figures
Review Reports Versions Notes

Abstract

:
Background: SAAP-148 is an antimicrobial peptide derived from LL-37. It exhibits excellent activity against drug-resistant bacteria and biofilms while resisting degradation in physiological conditions. Despite its optimal pharmacological properties, its mechanism of action at the molecular level has not been explored. Methods: The structural properties of SAAP-148 and its interaction with phospholipid membranes mimicking mammalian and bacterial cells were studied using liquid and solid-state NMR spectroscopy as well as molecular dynamics simulations. Results: SAAP-148 is partially structured in solution and stabilizes its helical conformation when interacting with DPC micelles. The orientation of the helix within the micelles was defined by paramagnetic relaxation enhancements and found similar to that obtained using solid-state NMR, where the tilt and pitch angles were determined based on 15N chemical shift in oriented models of bacterial membranes (POPE/POPG). Molecular dynamic simulations revealed that SAAP-148 approaches the bacterial membrane by forming salt bridges between lysine and arginine residues and lipid phosphate groups while interacting minimally with mammalian models containing POPC and cholesterol. Conclusions: SAAP-148 stabilizes its helical fold onto bacterial-like membranes, placing its helix axis almost perpendicular to the surface normal, thus probably acting by a carpet-like mechanism on the bacterial membrane rather than forming well-defined pores.

1. Introduction

The development of new antimicrobial agents is one of the main challenges of our times in view of the increasing threat of bacterial resistance. Multidrug-resistant (MDR) organisms expose a constantly increasing number of patients to untreatable infections [1,2,3]. Recalcitrant bacterial strains are often caused by nosocomial infections, many of them involving the so-called ESKAPE bacteria: Enterococcus faecium, Staphylococcus aureus, Klebsiella pneumoniae, Acinetobacter baumannii, Pseudomonas aeruginosa, and Enterobacter spp., some of which are able to form dangerous biofilms [3,4,5,6].
Antimicrobial peptides (AMPs) are promising molecules in the quest for new antibiotics. Produced by several organisms in nature, AMPs are relatively short (10–30 residues) and exhibit a wide range of antimicrobial or immunomodulating activities [7,8,9]. In contrast to standard antibiotics, many AMPs are able to quickly permeate cell membranes causing irreversible damage [7,10,11]. For this reason, they can be used in synergy with antibiotics to grant intracellular access to the latter, thus overcoming resistance based on reduced permeabilization of the bacterial membrane [12,13]. Some AMPs are able to bypass the bacterial membrane without damaging it and exert their antimicrobial action intracellularly [7,10,11,14,15,16,17,18]. While exceptions exist [19], AMPs are less prone to resistance due to multiple factors. First, their fast-killing rate (compared to antibiotics) does not favor bacterial replication, possibly leading to new resistance [7,20,21]. Second, the heterogeneity of their target (bacterial membrane) makes the interaction less specific, thus making it harder for bacteria to completely impair their action by single-point mutations [22]. Third, some AMPs act by multiple modes of action [23,24,25].
One of the main limitations of AMPs is their degradation or loss of activity in blood and tissues [3,26]. Several strategies have been used to overcome these problems such as the substitution of L-amino acids with their D-isomers or the introduction of non-natural amino acids [27,28,29,30]. Other strategies are based on the study of a huge portfolio of peptides designed from an existing one, aiming to optimize its activity and reduce its degradation properties. One of these peptides is SAAP-148 (Synthetic Antimicrobial and Antibiofilm Peptide), which derives from LL-37 and shows an excellent activity profile against drug-resistant bacteria and biofilms while resisting degradation when tested in physiological conditions [31,32,33]. Together with SAAP-148, another 25 peptides were inspired by LL-37, all designed using the substitution of anionic glutamine with cationic arginines or lysines, thus enhancing their cationicity and helicity in order to increase their antimicrobial activities. When tested against Staphylococcus aureus, SAAP-148 performed better than others with an LC99.9 of 1.6 μM in PBS and 12.8 μΜ in the presence of plasma in the medium [31]. Similar values were obtained using other resistant strains such as Pseudomonas aeruginosa, Enterococcus faecium, Klebsiella pneumoniae, Acinetobacter baumannii [31], or Enterococcus hirae [33]. Despite its optimal pharmacological properties (without remarkable toxic effects) [31,32,33], its mechanism of action remains poorly understood at the molecular level. The aim of this work is to elucidate its interaction with bacterial and mammalian membranes at the molecular level of detail.
The interaction of AMPs with membranes has been described in detail using the SMART model [34], which shows how amphipathic peptides tend to align to the surface and intercalate among the lipid headgroups, while apolar peptides tend to insert and span the thickness of the bilayer. The model suggests that bilayers adapt to membrane-inserted peptides until threshold concentrations are reached. In the presence of small headgroups, such as ethanolamine of PE, the amphipathic peptide may even stabilize the bilayer properties by intercalation. However, at high peptide-to-lipid ratios, the peptide tends to cover the surface (carpet model) causing the formation of transient micro-sized pores [35,36,37,38,39,40,41,42,43,44,45] that allow the transmembrane passage of molecules [11,46,47]. This is possible because of the “soft” nature of the bilayer which can deform, change its thickness, and adapt to different perturbations [48,49]. Depending on the nature of the peptide and the conditions (pH, buffer, temperature, etc.), different structures can be formed including barrel-stave or toroidal (also called “worm-hole”) pores, induction of non-lamellar phases, non-lytic depolarization, localized thinning bicelle formation, and detergent-type micellization which efficiently dissolve the membrane [7,11,35,48].

2. Materials and Methods

2.1. Peptide Synthesis

The peptide SAAP-148 (Ac-LKRVWKRVFKLLKRYWRQLKKPVR-NH2) was synthesized following the standard Fmoc-synthesis protocol using a Millipore 9050 automatic peptide synthesizer. The carboxyl amide C-terminus was archived using TentaGel S Ram resin (Rapp Polymeres, Tübingen Germany). Acetylation of the N-terminus was performed using 2% acetic anhydride. The leucine at position 11 or at position 12, one at a time, were 15N labeled using labeled Fmoc Leucine (Cortecnet, Voisins le Bretonneux, France). The peptides were purified with reverse phase HPLC (Gilson, Villiers-le-bel, France) using a preparative C18 column (Luna, C-18-100 Å-5 mm, Phenomenex, Le Pecq, France) and an acetonitrile/water gradient (Figure S1). Their identity and purity (>90%) were checked using MALDI mass spectrometry (MALDI-TOF Autoflex, Bruker Daltonics, Bremen, Germany) (Figure S2).

2.2. Liquid-State NMR Spectroscopy and Sample Preparation

2.2.1. Titration of SAAP-148 with DPC Micelles

SAAP-148 was solubilized at a concentration of 2.4 mM in 500 µL of 50 mM phosphate buffer at pH 6.6 containing 10% of D2O and 0.1 mM deuterated sodium 3-(trimethylsilyl)propionate:d4 (TSP:d4) as an internal reference. A few microliters of a 1 M solution of dodecylphosphocholine-d38 (DPC:d38) were gradually added to obtain a final concentration of 60 mM at 310 K.

2.2.2. Paramagnetic NMR Experiment

A measured amount of 16-doxyl stearic acid (16-DSA) was dissolved in deuterated methanol and added to samples of SAAP-148 in DPC micelles in microliter amounts until the paramagnetic effect was observed as a wider linewidth on selected peaks (typically in the ratio 16-DSA:DPC 1:10). Paramagnetic contributions were evaluated as percentage loss in the intensity of Hα/Cα cross peaks in 1H,13C-HMQC spectra. Errors were evaluated based on spectral noise using the relation for the propagation in ratios A = B C : ΔA = A( Δ B B + Δ C C ), where B is the cross-peak intensity without 16-DSA, C is the cross-peak intensity with 16-DSA, and ΔB and ΔC are the spectral noise. Experiments were run at 310 K.

2.2.3. Titration of SAAP-148 with Bicelles

SAAP-148 was solubilized in 500 µL of 50 mM phosphate buffer at pH 6.6 containing 10% of D2O and 0.1 mM TSP-d4 as an internal reference. Samples were typically at a concentration of 0.8 mM. A 1 M solution of isotropic bicelles was obtained by solubilizing in chloroform 1,2-dihexanoyl-sn-glycero-3-phosphocholine (DHPC), 1,2-dimyristoyl-sn-glycero-3- phosphocholine (DMPC), and 1,2-dimyristoyl-sn-glycero-3-phospho-(1′-rac-glycerol) (DMPG) at molar ratios 66.7:25.0:8.3. Few microliters of a 1 M solution of DHPC/DMPC/DMPG bicelles were gradually added to obtain a final concentration of 50 mM at 310 K.

2.2.4. NMR Acquisition and Processing

Liquid state NMR experiments were recorded using a 500 MHz Bruker spectrometer equipped with a 5 mm Broadband Inverse (BBI) probe. Backbone resonance assignments were achieved by 1H,13C-HMQC, 1H,1H-TOCSY (mixing time of 60 ms), and 1H,1H-NOESY (mixing time of 200 ms). Spectra were analyzed with TopSpin 4 (Bruker BioSpin). Secondary structure prediction was achieved by observing deviations from random coil values obtained with the POTENCI web server accessed on 21 October 2021 (https://st-protein02.chem.au.dk/potenci/) [50] at 310 K.

2.3. Solid-State NMR Spectroscopy and Sample Preparation

The sample preparation and NMR measurements here are described in detail and illustrated in [51]. In short, oriented samples for solid-state NMR were prepared by co-dissolving 3.25 mg of SAAP-148 powder in methanol/chloroform (1:1 by volume) with corresponding amounts of POPE (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanolamine) and POPG (1-palmitoyl-2-oleoyl-sn-glycero-3-phospho-(1′-rac-glycerol)) lipids to reach a 2:75:25 molar ratio. The sample solution was vortexed and sonicated using a bath sonicator for 5 min and the solvent was partially evaporated under a stream of nitrogen gas. The viscous mixture was spread on 22 ultrathin cover glass plates (8 × 12 mm2, Marienfeld, Lauda-Königshofen, Germany), and the samples were dried under vacuum overnight to remove the organic solvent. Thereafter, the glass plates were incubated at 96% relative humidity at 37 °C (i.e., well above the gel-to-fluid transition temperature) for two days, stacked on top of each other, tightly sealed and introduced in the NMR spectrometer with the glass plates normal being parallel to external magnetic field direction.
The preparation of non-oriented samples started by mixing POPE/POPG-d31 (3:1, 4 mg/1.5 mg) or POPE/POPE-d31/POPG (2:1:1, 2.9 mg/1.5 mg/1.5 mg) and the appropriate amount of peptide to reach a 2% mole ratio in methanol/chloroform 1/1 (v/v). The sample solution was vortexed and sonicated using a bath sonicator for 5 min and the solvent was evaporated using exposure to a stream of nitrogen and to a high vacuum overnight in such a manner to form a film along the walls of a glass tube. The sample was resuspended in 23.5 μL (water content h = 0.81) of 100 mM Tris buffer (pH 7.4) and involved vortexing and bath sonication, as well as 3 freeze/heat cycles at −20 °C and 40 °C. The glass tube (6 mm outer diameter) with the sample was inserted into the solenoidal coil of the solid-state NMR probe.
Proton-decoupled 31P solid-state NMR spectra were acquired at 121.575 MHz using a Bruker Avance wide-bore 300 solid-state NMR spectrometer equipped with a commercial double-resonance flat-coil probe (Bruker, Rheinstetten, Germany) [52]. A Hahn-echo pulse sequence [53] was used with a π/2 pulse of 5 µs, a spectral width of 100 kHz, an echo delay of 100 µs, an acquisition time of 10.2 ms, and a recycle delay of 3 s. External 85% H3PO4 at 0 ppm was used for calibration. The temperature was set to 310 K.
The proton-decoupled 15N cross-polarization spectra of static aligned samples were acquired at 30.43 MHz using a Bruker Avance wide bore 300 NMR spectrometer. An adiabatic CP pulse sequence was used with a spectral width, acquisition time, CP contact time, and recycle delay times of 25 kHz, 10.2 ms, 0.4 ms, and 3 s, respectively [51,54]. The 1H π/2 pulse and spinal-64 heteronuclear decoupling field strengths B1 corresponded to a nutation frequency of 31 kHz. Then, 40 k scans were accumulated, and the spectra were zero-filled to 4 k points. An exponential line broadening of 50 Hz was applied before the Fourier transformation. Spectra were externally referenced to 15NH4Cl at 39.3 ppm [55]. Samples were cooled with a stream of air at a temperature of 310 K.
Next, 2H solid-state NMR spectra of deuterated lipid samples were acquired using a quadrupole pulse-echo sequence [56] with a recycle delay of 0.3 s, an echo time of 100 μs, a dwell time of 0.5 μs, and a π/2 pulse of 5 µs. Before the Fourier transformation of the free induction decay, an exponential apodization with a line broadening of 100 Hz was applied. Spectra were externally referenced to 2H2O (0 Hz). The temperature was set to 310 K.
The deuterium order parameters (SCD) for each CD2 and CD3 group were determined according to S C D i = 4 3   h e 2 q Q Δ i υ , where Δiυ is the quadrupolar splitting of segment i and (e2qQ/h) is the static quadrupole coupling constant (167 kHz) observed for deuterons within C–D bonds [57].
Orientational restraints from oriented solid-state NMR spectra. A coordinate system was defined where the tilt angle is defined as the angle between the helix long axis and the membrane normal. The α-helical conformation had Ramachandran angles of (ϕ = −65°, ψ = −45°). The tilt and pitch angles were scanned, where at each orientation, the corresponding 15N chemical shift was calculated. The standard deviation (SD) of a Gaussian line shape takes into account orientational heterogeneity. Independent wobbling (10° Gaussian distribution) and azimuthal fluctuations around the helix long axis (SD 18°) were taken into consideration by averaging the resonance values on the ensemble of orientations with the corresponding Gaussian distributions, similar to the dynamics of surface-bound sequences that have been tested previously including N-terminal 17 residues domain of Huntingtin protein [58].

2.4. Circular Dichroism

The peptides were dissolved at 50 µM in a 10 mM phosphate buffer at pH 7 in the presence and in the absence of NaCl 100 mM or Tris buffer at pH 7. CD spectra were recorded at 37 °C using a J-810 spectropolarimeter (Jasco, Tokyo, Japan) using 2 scans, a step resolution of 1 nm, and an adaptive integration time between 1 and 8 s. The path length of the quartz cuvette was 1 mm. The spectra were recorded between λ = 190 and 250 nm.

2.5. Molecular Dynamics Simulations

The CHARMM-GUI web server [59,60,61] was used to prepare the bilayer systems for MD simulations. GROMACS software [62] was used for the calculations with the CHARMM36m force field [63]. Water molecules were described with the TIP3P model [64]. Bilayers contained 64 lipids in each leaflet, and peptide molecules (whose initial structure was calculated using I-Tasser [65]) were placed over the upper leaflet at a distance >10 Å, avoiding intermolecular contacts in the setups containing more than one molecule of the peptide. The C-terminus was amidated with the CHARMM terminal group functionality. Lysine and arginine side chains were protonated. A 50 Å thick water layer was added above and below the membrane resulting in about 15,000 water molecules (30,000 in the case of pure cardiolipin (CL) systems). A concentration of 150 mM of NaCl, CaCl2, or MgCl2 was used when appropriate.
Periodic boundary conditions were applied for all the simulations. The particle mesh Ewald (PME) method was used for long-range electrostatic interactions. Each system was energy-minimized with a 5000-step steepest-descent algorithm. Equilibration was performed with the Berendsen barostat [66] while a Parrinello–Rahman barostat [67] was used to maintain pressure (1 bar) semi-isotropically with a time constant of 5 ps and a compressibility of 4.5 × 10–5 bar–1. A Nose–Hoover thermostat [68,69] with a time constant of 1 ps was chosen to maintain the systems at 310.15 K. All bonds were constrained using the LINear Constraint Solver (LINCS) algorithm, allowing a 2 fs integration step. Then, 500 ns simulations were run (except when a different duration is specified in the text) after the standard CHARMM-GUI minimization and equilibration steps [70]. The whole process (minimization, equilibration, and production run) was performed twice in the absence of peptides and thrice in their presence. Convergence was assessed using the analysis of RMSD and polar contacts.
The following molar proportions were used for the mixed bilayers: POPE/POPG 70/30, POPE/POPG/TOCL2 67/27/6, and POPC/CHOL 70/30. MD trajectories were analyzed using GROMACS tools. Graphs and images were created with GNUplot [71] and PyMol [72]. MOLMOL [73] and VMD [74] were used for visualization.

3. Results and Discussion

3.1. Liquid-State NMR studies of SAAP-148 Interacting with Model Membranes

3.1.1. SAAP-148 in Aqueous Solution

The 1H NMR spectrum of SAAP-148 in phosphate buffer at pH 6.6 is shown in Figure 1A (bottom spectrum). The significant broadening of signals might be due to exchange with the water solvent (for exchangeable protons) and/or conformational exchange/aggregation equilibria that occur on times scales equivalent to the acquisition time (micro- to milliseconds). The latter hypothesis is supported by the fact that broadening is observed also in the aromatic and aliphatic region that does not contain exchangeable protons. The concentration dependence of the 1H NMR spectrum (Figure S3A) strongly suggests that SAAP-148 is aggregating or rather forming a multimer at concentrations larger than about 100 µM. In helical conformation, this peptide would form a perfectly amphipathic helix, suggesting that it might dimerize to avoid exposure of the apolar side chain to the water. At concentrations lower than 100 µM, the NMR lines are much narrower and CD spectra (Figure S3B) show the characteristic pattern of a random coil (minimum at about 200 nm).
Even though extreme broadening in the amide region prevents resonance assignment by 1H,1H-TOCSY and 1H,1H-NOESY at the concentration used for NMR, the signal dispersion appears compatible with residual structuring. A value of 0.88 ppm is in fact intermediate between a predicted dispersion of 0.48 ppm for a random coil (estimated usingPOTENCI web server [50], Table S1) and that of a full alpha-helical structure (1.05 ppm as estimated using SHIFTX2 web server [75], Table S1). However, the absence of HN-HN cross peaks in the NOESY spectrum indicates that such structures, if helical, are not well defined in solution.

3.1.2. SAAP-148 Interaction with Biomimetic Membrane Models

The addition of a concentrated solution of DPC-micelles significantly modifies the spectrum of SAAP-148. Many peaks re-emerge from the baseline at a protein:detergent molar ratio of 1:25 (see, for example, the two HN protons of the tryptophan side chain at about 10 ppm in Figure 1A), and the final spectrum in micelles become relatively well resolved to allow resonance assignment (Figure 1B). The gradual shift and broadening of SAAP-148 peaks as DPC micelles are added to the solution are indicative of fast exchange close to coalescence on the NMR time scale.
The main limitation of DPC resides in the fact that micelles do not form bilayers and that access to the lipid core is facilitated by their high curvature [76]. A better membrane model is provided with bicelles constituted by a mixture of DHPC and DMPC where the bilayer formed by DMPC is stabilized by the short chains of DHPC, which prevent direct exposure of the DMPC acyl chains to water [77]. Furthermore, these systems can incorporate to some extent different headgroups for modeling various types of biological membranes (e.g., PG for bacteria [78,79]). For solution-state NMR, isotropic bicelles are ideal for their fast-tumbling properties [77]. For this reason, we titrated a solution of SAAP-148 with a concentrated solution of isotropic bicelles constituted by a mixture of DHPC/DMPC/DMPG at a ratio of 2:0.25:0.75 (Figure 1C). The comparison of Figure 1A,C suggests that SAAP-148 interacts in a similar fashion with micelles and isotropic bicelles. Unfortunately, the size of the complex results in a shortening of the transversal relaxation time, which prevents the observation of Hα/Cα resonances in the 1H,13C-HSQC (or 1H,13C-HMQC) spectra and consequently the assignment of the bound form of SAAP-148. Interestingly, Figure 1C shows that the binding process takes place with a slower exchange than with DPC micelles (peaks do not shift gradually but tend to reappear at different positions along the titration), indicating a stronger interaction.
The full 1H backbone assignment was achieved for SAAP-148 in DPC micelles using the 2D 1H,1H-TOCSY, 1H,1H-NOESY, and 1H,13C-HMQC spectra (Figure 1B and Figure S4 and Table S1), allowing the calculation of the Hα and Cα resonance deviations from random coil values (Figure 2A,B). Persistent negative (Hα) and positive (Cα) deviations beyond the threshold limit (0.1 ppm for 1H and 0.7 ppm for 13C) [80] indicate the formation of a stable alpha helix interrupted only by proline 22 at the C-terminus [81,82,83].
By introducing a paramagnetic species in the system, the formation of a stable alpha-helix and its orientation in the micelles can be determined. Paramagnetic probes are commonly used in NMR measurements for their strong peak-broadening effect, which depends on the sixth power of the distance [84,85,86]. They can serve as a reference point to determine the location of atoms in the species under study. In our case, we added 16-doxyl stearate (16-DSA) at a molar ratio of 16-DSA:DPC 1:10. This is a stearic acid bearing a paramagnetic spin label at the end of its acyl chain. This location allows placing the paramagnetic probe in the core of the micelles once the molecule gets inserted [84]. Paramagnetic effects were quantified by the percentage loss in the intensity of Hα/Cα cross peaks in the 1H,13C-HMQC spectra. Figure 2C shows that the paramagnetic effect approximately follows the periodicity of a helix (3.6 amino acid per turn), thus confirming the results obtained with the analysis of chemical shift deviations (Figure 2A,B). Interestingly, the maxima in Figure 2C have comparable intensity along the sequence, an effect compatible with an orientation of the helix axis parallel to the membrane surface (or to the tangent touching it, given the high curvature of the micelle). A more pronounced paramagnetic effect at the N-terminus suggests its slightly deeper insertion in the micelle. Furthermore, positively charged Lys and Arg residues are clearly less affected than the others, indicating that the amphipathic helix (see insert in Figure 2C) exposes its apolar face towards the center of the micelle. Assuming that lysine side chains interact with phosphate groups of phospholipids (as often observed for antimicrobial peptides [87]), this implies that the peptide is well inserted in the micelle with its apolar face making important van der Waals contact with lipid chains.

3.2. Solid-State NMR Studies on SAAP-148 Peptide Interacting with Model Membranes

3.2.1. The Orientation of SAAP-148 in Bacterial Biomimetic Membranes

While NMR data in solution have provided important information on the interaction of SAAP-148 with DPC micelles, a better model of bacterial membranes can be obtained from phospholipid bilayers, which can be arranged either as liposomes or as planar structures on mechanical supports [51]. The large size of such systems requires the use of solid-state NMR spectroscopy in combination with isotopic labeling of either the peptide and/or the lipids. In order to confirm the results obtained with solution NMR spectroscopy also for more realistic systems, we have studied the interaction of SAAP-148 selectively labeled with 15N with POPE/POPG uniaxially oriented bilayers. Our choice is justified by the fact that bacterial cell membranes usually expose anionic lipids such as PG and CL and, in most gram-negative species, PE is the major phospholipid exposed to the outer leaflet of the cytoplasmic membrane [78,79].
The anisotropy of the 15N chemical shift is used to determine the pitch and tilt angles of the helix formed by SAAP-148 in oriented systems. These two parameters refer to the angles formed by the long (tilt) and short (pitch) cylindrical helix axes with the bilayer normal (which is assumed to be along the magnetic field direction in oriented systems). The homogeneity of the lipid orientation, parallel to the magnetic field, was ascertained using 31P solid-state NMR spectra (Figure 3A,B). The predominant intensity around 30 ppm is indicative that most lipids are well aligned along the surfaces of the mechanical support while some intensity extends up to −15 ppm, suggesting a global deformation of the lipid bilayer concomitant with other conformations and/or alignments of the lipid head group [89].
In a next step, peptides labeled with 15N in the central part of the peptide were investigated. This region adopts a stable helical conformation in the presence of detergent micelles (Figure 2) suggesting that Leu 11 and 12 are also part of this structured domain when interacting with lipid membranes. Indeed, when reconstituted into uniaxially oriented phospholipid bilayers the 15N spectra show well defined resonances indicative that both residues are part of a helical structure that is homogeneously aligned relative to the membrane normal (Figure 3C,D). The 15N chemical shifts are 80.0 ± 2.5 ppm and 70.0 ± 4.0 ppm, respectively, indicating a helix alignment approximately parallel to the membrane surface [90].
When analyzed in quantitative detail, these values are compatible with the combinations of tilt and pitch angles shown in Figure 3E as blue (Leu 11) and green lines (Leu 12), respectively. However, only two possible solutions (I and II) are compatible with both measurements and are highlighted in red in panel E of Figure 3. The two regions of tilt/pitch angular pairs correspond to the orientations represented in panel F. Despite the different conditions (DPC micelles vs. POPE/POPG bilayers), solution II is compatible with what was observed in the liquid state. Additionally, in that case, a slight inclination of the N-terminus towards the membrane core was observed (tilt), together with the orientation of apolar residues towards the interior of the bilayer (pitch).

3.2.2. The Effect of SAAP-148 on Biomimetic Membrane Dynamics

Once the orientation of the peptide was determined in our membrane models, the effect of SAAP-148 on the phospholipids and their order parameters was also investigated. This was achieved using 2H static solid-state NMR spectroscopy of liposomes where either POPG or POPE were 2H labeled along the palmitoyl chain (Figure 4A,B). The resulting 2H solid-state NMR spectra are characterized by superimposed quadrupolar splittings, which represent the order parameter of the various C–D bond vectors [91,92]. This implies that it is possible to calculate the order parameter of each position of the deuterated chain by measuring the corresponding quadrupolar splitting in the spectrum [91,92] (Figure 4C,D).
The presence of SAAP-148 clearly affects the spectrum of POPE/POPG liposomes, independently of the 2H labeled species (POPG in panels A, C, and E or POPE in panels B, D, and F), thus ruling out the specificity for potential inhomogeneous phospholipid domains in the membrane. When focusing on the effect on the order parameter (Figure 4B–E), SAAP-148 perturbs all the lipid palmitoyl chain (Figure 4C,D), with particular emphasis on the lipid tail (Figure 4E,F), reinforcing the hypothesis of deep insertion into the membrane interface as suggested by liquid state NMR data. Other linear cationic amphipathic peptides such as magainin, HB43, or designed model sequences have also shown considerable disordering of the fatty acyl chain packing when intercalating into the membrane interface [91,93,94,95]. SAAP-148 shows a particularly pronounced effect with an up to >30% reduction in the order parameter (Figure 4E,F). Furthermore, in these previous cases, a preferential disordering of the PG palmitoyl chain has been observed [91,93].

3.3. MD Simulations of SAAP-148 Interacting with Biomimetic Bilayers

In order to better understand the details of the interaction of SAAP-148 with its targeted membranes, we complemented our NMR results with several molecular dynamics simulations. To this end, we have simulated membrane systems containing phospholipid ratios representative of bacterial and mammalian cells [79], to distinguish between their bactericidal and possibly toxic effects, respectively.
One of the simplest models of a mammalian cell membrane consists of POPC phospholipids. Despite its simplicity, this system has been used in a wide range of studies in which it was able to successfully reproduce biological effects [16,96,97,98]. This is due to the fact that the PC moiety is the main headgroup exposed on mammalian cells, while other groups such as PE are often exposed to the inner leaflet [79,99,100,101]. In the present work, we have also considered the effect of cholesterol, another key lipid in mammalian cell membranes.
Regarding bacteria, as mentioned above, their cell membranes usually expose anionic lipids such as PG and CL [78,79] but also PE. As a consequence, we have extensively focused on PE/PG or PE/PG/CL bilayers and also on their pure components PE, PG, or CL to isolate their contribution to the interaction in the mixtures.
Since most media used for the determination of minimal inhibitory concentrations contain Na+, Mg2+, and Ca2+, simulations were performed with each of these ions one at a time. Calcium is known to affect the interaction of peptides and proteins with the bilayers by modulating the interplay between self-aggregation and membrane binding [102,103]. Calcium can either promote aggregation [102] or binding [104,105,106,107]. In the first case, Ca2+ binds to the phosphate groups of the bilayer and shields its negative charge, thus preventing the binding of positively charged residues of the peptide. Probably for the same reason, Mg2+ is known to impair the action of some AMPs at high concentrations [108]. Less commonly, Ca2+ forms salt bridges connecting charged glutamate residues of the peptide to phosphate groups of the membrane, facilitating the interaction [102,107].

3.3.1. Effect of SAAP-148 on Mammalian Model Membranes

In order to highlight the occurrence of key interatomic interactions triggering the association of SAAP-148 with membranes, we estimated the occurrence of H bonds, salt bridges, and van der Waals contacts along the trajectory (Figure 5A–C). This was quantified by calculating the radial distribution function [64] of all peptide atoms from each phospholipid atom. While polar contacts (H bonds and salt bridges) are particularly informative to study the first stages of the interaction, van der Waals contacts mostly inform on the final stages, in case the peptide has been inserted into the membrane.
While SAAP-148 is able to establish several polar contacts with POPC membranes (Figures S5–S7), its penetration in the membrane remains superficial as testified by the infrequent apolar contacts (Figures S6–S8). When the model is improved with the addition of cholesterol in the ratio POPC:CHO 70:30, the number of contacts is drastically reduced (Figures S5 and S8), and the preservation of the order parameter (Figures S9–S12) indicates that the membrane interior is not affected by the presence of the peptide (with one exception in the presence of Mg2+). Interestingly, in the presence of cholesterol, we observe a loss of the alpha-helix fold at both peptide termini, which might be due to the selectivity of SAAP-148 towards bacterial membranes (Figure S13).

3.3.2. Effect of SAAP-148 on Bacterial Model Membranes

With membrane models containing lipids characteristic of bacteria (PE, PG, CL), the occurrence of polar contacts increases by one order of magnitude as compared to mammalian models (compare Figures S5–S7 with Figures S13–S17). Arginines and lysines drive the interaction by forming salt bridges with phosphate oxygen atoms of phospholipids (Figure 5B) following a clear preference order CL > PG > PE (Figure S16), which is probably based on the increasing negative charge of the target (−2, −1, 0).
It should be noted that the insertion of SAAP-148 into bacterial models requires long simulation times (up to 6 μs) in order to observe penetration of the peptide, which is observed only in rare cases, mostly in the presence of calcium. Following the pattern observed for other AMPs [14,15,87,109,110], the peptide first sticks to the membrane from its exterior by the formation of the salt bridges described above and rigidifying the motion of acyl chains as observed by higher values of their order parameter (Figure S18). Subsequently, the peptide fluidifies the membrane (Figure S18) by deforming the surface but not necessarily being fully internalized, thus reproducing in some repetitions the experimental behavior of the order parameter (compare Figure 4C,D with Figures S11 or S18C). Although a significant reduction in the order parameter along the acyl chain of phospholipids can be obtained even without complete peptide internalization, we believe that SAAP-148 does get internalized in longer time scales (which might require enhanced sampling methods [111]). The orientation of the helix obtained by paramagnetic effects (where apolar side chains are closer to the center of micelles than charged ones) and the large effect on the order parameter of the very end of lipid tails (Figure 4E,F) suggest a deep internalization.
The significant perturbation of the surface can be monitored by changes in the electron density profile (Figures S19 and S20) and dipole potential profiles (Figure S21). The former (see Figure S19C) shows a clear inter-digitalization of the acyl chain extremities (increased density in the central part of the membrane) and a splitting of the phosphate maxima (at −2 and +2 nm from the center) in two populations due to disorder induced by the presence of the peptide and the entrance of water (whose density is not included in the graph). The latter (see Figure S21C) indicates that the presence of SAAP-148 lowers the dipole potential, thus facilitating the destabilizing transit of dipolar molecules across the bilayer [79,112,113].
Although a complete penetration of the peptide was observed only in simulations with a single peptide (Figure 5D), visual examination of simulations with multiple peptides suggests that dimerization could be an important phenomenon for the interaction (Figure 5E, Figures S22 and S23), although the dimer poorly penetrates the bilayers, at least in the time scale of our simulations (Figures S24–S26). Dimerization of SAAP-148 would not be surprising for a peptide that in helical conformation forms a perfect amphipathic helix and is perfectly consistent with what is observed experimentally in NMR spectra of the free peptide (see Section 3.1.1). Simulations show that, as long as the two hydrophobic faces shield each other from exposure to water, parallel or antiparallel alignment are both possible (Figure 5E).

4. Conclusions

SAAP-148 is an extremely promising antimicrobial peptide due to its capacity to kill MDR bacteria, prevent the formation of their biofilm, and even eradicate persister cells [31]. Furthermore, it has been shown to be resistant to degradation in physiological fluids [31]. The antimicrobial activity of SAAP-148 has been explained by its ability to interact and destabilize the bacterial membrane. This is reflected by the peptide’s ability to disrupt liposomes containing bacterial phospholipids [31,33] and potentiate the activity of other antibiotics [12]. In this work, we have elucidated its mechanism of action at the molecular level of detail. Experiments based on liquid and solid-state NMR spectroscopies show that SAAP-148 forms stable alpha-helical structures on the membrane surface and is able to become internalized in bacterial membrane mimetics. Molecular dynamics simulations require a long time to observe this internalization, indicating that the energy barrier is higher than for other AMPs [94]. Such a barrier could be lowered by the presence of metal ions such as Ca2+, although further studies are needed to verify this phenomenon.
Most importantly, using a variety of complementary approaches, we have consistently shown that the axis of the helix formed by SAAP-148 on bacterial-like membranes is almost perpendicular to the bilayer normal. Thus, SAAP-148 can act by breaking the membrane once a dense peptide carpet has been accumulated at the membrane interface, rather than forming transmembrane pores. This is consistent with the SMART model [34], in which amphipathic helices lay on the surface. At low peptide:lipid ratios, the soft bilayers can adjust to the disturbance caused by the interfacial peptide. According to the model, the bilayer packing can even be stabilized due the presence of peptide, e.g., in PE-rich membranes [114] (see the increase in the order parameter in Figure S18A,B). At higher peptide:lipid ratios, the changes in the order parameter are consistent with a disruption in the membrane order (Figure 4C–F and one repetition of simulations in Figure S18C) leading to the formation of transient pores [34]. This allows the internalization of the peptide and the subsequent formation of membrane-disassembling structures.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/pharmaceutics15030761/s1, Table S1. Chemical shift in NH, Hα, and Cα obtained experimentally (Spectra) for a random coil structure (POTENCI) or for an alpha-helix (ShiftX). Figure S1. Analytical purity of SAAP-148 peptide: HPLC chromatogram. Figure S2. Analytical purity of SAAP-148 peptide: MALDI-TOF mass (A) with zoom (B) in the relevant region. Figure S3. (A) Amino/aromatic region of 500 MHz 1H spectrum of SAAP-148 in phosphate buffer 50 mM at pH 6.6 (H2O/10%D2O) at 310 K as a function of concentration. (B) Normalized circular dichroism spectra of 50µM SAAP-148 in 10 mM Tris pH7 (black), 10 mM phosphate buffer pH7 (red), or 10mM phosphate buffer pH7 + 100mM NaCl (blue). The temperature was 310K. Figure S4. Hα/Cα region of 500 MHz 1H,13C-HSQC spectrum of SAAP-148 1.3 mM in DPC:d38 micelles. The assignment of non-glycine residues (whose signal is outside the chosen region) is shown on each peak. Figure S5. Polar contacts (hydrogen bonds and salt bridges) of SAAP-148 with different lipid bilayers. Sodium was used as a counterion. Figure S6. Polar (hydrogen bonds and salt bridges) and apolar (van der Waals interactions) contacts of SAAP-148 with POPC in the presence of calcium. Figure S7. Polar (hydrogen bonds and salt bridges) and apolar (van der Waals interactions) contacts of SAAP-148 with POPC in the presence of magnesium. Figure S8. Apolar contacts (van der Waals interactions) of SAAP-148 with different lipid bilayers. Sodium was used as a counterion. Figure S9. Order parameter for different bilayer models mimicking common bacterial and mammalian lipids with (in red) and without (in black) one molecule of SAAP-148 peptide. Sodium was used as a counterion. Figure S10. Order parameter for different bilayer models mimicking common bacterial and mammalian lipids with (in red) and without (in black) multiple SAAP-148 peptides. Sodium was used as a counterion. Figure S11. Order parameter for different bilayer models mimicking common bacterial and mammalian lipids with (in red) and without (in black) single (up) or multiple (down) SAAP-148 peptides in the presence of calcium. Figure S12. Order parameter for different bilayer models mimicking common bacterial and mammalian lipids with (in red) and without (in black) single (up) or multiple (down) SAAP-148 peptide in the presence of magnesium. Figure S13. DSSP secondary structures calculated along molecular dynamics (MD) simulations. Sodium was used as a counterion. Figure S14. Polar contacts (hydrogen bonds and salt bridges) of SAAP-148 with different lipid bilayers. Sodium was used as a counterion. Figure S15. Polar contacts (hydrogen bonds and salt bridges) of SAAP-148 with different lipid bilayers. Sodium was used as a counterion. Figure S16. Polar contacts (hydrogen bonds and salt bridges) of SAAP-148 with PE/PG (at 500 ns and 4us) and PE/PG/CL (500 ns) bilayers. Sodium was used as a counterion. Figure S17. Polar (hydrogen bonds and salt bridges) and apolar (van der Waals interactions) contacts of SAAP-148 with POPE/POPG in the presence of magnesium. Figure S18. Order parameter for POPE/POPG (70:30) bacterial mimicking bilayers with (in red) and without (in black) SAAP-148 multiple peptides. (A) Slight increase in the order parameter at 500 ns and (B) a decrease at 4 μs and (C) 6 μs. Sodium was used as a counterion. Figure S19. Electron density profiles with (in red) and without (in black) SAAP-148 peptide. (A) POPC model mimicking mammalian membranes. (B) POPE/POPG (70:30) model mimicking bacterial membranes at 500 ns and (C) 4 us. Sodium was used as a counterion. Figure S20. Electron density profiles with (in red) and without (in black) SAAP-148 peptide. For different membrane model bilayers. Sodium was used as a counterion. Figure S21. Dipole potential profiles with (in red) and without (in black) SAAP-148 peptide. For different membrane model bilayers: (A) POPC:CHO as a mammalian cell model, (B) POPE:POPG up to 500 ns, and (C) POPE:POPG up to 4us. Sodium was used as a counterion. Figure S22. Contact map (left) showing the dimerization effect of SAAP-128 (right). Sodium was used as a counterion. Figure S23. Contact maps showing the aggregation behavior in the presence of calcium (left) and magnesium (right). Figure S24. Apolar contacts (van der Waals interactions) of SAAP-148 with different lipid bilayers. Sodium was used as a counterion. Figure S25. Apolar contacts (van der Waals interactions) of SAAP-148 with different lipid bilayers. Sodium was used as a counterion. Figure S26. Apolar contacts (van der Waals interactions) of SAAP-148 with PE/PG (at 500 ns and 4 us) and PE/PG/CL (500 ns) bilayers. Sodium was used as a counterion.

Author Contributions

Conceptualization, C.S., B.B. and N.D.; methodology, E.S., F.R.-M., B.B. and N.D.; software, F.R.-M.; validation, F.R.-M., C.S., B.B. and N.D.; formal analysis, M.A., E.S., F.R.-M., and N.D.; investigation, M.A., E.S., F.R.-M. and C.A.; resources, C.S., B.B. and N.D.; data curation, M.A., E.S., F.R.-M. and C.A.; writing—original draft preparation, M.A., F.R.-M. and N.D.; writing—review and editing, M.A., E.S., F.R.-M., C.A., C.S., B.B. and N.D.; visualization, M.A., E.S., F.R.-M. and N.D.; supervision, C.S., B.B. and N.D.; project administration, C.S., B.B. and N.D.; funding acquisition, C.S., B.B. and N.D. All authors have read and agreed to the published version of the manuscript.

Funding

This work and the PhD fellowship of Morgane Adélaı̈de were funded by the ANR-BMBF NATURAL-ARSENAL project. Funding number: ANR-19-AMRB-0004.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

We would like to thank the Matrics platform at the University Picardie Jules Verne for providing computing resources. We would also like to thank Dominique Cailleu for his help with setting up NMR experiments.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Tacconelli, E.; Carrara, E.; Savoldi, A.; Harbarth, S.; Mendelson, M.; Monnet, D.L.; Pulcini, C.; Kahlmeter, G.; Kluytmans, J.; Carmeli, Y.; et al. Discovery, Research, and Development of New Antibiotics: The WHO Priority List of Antibiotic-Resistant Bacteria and Tuberculosis. Lancet Infect. Dis. 2018, 18, 318–327. [Google Scholar] [CrossRef] [PubMed]
  2. Pendleton, J.N.; Gorman, S.P.; Gilmore, B.F. Clinical Relevance of the ESKAPE Pathogens. Expert Rev. Anti. Infect. Ther. 2013, 11, 297–308. [Google Scholar] [CrossRef] [PubMed]
  3. Mulani, M.S.; Kamble, E.E.; Kumkar, S.N.; Tawre, M.S.; Pardesi, K.R. Emerging Strategies to Combat ESKAPE Pathogens in the Era of Antimicrobial Resistance: A Review. Front. Microbiol. 2019, 10, 539. [Google Scholar] [CrossRef] [PubMed]
  4. Santajit, S.; Indrawattana, N. Mechanisms of Antimicrobial Resistance in ESKAPE Pathogens. BioMed Res. Int. 2016, 2016, 1–8. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Lewis, K. Persister Cells and Infectious Disease; Springer: Berlin/Heidelberg, Germany, 2019; ISBN 9783030252410. [Google Scholar]
  6. Ramos-Vivas, J.; Chapartegui-González, I.; Fernández-Martínez, M.; González-Rico, C.; Fortún, J.; Escudero, R.; Marco, F.; Linares, L.; Montejo, M.; Aranzamendi, M.; et al. Biofilm Formation by Multidrug Resistant Enterobacteriaceae Strains Isolated from Solid Organ Transplant Recipients. Sci. Rep. 2019, 9, 8928. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Lei, J.; Sun, L.; Huang, S.; Zhu, C.; Li, P.; He, J.; Mackey, V.; Coy, D.H.; He, Q. The Antimicrobial Peptides and Their Potential Clinical Applications. Am. J. Transl. Res. 2019, 11, 3919–3931. [Google Scholar] [PubMed]
  8. Ramos-Martín, F.; Herrera-León, C.; Antonietti, V.; Sonnet, P.; Sarazin, C.; D’Amelio, N. The Potential of Antifungal Peptide Sesquin as Natural Food Preservative. Biochimie 2022, 203, 51–64. [Google Scholar] [CrossRef]
  9. Rizzetto, G.; Gambini, D.; Maurizi, A.; Candelora, M.; Molinelli, E.; Cirioni, O.; Brescini, L.; Giacometti, A.; Offidani, A.; Simonetti, O. Our Experience over 20 Years: Antimicrobial Peptides against Gram Positives, Gram Negatives, and Fungi. Pharmaceutics 2022, 15, 40. [Google Scholar] [CrossRef]
  10. Zhu, S.; Sani, M.-A.; Separovic, F. Interaction of Cationic Antimicrobial Peptides from Australian Frogs with Lipid Membranes. Pept. Sci. 2018, 110, e24061. [Google Scholar] [CrossRef]
  11. Nguyen, L.T.; Haney, E.F.; Vogel, H.J. The Expanding Scope of Antimicrobial Peptide Structures and Their Modes of Action. Trends Biotechnol. 2011, 29, 464–472. [Google Scholar] [CrossRef]
  12. van Gent, M.E.; van der Reijden, T.J.K.; Lennard, P.R.; de Visser, A.W.; Schonkeren-Ravensbergen, B.; Dolezal, N.; Cordfunke, R.A.; Drijfhout, J.W.; Nibbering, P.H. Synergism between the Synthetic Antibacterial and Antibiofilm Peptide (SAAP)-148 and Halicin. Antibiotics 2022, 11, 673. [Google Scholar] [CrossRef]
  13. Corbett, D.; Wise, A.; Langley, T.; Skinner, K.; Trimby, E.; Birchall, S.; Dorali, A.; Sandiford, S.; Williams, J.; Warn, P.; et al. Potentiation of Antibiotic Activity by a Novel Cationic Peptide: Potency and Spectrum of Activity of SPR741. Antimicrob. Agents Chemother. 2017, 61, e00200-17. [Google Scholar] [CrossRef] [Green Version]
  14. Ramos-Martín, F.; D’Amelio, N. Molecular Basis of the Anticancer and Antibacterial Properties of CecropinXJ Peptide: An In Silico Study. Int. J. Mol. Sci. 2021, 22, 691. [Google Scholar] [CrossRef]
  15. Ramos-Martín, F.; Herrera-León, C.; D’Amelio, N. Molecular Basis of the Anticancer, Apoptotic and Antibacterial Activities of Bombyx Mori Cecropin A. Arch. Biochem. Biophys. 2022, 715, 109095. [Google Scholar] [CrossRef] [PubMed]
  16. Ramos-Martín, F.; Herrera-León, C.; D’Amelio, N. Bombyx Mori Cecropin D Could Trigger Cancer Cell Apoptosis by Interacting with Mitochondrial Cardiolipin. Biochim. Biophys. Acta Biomembr. 2022, 1864, 184003. [Google Scholar] [CrossRef] [PubMed]
  17. Langel, Ü. CPP, Cell-Penetrating Peptides; Springer: Berlin/Heidelberg, Germany, 2019; ISBN 9789811387470. [Google Scholar]
  18. Lee, H.; Lim, S.I.; Shin, S.-H.; Lim, Y.; Koh, J.W.; Yang, S. Conjugation of Cell-Penetrating Peptides to Antimicrobial Peptides Enhances Antibacterial Activity. ACS Omega 2019, 4, 15694–15701. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Joo, H.-S.; Fu, C.-I.; Otto, M. Bacterial Strategies of Resistance to Antimicrobial Peptides. Philos. Trans. R. Soc. B Biol. Sci. 2016, 371, 20150292. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. Yu, G.; Baeder, D.Y.; Regoes, R.R.; Rolff, J. Predicting Drug Resistance Evolution: Insights from Antimicrobial Peptides and Antibiotics. Proc. Biol. Sci. 2018, 285, 20172687. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  21. Moretta, A.; Scieuzo, C.; Petrone, A.M.; Salvia, R.; Manniello, M.D.; Franco, A.; Lucchetti, D.; Vassallo, A.; Vogel, H.; Sgambato, A.; et al. Antimicrobial Peptides: A New Hope in Biomedical and Pharmaceutical Fields. Front. Cell. Infect. Microbiol. 2021, 11, 668632. [Google Scholar] [CrossRef]
  22. Rodríguez-Rojas, A.; Moreno-Morales, J.; James Mason, A.; Rolff, J. Cationic Antimicrobial Peptides Do Not Change Recombination Frequency in Escherichia Coli. Biol. Lett. 2018, 14, 20180006. [Google Scholar] [CrossRef] [Green Version]
  23. Huang, H.W.; Charron, N.E. Understanding Membrane-Active Antimicrobial Peptides. Q. Rev. Biophys. 2017, 50, e10. [Google Scholar] [CrossRef] [PubMed]
  24. Lee, T.-H.; Hofferek, V.; Separovic, F.; Reid, G.E.; Aguilar, M.-I. The Role of Bacterial Lipid Diversity and Membrane Properties in Modulating Antimicrobial Peptide Activity and Drug Resistance. Curr. Opin. Chem. Biol. 2019, 52, 85–92. [Google Scholar] [CrossRef] [PubMed]
  25. Zasloff, M. Antimicrobial Peptides of Multicellular Organisms: My Perspective. Adv. Exp. Med. Biol. 2019, 1117, 3–6. [Google Scholar] [PubMed]
  26. Vrancianu, C.O.; Gheorghe, I.; Czobor, I.B.; Chifiriuc, M.C. Antibiotic Resistance Profiles, Molecular Mechanisms and Innovative Treatment Strategies of. Microorganisms 2020, 8, 935. [Google Scholar] [CrossRef] [PubMed]
  27. Matthyssen, T.; Li, W.; Holden, J.A.; Lenzo, J.C.; Hadjigol, S.; O’Brien-Simpson, N.M. The Potential of Modified and Multimeric Antimicrobial Peptide Materials as Superbug Killers. Front. Chem. 2021, 9, 795433. [Google Scholar] [CrossRef]
  28. Gao, Y.; Fang, H.; Fang, L.; Liu, D.; Liu, J.; Su, M.; Fang, Z.; Ren, W.; Jiao, H. The Modification and Design of Antimicrobial Peptide. Curr. Pharm. Des. 2018, 24, 904–910. [Google Scholar] [CrossRef]
  29. Sandín, D.; Valle, J.; Chaves-Arquero, B.; Prats-Ejarque, G.; Larrosa, M.N.; González-López, J.J.; Jiménez, M.Á.; Boix, E.; Andreu, D.; Torrent, M. Rationally Modified Antimicrobial Peptides from the N-Terminal Domain of Human RNase 3 Show Exceptional Serum Stability. J. Med. Chem. 2021, 64, 11472–11482. [Google Scholar] [CrossRef]
  30. Han, Y.; Zhang, M.; Lai, R.; Zhang, Z. Chemical Modifications to Increase the Therapeutic Potential of Antimicrobial Peptides. Peptides 2021, 146, 170666. [Google Scholar] [CrossRef]
  31. de Breij, A.; Riool, M.; Cordfunke, R.A.; Malanovic, N.; de Boer, L.; Koning, R.I.; Ravensbergen, E.; Franken, M.; van der Heijde, T.; Boekema, B.K.; et al. The Antimicrobial Peptide SAAP-148 Combats Drug-Resistant Bacteria and Biofilms. Sci. Transl. Med. 2018, 10, eaan4044. [Google Scholar] [CrossRef] [Green Version]
  32. Scheper, H.; Wubbolts, J.M.; Verhagen, J.A.M.; de Visser, A.W.; van der Wal, R.J.P.; Visser, L.G.; de Boer, M.G.J.; Nibbering, P.H. SAAP-148 Eradicates MRSA Persisters Within Mature Biofilm Models Simulating Prosthetic Joint Infection. Front. Microbiol. 2021, 12, 625952. [Google Scholar] [CrossRef]
  33. Piller, P.; Wolinski, H.; Cordfunke, R.A.; Drijfhout, J.W.; Keller, S.; Lohner, K.; Malanovic, N. Membrane Activity of LL-37 Derived Antimicrobial Peptides against: Superiority of SAAP-148 over OP-145. Biomolecules 2022, 12, 523. [Google Scholar] [CrossRef] [PubMed]
  34. Bechinger, B. The SMART Model: Soft Membranes Adapt and Respond, Also Transiently, in the Presence of Antimicrobial Peptides. J. Pept. Sci. 2015, 21, 346–355. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Hammond, K.; Ryadnov, M.G.; Hoogenboom, B.W. Atomic Force Microscopy to Elucidate How Peptides Disrupt Membranes. Biochim. Biophys. Acta Biomembr. 2021, 1863, 183447. [Google Scholar] [CrossRef] [PubMed]
  36. Henderson, J.M.; Henderson, J.M.; Waring, A.J.; Separovic, F.; Lee, K.Y.C. Antimicrobial Peptides Share a Common Interaction Driven by Membrane Line Tension Reduction. Biophys. J. 2016, 111, 2176–2189. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Kim, C.; Spano, J.; Park, E.-K.; Wi, S. Evidence of Pores and Thinned Lipid Bilayers Induced in Oriented Lipid Membranes Interacting with the Antimicrobial Peptides, Magainin-2 and Aurein-3.3. Biochim. Biophys. Acta 2009, 1788, 1482–1496. [Google Scholar] [CrossRef] [Green Version]
  38. Zeth, K.; Sancho-Vaello, E. The Human Antimicrobial Peptides Dermcidin and LL-37 Show Novel Distinct Pathways in Membrane Interactions. Front. Chem. 2017, 5, 86. [Google Scholar] [CrossRef] [Green Version]
  39. Rakowska, P.D.; Jiang, H.; Ray, S.; Pyne, A.; Lamarre, B.; Carr, M.; Judge, P.J.; Ravi, J.; Gerling, U.I.M.; Koksch, B.; et al. Nanoscale Imaging Reveals Laterally Expanding Antimicrobial Pores in Lipid Bilayers. Proc. Natl. Acad. Sci. USA 2013, 110, 8918–8923. [Google Scholar] [CrossRef] [Green Version]
  40. Priyadarshini, D.; Ivica, J.; Separovic, F.; de Planque, M.R.R. Characterisation of Cell Membrane Interaction Mechanisms of Antimicrobial Peptides by Electrical Bilayer Recording. Biophys. Chem. 2022, 281, 106721. [Google Scholar] [CrossRef]
  41. Huang, Y.-T.; Kumar, S.R.; Chan, H.-C.; Jhan, Z.-H.; Chen, D.W.; Lue, S.J. Efficacy of Antimicrobial Peptides (AMPs) against Escherichia Coli and Bacteria Morphology Change after AMP Exposure. J. Taiwan Inst. Chem. Eng. 2021, 126, 307–312. [Google Scholar] [CrossRef]
  42. Lee, C.-C.; Sun, Y.; Qian, S.; Huang, H.W. Transmembrane Pores Formed by Human Antimicrobial Peptide LL-37. Biophys. J. 2011, 100, 1688–1696. [Google Scholar] [CrossRef] [Green Version]
  43. Savini, F.; Loffredo, M.R.; Troiano, C.; Bobone, S.; Malanovic, N.; Eichmann, T.O.; Caprio, L.; Canale, V.C.; Park, Y.; Mangoni, M.L.; et al. Binding of an Antimicrobial Peptide to Bacterial Cells: Interaction with Different Species, Strains and Cellular Components. Biochim. Biophys. Acta Biomembr. 2020, 1862, 183291. [Google Scholar] [CrossRef] [PubMed]
  44. Saigo, N.; Izumi, K.; Kawano, R. Electrophysiological Analysis of Antimicrobial Peptides in Diverse Species. ACS Omega 2019, 4, 13124–13130. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Kordi, M.; Borzouyi, Z.; Chitsaz, S.; Asmaei, M.H.; Salami, R.; Tabarzad, M. Antimicrobial Peptides with Anticancer Activity: Today Status, Trends and Their Computational Design. Arch. Biochem. Biophys. 2023, 733, 109484. [Google Scholar] [CrossRef] [PubMed]
  46. Malanovic, N.; Marx, L.; Blondelle, S.E.; Pabst, G.; Semeraro, E.F. Experimental Concepts for Linking the Biological Activities of Antimicrobial Peptides to Their Molecular Modes of Action. Biochim. Biophys. Acta Biomembr. 2020, 1862, 183275. [Google Scholar] [CrossRef]
  47. Huang, H.W. DAPTOMYCIN, Its Membrane-Active Mechanism vs. that of Other Antimicrobial Peptides. Biochim. Biophys. Acta Biomembr. 2020, 1862, 183395. [Google Scholar] [CrossRef]
  48. Marquette, A.; Bechinger, B. Biophysical Investigations Elucidating the Mechanisms of Action of Antimicrobial Peptides and Their Synergism. Biomolecules 2018, 8, 18. [Google Scholar] [CrossRef] [Green Version]
  49. Aisenbrey, C.; Marquette, A.; Bechinger, B. The Mechanisms of Action of Cationic Antimicrobial Peptides Refined by Novel Concepts from Biophysical Investigations. Adv. Exp. Med. Biol. 2019, 1117, 33–64. [Google Scholar]
  50. Nielsen, J.T.; Mulder, F.A.A. POTENCI: Prediction of Temperature, Neighbor and pH-Corrected Chemical Shifts for Intrinsically Disordered Proteins. J. Biomol. NMR 2018, 70, 141–165. [Google Scholar] [CrossRef]
  51. Aisenbrey, C.; Bertani, P.; Bechinger, B. Solid-State NMR Investigations of Membrane-Associated Antimicrobial Peptides. Methods Mol. Biol. 2010, 618, 209–233. [Google Scholar]
  52. Bechinger, B.; Opella, S.J. Flat-Coil Probe for NMR Spectroscopy of Oriented Membrane Samples. J. Magn. Reson. 1991, 95, 585–588. [Google Scholar] [CrossRef]
  53. Rance, M.; Byrd, R.A. Obtaining High-Fidelity Powder Spectra in Anisotropic Media: Phase-Cycled Hahn Echo Spectroscopy. J. Magn. Reson. 1983, 52, 221–240. [Google Scholar] [CrossRef]
  54. Pines, A.; Gibby, M.G.; Waugh, J.S. Proton-enhanced NMR of Dilute Spins in Solids. J. Chem. Phys. 1973, 59, 569–590. [Google Scholar] [CrossRef]
  55. Bertani, P.; Raya, J.; Bechinger, B. 15N Chemical Shift Referencing in Solid State NMR. Solid State Nucl. Magn. Reson. 2014, 61-62, 15–18. [Google Scholar] [CrossRef] [PubMed]
  56. Davis, J.H.; Jeffrey, K.R.; Bloom, M.; Valic, M.I.; Higgs, T.P. Quadrupolar Echo Deuteron Magnetic Resonance Spectroscopy in Ordered Hydrocarbon Chains. Chem. Phys. Lett. 1976, 42, 390–394. [Google Scholar] [CrossRef]
  57. Batchelder, L.S.; Niu, C.H.; Torchia, D.A. Methyl Reorientation in Polycrystalline Amino Acids and Peptides: A Deuteron NMR Spin-Lattice Relaxation Study. J. Am. Chem. Soc. 1983, 105, 2228–2231. [Google Scholar] [CrossRef]
  58. Michalek, M.; Salnikov, E.S.; Bechinger, B. Structure and Topology of the Huntingtin 1–17 Membrane Anchor by a Combined Solution and Solid-State NMR Approach. Biophys. J. 2013, 105, 699–710. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  59. Jo, S.; Lim, J.B.; Klauda, J.B.; Im, W. CHARMM-GUI Membrane Builder for Mixed Bilayers and Its Application to Yeast Membranes. Biophys. J. 2009, 97, 50–58. [Google Scholar] [CrossRef] [Green Version]
  60. Lee, J.; Cheng, X.; Swails, J.M.; Yeom, M.S.; Eastman, P.K.; Lemkul, J.A.; Wei, S.; Buckner, J.; Jeong, J.C.; Qi, Y.; et al. CHARMM-GUI Input Generator for NAMD, GROMACS, AMBER, OpenMM, and CHARMM/OpenMM Simulations Using the CHARMM36 Additive Force Field. J. Chem. Theory Comput. 2016, 12, 405–413. [Google Scholar] [CrossRef]
  61. Wu, E.L.; Cheng, X.; Jo, S.; Rui, H.; Song, K.C.; Dávila-Contreras, E.M.; Qi, Y.; Lee, J.; Monje-Galvan, V.; Venable, R.M.; et al. CHARMM-GUI Membrane Builder toward Realistic Biological Membrane Simulations. J. Comput. Chem. 2014, 35, 1997–2004. [Google Scholar] [CrossRef] [Green Version]
  62. Abraham, M.J.; Murtola, T.; Schulz, R.; Páll, S.; Smith, J.C.; Hess, B.; Lindahl, E. GROMACS: High Performance Molecular Simulations through Multi-Level Parallelism from Laptops to Supercomputers. SoftwareX 2015, 1-2, 19–25. [Google Scholar] [CrossRef] [Green Version]
  63. Huang, J.; Rauscher, S.; Nawrocki, G.; Ran, T.; Feig, M.; de Groot, B.L.; Grubmüller, H.; MacKerell, A.D. CHARMM36m: An Improved Force Field for Folded and Intrinsically Disordered Proteins. Nat. Methods 2016, 14, 71–73. [Google Scholar] [CrossRef] [Green Version]
  64. Berendsen, H.J.C.; Postma, J.P.M.; van Gunsteren, W.F.; Hermans, J. Interaction models for water in relation to protein hydration. In Intermolecular Forces, Proceedings of the Fourteenth Jerusalem Symposium on Quantum Chemistry and Biochemistry, Jerusalem, Israel, 13–16 April 1981; Springer: Dordrecht, The Netherlands; pp. 331–342. [CrossRef]
  65. Yang, J.; Yan, R.; Roy, A.; Xu, D.; Poisson, J.; Zhang, Y. The I-TASSER Suite: Protein Structure and Function Prediction. Nat. Methods 2015, 12, 7–8. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Berendsen, H.J.C.; Postma, J.P.M.; van Gunsteren, W.F.; DiNola, A.; Haak, J.R. Molecular Dynamics with Coupling to an External Bath. J. Chem. Phys. 1984, 81, 3684–3690. [Google Scholar] [CrossRef] [Green Version]
  67. Parrinello, M.; Rahman, A. Polymorphic Transitions in Single Crystals: A New Molecular Dynamics Method. J. Appl. Phys. 1981, 52, 7182–7190. [Google Scholar] [CrossRef]
  68. Nosé, S. A Unified Formulation of the Constant Temperature Molecular Dynamics Methods. J. Chem. Phys. 1984, 81, 511–519. [Google Scholar] [CrossRef] [Green Version]
  69. Hoover, W.G. Canonical Dynamics: Equilibrium Phase-Space Distributions. Phys. Rev. A Gen. Phys. 1985, 31, 1695–1697. [Google Scholar] [CrossRef] [Green Version]
  70. Jo, S.; Kim, T.; Im, W. Automated Builder and Database of Protein/membrane Complexes for Molecular Dynamics Simulations. PLoS ONE 2007, 2, e880. [Google Scholar] [CrossRef] [Green Version]
  71. Williams, T.; Kelley, C.; Merritt, E.A.; Bersch, C.; Bröker, H.B.; Campbell, J.; Cunningham, R.; Denholm, D.; Elber, E.; Fearick, R. Gnuplot 5.4. 4: An Interactive Plotting Program; 2022. Available online: http://www.gnuplot.info/documentation.html (accessed on 1 May 2022).
  72. DeLano, W.L. Pymol: An Open-Source Molecular Graphics Tool. CCP4 Newsl. Protein Crystallogr. 2002, 40, 82–92. [Google Scholar]
  73. Koradi, R.; Billeter, M.; Wüthrich, K. MOLMOL: A Program for Display and Analysis of Macromolecular Structures. J. Mol. Graph. 1996, 14, 51–55. [Google Scholar] [CrossRef]
  74. Humphrey, W.; Dalke, A.; Schulten, K. VMD: Visual Molecular Dynamics. J. Mol. Graph. 1996, 14, 33–38. [Google Scholar] [CrossRef]
  75. Han, B.; Liu, Y.; Ginzinger, S.W.; Wishart, D.S. SHIFTX2: Significantly Improved Protein Chemical Shift Prediction. J. Biomol. NMR 2011, 50, 43–57. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  76. Porcelli, F.; Ramamoorthy, A.; Barany, G.; Veglia, G. On the Role of NMR Spectroscopy for Characterization of Antimicrobial Peptides. Membr. Proteins 2013, 159–180. [Google Scholar] [CrossRef] [Green Version]
  77. Marcotte, I.; Auger, M. Bicelles as Model Membranes for Solid-and Solution-State NMR Studies of Membrane Peptides and Proteins. Concepts Magn. Reson. Part A Educ. J. 2005, 24, 17–37. [Google Scholar] [CrossRef]
  78. Randle, C.L.; Albro, P.W.; Dittmer, J.C. The Phosphoglyceride Composition of Gram-Negative Bacteria and the Changes in Composition during Growth. Biochim. Biophys. Acta 1969, 187, 214–220. [Google Scholar] [CrossRef] [PubMed]
  79. Ramos-Martín, F.; D’Amelio, N. Biomembrane Lipids: When Physics and Chemistry Join to Shape Biological Activity. Biochimie 2022, 203, 118–138. [Google Scholar] [CrossRef] [PubMed]
  80. Wang, Y.; Jardetzky, O. Probability-Based Protein Secondary Structure Identification Using Combined NMR Chemical-Shift Data. Protein Sci. 2002, 11, 852–861. [Google Scholar] [CrossRef]
  81. Wishart, D.S.; Sykes, B.D.; Richards, F.M. The Chemical Shift Index: A Fast and Simple Method for the Assignment of Protein Secondary Structure through NMR Spectroscopy. Biochemistry 1992, 31, 1647–1651. [Google Scholar] [CrossRef]
  82. Wishart, D.S.; Sykes, B.D. The 13 C Chemical-Shift Index: A Simple Method for the Identification of Protein Secondary Structure Using 13 C Chemical-Shift Data. J. Biomol. NMR 1994, 4, 171–180. [Google Scholar] [CrossRef]
  83. Wishart, D.S. Interpreting Protein Chemical Shift Data. Prog. Nucl. Magn. Reson. Spectrosc. 2011, 58, 62–87. [Google Scholar] [CrossRef]
  84. Hilty, C.; Wider, G.; Fernández, C.; Wüthrich, K. Membrane Protein-Lipid Interactions in Mixed Micelles Studied by NMR Spectroscopy with the Use of Paramagnetic Reagents. Chembiochem 2004, 5, 467–473. [Google Scholar] [CrossRef]
  85. Grisham, C.M. Paramagnetic Probes in NMR and EPR Studies of Membrane Enzymes. J. Biochem. Biophys. Methods 1980, 3, 39–59. [Google Scholar] [CrossRef] [PubMed]
  86. Sommer, L.A.M.; Joel Janke, J.; Drew Bennett, W.F.; Bürck, J.; Ulrich, A.S.; Peter Tieleman, D.; Dames, S.A. Characterization of the Immersion Properties of the Peripheral Membrane Anchor of the FATC Domain of the Kinase “Target of Rapamycin” by NMR, Oriented CD Spectroscopy, and MD Simulations. J. Phys. Chem. B 2014, 118, 4817–4831. [Google Scholar] [CrossRef] [PubMed]
  87. Ramos-Martín, F.; Herrera-León, C.; Antonietti, V.; Sonnet, P.; Sarazin, C.; D’Amelio, N. Antimicrobial Peptide K11 Selectively Recognizes Bacterial Biomimetic Membranes and Acts by Twisting Their Bilayers. Pharmaceuticals 2020, 14, 1. [Google Scholar] [CrossRef] [PubMed]
  88. Mól, A.R.; Castro, M.S.; Fontes, W. NetWheels: A Web Application to Create High Quality Peptide Helical Wheel and Net Projections. BioRxiv 2018. [Google Scholar] [CrossRef] [Green Version]
  89. Bechinger, B.; Salnikov, E.S. The Membrane Interactions of Antimicrobial Peptides Revealed by Solid-State NMR Spectroscopy. Chem. Phys. Lipids 2012, 165, 282–301. [Google Scholar] [CrossRef]
  90. Bechinger, B.; Sizun, C. Alignment and Structural Analysis of Membrane Polypeptides by15N and31P Solid-State NMR Spectroscopy. Concepts Magn. Reson. 2003, 18A, 130–145. [Google Scholar] [CrossRef]
  91. Salnikov, E.S.; James Mason, A.; Bechinger, B. Membrane Order Perturbation in the Presence of Antimicrobial Peptides by 2H Solid-State NMR Spectroscopy. Biochimie 2009, 91, 734–743. [Google Scholar] [CrossRef]
  92. Davis, J.H. The Description of Membrane Lipid Conformation, Order and Dynamics by 2H-NMR. Biochim. Biophys. Acta 1983, 737, 117–171. [Google Scholar] [CrossRef]
  93. Harmouche, N.; Bechinger, B. Lipid-Mediated Interactions between the Antimicrobial Peptides Magainin 2 and PGLa in Bilayers. Biophys. J. 2018, 115, 1033–1044. [Google Scholar] [CrossRef] [Green Version]
  94. Herrera-León, C.; Ramos-Martín, F.; Antonietti, V.; Sonnet, P.; D’Amelio, N. The Impact of Phosphatidylserine Exposure on Cancer Cell Membranes on the Activity of the Anticancer Peptide HB43. FEBS J. 2022, 289, 1984–2003. [Google Scholar] [CrossRef]
  95. Herrera-León, C.; Ramos-Martín, F.; El Btaouri, H.; Antonietti, V.; Sonnet, P.; Martiny, L.; Zevolini, F.; Falciani, C.; Sarazin, C.; D’Amelio, N. The Influence of Short Motifs on the Anticancer Activity of HB43 Peptide. Pharmaceutics 2022, 14, 1089. [Google Scholar] [CrossRef] [PubMed]
  96. Vahedi, A.; Bigdelou, P.; Farnoud, A.M. Quantitative Analysis of Red Blood Cell Membrane Phospholipids and Modulation of Cell-Macrophage Interactions Using Cyclodextrins. Sci. Rep. 2020, 10, 15111. [Google Scholar] [CrossRef]
  97. Li, G.; Kim, J.; Huang, Z.; St Clair, J.R.; Brown, D.A.; London, E. Efficient Replacement of Plasma Membrane Outer Leaflet Phospholipids and Sphingolipids in Cells with Exogenous Lipids. Proc. Natl. Acad. Sci. USA 2016, 113, 14025–14030. [Google Scholar] [CrossRef] [Green Version]
  98. Luchini, A.; Vitiello, G. Mimicking the Mammalian Plasma Membrane: An Overview of Lipid Membrane Models for Biophysical Studies. Biomimetics 2020, 6, 3. [Google Scholar] [CrossRef]
  99. Zachowski, A. Phospholipids in Animal Eukaryotic Membranes: Transverse Asymmetry and Movement. Biochem. J. 1993, 294, 1–14. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  100. Van Meer, G.; Voelker, D.R.; Feigenson, G.W. Membrane Lipids: Where They Are and How They Behave. Nat. Rev. Mol. Cell Biol. 2008, 9, 112–124. [Google Scholar] [CrossRef] [PubMed]
  101. Vance, J.E. Phospholipid Synthesis and Transport in Mammalian Cells. Traffic 2015, 16, 1–18. [Google Scholar] [CrossRef] [PubMed]
  102. Yang, Y.; Jalali, S.; Nilsson, B.L.; Dias, C.L. Binding Mechanisms of Amyloid-like Peptides to Lipid Bilayers and Effects of Divalent Cations. ACS Chem. Neurosci. 2021, 12, 2027–2035. [Google Scholar] [CrossRef]
  103. Sciacca, M.F.M.; Monaco, I.; La Rosa, C.; Milardi, D. The Active Role of Ca Ions in Aβ-Mediated Membrane Damage. Chem. Commun. 2018, 54, 3629–3631. [Google Scholar] [CrossRef]
  104. Liu, Y.; Ren, B.; Zhang, Y.; Sun, Y.; Chang, Y.; Liang, G.; Xu, L.; Zheng, J. Molecular Simulation Aspects of Amyloid Peptides at Membrane Interface. Biochim. Biophys. Acta Biomembr. 2018, 1860, 1906–1916. [Google Scholar] [CrossRef]
  105. Lockhart, C.; Klimov, D.K. Calcium Enhances Binding of Aβ Monomer to DMPC Lipid Bilayer. Biophys. J. 2015, 108, 1807–1818. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  106. Lemkul, J.A.; Bevan, D.R. Perturbation of Membranes by the Amyloid β-Peptide—A Molecular Dynamics Study. FEBS J. 2009, 276, 3060–3075. [Google Scholar] [CrossRef] [PubMed]
  107. Khelashvili, G.; Plante, A.; Doktorova, M.; Weinstein, H. Ca2-Dependent Mechanism of Membrane Insertion and Destabilization by the SARS-CoV-2 Fusion Peptide. Biophys. J. 2021, 120, 1105–1119. [Google Scholar] [CrossRef] [PubMed]
  108. Bowdish, D.M.E.; Davidson, D.J.; Elaine Lau, Y.; Lee, K.; Scott, M.G.; Hancock, R.E.W. Impact of LL-37 on Anti-Infective Immunity. J. Leukoc. Biol. 2005, 77, 451–459. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  109. Dufourc, E.J.; Smith, I.C.; Dufourcq, J. Molecular Details of Melittin-Induced Lysis of Phospholipid Membranes as Revealed by Deuterium and Phosphorus NMR. Biochemistry 1986, 25, 6448–6455. [Google Scholar] [CrossRef]
  110. Henzler-Wildman, K.A.; Martinez, G.V.; Brown, M.F.; Ramamoorthy, A. Perturbation of the Hydrophobic Core of Lipid Bilayers by the Human Antimicrobial Peptide LL-37. Biochemistry 2004, 43, 8459–8469. [Google Scholar] [CrossRef]
  111. Hénin, J.; Lelièvre, T.; Shirts, M.R.; Valsson, O.; Delemotte, L. Enhanced Sampling Methods for Molecular Dynamics Simulations [Article v1.0]. Living J. Comp. Mol. Sci. 2022, 4, 1583. [Google Scholar] [CrossRef]
  112. Pearlstein, R.A.; Dickson, C.J.; Hornak, V. Contributions of the Membrane Dipole Potential to the Function of Voltage-Gated Cation Channels and Modulation by Small Molecule Potentiators. Biochim. Et Biophys. Acta (BBA) Biomembr. 2017, 1859, 177–194. [Google Scholar] [CrossRef]
  113. Dreyer, J.; Zhang, C.; Ippoliti, E.; Carloni, P. Role of the Membrane Dipole Potential for Proton Transport in Gramicidin A Embedded in a DMPC Bilayer. J. Chem. Theory Comput. 2013, 9, 3826–3831. [Google Scholar] [CrossRef]
  114. Bechinger, B.; Juhl, D.W.; Glattard, E.; Aisenbrey, C. Revealing the Mechanisms of Synergistic Action of Two Magainin Antimicrobial Peptides. Front. Med. Technol. 2020, 2, 615494. [Google Scholar] [CrossRef]
Pharmaceutics 15 00761 g001 550
Figure 1. (A) The 500 MHz 1H NMR spectra (amide and aromatic region) recording during the titration of SAAP-148 with DPC:d38 micelles up to a concentration of 60 mM. (B) NMR assignment of the amide region of SAAP-148 in DPC micelles. (C) The 500 MHz 1H NMR spectra of the titration of SAAP-148 with DHPC/DMPC/DMPG isotropic bicelles up to a concentration of 50 mM. The temperature was 310 K.
Figure 1. (A) The 500 MHz 1H NMR spectra (amide and aromatic region) recording during the titration of SAAP-148 with DPC:d38 micelles up to a concentration of 60 mM. (B) NMR assignment of the amide region of SAAP-148 in DPC micelles. (C) The 500 MHz 1H NMR spectra of the titration of SAAP-148 with DHPC/DMPC/DMPG isotropic bicelles up to a concentration of 50 mM. The temperature was 310 K.
Pharmaceutics 15 00761 g001
Pharmaceutics 15 00761 g002 550
Figure 2. (A,B) Hα (A) and Cα (B) chemical shift deviations from random coil values (see Materials and Methods). (C) Percental intensity loss of Hα/Cα 1H,13C-HMQC peaks in the presence of paramagnetic 16-DSA. Insert shows a helical wheel for SAAP-148 peptide (created using Netwheels web server [88]).
Figure 2. (A,B) Hα (A) and Cα (B) chemical shift deviations from random coil values (see Materials and Methods). (C) Percental intensity loss of Hα/Cα 1H,13C-HMQC peaks in the presence of paramagnetic 16-DSA. Insert shows a helical wheel for SAAP-148 peptide (created using Netwheels web server [88]).
Pharmaceutics 15 00761 g002
Pharmaceutics 15 00761 g003 550
Figure 3. Proton-decoupled 31P (A,B) and 15N (C,D) spectra of 2 mol% SAAP-148 carrying 15N label in position Leu11 (A,C) and Leu12 (B,D) in POPE:POPG 3:1 at 310 K. The 15N signal from the PE headgroup at natural isotope abundance is indicated by asterisks in panels (C) and (D). (E) Restriction analysis and corresponding peptide alignments relative to the membrane normal (F) for 2 mol% SAAP-148 in POPE:POPG 3:1 bilayer. The resulting restrictions for SAAP-148 topology are shown in panel E: 15N-Leu11 is 80.0 ± 2.5 ppm in blue, and 15N-Leu12 is 70.0 ± 4.0 ppm in green. The restraints were obtained assuming rocking (18° Gaussian distribution) and wobbling (10°) motions of an ideal helix (φ = −65°, ψ = −45°). Leucines, phenylalanines, valines, tryptophan, and proline are shown in yellow, lysines in blue, glutamines in magenta, and arginines in cyan. Two possible intersections are obtained in (E) (orange ellipses), both are represented in panel (F) (note that the perfect helix shown in panel F does not represent an experimental structure but a way to show the axis orientation in the assumption of alpha helix conformation).
Figure 3. Proton-decoupled 31P (A,B) and 15N (C,D) spectra of 2 mol% SAAP-148 carrying 15N label in position Leu11 (A,C) and Leu12 (B,D) in POPE:POPG 3:1 at 310 K. The 15N signal from the PE headgroup at natural isotope abundance is indicated by asterisks in panels (C) and (D). (E) Restriction analysis and corresponding peptide alignments relative to the membrane normal (F) for 2 mol% SAAP-148 in POPE:POPG 3:1 bilayer. The resulting restrictions for SAAP-148 topology are shown in panel E: 15N-Leu11 is 80.0 ± 2.5 ppm in blue, and 15N-Leu12 is 70.0 ± 4.0 ppm in green. The restraints were obtained assuming rocking (18° Gaussian distribution) and wobbling (10°) motions of an ideal helix (φ = −65°, ψ = −45°). Leucines, phenylalanines, valines, tryptophan, and proline are shown in yellow, lysines in blue, glutamines in magenta, and arginines in cyan. Two possible intersections are obtained in (E) (orange ellipses), both are represented in panel (F) (note that the perfect helix shown in panel F does not represent an experimental structure but a way to show the axis orientation in the assumption of alpha helix conformation).
Pharmaceutics 15 00761 g003
Pharmaceutics 15 00761 g004 550
Figure 4. (A,B) Quadrupolar 2H solid-state NMR spectra of chain-deuterated POPG-d31 (A) or POPE-d31 (B) in the presence or absence of SAAP-148 peptide at 2 mol% concentration at 310 K. Samples were hydrated with Tris buffer pH 7.4 to reach a water content of h = 0.81. (C,D) Order parameters SCD of lipid palmitoyl chain C–D bonds in the absence and in the presence of SAAP-148 for POPG-d31 (C) and POPE-d31 (D). The effect of SAAP-148 is more clearly evidenced when plotting the ratio of SCD in the presence of SAAP-148 with that in its absence (E,F).
Figure 4. (A,B) Quadrupolar 2H solid-state NMR spectra of chain-deuterated POPG-d31 (A) or POPE-d31 (B) in the presence or absence of SAAP-148 peptide at 2 mol% concentration at 310 K. Samples were hydrated with Tris buffer pH 7.4 to reach a water content of h = 0.81. (C,D) Order parameters SCD of lipid palmitoyl chain C–D bonds in the absence and in the presence of SAAP-148 for POPG-d31 (C) and POPE-d31 (D). The effect of SAAP-148 is more clearly evidenced when plotting the ratio of SCD in the presence of SAAP-148 with that in its absence (E,F).
Pharmaceutics 15 00761 g004
Pharmaceutics 15 00761 g005 550
Figure 5. (A,B) Recurrence of polar atom contacts (H bonds and salt bridges) between SAAP-148 peptide and POPC (A) and POPE/POPG (B) bilayers calculated along MD simulation trajectories. (C) Recurrence of van der Waals contacts between SAAP-148 peptide and POPE/POPG bilayers calculated along MD simulation trajectories. (D) MD snapshot of SAAP-148 interacting with POPE/POPG bilayers in the presence of Ca2+ ion. Color code: phosphorus atom: yellow, POPE: dark green (body), turquoise (headgroup), light green (amine of the headgroup); POPG: dark violet (body), violet (headgroup), light violet (hydroxyls of the headgroup). For clarity, only functional groups of headgroups are shown (spheres) in the upper leaflet. SAAP-148 is shown as a ‘tube’ colored from blue (N-terminus) to red (C-terminus). Sidechains are shown as sticks with the following color code: positively charged (blue) and nonpolar (light gray). (E) SAAP-148 dimers found along the trajectory and interacting with POPE/POPG bilayers. In the dimer, the two units (colored in red and blue) can arrange in antiparallel (I) or parallel fashion (II).
Figure 5. (A,B) Recurrence of polar atom contacts (H bonds and salt bridges) between SAAP-148 peptide and POPC (A) and POPE/POPG (B) bilayers calculated along MD simulation trajectories. (C) Recurrence of van der Waals contacts between SAAP-148 peptide and POPE/POPG bilayers calculated along MD simulation trajectories. (D) MD snapshot of SAAP-148 interacting with POPE/POPG bilayers in the presence of Ca2+ ion. Color code: phosphorus atom: yellow, POPE: dark green (body), turquoise (headgroup), light green (amine of the headgroup); POPG: dark violet (body), violet (headgroup), light violet (hydroxyls of the headgroup). For clarity, only functional groups of headgroups are shown (spheres) in the upper leaflet. SAAP-148 is shown as a ‘tube’ colored from blue (N-terminus) to red (C-terminus). Sidechains are shown as sticks with the following color code: positively charged (blue) and nonpolar (light gray). (E) SAAP-148 dimers found along the trajectory and interacting with POPE/POPG bilayers. In the dimer, the two units (colored in red and blue) can arrange in antiparallel (I) or parallel fashion (II).
Pharmaceutics 15 00761 g005
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Adélaïde, M.; Salnikov, E.; Ramos-Martín, F.; Aisenbrey, C.; Sarazin, C.; Bechinger, B.; D’Amelio, N. The Mechanism of Action of SAAP-148 Antimicrobial Peptide as Studied with NMR and Molecular Dynamics Simulations. Pharmaceutics 2023, 15, 761. https://doi.org/10.3390/pharmaceutics15030761

AMA Style

Adélaïde M, Salnikov E, Ramos-Martín F, Aisenbrey C, Sarazin C, Bechinger B, D’Amelio N. The Mechanism of Action of SAAP-148 Antimicrobial Peptide as Studied with NMR and Molecular Dynamics Simulations. Pharmaceutics. 2023; 15(3):761. https://doi.org/10.3390/pharmaceutics15030761

Chicago/Turabian Style

Adélaïde, Morgane, Evgeniy Salnikov, Francisco Ramos-Martín, Christopher Aisenbrey, Catherine Sarazin, Burkhard Bechinger, and Nicola D’Amelio. 2023. "The Mechanism of Action of SAAP-148 Antimicrobial Peptide as Studied with NMR and Molecular Dynamics Simulations" Pharmaceutics 15, no. 3: 761. https://doi.org/10.3390/pharmaceutics15030761

APA Style

Adélaïde, M., Salnikov, E., Ramos-Martín, F., Aisenbrey, C., Sarazin, C., Bechinger, B., & D’Amelio, N. (2023). The Mechanism of Action of SAAP-148 Antimicrobial Peptide as Studied with NMR and Molecular Dynamics Simulations. Pharmaceutics, 15(3), 761. https://doi.org/10.3390/pharmaceutics15030761

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop