Theory of Low-Frequency Vibrations in DNA Macromolecules

Abstract

To describe low-frequency dynamics of DNA macromolecules a model is developed taking into account the hydrogen bond stretching in base pairs, the backbone flexibility and intranucleoside mobility. For double-stranded DNA the normal vibrations are found and the structure of low- frequency spectrum is determined. The agreement between theory and Raman spectroscopy data for B-DNA is demonstrated. Conformational dependences of vibration spectrum during the B→A and helix→coil DNA transitions are studied. The contribution coming from low-frequency mobility to the nucleic-protein recognition processes is discussed.     

Conformation oscillations of DNA

A conformational mobility model is proposed to describe the low-frequency DNA structure dynamics. The transverse oscillations of structure elements in a double-stranded chain are studied. The frequency of DNA long-wave conformational oscillations are estimated and the theory is shown to agree with experiment. DNA low-frequency Raman spectra are interpreted.     

Rate constants and activation parameters for the mobility of bulk and DNA-associated glycol-water solvents

The decrease in mobility of viscous glycol–water solvents when associated with native DNA is quantified from a study of the loss of the exciting-wavelength dependence of the phosphorescence spectrum of free and bound proflavin with increasing temperature. The data are interpreted in terms of a distribution of rate constants with an Arrhenius temperature dependence. Over the temperature range of the experiments a relative decrease of ～ 10⁴ in the average rate constant is observed for reorientation of the solvent when associated with DNA. The basis for this large reduction is found to derive from a large decrease in the pre-exponential factors (i.e., activation entropy) associated with the reorientation rate constants. The changes in the distribution of rate constants and the activation parameters for solvent mobility induced by DNA do not resemble the changes observed for any one of a number of small ion or molecule perturbations. The results suggest the presence of disorganized, relatively immobile solvent in association with DNA.     

Low-frequency vibrations of DNA molecules.

A model for calculating the low-frequency modes in DNA molecules is presented. The present model is associated with the 'breathing' of a DNA molecule as well as its complementary hydrogen bonds. The calculated results show excellent agreement with the observed low-frequency wavenumber (30 cm-1). Consequently, such an internal motion as reflected in the proposed model might be the origin of the observed low-frequency vibration in DNA molecules. This is helpful for investigating the relevant biological functions, which so far have been discussed by many scientists.     

Active miniature transposons from a plant genome and its nonrecombining Y chromosome

Mechanisms involved in eroding fitness of evolving Y chromosomes have been the focus of much theoretical and empirical work. Evolving Y chromosomes are expected to accumulate transposable elements (TEs), but it is not known whether such accumulation contributes to their genetic degeneration. Among TEs, miniature inverted-repeat transposable elements are nonautonomous DNA transposons, often inserted in introns and untranslated regions of genes. Thus, if they invade Y-linked genes and selection against their insertion is ineffective, they could contribute to genetic degeneration of evolving Y chromosomes. Here, we examine the population dynamics of active MITEs in the young Y chromosomes of the plant Silene latifolia and compare their distribution with those in recombining genomic regions. To isolate active MITEs, we developed a straightforward approach on the basis of the assumption that recent transposon insertions or excisions create singleton or low-frequency size polymorphisms that can be detected in alleles from natural populations. Transposon display was then used to infer the distribution of MITE insertion frequencies. The overall frequency spectrum showed an excess of singleton and low-frequency insertions, which suggests that these elements are readily removed from recombining chromosomes. In contrast, insertions on the Y chromosomes were present at high frequencies. Their potential contribution to Y degeneration is discussed.     

Interactions of two cytotoxic organotin(IV) compounds with calf thymus DNA

The reactions with DNA of two antitumor active organotin(IV) compounds, the dimer of bis[(di-n-butyl 3,6-dioxaheptanoato)tin] (C(52)H(108)Sn(4)O(1) x 2H(2)O), compound 1, and tri-n-butyltin 3,6,9-trioxodecanoate (C(19)H(40)SnO(5) x 1/2H(2)O), compound 2, were analysed by circular dichroism, DNA melting experiments and gel mobility shift assays. It is found that both complexes modify only slightly the B-type circular dichroism spectroscopy (CD) spectrum of calf thymus DNA. On the other hand, both complexes were found to affect significantly the parameters of the thermally induced helix-to-coil transition. Addition of 1 or 2 to calf thymus DNA samples does not favor DNA renaturation after melting ruling out formation of interstrand crosslinks. Moreover, the effects of both compounds on plasmid DNA gel mobility were investigated. From the analysis of the present results it is inferred that both organotin(IV) compounds do interact with DNA, probably at the level of the phosphate groups.     

Gel electrophoretic technique for separating crosslinked RNAs. Application to improved electron microscopic analysis of psoralen crosslinked 16 S ribosomal RNA

We have developed a gel electrophoresis technique for separating crosslinked RNA molecules into a series of discrete fractions. The gel used is polyacrylamide made in formamide and low salt designed to denature the RNA during electrophoresis. The mobility depends upon the position of crosslinking within each molecule, as demonstrated by electron microscopy of RNA eluted from the gel. In general, molecules with large loops electrophorese more slowly than molecules with small loops or uncrosslinked molecules. We have used this technique to re-examine the psoralen crosslinking pattern of Escherichia coli 16 S ribosomal RNA in inactivated 30 S ribosomal subunits. To determine the correct orientation of each type of crosslink, we have covalently attached DNA restriction fragments to the RNA so that the polarity of the RNA in the microscope would be known. Our previous major conclusions are confirmed: the predominant long-distance crosslink detected by gel electrophoresis involves a residue close to the 3′ end and a residue approximately 600 nucleotides away: the formamide/polyacrylamide gel is able to separate two closely spaced 1100-nucleotide interactions beginning close to the 3′ end, which were reported as one interaction before: and an interaction joining the ends is detected as before. However, one low-frequency crosslinked interaction, between positions 950 and 1400, and possibly another low-frequency interaction, between positions 550 and 870, are determined to be in the opposite polarity to that described previously.     

Structure and stability of complexes formed by nucleic acids. I. Binding of acridine to DNA

Nuclear magnetic resonance spectra of acridine have been measured in aqueous methanol solutions over a wide concentration range in the presence and absence of dissolved DNA. In solutions containing DNA the acridine spectra show a marked line broadening and intensity decrease at temperatures lower than 50°C. These line-shape changes can be associated with two types of binding interactions: (1) a tight, irrotational binding of the acridine at low acridine:phosphate ratios and (2) a weaker, rotationally less restrictive binding at high acridine concentrations. At temperatures above 50°C. a marked line narrowing is noted for the acridine spectrum and is attributed to an increase in mobility of the bound acridine as the DNA complex undergoes a helix–coil transition. A loose association of acridine molecules with the purine and pyrimidine bases in heat-denatured DNA is indicated by chemical shift changes in the acridine spectrum. The NMR measurements also show that the presence of acridine in denatured DNA solutions greatly reduces renaturation of the DNA.     

Studies with a fluorescent vital probe for DNA in mammalian cells

Olivomycin is taken up efficiently by HeLa cells and by rat fibroblast cells at 38.5 °C, but not by BHK cells. On irradiation with light of 425 nm wavelength, the nuclei of living cells that have taken up olivomycin fluoresce. When olivomycin complexes with DNA in solution, the emission spectrum broadens and shifts, the excitation wavelength maximum shifts up 15 nm, and the fluorescence polarization increases. In HeLa and fibroblast cells, the fluorescence characteristics indicate that olivomycin is entirely complexed to DNA, and its rotational mobility indicates that it is complexed to DNA in regions where other components of the chromatin offer no steric hindrance.     

Quasi-continuum models of twist-like and accordion-like low-frequency motions in DNA.

Formulae for calculating low-frequency twist-like and accordion-like modes of DNA molecules have been derived using a quasi-continuum model. The formulae can be employed in essentially all (viz. A, B, C, D, E, and Z) forms of DNA. Calculated results indicate that the experimentally observed low-frequency modes at 22 cm-1 for the A-form octanucleotide (d[CCCCGGGG]) and at 18 cm-1 for the B-form dodecanucleotide (d[CGCAA ATTTGCG]) may result from accordion-like motions, while those observed at 12 cm-1 and 15 cm-1 may result from combinations of twist-like oscillations excited in the intact segments of B- and A-DNA's, respectively. Frequency shifts in the low-frequency modes observed when DNA molecules undergo conformational changes among different forms are also discussed in terms of the current model.     

Normal vibrations of -glycylglycine

A complete normal coordinate analysis, infrared absorption, and inelastic neutron scattering studies on α-glycylglycine are reported. A Urey-Bradley force field which is also valid for N-deuterated sample is obtained. Assignments in the low-frequency region are based on inelastic neutron spectrum obtained in one phonon and cubic approximation.     

Low-frequency fourier transform infrared spectroscopy of the oxygen-evolving and quinone acceptor complexes in photosystem II

The low-frequency (<1000 cm-1) region of the IR spectrum has the potential to provide detailed structural and mechanistic insight into the photosystem II/oxygen evolving complex (PSII/OEC). A cluster of four manganese ions forms the core of the OEC and diagnostic manganese-ligand and manganese-substrate modes are expected to occur in the 200-900 cm-1 range. However, water also absorbs IR strongly in this region, which has limited previous Fourier transform infrared (FTIR) spectroscopic studies of the OEC to higher frequencies (>1000 cm-1). We have overcome the technical obstacles that have blocked FTIR access to low-frequency substrate, cofactor, and protein vibrational modes by using partially dehydrated samples, appropriate window materials, a wide-range MCT detector, a novel band-pass filter, and a closely regulated temperature control system. With this design, we studied PSII/OEC samples that were prepared by brief illumination of O2 evolving and Tris-washed preparations at 200 K or by a single saturating laser flash applied to O2 evolving and inhibited samples at 250 K. These protocols allowed us to isolate low-frequency modes that are specific to the QA-/QA and S2/S1 states. The high-frequency FTIR spectra recorded for these samples and parallel EPR experiments confirmed the states accessed by the trapping procedures we used. In the S2/S1 spectrum, we detect positive bands at 631 and 602 cm-1 and negative bands at 850, 679, 664, and 650 cm-1 that are specifically associated with these two S states. The possible origins of these IR bands are discussed. For the low-frequency QA-/QA difference spectrum, several modes can be assigned to ring stretching and bending modes from the neutral and anion radical states of the quinone acceptor. These results provide insight into the PSII/OEC and demonstrate the utility of FTIR techniques in accessing low-frequency modes in proteins.     

Diagonalization in a mixed basis: A method to compute low-frequency normal modes for large macromolecules

Low-frequency collective motions in proteins are generally very important for their biological functions. To study such motions, harmonic dynamics proved most useful since it is a straightforward method; it consists of the diagonalization of the Hessian matrix of the potential energy, yielding the vibrational spectrum and the directions of internal motions. Unfortunately, the diagonalization of this matrix requires a large computer memory, which is a limiting factor when the protein contains several thousand atoms. To circumvent this limitation we have developed three methods that enable us to diagonalize large matrices using much less computer memory than the usual harmonic dynamics. The first method is approximate; it consists of diagonalizing small blocks of the Hessian matrix, followed by the coupling of the low-frequency modes obtained for each block. It yields the low-frequency vibrational spectrum with a maximum error of 20%. The second method consists, after diagonalizing small blocks, of coupling the high- and low-frequency modes using an iterative procedure. It yields the exact low-frequency normal modes, but requires a long computational time with convergence problems. The third method, DIMB (Diagonalization in a Mixed Basis), which has the best performance, consists of coupling the approximate low-frequency modes with the mass-weighted cartesian coordinates, also using an iterative procedure. It reduces significantly the required computer memory and converges rapidly. The eigenvalues and eigenvectors obtained by this method are without significant error in the chosen frequency range. Moreover, it is a general method applicable to any problem of diagonalization of a large matrix. We report the application of these methods to a deca-alanine helix, trypsin inhibitor, a neurotoxin, and lysozyme. © 1993 John Wiley & Sons, Inc.     

Resonance Raman studies of lactoperoxidase

Resonance Raman (RR) spectra obtained at three excitation wavelengths are reported for various ferric, ferrous, and ferryl derivatives of bovine lactoperoxidase. The RR spectra of the ferric derivatives show the full complement of the vinyl stretching and scissor modes indicating that the two vinyls in the protoporphyrin IX prosthetic group are present in unmodified forms. The cysteine thiol complex exhibits a RR spectrum identical to that of the native enzyme, an observation which strongly suggests a nonheme binding site for the thiol substrates. The different ferrous complexes of lactoperoxidase which result from heme reduction at slightly alkaline and acidic pH gave identical low-frequency RR spectra. Differences are observed, however, in the high-frequency region. Reduction in the presence of cyanide, however, yields two time-resolved complexes. Changes in the ligand field during the conversion to the final form of the cyanoferrous complex are proposed based on the changes observed in the low-frequency vibrational spectrum. Comparisons are made between the low-frequency RR spectra of the limiting form of the cyanoferrous and the nitric oxide lactoperoxidase complexes. The similarity between the RR spectra of these two complexes in the 150-500 cm-1 region supports the assignment of structures for these complexes where the six-coordinate heme iron is displaced from the heme plane and away from the proximal histidine ligand.     

Effects of sequence selective drugs on the gel mobility of a bent DNA fragment.

The effects of various drugs on the structure of a bent DNA fragment have been investigated by studying DNA mobility in polyacrylamide gels. This DNA fragment has an anomalously slow rate of migration on account of its phased runs of adenines. Nogalamycin and echinomycin increase the gel mobility of kinetoplast DNA suggesting that the bending has been removed. Mithramycin, actinomycin, distamycin and ethidium have either no effect or cause a further reduction in mobility. These results are compared with other, non-bent DNA species which always show a decrease in gel mobility in the presence of DNA binding drugs.     

Topological effects on the electrophoretic mobility of rigid rodlike DNA in polyacrylamide gels

The electrophoretic migration of rigid rodlike DNA structures with well defined topologies has been investigated in polyacrylamide (PA) hydrogels prepared by copolymerization of acrylamide and N, N'-methylenebisacrylamide. Previous studies have reported structural and dynamic characteristics of linear and branched DNA during electrophoresis in PA gels using a variety of experimental parameters. However, a thorough investigation aimed at establishing specific relationships between topological features of rigid rodlike DNA structures and their electrophoretic behavior is still needed. In order to study these topological effects on mobility, an intensive examination of the electrophoretic mobility of small linear and starlike DNA was performed. A series of model DNA structures with well-defined branched topologies were synthesized with varying molecular parameters, such as number of arms surrounding the branch point and arm length. The electrophoretic mobility of these structures was then contrasted with a series of data obtained using linear DNA of comparable molecular size. When large DNA stars (M >/= 60 bp) were compared with linear DNA of identical molecular weight, the Ferguson plots were quite different. However, small DNA stars (24-32 bp) and linear analogues had identical Ferguson plots. This indicates that a different motional mode or greater interaction with the gel exists for the larger DNA stars. When the total molecular weight of the DNA stars was held constant and the number of arms varied, the Ferguson plots for all the stars were identical. Additionally, a critical pore size was reached when the ratio of linear DNA mobility to star DNA mobility increased dramatically. Thus, while the incorporation of a single branch point can produce a large reduction in mobility, above a critical molecular size, the incorporation of additional branch points does not appear to provide further reduction in mobility. This finding is consistent with the transport properties of large synthetic star polymers, where a large reduction in their diffusion coefficient is observed when a single branch is added. When additional arms are incorporated, large synthetic stars do not display an appreciable further reduction in diffusion coefficient. The effect of arm length on mobility for rigid rod DNA stars was also studied. For four-arm DNA stars, the mobility was found to scale as an exponential function of the arm length. Finally, a recently proposed phenomenological model was used to successfully fit the mobility data for linear rigid rod DNA at various concentrations of PA.