Search

The thermodynamic landscape of the CRISPR/Cas9 system plays a crucial role in understanding and optimizing the performance of this revolutionary genome-editing technology. In this research, we utilized isothermal titration calorimetry and microscale thermophoresis techniques to thoroughly investigate the thermodynamic properties governing CRISPR/Cas9 interactions. Our findings revealed that the binding between sgRNA and Cas9 is primarily governed by entropy, which compensates for an unfavorable enthalpy change. Conversely, the interaction between the CRISPR RNP complex and the target DNA is characterized by a favorable enthalpy change, offsetting an unfavorable entropy change. Notably, both interactions displayed negative heat capacity changes, indicative of potential hydration, ionization, or structural rearrangements. However, we noted that the involvement of water molecules and counterions in the interactions is minimal, suggesting that structural rearrangements play a significant role in influencing the binding thermodynamics. These results offer a nuanced understanding of the energetic contributions and structural dynamics underlying CRISPR-mediated gene editing. Such insights are invaluable for optimizing the efficiency and specificity of CRISPR-based genome editing applications, ultimately advancing our ability to precisely manipulate genetic material in various organisms for research, therapeutic, and biotechnological purposes.

A DNA nick, defined as a discontinuity in a double-stranded DNA molecule where the phosphodiester bond between adjacent nucleotides of one strand is absent due to enzyme action, serves as an effective mechanism to alleviate stress in supercoiled DNA. This stress release is essential for the smooth operation of transcriptional machinery. However, the underlying mechanisms and their impact on protein search dynamics, which are crucial for initiating transcription, remain unclear. Through extensive computer simulations, we unravel the molecular picture, demonstrating that intramolecular stress release due to a DNA nick is driven by a combination of writhing and twisting motions, depending on the nick's position. This stress release is quantitatively manifested as a step-like increase in the linking number. Furthermore, we elucidate that the nicked supercoiled minicircles exhibit enhanced torsional dynamics, promoting rapid conformational changes and frequent shifts in the identities of juxtaposed DNA sites on the plectoneme. The dynamics of the juxtaposition sites facilitates communication between protein and DNA, resulting in faster protein diffusion compared with native DNA with the same topology. Our findings highlight the mechanistic intricacies and underscore the importance of DNA nicks in facilitating transcription elongation by actively managing torsional stress during DNA unwinding by the RNA polymerase.

Performing spectroscopic measurements on biomolecules labeled with fluorescent probes is a powerful approach to locating the molecular behavior and dynamics of large systems at specific sites within their local environments. The indocarbocyanine dye Cy3 has emerged as one of the most commonly used chromophores. The incorporation of Cy3 dimers into DNA enhances experimental resolution owing to the spectral characteristics influenced by the geometric orientation of excitonically coupled monomeric units. Various theoretical models and simulations have been utilized to aid in the interpretation of the experimental spectra. In this study, we employ all-atom molecular dynamics simulations to study the structural dynamics of Cy3 dimers internally linked to the dsDNA backbone. We used quantum mechanical calculations to derive insights from both the linear absorption spectra and the circular dichroism data. Furthermore, we explore potential limitations within a commonly used force field for cyanine dyes. The molecular dynamics simulations suggest the presence of four possible Cy3 dimeric populations. The spectral simulations on the four populations show one of them to agree better with the experimental signatures, suggesting it to be the dominant population. The relative orientation of Cy3 in this population compares very well with previous predictions from the Holstein-Frenkel Hamiltonian model.

The DNA binding and cellular uptake of the lambda enantiomer of two bis-tetraazaphenanthrene (TAP) Ru(II) polypyridyl complexes containing either a linear dppn ( 1 ) or a hooked bdppz ( 2 ) benzodipyridophenazine ligand are reported, and the role of different charge-transfer states of the structural isomers in the photo-oxidation of guanine is explored. Both complexes possess characteristic metal-to-ligand charge-transfer (MLCT) bands between 400 and 500 nm and emission at ca. 630 nm in an aerated aqueous solution. Transient visible absorption (TrA) spectroscopy reveals that 400 nm excitation of 1 yields a dppn-based metal-to-ligand charge-transfer (MLCT) state, which in turn populates a dppn intraligand ( 3 IL) state. In contrast, photoexcitation of 2 results in an MLCT state on the TAP ligand and not the intercalating bdppz ligand. Both 1 and 2 bind strongly to double-stranded guanine-rich DNA with a loss of emission. Combined TrA and time-resolved infrared (TRIR) spectroscopy confirms formation of the guanine radical cation when 2 is bound to the d(G 5 C 5 ) 2 duplex, which is not the case when 1 is bound to the same duplex and indicates a different mechanism of action in DNA. Utilizing the long-lived triplet excited lifetime, we show good uptake and localization of 2 in live cells as well as isolated chromosomes. The observed shortening of the excited-state lifetime of 2 when internalized in cell chromosomes is consistent with DNA binding and luminescent quenching due to guanine photo-oxidation.

The protein-induced fluorescence change technique was employed to investigate the interactions between proteins and their DNA substrates modified with the Cy3 fluorophore. It has been reported that the human hepatoma-derived growth factor (HDGF), containing the chromatin-associated N-terminal proline-tryptophan-tryptophan-proline (PWWP) domain (the N-terminal 100 amino acids of HDGF) capable of binding the SMYD1 promoter, participates in various cellular processes and is involved in human cancer. This project investigated the specific binding behavior of HDGF, the PWWP domain, and the C140 domain (the C-terminal 140 amino acids of HDGF) sequentially using protein-induced fluorescence change. We found that the binding of HDGF and its related proteins on Cy3-labeled 15 bp SMYD1 dsDNA will cause a significant decrease in the recorded Cy3 fluorophore intensity, indicating the occurrence of protein-induced fluorescence quenching. The dissociation equilibrium constant was determined by fitting the bound fraction curve to a binding model. An approximate 10-time weaker SMYD1 binding affinity of the PWWP domain was found in comparison to HDGF. Moreover, the PWWP domain is required for DNA binding, and the C140 domain can enhance the DNA binding affinity. Furthermore, we found that the C140 domain can regulate the sequence-specific binding capability of HDGF on SMYD1 .

This study investigates the role of alkali cations in modulating the interaction between deoxyribonucleic acid (DNA) and Thioflavin T (ThT) in dilute and condensed phases. The emission characteristics of ThT were analyzed in the presence of double-stranded DNA and G-quadruplex structures with a focus on the effects of four cations: sodium, potassium, calcium, and magnesium. The ThT emission in double-stranded DNA was influenced by direct DNA binding and steric hindrance within the hydration shell of DNA, which was modulated by the presence of alkali cations. Lasing spectroscopy experiments further highlighted ThT sensitivity to the spatial arrangement of water molecules in the DNA hydration shell. Lasing was exclusively observed in the presence of Mg 2+ in the G-quadruplex structure, suggesting that the parallel propeller configuration of G4 provides an optimal environment for ThT light amplification. This study highlights the critical role of cations in DNA-dye interactions and reaffirms the significance of ThT in biophysical studies.

Tethered Particle Motion (TPM) is a single molecule technique, which consists in tracking the motion of a nanoparticle (NP) immersed in a fluid and tethered to a glass surface by a DNA molecule. The present work addresses the question of the applicability of TPM to fluids which contain crowders at volume fractions ranging from that of the nucleoid of living bacteria (around 30%) up to the jamming threshold (around 66%). In particular, we were interested in determining whether TPM can be used to characterize the compaction of DNA by globular crowders. To this end, extensive Brownian Dynamics simulations were performed with a specifically built coarse-grained model. Analysis of the simulations reveals several effects not observed in dilute media, which impose constraints on the TPM setup. In particular, the Tethered Fluorophore Motion (TFM) technique, which consists in replacing the NP by a much smaller fluorophore, is probably better suited than standard TPM. Moreover, a sample preparation technique which does not involve hydrophilic patches may be required. Finally, the use of a DNA brush may be needed to achieve DNA concentrations close to in vivo ones.

Many interactions between small molecules and particles occur in solutions. They are surrounded by other molecules that do not react, for example, biological processes in water, chemical reactions in gas or liquid solutions, and environmental reactions in air and water. However, predicting the rate of such diffusive interactions remains challenging, due to the random motion of molecules in solutions, as exampled by the famous Brownian motion of pollen particles. In this report, we experimentally confirmed that a disruptive rate equation we have published before can predict the association rate of typical adsorption at interfaces, which has a surprising fraction order of 4/3 that has not been considered before. This could be an important step toward a generalized method to predict the adsorption rate of many reactions.

The i-motif, a secondary structure of a four-helix formed by cytosine-rich DNA (i-DNA) through C-C + base pairing, is prevalent in human telomeres and promoters. This structure creates steric hindrance, thereby inhibiting both gene expression and protein coding. The conformation of i-DNA is intricately linked to the intracellular ionic environment. Hence, investigating its conformation under various ion conditions holds significant importance. In this study, we explored the impact of cations on the i-motif structure at the single-molecule level using the α-hemolysin (α-HL) nanochannel. Our findings reveal that the ability of i-DNA to fold into the i-motif structure follows the order Cs + > Na + > K + > Li + for monovalent cations. Furthermore, we observed the interconversion of single-stranded DNA (ss-DNA) and the i-motif structure at high and low concentrations of Mg 2+ and Ba 2+ electrolyte solutions. This study not only has the potential to extend the application of i-motif-based sensors in complex solution environments but also provides a new idea for the detection of metal ions.

Protein labeling through transient and repetitive hybridization of short, fluorophore-labeled DNA oligonucleotides has become widely applied in various optical super-resolution microscopy methods. The main advantages are multitarget imaging and molecular quantification. A challenge is the high background signal originating from the presence of unbound fluorophore-DNA labels in solution. Here, we report the self-quenching of fluorophore dimers conjugated to DNA oligonucleotides as a general concept to reduce the fluorescence background. Upon hybridization, the fluorescence signals of both fluorophores are restored. We expand the toolbox of fluorophores suitable for self-quenching and report their spectra and hybridization equilibria. We apply self-quenched fluorophore-DNA labels to stimulated emission depletion microscopy and single-molecule localization microscopy and report improved imaging performances.

Berubicin, a chemotherapy medication belonging to the class of anthracyclines, is simulated in double-stranded DNA sequences and cyclodextrins in an aqueous environment via full-atom molecular dynamics simulations on the time scale of microseconds. The drug is studied in both the neutral and protonated states so as to better comprehend the role of its charge in the formed complexes. The noncovalent berubicin-DNA and berubicin-cyclodextrin complexes are investigated in detail, paying special attention to their thermodynamic description by employing the double decoupling method, the solvent balance method, the weighted solvent accessible surface model, and the linear interaction energy method. A novel approach for extracting the desolvation thermodynamics of the binding process is also presented. Both the binding and desolvation Gibbs energies are decomposed into entropic and enthalpic contributions so as to elucidate the nature of complexation and its driving forces. Selected structural and geometrical properties of all the complexes, which are all stable, are analyzed. Both cyclodextrins under consideration are widely utilized for drug delivery purposes, and a comparative investigation between their bound states with berubicin is carried out.

Centromeric chromatin is a subset of chromatin structure and governs chromosome segregation. The centromere is composed of both CENP-A nucleosomes (CENP-A nuc ) and H3 nucleosomes (H3 nuc ) and is enriched with alpha-satellite (α-sat) DNA repeats. These CENP-A nuc have a different structure than H3 nuc , decreasing the base pairs (bp) of wrapped DNA from 147 bp for H3 nuc to 121 bp for CENP-A nuc . All these factors can contribute to centromere function. We investigated the interaction of H3 nuc and CENP-A nuc with NF-κB, a crucial transcription factor in regulating immune response and inflammation. We utilized atomic force microscopy (AFM) to characterize complexes of both types of nucleosomes with NF-κB. We found that NF-κB unravels H3 nuc , removing more than 20 bp of DNA, and that NF-κB binds to the nucleosomal core. Similar results were obtained for the truncated variant of NF-κB comprised only of the Rel homology domain and missing the transcription activation domain (TAD), suggesting that RelA TAD is not critical in unraveling H3 nuc . By contrast, NF-κB did not bind to or unravel CENP-A nuc . These findings with different affinities for two types of nucleosomes to NF-κB may have implications for understanding the mechanisms of gene expression in bulk and centromere chromatin.

Tethered motion is ubiquitous in nature, offering controlled movement and spatial constraints to otherwise chaotic systems. The enhanced functionality and practical utility of tethers has been exploited in biotechnology, catalyzing the design of novel biosensors and molecular assembly techniques. While notable technological advances incorporating tethered motifs have been made, a theoretical gap persists within the paradigm, hindering a comprehensive understanding of tethered-based technologies. In this work, we focus on the characterization of the binding kinetics of two tethered molecules functionalized to a hard surface. Using a mean-field approximation, the binding time of such bimolecular system is determined analytically. Furthermore, estimates of the grafting site separation and polymer lengths which expedite binding are provided. These estimates, along with the analytical theories and frameworks established here, have the potential to improve efficacy in self-assembly methods in DNA nanotechnology and can be extended to more biologically specific endeavors including targeted drug-delivery and molecular sensing.

The surface patterning in natural systems has exhibited appreciable functional advantages for life activities, which serve as inspiration for the design of artificial counterparts to achieve functions such as directional liquid transport at the nanoscale. Here, we propose a patterned two-dimensional (2D) in-plane heterostructure with a triangle-shaped hexagonal boron nitride (hBN) track embedded in graphene nanosheets, which can achieve unidirectional and self-propelled transport of nanodroplets carrying various biomolecules such as DNA, RNA, and peptides. Our extensive MD simulations show that the wettability gradient on the patterned heterostructure can drive the motion of nanodroplet with an instantaneous acceleration, which also permits long-distance transport (>100 nm) at the microsecond time scale. The different behaviors of various types of biomolecules have been further studied systematically within the transporting nanodroplets. These findings suggest that these specially designed, patterned heterostructures have the potential for spontaneous, directional transport of important biomolecules, which might be useful in biosensing, drug delivery, and biomedical nanodevices.

DNA origami is a pioneering approach for producing complex 2- or 3-D shapes for use in molecular electronics due to its inherent self-assembly and programmability properties. The electronic properties of DNA origami structures are not yet fully understood, limiting the potential applications. Here, we conduct a theoretical study with a combination of molecular dynamics, first-principles, and charge transmission calculations. We use four separate single strand DNAs, each having 8 bases (4 × G 4 C 4 and 4 × A 4 T 4 ), to form two different DNA nanostructures, each having two helices bundled together with one crossover. We also generated double-stranded DNAs to compare electronic properties to decipher the effects of crossovers and bundle formations. We demonstrate that density of states and band gap of DNA origami depend on its sequence and structure. The crossover regions could reduce the conductance due to a lack of available states near the HOMO level. Furthermore, we reveal that, despite having the same sequence, the two helices in the DNA origami structure could exhibit different electronic properties, and electrode position can affect the resulting conductance values. Our study provides better understanding of the electronic properties of DNA origamis and enables us to tune these properties for electronic applications such as nanowires, switches, and logic gates.

We have studied the excited states and structural properties for the complexes of cytosine (dC) 10 chains with silver ions (Ag + ) in a wide range of the Ag + to DNA ratio ( r ) and pH conditions using circular dichroism, steady-state absorption, and fluorescence spectroscopy along with the ultrafast fluorescence upconversion technique. We also calculated vertical electronic transition energies and determined the nature of the corresponding excited states in some models of the cytosine-Ag + complexes. We show that (dC) 10 chains in the presence of silver ions form a duplex stabilized by C-Ag + -C bonds. It is also shown that the i-motif structure formed by (dC) 10 chains is destabilized in the presence of Ag + ions. The excited-state properties in the studied complexes depend on the amount of binding ions and the binding sites, which is supported by the calculations. In particular, new low-lying excited states appear when the second Ag + ion interacts with the O atom of cytosine in the C-Ag + -C pairs. A similar picture is observed in the case when one Ag + ion interacts with one cytosine via the N7 atom.

Deoxyribonucleic acid (DNA) stability is a prerequisite in many applications, ranging from DNA-based vaccines and data storage to gene therapy. However, the strategies to enhance DNA stability are limited, and the underlying mechanisms are poorly understood. Ionic liquids (ILs), molten salts of organic cations and organic/inorganic anions, are showing tremendous prospects in myriads of applications. With a judicious choice of constituent ions, the protic nature of ILs can be tuned. In this work, we investigate the relative stability of full-length genomic DNA in aqueous IL solutions of increasing protic nature. Our experimental measurements show that the protic ionic liquids (PILs) enhance the DNA melting temperature significantly while unaltering its native B-conformation. Molecular dynamics simulations and quantum mechanical calculation results suggest that the intramolecular Watson-Crick H-bonding in DNA remains unaffected and, in addition, the PILs induce stronger H-bonding networks in solution through their ability to make multiple intermolecular H-bonds with the nucleobases and among its constituent ions, thus aiding greater DNA stability. The detailed understanding obtained from this study could bring about the much-awaited breakthrough in improved DNA stability for its sustained use in the aforesaid applications!

The migration of an electron-loss center (hole) in calf thymus DNA to bisbenzimidazole ligands bound in the minor groove is followed by pulse radiolysis combined with time-resolved spectrophotometry. The initially observed absorption spectrum upon oxidation of DNA by the selenite radical is consistent with spin on cytosine (C), as the GC • pair neutral radical, followed by the spectra of oxidized ligands. The rate of oxidation of bound ligands increased with an increase in the ratio ( r ) ligands per base pair from 0.005 to 0.04. Both the rate of ligand oxidation and the estimated range of hole transfer (up to 30 DNA base pairs) decrease with the decrease in one-electron reduction potential between the GC • pair neutral radical of ca. 1.54 V and that of the ligand radicals ( E 0 ', 0.90-0.99 V). Linear plots of log of the rate of hole transfer versus r give a common intercept at r = 0 and a free energy change of 12.2 ± 0.3 kcal mol -1 , ascribed to the GC • pair neutral radical undergoing a structural change, which is in competition to the observed hole transfer along DNA. The rate of hole transfer to the ligands at distance, R , from the GC • pair radical, k 2 , is described by the relationship k 2 = k 0 exp(constant/ R ), where k 0 includes the rate constant for surmounting a small barrier.

Low-energy (<20 eV) electrons (LEEs) can resonantly interact with DNA to form transient anions (TAs) of fundamental units, inducing single-strand breaks (SSBs), and cluster damage, such as double-strand breaks (DSBs). Shape resonances, which arise from electron capture in a previously unfilled orbital, can induce only a SSB, whereas a single core-excited resonance (i.e., two electrons in excited orbitals of the field of a hole) has been shown experimentally to cause cluster lesions. Herein, we show from time-dependent density functional theory (TDDFT) that a core-excited resonance can produce a DSB, i.e., a single 5 eV electron can induce two close lesions in DNA. We considered the nucleotide with the G-C base pair (ds[5'-G-3']) as a model for electron localization in the DNA double helix and calculated the potential energy surfaces (PESs) of excited states of the ground-state TA of ds[5'-G-3'], which correspond to shape and core-excited resonances. The calculations show that shape TAs start at ca. 1 eV, while core-excited TAs occur only above 4 eV. The energy profile of each excited state and the corresponding PES are obtained by simultaneously stretching both C5'-O5' bonds of ds[5'-G-3']. From the nature of the PES, we find two dissociative (σ*) states localized on the PO 4 groups at the C5' sites of ds[5'-G-3']. The first σ* state at 1 eV is due to a shape resonance, while the second σ* state is induced by a core-excited resonance at 5.4 eV. As the bond of the latter state stretches and arrives close to the dissociation limit, the added electron on C transfers to C5' phosphate, thus demonstrating the possibility of producing a DSB with only one electron of ca. 5 eV.

Elucidating how damage impacts DNA dynamics is essential for understanding the mechanisms of damage recognition and repair. Many DNA lesions alter their propensities to form low-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 1 H chemical exchange saturation transfer (CEST) NMR experiment in measuring the dynamics of oxidatively damaged 8-oxoguanine (8OG) in the mutagenic 8OG syn ·A anti 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) nonmutagenic 8OG anti ·A anti . These exchange parameters were in quantitative agreement with values from 13 C off-resonance R 1ρ and CEST on the labeled partner adenine. The Watson-Crick-like 8OG syn ·A anti mismatch also rescued the kinetics of Hoogsteen motions at distant A-T base pairs, which the G·A mismatch had slowed down. The results lend further support for 8OG anti ·A anti as a minor conformational state of 8OG·A, reveal that 8OG damage can impact Hoogsteen dynamics at a distance, and demonstrate the utility of 1 H CEST for measuring damage-dependent dynamics in unlabeled DNA.

Cas1 and Cas2 are highly conserved proteins among the clustered regularly interspaced short palindromic repeat Cas (CRISPR-Cas) systems and play a crucial role in protospacer selection and integration. According to the double-forked CRISPR Cas1-Cas2 complex, we conducted extensive all-atom molecular dynamics simulations to investigate the protospacer DNA binding and recognition. Our findings revealed that single-point amino acid mutations in Cas1 or in Cas2 had little impact on the binding of the protospacer, both in the binding and precatalytic states. In contrast, multiple-point amino acid mutations, particularly G74A, P80L, and V89A mutations on Cas2 and Cas2' proteins (m-multiple1 system), significantly affected the protospacer binding and selection. Notably, mutations on Cas2 and Cas2' led to an increased number of hydrogen bonds (#HBs) between Cas2&Cas2' and the dsDNA in the m-multiple1 system compared with the wild-type system. And the strong electrostatic interactions between Cas1-Cas2 and the protospacer DNA (psDNA) in the m-multiple1 system again suggested the increase in the binding affinity of protospacer acquisition. Specifically, mutations in Cas2 and Cas2' can remotely make the protospacer adjacent motif complementary (PAMc) sequences better in recognition by the two active sites, while multiple mutations K211E, P202Q, P212L, R138L, V134A, A286T, P282H, and P294H on Cas1a/Cas1b&Cas1a'/Cas1b' (m-multiple2 system) decrease the #HBs and the electrostatic interactions and make the PAMc worse in recognition compared with the wild-type system.

In nature, DNA exists primarily in a highly compacted form. The compaction of DNA in vivo is mediated by cationic proteins: histones in somatic nuclei and protamines in sperm chromatin. The extreme, nearly crystalline packaging of DNA by protamines in spermatozoa is thought to be essential for both efficient genetic delivery as well as DNA protection against damage by mutagens and oxidative species. The protective role of protamines is required in sperm, as they are sensitive to ROS damage due to the progressive loss of DNA repair mechanisms during maturation. The degree to which DNA packaging directly relates to DNA protection in the condensed state, however, is poorly understood. Here, we utilized different polycation condensing agents to achieve varying DNA packaging densities and quantify DNA damage by free radical oxidation within the condensates. Although we see that tighter DNA packaging generally leads to better protection, the length of the polycation also plays a significant role. Molecular dynamics simulations suggest that longer polyarginine chains offer increased protection by occupying more space on the DNA surface and forming more stable interactions. Taken together, our results suggest a complex interplay among polycation properties, DNA packaging density, and DNA protection against free radical damage within condensed states.

The basic packaging unit of eukaryotic chromatin is the nucleosome that contains 145-147 base pair duplex DNA wrapped around an octameric histone protein. While the DNA sequence plays a crucial role in controlling the positioning of the nucleosome, the molecular details behind the interplay between DNA sequence and nucleosome dynamics remain relatively unexplored. This study analyzes this interplay in detail by performing all-atom molecular dynamics simulations of nucleosomes, comparing the human α-satellite palindromic (ASP) and the strong positioning "Widom-601" DNA sequence at time scales of 12 μs. The simulations are performed at salt concentrations 10-20 times higher than physiological salt concentrations to screen the electrostatic interactions and promote unwrapping. These microsecond-long simulations give insight into the molecular-level sequence-dependent events that dictate the pathway of DNA unwrapping. We find that the "ASP" sequence forms a loop around SHL ± 5 for three sets of simulations. Coincident with loop formation is a cooperative increase in contacts with the neighboring N-terminal H2B tail and C-terminal H2A tail and the release of neighboring counterions. We find that the Widom-601 sequence exhibits a strong breathing motion of the nucleic acid ends. Coincident with the breathing motion is the collapse of the full N-terminal H3 tail and formation of an α-helix that interacts with the H3 histone core. We postulate that the dynamics of these histone tails and their modification with post-translational modifications (PTMs) may play a key role in governing this dynamics.

2-Aminopurine (2AP) is the most widely used fluorescent nucleobase analogue in DNA and RNA research. Its unique photophysical properties and sensitivity to environmental changes make it a useful tool for understanding nucleic acid dynamics and DNA-protein interactions. We studied the effect of ions present in commonly used buffer solutions on the excited-state photophysical properties of 2AP. Fluorescence quenching was negligible for tris(hydroxymethyl)aminomethane (TRIS), but significant for phosphate, carbonate, 3-(N-morpholino) propanesulfonic acid (MOPS), and 4-(2-hydroxyethyl)-1-piperazineethanesulfonic acid (HEPES) buffers. Results indicate that the two tautomers of 2AP (7H, 9H) are quenched by phosphate ions to different extents. Quenching by the H 2 PO 4 - ion is more pronounced for the 7H tautomer, while the opposite is true for the HPO 4 2- ion. For phosphate ions, the results of the time-resolved fluorescence study cannot be explained using a simple collisional quenching mechanism. Instead, results are consistent with transient interactions between 2AP and the phosphate ions. We postulate that excited-state interactions between the 2AP tautomers and an H-bond acceptor (phosphate and carbonate) result in significant quenching of the singlet-excited state of 2AP. Such interactions manifest in biexponential fluorescence intensity decays with pre-exponential factors that vary with quencher concentration, and downward curvatures of the Stern-Volmer plots.

The transport behavior of biomolecules at the confined nanoscale is very different from that of the bulk state. Numerous disease diagnostics and targeted drug treatments are performed based on nanochannels in cells. The specific structure and shape of nanochannels play an important role in the behavior and efficiency of substance transport. In this paper, we fabricated nanopores with different tilt angles and the same diameters using focused ion beam. The capture frequency and the blocking current amplitude of λ-DNA within large-angle nanopores decrease obviously, suggesting an increase in the energy barrier of large-angle nanopores and the fact that they stretch biomolecules to thinness. Most importantly, large-angle nanopores slow down λ-DNA transport by 2-4 times. MD simulations find that the sloped electroosmotic flow inside the tilted nanopores is the main factor contributing to the transport phenomena. The increase in the capture time of biomolecules by nanopores assists in obtaining more biological information from the current trajectories. Our study provides a new understanding of substance transport in specially shaped nanopores, which can be instrumental in providing fresh inspiration and approaches to the biomedical field.

Despite its rigid structure, DNA is a remarkably flexible molecule. Flexibility is essential for biological functions (such as transcription and gene repair), which require large-amplitude structural changes such as bubble formation. The bubbles thus formed are required to have a certain stability of their own and survive long on the time scale of molecular motions. A molecular understanding of fluctuations leading to quasi-stable structures is not available. Through extensive atomistic molecular dynamics simulations, we identify a sequence of microscopic events that culminate in local bubble formation, which is initiated by base-pair (BP) opening, resulting from the cleavage of native BP hydrogen bonds (HBs). This is followed by the formation of mismatched BPs with non-native contacts. These metastable structures can either revert to their original forms or undergo a flipping transition to form a local bubble that can span across 3-4 BPs. A substantial distortion of the DNA backbone and a disruption of BP stacking are observed because of the structural changes induced by these local perturbations. We also explored how water helps in the entire process. A small number of water molecules undergo rearrangement to stabilize the intermediate states by forming HBs with DNA bases. Water thus acts as a lubricant that counteracts the enthalpic penalty suffered from the loss of native BP contacts. Although the process of bubble formation is reversible, the sequence of steps involved poses an entropic barrier, preventing it from easily retracing the path to the native state.

RNA plays critical roles in the transmission and regulation of genetic information and is increasingly used in biomedical and biotechnological applications. Functional RNAs contain extended double-stranded regions, and the structure of double-stranded RNA (dsRNA) has been revealed at high resolution. However, the dependence of the properties of the RNA double helix on environmental effects, notably temperature, is still poorly understood. Here, we use single-molecule magnetic tweezer measurements to determine the dependence of the dsRNA twist on temperature. We find that dsRNA unwinds with increasing temperature, even more than DNA, with Δ T wRNA = -14.4 ± 0.7°/(°C·kbp), compared to Δ T wDNA = -11.0 ± 1.2°/(°C·kbp). All-atom molecular dynamics (MD) simulations using a range of nucleic acid force fields, ion parameters, and water models correctly predict that dsRNA unwinds with rising temperature but significantly underestimate the magnitude of the effect. These MD data, together with additional MD simulations involving DNA and DNA-RNA hybrid duplexes, reveal a linear correlation between the twist temperature decrease and the helical rise, in line with DNA but at variance with RNA experimental data. We speculate that this discrepancy might be caused by some unknown bias in the RNA force fields tested or by as yet undiscovered transient alternative structures in the RNA duplex. Our results provide a baseline to model more complex RNA assemblies and to test and develop new parametrizations for RNA simulations. They may also inspire physical models of the temperature-dependent dsRNA structure.

DNA nanotechnology has enabled the creation of supramolecular machines, whose shape and function are inspired from traditional mechanical engineering as well as from biological examples. As DNA inherently is a highly charged biopolymer, the external application of electric fields provides a versatile, computer-programmable way to control the movement of DNA-based machines. However, the details of the electrohydrodynamic interactions underlying the electrical manipulation of these machines are complex, as the influence of their intrinsic charge, the surrounding cloud of counterions, and the effect of electrokinetic fluid flow have to be taken into account. In this work, we identify the relevant effects involved in this actuation mechanism by determining the electric response of an established DNA-based nanorobotic arm to varying design and operation parameters. Borrowing an approach from single-molecule biophysics, we determined the electrical torque exerted on the nanorobotic arms by analyzing their thermal fluctuations when oriented in an electric field. We analyze the influence of various experimental and design parameters on the "actuatability" of the nanostructures and optimize the generated torque according to these parameters. Our findings give insight into the physical processes involved in the actuation mechanism and provide general guidelines that aid in designing and efficiently operating electrically driven nanorobotic devices made from DNA.

Nile Blue (NB) is a red-emissive dye that is well-known for imaging and staining applications. In this work, we describe the interaction of NB with various types of G-quadruplexes belonging to different topologies, molecularities, and conformations. Using spectroscopic techniques, we have determined the preferential binding of NB to c-Myc G-quadruplex and the other aspects of its binding. Concentration- and temperature-dependent studies showed that NB exists in a dynamic equilibrium between monomeric and H-aggregated states, which could be modulated by the addition of external agents such as anionic surfactants. NB displayed differential self-assembly with different types of G-quadruplex and duplex DNAs modulating its dynamic equilibrium between the monomeric and H-aggregated states. Fluorescence-based displacement studies revealed a 1:1 binding stoichiometry upon interaction with c- Myc G-quadruplex and an association constant of K app = 6.7 × 10 6 M -1 . Circular dichroism studies indicated that NB does not cause changes in the overall conformation of either G-quadruplexes or duplexes; however, it does indicate nucleic acid-dependent self-assembly at higher concentrations. Heat capacity measurement showed a more negative change when compared to that in DNA duplex, indicating more burial of the polar surface area by NB to the G-quadruplex host.

DNA is damaged through various exogenous sources (e.g., automobile exhaust, tobacco smoke, and processed foods), which can yield diverse C8-dG bulky aryl adducts. Adducts are known to induce structural changes to DNA that can lead to various biological outcomes, ranging from cell death to diseases such as cancer. Unfortunately, the relationship between the chemical composition of the damaged product, the adducted DNA structure, and the biological consequences is not well understood, which limits the development of disease detection and prevention strategies. The present study uses density functional theory (DFT) calculations and quintuplicate 1 μs molecular dynamics (MD) simulations to characterize the structure of DNA containing 21 model C8-dG adducts that systematically differ in size (phenyl to pyrenyl), shape (α (2,3), β (3,4) fusion, or ring substitution), and nucleobase-aryl group linkage (N, O, and C-linked). DFT calculations reveal that the inherent structural features of the G nucleobase adducts are impacted by linker type and bulky moiety shape, but not size, with the conformational flexibility reducing with α-ring fusion and linker composition as N > O > C. These structural properties are maintained in nucleoside models, which also reveal an increased propensity for anti -to- syn rotation about the glycosidic bond with N < O < C linker type. Although these diverse chemical features do not influence the global structure of adducted DNA, the adducts differentially impact the conformation local to the adducted site, including the relative populations of structures with the bulky moiety in the major groove (B conformer) and intercalated (stacked) into the helix (S conformer). Specifically, while the smallest phenyl adducts favor the B conformation and the largest pyrenyl-derived adducts stabilize the S conformation, the B/S ratio decreases with an increase in ring size and N > O > C linker composition. The shape and size (length) of the adduct can further finetune the B/S ratio, with β-fused naphthyl or α-fused phenanthryl N-linked adducts and O or C-linked adducts containing ring substitution increasing the prevalence of the S adducted DNA conformation. Overall, this work uncovers the significant effect of bulky moiety size and linker type, as well as the lesser impact of aryl group shape, on adducted DNA structure, which suggests differential replication and repair outcomes, and thereby represents an important step toward rationalizing connections between the structure and biological consequences of diverse DNA adducts.

Adenine DNA glycosylase (MutY) is a monofunctional glycosylase, removing adenines (A) misinserted opposite 8-oxo-7,8-dihydroguanine (OG), a common product of oxidative damage to DNA. Through multiscale calculations, we decipher a detailed adenine excision mechanism of MutY that is consistent with all available experimental data, involving an initial protonation step and two nucleophilic displacement steps. During the first displacement step, N-glycosidic bond cleavage is accompanied by the attack of the carboxylate group of residue Asp144 at the anomeric carbon (C1'), forming a covalent glycosyl-enzyme intermediate to stabilize the fleeting oxocarbenium ion. After departure of the excised base, water nucleophiles can be recruited to displace Asp144, completing the catalytic cycle with retention of stereochemistry at the C1' position. The two displacement reactions are found to mostly involve the movement of the oxocarbenium ion, occurring with large charge reorganization and thus sensitive to the internal electric field (IEF) exerted by the polar protein environment. Intriguingly, we find that the negatively charged carboxylate group is a good nucleophile for the oxocarbenium ion, yet an unactivated water molecule is not, and that the electric field catalysis strategy is used by the enzyme to enable its unique double-displacement reaction mechanism. A strong IEF, pointing toward 5' direction of the substrate sugar ring, greatly facilitates the second displacement reaction at the expense of elevating the barrier of the first one, thereby allowing both reactions to occur. These findings not only increase our understanding of the strategies used by DNA glycosylases to repair DNA lesions, but also have important implications for how internal/external electric field can be applied to modulate chemical reactions.

Protein-DNA interactions play an important role in numerous biological functions within the living cell. In many of these interactions, the DNA helix is significantly distorted upon protein-DNA complex formation. The Hha I restriction-modification system is one such system, where the methylation target is flipped out of the helix when bound to the methyltransferase. However, the base flipping mechanism is not well understood. The dynamics of the binding site of the Hha I methyltransferase and endonuclease (underlined) within the DNA oligomer [d(G 1 A 2 T 3 A 4G5C6G7C8 T 9 A 10 T 11 C 12 )] 2 are studied using deuterium solid-state NMR (SSNMR). SSNMR spectra obtained from DNAs deuterated on the base of nucleotides within and flanking the [5'-GCGC-3'] 2 sequence indicate that all of these positions are structurally flexible. Previously, conformational flexibility within the phosphodiester backbone and furanose ring within the target sequence has been observed and hypothesized to play a role in the distortion mechanism. However, whether that distortion was occurring through an active or passive mechanism remained unclear. These NMR data demonstrate that although the [5'-GCGC-3'] 2 sequence is dynamic, the target cytosine is not passively flipping out of the double-helix on the millisecond-picosecond time scale. Additionally, although previous studies have shown that both the furanose ring and phosphodiester backbone experience a change in dynamics upon methylation, which may play a role in recognition and cleavage by the endonuclease, our observations here indicate that methylation has no effect on the dynamics of the base itself.

Methylation at the C5 position of cytosine, a naturally occurring epigenetic modification on DNA, shows a high correlation with mutational hotspots in disease such as skin cancer. Due to its essential biological relevance, numerous studies were devoted to confirming that the methylated sites favor the formation of the cyclobutane pyrimidine dimer (CPD), a well-known UV-induced lesion. However, photophysical and photochemical properties of dinucleotides and polynucleotides containing 5-methylcytosine ( 5m C) remain elusive. Herein, a charge transfer (CT) triplet state, generated via intersystem crossing (ISC) from a CT singlet state that enhanced after methylation on cytosine, is directly observed by using femtosecond transient absorption (TA) and time-resolved mid-infrared (TRIR) spectroscopy together with quantum chemical calculations for the first time in the T 5m C dimer. Such an ISC process is quenched due to limitations of the ground-state geometries in 5m C-containing single-strand oligomer d(T 5m C) 9 . This mechanistic information is important for understanding the early stage of triplet state-induced CPD formation in 5m C containing DNA.

Protein-DNA interactions are fundamental to many biological processes. Proteins must find their target site on a DNA molecule to perform their function, and mechanisms for target search differ across proteins. Especially challenging phenomena to monitor and understand are transient binding events that occur across two DNA target sites, whether occurring in cis or trans. Type IIS restriction endonucleases rely on such interactions. They play a crucial role in safeguarding bacteria against foreign DNA, including viral genetic material. BfiI, a type IIS restriction endonuclease, acts upon a specific asymmetric sequence, 5-ACTGGG-3, and precisely cuts both upper and lower DNA strands at fixed locations downstream of this sequence. Here, we present two single-molecule Förster resonance energy-transfer-based assays to study such interactions in a BfiI-DNA system. The first assay focuses on DNA looping, detecting both "Phi"- and "U"-shaped DNA looping events. The second assay only allows in trans BfiI-target DNA interactions, improving the specificity and reducing the limits on observation time. With total internal reflection fluorescence microscopy, we directly observe on- and off-target binding events and characterize BfiI binding events. Our results show that BfiI binds longer to target sites and that BfiI rarely changes conformations during binding. This newly developed assay could be employed for other DNA-interacting proteins that bind two targets and for the dsDNA substrate BfiI-PAINT, a useful strategy for DNA stretch assays and other super-resolution fluorescence microscopy studies.

RNA duplexes are relatively rare but play very important biological roles. As an end-product of template-based RNA replication, they also have key implications for hypothetical primitive forms of life. Unless they are specifically separated by enzymes, these duplexes denature upon a temperature increase. However, mechanistic and kinetic aspects of RNA (and DNA) duplex thermal denaturation remain unclear at the microscopic level. We propose an in silico strategy that probes the thermal denaturation of RNA duplexes and allows for an extensive conformational space exploration along a wide temperature range with atomistic precision. We show that this approach first accounts for the strong sequence and length dependence of the duplexes melting temperature, reproducing the trends seen in the experiments and predicted by nearest-neighbor models. The simulations are then instrumental at providing a molecular picture of the temperature-induced strand separation. The textbook canonical "all-or-nothing" two-state model, very much inspired by the protein folding mechanism, can be nuanced. We demonstrate that a temperature increase leads to significantly distorted but stable structures with extensive base-fraying at the extremities, and that the fully formed duplexes typically do not form around melting. The duplex separation therefore appears as much more gradual than commonly thought.

Electron/hole transfer mechanisms in DNA and polynucleotide structures continue to garner considerable interest as emerging charge-transport systems and molecular electronics. To shed mechanistic insight into these electronic properties, we carried out large-scale density functional theory (DFT) calculations (up to 650 atoms) to systematically analyze the structural and electron/hole transport properties of fully periodic single- and double-stranded DNA. We examined the performance of various exchange-correlation functionals (LDA, BLYP, B3LYP, and B3LYP-D) and found that single-stranded thymine (T) and cytosine (C) are predominantly hole conductors, whereas single-stranded adenine (A) and guanine (G) are better electron conductors. For double-stranded DNA structures, the periodic A-T and G-C electronic band structures undergo a significant renormalization, which causes hole transport to only occur on the A and G nucleobases. Our calculations (1) provide new benchmarks for periodic nucleobase structures using dispersion-corrected hybrid functionals with large basis sets and (2) highlight the importance of dispersion effects for obtaining accurate geometries and electron/hole mobilities in these extended systems.

Most G-quadruplex (G4)-targeting ligands reported so far contain planar heteroaromatic groups and can intercalate into adjacent base pairs of double-stranded DNA (dsDNA). However, quantitative data on the binding number γ (ligands/bp) of G4 ligands that intercalate into long dsDNA remain lacking, which are essential for understanding the selectivity of G4 ligands. Here, using a single-molecule stretching assay based on the lengthening of dsDNA, we analyzed the dissociation constants and the binding number of eight most commonly used G4 ligands that intercalate into dsDNA. We showed that five ligands (CX-5461, BRACO-19, RHPS4, TrisQ, and Phen-DC3) intercalate into dsDNA avidly ( K d = 0.5-2.1 μM, saturated γ > 0.2 ligands/bp), which was similar to the typical dsDNA intercalator EB. Two bisquinolines, PDS and 360A, showed moderate intercalation ability ( K d = 22.5 and 48.7 μM) and γ < 0.01 ligands/bp in the presence of 1 μM ligands. Porphyrin NMM showed no intercalative binding even at 200 μM. Molecular docking and molecular dynamics simulations were carried out to further evaluate the intercalative binding of these G4 ligands with dsDNA by calculating the binding energies and π-π stacking probability.

The link between the chemical stability of G-quadruplex (qDNA) structures and their roles in eukaryotic genomic maintenance processes has been an area of interest now for several decades. This Review seeks to demonstrate how single-molecule force-based techniques can provide insight into the mechanical stabilities of a variety of qDNA structures as well as their ability to interconvert between different conformations under conditions of stress. Atomic force microscopy (AFM) and magnetic and optical tweezers have been the primary tools used in these investigations and have been used to examine both free and ligand-stabilized G-quadruplex structures. These studies have shown that the degree of stabilization of G-quadruplex structures has a significant effect on the ability of nuclear machinery to bypass these roadblocks on DNA strands. This Review will illustrate how various cellular components including replication protein A (RPA), Bloom syndrome protein (BLM), and Pif1 helicases are capable of unfolding qDNA. Techniques such as single-molecule fluorescence resonance energy transfer (smFRET), often in conjunction with the aforementioned force-based techniques, have proven extremely effective at elucidating the factors underpinning the mechanisms by which these proteins unwind qDNA structures. We will provide insight into how single-molecule tools have facilitated the direct visualization of qDNA roadblocks and also showcase results obtained from experiments designed to examine the ability of G-quadruplexes to limit the access of specific cellular proteins normally associated with telomeres.

With the development of solid-phase synthesis and DNA nanotechnology, DNA-based drug delivery systems have seen large advancements over the past decades. By combining various drugs (small-molecular drugs, oligonucleotides, peptides, and proteins) with DNA technology, drug-grafted DNA has demonstrated great potential as a promising platform in recent years, in which complementary properties of both components have been discovered; for instance, the synthesis of amphiphilic drug-grafted DNA has enabled the production of DNA nanomedicines for gene therapy and chemotherapy. Through the design of linkages between drug and DNA parts, stimuli-responsiveness can be instilled, which has boosted the application of drug-grafted DNA in various biomedical applications such as cancer therapy. This review discusses the progress of various drug-grafted DNA therapeutic agents, exploring the synthetic techniques and anticancer applications afforded through the combination of drug and nucleic acids.

A variety of anticancer and antibacterial drugs target DNA as one of their primary intracellular targets. Understanding ligand-DNA interactions and developing new, promising bioactive molecules for clinical use are greatly aided by elucidating the interaction between small molecules and natural polymeric DNAs. Small molecules' ability to attach to and inhibit DNA replication and transcription provides more information on how drugs impact the expression of genes. Yohimbine has been broadly studied in pharmacological properties, while its binding mode to DNA has not been explicated so far. In this study, an attempt was made to explore the interaction between Yohimbine (YH) and calf thymus (CT-DNA) by using varying thermodynamics and in silico approaches. Minor hypochromic and bathochromic shifts of fluorescence intensity were observed, suggesting the binding of YH to CT-DNA. The Scatchard plot analysis using the McGhee-von Hipple method revealed noncooperative binding and affinities in the range of 10 5 M -1 . The binding stoichiometry value is 2:1 (2 molecules of YH were span by 1 base pair) and was determined by Job's plot. The thermodynamic parameters suggested exothermic binding, which was favored by negative enthalpy and positive entropy changes from both isothermal titration calorimetry and temperature-dependent fluorescence experiment. Salt-dependent fluorescence suggested that the interaction between the ligand and DNA was governed by nonpolyelectrolytic forces. Kinetics experiment confirmed the static type of quenching. The results of iodide quenching, urea denaturation assay, dye displacement, DNA melting, and in silico molecular docking (MD) suggested groove binding of YH to CT-DNA. Circular dichroism spectra confirmed minimal perturbation of CT-DNA with YH binding via groove region. Therefore, the groove binding mechanism of interaction was validated by biophysical techniques and in silico , MD approaches. The findings supported here may contribute to the development of new YH therapeutics possessing better efficacy and lesser side effects.

Cyclic GMP-AMP synthase (cGAS) has been widely investigated as a drug target for its crucial role in innate immunity. However, the inhibitors designed using mouse model were often shown to be ineffective for humans. This outcome indicates that the activation mechanisms of human and mouse cGAS (mcGAS) are different. The cGAS activation is achieved by dimerization via binding to DNA, the detailed mechanism of which, however, is not entirely clear. To investigate these mechanisms, molecular dynamics (MD) simulations were performed on several states of four types of cGAS, namely, the mcGAS, the wild-type and A- and C-type mutations of human cGAS (hcGAS). We find that sequence differences between hcGAS and mcGAS can directly affect the protein structure stability, especially that of the siteB domain. The sequence and structural differences also contribute to DNA-binding differences. In addition, the conformational fluctuations of cGAS are found to correlate with the regulation of catalytic capacity. More importantly, we illustrate that dimerization enhances the correlation among distant residues and significantly reinforces the allosteric signal transmission among the DNA-binding interfaces and the catalytic pocket, which facilitates rapid immune response to cytosolic DNA. We conclude that siteB domain plays a prominent role in mcGAS activation, while siteA domain is key to hcGAS activation.

CRISPR (clustered regularly interspaced short palindromic repeats)/CRISPR-associated protein (Cas9) has been widely used for gene editing. Not all guide RNAs can cleave the DNA efficiently remains a major challenge to CRISPR/Cas9-mediated genome engineering. Therefore, understanding how the Cas9 complex successfully and efficiently identifies specific functional targets through base-pairing has great implications for such applications. The 10-nt seed sequence at the 3' end of the guide RNA is critical to target recognition and cleavage. Here, through stretching molecular dynamics simulation, we studied the thermodynamics and kinetics of the binding-dissociation process of the seed base and the target DNA base with the Cas9 protein. The results showed that in the presence of Cas9 protein, the enthalpy change and entropy change in binding-dissociation of the seed base with the target are smaller than those without the Cas9 protein. The reduction of entropy penalty upon association with the protein resulted from the pre-organization of the seed base in an A-form helix, and the reduction of enthalpy change was due to the electrostatic attraction of the positively charged channel with the negative target DNA. The binding barrier coming from the entropy loss and the dissociation barrier resulting from the destruction of the base pair in the presence of Cas9 protein were lower than those without protein, which indicates that the seed region is crucial for efficiently searching the correct target by accelerating the binding rate and dissociating fast from the wrong target.

A single nucleotide variant present between two otherwise identical nucleic acids will have unexpected functional consequences frequently. Here, a neoteric single nucleotide variation (SNV) detection assay that integrates two complementary nanotechnology systems, nanoassembly technology and an ingenious nanopore biosensing platform, has been applied to this research. Specifically, we set up a detection system to reflect the binding efficiency of the polymerase and nanoprobe through the difference of nanopore signals and then explore the effect of base mutation at the binding site. In addition, machine learning based on support vector machines is used to automatically classify characteristic events mapped by nanopore signals. Our system reliably discriminates single nucleotide variants at binding sites, even possessing the recognition among transitions, transversions, and hypoxanthine (base I). Our results demonstrate the potential of solid-state nanopore detection for SNV and provide some ideas for expanding solid-state nanopore detection platforms.

The denaturation of DNA is a critical process in biology and has many biotechnological applications. We investigated the compaction of locally denatured DNA by a chemical denaturation agent, dimethyl sulfoxide (DMSO), using magnetic tweezers (MTs), atomic force microscopy (AFM), and dynamic light scattering (DLS). Our results show that DMSO not only is capable of denaturing DNA but also able to compact DNA directly. When the DMSO concentration is above 10%, DNA condensation occurs due to the reduction in the persistence length of DNA and excluded volume interactions. Meanwhile, locally denatured DNA is easily condensed by divalent cations, such as magnesium ions (Mg 2+ ), contrasting with no native DNA condensation by the classical divalent cations. Specifically, the addition of more than 3 mM Mg 2+ to a 5% DMSO solution leads to DNA condensation. The critical condensing force ( F C ) increases from 6.4 to 9.5 pN when the Mg 2+ concentration grows from 3 to 10 mM. However, F C decreases gradually with a further increase in Mg 2+ concentration. For 3% DMSO solution, above 30 mM Mg 2+ is needed to compact DNA and a weaker condensing force was measured. With increasing Mg 2+ concentration, the morphology of the DMSO partially denatured DNA complex changes from loosely random coils to a dense network structure, even forming a spherical condensation nucleus, and finally to a partially disintegrated network. These findings show that the elasticity of DNA plays an important role in its denaturation and condensation behavior.

While the double helix is the most stable conformation of DNA inside cells, its transient unwinding and subsequent partial separation of the two complementary strands yields an intermediate single-stranded DNA (ssDNA). The ssDNA is involved in all major DNA transactions such as replication, transcription, recombination, and repair. The process of DNA unwinding and translocation is shouldered by helicases that transduce the chemical energy derived from nucleotide triphosphate (NTP) hydrolysis to mechanical energy and utilize it to destabilize hydrogen bonds between complementary base pairs. Consequently, a comprehensive understanding of the molecular mechanisms of these enzymes is essential. In the last few decades, a combination of single-molecule techniques (force-based manipulation and visualization) have been employed to study helicase action. These approaches have allowed researchers to study the single helicase-DNA complex in real-time and the free energy landscape of their interaction together with the detection of conformational intermediates and molecular heterogeneity during the course of helicase action. Furthermore, the unique ability of these techniques to resolve helicase motion at nanometer and millisecond spatial and temporal resolutions, respectively, provided surprising insights into their mechanism of action. This perspective outlines the contribution of single-molecule methods in deciphering molecular details of helicase activities. It also exemplifies how each technique was or can be used to study the helicase action of RecD2 in recombination DNA repair.

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.

It is generally believed that DNA nick is an effective way to release stress in supercoiled DNA, resulting from the twisting motion that individual strands rotate around the axis of the DNA helix. Here, we use MD simulations based on the oxDNA model to investigate the relaxation of 336 bp supercoiled minicircular DNA under DNA nick. Our simulations show that stress release, characterized by the abrupt decrease in linking number, may be induced by two types of DNA motion depending on the nick position. Except for the twisting motion, there is a writhing motion, that is, double strands collectively rotating with one plectoneme removal, which may occur in the process of DNA relaxation with the nick position in the loop region. Moreover, the writhing motion is more likely to occur in the DNA with relatively high hardness, such as C-G pairs. Our simulation results uncover the relationship between structural transformation, stress release, and DNA motion during the dynamic process under DNA nick, indicating the influence of nick position on the relaxation of the supercoiled DNA.

Antimalarial action of a drug is closely associated with the interaction with the parasite's DNA. Hence, in this study, the interaction of an important antimalarial drug, chloroquine (CLQ), has been investigated with six different sequences of DNA having pure adenine (A)-thymine (T) and pure cytosine (C)-guanine (G) as well as mixed nucleobases to achieve the nucleobase level of information in the binding of antimalarial drug with DNA along with binding induced stabilization/destabilization of DNA using different spectroscopic methods and molecular dynamics simulation technique. Further, the experiments have been also performed with 4-amino-7-chloroquinoline (7CLQ), an analogue of CLQ, to understand the role of the quinoline ring and side chain of CLQ in the binding with different sequences of DNA. The binding efficiency of CLQ with any sequence of DNA is higher than 7CLQ suggesting that the presence of charge on CLQ plays a prominent role in DNA binding. The data suggest that the binding of drug as well as induced stabilization of DNA depends significantly on the nature as well as the arrangement of the nucleobases. In general, the binding of CLQ with pure CG DNA is higher than with pure AT DNA; moreover, it prefers an alternate order of CG/AT than continual nucleobases in duplex DNA. CLQ predominately accommodates in the minor groove of AT DNA and prefers to form hydrogen bond mostly with the adenine nucleobase. In contrast to AT DNA, CLQ intrudes into the both major and minor grooves, but it is primarily accommodated into the major groove of CG DNA. CLQ forms a hydrogen bond mainly with guanine in the major groove and cytosine in the minor groove of CG DNA which enhances the binding of CLQ compared to AT DNA as well as induces higher stabilization in CG DNA. The molecular level information obtained about the functional group responsible for the interaction of CLQ as well as the role of chemical nature of nucleobases along with its ordering on the binding of CLQ with DNA may be useful in comprehensive understanding of its action mechanism.

The adsorption of oxidatively damaged DNA onto ferromagnetic substrates was investigated. Both confocal fluorescence microscopy and quartz crystal microbalance methods show that the adsorption rate and the coverage depend on the magnetization direction of the substrate and the position of the damage site on the DNA relative to the substrate. SQUID magnetometry measurements show that the subsequent magnetic susceptibility of the DNA-coated ferromagnetic film depends on the direction of the magnetic field that was applied to the ferromagnetic film as the molecules were adsorbed. This study reveals that (i) the spin and charge polarization in DNA molecules is changed significantly by oxidative damage in the G bases and (ii) the rate of adsorption on a ferromagnet, as a function of the direction of the magnetic dipole of the surface, can be used as an assay to detect oxidative damage in the DNA.

Single-molecule DNA studies have improved our understanding of the DNAs' structure and their interactions with other molecules. A variety of DNA labeling dyes are available for single-molecule studies, among which the bis-intercalating dye YOYO-1 and mono-intercalating dye YO-PRO-1 are widely used. They have an extraordinarily strong affinity toward DNA and are bright with a high quantum yield (>0.5) when bound to DNAs. However, it is still not clear how these dyes behave in DNA molecules under higher ionic strength and strong buffer flow. Here, we have studied the effect of ionic strength and flow rate of buffer on their binding in single DNA molecules. The larger the flow rate and the higher the ionic strength, the faster the intercalated dyes are washed away from the DNAs. In the buffer with 1 M ionic strength, YOYO-1 and YO-PRO-1 are mostly washed away from DNA within 2 min of moderate buffer flow.