Only subtle protein conformational adaptations are required for ligand binding to thyroid hormone receptors: Simulations using a novel multipoint steered molecular dynamics approach

Thyroid hormone receptors (TR) are hormone-dependent transcription regulators that play a major role in human health, development, and metabolic functions. The thyroid hormone resistance syndrome, diabetes, obesity, and some types of cancer are just a few examples of important diseases that are related to TR malfunctioning, particularly impaired hormone binding. Ligand binding to and dissociation from the receptor ultimately control gene transcription and, thus, detailed knowledge of binding and release mechanisms are fundamental for the comprehension of the receptor`s biological function and development of pharmaceuticals. In this work, we present the first computational study of ligand entry into the ligand binding domain (LBD) of a nuclear receptor. We report molecular dynamics simulations of ligand binding to TRs using a generalization of the steered molecular dynamics technique designed to perform single-molecule pulling simulations along arbitrarily nonlinear driving pathways. We show that only gentle protein movements and conformational adaptations are required for ligand entry into the LBDs and that the magnitude of the forces applied to assist ligand binding are of the order of the forces involved in ligand dissociation. Our simulations suggest an alternative view for the mechanisms ligand binding and dissociation of ligands from nuclear receptors in which ligands can simply diffuse through the protein surface to reach proper positioning within the binding pocket. The proposed picture indicates that the large-amplitude protein motions suggested by the apo- and holo-RXR alpha crystallographic structures are not required, reconciling conformational changes of LBDs required for ligand entry with other nuclear receptors apo-structures that resemble the ligand-bound LBDs.

A novel locus for split-hand/foot malformation associated with tibial hemimelia (SHFLD syndrome) maps to chromosome region 17p13.1-17p13.3

Split-hand/foot malformation (SHFM) associated with aplasia of long bones, SHFLD syndrome or Tibial hemimelia-ectrodactyly syndrome is a rare condition with autosomal dominant inheritance, reduced penetrance and an incidence estimated to be about 1 in 1,000,000 liveborns. To date, three chromosomal regions have been reported as strong candidates for harboring SHFLD syndrome genes: 1q42.2-q43, 6q14.1 and 2q14.2. We characterized the phenotype of nine affected individuals from a large family with the aim of mapping the causative gene. Among the nine affected patients, four had only SHFM of the hands and no tibial defects, three had both defects and two had only unilateral tibial hemimelia. In keeping with previous publications of this and other families, there was clear evidence of both variable expression and incomplete penetrance, the latter bearing hallmarks of anticipation. Segregation analysis and multipoint Lod scores calculations (maximum Lod score of 5.03 using the LINKMAP software) using all potentially informative family members, both affected and unaffected, identified the chromosomal region 17p13.1-17p13.3 as the best and only candidate for harboring a novel mutated gene responsible for the syndrome in this family. The candidate gene CRK located within this region was sequenced but no pathogenic mutation was detected.