Molecular dynamics simulations reveal initial structural and dynamic features for the A ₂AR as a result of ligand binding

G-protein-coupled receptors (GPCRs) are membrane proteins that have a wide variety of physiological roles. Adenosine receptors belong to the GPCR family. Adenosine receptors are implicated in many physiological disorders, such as Parkinson's disease, Huntington's disease, inflammatory and immune's disease and many others. Interestingly, crystal structures of the active and inactive conformations of the A₂-subtype adenosine receptor (A₂AR) have been solved. These two structures could be used to get insights about the conformational changes that occur during the process of activation/inactivation processes of this receptor. Therefore, two ligand-free simulations of the native active (PDB code: 3QAK) and inactive (PDB code: 3EML) conformations of the A₂AR and two halo-simulations were carried out to observe the initial conformational changes induced by coupling adenosine to the inactive conformation and caffeine to the active conformation. Furthermore, we constructed an A₂AR model that contained four thermostabilising mutations, L48A, T65A, Q89A and A54L, which had previously been determined to stabilise the bound conformation of the agonist, and we ran molecular dynamics simulations of this mutant to investigate how these point mutations might affect the inactive conformation of this receptor. This study provides insights about the initial structural and dynamic features that occur as a result of the binding of caffeine and adenosine in the active and inactive A₂AR structures, respectively, as well as the introduction of some mutations on the inactive structure of the A₂AR. Moreover, we provide useful and detailed information regarding structural features such as toggle switch and ionic lock during the activation/inactivation processes of this receptor.