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

Breaking Barriers in Ultrafast Spectroscopy and Imaging Using 100 kHz Amplified Yb-Laser Systems

ConspectusUltrafast spectroscopy and imaging have become tools utilized by a broad range of scientists involved in materials, energy, biological, and chemical sciences. Commercialization of ultrafast spectrometers including transient absorption spectrometers, vibrational sum frequency generation spectrometers, and even multidimensional spectrometers have put these advanced spectroscopy measurements into the hands of practitioners originally outside the field of ultrafast spectroscopy. There is now a technology shift occurring in ultrafast spectroscopy, made possible by new Yb-based lasers, that is opening exciting new experiments in the chemical and physical sciences. Amplified Yb-based lasers are not only more compact and efficient than their predecessors but also, most importantly, operate at many times the repetition rate with improved noise characteristics in comparison to the previous generation of Ti:sapphire amplifier technologies. Taken together, these attributes are enabling new experiments, generating improvements to long-standing techniques, and affording the transformation of spectroscopies to microscopies. This Account aims to show that the shift to 100 kHz lasers is a transformative step in nonlinear spectroscopy and imaging, much like the dramatic expansion that occurred with the commercialization of Ti:sapphire laser systems in the 1990s. The impact of this technology will be felt across a great swath of scientific communities. We first describe the technology landscape of amplified Yb-based laser systems used in conjunction with 100 kHz spectrometers operating with shot-to-shot pulse shaping and detection. We also identify the range of different parametric conversion and supercontinuum techniques which now provide a path to making pulses of light optimal for ultrafast spectroscopy. Second, we describe specific instances from our laboratories of how the amplified Yb-based light sources and spectrometers are transformative. For multiple probe time-resolved infrared and transient 2D IR spectroscopy, the gain in temporal span and signal-to-noise enables dynamical spectroscopy measurements from femtoseconds to seconds. These gains widen the applicability of time-resolved infrared techniques across a range of topics in photochemistry, photocatalysis, and photobiology as well as lower the technical barriers to implementation in a laboratory. For 2D visible spectroscopy and microscopy with white light, as well as 2D IR imaging, the high repetition rates of these new Yb-based light sources allow one to spatially map 2D spectra while maintaining high signal-to-noise in the data. To illustrate the gains, we provide examples of imaging applications in the study of photovoltaic materials and spectroelectrochemistry.

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.

The mutual interactions of RNA, counterions and water - quantifying the electrostatics at the phosphate-water interface

The structure and dynamics of polyanionic biomolecules, like RNA, are decisively determined by their electric interactions with the water molecules and the counterions in the environment. The solvation dynamics of the biomolecules involves a subtle balance of non-covalent and many-body interactions with structural fluctuations due to thermal motion occurring in a femto- to subnanosecond time range. This complex fluctuating many particle scenario is crucial in defining the properties of biological interfaces with far reaching significance for the folding of RNA structures and for facilitating RNA-protein interactions. Given the inherent complexity, suited model systems, carefully calibrated and benchmarked by experiments, are required to quantify the relevant interactions of RNA with the aqueous environment. In this feature article we summarize our recent progress in the understanding of the electrostatics at the biological interface of double stranded RNA (dsRNA) and transfer RNA (tRNA). Dimethyl phosphate (DMP) is introduced as a viable and rigorously accessible model system allowing the interaction strength with water molecules and counterions, their relevant fluctuation timescales and the spatial reach of interactions to be established. We find strong (up to ≈90 MV cm-1) interfacial electric fields with fluctuations extending up to ≈20 THz and demonstrate how the asymmetric stretching vibration νAS(PO2)- of the polarizable phosphate group can serve as the most sensitive probe for interfacial interactions, establishing a rigorous link between simulations and experiment. The approach allows for the direct interfacial observation of interactions of biologically relevant Mg2+ counterions with phosphate groups in contact pair geometries via the rise of a new absorption band imposed by exchange repulsion interactions at short interatomic distances. The systematic extension to RNA provides microscopic insights into the changes of the hydration structure that accompany the temperature induced melting of the dsRNA double helix and quantify the ionic interactions in the folded tRNA. The results show that pairs of negatively charged phosphate groups and Mg2+ ions represent a key structural feature of RNA embedded in water. They highlight the importance of binding motifs made of contact pairs in the electrostatic stabilization of RNA structures that have a strong impact on the surface potential and enable the fine tuning of the local electrostatic properties which are expected to be relevant for mediating the interactions between biomolecules.

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.

Phosphate Vibrations Probe Electric Fields in Hydrated Biomolecules: Spectroscopy, Dynamics, and Interactions

Electric interactions have a strong impact on the structure and dynamics of biomolecules in their native water environment. Given the variety of water arrangements in hydration shells and the femto- to subnanosecond time range of structural fluctuations, there is a strong quest for sensitive noninvasive probes of local electric fields. The stretching vibrations of phosphate groups, in particular the asymmetric (PO₂)⁻ stretching vibration ν_AS(PO₂)⁻, allow for a quantitative mapping of dynamic electric fields in aqueous environments via a field-induced redshift of their transition frequencies and concomitant changes of vibrational line shapes. We present a systematic study of ν_AS(PO₂)⁻ excitations in molecular systems of increasing complexity, including dimethyl phosphate (DMP), short DNA and RNA duplex structures, and transfer RNA (tRNA) in water. A combination of linear infrared absorption, two-dimensional infrared (2D-IR) spectroscopy, and molecular dynamics (MD) simulations gives quantitative insight in electric-field tuning rates of vibrational frequencies, electric field and fluctuation amplitudes, and molecular interaction geometries. Beyond neat water environments, the formation of contact ion pairs of phosphate groups with Mg²⁺ ions is demonstrated via frequency upshifts of the ν_AS(PO₂)⁻ vibration, resulting in a distinct vibrational band. The frequency positions of contact geometries are determined by an interplay of attractive electric and repulsive exchange interactions.

Change of Hydration Patterns upon RNA Melting Probed by Excitations of Phosphate Backbone Vibrations

The water hydration shell has a decisive impact on the structural and functional properties of RNA. Changes of the RNA structure upon melting and in biochemical processes are accompanied by a change of hydration patterns, a process that is barely characterized. To discern hydration geometries around the backbone phosphate groups of an RNA double helix at the molecular level, we combine two-dimensional infrared spectroscopy of phosphate vibrations with theoretical simulations. There are three distinct coexisting hydration motifs of the RNA A-helix: an ordered chain-like arrangement of water molecules with links between neighboring phosphate groups, separate local hydration shells of up to six water molecules, and hydrated phosphate/counterion contact pairs. RNA disordering upon melting is connected with a transition from predominantly ordered water structures to local hydration shells around phosphate units. Structural fluctuations are dominated by librational water motions occurring on a 300 fs time scale, without exchange between hydration motifs.

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.

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.

Central-metal effect on intramolecular vibrational energy transfer of M(CO)5Br (M = Mn, Re) probed by two-dimensional infrared spectroscopy

Vibrational energy transfer in transition metal complexes with flexible structures in condensed phases is of central importance to catalytical chemistry processes. In this work, two molecules with different metal atoms, M(CO)₅Br (where M = Mn, Re), were used as model systems, and their axial and radial carbonyl stretching modes as infrared probes. The central-metal effect on intramolecular vibrational energy redistribution (IVR) in M(CO)₅Br was investigated in polar and nonpolar solvents. The linear infrared (IR) peak splitting between carbonyl vibrations increases as the metal atom changes from Mn to Re. The waiting-time dependent two-dimensional infrared diagonal- and off-diagonal peak amplitudes reveal a faster IVR process in Re(CO)₅Br than in Mn(CO)₅Br. With the aid of density functional theory (DFT) calculations, the central-metal effect on IVR time linearly correlates with the vibrational coupling strength between the two involved modes. In addition, the polar solvent is found to accelerate the IVR process by affecting the anharmonic vibrational potentials of a solute vibration mode.

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.