Unraveling the Interplay of Extracellular Domain Conformational Changes and Parathyroid Hormone Type 1 Receptor Activation in Class B1 G Protein-Coupled Receptors: Integrating Enhanced Sampling Molecular Dynamics Simulations and Markov State Models

Parathyroid hormone (PTH) type 1 receptor (PTH1R), as a typical class B1 G protein-coupled receptor (GPCR), is responsible for regulating bone turnover and maintaining calcium homeostasis, and its dysregulation has been implicated in the development of several diseases. The extracellular domain (ECD) of PTH1R is crucial for the recognition and binding of ligands, and the receptor may exhibit an autoinhibited state with the closure of the ECD in the absence of ligands. However, the correlation between ECD conformations and PTH1R activation remains unclear. Thus, this study combines enhanced sampling molecular dynamics (MD) simulations and Markov state models (MSMs) to reveal the possible relevance between the ECD conformations and the activation of PTH1R. First, 22 intermediate structures are generated from the autoinhibited state to the active state and conducted for 10 independent 200 ns simulations each. Then, the MSM is constructed based on the cumulative 44 μs simulations with six identified microstates. Finally, the potential interplay between ECD conformational changes and PTH1R activation as well as cryptic allosteric pockets in the intermediate states during receptor activation is revealed. Overall, our findings reveal that the activation of PTH1R has a specific correlation with ECD conformational changes and provide essential insights for GPCR biology and developing novel allosteric modulators targeting cryptic sites.

Understanding the Allosteric Modulation of PTH1R by a Negative Allosteric Modulator

The parathyroid hormone type 1 receptor (PTH1R) acts as a canonical class B G protein-coupled receptor, regulating crucial functions including calcium homeostasis and bone formation. The identification and development of PTH1R non-peptide allosteric modulators have obtained widespread attention. It has been found that a negative allosteric modulator (NAM) could inhibit the activation of PTH1R, but the implied mechanism remains unclear. Herein, extensive molecular dynamics simulations together with multiple analytical approaches are utilized to unravel the mechanism of PTH1R allosteric inhibition. The results suggest that the binding of NAM destabilizes the structure of the PTH1R-PTH-spep/qpep (the C terminus of Gs/Gq proteins) complexes. Moreover, the presence of NAM weakens the binding of PTH/peps (spep and qpep) and PTH1R. The intra- and inter-molecular couplings are also weakened in PTH1R upon NAM binding. Interestingly, compared with our previous study of the positive allosteric effects induced by extracellular Ca²⁺, the enhanced correlation between the PTH and G-protein binding sites is significantly reduced by the replacement of this negative allosteric regulator. Our findings might contribute to the development of new therapeutic agents for diseases caused by the abnormal activation of PTH1R.

Endogenous ligand recognition and structural transition of a human PTH receptor

Endogenous parathyroid hormone (PTH) and PTH-related peptide (PTHrP) bind to the parathyroid hormone receptor 1 (PTH1R) and activate the stimulatory G-protein (Gs) signaling pathway. Intriguingly, the two ligands have distinct signaling and physiological properties: PTH evokes prolonged Gs activation, whereas PTHrP evokes transient Gs activation with reduced bone-resorption effects. The distinct molecular actions are ascribed to the differences in ligand recognition and dissociation kinetics. Here, we report cryoelectron microscopic structures of six forms of the human PTH1R-Gs complex in the presence of PTH or PTHrP at resolutions of 2.8 -4.1 Å. A comparison of the PTH-bound and PTHrP-bound structures reveals distinct ligand-receptor interactions underlying the ligand affinity and selectivity. Furthermore, five distinct PTH-bound structures, combined with computational analyses, provide insights into the unique and complex process of ligand dissociation from the receptor and shed light on the distinct durations of signaling induced by PTH and PTHrP.

Critical Extracellular Ca2+ Dependence of the Binding between PTH1R and a G-Protein Peptide Revealed by MD Simulations

The parathyroid hormone type 1 receptor (PTH1R), a canonical class B GPCR, is regulated by a positive allosteric modulator, extracellular Ca²⁺. Calcium ions prolong the residence time of PTH on the PTH1R, leading to increased receptor activation and duration of cAMP signaling. But the essential mechanism of the allosteric behavior of PTH1R is not fully understood. Here, extensive molecular dynamics (MD) simulations are performed for the PTH1R-G-protein combinations with and without Ca²⁺ to describe how calcium ions allosterically engage receptor-G-protein coupling. We find that the binding of Ca²⁺ stabilizes the conformation of the PTH1R-PTH-spep (the α5 helix of Gs protein) complex, especially the extracellular loop 1 (ECL1). Moreover, the MM-GBSA result indicates that Ca²⁺ allosterically promotes the interaction between PTH1R and spep, consistent with the observation of steered molecular dynamics (SMD) simulations. We further illuminate the possible allosteric signaling pathway from the stable Ca²⁺-coupling site to the intracellular G-protein binding site. These results unveil structural determinants for Ca²⁺ allosterism in the PTH1R-PTH-spep complex and give insights into pluridimensional GPCR signaling regulated by calcium ions.

Parathyroid hormone-like hormone plays a dual role in neuroblastoma depending on PTH1R expression

We have previously reported the expression of parathyroid hormone‐like hormone (PTHLH) in well‐differentiated, Schwannian stroma‐rich neuroblastic tumors. The aim of this study was to functionally assess the role of PTHLH and its receptor, PTH1R, in neuroblastoma. Stable knockdown of PTHLH and PTH1R was conducted in neuroblastoma cell lines to investigate the succeeding phenotype induced both in vitro and in vivo. Downregulation of PTHLH reduced MYCN expression and subsequently induced cell cycle arrest, senescence, and migration and invasion impairment in a MYCN‐amplified, TP53‐mutated neuroblastoma cell line. These phenotypes were associated with reduced tumorigenicity in a murine model. We also show that PTHLH expression is not under the control of the calcium‐sensing receptor in neuroblastoma. Conversely, its production is stimulated by epidermal growth factor receptor (EGFR). Accordingly, irreversible EGFR inhibition with canertinib abolished PTHLH expression. The oncogenic role of PTHLH appeared to be a consequence of its intracrine function, as downregulation of its receptor, PTH1R, increased anchorage‐independent growth and induced a more undifferentiated, invasive phenotype. Respectively, high PTH1R mRNA expression was found in MYCN nonamplified primary tumors and also significantly associated with other prognostic factors of good outcome. This study provides the first evidence of the dual role of PTHLH in the behavior of neuroblastomas. Moreover, the identification of EGFR as a transcriptional regulator of PTHLH in neuroblastoma provides a novel therapeutic opportunity to promote a less aggressive tumor phenotype through irreversible inhibition of EGFR tyrosine kinase activity.