Power-Law Solvation Dynamics in G-Quadruplex DNA: Role of Hydration Dynamics on Ligand Solvation inside DNA

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.

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.

The Role of Hydrogen Bonding in the Preferential Solvation of 5-Aminoquinoline in Binary Solvent Mixtures

5-Aminoquinoline (5AQ) has been used as a fluorescent probe of preferential solvation (PS) in binary solvent mixtures in which the nonpolar component is diethyl ether and the polar component is protic (methanol) or aprotic (acetonitrile). Hence, the roles of solvent polarity and solute-solvent hydrogen bonding have been delineated. Positive deviations of spectral shifts from a linear dependence on the concentration of the polar component, signifying PS, are markedly more pronounced in case of the protic solvent. Solvation dynamics on a nanosecond time scale mark the formation of the solvation shell around the fluorescent probe. Time-resolved area-normalized emission spectra indicate the occurrence of the continuous solvation of the excited state when the polar component is acetonitrile. In contrast, two distinct states were observed when the polar component was methanol, the second state being the hydrogen bonded one. Translational diffusion is the rate-determining step for formation of the solvation shell. The time constant associated with it has been estimated from rise times observed in fluorescence transients monitored at the red end of the fluorescence spectra and also from the time evolution of the spectral width of time-resolved emission spectra.

To unravel the connection between the non-equilibrium and equilibrium solvation dynamics of tryptophan: success and failure of the linear response theory of fluorescence Stokes shift

The connections between the non-equilibrium solvation dynamics upon optical transitions and the system's equilibrium fluctuations are explored in aqueous liquid. Linear response theory correlates time-dependent fluorescence with the equilibrium time correlation functions. In the previous work [T. Li, J. Chem. Theory Comput., 2017, 13, 1867], Stokes shift was explicitly decomposed into the contributions of various order time correlation functions on the excited state surface. Gaussian fluctuations of the solute-solvent interactions validate linear response theory. Correspondingly, the deviation of the Gaussian statistics causes the inefficiency of linear response evaluation. The above mechanism is thoroughly tested in this study. By employing molecular simulations, multiple non-equilibrium processes, not necessarily initiated from the ground state equilibrium minimum, were examined for tryptophan. Both the success and failure of linear response theory are found for this simple system and the mechanism is analyzed. These observations, assisted by the width dynamics, the initial state linear response approach, and the variation of the solvation structures, integrally verify the virtue of the excited state Gaussian statistics on the dynamics of Stokes shift.

Experimental and theoretical investigation of ground state intramolecular proton transfer (GSIPT) in salicylideneaniline Schiff base derivatives in polar protic medium

Ground state intramolecular proton transfer process has been comprehensively investigated in three salicylideneaniline Schiff base derivatives (SB1, SB2, and SB3) using experimental and theoretical methods. It has been confirmed that all the three Schiff base molecules in the ground electronic state exist in the enol form in non-polar and polar aprotic solvents. Keto form is being populated by the polar protic solvent through ground state intramolecular proton transfer (GSIPT) process. Ground state equilibrium between the enol and keto tautomers for SB1 and SB3 is mainly governed by the proton donating ability of the solvent. Ground state equilibria between the enol and keto tautomers of SB2 which is a positional isomer of SB3 is governed by the polarity and proton donating ability of the solvents. Excited state intramolecular proton transfer (ESIPT) process is also evidenced in all the three Schiff base molecules. Theoretical calculations at the B3LYP/cc-pVDZ level in the gas phase and in different solvents using polarisable continuum model (PCM) have failed to establish the GSIPT process. Microsolvation of individual enol and keto conformers has been investigated considering upto three solvent molecules. The energetics of the individual conformers together with the corresponding transition state have been calculated. It has been confirmed that the keto conformer is more stable compared to the enol conformer in microsolvated cluster of three methanol molecules. Lowering of activation energy for the enol to keto tautomerisation in the presence of methanol also supports the experimental observation for GSIPT process. TDDFT/B3LYP/cc-pVDZ single point calculations for microsolvated clusters of enol and keto form of the Schiff base molecules exhibit an excellent agreement with the experimentally obtained absorption spectra. Difference in spectral nature of the Schiff base molecules has been explained using natural bond orbital (NBO) analysis. Quantum theory of atoms in molecules (QTAIM) has also been utilised to understand the GSIPT process in detail.

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.

Live Cell Microscopy: A Physical Chemistry Approach

Probing dynamics of intracellular components using physical chemistry techniques is a remarkable bottom-up approach for understanding the structures and functions of a biological cell. In this "Feature Article", we give an overview on local polarity, solvation, viscosity, acid-base property, red-ox processes (thiol-disulfide exchange), and gene silencing at selected intracellular components inside a live cell. Significant differences have been observed between cancer cells and their noncancer counterparts. We demonstrate that thiol-disulfide exchange, calcium oscillation, and gene silencing are manifested in time dependence of fluorescence intensity. We show that fluorescent gold nanoclusters may be used in drug delivery (e.g., doxorubicin) and selective killing of cancer cells. Further, we discuss dynamics and structural changes of DNA quadruplexes and i-motifs, induced by different external conditions (e.g., pH) and additives (e.g., K⁺ and other target specific small molecules). We demonstrate that peptidomimetic analogues have high specificity over double-stranded DNA for binding with i-motifs and G-quadruplexes. These results may have significant biological implications.

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.

Efficient Criterion To Evaluate Linear Response Theory in Optical Transitions

The role of the Gaussian statistics on the solvation dynamics upon the photon excitation of the chromophore is deeply explored. The linear response theory for the fluorescence Stokes shift is investigated. An analytical formulism is presented to recast Stokes shift into the contributions of the equilibrium time correlation functions of the solute-solvent interactions on the excited-state surface, and the latter is further reformed and depicted by the time relaxation of the moment. As the first application of the formulism in the molecular dynamics simulations, it is verified that the efficiency of the linear response theory relies on the Gaussian characteristics of the dominant moments in terms of the Stokes shift, which is identified by the same relaxation dynamics between those moments and the linear order one. The comparisons between the above observations on the linearity of Stokes shift and the explanations in the literature are discussed. The key finding is the development of explicit criterion to measure the appropriateness of applying linear response theory.

Effect of water-DNA interactions on elastic properties of DNA self-assembled monolayers

DNA-water interactions have revealed as very important actor in DNA mechanics, from the molecular to the macroscopic scale. Given the particularly useful properties of DNA molecules to engineer novel materials through self-assembly and by bridging organic and inorganic materials, the interest in understanding DNA elasticity has crossed the boundaries of life science to reach also materials science and engineering. Here we show that thin films of DNA constructed through the self-assembly of sulfur tethered ssDNA strands demonstrate a Young's modulus tuning range of about 10 GPa by simply varying the environment relative humidity from 0% up to 70%. We observe that the highest tuning range occurs for ssDNA grafting densities of about 3.5 × 10¹³ molecules/cm ², where the distance between the molecules maximizes the water mediated interactions between the strands. Upon hybridization with the complementary strand, the DNA self-assembled monolayers significantly soften by one order of magnitude and their Young's modulus dependency on the hydration state drastically decreases. The experimental observations are in agreement with molecular dynamics simulations.

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

G-quadruplex DNA (GqDNA) structures play an important role in many specific cellular functions and are promising anti-tumor targets for small molecules (ligands). Here, we measured the dynamic Stokes shift of a ligand (Hoechst) bound to parallel c-Myc (mPu22) GqDNA over five decades of time from 100 fs to 10 ns, and compared it with the previously reported dynamics of DAPI bound to antiparallel human telomeric (hTelo22) GqDNA (Pal et al 2015 J. Phys. Chem. Lett. 6 1754). Stokes shift data from fluorescence up-conversion and time-correlated single photon counting experiments was combined to cover the broad dynamic range. The results show that the solvation dynamics of Hoechst in parallel mPu22 GqDNA follow a power law relaxation, added to fast 2 ps exponential relaxation, from 100 fs to 10 ns, with only a subtle difference of power law exponents in the two ligand-GqDNA systems (0.06 in Hoechst-mPu22 compared to 0.16 in DAPI-hTelo22). We measured steady-state fluorescence spectra and time-resolved anisotropy decays which confirm the tight binding of Hoechst to parallel mPu22 with a binding constant of ~1  ×  10⁵ M^-1. The molecular docking of Hoechst in parallel GqDNA followed by a 50 ns molecular dynamics (MD) simulation on a Hoechst-GqDNA complex reveals that Hoechst binds to one of the outer G-tetrads by end-stacking near G13 and G4, which is different from the binding site of DAPI inside a groove of antiparallel hTelo22 GqDNA. Reconciling previous experimental and simulation results, we assign the 2 ps component to the hydration dynamics of only weakly perturbed water near mPu22 and the power law relaxation to the coupled motion of water and DNA (i.e. DNA backbone, unpaired bases and loops connecting G-tetrads) which come near the Hoechst inside parallel GqDNA.

Solvatochromic fluorene-linked nucleoside and DNA as color-changing fluorescent probes for sensing interactions

A nucleoside bearing a solvatochromic push-pull fluorene fluorophore (dCFL ) was designed and synthesized by the Sonogashira coupling of alkyne-linked fluorene 8 with 5-iodo-2'-deoxycytidine. The fluorene building block 8 and labeled nucleoside dCFL exerted bright fluorescence with significant solvatochromic effect providing emission maxima ranging from 421 to 544 nm and high quantum yields even in highly polar solvents, including water. The solvatochromism of 8 was studied by DFT and ADC(2) calculations to show that, depending on the polarity of the solvent, emission either from the planar or the twisted conformation of the excited state can occur. The nucleoside was converted to its triphosphate variant dCFLTP which was found to be a good substrate for DNA polymerases suitable for the enzymatic synthesis of oligonucleotide or DNA probes by primer extension or PCR. The fluorene-linked DNA can be used as fluorescent probes for DNA-protein (p53) or DNA-lipid interactions, exerting significant color changes visible even to the naked eye. They also appear to be suitable for time-dependent fluorescence shift studies on DNA, yielding information on DNA hydration and dynamics.

Fluorescence Dynamics in the Endoplasmic Reticulum of a Live Cell: Time-Resolved Confocal Microscopy

Fluorescence dynamics in the endoplasmic reticulum (ER) of a live non-cancer lung cell (WI38) and a lung cancer cell (A549) are studied by using time-resolved confocal microscopy. To selectively study the organelle, ER, we have used an ER-Tracker dye. From the emission maximum (λmaxem) of the ER-Tracker dye, polarity (i.e. dielectric constant, ϵ) in the ER region of the cells (≈500 nm in WI38 and ≈510 nm in A549) is estimated to be similar to that of chloroform (λmaxem =506 nm, ϵ≈5). The red shift by 10 nm in λmaxem in the cancer cell (A549) suggests a slightly higher polarity compared to the non-cancer cell (WI38). The fluorescence intensity of the ER-Tracker dye exhibits prolonged intermittent oscillations on a timescale of 2-6 seconds for the cancer cell (A549). For the non-cancer cell (WI38), such fluorescence oscillations are much less prominent. The marked fluorescence intensity oscillations in the cancer cell are attributed to enhanced calcium oscillations. The average solvent relaxation time (<τs >) of the ER region in the lung cancer cell (A549, 250±50 ps) is about four times faster than that in the non-cancer cell (WI38, 1000±50 ps).

New insight into probe-location dependent polarity and hydration at lipid/water interfaces: comparison between gel- and fluid-phases of lipid bilayers

Location dependent polarity and hydration probed by a new series of 4-aminophthalimide-based fluorescent molecules (4AP-Cn;n= 2–10, 12) show different behaviour at gel- and fluid-phase lipid/water interfaces.

Solvation Dynamics in Different Phases of the Lyotropic Liquid Crystalline System

Reverse hexagonal (HII) liquid crystalline material based on glycerol monooleate (GMO) is considered as a potential carrier for drugs and other important biomolecules due to its thermotropic phase change and excellent morphology. In this work, the dynamics of encapsulated water, which plays important role in stabilization and formation of reverse hexagonal mesophase, has been investigated by time dependent Stokes shift method using Coumarin-343 as a solvation probe. The formation of the reverse hexagonal mesophase (HII) and transformation to the L2 phase have been monitored using small-angle X-ray scattering and polarized light microscopy experiments. REES studies suggest the existence of different polar regions in both HII and L2 systems. The solvation dynamics study inside the reverse hexagonal (HII) phase reveals the existence of two different types of water molecules exhibiting dynamics on a 120-900 ps time scale. The estimated diffusion coefficients of both types of water molecules obtained from the observed dynamics are in good agreement with the measured diffusion coefficient collected from the NMR study. The calculated activation energy is found to be 2.05 kcal/mol, which is associated with coupled rotational-translational water relaxation dynamics upon the transition from "bound" to "quasi-free" state. The observed ∼2 ns faster dynamics of the L2 phase compared to the HII phase may be associated with both the phase transformation as well as thermotropic effect on the relaxation process. Microviscosities calculated from time-resolved anisotropy studies infer that the interface is almost ∼22 times higher viscous than the central part of the cylinder. Overall, our results reveal the unique dynamical features of water inside the cylinder of reverse hexagonal and inverse micellar phases.