Abstract
The full activation process of G protein-coupled receptor (GPCR) plays an important role in cellular signal transduction. However, it remains challenging to simulate the whole process in which the GPCR is recognized and activated by a ligand and then couples to the G protein on a reasonable simulation timescale. Here, we developed a molecular dynamics (MD) approach named supervised (Su) Gaussian accelerated MD (GaMD) by incorporating a tabu-like supervision algorithm into a standard GaMD simulation. By using this Su-GaMD method, from the active and inactive structure of adenosine A(1) receptor (A(1)R), we successfully revealed the full activation mechanism of A(1)R, including adenosine (Ado)-A(1)R recognition, preactivation of A(1)R, and A(1)R-G protein recognition, in hundreds of nanoseconds of simulations. The binding of Ado to the extracellular side of A(1)R initiates conformational changes and the preactivation of A(1)R. In turn, the binding of G(i2) to the intracellular side of A(1)R causes a decrease in the volume of the extracellular orthosteric site and stabilizes the binding of Ado to A(1)R. Su-GaMD could be a useful tool to reconstruct or even predict ligand-protein and protein-protein recognition pathways on a short timescale. The intermediate states revealed in this study could provide more detailed complementary structural characterizations to facilitate the drug design of A(1)R in the future.