[Source: PLoS ONE, full text: (LINK). Abstract, edited.]
-------
Exploring the Mechanism of Zanamivir Resistance in a Neuraminidase Mutant: A Molecular Dynamics Study
Nanyu Han<SUP>1</SUP>, Xuewei Liu<SUP>2</SUP>, Yuguang Mu<SUP>1</SUP><SUP>*</SUP>
<SUP></SUP>
1 School of Biological Sciences, Nanyang Technological University, Singapore, Singapore, 2 School of Physical and Mathematical Sciences, Nanyang Technological University, Singapore, Singapore
Abstract
It is critical to understand the molecular basis of the drug resistance of influenza viruses to efficiently treat this infectious disease. Recently, H1N1 strains of influenza A carrying a mutation of Q136K in neuraminidase were found. The new strain showed a strong Zanamivir neutralization effect. In this study, normal molecular dynamics simulations and metadynamics simulations were employed to explore the mechanism of Zanamivir resistance. The wild-type neuraminidase contained a 3<SUB>10</SUB> helix before the 150 loop, and there was interaction between the 150 and 430 loops. However, the helix and the interaction between the two loops were disturbed in the mutant protein due to interaction between K136 and nearby residues. Hydrogen-bond network analysis showed weakened interaction between the Zanamivir drug and E276/D151 on account of the electrostatic interaction between K136 and D151. Metadynamics simulations showed that the free energy landscape was different in the mutant than in the wild-type neuraminidase. Conformation with the global minimum of free energy for the mutant protein was different from the wild-type conformation. While the drug fit completely into the active site of the wild-type neuraminidase, it did not match the active site of the mutant variant. This study indicates that the altered hydrogen-bond network and the deformation of the 150 loop are the key factors in development of Zanamivir resistance. Furthermore, the Q136K mutation has a variable effect on conformation of different N1 variants, with conformation of the 1918 N1 variant being more profoundly affected than that of the other N1 variants studied in this paper. This observation warrants further experimental investigation.
Citation: Han N, Liu X, Mu Y (2012) Exploring the Mechanism of Zanamivir Resistance in a Neuraminidase Mutant: A Molecular Dynamics Study. PLoS ONE 7(9): e44057. doi:10.1371/journal.pone.0044057
Editor: Yury E. Khudyakov, Centers for Disease Control and Prevention, United States of America
Received: December 23, 2011; Accepted: July 30, 2012; Published: September 6, 2012
Copyright: ? 2012 Han et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Funding: The support of research grants, NMRC/H1N1R/001/2009 and the IDA Cloud Computing Call for Project Proposals 2012 is gratefully acknowledged. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
* E-mail: ygmu@ntu.edu.sg
-Nanyu Han<SUP>1</SUP>, Xuewei Liu<SUP>2</SUP>, Yuguang Mu<SUP>1</SUP><SUP>*</SUP>
<SUP></SUP>
1 School of Biological Sciences, Nanyang Technological University, Singapore, Singapore, 2 School of Physical and Mathematical Sciences, Nanyang Technological University, Singapore, Singapore
Abstract
It is critical to understand the molecular basis of the drug resistance of influenza viruses to efficiently treat this infectious disease. Recently, H1N1 strains of influenza A carrying a mutation of Q136K in neuraminidase were found. The new strain showed a strong Zanamivir neutralization effect. In this study, normal molecular dynamics simulations and metadynamics simulations were employed to explore the mechanism of Zanamivir resistance. The wild-type neuraminidase contained a 3<SUB>10</SUB> helix before the 150 loop, and there was interaction between the 150 and 430 loops. However, the helix and the interaction between the two loops were disturbed in the mutant protein due to interaction between K136 and nearby residues. Hydrogen-bond network analysis showed weakened interaction between the Zanamivir drug and E276/D151 on account of the electrostatic interaction between K136 and D151. Metadynamics simulations showed that the free energy landscape was different in the mutant than in the wild-type neuraminidase. Conformation with the global minimum of free energy for the mutant protein was different from the wild-type conformation. While the drug fit completely into the active site of the wild-type neuraminidase, it did not match the active site of the mutant variant. This study indicates that the altered hydrogen-bond network and the deformation of the 150 loop are the key factors in development of Zanamivir resistance. Furthermore, the Q136K mutation has a variable effect on conformation of different N1 variants, with conformation of the 1918 N1 variant being more profoundly affected than that of the other N1 variants studied in this paper. This observation warrants further experimental investigation.
Citation: Han N, Liu X, Mu Y (2012) Exploring the Mechanism of Zanamivir Resistance in a Neuraminidase Mutant: A Molecular Dynamics Study. PLoS ONE 7(9): e44057. doi:10.1371/journal.pone.0044057
Editor: Yury E. Khudyakov, Centers for Disease Control and Prevention, United States of America
Received: December 23, 2011; Accepted: July 30, 2012; Published: September 6, 2012
Copyright: ? 2012 Han et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Funding: The support of research grants, NMRC/H1N1R/001/2009 and the IDA Cloud Computing Call for Project Proposals 2012 is gratefully acknowledged. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
* E-mail: ygmu@ntu.edu.sg
-------