Monomeric and dimeric bZIP transcription factor GCN4 bind at the same rate to their target DNA site

Domain tethering impacts dimerization and DNA-mediated allostery in the human transcription factor FoxP1

Transcription factors are multidomain proteins with specific DNA binding and regulatory domains. In the human FoxP subfamily (FoxP1, FoxP2, FoxP3, and FoxP4) of transcription factors, a 90 residue-long disordered region links a Leucine Zipper (ZIP)-known to form coiled-coil dimers-and a Forkhead (FKH) domain-known to form domain swapping dimers. We used replica exchange discrete molecular dynamics simulations, single-molecule fluorescence experiments, and other biophysical tools to understand how domain tethering in FoxP1 impacts dimerization at ZIP and FKH domains and how DNA binding allosterically regulates their dimerization. We found that domain tethering promotes FoxP1 dimerization but inhibits a FKH domain-swapped structure. Furthermore, our findings indicate that the linker mediates the mutual organization and dynamics of ZIP and FKH domains, forming closed and open states with and without interdomain contacts, thus highlighting the role of the linkers in multidomain proteins. Finally, we found that DNA allosterically promotes structural changes that decrease the dimerization propensity of FoxP1. We postulate that, upon DNA binding, the interdomain linker plays a crucial role in the gene regulatory function of FoxP1.

Human FoxP Transcription Factors as Tractable Models of the Evolution and Functional Outcomes of Three-Dimensional Domain Swapping

The association of two or more proteins to adopt a quaternary complex is one of the most widespread mechanisms by which protein function is modulated. In this scenario, three-dimensional domain swapping (3D-DS) constitutes one plausible pathway for the evolution of protein oligomerization that exploits readily available intramolecular contacts to be established in an intermolecular fashion. However, analysis of the oligomerization kinetics and thermodynamics of most extant 3D-DS proteins shows its dependence on protein unfolding, obscuring the elucidation of the emergence of 3D-DS during evolution, its occurrence under physiological conditions, and its biological relevance. Here, we describe the human FoxP subfamily of transcription factors as a feasible model to study the evolution of 3D-DS, due to their significantly faster dissociation and dimerization kinetics and lower dissociation constants in comparison to most 3D-DS models. Through the biophysical and functional characterization of FoxP proteins, relevant structural aspects highlighting the evolutionary adaptations of these proteins to enable efficient 3D-DS have been ascertained. Most biophysical studies on FoxP suggest that the dynamics of the polypeptide chain are crucial to decrease the energy barrier of 3D-DS, enabling its fast oligomerization under physiological conditions. Moreover, comparison of biophysical parameters between human FoxP proteins in the context of their minute sequence differences suggests differential evolutionary strategies to favor homoassociation and presages the possibility of heteroassociations, with direct impacts in their gene regulation function.

Group-C/S1 bZIP heterodimers regulate MdIPT5b to negatively modulate drought tolerance in apple species

Cytokinins play a central role in delaying senescence, reducing oxidative damage and maintaining plant growth during drought. This study showed that the ectopic expression of ProRE-deleted MdIPT5b, a key enzyme involved in cytokinin metabolism, increased the drought tolerance of transgenic Malus domestica (apple) callus and Solanum lycopersicum (tomato) seedlings by maintaining cytokinin homeostasis, and thus maintaining redox balance. Under restricted watering regimes, the yields of transgenic tomato plants were enhanced. Heterodimers of C/S1 bZIP are involved in the cytokinin-mediated drought response. The heterodimers bind the ProRE of MdIPT5b promoter, thus directly suppressing gene transcription. Single C/S1 bZIP members could not independently function as suppressors. However, specific paired members (heterodimers of MdbZIP80 with MdbZIP2 or with MdbZIP39) effectively suppressed transcription. The α-helical structure is essential for the heterodimerization of C/S1 bZIP members and for synergistic transcriptional suppression. As negative regulators of drought tolerance, suppressing either MdbZIP2 or MdbZIP39 alone does not improve the expression of MdIPT5b and did not increase the drought tolerance of transgenic apple callus. However, this could be achieved when they were co-suppressed. The suppression of MdbZIP80 alone could improve MdIPT5b expression and increase the drought tolerance of transgenic apple callus. However, these effects were reversed in response to the cosuppression of MdbZIP80 and MdIPT5b. Similar results were also observed during delayed dark-induced senescence in apple leaves. In conclusion, the apple C/S1 bZIP network (involving MdbZIP2, MdbZIP39 and MdbZIP80) directly suppressed the expression of MdIPT5b, thus negatively modulating drought tolerance and dark-induced senescence in a functionally redundant manner.

Selective antagonism of cJun for cancer therapy

The activator protein-1 (AP-1) family of transcription factors modulate a diverse range of cellular signalling pathways into outputs which can be oncogenic or anti-oncogenic. The transcription of relevant genes is controlled by the cellular context, and in particular by the dimeric composition of AP-1. Here, we describe the evidence linking cJun in particular to a range of cancers. This includes correlative studies of protein levels in patient tumour samples and mechanistic understanding of the role of cJun in cancer cell models. This develops an understanding of cJun as a focal point of cancer-altered signalling which has the potential for therapeutic antagonism. Significant work has produced a range of small molecules and peptides which have been summarised here and categorised according to the binding surface they target within the cJun-DNA complex. We highlight the importance of selectively targeting a single AP-1 family member to antagonise known oncogenic function and avoid antagonism of anti-oncogenic function.

Sequence-Specific DNA Binding by Noncovalent Peptide-Azocyclodextrin Dimer Complex as a Suitable Model for Conformational Fuzziness

Transcription factors are proteins lying at the endpoint of signaling pathways that control the complex process of DNA transcription. Typically, they are structurally disordered in the inactive state, but in response to an external stimulus, like a suitable ligand, they change their conformation, thereby activating DNA transcription in a spatiotemporal fashion. The observed disorder or fuzziness is functionally beneficial because it can add adaptability, versatility, and reversibility to the interaction. In this context, mimetics of the basic region of the GCN4 transcription factor (Tf) and their interaction with dsDNA sequences would be suitable models to explore the concept of conformational fuzziness experimentally. Herein, we present the first example of a system that mimics the DNA sequence-specific recognition by the GCN4 Tf through the formation of a non- covalent tetra-component complex: peptide–azoβ-CyD(dimer)–peptide–DNA. The non-covalent complex is constructed on the one hand by a 30 amino acid peptide corresponding to the basic region of GCN4 and functionalized with an adamantane moiety, and on the other hand an allosteric receptor, the azoCyDdimer, that has an azobenzene linker connecting two β-cyclodextrin units. The azoCyDdimer responds to light stimulus, existing as two photo-states: the first thermodynamically stable with an E:Z isomer ratio of 95:5 and the second obtained after irradiation with ultraviolet light, resulting in a photostationary state with a 60:40 E:Z ratio. Through electrophoretic shift assays and circular dichroism spectroscopy, we demonstrate that the E isomer is responsible for dimerization and recognition. The formation of the non-covalent tetra component complex occurs in the presence of the GCN4 cognate dsDNA sequence (′5-..ATGA cg TCAT..-3′) but not with (′5-..ATGA c TCAT..-3′) that differs in only one spacing nucleotide. Thus, we demonstrated that the tetra-component complex is formed in a specific manner that depends on the geometry of the ligand, the peptide length, and the ds DNA sequence. We hypothesized that the mechanism of interaction is sequential, and it can be described by the polymorphism model of static fuzziness. We argue that chemically modified peptides of the GCN4 Tf are suitable minimalist experimental models to investigate conformational fuzziness in protein–DNA interactions.

Insights from the crystal structure of the chicken CREB3 bZIP suggest that members of the CREB3 subfamily transcription factors may be activated in response to oxidative stress

cAMP response element binding Protein 3 (CREB3) is an endoplasmic reticulum (ER) membrane‐bound transcription factor, which belongs to the basic leucine zipper (bZIP) superfamily of eukaryotic transcription factors. CREB3 plays a role in the ER‐stress induced unfolded protein response (UPR) and is a multifunctional cellular factor implicated in a number of biological processes including cell proliferation and migration, tumor suppression, and immune‐related gene expression. To gain structural insights into the transcription factor, we determined the crystal structure of the conserved bZIP domain of chicken CREB3 (chCREB3) to a resolution of 3.95 Å. The X‐ray structure provides evidence that chCREB3 can form a stable homodimer. The chCREB3 bZIP has a structured, pre‐formed DNA binding region, even in the absence of DNA, a feature that could potentially enhance both the DNA binding specificity and affinity of chCREB3. Significantly, the homodimeric bZIP possesses an intermolecular disulfide bond that connects equivalent cysteine residues of the parallel helices in the leucine zipper region. This disulfide bond in the hydrophobic core of the bZIP may increase the stability of the homodimer under oxidizing conditions. Moreover, sequence alignment of bZIP sequences from chicken, human, and mouse reveals that only members of the CREB3 subfamily contain this cysteine residue, indicating that it could act as a redox‐sensor. Taken together, these results suggest that the activity of these transcription factors may be redox‐regulated and they may be activated in response to oxidative stress.

PDB Code(s): 6IAK

Quantitative analyses of the equilibria among DNA complexes of a blue-light-regulated bZIP module, Photozipper

Aureochrome1 is a blue-light-receptor protein identified in a stramenopile alga, Vaucheria frigida. Photozipper (PZ) is an N-terminally truncated, monomeric, V. frigida aureochrome1 fragment containing a basic leucine zipper (bZIP) domain and a light–oxygen–voltage (LOV)-sensing domain. PZ dimerizes upon photoexcitation and consequently increases its affinity for the target sequence. In the present study, to understand the equilibria among DNA complexes of PZ, DNA binding by PZ and mutational variants was quantitatively investigated by electrophoretic-mobility-shift assay and fluorescence-correlation spectroscopy in the dark and light states. DNA binding by PZ was sequence-specific and light-dependent. The half-maximal effective concentration of PZ for binding to the target DNA sequence was ~40 nM in the light, which was >10-fold less than the value in the dark. By contrast, the dimeric PZ-S₂C variant (with intermolecular disulfide bonds) had higher affinity for the target sequence, with dissociation constants of ~4 nM, irrespective of the light conditions. Substitutions of Glu159 and Lys164 in the leucine zipper region decreased the affinity of PZ for the target sequence, especially in the light, suggesting that these residues form inter-helical salt bridges between leucine zipper regions, stabilizing the dimer–DNA complex. Our quantitative analyses of the equilibria in PZ–DNA-complex formation suggest that the blue-light-induced dimerization of LOV domains and coiled-coil formation by leucine zipper regions are the primary determinants of the affinity of PZ for the target sequence.

A-ZIP53, a dominant negative reveals the molecular mechanism of heterodimerization between bZIP53, bZIP10 and bZIP25 involved in Arabidopsis seed maturation

In Arabidopsis, maturation phase, an intricate process in seed formation is tightly regulated by the DNA binding activity of protagonist basic leucine zipper 53 (bZIP53) transcription factor and its heterodimerizing partners, bZIP10 and bZIP25. Structural determinants responsible for heterodimerization specificity of bZIP53 are poorly understood. Analysis of amino acid sequences of three bZIPs does not identify interactions that may favor heterodimerization. Here, we describe a designed dominant negative termed A-ZIP53 that has a glutamic acid-rich amphipathic peptide sequence attached to N-terminal of bZIP53 leucine zipper. Circular dichroism (CD) and mass spectrometry studies with equimolar mixture of three bZIP proteins in pairs showed no heterodimer formation whereas A-ZIP53 interacted and formed stable heterodimers with bZIP53, bZIP10, and bZIP25. A-ZIP53 electrostatically mimics DNA and can overcome repulsion between basic DNA binding regions of three bZIP proteins. Gel shift experiments showed that A-ZIP53 can inhibit the DNA binding of three proteins. CD studies demonstrated the specificity of A-ZIP53 as it did not interact with bZIP39 and bZIP72. Transient co-transfections in Arabidopsis protoplasts showed that A-ZIP53 inhibited three bZIPs and their putative heterodimers-mediated transactivation of GUS reporter gene. Furthermore, four newly designed acidic extensions were evaluated for their ability to interact with three bZIPs.

Facilitated diffusion in the presence of obstacles on the DNA

Biological functions of DNA depend on the sequence-specific binding of DNA-binding proteins to their corresponding binding sites. Binding of these proteins to their binding sites occurs through a facilitated diffusion process that combines three-dimensional diffusion in the cytoplasm with one-dimensional diffusion (sliding) along the DNA. In this work, we use a lattice model of facilitated diffusion to study how the dynamics of binding of a protein to a specific site (e.g., binding of an RNA polymerase to a promoter or of a transcription factor to its operator site) is affected by the presence of other proteins bound to the DNA, which act as 'obstacles' in the sliding process. Different types of these obstacles with different dynamics are implemented. While all types impair facilitated diffusion, the extent of the hindrance depends on the type of obstacle. As a consequence of hindrance by obstacles, more excursions into the cytoplasm are required for optimal target binding compared to the case without obstacles.

The DNA target determines the dimerization partner selected by bHLHZ-like hybrid proteins AhRJun and ArntFos

The molecular basis of protein–partner selection and DNA binding of the basic helix–loop–helix (bHLH) and basic region-leucine zipper (bZIP) superfamilies of dimeric transcription factors is fundamental toward understanding gene regulation.

Mechanism of allele specific assembly and disruption of master regulator transcription factor complexes of NF-KBp50, NF-KBp65 and HIF1a on a non-coding FAS SNP

A challenging question in genetics is to understand the molecular function of non-coding variants of the genome. By using differential EMSA, ChIP and functional genome analysis, we have found that changes in transcription factors (TF) apparent binding affinity and dissociation rates are responsible for allele specific assembly or disruption of master TFs: we observed that NF-KBp50, NF-KBp65 and HIF1a bind with an affinity of up to 10 fold better to the C-allele than to the T-allele of rs7901656 both in vivo and in vitro. Furthermore, we showed that NF-KBp50, p65 and HIF1a form higher order heteromultimeric complexes overlapping rs7901656, implying synergism of action among TFs governing cellular response to infection and hypoxia. With rs7901656 on the FAS gene as a paradigm, we show how allele specific transcription factor complex assembly and disruption by a causal variant contributes to disease and phenotypic diversity. This finding provides the highly needed mechanistic insight into how the molecular etiology of regulatory SNPs can be understood in functional terms.

Retroactivity effects dependency on the transcription factors binding mechanisms

Downstream connection effects on transcription are caused by retroactivity. When biomolecular dynamical systems interconnect retroactivity is a property that becomes important. The biological functional meaning of these effects is increasingly becoming an area of interest. Downstream targets, which are operator binding sites in transcriptional networks, may induce behaviors such as ultrasensitive responses or even represent an undesired issue in regulation. To the best of our knowledge, the role of the binding mechanisms of transcription factors in relation to minimizing - or enhancing - retroactivity effects has not been previously addressed. Our aim is to evaluate retroactivity effects considering how the binding mechanism impacts the number of free functional transcription factor (FFTF) molecules using a simple model via deterministic and stochastic simulations. We study four transcription factor binding mechanisms (BM): simple monomer binding (SMB), dimer binding (DB), cooperative sequential binding (CSB) and cooperative sequential binding with dimerization (CSB_D). We consider weak and strong binding regimes for each mechanism, where we contrast the cases when the FFTF is bound or unbound to the downstream loads. Upon interconnection, the number of FFTF molecules changed less for the SMB mechanism while for DB they changed the most. Our results show that for the chosen mechanisms (in terms of the corresponding described dynamics), retroactivity effects depend on transcription binding mechanisms. This contributes to the understanding of how the transcription factor regulatory function-such as decision making-and its dynamic needs for the response, may determine the nature of the selected binding mechanism.

Time-Resolved Detection of Light-Induced Dimerization of Monomeric Aureochrome-1 and Change in Affinity for DNA

Aureochrome (Aureo) is a recently discovered blue light sensor protein initially from Vaucheria frigida, in which it controls blue light-dependent branch formation and/or development of a sex organ by a light-dependent change in the affinity for DNA. Although photochemical reactions of Aureo-LOV (LOV is a C-terminal light-oxygen-voltage domain) and the N-terminal truncated construct containing a bZIP (N-terminal basic leucine zipper domain) and a LOV domain have previously been reported, the reaction kinetics of the change in affinity for DNA have never been elucidated. The reactions of Aureo where the cysteines are replaced by serines (AureoCS) as well as the kinetics of the change in affinity for a target DNA are investigated in the time-domain. The dimerization rate constant is obtained as 2.8 × 10(4) M(-1) s(-1), which suggests that the photoinduced dimerization occurs in the LOV domain and the bZIP domain dimerizes using the interaction with DNA. Surprisingly, binding with the target DNA is completed very quickly, 7.7 × 10(4) M(-1) s(-1), which is faster than the protein dimerization rate. It is proposed that the nonspecific electrostatic interaction, which is observed as a weak binding with DNA, may play a role in the efficient searching for the target sequence within the DNA.

Tunable Helicity, Stability and DNA-Binding Properties of Short Peptides with Hybrid Metal Coordination Motifs

Given the prevalent role of α-helical motifs on protein surfaces in mediating protein-protein and protein-DNA interactions, there have been significant efforts to develop strategies to induce α-helicity in short, unstructured peptides to interrogate such interactions. Toward this goal, we have recently introduced hybrid metal coordination motifs (HCMs). HCMs combine a natural metal-binding amino acid side chain with a synthetic chelating group that are appropriately positioned in a peptide sequence to stabilize an α-helical conformation upon metal coordination. Here, we present a series of short peptides modified with HCMs consisting of a His and a phenanthroline group at i and i+7 positions that can induce α-helicity in a metal-tunable fashion as well as direct the formation of discrete dimeric architectures for recognition of biological targets. We show that the induction of α-helicity can be further modulated by secondary sphere interactions between amino acids at the i+4 position and the HCM. A frequently cited drawback of the use of peptides as therapeutics is their propensity to be quickly digested by proteases; here, we observe an enhancement of up to ∼100-fold in the half-lifes of the metal-bound HCM-peptides in the presence of trypsin. Finally, we show that an HCM-bearing peptide sequence, which contains the DNA-recognition domain of a bZIP protein but is devoid of the obligate dimerization domain, can dimerize with the proper geometry and in an α-helical conformation to bind a cognate DNA sequence with high affinities (K_d≥ 65 nM), again in a metal-tunable manner.

Inferring diffusion dynamics from FCS in heterogeneous nuclear environments

Fluorescence correlation spectroscopy (FCS) is a noninvasive technique that probes the diffusion dynamics of proteins down to single-molecule sensitivity in living cells. Critical mechanistic insight is often drawn from FCS experiments by fitting the resulting time-intensity correlation function, G(t), to known diffusion models. When simple models fail, the complex diffusion dynamics of proteins within heterogeneous cellular environments can be fit to anomalous diffusion models with adjustable anomalous exponents. Here, we take a different approach. We use the maximum entropy method to show-first using synthetic data-that a model for proteins diffusing while stochastically binding/unbinding to various affinity sites in living cells gives rise to a G(t) that could otherwise be equally well fit using anomalous diffusion models. We explain the mechanistic insight derived from our method. In particular, using real FCS data, we describe how the effects of cell crowding and binding to affinity sites manifest themselves in the behavior of G(t). Our focus is on the diffusive behavior of an engineered protein in 1) the heterochromatin region of the cell's nucleus as well as 2) in the cell's cytoplasm and 3) in solution. The protein consists of the basic region-leucine zipper (BZip) domain of the CCAAT/enhancer-binding protein (C/EBP) fused to fluorescent proteins.

Kinetically-defined component actions in gene repression

Gene repression by transcription factors, and glucocorticoid receptors (GR) in particular, is a critical, but poorly understood, physiological response. Among the many unresolved questions is the difference between GR regulated induction and repression, and whether transcription cofactor action is the same in both. Because activity classifications based on changes in gene product level are mechanistically uninformative, we present a theory for gene repression in which the mechanisms of factor action are defined kinetically and are consistent for both gene repression and induction. The theory is generally applicable and amenable to predictions if the dose-response curve for gene repression is non-cooperative with a unit Hill coefficient, which is observed for GR-regulated repression of AP1LUC reporter induction by phorbol myristate acetate. The theory predicts the mechanism of GR and cofactors, and where they act with respect to each other, based on how each cofactor alters the plots of various kinetic parameters vs. cofactor. We show that the kinetically-defined mechanism of action of each of four factors (reporter gene, p160 coactivator TIF2, and two pharmaceuticals [NU6027 and phenanthroline]) is the same in GR-regulated repression and induction. What differs is the position of GR action. This insight should simplify clinical efforts to differentially modulate factor actions in gene induction vs. gene repression.

A modular analysis of the auxin signalling network

Auxin is essential for plant development from embryogenesis onwards. Auxin acts in large part through regulation of transcription. The proteins acting in the signalling pathway regulating transcription downstream of auxin have been identified as well as the interactions between these proteins, thus identifying the topology of this network implicating 54 Auxin Response Factor (ARF) and Aux/IAA (IAA) transcriptional regulators. Here, we study the auxin signalling pathway by means of mathematical modeling at the single cell level. We proceed analytically, by considering the role played by five functional modules into which the auxin pathway can be decomposed: the sequestration of ARF by IAA, the transcriptional repression by IAA, the dimer formation amongst ARFs and IAAs, the feedback loop on IAA and the auxin induced degradation of IAA proteins. Focusing on these modules allows assessing their function within the dynamics of auxin signalling. One key outcome of this analysis is that there are both specific and overlapping functions between all the major modules of the signaling pathway. This suggests a combinatorial function of the modules in optimizing the speed and amplitude of auxin-induced transcription. Our work allows identifying potential functions for homo- and hetero-dimerization of transcriptional regulators, with ARF:IAA, IAA:IAA and ARF:ARF dimerization respectively controlling the amplitude, speed and sensitivity of the response and a synergistic effect of the interaction of IAA with transcriptional repressors on these characteristics of the signaling pathway. Finally, we also suggest experiments which might allow disentangling the structure of the auxin signaling pathway and analysing further its function in plants.

Exploiting anthracene photodimerization within peptides: light induced sequence-selective DNA binding

Here we detail the first example of anthracene photodimerisation in peptides, and use it to trigger a selective biomolecular recognition event.

Elucidation of the specific formation of homo- and heterodimeric forms of ThbZIP1 and its role in stress

Protein–protein interactions are important for the molecular understanding of the biological processes of proteins. The dimerization of bZIPs (basic leucine zipper proteins) is involved in modifying binding site specificities, altering dimer stability, and permitting a new set of specific protein-to-protein interactions to occur at the promoter. In the present study, we studied the whether ThbZIP1 form homo- and heterodimers using the yeast two-hybrid method. Five bZIP genes were cloned from Tamarix hispida to investigate their interaction with ThbZIP1. Our results showed that ThbZIP1 can form homodimers with itself, and three out of five bZIPs could interact with the ThbZIP1 protein to form heterodimers. Real-time RT-PCR results suggested that these ThbZIPs can all respond to abiotic stresses and abscisic acid (ABA), and shared very similar expression patterns in response to NaCl, ABA or PEG6000. Subcellular localization studies showed that all ThbZIPs are targeted to the nucleus. Our results showed that ThbZIP1 are dimeric proteins, which can form homo- or heterodimers.

Asymmetric DNA-search dynamics by symmetric dimeric proteins

We focus on dimeric DNA-binding proteins from two well-studied families: orthodox type II restriction endonucleases (REs) and transcription factors (TFs). Interactions of the protein's recognition sites with the DNA and, particularly, the contribution of each of the monomers to one-dimensional (1D) sliding along nonspecific DNA were studied using computational tools. Coarse-grained molecular dynamics simulations of DNA scanning by various TFs and REs provide insights into how the symmetry of a homodimer can be broken while they nonspecifically interact with DNA. The characteristics of protein sliding along DNA, such as the average sliding length, partitioning between 1D and 3D search, and the one-dimensional diffusion coefficient D1, strongly depend on the salt concentration, which in turn affects the probability of the two monomers adopting a cooperative symmetric sliding mechanism. Indeed, we demonstrate that maximal DNA search efficiency is achieved when the protein adopts an asymmetric search mode in which one monomer slides while its partner hops. We find that proteins classified as TFs have a higher affinity for the DNA, longer sliding lengths, and an increased probability of symmetric sliding in comparison with REs. Moreover, TFs can perform their biological function over a much wider range of salt concentrations than REs. Our results demonstrate that the different biological functions of DNA-binding proteins are related to the different nonspecific DNA search mechanisms they adopt.

The bZIP dimer localizes at DNA full-sites where each basic region can alternately translocate and bind to subsites at the half-site

Crystal structures of the GCN4 bZIP (basic region/leucine zipper) with the AP-1 or CRE site show how each GCN4 basic region binds to a 4 bp cognate half-site as a single DNA target; however, this may not always fully describe how bZIP proteins interact with their target sites. Previously, we showed that the GCN4 basic region interacts with all 5 bp in half-site TTGCG (termed 5H-LR) and that 5H-LR comprises two 4 bp subsites, TTGC and TGCG, which individually are also target sites of the basic region. In this work, we explore how the basic region interacts with 5H-LR when the bZIP dimer localizes to full-sites. Using AMBER molecular modeling, we simulated GCN4 bZIP complexes with full-sites containing 5H-LR to investigate in silico the interface between the basic region and 5H-LR. We also performed in vitro investigation of bZIP-DNA interactions at a number of full-sites that contain 5H-LR versus either subsite: we analyzed results from DNase I footprinting and electrophoretic mobility shift assay (EMSA) and from EMSA titrations to quantify binding affinities. Our computational and experimental results together support a highly dynamic DNA-binding model: when a bZIP dimer localizes to its target full-site, the basic region can alternately recognize either subsite as a distinct target at 5H-LR and translocate between the subsites, potentially by sliding and hopping. This model provides added insights into how α-helical DNA-binding domains of transcription factors can localize to their gene regulatory sequences in vivo.

Redox-mediated mechanisms regulate DNA binding activity of the G-group of basic region leucine zipper (bZIP) transcription factors in Arabidopsis

Plant genes that contain the G-box in their promoters are responsive to a variety of environmental stimuli. Bioinformatics analysis of transcriptome data revealed that the G-box element is significantly enriched in promoters of high light-responsive genes. From nuclear extracts of high light-treated Arabidopsis plants, we identified the AtbZIP16 transcription factor as a component binding to the G-box-containing promoter fragment of light-harvesting chlorophyll a/b-binding protein2.4 (LHCB2.4). AtbZIP16 belongs to the G-group of Arabidopsis basic region leucine zipper (bZIP) type transcription factors. Although AtbZIP16 and its close homologues AtbZIP68 and AtGBF1 bind the G-box, they do not bind the mutated half-sites of the G-box palindrome. In addition, AtbZIP16 interacts with AtbZIP68 and AtGBF1 in the yeast two-hybrid system. A conserved Cys residue was shown to be necessary for redox regulation and enhancement of DNA binding activity in all three proteins. Furthermore, transgenic Arabidopsis lines overexpressing the wild type version of bZIP16 and T-DNA insertion mutants for bZIP68 and GBF1 demonstrated impaired regulation of LHCB2.4 expression. Finally, overexpression lines for the mutated Cys variant of bZIP16 provided support for the biological significance of Cys(330) in redox regulation of gene expression. Thus, our results suggest that environmentally induced changes in the redox state regulate the activity of members of the G-group of bZIP transcription factors.

Structural basis of RNA binding by leucine zipper GCN4

Recently, we showed that leucine zipper (LZ) motifs of basic leucine zipper (bZIP) transcription factors GCN4 and c-Jun are capable of catalyzing degradation of RNA (Nikolaev et al., PLoS ONE 2010; 5:e10765). This observation is intriguing given the tight regulation of RNA turnover control and the antiquity of bZIP transcription factors. To support further mechanistic studies, herein, we elucidated RNA binding interface of the GCN4 leucine zipper motif from yeast. Solution NMR experiments showed that the LZ-RNA interaction interface is located in the first two heptads of LZ moiety, and that only the dimeric (coiled coil) LZ conformation is capable of binding RNA. Site-directed mutagenesis of the LZ-GCN4 RNA binding interface showed that substrate binding is facilitated by lysine and arginine side chains, and that at least one nucleophilic residue is located in proximity to the RNA phosphate backbone. Further studies in the context of full-length bZIP factors are envisaged to address the biological relevance of LZ RNase activity.

N-terminal segments modulate the α-helical propensities of the intrinsically disordered basic regions of bZIP proteins

Basic region leucine zippers (bZIPs) are modular transcription factors that play key roles in eukaryotic gene regulation. The basic regions of bZIPs (bZIP-bRs) are necessary and sufficient for DNA binding and specificity. Bioinformatic predictions and spectroscopic studies suggest that unbound monomeric bZIP-bRs are uniformly disordered as isolated domains. Here, we test this assumption through a comparative characterization of conformational ensembles for 15 different bZIP-bRs using a combination of atomistic simulations and circular dichroism measurements. We find that bZIP-bRs have quantifiable preferences for α-helical conformations in their unbound monomeric forms. This helicity varies from one bZIP-bR to another despite a significant sequence similarity of the DNA binding motifs (DBMs). Our analysis reveals that intramolecular interactions between DBMs and eight-residue segments directly N-terminal to DBMs are the primary modulators of bZIP-bR helicities. We test the accuracy of this inference by designing chimeras of bZIP-bRs to have either increased or decreased overall helicities. Our results yield quantitative insights regarding the relationship between sequence and the degree of intrinsic disorder within bZIP-bRs, and might have general implications for other intrinsically disordered proteins. Understanding how natural sequence variations lead to modulation of disorder is likely to be important for understanding the evolution of specificity in molecular recognition through intrinsically disordered regions (IDRs).

Direct measurement of DNA affinity landscapes on a high-throughput sequencing instrument

Several methods for characterizing DNA-protein interactions are available, but none have demonstrated both high throughput and quantitative measurement of affinity. Here we describe 'high-throughput sequencing'-'fluorescent ligand interaction profiling' (HiTS-FLIP), a technique for measuring quantitative protein-DNA binding affinity at unprecedented depth. In this approach, the optics built into a high-throughput sequencer are used to visualize in vitro binding of a protein to sequenced DNA in a flow cell. Application of HiTS-FLIP to the protein Gcn4 (Gcn4p), the master regulator of the yeast amino acid starvation response, yielded ~440 million binding measurements, enabling determination of dissociation constants for all 12-mer sequences having submicromolar affinity. These data revealed a complex interdependency between motif positions, allowed improved discrimination of in vivo Gcn4p binding sites and regulatory targets relative to previous methods and showed that sets of genes with different promoter affinities to Gcn4p have distinct functions and expression kinetics. Broad application of this approach should increase understanding of the interactions that drive transcription.

A kinetic model of TBP auto-regulation exhibits bistability

Background

TATA Binding Protein (TBP) is required for transcription initiation by all three eukaryotic RNA polymerases. It participates in transcriptional initiation at the majority of eukaryotic gene promoters, either by direct association to the TATA box upstream of the transcription start site or by indirectly localizing to the promoter through other proteins. TBP exists in solution in a dimeric form but binds to DNA as a monomer. Here, we present the first mathematical model for auto-catalytic TBP expression and use it to study the role of dimerization in maintaining the steady state TBP level.

Results

We show that the autogenous regulation of TBP results in a system that is capable of exhibiting three steady states: an unstable low TBP state, one stable state corresponding to a physiological TBP concentration, and another stable steady state corresponding to unviable cells where no TBP is expressed. Our model predicts that a basal level of TBP is required to establish the transcription of the TBP gene, and hence for cell viability. It also predicts that, for the condition corresponding to a typical mammalian cell, the high-TBP state and cell viability is sensitive to variation in DNA binding strength. We use the model to explore the effect of the dimer in buffering the response to changes in TBP levels, and show that for some physiological conditions the dimer is not important in buffering against perturbations.

Conclusions

Results on the necessity of a minimum basal TBP level support the in vivo observations that TBP is maternally inherited, providing the small amount of TBP required to establish its ubiquitous expression. The model shows that the system is sensitive to variations in parameters indicating that it is vulnerable to mutations in TBP. A reduction in TBP-DNA binding constant can lead the system to a regime where the unviable state is the only steady state. Contrary to the current hypotheses, we show that under some physiological conditions the dimer is not very important in restoring the system to steady state. This model demonstrates the use of mathematical modelling to investigate system behaviour and generate hypotheses governing the dynamics of such nonlinear biological systems.

Reviewers

This article was reviewed by Tomasz Lipniacki, James Faeder and Anna Marciniak-Czochra.

Arc-repressor dimerization on DNA: folding rate enhancement by colocalization

Multimeric proteins are ubiquitous in many cellular processes that require high levels of regulation. Eukaryotic gene expression is often regulated by a mechanism of combinatorial control that involves the binding of dimeric transcription factors to DNA together with the coordinated activity of additional proteins. In this study, we investigated the dimerization of the Arc-repressor on DNA with the aim of achieving microscopic insight into the possible advantages of interacting with DNA as a complex rather than as a monomeric single-domain protein. We used a computational coarse-grained model in which the protein dynamics was governed by native interactions and protein-DNA interactions were dictated by electrostatic forces. Inspired by previous experimental work that showed an enhanced refolding rate for the Arc-repressor in the presence of DNA and other polyanions, we focused on the mechanism and kinetics of the assembly of Arc monomers in the presence of single- (ssDNA) and double-stranded DNA (dsDNA) molecules in a low-salt concentration environment. The electrostatic interactions that attract the protein to the dsDNA were shown to be fundamental in colocalizing the unfolded Arc chains and in accelerating refolding. Arc monomers bind the dsDNA efficiently and nonspecifically, and search for each other via one-dimensional diffusion. The fastest folding of Arc is observed for DNA of 30 bp. Longer DNA is significantly less efficient in accelerating the Arc refolding rate, since the two subunits search distinct regions of the one-dimensional DNA and are therefore much less colocalized. The probability that the two unfolded chains will meet on 200 bp DNA is similar to that in the bulk. The colocalization of Arc subunits on ssDNA results in much faster folding compared to that obtained on dsDNA of the same length. Differences in the rate of Arc refolding, cooperativity, and the structure of its transition state ensemble introduced by ssDNA and dsDNA molecules demonstrate the important role of colocalization in biological self-assembly processes.

Modeling the dynamics of transcriptional gene regulatory networks for animal development

The dynamic process of cell fate specification is regulated by networks of regulatory genes. The architecture of the network defines the temporal order of specification events. To understand the dynamic control of the developmental process, the kinetics of mRNA and protein synthesis and the response of the cis-regulatory modules to transcription factor concentration must be considered. Here we review mathematical models for mRNA and protein synthesis kinetics which are based on experimental measurements of the rates of the relevant processes. The model comprises the response functions of cis-regulatory modules to their transcription factor inputs, by incorporating binding site occupancy and its dependence on biologically measurable quantities. We use this model to simulate gene expression, to distinguish between cis-regulatory execution of "AND" and "OR" logic functions, rationalize the oscillatory behavior of certain transcriptional auto-repressors and to show how linked subcircuits can be dealt with. Model simulations display the effects of mutation of binding sites, or perturbation of upstream gene expression. The model is a generally useful tool for understanding gene regulation and the dynamics of cell fate specification.

Elucidation of the structural determinants responsible for the specific formation of heterodimeric Mxd1/Max b-HLH-LZ and its binding to E-box sequences

The proteins of the Mxd family (formally known as Mad) are antagonists of the oncoprotein c-Myc. They compete with c-Myc for their obligate partner Max to prevent the c-Myc/Max heterodimer from binding to E-box sequences in the target gene promoters. In cancer cells, where Myc is overexpressed, the expression of Mxd proteins is usually insufficient or abrogated. However, the reintroduction of Mxd1 expression in these cells prevents growth and proliferation. While the antagonism of c-Myc functions by Mxd proteins is of potential relevance for the development of cancer treatment strategies, the structural determinants responsible for the specific heterodimerization between the Mxd and the Max b-helix-loop-helix-leucine zippers are not fully understood. Moreover, whether the heterodimer is assembled on DNA or in the nucleoplasm prior to DNA binding is under debate. In this article, we demonstrate that Mxd1 D112a and Max N78a and H81d, which are located in the leucine zippers of the proteins, can dictate the specificity of heterodimerization and whether or not the Mxd1/Max/DNA complex forms. Our results also indicate that additional specific determinants exist in the helix-loop-helix domains of Max and Mxd1. Finally, we provide evidence that heterodimerization must precede DNA binding in vivo.

Heme oxygenase-1 induction by NRF2 requires inactivation of the transcriptional repressor BACH1

Oxidative stress activates the transcription factor NRF2, which in turn binds cis-acting antioxidant response element (ARE) enhancers and induces expression of protective antioxidant genes. In contrast, the transcriptional repressor BACH1 binds ARE-like enhancers in cells naïve to oxidative stress and antagonizes NRF2 binding until it becomes inactivated by pro-oxidants. Here, we describe the dynamic roles of BACH1 and NRF2 in the transcription of the heme oxygenase-1 (HMOX1) gene. HMOX1 induction, elicited by arsenite-mediated oxidative stress, follows inactivation of BACH1 and precedes activation of NRF2. BACH1 repression is dominant over NRF2-mediated HMOX1 transcription and inactivation of BACH1 is a prerequisite for HMOX1 induction. In contrast, thioredoxin reductase 1 (TXNRD1) is regulated by NRF2 but not by BACH1. By comparing the expression levels of HMOX1 with TXNRD1, we show that nuclear accumulation of NRF2 is not necessary for HMOX1 induction; rather, BACH1 inactivation permits NRF2 already present in the nucleus at low basal levels to bind the HMOX1 promoter and elicit HMOX1 induction. Thus, BACH1 confers an additional level of regulation to ARE-dependent genes that reveals a new dimension to the oxidative stress response.

An analytical rate expression for the kinetics of gene transcription mediated by dimeric transcription factors

To model gene transcription kinetics, empirical fitting with the Hill function or S-system is often used. In this study, we derived an analytical expression for gene transcription rates in a manner similar to that developed for enzyme kinetics to describe the kinetics of gene transcription mediated by dimeric transcription factors (TFs) such as Gcn4p, a Saccharomyces cerevisiae master gene regulator. We showed that the analytical rate expression and its parameters estimated from several sets of experimental data could accurately reproduce the experimentally measured promoter-binding activity of Gcn4p. Furthermore, the analytical rate expression allowed us to derive analytically, rather than fit empirically, the parameters of the Hill function and S-system for use in modelling transcription kinetics. We found that a plot of gene transcription rate against Gcn4p concentration gave a sigmoidal dose-response curve with a positive co-operativity Hill coefficient (approximately 1.25), in accordance with previous experimental findings on the promoter binding of dimeric TFs. The characteristics of the dose-response curve around the estimated cellular Gcn4p concentration suggest that transcription regulation is efficiently controlled under physiological conditions. This work is a useful initial step towards analytically modelling and simulating complicated gene transcription networks.

The GCN4 bZIP targets noncognate gene regulatory sequences: quantitative investigation of binding at full and half sites

We previously reported that a basic region/leucine zipper (bZIP) protein, a hybrid of the GCN4 basic region and C/EBP leucine zipper, not only recognizes cognate target sites AP-1 (5'-TGACTCA-3') and cAMP-response element (CRE) (5'-TGACGTCA-3') but also binds selectively to noncognate DNA sites: C/EBP (CCAAT/enhancer binding protein, 5'-TTGCGCAA), XRE1 (xenobiotic response element, 5'-TTGCGTGA), HRE (HIF response element, 5'-GCACGTAG), and E-box (5'-CACGTG). In this work, we used electrophoretic mobility shift assay (EMSA) and circular dichroism (CD) for more extensive characterization of the binding of wt bZIP dimer to noncognate sites as well as full- and half-site derivatives, and we examined changes in flanking sequences. Quantitative EMSA titrations were used to measure dissociation constants of this hybrid, wt bZIP, to DNA duplexes: Full-site binding affinities gradually decrease from cognate sites AP-1 and CRE with Kd values of 13 and 12 nM, respectively, to noncognate sites with Kd values of 120 nM to low microM. DNA-binding selectivity at half sites is maintained; however, half-site binding affinities sharply decrease from the cognate half site (Kd = 84 nM) to noncognate half sites (all Kd values > 2 microM). CD shows that comparable levels of alpha-helical structure are induced in wt bZIP upon binding to cognate AP-1 or noncognate sites. Thus, noncognate sites may contribute to preorganization of stable protein structure before binding target DNA sites. This work demonstrates that the bZIP scaffold may be a powerful tool in the design of small, alpha-helical proteins with desired DNA recognition properties.

Combinatorial control of Arabidopsis proline dehydrogenase transcription by specific heterodimerisation of bZIP transcription factors

Proline metabolism has been implicated in plant responses to abiotic stresses. The Arabidopsis thaliana proline dehydrogenase (ProDH) is catalysing the first step in proline degradation. Transcriptional activation of ProDH by hypo-osmolarity is mediated by an ACTCAT cis element, a typical binding site of basic leucine zipper (bZIP) transcription factors. In this study, we demonstrate by gain-of-function and loss-of-function approaches, as well as chromatin immunoprecipitation (ChIP), that ProDH is a direct target gene of the group-S bZIP factor AtbZIP53. Dimerisation studies making use of yeast and Arabidopsis protoplast-based two-hybrid systems, as well as bimolecular fluorescence complementation (BiFC) reveal that AtbZIP53 does not preferentially form dimers with group-S bZIPs but strongly interacts with members of group-C. In particular, a synergistic interplay of AtbZIP53 and group-C AtbZIP10 was demonstrated by colocalisation studies, strong enhancement of ACTCAT-mediated transcription as well as complementation studies in atbzip53 atbzip10 T-DNA insertion lines. Heterodimer mediated activation of transcription has been found to operate independent of the DNA-binding properties and is described as a crucial mechanism to modulate transcription factor activity and function.

Two-hybrid protein-protein interaction analysis in Arabidopsis protoplasts: establishment of a heterodimerization map of group C and group S bZIP transcription factors

In vivo protein-protein interactions are frequently studied by means of yeast two-hybrid analysis. However, interactions detected in yeast might differ considerably in the plant system. Based on GAL4 DNA-binding (BD) and activation domains (AD) we established an Arabidopsis protoplast two-hybrid (P2H) system. The use of Gateway-compatible vectors enables the high-throughput screening of protein-protein interactions in plant cells. The efficiency of the system was tested by examining the homo- and heterodimerization properties of basic leucine zipper (bZIP) transcription factors. A comprehensive heterodimerization matrix of Arabidopsis thaliana group C and group S bZIP transcription factors was generated by comparing the results of yeast and protoplast two-hybrid experiments. Surprisingly, almost no homodimerization but rather specific and selective heterodimerization was detected. Heterodimers were preferentially formed between group C members (AtbZIP9, -10, -25, -63) and members of group S1 (AtbZIP1, -2, -11, -44, -53). In addition, significant but low-affinity interactions were detected inside group S1, S2 or C AtbZIPs, respectively. As a quantitative approach, P2H identified weak heterodimerization events which were not detected in the yeast system. Thus, in addition to cell biological techniques, P2H is a valuable tool for studying protein-protein interaction in living plant cells.

Modular enzyme design: regulation by mutually exclusive protein folding

A regulatory mechanism is introduced whereupon the catalytic activity of a given enzyme is controlled by ligand binding to a receptor domain of choice. A small enzyme (barnase) and a ligand-binding polypeptide (GCN4) are fused so that a simple topological constraint prevents them from existing simultaneously in their folded states. The two domains consequently engage in a thermodynamic tug-of-war in which the more stable domain forces the less stable domain to unfold. In the absence of ligand, the barnase domain is more stable and is therefore folded and active; the GCN4 domain is substantially unstructured. DNA binding induces folding of GCN4, forcibly unfolding and inactivating the barnase domain. Barnase-GCN4 is thus a "natively unfolded" protein that uses ligand binding to switch between partially folded forms. The key characteristics of each parent protein (catalytic efficiency of barnase, DNA binding affinity and sequence specificity of GCN4) are retained in the chimera. Barnase-GCN4 thus defines a modular approach for assembling enzymes with novel sensor capabilities from a variety of catalytic and ligand binding domains.

Slow Assembly and Disassembly of λ Cro Repressor Dimers

Dimers of Cro are required to recognize operator DNA and repress transcription, but dimerization is weak compared to DNA binding. Fluorophore-conjugated, single-cysteine variants of Cro have been used to investigate the equilibria and kinetics of dimer assembly. Equilibrium distributions of mixed dimers, monitored by fluorescence resonance energy transfer (FRET), confirm that labeled variants have equilibrium dimer dissociation constants in the micromolar concentration range. Subunit exchange experiments yield first order rate constants for dimer dissociation that range from 0.02 s(-1) to 0.04 s(-1). Association rate constants calculated from the ratios of dissociation equilibrium and rate constants range from 0.7x10(4) M(-1) s(-1) to 3x10(4) M(-1) s(-1), depending on the site of the fluorescent label. At nanomolar concentrations of subunits, assembly can be driven by addition of DNA. The bimolecular association rate constants measured under these conditions are not dramatically enhanced, ranging from 7x10(4) M(-1) s(-1) to 9x10(4) M(-1) s(-1). The association rate is second order in protein but independent of DNA concentration between 10 nM and 200 nM. The association of subunits under native conditions is more than four orders of magnitude slower than the fast assembly phase measured previously in refolding experiments, and is unaffected by peptidyl-prolyl isomerases. Stabilization of the folded structure of the protein by residue substitution in Cro F58W or reduced temperature increases the ratio of dimers to monomers and decreases the rate of subunit exchange. These data suggest that native monomers have compact structures with substantial barriers to unfolding and that unfolded or partially folded monomers are the preferred substrates for dimer assembly. Cro binding in vivo may be under kinetic rather than thermodynamic control. The slow assembly of Cro dimers demonstrated here provides a new perspective on the lysis/lysogeny switch of bacteriophage lambda.

DNA-binding domain of GCN4 induces bending of both the ATF/CREB and AP-1 binding sites of DNA

The interaction of proteins with DNA results, in some cases, in DNA bending, and this might have functional importance. However, when the protein-induced bending of DNA is small, its measurement presents a problem. It is shown that the fluorescence resonance energy transfer between fluorophores placed on the ends of the specially designed U-shaped DNA, which contains the DNA-binding sites at its central part, can be successfully used for this purpose. The lever effect of the arms of such U-shaped DNA ensures that the distance between the fluorophores is very sensitive to bending of the central part. Using this technique, it was shown that (i) the AP-1 and ATF/CREB binding sites of GCN4 transcription factor are pre-bent to the same extent (approximately 12 degrees toward the major groove) and (ii) binding of the GCN4 DNA-binding domain (GCN4-bZIP) results in additional bending of both these target sites but to a greater extent at the ATF/CREB site. In total, in the complex with GCN4-bZIP, the ATF/CREB site is bent by (25 +/- 2) degrees and the AP-1 site by (20 +/- 2) degrees toward the minor groove.

DNA-mediated assembly of weakly interacting DNA-binding protein subunits: in vitro recruitment of phage 434 repressor and yeast GCN4 DNA-binding domains

The specificity of DNA-mediated protein assembly was studied in two in vitro systems, based on (i) the DNA-binding domain of bacteriophage 434 repressor cI (amino acid residues 1-69), or (ii) the DNA-binding domain of the yeast transcription factor GCN4, (amino acids 1-34) and their respective oligonucleotide cognates. In vivo, both of these peptides are part of larger protein molecules that also contain dimerization domains, and the resulting dimers recognize cognate palindromic DNA sequences that contain two half-sites of 4 bp each. The dimerization domains were not included in the peptides tested, so in solution-in the presence or absence of non-cognate DNA oligonucleotides-these molecules did not show appreciable dimerization, as determined by pyrene excimer fluorescence spectroscopy and oxidative cross-linking monitored by mass spectrometry. Oligonucleotides with only one 4 bp cognate half-site were able to initiate measurable dimerization, and two half-sites were able to select specific dimers even from a heterogeneous pool of molecules of closely related specificity (such as DNA-binding domains of the 434 repressor and their engineered mutants that mimic the binding helix of the related P22 phage repressor). The fluorescent technique allowed us to separately monitor the unspecific, ionic interaction of the peptides with DNA which produced a roughly similar signal in the case of both cognate and non-cognate oligonucleotides. But in the former case, a concomitant excimer fluorescence signal showed the formation of correctly positioned dimers. The results suggest that DNA acts as a highly specific template for the recruitment of weakly interacting protein molecules that can thus build up highly specific complexes.

Gcn4p and novel upstream activating sequences regulate targets of the unfolded protein response

Eukaryotic cells respond to accumulation of unfolded proteins in the endoplasmic reticulum (ER) by activating the unfolded protein response (UPR), a signal transduction pathway that communicates between the ER and the nucleus. In yeast, a large set of UPR target genes has been experimentally determined, but the previously characterized unfolded protein response element (UPRE), an upstream activating sequence (UAS) found in the promoter of the UPR target gene KAR2, cannot account for the transcriptional regulation of most genes in this set. To address this puzzle, we analyzed the promoters of UPR target genes computationally, identifying as candidate UASs short sequences that are statistically overrepresented. We tested the most promising of these candidate UASs for biological activity, and identified two novel UPREs, which are necessary and sufficient for UPR activation of promoters. A genetic screen for activators of the novel motifs revealed that the transcription factor Gcn4p plays an essential and previously unrecognized role in the UPR: Gcn4p and its activator Gcn2p are required for induction of a majority of UPR target genes during ER stress. Both Hac1p and Gcn4p bind target gene promoters to stimulate transcriptional induction. Regulation of Gcn4p levels in response to changing physiological conditions may function as an additional means to modulate the UPR. The discovery of a role for Gcn4p in the yeast UPR reveals an additional level of complexity and demonstrates a surprising conservation of the signaling circuit between yeast and metazoan cells.

The yeast unfolded protein response activates a large set of target genes, but a characterized element found in the promoter of one target, KAR2, cannot account for most targets. Using computational and experimental methods, the authors identify additional elements, as well a role for GCN4p in the response

How do site-specific DNA-binding proteins find their targets?

Essentially all the biological functions of DNA depend on site-specific DNA-binding proteins finding their targets, and therefore 'searching' through megabases of non-target DNA. In this article, we review current understanding of how this sequence searching is done. We review how simple diffusion through solution may be unable to account for the rapid rates of association observed in experiments on some model systems, primarily the Lac repressor. We then present a simplified version of the 'facilitated diffusion' model of Berg, Winter and von Hippel, showing how non-specific DNA-protein interactions may account for accelerated targeting, by permitting the protein to sample many binding sites per DNA encounter. We discuss the 1-dimensional 'sliding' motion of protein along non-specific DNA, often proposed to be the mechanism of this multiple site sampling, and we discuss the role of short-range diffusive 'hopping' motions. We then derive the optimal range of sliding for a few physical situations, including simple models of chromosomes in vivo, showing that a sliding range of approximately 100 bp before dissociation optimizes targeting in vivo. Going beyond first-order binding kinetics, we discuss how processivity, the interaction of a protein with two or more targets on the same DNA, can reveal the extent of sliding and we review recent experiments studying processivity using the restriction enzyme EcoRV. Finally, we discuss how single molecule techniques might be used to study the dynamics of DNA site-specific targeting of proteins.