Vibrational dynamics of DNA. I. Vibrational basis modes and couplings

Carrying out density functional theory calculations of four DNA bases, base derivatives, Watson-Crick (WC) base pairs, and multiple-layer base pair stacks, we studied vibrational dynamics of delocalized modes with frequency ranging from 1400 to 1800 cm(-1). These modes have been found to be highly sensitive to structure fluctuation and base pair conformation of DNA. By identifying eight fundamental basis modes, it is shown that the normal modes of base pairs and multilayer base pair stacks can be described by linear combinations of these vibrational basis modes. By using the Hessian matrix reconstruction method, vibrational coupling constants between the basis modes are determined for WC base pairs and multilayer systems and are found to be most strongly affected by the hydrogen bonding interaction between bases. It is also found that the propeller twist and buckle motions do not strongly affect vibrational couplings and basis mode frequencies. Numerically simulated IR spectra of guanine-cytosine and adenine-thymine bases pairs as well as of multilayer base pair stacks are presented and described in terms of coupled basis modes. It turns out that, due to the small interlayer base-base vibrational interactions, the IR absorption spectrum of multilayer base pair system does not strongly depend on the number of base pairs.     

Interpreting DNA vibrational circular dichroism spectra using a coupling model from two-dimensional infrared spectroscopy

Two-dimensional infrared spectroscopy was recently used to measure the vibrational couplings between carbonyl bonds located on DNA nucleobases (Krummel, A. T.; Mukherjee, P.; Zanni, M. T. J. Phys. Chem. B 2003, 107, 9165 and Krummel, A. T.; Zanni, M. T. J. Phys. Chem. B 2006, 110, 13991). Here, we extend the coupling model derived from these 2D IR experiments to simulate the vibrational absorption and vibrational circular dichroism (VCD) spectra of three double-stranded DNA oligomers: poly(dG)-poly(dC), poly(dG-dC), and dGGCC. Using this model, we determine that the VCD spectrum of A-form poly(dG)-poly(dC) is dominated by interactions between stacked bases, whereas the coupling between base pairs and stacked bases carries equal importance in the VCD spectrum of B-form poly(dG-dC). We also simulate the absorption and VCD spectra of dGGCC, which is a combination of A- and B-form configurations. These simulations give insight into the structural interpretation of VCD and absorption spectroscopies that have long been used to monitor DNA secondary structure and kinetics.     

核酸碱基的非谐性振动:正规及错配碱基对的一维和二维红外光谱及其结构基础

用量子化学计算研究了正规Watson-Crick碱基对和四例典型的错配碱基对.对碱基单体和二聚体进行了详细的非谐性频率分析,以揭示其结构方面的一些振动特征.研究发现这些振动特征能在模拟的一维和二维红外光谱中很好地表现出来.利用势能分布研究了所选简正模式的离域化程度,发现从孤立的碱基单体到参与氢键的二聚体,模式的离域化程度变化很大;同时,这些模式的非谐性常数也发生了相应的改变.以人们通常认为的位于红外光谱中6-μm波长区域的羰基伸缩模式为例进行了探讨.     

Calculations of intermode coupling constants and simulations of amide I,II, and III vibrational spectra of dipeptides

Amide I, II, and III vibrations of polypeptides are important marker modes whose vibrational spectra can provide critical information on structure and dynamics of proteins in solution. The extent of delocalization and vibrational properties of amide normal mode can be described by the amide local mode frequencies and intermode coupling constants between a pair of amide local modes. To determine these fundamental quantities, the previous Hessian matrix reconstruction method has been generalized here and applied to the density functional theory results for various dipeptide conformers. The calculation results are then used to simulate IR absorption, vibrational circular dichroism, and 2D IR spectra of dipeptides. The relationships between dipeptide backbone conformations and these vibrational spectra are discussed. It is believed that the present computational method and results will be of use to quantitatively simulate vibrational spectra of complicated polypeptides beyond simple dipeptides     

DNA vibrational coupling revealed with two-dimensional infrared spectroscopy: insight into why vibrational spectroscopy is sensitive to DNA structure

Two-dimensional infrared (2D IR) spectroscopy was used to study the carbonyl vibrational modes of guanine and cytosine bases in A- and B-form DNA. Located between 1600 and 1700 cm(-1), these modes are often used to monitor DNA secondary structure with traditional infrared spectroscopies such as FTIR, but traditional spectroscopies lack the necessary observables to unravel the coupling mechanisms that make these modes sensitive to secondary structure. By using 2D IR spectroscopy and electronic structure calculations on d(G(5)C(5)) and d(GC)(8) model nucleic acids, we find that hydrogen-bonded guanine/cytosine base pairs are primarily electrostatically coupled and that the coupling between these modes can be modeled with a transition dipole density approach. In comparison, electrostatics is insufficient to model stacked bases because of cooperative charge-sharing effects, but the coupling can be accurately calculated using a finite difference method. We find that the coupling is very strong for both hydrogen-bonded and stacked base geometries, creating vibrational modes that extend both across the base pairs and along the lengths of the helices. Our results provide a physical basis for understanding how strong coupling gives rise to the empirically established relationship between infrared spectroscopy and DNA/RNA secondary structure.     

Vibrational dynamics of DNA: IV. Vibrational spectroscopic characteristics of A-, B-, and Z-form DNA's

Linear and nonlinear IR spectroscopic studies of nucleic acids can provide crucial information on solution conformations of DNA double helix and its complex with other molecules. Carrying out density functional theory calculations of A-, B-, and Z-form DNA's, the authors obtained vibrational spectroscopic properties as well as coupling constants between different basis modes. The vibrational couplings that determine the extent of exciton delocalization are strongly dependent on DNA conformation mainly because the interlayer distance between two neighboring base pairs changes with respect to the DNA conformation. The Z-DNA has comparatively small interlayer vibrational coupling constants so that its vibrational spectrum depends little on the number of base pairs, whereas the A-DNA shows a notable dependency on the size. Furthermore, it is shown that a few distinctively different line shape changes in both IR and two-dimensional IR spectra as the DNA conformation changes from B to A or from B to Z can be used as marker bands and characteristic features distinguishing different DNA conformations.     

Spectroscopic investigation of organosilicon polymers. I. Vibrational analysis of poly(tetramethyl-p-silphenylene-siloxane)

This article is a detailed, normal coordinate analysis of an isolated chain of poly(tetramethyl-p-silphenylene-siloxane), based on results obtained for the corresponding monomer and a recent x-ray structure. Predicted vibrational frequencies are compared with new experimental infrared (IR) and Raman data obtained from powder and well oriented films. Vibrational spectra were also obtained for a deuterio derivative in which all hydrogens in the benzene rings were replaced with deuterium atoms. The results indicate that an intermode mixing occurs between vibrations of the ring and the saturated part which lead to low frequency normal modes characteristic of the chain. The simplified valence force field obtained for this polymer should provide a useful tool for understanding the vibrational spectra of other members of the class of polysiloxanes.     

Noncovalent interactions between modified cytosine and guanine DNA base pair mimics investigated by terahertz spectroscopy and solid-state density functional theory

Modified cytosine and guanine nucleobases cocrystallize in a hydrogen bonding configuration similar to that observed in native DNA. The noncovalent interactions binding these base pairs in the crystalline solid were investigated using terahertz (THz) spectroscopy and solid-state density functional theory (DFT). While stronger hydrogen bonding interactions are responsible for the general molecular orientations in the crystalline state, it is the weaker dipole-dipole and dispersion forces that determine the overall packing arrangement. The inclusion of dispersion interactions in the DFT calculations was found to be necessary to accurately simulate the unit cell structure and THz vibrational spectrum. Using properly modeled intermolecular potentials, the lattice vibrational motions of the cytosine and guanine derivatives were calculated. The vibrational characters of the modes exhibited by the DNA base pair mimic in the THz region were primarily rotational motions and are indicative of the energies and the nature of vibrations that would likely be observed between similar base pairs in DNA molecules.     

核酸碱基的非谐性振动:Ⅲ,DNA叠加二聚体的一维和二维红外光谱

用量子化学计算方法研究了DNA碱基单体及其15个B型和2个G-四链体的叠加二聚体的非谐性振动光谱特征.研究发现,从碱基单体到二聚体,C=O伸缩模式之间的相互作用很强,表现在其非谐性振动频率和非谐性常数都发生了明显的改变,特别是在含碱基G的叠加体中.这些变化能在一维和二维红外光谱中很好地表现出来.利用振动模式的势能分布和非谐性常数的组成分析,讨论了叠加二聚体中C=O模式的离域化程度.     

Vibrational spectra of nucleic acid bases and their Watson–Crick pair complexes

The vibrational spectra of the nucleic acid bases adenine, thymine, guanine, and cytosine are calculated in the frame of density functional theory (DFT). In particular we use the Kohn–Sham scheme with gradient corrections for exchange and correlation to determine normal modes, frequencies, and intensities. The DFT results are found to be in good agreement with the experiment. Our computations provide assignments for IR, Raman, and neutron inelastic scattering spectroscopies; yield characteristic vibrational fingerprints of each compound for its identification in larger systems; and show general vibrational trends of nucleic acids. The Kohn–Sham scheme is further applied to obtain the spectra of the Watson–Crick pairs adenine-thymine and guanine-cytosine. A large number of monomeric vibrations are recognized in dimers; characteristic vibrations of pairs, which are mainly attributed to hydrogen bridges, are quantified according to changes in normal modes and frequency shifts. Binding and zero-point vibrational energies are analyzed to establish the stability of the complexes and discuss the quality of the energetic calculations. ©1999 John Wiley & Sons, Inc. J Comput Chem 20: 511–530, 1999     

Anharmonic Vibrational Signatures of DNA Bases and Watson-Crick Base Pairs

Changes of molecular structure and associated charge distributions, and changes of anharmonic vibrational parameters from DNA base monomers to the Watson-Crick base pairs, have been investigated at the density functional theory level. Through examination of the NH2, N H, and C=O stretching vibrational modes that are involved in the multiple H-bonds in the base pairs, sensitivity of their diagonal and off-diagonal anharmonicities, as well as anharmonic vibrational couplings, to the structure change are predicted. Our results reveal the intrinsic connection between the anharmonic vibrational potentials, H-bonding, and electrostatic interactions in DNA bases.     

Spectroscopic study of the interaction of pazelliptine with nucleic acids

The antitumor drug pazelliptine (PZE) binds to natural and synthetic DNA sequences at 100 mM NaCl, pH 7.0, as deduced from the absorption and fluorescence data. Scatchard plots constructed from the results obtained with poly(dG-dC)-poly(dG-dC) give binding constants of base pairs in the range (2–6) × 10⁵ M⁻¹. The modifications in the absorption and fluorescence spectra observed when PZE binds to various polynucleotides, namely poly(dA-dT)-poly(dA-dT), poly(dA)-poly(dT), poly(dG-dC)-poly(dG-dC) and calf thymus DNA. reveal a change in the protonation state of the drug upon binding, increasing the apparent pK_a of its 9-N nitrogen atom. The PZE excited state properties serve as a sensitive probe to distinguish between homo and hetero A-T sites as well as between AT and GC sites. Fluorescence studies reveal that energy transfer occurs from polynucleotide bases to the bound PZE chromophore, a result consistent with an intercalative mode of binding of the drug to DNA. The emission is enhanced when PZE is bound to A-T base pairs ( 30% increase of φ_F) whereas it is quenched in the vicinity of G-C base pairs ( 90% decrease of φ_F). Furthermore, the fluorescence spectrum obtained with calf thymus DNA is hardly distinguishable from that obtained with poly(dG-dC)-polu(dG-dC), suggesting a binding of PZE to G-C rich regions.     

Infrared Spectroscopy Evidence of Weak Interactions in Frustrated Lewis Pairs Formed by Tris(pentafluorophenyl)borane

Frustrated Lewis pairs (FLPs) have been widely investigated as promising catalysts due to their metal-free feature and ability to activate small molecules. Since their discovery, many works have been investigating how these Lewis pairs (intermolecular pairs) are held together in an encounter complex. This prompted several studies based on theoretical investigations, but experimental ones are limited yet. In this communication we show evidence of weak intermolecular interactions between Lewis acids and Lewis bases, distinguishing the Lewis adduct from FLPs, by probing fluorine-carbon vibrational modes using infrared spectroscopy. The main evidence is based on the band shifts occurring in FLPs due to weak hydrogen bonds between the hydrogen atoms of the Lewis base and the fluorine atoms of Lewis acid.     

Photoacoustic spectra and modes of vibration of TNT and RDX at CO2 laser wavelengths

IR spectra in the 'Finger Print' spectral range has great importance in qualitative and quantitative analysis of explosives like trinitrotoluene (TNT) and cyclotrimethyltrinitramine (RDX). Highly resolved IR bands of these compounds have been recorded in the 9.6 and 10.6 microm regions of CO2 laser. TNT and RDX are large molecules each having 21 atoms and it is very difficult to assign the modes of vibrations by comparison with those in other molecules making the vibrational assignments of observed bands a difficult task. The ab initio quantum chemical calculation is used for determining the molecular geometries and modes of vibration of these molecules with a view to assign their normal modes in the high resolution vibrational photoacoustic spectra. These assignments are very reliable in view of the good agreement between the observed and calculated frequencies of deuterated TNT.     

Ground- and excited-state properties of DNA base molecules from plane-wave calculations using ultrasoft pseudopotentials

We present equilibrium geometries, vibrational modes, dipole moments, ionization energies, electron affinities, and optical absorption spectra of the DNA base molecules adenine, thymine, guanine, and cytosine calculated from first principles. The comparison of our results with experimental data and results obtained by using quantum chemistry methods show that in specific cases gradient-corrected density-functional theory (DFT-GGA) calculations using ultrasoft pseudopotentials and a plane-wave basis may be a numerically efficient and accurate alternative to methods employing localized orbitals for the expansion of the electron wave functions.     

Terahertz absorption of DNA decamer duplex

This work combines experimental and theoretical approaches to investigate terahertz absorption spectra of the DNA formed by the sequence oligomer 5'-CCGGCGCCGG-3'. The three-dimensional structure of this self-complimentary DNA decamer has been well-studied, permitting us to perform direct identification of the low-frequency phonon modes associated with specific conformation and to conduct comprehensive computer simulations. Two modeling techniques, normal-mode analysis and nanosecond molecular dynamics with explicit solvent molecules, were employed to extract the low-frequency vibrational modes based on which the absorption spectra were calculated. The absorption spectra of the DNA decamer in aqueous solution were measured in the frequency range 10-25 cm(-1) using the terahertz Fourier transform infrared spectroscopy. Multiple well-resolved and reproducible resonance modes were observed. When calculated and experimental spectra were compared, the spectrum based on molecular dynamics simulations showed a better correlation with the experimental spectra than the one based on normal-mode analysis. These results demonstrate that there exist a considerable number of active low-frequency phonon modes in this short DNA duplex.     

Strong enhancement of vibrational relaxation by Watson-Crick base pairing

We have studied the ultrafast dynamics of NH-stretch vibrational excitations in Watson-Crick base pairs consisting of adenine and uracil derivatives. To estimate the influence of the A:U hydrogen bonding on the vibrational dynamics, we have also studied the uracil derivative in monomeric form. The vibrational relaxation of the NH-stretching mode is found to occur much faster in the Watson-Crick base pair than in monomeric uracil. From the delay dependence of the transient vibrational spectra, it can be concluded that both in base-paired and monomeric uracil, the energy relaxation takes place in two steps, the first step being a rapid transfer of energy from the NH-stretching mode to an accepting mode, the second step the relaxation of this accepting mode. The transient spectra show evidence that in the base pair the hydrogen bond between the nucleobases acts as the accepting mode, and that the hydrogen bonding between the bases is responsible for the extremely fast vibrational relaxation in this system.     

Infrared study on the interaction of poly(vinyl chloride) and ketones. I. PVC-MEK film system

In a previous article we presented preliminary information on IR spectrometric evidence for the interaction between poly(vinyl chloride) and ketones by using films with residual solvent as samples for spectrometry. This article reports detailed information on PVC–MEK film systems as compared to model compounds. The results indicate that there are “specific” interactions which were observed by a frequency shift of the carbonyl absorption of the ketone and changes on intensities of the C? Cl vibrational modes of the polymer, which are similar to the effects observed for nonpolymeric systems.     

Anharmonic vibrational modes of nucleic acid bases revealed by 2D IR spectroscopy

Polarization-dependent two-dimensional infrared (2D IR) spectra of the purine and pyrimadine base vibrations of five nucleotide monophosphates (NMPs) were acquired in D(2)O at neutral pH in the frequency range 1500-1700 cm(-1). The distinctive cross-peaks between the ring deformations and carbonyl stretches of NMPs indicate that these vibrational modes are highly coupled, in contrast with the traditional peak assignment, which is based on a simple local mode picture such as C═O, C═N, and C═C double bond stretches. A model of multiple anharmonically coupled oscillators was employed to characterize the transition energies, vibrational anharmonicities and couplings, and transition dipole strengths and orientations. No simple or intuitive structural correlations are found to readily assign the spectral features, except in the case of guanine and cytosine, which contain a single local CO stretching mode. To help interpret the nature of these vibrational modes, we performed density functional theory (DFT) calculations and found that multiple ring vibrations are coupled and delocalized over the purine and pyrimidine rings. Generally, there is close correspondence between the experimental and computational results, provided that the DFT calculations include explicit waters solvating hydrogen-bonding sites. These results provide direct experimental evidence of the delocalized nature of the nucleotide base vibrations via a nonperturbative fashion and will serve as building blocks for constructing a structure-based model of DNA and RNA vibrational spectroscopy.