Network-based prediction of drug combinations

Predicting anti-cancer drug combination responses with a temporal cell state network model

Cancer chemotherapy combines multiple drugs, but predicting the effects of drug combinations on cancer cell proliferation remains challenging, even for simple in vitro systems. We hypothesized that by combining knowledge of single drug dose responses and cell state transition network dynamics, we could predict how a population of cancer cells will respond to drug combinations. We tested this hypothesis here using three targeted inhibitors of different cell cycle states in two different cell lines in vitro. We formulated a Markov model to capture temporal cell state transitions between different cell cycle phases, with single drug data constraining how drug doses affect transition rates. This model was able to predict the landscape of all three different pairwise drug combinations across all dose ranges for both cell lines with no additional data. While further application to different cell lines, more drugs, additional cell state networks, and more complex co-culture or in vivo systems remain, this work demonstrates how currently available or attainable information could be sufficient for prediction of drug combination response for single cell lines in vitro.

A Science-Based Methodology Framework for the Assessment of Combination Safety Risks in Clinical Trials

Multiple components factor into the assessment of combination safety risks when two or more novel individual products are used in combination in clinical trials. These include, but are not limited to, biology, biochemistry, pharmacology, class effects, and preclinical and clinical findings (such as adverse drug reactions, drug target and mechanism of action, target expression, signaling, and drug-drug interactions). This paper presents a science-based methodology framework for the assessment of combination safety risks when two or more investigational products are used in clinical trials. The aim of this methodology framework is to improve prediction of the risks, to enable the appropriate safety risk mitigation and management to be put in place for the combination, and the development of the project combination safety strategy.

CancerGPT: Few-shot Drug Pair Synergy Prediction using Large Pre-trained Language Models

Large pre-trained language models (LLMs) have been shown to have significant potential in few-shot learning across various fields, even with minimal training data. However, their ability to generalize to unseen tasks in more complex fields, such as biology, has yet to be fully evaluated. LLMs can offer a promising alternative approach for biological inference, particularly in cases where structured data and sample size are limited, by extracting prior knowledge from text corpora. Our proposed few-shot learning approach uses LLMs to predict the synergy of drug pairs in rare tissues that lack structured data and features. Our experiments, which involved seven rare tissues from different cancer types, demonstrated that the LLM-based prediction model achieved significant accuracy with very few or zero samples. Our proposed model, the CancerGPT (with ~ 124M parameters), was even comparable to the larger fine-tuned GPT-3 model (with ~ 175B parameters). Our research is the first to tackle drug pair synergy prediction in rare tissues with limited data. We are also the first to utilize an LLM-based prediction model for biological reaction prediction tasks.

Study on the Metabolic Transformation Rule of Enrofloxacin Combined with Tilmicosin in Laying Hens

There is often abuse of drugs in livestock and poultry production, and the improper use of drugs leads to the existence of a low level of residues in eggs, which is a potential threat to human safety. Enrofloxacin (EF) and tilmicosin (TIM) are regularly combined for the prevention and treatment of poultry diseases. The current studies on EF or TIM mainly focus on a single drug, and the effects of the combined application of these two antibiotics on EF metabolism in laying hens are rarely reported. In this study, liquid chromatography-tandem mass spectrometry (LC-MS/MS) was used to determine the residual EF and TIM in laying hens and to investigate the effect of TIM on the EF metabolism in laying hens. In this paper, we first establish a method that can detect EF and TIM simultaneously. Secondly, the results showed that the highest concentration of EF in the egg samples was 974.92 ± 441.71 μg/kg on the 5th day of treatment. The highest concentration of EF in the egg samples of the combined administration group was 1256.41 ± 226.10 μg/kg on the 5th day of administration. The results showed that when EF and TIM were used in combination, the residue of EF in the eggs was increased, the elimination rate of EF was decreased, and the half-life of EF was increased. Therefore, the use of EF and TIM in combination should be treated with greater care and supervision should be strengthened to avoid risks to human health.

Lacking mechanistic disease definitions and corresponding association data hamper progress in network medicine and beyond

A long-term objective of network medicine is to replace our current, mainly phenotype-based disease definitions by subtypes of health conditions corresponding to distinct pathomechanisms. For this, molecular and health data are modeled as networks and are mined for pathomechanisms. However, many such studies rely on large-scale disease association data where diseases are annotated using the very phenotype-based disease definitions the network medicine field aims to overcome. This raises the question to which extent the biases mechanistically inadequate disease annotations introduce in disease association data distort the results of studies which use such data for pathomechanism mining. We address this question using global- and local-scale analyses of networks constructed from disease association data of various types. Our results indicate that large-scale disease association data should be used with care for pathomechanism mining and that analyses of such data should be accompanied by close-up analyses of molecular data for well-characterized patient cohorts.

SWEET: a single-sample network inference method for deciphering individual features in disease

Recently, extracting inherent biological system information (e.g. cellular networks) from genome-wide expression profiles for developing personalized diagnostic and therapeutic strategies has become increasingly important. However, accurately constructing single-sample networks (SINs) to capture individual characteristics and heterogeneity in disease remains challenging. Here, we propose a sample-specific-weighted correlation network (SWEET) method to model SINs by integrating the genome-wide sample-to-sample correlation (i.e. sample weights) with the differential network between perturbed and aggregate networks. For a group of samples, the genome-wide sample weights can be assessed without prior knowledge of intrinsic subpopulations to address the network edge number bias caused by sample size differences. Compared with the state-of-the-art SIN inference methods, the SWEET SINs in 16 cancers more likely fit the scale-free property, display higher overlap with the human interactomes and perform better in identifying three types of cancer-related genes. Moreover, integrating SWEET SINs with a network proximity measure facilitates characterizing individual features and therapy in diseases, such as somatic mutation, mut-driver and essential genes. Biological experiments further validated two candidate repurposable drugs, albendazole for head and neck squamous cell carcinoma (HNSCC) and lung adenocarcinoma (LUAD) and encorafenib for HNSCC. By applying SWEET, we also identified two possible LUAD subtypes that exhibit distinct clinical features and molecular mechanisms. Overall, the SWEET method complements current SIN inference and analysis methods and presents a view of biological systems at the network level to offer numerous clues for further investigation and clinical translation in network medicine and precision medicine.

Combined Multiomics and In Silico Approach Uncovers PRKAR1A as a Putative Therapeutic Target in Multi-Organ Dysfunction Syndrome

Despite all epidemiological, clinical, and experimental research efforts, therapeutic concepts in sepsis and sepsis-induced multi-organ dysfunction syndrome (MODS) remain limited and unsatisfactory. Currently, gene expression data sets are widely utilized to discover new biomarkers and therapeutic targets in diseases. In the present study, we analyzed MODS expression profiles (comprising 13 sepsis and 8 control samples) retrieved from NCBI-GEO and found 359 differentially expressed genes (DEGs), among which 170 were downregulated and 189 were upregulated. Next, we employed the weighted gene co-expression network analysis (WGCNA) to establish a MODS-associated gene co-expression network (weighted) and identified representative module genes having an elevated correlation with age. Based on the results, a turquoise module was picked as our hub module. Further, we constructed the PPI network comprising 35 hub module DEGs. The DEGs involved in the highest-confidence PPI network were utilized for collecting pathway and gene ontology (GO) terms using various libraries. Nucleotide di- and triphosphate biosynthesis and interconversion was the most significant pathway. Also, 3 DEGs within our PPI network were involved in the top 5 significantly enriched ontology terms, with hypercortisolism being the most significant term. PRKAR1A was the overlapping gene between top 5 significant pathways and GO terms, respectively. PRKAR1A was considered as a therapeutic target in MODS, and 2992 ligands were screened for binding with PRKAR1A. Among these ligands, 3 molecules based on CDOCKER score (molecular dynamics simulated-based score, which allows us to rank the binding poses according to their quality and to identify the best pose for each system) and crucial interaction with human PRKAR1A coding protein and protein kinase-cyclic nucleotide binding domains (PKA RI alpha CNB-B domain) via active site binding residues, viz. Val283, Val302, Gln304, Val315, Ile327, Ala336, Ala337, Val339, Tyr373, and Asn374, were considered as lead molecules.

Lower Concentrations of Amphotericin B Combined with Ent-Hardwickiic Acid Are Effective against Candida Strains

Life-threatening Candida infections have increased with the COVID-19 pandemic, and the already limited arsenal of antifungal drugs has become even more restricted due to its side effects associated with complications after SARS-CoV-2 infection. Drug combination strategies have the potential to reduce the risk of side effects without loss of therapeutic efficacy. The aim of this study was to evaluate the combination of ent-hardwickiic acid with low concentrations of amphotericin B against Candida strains. The minimum inhibitory concentration (MIC) values were determined for amphotericin B and ent-hardwickiic acid as isolated compounds and for 77 combinations of amphotericin B and ent-hardwickiic acid concentrations that were assessed by using the checkerboard microdilution method. Time–kill assays were performed in order to assess the fungistatic or fungicidal nature of the different combinations. The strategy of combining both compounds markedly reduced the MIC values from 16 µg/mL to 1 µg/mL of amphotericin B and from 12.5 µg/mL to 6.25 µg/mL of ent-hardwickiic acid, from isolated to combined, against C. albicans resistant to azoles. The combination of 1 µg/mL of amphotericin B with 6.25 µg/mL of ent-hardwickiic acid killed all the cells of the same strain within four hours of incubation.

Network-informed discovery of multidrug combinations for ERα+/HER2-/PI3Kα-mutant breast cancer

Breast cancer is a persistent threat to women worldwide. A large proportion of breast cancers are dependent on the estrogen receptor α (ERα) for tumor progression. Therefore, targeting ERα with antagonists, such as tamoxifen, or estrogen deprivation by aromatase inhibitors remain standard therapies for ERα + breast cancer. The clinical benefits of monotherapy are often counterbalanced by off-target toxicity and development of resistance. Combinations of more than two drugs might be of great therapeutic value to prevent resistance, and to reduce doses, and hence, decrease toxicity. We mined data from the literature and public repositories to construct a network of potential drug targets for synergistic multidrug combinations. With 9 drugs, we performed a phenotypic combinatorial screen with ERα + breast cancer cell lines. We identified two optimized low-dose combinations of 3 and 4 drugs of high therapeutic relevance to the frequent ERα + /HER2-/PI3Kα-mutant subtype of breast cancer. The 3-drug combination targets ERα in combination with PI3Kα and cyclin-dependent kinase inhibitor 1 (p21). In addition, the 4-drug combination contains an inhibitor for poly (ADP-ribose) polymerase 1 (PARP1), which showed benefits in long-term treatments. Moreover, we validated the efficacy of the combinations in tamoxifen-resistant cell lines, patient-derived organoids, and xenograft experiments. Thus, we propose multidrug combinations that have the potential to overcome the standard issues of current monotherapies.

A hybrid deep forest-based method for predicting synergistic drug combinations

Combination therapy is a promising approach in treating multiple complex diseases. However, the large search space of available drug combinations exacerbates challenge for experimental screening. To predict synergistic drug combinations in different cancer cell lines, we propose an improved deep forest-based method, ForSyn, and design two forest types embedded in ForSyn. ForSyn handles imbalanced and high-dimensional data in medium-/small-scale datasets, which are inherent characteristics of drug combination datasets. Compared with 12 state-of-the-art methods, ForSyn ranks first on four metrics for eight datasets with different feature combinations. We conduct a systematic analysis to identify the most appropriate configuration parameters. We validate the predictive value of ForSyn with cell-based experiments on several previously unexplored drug combinations. Finally, a systematic analysis of feature importance is performed on the top contributing features extracted by ForSyn. The resulting key genes may play key roles on corresponding cancers.

Discovery and Validation of Traditional Chinese and Western Medicine Combination Antirheumatoid Arthritis Drugs Based on Machine Learning (Random Forest Model)

The combination of traditional Chinese medicine (TCM) and Western medicine is a promising method for treating rheumatoid arthritis (RA). Combining the two fully exploits the advantages of Western and TCM to treat RA and has the potential to greatly improve the therapeutic effect on RA. In this study, we developed a combination drug training set by using 16 characteristic variables based on the characteristics of small molecules of TCM ingredients and Food and Drug Administration-certified combination drug data downloaded from the DrugCombDB database. Furthermore, we compared the prediction and classification abilities of five models: the k-nearest neighbors, naive Bayes, support vector machine, random forest, and AdaBoost algorithms. The random forest model was selected as the classification and prediction model for Western and TCM and Western combination drugs. We collected data for 41 small molecules of TCM ingredients from the Traditional Chinese Medicine Systems Pharmacology database and 10 small molecule drugs commonly used in anti-RA treatment from the DrugBank database. Combinations of Western and TCM for anti-RA treatment were screened. Finally, the CellTiter-Glo method was used to determine the synergy of these combinations, and the 15 most predicted drug combinations were carried out experimental verification. Myricetin, rhein, nobiletin, and fisetin had high synergy with celecoxib, and rhein had high synergy with hydroxychloroquine. The preliminary findings of this study can be further applied for practical clinical anti-RA combined treatment strategies and serve as a reference for clinical treatment of RA with integrated Western and TCM.

Gene Regulatory Network-Classifier: Gene Regulatory Network-Based Classifier and Its Applications to Gastric Cancer Drug (5-Fluorouracil) Marker Identification

The complex mechanisms of diseases involve the disturbance of the molecular network, rather than disorder in a single gene, implying that single gene-based analysis is insufficient to understand these mechanisms. Gene regulatory networks (GRNs) have attracted a lot of interest and various approaches have been developed for their statistical inference and gene network-based analysis. Although various computational methods have been developed, relatively little attention has been paid to incorporation of biological knowledge into the computational approaches. Furthermore, existing studies on network-based analysis perform prediction/classification of status of cell lines based on preconstructed GRNs, implying that we cannot extract prediction/classification-specific gene networks, leading to difficulty in interpretation of biological mechanisms and marker identification related to the status of cancer cell lines. We developed a novel strategy to build a GRN-based classifier, called a GRN-classifier. The proposed GRN-classifier estimates GRNs and classifies cell lines simultaneously, where the gene network is estimated to minimize error in gene network estimation and the negative log-likelihood for classifying cell lines. Thus, we can identify biological status-specific gene regulatory systems, enabling us to achieve biologically reliable interpretation of the classification. We also propose an algorithm to implement the GRN-classifier based on coordinate descent update. Monte Carlo simulations were conducted to examine performance of the GRN-classifier. Results: Our strategy provides effective results in feature selection in the classification model and edge selection in gene network estimation. The GRN-classifier also shows outstanding classification accuracy. We apply the GRN-classifier to classify cancer cell lines into anticancer drug-related status, that is, 5-fluorouracil (5-FU)-sensitive/resistant and 5-FU target/nontarget cancer cell lines. We then identified 5-FU markers based on 5-FU-related status classification-specific gene networks. The mechanisms of the identified markers were verified through literature survey. Our results suggest that the molecular interplay between MYOF and AHNAK2 may play a crucial role in drug resistance and can provide information on the chemotherapy efficiency of 5-FU. It is also suggested that suppression of the identified 5-FU markers, including MYOF/AHNAK2 and AKR1C1/AKR1C3 may improve 5-FU resistance of cancer cell lines.

Combination of Antidepressants and Antipsychotics as A Novel Treatment Option for Psychosis in Alzheimer's Disease

Background

Psychotic symptoms are reported as one of the most common complications of Alzheimer's disease (AD), affecting approximately half of AD patients, in whom they are associated with more rapid deterioration and increased mortality. Empiric treatments, namely first and second-generation antipsychotics, confer modest efficacy in AD patients with psychosis (AD+P) and themselves increase mortality. A recent genome-wide meta-analysis and early clinical trials suggest the use and beneficial effects of antidepressants among AD+P patients. This motivates our rationale for exploring their potential as a novel combination therapy option amongst these patients.

Methods

We included University of Pittsburgh Medical Center (UPMC) electronic medical records (EMRs) of 10,260 AD patients from January 2004 and October 2019 in our study. Survival analysis was performed to assess the effects of the combination of antipsychotics and antidepressants on the mortality of these patients. To provide more valuable insights on the hidden mechanisms of the combinatorial therapy, a protein-protein interaction (PPI) network representing AD+P was built, and network analysis methods were used to quantify the efficacy of these drugs on AD+P. An indicator score combining the measurements on the separation between drugs and the proximity between the drugs and AD+P was used to measure the effect of an antipsychotic-antidepressant drug pair against AD+P.

Results

Our survival analyses replicated that antipsychotic usage is strongly associated with increased mortality in AD patients while the co-administration of antidepressants with antipsychotics showed a significant beneficial effect in reducing mortality. Our network analysis showed that the targets of antipsychotics and antidepressants are well-separated, and antipsychotics and antidepressants have similar proximity scores to AD+P. Eight drug pairs, including some popular recommendations like Aripiprazole/Sertraline and other pairs not reported previously like Iloperidone/Maprotiline showed higher than average indicator scores which suggest their potential in treating AD+P via strong synergetic effects as seen in our study.

Conclusion

Our proposed combinations of antipsychotics and antidepressants therapy showed a strong superiority over current antipsychotics treatment for AD+P. The observed beneficial effects can be further strengthened by optimizing drug-pair selection based on our systems pharmacology analysis.

Drug repositioning via host-pathogen protein-protein interactions for the treatment of cervical cancer

Introduction

Integrating interaction data with biological knowledge can be a critical approach for drug development or drug repurposing. In this context, host-pathogen-protein-protein interaction (HP-PPI) networks are useful instrument to uncover the phenomena underlying therapeutic effects in infectious diseases, including cervical cancer, which is almost exclusively due to human papillomavirus (HPV) infections. Cervical cancer is one of the second leading causes of death, and HPV16 and HPV18 are the most common subtypes worldwide. Given the limitations of traditionally used virus-directed drug therapies for infectious diseases and, at the same time, recent cancer statistics for cervical cancer cases, the need for innovative treatments becomes clear.

Methods

Accordingly, in this study, we emphasize the potential of host proteins as drug targets and identify promising host protein candidates for cervical cancer by considering potential differences between HPV subtypes (i.e., HPV16 and HPV18) within a novel bioinformatics framework that we have developed. Subsequently, subtype-specific HP-PPI networks were constructed to obtain host proteins. Using this framework, we next selected biologically significant host proteins. Using these prominent host proteins, we performed drug repurposing analysis. Finally, by following our framework we identify the most promising host-oriented drug candidates for cervical cancer.

Results

As a result of this framework, we discovered both previously associated and novel drug candidates, including interferon alfacon-1, pimecrolimus, and hyaluronan specifically for HPV16 and HPV18 subtypes, respectively.

Discussion

Consequently, with this study, we have provided valuable data for further experimental and clinical efforts and presented a novel bioinformatics framework that can be applied to any infectious disease.

Synthesize heterogeneous biological knowledge via representation learning for Alzheimer's disease drug repurposing

Developing drugs for treating Alzheimer's disease has been extremely challenging and costly due to limited knowledge of underlying mechanisms and therapeutic targets. To address the challenge in AD drug development, we developed a multi-task deep learning pipeline that learns biological interactions and AD risk genes, then utilizes multi-level evidence on drug efficacy to identify repurposable drug candidates. Using the embedding derived from the model, we ranked drug candidates based on evidence from post-treatment transcriptomic patterns, efficacy in preclinical models, population-based treatment effects, and clinical trials. We mechanistically validated the top-ranked candidates in neuronal cells, identifying drug combinations with efficacy in reducing oxidative stress and safety in maintaining neuronal viability and morphology. Our neuronal response experiments confirmed several biologically efficacious drug combinations. This pipeline showed that harmonizing heterogeneous and complementary data/knowledge, including human interactome, transcriptome patterns, experimental efficacy, and real-world patient data shed light on the drug development of complex diseases.

SNRMPACDC: computational model focused on Siamese network and random matrix projection for anticancer synergistic drug combination prediction

Synergistic drug combinations can improve the therapeutic effect and reduce the drug dosage to avoid toxicity. In previous years, an in vitro approach was utilized to screen synergistic drug combinations. However, the in vitro method is time-consuming and expensive. With the rapid growth of high-throughput data, computational methods are becoming efficient tools to predict potential synergistic drug combinations. Considering the limitations of the previous computational methods, we developed a new model named Siamese Network and Random Matrix Projection for AntiCancer Drug Combination prediction (SNRMPACDC). Firstly, the Siamese convolutional network and random matrix projection were used to process the features of the two drugs into drug combination features. Then, the features of the cancer cell line were processed through the convolutional network. Finally, the processed features were integrated and input into the multi-layer perceptron network to get the predicted score. Compared with the traditional method of splicing drug features into drug combination features, SNRMPACDC improved the interpretability of drug combination features to a certain extent. In addition, the introduction of convolutional networks can better extract the potential information in the features. SNRMPACDC achieved the root mean-squared error of 15.01 and the Pearson correlation coefficient of 0.75 in 5-fold cross-validation of regression prediction for response data. In addition, SNRMPACDC achieved the AUC of 0.91 ± 0.03 and the AUPR of 0.62 ± 0.05 in 5-fold cross-validation of classification prediction of synergistic or not. These results are almost better than all the previous models. SNRMPACDC would be an effective approach to infer potential anticancer synergistic drug combinations.

A social theory-enhanced graph representation learning framework for multitask prediction of drug-drug interactions

Current machine learning-based methods have achieved inspiring predictions in the scenarios of mono-type and multi-type drug-drug interactions (DDIs), but they all ignore enhancive and depressive pharmacological changes triggered by DDIs. In addition, these pharmacological changes are asymmetric since the roles of two drugs in an interaction are different. More importantly, these pharmacological changes imply significant topological patterns among DDIs. To address the above issues, we first leverage Balance theory and Status theory in social networks to reveal the topological patterns among directed pharmacological DDIs, which are modeled as a signed and directed network. Then, we design a novel graph representation learning model named SGRL-DDI (social theory-enhanced graph representation learning for DDI) to realize the multitask prediction of DDIs. SGRL-DDI model can capture the task-joint information by integrating relation graph convolutional networks with Balance and Status patterns. Moreover, we utilize task-specific deep neural networks to perform two tasks, including the prediction of enhancive/depressive DDIs and the prediction of directed DDIs. Based on DDI entries collected from DrugBank, the superiority of our model is demonstrated by the comparison with other state-of-the-art methods. Furthermore, the ablation study verifies that Balance and Status patterns help characterize directed pharmacological DDIs, and that the joint of two tasks provides better DDI representations than individual tasks. Last, we demonstrate the practical effectiveness of our model by a version-dependent test, where 88.47 and 81.38% DDI out of newly added entries provided by the latest release of DrugBank are validated in two predicting tasks respectively.

CCSynergy: an integrative deep-learning framework enabling context-aware prediction of anti-cancer drug synergy

Combination therapy is a promising strategy for confronting the complexity of cancer. However, experimental exploration of the vast space of potential drug combinations is costly and unfeasible. Therefore, computational methods for predicting drug synergy are much needed for narrowing down this space, especially when examining new cellular contexts. Here, we thus introduce CCSynergy, a flexible, context aware and integrative deep-learning framework that we have established to unleash the potential of the Chemical Checker extended drug bioactivity profiles for the purpose of drug synergy prediction. We have shown that CCSynergy enables predictions of superior accuracy, remarkable robustness and improved context generalizability as compared to the state-of-the-art methods in the field. Having established the potential of CCSynergy for generating experimentally validated predictions, we next exhaustively explored the untested drug combination space. This resulted in a compendium of potentially synergistic drug combinations on hundreds of cancer cell lines, which can guide future experimental screens.

Optimizing Genomic Control in Hit Network-Target Set Model Associations with Lung Adenocarcinoma

Background: Hit network-target sets (HNSs), compiled sets of different network nodes of the same type, are available and play a significant role in cancer development but are notoriously more difficult to select than a single target. This is due to a combination of challenges attributed to the differential of node interactions, node heterogeneity, and the limitations of node-hit information. Methods: In this study, we constructed a lung adenocarcinoma regulatory network using TCGA data and obtained different HNSs of driver nodes (DNs), core modules (CMs) and core nodes (CNs) through three kinds of methods. Then, the optimized HNS (OHNS) was obtained by integrating CMs, CNs and DNs, and the performance of different HNSs was evaluated according to network structure importance, control capability, and clinical value. Results: We found that the OHNS has two main advantages, the central location of the network and the ability to control the network, and it plays an important role in the disease network through its multifaceted capabilities. Three unique pathways were discovered in the OHNS, which is consistent with previous experiments. Additionally, 13 genes were predicted to play roles in risk prognosis, disease drivers, and cell perturbation effects of lung adenocarcinoma, of which 12 may be candidates for new drugs and biomarkers of lung adenocarcinoma. Conclusion: This study can help us understand and control a network more effectively to determine the development trend of a disease, design effective multitarget drugs, and guide the therapeutic community to optimize appropriate strategies according to different research aims in cancer treatment.

Nardosinone and aurantio-obtusin, two medicine food homology natural compounds, are anti-influenza agents as indicated by transcriptome signature reversion

BACKGROUND

Medicine food homology (MFH) refers to food that can be used as medicine, and compounds isolated from MFH materials are valuable in novel drug discovery due to their good safety. Transcriptome signature reversion (TSR) is an attractive method for discovering drugs through transcriptional reverse matching; namely, the changes in transcriptional signatures induced by compounds are matched to a certain disease. This strategy can be used to discover anti-influenza agents among MFH natural compounds.

PURPOSE

MFH natural compounds with anti-influenza activities were identified through analyses of the reversal in the expression of multiple informative genes followed by in vitro evaluation of the cytopathic effect (CPE) caused by influenza infection and relative quantification of the nucleoprotein (NP) gene in viral RNA (vRNA). The combined effect of active compounds was determined through network-based separation score prediction followed by quantification of the viral hemagglutinin (HA) level.

METHODS

The transcriptome profiles of 4 lung or airway cell lines infected with 7 influenza virus strains were analyzed by robust rank aggregation (RRA) to identify informative genes in the signature of influenza virus infection. The identified informative genes were then matched to a transcriptomic profile library of MFH natural compounds. The anti-influenza activities of MFH natural compounds with negative enrichment scores (ESs) were evaluated in vitro using a CPE assay and relative quantification of the NP gene in the vRNA in the supernatant and cytoplasm to identify anti-influenza agents. The effects of combinations of active compounds were analyzed using network-based calculations followed by confirmation through bioassays for quantifying the viral HA levels.

RESULTS

Among the 159 MFH natural compounds, 54 compounds had negative ESs, as determined through TSR, and the anti-influenza activities of nardosinone and aurantio-obtusin were confirmed by bioassays. The half-maximal effective concentrations (EC₅₀) of nardosinone and aurantio-obtusin were 4.3-84.4 μM and 31.9-113.6 μM, respectively. The separation score between the informative genes with expression that was negatively regulated by nardosinone and aurantio-obtusin in the human protein-protein interaction (PPI) network was calculated to be 0.10, which indicated that the two compounds potentially exert a synergistic effect, and this effect was confirmed by the finding that the combination indexes (CIs) were calculated to equal 0.86 at inhibition level of 50% and 0.44 at inhibition level of 90%.

CONCLUSION

The TSR analysis and in vitro evaluation identified nardosinone and aurantio-obtusin as anti-influenza agents. Their antiviral activities were exerted by reversing the expression of multiple informative genes of the host cells. The separation analysis between the informative genes that were reversely regulated by nardosinone and aurantio-obtusin indicated that their combination may exert a synergistic effect, which was confirmed in vitro.

Roles of Skeletal Muscle in Development: A Bioinformatics and Systems Biology Overview

The ability to assess various cellular events consequent to perturbations, such as genetic mutations, disease states and therapies, has been recently revolutionized by technological advances in multiple "omics" fields. The resulting deluge of information has enabled and necessitated the development of tools required to both process and interpret the data. While of tremendous value to basic researchers, the amount and complexity of the data has made it extremely difficult to manually draw inference and identify factors key to the study objectives. The challenges of data reduction and interpretation are being met by the development of increasingly complex tools that integrate disparate knowledge bases and synthesize coherent models based on current biological understanding. This chapter presents an example of how genomics data can be integrated with biological network analyses to gain further insight into the developmental consequences of genetic perturbations. State of the art methods for conducting similar studies are discussed along with modern methods used to analyze and interpret the data.

A comprehensive SARS-CoV-2-human protein-protein interactome reveals COVID-19 pathobiology and potential host therapeutic targets

Studying viral-host protein-protein interactions can facilitate the discovery of therapies for viral infection. We use high-throughput yeast two-hybrid experiments and mass spectrometry to generate a comprehensive SARS-CoV-2-human protein-protein interactome network consisting of 739 high-confidence binary and co-complex interactions, validating 218 known SARS-CoV-2 host factors and revealing 361 novel ones. Our results show the highest overlap of interaction partners between published datasets and of genes differentially expressed in samples from COVID-19 patients. We identify an interaction between the viral protein ORF3a and the human transcription factor ZNF579, illustrating a direct viral impact on host transcription. We perform network-based screens of >2,900 FDA-approved or investigational drugs and identify 23 with significant network proximity to SARS-CoV-2 host factors. One of these drugs, carvedilol, shows clinical benefits for COVID-19 patients in an electronic health records analysis and antiviral properties in a human lung cell line infected with SARS-CoV-2. Our study demonstrates the value of network systems biology to understand human-virus interactions and provides hits for further research on COVID-19 therapeutics.

Identify Biomarkers and Design Effective Multi-Target Drugs in Ovarian Cancer: Hit Network-Target Sets Model Optimizing

Epithelial ovarian cancer (EOC) is highly aggressive with poor patient outcomes, and a deeper understanding of ovarian cancer tumorigenesis could help guide future treatment development. We proposed an optimized hit network-target sets model to systematically characterize the underlying pathological mechanisms and intra-tumoral heterogeneity in human ovarian cancer. Using TCGA data, we constructed an epithelial ovarian cancer regulatory network in this study. We use three distinct methods to produce different HNSs for identification of the driver genes/nodes, core modules, and core genes/nodes. Following the creation of the optimized HNS (OHNS) by the integration of DN (driver nodes), CM (core module), and CN (core nodes), the effectiveness of various HNSs was assessed based on the significance of the network topology, control potential, and clinical value. Immunohistochemical (IHC), qRT-PCR, and Western blotting were adopted to measure the expression of hub genes and proteins involved in epithelial ovarian cancer (EOC). We discovered that the OHNS has two key advantages: the network's central location and controllability. It also plays a significant role in the illness network due to its wide range of capabilities. The OHNS and clinical samples revealed the endometrial cancer signaling, and the PI3K/AKT, NER, and BMP pathways. MUC16, FOXA1, FBXL2, ARID1A, COX15, COX17, SCO1, SCO2, NDUFA4L2, NDUFA, and PTEN hub genes were predicted and may serve as potential candidates for new treatments and biomarkers for EOC. This research can aid in better capturing the disease progression, the creation of potent multi-target medications, and the direction of the therapeutic community in the optimization of effective treatment regimens by various research objectives in cancer treatment.

Machine learning prediction of side effects for drugs in clinical trials

Early and accurate detection of side effects is critical for the clinical success of drugs under development. Here, we aim to predict unknown side effects for drugs with a small number of side effects identified in randomized controlled clinical trials. Our machine learning framework, the geometric self-expressive model (GSEM), learns globally optimal self-representations for drugs and side effects from pharmacological graph networks. We show the usefulness of the GSEM on 505 therapeutically diverse drugs and 904 side effects from multiple human physiological systems. Here, we also show a data integration strategy that could be adopted to improve the ability of side effect prediction models to identify unknown side effects that might only appear after the drug enters the market.

Potential bioactive compounds and mechanisms of Fibraurea recisa Pierre for the treatment of Alzheimer's disease analyzed by network pharmacology and molecular docking prediction

Introduction

Heat-clearing and detoxifying Chinese medicines have been documented to have anti-Alzheimer's disease (AD) activities according to the accumulated clinical experience and pharmacological research results in recent decades. In this study, Fibraurea recisa Pierre (FRP), the classic type of Heat-clearing and detoxifying Chinese medicine, was selected as the object of research.

Methods

12 components with anti-AD activities were identified in FRP by a variety of methods, including silica gel column chromatography, multiple databases, and literature searches. Then, network pharmacology and molecular docking were adopted to systematically study the potential anti-AD mechanism of these compounds. Consequently, it was found that these 12 compounds could act on 235 anti-AD targets, of which AKT and other targets were the core targets. Meanwhile, among these 235 targets, 71 targets were identified to be significantly correlated with the pathology of amyloid beta (Aβ) and Tau.

Results and discussion

In view of the analysis results of the network of active ingredients and targets, it was observed that palmatine, berberine, and other alkaloids in FRP were the key active ingredients for the treatment of AD. Further, Kyoto encyclopedia of genes and genomes (KEGG) pathway enrichment analysis revealed that the neuroactive ligand-receptor interaction pathway and PI3K-Akt signaling pathway were the most significant signaling pathways for FRP to play an anti-AD role. Findings in our study suggest that multiple primary active ingredients in FRP can play a multitarget anti-AD effect by regulating key physiological processes such as neurotransmitter transmission and anti-inflammation. Besides, key ingredients such as palmatine and berberine in FRP are expected to be excellent leading compounds of multitarget anti-AD drugs.

DrugRepo: a novel approach to repurposing drugs based on chemical and genomic features

The drug development process consumes 9-12 years and approximately one billion US dollars in costs. Due to the high finances and time costs required by the traditional drug discovery paradigm, repurposing old drugs to treat cancer and rare diseases is becoming popular. Computational approaches are mainly data-driven and involve a systematic analysis of different data types leading to the formulation of repurposing hypotheses. This study presents a novel scoring algorithm based on chemical and genomic data to repurpose drugs for 669 diseases from 22 groups, including various cancers, musculoskeletal, infections, cardiovascular, and skin diseases. The data types used to design the scoring algorithm are chemical structures, drug-target interactions (DTI), pathways, and disease-gene associations. The repurposed scoring algorithm is strengthened by integrating the most comprehensive manually curated datasets for each data type. At DrugRepo score ≥ 0.4, we repurposed 516 approved drugs across 545 diseases. Moreover, hundreds of novel predicted compounds can be matched with ongoing studies at clinical trials. Our analysis is supported by a web tool available at: http://drugrepo.org/ .

Network-based prediction of the disclosure of ideation about self-harm and suicide in online counseling sessions

BACKGROUND

In psychological services, the transition to the disclosure of ideation about self-harm and suicide (ISS) is a critical point warranting attention. This study developed and tested a succinct descriptor to predict such transitions in an online synchronous text-based counseling service.

METHOD

We analyzed two years' worth of counseling sessions (N = 49,770) from Open Up, a 24/7 service in Hong Kong. Sessions from Year 1 (N = 20,618) were used to construct a word affinity network (WAN), which depicts the semantic relationships between words. Sessions from Year 2 (N = 29,152), including 1168 with explicit ISS, were used to train and test the downstream ISS prediction model. We divided and classified these sessions into ISS blocks (ISSBs), blocks prior to ISSBs (PISSBs), and non-ISS blocks (NISSBs). To detect PISSB, we adopted complex network approaches to examine the distance among different types of blocks in WAN.

RESULTS

Our analyses find that words within a block tend to form a module in WAN and that network-based distance between modules is a reliable indicator of PISSB. The proposed model yields a c-statistic of 0.79 in identifying PISSB.

CONCLUSIONS

This simple yet robust network-based model could accurately predict the transition point of suicidal ideation prior to its explicit disclosure. It can potentially improve the preparedness and efficiency of help-providers in text-based counseling services for mitigating self-harm and suicide.

The Use of Zidovudine Pharmacophore in Multi-Target-Directed Ligands for AIDS Therapy

The concept of polypharmacology embraces multiple drugs combined in a therapeutic regimen (drug combination or cocktail), fixed dose combinations (FDCs), and a single drug that binds to different targets (multi-target drug). A polypharmacology approach is widely applied in the treatment of acquired immunodeficiency syndrome (AIDS), providing life-saving therapies for millions of people living with HIV. Despite the success in viral load suppression and patient survival of combined antiretroviral therapy (cART), the development of new drugs has become imperative, owing to the emergence of resistant strains and poor adherence to cART. 3'-azido-2',3'-dideoxythymidine, also known as azidothymidine or zidovudine (AZT), is a widely applied starting scaffold in the search for new compounds, due to its good antiretroviral activity. Through the medicinal chemistry tool of molecular hybridization, AZT has been included in the structure of several compounds allowing for the development of multi-target-directed ligands (MTDLs) as antiretrovirals. This review aims to systematically explore and critically discuss AZT-based compounds as potential MTDLs for the treatment of AIDS. The review findings allowed us to conclude that: (i) AZT hybrids are still worth exploring, as they may provide highly active compounds targeting different steps of the HIV-1 replication cycle; (ii) AZT is a good starting point for the preparation of co-drugs with enhanced cell permeability.

DeepIDC: A Prediction Framework of Injectable Drug Combination Based on Heterogeneous Information and Deep Learning

BACKGROUND AND OBJECTIVE

In clinical practice, injectable drug combination (IDC) usually provides good therapeutic effects for patients. Numerous clinical studies have directly indicated that inappropriate IDC generates adverse drug events (ADEs). The clinical application of injections is increasing, and many injections lack relevant combination information. It is still a significant need for experienced clinical pharmacists to participate in evidence-based drug decision making, monitor medication safety, and manage drug interactions. Meanwhile, a large number of injection pairs and dosage combinations limit exhaustive screening. Here, we present a prediction framework, called DeepIDC, that can expediently screen the feasibility of IDCs using heterogeneous information with deep learning. This is the first specific prediction framework to identify IDCs.

METHODS

Since the interaction between the injected drugs may occur in the direct physical and chemical reactions at the time of mixing or may be the indirect interaction of their drug targets and pathways, we used molecular fingerprints, drug-target associations, and drug-pathway associations to convert injections into a string of digital vectors. Then, based on these injection vectors, we combined a bidirectional long short-term memory and a feed-forward neural network to build a prediction model for accurate and instructive prediction of IDC.

RESULTS

In three realistic evaluation scenarios, DeepIDC has achieved ideal prediction results. Furthermore, compared with the other five machine-learning methods, the proposed predictor is more efficient and robust. Among the top 30 potential IDCs of each IDC class predicted by DeepIDC, we found that 9 cases were experimentally verified in the literature or available on Drug.com.

CONCLUSION

The information we extracted in vivo and in vitro can effectively characterize injectable drugs. DeepIDC developed based on deep learning algorithm provides a valuable unified framework for new IDC discovery, which can make up for the lack of IDC information and predict potential IDC events.

Graph representation learning in biomedicine and healthcare

Networks-or graphs-are universal descriptors of systems of interacting elements. In biomedicine and healthcare, they can represent, for example, molecular interactions, signalling pathways, disease co-morbidities or healthcare systems. In this Perspective, we posit that representation learning can realize principles of network medicine, discuss successes and current limitations of the use of representation learning on graphs in biomedicine and healthcare, and outline algorithmic strategies that leverage the topology of graphs to embed them into compact vectorial spaces. We argue that graph representation learning will keep pushing forward machine learning for biomedicine and healthcare applications, including the identification of genetic variants underlying complex traits, the disentanglement of single-cell behaviours and their effects on health, the assistance of patients in diagnosis and treatment, and the development of safe and effective medicines.

Multi-omics peripheral and core regions of cancer

Thousands of genes are perturbed by cancer, and these disturbances can be seen in transcriptome, methylation, somatic mutation, and copy number variation omics studies. Understanding their connectivity patterns as an omnigenic neighbourhood in a molecular interaction network (interactome) is a key step towards advancing knowledge of the molecular mechanisms underlying cancers. Here, we introduce a unified connectivity line (CLine) to pinpoint omics-specific omnigenic patterns across 15 curated cancers. Taking advantage of the universality of CLine, we distinguish the peripheral and core genes for each omics aspect. We propose a network-based framework, multi-omics periphery and core (MOPC), to combine peripheral and core genes from different omics into a button-like structure. On the basis of network proximity, we provide evidence that core genes tend to be specifically perturbed in one omics, but the peripheral genes are diversely perturbed in multiple omics. And the core of one omics is regulated by multiple omics peripheries. Finally, we take the MOPC as an omnigenic neighbourhood, describe its characteristics, and explore its relative contribution to network-based mechanisms of cancer. We were able to present how multi-omics perturbations percolate through the human interactome and contribute to an integrated periphery and core.

Interpretable deep learning translation of GWAS and multi-omics findings to identify pathobiology and drug repurposing in Alzheimer's disease

Translating human genetic findings (genome-wide association studies [GWAS]) to pathobiology and therapeutic discovery remains a major challenge for Alzheimer's disease (AD). We present a network topology-based deep learning framework to identify disease-associated genes (NETTAG). We leverage non-coding GWAS loci effects on quantitative trait loci, enhancers and CpG islands, promoter regions, open chromatin, and promoter flanking regions under the protein-protein interactome. Via NETTAG, we identified 156 AD-risk genes enriched in druggable targets. Combining network-based prediction and retrospective case-control observations with 10 million individuals, we identified that usage of four drugs (ibuprofen, gemfibrozil, cholecalciferol, and ceftriaxone) is associated with reduced likelihood of AD incidence. Gemfibrozil (an approved lipid regulator) is significantly associated with 43% reduced risk of AD compared with simvastatin using an active-comparator design (95% confidence interval 0.51-0.63, p < 0.0001). In summary, NETTAG offers a deep learning methodology that utilizes GWAS and multi-genomic findings to identify pathobiology and drug repurposing in AD.

Deep learning and multi-omics approach to predict drug responses in cancer

BACKGROUND

Cancers are genetically heterogeneous, so anticancer drugs show varying degrees of effectiveness on patients due to their differing genetic profiles. Knowing patient's responses to numerous cancer drugs are needed for personalized treatment for cancer. By using molecular profiles of cancer cell lines available from Cancer Cell Line Encyclopedia (CCLE) and anticancer drug responses available in the Genomics of Drug Sensitivity in Cancer (GDSC), we will build computational models to predict anticancer drug responses from molecular features.

RESULTS

We propose a novel deep neural network model that integrates multi-omics data available as gene expressions, copy number variations, gene mutations, reverse phase protein array expressions, and metabolomics expressions, in order to predict cellular responses to known anti-cancer drugs. We employ a novel graph embedding layer that incorporates interactome data as prior information for prediction. Moreover, we propose a novel attention layer that effectively combines different omics features, taking their interactions into account. The network outperformed feedforward neural networks and reported 0.90 for [Formula: see text] values for prediction of drug responses from cancer cell lines data available in CCLE and GDSC.

CONCLUSION

The outstanding results of our experiments demonstrate that the proposed method is capable of capturing the interactions of genes and proteins, and integrating multi-omics features effectively. Furthermore, both the results of ablation studies and the investigations of the attention layer imply that gene mutation has a greater influence on the prediction of drug responses than other omics data types. Therefore, we conclude that our approach can not only predict the anti-cancer drug response precisely but also provides insights into reaction mechanisms of cancer cell lines and drugs as well.

Molecular Insights into the Heterotropic Allosteric Mechanism in Cytochrome P450 3A4-Mediated Midazolam Metabolism

Cytochrome P450 3A4 (CYP3A4) is the main P450 enzyme for drug metabolism and drug-drug interactions (DDIs), as it is involved in the metabolic process of approximately 50% of drugs. A detailed mechanistic elucidation of DDIs mediated by CYP3A4 is commonly believed to be critical for drug optimization and rational use. Here, two typical probes, midazolam (MDZ, substrate) and testosterone (TST, allosteric effector), are used to investigate the molecular mechanism of CYP3A4-mediated heterotropic allosteric interactions, through conventional molecular dynamics (cMD) and well-tempered metadynamics (WT-MTD) simulations. Distance monitoring shows that TST can stably bind in two potential peripheral sites (Site 1 and Site 2) of CYP3A4. The binding of TST at these two sites can induce conformational changes in CYP3A4 flexible loops on the basis of conformational analysis, thereby promoting the transition of the MDZ binding mode and affecting the ratio of MDZ metabolites. According to the results of the residue interaction network, multiple allosteric communication pathways are identified that can provide vivid and applicable insights into the heterotropic allostery of TST on MDZ metabolism. Comparing the regulatory effects and the communication pathways, the allosteric effect caused by TST binding in Site 2 seems to be more pronounced than in Site 1. Our findings could provide a deeper understanding of CYP3A4-mediated heterotropic allostery at the atomic level and would be helpful for rational drug use as well as the design of new allosteric modulators.

Computational Tactics for Precision Cancer Network Biology

Network biology has garnered tremendous attention in understanding complex systems of cancer, because the mechanisms underlying cancer involve the perturbations in the specific function of molecular networks, rather than a disorder of a single gene. In this article, we review the various computational tactics for gene regulatory network analysis, focused especially on personalized anti-cancer therapy. This paper covers three major topics: (1) cell line's (or patient's) cancer characteristics specific gene regulatory network estimation, which enables us to reveal molecular interplays under varying conditions of cancer characteristics of cell lines (or patient); (2) computational approaches to interpret the multitudinous and massive networks; (3) network-based application to uncover molecular mechanisms of cancer and related marker identification. We expect that this review will help readers understand personalized computational network biology that plays a significant role in precision cancer medicine.

Jha P, Lupieri A, Nakano T, Passos LS, Rogers MA, Becker-Greene D, Le TD, Decano JL, Ho Lee L, Guimaraes GC, Abdelhamid I, Halu A, Muscoloni A, V. Cannistraci C, Higashi H, Zhang H, Vromman A, Libby P, Keith Ozaki C, Sharma A, Singh SA, Aikawa E, Aikawa M. Proprotein Convertase Subtilisin/Kexin 9 (PCSK9) Promotes Macrophage Activation via LDL Receptor-Independent Mechanisms

BACKGROUND

Activated macrophages contribute to the pathogenesis of vascular disease. Vein graft failure is a major clinical problem with limited therapeutic options. PCSK9 (proprotein convertase subtilisin/kexin 9) increases low-density lipoprotein (LDL)-cholesterol levels via LDL receptor (LDLR) degradation. The role of PCSK9 in macrophage activation and vein graft failure is largely unknown, especially through LDLR-independent mechanisms. This study aimed to explore a novel mechanism of macrophage activation and vein graft disease induced by circulating PCSK9 in an LDLR-independent fashion.

METHODS

We used Ldlr^-/- mice to examine the LDLR-independent roles of circulating PCSK9 in experimental vein grafts. Adeno-associated virus (AAV) vector encoding a gain-of-function mutant of PCSK9 (rAAV8/D377Y-mPCSK9) induced hepatic PCSK9 overproduction. To explore novel inflammatory targets of PCSK9, we used systems biology in Ldlr^-/- mouse macrophages.

RESULTS

In Ldlr^-/- mice, AAV-PCSK9 increased circulating PCSK9, but did not change serum cholesterol and triglyceride levels. AAV-PCSK9 promoted vein graft lesion development when compared with control AAV. In vivo molecular imaging revealed that AAV-PCSK9 increased macrophage accumulation and matrix metalloproteinase activity associated with decreased fibrillar collagen, a molecular determinant of atherosclerotic plaque stability. AAV-PCSK9 induced mRNA expression of the pro-inflammatory mediators IL-1β (interleukin-1 beta), TNFα (tumor necrosis factor alpha), and MCP-1 (monocyte chemoattractant protein-1) in peritoneal macrophages underpinned by an in vitro analysis of Ldlr^-/- mouse macrophages stimulated with endotoxin-free recombinant PCSK9. A combination of unbiased global transcriptomics and new network-based hyperedge entanglement prediction analysis identified the NF-κB (nuclear factor-kappa B) signaling molecules, lectin-like oxidized LOX-1 (LDL receptor-1), and SDC4 (syndecan-4) as potential PCSK9 targets mediating pro-inflammatory responses in macrophages.

CONCLUSIONS

Circulating PCSK9 induces macrophage activation and vein graft lesion development via LDLR-independent mechanisms. PCSK9 may be a potential target for pharmacologic treatment for this unmet medical need.

Exploratory Analysis of the Sasang Constitution by Combining Network Analysis and Information Entropy

Sasang constitutional medicine is a unique concept in Korean medicine that can provide valuable insights into personalized healthcare and disease treatment. In this study, we combined network analysis and information entropy to systematically investigate the related information of Sasang constitutional (SC) types. A feature network was constructed using SC type and clinical information. The SC type-associated features and feature classes were identified using statistical analysis and entropy ranking. The patient network was constructed based on SC-type-associated features. We found that the feature network was closely connected within the features of the same classes and between several feature class pairs, including the symptom class. Most of the separation values between the feature classes, including the symptom class, were negative. In addition, we found 42 clinical features related to the SC type, and two important classes -personality and cold/heat- that increase the entropy ranking of the SC type. In the patient network, we found sparsely connected modules between SC types and a positive separation value between the Taeeumin–Soeumin and Taeeumin–Soyangin pairs. Our data-driven approach provides a deeper understanding of modernized forms of SC types and suggests that SC type is a practically useful concept for stratified healthcare and personalized medicine.

A network paradigm predicts drug synergistic effects using downstream protein-protein interactions

In some cases, drug combinations affect adverse outcome phenotypes by binding the same protein; however, drug-binding proteins are associated through protein-protein interaction (PPI) networks within the cell, suggesting that drug phenotypes may result from long-range network effects. We first used PPI network analysis to classify drugs based on proteins downstream of their targets and next predicted drug combination effects where drugs shared network proteins but had distinct binding proteins (e.g., targets, enzymes, or transporters). By classifying drugs using their downstream proteins, we had an 80.7% sensitivity for predicting rare drug combination effects documented in gold-standard datasets. We further measured the effect of predicted drug combinations on adverse outcome phenotypes using novel observational studies in the electronic health record. We tested predictions for 60 network-drug classes on seven adverse outcomes and measured changes in clinical outcomes for predicted combinations. These results demonstrate a novel paradigm for anticipating drug synergistic effects using proteins downstream of drug targets.

An optimized herbal combination for the treatment of liver fibrosis: Hub genes, bioactive ingredients, and molecular mechanisms

ETHNOPHARMACOLOGICAL RELEVANCE

Liver fibrosis is a chronic liver disease that can lead to cirrhosis, liver failure, and hepatocellular carcinoma, and it is associated with long-term adverse outcomes and mortality. As a primary resource for complementary and alternative medicine, traditional Chinese medicine (TCM) has accumulated a large number of effective formulas for the treatment of liver fibrosis in clinical practice. However, studies on how to systematically optimize TCM formulas are still lacking.

AIM OF THE REVIEW

To provide a methodological reference for the systematic optimization of TCM formulae against liver fibrosis and explored the underlying molecular mechanisms; To provide an efficient method for searching for lead compounds from natural sources and developing from herbal medicines; To enable clinicians and patients to make more reasonable choices and promote the effective treatment toward those patients with liver fibrosis.

MATERIALS AND METHODS

TCM formulas related to treating liver fibrosis were collected from the Web of Science, PubMed, the China National Knowledge Infrastructure (CNKI), Wan Fang, and the Chinese Scientific Journals Database (VIP). Furthermore, the TCM compatibility patterns were mined using association analysis. The core TCM combinations were found by designing an optimized formulas algorithm. Finally, the hub target proteins, potential molecular mechanisms, and active compounds were explored through integrative pharmacology and docking-based inverse virtual screening (IVS) approaches.

RESULTS

We found that the herbs for reinforcing deficiency, activating blood, removing blood stasis, and clearing heat were the basis of TCM formulae patterns. Furthermore, the combination of Salviae Miltiorrhizae (Salvia miltiorrhiza Bunge; Chinese salvia/Danshen), Astragali Radix (Astragalus membranaceus (Fisch.) Bunge; Astragalus/Huangqi), and Radix Bupleuri (Bupleurum chinense DC.; Bupleurum/Chaihu) was identified as core groups. A total of six targets (TNF, STAT3, EGFR, IL2, ICAM1, PTGS2) play a pivotal role in TCM-mediated liver fibrosis inhibition. (-)-Cryptotanshinone, Tanshinaldehyde, Ononin, Thymol, Daidzein, and Formononetin were identified as active compounds in TCM. And mechanistically, TCM could affect the development of liver fibrosis by regulating inflammation, immunity, angiogenesis, antioxidants, and involvement in TNF, MicroRNAs, Jak-STAT, NF-kappa B, and C-type lectin receptors (CLRs) signaling pathways. Molecular docking results showed that key components had good potential to bind to the target genes.

CONCLUSION

In summary, this study provides a methodological reference for the systematic optimization of TCM formulae and exploration of underlying molecular mechanisms.

PINet 1.0: A pathway network-based evaluation of drug combinations for the management of specific diseases

Drug combinations can increase the therapeutic effect by reducing the level of toxicity and the occurrence of drug resistance. Therefore, several drug combinations are often used in the management of complex diseases. However, due to the exponential growth in drug development, it would be impractical to evaluate all combinations through experiments. In view of this, we developed Pathway Interaction Network (PINet) biological model to estimate the optimal drug combinations for various diseases. The random walk with restart (RWR) algorithm was used to capture the "disease state" and "drug state," while PINet was used to evaluate the optimal drug combinations and the high-order drug combination. The model achieved a mean area under the curve of a receiver operating characteristic curve of 0.885. In addition, for some diseases, PINet predicted the optimal drug combination. For example, in the case of acute myeloid leukemia, PINet correctly predicted midostaurin and gemtuzumab as effective drug combinations, as demonstrated by the results of a Phase-I clinical trial. Moreover, PINet also correctly predicted the potential drug combinations for diseases that lacked a training dataset that could not be predicted using standard machine learning models.

Reprogramming of myeloid cells and their progenitors in patients with non-medullary thyroid carcinoma

Myeloid cells, crucial players in antitumoral defense, are affected by tumor-derived factors and treatment. The role of myeloid cells and their progenitors prior to tumor infiltration is poorly understood. Here we show single-cell transcriptomics and functional analyses of the myeloid cell lineage in patients with non-medullary thyroid carcinoma (TC) and multinodular goiter, before and after treatment with radioactive iodine compared to healthy controls. Integrative data analysis indicates that monocytes of TC patients have transcriptional upregulation of antigen presentation, reduced cytokine production capacity, and overproduction of reactive oxygen species. Interestingly, these cancer-related pathological changes are partially removed upon treatment. In bone marrow, TC patients tend to shift from myelopoiesis towards lymphopoiesis, reflected in transcriptional differences. Taken together, distinct transcriptional and functional changes in myeloid cells arise before their infiltration of the tumor and are already initiated in bone marrow, which suggests an active role in forming the tumor immune microenvironment.

Recent computational drug repositioning strategies against SARS-CoV-2

•

We performed a comprehensive review of computational drug repositioning methods applied to COVID-19 based on differing data types including sequence data, expression data, structure data and interaction data.

•

We found that graph theory and neural network were the most used strategies for drug repositioning in the case of COVID-19.

•

Integrating different levels of data may improve the success rate for drug repositioning.

Since COVID-19 emerged in 2019, significant levels of suffering and disruption have been caused on a global scale. Although vaccines have become widely used, the virus has shown its potential for evading immunities or acquiring other novel characteristics. Whether current drug treatments are still effective for people infected with Omicron remains unclear. Due to the long development cycles and high expense requirements of de novo drug development, many researchers have turned to consider drug repositioning in the search to find effective treatments for COVID-19. Here, we review such drug repositioning and combination efforts towards providing better handling. For potential drugs under consideration, aspects of both structure and function require attention, with specific categories of sequence, expression, structure, and interaction, the key parameters for investigation. For different data types, we show the corresponding differing drug repositioning methods that have been exploited. As incorporating drug combinations can increase therapeutic efficacy and reduce toxicity, we also review computational strategies to reveal drug combination potential. Taken together, we found that graph theory and neural network were the most used strategy with high potential towards drug repositioning for COVID-19. Integrating different levels of data may further improve the success rate of drug repositioning.

DrugMAP: molecular atlas and pharma-information of all drugs

The efficacy and safety of drugs are widely known to be determined by their interactions with multiple molecules of pharmacological importance, and it is therefore essential to systematically depict the molecular atlas and pharma-information of studied drugs. However, our understanding of such information is neither comprehensive nor precise, which necessitates the construction of a new database providing a network containing a large number of drugs and their interacting molecules. Here, a new database describing the molecular atlas and pharma-information of drugs (DrugMAP) was therefore constructed. It provides a comprehensive list of interacting molecules for >30 000 drugs/drug candidates, gives the differential expression patterns for >5000 interacting molecules among different disease sites, ADME (absorption, distribution, metabolism and excretion)-relevant organs and physiological tissues, and weaves a comprehensive and precise network containing >200 000 interactions among drugs and molecules. With the great efforts made to clarify the complex mechanism underlying drug pharmacokinetics and pharmacodynamics and rapidly emerging interests in artificial intelligence (AI)-based network analyses, DrugMAP is expected to become an indispensable supplement to existing databases to facilitate drug discovery. It is now fully and freely accessible at: https://idrblab.org/drugmap/.

CTpathway: a CrossTalk-based pathway enrichment analysis method for cancer research

BACKGROUND

Pathway enrichment analysis (PEA) is a common method for exploring functions of hundreds of genes and identifying disease-risk pathways. Moreover, different pathways exert their functions through crosstalk. However, existing PEA methods do not sufficiently integrate essential pathway features, including pathway crosstalk, molecular interactions, and network topologies, resulting in many risk pathways that remain uninvestigated.

METHODS

To overcome these limitations, we develop a new crosstalk-based PEA method, CTpathway, based on a global pathway crosstalk map (GPCM) with >440,000 edges by combing pathways from eight resources, transcription factor-gene regulations, and large-scale protein-protein interactions. Integrating gene differential expression and crosstalk effects in GPCM, we assign a risk score to genes in the GPCM and identify risk pathways enriched with the risk genes.

RESULTS

Analysis of >8300 expression profiles covering ten cancer tissues and blood samples indicates that CTpathway outperforms the current state-of-the-art methods in identifying risk pathways with higher accuracy, reproducibility, and speed. CTpathway recapitulates known risk pathways and exclusively identifies several previously unreported critical pathways for individual cancer types. CTpathway also outperforms other methods in identifying risk pathways across all cancer stages, including early-stage cancer with a small number of differentially expressed genes. Moreover, the robust design of CTpathway enables researchers to analyze both bulk and single-cell RNA-seq profiles to predict both cancer tissue and cell type-specific risk pathways with higher accuracy.

CONCLUSIONS

Collectively, CTpathway is a fast, accurate, and stable pathway enrichment analysis method for cancer research that can be used to identify cancer risk pathways. The CTpathway interactive web server can be accessed here http://www.jianglab.cn/CTpathway/ . The stand-alone program can be accessed here https://github.com/Bioccjw/CTpathway .

Transcriptomics-based network medicine approach identifies metformin as a repurposable drug for atrial fibrillation

Effective drugs for atrial fibrillation (AF) are lacking, resulting in significant morbidity and mortality. This study demonstrates that network proximity analysis of differentially expressed genes from atrial tissue to drug targets can help prioritize repurposed drugs for AF. Using enrichment analysis of drug-gene signatures and functional testing in human inducible pluripotent stem cell (iPSC)-derived atrial-like cardiomyocytes, we identify metformin as a top repurposed drug candidate for AF. Using the active compactor, a new design analysis of large-scale longitudinal electronic health record (EHR) data, we determine that metformin use is significantly associated with a reduced risk of AF (odds ratio = 0.48, 95%, confidence interval [CI] 0.36-0.64, p < 0.001) compared with standard treatments for diabetes. This study utilizes network medicine methodologies to identify repurposed drugs for AF treatment and identifies metformin as a candidate drug.

Combinatorial approaches for novel cardiovascular drug discovery: a review of the literature

INTRODUCTION

In this article, authors report an inclusive discussion about the combinatorial approach for the treatment of cardiovascular diseases (CVDs) and for counteracting the cardiovascular risk factors. The mentioned strategy was demonstrated to be useful for improving the efficacy of pharmacological treatments and in CVDs showed superior efficacy with respect to the classical monotherapeutic approach.

AREAS COVERED

According to this topic, authors analyzed the combinatorial treatments that are available on the market, highlighting clinical studies that demonstrated the efficacy of combinatorial drug strategies to cure CVDs and related risk factors. Furthermore, the review gives an outlook on the future perspective of this therapeutic option, highlighting novel drug targets and disease models that could help the future cardiovascular drug discovery.

EXPERT OPINION

The use of specifically designed and increasingly rational and effective drug combination therapies can therefore be considered the evolution of polypharmacy in cardiometabolic and CVDs. This approach can allow to intervene on multiple etiopathogenetic mechanisms of the disease or to act simultaneously on different pathologies/risk factors, using the combinations most suitable from a pharmacodynamic, pharmacokinetic, and toxicological perspective, thus finding the most appropriate therapeutic option.

A genetic map of the chromatin regulators to drug response in cancer cells

BACKGROUND

Diverse drug vulnerabilities owing to the Chromatin regulators (CRs) genetic interaction across various cancers, but the identification of CRs genetic interaction remains challenging.

METHODS

In order to provide a global view of the CRs genetic interaction in cancer cells, we developed a method to identify potential drug response-related CRs genetic interactions for specific cancer types by integrating the screen of CRISPR-Cas9 and pharmacogenomic response datasets.

RESULTS

Totally, 625 drug response-related CRs synthetic lethality (CSL) interactions and 288 CRs synthetic viability (CSV) interactions were detected. Systematically network analysis presented CRs genetic interactions have biological function relationship. Furthermore, we validated CRs genetic interactions induce multiple omics deregulation in The Cancer Genome Atlas. We revealed the colon adenocarcinoma patients (COAD) with mutations of a CRs set (EP300, MSH6, NSD2 and TRRAP) mediate a better survival with low expression of MAP2 and could benefit from taxnes. While the COAD patients carrying at least one of the CSV interactions in Vorinostat CSV module confer a poor prognosis and may be resistant to Vorinostat treatment.

CONCLUSIONS

The CRs genetic interaction map provides a rich resource to investigate cancer-associated CRs genetic interaction and proposes a powerful strategy of biomarker discovery to guide the rational use of agents in cancer therapy.