Insights into susceptibility of antiviral drugs against the E119G mutant of 2009 influenza A (H1N1) neuraminidase by molecular dynamics simulations and free energy calculations

Neuraminidase inhibitors (NAIs) play vital roles in controlling human influenza epidemics and pandemics. However, the emergence of new human influenza virus mutant strains resistant to existing antiviral drugs has been becoming a major challenge. Therefore, it is critical to uncover the mechanisms of drug resistance and seek alternative treatments to combat drug resistance. In this study, molecular dynamics (MD) simulations and Molecular Mechanics/Generalized Born Surface Area (MM/GBSA) were applied to investigate the different sensitivities of oseltamivir (OTV), zanamivir (ZNV), and peramivir (PRV) against the E119G mutant of 2009 A/H1N1 neuraminidase. The predicted binding free energies indicate that the E119G mutation in NA confers resistance to all of the three studied inhibitors. The ordering of the level of drug resistance predicted by the binding free energies for the three inhibitors is ZNV > PRV > OTV, which agrees well with the experimental data. Drug resistance arises primarily from the unfavorable shifts of the polar interactions between NA and the inhibitors. It comes as a surprise that the mutation of Glu119 that can form strong H-bonds with the inhibitors in the wild-type protein does not have direct impact on the binding affinities of both OTV and PRV due to the regulation of the strong unfavorable polar desolvation energies. The indirectly conformational variations of the inhibitors, which caused by the E119G mutation, are responsible for the loss of the binding free energies. However, for ZNV, the E119G mutation has both direct and indirect influences on the drug binding. The structural and quantitative viewpoint obtained from this study provides valuable information for the rational design of novel and effective drugs to combat drug resistance.