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

Beyond the "spine of hydration": Chiral SFG spectroscopy detects DNA first hydration shell and base pair structures

Experimental methods capable of selectively probing water at the DNA minor groove, major groove, and phosphate backbone are crucial for understanding how hydration influences DNA structure and function. Chiral-selective sum frequency generation spectroscopy (chiral SFG) is unique among vibrational spectroscopies because it can selectively probe water molecules that form chiral hydration structures around biomolecules. However, interpreting chiral SFG spectra is challenging since both water and the biomolecule can produce chiral SFG signals. Here, we combine experiment and computation to establish a theoretical framework for the rigorous interpretation of chiral SFG spectra of DNA. We demonstrate that chiral SFG detects the N-H stretch of DNA base pairs and the O-H stretch of water, exclusively probing water molecules in the DNA first hydration shell. Our analysis reveals that DNA transfers chirality to water molecules only within the first hydration shell, so they can be probed by chiral SFG spectroscopy. Beyond the first hydration shell, the electric field-induced water structure is symmetric and, therefore, precludes chiral SFG response. Furthermore, we find that chiral SFG can differentiate chiral subpopulations of first hydration shell water molecules at the minor groove, major groove, and phosphate backbone. Our findings challenge the scientific perspective dominant for more than 40 years that the minor groove "spine of hydration" is the only chiral water structure surrounding the DNA double helix. By identifying the molecular origins of the DNA chiral SFG spectrum, we lay a robust experimental and theoretical foundation for applying chiral SFG to explore the chemical and biological physics of DNA hydration.

Modeling the Infrared Spectroscopy of Oligonucleotides with 13C Isotope Labels

The carbonyl stretching modes have been widely used in linear and two-dimensional infrared (IR) spectroscopy to probe the conformation, interaction, and biological functions of nucleic acids. However, due to their universal appearance in nucleobases, the IR absorption bands of nucleic acids are often highly congested in the 1600-1800 cm^-1 region. Following the fruitful applications in proteins, ¹³C isotope labels have been introduced to the IR measurements of oligonucleotides to reveal their site-specific structural fluctuations and hydrogen bonding conditions. In this work, we combine recently developed frequency and coupling maps to develop a theoretical strategy that models the IR spectra of oligonucleotides with ¹³C labels directly from molecular dynamics simulations. We apply the theoretical method to nucleoside 5'-monophosphates and DNA double helices and demonstrate how elements of the vibrational Hamiltonian determine the spectral features and their changes upon isotope labeling. Using the double helices as examples, we show that the calculated IR spectra are in good agreement with experiments and the ¹³C isotope labeling technique can potentially be applied to characterize the stacking configurations and secondary structures of nucleic acids.

Sequence-Dependent Melting and Refolding Dynamics of RNA UNCG Tetraloops Using Temperature-Jump/Drop Infrared Spectroscopy

Time-resolved temperature-jump/drop infrared (IR) spectroscopy has been used to measure the impact of stem base sequence on the melting and refolding dynamics of ribonucleic acid (RNA) tetraloops. A series of three 12-nucleotide RNA hairpin sequences were studied, each featuring a UACG tetraloop motif and a double-stranded stem containing four base pairs. In each case, the stem comprised three GC pairs plus a single AU base pair inserted at the closing point of the loop (RNA_loop), in the middle of the stem (RNA_mid), or at the stem terminus (RNA_end). Results from analogous DNA tetraloop (TACG) sequences were also obtained. Inclusion of AU or AT base pairs in the stem leads to faster melting of the stem-loop structure compared to a stem sequence featuring four GC base pairs while refolding times were found to be slower, consistent with a general reduction in stem-loop stability caused by the AU/AT pair. Independent measurement of the dynamic timescales for melting and refolding of ring vibrational modes of guanine (G_R) and adenine (A_R) provided position-specific insight into hairpin dynamics. The G_R-derived data showed that DNA sequences melted more quickly (0.5 ± 0.1 to 0.7 ± 0.1 μs at 70 °C) than analogous RNA sequences (4.3 ± 0.4 to 4.4 ± 0.3 μs at 70 °C). Position-sensitive data from the A_R modes suggests that DNA hairpins begin melting from the terminal end of the stem toward the loop while RNA sequences begin melting from the loop. Refolding timescales for both RNA and DNA hairpins were found to be similar (250 ± 50 μs at 70 °C) except for RNA_end and DNA_loop which refolded much more slowly (746 ± 36 and 430 ± 31 μs, respectively), showing that the refolding pathway is significantly impaired by the placement of AU/AT pairs at different points in the stem. We conclude that conformational changes of analogous pairs of RNA and DNA tetraloops proceed by different mechanisms.

Measuring RNA UNCG Tetraloop Refolding Dynamics Using Temperature-Jump/Drop Infrared Spectroscopy

Determining the structural dynamics of RNA and DNA is essential to understanding their cellular function, but direct measurement of strand association or folding remains experimentally challenging. Here we illustrate a temperature-jump/drop method able to reveal refolding dynamics. Time-resolved temperature-jump/drop infrared spectroscopy is used to measure the melting and refolding dynamics of a 12-nucleotide RNA sequence comprising a UACG tetraloop and a four-base-pair double-stranded GC stem, comparing them to an equivalent DNA (TACG) sequence. Stem-loop melting occurred an order of magnitude more slowly in RNA than DNA (6.0 ± 0.1 μs versus 0.8 ± 0.1 μs at 70 °C). In contrast, the refolding dynamics of both sequences occurred on similar time scales (200 μs). While the melting and refolding dynamics of RNA and DNA hairpins both followed Arrhenius temperature dependences, refolding was characterized by an apparent negative activation energy, consistent with a mechanism involving multiple misfolded intermediates prior to zipping of the stem base pairs.

Diffusion mechanisms of DNA in agarose gels: NMR studies and Monte Carlo simulations

We report on the diffusion mechanism of short, single-stranded DNA molecules with up to 100 nucleobases in agarose gels with concentrations of up to 2.0% with the aim to characterize the DNA–agarose interaction. The diffusion coefficients were measured directly, i.e., without any model assumptions, by pulsed field gradient nuclear magnetic resonance (PFG-NMR). We find that the diffusion coefficient decreases, as expected, with an increase in both DNA strand length and gel concentration. In addition, we performed Monte Carlo simulations of particle diffusion in a model network of polymer chains, considering our experimental conditions. Together, the Monte Carlo simulations and the PFG-NMR results show that the decrease in diffusion coefficients in the presence of the agarose gel is due to a temporary adhesion of the DNA molecules to the surface of gel fibers. The average adhesion time to a given gel fiber increases with the length of the DNA strands but is independent of the number of gel fibers. The corresponding magnitude of the binding enthalpies of DNA strands to gel fibers indicates that a mixture of van der Waals interactions and hydrogen bonding contributes to the decreased diffusion of DNA in agarose gels.

The polarization dependence of 2D IR cross-peaks distinguishes parallel-stranded and antiparallel-stranded DNA G-quadruplexes

Guanine-rich nucleic acid sequences have a tendency to form four-stranded non-canonical motifs known as G-quadruplexes. These motifs may adopt a wide range of structures characterized by size, strand orientation, guanine base conformation, and fold topology. Using three K⁺-bound model systems, we show that vibrational coupling between guanine C6 = O and ring modes varies between parallel-stranded and antiparallel-stranded G-quadruplexes, and that such structures can be distinguished by comparison of the polarization dependences of cross-peaks in their two-dimensional infrared (2D IR) spectra. Combined with previously defined vibrational frequency trends, this analysis reveals key features of a 30-nucleotide unimolecular variant of the Bcl-2 proximal promoter that are consistent with its reported structure. This study shows that 2D IR spectroscopy is a convenient method for analyzing G-quadruplex structures that can be applied to complex sequences where traditional high-resolution methods are limited by solubility and disorder.

Probing the Hydrogen-Bonding Environment of Individual Bases in DNA Duplexes with Isotope-Edited Infrared Spectroscopy

Measuring the strength of the hydrogen bonds between DNA base pairs is of vital importance for understanding how our genetic code is physically accessed and recognized in cells, particularly during replication and transcription. Therefore, it is important to develop probes for these key hydrogen bonds (H-bonds) that dictate events critical to cellular function, such as the localized melting of DNA. The vibrations of carbonyl bonds are well-known probes of their H-bonding environment, and their signals can be observed with infrared (IR) spectroscopy. Yet, pinpointing a single bond of interest in the complex IR spectrum of DNA is challenging due to the large number of carbonyl signals that overlap with each other. Here, we develop a method using isotope editing and infrared (IR) spectroscopy to isolate IR signals from the thymine (T) C2=O carbonyl. We use solvatochromatic studies to show that the TC2=O signal’s position in the IR spectrum is sensitive to the H-bonding capacity of the solvent. Our results indicate that C2=O of a single T base within DNA duplexes experiences weak H-bonding interactions. This finding is consistent with the existence of a third, noncanonical CH···O H-bond between adenine and thymine in both Watson–Crick and Hoogsteen base pairs in DNA.

Hydrogen Bonding in Natural and Unnatural Base Pairs-A Local Vibrational Mode Study

In this work hydrogen bonding in a diverse set of 36 unnatural and the three natural Watson Crick base pairs adenine (A)-thymine (T), adenine (A)-uracil (U) and guanine (G)-cytosine (C) was assessed utilizing local vibrational force constants derived from the local mode analysis, originally introduced by Konkoli and Cremer as a unique bond strength measure based on vibrational spectroscopy. The local mode analysis was complemented by the topological analysis of the electronic density and the natural bond orbital analysis. The most interesting findings of our study are that (i) hydrogen bonding in Watson Crick base pairs is not exceptionally strong and (ii) the N-H⋯N is the most favorable hydrogen bond in both unnatural and natural base pairs while O-H⋯N/O bonds are the less favorable in unnatural base pairs and not found at all in natural base pairs. In addition, the important role of non-classical C-H⋯N/O bonds for the stabilization of base pairs was revealed, especially the role of C-H⋯O bonds in Watson Crick base pairs. Hydrogen bonding in Watson Crick base pairs modeled in the DNA via a QM/MM approach showed that the DNA environment increases the strength of the central N-H⋯N bond and the C-H⋯O bonds, and at the same time decreases the strength of the N-H⋯O bond. However, the general trends observed in the gas phase calculations remain unchanged. The new methodology presented and tested in this work provides the bioengineering community with an efficient design tool to assess and predict the type and strength of hydrogen bonding in artificial base pairs.

Tracing the Photoaddition of Pharmaceutical Psoralens to DNA

The psoralens 8-methoxypsoralen (8-MOP), 4,5',8-trimethylpsoralen (TMP) and 5-methoxypsoralen (5-MOP) find clinical application in PUVA (psoralen + UVA) therapy. PUVA treats skin diseases like psoriasis and atopic eczema. Psoralens target the DNA of cells. Upon photo-excitation psoralens bind to the DNA base thymine. This photo-binding was studied using steady-state UV/Vis and IR spectroscopy as well as nanosecond transient UV/Vis absorption. The experiments show that the photo-addition of 8-MOP and TMP involve the psoralen triplet state and a biradical intermediate. 5-MOP forms a structurally different photo-product. Its formation could not be traced by the present spectroscopic technique.

Vibrational Spectroscopic Map, Vibrational Spectroscopy, and Intermolecular Interaction

Vibrational spectroscopy is an essential tool in chemical analyses, biological assays, and studies of functional materials. Over the past decade, various coherent nonlinear vibrational spectroscopic techniques have been developed and enabled researchers to study time-correlations of the fluctuating frequencies that are directly related to solute-solvent dynamics, dynamical changes in molecular conformations and local electrostatic environments, chemical and biochemical reactions, protein structural dynamics and functions, characteristic processes of functional materials, and so on. In order to gain incisive and quantitative information on the local electrostatic environment, molecular conformation, protein structure and interprotein contacts, ligand binding kinetics, and electric and optical properties of functional materials, a variety of vibrational probes have been developed and site-specifically incorporated into molecular, biological, and material systems for time-resolved vibrational spectroscopic investigation. However, still, an all-encompassing theory that describes the vibrational solvatochromism, electrochromism, and dynamic fluctuation of vibrational frequencies has not been completely established mainly due to the intrinsic complexity of intermolecular interactions in condensed phases. In particular, the amount of data obtained from the linear and nonlinear vibrational spectroscopic experiments has been rapidly increasing, but the lack of a quantitative method to interpret these measurements has been one major obstacle in broadening the applications of these methods. Among various theoretical models, one of the most successful approaches is a semiempirical model generally referred to as the vibrational spectroscopic map that is based on a rigorous theory of intermolecular interactions. Recently, genetic algorithm, neural network, and machine learning approaches have been applied to the development of vibrational solvatochromism theory. In this review, we provide comprehensive descriptions of the theoretical foundation and various examples showing its extraordinary successes in the interpretations of experimental observations. In addition, a brief introduction to a newly created repository Web site (http://frequencymap.org) for vibrational spectroscopic maps is presented. We anticipate that a combination of the vibrational frequency map approach and state-of-the-art multidimensional vibrational spectroscopy will be one of the most fruitful ways to study the structure and dynamics of chemical, biological, and functional molecular systems in the future.

Modeling the vibrational couplings of nucleobases

Vibrational spectroscopy, in particular infrared spectroscopy, has been widely used to probe the three-dimensional structures and conformational dynamics of nucleic acids. As commonly used chromophores, the C=O and C=C stretch modes in the nucleobases exhibit distinct spectral features for different base pairing and stacking configurations. To elucidate the origin of their structural sensitivity, in this work, we develop transition charge coupling (TCC) models that allow one to efficiently calculate the interactions or couplings between the C=O and C=C chromophores based on the geometric arrangements of the nucleobases. To evaluate their performances, we apply the TCC models to DNA and RNA oligonucleotides with a variety of secondary and tertiary structures and demonstrate that the predicted couplings are in quantitative agreement with the reference values. We further elucidate how the interactions between the paired and stacked bases give rise to characteristic IR absorption peaks and show that the TCC models provide more reliable predictions of the coupling constants as compared to the transition dipole coupling scheme. The TCC models, together with our recently developed through-bond coupling constants and vibrational frequency maps, provide an effective theoretical strategy to model the vibrational Hamiltonian, and hence the vibrational spectra of nucleic acids in the base carbonyl stretch region directly from atomistic molecular simulations.

Electron Transfer in DNA at Electrified Interfaces

The ability of the DNA double helix to transport electrons underlies many life-centered biological processes and bio-electronic applications. However, there is little consensus on how efficiently the base pair π-stacks of DNA mediate electron transport. This minireview scrutinizes the current state-of-the-art knowledge on electron transfer (ET) properties of DNA and its long-range ability to transfer (mediate) electrical signals at electrified interfaces, without being oxidized or reduced. Complex changes an electric field induces in the DNA structure and its electronic properties govern the efficiency of DNA-mediated ET at electrodes and allow addressing the existing phenomenological riddles, while recently discovered rectifying properties of DNA contribute both to our understanding of DNA's ET in living systems and to advances in molecular bioelectronics.

Monitoring Base-Specific Dynamics during Melting of DNA-Ligand Complexes Using Temperature-Jump Time-Resolved Infrared Spectroscopy

Ultrafast time-resolved infrared spectroscopy employing nanosecond temperature-jump initiation has been used to study the melting of double-stranded (ds)DNA oligomers in the presence and absence of minor groove-binding ligand Hoechst 33258. Ligand binding to ds(5'-GCAAATTTCC-3'), which binds Hoechst 33258 in the central A-tract region with nanomolar affinity, causes a dramatic increase in the timescales for strand melting from 30 to ∼250 μs. Ligand binding also suppresses premelting disruption of the dsDNA structure, which takes place on 100 ns timescales and includes end-fraying. In contrast, ligand binding to the ds(5'-GCATATATCC-3') sequence, which exhibits an order of magnitude lower affinity for Hoechst 33258 than the A-tract motif, leads to an increase by only a factor of 5 in melting timescales and reduced suppression of premelting sequence perturbation and end-fraying. These results demonstrate a dynamic impact of the minor groove ligand on the dsDNA structure that correlates with binding strength and thermodynamic stabilization of the duplex. Moreover, the ability of the ligand to influence base pairs distant from the binding site has potential implications for allosteric communication mechanisms in dsDNA.

Development of Vibrational Frequency Maps for Nucleobases

Vibrational spectroscopy provides a powerful tool to probe the structure and dynamics of nucleic acids because specific normal modes, particularly the base carbonyl stretch modes, are highly sensitive to the hydrogen bonding patterns and stacking configurations in these biomolecules. In this work, we develop vibrational frequency maps for the C═O and C═C stretches in nucleobases that allow the calculations of their site frequencies directly from molecular dynamics simulations. We assess the frequency maps by applying them to nucleobase derivatives in aqueous solutions and nucleosides in organic solvents and demonstrate that the predicted infrared spectra are in good agreement with experimental measurements. The frequency maps can be readily used to model the linear and nonlinear vibrational spectroscopy of nucleic acids and elucidate the molecular origin of the experimentally observed spectral features.

Spectroelectrochemical Characterization of 1,2-Dipalmitoyl- sn-glycero-3-cytidine Diphosphate Nucleolipid Monolayer Supported on Gold (111) Electrode

The effect of the electrode potential on the orientation and conformation of the 1,2-dipalmitoyl- sn-glycero-3-cytidine monolayer deposited on a gold (111) electrode surface was described. The potential of zero free charge ( E_pzc) for the monolayer-covered electrode was determined to be -0.2 V vs SCE. The differential capacitance and charge density data indicated that the monolayer is stable at the electrode surface when ( E - E_pzc) > 0.0 V. At negative rational potentials, a progressive detachment (electrodewetting) of the monolayer occurs. The monolayer is fully detached from the electrode surface at ( E - E_pzc) < -0.6 V. The conformation and orientation of the acyl chains and the orientation of the cytosine moiety were determined with the help of photon polarization modulation infrared reflection absorption spectroscopy (PM-IRRAS). The IR measurements demonstrate that the acyl chains are predominantly in the gel phase in the adsorbed state and tilted at an angle of ∼30° with respect to the electrode surface normal. The tilt angle of the acyl chains increases when the film is detached from the gold surface, indicating that the monolayer becomes more disordered. At ( E - E_pzc) > 0.0 V, the plane of the cytosine moiety assumes a small angle of ∼20° with respect to the surface. At negative potentials, the tilt angle of the cytosine fragment increases and rotates. With the help of DFT calculations, these changes were explained by the repulsion of the positive pole of the cytosine permanent dipole moment by the positively charged gold surface and its attraction to the metal surface at negative electrode potentials. This work provides unique information for the future development of sensors based on the molecular recognition of nucleoside targets.

Submolecular Structure and Orientation of Oligonucleotide Duplexes Tethered to Gold Electrodes Probed by Infrared Reflection Absorption Spectroscopy: Effect of the Electrode Potentials

Unique electronic and ligand recognition properties of the DNA double helix provide basis for DNA applications in biomolecular electronic and biosensor devices. However, the relation between the structure of DNA at electrified interfaces and its electronic properties is still not well understood. Here, potential-driven changes in the submolecular structure of DNA double helices composed of either adenine-thymine (dAdT)₂₅ or cytosine-guanine (dGdC)₂₀ base pairs tethered to the gold electrodes are for the first time analyzed by in situ polarization modulation infrared reflection absorption spectroscopy (PM IRRAS) performed under the electrochemical control. It is shown that the conformation of the DNA duplexes tethered to gold electrodes via the C₆ alkanethiol linker strongly depends on the nucleic acid sequence composition. The tilt of purine and pyrimidine rings of the complementary base pairs (dAdT and dGdC) depends on the potential applied to the electrode. By contrast, neither the conformation nor orientation of the ionic in character phosphate-sugar backbone is affected by the electrode potentials. At potentials more positive than the potential of zero charge (pzc), a gradual tilting of the double helix is observed. In this tilted orientation, the planes of the complementary purine and pyrimidine rings lie ideally parallel to each other. These potentials do not affect the integral stability of the DNA double helix at the charged interface. At potentials more negative than the pzc, DNA helices adopt a vertical to the gold surface orientation. Tilt of the purine and pyrimidine rings depends on the composition of the double helix. In monolayers composed of (dAdT)₂₅ molecules the rings of the complementary base pairs lie parallel to each other. By contrast, the tilt of purine and pyrimidine rings in (dGdC)₂₀ helices depends on the potential applied to the electrode. Such potential-induced mobility of the complementary base pairs can destabilize the helix structure at a submolecular level. These pioneer results on the potential-driven changes in the submolecular structure of double stranded DNA adsorbed on conductive supports contribute to further understanding of the potential-driven sequence-specific electronic properties of surface-tethered oligonucleotides.

2D-IR Spectroscopy Shows that Optimized DNA Minor Groove Binding of Hoechst33258 Follows an Induced Fit Model

The induced fit binding model describes a conformational change occurring when a small molecule binds to its biomacromolecular target. The result is enhanced noncovalent interactions between the ligand and biomolecule. Induced fit is well-established for small molecule-protein interactions, but its relevance to small molecule-DNA binding is less clear. We investigate the molecular determinants of Hoechst33258 binding to its preferred A-tract sequence relative to a suboptimal alternating A-T sequence. Results from two-dimensional infrared spectroscopy, which is sensitive to H-bonding and molecular structure changes, show that Hoechst33258 binding results in loss of the minor groove spine of hydration in both sequences, but an additional perturbation of the base propeller twists occurs in the A-tract binding region. This induced fit maximizes favorable ligand-DNA enthalpic contributions in the optimal binding case and demonstrates that controlling the molecular details that induce subtle changes in DNA structure may hold the key to designing next-generation DNA-binding molecules.

Ultrafast 2D-IR and optical Kerr effect spectroscopy reveal the impact of duplex melting on the structural dynamics of DNA

Changes in the structural and solvation dynamics of DNA upon duplex melting are observed by 2D-IR and optical Kerr-effect spectroscopies.

Long-Range Vibrational Dynamics Are Directed by Watson-Crick Base Pairing in Duplex DNA

Ultrafast two-dimensional infrared (2D-IR) spectroscopy of a 15-mer A-T DNA duplex in solution has revealed structure-dependent vibrational coupling and energy transfer processes linking bases with the sugar-phosphate backbone. Duplex melting induces significant changes in the positions of off-diagonal peaks linking carbonyl and ring-stretching vibrational modes of the adenine and thymine bases with vibrations of the phosphate group and phosphodiester linkage. These indicate that Watson-Crick hydrogen bonding and helix formation lead to a unique vibrational coupling arrangement of base vibrational modes with those of the phosphate unit. On the basis of observations from time-resolved 2D-IR data, we conclude that rapid energy transfer processes occur between base and backbone, mediated by additional modes located on the deoxyribose moiety within the same nucleotide. These relaxation dynamics are insensitive to duplex melting, showing that efficient intramolecular energy relaxation to the solvent via the phosphate groups is the key to excess energy dissipation in both single- and double-stranded DNA.

Structural Characterization of Single-Stranded DNA Monolayers Using Two-Dimensional Sum Frequency Generation Spectroscopy

DNA-covered materials are important in technological applications such as biosensors and microarrays, but obtaining structural information on surface-bound biomolecules is experimentally challenging. In this paper, we structurally characterize single-stranded DNA monolayers of poly(thymine) from 10 to 25 bases in length with an emerging surface technique called two-dimensional sum frequency generation (2D SFG) spectroscopy. These experiments are carried out by adding a mid-IR pulse shaper to a femtosecond broad-band SFG spectrometer. Cross peaks and 2D line shapes in the 2D SFG spectra provide information about structure and dynamics. Because the 2D SFG spectra are heterodyne detected, the monolayer spectra can be directly compared to 2D infrared (2D IR) spectra of poly(thymine) in solution, which aids interpretation. We simulate the 2D SFG spectra using DFT calculations and an excitonic Hamiltonian that relates the molecular geometry to the vibrational coupling. Intrabase cross peaks help define the orientation of the bases and interbase cross peaks, created by coupling between bases, and resolves features not observed in 1D SFG spectra that constrain the relative geometries of stacked bases. We present a structure for the poly(T) oligomer that is consistent with the 2D SFG data. These experiments provide insight into the DNA monolayer structure and set precedent for studying complex biomolecules on surfaces with 2D SFG spectroscopy.

Characterizing Anharmonic Vibrational Modes of Quinones with Two-Dimensional Infrared Spectroscopy

Two-dimensional infrared (2D IR) spectroscopy was used to study the vibrational modes of three quinones--benzoquinone, naphthoquinone, and anthraquinone. The vibrations of interest were in the spectral range of 1560-1710 cm(-1), corresponding to the in-plane carbonyl and ring stretching vibrations. Coupling between the vibrational modes is indicated by the cross peaks in the 2D IR spectra. The diagonal and off-diagonal anharmonicities range from 4.6 to 17.4 cm(-1) for the quinone series. In addition, there is significant vibrational coupling between the in-plane carbonyl and ring stretching vibrations. The diagonal anharmonicity, off-diagonal anharmonicity, and vibrational coupling constants are reported for benzoquinone, naphthoquinone, and anthraquinone.

Applications of two-dimensional infrared spectroscopy

Two-dimensional infrared (2D IR) spectroscopy has recently emerged as a powerful tool with applications in many areas of scientific research. The inherent high time resolution coupled with bond-specific spatial resolution of IR spectroscopy enable direct characterization of rapidly interconverting species and fast processes, even in complex systems found in chemistry and biology. In this minireview, we briefly outline the fundamental principles and experimental procedures of 2D IR spectroscopy. Using illustrative example studies, we explain the important features of 2D IR spectra and their capability to elucidate molecular structure and dynamics. Primarily, this minireview aims to convey the scope and potential of 2D IR spectroscopy by highlighting select examples of recent applications including the use of innate or introduced vibrational probes for the study of nucleic acids, peptides/proteins, and materials.

Quantitative proteomics reveals middle infrared radiation-interfered networks in breast cancer cells

Breast cancer is one of the leading cancer-related causes of death worldwide. Treatment of triple-negative breast cancer (TNBC) is complex and challenging, especially when metastasis has developed. In this study, we applied infrared radiation as an alternative approach for the treatment of TNBC. We used middle infrared (MIR) with a wavelength range of 3-5 μm to irradiate breast cancer cells. MIR significantly inhibited cell proliferation in several breast cancer cells but did not affect the growth of normal breast epithelial cells. We performed iTRAQ-coupled LC-MS/MS analysis to investigate the MIR-triggered molecular mechanisms in breast cancer cells. A total of 1749 proteins were identified, quantified, and subjected to functional enrichment analysis. From the constructed functionally enriched network, we confirmed that MIR caused G2/M cell cycle arrest, remodeled the microtubule network to an astral pole arrangement, altered the actin filament formation and focal adhesion molecule localization, and reduced cell migration activity and invasion ability. Our results reveal the coordinative effects of MIR-regulated physiological responses in concentrated networks, demonstrating the potential implementation of infrared radiation in breast cancer therapy.

Mode-specific vibrational relaxation of photoexcited guanosine 5'-monophosphate and its acid form: a femtosecond broadband mid-IR transient absorption and theoretical study

UV-pump/broadband-mid-IR-probe transient absorption (TA) experiments and ab initio quantum mechanical (QM) calculations were used to investigate the photophysics in heavy water of the neutral and acid forms of guanosine 5'-monophosphate (GMP and GMPD(+), respectively). Excited GMP undergoes ultrafast internal conversion (IC) and returns to the electronic ground state in less than one picosecond with a large amount of excess vibrational energy. The spectroscopic signals are dominated by vibrational cooling - a process in which the solute dissipates vibrational energy to the solvent. For neutral GMP, cooling proceeds with a time constant of 3 ps. Following IC, at least some medium-frequency modes such as the carbonyl stretch and an in-plane ring vibration are excited, suggesting that the vibrational energy distribution is non-statistical. This is consistent with predicted structural changes upon passage through the S1/S0 conical intersection. GMPD(+) differs from GMP by a single deuteron at the N7 position, but has a dramatically longer lifetime of 200 ps. Vibrational cooling of the S1 state of GMPD(+) was monitored via several medium-frequency modes that were assigned using QM calculations. These medium-frequency modes are also vibrationally excited in a non-statistical fashion. Excitation of these modes is in line with the change in geometry at the S1 minimum of GMPD(+) predicted by QM calculations. Furthermore, these modes relax at different rates, fully consistent with QM calculations, which predict that excited vibrational states of the carbonyl stretch couple strongly to the D2O solvent and thus deactivate via intermolecular energy transfer (IET). In contrast, the ring stretch couples strongly to other ring modes of the guanine chromophore and appears to decay via intramolecular vibrational energy redistribution (IVR).

N-H stretching modes of adenosine monomer in solution studied by ultrafast nonlinear infrared spectroscopy and ab initio calculations

The N-H stretching vibrations of adenine, one of the building blocks of DNA, are studied by combining infrared absorption and nonlinear two-dimensional infrared spectroscopy with ab initio calculations. We determine diagonal and off-diagonal anharmonicities of N-H stretching vibrations in chemically modified adenosine monomer dissolved in chloroform. For the single-quantum excitation manifold, the normal mode picture with symmetric and asymmetric NH(2) stretching vibrations is fully appropriate. For the two-quantum excitation manifold, however, the interplay between intermode coupling and frequency shifts due to a large diagonal anharmonicity leads to a situation where strong mixing does not occur. We compare our findings with previously reported values obtained on overtone spectroscopy of coupled hydrogen stretching oscillators.

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.

ULTRA: A Unique Instrument for Time-Resolved Spectroscopy

We report the development of a high-sensitivity time-resolved infrared and Raman spectrometer with exceptional experimental flexibility based on a 10-kHz synchronized dual-arm femtosecond and picosecond laser system. Ultrafast high-average-power titanium sapphire lasers and optical parametric amplifiers provide wavelength tuning from the ultraviolet (UV) to the mid-infrared region. Customized silicon, indium gallium arsenide, and mercury cadmium telluride linear array detectors are provided to monitor the probe laser intensity in the UV to mid-infrared wavelength range capable of measuring changes in sample absorbance of ΔOD ~ 10(-5) in 1 second. The system performance is demonstrated for the time-resolved infrared, two-dimensional (2D) infrared, and femtosecond stimulated Raman spectroscopy techniques with organometallic intermediates, organic excited states, and the dynamics of the tertiary structure of DNA.

The effect of water on infrared spectra of DNA

The effect of water on infrared (IR) spectra of DNA is studied by using first-principles calculations based on the hybrid density functional theory. Our calculations of frequencies of C = O stretching modes in DNA without water do not reproduce the IR spectra, whereas calculations for DNA surrounded by water reproduce the IR spectra well. We also find that the energy of adsorption of one water molecule around a C = O site is large ([Formula: see text] eV) owing to hydrogen bond formation. We therefore conclude that the effect of water plays an important role as regards IR spectra of DNA in a water solution.

Structural and optical properties of DNA layers covalently attached to diamond surfaces

Label-free detection of DNA molecules on chemically vapor-deposited diamond surfaces is achieved with spectroscopic ellipsometry in the infrared and vacuum ultraviolet range. This nondestructive method has the potential to yield information on the average orientation of single as well as double-stranded DNA molecules, without restricting the strand length to the persistence length. The orientational analysis based on electronic excitations in combination with information from layer thicknesses provides a deeper understanding of biological layers on diamond. The pi-pi* transition dipole moments, corresponding to a transition at 4.74 eV, originate from the individual bases. They are in a plane perpendicular to the DNA backbone with an associated n-pi* transition at 4.47 eV. For 8-36 bases of single- and double-stranded DNA covalently attached to ultra-nanocrystalline diamond, the ratio between in- and out-of-plane components in the best fit simulations to the ellipsometric spectra yields an average tilt angle of the DNA backbone with respect to the surface plane ranging from 45 degrees to 52 degrees . We comment on the physical meaning of the calculated tilt angles. Additional information is gathered from atomic force microscopy, fluorescence imaging, and wetting experiments. The results reported here are of value in understanding and optimizing the performance of the electronic readout of a diamond-based label-free DNA hybridization sensor.

High-resolution waveguide THz spectroscopy of biological molecules

Low-frequency vibrational modes of biological molecules consist of intramolecular modes, which are dependent on the molecule as a whole, as well as intermolecular modes, which arise from hydrogen-bonding interactions and van der Waals forces. Vibrational modes thus contain important information about conformation dynamics of biological molecules, and can also be used for identification purposes. However, conventional Fourier transform infrared spectroscopy and terahertz time-domain spectroscopy (THz-TDS) often result in broad, overlapping features that are difficult to distinguish. The technique of waveguide THz-TDS has been recently developed, resulting in sharper features. For this technique, an ordered polycrystalline film of the molecule is formed on a metal sample plate. This plate is incorporated into a metal parallel-plate waveguide and probed via waveguide THz-TDS. The planar order of the film reduces the inhomogeneous broadening, and cooling of the samples to 77K reduces the homogenous broadening. This combination results in the line-narrowing of THz vibrational modes, in some cases to an unprecedented degree. Here, this technique has been demonstrated with seven small biological molecules, thymine, deoxycytidine, adenosine, D-glucose, tryptophan, glycine, and L-alanine. The successful demonstration of this technique shows the possibilities and promise for future studies of internal vibrational modes of large biological molecules.

Vibrational dynamics of DNA. III. Molecular dynamics simulations of DNA in water and theoretical calculations of the two-dimensional vibrational spectra

A theoretical description of the vibrational excitons in DNA is presented by using the vibrational basis mode theory developed in Papers I and II. The parameters obtained from the density functional theory calculations, such as vibrational coupling constants and basis mode frequencies, are used to numerically simulate two-dimensional (2D) IR spectra of dG(n):dC(n) and dA(n):dT(n) double helices with n varying from 1 to 10. From the molecular dynamics simulations of dG(5)C(5) and dA(5)T(5) double helices in D(2)O solution, it is found that the thermally driven internal motions of these systems in an aqueous solution do not induce strong fluctuations of basis mode frequencies nor vibrational couplings. In order to construct the two-exciton Hamiltonian, the vibrational anharmonicities of eight basis modes are obtained by carrying out B3LYP6-31G(*) calculations for the nine basis modes. The simulated 2D IR spectra of dG(n):dC(n) double helix in D(2)O solution are directly compared with closely related experimental results. The 2D IR spectra of dG(n):dC(n) and dA(n):dT(n) are found to be weakly dependent on the number of base pairs. The present work demonstrates that the computational procedure combining quantum chemistry calculation and molecular dynamics simulation methods can be of use to predict 2D IR spectra of nucleic acids in solutions.

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.