The speed of sound in DNA

We have used Brillouin scattering to determine the speed of sound in (and hence longitudinal modulus of) A- and B-DNA fibers. The speed of sound is very sensitive to the degree of hydration of the fibers, and measurements have to be made at laser powers below 5 mW to avoid local heating and dehydration. Under those conditions, we obtain sound speed perpendicular to the fiber axis of about 2.2 and 1.9 km/s in A- and B-DNA fibers, respectively. A-DNA fibers show a small anisotropy with sound speeds along the fiber axis higher by up to 10% B-DNA fibers appear to be isotropic.     

Characterization of DNA structures by laser Raman spectroscopy

Oriented fibers drawn from aqueous gels of calf-thymus DNA were maintained at constant relative humidites of 75 and 92% to yield canonical A-DNA and B-DNA structures, respectively. Raman spectra of the two forms of DNA were recorded over the spectral range 300–4000 cm^?1. The authenticated DNA fibers were deuterated in hygrostatic cells containing D₂O at appropriate relative humidities, and the corresponding spectra of deuterated DNAs were also obtained. The spectra reveal all of the Raman scattering frequencies and intensities characteristic of A- and B-DNA structures in both nondeuterated and deuterated froms, as well as the frequencies and intensities of adsorbed solvent molecules from which the hydration content of DNA fibers can be calculated. Numerous conformation-sensitive vibrational modes of DNA bases and phosphate groups have been identified throughout the 300–1700-cm^?1 interval. Evidence has also been obtained for conformation sensitivity of deoxyribosyl CH stretching modes in the 2800–3000-cm^?1 region. Raman lines of both the backbone and the bases are proposed as convenient indicators of A- and B-DNA structures. The results are extended to Z-DNA models investigated previously. Some implications of these findings for the determination of DNA or RNA structure from Raman spectra of nucleoproteins and viruses are considered.     

Mechanochemical study of NaDNA and NaDNA-netropsin fibers in ethanol-water and trifluoroethanol-water solutions.

Highly oriented calf-thymus NaDNA fibers, prepared by a wet-spinning method, were complexed with netropsin in ethanol-water and trifluoroethanol (TFE)-water solutions. The relative fiber length, L/L0, was measured at room temperature as a function of ethanol or TFE concentration to obtain information on the B-A conformational transition. The B-A transition point and transition cooperativity of the fibers were calculated. The binding of netropsin to NaDNA fibers was found to stabilize B form and to displace the B-A transition to higher ethanol concentration, as indicated by its elongational effect on the fiber bundles. An increased salt concentration was found to reduce netropsin binding. In netropsin-free ethanol solution, the dissociation of bound netropsin from the DNA fibers was observable. Pure B-NaDNA fibers were found to be more stable in TFE solution than in ethanol solution. This was interpreted as being due to a different steric factor and a larger polarity of TFE compared with ethanol, resulting in its smaller capacity to reduce the water activity and dielectric constant of the medium in the immediate vicinity of DNA fibers. Therefore, the effect of netropsin binding on the B-A transition of NaDNA fibers became less obvious in TFE solution. In another series of experiments, L/L0 was measured as a function of temperature to obtain information on the helix-coil transition, or melting, as well as the B-A transition of NaDNA and NaDNA-netropsin fibers. The melting temperature and helix-coil transition width were calculated from the melting curves. A phenomenological approach was used to describe the melting behavior of the fibers in and around the B-A transition region. The effect of netropsin on the melting of DNA fibers was attributed mainly to the stabilization of B-DNA and to a higher melting cooperativity in the B-DNA region.     

Modelling and refinement of the crystal structure of nucleoprotamine from Gibbula divaricata

The molecular structure of nucleoprotamine from Gibbula divaricata and its packing in oriented fibers has been modelled both to fit the X-ray diffraction pattern and to avoid steric compression. The representative model consists of 51 poly (dinucleotide) B-DNA helices with 51 poly(hexapeptide) chains associated with the major grooves. The prevailing peptide conformation is beta. The four arginine residues present are hydrogen-bonded to DNA phosphates while neutral peptides protrude into the minor grooves of neighboring nucleoprotamine molecules which are packed 2.61 nm apart in a screw-disordered, quasi-hexagonal lattice. This model reconciles a number of earlier, apparently conflicting experimental results and explains the remarkable stability of nucleoprotamines.     

Mechanochemical study of conformational transitions and melting of Li-, Na-, K-, and CsDNA fibers in ethanol–water solutions

Highly oriented fibers of Li-, Na-, K-, and CsDNA were prepared with a previously developed wet spinning method. The procedure gave a large number of equivalent fiber bundle samples (reference length, L₀, typically = 12–15 cm) for systematic measurements of the fiber length L in ethanol–water solutions, using a simple mechanochemical set up. The decrease in relative length L/L₀ with increasing ethanol concentration at room temperature gave evidence for the B-A transition centered at 76% (v/v) ethanol for NaDNA fibers and at 80 and 84% ethanol for K- and CsDNA fibers. A smaller decrease in L/L₀ of LiDNA fibers was attributed to the B-C transition centered at 80% ethanol. In a second type of experiment with DNA fibers in ethanol–water solutions, the heat-induced helix–coil transition, or melting, revealed itself in a marked contraction of the DNA fibers. The melting temperature T_m, decreased linearly with increasing ethanol concentration for fibers in the B-DNA ethanol concentration region. In the B-A transition region, Na- and KDNA fibers showed a local maximum in T_m. On further increase of the ethanol concentration, the A-DXA region followed with an even steeper linear decrease in T_m. The dependence on the identity of the counterion is discussed with reference to the model for groove binding of cations in B-DNA developed by Skuratovskii and co-workers and to the results from Raman studies of the interhelical bonds in A-DNA performed by Lindsay and co-workers. An attempt to apply the theory of Chogovadze and Frank-Kamenetskii on DNA melting in the B-A transition region to the curves failed. However, for Na- and KDNA the T_m dependence in and around the A-B transition region could be expressed as a weighted mean value of T_m of A- and B-DNA. On further increase of the ethanol concentration, above 84% ethanol for LiDNA and above about 90% ethanol for Na-, K-, and CsDNA, a drastic change occurred. T_m increased and a few percentages higher ethanol concentrations were found to stabilize the DNA fibers so that they did not melt at all, not even at the upper temperature limit of the experiments (～ 80°C). This is interpreted as being due to the strong aggregation induced by these high ethanol concentrations and to the formation of P-DNA. Many features of the results are compatible with the counterion–water affinity model. In another series of measurements, T_m of DNA fibers in 75% ethanol was measured at various salt concentrations. No salt effect was observed (with the exception of LiDNA at low salt concentrations). This result is supported by calculations within the Poisson–Boltzmann cylindrical cell model. © 1994 John Wiley & Sons, Inc.     

Crystal and solution structures of the B-DNA dodecamer d(CGCAAATTTGCG) probed by Raman spectroscopy: heterogeneity in the crystal structure does not persist in the solution structure

The self-complementary dodecamer d(CGCAAATTTGCG) crystallizes as a double helix of the B form and manifests a Raman spectrum with features not observed in Raman spectra of either DNA solutions or wet DNA fibers. A number of Raman bands are assigned to specific nucleoside sugar and phosphodiester conformations associated with this model B-DNA crystal structure. The Raman bands proposed as markers of the crystalline B-DNA structure are compared and contrasted with previously proposed markers of Z-DNA and A-DNA crystals. The results indicate that the three canonical forms of DNA can be readily distinguished by Raman spectroscopy. However, unlike Z-DNA and A-DNA, which retain their characteristic Raman fingerprints in aqueous solution, the B-DNA Raman spectrum is not completely conserved between crystal and solution states. The Raman spectra reveal greater heterogeneity of nucleoside conformations (sugar puckers) in the DNA molecules of the crystal structure than in those of the solution structure. The results are consistent with conversion of one-third of the dG residues from the C2'-endo/anti conformation in the solution structure to another conformation, deduced to be C1'-exo/anti, in the crystal. The dodecamer crystal also exhibits unusually broad Raman bands at 790 and 820 cm-1, associated with the geometry of the phosphodiester backbone and indicating a wider range of (alpha, zeta) backbone torsion angles in the crystal than in the solution structure. The results suggest that backbone torsion angles in the CGC and GCG sequences, which flank the central AAATTT sequence, are significantly different for crystal and solution structures, the former containing the greater diversity.(ABSTRACT TRUNCATED AT 250 WORDS)     

Internal correlations in minor groove profiles of experimental and computed B-DNA conformations.

It has been noticed that converged conformations of B-DNA oligomers obtained in MD calculations often have very small atom position rmsd values from the canonical B-DNA and all helical parameters close to the standard values, but their minor grooves tend to be somewhat narrower. This apparent bias disappears, however, when C5' rather than phosphorus atoms are used for measuring the groove width. At the origin of this effect is the specific orientation of phosphate groups in the canonical B-DNA model which maximizes their separation across the minor groove. When measured by C5' traces, minor groove profiles of experimental structures available in the Nucleic Acids Database show much less tendency to narrow below the canonical width. Correlation analysis reveals a high degree of correspondence in shapes of minor grooves of calculated and experimental single-crystal structures of B-DNA oligomers.     

Energetics of left and right handed models of DNA

It has been shown by model building studies that various right handed and left handed models are compatible with X-ray data of B-DNA and C-DNA. These models are also found to be in good agreement with infrared dichroism data. Detailed potential energy calculations have now been carried out for these models, viz., right and left handed B-DNA and right and left handed C-DNA. It is found that base sugar stacking and interactions involving the phosphate groups are the dominant forces for stabilizing a particular structure. For some sequences, viz., A-A, T-A and C-A, left handed stacking is quite favourable in both B and C structures. But intranucleotide interactions make the left B-DNA unfavourable while the left C-DNA structure is more stable, for all the sequences, than the right C-DNA structure, proposed from fibre data. For the hexanucleoside pentaphosphate fragments the same trend is observed, with the right handed B-DNA being the most stable of the four models studied. However, the left C-DNA structure is only marginally higher in energy, particularly if the shielding effect of the counter ions, on the phosphate group is taken into consideration.     

Predicting accurate ab initio DNA electron densities with equivariant neural networks

One of the fundamental limitations of accurately modeling biomolecules like DNA is the inability to perform quantum chemistry calculations on large molecular structures. We present a machine learning model based on an equivariant Euclidean neural network framework to obtain accurate ab initio electron densities for arbitrary DNA structures that are much too large for conventional quantum methods. The model is trained on representative B-DNA basepair steps that capture both base pairing and base stacking interactions. The model produces accurate electron densities for arbitrary B-DNA structures with typical errors of less than 1%. Crucially, the error does not increase with system size, which suggests that the model can extrapolate to large DNA structures with negligible loss of accuracy. The model also generalizes reasonably to other DNA structural motifs such as the A- and Z-DNA forms, despite being trained on only B-DNA configurations. The model is used to calculate electron densities of several large-scale DNA structures, and we show that the computational scaling for this model is essentially linear. We also show that this machine learning electron density model can be used to calculate accurate electrostatic potentials for DNA. These electrostatic potentials produce more accurate results compared with classical force fields and do not show the usual deficiencies at short range.     

Modelling and Refinement of the Crystal Structure of Nucleoprotamine from Gibbuta Divaricata

Abstract

The molecular structure of nucleoprotamine from Gibbula divaricata and its packing in oriented fibers has been modelled both to fit the X-ray diffraction pattern and to avoid steric compression. The representative model consists of 5₁ poly(dinucleotide) B-DNA helices with 5₁ poly(hexapeptide) chains associated with the major grooves. The prevailing peptide conformation is β, The four arginine residues present are hydrogen-bonded to DNA phosphates while neutral peptides protrude into the minor grooves of neighboring nucleoprotamine molecules which are packed 2.61 nm apart in a screw-disordered, quasi-hexagonal lattice. This model reconciles a number of earlier, apparently conflicting experimental results and explains the remarkable stability of nucleoprotamines.     

Approach to the limit of counterion condensation

According to counterion condensation theory, one of the contributions to the polyelectrolyte free energy is a pairwise sum of Debye-Hückel potentials between polymer charges that are reduced by condensed counterions. When the polyion model is taken as an infinitely long and uniformly spaced line of charges, a simple closed expression for the summation, combined with entropy-derived mixing contributions, leads to the central result of the theory, a condensed fraction of counterions dependent only on the linear charge density of the polyion and the valence of the counterion, stable against increases of salt up to concentrations in excess of 0.1 M. Here we evaluate the sum numerically for B-DNA models other than the infinite line of B-DNA charges. For a finite-length line there are end effects at low salt. The condensation limit is reached as a flat plateau by increasing the salt concentration. At a fixed salt concentration the condensation limit is reached by increasing the length of the line. At moderate salt even very short B-DNA line-model oligomers have condensed fractions not far from the infinite polymer limit. For a long double-helical array with charge coordinates at the phosphates of B-DNA, the limiting condensed fraction appears to be approached at low salt. In contrast to the results for the line of charges, however, the computed condensed fraction varies strongly with salt in the range of experimentally typical concentrations. Salt invariance is restored, in agreement with both the line model and experimental data, when dielectric saturation is considered by means of a distance-dependent dielectric function. For sufficiently long B-DNA line and helical models, as typical salt concentrations, the counterion binding fraction approaches the polymer limit as a linear function of 1/P, where P is the number of phosphate groups of B-DNA.     

Energetics of Left and Right Handed Models of DNA

Abstract

It has been shown by model building studies that various right handed and left handed models are compatible with X-ray data of B-DNA and C-DNA. These models are also found to be in good agreement with infrared dichroism data. Detailed potential energy calculations have now been carried out for these models, viz., right and left handed B-DNA and right and left handed C-DNA It is found that base sugar stacking and interactions involving the phosphate groups are the dominant force s for stabilizing a particular structure. For some sequences, viz., A-A, T-A and C-A, left handed stacking is quite favourable in both B and C structures. But intranucleotide interactions make the left B-DNA unfavourable while the left C-DNA structure is more stable, for all the sequences, than the right C-DNA structure, proposed from fibre data. For the hexanucleoside pentaphosphate fragments the same trend is observed, with the right handed B-DNA being the most stable of the four models studied. However, the left C-DNA structure is only marginally higher in energy, particularly if the shielding effect of the counter ions, on the phosphate group is taken into consideration.     

Simulations of the solvent structure for macromolecules. III. Determination of the Na+ counter ion structure

We report on a computer experiment in which, using Monte Carlo techniques, we considered a three-turn (30-base-pairs) B-DNA fragment as a solute and a set of 1200 water molecules and 60 sodium counterions (at a temperature of 300 K) as a solvent. From a statistical analysis of the Monte Carlo simulation (applied to the water molecules and counterions in the B-DNA field), we determined that the counterions themselves conform to two helical structures intertwined with the two strands. The strutures of the water molecules solvating both counterion helices and the two B-DNA strands are fully analyzed and described in detail. A model for base-pair recognition based on the above findings is proposed. Aspects of the unwinding mechanism are discussed.     

Conformational transitions of the phosphodiester backbone in native DNA: two-dimensional magic-angle-spinning 31P-NMR of DNA fibers.

Solid-state 31P-NMR is used to investigate the orientation of the phosphodiester backbone in NaDNA-, LiDNA-, MgDNA-, and NaDNA-netropsin fibers. The results for A- and B-DNA agree with previous interpretations. We verify that the binding of netropsin to NaDNA stabilizes the B form, and find that in NaDNA, most of the phosphate groups adopt a conformation typical of the A form, although there are minor components with phosphate orientations close to the B form. For LiDNA and MgDNA samples, on the other hand, we find phosphate conformations that are in variance with previous models. These samples display x-ray diffraction patterns that correspond to C-DNA. However, we find two distinct phosphate orientations in these samples, one resembling that in B-DNA, and one displaying a twist of the PO4 groups about the O3-P-O4 bisectors. The latter conformation is not in accordance with previous models of C-DNA structure.     

Structure of poly d(A).poly d(T).

On the basis of the x-ray data from polycrystalline and well oriented fibers of the sodium salt of poly d(A).poly d(T) (Arnott et al, Nucl. Acids Res. 11, 4141-4155 (1983), a revised B'-DNA model incorporating B-like adenine and thymine strands is shown to give a much better x-ray agreement (R = 0.25) than the previously assigned model consisting of mixed sugar conformations in the two strands. The narrowing of the minor and the widening of the major grooves are promiscuous features of B'-DNA, which are common to all poly d(purine).poly d(pyrimidine) duplexes with two hydrogen bonded base-pairs and are in marked contrast with classical B-DNA. Due to modest propeller (-15 degrees), the cross strand diagonal hydrogen bonds (0.37 nm) in this duplex are not as strong as those in A,T-rich oligonucleotide crystal structures.     

Structure of a bent DNA: two-dimensional NMR studies on d(GAAAATTTTC)2

Intrinsic DNA bending is caused by specific DNA sequences. The decamer d(GA4T4C)2, when it repeats in a synthetic polymer or in kinetoplast DNA, results in a macroscopic bending of the molecule as a whole. We employed high-resolution two-dimensional NMR methods to examine the intrinsic structural properties of the d(GA4T4C)2 duplex in solution. Examination of the NOESY data at 50- and 100-ms mixing times indicated that the kinds of observed NOEs can originate if each of the ten nucleotidyl residues belongs to the B-DNA family, i.e., C2'-endo,anti. However, the degree of observed NOE intensities from the A-T junction as well as the observed AH2-AH2 cross-peaks from adjacent AT pairs could not be rationalized on the basis of a straight B-DNA model but could be explained by only a B-DNA model with some structural discontinuity at the A-T junction--the site of 2-fold symmetry in the molecule. In view of the fact that the degree of observed NOE intensities can be complicated by spin diffusion and by fine structural distortion, we have resorted to the use of quantitative theoretical NOESY simulation (which takes into account primary, secondary, and higher orders of NOE) to delineate the structural discontinuity at the A-T junction and to arrive at a structure for the duplex d(GA4T4C)2. We propose a "junction B-DNA model" which can quantitatively explain the 2D NOESY data at 100- and 50-ms mixing times. In this model the two structural blocks in the molecule, i.e., d(GA4).d(T4C) and d(T4C).d(GA4), are conformationally equivalent and are connected at the A-T junction where the base pairs are stably stacked, but the two local structural frames do not coincide in space. This model can create an overall bending of 10 degrees with a center of curvature 50 A away from the center of the duplex. It is the thesis of this paper that the observed bending in polymers with a repeat of d(GA4T4C)2 and the bending in natural DNAs where AnTn.AnTn repeats are present originate at the oligonucleotide repeat level.     

Double helix conformation, groove dimensions and ligand binding potential of a G/C stretch in B-DNA.

The self-complementary DNA fragment CCGGCGCCGG crystallizes in the rhombohedral space group R3 with unit cell parameters a = 54.07 A and c = 44.59 A. The structure has been determined by X-ray diffraction methods at 2.2 A resolution and refined to an R value of 16.7%. In the crystal, the decamer forms B-DNA double helices with characteristic groove dimensions: compared with B-DNA of random sequence, the minor groove is wide and deep and the major groove is rather shallow. Local base pair geometries and stacking patterns are within the range commonly observed in B-DNA crystal structures. The duplex bears no resemblance to A-form DNA as might have been expected for a sequence with only GC base pairs. The shallow major groove permits an unusual crystal packing pattern with several direct intermolecular hydrogen bonds between phosphate oxygens and cytosine amino groups. In addition, decameric duplexes form quasi-infinite double helices in the crystal by end-to-end stacking. The groove geometries and accessibilities of this molecule as observed in the crystal may be important for the mode of binding of both proteins and drug molecules to G/C stretches in DNA.     

Simulations of the B-DNA molecular dynamics of d(CGCGAATTCGCG)2 and d(GCGCGCGCGC)2: an analysis of the role of initial geometry and a comparison of united and all-atom models

Molecular dynamics simulations on the sequence d(CGCGAATTCGCG)2 have been carried out using both united atom and all-atom representations, and starting the simulations both from a regular repeating B-DNA structure and from the x-ray single crystal B-DNA structure. An all-atom B-DNA simulation on the sequence d(GCGCGCGCGC)2 has also been carried out, in order to compare it with a previous united atom simulation. The helix repeats, H-bonding, sugar pucker profiles, and average torsional angles are all in the range observed in crystallographic and nmr studies for B-DNA helices. In some of the sequences, there is a significant bend in the DNA helices. The individual helix repeats, with focus on 3'CpG5' and 3'GpC5' units, show the opposite helix repeat to that suggested by Calladine's rules.     

A one- and two-dimensional NMR study of the B to Z transition of (m5dC-dG)3 in methanolic solution

The deoxyribose hexanucleoside pentaphosphate (m5dC-dG)3 has been studied by 500 MHz 1H NMR in D2O (0.1 M NaCl) and in D2O/deuterated methanol mixtures. Two conformations, in slow equilibrium on the NMR time scale, were detected in methanolic solution. Two-dimensional nuclear Overhauser effect (NOE) experiments were used to assign the base and many of the sugar resonances as well as to determine structural features for both conformations. The results were consistent with the an equilibrium in solution between B-DNA and Z-DNA. The majority of the molecules have a B-DNA structure in low-salt D2O and a Z-DNA structure at high methanol concentrations. A cross-strand NOE between methyl groups on adjacent cytosines is observed for Z-DNA but not B-DNA. The B-DNA conformation predominates at low methanol concentrations and is stabilized by increasing temperature, while the Z-DNA conformation predominates at high methanol concentrations and low temperatures. 31P NMR spectra gave results consistent with those obtained by 1H NMR. Comparison of the 31P spectra with those obtained on poly(dG-m5dC) allow assignment of the lower field resonances to GpC in the Z conformation.     

Elasticity theory of the B-DNA to S-DNA transition.

We propose in this note a simple model--the two-state Worm Like Chain--to describe the elasticity of the recently discovered stress-induced transformation from B-DNA to S-DNA. The model reduces for low tractions to the well-known Worm Like chain theory, which is used to describe the elastic properties of B-DNA, while in the limit of high chain-bending moduli it reduces to the two-state Ising model proposed by Cluzel et al. for the B-S transition [Cluzel, P., A. Lebrun, C. Heller, R. Lavery, J-L. Viovy, D. Chatenay, and F. Caron. 1996. DNA: an extensible molecule. Science. 271:792-794]. Our model can be treated analytically to produce an explicit form of the force-extension relationship which agrees reasonably with the observations. We use the model to show that conformational fluctuations of the chain play a role also for the B to S transformation.