Role of Monovalent Counterions in the Ultrafast Dynamics of DNA

Anomalous Concentration Dependence of Viscosity: Hidden Role of Cross-Correlations in Aqueous Electrolyte Solutions

The viscosity of aqueous electrolyte solutions exhibits well-known composition-dependent anomalies that show certain definitive trends and universal features. The viscosity of LiCl and NaCl solutions increases with concentration in a monotonic fashion, while solutions of KCl, RbCl, and CsCl exhibit a more complex behavior. Here, the viscosity first decreases and then increases with increasing concentration, with a rather broad minimum at intermediate concentrations (ca. 1-3 m). To unearth the origin of such puzzling behavior, we carried out detailed molecular-level analyses by interrogating the exact Green-Kubo expression of viscosity in terms of the stress-stress time correlation function (SS-TCF). The total SS-TCF can be decomposed into a collection of three self- and three cross-SS-TCFs arising from the three constituent components (water, cations, and anions). Mode coupling theory (MCT) analysis for the friction on ions and the viscosity of the solution suggests the possible importance of two-particle static and time-dependent cross-correlations between water and the ions. We calculate the viscosity and other dynamical properties for all five electrolyte (LiCl, NaCl, KCl, RbCl, and CsCl) solutions over a range of concentrations, using two models of water (SPC/E and TIP4P/2005). The total viscosity derives non-negligible contributions from all of the terms. The cross-correlations are found to be surprisingly large and seen to play a hidden role in the concentration dependence. However, the importance of cross-correlations is often not discussed. Our study leads to a theoretical understanding of the microscopic origin of the observed anomalies in the composition dependence of viscosity across all five electrolytes.

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.

Solvation and hydrogen bonding aided efficient non-radiative deactivation of polar excited state of 5-aminoquinoline

The mechanism of efficient fluorescence quenching of 5-aminoquinoline in alcoholic solvents.

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.

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.

A deeper insight into an intriguing acetonitrile–water binary mixture: synergistic effect, dynamic Stokes shift, fluorescence correlation spectroscopy, and NMR studies

An insight study reveals the strong synergistic solvation behaviours from reporter dye molecules within the acetonitrile (ACN)–water (WT) binary mixture.

Algebraic Statistics of Poincaré Recurrences in a DNA Molecule

The statistics of Poincaré recurrences is studied for the base-pair breathing dynamics of an all-atom DNA molecule in a realistic aqueous environment with thousands of degrees of freedom. It is found that at least over five decades in time the decay of recurrences is described by an algebraic law with the Poincaré exponent close to β=1.2. This value is directly related to the correlation decay exponent ν=β-1, which is close to ν≈0.15 observed in the time resolved Stokes shift experiments. By applying the virial theorem we analyze the chaotic dynamics in polynomial potentials and demonstrate analytically that an exponent β=1.2 is obtained assuming the dominance of dipole-dipole interactions in the relevant DNA dynamics. Molecular dynamics simulations also reveal the presence of strong low frequency noise with the exponent η=1.6. We trace parallels with the chaotic dynamics of symplectic maps with a few degrees of freedom characterized by the Poincaré exponent β~1.5.

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

G-quadruplex DNA (GqDNA) structures act as promising anticancer targets for small-molecules (ligands). Solvation dynamics of a ligand (DAPI: 4',6-diamidino-2-phenylindole) inside antiparallel-GqDNA is studied through direct comparison of time-resolved experiments to molecular dynamics (MD) simulation. Dynamic Stokes shifts of DAPI in GqDNA prepared in H2O buffer and D2O are compared to find the effect of water on ligand solvation. Experimental dynamics (in H2O) is then directly compared with the dynamics computed from 65 ns simulation on the same DAPI-GqDNA complex. Ligand solvation follows power-law relaxation (summed with fast exponential relaxation) from ~100 fs to 10 ns. Simulation results show relaxation below ~5 ps is dominated by water motion, while both water and DNA contribute comparably to dictate long-time power-law dynamics. Ion contribution is, however, found to be negligible. Simulation results also suggest that anomalous solvation dynamics may have origin in subdiffusive motion of perturbed water near GqDNA.

Probe Position-Dependent Counterion Dynamics in DNA: Comparison of Time-Resolved Stokes Shift of Groove-Bound to Base-Stacked Probes in the Presence of Different Monovalent Counterions

Time-resolved fluorescence Stokes shifts (TRFSS) of 4',6-diamidino-2-phenylindole (DAPI) inside the minor groove of DNA are measured in the presence of three different monovalent counterions: sodium (Na(+)), rubidium (Rb(+)), and tetrabutylammonium (TBA(+)). Fluorescence up-conversion and time-correlated single photon counting are combined to obtain the time-resolved emission spectra (TRES) of DAPI in DNA from 100 fs to 10 ns. Time-resolved Stokes shift data suggest that groove-bound DAPI can not sense the counterion dynamics because the ions are displaced by DAPI far from the probe-site. However, when these results are compared to the earlier base-stacked coumarin data, the same ions are found to affect the nanosecond dynamics significantly. This suggests that the ions come close to the probe-site, such that they can affect the dynamics when measured by base-stacked coumarin. These results support previous molecular dynamics (MD) simulation data of groove-bound and base-stacked probes inside DNA.

Effects of an unnatural base pair replacement on the structure and dynamics of DNA and neighboring water and ions

Incorporating small molecule probes into biomolecular systems to report on local structure and dynamics is a powerful strategy that underlies a wide variety of experimental techniques, including fluorescence, electron paramagnetic resonance (EPR), and Forster resonance energy transfer (FRET) measurements. When an unnatural probe is inserted into a protein or DNA, the degree to which the presence of the probe has perturbed the local structure and dynamics it was intended to study is always an important concern. Here, molecular dynamics (MD) simulations are used to systematically study the effect of replacing a DNA base pair with a fluorescent probe, coumarin 102 deoxyriboside, at six unique sites along an A-tract DNA dodecamer. While the overall structure of the DNA oligonucleotide remains intact, replacement of A*T base pairs leads to widespread structural and dynamic perturbations up to four base pairs away from the probe site, including widening of the minor groove and increased DNA flexibility. New DNA conformations, not observed in the native sequence, are sometimes found in the vicinity of the probe and its partner abasic site analog. Strong correlations are demonstrated between DNA surface topology and water mobility.

Computational exploration of mobile ion distributions around RNA duplex

Atomically detailed distributions of ions around an A-form RNA are computed. Different mixtures of monovalent and divalent ions are considered explicitly. Studies of tightly bound and of diffusive (but bound) ions around 25 base pairs RNA are conducted in explicit solvent. Replica exchange simulations provide detailed equilibrium distributions with moderate computing resources (20 ns of simulation using 64 replicas). The simulations show distinct behavior of single and double charged cations. Binding of Mg(2+) ion includes tight binding to specific sites while Na(+) binds only diffusively. The tight binding of Mg(2+) is with a solvation shell while Na(+) can bind directly to RNA. Negative mobile ions can be found near the RNA but must be assisted by proximate and mobile cations. At distances larger than 16 A from the RNA center, a model of RNA as charged rod in a continuum of ionic solution provides quantitative description of the ion density (the same as in atomically detailed simulation). At shorter distances, the structure of RNA (and ions) has a significant impact on the pair correlation functions. Predicted binding sites of Mg(2+) at the RNA surface are in accord with structures from crystallography. Electric field relaxation is investigated. The relaxation due to solution rearrangements is completed in tens of picoseconds, while the contribution of RNA tumbling continues to a few nanoseconds.

Dynamics of water and ions near DNA: comparison of simulation to time-resolved stokes-shift experiments

Time-resolved Stokes-shift experiments measure the dynamics of biomolecules and of the perturbed solvent near them on subnanosecond time scales, but molecular dynamics simulations are needed to provide a clear interpretation of the results. Here we show that simulations using standard methods quantitatively reproduce the main features of TRSS experiments in DNA and provide a molecular assignment for the dynamics. The simulations reproduce the magnitude and unusual power-law dynamics of the Stokes shift seen in recent experiments [ Andreatta, D., et al. J. Am. Chem. Soc. 2005, 127, 7270 ]. A polarization model is introduced to eliminate cross-correlations between the different components contributing to the signal. Using this model, well-defined contributions of the DNA, water, and counterion to the experimental signal are extracted. Water is found to have the largest contribution and to be responsible for the power-law dynamics. The counterions have a smaller, but non-negligible, contribution with a time constant of 220 ps. The contribution to the signal of the DNA itself is minor and fits a 30 ps stretched exponential. Both time-averaged and dynamic distributions are calculated. They show a small subset of ions with a different coupling but no other evidence of substates or rate heterogeneity.

Excitation Wavelength Dependence of Solvation Dynamics in a Supramolecular Assembly: PEO−PPO−PEO Triblock Copolymer and SDS

The triblock copolymer (PEO)20-(PPO)70-(PEO)20 (P123) forms a supramolecular aggregate with sodium dodecyl sulfate (SDS). The solvation dynamics and anisotropy decay of coumarin 480 (C480) in different regions of a P123-SDS aggregate are studied through variation of the excitation wavelength (lambdaex) using femtosecond upconversion. In a P123 micelle, because of the drastic differences in polarity between the hydrophilic corona region (PEO block) and the hydrophobic PPO core, C480 exhibits a pronounced red edge excitation shift (REES) of emission maximum by 24 nm. In the P123-SDS aggregate, SDS penetrates the core of the P123 micelle. This increases the polarity of the core and reduces the difference in the polarity between the core and the corona region. In a P123-SDS aggregate, the REES is much smaller (5 nm) which suggests a reduced difference between the core and the corona. Solvation dynamics in a P123 micelle displays a bulklike ultrafast component (<0.3 and 1 ps) in the PEO corona region, a 200 ps component arising from dynamics of polymer segments, and a very long component (5000 or 3000 ps) due to the highly restricted PPO core. In a P123-SDS aggregate, at lambdaex = 375 and 405 nm, the solvation dynamics is found to be faster than that in P123 micelle. In this case, the component (3000 ps) arising from the core region is faster than that (5000 ps) in P123 micelle. In both P123 micelle and P123-SDS aggregate, the relative contribution of the core region decreases and that of the corona region increases with an increase in lambdaex. At lambdaex = 435 nm, which probes the hydrophilic corona, the solvation dynamics for both P123 micelle and P123-SDS aggregate are almost similar.

Influence of the nature of quantum dot surface cations on interactions with DNA

Quantum dots are semiconductor nanoparticles that are approximately 1-10nm in diameter, similar to small proteins, and their photoluminescence is sensitive to the presence and nature of adsorbates. We have deployed these nanomaterials as luminescent probes of DNA structure. Sequence dependent conformational flexibility of DNA is of great interest due to its implications for drug-DNA and DNA-protein interactions. The counterion atmosphere surrounding DNA plays an important role in its structure, dynamics, and packaging. In this paper, we investigate the effect that various monovalent and divalent cations have on the binding of 4.5 nm CdS quantum dots to oligonucleotides that have sequence-directed intrinsic structure.