Origin of Slow Solvation Dynamics in DNA: DAPI in Minor Groove of Dickerson-Drew DNA

Melting and Bubble Formation in a Double-Stranded DNA: Microscopic Aspects of Early Base-Pair Opening Events and the Role of Water

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.

Unusual similarity of DNA solvation dynamics in high-salinity crowding with divalent cations of varying concentrations

Double-stranded DNA bears the highest linear negative charge density (2e- per base-pair) among all biopolymers, leading to strong interactions with cations and dipolar water, resulting in the formation of a dense 'condensation layer' around DNA. Interactions involving proteins and ligands binding to DNA are primarily governed by strong electrostatic forces. Increased salt concentrations impede such electrostatic interactions - a situation that prevails in oceanic species due to their cytoplasm being enriched with salts. Nevertheless, how these interactions' dynamics are affected in crowded hypersaline environments remains largely unexplored. Here, we employ steady-state and time-resolved fluorescence Stokes shifts (TRFSS) of a DNA-bound ligand (DAPI) to investigate the static and dynamic solvation properties of DNA in the presence of two divalent cations, magnesium (Mg2+), and calcium (Ca2+) at varying high to very-high concentrations of 0.15 M, 1 M and 2 M. We compare the results to those obtained in physiological concentrations (0.15 M) of monovalent Na+ ions. Combining data from fluorescence femtosecond optical gating (FOG) and time-correlated single photon counting (TCSPC) techniques, dynamic fluorescence Stokes shifts in DNA are analysed over a broad range of time-scales, from 100 fs to 10 ns. We find that while divalent cation crowding strongly influences the DNA stability and ligand binding affinity to DNA, the dynamics of DNA solvation remain remarkably similar across a broad range of five decades in time, even in a high-salinity crowded environment with divalent cations, as compared to the physiological concentration of the Na+ ion. Steady-state and time-resolved data of the DNA-groove-bound ligand are seemingly unaffected by ion-crowding in hypersaline solution, possibly due to ions being mostly displaced by the DNA-bound ligand. Furthermore, the dynamic coupling of cations with nearby water may possibly contribute to a net-neutral effect on the overall collective solvation dynamics in DNA, owing to the strong anti-correlation of their electrostatic interaction energy fluctuations. Such dynamic scenarios may persist within the cellular environment of marine life and other biological cells that experience hypersaline conditions.

Bimodal 1/f Noise and Anticorrelation between DNA-Water and DNA-Ion Energy Fluctuations

The coupling between the conformational fluctuations of DNA and its surrounding environment, consisting of water and ions in solution, remains poorly understood and relatively less investigated as compared to proteins. Here, with the help of molecular dynamics simulations and statistical mechanical analyses, we explore the dynamical coupling among DNA, water, and counterions through correlations among respective energy fluctuations in both double- (ds-) and single-stranded (ss-) DNA solutions. Fluctuations in the collective DNA-water and DNA-ion interaction energies are found to be strongly anticorrelated across all the systems. The fluctuations of DNA self-energy, however, are weakly coupled to DNA-water and DNA-ion interactions in ds-DNA. An enhancement of the DNA-water coupling is observed in ss-DNA, where the system is less rigid. All the interaction energies exhibit 1/f noise in their energy power spectra with surprisingly prominent bimodality in the DNA-water and DNA-ion fluctuations. The nature of the energy spectra appears to be indifferent to the relative rigidity of the DNA. We discuss the role of the observed correlations in ion-water motions on a DNA duplex in the experimentally observed anomalous slow dielectric relaxation and solvation dynamics and in furthering our understanding of the DNA energy landscape.

Interplay of conformational relaxation and hydrogen bond dynamics in the excited states of fluorescent Schiff base anions

Time resolved fluorescence spectroscopic investigation of four Schiff base anions has established that their excited state dynamics is governed by several solvent properties: polarity, viscosity and hydrogen bond donating ability. With viscous protic solvents like glycerol, fluorescence lifetimes of anions have been found to be markedly longer than those in ethanol, implying that conformational relaxation of molecules plays a key role in their nonradiative relaxation. Surprisingly, the lifetimes in less viscous aprotic solvents, like acetonitrile, are found to be even longer. The only plausible rationalization of this observation is in the light of hydrogen bond-assisted nonradiative phenomena that are operative in protic solvents. This contention draws support from a time evolution of the emission in the red end of the spectrum in low to moderately hydrogen bond donating protic solvents, with regard to an absence of such a rise time in aprotic solvents and strongly hydrogen bond donating solvents, viz., 2,2,2-trifluoroethanol. Rudimentary quantum chemical calculations provide a preliminary idea about the nature of excited state hydrogen bond redistribution involved in the process.

From structure and dynamics to biomolecular functions: The ubiquitous role of solvent in biology

Biological activity requires a solvent that can provide a suitable environment, which satisfies the twin need for stability and the ability to change. Among all the solvents water plays the most important role. We review, analyze, and comment on recent works on the structure and dynamics of water around biomolecules and their role in specific biological functions. While studies in the past have focused on understanding the biomolecule-water interactions through a hydration layer; recently the attention has shifted towards understanding functions at a molecular level. Such a microscopic understanding clearly requires elucidation of detailed dynamical processes where solvent molecules play an important role. Finally, we comment on the advances made in understanding the role of water inside a biological cell.

Knowledge-based prediction of DNA hydration using hydrated dinucleotides as building blocks

Database-derived water probability densities around structurally and sequentially distinct DNA dinucleotide fragments reproduce the known hydration motifs, which thus can be used as building blocks to predict DNA hydration.

Water plays an important role in stabilizing the structure of DNA and mediating its interactions. Here, the hydration of DNA was analyzed in terms of dinucleotide fragments from an ensemble of 2727 nonredundant DNA chains containing 41 853 dinucleotides and 316 265 associated first-shell water molecules. The dinucleotides were classified into categories based on their 16 sequences and the previously determined structural classes known as nucleotide conformers (NtCs). The construction of hydrated dinucleotide building blocks allowed dinucleotide hydration to be calculated as the probability of water density distributions. Peaks in the water densities, known as hydration sites (HSs), uncovered the interplay between base and sugar-phosphate hydration in the context of sequence and structure. To demonstrate the predictive power of hydrated DNA building blocks, they were then used to predict hydration in an independent set of crystal and NMR structures. In ten tested crystal structures, the positions of predicted HSs and experimental waters were in good agreement (more than 40% were within 0.5 Å) and correctly reproduced the known features of DNA hydration, for example the ‘spine of hydration’ in B-DNA. Therefore, it is proposed that hydrated building blocks can be used to predict DNA hydration in structures solved by NMR and cryo-EM, thus providing a guide to the interpretation of experimental data and computer models. The data for the hydrated building blocks and the predictions are available for browsing and visualization at the website https://watlas.datmos.org/watna/.

Molecular Picture of the Effect of Cosolvent Crowding on Ligand Binding and Dispersed Solvation Dynamics in G-Quadruplex DNA

Understanding molecular interactions and dynamics of proteins and DNA in a cell-like crowded environment is crucial for predicting their functions within the cell. Noncanonical G-quadruplex DNA (GqDNA) structures adopt various topologies that were shown to be strongly affected by molecular crowding. However, it is unknown how such crowding affects the solvation dynamics in GqDNA. Here, we study the effect of cosolvent (acetonitrile) crowding on ligand (DAPI) solvation dynamics within human telomeric antiparallel GqDNA through direct comparison of time-resolved fluorescence Stokes shift (TRFSS) experiments and molecular dynamics (MD) simulations results. We show that ligand binding affinity to GqDNA is drastically affected by acetonitrile (ACN). Solvation dynamics probed by DAPI in GqDNA groove show dispersed dynamics from ∼100 fs to 10 ns in the absence and presence of 20% and 40% (v/v) ACN. The nature of dynamics remain similar in buffer and 20% ACN, although in 40% ACN, distinct dynamics is observed in <100 ps. MD simulations performed on GqDNA/DAPI complex reveal preferential solvation of ligand by ACN, particularly in 40% ACN. Simulated solvation time-correlation functions calculated from MD trajectories compare very well to the overall solvation dynamics of DAPI in GqDNA, observed in experiments. Linear response decomposition of simulated solvation correlation functions unfolds the origin of dispersed dynamics, showing that the slower dynamics is dominated by DNA-motion in the presence of ACN (and also by the ACN dynamics at higher concentration). However, water-DNA coupled motion controls the slow dynamics in the absence of ACN. Our data, thus, unravel a detailed molecular picture showing that though ACN crowding affect ligand binding affinity to GqDNA significantly, the overall dispersed solvation dynamics in GqDNA remain similar in the absence and the presence of 20% ACN, albeit with a small effect on the dynamics in the presence of 40% ACN due to preferential solvation of ligand by ACN.

Long-range DNA-water interactions

DNA functions only in aqueous environments and adopts different conformations depending on the hydration level. The dynamics of hydration water and hydrated DNA leads to rotating and oscillating dipoles that, in turn, give rise to a strong megahertz to terahertz absorption. Investigating the impact of hydration on DNA dynamics and the spectral features of water molecules influenced by DNA, however, is extremely challenging because of the strong absorption of water in the megahertz to terahertz frequency range. In response, we have employed a high-precision megahertz to terahertz dielectric spectrometer, assisted by molecular dynamics simulations, to investigate the dynamics of water molecules within the hydration shells of DNA as well as the collective vibrational motions of hydrated DNA, which are vital to DNA conformation and functionality. Our results reveal that the dynamics of water molecules in a DNA solution is heterogeneous, exhibiting a hierarchy of four distinct relaxation times ranging from ∼8 ps to 1 ns, and the hydration structure of a DNA chain can extend to as far as ∼18 Å from its surface. The low-frequency collective vibrational modes of hydrated DNA have been identified and found to be sensitive to environmental conditions including temperature and hydration level. The results reveal critical information on hydrated DNA dynamics and DNA-water interfaces, which impact the biochemical functions and reactivity of DNA.

Preferential Binding of Epirubicin Hydrochloride with Single Nucleotide Mismatched DNA and Subsequent Sequestration by a Mixed Micelle

Targeting mismatched base pairs containing DNA using small molecules and exploring the underlying mechanism involved during the binding interactions is one of the fundamental aspects of drug design. These molecules in turn are used in nucleic acid targeted therapeutics and cancer diagnosis. In this work, we systematically delineate the binding of the anticancer drug, epirubicin hydrochloride (EPR) with 20-mer duplex DNA, having both natural nucleobase pairing and thermodynamically least stable non-Watson-Crick base pairing. From the thermal denaturation studies, we observed that EPR can remarkably enhance the thermal stability of cytosine-cytosine (CC) and cytosine-thymine (CT) mismatched (MM) DNA over other 20-mer duplex DNA. From steady-state fluorescence spectroscopy and isothermal titration calorimetry studies, we concluded that EPR binds strongly with the mismatched duplex DNA through the intercalation binding mode. The interaction of EPR and duplex DNA has also been monitored at a single molecular resolution using fluorescence correlation spectroscopy (FCS). Dynamic quantitates such as diffusion coefficients and hydrodynamic radii obtained from an FCS study along with association and dissociation rate constants estimated from intensity time trace analyses further substantiate the stronger binding affinity of EPR to the thermally less stable mismatched DNA, formed by the most discriminating nucleobase (viz. cytosine). Additionally, we have shown that EPR can be sequestered from nucleic acids using a mixed micellar system of an anionic surfactant and a triblock copolymer. From thermal denaturation studies and circular dichroism spectroscopy, we found that the extent of drug sequestration depends on the binding affinity of EPR to the duplex DNA, and this mixed micellar system can be employed for the removal of excess drug in the case of a drug overdose.

Solvation Dynamics in Water. 4. On the Initial Regime of Solvation Relaxation

It is shown, by means of numerical and analytic work, that initial molecular momenta play little significant role in the initial fast solvation relaxation that follows electronic excitation of, and charge creation for, a standard model system of a solute in water. Instead, the nonequilibrium dynamics are predominantly described by noninertial "steering" by the torques directly generated by the newly created charge distribution. It is this process that largely overcomes inertia and drives the relaxation dynamics on a time scale of a few tens of femtoseconds in the key initial regime of the dynamics. These results are discussed in the context of commonly employed descriptions such as inertial, Gaussian, and underdamped dynamical behavior.