Modeling interactions between Heparan sulfate and proteins based on the Heparan sulfate microarray analysis.

Heparan sulfate (HS), a sulfated polysaccharide abundant in the extracellular matrix, plays pivotal roles in various physiological and pathological processes by interacting with proteins. Investigating the binding selectivity of HS oligosaccharides to target proteins is essential, but the exhaustive inclusion of all possible oligosaccharides in microarray experiments is impractical. To address this challenge, we present a hybrid pipeline that integrates microarray and in silico techniques to design oligosaccharides with desired protein affinity. Using fibroblast growth factor 2 (FGF2) as a model protein, we assembled an in-house dataset of HS oligosaccharides on microarrays and developed two structural representations: a standard representation with all atoms explicit and a simplified representation with disaccharide units as "quasi-atoms." Predictive Quantitative Structure-Activity Relationship (QSAR) models for FGF2 affinity were developed using the Random Forest (RF) algorithm. The resulting models, considering the applicability domain, demonstrated high predictivity, with a correct classification rate of 0.81-0.80 and improved positive predictive values (PPV) up to 0.95. Virtual screening of 40 new oligosaccharides using the simplified model identified 15 computational hits, 11 of which were experimentally validated for high FGF2 affinity. This hybrid approach marks a significant step toward the targeted design of oligosaccharides with desired protein interactions, providing a foundation for broader applications in glycobiology.

STOPLIGHT: A Hit Scoring Calculator.

We introduce STOPLIGHT, a web portal to assist medicinal chemists in prioritizing hits from screening campaigns and the selection of compounds for optimization. STOPLIGHT incorporates services to assess 6 physiochemical and structural properties, 6 assay liabilities, and 11 pharmacokinetic properties, for any small molecule represented by its SMILES string. We briefly describe each service and illustrate the utility of this portal with a case study. The STOPLIGHT portal provides a user-friendly tool to guide hit selection in early drug discovery campaigns, whereby compounds with unfavorable properties can be quickly recognized and eliminated.

Pharmacokinetics Profiler (PhaKinPro): Model Development, Validation, and Implementation as a Web Tool for Triaging Compounds with Undesired Pharmacokinetics Profiles.

Computational models that predict pharmacokinetic properties are critical to deprioritize drug candidates that emerge as hits in high-throughput screening campaigns. We collected, curated, and integrated a database of compounds tested in 12 major end points comprising over 10,000 unique molecules. We then employed these data to build and validate binary quantitative structure-activity relationship (QSAR) models. All trained models achieved a correct classification rate above 0.60 and a positive predictive value above 0.50. To illustrate their utility in drug discovery, we used these models to predict the pharmacokinetic properties for drugs in the NCATS Inxight Drugs database. In addition, we employed the developed models to predict the pharmacokinetic properties of all compounds in the DrugBank. All models described in this paper have been integrated and made publicly available via the PhaKinPro Web-portal that can be accessed at https://phakinpro.mml.unc.edu/.

Cheminformatics-Guided Cell-Free Exploration of Peptide Natural Products.

There have been significant advances in the flexibility and power of in vitro cell-free translation systems. The increasing ability to incorporate noncanonical amino acids and complement translation with recombinant enzymes has enabled cell-free production of peptide-based natural products (NPs) and NP-like molecules. We anticipate that many more such compounds and analogs might be accessed in this way. To assess the peptide NP space that is directly accessible to current cell-free technologies, we developed a peptide parsing algorithm that breaks down peptide NPs into building blocks based on ribosomal translation logic. Using the resultant data set, we broadly analyze the biophysical properties of these privileged compounds and perform a retrobiosynthetic analysis to predict which peptide NPs could be directly synthesized in augmented cell-free translation reactions. We then tested these predictions by preparing a library of highly modified peptide NPs. Two macrocyclases, PatG and PCY1, were used to effect the head-to-tail macrocyclization of candidate NPs. This retrobiosynthetic analysis identified a collection of high-priority building blocks that are enriched throughout peptide NPs, yet they had not previously been tested in cell-free translation. To expand the cell-free toolbox into this space, we established, optimized, and characterized the flexizyme-enabled ribosomal incorporation of piperazic acids. Overall, these results demonstrate the feasibility of cell-free translation for peptide NP total synthesis while expanding the limits of the technology. This work provides a novel computational tool for exploration of peptide NP chemical space, that could be expanded in the future to allow design of ribosomal biosynthetic pathways for NPs and NP-like molecules.

ExEmPLAR (Extracting, Exploring, and Embedding Pathways Leading to Actionable Research): a user-friendly interface for knowledge graph mining.

SUMMARY: Knowledge graphs are being increasingly used in biomedical research to link large amounts of heterogenous data and facilitate reasoning across diverse knowledge sources. Wider adoption and exploration of knowledge graphs in the biomedical research community is limited by requirements to understand the underlying graph structure in terms of entity types and relationships, represented as nodes and edges, respectively, and learn specialized query languages for graph mining and exploration. We have developed a user-friendly interface dubbed ExEmPLAR (Extracting, Exploring, and Embedding Pathways Leading to Actionable Research) to aid reasoning over biomedical knowledge graphs and assist with data-driven research and hypothesis generation. We explain the key functionalities of ExEmPLAR and demonstrate its use with a case study considering the relationship of Trypanosoma cruzi, the etiological agent of Chagas disease, to frequently associated cardiovascular conditions. AVAILABILITY AND IMPLEMENTATION: ExEmPLAR is freely accessible at https://www.exemplar.mml.unc.edu/. For code and instructions for the using the application, see: https://github.com/beasleyjonm/AOP-COP-Path-Extractor.

Lies and Liabilities: Computational Assessment of High-Throughput Screening Hits to Identify Artifact Compounds.

Hits from high-throughput screening (HTS) of chemical libraries are often false positives due to their interference with assay detection technology. In response, we generated the largest publicly available library of chemical liabilities and developed "Liability Predictor," a free web tool to predict HTS artifacts. More specifically, we generated, curated, and integrated HTS data sets for thiol reactivity, redox activity, and luciferase (firefly and nano) activity and developed and validated quantitative structure-interference relationship (QSIR) models to predict these nuisance behaviors. The resulting models showed 58-78% external balanced accuracy for 256 external compounds per assay. QSIR models developed and validated herein identify nuisance compounds among experimental hits more reliably than do popular PAINS filters. Both the models and the curated data sets were implemented in "Liability Predictor," publicly available at https://liability.mml.unc.edu/. "Liability Predictor" may be used as part of chemical library design or for triaging HTS hits.

Multitask learning-driven identification of novel antitrypanosomal compounds.

Background: Chagas disease and human African trypanosomiasis cause substantial death and morbidity, particularly in low- and middle-income countries, making the need for novel drugs urgent. Methodology & results: Therefore, an explainable multitask pipeline to profile the activity of compounds against three trypanosomes (Trypanosoma brucei brucei, Trypanosoma brucei rhodesiense and Trypanosoma cruzi) were created. These models successfully discovered four new experimental hits (LC-3, LC-4, LC-6 and LC-15). Among them, LC-6 showed promising results, with IC50 values ranging 0.01-0.072 µM and selectivity indices >10,000. Conclusion: These results demonstrate that the multitask protocol offers predictivity and interpretability in the virtual screening of new antitrypanosomal compounds and has the potential to improve hit rates in Chagas and human African trypanosomiasis projects.

Cheminformatics-driven discovery of hit compounds against Paracoccidioides spp.

Aims: The development of safe and effective therapies for treating paracoccidioidomycosis using computational strategies were employed to discover anti-Paracoccidioides compounds. Materials & methods: We 1) collected, curated and integrated the largest library of compounds tested against Paracoccidioides spp.; 2) employed a similarity search to virtually screen the ChemBridge database and select nine compounds for experimental evaluation; 3) performed an experimental evaluation to determine the minimum inhibitory concentration and minimum fungicidal concentration as well as cytotoxicity; and 4) employed computational tools to identify potential targets for the most active compounds. Seven compounds presented activity against Paracoccidioides spp. Conclusion: These compounds are new hits with a predicted mechanisms of action, making them potentially attractive to develop new compounds.

Praemonitus praemunitus: can we forecast and prepare for future viral disease outbreaks?

Understanding the origins of past and present viral epidemics is critical in preparing for future outbreaks. Many viruses, including SARS-CoV-2, have led to significant consequences not only due to their virulence, but also because we were unprepared for their emergence. We need to learn from large amounts of data accumulated from well-studied, past pandemics and employ modern informatics and therapeutic development technologies to forecast future pandemics and help minimize their potential impacts. While acknowledging the complexity and difficulties associated with establishing reliable outbreak predictions, herein we provide a perspective on the regions of the world that are most likely to be impacted by future outbreaks. We specifically focus on viruses with epidemic potential, namely SARS-CoV-2, MERS-CoV, DENV, ZIKV, MAYV, LASV, noroviruses, influenza, Nipah virus, hantaviruses, Oropouche virus, MARV, and Ebola virus, which all require attention from both the public and scientific community to avoid societal catastrophes like COVID-19. Based on our literature review, data analysis, and outbreak simulations, we posit that these future viral epidemics are unavoidable, but that their societal impacts can be minimized by strategic investment into basic virology research, epidemiological studies of neglected viral diseases, and antiviral drug discovery.

Small molecule antiviral compound collection (SMACC): A comprehensive, highly curated database to support the discovery of broad-spectrum antiviral drug molecules.

Diseases caused by new viruses cost thousands if not millions of human lives and trillions of dollars. We have identified, collected, curated, and integrated all chemogenomics data from ChEMBL for 13 emerging viruses that hold the greatest potential threat to global human health. By identifying and solving several challenges related to data annotation accuracy, we developed a highly curated and thoroughly annotated database of compounds tested in both phenotypic and target-based assays for these viruses that we dubbed SMACC (Small Molecule Antiviral Compound Collection). The pilot version of the SMACC database contains over 32,500 entries for 13 viruses. By analyzing data in SMACC, we have identified ∼50 compounds with polyviral inhibition profile, mostly covering flavi- and coronaviruses. The SMACC database may serve as a reference for virologists and medicinal chemists working on the development of novel BSA agents in preparation for future viral outbreaks. SMACC is publicly available at https://smacc.mml.unc.edu.

Selene-Ethylenelacticamides and N-Aryl-Propanamides as Broad-Spectrum Leishmanicidal Agents.

The World Health Organization classifies Leishmania as one of the 17 "neglected diseases" that burden tropical and sub-tropical climate regions with over half a million diagnosed cases each year. Despite this, currently available anti-leishmania drugs have high toxicity and the potential to be made obsolete by parasite drug resistance. We chose to analyze organoselenides for leishmanicidal potential given the reduced toxicity inherent to selenium and the displayed biological activity of organoselenides against Leishmania. Thus, the biological activities of 77 selenoesters and their N-aryl-propanamide derivatives were predicted using robust in silico models of Leishmania infantum, Leishmania amazonensis, Leishmania major, and Leishmania (Viannia) braziliensis. The models identified 28 compounds with >60% probability of demonstrating leishmanicidal activity against L. infantum, and likewise, 26 for L. amazonesis, 25 for L. braziliensis, and 23 for L. major. The in silico prediction of ADMET properties suggests high rates of oral absorption and good bioavailability for these compounds. In the in silico toxicity evaluation, only seven compounds showed signs of toxicity in up to one or two parameters. The methodology was corroborated with the ensuing experimental validation, which evaluated the inhibition of the Promastigote form of the Leishmania species under study. The activity of the molecules was determined by the IC50 value (µM); IC50 values < 20 µM indicated better inhibition profiles. Sixteen compounds were synthesized and tested for their activity. Eight molecules presented IC50 values < 20 µM for at least one of the Leishmania species under study, with compound NC34 presenting the strongest parasite inhibition profile. Furthermore, the methodology used was effective, as many of the compounds with the highest probability of activity were confirmed by the in vitro tests performed.

Four diterpenes identified in silico were isolated from Hyptidinae and demonstrated in vitro activity against Mycobacterium tuberculosis.

Plants of Hyptidinae subtribe (Lamiaceae - family), as Mesosphaerum sidifolium, are a source of bioactive molecules. In the search for new drug candidates, we perform chemical characterization of diterpenes isolated from the aerial parts of M. sidifolium was carried out with uni- and bidimensional NMR spectral data, and evaluate in silico through the construction of a predictive model followed by in vitro testing Mycobacterium tuberculosis and Mycobacterium smegmatis. Resulted in the isolation of four components: Pomiferin D (1), Salviol (2), Pomiferin E (3) and 2α-hydroxysugiol (4), as well as two phenolic compounds, rosmarinic and caffeic acids. In silico model identified 48 diterpenes likely to have biological activity against M. tuberculosis. The diterpenes isolated were tested in vitro against M. tuberculosis demonstrating MIC = 125 µM for 4 and 1, while 2 and 3 -MIC = 250 µM. These compounds did not show biological activity at these concentrations for M. smegmatis.

Discovery of New Zika Protease and Polymerase Inhibitors through the Open Science Collaboration Project OpenZika.

The Zika virus (ZIKV) is a neurotropic arbovirus considered a global threat to public health. Although there have been several efforts in drug discovery projects for ZIKV in recent years, there are still no antiviral drugs approved to date. Here, we describe the results of a global collaborative crowdsourced open science project, the OpenZika project, from IBM's World Community Grid (WCG), which integrates different computational and experimental strategies for advancing a drug candidate for ZIKV. Initially, molecular docking protocols were developed to identify potential inhibitors of ZIKV NS5 RNA-dependent RNA polymerase (NS5 RdRp), NS3 protease (NS2B-NS3pro), and NS3 helicase (NS3hel). Then, a machine learning (ML) model was built to distinguish active vs inactive compounds for the cytoprotective effect against ZIKV infection. We performed three independent target-based virtual screening campaigns (NS5 RdRp, NS2B-NS3pro, and NS3hel), followed by predictions by the ML model and other filters, and prioritized a total of 61 compounds for further testing in enzymatic and phenotypic assays. This yielded five non-nucleoside compounds which showed inhibitory activity against ZIKV NS5 RdRp in enzymatic assays (IC50 range from 0.61 to 17 µM). Two compounds thermally destabilized NS3hel and showed binding affinity in the micromolar range (Kd range from 9 to 35 µM). Moreover, the compounds LabMol-301 inhibited both NS5 RdRp and NS2B-NS3pro (IC50 of 0.8 and 7.4 µM, respectively) and LabMol-212 thermally destabilized the ZIKV NS3hel (Kd of 35 µM). Both also protected cells from death induced by ZIKV infection in in vitro cell-based assays. However, while eight compounds (including LabMol-301 and LabMol-212) showed a cytoprotective effect and prevented ZIKV-induced cell death, agreeing with our ML model for prediction of this cytoprotective effect, no compound showed a direct antiviral effect against ZIKV. Thus, the new scaffolds discovered here are promising hits for future structural optimization and for advancing the discovery of further drug candidates for ZIKV. Furthermore, this work has demonstrated the importance of the integration of computational and experimental approaches, as well as the potential of large-scale collaborative networks to advance drug discovery projects for neglected diseases and emerging viruses, despite the lack of available direct antiviral activity and cytoprotective effect data, that reflects on the assertiveness of the computational predictions. The importance of these efforts rests with the need to be prepared for future viral epidemic and pandemic outbreaks.

Integrated approach to elucidate metal-implant related adverse outcome pathways.

Exogenous metal particles and ions from implant devices are known to cause severe toxic events with symptoms ranging from adverse local tissue reactions to systemic toxicities, potentially leading to the development of cancers, heart conditions, and neurological disorders. Toxicity mechanisms, also known as Adverse Outcome Pathways (AOPs), that explain these metal-induced toxicities are severely understudied. Therefore, we deployed in silico structure- and knowledge-based approaches to identify proteome-level perturbations caused by metals and pathways that link these events to human diseases. We captured 177 structure-based, 347 knowledge-based, and 402 imputed metal-gene/protein relationships for chromium, cobalt, molybdenum, nickel, and titanium. We prioritized 72 proteins hypothesized to directly contact implant surfaces and contribute to adverse outcomes. Results of this exploratory analysis were formalized as structured AOPs. We considered three case studies reflecting the following possible situations: (i) the metal-protein-disease relationship was previously known; (ii) the metal-protein, protein-disease, and metal-disease relationships were individually known but were not linked (as a unified AOP); and (iii) one of three relationships was unknown and was imputed by our methods. These situations were illustrated by case studies on nickel-induced allergy/hypersensitivity, cobalt-induced heart failure, and titanium-induced periprosthetic osteolysis, respectively. All workflows, data, and results are freely available in https://github.com/DnlRKorn/Knowledge_Based_Metallomics/. An interactive view of select data is available at the ROBOKOP Neo4j Browser at http://robokopkg.renci.org/browser/.

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/.

MolPredictX: Online Biological Activity Predictions by Machine Learning Models.

Here we report the development of MolPredictX, an innovate and freely accessible web interface for biological activity predictions of query molecules. MolPredictX utilizes in-house QSAR models to provide 27 qualitative predictions (active or inactive), and quantitative probabilities for bioactivity against parasitic (Trypanosoma and Leishmania), viral (Dengue, Sars-CoV and Hepatitis C), pathogenic yeast (Candida albicans), bacterial (Salmonella enterica and Escherichia coli), and Alzheimer disease enzymes. In this article, we introduce the methodology and usability of this webtool, highlighting its potential role in the development of new drugs against a variety of diseases. MolPredictX is undergoing continuous development and is freely available at https://www.molpredictx.ufpb.br/.

Allosteric Binders of ACE2 Are Promising Anti-SARS-CoV-2 Agents.

The COVID-19 pandemic has had enormous health, economic, and social consequences. Vaccines have been successful in reducing rates of infection and hospitalization, but there is still a need for acute treatment of the disease. We investigate whether compounds that bind the human angiotensin-converting enzyme 2 (ACE2) protein can decrease SARS-CoV-2 replication without impacting ACE2's natural enzymatic function. Initial screening of a diversity library resulted in hit compounds active in an ACE2-binding assay, which showed little inhibition of ACE2 enzymatic activity (116 actives, success rate ∼4%), suggesting they were allosteric binders. Subsequent application of in silico techniques boosted success rates to ∼14% and resulted in 73 novel confirmed ACE2 binders with K d values as low as 6 nM. A subsequent SARS-CoV-2 assay revealed that five of these compounds inhibit the viral life cycle in human cells. Further effort is required to completely elucidate the antiviral mechanism of these ACE2-binders, but they present a valuable starting point for both the development of acute treatments for COVID-19 and research into the host-directed therapy.

Small Molecule Antiviral Compound Collection (SMACC): a database to support the discovery of broad-spectrum antiviral drug molecules.

Diseases caused by new viruses costs thousands if not millions of human lives and trillions of dollars in damage to the global economy. Despite the rapid development of vaccines for SARS-CoV-2, the lack of small molecule antiviral drugs that work against multiple viral families (broad-spectrum antivirals; BSAs) has left the entire world’s human population vulnerable to the infection between the beginning of the outbreak and the widespread availability of vaccines. Developing BSAs is an attractive, yet challenging, approach that could prevent the next, inevitable, viral outbreak from becoming a global catastrophe. To explore whether historical medicinal chemistry efforts suggest the possibility of discovering novel BSAs, we (i) identified, collected, curated, and integrated all chemical bioactivity data available in ChEMBL for molecules tested in respective assays for 13 emerging viruses that, based on published literature, hold the greatest potential threat to global human health; (ii) identified and solved the challenges related to data annotation accuracy including assay description ambiguity, missing cell or target information, and incorrect BioAssay Ontology (BAO) annotations; (iii) developed a highly curated and thoroughly annotated database of compounds tested in both phenotypic (21,392 entries) and target-based (11,123 entries) assays for these viruses; and (iv) identified a subset of compounds showing BSA activity. For the latter task, we eliminated inconclusive and annotated duplicative entries by checking the concordance between multiple assay results and identified eight compounds active against 3-4 viruses from the phenotypic data, 16 compounds active against two viruses from the target-based data, and 35 compounds active in at least one phenotypic and one target-based assay. The pilot version of our SMACC (Small Molecule Antiviral Compound Collection) database contains over 32,500 entries for 13 viruses. Our analysis indicates that previous research yielded very small number of BSA compounds. We posit that focused and coordinated efforts strategically targeting the discovery of such agents must be established and maintained going forward. The SMACC database publicly available at https://smacc.mml.unc.edu may serve as a reference for virologists and medicinal chemists working on the development of novel BSA agents in preparation for future viral outbreaks.

Natural Products from Annonaceae as Potential Antichagasic Agents.

Chagas disease, a neglected tropical disease, is endemic in 21 Latin American countries and particularly prevalent in Brazil. Chagas disease has drawn more attention in recent years due to its expansion into non-endemic areas. The aim of this work was to computationally identify and experimentally validate the natural products from an Annonaceae family as antichagasic agents. Through the ligand-based virtual screening, we identified 57 molecules with potential activity against the epimastigote form of T. cruzi. Then, 16 molecules were analyzed in the in vitro study, of which, six molecules displayed previously unknown antiepimastigote activity. We also evaluated these six molecules for trypanocidal activity. We observed that all six molecules have potential activity against the amastigote form, but no molecules were active against the trypomastigote form. 13-Epicupressic acid seems to be the most promising, as it was predicted as an active compound in the in silico study against the amastigote form of T. cruzi, in addition to having in vitro activity against the epimastigote form.

Conserved coronavirus proteins as targets of broad-spectrum antivirals.

Coronaviruses are a class of single-stranded, positive-sense RNA viruses that have caused three major outbreaks over the past two decades: Middle East respiratory syndrome-related coronavirus (MERS-CoV), severe acute respiratory syndrome coronavirus (SARS-CoV), and severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). All outbreaks have been associated with significant morbidity and mortality. In this study, we have identified and explored conserved binding sites in the key coronavirus proteins for the development of broad-spectrum direct acting anti-coronaviral compounds and validated the significance of this conservation for drug discovery with existing experimental data. We have identified four coronaviral proteins with highly conserved binding site sequence and 3D structure similarity: PLpro, Mpro, nsp10-nsp16 complex(methyltransferase), and nsp15 endoribonuclease. We have compiled all available experimental data for known antiviral medications inhibiting these targets and identified compounds active against multiple coronaviruses. The identified compounds representing potential broad-spectrum antivirals include: GC376, which is active against six viral Mpro (out of six tested, as described in research literature); mycophenolic acid, which is active against four viral PLpro (out of four); and emetine, which is active against four viral RdRp (out of four). The approach described in this study for coronaviruses, which combines the assessment of sequence and structure conservation across a viral family with the analysis of accessible chemical structure - antiviral activity data, can be explored for the development of broad-spectrum drugs for multiple viral families.