Does TATA matter? A structural exploration of the selectivity determinants in its complexes with TATA box-binding protein

Tuning the Transcriptional Activity of the CaMV 35S Promoter in Plants by Single-Nucleotide Changes in the TATA Box

Synthetic biology uses genetically encoded devices and circuits to implement novel complex functions in living cells and organisms. A hallmark of these genetic circuits is the interaction among their individual parts, according to predefined rules, to process cellular information and produce a circuit output or response. As the number of individual components in a genetic circuit increases, so does the number of interactions needed to achieve the correct behavior, and hence, a greater need to fine-tune the levels of expression of each component. Transcriptional promoters play a key regulatory role in genetic circuits, as they influence the levels of RNA and proteins produced. In multicellular organisms, such as plants, they can also determine developmental, spatial, and tissue-specific patterns of gene expression. The 35S promoter from the Cauliflower Mosaic Virus (CaMV 35S) is widely used in plant biotechnology to direct high levels of gene expression in a variety of plant species. We produced a library of 21 variants of the CaMV 35S promoter by introducing all single nucleotide substitutions to the promoter's TATA box sequence. We then characterized the activity of all variants in homozygous transgenic plants and showed that some of these variants have lower activity than the wild type in plants. These promoter variants could be used to fine-tune the behavior of synthetic genetic circuits in plants.

Recent developments in empirical atomistic force fields for nucleic acids and applications to studies of folding and dynamics

Nucleic acids play critical roles in carrying genetic information, participating in catalysis, and preserving chromosomal structure. Despite over a century of study, efforts to understand the dynamics and structure-function relationships of DNA and RNA at the atomic level are still ongoing. Molecular dynamics (MD) simulations augment experiments by providing atomistic resolution and quantitative relationships between structure and conformational energy. Steady advancements in computer hardware, software, and atomistic force fields (FFs) over 40 years have facilitated new discoveries. Here, we review nucleic acid FF development with emphasis on recent refinements that have improved descriptions of important nucleic acid properties. We then discuss several key examples of successes and challenges in modeling nucleic acid structure and dynamics using the latest FFs.

Effect of temperature on the structure and hydration layer of TATA-box DNA: A molecular dynamics simulation study

DNA within the living cells experiences a diverse range of temperature, ranging from freezing condition to hot spring water. How the structure, the mechanical properties of DNA, and the solvation dynamics around DNA changes with the temperature is important to understand the functionality of DNA under those acute temperature conditions. In that notion, we have carried out molecular dynamics simulations of a DNA oligomer, containing TATA-box sequence for three different temperatures (250K, 300K and 350K). We observed that the structure of the DNA, in terms of backbone torsion angles, sugar pucker, base pair parameters, and base pair step parameters, did not show any unusual properties within the studied range of temperatures, but significant structural alteration was noticed between BI and BII forms at higher temperature. As expected, the flexibility of the DNA, in terms of the torsional rigidity and the bending rigidity is highly temperature dependent, confirming that flexibility increases with increase in temperature. Additionally, the groove widths of the studied DNA showed temperature sensitivity, specifically, the major groove width decreases and the minor groove width increases, respectively, with the increase in temperature. We observed that at higher temperature, water around both the major and the minor groove of the DNA is less structured. However, the water dynamics around the minor groove of the DNA is more restricted as compared to the water around the major groove throughout the studied range of temperatures, without any anomalous behavior.

Insights into protein-DNA interactions through structure network analysis

Protein–DNA interactions are crucial for many cellular processes. Now with the increased availability of structures of protein–DNA complexes, gaining deeper insights into the nature of protein–DNA interactions has become possible. Earlier, investigations have characterized the interface properties by considering pairwise interactions. However, the information communicated along the interfaces is rarely a pairwise phenomenon, and we feel that a global picture can be obtained by considering a protein–DNA complex as a network of noncovalently interacting systems. Furthermore, most of the earlier investigations have been carried out from the protein point of view (protein-centric), and the present network approach aims to combine both the protein-centric and the DNA-centric points of view. Part of the study involves the development of methodology to investigate protein–DNA graphs/networks with the development of key parameters. A network representation provides a holistic view of the interacting surface and has been reported here for the first time. The second part of the study involves the analyses of these graphs in terms of clusters of interacting residues and the identification of highly connected residues (hubs) along the protein–DNA interface. A predominance of deoxyribose–amino acid clusters in β-sheet proteins, distinction of the interface clusters in helix–turn–helix, and the zipper-type proteins would not have been possible by conventional pairwise interaction analysis. Additionally, we propose a potential classification scheme for a set of protein–DNA complexes on the basis of the protein–DNA interface clusters. This provides a general idea of how the proteins interact with the different components of DNA in different complexes. Thus, we believe that the present graph-based method provides a deeper insight into the analysis of the protein–DNA recognition mechanisms by throwing more light on the nature and the specificity of these interactions.

The interaction of proteins with DNA is crucial for several cellular processes. Some insights into the mode of interaction can be obtained from the analysis of the complexed structures. Conventional analyses are based on the identification of pairwise interactions. However, a collective representation of the network of interactions and the analyses of such networks provide valuable information, which is not easy to obtain from pairwise analyses. Although the protein structure networks have been described in the literature, this is the first time that a network representation of protein–DNA is described. Construction and analysis of such networks have given valuable information on protein–DNA interactions in terms of network parameters, such as clusters of interacting residues and hubs, which are highly connected residues. Furthermore, the results also represent both the protein- and the DNA-centric viewpoints, because the analysis is carried out on combined networks. The methodology developed here can lead to predictions, such as important residues responsible for stabilizing protein–DNA interactions, and will be of interest to experimentalists.

Sequence-dependence of the energetics of opening of at basepairs in DNA

Proton exchange and nuclear magnetic resonance spectroscopy are being used to characterize the energetics of opening of AT/TA basepairs in the DNA dodecamer 5'-d(GCTATAAAAGGG)-3'/5'-d(CCCTTTTATAGC)-3'. The dodecamer contains the TATA box of the adenovirus major late promoter. The equilibrium constants for opening of each basepair are measured from the dependence of the exchange rates of imino protons on ammonia concentration. The enthalpy, entropy, and free energy changes in the opening reaction of each basepair are determined from the temperature dependence of the exchange rates. The results reveal that the opening enthalpy changes encompass a wide range of values, namely, from 17 to 29 kcal/mol. The largest values are observed for the AT basepairs in 7th and 8th positions. These values, and the exchange rates of the corresponding imino protons, suggest that these two basepairs open in a single concerted reaction. The enthalpy changes for opening of the central six basepairs are correlated to the opening entropy changes. This enthalpy-entropy compensation minimizes the variations in the opening free energies among these central basepairs. Deviations from the enthalpy-entropy compensation pattern are observed for basepairs located close to the ends of the duplex structure, suggesting a different mode of opening for these basepairs.

DNA dynamically directs its own transcription initiation

It has long been known that double-stranded DNA is subject to temporary, localized openings of its two strands. Particular regions along a DNA polymer are destabilized structurally by available thermal energy in the system. The localized sequence of DNA determines the physical properties of a stretch of DNA, and that in turn determines the opening profile of that DNA fragment. We show that the Peyrard-Bishop nonlinear dynamical model of DNA, which has been used to simulate denaturation of short DNA fragments, gives an accurate representation of the instability profile of a defined sequence of DNA, as verified using S1 nuclease cleavage assays. By comparing results for a non-promoter DNA fragment, the adenovirus major late promoter, the adeno-associated viral P5 promoter and a known P5 mutant promoter that is inactive for transcription, we show that the predicted openings correlate almost exactly with the promoter transcriptional start sites and major regulatory sites. Physicists have speculated that localized melting of DNA might play a role in gene transcription and other processes. Our data link sequence-dependent opening behavior in DNA to transcriptional activity for the first time.

New Analogues of Amonafide and Elinafide, Containing Aromatic Heterocycles: Synthesis, Antitumor Activity, Molecular Modeling, and DNA Binding Properties

Amonafide- and elinafide-related mono and bisintercalators, modified by the introduction of a pi-excedent furan or thiophene ring fused to the naphthalimide moiety, have been synthesized. These compounds have shown an interesting antitumor profile. The best compound, 9, was 2.5-fold more potent than elinafide against human colon carcinoma cells (HT-29). Molecular dynamic simulations and physicochemical experiments have demonstrated that these compounds are capable of forming stable DNA complexes. These results, together with those previously reported by us for imidazo- and pyrazinonaphthalimide analogues, have prompted us to propose that the DNA binding process does not depend on the electronic nature of the fused heterocycle.

Sequence-dependent solution structure and motions of 13 TATA/TBP (TATA-box binding protein) complexes

The TATA element is a well-known example of a DNA promoter sequence recognized by the TATA box binding protein (TBP) through its intrinsic motion and deformability. Although TBP recognizes the TATA element octamer unusually (through the minor groove, which lacks the distinctive features of the major groove), single base-pair replacements alter transcriptional activity. Recent crystallographic experiments have suggested that TATA/TBP complexes differing by a single base pair retain substantial structural similarity despite their functional differences in activating transcription. To investigate the subtle role of sequence-dependent motion within the TATA element and certain aspects of its effect on assembly of the transcriptional complex, we examine 5-ns dynamics trajectories of 13 variant TATA/TBP complexes differing from each other by a single base pair. They include the wild-type (WT) adenovirus 2 major late promoter (AdMLP) TATA element, TATAAAAG (the octamer specifies positions -31 to -24 with respect to the transcription initiation site), and the variants A31 (i.e., AATAAAAG), T30, A29, C29, G28, T28, T27, G26, T26, C25, T25, and T24. Our simulated TATA/TBP complexes develop sequence-dependent structure and motion trends that may lead to favorable orientations for high-activity variants (with respect to binding TFIIA, TFIIB, and other transcription factors), while conversely, accelerate dissociation of low-activity TATA/TBP complexes. The motions that promote favorable geometries for preinitiation complexes include small rotations between TBP's N- and C-terminal domains, sense strand DNA backbone "slithering," and rotations in TBP's H2 and H2' helices. Low-activity variants tend to translate the H1 and H1' helices and withdraw the intercalating phenylalanines. These cumulative DNA and protein motions lead to a spatial spread of complex orientations up to 4 A; this is associated with an overall bend of the variant TATA/TBP complexes that spans 93 degrees to 110 degrees (107 degrees for the crystal reference). Taken together, our analyses imply larger differences when these local structural and bending changes are extended to longer DNA (upstream and downstream) and suggest that specific local TATA/TBP motions (e.g., shifts in TBP helices and TATA bases and backbone) play a role in modulating the formation and maintenance of the transcription initiation complex.

Comparison of TATA-binding protein recognition of a variant and consensus DNA promoters

Assembly of transcription pre-initiation complexes proceeds from the initial complex formed between "TATA" bearing promoter DNA and the TATA-binding protein (TBP). Our laboratory has been investigating the relationships among TATA sequence, TBP center dot TATA solution structure, recognition mechanisms, and transcription efficiency. TBP center dot TATA interactions have been modeled by global analysis of detailed kinetic and thermodynamic data obtained using fluorimetric and fluorometric techniques in conjunction with fluorescence resonance energy transfer. We have reported recently that TBP recognition of two consensus promoters, adenovirus major late (AdMLP: TATAAAAG) and E4 (TATATATA), is well described by a linear two-intermediate mechanism with simultaneous DNA binding and bending. Similar DNA geometries and high transcription efficiencies characterize these TBP x TATA complexes. Here we show that, in contrast to the consensus sequences, TBP recognition of a variant sequence (C7: TATAAACG) is described by a three-step model with two branching pathways. One pathway proceeds through an intermediate having severely bent DNA, reminiscent of the consensus interactions, with the other branch yielding a unique conformer with shallowly bent DNA. The resulting TBP x C7 complex has a dramatically different solution conformation than for TBP x DNA(CONSENSUS) and is correlated with diminished relative transcription activity. The temperature dependence of the TBP x C7 helical bend is postulated to derive from population shifts between the conformers with slightly and severely bent DNA.

ADAPT: a molecular mechanics approach for studying the structural properties of long DNA sequences

We describe an original approach to determining sequence-structure relationships for DNA. This approach, termed ADAPT, combines all-atom molecular mechanics with a multicopy algorithm to build nucleotides that contain all four standard bases in variable proportions. These nucleotides enable us to search very rapidly for base sequences that energetically favor chosen types of DNA deformation or chosen DNA-protein or DNA-ligand interactions. Sequences satisfying the chosen criteria can be found by energy minimization, combinatorial sequence searching, or genome scanning, in a manner similar to the threading approaches developed for protein structure prediction. In the latter case, we are able to analyze roughly 2000 base pairs per second. Applications of the method to DNA allomorphic transitions, DNA deformation, and specific DNA interactions are presented.

Molecular dynamics simulations of B '-DNA: sequence effects on A-tract-induced bending and flexibility

Molecular dynamics (MD) simulations including water and counterions are reported on five examples of A-tract DNA oligonucleotide dodecamer duplexes for which crystal structures are available, the homopolymeric duplex sequences poly(dA) and poly(dG), and two related sequences that serve as controls. MD was performed using the AMBER suite of programs for 3 ns on each sequence. These results, combined with previously reported MDs on 25-mer and 30-mer oligonucleotides on sequences with phased A-tracts carried out under a similar simulation protocol, are used to examine salient issues in the structural chemistry of ApA steps and A-tract induced axis bending. MD modeling successfully describes the distinctive B' structure of A-tracts in solution as essentially straight (wedge angles of <1 degrees ), more rigid than generic B-form DNA, with slight base-pair inclination, high propeller twist and a minor groove narrowing 5' to 3'. The MD structures in solution agree closely with corresponding crystal structures, supporting the idea that crystal structures provide a good model for A-tract DNA structure in solution. From the collective MD results, bending and flexibility are calculated by step. Pyrimidine-purine steps are predicted to be most intrinsically bent and also most bendable, i.e. susceptible to bending. Pyrimidine-pyrimidine ( approximately purine-purine) and purine-pyrimidine steps show less intrinsic deformation and deformability. The MD calculated flexibility correlates well with the protein-induced bendability derived independently from the protein DNA crystal structures. The MD results indicate that bending and flexibility of base-pair steps in DNA are highly correlated, i.e. steps that exhibit the most intrinsic deformation from B-form DNA turn are also the most dynamically deformable. The MD description of A-tract-induced axis bending shows most consistency with the non A-tract, general-sequence model, in which the sequence curvature originates primarily in base-pair roll towards the major groove in non-A-tract regions of the sequence, particularly pyrimidine-purine steps. The direction of curvature is towards the minor groove viewed from opposite the A-tracts, but the A-tracts per se exhibit only minor deformation. The MD results are found to be consistent with the directionality of bending inferred for DNA sequences from gel retardation and cyclization experiments.

Dynamic simulations of 13 TATA variants refine kinetic hypotheses of sequence/activity relationships 1 1Edited by B. Honig

The fundamental relationship between DNA sequence/deformability and biological function has attracted numerous experimental and theoretical studies. A classic prototype system used for such studies in eukaryotes is the complex between the TATA element transcriptional regulator and the TATA-box binding protein (TBP). The recent crystallographic study by Burley and co-workers demonstrated the remarkable structural similarity contrasted to different transcriptional activity of 11 TBP/DNA complexes in which the DNAs differed by single base-pairs. By simulating these TATA variants and two other single base-pair variants that were not crystallizable, we uncover sequence-dependent structural, energetic, and flexibility properties that tailor TATA elements to TBP interactions, complementing many previous studies by refining kinetic hypotheses on sequence/activity correlations. The factors that combine to produce favorable elements for TBP activity include overall flexibility; minor groove widening, as well as roll, rise, and shift increases at the ends of the TATA element; untwisting within the TATA element accompanied by large roll at the TATA element ends; and relatively low maximal water densities around the DNA. These features accompany the severe deformation induced by the minor-groove binding protein, which kinks the TATA element at the ends and displaces local water molecules to form stabilizing hydrophobic contacts. Interestingly, the preferred bending direction itself is not a significant predictor of activity disposition, although certain variants (such as wild-type AdMLP, 5'-TATA4G-3', and inactive A29, 5'-TA6G-3') exhibit large preferred bends in directions consistent with their activity or inactivity (major groove and minor groove bends, respectively). These structural, flexibility, and hydration preferences, identified here and connected to a new crystallographic study of a larger group of DNA variants than reported to date, highlight the profound influence of single base-pair DNA variations on DNA motion. Our refined kinetic hypothesis suggests the functional implications of these motions in a kinetic model of TATA/TBP recognition, inviting further theoretical and experimental research.

Molecular dynamics simulation of nucleic acids

We review molecular dynamics simulations of nucleic acids, including those completed from 1995 to 2000, with a focus on the applications and results rather than the methods. After the introduction, which discusses recent advances in the simulation of nucleic acids in solution, we describe force fields for nucleic acids and then provide a detailed summary of the published literature. We emphasize simulations of small nucleic acids ( approximately 6 to 24 mer) in explicit solvent with counterions, using reliable force fields and modern simulation protocols that properly represent the long-range electrostatic interactions. We also provide some limited discussion of simulation in the absence of explicit solvent. Absent from this discussion are results from simulations of protein-nucleic acid complexes and modified DNA analogs. Highlights from the molecular dynamics simulation are the spontaneous observation of A B transitions in duplex DNA in response to the environment, specific ion binding and hydration, and reliable representation of protein-nucleic acid interactions. We close by examining major issues and the future promise for these methods.

Site-specific cation binding mediates TATA binding protein-DNA interaction from a hyperthermophilic archaeon

Pyrococcus woesei (Pw) is a hyperthermophilic archaeal organism that exists under conditions of high salt and elevated temperature. In a previous study [O'Brien, R., DeDecker, B., Fleming, K., Sigler, P. B., and Ladbury, J. E., (1998) J. Mol. Biol. 279, 117-125], we showed that, despite the similarity of primary and secondary structure, the TATA box binding protein (TBP) from Pw binds thermodynamically in a fundamentally different way to its mesophilic counterparts. The affinity of the interaction increases as the salt concentration is increased. The formation of the protein-DNA complex involves the release of water and the uptake of ions, which were hypothesized to be cations. Here we test this hypothesis by selecting potential cation binding sites at negatively charged, acidic residues in the complex interface. These were substituted using site-directed mutagenesis of specific residues. Changes in the thermodynamic parameters on formation of the mutant protein-DNA complex were determined using isothermal titration calorimetry and compared to the wild type interaction. Removal of a glutamate residue from the binding site resulted in the uptake of one less cation on formation of the complex. This glutamate (E12) is directly involved in the binding of cations in the complex interface. Substitution of another acidic residue proximal to the DNA binding site (D101) had no effect on cation uptake, suggesting that the location of the amino acid on the protein surface is important in dictating the potential to coordinate cations. Removal of the cation binding site provided a more favorable entropy of binding; however, this effect is significantly reduced at higher salt concentrations. The removal of the cation binding site led to an increase in affinity with respect to the wild-type TBP at low salt concentrations.

A detailed interpretation of OH radical footprints in a TBP-DNA complex reveals the role of dynamics in the mechanism of sequence-specific binding

The hydroxyl radical footprint of the TATA-binding protein (TBP) bound to the high-affinity sequence TATAAAAG of the adenovirus 2 major late promoter has been quantitatively compared to a 2 ns molecular dynamics simulation of the complex in aqueous solution at room temperature using the CHARMM23 potential. The nucleotide-by-nucleotide analysis of the TBP-TATA hydroxyl radical footprint correlates with the solvent-accessible surface calculated from the dynamics simulation. The results suggest that local reactivity towards OH radicals results from the interplay between the local DNA geometry imposed by TBP binding, and the dynamics of the side-chains contacting the sugar hydrogen atoms. Analysis of the dynamics suggests that, over time, TBP forms stable interactions with the sugar-phosphate backbone through multiple contacts to different partners. This mechanism results in an enthalpic advantage to complex formation at a low entropic cost.

On the truncation of long-range electrostatic interactions in DNA

Long-range interactions are known to play an important role in highly polar biomolecules like DNA. In molecular dynamics simulations of nucleic acids and proteins, an accurate treatment of the long-range interactions are crucial for achieving stable nanosecond trajectories. In this report, we evaluate the structural and dynamic effects on a highly charged oligonucleotide in aqueous solution from different long-range truncation methods. Two group-based truncation methods, one with a switching function and one with a force-switching function were found to fail to give accurate stable trajectories close to the crystal structure. For these group-based truncation methods, large root mean square (rms) deviations from the initial structure were obtained and severe distortions of the oligonucleotide were observed. Another group-based truncation scheme, which used an abrupt truncation at 8. 0 A or at 12.0 A was also investigated. For the short cutoff distance, the conformations deviated far away from the initial structure and were significantly distorted. However, for the longer cutoff, where all necessary electrostatic interactions were included, the trajectory was quite stable. For the particle mesh Ewald (PME) truncation method, a stable DNA simulation with a heavy atom rms deviation of 1.5 A was obtained. The atom-based truncation methods also resulted in stable trajectories, according to the rms deviation from the initial B-DNA structure, of between 1.5 and 1.7 A for the heavy atoms. In these stable simulations, the heavy atom rms deviations were approximately 0.6-1.0 A lower for the bases than for the backbone. An increase of the cutoff radius from 8 to 12 A decreased the rms deviation by approximately 0.2 A for the atom-based truncation method with a force-shifting function, but increased the computational time by a factor of 2. Increasing the cutoff from 12 to 18 A for the atom-based truncation method with a force-shifting function requires 2-3 times more computational time, but did not significantly change the rms deviation. Similar rms deviations from the initial structure were found for the atom-based method with a force-shifting function and for the PME method. The computational cost was longer for the PME method with a cutoff of 12. 0 A for the direct space nonbonded calculations than for the atom-based truncation method with a force-shifting function and a cutoff of 12.0 A. If a nonperiodic boundary, e.g., a spherical boundary, was used, a considerable speedup could be achieved. From the rms fluctuations, the terminal nucleotides and especially the cytidines were found to be more flexible than the nonterminal nucleotides. The B-DNA form of the oligonucleotide was maintained throughout the simulations and is judged to depend on the parameters of the energy function and not on the truncation method used to handle the long-range electrostatic interactions. To perform accurate and stable simulations of highly charged biological macromolecules, we recommend that the atom-based force-shift method or the PME method should be used for the long-range electrostatics interactions.

The exchange of cognate TATA boxes results in a corresponding change in the strength of two HSV-1 early promoters

Previous analysis of two Herpes Simplex Virus Type-1 (HSV-1) promoters controlling expression of mRNA encoding early genes (U(L)37 and U(L)50) showed that the U(L)50 (dUTPase) promoter is at least 6-fold stronger both in its normal genomic location and in the non-essential gC locus. In the present report we demonstrate that the TATA element of either promoter is the major determinant of promoter strength. When the U(L)37 TATA element (CGTATAAC) was mutated with two base changes to the U(L)50 TATA sequence (CATAAAAC) in recombinant viruses, the activity of the U(L)37 promoter was increased to that of the U(L)50 promoter. Conversely, when the U(L)50 TATA element was changed to that of the U(L)37 promoter, U(L)50 promoter activity was reduced to the level of the U(L)37 promoter. In addition, we investigated the spacing of the TATA box with respect to upstream promoter elements. We found that re-positioning the U(L)37 TATA box to a location equivalent to that of the U(L)50 promoter relative to the transcript start site; i.e. three bases upstream of its cognate location, significantly diminished activity. Substitution of the U(L)50 TATA box at the new position could only partially restore promoter activity. Thus, we also conclude that the spacing of TATA elements vis-à-vis upstream promoter elements is also a critical determinant of promoter strength.

Nucleic acids: theory and computer simulation, Y2K

Molecular dynamics simulations on DNA and RNA that include solvent are now being performed under realistic environmental conditions of water activity and salt. Improvements to force-fields and treatments of long-range interactions have significantly increased the reliability of simulations. New studies of sequence effects, axis bending, solvation and conformational transitions have appeared.

Binding mechanisms of TATA box-binding proteins: DNA kinking is stabilized by specific hydrogen bonds

One of the common mechanisms of DNA bending by minor groove-binding proteins is the insertion of protein side chains between basepair steps, exemplified in TBP (TATA box-binding protein)/DNA complexes. At the central basepair step of the TATA box TBP produces a noticeable decrease in twist and an increase in roll, while engaging in hydrogen bonds with the bases and sugars. This suggests a mechanism for the stabilization of DNA kinks that was explored here with ab initio quantum mechanical calculations and molecular dynamics/potential of mean force calculations. The hydrogen bonds are found to contribute the energy necessary to drive the conformational transition at the central basepair step. The Asn, Thr, and Gly residues involved in hydrogen bonding to the DNA bases and sugar oxygens form a relatively rigid motif in TBP. The interaction of this motif with DNA is found to be responsible for inducing the untwisting and rolling of the central basepair step. Notably, direct readout is shown not to be capable of discriminating between AA and AT steps, as the strength of the hydrogen bonds between TBP and the DNA are the same for both sequences. Rather, the calculated free energy cost for an equivalent conformational transition is found to be sequence-dependent, and is calculated to be higher for AA steps than for AT steps.

How NF-kappaB can be attracted by its cognate DNA

NF-kappaB is involved in the transcriptional regulation of a large number of genes, in particular those of human immunodeficiency virus (HIV). Recently, we used NMR spectroscopy and molecular modelling to study the solution structure of a native duplex related to the HIV-1 kappaB site, together with a mutated duplex for which a three base-pair change abolishes NF-kappaB binding. The native duplex shows unusual dynamics of the four steps surrounding the kappaB site. Here, we explore the intrinsic properties of the NMR-refined structures of both duplexes in order to understand why the native sequence is recognised by NF-kappaB among other DNA sequences. We establish that only the native kappaB site can adopt a conformation where its structure (curvature and base displacement), the accessibility and the electrostatic potentials of key atoms become very favourable for binding the large loops of NF-kappaB, in contrast to the mutated duplex. Finally, we show that the neutralisation of phosphate groups contacted by NF-kappaB favours a more canonical DNA structure. These findings lead to a new hypothesis for specific recognition through the phosphodiester backbone dynamics of the sequences flanking a binding site. Such unusual behaviour confers upon the overall duplex properties that can be used by NF-kappaB to select its binding site. Thus, the selectivity determinants for NF-kappaB binding appear to depend on deformability of an "extended" consensus sequence.

TATA box DNA deformation with and without the TATA box-binding protein

DNA ring closure methods have been applied to TATA box DNA and its complex with the TATA box-binding protein (TBP). The J factors for cyclization (effective concentrations of one DNA end about the other) have been measured using cyclization kinetics, with and without bound TBP, for 18 DNA constructs containing the adenovirus major late promoter TATA box (TATAAAAG) separated by a variable helical phasing adapter from sequence-induced A-tract DNA bends. Six phasing lengths were used at three overall DNA lengths each. Cyclization kinetics were also measured in the absence of protein for the same set of molecules bearing a mutant TATA box (TACAAAAG). The results suggest that the TATA box DNA itself is strongly bent and anisotropically flexible, in a direction opposite to the bend induced by TBP, and that the mutant TACA box is much less bent/flexible. The bending and flexibility of the free DNA may govern the energetics of recognition of different DNA sequences by TBP, and the intrinsic bend may act to repress transcription complex assembly in the absence of TBP. The cyclization kinetics of TBP-DNA complexes in solution predict a geometry generally consistent with crystal structures, which show dramatic bending and unwinding. The novel observation of TBP-induced topoisomers suggests that this minicircle approach is able to distinguish TBP-induced unwinding from writhe (these cancel out in larger DNA), and this in turn suggests that changes in supercoiling in small topological domains can control TBP binding.

Modeling DNA deformations induced by minor groove binding proteins

Molecular modeling is used to demonstrate that the major structural deformations of DNA caused by four different minor groove binding proteins, TBP, SRY, LEF-1, and PurR, can all be mimicked by stretching the double helix between two 3'-phosphate groups flanking the binding region. This deformation reproduces the widening of the minor groove and the overall bending and unwinding of DNA caused by protein binding. It also reproduces the principal kinks associated with partially intercalated amino acid side chains, observed with such interactions. In addition, when protein binding involves a local transition to an A-like conformation, phosphate neutralization, via the formation of protein-DNA salt bridges, appears to favor the resulting deformation.

TIT for TAT: the properties of inosine and adenosine in TATA box DNA

The sequence dependent conformation, flexibility and hydration properties of DNA molecules constitute selectivity determinants in the formation of protein-DNA complexes. TATA boxes in which AT basepairs (bp) have been substituted by IC bp (TITI box) allow for probing these selectivity determinants for the complexation with the TATA box-binding protein (TBP) with different sequences but identical chemical surfaces. The reference promoter Adenovirus 2 Major Late Promoter (mlp) is formed by the apposition of two sequences with very different dynamic properties: an alternating TATA sequence and an A-tract. For a comparative study, we carried out molecular dynamics simulations of two DNA oligomers, one containing the mlp sequence (2 ns), and the other an analog where AT basepairs were substituted by IC basepairs (1 ns). The simulations, carried out with explicit solvent and counterinons, yield straight purine tracts, the A-tract being stiffer than the I-tract, an alternating structure for the YRYR tracts, and hydration patterns that differ between the purine tracts and the alternating sequence tracts. A detailed analysis of the proposed interactions responsible for the stiffness of the purine tracts indicates that the stacking between the bases bears the strongest correlation to stiffness. The hydration properties of the minor groove in the two oligomers are distinctly different. Such differences are likely to be responsible for the stronger binding of TBP to mlp over the inosine-substituted variant. The calculations were made possible by the development, described here, of a new set of forcefield parameters for inosine that complement the published CHARMM all-hydrogen nucleic acid parametrization.

Selective binding of the TATA box-binding protein to the TATA box-containing promoter: analysis of structural and energetic factors

We report the results of an energy-based exploration of the components of selective recognition of the TATA box-binding protein (TBP) to a TATA box sequence that includes 1) the interaction between the hydrophobic Leu, Pro, and Phe residues of TBP with the TA, AT, AA, TT, and CG steps, by ab initio quantum mechanical calculations; and 2) the free energy penalty, calculated from molecular dynamics/potential of mean force simulations, for the conformational transition from A-DNA and B-DNA into the TA-DNA form of DNA observed in a complex with TBP. The GTAT, GATT, GAAT, and GTTT tetramers were explored. The results show that 1) the discrimination of TA, AT, AA, TT, or CG steps by TBP cannot rest on their interaction with the inserting Phe side chains; 2) the steric clash between the bulky and hydrophobic Pro and Leu residues and the protruding -NH2 group of guanine is responsible for the observed selectivity against any Gua-containing basepair; 3) the Pro and Leu residues cannot selectively discriminate among TA, AT, AA, or TT steps; and 4) the calculated energy required to achieve the TA-DNA conformation of DNA that is observed in the complex with TBP appears to be a key determinant for the observed selectivity against the AT, AA, and TT steps. The simulations also indicate that only the TA step can form a very efficient interbase hydrogen bond network in the TA-DNA conformation. Such an energetically stabilizing network is not achievable in the AA and TT steps. While it is viable in the AT step, structural constraints render the hydrogen bonding network energetically ineffective there.

Inherent DNA curvature and flexibility correlate with TATA box functionality

Four 1.5 ns molecular dynamics (MD) simulations were performed on the d(GCTATAAAAGGG).d(CCCTTTTATAGC) double helix dodecamer bearing the Adenovirus major late promoter TATA element and three iso-composition mutants for which physical and biochemical data are available from the same laboratory. Three of these DNA sequences experimentally induce tight binding with the TATA box binding protein (TBP) and induce high transcription rates; the other DNA sequence induces much lower TBP binding and transcription. The x-ray crystal structures have previously shown that the duplex DNA in DNA-TBP complexes are highly bent. We performed and analyzed MD simulations for these four DNAs, whose experimental structures are not available, in order to address the issue of whether inherent DNA structure and flexibility play a role in establishing these observed preferences. A comparison of the experimental and simulated results demonstrated that DNA duplex sequence-dependent curvature and flexibility play a significant role in TBP recognition, binding, and transcriptional activation.

Sequence-dependent dynamics of TATA-Box binding sites

We have carried out two nanosecond-length molecular dynamics simulations on a DNA oligomer, d(GCGTAAAAAAAACGC)2, which contains a weak binding site for the TATA-box binding protein. An analysis of the resulting trajectories shows that this oligomer behaves differently from a related oligomer [d(GCGTATATAAAACGC)2] studied earlier using the same protocol (Flatters, D., M. Young, D. L. Beveridge, and R. Lavery. 1997. Conformational properties of the TATA-box binding sequence of DNA. J. Biomol. Struct. & Dyn. 14:757-765), and which contains a strong binding site for the same protein. The two basepair mutations that relate these oligomers lead to significant changes in time-averaged structure and in dynamic behavior, which extend over entire length of the oligomer and appear to be compatible with the experimentally observed decrease of binding and functional activity. These results suggest that molecular dynamics simulations, taking into account explicit solvent and counterions, and avoiding the truncation of electrostatic interactions, are a powerful tool for investigating the indirect aspects of protein-nucleic acid recognition.

Progressive DNA bending is made possible by gradual changes in the torsion angle of the glycosyl bond

Structural comparisons have led to the suggestion that the conformational rearrangement that would be required to change A-DNA into the TA-DNA form of DNA observed in the complex with the TATA box binding protein (TBP) could be completed by modifying only the value of the glycosyl bond chi by approximately 45 degrees. The lack of a high number of crystal structures of this type makes it difficult to conclude whether a smooth transition from A-DNA to TA-DNA can occur without disrupting at any point either the Watson-Crick base pairing or the A-DNA conformation of the backbone. To explore the possibility of such a smooth transition, constrained molecular dynamics simulations were carried out for the double-stranded dodecamer d(GGTATATAAAAC), in which a transition from A-DNA to TA-DNA was induced by modifying only the chi angle values. The results demonstrate the feasibility of a continuous path in the A-DNA to TA-DNA transition. Varying extents of DNA curvature are also attainable, by maintaining the A-DNA backbone structure and Watson-Crick hydrogen bonding while changing the chi angle value smoothly from that in A-DNA to one corresponding to B-DNA.