Base-Pair Opening in Double-Stranded Nucleic Acids

Investigation of changes in the arrangement of water molecules and salt ions surrounding different atoms of the DNA molecule during the melting process: a molecular dynamics simulation study

A molecular dynamics simulation was performed on a B-DNA duplex (CGCGAATTGCGC) at different temperatures. The DNA was immerged in a saltwater medium with 1 mol/L NaCl concentration. The arrangements of water molecules and cations around the different atoms of DNA on the melting pathway were investigated. Almost for all atoms of the DNA by double helix → single-stranded transition, the water molecules released from the DNA duplex and cations were close to single-stranded DNA, but this behavior was not clearly seen at melting temperatures. Therefore, release of water molecules and cations approaching the DNA by the increase of temperature does not have any effect on the sharpness of the transition curve. Most of the water molecules and cations were found to be around the negatively charged phosphate oxygen atoms. The number of water molecules released from the first shell hydration upon melting in the minor groove was higher than in the major groove, and intrusion of cations into the minor groove after melting was higher than into the major groove. The hydrations of imino protons were different from each other and were dependent on DNA bases.

Long lifetime of hydrogen-bonded DNA basepairs by force spectroscopy

Electron-tunneling data suggest that a noncovalently-bonded complex of three molecules, two recognition molecules that present hydrogen-bond donor and acceptor sites via a carboxamide group, and a DNA base, remains bound for seconds. This is surprising, given that imino-proton exchange rates show that basepairs in a DNA double helix open on millisecond timescales. The long lifetime of the three-molecule complex was confirmed using force spectroscopy, but measurements on DNA basepairs are required to establish a comparison with the proton-exchange data. Here, we report on a dynamic force spectroscopy study of complexes between the bases adenine and thymine (A-T, two-hydrogen bonds) and 2-aminoadenine and thymine (2AA-T, three-hydrogen bonds). Bases were tethered to an AFM probe and mica substrate via long, covalently linked polymer tethers. Data for bond-survival probability versus force and the rupture-force distributions were well fitted by the Bell model. The resulting lifetime of the complexes at zero pulling force was ~2 s for two-hydrogen bonds (A-T) and ~4 s for three-hydrogen bonds (2AA-T). Thus, DNA basepairs in an AFM pulling experiment remain bonded for long times, even without the stabilizing influence of base-stacking in a double helix. This result suggests that the pathways for opening, and perhaps the open states themselves, are very different in the AFM and proton-exchange measurements.

A pseudocatenane structure formed between DNA and A cyclic bisintercalator

Targeting double-stranded DNA with small molecules remains an active area of basic research. Herein is described a cyclic DNA bisintercalator that is based on two naphthalene diimide (NDI) intercalating units tethered by one linking element specific for binding in the minor groove and the other linking element specific for binding in the major groove. DNase I footprinting revealed a strong preference for binding the sequence 5'-GGTACC-3'. NMR structural studies of the complex with d(CGGTACCG)(2) verified a pseudocatenane structure in which the NDI units reside four base pairs apart, with one linker segment located in the minor groove and the other in the major groove consistent with the linker designs. To the best of our knowledge, this is the first structurally well-characterized pseudocatenane complex between a sequence specific cyclic bisintercalator and intact DNA.

Spectroscopic studies of position-specific DNA "breathing" fluctuations at replication forks and primer-template junctions

Junctions between ssDNA and dsDNA sequences are important in many cellular processes, including DNA replication, transcription, recombination, and repair. Significant transient conformational fluctuations ("DNA breathing") can occur at these ssDNA-dsDNA junctions. The involvement of such breathing in the mechanisms of macromolecular complexes that operate at these loci is not well understood, in part because these fluctuations have been difficult to measure in a position-specific manner. To address this issue we constructed forked or primer-template DNA constructs with 1 or 2 adjacent 2-aminopurine (2-AP) nucleotide residues (adenine analogues) placed at specific positions on both sides of the ssDNA-dsDNA junction. Unlike canonical DNA bases, 2-AP absorbs, fluoresces, and displays CD spectra at wavelengths >300 nm, where other nucleic acid and protein components are transparent. We used CD and fluorescence spectra and acrylamide quenching of these probes to monitor the extent and nature of DNA breathing of A-T base pairs at specific positions around the ssDNA-dsDNA junction. As expected, spectroscopically measurable unwinding penetrates approximately 2 bp into the duplex region of these junctions under physiological conditions for the constructs examined. Surprisingly, we found that 2-AP bases at ssDNA sites directly adjacent to ssDNA-dsDNA junctions are significantly more unstacked than those at more distant ssDNA positions. These local and transient DNA conformations on both sides of ssDNA-dsDNA junctions may serve as specific interaction targets for enzymes that manipulate DNA in the processes of gene expression.

DNA twisting flexibility and the formation of sharply looped protein-DNA complexes

Gene-regulatory complexes often require that pairs of DNA-bound proteins interact by looping-out short (often approximately 100-bp) stretches of DNA. The loops can vary in detailed length and sequence and, thus, in total helical twist, which radically alters their geometry. How this variability is accommodated structurally is not known. Here we show that the inherent twistability of 89- to 105-bp DNA circles exceeds theoretical expectation by up to 400-fold. These results can be explained only by greatly enhanced DNA flexibility, not by permanent bends. They invalidate the use of classic theories of flexibility for understanding sharp DNA looping but support predictions of two recent theories. Our findings imply an active role for DNA flexibility in loop formation and suggest that variability in the detailed helical twist of regulatory loops is accommodated naturally by the inherent twistability of the DNA.

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.

Alteration of A.T base-pair opening kinetics by the ammonium cation in DNA A-tracts

We have investigated by NMR the effects of NH(4)(+) on the chemical shifts, on the structure, and on the imino proton exchange kinetics of two duplexes containing an A-tract, [d(CGCGAATTCGCG)](2) and [d(GCA(4)T(4)GC)](2), and of a B-DNA duplex,[d(CGCGATCGCG)](2). Upon NH(4)(+) addition to [d(CGCGAATTCGCG)](2), the adenosine H2 protons, the thymidine imino protons, and the guanosine imino proton of the adjacent G.C pair show unambiguous chemical shifts. Similar shifts are observed in the A-tract of [d(GCA(4)T(4)GC)](2) and for the A5(H2) proton of the B DNA duplex [d(CGCGATCGCG)](2). The localization of the shifted protons suggests an effect related to NH(4)(+) binding in the minor groove. The cross-peak intensities of the NOESY spectra collected at low and high NH(4)(+) concentrations are comparable, and the COSY spectra do not show any change of the sugar pucker. This indicates a modest effect of ammonium binding on the duplex structures. Nevertheless, the imino proton exchange catalysis by ammonia provides evidence for a substantial effect of NH(4)(+) binding on the A.T base-pair kinetics in the A-tracts. Proton exchange experiments performed at high and low NH(4)(+) concentrations show the occurrence of two native conformations in proportions depending on the NH(4)(+) concentration. The base-pair lifetimes and the open-state lifetimes of each conformation are distinct. Exchange from each conformation proceeds via a single open state. But if, and only if, the NH(4)(+) concentration is kept larger than 1 M, the A.T imino proton exchange times of A-tract sequences exhibit a linear dependence versus the inverse of the NH(3) proton acceptor concentration. This had been interpreted as an indication for two distinct base-pair opening modes (Wärmländer, S., Sen, A., and Leijon, M. (2000) Biochemistry 39, 607-615).

A-tract DNA disfavours triplex formation

Optimisation of DNA triplex stability is of fundamental importance in the anti-gene strategy. In the present work, thermal denaturation studies by UV-spectrophotometry and structural and dynamical characterizations by NMR spectroscopy have been used systematically to investigate the effects on triplex stability of isolated insertions of different base triplets into an otherwise homogeneous 15-mer dT x dA-dT oligo-triplex. It is found that insertion of a single central C(+) x G-C or T x D-T triplet (D=2,6-diaminopurine) leads to a pronounced stabilization (up to 20 deg. C if the cytosine base is C5 methylated) at acidic as well as neutral pH. To a smaller degree, this is the case also for a C(+) x I-C triplet insertion. Using imino proton exchange measurements, it is shown that insertion of a DT base-pair in the underlying duplex perturbs the intrinsic A-tract structure in the same way as has been shown for a GC insert. We propose that the intrinsic properties of A-tract duplex DNA (e. g. high propeller twist and rigidity) are unfavourable for triplex formation and that GC- or DT-inserts stabilize the triplex by interfering with the A-tract features of the underlying duplex. The C(+) x I-C triplet without the N2 amino group in the minor groove is readily accommodated within the typical, highly propeller-twisted A-tract structure. This might be related to its smaller effect on the stability of the corresponding triplex. These results may be valuable for understanding DNA triplex formation in vivo as well as for the design of efficient triplex-forming oligonucleotides and in choosing suitable target sequences in the anti-gene strategy.

Additional hydrogen bonds and base-pair kinetics in the symmetrical AMP-DNA aptamer complex

The solution structure of an adenosine monophosphate (AMP)-DNA aptamer complex has been determined previously [Lin, C. H., and Patel, D. J. (1997) Chem. Biol. 4:817-832]. On a symmetrical aptamer complex containing the same binding loop, but with better resolved spectra, we have identified two additional hydrogen bond-mediated associations in the binding loop. One of these involves a rapidly exchanging G imino proton. The phosphate group of the AMP ligand was identified as the acceptor by comparison with other aptamer complexes. Imino proton exchange measurements also yielded the dissociation constants of the stem and binding loop base pairs. This study shows that nuclear magnetic resonance-based imino proton exchange is a good probe for detection of weak hydrogen-bond associations.

Influence of intervening mismatches on long-range guanine oxidation in DNA duplexes

A systematic investigation of the efficiency of oxidative damage at guanine residues through long-range charge transport was carried out as a function of intervening base mismatches. A series of DNA oligonucleotides were synthesized that incorporate a ruthenium intercalator linked covalently to the 5' terminus of one strand and containing two 5'-GG-3' sites in the complementary strand. Single base mismatches were introduced between the two guanine doublet steps, and the efficiency of transport through the mismatches was determined through measurements of the ratio of oxidative damage at the guanine doublets distal versus proximal to the intercalated ruthenium oxidant. Differing relative extents of guanine oxidation were observed for the different mismatches. The damage ratio of oxidation at the distal versus proximal site for the duplexes containing different mismatches varies in the order GC approximately GG approximately GT approximately GA > AA > CC approximately TT approximately CA approximately CT. For all assemblies, damage found with the Delta-Ru diastereomer was found to be greater than with the Lambda-diastereomer. The extent of distal/proximal guanine oxidation in different mismatch-containing duplexes was compared with the helical stability of the duplexes, electrochemical data for intercalator reduction on different mismatch-containing DNA films, and base-pair lifetimes for oligomers containing the different mismatches derived from 1H NMR measurements of the imino proton exchange rates. While a clear correlation is evident both with helix stability and electrochemical data monitoring reduction of an intercalator through DNA films, damage ratios correlate most closely with base-pair lifetimes. Competitive hole trapping at the mismatch site does not appear to be a key factor governing the efficiency of transport through the mismatch. These results underscore the importance of base dynamics in modulating long-range charge transport through the DNA base-pair stack.

Imino proton exchange and base-pair kinetics in RNA duplexes

Using NMR magnetization transfer from water and ammonia-catalyzed exchange of the imino proton, we have measured the base-pair lifetimes and the dissociation constants of six RNA duplexes: [r(CGCGAUCGCG)](2), [r(CGCGAAUUCGCG)](2), [r(CCUUUCGAAAGG)](2), [r(CGCACGUGCG)](2), [r(GGU(8)CC).r(GGA(8)CC)], and [poly(rA).poly(rU)], and we compare them with those of their DNA homologues. As predicted by a two-state (closed/open) model of the pair, the imino proton exchange times decrease linearly vs. the inverse of catalyst concentration. As in DNA duplexes, base pairs open one at a time, and the kinetics is in most cases insensitive to the nature of the adjacent residues. The lifetime of the r(G.C) pairs, 40 to 50 ms, is longer than that of the equivalent in the corresponding oligodeoxynucleotides, and the dissociation constants, about 10(-)(7), are slightly smaller. The r(A.U) opening and closing rates are much larger than those of the d(A.T) pairs, but the stabilities are comparable.

Acid-induced exchange of the imino proton in G.C pairs

Acid-induced catalysis of imino proton exchange in G.C pairs of DNA duplexes is surprisingly fast, being nearly as fast as for the isolated nucleoside, despite base-pair dissociation constants in the range of 10(-5) at neutral or basic pH. It is also observed in terminal G.C pairs of duplexes and in base pairs of drug-DNA complexes. We have measured imino proton exchange in deoxyguanosine and in the duplex (ATATAGATCTATAT) as a function of pH. We show that acid-induced exchange can be assigned to proton transfer from N7-protonated guanosine to cytidine in the open state of the pair. This is faster than transfer from neutral guanosine (the process of intrinsic catalysis previously characterized at neutral ph) due to the lower imino proton pK of the protonated form, 7.2 instead of 9.4. Other interpretations are excluded by a study of exchange catalysis by formiate and cytidine as exchange catalysts. The cross-over pH between the regimes of pH-independent and acid-induced exchange rates is more basic in the case of base pairs than in the mononucleoside, suggestive of an increase by one to two decades in the dissociation constant of the base pair upon N7 protonation of G. Acid-induced catalysis is much weaker in A.T base pairs, as expected in view of the low pK for protonation of thymidine.

NMR studies and restrained-molecular-dynamics calculations of a long A+T-rich stretch in DNA. Effects of phosphate charge and solvent approximations

The nonamer duplex d(GCAAAAACG).d(CGTTTTTGC) was studied by 1H-NMR at 500 MHz. With the exception of the H5' and H5" sugar protons, all protons were assigned by two-dimensional NMR experiments [NOE spectroscopy (NOESY), double-quantum-filtered J-correlated spectroscopy (DQF-COSY) and total correlation spectroscopy (TOCSY)]. The exchange kinetics of the imino protons of the Watson-Crick base pairing were studied at 15 degrees C by measuring inversion-recovery rates under conditions of extensive ammonia base catalysis. Extrapolation to infinite base concentration gave anomalous long lifetimes for the A-tract in accordance with previous results [Leroy, J.-L., Charettier, E., Kochoyan, M. & Guéron, M. (1988) Biochemistry 27, 8894-8898]. On average, 11 NOESY distance constraints/nucleotide were evaluated using the complete relaxation matrix approach. Deoxyribose coupling constants were obtained from simulations of the DQF-COSY cross-peaks, assuming a rapid two-state equilibrium between a C2'-endo and C3'-endo conformer. The sugars were found to be predominantly in the C2'-endo conformation. The NMR-derived distance and torsion constraints were implemented into three different restrained-molecular-dynamics (rMD) protocols, two in vacuo, with different charges on the phosphate group and the third with the solvent explicitly included. All protocols displayed good convergence from different starting structures. The structures derived from the three protocols satisfied experimental restraints equally well and had similar final energies. Although the overall pattern of sequence dependence of helical parameters shows some resemblance in all structures, we find that the absolute amplitudes of the parameters are largely dependent on the rMD protocols, particularly the twist parameters. The minor groove distance P(n + 2)-P(m + 2) varies from 0.7 nm to 1.2 nm in the three protocols. Still the NOESY-derived anomalously short distances AH2(n)-H1'(m + 1) and AH2(n)-H1'(n + 1), n and m denote complementary residues, which are assumed to be indicative of a compressed minor groove, are kept in all calculated structures.

Thermodynamic parameters from hydrogen exchange measurements

Just as exchangeable hydrogens that are controlled by global unfolding can be used to measure thermodynamic parameters at a global level, hydrogens that are exposed to exchange by local unfolding reactions may be used to obtain locally resolved energy parameters. Results with the hemoglobin system demonstrate the ability of HX methods to locate functionally important changes in a protein and to measure the energetic contribution of each. These results offer the promise that HX measurements may be used to delineate, in terms of definable bonds and their energies and interactions, the network of interactions that Hb and other proteins use to produce their various functions.

The interaction of ellipticine derivatives with nucleic acids studied by optical and 1H-nmr spectroscopy: effect of size of the heterocyclic ring system

The DNA interaction of derivatives of ellipticine with heterocyclic ring systems with three, four, or five rings and a dimethylaminoethyl side chain was studied. Optical spectroscopy of drug complexes with calf thymus DNA, poly[(dA-dT).(dA-dT)], or poly[(dG-dC).(dG-dC)] showed a 10 nm bathochromic shift of the light absorption bands of the pentacyclic and tetracyclic compounds upon binding to the nucleic acids, which indicates binding by intercalation. For the tricyclic compound a smaller shift of 1-3 nm was observed upon binding to the nucleic acids. Flow linear dichroism studies show that the geometry of all complexes is consistent with intercalation of the ring system, except for the DNA and poly[(dG-dC).(dG-dC)] complexes of the tricyclic compound, where the average angle between the drug molecular plane and the DNA helix axis was found to be 65 degrees. One-dimensional 1H-nmr spectroscopy was used to study complexes between d(CGCGATCGCG)2 and the tricyclic and pentacyclic compounds. The results on the pentacyclic compound show nonselective broadening due to intermediate chemical exchange of most oligonucleotide resonances upon drug binding. The imino proton resonances are in slow chemical exchange, and new resonances with upfield shifts approaching 1 ppm appear upon drug binding, which supports intercalative binding of the pentacyclic compound. The results on the tricyclic compound show more rapid binding kinetics and very selective broadening of resonances. The data suggest that the tricyclic compound is in an equilibrium between intercalation and minor groove binding, with a preference to bind close to the AT base pairs with the side chain residing in the minor groove.