Studying the collective motions of the adenosine A2A receptor as a result of ligand binding using principal component analysis

Adenosine receptors (ARs) belong to family A of GPCRs that are involved in many diseases, including cerebral and cardiac ischemic diseases, immune and inflammatory disorders, etc. Thus, they represent important therapeutic targets to treat these conditions. Computational techniques such as molecular dynamics (MD) simulations permit researchers to obtain structural information about these proteins, and principal component analysis (PCA) allows for the identification of collective motions. There are available structures for the active form (3QAK) and the inactive form (3EML) of A2AR which permit us to gain insight about their activation/inactivation mechanism. In this work, we have proposed an inverse strategy using MD simulations where the active form was coupled to the antagonist caffeine and the inactive form was coupled to adenosine agonist. Moreover, we have included four reported thermostabilizing mutations in the inactive form to study A2AR structural differences under different conditions. Some observations stand out from the PCA studies. For instance, the apo structures showed remarkable similarities, and the principal components (PCs) were rearranged in a ligand-dependent manner. Additionally, the active conformation was less stable compared to the inactive one. Some PCs inverted their direction in the presence of a ligand, and comparison of the PCs between 3EML and 3EML_ADN showed that adenosine induced major changes in the structure of A2AR. Rearrangement of PCs precedes and drives conformational changes that occur after ligand binding. Knowledge about these conformational changes provides important insights about the activity of A2AR.