Monte Carlo loop refinement and virtual screening of the thyroid-stimulating hormone receptor transmembrane domain

Functional Water Channels Within the TSH Receptor: A New Paradigm for TSH Action With Disease Implications

The thyroid-stimulating hormone receptor (TSHR) transmembrane domain (TMD) is found in the plasma membrane and consists of lipids and water molecules. To understand the role of TSHR-associated water molecules, we used molecular dynamic simulations of the TMD and identified a network of putative receptor-associated transmembrane water channels. This result was confirmed with extended simulations of the full-length TSHR with and without TSH ligand binding. While the transport time observed in the simulations via the TSHR protein was slower than via the lipid bilayer itself, we found that significantly more water traversed via the TSHR than via the lipid bilayer, which more than doubled with the binding of TSH. Using rat thyroid cells (FRTL-5) and a calcein fluorescence technique, we measured cell volumes after blockade of aquaporins 1 and 4, the major thyroid cell water transporters. TSH showed a dose-dependent ability to influence water transport, and similar effects were observed with stimulating TSHR autoantibodies. Small molecule TSHR agonists, which are allosteric activators of the TMD, also enhanced water transport, illustrating the role of the TMD in this phenomenon. Furthermore, the water channel pathway was also mapped across 2 activation motifs within the TSHR TMD, suggesting how water movement may influence activation of the receptor. In pathophysiological conditions such as hypothyroidism and hyperthyroidism where TSH concentrations are highly variable, this action of TSH may greatly influence water movement in thyroid cells and many other extrathyroidal sites where the TSHR is expressed, thus affecting normal cellular function.

Modeling TSH Receptor Dimerization at the Transmembrane Domain

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.

Fetoplacental vasculature as a model to study human cardiovascular endocrine disruption

Increasing evidence has associated the exposure of endocrine-disrupting chemicals (EDCs) with the cardiovascular (CV) system. This exposure is particularly problematic in a sensitive window of development, pregnancy. Pregnancy exposome can affect the overall health of the pregnancy by dramatic changes in vascular physiology and endocrine activity, increasing maternal susceptibility. Moreover, fetoplacental vascular function is generally altered, increasing the risk of developing pregnancy complications (including cardiovascular diseases, CVD) and predisposing the foetus to adverse health risks later in life. Thus, our review summarizes the existing literature on exposures to EDCs during pregnancy and adverse maternal health outcomes, focusing on the human placenta, vein, and umbilical artery associated with pregnancy complications. The purpose of this review is to highlight the role of fetoplacental vasculature as a model for the study of human cardiovascular endocrine disruption. Therefore, we emphasize that the placenta, together with the umbilical arteries and veins, allows a better characterization of the pregnant woman's exposome. Consequently, it contributes to the protection of the mother and foetus against CV disorders in life.

Thyrotropin, Hyperthyroidism, and Bone Mass

Thyrotropin (TSH), traditionally seen as a pituitary hormone that regulates thyroid glands, has additional roles in physiology including skeletal remodeling. Population-based observations in people with euthyroidism or subclinical hyperthyroidism indicated a negative association between bone mass and low-normal TSH. The findings of correlative studies were supported by small intervention trials using recombinant human TSH (rhTSH) injection, and genetic and case-based evidence. Genetically modified mouse models, which disrupt the reciprocal relationship between TSH and thyroid hormone, have allowed us to examine an independent role of TSH. Since the first description of osteoporotic phenotype in haploinsufficient Tshr +/- mice with normal thyroid hormone levels, the antiosteoclastic effect of TSH has been documented in both in vitro and in vivo studies. Further studies showed that increased osteoclastogenesis in Tshr-deficient mice was mediated by tumor necrosis factor α. Low TSH not only increased osteoclastogenesis, but also decreased osteoblastogenesis in bone marrow-derived primary osteoblast cultures. However, later in vivo studies using small and intermittent doses of rhTSH showed a proanabolic effect, which suggests that its action might be dose and frequency dependent. TSHR was shown to interact with insulin-like growth factor 1 receptor, and vascular endothelial growth factor and Wnt pathway might play a role in TSH's effect on osteoblasts. The expression and direct skeletal effect of a biologically active splice variant of the TSHβ subunit (TSHβv) in bone marrow-derived macrophage and other immune cells suggest a local skeletal effect of TSHR. Further studies of how locally secreted TSHβv and systemic TSHβ interact in skeletal remodeling through the endocrine, immune, and skeletal systems will help us better understand the hyperthyroidism-induced bone disease.

Implications of an Improved Model of the TSH Receptor Transmembrane Domain (TSHR-TMD-TRIO)

The thyroid-stimulating hormone receptor (TSHR) is a G-protein-coupled receptor group A family member with 7 transmembrane helices. We generated 3 new models of its entire transmembrane region using a 600 ns molecular simulation. The simulation started from our previously published model, which we have now revised by also modeling the intracellular loops and the C-terminal tail, adding internal waters and embedding it into a lipid bilayer with a water layer and with ions added to complete the system. We have named this model TSHR-TMD-TRIO since 3 representative dominant structures were then extracted from the simulation trajectory and compared with the original model. These structures each showed small but significant changes in the relative positions of the helices. The 3 models were also used as targets to dock a set of small molecules that are known active compounds including a new TSHR antagonist (BT362), which confirmed the appropriateness of the model with some small molecules showing significant preference for one or other of the structures.

UV-B Filter Octylmethoxycinnamate Alters the Vascular Contractility Patterns in Pregnant Women with Hypothyroidism

Increasing evidence relating the exposure and/or bioaccumulation of endocrine-disrupting compounds (EDCs) with cardiovascular system are arising. Octylmethoxycinnamate (OMC) is the most widely used UV-B filter and as EDC interacts with TH receptors. However, their effects on thyroid diseases during pregnancy remain unknown. The purpose of this work was to assess the short- and long-term effects of OMC on arterial tonus of pregnant women with hypothyroidism. To elucidate this, human umbilical artery (HUA) rings without endothelium were used to explore the vascular effects of OMC by arterial and cellular experiments. The binding energy and the modes of interaction of the OMC into the active center of the TSHR and THRα were analyzed by molecular docking studies. Our results indicated that OMC altered the contractility patterns of HUA contracted with serotonin, histamine and KCl, possibly due to an interference with serotonin and histamine receptors or an involvement of the Ca²⁺ channels. The molecular docking analysis show that OMC compete with T₃ for the binding center of THRα. Taken together, these findings pointed out to alterations in HUA reactivity as result of OMC-exposure, which may be involved in the development and increased risk of cardiovascular diseases.

A Gq Biased Small Molecule Active at the TSH Receptor

G protein coupled receptors (GPCRs) can lead to G protein and non-G protein initiated signals. By virtue of its structural property, the TSH receptor (TSHR) has a unique ability to engage different G proteins making it highly amenable to selective signaling. In this study, we describe the identification and characterization of a novel small molecule agonist to the TSHR which induces primary engagement with G_αq/11. To identify allosteric modulators inducing selective signaling of the TSHR we used a transcriptional-based luciferase assay system with CHO-TSHR cells stably expressing response elements (CRE, NFAT, SRF, or SRE) that were capable of measuring signals emanating from the coupling of G_αs , G_αq/11, G_βγ, and G_α12/13, respectively. Using this system, TSH activated G_αs , G_αq/11, and G_α12/13 but not G_βγ. On screening a library of 50K molecules at 0.1,1.0 and 10 μM, we identified a novel G_q/11 agonist (named MSq1) which activated G_q/11 mediated NFAT-luciferase >4 fold above baseline and had an EC₅₀= 8.3 × 10^-9 M with only minor induction of G_αs and cAMP. Furthermore, MSq1 is chemically and structurally distinct from any of the previously reported TSHR agonist molecules. Docking studies using a TSHR transmembrane domain (TMD) model indicated that MSq1 had contact points on helices H1, H2, H3, and H7 in the hydrophobic pocket of the TMD and also with the extracellular loops. On co-treatment with TSH, MSq1 suppressed TSH-induced proliferation of thyrocytes in a dose-dependent manner but lacked the intrinsic ability to influence basal thyrocyte proliferation. This unexpected inhibitory property of MSq1 could be blocked in the presence of a PKC inhibitor resulting in derepressing TSH induced protein kinase A (PKA) signals and resulting in the induction of proliferation. Thus, the inhibitory effect of MSq1 on proliferation resided in its capacity to overtly activate protein kinase C (PKC) which in turn suppressed the proliferative signal induced by activation of the predomiant cAMP-PKA pathway of the TSHR. Treatment of rat thyroid cells (FRTL5) with MSq1 did not show any upregulation of gene expression of the key thyroid specific markers such as thyroglobulin(Tg), thyroid peroxidase (Tpo), sodium iodide symporter (Nis), and the TSH receptor (Tshr) further suggesting lack of involvement of MSq1 and G_αq/11 activation with cellular differentation. In summary, we identified and characterized a novel G_αq/11 agonist molecule acting at the TSHR and which showed a marked anti-proliferative ability. Hence, Gq biased activation of the TSHR is capable of ameliorating the proliferative signals from its orthosteric ligand and may offer a therapeutic option for thyroid growth modulation.

Born This Way: Using Intrinsic Disorder to Map the Connections between SLITRKs, TSHR, and Male Sexual Orientation

Recently, genome-wide association study reveals a significant association between specific single nucleotide polymorphisms (SNPs) in men and their sexual orientation. These SNPs (rs9547443 and rs1035144) reside in the intergenic region between the SLITRK5 and SLITRK6 genes and in the intronic region of the TSHR gene and might affect functionality of SLITRK5, SLITRK6, and TSHR proteins that are engaged in tight control of key developmental processes, such as neurite outgrowth and modulation, cellular differentiation, and hormonal regulation. SLITRK5 and SLITRK6 are single-pass transmembrane proteins, whereas TSHR is a heptahelical G protein-coupled receptor (GPCR). Mutations in these proteins are associated with various diseases and are linked to phenotypes found at a higher rate in homosexual men. A bioinformatics analysis of SLITRK5, SLITRK6, and TSHR proteins is conducted to look at their structure, protein interaction networks, and propensity for intrinsic disorder. It is assumed that this information might improve understanding of the roles that SLITRK5, SLITRK6, and TSHR play within neuronal and thyroidal tissues and give insight into the phenotypes associated with male homosexuality.

Identification of a Novel Class of BRD4 Inhibitors by Computational Screening and Binding Simulations

Computational screening is a method to prioritize small-molecule compounds based on the structural and biochemical attributes built from ligand and target information. Previously, we have developed a scalable virtual screening workflow to identify novel multitarget kinase/bromodomain inhibitors. In the current study, we identified several novel N-[3-(2-oxo-pyrrolidinyl)phenyl]-benzenesulfonamide derivatives that scored highly in our ensemble docking protocol. We quantified the binding affinity of these compounds for BRD4(BD1) biochemically and generated cocrystal structures, which were deposited in the Protein Data Bank. As the docking poses obtained in the virtual screening pipeline did not align with the experimental cocrystal structures, we evaluated the predictions of their precise binding modes by performing molecular dynamics (MD) simulations. The MD simulations closely reproduced the experimentally observed protein-ligand cocrystal binding conformations and interactions for all compounds. These results suggest a computational workflow to generate experimental-quality protein-ligand binding models, overcoming limitations of docking results due to receptor flexibility and incomplete sampling, as a useful starting point for the structure-based lead optimization of novel BRD4(BD1) inhibitors.

Structural-Functional Features of the Thyrotropin Receptor: A Class A G-Protein-Coupled Receptor at Work

The thyroid-stimulating hormone receptor (TSHR) is a member of the glycoprotein hormone receptors, a sub-group of class A G-protein-coupled receptors (GPCRs). TSHR and its endogenous ligand thyrotropin (TSH) are of essential importance for growth and function of the thyroid gland and proper function of the TSH/TSHR system is pivotal for production and release of thyroid hormones. This receptor is also important with respect to pathophysiology, such as autoimmune (including ophthalmopathy) or non-autoimmune thyroid dysfunctions and cancer development. Pharmacological interventions directly targeting the TSHR should provide benefits to disease treatment compared to currently available therapies of dysfunctions associated with the TSHR or the thyroid gland. Upon TSHR activation, the molecular events conveying conformational changes from the extra- to the intracellular side of the cell across the membrane comprise reception, conversion, and amplification of the signal. These steps are highly dependent on structural features of this receptor and its intermolecular interaction partners, e.g., TSH, antibodies, small molecules, G-proteins, or arrestin. For better understanding of signal transduction, pathogenic mechanisms such as autoantibody action and mutational modifications or for developing new pharmacological strategies, it is essential to combine available structural data with functional information to generate homology models of the entire receptor. Although so far these insights are fragmental, in the past few decades essential contributions have been made to investigate in-depth the involved determinants, such as by structure determination via X-ray crystallography. This review summarizes available knowledge (as of December 2016) concerning the TSHR protein structure, associated functional aspects, and based on these insights we suggest several receptor complex models. Moreover, distinct TSHR properties will be highlighted in comparison to other class A GPCRs to understand the molecular activation mechanisms of this receptor comprehensively. Finally, limitations of current knowledge and lack of information are discussed highlighting the need for intensified efforts toward TSHR structure elucidation.

Structure and activation of the TSH receptor transmembrane domain

PURPOSE

The thyroid-stimulating hormone receptor (TSHR) is the target autoantigen for TSHR-stimulating autoantibodies in Graves' disease. The TSHR is composed of: a leucine-rich repeat domain (LRD), a hinge region or cleavage domain (CD) and a transmembrane domain (TMD). The binding arrangements between the TSHR LRD and the thyroid-stimulating autoantibody M22 or TSH have become available from the crystal structure of the TSHR LRD-M22 complex and a comparative model of the TSHR LRD in complex with TSH, respectively. However, the mechanism by which the TMD of the TSHR and the other glycoprotein hormone receptors (GPHRs) becomes activated is unknown.

METHODS

We have generated comparative models of the structures of the inactive (TMD_In) and active (TMD_Ac) conformations of the TSHR, follicle-stimulating hormone receptor (FSHR) and luteinizing hormone receptor (LHR) TMDs. The structures of TMD_Ac and TMD_In were obtained using class A GPCR crystal structures for which fully active and inactive conformations were available.

RESULTS

Most conserved motifs observed in GPCR TMDs are also observed in the amino acid sequences of GPHR TMDs. Furthermore, most GPCR TMD conserved helix distortions are observed in our models of the structures of GPHR TMDs. Analysis of these structures has allowed us to propose a mechanism for activation of GPHR TMDs.

CONCLUSIONS

Insight into the mechanism of activation of the TSHR by both TSH and TSHR autoantibodies is likely to be useful in the development of new treatments for Graves' disease.

TSH Receptor Signaling Abrogation by a Novel Small Molecule

Pathological activation of the thyroid-stimulating hormone receptor (TSHR) is caused by thyroid-stimulating antibodies in patients with Graves' disease (GD) or by somatic and rare genomic mutations that enhance constitutive activation of the receptor influencing both G protein and non-G protein signaling. Potential selective small molecule antagonists represent novel therapeutic compounds for abrogation of such abnormal TSHR signaling. In this study, we describe the identification and in vitro characterization of a novel small molecule antagonist by high-throughput screening (HTS). The identification of the TSHR antagonist was performed using a transcription-based TSH-inhibition bioassay. TSHR-expressing CHO cells, which also expressed a luciferase-tagged CRE response element, were optimized using bovine TSH as the activator, in a 384 well plate format, which had a Z score of 0.3-0.6. Using this HTS assay, we screened a diverse library of ~80,000 compounds at a final concentration of 16.7 μM. The selection criteria for a positive hit were based on a mean signal threshold of ≥50% inhibition of control TSH stimulation. The screening resulted in 450 positive hits giving a hit ratio of 0.56%. A secondary confirmation screen against TSH and forskolin - a post receptor activator of adenylyl cyclase - confirmed one TSHR-specific candidate antagonist molecule (named VA-K-14). This lead molecule had an IC₅₀ of 12.3 μM and a unique chemical structure. A parallel analysis for cell viability indicated that the lead inhibitor was non-cytotoxic at its effective concentrations. In silico docking studies performed using a TSHR transmembrane model showed the hydrophobic contact locations and the possible mode of inhibition of TSHR signaling. Furthermore, this molecule was capable of inhibiting TSHR stimulation by GD patient sera and monoclonal-stimulating TSHR antibodies. In conclusion, we report the identification of a novel small molecule TSHR inhibitor, which has the potential to be developed as a therapeutic antagonist for abrogation of TSHR signaling by TSHR autoantibodies in GD.

Transmembrane domains of attraction on the TSH receptor

The TSH receptor (TSHR) has the propensity to form dimers and oligomers. Our data using ectodomain-truncated TSHRs indicated that the predominant interfaces for oligomerization reside in the transmembrane (TM) domain. To map the potentially interacting residues, we first performed in silico studies of the TSHR transmembrane domain using a homology model and using Brownian dynamics (BD). The cluster of dimer conformations obtained from BD analysis indicated that TM1 made contact with TM4 and two residues in TM2 made contact with TM5. To confirm the proximity of these contact residues, we then generated cysteine mutants at all six contact residues predicted by the BD analysis and performed cysteine cross-linking studies. These results showed that the predicted helices in the protomer were indeed involved in proximity interactions. Furthermore, an alternative experimental approach, receptor truncation experiments and LH receptor sequence substitution experiments, identified TM1 harboring a major region involved in TSHR oligomerization, in agreement with the conclusion from the cross-linking studies. Point mutations of the predicted interacting residues did not yield a substantial decrease in oligomerization, unlike the truncation of the TM1, so we concluded that constitutive oligomerization must involve interfaces forming domains of attraction in a cooperative manner that is not dominated by interactions between specific residues.

New small molecule agonists to the thyrotropin receptor

BACKGROUND

Novel small molecular ligands (SMLs) to the thyrotropin receptor (TSHR) have potential as improved molecular probes and as therapeutic agents for the treatment of thyroid dysfunction and thyroid cancer.

METHODS

To identify novel SMLs to the TSHR, we developed a transcription-based luciferase-cAMP high-throughput screening system and we screened 48,224 compounds from a 100K library in duplicate.

RESULTS

We obtained 62 hits using the cut-off criteria of the mean±three standard deviations above the baseline. Twenty molecules with the greatest activity were rescreened against the parent CHO-luciferase cell for nonspecific activation, and we selected two molecules (MS437 and MS438) with the highest potency for further study. These lead molecules demonstrated no detectible cross-reactivity with homologous receptors when tested against luteinizing hormone (LH)/human chorionic gonadotropin receptor and follicle stimulating hormone receptor-expressing cells. Molecule MS437 had a TSHR-stimulating potency with an EC50 of 13×10(-8) M, and molecule MS438 had an EC50 of 5.3×10(-8) M. The ability of these small molecule agonists to bind to the transmembrane domain of the receptor and initiate signal transduction was suggested by their activation of a chimeric receptor consisting of an LHR ectodomain and a TSHR transmembrane. Molecular modeling demonstrated that these molecules bound to residues S505 and E506 for MS438 and T501 for MS437 in the intrahelical region of transmembrane helix 3. We also examined the G protein activating ability of these molecules using CHO cells co-expressing TSHRs transfected with luciferase reporter vectors in order to measure Gsα, Gβγ, Gαq, and Gα12 activation quantitatively. The MS437 and MS438 molecules showed potent activation of Gsα, Gαq, and Gα12 similar to TSH, but neither the small molecule agonists nor TSH showed activation of the Gβγ pathway. The small molecules MS437 and MS438 also showed upregulation of thyroglobulin (Tg), sodium iodine symporter (NIS), and TSHR gene expression.

CONCLUSIONS

Pharmacokinetic analysis of MS437 and MS438 indicated their pharmacotherapeutic potential, and their intraperitoneal administration to normal female mice resulted in significantly increased serum thyroxine levels, which could be maintained by repeated treatments. These molecules can therefore serve as lead molecules for further development of powerful TSH agonists.

Expanding the horizons of G protein-coupled receptor structure-based ligand discovery and optimization using homology models

We show the key role of structural homology models in GPCR structure-based lead discovery and optimization, highlighting methodological aspects, recent progress and future directions.