Structural differences in the NOE-derived structure of G-T mismatched DNA relative to normal DNA are correlated with differences in (13)C relaxation-based internal dynamics

DNA strand breaks and gaps target retroviral intasome binding and integration

Retrovirus integration into a host genome is essential for productive infections. The integration strand transfer reaction is catalyzed by a nucleoprotein complex (Intasome) containing the viral integrase (IN) and the reverse transcribed (RT) copy DNA (cDNA). Previous studies suggested that DNA target-site recognition limits intasome integration. Using single molecule Förster resonance energy transfer (smFRET), we show prototype foamy virus (PFV) intasomes specifically bind to DNA strand breaks and gaps. These break and gap DNA discontinuities mimic oxidative base excision repair (BER) lesion-processing intermediates that have been shown to affect retrovirus integration in vivo. The increased DNA binding events targeted strand transfer to the break/gap site without inducing substantial intasome conformational changes. The major oxidative BER substrate 8-oxo-guanine as well as a G/T mismatch or +T nucleotide insertion that typically introduce a bend or localized flexibility into the DNA, did not increase intasome binding or targeted integration. These results identify DNA breaks or gaps as modulators of dynamic intasome-target DNA interactions that encourage site-directed integration.

Backbone Conformational Equilibrium in Mismatched DNA Correlates with Enzyme Activity

T:G mismatches in mammals arise primarily from the deamination of methylated CpG sites or the incorporation of improper nucleotides. The process by which repair enzymes such as thymine DNA glycosylase (TDG) identify a canonical DNA base in the incorrect pairing context remains a mystery. However, the abundant contacts of the repair enzymes with the DNA backbone suggest a role for protein-phosphate interaction in the recognition and repair processes, where conformational properties may facilitate the proper interactions. We have previously used 31P NMR to investigate the energetics of DNA backbone BI-BII interconversion and the effect of a mismatch or lesion compared to canonical DNA and found stepwise differences in ΔG of 1-2 kcal/mol greater than equivalent steps in unmodified DNA. We have currently compared our results to substrate dependence for TDG, MBD4, M. HhaI, and CEBPβ, testing for correlations to sequence and base-pair dependence. We found strong correlations of our DNA phosphate backbone equilibrium (Keq) to different enzyme kinetics or binding parameters of these varied enzymes, suggesting that the backbone equilibrium may play an important role in mismatch recognition and/or conformational rearrangement and energetics during nucleotide flipping or other aspects of enzyme interrogation of the DNA substrate.

"Flexible hinge" dynamics in mismatched DNA revealed by fluorescence correlation spectroscopy

Altered unwinding/bending fluctuations at DNA lesion sites are implicated as plausible mechanisms for damage sensing by DNA-repair proteins. These dynamics are expected to occur on similar timescales as one-dimensional (1D) diffusion of proteins on DNA if effective in stalling these proteins as they scan DNA. We examined the flexibility and dynamics of DNA oligomers containing 3 base pair (bp) mismatched sites specifically recognized in vitro by nucleotide excision repair protein Rad4 (yeast ortholog of mammalian XPC). A previous Forster resonance energy transfer (FRET) study mapped DNA conformational distributions with cytosine analog FRET pair primarily sensitive to DNA twisting/unwinding deformations (Chakraborty et al. Nucleic Acids Res. 46: 1240-1255 (2018)). These studies revealed B-DNA conformations for nonspecific (matched) constructs but significant unwinding for mismatched constructs specifically recognized by Rad4, even in the absence of Rad4. The timescales of these unwinding fluctuations, however, remained elusive. Here, we labeled DNA with Atto550/Atto647N FRET dyes suitable for fluorescence correlation spectroscopy (FCS). With these probes, we detected higher FRET in specific, mismatched DNA compared with matched DNA, reaffirming unwinding/bending deformations in mismatched DNA. FCS unveiled the dynamics of these spontaneous deformations at ~ 300 µs with no fluctuations detected for matched DNA within the ~ 600 ns-10 ms FCS time window. These studies are the first to visualize anomalous unwinding/bending fluctuations in mismatched DNA on timescales that overlap with the < 500 µs "stepping" times of repair proteins on DNA. Such "flexible hinge" dynamics at lesion sites could arrest a diffusing protein to facilitate damage interrogation and recognition.

Unusual mammalian usage of TGA stop codons reveals that sequence conservation need not imply purifying selection

The assumption that conservation of sequence implies the action of purifying selection is central to diverse methodologies to infer functional importance. GC-biased gene conversion (gBGC), a meiotic mismatch repair bias strongly favouring GC over AT, can in principle mimic the action of selection, this being thought to be especially important in mammals. As mutation is GC→AT biased, to demonstrate that gBGC does indeed cause false signals requires evidence that an AT-rich residue is selectively optimal compared to its more GC-rich allele, while showing also that the GC-rich alternative is conserved. We propose that mammalian stop codon evolution provides a robust test case. Although in most taxa TAA is the optimal stop codon, TGA is both abundant and conserved in mammalian genomes. We show that this mammalian exceptionalism is well explained by gBGC mimicking purifying selection and that TAA is the selectively optimal codon. Supportive of gBGC, we observe (i) TGA usage trends are consistent at the focal stop codon and elsewhere (in UTR sequences); (ii) that higher TGA usage and higher TAA→TGA substitution rates are predicted by a high recombination rate; and (iii) across species the difference in TAA <-> TGA substitution rates between GC-rich and GC-poor genes is largest in genomes that possess higher between-gene GC variation. TAA optimality is supported both by enrichment in highly expressed genes and trends associated with effective population size. High TGA usage and high TAA→TGA rates in mammals are thus consistent with gBGC’s predicted ability to “drive” deleterious mutations and supports the hypothesis that sequence conservation need not be indicative of purifying selection. A general trend for GC-rich trinucleotides to reside at frequencies far above their mutational equilibrium in high recombining domains supports the generality of these results.

Kinetics and thermodynamics of BI-BII interconversion altered by T:G mismatches in DNA

T:G mismatches in DNA result in humans primarily from deamination of methylated CpG sites. They are repaired by redundant systems, such as thymine DNA glycosylase (TDG) and methyl-binding domain enzyme (MBD4), and maintenance of these sites has been implicated in epigenetic processes. The process by which these enzymes identify a canonical DNA base in the incorrect basepairing context remains a mystery. However, the conserved contacts of the repair enzymes with the DNA backbone suggests a role for protein-phosphate interaction in the recognition and repair processes. We have used 31P NMR to investigate the energetics of DNA backbone BI-BII interconversion, and for this work have focused on alterations to the activation barriers to interconversion and the effect of a mismatch compared with canonical DNA. We have found that alterations to the ΔG of interconversion for T:G basepairs are remarkably similar to U:G basepairs in the form of stepwise differences in ΔG of 1-2 kcal/mol greater than equivalent steps in unmodified DNA, suggesting a universality of this result for TDG substrates. Likewise, we see perturbations to the free energy (∼1 kcal/mol) and enthalpy (2-5 kcal/mol) of activation for the BI-BII interconversion localized to the phosphates flanking the mismatch. Overall our results strongly suggest that the perturbed backbone energetics in T:G basepairs play a significant role in the recognition process of DNA repair enzymes.

Single-Base Lesions and Mismatches Alter the Backbone Conformational Dynamics in DNA

DNA damage has been implicated in numerous human diseases, particularly cancer, and the aging process. Single-base lesions and mismatches in DNA can be cytotoxic or mutagenic and are recognized by a DNA glycosylase during the process of base excision repair. Altered local dynamics and conformational properties in damaged DNAs have previously been suggested to assist in recognition and specificity. Herein, we use solution nuclear magnetic resonance to quantify changes in BI-BII backbone conformational dynamics due to the presence of single-base lesions in DNA, including uracil, dihydrouracil, 1,N⁶-ethenoadenine, and T:G mismatches. Stepwise changes to the %BII and ΔG of the BI-BII dynamic equilibrium compared to those of unmodified sequences were observed. Additionally, the equilibrium skews toward endothermicity for the phosphates nearest the lesion/mismatched base pair. Finally, the phosphates with the greatest alterations correlate with those most relevant to the repair of enzyme binding. All of these results suggest local conformational rearrangement of the DNA backbone may play a role in lesion recognition by repair enzymes.

A Highly Fluorescent Nucleobase Molecular Rotor

Fluorescent base analogs (FBAs) are powerful probes of nucleic acids' structures and dynamics. However, previously reported FBAs exhibit relatively low brightness and therefore limited sensitivity of detection. Here we report the hitherto brightest FBA that has ideal molecular rotor properties for detecting local dynamic motions associated with base pair mismatches. The new trans-stilbene annulated uracil derivative "^tsT" exhibits bright fluorescence emissions in various solvents (ε × Φ = 3400-29 700 cm^-1 M^-1) and is highly sensitive to mechanical motions in duplex DNA (ε × Φ = 150-4250 cm^-1 M^-1). ^tsT is thereby a "smart" thymidine analog, exhibiting a 28-fold brighter fluorescence intensity when base paired with A as compared to T or C. Time-correlated single photon counting revealed that the fluorescence lifetime of ^tsT (τ = 4-11 ns) was shorter than its anisotropy decay in well-matched duplex DNA (θ = 20 ns), yet longer than the dynamic motions of base pair mismatches (0.1-10 ns). These properties enable unprecedented sensitivity in detecting local dynamics of nucleic acids.

Melting temperature measurement and mesoscopic evaluation of single, double and triple DNA mismatches

Unlike the canonical base pairs AT and GC, the molecular properties of mismatches such as hydrogen bonding and stacking interactions are strongly dependent on the identity of the neighbouring base pairs. As a result, due to the sheer number of possible combinations of mismatches and flanking base pairs, only a fraction of these have been studied in varying experiments or theoretical models. Here, we report on the melting temperature measurement and mesoscopic analysis of contiguous DNA mismatches in nearest-neighbours and next-nearest neighbour contexts. A total of 4032 different mismatch combinations, including single, double and triple mismatches were covered. These were compared with 64 sequences containing all combinations of canonical base pairs in the same location under the same conditions. For a substantial number of single mismatch configurations, 15%, the measured melting temperatures were higher than the least stable AT base pair. The mesoscopic calculation, using the Peyrard-Bishop model, was performed on the set of 4096 sequences, and resulted in estimates of on-site and nearest-neighbour interactions that can be correlated to hydrogen bonding and base stacking. Our results confirm many of the known properties of mismatches, including the peculiar sheared stacking of tandem GA mismatches. More intriguingly, it also reveals that a number of mismatches present strong hydrogen bonding when flanked on both sites by other mismatches. To highlight the applicability of our results, we discuss a number of practical situations such as enzyme binding affinities, thymine DNA glycosylase repair activity, and trinucleotide repeat expansions.

Mismatch Recognition by Saccharomyces cerevisiae Msh2-Msh6: Role of Structure and Dynamics

The mismatch repair (MMR) pathway maintains genome integrity by correcting errors such as mismatched base pairs formed during DNA replication. In MMR, Msh2–Msh6, a heterodimeric protein, targets single base mismatches and small insertion/deletion loops for repair. By incorporating the fluorescent nucleoside base analog 6-methylisoxanthopterin (6-MI) at or adjacent to a mismatch site to probe the structural and dynamic elements of the mismatch, we address how Msh2–Msh6 recognizes these mismatches for repair within the context of matched DNA. Fluorescence quantum yield and rotational correlation time measurements indicate that local base dynamics linearly correlate with Saccharomyces cerevisiae Msh2–Msh6 binding affinity where the protein exhibits a higher affinity (K_D ≤ 25 nM) for mismatches that have a significant amount of dynamic motion. Energy transfer measurements measuring global DNA bending find that mismatches that are both well and poorly recognized by Msh2–Msh6 experience the same amount of protein-induced bending. Finally, base-specific dynamics coupled with protein-induced blue shifts in peak emission strongly support the crystallographic model of directional binding, in which Phe 432 of Msh6 intercalates 3′ of the mismatch. These results imply an important role for local base dynamics in the initial recognition step of MMR.

Bending of DNA duplexes with mutation motifs

Mutations can be induced by environmental factors but also arise spontaneously during DNA replication or due to deamination of methylated cytosines at CpG dinucleotides. Sites where mutations occur with higher frequency than would be expected by chance are termed hotspots while sites that contain mutations rarely are termed coldspots. Mutations are permanently scanned and repaired by repair systems. Among them, the mismatch repair targets base pair mismatches, which are discriminated from canonical base pairs by probing altered elasticity of DNA. Using biased molecular dynamics simulations, we investigated the elasticity of coldspots and hotspots motifs detected in human genes associated with inherited disorders, and also of motifs with Czech population hotspots and de novo mutations. Main attention was paid to mutations leading to G/T and A+/C pairs. We observed that hotspots without CpG/CpHpG sequences are less flexible than coldspots, which indicates that flexible sequences are more effectively repaired. In contrary, hotspots with CpG/CpHpG sequences exhibited increased flexibility as coldspots. Their mutability is more likely related to spontaneous deamination of methylated cytosines leading to C > T mutations, which are primarily targeted by base excision repair. We corroborated conclusions based on computer simulations by measuring melting curves of hotspots and coldspots containing G/T mismatch.

Enhanced spontaneous DNA twisting/bending fluctuations unveiled by fluorescence lifetime distributions promote mismatch recognition by the Rad4 nucleotide excision repair complex

Rad4/XPC recognizes diverse DNA lesions including ultraviolet-photolesions and carcinogen-DNA adducts, initiating nucleotide excision repair. Studies have suggested that Rad4/XPC senses lesion-induced helix-destabilization to flip out nucleotides from damaged DNA sites. However, characterizing how DNA deformability and/or distortions impact recognition has been challenging. Here, using fluorescence lifetime measurements empowered by a maximum entropy algorithm, we mapped the conformational heterogeneities of artificially destabilized mismatched DNA substrates of varying Rad4-binding specificities. The conformational distributions, as probed by FRET between a cytosine-analog pair exquisitely sensitive to DNA twisting/bending, reveal a direct connection between intrinsic DNA deformability and Rad4 recognition. High-specificity CCC/CCC mismatch, free in solution, sampled a strikingly broad range of conformations from B-DNA-like to highly distorted conformations that resembled those observed with Rad4 bound; the extent of these distortions increased with bound Rad4 and with temperature. Conversely, the non-specific TAT/TAT mismatch had a homogeneous, B-DNA-like conformation. Molecular dynamics simulations also revealed a wide distribution of conformations for CCC/CCC, complementing experimental findings. We propose that intrinsic deformability promotes Rad4 damage recognition, perhaps by stalling a diffusing protein and/or facilitating ‘conformational capture’ of pre-distorted damaged sites. Surprisingly, even mismatched DNA specifically bound to Rad4 remains highly dynamic, a feature that may reflect the versatility of Rad4/XPC to recognize many structurally dissimilar lesions.

MutSγ-Induced DNA Conformational Changes Provide Insights into Its Role in Meiotic Recombination

In many organisms, MutSγ plays a role in meiotic recombination, facilitating crossover formation between homologous chromosomes. Failure to form crossovers leads to improper segregation of chromosomes and aneuploidy, which in humans result in infertility and birth defects. To improve current understanding of MutSγ function, this study investigates the binding affinities and structures of MutSγ in complex with DNA substrates that model homologous recombination intermediates. For these studies, we overexpressed and isolated from Escherichia coli the yeast MutSγ protein Saccharomyces cerevisiae (Sc) Msh4-Msh5. Sc Msh4-Msh5 binds Holliday junction (HJ)-like substrates, 3' overhangs, single-stranded (ss) forks, and the displacement loop with nanomolar affinity. The weakest binding affinities are detected for an intact duplex and open-junction construct. Similar to the human protein, Sc Msh4-Msh5 exhibits the highest affinity for the HJ with a K_d < 0.4 nM in solution. Energy-transfer experiments further demonstrate that DNA structure is modulated by the binding interaction with the largest changes associated with substrates containing an ss end. Upon binding, Sc Msh4-Msh5 displaces the ss away from the duplex in most of the ss-containing intermediates, potentially enabling the binding of RPA and other proteins. In the case of the junction-like intermediates, Msh4-Msh5 binding either stabilizes the existing stacked structure or induces formation of the stacked X conformation. Significantly, we find that upon binding, Msh4-Msh5 stacks an open-junction construct to the same extent as the standard junction. Stabilization of the junction in the stacked conformation is generally refractory to branch migration, which is consistent with a potential role for MutSγ to stabilize HJs and prevent branch migration until resolution by MutLγ. The different binding modalities observed suggest that Msh4-Msh5 not only binds to and stabilizes stacked junctions but also participates in meiotic recombination before junction formation through the stabilization of single-end invasion intermediates.

Active Site Interactions Impact Phosphoryl Transfer during Replication of Damaged and Undamaged DNA by Escherichia coli DNA Polymerase I

Replicative DNA polymerases are able to discriminate between very similar substrates with high accuracy. One mechanism by which E. coli DNA polymerase I checks for Watson-Crick geometry is through a hydrogen bonding fork between Arg668 and the incoming dNTP and the minor groove of the primer terminus. The importance of the Arg-fork was examined by disrupting it with either a guanine to 3-deazaguanine substitution at the primer terminus or the use of a carbocyclic deoxyribose analog of dUTP. Using thio-substituted dNTPs and differential quench techniques, we determined that when the Arg-fork was disrupted, the rate-limiting step changed from a conformational change to phosphodiester bond formation. This result indicates that Arg668 is involved in the phosphoryl transfer step. We examined the role of the Arg-fork in the replication of four DNA damaged templates, O⁶-methylguanine (O⁶-mG), 8-oxo-7,8-dihydroguanine (oxoG), O²-[4-(3-pyridyl)-4-oxobutyl]thymine (O²-POB-T), and N²-[(7S,8R,9S,10R)-7,8,9,10-tetrahydro-8,9,10-trihydroxybenzo[a]pyren-7-yl]-guanine (N²-BP-G). In general, the guanine to 3-deazaguanine substitution caused a decrease in k_pol that was proportional to k_pol over five orders of magnitude. The linear relationship indicates that the Arg668-fork helps catalyze phosphoryl transfer by the same mechanism with all the substrates. Exceptions to the linear relationship were the incorporations of dTTP opposite G, oxoG, and O⁶mG, which showed large decreases in k_pol, similar to that exhibited by the Watson-Crick base pairs. It was proposed that the incorporation of dTTP opposite G, oxoG, and O⁶mG occurred via Watson-Crick-like structures.

The structural impact of DNA mismatches

The structure and dynamics of all the transversion and transition mismatches in three different DNA environments have been characterized by molecular dynamics simulations and NMR spectroscopy. We found that the presence of mismatches produced significant local structural alterations, especially in the case of purine transversions. Mismatched pairs often show promiscuous hydrogen bonding patterns, which interchange among each other in the nanosecond time scale. This therefore defines flexible base pairs, where breathing is frequent, and where distortions in helical parameters are strong, resulting in significant alterations in groove dimension. Even if the DNA structure is plastic enough to absorb the structural impact of the mismatch, local structural changes can be propagated far from the mismatch site, following the expected through-backbone and a previously unknown through-space mechanism. The structural changes related to the presence of mismatches help to understand the different susceptibility of mismatches to the action of repairing proteins.

Dynamic hydrogen bonding and DNA flexibility in minor groove binders: molecular dynamics simulation of the polyamide f-ImPyIm bound to the Mlu1 (MCB) sequence 5'-ACGCGT-3' in 2:1 motif

Molecular dynamics simulations of the DNA 10-mer 5'-CCACGCGTGG-3' alone and complexed with the formamido-imidazole-pyrrole-imidazole (f-ImPyIm) polyamide minor groove binder in a 2:1 fashion were conducted for 50 ns using the pbsc0 parameters within the AMBER 12 software package. The change in DNA structure upon binding of f-ImPyIm was evaluated via minor groove width and depth, base pair parameters of Slide, Twist, Roll, Stretch, Stagger, Opening, Propeller, and x-displacement, dihedral angle distributions of ζ, ε, α, and γ determined using the Curves+ software program, and hydrogen bond formation. The dynamic hydrogen bonding between the f-ImPyIm and its cognate DNA sequence was compared to the static image used to predict sequence recognition by polyamide minor groove binders. Many of the predicted hydrogen bonds were present in less than 50% of the simulation; however, persistent hydrogen bonds between G5/15 and the formamido group of f-ImPyIm were observed. It was determined that the DNA is wider in the Complex than without the polyamide binder; however, there is flexibility in this particular sequence, even in the presence of the f-ImPyIm as evidenced by the range of minor groove widths the DNA exhibits and the dynamics of the hydrogen bonding that binds the two f-ImPyIm ions to the minor groove. The Complex consisting of the DNA and the 2 f-ImPyIm binders shows slight fraying of the 5' end of the 10-mer at the end of the simulation, but the portion of the oligomer responsible for recognition and binding is stable throughout the simulation. Several structural changes in the Complex indicate that minor groove binders may have a more active role in inhibiting transcription than just preventing binding of important transcription factors.

Molecular basis for chiral selection in RNA aminoacylation

The chiral-selective aminoacylation of an RNA minihelix is a potential progenitor to modern tRNA-based protein synthesis using l-amino acids. This article describes the molecular basis for this chiral selection. The extended double helical form of an RNA minihelix with a CCA triplet (acceptor of an amino acid), an aminoacyl phosphate donor nucleotide (mimic of aminoacyl-AMP), and a bridging nucleotide facilitates chiral-selective aminoacylation. Energetically, the reaction is characterized by a downhill reaction wherein an amino acid migrates from a high-energy acyl phosphate linkage to a lower-energy carboxyl ester linkage. The reaction occurs under the restriction that the nucleophilic attack of O, from 3′-OH in the terminal CCA, to C, from C=O in the acyl phosphate linkage, must occur at a Bürgi-Dunitz angle, which is defined as the O–C=O angle of approximately 105°. The extended double helical form results in a steric hindrance at the side chain of the amino acid leading to chiral preference combined with cation coordinations in the amino acid and the phosphate oxygen. Such a system could have developed into the protein biosynthetic system with an exclusively chiral component (l-amino acids) via (proto) ribosomes.

Comparative reactivity of mismatched and unpaired bases in relation to their type and surroundings. Chemical cleavage of DNA mismatches in mutation detection analysis

Systematic study of chemical reactivity of non-Watson-Crick base pairs depending on their type and microenvironment was performed on a model system that represents two sets of synthetic DNA duplexes with all types of mismatched and unmatched bases flanked by T.A or G.C pairs. Using comparative cleavage pattern analysis, we identified the main and additional target bases and performed quantitative study of the time course and efficacy of DNA modification caused by potassium permanganate or hydroxylamine. Potassium permanganate in combination with tetraethylammonium chloride was shown to induce DNA cleavage at all mismatched or bulged T residues, as well as at thymines of neighboring canonical pairs. Other mispaired (bulged) bases and thymine residues located on the second position from the mismatch site were not the targets for KMnO(4) attack. In contrast, hydroxylamine cleaved only heteroduplexes containing mismatched or unmatched C residues, and did not modify adjacent cytosines. However when G.C pairs flank bulged C residue, neighboring cytosines are also attacked by hydroxylamine due to defect migration. Chemical reactivity of target bases was shown to correlate strongly with the local disturbance of DNA double helix at mismatch or bulge site. With our model system, we were able to prove the absence of false-negative and false-positive results. Portion of heteroduplex reliably revealed in a mixture with corresponding homoduplex consists of 5% for bulge bases and "open" non-canonical pairs, and 10% for wobble base pairs giving minimal violations in DNA structure. This study provides a complete understanding of the principles of mutation detection methodology based on chemical cleavage of mismatches and clarifies the advantages and limitations of this approach in various biological and conformational studies of DNA.

Extensive backbone dynamics in the GCAA RNA tetraloop analyzed using 13C NMR spin relaxation and specific isotope labeling

Conformational dynamics play a key role in the properties and functions of proteins and nucleic acids. Heteronuclear NMR spin relaxation is a uniquely powerful site-specific probe of dynamics in proteins and has found increasing applications to nucleotide base side chains and anomeric sites in RNA. Applications to the nucleic acid ribose backbone, however, have been hampered by strong magnetic coupling among ring carbons in uniformly 13C-labeled samples. In this work, we apply a recently developed, metabolically directed isotope labeling scheme that places 13C with high efficiency and specificity at the nucleotide ribose C2' and C4' sites. We take advantage of this scheme to explore backbone dynamics in the well-studied GCAA RNA tetraloop. Using a combination of CPMG (Carr-Purcell-Meiboom-Gill) and R(1rho) relaxation dispersion spectroscopy to explore exchange processes on the microsecond to millisecond time scale, we find an extensive pattern of dynamic transitions connecting a set of relatively well-defined conformations. In many cases, the observed transitions appear to be linked to C3'-endo/C2'-endo sugar pucker transitions of the corresponding nucleotides, and may also be correlated across multiple nucleotides within the tetraloop. These results demonstrate the power of NMR spin relaxation based on alternate-site isotope labeling to open a new window into the dynamic properties of ribose backbone groups in RNA.

Sequence context effect for hMSH2-hMSH6 mismatch-dependent activation

Numerous DNA mismatches and lesions activate MutS homologue (MSH) ATPase activity that is essential for mismatch repair (MMR). We have found that a mismatch embedded in a nearest-neighbor sequence context containing symmetric 3'-purines (2 x 3'-purines) enhanced, whereas symmetric 3'-pyrimidines (2 x 3'-pyrimidines) reduced, hMSH2-hMSH6 ATPase activation. The 3'-purine/pyrimidine effect was most evident for G-containing mispairs. A similar trend pervaded mismatch binding (K(D)) and the melting of unbound oligonucleotides (T(m); DeltaG). However, these latter measures did not accurately predict the hierarchy of MSH ATPase activation. NMR studies of imino proton lifetime, solvent accessibility, and NOE connectivity suggest that sequence contexts that provoke improved MSH-activation displayed enhanced localized DNA flexibility: a dynamic DNA signature that may account for the wide range of lesions that activate MSH functions.

Mechanism of MutS searching for DNA mismatches and signaling repair

DNA mismatch repair is initiated by the recognition of mismatches by MutS proteins. The mechanism by which MutS searches for and recognizes mismatches and subsequently signals repair remains poorly understood. We used single-molecule analyses of atomic force microscopy images of MutS-DNA complexes, coupled with biochemical assays, to determine the distributions of conformational states, the DNA binding affinities, and the ATPase activities of wild type and two mutants of MutS, with alanine substitutions in the conserved Phe-Xaa-Glu mismatch recognition motif. We find that on homoduplex DNA, the conserved Glu, but not the Phe, facilitates MutS-induced DNA bending, whereas at mismatches, both Phe and Glu promote the formation of an unbent conformation. The data reveal an unusual role for the Phe residue in that it promotes the unbending, not bending, of DNA at mismatch sites. In addition, formation of the specific unbent MutS-DNA conformation at mismatches appears to be required for the inhibition of ATP hydrolysis by MutS that signals initiation of repair. These results provide a structural explanation for the mechanism by which MutS searches for and recognizes mismatches and for the observed phenotypes of mutants with substitutions in the Phe-Xaa-Glu motif.

Sequence specific and high affinity recognition of 5'-ACGCGT-3' by rationally designed pyrrole-imidazole H-pin polyamides: thermodynamic and structural studies

Imidazole (Im) and Pyrrole (Py)-containing polyamides that can form stacked dimers can be programmed to target specific sequences in the minor groove of DNA and control gene expression. Even though various designs of polyamides have been thoroughly investigated for DNA sequence recognition, the use of H-pin polyamides (covalently cross-linked polyamides) has not received as much attention. Therefore, experiments were designed to systematically investigate the DNA recognition properties of two symmetrical H-pin polyamides composed of PyImPyIm (5) or f-ImPyIm (3e, f=formamido) tethered with an ethylene glycol linker. These compounds were created to recognize the cognate 5'-ACGCGT-3' through an overlapped and staggered binding motif, respectively. Results from DNaseI footprinting, thermal denaturation, circular dichroism, surface plasmon resonance and isothermal titration microcalorimetry studies demonstrated that both H-pin polyamides bound with higher affinity than their respective monomers. The binding affinity of formamido-containing H-pin 3e was more than a hundred times greater than that for the tetraamide H-pin 5, demonstrating the importance of having a formamido group and the staggered motif in enhancing affinity. However, compared to H-pin 3e, tetraamide H-pin 5 demonstrated superior binding preference for the cognate sequence over its non-cognates, ACCGGT and AAATTT. Data from SPR experiments yielded binding constants of 1.6x10(8)M(-1) and 2.0x10(10)M(-1) for PyImPyIm H-pin 5 and f-ImPyIm H-pin 3e, respectively. Both H-pins bound with significantly higher affinity (ca. 100-fold) than their corresponding unlinked PyImPyIm 4 and f-ImPyIm 2 counterparts. ITC analyses revealed modest enthalpies of reactions at 298 K (DeltaH of -3.3 and -1.0 kcal mol(-1) for 5 and 3e, respectively), indicating these were entropic-driven interactions. The heat capacities (DeltaC(p)) were determined to be -116 and -499 cal mol(-1)K(-1), respectively. These results are in general agreement with DeltaC(p) values determined from changes in the solvent accessible surface areas using complexes of the H-pins bound to (5'-CCACGCGTGG)(2). According to the models, the H-pins fit snugly in the minor groove and the linker comfortably holds both polyamide portions in place, with the oxygen atoms pointing into the solvent. In summary, the H-pin polyamide provides an important molecular design motif for the discovery of future generations of programmable small molecules capable of binding to target DNA sequences with high affinity and selectivity.

Mutations in the signature motif in MutS affect ATP-induced clamp formation and mismatch repair

SUMMARY

MutS protein dimer recognizes and co-ordinates repair of DNA mismatches. Mismatch recognition by the N-terminal mismatch recognition domain and subsequent downstream signalling by MutS appear coupled to the C-terminal ATP catalytic site, Walker box, through nucleotide-mediated conformational transitions. Details of this co-ordination are not understood. The focus of this study is a conserved loop in Escherichia coli MutS that is predicted to mediate cross-talk between the two ATP catalytic sites in MutS homodimer. Mutagenesis was employed to assess the role of this loop in regulating MutS function. All mutants displayed mismatch repair defects in vivo. Biochemical characterization further revealed defects in ATP binding, ATP hydrolysis as well as effective mismatch recognition. The kinetics of initial burst of ATP hydrolysis was similar to wild type but the magnitude of the burst was reduced for the mutants. Given its proximity to the ATP bound in the opposing monomer in the crystal and its potential analogy with signature motif of ABC transporters, the results strongly suggest that the loop co-ordinates ATP binding/hydrolysis in trans by the two catalytic sites. Importantly, our data reveal that the loop plays a direct role in co-ordinating conformational changes involved in long-range communication between Walker box and mismatch recognition domains.

Cytosine ribose flexibility in DNA: a combined NMR 13C spin relaxation and molecular dynamics simulation study

Using (13)C spin relaxation NMR in combination with molecular dynamic (MD) simulations, we characterized internal motions within double-stranded DNA on the pico- to nano-second time scale. We found that the C-H vectors in all cytosine ribose moieties within the Dickerson-Drew dodecamer (5'-CGCGAATTCGCG-3') are subject to high amplitude motions, while the other nucleotides are essentially rigid. MD simulations showed that repuckering is a likely motional model for the cytosine ribose moiety. Repuckering occurs with a time constant of around 100 ps. Knowledge of DNA dynamics will contribute to our understanding of the recognition specificity of DNA-binding proteins such as cytosine methyltransferase.

Modifying the N-terminus of polyamides: PyImPyIm has improved sequence specificity over f-ImPyIm

Seven N-terminus modified derivatives of a previously published minor-groove binding polyamide (f-ImPyIm, 1) were synthesized and the biochemical and biophysical chemistry evaluated. These compounds were synthesized with the aim of attaining a higher level of sequence selectivity over f-ImPyIm (1), a previously published strong minor-groove binder. Two compounds possessing a furan or a benzofuran moiety at the N-terminus showed a footprint of 0.5microM at the cognate ACGCGT site (determined by DNase I footprinting); however, the specificity of these compounds was not improved. In contrast, PyImPyIm (4) produced a footprint of 0.5microM but showed a superior specificity using the same technique. When evaluated by thermal melting experiments and circular dichroism using ACGCGT and the non-cognate AAATTT sequence, all compounds were shown to bind in the minor-groove of DNA and stabilize the cognate sequence much better than the non-cognate (except for the non-amido-compound that did not bind either sequence, as expected). PyImPyIm (4) was interesting as the DeltaT(m) for this compound was only 4 degrees C but the footprint was very selective. No binding was observed for this compound with a third DNA (non-cognate, ACCGGT). ITC studies on compound 4 showed exothermic binding with ACGCGT and no heat change was observed for titrating the compound to the other two DNA sequences. The heat capacity (DeltaC(p)) of the PIPI/ACGCGT complex calculated from the hydrophobic interactions and SASA calculations was comparable to the experimental value obtained from ITC (-146calmol(-1)K(-1)). SPR results provided confirmation of the sequence specificity of PyImPyIm (4), with a K(eq) value determined to be 7.1x10(6) M(-1) for the cognate sequence and no observable binding to AAATTT and ACCGGT. Molecular dynamic simulations affirmed that PyImPyIm (4) binds as a dimer in an overlapped conformation, and it fits snugly in the minor-groove of the ACGCGT oligonucleotide. PyImPyIm (4) is an especially interesting molecule, because although the binding affinity is slightly reduced, the specificity with respect to f-ImPyIm (1) is significantly improved.

Magnitude and direction of DNA bending induced by screw-axis orientation: influence of sequence, mismatches and abasic sites

DNA-bending flexibility is central for its many biological functions. A new bending restraining method for use in molecular mechanics calculations and molecular dynamics simulations was developed. It is based on an average screw rotation axis definition for DNA segments and allows inducing continuous and smooth bending deformations of a DNA oligonucleotide. In addition to controlling the magnitude of induced bending it is also possible to control the bending direction so that the calculation of a complete (2-dimensional) directional DNA-bending map is now possible. The method was applied to several DNA oligonucleotides including A(adenine)-tract containing sequences known to form stable bent structures and to DNA containing mismatches or an abasic site. In case of G:A and C:C mismatches a greater variety of conformations bent in various directions compared to regular B-DNA was found. For comparison, a molecular dynamics implementation of the approach was also applied to calculate the free energy change associated with bending of A-tract containing DNA, including deformations significantly beyond the optimal curvature. Good agreement with available experimental data was obtained offering an atomic level explanation for stable bending of A-tract containing DNA molecules. The DNA-bending persistence length estimated from the explicit solvent simulations is also in good agreement with experiment whereas the adiabatic mapping calculations with a GB solvent model predict a bending rigidity roughly two times larger.

Structure and mechanism for DNA lesion recognition

A fundamental question in DNA repair is how a lesion is detected when embedded in millions to billions of normal base pairs. Extensive structural and functional studies reveal atomic details of DNA repair protein and nucleic acid interactions. This review summarizes seemingly diverse structural motifs used in lesion recognition and suggests a general mechanism to recognize DNA lesion by the poor base stacking. After initial recognition of this shared structural feature of lesions, different DNA repair pathways use unique verification mechanisms to ensure correct lesion identification and removal.

Altered dynamics of DNA bases adjacent to a mismatch: a cue for mismatch recognition by MutS

The structural deviations as well as the alteration in the dynamics of DNA at mismatch sites are considered to have a crucial role in mismatch recognition followed by its repair utilizing mismatch repair family proteins. To compare the dynamics at a mismatch and a non-mismatch site, we incorporated 2-aminopurine, a fluorescent analogue of adenine next to a G.T mismatch, a C.C mismatch, or an unpaired T, and at several other non-mismatch positions. Rotational diffusion of 2-aminopurine at these locations, monitored by time-resolved fluorescence anisotropy, showed distinct differences in the dynamics. This alteration in the motional dynamics is largely confined to the normally matched base-pairs that are immediately adjacent to a mismatch/ unpaired base and could be used by MutS as a cue for mismatch-specific recognition. Interestingly, the enhanced dynamics associated with base-pairs adjacent to a mismatch are significantly restricted upon MutS binding, perhaps "resetting" the cues for downstream events that follow MutS binding. Recognition of such details of motional dynamics of DNA for the first time in the current study enabled us to propose a model that integrates the details of mismatch recognition by MutS as revealed by the high-resolution crystal structure with that of observed base dynamics, and unveils a minimal composite read-out involving the base mismatch and its adjacent normal base-pairs.

Modulation of MutS ATP-dependent functional activities by DNA containing a cisplatin compound lesion (base damage and mismatch)

DNA damage-dependent signaling by the DNA mismatch repair (MMR) system is thought to mediate cytotoxicity of the anti-tumor drug cisplatin through molecular mechanisms that could differ from those required for normal mismatch repair. The present study investigated whether ATP-dependent biochemical properties of Escherichia coli MutS protein differ when the protein interacts with a DNA oligonucleotide containing a GT mismatch versus a unique site specifically placed cisplatin compound lesion, a cisplatin 1,2-d(GpG) intrastrand cross-link with a mispaired thymine opposite the 3' platinated guanine. MutS exhibited substantial affinity for this compound lesion in hydrolytic and in non-hydrolytic conditions of ATP, contrasting with the normal nucleotide inhibition effect of mispair binding. The cisplatin compound lesion was also shown to stimulate poorly MutS ATPase activity to approach the hydrolysis rate induced by nonspecific DNA. Moreover, MutS undergoes distinct conformation changes in the presence of the compound lesion and ATP under hydrolytic conditions as shown by limited proteolysis. In the absence of MutS, the cisplatin compound lesion was shown to induce a 39 degrees rigid bending of the DNA double helix contrasting with an unbent state for DNA containing a GT mispair. Furthermore, an unbent DNA substrate containing a monofunctional adduct mimicking a cisplatin residue failed to form a persistent nucleoprotein complex with MutS in the presence of adenine nucleotide. We propose that DNA bending could play a role in MutS biochemical modulations induced by a compound lesion and that cisplatin DNA damage signaling by the MMR system could be modulated in a direct mode.

NMR assessment of the global shape of a non-labelled DNA dodecamer containing a tandem of G-T mismatches

We have carried out a solution study of two non-labelled self-complementary DNA dodecamers d(GACTGTACAGTC)2 and d(GACTGTGCAGTC)2 by NMR, the second sequence composed of two G-T mismatches. Structures were determined using distances extracted from NOE effects alone or using both NOE and RDC constraints, measured in three different liquid crystalline media. We ensured that our data on the influence of the mesogen on the DNA structures, and the way in which the RDCs were incorporated as constraints in the protocol refinement, were consistent. We also tested the influence of different sets of RDCs and the best means of optimizing the calculation of D(a) and R. Resolution and accuracy of the ten best energy final structures were compared. The addition of a small set of RDC constraints significantly improves the final determined structures. We took advantage of the specificity of the RDC, i.e. it contains orientational information, and explored the global shape of the DNA duplexes; it was found that the duplexes do not have a large curvature. For the G-T base pair, we observed, in this particular sequence (tandem of G-T mismatches), a new pattern of base pairing, which involved the formation of a bifurcated hydrogen bond.

Physical and Structural Basis for the Strong Interactions of the -ImPy- Central Pairing Motif in the Polyamide f-ImPyIm

The polyamide f-ImPyIm has a higher affinity for its cognate DNA than either the parent analogue, distamycin A (10-fold), or the structural isomer, f-PyImIm (250-fold), has for its respective cognate DNA sequence. These findings have led to the formulation of a two-letter polyamide "language" in which the -ImPy- central pairings associate more strongly with Watson-Crick DNA than -PyPy-, -PyIm-, and -ImIm-. Herein, we further characterize f-ImPyIm and f-PyImIm, and we report thermodynamic and structural differences between -ImPy- (f-ImPyIm) and -PyIm- (f-PyImIm) central pairings. DNase I footprinting studies confirmed that f-ImPyIm is a stronger binder than distamycin A and f-PyImIm and that f-ImPyIm preferentially binds CGCG over multiple competing sequences. The difference in the binding of f-ImPyIm and f-PyImIm to their cognate sequences was supported by the Na(+)-dependent nature of DNA melting studies, in which significantly higher Na(+) concentrations were needed to match the ability of f-ImPyIm to stabilize CGCG with that of f-PyImIm stabilizing CCGG. The selectivity of f-ImPyIm beyond the four-base CGCG recognition site was tested by circular dichroism and isothermal titration microcalorimetry, which shows that f-ImPyIm has marginal selectivity for (A.T)CGCG(A.T) over (G.C)CGCG(G.C). In addition, changes adjacent to this 6 bp binding site do not affect f-ImPyIm affinity. Calorimetric studies revealed that binding of f-ImPyIm, f-PyImIm, and distamycin A to their respective hairpin cognate sequences is exothermic; however, changes in enthalpy, entropy, and heat capacity (DeltaC(p)) contribute differently to formation of the 2:1 complexes for each triamide. Experimental and theoretical determinations of DeltaC(p) for binding of f-ImPyIm to CGCG were in good agreement (-142 and -177 cal mol(-)(1) K(-)(1), respectively). (1)H NMR of f-ImPyIm and f-PyImIm complexed with their respective cognate DNAs confirmed positively cooperative formation of distinct 2:1 complexes. The NMR results also showed that these triamides bind in the DNA minor groove and that the oligonucleotide retains the B-form conformation. Using minimal distance restraints from the NMR experiments, molecular modeling and dynamics were used to illustrate the structural complementarity between f-ImPyIm and CGCG. Collectively, the NMR and ITC experiments show that formation of the 2:1 f-ImPyIm-CGCG complex achieves a structure more ordered and more thermodynamically favored than the structure of the 2:1 f-PyImIm-CCGG complex.

Chiral-selective aminoacylation of an RNA minihelix: Mechanistic features and chiral suppression

Aminoacylation of RNA minihelices is speculated to be a key step in the transition from the putative RNA world to the theater of proteins. This reaction affords the opportunity to make chiral selection of an l- or d-amino acid and thus determine the ultimate chirality that is incorporated into proteins. Previous work showed chiral preference of aminoacylation with a nonprotein, nonribozyme, RNA-directed aminoacylation system. This preference was, in turn, determined by the preexisting chirality of the RNA. The alpha-amino group attached to the asymmetric alpha-carbon of the amino acid was an obvious candidate to play a role in chiral selectivity through interactions with the RNA. Also not clear was whether a simple manipulation could change the chiral selectivity, thereby giving insight into the basis of chiral selection in the first place. Here we show, surprisingly, no role for the free alpha-amino group in chiral selection. However, by a sequence manipulation, chiral preference was suppressed and partly reversed. This result and those with further RNA constructs support the idea that the chiral preference for an l-amino acid in these constructs depends on avoiding a sugar-pucker-sensitive steric clash between a pendant group of a base with the amino acid side chain.

Alternate-site isotopic labeling of ribonucleotides for NMR studies of ribose conformational dynamics in RNA

Heteronuclear NMR spin relaxation studies of conformational dynamics are coming into increasing use to help understand the functions of ribozymes and other RNAs. Due to strong 13C-13C magnetic interactions within the ribose ring, however, these studies have thus far largely been limited to (13)C and (15)N resonances on the nucleotide base side chains. We report here the application of the alternate-site (13)C isotopic labeling scheme, pioneered by LeMaster for relaxation studies of amino acid side chains, to nucleic acid systems. We have used different strains of E. coli to prepare mononucleotides containing (13)C label in one of two patterns: Either C1' or C2' in addition to C4', termed (1'/2',4') labeling, or nearly complete labeling at the C2' and C4' sites only, termed (2',4') labeling. These patterns provide isolated 13C-1H spin systems on the labeled carbon atoms and thus allow spin relaxation studies without interference from 13C-13C scalar or dipolar coupling. Using relaxation studies of AMP dissolved in glycerol at varying temperature to produce systems with correlation times characteristic of different size RNAs, we demonstrate the removal of errors due to 13C-13C interaction in T (1) measurements of larger nucleic acids and in T (1rho) measurements in RNA molecules. By extending the applicability of spin relaxation measurements to backbone ribose groups, this technology should greatly improve the flexibility and completeness of NMR analyses of conformational dynamics in RNA.

Poor base stacking at DNA lesions may initiate recognition by many repair proteins

A fundamental question in DNA repair is how a mismatched or modified base is detected when embedded in millions to billions of normal base pairs. A survey of the literature and structural database reveals a common feature in all repair protein-DNA complexes: the DNA double helix is discontinuous at a lesion site due to base unstacking, kinking and/or nucleotide extrusion. Lesions induce destabilization and distortion of short linear DNAs, and underwinding in negatively supercoiled DNA presumably could compound the reduced stability caused by a lesion. A hypothesis is thus put forward that DNA lesion recognition occurs in two steps. Repair proteins initially recognize the weakened base stacking, and thus a flexible hinge at a DNA lesion. Sampling of flexible hinges rather than all DNA base pairs can reduce the task of finding a lesion by two to three orders of magnitude, from searching millions base pairs to thousands. After the initial encounter, a repair protein scrutinizes the shape, hydrogen bonding and electrostatic potentials of bases at the flexible hinge and dissociates if it is not a correct substrate. MutS, which has a broad range of substrates, actively dissociates from non-specific binding via an ATP-dependent proofreading mechanism. A single lesion may thus be sampled by BER, NER and MMR proteins until repaired. This proposition immediately suggests a mechanism for crosstalk between different repair and signaling pathways. It also raises the possibility that sampling of a lesion by one protein could facilitate loading of another by direct protein-protein or DNA mediated interactions.

The structure of a mixed LNA/DNA:RNA duplex is driven by conformational coupling between LNA and deoxyribose residues as determined from 13C relaxation measurements

A study of the internal dynamics of an LNA/DNA:RNA duplex has been performed to further characterize the conformational changes associated with the incorporation of locked nucleic acid (LNA) nucleotides in a DNA:RNA duplex. In general, it was demonstrated that the LNA/DNA:RNA duplex has a very high degree of order compared to dsDNA and dsRNA duplexes. The order parameters of the aromatic carbon atoms in the LNA/DNA strand are uniformly high, whereas a sharp drop in the degree of order was seen in the RNA strand in the beginning of the AUAU stretch in the middle of the strand. This can be related to a return to normal dsRNA dynamics for the central A:U base pair. The high order of the heteroduplex is consistent with preorganization of the chimera strand for an A-form duplex conformation. These results partly explain the dramatic increase in T(m) of the chimeric heteroduplex over dsDNA and DNA:RNA hybrids of the same sequence.

A simple motif for protein recognition in DNA secondary structures

DNA in a single-stranded form (ssDNA) exists transiently within the cell and comprises the telomeres of linear chromosomes and the genomes of some DNA viruses. As with RNA, in the single-stranded state, some DNA sequences are able to fold into complex secondary and tertiary structures that may be recognized by proteins and participate in gene regulation. To better understand how such DNA elements might fold and interact with proteins, and to compare recognition features to those of a structured RNA, we used in vitro selection to identify ssDNAs that bind an RNA-binding peptide from the HIV Rev protein with high affinity and specificity. The large majority of selected binders contain a non-Watson-Crick G.T base-pair and an adjacent C:G base-pair and both are essential for binding. This GT motif can be presented in different DNA contexts, including a nearly perfect duplex and a branched three-helix structure, and appears to be recognized in large part by arginine residues separated by one turn of an alpha-helix. Interestingly, a very similar GT motif is necessary also for protein binding and function of a well-characterized model ssDNA regulatory element from the proenkephalin promoter.

Enzymatic mutation detection technologies

Mutation is as necessary for life as fidelity is in DNA replication. The study of mutations reveals the normal functions of genes, messages, proteins, the causes of many diseases, and the variability of responses among individuals. Indeed, recent mutations that have not yet become polymorphisms are often deleterious and pertinent to the disease history of afflicted individuals. This review discusses the principles behind a variety of methods for the detection of mutations and factors that should be considered in future methods design. One enzymatic approach in particular using orthologs of the CEL I nuclease that show high specificity for all mismatches, appears to be easy and robust. Further developments of this and other methods will allow mutation detection to become an integral component of individualized medicine.

Extending the language of DNA molecular recognition by polyamides: unexpected influence of imidazole and pyrrole arrangement on binding affinity and specificity

Pyrrole (Py) and imidazole (Im) polyamides can be designed to target specific DNA sequences. The effect that the pyrrole and imidazole arrangement, plus DNA sequence, have on sequence specificity and binding affinity has been investigated using DNA melting (DeltaT(M)), circular dichroism (CD), and surface plasmon resonance (SPR) studies. SPR results obtained from a complete set of triheterocyclic polyamides show a dramatic difference in the affinity of f-ImPyIm for its cognate DNA (K(eq) = 1.9 x 10(8) M(-1)) and f-PyPyIm for its cognate DNA (K(eq) = 5.9 x 10(5) M(-1)), which could not have been anticipated prior to characterization of these compounds. Moreover, f-ImPyIm has a 10-fold greater affinity for CGCG than distamycin A has for its cognate, AATT. To understand this difference, the triamide dimers are divided into two structural groupings: central and terminal pairings. The four possible central pairings show decreasing selectivity and affinity for their respective cognate sequences: -ImPy > -PyPy- >> -PyIm- approximately -ImIm-. These results extend the language of current design motifs for polyamide sequence recognition to include the use of "words" for recognizing two adjacent base pairs, rather than "letters" for binding to single base pairs. Thus, polyamides designed to target Watson-Crick base pairs should utilize the strength of -ImPy- and -PyPy- central pairings. The f/Im and f/Py terminal groups yielded no advantage for their respective C/G or T/A base pairs. The exception is with the -ImPy- central pairing, for which f/Im has a 10-fold greater affinity for C/G than f/Py has for T/A.

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.

Molecular recognition of DNA base pairs by the formamido/pyrrole and formamido/imidazole pairings in stacked polyamides

Polyamides containing an N-terminal formamido (f) group bind to the minor groove of DNA as staggered, antiparallel dimers in a sequence-specific manner. The formamido group increases the affinity and binding site size, and it promotes the molecules to stack in a staggered fashion thereby pairing itself with either a pyrrole (Py) or an imidazole (Im). There has not been a systematic study on the DNA recognition properties of the f/Py and f/Im terminal pairings. These pairings were analyzed here in the context of f-ImPyPy, f-ImPyIm, f-PyPyPy and f-PyPyIm, which contain the central pairing modes, -ImPy- and -PyPy-. The specificity of these triamides towards symmetrical recognition sites allowed for the f/Py and f/Im terminal pairings to be directly compared by SPR, CD and DeltaT (M) experiments. The f/Py pairing, when placed next to the -ImPy- or -PyPy- central pairings, prefers A/T and T/A base pairs to G/C base pairs, suggesting that f/Py has similar DNA recognition specificity to Py/Py. With -ImPy- central pairings, f/Im prefers C/G base pairs (>10 times) to the other Watson-Crick base pairs; therefore, f/Im behaves like the Py/Im pair. However, the f/Im pairing is not selective for the C/G base pair when placed next to the -PyPy- central pairings.

Solution structure of a DNA duplex containing an alpha-anomeric adenosine: insights into substrate recognition by endonuclease IV

The cytotoxic alpha anomer of adenosine, generated in situ by radicals, must be recognized and repaired to maintain genomic stability. Endonuclease IV (Endo IV), a member of the base excision repair (BER) enzyme family, in addition to acting on abasic sites, has the auxiliary function of removing this mutagenic nucleotide in Escherichia coli. We have employed enzymatic, thermodynamic, and structural studies on DNA duplexes containing a central alpha-anomeric adenosine residue to characterize the role of DNA structure on recognition and catalysis by Endo IV. The enzyme recognizes and cleaves our alphaA-containing DNA duplexes at the site of the modification. The NMR solution structure of the DNA decamer duplex establishes that the single alpha-anomeric adenosine residue is intrahelical and stacks in a reverse Watson-Crick fashion consistent with the slight decrease in thermostability. However, the presence of this lesion confers significant changes to the global duplex conformation, resulting from a kink of the helical axis into the major groove and an opening of the minor groove emanating from the alpha-anomeric site. Interestingly, the conformation of the flanking base-paired segments is not greatly altered from a B-type conformation. The global structural changes caused by this lesion place the DNA along the conformational path leading to the DNA structure observed in the complex. Thus, it appears that the alpha-anomeric lesion facilitates recognition by Endo IV.

Insight into G[bond]T mismatch recognition using molecular dynamics with time-averaged restraints derived from NMR spectroscopy

Molecular dynamics (MD) simulations were conducted for a G[bond]T mismatch-containing DNA decamer, d(CCATGCGTGG)(2), and its Watson-Crick parent sequence, d(CCACGCGTGG)(2). Dynamics in unrestrained MD trajectories were in poor agreement with prior (13)C NMR studies. However, the accuracy of the trajectories was improved by the use of time-averaged interatomic distance restraints derived from (1)H NMR. Postprocess smoothing of the trajectories further improved accuracy. Comparison of restrained and smoothed trajectories of the two DNA molecules revealed distinct differences in dynamics. The major groove width of the mismatched oligomer was more variable over the course of the simulation compared to its parent sequence. Greater variability in helical parameters stretch and opening for the mismatches indicated less kinetically stable base pairing. Interbase helical parameters rise, roll, and tilt were also more variable in certain base steps involving mismatched bases. These dynamic differences between normal and G[bond]T mismatched DNA reflect differences in local flexibility that may play a role in mismatch recognition by the MutS. A potential alternate G[bond]T mismatch binding mode for MutS is also proposed.

DNA bending and unbending by MutS govern mismatch recognition and specificity

DNA mismatch repair is central to the maintenance of genomic stability. It is initiated by the recognition of base-base mismatches and insertion/deletion loops by the family of MutS proteins. Subsequently, ATP induces a unique conformational change in the MutS-mismatch complex but not in the MutS-homoduplex complex that sets off the cascade of events that leads to repair. To gain insight into the mechanism by which MutS discriminates between mismatch and homoduplex DNA, we have examined the conformations of specific and nonspecific MutS-DNA complexes by using atomic force microscopy. Interestingly, MutS-DNA complexes exhibit a single population of conformations, in which the DNA is bent at homoduplex sites, but two populations of conformations, bent and unbent, at mismatch sites. These results suggest that the specific recognition complex is one in which the DNA is unbent. Combining our results with existing biochemical and crystallographic data leads us to propose that MutS: (i) binds to DNA nonspecifically and bends it in search of a mismatch; (ii) on specific recognition of a mismatch, undergoes a conformational change to an initial recognition complex in which the DNA is kinked, with interactions similar to those in the published crystal structures; and (iii) finally undergoes a further conformational change to the ultimate recognition complex in which the DNA is unbent. Our results provide a structural explanation for the long-standing question of how MutS achieves mismatch repair specificity.