# Temporin L and aurein 2.5 have identical conformations but subtly distinct membrane and antibacterial activities

Giorgia Manzo, Philip M Ferguson, Charlotte K Hind, Melanie Clifford, V. Benjamin Benjamin Gustilo, Hind Ali, Sukhvinder S Bansal, Tam T Bui, Alex F Drake, Robert Andrew Atkinson, et al.

► To cite this version:

Giorgia Manzo, Philip M Ferguson, Charlotte K Hind, Melanie Clifford, V. Benjamin Benjamin Gustilo, et al. Temporin L and aurein 2.5 have identical conformations but subtly distinct membrane and antibacterial activities. Scientific Reports, 2019, 9 (1), 10.1038/s41598-019-47327-w. hal-03027493

HAL Id: hal-03027493

https://hal.science/hal-03027493

Submitted on 27 Nov 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
OPEN

# Temporin L and aurein 2.5 have identical conformations but subtly distinct membrane and antibacterial activities

Received: 13 August 2018

Accepted: 16 July 2019

Published online: 29 July 2019

Giorgia Manzo<fnref n="1" />, Philip M. Ferguson<fnref n="1" />, Charlotte K. Hind², Melanie Clifford², V. Benjamin Gustilo<fnref n="1" />, Hind Ali<fnref n="1" />, Sukhvinder S. Bansal<fnref n="1" />, Tam T. Bui³, Alex F. Drake³, R. Andrew Atkinson D<fnref n="1" />, J. Mark Sutton², Christian D. Lorenz⁴, David A. Phoenix⁵ & A. James Mason D<fnref n="1" />

Frogs such as *Rana temporaria* and *Litoria aurea* secrete numerous closely related antimicrobial peptides (AMPs) as an effective chemical dermal defence. Damage or penetration of the bacterial plasma membrane is considered essential for AMP activity and such properties are commonly ascribed to their ability to form secondary amphipathic, α-helix conformations in membrane mimicking milieu. Nevertheless, despite the high similarity in physical properties and preference for adopting such conformations, the spectrum of activity and potency of AMPs often varies considerably. Hence distinguishing apparently similar AMPs according to their behaviour in, and effects on, model membranes will inform understanding of primary-sequence-specific antimicrobial mechanisms. Here we use a combination of molecular dynamics simulations, circular dichroism and patch-clamp to investigate the basis for differing anti-bacterial activities in representative AMPs from each species; temporin L and aurein 2.5. Despite adopting near identical, α-helix conformations in the steady-state in a variety of membrane models, these two AMPs can be distinguished both *in vitro* and *in silico* based on their dynamic interactions with model membranes, notably their differing conformational flexibility at the N-terminus, ability to form higher order aggregates and the characteristics of induced ion conductance. Taken together, these differences provide an explanation of the greater potency and broader antibacterial spectrum of activity of temporin L over aurein 2.5. Consequently, while the secondary amphipathic, α-helix conformation is a key determinant of the ability of a cationic AMP to penetrate and disrupt the bacterial plasma membrane, the exact mechanism, potency and spectrum of activity is determined by precise structural and dynamic contributions from specific residues in each AMP sequence.

Antimicrobial peptides (AMPs) have remained an effective component of the innate immune system throughout evolutionary history and have been isolated from numerous and diverse organisms<fnref n="1" />⁻⁶. Naturally occurring AMPs have attracted considerable interest as a starting point for rational design processes aiming to enhance antimicrobial or immunomodulatory properties³⁻⁶. Such AMPs are often found to be part of a larger family of very close-related peptides which share substantial primary sequence and physico-chemical similarities and seemingly have very similar functions. The benefits of producing such a range of similar AMPs remain poorly understood and it is not yet clear the extent to which antimicrobial functions are duplicated by such AMPs. Further, since many AMPs that share the same physico-chemical properties and have the same antimicrobial function, have

<footnote n="1">¹Institute of Pharmaceutical Science, School of Cancer & Pharmaceutical Science, King's College London, Franklin-Wilkins Building, 150 Stamford Street, London, SE1 9NH, United Kingdom. ²Technology Development Group, National Infection Service, Public Health England, Salisbury, UK. ³Centre for Biomolecular Spectroscopy and Randall Division of Cell and Molecular Biophysics, King's College London, New Hunt's House, London, SE1 1UL, United Kingdom. ⁴Department of Physics, King's College London, London, WC2R 2LS, United Kingdom. ⁵School of Applied Science, London South Bank University, 103 Borough Road, London, SE1 0AA, United Kingdom. Correspondence and requests for materials should be addressed to A.J.M. (email: james.mason@kcl.ac.uk)</footnote>
<table><caption>**Table 1.** Peptides sequences, biophysical characteristics and concentration of peptide necessary to start membrane activity in electrophysiology experiments. Both peptides were amidated at the C-terminus. Hydrophobicity (H) and hydrophobic moment (µH) were calculated using HeliQuest⁵⁵.</caption><thead><tr><th rowspan="2">Peptide</th><th rowspan="2">Sequence</th><th rowspan="2">Length</th><th rowspan="2">Charge</th><th rowspan="2">H</th><th rowspan="2">µH</th><th colspan="2">Peptide concentration (µM)</th></tr><tr><th>DPhPE/DPhPG</th><th>DPhPG</th></tr></thead><tbody><tr><td>Aurein 2.5</td><td>GLFDIVKKVV GAFGSL</td><td>16</td><td>+2</td><td>0.622</td><td>0.609</td><td>7.69</td><td>5.75</td></tr><tr><td>Temporin L</td><td>FVQWFSKFLG RIL</td><td>13</td><td>+3</td><td>0.906</td><td>0.710</td><td>9.52</td><td>10</td></tr></tbody></table>

been described in unrelated organisms, it is tempting to consider such AMPs as an example of parallel evolution and that the ability of cationic AMPs to form, secondary amphipathic, α-helix conformation is largely sufficient for antimicrobial activity. Of fundamental interest, therefore, is to what extent such related AMPs can be considered alike or, alternatively, in what ways can their activities and functions be distinguished and explained at a molecular level with activity dependent on specific residues in the primary sequence. Although the effect of interactions between AMPs and numerous bacterial structures and process, notably bacterial proteases and efflux systems and bacterial cell walls⁷ will influence susceptibility to AMPs, the ability of AMPs to damage or cross the bacterial plasma membrane is considered to be a major determinant of antibacterial outcomes⁸. Therefore, the mechanism by which AMPs interact with models of such membranes is expected to be particularly informative.

Recently we have considered how relatively minor modifications to the primary sequence of temporin B, an AMP from *Rana temporaria*, can fundamentally alter, not only its potency and spectrum of antibacterial activity, but also its interaction and mode of disruption of membranes designed to mimic the plasma membrane of Gram-positive or Gram-negative bacteria⁹. Interestingly, this study showed that AMPs can be produced that share substantial sequence similarity and that have very similar antibacterial potencies but may nevertheless act with significantly different mechanisms of action. Notably, while biophysical studies conducted in the steady-state struggled to identify any difference in the conformation or membrane interaction of the three temporin B analogues that could be associated with differences in antibacterial activity, time-resolved techniques – notably all-atom MD simulations and patch-clamp – could effectively distinguish the behaviour of the three analogues and provide a rationale for the impact of sequence modification on membrane and antibacterial activities. Previously we used a similar approach to study and compare magainin 2 and pleurocidin. These AMPs have dissimilar primary sequences but share numerous physico-chemical properties, including a preference for α-helix conformation. Nevertheless, substantial differences in the ability of these two AMPs to penetrate model membranes were shown¹⁰. This was associated with a greater conformational flexibility of pleurocidin which is consistent with its ability to access intracellular, bactericidal targets and greater potency¹¹.

Temporin L is a further example of the ten temporin peptides originally isolated from the European red frog *Rana temporaria*¹² and has received considerable attention¹³. It is noticeably more potent than other temporin peptides against both Gram-positive and Gram-negative bacteria, including temporin B, but can act in synergy with at least some of these¹⁴,¹⁵. Its preference for adopting a highly stable α-helix conformation has been described¹⁶ – again more ordered than temporin B¹⁵ – and numerous structure-activity-relationship studies have been performed to enhance its antimicrobial¹⁷⁻²¹ or anti-endotoxin²²⁻²⁵ activities. AMPs have similarly been isolated from Australian frogs; seventeen aurein peptides were found in secretions from the granular dorsal glands of the Green and Golden Bell Frog *Litoria aurea*²⁶. Of these, aurein 2.5 has been studied in some detail²⁷⁻³¹. Aurein 2.5 also has a strong preference for adopting ordered α-helix conformations in a variety of membrane environments²⁸,³⁰,³¹ and strong surface activity was found in both anionic and zwitterionic model membranes²⁷.

*Rana temporaria* and *Litoria aurea* are examples of frogs from, respectively, the Ranidae and Hylidae families. Divergence from their common ancestor occurred during the Jurassic period, approximately 150 million years ago and³², although temporin L and aurein 2.5 are both AMPs derived from frog dorsal secretions¹²,²⁶, they share no sequence similarity. Nevertheless, although temporin L and aurein 2.5 are each members of much larger families of AMP produced by each frog, they represent interesting examples of parallel evolution as, despite the lack of sequence similarity, they do share numerous physico-chemical properties (Table 1); they are both relatively short and hydrophobic, they carry a modest positive charge and are known to adopt ordered α-helix conformations in model membranes or membrane mimicking milieu due to their secondary amphipathicity¹³,²⁸. Here we consider the antibacterial activity of temporin L and aurein 2.5 and investigate, using an approach that combines both time-resolved and steady-state biophysical methods, the extent to which their membrane interactions are similar or can be distinguished and whether such differences can be ascribed to the role of individual residues. The ability to describe the similarities and some key differences, between the two AMPs, at the molecular level, allows better understanding of the importance of specific molecular features – notably conformational flexibility and the role of individual hydrophobic amino acids – refining our appreciation of the requirements for potent antibacterial activity.

## Experimental Procedures

**Peptides and lipids.** Peptides (Table 1) were purchased from Pepceuticals Ltd (Enderby, UK) or Cambridge Research Biochemicals (Cleveland, UK) as desalted grade (crude) and were further purified using water/acetonitrile gradients using a Waters SymmetryPrep C8, 7 mm, 19 × 300 mm column. They were both amidated at the C-terminus. The lipids 1-palmitoyl-2-oleoyl-*sn*-glycero-3-phospho-(1'-*rac*-glycerol) (POPG), 1-palmitoyl*sn*-2-oleoyl-*sn*-glycero-3-phosphoethanolamine (POPE), 1-palmitoyl*sn*-2-oleoyl-*sn*-glycero-3-phosphoethanolamine (POPE-d31),
1,2-diphtyanoyl-*sn*-glycero-3-phospho-(1'-*rac*-glycerol) (DPhPG) and 1,2-diphtyanoyl-*sn*-glycero-3-phosphoethanolamine (DPhPE) were purchased from Avanti Polar Lipids, Inc. (Alabaster, AL) and used without any purification. All other reagents were used as analytical grade or better.

**Antibacterial activity assay.** The antibacterial activity of the peptides was assessed through a modified two-fold broth microdilution assay, as described before¹⁰. Briefly, a two-fold dilution of peptide stock solutions was performed in non-cation adjusted Müller Hinton broth in a 96-well polypropylene microtiter plate. Then, 100 µl of bacterial suspension (back diluted from overnight cultures to an A₆₀₀ of 0.01) were added to 100 µl of peptide solution in each well of the plate. Growth and sterility controls were present for each experiment. The plate was incubated at 37 °C for 18 hours without shaking and the MIC was defined as the lowest concentration where no visible growth could be detected by measuring the absorbance (A₆₀₀) on an Omega plate reader. Strains were selected to look at one susceptible strain and one multidrug resistant isolate of each bacterial species. All strains have been sequenced, to allow linkage of resistance phenotypes to known genetic traits and are readily available through the National Collection of Type Cultures (NCTC).

**NMR structure determination.** Liquid state NMR spectroscopy experiments and data analysis were performed as described previously¹⁰. The samples solution consisted of 2 mM peptide and 100 mM deuterated sodium dodecyl sulphate (SDS-d₂₅) micelles. 10% D₂O containing trimethylsilyl propanoic acid (TSP) was added for the lock signal and as internal chemical shift reference. NMR spectra were acquired on a Bruker Avance 500 MHz spectrometer (Bruker, Coventry, UK) equipped with a cryoprobe with samples maintained at 298 K. Standard Bruker TOCSY and NOESY pulse sequences were used, with water suppression using an excitation sculping sequence with gradients (mlevesgpph and noesyesgpph). The ¹H 90-degree pulse length was 8.0 µs. One TOCSY with a mixing time of 60 ms, and two NOESY with mixing time of 100 and 200 ms were acquired for each sample. The relaxation delay was 1 s. 2048 data points were recorded in the direct dimension, and either 256 or 512 data points in the indirect dimension. CARA (version 1.9.1.2) and Dynamo software³³,³⁴ were used for the structure calculation. Inter-proton NOEs interactions were used as distance restraints in the structure calculation. CARA software generated a total of 200 structures on UNIO'08 (version 1.0.4)³⁵ and XPLOR-NIH (version 2.40)³⁶,³⁷, after 7 iterations, using the simulated annealing protocol. 20 structures with the lowest energy were chosen to produce a final average structure. In the case of ambiguous NOEs assignments by CARA software, Dynamo software's annealing protocol was applied. Only unambiguous NOEs were used in this case after being classified as strong, medium and weak on the base of the relative intensity of the cross-peaks in the NOESY spectra. On the basis of this classification upper limits of 0.27, 0.33 and 0.50 nm have been applied, respectively, as restraint on the corresponding inter-proton distance, as described previously³⁸. One thousand structures were calculated and the 100 conformers with the lowest potential energy were selected for the analysis. The selected 100 conformers were aligned, and the root mean square deviation (RMSD) of the backbone heavy atoms was calculated with respect to their average structure. Solvent molecules were not included in the calculations. Structural coordinates were deposited in the Protein Data Bank (www.rcsb.org) and Biological Magnetic Resonance Bank (BMRB; www.bmrb.wisc.edu) under accession codes 6GS5 and 6GS9 (PDB) and 34283 and 34284 (BMRB) for temporin L and aurein 2.5 respectively. Fingerprint regions of TOCSY and NOESY spectra (Fig. S1) and NOE tables (Fig. S2) are provided as supplementary material.

**Molecular dynamics simulations.** Simulations were carried out on either a Dell Precision quad core T3400 or T3500 workstation with a 1 kW Power supply (PSU) and two NVIDA PNY GeForce GTX570 or GTX580 graphics cards using Gromacs³⁹. The CHARMM36 all-atom force field was used in all simulations⁴⁰,⁴¹. The initial bilayer configuration was built using CHARMM-GUI⁴². All membranes in this project contained a total of 512 lipids, composed either of POPE/POPG (75:25 mol:mol) or POPG to reflect the lipid charge ratios of the plasma membrane of Gram-negative and Gram-positive bacteria, respectively⁴³,⁴⁴. Eight peptides were inserted at least 30 Angstrom above the lipid bilayer in a random position and orientation at least 20 Angstrom apart. The starting structures were taken from the NMR calculation in SDS micelles. The system was solvated with TIP3P water and Na⁺ ions added to neutralize. Energy minimization was carried out at 310 K with the Nose-Hoover thermostat using the steepest descent algorithm until the maximum force was less than 1000.0 kJ/ml/nm (~3000-4000 steps). Equilibration was carried out using the NVT ensemble for 100 ps and then the NPT ensemble for 1000 ps with position restraints on the peptides. Hydrogen-containing bond angles were constrained with the LINCS algorithm. Final simulations were run in the NPT ensemble using 2 femtosecond intervals, with trajectories recorded every 2 picoseconds. All simulations were run for a total of 100 nanoseconds and repeated twice, with peptides inserted at different positions and orientations, giving a total of approximately 0.8 µs simulation. Trajectory RMSD plots (Fig. S3) are provided in the supplementary material. Torsion angles are circular quantities and the circular mean of psi or phi angles may be calculated as follows:

$$ \bar{\psi} = atan2\left(\frac{1}{n}\sum_{j=1}^{n} \sin \psi_j, \frac{1}{n}\sum_{j=1}^{n} \cos \psi_j\right) $$

Similarly, the associated circular variance for psi or phi angles is calculated as follows:

$$ \text{Var}(\psi) = 1 - R_{av} $$

with R being given by:
$$R^2 = \left[ \sum_{i=1}^{n} \cos \psi_i \right]^2 + \left[ \sum_{i=1}^{n} \sin \psi_i \right]^2$$

**Liposome preparation.** Small unilamellar vesicles (SUVs) and multi-lamellar vesicles (MLVs) were prepared for circular dichroism (CD) as previously described¹⁰. Lipid powders were solubilized in chloroform and dried under rotor-evaporation. To completely remove the organic solvent, the lipid films were left overnight under vacuum and hydrate in 5 mM Tris buffer with or without the addition of 100 mM NaCl (pH 7.0). Lipid suspension was subjected to 5 rapid freeze-thaw cycles for further sample homogenisation. POPE/POPG (75:25, mol:mol) and POPG SUVs were obtained by sonicating the lipid suspension on Soniprep 150 (Measuring and Scientific Equipment, London, UK) for 3 × 7 minutes with amplitude of 6 microns in the presence of ice to avoid lipid degradation. The SUVs were stored at 4 °C and used within 5 days of preparation.

**Circular dichroism spectroscopy.** Far-UV spectra of the peptides in the presence of SUVs and SDS micelles were acquired on a Chirascan Plus spectrometer (Applied Photophysics, Leatherhead, UK). Liposome samples were maintained at 310 K. Spectra were recorded from 260 to 190 nm. Lipid suspension was added to a 0.5 mm cuvette at a final concentration of 5.0 mM and then a few µl of a concentrated peptide solution were added and thoroughly mixed to give the indicated final peptide-to-lipid molar ratios. The same experimental conditions were used to investigate peptides secondary structure in SDS micelles. Final peptide concentration in the 0.5 mm cuvette was 40 µM, while SDS micelles concentration was 2 mM (L/P = 50). In processing, a spectrum of the peptide free suspension was subtracted and Savitsky-Golay smoothing with a convolution width of 5 points applied.

**Electrophysiology experiments (Patch-clamp).** Giant unilamellar vesicles (GUVs) composed of DPhPE/DPhPG (60:40, mol:mol) and DPhPG were prepared in the presence of 1 M sorbitol by the electroformation method in an indium-tin oxide (ITO) coated glass chamber connected to the Nanion Vesicle Prep Pro setup (Nanion Technologies GmbH, Munich, Germany) using a 3-V peak-to-peak AC voltage at a frequency of 5 Hz for 120 and 140 minutes, respectively, at 37 °C⁴⁵⁻⁴⁷. Bilayers were formed by adding the GUVs solution to a buffer containing 250 mM KCl, 50 mM MgCl₂ and 10 mM Hepes (pH 7.00) onto an aperture in a borosilicate chip (Port-a-Patch®, Nanion Technologies) and applying 70–90 mbar negative pressure resulting in a solvent-free membrane with a resistance in the GΩ range. Diphytanoyl chains are used here for practical reasons since, unlike lipids with mixed palmitoyl-oleoyl chains such as POPE, these lipids do not undergo the main, temperature dependent transition from disordered fluid into the all *trans* configuration and remain in the same phase between –120° and +120 °C⁴⁸,⁴⁹ while, crucially, the membranes composed of these lipids are mechanically stable and have high specific resistance⁴⁸, essential for electrophysiology experiments. After formation, a small amount of peptide stock solution (in water) was added to 50 µL of buffer solution in order to obtain its active concentration. All the experiments were carried on with a positive holding potential of 50 mV. The active concentration, the concentration at which the peptide first showed membrane activity, for each peptide was obtained through a titration performed in the same conditions. For all the experiments a minimum of 6 according repeats was done. Current traces were recorded at a sampling rate of 50 kHz using an EPC-10 amplifier from HEKA Elektronik (Lambrecht, Germany). The system was computer controlled by the PatchControl™ software (Nanion) and GePulse (Michael Pusch, Genoa, Italy, http://www.ge.cnr.it/ICB/conti_moran_pusch/programs-pusch/software-mik.htm). The data were filtered using the built-in Bessel filter of the EPC-10 at a cut-off frequency of 10 kHz. The experiments were performed at room temperature. Data analysis was performed with the pClamp 10 software package (Axon Instruments). Estimation of pore radii was performed as previously⁵⁰.

## Results

**Temporin L is a more potent antimicrobial than aurein 2.5.** In the present study, temporin L is shown to be more potent than aurein 2.5 against all strains included in both the Gram-positive and Gram-negative bacteria panels and Candida albicans (Table 2). For the Gram-positive strains, aurein 2.5 is outperformed by temporin L with, on average a 7.6-fold greater potency for the latter. For the Gram-negative strains, the difference between the potency of aurein 2.5 and temporin L is less, notably for both Acinetobacter baumannii isolates. However, the spectrum of activity for aurein 2.5 is narrower with almost no detectable activity against either Pseudomonas aeruginosa isolate suggesting the two peptides have differing abilities to penetrate its outer membrane, the potency of each AMP is likely affected by interactions with the differing cell walls in each pathogen and/or may differ in their susceptibilities to bacterial proteases and/or efflux systems. Nevertheless, if the antibacterial activity of these two AMPs is dependent on their membrane interactions, there is reason to expect differing interactions and behaviours of each AMP with model membranes designed to mimic the differing plasma membranes of these two groups.

**Both temporin L and aurein 2.5 adopt highly ordered α-helix conformations in membrane mimicking media.** Far-UV circular dichroism spectra were obtained in anionic micelles (Fig. 1C) but also in small unilamellar vesicles (SUVs) whose composition was designed to mimic the plasma membrane of Gram-positive or Gram-negative bacteria, and in aqueous solution (Fig. S4). The structure of temporin L has previously been determined both in SDS and dodecylphosphocholine (DPC) micelles²⁰ but coordinates are not available in the PDB, while no structure has been determined for aurein 2.5. Structures of both peptides were therefore also solved in the presence of SDS (Fig. 1A,B) and were used subsequently as starting structures for the molecular dynamics simulations. In SDS, both the far-UV CD and the NMR structures agree that both peptides adopt ordered α-helix conformations. The NMR structures and CD spectra can be compared with those obtained
<table><caption>**Table 2.** Antimicrobial activity. MS – methicillin sensitive; EMR – epidemic methicillin resistant; VS – vancomycin sensitive; VR – vancomycin resistant.</caption><thead><tr><th rowspan="2"></th><th rowspan="2">Isolate</th><th colspan="2">Peptide concentration (µg/ml)</th></tr><tr><th>Aurein 2.5</th><th>Temporin L</th></tr></thead><tbody><tr><td rowspan="7">Gram-negative</td><td>Klebsiella pneumoniae NCTC 13368</td><td>64</td><td>16</td></tr><tr><td>Klebsiella pneumoniae M6</td><td>32–64</td><td>16</td></tr><tr><td>Acinetobacter baumanii AYE</td><td>4–8</td><td>4</td></tr><tr><td>Acinetobacter baumanii ATCC 17978</td><td>8</td><td>4</td></tr><tr><td>Pseudomonas aeruginosa PAO1</td><td>128</td><td>16</td></tr><tr><td>Pseudomonas aeruginosa NCTC 13437</td><td>128</td><td>32–64</td></tr><tr><td>Escherichia coli NCTC 12923</td><td>16</td><td>4</td></tr><tr><td rowspan="5">Gram-positive</td><td>MS Staphylococcus aureus ATCC 9144</td><td>8</td><td>2</td></tr><tr><td>EMR Staphylococcus aureus-15</td><td>32</td><td>4</td></tr><tr><td>EMR Staphylococcus aureus-16</td><td>32–64</td><td>4</td></tr><tr><td>VS Enterococcus faecalis NCTC 775</td><td>32–64</td><td>4–8</td></tr><tr><td>VR Enterococcus faecalis NCTC 12201</td><td>64</td><td>8</td></tr><tr><td>Yeast</td><td>Candida albicans NCPF 3179</td><td>16</td><td>8</td></tr></tbody></table>

for temporin B and its analogues which describe peptides with a much greater degree of disorder⁹. The α-helix of aurein 2.5 extends almost throughout the length of the peptide while, for temporin L a little disorder is observable at the N-terminus only.

Both peptides display characteristic secondary amphipathicity when adopting this ordered conformation. After binding to lipid bilayers, in the steady-state far-UV CD again indicates similar ordered α-helix conformations are adopted by both peptides; in the absence of lipid both peptides adopt disordered conformations (Fig. S4). There are no apparent differences in conformation between either peptide in SUVs designed to mimic Gram-negative (Fig. S4A,B) or Gram-positive (Fig. S4C,D) plasma membranes and there was no significant effect of varying the peptide to lipid ratio.

**MD simulations reveal N-terminus conformational flexibility in temporin L but not aurein 2.5 on initial binding to model membranes.** Conformational flexibility has proven to be a key parameter capable of distinguishing the behaviour of AMPs with models of both Gram-negative and Gram-positive bacteria plasma membranes⁹,¹⁰. Therefore all-atom MD simulations were used to characterise the atomistic details of the interaction and conformational flexibility of aurein 2.5 and temporin L with model bacterial membranes. For each simulation, the time average of psi dihedral angles for each residue in the eight peptides is calculated, along with the circular variance as a measure of the spread of the dihedral angles⁵¹, over the duration of the simulation. The circular variance varies from 0 to 1 with a low value equating to low conformational flexibility and vice versa. For aurein 2.5 the average psi dihedral angles are similar in the presence of both POPE/POPG (Fig. 2A) and POPG (Fig. 2B) membranes and lie in the range expected for α-helix conformation. With the exception of Ser15, this is maintained throughout the length of the peptide. As evidenced by the relatively low values obtained for the circular variance throughout the length of the peptide, conformational flexibility is absent from the N-terminus but rises from Val10 towards the C-terminus. As in previous work⁹,¹⁰, phi dihedral angles remain relatively constant throughout the length of the peptide and, apart from the two glycine residues (Gly11, Gly 14), variance in phi was observed only for Ser15 (Fig. S5).

In the analogous simulations performed for temporin L however, features are revealed that now distinguish its behaviour in model membranes from aurein 2.5 (Fig. 2C,D). In both POPE/POPG (Fig. 2C) and POPG (Fig. 2D) the α-helix conformation is retained for the first ten (of thirteen) residues with the three residues that follow Gly10 (Arg11, Ile12 and Leu13) possibly adopting more extended conformations. Similar to aurein 2.5, the circular variance indicates that the four residues at the C-terminus experience substantial conformational flexibility. However, in contrast with aurein 2.5, further conformational flexibility is also apparent at the N-terminus with psi angles between the first three or four residues subject to substantial variance. Notably, between Trp4 and Leu9, there is a section of highly ordered α-helix (WFSKFL) with almost no variation in psi dihedral angle detected. There is little apparent difference between the secondary structure adopted by temporin L during the initial binding to either POPE/POPG or POPG bilayers apart from perhaps slightly greater conformational flexibility at the N-terminus when binding to POPG (Fig. 2D). Again, the phi dihedral angles are relatively constant throughout the length of the peptide, with variance in phi observed only for Gly10 and Arg11 (Fig. S5).

The topology of temporin B during the initial binding to lipid bilayers is sensitive to changes in primary sequence⁹. Here the ability of aurein 2.5 or temporin L to penetrate POPE/POPG or POPG bilayers can be visualized (Figs 3 and S6). The two peptides both attempt to penetrate either POPE/POPG or POPG bilayers via the N-terminus. For aurein 2.5, although the differing bilayer composition ensures the peptide penetrates deeper and more rapidly in POPG when compared with POPE/POPG, in both cases the first six to nine residues were consistently able to penetrate below the plane of the lipid phosphate groups (Fig. 3A,B). In contrast, temporin L penetrates POPE/POPG and POPG bilayers in subtly different ways (Fig. 3C,D). In both cases, penetration proceeds via the N-terminus. When binding to POPE/POPG, the first eight residues of temporin L penetrate to
<figure><img src="image_2.png" /><figcaption>**Figure 1.** Three-dimensional structure for aurein 2.5 and temporin L. Structures were determined through ¹H-NOESY NMR spectroscopy in SDS-d25 micelles (peptide/detergent equal to 1/50). The figure shows the structure with the minimum RMSD, used as starting point in the MD simulations. The hydrophobic and hydrophilic residues are shown in yellow and blue for aurein 2.5 (A) and temporin L (B). CD spectra of the peptides in the same conditions used for NMR experiments (C).</figcaption></figure>

the plane of the lipid phosphate and a certain periodicity, consistent with the ordered α-helix conformation, can be observed indicating an alignment roughly parallel to the bilayer surface (Fig. 3C). When binding to POPG however, a deeper penetration by the first four residues at the N-terminus is observed. Residues 5–8 (FSKF) penetrate more slowly and the periodicity is absent, indicating a steeper angle of attack on the bilayer (Fig. 3D) and reflecting the greater conformational flexibility in the first four residues (Fig. 2D).

**Hydrophobic amino acids mediate peptide aggregation at the bilayer surface.** Starting as randomly positioned monomers positioned at least 20 Å apart, both aurein 2.5 (Fig. 4) and temporin L (Fig. 5) were frequently observed to self-associate during the simulations (Fig. S7). In both POPE/POPG and POPG, aurein 2.5 was observed as a monomer and in dimer or trimer self-associations (Fig. 4B,E). For both dimer (Fig. 4A) and trimer (Fig. 4D) configurations, the self-association was mediated by hydrophobic groups, notably Phe3 and Phe13 which are located near to N- and C-termini respectively and the three valines that are located in the centre of the peptide (Fig. 4C,F). This configuration was qualitatively the same when inserting into both POPE/POPG and POPG membranes and was reproduced in both duplicates.

Temporin L also adopts oligomeric configurations in both bilayers (Fig. 5). In POPE/POPG, again populations of monomer, dimer and trimer aggregates are observed (Fig. 5B) with self-association mediated by contacts between hydrophobic residues in the centre of the primary sequence (Phe5, Phe8, Leu9). Again, this was reproduced in both duplicates. In POPG, while again both monomers and dimers are observed, a higher order aggregate was detected in one of the duplicated simulations. This comprised four (Fig. 5D) or sometimes even five peptides (Fig. 5E) with a series of aligned, hydrophobic, phenylalanine residues (Phe1, Phe5, Phe8) found to be mediating the self-association (Fig. 5F). As can be seen from the snapshot, within this aggregate, peptides are oriented with the α-helix long axis oriented both tilted towards and perpendicular to the bilayer surface. In the duplicate simulation, monomers and one trimer were observed, with self-association in that aggregate mediated by all residues between Phe8 and the C-terminus (FLGRIL).

**Aurein 2.5 and temporin L induce channel like activity in model membranes.** As in our previous study⁹, we made use of the Port-a-patch® automated patch-clamp system from Nanion Technologies (Munich, Germany) to determine whether the two peptides could be distinguished on the basis of their disruption of either DPhPE:DPhPG (60:40 mol:mol) or DPhPG bilayers, mimicking Gram-negative and Gram-positive bacteria cytoplasmic membranes, respectively⁴³,⁴⁴. For the temporin B analogues the technique determined substantial
<figure><img src="image_3.png" /><figcaption>**Figure 2.** Secondary structure analysis of aurein 2.5 and temporin L peptides from MD simulations. Mean dihedral angles (psi) (upper plot) and their circular variance (lower plot) are shown for each residue to enable identification of regiospecific changes in conformation and/or flexibility. Angles were averaged over 100 ns of simulation and eight aurein 2.5 (A,B) or temporin L (C,D) peptides when binding to POPE/POPG (A,C) or POPG (B,D) membranes. Backbone dihedral angles (psi) for residues in α-helix conformation are around -45°.</figcaption></figure>

<figure><img src="image_4.png" /><figcaption>**Figure 3.** Average aurein 2.5 and temporin L peptide topology during 100 ns MD simulation. The depth of aurein 2.5 (A,B) and temporin L (C,D) insertion into each membrane is shown as the Z-position for each residue, averaged over all eight peptides in each simulation, relative to the phosphate group plane in the upper bilayer leaflet. Positive or negative values indicate the peptides are below or above the phosphate group, respectively for peptides in POPE/POPG (A,C) or POPG (B,D).</figcaption></figure>
<figure><img src="image_5.png" /><figcaption>**Figure 4.** Interaction of aurein 2.5 with lipid bilayers *in silico*. Examples of peptides aggregating as a dimer or a trimer are shown respectively for POPE/POPG (A) and POPG (D) simulations. Basic, polar, acidic and hydrophobic residues are coloured in blue, green, red and white respectively with lipids and water removed for clarity. The top view snapshot shows the tendency to aggregate for the eight peptides in POPE/POPG (B) or POPG (E). Heatmaps of average inter-peptide $C_α-C_α$ distances, in Angstroms, for peptides in POPE/POPG (C) or POPG (F) show the stronger interactions that mediate assembly of aggregates.</figcaption></figure>

variation in the concentration of peptide required to induce conductance, the latency – the time taken for conductance to commence after addition of peptide –, the duration of activity and whether the membrane was broken through the action of the AMP. In no cases was any characteristic channel like activity – well-defined events with discrete opening levels – observed.

In contrast, in the present study both aurein 2.5 and temporin L were found to be capable of inducing channel like activity (Fig. 6). For aurein 2.5 this was observed in both membrane types (Fig. 6A,B) whereas for temporin L this was only found when challenging membranes formed from DPhPG (Fig. 6D). Similar concentrations of aurein 2.5 and temporin L were required to induce channel activity in both membranes (Table 1) and this parameter therefore does not explain any of the observed differences in antibacterial potency. Membrane activity was observed to commence a little earlier for aurein 2.5 in both membranes ($p < 0.05$) but then the activity continued for both peptides for some minutes before the membrane ruptured (Fig. S8). For aurein 2.5 the measurements in both DPhPE/DPhPG and DPhPG membranes adopted a similar pattern with very many low amplitude events and few higher amplitude events (Fig. S9). Opening levels varied within and between replicate traces indicating that pores that form are both transient and irregular. Nevertheless, a set of discrete levels can be defined based on the frequency of certain amplitude events and this enables a measure of the conductance – a function of dwell time – which in turn can be used to estimate the pore size (Table 3). Three discrete levels were consistently identified for aurein 2.5 in both membranes with a slightly higher amplitude and conductance found in DPhPG membranes. The pore size estimate suggests pores approximately equal to the size of a chloride ion are formed by the peptide.

Similar to aurein 2.5, temporin L has channel like activity when challenging DPhPG membranes. For this peptide only two levels could be consistently detected but the amplitude and conductance associated with the higher of the two levels were substantially greater than that found and calculated for aurein 2.5 in the same membrane (Table 3). The resulting estimated pore size indicates that temporin L is capable of creating pores that are at least twice the diameter of a chloride ion. When challenging DPhPE/DPhPG, temporin L induces numerous high amplitude conductance events (Fig. 6C). However, these events are not well-defined or discrete. Instead, bursts of current (fast events seen as spikes) with a wide distribution of amplitudes are detected. This is similar to the activity observed for temporin B and its analogues⁹, albeit with many events of much greater amplitude.
<figure><img src="image_6.png" /><figcaption>**Figure 5.** Interaction of temporin L with lipid bilayers *in silico*. Examples of peptides aggregating as a trimer or a tetramer are shown respectively for POPE/POPG (A) and POPG (D). Basic, polar, acidic and hydrophobic residues are coloured in blue, green, red and white respectively with lipids and water removed for clarity. The top view snapshot shows the tendency to aggregate for the eight peptides in POPE/POPG (B) or POPG (E). Heatmaps of average inter-peptide $C_α-C_α$ distances, in Angstroms, for peptides in POPE/POPG (C) or POPG (F) show the stronger interactions that mediate assembly of aggregates.</figcaption></figure>

<figure><img src="image_7.png" /><figcaption>**Figure 6.** Membrane activity of the temporin B peptides. Representative current traces and all-points histograms for aurein 2.5 (A,B) and temporin L (C,D) in DPhPE/DPhPG (A,C) or DPhPG (B,D) membranes. Experiments were acquired with a holding potential of +50 mV and in the presence of 250 mM KCl, 50 mM MgCl₂, 10 mM HEPES, pH 7.0. y-axis scaled accordingly to the different amplitudes for the different peptides.</figcaption></figure>
<table><caption>**Table 3.** Summary of channel-like activity detected at various opening levels. Pore size is estimated as previously⁵⁰.</caption><thead><tr><th rowspan="2">Peptide</th><th rowspan="2">Parameter</th><th colspan="3">Level</th></tr><tr><th>Level 1</th><th>Level 2</th><th>Level 3</th></tr></thead><tbody><tr><td>Aurein 2.5 (DPhPE/DPhPG)</td><td rowspan="3">Amplitude (pA)</td><td>0.58 ± 0.03</td><td>1.29 ± 0.07</td><td>3.21 ± 0.17</td></tr><tr><td>Aurein 2.5 (DPhPG)</td><td>0.83 ± 0.03</td><td>2.46 ± 0.06</td><td>5.43 ± 0.04</td></tr><tr><td>Temporin L (DPhPG)</td><td>0.89 ± 0.02</td><td><strong>25.38 ± 0.14</strong></td><td>—</td></tr><tr><td>Aurein 2.5 (DPhPE/DPhPG)</td><td rowspan="3">Conductance (pS)</td><td>11.54 ± 0.52</td><td>25.89 ± 1.36</td><td>64.24 ± 3.38</td></tr><tr><td>Aurein 2.5 (DPhPG)</td><td>16.8 ± 0.67</td><td>49.22 ± 1.30</td><td>108.56 ± 0.76</td></tr><tr><td>Temporin L (DPhPG)</td><td>17.89 ± 0.39</td><td><strong>507.54 ± 2.81</strong></td><td>—</td></tr><tr><td>Aurein 2.5 (DPhPE/DPhPG)</td><td rowspan="3">Estimated pore radius (nm)</td><td>0.06 ± 0.01</td><td>0.09 ± 0.02</td><td>0.14 ± 0.03</td></tr><tr><td>Aurein 2.5 (DPhPG)</td><td>0.07 ± 0.01</td><td>0.13 ± 0.02</td><td>0.19 ± 0.01</td></tr><tr><td>Temporin L (DPhPG)</td><td>0.08 ± 0.01</td><td><strong>0.43 ± 0.03</strong></td><td>—</td></tr></tbody></table>

## Discussion

**Time-resolved techniques distinguish mechanistic features between the frog antimicrobial peptides aurein 2.5 and temporin L.** Here we investigated whether two short, hydrophobic antimicrobial peptides, both derived from frog dorsal secretions and either known or predicted to adopt highly ordered α-helix conformations are identical or can be distinguished in their interactions with model membranes and whether such behaviour can be ascribed to specific amino acids in the primary sequences. This was done to better understand differences in their antimicrobial potencies and consequently the extent to which differences in sequence or similarities in physico-chemical properties in AMPs from two distantly related frog species determine antibacterial outcomes. Both aurein 2.5 and temporin L were confirmed as adopting highly ordered α-helix conformations in SDS, by CD and NMR, but also in more realistic membrane environments by CD. Both peptides were also capable of inducing channel-like conductance in model membranes and this should be contrasted with the membrane activity of temporin B⁹, a similarly short and hydrophobic AMP that nevertheless adopts much less ordered α-helix conformation, which is incapable of inducing channel-like activity in the exact same conditions used for the present study. More detailed analysis of the channel conductance measurements and comparison of the results obtained from MD simulations however allow key differences in behaviour of the peptides to be revealed. This suggests, not only that the two peptides operate with distinct mechanisms against the same membrane targets, but also that temporin L may adopt differing strategies for disrupting the plasma membranes of Gram-negative and Gram-positive bacteria. With the caveat that other interactions with the cell wall, bacterial proteases and efflux systems will influence activity, the qualitative differences in membrane activity may explain the observed differences in potency of the two AMPs against Gram-positive and Gram-negative bacteria and can be linked to the subtly differing roles of hydrophobic amino acids in each AMP. These two time-resolved techniques therefore appear useful in defining features of these AMPs that determine their potency. It is notable that the timescale of between 10 and 20 minutes for the appearance and duration of channel conductance is consistent with that observed for the attack of AMPs on bacteria reported elsewhere⁵²,⁵³. Although recent MD simulation studies of the behaviour of multiple AMPs binding to membranes in systems of comparable size explore up to 800 ns⁵⁴, it is currently not possible to sample such timescales using all-atom molecular dynamics simulations. The simulations here describe the initial binding of the two peptides to the model bilayers - 100 ns is sufficient to sample a period in which a local minimum is attained⁵⁴. Notably, in the case of temporin L, the importance of the conformational flexibility found in the N- and C-terminus regions and the role of key residues in the formation of higher order aggregates, described in the present study, is supported by a substantial body of previous work. We therefore consider, below, the extent to which differences in antimicrobial and membrane activities can be related to the apparent function of individual elements of the primary sequence whose behaviour has been partly revealed here at the molecular level.

**The high relative potency of temporin L against Gram-positive bacteria is associated with larger pores and higher conductance due to the formation of higher order aggregates.** Interpretation of the findings from the MD simulations is substantially aided by previous structure-activity relationship and biophysical studies, most notably of temporin L¹⁸⁻²¹,²³⁻²⁵. In particular, substitutions for Gln3, Gly10 and various replacements of the three phenylalanine residues have been made and their effect on antibacterial potency, cell cytotoxicity, haemolysis and anti-endotoxin activity determined¹⁸⁻²¹,²³⁻²⁵.

The C-terminal segment, Gly10-Leu13, and the N-terminus to Gln3 are shown here to experience substantial conformational flexibility. The C-terminal segment penetrates less into the hydrophobic core of model membranes than the remainder of the peptide while Gln3 is less buried when binding to POPE/POPG bilayers. The C-terminal segment and Gln3 were the most exposed to solvent in SDS micelles in a previous NMR study²⁰ hence there is agreement that temporin L is located at the interfacial region of the membrane with these sections of the sequence having specific roles. Further support for this comes from considering the previously reported effects of substituting Gly10 for either leucine or proline¹⁸. Replacement of Gly10 with leucine substantially increases the propensity to form α-helix and leads to abrogation of the antibacterial activity against both Gram-positive and Gram-negative species and drastically increased haemolytic potential. In contrast, replacement of Gly10 with proline produced a peptide with negligible haemolytic potential but at the costs of a significant reduction
in antibacterial activity. The conclusion that helical conformation at the C-terminus is essential for antibacterial activity¹⁸ can be refined to state that conformational flexibility in the C-terminal region is essential for antibacterial activity and this is afforded by glycine and not proline. As seen previously¹⁰ for magainin 2, glycine alone is insufficient to induce conformational flexibility in an AMP sequence and the residues located close to Gly10 can be expected to contribute.

Gln3 has also been replaced with proline¹⁸,²⁰, which has little effect on the antimicrobial activity but reduces a little the haemolytic potential, and lysine²³,²⁵, which increases the antiendotoxin capability of the peptide and produced a modest improvement in antibacterial activity. In these analogues the first two, hydrophobic residues were not altered, and the retention of activity is consistent with the essential role, shown here, for these residues in the initial binding as these penetrate the most into the hydrophobic core. Increasing the flexibility around the third residue may alter binding to membranes of differing composition while increasing the cationicity may enhance binding to the anionic surface of bacterial plasma membranes without compromising insertion; as was observed for analogues of temporin B⁹.

The importance of an identified phenylalanine zipper motif²³ can also now be considered in the context of its role, shown here, in mediating contacts between self-associated peptides in dimers and higher order aggregates. The hydrophobic residues, in particular the three phenylalanines, are involved in the initial binding to the membrane and the stabilisation of peptide aggregates in the membrane. Phe5 and Phe8 have been substituted for either alanines or leucines²⁴. Substitution of either residue with alanine substantially reduces antibacterial activity, cytotoxicity and haemolytic potential with substitution of Phe5 having greater impact. Substituting leucine for either Phe5 or Phe8 has little impact on antibacterial activity or cytotoxicity but reduces a little the haemolytic potential. Interestingly however, substitution of leucine for both Phe5 and Phe8 leads to total loss of antibacterial activity against S. aureus 25923 while activity against E. coli 25922 and P. aeruginosa BAA 427 is retained; haemolysis and cytotoxicity are increased. This indicates that both Phe5 and Phe8 are required for activity against Gram-positive S. aureus and that, combined, they have a specific role when temporin L challenges this, and related, species which cannot be substituted by other hydrophobic residues. In the present study, higher order aggregates and pores with diameters greater than twice the size of a chloride ion were observed only in membranes designed to model the Gram-positive plasma membrane i.e. uniformly negatively charged. Further, in membranes designed to model the Gram-negative plasma membrane, neither higher order aggregates nor regular channel like activity was observed. Therefore, although the phenylalanine zipper might play a role in determining the outcome of interactions between temporin L and the cell wall or other important events, we can hypothesise that in Gram-positive plasma membranes only, higher order temporin L aggregates can form and these lead to relatively larger channels and this underpins high antibacterial potency against such species. By extension, temporin L uses a different mechanism of permeabilising the plasma membrane of Gram-negative bacteria, which does not require the formation of higher order aggregates or regular channels, but which nevertheless induces very high amplitude conductance.

**Aurein 2.5 forms low order aggregates and acts via low conductance channel like permeabilisation of both Gram-negative and Gram-positive plasma membranes.** In comparison with temporin L, less is known about the role of individual residues in aurein 2.5. Further, although both peptides are confirmed as preferring ordered α-helix conformations, aurein 2.5 behaves differently to temporin L in models of both Gram-positive and Gram-negative membranes. Hydrophobic residues, which include two phenylalanines but also three valines and two leucines, mediate the formation of dimers or trimers but not higher order aggregates. Channels that are formed are small and induce relatively low conductance.

There are two glycines in the primary sequence of aurein 2.5, Gly14 and Gly11 which might be expected to promote conformational flexibility. Our recent comparison of conformational flexibility in pleurocidin and magainin 2 found substantial difference in conformational flexibility between the two peptides despite very similar positioning of glycine residues¹⁰ and, taking both studies together, we can conclude that this important property will be dependent both on positioning of the glycines themselves but also on contributions from residues in their vicinity. Gly11 is positioned to interrupt a long section of α-helix and, while preference for α-helix conformation is retained, conformational flexibility increases gradually from this residue towards the C-terminus until Gly14, where conformational flexibility is enhanced still further. In contrast with temporin L, where conformational flexibility is observed in segments located at both the N- or C-terminus, conformational flexibility is absent from the N-terminus in aurein 2.5. Thus, while temporin L differs in the way the N-terminus mediates penetration of the bilayer in models of Gram-negative and Gram-positive plasma membranes, aurein 2.5 penetrates both models in a similar fashion, albeit more readily into the latter.

The hydrophobic phenylalanines that are distributed in a heptad zipper motif and are responsible for mediating peptide self-association in temporin L are, in aurein 2.5, located close to either the N- or C-terminus. Instead, hydrophobic valine residues, Val 6, Val 9 and Val 10, are distributed throughout the aurein 2.5 primary sequence and perform the role of mediating self-association. The selection of valine over phenylalanine in these positions in aurein 2.5, coupled with the topology and/or geometry of the AMP, precludes the formation of higher order aggregates and highlights a specific function that can be fulfilled by phenylalanine but not other hydrophobic amino acids.

Taken together, the absence of conformational flexibility in the N-terminus and/or phenylalanines positioned to stabilise the assembly of higher order aggregates ensures there are fundamental differences in how aurein 2.5 interacts with models of bacterial plasma membranes, when compared with temporin L, despite the two cationic AMPs adopting very similar secondary amphipathic, α-helix conformations. The selection of an ability to form low or higher order aggregates may well impact on the ability to cause high conductance membrane permeabilisation but it may impact other parameters. Notably aurein 2.5 induces channel like activity much more rapidly than temporin L in DPhPG membranes and a little less peptide is needed. Future work may explore whether indeed higher order aggregates are needed for high conductance channel like activity and whether this is associated with a slower onset of channel activity.
Conclusion

Two AMPs from distantly related frogs were found to adopt seemingly identical, ordered α-helix conformations in model membranes. MD simulations and conductance measurements indicated however that there are fundamental differences in how the two AMPs bind to and disrupt these membranes. In particular, three phenylalanines in temporin L, which cannot be substituted by other hydrophobic amino acids, are suggested to be the origin of its ability to form higher order aggregates and relatively large ion conducting channels in models of Gram-positive bacteria plasma membranes and hence adopt a mechanism of action against such species that is distinct from that used against Gram-negative species. This shows that, while secondary amphipathic α-helix conformations are a vital element in the binding and disruption of bacterial plasma membranes, individual amino acids in the distinct structures make important structural and dynamic contributions that determine the precise mechanism of action, potency and species selectivity of AMPs.

References

1. Andersson, D. I., Hughes, D. & Kubicek-Sutherland, J. Z. Mechanisms and consequences of bacterial resistance to antimicrobial peptides. *Drug Resist. Update.* **26**, 43–57 (2016).

2. Zasloff, M. Antimicrobial peptides of multicellular organisms. *Nature*. **415**, 389–395 (2002).

3. Giuliani, A. et al. Antimicrobial peptides: natural templates for synthetic membrane-active compounds. *Cell. Mol. Life Sci.* **65**, 2450–2460 (2008).

4. Hancock, R. E. & Diamond, G. The role of cationic antimicrobial peptides in innate host defences. *Trends Microbiol.* **8**, 402–410 (2000).

5. Haney, E. F. & Hancock, R. E. W. Peptide design for antimicrobial and immunomodulatory applications. *Biopolymers.* **100**, 572–583 (2013).

6. Jenssen, H., Hamill, P. & Hancock, R. E. W. Peptide antimicrobial agents. *Clin. Microbiol. Rev.* **19**, 491–511 (2006).

7. Malanovic, N. & Lohner, K. Gram-positive bacterial cell envelopes: The impact on the activity of antimicrobial peptides. *Biochim. Biophys. Acta* **1858**, 936–946 (2016).

8. Nguyen, L. T., Haney, E. F. & Vogel, H. J. The expanding scope of antimicrobial peptide structures and their modes of action. *Trends Biotechnol.* **29**, 464–472 (2011).

9. Manzo, G. et al. Minor sequence modifications in temporin B cause drastic changes in antibacterial potency and selectivity by fundamentally altering membrane activity. *Scientific Reports*, https://doi.org/10.1038/s41598-018-37630-3 (2019).

10. Amos, S.-B. T. *et al.* Antimicrobial peptide potency is facilitated by greater conformational flexibility when binding to Gram-negative bacterial inner membranes. *Sci. Rep.* **6**, 37639–37651 (2016).

11. Kozłowska, J. et al. Combined systems approaches reveal highly plastic responses to antimicrobial peptide challenge in *Escherichia coli*. *PLoS Pathog* **10**(5), e1004104 (2014).

12. Simmaco, M. et al. Temporins, antimicrobial peptides from the European red frog *Rana temporaria*. *Eur. J. Biochem.* **242**, 788–792 (1996).

13. Mangoni, M. L., Di Grazia, A., Cappiello, F., Casciaro, B. & Luca, V. Naturally occurring peptides from *Rana temporaria*: Antimicrobial properties and more. *Curr. Top. Med. Chem.* **16**, 54–64 (2016).

14. Rinaldi, A. C. et al. Temporin L: antimicrobial, haemolytic and cytotoxic activities, and effects on membrane permeabilization in lipid vesicles. *Biochem. J.* **368**, 91–100 (2002).

15. Rosenfeld, Y., Barra, M., Simmaco, M., Shai, Y. & Mangoni, M. L. A synergism between temporins towards Gram-negative bacteria overcomes resistance imposed by the lipopolysaccharide protective layer. *J. Biol. Chem.* **281**, 28565–28574 (2006).

16. D'Abramo, M. et al. Conformational behaviour of temporin A and temporin L in aqueous solution: a computational/experimental study. *Biopolymers* **81**, 215–224 (2006).

17. Mangoni, M. L. et al. Effects of the antimicrobial peptide temporin L on cell morphology, membrane permeability and viability of *Escherichia coli*. *Biochem. J.* **380**, 859–865 (2004).

18. Mangoni, M. L. et al. Structure-activity relationship, conformational and biological studies of temporin L analogues. *J. Med. Chem.* **54**, 1298–1307 (2011).

19. Merlino, F. et al. Glycine-replaced derivatives of [Pro³,DLeu⁹]TL, a temporin L analogue: Evaluation of antimicrobial, cytotoxic and hemolytic activities. *Eur. J. Med. Chem.* **139**, 750–761 (2017).

20. Carotenuto, A. et al. A different molecular mechanism underlying antimicrobial and hemolytic actions of temporins A and L. *J. Med. Chem.* **51**, 2354–2362 (2008).

21. Saviello, M. R. et al. New insight into the mechanism of action of the temporin antimicrobial peptides. *Biochemistry* **49**, 1477–1485 (2010).

22. Giacometti, A. et al. Interaction of antimicrobial peptide temporin L with lipopolysaccharide *in vitro* and in experimental rat models of septic shock caused by Gram-negative bacteria. *Antimicrob. Agents Chemother.* **50**, 2478–2486 (2006).

23. Srivastava, S. & Ghosh, J. K. Introduction of a lysine residue promotes aggregation of temporin L in lipopolysaccharides and augmentation of its antiendotoxin property. *Antimicrob. Agents Chemother.* **57**, 2457–2466 (2013).

24. Srivastava, S., Kumar, A., Tripathi, A. K., Tandon, A. & Ghosh, J. K. Modulation of anti-endotoxin property of temporin L by minor amino acid substitution in identified phenylalanine zipper sequence. *Biochem. J.* **473**, 4045–4062 (2016).

25. Farrotti, A. et al. Molecular dynamics simulations of the host defense peptide temporin L and its Q3K derivative: an atomic level view from aggregation in water to bilayer perturbation. *Molecules* **22**, E1235 (2017).

26. Rozek, T. et al. The antibiotic and anticancer active aurein peptides from the Australian Bell Frogs *Litoria aurea* and *Litoria raniformis*. *Eur. J. Biochem.* **267**, 5330–5341 (2000).

27. Dennison, S. R., Morton, L. H. G., Shorrocks, A. J., Harris, F. H. & Phoenix, D. A. A study on the interactions of Aurein 2.5 with bacterial membranes. *Colloids Surf. B* *Biointerfaces* **68**, 225–230 (2009).

28. Dennison, S. R., Morton, L. H. G. & Phoenix, D. A. Role of molecular architecture on the relative efficacy of aurein 2.5 and modelin 5. *Biochim. Biophys. Acta* **1818**, 2094–2102 (2012).

29. Dennison, S. R., Harris, F., Morton, L. H. & Phoenix, D. A. Antimicrobial activity of aurein 2.5 against yeasts. *FEMS Microbiol. Lett.* **346**, 140–145 (2013).

30. Dennison, S. R., Morton, L. H., Harris, F. & Phoenix, D. A. The interaction of aurein 2.5 with fungal membranes. *Eur. Biophys. J.* **43**, 255–264 (2014).

31. Wang, J. et al. The cooperative behaviour of antimicrobial peptides in model membranes. *Biochim. Biophys. Acta.* **1838**, 2870–2881 (2014).

32. Zhang, P. et al. Efficient sequencing of Anuran mtDNAs and a mitogenomic exploration of the phylogeny and evolution of frogs. *Mol. Biol. Evol.* **30**, 1899–1915 (2013).

33. Keller, R. *The Computer Aided Resonance Assignment Tutorial*. (Cantana Verlag, Zurich, 2004).

34. Dynamo software: The NMR Molecular Dynamics and Analysis System. Available online at: http://spin.niddk.nih.gov/NMRPipe/dynamo.
35. Guerry, P. & Herrmann, T. Comprehensive Automation for NMR Structure Determination of Proteins. *Methods in Molecular Biology* **831**, 429–451 (2012).

36. Schwieters, C., Kuszewski, J., Tjandra, N. & Clore, G. M. The Xplor-NIH NMR molecular structure determination package. *Journal of Magnetic Resonance* **160**(1), 65–73 (2003).

37. Schwieters, C., Kuszewski, J. & Clore, G. M. Using Xplor-NIH for NMR molecular structure determination. *Prog. Nucl. Magn. Reson. Spectrosc.* **48**, 47–62 (2006).

38. Manzo, G. et al. The singular behavior of a β-type semi-synthetic two branched polypeptide: three-dimensional structure and mode of action. *Phys. Chem. Chem. Phys.* **18**, 30998–31011 (2016).

39. Abraham, M. J. et al. GROMACS: High performance molecular simulations through multi-level parallelism from laptops to supercomputers. *SoftwareX* **1**, 19–25 (2015).

40. Best, R. B. et al. Optimization of the additive CHARMM all-atom protein force field targeting improved sampling of the backbone φ, ψ and side-chain χ1 and χ2 dihedral angles. *J. Chem. Theory Comput.* **8**(9), 3257–3273 (2012).

41. Huang, J. & MacKerell, A. D. CHARMM36 all-atom additive protein force field: validation based on comparison to NMR data. *J. Comput. Chem.* **34**, 2135–2145 (2013).

42. Lee, J. 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.* **12**(1), 405–413 (2016).

43. Epand, R. M. & Epand, R. F. Bacterial membrane lipids in the action of antimicrobial agents. *J. Pept. Sci.* **17**, 298–305 (2011).

44. Sohlenkamp, C. & Geiger, O. Bacterial membrane lipids: diversity in structures and pathways. *FEMS Microbiology Reviews* **40**, 133–159 (2016).

45. Angelova, M. & Dimitrov, D. S. Liposome electroformation. *Faraday Discuss. Chem. Soc.* **81**, 303–311 (1986).

46. Angelova, M. & Dimitrov, D.S. A mechanism of liposome electroformation. Trends in Colloid and Interface Science II, ed. Degiorgio, V. (Springer, Berlin), 59–67 (1988).

47. Angelova, M. Giant vesicles. Perspectives in Supramolecular Chemistry, eds Luisi, P. L. & Walde, P. (Wiley-Interscience, Chichester, UK), 1st Ed, 27–36 (2000).

48. Lindsey, H., Petersen, N. O. & Chan, S. I. Physicochemical characterization of 1,2-diphytanoyl-sn-glycero-3-phosphocholine in model membrane systems. *Biochim. Biophys. Acta* **555**, 147–167 (1979).

49. Redwood, W. R., Pfeiffer, F. R., Weisbach, J. A. & Thompson, T. E. Physical properties of bilayer membranes formed from a synthetic saturated phospholipid in n-decane. *Biochim. Biophys. Acta* **233**, 1–6 (1971).

50. Tosatto, L. et al. Alpha-synuclein pore forming activity upon membrane association. *Biochim. Biophys. Acta* **1818**, 2876–2883 (2012).

51. Allen, F. H. & Johnson, O. Automated conformational analysis from crystallographic data. 4. Statistical descriptors for a distribution of torsion angles. *Acta Cryst.* **B47**, 62–67 (1991).

52. Sochacki, K. A., Barns, K. J., Bucki, R. & Weisshaar, J. C. Real-time attack on single *Escherichia coli* cells by the human antimicrobial peptide LL-37. *Proc. Natl. Acad. Sci. USA* **108**, E77–E81 (2011).

53. Choi, H., Yang, Z. & Weisshaar, J. C. Oxidative stress induced in *E. coli* by the human antimicrobial peptide LL-37. *PLoS Pathog.* **13**, e1006481 (2017).

54. Hsiao, Y.-W. et al. Cooperative modes of action of antimicrobial peptides characterized with atomistic simulations: a study on cecropin B. *J. Phys. Chem. B* **122**, 5908–5921 (2018).

55. Gautier, R., Douguet, D., Antonny, B. & Drin, G. HELIQUEST: a web server to screen sequences with specific α-helical properties. *Bioinformatics.* **24**, 2101–2102 (2008).

## Acknowledgements

NMR experiments described in this paper were produced using the facilities of the Centre for Biomolecular Spectroscopy, King's College London, acquired with a Multi-user Equipment Grant from the Wellcome Trust and an Infrastructure Grant from the British Heart Foundation. CDL acknowledges the stimulating research environment provided by the EPSRC Centre for Doctoral Training in Cross-Disciplinary Approaches to Non-Equilibrium Systems (CANES, EP/L015854/1). PMF is supported by a Health Schools Studentship funded by the EPSRC (EP/M50788X/1).

## Author Contributions

G.M., D.A.P. and A.J.M. designed the study. G.M., P.M.F., and A.J.M. wrote the main manuscript text and prepared all figures. P.M.F. and C.D.L. performed and/or analysed atomistic simulation data. S.S.B. provided materials. G.M., H.A. and R.A.A. performed structural studies. G.M., V.B.G., H.A., T.T.B. and A.F.D. performed and/or analysed CD measurements. G.M. and V.B.G. performed and analysed patch-clamp measurements. G.M., C.H., M.C. and M.J.S. conducted antimicrobial susceptibility testing. All authors approved the manuscript.

## Additional Information

Supplementary information accompanies this paper at https://doi.org/10.1038/s41598-019-47327-w.

**Competing Interests:** The authors declare no competing interests.

**Publisher's note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_8.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International

License, which permits use, sharing, adaptation, distribution and reproduction in any medium or
format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Cre-
ative Commons license, and indicate if changes were made. The images or other third party material in this
article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the
material. If material is not included in the article's Creative Commons license and your intended use is not per-
mitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the
copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2019