Homeodomain-DNA recognition

Genome-wide identification, evolutionary expansion, and expression profile of homeodomain-leucine zipper gene family in poplar (Populus trichocarpa)

Background

Homeodomain-leucine zipper (HD-ZIP) proteins are plant-specific transcriptional factors known to play crucial roles in plant development. Although sequence phylogeny analysis of Populus HD-ZIPs was carried out in a previous study, no systematic analysis incorporating genome organization, gene structure, and expression compendium has been conducted in model tree species Populus thus far.

Principal Findings

In this study, a comprehensive analysis of Populus HD-ZIP gene family was performed. Sixty-three full-length HD-ZIP genes were found in Populus genome. These Populus HD-ZIP genes were phylogenetically clustered into four distinct subfamilies (HD-ZIP I–IV) and predominately distributed across 17 linkage groups (LG). Fifty genes from 25 Populus paralogous pairs were located in the duplicated blocks of Populus genome and then preferentially retained during the sequential evolutionary courses. Genomic organization analyses indicated that purifying selection has played a pivotal role in the retention and maintenance of Populus HD-ZIP gene family. Microarray analysis has shown that 21 Populus paralogous pairs have been differentially expressed across different tissues and under various stresses, with five paralogous pairs showing nearly identical expression patterns, 13 paralogous pairs being partially redundant and three paralogous pairs diversifying significantly. Quantitative real-time RT-PCR (qRT-PCR) analysis performed on 16 selected Populus HD-ZIP genes in different tissues and under both drought and salinity stresses confirms their tissue-specific and stress-inducible expression patterns.

Conclusions

Genomic organizations indicated that segmental duplications contributed significantly to the expansion of Populus HD-ZIP gene family. Exon/intron organization and conserved motif composition of Populus HD-ZIPs are highly conservative in the same subfamily, suggesting the members in the same subfamilies may also have conservative functionalities. Microarray and qRT-PCR analyses showed that 89% (56 out of 63) of Populus HD-ZIPs were duplicate genes that might have been retained by substantial subfunctionalization. Taken together, these observations may lay the foundation for future functional analysis of Populus HD-ZIP genes to unravel their biological roles.

Cofactor binding evokes latent differences in DNA binding specificity between Hox proteins

Members of transcription factor families typically have similar DNA binding specificities yet execute unique functions in vivo. Transcription factors often bind DNA as multiprotein complexes, raising the possibility that complex formation might modify their DNA binding specificities. To test this hypothesis, we developed an experimental and computational platform, SELEX-seq, that can be used to determine the relative affinities to any DNA sequence for any transcription factor complex. Applying this method to all eight Drosophila Hox proteins, we show that they obtain novel recognition properties when they bind DNA with the dimeric cofactor Extradenticle-Homothorax (Exd). Exd-Hox specificities group into three main classes that obey Hox gene collinearity rules and DNA structure predictions suggest that anterior and posterior Hox proteins prefer DNA sequences with distinct minor groove topographies. Together, these data suggest that emergent DNA recognition properties revealed by interactions with cofactors contribute to transcription factor specificities in vivo.

Distinct DNA binding and transcriptional repression characteristics related to different ARX mutations

Mutations in the Aristaless-related homeobox gene (ARX) are associated with a wide variety of neurologic disorders including lissencephaly, hydrocephaly, West syndrome, Partington syndrome, and X-linked intellectual disability with or without epilepsy. A genotype-phenotype correlation exists for ARX mutations; however, the molecular basis for this association has not been investigated. To begin understanding the molecular basis for ARX mutations, we tested the DNA binding sequence preference and transcriptional repression activity for Arx, deletion mutants and mutants associated with various neurologic disorders. We found DNA binding preferences of Arx are influenced by the amino acid sequences adjacent to the homeodomain. Mutations in the homeodomain show a loss of DNA binding activity, while the T333N and P353R homeodomain mutants still possess DNA binding activities, although less than the wild type. Transcription repression activity, the primary function of ARX, is reduced in all mutants except the L343Q, which has no DNA binding activity and does not functionally repress Arx targets. These data indicate that mutations in the homeodomain result in not only a loss of DNA binding activity but also loss of transcriptional repression activity. Our results provide novel insights into the pathogenesis of ARX-related disorders and possible directions to pursue potential therapeutic interventions.

Learning generative models of molecular dynamics

We introduce three algorithms for learning generative models of molecular structures from molecular dynamics simulations. The first algorithm learns a Bayesian-optimal undirected probabilistic model over user-specified covariates (e.g., fluctuations, distances, angles, etc). L₁ reg-ularization is used to ensure sparse models and thus reduce the risk of over-fitting the data. The topology of the resulting model reveals important couplings between different parts of the protein, thus aiding in the analysis of molecular motions. The generative nature of the model makes it well-suited to making predictions about the global effects of local structural changes (e.g., the binding of an allosteric regulator). Additionally, the model can be used to sample new conformations. The second algorithm learns a time-varying graphical model where the topology and parameters change smoothly along the trajectory, revealing the conformational sub-states. The last algorithm learns a Markov Chain over undirected graphical models which can be used to study and simulate kinetics. We demonstrate our algorithms on multiple molecular dynamics trajectories.

Roles for intrinsic disorder and fuzziness in generating context-specific function in Ultrabithorax, a Hox transcription factor

Surprisingly few transcription factors drive animal development relative to the number and diversity of final tissues and body structures. Therefore, most transcription factors must function in more than one tissue. In a famous example, members of the Hox transcription factor family are expressed in contiguous stripes along the anterior/posterior axis during animal development. Individual Hox transcription factors specify all tissues within their expression domain and thus must respond to cellular cues to instigate the correct tissue-specific gene regulatory cascade. We describe how, in the Drosophila Hox protein Ultrabithorax, intrinsically disordered regions implement, regulate and co-ordinate multiple functions, potentially enabling context-specific gene regulation. The large N-terminal disordered domain encodes most of the transcription activation domain and directly impacts DNA binding affinity by the Ubx homeodomain. Similarly, the C-terminal disordered domain alters DNA binding affinity and specificity, interaction with a Hox binding protein and strongly influences both transcription activation and repression. Phosphorylation of the N-terminal disordered domain and alternative splicing of the C-terminal disordered domain could allow the cell to both regulate and co-ordinate DNA binding, protein interactions and transcription regulation. For regulatory mechanisms relying on disorder to continue to be available when Ubx is bound to other proteins or DNA, fuzziness would need to be preserved in these macromolecular complexes. The intrinsically disordered domains in Hox proteins are predicted to be on the very dynamic end of the disorder spectrum, potentially allowing disorder to persist when Ubx is bound to proteins or DNA to regulate the function of these "fuzzy" complexes. Because both intrinsically disordered regions within Ubx have multiple roles, each region may implement several different regulatory mechanisms identified in fuzzy complexes. These intrinsic disorder-based regulatory mechanisms are likely to be critical for allowing Ubx to sense tissue identity and respond by implementing a context-specific gene regulatory cascade.

Paracrine Pax6 activity regulates oligodendrocyte precursor cell migration in the chick embryonic neural tube

Homeoprotein transcription factors play fundamental roles in development, ranging from embryonic polarity to cell differentiation and migration. Research in recent years has underscored the physiological importance of homeoprotein intercellular transfer in eye field development, axon guidance and retino-tectal patterning, and visual cortex plasticity. Here, we have used the embryonic chick neural tube to investigate a possible role for homeoprotein Pax6 transfer in oligodendrocyte precursor cell (OPC) migration. We report the extracellular expression of Pax6 and the effects of gain and loss of extracellular Pax6 activity on OPCs. Open book cultures with recombinant Pax6 protein or Pax6 blocking antibodies, as well as in ovo gene transfer experiments involving expression of secreted Pax6 protein or secreted Pax6 antibodies, provide converging evidences that OPC migration is promoted by extracellular Pax6. The paracrine effect of Pax6 on OPC migration is thus a new example of direct non-cell autonomous homeoprotein activity.

Homeodomain transcription factors regulate BMP-2-induced osteoactivin transcription in osteoblasts

Osteoactivin (OA) is required for the differentiation of osteoblast cells. OA expression is stimulated by bone morphogenetic protein-2 (BMP-2). BMP-2 recruits homeodomain transcription factors Dlx3, Dlx5, and Msx2 to selectively activate or repress transcription of osteogenic genes and hence tightly regulate their transcription during osteoblast differentiation. Considering the key roles of Dlx3, Dlx5, and Msx2 in osteoblast differentiation, here we hypothesize that homeodomain proteins regulate BMP-2-induced OA transcription during osteoblast differentiation. Four classical homeodomain binding sites were identified in the proximal 0.96 kb region of rat OA promoter. Deletions and mutagenesis studies of the OA promoter region indicated that all four homeodomain binding sites are crucial for BMP-2-induced OA promoter activity. Simultaneous disruption of homeodomain binding sites at -852 and -843 of the transcription start site of OA gene significantly decreased the BMP-2-induced OA transcription and inhibited binding of Dlx3, Dlx5, and Msx2 proteins to the OA promoter. Dlx3 and Dlx5 proteins were found to activate the OA transcription, whereas, Msx2 suppressed BMP-2-induced OA transcription. Using chromatin immunoprecipitation assays, we demonstrated that the OA promoter is predominantly occupied by Dlx3 and Dlx5 during the proliferation and matrix maturation stages of osteoblast differentiation, respectively. During the matrix mineralization stage, BMP-2 robustly enhanced the recruitment of Dlx5 and to a lesser extent of Dlx3 and Msx2 to the OA promoter region. Collectively, our results show that the BMP-2-induced OA transcription is differentially regulated by Dlx3, Dlx5, and Msx2 during osteoblast differentiation.

Cooperative transcriptional activation by Klf4, Meis2, and Pbx1

The Kruppel-like factor Klf4 is implicated in tumorigenesis and maintaining stem cell pluripotency, and Klf4 can both activate and repress gene expression. We show that the Pbx1 and Meis2 homeodomain proteins interact with Klf4 and can be recruited to DNA elements comprising a Klf4 site or GC box, with adjacent Meis and Pbx sites. Meis2d and Pbx1a activate expression of p15(Ink4a) and E-cadherin, dependent on the Meis2d transcriptional activation domain. In HepG2 cells, reducing expression of endogenous Meis2 or Pbx1 decreases p15 gene expression and increases the number of cells entering S phase. Although DNA binding by all three proteins contributes to full cooperative activation, the sequence requirements for binding by Meis2 and Pbx1 are variable. In the E-cadherin promoter, a Pbx-like site is required for full activation, whereas in the p15 promoter, the Klf4 site appears to play the major role. Through a bioinformatics search we identified additional genes with conserved binding sites for Klf4, Meis2, and Pbx1 and show that at least some of these genes can be activated cooperatively by Klf4 and Meis2/Pbx1. We suggest a model in which genes with Klf4 sites can be cooperatively activated by Meis2/Pbx1 and Klf4, dependent primarily on recruitment by Klf4. This provides a mechanism to modulate transcriptional regulation by the multifunctional Klf4 transcription factor.

Use of a Drosophila genome-wide conserved sequence database to identify functionally related cis-regulatory enhancers

Background: Phylogenetic footprinting has revealed that cis-regulatory enhancers consist of conserved DNA sequence clusters (CSCs). Currently, there is no systematic approach for enhancer discovery and analysis that takes full-advantage of the sequence information within enhancer CSCs. Results: We have generated a Drosophila genome-wide database of conserved DNA consisting of >100,000 CSCs derived from EvoPrints spanning over 90% of the genome. cis-Decoder database search and alignment algorithms enable the discovery of functionally related enhancers. The program first identifies conserved repeat elements within an input enhancer and then searches the database for CSCs that score highly against the input CSC. Scoring is based on shared repeats as well as uniquely shared matches, and includes measures of the balance of shared elements, a diagnostic that has proven to be useful in predicting cis-regulatory function. To demonstrate the utility of these tools, a temporally-restricted CNS neuroblast enhancer was used to identify other functionally related enhancers and analyze their structural organization. Conclusions:cis-Decoder reveals that co-regulating enhancers consist of combinations of overlapping shared sequence elements, providing insights into the mode of integration of multiple regulating transcription factors. The database and accompanying algorithms should prove useful in the discovery and analysis of enhancers involved in any developmental process. Developmental Dynamics 241:169–189, 2012. © 2011 Wiley Periodicals, Inc.

How do Hox transcription factors find their target genes in the nucleus of living cells?

Homeotic mutations first found in Drosophila led to the identification of Hox genes in all bilateria. These genes are exceptional in that they are arranged in an ordered cluster, in which they are positioned in the same order along the chromosome as they are expressed along the antero-posterior axis to specify the corresponding body regions. They share a highly conserved DNA sequence of 180 bp, the homeobox which encodes the homeodomain, a 60 amino acid polypeptide involved in specific DNA and RNA binding and in protein-protein interactions. The discovery of the homeobox has uncovered for the first time a universal principle of specification of the body plan along the antero-posterior axis. The structure of the homeodomain has been determined by NMR spectroscopy and by X-ray crystallography. However, the mechanism by which the Hox proteins find their target genes in the nucleus of a living cell has been enigmatic. Transcriptome analysis indicates that there are hundreds of target genes to be regulated, both positively and negatively to ensure normal development. In the following, we show by Fluorescence Correlation Spectroscopy (FCS) and single molecule imaging in live salivary gland cells, that the mechanism of recognition is purely stochastic. The homeodomain associates and dissociates rapidly (in the ms range) with chromatin all along the chromosomes. If, however, it associates with a specific binding site in a puffed chromosome region, it remains bound for seconds or minutes to exert its function, by forming a complex with co-activators or co-repressors respectively. These direct measurements solve an old enigma of how Hox transcription factors find their target genes in the nucleus of live cells.

Good use and misuse of "genetic determinism"

After sequencing the human genome, scientists believed it would be possible to draw up a list of diseases, morphological characteristics and behavioral traits linked to each gene, but the post-genome era has shown that while links between genes and phenotypes, including behavioral phenotypes, do exist, they are more complex than was previously thought. There is no linear connection between genotype and brain and between brain and behavior; consequently, genomic and behavioral levels of organization are not isomorphous. There is no isomorphism because one gene plays many different roles, which means that the integrative processes needed for the development and functioning of an organism inevitably occurs in situations of non-linear multiple causality. Pleiotropy and epistasis, interactions between genes and the environment, alternative splicing and neuronal integration are all crucial mechanisms contributing to the many and varied aspects of brain-related genes.

Functional synthetic Antennapedia genes and the dual roles of YPWM motif and linker size in transcriptional activation and repression

Segmental identity along the anteroposterior axis of bilateral animals is specified by Hox genes. These genes encode transcription factors, harboring the conserved homeodomain and, generally, a YPWM motif, which binds Hox cofactors and increases Hox transcriptional specificity in vivo. Here we derive synthetic Drosophila Antennapedia genes, consisting only of the YPWM motif and homeodomain, and investigate their functional role throughout development. Synthetic peptides and full-length Antennapedia proteins cause head-to-thorax transformations in the embryo, as well as antenna-to-tarsus and eye-to-wing transformations in the adult, thus converting the entire head to a mesothorax. This conversion is achieved by repression of genes required for head and antennal development and ectopic activation of genes promoting thoracic and tarsal fates, respectively. Synthetic Antennapedia peptides bind DNA specifically and interact with Extradenticle and Bric-à-brac interacting protein 2 cofactors in vitro and ex vivo. Substitution of the YPWM motif by alanines abolishes Antennapedia homeotic function, whereas substitution of YPWM by the WRPW repressor motif, which binds the transcriptional corepressor Groucho, allows all proteins to act as repressors only. Finally, naturally occurring variations in the size of the linker between the homeodomain and YPWM motif enhance Antennapedia repressive or activating efficiency, emphasizing the importance of linker size, rather than sequence, for specificity. Our results clearly show that synthetic Antennapedia genes are functional in vivo and therefore provide powerful tools for synthetic biology. Moreover, the YPWM motif is necessary--whereas the entire N terminus of the protein is dispensable--for Antennapedia homeotic function, indicating its dual role in transcriptional activation and repression by recruiting either coactivators or corepressors.

Combinatorial regulation of photoreceptor differentiation factor, neural retina leucine zipper gene NRL, revealed by in vivo promoter analysis

Development and homeostasis require stringent spatiotemporal control of gene expression patterns that are established, to a large extent, by combinatorial action of transcription regulatory proteins. The bZIP transcription factor NRL (neural retina leucine zipper) is critical for rod versus cone photoreceptor cell fate choice during retinal development and acts as a molecular switch to produce rods from postmitotic precursors. Loss of Nrl in mouse leads to a cone-only retina, whereas ectopic expression of Nrl in photoreceptor precursors generates rods. To decipher the transcriptional regulatory mechanisms upstream of Nrl, we identified putative cis-control elements in the Nrl promoter/enhancer region by examining cross-species sequence conservation. Using in vivo transfection of promoter-reporter constructs into the mouse retina, we show that a 0.9-kb sequence upstream of the Nrl transcription initiation site is sufficient to drive reporter gene expression in photoreceptors. We further define a 0.3-kb sequence including a proximal promoter (cluster A1) and an enhancer (cluster B) that can direct rod-specific expression in vivo. Electrophoretic mobility shift assays using mouse retinal nuclear extracts, in combination with specific antibodies, demonstrate the binding of retinoid-related orphan nuclear receptor β (RORβ), cone rod homeobox, orthodenticle homolog 2, and cyclic AMP response element-binding protein to predicted consensus elements within clusters A and B. Our studies demonstrate Nrl as a direct transcriptional target of RORβ and suggest that combinatorial action of multiple regulatory factors modulates the expression of Nrl in developing and mature retina.

Modulation of plant growth by HD-Zip class I and II transcription factors in response to environmental stimuli

Plant development is adapted to changing environmental conditions for optimizing growth. This developmental adaptation is influenced by signals from the environment, which act as stimuli and may include submergence and fluctuations in water status, light conditions, nutrient status, temperature and the concentrations of toxic compounds. The homeodomain-leucine zipper (HD-Zip) I and HD-Zip II transcription factor networks regulate these plant growth adaptation responses through integration of developmental and environmental cues. Evidence is emerging that these transcription factors are integrated with phytohormone-regulated developmental networks, enabling environmental stimuli to influence the genetically preprogrammed developmental progression. Dependent on the prevailing conditions, adaptation of mature and nascent organs is controlled by HD-Zip I and HD-Zip II transcription factors through suppression or promotion of cell multiplication, differentiation and expansion to regulate targeted growth. In vitro assays have shown that, within family I or family II, homo- and/or heterodimerization between leucine zipper domains is a prerequisite for DNA binding. Further, both families bind similar 9-bp pseudopalindromic cis elements, CAATNATTG, under in vitro conditions. However, the mechanisms that regulate the transcriptional activity of HD-Zip I and HD-Zip II transcription factors in vivo are largely unknown. The in planta implications of these protein-protein associations and the similarities in cis element binding are not clear.

Differential Hematopoietic Supportive Potential and Gene Expression of Stroma Cell Lines from Midgestation Mouse Placenta and Adult Bone Marrow

During mouse embryogenesis, hematopoietic development takes place in several distinct anatomic locations. The microenvironment of different hematopoietic organs plays an important role in the proliferation and maturation of the hematopoietic cells. We hypothesized that fetal stromal cells would be distinct to adult bone marrow (BM)-derived stromal cells because the BM contributes mainly to the homeostasis of hematopoietic stem cells (HSCs), while extensive expansion of HSCs occurs during fetal development. Here we report the establishment of stromal cell lines from fetal hematopoietic organs, namely aorta-gonad-mesonephros (AGM), midgestation placenta (PL), and fetal liver (FL) together with adult bone marrow (BM). The growth patterns and hematopoietic supportive potential were studied. Their phenotypic and molecular gene expression profiles were also determined. Stromal cell lines from each tissue were able to support cobblestone area formation of BM c-Kit+Sca-1+ hematopoietic cells: 22 (22/47) from AGM, three (3/4) from PL, three (3/4) from FL, and three (3/3) from BM. There were similar levels of expansion of total mononuclear cells (TMNs) when HSCs were cocultured with fetal stroma and adult BM stroma. However, PL-derived stromal cells supported higher levels of generation of colony-forming progenitor cell (CFU-C), indicated by more colonies and colonies with significantly larger size. Flow cytometric analysis of the PL1 cells demonstrated a phenotype of CD45-, CD105+, Sca-1+, CD34+, and CD49d+, compared to adult BM1 cells, which were CD45-, CD105+, Sca-1+, CD34-, and CD49d-. Using Affymetrix microarray analysis, we identified that genes specifically express in endothelial cells, such as Tie1, Tek, Kdr, Flt4, Emcn, Pecam1, Icam2, Cdh5, Esam1, Prom1, Cd34, and Sele were highly expressed in stroma PL1, consistent with an endothelial phenotype, while BM1 expressed a mesenchymal stromal phenotype. In summary, these data demonstrate distinct characteristics of stromal cells that provide insights into the microenvironmental control of HSCs.

Genome-wide tissue-specific occupancy of the Hox protein Ultrabithorax and Hox cofactor Homothorax in Drosophila

The Hox genes are responsible for generating morphological diversity along the anterior-posterior axis during animal development. The Drosophila Hox gene Ultrabithorax (Ubx), for example, is required for specifying the identity of the third thoracic (T3) segment of the adult, which includes the dorsal haltere, an appendage required for flight, and the ventral T3 leg. Ubx mutants show homeotic transformations of the T3 leg towards the identity of the T2 leg and the haltere towards the wing. All Hox genes, including Ubx, encode homeodomain containing transcription factors, raising the question of what target genes Ubx regulates to generate these adult structures. To address this question, we carried out whole genome ChIP-chip studies to identify all of the Ubx bound regions in the haltere and T3 leg imaginal discs, which are the precursors to these adult structures. In addition, we used ChIP-chip to identify the sites bound by the Hox cofactor, Homothorax (Hth). In contrast to previous ChIP-chip studies carried out in Drosophila embryos, these binding studies reveal that there is a remarkable amount of tissue- and transcription factor-specific binding. Analyses of the putative target genes bound and regulated by these factors suggest that Ubx regulates many downstream transcription factors and developmental pathways in the haltere and T3 leg. Finally, we discovered additional DNA sequence motifs that in some cases are specific for individual data sets, arguing that Ubx and/or Hth work together with many regionally expressed transcription factors to execute their functions. Together, these data provide the first whole-genome analysis of the binding sites and target genes regulated by Ubx to specify the morphologies of the adult T3 segment of the fly.

Using the recognition code to swap homeodomain target specificity in cell culture

The homeodomain (HD) is a 60 amino acid-long DNA-binding domain. A large fraction of HDs binds with high affinity sequences containing the 5'-TAAT-3' core motif. However, NK-2 class HDs recognizes sequences containing the 5'-CAAG-3' core motif. By using a cell transfection approach, here we show that modification of residues located in the N-terminal arm (at positions 6, 7 and 8) and in the recognition helix (at position 54) is enough to swap the "in vivo" binding specificity of TTF-1 HD (which is a member of the NK-2 class HD) from 5'-CAAG-3' to 5'-TAAT-3'-containing targets. The role of residue at position 54 is also supported by data obtained with the HD of the Drosophila engrailed protein. These data support the notion that DNA-binding specificity "in vivo" is dictated by few critical residues.

The hemangioblast: from concept to authentication

The hemangioblast hypothesis has been hotly debated for over a century. Hemangioblasts are defined as multipotent cells that can give rise to both hematopoietic cells and endothelial cells. The existence of hemangioblasts has now been confirmed and many important molecules and several signaling pathways are involved in their generation and differentiation. Fibroblast growth factor, renin-angiotensin system and runt-related transcription factor 1 (Runx1) direct the formation of hemangioblasts through highly selective gene expression patterns. On the other hand, the hemogenic endothelium theory and a newly discovered pattern of hematopoietic/endothelial differentiation make the genesis of hemangioblasts more complicated. But how hemangioblasts are formed and how hematopoietic cells or endothelial cells are derived from remains largely unknown. Here we summarize the current knowledge of the signaling pathways and molecules involved in hemangioblast development and suggest some future clinical applications.

Functional analysis of conserved sequences within a temporally restricted neural precursor cell enhancer

Many of the key regulators of Drosophila CNS neural identity are expressed in defined temporal orders during neuroblast (NB) lineage development. To begin to understand the structural and functional complexity of enhancers that regulate ordered NB gene expression programs, we have undertaken the mutational analysis of the temporally restricted nerfin-1 NB enhancer. Our previous studies have localized the enhancer to a region just proximal to the nerfin-1 transcription start site. Analysis of this enhancer, using the phylogenetic footprint program EvoPrinter, reveals the presence of multiple sequence blocks that are conserved among drosophilids. cis-Decoder alignments of these conserved sequence blocks (CSBs) has identified shorter elements that are conserved in other Drosophila NB enhancers. Mutagenesis of the enhancer reveals that although each CSB is required for wild-type expression, neither position nor orientation of the CSBs within the enhancer is crucial for enhancer function; removal of less-conserved or non-conserved sequences flanking CSB clusters also does not significantly alter enhancer activity. While all three conserved E-box transcription factor (TF) binding sites (CAGCTG) are required for full function, adding an additional site at different locations within non-conserved sequences interferes with enhancer activity. Of particular note, none of the mutations resulted in ectopic reporter expression outside of the early NB expression window, suggesting that the temporally restricted pattern is defined by transcriptional activators and not by direct DNA binding repressors. Our work also points to an unexpectedly large number of TFs required for optimal enhancer function - mutant TF analysis has identified at least four that are required for full enhancer regulation.

Karyopherins in nuclear transport of homeodomain proteins during development

Homeodomain proteins are crucial transcription factors for cell differentiation, cell proliferation and organ development. Interestingly, their homeodomain signature structure is important for both their DNA-binding and their nucleocytoplasmic trafficking. The accurate nucleocytoplasmic distribution of these proteins is essential for their functions. We summarize information on (a) the roles of karyopherins for import and export of homeoproteins, (b) the regulation of their nuclear transport during development, and (c) the corresponding complexity of homeoprotein nucleocytoplasmic transport signals. This article is part of a Special Issue entitled: Regulation of Signaling and Cellular Fate through Modulation of Nuclear Protein Import.

Solution structures of the trihelix DNA-binding domains of the wild-type and a phosphomimetic mutant of Arabidopsis GT-1: mechanism for an increase in DNA-binding affinity through phosphorylation

GT-1 is a plant transcription factor that binds to one of the cis-acting elements, BoxII, which resides within the upstream promoter region of light-responsive genes. GT-1 was assumed to act as a molecular switch modulated through Ca(2+)-dependent phosphorylation/dephosphorylation in response to light signals. It was shown previously that the phosphorylation of threonine 133 in the DNA-binding domain (DBD) of GT-1 results in enhancement of the BoxII-binding activity. Interestingly, point mutation of Thr133 to Asp also enhances the BoxII-binding activity. Here, we report the solution structures of hypothetical trihelix DBDs of the wild-type (WT) and a phosphomimetic mutant (T133D) of GT-1. First, we demonstrated that the isolated DBD of GT-1 alone has the ability to bind to DNA, and that the T133D mutation of the isolated DBD can enhance the DNA-binding affinity. The structures of these DBDs turned out to be almost identical. The structural topology resembles that of Myb DBDs, but all α-helices are longer in GT-1. Our NMR titration experiments suggested that these longer α-helices yield an enlarged DNA-binding surface. The phosphorylation site is located at the N-terminus of the third α-helix. We built a structural model of the T133D DBD:BoxII complex with the program HADDOCK. The model resembles the structure of the TRF1 DBD:telomeric DNA complex. Interestingly, the model implies that the phosphorylated side chain may directly interact with the bases of DNA. On the basis of our findings, we propose a mechanism by which the DNA-binding activity toward BoxII of the phosphorylated GT-1 could be enhanced.

Suppression of colonic polyposis by homeoprotein CDX2 through its nontranscriptional function that stabilizes p27Kip1

Caudal-related homeoprotein CDX2 is expressed in intestinal epithelial cells, in which it is essential for their development and differentiation. A tumor suppressor function is suggested by evidence that CDX2 levels are decreased in human colon cancer specimens and that an inactivating mutation of Cdx2 in Apc(Δ716) mice markedly increases the incidence of colonic polyps. In this study, we investigated roles for transcriptional and nontranscriptional functions of CDX2 in suppression of colonic tumorigenesis. Mutagenic analysis of CDX2 revealed that loss of function stabilizes CDK inhibitor p27Kip1 by a nontranscriptional but homeodomain-dependent mechanism that inhibits cyclin E-CDK2 activity and blocks G0/G1-S progression in colon cancer cells. p27Kip1 stabilization was mediated by an inhibition of ubiquitylation-dependent proteolysis associated with decreased phosphorylation of Thr187 in p27Kip1. siRNA-mediated knockdown of p27Kip1 relieved the decrease in cyclin E-CDK2 activity and S-phase cell fraction elicited by CDX2 expression. Together, these results implicate a nontranscriptional function of CDX2 in tumor suppression mediated by p27Kip1 stabilization. Up to approximately 75% of low-CDX2 human colon cancer lesions show reduced levels of p27Kip1, whereas approximately 68% of high-CDX2 lesions retain expression of p27Kip1. These results show that low levels of CDX2 accelerate colon tumorigenesis by reducing p27Kip1 levels.

Penetratin story: an overview

Cell-penetrating peptides are short, often hydrophilic peptides that get access to the intracellular milieu. They have aroused great interest both in academic and applied research. First, cellular internalization of CPPs often involves the crossing of a biological membrane (plasma or vesicular), thus challenging the view of the nonpermeability of these structures to large hydrophilic molecules. Secondly, CPPs can drive the internalization of hydrophilic cargoes into cells, a rate-limiting step in the development of many therapeutic substances. Interestingly, the two mostly used CPPs, TAT and Penetratin peptides, are derived from natural proteins, HIV Tat and Antennapedia homeoprotein, respectively. The identification of the Penetratin peptide, summarized in this review, is intimately linked to the study of its parental natural protein.

Characterization of genomic structures and expression profiles of three tandem repeats of a mouse double homeobox gene: Duxbl

We identified and cloned a mouse double homeobox gene (Duxbl), which encodes two homeodomains. Duxbl gene, a tandem triplicate produces two major transcripts, Duxbl and Duxbl-s. The amino acid sequences of Duxbl homeodomains are most similar to those of human DUX4 protein, associated with facioscapulohumeral muscular dystrophy. In adult tissues, Duxbl is predominantly expressed in female reproductive organs and eyes, and slightly expressed in brain and testes. During gonad development, Duxbl is expressed from embryonic to adult stages and specifically expressed in oocytes and spermatogonia. During embryonic development, Duxbl is transcribed in limbs and tail. However, Duxbl proteins were only detected in trunk and limb muscles and in elongated myocytes and myotubes. In C2C12 muscle cell line, Duxbl expression pattern is similar to differentiated marker gene, Myogenin, increased in expression from 2 days onward in differentiating medium. We suggest that Duxbl proteins play regulatory roles during myogenesis and reproductive developments.

Comprehensive analysis of yeast surface displayed cDNA library selection outputs by exon microarray to identify novel protein-ligand interactions

Phosphatidylinositides are important signaling molecules that interact with a myriad of cellular proteins, many of which remain unidentified. We previously screened a yeast surface displayed human proteome library to identify protein fragments with affinity for the phosphatidylinositides, phosphatidylinositol-4,5-bisphosphate and phosphatidylinositol-3,4,5-trisphosphate. Much of the diversity in the screened selection outputs was represented by clones present at low frequencies, suggesting that a significant number of additional phosphatidylinositide-binding protein fragments might be present in the selection outputs. In the studies described in this report, we developed a novel cDNA library analysis method and comprehensively analyzed the polyclonal selection outputs from the phosphatidylinositol-4,5-bisphosphate and phosphatidylinositol-3,4,5-trisphosphate selections using a high-density exon microarray. In addition to the nine previously reported phosphatidylinositide-binding protein fragments, we identified 37 new phosphatidylinositide-binding candidates. Nine of 37 contain known phosphatidylinositide-binding domains, whereas the remaining 28 contain no known phosphatidylinositide-binding domain. We cloned and confirmed phosphatidylinositide binding by fluorescence-activated cell sorting for 17 of these novel candidate protein fragments. Our experiments suggest that phosphatidylinositide binding by these 17 novel protein fragments is dependent on both the inositol phosphate "headgroup" and the lipid "tail." This is in contrast with the PH domain containing fragments we tested, for which the inositol phosphate headgroup was sufficient for binding. The novel PtdIns-binding fragments come from a wide variety of proteins, including splicing factors, transcription factors, a kinase, and a polymerase. Intriguingly, 11 of the phosphatidylinositide-binding protein fragments are from nuclear proteins, including four containing homeobox domains. We found that phosphatidylinositides and double-stranded DNA oligonucleotides derived from homeobox domain target sequences compete for binding to homeobox domain-containing protein fragments, suggesting a possible mechanism for phospholipid-dependent transcriptional regulation. FACS enrichment of target-binding clones in yeast human cDNA display libraries coupled with comprehensive analysis of the selection output by DNA microarray analysis is an effective method for investigating common as well as rare protein interactions. In particular, this method is well suited for the study of small molecule/protein and drug/protein interactions.

Pax6 localizes to chromatin-rich territories and displays a slow nuclear mobility altered by disease mutations

The transcription factor Pax6 is crucial for the embryogenesis of multiple organs, including the eyes, parts of the brain and the pancreas. Mutations in one allele of PAX6 lead to eye diseases including Peter's anomaly and aniridia. Here, we use fluorescence recovery after photobleaching to show that Pax6 and also other Pax family proteins display a strikingly low nuclear mobility compared to other transcriptional regulators. For Pax6, the slow mobility is largely due to the presence of two DNA-binding domains, but protein-protein interactions also contribute. Consistently, the subnuclear localization of Pax6 suggests that it interacts preferentially with chromatin-rich territories. Some aniridia-causing missense mutations in Pax6 have impaired DNA-binding affinity. Interestingly, when these mutants were analyzed by FRAP, they displayed a pronounced increased mobility compared to wild-type Pax6. Hence, our results support the conclusion that disease mutations result in proteins with impaired function because of altered DNA- and protein-interaction capabilities.

Pentapeptide insertion mutagenesis of the Hoxa1 protein: mapping of transcription activation and DNA-binding regulatory domains

The mode of action of Hoxa1, like that of most Hox proteins, remains poorly characterized. In an effort to identify functional determinants contributing to the activity of Hoxa1 as a transcription factor, we generated 18 pentapeptide insertion mutants of the Hoxa1 protein and we assayed them in transfected cells for their activity on target enhancers from the EphA2 and Hoxb1 genes known to respond to Hoxa1 in the developing hindbrain. Only four mutants displayed a complete loss-of-function. Three of them contained an insertion in the homeodomain of Hoxa1, whereas the fourth loss-of-function mutant harbored an insertion in the very N-terminal end of the protein. Transcription activation assays in yeast further revealed that the integrity of both the N-terminal end and homeodomain is required for Hoxa1-mediated transcriptional activation. Furthermore, an insertion in the serine-threonine-proline rich C-terminal extremity of Hoxa1 induced an increase in activity in mammalian cells as well as in the yeast assay. The C-terminal extremity thus modulates the transcriptional activation capacity of the protein. Finally, electrophoretic mobility shift assays revealed that the N-terminal extremity of the protein also exerts a modulatory influence on DNA binding by Hoxa1-Pbx1a heterodimers.

Dissecting the functional specificities of two Hox proteins

Hox proteins frequently select and regulate their specific target genes with the help of cofactors like Extradenticle (Exd) and Homothorax (Hth). For the Drosophila Hox protein Sex combs reduced (Scr), Exd has been shown to position a normally unstructured portion of Scr so that two basic amino acid side chains can insert into the minor groove of an Scr-specific DNA-binding site. Here we provide evidence that another Drosophila Hox protein, Deformed (Dfd), uses a very similar mechanism to achieve specificity in vivo, thus generalizing this mechanism. Furthermore, we show that subtle differences in the way Dfd and Scr recognize their specific binding sites, in conjunction with non-DNA-binding domains, influence whether the target gene is transcriptionally activated or repressed. These results suggest that the interaction between these DNA-binding proteins and the DNA-binding site determines the architecture of the Hox-cofactor-DNA ternary complex, which in turn determines whether the complex recruits coactivators or corepressors.

An autoinhibitory effect of the homothorax domain of Meis2

Myeloid ecotropic insertion site (Meis)2 is a homeodomain protein containing a conserved homothorax (Hth) domain that is present in all Meis and Prep family proteins and in the Drosophila Hth protein. The Hth domain mediates interaction with Pbx homeodomain proteins, allowing for efficient DNA binding. Here we show that, like Meis1, Meis2 has a strong C-terminal transcriptional activation domain, which is required for full activation of transcription by homeodomain protein complexes composed of Meis2 and Pbx1. We also show that the activity of the activation domain is inhibited by the Hth domain, and that this autoinhibition can be partially relieved by the interaction of Pbx1 with the Hth domain of Meis2. Targeting of the Hth domain to DNA suggests that it is not a portable trans-acting repression domain. However, the Hth domain can inhibit a linked activation domain, and this inhibition is not limited to the Meis2 activation domain. Database searching reveals that the Meis3.2 splice variant, which is found in several vertebrate species, disrupts the Hth domain by removing 17 codons from the 5'-end of exon 6. We show that the equivalent deletion in Meis2 derepresses the C-terminal activation domain and weakens interaction with Pbx1. This work suggests that the transcriptional activity of all members of the Meis/Prep Hth protein family is subject to autoinhibition by their Hth domains, and that the Meis3.2 splice variant encodes a protein that bypasses this autoinhibitory effect.

Novel structural features in two ZHX homeodomains derived from a systematic study of single and multiple domains

BACKGROUND

Zhx1 to 3 (zinc-fingers and homeoboxes) form a set of paralogous genes encoding multi-domain proteins. ZHX proteins consist of two zinc fingers followed by five homeodomains. ZHXs have biological roles in cell cycle control by acting as co-repressors of the transcriptional regulator Nuclear Factor Y. As part of a structural genomics project we have expressed single and multi-domain fragments of the different human ZHX genes for use in structure determination.

RESULTS

A total of 30 single and multiple domain ZHX1-3 constructs selected from bioinformatics protocols were screened for soluble expression in E. coli using high throughput methodologies. Two homeodomains were crystallized leading to structures for ZHX1 HD4 and ZHX2 HD2. ZHX1 HD4, although closest matched to homeodomains from 'homez' and 'engrailed', showed structural differences, notably an additional C-terminal helix (helix V) which wrapped over helix I thereby making extensive contacts. Although ZHX2 HD2-3 was successfully expressed and purified, proteolysis occurred during crystallization yielding crystals of just HD2. The structure of ZHX2 HD2 showed an unusual open conformation with helix I undergoing 'domain-swapping' to form a homodimer.

CONCLUSIONS

Although multiple-domain constructs of ZHX1 selected by bioinformatics studies could be expressed solubly, only single homeodomains yielded crystals. The crystal structure of ZHX1 HD4 showed additional hydrophobic interactions relative to many known homeodomains via extensive contacts formed by the novel C-terminal helix V with, in particular, helix I. Additionally, the replacement of some charged covariant residues (which are commonly observed to form salt bridges in non-homeotherms such as the Drosophila 'engrailed' homeodomain), by apolar residues further increases hydrophobic contacts within ZHX1 HD4, and potentially stability, relative to engrailed homeodomain. ZHX1 HD4 helix V points away from the normally observed DNA major groove binding site on homeodomains and thus would not obstruct the putative binding of nucleic acid. In contrast, for ZHX2 HD2 the observed altered conformation involving rearrangement of helix I, relative to the canonical homeodomain fold, disrupts the normal DNA binding site, although protein-protein binding is possible as observed in homodimer formation.

PhenoFam-gene set enrichment analysis through protein structural information

BACKGROUND

With the current technological advances in high-throughput biology, the necessity to develop tools that help to analyse the massive amount of data being generated is evident. A powerful method of inspecting large-scale data sets is gene set enrichment analysis (GSEA) and investigation of protein structural features can guide determining the function of individual genes. However, a convenient tool that combines these two features to aid in high-throughput data analysis has not been developed yet. In order to fill this niche, we developed the user-friendly, web-based application, PhenoFam.

RESULTS

PhenoFam performs gene set enrichment analysis by employing structural and functional information on families of protein domains as annotation terms. Our tool is designed to analyse complete sets of results from quantitative high-throughput studies (gene expression microarrays, functional RNAi screens, etc.) without prior pre-filtering or hits-selection steps. PhenoFam utilizes Ensembl databases to link a list of user-provided identifiers with protein features from the InterPro database, and assesses whether results associated with individual domains differ significantly from the overall population. To demonstrate the utility of PhenoFam we analysed a genome-wide RNA interference screen and discovered a novel function of plexins containing the cytoplasmic RasGAP domain. Furthermore, a PhenoFam analysis of breast cancer gene expression profiles revealed a link between breast carcinoma and altered expression of PX domain containing proteins.

CONCLUSIONS

PhenoFam provides a user-friendly, easily accessible web interface to perform GSEA based on high-throughput data sets and structural-functional protein information, and therefore aids in functional annotation of genes.

Brassica GLABRA2 genes: analysis of function related to seed oil content and development of functional markers

Regulation of seed oil accumulation in oilseed rape (Brassica napus) has important economic significance. However, few genes have been characterized that affect final seed oil content. Through a mutant identification, the class IV homeodomain-ZIP transcription factor GLABRA2 (GL2) has been found to regulate seed oil accumulation in Arabidopsis, in addition to its role in trichome development. In this study, we isolated four distinct orthologues of GL2 from B. napus (AC-genome), B. rapa (A) and B. oleracea (C), using an overlapping-PCR strategy. The four GL2 orthologues were very similar, with 96.10-99.69% identity in exon regions, 75.45-93.84% in intron regions, 97.34-99.87% in amino acid sequences. Alignments of the four genes revealed that the A-genome sequences of BnaA.GL2.a from B. napus and BraA.GL2.a from B. rapa are more similar than the others, and likewise the C-genome sequences of BnaC.GL2.b from B. napus and BolC.GL2.a from B. oleracea are more similar. BnaA.GL2.a and BraA.GL2.a from the A-genome are highly expressed in roots, whilst BnaC.GL2.b and BolC.GL2.a from the C-genome are preferentially expressed in seeds. Transgenic ectopic overexpression and suppression of BnaC.GL2.b in Arabidopsis allowed further investigation of the effect on seed oil content. Overexpression generated two phenotypes: the wild-type-like and the gl2-mutant-like (an Arabidopsis glabrous mutant of gl2-2), with increases in seed oil content of 3.5-5.0% in the gl2-mutant-like transgenic plants. Suppression resulted in increases of 2.5-6.1% in seed oil content, and reduced trichome number at the leaf margins. These results suggest that BnaC.GL2.b can negatively regulate oil accumulation in Arabidopsis seeds. As a result of comparing the four GL2 genes, three A/C-genome-specific primer sets were developed and a C-genome-specific EcoRV cleavage site was identified, which can be used as functional markers to distinguish these orthologues within Brassica species. The genes identified and their molecular markers developed in this study will be valuable both for oilseed rape breeding focusing on improvement of seed oil content, and for detecting gene flow between populations.

DNA replication, development and cancer: a homeotic connection?

The homeotic proteins are transcription factors, highly conserved in metazoan organisms, exerting a pivotal role in development and differentiation. They individually display a loose specificity for the DNA sequence they can bind, but operate mainly in multi-molecular associations that assure their target and function specificity. Homeotic proteins are known to play a role in the positive or negative regulation of cell proliferation. Furthermore, many homeotic proteins are actually proto-oncogenes, since different translocations involving their genes cause tumors, particularly in the hematopoietic system. A one-hybrid screen to detect proteins with affinity for the lamin B2 replication origin identified three homeotic proteins, namely HoxA13, HoxC10 and HoxC13. Recent data demonstrate that the HoxC13 oncoprotein specifically associates with replication foci and binds in vitro and in vivo to several human DNA replication origins. Moreover, Hox proteins interact with geminin, a regulator of cell cycle progression, and control the interaction of this protein with the DNA replication licensing factor Ctd1. Thus, the homeotic proteins, by participating directly in the function of DNA replication origins, may provide a direct link between the accurate regulation of DNA replication required by the morphogenetic program and the deregulation of this process typical of cancer.

Cooperative DNA-binding and sequence-recognition mechanism of aristaless and clawless

To achieve accurate gene regulation, some homeodomain proteins bind cooperatively to DNA to increase those site specificities. We report a ternary complex structure containing two homeodomain proteins, aristaless (Al) and clawless (Cll), bound to DNA. Our results show that the extended conserved sequences of the Cll homeodomain are indispensable to cooperative DNA binding. In the Al-Cll-DNA complex structure, the residues in the extended regions are used not only for the intermolecular contacts between the two homeodomain proteins but also for the sequence-recognition mechanism of DNA by direct interactions. The residues in the extended N-terminal arm lie within the minor groove of DNA to form direct interactions with bases, whereas the extended conserved region of the C-terminus of the homeodomain interacts with Al to stabilize and localize the third alpha helix of the Cll homeodomain. This structure suggests a novel mode for the cooperativity of homeodomain proteins.

Autoregulatory binding sites in the zebrafish six3a promoter region define a new recognition sequence for Six3 proteins

The homeodomain (HD) transcription factor Six3, which is a member of the Six/Sine oculis family, is essential for development of the eyes and forebrain in vertebrates. It has recently been claimed that the HDs of Six3 and other members of the Six family have a common recognition sequence, TGATAC. However, a different recognition sequence including the typical TAAT core motif, which has not yet been fully defined, has also been proposed for the Six3 HD in mice. Our study of the zebrafish orthologue six3a, which has an identical HD, shows that it binds in vitro to multiple TAAT-containing sites within its promoter region. Comparison of the different binding affinities for these sequences identifies three high-affinity sites with a common TAATGTC motif. Notably, this new recognition sequence, which is supported by our analysis of the influence of single-nucleotide substitutions on the DNA-binding affinity, is distinct from all of the DNA-binding specificities previously described in surveys of HDs. In addition, our comparison of Six3a HD binding to the novel TAATGTC motif and the common recognition sequence of Six family HDs (TGATAC) shows very similar affinities, suggesting two distinct DNA-binding modes. Transient reporter assays of the six3a promoter in zebrafish embryos also indicate that the three high-affinity sites are involved in autoregulation. In support of this, chromatin immunoprecipitation experiments show enrichment of Six3a binding to a six3a promoter fragment containing two clustered high-affinity sites. These findings provide strong evidence that the TAATGTC motif is an important target sequence for vertebrate Six3 proteins in vivo.

Comparison of DNA binding across protein superfamilies

Specific protein-DNA interactions are central to a wide group of processes in the cell and have been studied both experimentally and computationally over the years. Despite the increasing collection of protein-DNA complexes, so far only a few studies have aimed at dissecting the structural characteristics of DNA binding among evolutionarily related proteins. Some questions that remain to be answered are: (a) what is the contribution of the different readout mechanisms in members of a given structural superfamily, (b) what is the degree of interface similarity among superfamily members and how this affects binding specificity, (c) how DNA-binding protein superfamilies distribute across taxa, and (d) is there a general or family-specific code for the recognition of DNA. We have recently developed a straightforward method to dissect the interface of protein-DNA complexes at the atomic level and here we apply it to study 175 proteins belonging to nine representative superfamilies. Our results indicate that evolutionarily unrelated DNA-binding domains broadly conserve specificity statistics, such as the ratio of indirect/direct readout and the frequency of atomic interactions, therefore supporting the existence of a set of recognition rules. It is also found that interface conservation follows trends that are superfamily-specific. Finally, this article identifies tendencies in the phylogenetic distribution of transcription factors, which might be related to the evolution of regulatory networks, and postulates that the modular nature of zinc finger proteins can explain its role in large genomes, as it allows for larger binding interfaces in a single protein molecule.

Quantitative study of synthetic Hox transcription factor-DNA interactions in live cells

Transcription factor-DNA interactions are life sustaining and therefore the subject of intensive research. In spite of vast effort, quantitative in vivo studies of the molecular mechanisms underlying these fundamental interactions remain challenging. In the preceding paper, we designed synthetic Sex combs reduced (Scr) peptides and validated genetically their function as transcriptional regulators. Here we present a controllable system for quantitative studies of protein-DNA interactions in live cells that enables us to "titrate" the concentration of the synthetic Scr peptides in a single cell. Using methods with single-molecule sensitivity, advanced fluorescence imaging and fluorescence correlation spectroscopy (FCS), we were able to study the kinetics of Scr-DNA interactions in live salivary gland cells, where Scr is normally expressed during development. We discerned freely moving Scr molecules, characterized the specific and nonspecific Scr peptide-DNA interactions, and estimated their corresponding dissociation constants (K(d)) in vivo. Our results suggest that the synthetic Scr transcription factors find their specific target sites primarily by multiple association/dissociation events, the rapidity of which is largely owed to electrostatic interactions. Based on these new findings, we formulate a model mechanism and emulate the kinetics of Scr homeodomain-DNA interactions in live cells using numerical simulations.

20ns molecular dynamics simulation of the antennapedia homeodomain-DNA complex: water interaction and DNA structure analysis

Homeodomains are one of the important families of eukaryotic DNA-binding motifs and provide an important model system for studying protein-DNA interactions. The crystal structure and NMR structure of the antennapedia homeodomain-DNA complex and comparison between them is available. Although earlier works have shown that the direct contacts and water mediated contacts are important for the binding and specificity. The detail dynamical structural characteristics of the complex, water mediating interactions in the complex and also the detail study of the free DNA and protein has not done. In the present paper we have reported the results of 20ns MD simulation of this complex with the presence of explicit water and also the 20ns MD simulation of the protein and the DNA separately in explicit water. The results show that the complex remains stable during the last 8ns of the simulation. The part of the protein which is interacting with the DNA has fewer fluctuations than other part of the protein. The pattern of water distribution around the interacting center has a typical pattern for this complex and it is quite different from the free protein and the free DNA. Water molecules penetrate into the interacting center during the simulation. Several water bridges have been identified which is responsible for recognition but not observed in the crystal structure. The recognized DNA sequence (14 mer) has been characterized by helical and step parameters. The correlated motions of the DNA and the protein in the complexed form and the free form has been analyzed.

The HMX/NKX homeodomain protein MLS-2 specifies the identity of the AWC sensory neuron type via regulation of the ceh-36 Otx gene in C. elegans

The differentiated features of postmitotic neurons are dictated by the expression of specific transcription factors. The mechanisms by which the precise spatiotemporal expression patterns of these factors are regulated are poorly understood. In C. elegans, the ceh-36 Otx homeobox gene is expressed in the AWC sensory neurons throughout postembryonic development, and regulates terminal differentiation of this neuronal subtype. Here, we show that the HMX/NKX homeodomain protein MLS-2 regulates ceh-36 expression specifically in the AWC neurons. Consequently, the AWC neurons fail to express neuron type-specific characteristics in mls-2 mutants. mls-2 is expressed transiently in postmitotic AWC neurons, and directly initiates ceh-36 expression. CEH-36 subsequently interacts with a distinct site in its cis-regulatory sequences to maintain its own expression, and also directly regulates the expression of AWC-specific terminal differentiation genes. We also show that MLS-2 acts in additional neuron types to regulate their development and differentiation. Our analysis describes a transcription factor cascade that defines the unique postmitotic characteristics of a sensory neuron subtype, and provides insights into the spatiotemporal regulatory mechanisms that generate functional diversity in the sensory nervous system.

Function and specificity of synthetic Hox transcription factors in vivo

Homeotic (Hox) genes encode transcription factors that confer segmental identity along the anteroposterior axis of the embryo. However the molecular mechanisms underlying Hox-mediated transcription and the differential requirements for specificity in the regulation of the vast number of Hox-target genes remain ill-defined. Here we show that synthetic Sex combs reduced (Scr) genes that encode the Scr C terminus containing the homedomain (HD) and YPWM motif (Scr-HD) are functional in vivo. Synthetic Scr-HD peptides can induce ectopic salivary glands in the embryo and homeotic transformations in the adult fly, act as transcriptional activators and repressors during development, and participate in protein-protein interactions. Their transformation capacity was found to be enhanced over their full-length counterpart and mutations known to transform the full-length protein into constitutively active or inactive variants behaved accordingly in the synthetic peptides. Our results show that synthetic Scr-HD genes are sufficient for homeotic function in Drosophila and suggest that the N terminus of Scr has a role in transcriptional potency, rather than specificity. We also demonstrate that synthetic peptides behave largely in a predictable way, by exhibiting Scr-specific phenotypes throughout development, which makes them an important tool for synthetic biology.

Cis-regulatory characterization of sequence conservation surrounding the Hox4 genes

Hox genes are key regulators of anterior-posterior axis patterning and have a major role in hindbrain development. The zebrafish Hox4 paralogs have strong overlapping activities in hindbrain rhombomeres 7 and 8, in the spinal cord and in the pharyngeal arches. With the aim to predict enhancers that act on the hoxa4a, hoxb4a, hoxc4a and hoxd4a genes, we used sequence conservation around the Hox4 genes to analyze all fish:human conserved non-coding sequences by reporter assays in stable zebrafish transgenesis. Thirty-four elements were functionally tested in GFP reporter gene constructs and more than 100 F1 lines were analyzed to establish a correlation between sequence conservation and cis-regulatory function, constituting a catalog of Hox4 CNEs. Sixteen tissue-specific enhancers could be identified. Multiple alignments of the CNEs revealed paralogous cis-regulatory sequences, however, the CNE sequence similarities were found not to correlate with tissue specificity. To identify ancestral enhancers that direct Hox4 gene activity, genome sequence alignments of mammals, teleosts, horn shark and the cephalochordate amphioxus, which is the most basal extant chordate possessing a single prototypical Hox cluster, were performed. Three elements were identified and two of them exhibited regulatory activity in transgenic zebrafish, however revealing no specificity. Our data show that the approach to identify cis-regulatory sequences by genome sequence alignments and subsequent testing in zebrafish transgenesis can be used to define enhancers within the Hox clusters and that these have significantly diverged in their function during evolution.

Regulation of Hox activity: insights from protein motifs

Deciphering the molecular bases of animal body plan construction is a central question in developmental and evolutionary biology. Genome analyses of a number of metazoans indicate that widely conserved regulatory molecules underlie the amazing diversity of animal body plans, suggesting that these molecules are reiteratively used for multiple purposes. Hox proteins constitute a good example of such molecules and provide the framework to address the mechanisms underlying transcriptional specificity and diversity in development and evolution. Here we examine the current knowledge of the molecular bases of Hox-mediated transcriptional control, focusing on how this control is encoded within protein sequences and structures. The survey suggests that the homeodomain is part of an extended multifunctional unit coordinating DNA binding and activity regulation and highlights the need for further advances in our understanding of Hox protein activity.