Rational Design of Novel Peptidomimetics against Influenza A Virus: Biological and Computational Studies

Effective therapy against the influenza virus is still an unmet goal. Drugs with antiviral effects exist, but the appearance of resistant viruses pushes towards the discovery of drugs with different mechanisms of action. New anti-influenza molecules should target a good candidate, as a new anti-influenza molecule could be an inhibitor of the influenza A virus hemagglutinin (HA), which plays a key role during the early phases of infection. In previous work, we identified two tetrapeptide sequences, SLDC (1) and SKHS (2), derived from bovine lactoferrin (bLf) C-lobe fragment 418–429, which were able to bind HA and inhibit cell infection at picomolar concentration. Considering the above, the aim of this study was to synthesize a new library of peptidomimetics active against the influenza virus. In order to test their ability to bind HA, we carried out a preliminary screening using biophysical assays such as surface plasmon resonance (SPR) and orthogonal immobilization-free microscale thermophoresis (MST). Biological and computational studies on the most interesting compounds were carried out. The methods applied allowed for the identification of a N-methyl peptide, S(N-Me)LDC, which, through high affinity binding of influenza virus hemagglutinin, was able to inhibit virus-induced hemagglutination and cell infection at picomolar concentration. This small sequence, with high activity, represents a good starting point for the design of new peptidomimetics and small molecules.


Introduction
Influenza A virus (IAV) is an extremely infectious pathogen endowed with zoonotic potential and represents one of the leading causes of morbidity and mortality worldwide [1,2].Over the decades, it has triggered several pandemics, causing major public health problems as well as severe economic losses [3,4].Flu therapy is based on vaccine prevention and antiviral drug treatment.However, vaccines have several limitations.They are not always protective against the circulating strains, and they cannot be administered to immunosuppressed patients or, in any case, to those with a weak immune system.
Antivirals are an alternative approach to combating the various subtypes of influenza viruses [5].Antiviral compounds approved by the Food and Drug Administration (FDA) for flu treatment can be classified into three groups: the neuroaminidase (NA) inhibitors oseltamivir phosphate, zanamivir, peramivir, and laninamivir octanoate hydrate; the M2 ion channel blockers amantadine and rimantadine; and the cap-dependent endonuclease protein inhibitor baloxavir marboxil [6][7][8][9][10].However, the emerging threat of new pandemic influenza strains spreading through the human population and the rise of resistance to licensed drugs have focused research into developing new treatments for influenza viruses [11][12][13][14].A viral component whose role in virus infection is crucial should be the ideal target for novel anti-influenza treatments.
The IAV envelope contains two glycoproteins: hemagglutinin (HA) and neuraminidase (NA).Both play pivotal roles in the infection cycle: HA for binding to specific cell receptors and mediating membrane fusion, and NA for releasing viral particles from the cell surface, but they are also involved in viral attachment and entry.The majority of current treatment approaches for IAV target NA.Oseltamivir, the only oral NA inhibitor, is the most common anti-flu therapy worldwide, so its broad use led to the selection of resistant viral strains [15].Frequently, when a tyrosine residue replaces the histidine residue 274 (position 275 of the N1 protein) on the NA gene, influenza A/H1N1 viruses develop oseltamivir resistance [16].For example, during the 2007-2008 influenza season in Europe, many seasonal A(H1N1) viruses emerged with the H275Y NA mutation, conferring a high resistance to oseltamivir [17].It is important to note that NA mutations responsible for oseltamivir resistance can emerge independently of the use of this drug [17].In fact, a high proportion of oseltamivir-resistant viruses have been isolated from patients who have not received treatment, demonstrating a high rate of transmission [6].Moreover, as these mutants are known to maintain fitness, they can easily spread globally, representing a serious public health risk.Furthermore, it is important to note that potential pandemic viruses, such as avian influenza A(H5N1) and A(H7N9), which cause severe disease in humans, are prone to the development of oseltamivir resistance in treated patients [18][19][20].Finally, it has been demonstrated that mallards can transmit an avian influenza virus that is oseltamivir-resistant to chickens.As this mutant is able to cross the species barrier and also maintain fitness, it could be able to infect humans [21].Finding alternative antiviral drugs with distinct modes of action is therefore urgently needed to fight infection by emerging viruses with pandemic potential, especially those carrying drug resistance markers.However, attempts to target HA are on the rise.The discovery of broadly neutralizing antibodies that target the HA stem has led to increased attention for HA as a target for broadly cross-protective antivirals.These antibodies are under investigation for direct use as therapeutic agents and have inspired the design of smaller cyclic peptides based on the antibody-binding interface.
Previously, Superti and coworkers demonstrated that bLf, in particular the bLf C-lobe, can bind and inhibit HA [22].Successively, through a truncation library, we identified the tetrapeptides SLDC (1) and SKHS (2), derived from bLf C-lobe fragment 418-429, which were able to bind HA and inhibit cell infection in a concentration range of picomolar [23].
It is well known that the direct application of peptides as medicinal entities has some severe limitations, including high degradation by proteolytic enzymes and poor cell membrane permeability.We designed new peptidomimetics derived from tetrapeptides 1 and 2 to increase metabolic stability and improve bioavailability, receptor affinity, and selectivity to overcome these drawbacks.Biological and computational assays on synthesized compounds were carried out.Direct binding assays, a bioactivity profile, and computational studies led to the identification of a very potent N-methyl peptide.

Design
Novel peptidomimetics were designed starting from the tetrapeptide sequences SLDC (1) and SKHS (2).Two different strategies were used to reduce the conformational freedom of the peptide backbone and increase the affinity towards the target: N-methylation and the synthesis of peptoids.

Design of N-Methyl Peptides
In peptide chemistry, N-methylation has been demonstrated to be a useful technique in structure-activity relationship studies.In physiologically active peptides, substitution of natural amino acids with N-methyl amino acids has enabled analogs with better pharmaco-logical properties, such as enzymatic stability [24,25], receptor selectivity [26,27], enhanced potency [28][29][30], and bioavailability [31][32][33][34].In fact, N-methylation of the peptide backbone gives peptides high affinity for their targets, proteolytic stability, membrane permeability, and conformational rigidity.
Hence, starting with tetrapeptides 1 and 2, we synthesized the corresponding Nmethyl peptides (3-10, Table 1).Therefore, we exchanged in a systematic manner Ser, Leu, Asp, Cys, Lys, and His with their N α -methylated analogs.Repeated combinations of poly-N-substituted glycine (NSG) form peptidomimetic compounds known as peptoids [35].They are a class of artificial, unnatural peptide mimics that have a modular structure that allows a wide variety of structural components to be easily incorporated.Peptoids have the functionality of their side chains linked to the nitrogen of the polyamide backbone instead of the α-carbon, creating an achiral and flexible oligomeric backbone without hydrogen bond donors.This N-alkyl backbone repeat motif allows peptoids greater stability against proteolytic degradation than analogous peptides [36].Additionally, NSGs may be more hydrophobic and possess higher cellular permeability [37][38][39].Figure 1 shows the similarities between peptides and peptoids in the spacing of the side chains and carbonyl groups and the variations in the chirality of the two monomers [40].Hence, we synthesized the peptoid analogues of 1 and 2 (11)(12)(13)(14)(15)(16)(17)(18), replacing in a systematic manner Ser, Leu, Asp, Cys, Lys, and His with their NGS.

Peptoids
In the solid-phase synthesis of peptoid, there are two distinct methods to introduce a N-alkylglycine (peptoid residue) on the growing peptide chain: (i) the N-substituted glycine derivative, appropriately protected at the tertiary nitrogen atom, can be prepared distinctly and used as a building block in the solid-phase procedure (monomer method) [43,44], or, alternatively, (ii) the peptoid residue is created during the growing peptide chain by combining two submonomers, a haloacetic acid and a primary amine (submonomer approach) [35].
We chose and improved the submonomer method, which allows direct assembly of the functionalized peptoid residue from commercially available chemicals, to speed up the synthesis of peptoids.

Chemistry 2.2.1. N-Methyl Peptides
Over the decades, several methods for solid-phase synthesis of N-methylated peptides have been described [41,42].The major limitation is the coupling of the N-methyl amino acid to the previous amino acid, which generally occurs in low yield and requires expensive coupling reagents and double coupling due to the sterically hindered N-methyl site.In this work, we planned to use N-methylated amino acids directly as building blocks, aided by microwave (MW) irradiation.Fmoc-Rink Amide Chem matrix resin was used as solid phase support (Scheme 1).Hence, we synthesized the peptoid analogues of 1 and 2 (11)(12)(13)(14)(15)(16)(17)(18), replacing in a systematic manner Ser, Leu, Asp, Cys, Lys, and His with their NGS.

N-Methyl Peptides
Over the decades, several methods for solid-phase synthesis of N-methylated peptides have been described [41,42].The major limitation is the coupling of the N-methyl amino acid to the previous amino acid, which generally occurs in low yield and requires expensive coupling reagents and double coupling due to the sterically hindered N-methyl site.In this work, we planned to use N-methylated amino acids directly as building blocks, aided by microwave (MW) irradiation.Fmoc-Rink Amide Chem matrix resin was used as solid phase support (Scheme 1).

Peptoids
In the solid-phase synthesis of peptoid, there are two distinct methods to introduce a N-alkylglycine (peptoid residue) on the growing peptide chain: (i) the N-substituted glycine derivative, appropriately protected at the tertiary nitrogen atom, can be prepared distinctly and used as a building block in the solid-phase procedure (monomer method) [43,44], or, alternatively, (ii) the peptoid residue is created during the growing peptide chain by combining two submonomers, a haloacetic acid and a primary amine (submonomer approach) [35].
We chose and improved the submonomer method, which allows direct assembly of the functionalized peptoid residue from commercially available chemicals, to speed up the synthesis of peptoids.
The peptoid residue is created during the peptide chain elongation for the synthesis of compounds 11-18 by mixing two submonomers, a haloacetic acid, and a primary amine (submonomer technique).Initially, we used a 25% piperidine solution in DMF for 30 min to deprotect the Rink Amide resin.

Peptoids
In the solid-phase synthesis of peptoid, there are two distinct methods to introduce a N-alkylglycine (peptoid residue) on the growing peptide chain: (i) the N-substituted glycine derivative, appropriately protected at the tertiary nitrogen atom, can be prepared distinctly and used as a building block in the solid-phase procedure (monomer method) [43,44], or, alternatively, (ii) the peptoid residue is created during the growing peptide chain by combining two submonomers, a haloacetic acid and a primary amine (submonomer approach) [35].
We chose and improved the submonomer method, which allows direct assembly of the functionalized peptoid residue from commercially available chemicals, to speed up the synthesis of peptoids.
The peptoid residue is created during the peptide chain elongation for the synthesis of compounds 11-18 by mixing two submonomers, a haloacetic acid, and a primary amine (submonomer technique).Initially, we used a 25% piperidine solution in DMF for 30 min to deprotect the Rink Amide resin.
Bromoacetic acid was then coupled to the NH 2 -peptide resin in the presence of N, N -diisopropylcarbodiimide (DIC), and the halogen was subsequently removed using an excessive amount of primary amine.Peptide chain elongation was carried out according to the standard protocol by adding the appropriate amine and bromoacetic acid.The last Nalkylglycine residue was acetylated as described in Section 3. Peptoids were simultaneously cleaved from the resin and deprotected.The synthesis is shown in Scheme 2.
Bromoacetic acid was then coupled to the NH2-peptide resin in the presence of N, N′diisopropylcarbodiimide (DIC), and the halogen was subsequently removed using an excessive amount of primary amine.Peptide chain elongation was carried out according to the standard protocol by adding the appropriate amine and bromoacetic acid.The last N-alkylglycine residue was acetylated as described in Section 3. Peptoids were simultaneously cleaved from the resin and deprotected.The synthesis is shown in Scheme 2. The crude compounds were purified by preparative RP-HPLC.Analytical HPLC indicated a purity greater than 98%, and molecular weights were confirmed by ESI-MS.

Microscale Thermophoresis (MST)
The 16 compounds were subjected to MST screening, as described in Materials and Methods (Section 3), and the results are shown in Table 1.According to the analysis, every peptide interacts with HA with different dissociation constants.Compounds 4, 9, 16, and 18 bind HA more effectively than compounds 1 and 2, showing an equilibrium dissociation constant (KD) value in a two-digit nanomolar concentration range (Figure 2).The crude compounds were purified by preparative RP-HPLC.Analytical HPLC indicated a purity greater than 98%, and molecular weights were confirmed by ESI-MS.

Direct Binding Assays 2.3.1. Microscale Thermophoresis (MST)
The 16 compounds were subjected to MST screening, as described in Materials and Methods (Section 3), and the results are shown in Table 1.According to the analysis, every peptide interacts with HA with different dissociation constants.Compounds 4, 9, 16, and 18 bind HA more effectively than compounds 1 and 2, showing an equilibrium dissociation constant (KD) value in a two-digit nanomolar concentration range (Figure 2).

Surface Plasmon Resonance (SPR)
Full-length recombinant HA (His-Tag) was immobilized up to about 11,000 response units (RU) on a sensor chip for the SPR study (Materials and Methods, Section 3).On the sensor surface, compound binding caused a change in refractive index [45].A regeneration step was required (glycine, pH 1.5).After injection, a running buffer was allowed to flow over the surface, and the dissociation of compounds from the surface was detected.On the other hand, no HA was immobilized on the control flow cell, so no appreciable signal alteration occurred.

Surface Plasmon Resonance (SPR)
Full-length recombinant HA (His-Tag) was immobilized up to about 11,000 response units (RU) on a sensor chip for the SPR study (Materials and Methods, Section 3).On the sensor surface, compound binding caused a change in refractive index [45].A regeneration step was required (glycine, pH 1.5).After injection, a running buffer was allowed to flow over the surface, and the dissociation of compounds from the surface was detected.On the other hand, no HA was immobilized on the control flow cell, so no appreciable signal alteration occurred.
The SPR study revealed that synthesized peptides effectively interact with the immobilized protein.The sensorgrams of the compounds 4, 9, 16, and 18 bound to HA in HBS-P buffer are shown in Figure 3. Interestingly, these peptidomimetics bind HA with higher efficiency with respect to 1 and 2, showing a K D value in a three-digit nanomolar concentration range.
Specific binding between compounds 4, 9, 16, and 18 and HA was shown by both direct binding assays.

Antiviral Activity 2.4.1. Toxicity of Compounds in MDCK cells
Two-fold serial dilutions of selected compounds (starting at 25 µM) were incubated for 72 h at 37 • C with semi-confluent MDCK cells grown in 96-well tissue culture microplates.Under these conditions, no peptide affected cell viability up to the highest concentration.
The SPR study revealed that synthesized peptides effectively interact with the immobilized protein.The sensorgrams of the compounds 4, 9, 16, and 18 bound to HA in HBS-P buffer are shown in Figure 3. Interestingly, these peptidomimetics bind HA with higher efficiency with respect to 1 and 2, showing a KD value in a three-digit nanomolar concentration range.Specific binding between compounds 4, 9, 16, and 18 and HA was shown by both direct binding assays.

Toxicity of Compounds in MDCK cells
Two-fold serial dilutions of selected compounds (starting at 25 µM) were incubated for 72 h at 37 °C with semi-confluent MDCK cells grown in 96-well tissue culture microplates.Under these conditions, no peptide affected cell viability up to the highest concentration.The ratio between the kinetic dissociation (k off ) and association (k on ) constants was used to determine the equilibrium dissociation constants (K D ).

Hemagglutination Inhibition Assay (HI)
The capability of compounds 4, 9, 16, and 18 to inhibit hemagglutinating activity was evaluated by HI.The inhibitory activity of peptides 1 and 2 is shown as a reference.The oseltamivir-resistant A/Parma/24/09 H1N1 virus strain was used.As shown in Table 2, compounds 4 and 16 prevented the hemagglutination of the oseltamivir-resistant strain A/Parma/24/09 H1N1 at quite similar concentrations to the reference peptides.

Neutralization Assay (NT)
Next, we evaluated the ability of the four compounds to inhibit viral infection by neutralization assay.All compounds were able to prevent viral infection in a dose range of roughly 50 to 140 pM, corresponding to selectivity indices ≈ 10 3 /10 5 .In particular, compound 4 inhibited virus replication at a concentration similar to that of the reference peptide (Table 3).

Computational Studies
The computational analysis of synthesized and tested compounds follows the previously generated protocol [46].Briefly, the hemagglutinin 3D structure was generated by homology modeling using the exact sequence of the assayed viral strains [46].The Receptor Binding Site, responsible for recognizing the sialic acid and highlighted by the SiteMap analysis on the HA surface, was selected to perform the docking calculations.Structure-based studies exploited the SP-peptide protocol available in Glide to explore the large conformational flexibility of the studied ligands.As already described, the 100 final poses were clustered to obtain each ligand's best representative binding geometry.
As all compounds showed activity toward A/Parma/24/09 H1N1, we focused our study on this viral strain's hemagglutinin.Additionally, we report the discussion about the most interesting compounds (4, 9, 16, and 18) binding to the RBS of A/Parma/24/09 H1N1.Docking score values for the four new ligands in comparison with previously studied tetrapeptides are reported in Table 4.The Glide calculated score is similar for all ligands and in line with the HI activity in Table 2. a Docking score values refer to the docking calculations performed in a previous article [46].
Compounds 16 and 18 are SKHS peptoids where the side chains of Lys2 and Ser4 residues, respectively, are moved on the peptide nitrogen atom.The docked binding mode follows the same profile we observed previously for the lead tetrapeptide SKHS (Figure S17): the ligands occupy the water-exposed RBS in an extended conformation, with Ser4 forming a network of hydrogen bonds with Gln210, Glu211, and Gly212.Differently from SKHS, the lysine side chain is able to form an ionic contact with Asp174 in both cases (Figure 4).The N-methyl derivative 9 presents reduced conformational flexibility that hampers the contact between the Lys2 protonated amine and Asp174 and causes a different ligand arrangement in the binding site.In fact, while the Ser4 OH occupies its conserved position in contact with Glu211, the serine residue in the opposite position reaches Asp174 to form an H-bond.Other H-bonds are formed with Thr177, Tyr80, Gln210, Glu211, and Tyr206 (Figure 5A,B).The N-methyl derivative 9 presents reduced conformational flexibility that hampers the contact between the Lys2 protonated amine and Asp174 and causes a different ligand arrangement in the binding site.In fact, while the Ser4 OH occupies its conserved position in contact with Glu211, the serine residue in the opposite position reaches Asp174 to form an H-bond.Other H-bonds are formed with Thr177, Tyr80, Gln210, Glu211, and Tyr206 (Figure 5A,B).Compound 4 is the derivative of SLDC where the Leu2 NH is methylated.Its docked interactions are driven by the hydrophobic contact of the Leu2 side chain in the most lipophilic region of the RBS.To satisfy this contact, the ligand inverts its binding mode with respect to the lead SLDC (Figure S18).The cysteine residue backbone forms H-bonds with key residues Gln210, Lys206, and Asp208, while Ser1 OH interacts with Val120 (Figure 5C,D).
Two treatments with piperidine/DMF (2:8, v/v) for 10 min and two extra 5 min treatments with piperidine/DBU/toluene/DMF (5:5:20:70, v/v) were used.Additionally, the chloranil test was carried out following each phase of deprotection and coupling to validate, respectively, the complete removal of the Fmoc protective group and to confirm that complete coupling has taken place on all the free amines on the resin.After removing the N-terminal Fmoc group, the peptides were acetylated by adding a solution of Ac 2 O/DCM (1:3) and agitating for 30 min.The final step was cleaving the peptides from the resin using the following cleavage cocktails: TFA/DCM (95:5 v/v) for 90 min (10 mL/g resin).By adding cold diethyl ether, peptides were precipitated, the solution was decanted, the solid was triturated, and the mixture was then again decanted.This procedure was carried out twice.The yield for the crude peptides was 30%.

Synthesis of Peptoids
A 2 M solution of bromoacetic acid in DMF (10 eq.) and DIC (10 eq.) was applied to the deprotected Rink Amide resin (100-200 mesh, 1% DVB, 0.59 mmol/g) in order to assemble N κ -protected N-N-aminoalkyl glycine residue on-resin.A 1 M solution of the selected amine in DMF was added after washing in DMF.In the Biotage Initiator + Alstra automated microwave synthesizer, the acylation reactions and nucleophilic displacement were carried out.
The beads were washed with DMF and prepared for the formation of the next residue.The finished peptoid resin was treated with a TFA-H 2 O-triisopropylsilane (TIS) mixture (95:2.5:2.5 by volume) for 90-120 min at rt to cleave the peptides from the resin and remove the acid-labile protecting groups.By adding cold diethyl ether, peptides were precipitated, then dried overnight under vacuum.A yield of 70-80% crude peptides was achieved.

Purification and Characterization
All crude peptides were purified by RP-HPLC on a preparative C18-bonded silica column (Phenomenex Kinetex AXIA 100 Å, 100 × 21.20 mm, 5 µm) using a Shimadzu SPD 20A UV/VIS detector, with detection at 214 and 254 nm.The column was perfused at a flow rate of 15 mL/min with solvent A (5%, v/v, water in 0.1% aqueous TFA), and a linear gradient from 5 to 90% of solvent B (85%, v/v, acetonitrile in 0.1% aqueous TFA) over 20 min was adopted for peptide elution.Analytical purity and retention time (tr) of each peptide were determined using HPLC conditions in the above solvent system (solvents A and B) programmed at a flow rate of 1.500 mL/min using a linear gradient from 5 to 90% B over 8 min, fitted with a C-18 column (Phenomenex, Aeris XB-C18 column; 150 mm × 4.60, 3.6 µm).All analogs showed >97% purity when monitored at 215 nm.Homogeneous fractions, as established using analytical HPLC, were pooled and lyophilized.Peptide molecular weights were determined by positive ESI infusion on a Q-Tof Premier Mass Spectrometer (Waters), equipped with the Xcalibur software (v4.2.47) for processing the data acquired.The sample was dissolved in a mixture of water and methanol (50/50) and injected directly into the electrospray source using a syringe pump at a constant flow (15 µL/min).Analytical data are shown in the Supplementary Materials (Table S1 and Figures S1-S16).

Microscale Thermophoresis (MST)
The Monolith NT.115pico (Nano-Temper Technologies GmbH, Munich, Germany) instrument was used for MST studies.His-HA was labeled using the Nanotemper His-Tag Labeling Kit RED-tris-NTA 2nd Generation as described elsewhere [47].His-HA protein (100 µL, 80 nM) in ddH 2 O was added to NT647-NHS fluorophore (100 µL, 40 nM, NanoTemper Technologies GmbH, Munich, Germany) in labeling buffer and incubated for 30 min at rt.Protein aggregates were removed from NT647-HA by centrifuging it at 15,000× g for 10 min at 4 • C. Pretests using standard-treated and premium-coated MST capillaries (NanoTemper Technologies GmbH, Munich, Germany) were performed as reported elsewhere.[46] In assay buffer (PBS), the protein did not adsorb to standardtreated capillary walls.Therefore, the following experiments were carried out on standardtreated capillaries.In order to determine the best state for MST signal repeatability and the suppression of unspecific adsorption to capillary walls, buffer conditions were then assessed.The assay buffer was used to dilute compound stocks (5 mM) in MST buffer to the greatest solubility concentration (50 µM).In MST tests, 16-fold 1:1 serial dilutions of each compound were added to NT647-HA (20 µL).The reaction solutions were loaded into standard-treated capillaries after 10 min of incubation at rt, and they were subsequently placed in the chip tray of Monolith NT.115 for thermophoresis analysis and the appraisal of K D values.Signals were captured using 10% LED power and high MST power.K D values were calculated from compound concentration-dependent changes in normalized fluorescence (F norm ) of HA after 21 s of thermophoresis.Each compound was tested in triplicate.Data examined using MO Affinity Analysis software (v2.3,NanoTemper Technologies GmbH, Munich, Germany).Confidence values (±) were indicated next to the K D value for each tested compound.

Surface Plasmon Resonance (SPR)
SPR binding studies were achieved at 25 • C using Biacore T200 (Cytiva, Uppsala, Sweden).His-HA was captured on the surface of the CM5 sensor chip using the antihistidine antibody provided in the His Capture Kit (Cytiva, Marlborough, MA, USA).The anti-histidine antibody was diluted to 50 µg/mL in the immobilization buffer included in the kit and covalently coupled to Sensor Chip CM5 by standard amine coupling to a level of approximately 11,000 RU.Then, His-HA was injected (21.5 µg/mL 10 mm acetate, pH 4.59) over the anti-histidine antibody surface for 60 s.On the reference surface, no protein was injected.The dissociation was monitored by injecting a running buffer for 600 s.Glycine buffer (10 mM, pH 1.5, 1 min) was injected to regenerate the surface.HBS-P+ buffer was used as the running buffer.After HA immobilization, the buffer was injected over the chip overnight (5 L/min).The association phase was carried out at 25 • C with a flow rate of 30 µL/min using a solution of peptide in HBS-P+ buffer at several concentrations (from 0.32 to 10 µM), and the dissociation phase was carried out with just the buffer.Using BIAevaluation software (v3.1) provided with the Biacore T200 instrument (Cytiva), the equilibrium dissociation constants (K D ), kinetic dissociation (k off ), and kinetic association (k on ) constants were calculated from the sensorgrams by global fitting of a 1:1 binding model.

Cells Line and VIRUS
Hemagglutination assays and/or plaque assays were used to determine the viral stocks' titers in accordance with accepted practices [48,49].The Madin-Darby canine kidney (MDCK, ATCC, CRL-2936) cell line was utilized.MDCK cells were grown as reported elsewhere [50].The Brisbane-like influenza A/H1N1 oseltamivir-resistant virus strain A/Parma/24/09 H1N1 was propagated in 80% confluent MDCK cells as previously described, and viral stocks were stored at −80 • C [51].

Cytotoxicity Assay
This assay was carried out on 80% confluent MDCK grown in 96-well tissue culture microplates in accordance with other reports [50].Two-fold serial dilutions of peptides (starting from 25 µM) were added to the cells and incubated at 37 • C. Light microscopy was used to analyze the cell morphology after 24, 48, and 72 h.Following this, cells were washed in phosphate-buffered saline (PBS, pH 7.4), and 100 µL of 0.5 mg/mL 3-(4,5dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (MTT) solution (0.5 mg/mL in PBS) was put to each well.After incubation at 37 • C for 3 h, the MTT solution was discarded and 100 µL of dimethyl sulfoxide (DMSO) were added to each well.Cells were incubated at room temperature for 15 min, and an ELISA plate reader (PerkinElmer Italia, Monza, Italy) with a 570 nm test wavelength and a 690 nm reference wavelength was used to read the microplates.Cytotoxicity was computed from an average of three replicates for each experiment, which was carried out in triplicate.The cytotoxicity was computed by comparing the average OD of wells containing compounds with the average OD of wells containing culture medium alone.Compound dilutions that did not affect cell viability compared to untreated cells were considered non-cytotoxic concentrations and utilized for antiviral assays.

Hemagglutination Inhibition Assay (HI)
First, the A/Parma/24/09 H1N1 virus was titrated using the hemagglutination test, exploiting its ability to agglutinate turkey red blood cells (RBC).To this end, 50 µL of the viral suspension were transferred to the 1st well of 24 of a 96-well U-bottom plate, and then 50 µL of PBS were added from well 2 up to well 24.Two-fold serial dilutions of the viral suspension (contained in the 1st well) were carried out in PBS from the 1st well to the 24th well.Then, 50 µL of turkey RBC suspension (0.5% in PBS) were added to each well, and the plate was incubated at room temperature (RT).After 1 h of incubation, the plate was tilted to let non-agglutinated RBC drip from the bottom of the wells.The reciprocal of the highest dilution of the viral suspension still capable of causing agglutination corresponded to the viral titer.This dilution represents 1 virus-agglutinating unit (HAU).The experiment was repeated six times, and the viral titer was obtained from the mean of six replicates.
A/Parma/24/09 H1N1 virus in phosphate-buffered saline (PBS, pH 7.4) was incubated for 1 h at 4 • C with serial dilutions of compounds in PBS.An equal volume of 0.5% turkey erythrocytes was then added and allowed to agglutinate for 60 min at RT. HI titers were expressed as the concentration of the compound producing 50% inhibition of hemagglutination by four virus-agglutinating units (HAUs).

Neutralization Assay (NT)
Cells were placed into wells of 96-well tissue culture microplates (20,000 cells/well) and incubated for 24 h at 37 • C in 5% CO 2 .NT was performed by incubating equal amounts of a viral suspension containing 105 plaque-forming units (p.f.u.) with serial two-fold compound dilutions, beginning at 25 µM, in culture medium for 1 h at 4 • C. A culture medium without compounds was used as a negative control.Cells were infected with 100 µL/well (three independent experiments were performed) of the virus-compound combinations.After 1 h of incubation, cells were washed carefully and incubated for 72 h at 37 • C. The viral cytopathic effect (c.p.e.) was measured by MTT assay as reported elsewhere by our laboratory [51].Results were expressed as a percentage of c.p.e.inhibition by comparison with untreated control cultures.EC 50 was defined as the compound concentration that generates 50% of the maximal protection.The percentage of protection was calculated by the equation (A570 compound-virus − A570 virus)/(A570 mock infected cells − A570 virus) × 100%.The data were expressed as the EC 50 relative to the value observed in the absence of the compound.
Neutralizing titers were expressed as EC 50 , that is the compound concentration providing 50% inhibition of infection.

Computational Studies
The in silico studies were carried out using the Schrodinger Suite, version 2023-1.The structure-based computational studies were performed on the HA homology models already generated and described in the previous work [46].
The studied peptide derivatives were built in Maestro using the Build facility and submitted to LigPrep calculations to obtain the correct protonation state at physiological pH.Prepared molecules were submitted to minimization using MacroModel, the OPLS4 force field, the GB/SA water solvation model, and the PRCG algorithm to a convergence of 0.01 kJÅ −1 mol −1 .
The SiteMap-calculated molecular interaction fields were used to define the Glide Grid for each HA in the Receptor Binding Site.Docking calculations were carried out using the Glide SP-peptide protocol and saving 100 poses for each ligand.Obtained poses were clustered using a Complete Linkage algorithm, defining the correct number of clusters suggested by the Kelley penalty evaluation.The most populated clusters and best docked compounds were selected for the following analysis.

Conclusions
In this paper, we synthesized a library of peptidomimetics, N-methyl peptides and peptoids, starting from the most promising tetrapeptides previously identified [23].Direct binding assays were carried out using two biophysical techniques (MST and SPR).Both demonstrated specific binding between compounds 4, 9, 16, and 18 and HA.Therefore, the antiviral activity of these compounds was assessed through biological assays, such as HI and NT.Results from the hemagglutination inhibition assay confirmed the specific interaction with HA.The neutralization test showed that S(N-Me)LDC (4) and SKHNhS (18) are the most active derivatives.Computational analysis allowed us to envision the putative interactions of these ligands with HA, where they establish conserved interactions with residues of the RBS, as already observed for peptides 1 and 2 [23].All applied methods agreed upon the identification of an N-methyl peptide, S(N-Me)LDC, able to bind HA with high affinity and inhibit influenza virus hemagglutination and cell infection at nanoand pico-molar concentrations, respectively.This small sequence can represent a valuable starting point for the design of small molecules.The work carried out opens the way to new perspectives for the development of new anti-influenza drugs, especially in a context in which the emergence of new and drug-resistant viruses highlights the need for new antiviral approaches and strategies.

Figure 1 .
Figure 1.Comparison of (a) a peptide and (b) a peptoid.

Figure 1 .
Figure 1.Comparison of (a) a peptide and (b) a peptoid.

Figure 1 .
Figure 1.Comparison of (a) a peptide and (b) a peptoid.

Figure 2 .
Figure 2. Direct binding measurements of the interaction between peptidomimetics 4, 9, 16, and 18 and HA (H1N1 hemagglutinin, Cat No.: Z03181, GenScript, NE).MST binding affinity assays were carried out as described in the Materials and Methods section (Section 3).The single dose-response curves for the binding of the peptidomimetics 4, 9, 16, and 18 to HA are depicted; Fnorm, normalized fluorescence.The experiments were repeated independently three times.The reported KD is the mean SD of three independent experiments.

Figure 2 .
Figure 2. Direct binding measurements of the interaction between peptidomimetics 4, 9, 16, and 18 and HA (H1N1 hemagglutinin, Cat No.: Z03181, GenScript, NE).MST binding affinity assays were carried out as described in the Materials and Methods section (Section 3).The single dose-response curves for the binding of the peptidomimetics 4, 9, 16, and 18 to HA are depicted; F norm , normalized fluorescence.The experiments were repeated independently three times.The reported K D is the mean SD of three independent experiments.

Figure 4 .
Figure 4. (Panels (A,C) Docked poses of peptoids 16 and 18 represented as sticks (orange and purple C atoms, respectively) in the RBS of A/Parma/24/09 H1N1 (deep-cyan cartoon).Just the HA residues that interact with the ligands are shown as lines for clarity's sake.H-bonds are represented as yellow dashed lines.(Panels (B,D)) Corresponding 2D ligand interaction diagrams; in these diagrams, residues close to the ligand are colored according to their properties (orange, negatively charged; blue, positively charged; green, hydrophobic; cyan, polar); H-bonds are shown as magenta arrows; and a gray shadow encircles atoms that are exposed to the solvent.

Figure 4 .
Figure 4. (Panels (A,C) Docked poses of peptoids 16 and 18 represented as sticks (orange and purple C atoms, respectively) in the RBS of A/Parma/24/09 H1N1 (deep-cyan cartoon).Just the HA residues that interact with the ligands are shown as lines for clarity's sake.H-bonds are represented as yellow dashed lines.(Panels (B,D)) Corresponding 2D ligand interaction diagrams; in these diagrams, residues close to the ligand are colored according to their properties (orange, negatively charged; blue, positively charged; green, hydrophobic; cyan, polar); H-bonds are shown as magenta arrows; and a gray shadow encircles atoms that are exposed to the solvent.

18 Figure 5 .
Figure 5. (Panels (A,C)) Docked poses of N-methyl peptides 9 and 4 represented as sticks (light orange and white C atoms, respectively) in the RBS of A/Parma/24/09 H1N1 (deep-cyan cartoon).Just the HA residues that interact with the ligands are shown as lines for clarity's sake.H-bonds are represented as yellow dashed lines.(Panels (B,D)) Corresponding 2D ligand interaction diagrams; in these diagrams, residues close to the ligand are colored according to their properties (orange, negatively charged; blue, positively charged; green, hydrophobic; cyan, polar); H-bonds are shown as magenta arrows; and a gray shadow encircles atoms that are exposed to the solvent.

Figure 5 .
Figure 5. (Panels (A,C)) Docked poses of N-methyl peptides 9 and 4 represented as sticks (light orange and white C atoms, respectively) in the RBS of A/Parma/24/09 H1N1 (deep-cyan cartoon).Just the HA residues that interact with the ligands are shown as lines for clarity's sake.H-bonds are represented as yellow dashed lines.(Panels (B,D)) Corresponding 2D ligand interaction diagrams; in these diagrams, residues close to the ligand are colored according to their properties (orange, negatively charged; blue, positively charged; green, hydrophobic; cyan, polar); H-bonds are shown as magenta arrows; and a gray shadow encircles atoms that are exposed to the solvent.

Table 1 .
Sequence and binding affinities between

Table 2 .
Sequence and HI activity of peptides 1

Table 3 .
Sequence and Nt activity of peptides and peptoids.
a EC 50 : the peptide or peptoid concentration that reduced cytopathic effect by 50%; b SI: the ratio of CC 50 (the peptide or peptoid concentration that reduced cell proliferation by 50%, equivalent to a concentration > 25 µM) and EC 50 ; the standard errors and means from three different experiments are displayed.The mean values of three independent experiments with standard errors are shown.c Scala et al. (2017) [23].