Sequence and Structural Determinants of Specificity Differences between Paralogous Transcription Factors

Transcription factors (TFs) control the temporal and spatial expression of target genes by interacting with DNA in a sequence-specific manner. Recent advances in high throughput experiments that measure TF-DNA interactions in vitro and in vivo have facilitated the identification of DNA binding sites for thousands of TFs. However, it remains unclear how each individual TF achieves its specificity, especially in the case of paralogous TFs that recognize distinct target genomic sites despite sharing very similar DNA binding motifs. In my work, I used a combination of high throughput in vitro protein-DNA binding assays and machine-learning algorithms to characterize and model the binding specificity of 11 paralogous TFs from 4 distinct structural families. My work proves that even very closely related paralogous TFs, with indistinguishable DNA binding motifs, oftentimes exhibit differential binding specificity for their genomic target sites, especially for sites with moderate binding affinity. Importantly, the differences I identify in vitro and through computational modeling help explain, at least in part, the differential in vivo genomic targeting by paralogous TFs. Future work will focus on in vivo factors that might also be important for specificity differences between paralogous TFs, such as DNA methylation, interactions with protein cofactors, or the chromatin environment. In this larger context, my work emphasizes the importance of intrinsic DNA binding specificity in targeting of paralogous TFs to the genome.

Genomic Analysis of Pseudomonas putida Phage tf with Localized Single-Strand DNA Interruptions

The complete sequence of the 46,267 bp genome of the lytic bacteriophage tf specific to Pseudomonas putida PpG1 has been determined. The phage genome has two sets of convergently transcribed genes and 186 bp long direct terminal repeats. The overall genomic architecture of the tf phage is similar to that of the previously described Pseudomonas aeruginosa phages PaP3, LUZ24 and phiMR299-2, and 39 out of the 72 products of predicted tf open reading frames have orthologs in these phages. Accordingly, tf was classified as belonging to the LUZ24-like bacteriophage group. However, taking into account very low homology levels between tf DNA and that of the other phages, tf should be considered as an evolutionary divergent member of the group. Two distinguishing features not reported for other members of the group were found in the tf genome. Firstly, a unique end structure - a blunt right end and a 4-nucleotide 3'-protruding left end - was observed. Secondly, 14 single-chain interruptions (nicks) were found in the top strand of the tf DNA. All nicks were mapped within a consensus sequence 5'-TACT/RTGMC-3'. Two nicks were analyzed in detail and were shown to be present in more than 90% of the phage population. Although localized nicks were previously found only in the DNA of T5-like and phiKMV-like phages, it seems increasingly likely that this enigmatic structural feature is common to various other bacteriophages.

Rapid evolution of the DNA-binding site in LAGLIDADG homing endonucleases

Sequence analysis of chloroplast and mitochondrial large subunit rRNA genes from over 75 green algae disclosed 28 new group I intron-encoded proteins carrying a single LAGLIDADG motif. These putative homing endonucleases form four subfamilies of homologous enzymes, with the members of each subfamily being encoded by introns sharing the same insertion site. We showed that four divergent endonucleases from the I-CreI subfamily cleave the same DNA substrates. Mapping of the 66 amino acids that are conserved among the members of this subfamily on the 3-dimensional structure of I-CreI bound to its recognition sequence revealed that these residues participate in protein folding, homodimerization, DNA recognition and catalysis. Surprisingly, only seven of the 21 I-CreI amino acids interacting with DNA are conserved, suggesting that I-CreI and its homologs use different subsets of residues to recognize the same DNA sequence. Our sequence comparison of all 45 single-LAGLIDADG proteins identified so far suggests that these proteins share related structures and that there is a weak pressure in each subfamily to maintain identical protein–DNA contacts. The high sequence variability we observed in the DNA-binding site of homologous LAGLIDADG endonucleases provides insight into how these proteins evolve new DNA specificity.

A new plasmid display technology for the in vitro selection of functional phenotype–genotype linked proteins

Background: Display technologies which allow peptides or proteins to be physically associated with the encoding DNA are central to procedures which involve screening of protein libraries in vitro for new or altered function. Here we describe a new system designed specifically for the display of libraries of diverse, functional proteins which utilises the DNA binding protein nuclear factor κB (NF-κB) p50 to establish a phenotype–genotype link between the displayed protein and the encoding gene. Results: A range of model fusion proteins to either the amino- or carboxy-terminus of NF-κB p50 have been constructed and shown to retain the picomolar affinity and DNA specificity of wild-type NF-κB p50. Through use of an optimal combination of binding buffer and DNA target sequence, the half-life of p50–DNA complexes could be increased to over 47 h, enabling the competitive selection of a variety of protein–plasmid complexes with enrichment factors of up to 6000-fold per round. The p50-based plasmid display system was used to enrich a maltose binding protein complex to homogeneity in only three rounds from a binary mixture with a starting ratio of 1:108 and to enrich to near homogeneity a single functional protein from a phenotype–genotype linked Escherichia coli genomic library using in vitro functional selections. Conclusions: A new display technology is described which addresses the challenge of functional protein display. The results demonstrate that plasmid display is sufficiently sensitive to select a functional protein from large libraries and that it therefore represents a useful addition to the repertoire of display technologies.

A new plasmid display technology for the in vitro selection of functional phenotype-genotype linked proteins

Background Display technologies which allow peptides or proteins to be physically associated with the encoding DNA are central to procedures which involve screening of protein libraries in vitro for new or altered function. Here we describe a new system designed specifically for the display of libraries of diverse, functional proteins which utilises the DNA binding protein nuclear factor κB (NF-κB) p50 to establish a phenotype-genotype link between the displayed protein and the encoding gene. Results A range of model fusion proteins to either the amino- or carboxy-terminus of NF-κB p50 have been constructed and shown to retain the picomolar affinity and DNA specificity of wild-type NF-κB p50. Through use of an optimal combination of binding buffer and DNA target sequence, the half-life of p50-DNA complexes could be increased to over 47 h, enabling the competitive selection of a variety of protein-plasmid complexes with enrichment factors of up to 6000-fold per round. The p50-based plasmid display system was used to enrich a maltose binding protein complex to homogeneity in only three rounds from a binary mixture with a starting ratio of 1:108 and to enrich to near homogeneity a single functional protein from a phenotype-genotype linked Escherichia coli genomic library using in vitro functional selections. Conclusions A new display technology is described which addresses the challenge of functional protein display. The results demonstrate that plasmid display is sufficiently sensitive to select a functional protein from large libraries and that it therefore represents a useful addition to the repertoire of display technologies.

Explicit DNase sequence bias modeling enables high-resolution transcription factor footprint detection.

DNaseI footprinting is an established assay for identifying transcription factor (TF)-DNA interactions with single base pair resolution. High-throughput DNase-seq assays have recently been used to detect in vivo DNase footprints across the genome. Multiple computational approaches have been developed to identify DNase-seq footprints as predictors of TF binding. However, recent studies have pointed to a substantial cleavage bias of DNase and its negative impact on predictive performance of footprinting. To assess the potential for using DNase-seq to identify individual binding sites, we performed DNase-seq on deproteinized genomic DNA and determined sequence cleavage bias. This allowed us to build bias corrected and TF-specific footprint models. The predictive performance of these models demonstrated that predicted footprints corresponded to high-confidence TF-DNA interactions. DNase-seq footprints were absent under a fraction of ChIP-seq peaks, which we show to be indicative of weaker binding, indirect TF-DNA interactions or possible ChIP artifacts. The modeling approach was also able to detect variation in the consensus motifs that TFs bind to. Finally, cell type specific footprints were detected within DNase hypersensitive sites that are present in multiple cell types, further supporting that footprints can identify changes in TF binding that are not detectable using other strategies.

The salt dependence of DNA recognition by NF-kappaB p50 : a detailed kinetic analysis of the effects on affinityand specificity

The binding kinetics of NF-kappaB p50 to the Ig-kappaB site and to a DNA duplex with no specific binding site were determined under varying conditions of potassium chloride concentration using a surface plasmonresonance biosensor. Association and dissociation rate constants were measured enabling calculation of the dissociation constants. Under previously established high affinity buffer conditions, the k a for both sequences was in the order of 10(7) M-1s-1whilst the k d values varied 600-fold in a sequence-dependent manner between 10(-1) and 10(-4 )s-1, suggesting that the selectivity of p50 for different sequences is mediated primarily through sequence-dependent dissociation rates. The calculated K D value for the Ig-kappaB sequence was 16 pM, whilst the K D for the non-specific sequence was 9.9 nM. As the ionic strength increased to levels which are closer to that of the cellular environment, the binding of p50 to the non-specific sequence was abolished whilst the specific affinity dropped to nanomolar levels. From these results, a mechanism is proposed in which p50 binds specific sequences with high affinity whilst binding non-specific sequences weakly enough to allow efficient searching of the DNA.

Substrate specificity plays an important role in uncoupling the catalytic and scaffolding activities of rat testis DNA topoisomerase IIalpha

Topoisomerase II (topo II) is a dyadic enzyme found in all eukaryotic cells. Topo II is involved in a number of cellular processes related to DNA metabolism, including DNA replication, recombination and the maintenance of genomic stability. We discovered a correlation between the development of postnatal testis and increased binding of topo IIalpha to the chromatin fraction. We used this observation to characterize DNA-binding specificity and catalytic properties of purified testis topo IIalpha. The results indicate that topo IIalpha binds a substrate containing the preferred site with greater affinity and, consequently, catalyzes the conversion of form I to form IV DNA more efficiently in contrast to substrates lacking such a site. Interestingly, topo IIalpha displayed high-affinity and cooperativity in binding to the scaffold associated region. In contrast to the preferred site, however, high-affinity binding of topo IIalpha to the scaffold-associated region failed to result in enhanced catalytic activity. Intriguingly, competition assays involving scaffold-associated region revealed an additional DNA-binding site within the dyadic topo IIalpha. These results implicate a dual role for topo IIalpha in vivo consistent with the notion that its sequestration to the chromatin might play a role in chromosome condensation and decondensation during spermatogenesis.

DNA Clasping by Mycobacterial HU: The C-Terminal Region of HupB Mediates Increased Specificity of DNA Binding

Background: HU a small, basic, histone like protein is a major component of the bacterial nucleoid. E. coli has two subunits of HU coded by hupA and hupB genes whereas Mycobacterium tuberculosis (Mtb) has only one subunit of HU coded by ORF Rv2986c (hupB gene). One noticeable feature regarding Mtb HupB, based on sequence alignment of HU orthologs from different bacteria, was that HupB(Mtb) bears at its C-terminal end, a highly basic extension and this prompted an examination of its role in Mtb HupB function. Methodology/Principal Findings: With this objective two clones of Mtb HupB were generated; one expressing full length HupB protein (HupB(Mtb)) and another which expresses only the N terminal region (first 95 amino acid) of hupB (HupB(MtbN)). Gel retardation assays revealed that HupBMtbN is almost like E. coli HU (heat stable nucleoid protein) in terms of its DNA binding, with a binding constant (K-d) for linear dsDNA greater than 1000 nM, a value comparable to that obtained for the HU alpha alpha and HU alpha beta forms. However CTR (C-terminal Region) of HupB(Mtb) imparts greater specificity in DNA binding. HupB(Mtb) protein binds more strongly to supercoiled plasmid DNA than to linear DNA, also this binding is very stable as it provides DNase I protection even up to 5 minutes. Similar results were obtained when the abilities of both proteins to mediate protection against DNA strand cleavage by hydroxyl radicals generated by the Fenton's reaction, were compared. It was also observed that both the proteins have DNA binding preference for A: T rich DNA which may occur at the regulatory regions of ORFs and the oriC region of Mtb. Conclusions/Significance: These data thus point that HupB(Mtb) may participate in chromosome organization in-vivo, it may also play a passive, possibly an architectural role.

Sequence specificity of Z-DNA formation in oligonucleotides

The sequence specific requirement for B----Z transition in solution was examined in d(CGTGCGCACG), d(CGTACGTACG), d(ACGTACGT) in presence of various Z-inducing factors. Conformational studies show that inspite of the alternating nature of purines and pyrimidines, the aforementioned sequences do not undergo B----Z transition under the influence of NaCl, hexamine cobalt chloride and ethanol. A comparison with the crystal structures of an assorted array of purine and pyrimidine sequences show that the sequence requirement for B----Z transition is much more stringent in solution as compared to the solid state. The disruptive influence of AT base pairs in B to Z transition is discussed.

Mutations in hpyAVIBM, C5 cytosine DNA methyltransferase from Helicobacter pylori result in relaxed specificity

The genome of Helicobacter pylori is rich in restrictionmodification (RM) systems. Approximately 4% of the genome codes for components of RM systems. hpyAVIBM, which codes for a phase-variable C5 cytosine methyltransferase (MTase) from H. pylori, lacks a cognate restriction enzyme. Over-expression of M.HpyAVIB in Escherichia coli enhances the rate of mutations. However, when the catalytically inactive F9N or C82W mutants of M.HpyAVIB were expressed in E. coli, mutations were not observed. The M.HpyAVIB gene itself was mutated to give rise to different variants of the MTase. M.HpyAVIB variants were purified and differences in kinetic properties and specificity were observed. Intriguingly, purified MTase variants showed relaxed substrate specificity. Homologues of hpyAVIBM homologues amplified and sequenced from different clinical isolates showed similar variations in sequence. Thus, hpyAVIBM presents an interesting example of allelic variations in H. pylori where changes in the nucleotide sequence result in proteins with new properties.

The zinc finger proteins Mxr1p and repressor of phosphoenolpyruvate carboxykinase (ROP) have the same DNA binding specificity but regulate Methanol Metabolism antagonistically in pichia pastoris

The methanol-inducible alcohol oxidase I (AOXI) promoter of the methylotrophic yeast, Pichia pastoris, is used widely for the production of recombinant proteins. AOXI transcription is regulated by the zinc finger protein Mxr1p (methanol expression regulator 1). ROP (repressor of phosphoenolpyruvate carboxykinase, PEPCK) is a methanol- and biotin starvation-inducible zinc finger protein that acts as a negative regulator of PEPCK in P. pastoris cultured in biotin-deficient, glucose-ammonium medium. The function of ROP during methanol metabolism is not known. In this study, we demonstrate that ROP represses methanol-inducible expression of AOXI when P. pastoris is cultured in a nutrient-rich medium containing yeast extract, peptone, and methanol (YPM). Deletion of the gene encoding ROP results in enhanced expression of AOXI and growth promotion whereas overexpression of ROP results in repression of AOXI and growth retardation of P. pastoris cultured in YPM medium. Surprisingly, deletion or overexpression of ROP has no effect on AOXI gene expression and growth of P. pastoris cultured in a minimal medium containing yeast nitrogen base and methanol (YNBM). Subcellular localization studies indicate that ROP translocates from cytosol to nucleus of cells cultured in YPM but not YNBM. In vitro DNA binding studies indicate that AOXI promoter sequences containing 5' CYCCNY 3' motifs serve as binding sites for Mxr1p as well as ROP. Thus, Mxr1p and ROP exhibit the same DNA binding specificity but regulate methanol metabolism antagonistically in P. pastoris. This is the first report on the identification of a transcriptional repressor of methanol metabolism in any yeast species.

Ca2+ Binding to the ExDxD Motif Regulates the DNA Cleavage Specificity of a Promiscuous Endonuclease

Most of the restriction endonucleases (REases) are dependent on Mg2+ for DNA cleavage, and in general, Ca2+ inhibits their activity. RKpnI, an HNH active site containing beta beta alpha-Me finger nuclease, is an exception. In presence of Ca2+, the enzyme exhibits high-fidelity DNA cleavage and complete suppression of Mg2+-induced promiscuous activity. To elucidate the mechanism of unusual Ca2+-mediated activity, we generated alanine variants in the putative Ca-2+ binding motif, E(132)xD(134)xD(136), of the enzyme. Mutants showed decreased levels of DNA cleavage in the presence of Ca2+. We demonstrate that ExDxD residues are involved in Ca2+ coordination; however, the invariant His of the catalytic HNH motif acts as a general base for nucleophile activation, and the other two active site residues, D148 and Q175, also participate in Ca2+-mediated cleavage. Insertion of a 10-amino acid linker to disrupt the spatial organization of the ExDxD and HNH motifs impairs Ca2+ binding and affects DNA cleavage by the enzyme. Although ExDxD mutant enzymes retained efficient cleavage at the canonical sites in the presence of Mg2+, the promiscuous activity was greatly reduced, indicating that the carboxyl residues of the acidic triad play an important role in sequence recognition by the enzyme. Thus, the distinct Ca2+ binding motif that confers site specific cleavage upon Ca2+ binding is also critical for the promiscuous activity of the Mg2+-bound enzyme, revealing its role in metal ion-mediated modulation of DNA cleavage.

Specificity of the ModA11, ModA12 and ModD1 epigenetic regulator N-6-adenine DNA methyltransferases of Neisseria meningitidis

Phase variation (random ON/OFF switching) of gene expression is a common feature of host-adapted pathogenic bacteria. Phase variably expressed N-6-adenine DNA methyltransferases (Mod) alter global methylation patterns resulting in changes in gene expression. These systems constitute phase variable regulons called phasevarions. Neisseria meningitidis phasevarions regulate genes including virulence factors and vaccine candidates, and alter phenotypes including antibiotic resistance. The target site recognized by these Type III N-6-adenine DNA methyltransferases is not known. Single molecule, real-time (SMRT) methylome analysis was used to identify the recognition site for three key N. meningitidis methyltransferases: ModA11 (exemplified by M.NmeMC58I) (5'-CGY(m6)AG-3'), ModA12 (exemplified by M.Nme77I, M.Nme18I and M.Nme579II) (5'-AC(m6)ACC-3') and ModD1 (exemplified by M.Nme579I) (5'-CC(m6)AGC-3'). Restriction inhibition assays and mutagenesis confirmed the SMRT methylome analysis. The ModA11 site is complex and atypical and is dependent on the type of pyrimidine at the central position, in combination with the bases flanking the core recognition sequence 5'-CGY(m6)AG-3'. The observed efficiency of methylation in the modA11 strain (MC58) genome ranged from 4.6% at 5'-GCGC(m6)AGG-3' sites, to 100% at 5'-ACGT(m6)AGG-3' sites. Analysis of the distribution of modified sites in the respective genomes shows many cases of association with intergenic regions of genes with altered expression due to phasevarion switching.