Molecular dynamics (MD) simulations and binding free energy calculation studies between inhibitors and type II dehydroquinase (DHQ2)

Type II dehydroquinase (DHQ2) is the third enzyme of the shikimic acid pathway, and it has been the effective target for tuberculosis (TB). So far, developing multiple potent inhibitors of the DHQ2 of Mycobacterium tuberculosis (DHQ2-Mt) has been considered to be the new therapy to TB. Molecular dynamics simulations followed by molecular mechanics-generalised Born surface area were carried out to calculate the free binding energy and to determine the affinity ability of the four chosen inhibitor molecules, L1, L2, L3 and L4. Energy decomposition analyses show the electrostatic interaction and van der Waals interaction of the ligands to every residue of the DHQ2-Mt. The results suggest that some important residues have different interactions with the four ligands, such as Arg19 and Tyr24. These interactions may have an effect on the ligand binding affinity. The binding affinity of monosubstituted inhibitor is higher than that of disubstituted inhibitor, due to some important interactions with the DHQ2-Mt residues. These computational works will be significant to the theoretical research in the future.