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.

Multiple structures for the 2-aminopurine-cytosine mispair.

The base pair formed between 2-aminopurine (2AP) and cytosine (C) is an intermediate in transition mutations generated by 2AP. To date, several structures have been proposed for the 2AP-C mispair, including those involving a rare tautomer, a protonated base pair, and a neutral wobble structure. In this paper, we describe a series of UV, fluorescence, and NMR studies which demonstrate that an equilibrium exists between the neutral wobble and the protonated Watson-Crick structures. The apparent pK value for the transition between the structures is 5.9-6.0. Formation of a Watson-Crick base pair is accomplished predominantly by protonation of the 2AP residue as indicated by UV spectral changes, fluorescence quenching, and changes in proton chemical shifts. Rapid transfer of the shared proton between the 2AP and cytosine residues is indicated by the rapid exchange of the cytosine amino protons of the protonated Watson-Crick configuration. The relative contribution of the neutral wobble and protonated Watson-Crick configurations to 2AP-induced transition mutations is discussed.

Solution structure by NMR and molecular dynamics of a duplex containing a guanine opposite a N-(2-deoxy-beta-D-erythro-pentofuranosyl)formamide lesion.

One- and two-dimensional NMR spectroscopy has been used combined with molecular dynamics to determine the fine structure of the DNA duplex 5'-d(AGGAGCCACG).d(CGTGGFTCCT) where F is the N-(2-deoxy-beta-D-erythro-pentofuranosyl)formamide residue which is a ring fragmentation product of thymine. The formamide deoxyribose exists as two isomers with respect to the orientation about the peptide bond. The two isomers (trans and cis) are observed in a ratio 3:2 in solution. For both species, the oligonucleotide adopts a globally B form structure although conformational changes are observed around the mismatch site. The formamide residue, whatever the isomer, is intrahelical and can pair with the guanine on the opposite strand with one hydrogen bond. For the cis isomer, the residue adopts a syn orientation and is able to form a second hydrogen bond with the guanine on the 5' side on the same strand. Off-resonance ROESY experiments have been used to investigate the chemical exchange observed at low temperature of the duplex. Conformational exchange has only been found for the oligonucleotide with the formamide residue in the trans conformation.

Structure of an oligonucleotide containing a N-(2-deoxy-beta-D-erythro-pentofuranosyl)formamide residue facing a guanine.

Formamide residue is a major oxidative DNA damage product from ionizing radiation on thymine residues in DNA. We report NMR and molecular modeling studies on a DNA duplex structure which contains guanine opposite formamide residue. Formamide residue exists as either the cis and trans isomer. For the trans and the cis isomers, we find that guanine and formamide are stacked inside the helix and are hydrogen bonded. The oligonucleotide adopts globally a B form structure for the two isomers. Conformational changes are observed between the two isomers.

Solution structures of a duplex containing an adenine opposite a gap (absence of one nucleotide). An NMR study and molecular dynamic simulations with explicit water molecules.

We investigated the behaviour of a 15mer DNA duplex, [5'd(CAGAGTCACTGGCTC)3']. [5'd(GAGCCAG)3' + 5'd(GACTCTG)3'] which contained an adenine opposite the gap. Analysis of the NMR data showed the existence of one major species, which was in equilibrium with two minor species. Their relative concentrations varied as a function of pH with a pKa of approximately 4.5. For the major species, the duplex was globally in B conformation with the central adenine stacked in the helix. The two G.C base pairs adjacent to the central adenine were well formed and a gap was present in front of this adenine. For the minor species, major structural perturbations occurred in the centre of the duplex. At neutral pH, the central adenine was involved in a G.A mismatch with G23 adjacent to the gap. Cytosine C7 was then extrahelical and no gap was observed. Under these conditions, the major neutral species corresponded to 70% of the total and the minor species to 30%. At acidic pH, the central adenine of the minor species was protonated and was involved in a G(syn).A+(anti) mismatch. The difference is that C9 is now extrahelical and G22 is implicated in the mispair. Three-dimensional models were built to initiate molecular dynamic simulations, which were in good agreement with the NMR data. Their structural stability in terms of hydrogen bonding and their flexibility are discussed and the biological significance for the interaction with DNA polymerase is evoked.

Conformations of nicked and gapped DNA structures by NMR and molecular dynamic simulations in water.

We have analyzed and compared the molecular structures and dynamics of DNA duplexes containing a nick or a gap of one nucleotide where the base in front of the gap is a guanine. The continuous strand has the sequence 5'(CAGAGTCXCTGGCTC) where the residue X is absent for the nick, 14-mer, and where it is a G residue for the gap. Duplexes were formed with the two corresponding 7-mers. Neither of these is phosphorylated adjacent at the nick site, but it is a good model for a single strand break. For the nick structure, the quantitative NMR data show that the global conformation is very close to canonical B-form DNA, but it displays enhanced local flexibility. For the gap structure, we observe only one species in which the extra G is well stacked into the helix. The two half-helices around this residue also show a B-form conformation. As with the nick duplex, the adjacent G imino protons show enhanced exchange with solvent. The gap does not close completely. Using distance constraints, MD calculations show that the nick conformation is very close to a duplex with no lesion but is indeed more flexible in the central part. The gapped structure shows two families of conformations. One is close to B-DNA, the other is significantly kinked at the gap which reduces the size of the cavity. We observe a spine of hydration within the cavities, similar, but of different geometry in the two cases.

Solution structure of N-(2-deoxy-D-erythro-pentofuranosyl)urea frameshifts, one intrahelical and the other extrahelical, by nuclear magnetic resonance and molecular dynamics.

The presence of a N-(2-deoxy-D-erythro pentofuranosyl)urea (henceforth referred to as deoxyribosylurea) residue, ring fragmentation product of a thymine, in a frameshift situation in the sequence 5'd(AGGACCACG).d(CGTGGurTCCT) has been studied by 1H and 31P nuclear magnetic resonance and molecular dynamics. At equilibrium, two species are found in slow exchange. We observe that the deoxyribosylurea residue can be either intra- or extrahelical within structures which otherwise do not deviate strongly from that of a B-DNA as observed by NMR. Our study suggests that this is determined by the nature and number of hydrogen bonds which this residue can form as a function of two possible isomers. There are two possible structures for the urea side chain, either cis or trans for the urido bond which significantly changes the hydrogen bonding geometry of the residue. In the intrahelical species, the cis isomer can form two good hydrogen bonds with the bases on the opposite strand in the intrahelical species, A4 and C5, which is not the case for the trans isomer. This results in a kink in the helical axis. For the major extrahelical species, the situation is reversed. The trans isomer is able to form two good hydrogen bonds, with G13 on the same strand and A7 on the opposite strand. For the extrahelical species, the cis isomer can form only one hydrogen bond. In this major structure the NMR data show that the bases which are on either side of the deoxyribosylurea residue in the sequence, G14 and T16, are stacked over each other in a way similar to a normal B-DNA structure. This requires the formation of a loop for the backbone between these two residues. This loop can belong to one of two families, right- or left-handed. In a previous study of an abasic frameshift [Cuniasse et al. (1989) Biochemistry 28, 2018-2026], a left-handed loop was observed, whereas in this study a right-handed loop is found for the first time in solution. The deoxyribosylurea residue lies in the minor groove and can form both an intra- and an interstrand hydrogen bond.

Solution structure as a function of pH of two central mismatches, C . T and C . C, in the 29 to 39 K-ras gene sequence, by nuclear magnetic resonance and molecular dynamics.

The DNA duplexes 5' d(GCCACCAGCTC) x d(GAGCTXGTGGC), where the base X is either cytosine or thymine, have been studied by one and two-dimensional nuclear magnetic resonance, energy minimization and molecular dynamics. The sequence studied corresponds to the region 29 to 39 of the K-ras gene and is a hot spot for mutations. The results show that both duplexes adopt a globally B-DNA-type structure. For the C x C mismatch, we observe a structural change as a function of pH with an apparent pK of 6.95. The neutral species has only one hydrogen bond between the two bases but shows two families of wobble structures where one base or the other is displaced in the major groove. The protonated species has two hydrogen bonds and two structures but of unequal populations. In both systems, the sugar puckers remain predominantly C2'-endo and no significant changes in the backbone structure are observed. The neutral C . T mismatch is stabilized by two hydrogen bonds but, surprisingly, it can also be protonated, although the apparent pK is much lower, 5.65. In this case, protonation does not result in an additional hydrogen bond but must be due to better base-stacking interactions for C+ x T. The NMR data show that the environment of the T imino proton is very similar for C x T and C+ x T, although the hydrogen bond acceptor would be expected to be a nitrogen atom in the former case and an oxygen atom in the latter. We propose that for both structures there is an intervening water molecule which in addition reduces backbone strain. We have also measured the fluctuations during molecular dynamics runs in these mismatches. All are greater than for Watson-Crick base-pairs and the C x C mismatch shows very pronounced mobility.

Structure of oligonucleotides with either a strand break or a bulged nucleotide.

We report NMR and molecular modelling studies on a DNA duplex structure which is composed of three oligonucleotides and mimics a strand break. Although it retains a B form conformation our model suggests that it is kinked at the strand break. In the same sequence with an extra bulged adenosine at the centre for the major species this residue is stacked in the helix and a kink is observed in the model.

Methoxyamine-induced mutagenesis of nucleic acids. A proton NMR study of oligonucleotides containing N4-methoxycytosine paired with adenine or guanine.

We report the solution structure of two heptanucleotides each containing a central N4-methoxycytosine, in one case paired with adenine on the opposite strand and the other with guanine. For the N4-methoxycytosine adenine pair, only the imino form of the N4-methoxycytosine residue is observed and base pairing is in Watson-Crick geometry. However, rotation of the methoxy group about the N-OCH3 bond is not constrained to a particular orientation although it must be anti to the N3 of N4-methoxycytosine. The slow exchange on a proton NMR time scale between the single strand and double strand forms is attributed to the strong preference of the cis conformation of the OCH3 group in the single strand, which inhibits base pair formation. For the N4-methoxycytosine that is base paired with guanine, we observe an amino form in Watson-Crick geometry in slow exchange with a base paired imino form in wobble geometry. The amino form is predominant at low temperature whereas the imino form predominates above 313 K. We have measured the exchange rate between the two forms at 303 K and observed a value of approximately 1 S-1. The relative ratio of amino and imino forms of N4-methoxycytosine is influenced by both the base that is in front and the temperature. Our results explain the preferential replacement of dTTP by N4-methoxycytosine in primer elongation.

Solution structure of two mismatches G.G and I.I in the K-ras gene context by nuclear magnetic resonance and molecular dynamics.

Two mismatches, G.G and I.I, have been incorporated at the central position of 5'd-(GCCACXAGCTC).d(GAGCTXGTGGC) in order to carry out NMR and molecular dynamics studies. These duplexes constitute the sequence 29-39 of the K-ras gene coding for the glycine 12, a hot spot for mutation. The NMR spectra show that the duplexes are not greatly distorted by the introduction of the mismatches and their global conformation is that of a canonical B-form double helix. For the duplex containing the G.G mismatch, we propose for the major species, a type of pairing involving one hydrogen bond between the imino group of one central guanine and the carbonyl group of the opposite guanine. Both bases are in an anti conformation. Two conformations, with the same donor and acceptor pattern can coexist, one is obtained from the other by a 180 degrees rotation about the pseudodyadic axis. Exchange between the two forms is observed by NMR at low temperature. A minor species involving hydrogen bonding between the guanine amino group and the carbonyl group of the guanine on the opposite strand may also exist as shown by the molecular dynamics calculations. For the I.I mismatch we observe the same major species, i.e., hydrogen bonding between an imino proton of one base and the carbonyl group of the base on the opposite strand with both bases in an anti conformation. Exchange between these two conformations is faster than for the G.G mismatch. Further, we observe that the I.I mismatch adopts a minor conformation, in which one or other of the bases is in the syn conformation.

Methoxyamine attack on cytosine produces ambivalent base pairing properties of the modified base.

We report the solution structure of two heptanucleotides each containing a central N4-methoxycytosine, in one case with adenine on the opposite strand and in the other with guanine. For the N4-methoxycytosine-adenine pair only the imino form of the N4-methoxycytosine residue is observed and base pairing is in Watson-Crick geometry. However, rotation of the methoxy group about the N-OCH3 bond is not constrained to a particular orientation although it must be anti to the N3 of N4-methoxycytosine. The slow exchange on a proton NMR time scale between the single strand and double strand forms is attributed to the strong preference of the syn conformation of the OCH3 group in the single strand which inhibits base pair formation. For N4-methoxycytosine base paired with guanosine we observe the N4-methoxycytosine base in the amino form in Watson-Crick geometry and a slow exchange of this species with an imino form base paired in wobble geometry. The amino form is predominant at low temperature whereas the imino form predominates above 40 degrees C. Our results point to preferential replacement of dTTP by N4-methoxycytosine in primer elongation.

Solution structure of two mismatches A.A and T.T in the K-ras gene context by nuclear magnetic resonance and molecular dynamics.

Two mismatches, one homopurine (A.A) and the other homopyrimidine (T.T), have been incorporated at the central position N of: 5'd(GCCACNAGCTC).d(GAGCTNGTGGC) in order to study nuclear magnetic resonance spectra and molecular dynamics. These duplexes constitute the sequence 29-39 of the K-ras gene coding for Gly12, a hot spot for mutation. The NMR spectra show that the duplexes are not greatly distorted by the introduction of the mismatches and their global conformation is that of a canonical B-form double helix. For both systems, no structural change is observed in the pH range 4.7-9. For the duplex containing the homopurine A.A mismatch, we propose a type of pairing involving one hydrogen bond between the amino group of one central adenine and the nitrogen N1 of the opposite adenine. For the duplex containing the mispaired T.T bases, NMR spectra recorded in H2O at 282 K indicate that these central bases are engaged in wobble pairing, involving two imino-carbonyl hydrogen bonds. For both systems two conformations with the same donor and acceptor pattern can coexist, one being obtained from the other by a 180 degrees rotation about the pseudodyadic axis. Exchange between the two forms is observed by NMR at low temperature for the T.T mispair and also inferred from NMR measurements on the A.A system. The presence of this exchange and its pathway has been investigated by molecular dynamics calculations on both systems. Distance restrained and unrestrained molecular dynamics are in very good agreement with the NMR data. The average structure for either mispair shows only small conformational change from normal B DNA. For each, a systematic pathway is observed for exchange between the two conformations.

Solution structure of an oncogenic DNA duplex, the K-ras gene and the sequence containing a central C.A or A.G mismatch as a function of pH: nuclear magnetic resonance and molecular dynamics studies.

The DNA duplex 5' d(GCCACCAGCTC)-d(GAGCTGGTGGC) corresponds to the sequence 29 to 39 of the K-ras gene, which contains a hot spot for mutations. This has been studied by one and two-dimensional nuclear magnetic resonance, energy minimization and molecular dynamics. The results show that it adopts a globally B-DNA type structure. We have introduced, at the central base-pair, the mismatches C.A and A.G. The mismatch position is that of the first base of the Gly12 codon, the hot spot. For the C.A mismatch we observe a structural change as a function of pH with an apparent pKa of 7.2. At low pH, the mismatch pair adopts a structure close to a classic wobble conformation with the cytidine residue displaced into the major groove. It is stabilised by two hydrogen bonds in which the adenosine residue is protonated and the cytidine residue has a significant C3'-endo population. At high pH, the mispair structure is in equilibrium between wobble and reverse wobble conformations. Similar studies are reported on the A.G mismatch, which also undergoes a transition as a function of pH. 31P spectra have been recorded on all systems and as a function of pH. No evidence for BII phosphodiester backbone conformations was found. The NMR results are well corroborated by molecular dynamics calculations performed with or without distance constraints. The dynamics at the mismatch sites have been examined. Although the overall structures are close to B-DNA, helical parameters fluctuate differently at these sites. Different hydrogen bonding alternatives in dynamic equilibrium that can involve three-centred hydrogen bonds are observed.

Helical parameters, fluctuations, alternative hydrogen bonding, and bending in oligonucleotides containing a mistmatched base-pair by NOESY distance restrained and distance free molecular dynamics.

We have analysed and compared the molecular structures and dynamics of the DNA duplex, that corresponds to the sequence 29 to 39 of the K-ras gene, where the central base-pair is the normal C.G base or a mismatch base-pair C.A, C.A+, A.G and A+.G. Molecular dynamics runs of 100 picoseconds without or with distance restraints derived from NOE measurements are analysed and compared in all cases. (1) The EMBO convention of helical parameters for nucleic acids is extended to account for the construction and the description of any DNA mismatched base-pair. (2) Both types of MD runs reproduce very well all NMR data, except the H8 H2 inter-residue distances where agreement is not as good. (3) Average parameter values and fluctuations are in good agreement with results derived from persistence length and torsion modulus measurements. (4) Our molecular dynamics suggest the presence, in certain cases, of three-centred hydrogen bonds, which should be viewed as an equilibrium between hydrogen-bonding alternatives. In the case of the C.A mismatch, we observe the correlation between transient DNA bending and possible hydrogen bonding between a base and its 5' neighbour on the opposite strand in the sequence CCA. (5) These molecular dynamics analyses and observations provide a coherent view for the results obtained from recent DNA crystal structures and for results derived from other techniques such as gel electrophoresis at C.A steps.

Structures of oligonucleotides containing 5-(methoxymethyl)-2'-deoxyuridine determined by NMR spectroscopy.

Base pairing of 5-(methoxymethyl)-2'-deoxyuridine (MMdU) opposite either adenine or guanine in a seven base pair oligonucleotide duplex has been studied by NMR spectroscopy. When paired with A, we observe that the MMdU.A base pair adopts Watson-Crick geometry. The methoxymethyl substituent is not held in a fixed conformation and may rotate around the C5-CH2 and CH2-O bonds. Examination of the potential energy as a function of rotation around these bonds indicates the presence of four low energy conformations. No hydrogen bonding is indicated for the methoxymethyl substituent, and the four potential minima result from reduced steric clash. For the MMdU.G base pair, the two bases adopt a wobble geometry which does not change with increasing solvent pH. Similarly, we find four low energy conformations for the methoxymethyl substituent in the major groove of the DNA helix.

Structures of base pairs with 5-(hydroxymethyl)-2'-deoxyuridine in DNA determined by NMR spectroscopy.

Base pairs with 5-(hydroxymethyl)-2'-deoxyuridine (HMdU) opposite either adenine or guanine in a seven-base oligonucleotide duplex have been studied by NMR spectroscopy. When paired with A, the HMdU-A base pair is in Watson-Crick geometry. The hydroxymethyl group maintains a fixed orientation in which the oxygen is on the 5' side of the base. The energy-minimized structure indicates the presence of a hydrogen bond between the hydroxymethyl group and the N7 of the 5' guanine residue. When paired with guanine, HMdU-G is in a wobble configuration at low pH. The hydroxymethyl group is on the 3' side of the base, positioned to form an intramolecular hydrogen bond with its own O4 carbonyl. With increasing pH, HMdU-G is observed to ionize with an apparent pK value of 9.7. The high-pH structure is in a Watson-Crick configuration, with the HMdU residue in a position similar to that observed for HMdU-A. It is proposed that interresidue hydrogen bonding of the HMdU residue may stabilize aberrant base-pair configurations.

Base-pair induced shifts in the tautomeric equilibrium of a modified DNA base.

We have examined the base-pairing properties of N4-methoxycytosine (mo4C), a mutagenic base analog, in DNA by nuclear magnetic resonance spectroscopy. Unlike standard bases, the tautomeric equilibrium of mo4C could be strongly influenced by base-pair formation. Paired with A, mo4C is found predominantly in the imino configuration in Watson-Crick geometry. However, when paired with G, two structurally distinct configurations are observed in equilibrium with one another. In one configuration, mo4C is in the amino form paired with G in Watson-Crick geometry. In the second species, mo4C is in the imino configuration paired with G in a wobble geometry. This is the first demonstration of basepair induced tautomeric shifts in DNA and supports the hypothesis that rare tautomeric forms may be involved in mutagenesis.

Solution conformation of an oligonucleotide containing a urea deoxyribose residue in front of a thymine.

Urea residues are produced by ionizing radiation on thymine residues in DNA. We have studied an oligodeoxynucleotide containing a thymine opposite the urea residue, by one and two dimensional NMR spectroscopy. The urea deoxyribose exists as two isomers with respect to the orientation about the peptide bond. For the trans isomer we find that the thymine and urea site are positioned within the helix and are probably hydrogen bonded. The oligonucleotide adopts a globally B form structure although conformational changes are observed around the mismatch site. A minor species is observed, in which the urea deoxyribose and the opposite base adopt an extrahelical position and this corresponds to the isomer cis for the peptide bond.