Free-energy analysis of lysozyme–triNAG binding modes with all-atom molecular dynamics simulation combined with the solution theory in the energy representation

Abstract We propose a method for calculating the binding free energy of protein-ligand complexes using all-atom molecular dynamics simulation combined with the solution theory in the energy representation. Four distinct modes for the binding of tri- N -acetyl- d -glucosamine (triNAG) to hen egg-white lysozyme were investigated, one from the crystal structure and three generated by docking predictions. The proposed method was demonstrated to be used to distinguish the most plausible binding mode (crystal model) as the lowest binding energy mode.

[1]  N. Matubayasi,et al.  Hydration property of globular proteins: An analysis of solvation free energy by energy representation method , 2010 .

[2]  David L Mobley,et al.  Alchemical free energy methods for drug discovery: progress and challenges. , 2011, Current opinion in structural biology.

[3]  J F Brandts,et al.  Rapid measurement of binding constants and heats of binding using a new titration calorimeter. , 1989, Analytical biochemistry.

[4]  David L. Mobley,et al.  Chapter 4 Alchemical Free Energy Calculations: Ready for Prime Time? , 2007 .

[5]  Nobuyuki Matubayasi,et al.  Theory of solutions in the energy representation. II. Functional for the chemical potential , 2002 .

[6]  Alexander D. MacKerell,et al.  CHARMM additive all-atom force field for carbohydrate derivatives and its utility in polysaccharide and carbohydrate-protein modeling. , 2011, Journal of chemical theory and computation.

[7]  Nobuyuki Matubayasi,et al.  Theory of solutions in the energy representation. III. Treatment of the molecular flexibility , 2003 .

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

[9]  J. Aqvist,et al.  A new method for predicting binding affinity in computer-aided drug design. , 1994, Protein engineering.

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

[11]  Laxmikant V. Kalé,et al.  Scalable molecular dynamics with NAMD , 2005, J. Comput. Chem..

[12]  N. Matubayasi,et al.  Free-energy analysis of the molecular binding into lipid membrane with the method of energy representation. , 2008, The Journal of chemical physics.

[13]  Nobuyuki Matubayasi,et al.  Evaluation of protein-protein docking model structures using all-atom molecular dynamics simulations combined with the solution theory in the energy representation. , 2012, The Journal of chemical physics.

[14]  Irwin Oppenheim,et al.  Statistical Mechanical Theory of Transport Processes. VII. The Coefficient of Thermal Conductivity of Monatomic Liquids , 1954 .

[15]  David S. Goodsell,et al.  Automated docking using a Lamarckian genetic algorithm and an empirical binding free energy function , 1998 .

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

[17]  René Thomsen,et al.  MolDock: a new technique for high-accuracy molecular docking. , 2006, Journal of medicinal chemistry.

[18]  P M Cullis,et al.  Affinities of amino acid side chains for solvent water. , 1981, Biochemistry.

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

[20]  Nobuyuki Matubayasi,et al.  End-point calculation of solvation free energy of amino-acid analogs by molecular theories of solution , 2010 .

[21]  Peter A. Kollman,et al.  FREE ENERGY CALCULATIONS : APPLICATIONS TO CHEMICAL AND BIOCHEMICAL PHENOMENA , 1993 .

[22]  Fumio Hirata,et al.  The effects of solvent on the conformation and the collective motions of protein: normal mode analysis and molecular dynamics simulations of melittin in water and in vacuum , 1991 .

[23]  Ioan Andricioaei,et al.  On the calculation of entropy from covariance matrices of the atomic fluctuations , 2001 .

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

[25]  N. Matubayasi,et al.  Communication: Free-energy analysis of hydration effect on protein with explicit solvent: equilibrium fluctuation of cytochrome c. , 2011, The Journal of chemical physics.

[26]  P Willett,et al.  Development and validation of a genetic algorithm for flexible docking. , 1997, Journal of molecular biology.

[27]  M. Vaney,et al.  High-resolution structure (1.33 A) of a HEW lysozyme tetragonal crystal grown in the APCF apparatus. Data and structural comparison with a crystal grown under microgravity from SpaceHab-01 mission. , 1996, Acta crystallographica. Section D, Biological crystallography.

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

[29]  Nobuyuki Matubayasi,et al.  Theory of solutions in the energetic representation. I. Formulation , 2000 .

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

[31]  Alexander D. MacKerell,et al.  Extending the treatment of backbone energetics in protein force fields: Limitations of gas‐phase quantum mechanics in reproducing protein conformational distributions in molecular dynamics simulations , 2004, J. Comput. Chem..

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