DNA bending propensity in the presence of base mismatches: implications for DNA repair

A high throughput single molecule platform to study DNA supercoiling effect on protein-DNA interactions

DNA supercoiling significantly influences DNA metabolic pathways. To examine its impact on DNA-protein interactions at the single-molecule level, we developed a highly efficient and reliable protocol to modify plasmid DNA at specific sites, allowing us to label plasmids with fluorophores and biotin. We then induced negative and positive supercoiling in these plasmids using gyrase and reverse gyrase, respectively. Comparing supercoiled DNA with relaxed circular DNA, we assessed the effects of supercoiling on CRISPR-Cas9 and mismatch repair protein MutS. We found that negative DNA supercoiling exacerbates off-target effects in DNA unwinding by Cas9. For MutS, we observed both negative and positive DNA supercoiling enhances the binding interaction between MutS and a mismatched base pair but does not affect the rate of ATP-induced sliding clamp formation. These findings not only underscore the versatility of our protocol but also opens new avenues for exploring the intricate dynamics of protein-DNA interactions under the influences of supercoiling.

Re-pairing DNA: binding of a ruthenium phi complex to a double mismatch

We report a crystal structure at atomic resolution (0.9 Å) of a ruthenium complex bound to a consecutive DNA double mismatch, which results in a TA basepair with flipped out thymine, together with the formation of an adenine bulge. The structure shows a form of metalloinsertion interaction of the Λ-[Ru(phen)2phi]2+ (phi = 9,10-phenanthrenediimine) complex at the bulge site. The metal complex interacts with the DNA via the major groove, where specific interactions between the adenines of the DNA and the phen ligands of the complex are formed. One Δ-[Ru(phen)2phi]2+ complex interacts via the minor groove, which shows sandwiching of its phi ligand between the phi ligands of the other two ruthenium complexes, and no interaction of its phen ligands with DNA. To our knowledge, this binding model represents a new form of metalloinsertion in showing major rather than minor groove insertion.

Molecular dynamics of mismatch detection-How MutS uses indirect readout to find errors in DNA

The mismatch repair protein MutS safeguards genomic integrity by finding and initiating repair of basepairing errors in DNA. Single-molecule studies show MutS diffusing on DNA, presumably scanning for mispaired/unpaired bases, and crystal structures show a characteristic "mismatch-recognition" complex with DNA enclosed within MutS and kinked at the site of error. But how MutS goes from scanning thousands of Watson-Crick basepairs to recognizing rare mismatches remains unanswered, largely because atomic-resolution data on the search process are lacking. Here, 10 μs all-atom molecular dynamics simulations of Thermus aquaticus MutS bound to homoduplex DNA and T-bulge DNA illuminate the structural dynamics underlying the search mechanism. MutS-DNA interactions constitute a multistep mechanism to check DNA over two helical turns for its 1) shape, through contacts with the sugar-phosphate backbone, 2) conformational flexibility, through bending/unbending engineered by large-scale motions of the clamp domain, and 3) local deformability, through basepair destabilizing contacts. Thus, MutS can localize a potential target by indirect readout due to lower energetic costs of bending mismatched DNA and identify a site that distorts easily due to weaker base stacking and pairing as a mismatch. The MutS signature Phe-X-Glu motif can then lock in the mismatch-recognition complex to initiate repair.

Computational Modeling of DNA 3D Structures: From Dynamics and Mechanics to Folding

DNA carries the genetic information required for the synthesis of RNA and proteins and plays an important role in many processes of biological development. Understanding the three-dimensional (3D) structures and dynamics of DNA is crucial for understanding their biological functions and guiding the development of novel materials. In this review, we discuss the recent advancements in computer methods for studying DNA 3D structures. This includes molecular dynamics simulations to analyze DNA dynamics, flexibility, and ion binding. We also explore various coarse-grained models used for DNA structure prediction or folding, along with fragment assembly methods for constructing DNA 3D structures. Furthermore, we also discuss the advantages and disadvantages of these methods and highlight their differences.

Noncanonical DNA structures are drivers of genome evolution

In addition to the canonical right-handed double helix, other DNA structures, termed 'non-B DNA', can form in the genomes across the tree of life. Non-B DNA regulates multiple cellular processes, including replication and transcription, yet its presence is associated with elevated mutagenicity and genome instability. These discordant cellular roles fuel the enormous potential of non-B DNA to drive genomic and phenotypic evolution. Here we discuss recent studies establishing non-B DNA structures as novel functional elements subject to natural selection, affecting evolution of transposable elements (TEs), and specifying centromeres. By highlighting the contributions of non-B DNA to repeated evolution and adaptation to changing environments, we conclude that evolutionary analyses should include a perspective of not only DNA sequence, but also its structure.

MutS recognition of mismatches within primed DNA replication intermediates

MutS initiates mismatch repair by recognizing mismatches in newly replicated DNA. Specific interactions between MutS and mismatches within double-stranded DNA promote ADP-ATP exchange and a conformational change into a sliding clamp. Here, we demonstrated that MutS from Pseudomonas aeruginosa associates with primed DNA replication intermediates. The predicted structure of this MutS-DNA complex revealed a new DNA binding site, in which Asn 279 and Arg 272 appeared to directly interact with the 3'-OH terminus of primed DNA. Mutation of these residues resulted in a noticeable defect in the interaction of MutS with primed DNA substrates. Remarkably, MutS interaction with a mismatch within primed DNA induced a compaction of the protein structure and impaired the formation of an ATP-bound sliding clamp. Our findings reveal a novel DNA binding mode, conformational change and intramolecular signaling for MutS recognition of mismatches within primed DNA structures.

"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.

Stretching of Long Double-Stranded DNA and RNA Described by the Same Approach

We propose an approach to help interpret polymer force-extension curves that exhibit plateau regimes. When coupled to a bead-spring dynamic model, the approach accurately reproduces a variety of experimental force-extension curves of long double-stranded DNA and RNA, including torsionally constrained and unconstrained DNA and negatively supercoiled DNA. A key feature of the model is a specific nonconvex energy function of the spring. We provide an algorithm to obtain the five required parameters of the model from experimental force-extension curves. The applicability of the approach to the force-extension curves of double-stranded (ds) DNA of variable GC content as well as to a DNA/RNA hybrid structure is explored and confirmed. We use the approach to explain counterintuitive sequence-dependent trends and make predictions. In the plateau region of the force-extension curves, our molecular dynamics simulations show that the polymer separates into a mix of weakly and strongly stretched states without forming macroscopically distinct phases. The distribution of these states is predicted to depend on the sequence.

Bending of Canonical and G/T Mismatched DNAs

Mismatched base pairs alter the flexibility and intrinsic curvature of DNA. The role of such DNA features is not fully understood in the mismatch repair pathway. MutS/DNA complexes exhibit DNA bending, PHE intercalation, and changes of base-pair parameters near the mismatch. Recently, we have shown that base-pair opening in the absence of MutS can discriminate mismatches from canonical base pairs better than DNA bending. However, DNA bending in the absence of MutS was found to be rather challenging to describe correctly. Here, we present a computational study on the DNA bending of canonical and G/T mismatched DNAs. Five types of geometric parameters covering template-based bending toward the experimental DNA structure, global, and local geometry parameters were employed in biased molecular dynamics in the absence of MutS. None of these parameters showed higher discrimination than the base-pair opening. Only roll could induce a sharply localized bending of DNA as observed in the experimental MutS/DNA structure. Further, we demonstrated that the intercalation of benzene mimicking PHE decreases the energetic cost of DNA bending without any effect on mismatch discrimination.

Molecular Mechanisms Associated with Clustered Lesion-Induced Impairment of 8-oxoG Recognition by the Human Glycosylase OGG1

The 8-oxo-7,8-dihydroguanine, referred to as 8-oxoG, is a highly mutagenic DNA lesion that can provoke the appearance of mismatches if it escapes the DNA Damage Response. The specific recognition of its structural signature by the hOGG1 glycosylase is the first step along the Base Excision Repair pathway, which ensures the integrity of the genome by preventing the emergence of mutations. 8-oxoG formation, structural features, and repair have been matters of extensive research; more recently, this active field of research expended to the more complicated case of 8-oxoG within clustered lesions. Indeed, the presence of a second lesion within 1 or 2 helix turns can dramatically impact the repair yields of 8-oxoG by glycosylases. In this work, we use μs-range molecular dynamics simulations and machine-learning-based postanalysis to explore the molecular mechanisms associated with the recognition of 8-oxoG by hOGG1 when embedded in a multiple-lesion site with a mismatch in 5′ or 3′. We delineate the stiffening of the DNA–protein interactions upon the presence of the mismatches, and rationalize the much lower repair yields reported with a 5′ mismatch by describing the perturbation of 8-oxoG structural features upon addition of an adjacent lesion.

DNA mechanics and its biological impact

Almost all nucleoprotein interactions and DNA manipulation events involve mechanical deformations of DNA. Extraordinary progresses in single-molecule, structural, and computational methods have characterized the average mechanical properties of DNA, such as bendability and torsional rigidity, in high resolution. Further, the advent of sequencing technology has permitted measuring, in high-throughput, how such mechanical properties vary with sequence and epigenetic modifications along genomes. We review these recent technological advancements, and discuss how they have contributed to the emerging idea that variations in the mechanical properties of DNA play a fundamental role in regulating, genome-wide, diverse processes involved in chromatin organization.

Coarse-grained modelling of DNA plectoneme pinning in the presence of base-pair mismatches

Damaged or mismatched DNA bases result in the formation of physical defects in double-stranded DNA. In vivo, defects in DNA must be rapidly and efficiently repaired to maintain cellular function and integrity. Defects can also alter the mechanical response of DNA to bending and twisting constraints, both of which are important in defining the mechanics of DNA supercoiling. Here, we use coarse-grained molecular dynamics (MD) simulation and supporting statistical-mechanical theory to study the effect of mismatched base pairs on DNA supercoiling. Our simulations show that plectoneme pinning at the mismatch site is deterministic under conditions of relatively high force (>2 pN) and high salt concentration (>0.5 M NaCl). Under physiologically relevant conditions of lower force (0.3 pN) and lower salt concentration (0.2 M NaCl), we find that plectoneme pinning becomes probabilistic and the pinning probability increases with the mismatch size. These findings are in line with experimental observations. The simulation framework, validated with experimental results and supported by the theoretical predictions, provides a way to study the effect of defects on DNA supercoiling and the dynamics of supercoiling in molecular detail.

Importance of base-pair opening for mismatch recognition

Mismatch repair is a highly conserved cellular pathway responsible for repairing mismatched dsDNA. Errors are detected by the MutS enzyme, which most likely senses altered mechanical property of damaged dsDNA rather than a specific molecular pattern. While the curved shape of dsDNA in crystallographic MutS/DNA structures suggests the role of DNA bending, the theoretical support is not fully convincing. Here, we present a computational study focused on a base-pair opening into the minor groove, a specific base-pair motion observed upon interaction with MutS. Propensities for the opening were evaluated in terms of two base-pair parameters: Opening and Shear. We tested all possible base pairs in anti/anti, anti/syn and syn/anti orientations and found clear discrimination between mismatches and canonical base-pairs only for the opening into the minor groove. Besides, the discrimination gap was also confirmed in hotspot and coldspot sequences, indicating that the opening could play a more significant role in the mismatch recognition than previously recognized. Our findings can be helpful for a better understanding of sequence-dependent mutability. Further, detailed structural characterization of mismatches can serve for designing anti-cancer drugs targeting mismatched base pairs.

Automated AFM analysis of DNA bending reveals initial lesion sensing strategies of DNA glycosylases

Base excision repair is the dominant DNA repair pathway of chemical modifications such as deamination, oxidation, or alkylation of DNA bases, which endanger genome integrity due to their high mutagenic potential. Detection and excision of these base lesions is achieved by DNA glycosylases. To investigate the remarkably high efficiency in target site search and recognition by these enzymes, we applied single molecule atomic force microscopy (AFM) imaging to a range of glycosylases with structurally different target lesions. Using a novel, automated, unbiased, high-throughput analysis approach, we were able to resolve subtly different conformational states of these glycosylases during DNA lesion search. Our results lend support to a model of enhanced lesion search efficiency through initial lesion detection based on altered mechanical properties at lesions. Furthermore, its enhanced sensitivity and easy applicability also to other systems recommend our novel analysis tool for investigations of diverse, fundamental biological interactions.

Investigating the pathogenic SNPs in BLM helicase and their biological consequences by computational approach

The BLM helicase protein plays a vital role in DNA replication and the maintenance of genomic integrity. Variation in the BLM helicase gene resulted in defects in the DNA repair mechanism and was reported to be associated with Bloom syndrome (BS) and cancer. Despite extensive investigation of helicase proteins in humans, no attempt has previously been made to comprehensively analyse the single nucleotide polymorphism (SNPs) of the BLM gene. In this study, a comprehensive analysis of SNPs on the BLM gene was performed to identify, characterize and validate the pathogenic SNPs using computational approaches. We obtained SNP data from the dbSNP database version 150 and mapped these data to the genomic coordinates of the "NM_000057.3" transcript expressing BLM helicase (P54132). There were 607 SNPs mapped to missense, 29 SNPs mapped to nonsense, and 19 SNPs mapped to 3'-UTR regions. Initially, we used many consensus tools of SIFT, PROVEAN, Condel, and PolyPhen-2, which together increased the accuracy of prediction and identified 18 highly pathogenic non-synonymous SNPs (nsSNPs) out of 607 SNPs. Subsequently, these 18 high-confidence pathogenic nsSNPs were analysed for BLM protein stability, structure-function relationships and disease associations using various bioinformatics tools. These 18 mutants of the BLM protein along with the native protein were further investigated using molecular dynamics simulations to examine the structural consequences of the mutations, which might reveal their malfunction and contribution to disease. In addition, 28 SNPs were predicted as "stop gained" nonsense SNPs and one SNP was predicted as "start lost". Two SNPs in the 3'UTR were found to abolish miRNA binding and thus may enhance the expression of BLM. Interestingly, we found that BLM mRNA overexpression is associated with different types of cancers. Further investigation showed that the dysregulation of BLM is associated with poor overall survival (OS) for lung and gastric cancer patients and hence led to the conclusion that BLM has the potential to be used as an important prognostic marker for the detection of lung and gastric cancer.

Key structural motifs in Thymine DNA glycosylase responsible for recognizing certain DNA bent conformation revealed by atomic simulations

Knowledge of how DNA bending facilitates the target-base searching by Thymine DNA glycosylase (TDG) is of major importance for unraveling the recognition mechanism between DNA and TDG in DNA repair process. An atomic-level understanding of the initial encounter between TDG and DNA before base-flipping, however, is still elusive. Here, we employ all-atom molecular dynamics (MD) simulations with an integrated simulation time of ∼3 μs to investigate how TDG responses to different DNA bending conformations. By constructing several TDG-DNA complexes with varied DNA bend angles (ranging from ∼0° to 60°), we pinpoint the key TDG motifs responsible for recognizing certain DNA bending conformations. Particularly, several positively charged residues, i.e., Lys232, Lys240, and Lys246, are critical for the tight binding with DNA backbones. Importantly, the roll-angle patterns, rather than the tilt and twist angles, are found to be strongly correlated with the extent of DNA bending, which in turn, governs the TDG recognition. Further comparisons between the naked and TDG-bound DNA conformations reveal that the TDG binding can impose a substantial DNA deformation, resulting in profound roll-angle alterations. Our studies warrant further experimental validations and provide deep structural insights into the recognition mechanism between TDG and DNA during their initial encounter.

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.

Cooperative recognition of T:T mismatch by echinomycin causes structural distortions in DNA duplex

Small-molecule compounds that target mismatched base pairs in DNA offer a novel prospective for cancer diagnosis and therapy. The potent anticancer antibiotic echinomycin functions by intercalating into DNA at CpG sites. Surprisingly, we found that the drug strongly prefers to bind to consecutive CpG steps separated by a single T:T mismatch. The preference appears to result from enhanced cooperativity associated with the binding of the second echinomycin molecule. Crystallographic studies reveal that this preference originates from the staggered quinoxaline rings of the two neighboring antibiotic molecules that surround the T:T mismatch forming continuous stacking interactions within the duplex. These and other associated changes in DNA conformation allow the formation of a minor groove pocket for tight binding of the second echinomycin molecule. We also show that echinomycin displays enhanced cytotoxicity against mismatch repair-deficient cell lines, raising the possibility of repurposing the drug for detection and treatment of mismatch repair-deficient cancers.

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.

Base-Pair Mismatch Can Destabilize Small DNA Loops through Cooperative Kinking

Base-pair mismatch can relieve mechanical stress in highly strained DNA molecules, but how it affects their kinetic stability is not known. Using single-molecule fluorescence resonance energy transfer, we measured the lifetimes of tightly bent DNA loops with and without base-pair mismatch. Surprisingly, for loops captured by stackable sticky ends which leave single-stranded DNA breaks (or nicks) upon annealing, the mismatch decreased the loop lifetime despite reducing the overall bending stress, and the decrease was largest when the mismatch was placed at the DNA midpoint. These findings suggest that base-pair mismatch increases bending stress at the opposite side of the loop through an allosteric mechanism known as cooperative kinking. Based on this mechanism, we present a three-state model that explains the apparent dichotomy between thermodynamic and kinetic stability.

Molecular Mechanisms of DNA Replication and Repair Machinery: Insights from Microscopic Simulations

Reproduction, the hallmark of biological activity, requires making an accurate copy of the genetic material to allow the progeny to inherit parental traits. In all living cells, the process of DNA replication is carried out by a concerted action of multiple protein species forming a loose protein-nucleic acid complex, the replisome. Proofreading and error correction generally accompany replication but also occur independently, safeguarding genetic information through all phases of the cell cycle. Advances in biochemical characterization of intracellular processes, proteomics and the advent of single-molecule biophysics have brought about a treasure trove of information awaiting to be assembled into an accurate mechanistic model of the DNA replication process. In this review, we describe recent efforts to model elements of DNA replication and repair processes using computer simulations, an approach that has gained immense popularity in many areas of molecular biophysics but has yet to become mainstream in the DNA metabolism community. We highlight the use of diverse computational methods to address specific problems of the fields and discuss unexplored possibilities that lie ahead for the computational approaches in these areas.

Photochemical analysis of structural transitions in DNA liquid crystals reveals differences in spatial structure of DNA molecules organized in liquid crystalline form

The anisotropic shape of DNA molecules allows them to form lyotropic liquid crystals (LCs) at high concentrations. This liquid crystalline arrangement is also found in vivo (e.g., in bacteriophage capsids, bacteria or human sperm nuclei). However, the role of DNA liquid crystalline organization in living organisms still remains an open question. Here we show that in vitro, the DNA spatial structure is significantly changed in mesophases compared to non-organized DNA molecules. DNA LCs were prepared from pBluescript SK (pBSK) plasmid DNA and investigated by photochemical analysis of structural transitions (PhAST). We reveal significant differences in the probability of UV-induced pyrimidine dimer photoproduct formation at multiple loci on the DNA indicative of changes in major groove architecture.

Supercoiling DNA Locates Mismatches

We present a method of detecting sequence defects by supercoiling DNA with magnetic tweezers. The method is sensitive to a single mismatched base pair in a DNA sequence of several thousand base pairs. We systematically compare DNA molecules with 0 to 16 adjacent mismatches at 1 M monovalent salt and 3.6 pN force and show that under these conditions, a single plectoneme forms and is stably pinned at the defect. We use these measurements to estimate the energy and degree of end-loop kinking at defects. From this, we calculate the relative probability of plectoneme pinning at the mismatch under physiologically relevant conditions. Based on this estimate, we propose that DNA supercoiling could contribute to mismatch and damage sensing in vivo.

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.

DNA mutation motifs in the genes associated with inherited diseases

Mutations in human genes can be responsible for inherited genetic disorders and cancer. Mutations can arise due to environmental factors or spontaneously. It has been shown that certain DNA sequences are more prone to mutate. These sites are termed hotspots and exhibit a higher mutation frequency than expected by chance. In contrast, DNA sequences with lower mutation frequencies than expected by chance are termed coldspots. Mutation hotspots are usually derived from a mutation spectrum, which reflects particular population where an effect of a common ancestor plays a role. To detect coldspots/hotspots unaffected by population bias, we analysed the presence of germline mutations obtained from HGMD database in the 5-nucleotide segments repeatedly occurring in genes associated with common inherited disorders, in particular, the PAH, LDLR, CFTR, F8, and F9 genes. Statistically significant sequences (mutational motifs) rarely associated with mutations (coldspots) and frequently associated with mutations (hotspots) exhibited characteristic sequence patterns, e.g. coldspots contained purine tract while hotspots showed alternating purine-pyrimidine bases, often with the presence of CpG dinucleotide. Using molecular dynamics simulations and free energy calculations, we analysed the global bending properties of two selected coldspots and two hotspots with a G/T mismatch. We observed that the coldspots were inherently more flexible than the hotspots. We assume that this property might be critical for effective mismatch repair as DNA with a mutation recognized by MutSα protein is noticeably bent.

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.

Single-molecule FRET unveils induced-fit mechanism for substrate selectivity in flap endonuclease 1

Human flap endonuclease 1 (FEN1) and related structure-specific 5’nucleases precisely identify and incise aberrant DNA structures during replication, repair and recombination to avoid genomic instability. Yet, it is unclear how the 5’nuclease mechanisms of DNA distortion and protein ordering robustly mediate efficient and accurate substrate recognition and catalytic selectivity. Here, single-molecule sub-millisecond and millisecond analyses of FEN1 reveal a protein-DNA induced-fit mechanism that efficiently verifies substrate and suppresses off-target cleavage. FEN1 sculpts DNA with diffusion-limited kinetics to test DNA substrate. This DNA distortion mutually ‘locks’ protein and DNA conformation and enables substrate verification with extreme precision. Strikingly, FEN1 never misses cleavage of its cognate substrate while blocking probable formation of catalytically competent interactions with noncognate substrates and fostering their pre-incision dissociation. These findings establish FEN1 has practically perfect precision and that separate control of induced-fit substrate recognition sets up the catalytic selectivity of the nuclease active site for genome stability.

DOI:http://dx.doi.org/10.7554/eLife.21884.001

When a cell divides it must copy its genetic information, which is found in the form of strands of DNA. Damage to the DNA may lead to cancer or a number of genetic diseases. However, every time a cell divides more than 10 million toxic “flaps” of excess DNA are generated. A protein called flap endonuclease 1 (FEN1) keeps the DNA in good repair by cutting off the flaps in a highly specific and selective manner.

Many proteins that interact with DNA are attracted to specific genetic sequences within the DNA strands. However, this is not the case for FEN1 and several other “structure-specific” proteins that help to repair and replicate DNA strands. So how do these proteins select the correct regions of DNA to interact with?

Rashid et al. used single-molecule fluorescence measurements to examine how purified FEN1 proteins interact with DNA flaps. The results show that FEN1 can perfectly recognize and correctly remove flaps through a process called “mutual-induced fit”, where the DNA and FEN1 are shaped by each other to produce a highly specific structure.

Further work is now needed to examine whether other proteins that are related to FEN1 use a similar process to ensure that they always cut DNA in the same way. More detailed and direct examination of the structure of FEN1 through other experimental methods may also help to reveal how the mutual-induced fit process enables FEN1 to achieve such high levels of precision. This could increase our understanding of how problems with FEN1 and similar proteins lead to different genetic diseases.

DOI:http://dx.doi.org/10.7554/eLife.21884.002

Long-Range Signaling in MutS and MSH Homologs via Switching of Dynamic Communication Pathways

Allostery is conformation regulation by propagating a signal from one site to another distal site. This study focuses on the long-range communication in DNA mismatch repair proteins MutS and its homologs where intramolecular signaling has to travel over 70 Å to couple lesion detection to ATPase activity and eventual downstream repair. Using dynamic network analysis based on extensive molecular dynamics simulations, multiple preserved communication pathways were identified that would allow such long-range signaling. The pathways appear to depend on the nucleotides bound to the ATPase domain as well as the type of DNA substrate consistent with previously proposed functional cycles of mismatch recognition and repair initiation by MutS and homologs. A mechanism is proposed where pathways are switched without major conformational rearrangements allowing for efficient long-range signaling and allostery.

Anisotropy of B-DNA Groove Bending

DNA bending is critical for DNA packaging, recognition, and repair, and occurs toward either the major or the minor groove. The anisotropy of B-DNA groove bending was quantified for eight DNA sequences by free energy simulations employing a novel reaction coordinate. The simulations show that bending toward the major groove is preferred for non-A-tracts while the A-tract has a high tendency of bending toward the minor groove. Persistence lengths were generally larger for bending toward the minor groove, which is thought to originate from differences in groove hydration. While this difference in stiffness is one of the factors determining the overall preference of bending direction, the dominant contribution is shown to be a free energy offset between major and minor groove bending. The data suggests that, for the A-tract, this offset is largely determined by inherent structural properties, while differences in groove hydration play a large role for non-A-tracts. By quantifying the energetics of DNA groove bending and rationalizing the origins of the anisotropy, the calculations provide important new insights into a key biological process.

A novel conception for spontaneous transversions caused by homo-pyrimidine DNA mismatches: a QM/QTAIM highlight

We have firstly shown that the T·T(w) and C·C(w) DNA mismatches with wobble (w) geometry stay in slow tautomeric equilibrium with short T·T*(WC) and C·C*(WC) Watson-Crick (WC) mispairs. These non-dissociative tautomeric rearrangements are controlled by the plane-symmetric, highly stable, highly polar and zwitterionic transition states. The obtained results allow us to understand in what way the T·T(w) and C·C(w) mismatches acquire enzymatically competent T·T*(WC) and C·C*(WC) conformations directly in the hydrophobic recognition pocket of a high-fidelity DNA-polymerase, thereby producing thermodynamically non-equilibrium spontaneous transversions. The simplest numerical estimation of the frequency ratio of the TT to CC spontaneous transversions satisfactorily agrees with experimental data.

Dynamic DNA binding licenses a repair factor to bypass roadblocks in search of DNA lesions

DNA-binding proteins search for specific targets via facilitated diffusion along a crowded genome. However, little is known about how crowded DNA modulates facilitated diffusion and target recognition. Here we use DNA curtains and single-molecule fluorescence imaging to investigate how Msh2-Msh3, a eukaryotic mismatch repair complex, navigates on crowded DNA. Msh2-Msh3 hops over nucleosomes and other protein roadblocks, but maintains sufficient contact with DNA to recognize a single lesion. In contrast, Msh2-Msh6 slides without hopping and is largely blocked by protein roadblocks. Remarkably, the Msh3-specific mispair-binding domain (MBD) licences a chimeric Msh2-Msh6(3MBD) to bypass nucleosomes. Our studies contrast how Msh2-Msh3 and Msh2-Msh6 navigate on a crowded genome and suggest how Msh2-Msh3 locates DNA lesions outside of replication-coupled repair. These results also provide insights into how DNA repair factors search for DNA lesions in the context of chromatin.

A bouquet of DNA structures: Emerging diversity

Structural polymorphism of DNA has constantly been evolving from the time of illustration of the double helical model of DNA by Watson and Crick. A variety of non-canonical DNA structures have constantly been documented across the globe. DNA attracted worldwide attention as a carrier of genetic information. In addition to the classical Watson–Crick duplex, DNA can actually adopt diverse structures during its active participation in cellular processes like replication, transcription, recombination and repair. Structures like hairpin, cruciform, triplex, G-triplex, quadruplex, i-motif and other alternative non-canonical DNA structures have been studied at length and have also shown their in vivo occurrence. This review mainly focuses on non-canonical structures adopted by DNA oligonucleotides which have certain prerequisites for their formation in terms of sequence, its length, number and orientation of strands along with varied solution conditions. This conformational polymorphism of DNA might be the basis of different functional properties of a specific set of DNA sequences, further giving some insights for various extremely complicated biological phenomena. Many of these structures have already shown their linkages with diseases like cancer and genetic disorders, hence making them an extremely striking target for structure-specific drug designing and therapeutic applications.

•

DNA can adopt diverse range of structures other than classical Watson–Crick duplex.

•

Discussion of alternate structures like hairpin, cruciform, triplex, quadruplex etc.

•

This review gives some insights for the biological relevance of DNA structures.

Tolerance of DNA Mismatches in Dmc1 Recombinase-mediated DNA Strand Exchange

Recombination between homologous chromosomes is required for the faithful meiotic segregation of chromosomes and leads to the generation of genetic diversity. The conserved meiosis-specific Dmc1 recombinase catalyzes homologous recombination triggered by DNA double strand breaks through the exchange of parental DNA sequences. Although providing an efficient rate of DNA strand exchange between polymorphic alleles, Dmc1 must also guard against recombination between divergent sequences. How DNA mismatches affect Dmc1-mediated DNA strand exchange is not understood. We have used fluorescence resonance energy transfer to study the mechanism of Dmc1-mediated strand exchange between DNA oligonucleotides with different degrees of heterology. The efficiency of strand exchange is highly sensitive to the location, type, and distribution of mismatches. Mismatches near the 3' end of the initiating DNA strand have a small effect, whereas most mismatches near the 5' end impede strand exchange dramatically. The Hop2-Mnd1 protein complex stimulates Dmc1-catalyzed strand exchange on homologous DNA or containing a single mismatch. We observed that Dmc1 can reject divergent DNA sequences while bypassing a few mismatches in the DNA sequence. Our findings have important implications in understanding meiotic recombination. First, Dmc1 acts as an initial barrier for heterologous recombination, with the mismatch repair system providing a second level of proofreading, to ensure that ectopic sequences are not recombined. Second, Dmc1 stepping over infrequent mismatches is likely critical for allowing recombination between the polymorphic sequences of homologous chromosomes, thus contributing to gene conversion and genetic diversity.

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.

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.

Differential mismatch recognition specificities of eukaryotic MutS homologs, MutSα and MutSβ

In eukaryotes, the recognition of the DNA postreplication errors and initiation of the mismatch repair is carried out by two MutS homologs: MutSα and MutSβ. MutSα recognizes base mismatches and 1 to 2 unpaired nucleotides whereas MutSβ recognizes longer insertion-deletion loops (IDLs) with 1 to 15 unpaired nucleotides as well as certain mismatches. Results from molecular dynamics simulations of native MutSβ:IDL-containing DNA and MutSα:mismatch DNA complexes as well as complexes with swapped DNA substrates provide mechanistic insight into how the differential substrate specificities are achieved by MutSα and MutSβ, respectively. Our simulations results suggest more extensive interactions between MutSβ and IDL-DNA and between MutSα and mismatch-containing DNA that suggest corresponding differences in stability. Furthermore, our simulations suggest more expanded mechanistic details involving a different degree of bending when DNA is bound to either MutSα or MutSβ and a more likely opening of the clamp domains when noncognate substrates are bound. The simulation results also provide detailed information on key residues in MutSβ and MutSα that are likely involved in recognizing IDL-DNA and mismatch-containing DNA, respectively.

Coupled binding-bending-folding: The complex conformational dynamics of protein-DNA binding studied by atomistic molecular dynamics simulations

BACKGROUND

Protein-DNA binding often involves dramatic conformational changes such as protein folding and DNA bending. While thermodynamic aspects of this behavior are understood, and its biological function is often known, the mechanism by which the conformational changes occur is generally unclear. By providing detailed structural and energetic data, molecular dynamics simulations have been helpful in elucidating and rationalizing protein-DNA binding.

SCOPE OF REVIEW

This review will summarize recent atomistic molecular dynamics simulations of the conformational dynamics of DNA and protein-DNA binding. A brief overview of recent developments in DNA force fields is given as well.

MAJOR CONCLUSIONS

Simulations have been crucial in rationalizing the intrinsic flexibility of DNA, and have been instrumental in identifying the sequence of binding events, the triggers for the conformational motion, and the mechanism of binding for a number of important DNA-binding proteins.

GENERAL SIGNIFICANCE

Molecular dynamics simulations are an important tool for understanding the complex binding behavior of DNA-binding proteins. With recent advances in force fields and rapid increases in simulation time scales, simulations will become even more important for future studies. This article is part of a Special Issue entitled Recent developments of molecular dynamics.