Structure-activity relationships of novel N-imidazoylpiperazines with potent anti-Trypanosoma cruzi activity.

Background: Chagas disease is caused by the parasite Trypanosoma cruzi, and the lack of effective and safe treatments makes identifying new classes of compounds with anti-T. cruzi activity of paramount importance. Methods: Hit-to-lead exploration of a metabolically stable N-imidazoylpiperazine was performed. Results: Compound 2, a piperazine derivative active against T. cruzi, was selected to perform the hit-to-lead exploration, which involved the design, synthesis and biological evaluation of 39 new derivatives. Conclusion: Compounds 6e and 10a were identified as optimized compounds with low micromolar in vitro activity, low cytotoxicity and suitable preliminary absorption, distribution, metabolism and excretion and physicochemical properties. Both compounds reduced parasitemia in mouse models of Chagas disease, providing a promising opportunity for further exploration of new antichagasic compounds.

Identification of a Prenyl Chalcone as a Competitive Lipoxygenase Inhibitor: Screening, Biochemical Evaluation and Molecular Modeling Studies.

Cyclooxygenase (COX) and lipoxygenase (LOX) are key targets for the development of new anti-inflammatory agents. LOX, which is involved in the biosynthesis of mediators in inflammation and allergic reactions, was selected for a biochemical screening campaign to identify LOX inhibitors by employing the main natural product library of Brazilian biodiversity. Two prenyl chalcones were identified as potent inhibitors of LOX-1 in the screening. The most active compound, (E)-2-O-farnesyl chalcone, decreased the rate of oxygen consumption to an extent similar to that of the positive control, nordihydroguaiaretic acid. Additionally, studies on the mechanism of the action indicated that (E)-2-O-farnesyl chalcone is a competitive LOX-1 inhibitor. Molecular modeling studies indicated the importance of the prenyl moieties for the binding of the inhibitors to the LOX binding site, which is related to their pharmacological properties.

Chalcones and their B-aryl analogues as myeloperoxidase inhibitors: In silico, in vitro and ex vivo investigations.

In the present study, a series of chalcones and their B-aryl analogues were prepared and evaluate as inhibitors of myeloperoxidase (MPO) chlorinating activity, using in vitro and ex vivo assays. Among these, B-thiophenyl chalcone (analogue 9) demonstrated inhibition of in vitro and ex vivo MPO chlorinating activity, exhibiting IC50 value of 0.53 and 19.2 µM, respectively. Potent ex vivo MPO inhibitors 5, 8 and 9 were not toxic to human neutrophils at 50 µM, as well as displayed weak 2,2-diphenyl-1-pycrylhydrazyl radical (DPPH•) and hypochlorous acid (HOCl) scavenger abilities. Docking simulations indicated binding mode of MPO inhibitors, evidencing hydrogen bonds between the amino group at 4'position (ring A) of chalcones with Gln91, Asp94, and Hys95 MPO residues. In this regard, the efficacy and low toxicity promoted aminochalcones and arylic analogues to the rank of hit compounds in the search for new non-steroidal anti-inflammatory compounds.

Multiparameter Optimization of Trypanocidal Cruzain Inhibitors With In Vivo Activity and Favorable Pharmacokinetics.

Cruzain, the main cysteine protease of Trypanosoma cruzi, plays key roles in all stages of the parasite's life cycle, including nutrition acquisition, differentiation, evasion of the host immune system, and invasion of host cells. Thus, inhibition of this validated target may lead to the development of novel drugs for the treatment of Chagas disease. In this study, a multiparameter optimization (MPO) approach, molecular modeling, and structure-activity relationships (SARs) were employed for the identification of new benzimidazole derivatives as potent competitive inhibitors of cruzain with trypanocidal activity and suitable pharmacokinetics. Extensive pharmacokinetic studies enabled the identification of metabolically stable and permeable compounds with high selectivity indices. CYP3A4 was found to be involved in the main metabolic pathway, and the identification of metabolic soft spots provided insights into molecular optimization. Compound 28, which showed a promising trade-off between pharmacodynamics and pharmacokinetics, caused no acute toxicity and reduced parasite burden both in vitro and in vivo.

Structure-Based and Molecular Modeling Studies for the Discovery of Cyclic Imides as Reversible Cruzain Inhibitors With Potent Anti-Trypanosoma cruzi Activity.

Chagas disease causes ~10,000 deaths each year, mainly in Latin America, where it is endemic. The currently available chemotherapeutic agents are ineffective in the chronic stage of the disease, and the lack of pharmaceutical innovation for Chagas disease highlights the urgent need for the development of new drugs. The enzyme cruzain, the main cysteine protease of Trypanosoma cruzi, has been explored as a validated molecular target for drug discovery. Herein, the design, molecular modeling studies, synthesis, and biological evaluation of cyclic imides as cruzain inhibitors are described. Starting with a micromolar-range cruzain inhibitor (3a, IC50 = 2.2 µM), this molecular optimization strategy resulted in the nanomolar-range inhibitor 10j (IC50 = 0.6 µM), which is highly active against T. cruzi intracellular amastigotes (IC50 = 1.0 µM). Moreover, most compounds were selective toward T. cruzi over human fibroblasts, which were used as host cells, and are less toxic to hepatic cells than the marketed drug benznidazole. This study enabled the discovery of novel chemical diversity and established robust structure-activity relationships to guide the design of optimized cruzain inhibitors as new trypanocidal agents.

Molecular modeling and structure-activity relationships for a series of benzimidazole derivatives as cruzain inhibitors.

AIM: Chagas disease is endemic in Latin America and no effective treatment is available. Efforts in drug research have focused on several enzymes from Trypanosoma cruzi, among which cruzain is a validated pharmacological target. METHODOLOGY: Chemometric analyses were performed on the data set using the hologram quantitative structure-activity relationship, comparative molecular field analysis and comparative molecular similarity index analysis methods. Docking simulations were executed using the crystallographic structure of cruzain in complex with a benzimidazole inhibitor. The top-scoring enzyme-inhibitor complexes were selected for the development of the 3D quantitative structure-activity relationship (QSAR) models and to assess the inhibitor binding modes and intermolecular interactions. RESULTS: Benzimidazole derivatives as cruzain inhibitors were used in molecular docking and QSAR studies. Significant statistical indicators were obtained, and the best models demonstrated high predictive ability for an external test set (r 2pred = 0.65, 0.94 and 0.82 for hologram QSAR, comparative molecular field analysis and comparative molecular similarity index analysis, respectively). Additionally, the graphical information of the chemometric analyses demonstrated substantial complementarity with the enzyme-binding site. CONCLUSION: These results demonstrate the relevance of the QSAR models to guide the design of structurally related benzimidazole derivatives with improved potency.

Functional, thermodynamics, structural and biological studies of in silico-identified inhibitors of Mycobacterium tuberculosis enoyl-ACP(CoA) reductase enzyme.

Novel chemotherapeutics agents are needed to kill Mycobacterium tuberculosis, the main causative agent of tuberculosis (TB). The M. tuberculosis 2-trans-enoyl-ACP(CoA) reductase enzyme (MtInhA) is the druggable bona fide target of isoniazid. New chemotypes were previously identified by two in silico approaches as potential ligands to MtInhA. The inhibition mode was determined by steady-state kinetics for seven compounds that inhibited MtInhA activity. Dissociation constant values at different temperatures were determined by protein fluorescence spectroscopy. van't Hoff analyses of ligand binding to MtInhA:NADH provided the thermodynamic signatures of non-covalent interactions (ΔH°, ΔS°, ΔG°). Phenotypic screening showed that five compounds inhibited in vitro growth of M. tuberculosis H37Rv strain. Labio_16 and Labio_17 compounds also inhibited the in vitro growth of PE-003 multidrug-resistant strain. Cytotoxic effects on Hacat, Vero and RAW 264.7 cell lines were assessed for the latter two compounds. The Labio_16 was bacteriostatic and Labio_17 bactericidal in an M. tuberculosis-infected macrophage model. In Zebrafish model, Labio_16 showed no cardiotoxicity whereas Labio_17 showed dose-dependent cardiotoxicity. Accordingly, a model was built for the MtInhA:NADH:Labio_16 ternary complex. The results show that the Labio_16 compound is a direct inhibitor of MtInhA, and it may represent a hit for the development of chemotherapeutic agents to treat TB.

Synthesis, biological evaluation, and structure-activity relationships of potent noncovalent and nonpeptidic cruzain inhibitors as anti-Trypanosoma cruzi agents.

The development of cruzain inhibitors has been driven by the urgent need to develop novel and more effective drugs for the treatment of Chagas' disease. Herein, we report the lead optimization of a class of noncovalent cruzain inhibitors, starting from an inhibitor previously cocrystallized with the enzyme (K(i) = 0.8 µM). With the goal of achieving a better understanding of the structure-activity relationships, we have synthesized and evaluated a series of over 40 analogues, leading to the development of a very promising competitive inhibitor (8r, IC50 = 200 nM, K(i) = 82 nM). Investigation of the in vitro trypanocidal activity and preliminary cytotoxicity revealed the potential of the most potent cruzain inhibitors in guiding further medicinal chemistry efforts to develop drug candidates for Chagas' disease.

Discovery of new inhibitors of Mycobacterium tuberculosis InhA enzyme using virtual screening and a 3D-pharmacophore-based approach.

Mycobacterium tuberculosis InhA (MtInhA) is an attractive enzyme to drug discovery efforts due to its validation as an effective biological target for tuberculosis therapy. In this work, two different virtual-ligand-screening approaches were applied in order to identify new InhA inhibitors' candidates from a library of ligands selected from the ZINC database. First, a 3-D pharmacophore model was built based on 36 available MtInhA crystal structures. By combining structure-based and ligand-based information, four pharmacophoric points were designed to select molecules able to satisfy the binding features of MtInhA substrate-binding cavity. The second approach consisted of using four well established docking programs, with different search algorithms, to compare the binding mode and score of the selected molecules from the aforementioned library. After detailed analyses of the results, six ligands were selected for in vitro analysis. Three of these molecules presented a satisfactory inhibitory activity with IC50 values ranging from 24 (±2) µM to 83 (±5) µM. The best compound presented an uncompetitive inhibition mode to NADH and 2-trans-dodecenoyl-CoA substrates, with Ki values of 24 (±3) µM and 20 (±2) µM, respectively. These molecules were not yet described as antituberculars or as InhA inhibitors, making its novelty interesting to start efforts on ligand optimization in order to identify new effective drugs against tuberculosis having InhA as a target. More studies are underway to dissect the discovered uncompetitive inhibitor interactions with MtInhA.

Conformational changes in 2-trans-enoyl-ACP (CoA) reductase (InhA) from M. tuberculosis induced by an inorganic complex: a molecular dynamics simulation study.

InhA, the NADH-dependent 2-trans-enoyl-ACP reductase enzyme from Mycobacterium tuberculosis (MTB), is involved in the biosynthesis of mycolic acids, the hallmark of mycobacterial cell wall. InhA has been shown to be the primary target of isoniazid (INH), one of the oldest synthetic antitubercular drugs. INH is a prodrug which is biologically activated by the MTB catalase-peroxidase KatG enzyme. The activation reaction promotes the formation of an isonicotinyl-NAD adduct which inhibits the InhA enzyme, resulting in reduction of mycolic acid biosynthesis. As a result of rational drug design efforts to design alternative drugs capable of inhibiting MTB's InhA, the inorganic complex pentacyano(isoniazid)ferrate(II) (PIF) was developed. PIF inhibited both wild-type and INH-resistant Ile21Val mutants of InhA and this inactivation did not require activation by KatG. Since no three-dimensional structure of the InhA-PIF complex is available to confirm the binding mode and to assess the molecular interactions with the protein active site residues, here we report the results of molecular dynamics simulations of PIF interaction with InhA. We found that PIF strongly interacts with InhA and that these interactions lead to macromolecular instabilities reflected in the long time necessary for simulation convergence. These instabilities were mainly due to perturbation of the substrate binding loop, particularly the partial denaturation of helices α6 and α7. We were also able to correlate the changes in the SASAs of Trp residues with the recent spectrofluorimetric investigation of the InhA-PIF complex and confirm their suggestion that the changes in fluorescence are due to InhA conformational changes upon PIF binding. The InhA-PIF association is very strong in the first 20.0 ns, but becomes very week at the end of the simulation, suggesting that the PIF binding mode we simulated may not reflect that of the actual InhA-PIF complex.

Crystal structure and molecular dynamics studies of human purine nucleoside phosphorylase complexed with 7-deazaguanine.

In humans, purine nucleoside phosphorylase (HsPNP) is responsible for degradation of deoxyguanosine, and genetic deficiency of this enzyme leads to profound T-cell mediated immunosuppression. HsPNP is a target for inhibitor development aiming at T-cell immune response modulation. Here we report the crystal structure of HsPNP in complex with 7-deazaguanine (HsPNP:7DG) at 2.75 A. Molecular dynamics simulations were employed to assess the structural features of HsPNP in both free form and in complex with 7DG. Our results show that some regions, responsible for entrance and exit of substrate, present a conformational variability, which is dissected by dynamics simulation analysis. Enzymatic assays were also carried out and revealed that 7-deazaguanine presents a lower inhibitory activity against HsPNP (K(i)=200 microM). The present structure may be employed in both structure-based design of PNP inhibitors and in development of specific empirical scoring functions.

Bioinformatics tools for screening of antiparasitic drugs.

Drug development has become the Holy Grail of many structural bioinformatics groups. The explosion of information about protein structures, ligand-binding affinity, parasite genome projects, and biological activity of millions of molecules opened the possibility to correlate this scattered information in order to generate reliable computational models to predict the likelihood of being able to modulate a target with a small-molecule drug. Computational methods have shown their potential in drug discovery and development allied with in vitro and in vivo methodologies. The present review discusses the main bioinformatics tools available for drug discovery and development.

Genomic databases and the search of protein targets for protozoan parasites.

The development of databases devoted to biological information opened the possibility to integrate, query and analyze biological data obtained from several sources that otherwise would be scattered through the web. Several issues arise in the handling of biological information, mainly due to the diversity of biological subject matter and the complexity of biological approaches towards phenomena of the living world. The integration of genomic data, three-dimensional structures of proteins, biological activity, and drugs availability allows a system approach to the study of the biology. Here we review the current status of these research efforts to develop genomic databases for protozoan parasites, such as the apicomplexan parasites, Trypanosoma cruzi and Leishmania spp. These databases may help in the discovery and development of new drugs against parasite-mediated diseases.

Protein-drug interaction studies for development of drugs against Plasmodium falciparum.

The study of protein-drug interaction is of pivotal importance to understand the structural features essential for ligand affinity. The explosion of information about protein structures has paved the way to develop structure-based virtual screening approaches. Parasitic protein kinases have been pointed out as potential targets for antiparasitic development. The identification of protein kinases in the Plasmodium falciparum genome has opened the possibility to test new families of inhibitors as potential antimalarial drugs. In addition, other key enzymes which play roles in biosynthetic pathways, such as enoyl reductase and chorismate synthase, can be valuable targets for drug development. This review is focused on these protein targets that may help to materialize new generations of antimalarial drugs.

Molecular modeling and dynamics studies of purine nucleoside phosphorylase from Bacteroides fragilis.

Purine nucleoside phosphorylase (PNP) catalyzes the reversible phosphorolysis of N-ribosidic bonds of purine nucleosides and deoxynucleosides, except adenosine, to generate ribose 1-phosphate and the purine base. This work describes for the first time a structural model of PNP from Bacteroides fragilis (Bf). We modeled the complexes of BfPNP with six different ligands in order to determine the structural basis for specificity of these ligands against BfPNP. Comparative analysis of the model of BfPNP and the structure of HsPNP allowed identification of structural features responsible for differences in the computationally determined ligand affinities. The molecular dynamics (MD) simulation was assessed to evaluate the overall stability of the BfPNP model. The superposition of the final onto the initial minimized structure shows that there are no major conformational changes from the initial model, which is consistent with the relatively low root mean square deviation (RMSD). The results indicate that the structure of the model was stable during MD, and does not exhibit loosely structured loop regions or domain terminals.

Molecular modeling and dynamics studies of Shikimate Kinase from Bacillus anthracis.

Bacillus anthracis has been used as weapon in bioterrorist activities, with high mortality, despite anti-microbial treatment, which strongly indicates a need of new drugs to treat anthrax. Shikimate Pathway is a seven-step biosynthetic route which generates chorismic acid. The shikimate pathway is essential for many pathological organisms, whereas it is absent in mammals. Therefore, these enzymes are potential targets for the development of non-toxic anti-microbial agents and herbicides and have been submitted to intensive structural studies. Shikimate Kinase is the fifth enzyme of shikimate pathway and catalyzes the specific phosphorylation of the 3-hydroxyl group of shikimate using ATP as a co-substrate, resulting in shikimate-3-phosphate and ADP. The present work describes for the first time a structural model for the Shikimate Kinase from B. anthracis using molecular modeling approach and molecular dynamics simulations. This study was able to identify the main residues of the ATP-binding and the shikimate pockets responsible for ligand affinities. Analysis of the molecular dynamics simulations indicates the structural features responsible for the stability of the structure. This study may help in the identification of new inhibitors for this enzyme.

In silico and in vitro: identifying new drugs.

Drug development is a high cost and laborious process, requiring a number of tests until a drug is made available in the market. Therefore, the use of methods to screen large number of molecules with less cost is crucial for faster identification of hits and leads. One strategy to identify drug-like molecules is the search for molecules able to interfere with a protein function, since protein interactions control most biological processes. Ideally the use of in silico screenings would make drug development faster and less expensive. Currently, however, the confirmation of biological activity is still needed. Due to the complexity of the task of drug discovery, an integrated and multi-disciplinary approach is ultimately required. Here we discuss examples of drugs developed through a combination of in silico and in vitro strategies. The potential use of these methodologies for the identification of active compounds as well as for early toxicity and bioavailability is also reviewed.

Molecular recognition models: a challenge to overcome.

Molecular recognition process describes the interaction involving two molecules. In the case of biomolecules, these pairs of molecules could be protein-protein, protein-ligand or protein-nucleic acid. The first model to capture the essential features, behind the molecular recognition problem, was the lock-and-key paradigm. The overall analysis protein-protein, protein-nucleic acid and protein-ligand interaction based on the three-dimensional structures and physicochemical parameters, such as binding affinity, opened the possibility to provide further insights in this basic phenomenon. The main ideas behind the molecular recognition are discussed in the present review.

Drug-binding databases.

Recent developments in computer power and chemoinformatics methodology make possible that a huge amount of data become available through internet. These databases are devoted to a wide spectrum of scientific fields. Here we are concerned with databases related to protein-drug interactions. More specifically, databases where potential new molecules could be accessed to be used in virtual screening initiatives. In the past decade several databases have been developed where molecules to be used in the virtual screening could be easily identified, downloaded and even purchased. This review describes and summarizes the recent advances in the development of these databases, and also the main applications related to virtual screening projects.

Molecular modeling as a tool for drug discovery.

With the progression of structural genomics projects, comparative modeling remains an increasingly important method of choice to obtain 3D structure of proteins. It helps to bridge the gap between the available sequence and structure information by providing reliable and accurate protein models. Comparative modeling based on more than 30% sequence identity is now approaching its natural template-based limits and further improvements require the development of effective refinement techniques capable of driving models toward native structure. For difficult targets, for which the most significant progress in recent years has been observed, optimal template selection and alignment accuracy are still the major problems. The past year has seen a maturation of molecular modeling, with an increasing number of comparative studies between established methods becoming possible, together with an explosion of new works especially in the areas of combinatorial chemistry and molecular diversity. To achieve this, knowledge about three-dimensional protein structures is crucial for the understanding of their functional mechanisms, and for a rational drug design. This review described recent progress in molecular modeling methodology.