Computational modelling reveals novel insights into GnRH receptor activation and binding dynamics

Gonadotrophin-releasing hormone (GnRH) regulates the mammalian reproductive system by binding to its receptor (GnRH1R) and is a target for treating reproductive hormone-dependent disorders and cancers. While the inactive structure of GnRH1R is known, the active conformation and GnRH binding mode that lead to receptor activation are not fully understood. The mechanism of GnRH-induced receptor activation remains poorly understood due to the absence of experimental structures of the active GnRH1R-GnRH complex. To address this gap, we employed computational docking simulations using Rosetta, coupled with a custom Python-based elimination protocol, to identify near-native binding poses. This approach yielded two top-ranked candidates, ROS-1 and ROS-2. Molecular dynamics simulations revealed that ROS-1 induced GnRH1R activation within 1.0 μs, characterised by a ≈ 4 Å outward shift of the cytoplasmic end of TM6. Key interactions included [Formula: see text] stacking between GnRH and GnRH1R (notably Y5 with Y283[Formula: see text], Y290[Formula: see text], and F309[Formula: see text]) and hydrogen bonds with L286[Formula: see text]. Intramolecular [Formula: see text] interactions within GnRH (Y5 and W3) also played a significant role. Two main communication pathways initiated by R8 of GnRH were identified. R8 formed cation-π interactions with W280[Formula: see text] and communicated with N87[Formula: see text] and the DPxxY motif via water-mediated hydrogen bonds. Additional interactions involved M125[Formula: see text] and the PAF and DRS motifs, which are critical for receptor activation. Key differences in [Formula: see text] interactions at they cytosolic end of TM7 between active and inactive states were identified due to the reorganisation of the DPxxY motif. Finally, GnRH1R communication with lipids through hydrogen bonds involving R240[Formula: see text], R75[Formula: see text], and S140[Formula: see text] was observed. This study provides insights into the active conformation and binding dynamics of the GnRH-GnRH1R complex, advancing our current understanding by providing a coherent picture that consolidates previous interpretations, thereby paving the way to better therapeutic applications.