Matrix metalloproteinases are zinc-containing enzymes capable of degrading all components of the extracellular matrix. Owing to their role in human disease, matrix metalloproteinase have been the subject of extensive study. A bioinorganic approach was recently used to identify novel inhibitors based on a maltol zinc-binding group, but accompanying molecular-docking studies failed to explain why one of these inhibitors, AM-6, had approximately 2500-fold selectivity for MMP-3 over MMP-2. A number of studies have suggested that the matrix-metalloproteinase active site is highly flexible, leading some to speculate that differences in active-site flexibility may explain inhibitor selectivity. To extend the bioinorganic approach in a way that accounts for MMP-2 and MMP-3 dynamics, we here investigate the predicted binding modes and energies of AM-6 docked into multiple structures extracted from matrix-metalloproteinase molecular dynamics simulations. Our findings suggest that accounting for protein dynamics is essential for the accurate prediction of binding affinity and selectivity. Additionally, AM-6 and other similar inhibitors likely select for and stabilize only a subpopulation of all matrix-metalloproteinase conformations sampled by the apo protein. Consequently, when attempting to predict ligand affinity and selectivity using an ensemble of protein structures, it may be wise to disregard protein conformations that cannot accommodate the ligand.

The nonmevalonate pathway is responsible for isoprenoid production in microbes, including H. pylori, M. tuberculosis and P. falciparum, but is nonexistent in humans, thus providing a desirable route for antibacterial and antimalarial drug discovery. We coordinate a structural study of IspH, a [4Fe-4S] protein responsible for converting HMBPP to IPP and DMAPP in the ultimate step in the nonmevalonate pathway. By performing accelerated molecular dynamics simulations on both substrate-free and HMBPP-bound [Fe4S4](2+) IspH, we elucidate how substrate binding alters the dynamics of the protein. Using principal component analysis, we note that while substrate-free IspH samples various open and closed conformations, the closed conformation observed experimentally for HMBPP-bound IspH is inaccessible in the absence of HMBPP. In contrast, simulations with HMBPP bound are restricted from accessing the open states sampled by the substrate-free simulations. Further investigation of the substrate-free simulations reveals large fluctuations in the HMBPP binding pocket, as well as allosteric pocket openings - both of which are achieved through the hinge motions of the individual domains in IspH. Coupling these findings with solvent mapping and various structural analyses reveals alternative druggable sites that may be exploited in future drug design efforts.