key: cord-0011061-2ie1v1wy authors: Ramírez-Salinas, G. Lizbeth; García-Machorro, Jazmín; Rojas-Hernández, Saúl; Campos-Rodríguez, Rafael; de Oca, Arturo Contis-Montes; Gomez, Miguel Medina; Luciano, Rocío; Zimic, Mirko; Correa-Basurto, José title: Bioinformatics design and experimental validation of influenza A virus multi-epitopes that induce neutralizing antibodies date: 2020-02-14 journal: Arch Virol DOI: 10.1007/s00705-020-04537-2 sha: f13a7b0b45330b36e40b686d6789935e9481bd2a doc_id: 11061 cord_uid: 2ie1v1wy Pandemics caused by influenza A virus (IAV) are responsible for the deaths of millions of humans around the world. One of these pandemics occurred in Mexico in 2009. Despite the impact of IAV on human health, there is no effective vaccine. Gene mutations and translocation of genome segments of different IAV subtypes infecting a single host cell make the development of a universal vaccine difficult. The design of immunogenic peptides using bioinformatics tools could be an interesting strategy to increase the success of vaccines. In this work, we used the predicted amino acid sequences of the neuraminidase (NA) and hemagglutinin (HA) proteins of different IAV subtypes to perform multiple alignments, epitope predictions, molecular dynamics simulations, and experimental validation. Peptide selection was based on the following criteria: promiscuity, protein surface exposure, and the degree of conservation among different medically relevant IAV strains. These peptides were tested using immunological assays to test their ability to induce production of antibodies against IAV. We immunized rabbits and mice and measured the levels of IgG and IgA antibodies in serum samples and nasal washes. Rabbit antibodies against the peptides P11 and P14 (both of which are hybrids of NA and HA) recognized HA from both group 1 (H1, H2, and H5) and group 2 (H3 and H7) IAV and also recognized the purified NA protein from the viral stock (influenza A Puerto Rico/916/34). IgG antibodies from rabbits immunized with P11 and P14 were capable of recognizing viral particles and inhibited virus hemagglutination. Additionally, intranasal immunization of mice with P11 and P14 induced specific IgG and IgA antibodies in serum and nasal mucosa, respectively. Interestingly, the IgG antibodies were found to have neutralizing capability. In conclusion, the peptides designed through in silico studies were validated in experimental assays. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (10.1007/s00705-020-04537-2) contains supplementary material, which is available to authorized users. Influenza A virus (IAV) is a lipid-enveloped, singlestranded, negative-sense RNA virus belonging to the family Orthomyxoviridae. The viral envelope contains three transmembrane proteins (NA [neuraminidase], HA [hemagglutinin] and M2 [proton channel]) on the viral surface and one protein (M1 [matrix protein]) below the membrane. The viral core contains the nucleoprotein (NP), viral RNA, and three polymerase proteins (PB1, PB2, and PA) [1] . IAV is classified into subtypes based on two major antigens: the surface spike glycoproteins NA and HA [2] . All IAV subtypes are known to cause infections in birds, which are their natural reservoir [3] . Humans are infected principally by the IAV subtypes H1N1, H2N2, H3N2, H7N9, and H5N1 [4] . Influenza pandemics have become serious socioeconomic and public-health problems worldwide. Moreover, seasonal flu causes approximately 250,000 to 500,000 deaths per year [5, 6] . IAV epidemics and pandemics are attributed to 1 3 mutations in the viral RNA genome. Mutations involving surface proteins (NA and HA) result in structural protein changes that cause a loss of antibody recognition against the virus. This is one reason why new flu vaccines need to be designed for each seasonal influenza or pandemic influenza strain. The development of vaccines is the major method used to prevent IAV infection and represents one of the most important contributions by the immunology field to public health [7] . An important strategy is to identify conserved epitopes that may be used to design new vaccines that are capable of conferring broad protection. Currently, the primary goal is to develop vaccines that protect by eliciting antibody responses against multiple subtypes and strains of influenza viruses [8] [9] [10] . These broadly neutralizing antibodies (bnAbs) generally target conserved and functional regions or epitopes on the major surface glycoproteins: hemagglutinin (head and stem), neuraminidase, and M2e [10] [11] [12] . The hemagglutinin (HA) is the main surface glycoprotein of influenza virus, which mediates the adsorption and penetration of the virus into host cells [13] . Each molecule of HA comprises a membrane distal globular head composed of HA1, which contains the receptorbinding site (RBS), and a stem region, which encompasses the fusion machinery [14] . Most bnAbs are directed against the HA protein. The receptor-binding site is a functionally conserved region on the HA1 globular head domain that is a target for bnAbs that inhibit viral entry by preventing HA binding to its host receptor [15, 16] . Since the stem region contains the most conserved epitopes for antibody recognition, antibodies produce against this region have a higher neutralization breadth than RBS-targeted bnAbs. These stem-binding bnAbs inhibit virus replication by blocking attachment and preventing conformational changes that are essential for membrane fusion [15] [16] [17] . NA is the second most abundant glycoprotein on the surface of influenza A and B viruses, and conserved domains or epitopes in NA induce bnAbs that protect against viruses of a single subtype [17] . Thus, NA epitopes may use in universal influenza vaccines [12, [17] [18] [19] [20] . Although NA-specific antibodies can control infection by several mechanisms, the main mechanism is the inhibition of enzyme activity [12, 18, 21] . Thus, "universal vaccines consisting of conserved domains or epitopes from HA and NA could be more broadly protective than single proteins. The rational design of peptide-based vaccines is based on computational procedures that employ knowledge of the antigen recognition of different protein targets, such as the major histocompatibility complex (MHC), T-cell receptor, and/or B-cell receptor. These procedures are performed at the molecular level [22] using information from the Protein Data Bank. This strategy offers many advantages, such as the production of short immunogenic peptides without the risk of infection and safe handling of biological samples under room-temperature conditions. Additionally, peptides can be synthesized in vitro to decrease production costs and can be handled without strict conditions due to their stability [23] . Peptide vaccines are classified according to the immunological cell type involved in the immune response (i.e., B-cell or T-cell epitopes) and are used to induce specific immune responses. T cells recognize peptides coupled to MHC from antigen-presenting cells, and the MHC is then coupled to the T-cell receptor (TCR). Following the formation of the peptide-MHC-TCR complex, T cells receive several biological signals that are capable of activating the immune system to act against microorganisms or cancer cells [24] . In addition, activation of B cells by epitopes can increase the protection spectrum of vaccines, making them more protective than those with limited immune system activation capacity [25] . In this study, we selected a set of peptides (P) from HA and NA (including hybrids) that were capable of generating antibodies against IAV. First, a protein sequence bioinformatics study was performed using the HA and NA proteins to obtain multi-epitope peptides characterized by conserved residues and promiscuous properties. The epitope selection process included criteria related to their protein surface localization. Then, three-dimensional (3D) models of the target epitopes were constructed and subjected to molecular dynamics (MD) simulations to explore their thermodynamic properties. Finally, the epitopes were used in experimental assays to explore their immunogenic properties in animal models. We retrieved the full primary sequences of the HA and NA proteins reported between 1918 and 2014 from the NCBI influenza virus resource [26] . We aligned the HA and NA sequences for each subtype using the Muscle/EBI server [27] and the Clustal X 2.0.11 program [28] . We constructed a three-dimensional (3D) structure model based on the NA consensus sequence (YP_009118627.1) of the H1N1 virus by employing the Swiss Model server [29] [30] [31] [32] . The NA monomer A (PDB ID: 3TI4) was used as a template. For the H1N1 HA consensus sequence (Table 1) , we used the Modeller 9.10 program to build 3D models of multi-chain HA by employing a multi-trimer template HA (PDB ID: 1RUY). The stereochemical quality of the models was evaluated using Ramachandran plots and PDBSUM [33] (data not show) and ERRAT servers [34] (data not shown). We submitted the NA (YP_009118627.1) and HA consensus sequences (Table 1) We selected epitopes based on promiscuity criteria among the IAV subtypes, the degree of conservation, and their surface localization on the quaternary structure models of the proteins that were built (Fig. 1S ), selecting six peptides (P1, P2, P3, P4, P5 and P6) without conjugation to keyhole limpet hemocyanin (KLH). P1 is located in the middle of HA, whereas P2-P4 and P6 are located in globular heads of HA. P5 is located in in globular head of NA (Fig. 1S) . We determined the degree of exposure of the regions by visual inspection, using the VMD program [35] . Once the most promising epitopes were identified, we added linkers of different length and structure (6-9 glycine residues) to control the distance between two selected epitopes and to maintain multi-epitope properties, discarding hybrids in which the added glycines or the neighboring peptides were predicted to be immunogenic. Additionally, epitopes previously identified by our research group were selected as candidates for the construction of hybrid peptides to evaluate their capacity to increase immunological responses [19, 36] . We predicted proteasome and immunoproteasome cleavage using the Proteasome Cleavage Prediction Server, PCPS (http://imed.med.ucm.es/Tools /pcps/index .html), to verify epitope structural stability. The PCPS analysis identified fragments of peptides ranging from 9 to 21 residues in length. These peptide sequences were submitted to the NetMHC 3.2 [37, 38] , NetMHCII 2.2 [37, 39] and PRED BALB/C (http://antig en.i2r.a-star.edu.sg/predB albc/.) servers to determine whether they maintained their immunogenic properties. We built 3D models of the selected multi-epitopes using the PEPstrMOD server (http://osddl inux.osdd.net/ragha va/ pepst rmod/, accessed on March 2018). We subjected these multi-epitopes to MD simulations using the NAMD 2.5 program [40] with CHARMM 27 as the force field [41] . The constant number-of-particles, pressure, and temperature (NPT) ensemble and the periodic boundary conditions were used. A constant temperature (310 K) was maintained using a Langevin thermostat set at a constant pressure (1 atm) maintained by an isotropic Langevin barostat [42] . These pressure and temperature values represented standard physiological conditions. Energy minimization was performed with 1000 steps in the conjugate-gradient algorithm with restraints on the peptide backbone, followed by 1000 steps without restraints. The system was heated for 20 ps and equilibrated for 60 ps with restraints to the α-carbon and then finished with no restraints and 310 K. These MD simulations were run for 500 ns. We calculated the root mean square deviation (RMSD), which measures whole-protein motions while generally considering the α-carbon coordinates in relation to MD simulation times, root mean square fluctuations (RMSF), which measure the α-carbon coordinates in relation to MD simulation time per residue, and radius of gyration (Rg), which measures the protein compactness from the center of the protein to the periphery according to the 3D coordinates in the structural analyses of multi-epitope MD simulations from the Carma64 program [43] . The secondary structure was determined using the Visual Molecular Dynamics (VMD) program [35] . Nine peptides were designed, but eight peptides (including P6 in P14) were selected (Table 1) to be synthesized (Peptide 2.0 Inc.) according to theoretical studies for experimental testing as described in the following experimental procedures ( Table 2 ). The conjugation of peptides with the KLH molecule was performed using succinimidyl 4-(N-maleimidomethyl)cyclohexane-1-carboxylate (SMCC), which is an heterobifunctional crosslinker that contains N-hydroxysuccinimide (NHS) ester and maleimide groups that allow covalent conjugation of amine molecules in KLH with the sulfhydryl groups of cysteine residues located in the middle of each hybrid peptide. We immunized eight rabbits (from Harlan Mexico) with peptides (P1, P2, P3, P4 and P5) without conjugation to KLH and hybrid peptides (P9, P11 and P14) conjugated to KLH at three time points with during seven-day intervals. For the primary immunization, 400 µg of peptide plus complete Freund's adjuvant (Sigma Chemical Co.) was administered by the subcutaneous route. The second immunization (7 days later) included 300 µg of peptide plus incomplete Freund's adjuvant administered by the subcutaneous route. The third immunization was administered by the intramuscular route 15 days after the primary immunization and contained 300 µg of peptide in 5 mL of saline solution. Seven days after the last immunization, we anesthetized the rabbits with pentobarbital via the intraperitoneal route. Serum samples were obtained from blood extracted by cardiac puncture and stored at -70 °C. Screening was carried out to select the peptides that induced the highest antibody titer in rabbit serum. For this purpose, five groups of peptides were used. Six BALB/c mice for each peptide were injected three times at seven-day intervals via the intranasal route with 100 µg of the peptide (P1, P4, P5, P11 and P14) plus 2 µg of cholera toxin (CT) mucosal adjuvant. Seven days after the last immunization, we anesthetized the mice with pentobarbital. Samples were obtained from blood extracted by cardiac puncture and washes of the nasal cavity with PBS, respectively. This study was carried out in strict accordance with the recommendations in the Guide for the Care and Use of Laboratory Animals, in accordance with Mexican Federal regulations for animal use and care (NOM-062-ZOO-1999) and was authorized by the Comité Interno del Cuidado y Uso de los Animales de Laboratorio (CICUAL) at ESM-IPN. The protocol was approved by the Committee on Ethics of the Escuela Nacional de Medicina y Homeopatia of Instituto Politecnico Naciocal, IPN (protocol number ENMH-CB-0095-2014). Mice were anesthetized with ether for intranasal immunization, and for biological assays, mice and rabbits were killed using pentobarbital. All efforts were made to minimize suffering and to use biological samples immediately. We determined anti-peptide antibody levels in serum samples from immunized rabbits and serum samples or nasal washes of mice immunized with peptides, using an indirect enzymelinked immunosorbent assay (ELISA). For biological assays, the mice or rabbits were killed using pentobarbital, and their biological samples were used immediately. Each sample was tested in duplicate. Briefly, 96-well plates were coated with 100 µL of different peptides (20 µg of peptide per mL in carbonate-bicarbonate buffer [15 mM Na 2 CO 3 and 35 mM NaHCO 3 , pH 9.6]). The plates were incubated overnight at 4 °C, followed by three washes with phosphate-buffered saline containing 0.05% Tween-20 (PBST). Blocking was performed using PBST plus 6% fat-free milk for 1 h, and the samples were then washed with PBST. Serial dilutions of serum samples from nine immunized rabbits (1:100 to 1:25000), serum from immunized mice (from 1:100 to 1:3200) and nasal washes (from 1:1 to 1:8) were prepared and added to the plates, which were incubated overnight at 4 °C and washed with PBST, after which 50 μL of a 1:3000 dilution of goat anti-rabbit IgG peroxidase (Thermo Scientific) for rabbit serum and rabbit anti-mouse IgG peroxidase (1:3000) or rabbit anti-mouse IgA (1:500) peroxidase (Thermo Scientific) for mouse samples was added to each well. The plates were incubated for 4 h at room temperature and washed with PBST. The enzymatic reactions The endpoint titer was defined as the reciprocal of the highest analyzed dilution that gave an absorbance value above 0.5. Influenza A virus strain Puerto Rico/916/34 was grown in Madin-Darby canine kidney (MDCK) cells. The cells were cultured in Dulbecco's modified Eagle medium nutrient mixture F-12 (DMEM/F-12) (Caisson) supplemented with 5% fetal bovine serum (SFB) (Mediatech) at 37 °C and 5% CO 2 . The viral stock was prepared in serum-free medium and frozen at -20 °C for later use as described previously [36] . Viral particles that had been isolated from infected MDCK cells and purified using serum-free medium were trapped in wells coated with concanavalin A (Con A), which immobilizes the detergent-solubilized viral glycoproteins of the viral envelope. Briefly, 96-well plates (Costar) were coated with 100 μL of Con A (Sigma-Aldrich) per well at a concentration of 50 μg/mL in PBS, pH 7.4, for 1 h. The wells were washed three times with PBS containing 0.1% Triton X-100 (PBS-TX) and incubated with solubilized and inactivated influenza A virus strain Puerto Rico/916/34 (serum-free virus stock) in PBS-TX for 1 h. After the wells were washed with PBS-TX, the unreacted ConA binding sites were blocked with RPMI medium 1640 containing 10% FBS (RPMI-10% SFB) for 1 h. Dilutions of the heat-inactivated rabbit serum samples were made in RPMI-10% SFB and incubated for 1 h at room temperature. The positive control was a 1:500 dilution of serum from a patient who was positive for influenza virus [21] . Rabbit pre-immune serum samples were used as negative controls. The wells were washed again and incubated with 100 μL of the appropriate peroxidase-conjugated anti-IgG (H+L) (Santa Cruz) as the secondary antibody. The wells were washed again and then incubated for 1 h with 100 μL of 2, 2'-azino-bis (3-ethylbenzothiazoline-6-sulfonic acid) (ABTS) (Sigma-Aldrich), and H 2 O 2 was added as a substrate. Absorbance values were determined at 405 nm [36] . Levels of antibody against HA and NA were determined using recombinant HA (rHA) proteins purchased from Serum samples from immunized rabbits were inactivated by heating at 56 ± 2 °C in a water bath for 30 min, and the total IgG were purified by affinity chromatography using immobilized protein A Sepharose CL-4B (Sigma-Aldrich) as described by Contis-Montes de Oca et al. [44] . Dilutions were made in PBS and added to a final volume of 25 μL in 96-well round-bottom plates using pre-immune serum as a control. The virus concentration was then adjusted to 4 hemagglutination units (HAU) per 25 μL in PBS and added to samples, except for the negative controls, which consisted of PBS without virus. The plates were incubated for 1 h at 37 °C while human red blood cells (HRBCs) of type "O" were washed twice with 1x PBS, pH 7.2, and centrifuged at 200 × g for 10 min at 4-8 °C. The cells were then resuspended in a final concentration of 0.75% in 1x PBS, and 50 μL of this suspension was added to each well. The plates were incubated again for 1 h at 37 °C, and the effects were observed [45] . The HI titer was the highest dilution that caused total inhibition of agglutination. HRBCs with PBS were used as a negative control, and HRBCs plus virus at 4 HAU were used as a positive control. Sera were tested at dilutions from 1:12 to 1:7680. Neutralizing antibodies were titrated as described previously [46, 47] Briefly, serial dilutions of the heat-inactivated test sera in duplicate starting from 1:100 were mixed with 20 PFU of influenza A virus strain Puerto Rico/916/34 for 1 h at 37 °C and added to MDCK cells at a density of 80,000 cells per well in 24-well plates. Positive and negative control sera and virus back titration to confirm the viral inoculum were included. At 1 h after infection, serum-free DMEM/F-12 medium with 2 µg of TPCK (L-1-tosylamide-2-phenylethyl chloromethyl ketone)-treated trypsin was added to each well, and the plates were incubated for 1 h at 37 °C. Then, 500 µL of overlay medium was added to each well, and the plates were incubated for 3 days under the same conditions. The supernatant was discarded, and plaques were visualized by staining with naphthol blue-black dye. The result is reported as the reciprocal of the dilution that protects 50% of the cells from a cytopathic effect due to infection. Three independent assays were conducted in duplicate, and the results are presented as the mean and standard deviation of one representative assay. The statistical significance of differences among groups was determined by ANOVA, followed by Tukey's test. Differences with p values less than 0.01 were considered significant. The analysis was performed using the PRISM program (GraphPad, San Diego, CA.). The half-maximal effective concentration (EC 50 ) was determined using the statistical program SigmaPlot for Windows version 11 (Systat Software Inc., San Jose, CA, USA). (Fig. 1A-D) of the influenza A H1N1 glycoproteins NA and HA (Fig. 1E) were obtained from multiple alignments. These consensus protein sequences were used to identify the immunogenic regions in HA and NA (Fig. 1) . Sequences that were conserved during evolution (data not shown), present in different influenza virus subtypes, or located on the surface (exposed to solvent, see Fig. 1S ) were predicted to be easily accessible to antibodies with neutralizing potential [48] , as demonstrated for HIV [49] . Six peptides (from H1N1 as reference), including four from HA (P1, P2, P3, P4 and P6), one from NA (P5), and one from NA reported previously by us [19] were selected as described above. Multiple sequence alignments of HA peptides from different subtypes of influenza A virus (H1N1, H7N9, H5N1, H2N2 and H3N2) showed a degree of conservation for each of the selected peptides when H1N1 (A/California/07/2009) was used as a reference (100% identity, Fig. 1 ). P1 was conserved to different degrees with respect to subtypes H1N1 (A/Puerto Rico/8/1934), H7N9 and H5N1 (Table 3) . The most similarity was observed with H1N1 (A/Puerto Rico/8/1934) and the lowest with H3N2 (see Fig. 1A ), suggesting that P1 could induce immunogenic responses against H1N1, H2N2, H3N2, H7N9 and H5N1 as a multi-vaccine. P2 and P4 were most similar to the corresponding sequences of pandemic H1N1 (A/Puerto Rico/8/1934), see Table 3 ), suggesting that they were potential candidates for specific vaccines and immunodiagnostics, whereas P3 was most similar to the corresponding sequence of H1N1 (A/Puerto Rico/8/1934) and had a low degree of similarity to H3N2 sequences. However, multiple alignment of NA sequences showed that P5 was more similar to the corresponding sequences of the H1N1 (A/Brisbane/59/2007 [H1N1]) and H5N1 subtypes than to those of other influenza virus subtypes (see Table 3 ). Furthermore, P1, P3 and P5 could work as a multi-vaccine epitope to provide broad protection against unspecified influenza A virus infections. One additional parameter that the epitopes needed to achieve for good recognition by MHC and TCR was low structural mobility [50] . This was investigated using MD simulations [50] measuring structural motions dependent on structural flexibility of hybrid peptides with long sequences (Fig. 2) . It is known that lower flexibility of peptides is associated with a better immune response [36] . As shown in Fig. 2 , P9 yielded larger RMSD values in MD simulations than compared to P11 and P14 (approximately 125-200 ns and higher fluctuation close to 225 ns) ( Fig. 2A) . From 275 ns to 500 ns, P11 showed less motion with minimal oscillations compared to P9 and P14, suggesting that it might be more immunogenic. Based on RMSF values, P9 exhibited more mobility than P11 and P14. Additionally, P9 had higher Rg values (approximately 175 ns) compared to P11 and P14. These results agree with some of the biological results described in the following paragraphs. The peptides that bind to HLA-DR, HLA-DQ, HLA-DP, human MHC class II alleles and mouse I-E d and I-A d molecules are shown in Table 4 . All of the peptides contain identical or very similar sequences (epitopes) that are pre- dicted to bind to both human and mouse MHC-II alleles. For example, P1 is an epitope that binds to human HLA-DRB1*0101, HLA-DRB1*0701, HLA-DRB1*1302 alleles and mouse I-A d allele. In summary, the selected peptide contains epitopes that are capable of interacting with human and mouse (BALB/c) MHC-II alleles. These results show that the predicted peptides contain helper T-cell epitopes, which could activate T-cell-dependent responses, such as IgG and IgA antibody responses. The peptides P2, P3 and P4 contain three epitopes, whereas P1 only contains one epitope. The hybrid peptides P9 and P11 contain four epitopes, and P14 contains only one that is likely to induce an immune response against HA and NA proteins [51, 52] . Table 5 . These peptides are predicted to contain several T-cell epitopes, some of which are predicted to bind to both human and mouse MHC-I alleles. For example, P1 contains the sequence FQNIHPITI, which binds to human HLA-B*3901 and mouse H2-D d alleles. The peptides were found to contain different number of epitopes that are recognizable by cytotoxic T cells. P1, P9 and P11 have four epitopes, and P2 and P14 have three. The epitopes that bind to MHC-I molecules could activate cytotoxic T cells that contribute to protection against invasion by influenza virus [51] [52] [53] . Using the B-cell epitope prediction program ABCpred, we identified several linear B-cell epitopes (Table 6 ). P1, P3 and P5 were predicted to contain two B-cell epitopes, whereas the hybrid peptides were predicted to contain 4-5 epitopes. Therefore, hybrid peptides could induce stronger mucosal and systemic humoral immune responses [54] . To confirm the immunogenicity of these peptides in vivo, we analyzed sera obtained from rabbits, and the peptides that induced the highest antibody titers in rabbits were used to immunize mice intranasally. IgG in sera and IgA in nasal washes, as well as neutralizing antibodies, were then detected. The unconjugated P1 induced the strongest IgG immune response in rabbits (Fig. 3A) . This result means that P1 could be immunogenic because it induced IgG antibodies in the absence of KLH [19] . According to epitope predictions, P1 contains T-cell (MHC II) and B-cell epitopes and is located in the middle of protein with 360º of exposure. P1 was included in peptides P9 and P11. P1 did not induce an immunogenic response when it was hybridized to P2 to yield P9; however, it induced a stronger IgG immune response when it was hybridized to P5 (NA) to yield P14. As expected, the peptides P11 and P14 induced the strongest immune response when conjugated to KLH. However, P9 induced a poor immune response (Fig. 3 ). P5 hybridized with P1 to yield P11 resulted in an immune response that could be attributed to P1. Although P1 is present in P9, it is located far from KLH. The position of the peptide in the hybrid could therefore explain the experimental results. P1 in the P9-KLH complex could be located at the KLH protein surface, and this location makes proteasome processing for presentation by MHC more difficult than the P11-KLH complex. The better immunogenic response obtained with P1 alone or P1 hybridized with P5 to yield P11 might have been due to the fact that P1 is very rich in basic residues, which have been shown to be important for immunogenicity [36] . P1 was capable of inducing an immune response despite having a lower molecular weight than P11 and P14. Similar results have been reported previously for another peptide [19] . As reported previously, the immunogenicity of P14 could be due to its NA epitope [19] . As mentioned above, we tested the immunogenicity of peptides conjugated to KLH to investigate the effect of KLH on the response to P9. Surprisingly, P9 still induced a weaker immune response than P11 despite the presence of P1 (Fig. 3, panel B) . These results confirmed that the peptide sequence was of key importance for immunogenicity. The carriers could play an important role due to their ability to affect the recognition pattern on MHC regardless of the presence of basic residues that are capable of recognizing important sites, as reported by Carrillo-Vazquez et al. [36] . The IgG immune response is a T-cell-dependent response because the class switch of the IgG isotype requires the interaction of peptides with MHC II molecules and T helper (CD4 T) cells [55, 56] . Thus, our results indirectly demonstrate that immunogenic peptides induce T helper cells restricted to MHC-II. Using ELISA, the immunogenicity of peptides in the mucosal compartment was examined in mice that were immunized intranasally with peptides mixed with CT (Fig. 4) . The IgG and IgA antibody responses to the peptides were similar in the serum samples. Intranasal immunization induced higher titers of IgG and IgA antibodies to P11, lower titers to P14, and very low titers to the unconjugated peptides P1, P4, and P5 (Fig. 4, panels A and B) . These results suggest that P1 requires other structural changes for proteasome protection or immunological transport in the P11 hybrid. To investigate the antibody responses in the mucosal compartment, the levels of IgG and IgA antibodies against the peptides were measured in nasal washes (Fig. 4 , panels C and D). We observed that P11, P1 and P14 induced higher titers of IgG antibody, whereas the IgA antibody response was highest against P1 and P4, followed by P11, P14 and P5. Surprisingly, all peptides (both conjugated [P11 and P14] and unconjugated [P1, P4, and P5] peptides) induced were added to microplates that had been coated with different peptides. IgG antibodies were detected using a secondary Ab specific for rabbit IgG (1:3000). Sera from rabbits immunized with either the unconjugated peptide 1 or with hybrid peptides P11 and P14 conjugated with KLH showed the highest titer (1:800), whereas sera against the remaining peptides had a lower titer (A 490 <0.5; 1:800). Samples of serum from a pre-immune rabbit were used as controls. No peptide-specific IgG responses were found in pre-immune sera. Individual samples were run in duplicate, and the data are shown as the mean ± SD high IgA antibody titers in nasal washes, but peptide 1 yielded the highest titers. The production of IgA antibodies in nasal mucosa requires that dendritic cells process and present peptides on class II MHC molecules (pMHC) to CD4 + T cells. In turn, these helper cells promote IgA class-switching recombination and affinity maturation of IgA-committed B-cells [57] . Thus, our results indicate that the peptides induced CD4 + T helper cells that recognize peptides in the form of an MHC-II-peptide complex [55, 57, 58] . In summary, intranasal immunization of mice with P11 and P14 along with CT induced specific IgG and IgA antibodies in serum and nasal mucosa. These antibodies could protect the host against virus invasion. According to the current concept, IgA is the dominant antibody involved in protection against infection in the nasal compartment (upper respiratory tract), whereas serum IgG, which diffuses into the lower respiratory tract, is predominantly involved in the protection of the lungs [59] [60] [61] . Finally, to test for substantial heterosubtypic binding activity, sera against P11 and P14 were tested for reactivity with purified hemagglutinin proteins from different subtypes and strains. We focused on testing P11 and P14 due to the promising results obtained with these peptides. Peptide-specific IgG or IgA antibody titers in sera and nasal washes of mice immunized with different peptides. Serial dilutions of sera (panels A-B) or nasal washes (panels C-D) from mice inoculated with different peptides were added to microplates that had been coated with specific peptides. Antibodies were detected using a sec-ondary Ab specific for mouse IgG (1:3000) or IgA (1:500). Samples of serum from pre-immune mice were used as controls, and no peptide-specific IgG or IgA responses were found in pre-immune sera or nasal washes, respectively. Individual samples were run in triplicate, and the data are shown as the mean ± SD As shown in Fig. 5A -G, the response patterns of the sera anti-P11 and anti-P14 against the hemagglutinins was similar for most hemagglutinins, with the exception of A/Brisbane/59/2007 (H1N1) and A/Puerto Rico/8/34 (H1N1), both of which are H1N1 strains (Fig. 5F and G) . Both sera reacted similarly against the H7N9 subtype. The responses of P11 and P14 against the H7N9 subtype were not clear from the bioinformatics studies (see Fig. 1 ). The linear epitopes could not be evaluated due to their low similarity to the H7N9 protein sequence, suggesting that a conformational epitope was involved. Conformational epitopes are capable of inducing immune responses despite being discontinuous epitopes and are recognized by antibodies as linear epitopes in a tertiary structure [62] . Antibodies produced in response to the A/California/04/09 pandemic virus cross-reacted with several influenza virus strains, including A/Brisbane/59/2007 (H1N1), A/Puerto Rico/8/34 (H1N1), A/Aichi/2/1968 (H3N2), and A/Anhui/1/13 (H7N9) [63] [64] [65] [66] . The origin of these crossreactive antibodies may be cross-reactive T and B cells specific for conserved epitopes in the hybrid peptides P11 and P14. Various studies have demonstrated that cross-protection against influenza virus strains may be due to cross-reactive T cells [67] [68] [69] . A previous in silico analysis provided a list of potential cross-reactive T-cell epitopes, including the sequences CPKYVKSTK and HAGAKSFYKNL, which were present in P1, P11 and P14 [70] . Also, as shown in Fig. 6 , P11 and P14 induced antibodies capable of . Serial dilutions of sera from rabbits inoculated with peptides con-jugated to KLH were added to microplates that had been coated with purified recombinant hemagglutinins (0.5 μg). Antibodies against P11 and P14 had similar binding activity to the majority of the rHAs, except for A/Brisbane/59/2007 and A/Puerto Rico/8/1934, both of which are H1N1 subtypes. As controls, serum samples from preimmunized rabbits were tested against the specific peptide recognizing purified NA protein from a viral stock in ELISA studies measuring IgG antibodies. According to ConA envELISA studies, P11 and P14 are capable of inducing antibodies in rabbits that recognize wild-type viral antigens (immobilized with ConA) from an influenza virus H1N1 strain at dilutions of 1:200 and 1:800 ( Fig. 7 ) with significant statistical differences (one-way ANOVA, p < 0.0001) compared to pre-immune sera at the same dilutions. These absorbance values are equal to values others have reported for patients who were previously infected with influenza virus [36] . In the hemagglutination inhibition assay, it was found that sera from rabbits that had been inoculated with P11 and P14 inhibited hemagglutination at a dilution of 1:384 and 1:960, respectively (Table 7) . There are reports in which viral antigens (surface glycoproteins) can be immobilized to be recognized by antibodies in an ELISA, as reported elsewhere for influenza virus and dengue virus [36, 71] . In this study, antibodies from rabbits immunized with P11 and P14 recognized viral particles of strain (A Puerto Rico/916/34) immobilized with ConA (Fig. 7) , and they also inhibited hemagglutination with titers higher than 40 (Fig. 7) , which decreases the risk of infection to 50% [72] . The above coincides with the neutralizing capacity of antibodies induced in mice immunized with peptides 11 and 14 (reciprocal of PRNT 50 titer of 380 and 1600, respectively; Table 7 ). These results indicate that in silico studies can be used to identify potential protective epitopes. In summary, we compared amino acid sequences ( Fig. 1) to provide a plausible explanation for the observed cross-protection between 2009 H1N1 and other influenza viruses [64] . Only P1 alone produced a response in vivo. The observed responses to P11 and P14 could have been due to coupling to KLH. Additionally, P11 contained P1 and P14 contained a recently reported peptide [19] . However, the responses against H7N9 were apparently not due to the presence of linear epitopes, suggesting that the peptides were capable of adopting 3D structures that formed a conformational epitope [62] . Rabbit antibodies to P11 and P14 recognized HA proteins from both group 1 (H1, H2, and H5) and group 2 (H3 and H7) influenza A viruses belonging to three different clades: the H1 clade (H1, H2, and H5), the H3 clade (H3), and the H7 clade (H7). Additionally, the antibodies also recognized purified neuraminidase from the viral stock and viral particles immobilized in the trapping ELISA, and they were able to inhibit hemagglutination and neutralize the virus. NA-specific IgG antibody titers in sera from rabbits immunized with different peptides. Serial dilutions of sera from rabbits inoculated with peptide 14 (P14), and peptide 11 (P11) were added to microplates previously coated with neuraminidase protein (0.5 µg). IgG antibodies were detected using a secondary Ab specific for rabbit IgG (1:3000). Individual samples were run in duplicate, and the data are shown as the mean ± SD Fig. 7 Trapping ELISA. Antibodies from rabbits immunized with P11 or P14 recognized influenza virus immobilized with ConA. Preimmune serum samples were used as negative controls, and a human serum was used as a positive control (p < 0.0001, pre-immune vs. immunized rabbit) The humoral immune responses against HA and NA are the primary means of resistance to influenza virus infection [68] ; therefore, hybrid peptides composed of HA and NA peptides can be used as potent synthetic vaccines. Because antibodies against the hybrid peptides P11 and P14 had broad heterosubtypic activity against the antigenically diverse H1, H2, H3, H5, and H9 influenza subtypes, these peptides might be used as broad-spectrum agents for prophylaxis (to induce neutralizing antibodies) and also for treatment of human or avian influenza infections, since the antibodies recognize the neuraminidase protein and could potentially prevent the release of new viral progeny during infection. The use of bioinformatics tools enabled the identification of potential epitopes for vaccine purposes. Some of these epitopes could act as immunogenic agents alone, offering several advantages, such as low cost and ease of handling without special treatment. Introduction to molecular biology of influenza a viruses The influenza virus: Antigenic composition and immune response Host adaptation and transmission of influenza A viruses in mammals Seasonal, avian, and novel H1N1 influenza: prevention and treatment modalities Seasonal influenza vaccine dose distribution in 157 countries Evolution and ecology of influenza A viruses Development of high-yield influenza A virus vaccine viruses Prospects of HA-based universal influenza vaccine Advances in universal influenza virus vaccine design and antibody mediated therapies based on conserved regions of the hemagglutinin Broadly neutralizing antibodies against influenza viruses Vaccine options for influenza: thinking small Neuraminidase as an influenza vaccine antigen: a low hanging fruit, ready for picking to improve vaccine effectiveness Structural basis of preexisting immunity to the 2009 H1N1 pandemic influenza virus Structure of the haemagglutinin membrane glycoprotein of influenza virus at 3 A resolution A perspective on the structural and functional constraints for immune evasion: insights from influenza virus Structural insights into the design of novel anti-influenza therapies NAction! How can neuraminidase-based immunity contribute to better influenza virus vaccines? MBio Influenza neuraminidase as a vaccine antigen Theoretical analysis of the neuraminidase epitope of the Mexican A H1N1 influenza strain, and experimental studies on its interaction with rabbit and human hosts Contribution of antibody production against neuraminidase to the protection afforded by influenza vaccines Influenza A: Biología, vacunas, y origen del virus pandémico A/H1N1 Design of an epitope-based peptide vaccine against spike protein of human coronavirus: an in silico approach Peptides for immunological purposes: design, strategies and applications The molecular bases of delta/alphabeta T cellmediated antigen recognition Identification of B cell and T cell epitopes using synthetic peptide combinatorial libraries The influenza virus resource at the National Center for Biotechnology Information MUSCLE: multiple sequence alignment with high accuracy and high throughput Clustal W and Clustal X version 2.0 SWISS-MODEL: modelling protein tertiary and quaternary structure using evolutionary information The SWISS-MODEL workspace: a web-based environment for protein structure homology modelling The SWISS-MODEL repository and associated resources Automated comparative protein structure modeling with SWISS-MODEL and Swiss-PdbViewer: a historical perspective PDBsum additions Knowledgebased validation of protein-structure coordinates derived by X-ray crystallography and nmr-spectroscopy VMD: visual molecular dynamics A continuous peptide epitope reacting with pandemic influenza AH1N1 predicted by bioinformatic approaches MHC class II epitope predictive algorithms Reliable prediction of T-cell epitopes using neural networks with novel sequence representations NN-align An artificial neural network-based alignment algorithm for MHC class II peptide binding prediction Scalable molecular dynamics with NAMD All-atom empirical potential for molecular modeling and dynamics studies of proteins Constant-pressure molecular-dynamics algorithms Software news and updates. Carma: a molecular dynamics analysis program Neutrophils extracellular traps damage Naegleria fowleri trophozoites opsonized with human IgG Influenza viruses: an introduction Simplified plaque reduction neutralization assay for dengue viruses by semimicro methods in BHK-21 cells: comparison of the BHK suspension test with standard plaque reduction neutralization Differences in antibody responses of individuals with natural infection and those vaccinated against pandemic H1N1 2009 influenza Immuno-fluorescence assay of leptospiral surface-exposed proteins Neutralization of virus infectivity by antibodies: old problems in new perspectives Structural and thermodynamic approach to peptide immunogenicity Virus-specific T cells as correlate of (cross-)protective immunity against influenza The effector T cell response to influenza infection Antiviral CD8+ T cell effector activities in situ are regulated by target cell type B cell responses to influenza infection and vaccination T cell activation Exploring the immunome: a brave new world for human vaccine development Bronchus-and nasal-associated lymphoid tissues Mechanism and regulation of class switch recombination Protection against influenza virus infection in polymeric Ig receptor knockout mice immunized intranasally with adjuvant-combined vaccines A nasally administered trivalent inactivated influenza vaccine is well tolerated, stimulates both mucosal and systemic immunity, and potentially protects against influenza illness Mucosal IgA responses in influenza virus infections; thoughts for vaccine design How a single amino acid change may alter the immunological information of a peptide Heterosubtypic neutralizing monoclonal antibodies cross-protective against H5N1 and H1N1 recovered from human IgM+ memory B cells Immunity to pre-1950 H1N1 influenza viruses confers cross-protection against the pandemic swineorigin 2009 A (H1N1) influenza virus Vaccination of monoglycosylated hemagglutinin induces cross-strain protection against influenza virus infections Preexisting human antibodies neutralize recently emerged H7N9 influenza strains Pre-existing immunity against swine-origin H1N1 influenza viruses in the general human population Assessment of seasonal influenza A virus-specific CD4 pandemic H1N1 swine-origin influenza A virus Seasonal influenza can poise hosts for CD4 T-cell immunity to H7N9 Avian influenza Immunoinformatic comparison of T-cell epitopes contained in novel swine-origin influenza A (H1N1) virus with epitopes in 2008-2009 conventional influenza vaccine DENV-2 subunit proteins fused to CR2 receptorbinding domain (P28)-induces specific and neutralizing antibodies to the Dengue virus in mice Hemagglutination inhibiting antibodies and protection against seasonal and pandemic influenza infection Publisher's Note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations Rafael Campos-Rodríguez citlicampos@gmail.com