OncoLoop: A Network-Based Precision Cancer Medicine Framework

Network analysis to identify driver genes and combination drugs in brain cancer

Brain cancer is one of the deadliest diseases, although many efforts have been made to treat it, there is no comprehensive and effective treatment approach yet. In recent years, the use of network-based analysis to identify important biological genes and pathways involved in various complex diseases, including brain cancer, has attracted the attention of researchers. The goal of this manuscript is to perform a comprehensive analysis of the various results presented related to brain cancer. For this purpose, firstly, based on the CORMINE medical database, collected all the genes related to brain cancer with a valid P-value. Then the structural and functional relationships between the above gene sets have been identified based on the STRING database. Next, in the PPI network, hub centrality analysis was performed to determine the proteins that have many connections with other proteins. After the modularization of the network, the module with the most hub vertices is considered as the most relevant module to the formation and progression of brain cancer. Since the driver vertices play an important role in biological systems, the edges of the selected module were oriented, and by analyzing the controllability of complex networks, a set of five proteins with the highest control power has been identified. Finally, based on the drug-gene interaction, a set of drugs effective on each of the driver genes has been obtained, which can potentially be used as new combination drugs. Validation of the hub and driver proteins shows that they are mainly essential proteins in the biological processes related to the various cancers and therefore the drugs that affect them can be considered as new combination therapy. The presented procedure can be used for any other complex disease.

Computer-aided drug discovery strategies for novel therapeutics for prostate cancer leveraging next-generating sequencing data

INTRODUCTION

Prostate cancer (PC) is the most common malignancy and accounts for a significant proportion of cancer deaths among men. Although initial therapy success can often be observed in patients diagnosed with localized PC, many patients eventually develop disease recurrence and metastasis. Without effective treatments, patients with aggressive PC display very poor survival. To curb the current high mortality rate, many investigations have been carried out to identify efficacious therapeutics. Compared to de novo drug designs, computational methods have been widely employed to offer actionable drug predictions in a fast and cost-efficient way. Particularly, powered by an increasing availability of next-generation sequencing molecular profiles from PC patients, computer-aided approaches can be tailored to screen for candidate drugs.

AREAS COVERED

Herein, the authors review the recent advances in computational methods for drug discovery utilizing molecular profiles from PC patients. Given the uniqueness in PC therapeutic needs, they discuss in detail the drug discovery goals of these studies, highlighting their translational values for clinically impactful drug nomination.

EXPERT OPINION

Evolving molecular profiling techniques may enable new perspectives for computer-aided approaches to offer drug candidates for different tumor microenvironments. With ongoing efforts to incorporate new compounds into large-scale high-throughput screens, the authors envision continued expansion of drug candidate pools.

Elucidating Compound Mechanism of Action and Polypharmacology with a Large-scale Perturbational Profile Compendium

The Mechanism of Action (MoA) of a drug is generally represented as a small, non-tissue-specific repertoire of high-affinity binding targets. Yet, drug activity and polypharmacology are increasingly associated with a broad range of off-target and tissue-specific effector proteins. To address this challenge, we have leveraged a microfluidics-based Plate-Seq technology to survey drug perturbational profiles representing >700 FDA-approved and experimental oncology drugs, in cell lines selected as high-fidelity models of 23 aggressive tumor subtypes. Built on this dataset, we implemented an efficient computational framework to define a tissue-specific protein activity landscape of these drugs and reported almost 50 million differential protein activities derived from drug perturbations vs. vehicle controls. These analyses revealed thousands of highly reproducible and novel, drug-mediated modulation of tissue-specific targets, leading to generation of a proteome-wide drug functional network, characterization of MoA-related drug clusters and off-target effects, dramatical expansion of druggable human proteome, and identification and experimental validation of novel, tissue-specific inhibitors of undruggable oncoproteins, most never reported before. The drug perturbation profile resource described here represents the first, large-scale, whole-genome-wide, RNA-Seq based dataset assembled to date, with the proposed framework, which is easily extended to elucidating the MoA of novel small-molecule libraries, facilitates mechanistic exploration of drug functions, supports systematic and quantitative approaches to precision oncology, and serves as a rich data foundation for drug discovery.

A systematic review on the culture methods and applications of 3D tumoroids for cancer research and personalized medicine

Cancer is a highly heterogeneous disease, and thus treatment responses vary greatly between patients. To improve therapy efficacy and outcome for cancer patients, more representative and patient-specific preclinical models are needed. Organoids and tumoroids are 3D cell culture models that typically retain the genetic and epigenetic characteristics, as well as the morphology, of their tissue of origin. Thus, they can be used to understand the underlying mechanisms of cancer initiation, progression, and metastasis in a more physiological setting. Additionally, co-culture methods of tumoroids and cancer-associated cells can help to understand the interplay between a tumor and its tumor microenvironment. In recent years, tumoroids have already helped to refine treatments and to identify new targets for cancer therapy. Advanced culturing systems such as chip-based fluidic devices and bioprinting methods in combination with tumoroids have been used for high-throughput applications for personalized medicine. Even though organoid and tumoroid models are complex in vitro systems, validation of results in vivo is still the common practice. Here, we describe how both animal- and human-derived tumoroids have helped to identify novel vulnerabilities for cancer treatment in recent years, and how they are currently used for precision medicine.

Identification and Pharmacological Targeting of Treatment-Resistant, Stem-like Breast Cancer Cells for Combination Therapy

Tumors frequently harbor isogenic yet epigenetically distinct subpopulations of multi-potent cells with high tumor-initiating potential-often called Cancer Stem-Like Cells (CSLCs). These can display preferential resistance to standard-of-care chemotherapy. Single-cell analyses can help elucidate Master Regulator (MR) proteins responsible for governing the transcriptional state of these cells, thus revealing complementary dependencies that may be leveraged via combination therapy. Interrogation of single-cell RNA sequencing profiles from seven metastatic breast cancer patients, using perturbational profiles of clinically relevant drugs, identified drugs predicted to invert the activity of MR proteins governing the transcriptional state of chemoresistant CSLCs, which were then validated by CROP-seq assays. The top drug, the anthelmintic albendazole, depleted this subpopulation in vivo without noticeable cytotoxicity. Moreover, sequential cycles of albendazole and paclitaxel-a commonly used chemotherapeutic -displayed significant synergy in a patient-derived xenograft (PDX) from a TNBC patient, suggesting that network-based approaches can help develop mechanism-based combinatorial therapies targeting complementary subpopulations.

Network-based elucidation of colon cancer drug resistance mechanisms by phosphoproteomic time-series analysis

Aberrant signaling pathway activity is a hallmark of tumorigenesis and progression, which has guided targeted inhibitor design for over 30 years. Yet, adaptive resistance mechanisms, induced by rapid, context-specific signaling network rewiring, continue to challenge therapeutic efficacy. Leveraging progress in proteomic technologies and network-based methodologies, we introduce Virtual Enrichment-based Signaling Protein-activity Analysis (VESPA)-an algorithm designed to elucidate mechanisms of cell response and adaptation to drug perturbations-and use it to analyze 7-point phosphoproteomic time series from colorectal cancer cells treated with clinically-relevant inhibitors and control media. Interrogating tumor-specific enzyme/substrate interactions accurately infers kinase and phosphatase activity, based on their substrate phosphorylation state, effectively accounting for signal crosstalk and sparse phosphoproteome coverage. The analysis elucidates time-dependent signaling pathway response to each drug perturbation and, more importantly, cell adaptive response and rewiring, experimentally confirmed by CRISPR knock-out assays, suggesting broad applicability to cancer and other diseases.

Subtype-selective prenylated isoflavonoids disrupt regulatory drivers of MYCN-amplified cancers

Transcription factors have proven difficult to target with small molecules because they lack pockets necessary for potent binding. Disruption of protein expression can suppress targets and enable therapeutic intervention. To this end, we developed a drug discovery workflow that incorporates cell-line-selective screening and high-throughput expression profiling followed by regulatory network analysis to identify compounds that suppress regulatory drivers of disease. Applying this approach to neuroblastoma (NBL), we screened bioactive molecules in cell lines representing its MYC-dependent (MYCNA) and mesenchymal (MES) subtypes to identify selective compounds, followed by PLATESeq profiling of treated cells. This revealed compounds that disrupt a sub-network of MYCNA-specific regulatory proteins, resulting in MYCN degradation in vivo. The top hit was isopomiferin, a prenylated isoflavonoid that inhibited casein kinase 2 (CK2) in cells. Isopomiferin and its structural analogs inhibited MYC and MYCN in NBL and lung cancer cells, highlighting the general MYC-inhibiting potential of this unique scaffold.

ERG activates a stem-like proliferation-differentiation program in prostate epithelial cells with mixed basal-luminal identity

To gain insight into how ERG translocations cause prostate cancer, we performed single cell transcriptional profiling of an autochthonous mouse model at an early stage of disease initiation. Despite broad expression of ERG in all prostate epithelial cells, proliferation was enriched in a small, stem-like population with mixed-luminal basal identity (called intermediate cells). Through a series of lineage tracing and primary prostate tissue transplantation experiments, we find that tumor initiating activity resides in a subpopulation of basal cells that co-express the luminal genes Tmprss2 and Nkx3.1 (called BasalLum) but not in the larger population of classical Krt8+ luminal cells. Upon ERG activation, BasalLum cells give rise to the highly proliferative intermediate state, which subsequently transitions to the larger population of Krt8+ luminal cells characteristic of ERG-positive human cancers. Furthermore, this proliferative population is characterized by an ERG-specific chromatin state enriched for NFkB, AP-1, STAT and NFAT binding, with implications for TF cooperativity. The fact that the proliferative potential of ERG is enriched in a small stem-like population implicates the chromatin context of these cells as a critical variable for unmasking its oncogenic activity.

Advances in Making Cancer Mouse Models More Accessible and Informative through Non-Germline Genetic Engineering

Genetically engineered mouse models (GEMMs) allow for modeling of spontaneous tumorigenesis within its native microenvironment in mice and have provided invaluable insights into mechanisms of tumorigenesis and therapeutic strategies to treat human disease. However, as their generation requires germline manipulation and extensive animal breeding that is time-, labor-, and cost-intensive, traditional GEMMs are not accessible to most researchers, and fail to model the full breadth of cancer-associated genetic alterations and therapeutic targets. Recent advances in genome-editing technologies and their implementation in somatic tissues of mice have ushered in a new class of mouse models: non-germline GEMMs (nGEMMs). nGEMM approaches can be leveraged to generate somatic tumors de novo harboring virtually any individual or group of genetic alterations found in human cancer in a mouse through simple procedures that do not require breeding, greatly increasing the accessibility and speed and scale on which GEMMs can be produced. Here we describe the technologies and delivery systems used to create nGEMMs and highlight new biological insights derived from these models that have rapidly informed functional cancer genomics, precision medicine, and immune oncology.

Elucidation and Pharmacologic Targeting of Master Regulator Dependencies in Coexisting Diffuse Midline Glioma Subpopulations

Diffuse Midline Gliomas (DMGs) are universally fatal, primarily pediatric malignancies affecting the midline structures of the central nervous system. Despite decades of clinical trials, treatment remains limited to palliative radiation therapy. A major challenge is the coexistence of molecularly distinct malignant cell states with potentially orthogonal drug sensitivities. To address this challenge, we leveraged established network-based methodologies to elucidate Master Regulator (MR) proteins representing mechanistic, non-oncogene dependencies of seven coexisting subpopulations identified by single-cell analysis-whose enrichment in essential genes was validated by pooled CRISPR/Cas9 screens. Perturbational profiles of 372 clinically relevant drugs helped identify those able to invert the activity of subpopulation-specific MRs for follow-up in vivo validation. While individual drugs predicted to target individual subpopulations-including avapritinib, larotrectinib, and ruxolitinib-produced only modest tumor growth reduction in orthotopic models, systemic co-administration induced significant survival extension, making this approach a valuable contribution to the rational design of combination therapy.

Network-Based Methods and Their Applications in Drug Discovery

Drug discovery is time-consuming, expensive, and predominantly follows the "one drug → one target → one disease" paradigm. With the rapid development of systems biology and network pharmacology, a novel drug discovery paradigm, "multidrug → multitarget → multidisease", has emerged. This new holistic paradigm of drug discovery aligns well with the essence of networks, leading to the emergence of network-based methods in the field of drug discovery. In this Perspective, we initially introduce the concept and data sources of networks and highlight classical methodologies employed in network-based methods. Subsequently, we focus on the practical applications of network-based methods across various areas of drug discovery, such as target prediction, virtual screening, prediction of drug therapeutic effects or adverse drug events, and elucidation of molecular mechanisms. In addition, we provide representative web servers for researchers to use network-based methods in specific applications. Finally, we discuss several challenges of network-based methods and the directions for future development. In a word, network-based methods could serve as powerful tools to accelerate drug discovery.

A Transcriptome-Based Precision Oncology Platform for Patient-Therapy Alignment in a Diverse Set of Treatment-Resistant Malignancies

Predicting in vivo response to antineoplastics remains an elusive challenge. We performed a first-of-kind evaluation of two transcriptome-based precision cancer medicine methodologies to predict tumor sensitivity to a comprehensive repertoire of clinically relevant oncology drugs, whose mechanism of action we experimentally assessed in cognate cell lines. We enrolled patients with histologically distinct, poor-prognosis malignancies who had progressed on multiple therapies, and developed low-passage, patient-derived xenograft models that were used to validate 35 patient-specific drug predictions. Both OncoTarget, which identifies high-affinity inhibitors of individual master regulator (MR) proteins, and OncoTreat, which identifies drugs that invert the transcriptional activity of hyperconnected MR modules, produced highly significant 30-day disease control rates (68% and 91%, respectively). Moreover, of 18 OncoTreat-predicted drugs, 15 induced the predicted MR-module activity inversion in vivo. Predicted drugs significantly outperformed antineoplastic drugs selected as unpredicted controls, suggesting these methods may substantively complement existing precision cancer medicine approaches, as also illustrated by a case study.

SIGNIFICANCE

Complementary precision cancer medicine paradigms are needed to broaden the clinical benefit realized through genetic profiling and immunotherapy. In this first-in-class application, we introduce two transcriptome-based tumor-agnostic systems biology tools to predict drug response in vivo. OncoTarget and OncoTreat are scalable for the design of basket and umbrella clinical trials. This article is highlighted in the In This Issue feature, p. 1275.

Systematic Pan-cancer Functional Inference and Validation of Hyper, Hypo and Neomorphic Mutations

While the functional effects of many recurrent cancer mutations have been characterized, the TCGA repository comprises more than 10M non-recurrent events, whose function is unknown. We propose that the context specific activity of transcription factor (TF) proteins-as measured by expression of their transcriptional targets-provides a sensitive and accurate reporter assay to assess the functional role of oncoprotein mutations. Analysis of differentially active TFs in samples harboring mutations of unknown significance-compared to established gain (GOF/hypermorph) or loss (LOF/hypomorph) of function-helped functionally characterize 577,866 individual mutational events across TCGA cohorts, including identification of mutations that are either neomorphic (gain of novel function) or phenocopy other mutations ( mutational mimicry ). Validation using mutation knock-in assays confirmed 15 out of 15 predicted gain and loss of function mutations and 15 of 20 predicted neomorphic mutations. This could help determine targeted therapy in patients with mutations of unknown significance in established oncoproteins.