Binding structures of tri‐N‐acetyl‐β‐glucosamine in hen egg white lysozyme using molecular dynamics with a polarizable force field

Lysozyme is a well‐studied enzyme that hydrolyzes the β‐(1,4)‐glycosidic linkage of N‐acetyl‐β‐glucosamine (NAG)n oligomers. The active site of hen egg‐white lysozyme (HEWL) is believed to consist of six subsites, A‐F that can accommodate six sugar residues. We present studies exploring the use of polarizable force fields in conjunction with all‐atom molecular dynamics (MD) simulations to analyze binding structures of complexes of lysozyme and NAG trisaccharide, (NAG)3. MD trajectories are applied to analyze structures and conformation of the complex as well as protein–ligand interactions, including the hydrogen‐bonding network in the binding pocket. Two binding modes (ABC and BCD) of (NAG)3 are investigated independently based on a fixed‐charge model and a polarizable model. We also apply molecular mechanics with generalized born and surface area (MM‐GBSA) methods based on MD using both nonpolarizable and polarizable force fields to compute binding free energies. We also study the correlation between root‐mean‐squared deviation and binding free energies of the wildtype and W62Y mutant; we find that for this prototypical system, approaches using the MD trajectories coupled with implicit solvent models are equivalent for polarizable and fixed‐charge models. © 2012 Wiley Periodicals, Inc.

[1]  Yang Zhong,et al.  Solvation properties of N‐acetyl‐β‐glucosamine: Molecular dynamics study incorporating electrostatic polarization , 2011, J. Comput. Chem..

[2]  Ray Luo,et al.  Virtual screening using molecular simulations , 2011, Proteins.

[3]  Georgios Archontis,et al.  Kinetics, in silico docking, molecular dynamics, and MM‐GBSA binding studies on prototype indirubins, KT5720, and staurosporine as phosphorylase kinase ATP‐binding site inhibitors: The role of water molecules examined , 2011, Proteins.

[4]  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..

[5]  Nohad Gresh,et al.  Polarizable water molecules in ligand-macromolecule recognition. Impact on the relative affinities of competing pyrrolopyrimidine inhibitors for FAK kinase. , 2010, Journal of the American Chemical Society.

[6]  Osamu Miyashita,et al.  Solution and crystal molecular dynamics simulation study of m4-cyanovirin-N mutants complexed with di-mannose. , 2009, Biophysical journal.

[7]  M. Matsushima,et al.  Importance of the hydrogen bonding network including Asp52 for catalysis, as revealed by Asn59 mutant hen egg-white lysozymes. , 2009, Journal of biochemistry.

[8]  Guohui Li,et al.  Trypsin‐ligand binding free energies from explicit and implicit solvent simulations with polarizable potential , 2009, J. Comput. Chem..

[9]  Jianpeng Ma,et al.  CHARMM: The biomolecular simulation program , 2009, J. Comput. Chem..

[10]  Alexander D. MacKerell,et al.  Additive empirical force field for hexopyranose monosaccharides , 2008, J. Comput. Chem..

[11]  Sandeep Patel,et al.  Revised charge equilibration potential for liquid alkanes. , 2008, The journal of physical chemistry. B.

[12]  Pengyu Y. Ren,et al.  Calculation of protein–ligand binding free energy by using a polarizable potential , 2008, Proceedings of the National Academy of Sciences.

[13]  D. F. Green.  Optimized parameters for continuum solvation calculations with carbohydrates. , 2008, Journal of Physical Chemistry B.

[14]  B. Berne,et al.  Role of the active-site solvent in the thermodynamics of factor Xa ligand binding. , 2008, Journal of the American Chemical Society.

[15]  Haruki Nakamura,et al.  Protein‐inhibitor flexible docking by a multicanonical sampling: Native complex structure with the lowest free energy and a free‐energy barrier distinguishing the native complex from the others , 2007, Proteins.

[16]  P. Konidala,et al.  Molecular dynamics simulations of pea (Pisum sativum) lectin structure with octyl glucoside detergents: the ligand interactions and dynamics. , 2007, Biophysical chemistry.

[17]  M. Gilson,et al.  Calculation of protein-ligand binding affinities. , 2007, Annual review of biophysics and biomolecular structure.

[18]  C. Laughton,et al.  A simple physical description of DNA dynamics: quasi-harmonic analysis as a route to the configurational entropy , 2007, Journal of physics. Condensed matter : an Institute of Physics journal.

[19]  Robert Abel,et al.  Motifs for molecular recognition exploiting hydrophobic enclosure in protein–ligand binding , 2007, Proceedings of the National Academy of Sciences.

[20]  N. Foloppe,et al.  Towards predictive ligand design with free-energy based computational methods? , 2006, Current medicinal chemistry.

[21]  Charles L. Brooks,et al.  Fluctuating charge force fields: recent developments and applications from small molecules to macromolecular biological systems , 2006 .

[22]  Mark A Olson,et al.  Calculation of absolute protein-ligand binding affinity using path and endpoint approaches. , 2006, Biophysical journal.

[23]  B. Kuhn,et al.  Validation and use of the MM-PBSA approach for drug discovery. , 2005, Journal of medicinal chemistry.

[24]  Alexander D. MacKerell,et al.  CHARMM fluctuating charge force field for proteins: II Protein/solvent properties from molecular dynamics simulations using a nonadditive electrostatic model , 2004, J. Comput. Chem..

[25]  Holger Gohlke,et al.  Converging free energy estimates: MM‐PB(GB)SA studies on the protein–protein complex Ras–Raf , 2004, J. Comput. Chem..

[26]  Charles L. Brooks,et al.  CHARMM fluctuating charge force field for proteins: I parameterization and application to bulk organic liquid simulations , 2004, J. Comput. Chem..

[27]  Charles L. Brooks,et al.  Generalized born model with a simple smoothing function , 2003, J. Comput. Chem..

[28]  Alexander D. MacKerell,et al.  A simple polarizable model of water based on classical Drude oscillators , 2003 .

[29]  Benoît Roux,et al.  Modeling induced polarization with classical Drude oscillators: Theory and molecular dynamics simulation algorithm , 2003 .

[30]  Pengyu Y. Ren,et al.  Polarizable Atomic Multipole Water Model for Molecular Mechanics Simulation , 2003 .

[31]  Pengyu Y. Ren,et al.  Consistent treatment of inter‐ and intramolecular polarization in molecular mechanics calculations , 2002, J. Comput. Chem..

[32]  P. Procacci,et al.  A transferable polarizable electrostatic force field for molecular mechanics based on the chemical potential equalization principle , 2002 .

[33]  Kevin J. Naidoo,et al.  Carbohydrate solution simulations: Producing a force field with experimentally consistent primary alcohol rotational frequencies and populations , 2002, J. Comput. Chem..

[34]  S. Withers,et al.  Catalysis by hen egg-white lysozyme proceeds via a covalent intermediate , 2001, Nature.

[35]  S. Rick.  Simulations of ice and liquid water over a range of temperatures using the fluctuating charge model , 2001 .

[36]  P. Kollman,et al.  Combined molecular mechanical and continuum solvent approach (MM-PBSA/GBSA) to predict ligand binding , 2000 .

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

[38]  Alexander D. MacKerell,et al.  All-atom empirical potential for molecular modeling and dynamics studies of proteins. , 1998, The journal of physical chemistry. B.

[39]  Peter Itskowitz,et al.  Chemical Potential Equalization Principle: Direct Approach from Density Functional Theory , 1997 .

[40]  D. Beglov,et al.  Atomic Radii for Continuum Electrostatics Calculations Based on Molecular Dynamics Free Energy Simulations , 1997 .

[41]  Bruce J. Berne,et al.  Dynamical Fluctuating Charge Force Fields: The Aqueous Solvation of Amides , 1996 .

[42]  K. Watanabe,et al.  Dissection of protein-carbohydrate interactions in mutant hen egg-white lysozyme complexes and their hydrolytic activity. , 1995, Journal of molecular biology.

[43]  Robert J. Woods,et al.  Molecular Mechanical and Molecular Dynamic Simulations of Glycoproteins and Oligosaccharides. 1. GLYCAM_93 Parameter Development , 1995 .

[44]  B. Berne,et al.  Dynamical fluctuating charge force fields: Application to liquid water , 1994, chem-ph/9406002.

[45]  M Karplus,et al.  The contribution of vibrational entropy to molecular association. The dimerization of insulin. , 1994, Journal of molecular biology.

[46]  J. Schlitter.  Estimation of absolute and relative entropies of macromolecules using the covariance matrix , 1993 .

[47]  M. Klein,et al.  Nosé-Hoover chains : the canonical ensemble via continuous dynamics , 1992 .

[48]  D. Phillips,et al.  Refinement of an enzyme complex with inhibitor bound at partial occupancy. Hen egg-white lysozyme and tri-N-acetylchitotriose at 1.75 A resolution. , 1992, Journal of molecular biology.

[49]  I. Kumagai,et al.  Redesign of the substrate-binding site of hen egg white lysozyme based on the molecular evolution of C-type lysozymes. , 1992, The Journal of biological chemistry.

[50]  M. James,et al.  Lysozyme revisited: crystallographic evidence for distortion of an N-acetylmuramic acid residue bound in site D. , 1991, Journal of molecular biology.

[51]  W. Goddard,et al.  Charge equilibration for molecular dynamics simulations , 1991 .

[52]  Car,et al.  Unified approach for molecular dynamics and density-functional theory. , 1985, Physical review letters.

[53]  M. Karplus,et al.  Method for estimating the configurational entropy of macromolecules , 1981 .

[54]  J. A. Rupley,et al.  STUDIES ON THE ENZYMIC ACTIVITY OF LYSOZYME, II. THE HYDROLYSIS AND TRANSFER REACTIONS OF N-ACETYLGLUCOSAMINE OLIGOSACCHARIDES*. , 1967 .

[55]  D. Phillips,et al.  THE HEN EGG-WHITE LYSOZYME MOLECULE , 1967 .

[56]  D. Phillips,et al.  The three-dimensional structure of an enzyme molecule. , 1966, Scientific American.

[57]  R. T. Sanderson,et al.  An Interpretation of Bond Lengths and a Classification of Bonds. , 1951, Science.

[58]  Hidalgo Raul Cruz,et al.  ロッド‐コイル ジブロック共重合体におけるスメクチック相 , 2007 .

[59]  R. Dreele.  Binding of N-acetylglucosamine oligosaccharides to hen egg-white lysozyme: a powder diffraction study. , 2005 .

[60]  C. Brooks.  Computer simulation of liquids , 1989 .

[61]  W. C. Still,et al.  A rapid approximation to the solvent accessible surface areas of atoms , 1988 .

[62]  L. Johnson,et al.  Protein-oligosaccharide interactions: lysozyme, phosphorylase, amylases. , 1988, Current topics in microbiology and immunology.

[63]  Wilfried J. Mortier,et al.  Electronegativity-equalization method for the calculation of atomic charges in molecules , 1986 .

[64]  Taiji Imoto,et al.  21 Vertebrate Lysozymes , 1972 .

[65]  M. Born.  Statistical Thermodynamics , 1944, Nature.