Molecular modeling of swine influenza A/H1N1, Spanish H1N1, and avian H5N1 flu N1 neuraminidases bound to Tamiflu and Relenza

A molecular model of the swine influenza A/H1N1 ( also called H1N1pdm) type-I neuraminidase was built using the pathogenic avian H5N1 type-I neuraminidase as a basis, due to the higher sequence identity between A/H1N1 and H5N1 (91.47%) compared to Spanish H1N1 (88.37%) neuraminidase. All-atom molecular dynamics (MD) simulations of all three neuraminidases were performed, either as apo-structures or with commercial antiviral drugs Tamiflu or Relenza separately bound; the simulations allowed for the identification of both conserved and unique drug-protein interactions across all three proteins. Specifically, conserved networks of hydrogen bonds stabilizing the drugs in the sialic acid binding site of the simulated neuraminidases are analyzed, providing insight into how disruption due to mutations may lead to increased drug resistance. In addition, a possible mechanism through which the residue 294 mutation acquires drug resistance is proposed by mapping the mutation site onto an electrostatic pathway which may play a role in controlling drug access to the binding pocket of neuraminidase, establishing a starting point for further investigations of neuraminidase drug resistance.