Classical and hologram QSAR studies on a series of inhibitors of trypanosomatid Glyceraldehyde-3-Phosphate Dehydrogenase

Leishmaniasis and trypanosomiasis are major causes of morbidity and mortality in both tropical and subtropical regions of the world. The current available drugs are limited, ineffective, and require long treatment regimens. Due to the high dependence of trypanosomatids on glycolysis as a source of energy, some glycolytic enzymes have been identified as attractive targets for drug design. In the present work, classical Two-Dimensional Quantitative Structure -Activity Relationships (2D QSAR) and Hologram QSAR (HQSAR) studies were performed on a series of adenosine derivatives as inhibitors of Leishmania mexicana Glyceraldehyde-3-Phosphate Dehydrogenase (LmGAPDH). Significant correlation coefficients (classical QSAR, r(2)=0.83 and q(2) =0.81; HQSAR, r(2)=0.91 and q(2) =0.86) were obtained for the 56 training set compounds, indicating the potential of the models for untested compounds. The models were then externally validated using a test set of 14 structurally related compounds and the predicted values were in good agreement with the experimental results (classical QSAR, r(pred)(2) = 0.94; HQSAR, r(pred)(2) = 0.92).

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.

Natural products biological screening and ligand-based virtual screening for the discovery of new antileishmanial agents

In the course of our research program to discover novel antileishmanial agents, a biological screening of natural products against Leishmania major promastigotes allowed the identification of a furoquinoline alkaloid (1) and a furanocoumarin (2) as new hits. Subsequently, an integrated ligand-based virtual screening approach was employed to search for new antileishmanial compounds using these naturally occurring molecules as templates. Fourteen out of 40 compounds selected from a database of about 800,000 compounds (extracted from ZINC, a free database for virtual screening) were experimentally confirmed to possess significant in vitro antileishmanial properties. The application of ligand-based virtual screening as a complementary approach to experimental natural product screening was a useful strategy to facilitate the identification of new promising lead candidates.

Fragment-based and classical quantitative structure-activity relationships for a series of hydrazides as antituberculosis agents

Worldwide, tuberculosis (TB) is the leading cause of death among curable infectious diseases. Multidrug-resistant Mycobacterium tuberculosis is an emerging problem of great importance to public health, and there is an urgent need for new anti-TB drugs. In the present work, classical 2D quantitative structure-activity relationships (QSAR) and hologram QSAR (HQSAR) studies were performed on a training set of 91 isoniazid derivatives. Significant statistical models (classical QSAR, q(2) = 0.68 and r(2) = 0.72; HQSAR, q(2) = 0.63 and r(2) = 0.86) were obtained, indicating their consistency for untested compounds. The models were then used to evaluate an external test set containing 24 compounds which were not included in the training set, and the predicted values were in good agreement with the experimental results (HQSAR, r(pred)(2) = 0.87; classical QSAR, r(pred)(2) = 0.75).

Classical and fragment-based hologram structure-activity relationships for a series of analgesic cyclic imides

Cyclic imides have been widely employed in drug design research due to their multiple pharmacological and biological properties. In the present study, two-dimensional quantitative structure-activity relationship (2D QSAR) studies were conducted on a series of potent analgesic cyclic imides using both classical and hologram QSAR (HQSAR) methods, yielding significant statistical models (classical QSAR, q(2) = 0.80; HQSAR, q(2) = 0.84). The models were then used to evaluate an external data test, and the predicted values were in good agreement with the experimental results, indicating their consistency for untested compounds.

Structural basis for selective inhibition of purine nucleoside phosphorylase from Schistosoma mansoni: Kinetic and structural studies

Selectivity plays a crucial role in the design of enzyme inhibitors as novel antiparasitic agents, particularly in cases where the target enzyme is also present in the human host. Purine nucleoside phosphorylase from Schistosoma mansoni (SmPNP) is an attractive target for the discovery of potential antischistosomal agents. In the present work, kinetic studies were carried out in order to determine the inhibitory potency, mode of action and enzyme selectivity of a series of inhibitors of SmPNP. In addition, crystallographic studies provided important structural insights for rational inhibitor design, revealing consistent structural differences in the binding mode of the inhibitors in the active sites of the SmPNP and human PNP (HsPNP) structures. The molecular information gathered in this work should be useful for future medicinal chemistry efforts in the design of new inhibitors of SmPNP having increased affinity and selectivity. (C) 2010 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.

Volsurf Descriptors to Analyse Anti-HCV and Cytotoxic Activities of Sesquiterpene Lactones from Asteraceae Family

Various significant anti-HCV and cytotoxic sesquiterpene lactones (SLs) have been characterized. In this work, the chemometric tool Principal Component Analysis (PCA) was applied to two sets of SLs and the variance of the biological activity was explored. The first principal component accounts for as much of the variability in the data as possible, and each succeeding component accounts for as much of the remaining variability as possible. The calculations were performed using VolSurf program. For anti-HCV activity, PC1 (First Principal Component) explained 30.3% and PC2 (Second Principal Component) explained 26.5% of matrix total variance, while for cytotoxic activity, PC1 explained 30.9% and PC2 explained 15.6% of the total variance. The formalism employed generated good exploratory and predictive results and we identified some structural features, for both sets, important to the suitable biological activity and pharmacokinetic profile.

Structural studies of prephenate dehydratase from Mycobacterium tuberculosis H37Rv by SAXS, ultracentrifugation, and computational analysis

Tuberculosis (TB) is one of the most common infectious diseases known to man and responsible for millions of human deaths in the world. The increasing incidence of TB in developing countries, the proliferation of multidrug resistant strains, and the absence of resources for treatment have highlighted the need of developing new drugs against TB. The shikimate pathway leads to the biosynthesis of chorismate, a precursor of aromatic amino acids. This pathway is absent from mammals and shown to be essential for the survival of Mycobacterium tuberculosis, the causative agent of TB. Accordingly, enzymes of aromatic amino acid biosynthesis pathway represent promising targets for structure-based drug design. The first reaction in phenylalanine biosynthesis involves the conversion of chorismate to prephenate, catalyzed by chorismate mutase. The second reaction is catalyzed by prephenate dehydratase (PDT) and involves decarboxylation and dehydratation of prephenate to form phenylpyruvate, the precursor of phenylalanine. Here, we describe utilization of different techniques to infer the structure of M. tuberculosis PDT (MtbPDT) in solution. Small angle X-ray scattering and ultracentrifugation analysis showed that the protein oligomeric state is a tetramer and MtbPDT is a flat disk protein. Bioinformatics tools were used to infer the structure of MtbPDT A molecular model for MtbPDT is presented and molecular dynamics simulations indicate that MtbPDT i.s stable. Experimental and molecular modeling results were in agreement and provide evidence for a tetrameric state of MtbPDT in solution.

CDK9 a potential target for drug development

The family of Cyclin-Dependent Kinases (CDKs) can be subdivided into two major functional groups based on their roles in cell cycle and/or transcriptional control. CDK9 is the catalytic subunit of positive transcription elongation factor b (P-TEFb). CDK9 is the kinase of the TAK complex (Tat-associated kinase complex), and binds to Tat protein of HIV, suggesting a possible role for CDK9 in AIDS progression. CDK9 complexed with its regulatory partner cyclin T1, serves as a cellular mediator of the transactivation function of the HIV Tat protein. P-TEFb is responsible for the phosphorylation of the carboxyl-terminal domain of RNA Pol II, resulting in stimulation of transcription. Furthermore, the complexes containing CDK9 induce the differentiation in distinct tissue. The CDK9/cyclin T1 complex is expressed at higher level in more differentiated primary neuroectodermal and neuroblastoma tumors, showing a correlation between the kinase expression and tumor differentiation grade. This may have clinical and therapeutical implications for these tumor types. Among the CDK inhibitors two have shown to be effective against CDK9: Roscovitine and Flavopiridol. These two inhibitors prevented the replication of human immunodeficiency virus (HIV) type 1 by blocking Tat transactivation of the HIV type 1 promoter. These compounds inhibit CDKs by binding to the catalytic domain in place of ATP, preventing transfer of a phosphate group to the substrate. More sensitive therapeutic agents of CDK9 can be designed, and structural studies can add information in the understanding of this kinase. The major features related to CDK9 inhibition will be reviewed in this article.

Estrogen receptor: Structural differences and potential implications on selectivity examined by the GRID/CPCA approach

Selective Estrogen Receptor Modulators ( SERMs) have been developed, but the selectivity towards the subtypes ( ER or ER is not well understood. Based on three-dimensional structural properties of ligand binding domains, a model that takes into account this aspect was developed via molecular interaction fields and consensus principal component analysis (GRID/CPCA).

2D QSAR and similarity studies on cruzain inhibitors aimed at improving selectivity over cathepsin L

Hologram quantitative structure-activity relationships (HQSAR) were applied to a data set of 41 cruzain inhibitors. The best HQSAR model (Q(2) = 0.77; R-2 = 0.90) employing Surflex-Sim, as training and test sets generator, was obtained using atoms, bonds, and connections as fragment distinctions and 4-7 as fragment size. This model was then used to predict the potencies of 12 test set compounds, giving satisfactory predictive R-2 value of 0,88. The contribution maps obtained from the best HQSAR model are in agreement with the biological activities of the study compounds. The Trypanosoma cruzi cruzain shares high similarity with the mammalian homolog cathepsin L. The selectivity toward cruzam was checked by a database of 123 compounds, which corresponds to the 41 cruzain inhibitors used in the HQSAR model development plus 82 cathepsin L inhibitors. We screened these compounds by ROCS (Rapid Overlay of Chemical Structures), a Gaussian-shape volume overlap filter that can rapidly identify shapes that match the query molecule. Remarkably, ROCS was able to rank the first 37 hits as being only cruzain inhibitors. In addition, the area under the curve (AUC) obtained with ROCS was 0.96, indicating that the method was very efficient to distinguishing between cruzain and cathepsin L inhibitors. (c) 2007 Elsevier Ltd. All rights reserved.

Artificial Neural Networks and the Study of the Psychoactivity of Cannabinoid Compounds

Cannabinoid compounds have widely been employed because of its medicinal and psychotropic properties. These compounds are isolated from Cannabis sativa (or marijuana) and are used in several medical treatments, such as glaucoma, nausea associated to chemotherapy, pain and many other situations. More recently, its use as appetite stimulant has been indicated in patients with cachexia or AIDS. In this work, the influence of several molecular descriptors on the psychoactivity of 50 cannabinoid compounds is analyzed aiming one obtain a model able to predict the psychoactivity of new cannabinoids. For this purpose, initially, the selection of descriptors was carried out using the Fisher`s weight, the correlation matrix among the calculated variables and principal component analysis. From these analyses, the following descriptors have been considered more relevant: E(LUMO) (energy of the lowest unoccupied molecular orbital), Log P (logarithm of the partition coefficient), VC4 (volume of the substituent at the C4 position) and LP1 (Lovasz-Pelikan index, a molecular branching index). To follow, two neural network models were used to construct a more adequate model for classifying new cannabinoid compounds. The first model employed was multi-layer perceptrons, with algorithm back-propagation, and the second model used was the Kohonen network. The results obtained from both networks were compared and showed that both techniques presented a high percentage of correctness to discriminate psychoactive and psychoinactive compounds. However, the Kohonen network was superior to multi-layer perceptrons.

The Relationship between Population Structure and Aluminum Tolerance in Cultivated Sorghum

Background: Acid soils comprise up to 50% of the world's arable lands and in these areas aluminum (Al) toxicity impairs root growth, strongly limiting crop yield. Food security is thereby compromised in many developing countries located in tropical and subtropical regions worldwide. In sorghum, SbMATE, an Al-activated citrate transporter, underlies the Alt(SB) locus on chromosome 3 and confers Al tolerance via Al-activated root citrate release. Methodology: Population structure was studied in 254 sorghum accessions representative of the diversity present in cultivated sorghums. Al tolerance was assessed as the degree of root growth inhibition in nutrient solution containing Al. A genetic analysis based on markers flanking Alt(SB) and SbMATE expression was undertaken to assess a possible role for Alt(SB) in Al tolerant accessions. In addition, the mode of gene action was estimated concerning the Al tolerance trait. Comparisons between models that include population structure were applied to assess the importance of each subpopulation to Al tolerance. Conclusion/Significance: Six subpopulations were revealed featuring specific racial and geographic origins. Al tolerance was found to be rather rare and present primarily in guinea and to lesser extent in caudatum subpopulations. Alt(SB) was found to play a role in Al tolerance in most of the Al tolerant accessions. A striking variation was observed in the mode of gene action for the Al tolerance trait, which ranged from almost complete recessivity to near complete dominance, with a higher frequency of partially recessive sources of Al tolerance. A possible interpretation of our results concerning the origin and evolution of Al tolerance in cultivated sorghum is discussed. This study demonstrates the importance of deeply exploring the crop diversity reservoir both for a comprehensive view of the dynamics underlying the distribution and function of Al tolerance genes and to design efficient molecular breeding strategies aimed at enhancing Al tolerance.

Chemometric Studies on Natural Products as Potential Inhibitors of the NADH Oxidase from Trypanosoma cruzi Using the VolSurf Approach

Natural products have widespread biological activities, including inhibition of mitochondrial enzyme systems. Some of these activities, for example cytotoxicity, may be the result of alteration of cellular bioenergetics. Based on previous computer-aided drug design (CADD) studies and considering reported data on structure-activity relationships (SAR), an assumption regarding the mechanism of action of natural products against parasitic infections involves the NADH-oxidase inhibition. In this study, chemometric tools, such as: Principal Component Analysis (PCA), Consensus PCA (CPCA), and partial least squares regression (PLS), were applied to a set of forty natural compounds, acting as NADH-oxidase inhibitors. The calculations were performed using the VolSurf+ program. The formalisms employed generated good exploratory and predictive results. The independent variables or descriptors having a hydrophobic profile were strongly correlated to the biological data.