Design and Characterization of Myristoylated and Non-Myristoylated Peptides Effective against Candida spp. Clinical Isolates

The increasing resistance of fungi to antibiotics is a severe challenge in public health, and newly effective drugs are required. Promising potential medications are lipopeptides, linear antimicrobial peptides (AMPs) conjugated to a lipid tail, usually at the N-terminus. In this paper, we investigated the in vitro and in vivo antifungal activity of three short myristoylated and non-myristoylated peptides derived from a mutant of the AMP Chionodracine. We determined their interaction with anionic and zwitterionic membrane-mimicking vesicles and their structure during this interaction. We then investigated their cytotoxic and hemolytic activity against mammalian cells. Lipidated peptides showed a broad spectrum of activity against a relevant panel of pathogen fungi belonging to Candida spp., including the multidrug-resistant C. auris. The antifungal activity was also observed vs. biofilms of C. albicans, C. tropicalis, and C. auris. Finally, a pilot efficacy study was conducted on the in vivo model consisting of Galleria mellonella larvae. Treatment with the most-promising myristoylated peptide was effective in counteracting the infection from C. auris and C. albicans and the death of the larvae. Therefore, this myristoylated peptide is a potential candidate to develop antifungal agents against human fungal pathogens.


Introduction
In recent years, extensive abuse of chemicals to control diseases has been associated with the increasing number of multidrug-resistant pathogens. The growth of fungi and bacteria resistant to conventional antimicrobial agents represents one of the biggest concerns worldwide [1][2][3][4].
Systemic fungal infections are challenging and represent a significant cause of morbidity and mortality, especially for immune-compromised people. Fungi infect billions of people every year, with a high rate of mortality of about 1.5 million [1]. Pathogenic fungi, like human cells, are eukaryotic, and potential selective drugs can only attack limited targets. Many antifungal drugs act on the biosynthetic pathway of ergosterol, the

Peptide Design and Physico-Chemical Properties
The first step in the activity and selectivity of cationic AMPs is the electrostatic interaction with the negatively charged bacterial membrane [23][24][25]. Thus, different strategies are used to improve the antibacterial activity with the aim to create more efficient potential agents in clinical. The most-used approach is to find the correct balance between the net positive charge and the amphipathic properties to obtain peptides that can fold in α-helices or β-structures [22,[26][27][28]. In the present study, we designed small peptides derived from the sequence of KHS-Cnd [22]. KHS-Cnd was chosen because of its high antimicrobial activity against Gram-positive, Gram-negative, and multidrug-resistant bacteria and relatively low cytotoxicity. We synthesized three small peptides of 11 amino acids from the parental peptide. The sequences and physicochemical properties of the three derivatives corresponding to the N-terminal, central, and C-terminal portions of KHS-Cnd, named Pep-A, Pep-B, and Pep-C, are reported in Table 1. Many studies report that the functionalization with fatty acid chains increases antifungal activity [20,21]. Thus, we modified the N-terminal by introducing a hydrophobic myristoyl moiety (14-carbon fatty acid) to favor the vehiculation and interaction of peptides with membranes and enhance their membrane-anchoring properties [29,30]. This modification increased the hydrophobicity of peptides without varying the peptide sequence. Physico-chemical properties of myristoylated peptides, called Myr-A, Myr-B, and Myr-C, are also reported in Table 1.

Critical Micellar Concentration (CMC) Determination
Lipopeptides usually self-assemble in organized structures that have attracted interest because of their enhanced bactericidal and fungicidal effects [31,32]. The self-assembly of myristoylated peptides was evaluated by determining the critical micellar concentration (CMC) that represents the concentration of spontaneous associations of monomers by a spectrofluorometric method. 8-anilinonaphtalene-1-sulfonate (ANS) was used as a fluorescent probe since its fluorescence is very sensitive to the microenvironment, showing a blue shift and intensity increase upon changing from a polar to non-polar medium [33]. Therefore, the fluorescence of ANS was monitored as a function of the lipopeptide concentration. Figure 1A shows the emission spectra of ANS at increasing concentrations of myristoylated peptides (Myr-C). Figure 1B shows the plot of fluorescence intensity of ANS (the maximum emission wavelength was 465 nm) vs. the logarithm of the lipidated peptide concentration. Fluorescence intensity slowly increases up to a certain lipopeptide concentration and then breaks abruptly, indicating a sudden change in the aggregation state of lipopeptides. For each peptide, the critical micellar concentration (CMC) was evaluated by the intersection of the two lines. The values of CMC are 4.7 ± 0.8 µM, 5.1 ± 0.8 µM, and 5.1 ± 0.4 µM for Myr-A, Myr-B, and Myr-C, respectively (data for Myr-A and Myr-B are reported in Supplementary Information Figures S1 and S2).  Non-lipidated peptides did not show a break in the fluorescence intensity upon titration, suggesting that they do not aggregate.

Membrane Partition Studies
To quantify the interaction of different lipidated and non-lipidated peptides with membrane mimetic systems, we measured the intrinsic fluorescence of Trp-1 upon partition with LUVs (large unilamellar vesicles) of different compositions. Briefly, the mixtures of POPC/POPG (70%/30% w/w) LUVs were used to mimic anionic membranes and POPC only to mimic zwitterionic membranes [34][35][36]. LUVs do not entirely describe the complexity of natural membranes; still, they are used here to evaluate the preference of different peptides for anionic (bacteria and fungi) and zwitterionic (mammalian cell) membranes. The mole fraction partition constants Kx were calculated from the partition isotherms as previously described [37]. Figure 2 reports the partition isotherms for the titration of myristoylated and non-myristoylated peptides with POPC and POPC/POPG vesicles.
The values of Kx for all peptides are reported in Table 2 together with G values, calculated using the expression G = −RT lnKx, and the selectivity ratio defined as the ratio between partition constants determined in anionic and zwitterionic membranes. Higher values of the selectivity ratio represent a preference of the peptide for anionic vesicles, while lower values indicate a preference for the zwitterionic one.
Partitioning data suggest that myristoylated peptides interact with LUVs with higher affinity. The mole fraction partition coefficients Kx range from 3.0 × 10 4 to 5.2 × 10 4 with calculated Gibbs free energies of the partition of ~ 26 kJ/mol for non-myristoylated peptides, suggesting poor binding to membranes. Myristoylated peptides showed larger Kx values ranging from 1.0 × 10 5 in POPC to 1.7 × 10 6 POPC/POPG vesicles, indicating a stronger interaction with LUVs. Non-lipidated peptides did not show a break in the fluorescence intensity upon titration, suggesting that they do not aggregate.

Membrane Partition Studies
To quantify the interaction of different lipidated and non-lipidated peptides with membrane mimetic systems, we measured the intrinsic fluorescence of Trp-1 upon partition with LUVs (large unilamellar vesicles) of different compositions. Briefly, the mixtures of POPC/POPG (70%/30% w/w) LUVs were used to mimic anionic membranes and POPC only to mimic zwitterionic membranes [34][35][36]. LUVs do not entirely describe the complexity of natural membranes; still, they are used here to evaluate the preference of different peptides for anionic (bacteria and fungi) and zwitterionic (mammalian cell) membranes. The mole fraction partition constants K x were calculated from the partition isotherms as previously described [37]. Figure 2 reports the partition isotherms for the titration of myristoylated and non-myristoylated peptides with POPC and POPC/POPG vesicles.    The values of K x for all peptides are reported in Table 2 together with ∆G values, calculated using the expression ∆G = −RT lnK x , and the selectivity ratio defined as the ratio between partition constants determined in anionic and zwitterionic membranes. Higher values of the selectivity ratio represent a preference of the peptide for anionic vesicles, while lower values indicate a preference for the zwitterionic one. Partitioning data suggest that myristoylated peptides interact with LUVs with higher affinity. The mole fraction partition coefficients K x range from 3.0 × 10 4 to 5.2 × 10 4 with calculated Gibbs free energies of the partition of~(−26) kJ/mol for non-myristoylated peptides, suggesting poor binding to membranes. Myristoylated peptides showed larger K x values ranging from 1.0 × 10 5 in POPC to 1.7 × 10 6 POPC/POPG vesicles, indicating a stronger interaction with LUVs.
Moreover, myristoylated peptides prefer to interact with anionic membranes as described by the values of the selectivity ratio ranging between~2 and 4. This feature should increase the antifungal activity, as described later. Mc Laughlin and colleagues stated that the presence of the myristoyl moiety is required to interact with membranes but is not always sufficient [38]. The higher propensity of myristoylated peptides to interact with membranes can be explained by the presence of the myristoyl group and charge effects. Mainly, the values of Gibbs free energy~(−33) kJ/mol are similar to that measured for the insertion of myristate into the bilayer [39]. Thus, the synergism between the myristate insertion into the bilayer and the electrostatic interaction between the lysine residues present in the peptides and the anionic lipids can be the driving force of the process [40].

Fluorescence Quenching Experiments
Peptides' topology in the absence and presence of lipid vesicles has been evaluated by monitoring the quenching of Trp-1 by acrylamide. Acrylamide is a small neutral molecule that quenches the tryptophan residue with high efficiency [41]. Quenching experiments were carried out in buffer and lipid vesicles as previously described [42]. Typical results are shown in Figure 3 for myristoylated and non-myristoylated peptides.

Fluorescence Quenching Experiments
Peptides' topology in the absence and presence of lipid vesicles has been evaluated by monitoring the quenching of Trp-1 by acrylamide. Acrylamide is a small neutral molecule that quenches the tryptophan residue with high efficiency [41]. Quenching experiments were carried out in buffer and lipid vesicles as previously described [42]. Typical results are shown in Figure 3 for myristoylated and non-myristoylated peptides. In the buffer, all non-myristoylated peptides and myristoylated peptides Myr-A and Myr-C showed typical linear behavior of the fluorescence intensity ratio in the absence (Fo) and in the presence (F) of the quencher vs. the quencher concentration. The same behavior was detected for myristoylated peptides in POPC and POPC/POPG lipid vesicles. The classical Stern-Volmer equation [43] 1 (1) was used to describe the fluorescence quenching, where KSV is the Stern-Volmer or collisional constant, while Fo and F are the fluorescence in the absence and presence of the acrylamide Q, respectively. Values of KSV are indicative of the dynamic quenching and are higher in buffer solution than in the presence of lipid vesicles ( Table 3), suggesting that Trp-1 is less accessible to the quencher in the presence of lipid vesicles and confirming the partition into the lipid bilayer. As reported in Figure 3, Myr-B in the buffer and non-myristoylated peptides in the presence of lipid vesicles showed non-linear behavior, with downward curvature at high quencher concentrations. This effect can be ascribed to selective quenching [44] due to the presence of two different fluorophore populations (a and b) with a diverse solvent exposition: one more accessible to the quencher and the other less. The fluorescence intensity can be written as: The fraction of total fluorescence has been determined using the following equation [43]: In the buffer, all non-myristoylated peptides and myristoylated peptides Myr-A and Myr-C showed typical linear behavior of the fluorescence intensity ratio in the absence (F o ) and in the presence (F) of the quencher vs. the quencher concentration. The same behavior was detected for myristoylated peptides in POPC and POPC/POPG lipid vesicles. The classical Stern-Volmer equation [43] was used to describe the fluorescence quenching, where K SV is the Stern-Volmer or collisional constant, while F o and F are the fluorescence in the absence and presence of the acrylamide Q, respectively. Values of K SV are indicative of the dynamic quenching and are higher in buffer solution than in the presence of lipid vesicles ( Table 3), suggesting that Trp-1 is less accessible to the quencher in the presence of lipid vesicles and confirming the partition into the lipid bilayer. As reported in Figure 3, Myr-B in the buffer and non-myristoylated peptides in the presence of lipid vesicles showed non-linear behavior, with downward curvature at high quencher concentrations. This effect can be ascribed to selective quenching [44] due to the presence of two different fluorophore populations (a and b) with a diverse solvent exposition: one more accessible to the quencher and the other less. The fluorescence intensity can be written as: The fraction of total fluorescence has been determined using the following equation [43]: 7 of 23 Table 3. Stern-Volmer constant (K sv and K a ), fraction accessible to the quencher (f a ), and Net Accessibility Factor (NAF) for myristoylated peptides in the presence/absence of vesicles of different compositions. NAF and fa are defined in the text. where f a is the fraction of the accessible fluorophore at an infinite concentration of the quencher and K a is an average Stern-Volmer constant for the accessible fluorophore both exposed and buried. A plot of, F o F−F o vs. 1 [Q] yields a straight line from which we determined the fitting parameters f a and K a ( Figures S3 and S4). Values of f a for non-myristoylated peptides in the presence of lipid vesicles are between 0.40 and 0.75, indicating two fluorophore populations, one accessible to the quencher and the other deeply buried in the phospholipid bilayer.

Peptide LUV Composition
Furthermore, we calculated the NAF, defined as the ratio between the Stern-Volmer constant in the presence and absence of LUV [45]: Low values of NAF suggest protection of the fluorophore from the quencher and strong interaction with lipids. NAF values reported in Table 3 are lower for myristoylated peptides than for non-myristoylated, suggesting a better interaction with the lipid vesicles. Peptide Myr-B has a lower value of NAF (0.07), indicating a strong interaction with the lipid and a preference for anionic lipid vesicles. Non-myristoylated peptides show high values of NAF, suggesting a weak interaction with lipid vesicles. All the data agree with the results obtained from the partition studies.

Circular Dichroism Experiments
Structural transitions of peptides upon the addition of lipid vesicles have been investigated by CD spectroscopy [42]. We explored structural changes of both myristoylated and non-myristoylated peptides upon interaction with anionic and zwitterionic membranemimicking systems. In Figure 4, typical titration curves for myristoylated peptides are reported at increasing lipid vesicle concentrations. The explored peptide-to-lipid ratios ranged from 1:0 to 1:40. CD data confirm that all peptides can interact with lipid vesicles of different compositions. Peptides Myr-A and Myr-C show similar CD spectra in buffer (black curves) with a minimum at~220 nm, typical of a β-conformation. This conformation seems to be stable, and it also persists upon the addition of POPC lipid vesicles. However, there is a slight transition towards a partially helical conformation in the presence of anionic lipid vesicles. The peptide Myr-B shows a typical spectrum of a random coil in buffer and in the presence of lipid vesicles of different compositions. are reported at increasing lipid vesicle concentrations. The explored peptide-to-lipid ratios ranged from 1:0 to 1:40. CD data confirm that all peptides can interact with lipid vesicles of different compositions. Peptides Myr-A and Myr-C show similar CD spectra in buffer (black curves) with a minimum at 220 nm, typical of a β-conformation. This conformation seems to be stable, and it also persists upon the addition of POPC lipid vesicles. However, there is a slight transition towards a partially helical conformation in the presence of anionic lipid vesicles. The peptide Myr-B shows a typical spectrum of a random coil in buffer and in the presence of lipid vesicles of different compositions. The fractional helical content fα for the peptides Myr-A and Myr-C in the presence of POPC/POPG lipid vesicles has been estimated using the formula [46]: (5) where is the observed ellipticity at 222 nm and and are the values for a completely random coil and an α-helical conformation, respectively. According to the Luo and Baldwin [47]   The fractional helical content f α for the peptides Myr-A and Myr-C in the presence of POPC/POPG lipid vesicles has been estimated using the formula [46]: where θ is the observed ellipticity at 222 nm and θ RC and θ H are the values for a completely random coil and an α-helical conformation, respectively. According to the Luo and Baldwin [47] formula, the limiting values are θ RC = −360 and θ H = −39.375. The helical fractions f α in the presence of POPC/POPG lipid vesicles are~0.2 and~0.1 at an L/P molar ratio of 40 for Myr-A and Myr-C, respectively.

Molecular Dynamics Simulations
To investigate more deeply, with an atomic resolution, the interaction between the peptides and the membrane models, we carried out a series of molecular dynamics simulations of peptides in different environments, including water, TFE/water (50% v/v), and in the presence of POPC and POPC/POPG lipid bilayers.
The initial structure of peptides used for the MD in the solution was obtained using the I-Tasser bioinformatics tool, which suggested, for all peptides, an α-helix conformation.
The MD simulations in water of all peptides showed a rapid loss of the initial α-helix conformation ( Figure S5). Only for Myr-B was the formation of a transient α-helix observed during 300 ns of simulation. Analysis of the peptide's secondary structure as a function of the simulation time suggests that all peptides are unstructured with high conformational flexibility in the water, in agreement with the experimental results obtained by CD spectra.
In the presence of TFE/water (50% v/v), a less polar solvent with respect to water, Pep-B, Pep-C, and Myr-C conserved the α-helix conformation for 300 ns with a helicity of about 45%, 64%, and 82%, respectively ( Figure S5). In the case of Pep-A and Myr-A, the analysis of the secondary structure showed an unfolding of peptides after 120 ns and 220 ns of simulation ( Figure S5A,B).
We also simulated the adsorption of peptides on the surface of anionic and zwitterionic membrane models. The structures obtained at the end of the MD simulations in TFE/water were used as the initial conformation of peptides for the simulations in the presence of POPC and POPC/POPG lipid bilayers. The peptides were positioned close to the surface of the membranes with the principal axis parallel to the lipid surface and a distance of 1.0 nm between the center of mass of the peptide and the phosphorous atoms of lipids. Figure 5 shows the representative snapshots of MD simulation at 1 µs of myristoylated peptides in the presence of POPC and POPC/POPG lipid bilayers. Myr-A and Myr-B interact with the surface of the zwitterionic and anionic membrane through an unfolded conformation, while Myr-C shows a higher tendency with respect to the other myristoylated peptide to fold in the presence of lipid bilayers. In particular, analysis of the time evolution of the secondary structure of Myr-C in POPC/POPG shows that the peptide conserves helicity (α-helix and π-helix) of about 45% (five residues) during the simulation duration ( Figure S6). On the other hand, in the presence of POPC, Myr-C interacts with the membrane, adopting, in the last 500 ns of simulation, a 3 10 -helix with helicity of about 45.5% ( Figure S6).
The adsorption of myristoylated peptides to the surface of the lipid bilayer is driven mainly by the alkyl chain, linked to the N-terminal residue of peptides, and by the presence of positively charged amino acids. The alkyl chain penetrates in the hydrophobic region of the lipid bilayer while the positively charged amino acids stabilize the peptide-membrane interaction by forming hydrogen bonds and electrostatic interactions with lipid molecules. The peptide Myr-A in POPC shows a deep penetration of the myristoyl chain and Trp-1 in the hydrophobic region of the membrane ( Figures 5A and S7A). In particular, the average position of Trp-1 is located slightly below the position of the carbonyl group (CO sn-1 and CO sn-2 ) of POPC, as evidenced by the analysis of the average density distribution calculated to the lipid bilayer center ( Figure S7A). In POPC/POPG ( Figures 5B and S7B), the residues of Trp-1 and Lys-4 are located close to the lipid/water surface. The residue Arg-7 interacts with phosphate groups forming an arginine-phosphate salt bridge. The different penetrations into the lipid bilayer of POPC and POPC/POPG observed by MD simulation for Trp-1 of the peptide Myr-A is in agreement with the NAF values obtained by fluorescence-quenching experiments. In POPC/POPG, Trp-1 is close to the lipid surface and then more exposed to the quencher molecules with respect to POPC, where the aromatic residue can penetrate the hydrophobic region of the membrane. The adsorption of myristoylated peptides to the surface of the lipid bilayer is driven mainly by the alkyl chain, linked to the N-terminal residue of peptides, and by the presence of positively charged amino acids. The alkyl chain penetrates in the hydrophobic region of the lipid bilayer while the positively charged amino acids stabilize the peptidemembrane interaction by forming hydrogen bonds and electrostatic interactions with lipid molecules. The peptide Myr-A in POPC shows a deep penetration of the myristoyl chain and Trp-1 in the hydrophobic region of the membrane (Figures 5A and S7A). In particular, the average position of Trp-1 is located slightly below the position of the carbonyl group (COsn-1 and COsn-2) of POPC, as evidenced by the analysis of the average density distribution calculated to the lipid bilayer center ( Figure S7A). In POPC/POPG (Figures 5B and S7B), the residues of Trp-1 and Lys-4 are located close to the lipid/water surface. The residue Arg-7 interacts with phosphate groups forming an arginine-phosphate salt bridge. The different penetrations into the lipid bilayer of POPC and POPC/POPG observed by MD simulation for Trp-1 of the peptide Myr-A is in agreement with the NAF values obtained by fluorescence-quenching experiments. In POPC/POPG, Trp-1 is close  (Figure 5C,D). Myr-C, which adopts a folding conformation (5-6 residues) when adsorbed by the membrane, in the presence of POPC shows residues of Trp-1 and all lysine residues except Lys-11 located inside the membrane close to the carbonyl group of lipids ( Figures 5E and S7E). On the contrary, in POPC/POPG, all charged amino acids of the peptide Myr-C are located close to the phosphate group (Figures 5F and S5F) and interact with lipid molecules via the salt bridge and hydrogen bonds. In the case of Myr-B and Myr-C, analysis of the density distribution calculated with respect to the lipid bilayer center suggests that in POPC/POPG, the residue of Trp-1 is located more inside the membrane with respect to POPC, in agreement with the factor NAF obtained by fluorescence experiments.
The analysis of MD trajectories of non-myristoylated peptides in POPC and POPC/POPG reveals a lower affinity of peptides for lipid bilayers with respect to myristoylated peptides, as shown by the partition constants K x reported in Table 2. In particular, the peptides Pep-B and Pep-C are adsorbed on the lipid surface and cannot deeply penetrate the membrane but lie on the surface, with the non-charged amino acids oriented toward the water phase ( Figures S8 and S9). Figure 6 shows, for example, the different penetrations of Pep-C and Myr-C adsorbed on the lipid surface of POPC/POPG. POPC shows residues of Trp-1 and all lysine residues except Lys-11 located inside the membrane close to the carbonyl group of lipids ( Figures 5E and S7E). On the contrary, in POPC/POPG, all charged amino acids of the peptide Myr-C are located close to the phosphate group (Figures 5F and S5F) and interact with lipid molecules via the salt bridge and hydrogen bonds. In the case of Myr-B and Myr-C, analysis of the density distribution calculated with respect to the lipid bilayer center suggests that in POPC/POPG, the residue of Trp-1 is located more inside the membrane with respect to POPC, in agreement with the factor NAF obtained by fluorescence experiments.
The analysis of MD trajectories of non-myristoylated peptides in POPC and POPC/POPG reveals a lower affinity of peptides for lipid bilayers with respect to myristoylated peptides, as shown by the partition constants Kx reported in Table 2. In particular, the peptides Pep-B and Pep-C are adsorbed on the lipid surface and cannot deeply penetrate the membrane but lie on the surface, with the non-charged amino acids oriented toward the water phase ( Figures S8 and S9). Figure 6 shows, for example, the different penetrations of Pep-C and Myr-C adsorbed on the lipid surface of POPC/POPG. The density distribution of Trp-1 atoms, calculated with respect to the lipid bilayer center ( Figure S8), shows that for all non-myristoylated peptides, Trp-1 does not penetrate the hydrophobic portion of the membrane but remains located close to phosphate groups of lipids at the lipid/water interface. At the lipid/water interface, the residue Trp-1 is more exposed to water molecules and all molecular species present in the water phase. MD simulation results suggest that Trp-1 of myristoylated peptides is deeply immersed in the The density distribution of Trp-1 atoms, calculated with respect to the lipid bilayer center ( Figure S8), shows that for all non-myristoylated peptides, Trp-1 does not penetrate the hydrophobic portion of the membrane but remains located close to phosphate groups of lipids at the lipid/water interface. At the lipid/water interface, the residue Trp-1 is more exposed to water molecules and all molecular species present in the water phase. MD simulation results suggest that Trp-1 of myristoylated peptides is deeply immersed in the hydrophobic region of the membrane and, therefore, poorly exposed to water, with respect to non-myristoylated peptides. These results are in good agreement with the fluorescence experiments.
The MD simulations evidenced that the presence of the myristoyl moiety linked to Pep-A, Pep-B, and Pep-C deeply influences the properties of peptides, particularly their structure and affinity for lipid bilayers.

Antifungal Activity
The MIC 90 values of non-myristoylated and myristoylated peptides against different Candida spp. are summarized in Table 4. All lipopeptides have much lower MIC 90 values than their unfunctionalized counterparts, with the sole exception of the Myr-A peptide that maintains an unchanged MIC 90 value towards C. auris. As already observed in the literature, these results confirm that increasing hydrophobicity by fatty acylation of AMPs with myristic acid enhances their antimicrobial activity and, in particular, the antifungal properties [48].  (Table S1). MIC values are higher than CMC, so MIC values are determined in the presence of lipopetide micellar aggregates and not the free monomer. Under these conditions, lipopetides are still soluble even if in the form of micellar aggregates. Thus, the activity of lipopetides might also depend on the so-called micellar mechanism, consisting of the solubilization of part of the host membrane into mixed micelles. However, further studies are needed to verify this hypothesis.

SEM Analysis
The antibiofilm activity of the Myr-B lipopeptide was tested on representative highbiofilm-producing strains of C. albicans, C. tropicalis, and C. auris. These strains were selected for their ability to form biofilms in vitro as assessed by the crystal violet assay (data not shown). SEM images clearly show that Myr-B is able to demolish the biofilm of all three fungal isolates at a concentration of 128 µg/mL (85 µM) (Figure 7). Untreated biofilms show a well-organized three-dimensional structure, with the presence of an extracellular matrix for the C. albicans and C. tropicalis strains and fungal hyphae, well-branched and adhered to the surface. The biofilm of C. auris, as already reported [49], is characterized by cellular aggregates with a poor extracellular matrix able to tolerate clinical concentrations of sodium hypochlorite. After treatment with the Myr-B lipopeptide, the biofilm of all fungal strains is strongly demolished. Few cellular hyphae fragments can be observed, but the morphological organization is completely lost. Biofilm demolition results validate the idea the peptide's myristoylation strongly enhances the antifungal properties also towards the biofilm form. These results have a significant impact on the search for new peptide drugs able to combat the rising antimicrobial resistance and also effective against biofilm-related infections.

In Vivo Toxicity and Efficacy Testing in Galleria Mellonella Model
Multidrug-resistant C. auris was chosen as a representative species for in vivo studies. After determining the infective fungal dose of 10 7 CFU/injection was capable of killing about 50% of the Galleria mellonella larvae 24 h after infection (data not shown), we evaluated the in vivo toxicity and efficacy of the Pep-B and Myr-B peptides (the Myr-B peptide was selected as the best active peptide due to the previously shown results). Groups of ten G. mellonella larvae were infected with C. auris and treated with the Pep-B and Myr-B peptides at a concentration of 640 µg/mL (about 400 µM). The results shown in Table 5 highlighted that the toxicity of Pep-B is slightly higher than that of the corresponding lipopeptide. Furthermore, Myr-B was able to contain the effects of fungal infection 24 h after the treatment, with a single addition of Myr-B: 8 out of 10 larvae survived, whereas, in the control group, only 4 survived. In the group of infected larvae treated with Pep-B, 7 out of 10 larvae died after 24 h, showing poor efficacy compared to the control. In these pilot and preliminary studies with a single drug administration, the data obtained after 24 h are the most significant, as in the following 48 h, the injected compound undergoes inactivation and elimination likely due to the action of internal proteases. Therefore, we could also confirm the antifungal activity of the Myr-B peptide in vivo.   Note. * The control group larvae were infected with the pathogen but not treated with peptides.

Cytotoxicity Assay against Human Cell Line
The cytotoxic effect of the six peptides has been investigated on a primary human fibroblast cell line (FB789). Four different peptide concentrations (5,10,20, and 50 µM) and two different time points (12 and 24 h) have been tested (see Figure 8)

Hemolytic Assay against Rabbit Erythrocytes.
Rabbit erythrocytes have been used to investigate the hemolytic effect of the six peptides (Figure 9). Hemolysis values from 5 to 20 μM are rather low, ranging from 0.1% (Peptide-A and Peptide Myr-A at 5 μM, see panel A) to 7% (Peptide Myr-B, see panel B). An increase in the membrane lysis induced by the peptides is visible at 50 μM, with an

Hemolytic Assay against Rabbit Erythrocytes
Rabbit erythrocytes have been used to investigate the hemolytic effect of the six peptides ( Figure 9). Hemolysis values from 5 to 20 µM are rather low, ranging from 0.1% (Peptide-A and Peptide Myr-A at 5 µM, see panel A) to 7% (Peptide Myr-B, see panel B). An increase in the membrane lysis induced by the peptides is visible at 50 µM, with an effect of 28.5% of hemolysis reached by Peptide Myr-C (see panel C). At the concentrations of Myr-B active against Candida species, no or very few hemolytic effects have been evidenced.  3). The asterisks indicate the significance level with respect to negative control (0% of toxicity): *** = p < 0.001; ** = p < 0.005; * = p < 0.05.

Hemolytic Assay against Rabbit Erythrocytes.
Rabbit erythrocytes have been used to investigate the hemolytic effect of the six peptides ( Figure 9). Hemolysis values from 5 to 20 μM are rather low, ranging from 0.1% (Peptide-A and Peptide Myr-A at 5 μM, see panel A) to 7% (Peptide Myr-B, see panel B). An increase in the membrane lysis induced by the peptides is visible at 50 μM, with an effect of 28.5% of hemolysis reached by Peptide Myr-C (see panel C). At the concentrations of Myr-B active against Candida species, no or very few hemolytic effects have been evidenced.

Conclusions
In recent years, fungal infections have been increasing worldwide, and they are usually associated with relatively high mortality rates. Still, the current antimycotic therapies are limited. Moreover, some Candida species resistant to antibiotics have emerged and new antifungal drugs need to be developed [50]. In this paper, we reported the design and characterization of three lipopeptides active in vitro against different multidrug-resistant Candida species. Lipopeptides were selective and active against fungal cells and not harmful against either fibroblast mammalian cells or erythrocytes. Finally, as the first proof of evidence, we tested the in vivo activity of the most-active peptide, Myr-B, against an infection of Candida auris in insect larvae of Galleria mellonella. The peptide was not toxic and could contain the effects of fungal infection 24 h after the treatment. Therefore, it could be considered a new promising drug candidate. These results could greatly impact the search for new antifungal peptide drugs able to combat the rising antibiotic resistance and against biofilm-related infections.

Materials
All the reagents were purchased from Sigma-Merck (Milan, Italy) unless otherwise specified. All peptides and lipo-peptides were purchased from CASLO ApS, c/o Scion Technical University of Denmark, with a purity >98%. Peptide concentrations were determined by UV spectroscopy at 280 nm before each sample preparation.

Lipid Vesicles Preparation
The lipids POPC (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine) and POPG (1palmitoyl-2-oleoyl-sn-glycero-3-phosphoglycerol) were purchased from Avanti Polar Lipids (Alabaster, AL, USA) and used without further purification. LUVs (Large Unilamellar Vesicles) were prepared as previously reported [37]. Briefly, appropriate amounts of lipids were dissolved in chloroform/methanol 9:1. The solvent was then removed by rotary evaporation and samples were placed overnight under a high vacuum. The lipid film was then hydrated by adding 1 mL of buffer (20 mM phosphate buffer at pH 7.4 with 150 mM NaCl and 0.8 mM EDTA) and then subjected to 5 freeze-thaw cycles and vortexed.
The multilamellar vesicles (MLVS) were extruded through a polycarbonate membrane with pores of 100 nm by an Avanti polar mini extruder, and the obtained LUVs were used within 48 h of preparation. LUVs were composed of 100% POPC and a 70%/30% (w/w) combination of POPC/POPG. The final concentration of the lipid was 10 mM.

Fluorescence Spectroscopy
All the steady-state fluorescence measurements were carried out using a Perkin Elmer LS55 fluorescence spectrophotometer. All the experiments were conducted at 298 K in a thermostatic cell holder and in 20 mM phosphate buffer at pH 7.4 containing 0.8 mM EDTA and 150 mM NaCl unless otherwise specified.

Critical Micellar Concentration (CMC) Determination
The CMC of the N-terminal myristoylated peptides was determined in water by fluorescence spectroscopy using 8-anilinonaphtalene-1-sulfonate (ANS) as the fluorescent probe [51]. The ANS assays were performed using peptide concentrations of 0.8-10.0 µM in 70 µM ANS solution.
The ANS was excited at λ ex = 356 nm, and the fluorescence spectrum was measured for λ em = 400-670 nm with excitation and emission bandwidths of 5 nm.

Partition Constant Determination
Since all peptides contain tryptophan, fluorescence spectroscopy is an appropriate tool for studying these molecules. Peptide partitioning between buffer and lipid vesicles was monitored by measuring the increase in tryptophan fluorescence upon the addition of LUVs. Trp fluorescence spectra were recorded with λ exc = 295 nm and scanning the emission between 305 and 400 nm. Measurements were carried out with a cross-oriented configuration of polarizers (pol em = 0 • and pol exc = 90 • ) [52]. Briefly, a 1.0 µM solution of the peptide was titrated with LUVs of different compositions (POPC 100% and 70%/30% POPC/POPG) as previously described [22,37]. The background effects of vesicles and buffer were subtracted from each spectrum. All the fluorescence emissions were corrected for dilution. Mole fraction partition coefficients K x were calculated from the fraction of the partitioned peptide f p as previously described [22,37]. The fraction of the partitioned peptide f p is given by the following equation [53] where K x is the mole fraction partition constant, and [L] and [W] are the molar concentrations of lipids and water, respectively.

Acrylamide Quenching Experiments
Quenching of tryptophan fluorescence both in the absence and presence of lipid vesicles was accomplished by the addition of small aliquots of acrylamide dissolved in buffer. The concentration was 2.0 µM for both myristoylated and non-myristoylated peptides whereas lipid vesicles were 400 µM in the experiments with myristoylated peptides and 1.0 mM in the experiments with non-myristoylated peptides. The excitation wavelength was set at 295 nm, and fluorescence emission spectra were recorded between 305 and 400 nm. Measured fluorescence was corrected for dilution of the sample, and data were fitted according to the Stern-Volmer or modified Stern-Volmer equation [43].

Circular Dichroism Spectroscopy
Circular Dichroism (CD) spectra were recorded in the 190-260 nm spectral range at 298 K on a J 715 JASCO spectropolarimeter with 1.0 nm intervals, using a 0.2 cm path-length quartz cuvette. The reported spectra are the average of 8 scans with an instrument scanning speed of 100 nm/min, a step size of 1 nm, and 1s response. Liposomes for CD experiments were prepared as described in Section 4.2 using a 10 mM phosphate buffer (pH 7.4 and 0.1 mM EDTA).
CD spectra of peptides were recorded in 10 mM phosphate buffer (pH 7.4 and 0.1 mM EDTA) in the absence and presence of LUVs of different compositions. The general procedure for non-myristoylated peptides was as follows: To a freshly prepared solution of 30 µM peptide, increasing aliquots of the proper liposome solution were added. After each addition, the solution was left for 5 min to equilibrate under stirring before recording the CD spectrum. The explored lipid/peptide ratios ranged from 1:0 to 1:60.

Simulation Systems
The initial structure of Pep-A, Pep-B, and Pep-C were obtained using the fold recognition algorithm implemented in I-TASSER (iterative threading assembly refinement algorithm) [54]. Arginine and lysine side chains were considered positively charged according to their protonation state at pH = 7. The structure of myristoylated peptides, Myr-A Myr-B, and Myr-C, were obtained by linking, with an amide bond, the acyl chain to the N-terminal tryptophan residue of Pep-A, Pep-B, and Pep-C using Avogadro software [55].
Four systems were constructed for each peptide. In the first system, the peptide was immersed in a cubic box 5.0 × 5.0 × 5.0 nm 3 and solvated with at list 4000 water molecules. In the second system, the peptide was immersed in a cubic box 5.0 × 5.0 × 5.0 nm 3 and solvated with a mixture of TFE/water (50% v/v). In the third system, the peptide was positioned close to the surface of the pre-equilibrated bilayer of POPC with the principal axis parallel to the lipid surface and a distance of 1.0 nm between the center of mass of the peptide and the phosphate atoms of lipids. A rectangular box of 6.5 × 6.5 × 9.0 nm 3 was used with at least 6500 water molecules. In the fourth system, each peptide was positioned above the pre-equilibrated bilayer surface of POPC/POPG with the principal axis parallel to the lipid surface and a distance of 1.0 nm between the center of mass of the peptide and the phosphate atoms of lipids. A rectangular box of 6.2 × 6.2 × 9.0 nm 3 was used with at least 6500 water molecules.
For the simulation in water and the mixture of TFE/water, the best-scored structures of each peptide obtained by the I-TASSER prediction were used. For simulations in the presence of the lipid bilayer of POPC and POPC/POPG, the structure of the peptide obtained at the end of the MD simulation in TFE/water was used to prepare the initial configuration of the system with the peptide close to the surface of the bilayer. The initial configuration of POPC and POPC/POPG (70:30) was built with 128 lipids (64 per leaflet). The mixed bilayer of POPC/POPG (70:30) was built using 92 molecules of POPC and 36 molecules of POPG (18 per leaflet).
Van der Waals and electrostatic interactions were calculated using a cut-off of 1.2 nm, and the particle mesh Ewald method (PME) [61] was applied to the long-range electrostatic interactions beyond 1.2 nm. Covalent bonds were constrained using the P-LINCS algorithm [62,63], whereas the geometry of water was constrained with the SETTLE algorithm [64]. Newton's equation of motion was integrated via the leap-frog algorithm using a time step of 2 fs.
The velocity-rescale thermostat [65] was used to control the temperature of the simulated systems by coupling the peptide, lipids, water, and ions separately to an external bath at 298 K with a coupling constant τ T of 0.1 ps. The pressure was controlled using the Berendsen barostat [66] (P = 1 bar, τ p = 1.0 ps). The pressure coupling was isotropic in the case of the simulations in water and in the mixture of TFE/water (50% v/v) and semi-isotropic in the presence of lipid bilayers. Periodic boundary conditions were applied in all three dimensions.
The trajectories obtained by MD simulation were analyzed with the GROMACS analysis tools, VMD 1.9.3 [67], and in-house scripts. Analysis of the secondary structure was performed with the gmx do_dssp tool of GROMACS using the program DSSP [68]. The antifungal activity of AMPs was determined against 50 clinical isolates, 10 for each of the following Candida species: C. albicans, C. tropicalis, C. glabrata, C. parapsilosis, and C. auris. The latter were derived from patients with positive blood cultures in four hospitals in India, with known resistance profiles to fluconazole, voriconazole, and echinocandins. Yeast cells were grown in Brain Heart Infusion (BHI) medium (Sigma-Aldrich, Gillingham, UK) for 16 h at 37 • C in a 150-rpm orbital shaker. Cells were then sub-inoculated in a fresh BHI medium and grown to an optical density (OD) of 0.3. The turbidity of the inoculum was adjusted to 0.5 McFarland and diluted 1:500 in RPMI 1640 broth, corresponding to around 2.5 × 10 5 CFU/mL.

MIC Evaluation
Minimal inhibitory concentrations (MICs) were determined by Broth microdilution susceptibility tests, according to The Clinical and Laboratory Standards Institute (CLSI) international guidelines [69]. A fungal suspension of 0.5 MacFarland was diluted 1:500 in RPMI 1640 broth and incubated in 96 multi-well plates containing AMPs at different concentrations ranging from 0.125 to 256 µg/mL. In this range of concentration, lipopeptides were completely soluble. MIC was defined as the lowest concentration of the peptides at which yeasts showed no visible growth compared to the control. The MIC 90 represents the value at which ≥90% of the strains within a test population were inhibited.

Biofilm Conditions and Biofilm Demolition Assays
Fungal strains of C. albicans, C. tropicalis, and C. auris were cultured in BHI and grown at 37 • C overnight, harvested by centrifugation (3000× g, 10 min), resuspended in sterile PBS, and adjusted to a density of 10 9 colony-forming units (CFU)/mL in RPMI media by measuring the absorbance at 600 nm. A 24-well plate with 13-mm-diameter glass coverslips (Thermanox NUNC, New York, NY, USA) was seeded with 500 mL and grown 3 days at 37 • C for biofilm formation. Once the biofilms were formed, after 3 days incubation, the culture medium was removed, and cells were washed three times with PBS to remove non-adherent cells. The fungal biofilms, highly visible on the disk surface, were then treated with a 200 µL final volume of peptide X in PBS at two different concentrations (128 and 16 µg/mL) for 16 h at 37 • C. The demolition activity was observed in scanning electron microscopy (SEM).

Galleria Mellonella Model
Ten larvae (250-320 mg each) were selected at random for each step in the procedure. Any larvae with darkening of the cuticle were discarded. The test compounds and yeast were injected into the hemocoel in water buffer through the last left proleg (Insumed 30G Pic, Italy, syringe volume 0.5 mL, needle size 30G) [70]. The larvae were incubated in the dark for 3 days, and mortality was recorded daily.
An infective dose of yeast Candida auris was determined by injecting groups of ten larvae with a suspension at 10 5 , 10 6 , 10 7 , and 10 8 colony-forming units (CFU) per injection. The larvae were incubated for 3 days at 30 • C. An infective dose was defined as the one that caused an immune response, recognizable by the darkening of the cuticle and able to kill approximately 50% of larvae after 24 h [70].
For the peptide efficacy experiment, ten larvae were injected into the last left proleg with a pre-determined infective dose of yeast and the dose of peptides at the same time in a final volume of 10 µL and incubated at 30 • C. Each step in the procedure included six groups of 10 larvae: C. auris-infected larvae (10 7 CFU), infected and treated larvae with peptide B (640 µg/mL), infected and treated larvae with peptide Myr-B (640 µg/mL), treated larvae with peptide B (640 µg/mL), treated larvae with peptide Myr-B (640 µg/mL), and vehicle-injected control (water). The mortality was recorded daily for 3 days.

Cytotoxicity Assay against Human Cell Line
The cytotoxicity of the six peptides was tested on a primary human fibroblast cell line (FB789) grown in Dulbecco's Modified Eagle Medium (DMEM) containing 10% fetal calf serum (FCS) and antibiotics (penicillin-streptomycin, Gibco) at 37 • C in a humidified 5% CO 2 atmosphere. Cells were seeded on 96-well microplates at a density of 3 × 10 3 cells per well in 100 µL of the medium. Four dilutions of the six peptides (from 5 µM to 50 µM) were added to the cells and maintained for 12 and 24 h. As a negative control, cells were grown in a normal medium without the peptides, while, as a positive control, cells were added to NaN 3 10% v/v. The cytotoxicity was determined by measuring the intracellular adenosine triphosphate (ATP) levels using the luciferase-based ATPlite assay (PerkinElmer, Waltham, MA, USA), according to the manufacturer's instructions. After 12 and 24 h, the cells were lysed, and the lysates were transferred into opaque well plates (OptiPlate-96, PerkinElmer). The amount of emitted light, linearly correlated with the ATP concentration, was measured with a microplate luminometer (Victor II PerkinElmer) for 10 min in the dark. Three replicates for each dilution were performed. Cell viability values were expressed as the mean + SD and calculated as the percent values of the treated samples with respect to the untreated cells (negative control). The statistical analysis was performed using one-way ANOVA analysis followed by the Dunnett post-test.

Hemolytic Activity Assay
The hemolytic assay was performed on rabbit erythrocytes (Rockland) maintained in Alsever's solution (Innovative Research). Briefly, after the removal of Alsever's solution, the erythrocytes were resuspended in PBS and counted in a hemocytometer. A suspension of 5,000,000 red blood cells was incubated with four dilutions (from 5 µm to 50 µm) of the six peptides in a 96-well microplate. As a negative control, erythrocytes were incubated without the peptides, while as a positive control, erythrocyte was added with triton 10% v/v. The plate was incubated at 37 • C for 2 h and subsequently centrifuged at 1200 rpm × 3 min to separate the pellet from the supernatant. The absorbance of the supernatant was measured at 492 nm. Each point was made in triplicate. The relative OD compared to the positive control defined the percentage of hemolysis [71]. Statistical analysis was performed using two-way ANOVA analysis followed by Bonferroni's post-test. Funding: This research was partially funded by Regione Lazio, grant number A0375-2020-36557, in the call "Gruppi di Ricerca 2020". Name of the project: "New Antimicrobial Peptides delivered in specific nanosystems against ESKAPE bacteria: a weapon against antibiotic resistance (NanopeptoBat)".

Informed Consent Statement: Not Applicable.
Data Availability Statement: The data are available on request from the corresponding author.

Conflicts of Interest:
The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results. Net accessibility factor POPC 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine POPG 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoglycerol SEM Scanning electron microscope TFE trifluoroethanol