Modeling TSH Receptor Dimerization at the Transmembrane Domain

Abstract

Biophysical studies have established that the thyrotropin (TSH) receptor (TSHR) undergoes posttranslational modifications including dimerization. Following our earlier simulation of a TSHR-transmembrane domain (TMD) monomer (called TSHR-TMD-TRIO) we have now proceeded with a molecular dynamics simulation (MD) of TSHR-TMD dimerization using this improved membrane-embedded model. The starting structure was the TMD protein with all extracellular and intracellular loops and internal waters, which was placed in the relative orientation of the model originally generated with Brownian dynamics. Furthermore, this model was embedded in a DPPC lipid bilayer further solvated with water and added salt. Data from the MD simulation studies showed that the dimeric subunits stayed in the same relative orientation and distance during the 1000 ns of study. Comparison of representative conformations of the individual monomers when dimerized with the conformations from the monomer simulation showed subtle differences as represented by the backbone root mean square deviations. Differences in the conformations of the ligand-binding sites, suggesting variable affinities for these "hot spots," were also revealed by comparing the docking scores of 46 small-molecule ligands that included known TSHR agonists and antagonists as well as their derivatives. These data add further insight into the tendency of the TSHR-TMD to form dimeric and oligomeric structures and show that the differing conformations influence small-molecule binding sites within the TMD.

Keywords: DPPC; GPCR; TSHR; dimerization; docking; molecular dynamics.