Effect of protein binding on ultrafast DNA dynamics: characterization of a DNA:APE1 complex

Environmentally sensitive fluorescent nucleoside analogues as probes for nucleic acid - protein interactions: molecular design and biosensing applications

Fluorescent nucleoside analogues (FNAs) are indispensable in studying the interactions of nucleic acids with nucleic acid-binding proteins. By replacing one of the poorly emissive natural nucleosides, FNAs enable real-time optical monitoring of the binding interactions in solutions, under physiologically relevant conditions, with high sensitivity. Besides that, FNAs are widely used to probe conformational dynamics of biomolecular complexes using time-resolved fluorescence methods. Because of that, FNAs are tools of high utility for fundamental biological research, with potential applications in molecular diagnostics and drug discovery. Here I review the structural and physical factors that can be used for the conversion of the molecular binding events into a detectable fluorescence output. Typical environmentally sensitive FNAs, their properties and applications, and future challenges in the field are discussed.

Structure and luminescence of DNA-templated silver clusters

DNA serves as a versatile template for few-atom silver clusters and their organized self-assembly. These clusters possess unique structural and photophysical properties that are programmed into the DNA template sequence, resulting in a rich palette of fluorophores which hold promise as chemical and biomolecular sensors, biolabels, and nanophotonic elements. Here, we review recent advances in the fundamental understanding of DNA-templated silver clusters (Ag N -DNAs), including the role played by the silver-mediated DNA complexes which are synthetic precursors to Ag N -DNAs, structure-property relations of Ag N -DNAs, and the excited state dynamics leading to fluorescence in these clusters. We also summarize the current understanding of how DNA sequence selects the properties of Ag N -DNAs and how sequence can be harnessed for informed design and for ordered multi-cluster assembly. To catalyze future research, we end with a discussion of several opportunities and challenges, both fundamental and applied, for the Ag N -DNA research community. A comprehensive fundamental understanding of this class of metal cluster fluorophores can provide the basis for rational design and for advancement of their applications in fluorescence-based sensing, biosciences, nanophotonics, and catalysis.

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.

Dynamics of water and ions around DNA: What is so special about them?

Water around biomolecules is special for behaving strangely - both in terms of structure and dynamics, while ions are found to control various interactions in biomolecules such as DNA, proteins and lipids. The questions that how water and ions around these biomolecules behave in terms of their structure and dynamics, and how they affect the biomolecular functions have triggered tremendous research activities worldwide. Such activities not only unfolded important static and dynamic properties of water and ions around these biomolecules, but also provoked heated debate regarding their explanation and role in biological functions. DNA, being negatively charged, interacts strongly with the surrounding dipolar water and positively charged counterions, leading to complex electrostatic coupling of water and ions with the DNA. Recent timeresolved fluorescence Stokes shift experiments and related computer simulation studies from our and other laboratories have unfolded some unique dynamic characteristics of water and ions near different structures of DNA. These results are discussed here to showcase the specialty of water and ion dynamics around DNA.

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.

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.

Design and activity of AP endonuclease-1 inhibitors

Apurinic/apyrimidinic endonuclease-1/redox effector factor-1 (APE-1) is a critical component of base excision repair that excises abasic lesions created enzymatically by the action of DNA glycosylases on modified bases and non-enzymatically by hydrolytic depurination/depyrimidination of nucleobases. Many anticancer drugs generate DNA adducts that are processed by base excision repair, and tumor resistance is frequently associated with enhanced APE-1 expression. Accordingly, APE-1 is a potential therapeutic target to treat cancer. Using computational approaches and the high resolution structure of APE-1, we developed a 5-point pharmacophore model for APE-1 small molecule inhibitors. One of the nM APE-1 inhibitors (AJAY-4) that was identified based on this model exhibited an overall median growth inhibition (GI50) of 4.19 μM in the NCI-60 cell line panel. The mechanism of action is shown to be related to the buildup of abasic sites that cause PARP activation and PARP cleavage, and the activation of caspase-3 and caspase-7, which is consistent with cell death by apoptosis. In a drug combination growth inhibition screen conducted in 10 randomly selected NCI-60 cell lines and with 20 clinically used non-genotoxic anticancer drugs, a synergy was flagged in the SK-MEL-5 melanoma cell line exposed to combinations of vemurafenib, which targets melanoma cells with V600E mutated BRAF, and AJAY-4, our most potent APE-1 inhibitor. The synergy between AJAY-4 and vemurafenib was not observed in cell lines expressing wild-type B-Raf protein. This synergistic combination may provide a solution to the resistance that develops in tumors treated with B-Raf-targeting drugs.

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.

Fluorescent DNA-based enzyme sensors

Fluorescent sensors that make use of DNA structures have become widely useful in monitoring enzymatic activities. Early studies focused primarily on enzymes that naturally use DNA or RNA as the substrate. However, recent advances in molecular design have enabled the development of nucleic acid sensors for a wider range of functions, including enzymes that do not normally bind DNA or RNA. Nucleic acid sensors present some potential advantages over classical small-molecule sensors, including water solubility and ease of synthesis. An overview of the multiple strategies under recent development is presented in this critical review, and expected future developments in microarrays, single molecule analysis, and in vivo sensing are discussed (160 references).

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.

Improved hydration-sensitive dual-fluorescence labels for monitoring peptide-nucleic acid interactions

Environmentally sensitive labels constitute a new, attractive tool for monitoring biomolecular interactions. 3-Hydroxychromone derivatives are of particular interest because they undergo excited-state intramolecular proton transfer (ESIPT) showing dual emission highly sensitive to environmental hydration. To overcome the drawbacks of the previously developed label for sensing protein-DNA interactions based on 2-furanyl-3-hydroxychromone (FC), a series of hydration-sensitive labels based on 3-hydroxy-4'-methoxyflavone have been synthesized. As compared to FC, the new labels display higher sensitivity of the ratio of their two emission bands (N*/T*) to solvent polarity and H-bond donor ability, as well as higher fluorescence quantum yields in water. Moreover, they show higher pK(a) values of their 3-hydroxyl group, allowing their application at neutral pH without interference of anionic forms. To illustrate the applications of these labels, we covalently coupled them to the N-terminus of the Tat(44-61) peptide that corresponds to the basic domain of the HIV-1 Tat protein. This coupling did not modify the nucleic acid chaperone properties of the peptide. Binding of oligonucleotides of varying length, sequence, and strandedness to the labeled peptides induced dramatic change in the N*/T* ratio of their two emission bands. This change indicated that the level of probe hydration in the peptide/oligonucleotide complexes decreases in the following order: short ssDNAs ≫ long ssDNAs > DNA hairpins > dsDNAs. The level of probe hydration was related to the ability of the probe to stack with the DNA bases or base pairs in the various complexes. The changes in the N*/T* ratio upon interaction of the labeled Tat peptides with DNA were about 3-fold larger with the new probes as compared to the parent FC label, in line with the higher sensitivity of the new probes to the environment. One of these labels, presenting the most compact geometry, showed the highest sensitivity, probably due to its optimal stacking with the DNA bases. Thus, the new hydration-sensitive labels appear as improved highly sensitive tools to site-selectively monitor the binding of peptides to oligonucleotides and nucleic acids.

Relaxation dynamics of nucleosomal DNA

Recent experimental and theoretical evidence demonstrates that proteins and water in the hydration layer can follow complex stretched exponential or power law relaxation dynamics. Here, we report on a 50 ns all atom molecular dynamics (MD) simulation of the yeast nucleosome, where the interactions between DNA, histones, surrounding water and ions are explicitly included. DNA interacts with the histone core in 14 locations, approximately every 10.4 base pairs. We demonstrate that all sites of interaction exhibit anomalously slow power law relaxation, extending up to 10 ns, while fast exponential relaxation dynamics of hundreds of picoseconds applies to DNA regions outside these locations. The appearance of 1/f(alpha) noise or pink noise in DNA dynamics is ubiquitous. For histone-bound nucleotide dynamics alpha --> 1 and is a signature of complexity of the protein-DNA interactions. For control purposes two additional DNA simulations free of protein are conducted. Both utilize the same sequence of DNA, as found the in the nucleosome. In one simulation the initial conformation of the double helix is a straight B-form. In the other, the initial conformation is super helical. Neither of these simulations exhibits the variation of alpha as a function of position, the measure of power law for dynamical behavior, which we observe in the nucleosome simulation. The unique correspondence (high alpha to DNA-histone interaction sites, low alpha to free DNA sites), suggests that alpha may be an important and new quantification of protein-DNA interactions for future experiments.

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.

Sensing peptide-oligonucleotide interactions by a two-color fluorescence label: application to the HIV-1 nucleocapsid protein

We present a new methodology for site-specific sensing of peptide–oligonucleotide (ODN) interactions using a solvatochromic fluorescent label based on 3-hydroxychromone (3HC). This label was covalently attached to the N-terminus of a peptide corresponding to the zinc finger domain of the HIV-1 nucleocapsid protein (NC). On interaction with target ODNs, the labeled peptide shows strong changes in the ratio of its two emission bands, indicating an enhanced screening of the 3HC fluorophore from the bulk water by the ODN bases. Remarkably, this two-color response depends on the ODN sequence and correlates with the 3D structure of the corresponding complexes, suggesting that the 3HC label monitors the peptide–ODN interactions site-specifically. By measuring the two-color ratio, we were also able to determine the peptide–ODN-binding parameters and distinguish multiple binding sites in ODNs, which is rather difficult using other fluorescence methods. Moreover, this method was found to be more sensitive than the commonly used steady-state fluorescence anisotropy, especially in the case of small ODNs. The described methodology could become a new universal tool for investigating peptide–ODN interactions.

Role of Monovalent Counterions in the Ultrafast Dynamics of DNA

This paper examines the contribution of counterion motion to the electric-field dynamics in the interior of DNA. The electric field is measured by a coumarin fluorophore that is synthetically incorporated into an oligonucleotide, where it replaces a native base pair. The DNA is a 17-base-pair oligomer with no A- or G-tracts. Time-resolved Stokes-shift measurements on the coumarin are made from 40 ps to 40 ns with each of the alkali ions and or one of several tetraalkylammonium ions as the DNA counterion. With the possible exception of rubidium, there are no indications of site-specific binding of the counterions. For sodium and other ions with a smaller hydrodynamic radius, the dynamics are identical and are fit to a power law. For larger ions, there is a progressive increase in the rate of shifting after 1 ns. This effect correlates with the hydrodynamic radius of the counterion. The lack of change in the spectral shape of the emission shows that neither the broadly distributed power-law relaxation nor the extra nanosecond dynamics are due to heterogeneity in the relaxation rates of different helices.

Ultrafast dynamics in DNA: "fraying" at the end of the helix

The dynamics of the electric fields in the interior of DNA are measured by using oligonucleotides in which a native base pair is replaced by a dye molecule (coumarin 102) whose emission spectrum is sensitive to the local electric field. Time-resolved measurements of the emission spectrum have been extended to a six decade time range (40 fs to 40 ns) by combining results from time-correlated photon counting, fluorescence up-conversion, and transient absorption. Recent results showed that when the reporter is placed in the center of the oligonucleotide, the dynamics are very broadly distributed over this entire time range and do not show specific time constants associated with individual processes (Andreatta, D.; et al. J. Am. Chem. Soc. 2005, 127, 7270). This paper examines an oligonucleotide with the reporter near its end. The broadly distributed relaxation seen before remains with little attenuation. In addition, a new relaxation with a well-defined relaxation time of 5 ps appears. This process is assigned to the rapid component of "fraying" at the end of the helix.

A new class of laser dyes, 2-oxa-bicyclo[3.3.0]octa-4,8-diene-3,6-diones, with unity fluorescence yield

A new class of highly fluorescent dyes, 4,8-diphenyl-2-oxa-bicyclo[3.3.0]octa-4,8-diene-3,6-diones (1a-c), have been synthesized, they all exhibit unity fluorescence quantum yield and short radiative lifetime (< 4 ns) in common organic solvents and have demonstrated remarkable amplified spontaneous emission with a gain efficiency of > 10.