Drugs targeting the four adenosine receptor (AR) subtypes can provide “soft" treatment of various significant diseases. Even for the two experimentally resolved AR subtypes the description of the orthosteric binding area and structure-activity relationships of ligands remains a demanding task due to the high similar amino acids sequence but also the broadness and flexibility of the ARs binding area. The identification of new pharmacophoric moieties and nanomolar leads and the exploration of their binding area with mutagenesis and state-of-the-art computational methods useful also for drug design purposes remains a challenging aim for all ARs. Here, we identified several low nanomolar ligands and potent competitive antagonists against A1R / A3R, containing the novel pyrazolo[3,4-c]pyridine pharmacophore for ARs, from a screen of an in-house library of only 52 compounds, originally designed for anti-proliferative activity. We identified L2-L10, A15, A17 with 3-aryl, 7-anilino and a electronegative group at 5-position as low micromolar to low nanomolar A1R / A3R antagonists. A17 has for A1R Kd = 5.62 nM and a residence time (RT) 41.33 min and for A3R Kd = 13.5 nM, RT = 47.23 min. The kinetic data showed that compared to the not potent or mediocre congeners the active compounds have similar association, for example at A1R Kon = 13.97 x106 M-1 (A17) vs Kon = 3.36 x106 M-1 (A26) but much lower dissociation rate Koff = 0.024 min-1 (A17) vs 0.134 min-1 (A26). Using molecular dynamics (MD) simulations and mutagenesis experiments we investigated the binding site of A17 showing that it can interact with an array of residues in transmembrane helix 5 (TM5), TM6, TM7 of A1R or A3R including residues E5.30, E5.28, T7.35 in A1R instead of Q5.28, V5.30 , L7.35 in A3R. A striking observation for drug design purposes is that for L2506.51A the binding affinity of A17 significantly increased at A1R. A17 provides a lead representative of a promising series and by means of the Thermodynamics Integration coupled with MD simulations (TI/MD) method, first applied here on whole GPCR- membrane system and showing a very good agreement between calculated and experimental relative binding free energies for A1R and A3R (spearman rank correlation p = 0.82 and 0.84, respectively), and kinetic experiments can lead to ligands with improved profile against ARs.