In situ identification of cellular drug targets in mammalian tissue.

The lack of tools to observe drug-target interactions at cellular resolution in intact tissue has been a major barrier to understanding in vivo drug actions. Here, we develop clearing-assisted tissue click chemistry (CATCH) to optically image covalent drug targets in intact mammalian tissues. CATCH permits specific and robust in situ fluorescence imaging of target-bound drug molecules at subcellular resolution and enables the identification of target cell types. Using well-established inhibitors of endocannabinoid hydrolases and monoamine oxidases, direct or competitive CATCH not only reveals distinct anatomical distributions and predominant cell targets of different drug compounds in the mouse brain but also uncovers unexpected differences in drug engagement across and within brain regions, reflecting rare cell types, as well as dose-dependent target shifts across tissue, cellular, and subcellular compartments that are not accessible by conventional methods. CATCH represents a valuable platform for visualizing in vivo interactions of small molecules in tissue.

Integrative Proteomic Characterization of Human Lung Adenocarcinoma.

Genomic studies of lung adenocarcinoma (LUAD) have advanced our understanding of the disease's biology and accelerated targeted therapy. However, the proteomic characteristics of LUAD remain poorly understood. We carried out a comprehensive proteomics analysis of 103 cases of LUAD in Chinese patients. Integrative analysis of proteome, phosphoproteome, transcriptome, and whole-exome sequencing data revealed cancer-associated characteristics, such as tumor-associated protein variants, distinct proteomics features, and clinical outcomes in patients at an early stage or with EGFR and TP53 mutations. Proteome-based stratification of LUAD revealed three subtypes (S-I, S-II, and S-III) related to different clinical and molecular features. Further, we nominated potential drug targets and validated the plasma protein level of HSP 90ß as a potential prognostic biomarker for LUAD in an independent cohort. Our integrative proteomics analysis enables a more comprehensive understanding of the molecular landscape of LUAD and offers an opportunity for more precise diagnosis and treatment.

Functional Profiling of a Plasmodium Genome Reveals an Abundance of Essential Genes.

The genomes of malaria parasites contain many genes of unknown function. To assist drug development through the identification of essential genes and pathways, we have measured competitive growth rates in mice of 2,578 barcoded Plasmodium berghei knockout mutants, representing >50% of the genome, and created a phenotype database. At a single stage of its complex life cycle, P. berghei requires two-thirds of genes for optimal growth, the highest proportion reported from any organism and a probable consequence of functional optimization necessitated by genomic reductions during the evolution of parasitism. In contrast, extreme functional redundancy has evolved among expanded gene families operating at the parasite-host interface. The level of genetic redundancy in a single-celled organism may thus reflect the degree of environmental variation it experiences. In the case of Plasmodium parasites, this helps rationalize both the relative successes of drugs and the greater difficulty of making an effective vaccine.

Unbiased interrogation of functional lysine residues in human proteome.

CRISPR screens have empowered the high-throughput dissection of gene functions; however, more explicit genetic elements, such as codons of amino acids, require thorough interrogation. Here, we establish a CRISPR strategy for unbiasedly probing functional amino acid residues at the genome scale. By coupling adenine base editors and barcoded sgRNAs, we target 215,689 out of 611,267 (35%) lysine codons, involving 85% of the total protein-coding genes. We identify 1,572 lysine codons whose mutations perturb human cell fitness, with many of them implicated in cancer. These codons are then mirrored to gene knockout screen data to provide functional insights into the role of lysine residues in cellular fitness. Mining these data, we uncover a CUL3-centric regulatory network in which lysine residues of CUL3 CRL complex proteins control cell fitness by specifying protein-protein interactions. Our study offers a general strategy for interrogating genetic elements and provides functional insights into the human proteome.

Pharmaceutical-Grade Rigosertib Is a Microtubule-Destabilizing Agent.

We recently used CRISPRi/a-based chemical-genetic screens and cell biological, biochemical, and structural assays to determine that rigosertib, an anti-cancer agent in phase III clinical trials, kills cancer cells by destabilizing microtubules. Reddy and co-workers (Baker et al., 2020, this issue of Molecular Cell) suggest that a contaminating degradation product in commercial formulations of rigosertib is responsible for the microtubule-destabilizing activity. Here, we demonstrate that cells treated with pharmaceutical-grade rigosertib (>99.9% purity) or commercially obtained rigosertib have qualitatively indistinguishable phenotypes across multiple assays. The two formulations have indistinguishable chemical-genetic interactions with genes that modulate microtubule stability, both destabilize microtubules in cells and in vitro, and expression of a rationally designed tubulin mutant with a mutation in the rigosertib binding site (L240F TUBB) allows cells to proliferate in the presence of either formulation. Importantly, the specificity of the L240F TUBB mutant for microtubule-destabilizing agents has been confirmed independently. Thus, rigosertib kills cancer cells by destabilizing microtubules, in agreement with our original findings.

Deciphering Drug Targets and Actions with Single-Cell and Spatial Resolution.

Recent advances in chemical, molecular, and genetic approaches have provided us with an unprecedented capacity to identify drug-target interactions across the whole proteome and genome. Meanwhile, rapid developments of single-cell and spatial omics technologies are revolutionizing our understanding of the molecular architecture of biological systems. However, a significant gap remains in how we align our understanding of drug actions, traditionally based on molecular affinities, with the in vivo cellular and spatial tissue heterogeneity revealed by these newer techniques. Here, we review state-of-the-art methods for profiling drug-target interactions and emerging multiomics tools to delineate the tissue heterogeneity at single-cell resolution. Highlighting the recent technical advances enabling high-resolution, multiplexable in situ small-molecule drug imaging (clearing-assisted tissue click chemistry, or CATCH), we foresee the integration of single-cell and spatial omics platforms, data, and concepts into the future framework of defining and understanding in vivo drug-target interactions and mechanisms of actions.

Uncovering dark matter in cancer by identifying epigenetic drivers.

The complex relationship between chromatin accessibility, transcriptional regulation, and cancer transitions presents a daunting puzzle. Terekhanova et al. created a pan-cancer epigenetic and transcriptomic atlas at single-cell resolution, yielding important insights into the underlying chromatin architecture of cancer transitions and novel discoveries with the potential to advance precision medicine.

MSI-DTI: predicting drug-target interaction based on multi-source information and multi-head self-attention.

Identifying drug-target interactions (DTIs) holds significant importance in drug discovery and development, playing a crucial role in various areas such as virtual screening, drug repurposing and identification of potential drug side effects. However, existing methods commonly exploit only a single type of feature from drugs and targets, suffering from miscellaneous challenges such as high sparsity and cold-start problems. We propose a novel framework called MSI-DTI (Multi-Source Information-based Drug-Target Interaction Prediction) to enhance prediction performance, which obtains feature representations from different views by integrating biometric features and knowledge graph representations from multi-source information. Our approach involves constructing a Drug-Target Knowledge Graph (DTKG), obtaining multiple feature representations from diverse information sources for SMILES sequences and amino acid sequences, incorporating network features from DTKG and performing an effective multi-source information fusion. Subsequently, we employ a multi-head self-attention mechanism coupled with residual connections to capture higher-order interaction information between sparse features while preserving lower-order information. Experimental results on DTKG and two benchmark datasets demonstrate that our MSI-DTI outperforms several state-of-the-art DTIs prediction methods, yielding more accurate and robust predictions. The source codes and datasets are publicly accessible at https://github.com/KEAML-JLU/MSI-DTI.

The current state of molecular testing in the treatment of patients with solid tumors, 2019.

The world of molecular profiling has undergone revolutionary changes over the last few years as knowledge, technology, and even standard clinical practice have evolved. Broad molecular profiling is now nearly essential for all patients with metastatic solid tumors. New agents have been approved based on molecular testing instead of tumor site of origin. Molecular profiling methodologies have likewise changed such that tests that were performed on patients a few years ago are no longer complete and possibly inaccurate today. As with all rapid change, medical providers can quickly fall behind or struggle to find up-to-date sources to ensure he or she provides optimum care. In this review, the authors provide the current state of the art for molecular profiling/precision medicine, practice standards, and a view into the future ahead.

Structural insights into the interaction of three Y-shaped ligands with PI3Kα.

Class IA phosphoinositide 3-kinase alpha (PI3Kα) is an important drug target because it is one of the most frequently mutated proteins in human cancers. However, small molecule inhibitors currently on the market or under development have safety concerns due to a lack of selectivity. Therefore, other chemical scaffolds or unique mechanisms of catalytic kinase inhibition are needed. Here, we report the cryo-electron microscopy structures of wild-type PI3Kα, the dimer of p110α and p85α, in complex with three Y-shaped ligands [cpd16 (compound 16), cpd17 (compound 17), and cpd18 (compound 18)] of different affinities and no inhibitory effect on the kinase activity. Unlike ATP-competitive inhibitors, cpd17 adopts a Y-shaped conformation with one arm inserted into a binding pocket formed by R770 and W780 and the other arm lodged in the ATP-binding pocket at an angle that is different from that of the ATP phosphate tail. Such a special interaction induces a conformation of PI3Kα resembling that of the unliganded protein. These observations were confirmed with two isomers (cpd16 and cpd18). Further analysis of these Y-shaped ligands revealed the structural basis of differential binding affinities caused by stereo- or regiochemical modifications. Our results may offer a different direction toward the design of therapeutic agents against PI3Kα.

Contrastive learning in protein language space predicts interactions between drugs and protein targets.

Sequence-based prediction of drug-target interactions has the potential to accelerate drug discovery by complementing experimental screens. Such computational prediction needs to be generalizable and scalable while remaining sensitive to subtle variations in the inputs. However, current computational techniques fail to simultaneously meet these goals, often sacrificing performance of one to achieve the others. We develop a deep learning model, ConPLex, successfully leveraging the advances in pretrained protein language models ("PLex") and employing a protein-anchored contrastive coembedding ("Con") to outperform state-of-the-art approaches. ConPLex achieves high accuracy, broad adaptivity to unseen data, and specificity against decoy compounds. It makes predictions of binding based on the distance between learned representations, enabling predictions at the scale of massive compound libraries and the human proteome. Experimental testing of 19 kinase-drug interaction predictions validated 12 interactions, including four with subnanomolar affinity, plus a strongly binding EPHB1 inhibitor (KD = 1.3 nM). Furthermore, ConPLex embeddings are interpretable, which enables us to visualize the drug-target embedding space and use embeddings to characterize the function of human cell-surface proteins. We anticipate that ConPLex will facilitate efficient drug discovery by making highly sensitive in silico drug screening feasible at the genome scale. ConPLex is available open source at https://ConPLex.csail.mit.edu.

Robust use of phenotypic heterogeneity at drug target genes for mechanistic insights: Application of cis-multivariable Mendelian randomization to GLP1R gene region.

Phenotypic heterogeneity at genomic loci encoding drug targets can be exploited by multivariable Mendelian randomization to provide insight into the pathways by which pharmacological interventions may affect disease risk. However, statistical inference in such investigations may be poor if overdispersion heterogeneity in measured genetic associations is unaccounted for. In this work, we first develop conditional F statistics for dimension-reduced genetic associations that enable more accurate measurement of phenotypic heterogeneity. We then develop a novel extension for two-sample multivariable Mendelian randomization that accounts for overdispersion heterogeneity in dimension-reduced genetic associations. Our empirical focus is to use genetic variants in the GLP1R gene region to understand the mechanism by which GLP1R agonism affects coronary artery disease (CAD) risk. Colocalization analyses indicate that distinct variants in the GLP1R gene region are associated with body mass index and type 2 diabetes (T2D). Multivariable Mendelian randomization analyses that were corrected for overdispersion heterogeneity suggest that bodyweight lowering rather than T2D liability lowering effects of GLP1R agonism are more likely contributing to reduced CAD risk. Tissue-specific analyses prioritized brain tissue as the most likely to be relevant for CAD risk, of the tissues considered. We hope the multivariable Mendelian randomization approach illustrated here is widely applicable to better understand mechanisms linking drug targets to diseases outcomes, and hence to guide drug development efforts.

Drug Target Identification in Tissues by Thermal Proteome Profiling.

Drug target deconvolution can accelerate the drug discovery process by identifying a drug's targets (facilitating medicinal chemistry efforts) and off-targets (anticipating toxicity effects or adverse drug reactions). Multiple mass spectrometry-based approaches have been developed for this purpose, but thermal proteome profiling (TPP) remains to date the only one that does not require compound modification and can be used to identify intracellular targets in living cells. TPP is based on the principle that the thermal stability of a protein can be affected by its interactions. Recent developments of this approach have expanded its applications beyond drugs and cell cultures to studying protein-drug interactions and biological phenomena in tissues. These developments open up the possibility of studying drug treatment or mechanisms of disease in a holistic fashion, which can result in the design of better drugs and lead to a better understanding of fundamental biology.

Identification and single-base gene-editing functional validation of a cis-EPO variant as a genetic predictor for EPO-increasing therapies.

Hypoxia-inducible factor prolyl hydroxylase inhibitors (HIF-PHIs) are currently under clinical development for treating anemia in chronic kidney disease (CKD), but it is important to monitor their cardiovascular safety. Genetic variants can be used as predictors to help inform the potential risk of adverse effects associated with drug treatments. We therefore aimed to use human genetics to help assess the risk of adverse cardiovascular events associated with therapeutically altered EPO levels to help inform clinical trials studying the safety of HIF-PHIs. By performing a genome-wide association meta-analysis of EPO (n = 6,127), we identified a cis-EPO variant (rs1617640) lying in the EPO promoter region. We validated this variant as most likely causal in controlling EPO levels by using genetic and functional approaches, including single-base gene editing. Using this variant as a partial predictor for therapeutic modulation of EPO and large genome-wide association data in Mendelian randomization tests, we found no evidence (at p < 0.05) that genetically predicted long-term rises in endogenous EPO, equivalent to a 2.2-unit increase, increased risk of coronary artery disease (CAD, OR [95% CI] = 1.01 [0.93, 1.07]), myocardial infarction (MI, OR [95% CI] = 0.99 [0.87, 1.15]), or stroke (OR [95% CI] = 0.97 [0.87, 1.07]). We could exclude increased odds of 1.15 for cardiovascular disease for a 2.2-unit EPO increase. A combination of genetic and functional studies provides a powerful approach to investigate the potential therapeutic profile of EPO-increasing therapies for treating anemia in CKD.

SAM-DTA: a sequence-agnostic model for drug-target binding affinity prediction.

Drug-target binding affinity prediction is a fundamental task for drug discovery and has been studied for decades. Most methods follow the canonical paradigm that processes the inputs of the protein (target) and the ligand (drug) separately and then combines them together. In this study we demonstrate, surprisingly, that a model is able to achieve even superior performance without access to any protein-sequence-related information. Instead, a protein is characterized completely by the ligands that it interacts. Specifically, we treat different proteins separately, which are jointly trained in a multi-head manner, so as to learn a robust and universal representation of ligands that is generalizable across proteins. Empirical evidences show that the novel paradigm outperforms its competitive sequence-based counterpart, with the Mean Squared Error (MSE) of 0.4261 versus 0.7612 and the R-Square of 0.7984 versus 0.6570 compared with DeepAffinity. We also investigate the transfer learning scenario where unseen proteins are encountered after the initial training, and the cross-dataset evaluation for prospective studies. The results reveals the robustness of the proposed model in generalizing to unseen proteins as well as in predicting future data. Source codes and data are available at https://github.com/huzqatpku/SAM-DTA.

Improved drug response prediction by drug target data integration via network-based profiling.

Drug response prediction (DRP) is important for precision medicine to predict how a patient would react to a drug before administration. Existing studies take the cell line transcriptome data, and the chemical structure of drugs as input and predict drug response as IC50 or AUC values. Intuitively, use of drug target interaction (DTI) information can be useful for DRP. However, use of DTI is difficult because existing drug response database such as CCLE and GDSC do not have information about transcriptome after drug treatment. Although transcriptome after drug treatment is not available, if we can compute the perturbation effects by the pharmacologic modulation of target gene, we can utilize the DTI information in CCLE and GDSC. In this study, we proposed a framework that can improve existing deep learning-based DRP models by effectively utilizing drug target information. Our framework includes NetGP, a module to compute gene perturbation scores by the network propagation technique on a network. NetGP produces genes in a ranked list in terms of gene perturbation scores and the ranked genes are input to a multi-layer perceptron to generate a fixed dimension vector for the integration with existing DRP models. This integration is done in a model-agnostic way so that any existing DRP tool can be incorporated. As a result, our framework boosts the performance of existing DRP models, in 64 of 72 comparisons. The performance gains are larger especially for test scenarios with samples with unseen drugs by large margins up to 34% in Pearson's correlation coefficient.

MHTAN-DTI: Metapath-based hierarchical transformer and attention network for drug-target interaction prediction.

Drug-target interaction (DTI) prediction can identify novel ligands for specific protein targets, and facilitate the rapid screening of effective new drug candidates to speed up the drug discovery process. However, the current methods are not sensitive enough to complex topological structures, and complicated relations between multiple node types are not fully captured yet. To address the above challenges, we construct a metapath-based heterogeneous bioinformatics network, and then propose a DTI prediction method with metapath-based hierarchical transformer and attention network for drug-target interaction prediction (MHTAN-DTI), applying metapath instance-level transformer, single-semantic attention and multi-semantic attention to generate low-dimensional vector representations of drugs and proteins. Metapath instance-level transformer performs internal aggregation on the metapath instances, and models global context information to capture long-range dependencies. Single-semantic attention learns the semantics of a certain metapath type, introduces the central node weight and assigns different weights to different metapath instances to obtain the semantic-specific node embedding. Multi-semantic attention captures the importance of different metapath types and performs weighted fusion to attain the final node embedding. The hierarchical transformer and attention network weakens the influence of noise data on the DTI prediction results, and enhances the robustness and generalization ability of MHTAN-DTI. Compared with the state-of-the-art DTI prediction methods, MHTAN-DTI achieves significant performance improvements. In addition, we also conduct sufficient ablation studies and visualize the experimental results. All the results demonstrate that MHTAN-DTI can offer a powerful and interpretable tool for integrating heterogeneous information to predict DTIs and provide new insights into drug discovery.

BindingSite-AugmentedDTA: enabling a next-generation pipeline for interpretable prediction models in drug repurposing.

While research into drug-target interaction (DTI) prediction is fairly mature, generalizability and interpretability are not always addressed in the existing works in this field. In this paper, we propose a deep learning (DL)-based framework, called BindingSite-AugmentedDTA, which improves drug-target affinity (DTA) predictions by reducing the search space of potential-binding sites of the protein, thus making the binding affinity prediction more efficient and accurate. Our BindingSite-AugmentedDTA is highly generalizable as it can be integrated with any DL-based regression model, while it significantly improves their prediction performance. Also, unlike many existing models, our model is highly interpretable due to its architecture and self-attention mechanism, which can provide a deeper understanding of its underlying prediction mechanism by mapping attention weights back to protein-binding sites. The computational results confirm that our framework can enhance the prediction performance of seven state-of-the-art DTA prediction algorithms in terms of four widely used evaluation metrics, including concordance index, mean squared error, modified squared correlation coefficient ($r^2_m$) and the area under the precision curve. We also contribute to three benchmark drug-traget interaction datasets by including additional information on 3D structure of all proteins contained in those datasets, which include the two most commonly used datasets, namely Kiba and Davis, as well as the data from IDG-DREAM drug-kinase binding prediction challenge. Furthermore, we experimentally validate the practical potential of our proposed framework through in-lab experiments. The relatively high agreement between computationally predicted and experimentally observed binding interactions supports the potential of our framework as the next-generation pipeline for prediction models in drug repurposing.

Metapath-aggregated heterogeneous graph neural network for drug-target interaction prediction.

Drug-target interaction (DTI) prediction is an essential step in drug repositioning. A few graph neural network (GNN)-based methods have been proposed for DTI prediction using heterogeneous biological data. However, existing GNN-based methods only aggregate information from directly connected nodes restricted in a drug-related or a target-related network and are incapable of capturing high-order dependencies in the biological heterogeneous graph. In this paper, we propose a metapath-aggregated heterogeneous graph neural network (MHGNN) to capture complex structures and rich semantics in the biological heterogeneous graph for DTI prediction. Specifically, MHGNN enhances heterogeneous graph structure learning and high-order semantics learning by modeling high-order relations via metapaths. Additionally, MHGNN enriches high-order correlations between drug-target pairs (DTPs) by constructing a DTP correlation graph with DTPs as nodes. We conduct extensive experiments on three biological heterogeneous datasets. MHGNN favorably surpasses 17 state-of-the-art methods over 6 evaluation metrics, which verifies its efficacy for DTI prediction. The code is available at https://github.com/Zora-LM/MHGNN-DTI.

Fine-grained selective similarity integration for drug-target interaction prediction.

The discovery of drug-target interactions (DTIs) is a pivotal process in pharmaceutical development. Computational approaches are a promising and efficient alternative to tedious and costly wet-lab experiments for predicting novel DTIs from numerous candidates. Recently, with the availability of abundant heterogeneous biological information from diverse data sources, computational methods have been able to leverage multiple drug and target similarities to boost the performance of DTI prediction. Similarity integration is an effective and flexible strategy to extract crucial information across complementary similarity views, providing a compressed input for any similarity-based DTI prediction model. However, existing similarity integration methods filter and fuse similarities from a global perspective, neglecting the utility of similarity views for each drug and target. In this study, we propose a Fine-Grained Selective similarity integration approach, called FGS, which employs a local interaction consistency-based weight matrix to capture and exploit the importance of similarities at a finer granularity in both similarity selection and combination steps. We evaluate FGS on five DTI prediction datasets under various prediction settings. Experimental results show that our method not only outperforms similarity integration competitors with comparable computational costs, but also achieves better prediction performance than state-of-the-art DTI prediction approaches by collaborating with conventional base models. Furthermore, case studies on the analysis of similarity weights and on the verification of novel predictions confirm the practical ability of FGS.