Innovative Multistage ML-QSAR Models for Malaria: From Data to Discovery.

Malaria presents a significant challenge to global public health, with around 247 million cases estimated to occur annually worldwide. The growing resistance of Plasmodium parasites to existing therapies underscores the urgent need for new and innovative antimalarial drugs. This study leveraged artificial intelligence (AI) to tackle this complex challenge. We developed multistage Machine Learning Quantitative Structure-Activity Relationship (ML-QSAR) models to effectively analyze large datasets and predict the efficacy of chemical compounds against multiple life cycle stages of Plasmodium parasites. We then selected 16 compounds for experimental evaluation, six of which showed at least dual-stage inhibitory activity and one inhibited all life cycle stages tested. Moreover, explainable AI (XAI) analysis provided insights into critical molecular features influencing model predictions, thereby enhancing our understanding of compound interactions. This study not only empowers the development of advanced predictive AI models but also accelerates the identification and optimization of potential antiplasmodial compounds.

Selective Bias Virtual Screening for Discovery of Promising Antimalarial Candidates targeting Plasmodium N-Myristoyltransferase.

Malaria remains a significant public health challenge, with Plasmodium vivax being the species responsible for the most prevalent form of the disease. Given the limited therapeutic options available, the search for new antimalarials against P. vivax is urgent. This study aims to identify new inhibitors for P. vivax N-myristoyltransferase (PvNMT), an essential drug target against malaria. Through a validated virtual screening campaign, we prioritized 23 candidates for further testing. In the yeast NMT system, seven compounds exhibit a potential inhibitor phenotype. In vitro antimalarial phenotypic assays confirmed the activity of four candidates while demonstrating an absence of cytotoxicity. Enzymatic assays reveal LabMol-394 as the most promising inhibitor, displaying selectivity against the parasite and a strong correlation within the yeast system. Furthermore, molecular dynamics simulations shed some light into its binding mode. This study constitutes a substantial contribution to the exploration of a selective quinoline scaffold and provides valuable insights into the development of new antimalarial candidates.

Identification of potential inhibitors of casein kinase 2 alpha of Plasmodium falciparum with potent in vitro activity.

Drug resistance to commercially available antimalarials is a major obstacle in malaria control and elimination, creating the need to find new antiparasitic compounds with novel mechanisms of action. The success of kinase inhibitors for oncological treatments has paved the way for the exploitation of protein kinases as drug targets in various diseases, including malaria. Casein kinases are ubiquitous serine/threonine kinases involved in a wide range of cellular processes such as mitotic checkpoint signaling, DNA damage response, and circadian rhythm. In Plasmodium, it is suggested that these protein kinases are essential for both asexual and sexual blood-stage parasites, reinforcing their potential as targets for multi-stage antimalarials. To identify new putative PfCK2α inhibitors, we utilized an in silico chemogenomic strategy involving virtual screening with docking simulations and quantitative structure-activity relationship predictions. Our investigation resulted in the discovery of a new quinazoline molecule (542), which exhibited potent activity against asexual blood stages and a high selectivity index (>100). Subsequently, we conducted chemical-genetic interaction analysis on yeasts with mutations in casein kinases. Our chemical-genetic interaction results are consistent with the hypothesis that 542 inhibits yeast Cka1, which has a hinge region with high similarity to PfCK2α. This finding is in agreement with our in silico results suggesting that 542 inhibits PfCK2α via hinge region interaction.

Transcriptomics-Guided In Silico Drug Repurposing: Identifying New Candidates with Dual-Stage Antiplasmodial Activity.

In tropical and subtropical areas, malaria stands as a profound public health challenge, causing an estimated 247 million cases worldwide annually. Given the absence of a viable vaccine, the timely and effective treatment of malaria remains a critical priority. However, the growing resistance of parasites to currently utilized drugs underscores the critical need for the identification of new antimalarial therapies. Here, we aimed to identify potential new drug candidates against Plasmodium falciparum, the main causative agent of malaria, by analyzing the transcriptomes of different life stages of the parasite and identifying highly expressed genes. We searched for genes that were expressed in all stages of the parasite's life cycle, including the asexual blood stage, gametocyte stage, liver stage, and sexual stages in the insect vector, using transcriptomics data from publicly available databases. From this analysis, we found 674 overlapping genes, including 409 essential ones. By searching through drug target databases, we discovered 70 potential drug targets and 75 associated bioactive compounds. We sought to expand this analysis to similar compounds to known drugs. So, we found a list of 1557 similar compounds, which we predicted as actives and inactives using previously developed machine learning models against five life stages of Plasmodium spp. From this analysis, two compounds were selected, and the reactions were experimentally evaluated. The compounds HSP-990 and silvestrol aglycone showed potent inhibitory activity at nanomolar concentrations against the P. falciparum 3D7 strain asexual blood stage. Moreover, silvestrol aglycone exhibited low cytotoxicity in mammalian cells, transmission-blocking potential, and inhibitory activity comparable to those of established antimalarials. These findings warrant further investigation of silvestrol aglycone as a potential dual-acting antimalarial and transmission-blocking candidate for malaria control.

Natural Compounds as Non-Nucleoside Inhibitors of Zika Virus Polymerase through Integration of In Silico and In Vitro Approaches.

Although the past epidemic of Zika virus (ZIKV) resulted in severe neurological consequences for infected infants and adults, there are still no approved drugs to treat ZIKV infection. In this study, we applied computational approaches to screen an in-house database of 77 natural and semi-synthetic compounds against ZIKV NS5 RNA-dependent RNA-polymerase (NS5 RdRp), an essential protein for viral RNA elongation during the replication process. For this purpose, we integrated computational approaches such as binding-site conservation, chemical space analysis and molecular docking. As a result, we prioritized nine virtual hits for experimental evaluation. Enzymatic assays confirmed that pedalitin and quercetin inhibited ZIKV NS5 RdRp with IC50 values of 4.1 and 0.5 µM, respectively. Moreover, pedalitin also displayed antiviral activity on ZIKV infection with an EC50 of 19.28 µM cell-based assays, with low toxicity in Vero cells (CC50 = 83.66 µM) and selectivity index of 4.34. These results demonstrate the potential of the natural compounds pedalitin and quercetin as candidates for structural optimization studies towards the discovery of new anti-ZIKV drug candidates.

PreS/MD: Predictor of Sensitization Hazard for Chemical Substances Released From Medical Devices.

In the United States, a pre-market regulatory submission for any medical device that comes into contact with either a patient or the clinical practitioner must include an adequate toxicity evaluation of chemical substances that can be released from the device during its intended use. These substances, also referred to as extractables and leachables, must be evaluated for their potential to induce sensitization/allergenicity, which traditionally has been done in animal assays such as the guinea pig maximization test (GPMT). However, advances in basic and applied science are continuously presenting opportunities to employ new approach methodologies, including computational methods which, when qualified, could replace animal testing methods to support regulatory submissions. Herein, we developed a new computational tool for rapid and accurate prediction of the GPMT outcome that we have named PreS/MD (predictor of sensitization for medical devices). To enable model development, we (1) collected, curated, and integrated the largest publicly available dataset for GPMT results; (2) succeeded in developing externally predictive (balanced accuracy of 70%-74% as evaluated by both 5-fold external cross-validation and testing of novel compounds) quantitative structure-activity relationships (QSAR) models for GPMT using machine learning algorithms, including deep learning; and (3) developed a publicly accessible web portal integrating PreS/MD models that can predict GPMT outcomes for any molecule of interest. We expect that PreS/MD will be used by both industry and regulatory scientists in medical device safety assessments and help replace, reduce, or refine the use of animals in toxicity testing. PreS/MD is freely available at https://presmd.mml.unc.edu/.

Update and elucidation of Plasmodium kinomes: Prioritization of kinases as potential drug targets for malaria.

Malaria is a tropical disease caused by Plasmodium spp. and transmitted by the bite of infected Anopheles mosquitoes. Protein kinases (PKs) play key roles in the life cycle of the etiological agent of malaria, turning these proteins attractive targets for antimalarial drug discovery campaigns. As part of an effort to understand parasite signaling functions, we report the results of a bioinformatics pipeline analysis of PKs of eight Plasmodium species. To date, no P. malariae and P. ovale kinome assemble has been conducted. We classified, curated and annotated predicted kinases to update P. falciparum, P. vivax, P. yoelii, P. berghei, P. chabaudi, and P. knowlesi kinomes published to date, as well as report for the first time the kinomes of P. malariae and P. ovale. Overall, from 76 to 97 PKs were identified among all Plasmodium spp. kinomes. Most of the kinases were assigned to seven of nine major kinase groups: AGC, CAMK, CMGC, CK1, STE, TKL, OTHER; and the Plasmodium-specific group FIKK. About 30% of kinases have been deeply classified into group, family and subfamily levels and only about 10% remained unclassified. Furthermore, updating and comparing the kinomes of P. vivax and P. falciparum allowed for the prioritization and selection of kinases as potential drug targets that could be explored for discovering new drugs against malaria. This integrated approach resulted in the selection of 37 protein kinases as potential targets and the identification of investigational compounds with moderate in vitro activity against asexual P. falciparum (3D7 and Dd2 strains) stages that could serve as starting points for the search of potent antimalarial leads in the future.

STopTox: An in Silico Alternative to Animal Testing for Acute Systemic and Topical Toxicity.

BACKGROUND: Modern chemical toxicology is facing a growing need to Reduce, Refine, and Replace animal tests (Russell 1959) for hazard identification. The most common type of animal assays for acute toxicity assessment of chemicals used as pesticides, pharmaceuticals, or in cosmetic products is known as a "6-pack" battery of tests, including three topical (skin sensitization, skin irritation and corrosion, and eye irritation and corrosion) and three systemic (acute oral toxicity, acute inhalation toxicity, and acute dermal toxicity) end points. METHODS: We compiled, curated, and integrated, to the best of our knowledge, the largest publicly available data sets and developed an ensemble of quantitative structure-activity relationship (QSAR) models for all six end points. All models were validated according to the Organisation for Economic Co-operation and Development (OECD) QSAR principles, using data on compounds not included in the training sets. RESULTS: In addition to high internal accuracy assessed by cross-validation, all models demonstrated an external correct classification rate ranging from 70% to 77%. We established a publicly accessible Systemic and Topical chemical Toxicity (STopTox) web portal (https://stoptox.mml.unc.edu/) integrating all developed models for 6-pack assays. CONCLUSIONS: We developed STopTox, a comprehensive collection of computational models that can be used as an alternative to in vivo 6-pack tests for predicting the toxicity hazard of small organic molecules. Models were established following the best practices for the development and validation of QSAR models. Scientists and regulators can use the STopTox portal to identify putative toxicants or nontoxicants in chemical libraries of interest. https://doi.org/10.1289/EHP9341.

Defining clinical outcome pathways.

Here, we propose a broad concept of 'Clinical Outcome Pathways' (COPs), which are defined as a series of key molecular and cellular events that underlie therapeutic effects of drug molecules. We formalize COPs as a chain of the following events: molecular initiating event (MIE) â†’ intermediate event(s) â†’ clinical outcome. We illustrate the concept with COP examples both for primary and alternative (i.e., drug repurposing) therapeutic applications. We also describe the elucidation of COPs for several drugs of interest using the publicly accessible Reasoning Over Biomedical Objects linked in Knowledge-Oriented Pathways (ROBOKOP) biomedical knowledge graph-mining tool. We propose that broader use of COP uncovered with the help of biomedical knowledge graph mining will likely accelerate drug discovery and repurposing efforts.

QSAR-Based Virtual Screening of Natural Products Database for Identification of Potent Antimalarial Hits.

With about 400,000 annual deaths worldwide, malaria remains a public health burden in tropical and subtropical areas, especially in low-income countries. Selection of drug-resistant Plasmodium strains has driven the need to explore novel antimalarial compounds with diverse modes of action. In this context, biodiversity has been widely exploited as a resourceful channel of biologically active compounds, as exemplified by antimalarial drugs such as quinine and artemisinin, derived from natural products. Thus, combining a natural product library and quantitative structure-activity relationship (QSAR)-based virtual screening, we have prioritized genuine and derivative natural compounds with potential antimalarial activity prior to in vitro testing. Experimental validation against cultured chloroquine-sensitive and multi-drug-resistant P. falciparum strains confirmed the potent and selective activity of two sesquiterpene lactones (LDT-597 and LDT-598) identified in silico. Quantitative structure-property relationship (QSPR) models predicted absorption, distribution, metabolism, and excretion (ADME) and physiologically based pharmacokinetic (PBPK) parameters for the most promising compound, showing that it presents good physiologically based pharmacokinetic properties both in rats and humans. Altogether, the in vitro parasite growth inhibition results obtained from in silico screened compounds encourage the use of virtual screening campaigns for identification of promising natural compound-based antimalarial molecules.

Chemogenomics and bioinformatics approaches for prioritizing kinases as drug targets for neglected tropical diseases.

Neglected tropical diseases (NTDs) are a group of twenty-one diseases classified by the World Health Organization that prevail in regions with tropical and subtropical climate and affect more than one billion people. There is an urgent need to develop new and safer drugs for these diseases. Protein kinases are a potential class of targets for developing new drugs against NTDs, since they play crucial role in many biological processes, such as signaling pathways, regulating cellular communication, division, metabolism and death. Bioinformatics is a field that aims to organize large amounts of biological data as well as develop and use tools for understanding and analyze them in order to produce meaningful information in a biological manner. In combination with chemogenomics, which analyzes chemical-biological interactions to screen ligands against selected targets families, these approaches can be used to stablish a rational strategy for prioritizing new drug targets for NTDs. Here, we describe how bioinformatics and chemogenomics tools can help to identify protein kinases and their potential inhibitors for the development of new drugs for NTDs. We present a review of bioinformatics tools and techniques that can be used to define an organisms kinome for drug prioritization, drug and target repurposing, multi-quinase inhibition approachs and selectivity profiling. We also present some successful examples of the application of such approaches in recent case studies.

Novel computational models offer alternatives to animal testing for assessing eye irritation and corrosion potential of chemicals.

Eye irritation and corrosion are fundamental considerations in developing chemicals to be used in or near the eye, from cleaning products to ophthalmic solutions. Unfortunately, animal testing is currently the standard method to identify compounds that cause eye irritation or corrosion. Yet, there is growing pressure on the part of regulatory agencies both in the USA and abroad to develop New Approach Methodologies (NAMs) that help reduce the need for animal testing and address unmet need to modernize safety evaluation of chemical hazards. In furthering the development and applications of computational NAMs in chemical safety assessment, in this study we have collected the largest expertly curated dataset of compounds tested for eye irritation and corrosion, and employed this data to build and validate binary and multi-classification Quantitative Structure-Activity Relationships (QSAR) models that can reliably assess eye irritation/corrosion potential of novel untested compounds. QSAR models were generated with Random Forest (RF) and Multi-Descriptor Read Across (MuDRA) machine learning (ML) methods, and validated using a 5-fold external cross-validation protocol. These models demonstrated high balanced accuracy (CCR of 0.68-0.88), sensitivity (SE of 0.61-0.84), positive predictive value (PPV of 0.65-0.90), specificity (SP of 0.56-0.91), and negative predictive value (NPV of 0.68-0.85). Overall, MuDRA models outperformed RF models and were applied to predict compounds' irritation/corrosion potential from the Inactive Ingredient Database, which contains components present in FDA-approved drug products, and from the Cosmetic Ingredient Database, the European Commission source of information on cosmetic substances. All models built and validated in this study are publicly available at the STopTox web portal (https://stoptox.mml.unc.edu/). These models can be employed as reliable tools for identifying potential eye irritant/corrosive compounds.

Pred-Skin: A Web Portal for Accurate Prediction of Human Skin Sensitizers.

Safety assessment is an essential component of the regulatory acceptance of industrial chemicals. Previously, we have developed a model to predict the skin sensitization potential of chemicals for two assays, the human patch test and murine local lymph node assay, and implemented this model in a web portal. Here, we report on the substantially revised and expanded freely available web tool, Pred-Skin version 3.0. This up-to-date version of Pred-Skin incorporates multiple quantitative structure-activity relationship (QSAR) models developed with in vitro, in chemico, and mice and human in vivo data, integrated into a consensus naïve Bayes model that predicts human effects. Individual QSAR models were generated using skin sensitization data derived from human repeat insult patch tests, human maximization tests, and mouse local lymph node assays. In addition, data for three validated alternative methods, the direct peptide reactivity assay, KeratinoSens, and the human cell line activation test, were employed as well. Models were developed using open-source tools and rigorously validated according to the best practices of QSAR modeling. Predictions obtained from these models were then used to build a naïve Bayes model for predicting human skin sensitization with the following external prediction accuracy: correct classification rate (89%), sensitivity (94%), positive predicted value (91%), specificity (84%), and negative predicted value (89%). As an additional assessment of model performance, we identified 11 cosmetic ingredients known to cause skin sensitization but were not included in our training set, and nine of them were accurately predicted as sensitizers by our models. Pred-Skin can be used as a reliable alternative to animal tests for predicting human skin sensitization.

Artificial Intelligence Applied to the Rapid Identification of New Antimalarial Candidates with Dual-Stage Activity.

Increasing reports of multidrug-resistant malaria parasites urge the discovery of new effective drugs with different chemical scaffolds. Protein kinases play a key role in many cellular processes such as signal transduction and cell division, making them interesting targets in many diseases. Protein kinase 7 (PK7) is an orphan kinase from the Plasmodium genus, essential for the sporogonic cycle of these parasites. Here, we applied a robust and integrative artificial intelligence-assisted virtual-screening (VS) approach using shape-based and machine learning models to identify new potential PK7 inhibitors with in vitro antiplasmodial activity. Eight virtual hits were experimentally evaluated, and compound LabMol-167 inhibited ookinete conversion of Plasmodium berghei and blood stages of Plasmodium falciparum at nanomolar concentrations with low cytotoxicity in mammalian cells. As PK7 does not have an essential role in the Plasmodium blood stage and our virtual screening strategy aimed for both PK7 and blood-stage inhibition, we conducted an in silico target fishing approach and propose that this compound might also inhibit P. falciparum PK5, acting as a possible dual-target inhibitor. Finally, docking studies of LabMol-167 with P. falciparum PK7 and PK5 proteins highlighted key interactions for further hit-to lead optimization.

Unveiling the Kinomes of Leishmania infantum and L. braziliensis Empowers the Discovery of New Kinase Targets and Antileishmanial Compounds.

Leishmaniasis is a neglected tropical disease caused by parasites of the genus Leishmania (NTD) endemic in 98 countries. Although some drugs are available, current treatments deal with issues such as toxicity, low efficacy, and emergence of resistance. Therefore, there is an urgent need to identify new targets for the development of new antileishmanial drugs. Protein kinases (PKs), which play an essential role in many biological processes, have become potential drug targets for many parasitic diseases. A refined bioinformatics pipeline was applied in order to define and compare the kinomes of L. infantum and L. braziliensis, species that cause cutaneous and visceral manifestations of leishmaniasis in the Americas, the latter being potentially fatal if untreated. Respectively, 224 and 221 PKs were identified in L. infantum and L. braziliensis overall. Almost all unclassified eukaryotic PKs were assigned to six of nine major kinase groups and, consequently, most have been classified into family and subfamily. Furthermore, revealing the kinomes for both Leishmania species allowed for the prioritization of potential drug targets that could be explored for discovering new drugs against leishmaniasis. Finally, we used a drug repurposing approach and prioritized seven approved drugs and investigational compounds to be experimentally tested against Leishmania. Trametinib and NMS-1286937 inhibited the growth of L. infantum and L. braziliensis promastigotes and amastigotes and therefore might be good candidates for the drug repurposing pipeline.

Computationally-guided drug repurposing enables the discovery of kinase targets and inhibitors as new schistosomicidal agents.

The development of novel therapeutics is urgently required for diseases where existing treatments are failing due to the emergence of resistance. This is particularly pertinent for parasitic infections of the tropics and sub-tropics, referred to collectively as neglected tropical diseases, where the commercial incentives to develop new drugs are weak. One such disease is schistosomiasis, a highly prevalent acute and chronic condition caused by a parasitic helminth infection, with three species of the genus Schistosoma infecting humans. Currently, a single 40-year old drug, praziquantel, is available to treat all infective species, but its use in mass drug administration is leading to signs of drug-resistance emerging. To meet the challenge of developing new therapeutics against this disease, we developed an innovative computational drug repurposing pipeline supported by phenotypic screening. The approach highlighted several protein kinases as interesting new biological targets for schistosomiasis as they play an essential role in many parasite's biological processes. Focusing on this target class, we also report the first elucidation of the kinome of Schistosoma japonicum, as well as updated kinomes of S. mansoni and S. haematobium. In comparison with the human kinome, we explored these kinomes to identify potential targets of existing inhibitors which are unique to Schistosoma species, allowing us to identify novel targets and suggest approved drugs that might inhibit them. These include previously suggested schistosomicidal agents such as bosutinib, dasatinib, and imatinib as well as new inhibitors such as vandetanib, saracatinib, tideglusib, alvocidib, dinaciclib, and 22 newly identified targets such as CHK1, CDC2, WEE, PAKA, MEK1. Additionally, the primary and secondary targets in Schistosoma of those approved drugs are also suggested, allowing for the development of novel therapeutics against this important yet neglected disease.

Chemoproteomic identification of molecular targets of antifungal prototypes, thiosemicarbazide and a camphene derivative of thiosemicarbazide, in Paracoccidioides brasiliensis.

Paracoccidioidomycosis (PCM) is a neglected human systemic disease caused by species of the genus Paracoccidioides. The disease attacks the host's lungs and may disseminate to many other organs. Treatment involves amphotericin B, sulfadiazine, trimethoprim-sulfamethoxazole, itraconazole, ketoconazole, or fluconazole. The treatment duration is usually long, from 6 months to 2 years, and many adverse effects may occur in relation to the treatment; co-morbidities and poor treatment adherence have been noted. Therefore, the discovery of more effective and less toxic drugs is needed. Thiosemicarbazide (TSC) and a camphene derivative of thiosemicarbazide (TSC-C) were able to inhibit P. brasiliensis growth at a low dosage and were not toxic to fibroblast cells. In order to investigate the mode of action of those compounds, we used a chemoproteomic approach to determine which fungal proteins were bound to each of these compounds. The compounds were able to inhibit the activities of the enzyme formamidase and interfered in P. brasiliensis dimorphism. In comparison with the transcriptomic and proteomic data previously obtained by our group, we determined that TSC and TSC-C were multitarget compounds that exerted effects on the electron-transport chain and cell cycle regulation, increased ROS formation, inhibited proteasomes and peptidases, modulated glycolysis, lipid, protein and carbohydrate metabolisms, and caused suppressed the mycelium to yeast transition.

The A-Z of Zika drug discovery.

Despite the recent outbreak of Zika virus (ZIKV), there are still no approved treatments, and early-stage compounds are probably many years away from approval. A comprehensive A-Z review of the recent advances in ZIKV drug discovery efforts is presented, highlighting drug repositioning and computationally guided compounds, including discovered viral and host cell inhibitors. Promising ZIKV molecular targets are also described and discussed, as well as targets belonging to the host cell, as new opportunities for ZIKV drug discovery. All this knowledge is not only crucial to advancing the fight against the Zika virus and other flaviviruses but also helps us prepare for the next emerging virus outbreak to which we will have to respond.

Computational drug discovery for the Zika virus

Few Zika virus (ZIKV) outbreaks had been reported since its first detection in 1947, until the recent epidemics occurred in South America (2014/2015) and expeditiously became a global public health emergency. This arbovirus reached 0.5-1.3 million cases of ZIKV infection in Brazil in 2015 and rapidly spread in new geographic areas such as the Americas. Despite the mild symptoms of the Zika fever, the major concern is related to the related severe neurological disorders, especially microcephaly in newborns. Advances in ZIKV drug discovery have been made recently and constitute promising approaches to ZIKV treatment. In this review, we summarize current computational drug discovery efforts and their applicability to discovery of anti-ZIKV drugs. Lastly, we present successful examples of the use of computational approaches to ZIKV drug discovery.