Conformation and membrane interaction studies of the potent antimicrobial and anticancer peptide palustrin-Ca

Palustrin-Ca (GFLDIIKDTGKEFAVKILNNLKCKLAGGCPP) is a host defence peptide with potent antimicrobial and anticancer activities, first isolated from the skin of the American bullfrog Lithobates catesbeianus. The peptide is 31 amino acid residues long, cationic and amphipathic. Two-dimensional NMR spectroscopy was employed to characterise its three-dimensional structure in a 50/50% water/2,2,2-trifluoroethanol-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$d_{3}$$\end{document}d3 mixture. The structure is defined by an \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha$$\end{document}α-helix that spans between Ile\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^{6}$$\end{document}6-Ala\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^{26}$$\end{document}26, and a cyclic disulfide-bridged domain at the C-terminal end of the peptide sequence, between residues 23 and 29. A molecular dynamics simulation was employed to model the peptide’s interactions with sodium dodecyl sulfate micelles, a widely used bacterial membrane-mimicking environment. Throughout the simulation, the peptide was found to maintain its \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha$$\end{document}α-helical conformation between residues Ile\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^{6}$$\end{document}6-Ala\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^{26}$$\end{document}26, while adopting a position parallel to the surface to micelle, which is energetically-favourable due to many hydrophobic and electrostatic contacts with the micelle.

www.nature.com/scientificreports/ While HDPs are typically positively charged, anionic molecules such as lipopolysaccharides (LPS), phospholipids and teichoic acids are prevalent in prokaryotic cell membranes, which enables attractive electrostatic interactions between HDPs and their target bacterial cell membranes 19,20 . While peptide cationicity is generally a requirement for attraction to membranes, the loss of N-terminal positive charge has been shown to allow a deeper peptide insertion into the lipid bilayer 21 , although a separate study has found that reduced N-terminal cationicity results in reduced antimicrobial activity 22 . Similarly to prokaryotic membranes, the cell membranes of cancer cells are enriched in negative charge, unlike non-cancerous cell membranes which are neutral, which allows anticancer peptides to be selectively electrostatically attracted to cancer cells.
HDPs usually exert their bactericidal activity via disruption of the bacterial cell membrane, for which three main mechanisms of action have been described: the toroidal pore model, barrel-stave model and carpet model. All three membranolytic mechanisms of action are relatively non-specific, which is the reason why few cases of bacterial resistance to HDPs have been described 23,24 . The initial step in all three mechanisms is the HDP binding of the target membrane, initially facilitated via the aforementioned electrostatic attraction, then via hydrophobic attraction. As the peptide concentration at the membrane increases, eventually a threshold concentration is reached, whereupon the peptides begin to induce changes in the membrane structure 25 . As the peptide concentration continues to increase, a second threshold concentration is reached, at which membrane lysis occurs 25 . This concentration may be determined by the organism's capacity to repair its damaged cell membranes 26 .
Generally, membrane lysis occurs either via the pore-forming toroidal pore or barrel-stave models, or the non-pore-forming carpet model. The latter requires that the bilayer curvature is disrupted, thereby disintegrating the bacterial membrane. Under the toroidal pore model, the bacterial membrane is bent from the outer leaflet inwards, resulting in the formation of a transmembrane toroidal pore lined by the peptides and the membrane's lipid head groups 27,28 . The barrel-stave model, meanwhile, involves the peptides forming a transmembrane pore; its lumen is lined by the peptides' hydrophobic residues, and its exterior is composed of hydrophilic residues, which are in contact with the membrane's lipid headgroups.
Where an amphipathic peptide structure is modified so that a residue side chain on its hydrophobic face has switched polarity, membrane lysis can also occur via a fractal rupture mechanism. This mechanism relies on apertures in the upper leaflet of the lipid bilayer being caused by chemically induced peptide-lipid interfaces, thereby exposing the membrane's hydrophobic layer and resulting in an energetically unfavourable water-lipid interface, and ultimately in structured defects in the membrane 29 .
Interestingly, some HDPs exert their activity without disruption of the lipid bilayer. Instead, these peptides translocate across the membrane without causing lysis, and instead exert antimicrobial activity by targetting an intracellular target 30 .
Our laboratory has conducted structural studies on various amphibian HDPs, including ranatuerin-2CSa 31 , XT-7 32 , alyteserin-1c 33 , brevinin-1BYa 34 and its analogues 35 , maximin 3 36 and maximin 1 37 . Palustrin-Ca (GFL-DIIKDTGKEFAVKILNNLKCKLAGGCPP) is a 31 amino acid residue HDP first isolated from the skin of the American bullfrog Lithobates catesbeianus 11 . In common with other HDPs isolated from frogs belonging to the Ranidae family, it contains a cyclic disulfide bridged domain at the C-terminal end of the peptide sequence, between residues 23 and 29. Palustrin-Ca is biologically interesting, as it is non-haemolytic, and exhibits potent broad-spectrum antibacterial activity, and as previously mentioned, very potent anticancer activity, with an IC 50 of 0.951 μg/ml against human gastric cancer SGC-7901 11 . In this work, palustrin-Ca's conformation and structural properties in a 50% H 2 O-trifluoroethanol-d 3 mixture are elucidated using NMR spectroscopy, with an ensemble of model structures being determined. Additionally, the peptide's interactions with a membrane-mimetic are simulated using an atomistic molecular dynamics simulation.

Materials and methods
Materials and NMR sample preparation. The palustrin-Ca peptide (MW = 3304 g/mol, purity > 95%) was purchased from ProteoGenix (Paris). The identity and purity of the peptide were confirmed using mass spectrometry and high performance liquid chromatography (Supplementary Information). 4.5 g of the peptide was dissolved in 0.6 mL of a 50% (v/v) TFE-d 3 /H 2 O solution, resulting in a final peptide concentration of 2.27 mM. 3-trimethylsilyl propionic acid (TSP) and 2,2,2-trifluoroethanol (TFE-d 3 ) of analytical grade were obtained from Sigma-Aldrich (Ireland). NMR spectroscopy. A Bruker Avance 600 NMR spectrometer with a 5 mm inverse probe head at a 1 H resonance frequency of 600.13 MHz was used to perform a number of NMR experiments at a temperature of 298 K. One-dimensional proton, 2D phase-sensitive total correlation spectroscopy (TOCSY) 38 , nuclear Overhauser effect spectroscopy (NOESY) 39 and natural abundance 1 H-13 C and 1 H-15 N heteronuclear single quantum coherence spectroscopy ( 1 H-13 C-HSQC) and 1 H-15 N-HSQC 40 spectra were acquired, with relaxation delays of 2.5 s, 2.0 s, 1.5 s, 1.0 s and 2.0 s, respectively, and acquisition times of 3.4 s, 340 ms, 280 ms, 170 ms and 100 ms, respectively. Mixing times of 60 ms and 200 ms were used for the TOCSY and NOESY, respectively.
The 1 H spectral widths were 6.0 kHz, 7.2 kHz, 6.0 kHz and 9.6 kHz, and the spectra were acquired with 4, 16, 32 and 64 transients for each of the 1024, 2048, 256 and 128 t1 increments for the TOCSY, NOESY, 1 H-13 C-HSQC and 1 H-15 N-HSQC, respectively. The 13 C spectral width was 21.1 kHz and the 15 N spectral width was 10.9 kHz. All two-dimensional spectra were processed with the Bruker TopSpin program, version 4.0.6 (Bruker BioSpin, Germany), using the sine squared window function, and the 1 H signal of TSP was used as the chemical shift reference. Structure calculation. The NMRFAM-SPARKY program, version 3.131 41 was used to analyse the acquired NMR spectra and integrate the NOESY cross-peaks. The integrated peak volumes were exported, and the nuclear www.nature.com/scientificreports/ Overhauser effect (NOE) cross-peak intensities were calibrated using the CALIBA 42 program, yielding a set of upper distance restraints. Protons that could not be stereospecifically assigned were treated as pseudoatoms. The set of obtained distance restraints was used as input to CYANA 43 , excluding those that represent fixed distances. A force constant of 1 kJ mol −1 Å −2 was used to weight the distance restraints, from which one hundred structures were generated by CYANA 44 , and submitted to 20,000 steps of simulated annealing and 20,000 of conjugate gradient minimization. The 20 structures with the lowest target function values were selected, and subjected to an additional 2000 steps of conjugate gradient energy minimization with restrained backbone atoms, using the CHARMM22 force field 45 in NAMD, version 2.12 46 . VMD (visual molecular dynamics), version 1.9.3 47 , was used to analyse the final ensemble of 20 structures. The geometry, stereochemical quality and structural statistics of the final ensemble of structures was assessed and validated using PROCHECK 48 and the wwPDB web service 49 .

Molecular dynamics simulation.
Structural model coordinates of an SDS micelle were obtained from Jakobtorweihen et al. 50 , and used in the construction of an SDS micelle-peptide system, which was constructed by aligning the centres of mass of the micelle and the medoid palustrin-Ca structure with VMD. To ensure that the system's charge remained neutral, chloride ions were added as counter-ions, and finally, TIP3P water was used to solvate the system. The system was energy minimized and equilibrated using the CHARMM22 all-atom forcefield 45,51 within NAMD version 2.12 46 . All calculations were conducted within the NPT ensemble, using the Langevin piston Nose-Hoover method 52,53 and periodic boundary conditions. Long-range non-bonded interactions were calculated up to a switching distance of 8.5 Å, beyond which a smooth switching function truncated the energy to a cut-off of 11 Å. The PME method 54 was used to calculate long-range electrostatic interactions at each time step. The non-bonded interaction list was also updated every step. All bonds to hydrogen atoms were constrained using the SHAKE algorithm 55 . A 2 fs timestep was employed.
The solvated peptide-micelle system was minimized for 2000 conjugate gradient steps with peptide backbone atoms fixed, and a further 2000 steps with its α-carbon atoms restrained. The system was then heated to a temperature of 310 K over 6000 steps using Langevin dynamics, with α-carbon remaining restrained. The system volume was equilibrated with the Langevin piston at 1 atm over 24,000 steps, with a further 24,000 steps without restraints. Finally, the system was simulated for 57.4 ns.
The radial distribution function g(r) was used for analysing the final 22.5 ns of simulation data. The measure gofr command within VMD was used to calculate the g(r) function between each residue's carbon atoms and the SDS micelle's non-sulfate atoms. The default values 0.1 and 25.0 Å were used for δr and max r, respectively.

Results and discussion
Conformational analysis by NMR. Two-dimensional NMR spectroscopy was used to determine the structure of palustrin-Ca in a 50:50 TFE-d 3 /water solvent mixture, which is a standard solvent mixture employed in the characterisation of peptide structures that promotes the formation of secondary structures, including α -helices 56 . Structural studies of another ranid frog-derived HDP, brevinin-1BYa, and its analogues in TFE and SDS micelles have shown that the peptide conformations in both media are comparable 34,35 . The effects that TFE exerts on peptide structures are well-studied, having been the subject of research for approximately 60 years 57,58 .
Peptides unrestrained by disulfide bonds typically are unstructured in pure aqueous solution [59][60][61][62] , as α-helical self-aggregation is hindered by the presence of cationic or sterically bulky residues. Furthermore, sequence enrichment in lysine, which is common among AMPs, can preclude the formation of β-sheets. Consequently, most linear cationic HDPs do not possess a defined secondary structure in pure aqueous solution 63,64 .
Nonetheless, the addition of TFE has been demonstrated to induce the formation of secondary structure by enhancing α-helical character in peptides 65 . This has been shown in studies that systematically increase TFE concentration 34,35,59,66 .
It is generally understood that the mechanism by which peptide α-helical character is enhanced commences with the preferential aggregation of TFE molecules around the peptide, facilitated by hydrophobic interaction between peptide's hydrophobic sidechains and the trifluoromethyl group, which allows the TFE molecules to substitute the peptide's hydration shell. This replacement removes alternative hydrogen-bonding partners, and results in a low dielectric environment, comparable to that experienced by peptides in biological membranes, which favours the formation of intramolecular hydrogen bonds between the backbone amide groups 67,68 .
NMR conformational studies of magainin 2 demonstrated that in a pure aqueous solution, the peptide possesses an extended random coil structure 69 . Similarly, the maximin 4 antimicrobial peptide does not have a well-defined structure in pure aqueous solution 70 . Previously, a 50% TFE concentration has been found to be appropriate for the characterisation of peptide structures, as low TFE concentrations result in poor amide chemical shift dispersion, which indicates a lack of defined secondary structure 71 .
The collected spectra were well-resolved, with well-dispersed peaks that facilitated unambiguous resonance assignment. Individual residue spin systems were identified with the aid of TOCSY, 1 H-13 C HSQC and 1 H-15 N HSQC spectra, while sequence-specific resonance assignment was conducted using the HN-HN and HN-Hα regions of the NOESY spectrum (Fig. 1). The first sequence-specific resonance assignments were made for the unique residues Thr 9 , Glu 12 and Val 15 , with the remainder being assigned through the use of a backbone-walk, and investigation of HN-Hα region peaks, where necessary. The Gly 1 amide proton chemical shift could not be determined due to chemical exchange. The presence of a disulfide bridge was confirmed by observation of clear long-range β i -β i+6 NOE peaks between residues Cys 23 -Cys 29 .
All the 1 H chemical shifts identified for palustrin-Ca are detailed in Table 1. Likewise, all of palustrin-Ca's intermolecular NOE connectivities are summarised in Fig. 2   Similarly, α-proton chemical shift deviations were calculated as the difference between the observed chemical shift and the random coil chemical shift ( �δ = δ obs − δ rc ). Four sequential α-proton chemical shift deviations at least 0.1 ppm lower than the expected random coil values indicate α-helical segments, while three or more deviations at least 0.1 ppm higher are indicative of β-strands; no change greater than 0.1 ppm is indicative of coiled regions 73 . Figure 4 shows the chemical shift deviations of palustrin-Ca, which shows a clear stretch of downfield shifts between residues Phe 2 -Asp 8 , and again between residues Lys 11 -Leu 25 . While these chemical shift deviations are indicative of α-helical structure, the former stretch of downfield shifts may instead be attributed to a turn conformation. While the α-proton chemical shift deviations for residues Ala 26 -Pro 30 are generally upfield, this trend is not sufficiently consistent to attribute any secondary structure to it. Molecular modelling. The NOESY cross-peaks were integrated, and the resultant volumes were converted into distance restraints. CYANA was employed to generate one hundred structures, of which the twenty possessing the lowest target function values were further energy-minimized. Table 2 summarizes the structural and energetic statistics of the twenty selected models. As depicted in Fig. 5, palustrin-Ca is predominantly α -helical in TFE-d 3 , which is known to promote the formation of the secondary structure. The ensemble's second- www.nature.com/scientificreports/   www.nature.com/scientificreports/ ary structural properties were analysed with the aid of STRIDE 74 . Most of the ensemble's structures are display α-helical structures between residues Ile 6 -Ala 26 . The terminal segments are predominantly defined as turns, but possess significant coil character as well. More specifically, the N-terminus exhibits mostly turn character, although the coil character of the C-terminus is greater. The peptide structure is amphipathic, with a clear sepa- Positively charged residues are coloured in cyan, negatively charged residues are coloured in pink, polar residues are coloured in purple, glycine residues are coloured in orange, cysteine residues are coloured in green, hydrophobic residues are coloured in slate. Table 2. Mean structural statistics of the twenty structural models comprising the palustrin-Ca ensemble. The peptide structure has been deposited in the PDB with deposition code 7P4X. Model 10 is the representative medoid structure. a RMSD values from VMD. b Based on PROCHECK.

Intraresidue 147
Interresidue, sequential 144 Interresidue, nonsequential 130 Ensemble RMSD values, Å, SA ± SD a www.nature.com/scientificreports/ ration between the polar and non-polar faces of the helix. This is unsurprising, considering that amphipathicity is a common requirement for antimicrobial activity 75 .

Molecular dynamics.
A molecular dynamics simulation of palustrin-Ca's medoid structure in an SDS micelle was conducted to characterise its behaviour in the prokaryotic membrane-mimetic environment. SDS micelles are often used as mimetics of the bacterial membrane in molecular dynamics simulations [76][77][78][79][80] . In the course of the simulation, the peptide translocated from its initial position aligned with the micelle's centre of mass, to the micelle's surface-water boundary, where it adopted a position parallel to the micelle surface (Fig. 6). The peptide's amphipathicity is clearly apparent; its hydrophobic residues remain closely associated with the micellar interior, while the hydrophilic residues are in contact with the aqueous solvent and the micelle's anionic sulfate headgroups. This amphipathic configuration is quite energetically favourable, and is contributes to the peptide's exertion of its antimicrobial activity 81 . The peptide's cationic residues are important for its bioactivity, as the electrostatic interaction between the target surface's anionic charges and a peptide's cationic residues drives the initial attraction of the peptide to the target 82 . Likewise, the hydrophobic residues are necessary for the peptide's insertion into the target hydrophobic core 83 . The peptide's interactions with the micelle's hydrophobic core were quantified using the radial distribution function g(r), which was calculated between the peptide's carbon atoms and the micelle's aliphatic chain, and plotted against the radius r, as shown in Fig. 7.
The radial distribution function plots illustrate clearly that the peptide's hydrophobic residues maintain a close association with the micelle, especially in the structured α-helical segment, where the hydrophobic Ile 6 , Phe 13 , Ala 14 , Ile 17 , Leu 18 , Leu 21 and Leu 25 residues are particularly closely associated to the micelle hydrophobic core. Conversely, hydrophilic residues, including Asp 8 , Glu 12 , Lys 16 , Asn 19 and Lys 22 are less associated with the micelle core, and in close contact with the aqueous solvent.
Interestingly, the N-terminal segment of the peptide is also found to maintain close association with the micelle, especially through the Gly 1 , Phe 2 and Ile 6 residues, despite the first five residues not forming part of a defined secondary structure. This is in contrast to the C-terminal segment, where beginning with Ala 26 onwards the peptide is only weakly associated with the micelle.
The peptide, which was experimentally determined to be α-helical between residues Ile 6 and Ala 26 , maintains this structure throughout the duration of the simulation. Despite the presence of a glycine residue at position 10, the α-helix remains relatively rigid, not exhibiting any significant flexibility during the simulation. This suggests that the glycine residue occurs too far from the helical segment's middle to impart flexibility. Although studies have previously demonstrated that a rigid α-helical structure can potentially result in an increased haemolytic activity 84 , this does not appear to apply in this case, which may be explained by the relatively flexible nonhelical N-terminal and C-terminal segments, which combined account for 10 of the sequences 31 residues. This Figure 6. Palustrin-Ca positioned at the SDS micelle-water boundary. Lysine residue side chains are shown interacting with the negatively charged sulfate headgroups. Positively charged residues are coloured in cyan, negatively charged residues are coloured in pink, polar residues are coloured in purple, glycine residues are coloured in orange, cysteine residues are coloured in green, hydrophobic residues are coloured in slate. SDS aliphatic chains are coloured in white, while the SDS headgroups are orange and red, and highlighted with dotted spheres. www.nature.com/scientificreports/ flexibility can be seen in the root-mean-square fluctuation (RMSF) plot, where the peptide's terminal residues have greater RMSF values than the residues located within the α-helical segment (Fig. 8).

Conclusion
A number of biologically active peptides have been isolated from the American bullfrog Lithobates catesbeianus. The isolated peptides are diverse, and not limited to a single peptide family; palustrin-Ca is the only palustrin peptide isolated from this amphibian. Palustrin-Ca can be considered to be a member of the palustrin-2 family www.nature.com/scientificreports/ of peptides, sharing an average 42.09% sequence identity with 18 peptides of the palustrin-2 family. Peptides that belong to the palustrin-2 family exhibit a variety of potent biological activities, including antimicrobial and anticancer activities, and therefore represent exciting candidates for novel drug development. Palustrin-Ca is 31 amino acid residues long, of which 5 are lysine, 4 are glycine and 4 are leucine. The peptide has a cyclic disulfide-bridged heptapeptide domain at its C-terminus, which is conserved with other peptides of the palustrin-2 family.
This study employed NMR spectroscopy and molecular modelling methods to characterise the peptide's threedimensional structure and obtain an ensemble of model structures. To date, palustrin-Ca is the only palustrin peptide to have had its three-dimensional structure elucidated. The results established that palustrin-Ca possesses an amphipathic, α-helical structure between residues Ile 6 -Ala 26 in a 50% TFE-H 2 O mixed solvent system, with the terminal segments predominantly defined as turns, with significant coil character as well.
A molecular dynamics simulation was conducted, whereby the medoid NMR-determined structure was simulated with an SDS micelle. The simulation results suggest that palustrin-Ca preferentially adopts a position parallel to the micelle's surface. This configuration is most energetically favourable as the peptide's hydrophobic residues can penetrate to the micelle's hydrophobic core, while the hydrophilic residues remain in contact with the aqueous solvent. The RDF plot in Fig. 7 illustrates this clearly; the greatest values for g(r) at small distances are observed for the hydrophobic Leu, Ile, Ala and Phe residues, demonstrating the importance of hydrophobicity for a peptide's ability to bind the target membrane. The observed preference for a position parallel to the micelle surface indicates that the peptide most probably exerts its antimicrobial activity through a non-pore-forming mechanism of action, such as the carpet model or the interfacial activity model. This study's results will facilitate future work focused on investigating HDP structures, and their interactions with zwitterionic lipid bilayers, in order to further our understanding of the relationships between HDP structure and function. This study has extended earlier work from our research group that led to the elucidation of a number of host-defence peptide structures, and also complements related bioinformatic studies on peptide structure and function 85-88 . Received: 3 August 2021; Accepted: 3 November 2021 Figure 8. Root-mean-square fluctuation (RMSF) calculated for the backbone atoms of each residue of palustrin-Ca over the course of the simulation. It is apparent that the peptide structure is generally rigid, with the exception of the N-and C-terminal residues, which are relatively flexible.