Large-Scale Analysis of B-Cell Epitopes on Influenza Virus Hemagglutinin – Implications for Cross-Reactivity of Neutralizing Antibodies

Influenza viruses continue to cause substantial morbidity and mortality worldwide. Fast gene mutation on surface proteins of influenza virus result in increasing resistance to current vaccines and available antiviral drugs. Broadly neutralizing antibodies (bnAbs) represent targets for prophylactic and therapeutic treatments of influenza. We performed a systematic bioinformatics study of cross-reactivity of neutralizing antibodies (nAbs) against influenza virus surface glycoprotein hemagglutinin (HA). This study utilized the available crystal structures of HA complexed with the antibodies for the analysis of tens of thousands of HA sequences. The detailed description of B-cell epitopes, measurement of epitope area similarity among different strains, and estimation of antibody neutralizing coverage provide insights into cross-reactivity status of existing nAbs against influenza virus. We have developed a method to assess the likely cross-reactivity potential of bnAbs for influenza strains, either newly emerged or existing. Our method catalogs influenza strains by a new concept named discontinuous peptide, and then provide assessment of cross-reactivity. Potentially cross-reactive strains are those that share 100% identity with experimentally verified neutralized strains. By cataloging influenza strains and their B-cell epitopes for known bnAbs, our method provides guidance for selection of representative strains for further experimental design. The knowledge of sequences, their B-cell epitopes, and differences between historical influenza strains, we enhance our preparedness and the ability to respond to the emerging pandemic threats.

Introduction

Influenza epidemics result in substantial morbidity and mortality (1). The World Health Organization (WHO) Global Influenza Network provides annual recommendations on antigenic variants to be included in the influenza vaccine formulations. Influenza virus has low-fidelity polymerases that result in high mutation rates (2). As a consequence, seasonal influenza viruses efficiently escape from acquired immunity in the human population through antigenic drift increasing the impact of seasonal influenza. The antigenic shift in influenza A viruses – the reassortment of multiple viral genomes resulting in new strains with recombined antigens – leads to occasional worldwide pandemics that result in significant morbidity and, usually, high mortality. High transmissibility of influenza combined with rapid mutation rates makes the discovery of novel influenza therapeutics an imperative (3). The main challenge in developing antibody-based prophylactics and therapeutic vaccine against influenza is to understand the variation generated by the virus and developing means to elicit broadly neutralizing antibody responses.

The majority of neutralizing antibodies (nAbs) generated during a normal immune response target hemagglutinin (HA) and block viral entry into host cells (4). However, significant sequence diversity among HA genes limits the protective breadth of these nAbs (5). This sequence diversity of influenza A virus is high – there are 17 HA serotypes that belong into two major groups called group 1 (Grp1: H1, H2, H5, H6, H8, H9, H11, H12, H13, H16, and H17), and group 2 (Grp2: H3, H4, H7, H10, H14, and H15) (6). C179, the first neutralizing antibody reported to neutralize strains from H1 and H2 of influenza A virus, was isolated from mice immunized with the A/Okuda/57 (H2N2) strain (7). Later it was found that C179 was able to cross-neutralize H1, H2, H5, H6, and H9 subtypes (8–11). The next major advance in the field came about 15 years later (12), a novel class of human antibodies encoded by the V_H1–69 gene were discovered. Among these antibodies, a series of broadly neutralizing antibodies (bnAbs) have been described, such as CR6261 and F10 (13). Most bnAbs that neutralize influenza A virus have been reported to neutralize strains from either exclusively Grp1 or Grp2. FI6v3 (14) and 39.29 (5) are the only antibodies reported to neutralize human influenza isolates from both Grp1 and Grp2. Influenza B viruses are classified within a single influenza type, with two antigenically and genetically distinct lineages that co-circulate (15), represented by the prototype viruses B/Victoria/2/1987 (Victoria lineage) and B/Yamagata/16/1988 (Yamagata lineage) (16). Antibody CR8071 (17) is a bnAb against influenza B viruses, with neutralizing ability for both Victoria and Yamagata lineages. bnAb CR9114 (17) binds a conserved epitope on the HA stem and was shown to neutralize all tested influenza A viruses. However, it did not show in vitro neutralizing activity against influenza B viruses at the tested concentrations (17).

Generally, the neutralizing effectiveness of these bnAbs was evaluated using representative strains from the subtypes of influenza A virus or lineages of influenza B virus. Because of the high variability of HA genes, such evaluation might result in a conclusion that is limited to the tested viral variants. To determine the landscape of nAbs and better understand their cross-reactivity properties, we performed a systematic study of B-cell epitopes of a selection of nAbs against influenza virus. Antibodies recognize discrete sites on the surface of macromolecule called B-cell epitopes (antigenic determinants). Some 10% of B-cell epitopes are linear peptides while 90% are formed from discontinuous amino acids that create surface patches through the three dimensional (3D) conformation of proteins (18). We defined a novel way of describing discontinuous motifs, using virtual peptides, to represent B-cell epitopes and further used this representation to estimate potential cross-reactivity and neutralizing coverage of these nAbs.

Functional characterization of the increasing number of nAbs and known crystal structures of these nAbs complexed with HA proteins enables us to precisely define their B-cell epitopes. A large number of sequences of influenza variants are available in public databases (19) enabling systematic bioinformatics analysis of cross-reactivity of nAbs against influenza virus. Such systematic analysis improves our understanding of antibody/antigen interactions, facilitates mapping of the known universe of target antigens, and allows the prediction of cross-reactivity. These methods and tools are useful for the design of broadly protective vaccines against emerging pathogens. This article describes a study of influenza HA cross-reactivity, but the method is applicable to any viral pathogen where information about nAbs and a collection of variant sequences of the target antigen are available.

Materials and Methods

Neutralizing Antibodies Against Hemagglutinin

The names and specificities of nAb against influenza virus HA were collected from published papers. Twenty-two nAbs against influenza virus with crystal structures available in PDB were collected from published articles (Table 1). Fifteen of these nAbs target at the globular head of HA, and for the other seven, the binding sites are located on HA stem region.

TABLE 1

Table 1. Summary of well-characterized neutralizing antibodies against influenza virus.

FIGURE 1

Figure 1. Phylogenetic tree of 17 influenza A virus subtypes. Representative sequences were selected for each subtype (34) and the phylogenetic tree was made with ClustalX (35) and TreeView (36). H1: A/California/04/2009(H1N1); H2: A/Singapore/1/1957(H2N2); H3: A/Aichi/2/1968(H3N2); H4: A/duck/Czechoslovakia/1956(H4N6); H5: A/VietNam/1203/2004(H5N1); H6: A/chicken/California/431/2000(H6N2); H7: A/Turkey/Italy/8458/2002(H7N3); H8: A/Turkey/Ontario/6118/1968(H8N4); H9: A/Swine/Hong Kong/9/98(H9N2); H10: A/chicken/Germany/N/1949(H10N7); H11: A/duck/England/1/1956(H11N6); H12: A/duck/Alberta/60/1976(H12N5); H13: A/gull/Maryland/704/1977(H13N6); H14: A/Mallard/Astrakhan/263/1982(H14N5); H15: A/shearwater/West Australia/2576/79(H15N9); H16: A/black-headed gull/Sweden/2/99(H16N3); and H17: A/little yellow-shouldered bat/Guatemala/060/2010 (H17N10).

The majority of these nAbs were observed to bind or neutralize influenza A virus isolated either from Grp1 or Grp2. Antibodies FI6v3, CR9114, and 39.29 were shown to neutralize influenza strains within both Grp1 and Grp2 (5, 14, 27). Antibodies CR8059 and CR8071 (17) were the only two nAbs for influenza B virus. CR8059 is a light chain D95aN variant of CR8071. Since the mutation on CR8059 is not present at the binding interface and does not affect the binding, only CR8071 was used in the following study (17). The majority of these nAbs were shown to neutralize more than one strain, some of them are broadly neutralizing across subtypes of influenza A virus or lineages of influenza B virus. The Abs BH151, HC19, HC45, and HC63 were shown to specifically neutralize HA from the A/X-31(H3N2) strain. The available structures of nAb/HA complexes were downloaded from PDB (37).

Validated Influenza Strains by Neutralizing Antibodies

Binding and neutralization assays were collected from published materials. Binding and non-binding strains were classified according to their affinity measurements. The thresholds used to discriminate binding and non-binding strains were inconsistent in different studies: the lowest affinity detectable values were set as 10⁻⁴ M (17), 10⁻⁵ M (33), and ~10⁻⁶ M (20). In some reports, nAbs showed positive binding results but did not display neutralization ability to the same strains [e.g., nAb CR9114 against strain B/Florida/4/2006 (Yamagata) (17)]. Because of the lack of standardized thresholds and ambiguous definition of binding, only results that indicate non-binding of antibodies were considered as useful information and were retained for the subsequent analysis as negatives.

The neutralized and the escape strains were detected using the microneutralization assay (38) or HA inhibition assay (39). Several measurements were suggested in these studies:

1. The lowest concentration of nAb that displayed inhibition of hemagglutination or microneutralizing activity were set as either 2.5 μg/mL (40) or 5 μg/mL (41).

2. The 50% inhibitory concentration was set to IC₅₀ = 50 μg/mL (17).

3. The effective concentration of antibody needed to inhibit at least 99% of viral infectivity was set as EC₉₉ = 100 μg/mL (24, 25).

The HA sequences of strains that were experimentally validated for neutralization by studied antibodies (“validated strains”) were retrieved from the literature. The influenza strains HA sequences were collected from the literature or, if absent, from the Influenza Knowledge Base (FLUKB)¹. All experimentally validated strains were grouped into either neutralized strains or escape strains. The neutralized strains were selected based on reported experimental evidence. The escape strains included true escape strains as well as strains that were reported not to bind nAbs. We did not find any discrepancies in reported neutralizing properties across different studies used to collect functional data.

Hemagglutinin Sequences

All HA sequences were downloaded from the Influenza Knowledge Base (FLUKB¹, dated August 26th, 2013). After removing the incomplete sequences (fragments), 45,812 full-length HA sequences were left in the data set (HA sequence dataset) for further analysis.

Generation of Multiple Sequence Alignment of Hemagglutinin Sequences

The HA sequences of influenza strains from FLUKB were aligned using the MAFFT tool (42). The resulting multiple sequence alignment (MSA) results provided a consistent numbering scheme for all the further analyses. MSA were generated for both experimentally validated strains of HA and for all entries from FLUKB. For each nAb, every HA sequence from the crystal structure and from the experimentally validated strains were searched individually within the FLUKB database to find a strain with highest similarity using BLAST (43). This procedure was done to ensure that residue position mapping in following steps is consistent with the numbering scheme.

Identification of B-Cell Epitopes

B-cell epitope were identified from antigen–antibody structure, using a formula with the combination of the measurements of accessible surface area (ASA) and atom distance. For each residue from HA antigen, the ASA value was calculated using Naccess software (44) for both free HA and for HA coupled with an antibody. Residues r_i with ASA loss more than 20% were selected as epitope residues,

$r_i\in \left\{\mathrm{epitope\hspace{0.17em}residues}\right\}\mathrm{if}\frac{{\mathrm{ASA}}_{\mathrm{free}}-{\mathrm{ASA}}_{\mathrm{coupled}}}{{\mathrm{ASA}}_{\mathrm{free}}}>0.2.$

The majority of contacts between two contacting atoms occur at distance smaller than 5 Å separation (45). Euclidean distance was calculated between atoms a_i and a_j using their coordinates a_i(x_i, y_i, z_i) and a_j(x_j, y_j, z_j) in PDB structure data,

$d_{ij}=\sqrt{{\left(x_i-x_j\right)}^2+{\left(y_i-y_j\right)}^2+{\left(z_i-z_j\right)}^2}.$

Hemagglutinin residues r_i whose minimum atom distance to the closest nAb atom was within 4 Å were also incorporated in the epitope. The minimal atom distance was defined as:

$\begin{array}{llll}\hfill & d_{\min }=\min \left\{d_{ij}\right\},a_i\in \mathrm{antigen\hspace{0.17em}residue}{r}_i,a_j\in \mathrm{antibody\hspace{0.17em}residue}{r}_j,& \hfill & \\ \hfill & r_i\in \left\{\mathrm{epitope\hspace{0.17em}residue}\right\}\mathrm{if}{d}_{\min }<4\AA .& \hfill & \end{array}$

The residues that satisfy either of these two conditions (ASA loss or minimum distance) are considered to constitute a B-cell epitope.

The specific residues on HA that form hydrogen bonds, salt bridges, disulfide bonds, and covalent bonds between the HA and nAb were considered to define a B-cell epitope. The antigen/antibody interaction was further analyzed using PISA tool (46). The analysis of HA structures showed that all the hydrogen bonds, salt bridges, disulfide bonds, and covalent bonds between HA and nAb in each studied structure were incorporated in B-cell epitopes defined in the previous step.

Extraction of Discontinuous Motifs from Validated Strains

For each nAb, using the MSA result and the standardized numbering, the residue positions of B-cell epitope identified from the HA/antibody crystal structure were mapped onto all HA sequence of validated strains. Then discontinuous motifs composed of mapped residues were extracted from these sequences. These discontinuous motifs were classified as either “neutralized” or “escape” motifs according to the experimental validation status of the corresponding strain.

Mapping of Discontinuous Motifs to HA Sequence Dataset

For each nAb, based on the MSA result, the residue positions of B-cell epitope identified from the HA/antibody crystal structure were mapped onto the HA sequence dataset. A “discontinuous peptide” composed of amino acids that form B-cell epitope, in order that they appear in the sequence, was extracted from each HA sequence. By comparing the discontinuous peptides to all validated neutralized and escape motifs from experimentally validated strains, each discontinuous peptide was classified as neutralized (if 100% matching a neutralized epitope motif), escape (if 100% matching an escape epitope motif), or non-validated (if 100% matching validation data are missing). The term “discontinuous motif” indicates positions that define each B-cell epitope extracted from experimentally validated strains collected from publications, while term “discontinuous peptide” represents specific B-cell epitopes extracted from the HA sequence dataset.

Results

B-Cell Epitope Regions

For each nAb, the B-cell epitope was identified from the crystal structure as described in Section “Materials and Methods.” The structure of nAb F10-H5 (13) and identified epitope are illustrated in Figure 2. After B-cell epitopes of all studied nAbs were mapped to the same template structure, the overlapping of binding sites were found among different nAbs, particularly at the receptor-binding site (RBS), which is the necessary structure for binding to the sialic acid receptors during virus infection.

FIGURE 2

Figure 2. B-cell epitope on the structure of neutralizing antibody F10 binding HA protein (PDB ID: 3FKU). (A) Complex of F10-HA [A/Vietnam/1203/04(H5N1)]. The structure is a HA trimer of three identical copies (one of them is colored as cyan and green; the other two are in gray). Each copy contains the HA1 (cyan) and HA2 (green) chain, also the heavy chain of F10 (red), the neutralized epitope is highlighted in pink; (B) Close-up view of neutralized epitope identified on the structure (highlighted as pink surface).

For cross-reactive nAbs against influenza A virus, four major binding locations on HA structure are apparent: two of them reside on the globular head of HA and the other two target the stem region of HA (Table 2; Figure 3). The RBS is a heavily targeted area, with overlapping epitopes defined by eight nAbs. The only nAb that binds HA head but not the RBS is 2D1 (21). The 2D1 recognizes the Sa site of A/South Carolina/1/1918(H1N1). Sa site is one of the earliest known antigenic sites (47), which is proximal to the receptor-binding pocket. The detailed comparison of epitope residue positions between 2D1 and the other HA head-targeted nAbs are listed in Table 3. In contrast to the Abs that interact with the HA head, a series of nAbs recognize another highly conserved helical region in the membrane-proximal HA stem. The epitopes on F subdomain (CR6261, 39.29, etc.) and stem base (CR8020) are adjacent to each other, with a small number of shared residues. The only broadly nAb neutralizing influenza B virus, CR8071 binds to the lower region of the globular head of HA – the “head base” (Figure 3C). All the remaining antibodies analyzed in our study bind specifically the HA on A/X-31(H3N2) strain. All X-31 specific nAbs complex with the membrane-distal domain of HA. NAbs BH151 and HC45 (22) recognize a single epitope located at the base of the eight-stranded antiparallel β-sheet structure. The HC19 binding site is adjacent to the RBS. The HC63 epitope shares several residues with HC19, thereby the antibody binding site overlaps the membrane-distal domains of two HA monomers.

TABLE 2

Table 2. B-cell epitope regions of the 22 neutralizing antibodies.

FIGURE 3

Figure 3. The distinct B-cell epitope regions recognized by representative nAbs. The B-cell epitope regions of (A) represent cross-reactive nAbs against influenza A virus; (B) represent strain-specific nAbs against X-31(H3N2); (C) represent broadly nAb CR8071 against influenza B virus. The epitope regions of nAbs target influenza A virus were mapped on the monomer (A) or trimer (B) HA from A/X-31(H3N2) (PDB ID: 1KEN). The structure of B/Brisbane/60/2008 HA (PDB ID: 4FQM) was used as a template structure for influenza B virus. Different colors here were used for distinguishing B-cell epitope regions.

TABLE 3

Table 3. B-cell epitope overlap for nAbs targeting HA head region.

Experimentally Validated Discontinuous Motifs

Discontinuous motifs were extracted from the validated sequences as described in Section “Materials and Methods,” and presented by WebLogo (48) and BlockLogo² [Ref. (49)]. WebLogo figures consist of stacks of amino acids, while the overall height of the stack indicates the sequence conservation at that position, and the height of symbols within the stack indicates the relative frequency of each amino or nucleic acid at that position. While BlockLogo is a web-based application for visualization of protein and nucleotide fragments, continuous protein sequence motifs, and discontinuous sequence motifs using calculation of block entropy from MSAs. The BlockLogo figures present the actual combinations of amino acids, and the height of each combination represents its relative frequency. In the nAb F10 as an example, the neutralized and escape discontinuous motifs are shown in Figures 4A,C (WebLogo figures), and Figures 4B,D (BlockLogo figures). WebLogos show a clear overall description of each residue conservation difference between individual neutralized and escape motifs. For example, 44N, 48T, 304R/D, 380L/Y, 391N, 394E/A/L on F10 epitope region are likely to contribute to the escape strains. In the BlockLogo figures, specific neutralized and escape B-cell epitopes of F10 were listed with their frequencies, which can be used for their direct comparison.

FIGURE 4

Figure 4. Neutralized and escape discontinuous motifs from experimentally validated sequences with nAb F10. The WebLogo shows global (A) neutralized motifs, and (C) escape motifs and BlockLogo shows individual (B) neutralized motifs, and (D) escape motifs. The extracted discontinuous motif extracted from the structure (PDB ID: 3FKU, chain A and B), corresponds to the positions of reference sequence [FLU0293715, A/Viet Nam/1203/2004(H5N1)]: 24, 44, 46, 48, 304, 305, 306, 331, 360, 361, 362, 363, 380, 383, 384, 387, 391, 394, 395, and 398.

Analysis of Variation of Discontinuous Peptides in HA Sequences Dataset

For each nAb, the residue positions of their B-cell epitopes were mapped on the complete HA sequences dataset collected from the FLUKB. Amino acid strings representing discontinuous peptides were extracted from the HA sequence of each strain. The variability of discontinuous peptides and validated discontinuous motif coverage were analyzed for each nAb.

For example, for the nAb F10, 589 different patterns of discontinuous peptides were generated among all 45,812 sequences in HA sequence dataset, using the F10 B-cell epitope identified from the crystal structure. In the next step, the discontinuous peptides were sorted according to their frequencies. The second most frequent peptide in FLUKB is identical an escape motif, while the 6th, 8th, and 19th are each identical to one of the neutralized motifs. However, the most frequent F10 discontinuous peptide in FLUKB (see text footnote 1) has not been experimentally tested (Figure 5), along with other 14 discontinuous peptides. The analysis of differences between the most frequent discontinuous peptide and neutralized or escape motifs was inconclusive. Therefore future experimental studies should include a representative sequence containing the discontinuous peptide HHVLSLPTVDGWLTQITVNI that is present in more than 10,000 entries in the FLUKB. We also recommend that motifs 1, 4, 5, 7, 9–18, and 20 are considered for the experimental validation. The remaining sequences are less common, each having <400 sequences in the data set.

FIGURE 5

Figure 5. Frequencies of top 20 discontinuous peptides (B-cell epitope of nAb F10) from the HA sequence dataset. FLU0243751 (A/Viet Nam/1203/2004) was used as reference HA sequence in the analysis of F10 B-cell epitopes. The corresponding positions of discontinuous peptides on FLU0243751 are: 24, 44, 46, 48, 304, 305, 306, 331, 360, 361, 362, 363, 380, 383, 384, 387, 391, 394, 395, and 398. Discontinuous peptides that were identical to neutralized motifs are shown in blue, while those identical to escape motifs are shown in red. The sequences of Top 20 most frequent discontinuous peptides are listed along with their validation status.

The discontinuous peptides were generated and the variability was investigated for all cross-reactive nAbs (Table 4). The B-cell epitope regions on the HA stem are less variable as compared to the epitopes on the HA head. The specific result generated within each subtype in HA sequence dataset show similar patterns as for all subtypes (data not shown). This conclusion is consistent with our previous knowledge that the globular head of HA1 has a higher mutation rate than the stem (29), making the stem a more conserved region for bnAbs targeting.

TABLE 4

Table 4. The number of different discontinuous peptides from B-cell epitopes of each nAb in the HA sequence dataset.

Discontinuous Motifs Coverage in HA Sequences Dataset

The neutralized and escape discontinuous motifs of nAb F10 have covered 19 and 17% of FLUKB, respectively, while the discontinuous peptides from 64% of the strains have not been validated (Figure 6A). Viewed by the subtype, F10 neutralized coverage of subtypes H5, H8, H9, and H11 are higher (50–90%) than of H1 and H2 (5–20%), while the coverage of subtype H6 is negligible (8 in 1708 H6 sequences) (Figure 6B).

FIGURE 6

Figure 6. Motif coverage for nAb F10. (A) The coverage of neutralized and escape discontinuous motifs, and non-validated discontinuous peptides within 45,812 HA sequences extracted from the FLUKB; (B) The motif coverage by subtype, the numbers in brackets indicate the number of sequence within the specific subtype (among 45,812 HA sequences).

The motif coverage analysis within the 45,812 HA sequences was performed for all nAbs. For the nAbs with available cross-reactivity data, the motif coverages were different between the nAbs targeting the HA globular head and those targeting the stem part. The nAbs that bind stem normally have higher neutralized motif coverage than those that bind the globular head (Figure 7).

FIGURE 7

Figure 7. Discontinuous motif coverage in the HA sequence dataset for cross-reactive nAbs. Only the coverage data for cross-reactive nAbs are shown here. The nAbs are grouped based on their binding locations and influenza types, from left to right: Sa site, near RBS, F subdomain, stem base on influenza A virus, and head base on influenza B virus.

The motif coverage is shown as heat map for each subtype and each nAb (Figure 8). The nAbs (such as CR6261, CR9114, F10, and FI6v3) that target stem region are more cross-reactive – they cover more strains, and also more subtypes of influenza.

FIGURE 8

Figure 8. The neutralized and escape motif coverage on each subtype within FLUKB of neutralizing antibodies. The heat map of (A) neutralized and (B) escape coverage are shown. Each cell on the grid represents coverage for a specific antibody (row name) on specific subtype (column name). The subtypes were sorted from Grp1 and Grp2 in influenza A virus, and influenza B virus. Different color schemes were used in order to differentiate neutralizing and escape coverage: from green to red/green to blue indicate rising neutralizing/escape coverage in HA sequence dataset. The boxes with symbol “X” indicate that no experimental validation data were available for this study.

Combining of Neutralizing Antibodies

For each sequence in the HA sequence dataset, 22 strings (discontinuous peptides) were extracted to represent 22 B-cell epitopes by all nAbs analyzed in this study. The majority (82.62%) of all strains in FLUKB have at least one discontinuous peptide that is identical to the validated neutralized motifs (Table 5). A small number (2.25%) of sequences can be neutralized by as many as seven nAbs.

TABLE 5

Table 5. Distribution of the number of neutralizing antibodies that share identical neutralized discontinuous motif with sequences within the HA sequence dataset.

Here, we propose a combination of nAbs, where a small number of nAbs can cover a large proportion of influenza strains. The nAbs FI6v3, F10, CR9114, and CR8071 (Figure 9A) were selected, and the neutralized coverage has increased from 18.91% (F10), 4.06% (CR8071), 43.89% (CR9114), and 58.44% (FI6v3) to 78.45% (Figure 9B) when these antibodies were combined. These nAbs also covered most subtypes of influenza A virus and both lineages in influenza B virus (Figure 9C).

FIGURE 9

Figure 9. A combination of four neutralizing subtype-diversified and potent neutralizing antibodies. (A) The heat map of the neutralized result for four nAbs, the color scheme is same as Figure 8A; (B) the neutralized coverage of four nAbs individually, and the combination of nAbs on HA dataset; and (C) the neutralizing result of combination of four nAbs by subtype. The subtypes were sorted from Grp1 and Grp2 in influenza A virus and influenza B virus. Only subtypes with neutralizing data are shown.

Discussion

This study presents an overview of binding specificities of reported nAbs, as well as an estimate of their neutralization and escape coverage (neutralization effectiveness) in more than 45,000 HA sequences available in FLUKB. The variety and frequency of discontinuous peptides within different B-cell epitopes have been analyzed in the HA data set. The results of the analysis of discontinuous peptides provide insights into further experimental design: strains with peptides that have high frequency among the strain populations should be given priority for experimental validation and their neutralizing status for specific nAbs.

Of note, additional sequence changes in HA outside the nAb epitope may result in either local or quaternary structural alterations that impacts antibody binding to the epitope per se. Likewise, modification of glycosylation sites through sequence change may impact accessibility of antibodies to the neutralization site, creating discordance between sequence identity of binding site shown in BlockLogo and neutralization outcome between two strains of viruses sharing the same epitope sequence. The frequency of such occurrences will be important to determine. Neutralization assays of strains with discontinuous epitopes identical to validated B-cell epitopes will provide a proof of cross-neutralization. Since the experimental validation is time and money consuming, the introduction of extended B-cell epitope (see Supplementary Material) aims to help select representative sequences that differ in extended B-cell epitopes. For each proposed neutralizing or escape peptide (actual B-cell epitope), a small number of variants defined by changes in its environment (extended B-cell epitopes) constitute the majority of strains with the proposed peptide.

On the other hand, before more experimental data generated to fill the existing “non-validated gap,” it will be meaningful to bring out some reasonable estimation. The assumption and methods in this paper are based on complete identity to discontinuous motifs on B-cell epitope (additionally extended B-cell epitope). To check the validity of this assumption, the similarity between discontinuous motifs and discontinuous peptides could be used to estimate and predict neutralization and binding results in the future. For example, a discontinuous peptide with mutated residues of similar feature to the neutralized motif would be considered as “possible neutralized peptide” against specific nAbs. These estimations could also be validated in experimental assays, and then be used to further experimental design iteratively.

Conclusion

Over the past few years, our understanding of nAbs and their responses against influenza HA have expanded tremendously. Besides the well-known HA head region interactions, an increasing number of characterized nAbs bind and neutralize influenza virus by targeting the more conserved stem regions. Among these stem-targeting nAbs, some show broadly neutralizing ability across subtypes/lineages, even across two groups in influenza A virus strains. However, the related experimental data for majority of nAbs are quite limited.

In sum, we have established a library of validated motifs (extracted from HA sequences in neutralized and escape strains) for each nAb. For any newly emerging strain, the cross-neutralization prediction can be made rapidly for existing nAbs and validation experiments can be designed judiciously. This study provides a method for investigation of cross-reactivity of nAbs against influenza viruses, but is directly applicable to any viral pathogen that has structurally characterized nAbs and a collection of variant sequences of the target antigen. Examples of such pathogens include orthomyxoviruses (influenza); flaviviruses such as dengue or West Nile; arenaviruses such as lymphocytic choriomeningitis virus and human immunodeficiency virus, among others. Insights from such bioinformatics analyses coupled with antibody antigenicity through crystallographic determinations will facilitate electronic neutralization profiling that can be tested empirically in subsequent laboratory neutralization assays.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

Jing Sun, Vladimir Brusic, and Ellis L. Reinherz acknowledge funding from NIH grant U01 AI 90043. Ulrich J. Kudahl was funded by Oticon Foundation, Otto Mønsted Foundation, Julie Dam’s Stipend, Henry Shaws Stipend, and Reinholdt Jorcks Stipend. Christian Simon was funded by the Novo Scholarship Programme, Direktør Ib Henriksens Fond, and Augustinus Fonden.

Supplementary Material

The Supplementary Material for this article can be found online at http://www.frontiersin.org/Journal/10.3389/fimmu.2014.00038/abstract

Footnotes

1. ^http://research4.dfci.harvard.edu/cvc/flukb
2. ^http://research4.dfci.harvard.edu/cvc/blocklogo

References

1. Thompson WW, Shay DK, Weintraub E, Brammer L, Cox N, Anderson LJ, et al. Mortality associated with influenza and respiratory syncytial virus in the United States. JAMA (2003) 289(2):179–86. doi: 10.1001/jama.289.2.179

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

2. Parvin JD, Moscona A, Pan WT, Leider JM, Palese P. Measurement of the mutation rates of animal viruses: influenza A virus and poliovirus type 1. J Virol (1986) 59(2):377–83.

Pubmed Abstract | Pubmed Full Text

3. Fraser C, Donnelly CA, Cauchemez S, Hanage WP, Van Kerkhove MD, Hollingsworth TD, et al. Pandemic potential of a strain of influenza A (H1N1): early findings. Science (2009) 324(5934):1557–61. doi:10.1126/science.1176062

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

4. Wiley DC, Skehel JJ. The structure and function of the hemagglutinin membrane glycoprotein of influenza virus. Annu Rev Biochem (1987) 56:365–94. doi:10.1146/annurev.bi.56.070187.002053

CrossRef Full Text

5. Nakamura G, Chai N, Park S, Chiang N, Lin Z, Chiu H, et al. An in vivo human-plasmablast enrichment technique allows rapid identification of therapeutic influenza a antibodies. Cell Host Microbe (2013) 14(1):93–103. doi:10.1016/j.chom.2013.06.004

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

6. Steel J, Lowen AC, Wang TT, Yondola M, Gao Q, Haye K, et al. Influenza virus vaccine based on the conserved hemagglutinin stalk domain. MBio (2010) 1(1):e18–10. doi:10.1128/mBio.00018-10

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

7. Okuno Y, Isegawa Y, Sasao F, Ueda S. A common neutralizing epitope conserved between the hemagglutinins of influenza A virus H1 and H2 strains. J Virol (1993) 67(5):2552–8.

Pubmed Abstract | Pubmed Full Text

8. Yoshida R, Igarashi M, Ozaki H, Kishida N, Tomabechi D, Kida H, et al. Cross-protective potential of a novel monoclonal antibody directed against antigenic site B of the hemagglutinin of influenza A viruses. PLoS Pathog (2009) 5(3):e1000350. doi:10.1371/journal.ppat.1000350

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

9. Ueda M, Maeda A, Nakagawa N, Kase T, Kubota R, Takakura H, et al. Application of subtype-specific monoclonal antibodies for rapid detection and identification of influenza A and B viruses. J Clin Microbiol (1998) 36(2):340–4.

Pubmed Abstract | Pubmed Full Text

10. Smirnov YA, Lipatov AS, Gitelman AK, Claas EC, Osterhaus AD. Prevention and treatment of bronchopneumonia in mice caused by mouse-adapted variant of avian H5N2 influenza A virus using monoclonal antibody against conserved epitope in the HA stem region. Arch Virol (2000) 145(8):1733–41. doi:10.1007/s007050070088

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

11. Sakabe S, Iwatsuki-Horimoto K, Horimoto T, Nidom CA, Le M, Takano R, et al. A cross-reactive neutralizing monoclonal antibody protects mice from H5N1 and pandemic (H1N1) 2009 virus infection. Antiviral Res (2010) 88(3):249–55. doi:10.1016/j.antiviral.2010.09.007

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

12. Ekiert DC, Wilson IA. Broadly neutralizing antibodies against influenza virus and prospects for universal therapies. Curr Opin Virol (2012) 2(2):134–41. doi:10.1016/j.coviro.2012.02.005

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

13. Sui J, Hwang WC, Perez S, Wei G, Aird D, Chen LM, et al. Structural and functional bases for broad-spectrum neutralization of avian and human influenza A viruses. Nat Struct Mol Biol (2009) 16(3):265–73. doi:10.1038/nsmb.1566

CrossRef Full Text

14. Corti D, Voss J, Gamblin SJ, Codoni G, Macagno A, Jarrossay D, et al. A neutralizing antibody selected from plasma cells that binds to group 1 and group 2 influenza A hemagglutinins. Science (2011) 333(6044):850–6. doi:10.1126/science.1205669

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

15. Yamashita M, Krystal M, Fitch WM, Palese P. Influenza B virus evolution: co-circulating lineages and comparison of evolutionary pattern with those of influenza A and C viruses. Virology (1988) 163(1):112–22. doi:10.1016/0042-6822(88)90238-3

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

16. Rota PA, Wallis TR, Harmon MW, Rota JS, Kendal AP, Nerome K. Cocirculation of two distinct evolutionary lineages of influenza type B virus since 1983. Virology (1990) 175(1):59–68. doi:10.1016/0042-6822(90)90186-U

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

17. Dreyfus C, Laursen NS, Kwaks T, Zuijdgeest D, Khayat R, Ekiert DC, et al. Highly conserved protective epitopes on influenza B viruses. Science (2012) 337(6100):1343–8. doi:10.1126/science.1222908

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

18. Huang J, Honda W. CED: a conformational epitope database. BMC Immunol (2006) 7:7. doi:10.1186/1471-2172-7-7

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

19. Squires RB, Noronha J, Hunt V, Garcia-Sastre A, Macken C, Baumgarth N, et al. Influenza research database: an integrated bioinformatics resource for influenza research and surveillance. Influenza Other Respir Viruses (2012) 6(6):404–16. doi:10.1111/j.1750-2659.2011.00331.x

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

20. Tsibane T, Ekiert DC, Krause JC, Martinez O, Crowe JE Jr, Wilson IA, et al. Influenza human monoclonal antibody 1F1 interacts with three major antigenic sites and residues mediating human receptor specificity in H1N1 viruses. PLoS Pathog (2012) 8(12):e1003067. doi:10.1371/journal.ppat.1003067

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

21. Xu R, Ekiert DC, Krause JC, Hai R, Crowe JE Jr, Wilson IA. Structural basis of preexisting immunity to the 2009 H1N1 pandemic influenza virus. Science (2010) 328(5976):357–60. doi:10.1126/science.1186430

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

22. Fleury D, Daniels RS, Skehel JJ, Knossow M, Bizebard T. Structural evidence for recognition of a single epitope by two distinct antibodies. Proteins (2000) 40(4):572–8. doi:10.1002/1097-0134(20000901)40:4<572::AID-PROT30>3.3.CO;2-E

CrossRef Full Text

23. Ekiert DC, Kashyap AK, Steel J, Rubrum A, Bhabha G, Khayat R, et al. Cross-neutralization of influenza A viruses mediated by a single antibody loop. Nature (2012) 489(7417):526–32. doi:10.1038/nature11414

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

24. Whittle JR, Zhang R, Khurana S, King LR, Manischewitz J, Golding H, et al. Broadly neutralizing human antibody that recognizes the receptor-binding pocket of influenza virus hemagglutinin. Proc Natl Acad Sci U S A (2011) 108(34):14216–21. doi:10.1073/pnas.1111497108

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

25. Schmidt AG, Xu H, Khan AR, O’Donnell T, Khurana S, King LR, et al. Preconfiguration of the antigen-binding site during affinity maturation of a broadly neutralizing influenza virus antibody. Proc Natl Acad Sci U S A (2013) 110(1):264–9. doi:10.1073/pnas.1218256109

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

26. Bizebard T, Gigant B, Rigolet P, Rasmussen B, Diat O, Bosecke P, et al. Structure of influenza virus haemagglutinin complexed with a neutralizing antibody. Nature (1995) 376(6535):92–4. doi:10.1038/376092a0

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

27. Fleury D, Barrere B, Bizebard T, Daniels RS, Skehel JJ, Knossow M. A complex of influenza hemagglutinin with a neutralizing antibody that binds outside the virus receptor binding site. Nat Struct Biol (1999) 6(6):530–4. doi:10.1038/9299

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

28. Barbey-Martin C, Gigant B, Bizebard T, Calder LJ, Wharton SA, Skehel JJ, et al. An antibody that prevents the hemagglutinin low pH fusogenic transition. Virology (2002) 294(1):70–4. doi:10.1006/viro.2001.1320

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

29. Lee PS, Yoshida R, Ekiert DC, Sakai N, Suzuki Y, Takada A, et al. Heterosubtypic antibody recognition of the influenza virus hemagglutinin receptor binding site enhanced by avidity. Proc Natl Acad Sci U S A (2012) 109(42):17040–5. doi:10.1073/pnas.1212371109

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

30. Dreyfus C, Ekiert DC, Wilson IA. Structure of a classical broadly neutralizing stem antibody in complex with a pandemic H2 influenza virus hemagglutinin. J Virol (2013) 87(12):7149–54. doi:10.1128/JVI.02975-12

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

31. Throsby M, van den Brink E, Jongeneelen M, Poon LL, Alard P, Cornelissen L, et al. Heterosubtypic neutralizing monoclonal antibodies cross-protective against H5N1 and H1N1 recovered from human IgM+ memory B cells. PLoS One (2008) 3(12):e3942. doi:10.1371/journal.pone.0003942

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

32. Ekiert DC, Bhabha G, Elsliger MA, Friesen RH, Jongeneelen M, Throsby M, et al. Antibody recognition of a highly conserved influenza virus epitope. Science (2009) 324(5924):246–51. doi:10.1126/science.1171491

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

33. Ekiert DC, Friesen RH, Bhabha G, Kwaks T, Jongeneelen M, Yu W, et al. A highly conserved neutralizing epitope on group 2 influenza A viruses. Science (2011) 333(6044):843–50. doi:10.1126/science.1204839

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

34. Sun X, Shi Y, Lu X, He J, Gao F, Yan J, et al. Bat-derived influenza hemagglutinin H17 does not bind canonical avian or human receptors and most likely uses a unique entry mechanism. Cell Rep (2013) 3(3):769–78. doi:10.1016/j.celrep.2013.01.025

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

35. Larkin MA, Blackshields G, Brown NP, Chenna R, McGettigan PA, McWilliam H, et al. Clustal W and Clustal X version 2.0. Bioinformatics (2007) 23(21):2947–8. doi:10.1093/bioinformatics/btm404

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

36. Page RD. TreeView: an application to display phylogenetic trees on personal computers. Comput Appl Biosci (1996) 12(4):357–8.

37. Berman HM, Westbrook J, Feng Z, Gilliland G, Bhat TN, Weissig H, et al. The protein data bank. Nucleic Acids Res (2000) 28(1):235–42. doi:10.1093/nar/28.1.235

CrossRef Full Text

38. Katz J, Hancock K, Veguilla V, Zhong W, Lu X, Sun H, et al. Serum cross-reactive antibody response to a novel influenza A (H1N1) virus after vaccination with seasonal influenza vaccine. MMWR Morb Mortal Wkly Rep (2009) 58(19):521–4.

Pubmed Abstract | Pubmed Full Text

39. Stavitsky AB. Micromethods for the study of proteins and antibodies II. Specific applications of hemagglutination and hemagglutination-inhibition reactions with tannic acid and protein-treated red blood cells. J Immunol (1954) 72(5):368–75.

40. Yu X, Tsibane T, McGraw PA, House FS, Keefer CJ, Hicar MD, et al. Neutralizing antibodies derived from the B cells of 1918 influenza pandemic survivors. Nature (2008) 455(7212):532–6. doi:10.1038/nature07231

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

41. Krause JC, Tumpey TM, Huffman CJ, McGraw PA, Pearce MB, Tsibane T, et al. Naturally occurring human monoclonal antibodies neutralize both 1918 and 2009 pandemic influenza A (H1N1) viruses. J Virol (2010) 84(6):3127–30. doi:10.1128/JVI.02184-09

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

42. Katoh K, Toh H. Recent developments in the MAFFT multiple sequence alignment program. Brief Bioinform (2008) 9(4):286–98. doi:10.1093/bib/bbn013

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

43. Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, et al. Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res (1997) 25(17):3389–402. doi:10.1093/nar/25.17.3389

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

44. Hubbard SJ, Thornton JM. Naccess. Computer Program. (Vol. 2). London: Department of Biochemistry and Molecular Biology, University College London (1993).

45. McConkey BJ, Sobolev V, Edelman M. Discrimination of native protein structures using atom-atom contact scoring. Proc Natl Acad Sci U S A (2003) 100(6):3215–20. doi:10.1073/pnas.0535768100

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

46. Krissinel E, Henrick K. Inference of macromolecular assemblies from crystalline state. J Mol Biol (2007) 372(3):774–97. doi:10.1016/j.jmb.2007.05.022

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

47. Wilson IA, Skehel JJ, Wiley DC. Structure of the haemagglutinin membrane glycoprotein of influenza virus at 3 A resolution. Nature (1981) 289(5796):366–73. doi:10.1038/289366a0

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

48. Crooks GE, Hon G, Chandonia JM, Brenner SE. WebLogo: a sequence logo generator. Genome Res (2004) 14(6):1188–90. doi:10.1101/gr.849004

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

49. Olsen LR, Kudahl UJ, Simon C, Sun J, Schönbach C, Reinherz EL, et al. BlockLogo: Visualization of peptide and sequence motif conservation. J Immunol Methods (2013) 400:37–44. doi:10.1016/j.jim.2013.08.014

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text