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.

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.

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.

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.

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.

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

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

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

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.

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.

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.

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 an acute treatment for the disease. We investigate whether compounds that bind the human ACE2 protein can interrupt SARS-CoV-2 replication without damaging ACE2’s natural enzymatic function. Initial compounds were screened for binding to ACE2 but little interruption of ACE2 enzymatic activity. This set of compounds was extended by application of quantitative structure-activity analysis, which resulted in 512 virtual hits for further confirmatory screening. A subsequent SARS-CoV-2 replication assay revealed that five of these compounds inhibit SARS-CoV-2 replication in human cells. Further effort is required to completely determine the antiviral mechanism of these compounds, but they serve as a strong starting point for both development of acute treatments for COVID-19 and research into the mechanism of infection.

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.