Many experimental studies have found that flavonoids can inhibit the activities of matrix metalloproteinases (MMPs), but the relevant mechanisms are still unclear. In this paper, the interaction mechanisms of MMP-9 with its five flavonoid inhibitors are investigated using a combination of molecular docking, hybrid quantum mechanical and molecular mechanical (QM/MM) calculations, and molecular dynamics simulations. The molecular dynamics simulation results show a good linear correlation between the calculated binding free energies of QM/MM−Poisson–Boltzmann surface area (PBSA) and the experimental −log(EC50) regarding the studied five flavonoids on MMP-9 inhibition in explicit solvent. It is found that compared with the MM−PBSA method, the QM/MM−PBSA method can obviously improve the accuracy for the calculated binding free energies. The predicted binding modes of the five flavonoid−MMP-9 complexes reveal that the different hydrogen bond networks can form besides producing the Zn−O coordination bonds, which can reasonably explain previous experimental results. The agreement between our calculated results and the previous experimental facts indicates that the force field parameters used here are effective and reliable for investigating the systems of flavonoid−MMP-9 interactions, and thus, these simulations and analyses could be reproduced for the other related systems involving protein−ligand interactions. This paper may be helpful for designing the new MMP-9 inhibitors having higher biological activities by carrying out the structural modifications of flavonoid molecules.
Protein-ligand free energies of binding from full-protein DFT calculations: convergence and choice of exchange-correlation functional