The recent occurrence of 2009 influenza A (H1N1) pandemic as well as others has raised concern of a far more dangerous outcome should this virus becomes resistant to current drug therapies. The number of clinical cases that are resistant to oseltamivir (Tamiflu®) is larger than the limited number of neuraminidase (NA) mutations (H275Y, N295S, and I223R) that have been identified at the active site and that are associated to oseltamivir resistance. In this study, we have performed a comparative analysis between a set of NAs that have the most representative mutations located outside the active site. The recently crystallized NA-oseltamivir complex (PDB ID: 3NSS) was used as a wild-type structure. After selecting the target NA sequences, their three-dimensional (3D) structure was built using 3NSS as a template by homology modeling. The 3D NA models were refined by molecular dynamics (MD) simulations. The refined models were used to perform a docking study, using oseltamivir as a ligand. Furthermore, the docking results were refined by free-energy analysis using the MM-PBSA method. The analysis of the MD simulation results showed that the NA models reached convergence during the first 10 ns. Visual inspection and structural measures showed that the mutated NA active sites show structural variations. The docking and MM-PBSA results from the complexes showed different binding modes and free energy values. These results suggest that distant mutations located outside the active site of NA affect its structure and could be considered to be a new source of resistance to oseltamivir, which agrees with reports in the clinical literature.
Copyright © 2012 Wiley Periodicals, Inc.