Urea Mimics Nucleobases by Preserving the Helical Integrity of B-DNA Duplexes via Hydrogen Bonding and Stacking Interactions

Systematic Selection of High-Affinity ssDNA Sequences to Carbon Nanotubes

Single-walled carbon nanotubes (SWCNTs) have gained significant interest for their potential in biomedicine and nanoelectronics. The functionalization of SWCNTs with single-stranded DNA (ssDNA) enables the precise control of SWCNT alignment and the development of optical and electronic biosensors. This study addresses the current gaps in the field by employing high-throughput systematic selection, enriching high-affinity ssDNA sequences from a vast random library. Specific base compositions and patterns are identified that govern the binding affinity between ssDNA and SWCNTs. Molecular dynamics simulations validate the stability of ssDNA conformations on SWCNTs and reveal the pivotal role of hydrogen bonds in this interaction. Additionally, it is demonstrated that machine learning could accurately distinguish high-affinity ssDNA sequences, providing an accessible model on a dedicated webpage (http://service.k-medai.com/ssdna4cnt). These findings open new avenues for high-affinity ssDNA-SWCNT constructs for stable and sensitive molecular detection across diverse scientific disciplines.

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.

Enhancement of the thermal stability of G-quadruplex structures by urea

The solute urea has been used extensively as a denaturant in protein folding studies; double-stranded nucleic acid structures are also destabilized by urea, but comparatively less than proteins. In previous research, the solute has been shown to strongly destabilize folded G-quadruplex DNA structures. This contribution demonstrates the stabilizing effect of urea on the G-quadruplex formed by the oligodeoxyribonucleotide (ODN), G3T (d[5'-GGGTGGGTGGGTGGG-3']), and related sequences in the presence of sodium or potassium cations. Stabilization is observed up to 7 M urea, which was the highest concentration we investigated. The folded structure of G3T has three G-tetrads and three loops that consist of single thymine residues. ODNs related to G3T, in which the thymine residues in the loop are substituted by adenosine residues, also exhibit enhanced stability in the presence of molar concentrations of urea. The circular dichroism (CD) spectra of these ODNs in the presence of urea are consistent with that of a G-quadruplex. As the urea concentration increases, the spectral intensities of the peaks and troughs change, while their positions change very little. The heat-induced transition from the folded to unfolded state, T_m, was measured by monitoring the change in the UV absorption as a function of temperature. G-quadruplex structures with loops containing single bases exhibited large increases in T_m with increasing urea concentrations. These data imply that the loop region play a significant role in the thermal stability of tetra-helical DNA structures in the presence of the solute urea.

DNA Materials Assembled from One DNA Strand

Due to the specific base-pairing recognition, clear nanostructure, programmable sequence and responsiveness of the DNA molecule, DNA materials have attracted extensive attention and been widely used in controlled release, drug delivery and tissue engineering. Generally, the strategies for preparing DNA materials are based on the assembly of multiple DNA strands. The construction of DNA materials using only one DNA strand can not only save time and cost, but also avoid defects in final assemblies generated by the inaccuracy of DNA ratios, which potentially promote the large-scale production and practical application of DNA materials. In order to use one DNA strand to form assemblies, the sequences have to be palindromes with lengths that need to be controlled carefully. In this review, we introduced the development of DNA assembly and mainly summarized current reported materials formed by one DNA strand. We also discussed the principle for the construction of DNA materials using one DNA strand.

Nucleic acid extraction from complex biofluid using toothpick-actuated over-the-counter medical-grade cotton

Nucleic acid amplification technique (NAAT)-assisted detection is the primary intervention for pathogen molecular diagnostics. However, NAATs such as quantitative real-time polymerase chain reaction (qPCR) require prior purification or extraction of target nucleic acid from the sample of interest since the latter often contains polymerase inhibitors. Similarly, genetic disease screening is also reliant on the successful extraction of pure patient genomic DNA from the clinical sample. However, such extraction techniques traditionally utilize spin-column techniques that in turn require centralized high-speed centrifuges. This hinders any potential deployment of qPCR- or PCR-like NAAT methods in resource-constrained settings. The development of instrument-free nucleic acid extraction methods, especially those utilizing readily available materials would be of great interest and benefit to NAAT-mediated molecular diagnosis workflows in resource-constrained settings. In this report, we screened medical-grade cotton, a readily available over-the-counter biomaterial to extract genomic DNA (gDNA) spiked in 30 %, 45 %, and 60 % serum or cell lysate. The extraction was carried out in a completely instrument-free manner using cotton and a sterilized toothpick and was completed in 30 min (with using chaotropic salt) or 10 min (without using chaotropic salt). The quality of the extracted DNA was then probed using PCR followed by agarose gel analysis for preliminary validation of the study. The qPCR experiments then quantitatively established the extraction efficiency (0.3-27 %, depending on serum composition). Besides, percent similarity score obtained from the Sanger sequencing experiments probed the feasibility of extracted DNA towards polymerase amplification with fluorescent nucleotide incorporation. Overall, our method demonstrated that DNA extraction could be performed utilizing toothpick-mounted cotton both with or without using a chaotropic salt, albeit with a difference in the quality of the extracted DNA.

Analysis of nucleotide insertion opposite urea and translesion synthesis across urea by DNA polymerases

Abstract

Urea (Ua) is produced in DNA as the result of oxidative damage to thymine and guanine. It was previously reported that Klenow fragment (Kf) exo⁻ incorporated dATP opposite Ua, and that DNA polymerase β was blocked by Ua. We report here the following nucleotide incorporations opposite Ua by various DNA polymerases: DNA polymerase α, dATP and dGTP (dATP > dGTP); DNA polymerase δ, dATP; DNA polymerase ζ, dATP; Kf exo⁻, dATP; Sulfolobus solfataricus P2 DNA polymerase IV (Dpo4), dGTP and dATP (dGTP > dATP); and DNA polymerase η, dCTP, dGTP, dATP, and dTTP (dCTP > dGTP > dATP > dTTP). DNA polymerases β and ε were blocked by Ua. Elongation by DNA polymerases δ and ζ stopped after inserting dATP opposite Ua. Importantly, the elongation efficiency to full-length beyond Ua using DNA polymerase η and Dpo4 were almost the same as that of natural DNA.

Graphical abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s41021-022-00236-3.

Sequence specific hydrogen bond of DNA with denaturants affects its stability: Spectroscopic and simulation studies

BACKGROUND

Several different small molecules have been used to target the DNA helix in order to treat the diseases caused by its mutation. Guanidinium(Gdm⁺) and urea based drugs have been used for the diseases related to central nervous system, also as the anti-inflammatory and chemotherapeutic agent. However, the role of Gdm⁺ and urea in the stabilization/destabilization of DNA is not well understood.

METHODS

Spectroscopic techniques along with molecular dynamics (MD) simulation have been performed on different sequences of DNA in the presence of guanidinium chloride (GdmCl) and urea to decode the binding of denaturants with DNA and the role of hydrogen bond with the different regions of DNA in its stability/destability.

RESULTS AND CONCLUSION

Our study reveals that, Gdm⁺ of GdmCl and urea both intrudes into the groove region of DNA along with the interaction with its phosphate backbone. However, interaction of Gdm⁺ and urea with the nucleobases in the groove region is different. Gdm⁺ forms the intra-strand hydrogen bond with the central region of the both sequences of DNA whereas inter-strand hydrogen bond along with water assisted hydrogen bond takes place in the case of urea. The intra-strand hydrogen bond formation capability of Gdm⁺ with the nucleobases in the minor groove of DNA decreases its groove width which probably causes the stabilization of B-DNA in GdmCl. In contrast, the propensity of the formation of inter-strand hydrogen bond of urea with the nucleobases in the groove region of DNA without affecting the groove width destabilizes B-DNA as compared to GdmCl. This study depicts that the opposite effect of GdmCl and urea on the stability is a general property of B-DNA. However, the extent of stabilization/destabilization of DNA in Gdm⁺ and urea depend on its sequence probably due to the difference in the intra/inter-strand hydrogen bonding with different bases present in both the sequences of DNA.

GENERAL SIGNIFICANCE

The information obtained from this study will be useful for the designing of Gdm⁺ based drug molecule which can target the DNA more specifically and selectively.

Urea-aromatic interactions in biology

Noncovalent interactions are key determinants in both chemical and biological processes. Among such processes, the hydrophobic interactions play an eminent role in folding of proteins, nucleic acids, formation of membranes, protein-ligand recognition, etc.. Though this interaction is mediated through the aqueous solvent, the stability of the above biomolecules can be highly sensitive to any small external perturbations, such as temperature, pressure, pH, or even cosolvent additives, like, urea-a highly soluble small organic molecule utilized by various living organisms to regulate osmotic pressure. A plethora of detailed studies exist covering both experimental and theoretical regimes, to understand how urea modulates the stability of biological macromolecules. While experimentalists have been primarily focusing on the thermodynamic and kinetic aspects, theoretical modeling predominantly involves mechanistic information at the molecular level, calculating atomistic details applying the force field approach to the high level electronic details using the quantum mechanical methods. The review focuses mainly on examples with biological relevance, such as (1) urea-assisted protein unfolding, (2) urea-assisted RNA unfolding, (3) urea lesion interaction within damaged DNA, (4) urea conduction through membrane proteins, and (5) protein-ligand interactions those explicitly address the vitality of hydrophobic interactions involving exclusively the urea-aromatic moiety.

Energetic, Structural and Dynamic Properties of Nucleobase-Urea Interactions that Aid in Urea Assisted RNA Unfolding

Understanding the structure-function relationships of RNA has become increasingly important given the realization of its functional role in various cellular processes. Chemical denaturation of RNA by urea has been shown to be beneficial in investigating RNA stability and folding. Elucidation of the mechanism of unfolding of RNA by urea is important for understanding the folding pathways. In addition to studying denaturation of RNA in aqueous urea, it is important to understand the nature and strength of interactions of the building blocks of RNA. In this study, a systematic examination of the structural features and energetic factors involving interactions between nucleobases and urea is presented. Results from molecular dynamics (MD) simulations on each of the five DNA/RNA bases in water and eight different concentrations of aqueous urea, and free energy calculations using the thermodynamic integration method are presented. The interaction energies between all the nucleobases with the solvent environment and the transfer free energies become more favorable with respect to increase in the concentration of urea. Preferential interactions of urea versus water molecules with all model systems determined using Kirkwood-Buff integrals and two-domain models indicate preference of urea by nucleobases in comparison to water. The modes of interaction between urea and the nucleobases were analyzed in detail. In addition to the previously identified hydrogen bonding and stacking interactions between urea and nucleobases that stabilize the unfolded states of RNA in aqueous solution, NH-π interactions are proposed to be important. Dynamic properties of each of these three modes of interactions have been presented. The study provides fundamental insights into the nature of interaction of urea molecules with nucleobases and how it disrupts nucleic acids.

Towards accurate infrared spectral density of weak H-bonds in absence of relaxation mechanisms

Following the previous theoretical developments to completely reproduce the IR spectra of weak hydrogen bond complexes within the framework of the linear response theory (LRT), the quantum theory of the high stretching mode spectral density (SD) of weak H-bonds is reconsidered. Within the LRT theory, the SD is the one sided Fourier transform of the autocorrelation function (ACF) of the high stretching mode dipole moment operator. In order to provide more accurate theoretical bandshapes, we have explored the equivalence between the SDs given in previous studies with respect to a new quantum one, and revealed that in place of the basic equations used in the precedent works for which the SD I_Old(ω)=2Re∫₀^∞G_Old(t)e^-iωtdt where the ACF G_Old(t) = ⟨μ(0)μ(t)⁺⟩ = tr {ρ {μ(0)} {μ(t)}⁺}, one can use a new expression for the SD, given by I_New(ω)=2ωRe∫₀^∞G_New(t)e^-iωtdt where G_New(t)=μ(0)μ(t)⁺=1βtrρ_B∫₀^βμ(0)μ(t+iλℏ)⁺dλ. Here ρ_B is the Boltzmann density operator, μ(0) the dipole moment operator at initial time and μ(t) the dipole moment operator at time t in the Heisenberg picture, ℏ is the Planck constant, β is the inverse of the Boltzmann factor k_BT where T is the absolute temperature and k_B the Boltzmann constant. Using this formalism, we demonstrated that the new quantum approach gives the same final SD as used by previous models, and reduces to the Franck-Condon progression appearing in the Maréchal and Witkowski's pioneering approach when the relaxation mechanisms are ignored. Results of this approach shed light on the equivalence between the quantum and classical IR SD approaches for weak H-bonds in absence of medium surroundings effect, which has been a subject of debate for decades.

Equivalence between the Classical and Quantum IR Spectral Density Approaches of Weak H-Bonds in the Absence of Damping

The aim of this paper is to overhaul the quantum elucidation of the spectral density (SD) of weak H-bonds treated without taking into account any of the damping mechanisms. The reconsideration of the SD is performed within the framework the linear response theory. Working in the setting of the strong anharmonic coupling theory and the adiabatic approximation, the simplified expression of the classical SD, in the absence of dampings, is equated to be I_Cl(ω) = Re[∫₀^∞G_Cl(t)e^-iΩt dt] in which the classical-like autocorrelation function (ACF), G_Cl(t), is given by G_Cl(t) = tr{ρ(β){μ(0)}{μ(t)}^†}. With this consideration, we have shown that the classical SD is equivalent to the line shape obtained by F(ω) = ΩI_Cl(ω), which in turn is equivalent to the quantum SD given by I_Qu(ω) = Re[∫₀^∞G_Qu(t)e^-iΩt dt], where G_Qu(t) is the corresponding quantum ACF having for expression G_Qu(t) = (1/β) tr{ρ∫₀^β[μ(0)}{μ(t + iλℏ)}^† dλ}. Thus, we have shown that for weak H-bonds dealt without dampings, the SDs obtained by the quantum approaches are equivalent to the SDs geted by the classical approach in which the incepation ACF is, however, of quantum nature and where the line shape is the Fourier transform of the ACF times the angular frequency. It is further shown that the classical approach dealing with the SD of weak H-bonds leads identically to the result found by Maréchal and Witkowski in their pioneering quantum treatment where they ignored the linear response theory and dampings.