Understanding the conformation transition in the activation pathway of β2 adrenergic receptor via a targeted molecular dynamics simulation

G protein coupled receptors (GPCRs) play a crucial role in regulating signal recognition and transduction through their activation. The conformation transition in the activation pathway is of particular importance for their function. However, it has been poorly elucidated due to experimental difficulties in determining the conformations and the time limitation of conventional molecular dynamics (CMD) simulation. Thus, in this work, we employed a targeted molecular dynamic (TMD) simulation to study the activation process from an inactive structure to a fully active one for β2 adrenergic receptor (β2AR). As a reference, 110 ns CMD simulations on wild β2AR and its D130N mutant were also carried out. TMD results show that there is at least an intermediate conformation cluster in the activation process, evidenced by the principal component analysis and the structural and dynamic differences of some important motifs. It is noteworthy that the activation of the ligand binding site lags the G-protein binding site, displaying uncoupled correlation. Comparisons between the CMD and TMD results show that the D130N mutation significantly speeds up ICL2 and key ionic lock to enter into the intermediate state, which to some extent facilitates the activation involved in the NPxxY, DRY region and the separation between TM3 and TM6. However, the contribution from the D130N mutation to the activation of the ligand binding site could not be observed within the scale of 110 ns time. These observations could provide novel insights into previous studies for better understanding of the activation mechanism for β2AR.