Guanine tracts enhance sequence directed DNA bends

Large-scale chromosome folding versus genomic DNA sequences: A discrete double Fourier transform technique

Using state-of-the-art techniques combining imaging methods and high-throughput genomic mapping tools leaded to the significant progress in detailing chromosome architecture of various organisms. However, a gap still remains between the rapidly growing structural data on the chromosome folding and the large-scale genome organization. Could a part of information on the chromosome folding be obtained directly from underlying genomic DNA sequences abundantly stored in the databanks? To answer this question, we developed an original discrete double Fourier transform (DDFT). DDFT serves for the detection of large-scale genome regularities associated with domains/units at the different levels of hierarchical chromosome folding. The method is versatile and can be applied to both genomic DNA sequences and corresponding physico-chemical parameters such as base-pairing free energy. The latter characteristic is closely related to the replication and transcription and can also be used for the assessment of temperature or supercoiling effects on the chromosome folding. We tested the method on the genome of E. coli K-12 and found good correspondence with the annotated domains/units established experimentally. As a brief illustration of further abilities of DDFT, the study of large-scale genome organization for bacteriophage PHIX174 and bacterium Caulobacter crescentus was also added. The combined experimental, modeling, and bioinformatic DDFT analysis should yield more complete knowledge on the chromosome architecture and genome organization.

Recent findings on the mechanisms involved in tenofovir resistance

Since its approval for clinical use in 2001, tenofovir (TFV) has become one of the most frequently prescribed nucleotide analogues used in combination with other antiretroviral agents against HIV-1 infection. Although reverse transcriptase inhibitors (RTIs) including TFV have been shown to be highly potent with reasonable safety profiles in the clinic, drug resistance hinders the effectiveness of current therapies and even causes treatment failure. Therefore, understanding the resistance mechanisms of RT and exploring the potential antiviral synergy between the different RTIs in combination therapies against the resistance mechanisms would greatly improve the long-term efficacy of existing and future regimens. We have studied the pyrophosphorolytic removal of TFV, a major resistance mechanism that RT utilizes, from two different viral sequences and observed interesting outcomes associated with the sequence context. Furthermore, addition of efavirenz, a non-nucleoside RTI, inhibits this removal process confirming the synergistic antiviral effects. This article highlights our recently published work on the viral sequence context contributing to the study of anti-HIV drug resistance in conjunction with the benefits of combining various RTIs that may have been neglected previously.

Monovalent cation binding by curved DNA molecules containing variable numbers of a-tracts

Monovalent cation binding by DNA A-tracts, runs of four or more contiguous adenine or thymine residues, has been determined for two curved approximately 200 basepair (bp) restriction fragments, one taken from the M13 origin of replication and the other from the VP1 gene of SV40. These two fragments have previously been shown to contain stable, centrally located bends of 44 degrees and 46 degrees , respectively, located within approximately 60 bp "curvature modules" containing four or five irregularly spaced A-tracts. Transient electric birefringence measurements of these two fragments, sequence variants containing reduced numbers of A-tracts in the SV40 curvature module or changes in the residues flanking the A-tracts in the M13 curvature module, have been combined with the free solution electrophoretic mobilities of the same fragments using known equations to estimate the effective charge of each fragment. The effective charge is reduced, on average, by one-third charge for each A-tract in the curvature module, suggesting that each A-tract binds a monovalent cation approximately one-third of the time. Monovalent cation binding to two or more A-tracts is required to observe significant curvature of the DNA helix axis.

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.

Signals for TBP/TATA box recognition

The TATA box-binding protein (TBP) recognizes its target sites (TATA boxes) by indirectly reading the DNA sequence through its conformation effects (indirect readout). Here, we explore the molecular mechanisms underlying indirect readout of TATA boxes by TBP by studying the binding of TBP to adenovirus major late promoter (AdMLP) sequence variants, including alterations inside as well as in the sequences flanking the TATA box. We measure here the dissociation kinetics of complexes of TBP with AdMLP targets and, by phase-sensitive assay, the intrinsic bending in the TATA box sequences as well as the bending of the same sequence induced by TBP binding. In these experiments we observe a correlation of the kinetic stability to sequence changes within the TATA recognition elements. Comparison of the kinetic data with structural properties of TATA boxes in known crystalline TBP/TATA box complexes reveals several "signals" for TATA box recognition, which are both on the single base-pair level, as well as larger DNA tracts within the TATA recognition element. The DNA bending induced by TBP on its binding sites is not correlated to the stability of TBP/TATA box complexes. Moreover, we observe a significant influence on the kinetic stability of alteration in the region flanking the TATA box. This effect is limited however to target sites with alternating TA sequences, whereas the AdMLP target, containing an A tract, is not influenced by these changes.

Compression of the DNA minor groove is responsible for termination of DNA synthesis by HIV-1 reverse transcriptase

HIV-1 reverse transcriptase (RT) generally terminates plus strand DNA synthesis at the centre of the viral genome. The central termination sequence (CTS) contains two termination sites which are located at the 3' end of AnTm motifs. These motifs generate a global curvature of the DNA helix which correlates with termination of DNA synthesis. Here, we have characterized HIV-1 RT termination sites on different DNA sequences. Again, they are located at the 3' end of A-tracts. Using hydroxyl radicals as a probe of the width of the DNA helix, we have shown that RT termination sites are always located a few base-pairs downstream of a compressed minor groove. Mutations which relieve these compressions also abolish the termination events. The replacement of the adenine tracts by 2,6-diaminopurine tracts has a similar effect. Moreover, no termination site is observed on DNA sequences containing phased GC-tracts which curve the DNA helix but compress the major groove. The compression of the DNA minor groove and not necessarily the curved trajectory of the DNA is, therefore, responsible for termination of DNA synthesis at the CTS by HIV-1 RT. This conclusion is consistent with interpretation of other biochemical data on the processivity of HIV-1 RT, based on the structure of a DNA-enzyme complex.

Small circles of helical DNA obtained on the basis of transcriptional pause sites sequence

The 21-base pair synthetic DNA duplexes with basic 'pause-motif site ('CATGC') were ligated head-to-tail to produce linear and circular multimers. This also was done from other closely related sequences. Electrophoretic mobilities of the linear multimers in polyacrylamide gels were determined under the standard and modified conditions. We revealed that small linear multimers (approximately 90 bp) were characterized by comparable value of gel retardation relative to the well known curved DNA, while longer multimers (130 to approximately 170 bp) had only slightly expressed mobility anomaly. Nevertheless these multimers containing nontruncated 'pause-motif were capable of cyclization, in particular, formation of unusually small circles while truncated ones were not. We conclude that basic 'pause-motif site increases the closure ability while the multimers based on truncated 'pause motif fail to curve into the small circles. We tend to explain this situation as a result of intrinsic bending as well as the influence of the thermal fluctuations of DNA, the latter most probably can be associated with 'pause motif'. We have estimated the equilibrial and maximal bend angles per 10.5 bp to be 12 degrees to approximately 16 degrees and 32 degrees accordingly under experimental conditions of our study.

DNA structural forms

In the years that have passed since the publication of Wolfram Saenger's classic book on nucleic acid structure (Saenger, 1984), a considerable amount of new data has been accumulated on the range of conformations which can be adopted by DNA. Many unusual species have joined the DNA zoo, including new varieties of two, three and four stranded helices. Much has been learnt about intrinsic DNA curvature, dynamics and conformational transitions and many types of damaged or deformed DNA have been investigated. In this article, we will try to summarise this progress, pointing out the scope of the various experimental techniques used to study DNA structure, and, where possible, trying to discern the rules which govern the behaviour of this subtle macromolecule. The article is divided into six major sections which begin with a general discussion of DNA structure and then present successively, B-DNA, DNA deformations, A-DNA, Z-DNA and DNARNA hybrids. An extensive set of references is included and should serve the reader who wishes to delve into greater detai.

The organic crystallizing agent 2-methyl-2,4-pentanediol reduces DNA curvature by means of structural changes in A-tracts

Contemporary predictive models for sequence-dependent DNA structure provide a good estimation of overall DNA curvature in most cases. However, the two current models differ fundamentally in their view of the origin of DNA curvature. An earlier model that associates DNA bending primarily, although not exclusively, with stretches of adenines (A-tracts) is based on results of comparative gel retardation, cyclization kinetics, hydroxyl radical cutting, and other solution measurements. It represents an intersection of wedge and junction models. More recently, a non-A-tract bending model has been proposed, built on structural results from x-ray crystallography and molecular modeling. In this view, A-tracts are proposed to be straight and rigid, whereas mixed sequence DNA is bent. Because a key premise of the non-A-tract bending model is the crystallographic observation that A-tracts are straight, we have examined the effect in solution of 2-methyl-2,4-pentanediol (MPD), an organic solvent used in crystal preparation for crystallographic DNA structure determinations. Using cyclization analysis, DNase I cutting, chemical probing, and electron microscopy on DNA oligomers with and without A-tracts, we show that the presence of MPD in solution dramatically affects A-tracts and that the effect is specific to these sequence elements. Combined with the previous observation that MPD affects gel mobility of curved sequences with A-tracts, our findings support the bent A-tract model and call for caution in the interpretation of crystallographic results on DNA structure as these are presently obtained.

Bending and torsional flexibility of G/C-rich sequences as determined by cyclization assays

The structural polymorphism of DNA is a vital aspect of its biological function. However, it has become increasingly apparent in recent years that DNA polymorphism is a complicated, multidimensional phenomenon that includes not only static sequence-directed structures but dynamic effects as well, including influences of counterions and sequence context. In order to address some of these additional factors that govern DNA conformation, we have used T4 ligase-mediated cyclization to investigate bending in a series of DNA sequences containing the GGGCCC.GGGCCC motif in different sequence contexts including various helical phasings with (A)5-tracts. We present evidence for curvature in GGGCCC.GGGCCC and (A)5-tract motifs in the presence of physiological levels of Mg2+ and show that these motifs curve through similar but oppositely directed bending angles under these ionic strength conditions. Although these two sequence motifs appear to bend similarly, our results suggest significant differences in stiffness and stability of curvature between them. We also show that under the same experimental conditions, the CTAG-CTAG sequence element possesses unusual torsional flexibility and that this appears to be associated with the central TA.TA dinucleotide. The results underscore the need to include sequence context and specific ion effects as well as a dynamic basis in more complete predictive models for functionally related DNA polymorphism.

A single-stranded DNA binding protein binds the origin of replication of the duplex kinetoplast DNA

Replication of the kinetoplast DNA (kDNA) minicircle of trypanosomatids initiates at a conserved 12-nt sequence, 5'-GGGGTTGGTGTA-3', termed the universal minicircle sequence (UMS). A sequence-specific single-stranded DNA-binding protein from Crithidia fasciculata binds the heavy strand of the 12-mer UMS. Whereas this UMS-binding protein (UMSBP) does not bind a duplex UMS dodecamer, it binds the double-stranded kDNA minicircle as well as a duplex minicircle fragment containing the origin-associated UMS. Binding of the minicircle origin region by the single-stranded DNA binding protein suggested the local unwinding of the DNA double helix at this site. Modification of thymine residues at this site by KMnO4 revealed that the UMS resides within an unwound or otherwise sharply distorted DNA at the minicircle origin region. Computer analysis predicts the sequence-directed curving of the minicircle origin region. Electrophoresis of a minicircle fragment containing the origin region in polyacrylamide gels revealed a significantly lower electrophoretic mobility than expected from its length. The fragment anomalous electrophoretic mobility is displayed only in its native conformation and is dependent on temperature and gel porosity, indicating the local curving of the DNA double helix. We suggest that binding of UMSBP at the minicircle origin of replication is possible through local unwinding of the DNA double helix at the UMS site. It is hypothesized here that this local melting is initiated through the untwisting of unstacked dinucleotide sequences at the bent origin site.

Oligo[d(C).(G)] runs exhibit a helical repeat of 11.1 bp in solution and cause slight DNA curvature when properly phased

We have inserted d(C)10 in a set of DNA fragments with bent segments on both ends, which are rotated with respect to each other by base pair wise increasing insertions. The electrophoretic mobilities on polyacrylamide gels of these DNA fragments were used to identify insertion sizes with cis conformations of the bent ends. These experiments revealed a helical repeat in solution of d(C).d(G) tracts of 11.1 +/- 0.08 bp. The electrophoretic mobilities of ligation ladders with properly phased d(C)5 and d(C)16 runs demonstrate a small but clearly detectable curvature of these fragments.

Strong sequence patterns in eukaryotic promoter regions: potential implications for DNA structure

1. Analysis of eukaryotic sequences reveals recurring trends in upstream regions. Oligomers composed of (G/C)n and (A/T)m blocks are preferentially flanked by (G/C)2 doublets on their 3' rather than on their 5' ends, that is (G/C)n(A/T)m(G/C)2 > (G/C)n+2(A/T)m. 2. These trends are stronger for larger n and smaller m. Additional trends are outlined below. 3. The trends are correlated with DNA structural parameters, in particular with twist and roll angles. 4. Generally, the trends hold if the base pair step joining the 5' (G/C)2 doublet to the (G/C)n (A/T)m oligomer is not undertwisted and is not strongly rolled into the major groove. 5. Other DNA parameters crucial for DNA-protein interactions are discussed as well.

Repeat units from a maize rDNA external spacer region exhibit DNA curvature and interact with high-mobility-group proteins

The 3-kb external spacer from a maize (Zea mays L. cv. A619) nuclear rRNA gene unit which contains nine highly homologous 200-bp repeat elements was found to include a region with DNA-curvature properties. The centre of curvature was localized within repeats 5 and 6 using a circular permutation assay. A 60-bp-long subfragment of this region was found to interact with nuclear proteins, including high-mobility-group (HMG) proteins, and with the maize HMGa protein synthesized in Escherichia coli from a recombinant plasmid. The potential influence of the binding of the HMG proteins on the conformation of this subfragment was studied with a permutation assay based on a bending vector.

In vitro evolution of intrinsically bent DNA

DNA fragments which are intrinsically bent or curved migrate anomalously during electrophoresis through polyacrylamide gels. Starting with an initial population of approximately 10(12) unique DNA sequences, DNA which exhibited the kind of anomalous mobility associated with DNA bending was selected and enriched using a variation of the SELEX procedure. After seven rounds of selection and amplification, the vast majority of the remaining population of DNA fragments migrated as bent DNA. Cloning and sequencing of 30 individual sequences from this population has yielded information regarding the relationship between DNA sequence and bending. Some of the previous conclusions on DNA bending have been confirmed while others have been modified, by the results presented here. In addition, the dinucleotide base step CA/TG, which had not been thought to be a major factor in DNA bending, appears to be important.

Trends in the 5' vs. 3' flanks of oligonucleotides in eukaryotic and prokaryotic genomes: the asymmetric roles played by cytosine and guanine

Studies of sequence context preferences of oligonucleotides composed of (G/C)n and (A/T)m blocks (n + m = 3,4,5) unravel strong patterns. Comparisons of the 5' and 3' nearest neighbor doublets flanking these oligomers reveal the preference of (G/C)2 to be positioned immediately next to the (A/T)m block, enclosing it by (G/C) nucleotides rather than extending the (G/C)n block. That is, for a (G/C)n(A/T)m oligomer and a (G/C)2 doublet, (G/C)n(A/T)m(G/C)2 greater than (G/C)n + 2 (A/T)m. Similarly for an (A/T)m(G/C)n oligomer, (G/C)2(A/T)m(G/C)n greater than (A/T)m(G/C)n + 2. In an analogous manner, (A/T)2 flanking doublets prefer enclosing the (G/C)n blocks, although these patterns are weaker. Here we show a strong, direct relationship between the magnitude of the trends and the presence of Cs in the (G/C)n block in the (G/C)n(A/T)m oligomer, and the presence of Gs in the complementary (A/T)m(G/C)n oligomers. The trends are stronger in eukaryotic than in prokaryotic sequences. They are stronger for longer (G/C)n and shorter (A/T)m blocks. We suggest that the preference for (A/T)m to be enclosed by (G/C) rather than be flanked by them on only one side is related to DNA structure and DNA-protein interaction. Sequences of the (G/C)(A/T)(G/C) type may have more homogeneous minor groove geometry. In particular, the strong G vs. C asymmetry in the trends may be related to pyrimidine-purine junctions, possibly to CG sequences.

Curved DNA without AA/TT dinucleotide step

The evidence is accumulating that dinucleotide steps other than AA/TT affect DNA flexure of AnTm (m + n greater than = 4) containing fragments. However, it is not clear whether macroscopic DNA flexure without AA/TT steps might occur. In this paper we demonstrate the anomaly in electrophoretic mobility of non AA/TT repetitive DNA sequences which is a function of sequence phasing. Therefore, our results show that PyPu (TA) and AG/CT steps, angulary separated by close to 180 degrees from Pu/Py (GC) and GG/CC steps, bend DNA, even in the absence of AnTm tracts.

Curved DNA without A-A: experimental estimation of all 16 DNA wedge angles

The principal sequence feature responsible for intrinsic DNA curvature is generally assumed to be runs of adenines. However, according to the wedge model of DNA curvature, each dinucleotide step is associated with a characteristic deflection of the local helix axis. Thus, an important test of a more general view of sequence-dependent DNA curvature is whether sequence elements other than A-A cause the DNA axis to deflect. To address this question, we have applied the wedge model to a large body of experimental data. The axial path of DNA can be described at each step by three Eulerian angles: the helical twist, the deflection angle (wedge angle), and the direction of the deflection. Circularization and gel electrophoretic mobility data on 54 synthetic DNA fragments, both from other laboratories and from our own, were used to compare the theoretical predictions of the wedge model with experiment. By minimizing misfit between calculated and observed DNA curvature, we have found that the stacks AG/CT, CG/CG, GA/TC, and GC/GC, in addition to AA/TT, have large wedge values. We have also synthesized seven sequences without AA/TT elements but with these other wedges correctly phased to cause appreciable predicted curvature. All appear curved as demonstrated by anomalous gel mobilities. The full set of 16 roll and tilt wedge angles is estimated and, together with the known 10 helical twists, these allow prediction of the general sequence-dependent trajectory of the DNA axis.

General nearest neighbor preferences in G/C oligomers interrupted by A/T: correlation with DNA structure

The frequencies of occurrence of the 5' and 3' nearest neighbor doublets of oligonucleotides containing (G/C) and (A/T) blocks show strong trends. Specifically, the following trends are observed. Given a (G/C)n (A/T)m oligomer (where G/C)n indicates a sequence of length n composed solely of Gs and/or Cs and (A/T)m is a sequence of length m composed solely of As and/or Ts, and n = 3,2,1; m = 1,2,3) and a (G/mC)2 doublet, (G/C)n (A/T)m (G/C)2 greater than (G/C)n + 2 (A/T)m. That is the (G/C)2 doublet is preferentially located 3' of the oligomer, enclosing the (A/T)m stretch. The trends are strongest for n = 3, m = 1 and gradually weaken as the size of the (mG/C)n block decreases (with a concomitant increase of (A/T)m). (A/T)2 nearest neighbor flank preferentially encloses the (G/C)n block (to produce (A/T)2 (G/C)n (A/T)m). The (A/T)2 flank trends are weaker than the (G/C)2 flank ones. The (A/T)2 flank trends also decrease in strength as the size of the (G/C)n block decreases. The statistical significance of these trends in eukaryotes is very high. A possible correlation with DNA structural parameters, in particular groove geometry, is discussed.