Structural basis for selective inhibition of trypanosomatid glyceraldehyde-3-phosphate dehydrogenase: Molecular docking and 3D QSAR studies

The glycolytic enzyme glyceraldehyde-3 -phosphate dehydrogenase (GAPDH) is as an attractive target for the development of novel antitrypanosomatid agents. In the present work, comparative molecular field analysis and comparative molecular similarity index analysis were conducted on a large series of selective inhibitors of trypanosomatid GAPDH. Four statistically significant models were obtained (r(2) > 0.90 and q(2) > 0.70), indicating their predictive ability for untested compounds. The models were then used to predict the potency of an external test set, and the predicted values were in good agreement with the experimental results. Molecular modeling studies provided further insight into the structural basis for selective inhibition of trypanosomatid GAPDH.

Docking and molecular dynamics simulation of quinone compounds with trypanocidal activity

In this work, two different docking programs were used, AutoDock and FlexX, which use different types of scoring functions and searching methods. The docking poses of all quinone compounds studied stayed in the same region in the trypanothione reductase. This region is a hydrophobic pocket near to Phe396, Pro398 and Leu399 amino acid residues. The compounds studied displays a higher affinity in trypanothione reductase (TR) than glutathione reductase (GR), since only two out of 28 quinone compounds presented more favorable docking energy in the site of human enzyme. The interaction of quinone compounds with the TR enzyme is in agreement with other studies, which showed different binding sites from the ones formed by cysteines 52 and 58. To verify the results obtained by docking, we carried out a molecular dynamics simulation with the compounds that presented the highest and lowest docking energies. The results showed that the root mean square deviation (RMSD) between the initial and final pose were very small. In addition, the hydrogen bond pattern was conserved along the simulation. In the parasite enzyme, the amino acid residues Leu399, Met400 and Lys402 are replaced in the human enzyme by Met406, Tyr407 and Ala409, respectively. In view of the fact that Leu399 is an amino acid of the Z site, this difference could be explored to design selective inhibitors of TR.

Role of Halogen Bonds in Thyroid Hormone Receptor Selectivity: Pharmacophore-Based 3D-QSSR Studies

Most physiological effects of thyroid hormones are mediated by the two thyroid hormone receptor subtypes, TR alpha and TR beta. Several pharmacological effects mediated by TR beta might be beneficial in important medical conditions such as obesity, hypercholesterolemia and diabetes, and selective TR beta activation may elicit these effects while maintaining an acceptable safety profile, To understand the molecular determinants of affinity and subtype selectivity of TR ligands, we have successfully employed a ligand- and structure-guided pharmacophore-based approach to obtain the molecular alignment of a large series of thyromimetics. Statistically reliable three-dimensional quantitative structure-activity relationship (3D-QSAR) and three-dimensional quantitative structure-selectivity relationship (3D-QSSR) models were obtained using the comparative molecular field analysis (CoMFA) method, and the visual analyses of the contour maps drew attention to a number of possible opportunities for the development of analogs with improved affinity and selectivity. Furthermore, the 3D-QSSR analysis allowed the identification of a novel and previously unmentioned halogen bond, bringing new insights to the mechanism of activity and selectivity of thyromimetics.

Convergent synthesis and cruzain inhibitory activity of novel 2-(N `-benzylidenehydrazino)-4-trifluoromethyl-pyrimidines

To search for new cruzain inhibitors, the synthesis of a series of novel 2-(N`-benzylidenehydrazino)-4-trifluoromethyl-pyrimidines in a convergent manner is presented. The cruzain inhibitory activity of some of these compounds was evaluated and a binding model was proposed. All derivatives tested were active and the most significant inhibitory effect (80% at 100 mu M) and IC(50) value (85 mu M) were obtained from the 2-(N`-4-chloro-benzylidenehydrazino)-4-trifluoromethyl-pyrimidine. Although further structural optimization to improve solubility is necessary, the molecular docking studies suggest that these inhibitors occupy the S2 pocket without irreversible enzyme inactivation, through hydrophobic interactions, thus, indicating a desirable mode of interaction for the design of novel inhibitors. (C) 2008 Elsevier Ltd. All rights reserved.

In silico screening of HIV-1 non-nucleoside reverse transcriptase and protease inhibitors

Two targets, reverse transcriptase (RT) and protease from HIV-1, were used during the past two decades to the discovery of non-nucleoside reverse transcriptase inhibitors (NNRTI) and protease inhibitors (PI) that belong to the arsenal of the antiretroviral therapy. Herein these enzymes were chosen as templates for conducting a computer-aided ligand design. Ligand and structure-based drug designs were the starting points to select compounds from a database bearing more than five million compounds by means of cheminformatic tools. New promising lead structures are retrieved from the database, which are open to acquisition and test. Classes of molecules already described as NNRTI or PI in the literature also came out and were useful to prove the reliability of the workflow, and thus validating the work carried out so far. (c) 2007 Elsevier Masson SAS. All rights reserved.

Structural modeling of high-affinity thyroid receptor-ligand complexes

Understanding the molecular basis of the binding modes of natural and synthetic ligands to nuclear receptors is fundamental to our comprehension of the activation mechanism of this important class of hormone regulated transcription factors and to the development of new ligands. Thyroid hormone receptors (TR) are particularly important targets for pharmaceuticals development because TRs are associated with the regulation of metabolic rates, body weight, and circulating levels of cholesterol and triglycerides in humans. While several high-affinity ligands are known, structural information is only partially available. In this work we obtain structural models of several TR-ligand complexes with unknown structure by docking high affinity ligands to the receptors` ligand binding domain with subsequent relaxation by molecular dynamics simulations. The binding modes of these ligands are discussed providing novel insights into the development of TR ligands. The experimental binding free energies are reasonably well-reproduced from the proposed models using a simple linear interaction energy free-energy calculation scheme.

Fragment-guided approach to incorporating structural information into a CoMFA study: BACE-1 as an example

Alzheimer`s disease is an ultimately fatal neurodegenerative disease, and BACE-1 has become an attractive validated target for its therapy, with more than a hundred crystal structures deposited in the PDB. In the present study, we present a new methodology that integrates ligand-based methods with structural information derived from the receptor. 128 BACE-1 inhibitors recently disclosed by GlaxoSmithKline R&D were selected specifically because the crystal structures of 9 of these compounds complexed to BACE-1, as well as five closely related analogs, have been made available. A new fragment-guided approach was designed to incorporate this wealth of structural information into a CoMFA study, and the methodology was systematically compared to other popular approaches, such as docking, for generating a molecular alignment. The influence of the partial charges calculation method was also analyzed. Several consistent and predictive models are reported, including one with r (2) = 0.88, q (2) = 0.69 and r (pred) (2) = 0.72. The models obtained with the new methodology performed consistently better than those obtained by other methodologies, particularly in terms of external predictive power. The visual analyses of the contour maps in the context of the enzyme drew attention to a number of possible opportunities for the development of analogs with improved potency. These results suggest that 3D-QSAR studies may benefit from the additional structural information added by the presented methodology.

The Guarguaraz Complex and the Neoproterozoic-Cambrian evolution of southwestern Gondwana: Geochemical signatures and geochronological constraints

The Guarguaraz Complex, in western Argentina, comprises a metasedimentary assemblage, associated with mafic sills and ultramafic bodies intruded by basaltic dikes, which are interpreted as Ordovician dismembered ophiolites. Two kinds of dikes are recognized, a group associated with the metasediments and the other ophiolite-related. Both have N-MORB signatures, with epsilon(Nd) between +3.5 and +8.2, indicating a depleted source, and Grenville model ages between 0.99 and 1.62 Ga. A whole-rock Sm-Nd isochron yielded an age of 655 +/- 76 Ma for these mafic rocks, which is compatible with cianobacteria and acritarchae recognized in the clastic metasedimentary platform sequences, that indicate a Neoproterozoic (Vendian)-Cambrian age of deposition. The Guarguaraz metasedimentary-ophiolitic complex represents, therefore, a remnant of an oceanic basin developed to the west of the Grenville-aged Cuyania terrane during the Neoproterozoic. The southernmost extension of these metasedimentary sequences in Cordon del Portillo might represent part of this platform and not fragments of the Chilenia terrane. An extensional event related to the fragmentation of Rodinia is represented by the mafic and ultramafic rocks. The Devonian docking of Chilenia emplaced remnants of ocean floor and slices of the Cuyania terrane (Las Yaretas Gneisses) in tectonic contact with the Neoproterozoic metasediments, marking the Devonian western border of Gondwana. (C) 2009 Elsevier Ltd. All rights reserved.

PACKMOL: A Package for Building Initial Configurations for Molecular Dynamics Simulations

Adequate initial configurations for molecular dynamics simulations consist of arrangements of molecules distributed in space in such a way to approximately represent the system`s overall structure. In order that the simulations are not disrupted by large van der Waals repulsive interactions, atoms from different molecules Must keep safe pairwise distances. Obtaining Such a molecular arrangement can be considered it packing problem: Each type molecule must satisfy spatial constraints related to the geometry of the system, and the distance between atoms of different molecules Must be greater than some specified tolerance. We have developed a code able to pack millions of atoms. grouped in arbitrarily complex molecules, inside a variety of three-dimensional regions. The regions may be intersections of spheres, ellipses, cylinders, planes, or boxes. The user must provide only the structure of one molecule of each type and the geometrical constraints that each type of molecule must satisfy. Building complex mixtures, interfaces, solvating biomolecules in water, other solvents, or mixtures of solvents, is straight forward. In addition. different atoms belonging to the same molecule may also be restricted to different spatial regions, in Such a way that more ordered molecular arrangements call be built, as micelles. lipid double-layers, etc. The packing time for state-of-the-art molecular dynamics systems varies front a few seconds to a few Minutes in a personal Computer. The input files are simple and Currently compatible with PDB, Tinker, Molden, or Moldy coordinate files. The package is distributed as free software and call be downloaded front http://www.ime.unicamp.br/similar to martinez/packmol/. (C) 2009 Wiley Periodicals. Inc. J Comput Chem 30: 2157-2164, 2009

Deciphering the role of the electrostatic interactions in the alpha-tropomyosin head-to-tail complex

Skeletal alpha-tropomyosin (Tm) is a dimeric coiled-coil protein that forms linear assemblies under low ionic strength conditions in vitro through head-to-tail interactions. A previously published NMR structure of the Tin head-to-tail complex revealed that it is formed by the insertion of the N-terminal coiled-coil of one molecule into a cleft formed by the separation of the helices at the C-terminus of a second molecule. To evaluate the contribution of charged residues to complex stability, we employed single and double-mutant Tm fragments in which specific charged residues were changed to alanine in head-to-tail binding assays, and the effects of the mutations were analyzed by thermodynamic double-mutant cycles and protein-protein docking. The results show that residues K5, K7, and D280 are essential to the stability of the complex. Though D2, K6, D275, and H276 are exposed to the solvent and do not participate in intermolecular contacts in the NMR structure, they may contribute to head-to-tail complex stability by modulating the stability of the helices at the Tm termini.

MPO Inhibitors Selected by Virtual Screening

The hemeprotein myeloperoxidase (MPO) participates in innate immune defense through its ability to generate potent microbicidal oxidants. However, these oxidants are also key mediators of the tissue damage associated with many inflammatory diseases. Thus, there is considerable interest in developing therapeutically useful MPO inhibitors. Here, we used structure-based drug design (SBDD) and ligand-based drug design (LBDD) to select for potentially new and selective MPO inhibitors. A pharmacophore model was developed based on the crystal structure of human MPO in complex with salicylhydroxamic acid (SHA), a known inhibitor of the enzyme. The pharmacophore model was used to screen the ZINC database for potential ligands, which were further filtered on the basis of their physical-chemical properties and docking score. The filtered compounds were visually inspected, and nine were purchased for experimental studies. Surprisingly, almost all of the selected compounds belonged to the aromatic hydrazide class, which had been previously described as MPO inhibitors. The compounds selected by virtual screening were shown to inhibit the chlorinating activity of MPO; the top four compounds displayed IC(50) values ranging from 1.0 to 2.8 mM. MPO inactivation by the most effective compound was shown to be irreversible. Overall, our results show that SBDD and LBDD may be useful for the rational development of new MPO inhibitors.

Flexibility and inhibitor binding in Cdc25 phosphatases

Cdc25 phosphatases involved in cell cycle checkpoints are now active targets for the development of anti-cancer therapies. Rational drug design would certainly benefit from detailed structural information for Cdc25s. However, only apo- or sulfate-bound crystal structures of the Cdc25 catalytic domain have been described so far. Together with previously available crystalographic data, results from molecular dynamics simulations, bioinformatic analysis, and computer-generated conformational ensembles shown here indicate that the last 30-40 residues in the C-terminus of Cdc25B are partially unfolded or disordered in solution. The effect of C-terminal flexibility upon binding of two potent small molecule inhibitors to Cdc25B is then analyzed by using three structural models with variable levels of flexibility, including an equilibrium distributed ensemble of Cdc25B backbone conformations. The three Cdc25B structural models are used in combination with flexible docking, clustering, and calculation of binding free energies by the linear interaction energy approximation to construct and validate Cdc25B-inhibitor complexes. Two binding sites are identified on top and beside the Cdc25B active site. The diversity of interaction modes found increases with receptor flexibility. Backbone flexibility allows the formation of transient cavities or compact hydrophobic units on the surface of the stable, folded protein core that are unexposed or unavailable for ligand binding in rigid and densely packed crystal structures. The present results may help to speculate on the mechanisms of small molecule complexation to partially unfolded or locally disordered proteins.

Integration of Ligand- and Target-Based Virtual Screening for the Discovery of Cruzain Inhibitors

A myriad of methods are available for virtual screening of small organic compound databases. In this study we have successfully applied a quantitative model of consensus measurements, using a combination of 3D similarity searches (ROCS and EON), Hologram Quantitative Structure Activity Relationships (HQSAR) and docking (FRED, FlexX, Glide and AutoDock Vina), to retrieve cruzain inhibitors from collected databases. All methods were assessed individually and then combined in a Ligand-Based Virtual Screening (LBVS) and Target-Based Virtual Screening (TBVS) consensus scoring, using Receiving Operating Characteristic (ROC) curves to evaluate their performance. Three consensus strategies were used: scaled-rank-by-number, rank-by-rank and rank-by-vote, with the most thriving the scaled-rank-by-number strategy, considering that the stiff ROC curve appeared to be satisfactory in every way to indicate a higher enrichment power at early retrieval of active compounds from the database. The ligand-based method provided access to a robust and predictive HQSAR model that was developed to show superior discrimination between active and inactive compounds, which was also better than ROCS and EON procedures. Overall, the integration of fast computational techniques based on ligand and target structures resulted in a more efficient retrieval of cruzain inhibitors with desired pharmacological profiles that may be useful to advance the discovery of new trypanocidal agents.

Stoichiometry and thermodynamics of the interaction between the C-terminus of human 90 kDa heat shock protein Hsp90 and the mitochondrial translocase of outer membrane Tom70

A large majority of the 1000-1500 proteins in the mitochondria are encoded by the nuclear genome, and therefore, they are translated in the cytosol in the form and contain signals to enable the import of proteins into the organelle. The TOM complex is the major translocase of the outer membrane responsible for preprotein translocation. It consists of a general import pore complex and two membrane import receptors, Tom20 and Tom70. Tom70 contains a characteristic TPR domain, which is a docking site for the Hsp70 and Hsp90 chaperones. These chaperones are involved in protecting cytosolic preproteins from aggregation and then in delivering them to the TOM complex. Although highly significant, many aspects of the interaction between Tom70 and Hsp90 are still uncertain. Thus, we used biophysical tools to study the interaction between the C-terminal domain of Hsp90 (C-Hsp90), which contains the EEVD motif that binds to TPR domains, and the cytosolic fragment of Tom70. The results indicate a stoichiometry of binding of one monomer of Tom70 per dimer of C-Hsp90 with a K(D) of 360 30 nM, and the stoichiometry and thermodynamic parameters obtained suggested that Tom70 presents a different mechanism of interaction with Hsp90 when compared with other TPR proteins investigated. (C) 2011 Elsevier Inc. All rights reserved.

Novel insights for dihydroorotate dehydrogenase class 1A inhibitors discovery

The enzyme dihydroorotate dehydrogenase (DHODH) has been suggested as a promising target for the design of trypanocidal agents. We report here the discovery of novel inhibitors of Trypanosoma cruzi DHODH identified by a combination of virtual screening and ITC methods. Monitoring of the enzymatic reaction in the presence of selected ligands together with structural information obtained from X-ray crystallography analysis have allowed the identification and validation of a novel site of interaction (S2 site). This has provided important structural insights for the rational design of T cruzi and Leishmania major DHODH inhibitors. The most potent compound (1) in the investigated series inhibits TcDHODH enzyme with K(i)(app) value of 19.28 mu M and possesses a ligand efficiency of 0.54 kcal mol(-1) per non-H atom. The compounds described in this work are promising hits for further development. (C) 2010 Elsevier Masson SAS. All rights reserved.