In silico identification of the potential drug resistance sites over 2009 influenza A (H1N1) virus neuraminidase

The outbreak and high speed global spread of the new strain of influenza A (H1N1) virus in 2009 poses a serious threat to the general population and governments. At present, the most effective drugs for the treatment of 2009 influenza A (H1N1) virus are neuraminidase inhibitors: mainly oseltamivir and zanamivir. The use of these two inhibitors will undoubtedly increase, and therefore it is more likely that drug-resistant influenza strains will arise. The identification of the potential resistance sites for these drugs in advance and the understanding of corresponding molecular basis to cause drug resistance are no doubt very important to fight against the new resistant influenza strains. In this study, first, the complexes of neuraminidase with the substrate sialic acid and two inhibitors oseltamivir and zanamivir were obtained by fitting them to the 3D structure of 2009 influenza A (H1N1) neuraminidase obtained by homology modeling. By using these complexes as the initial structures, molecular dynamics simulation and molecular mechanics generalized Born surface area (MM-GBSA) calculations were performed to identify the residues with significant contribution to the binding of substrate and inhibitors. By analyzing the difference of interaction profiles of substrate and inhibitors, the potential drug resistance sites for two inhibitors were identified. Parts of the identified sites have been verified to confer resistance to oseltamivir and zanamivir for influenza virus of the past flu epidemic. The identified potential resistance sites in this study will be useful for the development of new effective drugs against the drug resistance and avoid the situation of having no effective drugs to treat new mutant influenza strains.