Deep thermal profiling for detection of functional proteoform groups.

The complexity of the functional proteome extends considerably beyond the coding genome, resulting in millions of proteoforms. Investigation of proteoforms and their functional roles is important to understand cellular physiology and its deregulation in diseases but challenging to perform systematically. Here we applied thermal proteome profiling with deep peptide coverage to detect functional proteoform groups in acute lymphoblastic leukemia cell lines with different cytogenetic aberrations. We detected 15,846 proteoforms, capturing differently spliced, cleaved and post-translationally modified proteins expressed from 9,290 genes. We identified differential co-aggregation of proteoform pairs and established links to disease biology. Moreover, we systematically made use of measured biophysical proteoform states to find specific biomarkers of drug sensitivity. Our approach, thus, provides a powerful and unique tool for systematic detection and functional annotation of proteoform groups.

Thermal proteome profiling identifies PIP4K2A and ZADH2 as off-targets of Polo-like kinase 1 inhibitor volasertib.

Polo-like kinase 1 (PLK1) is an important cell cycle kinase and an attractive target for anticancer treatments. An ATP-competitive small molecular PLK1 inhibitor, volasertib, has reached phase III in clinical trials in patients with refractory acute myeloid leukemia as a combination treatment with cytarabine. However, severe side effects limited its use. The origin of the side effects is unclear and might be due to insufficient specificity of the drug. Thus, identifying potential off-targets to volasertib is important for future clinical trials and for the development of more specific drugs. In this study, we used thermal proteome profiling (TPP) to identify proteome-wide targets of volasertib. Apart from PLK1 and proteins regulated by PLK1, we identified about 200 potential volasertib off-targets. Comparison of this result with the mass-spectrometry analysis of volasertib-treated cells showed that phosphatidylinositol phosphate and prostaglandin metabolism pathways are affected by volasertib. We confirmed that PIP4K2A and ZADH2-marker proteins for these pathways-are, indeed, stabilized by volasertib. PIP4K2A, however, was not affected by another PLK1 inhibitor onvansertib, suggesting that PIP4K2A is a true off-target of volasertib. Inhibition of these proteins is known to impact both the immune response and fatty acid metabolism and could explain some of the side effects seen in volasertib-treated patients.

The transcriptome-wide landscape of molecular subtype-specific mRNA expression profiles in acute myeloid leukemia.

Molecular classification of acute myeloid leukemia (AML) aids prognostic stratification and clinical management. Our aim in this study is to identify transcriptome-wide mRNAs that are specific to each of the molecular subtypes of AML. We analyzed RNA-sequencing data of 955 AML samples from three cohorts, including the BeatAML project, the Cancer Genome Atlas, and a cohort of Swedish patients to provide a comprehensive transcriptome-wide view of subtype-specific mRNA expression. We identified 729 subtype-specific mRNAs, discovered in the BeatAML project and validated in the other two cohorts. Using unique proteomics data, we also validated the presence of subtype-specific mRNAs at the protein level, yielding a rich collection of potential protein-based biomarkers for the AML community. To enable the exploration of subtype-specific mRNA expression by the broader scientific community, we provide an interactive resource to the public.

Delineating functional and molecular impact of ex vivo sample handling in precision medicine.

Consistent handling of samples is crucial for achieving reproducible molecular and functional testing results in translational research. Here, we used 229 acute myeloid leukemia (AML) patient samples to assess the impact of sample handling on high-throughput functional drug testing, mass spectrometry-based proteomics, and flow cytometry. Our data revealed novel and previously described changes in cell phenotype and drug response dependent on sample biobanking. Specifically, myeloid cells with a CD117 (c-KIT) positive phenotype decreased after biobanking, potentially distorting cell population representations and affecting drugs targeting these cells. Additionally, highly granular AML cell numbers decreased after freezing. Secondly, protein expression levels, as well as sensitivity to drugs targeting cell proliferation, metabolism, tyrosine kinases (e.g., JAK, KIT, FLT3), and BH3 mimetics were notably affected by biobanking. Moreover, drug response profiles of paired fresh and frozen samples showed that freezing samples can lead to systematic errors in drug sensitivity scores. While a high correlation between fresh and frozen for the entire drug library was observed, freezing cells had a considerable impact at an individual level, which could influence outcomes in translational studies. Our study highlights conditions where standardization is needed to improve reproducibility, and where validation of data generated from biobanked cohorts may be particularly important.

Targeting autophagy as a therapeutic strategy in pediatric acute lymphoblastic leukemia.

Autophagy is activated in response to a variety of stress conditions including anti-cancer therapies, and tumors cells often depend on autophagy for survival. In this study, we have evaluated inhibition of autophagy as therapeutic strategy in acute lymphoblastic leukemia (ALL) in children, both as a single treatment and in combination with glucocorticoid (GC) Dexamethasone (Dexa). Analysis of proteomics and RNA-seq of ALL cell lines and primary samples identified an upregulation of Vps34 and ATG14 proteins and autophagy and lysosomal pathway enrichment in a genetic subgroup with a recurrent t(12;21) translocation. Cells from this sugbroup were also significantly more sensitive to the selective autophagy or lysosomal inhibitors than cells with other genetic rearrangements. Further, combination of Dexa with either lysosomal or autophagy inhibitors was either synergistic or additive in killing leukemic cells across various genetic and lineage backgrounds, for both cell lines and primary samples, as assessed using viability assays and SynergyFinder as well as apoptotic caspase 3/7-based live-cell assays. Our data demonstrate that targeting autophagy represents a promising strategy for the treatment of pediatric ALL, both as a selective modality for the t(12;21) pre-B-ALL subgroup, and in combination treatments to sensitize to GC-induced cytotoxicity.

FORALL: an interactive shiny/R web portal to navigate multi-omics high-throughput data of pediatric acute lymphoblastic leukemia.

Motivation: Pediatric acute lymphoblastic leukemia (ALL) is the most common cancer among children worldwide. The availability of easily accessible multi-omics data provides unprecedented resources and opportunities for discovering and refining disease biology, cancer biomarkers, and drug mechanisms of action. This has led to exponential growth of omics data available in public repositories. However, delivering the useful information and knowledge extraction from this data is one of the bottlenecks of multi-omics. Presenting, navigating, and downloading ALL omics data in a user-friendly interface provide a valuable platform for biologists and clinicians to get most of the omics data. Our in-house data provides in-depth mass spectrometry-based protein abundance data for a large panel of commercially available ALL cell lines. Providing this data to the scientific community in the form of a user-friendly web-portal allows for easy and detailed exploration of the data. Results: We have developed the Functional Omics Resource of Acute Lymphoblastic Leukemia (FORALL) web-portal. FORALL is a shiny-based web portal designed to navigate in-depth mass spectrometry-based proteomics data of 51 cell lines. The proteomics data can be navigated and visualized along with matched RNA expression data as well as drug sensitivity data of 528 investigational and approved drugs. Availability and implementation: FORALL is available at https://proteomics.se/forall/.

Prediction model for drug response of acute myeloid leukemia patients.

Despite some encouraging successes, predicting the therapy response of acute myeloid leukemia (AML) patients remains highly challenging due to tumor heterogeneity. Here we aim to develop and validate MDREAM, a robust ensemble-based prediction model for drug response in AML based on an integration of omics data, including mutations and gene expression, and large-scale drug testing. Briefly, MDREAM is first trained in the BeatAML cohort (n = 278), and then validated in the BeatAML (n = 183) and two external cohorts, including a Swedish AML cohort (n = 45) and a relapsed/refractory acute leukemia cohort (n = 12). The final prediction is based on 122 ensemble models, each corresponding to a drug. A confidence score metric is used to convey the uncertainty of predictions; among predictions with a confidence score >0.75, the validated proportion of good responders is 77%. The Spearman correlations between the predicted and the observed drug response are 0.68 (95% CI: [0.64, 0.68]) in the BeatAML validation set, -0.49 (95% CI: [-0.53, -0.44]) in the Swedish cohort and 0.59 (95% CI: [0.51, 0.67]) in the relapsed/refractory cohort. A web-based implementation of MDREAM is publicly available at https://www.meb.ki.se/shiny/truvu/MDREAM/ .

Preclinical efficacy of azacitidine and venetoclax for infant KMT2A-rearranged acute lymphoblastic leukemia reveals a new therapeutic strategy.

Infants with KMT2A-rearranged B-cell acute lymphoblastic leukemia (ALL) have a dismal prognosis. Survival outcomes have remained static in recent decades despite treatment intensification and novel therapies are urgently required. KMT2A-rearranged infant ALL cells are characterized by an abundance of promoter hypermethylation and exhibit high BCL-2 expression, highlighting potential for therapeutic targeting. Here, we show that hypomethylating agents exhibit in vitro additivity when combined with most conventional chemotherapeutic agents. However, in a subset of samples an antagonistic effect was seen between several agents. This was most evident when hypomethylating agents were combined with methotrexate, with upregulation of ATP-binding cassette transporters identified as a potential mechanism. Single agent treatment with azacitidine and decitabine significantly prolonged in vivo survival in KMT2A-rearranged infant ALL xenografts. Treatment of KMT2A-rearranged infant ALL cell lines with azacitidine and decitabine led to differential genome-wide DNA methylation, changes in gene expression and thermal proteome profiling revealed the target protein-binding landscape of these agents. The selective BCL-2 inhibitor, venetoclax, exhibited in vitro additivity in combination with hypomethylating or conventional chemotherapeutic agents. The addition of venetoclax to azacitidine resulted in a significant in vivo survival advantage indicating the therapeutic potential of this combination to improve outcome for infants with KMT2A-rearranged ALL.

Integrative multi-omics and drug response profiling of childhood acute lymphoblastic leukemia cell lines.

Acute lymphoblastic leukemia (ALL) is the most common childhood cancer. Although standard-of-care chemotherapeutics are sufficient for most ALL cases, there are subsets of patients with poor response who relapse in disease. The biology underlying differences between subtypes and their response to therapy has only partially been explained by genetic and transcriptomic profiling. Here, we perform comprehensive multi-omic analyses of 49 readily available childhood ALL cell lines, using proteomics, transcriptomics, and pharmacoproteomic characterization. We connect the molecular phenotypes with drug responses to 528 oncology drugs, identifying drug correlations as well as lineage-dependent correlations. We also identify the diacylglycerol-analog bryostatin-1 as a therapeutic candidate in the MEF2D-HNRNPUL1 fusion high-risk subtype, for which this drug activates pro-apoptotic ERK signaling associated with molecular mediators of pre-B cell negative selection. Our data is the foundation for the interactive online Functional Omics Resource of ALL (FORALL) with navigable proteomics, transcriptomics, and drug sensitivity profiles at https://proteomics.se/forall .

Proteogenomics refines the molecular classification of chronic lymphocytic leukemia.

Cancer heterogeneity at the proteome level may explain differences in therapy response and prognosis beyond the currently established genomic and transcriptomic-based diagnostics. The relevance of proteomics for disease classifications remains to be established in clinically heterogeneous cancer entities such as chronic lymphocytic leukemia (CLL). Here, we characterize the proteome and transcriptome alongside genetic and ex-vivo drug response profiling in a clinically annotated CLL discovery cohort (n = 68). Unsupervised clustering of the proteome data reveals six subgroups. Five of these proteomic groups are associated with genetic features, while one group is only detectable at the proteome level. This new group is characterized by accelerated disease progression, high spliceosomal protein abundances associated with aberrant splicing, and low B cell receptor signaling protein abundances (ASB-CLL). Classifiers developed to identify ASB-CLL based on its characteristic proteome or splicing signature in two independent cohorts (n = 165, n = 169) confirm that ASB-CLL comprises about 20% of CLL patients. The inferior overall survival in ASB-CLL is also independent of both TP53- and IGHV mutation status. Our multi-omics analysis refines the classification of CLL and highlights the potential of proteomics to improve cancer patient stratification beyond genetic and transcriptomic profiling.

Optimization of production of the anti-keratin 8 single-chain Fv TS1-218 in Pichia pastoris using design of experiments.

BACKGROUND: Optimization of conditions during recombinant protein production for improved yield is a major goal for protein scientists. Typically this is achieved by changing single crucial factor settings one at a time while other factors are kept fixed through trial-and-error experimentation. This approach may introduce larger bias and fail to identify interactions between the factors resulting in failure of finding the true optimal conditions. RESULTS: In this study we have utilized design of experiments in order to identify optimal culture conditions with the aim to improve the final yield of the anti-keratin 8 scFv TS1-218, during expression in P. pastoris in shake flasks. The effect of: pH, temperature and methanol concentration on the yield of TS1-218 using buffered minimal medium was investigated and a predictive model established. The results demonstrated that higher starting pH and lower temperatures during induction significantly increased the yield of TS1-218. Furthermore, the result demonstrated increased biomass accumulation and cell viability at lower temperatures which suggested that the higher yield of TS1-218 could be attributed to lower protease activity in the culture medium. The optimal conditions (pH 7.1, temperature 11°C and methanol concentration 1.2%) suggested by the predictive model yielded 21.4 mg TS1-218 which is a 21-fold improvement compared to the yield prior to optimization. CONCLUSION: The results demonstrated that design of experiments can be utilized for a rapid optimization of initial culture conditions and that P. pastoris is highly capable of producing and secreting functional single-chain antibody fragments at temperatures as low as 11°C.

Inhibition of the ubiquitin-proteasome system by an NQO1-activatable compound.

Malignant cells display an increased sensitivity towards drugs that reduce the function of the ubiquitin-proteasome system (UPS), which is the primary proteolytic system for destruction of aberrant proteins. Here, we report on the discovery of the bioactivatable compound CBK77, which causes an irreversible collapse of the UPS, accompanied by a general accumulation of ubiquitylated proteins and caspase-dependent cell death. CBK77 caused accumulation of ubiquitin-dependent, but not ubiquitin-independent, reporter substrates of the UPS, suggesting a selective effect on ubiquitin-dependent proteolysis. In a genome-wide CRISPR interference screen, we identified the redox enzyme NAD(P)H:quinone oxidoreductase 1 (NQO1) as a critical mediator of CBK77 activity, and further demonstrated its role as the compound bioactivator. Through affinity-based proteomics, we found that CBK77 covalently interacts with ubiquitin. In vitro experiments showed that CBK77-treated ubiquitin conjugates were less susceptible to disassembly by deubiquitylating enzymes. In vivo efficacy of CBK77 was validated by reduced growth of NQO1-proficient human adenocarcinoma cells in nude mice treated with CBK77. This first-in-class NQO1-activatable UPS inhibitor suggests that it may be possible to exploit the intracellular environment in malignant cells for leveraging the impact of compounds that impair the UPS.

Novel Broad-Spectrum Antiviral Inhibitors Targeting Host Factors Essential for Replication of Pathogenic RNA Viruses.

Recent RNA virus outbreaks such as Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and Ebola virus (EBOV) have caused worldwide health emergencies highlighting the urgent need for new antiviral strategies. Targeting host cell pathways supporting viral replication is an attractive approach for development of antiviral compounds, especially with new, unexplored viruses where knowledge of virus biology is limited. Here, we present a strategy to identify host-targeted small molecule inhibitors using an image-based phenotypic antiviral screening assay followed by extensive target identification efforts revealing altered cellular pathways upon antiviral compound treatment. The newly discovered antiviral compounds showed broad-range antiviral activity against pathogenic RNA viruses such as SARS-CoV-2, EBOV and Crimean-Congo hemorrhagic fever virus (CCHFV). Target identification of the antiviral compounds by thermal protein profiling revealed major effects on proteostasis pathways and disturbance in interactions between cellular HSP70 complex and viral proteins, illustrating the supportive role of HSP70 on many RNA viruses across virus families. Collectively, this strategy identifies new small molecule inhibitors with broad antiviral activity against pathogenic RNA viruses, but also uncovers novel virus biology urgently needed for design of new antiviral therapies.

Proteogenomics and Hi-C reveal transcriptional dysregulation in high hyperdiploid childhood acute lymphoblastic leukemia.

Hyperdiploidy, i.e. gain of whole chromosomes, is one of the most common genetic features of childhood acute lymphoblastic leukemia (ALL), but its pathogenetic impact is poorly understood. Here, we report a proteogenomic analysis on matched datasets from genomic profiling, RNA-sequencing, and mass spectrometry-based analysis of >8,000 genes and proteins as well as Hi-C of primary patient samples from hyperdiploid and ETV6/RUNX1-positive pediatric ALL. We show that CTCF and cohesin, which are master regulators of chromatin architecture, display low expression in hyperdiploid ALL. In line with this, a general genome-wide dysregulation of gene expression in relation to topologically associating domain (TAD) borders were seen in the hyperdiploid group. Furthermore, Hi-C of a limited number of hyperdiploid childhood ALL cases revealed that 2/4 cases displayed a clear loss of TAD boundary strength and 3/4 showed reduced insulation at TAD borders, with putative leukemogenic effects.

Functional mapping and single chain construction of the anti-cytokeratin 8 monoclonal antibody TS1.

The anti-cytokeratin (CK) 8 monoclonal antibody (mab) TS1 has been shown to efficiently bind to CK8 expressed in carcinomas in vivo. The anti-idiotypic antibody of TS1, alphaTS1, can be used to regulate the tumor:non-tumor ratio of TS1 by clearing non-tumor binding TS1 from the circulation. If the interaction of TS1 to CK8 and alphaTS1 is fully understood, mutations can be used to improve the tumor:non-tumor ratio. A scFv was made of the mab TS1 and residues earlier identified by Erlandsson et al. as important for the interaction with both its antigen CK8 and its anti-idiotype alphaTS1, were mutated to alanine or amides and expressed in E. coli. The effects of the mutations were studied by ELISA and residues important for the interactions to both CK8 and alphaTS1 were identified as mainly tyrosines, charged residues, a serine and a tryptophan. Altogether, nine amino acid residues in TS1 were found to be important in the interaction to alphaTS1 and six residues for the interaction to CK8. Important residues, clustered together in the modelled protein, were identified as residues from CDR 3 of the heavy chain and the unexpected participation of a residue in CDR 2 of the light chain. Some of the important residues are likely to be hotspots. Hotspots constitute a few residues in an interaction that contribute most to the binding, energetically. Amino acid residues in hotspots often cluster together in the center of the interaction interface, but can also be spread out to the periphery. The hotspots are often surrounded by hydrophobic patches, which are seen in the modelled TS1 protein used in this study. Amino acid residues that increased the affinity when mutated were also identified for both interactions. These residues are likely to be located outside the interacting interface. It can from this study be concluded that it is wise to precede the mutational procedure with experiments that can give guidelines for the selection of which amino acid residues to mutate. If the guidelines from the chemical modifications from Erlandsson et al. not had been used, this study would have left some residues unmutated and thereby missed important information.

CETSA screening identifies known and novel thymidylate synthase inhibitors and slow intracellular activation of 5-fluorouracil.

Target engagement is a critical factor for therapeutic efficacy. Assessment of compound binding to native target proteins in live cells is therefore highly desirable in all stages of drug discovery. We report here the first compound library screen based on biophysical measurements of intracellular target binding, exemplified by human thymidylate synthase (TS). The screen selected accurately for all the tested known drugs acting on TS. We also identified TS inhibitors with novel chemistry and marketed drugs that were not previously known to target TS, including the DNA methyltransferase inhibitor decitabine. By following the cellular uptake and enzymatic conversion of known drugs we correlated the appearance of active metabolites over time with intracellular target engagement. These data distinguished a much slower activation of 5-fluorouracil when compared with nucleoside-based drugs. The approach establishes efficient means to associate drug uptake and activation with target binding during drug discovery.

The cellular thermal shift assay for evaluating drug target interactions in cells.

Thermal shift assays are used to study thermal stabilization of proteins upon ligand binding. Such assays have been used extensively on purified proteins in the drug discovery industry and in academia to detect interactions. Recently, we published a proof-of-principle study describing the implementation of thermal shift assays in a cellular format, which we call the cellular thermal shift assay (CETSA). The method allows studies of target engagement of drug candidates in a cellular context, herein exemplified with experimental data on the human kinases p38α and ERK1/2. The assay involves treatment of cells with a compound of interest, heating to denature and precipitate proteins, cell lysis, and the separation of cell debris and aggregates from the soluble protein fraction. Whereas unbound proteins denature and precipitate at elevated temperatures, ligand-bound proteins remain in solution. We describe two procedures for detecting the stabilized protein in the soluble fraction of the samples. One approach involves sample workup and detection using quantitative western blotting, whereas the second is performed directly in solution and relies on the induced proximity of two target-directed antibodies upon binding to soluble protein. The latter protocol has been optimized to allow an increased throughput, as potential applications require large numbers of samples. Both approaches can be completed in a day.

Tracking cancer drugs in living cells by thermal profiling of the proteome.

The thermal stability of proteins can be used to assess ligand binding in living cells. We have generalized this concept by determining the thermal profiles of more than 7000 proteins in human cells by means of mass spectrometry. Monitoring the effects of small-molecule ligands on the profiles delineated more than 50 targets for the kinase inhibitor staurosporine. We identified the heme biosynthesis enzyme ferrochelatase as a target of kinase inhibitors and suggest that its inhibition causes the phototoxicity observed with vemurafenib and alectinib. Thermal shifts were also observed for downstream effectors of drug treatment. In live cells, dasatinib induced shifts in BCR-ABL pathway proteins, including CRK/CRKL. Thermal proteome profiling provides an unbiased measure of drug-target engagement and facilitates identification of markers for drug efficacy and toxicity.

Monitoring drug target engagement in cells and tissues using the cellular thermal shift assay.

The efficacy of therapeutics is dependent on a drug binding to its cognate target. Optimization of target engagement by drugs in cells is often challenging, because drug binding cannot be monitored inside cells. We have developed a method for evaluating drug binding to target proteins in cells and tissue samples. This cellular thermal shift assay (CETSA) is based on the biophysical principle of ligand-induced thermal stabilization of target proteins. Using this assay, we validated drug binding for a set of important clinical targets and monitored processes of drug transport and activation, off-target effects and drug resistance in cancer cell lines, as well as drug distribution in tissues. CETSA is likely to become a valuable tool for the validation and optimization of drug target engagement.

Construction of divalent anti-keratin 8 single-chain antibodies (sc(Fv)(2)), expression in Pichia pastoris and their reactivity with multicellular tumor spheroids.

Single-chain variable fragments (scFvs) are small monovalent recombinant antibody fragments that retain the specificity of their parent immunoglobulins. ScFvs are excellent building blocks for new and improved immunodiagnostic and therapeutic proteins. However, the monovalency and the rapid renal elimination of scFvs result in poor tumor accumulation and retention. Engineering divalent antibody fragments is an excellent way to address these shortcomings. In this study, covalent divalent single-chain variable fragments (sc(Fv)(2)s), were constructed from the monovalent anti-keratin 8 scFvs, TS1-218 and its mutant, HE1-Q. The scFvs and sc(Fv)(2)s were expressed in the methylotrophic yeast Pichia pastoris, utilizing the alpha-factor secretion signal (α-factor) for extracellular secretion. The immunoreactivity and specificity of the antibody fragments were analyzed with enzyme-linked immunosorbent assay (ELISA) and the uptake and retention of the (125)I labeled antibody fragments were evaluated using HeLa HEp-2 multicellular tumor spheroids (MCTSs). Analysis of the antibody fragments demonstrated that parts of the α-factor remained at the N-terminal of the antibody fragments. Despite incomplete processing of the α-factor, the antibody fragments were functional where the sc(Fv)(2)s gave a three-fold stronger signal in ELISA compared to their scFv counterparts and the mutant antibodies demonstrated a stronger signal than their initial wild types. In addition, the sc(Fv)(2)s DiTS1-218 and DiHE1-Q displayed an approximately two-fold higher uptake and were retained to a larger extent in the MCTS, demonstrating a 3.9 and 9.4-fold increase in half-life respectively compared to their corresponding scFvs. In conclusion, expression in P. pastoris improved the yield 20-fold and facilitated the purification of the antibody fragments. Furthermore, the sc(Fv)(2)s presented a higher functional affinity to K 8 both in ELISA and MCTS compared to the scFvs with DiHE1-Q being the best candidate for further studies.