Search

Childhood neuroblastoma with MYCN-amplification is classified as high-risk and often relapses after intensive treatments. Immune checkpoint blockade therapy against the PD-1/L1 axis shows limited efficacy in neuroblastoma patients and the cancer intrinsic immune regulatory network is poorly understood. Here, we leverage genome-wide CRISPR/Cas9 screens and identify H2AFY as a resistance gene to the clinically approved PD-1 blocking antibody, nivolumab. Analysis of single-cell RNA sequencing datasets reveals that H2AFY mRNA is enriched in adrenergic cancer cells and is associated with worse patient survival. Genetic deletion of H2afy in MYCN-driven neuroblastoma cells reverts in vivo resistance to PD-1 blockade by eliciting activation of the adaptive and innate immunity. Mapping of the epigenetic and translational landscape demonstrates that H2afy deletion promotes cell transition to a mesenchymal-like state. With a multi-omics approach, we uncover H2AFY-associated genes that are functionally relevant and prognostic in patients. Altogether, our study elucidates the role of H2AFY as an epigenetic gatekeeper for cell states and immunogenicity in high-risk neuroblastoma.

The development of functional neurons is a complex orchestration of multiple signaling pathways controlling cell proliferation and differentiation. Because the balance of antioxidants is important for neuronal survival and development, we hypothesized that ferroptosis must be suppressed to gain neurons. We find that removal of antioxidants diminishes neuronal development and laminar organization of cortical organoids, which is fully restored when ferroptosis is inhibited by ferrostatin-1 or when neuronal differentiation occurs in the presence of vitamin A. Furthermore, iron-overload-induced developmental growth defects in C. elegans are ameliorated by vitamin E and A. We determine that all-trans retinoic acid activates the Retinoic Acid Receptor, which orchestrates the expression of anti-ferroptotic genes. In contrast, retinal and retinol show radical-trapping antioxidant activity. Together, our study reveals an unexpected function of vitamin A in coordinating the expression of essential cellular gatekeepers of ferroptosis, and demonstrates that suppression of ferroptosis by radical-trapping antioxidants or by vitamin A is required to obtain mature neurons and proper laminar organization in cortical organoids., (© 2024. The Author(s).)

Immune checkpoint blockade therapy aims to activate the immune system to eliminate cancer cells. However, clinical benefits are only recorded in a subset of patients. Here, we leverage genome-wide CRISPR/Cas9 screens in a Tumor-Immune co-Culture System focusing on triple-negative breast cancer (TNBC). We reveal that NEDD8 loss in cancer cells causes a vulnerability to nivolumab (anti-PD-1). Genetic deletion of NEDD8 only delays cell division initially but cell proliferation is unaffected after recovery. Since the NEDD8 gene is commonly essential, we validate this observation with additional CRISPR screens and uncover enhanced immunogenicity in NEDD8 deficient cells using proteomics. In female immunocompetent mice, PD-1 blockade lacks efficacy against established EO771 breast cancer tumors. In contrast, we observe tumor regression mediated by CD8+ T cells against Nedd8 deficient EO771 tumors after PD-1 blockade. In essence, we provide evidence that NEDD8 is conditionally essential in TNBC and presents as a synergistic drug target for PD-1/L1 blockade therapy., (© 2024. The Author(s).)

Sebaceous glands drive acne, however, their role in other inflammatory skin diseases remains unclear. To shed light on their potential contribution to disease development, we investigated the spatial transcriptome of sebaceous glands in psoriasis and atopic dermatitis patients across lesional and non-lesional human skin samples. Both atopic dermatitis and psoriasis sebaceous glands expressed genes encoding key proteins for lipid metabolism and transport such as ALOX15B, APOC1, FABP7, FADS1/2, FASN, PPARG , and RARRES1. Also, inflammation-related SAA1 was identified as a common spatially variable gene. In atopic dermatitis, genes mainly related to lipid metabolism (e.g. ACAD8, FADS6 , or EBP) as well as disease-specific genes, i.e., Th2 inflammation-related lipid-regulating HSD3B1 were differentially expressed. On the contrary, in psoriasis, more inflammation-related spatially variable genes (e.g. SERPINF1 , FKBP5 , IFIT1/3, DDX58 ) were identified. Other psoriasis-specific enriched pathways included lipid metabolism (e.g. ACOT4, S1PR3) , keratinization (e.g. LCE5A, KRT5/7/16 ), neutrophil degranulation, and antimicrobial peptides (e.g. LTF, DEFB4A, S100A7-9 ). In conclusion, our results show that sebaceous glands contribute to skin homeostasis with a cell type-specific lipid metabolism, which is influenced by the inflammatory microenvironment. These findings further support that sebaceous glands are not bystanders in inflammatory skin diseases, but can actively and differentially modulate inflammation in a disease-specific manner., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2024 Seiringer, Hillig, Schäbitz, Jargosch, Pilz, Eyerich, Szegedi, Sochorová, Gruber, Zouboulis, Biedermann, Menden, Eyerich and Törőcsik.)

Introduction: The concept of Digital Twins (DTs) translated to drug development and clinical trials describes virtual representations of systems of various complexities, ranging from individual cells to entire humans, and enables in silico simulations and experiments. DTs increase the efficiency of drug discovery and development by digitalizing processes associated with high economic, ethical, or social burden. The impact is multifaceted: DT models sharpen disease understanding, support biomarker discovery and accelerate drug development, thus advancing precision medicine. One way to realize DTs is by generative artificial intelligence (AI), a cutting-edge technology that enables the creation of novel, realistic and complex data with desired properties., Areas Covered: The authors provide a brief introduction to generative AI and describe how it facilitates the modeling of DTs. In addition, they compare existing implementations of generative AI for DTs in drug discovery and clinical trials. Finally, they discuss technical and regulatory challenges that should be addressed before DTs can transform drug discovery and clinical trials., Expert Opinion: The current state of DTs in drug discovery and clinical trials does not exploit the entire power of generative AI yet and is limited to simulation of a small number of characteristics. Nonetheless, generative AI has the potential to transform the field by leveraging recent developments in deep learning and customizing models for the needs of scientists, physicians and patients.

Constitutive mucosa-associated lymphoid tissue lymphoma translocation protein 1 (MALT1) activity drives survival of malignant lymphomas addicted to chronic B-cell receptor signaling, oncogenic CARD11, or the API2-MALT1 (also BIRC3::MALT1) fusion oncoprotein. Although MALT1 scaffolding induces NF-κB-dependent survival signaling, MALT1 protease function is thought to augment NF-κB activation by cleaving signaling mediators and transcriptional regulators in B-cell lymphomas. However, the pathological role of MALT1 protease function in lymphomagenesis is not well understood. Here, we show that TRAF6 controls MALT1-dependent activation of NF-κB transcriptional responses but is dispensable for MALT1 protease activation driven by oncogenic CARD11. To uncouple enzymatic and nonenzymatic functions of MALT1, we analyzed TRAF6-dependent and -independent as well as MALT1 protease-dependent gene expression profiles downstream of oncogenic CARD11 and API2-MALT1. The data suggest that by cleaving and inactivating the RNA binding proteins Regnase-1 and Roquin-1/2, MALT1 protease induces posttranscriptional upregulation of many genes including NFKBIZ/IκBζ, NFKBID/IκBNS, and ZC3H12A/Regnase-1 in activated B-cell-like diffuse large B-cell lymphoma (ABC DLBCL). We demonstrate that oncogene-driven MALT1 activity in ABC DLBCL cells regulates NFKBIZ and NFKBID induction on an mRNA level via releasing a brake imposed by Regnase-1 and Roquin-1/2. Furthermore, MALT1 protease drives posttranscriptional gene induction in the context of the API2-MALT1 fusion created by the recurrent t(11;18)(q21;q21) translocation in MALT lymphoma. Thus, MALT1 paracaspase acts as a bifurcation point for enhancing transcriptional and posttranscriptional gene expression in malignant lymphomas. Moreover, the identification of MALT1 protease-selective target genes provides specific biomarkers for the clinical evaluation of MALT1 inhibitors., (© 2023 by The American Society of Hematology. Licensed under Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International (CC BY-NC-ND 4.0), permitting only noncommercial, nonderivative use with attribution. All other rights reserved.)

Mechanistic insight into ageing may empower prolonging the lifespan of humans; however, a complete understanding of this process is still lacking despite a plethora of ageing theories. In order to address this, we investigated the association of lifespan with eight phenotypic traits, that is, litter size, body mass, female and male sexual maturity, somatic mutation, heart, respiratory, and metabolic rate. In support of the somatic mutation theory, we analysed 15 mammalian species and their whole-genome sequencing deriving somatic mutation rate, which displayed the strongest negative correlation with lifespan. All remaining phenotypic traits showed almost equivalent strong associations across this mammalian cohort, however, resting heart rate explained additional variance in lifespan. Integrating somatic mutation and resting heart rate boosted the prediction of lifespan, thus highlighting that resting heart rate may either directly influence lifespan, or represents an epiphenomenon for additional lower-level mechanisms, for example, metabolic rate, that are associated with lifespan., (© 2023 The Authors. Aging Cell published by Anatomical Society and John Wiley & Sons Ltd.)

Precision medicine has revolutionised cancer treatments; however, actionable biomarkers remain scarce. To address this, we develop the Oncology Biomarker Discovery (OncoBird) framework for analysing the molecular and biomarker landscape of randomised controlled clinical trials. OncoBird identifies biomarkers based on single genes or mutually exclusive genetic alterations in isolation or in the context of tumour subtypes, and finally, assesses predictive components by their treatment interactions. Here, we utilise the open-label, randomised phase III trial (FIRE-3, AIO KRK-0306) in metastatic colorectal carcinoma patients, who received either cetuximab or bevacizumab in combination with 5-fluorouracil, folinic acid and irinotecan (FOLFIRI). We systematically identify five biomarkers with predictive components, e.g., patients with tumours that carry chr20q amplifications or lack mutually exclusive ERK signalling mutations benefited from cetuximab compared to bevacizumab. In summary, OncoBird characterises the molecular landscape and outlines actionable biomarkers, which generalises to any molecularly characterised randomised controlled trial., (© 2023. Springer Nature Limited.)

Aberrant DNA methylation accompanies genetic alterations during oncogenesis and tumour homeostasis and contributes to the transcriptional deregulation of key signalling pathways in cancer. Despite increasing efforts in DNA methylation profiling of cancer patients, there is still a lack of epigenetic biomarkers to predict treatment efficacy. To address this, we analyse 721 cancer cell lines across 22 cancer types treated with 453 anti-cancer compounds. We systematically detect the predictive component of DNA methylation in the context of transcriptional and mutational patterns, i.e., in total 19 DNA methylation biomarkers across 17 drugs and five cancer types. DNA methylation constitutes drug sensitivity biomarkers by mediating the expression of proximal genes, thereby enhancing biological signals across multi-omics data modalities. Our method reproduces anticipated associations, and in addition, we find that the NEK9 promoter hypermethylation may confer sensitivity to the NEDD8-activating enzyme (NAE) inhibitor pevonedistat in melanoma through downregulation of NEK9. In summary, we envision that epigenomics will refine existing patient stratification, thus empowering the next generation of precision oncology., (© 2023. Springer Nature Limited.)

The nasal epithelium is an important target for drug delivery to the nose and secondary organs such as the brain via the olfactory bulb. For both topical and brain delivery, the targeting of specific nasal regions such as the olfactory epithelium (brain) is essential, yet challenging. In this study, a numerical model was developed to predict the regional dose as mass per surface area (for an inhaled mass of 2.5 mg), which is the biologically most relevant dose metric for drug delivery in the respiratory system. The role of aerosol diameter (particle diameter: 1 nm to 30 µm) and inhalation flow rate (4, 15 and 30 L/min) in optimal drug delivery to the vestibule, nasal valve, olfactory and nasopharynx is assessed. To obtain the highest doses in the olfactory region, we suggest aerosols with a diameter of 20 µm and a medium inlet air flow rate of 15 L/min. High deposition on the olfactory epithelium was also observed for nanoparticles below 1 nm, as was high residence time (slow flow rate of 4 L/min), but the very low mass of 1 nm nanoparticles is prohibitive for most therapeutic applications. Moreover, high flow rates (30 L/min) and larger micro-aerosols lead to highest doses in the vestibule and nasal valve regions. On the other hand, the highest drug doses in the nasopharynx are observed for nano-aerosol (1 nm) and fine microparticles (1-20 µm) with a relatively weak dependence on flow rate. Furthermore, using the 45 different inhalation scenarios generated by numerical models, different machine learning models with five-fold cross-validation are trained to predict the delivered dose and avoid partial differential equation solvers for future predictions. Random forest and gradient boosting models resulted in R 2 scores of 0.89 and 0.96, respectively. The aerosol diameter and region of interest are the most important features affecting delivered dose, with an approximate importance of 42% and 47%, respectively.

Abundant heterogeneous immune cells infiltrate lesions in chronic inflammatory diseases and characterization of these cells is needed to distinguish disease-promoting from bystander immune cells. Here, we investigate the landscape of non-communicable inflammatory skin diseases (ncISD) by spatial transcriptomics resulting in a large repository of 62,000 spatially defined human cutaneous transcriptomes from 31 patients. Despite the expected immune cell infiltration, we observe rather low numbers of pathogenic disease promoting cytokine transcripts (IFNG, IL13 and IL17A), i.e. >125 times less compared to the mean expression of all other genes over lesional skin sections. Nevertheless, cytokine expression is limited to lesional skin and presented in a disease-specific pattern. Leveraging a density-based spatial clustering method, we identify specific responder gene signatures in direct proximity of cytokines, and confirm that detected cytokine transcripts initiate amplification cascades of up to thousands of specific responder transcripts forming localized epidermal clusters. Thus, within the abundant and heterogeneous infiltrates of ncISD, only a low number of cytokine transcripts and their translated proteins promote disease by initiating an inflammatory amplification cascade in their local microenvironment., (© 2022. The Author(s).)

Pharmacometrics (PM) and machine learning (ML) are both valuable for drug development to characterize pharmacokinetics (PK) and pharmacodynamics (PD). Pharmacokinetic/pharmacodynamic (PKPD) analysis using PM provides mechanistic insight into biological processes but is time- and labor-intensive. In contrast, ML models are much quicker trained, but offer less mechanistic insights. The opportunity of using ML predictions of drug PK as input for a PKPD model could strongly accelerate analysis efforts. Here exemplified by rifampicin, a widely used antibiotic, we explore the ability of different ML algorithms to predict drug PK. Based on simulated data, we trained linear regressions (LASSO), Gradient Boosting Machines, XGBoost and Random Forest to predict the plasma concentration-time series and rifampicin area under the concentration-versus-time curve from 0-24 h (AUC 0-24h ) after repeated dosing. XGBoost performed best for prediction of the entire PK series ( R 2 : 0.84, root mean square error (RMSE): 6.9 mg/L, mean absolute error (MAE): 4.0 mg/L) for the scenario with the largest data size. For AUC 0-24h prediction, LASSO showed the highest performance ( R 2 : 0.97, RMSE: 29.1 h·mg/L, MAE: 18.8 h·mg/L). Increasing the number of plasma concentrations per patient (0, 2 or 6 concentrations per occasion) improved model performance. For example, for AUC 0-24h prediction using LASSO, the R 2 was 0.41, 0.69 and 0.97 when using predictors only (no plasma concentrations), 2 or 6 plasma concentrations per occasion as input, respectively. Run times for the ML models ranged from 1.0 s to 8 min, while the run time for the PM model was more than 3 h. Furthermore, building a PM model is more time- and labor-intensive compared with ML. ML predictions of drug PK could thus be used as input into a PKPD model, enabling time-efficient analysis.

Hospital staff are at high risk for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection during the coronavirus disease (COVID-19) pandemic. This cross-sectional study aimed to determine the prevalence of SARS-CoV-2 infection in hospital staff at the University Hospital rechts der Isar in Munich, Germany, and identify modulating factors. Overall seroprevalence of SARS-CoV-2-IgG in 4,554 participants was 2.4%. Staff engaged in direct patient care, including those working in COVID-19 units, had a similar probability of being seropositive as non-patient-facing staff. Increased probability of infection was observed in staff reporting interactions with SARS-CoV-2‒infected coworkers or private contacts or exposure to COVID-19 patients without appropriate personal protective equipment. Analysis of spatiotemporal trajectories identified that distinct hotspots for SARS-CoV-2‒positive staff and patients only partially overlap. Patient-facing work in a healthcare facility during the SARS-CoV-2 pandemic might be safe as long as adequate personal protective equipment is used and infection prevention practices are followed inside and outside the hospital.

Background: Parkinson's disease (PD) is the second most common neurodegenerative disorder whose prevalence is rapidly increasing worldwide. The molecular mechanisms underpinning the pathophysiology of sporadic PD remain incompletely understood. Therefore, causative therapies are still elusive. To obtain a more integrative view of disease-mediated alterations, we investigated the molecular landscape of PD in human post-mortem midbrains, a region that is highly affected during the disease process., Methods: Tissue from 19 PD patients and 12 controls were obtained from the Parkinson's UK Brain Bank and subjected to multi-omic analyses: small and total RNA sequencing was performed on an Illumina's HiSeq4000, while proteomics experiments were performed in a hybrid triple quadrupole-time of flight mass spectrometer (TripleTOF5600+) following quantitative sequential window acquisition of all theoretical mass spectra. Differential expression analyses were performed with customized frameworks based on DESeq2 (for RNA sequencing) and with Perseus v.1.5.6.0 (for proteomics). Custom pipelines in R were used for integrative studies., Results: Our analyses revealed multiple deregulated molecular targets linked to known disease mechanisms in PD as well as to novel processes. We have identified and experimentally validated (quantitative real-time polymerase chain reaction/western blotting) several PD-deregulated molecular candidates, including miR-539-3p, miR-376a-5p, miR-218-5p and miR-369-3p, the valid miRNA-mRNA interacting pairs miR-218-5p/RAB6C and miR-369-3p/GTF2H3, as well as multiple proteins, such as CHI3L1, HSPA1B, FNIP2 and TH. Vertical integration of multi-omic analyses allowed validating disease-mediated alterations across different molecular layers. Next to the identification of individual molecular targets in all explored omics layers, functional annotation of differentially expressed molecules showed an enrichment of pathways related to neuroinflammation, mitochondrial dysfunction and defects in synaptic function., Conclusions: This comprehensive assessment of PD-affected and control human midbrains revealed multiple molecular targets and networks that are relevant to the disease mechanism of advanced PD. The integrative analyses of multiple omics layers underscore the importance of neuroinflammation, immune response activation, mitochondrial and synaptic dysfunction as putative therapeutic targets for advanced PD., (© 2022 The Authors. Clinical and Translational Medicine published by John Wiley & Sons Australia, Ltd on behalf of Shanghai Institute of Clinical Bioinformatics.)

The nasal olfactory region is a potential route for non-invasive delivery of drugs directly from the nasal epithelium to the brain, bypassing the often impermeable blood-brain barrier. However, efficient aerosol delivery to the olfactory region is challenging due to its location in the nose. Here we explore aerosol delivery with bi-directional pulsatile flow conditions for targeted drug delivery to the olfactory region using a computational fluid dynamics (CFD) model on the patient-specific nasal geometry. Aerosols with aerodynamic diameter of 1 µm, which is large enough for delivery of large enough drug doses and yet potentially small enough for non-inertial aerosol deposition due to, e.g., particle diffusion and flow oscillations, is inhaled for 1.98 s through one nostril and exhaled through the other one. The bi-directional aerosol delivery with steady flow rate of 4 L/min results in deposition efficiencies (DEs) of 50.9 and 0.48% in the nasal cavity and olfactory region, respectively. Pulsatile flow with average flow rate of 4 L/min (frequency: 45 Hz) reduces these values to 34.4 and 0.12%, respectively, and it mitigates the non-uniformity of right-left deposition in both the cavity (from 1.77- to 1.33-fold) and the olfactory region (from 624- to 53.2-fold). The average drug dose deposited in the nasal cavity and the olfactory epithelium region is very similar in the right nasal cavity independent of pulsation conditions (inhalation side). In contrast, the local aerosol dose in the olfactory region of the left side is at least 100-fold lower than that in the nasal cavity independent of pulsation condition. Hence, while pulsatile flow reduces the right-left (inhalation-exhalation) imbalance, it is not able to overcome it. However, the inhalation side (even with pulsation) allows for relatively high olfactory epithelium drug doses per area reaching the same level as in the total nasal cavity. Due to the relatively low drug deposition in olfactory region on the exhalation side, this allows either very efficient targeting of the inhalation side, or uniform drug delivery by performing bidirectional flow first from the one and then from the other side of the nose., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2021 Farnoud, Tofighian, Baumann, Martin, Rashidi, Menden and Schmid.)

Disparities between risk, treatment outcomes and survival rates in cancer patients across the world may be attributed to socioeconomic factors. In addition, the role of ancestry is frequently discussed. In preclinical studies, high-throughput drug screens in cancer cell lines have empowered the identification of clinically relevant molecular biomarkers of drug sensitivity; however, the genetic ancestry from tissue donors has been largely neglected in this setting. In order to address this, here, we show that the inferred ancestry of cancer cell lines is conserved and may impact drug response in patients as a predictive covariate in high-throughput drug screens. We found that there are differential drug responses between European and East Asian ancestries, especially when treated with PI3K/mTOR inhibitors. Our finding emphasizes a new angle in precision medicine, as cancer intervention strategies should consider the germline landscape, thereby reducing the failure rate of clinical trials.

Introduction : Precision medicine is the concept of treating diseases based on environmental factors, lifestyles, and molecular profiles of patients. This approach has been found to increase success rates of clinical trials and accelerate drug approvals. However, current precision medicine applications in early drug discovery use only a handful of molecular biomarkers to make decisions, whilst clinics gear up to capture the full molecular landscape of patients in the near future. This deep multi-omics characterization demands new analysis strategies to identify appropriate treatment regimens, which we envision will be pioneered by artificial intelligence. Areas covered : In this review, the authors discuss the current state of drug discovery in precision medicine and present our vision of how artificial intelligence will impact biomarker discovery and drug design. Expert opinion : Precision medicine is expected to revolutionize modern medicine; however, its traditional form is focusing on a few biomarkers, thus not equipped to leverage the full power of molecular landscapes. For learning how the development of drugs can be tailored to the heterogeneity of patients across their molecular profiles, artificial intelligence algorithms are the next frontier in precision medicine and will enable a fully personalized approach in drug design, and thus ultimately impacting clinical practice.

High-throughput testing of drugs across molecular-characterised cell lines can identify candidate treatments and discover biomarkers. However, the cells' response to a drug is typically quantified by a summary statistic from a best-fit dose-response curve, whilst neglecting the uncertainty of the curve fit and the potential variability in the raw readouts. Here, we model the experimental variance using Gaussian Processes, and subsequently, leverage uncertainty estimates to identify associated biomarkers with a new Bayesian framework. Applied to in vitro screening data on 265 compounds across 1074 cancer cell lines, our models identified 24 clinically established drug-response biomarkers, and provided evidence for six novel biomarkers by accounting for association with low uncertainty. We validated our uncertainty estimates with an additional drug screen of 26 drugs, 10 cell lines with 8 to 9 replicates. Our method is applicable to any dose-response data without replicates, and improves biomarker discovery for precision medicine., Competing Interests: DW, GK, TT, AG, JS, MG, MM No competing interests declared, JH James Hensman is an employee of Amazon.com. The author has no competing financial interests to declare. JD Jonathan Dry is affiliated with AstraZeneca and Tempus. The author has no competing financial interests to declare. FD Frank Dondelinger is an employee of Roche. The author has no competing financial interests to declare., (© 2020, Wang et al.)

High-throughput drug screens in cancer cell lines test compounds at low concentrations, thereby enabling the identification of drug-sensitivity biomarkers, while resistance biomarkers remain underexplored. Dissecting meaningful drug responses at high concentrations is challenging due to cytotoxicity, i.e., off-target effects, thus limiting resistance biomarker discovery to frequently mutated cancer genes. To address this, we interrogate subpopulations carrying sensitivity biomarkers and consecutively investigate unexpectedly resistant (UNRES) cell lines for unique genetic alterations that may drive resistance. By analyzing the GDSC and CTRP datasets, we find 53 and 35 UNRES cases, respectively. For 24 and 28 of them, we highlight putative resistance biomarkers. We find clinically relevant cases such as EGFR T790M mutation in NCI-H1975 or PTEN loss in NCI-H1650 cells, in lung adenocarcinoma treated with EGFR inhibitors. Interrogating the underpinnings of drug resistance with publicly available CRISPR phenotypic assays assists in prioritizing resistance drivers, offering hypotheses for drug combinations., Competing Interests: F.I. receives funding from Open Targets, a public-private initiative involving academia and industry, and performs consultancy for the joint 10.13039/501100000289CRUK-10.13039/100004325AstraZeneca Functional Genomics Center. B.S., J.R.D., and U.M. are employees of AstraZeneca., (© 2020 The Authors.)

In silico models to predict which tumors will respond to a given drug are necessary for Precision Oncology. However, predictive models are only available for a handful of cases (each case being a given drug acting on tumors of a specific cancer type). A way to generate predictive models for the remaining cases is with suitable machine learning algorithms that are yet to be applied to existing in vitro pharmacogenomics datasets. Here, we apply XGBoost integrated with a stringent feature selection approach, which is an algorithm that is advantageous for these high-dimensional problems. Thus, we identified and validated 118 predictive models for 62 drugs across five cancer types by exploiting four molecular profiles (sequence mutations, copy-number alterations, gene expression, and DNA methylation). Predictive models were found in each cancer type and with every molecular profile. On average, no omics profile or cancer type obtained models with higher predictive accuracy than the rest. However, within a given cancer type, some molecular profiles were overrepresented among predictive models. For instance, CNA profiles were predictive in breast invasive carcinoma (BRCA) cell lines, but not in small cell lung cancer (SCLC) cell lines where gene expression (GEX) and DNA methylation profiles were the most predictive. Lastly, we identified the best XGBoost model per cancer type and analyzed their selected features. For each model, some of the genes in the selected list had already been found to be individually linked to the response to that drug, providing additional evidence of the usefulness of these models and the merits of the feature selection scheme.

Drug combinations can expand therapeutic options and address cancer's resistance. However, the combinatorial space is enormous precluding its systematic exploration. Therefore, synergy prediction strategies are essential. We here present an approach to prioritise drug combinations in high-throughput screens and to stratify synergistic responses. At the core of our approach is the observation that the likelihood of synergy increases when targeting proteins with either strong functional similarity or dissimilarity. We estimate the similarity applying a multitask machine learning approach to basal gene expression and response to single drugs. We tested 7 protein target pairs (representing 29 combinations) and predicted their synergies in 33 breast cancer cell lines. In addition, we experimentally validated predicted synergy of the BRAF/insulin receptor combination (Dabrafenib/BMS-754807) in 48 colorectal cancer cell lines. We anticipate that our approaches can be used for prioritization of drug combinations in large scale screenings, and to maximize the efficacy of drugs already known to induce synergy, ultimately enabling patient stratification.

Synergistic drug combinations are commonly sought to overcome monotherapy resistance in cancer treatment. To identify such combinations, high-throughput cancer cell line combination screens are performed; and synergy is quantified using competing models based on fundamentally different assumptions. Here, we compare the behaviour of four synergy models, namely Loewe additivity, Bliss independence, highest single agent and zero interaction potency, using the Merck oncology combination screen. We evaluate agreements and disagreements between the models and investigate putative artefacts of each model's assumptions. Despite at least moderate concordance between scores (Pearson's r >0.32, Spearman's ρ>0.34), multiple instances of strong disagreement were observed. Those disagreements are driven by, among others, large differences in tested concentrations, maximum response values and median effective concentrations., (Copyright © 2019 Elsevier Ltd. All rights reserved.)

Here we describe a proteomic data resource for the NCI-60 cell lines generated by pressure cycling technology and SWATH mass spectrometry. We developed the DIA-expert software to curate and visualize the SWATH data, leading to reproducible detection of over 3,100 SwissProt proteotypic proteins and systematic quantification of pathway activities. Stoichiometric relationships of interacting proteins for DNA replication, repair, the chromatin remodeling NuRD complex, β-catenin, RNA metabolism, and prefoldins are more evident than that at the mRNA level. The data are available in CellMiner (discover.nci.nih.gov/cellminercdb and discover.nci.nih.gov/cellminer), allowing casual users to test hypotheses and perform integrative, cross-database analyses of multi-omic drug response correlations for over 20,000 drugs. We demonstrate the value of proteome data in predicting drug response for over 240 clinically relevant chemotherapeutic and targeted therapies. In summary, we present a novel proteome resource for the NCI-60, together with relevant software tools, and demonstrate the benefit of proteome analyses., (Copyright © 2019 The Author(s). Published by Elsevier Inc. All rights reserved.)

Personalised medicine has predominantly focused on genetically altered cancer genes that stratify drug responses, but there is a need to objectively evaluate differential pharmacology patterns at a subpopulation level. Here, we introduce an approach based on unsupervised machine learning to compare the pharmacological response relationships between 327 pairs of cancer therapies. This approach integrated multiple measures of response to identify subpopulations that react differently to inhibitors of the same or different targets to understand mechanisms of resistance and pathway cross-talk. MEK, BRAF, and PI3K inhibitors were shown to be effective as combination therapies for particular BRAF mutant subpopulations. A systematic analysis of preclinical data for a failed phase III trial of selumetinib combined with docetaxel in lung cancer suggests potential indications in pancreatic and colorectal cancers with KRAS mutation. This data-informed study exemplifies a method for stratified medicine to identify novel cancer subpopulations, their genetic biomarkers, and effective drug combinations., Competing Interests: Competing interestsNK is an employee of Constellation Analytics, LLC. The remaining authors declare no competing interests., (© The Author(s) 2019.)

The effectiveness of most cancer targeted therapies is short-lived. Tumors often develop resistance that might be overcome with drug combinations. However, the number of possible combinations is vast, necessitating data-driven approaches to find optimal patient-specific treatments. Here we report AstraZeneca's large drug combination dataset, consisting of 11,576 experiments from 910 combinations across 85 molecularly characterized cancer cell lines, and results of a DREAM Challenge to evaluate computational strategies for predicting synergistic drug pairs and biomarkers. 160 teams participated to provide a comprehensive methodological development and benchmarking. Winning methods incorporate prior knowledge of drug-target interactions. Synergy is predicted with an accuracy matching biological replicates for >60% of combinations. However, 20% of drug combinations are poorly predicted by all methods. Genomic rationale for synergy predictions are identified, including ADAM17 inhibitor antagonism when combined with PIK3CB/D inhibition contrasting to synergy when combined with other PI3K-pathway inhibitors in PIK3CA mutant cells.

Patients with seemingly the same tumour can respond very differently to treatment. There are strong, well-established effects of somatic mutations on drug efficacy, but there is at-most anecdotal evidence of a germline component to drug response. Here, we report a systematic survey of how inherited germline variants affect drug susceptibility in cancer cell lines. We develop a joint analysis approach that leverages both germline and somatic variants, before applying it to screening data from 993 cell lines and 265 drugs. Surprisingly, we find that the germline contribution to variation in drug susceptibility can be as large or larger than effects due to somatic mutations. Several of the associations identified have a direct relationship to the drug target. Finally, using 17-AAG response as an example, we show how germline effects in combination with transcriptomic data can be leveraged for improved patient stratification and to identify new markers for drug sensitivity.

Motivation: Large pharmacogenomic screenings integrate heterogeneous cancer genomic datasets as well as anti-cancer drug responses on thousand human cancer cell lines. Mining this data to identify new therapies for cancer sub-populations would benefit from common data structures, modular computational biology tools and user-friendly interfaces., Results: We have developed GDSCTools: a software aimed at the identification of clinically relevant genomic markers of drug response. The Genomics of Drug Sensitivity in Cancer (GDSC) database (www.cancerRxgene.org) integrates heterogeneous cancer genomic datasets as well as anti-cancer drug responses on a thousand cancer cell lines. Including statistical tools (analysis of variance) and predictive methods (Elastic Net), as well as common data structures, GDSCTools allows users to reproduce published results from GDSC and to implement new analytical methods. In addition, non-GDSC data resources can also be analysed since drug responses and genomic features can be encoded as CSV files., Contact: thomas.cokelaer@pasteur.fr or saezrodriguez.rwth-aachen.de or mg12@sanger.ac.uk., Supplementary Information: Supplementary data are available at Bioinformatics online.

Systematic studies of cancer genomes have provided unprecedented insights into the molecular nature of cancer. Using this information to guide the development and application of therapies in the clinic is challenging. Here, we report how cancer-driven alterations identified in 11,289 tumors from 29 tissues (integrating somatic mutations, copy number alterations, DNA methylation, and gene expression) can be mapped onto 1,001 molecularly annotated human cancer cell lines and correlated with sensitivity to 265 drugs. We find that cell lines faithfully recapitulate oncogenic alterations identified in tumors, find that many of these associate with drug sensitivity/resistance, and highlight the importance of tissue lineage in mediating drug response. Logic-based modeling uncovers combinations of alterations that sensitize to drugs, while machine learning demonstrates the relative importance of different data types in predicting drug response. Our analysis and datasets are rich resources to link genotypes with cellular phenotypes and to identify therapeutic options for selected cancer sub-populations., (Copyright © 2016 The Author(s). Published by Elsevier Inc. All rights reserved.)

Unlabelled: DREAM challenges are community competitions designed to advance computational methods and address fundamental questions in system biology and translational medicine. Each challenge asks participants to develop and apply computational methods to either predict unobserved outcomes or to identify unknown model parameters given a set of training data. Computational methods are evaluated using an automated scoring metric, scores are posted to a public leaderboard, and methods are published to facilitate community discussions on how to build improved methods. By engaging participants from a wide range of science and engineering backgrounds, DREAM challenges can comparatively evaluate a wide range of statistical, machine learning, and biophysical methods. Here, we describe DREAMTools, a Python package for evaluating DREAM challenge scoring metrics. DREAMTools provides a command line interface that enables researchers to test new methods on past challenges, as well as a framework for scoring new challenges. As of September 2015, DREAMTools includes more than 80% of completed DREAM challenges. DREAMTools complements the data, metadata, and software tools available at the DREAM website http://dreamchallenges.org and on the Synapse platform https://www.synapse.org., Availability: DREAMTools is a Python package. Releases and documentation are available at http://pypi.python.org/pypi/dreamtools. The source code is available at http://github.com/dreamtools.

The ability to computationally predict the effects of toxic compounds on humans could help address the deficiencies of current chemical safety testing. Here, we report the results from a community-based DREAM challenge to predict toxicities of environmental compounds with potential adverse health effects for human populations. We measured the cytotoxicity of 156 compounds in 884 lymphoblastoid cell lines for which genotype and transcriptional data are available as part of the Tox21 1000 Genomes Project. The challenge participants developed algorithms to predict interindividual variability of toxic response from genomic profiles and population-level cytotoxicity data from structural attributes of the compounds. 179 submitted predictions were evaluated against an experimental data set to which participants were blinded. Individual cytotoxicity predictions were better than random, with modest correlations (Pearson's r < 0.28), consistent with complex trait genomic prediction. In contrast, predictions of population-level response to different compounds were higher (r < 0.66). The results highlight the possibility of predicting health risks associated with unknown compounds, although risk estimation accuracy remains suboptimal.

Recent therapeutic successes have renewed interest in drug combinations, but experimental screening approaches are costly and often identify only small numbers of synergistic combinations. The DREAM consortium launched an open challenge to foster the development of in silico methods to computationally rank 91 compound pairs, from the most synergistic to the most antagonistic, based on gene-expression profiles of human B cells treated with individual compounds at multiple time points and concentrations. Using scoring metrics based on experimental dose-response curves, we assessed 32 methods (31 community-generated approaches and SynGen), four of which performed significantly better than random guessing. We highlight similarities between the methods. Although the accuracy of predictions was not optimal, we find that computational prediction of compound-pair activity is possible, and that community challenges can be useful to advance the field of in silico compound-synergy prediction.

Predicting the best treatment strategy from genomic information is a core goal of precision medicine. Here we focus on predicting drug response based on a cohort of genomic, epigenomic and proteomic profiling data sets measured in human breast cancer cell lines. Through a collaborative effort between the National Cancer Institute (NCI) and the Dialogue on Reverse Engineering Assessment and Methods (DREAM) project, we analyzed a total of 44 drug sensitivity prediction algorithms. The top-performing approaches modeled nonlinear relationships and incorporated biological pathway information. We found that gene expression microarrays consistently provided the best predictive power of the individual profiling data sets; however, performance was increased by including multiple, independent data sets. We discuss the innovations underlying the top-performing methodology, Bayesian multitask MKL, and we provide detailed descriptions of all methods. This study establishes benchmarks for drug sensitivity prediction and identifies approaches that can be leveraged for the development of new methods.

Predicting the response of a specific cancer to a therapy is a major goal in modern oncology that should ultimately lead to a personalised treatment. High-throughput screenings of potentially active compounds against a panel of genomically heterogeneous cancer cell lines have unveiled multiple relationships between genomic alterations and drug responses. Various computational approaches have been proposed to predict sensitivity based on genomic features, while others have used the chemical properties of the drugs to ascertain their effect. In an effort to integrate these complementary approaches, we developed machine learning models to predict the response of cancer cell lines to drug treatment, quantified through IC₅₀ values, based on both the genomic features of the cell lines and the chemical properties of the considered drugs. Models predicted IC₅₀ values in a 8-fold cross-validation and an independent blind test with coefficient of determination R² of 0.72 and 0.64 respectively. Furthermore, models were able to predict with comparable accuracy (R² of 0.61) IC50s of cell lines from a tissue not used in the training stage. Our in silico models can be used to optimise the experimental design of drug-cell screenings by estimating a large proportion of missing IC₅₀ values rather than experimentally measuring them. The implications of our results go beyond virtual drug screening design: potentially thousands of drugs could be probed in silico to systematically test their potential efficacy as anti-tumour agents based on their structure, thus providing a computational framework to identify new drug repositioning opportunities as well as ultimately be useful for personalized medicine by linking the genomic traits of patients to drug sensitivity.

Whereas genomic data are universally machine-readable, data from imaging, multiplex biochemistry, flow cytometry and other cell- and tissue-based assays usually reside in loosely organized files of poorly documented provenance. This arises because the relational databases used in genomic research are difficult to adapt to rapidly evolving experimental designs, data formats and analytic algorithms. Here we describe an adaptive approach to managing experimental data based on semantically typed data hypercubes (SDCubes) that combine hierarchical data format 5 (HDF5) and extensible markup language (XML) file types. We demonstrate the application of SDCube-based storage using ImageRail, a software package for high-throughput microscopy. Experimental design and its day-to-day evolution, not rigid standards, determine how ImageRail data are organized in SDCubes. We applied ImageRail to collect and analyze drug dose-response landscapes in human cell lines at single-cell resolution.