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.

Molecular Modeling Suggests Cruzain Specificity for Peptide Primaquine Prodrugs

Chagas` disease, infection caused by the protozoan Trypanosoma cruzi, is an important, social and medical ailment in the Latin America. This disease is endemic in 21 countries, mostly Latin America countries, with more than 300,000 new cases every year and about 16-18 million infected people. Current therapy is not effective in the chronic phase of the disease. Thus, new and better drugs are urgently needed. In this sense, the in vitro activity of primaquine (PQ) was reported. Based on this, peptide prodrugs of primaquine containing dipeptides - lysine-arginine (LysArg), phenylalanine-alanine (PheAla) and phenylalanine-arginine (PheArg) -- as carriers, were designed to be selectively cleaved by cruzain, a specific cysteine protease of T. cruzi. The prodrugs have shown to be active against tripomastigote forms according to this order: LysArg-PQ> PheAla-PQ> PheArg-PQ. The molecular mechanism of action considered a probable nucleophilic attack of the catalytic residue of cruzain (Cys25) on the respective prodrug amide carbonyl carbon, releasing PQ. In order to test this hypothesis, molecular modeling studies were performed, physicochemical parameters and stereoelectronic features calculated by using the AM1 semi-empirical method suggest that the amide carbonyl carbon is favorable for cleavage, where the LysArg showed the most electronic reactive and sterically disposable, leading to the prodrug release and action. In addition, the docking study indicates the occurrence of specific interactions between prodrugs and the pockets S1 and S2 of cruzain through the dipeptides carriers, being the distance between cruzain Cys25 and the amide carbonyl group related to the biological activity of the prodrugs.

Mechanism of Heparin Acceleration of Tissue Inhibitor of Metalloproteases-1 (TIMP-1) Degradation by the Human Neutrophil Elastase

Heparin has been shown to regulate human neutrophil elastase (HNE) activity. We have assessed the regulatory effect of heparin on Tissue Inhibitor of Metalloproteases-1 [TIMP-1] hydrolysis by HNE employing the recombinant form of TIMP-1 and correlated FRET-peptides comprising the TIMP-1 cleavage site. Heparin accelerates 2.5-fold TIMP-1 hydrolysis by HNE. The kinetic parameters of this reaction were monitored with the aid of a FRET-peptide substrate that mimics the TIMP-1 cleavage site in pre-steady-state conditionsby using a stopped-flow fluorescence system. The hydrolysis of the FRET-peptide substrate by HNE exhibits a pre-steady-state burst phase followed by a linear, steady-state pseudo-first-order reaction. The HNE acylation step (k(2)=21 +/- 1 s(-1)) was much higher than the HNE deacylation step (k(3)=0.57 +/- 0.05 s(-1)). The presence of heparin induces a dramatic effect in the pre-steady-state behavior of HNE. Heparin induces transient lag phase kinetics in HNE cleavage of the FRET-peptide substrate. The pre-steady-state analysis revealed that heparin affects all steps of the reaction through enhancing the ES complex concentration, increasing k(1) 2.4-fold and reducing k(-1) 3.1-fold. Heparin also promotes a 7.8-fold decrease in the k(2) value, whereas the k(3) value in the presence of heparin was increased 58-fold. These results clearly show that heparin binding accelerates deacylation and slows down acylation. Heparin shifts the HNE pH activity profile to the right, allowing HNE to be active at alkaline pH. Molecular docking and kinetic analysis suggest that heparin induces conformational changes in HNE structure. Here, we are showing for the first time that heparin is able to accelerate the hydrolysis of TIMP-1 by HNE. The degradation of TIMP-1is associated to important physiopathological states involving excessive activation of MMPs.

Sulfonyl-hydrazones of cyclic imides derivatives as potent inhibitors of the Mycobacterium tuberculosis protein tyrosine phosphatase B (PtpB)

Searching lead compounds for new antituberculosis drugs, the activity of synthetic sulfonamides and sulfonyl-hydrazones were assayed for their potential inhibitory activity towards a protein tyrosine phosphatase from Mycobacterium tuberculosis - PtpB. Four sulfonyl-hydrazones N-phenylmaleimide derivatives were active (compounds 14, 15, 19 and 21), and the inhibition of PtpB was found to be competitive with respect to the substrate p-nitrophenyl phosphate. Structure-based molecular docking simulations were performed and indicated that the new inhibitor candidates showed similar binding modes, filling the hydrophobic pocket of the protein by the establishment of van der Waals contacts, thereby contributing significantly to the complex stability.

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.

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.

Studies toward the structural optimization of novel thiazolylhydrazone-based potent antitrypanosomal agents

In previous studies, we identified promising anti-Trypanosoma cruzi cruzain inhibitors based on thiazolylhydrazones. To optimize this series, a number of medicinal chemistry directions were explored and new thiazolylhydrazones and thiosemicarbazones were thus synthesized. Potent cruzain inhibitors were identified, such as thiazolylhydrazones 3b and 3j, which exhibited IC(50) of 200-400 nM. Furthermore, molecular docking studies showed concordance with experimentally derived structure-activity relationships (SAR) data. In the course of this work, lead compounds exhibiting in vitro activity against both the epimastigote and trypomastigote forms of T. cruzi were identified and in vivo general toxicity analysis was subsequently performed. Novel SAR were documented, including the importance of the thiocarbonyl carbon attached to the thiazolyl ring and the direct comparison between thiosemicarbazones and thiazolylhydrazones. (C) 2010 Elsevier Ltd. All rights reserved.

Cloning, expression, molecular modelling and docking analysis of glutathione transferase from Saccharum officinarum

Sugarcane yield and quality are affected by a number of biotic and abiotic stresses. In response to such stresses, plants may increase the activities of some enzymes such as glutathione transferase (GST), which are involved in the detoxification of xenobiotics. Thus, a sugarcane GST was modelled and molecular docked using the program LIGIN to investigate the contributions of the active site residues towards the binding of reduced glutathione (GSH) and 1-chloro-2,4-dinitrobenzene (CDNB). As a result, W13 and I119 were identified as key residues for the specificity of sugarcane GSTF1 (SoGSTF1) towards CDNB. To obtain a better understanding of the catalytic specificity of sugarcane GST (SoGSTF1), two mutants were designed, W13L and I119F. Tertiary structure models and the same docking procedure were performed to explain the interactions between sugarcane GSTs with GSH and CDNB. An electron-sharing network for GSH interaction was also proposed. The SoGSTF1 and the mutated gene constructions were cloned and expressed in Escherichia coli and the expressed protein purified. Kinetic analyses revealed different Km values not only for CDNB, but also for GSH. The Km values were 0.2, 1.3 and 0.3 mM for GSH, and 0.9, 1.2 and 0.5 mM for CDNB, for the wild type, W13L mutant and I119F mutant, respectively. The V(max) values were 297.6, 224.5 and 171.8 mu mol min(-1) mg(-1) protein for GSH, and 372.3, 170.6 and 160.4 mu mol min(-1) mg(-1) protein for CDNB.

Use of virtual screening, flexible docking, and molecular interaction fields to design novel HMG-CoA reductase inhibitors for the treatment of hypercholesterolemia

Dietary changes associated with drug therapy can reduce high serum cholesterol levels and dramatically decrease the risk of coronary artery disease, stroke, and overall mortality. Statins are hypolipemic drugs that are effective in the reduction of cholesterol serum levels, attenuating cholesterol synthesis in liver by competitive inhibition regarding the substrate or molecular target HMG-CoA reductase. We have herewith used computer-aided molecular design tools, i.e., flexible docking, virtual screening in large data bases, molecular interaction fields to propose novel potential HMG-CoA reductase inhibitors that are promising for the treatment of hypercholesterolemia.

Molecular dynamics, flexible docking, virtual screening, ADMET predictions, and molecular interaction field studies to design novel potential MAO-B inhibitors

Monoamine oxidase is a flavoenzyme bound to the mitochondrial outer membranes of the cells, which is responsible for the oxidative deamination of neurotransmitter and dietary amines. It has two distinct isozymic forms, designated MAO-A and MAO-B, each displaying different substrate and inhibitor specificities. They are the well-known targets for antidepressant, Parkinson`s disease, and neuroprotective drugs. Elucidation of the x-ray crystallographic structure of MAO-B has opened the way for the molecular modeling studies. In this work we have used molecular modeling, density functional theory with correlation, virtual screening, flexible docking, molecular dynamics, ADMET predictions, and molecular interaction field studies in order to design new molecules with potential higher selectivity and enzymatic inhibitory activity over MAO-B.

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.

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

Molecular systematics of Thorea (Rhodophyta, Thoreales) species in Brazil

This study aimed to evaluate species level taxonomy and phylogenetic relationship among Thorea species in Brazil and other regions of the world using two molecular markers - RUBISCO large subunit plastid gene (rbcL) and nuclear small-subunit ribosomal DNA (SSU rDNA). Three samples of Thorea from Brazil (states of Mato Grosso do Sul and São Paulo) and one sample from Dominican Republic (DR) were sequenced. Analyses based on partial sequences of rbcL (1,282 bp) and complete sequences of SSU (1,752 bp) were essentially congruent and revealed that Thoreales formed a distinct monophyletic clade, which had two major branches with high support, representing the genera Thorea and Nemalionopsis. Thorea clade had four main branches with high support for all analyses, each one representing the species: 1) T. gaudichaudii C. Agardh from Asia (Japan and Philippines) - this clade occurred only in the rbcL analyses; 2) T. violacea Bory from Asia (Japan) and North America (U.S.A. and DR); 3) T. hispida (Thore) Desvaux from Europe (England) and Asia (Japan); 4) a distinct group with the three Brazilian samples (sequence identity: rbcL 97.2%, 1,246 bp; SSU 96.0-98.1%, 1,699-1,720 bp). The Brazilian samples clearly formed a monophyletic clade based on both molecular markers and was interpreted as a separate species, for which we resurrected the name T. bachmannii Pujals. Morphological and molecular evidences indicate that the Thoreales is well-resolved at ordinal and generic levels. In contrast, Thorea species recognized by molecular data require additional characters (e.g. reproductive and chromosome numbers) to allow consistent and reliable taxonomic circumscription aiming at a world revision based on molecular and morphological evidences.

Thaptomys Thomas 1915 (Rodentia, Sigmodontinae, Akodontini) with karyotypes 2n = 50, FN = 48, and 2n = 52, FN = 52: two monophyletic lineages recovered by molecular phylogeny

A novel karyotype with 2n = 50, FN = 48, was described for specimens of Thaptomys collected at Una, State of Bahia, Brazil, which are morphologically indistinguishable from Thaptomys nigrita, 2n = 52, FN = 52, found in other localities. It was hence proposed that the 2n = 50 karyotype could belong to a distinct species, cryptic of Thaptomys nigrita, once chromosomal rearrangements observed, along with the geographic distance, might represent a reproductive barrier between both forms. Phylogenetic analyses using maximum parsimony and maximum likelihood based on partial cytochrome b sequences with 1077 bp were performed, attempting to establish the relationships among the individuals with distinct karyotypes along the geographic distribution of the genus; the sample comprised 18 karyotyped specimens of Thaptomys, encompassing 15 haplotypes, from eight different localities of the Atlantic Rainforest. The intra-generic relationships corroborated the distinct diploid numbers, once both phylogenetic reconstructions recovered two monophyletic lineages, a northeastern clade grouping the 2n = 50 and a southeastern clade with three subclades, grouping the 2n = 52 karyotype. The sequence divergence observed between their individuals ranged from 1.9% to 3.5%.