All roads lead to Rome: Emergence of Mechanistic Divergence Shapes Activation of Class B1 G-Protein Coupled Receptors

Class B1 GPCRs are crucial to maintaining homeostasis along a multitude of vital biochemical pathways. Understanding the activation mechanism of these proteins at both a family and clade-specific level is particularly relevant for designing multi-target agonists, as exemplified by recently designed dual-agonists for GLP-1R and GIPR, for treating obesity. Here, we use 6 milliseconds of unbiased all-atom MD simulations of GCGR, GLP1R, PAC1R, SCTR, PTH1R and CALCR from the four different clades of Class B1 GPCRs to establish the universal mechanism of their activation. We show that the activation of Class B1 GPCRs involves a clade-independent intermediate state characterized by the outward movement of helix 6. We use a combination of Markov state models and transition path theory to show that the activation of these proteins occurs at a millisecond timescale. We identify characteristic molecular locks that are conserved at a clade-level, showcasing the uniqueness among the activation mechanisms of these proteins. We show that these proteins show similar inactive and active states, but show unique activation mechanisms at a residue level. These sites can be targeted directly or allosterically to design therapeutics targeting a specific clade of proteins. Thus, this study provides an integrated atomistic view of the activation for Class B1 GPCRs from a mechanistic, thermodynamic and kinetic perspective.