Complex local dynamics in DNA on the picosecond and nanosecond time scales

Interfacial Layer Breaker: A Violation of Stokes' Law in High-Speed Atomic Force Microscope Flows

Structured water near surfaces is important in nonclassical crystallization, biomineralization, and restructuring of cellular membranes. In addition to equilibrium structures, studied by atomic force microscopy (AFM), high-speed AFM (H-S AFM) can now detect piconewton forces in microseconds. With increasing speeds and decreasing tip diameters, there is a danger that continuum water models will not hold, and molecular dynamic (MD) simulations would be needed for accurate predictions. MD simulations, however, can only evolve over tens of nanoseconds due to memory and computational efficiency/speed limitations, so new methods are needed to bridge the gap. Here, we report a hybrid, multiscale simulation method, which can bridge the size and time scale gaps to existing experiments. Structured water is studied between a moving silica AFM colloidal tip and a cleaved mica surface. The computational domain includes 1,472,766 atoms. To mimic the effect of long-range hydrodynamic forces occurring in water, when moving the AFM tip at speeds from 5 × 10^-7 to 30 m/s, a hybrid multiscale method with local atomistic resolution is used, which serves as an effective open-domain boundary condition. The multiscale simulation is thus equivalent to using a macroscopically large computational domain with equilibrium boundary conditions. Quantification of the drag force shows the breaking of continuum behavior. Nonmonotonic dependence on both the tip speed and distance from the surface implies breaking of the hydration layer around the moving tip at time scales smaller than water cluster formation and strong water compressibility effects at the highest speeds.

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.

Weak nanoscale chaos and anomalous relaxation in DNA

Anomalous nonexponential relaxation in hydrated biomolecules is commonly attributed to the complexity of the free-energy landscapes, similarly to polymers and glasses. It was found recently that the hydrogen-bond breathing of terminal DNA base pairs exhibits a slow power-law relaxation attributable to weak Hamiltonian chaos, with parameters similar to experimental data. Here, the relationship is studied between this motion and spectroscopic signals measured in DNA with a small molecular photoprobe inserted into the base-pair stack. To this end, the earlier computational approach in combination with an analytical theory is applied to the experimental DNA fragment. It is found that the intensity of breathing dynamics is strongly increased in the internal base pairs that flank the photoprobe, with anomalous relaxation quantitatively close to that in terminal base pairs. A physical mechanism is proposed to explain the coupling between the relaxation of base-pair breathing and the experimental response signal. It is concluded that the algebraic relaxation observed experimentally is very likely a manifestation of weakly chaotic dynamics of hydrogen-bond breathing in the base pairs stacked to the photoprobe and that the weak nanoscale chaos can represent an ubiquitous hidden source of nonexponential relaxation in ultrafast spectroscopy.

Reaction Path Averaging: Characterizing the Structural Response of the DNA Double Helix to Electron Transfer

A polarizable environment, prominently the solvent, responds to electronic changes in biomolecules rapidly. The knowledge of conformational relaxation of the biomolecule itself, however, may be scarce or missing. In this work, we describe in detail the structural changes in DNA undergoing electron transfer between two adjacent nucleobases. We employ an approach based on averaging of tens to hundreds of thousands of nonequilibrium trajectories generated with molecular dynamics simulation, and a reduction of dimensionality suitable for DNA. We show that the conformational response of the DNA proceeds along a single collective coordinate that represents the relative orientation of two consecutive base pairs, namely, a combination of helical parameters shift and tilt. The structure of DNA relaxes on time scales reaching nanoseconds, contributing marginally to the relaxation of energies, which is dominated by the modes of motion of the aqueous solvent. The concept of reaction path averaging (RPA), conveniently exploited in this context, makes it possible to filter out any undesirable noise from the nonequilibrium data, and is applicable to any chemical process in general.

Observation of coherent delocalized phonon-like modes in DNA under physiological conditions

Underdamped terahertz-frequency delocalized phonon-like modes have long been suggested to play a role in the biological function of DNA. Such phonon modes involve the collective motion of many atoms and are prerequisite to understanding the molecular nature of macroscopic conformational changes and related biochemical phenomena. Initial predictions were based on simple theoretical models of DNA. However, such models do not take into account strong interactions with the surrounding water, which is likely to cause phonon modes to be heavily damped and localized. Here we apply state-of-the-art femtosecond optical Kerr effect spectroscopy, which is currently the only technique capable of taking low-frequency (GHz to THz) vibrational spectra in solution. We are able to demonstrate that phonon modes involving the hydrogen bond network between the strands exist in DNA at physiologically relevant conditions. In addition, the dynamics of the solvating water molecules is slowed down by about a factor of 20 compared with the bulk.

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.

Solvation dynamics and rotational relaxation of coumarin 153 in mixed micelles of Triton X-100 and cationic gemini surfactants: effect of composition and spacer chain length of gemini surfactants

To demonstrate simultaneously how the solvation dynamics and rotational relaxation in nonionic micelles change with the composition of a gemini surfactant and how this change depends on spacer chain length of gemini surfactants.

Modeling DNA Dynamics under Steady Deforming Forces and Torques

An algorithm is developed for modeling atom-level dynamics of DNA subjected to steady external torques. For completeness, simulations with steady stretching loads are also considered. The algorithms were tested in Brownian dynamics simulations of discrete wormlike chain models with calibrated elastic properties to confirm that the elastic responses induced are of desired type and magnitude and that no side effects appear. The same methods were next used in a series of 100-ns all-atom MD simulations of tetradecamer DNA fragments with explicit water and counterions. The results demonstrate the possibility of probing regular elastic responses in DNA under low, nearly physiological amplitudes of forces and torques.

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.

Direct Observation of Coupling between Structural Fluctuation and Ultrafast Hydration Dynamics of Fluorescent Probes in Anionic Micelles

The coupling of structural fluctuation and the dynamics of associated water molecules of biological macromolecules is vital for various biological activities. Although a number of molecular dynamics (MD) studies on proteins/DNA predicted the importance of such coupling, experimental evidence of variation of hydration dynamics with controlled structural fluctuation even in model macromolecule is sparse and raised controversies in the contemporary literature. Here, we have investigated dynamics of hydration at the surfaces of two similar anionic micelles sodium dodecyl sulfate (SDS) and sodium dodecylbenzenesulfonate (SDBS) as model macromolecules using coumarin 500 (C500) as spectroscopic probe with femtosecond to picosecond time resolution up to 20 ns time window. The constituting surfactants SDS and SDBS are structurally similar except one benzene moiety in the SDBS may offer additional rigidity to the SDBS micelles through π-stacking and added bulkiness. The structural integrity of the micelles in the aqueous medium is confirmed in dynamic light scattering (DLS) studies. A variety of studies including polarization gated fluorescence spectroscopy and quasielastic neutron scattering (QENS) have been used to confirm differential structural fluctuation of SDS and SDBS micelles. We have also employed femtosecond-resolved Förster resonance energy transfer (FRET) in order to study binding of a cationic organic ligand ethidium bromide (EtBr) salt at the micellar surfaces. The distance distribution of the donor (C500)-acceptor (EtBr) in the micellar media reveals the manifestation of the structural flexibility of the micelles. Our studies on dynamical coupling of the structural flexibility with surface hydration in the nanoscopic micellar media may find the relevance in the "master-slave" type water dynamics in biologically relevant macromolecules.

Investigation of changes in the arrangement of water molecules and salt ions surrounding different atoms of the DNA molecule during the melting process: a molecular dynamics simulation study

A molecular dynamics simulation was performed on a B-DNA duplex (CGCGAATTGCGC) at different temperatures. The DNA was immerged in a saltwater medium with 1 mol/L NaCl concentration. The arrangements of water molecules and cations around the different atoms of DNA on the melting pathway were investigated. Almost for all atoms of the DNA by double helix → single-stranded transition, the water molecules released from the DNA duplex and cations were close to single-stranded DNA, but this behavior was not clearly seen at melting temperatures. Therefore, release of water molecules and cations approaching the DNA by the increase of temperature does not have any effect on the sharpness of the transition curve. Most of the water molecules and cations were found to be around the negatively charged phosphate oxygen atoms. The number of water molecules released from the first shell hydration upon melting in the minor groove was higher than in the major groove, and intrusion of cations into the minor groove after melting was higher than into the major groove. The hydrations of imino protons were different from each other and were dependent on DNA bases.

On the absence of intrahelical DNA dynamics on the μs to ms timescale

DNA helices display a rich tapestry of motion on both short (<100 ns) and long (>1 ms) timescales. However, with the exception of mismatched or damaged DNA, experimental measures indicate that motions in the 1 μs to 1 ms range are effectively absent, which is often attributed to difficulties in measuring motions in this time range. We hypothesized that these motions have not been measured because there is effectively no motion on this timescale, as this provides a means to distinguish faithful Watson-Crick base-paired DNA from damaged DNA. The absence of motion on this timescale would present a 'static' DNA sequence-specific structure that matches the encounter timescales of proteins, thereby facilitating recognition. Here we report long-timescale (~10-44 μs) molecular dynamics simulations of a B-DNA duplex structure that addresses this hypothesis using both an 'Anton' machine and large ensembles of AMBER GPU simulations.

Convergence and reproducibility in molecular dynamics simulations of the DNA duplex d(GCACGAACGAACGAACGC)

BACKGROUND

The structure and dynamics of DNA are critically related to its function. Molecular dynamics simulations augment experiment by providing detailed information about the atomic motions. However, to date the simulations have not been long enough for convergence of the dynamics and structural properties of DNA.

METHODS

Molecular dynamics simulations performed with AMBER using the ff99SB force field with the parmbsc0 modifications, including ensembles of independent simulations, were compared to long timescale molecular dynamics performed with the specialized Anton MD engine on the B-DNA structure d(GCACGAACGAACGAACGC). To assess convergence, the decay of the average RMSD values over longer and longer time intervals was evaluated in addition to assessing convergence of the dynamics via the Kullback-Leibler divergence of principal component projection histograms.

RESULTS

These molecular dynamics simulations-including one of the longest simulations of DNA published to date at ~44μs-surprisingly suggest that the structure and dynamics of the DNA helix, neglecting the terminal base pairs, are essentially fully converged on the ~1-5μs timescale.

CONCLUSIONS

We can now reproducibly converge the structure and dynamics of B-DNA helices, omitting the terminal base pairs, on the μs time scale with both the AMBER and CHARMM C36 nucleic acid force fields. Results from independent ensembles of simulations starting from different initial conditions, when aggregated, match the results from long timescale simulations on the specialized Anton MD engine.

GENERAL SIGNIFICANCE

With access to large-scale GPU resources or the specialized MD engine "Anton" it is possible for a variety of molecular systems to reproducibly and reliably converge the conformational ensemble of sampled structures. This article is part of a Special Issue entitled: Recent developments of molecular dynamics.

Microscopic origin of the logarithmic time evolution of aging processes in complex systems

There exists compelling experimental evidence in numerous systems for logarithmically slow time evolution, yet its full theoretical understanding remains elusive. We here introduce and study a generic transition process in complex systems, based on nonrenewal, aging waiting times. Each state n of the system follows a local clock initiated at t = 0. The random time τ between clock ticks follows the waiting time density ψ(τ). Transitions between states occur only at local clock ticks and are hence triggered by the local forward waiting time, rather than by ψ(τ). For power-law forms ψ(τ) ≃ τ(-1-α) (0 < α < 1) we obtain a logarithmic time evolution of the state number ⟨n(t)⟩ ≃ log(t/t(0)), while for α > 2 the process becomes normal in the sense that ⟨n(t)⟩ ≃ t. In the intermediate range 1 < α < 2 we find the power-law growth ⟨n(t)⟩ ≃ t(α-1). Our model provides a universal description for transition dynamics between aging and nonaging states.

Hidden in plain sight: subtle effects of the 8-oxoguanine lesion on the structure, dynamics, and thermodynamics of a 15-base pair oligodeoxynucleotide duplex

The base lesion 8-oxoguanine is formed readily by oxidation of DNA, potentially leading to G → T transversion mutations. Despite the apparent similarity of 8-oxoguanine-cytosine base pairs to normal guanine-cytosine base pairs, cellular base excision repair systems effectively recognize the lesion base. Here we apply several techniques to examine a single 8-oxoguanine lesion at the center of a nonpalindromic 15-mer duplex oligonucleotide in an effort to determine what, if anything, distinguishes an 8-oxoguanine-cytosine (8oxoG-C) base pair from a normal base pair. The lesion duplex is globally almost indistinguishable from the unmodified parent duplex using circular dichroism spectroscopy and ultraviolet melting thermodynamics. The DNA mismatch-detecting photocleavage agent Rh(bpy)(2)chrysi(3+) cleaves only weakly and nonspecifically, revealing that the 8oxoG-C pair is locally stable at the level of the individual base pairs. Nuclear magnetic resonance spectra are also consistent with a well-conserved B-form duplex structure. In the two-dimensional nuclear Overhauser effect spectra, base-sugar and imino-imino cross-peaks are strikingly similar between parent and lesion duplexes. Changes in chemical shift due to the 8oxoG lesion are localized to its complementary cytosine and to the 2-3 bp immediately flanking the lesion on the lesion strand. Residues further removed from the lesion are shown to be unperturbed by its presence. Notably, imino exchange experiments indicate that the 8-oxoguanine-cytosine pair is strong and stable, with an apparent equilibrium constant for opening equal to that of other internal guanine-cytosine base pairs, on the order of 10(-6). This collection of experiments shows that the 8-oxoguanine-cytosine base pair is incredibly stable and similar to the native pair.

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.

Logarithmic based optical delay for time-resolved data collection

A method has been established that generates values spaced according to a mathematical function, specifically the logarithm function that can be applied to a stepper motor. Here, it is applied to yield logarithmically spaced time delay points for subnanosecond interferometric time-resolved experiments using a stepper motor controlled translation stage. Application of this method is discussed in terms of three input parameters: the optical delay stage time resolution, dt; the time of maximum delay, d(stop); and the desired number of data points, N. The method improves the efficiency of interferometric time-resolved data collection while providing data collection effective to determine decay parameters. In principle, this technique could be generalized to any mathematical function.

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.

Molecular recognition of plant DNA: does it differ from conventional animal DNA?

The recognition mechanism of DNA with small drugs/ligands is an important field of research from pharmacological point of view. Such studies are ample with DNAs extracted from animal cells, but are rare for those extracted from plant cells. However, such a study is strongly demanding for the formulation of pesticides and other agrochemicals. In this contribution, for the first time, we report the interaction of two well-known DNA binder ethidium bromide (EB) and Hoechst 33258 (H33258) with two genomic DNAs extracted from the leaves of Ricinus communis L. (castor bean) and Mangifera indica (mango) using steady-state and picosecond-resolved fluorescence spectroscopy. The purity of the extracted DNAs is confirmed from gel electrophoresis and optical absorption studies. As evidenced from the circular dichroism (CD) measurements the DNAs retain physiologically relevant B forms. The well-known DNA intercalator EB has been found to show an additional electrostatic mode of binding with the DNAs, which is not present in the conventional animal DNAs. The binding affinity of EB is found to be even weaker for the DNA extracted from M. indica compared to that in R. communis L. On the other hand, the binding affinity of H33258 with the plant DNAs is found to be comparable to that of animal DNAs. The difference in interaction could be rationalized from the possible differences in the base sequences.

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.

Conformational dynamics at the active site of alpha-chymotrypsin and enzymatic activity

The role of dynamical flexibility at the active site of a proteolytic enzyme alpha-chymotrypsin (CHT) has been correlated with its catalytic activity. The temperature-dependent efficiency of catalysis reveals a bell-shaped feature with a peak at 37 degrees C, the typical body temperature of homeothermal animals. The overall structural integrity of the enzyme in our experimental temperature range has been confirmed from dynamic light scattering (DLS) and circular dichroism (CD) studies. We have followed the dynamical evolution at the active site of CHT with temperature using picosecond-resolved fluorescence anisotropy of anthraniloyl probe (covalently attached to the serine-195 residue) and a substrate mimic (inhibitor) proflavin. The conformational dynamics at the active site is found to have a distinct connection with the enzyme functionality. The conformational flexibility of the enzyme is also evidenced from the compressibility studies on the enzyme. The site selective fluorescence detected circular dichroism (FDCD) studies reveal that the conformational flexibility of the enzyme has an effect on the structural perturbation at the active site. We have also proposed the possible implications of the dynamics in the associated energetics.

Dynamics in the DNA recognition by DAPI: exploration of the various binding modes

Two distinct modes of interaction of the fluorescent probe 4',6-diamidino-2-phenylindole (DAPI), depending on the sequence of DNA, have been reported in the literature. In the present study, the dynamics of solvation has been utilized to explore the binding interaction of DAPI to DNA oligomers of different sequences. Picosecond-resolved fluorescence and polarization-gated anisotropy have been used to characterize the binding of DAPI to the different oligomers. In the double-stranded dodecamer of sequence CGCGAATTCGCG (oligo1), the solvation relaxation dynamics of the probe reveals time constants of 0.130 ns (75%) and 2.35 ns (25%). Independent exploration of the minor-groove environment of oligo1 using another well-known minor-groove binder Hoechst 33258 (H258) shows similar timescales, further confirming minor-groove binding of DAPI to oligo1. In the double-stranded dodecamer (oligo2) having the sequence GCGCGCGCGCGC, where intercalation has been reported in the literature, no solvation is observed in our experimental window. DAPI bound to oligo2 shows quenching of fluorescence compared to that of DAPI in a buffer. The quenching of fluorescence of DAPI intercalated in DNA is also borne out by the appearance of a fast component of 30 ps in the fluorescence lifetime, revealing electron transfer to DAPI from GC base pairs, between which it intercalates. In addition to this, the excited-state lifetime of the probe in the DAPI-DNA complex also shows a time constant similar to that of the dye in a buffer, indicating that the excited-state photoprocesses associated with the free dye is also operative in this binding mode, consistent with the binding geometry of the DAPI in the DNA. The dynamics of DAPI in calf thymus DNA having a random sequence of base pairs is similar to that associated with the DNA minor groove. Our studies clearly explore the structure-dynamics correlation of the DAPI-DNA complex in the two distinct modes of interaction of DAPI with DNA.

Interaction of Hoechst 33258 and Ethidium with Histone1−DNA Condensates

We report structural and dynamical aspects of DNAs from various sources including synthetic oligonucleotides in bulk buffer and as a complex with histone1 (H1). High-resolution transmission electron microscopic (HRTEM) studies reveal the structural change of the DNAs upon complexation with H1 leading to formation of compact-globular and hollow-toroidal particles. In order to explore the functionality of ligand binding of the DNAs and their complexes with H1, we have used two biologically common fluorescent probes Hoechst 33258 (H33258) and Ethidium (EB) as model ligands. Picosecond resolved fluorescence and polarization gated anisotropy studies examined that the minor groove binding of H33258 to the genomic DNA-H1 complex remains almost unaltered. However, the intercalative interaction of EB with the DNA in the complex is severely perturbed compared to that with the DNA in bulk buffer. Time-dependent solvochromic effect of the probe H33258 further elucidates the dynamical solvation, which is reflective of the overall environmental relaxation of the DNAs upon condensation by H1. We have also performed circular dichroism (CD) studies on the DNAs and their complexes with H1, which reveal the change in conformation of the DNAs in the complexes. Our studies in the ligand-binding mechanisms of the DNA-H1 complex are important to understand the mechanism of drug binding to linker DNA in condensed chromatin.

Femtosecond time- and wavelength-resolved fluorescence and absorption spectroscopic study of the excited states of adenosine and an adenine oligomer

By employing broadband femtosecond Kerr-gated time-resolved fluorescence (KTRF) and transient absorption (TA) techniques, we report the first (to our knowledge) femtosecond combined time- and wavelength-resolved study on an ultraviolet-excited nucleoside and a single-stranded oligonucleotide (namely adenosine (Ado) and single-stranded adenine oligomer (dA)(20)) in aqueous solution. With the advantages of the ultrafast time resolution, the broad spectral and temporal probe window, and a high sensitivity, our KTRF and TA results enable the real time monitoring and spectral characterization of the excited-state relaxation processes of the Ado nucleoside and (dA)(20) oligonucleotide investigated. The temporal evolution of the 267 nm excited Ado KTRF spectra indicates there are two emitting components with lifetimes of approximately 0.13 ps and approximately 0.45 ps associated with the L(a) and L(b) pipi* excited states, respectively. These Ado results reveal no obvious evidence for the involvement of the npi* state along the irradiative internal conversion pathway. A distinct mechanism involving only the two pipi* states has been proposed for the ultrafast Ado deactivation dynamics in aqueous solution. The time dependence of the 267 nm excited (dA)(20) KTRF and TA spectra reveals temporal evolution from an ultrafast "A-like" state (with a approximately 0.39 ps decay time) to a relatively long-lived E(1) "excimer" (approximately 4.3 ps decay time) and an E(2) "excimer-like" (approximately 182 ps decay time) state. The "A-like" state has a spectral character closely resembling the excited state of Ado. Comparison of the spectral evolution between the results for Ado and (dA)(20) provides unequivocal evidence for the local excitation character of the initially photoexcited (dA)(20). The rapid transformation of the locally excited (dA)(20) component into the delocalized E(1) "excimer" state which then further evolves into the E(2) "excimer-like" state indicates that base stacking has a high ability to modify the excited-state deactivation pathway. This modification appears to occur by suppressing the internal conversion pathway of an individually excited base component where the stacking interaction mediates efficient interbase energy transfer and promotes formation of the collective excited states. This feature of the local excitation that is subsequently followed by rapid energy delocalization into nearby bases may occur in many base multimer systems. Our results provide an important new contribution to better understanding DNA photophysics.

Direct Observation of Essential DNA Dynamics: Melting and Reformation of the DNA Minor Groove

The dynamics of bound water and ions present in the minor groove of a dodecamer DNA has been decoupled from that of the long-range twisting/bending of the DNA backbone, using the minor groove binder Hoechst 33258 as a fluorescence reporter in the picosecond-resolved time window. The bound water and ions are essential structural components of the minor groove and are destroyed with the destruction of the minor groove when the dodecamer melts at high temperatures and reforms on subsequent cooling of the melted DNA. The melting and rehybridization of the DNA has been monitored by the changes in secondary structure using circular dichroism (CD) spectroscopy. The change in the relaxation dynamics of the DNA has been studied with picosecond resolution at different temperatures, following the temperature-dependent melting and rehybridization profile of the dodecamer, using time-resolved emission spectra (TRES). At room temperature, the relaxation dynamics of DNA is governed by a 40 ps (30%) and a 12.3 ns (70%) component. The dynamics of bound water and ions present in the minor groove is characterized by the 40 ps component in the relaxation dynamics of the probe bound in the minor groove of the dodecamer DNA. Analyses of the TRES taken at different temperatures show that the contribution of this component decreases and ultimately vanishes with the destruction of the minor groove and reappears again with the reformation of the groove. The dynamical behavior of bound water molecules and ions of a genomic DNA (from salmon testes) at different temperatures is also found to be consistent with that of the dodecamer. The longer component of approximately 10 ns in the DNA dynamics is found to be associated with the long-range bending/twisting of the DNA backbone and the associated counterions. The transition from bound water to free water at the DNA surface, indicative of the change in the hydration number associated with each base pair, has also been ascertained in the case of the genomic DNA at different temperatures by employing densimetric and acoustic techniques.

Ultrafast energy transfer and structural dynamics in DNA

Ultrafast structural dynamics concomitant to excitation energy transfer in DNA has been studied using a pair of pyrene-labeled DNA bases. The temporal evolution of the femtosecond pump-probe spectra reveals the existence of two electronic coupling pathways, through-base stack and through-space, which lead to excitation energy transfer and excimer formation even when the labeled DNA bases are separated by one AT base pair. The electronic coupling which mediates through-base stack energy transfer is so strong that a new absorption band arises in the excited-state absorption spectrum within 300 fs. From the analysis of time-dependent spectral shifts due to through-space excimer formation, the local structural dynamics and flexibility of DNA are characterized on the picosecond and nanosecond time scale.

Furan Decorated Nucleoside Analogues as Fluorescent Probes: synthesis, photophysical evaluation and site-specific incorporation

The synthesis and photophysical evaluation of modified nucleoside analogues in which a five-membered heterocycle (furan, thiophene, oxazole and thiazole) is attached to the 5 position of 2'-deoxyuridine are reported. The furan containing derivative is identified as the most promising responsive nucleoside of this family due to its emission quantum efficiency and degree of sensitivity to its microenvironment. The furan moiety was then attached to the 5 position of 2'-deoxycytidine as well as the 8 position of adenosine and guanosine. Photophysical evaluation of these four furan containing nucleoside analogues reveal distinct differences in the absorption, emission and quantum efficiency depending upon the class of nucleoside (pyrimidine or purine). Comparing the photophysical properties of all furan containing nucleosides, identifies the furan thymidine analogue, 5-(fur-2-yl)-2'-deoxyuridine, as the best candidate for use as a responsive fluorescent probe in nucleic acids. 5-(fur-2-yl)-2'-deoxyuridine was then converted to the corresponding phosphoramidite and site specifically incorporated into DNA oligonucleotides with greater than 88% coupling efficiency. Such furan-modified oligonucleotides form stable duplexes upon hybridization to their complementary DNA strands and display favorable fluorescent features.