Molecular dynamics simulation of the effects of single (S221P) and double (S221P and K216E) mutations in the hemagglutinin protein of influenza A H5N1 virus: a study on host receptor specificity

Avian influenza viruses of subtype H5N1 circulating in animals continue to pose threats to human health. The binding preference of the viral surface protein hemagglutinin (HA) to sialosaccharides of receptors is an important area for understanding mutations in the receptor binding site that could be the cause for avian-to-human transmission. In the present work, we studied the effect of two receptor binding site mutations, S221P singly and in combination with another mutation K216E in the HA protein of influenza A H5N1 viruses. Docking of sialic acid ligands corresponding to both avian and human receptors and molecular dynamics simulations of the complexes for wild and mutant strains of H5N1 viruses were carried out. The H5N1 strain possessing the S221P mutation indicated decreased binding to α2,3-linked sialic acids (avian receptor, SAα2,3Gal) when compared to the binding of the wild-type strain that did not possess the HA-221 mutation. The binding to α2,6-linked sialic acids (human receptor, SAα2,6Gal) was found to be comparable, indicating that the mutant strain shows limited dual receptor specificity. On the other hand, the S221P mutation in synergism with the K216E mutation in the binding site, resulted in increased binding affinity for SAα2,6Gal when compared to SAα2,3Gal, indicative of enhanced binding to human receptors. The in-depth study of the molecular interactions in the docked complexes could explain how co-occurring mutations in the HA viral protein can aid in providing fitness advantage to the virus, in the context of host receptor specificity in emerging variants of H5N1 influenza viruses.