Abstract:Glycoprotein hormone receptors [thyrotropin (TSHR), luteinizing hormone/chorionic gonadotropin (LHCGR), and follicle stimulating hormone (FSHR) receptors] are rhodopsin-like G protein-coupled receptors. These receptors display common structural features including a prominent extracellular domain with leucine-rich repeats (LRR) stabilized by β-sheets and a long and flexible loop known as the hinge region (HR), and a transmembrane (TM) domain with seven α-helices interconnected by intra- and extracellular loops. Binding of the ligand to the LRR resembles a hand coupling transversally to the α- and β-subunits of the hormone, with the thumb being the HR. The structure of the FSH-FSHR complex suggests an activation mechanism in which Y335 at the HR binds into a pocket between the α- and β-chains of the hormone, leading to an adjustment of the extracellular loops. In this study, we performed molecular dynamics (MD) simulations to identify the conformational changes of the FSHR and LHCGR. We set up a FSHR structure as predicted by AlphaFold (AF-P23945); for the LHCGR structure we took the cryo-electron microscopy structure for the active state (PDB:7FII) as initial coordinates. Specifically, the flexibility of the HR domain and the correlated motions of the LRR and TM domain were analyzed. From the conformational changes of the LRR, TM domain, and HR we explored the conformational landscape by means of MD trajectories in all-atom approximation, including a membrane of polyunsaturated phospholipids. The distances and procedures here defined may be useful to propose reaction coordinates to describe diverse processes, such as the active-to-inactive transition, and to identify intermediaries suited for allosteric regulation and biased binding to cellular transducers in a selective activation strategy. In the present study, we describe the results from a computational microscopy perspective (also known as molecular dynamics simulation) at the atomistic resolution for the two gonadotropin hormone receptors, the follicle-stimulant hormone receptor and the luteinizing/chorionic gonadotropin hormone receptor, which are essential for reproduction in humans. Several dysfunctional mutations in these receptors, leading to reproductive failure, have been detected in the clinical arena. To better understand the process whereby these two receptors perform their signaling tasks, we assembled the active state of the receptor structures in a membrane bilayer of phospholipids with water molecules as solvent at both sides of the membrane, and without any agonist bound to the receptor. The systems included nearly 200 thousand atoms, each moving around at 300 kelvin and 1 bar given the interactions (attractive or repulsive forces) from each other. As the motion equations are solved in each time step (at femtoseconds time scale), the system evolved over time during hundreds of nanoseconds (millions of time steps) for three independent replicates. The receptor conformation displayed non-random motions due to the stability of specific structures in the complex molecular environment, including the hydrophobic membrane core, the bilayer interfaces, and the aqueous medium. From the analysis of simulation trajectories and structural changes of the receptors, we could identify the main conformational changes exhibited by each receptor explored in a model cellular environment. We discussed the role of the hinge region at the extracellular domain in triggering the receptor conformational changes, as well as differences in the dynamics between these two receptors in terms of the flexibility of the structures. Importantly, we proposed relative distances among the different receptor domains as parameters to characterize conformational intermediaries along a transition of states. Understanding of the signaling process in gonadotropin hormone receptors might be useful to explore new strategies for the modulation of the receptor functions, the bias of signaling pathways, or the selective binding of agonists.

Impact of temperature and pH value on the stability of hGHRH: An MD approach

Comparison of the Ph- and Thermally-Induced Fluctuations of a Therapeutic Antibody Fab Fragment by Molecular Dynamics Simulation

Effect of pH on structural dynamics of HMG-CoA reductase and binding affinity to β-sitosterol

Dissecting the Effect of Temperature on Hyperthermophilic Pf2001 Esterase Dimerization by Molecular Dynamics

EFFECT OF pH ON STRUCTURE AND STABILITY OF COLLAGEN-LIKE PEPTIDE: INSIGHT FROM MOLECULAR DYNAMICS SIMULATION

Thermal stabilization of dihydrofolate reductase using monte carlo unfolding simulations and its functional consequences

Thermal Denaturing of Mutant Lysozyme with Both the OPLSAA and the CHARMM Force Fields

Molecular dynamics simulation of the effect of temperature on the conformation of ubiquitin protein

Effect of pH on the hinge region of influenza viral protein: a combined constant pH and well-tempered molecular dynamics study

Thermal Stabilities of the Microhydrated Zwitterionic Glycine: a Kinetics and Dynamics Study.

Application of conventional molecular dynamics simulation in evaluating the stability of apomyoglobin in urea solution

Constant-pH Molecular Dynamics Simulations of Closed and Open States of a Proton-gated Ion Channel

Constant pH Simulation with FMM Electrostatics in GROMACS. (A) Design and Applications

Simulation on the Ligand Leaving Process of Human Heat Shock Protein

Efffects of electrostatic and hydrophobic interaction on the stability of the TGEV main proteinase dimer

Recent Development and Application of Constant Ph Molecular Dynamics.

Temperature effects on structure and dynamics of the psychrophilic protease subtilisin S41 and its thermostable mutants in solution

Tracking conformational transitions of the gonadotropin hormone receptors in a bilayer of (SDPC) poly-unsaturated lipids from all-atom molecular dynamics simulations

A molecular dynamics study on pH response of protein adsorbed on peptide-modified polyvinyl alcohol hydrogel

Effects of the RGTFEGKF inhibitor on the structures of the transmembrane fragment 70-86 of glycophorin A: an all-atom molecular dynamics study.

Rational Design of Thermostable Carbonic Anhydrase Mutants Using Molecular Dynamics Simulations