Selective and potent adenosine A3 receptor antagonists by methoxyaryl substitution on the N-(2,6-diarylpyrimidin-4-yl)acetamide scaffold

Computational Workflow for Refining AlphaFold Models in Drug Design Using Kinetic and Thermodynamic Binding Calculations: A Case Study for the Unresolved Inactive Human Adenosine A3 Receptor

A structure-based drug design pipeline that considers both thermodynamic and kinetic binding data of ligands against a receptor will enable the computational design of improved drug molecules. For unresolved GPCR-ligand complexes, a workflow that can apply both thermodynamic and kinetic binding data in combination with alpha-fold (AF)-derived or other homology models and experimentally resolved binding modes of relevant ligands in GPCR-homologs needs to be tested. Here, as test case, we studied a congeneric set of ligands that bind to a structurally unresolved G protein-coupled receptor (GPCR), the inactive human adenosine A3 receptor (hA3R). We tested three available homology models from which two have been generated from experimental structures of hA1R or hA2AR and one model was a multistate alphafold 2 (AF2)-derived model. We applied alchemical calculations with thermodynamic integration coupled with molecular dynamics (TI/MD) simulations to calculate the experimental relative binding free energies and residence time (τ)-random accelerated MD (τ-RAMD) simulations to calculate the relative residence times (RTs) for antagonists. While the TI/MD calculations produced, for the three homology models, good Pearson correlation coefficients, correspondingly, r = 0.74, 0.62, and 0.67 and mean unsigned error (mue) values of 0.94, 1.31, and 0.81 kcal mol-1, the τ-RAMD method showed r = 0.92 and 0.52 for the first two models but failed to produce accurate results for the multistate AF2-derived model. With subsequent optimization of the AF2-derived model by reorientation of the side chain of R1735.34 located in the extracellular loop 2 (EL2) that blocked ligand's unbinding, the computational model showed r = 0.84 for kinetic data and improved performance for thermodynamic data (r = 0.81, mue = 0.56 kcal mol-1). Overall, after refining the multistate AF2 model with physics-based tools, we were able to show a strong correlation between predicted and experimental ligand relative residence times and affinities, achieving a level of accuracy comparable to an experimental structure. The computational workflow used can be applied to other receptors, helping to rank candidate drugs in a congeneric series and enabling the prioritization of leads with stronger binding affinities and longer residence times.

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

A library of potent and highly A₃AR selective pyrimidine-based compounds was designed to explore non-orthosteric interactions within this receptor. Starting from a prototypical orthosteric A₃AR antagonist (ISVY130), the structure-based design explored functionalized residues at the exocyclic amide L1 region and aimed to provide additional interactions outside the A₃AR orthosteric site. The novel ligands were assembled through an efficient and succinct synthetic approach, resulting in compounds that retain the A₃AR 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 A₃AR ligands, including covalent drugs and molecular probes for diagnostic and delivery purposes.

Adenosine receptor antagonists: Recent advances and therapeutic perspective

Adenosine is an endogenous purine-based nucleoside expressed nearly in all body tissues. It regulates various body functions by activating four G-protein coupled receptors, A₁, A_2A, A_2B, and A₃. These receptors are widely acknowledged as drug targets for treating different neurological, metabolic, and inflammatory diseases. Although numerous adenosine receptor inhibitors have been developed worldwide, achieving target selectivity is still a big hurdle in drug development. However, the identification of specific radioligands-based affinity assay, fluorescent ligands, and MS-based ligand assay have contributed to the development of selective and potent adenosine ligands. In recent years various small heterocyclic-based molecules have shown some promising results. Istradefylline has been approved for treating Parkinson's in Japan, while preladenant, tozadenant, CVT-6883, MRS-1523, and many more are under different phases of clinical development. The present review is focused on the quest to develop potent and selective adenosine inhibitors from 2013 to early 2021 by various research groups. The review also highlights their biological activity, selectivity, structure-activity relationship, molecular docking, and mechanistic studies. A special emphsesis on drug designing strategies has been also given the manuscript. The comprehensive compilation of research work carried out in the field will provide inevitable scope for designing and developing novel adenosine inhibitors with improved selectivity and efficacy.

Pharmacological characterisation of novel adenosine A3 receptor antagonists

The adenosine A3 receptor (A3R) belongs to a family of four adenosine receptor (AR) subtypes which all play distinct roles throughout the body. A3R antagonists have been described as potential treatments for numerous diseases including asthma. Given the similarity between (adenosine receptors) orthosteric binding sites, obtaining highly selective antagonists is a challenging but critical task. Here we screen 39 potential A3R, antagonists using agonist-induced inhibition of cAMP. Positive hits were assessed for AR subtype selectivity through cAMP accumulation assays. The antagonist affinity was determined using Schild analysis (pA2 values) and fluorescent ligand binding. Structure-activity relationship investigations revealed that loss of the 3-(dichlorophenyl)-isoxazolyl moiety or the aromatic nitrogen heterocycle with nitrogen at α-position to the carbon of carboximidamide group significantly attenuated K18 antagonistic potency. Mutagenic studies supported by molecular dynamic simulations combined with Molecular Mechanics-Poisson Boltzmann Surface Area calculations identified the residues important for binding in the A3R orthosteric site. We demonstrate that K18, which contains a 3-(dichlorophenyl)-isoxazole group connected through carbonyloxycarboximidamide fragment with a 1,3-thiazole ring, is a specific A3R (< 1 µM) competitive antagonist. Finally, we introduce a model that enables estimates of the equilibrium binding affinity for rapidly disassociating compounds from real-time fluorescent ligand-binding studies. These results demonstrate the pharmacological characterisation of a selective competitive A3R antagonist and the description of its orthosteric binding mode. Our findings may provide new insights for drug discovery.

Structure-Based Design of Potent and Selective Ligands at the Four Adenosine Receptors

The four receptors that signal for adenosine, A₁, A_2A, A_2B and A₃ ARs, belong to the superfamily of G protein-coupled receptors (GPCRs). They mediate a number of (patho)physiological functions and have attracted the interest of the biopharmaceutical sector for decades as potential drug targets. The many crystal structures of the A_2A, and lately the A₁ ARs, allow for the use of advanced computational, structure-based ligand design methodologies. Over the last decade, we have assessed the efficient synthesis of novel ligands specifically addressed to each of the four ARs. We herein review and update the results of this program with particular focus on molecular dynamics (MD) and free energy perturbation (FEP) protocols. The first in silico mutagenesis on the A₁AR here reported allows understanding the specificity and high affinity of the xanthine-antagonist 8-Cyclopentyl-1,3-dipropylxanthine (DPCPX). On the A_2AAR, we demonstrate how FEP simulations can distinguish the conformational selectivity of a recent series of partial agonists. These novel results are complemented with the revision of the first series of enantiospecific antagonists on the A_2BAR, and the use of FEP as a tool for bioisosteric design on the A₃AR.

Bis-pyrimidine acetamides: design, synthesis and biological evaluation

BACKGROUND

In the past few years, increased resistance of microorganisms towards antimicrobial agents become a serious health problem, so there is a need for the discovery of new antimicrobial agents. On the other hand, bis-pyrimidines possess various types of biological activity. In view of this, in the present study we have designed and synthesized a new series of bis-pyrimidine acetamides by Claisen-Schmidt condensation and screened for its in vitro antimicrobial and anticancer activities.

RESULTS

The synthesized bis-pyrimidine acetamide derivatives were confirmed by IR, 1H/13C-NMR, Mass spectral studies as well C, H, N analyses. The synthesized compounds were evaluated for their in vitro antimicrobial potential against Gram positive (Staphylococcus aureus and Bacillus subtilis); Gram negative (Escherichia coli, Pseudomonas aeruginosa and Salmonella enterica) bacterial and fungal (Candida albicans and Aspergillus niger) strains by tube dilution technique and the minimum inhibitory concentration (MIC) recorded in µmol/mL was comparable to reference drugs, cefadroxil (antibacterial) and fluconazole (antifungal). The in vitro anticancer activity (IC50 value) determined against human colorectal carcinoma (HCT116) cancer cell line by Sulforhodamine B (SRB) technique and 5-fluorouracil used as reference drug.

CONCLUSIONS

The biological study demonstrated that compounds 3, 13, 16, 17 and 18 were found to be most active antimicrobial agents with best MIC values than the cefadroxil (antibacterial) and fluconazole (antifungal) and compounds 12, 16 and 18 found to have better anticancer activity against human colorectal carcinoma (HCT116) cancer cell line with best IC50 value than the 5-fluorouracil (anticancer). Graphical abstract SAR of bis-pyrimidine acetamides.

Expanding the horizons of G protein-coupled receptor structure-based ligand discovery and optimization using homology models

We show the key role of structural homology models in GPCR structure-based lead discovery and optimization, highlighting methodological aspects, recent progress and future directions.

Structure-based drug design for G protein-coupled receptors

Our understanding of the structural biology of G protein-coupled receptors has undergone a transformation over the past 5 years. New protein-ligand complexes are described almost monthly in high profile journals. Appreciation of how small molecules and natural ligands bind to their receptors has the potential to impact enormously how medicinal chemists approach this major class of receptor targets. An outline of the key topics in this field and some recent examples of structure- and fragment-based drug design are described. A table is presented with example views of each G protein-coupled receptor for which there is a published X-ray structure, including interactions with small molecule antagonists, partial and full agonists. The possible implications of these new data for drug design are discussed.