Search

BackgroundThe thyroid hormone receptor-like (THR-like) family is the largest transcription factors family belonging to the nuclear receptor superfamily, which directly binds to DNA and regulates the gene expression and thereby controls various metabolic processes in a ligand-dependent manner. The THR-like family contains receptors THRs, RARs, VDR, PPARs, RORs, Rev-erbs, CAR, PXR, LXRs, and others. THR-like receptors are involved in many aspects of human health, including development, metabolism and homeostasis. Therefore, it is considered an important therapeutic target for various diseases such as osteoporosis, rickets, diabetes, etc.MethodsIn this study, we have performed an extensive sequence and structure analysis of the ligand-binding domain (LBD) of the THR-like family spanning multiple taxa. We have use different computational tools (information-theoretic measures; relative entropy) to predict the key residues responsible for fold and functional specificity in the LBD of the THR-like family. The MSA of THR-like LBDs was further used as input in conservation studies and phylogenetic clustering studies.ResultsPhylogenetic analysis of the LBD domain of THR-like proteins resulted in the clustering of eight subfamilies based on their sequence homology. The conservation analysis by relative entropy (RE) revealed that structurally important residues are conserved throughout the LBDs in the THR-like family. The multi-harmony conservation analysis further predicted specificity in determining residues in LBDs of THR-like subfamilies. Finally, fold and functional specificity determining residues (residues critical for ligand, DBD and coregulators binding) were mapped on the three-dimensional structure of thyroid hormone receptor protein. We then compiled a list of natural mutations in THR-like LBDs and mapped them along with fold and function-specific mutations. Some of the mutations were found to have a link with severe diseases like hypothyroidism, rickets, obesity, lipodystrophy, epilepsy, etc.ConclusionOur study identifies fold and function-specific residues in THR-like LBDs. We believe that this study will be useful in exploring the role of these residues in the binding of different drugs, ligands, and protein-protein interaction among partner proteins. So this study might be helpful in the rational design of either ligands or receptors.

Disorders of sex development (DSD) consist of a wide range of conditions involving numerous genes. Nevertheless, about half of 46,XY individuals remain genetically unsolved. GATA4 gene variants, mainly related to congenital heart defects (CHD), have also been recently associated with 46,XY DSD. In this study, we characterized three individuals presenting with 46,XY DSD with or without CHD and GATA4 variants in order to understand the phenotypical variability. We studied one patient presenting CHD and 46,XY gonadal dysgenesis, and two patients with a history of genetically unsolved 46,XY DSD, also known as male primary hypogonadism. Mutation analysis was carried out by candidate gene approach or targeted gene panel sequencing. Functional activity of GATA4 variants was tested in vitro on the CYP17 promoter involved in sex development using JEG3 cells. We found two novel and one previously described GATA4 variants located in the N-terminal zinc finger domain of the protein. Cys238Arg variant lost transcriptional activity on the CYP17 promoter reporter, while Trp228Cys and Pro226Leu behaved similar to wild type. These results were in line with bioinformatics simulation studies. Additional DSD variations, in the LRP4 and LHCGR genes, respectively, were identified in the two 46,XY individuals without CHD. Overall, our study shows that human GATA4 mutations identified in patients with 46,XY DSD may or may not be associated with CHD. Possible explanations for phenotypical variability may comprise incomplete penetrance, variable sensitivity of partner genes, and oligogenic mechanisms.