A proteome-wide atlas of drug mechanism of action

Polypharmacology prediction: the long road toward comprehensively anticipating small-molecule selectivity to de-risk drug discovery

INTRODUCTION

Small molecules often bind to multiple targets, a behavior termed polypharmacology. Anticipating polypharmacology is essential for drug discovery since unknown off-targets can modulate safety and efficacy - profoundly affecting drug discovery success. Unfortunately, experimental methods to assess selectivity present significant limitations and drugs still fail in the clinic due to unanticipated off-targets. Computational methods are a cost-effective, complementary approach to predict polypharmacology.

AREAS COVERED

This review aims to provide a comprehensive overview of the state of polypharmacology prediction and discuss its strengths and limitations, covering both classical cheminformatics methods and bioinformatic approaches. The authors review available data sources, paying close attention to their different coverage. The authors then discuss major algorithms grouped by the types of data that they exploit using selected examples.

EXPERT OPINION

Polypharmacology prediction has made impressive progress over the last decades and contributed to identify many off-targets. However, data incompleteness currently limits most approaches to comprehensively predict selectivity. Moreover, our limited agreement on model assessment challenges the identification of the best algorithms - which at present show modest performance in prospective real-world applications. Despite these limitations, the exponential increase of multidisciplinary Big Data and AI hold much potential to better polypharmacology prediction and de-risk drug discovery.

HDAC inhibitors activate lipid peroxidation and ferroptosis in gastric cancer

Gastric cancer remains among the deadliest neoplasms worldwide, with limited therapeutic options. Since efficacies of targeted therapies are unsatisfactory, drugs with broader mechanisms of action rather than a single oncogene inhibition are needed. Preclinical studies have identified histone deacetylases (HDAC) as potential therapeutic targets in gastric cancer. However, the mechanism(s) of action of HDAC inhibitors (HDACi) are only partially understood. This is particularly true with regard to ferroptosis as an emerging concept of cell death. In a panel of gastric cancer cell lines with different molecular characteristics, tumor cell inhibitory effects of different HDACi were studied. Lipid peroxidation levels were measured and proteome analysis was performed for the in-depth characterization of molecular alterations upon HDAC inhibition. HDACi effects on important ferroptosis genes were validated on the mRNA and protein level. Upon HDACi treatment, lipid peroxidation was found increased in all cell lines. Class I HDACi (VK1, entinostat) showed the same toxicity profile as the pan-HDACi vorinostat. Proteome analysis revealed significant and concordant alterations in the expression of proteins related to ferroptosis induction. Key enzymes like ACSL4, POR or SLC7A11 showed distinct alterations in their expression patterns, providing an explanation for the increased lipid peroxidation. Results were also confirmed in primary human gastric cancer tissue cultures as a relevant ex vivo model. We identify the induction of ferroptosis as new mechanism of action of class I HDACi in gastric cancer. Notably, these findings were independent of the genetic background of the cell lines, thus introducing HDAC inhibition as a more general therapeutic principle.

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.

Molecular glues and induced proximity: An evolution of tools and discovery

Small molecule molecular glues can nucleate protein complexes and rewire interactomes. Molecular glues are widely used as probes for understanding functional proximity at a systems level, and the potential to instigate event-driven pharmacology has motivated their application as therapeutics. Despite advantages such as cell permeability and the potential for low off-target activity, glues are still rare when compared to canonical inhibitors in therapeutic development. Their often simple structure and specific ability to reshape protein-protein interactions pose several challenges for widespread, designer applications. Molecular glue discovery and design campaigns can find inspiration from the fields of synthetic biology and biophysics to mine chemical libraries for glue-like molecules.

Unraveling the role of plasma proteins in dementia: insights from two cohort studies in the UK, with causal evidence from Mendelian randomization

Population-based proteomics offer a groundbreaking avenue to predict dementia onset. This study employed a proteome-wide, data-driven approach to investigate protein-dementia associations in 229 incident all-cause dementia (ACD) among 3,249 participants from the English Longitudinal Study of Ageing (ELSA) over a median 9.8-year follow-up, then validated in 1,506 incident ACD among 52,745 individuals from the UK Biobank (UKB) over median 13.7 years. NEFL and RPS6KB1 were robustly associated with incident ACD; MMP12 was associated with vascular dementia in ELSA. Additional markers EDA2R and KIM1 (HAVCR1) were identified from sensitivity analyses. Combining NEFL and RPS6KB1 with other factors yielded high predictive accuracy (area under the curve (AUC)=0.871) for incident ACD. Replication in the UKB confirmed associations between identified proteins with various dementia subtypes. Results from reverse Mendelian Randomization also supported the role of several proteins as early dementia biomarkers. These findings underscore proteomics' potential in identifying novel risk screening targets for dementia.

Proteomics analysis of histone deacetylase inhibitor-resistant solid tumors reveals resistant signatures and potential drug combinations

Histone deacetylase inhibitors (HDACis) are important drugs for cancer therapy, but the indistinct resistant mechanisms of solid tumor therapy greatly limit their clinical application. In this study we conducted HDACi-perturbated proteomics and phosphoproteomics analyses in HDACi-sensitive and -resistant cell lines using a tandem mass tag (TMT)-based quantitative proteomic strategy. We found that the ribosome biogenesis proteins MRTO4, PES1, WDR74 and NOP16 vital to tumorigenesis might regulate the tumor sensitivity to HDACi. By integrating HDACi-perturbated protein signature with previously reported proteomics and drug sensitivity data, we predicted and validated a series of drug combination pairs potentially to enhance the sensitivity of HDACi in diverse solid tumor. Functional phosphoproteomic analysis further identified the kinase PDK1 and ROCK as potential HDACi-resistant signatures. Overall, this study reveals the potential HDACi-resistant signatures and may provide promising drug combination strategies to attenuate the resistance of solid tumor to HDACi.

Assessing Target Specificity of the Small Molecule Inhibitor MARIMASTAT to Snake Venom Toxins: A Novel Application of Thermal Proteome Profiling

New treatments that circumvent the pitfalls of traditional antivenom therapies are critical to address the problem of snakebite globally. Numerous snake venom toxin inhibitors have shown promising cross-species neutralization of medically significant venom toxins in vivo and in vitro. The development of high-throughput approaches for the screening of such inhibitors could accelerate their identification, testing, and implementation and thus holds exciting potential for improving the treatments and outcomes of snakebite envenomation worldwide. Energetics-based proteomic approaches, including thermal proteome profiling and proteome integral solubility alteration (PISA) assays, represent "deep proteomics" methods for high throughput, proteome-wide identification of drug targets and ligands. In the following study, we apply thermal proteome profiling and PISA methods to characterize the interactions between venom toxin proteoforms in Crotalus atrox (Western Diamondback Rattlesnake) and the snake venom metalloprotease (SVMP) inhibitor marimastat. We investigate its venom proteome-wide effects and characterize its interactions with specific SVMP proteoforms, as well as its potential targeting of non-SVMP venom toxin families. We also compare the performance of PISA thermal window and soluble supernatant with insoluble precipitate using two inhibitor concentrations, providing the first demonstration of the utility of a sensitive high-throughput PISA-based approach to assess the direct targets of small molecule inhibitors for snake venom.

Lineage-specific proteome remodeling of diverse lung cancer cells by targeted epigenetic inhibitors

Epigenetic inhibitors exhibit powerful antiproliferative and anticancer activities. However, cellular responses to small-molecule epigenetic inhibition are heterogenous and dependent on factors such as the genetic background, metabolic state, and on-/off-target engagement of individual small-molecule drugs. To determine the mechanisms that drive these heterogeneous cellular responses, we quantified chromatin, proteome, and transcriptome remodeling due to histone deacetylase inhibitor (HDACi) -treated cells derived from diverse genetic backgrounds. We utilized high-throughput sample multiplexed proteomics and integrated intelligent data acquisition methods to map proteomes of cancer cell lines in response to HDACi. We determined cell type-specific and ubiquitous cellular responses based on the quantification of 10,621 total proteins. We then established how coordinated remodeling of the proteome, transcriptome and chromatin state of HDACi treated cancer cells revealed convergent (JUN, MAP2K3, CDKN1A) and divergent (CCND3, ASF1B, BRD7) molecular phenotypes. HDACi-regulated proteins differ greatly across cell lines owing to heterogeneous molecular states of these cell lines. Finally, we demonstrated that HDACi treatment drove a highly cell-type specific response that may in part be explained by cell line-specific off-target drug engagement.

Mass Spectrometry-Based Proteomics Analysis Unveils PTPRS Inhibits Proliferation and Inflammatory Response of Keratinocytes in Psoriasis

In this study, we used data-independent acquisition-mass spectrometry (DIA-MS) to analyze the serum proteome in psoriasis vulgaris (PsO). The serum proteomes of seven healthy controls and eight patients with PsO were analyzed using DIA-MS. Weighted gene co-expression network analysis was used to identify differentially expressed proteins (DEPs) that were closely related to PsO. Hub proteins of PsO were also identified. The Proteomics Drug Atlas 2023 was used to predict candidate hub protein drugs. To confirm the expression of the candidate factor, protein tyrosine phosphatase receptor S (PTPRS), in psoriatic lesions and the psoriatic keratinocyte model, immunohistochemical staining, quantitative real-time polymerase chain reaction, and western blotting were performed. A total of 129 DEPs were found to be closely related to PsO. The hub proteins for PsO were PVRL1, FGFR1, PTPRS, CDH2, CDH1, MCAM, and THY1. Five candidate hub protein drugs were identified: encorafenib, leupeptin, fedratinib, UNC 0631, and SCH 530348. PTPRS was identified as a common pharmacological target for these five drugs. PTPRS knockdown in keratinocytes promoted the proliferation and expression of IL-1α, IL-1β, IL-23A, TNF-α, MMP9, CXCL8, and S100A9. PTPRS expression was decreased in PsO, and PTPRS negatively regulated PsO. PTPRS may be involved in PsO pathogenesis through the inhibition of keratinocyte proliferation and inflammatory responses and is a potential treatment target for PsO.

Chromatin damage generated by DNA intercalators leads to degradation of RNA Polymerase II

In cancer therapy, DNA intercalators are mainly known for their capacity to kill cells by inducing DNA damage. Recently, several DNA intercalators have attracted much interest given their ability to inhibit RNA Polymerase I transcription (BMH-21), evict histones (Aclarubicin) or induce chromatin trapping of FACT (Curaxin CBL0137). Interestingly, these DNA intercalators lack the capacity to induce DNA damage while still retaining cytotoxic effects and stabilize p53. Herein, we report that these DNA intercalators impact chromatin biology by interfering with the chromatin stability of RNA polymerases I, II and III. These three compounds have the capacity to induce degradation of RNA polymerase II and they simultaneously enable the trapping of Topoisomerases TOP2A and TOP2B on the chromatin. In addition, BMH-21 also acts as a catalytic inhibitor of Topoisomerase II, resembling Aclarubicin. Moreover, BMH-21 induces chromatin trapping of the histone chaperone FACT and propels accumulation of Z-DNA and histone eviction, similarly to Aclarubicin and CBL0137. These DNA intercalators have a cumulative impact on general transcription machinery by inducing accumulation of topological defects and impacting nuclear chromatin. Therefore, their cytotoxic capabilities may be the result of compounding deleterious effects on chromatin homeostasis.

Multikinase inhibitors modulate non-constitutive proteasome expression in colorectal cancer cells

Introduction: Proteasomes are multi-subunit protein complexes responsible for protein degradation in cells. Immunoproteasomes and intermediate proteasomes (together non-constitutive proteasomes) are specific forms of proteasomes frequently associated with immune response, antigen presentation, inflammation and stress. Expression of non-constitutive proteasome subunits has a prognostic value in several types of cancer. Thus, factors that modulate non-constitutive proteasome expression in tumors are of particular interest. Multikinase inhibitors (MKIs) demonstrate promising results in treatment of cancer. At the same time, their immunomodulatory properties and effects on non-constitutive proteasome expression in colorectal cancer cells are poorly investigated. Methods: Proteasome subunit expression in colorectal cancer was evaluated by bioinformatic analysis of available datasets. Two colorectal cancer cell lines, expressing fluorescent non-constitutive proteasomes were treated with multikinase inhibitors: regorafenib and sorafenib. The proteasome subunit expression was assessed by real-time PCR, Western blotting and flow cytometry. The proteasome activity was studied using proteasome activity-based probe and fluorescent substrates. Intracellular proteasome localization was revealed by confocal microscopy. Reactive oxygen species levels following treatment were determined in cells. Combined effect of proteasome inhibition and treatment with MKIs on viability of cells was estimated. Results: Expression of non-constitutive proteasomes is increased in BRAF-mutant colorectal tumors. Regorafenib and sorafenib stimulated the activity and synthesis of non-constitutive proteasomes in examined cell lines. MKIs induced oxidative stress and redistribution of proteasomes within cells. Sorafenib stimulated formation of cytoplasmic aggregates, containing proteolyticaly active non-constitutive proteasomes, while regorafenib had no such effect. MKIs caused no synergistic action when were combined with the proteasome inhibitor. Discussion: Obtained results indicate that MKIs might affect the crosstalk between cancer cells and immune cells via modulation of intracellular proteasome pool. Observed phenomenon should be considered when MKI-based therapy is applied.

Pooled multicolour tagging for visualizing subcellular protein dynamics

Imaging-based methods are widely used for studying the subcellular localization of proteins in living cells. While routine for individual proteins, global monitoring of protein dynamics following perturbation typically relies on arrayed panels of fluorescently tagged cell lines, limiting throughput and scalability. Here, we describe a strategy that combines high-throughput microscopy, computer vision and machine learning to detect perturbation-induced changes in multicolour tagged visual proteomics cell (vpCell) pools. We use genome-wide and cancer-focused intron-targeting sgRNA libraries to generate vpCell pools and a large, arrayed collection of clones each expressing two different endogenously tagged fluorescent proteins. Individual clones can be identified in vpCell pools by image analysis using the localization patterns and expression level of the tagged proteins as visual barcodes, enabling simultaneous live-cell monitoring of large sets of proteins. To demonstrate broad applicability and scale, we test the effects of antiproliferative compounds on a pool with cancer-related proteins, on which we identify widespread protein localization changes and new inhibitors of the nuclear import/export machinery. The time-resolved characterization of changes in subcellular localization and abundance of proteins upon perturbation in a pooled format highlights the power of the vpCell approach for drug discovery and mechanism-of-action studies.

A heterogeneous pharmaco-transcriptomic landscape induced by targeting a single oncogenic kinase

Over-activation of the epidermal growth factor receptor (EGFR) is a hallmark of glioblastoma. However, EGFR-targeted therapies have led to minimal clinical response. While delivery of EGFR inhibitors (EGFRis) to the brain constitutes a major challenge, how additional drug-specific features alter efficacy remains poorly understood. We apply highly multiplex single-cell chemical genomics to define the molecular response of glioblastoma to EGFRis. Using a deep generative framework, we identify shared and drug-specific transcriptional programs that group EGFRis into distinct molecular classes. We identify programs that differ by the chemical properties of EGFRis, including induction of adaptive transcription and modulation of immunogenic gene expression. Finally, we demonstrate that pro-immunogenic expression changes associated with a subset of tyrphostin family EGFRis increase the ability of T-cells to target glioblastoma cells.

[Deciphering cellular processes responding to lethality of 17β-estradiol by quantitative phosphoproteomics]

17β-Estradiol (E2), an important endocrine hormone in the mammalian body, participates in the regulation of the physiological functions of the reproductive system, mammary glands, bone, and cardiovascular system, among others. Paradoxically, despite the physiological actions of endogenous E2 (0.2-1.0 nmol/L), numerous clinical and experimental studies have demonstrated that high-dose E2 treatment can cause tumor regression and exert pro-apoptotic actions in multiple cell types; however, the underlying mechanism remains undescribed. In particular, little information of the cellular processes responding to the lethality of E2 is available. In the present study, we attempted to characterize the cellular processes responding to high-dose (μmol/L) E2 treatment using quantitative phosphoproteomics to obtain a better understanding of the regulatory mechanism of E2-induced cell death. First, the cell phenotype induced by high-dose E2 was determined by performing Cell Counting Kit-8 assay (CCK8), cell cytotoxicity analysis by trypan blue staining, and microscopic imaging on HeLa cells treated with 1-10 μmol/L E2 or dimethyl sulfoxide (DMSO) for 1-3 d. E2 inhibited cell proliferation and induced cell death in a dose- and time-dependent manner. Compared with the DMSO-treated HeLa cells, the cells treated with 5 μmol/L E2 for 2 d demonstrated >74% growth inhibition and approximately 50% cell death. Thus, these cells were used for quantitative phosphoproteomic analysis. Next, a solid-phase extraction (SPE)-based immobilized titanium ion affinity chromatography (Ti4+-IMAC) phosphopeptide-enrichment method coupled with data-independent acquisition (DIA)-based quantitative proteomics was employed for the in-depth screening of high-dose E2-regulated phosphorylation sites to investigate the intracellular processes responding to high-dose E2 treatment. Liquid chromatography-tandem mass spectrometry (LC-MS/MS) identified over 10000 phosphorylation sites regulated by E2 and DMSO in HeLa cells. In comparison with the DMSO-treated cells, the cells treated with 5 μmol/L E2 showed 537 upregulated phosphorylation sites and 387 downregulated phosphorylation sites, with a threshold of p<0.01 and |log2(fold change)|≥1. A total of 924 phosphorylation sites on 599 proteins were significantly regulated by high-dose E2, and these sites were subjected to enrichment analysis. In addition, 453 differently regulated phosphorylation sites on 325 proteins were identified only in the E2- or DMSO-treated cell samples. These phosphorylation sites may be phosphorylated or dephosphorylated in response to high-dose E2 stimulation and were subjected to parallel enrichment analyses. Taken together, 1218 phosphorylation sites on 741 proteins were significantly regulated by high-dose E2 treatment. The functional phosphoproteins in these two groups were then analyzed using Gene Ontology (GO) and Gene Set Enrichment Analysis (GSEA) to determine the biological processes in which they participate and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway database. Consistent with the cell-phenotype data, cell cycle-related proteins were highly enriched in the two groups of E2-regulated phosphoproteins (p<0.05), indicating that high-dose E2 treatment can regulate cell proliferation. In addition, E2-regulated phosphoproteins were highly enriched in the cellular processes of ribosome biogenesis, nucleocytoplasmic transport, and messenger ribonucleic acid (mRNA) processing/splicing (p<0.05), indicating that the activation of these processes may contribute to high-dose E2-induced cell death. These results further confirm that high-dose E2 treatment inhibits protein translation and induces cell death. Furthermore, the significant upregulation of multiple phosphorylation sites associated with epidermal growth factor receptor (EGFR) and mitogen-activated protein kinases (MAPKs) MAPK1, MAPK4, and MAPK14 by high-dose E2 indicates that the EGFR and MAPK signaling pathways are likely involved in the regulation of E2-induced cell death. These phosphorylation sites likely play vital roles in E2-induced cell death in HeLa cells. Overall, our phosphoproteomic data could be a valuable resource for uncovering the regulatory mechanisms of E2 in the micromolar range.

Accelerating multiplexed profiling of protein-ligand interactions: High-throughput plate-based reactive cysteine profiling with minimal input

Chemoproteomics has made significant progress in investigating small-molecule-protein interactions. However, the proteome-wide profiling of cysteine ligandability remains challenging to adapt for high-throughput applications, primarily due to a lack of platforms capable of achieving the desired depth using low input in 96- or 384-well plates. Here, we introduce a revamped, plate-based platform which enables routine interrogation of either ∼18,000 or ∼24,000 reactive cysteines based on starting amounts of 10 or 20 μg, respectively. This represents a 5-10X reduction in input and 2-3X improved coverage. We applied the platform to screen 192 electrophiles in the native HEK293T proteome, mapping the ligandability of 38,450 reactive cysteines from 8,274 human proteins. We further applied the platform to characterize new cellular targets of established drugs, uncovering that ARS-1620, a KRASG12C inhibitor, binds to and inhibits an off-target adenosine kinase ADK. The platform represents a major step forward to high-throughput proteome-wide evaluation of reactive cysteines.

Mass spectrometry-based proteomic analysis to characterize cisplatin induced early signaling events in head and neck squamous cell carcinoma

Cisplatin is the commonly used chemotherapeutic drug in treatment of various cancers. However, development of resistance towards cisplatin results in tumor recurrence. Here, we aim to understand the mechanisms of action of cisplatin and emergence of resistance to cisplatin using mass spectrometry-based proteomic approach. A panel of head and neck squamous cell carcinoma (HNSCC) cell lines were treated with cisplatin at respective IC50 for 24 h and label-free mass spectrometry analysis was carried out. Proteomic analysis of A253, FaDu, Det562 and CAL27 cell lines upon cisplatin treatment resulted in the identification of 5,060, 4,816, 4,537 and 4,142 proteins, respectively. Bioinformatics analysis of differentially regulated proteins revealed proteins implicated in DNA damage bypass pathway, translation and mRNA splicing to be enriched. Further, proteins associated with cisplatin resistance exhibited alterations following short-term cisplatin exposure. Among these, class III tubullin protein (TUBB3) was found to be upregulated in cisplatin-treated cells compared to untreated cells. Western blot analysis confirmed the elevated expression of TUBB3 in cells treated with cisplatin for 24 h, and also in cisplatin resistant HNSCC cell lines. This study delineates the early signaling events that enable HNSCC cells to counteract the cytotoxic effects of cisplatin and facilitate the development of resistance.

Uncovering PROTAC Sensitivity and Efficacy by Multidimensional Proteome Profiling: A Case for STAT3

Proteolysis-targeting chimera (PROTAC) is a powerful technology that can effectively trigger the degradation of target proteins. The intricate interplay among various factors leads to a heterogeneous drug response, bringing about significant challenges in comprehending drug mechanisms. Our study applied data-independent acquisition-based mass spectrometry to multidimensional proteome profiling of PROTAC (DIA-MPP) to uncover the efficacy and sensitivity of the PROTAC compound. We profiled the signal transducer and activator of transcription 3 (STAT3) PROTAC degrader in six leukemia and lymphoma cell lines under multiple conditions, demonstrating the pharmacodynamic properties and downstream biological responses. Through comparison between sensitive and insensitive cell lines, we revealed that STAT1 can be regarded as a biomarker for STAT3 PROTAC degrader, which was validated in cells, patient-derived organoids, and mouse models. These results set an example for a comprehensive description of the multidimensional PROTAC pharmacodynamic response and PROTAC drug sensitivity biomarker exploration.

Intracellular proteomics and extracellular vesiculomics as a metric of disease recapitulation in 3D-bioprinted aortic valve arrays

In calcific aortic valve disease (CAVD), mechanosensitive valvular cells respond to fibrosis- and calcification-induced tissue stiffening, further driving pathophysiology. No pharmacotherapeutics are available to treat CAVD because of the paucity of (i) appropriate experimental models that recapitulate this complex environment and (ii) benchmarking novel engineered aortic valve (AV)-model performance. We established a biomaterial-based CAVD model mimicking the biomechanics of the human AV disease-prone fibrosa layer, three-dimensional (3D)-bioprinted into 96-well arrays. Liquid chromatography-tandem mass spectrometry analyses probed the cellular proteome and vesiculome to compare the 3D-bioprinted model versus traditional 2D monoculture, against human CAVD tissue. The 3D-bioprinted model highly recapitulated the CAVD cellular proteome (94% versus 70% of 2D proteins). Integration of cellular and vesicular datasets identified known and unknown proteins ubiquitous to AV calcification. This study explores how 2D versus 3D-bioengineered systems recapitulate unique aspects of human disease, positions multiomics as a technique for the evaluation of high throughput-based bioengineered model systems, and potentiates future drug discovery.

High-throughput drug target discovery using a fully automated proteomics sample preparation platform

Drug development is plagued by inefficiency and high costs due to issues such as inadequate drug efficacy and unexpected toxicity. Mass spectrometry (MS)-based proteomics, particularly isobaric quantitative proteomics, offers a solution to unveil resistance mechanisms and unforeseen side effects related to off-targeting pathways. Thermal proteome profiling (TPP) has gained popularity for drug target identification at the proteome scale. However, it involves experiments with multiple temperature points, resulting in numerous samples and considerable variability in large-scale TPP analysis. We propose a high-throughput drug target discovery workflow that integrates single-temperature TPP, a fully automated proteomics sample preparation platform (autoSISPROT), and data independent acquisition (DIA) quantification. The autoSISPROT platform enables the simultaneous processing of 96 samples in less than 2.5 hours, achieving protein digestion, desalting, and optional TMT labeling (requires an additional 1 hour) with 96-channel all-in-tip operations. The results demonstrated excellent sample preparation performance with >94% digestion efficiency, >98% TMT labeling efficiency, and >0.9 intra- and inter-batch Pearson correlation coefficients. By automatically processing 87 samples, we identified both known targets and potential off-targets of 20 kinase inhibitors, affording over a 10-fold improvement in throughput compared to classical TPP. This fully automated workflow offers a high-throughput solution for proteomics sample preparation and drug target/off-target identification.

Special Issue "Deployment of Proteomics Approaches in Biomedical Research"

Many angles of personalized medicine, such as diagnostic improvements, systems biology [...].

Highly specific intracellular ubiquitination of a small molecule

Ubiquitin is a small, highly conserved protein that acts as a posttranslational modification in eukaryotes. Ubiquitination of proteins frequently serves as a degradation signal, marking them for disposal by the proteasome. Here, we report a novel small molecule from a diversity-oriented synthesis library, BRD1732, that is directly ubiquitinated in cells, resulting in dramatic accumulation of inactive ubiquitin monomers and polyubiquitin chains causing broad inhibition of the ubiquitin-proteasome system. Ubiquitination of BRD1732 and its associated cytotoxicity are stereospecific and dependent upon two homologous E3 ubiquitin ligases, RNF19A and RNF19B. Our finding opens the possibility for indirect ubiquitination of a target through a ubiquitinated bifunctional small molecule, and more broadly raises the potential for posttranslational modification in trans .

Ligand Modification-Free Methods for the Profiling of Protein-Environmental Chemical Interactions

Adverse health outcomes caused by environmental chemicals are often initiated via their interactions with proteins. Essentially, one environmental chemical may interact with a number of proteins and/or a protein may interact with a multitude of environmental chemicals, forming an intricate interaction network. Omics-wide protein-environmental chemical interaction profiling (PECI) is of prominent importance for comprehensive understanding of these interaction networks, including the toxicity mechanisms of action (MoA), and for providing systematic chemical safety assessment. However, such information remains unknown for most environmental chemicals, partly due to their vast chemical diversity. In recent years, with the continuous efforts afforded, especially in mass spectrometry (MS) based omics technologies, several ligand modification-free methods have been developed, and new attention for systematic PECI profiling was gained. In this Review, we provide a comprehensive overview on these methodologies for the identification of ligand-protein interactions, including affinity interaction-based methods of affinity-driven purification, covalent modification profiling, and activity-based protein profiling (ABPP) in a competitive mode, physicochemical property changes assessment methods of ligand-directed nuclear magnetic resonance (ligand-directed NMR), MS integrated with equilibrium dialysis for the discovery of allostery systematically (MIDAS), thermal proteome profiling (TPP), limited proteolysis-coupled mass spectrometry (LiP-MS), stability of proteins from rates of oxidation (SPROX), and several intracellular downstream response characterization methods. We expect that the applications of these ligand modification-free technologies will drive a considerable increase in the number of PECI identified, facilitate unveiling the toxicological mechanisms, and ultimately contribute to systematic health risk assessment of environmental chemicals.

CDS-DB, an omnibus for patient-derived gene expression signatures induced by cancer treatment

Patient-derived gene expression signatures induced by cancer treatment, obtained from paired pre- and post-treatment clinical transcriptomes, can help reveal drug mechanisms of action (MOAs) in cancer patients and understand the molecular response mechanism of tumor sensitivity or resistance. Their integration and reuse may bring new insights. Paired pre- and post-treatment clinical transcriptomic data are rapidly accumulating. However, a lack of systematic collection makes data access, integration, and reuse challenging. We therefore present the Cancer Drug-induced gene expression Signature DataBase (CDS-DB). CDS-DB has collected 78 patient-derived, paired pre- and post-treatment transcriptomic source datasets with uniformly reprocessed expression profiles and manually curated metadata such as drug administration dosage, sampling time and location, and intrinsic drug response status. From these source datasets, 2012 patient-level gene perturbation signatures were obtained, covering 85 therapeutic regimens, 39 cancer subtypes and 3628 patient samples. Besides data browsing, download and search, CDS-DB also supports single signature analysis (including differential gene expression, functional enrichment, tumor microenvironment and correlation analyses), signature comparative analysis and signature connectivity analysis. This provides insights into drug MOA and its heterogeneity in patients, drug resistance mechanisms, drug repositioning and drug (combination) discovery, etc. CDS-DB is available at http://cdsdb.ncpsb.org.cn/.

Illuminating phenotypic drug responses of sarcoma cells to kinase inhibitors by phosphoproteomics

Kinase inhibitors (KIs) are important cancer drugs but often feature polypharmacology that is molecularly not understood. This disconnect is particularly apparent in cancer entities such as sarcomas for which the oncogenic drivers are often not clear. To investigate more systematically how the cellular proteotypes of sarcoma cells shape their response to molecularly targeted drugs, we profiled the proteomes and phosphoproteomes of 17 sarcoma cell lines and screened the same against 150 cancer drugs. The resulting 2550 phenotypic profiles revealed distinct drug responses and the cellular activity landscapes derived from deep (phospho)proteomes (9-10,000 proteins and 10-27,000 phosphorylation sites per cell line) enabled several lines of analysis. For instance, connecting the (phospho)proteomic data with drug responses revealed known and novel mechanisms of action (MoAs) of KIs and identified markers of drug sensitivity or resistance. All data is publicly accessible via an interactive web application that enables exploration of this rich molecular resource for a better understanding of active signalling pathways in sarcoma cells, identifying treatment response predictors and revealing novel MoA of clinical KIs.

A proteomic landscape of pharmacologic perturbations for functional relevance

Pharmacological perturbation studies based on protein-level signatures are fundamental for drug discovery. In the present study, we used a mass spectrometry (MS)-based proteomic platform to profile the whole proteome of the breast cancer MCF7 cell line under stress induced by 78 bioactive compounds. The integrated analysis of perturbed signal abundance revealed the connectivity between phenotypic behaviors and molecular features in cancer cells. Our data showed functional relevance in exploring the novel pharmacological activity of phenolic xanthohumol, as well as the noncanonical targets of clinically approved tamoxifen, lovastatin, and their derivatives. Furthermore, the rational design of synergistic inhibition using a combination of histone methyltransferase and topoisomerase was identified based on their complementary drug fingerprints. This study provides rich resources for the proteomic landscape of drug responses for precision therapeutic medicine.

Proteome-Wide Profiling Using Sample Multiplexing of a Human Cell Line Treated with Cannabidiol (CBD) and Tetrahydrocannabinol (THC)

Cannabis has been used historically for both medicinal and recreational purposes, with the most notable cannabinoids being cannabidiol (CBD) and tetrahydrocannabinol (THC). Although their therapeutic effects have been well studied and their recreational use is highly debated, the underlying mechanisms of their biological effects remain poorly defined. In this study, we use isobaric tag-based sample multiplexed proteome profiling to investigate protein abundance differences in the human neuroblastoma SH-SY5Y cell line treated with CBD and THC. We identified significantly regulated proteins by each treatment and performed a pathway classification and associated protein-protein interaction analysis. Our findings suggest that these treatments may lead to mitochondrial dysfunction and induce endoplasmic reticulum stress. These data can potentially be interrogated further to investigate the potential role of CBD and THC in various biological and disease contexts, providing a foundation for future studies.

Deep Proteomic Compound Profiling with the Orbitrap Ascend Tribrid Mass Spectrometer Using Tandem Mass Tags and Real-Time Search

Tandem mass tags (TMT) and tribrid mass spectrometers are a powerful combination for high-throughput proteomics with high quantitative accuracy. Increasingly, this technology is being used to map the effects of drugs on the proteome. However, the depth of proteomic profiling is still limited by sensitivity and speed. The new Orbitrap Ascend mass spectrometer was designed to address these limitations with a combination of hardware and software improvements. We evaluated the performance of the Ascend in multiple contexts including deep proteomic profiling. We found that the Ascend exhibited increased sensitivity, yielding higher signal-to-noise ratios than the Orbitrap Eclipse with shorter injection times. As a result, higher numbers of peptides and proteins were identified and quantified, especially with low sample input. TMT measurements had significantly improved signal-to-noise ratios, improving quantitative precision. In a fractionated 16plex sample that profiled proteomic differences across four human cell lines, the Ascend was able to quantify hundreds more proteins than the Eclipse, many of them low-abundant proteins, and the Ascend was able to quantify >8000 proteins in 30% less instrument time. We used the Ascend to analyze 8881 proteins in HCT116 cancer cells treated with covalent sulfolane/sulfolene inhibitors of peptidyl-prolyl cis-trans isomerase NIMA-interacting 1 (PIN1), a phosphorylation-specific peptidyl-prolyl cis-trans isomerase implicated in several cancers. We characterized these PIN1 inhibitors' effects on the proteome and identified discrepancies among the different compounds, which will facilitate a better understanding of the structure-activity relationship of this class of compounds. The Ascend was able to quantify statistically significant, potentially therapeutically relevant changes in proteins that the Eclipse could not detect.

Screening for chemicals with thyroid hormone-disrupting effects using zebrafish embryo

Thyroid disruption is an increasingly recognized issue in the use and development of chemicals and new drugs, especially to help toxicologist to complement the reproductive and developmental toxicology information of chemicals. Still, adequate assessment methods are scarce and often suffer a trade-off between physiological relevance and labor- and cost-intensive assays. Here, we present a tiered approach for a medium-throughput screening of chemicals to identify their thyroid disrupting potential in zebrafish embryos as a New Approach Methodology (NAM). After identifying the maximum tolerated concentrations, we exposed zebrafish larvae to sub-adverse effect levels of the reference compounds benzophenone-2, bisphenol A, phenylthiourea, potassium perchlorate, propylthiouracil, and phloroglucinol to exclude any systemic toxicity. Applying the transgenic zebrafish line that carries a gene for the red fluorescence protein (Tg(tg:mCherry)) under the thyroglobulin promoter, we could identify the thyroid disrupting effects of the chemicals by a time and cost-effective image analysis measuring the fluorescence levels in the thyroid glands. Our observations could be confirmed by altered expression patterns of genes involved in the hypothalamus-pituitary-thyroid (HPT) axis. Finally, to anchor the observed thyroid disruption, we determined some changes in the Thyroid hormone levels of triiodothyronine (T3) and Thyroxine (T4) using a newly developed liquid chromatography mass spectrometric (LCMS) method. The presented approach carries the potential to extend the toolbox for legislative authorities and chemical producers for the assessment of thyroid-specific endocrine disruption and to overcome current challenges in the evaluation of endocrine disruptors.

Real-Time Spectral Library Matching for Sample Multiplexed Quantitative Proteomics

Sample multiplexed quantitative proteomics assays have proved to be a highly versatile means to assay molecular phenotypes. Yet, stochastic precursor selection and precursor coisolation can dramatically reduce the efficiency of data acquisition and quantitative accuracy. To address this, intelligent data acquisition (IDA) strategies have recently been developed to improve instrument efficiency and quantitative accuracy for both discovery and targeted methods. Toward this end, we sought to develop and implement a new real-time spectral library searching (RTLS) workflow that could enable intelligent scan triggering and peak selection within milliseconds of scan acquisition. To ensure ease of use and general applicability, we built an application to read in diverse spectral libraries and file types from both empirical and predicted spectral libraries. We demonstrate that RTLS methods enable improved quantitation of multiplexed samples, particularly with consideration for quantitation from chimeric fragment spectra. We used RTLS to profile proteome responses to small molecule perturbations and were able to quantify up to 15% more significantly regulated proteins in half the gradient time compared to traditional methods. Taken together, the development of RTLS expands the IDA toolbox to improve instrument efficiency and quantitative accuracy for sample multiplexed analyses.

Enrichr-KG: bridging enrichment analysis across multiple libraries

Gene and protein set enrichment analysis is a critical step in the analysis of data collected from omics experiments. Enrichr is a popular gene set enrichment analysis web-server search engine that contains hundreds of thousands of annotated gene sets. While Enrichr has been useful in providing enrichment analysis with many gene set libraries from different categories, integrating enrichment results across libraries and domains of knowledge can further hypothesis generation. To this end, Enrichr-KG is a knowledge graph database and a web-server application that combines selected gene set libraries from Enrichr for integrative enrichment analysis and visualization. The enrichment results are presented as subgraphs made of nodes and links that connect genes to their enriched terms. In addition, users of Enrichr-KG can add gene-gene links, as well as predicted genes to the subgraphs. This graphical representation of cross-library results with enriched and predicted genes can illuminate hidden associations between genes and annotated enriched terms from across datasets and resources. Enrichr-KG currently serves 26 gene set libraries from different categories that include transcription, pathways, ontologies, diseases/drugs, and cell types. To demonstrate the utility of Enrichr-KG we provide several case studies. Enrichr-KG is freely available at: https://maayanlab.cloud/enrichr-kg.

Curcuminoids as Anticancer Drugs: Pleiotropic Effects, Potential for Metabolic Reprogramming and Prospects for the Future

The number of published studies on curcuminoids in cancer research, including its lead molecule curcumin and synthetic analogs, has been increasing substantially during the past two decades. Insights on the diversity of inhibitory effects they have produced on a multitude of pathways involved in carcinogenesis and tumor progression have been provided. As this wealth of data was obtained in settings of various experimental and clinical data, this review first aimed at presenting a chronology of discoveries and an update on their complex in vivo effects. Secondly, there are many interesting questions linked to their pleiotropic effects. One of them, a growing research topic, relates to their ability to modulate metabolic reprogramming. This review will also cover the use of curcuminoids as chemosensitizing molecules that can be combined with several anticancer drugs to reverse the phenomenon of multidrug resistance. Finally, current investigations in these three complementary research fields raise several important questions that will be put among the prospects for the future research related to the importance of these molecules in cancer research.

Comprehensive Profiling of Rapamycin Interacting Proteins with Multiple Mass Spectrometry-Based Omics Techniques

Profiling drug-protein interactions is critical for understanding a drug's mechanism of action and predicting the possible adverse side effects. However, to comprehensively profile drug-protein interactions remains a challenge. To address this issue, we proposed a strategy that integrates multiple mass spectrometry-based omics analysis to provided global drug-protein interactions, including physical interactions and functional interactions, with rapamycin (Rap) as a model. Chemoproteomics profiling reveals 47 Rap binding proteins including the known target protein FKBP12 with high confidence. Gen Ontology enrichment analysis suggested that the Rap binding proteins are implicated in several important cellular processes, such as DNA replication, immunity, autophagy, programmed cell death, aging, transcription modulation, vesicle-mediated transport, membrane organization, and carbohydrate and nucleobase metabolic processes. The phosphoproteomics profiling revealed 255 down-regulated and 150 up-regulated phosphoproteins responding to Rap stimulation; they mainly involve the PI3K-Akt-mTORC1 signaling axis. Untargeted metabolomic profiling revealed 22 down-regulated metabolites and 75 up-regulated metabolites responding to Rap stimulation; they are mainly associated with the synthesis processes of pyrimidine and purine. The integrative multiomics data analysis provides deep insight into the drug-protein interactions and reveals Rap's complicated mechanism of action.