Distinct Glycan Topology for Avian and Human Sialopentasaccharide Receptor Analogues upon Binding Different Hemagglutinins: A Molecular Dynamics Perspective

• Published in: Journal of molecular biology 2009-03, Vol.387 (2), p.465-491
• Main Authors: Xu, Dong, Newhouse, E. Irene, Amaro, Rommie E., Pao, Hsing C., Cheng, Lily S., Markwick, Phineus R.L., McCammon, J. Andrew, Li, Wilfred W., Arzberger, Peter W.
• Format: Article
• Language: English
• Subjects: Animals; avian influenza; Avian influenza virus; binding specificity; Birds - metabolism; Computer Simulation; Entropy; Galactose - metabolism; glycan topology; hemagglutinin; Hemagglutinins - metabolism; Humans; Hydrogen Bonding; Models, Molecular; molecular dynamics; N-Acetylneuraminic Acid - metabolism; Polysaccharides - chemistry; Protein Binding; Protein Structure, Secondary; Receptors, Cell Surface - chemistry; Receptors, Cell Surface - metabolism; Sequence Homology, Amino Acid; Solutions
• ISSN: 0022-2836; 1089-8638
• DOI: 10.1016/j.jmb.2009.01.040

Hemagglutinin (HA) binds to sialylated glycans exposed on the host cell surface in the initial stage of avian influenza virus infection. It has been previously hypothesized that glycan topology plays a critical role in the human adaptation of avian flu viruses, such as the potentially pandemic H5N1. Comparative molecular dynamics studies are complementary to experimental techniques, including glycan microarray, to understand the mechanism of species-specificity switch better. The examined systems comprise explicitly solvated trimeric forms of avian H3, H5, and swine H9 in complex with avian and human glycan receptor analogues—LSTa (α-2,3-linked lactoseries tetrasaccharide a) and LSTc (α-2,6-linked lactoseries tetrasaccharide c), respectively. The glycans adopted distinct topological profiles with inducible torsional angles when bound to different HAs. The corresponding receptor binding domain amino acid contact profiles were also distinct. Avian H5 was able to accommodate LSTc in a tightly “folded umbrella”-like topology through interactions with all five sugar residues. After considering conformational entropy, the relative binding free-energy changes, calculated using the molecular mechanics-generalized Born surface area technique, were in agreement with previous experimental findings and provided insights on electrostatic, van der Waals, desolvation, and entropic contributions to HA–glycan interactions. The topology profile and the relative abundance of free glycan receptors may influence receptor binding kinetics. Glycan composition and topological changes upon binding different HAs may be important determinants in species-specificity switch.