'Unconventional' coordination chemistry by metal chelating fragments in a metalloprotein active site.

The binding of three closely related chelators: 5-hydroxy-2-methyl-4H-pyran-4-thione (allothiomaltol, ATM), 3-hydroxy-2-methyl-4H-pyran-4-thione (thiomaltol, TM), and 3-hydroxy-4H-pyran-4-thione (thiopyromeconic acid, TPMA) to the active site of human carbonic anhydrase II (hCAII) has been investigated. Two of these ligands display a monodentate mode of coordination to the active site Zn(2+) ion in hCAII that is not recapitulated in model complexes of the enzyme active site. This unprecedented binding mode in the hCAII-thiomaltol complex has been characterized by both X-ray crystallography and X-ray spectroscopy. In addition, the steric restrictions of the active site force the ligands into a 'flattened' mode of coordination compared with inorganic model complexes. This change in geometry has been shown by density functional computations to significantly decrease the strength of the metal-ligand binding. Collectively, these data demonstrate that the mode of binding by small metal-binding groups can be significantly influenced by the protein active site. Diminishing the strength of the metal-ligand bond results in unconventional modes of metal coordination not found in typical coordination compounds or even carefully engineered active site models, and understanding these effects is critical to the rational design of inhibitors that target clinically relevant metalloproteins.

Utilizing a dynamical description of IspH to aid in the development of novel antimicrobial drugs.

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.

Guide to virtual screening: application to the Akt phosphatase PHLPP.

We present an example-based description of virtual screening (VS) techniques used to identify new regulators of the Akt phosphatase PHLPP (PH domain Leucine repeat Protein Phosphatase). This enzyme opposes the effects of two kinases, Akt and PKC, which play a major role in cell growth and survival. Therefore, PHLPP is a potential therapeutic target in pathophysiologies where these pathways are either repressed, such as in diabetes and cardiovascular diseases, or over-activated as in cancer. To the best of our knowledge, no PHLPP inhibitors have been reported so far in the literature. In this study, we used a combination of chemical and virtual screening techniques that led to the identification of a number of inhibiting compounds with diverse scaffolds. These compounds bind PHLPP and inhibit cell death when tested in cellular assays. We employed GLIDE docking software to screen a library of more than 40,000 compounds selected from the NCI open depository (250,000 compounds) by similarity searches. We compare the efficiency at which we determined binding compounds from the chemical screen, and compare enrichment factors of the virtually discovered compounds over chemical screening.

Pyrone-based inhibitors of metalloproteinase types 2 and 3 may work as conformation-selective inhibitors.

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.

E9-Im9 colicin DNase-immunity protein biomolecular association in water: a multiple-copy and accelerated molecular dynamics simulation study.

Protein-protein transient and dynamic interactions underlie all biological processes. The molecular dynamics (MD) of the E9 colicin DNase protein, its Im9 inhibitor protein, and their E9-Im9 recognition complex are investigated by combining multiple-copy (MC) MD and accelerated MD (aMD) explicit-solvent simulation approaches, after validation with crystalline-phase and solution experiments. Im9 shows higher flexibility than its E9 counterpart. Im9 displays a significant reduction of backbone flexibility and a remarkable increase in motional correlation upon E9 association. Im9 loops 23-31 and 54-64 open with respect to the E9-Im9 X-ray structure and show high conformational diversity. Upon association a large fraction (approximately 20 nm2) of E9 and Im9 protein surfaces become inaccessible to water. Numerous salt bridges transiently occurring throughout our six 50 ns long MC-MD simulations are not present in the X-ray model. Among these Im9 Glu31-E9 Arg96 and Im9 Glu41-Lys89 involve interface interactions. Through the use of 10 ns of Im9 aMD simulation, we reconcile the largest thermodynamic impact measured for Asp51Ala mutation with Im9 structure and dynamics. Lys57 acts as an essential molecular switch to shift Im9 surface loop towards an ideal configuration for E9 inhibition. This is achieved by switching Asp60-Lys57 and Asp62-Lys57 hydrogen bonds to Asp51-Lys57 salt bridge. E9-Im9 recognition involves shifts of conformational distributions, reorganization of intramolecular hydrogen bond patterns, and formation of new inter- and intramolecular interactions. The description of key transient biological interactions can be significantly enriched by the dynamic and atomic-level information provided by computer simulations.