Characterizing receptor flexibility to predict mutations that lead to human adaptation of influenza hemagglutinin

Abstract

A key step in the emergence of human pandemic influenza strains has been a switch in the binding preference of the viral glycoprotein hemagglutinin (HA) from avian to human sialic acid receptors (SAs). The conformation of the bound SA varies substantially with HA sequence, and crystallographic evidence suggests that the bound SA is flexible, so it is difficult to predict from crystal structures which mutations are responsible for the change in HA binding preference. We performed molecular dynamics (MD) simulations of SA analogs binding to various HAs, and observed a dynamic equilibrium among structurally diverse receptor conformations, including novel conformations that have not been experimentally observed. Using one such novel conformation, we predicted--and subsequently confirmed with microscale thermophoresis experiments--a set of mutations that substantially increased an HAs affinity for a human SA analog. This prediction could not have been inferred from existing crystal structures, suggesting that MD-generated HA-SA conformational ensembles could help researchers predict human-adaptive mutations, aiding in the surveillance of emerging pandemic threats.