Insights into the structure-affinity relationships and solvation effects between OfHex1 and inhibitors using molecular dynamics simulations.

OfHex1 is a potential target for the rational design of pesticides. TMG-chitotriomycin is one of the most highly specific known inhibitors of chitinolytic β-GlcNAcases from bacteria, fungi and insects. TMG-chitotriomycin and its analogues show different activities to OfHex1, dependent on the number of GlcNAc units. Subsequently, it is essential to explore how these GlcNAc unit number changes cause alterations in activity. In this study, we examined the free energy patterns and per residue decomposition of binding within the complexes of OfHex1 and a series of inhibitors, utilizing restricted molecular dynamics (MD) and water-mediated MM/GBSA calculations. The results indicated Glu328 could form a stronger polar interaction with OfHex1 inhibitors, while Trp448 and Trp490 had important non-polar contributions. Interestingly, the conformation of Trp448 was different in the open or closed state, when OfHex1 bound different inhibitors. Moreover, the water molecule that mediates the GlcNAc Ⅱ and Trp490 may be critical to stabilizing the hydrophobic interaction. Further study showed that isomerization of TMG-chitotriomycin analogs did not decrease binding affinity, however, there was a highly positive correlation between the calculated binding affinities and the experimental activity data (r2 = 0.92) when water molecules were explicitly taken into account. Moreover, the water molecules that mediated GlcNAc II and Trp490 might be critical to the stabilization of the hydrophobic interaction and cause the activity difference between TMG-(GlcNAc)2 and TMG-(GlcNAc).

[1]  T. Hatanaka,et al.  TMG-chitotriomycin as a probe for the prediction of substrate specificity of β-N-acetylhexosaminidases. , 2013, Carbohydrate research.

[2]  Guang-Fu Yang,et al.  Computational Design and Discovery of Conformationally Flexible Inhibitors of Acetohydroxyacid Synthase to Overcome Drug Resistance Associated with the W586L Mutation , 2008, ChemMedChem.

[3]  S. Withers,et al.  NAG-thiazoline, An N-Acetyl-β-hexosaminidase Inhibitor That Implicates Acetamido Participation , 1996 .

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

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

[6]  X. Qian,et al.  A crystal structure-guided rational design switching non-carbohydrate inhibitors' specificity between two β-GlcNAcase homologs , 2014, Scientific Reports.

[7]  Q. Yang,et al.  Exploring NAG‐thiazoline and its derivatives as inhibitors of chitinolytic β‐acetylglucosaminidases , 2015, FEBS letters.

[8]  D. V. van Aalten,et al.  GlcNAcstatins are nanomolar inhibitors of human O-GlcNAcase inducing cellular hyper-O-GlcNAcylation , 2009, The Biochemical journal.

[9]  G. Pavesi,et al.  Phylogenetic analyses suggest multiple changes of substrate specificity within the Glycosyl hydrolase 20 family , 2008, BMC Evolutionary Biology.

[10]  Christopher I. Bayly,et al.  Fast, efficient generation of high‐quality atomic charges. AM1‐BCC model: II. Parameterization and validation , 2002, J. Comput. Chem..

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

[12]  T. Aoyagi,et al.  The structure of nagstatin, a new inhibitor of N-acetyl-beta-D-glucosaminidase. , 1992, The Journal of antibiotics.

[13]  Qing Yang,et al.  Structural Determinants of an Insect β-N-Acetyl-d-hexosaminidase Specialized as a Chitinolytic Enzyme* , 2010, The Journal of Biological Chemistry.

[14]  A. Vasella,et al.  Synthesis of N‐Acetylglucosamine‐Derived Nagstatin Analogues and Their Evaluation as Glycosidase Inhibitors , 2005 .

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

[16]  Q. Ma,et al.  Structural insights into chitinolytic enzymes and inhibition mechanisms of selective inhibitors. , 2014, Current pharmaceutical design.

[17]  M. Horsch,et al.  N‐Acetylgucosaminono‐1,5‐lactone oxime and the corresponding (phenylcarbamoyl)oxime , 1991 .

[18]  Ge-Fei Hao,et al.  Understanding the mechanism of drug resistance due to a codon deletion in protoporphyrinogen oxidase through computational modeling. , 2009, The journal of physical chemistry. B.

[19]  J. Andrew McCammon,et al.  MM-PBSA Captures Key Role of Intercalating Water Molecules at a Protein−Protein Interface , 2009, Journal of chemical theory and computation.

[20]  Yong-Liang Zhu,et al.  Including Explicit Water Molecules as Part of the Protein Structure in MM/PBSA Calculations , 2014, J. Chem. Inf. Model..

[21]  D. Case,et al.  Exploring protein native states and large‐scale conformational changes with a modified generalized born model , 2004, Proteins.

[22]  Yonghua Wang,et al.  Identification of Novel Potential β-N-Acetyl-D-Hexosaminidase Inhibitors by Virtual Screening, Molecular Dynamics Simulation and MM-PBSA Calculations , 2012, International journal of molecular sciences.

[23]  T. Iwashita,et al.  TMG-chitotriomycin, an enzyme inhibitor specific for insect and fungal beta-N-acetylglucosaminidases, produced by actinomycete Streptomyces anulatus NBRC 13369. , 2008, Journal of the American Chemical Society.

[24]  Alessandro Contini,et al.  Explicit Ligand Hydration Shells Improve the Correlation between MM-PB/GBSA Binding Energies and Experimental Activities. , 2013, Journal of chemical theory and computation.

[25]  Chang-Guo Zhan,et al.  Modeling the catalysis of anti-cocaine catalytic antibody: competing reaction pathways and free energy barriers. , 2008, Journal of the American Chemical Society.

[26]  Tian Liu,et al.  Synthesis, Evaluation, and Mechanism of N,N,N‐Trimethyl‐D‐glucosamine‐(1→4)‐chitooligosaccharides as Selective Inhibitors of Glycosyl Hydrolase Family 20 β‐N‐Acetyl‐D‐hexosaminidases , 2011, Chembiochem : a European journal of chemical biology.

[27]  Chi‐Huey Wong,et al.  Novel Five‐Membered Iminocyclitol Derivatives as Selective and Potent Glycosidase Inhibitors: New Structures for Antivirals and Osteoarthritis , 2006, Chembiochem : a European journal of chemical biology.

[28]  T. Butters,et al.  Scalable syntheses of both enantiomers of DNJNAc and DGJNAc from glucuronolactone: the effect of N-alkylation on hexosaminidase inhibition. , 2012, Chemistry.

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

[30]  Dan Li,et al.  Importance of protein flexibility in ranking inhibitor affinities: modeling the binding mechanisms of piperidine carboxamides as Type I1/2 ALK inhibitors. , 2015, Physical chemistry chemical physics : PCCP.

[31]  R. Dror,et al.  Improved side-chain torsion potentials for the Amber ff99SB protein force field , 2010, Proteins.

[32]  Holger Gohlke,et al.  MMPBSA.py: An Efficient Program for End-State Free Energy Calculations. , 2012, Journal of chemical theory and computation.