Search

PRODORIC is worldwide one of the largest collections of prokaryotic transcription factor binding sites from multiple bacterial sources with corresponding interpretation and visualization tools. With the introduction of PRODORIC2 in 2017, the transition to a modern web interface and maintainable backend was started. With this latest PRODORIC release the database backend is now fully API-based and provides programmatical access to the complete PRODORIC data. The visualization tools Genome Browser and ProdoNet from the original PRODORIC have been reintroduced and were integrated into the PRODORIC website. Missing input and output options from the original Virtual Footprint were added again for position weight matrix pattern-based searches. The whole PRODORIC dataset was reannotated. Every transcription factor binding site was re-evaluated to increase the overall database quality. During this process, additional parameters, like bound effectors, regulation type and different types of experimental evidence have been added for every transcription factor. Additionally, 109 new transcription factors and 6 new organisms have been added. PRODORIC is publicly available at https://www.prodoric.de.

To investigate how exogenous DNA concatemerizes to form episomal artificial chromosomes (ACs), acquire equal segregation ability and maintain stable holocentromeres, we injected DNA sequences with different features, including sequences that are repetitive or complex, and sequences with different AT-contents, into the gonad of Caenorhabditis elegans to form ACs in embryos, and monitored AC mitotic segregation. We demonstrated that AT-poor sequences (26% AT-content) delayed the acquisition of segregation competency of newly formed ACs. We also co-injected fragmented Saccharomyces cerevisiae genomic DNA, differentially expressed fluorescent markers and ubiquitously expressed selectable marker to construct a less repetitive, more complex AC. We sequenced the whole genome of a strain which propagates this AC through multiple generations, and de novo assembled the AC sequences. We discovered CENP-AHCP-3 domains/peaks are distributed along the AC, as in endogenous chromosomes, suggesting a holocentric architecture. We found that CENP-AHCP-3 binds to the unexpressed marker genes and many fragmented yeast sequences, but is excluded in the yeast extremely high-AT-content centromeric and mitochondrial DNA (> 83% AT-content) on the AC. We identified A-rich motifs in CENP-AHCP-3 domains/peaks on the AC and on endogenous chromosomes, which have some similarity with each other and similarity to some non-germline transcription factor binding sites.

Quorum sensing gene expression in vibrios is regulated by the LuxR/HapR family of transcriptional factors, which includes Vibrio vulnificus SmcR. The consensus binding site of Vibrio LuxR/HapR/SmcR proteins is palindromic but highly degenerate with sequence variations at each promoter. To examine the mechanism by which SmcR recognizes diverse DNA sites, we generated SmcR separation-of-function mutants that either repress or activate transcription but not both. SmcR N55I is restricted in recognition of single base-pair variations in DNA binding site sequences and thus is defective at transcription activation but retains interaction with RNA polymerase (RNAP) alpha. SmcR S76A, L139R and N142D substitutions disrupt the interaction with RNAP alpha but retain functional DNA binding activity. X-ray crystallography and small angle X-ray scattering data show that the SmcR DNA binding domain exists in two conformations (wide and narrow), and the protein complex forms a mixture of dimers and tetramers in solution. The three RNAP interaction-deficient variants also have two DNA binding domain conformations, whereas SmcR N55I exhibits only the wide conformation. These data support a model in which two mechanisms drive SmcR transcriptional activation: interaction with RNAP and a multi-conformational DNA binding domain that permits recognition of variable DNA sites.

The glucocorticoid (GR) and androgen (AR) receptors execute unique functions in vivo, yet have nearly identical DNA binding specificities. To identify mechanisms that facilitate functional diversification among these transcription factor paralogs, we studied AR and GR in an equivalent cellular context. Analysis of chromatin and sequence features suggest that divergent binding, and corresponding gene regulation, are driven by different abilities of AR and GR to interact with relatively inaccessible chromatin. Divergent genomic binding patterns can also be the results of subtle differences in DNA binding preference between AR and GR. Furthermore, the sequence composition of large regions (>10 kb) surrounding selectively occupied binding sites differs significantly, indicating a role for the sequence environment in selectively guiding AR and GR to distinct binding sites. The comparison of binding sites that are shared between AR and GR shows that the specificity paradox can also be resolved by differences in the events that occur downstream of receptor binding. Specifically, we find that shared binding sites display receptor-specific enhancer activity, cofactor recruitment and changes in histone modifications. Genomic deletion of shared binding sites demonstrates their contribution to directing receptor-specific gene regulation. Together, these data suggest that differences in genomic occupancy as well as divergence in the events that occur downstream of receptor binding direct functional diversification among transcription factor paralogs.

Sequence-specific protein-DNA interactions are at the heart of the response of the tumor-suppressor p53 to numerous physiological and stress-related signals. Large variability has been previously reported in p53 binding to and transactivating from p53 response elements (REs) due, at least in part, to changes in direct (base) and indirect (shape) readouts of p53 REs. Here, we dissect p53 REs to decipher the mechanism by which p53 optimizes this highly regulated variable level of interaction with its DNA binding sites. We show that hemi-specific binding is more prevalent in p53 REs than previously envisioned. We reveal that sequences flanking the REs modulate p53 binding and activity and show that these effects extend to 4–5 bp from the REs. Moreover, we show here that the arrangement of p53 half-sites within its REs, relative to transcription direction, has been fine-tuned by selection pressure to optimize and regulate the response levels from p53 REs. This directionality in the REs arrangement is at least partly encoded in the structural properties of the REs. Furthermore, we show here that in the p21-5′ RE the orientation of the half-sites is such that the effect of the flanking sequences is minimized and we discuss its advantages.

We report an approach to predict DNA specificity of the tetracycline repressor (TetR) family transcription regulators (TFRs). First, a genome sequence-based method was streamlined with quantitative P-values defined to filter out reliable predictions. Then, a framework was introduced to incorporate structural data and to train a statistical energy function to score the pairing between TFR and TFR binding site (TFBS) based on sequences. The predictions benchmarked against experiments, TFBSs for 29 out of 30 TFRs were correctly predicted by either the genome sequence-based or the statistical energy-based method. Using P-values or Z-scores as indicators, we estimate that 59.6% of TFRs are covered with relatively reliable predictions by at least one of the two methods, while only 28.7% are covered by the genome sequence-based method alone. Our approach predicts a large number of new TFBs which cannot be correctly retrieved from public databases such as FootprintDB. High-throughput experimental assays suggest that the statistical energy can model the TFBSs of a significant number of TFRs reliably. Thus the energy function may be applied to explore for new TFBSs in respective genomes. It is possible to extend our approach to other transcriptional factor families with sufficient structural information.

Transcription factor decoy binding sites are short DNA sequences that can serve as “sponges” to titrate a transcription factor away from its natural binding site, therefore regulating gene expression. In this study, we harness decoy sites to develop synthetic transcription factor sponge systems to regulate gene expression for metabolic pathways in Escherichia coli. We show that transcription factor sponges can effectively regulate expression of native and heterologous genes. Tunability of the sponge can be engineered via changes in copy number or modifications to the DNA decoy site sequence. Using arginine biosynthesis as a showcase, we observe a 16-fold increase in arginine production when we introduce the sponge system to steer metabolic flux towards increased arginine biosynthesis, with negligible growth differences compared to the wild type strain. The sponge-based production strain shows high genetic stability; in contrast to a gene knock-out approach where mutations were common, we detected no mutations in the production system using the sponge-based strain. We further show that transcription factor sponges are amenable to multiplexed library screening by demonstrating enhanced tolerance to pinene with a combinatorial sponge library. Our study shows that transcription factor sponges are a powerful and compact tool for metabolic engineering.

Deregulation of microRNA (miRNA) expression plays a critical role in the transition from a physiological to a pathological state. The accurate miRNA promoter identification in multiple cell types is a fundamental endeavor towards understanding and characterizing the underlying mechanisms of both physiological as well as pathological conditions. DIANA-miRGen v4 (www.microrna.gr/mirgenv4) provides cell type specific miRNA transcription start sites (TSSs) for over 1500 miRNAs retrieved from the analysis of >1000 cap analysis of gene expression (CAGE) samples corresponding to 133 tissues, cell lines and primary cells available in FANTOM repository. MiRNA TSS locations were associated with transcription factor binding site (TFBSs) annotation, for >280 TFs, derived from analyzing the majority of ENCODE ChIP-Seq datasets. For the first time, clusters of cell types having common miRNA TSSs are characterized and provided through a user friendly interface with multiple layers of customization. DIANA-miRGen v4 significantly improves our understanding of miRNA biogenesis regulation at the transcriptional level by providing a unique integration of high-quality annotations for hundreds of cell specific miRNA promoters with experimentally derived TFBSs.

We describe an updated comprehensive database, LincSNP 3.0 (http://bioinfo.hrbmu.edu.cn/LincSNP), which aims to document and annotate disease or phenotype-associated variants in human long non-coding RNAs (lncRNAs) and circular RNAs (circRNAs) or their regulatory elements. LincSNP 3.0 has updated with several novel features, including (i) more types of variants including single nucleotide polymorphisms (SNPs), linkage disequilibrium SNPs (LD SNPs), somatic mutations and RNA editing sites have been expanded; (ii) more regulatory elements including transcription factor binding sites (TFBSs), enhancers, DNase I hypersensitive sites (DHSs), topologically associated domains (TADs), footprintss, methylations and open chromatin regions have been added; (iii) the associations among circRNAs, regulatory elements and variants have been identified; (iv) more experimentally supported variant-lncRNA/circRNA-disease/phenotype associations have been manually collected; (v) the sources of lncRNAs, circRNAs, SNPs, somatic mutations and RNA editing sites have been updated. Moreover, four flexible online tools including Genome Browser, Variant Mapper, Circos Plotter and Functional Annotation have been developed to retrieve, visualize and analyze the data. Collectively, LincSNP 3.0 provides associations among functional variants, regulatory elements, lncRNAs and circRNAs in diseases. It will serve as an important and continually updated resource for investigating functions and mechanisms of lncRNAs and circRNAs in diseases.

The CRISPR–Cas system provides adaptive immunity for bacteria and archaea to combat invading phages and plasmids. Phages evolved anti-CRISPR (Acr) proteins to neutralize the host CRISPR–Cas immune system as a counter-defense mechanism. AcrIF7 in Pseudomonas aeruginosa prophages strongly inhibits the type I-F CRISPR–Cas system. Here, we determined the solution structure of AcrIF7 and identified its target, Cas8f of the Csy complex. AcrIF7 adopts a novel β1β2α1α2β3 fold and interacts with the target DNA binding site of Cas8f. Notably, AcrIF7 competes with AcrIF2 for the same binding interface on Cas8f without common structural motifs. AcrIF7 binding to Cas8f is driven mainly by electrostatic interactions that require position-specific surface charges. Our findings suggest that Acrs of divergent origin may have acquired specificity to a common target through convergent evolution of their surface charge configurations.

Human Y-box binding protein 1 (YB-1) is a multifunctional protein and overexpressed in many types of cancer. It specifically recognizes DNA/RNA through a cold shock domain (CSD) and regulates nucleic acid metabolism. The C-terminal extension of CSD and the phosphorylation of S102 are indispensable for YB-1 function. Until now, the roles of the C-terminal extension and phosphorylation in gene transcription and translation are still largely unknown. Here, we solved the structure of human YB-1 CSD with a C-terminal extension sequence (CSDex). The structure reveals that the extension interacts with several residues in the conventional CSD and adopts a rigid structure instead of being disordered. Either deletion of this extension or phosphorylation of S102 destabilizes the protein and results in partial unfolding. Structural characterization of CSDex in complex with a ssDNA heptamer shows that all the seven nucleotides are involved in DNA–protein interactions and the C-terminal extension provides a unique DNA binding site. Our DNA-binding study indicates that CSDex can recognize more DNA sequences than previously thought and the phosphorylation reduces its binding to ssDNA dramatically. Our results suggest that gene transcription and translation can be regulated by changing the affinity of CSDex binding to DNA and RNA through phosphorylation, respectively.

G-quadruplexes are four-stranded nucleic acid structures involved in multiple cellular pathways including DNA replication and telomere maintenance. Such structures are formed by G-rich DNA sequences typified by telomeric DNA repeats. Whilst there is evidence for proteins that bind and regulate G-quadruplex formation, the molecular basis for this remains poorly understood. The budding yeast telomeric protein Rap1, originally identified as a transcriptional regulator functioning by recognizing double-stranded DNA binding sites, was one of the first proteins to be discovered to also bind and promote G-quadruplex formation in vitro. Here, we present the 2.4 Å resolution crystal structure of the Rap1 DNA-binding domain in complex with a G-quadruplex. Our structure not only provides a detailed insight into the structural basis for G-quadruplex recognition by a protein, but also gives a mechanistic understanding of how the same DNA-binding domain adapts to specifically recognize different DNA structures. The key observation is the DNA-recognition helix functions in a bimodal manner: In double-stranded DNA recognition one helix face makes electrostatic interactions with the major groove of DNA, whereas in G-quadruplex recognition a different helix face is used to make primarily hydrophobic interactions with the planar face of a G-tetrad. Ministry of Education (MOE) Published version Singapore Ministry of Education Academic Research Fund (AcRF) Tier 3 [MOE2012-T3-1-001]. Funding for open access charge: Singapore Ministry of Education Academic Research.

Modified DNA bases functionally distinguish the taxonomic forms of life—5-methylcytosine separates prokaryotes from eukaryotes and 5-hydroxymethylcytosine (5hmC) invertebrates from vertebrates. We demonstrate here that mouse endonuclease G (mEndoG) shows specificity for both 5hmC and Holliday junctions. The enzyme has higher affinity (>50-fold) for junctions over duplex DNAs. A 5hmC-modification shifts the position of the cut site and increases the rate of DNA cleavage in modified versus unmodified junctions. The crystal structure of mEndoG shows that a cysteine (Cys69) is positioned to recognize 5hmC through a thiol-hydroxyl hydrogen bond. Although this Cys is conserved from worms to mammals, a two amino acid deletion in the vertebrate relative to the invertebrate sequence unwinds an α-helix, placing the thiol of Cys69 into the mEndoG active site. Mutations of Cys69 with alanine or serine show 5hmC-specificity that mirrors the hydrogen bonding potential of the side chain (C–H < S–H < O–H). A second orthogonal DNA binding site identified in the mEndoG structure accommodates a second arm of a junction. Thus, the specificity of mEndoG for 5hmC and junctions derives from structural adaptations that distinguish the vertebrate from the invertebrate enzyme, thereby thereby supporting a role for 5hmC in recombination processes.

Studies of transcription regulation are often focused on binding of transcription factors (TFs) to a small number of promoters of interest. It is often assumed that TFs are in great excess to their binding sites (TFBSs) and competition for TFs between DNA sites is seldom considered. With increasing evidence that TFBSs are exceedingly abundant for many TFs and significant variations in TF and TFBS numbers occur during growth, the interplay between a TF and all TFBSs should not be ignored. Here, we use additional decoy DNA sites to quantitatively analyze how the relative abundance of a TF to its TFBSs impacts the steady-state level and onset time of gene expression for the auto-activated Escherichia coli PhoB response regulator. We show that increasing numbers of decoy sites progressively delayed transcription activation and lowered promoter activities. Perturbation of transcription regulation by additional TFBSs did not require extreme numbers of decoys, suggesting that PhoB is approximately at capacity for its DNA sites. Addition of decoys also converted a graded response to a bi-modal response. We developed a binding competition model that captures the major features of experimental observations, providing a quantitative framework to assess how variations in TFs and TFBSs influence transcriptional responses.

The most common form of DNA methylation involves the addition of a methyl group to a cytosine base in the context of a cytosine–phosphate–guanine (CpG) dinucleotide. Genomes from more primitive organisms are more abundant in CpG sites that, through the process of methylation, deamination and subsequent mutation to thymine–phosphate–guanine (TpG) sites, can produce new transcription factor binding sites. Here, we examined the evolutionary history of the over 36 000 glucocorticoid receptor (GR) consensus binding motifs in the human genome and identified a subset of them in regulatory regions that arose via a deamination and subsequent mutation event. GR can bind to both unmodified and methylated pre-GR binding sequences (GBSs) that contain a CpG site. Our structural analyses show that CpG methylation in a pre-GBS generates a favorable interaction with Arg447 mimicking that made with a TpG in a GBS. This methyl-specific recognition arose 420 million years ago and was conserved during the evolution of GR and likely helps fix the methylation on the relevant cytosines. Our study provides the first genetic, biochemical and structural evidence of high-affinity binding for the likely evolutionary precursor of extant TpG-containing GBS.

The identification of transcription factor binding sites and cis-regulatory motifs is a frontier whereupon the rules governing protein–DNA binding are being revealed. Here, we developed a new method (DEep Sequence and Shape mOtif or DESSO) for cis-regulatory motif prediction using deep neural networks and the binomial distribution model. DESSO outperformed existing tools, including DeepBind, in predicting motifs in 690 human ENCODE ChIP-sequencing datasets. Furthermore, the deep-learning framework of DESSO expanded motif discovery beyond the state-of-the-art by allowing the identification of known and new protein–protein–DNA tethering interactions in human transcription factors (TFs). Specifically, 61 putative tethering interactions were identified among the 100 TFs expressed in the K562 cell line. In this work, the power of DESSO was further expanded by integrating the detection of DNA shape features. We found that shape information has strong predictive power for TF–DNA binding and provides new putative shape motif information for human TFs. Thus, DESSO improves in the identification and structural analysis of TF binding sites, by integrating the complexities of DNA binding into a deep-learning framework.

A species-specific region, denoted SpG8-1b allowing hydroxycinnamic acids (HCAs) degradation is important for the transition between the two lifestyles (rhizospheric versus pathogenic) of the plant pathogen Agrobacterium fabrum. Indeed, HCAs can be either used as trophic resources and/or as induced-virulence molecules. The SpG8-1b region is regulated by two transcriptional regulators, namely, HcaR (Atu1422) and Atu1419. In contrast to HcaR, Atu1419 remains so far uncharacterized. The high-resolution crystal structures of two fortuitous citrate complexes, two DNA complexes and the apoform revealed that the tetrameric Atu1419 transcriptional regulator belongs to the VanR group of Pfam PF07729 subfamily of the large GntR superfamily. Until now, GntR regulators were described as dimers. Here, we showed that Atu1419 represses three genes of the HCAs catabolic pathway. We characterized both the effector and DNA binding sites and identified key nucleotides in the target palindrome. From promoter activity measurement using defective gene mutants, structural analysis and gel-shift assays, we propose N5,N10-methylenetetrahydrofolate as the effector molecule, which is not a direct product/substrate of the HCA degradation pathway. The Zn2+ ion present in the effector domain has both a structural and regulatory role. Overall, our work shed light on the allosteric mechanism of transcription employed by this GntR repressor.

Analysis of ENCODE long RNA-Seq and ChIP-seq (Chromatin Immunoprecipitation Sequencing) datasets for HepG2 and HeLa cell lines uncovered 1647 and 1958 transcripts that interfere with transcription factor binding to human enhancer domains. TFBSs (Transcription Factor Binding Sites) intersected by these ‘Enhancer Occlusion Transcripts’ (EOTrs) displayed significantly lower relative transcription factor (TF) binding affinities compared to TFBSs for the same TF devoid of EOTrs. Expression of most EOTrs was regulated in a cell line specific manner; analysis for the same TFBSs across cell lines, i.e. in the absence or presence of EOTrs, yielded consistently higher relative TF/DNA-binding affinities for TFBSs devoid of EOTrs. Lower activities of EOTr-associated enhancer domains coincided with reduced occupancy levels for histone tail modifications H3K27ac and H3K9ac. Similarly, the analysis of EOTrs with allele-specific expression identified lower activities for alleles associated with EOTrs. ChIA-PET (Chromatin Interaction Analysis by Paired-End Tag Sequencing) and 5C (Carbon Copy Chromosome Conformation Capture) uncovered that enhancer domains associated with EOTrs preferentially interacted with poised gene promoters. Analysis of EOTr regions with GRO-seq (Global run-on) data established the correlation of RNA polymerase pausing and occlusion of TF-binding. Our results implied that EOTr expression regulates human enhancer domains via transcriptional interference.

The rolling-circle replication is the most common mechanism for the replication of small plasmids carrying antibiotic resistance genes in Gram-positive bacteria. It is initiated by the binding and nicking of double-stranded origin of replication by a replication initiator protein (Rep). Duplex unwinding is then performed by the PcrA helicase, whose processivity is critically promoted by its interaction with Rep. How Rep and PcrA proteins interact to nick and unwind the duplex is not fully understood. Here, we have used magnetic tweezers to monitor PcrA helicase unwinding and its relationship with the nicking activity of Staphylococcus aureus plasmid pT181 initiator RepC. Our results indicate that PcrA is a highly processive helicase prone to stochastic pausing, resulting in average translocation rates of 30 bp s-1, while a typical velocity of 50 bp s-1 is found in the absence of pausing. Single-strand DNA binding protein did not affect PcrA translocation velocity but slightly increased its processivity. Analysis of the degree of DNA supercoiling required for RepC nicking, and the time between RepC nicking and DNA unwinding, suggests that RepC and PcrA form a protein complex on the DNA binding site before nicking. A comprehensive model that rationalizes these findings is presented., Ministerio de Economía y Competitividad (MINECO) [BFU2017-83794-P (AEI/FEDER, UE) to F.M.-H.]; European Research Council (ERC) under the European Union Horizon 2020 research and innovation grant agreemen t[681299 to F.M.-H.]; National Science Foundation [toS.H.L.]. C.C. was supported by a ‘Severo Ochoa’ postdoc-toral contract from the National Center of Biotechnology (CNB-CSIC); Funding for open access charge: EuropeanResearch Council [681299].

Bound by transcription factors, DNA motifs (i.e. transcription factor binding sites) are prevalent and important for gene regulation in different tissues at different developmental stages of eukaryotes. Although considerable efforts have been made on elucidating monomeric DNA motif patterns, our knowledge on heterodimeric DNA motifs are still far from complete. Therefore, we propose to develop a computational approach to synthesize a heterodimeric DNA motif from two monomeric DNA motifs. The approach is sequentially divided into two components (Phases A and B). In Phase A, we propose to develop the inference models on how two DNA monomeric motifs can be oriented and overlapped with each other at nucleotide level. In Phase B, given the two monomeric DNA motifs oriented, we further propose to develop DNA-binding family-specific input-output hidden Markov models (IOHMMs) to synthesize a heterodimeric DNA motif. To validate the approach, we execute and cross-validate it with the experimentally verified 618 heterodimeric DNA motifs across 49 DNA-binding family combinations. We observe that our approach can even “rescue" the existing heterodimeric DNA motif pattern (i.e. HOXB2_EOMES) previously published on Nature. Lastly, we apply the proposed approach to infer previously uncharacterized heterodimeric motifs. Their motif instances are supported by DNase accessibility, gene ontology, protein-protein interactions, in vivo ChIP-seq peaks, and even structural data from PDB. A public web-server is built for open accessibility and scientific impact. Its address is listed as follows: http://motif.cs.cityu.edu.hk/custom/MotifKirin.

The Cistrome Data Browser (DB) is a resource of human and mouse cis-regulatory information derived from ChIP-seq, DNase-seq and ATAC-seq chromatin profiling assays, which map the genome-wide locations of transcription factor binding sites, histone post-translational modifications and regions of chromatin accessible to endonuclease activity. Currently, the Cistrome DB contains approximately 47,000 human and mouse samples with about 24,000 newly collected datasets compared to the previous release two years ago. Furthermore, the Cistrome DB has a new Toolkit module with several features that allow users to better utilize the large-scale ChIP-seq, DNase-seq, and ATAC-seq data. First, users can query the factors which are likely to regulate a specific gene of interest. Second, the Cistrome DB Toolkit facilitates searches for factor binding, histone modifications, and chromatin accessibility in any given genomic interval shorter than 2Mb. Third, the Toolkit can determine the most similar ChIP-seq, DNase-seq, and ATAC-seq samples in terms of genomic interval overlaps with user-provided genomic interval sets. The Cistrome DB is a user-friendly, up-to-date, and well maintained resource, and the new tools will greatly benefit the biomedical research community. The database is freely available at http://cistrome.org/db, and the Toolkit is at http://dbtoolkit.cistrome.org.

The current version of the Gene Transcription Regulation Database (GTRD; http://gtrd.biouml.org) contains information about: (i) transcription factor binding sites (TFBSs) and transcription coactivators identified by ChIP-seq experiments for Homo sapiens, Mus musculus, Rattus norvegicus, Danio rerio, Caenorhabditis elegans, Drosophila melanogaster, Saccharomyces cerevisiae, Schizosaccharomyces pombe and Arabidopsis thaliana; (ii) regions of open chromatin and TFBSs (DNase footprints) identified by DNase-seq; (iii) unmappable regions where TFBSs cannot be identified due to repeats; (iv) potential TFBSs for both human and mouse using position weight matrices from the HOCOMOCO database. Raw ChIP-seq and DNase-seq data were obtained from ENCODE and SRA, and uniformly processed. ChIP-seq peaks were called using four different methods: MACS, SISSRs, GEM and PICS. Moreover, peaks for the same factor and peak calling method, albeit using different experiment conditions (cell line, treatment, etc.), were merged into clusters. To reduce noise, such clusters for different peak calling methods were merged into meta-clusters; these were considered to be non-redundant TFBS sets. Moreover, extended quality control was applied to all ChIP-seq data. Web interface to access GTRD was developed using the BioUML platform. It provides browsing and displaying information, advanced search possibilities and an integrated genome browser.

Super-enhancers are important for controlling and defining the expression of cell-specific genes. With research on human disease and biological processes, human H3K27ac ChIP-seq datasets are accumulating rapidly, creating the urgent need to collect and process these data comprehensively and efficiently. More importantly, many studies showed that super-enhancer-associated single nucleotide polymorphisms (SNPs) and transcription factors (TFs) strongly influence human disease and biological processes. Here, we developed a comprehensive human super-enhancer database (SEdb, http://www.licpathway.net/sedb) that aimed to provide a large number of available resources on human super-enhancers. The database was annotated with potential functions of super-enhancers in the gene regulation. The current version of SEdb documented a total of 331 601 super-enhancers from 542 samples. Especially, unlike existing super-enhancer databases, we manually curated and classified 410 available H3K27ac samples from >2000 ChIP-seq samples from NCBI GEO/SRA. Furthermore, SEdb provides detailed genetic and epigenetic annotation information on super-enhancers. Information includes common SNPs, motif changes, expression quantitative trait locus (eQTL), risk SNPs, transcription factor binding sites (TFBSs), CRISPR/Cas9 target sites and Dnase I hypersensitivity sites (DHSs) for in-depth analyses of super-enhancers. SEdb will help elucidate super-enhancer-related functions and find potential biological effects.

The REDfly database provides a comprehensive curation of experimentally-validated Drosophila transcriptional cis-regulatory elements and includes information on DNA sequence, experimental evidence, patterns of regulated gene expression, and more. Now in its thirteenth year, REDfly has grown to over 23 000 records of tested reporter gene constructs and 2200 tested transcription factor binding sites. Recent developments include the start of curation of predicted cis-regulatory modules in addition to experimentally-verified ones, improved search and filtering, and increased interaction with the authors of curated papers. An expanded data model that will capture information on temporal aspects of gene regulation, regulation in response to environmental and other non-developmental cues, sexually dimorphic gene regulation, and non-endogenous (ectopic) aspects of reporter gene expression is under development and expected to be in place within the coming year. REDfly is freely accessible at http://redfly.ccr.buffalo.edu, and news about database updates and new features can be followed on Twitter at @REDfly_database.

Genetic variation in cis-regulatory elements is thought to be a major driving force in morphological and physiological changes. However, identifying transcription factor binding events that code for complex traits remains a challenge, motivating novel means of detecting putatively important binding events. Using a curated set of 1154 high-quality transcription factor motifs, we demonstrate that independently eroded binding sites are enriched for independently lost traits in three distinct pairs of placental mammals. We show that these independently eroded events pinpoint the loss of hindlimbs in dolphin and manatee, degradation of vision in naked mole-rat and star-nosed mole, and the loss of external testes in white rhinoceros and Weddell seal. We additionally show that our method may also be utilized with more than two species. Our study exhibits a novel methodology to detect cis-regulatory mutations which help explain a portion of the molecular mechanism underlying complex trait formation and loss.

The accumulation of large epigenomics data consortiums provides us with the opportunity to extrapolate existing knowledge to new cell types and conditions. We propose Epitome, a deep neural network that learns similarities of chromatin accessibility between well characterized reference cell types and a query cellular context, and copies over signal of transcription factor binding and modification of histones from reference cell types when chromatin profiles are similar to the query. Epitome achieves state-of-the-art accuracy when predicting transcription factor binding sites on novel cellular contexts and can further improve predictions as more epigenetic signals are collected from both reference cell types and the query cellular context of interest.

The dynamics and mechanism of how site-specific DNA-bending proteins initially interrogate potential binding sites prior to recognition have remained elusive for most systems. Here we present these dynamics for Integration Host factor (IHF), a nucleoid-associated architectural protein, using a μs-resolved T-jump approach. Our studies show two distinct DNA-bending steps during site recognition by IHF. While the faster (∼100 μs) step is unaffected by changes in DNA or protein sequence that alter affinity by >100-fold, the slower (1–10 ms) step is accelerated ∼5-fold when mismatches are introduced at DNA sites that are sharply kinked in the specific complex. The amplitudes of the fast phase increase when the specific complex is destabilized and decrease with increasing [salt], which increases specificity. Taken together, these results indicate that the fast phase is non-specific DNA bending while the slow phase, which responds only to changes in DNA flexibility at the kink sites, is specific DNA kinking during site recognition. Notably, the timescales for the fast phase overlap with one-dimensional diffusion times measured for several proteins on DNA, suggesting that these dynamics reflect partial DNA bending during interrogation of potential binding sites by IHF as it scans DNA.

International audience; The Saccharomyces cerevisiae Pif1 protein (ScPif1p) is the prototypical member of the Pif1 family of DNA helicases. ScPif1p is involved in the maintenance of mitochondrial, ribosomal and telomeric DNA and suppresses genome instability at G-quadruplex motifs. Here, we report the crystal structures of a truncated ScPif1p (ScPif1p237−780) in complex with different ssDNAs. Our results have revealed that a yeast-specific insertion domain protruding from the 2B domain folds as a bundle bearing an α-helix, α16. The α16 helix regulates the helicase activities of ScPif1p through interactions with the previously identified loop3. Furthermore, a biologically relevant dimeric structure has been identified, which can be further specifically stabilized by G-quadruplex DNA. Basing on structural analyses and mutational studies with DNA binding and unwinding assays, a potential G-quadruplex DNA binding site in ScPif1p monomers is suggested. Our results also show that ScPif1p uses the Q-motif to preferentially hydrolyze ATP, and a G-rich tract is preferentially recognized by more residues, consistent with previous biochemical observations. These findings provide a structural and mechanistic basis for understanding the multifunctional ScPif1p.

Bacteria adapt to changes in their environment via differential gene expression mediated by DNA binding transcriptional regulators. The PRODORIC2 database hosts one of the largest collections of DNA binding sites for prokaryotic transcription factors. It is the result of the thoroughly redesigned PRODORIC database. PRODORIC2 is more intuitive and user-friendly. Besides significant technical improvements, the new update offers more than 1000 new transcription factor binding sites and 110 new position weight matrices for genome-wide pattern searches with the Virtual Footprint tool. Moreover, binding sites deduced from high-throughput experiments were included. Data for 6 new bacterial species including bacteria of the Rhodobacteraceae family were added. Finally, a comprehensive collection of sigma- and transcription factor data for the nosocomial pathogen Clostridium difficile is now part of the database. PRODORIC2 is publicly available at http://www.prodoric2.de.

Fanconi anemia (FA) is a recessive genetic disorder characterized by congenital abnormalities, progressive bone-marrow failure, and cancer susceptibility. The FA pathway consists of at least 21 FANC genes (FANCA-FANCV), and the encoded protein products interact in a common cellular pathway to gain resistance against DNA interstrand crosslinks. After DNA damage, FANCD2 is monoubiquitinated and accumulates on chromatin. FANCD2 plays a central role in the FA pathway, using yet unidentified DNA binding regions. By using synthetic peptide mapping and DNA binding screen by electromobility shift assays, we found that FANCD2 bears two major DNA binding domains predominantly consisting of evolutionary conserved lysine residues. Furthermore, one domain at the N-terminus of FANCD2 bears also nuclear localization sequences for the protein. Mutations in the bifunctional DNA binding/NLS domain lead to a reduction in FANCD2 monoubiquitination and increase in mitomycin C sensitivity. Such phenotypes are not fully rescued by fusion with an heterologous NLS, which enable separation of DNA binding and nuclear import functions within this domain that are necessary for FANCD2 functions. Collectively, our results enlighten the importance of DNA binding and NLS residues in FANCD2 to activate an efficient FA pathway.

Target search as performed by DNA-binding proteins is a complex process, in which multiple factors contribute to both thermodynamic discrimination of the target sequence from overwhelmingly abundant off-target sites and kinetic acceleration of dynamic sequence interrogation. TRF1, the protein that binds to telomeric tandem repeats, faces an intriguing variant of the search problem where target sites are clustered within short fragments of chromosomal DNA. In this study, we use extensive (>0.5 ms in total) MD simulations to study the dynamical aspects of sequence-specific binding of TRF1 at both telomeric and non-cognate DNA. For the first time, we describe the spontaneous formation of a sequence-specific native protein–DNA complex in atomistic detail, and study the mechanism by which proteins avoid off-target binding while retaining high affinity for target sites. Our calculated free energy landscapes reproduce the thermodynamics of sequence-specific binding, while statistical approaches allow for a comprehensive description of intermediate stages of complex formation.

AGRO100 (also known as AS1411) is a G-rich oligonucleotide that has long been established as a potent anti-cancer aptamer. However, the structure of AGRO100 remained unresolved, due to the co-existence of multiple different G-quadruplex conformations. We identified a DNA sequence named AT11, derived from AGRO100, which formed a single major G-quadruplex conformation and exhibited a similar anti-proliferative activity as AGRO100. The solution structure of AT11 revealed a four-layer G-quadruplex comprising of two propeller-type parallel-stranded subunits connected through a central linker. The stacking between the two subunits occurs at the 3΄-end of the first block and the 5΄-end of the second block. The structure of the anti-proliferative DNA sequence AT11 will allow greater understanding on the G-quadruplex folding principles and aid in structural optimization of anti-proliferative oligonucleotides. NRF (Natl Research Foundation, S’pore) MOE (Min. of Education, S’pore) Published version

Transcription factors are important gene regulators with distinctive roles in development, cell signaling and cell cycling, and they have been associated with many diseases. The ConTra v3 web server allows easy visualization and exploration of predicted transcription factor binding sites (TFBSs) in any genomic region surrounding coding or non-coding genes. In this updated version, with a completely re-implemented user interface using latest web technologies, users can choose from nine reference organisms ranging from human to yeast. ConTra v3 can analyze promoter regions, 5΄-UTRs, 3΄-UTRs and introns or any other genomic region of interest. Thousands of position weight matrices are available to choose from for detecting specific binding sites. Besides this visualization option, additional new exploration functionality is added to the tool that will automatically detect TFBSs having at the same time the highest regulatory potential, the highest conservation scores of the genomic regions covered by the predicted TFBSs and strongest co-localizations with genomic regions exhibiting regulatory activity. The ConTra v3 web server is freely available at http://bioit2.irc.ugent.be/contra/v3.

FadR is a fatty acyl-CoA dependent transcription factor that regulates genes encoding proteins involved in fatty-acid degradation and synthesis pathways. In this study, the crystal structures of Bacillus halodurans FadR, which belong to the TetR family, have been determined in three different forms: ligand-bound, ligand-free and DNA-bound at resolutions of 1.75, 2.05 and 2.80 Å, respectively. Structural and functional data showed that B. halodurans FadR was bound to its operator site without fatty acyl-CoAs. Structural comparisons among the three different forms of B. halodurans FadR revealed that the movement of DNA binding domains toward the operator DNA was blocked upon binding of ligand molecules. These findings suggest that the TetR family FadR negatively regulates the genes involved in fatty acid metabolism by binding cooperatively to the operator DNA as a dimer of dimers.

Isovaleryl coenzyme A (IV-CoA) is an important building block of iso-fatty acids. In myxobacteria, IV-CoA is essential for the formation of signaling molecules involved in fruiting body formation. Leucine degradation is the common source of IV-CoA, but a second, de novo biosynthetic route to IV-CoA termed AIB (alternative IV-CoA biosynthesis) was recently discovered in M. xanthus. The AIB-operon contains the TetR-like transcriptional regulator AibR, which we characterize in this study. We demonstrate that IV-CoA binds AibR with micromolar affinity and show by gelshift experiments that AibR interacts with the promoter region of the AIB-operon once IV-CoA is present. We identify an 18-bp near-perfect palindromic repeat as containing the AibR operator and provide evidence that AibR also controls an additional genomic locus coding for a putative acetyl–CoA acetyltransferase. To elucidate atomic details, we determined crystal structures of AibR in the apo, the IV-CoA- and the IV-CoA-DNA-bound state to 1.7 Å, 2.35 Å and 2.92 Å, respectively. IV-CoA induces partial unfolding of an α-helix, which allows sequence-specific interactions between AibR and its operator. This study provides insights into AibR-mediated regulation and shows that AibR functions in an unusual TetR-like manner by blocking transcription not in the ligand-free but in the effector-bound state.

Super-enhancers (SEs) are critical for the transcriptional regulation of gene expression. We developed the super-enhancer archive version 3.0 (SEA v. 3.0, http://sea.edbc.org) to extend SE research. SEA v. 3.0 provides the most comprehensive archive to date, consisting of 164 545 super-enhancers. Of these, 80 549 are newly identified from 266 cell types/tissues/diseases using an optimized computational strategy, and 52 have been experimentally confirmed with manually curated references. We now support super-enhancers in 11 species including 7 new species (zebrafish, chicken, chimp, rhesus, sheep, Xenopus tropicalis and stickleback). To facilitate super-enhancer functional analysis, we added several new regulatory datasets including 3 361 785 typical enhancers, chromatin interactions, SNPs, transcription factor binding sites and SpCas9 target sites. We also updated or developed new criteria query, genome visualization and analysis tools for the archive. This includes a tool based on Shannon Entropy to evaluate SE cell type specificity, a new genome browser that enables the visualization of SE spatial interactions based on Hi-C data, and an enhanced enrichment analysis interface that provides online enrichment analyses of SE related genes. SEA v. 3.0 provides a comprehensive database of all available SE information across multiple species, and will facilitate super-enhancer research, especially as related to development and disease.

TFBSshape (https://tfbsshape.usc.edu) is a motif database for analyzing structural profiles of transcription factor binding sites (TFBSs). The main rationale for this database is to be able to derive mechanistic insights in protein–DNA readout modes from sequencing data without available structures. We extended the quantity and dimensionality of TFBSshape, from mostly in vitro to in vivo binding and from unmethylated to methylated DNA. This new release of TFBSshape improves its functionality and launches a responsive and user-friendly web interface for easy access to the data. The current expansion includes new entries from the most recent collections of transcription factors (TFs) from the JASPAR and UniPROBE databases, methylated TFBSs derived from in vitro high-throughput EpiSELEX-seq binding assays and in vivo methylated TFBSs from the MeDReaders database. TFBSshape content has increased to 2428 structural profiles for 1900 TFs from 39 different species. The structural profiles for each TFBS entry now include 13 shape features and minor groove electrostatic potential for standard DNA and four shape features for methylated DNA. We improved the flexibility and accuracy for the shape-based alignment of TFBSs and designed new tools to compare methylated and unmethylated structural profiles of TFs and methods to derive DNA shape-preserving nucleotide mutations in TFBSs.

Escherichia coli single strand (ss) DNA binding (SSB) protein protects ssDNA intermediates and recruits at least 17 SSB interacting proteins (SIPs) during genome maintenance. The SSB C-termini contain a 9 residue acidic tip and a 56 residue intrinsically disordered linker (IDL). The acidic tip interacts with SIPs; however a recent proposal suggests that the IDL may also interact with SIPs. Here we examine the binding to four SIPs (RecO, PriC, PriA and χ subunit of DNA polymerase III) of three peptides containing the acidic tip and varying amounts of the IDL. Independent of IDL length, we find no differences in peptide binding to each individual SIP indicating that binding is due solely to the acidic tip. However, the tip shows specificity, with affinity decreasing in the order: RecO > PriA ∼ χ > PriC. Yet, RecO binding to the SSB tetramer and an SSB–ssDNA complex show significant thermodynamic differences compared to the peptides alone, suggesting that RecO interacts with another region of SSB, although not the IDL. SSB containing varying IDL deletions show different binding behavior, with the larger linker deletions inhibiting RecO binding, likely due to increased competition between the acidic tip interacting with DNA binding sites within SSB.

Super-enhancers are established through the interactions of several enhancers and a large number of proteins, including transcription factors and co-regulators; however, the formation of these interactions is poorly understood. By re-analysing previously published estrogen receptor alpha (ERα) ChIP-seq data sets derived from the MCF-7 cell line, we observed that in the absence of stimulation, future super-enhancers are represented by one or a few transcription factor binding event(s) and these extraordinary enhancers possess a response element largely specific to the ERα dimer. Upon hormonal stimulation, these primary binding sites are surrounded by a large amount of ERα and the critical components of active enhancers, such as P300 and MED1, and together with neighbouring sites bound by newly recruited ERα, they generate the functional super-enhancers. To further validate the role of canonical elements in super-enhancer formation, we investigated some additional signal-dependent transcription factors, confirming that certain, distinguished binding elements have a general organizer function. These results suggest that certain signal-specific transcription factors guide super-enhancer formation upon binding to strong response elements. These findings may reshape the current understanding of how these regulatory units assemble, highlighting the involvement of DNA elements instead of protein–protein interactions.

SNP2TFBS is a computational resource intended to support researchers investigating the molecular mechanisms underlying regulatory variation in the human genome. The database essentially consists of a collection of text files providing specific annotations for human single nucleotide polymorphisms (SNPs), namely whether they are predicted to abolish, create or change the affinity of one or several transcription factor (TF) binding sites. A SNP's effect on TF binding is estimated based on a position weight matrix (PWM) model for the binding specificity of the corresponding factor. These data files are regenerated at regular intervals by an automatic procedure that takes as input a reference genome, a comprehensive SNP catalogue and a collection of PWMs. SNP2TFBS is also accessible over a web interface, enabling users to view the information provided for an individual SNP, to extract SNPs based on various search criteria, to annotate uploaded sets of SNPs or to display statistics about the frequencies of binding sites affected by selected SNPs. Homepage: http://ccg.vital-it.ch/snp2tfbs/.

We describe LincSNP 2.0 (http://bioinfo.hrbmu.edu.cn/LincSNP), an updated database that is used specifically to store and annotate disease-associated single nucleotide polymorphisms (SNPs) in human long non-coding RNAs (lncRNAs) and their transcription factor binding sites (TFBSs). In LincSNP 2.0, we have updated the database with more data and several new features, including (i) expanding disease-associated SNPs in human lncRNAs; (ii) identifying disease-associated SNPs in lncRNA TFBSs; (iii) updating LD-SNPs from the 1000 Genomes Project; and (iv) collecting more experimentally supported SNP-lncRNA-disease associations. Furthermore, we developed three flexible online tools to retrieve and analyze the data. Linc-Mart is a convenient way for users to customize their own data. Linc-Browse is a tool for all data visualization. Linc-Score predicts the associations between lncRNA and disease. In addition, we provided users a newly designed, user-friendly interface to search and download all the data in LincSNP 2.0 and we also provided an interface to submit novel data into the database. LincSNP 2.0 is a continually updated database and will serve as an important resource for investigating the functions and mechanisms of lncRNAs in human diseases.

To understand the epigenetic regulation of transcriptional response of macrophages during early-stage M. tuberculosis (Mtb) infection, we performed ChIPseq analysis of H3K4 monomethylation (H3K4me1), a marker of poised or active enhancers. De novo H3K4me1 peaks in infected cells were associated with genes implicated in host defenses and apoptosis. Our analysis revealed that 40% of de novo regions contained human/primate-specific Alu transposable elements, enriched in the AluJ and S subtypes. These contained several transcription factor binding sites, including those for members of the MEF2 and ATF families, and LXR and RAR nuclear receptors, all of which have been implicated in macrophage differentiation, survival, and responses to stress and infection. Combining bioinformatics, molecular genetics, and biochemical approaches, we linked genes adjacent to H3K4me1-associated Alu repeats to macrophage metabolic responses against Mtb infection. In particular, we show that LXRα signaling, which reduced Mtb viability 18-fold by altering cholesterol metabolism and enhancing macrophage apoptosis, can be initiated at response elements present in Alu repeats. These studies decipher the mechanism of early macrophage transcriptional responses to Mtb, highlighting the role of Alu element transposition in shaping human transcription programs during innate immunity.

Tyrosyl DNA phosphodiesterase 2 (TDP2) is a multifunctional protein implicated in DNA repair, signal transduction and transcriptional regulation. In its DNA repair role, TDP2 safeguards genome integrity by hydrolyzing 5'-tyrosyl DNA adducts formed by abortive topoisomerase II (Top2) cleavage complexes to allow error-free repair of DNA double-strand breaks, thereby conferring cellular resistance against Top2 poisons. TDP2 consists of a C-terminal catalytic domain responsible for its phosphodiesterase activity, and a functionally uncharacterized N-terminal region. Here, we demonstrate that this N-terminal region contains a ubiquitin (Ub)-associated (UBA) domain capable of binding multiple forms of Ub with distinct modes of interactions and preference for either K48- or K63-linked polyUbs over monoUb. The structure of TDP2 UBA bound to monoUb shows a canonical mode of UBA-Ub interaction. However, the absence of the highly conserved MGF motif and the presence of a fourth α-helix make TDP2 UBA distinct from other known UBAs. Mutations in the TDP2 UBA-Ub binding interface do not affect nuclear import of TDP2, but severely compromise its ability to repair Top2-mediated DNA damage, thus establishing the importance of the TDP2 UBA-Ub interaction in DNA repair. The differential binding to multiple Ub forms could be important for responding to DNA damage signals under different contexts or to support the multi-functionality of TDP2.

Fragile X mental retardation protein (FMRP) is a multifunctional RNA-binding protein with crucial roles in neuronal development and function. Efforts aimed at elucidating how FMRP target mRNAs are selected have produced divergent sets of target mRNA and putative FMRP-bound motifs, and a clear understanding of FMRP's binding determinants has been lacking. To clarify FMRP's binding to its target mRNAs, we produced a shared dataset of FMRP consensus binding sequences (FCBS), which were reproducibly identified in two published FMRP CLIP sequencing datasets. This comparative dataset revealed that of the various sequence and structural motifs that have been proposed to specify FMRP binding, the short sequence motifs TGGA and GAC were corroborated, and a novel TAY motif was identified. In addition, the distribution of the FCBS set demonstrates that FMRP preferentially binds to the coding region of its targets but also revealed binding along 3′ UTRs in a subset of target mRNAs. Beyond probing these putative motifs, the FCBS dataset of reproducibly identified FMRP binding sites is a valuable tool for investigating FMRP targets and function.

Single-stranded DNA binding proteins (SSBs) play an important role in DNA processing events such as replication, recombination and repair. Human single-stranded DNA binding protein 1 (hSSB1/NABP2/OBFC2B) contains a single oligosaccharide/oligonucleotide binding (OB) domain followed by a charged C-terminus and is structurally homologous to the SSB from the hyperthermophilic crenarchaeote Sulfolobus solfataricus. Recent work has revealed that hSSB1 is critical to homologous recombination and numerous other important biological processes such as the regulation of telomeres, the maintenance of DNA replication forks and oxidative damage repair. Since the ability of hSSB1 to directly interact with single-stranded DNA (ssDNA) is paramount for all of these processes, understanding the molecular details of ssDNA recognition is essential. In this study, we have used solution-state nuclear magnetic resonance in combination with biophysical and functional experiments to structurally analyse ssDNA binding by hSSB1. We reveal that ssDNA recognition in solution is modulated by base-stacking of four key aromatic residues within the OB domain. This DNA binding mode differs significantly from the recently determined crystal structure of the SOSS1 complex containing hSSB1 and ssDNA. Our findings elucidate the detailed molecular mechanism in solution of ssDNA binding by hSSB1, a major player in the maintenance of genomic stability.

Transactivation by the ETS family of transcription factors, whose members share structurally conserved DNA-binding domains, is variably sensitive to methylation of their target genes. The mechanism by which DNA methylation controls ETS proteins remains poorly understood. Uncertainly also pervades the effects of hemi-methylated DNA, which occurs following DNA replication and in response to hypomethylating agents, on site recognition by ETS proteins. To address these questions, we measured the affinities of two sequence-divergent ETS homologs, PU.1 and Ets-1, to DNA sites harboring a hemi- and fully methylated CpG dinucleotide. While the two proteins bound unmethylated DNA with indistinguishable affinity, their affinities to methylated DNA are markedly heterogeneous and exhibit major energetic coupling between the two CpG methylcytosines. Analysis of simulated DNA and existing co-crystal structures revealed that hemi-methylation induced non-local backbone and groove geometries that were not conserved in the fully methylated state. Indirect readout of these perturbations was differentially achieved by the two ETS homologs, with the distinctive interfacial hydration in PU.1/DNA binding moderating the inhibitory effects of DNA methylation on binding. This data established a biophysical basis for the pioneering properties associated with PU.1, which robustly bound fully methylated DNA, but not Ets-1, which was substantially inhibited.

Gene expression is a multi-step process involving many layers of regulation. The main regulators of the pathway are DNA and RNA binding proteins. While over the years, a large number of DNA and RNA binding proteins have been identified and extensively studied, it is still expected that many other proteins, some with yet another known function, are awaiting to be discovered. Here we present a new web server, BindUP, freely accessible through the website http://bindup.technion.ac.il/, for predicting DNA and RNA binding proteins using a non-homology-based approach. Our method is based on the electrostatic features of the protein surface and other general properties of the protein. BindUP predicts nucleic acid binding function given the proteins three-dimensional structure or a structural model. Additionally, BindUP provides information on the largest electrostatic surface patches, visualized on the server. The server was tested on several datasets of DNA and RNA binding proteins, including proteins which do not possess DNA or RNA binding domains and have no similarity to known nucleic acid binding proteins, achieving very high accuracy. BindUP is applicable in either single or batch modes and can be applied for testing hundreds of proteins simultaneously in a highly efficient manner.

Uracil-DNA glycosylases (UDGs) are highly conserved proteins that can be found in a wide range of organisms, and are involved in the DNA repair and host defense systems. UDG activity is controlled by various cellular factors, including the uracil-DNA glycosylase inhibitors, which are DNA mimic proteins that prevent the DNA binding sites of UDGs from interacting with their DNA substrate. To date, only three uracil-DNA glycosylase inhibitors, phage UGI, p56, and Staphylococcus aureus SAUGI, have been determined. We show here that SAUGI has differential inhibitory effects on UDGs from human, bacteria, Herpes simplex virus (HSV; human herpesvirus 1) and Epstein-Barr virus (EBV; human herpesvirus 4). Newly determined crystal structures of SAUGI/human UDG and a SAUGI/HSVUDG complex were used to explain the differential binding activities of SAUGI on these two UDGs. Structural-based protein engineering was further used to modulate the inhibitory ability of SAUGI on human UDG and HSVUDG. The results of this work extend our understanding of DNA mimics as well as potentially opening the way for novel therapeutic applications for this kind of protein.

Non-LTR retrotransposons are an important class of mobile elements that insert into host DNA by target-primed reverse transcription (TPRT). Non-LTR retrotransposons must bind to their mRNA, recognize and cleave their target DNA, and perform TPRT at the site of DNA cleavage. As DNA binding and cleavage are such central parts of the integration reaction, a better understanding of the endonuclease encoded by non-LTR retrotransposons is needed. This paper explores the R2 endonuclease domain from Bombyx mori using in vitro studies and in silico modeling. Mutations in conserved sequences located across the putative PD-(D/E)XK endonuclease domain reduced DNA cleavage, DNA binding and TPRT. A mutation at the beginning of the first α-helix of the modeled endonuclease obliterated DNA cleavage and greatly reduced DNA binding. It also reduced TPRT when tested on pre-cleaved DNA substrates. The catalytic K was located to a non-canonical position within the second α-helix. A mutation located after the fourth β-strand reduced DNA binding and cleavage. The motifs that showed impaired activity form an extensive basic region. The R2 biochemical and structural data are compared and contrasted with that of two other well characterized PD-(D/E)XK endonucleases, restriction endonucleases and archaeal Holliday junction resolvases.

We determined the effects of DNA damage caused by ultraviolet radiation on gene expression in Leptospira interrogans using DNA microarrays. These data were integrated with DNA binding in vivo of LexA1, a regulator of the DNA damage response, assessed by chromatin immunoprecipitation and massively parallel DNA sequencing (ChIP-seq). In response to DNA damage, Leptospira induced expression of genes involved in DNA metabolism, in mobile genetic elements and defective prophages. The DNA repair genes involved in removal of photo-damage (e.g. nucleotide excision repair uvrABC, recombinases recBCD and resolvases ruvABC) were not induced. Genes involved in various metabolic pathways were down regulated, including genes involved in cell growth, RNA metabolism and the tricarboxylic acid cycle. From ChIP-seq data, we observed 24 LexA1 binding sites located throughout chromosome 1 and one binding site in chromosome 2. Expression of many, but not all, genes near those sites was increased following DNA damage. Binding sites were found as far as 550 bp upstream from the start codon, or 1 kb into the coding sequence. Our findings indicate that there is a shift in gene expression following DNA damage that represses genes involved in cell growth and virulence, and induces genes involved in mutagenesis and recombination., Conselho Nacional de Pesquisas (Brazil), Fundacao de Amparo a Pesquisa do Estado de Sao Paulo, Fundacao Butantan, National Institutes of Health (U.S.) (National Institute of General Medical Sciences (U.S.) Award GM41934)