- Open Access
HASP server: a database and structural visualization platform for comparative models of influenza A hemagglutinin proteins
BMC Bioinformatics volume 14, Article number: 197 (2013)
Influenza A viruses possess RNA genomes that mutate frequently in response to immune pressures. The mutations in the hemagglutinin genes are particularly significant, as the hemagglutinin proteins mediate attachment and fusion to host cells, thereby influencing viral pathogenicity and species specificity. Large-scale influenza A genome sequencing efforts have been ongoing to understand past epidemics and pandemics and anticipate future outbreaks. Sequencing efforts thus far have generated nearly 9,000 distinct hemagglutinin amino acid sequences.
Comparative models for all publicly available influenza A hemagglutinin protein sequences (8,769 to date) were generated using the Rosetta modeling suite. The C-alpha root mean square deviations between a randomly chosen test set of models and their crystallographic templates were less than 2 Å, suggesting that the modeling protocols yielded high-quality results. The models were compiled into an online resource, the Hemagglutinin Structure Prediction (HASP) server. The HASP server was designed as a scientific tool for researchers to visualize hemagglutinin protein sequences of interest in a three-dimensional context. With a built-in molecular viewer, hemagglutinin models can be compared side-by-side and navigated by a corresponding sequence alignment. The models and alignments can be downloaded for offline use and further analysis.
The modeling protocols used in the HASP server scale well for large amounts of sequences and will keep pace with expanded sequencing efforts. The conservative approach to modeling and the intuitive search and visualization interfaces allow researchers to quickly analyze hemagglutinin sequences of interest in the context of the most highly related experimental structures, and allow them to directly compare hemagglutinin sequences to each other simultaneously in their two- and three-dimensional contexts. The models and methodology have shown utility in current research efforts and the ongoing aim of the HASP server is to continue to accelerate influenza A research and have a positive impact on global public health.
Influenza A viruses (IAV) are among the most common causes of human respiratory infections and among the most significant because they cause high morbidity and mortality, both in annual epidemics and in unpredictable pandemics . IAVs are enveloped negative-strand RNA viruses with segmented genomes containing 8 gene segments encoding at least 11 open reading frames . IAVs are covered with proteins, notably hemagglutinin (HA) and neuraminidase (NA). The combination of alleles of HA and NA define strain nomenclature because of their variability and importance in pathogenicity, in host immunity, and in host species specificity. Currently, there are 17 subtypes of HA and 10 subtypes of NA.
HA is a glycosylated type I integral membrane protein that functions both as the viral receptor-binding protein and fusion protein. HA recognizes sialic-acid (SA) bound glycans with variable specificity for SA with α2-3 or α2-6 glycosidic linkages; these linkages are critical in determining host species specificity . The HA protein is a trimer, with each monomer composed of a heavy (~40 kDa) and light (~20 kDa) chain cleaved from a single precursor. Given the HA diversity of IAVs, structural predictions are an important tool to map receptor-binding and antigenic regions.
To date, nearly 9,000 distinct HA alleles have been sequenced , which is roughly three orders of magnitude larger than the number of experimentally determined structures (Table 1). Given the large scale of HA sequencing efforts and the experimental barriers in determining HA structures, the gap between the experimental characterization of primary and tertiary structures is likely to persist. Despite the many changes HAs accumulate in their primary amino acid sequences as a result of antigenic shift and drift, the tertiary structures of even highly divergent HAs are largely conserved. For example, the HA structures of a group 1 H9 subtype (PDB ID: 1JSD ) and a group 2 H7 subtype (PDB ID: 1TI8 ) can be superposed with a C-alpha Root Mean Square Deviation (RMSD)of 2.8 Å over 467 amino acids, despite only 40% sequence identity. With this structural similarity between divergent HAs, comparative models of HA proteins may be of reasonable quality and hence useful in filling the information gap between HA sequence and structure space. Here we present the Hemagglutinin Structure Prediction (HASP) server, a database of HA comparative models of all publicly available HA sequences. The database is provided with an integrated molecular viewer to accelerate model analysis and hypothesis generation.
Construction and content
All crystal structures of influenza HAs were retrieved from the Protein Data Bank (PDB) and filtered through the PISCES server  using a 95% sequence identity cutoff. Additional criteria included having an R-factor below 0.30 and a resolution better than 4.0 Å, resulting in a set of 12 structures (Table 1). Each of these structures was processed with the Molprobity server  to optimize hydrogen positions and rotamers for asparagine, glutamine, and histidine residues. A symmetric trimer for each model was then generated from the first (or only) monomer in the model by applying the crystallographic or non-crystallographic symmetry operators using the symmetry functionality of ROSETTA3, v3.2 . The resulting models and their amino acid sequences were used for comparative modeling.
All unique HA sequences were downloaded from the Influenza Research Database (8,769 as of June 2011) . The template sequences were searched with each sequence using BLAST . The top scoring high-scoring segment pair (hsp) was used as the alignment for comparative modeling.
A conservative approach was taken in the modeling, where only amino acid types, side-chain conformations, and dihedrals were allowed to change, with insertions and deletions omitted from the model. This approach was chosen because of the high sequence similarity between the query and template HA sequences and the desire to retain as much information from the crystal structure as possible.
The models were created by mapping the query sequence onto the template sequence from the hsp alignment using the fixed-backbone design functionality of ROSETTA3, v3.2 . Rotamers with χ1 angles ±1 standard deviation from canonical Dunbrack rotamers were included in modeling. Energy minimization of the side chain dihedral angles was also performed. For positions with identities between the query and template sequences in the hsp alignments, the side-chain conformation of the template structure was retained. Only aligned residues in the hsp alignment are present in the final models.
Model quality estimates
To derive semi-quantitative estimates of the quality that can be expected for modeling a given HA sequence, the sequences of the twelve HA crystal structures used as templates were modeled using the server, omitting their structures from the possible templates, and the models were compared to the crystal structures by various criteria (Figure 1, Table 2). The recovery rate of side chain conformations between the HASP-generated model and the crystal structure were calculated using the ROSETTA3 scientific benchmarking functions (Leaver-Fay A, O’Meara MJ, et al.; unpublished). The C-alpha RMSD between the template structure and the crystal structure were calculated. The rotamer recovery rate and C-alpha RMSD were evaluated as a function of sequence identity between the query sequences and sequences of the templates, the only information a user of the HASP server would have a priori. These relations could be fit by linear regression such that the rotamer recovery rate can be approximated as 0.32 * identities + 0.41, with an R2 of 0.88, and the C-alpha RMSD can be approximated as −2.52 * identities + 3.16, with an R2 of 0.62 (Figure 1). Given the limited dataset (n = 12), these equations are intended only to serve as an estimate of the quality of a model that can be expected for an HA query-template pair with a given sequence identity.
On a finer scale, the rotamer recovery rate was evaluated as a function of amino acid type (Table 2). Rotamer recovery rates were roughly proportional to hydrophobicity and the number of dihedrals for a given amino acid type. Large hydrophobic amino acids with fewer dihedrals (F, W, Y) were recovered at rates of approximately >90%. Amino acids with many dihedrals (K, R, M) were recovered at the lowest rates, 35-37%, followed by hydrophilic amino acids with fewer dihedrals (D, E, N, Q) with rates between 40-50%.
Utility and discussion
Main features of the HASP server
The HASP server interface has two primary components, the Search tab (buttons or tabs are identified hereafter in italics) to identify HA sequences of interest and the Viewer tab to display both the sequences and structures of the selected HA proteins. In the Search tab, HA sequences of interest can be selected from the database based on the H/N subtype of the strain, the geographical location of strain collection, by keyword, and other strain features.
To start, a user chooses the parameters of the query in the Search tab. By clicking Toggle Map Viewer, results can be narrowed by geographical location of strain collection through a color-coded world map; countries labeled in dark and light shades of green have the highest and lowest number of cases, respectively. Clicking Go updates the results, and after toggling out of the map viewer, a list of HA sequences is displayed. Information provided for each sequence includes the EMBL ID, Subtype, Strain, Year, Location (City, State/Country), and Species (host).
Sequences of interest are selected by checking the box in the right-hand column of the search results. In the Viewer tab, up to two of these models can be displayed simultaneously within the interface using the built-in molecular viewer. Residues can be easily identified within the model through interactive alignment functionality. In the viewer, preset views of the models may be selected through buttons (Hydrogens On, All Atoms On, Single Chain) with options for displaying models as cartoons, or with or without labels. For further analysis, a drop-down menu allows the user to export sequences or structures in FASTA or PDB format, respectively.
Discussion of HASP server use through a case study
In the World Health Organization’s Human Animal Interface online database for H5N1 Avian Influenza in Humans, it is reported that there were three deaths from H5N1 in Vietnam in 2003 and twenty deaths in 2004 . The differences between a pair of HA proteins sequenced from infected humans and chickens in Vietnam in 2003 and 2004 are explored below with the HASP server as a case study.
To narrow the search to HA sequences from H5N1 strains, H5 and N1 are selected from the subtype menu below the Search box (Figure 2A). To refine the search to H5N1 strains by country, Vietnam is selected on the map. The search term, “2003 OR 2004”, is added to the search box to restrict the search to strains sampled in those years (Figure 2B). Clicking the Go button updates the query, and toggling out of the map viewer displays the results page, where twenty-four H5 sequences are found for the H5N1 subtype in Vietnam between 2003-2004 (Figure 3). To compare the HA sequences sequenced from human and chicken hosts, boxes are checked next to “A/Viet Nam/1203/2004 (AY818135) (human hosts)” and “A/chicken/Viet Nam/Ncvd8/2003 (EF541407) (chicken hosts),” Clicking Update View displays the two sequences in the Viewer tab.
In the viewer, a model is displayed for the selected 2003-chicken H5. To compare the two models, the identifiers for both sequences are selected at the left of the sequence alignment, highlighting them in purple and instantiating a side-by-side viewer (Figure 4). In order to clarify the view, the Single Chain button is clicked, which displays only one chain of each trimer. Scrolling through the alignment, a non-conservative substitution at position 145 adjacent the receptor-binding domain is observed. At this position, there is a serine in the chicken host H5 and a leucine in the human host H5. Clicking on this residue in the sequence highlights in cyan the corresponding residue in each model with neighboring residues highlighted in yellow in the alignment. The side-by-side view of each model is also centered on the residue of interest and it and its neighbors are rendered as sticks on top of the ribbon diagram. From the viewer, the leucine appears exposed in the human host H5 (Figure 4). To investigate this observation further through additional modeling and sequence analysis, the models and sequences are downloaded by pressing the Download button, and selecting Models and Sequences, respectively, from the drop down list.
Because HA regulates IAV entry into the cell  and strongly impacts both pathogenicity and interspecies transmission [16, 17], new insights into HA structure and function are critical contributions to the study of influenza A. The HASP server was designed to enable researchers to visualize their HA sequences of interest in the three-dimensional context of the most related HA crystal structure. It employs computationally fast protocols for model generation, database navigation, and visualization appropriate for the large number of HA sequences and scalable with increasing sequence data. These protocols will be used in planned implementations of the HASP server to allow for models to be generated in real-time for user submitted sequences. The HASP server makes structural information on HA sequences quickly and easily accessible to all researchers, providing a valuable aid for interpreting data and generating new hypotheses. For example, structural information derived from the HASP pipeline was recently utilized in a study of a HA variant from the 1918 influenza pandemic  that resulted in nearly 50 million deaths worldwide . In that study, the models generated were used in docking studies with receptor analogs to assess potential changes in receptor binding specificity resulting from point mutations .
The primary aim of the HASP server is to provide a tool for viewing and comparing amino acid changes in HA subtypes at all levels of protein structure, from primary to quaternary, giving a complete and integrated view of those changes to facilitate understanding. Understanding these changes, which determine the efficiency of transmission, pathology, and ecology of these viruses, is of critical and vital importance to global public health.
Availability and requirements
The HASP server is available free of charge as a web application at: http://exon.niaid.nih.gov/HASP.html
Hemagglutinin Structure Prediction (server)
Root Mean Square Deviation
High-scoring segment pair
Influenza A virus
Protein Data Bank
Google Web Toolkit.
Clark NM, Lynch JP: Influenza: epidemiology, clinical features, therapy, and prevention. Semin Respir Crit Care Med. 2011, 32 (4): 373-392. 10.1055/s-0031-1283278.
Taubenberger JK, Kash JC: Influenza virus evolution, host adaptation, and pandemic formation. Cell Host Microbe. 2010, 7 (6): 440-451. 10.1016/j.chom.2010.05.009.
Suzuki Y, Ito T, Suzuki T, Holland RE, Chambers TM, Kiso M, Ishida H, Kawaoka Y: Sialic acid species as a determinant of the host range of influenza A viruses. J Virol. 2000, 74 (24): 11825-11831. 10.1128/JVI.74.24.11825-11831.2000.
Bao Y, Bolotov P, Dernovoy D, Kiryutin B, Zaslavsky L, Tatusova T, Ostell J, Lipman D: The influenza virus resource at the National Center for Biotechnology Information. J Virol. 2008, 82 (2): 596-601. 10.1128/JVI.02005-07.
Ha Y, Stevens DJ, Skehel JJ, Wiley DC: H5 avian and H9 swine influenza virus haemagglutinin structures: possible origin of influenza subtypes. EMBO J. 2002, 21 (5): 865-875. 10.1093/emboj/21.5.865.
Russell RJ, Gamblin SJ, Haire LF, Stevens DJ, Xiao B, Ha Y, Skehel JJ: H1 and H7 influenza haemagglutinin structures extend a structural classification of haemagglutinin subtypes. Virology. 2004, 325 (2): 287-296. 10.1016/j.virol.2004.04.040.
Wang G, Dunbrack RL: PISCES: a protein sequence culling server. Bioinformatics. 2003, 19 (12): 1589-1591. 10.1093/bioinformatics/btg224.
Chen VB, Arendall WB, Headd JJ, Keedy DA, Immormino RM, Kapral GJ, Murray LW, Richardson JS, Richardson DC: MolProbity: all-atom structure validation for macromolecular crystallography. Acta Crystallogr D Biol Crystallogr. 2010, 66 (Pt 1): 12-21.
Leaver-Fay A, Tyka M, Lewis SM, Lange OF, Thompson J, Jacak R, Kaufman K, Renfrew PD, Smith CA, Sheffler W, et al: ROSETTA3: an object-oriented software suite for the simulation and design of macromolecules. Methods Enzymol. 2011, 487: 545-574.
Squires B, Macken C, Garcia-Sastre A, Godbole S, Noronha J, Hunt V, Chang R, Larsen CN, Klem E, Biersack K, et al: BioHealthBase: informatics support in the elucidation of influenza virus hostpathogen interactions and virulence. Nucleic Acids Res. 2008, 36: D497-D503.
Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol. 1990, 215 (3): 403-410.
Benson DA, Karsch-Mizrachi I, Lipman DJ, Ostell J, Sayers EW: GenBank. Nucleic Acids Res. 2011, 39: D32-37. 10.1093/nar/gkq1079.
Katoh K, Toh H: Recent developments in the MAFFT multiple sequence alignment program. Brief Bioinform. 2008, 9 (4): 286-298. 10.1093/bib/bbn013.
Cumulative number of confirmed human cases for avian influenza A (H5N1) reported to WHO. 2003, [http://www.who.int/influenza/human_animal_interface/H5N1_cumulative_table_archives/en/index.html], -2012
Jiang S, Li R, Du L, Liu S: Roles of the hemagglutinin of influenza A virus in viral entry and development of antiviral therapeutics and vaccines. Protein Cell. 2010, 1 (4): 342-354. 10.1007/s13238-010-0054-6.
Stevens J, Donis RO: Influenza virus hemagglutinin-structural studies and their implications for the development of therapeutic approaches. Infect Disord Drug Targets. 2007, 7 (4): 329-335. 10.2174/187152607783018727.
Rungrotmongkol T, Yotmanee P, Nunthaboot N, Hannongbua S: Computational studies of influenza A virus at three important targets: hemagglutinin, neuraminidase and M2 protein. Curr Pharm Des. 2011, 17 (17): 1720-1739. 10.2174/138161211796355083.
Sheng ZM, Chertow DS, Ambroggio X, McCall S, Przygodzki RM, Cunningham RE, Maximova OA, Kash JC, Morens DM, Taubenberger JK: Autopsy series of 68 cases dying before and during the 1918 influenza pandemic peak. Proc Natl Acad Sci USA. 2011, 108 (39): 16416-16421. 10.1073/pnas.1111179108.
Taubenberger JK, Morens DM: 1918 Influenza: the mother of all pandemics. Emerg Infect Dis. 2006, 12 (1): 15-22.
We are grateful to Dr. Yentram Huyen and Mike Tartakovsky for their support of the project. We also thank Dr. Michael Dolan for helpful comments on HASP, and Dr. Meghan Coakley for assistance in writing this manuscript. This research was supported in part by the Office of Science Management and Operations of the National Institute of Allergy and Infectious Diseases. This study utilized the high-performance computational capabilities of the Biowulf Linux cluster at the National Institutes of Health, Bethesda, Maryland. (http://biowulf.nih.gov).
The authors declare that they have no competing interests.
XIA, EJD, JKT, and DEH conceived the HASP server. XIA designed the HASP server, created the modeling pipeline, performed the model quality assessment, and authored the manuscript. JD led the development and implementation of the HASP server. VG assisted JD with the development and implementation of the HASP server. All authors contributed to and edited the manuscript, and approved the final version.
Xavier I Ambroggio, Jennifer Dommer contributed equally to this work.
About this article
Cite this article
Ambroggio, X.I., Dommer, J., Gopalan, V. et al. HASP server: a database and structural visualization platform for comparative models of influenza A hemagglutinin proteins. BMC Bioinformatics 14, 197 (2013). https://doi.org/10.1186/1471-2105-14-197
- Influenza A
- Receptor binding
- Membrane fusion
- Sialic acid
- Homology modeling
- Molecular visualization