The MM2QM tool for combining docking, molecular dynamics, molecular mechanics, and quantum mechanics

The use of the MM2QM tool in a combined docking + molecular dynamics (MD) + molecular mechanics (MM) + quantum mechanical (QM) binding affinity prediction study is presented, and the tool itself is discussed. The system of interest is Mycobacterium tuberculosis (MTB) pantothenate synthetase in complexes with three highly similar sulfonamide inhibitors, for which crystal structures are available. Starting from the structure of MTB pantothenate synthetase in the "open" conformation and following the combined docking + MD + MM + QM procedure, we were able to capture the closing of the enzyme binding pocket and to reproduce the position of the ligands with an average root mean square deviation of 1.6 Å. Protein-ligand interaction energies were reproduced with an average error lower than 10%. The discussion on the MD part and a protein flexibility importance is carried out. The presented approach may be useful especially for finding analog inhibitors or improving drug candidates.