Local mobility of nucleic acids as determined from crystallographic data. II. Z-form DNA

Nucleic-acid structural deformability deduced from anisotropic displacement parameters

The growing numbers of very well resolved nucleic-acid crystal structures with anisotropic displacement parameters provide an unprecedented opportunity to learn about the natural motions of DNA and RNA. Here we report a new Monte-Carlo approach that takes direct account of this information to extract the distortions of covalent structure, base pairing, and dinucleotide geometry intrinsic to regularly organized double-helical molecules. We present new methods to test the validity of the anisotropic parameters and examine the apparent deformability of a variety of structures, including several A, B, and Z DNA duplexes, an AB helical intermediate, an RNA, a ligand-DNA complex, and an enzyme-bound DNA. The rigid-body parameters characterizing the positions of the bases in the structures mirror the mean parameters found when atomic motion is taken into account. The base-pair fluctuations intrinsic to a single structure, however, differ from those extracted from collections of nucleic-acid structures, although selected base-pair steps undergo conformational excursions along routes suggested by the ensembles. The computations reveal surprising new molecular insights, such as the stiffening of DNA and concomitant separation of motions of contacted nucleotides on opposite strands by the binding of Escherichia coli endonuclease VIII, which suggest how the protein may direct enzymatic action.

DNA structural forms

In the years that have passed since the publication of Wolfram Saenger's classic book on nucleic acid structure (Saenger, 1984), a considerable amount of new data has been accumulated on the range of conformations which can be adopted by DNA. Many unusual species have joined the DNA zoo, including new varieties of two, three and four stranded helices. Much has been learnt about intrinsic DNA curvature, dynamics and conformational transitions and many types of damaged or deformed DNA have been investigated. In this article, we will try to summarise this progress, pointing out the scope of the various experimental techniques used to study DNA structure, and, where possible, trying to discern the rules which govern the behaviour of this subtle macromolecule. The article is divided into six major sections which begin with a general discussion of DNA structure and then present successively, B-DNA, DNA deformations, A-DNA, Z-DNA and DNARNA hybrids. An extensive set of references is included and should serve the reader who wishes to delve into greater detai.

How long does DNA keep the memory of its conformation? A time-dependent canonical correlation analysis of molecular dynamics simulation

The time dependence of the correlation between motions of different parts of DNA is analyzed from a 200 ps molecular dynamics simulation of the double-stranded self-complementary d(CTGATCAG) in the B form. Each nucleotide is decomposed into three subunits corresponding to the furanose ring (SU), the base (BA), and the backbone (SK). The motion of each subunit is considered as the superimposition of rigid body translation, rigid body rotation, and internal deformation. Canonical time-dependent correlation functions calculated with coordinates describing the different components of the subunits motion are defined and computed. This allows us to probe how long a particular type of motion of one subunit influences the other types of motions of other subunits (cross correlation functions) or how long a particular subunit keeps the memory of its own conformation or location (autocorrelation functions). From autocorrelation analysis it is found that deformation decorrelates within a few tenths of picoseconds, rotational correlation times are on the order of 8 ps, while translational motions are long-time correlated. The deformation of a subunit is not correlated to the deformation of another one (at the 200 ps time scale of our simulation), but influences slightly their translation and orientation as time increases.

Rigid-body motions of sub-units in DNA: a correlation analysis of a 200 ps molecular dynamics simulation

A 200 ps molecular dynamics simulation of the B-form double stranded self-complementary octanucleotide d(CTGATCAG) is analyzed in terms of correlated motions using the canonical analysis approach. Each nucleotide is decomposed in three sub-units corresponding to the base, the sugar ring and the backbone respectively. The correlation between the full dynamics of two sub-units was found to decrease as their mutual distance increases. The interpretation of the full dynamics of sub-units as the superimposition of rigid-body motions (translation and orientation) and deformation shows that the main source of correlation is rigid-body motions. Correlation between sub-units deformation is weak and practically vanishes for sub-units belonging to non-adjacent nucleotides. It is also shown that the correlation is much more important for sub-units of the same strand than of opposite strands. We conclude that the internal dynamics of the octanucleotide may be well described by rigid-body motions, the sub-units deformation having only local influence whereas sub-units translation and rotation have repercussion to long distances. The results presented in this study suggest how the number of degrees of freedom may be reduced for simulating long-time dynamics of oligonucleotides.

Local mobility of nucleic acids as determined from crystallographic data. III. A daunomycin-DNA complex

The local mobility of the complex between the anti-tumor drug daunomycin and a DNA hexanucleotide duplex of sequence d(CpGpTpApCpG)2 has been determined by anisotropic refinement of single crystal X-ray diffraction data of 1.2 A resolution (1 A = 0.1 nm). The directions and amplitudes of the local motion indicate that changes in mobility of DNA due to daunomycin binding are primarily limited to the residues forming the intercalation site and do not propagate to the neighboring residues. The intercalated daunomycin ring system (aglycone) is rigidly fixed in the base stack, apparently serving as an anchor for the amino sugar segment of the drug which is one of the most mobile regions of the entire complex. The high flexibility of this amino sugar may be important for inhibition of replication and transcription not only by sterically blocking the minor groove, but also by allowing nonproductive interactions to be formed with various polymerases or other DNA-binding proteins. The crystallographic model is improved sufficiently by the rigid group anisotropic refinement to allow additional bound water molecules to be located.

Mobility of individual 5-fluorouridine residues in 5-fluorouracil-substituted Escherichia coli valine transfer RNA. A 19F nuclear magnetic resonance relaxation study

19F nuclear magnetic resonance (n.m.r.) relaxation parameters of 5-fluorouracil-substituted Escherichia coli tRNA(Val)1 were measured and used to characterize the internal mobility of individual 5-fluorouridine (FUrd) residues in terms of several models of molecular motion. Measured relaxation parameters include the spin-lattice (T1) relaxation time at 282 MHz, the 19F(1H) NOE at 282 MHz, and the spin-spin (T2) relaxation time, estimated from linewidth data at 338 MHz, 282 MHz and 84 MHz. Dipolar and chemical shift anisotropy contributions to the 19F relaxation parameters were determined from the field-dependence of T2. The results demonstrate a large chemical shift anisotropy contribution to the 19F linewidths at 282 and 338 MHz. Analysis of chemical shift anisotropy relaxation data shows that, relative to overall tumbling of the macromolecule, negligible torsional motion occurs about the glycosidic bond of FUrd residues in 19F-labeled tRNA(Val)1, consistent with the maintenance of base-base hydrogen-bond and/or stacking interactions at all fluorouracil residues in the molecule. The dipolar relaxation data are analyzed by using the "two-state jump" and "diffusion in a cone" formalisms. Motional amplitudes (theta) are interpreted as being due to pseudorotational fluctuations within the ribose ring of the fluorinated nucleoside. These amplitudes range from approximately 30 degrees to 60 degrees, assuming a correlation time (tau i,2) of 1.6 ns. By using available 19F n.m.r. assignment data for the 14 FUrd residues in 5-fluorouracil-substituted tRNA(Val)1, these motional amplitudes can be correlated directly with the environmental domain of the residue. Residues located in tertiary and helical structural domains show markedly less motion (theta approximately equal to 30 to 35 degrees) than residues located in loops (theta approximately equal to 45 to 60 degrees). A correlation between residue mobility and solvent exposure is also demonstrated. The amplitudes of internal motion for specific residues agree quite well with those derived from X-ray diffraction and molecular dynamics data for yeast tRNA(Phe).

Comparison of the reactivity of B-DNA and Z-DNA with two isosteric chemical carcinogens: 2-N,N-acetoxyacetylaminofluorene and 3-N,N-acetoxyacetylamino-4,6-dimethyldipyrido-[1,2-a:3',2' -d] imidazole

The reactivity of nucleic acids in various conformations and two isosteric chemical carcinogens 2-N,N-acetoxyacetylaminofluorene (N-AcO-AAF) and 3-N,N-acetoxyacetylamino-4,6-dimethyldipyrido [1,2-a:3',2'-d] imidazole (N-AcO-AGlu-P-3) have been studied. Both carcinogens bind covalently to poly(dG-dC).poly(dG-dC) (B form) and to poly(dG-br5C).poly(dG-br5dC) (Z form). They also bind covalently to (dC-dG)16 and to (dG-dT)15 sequences inserted in plasmids when the inserts are in the B form but they do not bind to the inserts in the Z form. The reactivity of guanine residues at the B-Z junctions depends upon the superhelical density of the plasmids and upon the base sequences at the junction. The distribution of AGlu-P-3 modified guanines in a restriction fragment of pBR322 is not uniform and is different from that of AAF-modified guanines. The conclusion is that N-AcO-Glu-P-3 as N-AcO-AAF can probe at the nucleotide level the polymorphism of DNA. On the other hand, the non-reactivity of both chemical carcinogens and Z-DNA and the hyperreactivity of some junctions might have some importance in the understanding of chemical carcinogenesis.

Temperature-dependent molecular dynamics and restrained X-ray refinement simulations of a Z-DNA hexamer

Molecular dynamics simulations of the Z-DNA hexamer 5BrdC-dG-5BrdC-dG-5BrdC-dG were performed at several temperatures between 100 K and 300 K. Above 250 K, a strong sequence-dependent flexibility in the nucleic acid is observed, with the guanine sugar and the phosphate of GpC sequences much more mobile than the cytosine sugar and phosphate of CpG sequences. At 300 K, the hexamer is in dynamic equilibrium between several Z forms, including the crystallographically determined ZI and ZII forms. The local base-pair geometry, however, is not very variable, except for the roll of the base-pairs. The hexamer molecular dynamics trajectories have been used to test the restrained parameter crystallographic refinement model for nucleic acids. X-ray diffraction intensities corresponding to observed diffraction data were computed. The average structures obtained from the simulations were then refined against the calculated intensities, using a restrained least-squares program developed for nucleic acids in order to analyse the effects of the refinement model on the derived quantities. In general, the temperature dependence of the atomic fluctuations determined directly from the refined Debye-Waller factors is in reasonably good agreement with the results obtained by calculating the atomic fluctuations directly from the Z-DNA molecular dynamics trajectories. The agreement is best for refinement of temperature factors without restraints. At the highest temperature studied (300 K), the effect of the refinement on the most mobile atoms (phosphates) is to significantly reduce the mean-square atomic fluctuations estimated from the refined Debye-Waller factors below the actual values (less than (delta r)2 greater than congruent to 0.5 A2). Analysis of the temperature-dependence of the mean-square atomic fluctuations provides information concerning the conformational potential within which the atoms move. The calculated temperature-dependence and anharmonicity of the Z-DNA helix are compared with the results observed for proteins. The average structures from the simulations were refined against the experimental X-ray intensities. It is found that low-temperature molecular dynamics simulations provide a useful tool for optimizing the refinement of X-ray structures.