Structure of a HoxB1-Pbx1 heterodimer bound to DNA: role of the hexapeptide and a fourth homeodomain helix in complex formation

The intrinsically disordered transactivation region of HOXA9 regulates its function by auto-inhibition of its DNA-binding activity

HOXA9 transcription factor is expressed in hematopoietic stem cells and is involved in the regulation of their differentiation and maturation to various blood cells. HOXA9 is linked to various leukemia and is a marker for poor prognosis of acute myeloid leukemia (AML). This protein has a conserved DNA-binding homeodomain and a transactivation domain. We show that this N-terminal transactivation domain is intrinsically disordered and inhibits DNA-binding by the homeodomain. Using NMR spectroscopy and molecular dynamics simulation, we show that the hexapeptide 197AANWLH202 in the disordered region transiently occludes the DNA-binding interface. The hexapeptide also forms a rigid segment, as determined by NMR dynamics, in an otherwise flexible disordered region. Interestingly, this hexapeptide is known to mediate the interaction of HOXA9 and its TALE partner proteins, such as PBX1, and help in cooperative DNA binding. Mutation of tryptophan to alanine in the hexapeptide abrogates the DNA-binding auto-inhibition. We propose that the disordered transactivation region plays a dual role in the regulation of HOXA9 function. In the absence of TALE partners, it inhibits DNA binding, and in the presence of TALE partners it interacts with the TALE protein and facilitates the cooperative DNA binding by the HOX-TALE complex.

TALE transcription factors: Cofactors no more

Exd/PBX, Hth/MEIS and PREP proteins belong to the TALE (three-amino-acid loop extension) superclass of transcription factors (TFs) with an atypical homedomain (HD). Originally discovered as "cofactors" to HOX proteins, revisiting their traditional role in light of genome-wide experiments reveals a strong and reproducible pattern of HOX and TALE co-occupancy across diverse embryonic tissues. While confirming that TALE increases HOX specificity and selectivity in vivo, this wider outlook also reveals novel aspects of HOX:TALE collaboration, namely that HOX TFs generally require pre-bound TALE factors to access their functional binding sites in vivo. In contrast to the restricted expression domains of HOX TFs, TALE factors are largely ubiquitous, and PBX and PREP are expressed at the earliest developmental stages. PBX and MEIS control development of many organs and tissues and their dysregulation is associated with congenital disease and cancer. Accordingly, many instances of TALE cooperation with non HOX TFs have been documented in various systems. The model that emerges from these studies is that TALE TFs create a permissive chromatin platform that is selected by tissue-restricted TFs for binding. In turn, HOX and other tissue-restricted TFs selectively convert a ubiquitous pool of low affinity TALE binding events into high confidence, tissue-restricted binding events associated with transcriptional activation. As a result, TALE:TF complexes are associated with active chromatin and domain/lineage-specific gene activity. TALE ubiquitous expression and broad genomic occupancy, as well as the increasing examples of TALE tissue-specific partners, reveal a universal and obligatory role for TALE in the control of tissue and lineage-specific transcriptional programs, beyond their initial discovery as HOX co-factors.

The advances of E2A-PBX1 fusion in B-cell acute lymphoblastic Leukaemia

The E2A-PBX1 gene fusion is a common translocation in B-cell acute lymphoblastic leukaemia. Patients harbouring the E2A-PBX1 fusion gene typically exhibit an intermediate prognosis. Furthermore, minimal residual disease has unsatisfactory prognostic value in E2A-PBX1 B-cell acute lymphoblastic leukaemia. However, the mechanism of E2A-PBX1 in the occurrence and progression of B-cell acute lymphoblastic leukaemia is not well understood. Here, we mainly review the roles of E2A and PBX1 in the differentiation and development of B lymphocytes, the mechanism of E2A-PBX1 gene fusion in B-cell acute lymphoblastic leukaemia, and the potential therapeutic approaches.

Ramifications of POU4F3 variants associated with autosomal dominant hearing loss in various molecular aspects

POU4F3, a member of the POU family of transcription factors, commonly causes autosomal dominant deafness. Exome sequencing was used to identify four novel variants in POU4F3 (NM_002700.2), including c.564dupA: p.Ala189SerfsTer26, c.743T > C:p.Leu248Pro, c.879C > A:p.Phe293Leu, and c.952G > A:p.Val318Met, and diverse aspects of the molecular consequences of their protein expression, stability, subcellular localization, and transcriptional activity were investigated. The expression of three mutant proteins, encoded by missense variants, was reduced compared to the wild-type protein, demonstrating that the mutants were unstable and vulnerable to degradation. Additionally, all the mutant proteins had distinct subcellular localization patterns. A mutant protein carrying p.Ala189SerfsTer26, in which both mono- and bi-partite nuclear localization signals were disrupted, showed abnormal subcellular localization. Resultantly, all the mutant proteins significantly reduced the transcriptional activity required to regulate the downstream target gene expression. Furthermore, we identified the altered expression of 14 downstream target genes associated with inner ear development using patient-derived lymphoblastoid cell lines. There was a significant correlation of the expression profile between patient-derived cells and the cochlear hair cells, which provided a breakthrough for cases where the collection of human cochlear samples for transcriptome studies was unfeasible. This study expanded the genotypic spectrum of POU4F3 in DFNA15, and further refined the molecular mechanisms underlying POU4F3-associated DFNA15.

NANOGP1, a tandem duplicate of NANOG, exhibits partial functional conservation in human naïve pluripotent stem cells

Gene duplication events can drive evolution by providing genetic material for new gene functions, and they create opportunities for diverse developmental strategies to emerge between species. To study the contribution of duplicated genes to human early development, we examined the evolution and function of NANOGP1, a tandem duplicate of the transcription factor NANOG. We found that NANOGP1 and NANOG have overlapping but distinct expression profiles, with high NANOGP1 expression restricted to early epiblast cells and naïve-state pluripotent stem cells. Sequence analysis and epitope-tagging revealed that NANOGP1 is protein coding with an intact homeobox domain. The duplication that created NANOGP1 occurred earlier in primate evolution than previously thought and has been retained only in great apes, whereas Old World monkeys have disabled the gene in different ways, including homeodomain point mutations. NANOGP1 is a strong inducer of naïve pluripotency; however, unlike NANOG, it is not required to maintain the undifferentiated status of human naïve pluripotent cells. By retaining expression, sequence and partial functional conservation with its ancestral copy, NANOGP1 exemplifies how gene duplication and subfunctionalisation can contribute to transcription factor activity in human pluripotency and development.

Novel Molecular Genetic Etiology of Asymmetric Hearing Loss: Autosomal-Dominant LMX1A Variants

INTRODUCTION

Sensorineural hearing loss is the most common sensory disorder in humans. Genetic analyses have greatly increased our understanding of the pathogenic mechanisms in play. Thus, characterization of audiologic phenotypes by the genetic etiology may aid elucidation of the etiologies of certain types of inherited hearing loss. Further, delineation of specific audiologic phenotypes based on the genetic etiology aids our understanding of some types of inherited hearing loss in terms of the prediction of clinical course, revelation of genotype-phenotype correlations, and application of appropriate audiologic rehabilitation. Here, we describe the interesting audiologic characteristics of LMX1A -associated deafness, which revealed significant asymmetry between two ears.

METHODS

Among 728 probands of which genomic DNA went through exome sequencing regardless of any specific audiologic phenotypes, probands for which exome sequencing was performed and a causative LMX1A variant was found were all included. Five LMX1A -associated DFNA7 families (approximately 0.7%), the pedigrees of whom indicated autosomal-dominant hearing loss, were identified, and segregation was studied using Sanger sequencing. The affected individuals underwent comprehensive evaluations, including medical history reviews, physical examinations, imaging, and auditory phenotyping. We functionally characterized the novel LMX1A variants via computational structural modeling and luciferase reporter assays.

RESULTS

Among 728 probands of which genomic DNA went through exome sequencing, we identified four novel LMX1A heterozygous variants related to DFNA7 (c.622C>T:p.Arg208*, c.719A>G:p.Gln240Arg, c.721G>A:p.Val241Met, and c.887dup:p.Gln297Thrfs*41) and one harboring a de novo heterozygous missense LMX1A variant (c.595A>G;p.Arg199Gly) previously reported. It is important to note that asymmetric hearing loss was identified in all probands and most affected individuals, although the extent of asymmetry varied. Structural modeling revealed that the two missense variants, p.Gln240Arg and p.Val241Met, affected conserved residues of the homeodomain, thus attenuating LMX1A-DNA interaction. In addition, Arg208*-induced premature termination of translation destroyed the structure of the LMX1A protein, including the DNA-binding homeodomain, and p.Gln297Thrfs*41 led to the loss of the C-terminal helix involved in LIM2 domain interaction. Compared with the wild-type protein, all mutant LMX1A proteins had significantly reduced transactivation efficiency, indicating that the ability to elicit transcription of the downstream target genes of LMX1A was severely compromised. Thus, in line with the American College of Medical Genetics and Genomics guideline specified to genetic hearing loss, the four novel LMX1A variants were identified as "pathogenic" (p.Arg208* and p.Gln297Thrfs*41), "likely pathogenic" (p.Val241Met), and as a "variant of uncertain significance'' (p.Gln240Arg).

CONCLUSION

For the first time, we suggest that LMX1A is one of the candidate genes which, if altered, could be associated with dominantly inherited asymmetric hearing loss. We also expand the genotypic spectrum of disease-causing variants of LMX1A causing DFNA7 by doubling the number of LMX1A variants reported thus far in the literature.

HOXA1, a breast cancer oncogene

More than 25 years ago, the first literature records mentioned HOXA1 expression in human breast cancer. A few years later, HOXA1 was confirmed as a proper oncogene in mammary tissue. In the following two decades, molecular data about the mode of action of the HOXA1 protein, the factors contributing to activate and maintain HOXA1 gene expression and the identity of its target genes have accumulated and provide a wider view on the association of this transcription factor to breast oncogenesis. Large-scale transcriptomic data gathered from wide cohorts of patients further allowed refining the relationship between breast cancer type and HOXA1 expression. Several recent reports have reviewed the connection between cancer hallmarks and the biology of HOX genes in general. Here we take HOXA1 as a paradigm and propose an extensive overview of the molecular data centered on this oncoprotein, from what its expression modulators, to the interactors contributing to its oncogenic activities, and to the pathways and genes it controls. The data converge to an intricate picture that answers questions on the multi-modality of its oncogene activities, point towards better understanding of breast cancer aetiology and thereby provides an appraisal for treatment opportunities.

Diversification and Functional Evolution of HOX Proteins

Gene duplication and divergence is a major contributor to the generation of morphological diversity and the emergence of novel features in vertebrates during evolution. The availability of sequenced genomes has facilitated our understanding of the evolution of genes and regulatory elements. However, progress in understanding conservation and divergence in the function of proteins has been slow and mainly assessed by comparing protein sequences in combination with in vitro analyses. These approaches help to classify proteins into different families and sub-families, such as distinct types of transcription factors, but how protein function varies within a gene family is less well understood. Some studies have explored the functional evolution of closely related proteins and important insights have begun to emerge. In this review, we will provide a general overview of gene duplication and functional divergence and then focus on the functional evolution of HOX proteins to illustrate evolutionary changes underlying diversification and their role in animal evolution.

Identification and functional characterization of transcriptional activators in human cells

Transcription is orchestrated by thousands of transcription factors (TFs) and chromatin-associated proteins, but how these are causally connected to transcriptional activation is poorly understood. Here, we conduct an unbiased proteome-scale screen to systematically uncover human proteins that activate transcription in a natural chromatin context. By combining interaction proteomics and chemical inhibitors, we delineate the preference of these transcriptional activators for specific co-activators, highlighting how even closely related TFs can function via distinct cofactors. We also identify potent transactivation domains among the hits and use AlphaFold2 to predict and experimentally validate interaction interfaces of two activation domains with BRD4. Finally, we show that many novel activators are partners in fusion events in tumors and functionally characterize a myofibroma-associated fusion between SRF and C3orf62, a potent p300-dependent activator. Our work provides a functional catalog of potent transactivators in the human proteome and a platform for discovering transcriptional regulators at genome scale.

PBX1: a key character of the hallmarks of cancer

Pre-B-cell leukemia homeobox transcription factor 1 (PBX1) was first identified as part of a fusion protein resulting from the chromosomal translocation t(1;19) in pre-B cell acute lymphoblastic leukemias. Since then, PBX1 has been associated with important developmental programs, and its expression dysregulation has been related to multifactorial disorders, including cancer. As PBX1 overexpression in many cancers is correlated to poor prognosis, we sought to understand how this transcription factor contributes to carcinogenesis, and to organize PBX1's roles in the hallmarks of cancer. There is enough evidence to associate PBX1 with at least five hallmarks: sustaining proliferative signaling, activating invasion and metastasis, inducing angiogenesis, resisting cell death, and deregulating cellular energetics. The lack of studies investigating a possible role for PBX1 on the remaining hallmarks made it impossible to defend or refute its contribution on them. However, the functions of some of the PBX1's transcription targets indicate a potential engagement of PBX1 in the avoidance of immune destruction and in the tumor-promoting inflammation hallmarks. Interestingly, PBX1 might be a player in tumor suppression by activating the transcription of some DNA damage response genes. This is the first review organizing PBX1 roles into the hallmarks of cancer. Thus, we encourage future studies to uncover the PBX1's underlying mechanisms to promote carcinogenesis, for it is a promising diagnostic and prognostic biomarker, as well as a potential target in cancer treatment.

Development of small molecule inhibitors targeting PBX1 transcription signaling as a novel cancer therapeutic strategy

PBX1 is a transcription factor involved in diverse cellular functions including organ development, stem cell renewal, and tumorigenesis. PBX1 is localized at chr1q23.3, a frequently amplified chromosomal region, and it is overexpressed in many human malignancies. Cancer cells with elevated PBX1 signaling are particularly vulnerable to PBX1 withdrawal. We designed a series of small molecule compounds capable of docking to the interface between PBX1 and its cognate DNA target sequence. Among them, T417 is found to be a lead compound. In cell-based assays, T417 significantly suppressed self-renewal and proliferation of cancer cells expressing high levels of PBX1. T417 also re-sensitized platinum-resistant ovarian tumors to carboplatin. T417 did not affect healthy tissues likely due to their lower PBX1 expression levels. Therefore, targeting PBX-DNA interface can be a promising strategy for treating human tumors reliant on PBX1 for survival.

The Hox protein conundrum: The "specifics" of DNA binding for Hox proteins and their partners

Homeotic genes (Hox genes) are homeodomain-transcription factors involved in conferring segmental identity along the anterior-posterior body axis. Molecular characterization of HOX protein function raises some interesting questions regarding the source of the binding specificity of the HOX proteins. How do HOX proteins regulate common and unique target specificity across space and time? This review attempts to summarize and interpret findings in this area, largely focused on results from in vitro and in vivo studies in Drosophila and mouse systems. Recent studies related to HOX protein binding specificity compel us to reconsider some of our current models for transcription factor-DNA interactions. It is crucial to study transcription factor binding by incorporating components of more complex, multi-protein interactions in concert with small changes in binding motifs that can significantly impact DNA binding specificity and subsequent alterations in gene expression. To incorporate the multiple elements that can determine HOX protein binding specificity, we propose a more integrative Cooperative Binding model.

A Computer-Based Methodology to Design Non-Standard Peptides Potentially Able to Prevent HOX-PBX1-Associated Cancer Diseases

In the last decades, HOX proteins have been extensively studied due to their pivotal role in transcriptional events. HOX proteins execute their activity by exploiting a cooperative binding to PBX proteins and DNA. Therefore, an increase or decrease in HOX activity has been associated with both solid and haematological cancer diseases. Thus, inhibiting HOX-PBX interaction represents a potential strategy to prevent these malignancies, as demonstrated by the patented peptide HTL001 that is being studied in clinical trials. In this work, a computational study is described to identify novel potential peptides designed by employing a database of non-natural amino acids. For this purpose, residue scanning of the HOX minimal active sequence was performed to select the mutations to be further processed. According to these results, the peptides were point-mutated and used for Molecular Dynamics (MD) simulations in complex with PBX1 protein and DNA to evaluate complex binding stability. MM-GBSA calculations of the resulting MD trajectories were exploited to guide the selection of the most promising mutations that were exploited to generate twelve combinatorial peptides. Finally, the latter peptides in complex with PBX1 protein and DNA were exploited to run MD simulations and the ΔGbinding average values of the complexes were calculated. Thus, the analysis of the results highlighted eleven combinatorial peptides that will be considered for further assays.

How structural biology transformed studies of transcription regulation

The past 4 decades have seen remarkable advances in our understanding of the structural basis of gene regulation. Technological advances in protein expression, nucleic acid synthesis, and structural biology made it possible to study the proteins that regulate transcription in the context of ever larger complexes containing proteins bound to DNA. This review, written on the occasion of the 50th anniversary of the founding of the Protein Data Bank focuses on the insights gained from structural studies of protein-DNA complexes and the role the PDB has played in driving this research. I cover highlights in the field, beginning with X-ray crystal structures of the first DNA-binding domains to be studied, through recent cryo-EM structures of transcription factor binding to nucleosomal DNA.

Posttranslational Modifications in Conserved Transcription Factors: A Survey of the TALE-Homeodomain Superclass in Human and Mouse

Transcription factors (TFs) guide effector proteins like chromatin-modifying or -remodeling enzymes to distinct sites in the genome and thereby fulfill important early steps in translating the genome’s sequence information into the production of proteins or functional RNAs. TFs of the same family are often highly conserved in evolution, raising the question of how proteins with seemingly similar structure and DNA-binding properties can exert physiologically distinct functions or respond to context-specific extracellular cues. A good example is the TALE superclass of homeodomain-containing proteins. All TALE-homeodomain proteins share a characteristic, 63-amino acid long homeodomain and bind to similar sequence motifs. Yet, they frequently fulfill non-redundant functions even in domains of co-expression and are subject to regulation by different signaling pathways. Here we provide an overview of posttranslational modifications that are associated with murine and human TALE-homeodomain proteins and discuss their possible importance for the biology of these TFs.

The Increased Expression of an Engrailed to Sustain Shell Formation in Response to Ocean Acidification

Engrailed is a transcription factor required in numerous species for important developmental steps such as neurogenesis, segment formation, preblastoderm organization, and compartment formation. Recent study has proved that engrailed is also a key gene related to shell formation in marine bivalves. In the present study, the expression pattern of an engrailed gene (Cgengrailed-1) in Pacific oyster Crassostrea gigas under CO₂-driven acidification was investigated to understand its possible role in the regulation of shell formation and adaptation to ocean acidification (OA). The open reading frame (ORF) of Cgengrailed-1 was obtained, which was of 690 bp encoding a polypeptide of 229 amino acids with a HOX domain. Phylogenetic analysis indicated that the deduced amino acid sequence of Cgengrailed-1 shared high homology with other engraileds from Drosophila melanogaster, Mizuhopecten yessoensi, and Crassostrea virginica. The mRNA transcripts of Cgengrailed-1 were constitutively expressed in various tissues with the highest expression levels detected in labial palp and mantle, which were 86.83-fold (p < 0.05) and 75.87-fold (p < 0.05) higher than that in hepatopancreas. The mRNA expression of Cgengrailed-1 in mantle decreased dramatically after moderate (pH 7.8) and severe (pH 7.4) acidification treatment (0.75- and 0.15-fold of that in control group, p < 0.05). The results of immunofluorescence assay demonstrated that the expression level of Cgengrailed-1 in the middle fold of mantle increased significantly upon moderate and severe acidification treatment. Moreover, after the oyster larvae received acidification treatment at trochophore stage, the mRNA expression levels of Cgengrailed-1 increased significantly in D-shape larvae stages, which was 3.11- (pH 7.8) and 4.39-fold (pH 7.4) of that in control group (p < 0.05). The whole-mount immunofluorescence assay showed that Cgengrailed-1 was mainly expressed on the margin of shell gland, and the periostracum in trochophore, early D-shape larvae and D-shape larvae in both control and acidification treatment groups, and the intensity of positive signals in early D-shape larvae and D-shape larvae increased dramatically under acidification treatment. These results collectively suggested that the expression of Cgengrailed-1 could be triggered by CO₂-driven acidification treatment, which might contribute to induce the initial shell formation in oyster larvae and the formation of periostracum in adult oyster to adapt to the acidifying marine environment.

A six-amino-acid motif is a major determinant in functional evolution of HOX1 proteins

Gene duplication and divergence is a major driver in the emergence of evolutionary novelties. How variations in amino acid sequences lead to loss of ancestral activity and functional diversification of proteins is poorly understood. We used cross-species functional analysis of Drosophila Labial and its mouse HOX1 orthologs (HOXA1, HOXB1, and HOXD1) as a paradigm to address this issue. Mouse HOX1 proteins display low (30%) sequence similarity with Drosophila Labial. However, substituting endogenous Labial with the mouse proteins revealed that HOXA1 has retained essential ancestral functions of Labial, while HOXB1 and HOXD1 have diverged. Genome-wide analysis demonstrated similar DNA-binding patterns of HOXA1 and Labial in mouse cells, while HOXB1 binds to distinct targets. Compared with HOXB1, HOXA1 shows an enrichment in co-occupancy with PBX proteins on target sites and exists in the same complex with PBX on chromatin. Functional analysis of HOXA1-HOXB1 chimeric proteins uncovered a novel six-amino-acid C-terminal motif (CTM) flanking the homeodomain that serves as a major determinant of ancestral activity. In vitro DNA-binding experiments and structural prediction show that CTM provides an important domain for interaction of HOXA1 proteins with PBX. Our findings show that small changes outside of highly conserved DNA-binding regions can lead to profound changes in protein function.

Mapping the native interaction surfaces of PREP1 with PBX1 by cross-linking mass-spectrometry and mutagenesis

Both onco-suppressor PREP1 and the oncogene MEIS1 bind to PBX1. This interaction stabilizes the two proteins and allows their translocation into the nucleus and thus their transcriptional activity. Here, we have combined cross-linking mass-spectrometry and systematic mutagenesis to detail the binding geometry of the PBX1-PREP1 (and PBX1-MEIS1) complexes, under native in vivo conditions. The data confirm the existence of two distinct interaction sites within the PBC domain of PBX1 and unravel differences among the highly similar binding sites of MEIS1 and PREP1. The HR2 domain has a fundamental role in binding the PBC-B domain of PBX1 in both PREP1 and MEIS1. The HR1 domain of MEIS1, however, seem to play a less stringent role in PBX1 interaction with respect to that of PREP1. This difference is also reflected by the different binding affinity of the two proteins to PBX1. Although partial, this analysis provides for the first time some ideas on the tertiary structure of the complexes not available before. Moreover, the extensive mutagenic analysis of PREP1 identifies the role of individual hydrophobic HR1 and HR2 residues, both in vitro and in vivo.

Identification and characterization of cytotoxic amyloid-like regions in human Pbx-regulating protein-1

The ability of many proteins to fold into well-defined structures has been traditionally considered a prerequisite for fulfilling their functions. Protein folding is also regarded as a valuable loophole to escape uncontrolled and harmful aggregations. Here we show that the PBX-regulating protein-1 (PREP1), an important homeodomain transcription factor involved in cell growth and differentiation during embryogenesis, is endowed with an uncommon thermostability. Indeed, circular dichroism analyses indicate that it retains most of its secondary structure at very high temperatures. These findings have important implications for PREP1 functions since it is a stabilizing factor of its partner PBX1. Predictive analyses suggest that the observed PREP1 thermostability could be related to the presence of aggregation-prone regions. Interestingly, synthetic peptides corresponding to these regions exhibit a remarkable propensity to form toxic β-rich amyloid-like aggregates in physiological conditions. On this basis, we suggest that PREP1 stability is an effective way to prevent or limit the formation of harmful aggregates. Notably, one of these PREP1 fragments (residues 117-132) is able to reversibly switch from α-helical to β-rich states depending on the environmental conditions. The chameleon conformational behavior of this peptide makes it an ideal system to study this intriguing and widespread structural transition.

A zebrafish hox gene acts before gastrulation to specify the hemangioblast

Hox genes encode transcription factors that have been implicated in embryonic, adult and disease processes. The earliest developmental program known to be directed by Hox genes is the timing of ingression of presumptive axial mesoderm during gastrulation. We previously used morpholino (MO)-based knockdown to implicate the zebrafish hoxd4a gene in the specification of the hemangioblast, an event occurring at pre-gastrulation stages, well before the earliest known Hox gene function. The precise time at which hoxd4a function is required for this specification is not defined. We therefore fused the hoxd4a coding region to the human estrogen receptor (hER^T2 ). Following co-injection of anti-hoxd4a MO with mRNA encoding the Hoxd4a-ER^T2 fusion protein, hemangioblast specification was fully rescued when embryos were exposed to the estrogen analog 4-hydroxy-tamoxifen (4-OHT) at 4 hr post-fertilization (hpf), but only poorly at 6 hpf and not at all at 8 hpf, thereby defining a pre-gastrulation role for Hoxd4a, the earliest developmental function of a vertebrate Hox gene so far described. Both DNA binding and interaction with cofactor Pbx were further shown to be required for rescue of the morphant phenotype. Confirmation of the morphant phenotype was sought via the generation of hoxd4a null mutants using CRISPR/Cas9 technology. Null mutants of hoxd4a up to the third generation (F₃ ) failed to recapitulate the morphant phenotype, and were largely refractory to the effects of injected anti-hoxd4a MO suggesting the action of genetic compensation.

A Self-Activation Loop Maintains Meristematic Cell Fate for Branching

In plants and animals, self-renewing stem cell populations play fundamental roles in many developmental contexts. Plants differ from most animals in their retained ability to initiate new cycles of growth and development, which relies on the establishment and activity of branch meristems. In seed plants, branching is achieved by stem-cell-containing axillary meristems, which are initiated from a leaf axil meristematic cell population originally detached from the shoot apical meristem. It remains unclear how the meristematic cell fate is maintained. Here, we show that ARABIDOPSISTHALIANAHOMEOBOXGENE1 (ATH1) maintains the meristem marker gene SHOOT MERISTEMLESS (STM) expression in the leaf axil to enable meristematic cell fate maintenance. Furthermore, ATH1 protein interacts with STM protein to form a STM self-activation loop. Genetic and biochemical data suggest that ATH1 anchors STM to activate STM as well as other axillary meristem regulatory genes. This auto-regulation allows the STM locus to remain epigenetically active. Taken together, our findings provide a striking example of a self-activation loop that maintains the flexibility required for stem cell niche re-establishment during organogenesis.

PBX1 promotes development of natural killer cells by binding directly to the Nfil3 promoter

The transcription factor nuclear factor interleukin-3-regulated protein (NFIL3, also called E4BP4) is crucial for commitment of natural killer (NK) cells from common lymphoid progenitors (CLPs). However, the identity of the factor that can regulate NFIL3 directly during the NK-cell development is not known. Here, we reveal that pre-B-cell leukemia transcription factor 1 (PBX1) can upregulate the NFIL3 expression directly. We used conditional knockout mice in which PBX1 in hematopoietic cells was specifically absent. The number of NK-committed progenitor pre-NKP cells and rNKP cells was reduced significantly in the absence of PBX1, which was consistent with NFIL3 deficiency. Also, the NFIL3 expression in NK cells was decreased if PBX1 was absent. We demonstrated that PBX1 was bound directly to the promoter of Nfil3 and facilitated transcription. Upon knockout of the binding site of PBX1 in the Nfil3 promoter, mice showed fewer NK-precursor cells and NK cells, just like that observed in Nfil3 knockout mice. Furthermore, asparagine N286 in the homeodomain of PBX1 controlled the binding of PBX1 to the Nfil3 promoter. Collectively, these findings demonstrate that the transcription factor PBX1 promotes the early development of NK cells by upregulating the Nfil3 expression directly.

Intrinsic DNA Shape Accounts for Affinity Differences between Hox-Cofactor Binding Sites

Transcription factors bind to their binding sites over a wide range of affinities, yet how differences in affinity are encoded in DNA sequences is not well understood. Here, we report X-ray crystal structures of four heterodimers of the Hox protein AbdominalB bound with its cofactor Extradenticle to four target DNA molecules that differ in affinity by up to ∼20-fold. Remarkably, despite large differences in affinity, the overall structures are very similar in all four complexes. In contrast, the predicted shapes of the DNA binding sites (i.e., the intrinsic DNA shape) in the absence of bound protein are strikingly different from each other and correlate with affinity: binding sites that must change conformations upon protein binding have lower affinities than binding sites that have more optimal conformations prior to binding. Together, these observations suggest that intrinsic differences in DNA shape provide a robust mechanism for modulating affinity without affecting other protein-DNA interactions.

Crystal Structure of the Double Homeodomain of DUX4 in Complex with DNA

Double homeobox (DUX) transcription factors are unique to eutherian mammals. DUX4 regulates expression of repetitive elements during early embryogenesis, but misexpression of DUX4 causes facioscapulohumeral muscular dystrophy (FSHD) and translocations overexpressing the DUX4 double homeodomain cause B cell leukemia. Here, we report the crystal structure of the tandem homeodomains of DUX4 bound to DNA. The homeodomains bind DNA in a head-to-head fashion, with the linker making anchoring DNA minor-groove interactions and unique protein contacts. Remarkably, despite being tandem duplicates, the DUX4 homeodomains recognize different core sequences. This results from an arginine-to-glutamate mutation, unique to primates, causing alternative positioning of a key arginine side chain in the recognition helix. Mutational studies demonstrate that this primate-specific change is responsible for the divergence in sequence recognition that likely drove coevolution of embryonically regulated repeats in primates. Our work provides a framework for understanding the endogenous function of DUX4 and its role in FSHD and cancer.

Genome-wide identification and characterization of TALE superfamily genes in cotton reveals their functions in regulating secondary cell wall biosynthesis

BACKGROUND

Cotton fiber length and strength are both key traits of fiber quality, and fiber strength (FS) is tightly correlated with secondary cell wall (SCW) biosynthesis. The three-amino-acid-loop-extension (TALE) superclass homeoproteins are involved in regulating diverse biological processes in plants, and some TALE members has been identified to play a key role in regulating SCW formation. However, little is known about the functions of TALE members in cotton (Gossypium spp.).

RESULTS

In the present study, based on gene homology, 46, 47, 88 and 94 TALE superfamily genes were identified in G. arboreum, G. raimondii, G. barbadense and G. hirsutum, respectively. Phylogenetic and evolutionary analysis showed the evolutionary conservation of two cotton TALE families (including BEL1-like and KNOX families). Gene structure analysis also indicated the conservation of GhTALE members under selection. The analysis of promoter cis-elements and expression patterns suggested potential transcriptional regulation functions in fiber SCW biosynthesis and responses to some phytohormones for GhTALE proteins. Genome-wide analysis of colocalization of TALE transcription factors with SCW-related QTLs revealed that some BEL1-like genes and KNAT7 homologs may participate in the regulation of cotton fiber strength formation. Overexpression of GhKNAT7-A03 and GhBLH6-A13 significantly inhibited the synthesis of lignocellulose in interfascicular fibers of Arabidopsis. Yeast two-hybrid (Y2H) experiments showed extensive heteromeric interactions between GhKNAT7 homologs and some GhBEL1-like proteins. Yeast one-hybrid (Y1H) experiments identified the upstream GhMYB46 binding sites in the promoter region of GhTALE members and defined the downstream genes that can be directly bound and regulated by GhTALE heterodimers.

CONCLUSION

We comprehensively identified TALE superfamily genes in cotton. Some GhTALE members are predominantly expressed during the cotton fiber SCW thicking stage, and may genetically correlated with the formation of FS. Class II KNOX member GhKNAT7 can interact with some GhBEL1-like members to form the heterodimers to regulate the downstream targets, and this regulatory relationship is partially conserved with Arabidopsis. In summary, this study provides important clues for further elucidating the functions of TALE genes in regulating cotton growth and development, especially in the fiber SCW biosynthesis network, and it also contributes genetic resources to the improvement of cotton fiber quality.

HOX genes as transcriptional and epigenetic regulators during tumorigenesis and their value as therapeutic targets

Several HOX genes are aberrantly expressed in a wide range of cancers interfering with their development and resistance to treatment. This seems to be often caused by alterations in the methylation profiles of their promoters. The role of HOX gene products in cancer is highly 'tissue specific', relying ultimately on their ability to regulate oncogenes or tumor-suppressor genes, directly as transcriptional regulators or indirectly interfering with the levels of epigenetic regulators. Nowadays, different strategies have been tested the use of HOX genes as therapeutic targets for cancer diagnosis and treatment. Here, we trace the history of the research concerning the involvement of HOX genes in cancer, their connection with epigenetic regulation and their potential use as therapeutic targets.

MEIS transcription factors in development and disease

MEIS transcription factors are key regulators of embryonic development and cancer. Research on MEIS genes in the embryo and in stem cell systems has revealed novel and surprising mechanisms by which these proteins control gene expression. This Primer summarizes recent findings about MEIS protein activity and regulation in development, and discusses new insights into the role of MEIS genes in disease, focusing on the pathogenesis of solid cancers.

Hox genes: Downstream "effectors" of retinoic acid signaling in vertebrate embryogenesis

One of the major regulatory challenges of animal development is to precisely coordinate in space and time the formation, specification, and patterning of cells that underlie elaboration of the basic body plan. How does the vertebrate plan for the nervous and hematopoietic systems, heart, limbs, digestive, and reproductive organs derive from seemingly similar population of cells? These systems are initially established and patterned along the anteroposterior axis (AP) by opposing signaling gradients that lead to the activation of gene regulatory networks involved in axial specification, including the Hox genes. The retinoid signaling pathway is one of the key signaling gradients coupled to the establishment of axial patterning. The nested domains of Hox gene expression, which provide a combinatorial code for axial patterning, arise in part through a differential response to retinoic acid (RA) diffusing from anabolic centers established within the embryo during development. Hence, Hox genes are important direct effectors of retinoid signaling in embryogenesis. This review focuses on describing current knowledge on the complex mechanisms and regulatory processes, which govern the response of Hox genes to RA in several tissue contexts including the nervous system during vertebrate development.

Targeting HOX/PBX dimer formation as a potential therapeutic option in esophageal squamous cell carcinoma

Homeobox genes are known to be classic examples of the intimate relationship between embryogenesis and tumorigenesis, which are a family of transcriptional factors involved in determining cell identity during early development, and also dysregulated in many malignancies. Previously, HOXB7, HOXC6 and HOXC8 were found abnormally upregulated in esophageal squamous cell carcinoma (ESCC) tissues compared with normal mucosa and seen as poor prognostic predictors for ESCC patients, and were shown to promote cell proliferation and anti‐apoptosis in ESCC cells. These three HOX members have a high level of functional redundancy, making it difficult to target a single HOX gene. The aim of the present study was to explore whether ESCC cells are sensitive to HXR9 disrupting the interaction between multiple HOX proteins and their cofactor PBX, which is required for HOX functions. ESCC cell lines (KYSE70, KYSE150, KYSE450) were treated with HXR9 or CXR9, and coimmunoprecipitation and immunofluorescent colocalization were carried out to observe HOX/PBX dimer formation. To further investigate whether HXR9 disrupts the HOX pro‐oncogenic function, CCK‐8 assay and colony formation assay were carried out. Apoptosis was assessed by flow cytometry, and tumor growth in vivo was investigated in a xenograft model. RNA‐seq was used to study the transcriptome of HXR9‐treated cells. Results showed that HXR9 blocked HOX/PBX interaction, leading to subsequent transcription alteration of their potential target genes, which are involved in JAK‐signal transducer and activator of transcription (STAT) activation and apoptosis inducement. Meanwhile, HXR9 showed an antitumor phenotype, such as inhibiting cell proliferation, inducing cell apoptosis and significantly retarding tumor growth. Therefore, it is suggested that targeting HOX/PBX may be a novel effective treatment for ESCC.

The human HOXA9 protein uses paralog-specific residues of the homeodomain to interact with TALE-class cofactors

HOX proteins interact with PBX and MEIS cofactors, which belong to the TALE-class of homeodomain (HD)-containing transcription factors. Although the formation of HOX-PBX complexes depends on a unique conserved HOX motif called hexapeptide (HX), the additional presence of MEIS induces a remodeling of the interaction, leading to a global dispensability of the HX motif for trimeric complex formation in the large majority of HOX proteins. In addition, it was shown that the anterior HOXB3 and central HOXA7 and HOXC8 proteins could use different alternative TALE interaction motifs, with or without the HX motif, depending on the DNA-binding site and cell context. Here we dissected the molecular interaction properties of the human posterior HOXA9 protein with its TALE cofactors, PBX1 and MEIS1. Analysis was performed on different DNA-binding sites in vitro and by doing Bimolecular Fluorescence Complementation (BiFC) in different cell lines. Notably, we observed that the HOXA9-TALE interaction relies consistently on the redundant activity of the HX motif and two paralog-specific residues of the HOXA9 HD. Together with previous work, our results show that HOX proteins interact with their generic TALE cofactors through various modalities, ranging from unique and context-independent to versatile and context-dependent TALE binding interfaces.

The Enigmatic HOX Genes: Can We Crack Their Code?

Homeobox genes (HOX) are a large family of transcription factors that direct the formation of many body structures during early embryonic development. There are 39 genes in the subgroup of homeobox genes that constitute the human HOX gene family. Correct embryonic development of flies and vertebrates is, in part, mediated by the unique and highly regulated expression pattern of the HOX genes. Disruptions in these fine-tuned regulatory mechanisms can lead to developmental problems and to human diseases such as cancer. Unfortunately, the molecular mechanisms of action of the HOX family of transcription factors are severely under-studied, likely due to idiosyncratic details of their structure, expression, and function. We suggest that a concerted and collaborative effort to identify interacting protein partners, produce genome-wide binding profiles, and develop HOX network inhibitors in a variety of human cell types will lead to a deeper understanding of human development and disease. Within, we review the technological challenges and possible approaches needed to achieve this goal.

The Generic Facet of Hox Protein Function

Hox transcription factors are essential to promote morphological diversification of the animal body. A substantial number of studies have focused on how Hox proteins reach functional specificity, an issue that arises from the fact that these transcription factors control distinct developmental functions despite sharing similar molecular properties. In this review, we highlight that, besides specific functions, for which these transcription factors are renowned, Hox proteins also often have nonspecific functions. We next discuss some emerging principles of these generic functions and how they relate to specific functions and explore our current grasp of the underlying molecular mechanisms.

Characterization of the interaction interface and conformational dynamics of human TGIF1 homeodomain upon the binding of consensus DNA

The TG interacting factor-1 homeodomain (TGIF1-HD) binds with the consensus DNA motif 5'-TGTCA-3' in gene promoters through its three-amino acid loop extension (TALE) type homeodomain, and then recruits co-regulators to regulate gene expression. Although the solution NMR structure of human TGIF1-HD has been reported previously, little is known about its DNA binding mechanism. NMR titrations have been extensively used to study mechanisms of ligand binding to target proteins; however, an intermediate exchange occurred predominantly between TGIF1-HD in the free and bound states when titrated with the consensus DNA, which resulted in poor-quality NMR spectra and precluded further exploration of its interaction interface and conformational dynamics. Here, the helix α3 of TGIF1-HD was speculated as the specific DNA binding interface by hydrogen-deuterium exchange mass spectrometry (HDX-MS) experiments, and subsequently confirmed by chemical exchange saturation transfer (CEST) spectroscopy. In addition, simultaneous conformational changes in other regions, including α1 and α2, were induced by DNA binding, explaining the observation of chemical shift perturbations from extensive residues besides those located in α3. Further, low-populated DNA-bound TGIF1-HD undergoing a slow exchange at a rate of 130.2 ± 3.6 s^-1 was derived from the analysis of the CEST data, and two residues, R220 and R221, located in the middle of α3 were identified to be crucial for DNA binding. Our study provides structural and dynamic insights into the mechanisms of TGIF1-HD recognition of extensive promoter DNA.

Detecting and interpreting DNA methylation marks

The generation, alteration, recognition, and erasure of epigenetic modifications of DNA are fundamental to controlling gene expression in mammals. These covalent DNA modifications include cytosine methylation by AdoMet-dependent methyltransferases and 5-methylcytosine oxidation by Fe(II)-dependent and α-ketoglutarate-dependent dioxygenases. Sequence-specific transcription factors are responsible for interpreting the modification status of specific regions of chromatin. This review focuses on recent developments in characterizing the functional and structural links between the modification status of two DNA bases: 5-methylcytosine and 5-methyluracil (thymine).

Functions of TGIF homeodomain proteins and their roles in normal brain development and holoprosencephaly

Holoprosencephaly (HPE) is a frequent human forebrain developmental disorder with both genetic and environmental causes. Multiple loci have been associated with HPE in humans, and potential causative genes at 14 of these loci have been identified. Although TGIF1 (originally TGIF, for Thymine Guanine-Interacting Factor) is among the most frequently screened genes in HPE patients, an understanding of how mutations in this gene contribute to the pathogenesis of HPE has remained elusive. However, mouse models based on loss of function of Tgif1, and the related Tgif2 gene, have shed some light on how human TGIF1 variants might cause HPE. Functional analyses of TGIF proteins and of TGIF1 single nucleotide variants from HPE patients, combined with analysis of forebrain development in mouse embryos lacking both Tgif1 and Tgif2, suggest that TGIFs regulate the transforming growth factor ß/Nodal signaling pathway and sonic hedgehog (SHH) signaling independently. Although, some developmental processes that are regulated by TGIFs may be Nodal-dependent, it appears that the forebrain patterning defects and HPE in Tgif mutant mouse embryos is primarily due to altered signaling via the Shh pathway.

Atomistic molecular dynamics simulations of bioactive engrailed 1 interference peptides (EN1-iPeps)

The neural-specific transcription factor Engrailed 1 - is overexpressed in basal-like breast tumours. Synthetic interference peptides - comprising a cell-penetrating peptide/nuclear localisation sequence and the Engrailed 1-specific sequence from the N-terminus have been engineered to produce a strong apoptotic response in tumour cells overexpressing EN1, with no toxicity to normal or non Engrailed 1-expressing cells. Here scaled molecular dynamics simulations were used to study the conformational dynamics of these interference peptides in aqueous solution to characterise their structure and dynamics. Transitions from disordered to α-helical conformation, stabilised by hydrogen bonds and proline-aromatic interactions, were observed throughout the simulations. The backbone of the wild-type peptide folds to a similar conformation as that found in ternary complexes of anterior Hox proteins with conserved hexapeptide motifs important for recognition of pre-B-cell leukemia Homeobox 1, indicating that the motif may possess an intrinsic preference for helical structure. The predicted NMR chemical shifts of these peptides are consistent with the Hox hexapeptides in solution and Engrailed 2 NMR data. These findings highlight the importance of aromatic residues in determining the structure of Engrailed 1 interference peptides, shedding light on the rational design strategy of molecules that could be adopted to inhibit other transcription factors overexpressed in other cancer types, potentially including other transcription factor families that require highly conserved and cooperative protein-protein partnerships for biological activity.

A hydrophobic residue in the TALE homeodomain of PBX1 promotes epithelial-to-mesenchymal transition of gastric carcinoma

Pre-B-cell leukemia homeobox 1 (PBX1) was originally identified as a proto-oncogene in human leukemia. Although this protein has been shown to contribute to cellular development and tumorigenesis, the role of PBX1 in gastric carcinoma (GC) remains unclear. In this study, we observed increased expression of PBX1 in GC tissues compared with adjacent normal tissues. This increase in PBX1 expression levels negatively correlated with HOXB9 mRNA expression and was also associated with malignancy and metastasis. PBX1 promoted proliferation and metastasis of GC cells both in vitro and in vivo. These phenomena were also accompanied by epithelial-to-mesenchymal transition (EMT). Additionally, we observed that PBX1 promotes the expression of tumor growth and angiogenic factors. A structural model of the PBX1-HOX complex revealed that hydrophobic binding between PBX1 and the hexapeptide motif might be required for EMT induction. This analysis also demonstrated that the Phe²⁵² residue in the first helix of the TALE homeodomain is involved in the latter hydrophobic binding reaction. In vitro data from PBX1 mutants suggest that PBX1 cannot promote tumorigenesis of GC cells via EMT induction when Phe²⁵² residues lose hydrophobicity. It is likely that the presence of this residue is essential in facilitating hydrophobic binding with the hexapeptide motif. These findings suggest that PBX1 may be a potential target for GC treatment and this study provides a platform to elucidate the molecular mechanisms that underpin the role of PBX1 in GC tumorigenesis.

Two distinct DNA sequences recognized by transcription factors represent enthalpy and entropy optima

Most transcription factors (TFs) can bind to a population of sequences closely related to a single optimal site. However, some TFs can bind to two distinct sequences that represent two local optima in the Gibbs free energy of binding (ΔG). To determine the molecular mechanism behind this effect, we solved the structures of human HOXB13 and CDX2 bound to their two optimal DNA sequences, CAATAAA and TCGTAAA. Thermodynamic analyses by isothermal titration calorimetry revealed that both sites were bound with similar ΔG. However, the interaction with the CAA sequence was driven by change in enthalpy (ΔH), whereas the TCG site was bound with similar affinity due to smaller loss of entropy (ΔS). This thermodynamic mechanism that leads to at least two local optima likely affects many macromolecular interactions, as ΔG depends on two partially independent variables ΔH and ΔS according to the central equation of thermodynamics, ΔG = ΔH - TΔS.

Genes are sections of DNA that carry the instructions needed to build other molecules including all the proteins that the cell needs to fulfill its role. The information in the DNA is stored as a code consisting of four chemical bases, often referred to simply as “A”, “C”, “G” and “T”. The order or sequence of these bases determines the role of a protein. Many organisms – including humans – are built of many different types of cells that perform unique roles. Almost all cells carry the same genetic information, but proteins called transcription factors can regulate the activity of genes so that only a relevant subset of genes is switched on at a particular time.

Transcription factors glide along DNA and bind to short DNA sequences by attaching to the DNA bases directly or through bridges made up of water molecules. Two physical concepts known as enthalpy and entropy determine the strength of the connection. Enthalpy relates to how strong the chemical bonds that form between the transcription factors and the DNA bases are, compared to a situation where the transcription factor and DNA do not form a complex and bind to water molecules around them. Entropy measures the disorder of the system – the more disordered the solvent and protein-DNA complex are compared to solvent-containing free DNA and protein, the stronger the binding. A water molecule that bridges a DNA base with an amino-acid of a protein contributes to enthalpy, but results in loss of entropy, because the system becomes more ordered since the water molecule can no longer move freely.

Most transcription factors can only bind to DNA sequences that are very similar to each other, but some transcription factors can recognize several different kinds of sequences, and until now it was not clear how they could do this. Morgunova et al. studied four different human transcription factors that can each bind to two distinct DNA sequences. The results showed that the transcription factors bound to both DNA sequences with similar strength, but via different mechanisms. For one DNA sequence, an enthalpy-based mechanism essentially ‘froze’ the transcription factor to the DNA through rigid water bridges. The other DNA sequence was bound equally strongly but through moving water molecules, because this increased the entropy of the system. It is possible that these mechanisms could also apply to many other molecules that interact with each other through water-molecule bridges.

A better knowledge of the chemical bonds between transcription factors and DNA bases may in future help efforts to develop new treatments that depend on molecules being able to bind to other molecules. In addition, these findings may one day help scientists to predict how strongly two molecules will interact simply by knowing the structures of the molecules involved.

The O-GlcNAc transferase OGT interacts with and post-translationally modifies the transcription factor HOXA1

HOXA1 belongs to the HOX family of transcription factors which are key regulators of animal development. Little is known about the molecular pathways controlling HOXA1. Recent data from our group revealed distinct partner proteins interacting with HOXA1. Among them, OGT is an O-linked N-acetylglucosamine (O-GlcNAc) transferase modifying a variety of proteins involved in different cellular processes including transcription. Here, we confirm OGT as a HOXA1 interactor, we characterise which domains of HOXA1 and OGT are required for the interaction, and we provide evidence that OGT post-translationally modifies HOXA1. Mass spectrometry experiments indeed reveal that HOXA1 can be phosphorylated on the AGGTVGSPQYIHHSY peptide and that upon OGT expression, the phosphate adduct is replaced by an O-GlcNAc group.

Targeting HOX/PBX dimers in cancer

The HOX and PBX gene families encode transcription factors that have key roles in establishing the identity of cells and tissues in early development. Over the last 20 years it has become apparent that they are also dysregulated in a wide range of solid and haematological malignancies and have a predominantly pro-oncogenic function. A key mode of transcriptional regulation by HOX and PBX proteins is through their interaction as a heterodimer or larger complex that enhances their binding affinity and specificity for DNA, and there is growing evidence that this interaction is a potential therapeutic target in malignancies that include prostate, breast, renal, ovarian and lung cancer, melanoma, myeloma, and acute myeloid leukaemia. This review summarizes the roles of HOX and PBX genes in cancer and assesses the therapeutic potential of HOX/PBX dimer inhibition, including the availability of biomarkers for its application in precision medicine.

PRDM14, a putative histone methyl-transferase, interacts with and decreases the stability and activity of the HOXA1 transcription factor

Understanding how the activity of transcription factors like HOX proteins is regulated remains a widely open question. In a recent screen for proteins interacting with HOXA1, we identified a PRDM protein family member, PRDM14, which is known to be transiently co-expressed with HOXA1 in epiblast cells before their specification towards somatic versus germ cell fate. Here, we confirm PRDM14 is an interactor of HOXA1 and we identify the homeodomain of HOXA1 as well as the PR domain and Zinc fingers of PRDM14 to be required for the interaction. An 11-His repeat of HOXA1 previously highlighted to contribute to HOXA1-mediated protein-protein interactions is also involved. At a functional level, we provide evidence that HOXA1 displays an unexpectedly long half-life and demonstrate that PRDM14 can reduce the stability and affect the transcriptional activity of HOXA1.

Impact of cytosine methylation on DNA binding specificities of human transcription factors

The majority of CpG dinucleotides in the human genome are methylated at cytosine bases. However, active gene regulatory elements are generally hypomethylated relative to their flanking regions, and the binding of some transcription factors (TFs) is diminished by methylation of their target sequences. By analysis of 542 human TFs with methylation-sensitive SELEX (systematic evolution of ligands by exponential enrichment), we found that there are also many TFs that prefer CpG-methylated sequences. Most of these are in the extended homeodomain family. Structural analysis showed that homeodomain specificity for methylcytosine depends on direct hydrophobic interactions with the methylcytosine 5-methyl group. This study provides a systematic examination of the effect of an epigenetic DNA modification on human TF binding specificity and reveals that many developmentally important proteins display preference for mCpG-containing sequences.

Germ cell-intrinsic requirement for the homeodomain transcription factor PKnox1/Prep1 in adult spermatogenesis

PKnox1 (also known as Prep1) belongs to the TALE family of homeodomain transcription factors that are critical for regulating growth and differentiation during embryonic and postnatal development in vertebrates. We demonstrate here that PKnox1 is required for adult spermatogenesis in a germ cell-intrinsic manner. Tamoxifen-mediated PKnox1 loss in the adult testes, as well as its germ cell-specific ablation, causes testis hypotrophy with germ cell apoptosis and, as a consequence, compromised spermatogenesis. In PKnox1-deficient testes, spermatogenesis was arrested at the c-Kit⁺ spermatogonia stage, with a complete loss of the meiotic spermatocytes, and was accompanied by compromised differentiation of the c-Kit⁺ spermatogonia. Taken together, these results indicate that PKnox1 is a critical regulator of maintenance and subsequent differentiation of the c-Kit⁺ stage of spermatogonia in the adult testes.

HOXA1 and TALE proteins display cross-regulatory interactions and form a combinatorial binding code on HOXA1 targets

Hoxa1 has diverse functional roles in differentiation and development. We identify and characterize properties of regions bound by HOXA1 on a genome-wide basis in differentiating mouse ES cells. HOXA1-bound regions are enriched for clusters of consensus binding motifs for HOX, PBX, and MEIS, and many display co-occupancy of PBX and MEIS. PBX and MEIS are members of the TALE family and genome-wide analysis of multiple TALE members (PBX, MEIS, TGIF, PREP1, and PREP2) shows that nearly all HOXA1 targets display occupancy of one or more TALE members. The combinatorial binding patterns of TALE proteins define distinct classes of HOXA1 targets, which may create functional diversity. Transgenic reporter assays in zebrafish confirm enhancer activities for many HOXA1-bound regions and the importance of HOX-PBX and TGIF motifs for their regulation. Proteomic analyses show that HOXA1 physically interacts on chromatin with PBX, MEIS, and PREP family members, but not with TGIF, suggesting that TGIF may have an independent input into HOXA1-bound regions. Therefore, TALE proteins appear to represent a wide repertoire of HOX cofactors, which may coregulate enhancers through distinct mechanisms. We also discover extensive auto- and cross-regulatory interactions among the Hoxa1 and TALE genes, indicating that the specificity of HOXA1 during development may be regulated though a complex cross-regulatory network of HOXA1 and TALE proteins. This study provides new insight into a regulatory network involving combinatorial interactions between HOXA1 and TALE proteins.

Structural perspective of cooperative transcription factor binding

In prokaryotes, individual transcription factors (TFs) can recognize long DNA motifs that are alone sufficient to define the genes that they induce or repress. In contrast, in higher organisms that have larger genomes, TFs recognize sequences that are too short to define unique genomic positions. In addition, development of multicellular organisms requires molecular systems that are capable of executing combinatorial logical operations. Co-operative recognition of DNA by multiple TFs allows both definition of unique genomic positions in large genomes, and complex information processing at the level of individual regulatory elements. The TFs can co-operate in multiple different ways, and the precise mechanism used for co-operation determines important features of the regulatory interactions. Here, we present an overview of the structural basis of the different mechanisms by which TFs can cooperate, focusing on insight from recent functional studies and structural analyses of specific TF-TF-DNA complexes.

New Insights into Cooperative Binding of Homeodomain Transcription Factors PREP1 and PBX1 to DNA

PREP1 and PBX1 are homeodomain (HD) transcription factors that play crucial roles in embryonic development. Here, we present the first biophysical characterization of a PREP1 HD, and the NMR spectroscopic study of its DNA binding pocket. The data show that residues flanking the HD participate in DNA binding. The kinetic parameters for DNA binding of individual PREP1 and PBX1 HDs, and of their combination, show that isolated PREP1 and PBX1 HDs bind to DNA in a cooperative manner. A novel PREP1 motif, flanking the HD at the C-terminus, is required for cooperativity.

Hox functional diversity: Novel insights from flexible motif folding and plastic protein interaction

How the formidable diversity of forms emerges from developmental and evolutionary processes is one of the most fascinating questions in biology. The homeodomain-containing Hox proteins were recognized early on as major actors in diversifying animal body plans. The molecular mechanisms underlying how this transcription factor family controls a large array of context- and cell-specific biological functions is, however, still poorly understood. Clues to functional diversity have emerged from studies exploring how Hox protein activity is controlled through interactions with PBC class proteins, also evolutionary conserved HD-containing proteins. Recent structural data and molecular dynamic simulations add further mechanistic insights into Hox protein mode of action, suggesting that flexible folding of protein motifs allows for plastic protein interaction. As we discuss in this review, these findings define a novel type of Hox-PBC interaction, weak and dynamic instead of strong and static, hence providing novel clues to understanding Hox transcriptional specificity and diversity.