Molecular investigation of TSHR gene in Bangladeshi congenital hypothyroid patients

The disorder of thyroid gland development or thyroid dysgenesis accounts for 80–85% of congenital hypothyroidism (CH) cases. Mutations in the TSHR gene are mostly associated with thyroid dysgenesis, and prevent or disrupt normal development of the gland. There is limited data available on the genetic spectrum of congenital hypothyroid children in Bangladesh. Thus, an understanding of the molecular aetiology of thyroid dysgenesis is a prerequisite. The aim of the study was to investigate the effect of mutations in the TSHR gene on the small molecule thyrogenic drug-binding site of the protein. We identified two nonsynonymous mutations (p.Ser508Leu, p.Glu727Asp) in the exon 10 of the TSHR gene in 21 patients with dysgenesis by sequencing-based analysis. Later, the TSHR(368-764) protein was modeled by the I-TASSER server for wild-type and mutant structures. The model proteins were targeted by thyrogenic drugs, MS437 and MS438 to perceive the effect of mutations. The damaging effect in drug-protein complexes of mutants was explored by molecular docking and molecular dynamics simulations. The binding affinity of wild-type protein was much higher than the mutant cases for both of the drug ligands (MS437 and MS438). Molecular dynamics simulates the dynamic behavior of wild-type and mutant complexes. MS437-TSHR(368-764)MT2 and MS438-TSHR(368-764)MT1 showed stable conformations in biological environments. Finally, Principle Component Analysis revealed structural and energy profile discrepancies. TSHR(368-764)MT1 exhibited much more variations than TSHR(368-764)WT and TSHR(368-764)MT2, emphasizing a more damaging pattern in TSHR(368-764)MT1. This genetic study might be helpful to explore the mutational impact on drug binding sites of TSHR protein which is important for future drug design and selection for the treatment of congenital hypothyroid children with dysgenesis.