Dynamics of spontaneous flipping of a mismatched base in DNA duplex

Selective tumor cell death induced by irradiated riboflavin through recognizing DNA G-T mismatch

Riboflavin (vitamin B2) has been thought to be a promising antitumoral agent in photodynamic therapy, though the further application of the method was limited by the unclear molecular mechanism. Our work reveals that riboflavin was able to recognize G–T mismatch specifically and induce single-strand breaks in duplex DNA targets efficiently under irradiation. In the presence of riboflavin, the photo-irradiation could induce the death of tumor cells that are defective in mismatch repair system selectively, highlighting the G–T mismatch as potential drug target for tumor cells. Moreover, riboflavin is a promising leading compound for further drug design due to its inherent specific recognition of the G–T mismatch.

Effect of T·T Mismatch on DNA Dynamics Probed by Minor Groove Binders: Comparison of Dynamic Stokes Shifts of Hoechst and DAPI

Recognition of DNA base mismatches and their subsequent repair by enzymes is vital for genomic stability. However, it is difficult to comprehend such a process in which enzymes sense and repair different types of mismatches with different ability. It has been suggested that the differential structural changes of mismatched bases act as cues to the repair enzymes, although the effect of such DNA structural changes on surrounding water and ion dynamics is inevitable due to strong electrostatic coupling among them. Thus, collective dynamics of DNA, water, and ions near the mismatch site is believed to be important for mismatch recognition and repair mechanism. Here we show that introduction of a T·T mismatch in the minor groove of DNA induces dispersed (collective) power-law solvation dynamics (of exponent ∼0.24), measured by monitoring the time-resolved fluorescence Stokes shifts (TRFSS) of two popular minor groove binders (Hoechst 33258 and DAPI) over five decades of time from 100 fs to 10 ns. The same ligands however sense different dynamics (power-law of exponent ∼0.15 or power-law multiplied with biexponential relaxation) in the minor groove of normal-DNA. The similar fluorescence anisotropy decays of ligands measured in normal- and T·T-DNA suggest that Stokes shift dynamics and their changes in T·T-DNA purely originate from the solvation process, and not from any internal rotational motion of probe-ligands. The dispersed power-law solvation dynamics seen in T·T-DNA indicate that the ligands do not sense any particular (exponential) relaxation specific to T·T wobbling and/or other conformational changes. This could be the reason why T·T mismatch is recognized by enzymes with lower efficiency compared to purine-pyrimidine and purine-purine mismatches.

Free Energy Coupling between DNA Bending and Base Flipping

Free energy simulations are presented to probe the energetic coupling between DNA bending and the flipping of a central thymine in double stranded DNA 13mers. The energetics are shown to depend on the neighboring base pairs, and upstream C or T or downstream C tended to make flipping more costly. Flipping to the major groove side was generally preferred. Bending aids flipping, by pushing the system up in free energy, but for small and intermediate bending angles the two were uncorrelated. At higher bending angles, bending and flipping became correlated, and bending primed the system for base flipping toward the major groove. Flipping of the 6-4 pyrimidine-pyrimidone and pyrimidine dimer photoproducts is shown to be more facile than for undamaged DNA. For the damages, major groove flipping was preferred, and DNA bending was much facilitated in the 6-4 pyrimidine-pyrimidone damaged system. Aspects of the calculations were verified by structural analyses of protein-DNA complexes with flipped bases.

Direct observation of breathing dynamics at the mismatch induced DNA bubble with nanometre accuracy: a smFRET study

The detailed conformational dynamics of the melted region in double-stranded DNA has been studied using a combination of ensemble and single-molecule FRET techniques. We monitored the millisecond time scale fluctuation kinetics of the two strands at the bubble region that varies with the size of the bubble. As the individual strands at the melting bubble behave as single-stranded DNA, and hence fluctuate dynamically to attain energetically favored configurations, the rates of these fluctuations increase with increase in the bubble size. In different short DNAs under investigation, the two strands never cross each other to form a knot, irrespective of the number of base pair mismatches present. Rather, they prefer to stay apart from each other, as the size of the bubble increases and follow exactly an opposite trend for bubbles of smaller size. The range within which the bubble strands fluctuate are monitored with great accuracy in the nanometre resolution from the single-molecule FRET measurements. The shape of the bubble that plays a crucial role in determining the activity of the DNA was speculated. These results shall be useful in quantifying the chemical processes within DNA as well as to develop a deeper understanding of the activity of the DNA due to induced mismatches.

Sequence dependency of canonical base pair opening in the DNA double helix

The flipping-out of a DNA base from the double helical structure is a key step of many cellular processes, such as DNA replication, modification and repair. Base pair opening is the first step of base flipping and the exact mechanism is still not well understood. We investigate sequence effects on base pair opening using extensive classical molecular dynamics simulations targeting the opening of 11 different canonical base pairs in two DNA sequences. Two popular biomolecular force fields are applied. To enhance sampling and calculate free energies, we bias the simulation along a simple distance coordinate using a newly developed adaptive sampling algorithm. The simulation is guided back and forth along the coordinate, allowing for multiple opening pathways. We compare the calculated free energies with those from an NMR study and check assumptions of the model used for interpreting the NMR data. Our results further show that the neighboring sequence is an important factor for the opening free energy, but also indicates that other sequence effects may play a role. All base pairs are observed to have a propensity for opening toward the major groove. The preferred opening base is cytosine for GC base pairs, while for AT there is sequence dependent competition between the two bases. For AT opening, we identify two non-canonical base pair interactions contributing to a local minimum in the free energy profile. For both AT and CG we observe long-lived interactions with water and with sodium ions at specific sites on the open base pair.

The DNA double helix, a molecule that stores biological information, has become an iconic image of biomedical research. In order to use or repair the information it carries, the bases that are stacked in the helix need to be chemically exposed. This can happen either by separating the two strands in the helix or by flipping out individual bases. Here, we focus on the latter process. Usually proteins are involved in interactions with bases, but it is still unclear if bases are pulled out actively by proteins or if they act on spontaneously flipped bases. Although experiments can detect base pair opening, it is difficult to detect which base moves in which direction. Here, we present results from molecular dynamics simulations using a recently developed sampling method which improves the statistics in the simulations by enhancing the probability of the base pair opening event. We observe differences in probability, modes and mechanism of opening that depend not only on the types of the bases in the pair, but also strongly on their neighbors. This provides essential information for understanding how DNA functions.

Single-molecule dissection of stacking forces in DNA

We directly measured at the single-molecule level the forces and lifetimes of DNA base-pair stacking interactions for all stack sequence combinations. Our experimental approach combined dual-beam optical tweezers with DNA origami components to allow positioning of blunt-end DNA helices so that the weak stacking force could be isolated. Base-pair stack arrays that lacked a covalent backbone connection spontaneously dissociated at average rates ranging from 0.02 to 500 per second, depending on the sequence combination and stack array size. Forces in the range from 2 to 8 piconewtons that act along the helical direction only mildly accelerated the stochastic unstacking process. The free-energy increments per stack that we estimate from the measured forward and backward kinetic rates ranged from -0.8 to -3.4 kilocalories per mole, depending on the sequence combination. Our data contributes to understanding the mechanics of DNA processing in biology, and it is helpful for designing the kinetics of DNA-based nanoscale devices according to user specifications.

Efficient sampling over rough energy landscapes with high barriers: A combination of metadynamics with integrated tempering sampling

In order to efficiently overcome high free energy barriers embedded in a complex energy landscape and calculate overall thermodynamics properties using molecular dynamics simulations, we developed and implemented a sampling strategy by combining the metadynamics with (selective) integrated tempering sampling (ITS/SITS) method. The dominant local minima on the potential energy surface (PES) are partially exalted by accumulating history-dependent potentials as in metadynamics, and the sampling over the entire PES is further enhanced by ITS/SITS. With this hybrid method, the simulated system can be rapidly driven across the dominant barrier along selected collective coordinates. Then, ITS/SITS ensures a fast convergence of the sampling over the entire PES and an efficient calculation of the overall thermodynamic properties of the simulation system. To test the accuracy and efficiency of this method, we first benchmarked this method in the calculation of ϕ - ψ distribution of alanine dipeptide in explicit solvent. We further applied it to examine the design of template molecules for aromatic meta-C-H activation in solutions and investigate solution conformations of the nonapeptide Bradykinin involving slow cis-trans isomerizations of three proline residues.

Structure of the EndoMS-DNA Complex as Mismatch Restriction Endonuclease

Archaeal NucS nuclease was thought to degrade the single-stranded region of branched DNA, which contains flapped and splayed DNA. However, recent findings indicated that EndoMS, the orthologous enzyme of NucS, specifically cleaves double-stranded DNA (dsDNA) containing mismatched bases. In this study, we determined the structure of the EndoMS-DNA complex. The complex structure of the EndoMS dimer with dsDNA unexpectedly revealed that the mismatched bases were flipped out into binding sites, and the overall architecture most resembled that of restriction enzymes. The structure of the apo form was similar to the reported structure of Pyrococcus abyssi NucS, indicating that movement of the C-terminal domain from the resting state was required for activity. In addition, a model of the EndoMS-PCNA-DNA complex was preliminarily verified with electron microscopy. The structures strongly support the idea that EndoMS acts in a mismatch repair pathway.

Scanning Single-Molecule Fluorescence Correlation Spectroscopy Enables Kinetics Study of DNA Hairpin Folding with a Time Window from Microseconds to Seconds

Single-molecule fluorescence measurements have been widely used to explore kinetics and dynamics of biological systems. Among them, single-molecule imaging (SMI) is good at tracking processes slower than tens of milliseconds, whereas fluorescence correlation spectroscopy (FCS) is good at probing processes faster than submilliseconds. However, there is still shortage of simple yet effective single-molecule fluorescence method to cover the time-scale between submilliseconds and tens of milliseconds. To effectively bridge this millisecond gap, we developed a single-molecule fluorescence correlation spectroscopy (smFCS) method that works on surface-immobilized single molecules through surface scanning. We validated it by monitoring the classical DNA hairpin folding process. With a wide time window from microseconds to seconds, the experimental data are well fitted to the two-state folding model. All relevant molecular parameters, including the relative fluorescence brightness, equilibrium constant, and reaction rate constants, were uniquely determined.

Base Flipping within the α-Hemolysin Latch Allows Single-Molecule Identification of Mismatches in DNA

A method for identifying and differentiating DNA duplexes containing the mismatched base pairs CC and CA at single molecule resolution with the protein pore α-hemolysin (αHL) is presented. Unique modulating current signatures are observed for duplexes containing the CC and CA mismatches when the mismatch site in the duplex is situated in proximity to the latch constriction of αHL during DNA residence inside the pore. The frequency and current amplitude of the modulation states are dependent on the mismatch type (CC or CA) permitting easy discrimination of these mismatches from one another, and from a fully complementary duplex that exhibits no modulation. We attribute the modulating current signatures to base flipping and subsequent interaction with positively charged lysine residues at the latch constriction of αHL. Our hypothesis is supported by the extended residence times of DNA duplexes within the pore when a mismatch is in proximity to the latch constriction, and by the loss of the two-state current signature in low pH buffers (<6.3), where the protonation of one of the cytosine bases increases the stability of the intrahelical state.

Conformational stabilities of iminoallantoin and its base pairs in DNA: implications for mutagenicity

Under acidic conditions, insertion of G opposite Ia may lead to G to C mutations in DNA.

Formaldehyde crosslinking: a tool for the study of chromatin complexes

Formaldehyde has been used for decades to probe macromolecular structure and function and to trap complexes, cells, and tissues for further analysis. Formaldehyde crosslinking is routinely employed for detection and quantification of protein-DNA interactions, interactions between chromatin proteins, and interactions between distal segments of the chromatin fiber. Despite widespread use and a rich biochemical literature, important aspects of formaldehyde behavior in cells have not been well described. Here, we highlight features of formaldehyde chemistry relevant to its use in analyses of chromatin complexes, focusing on how its properties may influence studies of chromatin structure and function.

Effects of non-CpG site methylation on DNA thermal stability: a fluorescence study

Cytosine methylation is a widespread epigenetic regulation mechanism. In healthy mature cells, methylation occurs at CpG dinucleotides within promoters, where it primarily silences gene expression by modifying the binding affinity of transcription factors to the promoters. Conversely, a recent study showed that in stem cells and cancer cell precursors, methylation also occurs at non-CpG pairs and involves introns and even gene bodies. The epigenetic role of such methylations and the molecular mechanisms by which they induce gene regulation remain elusive. The topology of both physiological and aberrant non-CpG methylation patterns still has to be detailed and could be revealed by using the differential stability of the duplexes formed between site-specific oligonucleotide probes and the corresponding methylated regions of genomic DNA. Here, we present a systematic study of the thermal stability of a DNA oligonucleotide sequence as a function of the number and position of non-CpG methylation sites. The melting temperatures were determined by monitoring the fluorescence of donor-acceptor dual-labelled oligonucleotides at various temperatures. An empirical model that estimates the methylation-induced variations in the standard values of hybridization entropy and enthalpy was developed.

Visualizing transient Watson-Crick-like mispairs in DNA and RNA duplexes

Rare tautomeric and anionic nucleobases are believed to have fundamental biological roles, but their prevalence and functional importance has remained elusive because they exist transiently, in low abundance, and involve subtle movements of protons that are difficult to visualize. Using NMR relaxation dispersion, we show here that wobble dG•dT and rG•rU mispairs in DNA and RNA duplexes exist in dynamic equilibrium with short-lived, low-populated Watson-Crick-like mispairs that are stabilized by rare enolic or anionic bases. These mispairs can evade Watson-Crick fidelity checkpoints and form with probabilities (10(-3) to 10(-5)) that strongly imply a universal role in replication and translation errors. Our results indicate that rare tautomeric and anionic bases are widespread in nucleic acids, expanding their structural and functional complexity beyond that attainable with canonical bases.

Lesion search and recognition by thymine DNA glycosylase revealed by single molecule imaging

The ability of DNA glycosylases to rapidly and efficiently detect lesions among a vast excess of nondamaged DNA bases is vitally important in base excision repair (BER). Here, we use single molecule imaging by atomic force microscopy (AFM) supported by a 2-aminopurine fluorescence base flipping assay to study damage search by human thymine DNA glycosylase (hTDG), which initiates BER of mutagenic and cytotoxic G:T and G:U mispairs in DNA. Our data reveal an equilibrium between two conformational states of hTDG–DNA complexes, assigned as search complex (SC) and interrogation complex (IC), both at target lesions and undamaged DNA sites. Notably, for both hTDG and a second glycosylase, hOGG1, which recognizes structurally different 8-oxoguanine lesions, the conformation of the DNA in the SC mirrors innate structural properties of their respective target sites. In the IC, the DNA is sharply bent, as seen in crystal structures of hTDG lesion recognition complexes, which likely supports the base flipping required for lesion identification. Our results support a potentially general concept of sculpting of glycosylases to their targets, allowing them to exploit the energetic cost of DNA bending for initial lesion sensing, coupled with continuous (extrahelical) base interrogation during lesion search by DNA glycosylases.