Computational Characterization of Membrane Proteins as Anticancer Targets: Current Challenges and Opportunities.

Cancer remains a leading cause of mortality worldwide and calls for novel therapeutic targets. Membrane proteins are key players in various cancer types but present unique challenges compared to soluble proteins. The advent of computational drug discovery tools offers a promising approach to address these challenges, allowing for the prioritization of "wet-lab" experiments. In this review, we explore the applications of computational approaches in membrane protein oncological characterization, particularly focusing on three prominent membrane protein families: receptor tyrosine kinases (RTKs), G protein-coupled receptors (GPCRs), and solute carrier proteins (SLCs). We chose these families due to their varying levels of understanding and research data availability, which leads to distinct challenges and opportunities for computational analysis. We discuss the utilization of multi-omics data, machine learning, and structure-based methods to investigate aberrant protein functionalities associated with cancer progression within each family. Moreover, we highlight the importance of considering the broader cellular context and, in particular, cross-talk between proteins. Despite existing challenges, computational tools hold promise in dissecting membrane protein dysregulation in cancer. With advancing computational capabilities and data resources, these tools are poised to play a pivotal role in identifying and prioritizing membrane proteins as personalized anticancer targets.

Still in Search for an EAAT Activator: GT949 Does Not Activate EAAT2, nor EAAT3 in Impedance and Radioligand Uptake Assays.

Excitatory amino acid transporters (EAATs) are important regulators of amino acid transport and in particular glutamate. Recently, more interest has arisen in these transporters in the context of neurodegenerative diseases. This calls for ways to modulate these targets to drive glutamate transport, EAAT2 and EAAT3 in particular. Several inhibitors (competitive and noncompetitive) exist to block glutamate transport; however, activators remain scarce. Recently, GT949 was proposed as a selective activator of EAAT2, as tested in a radioligand uptake assay. In the presented research, we aimed to validate the use of GT949 to activate EAAT2-driven glutamate transport by applying an innovative, impedance-based, whole-cell assay (xCELLigence). A broad range of GT949 concentrations in a variety of cellular environments were tested in this assay. As expected, no activation of EAAT3 could be detected. Yet, surprisingly, no biological activation of GT949 on EAAT2 could be observed in this assay either. To validate whether the impedance-based assay was not suited to pick up increased glutamate uptake or if the compound might not induce activation in this setup, we performed radioligand uptake assays. Two setups were utilized; a novel method compared to previously published research, and in a reproducible fashion copying the methods used in the existing literature. Nonetheless, activation of neither EAAT2 nor EAAT3 could be observed in these assays. Furthermore, no evidence of GT949 binding or stabilization of purified EAAT2 could be observed in a thermal shift assay. To conclude, based on experimental evidence in the present study GT949 requires specific assay conditions, which are difficult to reproduce, and the compound cannot simply be classified as an activator of EAAT2 based on the presented evidence. Hence, further research is required to develop the tools needed to identify new EAAT modulators and use their potential as a therapeutic target.

Molecular insights into disease-associated glutamate transporter (EAAT1 / SLC1A3) variants using in silico and in vitro approaches.

Glutamate is an essential excitatory neurotransmitter and an intermediate for energy metabolism. Depending on the tumor site, cancer cells have increased or decreased expression of excitatory amino acid transporter 1 or 2 (EAAT1/2, SLC1A3/2) to regulate glutamate uptake for the benefit of tumor growth. Thus, EAAT1/2 may be an attractive target for therapeutic intervention in oncology. Genetic variation of EAAT1 has been associated with rare cases of episodic ataxia, but the occurrence and functional contribution of EAAT1 mutants in other diseases, such as cancer, is poorly understood. Here, 105 unique somatic EAAT1 mutations were identified in cancer patients from the Genomic Data Commons dataset. Using EAAT1 crystal structures and in silico studies, eight mutations were selected based on their close proximity to the orthosteric or allosteric ligand binding sites and the predicted change in ligand binding affinity. In vitro functional assessment in a live-cell, impedance-based phenotypic assay demonstrated that these mutants differentially affect L-glutamate and L-aspartate transport, as well as the inhibitory potency of an orthosteric (TFB-TBOA) and allosteric (UCPH-101) inhibitor. Moreover, two episodic ataxia-related mutants displayed functional responses that were in line with literature, which confirmed the validity of our assay. Of note, ataxia-related mutant M128R displayed inhibitor-induced functional responses never described before. Finally, molecular dynamics (MD) simulations were performed to gain mechanistic insights into the observed functional effects. Taken together, the results in this work demonstrate 1) the suitability of the label-free phenotypic method to assess functional variation of EAAT1 mutants and 2) the opportunity and challenges of using in silico techniques to rationalize the in vitro phenotype of disease-relevant mutants.

Artificial intelligence for natural product drug discovery.

Developments in computational omics technologies have provided new means to access the hidden diversity of natural products, unearthing new potential for drug discovery. In parallel, artificial intelligence approaches such as machine learning have led to exciting developments in the computational drug design field, facilitating biological activity prediction and de novo drug design for molecular targets of interest. Here, we describe current and future synergies between these developments to effectively identify drug candidates from the plethora of molecules produced by nature. We also discuss how to address key challenges in realizing the potential of these synergies, such as the need for high-quality datasets to train deep learning algorithms and appropriate strategies for algorithm validation.

3DDPDs: describing protein dynamics for proteochemometric bioactivity prediction. A case for (mutant) G protein-coupled receptors.

Proteochemometric (PCM) modelling is a powerful computational drug discovery tool used in bioactivity prediction of potential drug candidates relying on both chemical and protein information. In PCM features are computed to describe small molecules and proteins, which directly impact the quality of the predictive models. State-of-the-art protein descriptors, however, are calculated from the protein sequence and neglect the dynamic nature of proteins. This dynamic nature can be computationally simulated with molecular dynamics (MD). Here, novel 3D dynamic protein descriptors (3DDPDs) were designed to be applied in bioactivity prediction tasks with PCM models. As a test case, publicly available G protein-coupled receptor (GPCR) MD data from GPCRmd was used. GPCRs are membrane-bound proteins, which are activated by hormones and neurotransmitters, and constitute an important target family for drug discovery. GPCRs exist in different conformational states that allow the transmission of diverse signals and that can be modified by ligand interactions, among other factors. To translate the MD-encoded protein dynamics two types of 3DDPDs were considered: one-hot encoded residue-specific (rs) and embedding-like protein-specific (ps) 3DDPDs. The descriptors were developed by calculating distributions of trajectory coordinates and partial charges, applying dimensionality reduction, and subsequently condensing them into vectors per residue or protein, respectively. 3DDPDs were benchmarked on several PCM tasks against state-of-the-art non-dynamic protein descriptors. Our rs- and ps3DDPDs outperformed non-dynamic descriptors in regression tasks using a temporal split and showed comparable performance with a random split and in all classification tasks. Combinations of non-dynamic descriptors with 3DDPDs did not result in increased performance. Finally, the power of 3DDPDs to capture dynamic fluctuations in mutant GPCRs was explored. The results presented here show the potential of including protein dynamic information on machine learning tasks, specifically bioactivity prediction, and open opportunities for applications in drug discovery, including oncology.

Cancer-Associated Mutations of the Adenosine A2A Receptor Have Diverse Influences on Ligand Binding and Receptor Functions.

The adenosine A2A receptor (A2AAR) is a class A G-protein-coupled receptor (GPCR). It is an immune checkpoint in the tumor micro-environment and has become an emerging target for cancer treatment. In this study, we aimed to explore the effects of cancer-patient-derived A2AAR mutations on ligand binding and receptor functions. The wild-type A2AAR and 15 mutants identified by Genomic Data Commons (GDC) in human cancers were expressed in HEK293T cells. Firstly, we found that the binding affinity for agonist NECA was decreased in six mutants but increased for the V275A mutant. Mutations A165V and A265V decreased the binding affinity for antagonist ZM241385. Secondly, we found that the potency of NECA (EC50) in an impedance-based cell-morphology assay was mostly correlated with the binding affinity for the different mutants. Moreover, S132L and H278N were found to shift the A2AAR towards the inactive state. Importantly, we found that ZM241385 could not inhibit the activation of V275A and P285L stimulated by NECA. Taken together, the cancer-associated mutations of A2AAR modulated ligand binding and receptor functions. This study provides fundamental insights into the structure-activity relationship of the A2AAR and provides insights for A2AAR-related personalized treatment in cancer.

Oncological drug discovery: AI meets structure-based computational research.

The integration of machine learning and structure-based methods has proven valuable in the past as a way to prioritize targets and compounds in early drug discovery. In oncological research, these methods can be highly beneficial in addressing the diversity of neoplastic diseases portrayed by the different hallmarks of cancer. Here, we review six use case scenarios for integrated computational methods, namely driver prediction, computational mutagenesis, (off)-target prediction, binding site prediction, virtual screening, and allosteric modulation analysis. We address the heterogeneity of integration approaches and individual methods, while acknowledging their current limitations and highlighting their potential to bring drugs for personalized oncological therapies to the market faster.

Structural Insights into the Mechanisms of Action of Functionally Distinct Classes of Chikungunya Virus Nonstructural Protein 1 Inhibitors.

Chikungunya virus (CHIKV) nonstructural protein 1 (nsP1) harbors the methyltransferase (MTase) and guanylyltransferase (GTase) activities needed for viral RNA capping and represents a promising antiviral drug target. We compared the antiviral efficacies of nsP1 inhibitors belonging to the MADTP, CHVB, and FHNA series (6'-fluoro-homoneplanocin A [FHNA], its 3'-keto form, and 6'-ß-fluoro-homoaristeromycin). Cell-based phenotypic cross-resistance assays revealed that the CHVB and MADTP series had similar modes of action that differed from that of the FHNA series. In biochemical assays with purified Semliki Forest virus and CHIKV nsP1, CHVB compounds strongly inhibited MTase and GTase activities, while MADTP-372 had a moderate inhibitory effect. FHNA did not directly inhibit the enzymatic activity of CHIKV nsP1. The first-of-their-kind molecular-docking studies with the cryo-electron microscopy (cryo-EM) structure of CHIKV nsP1, which is assembled into a dodecameric ring, revealed that the MADTP and CHVB series bind at the S-adenosylmethionine (SAM)-binding site in the capping domain, where they would function as competitive or noncompetitive inhibitors. The FHNA series was predicted to bind at the secondary binding pocket in the ring-aperture membrane-binding and oligomerization (RAMBO) domain, potentially interfering with the membrane binding and oligomerization of nsP1. Our cell-based and enzymatic assays, in combination with molecular docking and mapping of compound resistance mutations to the nsP1 structure, allowed us to group nsP1 inhibitors into functionally distinct classes. This study identified druggable pockets in the nsP1 dodecameric structure and provides a basis for the rational design, optimization, and combination of inhibitors of this unique and promising antiviral drug target.