Visualization of drug target interactions in the contexts of pathways and networks with ReactomeFIViz

scRank infers drug-responsive cell types from untreated scRNA-seq data using a target-perturbed gene regulatory network

Cells respond divergently to drugs due to the heterogeneity among cell populations. Thus, it is crucial to identify drug-responsive cell populations in order to accurately elucidate the mechanism of drug action, which is still a great challenge. Here, we address this problem with scRank, which employs a target-perturbed gene regulatory network to rank drug-responsive cell populations via in silico drug perturbations using untreated single-cell transcriptomic data. We benchmark scRank on simulated and real datasets, which shows the superior performance of scRank over existing methods. When applied to medulloblastoma and major depressive disorder datasets, scRank identifies drug-responsive cell types that are consistent with the literature. Moreover, scRank accurately uncovers the macrophage subpopulation responsive to tanshinone IIA and its potential targets in myocardial infarction, with experimental validation. In conclusion, scRank enables the inference of drug-responsive cell types using untreated single-cell data, thus providing insights into the cellular-level impacts of therapeutic interventions.

Enhancing the anticancer potential of metformin: fabrication of efficient nanospanlastics, in vitro cytotoxic studies on HEP-2 cells and reactome enhanced pathway analysis

Metformin (MET), an oral antidiabetic drug, was reported to possess promising anticancer effects. We hypothesized that MET encapsulation in unique nanospanlastics would enhance its anticancer potential against HEP-2 cells. Our results showed the successful fabrication of Nano-MET spanlastics (d = 232.10 ± 0.20 nm; PDI = 0.25 ± 0.11; zeta potential = (-) 44.50 ± 0.96; drug content = 99.90 ± 0.11 and entrapment efficiency = 88.01 ± 2.50%). MTT assay revealed the enhanced Nano-MET cytotoxicity over MET with a calculated IC50 of 50 μg/mL and > 500 μg/mL, respectively. Annexin V/PI apoptosis assay showed that Nano-MET significantly decreased the percentage of live cells from 95.49 to 93.70 compared to MET and increased the percentage of cells arrested in the G0/G1 phase by 8.38%. Moreover, Nano-MET downregulated BCL-2 and upregulated BAX protein levels by 1.57 and 1.88 folds, respectively. RT-qPCR revealed that Nano-MET caused a significant 13.75, 4.15, and 2.23-fold increase in caspase-3, -8, and - 9 levels as well as a 100 and 43.47-fold decrease in cyclin D1 and mTOR levels, respectively. The proliferation marker Ki67 immunofluorescent staining revealed a 3-fold decrease in positive cells in Nano-MET compared to the control. Utilizing the combined Pathway-Enrichment Analysis (PEA) and Reactome analysis indicated high enrichment of certain pathways including nucleotides metabolism, Nudix-type hydrolase enzymes, carbon dioxide hydration, hemostasis, and the innate immune system. In summary, our results confirm MET cytotoxicity enhancement by its encapsulation in nanospanlastics. We also highlight, using PEA, that MET can modulate multiple pathways implicated in carcinogenesis.

Illuminating Dark Proteins using Reactome Pathways

Limited knowledge about a substantial portion of protein coding genes, known as "dark" proteins, hinders our understanding of their functions and potential therapeutic applications. To address this, we leveraged Reactome, the most comprehensive, open source, open-access pathway knowledgebase, to contextualize dark proteins within biological pathways. By integrating multiple resources and employing a random forest classifier trained on 106 protein/gene pairwise features, we predicted functional interactions between dark proteins and Reactome-annotated proteins. We then developed three scores to measure the interactions between dark proteins and Reactome pathways, utilizing enrichment analysis and fuzzy logic simulations. Correlation analysis of these scores with an independent single-cell RNA sequencing dataset provided supporting evidence for this approach. Furthermore, systematic natural language processing (NLP) analysis of over 22 million PubMed abstracts and manual checking of the literature associated with 20 randomly selected dark proteins reinforced the predicted interactions between proteins and pathways. To enhance the visualization and exploration of dark proteins within Reactome pathways, we developed the Reactome IDG portal, deployed at https://idg.reactome.org, a web application featuring tissue-specific protein and gene expression overlay, as well as drug interactions. Our integrated computational approach, together with the user-friendly web platform, offers a valuable resource for uncovering potential biological functions and therapeutic implications of dark proteins.

Using the Reactome Database

Pathway databases provide descriptions of the roles of proteins, nucleic acids, lipids, carbohydrates, and other molecular entities within their biological cellular contexts. Pathway-centric views of these roles may allow for the discovery of unexpected functional relationships in data such as gene expression profiles and somatic mutation catalogues from tumor cells. For this reason, there is a high demand for high-quality pathway databases and their associated tools. The Reactome project (a collaboration between the Ontario Institute for Cancer Research, New York University Langone Health, the European Bioinformatics Institute, and Oregon Health & Science University) is one such pathway database. Reactome collects detailed information on biological pathways and processes in humans from the primary literature. Reactome content is manually curated, expert-authored, and peer-reviewed and spans the gamut from simple intermediate metabolism to signaling pathways and complex cellular events. This information is supplemented with likely orthologous molecular reactions in mouse, rat, zebrafish, worm, and other model organisms. © 2023 The Authors. Current Protocols published by Wiley Periodicals LLC. Basic Protocol 1: Browsing a Reactome pathway Basic Protocol 2: Exploring Reactome annotations of disease and drugs Basic Protocol 3: Finding the pathways involving a gene or protein Alternate Protocol 1: Finding the pathways involving a gene or protein using UniProtKB (SwissProt), Ensembl, or Entrez gene identifier Alternate Protocol 2: Using advanced search Basic Protocol 4: Using the Reactome pathway analysis tool to identify statistically overrepresented pathways Basic Protocol 5: Using the Reactome pathway analysis tool to overlay expression data onto Reactome pathway diagrams Basic Protocol 6: Comparing inferred model organism and human pathways using the Species Comparison tool Basic Protocol 7: Comparing tissue-specific expression using the Tissue Distribution tool.

The potential of a data centred approach & knowledge graph data representation in chemical safety and drug design

Big Data pervades nearly all areas of life sciences, yet the analysis of large integrated data sets remains a major challenge. Moreover, the field of life sciences is highly fragmented and, consequently, so is its data, knowledge, and standards. This, in turn, makes integrated data analysis and knowledge gathering across sub-fields a demanding task. At the same time, the integration of various research angles and data types is crucial for modelling the complexity of organisms and biological processes in a holistic manner. This is especially valid in the context of drug development and chemical safety assessment where computational methods can provide solutions for the urgent need of fast, effective, and sustainable approaches. At the same time, such computational methods require the development of methodologies suitable for an integrated and data centred Big Data view. Here we discuss Knowledge Graphs (KG) as a solution to a data centred analysis approach for drug and chemical development and safety assessment. KGs are knowledge bases, data analysis engines, and knowledge discovery systems all in one, allowing them to be used from simple data retrieval, over meta-analysis to complex predictive and knowledge discovery systems. Therefore, KGs have immense potential to advance the data centred approach, the re-usability, and informativity of data. Furthermore, they can improve the power of analysis, and the complexity of modelled processes, all while providing knowledge in a natively human understandable network data model.

Integrative analysis of drug response and clinical outcome in acute myeloid leukemia

Acute myeloid leukemia (AML) is a cancer of myeloid-lineage cells with limited therapeutic options. We previously combined ex vivo drug sensitivity with genomic, transcriptomic, and clinical annotations for a large cohort of AML patients, which facilitated discovery of functional genomic correlates. Here, we present a dataset that has been harmonized with our initial report to yield a cumulative cohort of 805 patients (942 specimens). We show strong cross-cohort concordance and identify features of drug response. Further, deconvoluting transcriptomic data shows that drug sensitivity is governed broadly by AML cell differentiation state, sometimes conditionally affecting other correlates of response. Finally, modeling of clinical outcome reveals a single gene, PEAR1, to be among the strongest predictors of patient survival, especially for young patients. Collectively, this report expands a large functional genomic resource, offers avenues for mechanistic exploration and drug development, and reveals tools for predicting outcome in AML.

A systems biology approach identifies candidate drugs to reduce mortality in severely ill patients with COVID-19

Despite the availability of highly efficacious vaccines, coronavirus disease 2019 (COVID-19) caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) lacks effective drug treatment, which results in a high rate of mortality. To address this therapeutic shortcoming, we applied a systems biology approach to the study of patients hospitalized with severe COVID. We show that, at the time of hospital admission, patients who were equivalent on the clinical ordinal scale displayed significant differential monocyte epigenetic and transcriptomic attributes between those who would survive and those who would succumb to COVID-19. We identified messenger RNA metabolism, RNA splicing, and interferon signaling pathways as key host responses overactivated by patients who would not survive. Those pathways are prime drug targets to reduce mortality of critically ill patients with COVID-19, leading us to identify tacrolimus, zotatifin, and nintedanib as three strong candidates for treatment of severely ill patients at the time of hospital admission.

Evaluating the predictive accuracy of curated biological pathways in a public knowledgebase

Abstract

Reactome is a database of human biological pathways manually curated from the primary literature and peer-reviewed by experts. To evaluate the utility of Reactome pathways for predicting functional consequences of genetic perturbations, we compared predictions of perturbation effects based on Reactome pathways against published empirical observations. Ten cancer-relevant Reactome pathways, representing diverse biological processes such as signal transduction, cell division, DNA repair and transcriptional regulation, were selected for testing. For each pathway, root input nodes and key pathway outputs were defined. We then used pathway-diagram-derived logic graphs to predict, either by inspection by biocurators or using a novel algorithm MP-BioPath, the effects of bidirectional perturbations (upregulation/activation or downregulation/inhibition) of single root inputs on the status of key outputs. These predictions were then compared to published empirical tests. In total, 4968 test cases were analyzed across 10 pathways, of which 847 were supported by published empirical findings. Out of the 847 test cases, curators’ predictions agreed with the experimental evidence in 670 and disagreed in 177 cases, resulting in ∼81% overall accuracy. MP-BioPath predictions agreed with experimental evidence for 625 and disagreed for 222 test cases, resulting in ∼75% overall accuracy. The expected accuracy of random guessing was 33%. Per-pathway accuracy did not correlate with the number of pathway edges nor the number of pathway nodes but varied across pathways, ranging from 56% (curator)/44% (MP-BioPath) for ‘Mitotic G1 phase and G1/S transition’ to 100% (curator)/94% (MP-BioPath) for ‘RAF/MAP kinase cascade’. This study highlights the potential of pathway databases such as Reactome in modeling genetic perturbations, promoting standardization of experimental pathway activity readout and supporting hypothesis-driven research by revealing relationships between pathway inputs and outputs that have not yet been directly experimentally tested.

Database URL

www.reactome.org

Exploring Toxins for Hunting SARS-CoV-2 Main Protease Inhibitors: Molecular Docking, Molecular Dynamics, Pharmacokinetic Properties, and Reactome Study

The main protease (Mpro) is a potential druggable target in SARS-CoV-2 replication. Herein, an in silico study was conducted to mine for Mpro inhibitors from toxin sources. A toxin and toxin-target database (T3DB) was virtually screened for inhibitor activity towards the Mpro enzyme utilizing molecular docking calculations. Promising toxins were subsequently characterized using a combination of molecular dynamics (MD) simulations and molecular mechanics-generalized Born surface area (MM-GBSA) binding energy estimations. According to the MM-GBSA binding energies over 200 ns MD simulations, three toxins—namely philanthotoxin (T3D2489), azaspiracid (T3D2672), and taziprinone (T3D2378)—demonstrated higher binding affinities against SARS-CoV-2 Mpro than the co-crystalized inhibitor XF7 with MM-GBSA binding energies of −58.9, −55.9, −50.1, and −43.7 kcal/mol, respectively. The molecular network analyses showed that philanthotoxin provides a ligand lead using the STRING database, which includes the biochemical top 20 signaling genes CTSB, CTSL, and CTSK. Ultimately, pathway enrichment analysis (PEA) and Reactome mining results revealed that philanthotoxin could prevent severe lung injury in COVID-19 patients through the remodeling of interleukins (IL-4 and IL-13) and the matrix metalloproteinases (MMPs). These findings have identified that philanthotoxin—a venom of the Egyptian solitary wasp—holds promise as a potential Mpro inhibitor and warrants further in vitro/in vivo validation.

Leveraging biochemical reactions to unravel functional impacts of cancer somatic variants affecting protein interaction interfaces

Background: Considering protein mutations in their biological context is essential for understanding their functional impact, interpretation of high-dimensional datasets and development of effective targeted therapies in personalized medicine. Methods: We combined the curated knowledge of biochemical reactions from Reactome with the analysis of interaction-mediating 3D interfaces from Mechismo. In addition, we provided a software tool for users to explore and browse the analysis results in a multi-scale perspective starting from pathways and reactions to protein-protein interactions and protein 3D structures. Results: We analyzed somatic mutations from TCGA, revealing several significantly impacted reactions and pathways in specific cancer types. We found examples of genes not yet listed as oncodrivers, whose rare mutations were predicted to affect cancer processes similarly to known oncodrivers. Some identified processes lack any known oncodrivers, which suggests potentially new cancer-related processes (e.g. complement cascade reactions). Furthermore, we found that mutations perturbing certain processes are significantly associated with distinct phenotypes (i.e. survival time) in specific cancer types (e.g. PIK3CA centered pathways in LGG and UCEC cancer types), suggesting the translational potential of our approach for patient stratification. Our analysis also uncovered several druggable processes (e.g. GPCR signalling pathways) containing enriched reactions, providing support for new off-label therapeutic options. Conclusions: In summary, we have established a multi-scale approach to study genetic variants based on protein-protein interaction 3D structures. Our approach is different from previously published studies in its focus on biochemical reactions and can be applied to other data types (e.g. post-translational modifications) collected for many types of disease.

Leveraging biochemical reactions to unravel functional impacts of cancer somatic variants affecting protein interaction interfaces

Background: Considering protein mutations in their biological context is essential for understanding their functional impact, interpretation of high-dimensional datasets and development of effective targeted therapies in personalized medicine. Methods: We combined the curated knowledge of biochemical reactions from Reactome with the analysis of interaction-mediating 3D interfaces from Mechismo. In addition, we provided a software tool for users to explore and browse the analysis results in a multi-scale perspective starting from pathways and reactions to protein-protein interactions and protein 3D structures. Results: We analyzed somatic mutations from TCGA, revealing several significantly impacted reactions and pathways in specific cancer types. We found examples of genes not yet listed as oncodrivers, whose rare mutations were predicted to affect cancer processes similarly to known oncodrivers. Some identified processes lack any known oncodrivers, which suggests potentially new cancer-related processes (e.g. complement cascade reactions). Furthermore, we found that mutations perturbing certain processes are significantly associated with distinct phenotypes (i.e. survival time) in specific cancer types (e.g. PIK3CA centered pathways in LGG and UCEC cancer types), suggesting the translational potential of our approach for patient stratification. Our analysis also uncovered several druggable processes (e.g. GPCR signalling pathways) containing enriched reactions, providing support for new off-label therapeutic options. Conclusions: In summary, we have established a multi-scale approach to study genetic variants based on protein-protein interaction 3D structures. Our approach is different from previously published studies in its focus on biochemical reactions and can be applied to other data types (e.g. post-translational modifications) collected for many types of disease.

Leveraging biochemical reactions to unravel functional impacts of cancer somatic variants affecting protein interaction interfaces

Background: Considering protein mutations in their biological context is essential for understanding their functional impact, interpretation of high-dimensional datasets and development of effective targeted therapies in personalized medicine. Methods: We combined the curated knowledge of biochemical reactions from Reactome with the analysis of interaction-mediating 3D interfaces from Mechismo. In addition, we provided a software tool for users to explore and browse the analysis results in a multi-scale perspective starting from pathways and reactions to protein-protein interactions and protein 3D structures. Results: We analyzed somatic mutations from TCGA, revealing several significantly impacted reactions and pathways in specific cancer types. We found examples of genes not yet listed as oncodrivers, whose rare mutations were predicted to affect cancer processes similarly to known oncodrivers. Some identified processes lack any known oncodrivers, which suggests potentially new cancer-related processes (e.g. complement cascade reactions). Furthermore, we found that mutations perturbing certain processes are significantly associated with distinct phenotypes (i.e. survival time) in specific cancer types (e.g. PIK3CA centered pathways in LGG and UCEC cancer types), suggesting the translational potential of our approach for patient stratification. Our analysis also uncovered several druggable processes (e.g. GPCR signalling pathways) containing enriched reactions, providing support for new off-label therapeutic options. Conclusions: In summary, we have established a multi-scale approach to study genetic variants based on protein-protein interaction 3D structures. Our approach is different from previously published studies in its focus on biochemical reactions and can be applied to other data types (e.g. post-translational modifications) collected for many types of disease.

Integrated Genomic Profiling and Drug Screening of Patient-Derived Cultures Identifies Individualized Copy Number-Dependent Susceptibilities Involving PI3K Pathway and 17q Genes in Neuroblastoma

Neuroblastoma is the commonest extracranial pediatric malignancy. With few recurrent single nucleotide variations (SNVs), mutation-based precision oncology approaches have limited utility, but its frequent and heterogenous copy number variations (CNVs) could represent genomic dependencies that may be exploited for personalized therapy. Patient-derived cell culture (PDC) models can facilitate rapid testing of multiple agents to determine such individualized drug-responses. Thus, to study the relationship between individual genomic aberrations and therapeutic susceptibilities, we integrated comprehensive genomic profiling of neuroblastoma tumors with drug screening of corresponding PDCs against 418 targeted inhibitors. We quantified the strength of association between copy number and cytotoxicity, and validated significantly correlated gene-drug pairs in public data and using machine learning models. Somatic mutations were infrequent (3.1 per case), but copy number losses in 1p (31%) and 11q (38%), and gains in 17q (69%) were prevalent. Critically, in-vitro cytotoxicity significantly correlated only with CNVs, but not SNVs. Among 1278 significantly correlated gene-drug pairs, copy number of GNA13 and DNA damage response genes CBL, DNMT3A, and PPM1D were most significantly correlated with cytotoxicity; the drugs most commonly associated with these genes were PI3K/mTOR inhibitor PIK-75, and CDK inhibitors P276-00, SNS-032, AT7519, flavopiridol and dinaciclib. Predictive Markov random field models constructed from CNVs alone recapitulated the true z-score-weighted associations, with the strongest gene-drug functional interactions in subnetworks involving PI3K and JAK-STAT pathways. Together, our data defined individualized dose-dependent relationships between copy number gains of PI3K and STAT family genes particularly on 17q and susceptibility to PI3K and cell cycle agents in neuroblastoma. Integration of genomic profiling and drug screening of patient-derived models of neuroblastoma can quantitatively define copy number-dependent sensitivities to targeted inhibitors, which can guide personalized therapy for such mutationally quiet cancers.

Gastroprotection against Rat Ulcers by Nephthea Sterol Derivative

Different species belonging to the genus Nephthea (Acyonaceae) are a rich resource for bioactive secondary metabolites. The literature reveals that the gastroprotective effects of marine secondary metabolites have not been comprehensively studied in vivo. Hence, the present investigation aimed to examine and determine the anti-ulcer activity of 4α,24-dimethyl-5α-cholest-8β,18-dihydroxy,22E-en-3β-ol (ST-1) isolated from samples of a Nephthea species. This in vivo study was supported by in silico molecular docking and protein-protein interaction techniques. Oral administration of ST-1 reduced rat stomach ulcers with a concurrent increase in gastric mucosa. Molecular docking calculations against the H+/K+-ATPase transporter showed a higher binding affinity of ST-1, with a docking score value of -9.9 kcal/mol and a pKi value of 59.7 nM, compared to ranitidine (a commercial proton pump inhibitor, which gave values of -6.2 kcal/mol and 27.9 µM, respectively). The combined PEA-reactome analysis results revealed promising evidence of ST-1 potency as an anti-ulcer compound through significant modulation of the gene set controlling the PI3K signaling pathway, which subsequently plays a crucial role in signaling regarding epithelialization and tissue regeneration, tissue repairing and tissue remodeling. These results indicate a probable protective role for ST-1 against ethanol-induced gastric ulcers.

Blue Biotechnology: Computational Screening of Sarcophyton Cembranoid Diterpenes for SARS-CoV-2 Main Protease Inhibition

The coronavirus pandemic has affected more than 150 million people, while over 3.25 million people have died from the coronavirus disease 2019 (COVID-19). As there are no established therapies for COVID-19 treatment, drugs that inhibit viral replication are a promising target; specifically, the main protease (Mpro) that process CoV-encoded polyproteins serves as an Achilles heel for assembly of replication-transcription machinery as well as down-stream viral replication. In the search for potential antiviral drugs that target Mpro, a series of cembranoid diterpenes from the biologically active soft-coral genus Sarcophyton have been examined as SARS-CoV-2 Mpro inhibitors. Over 360 metabolites from the genus were screened using molecular docking calculations. Promising diterpenes were further characterized by molecular dynamics (MD) simulations based on molecular mechanics-generalized Born surface area (MM-GBSA) binding energy calculations. According to in silico calculations, five cembranoid diterpenes manifested adequate binding affinities as Mpro inhibitors with ΔGbinding < -33.0 kcal/mol. Binding energy and structural analyses of the most potent Sarcophyton inhibitor, bislatumlide A (340), was compared to darunavir, an HIV protease inhibitor that has been recently subjected to clinical-trial as an anti-COVID-19 drug. In silico analysis indicates that 340 has a higher binding affinity against Mpro than darunavir with ΔGbinding values of -43.8 and -34.8 kcal/mol, respectively throughout 100 ns MD simulations. Drug-likeness calculations revealed robust bioavailability and protein-protein interactions were identified for 340; biochemical signaling genes included ACE, MAPK14 and ESR1 as identified based on a STRING database. Pathway enrichment analysis combined with reactome mining revealed that 340 has the capability to re-modulate the p38 MAPK pathway hijacked by SARS-CoV-2 and antagonize injurious effects. These findings justify further in vivo and in vitro testing of 340 as an antiviral agent against SARS-CoV-2.

Two novel oxetane containing lignans and a new megastigmane from Paronychia arabica and in silico analysis of them as prospective SARS-CoV-2 inhibitors

The chemical characterization of the extract of the aerial parts of Paronychia arabica afforded two oxetane containing lignans, paronychiarabicine A (1) and B (2), and one new megastigmane, paronychiarabicastigmane A (3), alongside a known lignan (4), eight known phenolic compounds (5–12), one known elemene sesquiterpene (13) and one steroid glycoside (14). The chemical structures of the isolated compounds were constructed based upon the HRMS, 1D, and 2D-NMR results. The absolute configurations were established via NOESY experiments as well as experimental and TDDFT-calculated electronic circular dichroism (ECD). Utilizing molecular docking, the binding scores and modes of compounds 1–3 towards the SARS-CoV-2 main protease (M^pro), papain-like protease (PL^pro), and RNA-dependent RNA polymerase (RdRp) were revealed. Compound 3 exhibited a promising docking score (−9.8 kcal mol⁻¹) against SARS-CoV-2 M^pro by forming seven hydrogen bonds inside the active site with the key amino acids. The reactome pathway enrichment analysis revealed a correlation between the inhibition of GSK3 and GSK3B genes (identified as the main targets of megastigmane treatment) and significant inhibition of SARS-CoV-1 viral replication in infected Vero E6 cells. Our results manifest a novel understanding of genes, proteins and corresponding pathways against SARS-CoV-2 infection and could facilitate the identification and characterization of novel therapeutic targets as treatments of SARS-CoV-2 infection.

The hydromethanolic extract of Paronychia arabica aerial parts afforded two oxetane containing lignans, paronychiarabicine A (1) and B (2), and one new megastigmane, paronychiarabicastigmane A (3), alongside a known secondary metabolites (4–14).

In Silico Mining of Terpenes from Red-Sea Invertebrates for SARS-CoV-2 Main Protease (Mpro) Inhibitors

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the causative agent for the COVID-19 pandemic, which generated more than 1.82 million deaths in 2020 alone, in addition to 83.8 million infections. Currently, there is no antiviral medication to treat COVID-19. In the search for drug leads, marine-derived metabolites are reported here as prospective SARS-CoV-2 inhibitors. Two hundred and twenty-seven terpene natural products isolated from the biodiverse Red-Sea ecosystem were screened for inhibitor activity against the SARS-CoV-2 main protease (Mpro) using molecular docking and molecular dynamics (MD) simulations combined with molecular mechanics/generalized Born surface area binding energy calculations. On the basis of in silico analyses, six terpenes demonstrated high potency as Mpro inhibitors with ΔGbinding ≤ -40.0 kcal/mol. The stability and binding affinity of the most potent metabolite, erylosides B, were compared to the human immunodeficiency virus protease inhibitor, lopinavir. Erylosides B showed greater binding affinity towards SARS-CoV-2 Mpro than lopinavir over 100 ns with ΔGbinding values of -51.9 vs. -33.6 kcal/mol, respectively. Protein-protein interactions indicate that erylosides B biochemical signaling shares gene components that mediate severe acute respiratory syndrome diseases, including the cytokine- and immune-signaling components BCL2L1, IL2, and PRKC. Pathway enrichment analysis and Boolean network modeling were performed towards a deep dissection and mining of the erylosides B target-function interactions. The current study identifies erylosides B as a promising anti-COVID-19 drug lead that warrants further in vitro and in vivo testing.

Biological characteristics of aging in human acute myeloid leukemia cells: the possible importance of aldehyde dehydrogenase, the cytoskeleton and altered transcriptional regulation

Patients with acute myeloid leukemia (AML) have a median age of 65-70 years at diagnosis. Elderly patients have more chemoresistant disease, and this is partly due to decreased frequencies of favorable and increased frequencies of adverse genetic abnormalities. However, aging-dependent differences may also contribute. We therefore compared AML cell proteomic and phosphoproteomic profiles for (i) elderly low-risk and younger low-risk patients with favorable genetic abnormalities; and (ii) high-risk patients with adverse genetic abnormalities and a higher median age against all low-risk patients with lower median age. Elderly low-risk and younger low-risk patients showed mainly phosphoproteomic differences especially involving transcriptional regulators and cytoskeleton. When comparing high-risk and low-risk patients both proteomic and phosphoproteomic studies showed differences involving cytoskeleton and immunoregulation but also transcriptional regulation and cell division. The age-associated prognostic impact of cyclin-dependent kinases was dependent on the cellular context. The protein level of the adverse prognostic biomarker mitochondrial aldehyde dehydrogenase (ALDH2) showed a similar significant upregulation both in elderly low-risk and elderly high-risk patients. Our results suggest that molecular mechanisms associated with cellular aging influence chemoresistance of AML cells, and especially the cytoskeleton function may then influence cellular hallmarks of aging, e.g. mitosis, polarity, intracellular transport and adhesion.

Illuminating biological pathways for drug targeting in head and neck squamous cell carcinoma

Head and neck squamous cell carcinoma (HNSCC) remains a morbid disease with poor prognosis and treatment that typically leaves patients with permanent damage to critical functions such as eating and talking. Currently only three targeted therapies are FDA approved for use in HNSCC, two of which are recently approved immunotherapies. In this work, we identify biological pathways involved with this disease that could potentially be targeted by current FDA approved cancer drugs and thereby expand the pool of potential therapies for use in HNSCC treatment. We analyzed 508 HNSCC patients with sequencing information from the Genomic Data Commons (GDC) database and assessed which biological pathways were significantly enriched for somatic mutations or copy number alterations. We then further classified pathways as either “light” or “dark” to the current reach of FDA-approved cancer drugs using the Cancer Targetome, a compendium of drug-target information. Light pathways are statistically enriched with somatic mutations (or copy number alterations) and contain one or more targets of current FDA-approved cancer drugs, while dark pathways are enriched with somatic mutations (or copy number alterations) but not currently targeted by FDA-approved cancer drugs. Our analyses indicated that approximately 35–38% of disease-specific pathways are in scope for repurposing of current cancer drugs. We further assess light and dark pathways for subgroups of patient tumor samples according to HPV status. The framework of light and dark pathways for HNSCC-enriched biological pathways allows us to better prioritize targeted therapies for further research in HNSCC based on the HNSCC genetic landscape and FDA-approved cancer drug information. We also highlight the importance in the identification of sub-pathways where targeting and cross targeting of other pathways may be most beneficial to predict positive or negative synergy with potential clinical significance. This framework is ideal for precision drug panel development, as well as identification of highly aberrant, untargeted candidates for future drug development.