Altering the DNA-binding specificity of the yeast Matalpha 2 homeodomain protein

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.

Determination and inference of eukaryotic transcription factor sequence specificity

Transcription factor (TF) DNA sequence preferences direct their regulatory activity, but are currently known for only ∼1% of eukaryotic TFs. Broadly sampling DNA-binding domain (DBD) types from multiple eukaryotic clades, we determined DNA sequence preferences for >1,000 TFs encompassing 54 different DBD classes from 131 diverse eukaryotes. We find that closely related DBDs almost always have very similar DNA sequence preferences, enabling inference of motifs for ∼34% of the ∼170,000 known or predicted eukaryotic TFs. Sequences matching both measured and inferred motifs are enriched in chromatin immunoprecipitation sequencing (ChIP-seq) peaks and upstream of transcription start sites in diverse eukaryotic lineages. SNPs defining expression quantitative trait loci in Arabidopsis promoters are also enriched for predicted TF binding sites. Importantly, our motif "library" can be used to identify specific TFs whose binding may be altered by human disease risk alleles. These data present a powerful resource for mapping transcriptional networks across eukaryotes.

Changing the recognition site of a conjugative relaxase by rational design

TrwC is a relaxase protein, which starts and finishes DNA processing during bacterial conjugation in plasmid R388. TrwC recognizes a specific sequence of DNA (25 nucleotides) in the donor cell: the nic-site. As a model example, a single transversion C24G in nic avoids DNA processing by TrwC. Using this simple model, our objective was to obtain a proof of principle that TrwC specificity can be changed. Several structures of DNA-TrwC complexes were used as reference to design a focused saturation mutagenesis library (NNK) randomizing amino acid Lys262, since its side chain seems to sterically hinder the recognition of the C24G nic mutation by wild-type TrwC. Using bacterial conjugation as an in vivo selection system, several TrwC variants were found that show changes in substrate specificity. These variants were also tested in a competitive assay to evaluate their conjugation efficiency.

Analysis of homeodomain specificities allows the family-wide prediction of preferred recognition sites

We describe the comprehensive characterization of homeodomain DNA-binding specificities from a metazoan genome. The analysis of all 84 independent homeodomains from D. melanogaster reveals the breadth of DNA sequences that can be specified by this recognition motif. The majority of these factors can be organized into 11 different specificity groups, where the preferred recognition sequence between these groups can differ at up to four of the six core recognition positions. Analysis of the recognition motifs within these groups led to a catalog of common specificity determinants that may cooperate or compete to define the binding site preference. With these recognition principles, a homeodomain can be reengineered to create factors where its specificity is altered at the majority of recognition positions. This resource also allows prediction of homeodomain specificities from other organisms, which is demonstrated by the prediction and analysis of human homeodomain specificities.

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.

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.

Combined analysis of expression data and transcription factor binding sites in the yeast genome

BACKGROUND

The analysis of gene expression using DNA microarrays provides genome wide profiles of the genes controlled by the presence or absence of a specific transcription factor. However, the question arises of whether a change in the level of transcription of a specific gene is caused by the transcription factor acting directly at the promoter of the gene or through regulation of other transcription factors working at the promoter.

RESULTS

To address this problem we have devised a computational method that combines microarray expression and site preference data. We have tested this approach by identifying functional targets of the a1-alpha2 complex, which represses haploid-specific genes in the yeast Saccharomyces cerevisiae. Our analysis identified many known or suspected haploid-specific genes that are direct targets of the a1-alpha2 complex, as well as a number of previously uncharacterized targets. We were also able to identify a number of haploid-specific genes which do not appear to be direct targets of the a1-alpha2 complex, as well as a1-alpha2 target sites that do not repress transcription of nearby genes. Our method has a much lower false positive rate when compared to some of the conventional bioinformatic approaches.

CONCLUSIONS

These findings show advantages of combining these two forms of data to investigate the mechanism of co-regulation of specific sets of genes.

Dominant negative dimerization of a mutant homeodomain protein in Axenfeld-Rieger syndrome

Axenfeld-Rieger syndrome is an autosomal-dominant disorder caused by mutations in the PITX2 homeodomain protein. We have studied the mechanism underlying the dominant negative K88E mutation, which occurs at position 50 of the homeodomain. By using yeast two-hybrid and in vitro pulldown assays, we have documented that PITX2a can form homodimers in the absence of DNA. Moreover, the K88E mutant had even stronger dimerization ability, primarily due to interactions involving the C-terminal region. Dimerization allowed cooperative binding of wild-type (WT) PITX2a to DNA containing tandem bicoid sites in a head-to-tail orientation (Hill coefficient, 1.73). In contrast, the WT-K88E heterodimer bound the tandem sites with greatly reduced cooperativity and decreased transactivation activity. To further explore the role of position 50 in PITX2a dimerization, we introduced a charge-conservative mutation of lysine to arginine (K88R). The K88R protein had greatly reduced binding to a TAATCC element and did not specifically bind any other TAATNN motif. Like K88E, K88R formed relatively stronger dimers with WT. As predicted by our model, the K88R protein acted in a dominant negative manner to suppress WT PITX2a activity. These results suggest that the position 50 residue in the PITX2 homeodomain plays an important role in both DNA binding and dimerization activities.

Dlx proteins position the neural plate border and determine adjacent cell fates

The lateral border of the neural plate is a major source of signals that induce primary neurons, neural crest cells and cranial placodes as well as provide patterning cues to mesodermal structures such as somites and heart. Whereas secreted BMP, FGF and Wnt proteins influence the differentiation of neural and non-neural ectoderm, we show here that members of the Dlx family of transcription factors position the border between neural and non-neural ectoderm and are required for the specification of adjacent cell fates. Inhibition of endogenous Dlx activity in Xenopus embryos with an EnR-Dlx homeodomain fusion protein expands the neural plate into non-neural ectoderm tissue whereas ectopic activation of Dlx target genes inhibits neural plate differentiation. Importantly, the stereotypic pattern of border cell fates in the adjacent ectoderm is re-established only under conditions where the expanded neural plate abuts Dlx-positive non-neural ectoderm. Experiments in which presumptive neural plate was grafted to ventral ectoderm reiterate induction of neural crest and placodal lineages and also demonstrate that Dlx activity is required in non-neural ectoderm for the production of signals needed for induction of these cells. We propose that Dlx proteins regulate intercellular signaling across the interface between neural and non-neural ectoderm that is critical for inducing and patterning adjacent cell fates.

Engineered improvements in DNA-binding function of the MATa1 homeodomain reveal structural changes involved in combinatorial control

We have engineered enhanced DNA-binding function into the a1 homeodomain by making changes in a loop distant from the DNA-binding surface. Comparison of the free and bound a1 structures suggested a mechanism linking van der Waals stacking changes in this loop to the ordering of a final turn in the DNA-binding helix of a1. Inspection of the protein sequence revealed striking differences in amino acid identity at positions 24 and 25 compared to related homeodomain proteins. These positions lie in the loop connecting helix-1 and helix-2, which is involved in heterodimerization with the alpha 2 protein. A series of single and double amino acid substitutions (a1-Q24R, a1-S25Y, a1-S25F and a1-Q24R/S25Y) were engineered, expressed and purified for biochemical and biophysical study. Calorimetric measurements and HSQC NMR spectra confirm that the engineered variants are folded and are equally or more stable than the wild-type a1 homeodomain. NMR analysis of a1-Q24R/S25Y demonstrates that the DNA recognition helix (helix-3) is extended by at least one turn as a result of the changes in the loop connecting helix-1 and helix-2. As shown by EMSA, the engineered variants bind DNA with enhanced affinity (16-fold) in the absence of the alpha 2 cofactor and the variant alpha 2/a1 heterodimers bind cognate DNA with specificity and affinity reflective of the enhanced a1 binding affinity. Importantly, in vivo assays demonstrate that the a1-Q24R/S25Y protein binds with fivefold greater affinity than wild-type a1 and is able to partially suppress defects in repression by alpha 2 mutants. As a result of these studies, we show how subtle differences in residues at a surface distant from the functional site code for a conformational switch that allows the a1 homeodomain to become active in DNA binding in association with its cofactor alpha 2.