Understanding the homeodomain

HOPX: A Unique Homeodomain Protein in Development and Tumor Suppression

Simple Summary

Homeobox (HOX) genes encode homeodomain proteins that regulate a wide range of molecular pathways. The homeodomain is highly conserved and binds to DNA. One exception is homeodomain-only protein (HOPX) that lacks DNA-binding capacity. HOPX plays a crucial role in development and its functional impairment is associated with a variety of diseases, including cancer. Loss of HOPX function occurs in a wide range of cancer types, where it functions as a tumor suppressor gene. Understanding the molecular mechanisms by which HOPX regulates carcinogenesis will likely lead to the development of new therapeutic approaches.

Abstract

Homeobox genes are master regulators of morphogenesis and differentiation by acting at the top of genetic hierarchies and their deregulation is associated with a variety of human diseases. They usually contain a highly conserved sequence that codes for the homeodomain of the protein, a specialized motif with three α helices and an N-terminal arm that aids in DNA binding. However, one homeodomain protein, HOPX, is unique among its family members in that it lacks the capacity to bind DNA and instead functions by interacting with transcriptional regulators. HOPX plays crucial roles in organogenesis and is expressed in both embryonic and adult stem cells. Loss of HOPX expression is common in cancer, where it functions primarily as a tumor suppressor gene. In this review, we describe the function of HOPX in development and discuss its role in carcinogenesis.

Evolutionary Analysis of the Zinc Finger and Homeoboxes Family of Proteins Identifies Multiple Conserved Domains and a Common Early Chordate Ancestor

The Zinc Fingers and Homeoboxes (Zhx) proteins, Zhx1, Zhx2, and Zhx3, comprise a small family of proteins containing two amino-terminal C₂–H₂ zinc fingers and four or five carboxy-terminal homeodomains. These multiple homeodomains make Zhx proteins unusual because the majority of homeodomain-containing proteins contain a single homeodomain. Studies in cultured cells and mice suggest that Zhx proteins can function as positive or negative transcriptional regulators. Zhx2 regulates numerous hepatic genes, and all three Zhx proteins have been implicated in different cancers. Because Zhx proteins contain multiple predicted homeodomains, are associated with interesting physiological traits, and seem to be only present in the vertebrate lineage, we investigated the evolutionary history of this small family by comparing Zhx homologs from a wide range of chordates. This analysis indicates that the zinc finger motifs and homeodomains are highly similar among all Zhx proteins and also identifies additional Zhx-specific conserved regions, including a 13 amino acid amino-terminal motif that is nearly identical among all gnathostome Zhx proteins. We found single Zhx proteins in the sea lamprey (Petromyzon marinus) and in the nonvertebrate chordates sea squirt (Ciona intestinalis) and lancelet (Branchiostoma floridae); these Zhx proteins are most similar to gnathostome Zhx3. Based on our analyses, we propose that a duplication of the primordial Zhx gene gave rise to Zhx3 and the precursor to Zhx1 and Zhx2. A subsequent tandem duplication of this precursor generated Zhx1 and Zhx2 found in gnathostomes.

Upregulation of MGP by HOXC8 promotes the proliferation, migration, and EMT processes of triple-negative breast cancer

Triple-negative breast cancer (TNBC) is the most aggressive breast cancer subtype which accounts for 15%-20% of all breast cancer cases. The management of TNBC has remained a challenge due to its lack of targeted therapy. Previously, we reported that homeobox C8 (HOXC8) was involved in metastasis and migration of breast cancer cells. By chromatin immunoprecipitation and luciferase assays, we found that HOXC8 functioned as a transcription factor to activate the transcription of matrix Gla protein (MGP) gene, leading to an increase in the proliferation, anchorage-independent growth, and migration of TNBC cells. We further demonstrated that MGP expression promoted the epithelial-mesenchymal transition (EMT) process of TNBC cells, but not the other subtypes of breast cancer, suggesting that MGP induced EMT to promote proliferation and migration of TNBC cells. Moreover, we found that MGP was upregulated in clinical breast specimens compared to normal breast tissues and high MGP expression was statistically associated with poor, relapse-free survival for TNBC patients, indicating that MGP is probably a novel biomarker or therapeutic target for TNBC patients. Together, our results showed that the HOXC8-MGP axis played an important role in the tumorigenesis of TNBC and might be a promising therapeutic target for TNBC treatment.

Lysine Side-Chain Dynamics in the Binding Site of Homeodomain/DNA Complexes As Observed by NMR Relaxation Experiments and Molecular Dynamics Simulations

An important but poorly characterized contribution to the thermodynamics of protein-DNA interactions is the loss of entropy that occurs from restricting the conformational freedom of amino acid side chains. The effect of restricting the flexibility of several side chains at a protein-DNA interface may be comparable in many cases to the other factors that determine the binding thermodynamics and may, therefore, play a key role in dictating the binding affinity and/or specificity. Because the entropic contributions, including the presence and influence of side-chain dynamics, are especially difficult to estimate based on structural information, it is important to pursue experimental and theoretical studies that can provide direct information regarding these issues. We report on studies of a model system, the homeodomain/DNA complex, focusing on the Lys50 class of homeodomains where a key lysine residue in position 50 was shown previously to be critical for binding site specificity. NMR methodology was employed for determining the dynamics of lysine side-chain amino groups via 15N relaxation measurements in the Lys50-class homeodomains from the Drosophila protein Bicoid and the human protein Pitx2. In the case of Pitx2, complexes with both a consensus and a nonconsensus DNA binding site were examined. NMR-derived order parameters indicated moderate to substantial conformational freedom for the lysine NH3+ group in the complexes studied. To complement the experimental NMR measurements, molecular dynamics simulations were performed for the consensus complexes to gain further, detailed insights regarding the dynamics of the Lys50 side chain and other important residues in the protein-DNA interface.

Accurate prediction of functional, structural, and stability changes in PITX2 mutations using in silico bioinformatics algorithms

Mutations in PITX2 have been implicated in several genetic disorders, particularly Axenfeld-Rieger syndrome. In order to determine the most reliable bioinformatics tools to assess the likely pathogenicity of PITX2 variants, the results of bioinformatics predictions were compared to the impact of variants on PITX2 structure and function. The MutPred, Provean, and PMUT bioinformatic tools were found to have the highest performance in predicting the pathogenicity effects of all 18 characterized missense variants in PITX2, all with sensitivity and specificity >93%. Applying these three programs to assess the likely pathogenicity of 13 previously uncharacterized PITX2 missense variants predicted 12/13 variants as deleterious, except A30V which was predicted as benign variant for all programs. Molecular modeling of the PITX2 homoedomain predicts that of the 31 known PITX2 variants, L54Q, F58L, V83F, V83L, W86C, W86S, and R91P alter PITX2's structure. In contrast, the remaining 24 variants are not predicted to change PITX2's structure. The results of molecular modeling, performed on all the PITX2 missense mutations located in the homeodomain, were compared with the findings of eight protein stability programs. CUPSAT was found to be the most reliable in predicting the effect of missense mutations on PITX2 stability. Our results showed that for PITX2, and likely other members of this homeodomain transcription factor family, MutPred, Provean, PMUT, molecular modeling, and CUPSAT can reliably be used to predict PITX2 missense variants pathogenicity.

HOXB7 overexpression in lung cancer is a hallmark of acquired stem-like phenotype

HOXB7 is a homeodomain (HOX) transcription factor involved in regional body patterning of invertebrates and vertebrates. We previously identified HOXB7 within a ten-gene prognostic signature for lung adenocarcinoma, where increased expression of HOXB7 was associated with poor prognosis. This raises the question of how HOXB7 overexpression can influence the metastatic behavior of lung adenocarcinoma. Here, we analyzed publicly available microarray and RNA-seq lung cancer expression datasets and found that HOXB7-overexpressing tumors are enriched in gene signatures characterizing adult and embryonic stem cells (SC), and induced pluripotent stem cells (iPSC). Experimentally, we found that HOXB7 upregulates several canonical SC/iPSC markers and sustains the expansion of a subpopulation of cells with SC characteristics, through modulation of LIN28B, an emerging cancer gene and pluripotency factor, which we discovered to be a direct target of HOXB7. We validated this new circuit by showing that HOXB7 enhances reprogramming to iPSC with comparable efficiency to LIN28B or its target c-MYC, which is a canonical reprogramming factor.

Protein environment affects the water–tryptophan binding mode. MD, QM/MM, and NMR studies of engrailed homeodomain mutants

Water molecules can interact with the π-face of tryptophan either forming an O–H⋯π hydrogen bond or by a lone-pair⋯π interaction. Surrounding amino acids can favor the one or the other interaction type.

Structural effect of the L16Q, K50E, and R53P mutations on homeodomain of pituitary homeobox protein 2

The transcription factor pituitary homeobox protein 2 (PITX2) is involved in genetic control of development. Mutations in PITX2, most in the homeodomain, cause the autosomal-dominant disorder Rieger syndrome. The mutants L16Q, K50E and R53P destabilize the structure and disrupt DNA-binding activity. The biological functions of these mutants have been characterized but not the structural basis behind the loss of DNA-binding activity. We performed multiple molecular dynamics simulations at 37°C to investigate the structural and dynamic effects of the 3 PITX2 homeodomain mutants. Compared with the wild type (WT), the L16Q mutant induces a kink in the α3 helix, which is stabilized by the hydrogen bond of Q21-R59. The disruption in backbone hydrogen bonds of V47-N51 and W48-R52 leads to a kink formation in the α3 helix of K50E. The R53P mutant alters the relative orientation of helices, which is apparently stabilized by the formation of new hydrogen bonds of T38-Q11, T38-Q12, T38-R2, N39-R2, L40-Q1, L40-R2, and T41-Q4. The hydrophobic core residues F8, L13, L40 and V45 change their positions in all mutants to break the hydrophobic core. Thus, changes in helical orientations and hydrophobic core cause rearrangement of the DNA-binding surface and disrupt DNA-binding activity in the mutants. The structural and molecular dynamics properties of 3 PITX2 homeodomain mutants differ from those of the WT, especially in formation of a kink in the recognition helix, change in the packing of helices and disruption of the hydrophobic core. This structural basis for the loss of DNA-binding activity for these polymorphisms may help in understanding the effect of mutations on other homeodomains with other diseases.

Exploring the DNA-recognition potential of homeodomains

The recognition potential of most families of DNA-binding domains (DBDs) remains relatively unexplored. Homeodomains (HDs), like many other families of DBDs, display limited diversity in their preferred recognition sequences. To explore the recognition potential of HDs, we utilized a bacterial selection system to isolate HD variants, from a randomized library, that are compatible with each of the 64 possible 3' triplet sites (i.e., TAANNN). The majority of these selections yielded sets of HDs with overrepresented residues at specific recognition positions, implying the selection of specific binders. The DNA-binding specificity of 151 representative HD variants was subsequently characterized, identifying HDs that preferentially recognize 44 of these target sites. Many of these variants contain novel combinations of specificity determinants that are uncommon or absent in extant HDs. These novel determinants, when grafted into different HD backbones, produce a corresponding alteration in specificity. This information was used to create more explicit HD recognition models, which can inform the prediction of transcriptional regulatory networks for extant HDs or the engineering of HDs with novel DNA-recognition potential. The diversity of recovered HD recognition sequences raises important questions about the fitness barrier that restricts the evolution of alternate recognition modalities in natural systems.

Ceramide synthases in mammalians, worms, and insects: emerging schemes

The ceramide synthase (CerS) gene family comprises a group of highly conserved transmembrane proteins, which are found in all studied eukaryotes. The key feature of the CerS proteins is their role in ceramide synthase activity. Therefore, their original name ‘longevity assurance gene (Lass) homologs’, after the founding member, the yeast longevity assurance gene lag1, was altered to ‘CerS’. All CerS have high sequence similarity in a domain called LAG1 motif and a subset of CerS proteins is predicted to contain a Homeobox (Hox) domain. These domains could be the key to the multiple roles CerS have. CerS proteins play a role in diverse biological processes such as proliferation, differentiation, apoptosis, stress response, cancer, and neurodegeneration. In this review, we focus on CerS structure and biological function with emphasis of biological functions in the widely used model systems Caenorhabditis elegans and Drosophila melanogaster. Also, we focus on the accumulating data suggesting a role for CerS in lipid homeostasis.

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 Oct4 and Nanog transcription network regulates pluripotency in mouse embryonic stem cells

Oct4 and Nanog are transcription factors required to maintain the pluripotency and self-renewal of embryonic stem (ES) cells. Using the chromatin immunoprecipitation paired-end ditags method, we mapped the binding sites of these factors in the mouse ES cell genome. We identified 1,083 and 3,006 high-confidence binding sites for Oct4 and Nanog, respectively. Comparative location analyses indicated that Oct4 and Nanog overlap substantially in their targets, and they are bound to genes in different configurations. Using de novo motif discovery algorithms, we defined the cis-acting elements mediating their respective binding to genomic sites. By integrating RNA interference-mediated depletion of Oct4 and Nanog with microarray expression profiling, we demonstrated that these factors can activate or suppress transcription. We further showed that common core downstream targets are important to keep ES cells from differentiating. The emerging picture is one in which Oct4 and Nanog control a cascade of pathways that are intricately connected to govern pluripotency, self-renewal, genome surveillance and cell fate determination.

Homology modeling based solution structure of Hoxc8-DNA complex: role of context bases outside TAAT stretch

The 3D structure of neither Hoxc8 nor Hoxc8-DNA complex is known. The repressor protein Hoxc8 binds to the TAAT stretch of the promoter of the osteopontin gene and modulates its expression. Over expression of the osteopontin gene is related to diseases like osteoporosis, multiple sclerosis, cancer et cetera. In this paper we have proposed a 3D structure of Hoxc8-DNA complex obtained by Homology modeling and molecular dynamics (MD) simulation in explicit water. The crystal structure (9ant.pdb) of Antennapedia homeodomain in complex with its DNA sequence was chosen as the template based on (i) high sequence identity (85% for the protein and 60% for the DNA) and (ii) the presence of the TAAT stretch in interaction with the protein. The resulting model was refined by MD simulation for 2.0ns in explicit water. This refined model was then characterized in terms of the structural and the interactional features to improve our understanding of the mechanism of Hoxc8-DNA recognition. The interaction pattern shows that the residues Ile(195), Gln(198), and Asn(199), and the bases S2-(4)TAATG(8) are most important for recognition suggesting the stretch TAATG as the "true recognition element" in the present case. A strong and long-lived water bridge connecting Gln(198) and the base of S1-C(7) complementary to S2-G(8) was observed. Our predicted model of Hoxc8-DNA complex provides us with features that are consistent with the available experimental data on Hoxc8 and the general features of other homeodomain-DNA complexes. The predictions based on the model are also amenable to experimental verification.

Structural basis of Prospero-DNA interaction: implications for transcription regulation in developing cells

The crystal structure of a complex between the novel homeodomain of the neural transcription factor Prospero and DNA shows that the invariant residues Lys1290, Asn1294, and Asp1297 make specific contacts with the noncanonical DNA binding site. The overall structure includes the homeodomain and the adjacent Prospero domain and confirms that they act as a single structural unit, a Homeo-Prospero domain. The Prospero domain facilitates the proper alignment of the protein on the DNA. Knowledge of the structure reconciles two different DNA sequences that have been proposed as transcriptional targets for Prospero. As in the apo structure, the C terminus of the Prospero domain shields a short helix within the homeodomain that includes a nuclear export signal (NES). The structural results suggest that exposure of the NES is not coupled directly to DNA binding. We propose a DNA recognition mechanism specific to Prospero-type homeodomains in developing cells.

Solution structure of the K50 class homeodomain PITX2 bound to DNA and implications for mutations that cause Rieger syndrome

We have determined the solution structure of a complex containing the K50 class homeodomain Pituitary homeobox protein 2 (PITX2) bound to its consensus DNA site (TAATCC). Previous studies have suggested that residue 50 is an important determinant of differential DNA-binding specificity among homeodomains. Although structures of several homeodomain-DNA complexes have been determined, this is the first structure of a native K50 class homeodomain. The only K50 homeodomain structure determined previously is an X-ray crystal structure of an altered specificity mutant, Engrailed Q50K (EnQ50K). Analysis of the NMR structure of the PITX2 homeodomain indicates that the lysine at position 50 makes contacts with two guanines on the antisense strand of the DNA, adjacent to the TAAT core DNA sequence, consistent with the structure of EnQ50K. Our evidence suggests that this side chain may make fluctuating interactions with the DNA, which is complementary to the crystal data for EnQ50K. There are differences in the tertiary structure between the native K50 structure and that of EnQ50K, which may explain differences in affinity and specificity between these proteins. Mutations in the human PITX2 gene are responsible for Rieger syndrome, an autosomal dominant disorder. Analysis of the residues mutated in Rieger syndrome indicates that many of these residues are involved in DNA binding, while others are involved in formation of the hydrophobic core of the protein. Overall, the role of K50 in homeodomain recognition is further clarified, and the results indicate that native K50 homeodomains may exhibit differences from altered specificity mutants.

Homeobox gene expression profile in human hematopoietic multipotent stem cells and T-cell progenitors: implications for human T-cell development

Class I homeobox (HOX) genes comprise a large family of transcription factors that have been implicated in normal and malignant hematopoiesis. However, data on their expression or function during T-cell development is limited. Using degenerated RT-PCR and Affymetrix microarray analysis, we analyzed the expression pattern of this gene family in human multipotent stem cells from fetal liver (FL) and adult bone marrow (ABM), and in T-cell progenitors from child thymus. We show that FL and ABM stem cells are similar in terms of HOX gene expression, but significant differences were observed between these two cell types and child thymocytes. As the most immature thymocytes are derived from immigrated FL and ABM stem cells, this indicates a drastic change in HOX gene expression upon entry into the thymus. Further analysis of HOX-A7, HOX-A9, HOX-A10, and HOX-A11 expression with specific RT-PCR in all thymocyte differentiation stages showed a sequential loss of 3' region HOX-A cluster genes during intrathymic T-cell development and an unexpected expression of HOX-A11, previously not recognized to play a role in hematopoiesis. Also HOX-B3 and HOX-C4 were expressed throughout thymocyte development. Overall, these data provide novel evidence for an important role of certain HOX genes in human T-cell development.

Down-regulation of the myeloid homeobox protein Hex is essential for normal T-cell development

The haematopoietic homeobox gene Hex (also called Prh) is expressed in myeloid cells and B cells but not T cells. To investigate whether Hex levels might play a role in myeloid versus T-cell development, two types of transgenic mouse lines were constructed, each with ectopic expression of Hex in T cells (CD11a/Hex and Lck/Hex). Both these types of transgenic mouse had the same defects in T-cell maturation, indicating that proper T-cell development may be dependent not just on the up-regulation of lymphoid-specific transcriptional regulators but also on the co-ordinated down-regulation of myeloid-specific transcriptional regulators such as Hex. In addition, Hex over-expression significantly increased myeloid progenitor cycling, which may explain its role in retrovirally induced murine leukaemia.

New insights in the transcriptional activity and coregulator molecules in the arterial wall

A number of vascular diseases are associated with abnormal expression of genes that contribute to their pathophysiological and clinical manifestations, but at the same time offer potential therapeutic targets. One of the promising therapeutic approaches targets the pathophysiological pathways leading to aberrant gene activation, namely transcriptional activity and its molecular modulators (agonists, antagonists, coregulators, and nuclear receptors). The transcription factors can be divided into four classes (I-IV) classified by structural elements, like basic leucine zipper (bZIP) or basic helix-loop-helix (bHLH), which mediate their DNA binding activity but also determine the classes of drugs that can affect their activity. For example, statins modulate activation of the class-I transcription factor sterol responsive element-binding protein (SREBP), whose target genes including hydroxyl-methyl-glutaryl acetyl Coenzyme-A (HMG-CoA) reductase, HMG-CoA synthase, and the low-density lipoprotein receptor, all of which are involved in cholesterol and fatty acid metabolism. Similarly, insulin-like drugs target the nuclear receptor peroxisome-proliferator-activator-receptor (PPAR)-gamma (class-II), several anti-inflammatory drugs inhibit activation of nuclear factor kappa B (NFkappaB) (class-IV), while others (e.g. flavopiridol, rapamycin, and paclitaxel) target regulation of cell-cycle proteins. Increased understanding of the genetic and molecular basis of disease (e.g. transcriptional activity and its coregulation) will potentially enhance future diagnosis, treatment, and prevention of vascular diseases.

Modulation of Cardiac Growth and Development by HOP, an Unusual Homeodomain Protein

We have discovered an unusual homeodomain protein, called HOP, which is comprised simply of a homeodomain. HOP is highly expressed in the developing heart where its expression is dependent on the cardiac-restricted homeodomain protein Nkx2.5. HOP does not bind DNA and acts as an antagonist of serum response factor (SRF), which regulates the opposing processes of proliferation and myogenesis. Mice homozygous for a HOP null allele segregate into two phenotypic classes characterized by an excess or deficiency of cardiac myocytes. We propose that HOP modulates SRF activity during heart development; its absence results in an imbalance between cardiomyocyte proliferation and differentiation with consequent abnormalities in cardiac morphogenesis.

The binding of Ku antigen to homeodomain proteins promotes their phosphorylation by DNA-dependent protein kinase

The Ku antigen (70- and 80-kDa subunits) is a regulatory subunit of DNA-dependent protein kinase (DNA-PK) that promotes the recruitment of the catalytic subunit of DNA-PK (DNA-PKcs) to DNA ends and to specific DNA sequences from which the kinase is activated. Ku and DNA-PKcs plays essential roles in double-stranded DNA break repair and V(D)J recombination and have been implicated in the regulation of specific gene transcription. In a yeast two-hybrid screen of a Jurkat T cell cDNA library, we have identified a specific interaction between the 70-kDa subunit of Ku heterodimer and the homeodomain of HOXC4, a homeodomain protein expressed in the hematopoietic system. Unexpectedly, a similar interaction with Ku was observed for several additional homeodomain proteins including octamer transcription factors 1 and 2 and Dlx2, suggesting that specific binding to Ku may be a property shared by many homeodomain proteins. Ku-homeodomain binding was mediated through the extreme C terminus of Ku70 and was abrogated by amino acid substitutions at Lys595/Lys596. Ku binding allowed the recruitment of the homeodomain to DNA ends and dramatically enhanced the phosphorylation of homeodomain-containing proteins by DNA-PK. These results suggest that Ku functions as a substrate docking protein for signaling by DNA-PK to homeodomain proteins from DNA ends.

Identification of an optimal Ncx binding sequence required for transcriptional activation

The Ncx gene encodes a homeobox containing transcription factor that belongs to the Hox11 gene family. We determined specific Ncx protein binding consensus DNA sequences. Optimal Ncx binding sequences were 5'-CGGTAATTGG-3' (TAAT core) and 5'-CGGTAAGTGG-3' (TAAG core), which coincided with the Hox11 binding sequence. Both Ncx and Hox11 could bind to the TAAT and the TAAG core oligonucleotide in vitro. However, they could efficiently transactivate the reporter plasmid linked to the TAAT core sequence but not to the TAAG core sequence. Thus, Ncx and Hox11 act as transcriptional activators via their target sequence, 5'-CGGTAATTGG-3'.

Developmental patterning genes and their conserved functions: from model organisms to humans

Molecular and genetic evidence accumulated during the past 20 years in the field of developmental biology indicates that different animals possess many common genetic systems for embryonic patterning. In this review we describe the conserved functions of such developmental patterning genes and their relevance for human pathological conditions. Special attention is given to the Hox genetic system, involved in establishing cell identities along the anterior-posterior axis of all higher metazoans. We also describe other conserved genetic systems, such as the involvement of Pax6 genes in eye development and the role of Nkx2.5-type proteins in heart development. Finally, we outline some fascinating problems at the forefront of the studies of developmental patterning genes and show how knowledge obtained from model genetic organisms such as Drosophila helps to explain normal human morphogenesis and the genetic basis of some birth defects.

The leukemia-associated protein Btg1 and the p53-regulated protein Btg2 interact with the homeoprotein Hoxb9 and enhance its transcriptional activation

BTG1 and BTG2 belong to a family of functionally related genes involved in the control of the cell cycle. As part of an ongoing attempt to understand their biological functions, we used a yeast two-hybrid screening to look for possible functional partners of Btg1 and Btg2. Here we report the physical and functional association between these proteins and the homeodomain protein Hoxb9. We further show that Btg1 and Btg2 enhance Hoxb9-mediated transcription in transfected cells, and we report the formation of a Hoxb9.Btg2 complex on a Hoxb9-responsive target, and the fact that this interaction facilitates the binding of Hoxb9 to DNA. The transcriptional activity of the Hoxb9.Btg complex is essentially dependent on the activation domain of Hoxb9, located in the N-terminal portion of the protein. Our data indicate that Btg1 and Btg2 act as transcriptional cofactors of the Hoxb9 protein, and suggest that this interaction may mediate their antiproliferative function.

Structural study of homeodomain protein-DNA complexes using a homology modeling approach

The homeodomain is a conserved protein motif that binds to DNA and plays a central role in gene regulation. We use homeodomain as a model system to study the specific interactions between protein and DNA in a complex. Following the fundamental concept of homology modeling, we have developed an algorithm for predicting structures of both protein and DNA using the known structure of a similar complex as the template. The accuracies of the algorithm in predicting the complex structures are evaluated when two of the homeodomain protein-DNA complexes with known structures (antennapedia and MATalpha2) are selected as test systems. This algorithm allows structural studies of homeodomain binds to DNA with different sequences.

Activity regulation of Hox proteins, a mechanism for altering functional specificity in development and evolution

The closely related Hox transcription factors Ultrabithorax (Ubx) and Antennapedia (Antp) respectively direct first abdominal (A1) and second thoracic (T2) segment identities in Drosophila. It has been proposed that their functional differences derive from their differential occupancy of DNA target sites. Here we show that a hybrid version of Ubx (Ubx-VP16), which possesses an enhanced transcriptional activation function, no longer directs A1 denticle pattern in embryonic epidermal cells. Instead, it mimics Antp in directing T2 denticle pattern, and it can rescue the cuticular loss-of-function phenotype of Antp mutants. In cells that do not produce denticles, Ubx-VP16 appears to have largely retained its normal repressive regulatory functions. These results suggest that the modulation of Hox activation and repression functions can account for segment-specific morphological differences that are controlled by different members of the Hox family. Our results also are consistent with the idea that activity regulation underlies the phenotypic suppression phenomenon in which a more posterior Hox protein suppresses the function of a more anterior member of the Hox cluster. The acquisition of novel activation and repression potentials in Hox proteins may be an important mechanism underlying the generation of subtle morphological differences during evolution.

The expression of tobacco knotted1-type class 1 homeobox genes correspond to regions predicted by the cytohistological zonation model

We have isolated and characterized four tobacco homeobox genes, NTH1, NTH9, NTH20, NTH22 (Nicotiana tabacum homeobox) which belong to the class 1 knotted1-type family of homeobox genes. Comparison of the inferred amino acid sequences of the ELK homeodomains of these genes and previously reported kn1-type class 1 proteins has revealed that the four new tobacco genes belong to distinct subclasses, suggesting that each NTH gene may have distinct functions. Using in situ hybridization and by analysing the distribution of GUS activity in tobacco plants transformed with NTH promoter::GUS constructs, localized expression of the three NTH genes was observed in the shoot apical meristem (SAM). In the vegetative SAM, NTH1 and NTH15 showed overlapping expression in the corpus, NTH20 was expressed in the peripheral zone, and NTH9 was predominantly expressed in the rib zone. The expression patterns of the different NTH genes correspond to regions predicted by the cytohistological zonation model, suggesting that each NTH gene specifies the function of the SAM zone with which it is associated.

Transcriptional antagonism between Hmx1 and Nkx2.5 for a shared DNA-binding site

The recently described Hmx family of homeodomain proteins is predominately expressed in discrete regions of developing sensory tissues. In this report, we have identified the preferred DNA-binding site of the murine Hmx3 homeodomain protein by the selection and amplification binding (SAAB) technique. The consensus Hmx-binding site contained the sequence 5'-CAAGTG-3', which differs from the 5'-TAAT-3' motif commonly associated with homeodomain proteins. Instead, the Hmx consensus is similar to the 5'-CAAGTG-3'-binding sites of Nkx2.1 and Nkx2.5 homeodomain proteins. Based on mutation studies, both the 5'-CAAG-3' core and the 3'-TG dinucleotide are required for high affinity binding by Hmx3 and the homologous Hmx1 protein. A critical determinant of this specificity is the glutamine at position 50 in the third helix of the Hmx homeodomain. Hmx1 binds to the 5'-CAAGTG-3' element with an apparent dissociation constant of 20 nM. Unexpectedly, the human Hmx1 protein specifically repressed transcription from a luciferase reporter gene containing 3 copies of the 5'-CAAGTG-3' sequence. In contrast, the Nkx2.5 protein transactivated this luciferase reporter. Interestingly, co-expression of Hmx1 and Nkx2.5 attenuated each others activity, suggesting that genes containing the CAAGTG element can integrate signals from these proteins. Therefore, Hmx1 and Nkx2. 5 proteins bind a unique DNA sequence and act as transcriptional antagonists.

OC-2, a novel mammalian member of the ONECUT class of homeodomain transcription factors whose function in liver partially overlaps with that of hepatocyte nuclear factor-6

Transcription factors of the ONECUT class, whose prototype is hepatocyte nuclear factor (HNF)-6, are characterized by the presence of a single cut domain and by a peculiar homeodomain (Lannoy, V. J., Bürglin, T. R., Rousseau, G. G., and Lemaigre, F. P. (1998) J. Biol. Chem. 273, 13552-13562). We report here the identification and characterization of human OC-2, the second mammalian member of this class. The OC-2 gene is located on human chromosome 18. The distribution of OC-2 mRNA in humans is tissue-restricted, the strongest expression being detected in the liver and skin. The amino acid sequence of OC-2 contains several regions of high similarity to HNF-6. The recognition properties of OC-2 for binding sites present in regulatory regions of liver-expressed genes differ from, but overlap with, those of HNF-6. Like HNF-6, OC-2 stimulates transcription of the hnf-3beta gene in transient transfection experiments, suggesting that OC-2 participates in the network of transcription factors required for liver differentiation and metabolism.

A gamete-specific, sex-limited homeodomain protein in Chlamydomonas

During fertilization in Chlamydomonas, gametes of opposite mating types interact with each other through sex-specific adhesion molecules on their flagellar surfaces. Flagellar adhesion brings the cell bodies of the gametes into close contact and initiates a signal transduction pathway in preparation for cell-cell fusion. We have identified a cDNA, gsp1, whose transcript levels are upregulated during flagellar adhesion. The GSP1 polypeptide is a novel, gamete-specific homeodomain protein, the first to be identified in an alga. Its homeodomain shows significant identity with several higher plant homeodomain proteins. Although encoded by a single copy gene present in cells of both mating types, immunoblot analysis showed that GSP1 was expressed in mating type (mt)+ gametes, but was not detectable in mt- gametes or in vegetative cells of either mating type. Moreover, GSP1 appeared late during gametogenesis, suggesting that it may function during adhesion with mt- gametes or after zygote formation. GSP1 is expressed in imp11, mt- mutant gametes, which have a lesion in the mid gene involved in sex determination and exhibit many phenotypic characteristics of mt+ gametes. Thus, gsp1 is negatively regulated by mid and is the first molecule to be identified in Chlamydomonas that shows sex-limited expression.

Role of WUSCHEL in regulating stem cell fate in the Arabidopsis shoot meristem

The shoot meristem gives rise to the aerial parts of higher plants by continuously initiating new organs. The basis of this activity is its ability to maintain a pool of pluripotent stem cells, which are the ultimate source of all tissues of the shoot. In Arabidopsis plants mutant for the WUSCHEL (WUS) gene, the stem cells are misspecified and appear to undergo differentiation. Here, we show that WUS encodes a novel homeodomain protein which presumably acts as a transcriptional regulator. The pattern of WUS expression suggests that stem cells in the shoot meristem are specified by an underlying cell group which is established in the 16-cell embryo and becomes localized to its prospective domain of function by asymmetric cell divisions.

Engrailed homeodomain-DNA complex at 2.2 A resolution: a detailed view of the interface and comparison with other engrailed structures

We report the 2.2 A resolution structure of the Drosophila engrailed homeodomain bound to its optimal DNA site. The original 2.8 A resolution structure of this complex provided the first detailed three-dimensional view of how homeodomains recognize DNA, and has served as the basis for biochemical studies, structural studies and molecular modeling. Our refined structure confirms the principal conclusions of the original structure, but provides important new details about the recognition interface. Biochemical and NMR studies of other homeodomains had led to the notion that Gln50 was an especially important determinant of specificity. However, our refined structure shows that this side-chain makes no direct hydrogen bonds to the DNA. The structure does reveal an extensive network of ordered water molecules which mediate contacts to several bases and phosphates (including contacts from Gln50), and our model provides a basis for detailed comparison with the structure of an engrailed Q50K altered-specificity variant. Comparing our structure with the crystal structure of the free protein confirms that the N and C termini of the homeodomain become ordered upon DNA-binding. However, we also find that several key DNA contact residues in the recognition helix have the same conformation in the free and bound protein, and that several water molecules also are "preorganized" to contact the DNA. Our structure helps provide a more complete basis for the detailed analysis of homeodomain-DNA interactions.

The molecular basis of Rieger syndrome. Analysis of Pitx2 homeodomain protein activities

Rieger syndrome is an autosomal-dominant developmental disorder that includes glaucoma and mild craniofacial dysmorphism in humans. Mutations in the Pitx2 homeobox gene have been linked to Rieger syndrome. We have characterized wild type and mutant Pitx2 activities using electrophoretic mobility shift assays, protein binding, and transient transfection assays. Pitx2 preferentially binds the bicoid homeodomain binding site and transactivates reporter genes containing this site. The combination of Pitx2 and another homeodomain protein, Pit-1, yielded a synergistic 55-fold activation of the prolactin promoter in transfection assays. Addition of Pit-1 increased Pitx2 binding to the bicoid element in electrophoretic mobility shift assays. Furthermore, we demonstrate specific binding of Pit-1 to Pitx2 in vitro. Thus, wild type Pitx2 DNA binding activity is modulated by protein-protein interactions. We next studied two Rieger mutants. A threonine to proline mutation (T68P) in the second helix of the homeodomain retained DNA binding activity with the same apparent KD and only about a 2-fold reduction in the Bmax. However, this mutant did not transactivate reporter genes containing the bicoid site. The mutant Pitx2 protein binds Pit-1, but there was no detectable synergism on the prolactin promoter. A second mutation (L54Q) in a highly conserved residue in helix 1 of the homeodomain yielded an unstable protein. Our results provide insights into the potential mechanisms underlying the developmental defects in Rieger syndrome.

Transcription factors: an overview

This special issue of the International Journal of Biochemistry and Cell Biology contains a series of review articles and original papers dealing with the topic of transcription factors. The purpose of this introductory article is to provide an overview of these factors, their mechanism of action, their regulation and the manner in which alterations in them can result in disease.

Evidence that the DNA binding specificity of winged helix proteins is mediated by a structural change in the amino acid sequence adjacent to the principal DNA binding helix

We present the first structural evidence supporting the hypothesis that the binding specificity of the winged helix DNA binding motif is mediated by residues adjacent to the alpha-helix (H3), the moiety which is primarily involved in the interaction with DNA. Using NMR to determine secondary structural elements of a winged helix family member, Genesis (formerly HFH-2), and comparing these with those found in the X-ray crystal structure of the HNF-3gamma/DNA complex [Clark, K. L., Halay, E. D., Lai, E., & Burley, S. K. (1993) Nature 364, 412-420], we show that the major differences observed occur for H3 and the region immediately prior to this. H3 in Genesis is slightly shorter than in HNF-3gamma and, in addition, we observe an extra small helix (H4) in the region between H2 and H3 which is not found in the HNF-3gamma/DNA complex. This is significant as it has been shown previously [Overdier, D. G., Porcella, A., & Costa R. H. (1994) Mol. Cell. Biol. 14, 2755-2766] that the DNA-binding specificity is influenced by amino acid residues in this region.

Homology modeling using simulated annealing of restrained molecular dynamics and conformational search calculations with CONGEN: application in predicting the three-dimensional structure of murine homeodomain Msx-1

We have developed an automatic approach for homology modeling using restrained molecular dynamics and simulated annealing procedures, together with conformational search algorithms available in the molecular mechanics program CONGEN (Bruccoleri RE, Karplus M, 1987, Biopolymers 26:137-168). The accuracy of the method is validated by "predicting" structures of two homeodomain proteins with known three-dimensional structures, and then applied to predict the three-dimensional structure of the homeodomain of the murine Msx-1 transcription factor. Regions of the unknown protein structure that are highly homologous to the known template structure are constrained by "homology distance constraints," whereas the conformations of nonhomologous regions of the unknown protein are defined only by the potential energy function. A full energy function (excluding explicit solvent) is employed to ensure that the calculated structures have good conformational energies and are physically reasonable. As in NMR structure determinations, information on the consistency of the structure prediction is obtained by superposition of the resulting family of protein structures. In this paper, our homology modeling algorithm is described and compared with related homology modeling methods using spatial constraints derived from the structures of homologous proteins. The software is then used to predict the DNA-bound structures of three homeodomain proteins from the X-ray crystal structure of the engrailed homeodomain protein (Kissinger CR et al., 1990, Cell 63:579-590). The resulting backbone and side-chain conformations of the modeled yeast Mat alpha 2 and D. melanogaster Antennapedia homeodomains are excellent matches to the corresponding published X-ray crystal (Wolberger C et al., 1991, Cell 67:517-528) and NMR (Billeter M et al., 1993, J Mol Biol 234:1084-1097) structures, respectively. Examination of these structures of Msx-1 reveals a network of highly conserved surface salt bridges that are proposed to play a role in regulating protein-protein interactions of homeodomains in transcription complexes.

The yeast homeodomain protein MATalpha2 shows extended DNA binding specificity in complex with Mcm1

The MATalpha2 (alpha2) repressor interacts with the Mcm1 protein to turn off a-cell type-specific genes in the yeast Saccharomyces cerevisiae. We compared five natural alpha2-Mcm1 sites with an alpha2-Mcm1 symmetric consensus site (AMSC) for their relative strength of repression and found that the AMSC functions slightly better than any of the natural sites. To further investigate the DNA binding specificity of alpha2 in complex with Mcm1, symmetric substitutions at each position in the alpha2 half-sites of AMSC were constructed and assayed for their effect on repression in vivo and DNA binding affinity in vitro. As expected, substitutions at positions in which there are base-specific contacts decrease the level of repression. Interestingly, substitutions at other positions, in which there are no apparent base-specific contacts made by the protein in the alpha2-DNA co-crystal structure, also significantly decrease repression. As an alternative method to examining the DNA binding specificity of alpha2, we performed in vitro alpha2 binding site selection experiments in the presence and absence of Mcm1. In the presence of Mcm1, the consensus sequences obtained were extended and more closely related to the natural alpha2 sites than the consensus sequence obtained in the absence of Mcm1. These results demonstrate that in the presence of Mcm1 the sequence specificity of alpha2 is extended to these positions.

Site-specific recognition by an isolated DNA-binding domain of the sine oculis protein

The sine oculis (so) gene is required for the development of the Drosophila visual system. The 416 amino acid SO protein contains a 40 amino acid region homologous to the helix-turn-helix (HtH) region of the homeodomain. Three HtH-containing peptides ranging in size from 63 to 93 amino acids (SO(218-279), SO(204-279), and SO(188-279)) were expressed in Escherichia coli and characterized in vitro. These fragments show circular dichroism spectra characteristic of helical proteins and cooperative unfolding transitions. Derivatization of these three peptides with the chemical nuclease 1,10-phenanthroline:copper (OP-Cu) allowed the identification of specific DNA-binding sites within the 3.1 kb pUC119 plasmid. Similar cleavage patterns with similar relative affinities were obtained for all three peptides. Nucleotide resolution mapping of the predominant cleavage area identified two primary cleavage sites with a similar core sequence. The DNA cleavage sites were confirmed by DNase I footprinting with both native and OP-Cu-conjugated SO HtH peptides. This study identifies a 63 amino acid peptide as sufficient for specific DNA binding.

Three knotted1-like homeobox genes in Arabidopsis

Five arabidopsis kn1-like homeobox genes were cloned through low-stringency screening of Arabidopsis cDNA libraries with the kn1 homeobox from maize. These five genes were named KNAT1-5 (for kn1-like Arabidopsis thaliana). An analysis of KNAT1 and 2 has been presented previously [19]. Here we present an analysis of the genes KNAT3, 4 and 5. On the basis of sequence and expression patterns, these three genes belong to the class II subfamily of kn1-like homeobox genes [16]. Low-stringency Southern analysis suggests several additional members of the class II genes exist in the Arabidopsis genome. The predicted amino acid sequences of the three genes share extensive homology outside of the homeodomain, including 84% between KNAT3 and KNAT4. Northern analysis shows that although all three genes are expressed in all tissues examined, the level of KNAT3 RNA is highest in young siliques, inflorescences and roots, KNAT4 RNA level is strongest in leaves and young siliques, and KNAT5 RNA level is highest in roots. The specificity of these patterns was confirmed by RNA fingerprint analysis. KNAT3 and 4 are light-regulated as they show reduced expression in etiolated seedlings and also in hy3, cop1 and det1 mutant backgrounds.

Isolation of an Arabidopsis homologue of the maize homeobox Knotted-1 gene

The shoot apical meristem (SAM) is responsible for forming most of the above-ground portion of the plant. We sought to isolate regulatory genes expressed in the Arabidopsis SAM by screening a Brassica oleracea (cauliflower) meristem cDNA library with the homeobox fragment from the maize Knotted-1 (Kn1) gene. We isolated and characterized the corresponding clone, Merihb1, from Arabidopsis. Analysis shows that the predicted MERIHB1 protein exhibits strong homology to KN1 and RS1 from maize, SBH1 from soybean, and KNAT1 and KNAT2 from Arabidopsis. Merihb1 is highly expressed in mRNA from cauliflower meristems and also accumulates in stem and flower mRNA. Based on the similarity of the Merihb1 and Kn1 sequences, expression patterns, and in situ hybridizations, we suggest that Merihb1 represents an Arabidopsis homologue of the maize Kn1 gene.

Drosophila engrailed-1,10-phenanthroline chimeras as probes of homeodomain-DNA complexes

We have converted the Drosophila engrailed homeodomain into a sequence-specific nuclease by linking the protein to the chemical nuclease 1,10-phenanthroline-copper (OP-Cu). Unique cysteines were introduced at six positions into the homeodomain by site-directed mutagenesis for the covalent attachment of OP-Cu. The varied DNA-binding affinity and specificity of these mutants and the DNA cleavage pattern of their OP-Cu derivatives allowed us to assess the crystal structure of the engrailed homeodomain-DNA complex. We have also achieved site-specific double-stranded DNA scission with one of the homeodomain mutants, E28C, which has the potential of being used to identify engrailed binding sites in the genome. Because the homeodomain is so well conserved among members of the homeodomain-containing protein family, other homeodomain proteins can be converted into nucleases by attaching OP-Cu at position 28 of their homeodomains.

Phosphorylation of the Drosophila engrailed protein at a site outside its homeodomain enhances DNA binding

The engrailed gene encodes a homeodomain-containing phosphoprotein that binds DNA. Here, we show that engrailed protein is posttranslationally modified in embryos and in embryo-derived cultured cells but is essentially unmodified when expressed in Escherichia coli. Engrailed protein produced by bacteria can be phosphorylated in nuclear extracts prepared from Drosophila embryos, and phosphotryptic peptides from this modified protein partly reproduce two-dimensional maps of phosphotryptic fragments obtained from metabolically labeled engrailed protein. The primary embryonic protein kinase modifying engrailed protein is casein kinase II (CK-II). Analysis of mutant proteins revealed that the in vitro phosphoacceptors are mainly clustered in a region outside the engrailed homeodomain and identified serines 394, 397, 401, and 402 as the targets for CK-II phosphorylation. CK-II-dependent phosphorylation of an N-truncated derivative of engrailed protein purified from bacteria increased its DNA binding 2-4-fold.

Mutually exclusive interactions between factors binding to adjacent Sp1 and AT-rich elements regulate gastrin gene transcription in insulinoma cells

The gastrin gene is transiently expressed in fetal pancreatic islets during islet neogenesis but then switched off after birth when islet cells become fully differentiated. Previous studies identified a cis-regulatory sequence between -109 and -75 in the human gastrin promoter which binds islet cell-specific activators and a nonspecific repressor and thus may act as a molecular switch. The present study identified another cis-regulatory sequence (-163ACACTAAATGAAAGGGCGGGGCAG-140) which bound two islet nuclear proteins in a mutually exclusive manner, as defined by gel shift competition, methylation interference, and DNase I foot-printing assays. The general transactivator Sp1 recognized the downstream GGGCGGGG sequence, but Sp1 binding was prevented when another islet factor bound to the adjacent AT-rich sequence (CTAAATGA). This gastrin AT-rich element is nearly identical to the binding site (ATAAATGA) for the islet-specific transcription factor beta TF-1. However, the gastrin AT-binding factor appeared to differ from beta TF-1 in its gel mobility shift pattern. Transfections of rat insulinoma cells revealed that mutations which blocked binding to the AT-rich element but allowed Sp1 binding up-regulated transcriptional activity. These results suggest that the gastrin AT-binding factor blocks transactivation by Sp1 and may have a role in the repression of gastrin transcription seen at the end of islet differentiation.

Sequence, genomic organization, and expression of the novel homeobox gene Hesx1

Extensive analyses of homeobox gene expression and function during murine embryogenesis have demonstrated that homeobox gene products are key components in the establishment of pattern formation and regional identity during development. In this paper we report the molecular characterization and expression of a novel murine homeobox sequence, Hesx1, isolated from pluripotent embryonic stem cells. Hesx1 is expressed as two transcripts of 1.0 and 1.2 kilobases which encode an identical 185 amino acid open reading frame. The transcripts differ in the 3'-untranslated region due to the differential utilization of a weak splice donor site located immediately downstream of the translation termination codon. The Hesx1 homeodomain shared 80% identity with the Xenopus homeoprotein XANF-1 and was less than 50% related to other homeodomain sequences. Hesx1 and XANF-1 therefore constitute the founder members of a new homeodomain class. Hesx1 expression was down-regulated during embryonic stem cell differentiation and was detected in tissue-specific RNA samples derived from the embryonic liver, and at lower levels in viscera, amnion, and yolk sac. Expression in adult mice was not detected. These sites of expression are consistent with a role for Hesx1 in the regulation of developmental decisions in the early mouse embryo and during fetal hematopoiesis.

Structural studies of the engrailed homeodomain

The structure of the Drosophila engrailed homeodomain has been solved by molecular replacement and refined to an R-factor of 19.7% at a resolution of 2.1 A. This structure offers a high-resolution view of an important family of DNA-binding proteins and allows comparison to the structure of the same protein bound to DNA. The most significant difference between the current structure and that of the 2.8-A engrailed-DNA complex is the close packing of an extended strand against the rest of the protein in the unbound protein. Structural features of the protein not previously noted include a "herringbone" packing of 4 aromatic residues in the core of the protein and an extensive network of salt bridges that covers much of the helix 1-helix 2 surface. Other features that may play a role in stabilizing the native state include the interaction of buried carbonyl oxygen atoms with the edge of Phe 49 and a bias toward statistically preferred side-chain dihedral angles. There is substantial disorder at both ends of the 61 amino acid protein. A 51-amino acid variant of engrailed (residues 6-56) was synthesized and shown by CD and thermal denaturation studies to be structurally and thermodynamically similar to the full-length domain.

A homology-based molecular model of the proline-rich homeodomain protein Prh, from haematopoietic cells

A molecular structural model for the homeodomain of the haematopoietic protein Prh together with its DNA recognition sequence, has been built using the known crystal structure of the MAT alpha 2 homeodomain as a starting-point. The modelling procedure used main and side-chain optimisations by means of molecular mechanics/simulated annealing procedures to obtain stereochemically plausible geometries. The resulting structure has a number of specific interactions in both major and minor grooves of the DNA that serve to define the consensus binding sequence for Prh. In particular, the side-chain of glutamine 50 is postulated to be involved in hydrogen bonds to adjacent adenine and cytosine bases within the consensus sequence.

Homeodomain proteins in development and therapy

Homeobox genes encode transcriptional regulators found in all organisms ranging from yeast to humans. In Drosophila, a specific class of homeobox genes, the homeotic genes, specifies the identity of certain spatial units of development. Their genomic organization, in Drosophila, as well as in vertebrates, is uniquely connected with their expression which follows a 5'-posterior-3'-anterior rule along the longitudinal body axis. The 180-bp homeobox is part of the coding sequence of these genes, and the sequence of 60 amino acids it encodes is referred to as the homeodomain. Structural analyses have shown that homeodomains consist of a helix-turn-helix motif that binds the DNA by inserting the recognition helix into the major groove of the DNA and its amino-terminal arm into the adjacent minor groove. Developmental as well as gene regulatory functions of homeobox genes are discussed, with special emphasis on one group, the Antennapedia (Antp) class homeobox genes and a representative 60-amino acid Antennapedia peptide (pAntp). In cultured neuronal cells, pAntp translocates through the membrane specifically and efficiently and accumulates in the nucleus. The internalization process is followed by a strong induction of neuronal morphological differentiation, which raises the possibility that motoneuron growth is controlled by homeodomain proteins. It has been demonstrated that chimeric peptide molecules encompassing pAntp are also captured by cultured neurons and conveyed to their nuclei. This may be of enormous interest for the internalization of drugs.