Exploring Biginelli-based scaffolds as A2B adenosine receptor antagonists: Unveiling novel structure-activity relationship trends, lead compounds, and potent colorectal anticancer agents.

Antagonists of the A2B adenosine receptor have recently emerged as targeted anticancer agents and immune checkpoint inhibitors within the realm of cancer immunotherapy. This study presents a comprehensive evaluation of novel Biginelli-assembled pyrimidine chemotypes, including mono-, bi-, and tricyclic derivatives, as A2BAR antagonists. We conducted a comprehensive examination of the adenosinergic profile (both binding and functional) of a large compound library consisting of 168 compounds. This approach unveiled original lead compounds and enabled the identification of novel structure-activity relationship (SAR) trends, which were supported by extensive computational studies, including quantum mechanical calculations and free energy perturbation (FEP) analysis. In total, 25 molecules showed attractive affinity (Ki < 100 nM) and outstanding selectivity for A2BAR. From these, five molecules corresponding to the new benzothiazole scaffold were below the Ki < 10 nM threshold, in addition to a novel dual A2A/A2B antagonist. The most potent compounds, and the dual antagonist, showed enantiospecific recognition in the A2BAR. Two A2BAR selective antagonists and the dual A2AAR/A2BAR antagonist reported in this study were assessed for their impact on colorectal cancer cell lines. The results revealed a significant and dose-dependent reduction in cell proliferation. Notably, the A2BAR antagonists exhibited remarkable specificity, as they did not impede the proliferation of non-tumoral cell lines. These findings support the efficacy and potential that A2BAR antagonists as valuable candidates for cancer therapy, but also that they can effectively complement strategies involving A2AAR antagonism in the context of immune checkpoint inhibition.

Visualizing the impact of disease-associated mutations on G protein-nucleotide interactions.

Activation of G proteins stimulates ubiquitous intracellular signaling cascades essential for life processes. Under normal physiological conditions, nucleotide exchange is initiated upon the formation of complexes between a G protein and G protein-coupled receptor (GPCR), which facilitates exchange of bound GDP for GTP, subsequently dissociating the trimeric G protein into its Gα and Gßγ subunits. However, single point mutations in Gα circumvent nucleotide exchange regulated by GPCR-G protein interactions, leading to either loss-of-function or constitutive gain-of-function. Mutations in several Gα subtypes are closely linked to the development of multiple diseases, including several intractable cancers. We leveraged an integrative spectroscopic and computational approach to investigate the mechanisms by which seven of the most frequently observed clinically-relevant mutations in the α subunit of the stimulatory G protein result in functional changes. Variable temperature circular dichroism (CD) spectroscopy showed a bimodal distribution of thermal melting temperatures across all GαS variants. Modeling from molecular dynamics (MD) simulations established a correlation between observed thermal melting temperatures and structural changes caused by the mutations. Concurrently, saturation-transfer difference NMR (STD-NMR) highlighted variations in the interactions of GαS variants with bound nucleotides. MD simulations indicated that changes in local interactions within the nucleotide-binding pocket did not consistently align with global structural changes. This collective evidence suggests a multifaceted energy landscape, wherein each mutation may introduce distinct perturbations to the nucleotide-binding site and protein-protein interaction sites. Consequently, it underscores the importance of tailoring therapeutic strategies to address the unique challenges posed by individual mutations.

N-(Heteroaryl)thiophene sulfonamides as angiotensin AT2 receptor ligands.

Two series of N-(heteroaryl)thiophene sulfonamides, encompassing either a methylene imidazole group or a tert-butylimidazolylacetyl group in the meta position of the benzene ring, have been synthesized. An AT2R selective ligand with a Ki of 42 nM was identified in the first series and in the second series, six AT2R selective ligands with significantly improved binding affinities and Ki values of <5 nM were discovered. The binding modes to AT2R were explored by docking calculations combined with molecular dynamics simulations. Although some of the high affinity ligands exhibited fair stability in human liver microsomes, comparable to that observed with C21 undergoing clinical trials, most ligands displayed a very low metabolic stability with t½ of less than 10 min in human liver microsomes. The most promising ligand, with an AT2R Ki value of 4.9 nM and with intermediate stability in human hepatocytes (t½ = 77 min) caused a concentration-dependent vasorelaxation of pre-contracted mouse aorta.

Exploring the Effect of Halogenation in a Series of Potent and Selective A2B Adenosine Receptor Antagonists.

The modulation of the A2B adenosine receptor is a promising strategy in cancer (immuno) therapy, with A2BAR antagonists emerging as immune checkpoint inhibitors. Herein, we report a systematic assessment of the impact of (di- and mono-)halogenation at positions 7 and/or 8 on both A2BAR affinity and pharmacokinetic properties of a collection of A2BAR antagonists and its study with structure-based free energy perturbation simulations. Monohalogenation at position 8 produced potent A2BAR ligands irrespective of the nature of the halogen. In contrast, halogenation at position 7 and dihalogenation produced a halogen-size-dependent decay in affinity. Eight novel A2BAR ligands exhibited remarkable affinity (Ki < 10 nM), exquisite subtype selectivity, and enantioselective recognition, with some eutomers eliciting sub-nanomolar affinity. The pharmacokinetic profile of representative derivatives showed enhanced solubility and microsomal stability. Finally, two compounds showed the capacity of reversing the antiproliferative effect of adenosine in activated primary human peripheral blood mononuclear cells.

A2B adenosine receptor antagonists rescue lymphocyte activity in adenosine-producing patient-derived cancer models.

BACKGROUND: Adenosine is a metabolite that suppresses antitumor immune response of T and NK cells via extracellular binding to the two subtypes of adenosine-2 receptors, A2ARs. While blockade of the A2AARs subtype effectively rescues lymphocyte activity, with four A2AAR antagonists currently in anticancer clinical trials, less is known for the therapeutic potential of the other A2BAR blockade within cancer immunotherapy. Recent studies suggest the formation of A2AAR/A2BAR dimers in tissues that coexpress the two receptor subtypes, where the A2BAR plays a dominant role, suggesting it as a promising target for cancer immunotherapy. METHODS: We report the synthesis and functional evaluation of five potent A2BAR antagonists and a dual A2AAR/A2BAR antagonist. The compounds were designed using previous pharmacological data assisted by modeling studies. Synthesis was developed using multicomponent approaches. Flow cytometry was used to evaluate the phenotype of T and NK cells on A2BAR antagonist treatment. Functional activity of T and NK cells was tested in patient-derived tumor spheroid models. RESULTS: We provide data for six novel small molecules: five A2BAR selective antagonists and a dual A2AAR/A2BAR antagonist. The growth of patient-derived breast cancer spheroids is prevented when treated with A2BAR antagonists. To elucidate if this depends on increased lymphocyte activity, immune cells proliferation, and cytokine production, lymphocyte infiltration was evaluated and compared with the potent A2AAR antagonist AZD-4635. We find that A2BAR antagonists rescue T and NK cell proliferation, IFNγ and perforin production, and increase tumor infiltrating lymphocytes infiltration into tumor spheroids without altering the expression of adhesion molecules. CONCLUSIONS: Our results demonstrate that A2BAR is a promising target in immunotherapy, identifying ISAM-R56A as the most potent candidate for A2BAR blockade. Inhibition of A2BAR signaling restores T cell function and proliferation. Furthermore, A2BAR and dual A2AAR/A2BAR antagonists showed similar or better results than A2AAR antagonist AZD-4635 reinforcing the idea of dominant role of the A2BAR in the regulation of the immune system.

Exploring Non-orthosteric Interactions with a Series of Potent and Selective A3 Antagonists.

A library of potent and highly A3AR selective pyrimidine-based compounds was designed to explore non-orthosteric interactions within this receptor. Starting from a prototypical orthosteric A3AR antagonist (ISVY130), the structure-based design explored functionalized residues at the exocyclic amide L1 region and aimed to provide additional interactions outside the A3AR orthosteric site. The novel ligands were assembled through an efficient and succinct synthetic approach, resulting in compounds that retain the A3AR potent and selective profile while improving the solubility of the original scaffold. The experimentally demonstrated tolerability of the L1 region to structural functionalization was further assessed by molecular dynamics simulations, giving hints of the non-orthosteric interactions explored by these series. The results pave the way to explore newly functionalized A3AR ligands, including covalent drugs and molecular probes for diagnostic and delivery purposes.

Optimization of 2-Amino-4,6-diarylpyrimidine-5-carbonitriles as Potent and Selective A1 Antagonists.

We herein document a large collection of 108 2-amino-4,6-disubstituted-pyrimidine derivatives as potent, structurally simple, and highly selective A1AR ligands. The most attractive ligands were confirmed as antagonists of the canonical cyclic adenosine monophosphate pathway, and some pharmacokinetic parameters were preliminarilly evaluated. The library, built through a reliable and efficient three-component reaction, comprehensively explored the chemical space allowing the identification of the most prominent features of the structure-activity and structure-selectivity relationships around this scaffold. These included the influence on the selectivity profile of the aromatic residues at positions R4 and R6 of the pyrimidine core but most importantly the prominent role to the unprecedented A1AR selectivity profile exerted by the methyl group introduced at the exocyclic amino group. The structure-activity relationship trends on both A1 and A2AARs were conveniently interpreted with rigorous free energy perturbation simulations, which started from the receptor-driven docking model that guided the design of these series.

Slow conformational dynamics of the human A2A adenosine receptor are temporally ordered.

A more complete depiction of protein energy landscapes includes the identification of different function-related conformational states and the determination of the pathways connecting them. We used total internal reflection fluorescence (TIRF) imaging to investigate the conformational dynamics of the human A2A adenosine receptor (A2AAR), a class A G protein-coupled receptor (GPCR), at the single-molecule level. Slow, reversible conformational exchange was observed among three different fluorescence emission states populated for agonist-bound A2AAR. Transitions among these states predominantly occurred in a specific order, and exchange between inactive and active-like conformations proceeded through an intermediate state. Models derived from molecular dynamics simulations with available A2AAR structures rationalized the relative fluorescence emission intensities for the highest and lowest emission states but not the transition state. This suggests that the functionally critical intermediate state required to achieve activation is not currently visualized among available A2AAR structures.

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.

Potent Inhibition of Nicotinamide N-Methyltransferase by Alkene-Linked Bisubstrate Mimics Bearing Electron Deficient Aromatics.

Nicotinamide N-methyltransferase (NNMT) methylates nicotinamide (vitamin B3) to generate 1-methylnicotinamide (MNA). NNMT overexpression has been linked to a variety of diseases, most prominently human cancers, indicating its potential as a therapeutic target. The development of small-molecule NNMT inhibitors has gained interest in recent years, with the most potent inhibitors sharing structural features based on elements of the nicotinamide substrate and the S-adenosyl-l-methionine (SAM) cofactor. We here report the development of new bisubstrate inhibitors that include electron-deficient aromatic groups to mimic the nicotinamide moiety. In addition, a trans-alkene linker was found to be optimal for connecting the substrate and cofactor mimics in these inhibitors. The most potent NNMT inhibitor identified exhibits an IC50 value of 3.7 nM, placing it among the most active NNMT inhibitors reported to date. Complementary analytical techniques, modeling studies, and cell-based assays provide insights into the binding mode, affinity, and selectivity of these inhibitors.

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.

Design, synthesis, HER2 inhibition and anticancer evaluation of new substituted 1,5-dihydro-4,1-benzoxazepines.

A series of 11 new substituted 1,5-dihydro-4,1-benzoxazepine derivatives was synthesised to study the influence of the methyl group in the 1-(benzenesulphonyl) moiety, the replacement of the purine by the benzotriazole bioisosteric analogue, and the introduction of a bulky substituent at position 6 of the purine, on the biological effects. Their inhibition against isolated HER2 was studied and the structure-activity relationships have been confirmed by molecular modelling studies. The most potent compound against isolated HER2 is 9a with an IC50 of 7.31 µM. We have investigated the effects of the target compounds on cell proliferation. The most active compound (7c) against all the tumour cell lines studied (IC50 0.42-0.86 µM) does not produce any modification in the expression of pro-caspase 3, but increases the caspase 1 expression, and promotes pyroptosis.

Potent and Subtype-Selective Dopamine D2 Receptor Biased Partial Agonists Discovered via an Ugi-Based Approach.

Using a previously unexplored, efficient, and versatile multicomponent method, we herein report the rapid generation of novel potent and subtype-selective DRD2 biased partial agonists. This strategy exemplifies the search for diverse and previously unexplored moieties for the secondary/allosteric pharmacophore of the common phenyl-piperazine scaffold. The pharmacological characterization of the new compound series led to the identification of several ligands with excellent DRD2 affinity and subtype selectivity and remarkable functional selectivity for either the cAMP (22a and 24d) or the ß-arrestin (27a and 29c) signaling pathways. These results were further interpreted on the basis of molecular models of these ligands in complex with the recent DRD2 crystal structures, highlighting the critical role of the secondary/allosteric pharmacophore in modulating the functional selectivity profile.

Structural Basis of Inhibition of Human Insulin-Regulated Aminopeptidase (IRAP) by Benzopyran-Based Inhibitors.

Inhibition of the insulin-regulated aminopeptidase (IRAP) improves memory and cognition in animal models. The enzyme has recently been crystallized and several series of inhibitors reported. We herein focused on one series of benzopyran-based inhibitors of IRAP known as the HFI series, with unresolved binding mode to IRAP, and developed a robust computational model to explain the structure-activity relationship (SAR) and potentially guide their further optimization. The binding model here proposed places the benzopyran ring in the catalytic binding site, coordinating the Zn2+ ion through the oxygen in position 3, in contrast to previous hypothesis. The whole series of HFI compounds was then systematically simulated, starting from this binding mode, using molecular dynamics and binding affinity estimated with the linear interaction energy (LIE) method. The agreement with experimental affinities supports the binding mode proposed, which was further challenged by rigorous free energy perturbation (FEP) calculations. Here, we found excellent correlation between experimental and calculated binding affinity differences, both between selected compound pairs and also for recently reported experimental data concerning the site directed mutagenesis of residue Phe544. The computationally derived structure-activity relationship of the HFI series and the understanding of the involvement of Phe544 in the binding of this scaffold provide valuable information for further lead optimization of novel IRAP inhibitors.

Free Energy Calculations for Protein-Ligand Binding Prediction.

Computational prediction of protein-ligand binding involves initial determination of the binding mode and subsequent evaluation of the strength of the protein-ligand interactions, which directly correlates with ligand binding affinities. As a consequence of increasing computer power, rigorous approaches to calculate protein-ligand binding affinities, such as free energy perturbation (FEP) methods, are becoming an essential part of the toolbox of computer-aided drug design. In this chapter, we provide a general overview of these methods and introduce the QFEP modules, which are open-source API workflows based on our molecular dynamics (MD) package Q. The module QligFEP allows estimation of relative binding affinities along ligand series, while QresFEP is a module to estimate binding affinity shifts caused by single-point mutations of the protein. We herein provide guidelines for the use of each of these modules based on data extracted from ligand-design projects. While these modules are stand-alone, the combined use of the two workflows in a drug-design project yields complementary perspectives of the ligand binding problem, providing two sides of the same coin. The selected case studies illustrate how to use QFEP to approach the two key questions associated with ligand binding prediction: identifying the most favorable binding mode from different alternatives and establishing structure-affinity relationships that allow the rational optimization of hit compounds.

3,4-Dihydropyrimidin-2(1H)-ones as Antagonists of the Human A2B Adenosine Receptor: Optimization, Structure-Activity Relationship Studies, and Enantiospecific Recognition.

We present and thoroughly characterize a large collection of 3,4-dihydropyrimidin-2(1H)-ones as A2BAR antagonists, an emerging strategy in cancer (immuno) therapy. Most compounds selectively bind A2BAR, with a number of potent and selective antagonists further confirmed by functional cyclic adenosine monophosphate experiments. The series was analyzed with one of the most exhaustive free energy perturbation studies on a GPCR, obtaining an accurate model of the structure-activity relationship of this chemotype. The stereospecific binding modeled for this scaffold was confirmed by resolving the two most potent ligands [(±)-47, and (±)-38 Ki = 10.20 and 23.6 nM, respectively] into their two enantiomers, isolating the affinity on the corresponding (S)-eutomers (Ki = 6.30 and 11.10 nM, respectively). The assessment of the effect in representative cytochromes (CYP3A4 and CYP2D6) demonstrated insignificant inhibitory activity, while in vitro experiments in three prostate cancer cells demonstrated that this pair of compounds exhibits a pronounced antimetastatic effect.

N-(Methyloxycarbonyl)thiophene sulfonamides as high affinity AT2 receptor ligands.

A series of meta-substituted acetophenone derivatives, encompassing N-(alkyloxycarbonyl)thiophene sulfonamide fragments have been synthesized. Several selective AT2 receptor ligands were identified, among those a tert-butylimidazole derivative (20) with a Ki of 9.3 nM, that demonstrates a high stability in human liver microsomes (t½ = 62 min) and in human hepatocytes (t½ = 194 min). This methyloxycarbonylthiophene sulfonamide is a 20-fold more potent binder to the AT2 receptor and is considerably more stable in human liver microsomes, than a previously reported and broadly studied structurally related AT2R prototype antagonist 3 (C38). Ligand 20 acts as an AT2R agonist and caused an AT2R mediated concentration-dependent vasorelaxation of pre-contracted mouse aorta. Furthermore, in contrast to imidazole derivative C38, the tert-butylimidazole derivative 20 is a poor inhibitor of CYP3A4, CYP2D6 and CYP2C9. It is demonstrated herein that smaller alkyloxycarbonyl groups make the ligands in this series of AT2R selective compounds less prone to degradation and that a high AT2 receptor affinity can be retained after truncation of the alkyloxycarbonyl group. Binding modes of the most potent AT2R ligands were explored by docking calculations combined with molecular dynamics simulations.

Publisher Correction: GPCRmd uncovers the dynamics of the 3D-GPCRome.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

GPCRmd uncovers the dynamics of the 3D-GPCRome.

G-protein-coupled receptors (GPCRs) are involved in numerous physiological processes and are the most frequent targets of approved drugs. The explosion in the number of new three-dimensional (3D) molecular structures of GPCRs (3D-GPCRome) over the last decade has greatly advanced the mechanistic understanding and drug design opportunities for this protein family. Molecular dynamics (MD) simulations have become a widely established technique for exploring the conformational landscape of proteins at an atomic level. However, the analysis and visualization of MD simulations require efficient storage resources and specialized software. Here we present GPCRmd (http://gpcrmd.org/), an online platform that incorporates web-based visualization capabilities as well as a comprehensive and user-friendly analysis toolbox that allows scientists from different disciplines to visualize, analyze and share GPCR MD data. GPCRmd originates from a community-driven effort to create an open, interactive and standardized database of GPCR MD simulations.