Abstract
The parathyroid hormone type 1 receptor (PTH1R) plays a crucial role in modulating various physiological functions and is considered an effective therapeutic target for osteoporosis. However, a lack of detailed molecular and energetic information about PTH1R limits our comprehensive understanding of its activation process. In this study, we performed computational simulations to explore key events in the activation process, such as conformational changes in PTH1R, Gs protein coupling, and the release of guanosine diphosphate (GDP). Our analysis identified kinetic information, including the rate-determining step, transition state, and energy barriers. Free-energy and structural analyses revealed that GDP could be released from the Gs protein when the binding cavity is partially open. Additionally, we predicted important residues, including potential pathogenic mutations, and verified their significance through site-directed mutations. These findings enhance our understanding of class B GPCR activation mechanisms. Furthermore, the methodology employed in this study can be applied to other biophysical systems.