Helical periodicity of DNA, Poly(dA) . poly(dT) and poly(dA-dT). poly(dA-dT) in solution

Poly(dA:dT)-rich DNAs are highly flexible in the context of DNA looping

Large-scale DNA deformation is ubiquitous in transcriptional regulation in prokaryotes and eukaryotes alike. Though much is known about how transcription factors and constellations of binding sites dictate where and how gene regulation will occur, less is known about the role played by the intervening DNA. In this work we explore the effect of sequence flexibility on transcription factor-mediated DNA looping, by drawing on sequences identified in nucleosome formation and ligase-mediated cyclization assays as being especially favorable for or resistant to large deformations. We examine a poly(dA:dT)-rich, nucleosome-repelling sequence that is often thought to belong to a class of highly inflexible DNAs; two strong nucleosome positioning sequences that share a set of particular sequence features common to nucleosome-preferring DNAs; and a CG-rich sequence representative of high G+C-content genomic regions that correlate with high nucleosome occupancy in vivo. To measure the flexibility of these sequences in the context of DNA looping, we combine the in vitro single-molecule tethered particle motion assay, a canonical looping protein, and a statistical mechanical model that allows us to quantitatively relate the looping probability to the looping free energy. We show that, in contrast to the case of nucleosome occupancy, G+C content does not positively correlate with looping probability, and that despite sharing sequence features that are thought to determine nucleosome affinity, the two strong nucleosome positioning sequences behave markedly dissimilarly in the context of looping. Most surprisingly, the poly(dA:dT)-rich DNA that is often characterized as highly inflexible in fact exhibits one of the highest propensities for looping that we have measured. These results argue for a need to revisit our understanding of the mechanical properties of DNA in a way that will provide a basis for understanding DNA deformation over the entire range of biologically relevant scenarios that are impacted by DNA deformability.

Local elasticity of strained DNA studied by all-atom simulations

Genomic DNA is constantly subjected to various mechanical stresses arising from its biological functions and cell packaging. If the local mechanical properties of DNA change under torsional and tensional stress, the activity of DNA-modifying proteins and transcription factors can be affected and regulated allosterically. To check this possibility, appropriate steady forces and torques were applied in the course of all-atom molecular dynamics simulations of DNA with AT- and GC-alternating sequences. It is found that the stretching rigidity grows with tension as well as twisting. The torsional rigidity is not affected by stretching, but it varies with twisting very strongly, and differently for the two sequences. Surprisingly, for AT-alternating DNA it passes through a minimum with the average twist close to the experimental value in solution. For this fragment, but not for the GC-alternating sequence, the bending rigidity noticeably changes with both twisting and stretching. The results have important biological implications and shed light on earlier experimental observations.

Sequence dependence of DNA bending rigidity

For many aspects of DNA-protein interaction, it is vital to know how DNA bending rigidity (or persistence length, a) depends on its sequence. We addressed this problem using the method based on cyclization of short DNA fragments, which allows very accurate determination of a. Our approach was based on assigning specific values of a to each of 10 distinct dinucleotide steps. We prepared DNA fragments, each about 200 bp in length, with various quasi-periodic sequences, measured their cyclization efficiencies (j factors), and fitted the data by the theoretical equation to obtain the values of a for each fragment. From these data, we obtained a set of a for the dinucleotide steps. To test this set, we used it to design DNA sequences that should correspond to very low and very high values of a, prepared the corresponding fragments, and determined their values of a experimentally. The measured and calculated values of a were very close to one another, confirming that we have found the correct solution to this long-standing problem. The same experimental data also allowed us to determine the sequence dependence of DNA helical repeat.

High pressure fourier transform infrared spectroscopy of poly(dA)poly(dT), poly(dA) and poly(dT)

The effect of hydrostatic pressure upon the DNA duplex, poly(dA)poly(dT), and its component single strands, poly(dA) and poly(dT) has been studied by fourier-transform infrared spectroscopy (FT-IR). The spectral data indicate that at 28 degrees C and pressures up to 12 kbar (1200 MPa) all three polymers retain the B conformation. Pressure causes the band at 967 cm(-1), arising from water-deoxyribose interactions, to shift to higher frequencies, a result consistent with increased hydration at elevated pressures. A larger pressure-induced frequency shift in this band is observed in the single stranded polymers than in the double stranded molecule, suggesting that the effect of pressure on the hydration of single strands may be greater than upon a double stranded complex. A pressure-dependent hypochromicity in the bands attributed to base stacking indicates that pressure facilitates the base stacking in the three polymers, in agreement with previous assessments of the importance of stacking in the stabilization of DNA secondary structure at ambient and high pressures.

Three Sequence Rules for Chromatin

Extensive DNA sequence analysis of three eukaryotes, S. cerevisiae, C. elegans, and D. melanogaster, reveals two different AA/TT periodical patterns associated with the nucleosome positioning. The first pattern is the counter-phase oscillation of AA and TT dinucleotides, which has been frequently considered as the nucleosome DNA pattern. This represents the sequence rule I for chromatin structure. The second pattern is the in-phase oscillation of the AA and TT dinucleotides with the same nucleosome DNA period, 10.4 bases. This pattern apparently corresponds to curved DNA, that also participates in the nucleosome formation, and represents the sequence rule II for chromatin. The positional correlations of AA and TT dinucleotides also indicate that the nucleosomes are separated by specific linker sizes (preferably 8, 18, ... bases), dictated by the steric exclusion rules. Thus, the sequence positions of the neighboring nucleosomes are correlated, and this represents the sequence rule III.

Sequence periodicity of Escherichia coli is concentrated in intergenic regions

BACKGROUND

Sequence periodicity with a period close to the DNA helical repeat is a very basic genomic property. This genomic feature was demonstrated for many prokaryotic genomes. The Escherichia coli sequences display the period close to 11 base pairs.

RESULTS

Here we demonstrate that practically only ApA/TpT dinucleotides contribute to overall dinucleotide periodicity in Escherichia coli. The noncoding sequences reveal this periodicity much more prominently compared to protein-coding sequences. The sequence periodicity of ApC/GpT, ApT and GpC dinucleotides along the Escherichia coli K-12 is found to be located as well mainly within the intergenic regions.

CONCLUSIONS

The observed concentration of the dinucleotide sequence periodicity in the intergenic regions of E. coli suggests that the periodicity is a typical property of prokaryotic intergenic regions. We suppose that this preferential distribution of dinucleotide periodicity serves many biological functions; first of all, the regulation of transcription.

Conformational characteristics and correlations in crystal structures of nucleic acid oligonucleotides: evidence for sub-states

Sugar phosphate backbone conformations are a structural element inextricably involved in a complete understanding of specific recognition nucleic acid ligand interactions, from early stage discrimination of the correct target to complexation per se, including any structural adaptation on binding. The collective results of high resolution DNA, RNA and protein/DNA crystal structures provide an opportunity for an improved and enhanced statistical analysis of standard and unusual sugar-phosphate backbone conformations together with corresponding dinucleotide sequence effects as a basis for further exploration of conformational effects on binding. In this study, we have analyzed the conformations of all relevant crystal structures in the nucleic acids data base, determined the frequency distribution of all possible epsilon, zeta, alpha, beta and gamma backbone angle arrangements within four nucleic acid categories (A-RNA and A-DNA, free and bound B-DNA) and explored the relationships between backbone angles, sugar puckers and selected helical parameters. The trends in the correlations are found to be similar regardless of the nucleic acid category. It is interesting that specific structural effects exhibited by the different unusual backbone sub-states are in some cases contravariant. Certain alpha/gamma changes are accompanied by C3' endo (north) sugars, small twist angles and positive values of base pair roll, and favor a displacement of nucleotide bases towards the minor groove compared to that of canonical B form structures. Unusual epsilon/zeta combinations occur with C2' (south) sugars, high twist angles, negative values of base pair roll, and base displacements towards the major groove. Furthermore, any unusual backbone correlates with a reduced dispersion of equilibrium structural parameters of the whole double helix, as evidenced by the reduced standard deviations of almost all conformational parameters. Finally, a strong sequence effect is displayed in the free oligomers, but reduced somewhat in the ligand bound forms. The most variable steps are GpA and CpA, and, to a lesser extent, their partners TpC and TpG. The results provide a basis for considering if the variable and non-variable steps within a biological active sequence precisely determine morphological structural features as the curvature direction, the groove depth, and the accessibility of base pair for non covalent associations.

Preferential binding of quinolones to DNA with alternating G, C / A, T sequences: a spectroscopic study

The binding of quinolones, nalidixic acid (Nal), oxolinic acid (Oxo) with double stranded polynucleotides was undertaken by using UV-melting, UV-Vis absorption, fluorescence and CD spectroscopic techniques. The binding of Nal or Oxo to the polynucleotides under low-salt buffer conditions were determined for poly (dA).(dT), poly [d(A-T)], poly (dG).(dC), poly [d(G-C)] and E. coli DNA. The fluorescence data were analyzed using a previously established two step mechanism with two different DNA-Drug complexes [Rajeswari et al., Biochemistry 26, 6825-31 (1987)]. The first complex [DN](1) with a binding constant K(1), is formed where the interactions are 'nonspecific' and complex [DN](2) with a binding constant K(2), is formed where the interactions are "specific" which involve (additional) hydrophobic type of interactions like 'stacking' of the drug and the overall association constant is represented as K(=K(1)K(2)). The order of binding for Nal and Oxo is: poly [d(G-C)] > poly [d(A- T)] > E. coli > poly (dG).(dC) > poly (dA).(dT). Interaction of quinolones seems to be preferential in the alternating G, C or A, T stretches of DNA than those of non-alternating. Within any alternating or non-alternating in DNA sequences the G, C rich sequences have distinctly greater binding than A, T sequences. The overall association constant data (K) reveal higher binding of Oxo to DNA compared to Nal to any given polynucleotide investigated; which also explains the higher antibacterial potency of Oxo. Changes in the absorption difference spectra and in circular dichroic spectra also manifest these results. As the melting temperatures of the polynucleotides were only marginally raised in presence of the quinolone, we rule out the possibility of 'classical intercalation' of the drug. Amino group of guanine facilitates the binding of quinolones and therefore has the greater binding with the DNA. However, poly (dG).(dC) is known to exist in 'A' conformation which is not adopted by quinolones as in the case of poly (dA).(dT). Present results suggest that Nal or Oxo bind to DNA in a non-classical fashion which is partially stacking in nature.

Critical effect of the N2 amino group on structure, dynamics, and elasticity of DNA polypurine tracts

Unrestrained 5-20-ns explicit-solvent molecular dynamics simulations using the Cornell et al. force field have been carried out for d[GCG(N)11GCG]2 (N, purine base) considering guanine*cytosine (G*C), adenine*thymine (A*T), inosine*5-methyl-cytosine (I*mC), and 2-amino-adenine*thymine (D*T) basepairs. The simulations unambiguously show that the structure and elasticity of N-tracts is primarily determined by the presence of the amino group in the minor groove. Simulated A-, I-, and AI-tracts show almost identical structures, with high propeller twist and minor groove narrowing. G- and D-tracts have small propeller twisting and are partly shifted toward the A-form. The elastic properties also differ between the two groups. The sequence-dependent electrostatic component of base stacking seems to play a minor role. Our conclusions are entirely consistent with available experimental data. Nevertheless, the propeller twist and helical twist in the simulated A-tract appear to be underestimated compared to crystallographic studies. To obtain further insight into the possible force field deficiencies, additional multiple simulations have been made for d(A)10, systematically comparing four major force fields currently used in DNA simulations and utilizing B and A-DNA forms as the starting structure. This comparison shows that the conclusions of the present work are not influenced by the force field choice.

Invariant and variable base stacking geometries in B-DNA and A-DNA

We calculated the interatomic distances between all couples of non-hydrogen atoms belonging to the neighboring Watson-Crick base pairs in the available crystal structures of DNA. Their standard deviations revealed remarkably large differences in the variability of the base stacking geometries of the particular steps. In line with experimental studies in solution, (CpA)-(TpG) and (TpA).(TpA) were identified as the most variable or flexible steps in the crystal structures of B-DNA. On the other hand, base stacking geometries of the (ApT).(ApT) steps were the most invariant, which was very surprising because all three steps composed only of C and G were much more flexible. This finding suggests that conformational stability of DNA and the rigidity have different origins. Furthermore, the nucleotide sequence dependence of the flexibility was almost reversed in A-DNA because the most flexible steps in B-DNA were the least flexible in A-DNA. The most invariant steps of B-DNA were variable in A-DNA. The (ApT).(ApT) step was a notable exception to this rule because it belonged to the most rigid steps in both B-DNA and A-DNA. The present results are fully consistent with the properties that poly(dA-dT).poly(dA-dT), poly(dA).poly(dT), poly(dAdC).poly(dG-dT) and poly(dA-dG).poly(dC-dT) exhibit in solution.

Molecular dynamics and continuum solvent studies of the stability of polyG-polyC and polyA-polyT DNA duplexes in solution

Molecular dynamics simulation in explicit solvent and continuum solvent models are applied to investigate the relative stability of A- and B-form helices for two DNA sequences, dA10-dT10 and dG10-dC10 in three structural forms. One structural form is based on an unrestrained molecular dynamics (MD) trajectory starting from a canonical B-DNA structure, the second is based on a MD trajectory starting in a canonical B-DNA structure with the sugars constrained to be C2'-endo and the third simulation started from a canonical A-DNA structure with the sugars constrained to C3'-endo puckers. For the energetic analysis, structures were taken as snapshots from nanosecond length molecular dynamics simulations computed in a consistent fashion in explicit solvent, applying the particle mesh Ewald method and the Cornell et al. force field. The electrostatic contributions to solvation free energies are computed using both a finite-difference Poisson-Boltzmann model and a pairwise Generalized Born model. The non-electrostatic contributions to the solvation free energies are estimated with a solvent accessible surface area dependent term. To estimate the gas phase component of the relative free energy between the various structures, the mean solute internal energies (determined with the Cornell et al. molecular mechanics potential including all pairwise interactions within the solute) and estimates of the solute entropy (using a harmonic approximation) were used. Consistent with experiment, the polyG-polyC (GC) structures are found to be much more A-phillic than the polyA-polyT (AT) structures, the latter being quite A-phobic. The dominant energy components responsible for this difference comes from the internal and van der Waal energies. A perhaps less appreciated difference between the GC and AT rich sequences is suggested by the calculated salt dependence which demonstrates a significantly enhanced ability to drive GC rich sequences towards an A-form structure compared to AT rich sequences. In addition to being A-phobic, the AT structure also has a noticably larger helical repeat than GC and other mixed sequence duplexes, consistent with experiment. Analysis of the average solvent density from the trajectories shows hydration patterns in qualitative agreement with experiment and previous theoretical treatments.

On the recombinational origin of protein-sequence-subunit structure

Since 1929 the concept that proteins are built from subunits of certain standard size (Svedberg 1929) has been revisited several times, each time with a new demonstration that, indeed, there are certain preferred protein sizes. According to recent estimates the overrepresented sizes are close to multiples of 125 amino acid (aa) residues for eukaryotes and 150 residues for prokaryotes. To explain these preferences, a hypothesis is suggested, and quantitatively developed, on the recombinational nature of this regularity. The protein-coding sequences are assumed to evolve at some early stage via recombinational events--insertions of DNA circles of a certain optimal size. The contour lengths of the protein-coding DNA circles had to be simultaneously divisible by three and, to minimize torsional constraint, by the DNA helical repeat. With these two conditions satisfied, the calculated contour lengths of the DNA circles, 250-500 base pairs (bp), turn out to correspond well to known optimal DNA circularization sizes and to the predicted range of the protein sequence subunit sizes: 80-170 aa residues, which covers experimentally observed values. The subunit size is found to be strongly influenced by the helical repeat of DNA. The sizes 125 and 150 aa are derived when the corresponding helical repeats of DNA are set within fractions of promilles from the 10.54 bp/turn value. This fits to the experimentally estimated mean for natural mixed DNA sequences, 10.53-10.57 bp/turn.(ABSTRACT TRUNCATED AT 250 WORDS)

Systematic study of nuclear Overhauser effects vis-à-vis local helical parameters, sugar puckers, and glycosidic torsions in B DNA: insensitivity of NOE to local transitions in B DNA oligonucleotides due to internal structural compensations

A method has been developed to solve structures of DNA oligomers in solution from the experimental NOE data. The method is a combination of two approaches: (1) full matrix NOESY simulations and (2) conformational calculations of DNA double helix based on generalized helical parameters. The process of the refinement of a solution structure does not involve NMR-derived interproton distance constraints; rather it consists of a direct fitting of a structure to the experimental NOE data, a weighted sum of energy, and R factor being under minimization. A helical parameters-based generation of DNA forms makes it possible to organize the search for the optimal structure more effectively, systematically varying starting conformations. The method has been used to calculate a structure for the self-complementary DNA hexamer GGATCC, which is consistent with the available experimental data. The structure belongs to the B family of forms, although the local structural heterogeneity is very strong. Sugar puckers vary from O4'-exo to C3'-exo; helical steps are open with different magnitudes toward the minor groove. Next, we have addressed the question of how uniquely the structure is defined by the existing NMR data. Different structural parameters have been systematically varied, and their effect on individual NOE's and the R factor has been studied. Two energetically conjugated parameters, sugar puckers and glycosidic angles, can be determined very reliably, because of the strong dependences of the intraresidue H6/H8 to H2'/H2''/H3' NOE's. In contrast, the local helical conformation of DNA and the geometry of base pairs proved to be underdetermined by the existing NOE information, because the effect of any helical parameter on interproton distances can be compensated by the concerted changes in other parameters.

Structural alteration in alternating adenine-thymine sequences in positively supercoiled DNA

An alternating adenine-thymine tract in a relaxed closed circular plasmid was found to become strongly reactive to osmium tetroxide in the presence of actinomycin D. We suggest that this is due to a local overwinding of the alternating tract as a result of positive supercoiling induced by intercalation of the antibiotic at GpC sequences elsewhere in the DNA. We have previously shown that (A.T)n sequences undergo a local underwinding in response to negative supercoiling, and it appears that such sequences are especially torsionally deformable in both directions.

Actinomycin D induced DNase I hypersensitivity and asymmetric structure transmission in a DNA hexadecamer

DNase I cleavage rates and nmr chemical shifts are shown to change for DNA sequences distal to an intercalated actinomycin D molecule in a duplex hexadecamer upon drug binding. Both sets of observations suggest that the source of these changes is a DNA-mediated structural response. The nmr results imply the response is transmitted preferentially in a 5'-to-3' direction from the drug binding site. An inequivalent response of the two strands to a ligand-induced conformational change immediately suggests a mechanism for distinguishing the sense and antisense strands of DNA.

DNA-sequence and metal-ion specificity of the formation of *H-DNA

The homopyrimidine-homopurine sequence d(CT/GA)22 undergoes, in the presence of zinc ions, transition to an altered DNA conformation (*H-DNA) which is neither H-DNA nor B-DNA. *H-DNA is characterized by a peculiar chemical reactivity pattern in which most of the polypyrimidine strand is hyperreactive to osmium tetroxide and the central part of the polypurine strand is sensitive to diethylpyrocarbonate. Formation of *H-DNA is specific of metal-ion. *H-DNA is detected in the presence of Zn++, Cd++ and Mn++. The efficiency on promoting the transition is in the order of Zn++ greater than Cd++ much greater than Mn++. Formation of *H-DNA is also specific of nucleotide sequence. From all the different homopolymeric sequences tested only the d(CT/GA)22 sequence showed the zinc-induced transition to *H-DNA. These results suggest that stabilization of *H-DNA involves the formation of a specific complex between the metal-ion and the nucleotide sequence. The biological relevance of these results is discussed in view of the important role that zinc ions play on many nucleic acids processes.

DNA B to D transition can be explained in terms of hydration economy of the minor groove atoms

Adjacent phosphate oxygen atoms in A and Z-DNA are located much closer together than in the B form and can be hydrated more economically due to the formation of water bridges between them, whereas in the B form phosphates are hydrated individually. This principle of hydration economy of phosphate groups discovered by Saenger and colleagues could not be applied to the B-D transition, which, like the B-A and B-Z transitions, occurs in a situation of water deficiency, because the distances between adjacent phosphates of individual polynucleotide chains in the D form are not much different from B-DNA. It follows from our calculations of B and D-DNA accessibility to solvent performed by the method of Lee & Richards, and from a simulation of solvent structure near DNA, that there is an economy of hydration only for the minor groove atoms. This feature and some experimental data can explain why only a limited range of sequences consisting of A.T or I.C pairs undergo the transition to the D form. The conformational transition in DNAs with such sequences to a poly[d(A]).poly[d(T])-like conformation (Bh-DNA), which is accompanied by a narrowing of the minor groove, can be explained in the same way. Calculations suggest that in the D-form minor groove of different A-T or I-C DNAs there is a double-layer hydration spine similar to that observed by Drew & Dickerson in the A-T tract of the d(C-G-C-G-A-A-T-T-C-G-C-G) dodecamer. The B-D and B-Bh transitions in A + T-rich DNAs can have biological implications, e.g. they can facilitate DNA bending upon the interaction with proteins.

Structure of the beta-form of poly d(A).poly d(U)

The crystalline beta-form of the sodium salt of poly d(A).poly d(U) trapped in oriented fibers forms a Watson-Crick base-paired, 10(1) double-helix of pitch 3.2 nm. Two molecules are present in a monoclinic unit cell apparently isomorphous with beta-poly d(A).poly d(T). The two chains in each molecule both carry C2'-endo puckered furanose rings but are conformationally not identical. The orientations of the A:U base-pairs relative to the helix-axis are distinctly different from those in classical B-DNA and the overall morphology of the duplex in which they reside resembles that of the alpha-forms of poly (purine).poly (pyrimidine) DNA duplexes previously reported.

Sequence and structural requirements of a herpes simplex viral DNA replication origin

Herpes simplex virus (HSV) types 1 and 2 contain two classes of origins of DNA replication, oriS and oriL, which are closely related. A series of plasmids was constructed which contained specifically altered versions of the HSV type 2 oriS replication origin. Their ability to replicate in an in vivo replicon assay allowed a core origin of 75 base pairs (bp) to be defined. It included both arms of a 56-bp palindrome and from 13 to 20 bp of sequence leftward of the palindrome. The AT-rich sequence at the center of the palindrome was essential. Sequences on either side of the core origin enhanced replication. When additional copies of the -AT-dinucleotide were introduced progressively into the center of the palindrome, an oscillating effect on origin function was observed. These and other data implicate a linear rather than a cruciform conformation of the oriS palindrome in the initiation of HSV replication.

Actinomycin D induced DNase I cleavage enhancement caused by sequence specific propagation of an altered DNA structure

Two DNA hexadecamers containing one central 5'-GC-3' base step have been examined by footprinting methodology in the presence and absence of actinomycin D. The results of these studies, coupled with imino proton NMR measurements indicate that the antitumor drug causes a change in DNA conformation at a distance from the actinomycin intercalation site in a molecule of sequence d[ATATATAGCTATATAT] that does not occur in d[AAAAAAAGCTTTTTTT]. The experiments demonstrate that DNase I rate enhancements associated with actinomycin D binding are caused by ligand alteration of equilibrium DNA structure.

A rigorous basepair oriented description of DNA structures

We propose new, rigorous definitions for (i) basepair fixed coordinate systems and (ii) the twist, tilt, and roll angles (called tau, t, rho) describing the relative orientation of adjacent basepairs and bases in a pair, in arbitrary DNA structures obtained from x-ray diffraction, 2D NMR, or energy calculations. In contrast to the corresponding angular parameters (tg, theta T, theta R) and coordinate systems introduced by Dickerson and co-workers and currently in use, our angular parameters and coordinate systems, together with a set of three displacement parameters, dx, dy, dz, provide a mathematically correct and general description of DNA conformations at the basepairs and/or base level. For instance, our description is applicable when the DNA structure considered is inherently curved, irregular, and/or does not possess dyad (or pseudodyad) axes. We develop a computationally convenient algorithm for rigorous DNA conformational analysis and apply it to some of the known crystal structures. We establish the connection to the currently used parameters and test the consistency and efficiency of our methodology by reconstructing the Dickerson B dodecamer using only the sequence and the set of parameters obtained from the atomic coordinates. The six parameter (tau, t, rho, dx, dy, dz) basepair level reconstruction is good but not perfect. Perfect reconstruction is obtained when one also considers each base in a basepair (consideration of propeller twist alone is not sufficient). The variation of the rigorous parameters proposed along the sequence is much larger, but their average values agree with fiber and solution data much better than in the case of the currently used set. The results of our analysis do not support Trifonov's AA.TT wedge model for DNA curvature but provide some evidence in favor of the Crothers junction-bend model. We point out some of the limitations of basepair level approaches when applied to DNA structure prediction and quantitative understanding of sequence-dependent variations in structure.

Sequence and structural requirements of a herpes simplex viral DNA replication origin

Herpes simplex virus (HSV) types 1 and 2 contain two classes of origins of DNA replication, oriS and oriL, which are closely related. A series of plasmids was constructed which contained specifically altered versions of the HSV type 2 oriS replication origin. Their ability to replicate in an in vivo replicon assay allowed a core origin of 75 base pairs (bp) to be defined. It included both arms of a 56-bp palindrome and from 13 to 20 bp of sequence leftward of the palindrome. The AT-rich sequence at the center of the palindrome was essential. Sequences on either side of the core origin enhanced replication. When additional copies of the -AT-dinucleotide were introduced progressively into the center of the palindrome, an oscillating effect on origin function was observed. These and other data implicate a linear rather than a cruciform conformation of the oriS palindrome in the initiation of HSV replication.

Structural basis of stable bending in DNA containing An tracts. Different types of bending

Structural determinants of DNA bending of different types have been studied by theoretical conformational analysis of duplexes. Their terminal parts were fixed either in an ordinary low-energy B-like conformation or in "anomalous" conformations with a narrowed minor groove typical of An tracts. The anomalous conformations had different negative tilt angles (up to about zero), different propeller twists and minor groove widths. Calculations have been performed for DNA fragments AnTm, TnAm, AnGCTm, AnCGTm, TmGCAn, TmCGAn which are the models of the junction of two anomalous structures on An and Tm tracts. On the AT step of the AnTm fragment the minor groove can be easily narrowed so that a whole unbent fragment of anomalous structure is formed on AnTm. According to our energy estimates, there should not be any reliable bending on AnTm. In contrast, in all other cases there was a pronounced roll-like bending into the major groove in the chemical symmetry region. Calculations of the junction between the anomalous and ordinary B-like structure for GnTm and CnAm have shown that there is an equilibrium bending with a tilt component towards the chain having the anomalous structure at the 5'-end. From our calculations it is impossible to determine precisely the direction of bending, though it can be suggested that the roll component of bending might be directed towards the major groove. The anomalous structure is the main reason of bending; alternations of pyrimidines and purines can modulate the value and the direction of equilibrium bending (only the value in the case of self-complementary fragments).(ABSTRACT TRUNCATED AT 250 WORDS)

Chromatin reconstitution on small DNA rings. I

Chromatin was reconstituted using the four core histones on 359 base-pair nicked and closed rings by salt dialysis and/or at physiological ionic strength by means of polyglutamic acid. The products, which consisted of mono- and dinucleosomes, were characterized by gel electrophoresis, sedimentation in sucrose gradients and high-resolution electron microscopy. The results were as follows. (1) The efficiency of the reconstitution was found first to increase with the negative linking difference of the closed rings relative to their relaxed configuration to reach a maximum for -2 turns, and then to decrease for the largest difference of -3 turns. Discrepancies between topoisomers were also observed with regard to differential formation of mono- and dinucleosomes. Topoisomer -1 reconstituted monomers easily but reconstituted dimers with difficulty, whilst this discrimination was virtually absent in the case of topoisomers -2 and -3. Moreover, mononucleosomes on the nicked ring were, with respect to their electrophoretic mobility, similar to mononucleosomes formed on topoisomer -1 but not to those on the other topoisomers, whose mobilities were greater. These features were interpreted in terms of the linking number change associated with the formation of a nucleosome monomer and dimer, approximately -1 and -2 turns, respectively. (2) Two dinucleosome subtypes were found to form in a sequential manner. Their different electrophoretic mobilities and sedimentation coefficients suggested that the early subtype is lighter, probably because of an incomplete histone complement in the second nucleosome of that subtype as a result of an impaired co-operativity in octamer assembly due to the small ring size. (3) An electron microscopic examination of the chromatin reconstituted on topoisomer -2 revealed that both mono- and dinucleosomes adopt two different, salt-dependent, morphologies each: in type I, entering and exiting DNAs do not cross, whilst they do in type II. Type I configuration is favoured in lower salt, whereas type II is favoured in higher salt. Such behaviour explains why nucleosomes in dimers were found to be always diametrically opposed on the rings rather than sometimes apposed, as would have been expected from a random deposition of the histone cores.

Calibration of DNA curvature and a unified description of sequence-directed bending

Chemically synthesized duplex oligodeoxynucleotides having different average numbers of adenine tracts (A6) per helix turn were ligated into multimers and analyzed by electrophoresis on polyacrylamide gels. The magnitude of the anomaly in gel mobility is found to be a quadratic function of the curvature of the DNA molecule. Parameters that describe intrinsic DNA bending, expressed as the tilt and roll components of the helix-axis deflection at the junctions between the adenine tract and adjacent B-DNA, were adjusted to fit the measured relative curvature of regularly repeated DNA bending sequences known from other studies and synthesized for this study. The model developed here retains the predominance of bending in the direction of tilt at the junctions but incorporates an appreciable roll component at the 5' end of an adenine tract, opening the minor groove there. This feature is consistent with chemical "footprinting" experiments on molecules containing adenine tracts. The overall direction of bending is effectively toward the minor groove, viewed from the center of an A5 or A6 tract. A possible underlying structure, which can also be described by a wedge bending model, is that derived from fiber diffraction studies of poly(dA).poly(dT). However, alternative models for the adenine tract, such as propeller twisted DNA, cannot be eliminated, although they do not lead to the correct direction of bending. The results permit calculation of the helix-axis trajectory of natural DNA molecules containing adenine-tract bends.

DNA supercoiling in vivo

DNA topoisomerase mutants of Escherichia coli and Saccharomyces cerevisiae were used to study the topological state of intracellular DNA. In E. coli, it is shown that switching off the gene topA encoding DNA topoisomerase I leads to an increase in the degree of negative supercoiling of intracellular DNA and inhibition of the growth of the cells: a d(pCpG)16.d(pCpG)16 sequence on a plasmid is also shown to flip from a right-handed B-helical structure to a left-handed Z-helical structure in vivo when topA is switched off. In S. cerevisiae, the topological state of intracellular DNA is little affected by the cellular levels of the topoisomerases.

Base sequence effects in double helical DNA. I. Potential energy estimates of local base morphology

A series of potential energy calculations have been carried out to estimate base sequence dependent structural differences in B-DNA. Attention has been focused on the simplest dimeric fragments that can be used to build long chains, computing the energy as a function of the orientation and displacement of the 16 possible base pair combinations within the double helix. Calculations have been performed, for simplicity, on free base pairs rather than complete nucleotide units. Conformational preferences and relative flexibilities are reported for various combinations of the roll, tilt, twist, lateral displacement, and propeller twist of individual residues. The predictions are compared with relevant experimental measures of conformation and flexibility, where available. The energy surfaces are found to fit into two distinct categories, some dimer duplexes preferring to bend in a symmetric fashion and others in a skewed manner. The effects of common chemical substitutions (uracil for thymine, 5-methyl cytosine for cytosine, and hypoxanthine for guanine) on the preferred arrangements of neighboring residues are also examined, and the interactions of the sugar-phosphate backbone are included in selected cases. As a first approximation, long range interactions between more distant neighbors, which may affect the local chain configuration, are ignored. A rotational isomeric state scheme is developed to describe the average configurations of individual dimers and is used to develop a static picture of overall double helical structure. The ability of the energetic scheme to account for documented examples of intrinsic B-DNA curvature is presented, and some new predictions of sequence directed chain bending are offered.

Structure of the alpha-form of poly[d(A)].poly[d(T)] and related polynucleotide duplexes

The alpha-form of poly[d(A)].poly[d(T)], observed in fibers at high (greater than 80%) relative humidity, is a 10-fold double-helical structure of pitch 3.2 nm. This new X-ray analysis shows that the two strands of the double helix are of the same kind conformationally and both B-like in containing C-2'-endo-puckered deoxyribose rings. Nevertheless, the two strands are different enough for the overall morphology of the duplex to resemble that of the heteromerous model for the drier (beta) form of poly[d(A)].poly[d(T)] in which one strand has C-2'-endo rings and the other C-3'-endo. Since the orientations of the bases in poly[d(A)].poly[d(T)] are persistently different from those of classical B-DNA it is likely that there will be local bending (about 10 degrees) at the junctions between general sequence tracts and the oligo[d(A)].oligo[d(T)] tracts that occur in some native DNAs. The conclusions about the structure of alpha-poly[d(A)].poly[d(T)] are reinforced by independent analyses of similar X-ray diffraction patterns from poly[d(A)].poly[d(U)] and poly[d(A-I)].poly[d(C-T)].

Molecular mechanics calculations of dA12.dT12 and of the curved molecule d(GCTCGAAAAA)4.d(TTTTTCGAGC)4

Using the AMBER software package (Weiner and Kollman 1981) substantially modified for electrostatic contributions, the structural energies of the double-stranded oligonucleotides dA12.dT12 and d(GCTCGAAAAA)4.d(TTTTTCGAGC)4 were minimized. Using various starting structures for the molecule dA12.dT12, one final structure is obtained which possesses the experimentally determined properties of poly(dA).poly(dT). This structure is an A-form-B-form-hybrid structure similar to that of Arnott et al. (1983). The dA-strand is similar to an A-form while the dT-strand is similar to normal B-form. This structure and separately optimized B-form sequence stretches were used to construct the double-stranded fragment d(GCTCGAAAAA)4 which again was optimized. This sequence, when imbedded in a DNA fragment as contiguous repeats, shows a gel migration anomaly which has been interpreted as stable curvature of the DNA (Diekmann 1986). The calculated structure of this sequence indeed has a curved helix axis and is discussed as a model for curved DNA. A theoretical formalism is presented which allows one to calculate the structural parameters of any nucleic acid double helix in two different geometrical representations. This formalism is used to determine the parameters of the base-pair orientations of the curved structure in terms of wedge as well as cylindrical parameters. In the structural model presented here, the curvature of the helix axis results from an alternation of two different DNA structures in which the base-pairs possess different angles with the helix axis ('cylinder tilt'). Resulting from geometric restraints, a negative cylinder tilt angle correlates strongly with the closing of the minor groove ('wedge roll'). The blocks with different structure are not exactly coincident with the dA5-blocks and the B-DNA stretches. Within the dA5 block, base-pair tilt and wedge roll adopt large values which proceed into the 3' flanking B-DNA sequence by about one base-pair. These properties of the structure calculated here are discussed in terms of different models explaining DNA curvature.

The structure of poly(dA):poly(dT) in a condensed state and in solution

New X-ray and energetically optimal models of poly(dA):poly(dT) with the hydration spine in the minor groove have been compared with the NMR data in solution (Behling, R.W. and Kearns, D.R. (1986) Biochemistry 25, 3335-3346). These models have been refined to achieve a better fit with the NMR data. The obtained results suggest that the poly(dA):poly(dT) structure in a condensed state is similar to that in solution. The proposed conformations of poly(dA):poly(dT), unlike the classic B form, satisfy virtually all geometrical requirements which follow from the NMR data. Thus, the X-ray and energetically optimal poly(dA):poly(dT) structures (or those with slight modifications) can be considered as credible models of the poly(dA):poly(dT) double helix in solution. One of the features distinguishing these models from the classic B form is a narrowed minor groove.

Helical repeat of DNA in solution. The V curve method

The V-like dependence of the electrophoretic mobility of small DNA rings on their topological constraint, which has been documented in a recent paper [Zivanovic et al. (1986), J. Mol. Biol., 192, 645-660], has been explored as a tool to measure the helical twist of the torsionally unstressed duplex. The method was applied to single mixed sequence fragments approximately 350 to 1400 base pairs in length, providing estimates of their average helical periodicity. It was also used to compare two DNA fragments, and to evaluate the helical repeat of poly(dA.dT).poly(dA.dT) and poly(dA).poly(dT) inserts, and the helix unwindings associated with dA and dC methylations by dam and Hhal methylases, respectively. Data were found to be highly reproducible and helical repeat estimates were in good agreement with those obtained from previous techniques.

Anomalous structure and properties of poly (dA).poly(dT). Computer simulation of the polynucleotide structure with the spine of hydration in the minor groove

The results of the search for low-energy conformations of poly(dA).poly(dT) and of the poly(dA).poly(dT) "complex" with the spine of hydration similar to that found by Dickerson and co-workers (Kopka, M.L., Fratini, A.V., Drew, H.R. and Dickerson, R.E. (1983) J. Mol. Biol. 163, 129-146) in the minor groove of the CGCGAATTCGCG crystals are described. It is shown that the existence of such a spine in the minor groove of poly(dA).poly(dT) is energetically favourable. Moreover, the spine of hydration makes the polynucleotide conformation similar to the poly(dA).poly(dT) structure in fibers and to the conformation of the central part of CGCGAATTCGCG in crystals; it also acquires features characteristic of the structure of poly(dA).poly(dT) and DNA oligo(dA)-tracts in solution. It is shown that the existence of the TpA step in conformations characteristic of the poly(dA).poly(dT) complex with the spine of hydration is energetically unfavourable (in contrast to the ApT step) and therefore this step should result in destabilization of the spine of hydration in the DNA minor groove. Thus, it appears that the spine of hydration as described by Dickerson and co-workers is unlikely to exist in the poly d(A-T).poly d(A-T) structure. The data obtained permit us to interpret a large body of experimental facts concerning the unusual structure and properties of poly(dA).poly(dT) and oligo(dA)-tracts in DNA both in fibers and in solution. The results provide evidence of the existence of the minor groove spine of hydration both in fibers and in solution on A/T tracts of DNA which do not contain the TpA step. The spine plays an active role in the formation of the anomalous conformation of these tracts.

(A-T)n tracts embedded in random sequence DNA--formation of a structure which is chemically reactive and torsionally deformable

Alternating d(A-T)n sequences which are contiguous with DNA of effectively random sequence have an abnormal conformation in linear DNA molecules. These regions are strongly reactive towards chemical modification by osmium tetroxide, and are preferentially cleaved by micrococcal nuclease. Both the chemical modification and the enzymic cutting occur uniformly through the alternating tract, and there is no evidence for enzyme or chemical sensitivity in the interfaces between the tract and DNA of normal conformation. These reactivities have a requirement for an alternating sequence. In addition to chemical reactivity, alternating (A-T)n sequences exhibit anomalously small twist changes on cruciform formation, suggesting that the pre-extruded DNA is underwound. We propose that the alternating sequences adopt an altered conformation which is subject to easy torsional deformation.

Properties of supercoiled DNA in gel electrophoresis. The V-like dependence of mobility on topological constraint. DNA-matrix interactions

The dependence of the electrophoretic mobility of small DNA rings on topological constraint was investigated in acrylamide or agarose gels as a function of DNA size (from approximately 350 to 1400 base-pairs), gel concentration and nucleotide sequence. Under appropriate adjustment between the size of the DNA and the gel concentration, this dependence was found to be V-shaped in a limited interval around constraint O, the minimum mobility at the apex of the V being obtained for relaxed DNA. Analysis of the DNA size dependence of the V suggests that it is the result of a modulated compaction of the DNA rings by the gel matrix. Compaction appears to be maximum upon relaxation, and to decrease with increase in supercoiling. Consistent with this interpretation, gels were found to oppose structural departures from the B helix, such as Z transition and cruciform extrusion, which tend to relax the DNA molecule and make it more expanded. In contrast, when DNA size or gel concentration are large enough relative to one another, U shapes are observed instead of Vs, as a consequence of an increase in the mobility of the rings closer to relaxation. The relevance of these results to the situation of superhelical DNA in vivo is discussed. Application of the V to the measurement of the DNA helical twist is mentioned.

Hydrophobic interactions in the major groove can influence DNA local structure

We have considered hydrophobic interactions among aliphatic hydrocarbon groups in A/T sequences. The slightly overwound sequences (T)n.(A)n yield structures with tightly stacked methyl groups along one side of the major groove. The sequence TTAA may yield a sharp bend by folding together the two pairs of stacked methyls on the opposite sides of the major groove. Thus the sequence can affect the formation of either a smooth bend or a sharp kink. These sequence dependent local conformations may be related to a number of biological results.

Nuclease recognition of an alternating structure in a d(AT)14 plasmid insert

The nuclease reactivity and specificity of a cloned tract of poly X (dA-dT) X poly(dA-dT) has been explored. Digestion with DNAse I, Mung Bean nuclease, S1 nuclease, DNAse II, and copper (1,10-phenanthroline)2 on a 256 base pair restriction fragment containing d(AT)14A revealed a dinucleotide repeat structure for the alternating sequence. Furthermore, conditions which wind or unwind the linear DNA had little effect on the reactivity of the AT insert. These preferred cleavages offer insights to structural alterations within the DNA helix which differ from A, B, or Z-DNA. Nucleation into flanking sequences by this structural alteration was not observed.

DNA bending at adenine . thymine tracts

Intrinsic bending of DNA molecules results from local structural polymorphism in regions of homopolymeric dA . dT which are at least 4 base pairs long; the A . T tracts must be repeated in phase with the helix screw. Bending, in the direction of base-pair tilt rather than roll, occurs at the junctions between the A . T tract and adjacent B-DNA, with a larger angle at the 3' than at the 5' end of the A tract.

Influence of amino acid and peptide counterions on the conformation of DNA

We describe fibre diffraction studies on the interaction of DNA with different amino acids and peptides. The B form of DNA, with ten base-pairs per turn, is always found at high levels of humidity. We suggest that this pitch is observed because the DNA molecules are maintained in a straight position. In solution, the DNA molecules are bent and may have a larger pitch. The A form of DNA is never found upon dehydration. Instead, the B form may be either stabilized by the counterions or altered so that the number of base-pairs per helical turn decreases upon dehydration. Alteration is favoured either by small counterions that have a single charge or by large basic polypeptides and proteins. Stabilization is favoured by small counterions that have several charged groups. A third type of behaviour is found with some amino acids that contain hydrophobic groups, which destabilize the secondary structure of DNA, probably due to a modification of its intramolecular interactions. We have not detected any specific effect of amino acid side-chains, although the amino acid sequence has a clear influence on the interaction. We think that these observations are of interest in the pursuit of more detailed crystallographic studies on protein-DNA interactions.

Curved DNA: design, synthesis, and circularization

Curved DNA molecules and unusually small circles have been obtained by ligation of synthetic 21-base precursors: (sequence in text). The ligation resulted in the formation of double-stranded oligo-(precursor)s possessing a strong 10.5-base-pair (bp) periodicity of the runs of adenines. Two-dimensional polyacrylamide gel electrophoresis of the ligation products showed two distinct families of spots: (i) noncircular oligo(precursor)s of 21 to 231 bp (1- to 11-mers) and (ii) four circles from 105 to 168 bp (eluted and analyzed by denaturing gel electrophoresis). The noncircular oligomers exhibited anomalously slow migration, as if they were as much as three times longer than they actually are. The amount of circular products peaked sharply at approximately equal to 126 bp, near which size the circles have been estimated to be nonconstrained both torsionally and in terms of bending. The nonconstrained circularization provides a technique for the direct measurement of the inherent curvature of DNA in solution. From the size of the circles, an estimate of 8.7 degrees is obtained for the absolute value of the AA X TT wedge angle (roll and tilt combined).

On the sequence determinants and flexibility of the kinetoplast DNA fragment with abnormal gel electrophoretic mobilities

A 410 base-pair (bp) Sau3A restriction fragment derived from a Leishmania tarentolae kinetoplast DNA minicircle, which is known to have slower than expected electrophoretic mobilities in polyacrylamide gels, has been cloned in a plasmid and deletions from one end of the cloned segment have been constructed. Analysis of the gel electrophoretic mobility data of a large number of restriction fragments derived from the kinetoplast DNA clone and its deletion subclones has led to the conclusion that two sequences, one in the region bp 100 to 170 and the other bp 190 to 250, both numbered from one end of the 410 bp kinetoplast DNA segment, are important for the abnormal gel electrophoretic behavior of the kinetoplast DNA fragment. One common feature of these sequences is the periodic presence of short runs of A residues (3 to 6 As in each); auto-correlation analysis of these runs of A residues shows a strong harmonic component with a period around 11 bp. These results support and extend the previous analysis of Wu & Crothers (1984). The abnormal electrophoretic behavior is accentuated at low temperature and by the addition of Mg2+ to the electrophoresis buffer; addition of Na+ has the opposite effect. Insertion of sequences derived from the kinetoplast DNA fragment into nicked circular DNA causes no unexpected change in its electrophoretic mobility in agarose gel, suggesting that the 410 bp sequence, or segments of it, has no significant spatial writhe. Abnormal shifts in agarose gel mobilities are observed, however, when certain segments of the kinetoplast DNA are inserted into positively or negatively supercoiled DNA topoisomers. These results are consistent with a bent structure of the kinetoplast DNA in which the bend has zero writhe in its undistorted form but is easily distorted.

Facile cruciform formation by an (A-T)34 sequence from a Xenopus globin gene

We have studied the structure adopted by an (A-T)34 sequence from a Xenopus globin gene when present in a negatively supercoiled plasmid. A variety of enzyme and chemical probing experiments and electrophoretic migration shift methods reveal that the sequence adopts cruciform geometry at moderate levels of supercoiling. The structure has the lowest free energy of formation yet observed for a cruciform, and no detectable kinetic barrier preventing rapid interconversion between extruded and unextruded conformations. Analysis of band-shift experiments reveals a twist change on cruciform formation of -5.8, slightly smaller than the -6.5 we would predict on the basis of a transition from B DNA. An attractive explanation consistent with this discrepancy is that the (A-T)34 stretch is locally underwound to about 11.7 base-pairs/helical turn at low levels of supercoiling. This calculation is made on the assumption that the cruciform junction is structurally similar to those examined previously, which is supported by the nuclease digestion results. This perturbed helical structure could be of considerable biological significance.