Abstract:The binding of a series of hydroxamate inhibitors with gelatinase-A is examined to evaluate the viability of calculating free energies of binding, DeltaG(b), utilizing molecular dynamics (MD) simulations with a linear interaction energy approach. In our simulations, a bonded model was used to represent the potentials of the catalytic zinc center. The electrostatic distribution of this model was derived using a two-stage electrostatic potential fitting calculations. The resulting bonded model was then used to generate the MD trajectories. Coulombic, van der Waals, and coordinate bond energy components determined from MD simulations of the bound and unbound inhibitors solvated in water were correlated with the free energies of binding for the 15 hydroxamate inhibitors. In the correlation process, several linear models consisted of different energy components were tested. We found that besides the usually used Coulombic and van der Waals energy terms, the introduction of a constant term could significantly improve the correlation. The best model yields an average error of 0.6 kcal/mol for the 15 binding affinities, which cover an observed range of 7.2 kcal/mol. The predictive ability of the best model was revealed by the high value of q(2) (0.854) from the leave-one-out cross-validation. To this series of inhibitors, the constant term can be treated as effective adjustment to the entropy contribution in the binding free energies. The MD simulations predicted the binding mode of the gelatinase-A with the studied inhibitors, and also provided insights into the interactions occurring in the active site and the origins of variations in DeltaG(b) The P1 ' groups of inhibitors make extensive van der Waals and hydrophobic contacts with the nonpolar side chains of four residues in the S1 ' subsite, including Leu 197, Val 198, Leu 218, and Tyr 223, which directly influence the ligand binding. Hydrogen bonds between hydroxamates and gelatinase-A are very important to stabilize the inhibitors in the active site. The hydrogen bonds between the P-3' group and gelatinase-A can produce more favorable electrostatic interactions.

The Prediction of Binding Free Energies for a Group of Trypsin Inhibitors

Prediction of binding free energies between benzamides and trypsin

Predictions Of Binding Of A Diverse Set Of Ligands To Gelatinase-A By A Combination Of Molecular Dynamics And Continuum Solvent Models

Binding Affinities for A Series of Selective Inhibitors of Gelatinase-A Using Molecular Dynamics with A Linear Interaction Energy Approach

Predicted Par1 Inhibitors from Multiple Computational Methods

Insights into the binding mechanisms of inhibitors of MDM2 based on molecular dynamics simulations and binding free energy calculations

Binding Modes and Conformational Changes of FK506-binding Protein 51 Induced by Inhibitor Bindings: Insight into Molecular Mechanisms Based on Multiple Simulation Technologies

Accurate Prediction for Protein-Peptide Binding Based on High-Temperature Molecular Dynamics Simulations

Binding Mechanism of Inhibitors to Heat Shock Protein 90 Investigated by Multiple Independent Molecular Dynamics Simulations and Prediction of Binding Free Energy

Free Energy Profile and Kinetics of Coupled Folding and Binding of the Intrinsically Disordered Protein p53 with MDM2

Binding Free Energy Calculations For Mmp2-Hydroxamate Complexes

Rapid and Reliable Binding Affinity Prediction of Bromodomain Inhibitors: A Computational Study

Fast and accurate predictions of binding free energies using MM-PBSA and MM-GBSA

A Computational Analysis of the Binding Affinities of FKBP12 Inhibitors Using the MM‐PB/SA Method

Non-bonded force field parameters from MBIS partitioning of the molecular electron density improve binding affinity predictions of the T4-lysozyme double mutant

Revealing The Binding Modes And The Unbinding Of 14-3-3 Sigma Proteins And Inhibitors By Computational Methods

Binding free energy of protein/ligand complexes calculated using dissociation Parallel Cascade Selection Molecular Dynamics and Markov state model

Investigate the Binding of Catechins to Trypsin Using Docking and Molecular Dynamics Simulation

Predicting the Protein-Ligand Affinity from Molecular Dynamics Trajectories

Nonbonded Force Field Parameters from MBIS Partitioning of the Molecular Electron Density Improve Binding Affinity Predictions of the T4-Lysozyme Double Mutant

Molecular Dynamic Simulations Give Insight into the Mechanism of Binding Between 2-Aminothiazole Inhibitors and CDK5