Molecular Dynamics Simulation Study of DNA Triplex Formed by Mixed Sequences in Solution

Abstract

The unrestrained molecular dynamics simulation of the triple helical DNA with mix sequences d(GACTGGTGAC)?d(CTGACCACTG)*d (GACTGGTGAC), using the particle mesh Ewald sum, is presented here. The Ewald summation method effectively eliminates the usual “cut-off” of the long—range interactions and allowed us to evaluate the full effect of the electrostatic forces. The AMBER5.0 force field has been used during the simulation in solvent. The MD results support a dynamically stable model of DNA triplex over the entire length of the trajectory. The duplex structure assumes the conformation, which is very close to B-DNA. In mixed sequences the purine bases occurs in both strand of DNA duplex. The bases of third strand do not favor the Hoogsteen or/and reverse Hoogsteen type of Hydrogen bonding but they form hydrogen bonds with the bases of both the strand of DNA duplex. The orientation of the third strand is parallel to one of the strand of duplex and all nucleotides (C, A, G & T) show isomorphic behavior with respect to the DNA duplex. The conformation of all the three strands is almost same except few exceptions. Due to interaction of third strand the conformational change in the duplex structure and a finite amount of displacement in the W-C base pairs have been observed. The conformational variation of the back bone torsion angles and helicoidal parameters, groove widths have been discussed. The sequence—dependent effects on local conformation, helicoidal and morphological structure, width of the grooves of DNA helix may have important implication for understanding the functional energetics and specificity of interactions of DNA and its triplexes with proteins, pharmaceutical agents and other legends.     

Triplex hydration: nanosecond molecular dynamics simulation of the solvated triplex formed by mixed sequences

A theoretical model for the hydration pattern and motion of ions around the triple helical DNA with mixed sequences d(GACTGGTGAC)d(GTCACCAGTC)*d(GACTGGTGAC) in solution, during MD simulation, using the particle mesh Ewald sum method, is elaborated here. The AMBER 5.0 force field has been used during the simulation in solvent. The simulation studies support a dynamically stable atmosphere around the DNA triplex in solution over the entire length of the trajectory. The results have been compared with Hoogsteen triplexes and examined in the context of the observed behaviour of hydration in crystallographic data of duplexes. The dynamical organization of counterions and water molecules around the triplex formed by mixed sequences is described here. It has been observed that cations prefer to bind between two adjoining purines of the second and the third strands. The idea of localized complexes (mobile counterions in unspecific electronegative pockets around the DNA triplex with water molecules) may have important implications for understanding the specificity of the interactions of nucleic acids with proteins and other ligands.     

Prediction of the structure of the Y+.R-.R(+)-type DNA triple helix by molecular modelling.

Molecular mechanics has been used to predict the structure of the Y+.R-.R(+)-type DNA triple helix, in which a second polypurine strand binds antiparallel to the homopurine strand of a homopurine/homopyrimidine stretch of duplex DNA. From calculations on the sequence d(C)10.d(G)10.d(G)10, two likely structures emerge. One has the glycosidic torsions of the third strand bases in the anti-conformation and Hoogsteen hydrogen-bonds to the purine strand of the duplex, the other has the third strand purines in the syn orientation and uses a reverse-Hoogsteen hydrogen-bonding pattern. Despite the large structural differences between these two types of triplex, calculations performed in vacuo with a distance-dependent dielectric constant to mimic the shielding effect of solvent show them to be energetically very similar, with the latter (syn) slightly preferred. However, if explicit solvent molecules are included in the calculation, the anti conformation is found to be much preferred. This difference in the results seems to stem from an underestimation of short-range electrostatic interactions in the in vacuo simulations. When TAA or TAT base triples are substituted for the sixth CGG triple in the sequence, it is found that, for the solvated model, the third strand base of the TAA triple prefers the syn orientation while that in the TAT triple retains a preference, though reduced, for the anti conformation.     

Structure of a triple helical DNA with a triplex-duplex junction

Extended purine sequences on a DNA strand can lead to the formation of triplex DNA in which the third strand runs parallel to the purine strand. Triplex DNA structures have been proposed to play a role in gene expression and recombination and also have potential application as antisense inhibitors of gene expression. Triplex structures have been studied in solution by NMR, but have hitherto resisted attempts at crystallization. Here, we report a novel design of DNA sequences, which allows the first crystallographic study of DNA segment containing triplexes and its junction with a duplex. In the 1.8 A resolution structure, the sugar-phosphate backbone of the third strand is parallel to the purine-rich strand. The bases of the third strand associate with the Watson and Crick duplex via Hoogsteen-type interactions, resulting in three consecutive C(+).GC, BU.ABU (BU = 5-bromouracil), and C(+).GC triplets. The overall conformation of the DNA triplex has some similarity to the B-form, but is distinct from both A- and B-forms. There are large changes in the phosphate backbone torsion angles (particularly gamma) of the purine strand, probably due to the electrostatic interactions between the phosphate groups and the protonated cytosine. These changes narrow the minor groove width of the purine-Hoogsteen strands and may represent sequence-specific structural variations of the DNA triplex.     

Structural properties of hybrid triplex of polycation deoxyribonucleic S-methylthiourea (DNmt) strands with a complementary DNA strand, probed by nanosecond molecular dynamics

The replacement of phosphodiester linkages of the polyanion DNA with S-methylthiourea linkers provides the polycation deoxyribonucleic S-methylthiourea (DNmt). Molecular dynamics studies to 1,220 ps of the hybrid triplex formed from octameric DNmt strands d(Tmt)8 with a complementary DNA oligomer strand d(Ap)8 have been carried out with explicit water solvent and Na+Cl- counterions under periodic boundary conditions using the CHARMM force field and the Ewald summation method. The Watson-Crick and Hoogsteen hydrogen-bonding patterns of the A/T tracts remained intact without any structural restraints for triplex structures throughout the simulation. The duplex portion of the triplex structure equilibrated at a B-DNA conformation in terms of the helical rise and other helical parameters. The dynamic structures of the DNmt x DNA x DNmt triplex were determined by examining histograms from the last 800 ps of the dynamics run. These included the hydrogen-bonding pattern (sequence recognition), three-centered bifurcating occurrences, minor groove width variations, and bending of tracts for the hybrid triplex structures. Together with the Watson-Crick hydrogen-bondings, the strong Hoogsteen hydrogen-bondings, the partially maintained three-centered bifurcatings in the Watson-Crick pair, and the medium-strength three-centered bifurcatings in the Hoogsteen pair suggest that the hybrid triplex is energetically favorable as compared to a duplex with similar base stacking, van der Waals interactions, and helical parameters. This is in agreement with our previously reported thermodynamic study, in which only triplex structures were observed in solution. The bending angle measured between the local axis vectors of the first and last helical axis segments is about 20 degrees for the Watson-Crick portion of the averaged structure. Propeller twist (associated with three-centered hydrogen-bonding) up to -30 degrees, native to DNA AT base pairing, was also observed for the triplex structure. The sugar pseudorotation phase angles and the ring rotation angles for the DNA strand are within the C3'-endo domain and C2'-endo domain for the DNmt strand. Water spines are observed in both minor and major grooves throughout the dynamics run. The molecular dynamics simulations of the structural properties of DNmt x DNA x DNmt hybrid triplex is compared to the DNG x DNA x DNG hybrid triplex (In DNG the -O-(PO2-)-O- linkers in DNA is replaced by -NH-C(=N+H2)-NH-).     

Structural Properties of Hybrid Triplex of Polycation Deoxyribonucleic S-methylthiourea (DNmt) Strands with a Complementary DNA Strand,Probed by Nanosecond Molecular Dynamics

Abstract

The replacement of phosphodiester linkages of the polyanion DNA with S-methylthiourea linkers provides the polycation deoxyribonucleic S-methylthiourea (DNmt). Molecular dynamics studies to 1,220 ps of the hybrid triplex formed from octameric DNmt strands d(Tmt)₈ with a complementary DNA oligomer strand d(Ap)₈ have been carried out with explicit water solvent and Na counterions under periodic boundary conditions using the CHARMM force field and the Ewald summation method. The Watson-Crick and Hoogsteen hydrogen-bonding patterns of the A/T tracts remained intact without any structural restraints for triplex structures throughout the simulation. The duplex portion of the triplex structure equilibrated at a B-DNA conformation in terms of the helical rise and other helical parameters. The dynamic structures of the DNmt·DNA·DNmt triplex were determined by examining histograms from the last 800 ps of the dynamics run. These included the hydrogen-bonding pattern (sequence recognition), three-centered bifurcating occurrences, minor groove width variations, and bending of tracts for the hybrid triplex structures. Together with the Watson-Crick hydrogen-bondings, the strong Hoogsteen hydrogen-bondings, the partially maintained three-centered bifurcatings in the Watson-Crick pair, and the medium-strength three-centered bifurcatings in the Hoogsteen pair suggest that the hybrid triplex is energetically favorable as compared to a duplex with similar base stacking, van der Waals interactions, and helical parameters. This is in agreement with our previously reported thermody- namic study, in which only triplex structures were observed in solution. The bending angle measured between the local axis vectors of the first and last helical axis segments is about 20° for the Watson-Crick portion of the averaged structure. Propeller twist (associated with three-centered hydrogen-bonding) up to ?30°, native to DNA AT base pairing, was also observed for the triplex structure. The sugar pseudorotation phase angles and the ring rotation angles for the DNA strand are within the C3′-endo domain and C2′-endo domain for the DNmt strand. Water spines are observed in both minor and major grooves throughout the dynamics run. The molecular dynamics simulations of the structural properties of DNmt·DNA·DNmt hybrid triplex is compared to the DNG·DNA·DNG hybrid triplex (In DNG the -O-(PO₂-)-O- linkers in DNA is replaced by -NH-C(=N₂)-NH-).     

Novel Hoogsteen-like bases for configurational recognition of the T-A base pair by DNA triplex formation

Effective sequence-specific recognition of duplex DNA is possible by triplex formation with natural oligonucleotides via Hoogsteen H-bonding. However, triplex formation is in practice limited to pyrimidine oligonucleotides binding duplex A-T or G-C base-pair DNA sequences specifically at homopurine sites in the major groove as T·A-T and C⁺·G-C triplets. Here we report the successful modeling of novel unnatural nucleosides that recognize the T-A DNA base pair by Hoogsteen interaction. Since the DNA triplex can be considered to assume an A-type or B-type conformation, these novel Hoogsteen nucleotides are tested within model A-type and B-type conformation triplex structures. A triplet consisting of the T-A base pair and one of the novel Hoogsteen nucleotides replaces the central T·A-T triplet in the triplex using the same deoxyribose-phosphodiester and base-deoxyribose dihedral angle configuration. The entire triplex is energy minimized and the presence of any structural or energetic perturbations due to the central triplet is assessed with respect to the unmodified energy-minimized (T·A-T)₁₁ proposed starting structures. Incorporation of these novel triplets into both A-type and B-type natural triplex structures provokes minimal change in the configuration of the central and adjacent triplets. The plan is to produce a series of Hoogsteen-like bases that preferentially bind the T-A major groove in either an A-type or B-type conformation. Selective recognition of the T-A major groove with respect to the G-C major groove, which presents similar keto and amine placement, is also assessed with configurational preference. Evaluation of the triplex solution structure by using these unnatural bases as binding conformational probes is a prerequisite to the further design of triplet forming bases. © 1996 John Wiley & Sons, Inc.     

Theoretical study of the Hoogsteen-Watson-Crick junctions in DNA

A series of d (AT)(n) oligonucleotides containing mixtures of normal B-type Watson-Crick and antiparallel Hoogsteen helices have been studied using molecular dynamics simulation techniques to analyze the structural and thermodynamic impact of the junction between Watson-Crick and antiparallel Hoogsteen structures. Analysis of molecular dynamics simulations strongly suggests that for all oligonucleotides studied the antiparallel Hoogsteen appears as a reasonable conformation, only slightly less stable than the canonical B-type Watson-Crick one. The junctions between the Watson-Crick and Hoogsteen structures introduces a priori a sharp discontinuity in the helix, because the properties of each type of conformation are very well preserved in the corresponding fragments. However, and quite counterintuitively, junctions do not largely distort the duplex in structural, dynamics or energetic terms. Our results strongly support the possibility that small fragments of antiparallel Hoogsteen duplex might be embedded into large fragments of B-type Watson-Crick helices, making possible protein-DNA interactions that are specific of the antiparallel Hoogsteen conformation.     

Three-dimensional structures of bulge-containing DNA fragments.

The three-dimensional structure of a DNA tridecamer d(CGCAGAATTCGCG)2 containing bulged adenine bases was determined by single crystal X-ray diffraction methods, at 120 K, to 2.6 A resolution. The structure is a B-DNA type double helix with a single duplex in the asymmetric unit. One of the bulged adenine bases loops out from the double helix, while the other stacks in to it. This is in contrast to our preliminary finding, which indicated that both adenine bases were looped out. This revised model was confirmed by the use of a covalently bound heavy-atom derivative. The conformation of the looped-out bulge hardly disrupts base stacking interactions of the bases flanking it. This is achieved by the backbone making a "loop-the-loop" curve with the extra adenine flipping over with respect to the other nucleotides in the strand. The looped-out base intercalates into the stacked-in bulge site of a symmetrically related duplex. The looped-out and stacked-in bases form an A.A reversed Hoogsteen base-pair that stacks between the surrounding base-pairs, thus stabilizing both bulges. The double helix is frayed at one end with the two "melted" bases participating in intermolecular interactions. A related structure, of the same tridecamer, after soaking the crystals with proflavin, was determined to 3.2 A resolution. The main features of this B-DNA duplex are basically similar to the native tridecamer but differ in detail especially in the conformation of the bulged-out base. Accommodation of a large perturbation such as that described here with minimal disruption of the double helix shows both the flexibility and resiliency of the DNA molecule.     

FTIR and UV spectroscopy studies of triplex formation between alpha-oligonucleotides with non-ionic phoshoramidate linkages and DNA targets

The triplexes formed by pyrimidine alpha-oligodeoxynucleotides, 15mers alpha dT(15) or 12mers alpha dCT having dimethoxyethyl (PNHdiME), morpholino (PMOR) or propyl (PNHPr) non-ionic phosphoramidate linkages with DNA duplex targets have been investigated by UV and FTIR spectroscopy. Due to the decrease in the electrostatic repulsion between partner strands of identical lengths all modifications result in triplexes more stable than those formed with unmodified phosphodiester beta-oligodeoxynucleotides (beta-ODNs). Among the alpha-ODN third strands having C and T bases and non-ionic phosphoramidate linkages (alpha dCTPN) the most efficient modification is (PNHdiME). The enhanced third strand stability of the alpha dCTPN obtained as diastereoisomeric mixtures is attenuated by the steric hindrance of the PMOR linkages or by the hydrophobicity of the PNHPr linkages. All alpha dCTPN strands form triplexes even at neutral pH. In the most favorable case (PNHdiME), we show by FTIR spectroscopy that the triplex formed at pH 7 is held by Hoogsteen T*A.T triplets and in addition by an hydrogen bond between O6 of G and C of the third strand (Tm = 30 degrees C). The detection of protonated cytosines is correlated at pH 6 with a high stabilization of the triplex (Tm = 65 degrees C). While unfavorable steric effects are overcome with alpha anomers, the limitation of the pH dependence is not completely suppressed. Different triplexes are evidenced for non pH dependent phosphoramidate alpha-thymidilate strands (alpha dT(15)PN) interacting with a target duplex of identical length. At low ionic strength and DNA concentration we observe the binding to beta dA(15) either of alpha dT(15)PN as duplex strand and beta dT(15) as third strand, or of two hydrophobic alpha dT(15)PNHPr strands. An increase in the DNA and counterion concentration stabilizes the anionic target duplex and then the alpha dT(15)PN binds as Hoogsteen third strand.     

Sequence-dependent elastic properties of DNA

Harmonic elastic constants of 3-11 bp duplex DNA fragments were evaluated using four 5 ns unrestrained molecular dynamics simulation trajectories of 17 bp duplexes with explicit inclusion of solvent and counterions. All simulations were carried out with the Cornell et al. force-field and particle mesh Ewald method for long-range electrostatic interactions. The elastic constants including anisotropic bending and all coupling terms were derived by analyzing the correlations of fluctuations of structural properties along the trajectories. The following sequences have been considered: homopolymer d(ApA)(n) and d(GpG)(n), and alternating d(GPC)(n) and d(APT)(n). The calculated values of elastic constants are in very good overall agreement with experimental values for random sequences. The atomic-resolution molecular dynamics approach, however, reveals a pronounced sequence-dependence of the stretching and torsional rigidity of DNA, while sequence-dependence of the bending rigidity is smaller for the sequences considered. The earlier predicted twist-bend coupling emerged as the most important cross-term for fragments shorter than one helical turn. The calculated hydrodynamic relaxation times suggest that damping of bending motions may play a role in molecular dynamics simulations of long DNA fragments. A comparison of elasticity calculations using global and local helicoidal analyses is reported. The calculations reveal the importance of the fragment length definition. The present work shows that large-scale molecular dynamics simulations represent a unique source of data to study various aspects of DNA elasticity including its sequence-dependence.     

UV spectroscopic identification and thermodynamic analysis of protonated third strand deoxycytidine residues at neutrality in the triplex d(C(+)-T)6:[d(A-G)6.d(C-T)6]; evidence for a proton switch.

Near-UV difference spectral analysis of the triplex formed from d(C-T)6 and d(A-G)6.d(C-T)6 in neutral and acidic solution shows that the third strand dC residues are protonated at pH 7.0, far above their intrinsic pKa. Additional support for ion-dipole interactions between the third strand dC residues and the G.C target base pairs comes from reduced positive dependence of triplet stability on ionic strength below 0.9 M Na+, inverse dependence above 0.9 M Na+ and strong positive dependence on hydrogen ion concentration. Molecular modeling (AMBER) of C:G.C and C+:G.C base triplets with the third strand base bound in the Hoogsteen geometry shows that only the C+:G.C triplet is energetically feasible. van't Hoff analysis of the melting of the triplex and target duplex shows that between pH 5.0 and 8.5 in 0.15 M NaCl/0.005 M MgCl2 the enthalpy of melting (delta H degree obs) varies from 5.7 to 6.6 kcal.mol-1 for the duplex in a duplex mixture and from 7.3 to 9.7 kcal.mol-1 for third strand dissociation in the triplex mixture. We have extended the condensation-screening theory of Manning to pH-dependent third strand binding. In this development we explicitly include the H+ contribution to the electrostatic free energy and obtain [formula: see text]. The number of protons released in the dissociation of the third strand from the target duplex at pH 7.0, delta n2, is thereby calculated to be 5.5, in good agreement with approximately six third strand dc residues per mole of triplex. This work shows that when third strand binding requires protonated residues that would otherwise be neutral, triplex formation and dissociation are mediated by proton uptake and release, i.e., a proton switch. As a by-product of this study, we have found that at low pH the Watson-Crick duplex d(A-G)6.d(C-T)6 undergoes a transition to a parallel Hoogsteen duplex d(A-G)6.d(C(+)-T)6.     

Proton nuclear magnetic resonance assignments and structural characterization of an intramolecular DNA triplex.

Two-dimensional 1H n.m.r. spectroscopy has been used to study the 31-base DNA oligonucleotide 5'-dAGAGAGAACCCCTTCTCTCTTTTTCTCTCTT-3', which folds to form a stable intramolecular triplex in solution at acidic pH. This structure is considerably more difficult to assign than short B-DNA duplexes and requires new assignment methods. The assignment strategy and assignments of almost all of the exchangeable and nonexchangeable resonances are presented. Seven base triplets and one Watson-Crick base-pair form the core of the structure and are connected by a four C and four T loop at either end. The second pyrimidine "strand" (bases 24 to 31) in this intramolecular pyrimidine-purine-pyrimidine triplex binds via Hoogsteen base-pairs in the major groove and is parallel to the purine "strand" (bases 1 to 8). Analysis of the sugar puckers reveals that, contrary to widely accepted belief, the triplex sugars are not predominantly in the N-type (close to C3'-endo) conformation. Except for some of the C nucleotides, all sugars are predominantly S-type (close to C2'-endo). Thus, the duplex DNA does not assume N-type sugar conformations to accommodate a third strand in the major groove. A preliminary model of the triplex structure is presented.     

The molecular structure of Okazaki fragment r(CGCA)d(AAAAAGCG):d(CGCTTTTTTGCG) in solution

The hetero duplex molecule, r(CGCA)d(AAAAAGCG):d(CGCTTTTTTGCG) which corresponds to Okazaki fragment was synthesized and its molecular structure has been analyzed by NMR study. The RNA strand of RNA-DNA hybrid region adopts A-form and DNA strand of the same region deviates from the standard B-form. The conformation of DNA-DNA duplex segment belongs to B-form. The hybrid-DNA duplex junction shows a structural discontinuities, A-B junction. The same conformational characteristic of oligo(dA): oligo(dT) tract as that of DNA oligomer which has same base sequence has been observed.     

Thermodynamics of triple helix formation: spectrophotometric studies on the d(A)10.2d(T)10 and d(C+3T4C+3).d(G3A4G3).d(C3T4C3) triple helices.

We have stabilized the d(A)10.2d(T)10 and d(C+LT4C+3).d(G3A4G3).d(C3T4C3) triple helices with either NaCl or MgCl2 at pH 5.5. UV mixing curves demonstrate a 1:2 stoichiometry of purine to pyrimidine strands under the appropriate conditions of pH and ionic strength. Circular dichroic titrations suggest a possible sequence-independent spectral signature for triplex formation. Thermal denaturation profiles indicate the initial loss of the third strand followed by dissociation of the underlying duplex with increasing temperature. Depending on the base sequence and ionic conditions, the binding affinity of the third strand for the duplex at 25 degrees C is two to five orders of magnitude lower than that of the two strands forming the duplex. Thermodynamic parameters for triplex formation were determined for both sequences in the presence of 50 mM MgCl2 and/or 2.0 M NaCl. Hoogsteen base pairs are 0.22-0.64 kcal/mole less stable than Watson-Crick base pairs, depending on ionic conditions and base composition. C+.G and T.A Hoogsteen base pairs appear to have similar stability in the presence of Mg2+ ions at low pH.     

Homopurine and homopyrimidine strands complementary in parallel orientation form an antiparallel duplex at neutral pH with A-C,G-T,and T-C mismatched base pairs

DNA sequences d-TGAGGAAAGAAGGT (a 14-mer) and d-CTCCTTTCTTCC (a 12-mer) are complementary in parallel orientation forming either Donahue (reverse Watson-Crick) base pairing at neutral pH or Hoogsteen base pairing at slightly acidic pH. The structure of the complex formed by dissolving the two strands in equimolar ratio in water has been investigated by nmr. At neutral pH, the system forms an ordered antiparallel duplex with five A : T and four G : C Watson-Crick base pairs and three mismatches, namely G-T, A-C, and T-C. The nuclear Overhauser effect cross-peak pattern suggests an overall B-DNA conformation with major structural perturbations near the mismatches. The duplex has a low melting point and dissociates directly into single strands with a broad melting profile. The hydrogen-bonding schemes in the mismatched base pairs have been investigated. It has been shown earlier that in acidic pH, the system prefers a triple-stranded structure with two pyrimidine strands and one purine strand. One of the pyrimidine strands has protonated cytosines, forms Hoogsteen base pairing, and is aligned parallel to the purine strand; the other has nonprotonated cytosines and has base-pairing scheme similar to the one discussed in this paper. The parallel duplex is therefore less stable than either the antiparallel duplex or the triplex, in spite of its perfect complementarity. © 1997 John Wiley & Sons, Inc. Biopoly 41: 773–784, 1997     

Site-specific intercalation at the triplex-duplex junction induces a conformational change which is detectable by hypersensitivity to diethylpyrocarbonate.

Using site-specific intercalation directed by intermolecular triplex formation, the conformation of an intercalation site in DNA was examined by footprinting with the purine-specific (A much greater than G) reagent diethylpyrocarbonate. Site specific intercalation was achieved by covalently linking an intercalator to the 5' end of a homopyrimidine oligodeoxynucleotide, which bound to a homopurinehomopyrimidine stretch in a recombinant plasmid via intermolecular triplex formation. This directs intercalation to a single site in 3kb of DNA at the 5' triplex-duplex junction. Footprinting with diethylpyrocarbonate and dimethylsulphate revealed strong protection from modification of adenine residues within the triple-helix in concordance with their Hoogsteen pairing with the third strand, and a strong hypersensitivity to diethylpyrocarbonate at the first adenine of the duplex. This result indicates that intercalation at this site induces a conformational change at the 5' triplex-duplex junction. Furthermore, the same diethlypyrocarbonate hypersensitivity was observed with an unmodified triple-strand forming oligonucleotide and a range of intercalating molecules present in solution. Thus the 5' triplex-duplex junction is a strong binding site for some intercalating molecules and the junction undergoes a conformational change which is sensitive to diethylpyrocarbonate upon insertion of the planar aromatic chromophore. This conformational change can be used to direct a single-strand cut in duplex DNA to a defined site.     

Calculated distortions of duplex DNA by a cis, syn cyclobutane thymine dimer are unaffected by a 3' TpA step.

Molecular dynamics simulations were performed on the duplex DNA dodecamers d(CGCGAA TT CGCG): d(CGCGAATTCGCG) and d(GCACGAA TT AAG): d(CTTAATTCGTGC), where TT denotes a cis, syn cyclobutane thymine dimer. The constant temperature and pressure algorithm of the AMBER 4.1 molecular-modeling package was used with explicit water and counterions, periodic boundary conditions and electrostatic interactions evaluated by the particle-mesh Ewald method. Results were analyzed by the CURVES algorithm and its implementation in DIALS and WINDOWS. Calculated distortions of DNA structure by the thymine dimer were qualitatively and quantitatively similar for the two sequences. Despite the enhanced flexibility of the native TpA dinucleotide step, major deviations from the B-DNA values of helicoidal parameters were found only at the Ap and p dinucleotide steps in both sequences. Only the AT base pairs of the two sequences that contain the 5' thymine of the dimers exhibited weakened Watson-Crick hydrogen bonds and anomalous stretching. Hence, we conclude that the pattern of structural perturbations responsible for recognition of cis, syn thymine dimers by repair enzymes is not sensitive to their sequence context.     

Base-pair dynamics in an antiparallel DNA triplex measured by catalyzed imino proton exchange monitored via 1H NMR spectroscopy

Using (1)H NMR spectroscopy, the base-pair opening dynamics of an antiparallel foldback DNA triplex and the corresponding duplex has been characterized via catalyzed imino proton exchange. The triplex system was found to be in an equilibrium between a duplex and a triplex form. The exchange rate between the two forms (i.e., the on/off-rate of the third strand) was measured to be 5 s(-1) at 1 degrees C, and the base-pair dynamics of both forms were investigated separately. Both Watson-Crick and reverse Hoogsteen base pairs were found to have base-pair lifetimes in the order of milliseconds. The stability of the Watson-Crick base pairs was, however, substantially increased in the presence of the third strand. In the DNA triplex, the opening dynamics of the reverse Hoogsteen base pairs was significantly faster than the dynamics of the Watson-Crick pairs. We were able to conclude that, for both Watson-Crick and reverse Hoogsteen base pairs, spontaneous and individual opening from within the closed base triplet is the dominating opening pathway.     

DNA conformation and dynamics

Nucleotide conformation and dynamics are important for the study of radiation damage to DNA at the atomic level. It is necessary to study not only normal oligonucleotide structure but also those containing modified bases which result from interaction with OH-radicals. There are now over 8000 atomic coordinate entries in the Brookhaven Protein Data Bank, of which over 900 relate to experimentally determined structures of nucleic acids and nucleic acid/protein complexes. We review some of these data which have led to the elucidation of novel DNA conformations, insight into DNA sequence specificity and knowledge of protein/DNA interactions. Further understanding of the conformation, stability and dynamics of nucleic acids has come from molecular modelling. We have used such techniques to study chemical modifications to bases such as alkylation of thymine and guanine and the effects of curvature in longer sequences. Recent improvements in this area include the inclusions of explicit counter-ions and solvent molecules, the use of Particle Mesh Ewald methods to incorporate the long-range electrostatic interactions and the use of longer time scale simulations. We have employed these methods to analyse the effects of incorporation of 8-oxodeoxyguanosine into duplex DNA. This lesion is a common result of radiation damage and is known to have important effects in mutagenesis, cancer and ageing. Received: 7 October 1998 / Accepted in revised form: 18 January 1999