Molecular-mechanical studies of the mismatched base analogs of d(CGCGAATTCGCG)2:d(CGTGAATTCGCG)2, d(CGAGAATTCGCG)2, d(CGCGAATTCACG)2, d(CGCGAATTCTCG)2, and d(CGCAGAATTCGCG).d(CGCGAATTCGCG)

Factors determining the deriving force of DNA formation: geometrical differences of base pairs, dehydration of bases, and the arginine assisting

The mechanism of the fidelity synthesis of DNA associated with the process of dGTP combination to the DNA template was explored. The exclusion of water molecules from the hydrated DNA bases can amplify the energy difference between the correct and incorrect base pairs, but the effect of the water molecules on the Gibbs free energy of formation is dependent on the binding sites for the water molecules. The water detachment from the incoming dNTP is not the only factor but the first step for the successful replication of DNA. The second step is the selection of the DNA polymerase on the DNA base pair through the comparison between the correct DNA base and the incorrect DNA base. The bonding of the Arg668 with the incoming dNTP can enlarge the Gibbs free energies of formation of the base pairs, especially the correct base pairs, thus increasing the driving force of DNA formation. When the DNA base of the primer terminus is correct, the extension of the guanine and the adenine is quicker than that of the cytosine and the thymine because of the hydrogen bonding fork formation of Arg668 with the minor groove of the primer terminus and the ring oxygen of the deoxyribose moiety of the incoming dNTP. Because of the geometry differences of the incorrect base pairs with the correct base pairs, the effect from the DNA polymerase is smaller on the incorrect base pair than on the correct base pair, and the extension of a mispair is slower than that of a correct base pair. This decreases the extension rate of the base pair and thus allows proofreading exonuclease activity to excise the incorrect base pair. Arg668 cannot prevent the extension of the GT mispair, as well as the GC correct base pair, and GA and GG mispairs. This may be attributed to the small geometry difference between the GT base pair and the correct AT base pair.

The role of Arg165 towards base flipping, base stabilization and catalysis in M.HhaI

Arg165 forms part of a previously identified base flipping motif in the bacterial DNA cytosine methyltransferase, M.HhaI. Replacement of Arg165 with Ala has no detectable effect on either DNA or AdoMet affinity, yet causes the base flipping and restacking transitions to be decreased approximately 16 and 190-fold respectively, thus confirming the importance of this motif. However, these kinetic changes cannot account for the mutant's observed 10(5)-fold decreased catalytic rate. The mutant enzyme/cognate DNA cocrystal structure (2.79 A resolution) shows the target cytosine to be positioned approximately 30 degrees into the major groove, which is consistent with a major groove pathway for nucleotide flipping. The pyrimidine-sugar chi angle is rotated to approximately +171 degrees, from a range of -95 degrees to -120 degrees in B DNA, and -77 degrees in the WT M.HhaI complex. Thus, Arg165 is important for maintaining the cytosine positioned for nucleophilic attack by Cys81. The cytosine sugar pucker is in the C2'-endo-C3'-exo (South conformation), in contrast to the previously reported C3'-endo (North conformation) described for the original 2.70 A resolution cocrystal structure of the WT M.HhaI/DNA complex. We determined a high resolution structure of the WT M.HhaI/DNA complex (1.96 A) to better determine the sugar pucker. This new structure is similar to the original, lower resolution WT M.HhaI complex, but shows that the sugar pucker is O4'-endo (East conformation), intermediate between the South and North conformers. In summary, Arg165 plays significant roles in base flipping, cytosine positioning, and catalysis. Furthermore, the previously proposed M.HhaI-mediated changes in sugar pucker may not be an important contributor to the base flipping mechanism. These results provide insights into the base flipping and catalytic mechanisms for bacterial and eukaryotic DNA methyltransferases.

Nucleobase pairing in expanded Watson-Crick-like genetic information systems

To guide the design of alternative genetic systems, we measured melting temperatures of DNA duplexes containing matched and mismatched nucleobase pairs from natural and unnatural structures. The pairs were analyzed in terms of structural features, including nucleobase size, number of hydrogen bonds formed, the presence of uncompensated hydrogen bonding functional groups, the nature of the bond joining the nucleobase to the sugar, and nucleobase charge. The results suggest that stability of nucleobase pairs correlates with the number of H-bonds, size complementarity, the presence of uncompensated functional groups, and the presence of charge on a nucleobase. Each of these properties appear to be more significant than the nature of the glycosidic bond and sequence context. The results provide guidelines for constructing stable Watson-Crick like nucleobase pairs with unnatural nucleobases. The experiments also demonstrate that expanded genetic systems can be constructed using size complementary nucleobase pairs that contain three hydrogen bonds.

New insights into the structure of abasic DNA from molecular dynamics simulations

Abasic (AP) sites constitute a common form of DNA damage, arising from the spontaneous or enzymatic breakage of the N-glycosyl bond and the loss of a nucleotide base. To examine the effects of such damage on DNA structure, especially in the vicinity of the abasic sugar, four 1.5 ns molecular dynamics simulations of double-helical DNA dodecamers with and without a single abasic (tetrahydrofuran, X) lesion in a 5'-d(CXT) context have been performed and analyzed. The results indicate that the abasic site does not maintain a hole or gap in the DNA, but instead perturbs the canonical structure and induces additional flexibility close to the abasic site. In the apurinic simulations (i.e., when a pyrimidine is opposite the AP site), the abasic sugar flipped in and out of the minor groove, and the gap was water filled, except during the occurrence of a novel non-Watson-Crick C-T base pair across the abasic site. The apyrimidinic gap was not penetrated by water until the abasic sugar flipped out and remained extrahelical. Both AP helices showed kinks of 20-30 degrees at the abasic site. The Watson-Crick hydrogen bonds are more transient throughout the DNA double helices containing an abasic site. The abasic sugar displayed an unusually broad range of sugar puckers centered around the northern pucker. The increased motion of the bases and backbone near the abasic site appear to correlate with sequence-dependent helical stability. The data indicate that abasic DNA contorts more easily and in specific ways relative to unmodified DNA, an aspect likely to be important in abasic site recognition and hydrolysis.

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.

Synthesis and characterization of 8-methoxy-2'- deoxyadenosine-containing oligonucleotides to probe the syn glycosidic conformation of 2'-deoxyadenosine within DNA

The synthesis of 8-methoxy-2'-deoxyadenosine (moA) protected at N6 as an N,N-dimethylformamidine derivative and incorporation of the modified nucleoside into oligodeoxynucleotides via the phosphoramidite method are described. UV thermal denaturation studies were conducted on duplexes containing moA:G, moA:C and moA:T base pairs to determine the thermodynamic stability of duplexes containing moA relative to their adenosine (A)-containing counterparts. In the case of moA:G base pairs the effect of moA substitution is sequence dependent. In A:G mismatch-containing sequences, which have been shown by structural characterization to have a syn conformational preference at the glycosidic bond of A, moA substitution results in stabilization of the duplex. In contrast, in sequences where the A in the A:G mismatch has been shown to prefer the anti conformation moA substitution is destabilizing to the duplex. Thus moA may be a useful probe for investigating the conformational preferences of the N-glycosidic bond of adenosine within DNA. In addition, moA nucleoside is more resistant to acid-catalyzed depurination than previously described 8-bromo-2'-deoxyadenosine, allowing for facile incorporation into oligonucleotides via automated solid phase DNA synthesis.

Molecular dynamics simulations suggest that the Eco RI kink is an example of molecular strain

The energy surface in the vicinity of the "Eco RI kink" was investigated by conducting both in vacuo molecular dynamics simulations as well as a simulation with explicit solvent. The in vacuo simulations used the "all atom" AMBER 3.0 force field with a distant dependent dielectric function and "hydrated" counter ions while the simulation with explicit solvent used the AMBER 4.0 force field, fully charged phosphates and counter ions and a dielectric constant of 1.0. The thrust of the simulations was to discriminate between two models of the energy surface of the deformed DNA as found in the recognition complex with Eco RI endonuclease. In the intrinsic model, the kinked DNA is a local minimum of the energy surface intrinsic to the DNA itself while in the strained model there is no significant energy barrier separating kinked and regular B-DNA. The two models have significant implications for theories of indirect recognition of DNA based on sequence-dependent deformability. The simulations suggest that the Eco RI-kinked structure is an example of molecular strain because it is not near a minimum of any of the potential energy functions examined. The simulations leave the question of an energy barrier somewhat open and raise the possibility that the Eco RI kink is at (or near) a point of dynamic instability of the energy surface (either a true maximum or a saddle point).

N.m.r. determination of the solution conformation and dynamics of the A.G mismatch in the d(CGCAAATTGGCG)2 dodecamer

A.G base-paired mismatches that occur during replication are among the most difficult to detect by repair enzymes. Such purine.purine mispairs can exist in two conformations, one of which is stabilized by protons [Gao & Patel (1988) J. Am. Chem. Soc. 110, 5178-5182]. We have undertaken a 1H-n.m.r. and 31P-n.m.r. study of the mismatched dodecamer d(CGCAAATTGGCG)2 as a function of both temperature and pH to determine the conformational features of the A.G mismatch. At pH greater than 7 the mispaired bases are each in the anti conformation and are stacked in the B-like helix. As the pH is decreased, a second conformation becomes populated (apparent pKa approx. 5.9) with concomitant changes in the chemical shifts of protons of the mispaired bases and their nearest neighbours. Data from two-dimensional nuclear-Overhauser-enhancement spectroscopy show unequivocally that, at low pH, the dominant conformation is one in which the mismatched G residues are in the syn conformation and are hydrogen-bonded to the A residues that remain in the anti conformation. Residues not adjacent to the A.G sites are almost unaffected by the transition or the mispairing, suggesting considerable local flexibility of the unconstrained duplexes. Despite the bulging of the mispaired bases, the conformation of the A(anti).G(anti) duplex is very similar to the native dodecamer, whereas the AH+(anti).G(syn) duplex shows a greater variation in the backbone conformation at the mismatched site. According to the chemical shifts, the duplex retains twofold symmetry in solution. The equilibrium between the syn and anti conformations of G9/G21 is strongly dependent on pH, but only weakly dependent on temperature (delta H approx. 16 kJ.mol-1). The first-order rate constant for the transition is approx. 9 s-1 at 283 K and approx. 60 s-1 at 298 K, with an activation enthalpy of approx. 100 kJ.mol-1. The stabilization of the A(anti).G(syn) conformation by protons is consistent with models invoking N1 protonation of adenine. Using the derived glycosidic torsion angles we have used restrained molecular dynamics to build models of the neutral and protonated d(CGCAAATTGGCG)2 oligomers. The results confirm that the A(anti).G(anti) and AH+(anti).G(syn) conformations are favoured at high pH and low pH respectively, in accord with n.m.r. and single-crystal X-ray data.

Theoretical and experimental measures of DNA helix stability and their relation to sequence specific repair of O6-ethylguanine lesions

Recent work (Breslauer et al. (1986) Proc. Natl. Acad. Sci. (U.S.A.), 83, 3746) has provided a method for calculating empirical thermodynamic quantities for helix to coil transitions from the base sequence of any oligomer. It is shown in this work that the DNA helix binding energy, calculated with the AMBER force field, for 9-mers of the type 5'-GGGXGeYGGG-3', where X and Y are any base and the central Ge is O6-ethylguanine, correlates well with the empirical delta G for helix to strand transitions. The mutation spectrum of ethane methylsulfonate (EMS) in the lacI gene of Escherichia coli can be modeled using the calculated local binding energy but the empirical free energies, enthalpies and melting temperatures predict these levels of repair less well. The relation of the binding energy to the mutation spectrum can be somewhat improved by including entropic effects in a theoretical free energy of binding as given by delta G theoretical identical to delta E binding - T delta S.

Block-units method for conformational calculations of large nucleic acid chains. II. The two-hierarchical approach and its application to conformational arrangement of the unusual T psi C loop of rabbit tRNA(Val)

The two-level hierarchical methodology is suggested for conformational calculations of large fragments of nucleic acids. The method of the first level is intended for performing a fast screening of the conformational phase space. The high-level method may be used to refine structurally important conformations. The method of the first level is the block-units method, which has been developed specially for these purposes (see part I). It has been shown that the block-units method allows the satisfactory calculation of the structure parameters of the optimal conformations of polynucleotides. The results of the conformational rearrangement calculations of the T psi C loop of the tRNA(Phe) after modification of its sequence are represented.

Crystal structure and sequence-dependent conformation of the A.G mispaired oligonucleotide d(CGCAAGCTGGCG)

The crystal structure of the dodecanucleotide d(CGCAAGCTGGCG) has been determined to a resolution of 2.5 A and refined to an R factor of 19.3% for 1710 reflections. The sequence crystallizes as a B-type double helix, with two G(anti).A(syn) base pairs. These are stabilized by three-center hydrogen bonds to pyrimidines that induce perturbations in base-pair geometry. The central AGCT region of the helix has a wide (greater than 6 A) minor groove.

Crystal structure and stability of a DNA duplex containing A(anti).G(syn) base-pairs

The synthetic dodecanucleotide d(CGCAAATTGGCG) has been analysed by single-crystal X-ray diffraction techniques and the structure refined to R = 0.16 and 2.25 A resolution, with the location of 94 solvent molecules. The sequence crystallizes as a full turn of a B-DNA helix with ten Watson-Crick base-pairs and two adenine-guanine mispairs. The analysis clearly shows that the mismatches are of the form A(anti).G(syn). Thermal denaturation studies indicate that the stability of the duplex is strongly pH dependent, with a maximum at pH 5.0, suggesting that the base-pair is stabilized by protonation. Three different arrangements have been observed for base-pairs between guanine and adenine and it is likely that A.G mismatch conformation is strongly influenced by dipole-dipole interactions with adjacent base-pairs.

Mispair formation in DNA can involve rare tautomeric forms in the template

The formation of pyridine-pyrimidine- and pyrimidine-pyrimidine base pairs after in vitro DNA replication with the large fragment of Escherichia coli DNA polymerase I indicates that Watson-Crick-like base pairing between pyrimidine bases can occur in the enzyme due to the presence of the rare tautomers of deoxycytidylate and thymidylate in the template strand. The implications to mispair formation in DNA, such as the difference between the structures of the mispairs during and after replication, are discussed and the possible action of mutagenic DNA protonating and deprotonating agents in vivo is considered.

Pyrimidine.pyrimidine base-pair mismatches in DNA. A nuclear magnetic resonance study of T.T pairing at neutral pH and C.C pairing at acidic pH in dodecanucleotide duplexes

Structural features of pyrimidine.pyrimidine mismatches in the interior of oligonucleotide duplexes have been investigated by high resolution two-dimensional proton nuclear magnetic resonance (n.m.r.) spectroscopy. These studies were conducted on the self-complementary d(C-G-C-T-A-G-C-T-T-G-C-G) duplex (designated T.T 12-mer) and the self-complementary d(C-G-C-C-A-G-C-T-C-G-C-G) duplex (designated C.C 12-mer) containing T.T and C.C pairs located at identical positions four base-pairs from either end of the duplex. Proton n.m.r. studies on the T.T 12-mer duplex were undertaken in the neutral pH range, while studies on the C.C 12-mer duplex were recorded at acidic pH. The proton spectra narrowed considerably on lowering the pH below neutrality for the C.C 12-mer duplex. Two-dimensional nuclear Overhauser enhancement spectroscopy (NOESY) data sets have been recorded on the T.T 12-mer and C.C 12-mer duplexes in high salt H2O and D2O solution. The magnitude of the NOE crosspeaks and the directionality of the NOE connectivities demonstrate that both duplexes are right-handed with all bases, including those at the mismatch site, adopting an anti configuration about the glycosidic bond. The observed base and sugar proton chemical shifts suggest structural similarities for the trinucleotide segments centered about the T.T and C.C mismatches. A NOE is detected between the resolved imino protons of T4 and T9 at the mismatch site, consistent with formation of a stacked "wobble" T4(anti).T9(anti) pair in the T.T 12-mer duplex. A comparison of the imino proton chemical shift and NOE data suggests that the imino-carbonyl hydrogen bonds in the wobble T.T mismatch are weaker than the corresponding imino-carbonyl hydrogen bonds in the wobble G.T mismatch. The 4-amino protons of C4 and C9 at the mismatch site in the C.C 12-mer duplex do not exhibit the pattern of hydrogen-bonded and exposed protons separated by approximately 1.5 parts per million characteristic of cytidine amino protons involved in Watson-Crick G.C pairing. The experimental data are insufficient to differentiate between wobble C(anti).C+(anti) and other pairing possibilities for the mismatch in the C.C 12-mer duplex at acidic pH.

A kinked model for the solution structure of DNA tridecamers with inserted adenosines: energy minimization and molecular dynamics

Structural modelling techniques using energy minimization and molecular dynamics have been employed to generate kinked models for the solution structure of two DNA tridecamer sequences containing inserted adenosines: d(CGCAGAATTCGCG)2 and d(CGCAGAGCTCGCG)2. These models are consistent with NMR studies of these sequences in solution. The overall shapes of the two models are similar, consisting of three B-DNA sections: two outer segments on the same side of the central portion, with the additional adenosines acting as wedges to kink the structure. An alternative scheme for the hydrogen bond pairing at the kink site is suggested as a way for the additional adenosines to be stabilized in the duplex.

Nuclear magnetic spectra of self complementary decanucleotides in solution; base sequence effect on the chemical shifts of nonexchangeable protons

The aim of this study was to attempt to determine the extent to which the chemical shifts of the nonexchangeable base protons of a DNA helix depend upon the base sequence. We measured the proton NMR spectra of twelve decadeoxynucleotides in order to carry out a "statistical" treatment. In the helices, the chemical shifts were found to be determined within +/- 0.04 ppm, largely by the nearest neighbor residues on the 5'-side, and to a smaller extent by the residue on the 3'-side. The theoretical chemical shift calculations reproduced very well the polymerization shifts measured for H2 protons of adenosines if the electrostatic field effect was taken into account. A fair agreement was also obtained for H8 protons of the adenosine and guanosine residues. However, theory underestimates the polarization effects of the base protons of cytidine. This discrepancy suggests that the conformation of this residue is different in the mononucleotides relative to double helices.

One- and two-dimensional 1H-NMR investigations of two 19-base-pair analogues of the tet operator

Two 19-base-pair oligodeoxynucleotides, analogues of one of the operators which specifically bind the repressor protein in the regulatory part of the transposon Tn10 tetracycline-resistance (tet) determinant, have been studied by 1H-NMR spectroscopy. The analogues contain a mismatch in the central base pair of the double helix (T.T or A.A). The imino protons have been assigned to the base pairs by one-dimensional NOE measurements, and the thermally induced transition from the duplex to the single strand has been followed. The cytidine amino resonances have been assigned by means of two-dimensional NOE spectroscopy in H2O. Two-dimensional phase-sensitive NOE and magnitude-correlated spectra have been recorded in 2H2O; all nonexchangeable protons, with the exception of some of H5', H5" protons, have been assigned. The NMR data made it possible to carry out a qualitative analysis of the structures of both oligodeoxynucleotides. The general structures close to B-DNA, show irregularities in the mismatch areas.

Helix geometry, hydration, and G.A mismatch in a B-DNA decamer

The DNA double helix is not a regular, featureless barberpole molecule. Different base sequences have their own special signature, in the way that they influence groove width, helical twist, bending, and mechanical rigidity or resistance to bending. These special features probably help other molecules such as repressors to read and recognize one base sequence in preference to another. Single crystal x-ray structure analysis is beginning to show us the various structures possible in the B-DNA family. The DNA decamer C-C-A-A-G-A-T-T-G-G appears to be a better model for mixed-sequence B-DNA than was the earlier C-G-C-G-A-A-T-T-C-G-C-G, which is more akin to regions of poly(dA).poly(dT). The G.A mismatch base pairs at the center of the decamer are in the anti-anti conformation about their bonds from base to sugar, in agreement with nuclear magnetic resonance evidence on this and other sequences, and in contrast to the anti-syn geometry reported for G.A pairs in C-G-C-G-A-A-T-T-A-G-C-G. The ordered spine of hydration seen earlier in the narrow-grooved dodecamer has its counterpart, in this wide-grooved decamer, in two strings of water molecules lining the walls of the minor groove, bridging from purine N3 or pyrimidine O2, to the following sugar O4'. The same strings of hydration are present in the phosphorothioate analog of G-C-G-C-G-C. Unlike the spine, which is broken up by the intrusion of amine groups at guanines, these water strings are found in general, mixed-sequence DNA because they can pass by unimpeded to either side of a guanine N2 amine. The spine and strings are perceived as two extremes of a general pattern of hydration of the minor groove, which probably is the dominant factor in making B-DNA the preferred form at high hydration.

Inosine.adenine base pairs in a B-DNA duplex

The structure of the synthetic deoxydodecamer d(C-G-C-I-A-A-T-T-A-G-C-G) has been determined by single crystal X-ray diffraction techniques at 2.5A resolution. The refinement converged with a crystallographic residual, R = 0.19 and the location of 64 solvent molecules. The sequence crystallises as a B-DNA helix with 10 Watson-Crick base-pairs (4 A.T. and 6 G.C) and 2 inosine.adenine (I.A) pairs. The present work shows that in the purine.purine base-pairs the adenine adopts syn orientation with respect to the furanose moiety while the inosine is in the trans (anti) orientation. Two hydrogen bonds link the I.A. base-pair, one between N-1(I) and N-7(A), the other between O-6(I) and N-6(A). This bulky purine.purine base-pair is incorporated in the double helix at two positions with little distortion of either local or global conformation. The pairing observed in this study is presented as a model for I.A base-pairs in RNA codon-anticodon interactions and may help explain the thermodynamic stability of inosine containing base-pairs. Conformational parameters and base stacking interactions are presented and where appropriate compared with those of the native compound, d(C-G-C-G-A-A-T-T-C-G-C-G) and with other studies of oligonucleotides containing purine.purine base-pairs.

Variation of nonexchangeable proton resonance chemical shifts as a probe of aberrant base pair formation in DNA

Variation of nonexchangeable proton resonance chemical shifts for deoxycytidine and deoxy-adenosine as a function of protonation and imino tautomer formation has been determined. Protonation induces downfield shifts of proton resonances whereas formation of the rare imino tautomer induces upfield shifts. Titration curves are constructed on the basis of spectrophotometrically determined pK values. Excellent fit is obtained between theoretical titration curves and experimental data, which indicates that chemical shifts of base protons may be used to quantitatively determine the relative concentrations of either rare imino tautomeric conformations or protonated base forms. These data may be utilized as an aid in the elucidation of the nature of hydrogen bonding between mismatched base pairs in DNA oligomers containing cytosine or adenine residues. These data, in conjunction with the oligonucleotide study of Patel et al. [Patel, D. J., Kozlowski, S.A., Ikuta, S., & Itakura, K. (1984) Biochemistry 23, 3218-3226], have been used to rigorously argue the existence of a "protonated" adenine residue in the A-C mismatch. This structure allows reconciliation of the NMR solution data with crystallographic data [Hunter, W.N., Brown, T., Anand, N.N., & Kennard, O. (1986) Nature (London) 320, 552-555], which support the protonated base pair.

Mismatches in DNA double strands: thermodynamic parameters and their correlation to repair efficiencies

The helix-coil transitions of the 16 octadecameric DNA duplexes dCGTCGTTTXACAACGTCG X dCGACGTTGTX1AAACGACG with A, T, G, and C for X and X1 were measured by UV-absorption. This sequence was taken from former studies of in vivo determination of efficiencies of mismatch repair (Kramer, Kramer, and Fritz (1984) Cell 38, 879-887). The thermodynamic parameters for double strand and mismatch formation have been obtained by evaluating the partition function of a stack model which allowed for loop formation. As a result the mismatches could be classified into wobble base pairs (T/G, G/G, C/A, A/A, A/G), open base pairs, i.e. permanent loops (T/T, C/T, T/C, C/C), and intermediate or weak base pairs (G/T, A/C, G/A). There is no correlation between Tm and the biological repair efficiency of X/X1. The structure classes, however, as described above show a close correlation: Open base pairs show the lowest repair efficiencies, whereas mismatches with high repair efficiency always belong to the structural class of wobble base pairs. Because of the palindromic nearest neighbors of the variation site X/X1, the influence of next-nearest neighbor interactions could be detected and be estimated to about 1 kJ/mol for one stack.

Multistep modeling (MSM) of biomolecular structure application to the A-G mispair in the B-DNA environment

A multistep modeling procedure has been evolved to study the structural changes introduced by lesions in DNA. We report here the change in the structure of regular B-DNA geometry due to the incorporation of Ganti-Aanti mispair in place of a regular G-C pair, preserving the helix continuity. The energetics of the structure so obtained is compared with the Ganti-Asyn configuration under similar constrained conditions. We present the methodology adopted and discuss the results.

Nuclear magnetic resonance and circular dichroism studies of a duplex--single-stranded hairpin loop equilibrium for the oligodeoxyribonucleotide sequence d(CGCGATTCGCG)

Nuclear magnetic resonance (NMR) and circular dichroism (CD) studies have been carried out with the oligodeoxyribonucleotide mismatch sequence, d(CGCGATTCGCG), 1. It has been found that 1 exists, in solution, as an equilibrium mixture of slowly interconverting, structured conformational isomers, 1a and 1b. On the basis of the concentration dependence of the 1a-1b equilibrium, the 1H NMR spectrum of the imino protons of the nucleotide bases, and the individual CD spectra of 1a and 1b, it is suggested that the two species correspond to a B-type DNA duplex and a single-stranded, hairpin-loop structure; the portion of the single-stranded species not involved in the loop appears to have a B-type DNA structure (on the basis of the CD measurements). To facilitate 1H NMR resonance assignments, the two possible des-methyl thymidine derivatives of 1 were synthesized; the effect of this substitution on the physical chemical properties of 1 was explored. The 1H NMR spectra of 1, as a function of temperature, showed that, under conditions wherein both species were present to a significant extent, the duplex form melted at a lower temperature than the single-stranded, hairpin loop structure.

Base-base mismatches. Thermodynamics of double helix formation for dCA3XA3G + dCT3YT3G (X, Y = A,C,G,T)

Thermodynamic parameters for double strand formation have been measured for the sixteen double helices of the sequence dCA3XA3G.dCT3YT3G, with each of the bases A, C, G and T at the positions labelled X and Y. The results are analyzed in terms of nearest-neighbors and are compared with thermodynamic parameters for RNA secondary structure. At room temperature the sequence (Formula: see text) is more stable than (Formula: see text) and is similar in stability to (Formula: see text) and (Formula: see text) are least stable. At higher temperatures the sequences containing a G.C base pair become more stable than those containing only A.T. All molecules containing mismatches are destabilized with respect to those with only Watson-Crick pairing, but there is a wide range of destabilization. At room temperature the most stable mismatches are those containing guanine (G.T, G.G, G.A); the least stable contain cytosine (C.A, C.C). At higher temperatures pyrimidine-pyrimidine mismatches become the least stable.