Probing multiple core-hole interactions in the nitrogen K-edge of DNA base pairs by multidimensional attosecond X-ray spectroscopy. A simulation study

Two-dimensional X-ray correlation spectroscopy (2DXCS) signals of the isolated DNA bases and Watson-Crick base pairs which contain multiple absorbing nitrogen atoms are calculated. Core-hole excited states are calculated using density functional theory with the B3LYP functional and 6-311G** basis set. Sum over states calculations of the signals reveal changes in cross-peak intensities between hydrogen-bonded and stacked base pairs. Nucleobase analogues are proposed for investigating base-stacking and hydrogen-bonding interactions.     

A TDDFT study of the excited states of DNA bases and their assemblies

We present a detailed study of the optical absorption spectra of DNA bases and base pairs, carried out by means of time dependent density functional theory. The spectra for the isolated bases are compared to available theoretical and experimental data and used to assess the accuracy of the method and the quality of the exchange-correlation functional. Our approach turns out to be a reliable tool to describe the response of the nucleobases. Furthermore, we analyze in detail the impact of hydrogen bonding and pi-stacking in the calculated spectra for both Watson-Crick base pairs and Watson-Crick stacked assemblies. We show that the reduction of the UV absorption intensity (hypochromicity) for light polarized along the base-pair plane depends strongly on the type of interaction. For light polarized perpendicular to the basal plane, the hypochromicity effect is reduced, but another characteristic is found, namely a blue shift of the optical spectrum of the base-assembly compared to that of the isolated bases. The use of optical tools as fingerprints for the characterization of the structure (and type of interaction) is extensively discussed.     

Hydrogen bonding effect on the structure and vibrational spectra of complementary pairs of nucleic acid bases. I. Adenine-uracil

Vibrational spectra of the adenine—uracil complementary pair of nucleic acid bases (NAB) in the isolated state are calculated in the B3LYP/6-311++G(d,p) approximation and analyzed. The hydrogen bonding effect on frequency positions and intensities of normal vibrations of NAB pairs is shown in comparison with isolated uracil and adenine molecules.     

Electronic excitations and spectra in single-stranded DNA

Using density functional theory and molecular dynamics simulations, we show that delocalized states extending over three bases can be directly excited in single-stranded poly(A) DNA. The results are in semiquantitative agreement with recent experimental results for the delocalization length of these states in single- and double-stranded DNA. The structures used in these molecular dynamics calculations are validated by comparing calculated circular dichroic spectra for d(A)2 and d(A)4 with experiment. These spectra, which arise from highly stacked structures, are in good agreement with experiment, suggesting that the short delocalization in ssDNA arises in spite of strong stacking.     

Vibrational dynamics of DNA. II. Deuterium exchange effects and simulated IR absorption spectra

In Paper I, we studied vibrational properties of normal bases, base derivatives, Watson-Crick base pairs, and multiple layer base pair stacks in the frequency range of 1400-1800 cm(-1). However, typical IR absorption spectra of single- and double-stranded DNA have been measured in D(2)O solution. Consequently, the more relevant bases and base pairs are those with deuterium atoms in replacement with labile amino hydrogen atoms. Thus, we have carried out density functional theory vibrational analyses of properly deuterated bases, base pairs, and stacked base pair systems. In the frequency range of interest, both aromatic ring deformation modes and carbonyl stretching modes appear to be strongly IR active. Basis mode frequencies and vibrational coupling constants are newly determined and used to numerically simulate IR absorption spectra. It turns out that the hydration effects on vibrational spectra are important. The numerically simulated vibrational spectra are directly compared with experiments. Also, the (18)O-isotope exchange effect on the poly(dG):poly(dC) spectrum is quantitatively described. The present calculation results will be used to further simulate two-dimensional IR photon echo spectra of DNA oligomers in the companion Paper III.     

Vibronic Structures in Absorption and Fluorescence Spectra of Firefly Oxyluciferin in Aqueous Solutions

To elucidate the factors determining the spectral shapes and widths of the absorption and fluorescence spectra for keto and enol oxyluciferin and their conjugate bases in aqueous solutions, the intensities of vibronic transitions between their ground and first electronic excited states were calculated for the first time via estimation of the vibrational Franck–Condon factors. The major normal modes, overtones and combination tones in absorption and fluorescence spectra are similar for all species. The theoretical full widths at half maximum of absorption spectra are 0.4–0.7 eV and those for the fluorescence spectra are 0.4–0.5 eV, except for phenolate‐keto that exhibits exceptionally sharp peak widths due to the dominance of the 0–0′ or 0′–0 band. These spectral shapes and widths explain many relevant features of the experimentally observed spectra.     

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     

Comparison of the secondary-ion and electron-impact mass spectra of alkaloids

An advantage of using secondary-ion spectra for determining the molecular masses of quaternary bases and bases the molecular ions of which are unstable to electron impact is shown. The LSIMS spectra of four groups of diterpene alkaloids include the 100% peaks of the (M + H)⁺ ions and the peaks of the main fragments present in the EI spectra. The B/E = const. spectra of the key ions obtained by the different methods are compared. The closeness of the values of A of analogous transitions calculated from the B/E and MD spectra is confirmed.     

Theoretical calculation of the photodetachment spectra of XAuY- (X, Y = Cl, Br, and I)

The photodetachment spectra of the title molecules have been calculated, taking electron correlation and spin-orbit coupling into account and employing improved relativistic effective core potentials for gold and the halogen atoms. The calculated spectra have been compared with existing experimental spectra. The spin-orbit splitting of several degenerate electronic states has been calculated. The composition of the spin-orbit eigenstates are analyzed in terms of scalar relativistic electronic states. A comparison of the relative position of peaks in the calculated photodetachment spectra of the title molecules has been made.     

Calculation of the stabilization energies of oxidatively damaged guanine base pairs with guanine

DNA is constantly exposed to endogenous and exogenous oxidative stresses. Damaged DNA can cause mutations, which may increase the risk of developing cancer and other diseases. G:C-C:G transversions are caused by various oxidative stresses. 2,2,4-Triamino-5(2H)-oxazolone (Oz), guanidinohydantoin (Gh)/iminoallantoin (Ia) and spiro-imino-dihydantoin (Sp) are known products of oxidative guanine damage. These damaged bases can base pair with guanine and cause G:C-C:G transversions. In this study, the stabilization energies of these bases paired with guanine were calculated in vacuo and in water. The calculated stabilization energies of the Ia:G base pairs were similar to that of the native C:G base pair, and both bases pairs have three hydrogen bonds. By contrast, the calculated stabilization energies of Gh:G, which form two hydrogen bonds, were lower than the Ia:G base pairs, suggesting that the stabilization energy depends on the number of hydrogen bonds. In addition, the Sp:G base pairs were less stable than the Ia:G base pairs. Furthermore, calculations showed that the Oz:G base pairs were less stable than the Ia:G, Gh:G and Sp:G base pairs, even though experimental results showed that incorporation of guanine opposite Oz is more efficient than that opposite Gh/Ia and Sp.     

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.     

Simulation of X-ray absorption near edge spectra of organometallic compounds in the ground and optically excited states

The Mg K-edge and Zn K- and L3-edge X-ray absorption near edge spectra of Mg and Zn porphyrins in the ground state and low-lying optically excited states are calculated. Also computed are X-ray absorption near edge spectra of Fe(II) spin crossover compound in its ground and low-lying optically excited states, motivated by a recent experiment (J. Phys. Chem. A 2006, 110, 38). The calculated absorption spectra of optically excited states can be used to simulate ultrafast optical pump/X-ray probe experiments.     

Absorption and fluorescence emission spectroscopic characters of naphtho‐homologated yy‐DNA bases and effect of methanol solution and base pairing

A comprehensive theoretical study of electronic transitions of naphtho‐homologated base analogs, namely, yy‐T , yy‐C , yy‐A , and yy‐G , was performed. The nature of the low‐lying excited states is discussed, and the results are compared with those from experiment and also with those of y‐bases. Geometrical characteristics of the lowest excited singlet ππ* and nπ* states were explored using the CIS method, and the effects of methanol solution and paring with their complementary natural bases on the relevant absorption and emission spectra of these modified bases were examined. The calculated excitation and emission energies agree well with the measured data, where experimental results are available. In methanol solution, the fluorescence from yy‐A and yy‐G would be expected to occur around 539 and 562 nm, respectively, suggesting that yy‐A is a green‐colored fluorophore, whereas yy‐G is a yellow‐colored fluorophore. The methanol solution was found to red‐shift both the absorption and emission maxima of yy‐A , yy‐T , and yy‐C , but blue‐shift those for yy‐G . Generally, though base pairing has no significant effects on the absorption and fluorescence maxima of yy‐A , yy‐C , and yy‐T , it blue‐shifts those for yy‐G . © 2009 Wiley Periodicals, Inc. J Comput Chem 2010     

Synthesis and characterization of 3-amino-5-methylisoxazole Schiff bases and their complexes with copper(II), nickel(II) and cobalt(II)

Schiff bases obtained by condensing 3-amino-5-methylisoxazole with salicylaldehyde, 2,3-dihydroxybenzaldehyde, 2,4-dihydroxy-benzaldehyde, 2,5-dihydroxybenzaldehyde or o-hydroxynaphthaldehyde were obtained and characterized by C, H, N analysis, mass, NMR and IR spectra. Copper, nickel and cobalt complexes of the Schiff bases were prepared and characterized using elemental analysis, conductivity measurements, magnetic moments, IR, UV-VIS and ESR spectra, X-ray diffraction, TGA, DTA and DSC thermal analysis. All the complexes are non-electrolytes. ESR spectra show isotropic as well as axial symmetry for the copper complexes. Thermal studies support the formulation of these complexes and showed that they decompose in two or three steps depending on the metal used. Activation energy E _a and enthalpies ΔH associated with the decomposition process were calculated and correlated with the complexed metal used.     

The dynamics of intermittent strand separation in double-stranded DNA

The transient rupture and reformation of hydrogen bonds between base pairs on distinct chains of double-stranded DNA ("bubble" dynamics) is modeled in terms of the fluctuating distance between the bases. The fluctuations in the distance are assumed to be governed by a simple Langevin equation with a quadratic potential under conditions of high friction. A critical distance of separation L must be achieved before a bubble defect is considered to have been formed. The decay of the dynamic correlations between states of the DNA that have such defects and those that do not has been calculated from the above model and has been found to reproduce the trends in experimental measurements of the same quantity.     

UV spectra and excitation delocalization in DNA: influence of the spectral width.

The singlet excited states of the model DNA duplex (dA)10.(dT)10 are studied. Calculations are performed in the exciton theory framework. Molecular dynamics calculations provide the duplex geometry. The dipolar coupling is determined using atomic transition charges. The monomer transition energies are simulated by Gaussian functions resembling the absorption bands of nucleosides in aqueous solutions. Most of the excited states are found to be delocalized over at least two bases and result from the mixing of different monomer states. Their properties are only weakly affected by conformational changes of the double helix. On average, the highest oscillator strength is carried by the upper eigenstates. The duplex absorption spectra are shifted a few nanometers to higher energies with respect to the spectra of noninteracting monomers. The states with larger spatial extent are located close to the maximum of the absorption spectrum.     

Estimates of electronic coupling for excess electron transfer in DNA

Electronic coupling V(da) is one of the key parameters that determine the rate of charge transfer through DNA. While there have been several computational studies of V(da) for hole transfer, estimates of electronic couplings for excess electron transfer (ET) in DNA remain unavailable. In the paper, an efficient strategy is established for calculating the ET matrix elements between base pairs in a pi stack. Two approaches are considered. First, we employ the diabatic-state (DS) method in which donor and acceptor are represented with radical anions of the canonical base pairs adenine-thymine (AT) and guanine-cytosine (GC). In this approach, similar values of V(da) are obtained with the standard 6-31G(*) and extended 6-31+ +G(**) basis sets. Second, the electronic couplings are derived from lowest unoccupied molecular orbitals (LUMOs) of neutral systems by using the generalized Mulliken-Hush or fragment charge methods. Because the radical-anion states of AT and GC are well reproduced by LUMOs of the neutral base pairs calculated without diffuse functions, the estimated values of V(da) are in good agreement with the couplings obtained for radical-anion states using the DS method. However, when the calculation of a neutral stack is carried out with diffuse functions, LUMOs of the system exhibit the dipole-bound character and cannot be used for estimating electronic couplings. Our calculations suggest that the ET matrix elements V(da) for models containing intrastrand thymine and cytosine bases are essentially larger than the couplings in complexes with interstrand pyrimidine bases. The matrix elements for excess electron transfer are found to be considerably smaller than the corresponding values for hole transfer and to be very responsive to structural changes in a DNA stack.     

The Generalized Sturmian Method for Calculating Spectra of Atoms and Ions

The properties of generalized Sturmian basis sets are reviewed, and functions of this type are used to perform direct configuration interaction calculations on the spectra of atoms and ions. Singlet excited states calculated in this way show good agreement with experimentally measured spectra. When the generalized Sturmian method is applied to atoms, the configurations are constructed from hydrogenlike atomic orbitals with an effective charge which is characteristic of the configuration. Thus, orthonormality between the orbitals of different configurations cannot be assumed, and the generalized Slater–Condon rules must be used. This aspect of the problem is discussed in detail. Finally spectra are calculated in the presence of a strong external electric field. In addition to the expected Stark effect, the calculated spectra exhibit anomalous states. These are shown to be states where one of the electrons is primarily outside the atom or ion, with only a small amplitude inside.     

Structure, stability, and dynamics of canonical and noncanonical base pairs: quantum chemical studies

The importance of non-Watson-Crick base pairs in the three-dimensional structure of RNA is now well established. The structure and stability of these noncanonical base pairs are, however, poorly understood. We have attempted to understand structural features of 33 frequently occurring base pairs using density functional theory. These are of three types, namely (i) those stabilized by two or more polar hydrogen bonds between the bases, (ii) those having one polar and another C-H...O/N type interactions, and (iii) those having one H-bond between the bases and another involving one of the sugars linked to the bases. We found that the base pairs having two polar H-bonds are very stable as compared to those having one C-H...O/N interaction. Our quantitatively analysis of structures of these optimized base pairs indicates that they possess a different amount of nonplanarity with large propeller or buckle values as also observed in the crystal structures. We further found that geometry optimization does not modify the hydrogen-bonding pattern, as values of shear and open angle of the base pairs remain conserved. The structures of initial crystal geometry and final optimized geometry of some base pairs having only one polar H-bond and a C-H...O/N interaction, however, are significantly different, indicating the weak nature of the nonpolar interaction. The base pair flexibility, as measured from normal-mode analysis, in terms of the intrinsic standard deviations of the base pair structural parameters are in conformity with those calculated from RNA crystal structures. We also noticed that deformation of a base pair along the stretch direction is impossible for all of the base pairs, and movements of the base pairs along shear and open are also quite restricted. The base pair opening mode through alteration of propeller or buckle is considerably less restricted for most of the base pairs.     

Photophysical characters of rationally designed hetero-ring-expanded guanine analogues and effect of cytosine pairing

We present the results of the CIS and TDB3LYP calculations of the optical absorption and emission spectra of some newly designed guanine (G) analogues and their Watson-Crick base pairs. Compared with natural G, the onset absorption peaks of these newly designed analogues are red-shifted, while the fluorescence peaks are blue-shifted. In general, the first excited singlet states (pipi*) of these analogues are nonplanar for all bases considered here. But, the Stokes shifts for the designed G-analogues are much smaller than that of natural G, suggesting that they have stronger molecular rigidity and higher fluorescence quantum yields than those of natural G. The first excited states of these Watson-Crick base pairs essentially originate from those of their isolated purine moieties, as demonstrated from the S1 geometries of their Watson-Crick base pairs. For G and its analogues, A1 and A2 (they are ring-expanded with one-bond intercalation at the C5 site), the pairing with cytosine reduces the oscillator strengths of both the first absorption peak (by 27%-60%) and the fluorescent emission (by 19%-23%), while for the analogues A3, A4, and xG in which G is ring-expanded with a two-bond intercalation at the C5 site, the pairing, in contrast, increases the oscillator strengths of both the first absorption peak (by 11%-15%) and the fluorescent emission (by 3%-20%). These observations indicate that the pairing with cytosine can quench the fluorescence for G, A1, and A2 but enhance the fluorescence quantum yields for A3, A4, and xG. The significant shifts induced by ring-expansion in the ring-expanded G with a two-bond intercalation at the C5 site reveal a possibility for their fluorescent detections.