Dispersed dynamics of solvation in G-quadruplex DNA: comparison of dynamic Stokes shifts of probes in parallel and antiparallel quadruplex structures

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.

G-Tetrad-Selective Ligand Binding Kinetics in G-Quadruplex DNA Probed with Fluorescence Correlation Spectroscopy

Probing the kinetics of ligand binding to biomolecules is of paramount interest in biology and pharmacology. Measurements of such kinetic processes provide information on the rate-determining steps that control the binding affinity of ligands to biomolecules, thereby predicting the mechanism of the molecular interaction. In this context, ligand binding to G-quadruplex DNA (GqDNA) structures has attracted tremendous attention primarily because of their use in possible anticancer therapy. Although a large number of G-quadruplex-specific ligands have been proposed, probing the kinetics of G-tetrad-selective binding of (multiple) ligands within a G-quadruplex DNA (GqDNA) structure remains challenging. Most of the earlier studies focused on the thermodynamics of ligand binding; however, the kinetics of ligand association and dissociation with GqDNA, particularly binding of multiple ligands within a GqDNA structure, have not been explored. Here, we propose a simple fluorescence correlation spectroscopy-based method that measures the G-tetrad-selective association and dissociation rates of ligands within a GqDNA structure by correlating the fluorescence fluctuations of a site-specific (5' or 3' end-labeled) fluorophore (Cy3) in GqDNA due to quenching of Cy3 fluorescence, induced by the ligand binding to the G-tetrads. We show that well-known GqDNA ligands, BRACO19, TMPyP4, Hoechst 33258, and Hoechst 33342, have G-tetrad-selective association and dissociation rates, which suggest site-dependent variation of free energy barriers for binding/unbinding of the ligands with GqDNA. We also show that the measured kinetic rates depend not only on the G-tetrad site (5' vs 3' end) but also on the ligand and GqDNA structures.

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.

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

The measurement and understanding of collective solvation dynamics in DNA have vital biological implications, as protein and ligand binding to DNA can be directly controlled by complex electrostatic interactions of anionic DNA and surrounding dipolar water, and ions. Time-resolved fluorescence Stokes shift (TRFSS) experiments revealed anomalously slow solvation dynamics in DNA much beyond 100 ps that follow either power-law or slow multiexponential decay over several nanoseconds. The origin of such dispersed dynamics remains difficult to understand. Here we compare results of TRFSS experiments to molecular dynamics (MD) simulations of well-known 4',6-diamidino-2-phenylindole (DAPI)/Dickerson-Drew DNA complex over five decades of time from 100 fs to 10 ns to understand the origin of such dispersed dynamics. We show that the solvation time-correlation function (TCF) calculated from 200 ns simulation trajectory (total 800 ns) captures most features of slow dynamics as measured in TRFSS experiments. Decomposition of TCF into individual components unravels that slow dynamics originating from dynamically coupled DNA-water motion, although contribution from coupled water-Na⁺ motion is non-negligible. The analysis of residence time of water molecules around the probe (DAPI) reveals broad distribution from ∼6 ps to ∼3.5 ns: Several (49 nos.) water molecules show residences time greater than 500 ps, of which at least 14 water molecules show residence times of more than 1 ns in the first solvation shell of DAPI. Most of these slow water molecules are found to occupy two hydration sites in the minor groove near DAPI binding site. The residence time of Na⁺, however, is found to vary within ∼17-120 ps. Remarkably, we find that freezing the DNA fluctuations in simulation eliminates slower dynamics beyond ∼100 ps, where water and Na⁺ dynamics become faster, although strong anticorrelation exists between them. These results indicate that primary origin of slow dynamics lies within the slow fluctuations of DNA parts that couple with nearby slow water and ions to control the dispersed collective solvation dynamics in DNA minor groove.

Effect of T·T Mismatch on DNA Dynamics Probed by Minor Groove Binders: Comparison of Dynamic Stokes Shifts of Hoechst and DAPI

Recognition of DNA base mismatches and their subsequent repair by enzymes is vital for genomic stability. However, it is difficult to comprehend such a process in which enzymes sense and repair different types of mismatches with different ability. It has been suggested that the differential structural changes of mismatched bases act as cues to the repair enzymes, although the effect of such DNA structural changes on surrounding water and ion dynamics is inevitable due to strong electrostatic coupling among them. Thus, collective dynamics of DNA, water, and ions near the mismatch site is believed to be important for mismatch recognition and repair mechanism. Here we show that introduction of a T·T mismatch in the minor groove of DNA induces dispersed (collective) power-law solvation dynamics (of exponent ∼0.24), measured by monitoring the time-resolved fluorescence Stokes shifts (TRFSS) of two popular minor groove binders (Hoechst 33258 and DAPI) over five decades of time from 100 fs to 10 ns. The same ligands however sense different dynamics (power-law of exponent ∼0.15 or power-law multiplied with biexponential relaxation) in the minor groove of normal-DNA. The similar fluorescence anisotropy decays of ligands measured in normal- and T·T-DNA suggest that Stokes shift dynamics and their changes in T·T-DNA purely originate from the solvation process, and not from any internal rotational motion of probe-ligands. The dispersed power-law solvation dynamics seen in T·T-DNA indicate that the ligands do not sense any particular (exponential) relaxation specific to T·T wobbling and/or other conformational changes. This could be the reason why T·T mismatch is recognized by enzymes with lower efficiency compared to purine-pyrimidine and purine-purine mismatches.