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.

Kinetic profiling of novel spirobenzo-oxazinepiperidinone derivatives as equilibrative nucleoside transporter 1 inhibitors.

Evaluation of kinetic parameters of drug-target binding, kon, koff, and residence time (RT), in addition to the traditional in vitro parameter of affinity is receiving increasing attention in the early stages of drug discovery. Target binding kinetics emerges as a meaningful concept for the evaluation of a ligand's duration of action and more generally drug efficacy and safety. We report the biological evaluation of a novel series of spirobenzo-oxazinepiperidinone derivatives as inhibitors of the human equilibrative nucleoside transporter 1 (hENT1, SLC29A1). The compounds were evaluated in radioligand binding experiments, i.e., displacement, competition association, and washout assays, to evaluate their affinity and binding kinetic parameters. We also linked these pharmacological parameters to the compounds' chemical characteristics, and learned that separate moieties of the molecules governed target affinity and binding kinetics. Among the 29 compounds tested, 28 stood out with high affinity and a long residence time of 87 min. These findings reveal the importance of supplementing affinity data with binding kinetics at transport proteins such as hENT1.

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.

Dual allosteric and orthosteric pharmacology of synthetic analog cannabidiol-dimethylheptyl, but not cannabidiol, on the cannabinoid CB2 receptor.

Cannabinoid CB2 receptor (CB2R) is a class A G protein-coupled receptor (GPCR) involved in a broad spectrum of physiological processes and pathological conditions. For that reason, targeting CB2R might provide therapeutic opportunities in neurodegenerative disorders, neuropathic pain, inflammatory diseases, and cancer. The main components from Cannabis sativa, such as Δ9-tetrahydrocannabinol (Δ9-THC) and cannabidiol (CBD), have been therapeutically exploited and synthetically-derived analogs have been generated. One example is cannabidiol-dimethylheptyl (CBD-DMH), which exhibits anti-inflammatory effects. Nevertheless, its pharmacological mechanism of action is not yet fully understood and is hypothesized for multiple targets, including CB2R. The aim of this study was to further investigate the molecular pharmacology of CBD-DMH on CB2R while CBD was taken along as control. These compounds were screened in equilibrium and kinetic radioligand binding studies and various functional assays, including G protein activation, inhibition of cAMP production and ß-arrestin-2 recruitment. In dissociation studies, CBD-DMH allosterically modulated the radioligand binding. Furthermore, CBD-DMH negatively modulated the G protein activation of reference agonists CP55,940, AEA and 2-AG, but not the agonist-induced ß-arrestin-2 recruitment. Nevertheless, CBD-DMH also displayed competitive binding to CB2R and partial agonism on G protein activation, inhibition of cAMP production and ß-arrestin-2 recruitment. CBD did not exhibit such allosteric behavior and only very weakly bound CB2R without activation. This study shows a dual binding mode of CBD-DMH, but not CBD, to CB2R with the suggestion of two different binding sites. Altogether, it encourages further research into this dual mechanism which might provide a new class of molecules targeting CB2R.

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.

Development of an Affinity-Based Probe to Profile Endogenous Human Adenosine A3 Receptor Expression.

The adenosine A3 receptor (A3AR) is a G protein-coupled receptor (GPCR) that exerts immunomodulatory effects in pathophysiological conditions such as inflammation and cancer. Thus far, studies toward the downstream effects of A3AR activation have yielded contradictory results, thereby motivating the need for further investigations. Various chemical and biological tools have been developed for this purpose, ranging from fluorescent ligands to antibodies. Nevertheless, these probes are limited by their reversible mode of binding, relatively large size, and often low specificity. Therefore, in this work, we have developed a clickable and covalent affinity-based probe (AfBP) to target the human A3AR. Herein, we show validation of the synthesized AfBP in radioligand displacement, SDS-PAGE, and confocal microscopy experiments as well as utilization of the AfBP for the detection of endogenous A3AR expression in flow cytometry experiments. Ultimately, this AfBP will aid future studies toward the expression and function of the A3AR in pathologies.

Microwave-Assisted Synthesis, Structure, and Preliminary Biological Evaluation of Novel 6-Methoxy-5,6-dihydro-5-azapurines.

We herein disclose the microwave-assisted synthesis of previously unreported 6-methoxy-5,6-dihydro-5-azapurines, whose purine-like scaffold is promising for drug discovery. The method is simple, fast, and relies on easily accessible reagents such as trimethyl orthoformate, acetic acid, and aminotriazole-derived N,N'-disubstituted formamidines. The preliminary biological evaluation revealed that selected representatives of synthesized 6-methoxy-5,6-dihydro-5-azapurines dose-dependently reduce the viability of HepG2 and A549 cancer cells having little to no influence on five tested purinergic receptors.

A Chemical Biological Approach to Study G Protein-Coupled Receptors: Labeling the Adenosine A1 Receptor Using an Electrophilic Covalent Probe.

G protein-coupled receptors (GPCRs) have been known for decades as attractive drug targets. This has led to the development and approval of many ligands targeting GPCRs. Although ligand binding effects have been studied thoroughly for many GPCRs, there are multiple aspects of GPCR signaling that remain poorly understood. The reasons for this are the difficulties that are encountered upon studying GPCRs, for example, a poor solubility and low expression levels. In this work, we have managed to overcome some of these issues by developing an affinity-based probe for a prototypic GPCR, the adenosine A1 receptor (A1AR). Here, we show the design, synthesis, and biological evaluation of this probe in various biochemical assays, such as SDS-PAGE, confocal microscopy, and chemical proteomics.

Development of subtype-selective covalent ligands for the adenosine A2B receptor by tuning the reactive group.

Signalling through the adenosine receptors (ARs), in particular through the adenosine A2B receptor (A2BAR), has been shown to play a role in a variety of pathological conditions, ranging from immune disorders to cancer. Covalent ligands for the A2BAR have the potential to irreversibly block the receptor, as well as inhibit all A2BAR-induced signalling pathways. This will allow a thorough investigation of the pathophysiological role of the receptor. In this study, we synthesized and evaluated a set of potential covalent ligands for the A2BAR. The ligands all contain a core scaffold consisting of a substituted xanthine, varying in type and orientation of electrophilic group (warhead). Here, we find that the right combination of these variables is necessary for a high affinity, irreversible mode of binding and selectivity towards the A2BAR. Altogether, this is the case for sulfonyl fluoride 24 (LUF7982), a covalent ligand that allows for novel ways to interrogate the A2BAR.

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.

Cancer-Related Somatic Mutations in Transmembrane Helices Alter Adenosine A1 Receptor Pharmacology.

Overexpression of the adenosine A1 receptor (A1AR) has been detected in various cancer cell lines. However, the role of A1AR in tumor development is still unclear. Thirteen A1AR mutations were identified in the Cancer Genome Atlas from cancer patient samples. We have investigated the pharmacology of the mutations located at the 7-transmembrane domain using a yeast system. Concentration-growth curves were obtained with the full agonist CPA and compared to the wild type hA1AR. H78L3.23 and S246T6.47 showed increased constitutive activity, while only the constitutive activity of S246T6.47 could be reduced to wild type levels by the inverse agonist DPCPX. Decreased constitutive activity was observed on five mutant receptors, among which A52V2.47 and W188C5.46 showed a diminished potency for CPA. Lastly, a complete loss of activation was observed in five mutant receptors. A selection of mutations was also investigated in a mammalian system, showing comparable effects on receptor activation as in the yeast system, except for residues pointing toward the membrane. Taken together, this study will enrich the view of the receptor structure and function of A1AR, enlightening the consequences of these mutations in cancer. Ultimately, this may provide an opportunity for precision medicine for cancer patients with pathological phenotypes involving these mutations.

Impedance-Based Phenotypic Readout of Transporter Function: A Case for Glutamate Transporters.

Excitatory amino acid transporters (EAAT/SLC1) mediate Na+-dependent uptake of extracellular glutamate and are potential drug targets for neurological disorders. Conventional methods to assess glutamate transport in vitro are based on radiolabels, fluorescent dyes or electrophysiology, which potentially compromise the cell's physiology and are generally less suited for primary drug screens. Here, we describe a novel label-free method to assess human EAAT function in living cells, i.e., without the use of chemical modifications to the substrate or cellular environment. In adherent HEK293 cells overexpressing EAAT1, stimulation with glutamate or aspartate induced cell spreading, which was detected in real-time using an impedance-based biosensor. This change in cell morphology was prevented in the presence of the Na+/K+-ATPase inhibitor ouabain and EAAT inhibitors, which suggests the substrate-induced response was ion-dependent and transporter-specific. A mechanistic explanation for the phenotypic response was substantiated by actin cytoskeleton remodeling and changes in the intracellular levels of the osmolyte taurine, which suggests that the response involves cell swelling. In addition, substrate-induced cellular responses were observed for cells expressing other EAAT subtypes, as well as in a breast cancer cell line (MDA-MB-468) with endogenous EAAT1 expression. These findings allowed the development of a label-free high-throughput screening assay, which could be beneficial in early drug discovery for EAATs and holds potential for the study of other transport proteins that modulate cell shape.

Cancer-related somatic mutations alter adenosine A1 receptor pharmacology-A focus on mutations in the loops and C-terminus.

G protein-coupled receptors (GPCRs) are known to be involved in tumor progression and metastasis. The adenosine A1 receptor (A1 AR) has been detected to be over-expressed in various cancer cell lines. However, the role of A1 AR in tumor development is not yet well characterized. A series of A1 AR mutations were identified in the Cancer Genome Atlas from cancer patient samples. In this study, we have investigated the pharmacology of mutations located outside of the 7-transmembrane domain by using a "single-GPCR-one-G protein" yeast system. Concentration-growth curves were obtained with the full agonist CPA for 12 mutant receptors and compared to the wild-type hA1 AR. Most mutations located at the extracellular loops (EL) reduced the levels of constitutive activity of the receptor and agonist potency. For mutants at the intracellular loops (ILs) of the receptor, an increased constitutive activity was found for mutant receptor L211R5.69 , while a decreased constitutive activity and agonist response were found for mutant receptor L113F34.51 . Lastly, mutations identified on the C-terminus did not significantly influence the pharmacological function of the receptor. A selection of mutations was also investigated in a mammalian system. Overall, similar effects on receptor activation compared to the yeast system were found with mutations located at the EL, but some contradictory effects were observed for mutations located at the IL. Taken together, this study will enrich the insight of A1 AR structure and function, enlightening the consequences of these mutations in cancer. Ultimately, this may provide potential precision medicine in cancer treatment.

Kinetic profiling and functional characterization of 8-phenylxanthine derivatives as A2B adenosine receptor antagonists.

A2B adenosine receptor (A2BAR) antagonists have therapeutic potential in inflammation-related diseases such as asthma, chronic obstructive pulmonary disease and cancer. However, no drug is currently clinically approved, creating a demand for research on novel antagonists. Over the last decade, the study of target binding kinetics, along with affinity and potency, has been proven valuable in early drug discovery stages, as it is associated with improved in vivo drug efficacy and safety. In this study, we report the synthesis and biological evaluation of a series of xanthine derivatives as A2BAR antagonists, including an isothiocyanate derivative designed to bind covalently to the receptor. All 28 final compounds were assessed in radioligand binding experiments, to evaluate their affinity and for those qualifying, kinetic binding parameters. Both structure-affinity and structure-kinetic relationships were derived, providing a clear relationship between affinity and dissociation rate constants. Two structurally similar compounds, 17 and 18, were further evaluated in a label-free assay due to their divergent kinetic profiles. An extended cellular response was associated with long A2BAR residence times. This link between a ligand's A2BAR residence time and its functional effect highlights the importance of binding kinetics as a selection parameter in the early stages of drug discovery.

Development of Bicyclo[3.1.0]hexane-Based A3 Receptor Ligands: Closing the Gaps in the Structure-Affinity Relationships.

The adenosine A3 receptor is a promising target for treating and diagnosing inflammation and cancer. In this paper, a series of bicyclo[3.1.0]hexane-based nucleosides was synthesized and evaluated for their P1 receptor affinities in radioligand binding studies. The study focused on modifications at 1-, 2-, and 6-positions of the purine ring and variations of the 5'-position at the bicyclo[3.1.0]hexane moiety, closing existing gaps in the structure-affinity relationships. The most potent derivative 30 displayed moderate A3AR affinity (Ki of 0.38 µM) and high A3R selectivity. A subset of compounds varied at 5'-position was further evaluated in functional P2Y1R assays, displaying no off-target activity.

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.

Deciphering conformational selectivity in the A2A adenosine G protein-coupled receptor by free energy simulations.

Transmembranal G Protein-Coupled Receptors (GPCRs) transduce extracellular chemical signals to the cell, via conformational change from a resting (inactive) to an active (canonically bound to a G-protein) conformation. Receptor activation is normally modulated by extracellular ligand binding, but mutations in the receptor can also shift this equilibrium by stabilizing different conformational states. In this work, we built structure-energetic relationships of receptor activation based on original thermodynamic cycles that represent the conformational equilibrium of the prototypical A2A adenosine receptor (AR). These cycles were solved with efficient free energy perturbation (FEP) protocols, allowing to distinguish the pharmacological profile of different series of A2AAR agonists with different efficacies. The modulatory effects of point mutations on the basal activity of the receptor or on ligand efficacies could also be detected. This methodology can guide GPCR ligand design with tailored pharmacological properties, or allow the identification of mutations that modulate receptor activation with potential clinical implications.

Piperazine squaric acid diamides, a novel class of allosteric P2X7 receptor antagonists.

The P2X7 receptor (P2X7R) stands out among the purinergic receptors due to its strong involvement in the regulation of tumor growth and metastasis formation as well as in innate immune responses and afferent signal transmission. Numerous studies have pointed out the beneficial effects of P2X7R antagonism for the treatment of a variety of cancer types, inflammatory diseases, and chronic pain. Herein we describe the development of novel P2X7R antagonists, incorporating piperazine squaric diamides as a central element. Besides improving the antagonists' potency from pIC50 values of 5.7-7.6, ADME properties (logD7.4 value, plasma protein binding, in vitro metabolic stability) of the generated compounds were investigated and optimized to provide novel P2X7R antagonists with drug-like properties. Furthermore, docking studies revealed the antagonists binding to the allosteric binding pocket in two distinct binding poses, depending on the substitution of the central piperazine moiety.

Anticancer opportunities at every stage of chemokine function.

The chemokine system, comprising 48 chemokines and 23 receptors, is critically involved in several hallmarks of cancer. Yet, despite extensive efforts from the pharmaceutical sector, only two drugs aimed at this system are currently approved for clinical use against cancer. To date, numerous pharmacological approaches have been developed to successfully intervene at different stages of chemokine function: (i) chemokine availability; (ii) chemokine-glycosaminoglycan binding; and (iii) chemokine receptor binding. Many of these strategies have been tested in preclinical cancer models, and some have advanced to clinical trials as potential anticancer therapies. Here we will review the strategies and growing pharmacological toolbox for manipulating the chemokine system in cancer, and address novel methods poised for future (pre)clinical testing.

Identification of V6.51L as a selectivity hotspot in stereoselective A2B adenosine receptor antagonist recognition.

The four adenosine receptors (ARs) A1AR, A2AAR, A2BAR, and A3AR are G protein-coupled receptors (GPCRs) for which an exceptional amount of experimental and structural data is available. Still, limited success has been achieved in getting new chemical modulators on the market. As such, there is a clear interest in the design of novel selective chemical entities for this family of receptors. In this work, we investigate the selective recognition of ISAM-140, a recently reported A2BAR reference antagonist. A combination of semipreparative chiral HPLC, circular dichroism and X-ray crystallography was used to separate and unequivocally assign the configuration of each enantiomer. Subsequently affinity evaluation for both A2A and A2B receptors demonstrate the stereospecific and selective recognition of (S)-ISAM140 to the A2BAR. The molecular modeling suggested that the structural determinants of this selectivity profile would be residue V2506.51 in A2BAR, which is a leucine in all other ARs including the closely related A2AAR. This was herein confirmed by radioligand binding assays and rigorous free energy perturbation (FEP) calculations performed on the L249V6.51 mutant A2AAR receptor. Taken together, this study provides further insights in the binding mode of these A2BAR antagonists, paving the way for future ligand optimization.