Theoretical Studies on the Susceptibility of Oseltamivir against Variants of 2009 A/H1N1 Influenza Neuraminidase

The outbreak and high speed global spread of the new strain of influenza A/H1N1 virus in 2009 posed a serious threat to global health. It is more likely that drug-resistant influenza strains will arise after the extensive use of anti-influenza drugs. Consequently, the identification of the potential resistant sites for drugs in advance and the understanding of the corresponding molecular mechanisms that cause drug resistance are quite important in the design of new drug candidates with better potency to combat drug resistance. Here, we performed molecular simulations to evaluate the potency of oseltamivir to combat drug resistance caused by the mutations in 2009 A/H1N1 neuraminidase (NA). We examined three representative drug-resistant mutations in NA, consisting of H274Y, N294S, and Y252H. First, a theoretical structure of A/H1N1 NA in complex with oseltamivir was constructed using homology modeling. Then, molecular dynamics (MD) simulations, molecular mechanics/Poisson-Boltzmann surface area (MM/PBSA) calculations, and MM/GBSA free energy decomposition were used to characterize the binding of oseltamivir with the wild type (WT) and three mutated NAs. Our predictions show that N294S and H274Y, two popular drug-resistant mutations in different variants of NA, still cause significant resistance to oseltamivir. However, the Y252H mutation does not impair the interactions between oseltamivir and A/H1N1 NA. An examination of individual energy components shows that the loss of polar interactions is the key source for the resistance of the studied mutations to oseltamivir. Moreover, free energy decomposition analysis and structural analysis reveal that the N294S or H274Y mutation triggers the large-scale conformational changes of the binding pocket and then impairs the affinity of oseltamivir. We expect that our results will be useful for the rational design of NA inhibitors with high potency against drug-resistant A/H1N1 mutants.

[1]  Rommie E. Amaro,et al.  Ensemble-Based Virtual Screening Reveals Potential Novel Antiviral Compounds for Avian Influenza Neuraminidase , 2008, Journal of medicinal chemistry.

[2]  Zhi-Yuan Su,et al.  Homology modeling, docking, and molecular dynamics reveal HR1039 as a potent inhibitor of 2009 A(H1N1) influenza neuraminidase. , 2010, Biophysical chemistry.

[3]  R. Zwanzig High‐Temperature Equation of State by a Perturbation Method. I. Nonpolar Gases , 1954 .

[4]  Jan H. Jensen,et al.  Very fast empirical prediction and rationalization of protein pKa values , 2005, Proteins.

[5]  T. Darden,et al.  A smooth particle mesh Ewald method , 1995 .

[6]  Wei Wang,et al.  Characterization of domain-peptide interaction interface: prediction of SH3 domain-mediated protein-protein interaction network in yeast by generic structure-based models. , 2012, Journal of proteome research.

[7]  Junmei Wang,et al.  Development and testing of a general amber force field , 2004, J. Comput. Chem..

[8]  D. Case,et al.  Characterization of domain-peptide interaction interface: a case study on the amphiphysin-1 SH3 domain. , 2008, Journal of molecular biology.

[9]  P. Kollman,et al.  Calculating structures and free energies of complex molecules: combining molecular mechanics and continuum models. , 2000, Accounts of chemical research.

[10]  Hiroshi Chuman,et al.  Correlation Analyses on Binding Affinity of Sialic Acid Analogues with Influenza Virus Neuraminidase-1 Using ab Initio MO Calculations on Their Complex Structures , 2010, J. Chem. Inf. Model..

[11]  Jie J. Zheng,et al.  Computational studies of H5N1 influenza virus resistance to oseltamivir , 2009, Protein science : a publication of the Protein Society.

[12]  Janet Daly,et al.  Mutations conferring zanamivir resistance in human influenza virus N2 neuraminidases compromise virus fitness and are not stably maintained in vitro. , 2006, The Journal of antimicrobial chemotherapy.

[13]  Hideo Goto,et al.  Avian flu: Isolation of drug-resistant H5N1 virus , 2005, Nature.

[14]  G. Klebe,et al.  Approaches to the Description and Prediction of the Binding Affinity of Small-Molecule Ligands to Macromolecular Receptors , 2002 .

[15]  Yi Guan,et al.  Oseltamivir resistance during treatment of influenza A (H5N1) infection. , 2005, The New England journal of medicine.

[16]  Huanxiang Liu,et al.  Prediction of zanamivir efficiency over the possible 2009 Influenza A (H1N1) mutants by multiple molecular dynamics simulations and free energy calculations , 2011, Journal of molecular modeling.

[17]  R. Dwek,et al.  Structural Characterization of the 1918 Influenza Virus H1N1 Neuraminidase , 2008, Journal of Virology.

[18]  Urban Bren,et al.  Do all pieces make a whole? Thiele cumulants and the free energy decomposition , 2007 .

[19]  J. Andrew McCammon,et al.  Independent-Trajectories Thermodynamic-Integration Free-Energy Changes for Biomolecular Systems: Determinants of H5N1 Avian Influenza Virus Neuraminidase Inhibition by Peramivir , 2009, Journal of chemical theory and computation.

[20]  Mark von Itzstein,et al.  The war against influenza: discovery and development of sialidase inhibitors , 2007, Nature Reviews Drug Discovery.

[21]  Susana López,et al.  Molecular anatomy of 2009 influenza virus A (H1N1). , 2009, Archives of medical research.

[22]  Ricky Chachra,et al.  Origins of Resistance Conferred by the R292K Neuraminidase Mutation via Molecular Dynamics and Free Energy Calculations. , 2008, Journal of chemical theory and computation.

[23]  W. L. Jorgensen,et al.  Comparison of simple potential functions for simulating liquid water , 1983 .

[24]  Tingjun Hou,et al.  Assessing the Performance of the MM/PBSA and MM/GBSA Methods. 1. The Accuracy of Binding Free Energy Calculations Based on Molecular Dynamics Simulations , 2011, J. Chem. Inf. Model..

[25]  Xiaowu Chen,et al.  Discovery and development of GS 4104 (oseltamivir): an orally active influenza neuraminidase inhibitor. , 2000, Current medicinal chemistry.

[26]  Thanyada Rungrotmongkol,et al.  Molecular prediction of oseltamivir efficiency against probable influenza A (H1N1-2009) mutants: molecular modeling approach , 2010, Amino Acids.

[27]  Irwin D Kuntz,et al.  Investigation of neuraminidase-substrate recognition using molecular dynamics and free energy calculations. , 2003, Journal of medicinal chemistry.

[28]  P A Kollman,et al.  Absolute and relative binding free energy calculations of the interaction of biotin and its analogs with streptavidin using molecular dynamics/free energy perturbation approaches , 1993, Proteins.

[29]  Ron A M Fouchier,et al.  Antigenic and Genetic Characteristics of Swine-Origin 2009 A(H1N1) Influenza Viruses Circulating in Humans , 2009, Science.

[30]  Yoshihiro Kawaoka,et al.  oseltamivir: descriptive study , 2022 .

[31]  Vasiliy P. Mishin,et al.  Susceptibilities of Antiviral-Resistant Influenza Viruses to Novel Neuraminidase Inhibitors , 2005, Antimicrobial Agents and Chemotherapy.

[32]  J. Andrew McCammon,et al.  Role of Secondary Sialic Acid Binding Sites in Influenza N1 Neuraminidase , 2010, Journal of the American Chemical Society.

[33]  J. Kirkwood Statistical Mechanics of Fluid Mixtures , 1935 .

[34]  Jan Balzarini,et al.  The molecular basis of resilience to the effect of the Lys103Asn mutation in non-nucleoside HIV-1 reverse transcriptase inhibitors studied by targeted molecular dynamics simulations. , 2005, Journal of the American Chemical Society.

[35]  Parimal Kar,et al.  Mutation-induced loop opening and energetics for binding of tamiflu to influenza N8 neuraminidase. , 2012, The journal of physical chemistry. B.

[36]  Pascal Bonnet,et al.  Molecular dynamics and free energy analysis of neuraminidase–ligand interactions , 2004, Protein science : a publication of the Protein Society.

[37]  Yoshihiro Kawaoka,et al.  Amino Acid Residues Contributing to the Substrate Specificity of the Influenza A Virus Neuraminidase , 1999, Journal of Virology.

[38]  Mark von Itzstein,et al.  The war against influenza: discovery and development of sialidase inhibitors. , 2007, Nature reviews. Drug discovery.

[39]  Klaus Schulten,et al.  Molecular Dynamics Simulations Suggest that Electrostatic Funnel Directs Binding of Tamiflu to Influenza N1 Neuraminidases , 2010, PLoS Comput. Biol..

[40]  R. Webster,et al.  Neuraminidase Inhibitor-Resistant Recombinant A/Vietnam/1203/04 (H5N1) Influenza Viruses Retain Their Replication Efficiency and Pathogenicity In Vitro and In Vivo , 2007, Journal of Virology.

[41]  H. Berendsen,et al.  Molecular dynamics with coupling to an external bath , 1984 .

[42]  Jung-Hsin Lin,et al.  Remarkable loop flexibility in avian influenza N1 and its implications for antiviral drug design. , 2007, Journal of the American Chemical Society.

[43]  Supot Hannongbua,et al.  Long time scale GPU dynamics reveal the mechanism of drug resistance of the dual mutant I223R/H275Y neuraminidase from H1N1-2009 influenza virus. , 2012, Biochemistry.

[44]  Thanyada Rungrotmongkol,et al.  Dynamic Behavior of Avian Influenza A Virus Neuraminidase Subtype H5N1 in Complex with Oseltamivir, Zanamivir, Peramivir, and Their Phosphonate Analogues , 2009, J. Chem. Inf. Model..

[45]  Gerhard Klebe,et al.  PDB2PQR: expanding and upgrading automated preparation of biomolecular structures for molecular simulations , 2007, Nucleic Acids Res..

[46]  Huanxiang Liu,et al.  Molecular modeling study on the resistance mechanism of HCV NS3/4A serine protease mutants R155K, A156V and D168A to TMC435. , 2012, Antiviral research.

[47]  William A. McLaughlin,et al.  Evaluating the potency of HIV‐1 protease drugs to combat resistance , 2008, Proteins.

[48]  A. Moscona,et al.  Oseltamivir resistance--disabling our influenza defenses. , 2005, The New England journal of medicine.

[49]  A. Warshel,et al.  Calculations of antibody-antigen interactions: microscopic and semi-microscopic evaluation of the free energies of binding of phosphorylcholine analogs to McPC603. , 1992, Protein engineering.

[50]  Wei Zhang,et al.  Predicting drug resistance of the HIV‐1 protease using molecular interaction energy components , 2009, Proteins.

[51]  Guy Boivin,et al.  Generation and characterization of recombinant pandemic influenza A(H1N1) viruses resistant to neuraminidase inhibitors. , 2011, The Journal of infectious diseases.

[52]  J. Andrew McCammon,et al.  Characterizing Loop Dynamics and Ligand Recognition in Human- and Avian-Type Influenza Neuraminidases via Generalized Born Molecular Dynamics and End-Point Free Energy Calculations , 2009, Journal of the American Chemical Society.

[53]  Xiaojie Xu,et al.  Binding Affinities for a Series of Selective Inhibitors of Gelatinase-A Using Molecular Dynamics with a Linear Interaction Energy Approach , 2001 .

[54]  J. C. Dyason,et al.  A study of the active site of influenza virus sialidase: an approach to the rational design of novel anti-influenza drugs. , 1996, Journal of medicinal chemistry.

[55]  W G Laver,et al.  Structures of aromatic inhibitors of influenza virus neuraminidase. , 1995, Biochemistry.

[56]  Tingjun Hou,et al.  Molecular dynamics and free energy studies on the wild-type and double mutant HIV-1 protease complexed with amprenavir and two amprenavir-related inhibitors: mechanism for binding and drug resistance. , 2007, Journal of medicinal chemistry.

[57]  William L Jorgensen,et al.  Perspective on Free-Energy Perturbation Calculations for Chemical Equilibria. , 2008, Journal of chemical theory and computation.

[58]  J. N. Varghese,et al.  Structure of the catalytic and antigenic sites in influenza virus neuraminidase , 1983, Nature.

[59]  Urban Bren,et al.  Decomposition of the solvation free energies of deoxyribonucleoside triphosphates using the free energy perturbation method. , 2006, The journal of physical chemistry. B.

[60]  J. Åqvist,et al.  Ion-water interaction potentials derived from free energy perturbation simulations , 1990 .

[61]  P. Kollman,et al.  A well-behaved electrostatic potential-based method using charge restraints for deriving atomic char , 1993 .

[62]  Herman Tse,et al.  Two Years after Pandemic Influenza A/2009/H1N1: What Have We Learned? , 2012, Clinical Microbiology Reviews.

[63]  G. Ciccotti,et al.  Numerical Integration of the Cartesian Equations of Motion of a System with Constraints: Molecular Dynamics of n-Alkanes , 1977 .

[64]  D. Case,et al.  Modification of the Generalized Born Model Suitable for Macromolecules , 2000 .

[65]  Jing Zhang,et al.  Detecting and understanding combinatorial mutation patterns responsible for HIV drug resistance , 2010, Proceedings of the National Academy of Sciences.

[66]  David Ozonoff,et al.  Novel Druggable Hot Spots in Avian Influenza Neuraminidase H5N1 Revealed by Computational Solvent Mapping of a Reduced and Representative Receptor Ensemble , 2008, Chemical biology & drug design.

[67]  Mai Suan Li,et al.  Study of Tamiflu Sensitivity to Variants of A/H5N1 Virus Using Different Force Fields , 2011, J. Chem. Inf. Model..

[68]  Martin Almlöf,et al.  Computational study of the binding affinity and selectivity of the bacterial ammonium transporter AmtB. , 2006, Biochemistry.

[69]  Brian J. Smith,et al.  Analysis of inhibitor binding in influenza virus neuraminidase , 2001, Protein science : a publication of the Protein Society.

[70]  D. Case,et al.  Insights into protein-protein binding by binding free energy calculation and free energy decomposition for the Ras-Raf and Ras-RalGDS complexes. , 2003, Journal of molecular biology.

[71]  Guy Boivin,et al.  Susceptibility of recent Canadian influenza A and B virus isolates to different neuraminidase inhibitors. , 2002, Antiviral research.

[72]  T. Šolmajer,et al.  Binding free energy calculations of N-sulphonyl-glutamic acid inhibitors of MurD ligase , 2009, Journal of molecular modeling.

[73]  Jin Woo Park,et al.  Infiltration of Water Molecules into the Oseltamivir-Binding Site of H274Y Neuraminidase Mutant Causes Resistance to Oseltamivir , 2009, J. Chem. Inf. Model..

[74]  Thanyada Rungrotmongkol,et al.  How does each substituent functional group of oseltamivir lose its activity against virulent H5N1 influenza mutants? , 2009, Biophysical chemistry.

[75]  David J. Stevens,et al.  The structure of H5N1 avian influenza neuraminidase suggests new opportunities for drug design , 2006, Nature.

[76]  F. Hayden,et al.  Selection of influenza virus mutants in experimentally infected volunteers treated with oseltamivir. , 2001, The Journal of infectious diseases.

[77]  Thanyada Rungrotmongkol,et al.  Susceptibility of antiviral drugs against 2009 influenza A (H1N1) virus. , 2009, Biochemical and biophysical research communications.

[78]  Xiaojie Xu,et al.  Recent Advances in Free Energy Calculations with a Combination of Molecular Mechanics and Continuum Models , 2006 .

[79]  V. Shanthi,et al.  Exploring the Cause of Oseltamivir Resistance Against Mutant H274Y Neuraminidase by Molecular Simulation Approach , 2012, Applied Biochemistry and Biotechnology.

[80]  Huanxiang Liu,et al.  Exploring the molecular basis of dsRNA recognition by NS1 protein of influenza A virus using molecular dynamics simulation and free energy calculation. , 2011, Antiviral research.

[81]  L. Kit,et al.  A revision of the system of nomenclature for influenza viruses: a WHO memorandum. , 1980, Bulletin of the World Health Organization.

[82]  P. Kollman,et al.  Continuum Solvent Studies of the Stability of DNA, RNA, and Phosphoramidate−DNA Helices , 1998 .

[83]  Tingjun Hou,et al.  Structural Bioinformatics Prediction of Peptides Binding to the Pka Riiα Subunit Using a Hierarchical Strategy , 2022 .

[84]  Huanxiang Liu,et al.  In silico identification of the potential drug resistance sites over 2009 influenza A (H1N1) virus neuraminidase. , 2010, Molecular pharmaceutics.

[85]  M. Zambon,et al.  Emergence of resistance to oseltamivir among influenza A(H1N1) viruses in Europe. , 2008, Euro surveillance : bulletin Europeen sur les maladies transmissibles = European communicable disease bulletin.

[86]  Klaus Schulten,et al.  Molecular Dynamics Simulations Suggest that Electrostatic Funnel Directs Binding of Tamiflu to Influenza N1 Neuraminidases , 2010, PLoS Comput. Biol..

[87]  Wei Zhang,et al.  A point‐charge force field for molecular mechanics simulations of proteins based on condensed‐phase quantum mechanical calculations , 2003, J. Comput. Chem..

[88]  J. McKimm-Breschkin,et al.  Resistance of influenza viruses to neuraminidase inhibitors--a review. , 2000, Antiviral research.

[89]  Alba Chimirri,et al.  Molecular dynamics studies of the wild-type and double mutant HIV-1 integrase complexed with the 5CITEP inhibitor: mechanism for inhibition and drug resistance. , 2003, Biophysical journal.

[90]  Tingjun Hou,et al.  Assessing the performance of the molecular mechanics/Poisson Boltzmann surface area and molecular mechanics/generalized Born surface area methods. II. The accuracy of ranking poses generated from docking , 2011, J. Comput. Chem..

[91]  Thanyada Rungrotmongkol,et al.  Understanding of known drug‐target interactions in the catalytic pocket of neuraminidase subtype N1 , 2008, Proteins.

[92]  H. Gohlke,et al.  Free Energy Calculations by the Molecular Mechanics Poisson−Boltzmann Surface Area Method , 2012, Molecular informatics.

[93]  Robert V. Swift,et al.  Mechanism of 150-cavity formation in influenza neuraminidase , 2011, Nature communications.

[94]  K. Goa,et al.  Zanamivir , 1999, Drugs.

[95]  Huanxiang Liu,et al.  Understanding the effect of drug-resistant mutations of HIV-1 intasome on raltegravir action through molecular modeling study. , 2012, Molecular bioSystems.

[96]  Alan J. Hay,et al.  Crystal structures of oseltamivir-resistant influenza virus neuraminidase mutants , 2008, Nature.

[97]  Lieve Naesens,et al.  Design and synthesis of bioactive adamantane spiro heterocycles. , 2007, Bioorganic & medicinal chemistry letters.

[98]  Ray Luo,et al.  Accelerated Poisson–Boltzmann calculations for static and dynamic systems , 2002, J. Comput. Chem..

[99]  Tatiana Tatusova,et al.  NCBI Reference Sequence (RefSeq): a curated non-redundant sequence database of genomes, transcripts and proteins , 2004, Nucleic Acids Res..

[100]  Holger Gohlke,et al.  The Amber biomolecular simulation programs , 2005, J. Comput. Chem..

[101]  Thanyada Rungrotmongkol,et al.  Source of oseltamivir resistance in avian influenza H5N1 virus with the H274Y mutation , 2009, Amino Acids.

[102]  Urban Bren,et al.  DNA duplex stability: the role of preorganized electrostatics. , 2010, The journal of physical chemistry. B.

[103]  A. Mankertz,et al.  Measles outbreak in Switzerland--an update relevant for the European football championship (EURO 2008). , 2008, Euro surveillance : bulletin Europeen sur les maladies transmissibles = European communicable disease bulletin.

[104]  T. B. Usman,et al.  Reduced Sensitivity of Influenza A (H5N1) to Oseltamivir , 2007, Emerging Infectious Diseases.

[105]  Wei Zhang,et al.  Characterization of Domain-Peptide Interaction Interface , 2009, Molecular & Cellular Proteomics.

[106]  Zhihua Lin,et al.  Molecular dynamics simulation of oseltamivir resistance in neuraminidase of avian influenza H5N1 virus , 2011, Journal of molecular modeling.