The role of response elements organization in transcription factor selectivity: the IFN-β enhanceosome example

Transcriptional kinetic synergy: A complex landscape revealed by integrating modeling and synthetic biology

Transcription factors (TFs) control gene expression, often acting synergistically. Classical thermodynamic models offer a biophysical explanation for synergy based on binding cooperativity and regulated recruitment of RNA polymerase. Because transcription requires polymerase to transition through multiple states, recent work suggests that "kinetic synergy" can arise through TFs acting on distinct steps of the transcription cycle. These types of synergy are not mutually exclusive and are difficult to disentangle conceptually and experimentally. Here, we model and build a synthetic circuit in which TFs bind to a single shared site on DNA, such that TFs cannot synergize by simultaneous binding. We model mRNA production as a function of both TF binding and regulation of the transcription cycle, revealing a complex landscape dependent on TF concentration, DNA binding affinity, and regulatory activity. We use synthetic TFs to confirm that the transcription cycle must be integrated with recruitment for a quantitative understanding of gene regulation.

Proteome-centric cross-omics characterization and integrated network analyses of triple-negative breast cancer

We report a comprehensive proteomic study of a 90-case cohort of paired samples of triple-negative breast cancer (TNBC) in quantification, phosphorylation, and DNA-binding capacity. Four integrative subtypes (iP-1-4) are stratified on the basis of global proteome and phosphoproteome, each of which exhibits distinct molecular and pathway features. Scaffold and co-expression network analyses of three proteomic datasets, integrated with those from genome and transcriptome of the same cohort, reveal key pathways and master regulators that, characteristic of TNBC subtypes, play important regulatory roles within and between scaffold sub-structures and co-expression communities. We find that NAE1 is a potential drug target for subtype iP-1, and a series of key molecules in fatty acid metabolism, such as AKT1/FASN, are plausible targets for subtype iP-2. Libraries of proteins, pathways and networks of TNBC provide a valuable molecular infrastructure for further clinical exploration and in-depth studies of the molecular mechanisms of the disease.

Of Keeping and Tipping the Balance: Host Regulation and Viral Modulation of IRF3-Dependent IFNB1 Expression

The type I interferon (IFN) response is a principal component of our immune system that allows to counter a viral attack immediately upon viral entry into host cells. Upon engagement of aberrantly localised nucleic acids, germline-encoded pattern recognition receptors convey their find via a signalling cascade to prompt kinase-mediated activation of a specific set of five transcription factors. Within the nucleus, the coordinated interaction of these dimeric transcription factors with coactivators and the basal RNA transcription machinery is required to access the gene encoding the type I IFN IFNβ (IFNB1). Virus-induced release of IFNβ then induces the antiviral state of the system and mediates further mechanisms for defence. Due to its key role during the induction of the initial IFN response, the activity of the transcription factor interferon regulatory factor 3 (IRF3) is tightly regulated by the host and fiercely targeted by viral proteins at all conceivable levels. In this review, we will revisit the steps enabling the trans-activating potential of IRF3 after its activation and the subsequent assembly of the multi-protein complex at the IFNβ enhancer that controls gene expression. Further, we will inspect the regulatory mechanisms of these steps imposed by the host cell and present the manifold strategies viruses have evolved to intervene with IFNβ transcription downstream of IRF3 activation in order to secure establishment of a productive infection.

Complex sex-biased antibody responses: estrogen receptors bind estrogen response elements centered within immunoglobulin heavy chain gene enhancers

Nuclear hormone receptors including the estrogen receptor (ERα) and the retinoic acid receptor regulate a plethora of biological functions including reproduction, circulation and immunity. To understand how estrogen and other nuclear hormones influence antibody production, we characterized total serum antibody isotypes in female and male mice of C57BL/6J, BALB/cJ and C3H/HeJ mouse strains. Antibody levels were higher in females compared to males in all strains and there was a female preference for IgG2b production. Sex-biased patterns were influenced by vitamin levels, and by antigen specificity toward influenza virus or pneumococcus antigens. To help explain sex biases, we examined the direct effects of estrogen on immunoglobulin heavy chain sterile transcript production among purified, lipopolysaccharide-stimulated B cells. Supplemental estrogen in B-cell cultures significantly increased immunoglobulin heavy chain sterile transcripts. Chromatin immunoprecipitation analyses of activated B cells identified significant ERα binding to estrogen response elements (EREs) centered within enhancer elements of the immunoglobulin heavy chain locus, including the Eµ enhancer and hypersensitive site 1,2 (HS1,2) in the 3' regulatory region. The ERE in HS1,2 was conserved across animal species, and in humans marked a site of polymorphism associated with the estrogen-augmented autoimmune disease, lupus. Taken together, the results highlight: (i) the important targets of ERα in regulatory regions of the immunoglobulin heavy chain locus that influence antibody production, and (ii) the complexity of mechanisms by which estrogen instructs sex-biased antibody production profiles.

Molecular mechanisms of the protein-protein interaction-regulated binding specificity of basic-region leucine zipper transcription factors

It is well known that the DNA-binding specificity of transcription factors (TFs) is influenced by protein-protein interactions (PPIs). However, the underlying molecular mechanisms remain largely unknown. In this work, we adopted the cAMP-response element-binding protein (CREB) of the basic leucine zipper (bZIP) TF family as a model system, and a workflow of combined bioinformatics and molecular modeling analysis of protein-DNA interaction was tested. First, the multiple sequence alignment and SDPsite method were used to find potential bZIP family binding specificity determining positions (SDPs) within the protein-protein interaction region. Second, the mutation system was analyzed using molecular dynamics simulation. Molecular mechanics Poisson-Boltzmann surface area (MM/PBSA) free energy calculations confirmed the enhancement of the binding affinity of the mutation, which was in agreement with experimental results. The root mean square fluctuation (RMSF) and hydrogen bonding changes suggested an open and close protein dimerization process after the system was mutated, which resulted in the change of the hydrogen bonding of the protein-DNA interface and a slight conformational change. We believe that this work will contribute to understanding the protein-protein interaction-regulated binding specificity of bZIP transcription factors.

The Functionality and Evolution of Eukaryotic Transcriptional Enhancers

Enhancers regulate precise spatial and temporal patterns of gene expression in eukaryotes and, moreover, evolutionary changes in these modular cis-regulatory elements may represent the predominant genetic basis for phenotypic evolution. Here, we review approaches to identify and functionally analyze enhancers and their transcription factor binding sites, including assay for transposable-accessible chromatin-sequencing (ATAC-Seq) and clustered regularly interspaced short palindromic repeats (CRISPR)/Cas9, respectively. We also explore enhancer functionality, including how transcription factor binding sites combine to regulate transcription, as well as research on shadow and super enhancers, and how enhancers can act over great distances and even in trans. Finally, we discuss recent theoretical and empirical data on how transcription factor binding sites and enhancers evolve. This includes how the function of enhancers is maintained despite the turnover of transcription factor binding sites as well as reviewing studies where mutations in enhancers have been shown to underlie morphological change.

Enhanceosome transcription factors preferentially dimerize with high mobility group proteins

Background

The enhanceosome is an enhancer located upstream of the human interferon β gene, bound by transcription factor (TF) complex of extremely rigid structure. Within these rigid constraints, even a slight change of distances between transcription factor binding sites (TFBS) results in loss of functionality of the enhanceosome. We hypothesized that smaller subunits of the enhanceosome may entail TF complex formation in other regulatory regions.

Results

In order to verify this hypothesis we systematically searched for dimerization preferences of the TFs that have TFBS in the enhanceosome. For this we utilized our recently developed tool, TACO. We performed this computational experiment in a cell-type–specific manner by utilizing cell-type–specific DNase-seq data for 105 human cell types. We also used 20 TRANSFAC motifs comprising not only the usual TFs constituting the enhanceosome but also the architectural proteins of High Mobility Group I(Y) (HMG I). A similar experiment used 42 DNase-seq data sets for mouse cell types. We found 137 statistically significant dimer predictions in the human genome, and 37 predictions in the mouse genome, that matched the positioning on the enhanceosome with ±2 bp tolerance. To characterize these predicted TF dimers, we performed functional analysis (Gene Ontology enrichment) for sets of genes which were in the neighbourhood of predicted dimer instances. A notable feature of these instances is that (1) most of them are located in introns of genes, (2) they are enriched in regulatory states, and (3) those instances that are located near transcription start sites are enriched for inclusion in computationally predicted enhancers. We also investigated similarity of dimer predictions between human and mouse.

Conclusions

It follows from our experiments that, except for homodimer formed by IRF proteins, the rest of the dimers were formed exclusively between one of the transcriptional activators (ATF-2/c-Jun and IRF) and a HMG I protein. NF- κB did not participate in forming dimers with other proteins. Dimers predicted in mouse were fully contained in those predicted in human, with exactly the same spacing and orientation. Intriguingly, in most of the cases the enhanceosome motifs have 1 bp wider spacing than the corresponding dimers predicted genome-wide, which is likely caused by the overall 3D structure constraints of the enhanceosome-bound complex.

Electronic supplementary material

The online version of this article (doi:10.1186/s12918-016-0258-3) contains supplementary material, which is available to authorized users.

Cooperative DNA Recognition Modulated by an Interplay between Protein-Protein Interactions and DNA-Mediated Allostery

Highly specific transcriptional regulation depends on the cooperative association of transcription factors into enhanceosomes. Usually, their DNA-binding cooperativity originates from either direct interactions or DNA-mediated allostery. Here, we performed unbiased molecular simulations followed by simulations of protein-DNA unbinding and free energy profiling to study the cooperative DNA recognition by OCT4 and SOX2, key components of enhanceosomes in pluripotent cells. We found that SOX2 influences the orientation and dynamics of the DNA-bound configuration of OCT4. In addition SOX2 modifies the unbinding free energy profiles of both DNA-binding domains of OCT4, the POU specific and POU homeodomain, despite interacting directly only with the first. Thus, we demonstrate that the OCT4-SOX2 cooperativity is modulated by an interplay between protein-protein interactions and DNA-mediated allostery. Further, we estimated the change in OCT4-DNA binding free energy due to the cooperativity with SOX2, observed a good agreement with experimental measurements, and found that SOX2 affects the relative DNA-binding strength of the two OCT4 domains. Based on these findings, we propose that available interaction partners in different biological contexts modulate the DNA exploration routes of multi-domain transcription factors such as OCT4. We consider the OCT4-SOX2 cooperativity as a paradigm of how specificity of transcriptional regulation is achieved through concerted modulation of protein-DNA recognition by different types of interactions.

Pluripotent stem cells can give rise to all somatic lineages. When taken out of the context of the embryo they can be maintained and for this a core transcriptional regulatory circuitry is crucial. OCT4 and SOX2, two factors of this network, are also critical for the induction of pluripotency in somatic cells. In pluripotent cells, OCT4 and SOX2 associate on DNA regulatory regions, enhancing or modifying each other's sequence specificity. In contrast, in the early stages during induction of pluripotency, it was proposed that OCT4 explores the genome independent of SOX2. Here we report the mechanism by which SOX2 influences the orientation, dynamics, and unbinding free energy profile of OCT4. This involves an interplay of protein-protein interactions and DNA-mediated allostery. We consider that this mechanism enables OCT4 to use its DNA binding domains and the interaction partners available in a certain biological context to access alternative genome exploration routes. This study enhances the understanding of the context specific function of OCT4 and provides a general perspective on how DNA-binding cooperativity is modulated by different types of interactions.

Structure-aided prediction of mammalian transcription factor complexes in conserved non-coding elements

Mapping the DNA-binding preferences of transcription factor (TF) complexes is critical for deciphering the functions of cis-regulatory elements. Here, we developed a computational method that compares co-occurring motif spacings in conserved versus unconserved regions of the human genome to detect evolutionarily constrained binding sites of rigid TF complexes. Structural data were used to estimate TF complex physical plausibility, explore overlapping motif arrangements seldom tackled by non-structure-aware methods, and generate and analyse three-dimensional models of the predicted complexes bound to DNA. Using this approach, we predicted 422 physically realistic TF complex motifs at 18% false discovery rate, the majority of which (326, 77%) contain some sequence overlap between binding sites. The set of mostly novel complexes is enriched in known composite motifs, predictive of binding site configurations in TF-TF-DNA crystal structures, and supported by ChIP-seq datasets. Structural modelling revealed three cooperativity mechanisms: direct protein-protein interactions, potentially indirect interactions and 'through-DNA' interactions. Indeed, 38% of the predicted complexes were found to contain four or more bases in which TF pairs appear to synergize through overlapping binding to the same DNA base pairs in opposite grooves or strands. Our TF complex and associated binding site predictions are available as a web resource at http://bejerano.stanford.edu/complex.

Combinatorial control of gene expression

The complexity and diversity of eukaryotic organisms are a feat of nature's engineering. Pulling the strings of such an intricate machinery requires an even more masterful and crafty approach. Only the number and type of responses that they generate exceed the staggering proportions of environmental signals perceived and processed by eukaryotes. Hence, at first glance, the cell's sparse stockpile of controlling factors does not seem remotely adequate to carry out this response. The question as to how eukaryotes sense and respond to environmental cues has no single answer. It is an amalgamation, an interplay between several processes, pathways, and factors—a combinatorial control. A short description of some of the most important elements that operate this entire conglomerate is given in this paper.

Taking into account nucleosomes for predicting gene expression

The eukaryotic genome is organized in a chain of nucleosomes that consist of 145-147 bp of DNA wrapped around a histone octamer protein core. Binding of transcription factors (TF) to nucleosomal DNA is frequently impeded, which makes it a challenging task to calculate TF occupancy at a given regulatory genomic site for predicting gene expression. Here, we review methods to calculate TF binding to DNA in the presence of nucleosomes. The main theoretical problems are (i) the computation speed that is becoming a bottleneck when partial unwrapping of DNA from the nucleosome is considered, (ii) the perturbation of the binding equilibrium by the activity of ATP-dependent chromatin remodelers, which translocate nucleosomes along the DNA, and (iii) the model parameterization from high-throughput sequencing data and fluorescence microscopy experiments in living cells. We discuss strategies that address these issues to efficiently compute transcription factor binding in chromatin.

Disentangling the many layers of eukaryotic transcriptional regulation

Regulation of gene expression in eukaryotes is an extremely complex process. In this review, we break down several critical steps, emphasizing new data and techniques that have expanded current gene regulatory models. We begin at the level of DNA sequence where cis-regulatory modules (CRMs) provide important regulatory information in the form of transcription factor (TF) binding sites. In this respect, CRMs function as instructional platforms for the assembly of gene regulatory complexes. We discuss multiple mechanisms controlling complex assembly, including cooperative DNA binding, combinatorial codes, and CRM architecture. The second section of this review places CRM assembly in the context of nucleosomes and condensed chromatin. We discuss how DNA accessibility and histone modifications contribute to TF function. Lastly, new advances in chromosomal mapping techniques have provided increased understanding of intra- and interchromosomal interactions. We discuss how these topological maps influence gene regulatory models.

Two polymorphic residues account for the differences in DNA binding and transcriptional activation by NF-κB proteins encoded by naturally occurring alleles in Nematostella vectensis

The NF-κB family of transcription factors is activated in response to many environmental and biological stresses, and plays a key role in innate immunity across a broad evolutionary expanse of animals. A simple NF-κB pathway is present in the sea anemone Nematostella vectensis, an important model organism in the phylum Cnidaria. Nematostella has previously been shown to have two naturally occurring NF-κB alleles (Nv-NF-κB-C and Nv-NF-κB-S) that encode proteins with different DNA-binding and transactivation abilities. We show here that polymorphic residues 67 (Cys vs. Ser) and 269 (Ala vs. Glu) play complementary roles in determining the DNA-binding activity of the NF-κB proteins encoded by these two alleles and that residue 67 is primarily responsible for the difference in their transactivation ability. Phylogenetic analysis indicates that Nv-NF-κB-S is the derived allele, consistent with its restricted geographic distribution. These results define polymorphic residues that are important for the DNA-binding and transactivating activities of two naturally occurring variants of Nv-NF-κB. The implications for the appearance of the two Nv-NF-κB alleles in natural populations of sea anemones are discussed.