NMR investigation of DNA primer-template models: structural insights into dislocation mutagenesis in DNA replication

NMR proton chemical shift prediction of C·C mismatches in B-DNA

Accurate prediction of DNA chemical shifts facilitates resonance assignment and allows recognition of different conformational features. Based on the nearest neighbor model and base pair replacement approach, we have determined a set of triplet chemical shift values and correction factors for predicting the proton chemical shifts of B-DNA containing an internal C·C mismatch. Our results provide a reliable chemical shift prediction with an accuracy of 0.07 ppm for non-labile protons and 0.09 ppm for labile protons. In addition, we have also shown that the correction factors for C·C mismatches can be used interchangeably with those for T·T mismatches. As a result, we have generalized a set of correction factors for predicting the flanking residue chemical shifts of pyrimidine·pyrimidine mismatches.

Cytosine unstacking and strand slippage at an insertion-deletion mutation sequence in an overhang-containing DNA duplex

Base unstacking in template strands, when accompanied by strand slippage, can result in deletion mutations during strand extension by nucleic acid polymerases. In a GCCC mutation hot-spot sequence, which was previously identified to have a 50% probability of causing such mutations during DNA replication by a Y-family polymerase, a single-base deletion mutation could result from such unstacking of any one of its three template cytosines. In this study, the intrinsic energetic differences in unstacking among these three cytosines in a solvated DNA duplex overhang model were examined using umbrella sampling molecular dynamics simulations. The free energy profiles obtained show that cytosine unstacking grows progressively more unfavorable as one moves inside the duplex from the 5′-end of the overhang template strand. Spontaneous strand slippage occurs in response to such base unstacking in the direction of both the major and minor grooves for all three cytosines. Unrestrained simulations run from three distinct strand-slipped states and one non-strand-slipped state suggest that a more duplexlike environment can help stabilize strand slippage. The possible underlying reasons and biological implications of these observations are discussed in the context of nucleic acid replication active site dynamics.

Analyzing the relationship between single base flipping and strand slippage near DNA duplex termini

Insertion-deletion (indel) mutations are caused by strand slippage between pairing primer and template strands during nucleic acid strand extension. A possible causative factor for such strand slippage is base flipping in the primer strand or template strand, for insertion or deletion mutations, respectively. A simple mechanistic description is that the "hole" in the nucleic acid duplex left behind by a flipping base is occupied by a neighboring base on the same strand, resulting in slippage with respect to its paired strand. The extent of single base flipping required for occupation of its former place in the double helix by a neighboring base is not fully understood. The present study uses restrained molecular dynamics (MD) simulations along a pseudohedihedral base flipping parameter to construct two-dimensional free energy profiles along base flipping and strand slippage geometric parameters. These profiles, generated for both cytosine and guanine single base flipping in a short repetitive indel mutation hot-spot DNA sequence, illustrate the extent of single base flipping that can allow strand slippage by one base position. Relatively minor base flipping into both the major and minor grooves can result in strand slippage. Deconstruction of the collective variable strand slippage geometric parameter into its component distances illustrates the details of how strand slippage can accompany base flipping. The trans Watson-Crick:sugar edge interaction that stabilizes cytosine flipping in this hot-spot sequence is also characterized energetically. The impact of these results on understanding sequence dependence of indel errors in nucleic acid strand extension is discussed, along with a suggestion for future studies that can generalize the present findings to all nearest-neighbor sequence contexts.

NMR proton chemical shift prediction of T·T mismatches in B-DNA duplexes

A proton chemical shift prediction scheme for B-DNA duplexes containing a T·T mismatch has been established. The scheme employs a set of T·T mismatch triplet chemical shift values, 5'- and 3'-correction factors extracted from reference sequences, and also the B-DNA chemical shift values predicted by Altona et al. The prediction scheme was tested by eight B-DNA duplexes containing T·T mismatches. Based on 560 sets of predicted and experimental proton chemical shift values, the overall prediction accuracy for non-labile protons was determined to be 0.07 ppm with an excellent correlation coefficient of 0.9996. In addition, the prediction accuracy for 96 sets of labile protons was found to be 0.22 ppm with a correlation coefficient of 0.9961. The prediction scheme developed herein can facilitate resonance assignments of B-DNA duplexes containing T·T mismatches and be generalized for the chemical shift prediction of other DNA mismatches. Our chemical shift data will also be useful for establishing structure-chemical shift information in B-DNA containing mismatches.

Partial base flipping is sufficient for strand slippage near DNA duplex termini

Strand slippage is a structural mechanism by which insertion-deletion (indel) mutations are introduced during replication by polymerases. Three-dimensional atomic-resolution structural pathways are still not known for the decades-old template slippage description. The dynamic nature of the process and the higher energy intermediates involved increase the difficulty of studying these processes experimentally. In the present study, restrained and unrestrained molecular dynamics simulations, carried out using multiple nucleic acid force fields, are used to demonstrate that partial base-flipping can be sufficient for strand slippage at DNA duplex termini. Such strand slippage can occur in either strand, i.e. near either the 3' or the 5' terminus of a DNA strand, which suggests that similar structural flipping mechanisms can cause both primer and template slippage. In the repetitive mutation hot-spot sequence studied, non-canonical base-pairing with exposed DNA groove atoms of a neighboring G:C base-pair stabilizes a partially flipped state of the cytosine. For its base-pair partner guanine, a similar partially flipped metastable intermediate was not detected, and the propensity for sustained slippage was also found to be lower. This illustrates that a relatively small metastable DNA structural distortion in polymerase active sites could allow single base insertion or deletion mutations to occur, and stringent DNA groove molecular recognition may be required to maintain intrinsic DNA polymerase fidelity. The implications of a close relationship between base-pair dissociation, base unstacking, and strand slippage are discussed in the context of sequence dependence of indel mutations.

Effect of an abasic site on strand slippage in DNA primer-templates

An abasic site is the most common lesion in DNA. It is also an intermediate product formed during base excision repair. Previously, we demonstrated that strand slippage can occur in primer-template model systems containing any kind of natural templating bases, suggesting deletion and expansion errors are possible in any kind of sequences during DNA replication. In this study, nuclear magnetic resonance spectroscopic investigations have been performed to study the intrinsic effect of a templating abasic residue on strand slippage in primer-template models. A DNA hairpin model system containing an abasic site and a 5'-overhang region was used to mimic the situation that a dNTP has just been incorporated opposite the abasic site. Our results show that, after dNTP incorporation, strand slippage occurs regardless of the type of terminal base pair formed. Compared to natural templating bases, abasic sites possess a higher slippage propensity, implicating a higher chance of expansion or deletion errors during DNA replication.

Sequence context effect on strand slippage in natural DNA primer-templates

Strand slippage has been found to occur in primer-templates containing a templating thymine, cytosine, and guanine, leading to the formation of misaligned structures with a single-nucleotide bulge. If remained in the active site of low-fidelity polymerases during DNA replication, these misaligned structures can ultimately bring about deletion mutations. In this study, we performed NMR investigations on primer-template models containing a templating adenine. Similar to our previous results on guanine, adenine templates are also less prone to strand slippage than pyrimidine templates. Misalignment occurs only in primer-templates that form a terminal C·G or G·C base pair. Together with our previous findings on thymine, cytosine, and guanine templates, the present study reveals strand slippage can occur in any kind of natural templating bases during DNA replication, providing insights into the origin of mutation hotspots in natural DNA sequences. In addition to the type of incoming base upon misincorporation, the propensity of strand slippage in primer-templates depends also on the type of templating base, its upstream and downstream bases.

Preferential base pairing modes of T·T mismatches

It has long been recognized that T·T mismatches can adopt two different modes of exchangeable wobble base pairs in which no preferential pairing mode has been observed. In this study, we have performed a systematic nuclear magnetic resonance (NMR) investigation to study the sequence context effect on the pairing modes of T·T mismatches. Our results reveal for the first time that preferential pairing mode does exist in T·T mismatches with specific type of flanking base pairs.

Impact of DNA physical properties on local sequence bias of human mutation

In selectively neutral regions of the human genome, nucleotide substitutions do not occur at random with respect to the local DNA sequence neighborhood. However, apart from the hypermutability of methylated CpG dinucleotides, which can explain the overrepresentation of nucleotide transitions in this context, the sequence-specific factors underlying point mutation bias remain largely to be determined, both in nature and in quantitative impact. One hypothesis suggests that the physical characteristics of a DNA context could have a modulating effect on its mutability, adjusting the impact of damage or the efficiency of repair. Here, we report a genome-wide computational test of this hypothesis, in which we utilize a constrained set of human non-CpG SNPs as the source of selectively neutral germline mutations. Interestingly, we observe that the quantitative context-dependencies of some substitution types display significant associations to measures of local structural topography and helix stability in DNA. Most prominently, we find that the local sequence bias of transition mutations is significantly associated with the sequence-dependent level of helix instability imposed by the potentially underlying DNA mismatches. The results of our work indicate the extent to which DNA physical properties could have shaped the recent point mutational spectrum in the human genome.

NMR investigation of DNA primer-template models: guanine templates are less prone to strand slippage upon misincorporation

Misaligned structures can result from strand slippage during DNA replication and, if not repaired, would lead to mutations. Previously, we showed that strand slippage can occur upon misincorporation of a dNTP opposite thymine and cytosine templates, resulting in a misaligned structure with a T- or C-bulge. The formation propensity for misaligned structures was found to depend on the type of terminal base pair. In this study, we performed NMR investigations on primer-template models containing a guanine template. Our results reveal guanine templates are less prone to strand slippage than pyrimidine templates. Misalignment was found to occur only in 5'-CG templates with a downstream purine. In addition to the significance of terminal base pair and upstream nucleotide, the present study reveals the importance of the templating base and its downstream nucleotide, which also determine the propensity of strand slippage in primer-templates.

Effect of hyperoxidized guanine on DNA primer-template structures: spiroiminodihydantoin leads to strand slippage

Oxidation of guanine in DNA can lead to mutagenic lesions such as 7-hydro-8-oxoguanine (oG). Upon further oxidation, a more mutagenic lesion, spirominodihydantoin (Sp), can occur. In this study, nuclear magnetic resonance (NMR) investigations were performed to determine the structural features of DNA primer-template models with 5'-GG, 5'-G(oG), 5'-G(Sp) and 5'-T(Sp) templates, that mimic the situation in which the downstream G of the template has been oxidized to oG or hyperoxidized to Sp. Our results show that misalignment occurs only in the 5'-G(Sp) and 5'-T(Sp) templates, providing structural insights into the observed differences in mutagenicity of Sp and oG during DNA replication.

Nuclear magnetic resonance investigation of primer--template models: formation of a pyrimidine bulge upon misincorporation

Our previous studies have shown that misaligned structures can occur upon misincorporation of a dNTP opposite thymine templates. The formation of misaligned structures during DNA replication, if not repaired properly, can be bypassed and extended by low-fidelity polymerases and ultimately lead to mutations. In this study, the base pair structures at the replicating sites of a set of primer-template models which mimic the situation upon misincorporation of a dNTP opposite cytosine templates have been determined. High-resolution NMR structural results show that misaligned structures with a C-bulge can be formed upon incorporation of dCTP, dTTP, and dATP opposite 5'-GC, 5'-AC, and 5'-TC templates, respectively. The stabilities of misaligned structures depend on the types of terminal base pairs at the replicating sites. Together with the structural findings in thymine templates, we conclude that terminal G.C and C.G base pairs always contribute a larger stabilizing effect to the misaligned structures containing a pyrimidine bulge than terminal A.T and T.A base pairs. Misalignment and thus deletion mutation are more likely to occur if misincorporation of a nucleotide opposite a pyrimidine template can cause template slippage to form a terminal G.C or C.G base pair. Although misalignment also occurs when the newly formed terminal base pair is an A.T base pair or a T.A base pair, both misaligned and mismatched conformers coexist, which can lead to deletion and substitution mutations, respectively.

A mutation in KIR3DS1 that results in truncation and lack of cell surface expression

The KIR gene cluster exhibits a high degree of polymorphism in terms of gene content as well as allelic polymorphism, and data suggest that it is evolving rapidly. The KIR3DL1 locus is one of the most polymorphic loci within this cluster and is unique in that it encodes an activating receptor KIR3DS1, as well as multiple inhibitory KIR3DL1 allotypes. Because KIR3DS1 has been implicated in a number of diseases, we tested for the presence of KIR3DS1 variants that might affect its expression and activating capacity. Preliminary FACS analysis indicated that indeed some individuals with the KIR3DS1 allele showed no cell surface expression of the molecule. Sequencing analysis identified a variant with a complex deletion/substitution mutation in exon 4 (which encodes the D1 extracellular domain), resulting in a premature stop codon. We subsequently genotyped 3,960 unrelated individuals and determined the frequencies of this allele across geographically distinct world populations. The data indicate that the null KIR3DS1 allele is uncommon, arose on a single haplotype, and spread across geographically distinct populations.

Proton chemical shift prediction of A.A mismatches in B-DNA duplexes

A proton chemical shift prediction method has been developed for double helical DNAs containing A.A mismatches. This method makes use of the chemical shift prediction scheme for normal B-DNA duplexes developed by Altona and co-workers and a set of A.A mismatch triplet chemical shift values and corrections factors extracted from reference sequences. The triplet values are used for predicting chemical shifts of A.A mismatches whereas the normal B-DNA chemical shifts and correction factors are used for the flanking residues of A.A mismatches. Both 5'- and 3'-correction factors have been determined from the chemical shift differences upon replacing the A.A mismatch in a duplex with an A.T base pair. Based on 560 sets of predicted and experimental chemical shifts, the overall prediction accuracy for various types of protons has been determined to be 0.07 ppm with an excellent correlation coefficient of 0.9996.