Watson-Crick Face Research Articles

WITHDRAWN: Detection of N3-methylcytidine in RNA at single base resolution using Protoscript II to read through and generate high mutation rates

N3-methylcytidine (m3C), an RNA modification involving methylation of cytidine at the Watson-Crick face, has been identified in human tRNA, mitochondrial tRNA (mt-tRNA), and mRNA. In tRNA and mt-tRNA, it has been shown to play important roles in tRNA structure and translation. m3C remains understudied in mRNA and other low-abundance RNAs due to its markedly lower levels, the lack of high-resolution techniques to reproducibly map its location, and under-characterization of its effector proteins on these RNA types. Here we describe a new approach which utilizes the high readthrough and mutation rates of Protoscript II (PSII) at m3C sites in combination with AlkB demethylation to detect m3C at single-base resolution. We first validate the method on synthetic probes. Then, using this method, we successfully identify the known m3C sites on tRNA and mt-tRNA and show that none of the m3C sites identified in cytosolic tRNA are likely METTL8 substrates. In contrast, the two m3C sites detected in mt-tRNA are installed by METTL8 since they are sensitive to METTL8 knockdown (KD) and over expression (OE). These results demonstrate that the method can be used to investigate changes in the m3C level in many sites using an unbiased approach. Overall, this study introduces a convenient high-resolution method for m3C mapping with high readthrough rates, paving the way for future mechanistic studies on the biological function of m3C and its effector proteins.

A new reagent for in vivo structure probing of RNA G and U residues that improves RNA structure prediction alone and combined with DMS.

A key to understanding the roles of RNA in regulating gene expression is knowing their structures in vivo. One way to obtain this information is through probing the structures of RNA with chemicals. To probe RNA structure directly in cells, membrane-permeable reagents that modify the Watson-Crick (WC) face of unpaired nucleotides can be used. Although dimethyl sulfate (DMS) has led to substantial insight into RNA structure, it has limited nucleotide specificity in vivo, with WC face reactivity only at adenine (A) and cytosine (C) at neutral pH. The reagent 1-ethyl-3-(3-dimethylaminopropyl) carbodiimide (EDC) was recently shown to modify the WC face of guanine (G) and uracil (U). Although useful at lower concentrations in experiments that measure chemical modifications by reverse transcription (RT) stops, at higher concentrations necessary for detection by mutational profiling (MaP), EDC treatment leads to degradation of RNA. Here, we demonstrate EDC-stimulated degradation of RNA in Gram-negative and Gram-positive bacteria. In an attempt to overcome these limitations, we developed a new carbodiimide reagent, 1-ethyl-3-(3-dimethylaminopropyl) carbodiimide methiodide (ETC), which we show specifically modifies unpaired Gs and Us in vivo without substantial degradation of RNA. We establish ETC as a probe for MaP and optimize the RT conditions and computational analysis in Escherichia coli Importantly, we demonstrate the utility of ETC as a probe for improving RNA structure prediction both alone and with DMS.

Replication Bypass of the N-(2-Deoxy-d-erythro-pentofuranosyl)-urea DNA Lesion by Human DNA Polymerase η.

Urea lesions in DNA arise from thymine glycol (Tg) or 8-oxo-dG; their genotoxicity is thought to arise in part due to their potential to accommodate the insertion of all four dNTPs during error-prone replication. Replication bypass with human DNA polymerase η (hPol η) confirmed that all four dNTPs were inserted opposite urea lesions but with purines exhibiting greater incorporation efficiency. X-ray crystal structures of ternary replication bypass complexes in the presence of Mg2+ ions with incoming dNTP analogs dAMPnPP, dCMPnPP, dGMPnPP, and dTMPnPP bound opposite urea lesions (hPol η·DNA·dNMPnPP complexes) revealed all were accommodated by hPol η. In each, the Watson-Crick face of the dNMPnPP was paired with the urea lesion, exploiting the ability of the amine and carbonyl groups of the urea to act as H-bond donors or acceptors, respectively. With incoming dAMPnPP or dGMPnPP, the distance between the imino nitrogen of urea and the N9 atoms of incoming dNMPnPP approximated the canonical distance of 9 Å in B-DNA. With incoming dCMPnPP or dTMPnPP, the corresponding distance of about 7 Å was less ideal. Improved base-stacking interactions were also observed with incoming purines vs pyrimidines. Nevertheless, in each instance, the α-phosphate of incoming dNMPnPPs was close to the 3'-hydroxyl group of the primer terminus, consistent with the catalysis of nucleotidyl transfer and the observation that all four nucleotides could be inserted opposite urea lesions. Preferential insertion of purines by hPol η may explain, in part, why the urea-directed spectrum of mutations arising from Tg vs 8-oxo-dG lesions differs.

Global RNA modifications to the MALAT1 triple helix differentially affect thermostability and weaken binding to METTL16

Therapeutic mRNAs are generated using modified nucleotides, namely N1-methylpseudouridine (m1Ψ) triphosphate, so that the mRNA evades detection by the immune system. RNA modifications, even at a single-nucleotide position, perturb RNA structure, although it is not well understood how structure and function is impacted by globally modified RNAs. Therefore, we examined the metastasis-associated lung adenocarcinoma transcript 1 triple helix, a highly structured stability element that includes single-, double-, and triple-stranded RNA, globally modified with N6-methyladenosine (m6A), pseudouridine (Ψ), or m1Ψ. UV thermal denaturation assays showed that m6A destabilizes both the Hoogsteen and Watson-Crick faces of the RNA by ∼20 °C, Ψ stabilizes the Hoogsteen and Watson-Crick faces of the RNA by ∼12 °C, and m1Ψ has minimal effect on the stability of the Hoogsteen face of the RNA but increases the stability of the Watson-Crick face by ∼9 °C. Native gel-shift assays revealed that binding of the methyltransferase-like protein 16 to the metastasis-associated lung adenocarcinoma transcript 1 triple helix was weakened by at least 8-, 99-, and 23-fold, respectively, when RNA is globally modified with m6A, Ψ, or m1Ψ. These results demonstrate that a more thermostable RNA structure does not lead to tighter RNA-protein interactions, thereby highlighting the regulatory power of RNA modifications by multiple means.

6-Pyrazolylpurine and its deaza derivatives as nucleobases for silver(I)-mediated base pairing with pyrimidines

The artificial nucleobase 6-pyrazolylpurine (6PP) and its deaza derivatives 1-deaza-6-pyrazolylpurine (1D6PP), 7-deaza-6-pyrazolylpurine (7D6PP), and 1,7-dideaza-6-pyrazolylpurine (1,7D6PP) were investigated with respect to their ability to differentiate between the canonical nucleobases cytosine and thymine by means of silver(I)-mediated base pairing. As shown by temperature-dependent UV spectroscopy and by circular dichroism spectroscopy, 6PP and (to a lesser extent) 7D6PP form stable silver(I)-mediated base pairs with cytosine, but not with thymine. 1D6PP and 1,7D6PP do not engage in the formation of stabilizing silver(I)-mediated base pairs with cytosine or thymine. The different behavior of 1D6PP, 7D6PP, and 1,7D6PP indicates that silver(I) binding occurs via the N1 position of the purine derivative, i.e. via the Watson–Crick face. The data show that 6PP is capable of differentiating between cytosine and thymine, which is potentially relevant in the context of detecting single-nucleotide polymorphisms.Graphical abstract

5-Formylcytosine weakens the G-C pair and imparts local conformational fluctuations to DNA duplexes.

DNA epigenetic modifications such as 5-methyl (5mC), 5-hydroxymethyl (5hmC), 5-formyl (5fC) and 5-carboxyl (5caC) cytosine have unique and specific biological roles. Crystallographic studies of 5mC containing duplexes were conducted in the A-, B- or the intermediate E-DNA polymorphic forms. 5fC-modified duplexes initially observed in the disputed F-DNA architecture were subsequently crystallized in the A-form, suggesting that epigenetic modifications enable DNA sequences to adopt diverse conformational states that plausibly contribute to their function. Solution-state studies of these modifications were found in the B-DNA form, with marked differences in the conformational flexibility of 5fC containing duplexes in comparison to C/5mC containing duplexes, compromising the DNA duplex's stability. Herein, we systematically evaluate sensitive and commonly inaccessible NMR parameters to map the subtle differences between C, 5mC, and their oxidized (5hmC/5fC) counterparts. We observe that 15N/1H chemical shifts effectively report on the weakening of 5fC-G Watson-Crick base-pair H-bonding, extending the instability beyond any achievable within the sequence-specific changes in DNA. Triple 5fC containing sequences propagate the destabilization farther from the site of modifications, explaining reduced duplex stability upon multiple modifications. Additionally, scalar and residual dipolar coupling measurements unravel local sugar pucker fluctuations. One-bond 13C-1H scalar coupling measurements point towards a significant deviation away from the anticipated C2'-endo pucker for the 5fC modified nucleotide. Structural models obtained employing 13C-1H residual dipolar couplings and inter-proton distances corroborate the sugar pucker's deviation for 5fC modified DNA duplexes. The changes in the sugar pucker equilibria remain local to the 5fC modified nucleotide sans additive/long-range effects arising from multiple contiguous modifications. These observations highlight the impact of a major groove modification that alters the physical properties of DNA duplex without disturbing the Watson-Crick face. The changes observed in our studies for the 5fC containing DNA contrast with the perturbations induced by damage/lesion highlight the varied conformational preferences that modified nucleobases impart to the DNA duplex. As sequence-specific DNA transactions are rooted in the base-pair stability and pucker deviations, the observed structural perturbations for 5fC-modified DNA potentially play critical functional roles, such as protein-DNA recognition and interactions.

Probing RNA Structure with Dimethyl Sulfate Mutational Profiling with Sequencing <em>In Vitro</em> and in Cells

The role of RNA structure in virtually any biological process has become increasingly evident, especially in the past decade. However, classical approaches to solving RNA structure, such as RNA crystallography or cryo-EM, have failed to keep up with the rapidly evolving field and the need for high-throughput solutions. Mutational profiling with sequencing using dimethyl sulfate (DMS) MaPseq is a sequencing-based approach to infer the RNA structure from a base's reactivity with DMS. DMS methylates the N1 nitrogen in adenosines and the N3 in cytosines at their Watson-Crick face when the base is unpaired. Reverse-transcribing the modified RNA with the thermostable group II intron reverse transcriptase (TGIRT-III) leads to the methylated bases being incorporated as mutations into the cDNA. When sequencing the resulting cDNA and mapping it back to a reference transcript, the relative mutation rates for each base are indicative of the base's "status" as paired or unpaired. Even though DMS reactivities have a high signal-to-noise ratio both in vitro and in cells, this method is sensitive to bias in the handling procedures. To reduce this bias, this paper provides a protocol for RNA treatment with DMS in cells and with in vitro transcribed RNA.

G4Atlas: a comprehensive transcriptome-wide G-quadruplex database.

RNA G-quadruplex (rG4) is a vital RNA tertiary structure motif that involves the base pairs on both Hoogsteen and Watson-Crick faces of guanines. rG4 is of great importance in the post-transcriptional regulation of gene expression. Experimental technologies have advanced to identify in vitro and in vivo rG4s across diverse transcriptomes. Building on these recent advances, here we present G4Atlas, the first transcriptome-wide G-quadruplex database, in which we have collated, classified, and visualized transcriptome rG4 experimental data, generated from rG4-seq, chemical profiling and ligand-binding methods. Our comprehensive database includes transcriptome-wide rG4s generated from 82 experimental treatments and 238 samples across ten species. In addition, we have also included RNA secondary structure prediction information across both experimentally identified and unidentified rG4s to enable users to display any potential competitive folding between rG4 and RNA secondary structures. As such, G4Atlas will enable users to explore the general functions of rG4s in diverse biological processes. In addition, G4Atlas lays the foundation for further data-driven deep learning algorithms to examine rG4 structural features.

Web-based platform for analysis of RNA folding from high throughput chemical probing data.

RNA structures play critical roles in regulating gene expression across all domains of life and viruses. Chemical probing methods coupled with massively parallel sequencing have revolutionized the RNA structure field by enabling the assessment of many structures in their native, physiological context. Previously, we developed Dimethyl-Sulfate-based Mutational Profiling and Sequencing (DMS-MaPseq), which uses DMS to label the Watson-Crick face of open and accessible adenine and cytosine bases in the RNA. We used this approach to determine the genome-wide structures of HIV-1 and SARS-CoV-2 in infected cells, which permitted uncovering new biology and identifying therapeutic targets. Due to the simplicity and ease of the experimental procedure, DMS-MaPseq has been adopted by labs worldwide. However, bioinformatic analysis remains a substantial hurdle for labs that often lack the necessary infrastructure and computational expertise. Here we present a modern web-based interface that automates the analysis of chemical probing profiles from raw sequencing files (http://rnadreem.org). The availability of a simple web-based platform for DMS-MaPseq analysis will dramatically expand studies of RNA structure and aid in the design of structure-based therapeutics.

DNA-Mimicking Metal-Organic Frameworks with Accessible Adenine Faces for Complementary Base Pairing.

Biologically derived metal–organic frameworks (Bio-MOFs) are significant, as they can be used in cutting-edge biomedical applications such as targeted gene delivery. Herein, adenine (Ade) and unnatural amino acids coordinate with Zn2+ to produce biocompatible frameworks, KBM-1 and KBM-2, with extremely defined porous channels. They feature an accessible Watson–Crick Ade face that is available for further hydrogen bonding and can load single-stranded DNA (ssDNA) with 13 and 41% efficiency for KBM-1 and KBM-2, respectively. Treatment of these frameworks with thymine (Thy), as a competitive guest for base pairing with the Ade open sites, led to more than 50% reduction of ssDNA loading. Moreover, KBM-2 loaded Thy-rich ssDNA more efficiently than Thy-free ssDNA. These findings support the role of the Thy-Ade base pairing in promoting ssDNA loading. Furthermore, theoretical calculations using the self-consistent charge density functional tight-binding (SCC-DFTB) method verified the role of hydrogen bonding and van der Waals type interactions in this host–guest interface. KBM-1 and KBM-2 can protect ssDNA from enzymatic degradation and release it at acidic pH. Most importantly, these biocompatible frameworks can efficiently deliver genetic cargo with retained activity to the cell nucleus. We envisage that this class of Bio-MOFs can find immediate applicability as biomimics for sensing, stabilizing, and delivering genetic materials.

Structural characterization of a new subclass of panicum mosaic virus-like 3' cap-independent translation enhancer.

Canonical eukaryotic mRNA translation requires 5′cap recognition by initiation factor 4E (eIF4E). In contrast, many positive-strand RNA virus genomes lack a 5′cap and promote translation by non-canonical mechanisms. Among plant viruses, PTEs are a major class of cap-independent translation enhancers located in/near the 3′UTR that recruit eIF4E to greatly enhance viral translation. Previous work proposed a single form of PTE characterized by a Y-shaped secondary structure with two terminal stem-loops (SL1 and SL2) atop a supporting stem containing a large, G-rich asymmetric loop that forms an essential pseudoknot (PK) involving C/U residues located between SL1 and SL2. We found that PTEs with less than three consecutive cytidylates available for PK formation have an upstream stem-loop that forms a kissing loop interaction with the apical loop of SL2, important for formation/stabilization of PK. PKs found in both subclasses of PTE assume a specific conformation with a hyperreactive guanylate (G*) in SHAPE structure probing, previously found critical for binding eIF4E. While PTE PKs were proposed to be formed by Watson–Crick base-pairing, alternative chemical probing and 3D modeling indicate that the Watson–Crick faces of G* and an adjacent guanylate have high solvent accessibilities. Thus, PTE PKs are likely composed primarily of non-canonical interactions.

Structural Rationalization for the Nonmutagenic and Mutagenic Bypass of the Tobacco-Derived O4-4-(3-Pyridyl)-4-oxobut-1-yl-thymine Lesion by Human Polymerase η: A Multiscale Computational Study.

Tobacco-derived pyridyloxobutyl (POB) DNA adducts are unique due to the large size and flexibility of the alkyl chain connecting the pyridyl ring to the nucleobase. Recent experimental work suggests that the O4-4-(3-pyridyl)-4-oxobut-1-yl-T (O4-POB-T) lesion can undergo both nonmutagenic (dATP) and mutagenic (dGTP) insertion by the translesion synthesis (TLS) polymerase (pol) η in human cells. Interestingly, the mutagenic rate for O4-POB-T replication is reduced compared to that for the smaller O4-methylthymine (O4-Me-T) lesion, and O4-POB-T yields a different mutagenic profile than the O2-POB-T variant (dTTP insertion). The present work uses a combination of density functional theory calculations and molecular dynamics simulations to probe the impact of the size and flexibility of O4-POB-T on pol η replication outcomes. Due to changes in the Watson-Crick binding face upon damage of canonical T, O4-POB-T does not form favorable hydrogen-bonding interactions with A. Nevertheless, dATP is positioned for insertion in the pol η active site by a water chain to the template strand, which suggests a pol η replication pathway similar to that for abasic sites. Although a favorable O4-POB-T:G mispair forms in the pol η active site and DNA duplexes, the inherent dynamical nature of O4-POB-T periodically disrupts interstrand hydrogen bonding that would otherwise facilitate dGTP insertion and stabilize damaged DNA duplexes. In addition to explaining the origin of the experimentally reported pol η outcomes associated with O4-POB-T replication, comparison to structural data for the O4-Me-T and O2-POB-T adducts highlights an emerging common pathway for the nonmutagenic replication of thymine alkylated lesions by pol η, yet underscores the broader impacts of bulky moiety size, flexibility, and position on the associated mutagenic outcomes.

Hoogsteen Base Pairs Increase the Susceptibility of Double-Stranded DNA to Cytotoxic Damage

Watson-Crick base pairing is believed to protect the Watson-Crick face of nucleobases in double-stranded DNA (dsDNA) from cytotoxic damage. Lesions on the Watson-Crick face are believed to occur when the DNA is transiently single stranded during replication and transcription. However, damage to the Watson-Crick face in dsDNA has been reported in vitro, but the mechanisms rendering dsDNA susceptible to damage remain poorly understood. in this work, we explored the role of Hoogsteen base pairs (bps) in mediating cytotoxic damage to dsDNA.

7-Connected FeIII3-Based Bio-MOF: Pore Space Partition and Gas Separations.

We reported herein a new 3D bio-MOF (NbU-12) using a pore space partition strategy: MIL-88D was selected as a primary framework, and adenine connected two independent MIL-88D to form a self-interpenetrated structure. Because of this, the hexagonal channel in MIL-88D split into two small rectangular channels. Different from the reported series CPM-35 materials, NbU-12 simultaneously maximized the retention of open metal sites from MIL-88D and introduced a Watson-Crick face to the pore surface of NbU-12. Remarkably, NbU-12 exhibits an excellent selectivity performance toward C2H6/C2H4 and C2H6/CH4, which was proven by ideal adsorbed solution theory calculation and breakthrough experiments.

Hoogsteen base pairs increase the susceptibility of double-stranded DNA to cytotoxic damage

As the Watson-Crick faces of nucleobases are protected in dsDNA, it is commonly assumed that deleterious alkylation damage to the Watson-Crick faces of nucleobases predominantly occurs when DNA becomes single-stranded during replication and transcription. However, damage to the Watson-Crick faces of nucleobases has been reported in dsDNA in vitro through mechanisms that are not understood. In addition, the extent of protection from methylation damage conferred by dsDNA relative to ssDNA has not been quantified. Watson-Crick base pairs in dsDNA exist in dynamic equilibrium with Hoogsteen base pairs that expose the Watson-Crick faces of purine nucleobases to solvent. Whether this can influence the damage susceptibility of dsDNA remains unknown. Using dot-blot and primer extension assays, we measured the susceptibility of adenine-N1 to methylation by dimethyl sulfate (DMS) when in an A-T Watson-Crick versus Hoogsteen conformation. Relative to unpaired adenines in a bulge, Watson-Crick A-T base pairs in dsDNA only conferred ∼130-fold protection against adenine-N1 methylation, and this protection was reduced to ∼40-fold for A(syn)-T Hoogsteen base pairs embedded in a DNA-drug complex. Our results indicate that Watson-Crick faces of nucleobases are accessible to alkylating agents in canonical dsDNA and that Hoogsteen base pairs increase this accessibility. Given the higher abundance of dsDNA relative to ssDNA, these results suggest that dsDNA could be a substantial source of cytotoxic damage. The work establishes DMS probing as a method for characterizing A(syn)-T Hoogsteen base pairs in vitro and also lays the foundation for a sequencing approach to map A(syn)-T Hoogsteen and unpaired adenines genome-wide in vivo.

Manganese Ions Individually Alter the Reverse Transcription Signature of Modified Ribonucleosides.

Reverse transcription of RNA templates containing modified ribonucleosides transfers modification-related information as misincorporations, arrest or nucleotide skipping events to the newly synthesized cDNA strand. The frequency and proportion of these events, merged from all sequenced cDNAs, yield a so-called RT signature, characteristic for the respective RNA modification and reverse transcriptase (RT). While known for DNA polymerases in so-called error-prone PCR, testing of four different RTs by replacing Mg2+ with Mn2+ in reaction buffer revealed the immense influence of manganese chloride on derived RT signatures, with arrest rates on m1A positions dropping from 82% down to 24%. Additionally, we observed a vast increase in nucleotide skipping events, with single positions rising from 4% to 49%, thus implying an enhanced read-through capability as an effect of Mn2+ on the reverse transcriptase, by promoting nucleotide skipping over synthesis abortion. While modifications such as m1A, m22G, m1G and m3C showed a clear influence of manganese ions on their RT signature, this effect was individual to the polymerase used. In summary, the results imply a supporting effect of Mn2+ on reverse transcription, thus overcoming blockades in the Watson-Crick face of modified ribonucleosides and improving both read-through rate and signal intensity in RT signature analysis.

Machine learning of reverse transcription signatures of variegated polymerases allows mapping and discrimination of methylated purines in limited transcriptomes.

Reverse transcription (RT) of RNA templates containing RNA modifications leads to synthesis of cDNA containing information on the modification in the form of misincorporation, arrest, or nucleotide skipping events. A compilation of such events from multiple cDNAs represents an RT-signature that is typical for a given modification, but, as we show here, depends also on the reverse transcriptase enzyme. A comparison of 13 different enzymes revealed a range of RT-signatures, with individual enzymes exhibiting average arrest rates between 20 and 75%, as well as average misincorporation rates between 30 and 75% in the read-through cDNA. Using RT-signatures from individual enzymes to train a random forest model as a machine learning regimen for prediction of modifications, we found strongly variegated success rates for the prediction of methylated purines, as exemplified with N1-methyladenosine (m1A). Among the 13 enzymes, a correlation was found between read length, misincorporation, and prediction success. Inversely, low average read length was correlated to high arrest rate and lower prediction success. The three most successful polymerases were then applied to the characterization of RT-signatures of other methylated purines. Guanosines featuring methyl groups on the Watson-Crick face were identified with high confidence, but discrimination between m1G and m22G was only partially successful. In summary, the results suggest that, given sufficient coverage and a set of specifically optimized reaction conditions for reverse transcription, all RNA modifications that impede Watson-Crick bonds can be distinguished by their RT-signature.

Computational Insight into the Differential Mutagenic Patterns of O-Methylthymine Lesions.

Differential mutagenic patterns were recently reported for O-methylated thymine lesions, which indicate that O4-methylthymine (O4-Me-T) frequently leads to G misinsertions, whereas O2-methylthymine (O2-Me-T) is primarily nonmutagenic. The reasons for these differences are unclear since both lesions similarly alter the Watson-Crick binding face of T. To rationalize these replication outcomes at a molecular level, this work uses density functional theory calculations and molecular dynamics simulations to probe the lesion base-pairing properties as well as lesion accommodation by human polymerase η (pol η) and post-extension DNA duplexes. O4-Me-T forms two strong hydrogen bonds with an opposing G in the active site of pol η, which rationalizes the observed lesion mutagenicity. Nevertheless, dATP insertion opposite O4-Me-T can proceed through water-mediated hydrogen bonding, which is similar to the pathway previously proposed for pol η bypass of abasic sites and other T alkylation lesions. In contrast, the position of O2-Me-T in the pol η active site is dynamic due to the presence of the aberrant methyl group on the minor groove side of DNA. In fact, the experimental replication outcomes can only be rationalized when the syn glycosidic orientation of O2-Me-T is considered, which stabilizes the pre-insertion complex by placing the damage in the polymerase open pocket on the major groove side of DNA. Although dATP insertion can occur opposite syn-O2-Me-T through a water-mediated pathway similar to O4-Me-T replication, rotation about the glycosidic bond precludes a stable pol η ternary complex corresponding to dGTP insertion, which correlates with the reported nonmutagenic bypass of O2-Me-T. In addition to providing structural insights into the differential mutagenicity of methylated T adducts, our data highlight an emerging theme in the literature for the replication of pyrimidine alkylation products in noncanonical glycosidic orientations and sets the stage for future work on the replication of other alkylated lesions by TLS polymerases.

DFT and MD Studies of Formaldehyde-Derived DNA Adducts: Molecular-Level Insights into the Differential Mispairing Potentials of the Adenine, Cytosine, and Guanine Lesions.

4-(Methylnitrosamino)-1-(3-pyridyl)-1-butanone is a potent nicotine-based carcinogen that generates many DNA lesions, including the HOCH2-C, HOCH2-G, and HOCH2-A hydroxymethyl adducts. Despite all lesions containing an altered exocyclic amino group, which allows the hydroxymethyl group to be directed away from the Watson-Crick binding face, only the most persistent adenine adduct is mutagenic. As a first step toward understanding this differential mutagenicity, density functional theory (DFT) and molecular dynamics (MD) simulations were used to gain atomic-level structural details of these DNA damage products. DFT calculations reveal that all three lesions exhibit conformational diversity. However, regardless of the hydroxymethyl-nucleobase orientation, both DFT and MD simulations highlight that HOCH2-C and HOCH2-G form pairs with the canonical complementary base (G and C, respectively) that are structural and energetically preferred over mispairs. In contrast, depending on the hydroxymethyl-nucleobase orientation, the Watson-Crick HOCH2-A:T pair can become significantly destabilized relative to undamaged A:T. As a result, HOCH2-A mispairs with G, C, and A are energetically accessible and maintain key geometrical features of canonical DNA. Overall, our data directly correlate with the reported differential mutagenicity of the hydroxylmethyl lesions and will encourage future studies to further uncover the cellular impact of the most persistent adenine lesion.

Nucleobase pairing and photodimerization in a biologically derived metal-organic framework nanoreactor

Biologically derived metal-organic frameworks (bio-MOFs) are of great importance as they can be used as models for bio-mimicking and in catalysis, allowing us to gain insights into how large biological molecules function. Through rational design, here we report the synthesis of a novel bio-MOF featuring unobstructed Watson-Crick faces of adenine (Ade) pointing towards the MOF cavities. We show, through a combined experimental and computational approach, that thymine (Thy) molecules diffuse through the pores of the MOF and become base-paired with Ade. The Ade-Thy pair binding at 40–45% loading reveals that Thy molecules are packed within the channels in a way that fulfill both the Woodward-Hoffmann and Schmidt rules, and upon UV irradiation, Thy molecules dimerize into Thy<>Thy. This study highlights the utility of accessible functional groups within the pores of MOFs, and their ability to ‘lock’ molecules in specific positions that can be subsequently dimerized upon light irradiation, extending the use of MOFs as nanoreactors for the synthesis of molecules that are otherwise challenging to isolate.