Triclabendazole and clofazimine reduce replication and spermine uptake in vitro in Toxoplasma gondii.

Toxoplasmosis is a worldwide zoonosis caused by the protozoan parasite Toxoplasma gondii. Although this infection is generally asymptomatic in immunocompetent individuals, it can cause serious clinical manifestations in newborns with congenital infection or in immunocompromised patients. As current treatments are not always well tolerated, there is an urgent need to find new drugs against human toxoplasmosis. Drug repurposing has gained considerable momentum in the last decade and is a particularly attractive approach for the search of therapeutic alternatives to treat rare and neglected diseases. Thus, in this study, we investigated the antiproliferative effect of several repurposed drugs. Of these, clofazimine and triclabendazole displayed a higher selectivity against T. gondii, affecting its replication. Furthermore, both compounds inhibited spermine incorporation into the parasite, which is necessary for the formation of other polyamines. The data reported here indicate that clofazimine and triclabendazole could be used for the treatment of human toxoplasmosis and confirms that drug repurposing is an excellent strategy to find new therapeutic targets of intervention.

Garbage in, garbage out: how reliable training data improved a virtual screening approach against SARS-CoV-2 MPro.

Introduction: The identification of chemical compounds that interfere with SARS-CoV-2 replication continues to be a priority in several academic and pharmaceutical laboratories. Computational tools and approaches have the power to integrate, process and analyze multiple data in a short time. However, these initiatives may yield unrealistic results if the applied models are not inferred from reliable data and the resulting predictions are not confirmed by experimental evidence. Methods: We undertook a drug discovery campaign against the essential major protease (MPro) from SARS-CoV-2, which relied on an in silico search strategy -performed in a large and diverse chemolibrary- complemented by experimental validation. The computational method comprises a recently reported ligand-based approach developed upon refinement/learning cycles, and structure-based approximations. Search models were applied to both retrospective (in silico) and prospective (experimentally confirmed) screening. Results: The first generation of ligand-based models were fed by data, which to a great extent, had not been published in peer-reviewed articles. The first screening campaign performed with 188 compounds (46 in silico hits and 100 analogues, and 40 unrelated compounds: flavonols and pyrazoles) yielded three hits against MPro (IC50 ≤ 25 µM): two analogues of in silico hits (one glycoside and one benzo-thiazol) and one flavonol. A second generation of ligand-based models was developed based on this negative information and newly published peer-reviewed data for MPro inhibitors. This led to 43 new hit candidates belonging to different chemical families. From 45 compounds (28 in silico hits and 17 related analogues) tested in the second screening campaign, eight inhibited MPro with IC50 = 0.12-20 µM and five of them also impaired the proliferation of SARS-CoV-2 in Vero cells (EC50 7-45 µM). Discussion: Our study provides an example of a virtuous loop between computational and experimental approaches applied to target-focused drug discovery against a major and global pathogen, reaffirming the well-known "garbage in, garbage out" machine learning principle.

Making the most effective use of available computational methods for drug repositioning.

INTRODUCTION: Over the last decades, there has been substantial debate around the apparent drop in productivity in the pharmaceutical sector. The development of second or further medical uses for known drugs is a possible answer to expedite the development of new therapeutic solutions. Computational methods are among the main strategies for exploring drug repurposing opportunities in a systematic manner. AREAS COVERED: This article reviews three general approximations to systematically discover new therapeutic uses for existing drugs: disease-, target-, and drug-centric approaches, along with some recently reported computational methods associated with them. EXPERT OPINION: Computational methods are essential for organizing and analyzing the large volume of available biomedical data, which has grown exponentially in the era of big data. The clearest trend in the field involves the use of integrative approaches where different types of data are combined into multipartite networks. Every aspect of computer-guided drug repositioning has currently incorporated state-of-the-art machine learning tools to boost their pattern recognition and predictive capabilities. Remarkably, a majority of the recently reported platforms are publicly available as web apps or open-source software. The introduction of nationwide electronic health records provides invaluable real-world data to detect unknown relationships between approved drug treatments and diseases.

A combined ligand and target-based virtual screening strategy to repurpose drugs as putrescine uptake inhibitors with trypanocidal activity.

Chagas disease, also known as American trypanosomiasis, is a neglected tropical disease caused by the protozoa Trypanosoma cruzi, affecting nearly 7 million people only in the Americas. Polyamines are essential compounds for parasite growth, survival, and differentiation. However, because trypanosomatids are auxotrophic for polyamines, they must be obtained from the host by specific transporters. In this investigation, an ensemble of QSAR classifiers able to identify polyamine analogs with trypanocidal activity was developed. Then, a multi-template homology model of the dimeric polyamine transporter of T. cruzi, TcPAT12, was created with Rosetta, and then refined by enhanced sampling molecular dynamics simulations. Using representative snapshots extracted from the trajectory, a docking model able to discriminate between active and inactive compounds was developed and validated. Both models were applied in a parallel virtual screening campaign to repurpose known drugs as anti-trypanosomal compounds inhibiting polyamine transport in T. cruzi. Montelukast, Quinestrol, Danazol, and Dutasteride were selected for in vitro testing, and all of them inhibited putrescine uptake in biochemical assays, confirming the predictive ability of the computational models. Furthermore, all the confirmed hits proved to inhibit epimastigote proliferation, and Quinestrol and Danazol were able to inhibit, in the low micromolar range, the viability of trypomastigotes and the intracellular growth of amastigotes.

iRaPCA and SOMoC: Development and Validation of Web Applications for New Approaches for the Clustering of Small Molecules.

The clustering of small molecules implies the organization of a group of chemical structures into smaller subgroups with similar features. Clustering has important applications to sample chemical datasets or libraries in a representative manner (e.g., to choose, from a virtual screening hit list, a chemically diverse subset of compounds to be submitted to experimental confirmation, or to split datasets into representative training and validation sets when implementing machine learning models). Most strategies for clustering molecules are based on molecular fingerprints and hierarchical clustering algorithms. Here, two open-source in-house methodologies for clustering of small molecules are presented: iterative Random subspace Principal Component Analysis clustering (iRaPCA), an iterative approach based on feature bagging, dimensionality reduction, and K-means optimization; and Silhouette Optimized Molecular Clustering (SOMoC), which combines molecular fingerprints with the Uniform Manifold Approximation and Projection (UMAP) and Gaussian Mixture Model algorithm (GMM). In a benchmarking exercise, the performance of both clustering methods has been examined across 29 datasets containing between 100 and 5000 small molecules, comparing these results with those given by two other well-known clustering methods, Ward and Butina. iRaPCA and SOMoC consistently showed the best performance across these 29 datasets, both in terms of within-cluster and between-cluster distances. Both iRaPCA and SOMoC have been implemented as free Web Apps and standalone applications, to allow their use to a wide audience within the scientific community.

Tetracycline Derivatives Inhibit Plasmodial Cysteine Protease Falcipain-2 through Binding to a Distal Allosteric Site.

Allosteric inhibitors regulate enzyme activity from remote and usually specific pockets. As they promise an avenue for less toxic and safer drugs, the identification and characterization of allosteric inhibitors has gained great academic and biomedical interest in recent years. Research on falcipain-2 (FP-2), the major papain-like cysteine hemoglobinase of Plasmodium falciparum, might benefit from this strategy to overcome the low selectivity against human cathepsins shown by active site-directed inhibitors. Encouraged by our previous finding that methacycline inhibits FP-2 noncompetitively, here we assessed other five tetracycline derivatives against this target and characterized their inhibition mechanism. As previously shown for methacycline, tetracycline derivatives inhibited FP-2 in a noncompetitive fashion, with Ki values ranging from 121 to 190 µM. A possible binding to the S' side of the FP-2 active site, similar to that described by X-ray crystallography (PDB: 6SSZ) for the noncompetitive inhibitor E-chalcone 48 (EC48), was experimentally discarded by kinetic analysis using a large peptidyl substrate spanning the whole active site. By combining lengthy molecular dynamics (MD) simulations that allowed methacycline to diffuse from solution to different FP-2 surface regions and free energy calculations, we predicted the most likely binding mode of the ligand. Of note, the proposed binding pose explains the low differences in Ki values observed for the tested tetracycline derivatives and the calculated binding free energies match the experimental values. Overall, this study has implications for the design of novel allosteric inhibitors against FP-2 and sets the basis for further optimization of the tetracycline scaffold to produce more potent and selective inhibitors.

Homology Modeling and Molecular Dynamics Simulations of Trypanosoma cruzi Phosphodiesterase b1.

Cyclic nucleotide phosphodiesterases have been implicated in the proliferation, differentiation and osmotic regulation of trypanosomatids; in some trypanosomatid species, they have been validated as molecular targets for the development of new therapeutic agents. Because the experimental structure of Trypanosoma cruzi PDEb1 (TcrPDEb1) has not been solved so far, an homology model of the target was created using the structure of Trypanosoma brucei PDEb1 (TbrPDEb1) as a template. The model was refined by extensive enhanced sampling molecular dynamics simulations, and representative snapshots were extracted from the trajectory by combined clustering analysis. This structural ensemble was used to develop a structure-based docking model of the target. The docking accuracy of the model was validated by redocking and cross-docking experiments using all available crystal structures of TbrPDEb1, whereas the scoring accuracy was validated through a retrospective screen, using a carefully curated dataset of compounds assayed against TbrPDEb1 and/or TcrPDEb1. Considering the results from in silico validations, the model may be applied in prospective virtual screening campaigns to identify novel hits, as well as to guide the rational design of potent and selective inhibitors targeting this enzyme.

Strengths and Weaknesses of Docking Simulations in the SARS-CoV-2 Era: the Main Protease (Mpro) Case Study.

The scientific community is working against the clock to arrive at therapeutic interventions to treat patients with COVID-19. Among the strategies for drug discovery, virtual screening approaches have the capacity to search potential hits within millions of chemical structures in days, with the appropriate computing infrastructure. In this article, we first analyzed the published research targeting the inhibition of the main protease (Mpro), one of the most studied targets of SARS-CoV-2, by docking-based methods. An alarming finding was the lack of an adequate validation of the docking protocols (i.e., pose prediction and virtual screening accuracy) before applying them in virtual screening campaigns. The performance of the docking protocols was tested at some level in 57.7% of the 168 investigations analyzed. However, we found only three examples of a complete retrospective analysis of the scoring functions to quantify the virtual screening accuracy of the methods. Moreover, only two publications reported some experimental evaluation of the proposed hits until preparing this manuscript. All of these findings led us to carry out a retrospective performance validation of three different docking protocols, through the analysis of their pose prediction and screening accuracy. Surprisingly, we found that even though all tested docking protocols have a good pose prediction, their screening accuracy is quite limited as they fail to correctly rank a test set of compounds. These results highlight the importance of conducting an adequate validation of the docking protocols before carrying out virtual screening campaigns, and to experimentally confirm the predictions made by the models before drawing bold conclusions. Finally, successful structure-based drug discovery investigations published during the redaction of this manuscript allow us to propose the inclusion of target flexibility and consensus scoring as alternatives to improve the accuracy of the methods.

Application of target repositioning and in silico screening to exploit fatty acid binding proteins (FABPs) from Echinococcus multilocularis as possible drug targets.

Fatty acid binding proteins (FABPs) are small intracellular proteins that reversibly bind fatty acids and other hydrophobic ligands. In cestodes, due to their inability to synthesise fatty acids and cholesterol de novo, FABPs, together with other lipid binding proteins, have been proposed as essential, involved in the trafficking and delivery of such lipophilic metabolites. Pharmacological agents that modify specific parasite FABP function may provide control of lipid signalling pathways, inflammatory responses and metabolic regulation that could be of crucial importance for the parasite development and survival. Echinococcus multilocularis and Echinococcus granulosus are, respectively, the causative agents of alveolar and cystic echinococcosis (or hydatidosis). These diseases are included in the World Health Organization's list of priority neglected tropical diseases. Here, we explore the potential of FABPs from cestodes as drug targets. To this end, we have applied a target repurposing approach to identify novel inhibitors of Echinococcus spp. FABPs. An ensemble of computational models was developed and applied in a virtual screening campaign of DrugBank library. 21 hits belonging to the applicability domain of the ensemble models were identified, and 3 of the hits were assayed against purified E. multilocularis FABP, experimentally confirming the model's predictions. Noteworthy, this is to our best knowledge the first report on isolation and purification of such four FABP, for which initial structural and functional characterization is reported here.

Positive Predictive Value Surfaces as a Complementary Tool to Assess the Performance of Virtual Screening Methods.

BACKGROUND: Since their introduction in the virtual screening field, Receiver Operating Characteristic (ROC) curve-derived metrics have been widely used for benchmarking of computational methods and algorithms intended for virtual screening applications. Whereas in classification problems, the ratio between sensitivity and specificity for a given score value is very informative, a practical concern in virtual screening campaigns is to predict the actual probability that a predicted hit will prove truly active when submitted to experimental testing (in other words, the Positive Predictive Value - PPV). Estimation of such probability is however, obstructed due to its dependency on the yield of actives of the screened library, which cannot be known a priori. OBJECTIVE: To explore the use of PPV surfaces derived from simulated ranking experiments (retrospective virtual screening) as a complementary tool to ROC curves, for both benchmarking and optimization of score cutoff values. METHODS: The utility of the proposed approach is assessed in retrospective virtual screening experiments with four datasets used to infer QSAR classifiers: inhibitors of Trypanosoma cruzi trypanothione synthetase; inhibitors of Trypanosoma brucei N-myristoyltransferase; inhibitors of GABA transaminase and anticonvulsant activity in the 6 Hz seizure model. RESULTS: Besides illustrating the utility of PPV surfaces to compare the performance of machine learning models for virtual screening applications and to select an adequate score threshold, our results also suggest that ensemble learning provides models with better predictivity and more robust behavior. CONCLUSION: PPV surfaces are valuable tools to assess virtual screening tools and choose score thresholds to be applied in prospective in silico screens. Ensemble learning approaches seem to consistently lead to improved predictivity and robustness.

In Silico Drug Repositioning for Chagas Disease.

Chagas disease is an infectious tropical disease included within the group of neglected tropical diseases. Though historically endemic to Latin America, it has lately spread to high-income countries due to human migration. At present, there are only two available drugs, nifurtimox and benznidazole, approved for this treatment, both with considerable side-effects (which often result in treatment interruption) and limited efficacy in the chronic stage of the disease in adults. Drug repositioning involves finding novel therapeutic indications for known drugs, including approved, withdrawn, abandoned and investigational drugs. It is today a broadly applied approach to develop innovative medications, since indication shifts are built on existing safety, ADME and manufacturing information, thus greatly shortening development timeframes. Drug repositioning has been signaled as a particularly interesting strategy to search for new therapeutic solutions for neglected and rare conditions, which traditionally present limited commercial interest and are mostly covered by the public sector and not-for-profit initiatives and organizations. Here, we review the applications of computer-aided technologies as systematic approaches to drug repositioning in the field of Chagas disease. In silico screening represents the most explored approach, whereas other rational methods such as network-based and signature-based approximations have still not been applied.

In silico Guided Drug Repurposing: Discovery of New Competitive and Non-competitive Inhibitors of Falcipain-2.

Malaria is among the leading causes of death worldwide. The emergence of Plasmodium falciparum resistant strains with reduced sensitivity to the first line combination therapy and suboptimal responses to insecticides used for Anopheles vector management have led to renewed interest in novel therapeutic options. Here, we report the development and validation of an ensemble of ligand-based computational models capable of identifying falcipain-2 inhibitors, and their subsequent application in the virtual screening of DrugBank and Sweetlead libraries. Among four hits submitted to enzymatic assays, two (odanacatib, an abandoned investigational treatment for osteoporosis and bone metastasis, and the antibiotic methacycline) confirmed inhibitory effects on falcipain-2, with Ki of 98.2 nM and 84.4 µM. Interestingly, Methacycline proved to be a non-competitive inhibitor (α = 1.42) of falcipain-2. The effects of both hits on falcipain-2 hemoglobinase activity and on the development of P. falciparum were also studied.

Molecular Topology and Other Promiscuity Determinants as Predictors of Therapeutic Class - A Theoretical Framework to Guide Drug Repositioning?

Much interest has been paid in the last decade on molecular predictors of promiscuity, including molecular weight, log P, molecular complexity, acidity constant and molecular topology, with correlations between promiscuity and those descriptors seemingly being context-dependent. It has been observed that certain therapeutic categories (e.g. mood disorders therapies) display a tendency to include multi-target agents (i.e. selective non-selectivity). Numerous QSAR models based on topological descriptors suggest that the topology of a given drug could be used to infer its therapeutic applications. Here, we have used descriptive statistics to explore the distribution of molecular topology descriptors and other promiscuity predictors across different therapeutic categories. Working with the publicly available ChEMBL database and 14 molecular descriptors, both hierarchical and non-hierchical clustering methods were applied to the descriptors mean values of the therapeutic categories after the refinement of the database (770 drugs grouped into 34 therapeutic categories). On the other hand, another publicly available database (repoDB) was used to retrieve cases of clinically-approved drug repositioning examples that could be classified into the therapeutic categories considered by the aforementioned clusters (111 cases), and the correspondence between the two studies was evaluated. Interestingly, a 3- cluster hierarchical clustering scheme based on only 14 molecular descriptors linked to promiscuity seem to explain up to 82.9% of approved cases of drug repurposing retrieved of repoDB. Therapeutic categories seem to display distinctive molecular patterns, which could be used as a basis for drug screening and drug design campaigns, and to unveil drug repurposing opportunities between particular therapeutic categories.

Cascade Ligand- and Structure-Based Virtual Screening to Identify New Trypanocidal Compounds Inhibiting Putrescine Uptake.

Chagas disease is a neglected tropical disease endemic to Latin America, though migratory movements have recently spread it to other regions. Here, we have applied a cascade virtual screening campaign combining ligand- and structure-based methods. In order to find novel inhibitors of putrescine uptake in Trypanosoma cruzi, an ensemble of linear ligand-based classifiers obtained by has been applied as initial screening filter, followed by docking into a homology model of the putrescine permease TcPAT12. 1,000 individual linear classifiers were inferred from a balanced dataset. Subsequently, different schemes were tested to combine the individual classifiers: MIN operator, average ranking, average score, average voting, with MIN operator leading to the best performance. The homology model was based on the arginine/agmatine antiporter (AdiC) from Escherichia coli as template. It showed 64% coverage of the entire query sequence and it was selected based on the normalized Discrete Optimized Protein Energy parameter and the GA341 score. The modeled structure had 96% in the allowed area of Ramachandran's plot, and none of the residues located in non-allowed regions were involved in the active site of the transporter. Positivity Predictive Value surfaces were applied to optimize the score thresholds to be used in the ligand-based virtual screening step: for that purpose Positivity Predictive Value was charted as a function of putative yields of active in the range 0.001-0.010 and the Se/Sp ratio. With a focus on drug repositioning opportunities, DrugBank and Sweetlead databases were subjected to screening. Among 8 hits, cinnarizine, a drug frequently prescribed for motion sickness and balance disorder, was tested against T. cruzi epimastigotes and amastigotes, confirming its trypanocidal effects and its inhibitory effects on putrescine uptake. Furthermore, clofazimine, an antibiotic with already proven trypanocidal effects, also displayed inhibitory effects on putrescine uptake. Two other hits, meclizine and butoconazole, also displayed trypanocidal effects (in the case of meclizine, against both epimastigotes and amastigotes), without inhibiting putrescine uptake.

Discovery of novel polyamine analogs with anti-protozoal activity by computer guided drug repositioning.

Chagas disease is a parasitic infection caused by the protozoa Trypanosoma cruzi that affects about 6 million people in Latin America. Despite its sanitary importance, there are currently only two drugs available for treatment: benznidazole and nifurtimox, both exhibiting serious adverse effects and limited efficacy in the chronic stage of the disease. Polyamines are ubiquitous to all living organisms where they participate in multiple basic functions such as biosynthesis of nucleic acids and proteins, proliferation and cell differentiation. T. cruzi is auxotroph for polyamines, which are taken up from the extracellular medium by efficient transporters and, to a large extent, incorporated into trypanothione (bis-glutathionylspermidine), the major redox cosubstrate of trypanosomatids. From a 268-compound database containing polyamine analogs with and without inhibitory effect on T. cruzi we have inferred classificatory models that were later applied in a virtual screening campaign to identify anti-trypanosomal compounds among drugs already used for other therapeutic indications (i.e. computer-guided drug repositioning) compiled in the DrugBank and Sweetlead databases. Five of the candidates identified with this strategy were evaluated in cellular models from different pathogenic trypanosomatids (T. cruzi wt, T. cruzi PAT12, T. brucei and Leishmania infantum), and in vitro models of aminoacid/polyamine transport assays and trypanothione synthetase inhibition assay. Triclabendazole, sertaconazole and paroxetine displayed inhibitory effects on the proliferation of T. cruzi (epimastigotes) and the uptake of putrescine by the parasite. They also interfered with the uptake of others aminoacids and the proliferation of infective T. brucei and L. infantum (promastigotes). Trypanothione synthetase was ruled out as molecular target for the anti-parasitic activity of these compounds.

Identification of levothyroxine antichagasic activity through computer-aided drug repurposing.

Cruzipain (Cz) is the major cysteine protease of the protozoan Trypanosoma cruzi, etiological agent of Chagas disease. A conformation-independent classifier capable of identifying Cz inhibitors was derived from a 163-compound dataset and later applied in a virtual screening campaign on the DrugBank database, which compiles FDA-approved and investigational drugs. 54 approved drugs were selected as candidates, 3 of which were acquired and tested on Cz and T. cruzi epimastigotes proliferation. Among them, levothyroxine, traditionally used in hormone replacement therapy in patients with hypothyroidism, showed dose-dependent inhibition of Cz and antiproliferative activity on the parasite.

Application of computer-aided drug repurposing in the search of new cruzipain inhibitors: discovery of amiodarone and bromocriptine inhibitory effects.

Cruzipain (Cz) is the major cystein protease of the protozoan Trypanosoma cruzi , etiological agent of Chagas disease. From a 163 compound data set, a 2D-classifier capable of identifying Cz inhibitors was obtained and applied in a virtual screening campaign on the DrugBank database, which compiles FDA-approved and investigational drugs. Fifty-four approved drugs were selected as candidates, four of which were acquired and tested on Cz and T. cruzi epimastigotes. Among them, the antiparkinsonian and antidiabetic drug bromocriptine and the antiarrhythmic amiodarone showed dose-dependent inhibition of Cz and antiproliferative activity on the parasite.