Molecular mechanics and dynamics of an abasic frameshift in DNA and comparison to NMR data

Dynamic conformational states of DNA containing T.T or BrdU.T mispaired bases: wobble H-bond pairing versus cross-strand inter-atomic contacts

The dynamic structure of 11-mer DNA duplexes of different sequences with or without homopyrimidine (T.T, or BrdU.T) mismatches was studied by molecular dynamics (MD) simulations on a time scale from 200 ps to 1 ns. The conformational analysis suggests that in mismatched duplexes the formation of classical T.T wobble H-bonding pairing is nearest-neighbor sequence-dependent and, in most cases, three-centered H-bonds and numerous alternative close cross-strand interatomic contacts exist. Thus, in duplex W1, where the central triplet is 5'd(CTA).d(TTG), two wobble conformations W upward arrow (alphabeta) and W downward arrow (betaalpha) are formed and exchange rapidly at 300 K. In contrast, when the central triplet is 5'd(TTT).d(ATA) (W2 duplex) wobble conformations are rarely observed at 300 K, and the T.T mispair most often adopts a "twisted" conformation with one largely persistent normal H-bond, plus a stable cross-strand contact involving a T flanking base. However, at elevated temperature (400 K) the same W2 duplex shows frequent exchange between the two classical wobble conformations (alphabeta<-->betaalpha), as is in the case when the central triplet is 5'd(TBrdUT).d(ATA) (W3 duplex at 300 K). It is suggested that in the W2 sequence, restrictions due to thymine-methyl/pi interactions prevent the formation of wobble pairing and thermal activation energy, and/or the chemical replacement of T by BrdU are required in order for the T(BrdU).T mismatch to adopt and exchange between wobble conformations. The specific short and/or long-lived (double/triple) cross-strand dynamic interactions in W1, W2 and W3 duplexes are throughout characterized. These frequent atomic encounters exemplify possible inter-strand charge transfer pathways in the studied DNA molecules.

Structural study of DNA duplex containing an N-(2-deoxy-beta-D-erythro-pentofuranosyl) formamide frameshift by NMR and restrained molecular dynamics

The presence of an N-(2-deoxy-beta-D-erythro-pentofuranosyl) formamide (F) residue, a ring fragmentation product of thymine, in a frameshift context in the sequence 5'-d-(AGGACCACG)*d(CGTGGFTCCT) has been studied by 1H and 31P nuclear magnetic resonance (NMR) and molecular dynamics. Two-dimensional NMR studies show that the formamide residue, whether the cis or trans isomer, is rotated out of the helix and that the bases on either side of the formamide residue in the sequence, G14 and T16, are stacked over each other in a way similar to normal B-DNA. The cis and trans isomers were observed in the ratio 3:2 in solution. Information extracted from 31P NMR data reveal a modification of the phosphodiester backbone conformation at the extrahelical site, which is also observed during the molecular dynamics simulations.

Biopolymer Chain Elasticity: A novel concept and a least deformation energy principle predicts backbone and overall folding of DNA TTT hairpins in agreement with NMR distances

A new molecular modelling methodology is presented and shown to apply to all published solution structures of DNA hairpins with TTT in the loop. It is based on the theory of elasticity of thin rods and on the assumption that single-stranded B-DNA behaves as a continuous, unshearable, unstretchable and flexible thin rod. It requires four construction steps: (i) computation of the tri-dimensional trajectory of the elastic line, (ii) global deformation of single-stranded helical DNA onto the elastic line, (iii) optimisation of the nucleoside rotations about the elastic line, (iv) energy minimisation to restore backbone bond lengths and bond angles. This theoretical approach called 'Biopolymer Chain Elasticity' (BCE) is capable of reproducing the tri-dimensional course of the sugar-phosphate chain and, using NMR-derived distances, of reproducing models close to published solution structures. This is shown by computing three different types of distance criteria. The natural description provided by the elastic line and by the new parameter, Omega, which corresponds to the rotation angles of nucleosides about the elastic line, offers a considerable simplification of molecular modelling of hairpin loops. They can be varied independently from each other, since the global shape of the hairpin loop is preserved in all cases.

Solution structure of a DNA duplex containing a formamide-adenine base pair

The N-(2-deoxy-beta3-D-erythro-pentofuranosyl) formamide residue results from a ring fragmentation product of thymine or cytosine. The presence of a formamide-adenine base pair in the sequence 5'd(AGGAACCACG).d(CGTGGFTCCT) has been studied by 1H and 31P nuclear magnetic resonance (NMR) and molecular dynamics. There are two possible isomers for the formamide side chain, either cis or trans. For each isomer, we observed an equilibrium in solution between two forms. First, a species where the formamide is intrahelical and paired with the facing adenine. For the cis isomer, the formamide is in a syn conformation and two hydrogen bonds with adenine are formed. The trans isomer is in an anti conformation and a single hydrogen bond is observed. In the second form, whatever the isomer, the formamide is rejected outside the helix, whereas the adenine remains inside.

Abasic DNA structure, reactivity, and recognition

Loss of a base in DNA, i.e., creation of an abasic site leaving a deoxyribose residue in the strand, is a frequent lesion that may occur spontaneously, or under the action of radiations and alkylating agents, or enzymatically as an intermediate in the repair of modified or abnormal bases. The abasic site lesion is mutagenic or lethal if not repaired. From a chemical point of view,the abasic site is an alkali-labile residue that leads to strand breakage through beta- and delta- elimination. Progress in the understanding of the chemistry and enzymology of abasic DNA largely relies upon the study of synthetic abasic duplexes. Several efficient synthetic methods have thus been developed to introduce the lesion (or a stable analogue) at defined position in the sequence. Physicochemical and spectroscopic examination of such duplexes, including calorimetry, melting temperature, high-field nmr and molecular modeling indicate that the lesion strongly destabilizes the duplex, although remaining in the canonical B-form with structural modifications strictly located at the site of the lesion. Probes have been developed to titrate the damage in DNA in vitro. Series of molecules have been devised to recognize specifically the abasic site, exhibiting a cleavage activity and mimicking the AP nucleases. Others have been prepared that bind strongly to the abasic site and show promise in potentiating the cytotoxic and antitumor activity of the clinically used nitrosourea (bis-chloroethylnitrosurea).

DNA structural forms

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

A database of 32 DNA triplets to study triple helices by molecular mechanics and dynamics

We present here a database of 32 deoxyribonucleotide triplets, that can be used as building blocks of triple helix forming deoxyribonucleotides on a computer. This database is made of all the pairing schemes of the triplets ATT, GCC+, ATA and GCG where the third base forms two hydrogen bonds with the purine of the first two Watson-Crick strands. The essential features of the known triple helices were preserved in the resulting structures. A triple helix can be easily built from any combination of these basic triplets. Four homogeneous and alternate triple helices thus obtained were studied by molecular mechanics and dynamics in vacuo. The results are in agreement with known experimental observations for ATT and suggest a possible structure for the GCG triple helix. In order to characterize the geometry of the structures obtained, the definitions of nucleic acid structure parameters (R.E. Dickerson et al., EMBO J. 8 (1989) 1-4) have been extended to triple helical polynucleotides.

Triple helix formation with short oligonucleotide-intercalator conjugates matching the HIV-1 U3 LTR end sequence

In an attempt to target short purine sequences in view of pharmacological application, we have synthesized three new TFO (triple-helix-forming oligonucleotide) conjugates in which an intercalating oxazolopyridocarbazole (OPC) chromophore is linked by a pentamethylene linker to a 7-mer oligonucleotide matching the polypurine/polypyrimidine sequence located in the HIV-1 U3 LTR end region. The TFO moiety of conjugates are 5'CCTTCCC, 5'GGGAAGG, and 5'GGGTTGG. Their ability to bind to double-stranded DNA targets was examined. This binding is demonstrated by a footprinting technique using DNase I as a cleaving agent. The complex involved intermolecular pyr-pur*pyr or pur-pur*pyr triple helix. Pyrimidine TFO-OPC binds in a pH-dependent manner, whereas the others do not. The formation of the complex has been investigated at neutral pH and increasing temperature. We observed that the protection due to the purine and mixed TFO-OPC was pH independent and remained identical up to 40 degrees C. To determine the position of the OPC chromophore, molecular modeling was undertaken on the purine-conjugate/target complex. It has been suggested that the complex involved the intercalation of the OPC at the triplex-duplex junction with a small unwinding at the next excluded site.

The pH dependent configurations of the C.A mispair in DNA

The structure of the cytosine-adenine mispair in a 7 base pair duplex has been investigated by proton NMR spectroscopy. At low pH, the predominant structure is protonated on the A residue and assumes a wobble conformation consistent with previous findings. The C residue of the mispair is found in a C2'-C3' endo equilibrium. This is confirmed by molecular dynamics calculations which suggest that the conformation of the protonated wobble is flexible and not as rigid as a normal base pair. As the solution pH is raised, a structural transition is observed with an apparent pK of 7.54 at 23 degrees C. At higher pH the predominant structure is one in which both the C and A residues are intrahelical. Evidence is presented that this structure corresponds to a reverse wobble in which the two bases are held together by one hydrogen bond. This structure is much less stable than the protonated form and even at low temperature single strands are observed in slow exchange with the neutral duplex form.

Solution conformation of an oligonucleotide containing a G.G mismatch determined by nuclear magnetic resonance and molecular mechanics

We have determined by two-dimensional nuclear magnetic resonance studies and molecular mechanics calculations the three dimensional solution structure of the non-selfcomplementary oligonucleotide, d(GAGGAGGCACG). d(CGTGCGTCCTC) in which the central base pair is G.G. This is the first structural determination of a G.G mismatch in a oligonucleotide. Two dimensional nuclear magnetic resonance spectra show that the bases of the mismatched pair are stacked into the helix and that the helix adopts a classical B-DNA form. Spectra of the exchangeable protons show that the two guanosines are base paired via their imino protons. For the non-exchangeable protons and for some of the exchangeable protons nuclear Overhauser enhancement build up curves at short mixing times have been measured. These give 84 proton-proton distances which are sensitive to the helix conformation. One of the guanosines adopts a normal anti conformation while the other is syn or close to syn. All non-terminal sugars are C2' endo. These data sets were incorporated into the refinement of the oligonucleotide structure by molecular mechanics calculations. The G.G mismatch shows a symmetrical base pairing structure. Although the mismatch is very bulky many of its features are close to that of normal B-DNA. The mismatch induces a small lateral shift in the helix axis and the sum of the helical twist above and below the mismatch is close to that of B-DNA.