Interpretation of DNA Vibration Modes: III - The Behaviour of the Sugar Pucker Vibration Modes as a Function of its Pseudorotation Parameters

Abstract

A systematic study of the sugar pucker characteristic vibration modes as a function of its geometrical conformations, has been performed. The present investigation is based on the Wilson GF method and a non-redundant valence force field. The calculated results allow to assign the modes arising mainly from the sugar motions and present in quasi whole vibrational spectra related to the right or left-handed double-helices (i.e., 1050 cm^?1,960 cm^?1 and 890 cm^?1). Moreover, the conformation dependent modes as those at 860 cm^?1 and around 810 cm^?1 (A form) as well as the one located around 830 cm^?1 (B form) are interpreted by the present investigation. The possibility of the interaction of the latter modes with the phosphate group motions along the DNA double-helical chains are also discussed.     

Normal modes of vibration in bovine pancreatic trypsin inhibitor and its mechanical property

The normal mode analysis of conformational fluctuation is carried out for a small globular protein, bovine pancreatic trypsin inhibitor. Results are analyzed mainly to reveal the mechanical construction of the protein molecule. We take dihedral angles, including peptide omega angles, as independent variables for the normal mode analysis. There are 306 such angles in this molecule. Motions in modes with frequencies lower than 120 cm-1 are shown to involve atoms in the whole protein molecule, and spatial change of displacement vectors is continuous, i.e., those of atoms near in space are similar. To quantitate the observation of the continuity, a correlation function of direction vectors of atomic displacements is calculated. From this function we define a quantity that is interpreted as the wave length of an equivalent elastic plane wave. From this quantity we deduce effective Young's modulus for each mode. For the mode with the lowest frequency 4.4 cm-1, it turned out to be 0.8 x 10(9) dyn cm-2, the value two orders of magnitude softer than, for instance, alpha-helices. Prompted by this observation, the four lowest frequency modes and also the harmonic motions in the thermal equilibrium are analyzed further mainly to detect relatively rigid structural elements in the molecule. From this analysis emerges a mechanical picture of the protein molecule that is made up of relatively rigid elements held together by very soft parts.     

Interpretation of DNA vibration modes: I--The guanosine and cytidine residues involved in poly(dG-dC).poly(dG-dC) and d(CG)3.d(CG)3

A normal coordinate analysis has been carried out on guanosine and cytidine residues appearing in oligo and polynucleotides by using a simplified valence force field that allows the vibrational spectra of 5'-dGMP and 2'-deoxycytidine molecules to be reproduced. The role of both C2'-endo and C3'-endo conformations on sugar pucker, as well as that of glycosidic torsion angle (X), on several characteristic vibration modes of these residues have been studied. The present calculations based on a non-redundant set of internal coordinates preserving the harmonic approximation of the potential field, allows us to explain quite satisfactorily the modifications of the vibrational spectra in the 1550-1250 cm-1 and 785-500 cm-1 regions, when the right----left-handed conformational transition occurs.     

Interpretation of DNA vibration modes: IV--A single-helical approach to assign the phosphate-backbone contribution to the vibrational spectra in A and B conformations

A calculated approach based on the Higgs method for assigning the vibration modes of an infinite helicoidal polymeric chain has been performed on the basis of a reliable valence force field. The calculated results allowed the phosphate-backbone marker modes of the A and B forms, to be interpreted. In the dynamic models used, the bases have been omitted and no interchain interaction was considered. The calculation can also interprete quite satisfactorily the characteristic Raman peaks and infrared bands in the 1250-700 cm-1 spectral region arising from the sugar or sugar-phosphate association and reproduce their evolution upon the B----A DNA conformational transition. They clearly show that the phosphate-backbone modes in the above mentioned spectral region constitute the optical branches of the phonon dispersion curves with no detectable variation in the first Brillouin-zone.     

Interpretation of DNA vibration modes. II--The adenosine and thymidine residues involved in oligonucleotides and polynucleotides

Normal coordinate analysis of the adenosine and thymidine residues involved in the right- and left-handed conformations of oligonucleotides and polynucleotides has been performed. The valence force field, employed in this work, allowed recently to reproduce the vibrational spectra of 2'-deoxythymidine and 2'-deoxyadenosine. The calculated wavenumbers based on a non-redundant set of internal coordinates have been compared to the Raman and infrared peak positions arising from A, B, C, D and Z conformations, in the 1550-1250 cm-1 and 800-600 cm-1 spectral regions: i.e. characteristic of adenosine and thymidine residues. Moreover, a systematic study has been performed on the evolution of the vibrational wavenumbers as a function of the glycosidic angle (chi) and the sugar pucker conformation.     

Identification of low-frequency modes in protein molecules.

It is demonstrated that the observed low-frequency motions with wave numbers of 22 cm-1 and 25 cm-1 for insulin and lysozyme respectively originate from the accordion-like motions of the principal helices therein. The calculated results based on such a model are in good agreement with the observed values. During calculations the role of the internal microenvironment upon the low-frequency motion is naturally revealed, so as to elucidate as well why this kind of low-frequency motion is so sensitive to the conformations of proteins observed.     

Resonant and localized breathing modes in terminal regions of the DNA double helix.

A Green's function approach is used in constructing a dynamic model of a semi-infinite length of the DNA homopolymer B poly(d) . poly(d). Considerable attention is focused on the hydrogen bond stretching close to the terminus. A melting (or breathing) coordinate (M) is defined as an average over the three linking hydrogen bond stretches in a unit cell. The thermal mean squared amplitude of (M) is enhanced at the chain end compared with the interior. Spectral branches at 69, 80 and 105 cm-1, as well as a local mode at 75 cm-1, are primary contributors to the enhancement. We suggest that this fact can affect the thermal melting of a DNA double helical homopolymer, enhancing the tendency to start from an end (if one is available). We show how certain infinite chain modes with small (M) amplitude can turn into breathing modes near the terminus, and suggest that the same phenomenon may occur near other specific base-pair sequences. There is also considerable attention paid to the low microwave region from approximately 0 to 1.75 cm-1. The thermally activated modes in this frequency region contribute approximately (0.02 A)2 to [M2(0)] at 40 K, approximately two orders of magnitude greater than for [M2(infinity)]. Most important however, is the existence of narrow resonant modes in this frequency region. Particularly pronounced resonances near 0.03 cm-1 and 0.08 cm-1 (approximately 0.9 and 2.4 GHz) amplify M2(0) at the terminus by about for orders of magnitude over the infinite chain value M2(infinity).     

Resonance Raman enhancement of phenyl ring vibrational modes in phenyl iron complex of myoglobin.

Resonance Raman spectra are reported for the organometallic phenyl-FeIII complexes of horse heart myoglobin. We observed the resonance enhancement of the ring vibrational modes of the bound phenyl group. They were identified at 642, 996, 1,009, and 1,048 cm-1, which shift to 619, 961, 972, and 1,030 cm-1, respectively, upon phenyl 13C substitution. The lines at 642 and 996 cm-1 are assigned, respectively, as in-plane phenyl ring deformation mode (derived from benzene vibration No. 6a at 606 cm-1) and out-of-plane CH deformation (derived from benzene vibration No. 5 at 995 cm-1). The frequencies of the ring "breathing" modes at 1,009 and 1,048 cm-1 are higher than the corresponding ones in phenylalanine (at 1,004 and 1,033 cm-1) and benzene (at 992 and 1,010 cm-1), indicating that the ring C--C bonds are strengthened (or shortened) when coordinated to the heme iron. The excitation profiles of these phenyl ring modes and a porphyrin ring vibrational mode at 674 cm-1 exhibit peaks near its Soret absorption maximum at 431 nm. This appears to indicate that these phenyl ring modes may be enhanced via resonance with the Soret pi-pi transition. The FeIII--C bond stretching vibration has not been detected with excitation wavelengths in the 406.7-457.9-nm region.     

Dynamics of transfer RNAs analyzed by normal mode calculation.

Normal mode calculation is applied to tRNAPhe and tRNAAsp, and their structural and vibrational aspects are analyzed. Dihedral angles along the phosphate-ribose backbone (alpha, beta, gamma, epsilon, zeta) and dihedral angles of glycosyl bonds (chi) are selected as movable parameters. The calculated displacement of each atom agrees with experimental data. In modes with frequencies higher than 130 cm-1, the motions are localized around each stem and the elbow region of the L-shape. On the other hand, collective motions such as bending or twisting of arms are seen in modes with lower frequencies. Hinge axes and bend angles are calculated without prior knowledge. Movements in modes with very low frequencies are combinations of hinge bending motions with various hinge axes and bend angles. The thermal fluctuations of dihedral angles well reflect the structural characters of transfer RNAs. There are some dihedral angles of nucleotides located around the elbow region of L-shape, which fluctuate about five to six times more than the average value. Nucleotides in the position seem to be influential in the dynamics of the entire structure. The normal mode calculation seems to provide much information for the study of conformational changes of transfer RNAs induced by aminoacyl-tRNA synthetase or codon during molecular recognition.     

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.     

Low-frequency dynamics and Raman scattering of crystals, of B-, A-, and Z-DNA, and fibers of C-DNA

Normal modes of vibration of DNA in the low-frequency region (10-300 cm-1 interval) have been identified from Raman spectra of crystals of B-DNA [d(CGCAAATTTGCG)], A-DNA [r(GCG)d(CGC) and d(CCCCGGGG)], and Z-DNA [d(CGCGCG) and d(CGCGTG)]. The lowest vibrational frequencies detected in the canonical DNA structures--at 18 +/- 2 cm-1 in the B-DNA crystal, near 24 +/- 2 cm-1 in A-DNA crystals, and near 30 +/- 2 cm-1 in Z-DNA crystals--are shown to correlate well with the degree of DNA hydration in the crystal structures, as well as with the level of hydration in calf thymus DNA fibers. These findings support the assignment [H. Urabe et al. (1985) J. Chem. Phys. 82, 531-535; C. Demarco et al. (1985) Biopolymers 24, 2035-2040] of the lowest frequency Raman band of each DNA to a helix mode, which is dependent primarily upon the degree of helix hydration, rather than upon the intrahelical conformation. The present results show also that B-, A-, C-, and Z-DNA structures can be distinguished from one another on the basis of their characteristic Raman intensity profiles in the region of 40-140 cm-1, even though all structures display two rather similar and complex bands centered within the intervals of 66-72 and 90-120 cm-1. The similarity of Raman frequencies for B-, A-, C-, and Z-DNA suggests that these modes originate from concerted motions of the bases (librations), which are not strongly dependent upon helix backbone geometry or handedness. Correlation of the Raman frequencies and intensities with the DNA base compositions suggests that the complex band near 90-120 cm-1 in all double-helix structures is due to in-plane librational motions of the bases, which involve stretching of the purine-pyrimidine hydrogen bonds. This would explain the centering of the band at higher frequencies in structures containing G.C pairs (greater than 100 cm-1) than in structures containing A.T pairs (less than 100 cm-1), consistent with the strengths of G.C and A.T hydrogen bonding.     

A Raman study of low frequency intrahelical modes in A-, B-, and C-DNA.

We have obtained low frequency (less than 200 cm-1) Raman spectra of calf-thymus DNA and poly(rI).poly(rC) as a function of water content and counterion species and of d(GGTATACC)2 and d(CGCGAATTCGCG)2 crystals. We have found that the Raman scattering from water in the first and second hydration shells does not contribute directly to the Raman spectra of DNA. We have determined the number of strong Raman active modes by comparing spectra for different sample orientations and polarizations and by obtaining fits to the spectra. We have found at least five Raman active modes in the spectra of A- and B-DNA. The frequencies of the modes above 40 cm-1 do not vary with counterion species, and there are only relatively small changes upon hydration. These modes are, therefore, almost completely internal. The mode near 34 cm-1 in A-DNA is mostly internal, whereas the mode near 25 cm-1 is dominated by interhelical interactions. The observed intensity changes upon dehydration were found to be due to the decrease in interhelical distance. Polymer length appears to play a role in the lowest frequency modes.     

Resonance Raman studies of Escherichia coli sulfite reductase hemoprotein. 3. Bound ligand vibrational modes

The vibrations of the bound diatomic heme ligands CO, CN-, and NO are investigated by resonance Raman spectroscopy in various redox states of Escherichia coli sulfite reductase hemoprotein, and assignments are generated by use of isotopically labeled ligands. For the fully reduced CO complex (ferrous siroheme, reduced Fe4S4 cluster) at room temperature, nu CO is observed at 1904 cm-1, shifting to 1920 cm-1 upon oxidation of the cluster. The corresponding delta FeCO modes are identified at 574 and 566 cm-1, respectively, by virtue of the zigzag pattern of their isotopic shifts. In frozen solution, two species are observed for the cluster-oxidized state, with nu CO at 1910 and 1936 cm-1 and nu FeC at 532 and 504 cm-1, respectively; nu FeC for the fully reduced species is identified at 526 cm-1 in the frozen state. For the ferrous siroheme-NO complex (cluster oxidized), nu NO is identified at 1555 cm-1 in frozen solution and a low-frequency mode is identified at 558 cm-1; this stretching mode is significantly lower than that observed in Mb-NO. For the ferric siroheme cyanide complexes evidence of two ligand-bonding forms is observed, with modes at 451/390 and 451/352 cm-1; they are distinguished by a reversal of the isotopic shift patterns of the upper and lower modes and could arise from a linear and a bent Fe-C unit, respectively. For the ferrous siroheme cyanide complex isotope-sensitive modes observed at 495 and 452 cm-1 are assigned to the FeCN- bending and FeC stretching vibrations, respectively.(ABSTRACT TRUNCATED AT 250 WORDS)     

Stability of triple-helical poly(dT)-poly(dA)-poly(dT) DNA with counterions.

Structural conformation of triple-helical poly(dT)-poly(dA)-poly(dT) has been a very controversial issue recently. Earlier investigations, based on fiber diffraction data and molecular modeling, indicated an A-form conformation with C'3-endo sugar pucker. On the other hand, Raman, solution infrared spectral, and NMR studies show a B-form structure with C'2-endo sugars. In accordance with these experimental results, a theoretical model with B-form, C'2-endo sugars was proposed in 1993. In the present work we investigate the dynamics and stability of the two conformations within the effective local field approach applied to the normal mode calculations for the system. The presence of counterions was explicitly taken into account. Stable equilibrium positions for the counterions were calculated by analyzing the normal mode dynamics and free energy of the system. The breathing modes of the triple helix are shifted to higher frequencies over those of the double helix by 4-16 cm-1. The characteristic marker band for the B conformation at 835 cm-1 is split up into two marker bands at 830 and 835 cm-1. A detailed comparison of the normal modes and the free energies indicates that the B-form structure, with C'2-endo sugar pucker, is more stable than the A-form structure. The normal modes and the corresponding dipole moments are found to be in close agreement with recent spectroscopic findings.     

Local modes in a DNA polymer with hydrogen bond defect.

Vibrations of a homopolymer DNA with localized hydrogen bond defects have been examined using the recently developed decaying mode theory for long-chain polymers with local structural defects. For a poly(dA)-poly(dT) homopolymer having perturbed hydrogen bonds in one base pair, a localized mode at 63.2 cm-1 has been found. This mode has a very nearly pure H-bond stretch or "breathing" character, although the backbones do not separate. This agrees in frequency with a similar result found by other authors using a different approach. We search the full microwave frequency range for other local modes for several models of weakened H bonds. Besides the local mode with breathing characteristics, local modes with other characteristic motions were found, but only for asymmetrically perturbed bonds. We find in general that local modes are not very robust, requiring quite specific, narrow ranges in parameter space. They are also not abundant, there being only three in our most prolific model.     

Breathing mode of conformational fluctuations in globular proteins.

Globular proteins in the native state are assumed to behave as continuous elastic spheres in the low frequency breathing motions. Reasonable values of Young's modulus E = 10(11) dyne/cm2 and the radius of the sphere ro = 20 A, yield a wave number of 26 cm-1 for the fundamental vibration of the sphere. The peak at around 30 cm-1 in the laser Raman spectra of native alpha-chymotrypsin and pepsin observed by Brown et al. might be assigned to the breathing motion which the native proteins undergo as continuous elastic bodies.     

Terahertz vibration modes in Na/K-ATPase

Mechanical vibration in the Terahertz range is believed to be connected with protein functions. In this paper, we present the results of a normal-mode analysis (modal analysis) of a Na/K-ATPase all-atom model, focusing the attention on low-frequency vibration modes. The numerical model helps in the interpretation of experimental results previously obtained by the authors via Raman spectroscopy of Na/K-ATPase samples, where several unassigned peaks were found in the sub-500 cm^?1 range. In particular, vibration modes corresponding to peaks at 27, 190 and 300 cm^?1, found experimentally, are confirmed here numerically, together with some other modes at lower frequencies (wavenumbers) that were not possible to observe in the experimental test. All the aforementioned modes correspond to vibrations involving the protein ends, i.e. portions directly related to the operating mechanism of the sodium-potassium pump.     

Investigations of optical line shapes and kinetic hole burning in myoglobin.

We present the results of an extensive investigation of the optical line shapes of deoxymyoglobin (Mb), the ligand-bound form (MbCO), and the low-temperature photoproduct (Mb*). The thermal properties and the pH dependence of the Soret band and the near infrared band III (approximately 760 nm) are analyzed, taking into account the underlying vibrational properties of the absorption bands. The strong temperature dependence associated with the Soret band of MbCO and band III of Mb indicates significant coupling to low-frequency modes that may not be directly observed in the resonance Raman spectra. On the basis of analogous line-shape studies in a variety of heme systems, we assign the low-frequency coupling in MbCO to torsional motions of the CO molecule. The low-frequency mode coupled to band III (approximately 70 cm-1) is found to lie quite close to the value for the heme-doming motion (approximately 50 cm-1) calculated by using the kinetically determined value of the force constant (17 N/m). Significant inhomogeneous broadening in the Soret region of Mb and Mb* is found to be due to a "nonkinetic" coordinate that we associate with the orientation of the proximal histidine. A "kinetic" coordinate, associated with the equilibrium displacement of the iron atom from the porphyrin plane (a) is found to contribute to the inhomogeneous broadening of both the Soret band and band III. The relaxation of the heme as the system evolves from from Mb* to Mb is followed optically as a function of temperature, and a sharp transition temperature is found at 185 K. The blue shifts of the Soret band and band III as Mb* evolves to Mb are found to be nearly identical (delta v*ABS approximately 140 cm-1) and attributed to changes in the mean value of a between Mb* (a*0) and Mb (a0 = 0.45 A). A simple quadratic model for the coordinate coupling that simultaneously accounts for the observed shift, delta v*ABS, the low-temperature kinetics and the kinetic hole burning predicts a*0 = 0.2 +/- 0.05 A and EA = 16 +/- 2 kJ/mol for the room temperature Arrhenius barrier height at the heme. A simple quantitative method for the analysis of kinetic hole-burning experiments is also developed and applied to recent studies involving quaternary and subunit-specific hemoglobin structures.     

Direct measurement of hydration-related dynamic changes in lysozyme using inelastic neutron scattering spectroscopy

Inelastic neutron scattering spectroscopy is used to investigate dynamic changes in lysozyme powder at two different low D2O hydrations (0.07g D2O/g protein and 0.20 g D2O/g protein). In the higher hydration sample, the inelastic scattering between 0.8 and 4.0 cm-1 energy transfer is increased and the elastic scattering is decreased. The decreased elastic scattering suggests increased atomic amplitudes of motion and the increased 0.8 to 4.0 cm-1 scattering suggests increased motions in this frequency range. Comparison with normal mode models of lysozyme dynamics shows that the inelastic difference occurs in the frequency region predicted for the lowest frequency, largest amplitude, global modes of the molecular [M. Levitt, C. Sander and P.S. Stern, J. Mol. Biol. 181, 423 (1985). B. Brooks and M. Karplus, Proc. Natl. Acad. Sci (U.S.A) 82, 4995 (1985), R.E. Bruccoleri, M. Karplus and J.A. McCammon, Biopolymers 25 1767 (1986)]. Our results are consistent with a model in which an increased number of low frequency global modes are present in the higher hydrated sample.     

Resonance Raman studies of the HOOP modes in octopus bathorhodopsin with deuterium-labeled retinal chromophores

Resonance Raman spectra of the hydrogen out-of-plane (HOOP) vibrational modes in the retinal chromophore of octopus bathorhodopsin with deuterium label(s) along the polyene chain have been obtained. In clear contrast with bovine bathorhodopsin's HOOP modes, there are only two major HOOP bands at 887 and 940 cm-1 for octopus bathorhodopsin. On the basis of their isotopic shifts upon deuterium labeling, we have assigned the band at 887 cm-1 to C10H and C14H HOOP modes, and the band at 940 cm-1 to C11H = C12H Au-like HOOP mode. Except for a 26 cm-1 downward shift, the C11H = C12H Au-like wag appears to be little disturbed in octopus bathorhodopsin from the chromophore in solution since its changes upon deuterium labeling are close to those found in solution model-compound studies. We found also that the C10H and C14H HOOP wags are also similar to those in the model-compound studies. However, we have found that the interaction between the C7H and C8H HOOP internal coordinates of the chromophore in octopus bathorhodopsin is different from that of the chromophore in solution. The intensity of the C11H = C12H and the other HOOP modes suggests that the chromophore of octopus bathorhodopsin is somewhat torsionally distorted from a planar trans geometry. Importantly, a twist about C11 = C12 double bond is inferred. Such a twist breaks the local symmetry, resulting in the observation of the normally Raman-forbidden C11H = C12H Au-like HOOP mode. The twisted nature of the chromophore, semiquantitatively discussed here, likely affects the lambda max of the chromophore and its enthalpy.(ABSTRACT TRUNCATED AT 250 WORDS)