DNA and its counterions: a molecular dynamics study

The behaviour of mobile counterions, Na⁺ and K⁺, was analysed around a B-DNA double helix with the sequence CCATGCGCTGAC in aqueous solution during two 50 ns long molecular dynamics trajectories. The movement of both monovalent ions remains diffusive in the presence of DNA. Ions sample the complete space available during the simulation time, although individual ions sample only about one-third of the simulation box. Ions preferentially sample electronegative sites around DNA, but direct binding to DNA bases remains a rather rare event, with highest site occupancy values of <13%. The location of direct binding sites depends greatly on the nature of the counterion. While Na⁺ binding in both grooves is strongly sequence-dependent with the preferred binding site in the minor groove, K⁺ mainly visits the major groove and binds close to the centre of the oligomer. The electrostatic potential of an average DNA structure therefore cannot account for the ability of a site to bind a given cation; other factors must also play a role. An extensive analysis of the influence of counterions on DNA conformation showed no evidence of minor groove narrowing upon ion binding. A significant difference between the conformations of the double helix in the different simulations can be attributed to extensive α/γ transitions in the phosphate backbone during the simulation with Na⁺. These transitions, with lifetimes over tens of nanoseconds, however, appear to be correlated with ion binding to phosphates. The ion-specific conformational properties of DNA, hitherto largely overlooked, may play an important role in DNA recognition and binding.     

Water permeation through gramicidin A: desformylation and the double helix: a molecular dynamics study

Multinanosecond molecular dynamics simulations of gramicidin A embedded in a dimyristoylphosphatidylcholine bilayer show a remarkable structural stability for both experimentally determined conformations: the head-to-head helical dimer and the double helix. Water permeability was found to be much higher in the double helical conformation, which is explained by lower hydrogen bond-mediated enthalpic barriers at the channel entrance and its larger pore size. Free-energy perturbation calculations show that the double helical structure is stabilized by the positive charges at the N termini introduced by the desformylation, whereas the helical dimer is destabilized. Together with the recent experimental observation that desformyl gramicidin conducts water hundredfold better than gramicidin, this suggests that desformyl gramicidin A predominantly occurs in the double helical conformation.     

Rotational dynamics of adenine amino groups in a DNA double helix

Exocyclic amino groups of the bases undergo conformational fluctuations that affect the recognition and reactivity of nucleic acid molecules. Among these fluctuations, rotation of amino groups around C-N bonds is of special interest. In the present paper, we report the first determination of the rates and energetic parameters for rotation of the N6-amino group of adenine in a DNA double helix. The DNA molecule studied is the dodecamer [d(CGCGAGCTCGCG)]2. The adenine in each A. T basepair of the dodecamer was labeled with 15N at the N6 position, and the NMR resonances of the two protons in the adenine amino group were selectively observed by 15N-editing methods. The rates of rotation of the amino group were obtained from experiments of transfer of magnetization between the two protons in the same group and from lineshape analysis of 15N-edited amino proton NMR resonances. The results show that, over the temperature range from 0 to 70 degrees C, the rates of rotation of adenine amino groups range from 60 to 24,000 s-1. Formation of the activated state during rotation has a standard enthalpy change of 15.3 +/- 0.2 kcal/mol and a standard entropy change of 6.0 +/- 0.7 cal/(mol. K). Analysis of the results suggests that rotation of the amino group occurs in the paired, closed state of the adenine in the A. T basepair of the double-helical DNA structure.     

Hydration structure and dynamics of B- and Z-DNA in the presence of counterions via molecular dynamics simulations

Following our previous attempts at understanding the structural and dynamical properties of water and counterions hydrating nucleic acids, we have performed molecular dynamics simulations for B- and Z-DNA. In these simulations, the nucleic acids were held rigid. In the case of B-DNA, one turn of B-DNA double helix was considered in the presence of 1500 water molecules and 20 counterions (K⁺). The simulations were performed for 4.0 ps after equilibrating the system. For Z-DNA, we considered one turn of the double helix in the presence of 1851 water molecules and 24 counterions (K⁺). The simulations were carried out for 3.5 ps after equilibration. The average temperature of these simulations was ～ 360 K for Z-DNA and ～ 345 K for B-DNA. In these simulations the hydrogen atoms were explicitly taken into account. For both simulations, a fifth-order predictor-corrector was used for solving the translational equations of motion. The rotational motion of the water molecules was represented in terms of quaternion algebra and the rotational equations of motion were solved with a second-order quaternion method using a sixth-order predictor-corrector method. A time step of 0.5 · 10^?15 s was used in these simulations. The structural and the dynamical properties of water solvating the counterions, and the phosphate groups of the DNA, were computed to understand the hydration structure. Diffusion coefficients and velocity correlation functions were calculated for both ions and the water molecules. The velocity correlation functions for the ions exhibit a caged behavior. The dipole correlation functions for the water molecules indicate that the water molecules close to the helix retain the memory of their initial orientations for longer periods of time than those away from the helix. During the time period of our simulation (3–4 ps) the ion probability distributions show a well-defined pattern and suggest limited mobility for the ions, being close to the helix.     

Localisation and dynamics of sodium counterions around DNA in solution from molecular dynamics simulation

The localisation and dynamics of sodium counterions around the DNA duplex d(AGCGTACTAGTACGCT)₂ corresponding to the trp operator fragment used in the crystal structure of the half site complex has been studied by a 1.4 ns molecular dynamics simulation in explicit solvent. A continuous and well-defined counterion density is shown to be present around the minor groove, while density patches are found in the major groove in regions where DNA bending is observed. A residence time analysis reveals the dynamic nature of these distributions. The resulting picture agrees with previous theoretical and experimental studies of A-tract DNA sequences, and is consistent with the polyelectrolyte condensation model. Received: 12 August 1999 / Revised version: 11 November 1999 / Accepted: 7 December 1999     

Hydration pattern of A4T4 and T4A4 DNA: a molecular dynamics study

Hydration pattern and energetics of 'A-tract' containing duplexes have been studied using molecular dynamics on 12-mer self-complementary sequences 5'-d(GCA4T4GC)-3' and 5'-d(CGT4A4CG)-3'. The structural features for the simulated duplexes showed correlation with the corresponding experimental structures. Analysis of the hydration pattern confirmed that water network around the simulated duplexes is more conformation specific rather than sequence specific. The calculated heat capacity change upon duplex formation showed that the process is entropically driven for both the sequences. Furthermore, the theoretical free energy estimates calculated using MMPBSA approach showed a higher net electrostatic contribution for A4T4 duplex formation than for T4A4, however, energetically both the duplexes are observed to be equally stable.     

Perpetuating the double helix: molecular machines at eukaryotic DNA replication origins

The hardest part of replicating a genome is the beginning. The first step of DNA replication (called "initiation") mobilizes a large number of specialized proteins ("initiators") that recognize specific sequences or structural motifs in the DNA, unwind the double helix, protect the exposed ssDNA, and recruit the enzymatic activities required for DNA synthesis, such as helicases, primases and polymerases. All of these components are orderly assembled before the first nucleotide can be incorporated. On the occasion of the 50th anniversary of the discovery of the DNA structure, we review our current knowledge of the molecular mechanisms that control initiation of DNA replication in eukaryotic cells, with particular emphasis on the recent identification of novel initiator proteins. We speculate how these initiators assemble molecular machines capable of performing specific biochemical tasks, such as loading a ring-shaped helicase onto the DNA double helix.     

Hydration of the dTn.dAn x dTn parallel triple helix: a Fourier transform infrared and gravimetric study correlated with molecular dynamics simulations.

We present a comparative analysis of the water organization around the dTn.dAn x dTn triple helix and the Watson-Crick double helix dTn.dAn respectively by means of gravimetric measurements, infrared spectroscopy and molecular dynamics simulations. The hydration per nucleotide determined by gravimetric and spectroscopic methods correlated with the molecular dynamics simulations shows that at high relative humidity (98% RH) the triple helix is less solvated than the duplex (17 +/- 2 water molecules per nucleotide instead of 21 +/-1). The experimental desorption curves are different for both structures and indicate that below 81% RH the triplex becomes more hydrated than the duplex. At this RH the FTIR spectra show the emergence of N-type sugars in the adenosine strand of the triplex. When the third strand is bound in the major groove of the Watson-Crick duplex molecular dynamics simulations show the formation of a spine of water molecules between the two thymidine strands.     

A molecular dynamics simulation study of oriented DNA with polyamine and sodium counterions: diffusion and averaged binding of water and cations

Four different molecular dynamics (MD) simulations have been performed for ordered DNA decamers, d(5′-ATGCAGTCAG)·d(5′-TGACTGCATC). The counterions were the two natural polyamines spermidine³⁺ (Spd³⁺) and putrescine²⁺ (Put²⁺), the synthetic polyamine diaminopropane²⁺ (DAP²⁺) and Na⁺. The simulation set-up corresponds to an infinite array of parallel DNA mimicking the state in oriented DNA fibers or crystals. This work describes general properties of polyamine and Na⁺ binding to DNA. Simulated diffusion coefficients show satisfactory agreement with experimental NMR diffusion data of comparable systems. The interaction of the polyamines with DNA is dynamic in character and the cations mostly form short-lived contacts with the electronegative binding sites of DNA. Polyamines, Na⁺ and water interact most frequently with the charged phosphate atoms with preference for association from the minor groove side with O1P over O2P. There is a strong anti-correlation in the cation binding to the electronegative groups of DNA, i.e. the presence of a cation near one of the DNA sites repels other cations from binding to this and to the other sites separated by <7.5 Å from each other. In contrast to the other polyamines, DAP²⁺ is able to form ‘bridges’ connecting neighboring phosphate groups along the DNA strand. A small fraction of DAP²⁺ and Put²⁺ can be found in the major grooves, while Spd³⁺ is absent there. The results of the MD simulations reveal principal differences in the polyamine–DNA interactions between the natural (Spd³⁺, Put²⁺ and spermine⁴⁺) and synthetic (DAP²⁺) polyamines.     

Ethidium ion binds more strongly to a DNA double helix with a bulged cytosine than to a regular double helix

Thermodynamic parameters for ethidium intercalation were determined for the double helices formed by the oligonucleotides dCA6G + dCT6G, which form a normal helix, and dCA3CA3G + dCT6G, which form a double helix with the middle cytosine bulged outside of the helix. Ethidium intercalation was measured by monitoring the absorbance at 260 and 283 nm as a function of temperature for a number of concentrations of ethidium. The binding to the normal helix occurs equally at all the intercalation sites, with an enthalpy of binding of -8 kcal mol-1, an entropy of binding of -6 eu, and an equilibrium constant at 25 degrees C of 2.2 X 10(4) M-1. The binding to the bulged double helix was considerably stronger and is consistent with a model in which the intercalation sites on either side of the bulged base bind 10 times stronger than the other sites. Thus, there are two strong binding sites on the perturbed helix with equilibrium constants for binding of 2 X 10(5) M-1 at 25 degrees C in addition to five normal sites. Several other binding models were tested but did not fit the data satisfactorily.     

Hydration dynamics of the collagen triple helix by NMR

The hydration of the collagen-like Ac-(Gly-Pro-Hyp)(6)-NH(2) triple-helical peptide in solution was investigated using an integrated set of high-resolution NMR hydration experiments, including different recently developed exchange-network editing methods. This approach was designed to explore the hydration dynamics in the proximity of labile groups, such as the hydroxyproline hydroxyl group, and revealed that the first shell of hydration in collagen-like triple helices is kinetically labile with upper limits for water molecule residence times in the nanosecond to sub-nanosecond range. This result is consistent with a "hopping" hydration model in which solvent molecules are exchanged in and out of solvation sites at a rate that is not directly correlated to the degree of site localization. The hopping model thus reconciles the dynamic view of hydration revealed by NMR with the previously suggested partially ordered semi-clathrate-like cylinder of hydration. In addition, the nanosecond to sub-nanosecond upper limits for water molecule residence times imply that hydration-dehydration events are not likely to be the rate-limiting step for triple helix self-recognition, complementing previous investigations on water dynamics in collagen fibers. This study has also revealed labile proton features expected to facilitate the characterization of the structure and folding of triple helices in collagen peptides.     

Formation of the double helix: a mutational study

To investigate the mechanisms by which oligonucleotides hybridize to target molecules, the binding of two oligodeoxynucleotide probes to RNA targets was measured over a broad range of temperatures. Mutations were then scanned across each DNA/RNA hybrid to map, at single base resolution, sequences important for hybridization. Despite being unrelated in sequence, each hybrid formed by a similar mechanism. In the absence of secondary structure, two stretches of bases, termed nucleation regions, cooperated with one another by a looping mechanism to nucleate hybridization. Mutations inside each nucleation region strongly decreased hybridization rates, even at temperatures well below the melting temperature (T_m) of the hybridized duplex. Surprisingly, nucleation regions were detected in a RNA target but not a corresponding DNA target. When either nucleation region was sequestered in secondary structure, the hybridization rate fell and the mechanism of hybridization changed. Single-stranded bases within the nucleation region of the probe and target first collided to form a double helix. If sufficiently G + C rich, the double helix then propagated throughout the oligonucleotide by a strand invasion process. On the basis of these results, general mechanisms for the hybridization of oligonucleotides to complementary and mutant targets are proposed.     

Thermodynamics of site-specific small molecular ion interactions with DNA duplex: a molecular dynamics study

The stability and dynamics of a double-stranded DNA (dsDNA) is affected by the preferential occupancy of small monovalent molecular ions. Small metal and molecular ions such as sodium and alkyl ammonium have crucial biological functions in human body, affect the thermodynamic stability of the duplex DNA and exhibit preferential binding. Here, using atomistic molecular dynamics simulations, we investigate the preferential binding of metal ion such as Na⁺ and molecular ions such as tetramethyl ammonium (TMA⁺) and 2-hydroxy-N,N,N-trimethylethanaminium (CHO⁺) to double-stranded DNA. The thermodynamic driving force for a particular molecular ion-DNA interaction is determined by decomposing the free energy of binding into its entropic and enthalpic contributions. Our simulations show that each of these molecular ions preferentially binds to the minor groove of the DNA and the extent of binding is highest for CHO⁺. The ion binding processes are found to be entropically favourable. In addition, the contribution of hydrophobic effects towards the entropic stabilisation (in case of TMA⁺) and the effect of hydrogen bonding contributing to enthalpic stabilisation (in case of CHO⁺) have also been investigated.     

A hydrogen exchange study of the open segment in a DNA double helix

The kinetics of the hydrogen-deuterium exchange reactions of double-helical poly[d(A-T)]·poly[d(A-T)], poly(dA)·poly(dT), and constituent nucleosides (deoxyadenosine and thymidine) have been examined at various temperatures by stopped-flow ultraviolet spectrophotometry, in the spectral region 240–300 nm. The results were interpreted on the basis of a mechanism of the hydrogen exchange reaction of a helical polynucleotide, proposed by Englander and colleagues as well as by the Tsuboi and Nakanishi group. It was concluded that the rates of the base-pair opening reactions are nearly equal to one another in double-helical DNAs, irrespective of the base sequence. On the other hand, the free energy required for bringing the open segment at a particular base-pair was found to be much greater for poly(dA)·poly(dT) than for poly[d(A-T)]· poly[d(A-T)].     

Conjugated double bonds in lipid bilayers: a molecular dynamics simulation study

Conjugated linoleic acids (CLA) are found naturally in dairy products. Two isomers of CLA, that differ only in the location of cis and trans double bonds, are found to have distinct and different biological effects. The cis 9 trans 11 (C9T11) isomer is attributed to have the anti-carcinogenic effects, while the trans 10 cis 12 (T10C12) isomer is believed to be responsible for the anti-obesity effects. Since dietary CLA are incorporated into membrane phospholipids, we have used Molecular Dynamics (MD) simulations to investigate the comparative effects of the two isomers on lipid bilayer structure. Specifically, simulations of phosphatidylcholine lipid bilayers in which the sn-2 chains contained one of the two isomers of CLA were performed. Force field parameters for the torsional potential of double bonds were obtained from ab initio calculations. From the MD trajectories we calculated and compared structural properties of the two lipid bilayers, including areas per molecule, density profiles, thickness of bilayers, tilt angle of tail chains, order parameters profiles, radial distribution function (RDF) and lateral pressure profiles. The main differences found between bilayers of the two CLA isomers, are (1) the order parameter profile for C9T11 has a dip in the middle of sn-2 chain while the profile for T10C12 has a deeper dip close to terminal of sn-2 chain, and (2) the lateral pressure profiles show differences between the two isomers. Our simulation results reveal localized physical structural differences between bilayers of the two CLA isomers that may contribute to different biological effects through differential interactions with membrane proteins or cholesterol.     

Relative flexibility of DNA and RNA: a molecular dynamics study

State of the art molecular dynamics simulations are used to study the structure, dynamics, molecular interaction properties and flexibility of DNA and RNA duplexes in aqueous solution. Special attention is paid to the deformability of both types of structures, revisiting concepts on the relative flexibility of DNA and RNA duplexes. Our simulations strongly suggest that the concepts of flexibility, rigidity and deformability are much more complex than usually believed, and that it is not always true that DNA is more flexible than RNA.