Novel aryl beta-aminocarbonyl derivatives as inhibitors of Trypanosoma cruzi trypanothione reductase: binding mode revised by docking and GRIND2-based 3D-QSAR procedures

Trypanothione reductase has long been investigated as a promising target for chemotherapeutic intervention in Chagas disease, since it is an enzyme of a unique metabolic pathway that is exclusively present in the pathogen but not in the human host, which has the analog Glutathione reductase. In spite of the present data-set includes a small number of compounds, a combined use of flexible docking, pharmacophore perception, ligand binding site prediction, and Grid-Independent Descriptors GRIND2-based 3D-Quantitative Structure-Activity Relationships (QSAR) procedures allowed us to rationalize the different biological activities of a series of 11 aryl beta-aminocarbonyl derivatives, which are inhibitors of Trypanosoma cruzi trypanothione reductase (TcTR). Three QSAR models were built and validated using different alignments, which are based on docking with the TcTR crystal structure, pharmacophore, and molecular interaction fields. The high statistical significance of the models thus obtained assures the robustness of this second generation of GRIND descriptors here used, which were able to detect the most important residues of such enzyme for binding the aryl beta-aminocarbonyl derivatives, besides to rationalize distances among them. Finally, a revised binding mode has been proposed for our inhibitors and independently supported by the different methodologies here used, allowing further optimization of the lead compounds with such combined structure- and ligand-based approaches in the fight against the Chagas disease.

Automatic design of decision-tree induction algorithms tailored to flexible-receptor docking data

Background: This paper addresses the prediction of the free energy of binding of a drug candidate with enzyme InhA associated with Mycobacterium tuberculosis. This problem is found within rational drug design, where interactions between drug candidates and target proteins are verified through molecular docking simulations. In this application, it is important not only to correctly predict the free energy of binding, but also to provide a comprehensible model that could be validated by a domain specialist. Decision-tree induction algorithms have been successfully used in drug-design related applications, specially considering that decision trees are simple to understand, interpret, and validate. There are several decision-tree induction algorithms available for general-use, but each one has a bias that makes it more suitable for a particular data distribution. In this article, we propose and investigate the automatic design of decision-tree induction algorithms tailored to particular drug-enzyme binding data sets. We investigate the performance of our new method for evaluating binding conformations of different drug candidates to InhA, and we analyze our findings with respect to decision tree accuracy, comprehensibility, and biological relevance. Results: The empirical analysis indicates that our method is capable of automatically generating decision-tree induction algorithms that significantly outperform the traditional C4.5 algorithm with respect to both accuracy and comprehensibility. In addition, we provide the biological interpretation of the rules generated by our approach, reinforcing the importance of comprehensible predictive models in this particular bioinformatics application. Conclusions: We conclude that automatically designing a decision-tree algorithm tailored to molecular docking data is a promising alternative for the prediction of the free energy from the binding of a drug candidate with a flexible-receptor.

Synthesis, anti-inflammatory activity and molecular docking studies of 2,5-diarylfuran amino acid derivatives

A series of 2,5-diaryl substituted furans functionalized with several amino acids were synthesized and evaluated as the cyclooxygenases COX-1 and COX-2 enzymes inhibitors. The proline-substituted compound inhibited PGE(2) secretion by LPS-stimulated neutrophils, suggesting selectivity for COX-2. Molecular docking studies in the binding site of COX-2 were performed. (C) 2011 Elsevier Masson SAS. All rights reserved.

Synthesis, biological evaluation and molecular docking studies of 3-(triazolyl)coumarin derivatives: Effect on inducible nitric oxide synthase

A series of 3-(triazolyl)-coumarins were synthesized and tested as anti-inflammatory agents. It was possible to infer that these compounds do not alter the interaction of LPS with TLR-4 or TLR-2, as the intracellular pathways involved in the TNF-alpha secretion and COX-2 activity were not affected. Nevertheless, the compounds inhibited iNOS-derived NO production, without affecting the eNOS activity. The outcome of the docking studies showed that it pi center dot center dot center dot pi interactions with the heme group are important for the iNOS inhibition, thus making compound 3c a promising lead. Moreover, the efficacy of this compound was visualized by the reduced number of neutrophils in the LPS-inflamed subcutaneous tissue. Together, biological and docking data show that triazolyl-substituted coumarins, that can act on iNOS, are a good scaffold to be explored. (C) 2012 Elsevier Masson SAS. All rights reserved.

A tellurium-based cathepsin B inhibitor: Molecular structure, modelling, molecular docking and biological evaluation

The crystallographically determined structure of biologically active 4,4-dichloro-1,3-diphenyl-4-telluraoct-2-en-1-one, 3, shows the coordination geometry for Te to be distorted psi-pentagonal bipyramidal based on a C2OCl3(lone pair) donor set. Notable is the presence of an intramolecular axial Te center dot center dot center dot O (carbonyl) interaction, a design element included to reduce hydrolysis. Raman and molecular modelling studies indicate the persistence of the Te center dot center dot center dot O(carbonyl) interaction in the solution (CHCl3) and gasphases, respectively. Docking studies of 3' (i.e. original 3 less one chloride) with Cathepsin B reveals a change in the configuration about the vinyl C = C bond. i.e. to E from Z (crystal structure). This isomerism allows the optimisation of interactions in the complex which features a covalent Te-SGCys29 bond. Crucially, the E configuration observed for 3' allows for the formation of a hypervalent Te center dot center dot center dot O interaction as well as an O center dot center dot center dot H-O hydrogen bond with the Gly27 and Glu122 residues, respectively. Additional stabilisation is afforded by a combination of interactions spanning the S1, S2, S1' and S2' sub-sites of Cathepsin B. The greater experimental inhibitory activity of 3 compared with analogues is rationalised by the additional interactions formed between 3' and the His110 and His111 residues in the occluding loop, which serve to hinder the entrance to the active site. (C) 2012 Elsevier B.V. All rights reserved.

Docking, synthesis and pharmacological activity of novel urea-derivatives designed as p38 MAPK inhibitors

p38 mitogen-activated protein kinase (p38 MAPK) is an important signal transducing enzyme involved in many cellular regulations, including signaling pathways, pain and inflammation. Several p38 MAPK inhibitors have been developed as drug candidates to treatment of autoimmune disorders, such as rheumatoid arthritis. In this paper we reported the docking, synthesis and pharmacological activity of novel urea-derivatives (4a-e) designed as p38 MAPK inhibitors. These derivatives presented good theoretical affinity to the target p38 MAPK, standing out compound 4e (LASSBio-998), which showed a better score value compared to the prototype GK-00687. This compound was able to reduce in vitro TNF-alpha production and was orally active in a hypernociceptive murine model sensible to p38 MAPK inhibitors. Otherwise, compound 4e presented a dose-dependent analgesic effect in a model of antigen (mBSA)-induced arthritis and anti-inflammatory profile in carrageenan induced paw edema, indicating its potential as a new antiarthritis prototype. (c) 2012 Elsevier Masson SAS. All rights reserved.

Structural insights regarding an insecticidal Talisia esculenta protein and its biotechnological potential for Diatraea saccharalis larval control

Talisin is a seed-storage protein from Talisia esculenta that presents lectin-like activities, as well as proteinase-inhibitor properties. The present study aims to provide new in vitro and in silico biochemical information about this protein, shedding some light on its mechanistic inhibitory strategies. A theoretical three-dimensional structure of Talisin bound to trypsin was constructed in order to determine the relative interaction mode. Since the structure of non-competitive inhibition has not been elucidated, Talisin-trypsin docking was carried out using Hex v5.1, since the structure of non-competitive inhibition has not been elucidated. The predicted non-coincidence of the trypsin binding site is completely different from that previously proposed for Kunitz-type inhibitors, which demonstrate a substitution of an Arg(64) for the Glu(64) residue. Data, therefore, provide more information regarding the mechanisms of non-competitive plant proteinase inhibitors. Bioassays with Talisin also presented a strong insecticide effect on the larval development of Diatraea saccharalis, demonstrating LD50 and ED50 of ca. 2.0% and 1.5%, respectively. (C) 2011 Elsevier Inc. All rights reserved.

The Cysteine-Rich Protein Thimet Oligopeptidase as a Model of the Structural Requirements for S-glutathiolation and Oxidative Oligomerization

Thimet oligopeptidase (EP24.15) is a cysteine-rich metallopeptidase containing fifteen Cys residues and no intra-protein disulfide bonds. Previous work on this enzyme revealed that the oxidative oligomerization of EP24.15 is triggered by S-glutathiolation at physiological GSSG levels (10-50 mu M) via a mechanism based on thiol-disulfide exchange. In the present work, our aim was to identify EP24.15 Cys residues that are prone to S-glutathiolation and to determine which structural features in the cysteinyl bulk are responsible for the formation of mixed disulfides through the reaction with GSSG and, in this particular case, the Cys residues within EP24.15 that favor either S-glutathiolation or inter-protein thiol-disulfide exchange. These studies were conducted by in silico structural analyses and simulations as well as site-specific mutation. S-glutathiolation was determined by mass spectrometric analyses and western blotting with anti-glutathione antibody. The results indicated that the stabilization of a thiolate sulfhydryl and the solvent accessibility of the cysteines are necessary for S-thiolation. The Solvent Access Surface analysis of the Cys residues prone to glutathione modification showed that the S-glutathiolated Cys residues are located inside pockets where the sulfur atom comes into contact with the solvent and that the positively charged amino acids are directed toward these Cys residues. The simulation of a covalent glutathione docking onto the same Cys residues allowed for perfect glutathione posing. A mutation of the Arg residue 263 that forms a saline bridge to the Cys residue 175 significantly decreased the overall S-glutathiolation and oligomerization of EP24.15. The present results show for the first time the structural requirements for protein S-glutathiolation by GSSG and are consistent with our previous hypothesis that EP24.15 oligomerization is dependent on the electron transfer from specific protonated Cys residues of one molecule to previously S-glutathionylated Cys residues of another one.

Peptide fingerprinting of the neurotoxic fractions isolated from the secretions of sea anemones Stichodactyla helianthus and Bunodosoma granulifera. New members of the APETx-like family identified by a 454 pyrosequencing approach

Sea anemones are known to contain a wide diversity of biologically active peptides, mostly unexplored according to recent peptidomic and transcriptomic studies. In the present work, the neurotoxic fractions from the exudates of Stichodactyla helianthus and Bunodosoma granulifera were analyzed by reversed-phase chromatography and mass spectrometry. The first peptide fingerprints of these sea anemones were assessed, revealing the largest number of peptide components (156) so far found in sea anemone species, as well as the richer peptide diversity of B. granulifera in relation to S. helianthus. The transcriptomic analysis of B. granulifera, performed by massive cDNA sequencing with 454 pyrosequencing approach allowed the discovery of five new APETx-like peptides (U-AITX-Bg1a-e - including the full sequences of their precursors for four of them), which together with type 1 sea anemone sodium channel toxins constitute a very distinguishable feature of studied sea anemone species belonging to genus Bunodosoma. The molecular modeling of these new APETx-like peptides showed a distribution of positively charged and aromatic residues in putative contact surfaces as observed in other animal toxins. On the other hand, they also showed variable electrostatic potentials, thus suggesting a docking onto their targeted channels in different spatial orientations. Moreover several crab paralyzing toxins (other than U-AITX-Bg1a-e), which induce a variety of symptoms in crabs, were isolated. Some of them presumably belong to new classes of crab-paralyzing peptide toxins, especially those with molecular masses below 2 kDa, which represent the smallest peptide toxins found in sea anemones. (C) 2011 Elsevier Inc. All rights reserved.

Advanced QSAR Studies on PPAR delta Ligands Related to Metabolic Diseases

PPAR delta is a nuclear receptor that, when activated, regulates the metabolism of carbohydrates and lipids and is related to metabolic syndrome and type 2 diabetes. To understand the main interactions between ligands and PPAR delta, we have constructed 2D and 3D QSAR models and compared them with HOMO, LUMO and electrostatic potential maps of the compounds studied, as well as docking results. All QSAR models showed good statistical parameters and prediction outcomes. The QSAR models were used to predict the biological activity of an external test set, and the predicted values are in good agreement with the experimental results. Furthermore, we employed all maps to evaluate the possible interactions between the ligands and PPAR delta. These predictive QSAR models, along with the HOMO, LUMO and MEP maps, can provide insights into the structural and chemical properties that are needed in the design of new PPAR delta ligands that have improved biological activity and can be employed to treat metabolic diseases.

FERM domain interaction with myosin negatively regulates FAK in cardiomyocyte hypertrophy

Focal adhesion kinase (FAK) regulates cellular processes that affect several aspects of development and disease. The FAK N-terminal FERM (4.1 protein-ezrin-radixin-moesin homology) domain, a compact clover-leaf structure, binds partner proteins and mediates intramolecular regulatory interactions. Combined chemical cross-linking coupled to MS, small-angle X-ray scattering, computational docking and mutational analyses showed that the FAK FERM domain has a molecular cleft (similar to 998 angstrom(2)) that interacts with sarcomeric myosin, resulting in FAK inhibition. Accordingly, mutations in a unique short amino acid sequence of the FERM myosin cleft, FP-1, impaired the interaction with myosin and enhanced FAK activity in cardiomyocytes. An FP-1 decoy peptide selectively inhibited myosin interaction and increased FAK activity, promoting cardiomyocyte hypertrophy through activation of the AKT-mammalian target of rapamycin pathway. Our findings uncover an inhibitory interaction between the FAK FERM domain and sarcomeric myosin that presents potential opportunities to modulate the cardiac hypertrophic response through changes in FAK activity.

Ligand- and Structure-Based Drug Design Strategies and PPAR delta/alpha Selectivity

Peroxisome-proliferator-activated receptors are a class of nuclear receptors with three subtypes: a, ? and d. Their main function is regulating gene transcription related to lipid and carbohydrate metabolism. Currently, there are no peroxisome-proliferator-activated receptors d drugs being marketed. In this work, we studied a data set of 70 compounds with a and d activity. Three partial least square models were created, and molecular docking studies were performed to understand the main reasons for peroxisome-proliferator-activated receptors d selectivity. The obtained results showed that some molecular descriptors (log P, hydration energy, steric and polar properties) are related to the main interactions that can direct ligands to a particular peroxisome-proliferator-activated receptors subtype.

The leishmanicidal flavonols quercetin and quercitrin target Leishmania (Leishmania) amazonensis arginase

Polyamine biosynthesis enzymes are promising drug targets for the treatment of leishmaniasis, Chagas' disease and African sleeping sickness. Arginase, which is a metallohydrolase, is the first enzyme involved in polyamine biosynthesis and converts arginine into ornithine and urea. Ornithine is used in the polyamine pathway that is essential for cell proliferation and ROS detoxification by trypanothione. The flavonols quercetin and quercitrin have been described as antitrypanosomal and antileishmanial compounds, and their ability to inhibit arginase was tested in this work. We characterized the inhibition of recombinant arginase from Leishmania (Leishmania) amazonensis by quercetin, quercitrin and isoquercitrin. The IC50 values for quercetin, quercitrin and isoquercitrin were estimated to be 3.8, 10 and 4.3 mu M, respectively. Quercetin is a mixed inhibitor, whereas quercitrin and isoquercitrin are uncompetitive inhibitors of L. (L.) amazonensis arginase. Quercetin interacts with the substrate L-arginine and the cofactor Mn2+ at pH 9.6, whereas quercitrin and isoquercitrin do not interact with the enzyme's cofactor or substrate. Docking analysis of these flavonols suggests that the cathecol group of the three compounds interact with Asp129, which is involved in metal bridge formation for the cofactors Mn-A(2+) and Mn-B(2+) in the active site of arginase. These results help to elucidate the mechanism of action of leishmanicidal flavonols and offer new perspectives for drug design against Leishmania infection based on interactions between arginase and flavones. (C) 2012 Elsevier Inc. All rights reserved.

Characterization of Heparin-induced Glyceraldehyde-3-phosphate Dehydrogenase Early Amyloid-like Oligomers and Their Implication in alpha-Synuclein Aggregation

Lewy bodies and Lewy neurites, neuropathological hallmarks of several neurological diseases, are mainly made of filamentous assemblies of alpha-synuclein. However, other macromolecules including Tau, ubiquitin, glyceraldehyde-3-phosphate dehydrogenase, and glycosaminoglycans are routinely found associated with these amyloid deposits. Glyceraldehyde-3-phosphate dehydrogenase is a glycolytic enzyme that can form fibrillar aggregates in the presence of acidic membranes, but its role in Parkinson disease is still unknown. In this work, the ability of heparin to trigger the amyloid aggregation of this protein at physiological conditions of pH and temperature is demonstrated by infrared and fluorescence spectroscopy, dynamic light scattering, small angle x-ray scattering, circular dichroism, and fluorescence microscopy. Aggregation proceeds through the formation of short rod-like oligomers, which elongates in one dimension. Heparan sulfate was also capable of inducing glyceraldehyde-3-phosphate dehydrogenase aggregation, but chondroitin sulfates A, B, and C together with dextran sulfate had a negligible effect. Aided with molecular docking simulations, a putative binding site on the protein is proposed providing a rational explanation for the structural specificity of heparin and heparan sulfate. Finally, it is demonstrated that in vitro the early oligomers present in the glyceraldehyde-3-phosphate dehydrogenase fibrillation pathway promote alpha-synuclein aggregation. Taking into account the toxicity of alpha-synuclein prefibrillar species, the heparin-induced glyceraldehyde-3-phosphate dehydrogenase early oligomers might come in useful as a novel therapeutic strategy in Parkinson disease and other synucleinopathies.

Looking over Toxin-K+ Channel Interactions. Clues from the Structural and Functional Characterization of alpha-KTx Toxin Tc32, a Kv1.3 Channel Blocker

alpha-KTx toxin Tc32, from the Amazonian scorpion Tityus cambridgei, lacks the dyad motif; including Lys27, characteristic of the family and generally associated with channel blockage. The toxin has been cloned and expressed for the first time. Electrophysiological experiments, by showing that the recombinant form blocks Kv1.3 channels of olfactory bulb periglomerular cells like the natural Tc32 toxin, when tested on the Kv1.3 channel of human T lymphocytes, confirmed it is in an active fold. The nuclear magnetic resonance-derived structure revealed it exhibits an alpha/beta scaffold typical of the members of the alpha-KTx family. TdK2 and TdK3, all belonging to the same alpha-KTx 18 subfamily, share significant sequence identity with Tc32 but diverse selectivity and affinity for Kv1.3 and Kv1.1 channels. To gain insight into the structural features that may justify those differences, we used the recombinant Tc32 nuclear magnetic resonance-derived structure to model the other two toxins, for which no experimental structure is available. Their interaction with Kv1.3 and Kv1.1 has been investigated by means of docking simulations. The results suggest that differences in the electrostatic features of the toxins and channels, in their contact surfaces, and in their total dipole moment orientations govern the affinity and selectivity of toxins. In addition, we found that, regardless of whether the dyad motif is present, it is always a Lys side chain that physically blocks the channels, irrespective of its position in the toxin sequence.