Nuclear magnetic resonance spectroscopy of a DNA complex with the uniformly 13C-labeled Antennapedia homeodomain and structure determination of the DNA-bound homeodomain

Of numbers and movement - understanding transcription factor pathogenesis by advanced microscopy

Transcription factors (TFs) are life-sustaining and, therefore, the subject of intensive research. By regulating gene expression, TFs control a plethora of developmental and physiological processes, and their abnormal function commonly leads to various developmental defects and diseases in humans. Normal TF function often depends on gene dosage, which can be altered by copy-number variation or loss-of-function mutations. This explains why TF haploinsufficiency (HI) can lead to disease. Since aberrant TF numbers frequently result in pathogenic abnormalities of gene expression, quantitative analyses of TFs are a priority in the field. In vitro single-molecule methodologies have significantly aided the identification of links between TF gene dosage and transcriptional outcomes. Additionally, advances in quantitative microscopy have contributed mechanistic insights into normal and aberrant TF function. However, to understand TF biology, TF-chromatin interactions must be characterised in vivo, in a tissue-specific manner and in the context of both normal and altered TF numbers. Here, we summarise the advanced microscopy methodologies most frequently used to link TF abundance to function and dissect the molecular mechanisms underlying TF HIs. Increased application of advanced single-molecule and super-resolution microscopy modalities will improve our understanding of how TF HIs drive disease.

Racemic phosphorothioate as a tool for NMR investigations of protein-DNA complexes

A major driving force for protein-nucleic acid association is electrostatic interactions via ion pairs of the positively charged basic side chains and negatively charged phosphates. For a better understanding of how proteins scan DNA and recognize particular signatures, it is important to gain atomic-level insight into the behavior of basic side chains at the protein-DNA interfaces. NMR spectroscopy is a powerful tool for investigating the structural, dynamic, and kinetic aspects of protein-DNA interactions. However, resonance assignment of basic side-chain cationic moieties at the molecular interfaces remains to be a major challenge. Here, we propose a fast, robust, and inexpensive approach that greatly facilitates resonance assignment of interfacial moieties and also allows for kinetic measurements of protein translocation between two DNA duplexes. This approach utilizes site-specific incorporation of racemic phosphorothioate at the position of a phosphate that interacts with a protein side chain. This modification retains the electric charge of phosphate and therefore is mild, but causes significant chemical shift perturbations for the proximal protein side chains, which facilitates resonance assignment. Due to the racemic nature of the modification, two different chemical shifts are observed for the species with different diastereomers R_P and S_P of the incorporated phosphorothioate group. Kinetic information on the exchange of the protein molecule between R_P and S_P DNA duplexes can be obtained by ¹⁵N_z exchange spectroscopy. We demonstrate the applications of this approach to the Antennapedia homeodomain-DNA complex and the CREB1 basic leucine-zipper (bZIP)-DNA complex.

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.

Characterization and target genes of nine human PRD-like homeobox domain genes expressed exclusively in early embryos

PAIRED (PRD)-like homeobox genes belong to a class of predicted transcription factor genes. Several of these PRD-like homeobox genes have been predicted in silico from genomic sequence but until recently had no evidence of transcript expression. We found recently that nine PRD-like homeobox genes, ARGFX, CPHX1, CPHX2, DPRX, DUXA, DUXB, NOBOX, TPRX1 and TPRX2, were expressed in human preimplantation embryos. In the current study we characterized these PRD-like homeobox genes in depth and studied their functions as transcription factors. We cloned multiple transcript variants from human embryos and showed that the expression of these genes is specific to embryos and pluripotent stem cells. Overexpression of the genes in human embryonic stem cells confirmed their roles as transcription factors as either activators (CPHX1, CPHX2, ARGFX) or repressors (DPRX, DUXA, TPRX2) with distinct targets that could be explained by the amino acid sequence in homeodomain. Some PRD-like homeodomain transcription factors had high concordance of target genes and showed enrichment for both developmentally important gene sets and a 36 bp DNA recognition motif implicated in Embryo Genome Activation (EGA). Our data implicate a role for these previously uncharacterized PRD-like homeodomain proteins in the regulation of human embryo genome activation and preimplantation embryo development.

Comprehensive comparative homeobox gene annotation in human and mouse

Homeobox genes are a group of genes coding for transcription factors with a DNA-binding helix-turn-helix structure called a homeodomain and which play a crucial role in pattern formation during embryogenesis. Many homeobox genes are located in clusters and some of these, most notably the HOX genes, are known to have antisense or opposite strand long non-coding RNA (lncRNA) genes that play a regulatory role. Because automated annotation of both gene clusters and non-coding genes is fraught with difficulty (over-prediction, under-prediction, inaccurate transcript structures), we set out to manually annotate all homeobox genes in the mouse and human genomes. This includes all supported splice variants, pseudogenes and both antisense and flanking lncRNAs. One of the areas where manual annotation has a significant advantage is the annotation of duplicated gene clusters. After comprehensive annotation of all homeobox genes and their antisense genes in human and in mouse, we found some discrepancies with the current gene set in RefSeq regarding exact gene structures and coding versus pseudogene locus biotype. We also identified previously un-annotated pseudogenes in the DUX, Rhox and Obox gene clusters, which helped us re-evaluate and update the gene nomenclature in these regions. We found that human homeobox genes are enriched in antisense lncRNA loci, some of which are known to play a role in gene or gene cluster regulation, compared to their mouse orthologues. Of the annotated set of 241 human protein-coding homeobox genes, 98 have an antisense locus (41%) while of the 277 orthologous mouse genes, only 62 protein coding gene have an antisense locus (22%), based on publicly available transcriptional evidence.

Binding polymorphism in the DNA bound state of the Pdx1 homeodomain

The subtle effects of DNA-protein recognition are illustrated in the homeodomain fold. This is one of several small DNA binding motifs that, in spite of limited DNA binding specificity, adopts crucial, specific roles when incorporated in a transcription factor. The homeodomain is composed of a 3-helix domain and a mobile N-terminal arm. Helix 3 (the recognition helix) interacts with the DNA bases through the major groove, while the N-terminal arm becomes ordered upon binding a specific sequence through the minor groove. Although many structural studies have characterized the DNA binding properties of homeodomains, the factors behind the binding specificity are still difficult to elucidate. A crystal structure of the Pdx1 homeodomain bound to DNA (PDB 2H1K) obtained previously in our lab shows two complexes with differences in the conformation of the N-terminal arm, major groove contacts, and backbone contacts, raising new questions about the DNA recognition process by homeodomains. Here, we carry out fully atomistic Molecular Dynamics simulations both in crystal and aqueous environments in order to elucidate the nature of the difference in binding contacts. The crystal simulations reproduce the X-ray experimental structures well. In the absence of crystal packing constraints, the differences between the two complexes increase during the solution simulations. Thus, the conformational differences are not an artifact of crystal packing. In solution, the homeodomain with a disordered N-terminal arm repositions to a partially specific orientation. Both the crystal and aqueous simulations support the existence of different stable binding conformers identified in the original crystallographic data with different degrees of specificity. We propose that protein-protein and protein-DNA interactions favor a subset of the possible conformations. This flexibility in DNA binding may facilitate multiple functions for the same transcription factor.

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.

Comparing anterior and posterior Hox complex formation reveals guidelines for predicting cis-regulatory elements

Hox transcription factors specify numerous cell fates along the anterior-posterior axis by regulating the expression of downstream target genes. While expression analysis has uncovered large numbers of de-regulated genes in cells with altered Hox activity, determining which are direct versus indirect targets has remained a significant challenge. Here, we characterize the DNA binding activity of Hox transcription factor complexes on eight experimentally verified cis-regulatory elements. Hox factors regulate the activity of each element by forming protein complexes with two cofactor proteins, Extradenticle (Exd) and Homothorax (Hth). Using comparative DNA binding assays, we found that a number of flexible arrangements of Hox, Exd, and Hth binding sites mediate cooperative transcription factor complexes. Moreover, analysis of a Distal-less regulatory element (DMXR) that is repressed by abdominal Hox factors revealed that suboptimal binding sites can be combined to form high affinity transcription complexes. Lastly, we determined that the anterior Hox factors are more dependent upon Exd and Hth for complex formation than posterior Hox factors. Based upon these findings, we suggest a general set of guidelines to serve as a basis for designing bioinformatics algorithms aimed at identifying Hox regulatory elements using the wealth of recently sequenced genomes.

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 solution structure of the native K50 Bicoid homeodomain bound to the consensus TAATCC DNA-binding site

The solution structure of the homeodomain of the Drosophila morphogenic protein Bicoid (Bcd) complexed with a TAATCC DNA site is described. Bicoid is the only known protein that uses a homeodomain to regulate translation, as well as transcription, by binding to both RNA and DNA during early Drosophila development; in addition, the Bcd homeodomain can recognize an array of different DNA sites. The dual functionality and broad recognition capabilities signify that the Bcd homeodomain may possess unique structural/dynamic properties. Bicoid is the founding member of the K50 class of homeodomain proteins, containing a lysine residue at the critical 50th position (K50) of the homeodomain sequence, a residue required for DNA and RNA recognition; Bcd also has an arginine residue at the 54th position (R54), which is essential for RNA recognition. Bcd is the only known homeodomain with the K50/R54 combination of residues. The Bcd structure indicates that this homeodomain conforms to the conserved topology of the homeodomain motif, but exhibits a significant variation from other homeodomain structures at the end of helix 1. A key result is the observation that the side-chains of the DNA-contacting residues K50, N51 and R54 all show strong signs of flexibility in the protein-DNA interface. This finding is supportive of the adaptive-recognition theory of protein-DNA interactions.

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.

Conformational analysis by NMR and distance geometry techniques of a peptide mimetic of the third helix of the Antennapedia homeodomain

The Antennapedia homeodomain structure consists of four helices. The helices II and III are connected by a tripeptide that forms a turn, and constitute the well-known helix-turn-helix motif. The recognition helix penetrates the DNA major groove, gives specific protein-DNA contacts and forms direct, or water-mediated, intermolecular hydrogen bonds. It was suggested that helix III (and perhaps also helix IV) might represent the recognition helix of Antennapedia homeodomain, which makes contact with the surface of the major groove of the DNA. In an attempt to clarify the helix III capabilities of assuming an helical conformation when separated from the rest of the protein, we carried out the structural determination of the recognition helix III in different solvent media. The conformational study of fragments 42-53, where residues W48 and F49, not involved in the protein-DNA interaction, were substituted by two alanines, was conducted in sodium dodecyl sulfate (SDS), trifluoroethanol (TFE) and TFE/water, using circular dichroism, nuclear magnetic resonance (NMR) and distance geometry (DG) techniques. The fragment assumes a well-defined secondary structure in TFE and in TFE/water (90/10, v/v) with an alpha-helix encompassing residues 4-9, while in TFE/water (70/30, v/v) a less regular structure was found. The DG results in the micellar system evidence the presence of a distorted alpha-helical conformation involving residues 4-8. Our results reveal that the isolated Antennapedia recognition helix III tend to preserve in solution the alpha-helical conformation even if separated from the rest of the molecule.

Specific DNA recognition by theAntp homeodomain: MD simulations of specific and nonspecific complexes

Four molecular dynamics simulation trajectories of complexes between the wild-type or a mutant Antennapedia homeodomain and 2 DNA sequences were generated in order to probe the mechanisms governing the specificity of DNA recognition. The starting point was published affinity measurements showing that a single protein mutation combined with a replacement of 2 base pairs yields a new high-affinity complex, whereas the other combinations, with changes on only 1 macromolecule, exhibited lower affinity. The simulations of the 4 complexes yielded fluctuating networks of interaction. On average, these networks differ significantly, explaining the switch of affinity caused by the alterations in the macromolecules. The network of mostly hydrogen-bonding interactions involving several water molecules, which was suggested both by X-ray and NMR structures of the wild-type homeodomain and its DNA operator sequence, could be reproduced in the trajectory. More interestingly, the high-affinity complex with alterations in both the protein and the DNA yielded again a dynamic but very tight network of intermolecular interactions, however, attributing a significantly stronger role to direct hydrophobic interactions at the expense of water bridges. The other 2 homeodomain-DNA complexes, with only 1 molecule altered, show on average over the trajectories a clearly reduced number of protein-DNA interactions. The observations from these simulations suggest specific experiments and thus close the circle formed by biochemical, structural, and computational studies. The shift from a water-dominated to a more "dry" interface may prove important in the design of proteins binding DNA in a specific manner.

Insights into binding cooperativity of MATa1/MATalpha2 from the crystal structure of a MATa1 homeodomain-maltose binding protein chimera

The Yeast MATa1 and MATalpha2 are homeodomain proteins that bind DNA cooperatively to repress transcription of cell type specific genes. The DNA affinity and specificity of MATa1 in the absence of MATalpha2, however, is very low. MATa1 is converted to a higher affinity DNA-binding protein by its interaction with the C-terminal tail of MATalpha2. To understand why MATa1 binds DNA weakly by itself, and how the MATalpha2 tail affects the affinity of MATa1 for DNA, we determined the crystal structure of a maltose-binding protein (MBP)-a1 chimera whose DNA binding behavior is similar to MATa1. The overall MATa1 conformation in the MBP-a1 structure, which was determined in the absence of alpha2 and DNA, is similar to that in the a1/alpha2/DNA structure. The sole difference is in the C-terminal portion of the DNA recognition helix of MATa1, which is flexible in the present structure. However, these residues are not in a location likely to be affected by binding of the MATalpha2 tail. The results argue against conformational changes in a1 induced by the tail of MATalpha2, suggesting instead that the MATalpha2 tail energetically couples the DNA binding of MATalpha2 and MATa1.

Investigation of penetratin peptides. Part 1. The environment dependent conformational properties of penetratin and two of its derivatives

The homeodomain, the DNA-binding domain of Antennapedia homeoprotein, is composed of three alpha-helices and one beta-turn between helices II and III. Its third helix from the N-terminal (helix III) can translocate through the cell membrane into the nucleus and can be used as an intracellular vehicle for the delivery of oligopeptides and oligonucleotides. To the best of our knowledge, this helix III, called penetratin, which consists of 16 amino acids, is internalized by cells in a specific, non-receptor-mediated manner. For a better understanding of the mechanism of the transfer, the structure of penetratin was examined in both extracellular matrix-mimetic and membrane-mimetic environments: 1H-NMR and CD spectroscopic measurements were performed in mixtures of TFE/water with different ratios. The molecular conformations of two analogue peptides [(6,14-Phe)-penetratin and a 12 amino acid penetratin derivative (peptide 3)] were also studied. An atomic level comprehensive analysis of penetratin and its two analogues was performed. In a membrane-mimetic solvent system (TFEd2/water = 9: 1), on the basis of 553 distance restraints, the 4-12 region of penetratin exhibits a bent, irregular helical structure on NMR examination. Interactions between hydrophobic amino acid residues in conjunction with H-bonds stabilize the secondary structure of the molecule. Thus, both derivatives adopt a helix-like conformation. However, while (6,14-Phe)-penetratin displays both alpha-helical and 310-helical features, the structure of peptide 3 is predominantly a 310-helix. Of the three peptides, surprisingly (6,14-Phe)-penetratin has the largest helical content. An increase in the polarity of the molecular environment gradually disintegrates these helix-like secondary structures. In a highly aqueous molecular system (TFEd2/water = 1 : 9), the fast exchange of multiple conformers leads to too few distance restraints being extracted, therefore the NMR structures can no longer be determined. The NMR data show that only short-range order can be traced in these peptides. Under these conditions, the molecules adopt nascent helix-like structures. On the other hand, CD spectra could be recorded at any TFE/water ratio and the conformational interconversion could therefore be monitored as a function of the polarity of the molecular environment. The CD data were analysed comprehensively by the quantitative deconvolution method (CCA+). All three penetratin peptides display helical conformational features in a low dielectric medium, with significant differences as a function of their amino acid composition. However, these conformational features are gradually lost during the shift from an apolar to a polar molecular environment.

Identification of a dominant negative homeodomain mutation in Rieger syndrome

Mutations in the PITX2 bicoid-like homeobox gene cause Rieger syndrome. Rieger syndrome is an autosomal-dominant human disorder characterized by glaucoma as well as dental hypoplasia, mild craniofacial dysmorphism, and umbilical stump abnormalities. PITX2 has also been implicated in the development of multiple organs and left-right asymmetry in the body plan. The PITX2 homeodomain has a lysine at position 50, which has been shown to impart the bicoid-type (TAATCC) DNA binding specificity to other homeodomain proteins. A mutation (K88E), found in a Rieger syndrome patient, changes this lysine to glutamic acid. We were intrigued by the relatively pronounced phenotypic consequences of this K88E mutation. In the initial analyses, the mutant protein appeared to simply be inactive, with essentially no DNA binding and transactivation activities and, unlike the wild type protein, with an inability to synergize with another transcription factor, Pit-1. However, when the K88E DNA was cotransfected with wild type PITX2, analogous to the patient genotype, the K88E mutant suppressed the synergism of wild type PITX2 with Pit-1. In contrast, a different PITX2 homeodomain mutant, T68P, which is also defective in DNA binding, transactivation, and Pit-1 synergism activities, did not suppress the wild type synergism with Pit-1. These results describe the first dominant negative missense mutation in a homeodomain and support a model that may partially explain the phenotypic variation within Rieger syndrome.

Time-shared X(omega(1))-half-filter for improved sensitivity in subspectral editing

Experiments with X-half-filter elements allow the separation of the resonances from protons bound and unbound to a spin X into different subspectra. This Communication presents a modified half-filter element where the filter delay is simultaneously used for chemical shift labeling and scalar coupling evolution in a semi-constant time experiment. The filter element is demonstrated with a (1)H NOESY spectrum of a 28.5-kDa 2:1 complex between the uniformly (13)C-labeled N-terminal domain of Escherichia coli arginine repressor and operator DNA.

Comparative stability analysis of the thyroid transcription factor 1 and Antennapedia homeodomains: evidence for residue 54 in controlling the structural stability of the recognition helix

The thyroid transcription factor 1 homeodomain (TTF-1 HD) shows a peculiar DNA-binding specificity which is partially dictated by several amino acids of the recognition helix. TTF-1 preferentially recognizes sequences containing the 5'-CAAG-3' core motif while most other homeodomains, such as Antennapedia (Antp), recognizes sites containing the 5'-TAAT-3' core motif. Since phenomena of 'induced fit' may occur during protein/DNA interaction, a primary role for high affinity binding and target discrimination has to be searched in the effect played by subtle structural determinants in these proteins. By using spectroscopic analysis in aqueous solution, we compared the structural stability of TTF-1 and Antp homeodomains. Although the three-dimensional structural architecture of homeodomains is conserved, some differences are detectable in terms of their structural stability. At 24 degrees C the TTF-1 HD is less structured than the Antp HD with 24 and 34% of the residues in the alpha-helical conformation, respectively. This poor folded structure reflects into different thermal and isothermal stability between the two homeodomains. TTF-1 HD exhibits a Tm of 39 degrees C and is stabilized by a delta GDH2O of +1487 cal/mol, calculated by Urea unfolding, while Antp HD exhibits a Tm of 48 degrees C and is stabilized by a delta GDH2O of +2742 cal/mol. By using mutants of both TTF-1 and Antp HDs we demonstrate that one of the major determinants in controlling the structural stability of the recognition helix is the residue at position 54. Since previous studies have shown that also residue at position 56 is involved in stabilization of the recognition helix, we conclude that the structure of this critical element is controlled by an interplay between residues at position 54 and 56 of the homeodomain.

Conformational changes of the BS2 operator DNA upon complex formation with the Antennapedia homeodomain studied by NMR with 13C/15N-labeled DNA

The NMR structures have been determined for a 13C/15N doubly labeled 14 base-pair DNA duplex comprising the BS2 operator sequence both free in solution and in the complex with the Antennapedia homeodomain. The impact of the DNA labeling is assessed from comparison with a previous structure of the same complex that was determined using isotope labeling only for the protein. Differences between the two structure determinations are nearly completely limited to the DNA, which retains the global B -conformation of the free DNA also in the complex. Local protein-induced conformational changes are a narrowing of the minor groove due to the interaction with the N-terminal arm of the homeodomain, and changes of the sugar puckers of the deoxyriboses G5 and C6, which are apparently induced by van der Waals interactions with Tyr25, and with Gln50 and Arg53, respectively. The high conservation of these amino acid residues in homeodomains suggests that protein-induced shifts in some sugar puckers contribute to the affinity of homeodomains to their cognate DNA. The data obtained here with the Antennapedia homeodomain-DNA complex clearly show that nucleic acid isotope-labeling can support detailed conformational characterization of DNA in complexes with proteins, which will be indispensable for structure determinations of complexes containing globally distorted DNA conformations.

The solution structure of the homeodomain of the rat insulin-gene enhancer protein isl-1. Comparison with other homeodomains

Homeodomains are one of the key families of eukaryotic DNA-binding motifs and provide an important model system for DNA recognition. We have determined a high-quality nuclear magnetic resonance (NMR) structure of the DNA-binding homeodomain of the insulin gene enhancer protein Isl-1 (Isl-1-HD). It forms the first solution structure of a homeodomain from the LIM family. It contains a well-defined inner core (residues 12-55) consisting of the classical three-helix structure observed in other homeodomains. The N terminus is unstructured up to residue 8, while the C terminus gradually becomes unstructured from residue 55 onwards. Some flexibility is evident in the loop parts of the inner core. Isl-1-HD has, despite its low sequence identity (23-34 %), a structure that is strikingly similar to that of the other homeodomains with known three-dimensional structures. Detailed analysis of Isl-1-HD and the other homeodomains rationalizes the differences in their temperature stability and explains the low stability of the Isl-1-HD in the free state (tm 22-30 degrees C). Upon DNA binding, a significant stabilization occurs (tm>55 degrees C). The low stability of Isl-1-HD (and other mammalian homeodomains) suggests that in vivo Isl-1-HD recognizes its cognate DNA from its unfolded state.

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

Hox homeodomain proteins are developmental regulators that determine body plan in a variety of organisms. A majority of the vertebrate Hox proteins bind DNA as heterodimers with the Pbx1 homeodomain protein. We report here the 2.35 A structure of a ternary complex containing a human HoxB1-Pbx1 heterodimer bound to DNA. Heterodimer contacts are mediated by the hexapeptide of HoxB1, which binds in a pocket in the Pbx1 protein formed in part by a three-amino acid insertion in the Pbx1 homeodomain. The Pbx1 DNA-binding domain is larger than the canonical homeodomain, containing an additional alpha helix that appears to contribute to binding of the HoxB1 hexapeptide and to stable binding of Pbx1 to DNA. The structure suggests a model for modulation of Hox DNA binding activity by Pbx1 and related proteins.

Mutation analysis of LMX1B gene in nail-patella syndrome patients

Nail-patella syndrome (NPS), a pleiotropic disorder exhibiting autosomal dominant inheritance, has been studied for >100 years. Recent evidence shows that NPS is the result of mutations in the LIM-homeodomain gene LMX1B. To determine whether specific LMX1B mutations are associated with different aspects of the NPS phenotype, we screened a cohort of 41 NPS families for LMX1B mutations. A total of 25 mutations were identified in 37 families. The nature of the mutations supports the hypothesis that NPS is the result of haploinsufficiency for LMX1B. There was no evidence of correlation between aspects of the NPS phenotype and specific mutations.

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.

15N-edited three-dimensional NOESY-HMQC with water flipback: enhancement of weak labile 1H resonances of protein side chains contacting DNA

Two pulse sequences are described that employ a modified water flipback technique to enhance the signal intensity of weak side chain resonances at the protein-DNA interface of the vnd/NK-2 homeodomain/DNA complex in an 15N-edited three-dimensional NOESY-HMQC spectrum. The pulse sequences presented employ water flipback pulses at the beginning of the NOESY mixing time, optimizing the direct NOE transfer of magnetization from the water to the protein by maximizing the z-component of the water magnetization. In one of the pulse sequences, radiation damping during the the indirect 1H and 15N evolution times is suppressed. A modified version of the WATERGATE water suppression technique is employed during the HMQC portion of the experiment. The signal enhancement is demonstrated for the resonances of the side chain amide of Asn51, an invariant homeodomain residue whose contact with the DNA is critical for binding. An ancillary advantage of the experiment is the ability to observe NOE transfer of magnetization from water. The information present in the water resonance plane of the three-dimensional spectrum is illustrated in a comparison with the corresponding HMQC spectrum of the protein/DNA complex.

Comparison of X-ray and NMR structures for the Antennapedia homeodomain-DNA complex

Homeodomains are one of the key families of eukaryotic DNA-binding motifs and provide an important model system for studying protein-DNA interactions. We have crystallized the Antennapedia homeodomain-DNA complex and solved this structure at 2.4 A resolution. NMR and molecular dynamics studies had implied that this homeodomain achieves specificity through an ensemble of rapidly fluctuating DNA contacts. The crystal structure is in agreement with the underlying NMR data, but our structure reveals a well-defined set of contacts and also reveals the locations and roles of water molecules at the protein-DNA interface. The synthesis of X-ray and NMR studies provides a unified, general model for homeodomain-DNA interactions.

A conserved cluster of homeodomain binding sites in the mouse Hoxa-4 intron functions in Drosophila embryos as an enhancer that is directly regulated by Ultrabithorax

The evolutionary conservation of the homeodomains suggests that their in vivo DNA binding sites may also be conserved between vertebrates and invertebrates. The regulatory function of the mouse Hoxa-4 and Hoxb-4 introns were analyzed in Drosophila since they both contain a cluster of three homeodomain binding sites, the HB1 element, which was also found in the introns of other Hox genes ranging from fish to humans as well as in the Ultrabithorax (Ubx) and decapentaplegic (dpp) genes of Drosophila. The enhancer of the Hoxa-4 intron was found to respond to several homeobox genes activating a lacZ reporter gene in particular cells of the epidermis in Drosophila embryos. The enhancer activity was found to be similar to previously described autoregulatory elements of Deformed (Dfd), the Drosophila homolog of Hoxa-4, but additional expression was observed in more posterior segments activated by Ubx and repressed by abdominal-A (abd-A). Point mutations in the homeodomain binding sites in HB1 abolished the enhancer activity. A second site suppression experiment showed that UBX interacts directly with the HB1 element. When the HB1 element in the Hoxa-4 intron was replaced by that of the mesodermal enhancer of dpp, which was previously shown to be directly controlled by Ubx, Ubx-dependent activation was retained, but repression by abd-A was lost. The same result was obtained when the third binding site of HB1 was altered, suggesting that this site is responsible for abd-A-dependent repression. Finally, deletion of potential cofactor binding sites flanking the HB1 element that are also conserved in the medaka, chicken, and mouse genes revealed that they are important for enhancer function in Drosophila and that the Dfd-dependent and the Ubx-dependent expression requires different sites. The evolutionary and functional conservation of the HB1 elements indicates that not only the homeodomains but also some of their in vivo binding sites are conserved between vertebrates and invertebrates.

Electrostatics and hydration at the homeodomain-DNA interface: chemical probes of an interfacial water cavity

Electrostatics and hydration of a homeodomain-DNA complex are dissected by chemical modification. Selective neutralization of phosphate charges by methylphosphonate substitution demonstrates the differential importance of short- and long-range electrostatic interactions. Whereas the footprint of direct contacts is in accord with crystal structures, interference is also observed at non-contacted sites. Such sites adjoin a novel interfacial water cavity in the major groove. Non-contacted phosphodiester groups in the cavity are proposed to contribute to long-range ordering of an extended protein-water-DNA interface. Use of isolated S(p) and R(p) methylphosphonate diastereomers demonstrates that interference at this extended interface is stereoselective and charge-independent. Attenuation of protein binding presumably reflects groove-specific reorganization of bound water. Surprisingly, such attenuation can exceed that due to neutralization of a direct phosphate-side-chain salt bridge. These results support the hypothesis that hydration of an interfacial cavity functions as a non-covalent extension of the DNA surface. Stereo-specific interrogation of bound water by chemical synthesis provides a general method to assess the coupling between solvation and DNA recognition.

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.

Solution structure of the Grb2 N-terminal SH3 domain complexed with a ten-residue peptide derived from SOS: direct refinement against NOEs, J-couplings and 1H and 13C chemical shifts

Refined ensembles of solution structures have been calculated for the N-terminal SH3 domain of Grb2 (N-SH3) complexed with the ac-VPPPVPPRRR-nh2 peptide derived from residues 1135 to 1144 of the mouse SOS-1 sequence. NMR spectra obtained from different combinations of both 13C-15N-labeled and unlabeled N-SH3 and SOS peptide fragment were used to obtain stereo-assignments for pro-chiral groups of the peptide, angle restraints via heteronuclear coupling constants, and complete 1H, 13C, and 15N resonance assignments for both molecules. One ensemble of structures was calculated using conventional methods while a second ensemble was generated by including additional direct refinements against both 1H and 13C(alpha)/13C(beta) chemical shifts. In both ensembles, the protein:peptide interface is highly resolved, reflecting the inclusion of 110 inter-molecular nuclear Overhauser enhancement (NOE) distance restraints. The first and second peptide-binding sub-sites of N-SH3 interact with structurally well-defined portions of the peptide. These interactions include hydrogen bonds and extensive hydrophobic contacts. In the third highly acidic sub-site, the conformation of the peptide Arg8 side-chain is partially ordered by a set of NOE restraints to the Trp36 ring protons. Overall, several lines of evidence point to dynamical averaging of peptide and N-SH3 side-chain conformations in the third subsite. These conformations are characterized by transient charge stabilized hydrogen bond interactions between the peptide arginine side-chain hydrogen bond donors and either single, or possibly multiple, acceptor(s) in the third peptide-binding sub-site.

Eukaryotic transcription factor-DNA complexes

Eukaryotes have three distinct RNA polymerases that catalyze transcription of nuclear genes. RNA polymerase II is responsible for transcribing nuclear genes encoding the messenger RNAs and several small nuclear RNAs. Like RNA polymerases I and III, polymerase II cannot recognize its target promoter directly and initiate transcription without accessory factors. Instead, this large multisubunit enzyme relies on general transcription factors and transcriptional activators and coactivators to regulate transcription from class II promoters. X-ray crystallography and nuclear magnetic resonance spectroscopy have been used to study complexes of general transcription factors and transcriptional activators with their specific DNA targets. This work has provided important structural insights into transcription initiation by polymerase II and the more general problem of DNA sequence recognition.

Kinetics of binding of Antp homeodomain to DNA analyzed by measurements of surface plasmon resonance

The kinetics of binding of the Antp homeodomain to three kinds of DNA fragment were analyzed by measurements of surface plasmon resonance at various temperatures. Non-specific and specific binding of the homeodomain to DNA was examined. In the case of non-specific binding, the association rate constant (k(ass)) was estimated to be 1.41-2.62 x 10(5) M(-1) s(-1) and the dissociation rate constant (k(diss)) was 1.36-3.10 x 10(-2) s(-1), thus, the dissociation constant (KD) was 0.847-1.72 x 10(-7) M. The association seemed to be driven by entropy. In the case of specific binding, by contrast, the enthalpy term seemed to contribute more to the binding than did the entropy term. The k(ass) was 2.04-2.59 x 10(5) M(-1) s(-1) and the k(diss) was 0.759-1.16 x 10(-3) s(-1), thus, the KD was 2.93-5.69 X 10(-9) M. These values were measured under the condition of 150 mM NaCl. Both interactions were strongly dependent on the concentration of NaCl. The KD at 50 mM NaCl became several tens of times smaller than those at 150 mM. Possible reasons for the differences between non-specific and specific interactions are discussed.

Conformational and associative behaviours of the third helix of antennapedia homeodomain in membrane-mimetic environments

The third helix of antennapedia homeodomain pAntp-(43-58) can translocate through cell membrane and has been used as an intracellular vehicle for delivering peptides and oligonucleotides. The conformational and associative behaviour of two peptidic vectors pAntp-(43-58) and [Pro50] pAntp-(43-58) has been analyzed by different biophysical methods. pAntp-(43-58) adopts an amphipathic helical structure in 30% (by vol.) hexafluoroisopropanol, in perfluoro-tert-butanol and in the presence of SDS micelles. CD spectra indicate that the conformation of [Pro50]pAntp-(43-58) in contrast to pAntp-(43-58) is independent of the media used. 1H-NMR spectroscopy in SDS micelles or in perfluoro-tert-butanol allows detection of aggregated peptides probably in a ribbon 2(7) type conformation. These conformations became the predominant structure when Gln50 was replaced by Pro50. Interproton-distance restraints derived from NOE measurements have been classified in two groups corresponding to two types of structures: alpha-helix and essentially extended structures. Consecutive CH alpha (i)/ CH alpha (i + 1) NOEs are only compatible with aggregates. Simulated annealing calculation of dimeric structure agrees with phi and psi angles in the beta-sheet and gamma-turn regions. Fluorescence spectroscopy analysis has shown that the indole groups of both peptides penetrate into SDS micelles; both peptides also induce the formation of micelles at very low concentration of SDS (20 microM). Similar interaction was observed with reverse-phase micelles made of bis(2-ethyhexyl) sodium sulfosuccinate and small unilamellar vesicles (SUV) made of a mixture of phosphatidylcholine/phosphatidylserine. 31P-NMR of vesicles (SUV and large unilamellar vesicles) indicated that the addition of pAntp analogues did not affect the size of phosphatidylcholine/phosphatidylserine vesicles. The addition of pAntp analogues to lipidic dispersions modulates lipid polymorphism in different ways depending on the mixtures of acidic lipids.

Intron of the mouse Hoxa-7 gene contains conserved homeodomain binding sites that can function as an enhancer element in Drosophila

The 5' flanking sequences and the intron of the mouse Hoxa-7 gene were searched for regulatory elements that can function in Drosophila. Only the intron is able to activate a lacZ fusion gene in various tissues of Drosophila embryos. This enhancer function requires a cluster of three homeodomain binding sites (HB1-element) that are also found in the introns of other Hox genes as well as in a putative autoregulatory element of Ultrabithorax (Ubx), the Drosophila homolog of Hoxa-7. If a single binding site in the autoregulatory element of fushi tarazu (ftz) is replaced by the HB1-element of Hoxa-7, the expression pattern is altered and newly controlled by the homeotic gene caudal (cad). These data suggest that HB1 is a potential target for different homeodomain proteins of both vertebrates and invertebrates.

Analysis of the solution structure of the homeodomain of rat thyroid transcription factor 1 by 1H-NMR spectroscopy and restrained molecular mechanics

The solution structure of the rat thyroid transcription factor 1 (TTF-1) homeodomain has been elucidated by 1H-NMR and restrained modeling. The TTF-1 homeodomain folds in the same manner as classical homeodomains, with three helices, a loose loop between the first two helices, and a tight turn between helix II and helix III. The typical assembly of the hydrophobic core is maintained and N-capping motifs are identified in helix I and helix III. The N-terminal stretch of helix II exhibits some mobility, similar to the preceding loop region, which may be related to its anomalous capping. The N-terminal decapeptide and the C-terminal octapeptide of the molecule (68 residues long) are disordered. All the previous characteristics are shared by all known isolated homeodomain structures. An important difference among these structures occurs at the C-terminal extension of helix III, which is either disordered or helically folded. In the TTF-1 homeodomain, the C-terminal extension of helix III (residues 51-59) appears structured, albeit not as rigidly as the preceding portion. Analysis of the NOEs and hydrogendeuterium exchange of backbone amides provides evidence for discontinuity between the two moieties of helix III, which is introduced by a tightening or a kink of residues 51-53.

Drosophila development: homeodomains and translational control

Translation of the transcription factor caudal is repressed at the anterior end of the Drosophila embryo. Surprisingly, the DNA-binding homeodomain of the transcription factor Bicoid mediates this repression by binding caudal mRNA.

Hydration and DNA recognition by homeodomains

A 2-nanosecond molecular dynamics (MD) simulation of an Antennapedia homeodomain-DNA complex in explicit solvent water at ambient temperature and pressure was performed to supplement experimental nuclear magnetic resonance (NMR) data on the structure and dynamics of this complex. In addition to direct protein-DNA contacts, the MD trajectory attributes an essential role for specific DNA recognition to hydration water molecules that mediate intermolecular contacts. The simulation provides a detailed description of the pathways of hydration water molecules exchanging in and out of the protein-DNA interface and indicates that the residence times of these "interior" waters are on the nanosecond time scale, near the lower end of the range determined by NMR.

The crystal structure of the DNA-binding domain of yeast RAP1 in complex with telomeric DNA

Telomeres, the nucleoprotein complexes at the ends of eukaryotic chromosomes, are essential for chromosome stability. In the yeast S. cerevisiae, telomeric DNA is bound in a sequence-specific manner by RAP1, a multifunctional protein also involved in transcriptional regulation. Here we report the crystal structure of the DNA-binding domain of RAP1 in complex with telomeric DNA site at 2.25 A resolution. The protein contains two similar domains that bind DNA in a tandem orientation, recognizing a tandemly repeated DNA sequence. The domains are structurally related to the homeodomain and the proto-oncogene Myb, but show novel features in their DNA-binding mode. A structured linker between the domains and a long C-terminal tail contribute to the binding specificity. This structure provides insight into the recognition of the conserved telomeric DNA sequences by a protein.

Homeodomain interactions

Homeodomain proteins play key roles in development and gene regulation in eukaryotes. Past structural studies have focused on the binding of monomeric homeodomains to DNA, but two recent structures have revealed how homeodomains bind DNA as multimers. The structures of the Drosophila Paired homodimer and the yeast a1/alpha2 heterodimer bound to DNA, along with a high-resolution study of a Drosophila eve-DNA complex, have deepened our understanding of how homeodomains locate their DNA targets.

Solution structure of the Oct-1 POU homeodomain determined by NMR and restrained molecular dynamics

The POU homeodomain (POUhd), a divergent member of the well-studied class of homeodomain proteins, is the C-terminal part of the bipartite POU domain, the conserved DNA-binding domain of the POU proteins. In this paper we present the solution structure of POUhd of the human Oct-1 transcription factor. This fragment was overexpressed in Escherichia coli and studied by two- and three-dimensional homo- and heteronuclear NMR techniques, resulting in virtually complete 1H and 15N resonance assignments for residues 2-60. Using distance and dihedral constraints derived from the NMR data, 50 distance geometry structures were calculated, which were refined by means of restrained molecular dynamics. From this set a total of 31 refined structures were selected, having low constraint energy and few constraint violations. The ensemble of 31 structures displays a root-mean-square deviation of the coordinates of 0.59 A with respect to the average structure, calculated over the backbone atoms of residues 6 to 54. The fold of POUhd is very similar to that of the canonical homeodomains. Interestingly, the recognition helix of the free POUhd ends at residue 53, while in the cocrystal structure of the intact POU domain with the DNA octamer motif [Klemm, J.D., Rould, M.A., Aurora, R., Herr, W. and Pabo, C.O. (1994) Cell, 77, 21-32] this helix in the POUhd subdomain is extended as far as residue 60.

Structure of the Oct-3 POU-homeodomain in solution, as determined by triple resonance heteronuclear multidimensional NMR spectroscopy

The POU-homeodomain (POUH) forms the bipartite DNA-binding POU domain in association with the POU-specific domain. The 1H, 15N, and 13C magnetic resonances of the 67-amino acid long POUH of mouse Oct-3 have almost completely been assigned, mainly through the combined use of three-dimensional triple resonance NMR methods. Based on the distance and dihedral angle constraints derived from the NMR data, the solution structure of the POUH domain has been calculated by the ab initio simulated annealing method. The average RMS deviation for all backbone heavy atoms of the 20 best calculated structures for residues 9-53 of the total 67 amino acid residues is 0.44 A. The POUH domain consists of three alpha-helices (helix-I, 10-20; helix-II, 28-38; and helix-III, 42-53), and helices-II and -III form a helix-turn-helix motif. In comparison with other classical homeodomains, the folding of the three helices is quite similar. However, the length of helix-III is fairly short. In the complex of the Oct-1 POU domain with an octamer site (Klemm JD, et al., 1994, Cell 77:21-32), the corresponding region is involved in helix-III. The structural difference between these two cases will be discussed.

The des(1-6)antennapedia homeodomain: comparison of the NMR solution structure and the DNA-binding affinity with the intact Antennapedia homeodomain

The nuclear magnetic resonance (NMR) solution structure of an N-terminally truncated mutant Antennapedia homeodomain, des(1-6)Antp(C39S), has been determined from 935 nuclear Overhauser effect upper distance constraints and 148 dihedral angle constraints by using the programs DIANA and OPAL. Twenty conformers representing the solution structure of des(1-6)Antp(C39S) have an average root-mean-square distance relative to the mean coordinates of 0.56 A for the backbone atoms of residues 8-59. Comparison with the intact Antp(C39S) homeodomain shows that the two proteins have identical molecular architectures. The removal of the N-terminal residues 1-6, which are flexibly disordered in the intact homeodomain, causes only strictly localized structure variations and does not noticeably affect the adjoining helix I from residues 10-21. The DNA-binding constant of des(1-6)Antp(C39S) is approximately 10-fold reduced relative to the intact Antp(C39S) homeodomain, which can now be attributed to the absence of the previously reported contacts of the N-terminal polypeptide segment of the intact Antp(C39S) homeodomain with the minor groove of the DNA duplex.

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.