3D-Pharmacophore mapping of thymidine-based inhibitors of TMPK as potential antituberculosis agents

Tuberculosis (TB) is the primary cause of mortality among infectious diseases. Mycobacterium tuberculosis monophosphate kinase (TMPKmt) is essential to DNA replication. Thus, this enzyme represents a promising target for developing new drugs against TB. In the present study, the receptor-independent, RI, 4D-QSAR method has been used to develop QSAR models and corresponding 3D-pharmacophores for a set of 81 thymidine analogues, and two corresponding subsets, reported as inhibitors of TMPKmt. The resulting optimized models are not only statistically significant with r (2) ranging from 0.83 to 0.92 and q (2) from 0.78 to 0.88, but also are robustly predictive based on test set predictions. The most and the least potent inhibitors in their respective postulated active conformations, derived from each of the models, were docked in the active site of the TMPKmt crystal structure. There is a solid consistency between the 3D-pharmacophore sites defined by the QSAR models and interactions with binding site residues. Moreover, the QSAR models provide insights regarding a probable mechanism of action of the analogues.

Rational Design and 3D-Pharmacophore Mapping of 5 `-Thiourea-Substituted alpha-Thymidine Analogues as Mycobacterial TMPK Inhibitors

Thymidine monophosphate kinase (TMPK) has emerged as an attractive target for developing inhibitors of Mycobacterium tuberculosis growth. In this study the receptor-independent (RI) 4D-QSAR formalism has been used to develop QSAR models and corresponding 3D-pharmacophores for a set of 5`-thiourea-substituted alpha-thymidine inhibitors. Models were developed for the entire training set and for a subset of the training set consisting of the most potent inhibitors. The optimized (RI) 4D-QSAR models are statistically significant (r(2) = 0.90, q(2) = 0.83 entire set, r(2) = 0.86, q(2) = 0.80 high potency subset) and also possess good predictivity based on test set predictions. The most and least potent inhibitors, in their respective postulated active conformations derived from the models, were docked in the active site of the TMPK crystallographic structure. There is a solid consistency between the 3D-pharmacophore sites defined by the QSAR models and interactions with binding site residues. This model identifies new regions of the inhibitors that contain pharmacophore sites, such as the sugar-pyrimidine ring structure and the region of the 5`-arylthiourea moiety. These new regions of the ligands can be further explored and possibly exploited to identify new, novel, and, perhaps, better antituberculosis inhibitors of TMPKmt. Furthermore, the 3D-pharmacophores defined by these models can be used as a starting point for future receptor-dependent antituberculosis drug design as well as to elucidate candidate sites for substituent addition to optimize ADMET properties of analog inhibitors.

Inhibition of immune complex-mediated neutrophil oxidative metabolism: A pharmacophore model for 3-phenylcoumarin derivatives using GRIND-based 3D-QSAR and 2D-QSAR procedures

In this study, twenty hydroxylated and acetoxylated 3-phenylcoumarin derivatives were evaluated as inhibitors of immune complex-stimulated neutrophil oxidative metabolism and possible modulators of the inflammatory tissue damage found in type III hypersensitivity reactions. By using lucigenin- and luminol-enhanced chemiluminescence assays (CL-luc and CL-lum, respectively), we found that the 6,7-dihydroxylated and 6,7-diacetoxylated 3-phenylcoumarin derivatives were the most effective inhibitors. Different structural features of the other compounds determined CL-luc and/or CL-lum inhibition. The 2D-QSAR analysis suggested the importance of hydrophobic contributions to explain these effects. In addition, a statistically significant 3D-QSAR model built applying GRIND descriptors allowed us to propose a virtual receptor site considering pharmacophoric regions and mutual distances. Furthermore, the 3-phenylcoumarins studied were not toxic to neutrophils under the assessed conditions. (C) 2007 Elsevier Masson SAS. All rights reserved.

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.

Discovery of New Inhibitors of Schistosoma mansoni PNP by Pharmacophore-Based Virtual Screening

Schistosomiasis is considered the second most important tropical parasitic disease, with severe socioeconomic consequences for millions of people worldwide. Schistosoma monsoni, one of the causative agents of human schistosomiasis, is unable to synthesize purine nucleotides de novo, which makes the enzymes of the purine salvage pathway important targets for antischistosomal drug development. In the present work, we describe the development of a pharmacophore model for ligands of S. mansoni purine nucleoside phosphorylase (SmPNP) as well as a pharmacophore-based virtual screening approach, which resulted in the identification of three thioxothiazolidinones (1-3) with substantial in vitro inhibitory activity against SmPNP. Synthesis, biochemical evaluation, and structure activity relationship investigations led to the successful development of a small set of thioxothiazolidinone derivatives harboring a novel chemical scaffold as new competitive inhibitors of SmPNP at the low-micromolar range. Seven compounds were identified with IC(50) values below 100 mu M. The most potent inhibitors 7, 10, and 17 with 1050 of 2, 18, and 38 mu M, respectively, could represent new potential lead compounds for further development of the therapy of schistosomiasis.

Pharmacophore-based 3D QSAR studies on a series of high affinity 5-HT(1A) receptor ligands

5-HT(1A) receptor antagonists have been employed to treat depression, but the lack of structural information on this receptor hampers the design of specific and selective ligands. In this study, we have performed CoMFA studies on a training set of arylpiperazines (high affinity 5-HT(1A) receptor ligands) and to produce an effective alignment of the data set, a pharmacophore model was produced using Galahad. A statistically significant model was obtained, indicating a good internal consistency and predictive ability for untested compounds. The information gathered from our receptor-independent pharmacophore hypothesis is in good agreement with results from independent studies using different approaches. Therefore, this work provides important insights on the chemical and structural basis involved in the molecular recognition of these compounds. (C) 2010 Elsevier Masson SAS. All rights reserved.

Pharmacophore Model of Cruzain Inhibitors

Chagas disease, caused by the protozoan Trypanosoma cruzi, is one of the most serious amongst the so-called neglected diseases in Latin America, specially in Brazil. So far there has been no effective treatment for the chronic phase of this disease. Cruzain is a major cysteine protease of T cruzi and it is recognized as a valid target for Chagas disease chemotherapy. The mechanism of cruzain action is associated with the nucleophilic attack of an activated sulfur atom towards electrophilic groups. In this report, features of a putative pharmacophore model of the enzyme, developed as a virtual screening tool for the selection of potential cruzain inhibitors, are described. The final proposed model was applied to the ZINC v.7 database and afterwards experimentally validated by an enzymatic inhibition assay. One of the compounds selected by the model showed cruzain inhibition in the low micromolar range.

The discovery of new 11beta-hydroxysteroid dehydrogenase type 1 inhibitors by common feature pharmacophore modeling and virtual screening

11beta-Hydroxysteroid dehydrogenase (11beta-HSD) enzymes catalyze the conversion of biologically inactive 11-ketosteroids into their active 11beta-hydroxy derivatives and vice versa. Inhibition of 11beta-HSD1 has considerable therapeutic potential for glucocorticoid-associated diseases including obesity, diabetes, wound healing, and muscle atrophy. Because inhibition of related enzymes such as 11beta-HSD2 and 17beta-HSDs causes sodium retention and hypertension or interferes with sex steroid hormone metabolism, respectively, highly selective 11beta-HSD1 inhibitors are required for successful therapy. Here, we employed the software package Catalyst to develop ligand-based multifeature pharmacophore models for 11beta-HSD1 inhibitors. Virtual screening experiments and subsequent in vitro evaluation of promising hits revealed several selective inhibitors. Efficient inhibition of recombinant human 11beta-HSD1 in intact transfected cells as well as endogenous enzyme in mouse 3T3-L1 adipocytes and C2C12 myotubes was demonstrated for compound 27, which was able to block subsequent cortisol-dependent activation of glucocorticoid receptors with only minor direct effects on the receptor itself. Our results suggest that inhibitor-based pharmacophore models for 11beta-HSD1 in combination with suitable cell-based activity assays, including such for related enzymes, can be used for the identification of selective and potent inhibitors.

Discovery of novel CB2 receptor ligands by a pharmacophore-based virtual screening workflow

Cannabinoid receptor 2 (CB(2) receptor) ligands are potential candidates for the therapy of chronic pain, inflammatory disorders, atherosclerosis, and osteoporosis. We describe the development of pharmacophore models for CB(2) receptor ligands, as well as a pharmacophore-based virtual screening workflow, which resulted in 14 hits for experimental follow-up. Seven compounds were identified with K(i) values below 25 microM. The CB(2) receptor-selective pyridine tetrahydrocannabinol analogue 8 (K(i) = 1.78 microM) was identified as a CB(2) partial agonist. Acetamides 12 (K(i) = 1.35 microM) and 18 (K(i) = 2.1 microM) represent new scaffolds for CB(2) receptor-selective antagonists and inverse agonists, respectively. Overall, our pharmacophore-based workflow yielded three novel scaffolds for the chemical development of CB(2) receptor ligands.

Library of libraries: approach to synthetic combinatorial library design and screening of "pharmacophore" motifs.

Construction of synthetic combinatorial libraries is described that allows for the generation of a library of motifs rather than a library of compounds. Peptide libraries based on this strategy were synthesized and screened with model targets streptavidin and anti-beta-endorphin antibody. The screens resulted in observation of expected motifs providing evidence of the effectiveness of the suggested approach.

Computer-aided Drug Design of Novel PLA(2) Inhibitor Candidates for Treatment of Snakebite

Phospholipases A(2) (PLA(2)) are enzymes commonly found in snake venoms from Viperidae and Elaphidae families, which are major components thereof. Many plants are used in traditional medicine its active agents against various effects induced by snakebite. This article presents the PLA(2) BthTX-I structure prediction based on homology modeling. In addition, we have performed virtual screening in a large database yielding a set of potential bioactive inhibitors. A flexible docking program was used to investigate the interactions between the receptor and the new ligands. We have performed molecular interaction fields (MIFs) calculations with the phospholipase model. Results confirm the important role of Lys49 for binding ligands and suggest three additional residues as well. We have proposed a theoretically nontoxic, drug-like, and potential novel BthTX-I inhibitor. These calculations have been used to guide the design of novel phospholipase inhibitors as potential lead compounds that may be optimized for future treatment of snakebite victims as well as other human diseases in which PLA(2) enzymes are involved.

Small molecular probes for G-protein-coupled C5a receptors: Conformationally constrained antagonists derived from the C terminus of the human plasma protein C5a

Activation of the human complement system of plasma proteins in response to infection or injury produces a 4-helix bundle glycoprotein (74 amino acids) known as C5a. C5a binds to G-protein-coupled receptors on cell surfaces triggering receptor-ligand internalization, signal transduction, and powerful inflammatory responses. Since excessive levels of C5a are associated with autoimmune and chronic inflammatory disorders, inhibitors of receptor activation may have therapeutic potential. We now report solution structures and receptor-binding and antagonist activities for some of the first small molecule antagonists of C5a derived from its hexapeptide C terminus. The antagonist NMe-Phe-Lys-Pro-D-Cha-Trp-D-Arg-CO2H (1) surprisingly shows an unusually well-defined solution structure as determined by H-1 NMR spectroscopy. This is one of the smallest acyclic peptides found to possess a defined solution conformation, which can be explained by the constraining role of intramolecular hydrogen bonding. NOE and coupling constant data, slow deuterium exchange, and a low dependence on temperature for the chemical shift of the D-Cha-NH strongly indicate an inverse gamma turn stabilized by a D-Cha-NH ... OC-Lys hydrogen bond. Smaller conformational populations are associated with a hydrogen bond between Trp-NH ... OC-Lys, defining a type II beta turn distorted by the inverse gamma turn incorporated within it. An excellent correlation between receptor-affinity and antagonist activity is indicated for a limited set of synthetic peptides. Conversion of the C-terminal carboxylate of 1 to an amide decreases antagonist potency 5-fold, but potency is increased up to 10-fold over 1 if the amide bond is made between the C-terminal carboxylate and a Lys/Orn side chain to form a cyclic analogue. The solution structure of cycle 6 also shows gamma and beta turns; however, the latter occurs in a different position, and there are clear conformational changes in 6 vs 1 that result in enhanced activity. These results indicate that potent C5a antagonists can be developed by targeting site 2 alone of the C5a receptor and define a novel pharmacophore for developing powerful receptor probes or drug candidates.

Structure-activity relationships of omega-conotoxins MVIIA, MVIIC and 14 loop splice hybrids at N and P/Q-type calcium channels

The omega-conotoxins are a set of structurally related, four-loop, six cysteine containing peptides, that have a range of selectivities for different subtypes of the voltage-sensitive calcium channel (VSCC). To investigate the basis of the selectivity displayed by these peptides, we have studied the binding affinities of two naturally occurring omega-conotoxins, MVIIA and MVIIC and a series of 14 MVIIA/MVIIC loop hybrids using radioligand binding assays for N and P/Q-type Ca2+ channels in rat brain tissue. A selectivity profile was developed from the ratio of relative potencies at N-type VSCCs (using [I-125]GVIA radioligand binding assays) and P/Q-type VSCCs (using [I-125]MVIIC radioligand binding assays). in these peptides, loops 2 and 4 make the greatest contribution to VSCC subtype selectivity, while the effects of loops 1 and 3 are negligible. Peptides with homogenous combinations of loop 2 and 4 display clear selectivity preferences, while those with heterogeneous combinations of loops 2 and 4 are less discriminatory. H-1 NMR spectroscopy revealed that the global folds of MVIIA, MVIIC and the 14 loop hybrid peptides were similar; however, several differences in local structure were identified. Based on the binding data and the 3D structures of MVIIA, GVIA and MVIIC, we have developed a preliminary pharmacophore based on the omega-conotoxin residues most Likely to interact with the N-type VSCC. (C) 1999 Academic Press.