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

Interaction-Dependent Secondary Structure of Peptides in Biomolecular Condensates

Biomolecular condensates provide a mechanism for compartmentalization of biomolecules in eukaryotic cells. These liquid-like condensates are formed via liquid-liquid phase separation, by a plethora of interactions, and can mediate several biological processes in healthy cells. Expansions of dipeptide repeat proteins, DPRs, in which arginine rich DPRs like poly-proline-arginine (PR), and poly-glycine-arginine (GR), partition RNA into condensates can however induce cell toxicity. Here, we use (GR)20 as a model for biological poly-GR and condense it using either excluded volume interactions with polyethylene glycol (PEG) as a crowder or direct electrostatic interactions with RNA oligomers. Using two-dimensional infrared (2D IR) spectroscopy, we observe that (GR)20 condensed through an excluded volume forms β-sheet structures, whereas (GR)20 condensed with RNA forms loops. We also investigate local hydrogen-bond dynamics in the condensate and compare the measurements with molecular dynamics simulations. Hydrogen bond lifetimes undergo a marked slowdown compared to dynamics in the dilute phase. This is representative of confined water within the percolated networks inside the condensate due to the interaction present in the condensate disrupting H-bond networks. Overall, our results show that both protein structure and dynamics are inherently dependent on the type of interactions that stabilize the condensates.

Modeling Infrared Spectroscopy of Nucleic Acids: Integrating Vibrational Non-Condon Effects with Machine Learning Schemes

Vibrational non-Condon effects, which describe how molecular vibrational transitions are influenced by a system's rotational and translational degrees of freedom, are often overlooked in spectroscopy studies of biological macromolecules. In this work, we explore these effects in the modeling of infrared (IR) spectra for nucleic acids in the 1600-1800 cm-1 region. Through electronic structure calculations, we reveal that the transition dipole moments of the C═O and C═C stretching modes in nucleobases are highly sensitive to solvation, hydrogen bonding, and base stacking conditions. To incorporate vibrational non-Condon effects into spectroscopy modeling, we use local electric fields on chromophore atoms as collective coordinates and leverage experimental IR spectra of oligonucleotides to develop deep neural network-based transition dipole strength (TDS) maps for the C═O and C═C chromophores. By integrating molecular dynamics simulations with a mixed quantum/classical treatment of the line shape theory, we apply the TDS maps to calculate the IR spectra of nucleoside 5'-monophosphates, DNA double helices and yeast phenylalanine tRNA. The resulting theoretical spectra show quantitative agreement with experimental measurements. While the predictions for nucleoside 5'-monophosphates are comparable to baseline performance, the TDS maps yield significantly improved IR peak intensities across all oligonucleotides. This theoretical framework effectively bridges atomistic simulations and IR spectroscopy experiments, offering molecular insights into how vibrational non-Condon effects impact the observed spectral features.

NAIRDB: a database of Fourier transform infrared (FTIR) data for nucleic acids

The Nucleic Acid InfraRed Data Bank (NAIRDB) serves as a comprehensive public repository dedicated to the archival and free distribution of Fourier transform infrared (FTIR) spectral data specific to nucleic acids. This database encompasses a collection of FTIR spectra covering diverse nucleic acid molecules, including DNA, RNA, DNA/RNA hybrids and their various derivatives. NAIRDB covers details of the experimental conditions for FTIR measurements, literature links, primary sequence data, information about experimentally determined structures for related nucleic acid molecules and/or computationally modeled 3D structures. All entries undergo expert validation and curation to maintain the completeness, consistency and quality of the data. NAIRDB can be searched by similarity of nucleic acid sequences or by direct comparison of spectra. The database is open for the submission of the FTIR data for nucleic acids. NAIRDB is available at https://nairdb.genesilico.pl.

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.

Unraveling the Structure-Spectrum Relationship of Yeast Phenylalanine Transfer RNA: Insights from Theoretical Modeling of Infrared Spectroscopy

Yeast phenylalanine tRNA (tRNAphe) is a paradigmatic model in structural biology. In this work, we combine molecular dynamics simulations and spectroscopy modeling to establish a direct link between its structure, conformational dynamics, and infrared (IR) spectra. Employing recently developed vibrational frequency maps and coupling models, we apply a mixed quantum/classical treatment of the line shape theory to simulate the IR spectra of tRNAphe in the 1600-1800 cm-1 region across its folded and unfolded conformations and under varying concentrations of Mg2+ ions. The predicted IR spectra of folded and unfolded tRNAphe are in good agreement with experimental measurements, validating our theoretical framework. We then elucidate how the characteristic L-shaped tertiary structure of the tRNA and its modulation in response to diverse chemical environments give rise to distinct IR absorption peaks and line shapes. These calculations effectively bridge IR spectroscopy experiments and atomistic molecular simulations, unraveling the molecular origins of the observed IR spectra of tRNAphe. This work presents a robust theoretical protocol for modeling the IR spectroscopy of nucleic acids, which will facilitate its application as a sensitive probe for detecting the fluctuating secondary and tertiary structures of these essential biological macromolecules.

Structures and dynamics of 8-oxo-7,8-dihydro-2'-deoxyguanosine in neutral and basic aqueous solutions by spectroscopy

8-oxo-7,8-dihydro-2'-dexyoguanine (8-oxo-dG) can be tautomerized to a 6-enolate,8-keto tautomer through nearby-NH deprotonation at elevated pH. In this work, the N3-protonated 8-oxo-dG tautomers in deuterated pH-buffer solutions were studied using steady-state UV/Vis, FTIR, and ultrafast two-dimensional IR spectroscopies. The presence of 6,8-diketo and C6-anionic tautomers at neutral to basic conditions (pD = 7.4-12.0) was revealed by UV/Vis and FTIR results and was further confirmed by 2D IR signals in both diagonal and off-diagonal regions. However, the C6-enol tautomer, which may be an intermediate during the transition from 6,8-diketo to C6-enolate,C8-keto, was not observed appreciably due to its extreme low population. Furthermore, the neutral-to-anionic tautomeric transition of N3H-8-oxo-dG studied in this work occurs under more basic conditions than the N1H-8-oxo-dG reported previously, showing a higher pKa value for N3H than N1H. Finally, vibrational relaxation of the carbonyl stretching mode was found to be both molecular site dependent and pD dependent for 8oxo-dG. Taken together, this work shows that the ultrafast infrared spectroscopic method is effective for examining tautomers and their dynamics in nucleic acids.

Isotope-Edited Variable Temperature Infrared Spectroscopy for Measuring Transition Temperatures of Single A-T Watson-Crick Base Pairs in DNA Duplexes

Experimental methods to determine transition temperatures for individual base pair melting events in DNA duplexes are lacking despite intense interest in these thermodynamic parameters. Here, we determine the dimensions of the thymine (T) C2═O stretching vibration when it is within the DNA duplex via isotopic substitutions at other atomic positions in the structure. First, we determined that this stretching state was localized enough to specific atoms in the molecule to make submolecular scale measurements of local structure and stability in high molecular weight complexes. Next, we develop a new isotope-edited variable temperature infrared method to measure melting transitions at various locations in a DNA structure. As an initial test of this "sub-molecular scale thermometer", we applied our T13C2 difference infrared signal to measure location-dependent melting temperatures (TmL) in a DNA duplex via variable temperature attenuated total reflectance Fourier transform infrared (VT-ATR-FTIR) spectroscopy. We report that the TmL of a single Watson-Crick A-T base pair near the end of an A-T rich sequence (poly T) is ∼34.9 ± 0.7°C. This is slightly lower than the TmL of a single base pair near the middle position of the poly T sequence (TmL ∼35.6±0.2°C). In addition, we also report that the TmL of a single Watson-Crick A-T base pair near the end of a 50% G-C sequence (12-mer) is ∼52.5 ± 0.3°C, which is slightly lower than the global melting Tm of the 12-mer sequence (TmL ∼54.0±0.9°C). Our results provide direct physical evidence for end fraying in DNA sequences with our novel spectroscopic methods.

Triplet Excimer Formation in a DNA Duplex with Silver Ion-Mediated Base Pairs

The dynamics of excited electronic states in self-assembled structures formed between silver(I) ions and cytosine-containing DNA strands or monomeric cytosine derivatives were investigated by time-resolved infrared (TRIR) spectroscopy and quantum mechanical calculations. The steady-state and time-resolved spectra depend sensitively on the underlying structures, which change with pH and the nucleobase and silver ion concentrations. At pH ∼ 4 and low dC20 strand concentration, an intramolecularly folded i-motif is observed, in which protons, and not silver ions, mediate C-C base pairing. However, at the higher strand concentrations used in the TRIR measurements, dC20 strands associate pairwise to yield duplex structures containing C-Ag+-C base pairs with a high degree of propeller twisting. UV excitation of the silver ion-mediated duplex produces a long-lived excited state, which we assign to a triplet excimer state localized on a pair of stacked cytosines. The computational results indicate that the propeller-twisted motifs induced by metal-ion binding are responsible for the enhanced intersystem crossing that populates the triplet state and not a generic heavy atom effect. Although triplet excimer states have been discussed frequently as intermediates in the formation of cyclobutane pyrimidine dimers, we find neither computational nor experimental evidence for cytosine-cytosine photoproduct formation in the systems studied. These findings provide a rare demonstration of a long-lived triplet excited state that is formed in a significant yield in a DNA duplex, demonstrating that supramolecular structural changes induced by metal ion binding profoundly affect DNA photophysics.

Biomolecular infrared spectroscopy: making time for dynamics

Time resolved infrared spectroscopy of biological molecules has provided a wealth of information relating to structural dynamics, conformational changes, solvation and intermolecular interactions. Challenges still exist however arising from the wide range of timescales over which biological processes occur, stretching from picoseconds to minutes or hours. Experimental methods are often limited by vibrational lifetimes of probe groups, which are typically on the order of picoseconds, while measuring an evolving system continuously over some 18 orders of magnitude in time presents a raft of technological hurdles. In this Perspective, a series of recent advances which allow biological molecules and processes to be studied over an increasing range of timescales, while maintaining ultrafast time resolution, will be reviewed, showing that the potential for real-time observation of biomolecular function draws ever closer, while offering a new set of challenges to be overcome.

Ultrafast Spectroscopy Reveals Slow Water Dynamics in Biocondensates

Cells achieve high spatiotemporal control over biochemical processes through compartmentalization to membrane-bound as well as membraneless organelles that assemble by liquid-liquid phase separation. Characterizing the balance of forces within these environments is essential to understanding their stability and function, and water is an integral part of the condensate, playing an important role in mediating electrostatic and hydrogen-bonding interactions. Here, we investigate the ultrafast, picosecond hydrogen-bond dynamics of a model biocondensate consisting of a peptide poly-l-arginine (Poly-R) and the nucleic acid adenosine monophosphate (AMP) using coherent two-dimensional infrared (2D IR) spectroscopy. We investigated three vibrational modes: the arginine side-chain C═N stretches, an AMP ring mode, and the amide backbone carbonyl stretching modes. Dynamics slow considerably between the dilute phase and the condensate phase for each vibrational probe. For example, the arginine side-chain C═N modes slow from 0.38 to 2.26 ps due to strong electrostatic interactions. All-atom molecular dynamics simulations provide an atomistic interpretation of the H-bond network disruption resulting from electrostatic contributions as well as collapse within the condensate. Simulations predict that a fraction of water molecules are highly constrained within the condensate, explaining the observed slowdown in the H-bond dynamics.

Unraveling Excited State Dynamics of a Single-Stranded DNA-Assembled Conjugated Polyelectrolyte

Conformational templating of conjugated polyelectrolytes with single-stranded DNAs (ssDNAs) has the prospect of tailoring excited state dynamics for specific optoelectronic applications. We use ultrafast time-resolved infrared spectroscopy to study the photophysics of a cationic polythiophene assembled with different ssDNAs, inducing distinct conformations (flexible disordered structures vs more rigid complexes with increased backbone planarity). Intrachain polarons are always produced upon selective excitation of the polymer, the extent being dependent on backbone torsional order. Polaron formation and decay were monitored through evolution of IR-active vibrational modes that interfere with mid-IR polaron electronic absorption giving rise to Fano-antiresonances. Selective UV excitation of ssDNAs revealed that stacking interactions between thiophene rings and nucleic acid bases can promote the formation of an intermolecular charge transfer complex. The findings inform designers of functional conjugated polymers by identifying that involvement of the scaffold in the photophysics needs to be considered when developing such structures for optoelectronic applications.

Thermodynamics and kinetics of DNA and RNA dinucleotide hybridization to gaps and overhangs

Hybridization of short nucleic acid segments (<4 nt) to single-strand templates occurs as a critical intermediate in processes such as nonenzymatic nucleic acid replication and toehold-mediated strand displacement. These templates often contain adjacent duplex segments that stabilize base pairing with single-strand gaps or overhangs, but the thermodynamics and kinetics of hybridization in such contexts are poorly understood because of the experimental challenges of probing weak binding and rapid structural dynamics. Here we develop an approach to directly measure the thermodynamics and kinetics of DNA and RNA dinucleotide dehybridization using steady-state and temperature-jump infrared spectroscopy. Our results suggest that dinucleotide binding is stabilized through coaxial stacking interactions with the adjacent duplex segments as well as from potential noncanonical base-pairing configurations and structural dynamics of gap and overhang templates revealed using molecular dynamics simulations. We measure timescales for dissociation ranging from 0.2-40 μs depending on the template and temperature. Dinucleotide hybridization and dehybridization involve a significant free energy barrier with characteristics resembling that of canonical oligonucleotides. Together, our work provides an initial step for predicting the stability and kinetics of hybridization between short nucleic acid segments and various templates.

Ultrafast transient vibrational action spectroscopy of cryogenically cooled ions

Ultrafast transient vibrational action spectra of cryogenically cooled Re(CO)3(CH3CN)3+ ions are presented. Nonlinear spectra were collected in the time domain by monitoring the photodissociation of a weakly bound N2 messenger tag as a function of delay times and phases between a set of three infrared pulses. Frequency-resolved spectra in the carbonyl stretch region show relatively strong bleaching signals that oscillate at the difference frequency between the two observed vibrational features as a function of the pump-probe waiting time. This observation is consistent with the presence of nonlinear pathways resulting from underlying cross-peak signals between the coupled symmetric-asymmetric C≡O stretch pair. The successful demonstration of frequency-resolved ultrafast transient vibrational action spectroscopy of dilute molecular ion ensembles provides an exciting, new framework for the study of molecular dynamics in isolated, complex molecular ion systems.

Direct monitoring of the thermodynamics and kinetics of DNA and RNA dinucleotide dehybridization from gaps and overhangs

Hybridization of short nucleic acid segments (<4 nucleotides) to single-strand templates occurs as a critical intermediate in processes such as non-enzymatic nucleic acid replication and toehold-mediated strand displacement. These templates often contain adjacent duplex segments that stabilize base pairing with single-strand gaps or overhangs, but the thermodynamics and kinetics of hybridization in such contexts are poorly understood due to experimental challenges of probing weak binding and rapid structural dynamics. Here we develop an approach to directly measure the thermodynamics and kinetics of DNA and RNA dinucleotide dehybridization using steady-state and temperature-jump infrared spectroscopy. Our results suggest that dinucleotide binding is stabilized through coaxial stacking interactions with the adjacent duplex segments as well as from potential non-canonical base pairing configurations and structural dynamics of gap and overhang templates revealed using molecular dynamics simulations. We measure timescales for dissociation ranging from 0.2 to 40 µs depending on the template and temperature. Dinucleotide hybridization and dehybridization involves a significant free energy barrier with characteristics resembling that of canonical oligonucleotides. Together, our work provides an initial step for predicting the stability and kinetics of hybridization between short nucleic acid segments and various templates.

Disruption of energetic and dynamic base pairing cooperativity in DNA duplexes by an abasic site

DNA duplex stability arises from cooperative interactions between multiple adjacent nucleotides that favor base pairing and stacking when formed as a continuous stretch rather than individually. Lesions and nucleobase modifications alter this stability in complex manners that remain challenging to understand despite their centrality to biology. Here, we investigate how an abasic site destabilizes small DNA duplexes and reshapes base pairing dynamics and hybridization pathways using temperature-jump infrared spectroscopy and coarse-grained molecular dynamics simulations. We show how an abasic site splits the cooperativity in a short duplex into two segments, which destabilizes small duplexes as a whole and enables metastable half-dissociated configurations. Dynamically, it introduces an additional barrier to hybridization by constraining the hybridization mechanism to a step-wise process of nucleating and zipping a stretch on one side of the abasic site and then the other.

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.

Vibrational relaxation by methylated xanthines in solution: Insights from 2D IR spectroscopy and calculations

Two-dimensional infrared (2D IR) spectroscopy, infrared pump-infrared probe spectroscopy, and density functional theory calculations were used to study vibrational relaxation by ring and carbonyl stretching modes in a series of methylated xanthine derivatives in acetonitrile and deuterium oxide (heavy water). Isotropic signals from the excited symmetric and asymmetric carbonyl stretch modes decay biexponentially in both solvents. Coherent energy transfer between the symmetric and asymmetric carbonyl stretching modes gives rise to a quantum beat in the time-dependent anisotropy signals. The damping time of the coherent oscillation agrees with the fast decay component of the carbonyl bleach recovery signals, indicating that this time constant reflects intramolecular vibrational redistribution (IVR) to other solute modes. Despite their similar frequencies, the excited ring modes decay monoexponentially with a time constant that matches the slow decay component of the carbonyl modes. The slow decay times, which are faster in heavy water than in acetonitrile, approximately match the ones observed in previous UV pump-IR probe measurements on the same compounds. The slow component is assigned to intermolecular energy transfer to solvent bath modes from low-frequency solute modes, which are populated by IVR and are anharmonically coupled to the carbonyl and ring stretch modes. 2D IR measurements indicate that the carbonyl stretching modes are weakly coupled to the delocalized ring modes, resulting in slow exchange that cannot explain the common solvent-dependence. IVR is suggested to occur at different rates for the carbonyl vs ring modes due to differences in mode-specific couplings and not to differences in the density of accessible states.

Excellent removal performance of 4,4'-biphenyldicarboxaldehyde m-phenylenediamine Schiff base magnetic polymer towards phenanthrene and 9-phenanthrol: Experimental, modeling and DFT calculations studies

Phenanthrene (PTH) and 9-phenanthrol (9-PTH) exhibited severe health threats and ecological hazards, for this reason, exploring a high-efficient removing strategy for PTH and 9-PTH could be considered of great urgency. Herein the 4,4'-biphenyldicarboxaldehyde m-phenylenediamine Schiff base magnetic polymer (magnetic BIPH-PHEN) was successfully fabricated via Schiff base polycondensation reaction and the subsequently one-pot embedded method. The mutual aromatic nucleus of BIPH-PHEN polymer and PTH/9-PTH could form π-π interaction, thus improving the capture ability, the embedded Fe₃O₄ nanoparticles provided the possibility for rapid separation. The physical and chemical properties of the magnetic BIPH-PHEN were systematically characterized. The removal rate of magnetic BIPH-PHEN towards PTH and 9-PTH was 85.65 % and 98.52 %, respectively (PTH or 9-PTH: 8 mg/L; Adsorbent: 0.2 g/L). The DFT calculations including energy calculations and electrostatic potential distribution analyzed the different bonding modes and proposed the most possible bonding modes in the adsorbent/adsorbate system. Moreover, the LUMO and HOMO orbits combined with energy gaps analysis proved the existence and specific types of the π-π interaction. The monolayer adsorption occurred on the homogeneous magnetic BIPH-PHEN surface, simultaneously the chemisorption was dominant. This work not only proposed new sights on assembling magnetic Schiff base polymer for removing polycyclic aromatic hydrocarbons, but also provided a deeper understanding of intramolecular interactions.

Quantum-chemical calculation of two-dimensional infrared spectra using localized-mode VSCF/VCI

Computational protocols for the simulation of two-dimensional infrared (2D IR) spectroscopy usually rely on vibrational exciton models which require an empirical parameterization. Here, we present an efficient quantum-chemical protocol for predicting static 2D IR spectra that does not require any empirical parameters. For the calculation of anharmonic vibrational energy levels and transition dipole moments, we employ the localized-mode vibrational self-consistent field (L-VSCF)/vibrational configuration interaction (L-VCI) approach previously established for (linear) anharmonic theoretical vibrational spectroscopy [P. T. Panek and C. R. Jacob, ChemPhysChem 15, 3365-3377 (2014)]. We demonstrate that with an efficient expansion of the potential energy surface using anharmonic one-mode potentials and harmonic two-mode potentials, 2D IR spectra of metal carbonyl complexes and dipeptides can be predicted reliably. We further show how the close connection between L-VCI and vibrational exciton models can be exploited to extract the parameters of such models from those calculations. This provides a novel route to the fully quantum-chemical parameterization of vibrational exciton models for predicting 2D IR spectra.

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.

Multidimensional insights into the repeated electromagnetic field stimulation and biosystems interaction in aging and age-related diseases

We provide a multidimensional sequence of events that describe the electromagnetic field (EMF) stimulation and biological system interaction. We describe this process from the quantum to the molecular, cellular, and organismal levels. We hypothesized that the sequence of events of these interactions starts with the oscillatory effect of the repeated electromagnetic stimulation (REMFS). These oscillations affect the interfacial water of an RNA causing changes at the quantum and molecular levels that release protons by quantum tunneling. Then protonation of RNA produces conformational changes that allow it to bind and activate Heat Shock Transcription Factor 1 (HSF1). Activated HSF1 binds to the DNA expressing chaperones that help regulate autophagy and degradation of abnormal proteins. This action helps to prevent and treat diseases such as Alzheimer's and Parkinson's disease (PD) by increasing clearance of pathologic proteins. This framework is based on multiple mathematical models, computer simulations, biophysical experiments, and cellular and animal studies. Results of the literature review and our research point towards the capacity of REMFS to manipulate various networks altered in aging (Reale et al. PloS one 9, e104973, 2014), including delay of cellular senescence (Perez et al. 2008, Exp Gerontol 43, 307-316) and reduction in levels of amyloid-β peptides (Aβ) (Perez et al. 2021, Sci Rep 11, 621). Results of these experiments using REMFS at low frequencies can be applied to the treatment of patients with age-related diseases. The use of EMF as a non-invasive therapeutic modality for Alzheimer's disease, specifically, holds promise. It is also necessary to consider the complicated and interconnected genetic and epigenetic effects of the REMFS-biological system's interaction while avoiding any possible adverse effects.

Chiral Sum Frequency Generation Spectroscopy Detects Double-Helix DNA at Interfaces

Many DNA-based technologies involve the immobilization of DNA and therefore require a fundamental understanding of the DNA structure-function relationship at interfaces. We present three immobilization methods compatible with chiral sum frequency generation (SFG) spectroscopy at interfaces. They are the "anchor" method for covalently attaching DNA on a glass surface, the "island" method for dropcasting DNA on solid substrates, and the "buoy" method using a hydrocarbon moiety for localizing DNA at the air-water interface. Although SFG was previously used to probe DNA, the chiral and achiral SFG responses of single-stranded and double-stranded DNA have not been compared systemically. Using the three immobilization methods, we obtain the achiral and chiral C-H stretching spectra. The results introduce four potential applications of chiral SFG. First, chiral SFG gives null response from single-stranded DNA but prominent signals from double-stranded DNA, providing a simple binary readout for label-free detection of DNA hybridization. Second, with heterodyne detection, chiral SFG gives an opposite-signed spectral response useful for distinguishing native (D-) right-handed double helix from non-native (L-) left-handed double helix. Third, chiral SFG captures the aromatic C-H stretching modes of nucleobases that emerge upon hybridization, revealing the power of chiral SFG to probe highly localized molecular structures within DNA. Finally, chiral SFG is sensitive to macroscopic chirality but not local chiral centers and thus can detect not only canonical antiparallel double helix but also other DNA secondary structures, such as a poly-adenine parallel double helix. Our work benchmarks the SFG responses of DNA immobilized by the three distinct methods, building a basis for new chiral SFG applications to solve fundamental and biotechnological problems.

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.

Structural Insights Into Tautomeric Dynamics in Nucleic Acids and in Antiviral Nucleoside Analogs

DNA (2'-deoxyribonucleic acid) and RNA (ribonucleic acid) play diverse functional roles in biology and disease. Despite being comprised primarily of only four cognate nucleobases, nucleic acids can adopt complex three-dimensional structures, and RNA in particular, can catalyze biochemical reactions to regulate a wide variety of biological processes. Such chemical versatility is due in part to the phenomenon of nucleobase tautomerism, whereby the bases can adopt multiple, yet distinct isomeric forms, known as tautomers. For nucleobases, tautomers refer to structural isomers that differ from one another by the position of protons. By altering the position of protons on nucleobases, many of which play critical roles for hydrogen bonding and base pairing interactions, tautomerism has profound effects on the biochemical processes involving nucleic acids. For example, the transient formation of minor tautomers during replication could generate spontaneous mutations. These mutations could arise from the stabilization of mismatches, in the active site of polymerases, in conformations involving minor tautomers that are indistinguishable from canonical base pairs. In this review, we discuss the evidence for tautomerism in DNA, and its consequences to the fidelity of DNA replication. Also reviewed are RNA systems, such as the riboswitches and self-cleaving ribozymes, in which tautomerism plays a functional role in ligand recognition and catalysis, respectively. We also discuss tautomeric nucleoside analogs that are efficacious as antiviral drug candidates such as molnupiravir for coronaviruses and KP1212 for HIV. The antiviral efficacy of these analogs is due, in part, to their ability to exist in multiple tautomeric forms and induce mutations in the replicating viral genomes. From a technical standpoint, minor tautomers of nucleobases are challenging to identify directly because they are rare and interconvert on a fast, millisecond to nanosecond, time scale. Nevertheless, many approaches including biochemical, structural, computational and spectroscopic methods have been developed to study tautomeric dynamics in RNA and DNA systems, and in antiviral nucleoside analogs. An overview of these methods and their applications is included here.

Advances, challenges and perspectives of quantum chemical approaches in molecular spectroscopy of the condensed phase

The purpose of this review is to demonstrate advances, challenges and perspectives of quantum chemical approaches in molecular spectroscopy of the condensed phase. Molecular spectroscopy, particularly vibrational spectroscopy and electronic spectroscopy, has been used extensively for a wide range of areas of chemical sciences and materials science as well as nano- and biosciences because it provides valuable information about structure, functions, and reactions of molecules. In the meantime, quantum chemical approaches play crucial roles in the spectral analysis. They also yield important knowledge about molecular and electronic structures as well as electronic transitions. The combination of spectroscopic approaches and quantum chemical calculations is a powerful tool for science, in general. Thus, our article, which treats various spectroscopy and quantum chemical approaches, should have strong implications in the wider scientific community. This review covers a wide area of molecular spectroscopy from far-ultraviolet (FUV, 120-200 nm) to far-infrared (FIR, 400-10 cm-1)/terahertz and Raman spectroscopy. As quantum chemical approaches, we introduce several anharmonic approaches such as vibrational self-consistent field (VSCF) and the combination of periodic harmonic calculations with anharmonic corrections based on finite models, grid-based techniques like the Numerov approach, the Cartesian coordinate tensor transfer (CCT) method, Symmetry-Adapted Cluster Configuration-Interaction (SAC-CI), and the ZINDO (Semi-empirical calculations at Zerner's Intermediate Neglect of Differential Overlap). One can use anharmonic approaches and grid-based approaches for both infrared (IR) and near-infrared (NIR) spectroscopy, while CCT methods are employed for Raman, Raman optical activity (ROA), FIR/terahertz and low-frequency Raman spectroscopy. Therefore, this review overviews cross relations between molecular spectroscopy and quantum chemical approaches, and provides various kinds of close-reality advanced spectral simulation for condensed phases.

Combining steady state and temperature jump IR spectroscopy to investigate the allosteric effects of ligand binding to dsDNA

Changes in the structural dynamics of double stranded (ds)DNA upon ligand binding have been linked to the mechanism of allostery without conformational change, but direct experimental evidence remains elusive. To address this, a combination of steady state infrared (IR) absorption spectroscopy and ultrafast temperature jump IR absorption measurements has been used to quantify the extent of fast (∼100 ns) fluctuations in (ds)DNA·Hoechst 33258 complexes at a range of temperatures. Exploiting the direct link between vibrational band intensities and base stacking shows that the absolute magnitude of the change in absorbance caused by fast structural fluctuations following the temperature jump is only weakly dependent on the starting temperature of the sample. The observed fast dynamics are some two orders of magnitude faster than strand separation and associated with all points along the 10-base pair duplex d(GCATATATCC). Binding the Hoechst 33258 ligand causes a small but consistent reduction in the extent of these fast fluctuations of base pairs located outside of the ligand binding region. These observations point to a ligand-induced reduction in the flexibility of the dsDNA near the binding site, consistent with an estimated allosteric propagation length of 15 Å, about 5 base pairs, which agrees well with both molecular simulation and coarse-grained statistical mechanics models of allostery leading to cooperative ligand binding.

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.

Advanced DNA Detection via Multispectral Plasmonic Metasurfaces

We propose and demonstrate a sensing platform based on plasmonic metasurfaces for the detection of very low concentrations of deoxyribo-nucleic acid (DNA) fragments. The platform relies on surface-enhanced infrared absorption spectroscopy, implemented via a multispectral metasurface. Specifically, different regions (“pixels”) are engineered so as to separately cover the medium-infrared range of the electromagnetic spectrum extending from the functional-groups to the fingerprint region of a single analyte. In conjunction with a suitable bio-functionalization, this enables univocal and label-free recognition of specific molecules. For experimental validation, we fabricate a large-area gold metasurface on a silicon chip, and functionalize it with a recognition layer of peptide nucleic acid (PNA). Our experimental results indicate the possibility to detect complementary DNA fragments in concentrations as low as 50 fM, i.e., well below the value attained by standard methods, with additional advantages in terms of processing time, versatility and ease of implementation/operation.

Computing the frequency fluctuation dynamics of highly coupled vibrational transitions using neural networks

The description of frequency fluctuations for highly coupled vibrational transitions has been a challenging problem in physical chemistry. In particular, the complexity of their vibrational Hamiltonian does not allow us to directly derive the time evolution of vibrational frequencies for these systems. In this paper, we present a new approach to this problem by exploiting the artificial neural network to describe the vibrational frequencies without relying on the deconstruction of the vibrational Hamiltonian. To this end, we first explored the use of the methodology to predict the frequency fluctuations of the amide I mode of N-methylacetamide in water. The results show good performance compared with the previous experimental and theoretical results. In the second part, the neural network approach is used to investigate the frequency fluctuations of the highly coupled carbonyl stretch modes for the organic carbonates in the solvation shell of the lithium ion. In this case, the frequency fluctuation predicted by the neural networks shows a good agreement with the experimental results, which suggests that this model can be used to describe the dynamics of the frequency in highly coupled transitions.

Large-scale investigation of the effects of nucleobase sequence on fluorescence excitation and Stokes shifts of DNA-stabilized silver clusters

DNA-stabilized silver clusters (AgN-DNAs) exhibit diverse sequence-programmed fluorescence, making these tunable nanoclusters promising sensors and bioimaging probes. Recent advances in the understanding of AgN-DNA structures and optical properties have largely relied on detailed characterization of single species isolated by chromatography. Because most AgN-DNAs are unstable under chromatography, such studies do not fully capture the diversity of these clusters. As an alternative method, we use high-throughput synthesis and spectroscopy to measure steady state Stokes shifts of hundreds of AgN-DNAs. Steady state Stokes shift is of interest because its magnitude is determined by energy relaxation processes which may be sensitive to specific cluster geometry, attachment to the DNA template, and structural engagement of solvent molecules. We identify 305 AgN-DNA samples with single-peaked emission and excitation spectra, a characteristic of pure solutions and single emitters, which thus likely contain a dominant emissive AgN-DNA species. Steady state Stokes shifts of these samples vary widely, are in agreement with values reported for purified clusters, and are several times larger than for typical organic dyes. We then examine how DNA sequence selects AgN-DNA excitation energies and Stokes shifts, comment on possible mechanisms for energy relaxation processes in AgN-DNAs, and discuss how differences in AgN-DNA structure and DNA conformation may result in the wide distribution of optical properties observed here. These results may aid computational studies seeking to understand the fluorescence process in AgN-DNAs and the relations of this process to AgN-DNA structure.

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.

Time-Resolved Vibrational Fingerprints for Two Silver Cluster-DNA Fluorophores

DNA-templated silver clusters are chromophores in which the nucleobases encode the cluster spectra and brightness. We describe the coordination environments of two nearly identical Ag₁₀⁶⁺ clusters that form with 18-nucleotide strands CCCCA CCCCT CCCX TTTT, with X = guanosine and inosine. For the first time, femtosecond time-resolved infrared (TRIR) spectroscopy with visible excitation and mid-infrared probing is used to correlate the response of nucleobase vibrational modes to electronic excitation of the metal cluster. A rich pattern of transient TRIR peaks in the 1400-1720 cm^-1 range decays synchronously with the visible emission. Specific infrared signatures associated with the single guanosine/inosine along with a subset of cytidines, but not the thymidines, are observed. These fingerprints suggest that the network of bonds between a silver cluster adduct and its polydentate DNA ligands can be deciphered to rationally tune the coordination and thus spectra of molecular silver chromophores.

DNA minor-groove binder Hoechst 33258 destabilizes base-pairing adjacent to its binding site

Understanding the dynamic interactions of ligands to DNA is important in DNA-based nanotechnologies. By structurally tracking the dissociation of Hoechst 33258-bound DNA (d(CGCAAATTTGCG)₂) complex (H-DNA) with T-jump 2D-IR spectroscopy, the ligand is found to strongly disturb the stability of the three C:G base pairs adjacent to A:T the binding site, with the broken base pairs being more than triple at 100 ns. The strong stabilization effect of the ligand on DNA duplex makes this observation quite striking, which dramatically increases the melting temperature and dissociation time. MD simulations demonstrate an important role of hydration water and counter cations in maintaining the separation of terminal base pairs. The hydrogen bonds between the ligand and thymine carbonyls are crucial in stabilizing H-DNA, whose breaking signal appearing prior to the complete dissociation. Thermodynamic analysis informs us that H-DNA association is a concerted process, where H cooperates with DNA single strands in forming H-DNA.

Temperature-Jump 2D IR Spectroscopy with Intensity-Modulated CW Optical Heating

Pulsed temperature-jump (T-jump) spectroscopy with infrared (IR) detection has been widely used to study biophysical processes occurring from nanoseconds to ∼1 ms with structural sensitivity. However, many systems exhibit structural dynamics on time scales longer than the millisecond barrier that is set by the time scale for thermal relaxation of the sample. We developed a linear and nonlinear infrared spectrometer coupled to an intensity-modulated continuous wave (CW) laser to probe T-jump-initiated chemical reactions from <1 ms to seconds. Time-dependent modulation of the CW laser leads to a <1 ms heating time as well as a constant final temperature (±3%) for the duration of the heating time. Temperature changes of up to 75 °C in D₂O are demonstrated, allowing for nonequilibrium measurements inaccessible to standard pulsed optical T-jump setups. T-jump linear absorption, pump-probe, and two-dimensional IR (2D IR) spectroscopy are applied to the unfolding and refolding of ubiquitin and a model intercalated motif (i-motif) DNA sequence, and analysis of the observed signals is used to demonstrate the limits and utility of each method. Overall, the ability to probe temperature-induced chemical processes from <1 ms to many seconds with 2D IR spectroscopy provides multiple new avenues for time-dependent spectroscopy in chemistry and biophysics.

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.

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.

Biomolecular and bioanalytical applications of infrared spectroscopy - A review

Infrared (IR; or mid-infrared, MIR; 4000-400 cm^-1; 2500-25,000 nm) spectroscopy has become one of the most powerful and versatile tools at the disposal of modern bioscience. Because of its high molecular specificity, applicability to wide variety of samples, rapid measurement and non-invasivity, IR spectroscopy forms a potent approach to elucidate qualitative and quantitative information from various kinds of biological material. For these reasons, it became an established bioanalytical technique with diverse applications. This work aims to be a comprehensive and critical review of the recent accomplishments in the field of biomolecular and bioanalytical IR spectroscopy. That progress is presented on a wider background, with fundamental characteristics, the basic principles of the technique outlined, and its scientific capability directly compared with other methods being used in similar fields (e.g. near-infrared, Raman, fluorescence). The article aims to present a complete examination of the topic, as it touches the background phenomena, instrumentation, spectra processing and data analytical methods, spectra interpretation and related information. To suit this goal, the article includes a tutorial information essential to obtain a thorough perspective of bio-related applications of the reviewed methodologies. The importance of the fundamental factors to the final performance and applicability of IR spectroscopy in various areas of bioscience is explained. This information is interpreted in critical way, with aim to gain deep understanding why IR spectroscopy finds extraordinarily intensive use in this remarkably diverse and dynamic field of research and utility. The major focus is placed on the diversity of the applications in which IR biospectroscopy has been established so far and those onto which it is expanding nowadays. This includes qualitative and quantitative analytical spectroscopy, spectral imaging, medical diagnosis, monitoring of biophysical processes, and studies of physicochemical properties and dynamics of biomolecules. The application potential of IR spectroscopy in light of the current accomplishments and the future prospects is critically evaluated and its significance in the progress of bioscience is comprehensively presented.

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.

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.

2D-IR studies of cyanamides (NCN) as spectroscopic reporters of dynamics in biomolecules: Uncovering the origin of mysterious peaks

Cyanamides (NCN) have been shown to have a larger transition dipole strength than cyano-probes. In addition, they have similar structural characteristics and vibrational lifetimes to the azido-group, suggesting their utility as infrared (IR) spectroscopic reporters for structural dynamics in biomolecules. To access the efficacy of NCN as an IR probe to capture the changes in the local environment, several model systems were evaluated via 2D IR spectroscopy. Previous work by Cho [G. Lee, D. Kossowska, J. Lim, S. Kim, H. Han, K. Kwak, and M. Cho, J. Phys. Chem. B 122(14), 4035-4044 (2018)] showed that phenylalanine analogues containing NCN show strong anharmonic coupling that can complicate the interpretation of structural dynamics. However, when NCN is embedded in 5-membered ring scaffolds, as in N-cyanomaleimide and N-cyanosuccinimide, a unique band structure is observed in the 2D IR spectrum that is not predicted by simple anharmonic frequency calculations. Further investigation indicated that electron delocalization plays a role in the origins of the band structure. In particular, the origin of the lower frequency transitions is likely a result of direct interaction with the solvent.

Oxidized Derivatives of 5-Methylcytosine Alter the Stability and Dehybridization Dynamics of Duplex DNA

The naturally occurring nucleobase 5-methylcytosine (mC) and its oxidized derivatives 5-hydroxymethylcytosine (hmC), 5-formylcytosine (fC), and 5-carboxylcytosine (caC) play important roles in epigenetic regulation and, along with cytosine (C), represent nucleobases currently implicated in the active cytosine demethylation pathway. Despite considerable interest in these modified bases, their impact on the thermodynamic stability of double-stranded DNA (dsDNA) remains ambiguous and their influence on hybridization kinetics and dynamics is even less well-understood. To address these unknowns, we employ steady-state and time-resolved infrared spectroscopy to measure the influence of cytosine modification on the thermodynamics and kinetics of hybridization by assessing the impact on local base pairing dynamics, shifts in the stability of the duplex state, and changes to the hybridization transition state. Modification with mC leads to more tightly bound base pairing below the melting transition and stabilizes the duplex relative to canonical DNA, but the free energy barrier to dehybridization at physiological temperature is nevertheless reduced slightly. Both hmC and fC lead to an increase in local base pair fluctuations, a reduction in the cooperativity of duplex melting, and a lowering of the dissociation barrier, but these effects are most pronounced when the 5-position is formylated. The caC nucleobase demonstrates little impact on dsDNA under neutral conditions, but we find that this modification can dynamically switch between C-like and fC-like behavior depending on the protonation state of the 5-position carboxyl group. Our results provide a consistent thermodynamic and kinetic framework with which to describe the modulation of the physical properties of double-stranded DNA containing these modified nucleobases.

5-Carboxylcytosine and Cytosine Protonation Distinctly Alter the Stability and Dehybridization Dynamics of the DNA Duplex

Applications associated with nucleobase protonation events are grounded in their fundamental impact on DNA thermodynamics, structure, and hybridization dynamics. Of the canonical nucleobases, N3 protonation of cytosine (C) is the most widely utilized in both biology and nanotechnology. Naturally occurring C derivatives that shift the N3 pK_a introduce an additional level of tunability. The epigenetic nucleobase 5-carboxylcytosine (caC) presents a particularly interesting example since this derivative forms Watson-Crick base pairs of similar stability and displays pH-dependent behavior over the same range as the canonical nucleobase. However, the titratable group in caC corresponds to the exocyclic carboxyl group rather than N3, and the implications of these divergent protonation events toward DNA hybridization thermodynamics, kinetics, and base pairing dynamics remain poorly understood. Here, we study the pH dependence of these physical properties using model oligonucleotides containing C and caC with FTIR and temperature-jump IR spectroscopy. We demonstrate that N3 protonation of C completely disrupts duplex stability, leading to large shifts in the duplex/single-strand equilibrium, a reduction in the cooperativity of melting, and an acceleration in the rate of duplex dissociation. In contrast, while increasing 5-carboxyl protonation in caC-containing duplexes induces an increase in base pair fluctuations, the DNA duplex can tolerate substantial protonation without significant perturbation to the duplex/single-strand equilibrium. However, 5-carboxyl protonation has a large impact on hybridization kinetics by reducing the transition state free energy. Our thermodynamic and kinetic analysis provides new insight on the impact of two divergent protonation mechanisms in naturally occurring nucleobases on the biophysical properties of DNA.

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.