Abstract
Influenza A(H3N2) has been a major cause of seasonal influenza in humans since 1968, and has evolved by antigenic drift under the constantly changing human herd immunity. Increasing evidence suggests that the antigenic change occasionally occurred concomitant with the alterations of the N-glycosylation site profile and hemagglutination activity of the virion surface protein hemagglutinin (HA). However, the structural basis of these changes remains largely unclear. To address this issue, we performed molecular dynamics simulations of the glycosylated HA trimers of the A(H3N2), which has a novel pattern of Asn-X-Ser/Thr sequons unique in the new A(H3N2) epidemic clade 3C.2a and is characterized by attenuated ability to agglutinate nonhuman erythrocytes. Comparison of the equilibrated structures of the glycosylated HA trimers with and without the 3C.2a-specific mutations reveals that the mutations could induce a drastic reduction in the apical space for the ligand binding via glycan-shield rearrangement. The results suggest that the 3C.2a strain has evolved an HA structure that is advantageous for evading pre-existing antibodies, while also increasing the ligand binding specificity. These findings have structural implications for our understanding of the phenotypic changes, evolution, and fate of influenza A(H3N2).
Introduction
The hemagglutinin (HA) protein of influenza virus is a glycosylated type I integral membrane protein that protrudes from the mature virion surface and plays critical roles in viral interactions with hosts. The HA protein is synthesized in infected cells as a precursor HA0, and is subsequently cleaved by cellular proteases to HA1 and HA2 subunits that are covalently attached by a disulfide bond. The mature HA protein on the virion is composed of three pairs of the HA1/HA2 subunits (). The tip of the HA protein forms a globular structure, termed the globular head, and confers on the virus an ability to attach cells via interactions with the sialic acid-containing glycan moiety on the target cell surface (). Meanwhile, the HA globular head constitutes the major viral antigenic sites that induce neutralization antibodies in infected hosts. These functional and antigenic features drive sequence and structural variations, particularly near the receptor-binding site in the globular head, according to specific rules (; ). Importantly, the sequence variation on the globular head causes various phenotypic changes of viruses, including changes in antigenicity and receptor specificity. Therefore, it is critical to determine the structural changes in the HA globular head in order to understand the viral interplay with the hosts and evolution. Unfortunately, however, it is usually time consuming to characterize mutation-induced structural changes by experimental approaches alone.
Computational science is a rapidly growing area that now successfully complements the experimental and theoretical sciences in various fields, including life science. For example, recent advances in molecular dynamics (MD) simulation enable us to characterize changes in the three-dimensional structures of the mutated proteins in relatively short timescales compared with the experimental approaches (; ). The MD simulations have been used to disclose the structural basis of the adaptation and evolution of the highly mutable human immunodeficiency virus (HIV). This includes elucidation of the HIV structural changes associated with the phenotypic changes in viral neutralization sensitivity and receptor tropism (; , ; ), viral sensitivity to antiviral protein (), viral drug sensitivity (), viral growth in non-natural host cells (), and viral sensitivities to antibodies by drug-resistance mutations (; ).
In this study, we used the MD simulation to gain new insights into the roles of mutations in a recent epidemic variant of the influenza A(H3N2) viruses. The A(H3N2) viruses have emerged on 1968 in humans of southern Asia and were soon widespread in the world. Thereafter, the A(H3N2) has been a major cause of seasonal influenza in humans to date. During the 2014/15 epidemic season of influenza, a new A(H3N2) substrain had rapidly predominated in humans worldwide (). Notably, the hemagglutination activity of this substrain somehow could be measured with only a small portion of the viral population using a conventional hemagglutination assay with nonhuman erythrocytes (). The A(H3N2) substrain is characterized by alterations of the N-glycosylation sequons on the globular heads of the HA protein as compared with other A(H3N2) clades () and is now referred to as 3C.2a. The oligosaccharides on the HA protein play key roles in viral antigenicity (; ; ; ; ; ; ) and binding specificity/affinity to the cellular receptor (; ; ; ; ; ; ; ). However, it remains unclear how the 3C.2a mutations altered the HA structure and attenuated the hemagglutination activity with nonhuman erythrocytes. To address this issue, we here examined the structural effects of the four mutations in the globular heads using MD simulations. The obtained results predicted that the mutations could induce rearrangement of the glycan shield around the receptor-binding surface of the HA protein, leading to shrinkage of the ligand-accessible space.
Materials and Methods
Genetic Clades Determination
Genetic clades determination of the influenza virus in Japan has been performed routinely as the part of the work of the National Epidemiological Surveillance of Infectious Diseases in Japan1 and the Global Surveillance of Influenza in the WHO Reference Laboratories2. Information on the specimen collection and weekly report of HA type is available on the National Epidemiological Surveillance of Infectious Diseases3. Briefly, RNAs were extracted from viruses by using QIAamp viral RNA kit (QIAGEN, Dusseldorf, Germany). The HA genes were amplified from extracted RNAs by RT-PCRs using gene-specific primers (the primer sequences are available upon request) and SuperScript III One-step RT-PCR system with Platinum Taq (Invitrogen, Carlsbad, CA, USA). Sequencing reactions were performed with BigDye terminator kit (Applied Biosystems) and sequences were determined using 3730xl DNA analyzer (Applied Biosystems, Foster City, CA, USA). The genetic clades were determined by the generation of phylogenetic trees of the HA genes. The phylogenetic trees were constructed using MEGA 6 software () with the neighbor-joining method. The numbers of HA sequences obtained in the individual periods are 39, 83, 125, 108, 70, and 77 for the September 2013 – Jan. 2014, February 2014 – August 2014, September 2014 – January 2015, February 2015 – August 2015, September 2015 – January 2016, and February 2016 – August 2016, respectively. The nucleotide sequences used in this study are registered at GISAID, a publicly accessible influenza virus database4. Accession number at GISAID for the HA sequence used for the molecular modeling in this study is EPI543763 (A/Switzerland/9715293/2013 strain).
Molecular Modeling of a Glycosylated HA Trimer in the Ligand-free State
Three-dimensional (3-D) models for glycosylated extracellular domains of HA trimers of influenza A (H3N2) in the ligand-free state were constructed by homology modeling with Molecular Operating Environment (MOE) (Chemical Computing Group Inc., Montreal, QC, Canada). The crystal structure of the HA trimer of the influenza A/Victoria/361/2011 (H3N2) virus (PDB code: 4O5N; resolution: 1.75 Å; amino acid residues 4–325 and 330–502 for HA1 and HA2 peptides, respectively) was used as the modeling template. Obtained models were optimized by energy minimization using MOE and an Amber10: Extended Huckel Theory (EHT) force field implemented in MOE, which combines Amber10 and EHT bonded parameters for the large-scale energy minimization (; ). The high-mannose oligosaccharide Man5GlcNAc2 was added to potential N-glycosylation sites in HA using Online Glycoprotein Builder5.
MD Simulation of Glycosylated HA Trimer Models
Glycosylated HA trimer models in a ligand-free state were subjected to MD simulation essentially as described for simulations of HIV-1 gp120 (). MD simulations were performed by the PMEMD (Particle Mesh Ewald Molecular Dynamics) module in the AMBER 14 program package (), employing the Amber ff99SB-ILDN force field, a protein force field with improved side-chain torsion potentials (), the GLYCAM06 force field, a biomolecular force field for glycans (), and the TIP3P water model for simulations of aqueous solutions (). Bond lengths involving hydrogen were constrained with SHAKE, a constraint algorithm to satisfy a Newtonian motion (), and the time step for all MD simulations was set to 2 fs. A non-bonded cutoff of 10 Å was used. After heating calculations for 20 ps until 310 K using the NVT ensemble for the constant volume, temperature, and numbers of particles in the system, simulations were executed using the NPT ensemble for the constant pressure, temperature, and numbers of particles in the system at 1 atm, at 310 K, and in 150 mM NaCl for 100 ns. Root mean square deviations (RMSDs) between the heavy atoms of the two superposed proteins were used to measure the overall structural differences between the two proteins (). The RMSD was calculated using the cpptraj module in AmberTools 14, a trajectory analysis tool (). We used “Computer System for the Prediction of Mutations of Pathogens” at Research Center for Zoonosis Control, Hokkaido University for the MD simulations.
Calculation of Root Mean Square Fluctuation (RMSF)
We calculated RMSFs of individual components of the high mannose oligosaccharides Man5GlcNAc2 around the receptor binding site between 50 to 100 ns of MD simulations to quantify structural dynamics of glycans during the MD simulations. The average structures during the last 50 ns of MD simulations were used as reference structures for RMSF calculation. RMSFs were calculated as previously described (; , ; ) by using the ptraj module in Amber, a trajectory analysis tool ().
Results
Temporal Dynamics of Clade Populations of A(H3N2) Viruses Since September 2013
Over the 2014/15 influenza season, influenza surveillance reports from reference laboratories warned of a global epidemic of a new A(H3N2) variant population, termed 3C.2a. In Japan, the 3C.2a was first detected as a relatively minor population over the season from February to August in 2014, during which four A(H3N2) clades co-existed (Figure 1, upper panel). However, the 3C.2a rapidly became dominant at the beginning of the winter season in 2015, displaced pre-existing clades, and has continually predominated at the collection sites in Japan ever since, representing 76 of 77 (98.7%) clades in the season from February to August in 2016 (Figure 1, lower panel). In parallel, another A (H3N2) clade, 3C.3a, which co-existed in the same period as the 3C. 2a over the February to August 2014 season, became minor after February 2015 (Figure 1). These results are consistent with the reports on the global epidemic of the 3C.2a from the WHO Reference Laboratories6, and suggest the selective advantage of the 3C.2a for human-to-human transmission during the study period. Notably, hemagglutination activity of the 3C.2a was significantly attenuated when measured with a conventional hemagglutination assay (). Together, these data suggest that certain structural changes occurred in the HA protein to confer selective advantages on the 3C.2a.
FIGURE 1
Amino Acid Signatures of the 3C.2a HA Protein
Molecular epidemiological data from above study and the global surveillance of influenza in the WHO Reference Laboratories suggest that the A (H3N2) clades, 3C.2a and 3C.3a, were diversified from 3C.2 and 3C.3, respectively and thereafter the 3C.2a had dominated over the 3C.3a. Therefore, we examined differences in HA proteins between the clades 3C.2a and 3C.3a. The HA protein of the 3C.2a population has multiple amino acid substitutions as compared with the clade 3C.3a (amino acid numbers 3, 128, 138, 142, 144, 159, 160, 311, 326, and 489). In this study, we focused on the four substitutions around the receptor-binding surface on the tip of the HA protein, i.e., Ala128Thr, Asn144Ser, Ser159Tyr, and Lys160Thr (Figure 2). The Ala128Thr and Lys160Thr create new potential N-glycosylation sites, Asn-X-Ser/Thr, whereas Asn144Ser results in the loss of a single N-glycosylation site. Ser159Tyr is located adjacent to the Lys160Thr. The Ala128Thr substitution initially detected in the 3C.2 and had been preserved in its descendent, 3C.2a, whereas the other three substitutions newly emerged in the 3C.2a (Figure 2A). All the mutations are placed at or near the major antigenic elements for H3 HAs, such as antigen sites A (amino acids 140-146) and B (amino acids 155-160 and 188-198) (; ; ; ; ) (Figure 2B).
FIGURE 2
MD Simulations of Glycosylated HA Trimers in the Ligand-free State
To understand the structural impacts of the four mutations on the tip of the 3C.2a HA protein, we constructed molecular models of the glycosylated HA trimers with and without the mutations. Two structural models of the glycosylated HA trimers in a ligand-free state were constructed by homology modeling: a model for the A(H3N2) 2015/16 vaccine strain (A/Switzerland/9715293/2013) that belongs to the perishing clade 3C.3a (Figures 1, 2) and a model for A/Switzerland/9715293/2013 possessing the four 3C.2a-specific mutations on the HA globular heads (Ala128Thr, Asn144Ser, Ser159Tyr, and Lys160Thr). The obtained models were optimized by energy minimization using the MOE:EHT force field (; ) and were subjected to the MD simulations using the Amber ff99SB-ILDN force field () and the GLYCAM06 force field () as described previously ().
The structural dynamics of the glycosylated HA proteins in solution were monitored by RMSD between the initial model structure and the structures at given time points of the MD simulation (Figure 3). For the protein portions of the glycosylated HA molecules, the RMSDs sharply increased in the beginning and reached a near plateau after 20 ns of the MD simulations in both the 3C.3a and its mutant with 3C.2a mutations (Figures 3A,B). This RMSD profile was similar to that obtained with the glycosylated HIV-1 gp120 protein (). In addition, we monitored the RMSDs of the glycan portions of the glycosylated HA molecules. We obtained an RMSD profile that was similar to that of the protein portions: there was a sharp increase soon after the MD simulation onset followed by a near plateau after 20 ns (Figures 3C,D). These results suggest that the structural distortions of the amino acid residues and glycans of the initial models were relieved shortly after the start of MD simulation under thermodynamic driving forces in solution. The data also predict that the HA structure can reach a state of thermodynamic equilibrium in solution.
FIGURE 3
Finally, we compared the 3-D structure of the HA trimer in a state of thermodynamic equilibrium during the MD simulations (100 ns) between the 3C.3a and its mutant (Figure 4A). A marked structural difference between the 3C.3a and its mutants was detected on the globular heads of the HA protein at 100 ns of two MD simulations: the glycan moieties were placed more compactly in the mutant than in the 3C.3a (Figures 4A,B). Consequently, the space for the access of ligands around the receptor-binding site was more confined in the mutant (Figure 4C). The apical spaces for ligand binding could be influenced not only by the arrangements of glycans but also by fluctuation of the glycans. Therefore, we examined structural fluctuations of the individual glycans around the receptor binding site by calculating RMSF using snapshots of structures from 50 to 100 ns of each MD simulation. In the present HA models, the high mannose oligosaccharides Man5GlcNAc2 were attached on the asparagine residues of potential N-glycosylation sites around the receptor binding site (Figure 5A). We calculated RMSFs of individual components of the oligosaccharides, such as GlcNAc and Man, at positions 1 to 7 (Figure 5B). Notably, RMSFs of the glycan components at given positions, and the profiles of RMSFs at the seven positions were similar among the four N-glycosylation sites. The data indicate that all glycans fluctuated in a similar fashion around the receptor binding site during the MD simulations. These results suggest that the apical spaces for ligand binding are mainly influenced by arrangements of glycans and that the spaces are constantly different between the two models during the MD simulations.
FIGURE 4
FIGURE 5
Discussion
In this report, we studied the structural impact of mutations residing in the HA globular head of a recent epidemic clade population of the A(H3N2). This clade, termed 3C.2a, predominated over other co-existing A(H3N2) clades including the clade 3C.3a during the 2014/15 and 2015/16 seasons, and retained a unique set of sequon mutations in the HA globular heads (Figures 1, 2). To analyze the HA structures under near physiological conditions, we performed MD simulations of the glycosylated HA trimers and obtained structures in a state of thermodynamic equilibrium in solution (Figure 3). Comparisons of the obtained structures revealed marked changes in the glycan shield around the receptor-binding site (Figures 4, 5). This finding has important implications for our understanding of phenotypic changes, evolution, and fate of the influenza virus A(H3N2).
First, the present results indicate that the 3C.2a-specific mutations have an impact on the immunological features of the HA protein. Our MD simulations show that the glycans on the HA globular head are rearranged by the mutations into a configuration markedly distinct from that of the 2015/16 vaccine strain (Figures 4, 5). Therefore, these mutations could induce changes in the HA antigenicity. These structural findings are consistent with previous reports demonstrating critical roles of the oligosaccharides on the HA protein in the viral antigenicity of influenza viruses (; ; ; ; ; ; ). Moreover, we found that the glycan rearrangement resulted in shrinkage of the access space on the top of the HA protein near the receptor-binding site of the globular head (Figure 4). This structural change could cause steric hindrance for binding of the antibodies directed to the antigen site B. Thus, it is likely that the HA protein of the 3C.2a had a selective advantage in evading antibodies against receptor-binding site as compared with the 3C.3a HA. These possibilities are consistent with the rapid spread and predominance of the clade C3.2a over the clade 3C.3a in human populations during the study period (Figure 1) (). Antigenic analysis of these viruses should be done to make sure these possibilities.
Secondly, the present results indicate that the 3C.2a-specific mutations impact the receptor specificity of the HA protein. Our MD simulations showed that the new arrangement of glycans on the HA globular head could shrink the space for the access of the sialic acid-containing glycan moiety on the target cell surface (Figure 4). This structural change could cause an increase in the ligand specificity and/or affinity of the HA protein, and thereby lead to a preference for receptors on the human respiratory organs, but not for receptors on the nonhuman erythrocytes. These structural findings are consistent with the attenuation of hemagglutination activity of the 3C.2a HA protein when assessed with a conventional hemagglutination assay using nonhuman erythrocytes (). Similarly, previous studies have demonstrated that the oligosaccharides on the HA protein surface play key roles in the binding specificity and the affinity to infection receptors (; ; ; ; ; ; ; ).
Thirdly, the present findings have implications in terms of the adaptive evolution of A (H3N2). Previous reports have highlighted the importance of changes in the glycan shield for viral adaptation (). Interestingly, the numbers of potential N-glycosylation sites in the HA protein have been continuously increasing since 1968: only two sequons existed –in the initial strains, whereas more than 7–10 sequons are common in the present epidemic strains (). This change is likely to be a basic strategy by which A (H3N2) has maintained its presence in human populations with changing herd immunity over the last 48 years. However, the change also runs the risk of creating an evolutionary “dead-end” for the virus, because the acquisition of new glycans on the HA globular head increases the chances for steric hindrance during receptor binding. Notably, both the present and recent studies suggest that the HA protein possessing the 3C.2a-type glycan shield would significantly affect ligand specificity (Figure 4) (). Thus, the 3C.2a has evolved an HA structure that is advantageous for evading pre-existing antibodies, while also reducing ligand affinity to nonhuman glycan moieties.
Finally, the present have structural implications in the fate of influenza A (H3N2) viruses. Our study indicates that an increase in the number of the glycan moieties around the receptor-binding site could induce a drastic reduction in the apical space for ligand binding (Figure 4). As discussed above, continuous interactions between A (H3N2) and human immunity seem to force the virus to serially increase the numbers of glycans on the HA protein. This eventually may reduce the HA binding affinity via an increase in steric hindrance even to the human ligands, and thereby decrease the replication fitness of the virus for spread among humans. Thus, if the overall tendency of HA to increase the numbers of glycosylation sites () continues in the future, it may eventually put the influenza A (H3N2) viruses in an evolutionary cul-de-sac.
Statements
Author contributions
MY, SW, TO, and HS conceived and designed the study. SF and MS performed sequencing and phylogenetic analysis. KI and HS prepared the computing environment. MY performed MD simulations. HS prepared the manuscript. All authors read and approved the final manuscript.
Funding
This study was supported by a Grant-in-Aid for emerging and reemerging infectious diseases to HS (Grant Number: H27-Shinkogyousei-shitei-002) from the Ministry of Health, Labour and Welfare of Japan, a financial support from the Joint Usage/Research Center Program at Research Center for Zoonosis Control, Hokkaido University, a Grant-in-Aid for Scientific Research (C) to MY (Grant Number:16K08824) from the Ministry of Education, Culture, Sports, Science and Technology, and a Grant-in-Aid for the Research Program on Re-emerging Infectious Diseases to TO from the Japan Agency for Medical Research and development, AMED.
Acknowledgments
We thank Ms. Hiromi Nakamura for technical assistance.
Conflict of interest
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.
Footnotes
1.^http://www.nih.go.jp/niid/en/influenza-e.html
2.^http://www.who.int/influenza/en/
3.^http://www.nih.go.jp/niid/en/influenza-e/2099-idsc/iasr-flu-e/6791-iasr-inf-e20160925.html
References
1
AbeY.TakashitaE.SugawaraK.MatsuzakiY.MurakiY.HongoS. (2004). Effect of the addition of oligosaccharides on the biological activities and antigenicity of influenza A/H3N2 virus hemagglutinin.J. Virol.789605–9611. 10.1128/JVI.78.18.9605-9611.2004
2
AlamM.KuwataT.ShimuraK.YokoyamaM.Ramirez ValdezK. P.TanakaK.et al (2016). Enhanced antibody-mediated neutralization of HIV-1 variants that are resistant to fusion inhibitors.Retrovirology13:70. 10.1186/s12977-016-0304-7
3
AytayS.SchulzeI. T. (1991). Single amino acid substitutions in the hemagglutinin can alter the host range and receptor binding properties of H1 strains of influenza A virus.J. Virol.653022–3028.
4
BlackburneB. P.HayA. J.GoldsteinR. A. (2008). Changing selective pressure during antigenic changes in human influenza H3.PLoS Pathog.4:e1000058. 10.1371/journal.ppat.1000058
5
CaseD. A.BabinV.BerrymanJ. T.BetzR. M.CaiQ.CeruttiD. S.et al (2014). AMBER 14.San Francisco, CA: University of California.
6
CaseD. A.CheathamT. E.IIIDardenT.GohlkeH.LuoR.MerzK. M.et al (2005). The Amber biomolecular simulation programs.J. Comput. Chem.261668–1688. 10.1002/jcc.20290
7
DasS. R.PuigboP.HensleyS. E.HurtD. E.BenninkJ. R.YewdellJ. W. (2010). Glycosylation focuses sequence variation in the influenza A virus H1 hemagglutinin globular domain.PLoS Pathog.6:e1001211. 10.1371/journal.ppat.1001211
8
de VriesR. P.de VriesE.BoschB. J.de GrootR. J.RottierP. J.de HaanC. A. (2010). The influenza A virus hemagglutinin glycosylation state affects receptor-binding specificity.Virology40317–25. 10.1016/j.virol.2010.03.047
9
DeomC. M.CatonA. J.SchulzeI. T. (1986). Host cell-mediated selection of a mutant influenza A virus that has lost a complex oligosaccharide from the tip of the hemagglutinin.Proc. Natl. Acad. Sci. U.S.A.833771–3775. 10.1073/pnas.83.11.3771
10
GambaryanA. S.MarininaV. P.TuzikovA. B.BovinN. V.RudnevaI. A.SinitsynB. V.et al (1998). Effects of host-dependent glycosylation of hemagglutinin on receptor-binding properties on H1N1 human influenza A virus grown in MDCK cells and in embryonated eggs.Virology247170–177. 10.1006/viro.1998.9224
11
GerberP. R.MullerK. (1995). MAB, a generally applicable molecular force field for structure modelling in medicinal chemistry.J. Comput. Aided Mol. Des.9251–268. 10.1007/BF00124456
12
GerhardW.YewdellJ.FrankelM. E.WebsterR. (1981). Antigenic structure of influenza virus haemagglutinin defined by hybridoma antibodies.Nature290713–717. 10.1038/290713a0
13
GuntherI.GlatthaarB.DollerG.GartenW. (1993). A H1 hemagglutinin of a human influenza A virus with a carbohydrate-modulated receptor binding site and an unusual cleavage site.Virus Res.27147–160. 10.1016/0168-1702(93)90078-2
14
HaY.StevensD. J.SkehelJ. J.WileyD. C. (2003). X-ray structure of the hemagglutinin of a potential H3 avian progenitor of the 1968 Hong Kong pandemic influenza virus.Virology309209–218. 10.1016/S0042-6822(03)00068-0
15
HikichiY.YokoyamaM.TakemuraT.FujinoM.KumakuraS.MaedaY.et al (2016). Increased HIV-1 sensitivity to neutralizing antibodies by mutations in the Env V3-coding region for resistance to CXCR4 antagonists.J. Gen. Virol.972427–2440. 10.1099/jgv.0.000536
16
JorgensenW. L.ChandrasekharJ.MaduraJ. D.ImpeyR. W.KleinM. L. (1983). Comparison of simple potential functions for simulating liquid water.J. Chem. Phys.79926–935. 10.1063/1.445869
17
KirschnerK. N.YongyeA. B.TschampelS. M.Gonzalez-OuteirinoJ.DanielsC. R.FoleyB. L.et al (2008). GLYCAM06: a generalizable biomolecular force field. Carbohydrates.J. Comput. Chem.29622–655. 10.1002/jcc.20820
18
KnossowM.DanielsR. S.DouglasA. R.SkehelJ. J.WileyD. C. (1984). Three-dimensional structure of an antigenic mutant of the influenza virus haemagglutinin.Nature311678–680. 10.1038/311678a0
19
KoelB. F.BurkeD. F.BestebroerT. M.van der VlietS.ZondagG. C.VervaetG.et al (2013). Substitutions near the receptor binding site determine major antigenic change during influenza virus evolution.Science342976–979. 10.1126/science.1244730
20
KuwataT.TakakiK.YoshimuraK.EnomotoI.WuF.OurmanovI.et al (2013). Conformational epitope consisting of the V3 and V4 loops as a target for potent and broad neutralization of simian immunodeficiency viruses.J. Virol.875424–5436. 10.1128/JVI.00201-13
21
LiaoH. Y.HsuC. H.WangS. C.LiangC. H.YenH. Y.SuC. Y.et al (2010). Differential receptor binding affinities of influenza hemagglutinins on glycan arrays.J. Am. Chem. Soc.13214849–14856. 10.1021/ja104657b
22
Lindorff-LarsenK.PianaS.PalmoK.MaragakisP.KlepeisJ. L.DrorR. O.et al (2010). Improved side-chain torsion potentials for the Amber ff99SB protein force field.Proteins781950–1958. 10.1002/prot.22711
23
MatrosovichM.ZhouN.KawaokaY.WebsterR. (1999). The surface glycoproteins of H5 influenza viruses isolated from humans, chickens, and wild aquatic birds have distinguishable properties.J. Virol.731146–1155.
24
MiyamotoT.NakayamaE. E.YokoyamaM.IbeS.TakeharaS.KonoK.et al (2012). The carboxyl-terminus of human immunodeficiency virus type 2 circulating recombinant form 01_AB capsid protein affects sensitivity to human TRIM5alpha.PLoS ONE7:e47757. 10.1371/journal.pone.0047757
25
NaganawaS.YokoyamaM.ShiinoT.SuzukiT.IshigatsuboY.UedaA.et al (2008). Net positive charge of HIV-1 CRF01_AE V3 sequence regulates viral sensitivity to humoral immunity.PLoS ONE3:e3206. 10.1371/journal.pone.0003206
26
OdeH.NakashimaM.KitamuraS.SugiuraW.SatoH. (2012). Molecular dynamics simulation in virus research.Front. Microbiol.3:258. 10.3389/fmicb.2012.00258
27
OhuchiM.OhuchiR.FeldmannA.KlenkH. D. (1997). Regulation of receptor binding affinity of influenza virus hemagglutinin by its carbohydrate moiety.J. Virol.718377–8384.
28
PingJ.LiC.DengG.JiangY.TianG.ZhangS.et al (2008). Single-amino-acid mutation in the HA alters the recognition of H9N2 influenza virus by a monoclonal antibody.Biochem. Biophys. Res. Commun.371168–171. 10.1016/j.bbrc.2008.04.045
29
RyckaertJ.-P.CiccottiG.BerendsenH. J. C. (1977). Numerical integration of the cartesian equations of motion of a system with constraints: molecular dynamics of n-alkanes.J. Comput. Phys.23327–341. 10.1016/0021-9991(77)90098-5
30
SaitoT.NakayaY.SuzukiT.ItoR.SaitoT.SaitoH.et al (2004). Antigenic alteration of influenza B virus associated with loss of a glycosylation site due to host-cell adaptation.J. Med. Virol.74336–343. 10.1002/jmv.20178
31
SatoH.YokoyamaM.TohH. (2013). Genomics and computational science for virus research.Front. Microbiol.4:42. 10.3389/fmicb.2013.00042
32
SkowronskiD. M.SabaiducS.ChambersC.EshaghiA.GubbayJ. B.KrajdenM.et al (2016). Mutations acquired during cell culture isolation may affect antigenic characterisation of influenza A(H3N2) clade 3C.2a viruses.Euro Surveill.21:30112. 10.2807/1560-7917.ES.2016.21.3.30112
33
SmithD. J.LapedesA. S.de JongJ. C.BestebroerT. M.RimmelzwaanG. F.OsterhausA. D.et al (2004). Mapping the antigenic and genetic evolution of influenza virus.Science305371–376. 10.1126/science.1097211
34
TamuraK.StecherG.PetersonD.FilipskiA.KumarS. (2013). MEGA6: molecular evolutionary genetics analysis version 6.0.Mol. Biol. Evol.302725–2729. 10.1093/molbev/mst197
35
TsuchiyaE.SugawaraK.HongoS.MatsuzakiY.MurakiY.LiZ. N.et al (2002). Effect of addition of new oligosaccharide chains to the globular head of influenza A/H2N2 virus haemagglutinin on the intracellular transport and biological activities of the molecule.J. Gen. Virol.83(Pt 5)1137–1146. 10.1099/0022-1317-83-5-1137
36
WangC. C.ChenJ. R.TsengY. C.HsuC. H.HungY. F.ChenS. W.et al (2009). Glycans on influenza hemagglutinin affect receptor binding and immune response.Proc. Natl. Acad. Sci. U.S.A.10618137–18142. 10.1073/pnas.0909696106
37
WangW.LuB.ZhouH.SuguitanALJrChengX.SubbaraoK.et al (2010). Glycosylation at 158N of the hemagglutinin protein and receptor binding specificity synergistically affect the antigenicity and immunogenicity of a live attenuated H5N1 A/Vietnam/1203/2004 vaccine virus in ferrets.J. Virol.846570–6577. 10.1128/JVI.00221-10
38
WanzeckK.BoydK. L.McCullersJ. A. (2011). Glycan shielding of the influenza virus hemagglutinin contributes to immunopathology in mice.Am. J. Respir. Crit. Care Med.183767–773. 10.1164/rccm.201007-1184OC
39
WebsterR. G.LaverW. G. (1980). Determination of the number of nonoverlapping antigenic areas on Hong Kong (H3N2) influenza virus hemagglutinin with monoclonal antibodies and the selection of variants with potential epidemiological significance.Virology104139–148. 10.1016/0042-6822(80)90372-4
40
WileyD. C.SkehelJ. J. (1987). The structure and function of the hemagglutinin membrane glycoprotein of influenza virus.Annu. Rev. Biochem.56365–394. 10.1146/annurev.bi.56.070187.002053
41
WileyD. C.WilsonI. A.SkehelJ. J. (1981). Structural identification of the antibody-binding sites of Hong Kong influenza haemagglutinin and their involvement in antigenic variation.Nature289373–378. 10.1038/289373a0
42
YokoyamaM.NaganawaS.YoshimuraK.MatsushitaS.SatoH. (2012). Structural dynamics of HIV-1 envelope Gp120 outer domain with V3 loop.PLoS ONE7:e37530. 10.1371/journal.pone.0037530
43
YokoyamaM.NomaguchiM.DoiN.KandaT.AdachiA.SatoH. (2016). In silico analysis of HIV-1 Env-gp120 reveals structural bases for viral adaptation in growth-restrictive cells.Front. Microbiol.7:110. 10.3389/fmicb.2016.00110
44
YuanY.YokoyamaM.MaedaY.TerasawaH.HaradaS.SatoH.et al (2013). Structure and dynamics of the gp120 V3 loop that confers noncompetitive resistance in R5 HIV-1(JR-FL) to maraviroc.PLoS ONE8:e65115. 10.1371/journal.pone.0065115
Summary
Keywords
MD simulation, influenza A(H3N2), HA protein, N-linked glycans, mutations, structural change
Citation
Yokoyama M, Fujisaki S, Shirakura M, Watanabe S, Odagiri T, Ito K and Sato H (2017) Molecular Dynamics Simulation of the Influenza A(H3N2) Hemagglutinin Trimer Reveals the Structural Basis for Adaptive Evolution of the Recent Epidemic Clade 3C.2a. Front. Microbiol. 8:584. doi: 10.3389/fmicb.2017.00584
Received
05 January 2017
Accepted
21 March 2017
Published
10 April 2017
Volume
8 - 2017
Edited by
Aeron Hurt, WHO Collaborating Centre for Reference and Research on Influenza, Australia
Reviewed by
Hirotaka Ode, National Hospital Organization Nagoya Medical Center, Japan; Hidekatsu Iha, Oita University, Japan
Updates
Copyright
© 2017 Yokoyama, Fujisaki, Shirakura, Watanabe, Odagiri, Ito and Sato.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Hironori Sato, hirosato@nih.go.jp
This article was submitted to Virology, a section of the journal Frontiers in Microbiology
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.