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.     

Solvent effects on the conformation of DNA dodecamer segment: a simulation study

Different solvent temperatures with five kinds of counterions are used to investigate solvent effects on the DNA microscopic structure. The dodecamer d (CGCGAATTCGCG) DNA segment is merged into the solvents and its conformation transition is studied with the molecular dynamics simulations in detail. For the simple point charge model of water molecule with Na(+) counterions, as temperature increases from 200 K to 343 K, the duplex DNA changes from stiff B form to a state between A form and B form, which we define as mixed (A-B) structure, with a double helix unwinding. To study the counterions effects, other four alkali cations, Li(+), K(+), Rb(+), or Cs(+) ions, are substituted for Na(+) ions at 298 K and 343 K, respectively. For the cases of Li(+), Rb(+), and Cs(+) ions, the duplex DNA becomes more flexible with sugar configuration changing form C2'-endo to C1'-endo type and the width and depth of minor groove at CpG and GpC steps moving towards A values, as the mass of the counterions decreasing. For the case of K(+) ions, DNA-K(+) interaction widens the width of minor and major grooves at ApA steps and TpT steps, respectively. It seems that the light ions (Li(+) or Na(+)) prefer to interact with the free phosphate oxygen atoms while the heavier ions (Rb(+) and Cs(+)) strongly interact with the base pairs.     

DNA dynamics in a water drop without counterions

Because of its polyionic character, the DNA double helix is stable and biologically active only in salty aqueous media where its charge is compensated by solvent counterions. Monovalent metal ions are ubiquitous in DNA environment, and they are usually considered as the possible driving force of sequence-dependent modulations of DNA structure that make it recognizable by proteins. In an effort to directly examine this hypothesis, MD simulations of DNA in a water drop surrounded by vacuum were carried out, which relieves the requirement of charge neutrality. Surprisingly, with zero concentration of counterions, a dodecamer DNA duplex appears metastable, and its structure remains similar to that observed in experiment, including the minor groove narrowing in the dodecamer d(CGCGAATTCGCG)(2) often considered as the most evident cation effect. It is suggested that the same computational approach will allow one to simulate dynamics of long DNA chains more efficiently than with periodical boundary conditions.     

The dynamics of water at DNA interfaces: computational studies of Hoechst 33258 bound to DNA

Together, spectroscopy combined with computational studies that relate directly to the experimental measurements have the potential to provide unprecedented insight into the dynamics of important biological processes. Recent time-resolved fluorescence experiments have shown that the time scales for collective reorganization at the interface of proteins and DNA with water are more than an order of magnitude slower than in bulk aqueous solution. The molecular interpretation of this change in the collective response is somewhat controversial some attribute the slower reorganization to dramatically retarded water motion, while others describe rapid water dynamics combined with a slower biomolecular response. To connect directly to solvation dynamics experiments of the fluorescent probe Hoechst 33258 (H33258) bound to DNA, we have generated 770 ns of molecular dynamics (MD) simulations and calculated the equilibrium and nonequilibrium solvation response to excitation of the probe. The calculated time scales for the solvation response of H33258 free in solution (0.17 and 1.4 ps) and bound to DNA (1.5 and 20 ps) are highly consistent with experiment (0.2 and 1.2 ps, 1.4 and 19 ps, respectively). Decomposition of the calculated response revealed that water solvating the probe bound to DNA was still relatively mobile, only slowing by a factor of 2-3, while DNA motion was responsible for the long-time component (approximately 20 ps).     

Influence of the dynamic positions of cations on the structure of the DNA minor groove: sequence-dependent effects

Different models for minor groove structures predict that the conformation is essentially fixed by sequence and has an influence on local ion distribution or alternatively that temporal positions of ions around the minor groove can affect the structure if they neutralize cross-strand phosphate charges. Our previous studies show that the minor groove in an AATT dodecamer responds to local sodium ion positions and is narrow when ions neutralize cross-strand phosphate-phosphate charges [J. Am. Chem. Soc. 2000, 122, 10513-10520]. Previous results from a number of laboratories have shown that G-tracts often have a wider minor groove than A-tracts, but they do not indicate whether this is due to reduced flexibility or differences in ion interactions. We have undertaken a molecular dynamics study of a d(TATAGGCCTATA) duplex to answer this question. The results show that the G-tract has the same amplitude of minor groove fluctuations as the A-tract sequence but that it has fewer ion interactions that neutralize cross-strand phosphate charges. These results demonstrate that differences in time-average groove width between A- and G-tracts are due to differences in ion interactions at the minor groove. When ions neutralize the cross-strand phosphates, the minor groove is narrow. When there are no neutralizing ion interactions, the minor groove is wide. The population of structures with no ion interactions is larger with the GGCC than with the AATT duplex, and GGCC has a wider time-average minor groove in agreement with experiment.     

Nature of minor-groove binders-DNA complexes in the gas phase

The structure of noncovalent complexes of DNA duplex with minor groove binders (mG-binders) has been analyzed by state of the art molecular dynamics (MD) simulations. More than 3.3 micros of MD trajectories (including 4 x 0.5 micros trajectories) were collected for the Dickerson's dodecamer bound to DAPI, Hoechst 33258, and Netropsin. Comparison of these trajectories with control simulations in water allowed us to determine that the extreme dehydration and partial neutralization occurring during electrospray experiments does not produce the disruption of the DNA:mG-binder complexes or the dissociation of the two strands of the duplex. Irrespective of the drug and the simulation conditions the mG-binders remains bound to the DNA near the preferential binding position in aqueous conditions. Large distortions appear in the two DNA strands, which maintain however a memory of the original DNA duplex structure in water, and a general helical-like conformation.     

Atomic Resolution Insight into Sac7d Protein Binding to DNA and Associated Global Changes by Molecular Dynamics Simulations

Sac7d is a small, thermostable protein that induces large helical deformations in DNA upon association. Starting from multiple initial placements of the unbound Sac7d structure relative to a B‐DNA oligonucleotide, molecular dynamics (MD) simulations were employed to directly follow several successful binding events at atomic resolution that resulted in structures in close agreement with the native complex geometry. The final native complex formed rapidly within tenths of nanoseconds and included simultaneous large‐scale kinking, groove opening, twisting, and intercalation in the target DNA. The simulations indicate that the complex formation process involves initial non‐native contacts that helped in reaching the final bound state, with residues intercalated at the center of the kinked DNA. It was also possible to identify several long‐lived trapped intermediate states of the binding process and to follow sliding processes of Sac7d along the DNA minor groove.     

Multiple time scales in solvation dynamics of DNA in aqueous solution: the role of water, counterions, and cross-correlations

Recent time domain experiments have explored solvation dynamics of a probe located inside a DNA duplex, in an effort to gain information, e.g., on the dynamics of water molecules in the DNA major and minor grooves and their environment. Multiple time constants in the range of a few picoseconds to several nanoseconds were obtained. We have carried out 15 ns long atomistic molecular dynamics simulations to study the solvation dynamics of bases of a 38 base-pair long DNA duplex in an aqueous solution containing counterions. We have computed the energy-energy time correlation function (TCF) of the four individual bases (A, T, G, and C) to characterize the solvation dynamics. All the TCFs display highly nonexponential decay with time. When the trajectories are analyzed with 100 fs time resolution, the TCF of each base shows initial ultrafast decay (with tau1 approximately equal 60-80 fs) followed by two intermediate components (tau2 approximately equal 1 ps, tau3 approximately equal 20-30 ps), in near complete agreement with a recent time domain experiment on DNA solvation. Interestingly, the solvation dynamics of each of the four different nucleotide bases exhibit rather similar time scales. To explore the existence of slow relaxation at longer times reported recently in a series of experiments, we also analyzed the solvation TCFs calculated with longer time trajectories and with a larger time resolution of 1 ps. In this case, an additional slow component with a time constant of the order of 250 ps is observed. Through an analysis of partial solvation TCFs, we find that the slow decay originates mainly from the interaction of the nucleotides with the dipolar water molecules and the counterions. An interesting negative cross-correlation between water and counterions is observed, which makes an important contribution to relaxation at intermediate to longer times.     

Theoretical study of the interaction between a high-valent manganese porphyrin oxyl-(hydroxo)-Mn(IV)-TMPyP and double-stranded DNA

Cationic porphyrin derivatives such as meso-tetrakis(4-N-methylpyridinium)porphyrin, TMPyP, have been shown to interact with double-stranded DNA. The manganese derivative, Mn(III)-TMPyP, activated by an oxygen donor like potassium monopersulfate, provides an efficient DNA-cleaving system. Previous experimental work1 has shown that DNA cleavage by the Mn(III)-TMPyP/KHSO(5) system was due to an oxidative attack, within the minor groove of B-DNA, at the C5' or C1' carbons of deoxyribose units. The aim of this study was to use molecular modeling to elucidate the specificity of the interactions between the transient active species oxyl-Mn(IV)-TMPyP and the DNA target. Geometric parameters, charges, and force field constants consistent with the AMBER 98 force field were calculated by DFT methods. Molecular modeling (mechanics and dynamic simulations) were performed for oxyl-(hydroxo)-Mn(IV)-TMPyP bound in the minor groove of the dodecamer d(5'-TCGTCAAACCGC)-d(5'-GCGGTTTGACGA). Geometry, interactions, and binding energy of the metalloporphyrin located at the A.T triplet region of the dodecamer were analyzed. These studies show no significant structural change of the DNA structure upon ligand binding. Mobility of the metalloporphyrin in the minor groove was restrained by the formation of a hydrogen bond between the hydroxo ligand trans to the metal-oxyl and a DNA phosphate, restricting the access of the oxyl group to the (pro-S) H atom at C5'.     

Exploring DNA groove water dynamics through hydrogen bond lifetime and orientational relaxation

Dynamics of water molecules in the grooves of DNA are of great interest both for practical (functionality of DNA) and fundamental (as examples of confined systems) interest. Here the authors employ atomistic molecular dynamics simulations to understand varying water dynamics at the minor and the major grooves of a 38 base-pair long DNA duplex in water. In order to understand and quantify the diversity in the nature of hydrogen bond due to many hydrogen bond donors and acceptors present in the four bases, they have undertaken study of hydrogen bond lifetime (HBLT) correlation functions of all the specific hydrogen bonds between the base atoms and water molecules. They find that the HBLT correlation functions are in general multiexponential, with the average lifetime depending significantly on the specificity and may thus be biologically relevant. The average hydrogen bond lifetime is longer in the minor groove than that in the major groove by almost a factor of 2. Analysis further shows that water hydrogen bonds with phosphate oxygen have substantially shorter lifetimes than those with the groove atoms. They also compute two different orientational time correlation functions (OTCFs) of the water molecules present at the major and the minor grooves and attempt to correlate OTCF with HBLT correlation function. The OTCFs in the minor groove exhibit three time scales, with the time constant of the slowest component one to two orders of magnitude longer than what is observed for bulk water. A slow component is also present for the major groove water but with shorter time constant. Interestingly, correlation between reformations allowed HBLT correlation function [C(HB)(t)] and the OTCF markedly deviates from each other in the grooves, indicating enhanced rigidity of water molecules in the grooves.     

31P NMR investigation of backbone dynamics in DNA binding sites

The backbone conformation of DNA plays an important role in the indirect readout mechanisms for protein--DNA recognition events. Thus, investigating the backbone dynamics of each step in DNA binding sequences provides useful information necessary for the characterization of these interactions. Here, we use 31P dynamic NMR to characterize the backbone conformation and dynamics in the Dickerson dodecamer, a sequence containing the EcoRI binding site, and confirm solid-state 2H NMR results showing that the C3pG4 and C9pG10 steps experience unique dynamics and that these dynamics are quenched upon cytosine methylation. In addition, we show that cytosine methylation affects the conformation and dynamics of neighboring nucleotide steps, but this effect is localized to only near neighbors and base-pairing partners. Last, we have been able to characterize the percent BII in each backbone step and illustrate that the C3pG4 and C9pG10 favor the noncanonical BII conformation, even at low temperatures. Our results demonstrate that 31P dynamic NMR provides a robust and efficient method for characterizing the backbone dynamics in DNA. This allows simple, rapid determination of sequence-dependent dynamical information, providing a useful method for studying trends in protein-DNA recognition events.     

Structural and dynamic variations in DNA hexamers containing T-T and F-F single and tandem internal mispairs

Molecular dynamics simulations of double-helical DNA oligomers have been performed to investigate differences in the structure, dynamics, and hydration of F-F and T-T mispairs. Hexamers containing F-F pairs were found to be more dynamic, especially in the region of the mispair itself. This dynamic variability derives from greater flexibility of F-F pairs. The T-T mispairs, on the other hand, were found to be comparatively tightly bound as wobble pairs. The major and minor groove edges of the T-T pairs were observed to be solvated at exposed carbonyl positions by at least one water molecule, while F-F pairs lacked solvating waters. Stacking interactions were nearly identical for T-T and F-F pairs, leading to similar average structures, even though F stacking was more dynamically variable. Solvation differences between F-F and T-T therefore support the steric exclusion model for nucleotide incorporation in DNA replication. Large differences in the orientation of minor groove functional groups, in addition to differences in solvation, further rationalize why F bases present during DNA extension events induce stalls. Two novel nucleotides are proposed to further elucidate minor groove interactions of DNA with polymerase molecules.Electronic Supplementary Material This Material consists of equilibration protocol, plots of center-of-mass stacking, water radial distribution functions, helical parameter dynamics, and dynamics data for a control AT sequence. Supplementary material is available in the online version of this article at Contribution to the Jacopo Tomasi Honorary Issue     

Competitive Na(+) and Rb(+) binding in the minor groove of DNA

Sequence-dependent coordination of alkali ions to the nucleotide bases in the minor groove of AT-tract B-DNA has recently been inferred from X-ray crystallography, solution NMR and computer simulations. Here, we present new (23)Na and (87)Rb magnetic relaxation dispersion (MRD) data that demonstrate competitive and long-lived binding of Na(+) and Rb(+) ions in the minor groove of the B-DNA duplex [d(CGCGAATTCGCG)](2). The Na(+)/Rb(+) selectivity of the minor groove is found to be weak, consistent with local structural flexibility. The ion occupancies derived from the MRD data are substantially higher than previously reported, suggesting that groove-bound ions significantly influence the energetics and structural polymorphism of DNA in vivo. For example, in the presence of 0.20 M Na(+) and 0.56 M Rb(+) at 4 degrees C, the ApT site in the minor groove is occupied by a Rb(+) ion, a Na(+) ion, or a water molecule 40, 10, and 50% of the time, respectively. In the absence of Rb(+), the Na(+) occupancy increases to 50%. At 4 degrees C, the mean residence time of groove-bound ions is 0.2 +/- 0.1 micros for Rb(+) and 10 ns to 100 micros for Na(+). A shorter correlation time of 2 ns is attributed to counterions bridging cross-strand phosphate groups.     

Molecular dynamics of minimal B‐DNA

Stable and accurate molecular dynamics (MD) of B‐DNA duplexes can be obtained in inexpensive computational conditions where only the minor groove is filled with water while the bulk solvent is represented implicitly. This model system presents significant theoretical as well as practical interest because, due to its simplicity and exceptional computational performance, it can be employed in simulations of very long DNA fragments. To better understand its properties and clarify the physical background of the effects produced by the limited water shell, dynamics of several different DNA oligomers was studied. It is found that optimal simulation conditions are reached when the explicit water is confined within the minor groove while the major groove is cleaned periodically. The internal solvent mobility appears high enough to observe in the nanosecond time scale spontaneous formation of sequence‐specific hydration patterns known from experiments. It is shown that the model produces stable MD trajectories close to the B‐DNA form regardless of the base pair sequence and that, on the other hand, the dynamics are strongly sequence dependent. Independent observations suggest that B‐DNA with only minor groove hydrated resembles its natural thermodynamic state at low water concentration; therefore, this model system can be tentatively called “minimal B‐DNA.” © 2001 John Wiley & Sons, Inc. J Comput Chem 22: 457–467, 2001     

Anomalous dielectric relaxation of water molecules at the surface of an aqueous micelle

Dielectric relaxation of aqueous solutions of micelles, proteins, and many complex systems shows an anomalous dispersion at frequencies intermediate between those corresponding to the rotational motion of bulk water and that of the organized assembly or macromolecule. The precise origin of this anomalous dispersion is not well-understood. In this work we employ large scale atomistic molecular dynamics simulations to investigate the dielectric relaxation (DR) of water molecules in an aqueous micellar solution of cesium pentadecafluorooctanoate. The simulations clearly show the presence of a slow component in the moment-moment time correlation function [PhiMW(t)] of water molecules, with a time constant of about 40 ps, in contrast to only 9 ps for bulk water. Interestingly, the orientational time correlation function [Cmu(t)] of individual water molecules at the surface exhibits a component with a time constant of about 19 ps. We show that these two time constants can be related by the well-known micro-macrorelations of statistical mechanics. In addition, the reorientation of surface water molecules exhibits a very slow component that decays with a time constant of about 500 ps. An analysis of hydrogen bond lifetime and of the rotational relaxation in the coordinate frame fixed on the micellar body seems to suggest that the 500 ps component owes its origin to the existence of an extended hydrogen bond network of water molecules at the surface. However, this ultraslow component is not found in the total moment-moment time correlation function of water molecules in the solution. The slow DR of hydration water is found to be well correlated with the slow solvation dynamics of cesium ions at the water-micelle interface.     

Power-law solvation dynamics in DNA over six decades in time

Time-resolved Stokes shifts in a dye-containing oligonucleotide have been observed over the entire time range from 40 fs to 40 ns. The dynamics fit to a power law with a small exponent of 0.15. Similar relaxation has been seen in proteins but has not been anticipated in DNA. Distinct relaxation components due to specific subcomponents of the system, bulk water, bound water, counterions, backbone, bases, and so on, are not found. The various subcomponents may be so strongly coupled that their motions cannot be treated separately.     

Simulation of adsorption of DNA on carbon nanotubes

We report molecular dynamics simulations of DNA adsorption on a single-walled carbon nanotube (SWNT) in an aqueous environment. We have modeled a DNA segment with 12 base pairs (Dickerson dodecamer) and a (8,8) SWNT in water, with counterions to maintain total charge neutrality. Simulations show that DNA binds to the external surface of an uncharged or positively charged SWNT on a time scale of a few hundred picoseconds. The hydrophobic end groups of DNA are attracted to the hydrophobic SWNT surface of uncharged SWNTs, while the hydrophilic backbone of DNA does not bind to the uncharged SWNT. The binding mode of DNA to charged SWNTs is qualitatively different from uncharged SWNTs. The phosphodiester groups of the DNA backbone are attracted to a positively charged SWNT surface while DNA does not adsorb on negatively charged SWNTs. There is no evidence for canonical double-stranded DNA wrapping around either charged or uncharged SWNTs on the very short time scales of the simulations. The adsorption process appears to have negligible effect on the internal stacking structure of the DNA molecule but significantly affects the A to B form conversion of A-DNA. The adsorption of A-DNA onto an uncharged SWNT inhibits the complete relaxation of A-DNA to B-DNA within the time scale of the simulations. In contrast, binding of the A-DNA onto a positively charged SWNT may promote slightly the A to B conversion.     

Stability of DNA containing a structural water mimic in an A-T rich sequence

We describe here the synthesis and properties of A-T rich DNA containing covalently bound water mimics located in the DNA minor groove.     

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.