Molecular dynamics simulation of the follicle-stimulating hormone receptor. Understanding the conformational dynamics of receptor variants at positions N680 and D408 from in silico analysis

Follicle-stimulating hormone receptor (FSHR) is a G-protein coupled receptor (GPCR) and a prototype of the glycoprotein hormone receptors subfamily of GPCRs. Structural data of the FSHR ectodomain in complex with follicle-stimulating hormone suggests a "pull and lift" activation mechanism that triggers a conformational change on the seven α-helix transmembrane domain (TMD). To analyze the conformational changes of the FSHR TMD resulting from sequence variants associated with reproductive impairment in humans, we set up a computational approach combining helix modeling and molecular simulation methods to generate conformational ensembles of the receptor at room (300 K) and physiological (310 K) temperatures. We examined the receptor dynamics in an explicit membrane environment of polyunsaturated phospholipids and solvent water molecules. The analysis of the conformational dynamics of the functional (N680 and S680) and dysfunctional (mutations at D408) variants of the FSHR allowed us to validate the FSHR-TMD model. Functional variants display a concerted motion of flexible intracellular regions at TMD helices 5 and 6. Disruption of side chain interactions and conformational dynamics were detected upon mutation at D408 when replaced with alanine, arginine, or tyrosine. Dynamical network analysis confirmed that TMD helices 2 and 5 may share communication pathways in the functional FSHR variants, whereas no connectivity was detected in the dysfunctional mutants, indicating that the global dynamics of the FSHR was sensitive to mutations at amino acid residue 408, a key position apparently linked to misfolding and variable cell surface plasma membrane expression of FSHRs with distinct mutations at this position.