HA198 Mutations in H9N2 Avian Influenza: Molecular Dynamics Insights into Receptor Binding

doi:10.21203/rs.3.rs-5297856/v1

Download PDF

Research Article

HA198 Mutations in H9N2 Avian Influenza: Molecular Dynamics Insights into Receptor Binding

https://doi.org/10.21203/rs.3.rs-5297856/v1

This work is licensed under a CC BY 4.0 License

Version 1

posted

You are reading this latest preprint version

The H9N2 avian influenza virus is widely disseminated in poultry and poses a zoonotic threat, despite vaccination efforts. Mutations at residue 198 of hemagglutinin (HA) are critical for antigenic variation and receptor-binding specificity, but the underlying molecular mechanisms remain unclear. This study explores the molecular mechanisms by which mutations at the HA 198 site affect the antigenicity, receptor specificity, and binding affinity of the H9N2 virus.

Using the sequence of the A/Chicken/Jiangsu/WJ57/2012 strain, we constructed recombinant H9N2 viruses, including rWJ57, rWJ57/HA_198A, and rWJ57/HA_198T, using reverse genetics. These variants were analyzed through hemagglutination inhibition (HI) assays, solid-phase receptor binding assays, and enzyme-linked immunosorbent assays (ELISA). Additionally, molecular dynamics (MD) simulations were performed to further dissect the atomic-level interactions between HA and sialic acids (SA).

The results showed that HA_198V and HA_198T mutations enhanced viral binding to human α2,6-linked SA, while HA_198A exhibited a preference for avian α2,3 SA. These mutations altered interactions with non-neutralizing antibodies but did not affect neutralizing antibodies, potentially contributing to immune evasion. MD simulations revealed that HA_198V/T formed stable complexes with the α2,6 SA, mediated by specific residues and water bridges, whereas HA_198A formed stable complexes with the α2,3 SA. Interestingly, residue 198 interacted with the α2,6 SA via water bridges but had minimal direct interaction with the α2,3 SA.

This study provides new insights into the molecular basis of receptor specificity, binding affinity, and antigenic drift in H9N2 viruses, highlighting the critical role of HA 198 mutations in regulating host adaptation. These findings are of great significance for H9N2 virus surveillance, vaccine development, and zoonotic transmission risk assessment.

H9N2 avian influenza virus

mutations at residue 198

hemagglutinin

receptor binding

molecular dynamics simulations

H9N2 avian influenza virus (AIV) has been identified as the predominant low-pathogenic strain in China [1], posing significant challenges to both the poultry industry and public health. This virus not only co-infects with bacterial pathogens such as Escherichia coli, Staphylococcus aureus, and Haemophilus paragallinarum [2–4], but also serves as a genetic reservoir for human-infecting influenza subtypes, including H3N8, H7N9, H10N3, and H10N8 [5–9]. Despite China's ongoing H9N2 vaccination efforts for more than two decades, the virus continues to circulate widely, driven by rapid antigenic variation [1]. The hemagglutinin (HA) protein, the virus's main antigen, is subject to mutations in its receptor-binding domain (RBD), which lead to significant changes to the virus’s antigenic properties [10]. One key site, amino acid 198, is especially important due to its frequent non-conservative mutations [11].

Previous research has shown that H9N2 vaccine strains commonly acquire the A198V mutation after multiple passages in specific-pathogen-free (SPF) chickens, leading to a significant reduction in hemagglutinin inhibition (HI) titers and an increase in receptor-binding activity [12–14]. This mutation is also linked to the virus's ability to evade antibody-mediated immunity, further complicating disease control efforts [15, 16]. The mutation at amino acid 198 within the HA protein, which is often found at this receptor-binding site, can also influence the binding affinity of H9N2 to sialic acid receptors [17]. Despite these findings, the molecular mechanisms by which mutations at the 198 site affect receptor-binding affinity remain poorly understood, underscoring the need for deeper investigation.

With recent advances in biological detection technologies, the precision of methods used to study virus-receptor interactions has significantly improved. Traditional techniques such as erythrocyte adsorption assays have been replaced by more advanced methods like solid-phase direct binding assays, glycan microarrays, surface plasmon resonance, and X-ray crystallography [18–20]. However, while these biological detection methods provide valuable insights, they are often complemented by computational approaches such as molecular docking and molecular dynamics (MD) simulations, which can simulate interactions between viral proteins and receptors under natural conditions [21, 22]. Although these computational techniques are widely used in drug discovery, they are underutilized in research on virus-host interactions, particularly in influenza viruses.

This study aims to investigate the molecular basis of receptor-binding alterations caused by mutations at position 198 in the HA protein of H9N2 AIV. Using the epidemic strain A/Chicken/Jiangsu/WJ57 (H9N2, WJ57) with 198V as the backbone [23], we employed reverse genetics to generate recombinant viruses carrying V198A and V198T mutations. The effects of these mutations on receptor-binding activity were assessed experimentally, while molecular docking and MD simulations were used to analyze the interactions between the HA protein and sialic acid receptors. Our findings reveal the molecular mechanisms by which mutations at the 198 site affect receptor-binding affinity and sialic acid specificity, contributing to a better understanding of the genetic evolution of H9N2 AIV and its antigenic variation.

Ethical compliance

The SPF chicken embryos were purchased from Nanjing Biology Medical Factory Qian Yuan-hao Biological Co., Ltd. All experimental procedures were approved by the Jiangsu Administrative Committee for Laboratory (SYXK 2016-0020) and were performed by the Jiangsu Laboratory Animal Welfare and Ethics guidelines of the Department of Science and Technology of Jiangsu Province.

Viruses and cells culture

The sequences of the H9N2 vaccine strain A/chicken/Jiangsu/WJ57/2012 (WJ57), published in the GenBank database under accession numbers KP893703 to KP893708 and KJ000709 to KJ000710, were synthesized by SYNBIO Technologies (Suzhou, China). The recombinant WJ57 strain (rWJ57) was generated using reverse genetics and stored at the Jiangsu Key Laboratory for High-Tech Research and Development of Veterinary Biopharmaceuticals. COS-1 cells were purchased from ATCC (Manassas, VA, USA), cultured in Dulbecco's Modified Eagle's Medium (DMEM) (Gibco, BRL, Grand Island, USA) supplemented with 10% fetal calf serum (Gemini, Woodland, CA, USA), and incubated at 37°C with 5% CO₂.

HA gene sequencing analysis

The frequency and prevalence of residues at amino acid position 198 in HA were analyzed by comparing 4,093 full-length H9N2 HA sequences from Chinese isolates deposited in GenBank from 2013 to 2022 (up to January 1, 2023).

Generation of recombinant H9N2 viruses

Primers for site-directed mutagenesis and homologous recombination were designed based on the HA sequence of the WJ57 virus using Primer 5.0 software (Primer-E Ltd., Plymouth, UK). Overlap PCR was used to amplify the mutant HA genes, which were then inserted into the transcriptional/expression vector pHW2000 using the ClonExpress II One Step Cloning Kit (Vazyme Biotech, Nanjing, China). After transformation into DH5α cells, ten colonies were randomly selected, cultured, and subjected to plasmid DNA extraction (TianGen, Beijing, China). The constructs were then sequenced (GENEWIZ, Suzhou, China) to confirm their accuracy. The primers used in the amplification process are listed in Table 1.

Table 1. Primers used for construction of recombinant plasmids

Primer	Primer sequence（5’-3’）
^apHW-HA-F	gtcgacctccgaagttgggggggAGCAAAAGCAGGGGAATTTC
^apHW-HA-R	ggcattttgggccgccgggttattAGTAGAAACAAGGGTGTTTTTGCC
^b V198A-F	CCACCGATACTGCACAGACAAATCT
^b V198A-R	AGATTTGTCTGTGCAGTATCGGTGG
^b V198T-F	CCACCGATACTACGCAGACAAATCT
^b V198T-R	AGATTTGTCTGCGTAGTATCGGTGG
^c pHWvector-F	ACCCGGCGGCCCAAAATGCC
^c pHWvector-R	CCCCAACTTCGGAGGTCGAC

^a Primers employed for pHW2000 construction using a binucleotide cloning strategy to amplify the full-length cDNA of the HA gene. Lowercase letters denote the homologous arm sequences from the pHW2000 vector; b Primers used for site-specific mutations in the HA gene, with sequences marked by a single underline indicating the nucleotides subjected to mutation; c Primers utilized for amplifying the linearized pHW2000 vector

As described by Zhu et al.[14], the H9N2 recombinant viruses, rWJ/HA_198A and rWJ/HA_198T, were identified and rescued using the rWJ57 strain backbone. Briefly, a total of 2.4 µg of eight plasmids, mixed in equal amounts, was added to 100 µL of Opti-MEM medium (Gibco, BRL, Grand Island, USA). Then, 5 µL of X-treme GENE HP DNA Transfection Reagent (Roche, Basel, Switzerland) was introduced to the mixture, which was incubated at room temperature for 15 minutes. The mixture was then added to COS-1 cells cultured in 12-well plates at 70–80% confluency. The cells were incubated at 37°C with 5% CO₂ for 24 hours, after which 2 µg/mL of TPCK-trypsin (Sigma, St. Louis, MO, USA) was added to each well. Following a 30-hour transfection period, the supernatants were collected and used to propagate the virus in 10-day-old specific pathogen-free (SPF) embryonated chicken eggs. Finally, the rescued virus was subjected to an HA assay, and the HA genes were sequenced at GENEWIZ (Suzhou, China) to confirm the accuracy of the intended mutations.

Antisera preparation

Antisera for the rWJ57 strain were prepared as previously described [11, 14] . Briefly, 3-week-old SPF chickens were immunized twice by subcutaneous injection of 0.3 mL of an oil-emulsified, inactivated whole-virus vaccine of the rWJ57 strain, which was inactivated by adding 0.2% formalin (v/v) for 24 hours at 37°C. Antisera were collected from the chickens two weeks after the booster vaccination.

Hemagglutination inhibition (HI) assay

Antisera for the rWJ57 strain were treated with cholera filtrate (Sigma-Aldrich, St. Louis, MO, USA) to remove nonspecific hemagglutination inhibitors before conducting the HI assay on three recombinant viruses: rWJ57, rWJ/HA_198A and rWJ/HA_198T. The assay was performed using 4 hemagglutination units (HAU) of H9N2 and 1% (v/v) chicken erythrocytes [14, 24].

Enzyme‑linked immunosorbent assay (ELISA)

The ELISA assay was conducted according to previously described methods [14, 24]. Sucrose gradient-purified viruses—rWJ57, rWJ/HA_198A and rWJ/HA_198T—were diluted in PBS and seeded into Nunc-Immuno MaxiSorp 96-well plates (Corning, NY, USA) at a concentration of 16 HAU per well. Following overnight incubation at 4°C, the wells were blocked with PBS-nonfat dry milk (Beyotime, Shanghai, China). Serial two-fold dilutions of antisera against the rWJ57 virus in chickens were prepared using PBS containing 0.05% Tween 20 and then added to the wells, followed by a 2-hour incubation at 37°C. After washing, a goat anti-chicken horseradish peroxidase antibody (Bioss, Beijing, China) was added, and the plates were incubated for 1.5 hours at 37°C. Finally, TMB (3,3′,5,5′-tetramethylbenzidine) substrate (Beyotime, Shanghai, China) was added, and the reaction was stopped using TMB stop substrate (Beyotime, Shanghai, China).

Absorbance was measured at 450 nm using an automated ELISA plate reader (model EL311SX; Biotek, Winooski, VT, USA). The area under the curve (AUC) for each virus was calculated using GraphPad Prism 9.5 software (San Diego, CA, USA), based on virus-antibody binding above the corresponding negative control.

Receptor-destroying enzyme (RDE) assay and Solid-phase binding assay

The RDE assay was performed as previously described [14, 25]. Chicken erythrocytes were treated with varying concentrations of α2-3,6,8 neuraminidase (New England Biolabs, Beverly, MA, USA) for 1 hour at 37°C. After washing with PBS, 1% (v/v) solutions of the treated erythrocytes were added to 4 HAU of each virus (determined using untreated chicken erythrocytes). Agglutination was measured after 1 hour of incubation.

To further evaluate the receptor-binding properties of the H9N2 viruses, we performed solid-phase binding assays using the synthetic sialyl glycopolymers Neu5Ac 2-3 Galβ1-4 GlcNAc (3ʹSLN)-PAA-biotin and Neu5Ac 2-6 Galβ1-4 GlcNAc (6ʹSLN)-PAA-biotin (GlycoTech), based on the method described by Zhang et al. [26, 27]. The glycan analogs were serially diluted in PBS and added to streptavidin-coated 96-well plates (Pierce). Plates were blocked with PBS containing 2% skim milk, and 64 HA units of virus were applied per well, followed by chicken antiserum as the primary antibody. Detection was done by adding HRP-conjugated rabbit anti-chicken IgG and tetramethylbenzidine substrate, with the reaction stopped by 1 M H₂SO₄. Absorbance at 450 nm was measured, and all samples were tested in triplicate.

Molecular docking and dynamics simulation

The structures of HA were predicted using the Swiss-Model. The protonation state of all compounds was set to pH 7.4, and 3D structures were generated using Open Babel [28]. AutoDock Tools (ADT3) were used to prepare and parameterize the receptor protein and ligands. Docking grid files were constructed with AutoGrid from the sitemap, and docking simulations were performed using AutoDock Vina (1.2.0) [29, 30]. The optimal pose was selected for protein-ligand interaction analysis. Finally, protein-ligand interaction figures were generated using PyMOL. The ligands used in this study were α2,6 sialic acid (SA) (6ʹSLN, CAS no. 501427-93-6) and α2,3 SA (3ʹSLN, CAS no. 501427-92-5).

MD simulations were conducted using the non-commercial version 2022.1 of Desmond/Maestro software (https://www.deshawresearch.com). TIP3P water molecules were added to the systems, which were neutralized with a 0.15 M NaCl solution. After system minimization and relaxation, a 100-ns production simulation was performed using an isothermal-isobaric ensemble at 300 K and 1 bar. Trajectory coordinates were recorded every 100 ps throughout the simulation. Molecular dynamics analysis was performed using the Simulation Interaction Diagram feature of Desmond software.

Statistical analysis

All data were expressed as mean ± SD. Statistical comparisons between groups were made using Student’s t-test in GraphPad Prism 9.5 (GraphPad, USA), with P-values less than 0.05 considered statistically significant.

Analysis of the natural mutation rate of residue 198

An analysis of 4,093 avian H9N2 AIV HA protein sequences at position 198, published in GenBank from 2013 to 2022, revealed three prevalent mutations: 198T, 198A, and 198V. The most frequent was 198T, representing 72.4% of the sequences, followed by 198A (20.7%), while 198V was the least common (6.9%) (Figure 1A). Between 2013 and 2015, the frequency of the 198T mutation steadily increased, reaching over 80% from 2015 to 2017. However, from 2017 to 2020, the prevalence of 198T steadily declined, dropping to 70% in 2021, while the 198V mutation showed an inverse trend, rising in frequency and becoming dominant by 2022 (Figure 1B). The 198A mutation frequency remained stable (10%-15%) throughout the study period, with minor fluctuations.

Rescue and identification of recombinant viruses

Following transfection of COS-1 cells with eight plasmids encoding the full H9N2 genome, the supernatants were inoculated into SPF chicken embryos. After 90 hours, HA titers were measured in the allantoic fluid, and recombinant viruses with HA titers ≥2² were considered positive. The HA titers for rWJ/HA_198A and rWJ/HA_198T viruses were 2⁶ and 2⁸, respectively. Sequencing confirmed that the HA genes contained the desired mutations at positions 591–594, while all other gene segments (PB2, PB1, PA, NP, NA, M, and NS) were identical to the parental rWJ57 strain.

Mutations at HA position 198 as a key modulator of HI titers

To investigate how mutations at residue 198 influence the antigenic properties of the virus, HI assays were performed using sera from chickens infected with rWJ57. The rWJ57 virus carrying the 198V mutation showed a two-fold lower HI titer compared to rWJ/HA_198A and rWJ/HA_198T (Figure 1C), suggesting that residue 198 plays a critical role in modulating viral antigenicity.

Mutations at HA position 198 affects the binding affinities to non-neutralizing antibodies

To explore the underlying mechanisms, ELISA assays were conducted to measure the binding affinities of neutralizing and non-neutralizing antibodies. The results demonstrated that the rWJ57 virus had significantly lower affinity for non-neutralizing antibodies compared to rWJ/HA_198A (P < 0.001), while rWJ/HA_198T displayed significantly higher affinity than both rWJ57 and rWJ/HA_198A (P < 0.001) (Figures 2A and 2B). However, no significant differences were observed in the binding affinities of the three viruses to neutralizing antibodies (Figure 2A and 2C).

Mutations at HA position 198 influences the receptor binding activity and specificity.

To assess how receptor binding activity affects HI titers, chicken erythrocytes were treated with α2-3,6,8 neuraminidase to modify sialic acid receptors. Stronger receptor binding was reflected by the virus's ability to agglutinate erythrocytes treated with higher concentrations of neuraminidase. Results showed that rWJ57 exhibited receptor binding activity four-fold higher than rWJ/HA_198A and 16-fold higher than rWJ/HA_198T (Figure 2D). These findings suggest that the mutation at position 198 alters receptor binding, which primarily drives HI titer fluctuations in the H9N2 virus, rather than neutralizing antibodies.

While the RDE assay demonstrates the virus's binding activity for sialic acids (both α2-3 SA and α2-6 SA), it cannot differentiate the virus's affinity for specific sialic acid types. To further explore the effect of mutations at HA position 198 on the affinity for specific sialic acids, we conducted a solid-phase binding assay. The results revealed distinct receptor preferences among H9N2 virus variants with mutations at HA position 198. HA_198V exhibited the highest affinity for 6'SLN, while HA_198A showed strongest binding to 3'SLN (Figure 2E). HA_198T demonstrated intermediate binding to 6'SLN but the lowest affinity for 3'SLN (Figure 2F). These findings indicate that the amino acid at position 198 of the HA protein plays a crucial role in determining receptor binding specificity.

Mutations at HA position 198 alter in HA protein-sialic acid receptor interactions

The α2,6 SA formed five hydrogen bond interactions with the HA protein of the rWJ57 strain (HA_198V protein) at amino acid positions 109, 147, 201, 235, and 236, with a binding energy of -6.1 kcal/mol (Figure 3A). It also formed five hydrogen bond interactions with the rWJ/HA_198T HA protein (HA_198T protein) at positions 109, 147, 148, 233, and 236, with a binding energy of -5.9 kcal/mol (Figure 3B), and four interactions with the rWJ/HA_198A HA protein (HA_198A protein) at positions 145, 147, 235, and 236, with a binding energy of -6.0 kcal/mol (Figure 3C).

The α2,3 SA formed six hydrogen bond interactions with the HA_198V protein at amino acid positions 109, 147, 148, 233, 235, and 236, with a binding energy of -5.6 kcal/mol (Figure 3D). It also formed six interactions with the HA_198T protein at positions 109, 147, 198, 233, 235, and 236, with a binding energy of -5.5 kcal/mol (Figure 3E), and three interactions with the HA_198A protein at positions 109, 145, and 147, with a binding energy of -5.5 kcal/mol (Figure 3F).

These results indicate that alterations at amino acid position 198 of the HA protein directly modify the microstructure of the receptor-binding pocket, leading to changes in the binding sites, interaction forces, torque, and binding energy with sialic acid receptors, ultimately influencing the interactions between the HA protein and the receptors.

MD simulation

To explore the molecular mechanisms by which mutations at residue 198 influence receptor binding affinity and viral tropism, molecular dynamics (MD) simulations were employed to evaluate the structural and dynamic behavior of the HA-sialic acid receptor complexes. The simulation focused on how the HA_198V, HA_198T, and HA_198A mutations alter the interaction between the HA protein and both α2,6 SA and α2,3 SA. The key metrics analyzed were root mean square deviation (RMSD), root mean square fluctuation (RMSF), hydrogen bonding patterns, and binding stability, which provided insights into conformational changes and receptor affinity.

MD simulation demonstrates stable binding stability between HA_198V and α2,6 SA

During the simulation, the HAprotein maintained a stable interaction with the α2,6 SA. The RMSD values for the HA_198V protein and the α2,6 SA stabilized around 7.0 Å and 9.0 Å, respectively, after initial fluctuations (Figure 4A), suggesting that the complex underwent minor adjustments before reaching a stable conformation. RMSF analysis (Figures 4B and 4C) revealed minimal fluctuations in key binding residues, particularly Tyr109, Thr147, Met235, and Gly236, highlighting their crucial role in maintaining the interaction’s stability.

Thr147 played a pivotal role in the binding process, forming hydrogen bonds with a high formation frequency of 118%, supported by secondary water bridge interactions with a frequency of 30% (Figures 4D and 4E). Other residues, including Tyr109, Asn145, Trp161, and Gly236, also contributed through hydrophobic and ionic interactions, further stabilizing the complex. The α2,6 SA itself formed several intramolecular hydrogen bonds, which contributed to its conformational stability (Figure 4E).

These findings suggest that the HA_198V mutation significantly enhances the binding affinity between the HA protein and α2,6 SA. The combination of hydrogen bonds, hydrophobic interactions, and water bridges forms a robust network that stabilizes the HA -α2,6 SA complex.

MD simulation reveals enhanced binding affinity between HA_198T and α2,6 SA

In the MD simulation of the HA_198Tprotein complexed with the α2,6 SA, the RMSD of the HA_198T protein stabilized at approximately 9.0 Å, while the α2,6 SA's RMSD stabilized at around 16 Å after initial fluctuations (Figure 5A). These RMSD values indicate that both the protein and α2,6 SA underwent structural adjustments before reaching a stable conformation.

Further analysis using RMSF revealed that the HA_198T protein exhibited relatively low fluctuations, indicating minimal conformational changes during the binding process (Figure 5B). In contrast, the α2,6 SA displayed higher RMSF values, highlighting its dynamic and flexible nature (Figure 5C). This flexibility likely facilitated binding interactions with the protein, contributing to the formation of a more stable complex.

A total of 30 interactions were identified between the HA_198T protein and the α2,6 SA, with Thr147 playing a critical role, showing a hydrogen bond formation frequency of 92% (Figures 5D and 5E). Other residues, such as Asn145 and Gly233, also contributed significantly to the stability of the complex, reinforcing the hydrogen bonding network. Additionally, the α2,6 SA formed several intramolecular hydrogen bonds, further stabilizing its conformation and strengthening the overall α2,6 SA-protein complex (Figure 5E).

Overall, the HA_198T mutation enhances the binding affinity between the HA protein and the α2,6 SA. The hydrogen bond network, particularly involving Thr147, is crucial in stabilizing the complex. The flexibility of the α2,6 SA, along with its intramolecular interactions, plays a significant role in the stability and dynamic behavior of the α2,6 SA-protein interaction. These findings provide valuable insights into the molecular mechanisms underlying the increased receptor binding affinity due to the HA_198T mutation, with potential implications for future vaccine or antiviral drug development.

MD simulation reveals unstable binding and dissociation between HA_198A and α2,6 SA

During the molecular dynamics simulation, the HA_198A protein exhibited significant fluctuations. The RMSD of HA_198A eventually stabilized at approximately 4.0 Å, whereas the RMSD of the α2,6 SA fluctuated substantially, finally stabilizing near 90 Å (Figure 6A).

While RMSF of the HA_198A protein remained relatively small, the RMSF of the α2,6 SA was notably larger (Figures 6B and 6C), indicating greater movement of the α2,6 SA during the simulation. Despite the formation of 29 interactions, the α2,6 SA eventually moved away from the active binding pocket (Figures 6D and 6E).

These results suggest that the initial binding conformation between the α2,6 SA and the HA_198A protein was unstable, characterized by weak binding affinity. During the simulation, the α2,6 SA ultimately dissociated from the binding pocket, failing to maintain a stable interaction with the HA protein.

MD simulation reveals moderate binding affinity between HA_198V and α2,3 SA after initial instability

During the MD simulation, the RMSD of the HA_198V protein stabilized at at approximately 8.0 Å, while the RMSD of the α2,3 SA, after some fluctuations, settled around 10.0 Å (Figure 7A). The RMSF values for the amino acid residues interacting with the α2,3 SA, highlighted by green lines in Figure 7B, were relatively low, indicating minimal conformational changes in these regions during the binding process. In contrast, the RMSF values for the atoms in the α2,3 SA were relatively high, reflecting substantial fluctuations in the α2,3 SA 's structure (Figure 7C). These findings suggest that the initial binding conformation was unstable, though the system eventually achieved a more stable conformation after significant fluctuations.

The HA_198V protein and the α2,3 SA formed 41 interactions, with Phe408 playing a pivotal role by forming hydrogen bonds with a 31% frequency, highlighting its importance in the binding process. Other interaction sites contributed via hydrogen bonds and water bridges, although these interactions were comparatively weak (Figures 7D and 7E). Additionally, the α2,3 SA formed several intramolecular hydrogen bonds that helped stabilize its conformation during binding (Figure 7E).

Overall, these results indicate that the HA_198V protein and α2,3 SA optimized their binding conformation through the simulation, achieving a more stable structure with moderate binding affinity.

MD simulation confirms weak binding characteristics between HA_198T and α2,3 SA

During the MD simulation, the HA_198T protein exhibited minimal fluctuations, with its RMSD stabilizing at approximately 6 Å. In contrast, the RMSD of the α2,3 SA displayed significant variability before eventually stabilizing at around 30 Å (Figure 8A). The RMSF values for the HA_198T protein were relatively low, indicating limited conformational changes, while the α2,3 SA exhibited much higher RMSF values, suggesting considerable fluctuations in the α2,3 SA's structure throughout the simulation (Figures 8B and 8C). This pattern suggests that the initial conformation of the HA_198T protein and α2,3 SA complex was unstable, with the α2,3 SA gradually moving away from the binding pocket during the course of the simulation.

Although a total of 91 interactions were formed between the protein and α2,3 SA, none of these interactions had frequencies exceeding 30% (Figure 8D), and only those above this threshold are shown in Figures 8D and 8E. The absence of strong, consistent interactions indicates a weak binding affinity between the HA_198T protein and the α2,3 SA. Consequently, the α2,3 SA was unable to maintain stable binding within the protein's receptor pocket, which is consistent with its detachment from the binding site.

MD simulation reveals high-affinity stable binding between HA_198A and α2,3 SA

During the MD simulation, both the HA_198A protein and the α2,3 SA demonstrated high stability. The RMSD of the HA_198A protein stabilized at approximately 6.0 Å, and after initial fluctuations, the RMSD of the α2,3 SA also stabilized at around 6.0 Å (Figure 9A). The RMSF values for the interacting amino acid residues of HA_198A (indicated in green in Figure 9B) and the α2,3 SA atoms showed minimal fluctuations overall (Figures 9B and 9C), indicating that both the protein and the α2,3 SA had relatively stable initial conformations. After some dynamic movements, a more stable binding conformation was formed.

The HA_198A protein and the α2,3 SA formed a total of 33 interactions. Notably, the hydrogen bond and water bridge interaction frequencies of Glu407 reached 159% and 99%, respectively, underscoring the critical role of Glu407 in the binding process (Figures 9D and 9E). Other important residues, including Glu99, Arg275, Ile276, Lys278, and Asp405, also contributed significantly to maintaining the complex's stability (Figures 9D and 9E). Additionally, several intramolecular hydrogen bonds were formed within the α2,3 SA, further stabilizing its binding conformation (Figure 9E).

These results suggest that the α2,3 SA and HA_198A protein optimized their initial binding conformations during molecular dynamics simulations, resulting in a more stable complex with high binding affinity.

In summary, MD simulations provided insights into the molecular interactions and conformational changes governing the binding affinity between the HA protein and sialic acids, which are critical for viral receptor specificity. The binding affinity between the 198th residue of the HA protein and the α2,6 SA decreased in the order of HA_198V, HA_198T, and HA_198A. Conversely, for the α2,3 SA, the binding affinity increased in the order of HA_198T, HA_198V, and HA_198A. Furthermore, when interacting with the α2,6 SA, the 198th residue primarily engaged through the formation of water bridges. In contrast, the 198th residue did not directly participate in interactions with the α2,3 SA, suggesting a differential role of this residue in modulating receptor specificity based on SA type.

The ongoing evolution of the H9N2 AIV has led to a mismatch between epidemic and vaccine strains, contributing significantly to immune failure and periodic outbreaks. Recent years have seen extensive documentation of the diversity and high mutability of amino acids at position 198 on the HA protein, closely associated with the adaptive evolution of H9N2 AIV [11, 31–34]. This study aimed to elucidate the molecular mechanisms by which high-frequency amino acids at position 198 on the HA protein influence receptor binding affinity and specificity, providing new insights into the genetic evolution patterns, antigenic variation mechanisms, and cross-species transmission potential of H9N2 AIV.

Our analysis revealed that between 2013 and 2023, 79.36% of H9N2 AIV strains in China contained the 190T/V mutation, corroborating and extending the findings of Sun et al. [35], who reported an increase from 22.9% before 1999 to 57.2% between 2013 and 2019. This significant increase could explain the rising number of H9N2 avian influenza viruses with enhanced mammalian infection potential, potentially contributing to the observed increase in human infections. The public health implications of this trend are considerable, highlighting the need for enhanced surveillance strategies and the development of more effective vaccines that can address the evolving nature of these viruses. Furthermore, we observed a negative correlation between the prevalence trends of the 198V and 198T mutations. This intriguing finding warrants further investigation, as it may provide crucial insights into the genetic evolution and adaptation strategies of H9N2 AIV. The mechanisms driving this negative correlation could be related to the differential effects of these mutations on viral fitness in various host environments or to complex interactions with other viral proteins. Understanding these dynamics could be key to predicting future evolutionary trajectories of H9N2 AIV and improving our preparedness for potential outbreaks.

Our findings demonstrate that mutations at HA position 198 significantly influence receptor binding preferences and affinities, subsequently impacting HI titers and contributing to antigenic drift. Specifically, HA_198V and HA_198T variants exhibit higher affinity for 6'SLN, a receptor more prevalent in the upper respiratory tract of mammals [36], compared to HA_198A. Conversely, the HA_198A variant shows a stronger preference for 3'SLN. These differential binding patterns, consistently demonstrated in both RDE and solid-phase binding assays, align with previous reports of HA_198V/T mutations promoting mammalian adaptation [31]. Moreover, our study uncovered novel differences in non-neutralizing antibody binding among HA_198V, HA_198T, and HA_198A variants, suggesting a complex interplay between receptor binding and antibody recognition. These findings have significant implications for vaccine design and efficacy, particularly in light of the ongoing challenge of antigenic drift in H9N2 AIV. The altered antibody binding profiles associated with different HA₁₉₈ variants underscore the need for vaccine strategies that can induce broadly neutralizing antibodies or that can be quickly updated to match circulating strains.

To further elucidate the mechanisms underlying changes in receptor-binding strength and specificity caused by HA198 amino acid mutations, we employed molecular docking and molecular dynamics simulations. These computational approaches revealed that mutations such as HA_198V, HA_198T, and HA_198A induce conformational changes in the receptor-binding pocket, significantly altering the interactions between the H9N2 virus and both α2,6 SA and α2,3 SA receptors. Specifically, during MD simulations, the HA_198V and HA_198T mutations enhanced binding stability with 6'SLN, while HA_198A formed a more stable complex with 3'SLN. We also identified key residues, such as Thr147, Glu199, and Tyr109, that play crucial roles in receptor binding by forming stable hydrogen bonds and water bridges with sialic acids. Notably, residue 198 primarily interacts with α2,6 SA through water bridges, while showing minimal direct involvement in binding with α2,3 SA. These molecular-level changes could translate to broader phenotypic effects, such as altered transmissibility or virulence in different host species. For instance, the enhanced stability of HA_198V and HA_198T with 6'SLN might contribute to more efficient viral entry and replication in mammalian cells, potentially increasing the risk of zoonotic transmission. Conversely, the preference of HA_198A for 3'SLN might optimize the virus for replication in avian hosts. Understanding these structure-function relationships at the molecular level provides valuable insights into the mechanisms of host adaptation and could inform the development of targeted antiviral strategies.

This study underscores the advantages of molecular docking and MD simulations in analyzing virus-receptor interactions. Compared to traditional biological assays, these computational techniques offer more direct visualization of the dynamic interactions between viral proteins and receptors, revealing molecular changes and their biological significance. These tools provide a powerful means to study the adaptive evolution and antigenic drift of the H9N2 virus, allowing for rapid assessment of the potential impact of newly identified HA mutations on receptor binding. However, it's important to acknowledge the limitations of these computational approaches. While they provide valuable insights into molecular interactions, they may not fully capture the complexity of the cellular environment or the potential influence of other viral and host factors. Therefore, these computational studies should be viewed as complementary to experimental approaches. Integrating computational predictions with in vitro and in vivo experimental validation can provide a more comprehensive understanding of virus-host interactions and enhance the reliability of our findings.

In conclusion, this study combined experimental data with computational simulations to elucidate how mutations at residue 198 of the HA protein regulate receptor-binding affinity, specificity, and antibody recognition, driving viral antigenic drift and host adaptation. These findings offer valuable insights into the underlying molecular mechanisms of H9N2 AIV evolution and adaptation.

HA_198V, an alanine residue at position 198 in HA; HA_198V; a valine residue at position 198 in HA; HA_198T; a threonine residue at position 198 in HA; TPCK- trypsin: tosyl-phenylalanine chloromethyl-ketone treated typsin.

Authors’ contributions

RZ and SC participated in the design of the study. RZ, JW and RC performed the experiments. MZ and ZW provide intellectual input. RZ and LW analysed the data. RZ wrote the manuscript. LZ and SZ reviewed and prepared the version of the manuscript. All authors read and approved the final manuscript.

Funding

This research was supported by Jiangsu Agri-Animal Husbandry Vocational College Natural Science Project (NSF2022CB15), Jiangsu University Blue Project [(2023) NO. 27)], Enterprise Practice Training Program for Young Teachers in Higher Vocational Colleges of Jiangsu Province (2023QYSJ037), and Taizhou 311 High-Level Talent Training Project [(2024) NO. 1)].

Conflicts of Interest

All the authors have declared that no conflicts of interest exist regarding to the publication of the data in this manuscript.

Sun Y, Liu J (2015) H9N2 influenza virus in China: A cause of concern. Protein Cell 6:18–25. https://doi.org/10.1007/s13238-014-0111-7
Kishida N, Sakoda Y, Eto M, et al (2004) Co-infection of staphylococcus aureus or haemophilus paragallinarum exacerbates H9N2 influenza a virus infection in chickens. Arch Virol 149:2095–2104. https://doi.org/10.1007/s00705-004-0372-1
Wang J, Li Y, Yin Y (2018) Respiratory phagocytes are implicated in enhanced colibacillosis in chickens co-infected with influenza virus H9N2 and escherichia coli. Br Poult Sci 59:160–165. https://doi.org/10.1080/00071668.2017.1406061
Wang X, Wang H, Zhang S, et al (2023) The role of transforming growth factor beta-1 protein in escherichia coli secondary infection induced by H9N2 avian influenza virus in chickens. Microb Pathog 175:105983. https://doi.org/10.1016/j.micpath.2023.105983
Pu J, Wang S, Yin Y, et al (2015) Evolution of the H9N2 influenza genotype that facilitated the genesis of the novel H7N9 virus. Proc Natl Acad Sci U S A 112:548–553. https://doi.org/10.1073/pnas.1422456112
Ding S, Zhou J, Xiong J, et al (2024) Continued evolution of H10N3 influenza virus with adaptive mutations poses an increased threat to mammals. Virol Sin 39:546–555. https://doi.org/10.1016/j.virs.2024.06.005
Li B, Su G, Xiao C, et al (2022) The PB2 co-adaptation of H10N8 avian influenza virus increases the pathogenicity to chickens and mice. Transbound Emerg Dis 69:1794–1803. https://doi.org/10.1111/tbed.14157
Q M, S Z, S L, et al (2024) Emergence of novel reassortant H3N3 avian influenza viruses with increased pathogenicity in chickens in 2023. Emerging microbes & infections 13:. https://doi.org/10.1080/22221751.2023.2287683
Li P, Niu M, Li Y, et al (2022) Human infection with H3N8 avian influenza virus: A novel H9N2-original reassortment virus. J Infect 85:e187–e189. https://doi.org/10.1016/j.jinf.2022.08.033
Koel BF, Burke DF, Bestebroer TM, et al (2013) Substitutions near the receptor binding site determine major antigenic change during influenza virus evolution. Science 342:976–979. https://doi.org/10.1126/science.1244730
Zhu R, Xu D, Yang X, et al (2018) Genetic and biological characterization of H9N2 avian influenza viruses isolated in China from 2011 to 2014. PLoS One 13:e0199260. https://doi.org/10.1371/journal.pone.0199260
Su H, Zhao Y, Zheng L, et al (2020) Effect of the selection pressure of vaccine antibodies on evolution of H9N2 avian influenza virus in chickens. AMB Express 10:98. https://doi.org/10.1186/s13568-020-01036-0
Jin H, Wang W, Yang X, et al (2018) Evolution of H9N2 avian influenza virus in embryonated chicken eggs with or without homologous vaccine antibodies. BMC Vet Res 14:71. https://doi.org/10.1186/s12917-018-1391-6
Zhu R, Xu S, Sun W, et al (2022) HA gene amino acid mutations contribute to antigenic variation and immune escape of H9N2 influenza virus. Vet Res 53:43. https://doi.org/10.1186/s13567-022-01058-5
Hensley SE, Das SR, Bailey AL, et al (2009) Hemagglutinin receptor binding avidity drives influenza a virus antigenic drift. Science 326:734–736. https://doi.org/10.1126/science.1178258
Sealy JE, Yaqub T, Peacock TP, et al (2018) Association of increased receptor-binding avidity of influenza a(H9N2) viruses with escape from antibody-based immunity and enhanced zoonotic potential. Emerg Infect Dis 25:63–72. https://doi.org/10.3201/eid2501.180616
Peacock TP, Sealy JE, Harvey WT, et al (2021) Genetic determinants of receptor-binding preference and zoonotic potential of H9N2 avian influenza viruses. J Virol 95:e01651-20. https://doi.org/10.1128/JVI.01651-20
Koehler M, Delguste M, Sieben C, et al (2020) Initial step of virus entry: Virion binding to cell-surface glycans. Annu Rev Virol 7:143–165. https://doi.org/10.1146/annurev-virology-122019-070025
Overeem NJ, van der Vries E, Huskens J (2021) A dynamic, supramolecular view on the multivalent interaction between influenza virus and host cell. Small 17:2007214. https://doi.org/10.1002/smll.202007214
Zhao C, Pu J (2022) Influence of host sialic acid receptors structure on the host specificity of influenza viruses. Viruses 14:2141. https://doi.org/10.3390/v14102141
Vora J, Athar M, Sinha S, et al (2020) Binding insight of anti-HIV phytocompounds with prime targets of HIV: A molecular dynamics simulation analysis. Curr HIV Res 18:132–141. https://doi.org/10.2174/1570162X18666200129112509
Durrant JD, McCammon JA (2011) Molecular dynamics simulations and drug discovery. BMC Biol 9:71. https://doi.org/10.1186/1741-7007-9-71
Hao X, Hu J, Wang J, et al (2016) Reassortant H5N1 avian influenza viruses containing PA or NP gene from an H9N2 virus significantly increase the pathogenicity in mice. Vet Microbiol 192:95–101. https://doi.org/10.1016/j.vetmic.2016.07.002
Peng Q, Zhu R, Wang X, et al (2019) Impact of the variations in potential glycosylation sites of the hemagglutinin of H9N2 influenza virus. Virus Genes 55:182–190. https://doi.org/10.1007/s11262-018-1623-7
Chambers BS, Parkhouse K, Ross TM, et al (2015) Identification of hemagglutinin residues responsible for H3N2 antigenic drift during the 2014–2015 influenza season. Cell Rep 12:1–6. https://doi.org/10.1016/j.celrep.2015.06.005
Zhang N, Quan K, Chen Z, et al (2023) The emergence of new antigen branches of H9N2 avian influenza virus in China due to antigenic drift on hemagglutinin through antibody escape at immunodominant sites. Emerg Microbes Infect 12:2246582. https://doi.org/10.1080/22221751.2023.2246582
Gao R, Gu M, Liu K, et al (2018) T160A mutation-induced deglycosylation at site 158 in hemagglutinin is a critical determinant of the dual receptor binding properties of clade 2.3.4.4 H5NX subtype avian influenza viruses. Veterinary Microbiology 217:158–166. https://doi.org/10.1016/j.vetmic.2018.03.018
O’Boyle NM, Banck M, James CA, et al (2011) Open babel: An open chemical toolbox. J Cheminform 3:33. https://doi.org/10.1186/1758-2946-3-33
Eberhardt J, Santos-Martins D, Tillack AF, Forli S (2021) AutoDock vina 1.2.0: New docking methods, expanded force field, and python bindings. J Chem Inf Model 61:3891–3898. https://doi.org/10.1021/acs.jcim.1c00203
Trott O, Olson AJ (2010) AutoDock vina: Improving the speed and accuracy of docking with a new scoring function, efficient optimization and multithreading. J Comput Chem 31:455–461. https://doi.org/10.1002/jcc.21334
Amin F, Mukhtar N, Ali M, et al (2024) Mapping genetic markers associated with antigenicity and host range in H9N2 influenza a viruses infecting poultry in pakistan. Avian Dis 68:43–51. https://doi.org/10.1637/aviandiseases-D-23-00029
Gao X, Wang N, Chen Y, et al (2021) Sequence characteristics and phylogenetic analysis of H9N2 subtype avian influenza a viruses detected from poultry and the environment in China, 2018. PeerJ 9:e12512. https://doi.org/10.7717/peerj.12512
Yang F, Xiao Y, Liu F, et al (2021) Molecular characterization and antigenic analysis of reassortant H9N2 subtype avian influenza viruses in eastern China in 2016. Virus Research 306:198577. https://doi.org/10.1016/j.virusres.2021.198577
Xia J, Li Y-X, Dong M-Y, et al (2023) Evolution of prevalent H9N2 subtype of avian influenza virus during 2019 to 2022 for the development of a control strategy in China. Poult Sci 102:102957. https://doi.org/10.1016/j.psj.2023.102957
Sun X, Belser JA, Maines TR (2020) Adaptation of H9N2 influenza viruses to mammalian hosts: A review of molecular markers. Viruses 12:541. https://doi.org/10.3390/v12050541
Sriwilaijaroen N, Kondo S, Yagi H, et al (2011) N-glycans from porcine trachea and lung: Predominant NeuAcα2-6Gal could be a selective pressure for influenza variants in favor of human-type receptor. PLoS One 6:e16302. https://doi.org/10.1371/journal.pone.0016302

Download PDF

Version 1

posted

You are reading this latest preprint version

HA198 Mutations in H9N2 Avian Influenza: Molecular Dynamics Insights into Receptor Binding

Status:

Version 1

Abstract

Figures

Introduction

Materials and methods

Results

Discussion

Abbreviations

Declarations

References

Status:

Version 1