Ultrafast two-dimensional infrared spectroscopy of guanine-cytosine base pairs in DNA oligomers

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.

Probing the local structure and dynamics of nucleotides using vibrationally enhanced alkynyl stretching

Monitoring the site-specific local structure and dynamics of polynucleotides and DNA is important for understanding their biological functions. However, structurally characterizing these biomolecules with high time resolution has been known to be experimentally challenging. In this work, several 5-silylethynyl-2'-deoxynucleosides and 5-substituted phenylethynyl-2'-deoxynucleosides on the basis of deoxycytidine (dC) and deoxythymidine (dT) were synthesized, in which the alkynyl group shows intensified CC stretching vibration with infrared transition dipole moment magnitude close to that of typical CO stretching, and exhibits structural sensitivities in both vibrational frequency and spectral width. In particular, 5-trimethylsilylethynyl-2'-dC (^TMSEdC, molecule 1a) was examined in detail using femtosecond nonlinear IR spectroscopy. The solvent dependent CC stretching frequency of 1a can be reasonably interpreted mainly as the hydrogen-bonding effect between the solvent and cytosine base ring structure. Transient 2D IR and pump-probe IR measurements of 1a carried out comparatively in two aprotic solvents (DMSO and THF) and one protic solvent (MeOH) further reveal solvent dependent ultrafast vibrational properties, including diagonal anharmonicity, spectral diffusion, vibrational relaxation and anisotropy dynamics. These observed sensitivities are rooted in an extended π-conjugation of the base ring structure in which the CC group is actively involved. Our results show that the intensified CC stretching vibration can potentially provide a site-specific IR probe for monitoring the equilibrium and ultrafast structural dynamics of polynucleotides.

Vibrational Spectra of Nucleotides in the Presence of the Au Cluster Enhancer in MD Simulation of a SERS Sensor

Surface-enhanced Raman scattering (SERS) nanoprobes have shown tremendous potential in in vivo imaging. The development of single oligomer resolution in the SERS promotes experiments on DNA and protein identification using SERS as a nanobiosensor. As Raman scanners rely on a multiple spectrum acquisition, faster imaging in real-time is required. SERS weak signal requires averaging of the acquired spectra that erases information on conformation and interaction. To build spectral libraries, the simulation of measurement conditions and conformational variations for the nucleotides relative to enhancer nanostructures would be desirable. In the molecular dynamic (MD) model of a sensing system, we simulate vibrational spectra of the cytosine nucleotide in FF2/FF3 potential in the dynamic interaction with the Au20 nanoparticles (NP) (EAM potential). Fourier transfer of the density of states (DOS) was performed to obtain the spectra of bonds in reaction coordinates for nucleotides at a resolution of 20 to 40 cm-1. The Au20 was optimized by ab initio density functional theory with generalized gradient approximation (DFT GGA) and relaxed by MD. The optimal localization of nucleotide vs. NP was defined and the spectral modes of both components vs. interaction studied. Bond-dependent spectral maps of nucleotide and NP have shown response to interaction. The marker frequencies of the Au20-nucleotide interaction have been evaluated.

Infrared Spectroscopy Reveals the Preferred Motif Size and Local Disorder in Parallel Stranded DNA G-Quadruplexes

Infrared spectroscopy detects the formation of G-quadruplexes in guanine-rich nucleic acid sequences through shifts in the guanine C=O stretch mode. Here, we use ultrafast 2D infrared (IR) spectroscopy and isotope substitution to show that these shifts arise from vibrational delocalization among stacked G-quartets. This provides a direct measure of the sizes of locally ordered motifs in heterogeneous samples with substantial disordered regions. We find that parallel-stranded, potassium-bound DNA G-quadruplexes are limited to five consecutive G-quartets and 3-4 consecutive layers are preferred for longer polyguanine tracts. The resulting potassium-dependent G-quadruplex assembly landscape reflects the polyguanine tract lengths found in genomes, the ionic conditions prevalent in healthy mammalian cells, and the onset of structural disorder in disease states. Our study describes spectral markers that can be used to probe other G-quadruplex structures and provides insight into the fundamental limits of their formation in biological and artificial systems.

Vibrations of the guanine-cytosine pair in chloroform: an anharmonic computational study

We compute at the anharmonic level the vibrational spectra of the Watson-Crick dimer formed by guanosine (G) and cytidine (C) in chloroform, together with those of G, C and the most populated GG dimer. The spectra for deuterated and partially deuterated GC are also computed. We use DFT calculations, with B3LYP and CAM-B3LYP as reference functionals. Solvent effects from chloroform are included via the Polarizable Continuum Model (PCM), and by performing tests on models including up two chloroform molecules. Both B3LYP and CAM-B3LYP calculations reproduce the shape of the experimental spectra well in the fingerprint region (1500-1700 cm-1) and in the N-H stretching region (2800-3600 cm-1), with B3LYP providing better quantitative agreement with experiments. According to our calculations, the N-H amido streching mode of G falls at ∼2900 cm-1, while the N-H amino of G and C falls at ∼3100 cm-1 when hydrogen-bonded, or ∼3500 cm-1 when free. Overtone and combination bands strongly contribute to the absorption band at ∼3300 cm-1. Inclusion of bulk solvent effects significantly increases the accuracy of the computed spectra, while solute-solvent interactions have a smaller, though still noticeable, effect. Some key aspects of the anharmonic treatment of strongly vibrationally coupled supermolecular systems and the related methodological issues are also discussed.

Watching Proteins Wiggle: Mapping Structures with Two-Dimensional Infrared Spectroscopy

Proteins exhibit structural fluctuations over decades of time scales. From the picosecond side chain motions to aggregates that form over the course of minutes, characterizing protein structure over these vast lengths of time is important to understanding their function. In the past 15 years, two-dimensional infrared spectroscopy (2D IR) has been established as a versatile tool that can uniquely probe proteins structures on many time scales. In this review, we present some of the basic principles behind 2D IR and show how they have, and can, impact the field of protein biophysics. We highlight experiments in which 2D IR spectroscopy has provided structural and dynamical data that would be difficult to obtain with more standard structural biology techniques. We also highlight technological developments in 2D IR that continue to expand the scope of scientific problems that can be accessed in the biomedical sciences.

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.

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.

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.

Proton transfer rates in ionized water clusters (H2O)n (n = 2–4)

A proton transfer process is usually dominant in several biological phenomena such as the energy relaxation of photo-excited DNA base pairs and a charge relay process in Ser-His-Glu.

On the absence of intrahelical DNA dynamics on the μs to ms timescale

DNA helices display a rich tapestry of motion on both short (<100 ns) and long (>1 ms) timescales. However, with the exception of mismatched or damaged DNA, experimental measures indicate that motions in the 1 μs to 1 ms range are effectively absent, which is often attributed to difficulties in measuring motions in this time range. We hypothesized that these motions have not been measured because there is effectively no motion on this timescale, as this provides a means to distinguish faithful Watson-Crick base-paired DNA from damaged DNA. The absence of motion on this timescale would present a 'static' DNA sequence-specific structure that matches the encounter timescales of proteins, thereby facilitating recognition. Here we report long-timescale (~10-44 μs) molecular dynamics simulations of a B-DNA duplex structure that addresses this hypothesis using both an 'Anton' machine and large ensembles of AMBER GPU simulations.