Unveiling the solution structure of a DNA duplex with continuous silver-modified Watson-Crick base pairs.

The challenge of transforming organized DNA structures into their metallized counterparts persists in the scientific field. In this context, utilizing DNA molecules modified with 7-deazapurine, provides a transformative solution. In this study, we present the solution structure of a DNA duplex that can be transformed into its metallized equivalent while retaining the natural base pairing arrangement through the creation of silver-modified Watson-Crick base pairs. Unlike previously documented X-ray structures, our research demonstrates the feasibility of preserving the intrinsic DNA self-assembly while incorporating AgI into the double helix, illustrating that the binding of silver does not disrupt the canonical base-pairing organization. Moreover, in our case, the uninterrupted AgI chain deviates from forming conventional straight linear chains; instead, it adheres to a helical arrangement dictated by the underlying DNA structure. This research challenges conventional assumptions and opens the door to precisely design structures based on the organization of highly stable Ag-DNA assemblies.

Investigation of the Nonradiative Photoprocesses of Unnatural DNA Base: 7-(2-Thienyl)-imidazo[4,5-b]pyridine (Ds)─A Computational Study.

7-(2-Thienyl)-imidazo[4,5-b]pyridine (Ds) is an unnatural nucleic acid that forms a stable pair with pyrrole-2-carbaldehyde (Pa) in DNA. This Ds-Pa pair gets stabilized via van der Waals interaction and shape fitting. In our previous study [Ghosh, P. J. Phys. Chem. A 2021, 125, 5556-5561], we investigated the nonradiative photoprocesses of the unnatural DNA base Pa, and also there are some studies on its stability and reactivity in the ground state. But, to consider it as a good unnatural base pair, one has to understand its stability not only in the ground state but also in the excited states after absorbing ultraviolet (UV) radiation. Therefore, in this study, the excited-state photoprocesses of Ds on UV irradiation and its nonradiative decay channels have been investigated using state-of-the-art multireference methods, and this investigation finally leads the molecule to access the minimum energy crossing point (MECP) via a downhill pathway.

Beyond the "spine of hydration": Chiral SFG spectroscopy detects DNA first hydration shell and base pair structures.

Experimental methods capable of selectively probing water at the DNA minor groove, major groove, and phosphate backbone are crucial for understanding how hydration influences DNA structure and function. Chiral-selective sum frequency generation spectroscopy (chiral SFG) is unique among vibrational spectroscopies because it can selectively probe water molecules that form chiral hydration structures around biomolecules. However, interpreting chiral SFG spectra is challenging since both water and the biomolecule can produce chiral SFG signals. Here, we combine experiment and computation to establish a theoretical framework for the rigorous interpretation of chiral SFG spectra of DNA. We demonstrate that chiral SFG detects the N-H stretch of DNA base pairs and the O-H stretch of water, exclusively probing water molecules in the DNA first hydration shell. Our analysis reveals that DNA transfers chirality to water molecules only within the first hydration shell, so they can be probed by chiral SFG spectroscopy. Beyond the first hydration shell, the electric field-induced water structure is symmetric and, therefore, precludes chiral SFG response. Furthermore, we find that chiral SFG can differentiate chiral subpopulations of first hydration shell water molecules at the minor groove, major groove, and phosphate backbone. Our findings challenge the scientific perspective dominant for more than 40 years that the minor groove "spine of hydration" is the only chiral water structure surrounding the DNA double helix. By identifying the molecular origins of the DNA chiral SFG spectrum, we lay a robust experimental and theoretical foundation for applying chiral SFG to explore the chemical and biological physics of DNA hydration.

Beyond the base pairs: comparative genome-wide DNA methylation profiling across sequencing technologies.

Deoxyribonucleic acid (DNA) methylation plays a key role in gene regulation and is critical for development and human disease. Techniques such as whole-genome bisulfite sequencing (WGBS) and reduced representation bisulfite sequencing (RRBS) allow DNA methylation analysis at the genome scale, with Illumina NovaSeq 6000 and MGI Tech DNBSEQ-T7 being popular due to their efficiency and affordability. However, detailed comparative studies of their performance are not available. In this study, we constructed 60 WGBS and RRBS libraries for two platforms using different types of clinical samples and generated approximately 2.8 terabases of sequencing data. We systematically compared quality control metrics, genomic coverage, CpG methylation levels, intra- and interplatform correlations, and performance in detecting differentially methylated positions. Our results revealed that the DNBSEQ platform exhibited better raw read quality, although base quality recalibration indicated potential overestimation of base quality. The DNBSEQ platform also showed lower sequencing depth and less coverage uniformity in GC-rich regions than did the NovaSeq platform and tended to enrich methylated regions. Overall, both platforms demonstrated robust intra- and interplatform reproducibility for RRBS and WGBS, with NovaSeq performing better for WGBS, highlighting the importance of considering these factors when selecting a platform for bisulfite sequencing.

Investigating the Effect of Chemical Modifications on the Ribose Sugar Conformation, Watson-Crick Base Pairing, and Intrastrand Stacking Interactions: A Theoretical Approach.

Over the last few decades, chemically modified sugars have been incorporated into nucleic acid-based therapeutics to improve their pharmacological potential. Chemical modification can influence the sugar conformation, Watson-Crick hydrogen (W-C) bonding, and nucleobase stacking interactions, which play major roles in the structural integrity and dynamic properties of nucleic acid duplexes. In this study, we categorized 33 uridine (U*) and cytidine (C*) sugar modifications and calculated their sugar conformational parameters. We also calculated the Watson-Crick hydrogen bond energies of the modified RNA-type base pairs (U*:A and C*:G) using DFT and sSAPT0 methods. The W-C base pairing energy calculations suggested that the South-type modified sugar strengthens the C*:G base pair and weakens the U*:A base pair compared to the unmodified one. In contrast, the North-type sugar modifications form weaker C*:G base pair and marginally stronger U*:A base pair compared to the South-type modified sugars. Moreover, intrastrand base stacking energies were calculated for 15 modifications incorporated at the fourth position in 7-mer non-self-complementary RNA duplexes [(GCAU*GAC)2 and (GCAC*GAC)2], utilizing molecular dynamics simulation and quantum mechanical (DFT and sSAPT0) methods. The sugar modifications were found to have minimal effect on the intrastrand base-stacking interactions. However, the glycol nucleic acid modification disturbs the intrastrand base-stacking significantly, corroborating the experimental data.

Insights into bacterial metabolism from small RNAs.

The study of small, regulatory RNAs (sRNA) that act by base-pairing with target RNAs in bacteria has been steadily advancing, particularly with the availability of more and more transcriptome and RNA-RNA interactome datasets. While the characterization of multiple sRNAs has helped to elucidate their mechanisms of action, these studies also are providing insights into protein function, control of metabolic flux, and connections between metabolic pathways as we will discuss here. In describing several examples of the metabolic insights gained, we will summarize the different types of base-pairing sRNAs including mRNA-derived sRNAs, sponge RNAs, RNA mimics, and dual-function RNAs as well as suggest how information about sRNAs could be exploited in the future.

7-Deaza-2'-deoxyisoguanosine, a Noncanonical Nucleoside for Nucleic Acid Code Expansion and New DNA Constructs: Nucleobase Functionalization of Inverse Watson-Crick and Purine-Purine Base Pairs.

7-Deaza-2'-deoxyisoguanosine forms stable inverse Watson-Crick base pairs with 5-methyl-2'-deoxyisocytidine and purine-purine base pairs with 2'-deoxyguanosine or 5-aza-7-deaza-2'-deoxyguanosine. Both base pairs expand the genetic coding system. The manuscript reports on the functionalization of these base pairs with halogen atoms and clickable side chains introduced at 7-position of the 7-deazapurine base. Oligonucleotides containing the functionalized base pairs were prepared by solid-phase synthesis. To this end, a series of phosphoramidites were synthesized and clickable side chains with short and long linkers were incorporated in oligonucleotides. Fluorescent pyrene conjugates were obtained by postmodification. Functionalization of DNA with a single inverse Watson-Crick base pair by halogens or clickable residues has only a minor impact on duplex stability. Pyrene click adducts increase (long linker) or decrease (short linker) the double helix stability. Stable hybrid duplexes were constructed containing three consecutive purine-purine pairs of 7-functionalized 7-deaza-2'-deoxyisoguanine with guanine or 5-aza-7-deazaguanine in the center and Watson-Crick pairs at both ends. The incorporation of a hybrid base pair tract of 7-deaza-2'-deoxyisoguanosine/5-aza-7-deaza-2'-deoxyguanosine pairs stabilizes the double helix strongly. Fluorescence intensity of pyrene short linker adducts increased when the 7-deazapurine base was positioned opposite to 5-methylisocytosine (inverse base pair) compared to purine-purine base pairs with guanine or 5-aza-7-deazaguanine in opposite positions. For long liker adducts, the situation is more complex. Circular dichroism (CD) spectra of purine DNA differ to those of Watson-Crick double helices and are indicative for the new DNA constructs. The impact of 7-deaza-2'-deoxyisoguanine base pair functionalization is studied for the first time and all experimental details are reported to prepare DNA functionalized at the 7-deazaisoguanine site. The influence of single and multiple incorporations on DNA structure and stability is shown. Clickable residues introduced at the 7-position of the 7-deazaisoguanine base provide handles for Huisgen-Sharpless-Meldal click cycloadditions without harming the stability of purine-pyrimidine and purine-purine base pairs. Other chemistries might be used for bioconjugation. Our investigation paves the way for the functionalization of a new DNA related recognition system expanding the common Watson-Crick regime.

The effect of the loop on the thermodynamic and kinetic of single base pair in pseudoknot.

RNA pseudoknots are RNA molecules with specialized three-dimensional structures that play important roles in various biological processes. To understand the functions and mechanisms of pseudoknots, it is essential to elucidate their structures and folding pathways. The most fundamental step in RNA folding is the opening and closing of a base pair. The effect of flexible loops on the base pair in pseudoknots remains unclear. In this work, we use molecular dynamics simulations and Markov state model to study the configurations, thermodynamic and kinetic of single base pair in pseudoknots. We find that the presence of the loop leads to a trap state. In addition, the rate-limiting step for the formation of base pair is the disruption of the trap state, rather than the open state to the closed state, which is quite different from the previous studies on non-pseudoknot RNA. For the thermodynamic parameters in pseudoknots, we find that the entropy difference upon opening the base pair between this simulation and the nearest-neighbor model results from the different entropy of different lengths of loop in solution. The thermodynamic parameters of the stack in pseudoknot are close to the nearest-neighbor parameters. The bases on the loop have different distribution patterns in different states, and the slow transition states of the loop are determined by the orientation of the bases.

Universal cold RNA phase transitions.

RNA's diversity of structures and functions impacts all life forms since primordia. We use calorimetric force spectroscopy to investigate RNA folding landscapes in previously unexplored low-temperature conditions. We find that Watson-Crick RNA hairpins, the most basic secondary structure elements, undergo a glass-like transition below [Formula: see text]C where the heat capacity abruptly changes and the RNA folds into a diversity of misfolded structures. We hypothesize that an altered RNA biochemistry, determined by sequence-independent ribose-water interactions, outweighs sequence-dependent base pairing. The ubiquitous ribose-water interactions lead to universal RNA phase transitions below TG, such as maximum stability at [Formula: see text]C where water density is maximum, and cold denaturation at [Formula: see text]C. RNA cold biochemistry may have a profound impact on RNA function and evolution.

The Influence of 2'-Deoxyguanosine Lesions on the Electronic Properties of OXOG:::C Base Pairs in Ds-DNA: A Comparative Analysis of Theoretical Studies.

DNA is continuously exposed to a variety of harmful factors, which, on the one hand, can force undesirable processes such as ageing, carcinogenesis and mutagenesis, while on the other hand, can accelerate evolutionary changes. Of all the canonical nucleosides, 2'-deoxyguanosine (dG) exhibits the lowest ionization potential, making it particularly prone to the one-electron oxidizing process. The most abundant type of nucleobase damage is constituted by 7,8-dihydro-8-oxo-2'-deoxyguanosine (OXOdG), with an oxidation potential that is 0.56 V lower than that of canonical dG. All this has led to OXOdG, as an isolated lesion, being perceived as a sink for radical cations in the genome. In this paper, a comparative analysis of the electronic properties of an OXOGC base pair within the context of a clustered DNA lesion (CDL) has been conducted. It is based on previous DFT studies that were carried out at the M06-2x/6-31++G** level of theory in non-equilibrated and equilibrated condensed phases. The results of the comparative analysis presented here reveal the following: (A) The ionization potentials of OXOG4C2 were largely unaffected by a second lesion. (B) The positive charge and spin were found predominantly on the OXOG4C2 moiety. (C) The electron-hole transfers A3T3→G4C2 and G4C2←A5T1 were found in the Marcus inverted region and were resistant to the presence of a second DNA lesion in close proximity. It can therefore be reasonably postulated that OXOGC becomes the sink for a radical cation migrating through the double helix, irrespective of the presence of other 2'-deoxyguanosine lesions in the CDL structure.

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.

DNA Hierarchical Superstructures from Micellar Units: Stiff Hydrogels and Anisotropic Nanofibers.

Supramolecular materials have been assembled using a wide range of interactions, including the hydrophobic effect, DNA base-pairing, and hydrogen bonding. Specifically, DNA amphiphiles with a hydrophobic building block self-assemble into diverse morphologies depending on the length and composition of both blocks. Herein, we take advantage of the orthogonality of different supramolecular interactions - the hydrophobic effect, Watson-Crick-Franklin base pairing and RNA kissing loops - to create hierarchical self-assemblies with controlled morphologies on both the nanometer and the micrometer scales. Assembly through base-pairing leads to the formation of hybrid, multi-phasic hydrogels with high stiffness and self-healing properties. Assembly via hydrophobic core interactions gives anisotropic, discrete assemblies, where DNA fibers with one sequence are terminated with DNA spheres bearing different sequences. This work opens new avenues for the bottom-up construction of DNA-based materials, with promising applications in drug delivery, tissue engineering, and the creation of complex DNA structures from a minimum array of components.

Excited-Ground-State Transition of the RNA Strand Slippage Mechanism Captured by the Base-Specific Force Field.

Excited-ground-state transition and strand slippage of RNA play key roles in transcription and translation of central dogma. Due to limitation of current experimental techniques, the dynamic structure ensembles of RNA remain inadequately understood. Molecular dynamics simulations offer a promising complementary approach, whose accuracy depends on the force field. Here, we develop the new version of RNA base-specific force field (BSFF2) to address underestimation of base pairing stability and artificial backbone conformations. Extensive evaluations on typical RNA systems have comprehensively confirmed the accuracy of BSFF2. Furthermore, BSFF2 demonstrates exceptional efficiency in de novo folding of tetraloops and reproducing base pair reshuffling transition between RNA excited and ground states. Then, we explored the RNA strand slippage mechanism with BSFF2. We conducted a comprehensive three-dimensional structural investigation into the strand slippage of the most complex r(G4C2)9 repeat element and presented the molecular details in the dynamic transition along with the underlying mechanism. Our results of capturing the strand slippage, excited-ground transition, de novo folding, and simulations for various typical RNA motifs indicate that BSFF2 should be one of valuable tools for dynamic conformation research and structure prediction of RNA, and a future contribution to RNA-targeted drug design as well as RNA therapy development.

GeF-seq: A Simple Procedure for Base-Pair Resolution ChIP-seq.

Nucleotide sequences recognized and bound by DNA-binding proteins (DBPs) are critical to controlling and maintaining gene expression, replication, chromosome segregation, cell division, and nucleoid structure in bacterial cells. Therefore, determination of the binding sequences of DBPs is important not only to study DBP recognition mechanisms but also to understand the fundamentals of cell homeostasis. While ChIP-seq analysis appears to be an effective way to determine DBP binding sites on the genome, the resolution is sometimes not sufficient to identify the sites precisely. Here we introduce a simple and effective method named Genome Footprinting with high-throughput sequencing (GeF-seq) to determine binding sites of DBPs with single base-pair resolution. GeF-seq detects binding sites of DBPs as sharp peaks and thus makes it possible to identify the recognition sequence in each "binding peak" more easily and accurately compared to the common ChIP-seq.

Unexpected large charge transfer rate mediated by adenine in twisted DNA structure.

DNA exhibits remarkable charge transfer ability, which is crucial for its biological functions and potential electronic applications. The charge transfer process in DNA is widely recognized as primarily mediated by guanine, while the contribution of other nucleobases is negligible. Using the tight-binding models in conjunction with first-principles calculations, we investigated the charge transfer behavior of homogeneous GC and AT pairs. We found that the charge transfer rate of adenine significantly changes. With overstretching, the charge transfer rate of adenine can even surpass that of guanine, by as much as five orders of magnitude at a twist angle of around 26°. Further analysis reveals that it is attributed to the turnover of the relative coupling strength between homogeneous GC and AT base pairs, which is caused by the symmetry exchange between the two highest occupied molecular orbitals of base pairs occurring at different twist angles. Given the high degree of flexibility of DNA in vivo and in vitro conditions, these findings prompt us to reconsider the mechanism of biological functions concerning the charge transfer in DNA molecules and further open the potential of DNA as a biomaterial for electronic applications.

Phase Separation Modulates the Thermodynamics and Kinetics of RNA Hybridization.

Biomolecular condensates can influence cellular function in a number of ways, including by changing the structural dynamics and conformational equilibria of the molecules partitioned within them. Here we use methyl transverse relaxation optimized spectroscopy (methyl-TROSY) NMR in conjunction with 2'-O-methyl labeling of RNA to characterize the thermodynamics and kinetics of RNA-RNA base pairing in condensates formed by the C-terminal intrinsically disordered region of CAPRIN1, an RNA-binding protein involved in RNA transport, translation, and stability. CAPRIN1 condensates destabilize RNA-RNA base pairing, resulting from a ∼270-fold decrease and a concomitant ∼15-fold increase in the on- and off-rates for duplex formation, respectively. The ∼30-fold slower diffusion of RNA single strands within the condensed phase partially accounts for the reduced on-rate, but the further ∼9-fold reduction likely reflects shedding of CAPRIN1 chains that are interacting with the RNA prior to hybridization. Our study emphasizes the important role of protein solvation in modulating nucleic acid recognition processes inside condensates.

Predicting the effect of binding molecules on the shape and mechanical properties of structured DNA assemblies.

Chemo-mechanical deformation of structured DNA assemblies driven by DNA-binding ligands has offered promising avenues for biological and therapeutic applications. However, it remains elusive how to effectively model and predict their effects on the deformation and mechanical properties of DNA structures. Here, we present a computational framework for simulating chemo-mechanical change of structured DNA assemblies. We particularly quantify the effects of ethidium bromide (EtBr) intercalation on the geometry and mechanical properties of DNA base-pairs through molecular dynamics simulations and integrated them into finite-element-based structural analysis to predict the shape and properties of DNA objects. The proposed model captures various structural changes induced by EtBr-binding such as shape variation, flexibility modulation, and supercoiling instability. It enables a rational design of structured DNA assemblies with tunable shapes and mechanical properties by binding molecules.

7-Deazapurine and Pyrimidine Nucleoside and Oligonucleotide Cycloadducts Formed by Inverse Diels-Alder Reactions with 3,6-Di(pyrid-2-yl)-1,2,4,5-tetrazine: Ethynylated and Vinylated Nucleobases for Functionalization and Impact of Pyridazine Adducts on DNA Base Pair Stability and Mismatch Discrimination.

The manuscript reports on 7-deazapurine and pyrimidine nucleoside and oligonucleotide cycloadducts formed by the inverse electron demand Diels-Alder (iEDDA) reaction with 3,6-di(pyrid-2-yl)-1,2,4,5-tetrazine. Cycloadducts were constructed from ethynylated and vinylated nucleobases. Oligonucleotides were synthesized containing iEDDA modifications, and the impact on duplex stability was investigated. iEDDA reactions were performed on nucleoside triple bond side chains. Oxidation was not required in these cases as dihydropyridazine intermediates are not formed. In contrast, oxidation is necessary for reactions performed on alkenyl compounds. This was verified on 5-vinyl-2'-deoxyuridine. A diastereomeric mixture of 1,2-dihydropyridazine cycloadduct intermediates was isolated, characterized, and later oxidized. 12-mer oligonucleotides containing 1,2-pyridazine inverse Diels-Alder cycloadducts and their precursors were hybridized to short DNA duplexes. For that, a series of phosphoramidites was prepared. DNA duplexes with 7-functionalized 7-deazaadenines and 5-functionalized pyrimidines display high duplex stability when spacer units are present between nucleobases and pyridazine cycloadducts. A direct connectivity of the pyridazine moiety to nucleobases as reported for metabolic labeling of vinyl nucleosides reduced duplex stability strongly. Oligonucleotides bearing linkers with and without pyridazine cycloadducts attached to the 7-deazaadenine nucleobase significantly reduced mismatch formation with dC and dG.

Comprehensive Assessment of Force-Field Performance in Molecular Dynamics Simulations of DNA/RNA Hybrid Duplexes.

Mixed double helices formed by RNA and DNA strands, commonly referred to as hybrid duplexes or hybrids, are essential in biological processes like transcription and reverse transcription. They are also important for their applications in CRISPR gene editing and nanotechnology. Yet, despite their significance, the hybrid duplexes have been seldom modeled by atomistic molecular dynamics methodology, and there is no benchmark study systematically assessing the force-field performance. Here, we present an extensive benchmark study of polypurine tract (PPT) and Dickerson-Drew dodecamer hybrid duplexes using contemporary and commonly utilized pairwise additive and polarizable nucleic acid force fields. Our findings indicate that none of the available force-field choices accurately reproduces all the characteristic structural details of the hybrid duplexes. The AMBER force fields are unable to populate the C3'-endo (north) pucker of the DNA strand and underestimate inclination. The CHARMM force field accurately describes the C3'-endo pucker and inclination but shows base pair instability. The polarizable force fields struggle with accurately reproducing the helical parameters. Some force-field combinations even demonstrate a discernible conflict between the RNA and DNA parameters. In this work, we offer a candid assessment of the force-field performance for mixed DNA/RNA duplexes. We provide guidance on selecting utilizable force-field combinations and also highlight potential pitfalls and best practices for obtaining optimal performance.