Computational Study of Drug Binding Affinity to Influenza A Neuraminidase Using Smooth Reaction Path Generation (SRPG) Method

Assessment of accurate drug binding affinity to a protein remains a challenge for in silico drug development. In this research, we used the smooth reaction path generation (SRPG) method to calculate binding free energies and determine potential of mean forces (PMFs) along the smoothed dissociation paths of influenza A neuraminidase and its variants with oseltamivir (Tamiflu) and zanamivir (Relenza) inhibitors. With the gained results, we found that the binding free energies of neuraminidase A/H5N1 in WT and two mutants (including H274Y and N294S) with oseltamivir and zanamivir show good agreement with experimental results. Additionally, the thermodynamic origin of the drug resistance of the mutants was also discussed from the PMF profiles.

[1]  David S. Goodsell,et al.  AutoDock4 and AutoDockTools4: Automated docking with selective receptor flexibility , 2009, J. Comput. Chem..

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

[3]  A. Laio,et al.  Flexible docking in solution using metadynamics. , 2005, Journal of the American Chemical Society.

[4]  J M Thornton,et al.  LIGPLOT: a program to generate schematic diagrams of protein-ligand interactions. , 1995, Protein engineering.

[5]  Ly Le,et al.  Steered molecular dynamics approach for promising drugs for influenza A virus targeting M2 channel proteins , 2015, European Biophysics Journal.

[6]  Ly Le,et al.  Discover binding pathways using the sliding binding-box docking approach: application to binding pathways of oseltamivir to avian influenza H5N1 neuraminidase , 2013, Journal of Computer-Aided Molecular Design.

[7]  Lin Li,et al.  Insights into susceptibility of antiviral drugs against the E119G mutant of 2009 influenza A (H1N1) neuraminidase by molecular dynamics simulations and free energy calculations. , 2013, Antiviral research.

[8]  Haruki Nakamura,et al.  Protein-Ligand Binding Free Energy Calculation by the Smooth Reaction Path Generation (SRPG) Method , 2009, J. Chem. Inf. Model..

[9]  Tingjun Hou,et al.  Development and Evaluation of an Integrated Virtual Screening Strategy by Combining Molecular Docking and Pharmacophore Searching Based on Multiple Protein Structures , 2013, J. Chem. Inf. Model..

[10]  H. Fujitani,et al.  Massively parallel computation of absolute binding free energy with well-equilibrated states. , 2009, Physical review. E, Statistical, nonlinear, and soft matter physics.

[11]  Y. Fukunishi,et al.  The Filling Potential Method: A Method for Estimating the Free Energy Surface for Protein−Ligand Docking , 2003 .

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

[13]  Francesco Luigi Gervasio,et al.  From A to B in free energy space. , 2007, The Journal of chemical physics.

[14]  Serdar Kuyucak,et al.  Potential of mean force calculations of ligand binding to ion channels from Jarzynski's equality and umbrella sampling. , 2008, The Journal of chemical physics.

[15]  J. Tam,et al.  Influenza A (H5N1) in Hong Kong: an overview. , 2002, Vaccine.

[16]  Tingjun Hou,et al.  Insight into Crizotinib Resistance Mechanisms Caused by Three Mutations in ALK Tyrosine Kinase using Free Energy Calculation Approaches , 2013, J. Chem. Inf. Model..

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

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

[19]  Mai Mai,et al.  Steered Molecular Dynamics-A Promising Tool for Drug Design , 2012 .

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

[21]  Roman A. Laskowski,et al.  LigPlot+: Multiple Ligand-Protein Interaction Diagrams for Drug Discovery , 2011, J. Chem. Inf. Model..

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

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

[24]  D. Baker,et al.  A simple physical model for binding energy hot spots in protein–protein complexes , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[25]  K Schulten,et al.  VMD: visual molecular dynamics. , 1996, Journal of molecular graphics.

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

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

[28]  Leslie Greengard,et al.  A fast algorithm for particle simulations , 1987 .

[29]  Tingjun Hou,et al.  P-loop Conformation Governed Crizotinib Resistance in G2032R-Mutated ROS1 Tyrosine Kinase: Clues from Free Energy Landscape , 2014, PLoS Comput. Biol..

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

[31]  Tingjun Hou,et al.  Assessing an Ensemble Docking-Based Virtual Screening Strategy for Kinase Targets by Considering Protein Flexibility , 2014, J. Chem. Inf. Model..

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

[33]  R. Webster,et al.  Pandemic influenza as a current threat. , 2009, Current topics in microbiology and immunology.

[34]  Lin Li,et al.  Theoretical Studies on the Susceptibility of Oseltamivir against Variants of 2009 A/H1N1 Influenza Neuraminidase , 2012, J. Chem. Inf. Model..

[35]  Tingjun Hou,et al.  Assessing the performance of MM/PBSA and MM/GBSA methods. 5. Improved docking performance using high solute dielectric constant MM/GBSA and MM/PBSA rescoring. , 2014, Physical chemistry chemical physics : PCCP.

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

[37]  Jaques Reifman,et al.  Quantitative Predictions of Binding Free Energy Changes in Drug-Resistant Influenza Neuraminidase , 2012, PLoS Comput. Biol..