The PAX8/PPARγ rearrangement, producing the PAX8–PPARγ fusion protein (PPFP), is thought to play an essential role in the oncogenesis of thyroid follicular tumors. To identify PPFP-targeted drug candidates and establish an early standard of care for thyroid tumors, we performed ensemble docking-based compound screening. Specifically, we investigated which pocket structure should be adopted to search for a promising ligand compound for the PPFP; the position of the ligand binding pocket on the PPARγ side of the PPFP is common to that of PPARγ; however, the shape is slightly different between them due to environmental factors. We developed a method for selecting a PPFP structure with a relevant pocket and high prediction accuracy for ligand binding. This method was validated using PPARγ, whose structure and activity values are known for many compounds. Then, we performed docking calculations to PPFP for 97 drug or drug-like compounds with a thiazolidine backbone, which is one of the characteristics of ligands that bind well to PPARγ, registered in the DrugBank database. Furthermore, the binding affinities of promising ligand candidates were estimated more reliably using the molecular mechanics-Poisson–Boltzmann surface area method. Thus, we propose promising drug candidates for the PPFP with a thiazolidine backbone.

Mexico shows a high birth prevalence of congenital hypothyroidism (CH) due to thyroid dysgenesis (TD). PAX8 defects underlie only 1% of these cases and NKX2-1 does not seem to be involved. Here, we analyzed other TD-related genes in 128 non-related Mexican patients (females 77.3%; 6 months to 16.6 years) with non-syndromic CH-TD diagnosis established by clinical evaluation, thyroid hormone serum profiling, and scintigraphy (74%) or ultrasonography (26%). We performed Sanger sequencing of FOXE1, NKX2-5, and TSHR and evaluated copy number variations (CNVs) in TSHR, FOXE1, PAX8, and NKX2-1 by multiplex ligation-dependent probe amplification. Odds ratios for TD risk were explored for FOXE1 polyalanine stretches [polyAla-rs71369530] in cases and controls (N=116). Five rare missense changes cataloged as benign (NKX2-5:p.(Ala119Ser)-rs137852684), of unknown significance (FOXE1:p.(Ala335Gly)-rs543372757; TSHR:p.(Asp118Asn)-rs1414102266), and likely pathogenic (FOXE1:p.(Gly124Arg)-rs774035532; TSHR:p.(Trp422Arg)-rs746029360) accounted for 1.5% (N=2/128) of clinically relevant genotypes (supported in part by protein modeling) in CH-TD. No CNVs were identified, nor did polyAla >14 alanines in FOXE1 significantly protect against TD. The present and previously published data collectively show that small clinically relevant germline variants in PAX8, FOXE1, and TSHR are found in only a very small proportion (2.5%) of isolated CH-TD Mexican patients.