Dual A1/A3 Adenosine Receptor Antagonists: Binding Kinetics and Structure-Activity Relationship Studies Using Mutagenesis and Alchemical Binding Free Energy Calculations.

Drugs targeting adenosine receptors (AR) can provide treatment for diseases. We report the identification of 7-(phenylamino)-pyrazolo[3,4-c]pyridines L2-L10, A15, and A17 as low-micromolar to low-nanomolar A1R/A3R dual antagonists, with 3-phenyl-5-cyano-7-(trimethoxyphenylamino)-pyrazolo[3,4-c]pyridine (A17) displaying the highest affinity at both receptors with a long residence time of binding, as determined using a NanoBRET-based assay. Two binding orientations of A17 produce stable complexes inside the orthosteric binding area of A1R in molecular dynamics (MD) simulations, and we selected the most plausible orientation based on the agreement with alanine mutagenesis supported by affinity experiments. Interestingly, for drug design purposes, the mutation of L2506.51 to alanine increased the binding affinity of A17 at A1R. We explored the structure-activity relationships against A1R using alchemical binding free energy calculations with the thermodynamic integration coupled with the MD simulation (TI/MD) method, applied on the whole G-protein-coupled receptor-membrane system, which showed a good agreement (r = 0.73) between calculated and experimental relative binding free energies.

Discovery of a High Affinity Adenosine A1/A3 Receptor Antagonist with a Novel 7-Amino-pyrazolo[3,4-d]pyridazine Scaffold.

Here we describe the design and synthesis of pyrazolo[3,4-d]pyridazines as adenosine receptor (AR) ligands. We demonstrate that the introduction of a 3-phenyl group, together with a 7-benzylamino and 1-methyl group at the pyrazolopyridazine scaffold, generated the antagonist compound 10b, which displayed 21 nM affinity and a residence time of ∼60 min, for the human A1R, 55 nM affinity and a residence time of ∼73 min, for the human A3R and 1.7 µΜ affinity for the human A2BR while not being toxic. Strikingly, the 2-methyl analog of 10b, 15b, had no significant affinity. Docking calculations and molecular dynamics simulations of the ligands inside the orthosteric binding area suggested that the 2-methyl group in 15b hinders the formation of hydrogen bonding interactions with N6.55 which are considered critical for the stabilization inside the orthosteric binding cavity. We, therefore, demonstrate that 10a is a novel scaffold for the development of high affinity AR ligands. From the mutagenesis experiments the biggest effect was observed for the Y2717.46A mutation which caused an ∼10-fold reduction in the binding affinity of 10b.