Pre-bending of a promoter sequence enhances affinity for the TATA-binding factor

Cell phenotypes can be predicted from propensities of protein conformations

Proteins exist as dynamic conformational ensembles. Here we suggest that the propensities of the conformations can be predictors of cell function. The conformational states that the molecules preferentially visit can be viewed as phenotypic determinants, and their mutations work by altering the relative propensities, thus the cell phenotype. Our examples include (i) inactive state variants harboring cancer driver mutations that present active state-like conformational features, as in K-Ras4BG12V compared to other K-Ras4BG12X mutations; (ii) mutants of the same protein presenting vastly different phenotypic and clinical profiles: cancer and neurodevelopmental disorders; (iii) alterations in the occupancies of the conformational (sub)states influencing enzyme reactivity. Thus, protein conformational propensities can determine cell fate. They can also suggest the allosteric drugs efficiency.

RNA conformational propensities determine cellular activity

Cellular processes are the product of interactions between biomolecules, which associate to form biologically active complexes1. These interactions are mediated by intermolecular contacts, which if disrupted, lead to alterations in cell physiology. Nevertheless, the formation of intermolecular contacts nearly universally requires changes in the conformations of the interacting biomolecules. As a result, binding affinity and cellular activity crucially depend both on the strength of the contacts and on the inherent propensities to form binding-competent conformational states2,3. Thus, conformational penalties are ubiquitous in biology and must be known in order to quantitatively model binding energetics for protein and nucleic acid interactions4,5. However, conceptual and technological limitations have hindered our ability to dissect and quantitatively measure how conformational propensities affect cellular activity. Here we systematically altered and determined the propensities for forming the protein-bound conformation of HIV-1 TAR RNA. These propensities quantitatively predicted the binding affinities of TAR to the RNA-binding region of the Tat protein and predicted the extent of HIV-1 Tat-dependent transactivation in cells. Our results establish the role of ensemble-based conformational propensities in cellular activity and reveal an example of a cellular process driven by an exceptionally rare and short-lived RNA conformational state.

Deciphering the mechanical code of the genome and epigenome

Diverse DNA-deforming processes are impacted by the local mechanical and structural properties of DNA, which in turn depend on local sequence and epigenetic modifications. Deciphering this mechanical code (that is, this dependence) has been challenging due to the lack of high-throughput experimental methods. Here we present a comprehensive characterization of the mechanical code. Utilizing high-throughput measurements of DNA bendability via loop-seq, we quantitatively established how the occurrence and spatial distribution of dinucleotides, tetranucleotides and methylated CpG impact DNA bendability. We used our measurements to develop a physical model for the sequence and methylation dependence of DNA bendability. We validated the model by performing loop-seq on mouse genomic sequences around transcription start sites and CTCF-binding sites. We applied our model to test the predictions of all-atom molecular dynamics simulations and to demonstrate that sequence and epigenetic modifications can mechanically encode regulatory information in diverse contexts.

OUP accepted manuscript

DNA mechanical properties play a critical role in every aspect of DNA-dependent biological processes. Recently a high throughput assay named loop-seq has been developed to quantify the intrinsic bendability of a massive number of DNA fragments simultaneously. Using the loop-seq data, we develop a software tool, DNAcycP, based on a deep-learning approach for intrinsic DNA cyclizability prediction. We demonstrate DNAcycP predicts intrinsic DNA cyclizability with high fidelity compared to the experimental data. Using an independent dataset from in vitro selection for enrichment of loopable sequences, we further verified the predicted cyclizability score, termed C-score, can well distinguish DNA fragments with different loopability. We applied DNAcycP to multiple species and compared the C-scores with available high-resolution chemical nucleosome maps. Our analyses showed that both yeast and mouse genomes share a conserved feature of high DNA bendability spanning nucleosome dyads. Additionally, we extended our analysis to transcription factor binding sites and surprisingly found that the cyclizability is substantially elevated at CTCF binding sites in the mouse genome. We further demonstrate this distinct mechanical property is conserved across mammalian species and is inherent to CTCF binding DNA motif.

Mechanical properties of DNA and DNA nanostructures: comparison of atomistic, Martini and oxDNA models

The flexibility and stiffness of small DNA molecules play a fundamental role ranging from several biophysical processes to nano-technological applications. Here, we estimate the mechanical properties of short double-stranded DNA (dsDNA) with lengths ranging from 12 base-pairs (bp) to 56 bp, paranemic crossover (PX) DNA and hexagonal DNA nanotubes (DNTs) using two widely used coarse-grained models - Martini and oxDNA. To calculate the persistence length (Lp) and the stretch modulus (γ) of the dsDNA, we incorporate the worm-like chain and elastic rod model, while for the DNTs, we implement our previously developed theoretical framework. We compare and contrast all of the results with previously reported all-atom molecular dynamics (MD) simulations and experimental results. The mechanical properties of dsDNA (Lp ∼ 50 nm, γ ∼ 800-1500 pN), PX DNA (γ ∼ 1600-2000 pN) and DNTs (Lp ∼ 1-10 μm, γ ∼ 6000-8000 pN) estimated using the Martini soft elastic network and oxDNA are in very good agreement with the all-atom MD and experimental values, while the stiff elastic network Martini reproduces values of Lp and γ which are an order of magnitude higher. The high flexibility of small dsDNA is also depicted in our calculations. However, Martini models proved inadequate to capture the salt concentration effects on the mechanical properties with increasing salt molarity. oxDNA captures the salt concentration effect on the small dsDNA mechanics. But it is found to be ineffective for reproducing the salt-dependent mechanical properties of DNTs. Also, unlike Martini, the time evolved PX DNA and DNT structures from the oxDNA models are comparable to the all-atom MD simulated structures. Our findings provide a route to study the mechanical properties of DNA and DNA based nanostructures with increased time and length scales and has a remarkable implication in the context of DNA nanotechnology.

DNA mismatches reveal conformational penalties in protein-DNA recognition

Transcription factors recognize specific genomic sequences to regulate complex gene-expression programs. Although it is well-established that transcription factors bind to specific DNA sequences using a combination of base readout and shape recognition, some fundamental aspects of protein-DNA binding remain poorly understood1,2. Many DNA-binding proteins induce changes in the structure of the DNA outside the intrinsic B-DNA envelope. However, how the energetic cost that is associated with distorting the DNA contributes to recognition has proven difficult to study, because the distorted DNA exists in low abundance in the unbound ensemble3-9. Here we use a high-throughput assay that we term SaMBA (saturation mismatch-binding assay) to investigate the role of DNA conformational penalties in transcription factor-DNA recognition. In SaMBA, mismatched base pairs are introduced to pre-induce structural distortions in the DNA that are much larger than those induced by changes in the Watson-Crick sequence. Notably, approximately 10% of mismatches increased transcription factor binding, and for each of the 22 transcription factors that were examined, at least one mismatch was found that increased the binding affinity. Mismatches also converted non-specific sites into high-affinity sites, and high-affinity sites into 'super sites' that exhibit stronger affinity than any known canonical binding site. Determination of high-resolution X-ray structures, combined with nuclear magnetic resonance measurements and structural analyses, showed that many of the DNA mismatches that increase binding induce distortions that are similar to those induced by protein binding-thus prepaying some of the energetic cost incurred from deforming the DNA. Our work indicates that conformational penalties are a major determinant of protein-DNA recognition, and reveals mechanisms by which mismatches can recruit transcription factors and thus modulate replication and repair activities in the cell10,11.

A dynamic view of DNA structure within the nucleosome: Biological implications

Most of eukaryotic cellular DNA is packed in nucleosome core particles (NCPs), in which the DNA (DNA_NCP) is wrapped around histones. The influence of this organization on the intrinsic local dynamics of DNA is largely unknown, in particular because capturing such information from experiments remains notoriously challenging. Given the importance of dynamical properties in DNA functions, we addressed this issue using CHARMM36 MD simulations of a nucleosome containing the NCP positioning 601 sequence and four related free dodecamers. Comparison between DNA_NCP and free DNA reveals a limited impact of the dense DNA-histone interface on correlated motions of dinucleotide constituents and on fluctuations of inter base pair parameters. A characteristic feature intimately associated with the DNA_NCP super-helical path is a set of structural periodicities that includes a marked alternation of regions enriched in backbone BI and BII conformers. This observation led to uncover a convincing correspondence between the sequence effect on BI/BII propensities in both DNA_NCP and free DNA, strengthening the idea that the histone preference for particular DNA sequences relies on those intrinsic structural properties. These results offer for the first time a detailed view of the DNA dynamical behavior within NCP. They show in particular that the DNA_NCP dynamics is substantial enough to preserve the ability to structurally adjust to external proteins, for instance remodelers. Also, fresh structural arguments highlight the relevance of relationships between DNA sequence and structural properties for NCP formation. Overall, our work offers a more rational framework to approach the functional, biological roles of NCP.

The Sox2 transcription factor binds RNA

Certain transcription factors are proposed to form functional interactions with RNA to facilitate proper regulation of gene expression. Sox2, a transcription factor critical for maintenance of pluripotency and neurogenesis, has been found associated with several lncRNAs, although it is unknown whether these interactions are direct or via other proteins. Here we demonstrate that human Sox2 interacts directly with one of these lncRNAs with high affinity through its HMG DNA-binding domain in vitro. These interactions are primarily with double-stranded RNA in a non-sequence specific fashion, mediated by a similar but not identical interaction surface. We further determined that Sox2 directly binds RNA in mouse embryonic stem cells by UV-cross-linked immunoprecipitation of Sox2 and more than a thousand Sox2-RNA interactions in vivo were identified using fRIP-seq. Together, these data reveal that Sox2 employs a high-affinity/low-specificity paradigm for RNA binding in vitro and in vivo.

Some transcription factors have been proposed to functionally interact with RNA to facilitate proper regulation of gene expression. Here the authors demonstrate that human Sox2 interact directly and with high affinity to RNAs through its HMG DNA-binding domain.

Decomposing protein-DNA binding and recognition using simplified protein models

We analyze the role of different physicochemical factors in protein/DNA binding and recognition by comparing the results from all-atom molecular dynamics simulations with simulations using simplified protein models. These models enable us to separate the role of specific amino acid side chains, formal amino acid charges and hydrogen bonding from the effects of the low-dielectric volume occupied by the protein. Comparisons are made on the basis of the conformation of DNA after protein binding, the ionic distribution around the complex and the sequence specificity. The results for four transcription factors, binding in either the minor or major grooves of DNA, show that the protein volume and formal charges, with one exception, play a predominant role in binding. Adding hydrogen bonding and a very small number of key amino acid side chains at the all-atom level yields results in DNA conformations and sequence recognition close to those seen in the reference all-atom simulations.

Protein-DNA interfaces: a molecular dynamics analysis of time-dependent recognition processes for three transcription factors

We have studied the dynamics of three transcription factor-DNA complexes using all-atom, microsecond-scale MD simulations. In each case, the salt bridges and hydrogen bond interactions formed at the protein-DNA interface are found to be dynamic, with lifetimes typically in the range of tens to hundreds of picoseconds, although some interactions, notably those involving specific binding to DNA bases, can be a hundred times longer lived. Depending on the complex studied, this dynamics may or may not lead to the existence of distinct conformational substates. Using a sequence threading technique, it has been possible to determine whether DNA sequence recognition is sensitive or not to such conformational changes, and, in one case, to show that recognition appears to be locally dependent on protein-mediated cation distributions.

Prediction of fine-tuned promoter activity from DNA sequence

The quantitative prediction of transcriptional activity of genes using promoter sequence is fundamental to the engineering of biological systems for industrial purposes and understanding the natural variation in gene expression. To catalyze the development of new algorithms for this purpose, the Dialogue on Reverse Engineering Assessment and Methods (DREAM) organized a community challenge seeking predictive models of promoter activity given normalized promoter activity data for 90 ribosomal protein promoters driving expression of a fluorescent reporter gene. By developing an unbiased modeling approach that performs an iterative search for predictive DNA sequence features using the frequencies of various k-mers, inferred DNA mechanical properties and spatial positions of promoter sequences, we achieved the best performer status in this challenge. The specific predictive features used in the model included the frequency of the nucleotide G, the length of polymeric tracts of T and TA, the frequencies of 6 distinct trinucleotides and 12 tetranucleotides, and the predicted protein deformability of the DNA sequence. Our method accurately predicted the activity of 20 natural variants of ribosomal protein promoters (Spearman correlation r = 0.73) as compared to 33 laboratory-mutated variants of the promoters (r = 0.57) in a test set that was hidden from participants. Notably, our model differed substantially from the rest in 2 main ways: i) it did not explicitly utilize transcription factor binding information implying that subtle DNA sequence features are highly associated with gene expression, and ii) it was entirely based on features extracted exclusively from the 100 bp region upstream from the translational start site demonstrating that this region encodes much of the overall promoter activity. The findings from this study have important implications for the engineering of predictable gene expression systems and the evolution of gene expression in naturally occurring biological systems.

Role of microscopic flexibility in tightly curved DNA

The genetic material in living cells is organized into complex structures in which DNA is subjected to substantial contortions. Here we investigate the difference in structure, dynamics, and flexibility between two topological states of a short (107 base pair) DNA sequence in a linear form and a covalently closed, tightly curved circular DNA form. By employing a combination of all-atom molecular dynamics (MD) simulations and elastic rod modeling of DNA, which allows capturing microscopic details while monitoring the global dynamics, we demonstrate that in the highly curved regime the microscopic flexibility of the DNA drastically increases due to the local mobility of the duplex. By analyzing vibrational entropy and Lipari–Szabo NMR order parameters from the simulation data, we propose a novel model for the thermodynamic stability of high-curvature DNA states based on vibrational untightening of the duplex. This novel view of DNA bending provides a fundamental explanation that bridges the gap between classical models of DNA and experimental studies on DNA cyclization, which so far have been in substantial disagreement.

Coupled binding-bending-folding: The complex conformational dynamics of protein-DNA binding studied by atomistic molecular dynamics simulations

BACKGROUND

Protein-DNA binding often involves dramatic conformational changes such as protein folding and DNA bending. While thermodynamic aspects of this behavior are understood, and its biological function is often known, the mechanism by which the conformational changes occur is generally unclear. By providing detailed structural and energetic data, molecular dynamics simulations have been helpful in elucidating and rationalizing protein-DNA binding.

SCOPE OF REVIEW

This review will summarize recent atomistic molecular dynamics simulations of the conformational dynamics of DNA and protein-DNA binding. A brief overview of recent developments in DNA force fields is given as well.

MAJOR CONCLUSIONS

Simulations have been crucial in rationalizing the intrinsic flexibility of DNA, and have been instrumental in identifying the sequence of binding events, the triggers for the conformational motion, and the mechanism of binding for a number of important DNA-binding proteins.

GENERAL SIGNIFICANCE

Molecular dynamics simulations are an important tool for understanding the complex binding behavior of DNA-binding proteins. With recent advances in force fields and rapid increases in simulation time scales, simulations will become even more important for future studies. This article is part of a Special Issue entitled Recent developments of molecular dynamics.

Micro-minicircle Gene Therapy: Implications of Size on Fermentation, Complexation, Shearing Resistance, and Expression

The minicircle (MC), composed of eukaryotic sequences only, is an interesting approach to increase the safety and efficiency of plasmid-based vectors for gene therapy. In this paper, we investigate micro-MC (miMC) vectors encoding small regulatory RNA. We use a construct encoding a splice-correcting U7 small nuclear RNA, which results in a vector of 650 base pairs (bp), as compared to a conventional 3600 bp plasmid carrying the same expression cassette. Furthermore, we construct miMCs of varying sizes carrying different number of these cassettes. This allows us to evaluate how size influences production, super-coiling, stability and efficiency of the vector. We characterize coiling morphology by atomic force microscopy and measure the resistance to shearing forces caused by an injector device, the Biojector. We compare the behavior of miMCs and plasmids in vitro using lipofection and electroporation, as well as in vivo in mice. We here show that when the size of the miMC is reduced, the formation of dimers and trimers increases. There seems to be a lower size limit for efficient expression. We demonstrate that miMCs are more robust than plasmids when exposed to shearing forces, and that they show extended expression in vivo.

Structural determinants of DNA recognition by plant MADS-domain transcription factors

Plant MADS-domain transcription factors act as key regulators of many developmental processes. Despite the wealth of information that exists about these factors, the mechanisms by which they recognize their cognate DNA-binding site, called CArG-box (consensus CCW₆GG), and how different MADS-domain proteins achieve DNA-binding specificity, are still largely unknown. We used information from in vivo ChIP-seq experiments, in vitro DNA-binding data and evolutionary conservation to address these important questions. We found that structural characteristics of the DNA play an important role in the DNA binding of plant MADS-domain proteins. The central region of the CArG-box largely resembles a structural motif called ‘A-tract’, which is characterized by a narrow minor groove and may assist bending of the DNA by MADS-domain proteins. Periodically spaced A-tracts outside the CArG-box suggest additional roles for this structure in the process of DNA binding of these transcription factors. Structural characteristics of the CArG-box not only play an important role in DNA-binding site recognition of MADS-domain proteins, but also partly explain differences in DNA-binding specificity of different members of this transcription factor family and their heteromeric complexes.

Winding single-molecule double-stranded DNA on a nanometer-sized reel

A molecular system of a nanometer-sized reel was developed from F₁–ATPase, a rotary motor protein. By combination with magnetic tweezers and optical tweezers, single-molecule double-stranded DNA (dsDNA) was wound around the molecular reel. The bending stiffness of dsDNA was determined from the winding tension (0.9–6.0 pN) and the diameter of the wound loop (21.4–8.5 nm). Our results were in good agreement with the conventional worm-like chain model and a persistence length of 54 ± 9 nm was estimated. This molecular reel system offers a new platform for single-molecule study of micromechanics of sharply bent DNA molecules and is expected to be applicable to the elucidation of the molecular mechanism of DNA-associating proteins on sharply bent DNA strands.

A flexible integrative approach based on random forest improves prediction of transcription factor binding sites

Transcription factor binding sites (TFBSs) are DNA sequences of 6–15 base pairs. Interaction of these TFBSs with transcription factors (TFs) is largely responsible for most spatiotemporal gene expression patterns. Here, we evaluate to what extent sequence-based prediction of TFBSs can be improved by taking into account the positional dependencies of nucleotides (NPDs) and the nucleotide sequence-dependent structure of DNA. We make use of the random forest algorithm to flexibly exploit both types of information. Results in this study show that both the structural method and the NPD method can be valuable for the prediction of TFBSs. Moreover, their predictive values seem to be complementary, even to the widely used position weight matrix (PWM) method. This led us to combine all three methods. Results obtained for five eukaryotic TFs with different DNA-binding domains show that our method improves classification accuracy for all five eukaryotic TFs compared with other approaches. Additionally, we contrast the results of seven smaller prokaryotic sets with high-quality data and show that with the use of high-quality data we can significantly improve prediction performance. Models developed in this study can be of great use for gaining insight into the mechanisms of TF binding.

Generation of DNA nanocircles containing mismatched bases

The DNA mismatch repair (MMR) system recognizes and repairs errors that escaped the proofreading function of DNA polymerases. To study molecular details of the MMR mechanism, in vitro biochemical assays require specific DNA substrates carrying mismatches and strand discrimination signals. Current approaches used to generate MMR substrates are time-consuming and/or not very flexible with respect to sequence context. Here we report an approach to generate small circular DNA containing a mismatch (nanocircles). Our method is based on the nicking of PCR products resulting in single-stranded 3' overhangs, which form DNA circles after annealing and ligation. Depending on the DNA template, one can generate mismatched circles containing a single hemimethylated GATC site (for use with the bacterial system) and/or nicking sites to generate DNA circles nicked in the top or bottom strand (for assays with the bacterial or eukaryotic MMR system). The size of the circles varied (323 to 1100 bp), their sequence was determined by the template DNA, and purification of the circles was achieved by ExoI/ExoIII digestion and/or gel extraction. The quality of the nanocircles was assessed by scanning-force microscopy and their suitability for in vitro repair initiation was examined using recombinant Escherichia coli MMR proteins.

Basic statistical recipes for the emergence of biochemical discernment

An essential step towards understanding life would be to identify the very basic mechanisms responsible for the discerning behaviour of living biochemical systems, absent from randomly reacting chemical soups. One intuitively feels that this question goes beyond the particular nature of the biological molecules and should relate to general physical principles. The pre-eminent physicist Ludwig Boltzmann early envisioned life as a struggle for entropy, in concordance with the subsequent principle of self-organization out of equilibrium. Re-examination of elementary steady state biochemical systems from a statistical perspective supports this view and shows that sigmoidal responses arising from microstates elimination, are sufficient to explain innermost characteristics of life, including its capacity to convert random molecular interactions into accurate biological reactions. A primary operating strategy to achieve this goal is the introduction of time-irreversible transitions in molecular state conversion cycles by injection of free energy, which confers decisional capacity to single macromolecules. Selected examples from various fields of molecular biology such as enzymology and gene expression, are provided to show that these non-equilibrium steady state mechanisms remain important in contemporary biochemical systems. But in addition, information archiving allowed the emergence of the time-reversible counterparts of these mechanisms, mediated by evolutionary pre-organized macromolecular complexes capable of generating discernment in a non-dissipative manner.

Pyrosequencing enhancement for better detection limit and sequencing homopolymers

Pyrosequencing is a DNA sequencing technique based on sequencing-by-synthesis enabling rapid and real-time sequence determination. Although ample genomic research has been undertaken using pyrosequencing, the requirement of relatively high amount of DNA template and the difficulty in sequencing the homopolymeric regions limit its key advantages in the applications directing towards clinical research. In this study, we demonstrate that pyrosequencing on homopolymeric regions with 10 identical nucleotides can be successfully performed with optimal amount of DNA (0.3125-5 pmol) immobilized on conventional non-porous Sepharose beads. We also validate that by using porous silica beads, the sequencing signal increased 3.5-folds as compared to that produced from same amount of DNA immobilized on solid Sepharose beads. Our results strongly indicate that with optimized quantity of DNA and suitable solid support, the performance of pyrosequencing on homopolymeric regions and its detection limit has been significantly improved.

Counterion and pH-Mediated Structural Changes in Charged Biopolymer Gels

DNA solutions and gels exhibit a wide range of phenomena, many of which have not yet been fully understood. In the presence of multivalent counterions, attraction between charged DNA strands occurs. Increasing the concentration of multivalent ions leads to a decrease of the osmotic pressure, and a sufficiently high ion concentration results in the precipitation of the polymer. Replacing the monovalent counterions by hydrogen ions (decreasing the pH) also produces a marked decrease of the osmotic pressure, and at low pH a phase transition takes place. In this paper we analyze osmotic swelling pressure measurements and small-angle neutron scattering (SANS) measurements made on chemically cross-linked DNA gels swollen in near physiological salt solutions. The effect of calcium ions is compared with that of decreasing the pH of the equilibrium salt solution. We demonstrate that both the concentration dependence of the osmotic pressure and the SANS response of DNA gels display significant differences in the two cases.

Effects of the nucleoid protein HU on the structure, flexibility, and ring-closure properties of DNA deduced from Monte Carlo simulations

The histone-like HU (heat unstable) protein plays a key role in the organization and regulation of the Escherichia coli genome. The nonspecific nature of HU binding to DNA complicates analysis of the mechanism by which the protein contributes to the looping of DNA. Conventional models of the looping of HU-bound duplexes attribute the changes in biophysical properties of DNA brought about by the random binding of protein to changes in the effective parameters of an ideal helical wormlike chain. Here, we introduce a novel Monte Carlo approach to study the effects of nonspecific HU binding on the configurational properties of DNA directly. We randomly decorated segments of an ideal double-helical DNA with HU molecules that induce the bends and other structural distortions of the double helix find in currently available X-ray structures. We find that the presence of HU at levels approximating those found in the cell reduces the persistence length by roughly threefold compared with that of naked DNA. The binding of protein has particularly striking effects on the cyclization properties of short duplexes, altering the dependence of ring closure on chain length in a way that cannot be mimicked by a simple wormlike model and accumulating at higher-than-expected levels on successfully closed chains. Moreover, the uptake of protein on small minicircles depends on chain length, taking advantage of the HU-induced deformations of DNA structure to facilitate ligation. Circular duplexes with bound HU show much greater propensity than protein-free DNA to exist as negatively supercoiled topoisomers, suggesting a potential role of HU in organizing the bacterial nucleoid. The local bending and undertwisting of DNA by HU, in combination with the number of bound proteins, provide a structural rationale for the condensation of DNA and the observed expression levels of reporter genes in vivo.

Using FRET to measure the angle at which a protein bends DNA: TBP binding a TATA box as a model system

An undergraduate biochemistry laboratory experiment that will teach the technique of fluorescence resonance energy transfer (FRET) while analyzing protein-induced DNA bending is described. The experiment uses the protein TATA binding protein (TBP), which is a general transcription factor that recognizes and binds specific DNA sequences known as TATA boxes. When TBP binds to a TATA box, it bends the DNA. Such bending will be detected using FRET to measure the distance between two fluorophores located on the ends of the DNA. When TBP binds and bends the DNA, the fluorophores move closer together, reflected by an increase in FRET. At the completion of the experiment, three parameters will be determined: 1) the efficiency of the FRET, 2) the end-to-end distance between the fluorophores, and 3) the angle at which TBP bends the DNA. In performing this experiment, students will be introduced to FRET, gain experience in quantitative biophysical measurements, and appreciate how a protein can induce a dramatic change in DNA conformation.

Structure and mechanism for DNA lesion recognition

A fundamental question in DNA repair is how a lesion is detected when embedded in millions to billions of normal base pairs. Extensive structural and functional studies reveal atomic details of DNA repair protein and nucleic acid interactions. This review summarizes seemingly diverse structural motifs used in lesion recognition and suggests a general mechanism to recognize DNA lesion by the poor base stacking. After initial recognition of this shared structural feature of lesions, different DNA repair pathways use unique verification mechanisms to ensure correct lesion identification and removal.

FOXM1, a typical proliferation-associated transcription factor

FOXM1 is a typical proliferation-associated transcription factor: it stimulates proliferation by promoting S-phase entry as well as M-phase entry and is involved in proper execution of mitosis. Accordingly, FOXM1 regulates genes that control G1/S-transition, S-phase progression, G2/M-transition and M-phase progression. Consistently, its expression and its activity are antagonistically regulated by many important proliferation and anti-proliferation signals. Furthermore, FOXM1 is implicated in tumorigenesis and contributes to both tumor initiation and progression. In addition to its function as a conventional transcription factor, FOXM1 transactivates the human c-myc P1 and P2 promoters directly via their TATA-boxes by a new transactivation mechanism, which it also employs for transactivation of the human c-fos, hsp70 and histone H2B/a promoters. This review summarizes the current knowledge on FOXM1, in particular its two different transactivation mechanisms, the regulation of its transcriptional activity by proliferation versus anti-proliferation signals and its function in normal cell cycle progression and tumorigenesis.

DNA bending induced by carbocyclic sugar analogs constrained to the north conformation

DNA bending caused by introduction of carbocyclic sugars constrained to the north conformation was studied, using explicit solvent molecular dynamic (MD) simulations. The native Drew-Dickerson (DD) dodecamer and its three modifications containing north carbocyclic sugars in the 7th (T7*), 8th (T8*) or both 7th and 8th (T7T8*) nucleotide positions were examined. Introduction of the carbocyclic sugar results in A-form conformations for the alpha, beta, chi, zeta, and sugar pucker backbone parameters in the modified nucleotides. Increased steric repulsion between the sugar and its parent base in the modified oligonucleotides impacts the roll and cup dinucleotide step parameters, increasing the bending of the oligomer axis. Increased buckling of the substituted nucleotides disrupts the usual stabilizing base stacking interactions. The level of overall bending depends on the number and position of carbocyclic sugars introduced in the DNA sequence. Single sugar substitutions are unable to induce substantial bending due to the neighboring unmodified nucleotides counterbalancing the distortion. Significant bending can, however, be induced by two consecutive north sugars (T7T8*), which is in agreement with experimental results. The modified oligomers populate a wide range of bend angles, indicating that they maintain flexibility in the bent state. The present results suggest that insertion of carbocyclic sugars into DNA or RNA duplexes can be used to engineer bending of the duplexes without impacting the electrostatic or chemical properties of the phosphodiester backbone, thereby serving as excellent tools for experimental elucidation of nucleic acid structure-function relationships.

TFIIA changes the conformation of the DNA in TBP/TATA complexes and increases their kinetic stability

Eukaryotic mRNA transcription by RNA polymerase II is a highly regulated complex reaction involving numerous proteins. In order to control tissue and promoter specific gene expression, transcription factors must work in concert with each other and with the promoter DNA to form the proper architecture to activate the gene of interest. The TATA binding protein (TBP) binds to TATA boxes in core promoters and bends the TATA DNA. We have used quantitative solution fluorescence resonance energy transfer (FRET) and gel-based FRET (gelFRET) to determine the effect of TFIIA on the conformation of the DNA in TBP/TATA complexes and on the kinetic stability of these complexes. Our results indicate that human TFIIA decreases the angle to which human TBP bends consensus TATA DNA from 104 degrees to 80 degrees when calculated using a two-kink model. The kinetic stability of TBP/TATA complexes was greatly reduced by increasing the KCl concentration from 50 mM to 140 mM, which is more physiologically relevant. TFIIA significantly enhanced the kinetic stability of TBP/TATA complexes, thereby attenuating the effect of higher salt concentrations. We also found that TBP bent non-consensus TATA DNA to a lesser degree than consensus TATA DNA and complexes between TBP and a non-consensus TATA box were kinetically unstable even at 50 mM KCl. Interestingly, TFIIA increased the calculated bend angle and kinetic stability of complexes on a non-consensus TATA box, making them similar to those on a consensus TATA box. Our data show that TFIIA induces a conformational change within the TBP/TATA complex that enhances its stability under both in vitro and physiological salt conditions. Furthermore, we present a refined model for the effect that TFIIA has on DNA conformation that takes into account potential changes in bend angle as well as twist angle.

Regulation of transcription by synthetic DNA-bending agents

Gene expression is regulated by a complex interplay between binding and the three-dimensional arrangement of transcription factors with RNA polymerase and DNA. Previous studies have supported a direct role for DNA bending and conformation in gene expression, which suggests that agents that induce bends in DNA might be able to control gene expression. To test this hypothesis, we examined the effect of triple-helix-forming oligonucleotide (TFO) bending agents on the transcription of luciferase in an in vitro transcriptional/translational system. We find that transcription is regulated only by a TFO that induces a bend in the DNA. Related TFOs that do not induce bends in DNA have no effect on transcription. Reporter expression can be increased by as much as 80 % or decreased by as much as 50 % depending on the phasing of the upstream bend relative to the promoter. We interpret the results as follows: when the bend is positioned such that the upstream DNA is curved toward the RNA polymerase on the same DNA face, transcription is enhanced. When the upstream DNA is curved away, transcription is attenuated. These results support the hypothesis that DNA-bending agents might have the capability to regulate gene expression, thereby opening up a previously undervalued avenue in research on the artificial control of gene expression.

The Saccharomyces cerevisiae RNA polymerase III recruitment factor subunits Brf1 and Bdp1 impose a strict sequence preference for the downstream half of the TATA box

Association of the TATA-binding protein (TBP) with its cognate site within eukaryotic promoters is key to accurate and efficient transcriptional initiation. To achieve recruitment of Saccharomyces cerevisiae RNA polymerase III, TBP is associated with two additional factors, Brf1 and Bdp1, to form the initiation factor TFIIIB. Previous data have suggested that the structure or dynamics of the TBP–DNA complex may be altered upon entry of Brf1 and Bdp1 into the complex. We show here, using the altered specificity TBP mutant TBPm3 and an iterative in vitro selection assay, that entry of Brf1 and Bdp1 into the complex imposes a strict sequence preference for the downstream half of the TATA box. Notably, the selected sequence (TGTAAATA) is a perfect match to the TATA box of the RNA polymerase III-transcribed U6 small nuclear RNA (SNR6) gene. We suggest that the selected T•A base pair step at the downstream end of the 8 bp TBP site may provide a DNA flexure that promotes TFIIIB-DNA complex formation.

Gcn5p contributes to the bidirectional character of the UGA3-GLT1 yeast promoter

Analysis of the UGA3-GLT1 bidirectional promoter has indicated that its transcriptional activation is determined by the combined action of Gcn4p and Gln3p, and that its bidirectional character is influenced by chromatin organization, through the action of an Abf1p binding site and a polydAdTtract. Results presented in this paper show that lack of Gcn5p impairs histone acetylation and nucleosomal organization of the UGA3-GLT1 promoter, resulting in an asymmetrical transcriptional activation response of UGA3 and GLT1. The phenotype displayed by a double mutant impaired in GCN5 and in the Abf1p binding site indicates that the combined action of these two elements determines the bidirectional capacity of the UGA3-GLT1 intergenic region.

3D reconstruction and comparison of shapes of DNA minicircles observed by cryo-electron microscopy

We use cryo-electron microscopy to compare 3D shapes of 158 bp long DNA minicircles that differ only in the sequence within an 18 bp block containing either a TATA box or a catabolite activator protein binding site. We present a sorting algorithm that correlates the reconstructed shapes and groups them into distinct categories. We conclude that the presence of the TATA box sequence, which is believed to be easily bent, does not significantly affect the observed shapes.

FOXM1c transactivates the human c-myc promoter directly via the two TATA boxes P1 and P2

FOXM1c transactivates the c-myc promoter via the P1 and P2 TATA boxes using a new mechanism. Whereas the P1 TATA box TATAATGC requires its sequence context to be FOXM1c responsive, the P2 TATA box TATAAAAG alone is sufficient to confer FOXM1c responsiveness to any minimal promoter. FOXM1c transactivates by binding to the TATA box as well as directly to TATA-binding protein, transcription factor IIB and transcription factor IIA. This new transactivation mechanism is clearly distinguished from the function of FOXM1c as a conventional transcription factor. The central domain of FOXM1c functions as an essential domain for activation via the TATA box, but as an inhibitory domain (retinoblastoma protein-independent transrepression domain and retinoblastoma protein-recruiting negative regulatory domain) for transactivation via conventional FOXM1c-binding sites. Each promoter with the P2 TATA box TATAAAAG is postulated to be transactivated by FOXM1c. This was demonstrated for the promoters of c-fos, hsp70 and histone H2B/a. A database search revealed almost 300 probable FOXM1c target genes, many of which function in proliferation and tumorigenesis. Accordingly, dominant-negative FOXM1c proteins reduced cell growth approximately threefold, demonstrating a proliferation-stimulating function for wild-type FOXM1c.

The UGA3-GLT1 intergenic region constitutes a promoter whose bidirectional nature is determined by chromatin organization in Saccharomyces cerevisiae

Transcription of an important number of divergent genes of Saccharomyces cerevisiae is controlled by intergenic regions, which constitute factual bidirectional promoters. However, few of such promoters have been characterized in detail. The analysis of the UGA3-GLT1 intergenic region has provided an interesting model to study the joint action of two global transcriptional activators that had been considered to act independently. Our results show that Gln3p and Gcn4p exert their effect upon cis-acting elements, which are shared in a bidirectional promoter. Accordingly, when yeast is grown on a low-quality nitrogen source, or under amino acid deprivation, the expression of both UGA3 and GLT1 is induced through the action of both these global transcriptional modulators that bind to a region of the bidirectional promoter. In addition, we demonstrate that chromatin organization plays a major role in the bidirectional properties of the UGA3-GLT1 promoter, through the action of an upstream Abf1p-binding consensus sequence and a polydAdT(tract). Mutations in these cis-elements differentially affect transcription of UGA3 and GLT1, and thus alter the overall relative expression. This is the first example of an intergenic region constituting a promoter whose bidirectional character is determined by chromatin organization.

Poor base stacking at DNA lesions may initiate recognition by many repair proteins

A fundamental question in DNA repair is how a mismatched or modified base is detected when embedded in millions to billions of normal base pairs. A survey of the literature and structural database reveals a common feature in all repair protein-DNA complexes: the DNA double helix is discontinuous at a lesion site due to base unstacking, kinking and/or nucleotide extrusion. Lesions induce destabilization and distortion of short linear DNAs, and underwinding in negatively supercoiled DNA presumably could compound the reduced stability caused by a lesion. A hypothesis is thus put forward that DNA lesion recognition occurs in two steps. Repair proteins initially recognize the weakened base stacking, and thus a flexible hinge at a DNA lesion. Sampling of flexible hinges rather than all DNA base pairs can reduce the task of finding a lesion by two to three orders of magnitude, from searching millions base pairs to thousands. After the initial encounter, a repair protein scrutinizes the shape, hydrogen bonding and electrostatic potentials of bases at the flexible hinge and dissociates if it is not a correct substrate. MutS, which has a broad range of substrates, actively dissociates from non-specific binding via an ATP-dependent proofreading mechanism. A single lesion may thus be sampled by BER, NER and MMR proteins until repaired. This proposition immediately suggests a mechanism for crosstalk between different repair and signaling pathways. It also raises the possibility that sampling of a lesion by one protein could facilitate loading of another by direct protein-protein or DNA mediated interactions.

The importance and identification of regulatory polymorphisms and their mechanisms of action

The search for the genetic variations underlying all human phenotypes is in its infancy but must be one of the long term goals of the scientific community. There is evidence that most, if not all human phenotypes, including illnesses are influenced by the genetic makeup of the individual. There are an estimated 11 million human genetic polymorphisms with a minor allele frequency >1% and possibly many times that number of rare sequence variants. The proportion of these sequence variants which have any functional effect is unknown but it is likely that the majority of those which influence illness lie outside of the amino acid coding regions of genes, and affect the regulation of gene expression--these are called rSNPs. Recent research suggests that about 50% of genes have one or more common rSNPs associated with them and probably most if not all genes have an rSNP within the human population. In the long term, determining which polymorphisms are potentially functional must be done bio-informatically using algorithms based upon experimental data. However, at the current time, the limited data that has been obtained does not allow the creation of such an algorithm. In vitro studies suggest that a large proportion of rSNPs lie within the core and proximal promoter regions of genes but it is not clear how the majority of these influence transcription, as they do not appear to be within any known transcription factor binding sites. However, promoter regions possess a number of sequence-dependent characteristics which make them distinct from the rest of the genome, namely stability, curvature and flexibility. Subtle changes to these features may underlie the mechanisms by which many polymorphisms exert their function.

DNA repair in a protein-DNA complex: searching for the key to get in

An obstacle encountered by nucleotide excision repair (NER) proteins during repair of the genome is the masking of bulky lesions by DNA binding proteins. For example, certain transcription factors are known to be impediments, and suppress damage removal at their recognition sequences. We have used well-defined protein-DNA complexes to study the molecular mechanism(s) used by repair proteins in gaining access to DNA lesions in chromatin. Using transcription factor IIIA (TFIIIA) and the 5S ribosomal RNA gene (5S rDNA), we previously measured position-dependent effects of cyclobutane pyrimidine dimers (CPDs) at five different sites within the internal control region (ICR) on complex formation [Y. Kwon, M.J. Smerdon, Binding of zinc finger protein transcription factor IIIA to its cognate DNA sequence with single UV photoproducts at specific sites and its effect on DNA repair, J. Biol. Chem. 278 (2003) 45451-45459]. We found that CPDs at two of these sites enhance the TFIIIA-rDNA dissociation rate, which correlates with enhanced repair at these two sites. Here, we used a novel approach to directly compare dissociation of randomly damaged rDNA with NER. We refined the relationship between dissociation and repair of the complex by examining all CPD sites in the transcribed strand. A 214 bp 5S rDNA fragment was irradiated with UV light to produce CPDs at dipyrimidine sites and approximately 1 CPD per fragment. Positions of CPDs that alter binding of TFIIIA were determined by T4 endonuclease V mapping of TFIIIA-bound and unbound fractions of UV-irradiated DNA. As expected, the results reveal that dissociation of TFIIIA from the complex is significantly enhanced by CPDs within the ICR. Moreover, the levels of dissociation induced by CPDs were quantitatively compared with their repair efficiency, and indicate that repair rates of most CPDs in the complex closely correlate with the dissociation rates. In addition, changes in dissociation rate are similar to changes in CPD frequency induced by TFIIIA binding. These findings indicate that structural compatibility of a DNA lesion within a protein-DNA complex can determine both lesion frequency and repair efficiency.

DNA bending by M.EcoKI methyltransferase is coupled to nucleotide flipping

The maintenance methyltransferase M.EcoKI recognizes the bipartite DNA sequence 5′-AACNNNNNNGTGC-3′, where N is any nucleotide. M.EcoKI preferentially methylates a sequence already containing a methylated adenine at or complementary to the underlined bases in the sequence. We find that the introduction of a single-stranded gap in the middle of the non-specific spacer, of up to 4 nt in length, does not reduce the binding affinity of M.EcoKI despite the removal of non-sequence-specific contacts between the protein and the DNA phosphate backbone. Surprisingly, binding affinity is enhanced in a manner predicted by simple polymer models of DNA flexibility. However, the activity of the enzyme declines to zero once the single-stranded region reaches 4 nt in length. This indicates that the recognition of methylation of the DNA is communicated between the two methylation targets not only through the protein structure but also through the DNA structure. Furthermore, methylation recognition requires base flipping in which the bases targeted for methylation are swung out of the DNA helix into the enzyme. By using 2-aminopurine fluorescence as the base flipping probe we find that, although flipping occurs for the intact duplex, no flipping is observed upon introduction of a gap. Our data and polymer model indicate that M.EcoKI bends the non-specific spacer and that the energy stored in a double-stranded bend is utilized to force or flip out the bases. This energy is not stored in gapped duplexes. In this way, M.EcoKI can determine the methylation status of two adenine bases separated by a considerable distance in double-stranded DNA and select the required enzymatic response.

Cytosine arabinoside substitution decreases transcription factor-DNA binding element complex formation

CONTEXT

The pyrimidine nucleoside analog, cytosine arabinoside (Ara-C), is an effective therapeutic agent for acute leukemia. The phosphorylated triphosphate, cytosine arabinoside triphosphate, competes with deoxycytosine triphosphate as a substrate for incorporation into DNA. Once incorporated into DNA, it inhibits DNA polymerase and topoisomerase I and modifies the tertiary structure of DNA.

OBJECTIVE

To determine if the substitution of Ara-C for cytosine in double-stranded oligonucleotides that contain 4 specific transcription factor binding sites (TATA, GATA, C/EBP, and AP-2alpha) alters transcription factor binding to their respective DNA binding elements.

DESIGN

Transcription factors were obtained from nuclear extracts from human promyelocytic leukemia HL-60 cells. [32P]-end-labeled double-stranded oligonucleotides that contained 1 or 2 specific transcription factor binding sites with or without Ara-C substitution for cytosine were used to assess transcription factor binding by electrophoretic mobility shift assay.

RESULTS

The substitution of Ara-C for cytosine within and outside the transcription factor binding element (AP-2alpha, C/EBP), outside the binding element only (GATA, TATA), or within the binding element only (AP-2alpha) all result in a reduction in transcription factor binding to their respective DNA binding element.

CONCLUSION

The reduction of the binding capacity of transcription factors with their respective DNA binding elements may depend on structural changes within oligonucleotides induced by Ara-C incorporation. This altered binding capacity of transcription factors to their DNA binding elements may represent one mechanism for Ara-C cytotoxicity secondary to inhibition of transcription of new messenger RNAs and, subsequently, translation of new proteins.

BMP4-dependent expression of Xenopus Grainyhead-like 1 is essential for epidermal differentiation

Morphogen-dependent epidermal-specific transacting factors have not been defined in vertebrates. We demonstrate that a member of the grainyhead transcription factor family, Grainyhead-like 1 (XGrhl1) is essential for ectodermal ontogeny in Xenopus laevis. Expression of this factor is restricted to epidermal cells. Moreover, XGrhl1 is regulated by the BMP4 signaling cascade. Disruption of XGrhl1 activity in vivo results in a severe defect in terminal epidermal differentiation, with inhibition of XK81A1 epidermal keratin gene expression, a key target of BMP4 signaling. Furthermore, transcription of the XK81A1 gene is modulated directly by binding of XGRHL1 to a promoter-localized binding motif that is essential for high-level expression. These results establish a novel developmental role for XGrhl1 as a crucial tissue-specific regulator of vertebrate epidermal differentiation.

The contribution of phosphate-phosphate repulsions to the free energy of DNA bending

DNA bending is important for the packaging of genetic material, regulation of gene expression and interaction of nucleic acids with proteins. Consequently, it is of considerable interest to quantify the energetic factors that must be overcome to induce bending of DNA, such as base stacking and phosphate–phosphate repulsions. In the present work, the electrostatic contribution of phosphate–phosphate repulsions to the free energy of bending DNA is examined for 71 bp linear and bent-form model structures. The bent DNA model was based on the crystallographic structure of a full turn of DNA in a nucleosome core particle. A Green's function approach based on a linear-scaling smooth conductor-like screening model was applied to ascertain the contribution of individual phosphate–phosphate repulsions and overall electrostatic stabilization in aqueous solution. The effect of charge neutralization by site-bound ions was considered using Monte Carlo simulation to characterize the distribution of ion occupations and contribution of phosphate repulsions to the free energy of bending as a function of counterion load. The calculations predict that the phosphate–phosphate repulsions account for ∼30% of the total free energy required to bend DNA from canonical linear B-form into the conformation found in the nucleosome core particle.

Transcriptional frequency and cell determination

The relative base composition of DNA regulatory sequences of certain genes of undetermined multipotent progenitor cells may account for the frequency of transcription of these genes in cell determination. The sequences of these regulatory regions of cell determination genes that are more AT-rich would create the potential for transcription at a higher frequency due to their lower melting temperature, as well as propensity to bend. An increase of one or more of the high mobility group (HMG) chromatin proteins would preferentially bind the more AT-rich regulatory sequences, thereby increasing the rate of transcription. The amount of unphosphorylated H1 histone reacting with these same regulatory sites may decrease transcription frequency. The level of cell growth, i.e. total protein synthesis of a cell, is correlated positively with the synthesis of HMG proteins. H1 histone synthesis is linked to DNA replication. Unbalanced growth would alter the amounts of HMG proteins and H1 histone, thus changing transcriptional frequency. The greater the enrichment of AT sequences in the regulatory regions of the cell determination genes, the greater may be the extent of evolutionary conservation. Higher frequency of transcription of the cell determination genes with the more AT-rich regulatory sequences could account for the earlier expression of the more conserved cell determination genes during embryonic development. Preferential binding of H1 histone to the more AT-rich regulatory sequences would subsequently restrict their transcription before that of less conserved cell determination genes.

Looking into DNA recognition: zinc finger binding specificity

We present a quantitative, theoretical analysis of the recognition mechanisms used by two zinc finger proteins: Zif268, which selectively binds to GC-rich sequences, and a Zif268 mutant, which binds to a TATA box site. This analysis is based on a recently developed method (ADAPT), which allows binding specificity to be analyzed via the calculation of complexation energies for all possible DNA target sequences. The results obtained with the zinc finger proteins show that, although both mainly select their targets using direct, pairwise protein-DNA interactions, they also use sequence-dependent DNA deformation to enhance their selectivity. A new extension of our methodology enables us to determine the quantitative contribution of these two components and also to measure the contributions of individual residues to overall specificity. The results show that indirect recognition is particularly important in the case of the TATA box binding mutant, accounting for 30% of the total selectivity. The residue-by-residue analysis of the protein-DNA interaction energy indicates that the existence of amino acid-base contacts does not necessarily imply sequence selectivity, and that side chains without contacts can nevertheless contribute to defining the protein's target sequence.