Mutations that affect the ability of the vnd/NK-2 homeoprotein to regulate gene expression: transgenic alterations and tertiary structure

Major β cell-specific functions of NKX2.2 are mediated via the NK2-specific domain

The consolidation of unambiguous cell fate commitment relies on the ability of transcription factors (TFs) to exert tissue-specific regulation of complex genetic networks. However, the mechanisms by which TFs establish such precise control over gene expression have remained elusive-especially in instances in which a single TF operates in two or more discrete cellular systems. In this study, we demonstrate that β cell-specific functions of NKX2.2 are driven by the highly conserved NK2-specific domain (SD). Mutation of the endogenous NKX2.2 SD prevents the developmental progression of β cell precursors into mature, insulin-expressing β cells, resulting in overt neonatal diabetes. Within the adult β cell, the SD stimulates β cell performance through the activation and repression of a subset of NKX2.2-regulated transcripts critical for β cell function. These irregularities in β cell gene expression may be mediated via SD-contingent interactions with components of chromatin remodelers and the nuclear pore complex. However, in stark contrast to these pancreatic phenotypes, the SD is entirely dispensable for the development of NKX2.2-dependent cell types within the CNS. Together, these results reveal a previously undetermined mechanism through which NKX2.2 directs disparate transcriptional programs in the pancreas versus neuroepithelium.

Cognate DNA Recognition by Engrailed Homeodomain Involves a Conformational Change Controlled via an Electrostatic-Spring-Loaded Latch

Transcription factors must scan genomic DNA, recognize the cognate sequence of their control element(s), and bind tightly to them. The DNA recognition process is primarily carried out by their DNA binding domains (DBD), which interact with the cognate site with high affinity and more weakly with any other DNA sequence. DBDs are generally thought to bind to their cognate DNA without changing conformation (lock-and-key). Here, we used nuclear magnetic resonance and circular dichroism to investigate the interplay between DNA recognition and DBD conformation in the engrailed homeodomain (enHD), as a model case for the homeodomain family of eukaryotic DBDs. We found that the conformational ensemble of enHD is rather flexible and becomes gradually more disordered as ionic strength decreases following a Debye–Hückel’s dependence. Our analysis indicates that enHD’s response to ionic strength is mediated by a built-in electrostatic spring-loaded latch that operates as a conformational transducer. We also found that, at moderate ionic strengths, enHD changes conformation upon binding to cognate DNA. This change is of larger amplitude and somewhat orthogonal to the response to ionic strength. As a consequence, very high ionic strengths (e.g., 700 mM) block the electrostatic-spring-loaded latch and binding to cognate DNA becomes lock-and-key. However, the interplay between enHD conformation and cognate DNA binding is robust across a range of ionic strengths (i.e., 45 to 300 mM) that covers the physiologically-relevant conditions. Therefore, our results demonstrate the presence of a mechanism for the conformational control of cognate DNA recognition on a eukaryotic DBD. This mechanism can function as a signal transducer that locks the DBD in place upon encountering the cognate site during active DNA scanning. The electrostatic-spring-loaded latch of enHD can also enable the fine control of DNA recognition in response to transient changes in local ionic strength induced by variate physiological processes.

The alternative protein isoform NK2B, encoded by the vnd/NK-2 proneural gene, directly activates transcription and is expressed following the start of cells differentiation

NK-2 is a homeodomain protein essential for the development of the central nervous system in the Drosophila embryo. Here, we show that the vnd/NK-2 gene encodes an additional protein isoform (NK-2B) that differs from the known one (NK-2A) in its N-terminal domain. While NK-2A is a transcription repressor, NK-2B directly activates transcription from promoters containing NK-2 binding sites, with its N-terminal domain possessing a strong transcription activation potency. The transcription of NK-2B starts at the onset of metamorphosis. Its expression is observed in precursors of differentiating photoreceptors and in photoreceptors of the adult eye. Both NK-2B and NK-2A are expressed in the lamina. However, the expression of NK-2A is mostly associated with the undifferentiated state of nervous cells.

Kinetic analysis of Drosophila Vnd protein containing homeodomain with its target sequence

Homeodomain (HD) is a highly conserved DNA-binding domain composed of helix-turn-helix motif. Drosophila Vnd (Ventral nervous system defective) containing HD acts as a regulator to either enhance or suppress gene expression upon binding to its target sequence. In this study, kinetic analysis of Vnd binding to DNA was performed. The result demonstrates that DNA-binding affinity of the recombinant protein containing HD and NK2-specific domain (NK2-SD) was higher than that of the full-length Vnd. To access whether phosphorylation sites within HD and NK2-SD affect the interaction of the protein with the target sequence, alanine substitutions were introduced. The result shows that S631A mutation within NK2-SD does not contribute significantly to the DNA-binding affinity. However, S571A and T600A mutations within HD showed lower affinity for DNA binding. In addition, DNA-binding analysis using embryonic nuclear protein also demonstrates that Vnd interacts with other nuclear proteins, suggesting the existence of Vnd as a complex.

Interactions of the acidic domain and SRF interacting motifs with the NKX3.1 homeodomain

NKX3.1 is a prostate tumor suppressor belonging to the NK-2 family of homeodomain (HD) transcription factors. NK-2 family members often possess a stretch of 10-15 residues enriched in acidic amino acids, the acidic domain (AD), in the flexible, disordered region N-terminal to the HD. Interactions between the N-terminal region of NKX3.1 and its homeodomain affect protein stability and DNA binding. CD spectroscopy measuring the thermal unfolding of NKX3.1 constructs showed a 2 degrees C intramolecular stabilization of the HD by the N-terminal region containing the acidic domain (residues 85-96). CD of mixtures of various N-terminal peptides with a construct containing just the HD showed that the acidic domain and the following region, the SRF interacting (SI) motif (residues 99-105), was necessary for this stabilization. Phosphorylation of the acidic domain is known to slow proteasomal degradation of NKX3.1 in prostate cells, and NMR spectroscopy was used to measure and map the interaction of the HD with phosphorylated and nonphosphorylated forms of the AD peptide. The interaction with the phosphorylated AD peptide was considerably stronger (K(d) = 0.5 +/- 0.2 mM), resulting in large chemical shift perturbations for residues Ser150 and Arg175 in the HD, as well as a 2 degrees C increase in the HD thermal stability compared to that of the nonphosphorylated form. NKX3.1 constructs with AD phosphorylation site threonine residues (89 and 93) mutated to glutamate were 4 degrees C more stable than HD alone. Using polymer theory, effective concentrations for interactions between domains connected by flexible linkers are predicted to be in the millimolar range, and thus, the weak intramolecular interactions observed here could conceivably modulate or compete with stronger, intermolecular interactions with the NKX3.1 HD.

How the Dorsal gradient works: insights from postgenome technologies

Gradients of extracellular signaling molecules and transcription factors are used in a variety of developmental processes, including the patterning of the Drosophila embryo, the establishment of diverse neuronal cell types in the vertebrate neural tube, and the anterior-posterior patterning of vertebrate limbs. Here, we discuss how a gradient of the maternal transcription factor Dorsal produces complex patterns of gene expression across the dorsal-ventral (DV) axis of the early Drosophila embryo. The identification of 60-70 Dorsal target genes, along with the characterization of approximately 35 associated regulatory DNAs, suggests that there are at least six different regulatory codes driving diverse DV expression profiles.

NOBOX homeobox mutation causes premature ovarian failure

NOBOX (newborn ovary homeobox gene) is an oocyte-specific homeobox gene that plays a critical role in early folliculogenesis and represents a candidate gene for nonsyndromic ovarian failure. We investigated whether mutations in the NOBOX gene cause premature ovarian failure (POF). We sequenced the NOBOX gene in 96 white women with POF and discovered seven known single-nucleotide polymorphisms and four novel variations, two of which, p.Arg355His and p.Arg360Gln, cause missense mutations in the homeobox domain. Electrophoretic mobility shift assay (EMSA) confirmed that the missense mutation, p.Arg355His, disrupted NOBOX homeodomain binding to NOBOX DNA-binding element (NBE) and had a dominant negative effect on the binding of wild-type NOBOX to DNA. Our findings demonstrate that NOBOX mutations can cause POF.

The Nk-2 box of the Drosophila homeodomain protein, Vnd, contributes to its repression activity in a Groucho-dependent manner

The transcription factor, Vnd, is a dual regulator that specifies ventral neuroblast identity in Drosophila by both repressing and activating target genes. Vnd and its homologues have a conserved amino acid sequence, the Nk-2 box or Nk specific domain, as well a conserved DNA-binding homeodomain and an EhI-type Groucho interaction domain. However, the function of the conserved Nk-2 box has not been fully defined. To explore its function, we deleted the Nk-2 box and compared the regulatory activity of mutant Vnd in transgenic over-expression assays to that of the wild-type protein. We were unable to assign regulatory activity to the Nk-2 box using an over-expression assay, because the mutant protein activated expression of endogenous Vnd, masking a requirement for the Nk-2 box. However, in transgenic rescue assays, Vnd lacking the Nk-2 box repressed ind expression at 30% lower levels than the wild-type protein. Moreover, in transient transfection assays using Gal4 DNA-binding domain-Vnd chimeras, the repression activity of Vnd lacking the Nk-2 box was compromised. Because Vnd represses target gene expression in conjunction with Groucho, we asked whether the Nk-2 box affects Vnd's ability to interact with this co-repressor. Vnd lacking the Nk-2 box binds Groucho 30% less efficiently than wild-type Vnd in co-immunoprecipitations. These data suggest that the Nk-2 box contributes to the repression activity of Vnd by stabilizing its interaction with the co-repressor, Groucho.

Physical and functional interactions between the prostate suppressor homeoprotein NKX3.1 and serum response factor

The NKX3.1 transcription factor is an NK family homeodomain protein and a tumor suppressor gene that is haploinsufficient and down-regulated in the early phases of prostate cancer. Like its cardiac homolog, NKX2.5, NKX3.1 acts synergistically with serum response factor (SRF) to activate expression from the smooth muscle gamma-actin (SMGA) gene promoter. Using NMR spectroscopy, three conserved motifs in a construct containing the N-terminal region and homeodomain of NKX3.1 were observed to interact with the MADS box domain of SRF. These motifs interacted both in the absence of DNA and when both proteins were bound to a SMGA promoter DNA sequence. No significant interaction was seen between the homeodomain and SRF MADS box. One of the SRF-interacting regions was the tinman (TN) or engrailed homology-1 motif (EH-1), residues 29-35 (FLIQDIL), which for other NK proteins is the site of interaction with the repressor protein Groucho. A second hydrophobic interacting region was designated the SRF-interacting (SI) motif and included residues 99-105 (LGSYLLD). A third interacting motif was the acidic region adjacent to the SI motif including residues 88-96 (ETLAETEPE). The acidic domain (AD) motif signals also showed strengthening upon the NKX3.1 homeodomain binding to DNA in the absence of SRF, consistent with the acidic region weakly interacting with the homeodomain in the unbound state. The importance of these linear motifs in the transcriptional interaction of NKX3.1 and SRF was demonstrated by targeted mutagenesis of an NKX3.1 expression vector in a SMGA reporter assay. The results implicate the NKX3.1 N-terminal region in regulation of transcriptional activity of this tumor suppressor.

An eh1-like motif in odd-skipped mediates recruitment of Groucho and repression in vivo

Drosophila Groucho, like its vertebrate Transducin-like Enhancer-of-split homologues, is a corepressor that silences gene expression in numerous developmental settings. Groucho itself does not bind DNA but is recruited to target promoters by associating with a large number of DNA-binding negative transcriptional regulators. These repressors tether Groucho via short conserved polypeptide sequences, of which two have been defined. First, WRPW and related tetrapeptide motifs have been well characterized in several repressors. Second, a motif termed Engrailed homology 1 (eh1) has been found predominantly in homeodomain-containing transcription factors. Here we describe a yeast two-hybrid screen that uncovered physical interactions between Groucho and transcription factors, containing eh1 motifs, with different types of DNA-binding domains. We show that one of these, the zinc finger protein Odd-skipped, requires its eh1-like sequence for repressing specific target genes in segmentation. Comparison between diverse eh1 motifs reveals a bias for the phosphoacceptor amino acids serine and threonine at a fixed position, and a mutational analysis of Odd-skipped indicates that these residues are critical for efficient interactions with Groucho and for repression in vivo. Our data suggest that phosphorylation of these phosphomeric residues, if it occurs, will down-regulate Groucho binding and therefore repression, providing a mechanism for posttranslational control of Groucho-mediated repression.

Automatic discovery of cross-family sequence features associated with protein function

Background

Methods for predicting protein function directly from amino acid sequences are useful tools in the study of uncharacterised protein families and in comparative genomics. Until now, this problem has been approached using machine learning techniques that attempt to predict membership, or otherwise, to predefined functional categories or subcellular locations. A potential drawback of this approach is that the human-designated functional classes may not accurately reflect the underlying biology, and consequently important sequence-to-function relationships may be missed.

Results

We show that a self-supervised data mining approach is able to find relationships between sequence features and functional annotations. No preconceived ideas about functional categories are required, and the training data is simply a set of protein sequences and their UniProt/Swiss-Prot annotations. The main technical aspect of the approach is the co-evolution of amino acid-based regular expressions and keyword-based logical expressions with genetic programming. Our experiments on a strictly non-redundant set of eukaryotic proteins reveal that the strongest and most easily detected sequence-to-function relationships are concerned with targeting to various cellular compartments, which is an area already well studied both experimentally and computationally. Of more interest are a number of broad functional roles which can also be correlated with sequence features. These include inhibition, biosynthesis, transcription and defence against bacteria. Despite substantial overlaps between these functions and their corresponding cellular compartments, we find clear differences in the sequence motifs used to predict some of these functions. For example, the presence of polyglutamine repeats appears to be linked more strongly to the "transcription" function than to the general "nuclear" function/location.

Conclusion

We have developed a novel and useful approach for knowledge discovery in annotated sequence data. The technique is able to identify functionally important sequence features and does not require expert knowledge. By viewing protein function from a sequence perspective, the approach is also suitable for discovering unexpected links between biological processes, such as the recently discovered role of ubiquitination in transcription.

Hox genes and their candidate downstream targets in the developing central nervous system

1. Homeobox (Hox) genes were originally discovered in the fruit fly Drosophila, where they function through a conserved homeodomain as transcriptional regulators to control embryonic morphogenesis. Since then over 1000 homeodomain proteins have been identified in several species. In vertebrates, 39 Hox genes have been identified as homologs of the original Drosophila complex, and like their Drosophila counterparts they are organized within chromosomal clusters. Vertebrate Hox genes have also been shown to play a critical role in embryonic development as transcriptional regulators. 2. Both the Drosophila and vertebrate Hox genes have been shown to interact with various cofactors, such as the TALE homeodomain proteins, in recognition of consensus sequences within regulatory elements of their target genes. These protein-protein interactions are believed to contribute to enhancing the specificity of target gene recognition in a cell-type or tissue- dependent manner. The regulatory activity of a particular Hox protein on a specific regulatory element is highly variable and dependent on its interacting partners within the transcriptional complex. 3. In vertebrates, Hox genes display spatially restricted patterns of expression within the developing CNS, both along the anterioposterior and dorsoventral axis of the embryo. Their restricted gene expression is suggestive of a regulatory role in patterning of the CNS, as well as in cell specification. Determining the precise function of individual Hox genes in CNS morphogenesis through classical mutational analyses is complicated due to functional redundancy between Hox genes. 4. Understanding the precise mechanisms through which Hox genes mediate embryonic morphogenesis requires the identification of their downstream target genes. Although Hox genes have been implicated in the regulation of several pathways, few target genes have been shown to be under their direct regulatory control. Development of methodologies used for the isolation of target genes and for the analysis of putative targets will be beneficial in establishing the genetic pathways controlled by Hox factors. 5. Within the developing CNS various cell adhesion molecules and signaling molecules have been identified as candidate downstream target genes of Hox proteins. These targets play a role in processes such as cell migration and differentiation, and are implicated in contributing to neuronal processes such as plasticity and/or specification. Hence, Hox genes not only play a role in patterning of the CNS during early development, but may also contribute to cell specification and identity.

Identification and analysis of vnd/NK-2 homeodomain binding sites in genomic DNA

Vnd/NK-2 homeodomain affinity column chromatography was used to purify Drosophila DNA fragments bound by the vnd/NK-2 homeodomain. Sequencing the selected genomic DNA fragments led to the identification of 77 Drosophila DNA fragments that were grouped into 42 vnd/NK-2 homeodomain-binding loci. Most loci were within upstream or intronic regions, especially first introns. Nineteen of the Drosophila DNA fragments cloned correspond to one locus, termed Clone A, which is 312 bp in length and contains five vnd/NK-2 homeodomain core consensus binding sites, 5'-AAGTG, and is part of the first intron of the Beadex gene. We further analyzed the interactions between Clone A and vnd/NK-2 homeodomain protein by mobility-shift assay, DNase I footprinting, methylation interference, and ethylation interference. The DNase I footprinting analysis of Clone A with vnd/NK-2 homeodomain protein revealed three strong binding sites and one weak binding site between 15 and 130 bp of Clone A. We also analyzed binding of the vnd/NK-2 homeodomain to the 5'-flanking sequence of vnd/NK-2 genomic DNA. The DNase I footprinting result showed that there are two strong binding sites and five weak binding sites in the fragment between -385 and -675 bp from the transcription start site of the vnd/NK-2 gene.

Contextual interactions determine whether the Drosophila homeodomain protein, Vnd, acts as a repressor or activator

At the molecular level, members of the NKx2.2 family of transcription factors establish neural compartment boundaries by repressing the expression of homeobox genes specific for adjacent domains [Muhr et al. (2001) Cell, 104, 861-873; Weiss et al. (1998) Genes Dev., 12, 3591-3602]. The Drosophila homologue, vnd, interacts genetically with the high-mobility group protein, Dichaete, in a manner suggesting co-operative activation [Zhao and Skeath (2002) Development, 129, 1165-1174]. However, evidence for direct interactions and transcriptional activation is lacking. Here, we present molecular evidence for the interaction of Vnd and Dichaete that leads to the activation of target gene expression. Two-hybrid interaction assays indicate that Dichaete binds the Vnd homeodomain, and additional Vnd sequences stabilize this interaction. In addition, Vnd has two activation domains that are typically masked in the intact protein. Whether vnd can activate or repress transcription is context-dependent. Full-length Vnd, when expressed as a Gal4 fusion protein, acts as a repressor containing multiple repression domains. A divergent domain in the N-terminus, not found in vertebrate Vnd-like proteins, causes the strongest repression. The co-repressor, Groucho, enhances Vnd repression, and these two proteins physically interact. The data presented indicate that the activation and repression domains of Vnd are complex, and whether Vnd functions as a transcriptional repressor or activator depends on both intra- and inter-molecular interactions.

NMR structure of the noncytotoxic alpha-sarcin mutant Delta(7-22): the importance of the native conformation of peripheral loops for activity

The deletion mutant Delta(7-22) of alpha-sarcin, unlike its wild-type protein counterpart, lacks the specific ability to degrade rRNA in intact ribosomes and exhibits an increased unspecific ribonuclease activity and decreased interaction with lipid vesicles. In trying to shed light on these differences, we report here on the three-dimensional structure of the Delta(7-22) alpha-sarcin mutant using NMR methods. We also evaluated its dynamic properties on the basis of theoretical models and measured its correlation time (6.2 nsec) by time-resolved fluorescence anisotropy. The global fold characteristic of ribotoxins is preserved in the mutant. The most significant differences with respect to the alpha-sarcin structure are concentrated in (1) loop 2, (2) loop 3, which adopts a new orientation, and (3) loop 5, which shows multiple conformations and an altered dynamics. The interactions between loop 5 and the N-terminal hairpin are lost in the mutant, producing increased solvent accessibility of the active-site residues. The degree of solvent exposure of the catalytic His 137 is similar to that shown by His 92 in RNase T1. Additionally, the calculated order parameters of residues belonging to loop 5 in the mutant correspond to an internal dynamic behavior more similar to RNase T1 than alpha-sarcin. On the other hand, changes in the relative orientation of loop 3 move the lysine-rich region 111-114, crucial for substrate recognition, away from the active site. All of the structural and dynamic data presented here reveal that the mutant is a hybrid of ribotoxins and noncytotoxic ribonucleases, consistent with its biological properties.

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.