Abstract
We present a robust protocol based on iterations of free energy perturbation (FEP) calculations, chemical synthesis, biophysical mapping and X-ray crystallography to reveal the binding mode of an antagonist series to the A(2A) adenosine receptor (AR). Eight A(2A) AR binding site mutations from biophysical mapping experiments were initially analyzed with sidechain FEP simulations, performed on alternate binding modes. The results distinctively supported one binding mode, which was subsequently used to design new chromone derivatives. Their affinities for the A(2A) AR were experimentally determined and investigated through a cycle of ligand-FEP calculations, validating the binding orientation of the different chemical substituents proposed. Subsequent X-ray crystallography of the A(2A) AR with a low and a high affinity chromone derivative confirmed the predicted binding orientation. The new molecules and structures here reported were driven by free energy calculations, and provide new insights on antagonist binding to the A(2A) AR, an emerging target in immuno-oncology.