Structure- and ligand-based drug design approaches for neglected tropical diseases

Drug discovery has moved toward more rational strategies based on our increasing understanding of the fundamental principles of protein-ligand interactions. Structure( SBDD) and ligand-based drug design (LBDD) approaches bring together the most powerful concepts in modern chemistry and biology, linking medicinal chemistry with structural biology. The definition and assessment of both chemical and biological space have revitalized the importance of exploring the intrinsic complementary nature of experimental and computational methods in drug design. Major challenges in this field include the identification of promising hits and the development of high-quality leads for further development into clinical candidates. It becomes particularly important in the case of neglected tropical diseases (NTDs) that affect disproportionately poor people living in rural and remote regions worldwide, and for which there is an insufficient number of new chemical entities being evaluated owing to the lack of innovation and R&D investment by the pharmaceutical industry. This perspective paper outlines the utility and applications of SBDD and LBDD approaches for the identification and design of new small-molecule agents for NTDs.

Extensive structural change of the envelope protein of dengue virus induced by a tuned ionic strength: conformational and energetic analyses

The Dengue has become a global public health threat, with over 100 million infections annually; to date there is no specific vaccine or any antiviral drug. The structures of the envelope (E) proteins of the four known serotype of the dengue virus (DENV) are already known, but there are insufficient molecular details of their structural behavior in solution in the distinct environmental conditions in which the DENVs are submitted, from the digestive tract of the mosquito up to its replication inside the host cell. Such detailed knowledge becomes important because of the multifunctional character of the E protein: it mediates the early events in cell entry, via receptor endocytosis and, as a class II protein, participates determinately in the process of membrane fusion. The proposed infection mechanism asserts that once in the endosome, at low pH, the E homodimers dissociate and insert into the endosomal lipid membrane, after an extensive conformational change, mainly on the relative arrangement of its three domains. In this work we employ all-atom explicit solvent Molecular Dynamics simulations to specify the thermodynamic conditions in that the E proteins are induced to experience extensive structural changes, such as during the process of reducing pH. We study the structural behavior of the E protein monomer at acid pH solution of distinct ionic strength. Extensive simulations are carried out with all the histidine residues in its full protonated form at four distinct ionic strengths. The results are analyzed in detail from structural and energetic perspectives, and the virtual protein movements are described by means of the principal component analyses. As the main result, we found that at acid pH and physiological ionic strength, the E protein suffers a major structural change; for lower or higher ionic strengths, the crystal structure is essentially maintained along of all extensive simulations. On the other hand, at basic pH, when all histidine residues are in the unprotonated form, the protein structure is very stable for ionic strengths ranging from 0 to 225 mM. Therefore, our findings support the hypothesis that the histidines constitute the hot points that induce configurational changes of E protein in acid pH, and give extra motivation to the development of new ideas for antivirus compound design.

Self-Organizing Maps of Molecular Descriptors for Sesquiterpene Lactones and Their Application to the Chemotaxonomy of the Asteraceae Family

The Asteraceae, one of the largest families among angiosperms, is chemically characterised by the production of sesquiterpene lactones (SLs). A total of 1,111 SLs, which were extracted from 658 species, 161 genera, 63 subtribes and 15 tribes of Asteraceae, were represented and registered in two dimensions in the SISTEMATX, an in-house software system, and were associated with their botanical sources. The respective 11 block of descriptors: Constitutional, Functional groups, BCUT, Atom-centred, 2D autocorrelations, Topological, Geometrical, RDF, 3D-MoRSE, GETAWAY and WHIM were used as input data to separate the botanical occurrences through self-organising maps. Maps that were generated with each descriptor divided the Asteraceae tribes, with total index values between 66.7% and 83.6%. The analysis of the results shows evident similarities among the Heliantheae, Helenieae and Eupatorieae tribes as well as between the Anthemideae and Inuleae tribes. Those observations are in agreement with systematic classifications that were proposed by Bremer, which use mainly morphological and molecular data, therefore chemical markers partially corroborate with these classifications. The results demonstrate that the atom-centred and RDF descriptors can be used as a tool for taxonomic classification in low hierarchical levels, such as tribes. Descriptors obtained through fragments or by the two-dimensional representation of the SL structures were sufficient to obtain significant results, and better results were not achieved by using descriptors derived from three-dimensional representations of SLs. Such models based on physico-chemical properties can project new design SLs, similar structures from literature or even unreported structures in two-dimensional chemical space. Therefore, the generated SOMs can predict the most probable tribe where a biologically active molecule can be found according Bremer classification.

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.

Design, synthesis and the effect of 1,2,3-triazole sialylmimetic neoglycoconjugates on Trypanosoma cruzi and its cell surface trans-sialidase

This work describes the synthesis of a series of sialylmimetic neoglycoconjugates represented by 1,4-disubstituted 1,2,3-triazole-sialic acid derivatives containing galactose modified at either C-1 or C-6 positions, glucose or gulose at C-3 position, and by the amino acid derivative 1,2,3-triazole fused threonine-3-O-galactose as potential TcTS inhibitors and anti-trypanosomal agents. This series was obtained by Cu(I)-catalysed azide-alkyne cycloaddition reaction ('click chemistry') between the azido-functionalized sugars 1-N(3)-Gal (commercial), 6-N(3)-Gal, 3-N(3)-Glc and 3-N(3)-Gul with the corresponding alkyne-based 2-propynyl-sialic acid, as well as by click chemistry reaction between the amino acid N(3)-ThrOBn with 3-O-propynyl-GalOMe. The 1,2,3-triazole linked sialic acid-6-O-galactose and the sialic acid-galactopyranoside showed high Trypanosoma cruzi trans-sialidase (TcTS) inhibitory activity at 1.0 mM (approx. 90%), whilst only the former displayed relevant trypanocidal activity (IC(50) 260 mu M). These results highlight the 1,2,3-triazole linked sialic acid-6-O-galactose as a prototype for further design of new neoglycoconjugates against Chagas' disease. (C) 2011 Elsevier Ltd. All rights reserved.

Molecular detection of enteropathogenic Escherichia coli in asymptomatic captive psittacines

Psittaciformes are one of the most endangered groups of birds, and several Brazilian species are classified between vulnerable and critically endangered. It is thus necessary to identify agents that cause infections in captive wild animals and to assess the risks posed thereof and to design interventions to minimize the possibility of disease outbreaks, leading to the conservation of endangered species. The purpose of this study was to identify enteropathogenic Escherichia coli (EPEC) cloacal isolates from asymptomatic psittacines in captivity and evaluate the distribution of the EPEC pathotype. Cloacal swabs were obtained from 46 asymptomatic birds, and resulting isolates were tested by polymerase chain reaction (PCR) for the presence of the attaching and effacing gene (eae) and bundle-forming pilus structural gene (bfpA) of EPEC. Samples from several species were tested, and three samples were found to be positive for the eae and bfpA genes and characterized as typical EPEC. This is the first report of this pathotype in asymptomatic psittacines. Although certain E. coli strains are more pathogenic than others, various factors should be considered when determining the potential of E. coli isolates to cause disease in captive psittacines. Birds that are positive for the EPEC (typical) strain could be zoonotic sources of infection, and may have acquired these strains through contact with humans or domestic animals. These findings may also be valuable for the long-term management of endangered species ex situ as one EPEC sample was isolated from a Red-tailed Amazon (Amazona brasiliensis).

Clinical, epidemiological and molecular features of the HIV-1 subtype C and recombinant forms that are circulating in the city of Sao Paulo, Brazil

Background: The city of Sao Paulo has the highest AIDS case rate, with nearly 60% in Brazil. Despite, several studies involving molecular epidemiology, lack of data regarding a large cohort study has not been published from this city. Objectives: This study aimed to describe the HIV-1 subtypes, recombinant forms and drug resistance mutations, according to subtype, with emphasis on subtype C and BC recombinants in the city of Sao Paulo, Brazil. Study design: RNA was extracted from the plasma samples of 302 HIV-1-seropositive subjects, of which 211 were drug-naive and 82 were exposed to ART. HIV-1 partial pol region sequences were used in phylogenetic analyses for subtyping and identification of drug resistance mutations. The envelope gene of subtype C and BC samples was also sequenced. Results: From partial pol gene analyses, 239 samples (79.1%) were assigned as subtype B, 23 (7.6%) were F1, 16 (5.3%) were subtype C and 24 (8%) were mosaics (3 CRF28/CRF29-like). The subtype C and BC recombinants were mainly identified in drug-naive patients (72.7%) and the heterosexual risk exposure category (86.3%), whereas for subtype B, these values were 69.9% and 57.3%, respectively (p = 0.97 and p = 0.015, respectively). An increasing trend of subtype C and BC recombinants was observed (p < 0.01). Conclusion: The HIV-1 subtype C and CRFs seem to have emerged over the last few years in the city of Sao Paulo, principally among the heterosexual population. These findings may have an impact on preventive measures and vaccine development in Brazil.

Genetic analysis of kernel oil content in tropical maize with design III and QTL mapping

Oil content and grain yield in maize are negatively correlated, and so far the development of high-oil high-yielding hybrids has not been accomplished. Then a fully understand of the inheritance of the kernel oil content is necessary to implement a breeding program to improve both traits simultaneously. Conventional and molecular marker analyses of the design III were carried out from a reference population developed from two tropical inbred lines divergent for kernel oil content. The results showed that additive variance was quite larger than the dominance variance, and the heritability coefficient was very high. Sixteen QTL were mapped, they were not evenly distributed along the chromosomes, and accounted for 30.91% of the genetic variance. The average level of dominance computed from both conventional and QTL analysis was partial dominance. The overall results indicated that the additive effects were more important than the dominance effects, the latter were not unidirectional and then heterosis could not be exploited in crosses. Most of the favorable alleles of the QTL were in the high-oil parental inbred, which could be transferred to other inbreds via marker-assisted backcross selection. Our results coupled with reported information indicated that the development of high-oil hybrids with acceptable yields could be accomplished by using marker-assisted selection involving oil content, grain yield and its components. Finally, to exploit the xenia effect to increase even more the oil content, these hybrids should be used in the Top Cross((TM)) procedure.

Sesquiterpene Lactones with Anti-Hepatitis C Virus Activity Using Molecular Descriptors

Hepatitis C is a worldwide public health problem. The available therapies are limited by their partial effectiveness and with meaningful side-effects. Sesquiterpene lactones (SLs) are a group of natural products with a wide variety of chemical structures and biological activities associated. There are few studies about the influence of the molecular structure of SLs for the anti-hepatitis C virus activity. In the present work, SLs are investigated in a subgenomic RNA replicon assay system and were analyzed using multiple linear regression along with self-organizing maps with DRAGON descriptors in order to identify the structural requirements for their biological activity and to predict the inhibitory potency of SLs. Characteristics such as stereochemistry and electronic effects demonstrated to be important for their anti-HCV activity, and the SOM produced a clear separation betwenn active and inactive compounds. Therefore, it is possible to use this map as a filter for virtual screening to predict the anti-HCV activity of SLs.

Molecular organization and doping in poly(2-methoxyaniline)/Ni(dmit)(2) films obtained with the Langmuir-Blodgett technique

The control of the properties of materials at the molecular level is pursued for many applications, especially those associated with nanostructures. In this paper, we show that the coordination compound [Ni(dmit)(2)], where (dmit) is the 1,3-dithiole-2-thione-4,5-dithiolate ligand, can induce doping of poly(2-methoxyaniline) (POMA) in molecularly ordered Langmuir and Langmuir-Blodgett (LB) films. Doping was associated with interactions between the components and the compression of the Langmuir film at the air-water interface, according to polarization-modulated infrared reflection-absorption spectroscopy (PM-IRRAS) data. Taking these results together with in situ UV-Vis absorption measurements, we could identify the molecular groups involved in the interaction, including the way they were reoriented upon film compression. The Langmuir films were sufficiently stable to be transferred as Y-type LB films, while the hybrid POMA/[Ni(dmit)(2)] films remain doped in the solid state. As expected, the molecular charges affected the film morphology, as observed from combined atomic and electric force microscopy measurements. In summary, with adequate spectroscopy and microscopy tools we characterized molecular-level interactions, which may allow one to design molecular electronic devices with controlled electrical properties.

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.

Synthesis, Biological Evaluation, And Molecular Modeling of Chalcone Derivatives As Potent Inhibitors of Mycobacterium tuberculosis Protein Tyrosine Phosphatases (PtpA and PtpB)

Tuberculosis (TB) is a major infectious disease caused by Mycobacterium tuberculosis (Mtb). According to the World Health Organization (WHO), about 1.8 million people die from TB and 10 million new cases are recorded each year. Recently, a new series of naphthylchalcones has been identified as inhibitors of Mtb protein tyrosine phosphatases (PTPs). In this work, 100 chalcones were designed, synthesized, and investigated for their inhibitory properties against MtbPtps. Structure-activity relationships (SAR) were developed, leading to the discovery of new potent inhibitors with IC50 values in the low-micromolar range. Kinetic studies revealed competitive inhibition and high selectivity toward the Mtb enzymes. Molecular modeling investigations were carried out with the aim of revealing the most relevant structural requirements underlying the binding affinity and selectivity of this series of inhibitors as potential anti-TB drugs.

Molecular mapping of the regenerative niche in a murine model of myocardial infarction

Adult stem cells are distributed through the whole organism, and present a great potential for the therapy of different types of disease. For the design of efficient therapeutic strategies, it is important to have a more detailed understanding of their basic biological characteristics, as well as of the signals produced by damaged tissues and to which they respond. Myocardial infarction (MI), a disease caused by a lack of blood flow supply in the heart, represents the most common cause of morbidity and mortality in the Western world. Stem cell therapy arises as a promising alternative to conventional treatments, which are often ineffective in preventing loss of cardiomyocytes and fibrosis. Cell therapy protocols must take into account the molecular events that occur in the regenerative niche of MI. In the present study, we investigated the expression profile of ten genes coding for chemokines or cytokines in a murine model of MI, aiming at the characterization of the regenerative niche. MI was induced in adult C57BL/6 mice and heart samples were collected after 24 h and 30 days, as well as from control animals, for quantitative RT-PCR. Expression of the chemokine genes CCL2, CCL3, CCL4, CCL7, CXCL2 and CXCL10 was significantly increased 24 h after infarction, returning to baseline levels on day 30. Expression of the CCL8 gene significantly increased only on day 30, whereas gene expression of CXCL12 and CX3CL1 were not significantly increased in either ischemic period. Finally, expression of the IL-6 gene increased 24 h after infarction and was maintained at a significantly higher level than control samples 30 days later. These results contribute to the better knowledge of the regenerative niche in MI, allowing a more efficient selection or genetic manipulation of cells in therapeutic protocols.

A Combined Study Using Ligand-Based Design, Synthesis, and Pharmacological Evaluation of Analogues of the Acetaminophen Ortho-Regioisomer with Potent Analgesic Activity

A ligand-based drug design study was performed to acetaminophen regioisomers as analgesic candidates employing quantum chemical calculations at the DFT/B3LYP level of theory and the 6-31G* basis set. To do so, many molecular descriptors were used such as highest occupied molecular orbital, ionization potential, HO bond dissociation energies, and spin densities, which might be related to quench reactivity of the tyrosyl radical to give N-acetyl-p-benzosemiquinone-imine through an initial electron withdrawing or hydrogen atom abstraction. Based on this in silico work, the most promising molecule, orthobenzamol, was synthesized and tested. The results expected from the theoretical prediction were confirmed in vivo using mouse models of nociception such as writhing, paw licking, and hot plate tests. All biological results suggested an antinociceptive activity mediated by opioid receptors. Furthermore, at 90 and 120 min, this new compound had an effect that was comparable to morphine, the standard drug for this test. Finally, the pharmacophore model is discussed according to the electronic properties derived from quantum chemistry calculations.

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.