Influence of dangling thymidine residues on the stability and structure of two DNA duplexes

We have employed temperature-dependent UV spectroscopy, circular dichroism (CD), 400-MHz proton nuclear magnetic resonance (NMR), and computer modeling to characterize both structurally and thermodynamically the influence of unpaired, dangling thymidine residues (T) on the thermal stability and melting behavior of two DNA core duplexes. The specific DNA double helices that we have investigated in this work are core duplexes [d(GC)3]2 (I) and [d(CG)3]2 (IV), 3' dangling T derivatives [d[(GC)3TT]]2 (II) and [d[(CG)3TT]]2 (V), and 5' dangling T derivatives [d[TT(GC)3]]2 (III) and [d[TT(CG)3]]2 (VI). Our experimental data allow us to reach the following conclusions: (1) For both core duplexes (I and IV), the addition of dangling T residues on either the 5' or 3' end causes an increase in the optical melting temperature tm. (2) For both core duplexes, 5' dangling T residues induce a greater increase in the optical tm's than 3' dangling T residues. (3) For both cores duplexes, the increase in tm induced by the addition of dangling T residues is enthalpic in origin, with 5' dangling T residues inducing a greater increase in the van't Hoff transition enthalpy than 3' dangling T's. (4) Dangling T residues cause downfield shifts in all of the nonexchangeable aromatic protons of the [d(GC)3]2 core duplex (I), with the 5' T residues inducing the largest shifts. For the most part, this trend does not hold with the [d(CG)3]2 core duplex (IV). (5) For both core duplexes, the addition of dangling T residues causes an increase in the NMR tm's of almost all the nonexchangeable aromatic protons of the core duplex.(ABSTRACT TRUNCATED AT 250 WORDS)     

Parallel-stranded DNA with mixed AT/GC composition: role of trans G.C base pairs in sequence dependent helical stability

Parallel-stranded (ps) DNAs with mixed AT/GC content comprising G.C pairs in a varying sequence context have been investigated. Oligonucleotides were devised consisting of two 10-nt strands complementary either in a parallel or in an antiparallel orientation and joined via nonnucleotide linkers so as to form 10-bp ps or aps hairpins. A predominance of intramolecular hairpins over intermolecular duplexes was achieved by choice of experimental conditions and verified by fluorescence determinations yielding estimations of rotational relaxation times and fractional base pairing. A multistate mode of ps hairpin melting was revealed by temperature gradient gel electrophoresis (TGGE). The thermal stability of the ps hairpins with mixed AT/GC content depends strongly on the specific sequence in a manner peculiar to the ps double helix. The thermodynamic effects of incorporating trans G.C base pairs into an AT sequence are context-dependent: an isolated G. C base pair destabilizes the duplex whereas a block of > or =2 consecutive G.C base pairs exerts a stabilizing effect. A multistate heterogeneous zipper model for the thermal denaturation of the hairpins was derived and used in a global minimization procedure to compute the thermodynamic parameters of the ps hairpins from experimental melting data. In 0.1 M LiCl at 3 degrees C, the formation of a trans G.C pair in a GG/CC sequence context is approximately 3 kJ mol(-)(1) more favorable than the formation of a trans A.T pair in an AT/TA sequence context. However, GC/AT contacts contribute a substantial unfavorable free energy difference of approximately 2 kJ mol(-)(1). As a consequence, the base composition and fractional distribution of isolated and clustered G.C base pairs determine the overall stability of ps-DNA with mixed AT/GC sequences. Thus, the stability of ps-DNA comprising successive > or =2 G.C base pairs is greater than that of ps-DNA with an alternating AT sequence, whereas increasing the number of AT/GC contacts by isolating G.C base pairs exerts a destabilizing effect on the ps duplex. Molecular modeling of the various helices by force field techniques provides insight into the structural basis for these distinctions.     

Thermodynamic analysis of monoclonal antibody binding to duplex DNA.

A technique based on fluorescence polarization (anisotropy) was used to measure the binding of antibodies to DNA under a variety of conditions. Fluorescein-labeled duplexes of 20 bp in length were employed as the standard because they are stable even at low ionic strength yet sufficiently short so that both arms of an IgG cannot bind to the same duplex. IgG Jel 274 binds duplexes in preference to single-stranded DNA; in 80 mM NaCl Kobs for (dG)20.(dC)20 is 4.1x10(7) M-1 compared with 6.4x10(5) M-1 for d(A5C10A5). There is little sequence specificity, but the interaction is very dependent on ionic strength. From plots of log Kobs against log[Na+] it was deduced that five or six ion pairs are involved in complex formation. At low ionic strength,Kobs is independent of temperature and complex formation is entropy driven with DeltaH degrees obs and DeltaC degrees p,obs both zero. In contrast, in 80 mM NaCl DeltaC degrees p,obs is -630 and -580 cal mol-1K-1 for [d(TG)]10.[d(CA)]10 and (dG)20.(dC)20 respectively. IgG Jel 241 also binds more tightly to duplexes than single-stranded DNA, but sequence preferences were apparent. The values for Kobs to [d(AT)]20 and [d(GC)]20 are 2.7x10(8) and 1.3x10(8) M-1 respectively compared with 5.7x10(6) M-1 for both (dA)20. (dT)20 and (dG)20.(dC)20. As with Jel 274, the binding of Jel 241 is very dependent on ionic strength and four or five ionic bonds are involved in complex formation with all the duplex DNAs which were tested. DeltaC degrees p,obs for Jel 241 binding to [d(AT)]20 was negative (-87 cal mol-1K-1) in 80 mM NaCl but was zero at high ionic strength (130 mM NaCl). Therefore, for duplex-specific DNA binding antibodies DeltaC degrees p,obs is dependent on [Na+] and a large negative value does not correlate with sequence-specific interactions.     

1H NMR determination of base-pair lifetimes in oligonucleotides containing single base mismatches

Proton nuclear magnetic resonance (NMR) spectroscopy is employed to characterize the kinetics of base-pair opening in a series of 9mer duplexes containing different single base mismatches. The imino protons from the different mismatched, as well as fully matched, duplexes are assigned from the imino-imino region in the WATERGATE NOESY spectra. The exchange kinetics of the imino protons are measured from selective longitudinal relaxation times. In the limit of infinite exchange catalyst concentration, the exchange times of the mismatch imino protons extrapolate to much shorter lifetimes than are commonly observed for an isolated GC base pair. Different mismatches exhibit different orders of base-pair lifetimes, e.g. a TT mismatch has a shorter base-pair lifetime than a GG mismatch. The effect of the mismatch was observed up to a distance of two neighboring base pairs. This indicates that disruption in the duplex caused by the mismatch is quite localized. The overall order of base-pair lifetimes in the selected sequence context of the base pair is GC > GG > AA > CC > AT > TT. Interestingly, the fully matched AT base pair has a shorter base-pair lifetime relative to many of the mismatches. Thus, in any given base pair, the exchange lifetime can exhibit a strong dependence on sequence context. These findings may be relevant to the way mismatch recognition is accomplished by proteins and small molecules.     

A nonlinear dynamic model of DNA with a sequence-dependent stacking term

No simple model exists that accurately describes the melting behavior and breathing dynamics of double-stranded DNA as a function of nucleotide sequence. This is especially true for homogenous and periodic DNA sequences, which exhibit large deviations in melting temperature from predictions made by additive thermodynamic contributions. Currently, no method exists for analysis of the DNA breathing dynamics of repeats and of highly G/C- or A/T-rich regions, even though such sequences are widespread in vertebrate genomes. Here, we extend the nonlinear Peyrard–Bishop–Dauxois (PBD) model of DNA to include a sequence-dependent stacking term, resulting in a model that can accurately describe the melting behavior of homogenous and periodic sequences. We collect melting data for several DNA oligos, and apply Monte Carlo simulations to establish force constants for the 10 dinucleotide steps (CG, CA, GC, AT, AG, AA, AC, TA, GG, TC). The experiments and numerical simulations confirm that the GG/CC dinucleotide stacking is remarkably unstable, compared with the stacking in GC/CG and CG/GC dinucleotide steps. The extended PBD model will facilitate thermodynamic and dynamic simulations of important genomic regions such as CpG islands and disease-related repeats.     

Sequence dependence of DNA conformational flexibility

By using conformational free energy calculations, we have studied the sequence dependence of flexibility and its anisotropy along various conformational variables of DNA base pairs. The results show the AT base step to be very flexible along the twist coordinate. On the other hand, homonucleotide steps, GG(CC) and AA(TT), are among the most rigid sequences. For the roll motion that would correspond to a bend, the TA step is most flexible, while the GG(CC) step is least flexible. The flexibility of roll is quite anisotropic; the ratio of fluctuations toward the major and minor grooves is the largest for the GC step and the smallest for the AA(TT) and CG steps. Propeller twisting of base pairs is quite flexible, especially of A.T base pairs; propeller twist can reach 19 degrees by thermal fluctuation. We discuss the effect of electrostatic parameters, comparison with available experimental results, and biological relevance of these results.     

Energetics of sequence-specific protein-DNA association: conformational stability of the DNA binding domain of integrase Tn916 and its cognate DNA duplex

Sequence-specific DNA recognition by bacterial integrase Tn916 involves structural rearrangements of both the protein and the DNA duplex. Energetic contributions from changes of conformation, thermal motions and soft vibrational modi of the protein, the DNA, and the complex significantly influence the energetic profile of protein-DNA association. Understanding the energetics of such a complicated system requires not only a detailed calorimetric investigation of the association reaction but also of the components in isolation. Here we report on the conformational stability of the integrase Tn916 DNA binding domain and its cognate 13 base pair target DNA duplex. Using a combination of temperature and denaturant induced unfolding experiments, we find that the 74-residue DNA binding domain is compact and unfolds cooperatively with only small deviation from two-state behavior. Scanning calorimetry reveals an increase of the heat capacity of the native protein attributable to increased thermal fluctuations. From the combined calorimetric and spectroscopic experiments, the parameters of protein unfolding are T(m) = 43.8 +/- 0.3 degrees C, DeltaH(m) = 255 +/- 18 kJ mol(-1), DeltaS(m) = 0.80 +/- 0.06 kJ mol(-1), and DeltaC(p) = 5.0 +/- 0.8 kJ K(-1) mol(-1). The DNA target duplex displays a thermodynamic signature typical of short oligonucleotide duplexes: significant heat absorption due to end fraying and twisting precedes cooperative unfolding and dissociation. The parameters for DNA unfolding and dissociation are DeltaH(m) = 335 +/- 4 kJ mol(-1) and DeltaC(p) = 2.7 +/- 0.9 kJ K(-(1) mol(-1). The results reported here have been instrumental in interpreting the thermodynamic features of the association reaction of the integrase with its 13 base pair target DNA duplex reported in the accompanying paper [Milev et al. (2003) Biochemistry 42, 3481-3491].     

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.     

The Intramolecular Impact to the Sequence Specificity of B→A Transition: Low Energy Conformational Variations in AA/TT and GG/CC Steps

Abstract It is well known, that local B→A transformation in DNA is involved in several biological processes. In vitro B?A transition is sequence-specific. The physical basis of this specificity is not known yet. Here we analyze the effect of intramolecular interactions on the structural behavior of the GG/CC and AA/TT steps. These steps exemplify sequence specific bias to the B- or A-form structure. Optimization of potential energy of the molecular systems composed of an octanucle-otide, neutralized by Na(+) and solvated with TIP3P water molecules in rectangular box with periodic boundary conditions gives the statistically representative sets of low energy structures for GG/CC and AA/TT steps in the middle of the diverse flanking sequences. Permissible 3D variations of GG/CC and AA/TT, and correlation of the relative motion of base pairs in these steps were analyzed. AA/TT step permits high variability for low energy conformers in the B-form DNA and small variability for low energy conformers in the A-form DNA. In contrast GG/CC step permits high variability for low energy conformers in the A-form DNA and small variability for low energy conformers in the B-form DNA. The relative motion of base pairs in GG/CC step is high correlated, while in AA/TT step this correlation is notably less. Atom-atom interactions inside-the-step always favors the B-form and their component - stacking interactions (atomatom interactions between nucleic bases) is crucial for the duplex stabilization. Formation of the A-form for both steps is a result of interactions with the flanking sequences and water-cation environment in the box. The average energy difference between conformations presenting B-form and A-form for the GG/CC step is high, while for the AA/TT step it is rather low. Thus, intramolecular interactions in GG/CC and AA/TT steps affect the possible conformational diversity ("conformational entropy") of the A- and B- type structures of DNA step. This determines the known bias of the A-form DNA depending on the enrichment of sequences with GG/CC. If structural tuning during the process of protein-DNA complex formation lead to the local B→A transformation of DNA, it is largely directed by high conformational diversity of GG/CC step in the A-form. In such a case the presence in the target site of both kinds of examined steps ensures the reversible character of ligand binding.     

Thermodynamic analysis of stacking hybridization of oligonucleotides with DNA template.

Contiguous stacking hybridization of oligodeoxyribonucleotides with DNA as template was investigated using three types of complexes: oligonucleotide contiguously stacked with the stem of the preformed minihairpin (complexes I), oligonucleotide tandems containing two (complexes II) or three (complexes III) short oligomers with a common DNA template. Enthalpy Delta H degrees and entropy Delta S degrees of the coaxial stacking of adjacent duplexes were determined for GC/G*pC, GT/A*pC, AC/G*pT, AT/A*pT, CT/A*pG, AG/C*pT, AA/T*pT and TT/A*pA nicked (*) dinucleotide base pairs. The maximal efficiency of co-operative interaction was found for the GC/G*pC interface (Delta G degrees(NN/N*pN)=-2.7 kcal/mol) and the minimal one for the AA/T*pT interface (Delta G degrees(NN/N*pN)=-1.2 kcal/mol) at 37 degrees C. As a whole, the efficiency of the base pairs interaction Delta G degrees(NN/N*pN) in the nick is not lower than that within the intact DNA helix (Delta G degrees(NN/NN)).These observed Delta G degrees(NN/N*pN) values are proposed may include the effect of the partial removal of fraying at the adjacent helix ends additionally to the effect of the direct stacking of the terminal base pairs in the duplex junction (Delta G degrees(NN/NN). The thermodynamic parameters have been found to describe adequately the formation of all tandem complexes of the II and III types with oligonucleotides of various length and hybridization properties. The performed thermodynamic analysis reveals features of stacking oligonucleotide hybridization which allow one to predict the temperature dependence of association of oligonucleotides and the DNA template within tandem complexes as well as to determine optimal concentration for formation of these complexes characterized by high co-operativity level.     

Thermodynamics of RNA duplexes with tandem mismatches containing a uracil-uracil pair flanked by C.G/G.C or G.C/A.U closing base pairs

The thermodynamics governing the denaturation of RNA duplexes containing 8 bp and a central tandem mismatch or 10 bp were evaluated using UV absorbance melting curves. Each of the eight tandem mismatches that were examined had one U-U pair adjacent to another noncanonical base pair. They were examined in two different RNA duplex environments, one with the tandem mismatch closed by G.C base pairs and the other with G.C and A.U closing base pairs. The free energy increments (Delta Gdegrees(loop)) of the 2 x 2 loops were positive, and showed relatively small differences between the two closing base pair environments. Assuming temperature-independent enthalpy changes for the transitions, (Delta Gdegrees(loop)) for the 2 x 2 loops varied from 0.9 to 1.9 kcal/mol in 1 M Na(+) at 37 degrees C. Most values were within 0.8 kcal/mol of previously estimated values; however, a few sequences differed by 1.2-2.0 kcal/mol. Single strands employed to form the RNA duplexes exhibited small noncooperative absorbance increases with temperature or transitions indicative of partial self-complementary duplexes. One strand formed a partial self-complementary duplex that was more stable than the tandem mismatch duplexes it formed. Transitions of the RNA duplexes were analyzed using equations that included the coupled equilibrium of self-complementary duplex and non-self-complementary duplex denaturation. The average heat capacity change (DeltaC(p)) associated with the transitions of two RNA duplexes was estimated by plotting DeltaH degrees and DeltaS degrees evaluated at different strand concentrations as a function of T(m) and ln T(m), respectively. The average DeltaC(p) was 70 +/- 5 cal K(-)(1) (mol of base pairs)(-)(1). Consideration of this heat capacity change reduced the free energy of formation at 37 degrees C of the 10 bp control RNA duplexes by 0.3-0.6 kcal/mol, which may increase Delta Gdegrees(loop) values by similar amounts.     

Influence of an exocyclic guanine adduct on the thermal stability, conformation, and melting thermodynamics of a DNA duplex.

As part of an overall program to characterize the impact of mutagenic lesions on the physiochemical properties of DNA, we report here the results of a comparative spectroscopic study on pairs of DNA duplexes both with and without an exocyclic guanine lesion. Specifically, we have studied a family of four 13-mer duplexes of the form d(CGCATGYGTACGC).d(GCGTACZCATGCG) in which Y is either the normal deoxyguanosine residue (G) or the exocyclic guanine adduct 1,N2-propanodeoxyguanosine (X), while Z is either deoxycytosine (C) or deoxyadenosine (A). Thus, the four duplexes studied, which can be designated by the identity of their central Y.Z base pair, are a Watson-Crick duplex (GC), a duplex with a central mismatch (GA), and two duplexes with exocyclic guanine lesions (X), that differ only by the base opposite the lesion (XC and XA). The data derived from our spectroscopic measurements on these four duplexes have allowed us to evaluate the influence of the exocyclic guanine lesion, as well as the base opposite the lesion, on the conformation, thermal stability, and melting energetics of the host DNA duplex. To be specific, our circular dichroism (CD) spectra show that the exocyclic guanine lesion induces alterations in the duplex structure, while our temperature-dependent optical measurements reveal that these lesion-induced structural alterations reduce the thermal stability, the transition enthalpy, and the transition free energy of the duplex.(ABSTRACT TRUNCATED AT 250 WORDS)     

Structure and stability of the consecutive stereoregulated chiral phosphorothioate DNA duplex

The duplex structures of the stereoregulated phosphorothioate DNAs, [R(p),R(p)]- and [S(p),S(p)]-[d(GC(ps)T(ps)ACG)] (ps, phosphorothioate; PS-DNA), with their complementary RNA have been investigated by combined use of (1)H NMR and restrained molecular dynamics calculation. Compared to those obtained for the unmodified duplex structures (PO-DNA.RNA), the NOE cross-peak intensities are virtually identical for the PS-DNA.RNA hybrid duplexes. The structural analysis on the basis of the NOE restraints reveals that all of the three DNA.RNA duplexes take a A-form conformation and that there is no significant difference in the base stacking for the DNA.RNA hybrid duplexes. On the other hand, the NOE cross-peak intensities of the protons around the central T(ps)A step of the PS-DNA.DNA duplexes are apparently different from those of PO-DNA. DNA. The chemical shifts of H8/6 and H1' at the T(ps)A step are also largely different among PS-DNA.DNAs and PO-DNA.DNA, suggesting that the DNA.DNA structure is readily changed by the introduction of the phosphorothioate groups to the central T(p)A step. The structure calculations indicate that all of these DNA.DNA duplexes are B-form although there exist some small differences in helical parameters between the [R(p),R(p)]- and [S(p),S(p)]PS-DNA.DNA duplexes. The melting temperatures (T(m)) were determined for all of the duplexes by plotting the chemical shift change of isolated peaks as a function of temperature. For the PS-DNA.RNA hybrid duplexes, the [S(p),S(p)] isomer is less stable than the [R(p),R(p)] isomer while this trend is reversed for the PS-DNA.DNA duplexes. Consequently, although the PS-DNA.RNA duplexes take the similar A-form structure, the duplex stability is different between PS-DNA.RNA duplexes. The stability of the DNA.RNA duplexes may not be governed by the A-form structure itself but by some other factors such as the hydration around the phosphorothioate backbone, although the T(m) difference of the DNA.DNA duplexes could be explained by the structural factor.     

Sequence specificity of DNA cleavage by bis(1,10-phenanthroline)copper(I): effects of single base pair transitions on the cleavage of preferred pyrimidine-purine-pyrimidine triplets

The cleavage of DNA restriction fragments by bis(1,10-phenanthroline)copper(I) [[(OP)2CuI]+] is sequence dependent: the trimer TAT is most strongly preferred, while the trimer TGT and tetramers TAAT, TAGT, and CAGT are strongly to moderately preferred [Veal, J. M., & R. L. (1988) Biochemistry 27, 1822-1827]. [(OP)2CuI]+ cleavage of a series of oligonucleotide duplexes of the type 5'-CCCTPyPuPyCCCC-3'/3'-GGGAPuPyPuGGGG-5' (Py = pyrimidine; Pu = purine) was examined to determine the effects of purine substituents in the central triplet on specificity. The relative cleavage rates of different PyPuPy triplets in oligomers were similar to those observed for restriction fragments. The undecamer duplex containing the trimer TAT (TTATC) was most preferentially cleaved, predominantly at the central adenosine and the adjacent 3'-thymidine. Duplexes differing from TTATC by a single A.T----G.C transition in the central triplet were cleaved at significantly reduced rates relative to TTATC, the order of preference being TAT greater than TGT greater than TAC greater than CAT. By contrast, duplexes differing from TTATC by a single A.T----I.C transition were cleaved at rates similar to those for TTATC when the transition occurred at the 5'-pyrimidine or central purine [i.e., C(.I)AT and TIT]. A duplex containing the trimer TAC(.I) was cleaved at a reduced rate similar to the duplex containing TAC(.G). The guanine 2-amino group at positions 1 and 2, but not position 3, of a 5'-PyPuPy-3' trimer is therefore implicated as a strong inhibitor of DNA binding by the copper-phenanthroline complex.(ABSTRACT TRUNCATED AT 250 WORDS)     

Stabilities of consecutive A.C, C.C, G.G, U.C, and U.U mismatches in RNA internal loops: Evidence for stable hydrogen-bonded U.U and C.C.+ pairs.

The stability and structure of RNA duplexes with consecutive A.C, C.A, C.C, G.G, U.C, C.U, and U.U mismatches were studied by UV melting, CD, and NMR. The results are compared to previous results for GA and AA internal loops [SantaLucia, J., Kierzek, R., & Turner, D. H. (1990) Biochemistry 29, 8813-8819; Peritz, A., Kierzek, R., & Turner, D.H. (1991) Biochemistry 30, 6428-6436)]. The observed order for stability increments of internal loop formation at pH 7 is AG = GA approximately UU greater than GG greater than or equal to CA greater than or equal to AA = CU = UC greater than or equal to CC greater than or equal to AC. The results suggest two classes for internal loops with consecutive mismatches: (1) loops that stabilize duplexes and have strong hydrogen bonding and (2) loops that destabilize duplexes and may not have strong hydrogen bonding. Surprisingly, rCGCUUGCG forms a very stable duplex at pH 7 in 1 M NaCl with a TM of 44.8 degrees C at 1 x 10(-4) M and a delta G degrees 37 of -7.2 kcal/mol. NOE studies of the imino protons indicate hydrogen bonding within the U.U mismatches in a wobble-type structure. Resonances corresponding to the hydrogen-bonded uridines are located at 11.3 and 10.4 ppm. At neutral pH, rCGCCCGCG is one of the least stable duplexes with a TM of 33.2 degrees C and delta G degrees 37 of -5.1 kcal/mol. Upon lowering the pH to 5.5, however, the TM increases by 12 degrees C, and delta G degrees 37 becomes more favorable by 2.5 kcal/mol. The pH dependence of rCGCCCGCG may be due to protonation of the internal loop C's, since no changes in thermodynamic parameters are observed for rCGCUUGCG between pH 7 and 5.5. Furthermore, two broad imino proton resonances are observed at 10.85 and 10.05 ppm for rCGCCCGCG at pH 5.3, but not at pH 6.5. This is also consistent with C.C+ base pairs forming at pH 5.5. rCGCCAGCG and rGGCACGCC have a small pH dependence, with TM increases of 5 and 3 degrees C, respectively, upon lowering the pH from 7 to 5.5. rCGCCUGCG and rCGCUCGCG also show little pH dependence, with TM increases of 0.8 and 1.4 degrees C, respectively, upon lowering the pH to 5.5.(ABSTRACT TRUNCATED AT 400 WORDS)     

Thermodynamic, counterion, and hydration effects for the incorporation of locked nucleic acid nucleotides into DNA duplexes

A locked nucleic acid (LNA) monomer is a conformationally restricted nucleotide analogue with an extra 2'-O, 4'-C-methylene bridge added to the ribose ring. LNA-modified oligonucleotides are known to exhibit enhanced hybridization affinity toward complementary DNA and RNA. In this work, we have evaluated the hybridization thermodynamics of a series of LNA-substituted DNA octamers, modified to various extents by one to three LNA substitutions, introduced at either adenine (5'-AGCACCAG) or thymine (5'-TGCTCCTG) nucleotides. To understand the energetics, counterion effects, and the hydration contribution of the incorporation of LNA modification, a combination of spectroscopic and calorimetric techniques was used. The CD spectra of the corresponding duplexes showed that the modified duplexes adopt an A-type conformation. UV and DSC melting studies revealed that each type of duplex unfolds in a two-state transition. A complete thermodynamic profile at 5 degrees C indicated that the net effect of modification on thermodynamic parameters might be positional and that the neighboring bases flanking the modification might influence the favorable formation of the modified duplexes. Furthermore, relative to the formation of the unmodified reference duplexes, the formation of modified duplexes is accompanied by a higher uptake of counterions and a lower uptake of water molecules.     

Conformational properties of Z-forming DNA oligomers bearing terminal unpaired bases.

Three sets of semi-self-complementary deoxyribonucleotide decamers with the sequence XX-(5meCG)4, (5meCG)4-XX, or Y-(5meCG)4-Y, where XX = AA, CC, GG, or TT and Y = A, C, G, or T, were synthesized along with the self-complementary octamer (5meCG)4. The 8-mer duplex readily undergoes a B-to-Z conformational conversion upon increasing the NaCl concentration with a transitional midpoint of approximately 1.1 M NaCl. The 10-mers should form 8-bp duplexes a with core sequence of [(5meCG)4]2 with 5'-XX overhangs, 3'-XX overhangs, or 5',3'-Y/Y mismatches. Circular dichroism was employed to determine the conformations of all oligomers. Salt titrations were performed to measure the effect of overhangs and terminal mismatches on the B-to-Z conversion. In general, the presence of 5'-XX overhangs results in a transition midpoint equal to or slightly higher than the control, whereas the presence of 3'-XX overhangs results in a transition midpoint slightly lower than the control. The 3'-CC and 5'-GG overhangs are exceptions, with transition midpoints much higher than the control. These oligomers apparently form duplexes with 5',3'-C/C or 5',3'-G/G mismatches abutting a [(G5meC)4]2 duplex core. The presence of terminal mismatches in the third set of oligomers results in transition midpoints higher than the control. Ultraviolet absorbance methods were used to evaluate the effect of the various stacking motifs of the 10-mers on the thermodynamics of melting relative to the 8-mer for both B and Z conformations. We found that in both the B and Z conformations, the presence of an overhang stabilizes the [(5meCG)4]2 duplex, with the 5' overhangs having a greater stabilizing effect relative to the 3' overhangs. The presence of 5',3'-Y/Y mismatches also imparts a stabilizing effect on the control 8-mer in both the B and Z conformations. These results are discussed in terms of stacking interactions of the terminal unpaired bases.     

Sequence structure of human nucleosome DNA

Positional distributions of various dinucleotides in experimentally derived human nucleosome DNA sequences are analyzed. Nucleosome positioning in this species is found to depend largely on GG and CC dinucleotides periodically distributed along the nucleosome DNA sequence, with the period of 10.4 bases. The GG and CC dinucleotides oscillate counterphase, i.e., their respective preferred positions are shifted about a half-period from one another, as it was observed earlier for AA and TT dinucleotides. Other purine-purine and pyrimidine-pyrimidine dinucleotides (RR and YY) display the same periodical and counterphase pattern. The dominance of oscillating GG and CC dinucleotides in human nucleosomes and the contribution of AG(CT), GA(TC), and AA(TT) suggest a general nucleosome DNA sequence pattern - counterphase oscillation of RR and YY dinucleotides. AA and TT dinucleotides, commonly accepted as major players, are only weak contributors in the case of human nucleosomes.     

Relative specificities in binding of Watson-Crick base pairs by third strand residues in a DNA pyrimidine triplex motif.

The specificity of binding of Watson-Crick base pairs by third strand nucleic acid residues via triple helix formation was investigated in a DNA pyrimidine triplex motif by thermal melting experiments. The host duplex was of the type A10-X-A10: T10-Y-T10, and the third strand T10-Z-T10, giving rise to 16 possible triplexes with Z:XY inserts, 4 duplexes with the Watson-Crick base pairs (XY) and 12 duplexes with mismatch pairs (XZ), all of whose stabilities were compared. Two Z:XY combinations confirm the primary binding of AT and GC target pairs in homopurine.homopyrimidine sequences by T and C residues, respectively. All other Z:XY combinations in the T:AT environment result in triplex destabilization. While some related observations have been reported, the present experiments differ importantly in that they were performed in a T:AT nearest neighbor environment and at physiological ionic strength and pH, all of which were previously untested. The conclusions now drawn also differ substantially from those in previous studies. Thus, by evaluating the depression in Tm due to base triplet mismatches strictly in terms of third strand residue affinity and specificity for the target base pair, it is shown that none of the triplet combinations that destabilize qualify for inclusion in the third strand binding code for the pyrimidine triplex motif. Hence, none of the mismatch triplets afford a general way of circumventing the requirement for homopurine.homopyrimidine targets when third strands are predominated by pyrimidines, as others have suggested. At the same time, the applicability of third strand binding is emphasized by the finding that triplexes are equally or much more sensitive to base triplet mismatches than are Watson-Crick duplexes to base pair mismatches.