Dynamic basis for dG•dT misincorporation via tautomerization and ionization

Kinetic Dissection of Proton-Coupled Conformational Transitions in Nucleic Acids by Integrating pH-Dependent NMR and Chemical Modifications

Proton-coupled conformational transitions play fundamental roles in nucleic acid recognition, catalysis, and folding, yet the kinetic mechanisms underlying these multistep protonation reactions remain unknown. Here, we present an approach to resolve the dominant kinetic pathway and rate-limiting step, which combines NMR chemical exchange measurements with chemical perturbations that shift pKa or modulate conformational equilibria. Applying the approach to three nucleic acid systems, we find the microscopic protonation step to be a diffusion-limited proton transfer reaction (kprot ∼ 1011 M-1 s-1), 2 orders of magnitude faster than diffusion-limited ligand-binding. For an A+-C mismatch in duplex DNA, protonation was the rate-limiting step occurring after the conformational change at a diffusion-limited kon ∼ 1011 M-1 s-1 via conformational selection of the wobble conformation, which forms rapidly and in significant abundance in the neutral ensemble. In RNA, the A-C wobble was sparsely populated in the neutral ensemble. The apparent kon was 2 orders of magnitude slower, and the reaction followed an induced-fit mechanism, where the unpaired adenine was initially protonated, followed by rate-limiting intrahelical flipping. The apparent kon was 5 orders of magnitude slower for the protonated G(syn)-C+ Hoogsteen conformation in duplex DNA in which cytosine protonation was rate-limiting occurring after the conformational change via conformational selection of an energetically disfavored G(syn)-C intermediate. These kinetic models quantitatively predicted the impact of pH shifts and chemical modifications on reaction kinetics. Our findings reveal how differences in nucleic acid conformational ensembles can drive diverse kinetic responses to pH changes and chemical modifications, even in binding reactions involving the simplest ligand: the proton.

Structure of a DNA Glycosylase Bound to a Nicked T:G Mismatch-Containing DNA

Mismatched T:G base pairs can arise during de novo replication as well as base excision repair (BER). In particular, the action of the gap-filling polymerase β (Polβ) can generate a T:G pair as well as a nick in the DNA backbone. The processing of a nicked T:G mispair is poorly understood. We are interested in understanding whether the T:G-specific DNA glycosylase MBD4 can recognize and process nicked T:G mismatches. We have discovered that MBD4 binds a nicked T:G-containing DNA, but does not cleave thymine opposite guanine. To gain insight into this, we have determined a crystal structure of human MBD4 bound to a nicked T:G-containing DNA. This structure displayed the full insertion of thymine into the catalytic site and the recognition of thymine based on the catalytic site's amino acid residues. However, thymine excision did not occur, presumably due to the inactivation of the catalytic D560 carboxylate nucleophile via a polar interaction with the 5'-hydrogen phosphate of the nicked DNA. The nicked complex was greatly stabilized by an ordered water molecule that formed four hydrogen bonds with the nicked DNA and MBD4. Interestingly, the arginine finger R468 did not engage in the phosphate pinching that is commonly observed in T:G mismatch recognition complex structures. Instead, the guanidinium moiety of R468 made bifurcated hydrogen bonding interactions with O6 of guanine, thereby stabilizing the estranged guanine. These observations suggest that R468 may sense and disrupt T:G pairs within the DNA duplex and stabilize the flipped-out thymine. The structure described here would be a close mimic of an intermediate in the base extrusion pathway induced by DNA glycosylase.

Insight into the Conformational Ensembles Formed by U-U and T-T Mismatches in RNA and DNA Duplexes From a Structure-based Survey, NMR, and Molecular Dynamics Simulations

Nucleic acid base pairs interconvert between alternative conformations on a free energy landscape, and these dynamics play critical roles in recognition, folding, and catalysis. U-U and T-T mismatches can adopt two nearly isoenergetic wobble conformations, distinguished by their relative shearing displacements. Experimental NMR evidence suggests that these conformations dynamically interconvert in RNA motifs containing tandem U-U mismatches. However, whether such motions occur ubiquitously across U-U and T-T mismatches remains unknown, as high-resolution nucleic acid structures typically report only a single conformation. Here, we used NMR spectroscopy, a structure-based survey of the Protein Data Bank, and molecular dynamics (MD) simulations to investigate wobble dynamics in U-U and T-T mismatches when flanked by canonical Watson-Crick base pairs in RNA and DNA duplexes. The structure-based survey revealed that U-U mismatches have propensities to adopt alternative wobble conformations even when controlling for sequence and identified potential intermediates along the wobble transition. Off-resonance R1ρ relaxation dispersion experiments detected no micro- to millisecond dynamics for U-U mismatches in duplex RNA and T-T mismatches in duplex DNA. However, alternative conformer refinement of the electron density in X-ray structures, inter-proton NOEs, carbonyl carbon chemical shifts, an RDC-derived conformational ensemble, and MD simulations indicated that U-U and T-T mismatches exist in a dynamic equilibrium between two wobble conformations, with the minor state exceeding 30% and the transitions occurring on the nanosecond timescale. Our findings suggest that U-U and T-T ubiquitously undergo sub-microsecond wobble motions, contributing to the energetic landscape and dynamic plasticity of nucleic acids, with important implications for processes that generate and act on these mismatches.

Quantum life science: biological nano quantum sensors, quantum technology-based hyperpolarized MRI/NMR, quantum biology, and quantum biotechnology

The emerging field of quantum life science combines principles from quantum physics and biology to study fundamental life processes at the molecular level. Quantum mechanics, which describes the properties of small particles, can help explain how quantum phenomena such as tunnelling, superposition, and entanglement may play a role in biological systems. However, capturing these effects in living systems is a formidable challenge, as it involves dealing with dissipation and decoherence caused by the surrounding environment. We overview the current status of the quantum life sciences from technologies and topics in quantum biology. Technologies such as biological nano quantum sensors, quantum technology-based hyperpolarized MRI/NMR, high-speed 2D electronic spectrometers, and computer simulations are being developed to address these challenges. These interdisciplinary fields have the potential to revolutionize our understanding of living organisms and lead to advancements in genetics, molecular biology, medicine, and bioengineering.

Ribose Sugar Alters Conformational Sampling of G⋅T Mismatched Duplex DNA

Polymerases erroneously incorporate Guanine-Thymine (dG⋅dT) mismatches in genomic DNA that further evades repair by transient sampling of tautomeric/ionic states compromising fidelity of repairing dG⋅dT mismatches. In conjunction, significant frequency of ribose (mis)incorporation in duplex DNA permits for misincorporated-mismatch in the genome. Ribose incorporated G (rG) mismatched with T (rG⋅dT) is the most stable across all misincorporated-mismatch calling into question the conformational consequences of the ribose sugar in addition to the mismatch. In this work, the effects of single rG⋅dT is investigated within a dodecamer DNA duplex employing solution-state NMR spectroscopy, partial anisotropic measurements in conjunction with molecular dynamics simulations to evaluate the impact on base pairs and the overall duplex structure. It is observed that rG⋅dT pairs exhibit enhanced flexibility in both base-pair and sugar dynamics compared to dG⋅dT, and the perturbations are enhanced in comparison to a ribose incorporated adenine-thymine (rA-dT) pair. The structural perturbations compared between rG⋅dT and dG⋅dT provides clues on plausible recognition modes of ribonucleotide excision repair (RER) pathway that looks for misincorporated ribose and mismatch repair (MMR) enzymes that scout for a mismatch.

Quantitative and Systematic NMR Measurements of Sequence-Dependent A-T Hoogsteen Dynamics in the DNA Double Helix

The dynamic properties of DNA depend on the sequence, providing an important source of sequence-specificity in biochemical reactions. However, comprehensively measuring how these dynamics vary with sequence is challenging, especially when they involve lowly populated and short-lived conformational states. Using 1H CEST supplemented by targeted 13C R1ρ NMR experiments, we quantitatively measured Watson-Crick to Hoogsteen dynamics for an A-T base pair in 13 trinucleotide sequence contexts. The Hoogsteen population and exchange rate varied 4-fold and 16-fold, respectively, and were dependent on both the 3'- and 5'-neighbors but only weakly dependent on monovalent ion concentration (25 versus 100 mM NaCl) and pH (6.8 versus 8.0). Flexible TA and CA dinucleotide steps exhibited the highest Hoogsteen populations, and their kinetics rates strongly depended on the 3'-neighbor. In contrast, the stiffer AA and GA steps had the lowest Hoogsteen population, and their kinetics were weakly dependent on the 3'-neighbor. The Hoogsteen lifetime was especially short when G-C neighbors flanked the A-T base pair. Our results uncover a unique conformational basis for sequence-specificity in the DNA double helix and establish the utility of NMR to quantitatively and comprehensively measure sequence-dependent DNA dynamics.

Structural basis of error-prone DNA synthesis by DNA polymerase θ

DNA polymerase θ (Pol θ) is an A-family DNA polymerase specialized in DNA double-strand breaks repair and translesion synthesis. Distinct from its high-fidelity homologs in DNA replication, Pol θ catalyzes template-dependent DNA synthesis with an inherent propensity for error incorporation. However, the structural basis of Pol θ's low-fidelity DNA synthesis is not clear. Here, we present cryo-electron microscopy structures detailing the polymerase domain of human Pol θ in complex with a cognate C:G base pair (bp), a mismatched T:G bp, or a mismatched T:T bp. Our structures illustrate that Pol θ snugly accommodates the mismatched nascent base pairs within its active site with the finger domain well-closed, consistent with our in-solution fluorescence measurement but in contrast to its high-fidelity homologs. In addition, structural examination and mutagenesis study show that unique residues surrounding the active site contribute to the stabilization of the mismatched nascent base pair. Furthermore, Pol θ can efficiently extend from the misincorporated T:G or T:T mismatches, yet with a preference for template or primer looping-out, resulting in insertions and deletions. Collectively, our results elucidate how an A-family polymerase is adapted for error-prone DNA synthesis.

Lifetime of ground conformational state determines the activity of structured RNA

Biomolecules continually sample alternative conformations. Consequently, even the most energetically favored ground conformational state has a finite lifetime. Here, we show that, in addition to the three-dimensional (3D) structure, the lifetime of a ground conformational state determines its biological activity. Using hydrogen-deuterium exchange nuclear magnetic resonance spectroscopy, we found that Zika virus exoribonuclease-resistant RNA (xrRNA) encodes a ground conformational state with a lifetime that is ~105-107 longer than that of canonical base pairs. Mutations that shorten the apparent lifetime of the ground state without affecting its 3D structure decreased exoribonuclease resistance in vitro and impaired virus replication in cells. Additionally, we observed this exceptionally long-lived ground state in xrRNAs from diverse infectious mosquito-borne flaviviruses. These results demonstrate the biological importance of the lifetime of a preorganized ground state and further suggest that elucidating the lifetimes of dominant 3D structures of biomolecules may be crucial for understanding their behaviors and functions.

A Single Multiplex PCR and Single-Nucleotide Extension Assay for the Detection of Common Thanatophoric Dysplasia I and II Mutations

Mutation analysis provides confirmation of a clinical and radiological diagnosis of thanatophoric dysplasia types I and II (TD I and II). We developed a single multiplexed PCR and a single-nucleotide extension (SNE) assay to identify 14 common mutations causing 99% of TD I and TD II, including the challenging three adjacent mutations in the stop codon of exon 18 of the FGFR3 gene. The assay design also provides a solution for resolving SNE PCR product sizing using performance optimized polymer-7. The assay was validated using 37 previously characterized, de-identified patient samples representing the nine wild-types and 10 of 14 mutant genotypes. Four artificial templates were synthesized to mimic four TD I mutations not represented in the available patient samples. Fragment size and fluorophore channel for each SNE product from 10 samples and the four artificial templates were used to define bins and panels for analysis with GeneMarker version 3.0 and GeneMapper version 6.0 software. Allele calls (bin placement within the panels) were verified using the remaining 27 previously characterized samples. This TD I and II PCR and SNE assay is a robust multiplexed assay, streamlined, to identify 14 mutations in one single reaction. This assay has a shorter turnaround time in comparison to traditional Sanger or next-generation sequencing.

Probing the nucleobase selectivity of RNA polymerases with dual-coding substrates

Formycin A (FOR) and pyrazofurin A (PYR) are nucleoside analogs with antiviral and antitumor properties. They are known to interfere with nucleic acid metabolism, but their direct effect on transcription is less understood. We explored how RNA polymerases (RNAPs) from bacteria, mitochondria, and viruses utilize FOR, PYR, and oxidized purine nucleotides. All tested polymerases incorporated FOR in place of adenine and PYR in place of uridine. FOR also exhibited surprising dual-coding behavior, functioning as a cytosine substitute, particularly for viral RNAP. In contrast, 8-oxoadenine and 8-oxoguanine were incorporated in place of uridine in addition to their canonical Watson-Crick codings. Our data suggest that the interconversion of canonical anti and alternative syn conformers underlies dual-coding abilities of FOR and oxidized purines. Structurally distinct RNAPs displayed varying abilities to utilize syn conformers during transcription. By examining base pairings that led to substrate incorporation and the entire spectrum of geometrically compatible pairings, we have gained new insights into the nucleobase selection processes employed by structurally diverse RNAPs. These insights may pave the way for advancements in antiviral therapies.

A review on point mutations via proton transfer in DNA base pairs in the absence and presence of electric fields

This comprehensive review focuses on spontaneous mutations that may occur during DNA replication, the fundamental process responsible for transferring genetic information. In 1963, Löwdin postulated that these mutations are primarily a result of proton transfer reactions within the hydrogen-bonded DNA base pairs. The single and double proton transfer reactions within the base pairs in DNA result in zwitterions and rare tautomers, respectively. For persistent mutations, these products must be generated at high rates and should be thermodynamically stable. This review covers the proton transfer reactions studied experimentally and computationally. The review also examines the influence of externally applied electric fields on the thermodynamics and kinetics of proton transfer reactions within DNA base pairs, and their biological implications.

Insights into the A-C Mismatch Conformational Ensemble in Duplex DNA and its Role in Genetic Processes through a Structure-based Review

Knowing the conformational ensembles formed by mismatches is crucial for understanding how they are generated and repaired and how they contribute to genomic instability. Here, we review structural and energetic studies of the A-C mismatch in duplex DNA and use the information to identify critical conformational states in its ensemble and their significance in genetic processes. In the 1970s, Topal and Fresco proposed the A-C wobble stabilized by two hydrogen bonds, one requiring protonation of adenine-N1. Subsequent NMR and X-ray crystallography studies showed that the protonated A-C wobble was in dynamic equilibrium with a neutral inverted wobble. The mismatch was shown to destabilize duplex DNA in a sequence- and pH-dependent manner by 2.4-3.8 kcal/mol and to have an apparent pKa ranging between 7.2 and 7.7. The A-C mismatch conformational repertoire expanded as structures were determined for damaged and protein-bound DNA. These structures included Watson-Crick-like conformations forming through tautomerization of the bases that drive replication errors, the reverse wobble forming through rotation of the entire nucleotide proposed to increase the fidelity of DNA replication, and the Hoogsteen base-pair forming through the flipping of the adenine base which explained the unusual specificity of DNA polymerases that bypass DNA damage. Thus, the A-C mismatch ensemble encompasses various conformational states that can be selectively stabilized in response to environmental changes such as pH shifts, intermolecular interactions, and chemical modifications, and these adaptations facilitate critical biological processes. This review also highlights the utility of existing 3D structures to build ensemble models for nucleic acid motifs.

Amber free energy tools: Interoperable software for free energy simulations using generalized quantum mechanical/molecular mechanical and machine learning potentials

We report the development and testing of new integrated cyberinfrastructure for performing free energy simulations with generalized hybrid quantum mechanical/molecular mechanical (QM/MM) and machine learning potentials (MLPs) in Amber. The Sander molecular dynamics program has been extended to leverage fast, density-functional tight-binding models implemented in the DFTB+ and xTB packages, and an interface to the DeePMD-kit software enables the use of MLPs. The software is integrated through application program interfaces that circumvent the need to perform "system calls" and enable the incorporation of long-range Ewald electrostatics into the external software's self-consistent field procedure. The infrastructure provides access to QM/MM models that may serve as the foundation for QM/MM-ΔMLP potentials, which supplement the semiempirical QM/MM model with a MLP correction trained to reproduce ab initio QM/MM energies and forces. Efficient optimization of minimum free energy pathways is enabled through a new surface-accelerated finite-temperature string method implemented in the FE-ToolKit package. Furthermore, we interfaced Sander with the i-PI software by implementing the socket communication protocol used in the i-PI client-server model. The new interface with i-PI allows for the treatment of nuclear quantum effects with semiempirical QM/MM-ΔMLP models. The modular interoperable software is demonstrated on proton transfer reactions in guanine-thymine mispairs in a B-form deoxyribonucleic acid helix. The current work represents a considerable advance in the development of modular software for performing free energy simulations of chemical reactions that are important in a wide range of applications.

Electronic and Nuclear Quantum Effects on Proton Transfer Reactions of Guanine-Thymine (G-T) Mispairs Using Combined Quantum Mechanical/Molecular Mechanical and Machine Learning Potentials

Rare tautomeric forms of nucleobases can lead to Watson-Crick-like (WC-like) mispairs in DNA, but the process of proton transfer is fast and difficult to detect experimentally. NMR studies show evidence for the existence of short-time WC-like guanine-thymine (G-T) mispairs; however, the mechanism of proton transfer and the degree to which nuclear quantum effects play a role are unclear. We use a B-DNA helix exhibiting a wGT mispair as a model system to study tautomerization reactions. We perform ab initio (PBE0/6-31G*) quantum mechanical/molecular mechanical (QM/MM) simulations to examine the free energy surface for tautomerization. We demonstrate that while the ab initio QM/MM simulations are accurate, considerable sampling is required to achieve high precision in the free energy barriers. To address this problem, we develop a QM/MM machine learning potential correction (QM/MM-ΔMLP) that is able to improve the computational efficiency, greatly extend the accessible time scales of the simulations, and enable practical application of path integral molecular dynamics to examine nuclear quantum effects. We find that the inclusion of nuclear quantum effects has only a modest effect on the mechanistic pathway but leads to a considerable lowering of the free energy barrier for the GT*⇌G*T equilibrium. Our results enable a rationalization of observed experimental data and the prediction of populations of rare tautomeric forms of nucleobases and rates of their interconversion in B-DNA.

Correction of non-random mutational biases along a linear bacterial chromosome by the mismatch repair endonuclease NucS

The linear chromosome of Streptomyces exhibits a highly compartmentalized structure with a conserved central region flanked by variable arms. As double strand break (DSB) repair mechanisms play a crucial role in shaping the genome plasticity of Streptomyces, we investigated the role of EndoMS/NucS, a recently characterized endonuclease involved in a non-canonical mismatch repair (MMR) mechanism in archaea and actinobacteria, that singularly corrects mismatches by creating a DSB. We showed that Streptomyces mutants lacking NucS display a marked colonial phenotype and a drastic increase in spontaneous mutation rate. In vitro biochemical assays revealed that NucS cooperates with the replication clamp to efficiently cleave G/T, G/G and T/T mismatched DNA by producing DSBs. These findings are consistent with the transition-shifted mutational spectrum observed in the mutant strains and reveal that NucS-dependent MMR specific task is to eliminate G/T mismatches generated by the DNA polymerase during replication. Interestingly, our data unveil a crescent-shaped distribution of the transition frequency from the replication origin towards the chromosomal ends, shedding light on a possible link between NucS-mediated DSBs and Streptomyces genome evolution.

Quantitative and systematic NMR measurements of sequence-dependent A-T Hoogsteen dynamics uncovers unique conformational specificity in the DNA double helix

The propensities to form lowly-populated short-lived conformations of DNA could vary with sequence, providing an important source of sequence-specificity in biochemical reactions. However, comprehensively measuring how these dynamics vary with sequence is challenging. Using 1H CEST and 13C R 1 ρ NMR, we measured Watson-Crick to Hoogsteen dynamics for an A-T base pair in thirteen trinucleotide sequence contexts. The Hoogsteen population and exchange rate varied 4-fold and 16-fold, respectively, and were dependent on both the 3'- and 5'-neighbors but only weakly dependent on monovalent ion concentration (25 versus 100 mM NaCl) and pH (6.8 versus 8.0). Flexible TA and CA dinucleotide steps exhibited the highest Hoogsteen populations, and their kinetics rates strongly depended on the 3'-neighbor. In contrast, the stiffer AA and GA steps had the lowest Hoogsteen population, and their kinetics were weakly dependent on the 3'-neighbor. The Hoogsteen lifetime was especially short when G-C neighbors flanked the A-T base pair. The Hoogsteen dynamics had a distinct sequence-dependence compared to duplex stability and minor groove width. Thus, our results uncover a unique source of sequence-specificity hidden within the DNA double helix in the form of A-T Hoogsteen dynamics and establish the utility of 1H CEST to quantitively measure sequence-dependent DNA dynamics.

NMR measurements of transient low-populated tautomeric and anionic Watson-Crick-like G·T/U in RNA:DNA hybrids: implications for the fidelity of transcription and CRISPR/Cas9 gene editing

Many biochemical processes use the Watson-Crick geometry to distinguish correct from incorrect base pairing. However, on rare occasions, mismatches such as G·T/U can transiently adopt Watson-Crick-like conformations through tautomerization or ionization of the bases, giving rise to replicative and translational errors. The propensities to form Watson-Crick-like mismatches in RNA:DNA hybrids remain unknown, making it unclear whether they can also contribute to errors during processes such as transcription and CRISPR/Cas editing. Here, using NMR R1ρ experiments, we show that dG·rU and dT·rG mismatches in two RNA:DNA hybrids transiently form tautomeric (Genol·T/U $ \mathbin{\lower.3ex\hbox{$\buildrel\textstyle\rightarrow\over {\smash{\leftarrow}\vphantom{_{\vbox to.5ex{\vss}}}}$}}$ G·Tenol/Uenol) and anionic (G·T-/U-) Watson-Crick-like conformations. The tautomerization dynamics were like those measured in A-RNA and B-DNA duplexes. However, anionic dG·rU- formed with a ten-fold higher propensity relative to dT-·rG and dG·dT- and this could be attributed to the lower pKa (ΔpKa ∼0.4-0.9) of U versus T. Our findings suggest plausible roles for Watson-Crick-like G·T/U mismatches in transcriptional errors and CRISPR/Cas9 off-target gene editing, uncover a crucial difference between the chemical dynamics of G·U versus G·T, and indicate that anionic Watson-Crick-like G·U- could play a significant role evading Watson-Crick fidelity checkpoints in RNA:DNA hybrids and RNA duplexes.

Surface-Accelerated String Method for Locating Minimum Free Energy Paths

We present a surface-accelerated string method (SASM) to efficiently optimize low-dimensional reaction pathways from the sampling performed with expensive quantum mechanical/molecular mechanical (QM/MM) Hamiltonians. The SASM accelerates the convergence of the path using the aggregate sampling obtained from the current and previous string iterations, whereas approaches like the string method in collective variables (SMCV) or the modified string method in collective variables (MSMCV) update the path only from the sampling obtained from the current iteration. Furthermore, the SASM decouples the number of images used to perform sampling from the number of synthetic images used to represent the path. The path is optimized on the current best estimate of the free energy surface obtained from all available sampling, and the proposed set of new simulations is not restricted to being located along the optimized path. Instead, the umbrella potential placement is chosen to extend the range of the free energy surface and improve the quality of the free energy estimates near the path. In this manner, the SASM is shown to improve the exploration for a minimum free energy pathway in regions where the free energy surface is relatively flat. Furthermore, it improves the quality of the free energy profile when the string is discretized with too few images. We compare the SASM, SMCV, and MSMCV using 3 QM/MM applications: a ribozyme methyltransferase reaction using 2 reaction coordinates, the 2'-O-transphosphorylation reaction of Hammerhead ribozyme using 3 reaction coordinates, and a tautomeric reaction in B-DNA using 5 reaction coordinates. We show that SASM converges the paths using roughly 3 times less sampling than the SMCV and MSMCV methods. All three algorithms have been implemented in the FE-ToolKit package made freely available.

Proton-polarized states in DNA

We observed signatures of a phase transition in the double-stranded DNA fragment of known length and sequences using a non-invasive semiconductor-electrolyte interface technique and statistical physics methods. Observations revealed a coherence peak in the electromotive force and a significant decline in calculated dynamic entropy at a critical temperature and pH. This behavior may arise from the dynamic interaction of proton (H+) pairs with opposite momentum and spin, carrying a charge q=2+ under critical conditions.

Direct Measurement of 8OG syn-anti Flips in Mutagenic 8OG•A and Long-Range Damage-Dependent Hoogsteen Breathing Dynamics Using 1H CEST NMR

Elucidating how damage impacts DNA dynamics is essential for understanding the mechanisms of damage recognition and repair. Many DNA lesions alter the propensities to form lowly-populated and short-lived conformational states. However, NMR methods to measure these dynamics require isotopic enrichment, which is difficult for damaged nucleotides. Here, we demonstrate the utility of the 1H chemical exchange saturation transfer (CEST) NMR experiment in measuring the dynamics of oxidatively damaged 8-oxoguanine (8OG) in the mutagenic 8OGsyn•Aanti mismatch. Using 8OG-H7 as an NMR probe of the damaged base, we directly measured 8OG syn-anti flips to form a lowly-populated (pop. ~ 5%) and short-lived (lifetime ~ 50 ms) non-mutagenic 8OGanti•Aanti. These exchange parameters were in quantitative agreement with values from 13C off-resonance R1ρ and CEST on a labeled partner adenine. The Watson-Crick-like 8OGsyn•Aanti mismatch also rescued the kinetics of Hoogsteen motions at distance A-T base pairs, which the G•A mismatch had slowed down. The results lend further support for 8OGanti•Aanti as a minor conformational state of 8OG•A, reveal that 8OG damage can impact Hoogsteen dynamics at a distance, and demonstrate the utility of 1H CEST for measuring damage-dependent dynamics in unlabeled DNA.

Molecular insight into how the position of an abasic site modifies DNA duplex stability and dynamics

Local perturbations to DNA base-pairing stability from lesions and chemical modifications can alter the stability and dynamics of an entire oligonucleotide. End effects may cause the position of a disruption within a short duplex to influence duplex stability and structural dynamics, yet this aspect of nucleic acid modifications is often overlooked. We investigate how the position of an abasic site (AP site) impacts the stability and dynamics of short DNA duplexes. Using a combination of steady-state and time-resolved spectroscopy and molecular dynamics simulations, we unravel an interplay between AP-site position and nucleobase sequence that controls energetic and dynamic disruption to the duplex. The duplex is disrupted into two segments by an entropic barrier for base-pairing on each side of the AP site. The barrier induces fraying of the short segment when an AP site is near the termini. Shifting the AP site inward promotes a transition from short-segment fraying to fully encompassing the barrier into the thermodynamics of hybridization, leading to further destabilization of the duplex. Nucleobase sequence determines the length scale for this transition by tuning the barrier height and base-pair stability of the short segment, and certain sequences enable out-of-register base-pairing to minimize the barrier height.

A unified Watson-Crick geometry drives transcription of six-letter expanded DNA alphabets by E. coli RNA polymerase

Artificially Expanded Genetic Information Systems (AEGIS) add independently replicable unnatural nucleotide pairs to the natural G:C and A:T/U pairs found in native DNA, joining the unnatural pairs through alternative modes of hydrogen bonding. Whether and how AEGIS pairs are recognized and processed by multi-subunit cellular RNA polymerases (RNAPs) remains unknown. Here, we show that E. coli RNAP selectively recognizes unnatural nucleobases in a six-letter expanded genetic system. High-resolution cryo-EM structures of three RNAP elongation complexes containing template-substrate UBPs reveal the shared principles behind the recognition of AEGIS and natural base pairs. In these structures, RNAPs are captured in an active state, poised to perform the chemistry step. At this point, the unnatural base pair adopts a Watson-Crick geometry, and the trigger loop is folded into an active conformation, indicating that the mechanistic principles underlying recognition and incorporation of natural base pairs also apply to AEGIS unnatural base pairs. These data validate the design philosophy of AEGIS unnatural basepairs. Further, we provide structural evidence supporting a long-standing hypothesis that pair mismatch during transcription occurs via tautomerization. Together, our work highlights the importance of Watson-Crick complementarity underlying the design principles of AEGIS base pair recognition.

Probing excited state 1Hα chemical shifts in intrinsically disordered proteins with a triple resonance-based CEST experiment: Application to a disorder-to-order switch

Over 40% of eukaryotic proteomes and 15% of bacterial proteomes are predicted to be intrinsically disordered based on their amino acid sequence. Intrinsically disordered proteins (IDPs) exist as heterogeneous ensembles of interconverting conformations and pose a challenge to the structure-function paradigm by apparently functioning without possessing stable structural elements. IDPs play a prominent role in biological processes involving extensive intermolecular interaction networks and their inherently dynamic nature facilitates their promiscuous interaction with multiple structurally diverse partner molecules. NMR spectroscopy has made pivotal contributions to our understanding of IDPs because of its unique ability to characterize heterogeneity at atomic resolution. NMR methods such as Chemical Exchange Saturation Transfer (CEST) and relaxation dispersion have enabled the detection of 'invisible' excited states in biomolecules which are transiently and sparsely populated, yet central for function. Here, we develop a 1Hα CEST pulse sequence which overcomes the resonance overlap problem in the 1Hα-13Cα plane of IDPs by taking advantage of the superior resolution in the 1H-15N correlation spectrum. In this sequence, magnetization is transferred after 1H CEST using a triple resonance coherence transfer pathway from 1Hα (i) to 1HN(i + 1) during which the 15N(t1) and 1HN(t2) are frequency labelled. This approach is integrated with spin state-selective CEST for eliminating spurious dips in CEST profiles resulting from dipolar cross-relaxation. We apply this sequence to determine the excited state 1Hα chemical shifts of the intrinsically disordered DNA binding domain (CytRN) of the bacterial cytidine repressor (CytR), which transiently acquires a functional globally folded conformation. The structure of the excited state, calculated using 1Hα chemical shifts in conjunction with other excited state NMR restraints, is a three-helix bundle incorporating a helix-turn-helix motif that is vital for binding DNA.

A quantum physics layer of epigenetics: a hypothesis deduced from charge transfer and chirality-induced spin selectivity of DNA

BACKGROUND

Epigenetic mechanisms are informational cellular processes instructing normal and diseased phenotypes. They are associated with DNA but without altering the DNA sequence. Whereas chemical processes like DNA methylation or histone modifications are well-accepted epigenetic mechanisms, we herein propose the existence of an additional quantum physics layer of epigenetics.

RESULTS

We base our hypothesis on theoretical and experimental studies showing quantum phenomena to be active in double-stranded DNA, even under ambient conditions. These phenomena include coherent charge transfer along overlapping pi-orbitals of DNA bases and chirality-induced spin selectivity. Charge transfer via quantum tunneling mediated by overlapping orbitals results in charge delocalization along several neighboring bases, which can even be extended by classical (non-quantum) electron hopping. Such charge transfer is interrupted by flipping base(s) out of the double-strand e.g., by DNA modifying enzymes. Charge delocalization can directly alter DNA recognition by proteins or indirectly by DNA structural changes e.g., kinking. Regarding sequence dependency, charge localization, shown to favor guanines, could influence or even direct epigenetic changes, e.g., modification of cytosines in CpG dinucleotides. Chirality-induced spin selectivity filters electrons for their spin along DNA and, thus, is not only an indicator for quantum coherence but can potentially affect DNA binding properties.

CONCLUSIONS

Quantum effects in DNA are prone to triggering and manipulation by external means. By the hypothesis put forward here, we would like to foster research on "Quantum Epigenetics" at the interface of medicine, biology, biochemistry, and physics to investigate the potential epigenetic impact of quantum physical principles on (human) life.

NMR measurements of transient low-populated tautomeric and anionic Watson-Crick-like G·T/U in RNA:DNA hybrids: Implications for the fidelity of transcription and CRISPR/Cas9 gene editing

Many biochemical processes use the Watson-Crick geometry to distinguish correct from incorrect base pairing. However, on rare occasions, mismatches such as G•T/U can transiently adopt Watson-Crick-like conformations through tautomerization or ionization of the bases, giving rise to replicative and translational errors. The propensities to form Watson-Crick-like mismatches in RNA:DNA hybrids remain unknown, making it unclear whether they can also contribute to errors during processes such as transcription and CRISPR/Cas editing. Here, using NMR R 1ρ experiments, we show that dG•rU and dT•rG mismatches in two RNA:DNA hybrids transiently form tautomeric (G enol •T/U ⇄G•T enol /U enol ) and anionic (G•T - /U - ) Watson-Crick-like conformations. The tautomerization dynamics were like those measured in A-RNA and B-DNA duplexes. However, anionic dG•rU - formed with a ten-fold higher propensity relative to dT - •rG and dG•dT - and this could be attributed to the lower pK a (Δ pK a ∼0.4-0.9) of U versus T. Our findings suggest plausible roles for Watson-Crick-like G•T/U mismatches in transcriptional errors and CRISPR/Cas9 off-target gene editing, uncover a crucial difference between the chemical dynamics of G•U versus G•T, and indicate that anionic Watson-Crick-like G•U - could play a significant role evading Watson-Crick fidelity checkpoints in RNA:DNA hybrids and RNA duplexes.

Molecular insight into how the position of an abasic site and its sequence environment influence DNA duplex stability and dynamics

Local perturbations to DNA base-pairing stability from lesions and chemical modifications can alter the stability and dynamics of an entire oligonucleotide. End effects may cause the position of a disruption within a short duplex to influence duplex stability and structural dynamics, yet this aspect of nucleic acid modifications is often overlooked. We investigate how the position of an abasic site (AP site) impacts the stability and dynamics of short DNA duplexes. Using a combination of steady-state and time-resolved spectroscopy and molecular dynamics simulations, we unravel an interplay between AP-site position and nucleobase sequence that controls energetic and dynamic disruption to the duplex. The duplex is disrupted into two segments by an entropic barrier for base pairing on each side of the AP site. The barrier induces fraying of the short segment when an AP site is near the termini. Shifting the AP site inward promotes a transition from short-segment fraying to fully encompassing the barrier into the thermodynamics of hybridization, leading to further destabilization the duplex. Nucleobase sequence determines the length scale for this transition by tuning the barrier height and base-pair stability of the short segment, and certain sequences enable out-of-register base pairing to minimize the barrier height.

Dynamic basis for dA•dGTP and dA•d8OGTP misincorporation via Hoogsteen base pairs

Replicative errors contribute to the genetic diversity needed for evolution but in high frequency can lead to genomic instability. Here, we show that DNA dynamics determine the frequency of misincorporating the A•G mismatch, and altered dynamics explain the high frequency of 8-oxoguanine (8OG) A•8OG misincorporation. NMR measurements revealed that Aanti•Ganti (population (pop.) of >91%) transiently forms sparsely populated and short-lived Aanti+•Gsyn (pop. of ~2% and kex = kforward + kreverse of ~137 s-1) and Asyn•Ganti (pop. of ~6% and kex of ~2,200 s-1) Hoogsteen conformations. 8OG redistributed the ensemble, rendering Aanti•8OGsyn the dominant state. A kinetic model in which Aanti+•Gsyn is misincorporated quantitatively predicted the dA•dGTP misincorporation kinetics by human polymerase β, the pH dependence of misincorporation and the impact of the 8OG lesion. Thus, 8OG increases replicative errors relative to G because oxidation of guanine redistributes the ensemble in favor of the mutagenic Aanti•8OGsyn Hoogsteen state, which exists transiently and in low abundance in the A•G mismatch.

Turing, von Neumann, and the computational architecture of biological machines

In the mid-1930s, the English mathematician and logician Alan Turing invented an imaginary machine which could emulate the process of manipulating finite symbolic configurations by human computers. His machine launched the field of computer science and provided a foundation for the modern-day programmable computer. A decade later, building on Turing's machine, the American-Hungarian mathematician John von Neumann invented an imaginary self-reproducing machine capable of open-ended evolution. Through his machine, von Neumann answered one of the deepest questions in Biology: Why is it that all living organisms carry a self-description in the form of DNA? The story behind how two pioneers of computer science stumbled on the secret of life many years before the discovery of the DNA double helix is not well known, not even to biologists, and you will not find it in biology textbooks. Yet, the story is just as relevant today as it was eighty years ago: Turing and von Neumann left a blueprint for studying biological systems as if they were computing machines. This approach may hold the key to answering many remaining questions in Biology and could even lead to advances in computer science.

19F NMR Untersuchung des Konformationsaustauschs mehrerer Zustände im synthetischen Neomycin-bindenden Riboschalter

Der synthetische Neomycin‐bindende Riboschalter interagiert mit seinem Liganden Neomycin sowie mit den verwandten Antibiotika Ribostamycin und Paromomycin. Die Bindung dieser Aminoglykoside induziert sehr ähnliche Grundzustandsstrukturen in der RNA, allerdings kann nur Neomycin die Initiierung der Translation effizient unterdrücken. Der molekulare Ursprung dieser Unterschiede wurde auf Unterschiede in der Dynamik der Ligand‐Riboschalter‐Komplexe zurückgeführt. In diesem Artikel kombinieren wir fünf komplementäre fluorbasierte NMR‐Methoden, um die Dynamik der drei Riboschalter‐Komplexe im Sekunden‐ bis Mikrosekundenbereich genau zu quantifizieren. Unsere Daten offenbaren komplexe Austauschprozesse mit bis zu vier strukturell unterschiedlichen Zuständen. Wir interpretieren unsere Ergebnisse in einem Modell, das ein Zusammenspiel zwischen verschiedenen chemischen Gruppen in den Antibiotika und spezifischen Basen im Riboschalter zeigt. Allgemeiner unterstreichen unsere Daten das Potenzial von 19F NMR‐Methoden, komplexe Austauschprozesse mit mehreren angeregten Zuständen zu charakterisieren.

Multi-Site Conformational Exchange in the Synthetic Neomycin-Sensing Riboswitch Studied by 19 F NMR

The synthetic neomycin-sensing riboswitch interacts with its cognate ligand neomycin as well as with the related antibiotics ribostamycin and paromomycin. Binding of these aminoglycosides induces a very similar ground state structure in the RNA, however, only neomycin can efficiently repress translation initiation. The molecular origin of these differences has been traced back to differences in the dynamics of the ligand:riboswitch complexes. Here, we combine five complementary fluorine based NMR methods to accurately quantify seconds to microseconds dynamics in the three riboswitch complexes. Our data reveal complex exchange processes with up to four structurally different states. We interpret our findings in a model that shows an interplay between different chemical groups in the antibiotics and specific bases in the riboswitch. More generally, our data underscore the potential of 19 F NMR methods to characterize complex exchange processes with multiple excited states.

Polγ coordinates DNA synthesis and proofreading to ensure mitochondrial genome integrity

Accurate replication of mitochondrial DNA (mtDNA) by DNA polymerase γ (Polγ) is essential for maintaining cellular energy supplies, metabolism, and cell cycle control. To illustrate the structural mechanism for Polγ coordinating polymerase (pol) and exonuclease (exo) activities to ensure rapid and accurate DNA synthesis, we determined four cryo-EM structures of Polγ captured after accurate or erroneous incorporation to a resolution of 2.4-3.0 Å. The structures show that Polγ employs a dual-checkpoint mechanism to sense nucleotide misincorporation and initiate proofreading. The transition from replication to error editing is accompanied by increased dynamics in both DNA and enzyme, in which the polymerase relaxes its processivity and the primer-template DNA unwinds, rotates, and backtracks to shuttle the mismatch-containing primer terminus 32 Å to the exo site for editing. Our structural and functional studies also provide a foundation for analyses of Polγ mutation-induced human diseases and aging.

Measuring the Enthalpy of an Individual Hydrogen Bond in a DNA Duplex with Nucleobase Isotope Editing and Variable-Temperature Infrared Spectroscopy

The level of interest in probing the strength of noncovalent interactions in DNA duplexes is high, as these weak forces dictate the range of suprastructures the double helix adopts under different conditions, in turn directly impacting the biological functions and industrial applications of duplexes that require making and breaking them to access the genetic code. However, few experimental tools can measure these weak forces embedded within large biological suprastructures in the native solution environment. Here, we develop experimental methods for detecting the presence of a single noncovalent interaction [a hydrogen bond (H-bond)] within a large DNA duplex in solution and measure its formation enthalpy (ΔH_f). We report that introduction of a H-bond into the TC2═O group from the noncanonical nucleobase 2-aminopurine produces an expected decrease ∼10 ± 0.76 cm^-1 (from ∼1720 cm^-1 in Watson-Crick to ∼1710 cm^-1 in 2-aminopurine), which correlates with an enthalpy of ∼0.93 ± 0.066 kcal/mol for this interaction.

Complementary base pair interactions between different rare tautomers of the second-generation artificial genetic alphabets

The functionality of a semisynthetic DNA in the biological environment will depend on the base pair nature of its complementary base pairs. To understand this, base pair interactions between complementary bases of recently proposed eight second-generation artificial nucleobases are studied herein by considering their rare tautomeric conformations and a dispersion-corrected density functional theoretic method. It is found that the binding energies of two hydrogen-bonded complementary base pairs are more negative than those of the three hydrogen-bonded base pairs. However, as the former base pairs are endothermic, the semisynthetic duplex DNA would involve the latter base pairs.

Tautomerisation Mechanisms in the Adenine-Thymine Nucleobase Pair during DNA Strand Separation

The adenine-thymine tautomer (A*-T*) has previously been discounted as a spontaneous mutagenesis mechanism due to the energetic instability of the tautomeric configuration. We study the stability of A*-T* while the nucleobases undergo DNA strand separation. Our calculations indicate an increase in the stability of A*-T* as the DNA strands unzip and the hydrogen bonds between the bases stretch. Molecular Dynamics simulations reveal the time scales and dynamics of DNA strand separation and the statistical ensemble of opening angles present in a biological environment. Our results demonstrate that the unwinding of DNA, an inherently out-of-equilibrium process facilitated by helicase, will change the energy landscape of the adenine-thymine tautomerization reaction. We propose that DNA strand separation allows the stable tautomerization of adenine-thymine, providing a feasible pathway for genetic point mutations via proton transfer between the A-T bases.

Tunneling Effect in Proton Transfer: Transfer Matrix Approach

The transfer matrix (TM) method was applied to calculate the transmission probability (TP) for proton transfer reactions. The tunneling factors in the reaction rate constants were also evaluated using the TPs. To test this method, TPs for the Eckart potentials modeled as a guanine-cytosine base pair were calculated by the TM method and compared to TPs by the analytical solution. As a result, the errors in the TPs by the TM method were quite small. The tunneling factors for the guanine-thymine (G-T) and adenine-cytosine (A-C) mispair reactions were then evaluated by the TM method. A shoulder appears on each potential for these reactions [Odai, K.; Umesaki,K. J. Phys. Chem. A. 2021, 125, 8196-8204]. As a result, the shoulder for the G-T mispair reaction contributes significantly to the tunneling, while the shoulder for the A-C mispair reaction contributes little to the tunneling. These results are difficult to obtain with methods such as Wigner's tunneling factor.

Quantum Tunnelling Effects in the Guanine-Thymine Wobble Misincorporation via Tautomerism

The misincorporation of a noncomplementary DNA base in the polymerase active site is a critical source of replication errors that can lead to genetic mutations. In this work, we model the mechanism of wobble mispairing and the subsequent rate of misincorporation errors by coupling first-principles quantum chemistry calculations to an open quantum systems master equation. This methodology allows us to accurately calculate the proton transfer between bases, allowing the misincorporation and formation of mutagenic tautomeric forms of DNA bases. Our calculated rates of genetic error formation are in excellent agreement with experimental observations in DNA. Furthermore, our quantum mechanics/molecular mechanics model predicts the existence of a short-lived "tunnelling-ready" configuration along the wobble reaction pathway in the polymerase active site, dramatically increasing the rate of proton transfer by a hundredfold, demonstrating that quantum tunnelling plays a critical role in determining the transcription error frequency of the polymerase.

Mutagenesis induced by protonation of single-stranded DNA is linked to glycolytic sugar metabolism

Mutagenesis can be thought of as random, in the sense that the occurrence of each mutational event cannot be predicted with precision in space or time. However, when sufficiently large numbers of mutations are analyzed, recurrent patterns of base changes called mutational signatures can be identified. To date, some 60 single base substitution or SBS signatures have been derived from analysis of cancer genomics data. We recently reported that the ubiquitous signature SBS5 matches the pattern of single nucleotide polymorphisms (SNPs) in humans and has analogs in many species. Using a temperature-sensitive single-stranded DNA (ssDNA) mutation reporter system, we also showed that a similar mutational pattern in yeast is dependent on error-prone translesion DNA synthesis (TLS) and glycolytic sugar metabolism. Here, we further investigated mechanisms that are responsible for this form of mutagenesis in yeast. We first confirmed that excess sugar metabolism leads to increased mutation rate, which was detectable by fluctuation assay. Since glycolysis is known to produce excess protons, we then investigated the effects of experimental manipulations on pH and mutagenesis. We hypothesized that yeast metabolizing 8% glucose would produce more excess protons than cells metabolizing 2% glucose. Consistent with this, cells metabolizing 8% glucose had lower intracellular and extracellular pH values. Similarly, deletion of vma3 (encoding a vacuolar H⁺-ATPase subunit) increased mutagenesis. We also found that treating cells with edelfosine (which renders membranes more permeable, including to protons) or culturing in low pH media increased mutagenesis. Analysis of the mutational pattern attributable to 20 µM edelfosine treatment revealed similarity to the SBS5-like TLS- and glycolysis-dependant mutational patterns previously observed in ssDNA. Altogether, our results agree with multiple biochemical studies showing that protonation of nitrogenous bases can alter base pairing so as to stabilize some mispairs, and shed new light on a common form of intrinsic mutagenesis.

Optimizing frequency sampling in CEST experiments

For the past decade chemical exchange saturation transfer (CEST) experiments have been successfully applied to study exchange processes in biomolecules involving sparsely populated, transiently formed conformers. Initial implementations focused on extensive sampling of the CEST frequency domain, requiring significant measurement times. Here we show that the lengthy sampling schemes often used are not optimal and that reduced frequency sampling schedules can be developed without a priori knowledge of the exchange parameters, that only depend on the chosen B₁ field, and, to a lesser extent, on the intrinsic transverse relaxation rates of ground state spins. The reduced sampling approach described here can be used synergistically with other methods for reducing measurement times such as those that excite multiple frequencies in the CEST dimension simultaneously, or make use of non-uniform sampling of indirectly detected time domains, to further decrease measurement times. The proposed approach is validated by analysis of simulated and experimental datasets.

Laplace's demon in biology: Models of evolutionary prediction

Our ability to predict natural phenomena can be limited by incomplete information. This issue is exemplified by "Laplace's demon," an imaginary creature proposed in the 18th century, who knew everything about everything, and thus could predict the full nature of the universe forward or backward in time. Quantum mechanics, among other things, has cast doubt on the possibility of Laplace's demon in the full sense, but the idea still serves as a useful metaphor for thinking about the extent to which prediction is limited by incomplete information on deterministic processes versus random factors. Here, we use simple analytical models and computer simulations to illustrate how data limits can be captured in a Bayesian framework, and how they influence our ability to predict evolution. We show how uncertainty in measurements of natural selection, or low predictability of external environmental factors affecting selection, can greatly reduce predictive power, often swamping the influence of intrinsic randomness caused by genetic drift. Thus, more accurate knowledge concerning the causes and action of natural selection is key to improving prediction. Fortunately, our analyses and simulations show quantitatively that reasonable improvements in data quantity and quality can meaningfully increase predictability.

Intrinsic base substitution patterns in diverse species reveal links to cancer and metabolism

Analyses of large-scale cancer sequencing data have revealed that mutagenic processes can create distinctive patterns of base substitutions, called mutational signatures. Interestingly, mutational patterns resembling some of these signatures can also be observed in normal cells. To determine whether similar patterns exist more generally, we analyzed large data sets of genetic variation, including mutations from 7 model species and single nucleotide polymorphisms in 42 species, totaling >1.9 billion variants. We found that base substitution patterns for most species closely match single base substitution (SBS) mutational signature 5 in the Catalog of Somatic Mutations in Cancer (COSMIC) database. SBS5 is ubiquitous in cancers and also present in normal human cells, suggesting that similar patterns of genetic variation across so many species are likely due to conserved biochemistry. We investigated the mechanistic origins of the SBS5-like mutational pattern in Saccharomyces cerevisiae, and show that translesion DNA synthesis and sugar metabolism are directly linked to this form of mutagenesis. We propose that conserved metabolic processes in cells are coupled to continuous generation of genetic variants, which can be acted upon by selection to drive the evolution of biological entities.

Structural basis for Cas9 off-target activity

The target DNA specificity of the CRISPR-associated genome editor nuclease Cas9 is determined by complementarity to a 20-nucleotide segment in its guide RNA. However, Cas9 can bind and cleave partially complementary off-target sequences, which raises safety concerns for its use in clinical applications. Here, we report crystallographic structures of Cas9 bound to bona fide off-target substrates, revealing that off-target binding is enabled by a range of noncanonical base-pairing interactions within the guide:off-target heteroduplex. Off-target substrates containing single-nucleotide deletions relative to the guide RNA are accommodated by base skipping or multiple noncanonical base pairs rather than RNA bulge formation. Finally, PAM-distal mismatches result in duplex unpairing and induce a conformational change in the Cas9 REC lobe that perturbs its conformational activation. Together, these insights provide a structural rationale for the off-target activity of Cas9 and contribute to the improved rational design of guide RNAs and off-target prediction algorithms.

Analyses of mutational patterns induced by formaldehyde and acetaldehyde reveal similarity to a common mutational signature

Formaldehyde and acetaldehyde are reactive small molecules produced endogenously in cells as well as being environmental contaminants. Both of these small aldehydes are classified as human carcinogens, since they are known to damage DNA and exposure is linked to cancer incidence. However, the mutagenic properties of formaldehyde and acetaldehyde remain incompletely understood, at least in part because they are relatively weak mutagens. Here, we use a highly sensitive yeast genetic reporter system featuring controlled generation of long single-stranded DNA regions to show that both small aldehydes induced mutational patterns characterized by predominantly C/G → A/T, C/G → T/A, and T/A → C/G substitutions, each in similar proportions. We observed an excess of C/G → A/T transversions when compared to mock-treated controls. Many of these C/G → A/T transversions occurred at TC/GA motifs. Interestingly, the formaldehyde mutational pattern resembles single base substitution signature 40 from the Catalog of Somatic Mutations in Cancer. Single base substitution signature 40 is a mutational signature of unknown etiology. We also noted that acetaldehyde treatment caused an excess of deletion events longer than 4 bases while formaldehyde did not. This latter result could be another distinguishing feature between the mutational patterns of these simple aldehydes. These findings shed new light on the characteristics of 2 important, commonly occurring mutagens.

Measuring thermodynamic preferences to form non-native conformations in nucleic acids using ultraviolet melting

Thermodynamic preferences to form non-native conformations are crucial for understanding how nucleic acids fold and function. However, they are difficult to measure experimentally because this requires accurately determining the population of minor low-abundance (<10%) conformations in a sea of other conformations. Here, we show that melting experiments enable facile measurements of thermodynamic preferences to adopt nonnative conformations in DNA and RNA. The key to this "delta-melt" approach is to use chemical modifications to render specific minor non-native conformations the major state. The validity and robustness of delta-melt is established for four different non-native conformations under various physiological conditions and sequence contexts through independent measurements of thermodynamic preferences using NMR. Delta-melt is faster relative to NMR, simple, and cost-effective and enables thermodynamic preferences to be measured for exceptionally low-populated conformations. Using delta-melt, we obtained rare insights into conformational cooperativity, obtaining evidence for significant cooperativity (1.0 to 2.5 kcal/mol) when simultaneously forming two adjacent Hoogsteen base pairs. We also measured the thermodynamic preferences to form G-C⁺ and A-T Hoogsteen and A-T base open states for nearly all 16 trinucleotide sequence contexts and found distinct sequence-specific variations on the order of 2 to 3 kcal/mol. This rich landscape of sequence-specific non-native minor conformations in the DNA double helix may help shape the sequence specificity of DNA biochemistry. Thus, melting experiments can now be used to access thermodynamic information regarding regions of the free energy landscape of biomolecules beyond the native folded and unfolded conformations.

Watching right and wrong nucleotide insertion captures hidden polymerase fidelity checkpoints

Efficient and accurate DNA synthesis is enabled by DNA polymerase fidelity checkpoints that promote insertion of the right instead of wrong nucleotide. Erroneous X-family polymerase (pol) λ nucleotide insertion leads to genomic instability in double strand break and base-excision repair. Here, time-lapse crystallography captures intermediate catalytic states of pol λ undergoing right and wrong natural nucleotide insertion. The revealed nucleotide sensing mechanism responds to base pair geometry through active site deformation to regulate global polymerase-substrate complex alignment in support of distinct optimal (right) or suboptimal (wrong) reaction pathways. An induced fit during wrong but not right insertion, and associated metal, substrate, side chain and pyrophosphate reaction dynamics modulated nucleotide insertion. A third active site metal hastened right but not wrong insertion and was not essential for DNA synthesis. The previously hidden fidelity checkpoints uncovered reveal fundamental strategies of polymerase DNA repair synthesis in genomic instability.

DNA polymerase (pol) λ performs DNA synthesis in base excision and double strand break repair. How pol λ accomplishes nucleotide insertion that can lead to mutagenesis and genomic instability was unclear. Here the authors employ time-lapse crystallography to reveal hidden polymerase checkpoints that enable right and wrong natural nucleotide insertion by pol λ.

Isotope Labels Combined with Solution NMR Spectroscopy Make Visible the Invisible Conformations of Small-to-Large RNAs

RNA is central to the proper function of cellular processes important for life on earth and implicated in various medical dysfunctions. Yet, RNA structural biology lags significantly behind that of proteins, limiting mechanistic understanding of RNA chemical biology. Fortunately, solution NMR spectroscopy can probe the structural dynamics of RNA in solution at atomic resolution, opening the door to their functional understanding. However, NMR analysis of RNA, with only four unique ribonucleotide building blocks, suffers from spectral crowding and broad linewidths, especially as RNAs grow in size. One effective strategy to overcome these challenges is to introduce NMR-active stable isotopes into RNA. However, traditional uniform labeling methods introduce scalar and dipolar couplings that complicate the implementation and analysis of NMR measurements. This challenge can be circumvented with selective isotope labeling. In this review, we outline the development of labeling technologies and their application to study biologically relevant RNAs and their complexes ranging in size from 5 to 300 kDa by NMR spectroscopy.

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.