Observation of coherent delocalized phonon-like modes in DNA under physiological conditions

Low-Frequency Spectra of Hydrated Ionic Liquids with Kosmotropic and Chaotropic Anions

In this study, we investigated the water concentration dependence of the intermolecular vibrations of two hydrated ionic liquids (ILs), cholinium dihydrogen phosphate ([ch][dhp]) and cholinium bromide ([ch]Br), using femtosecond Raman-induced Kerr effect spectroscopy (fs-RIKES). The anions of the former and latter hydrated ILs are kosmotropic and chaotropic, respectively. We found that the spectral peak of ∼50 cm-1 shifted to the low-frequency side in hydrated [ch][dhp], indicating the weakening of its intermolecular interactions. In contrast, no change in the peak frequency of the low-frequency band at ∼50 cm-1 was observed with increasing water concentration in hydrated [ch]Br. The vibrational density of states (VDOS) spectra generated from molecular dynamics (MD) simulations were in qualitative agreement with the experimental results. Decomposition analysis of the VDOS spectra for each component revealed that the red shift of the low-frequency band in the hydrated [ch][dhp] upon water addition was essentially due to the contributions of anions and water rather than that of the cholinium cation. We also found from the low-frequency spectra of the two hydrated ILs that they differed in the concentration dependence of the 180 cm-1 band, which is assigned as a hindered translational motion of water molecules combined to form O···O stretching motions. From the relationship between the peak frequency of the low-frequency band and the bulk parameter, which is the square root of the surface tension divided by the density, we found that the peak frequency in the hydrated IL with kosmotropic [dhp]- depends on the bulk parameter, similar to the case for an aqueous solution of the typical deep eutectic solvent reline. However, the peak frequency of the hydrated IL with chaotropic Br- is constant with the bulk parameter.

Vibration spectra of DNA and RNA segments

The dispersion curves and density of states are used to analyze the vibrational characteristics of DNA and RNA segments. This is done using a harmonic Hamiltonian and the Green's function technique. Two configurations of DNA and RNA, finite and cyclic, have been investigated and compared to their infinite counterparts. For the DNA molecule, three models, including a fishbone model, a ldder model, and a fishbone ladder model, have been employed, while the RNA molecule has been represented using a half fishbone model. To enhance the realism of DNA and RNA simulations, the unit cells within each infinite system as well as the length of the finite and cyclic cases are gradually enlarged. The connections between the sub-sites have been modeled using linear springs, where the stiffness of the vertical springs exhibits random variations throughout the length of the DNA and RNA models. Shorter DNA and RNA segments exhibit additional peaks in their density of states, resulting in more bands in dispersion curves. This indicates that as the number of building blocks grows in these segments, their curves resemble those of infinite systems. These findings have practical implications for studying the vibration characteristics of similar macro-systems.

Harmonic and anharmonic studies on THz spectra of two vanillin polymorphs

Polymorphism commonly exists in organic molecular crystals. The fingerprint features in low-frequency vibrational range are important information reflecting different intermolecular interactions of polymorphs. Interpreting these features is very helpful to understand vibrational property of polymorphs and reveal the thermodynamic stability. In this work, the low-frequency vibrations of form I and II of vanillin are investigated using terahertz time-domain spectroscopy. Static DFT calculation and ab initio molecular dynamics (AIMD) are employed to interpret their low-frequency vibrations of both forms in harmonic and anharmonic ways, respectively. Their low-frequency vibration characteristics in harmonic calculations are discussed, and anharmonic mode couplings between OH bond stretch and the stretching and bending motion of hydrogen bonds are uncovered. Moreover, the thermodynamic energies including electronic potential energy and vibrational/kinetic energy arising from nuclear motions are calculated. The result reveals that the stability order of the two forms is mainly dependent on their electric potential energy difference.

Lifting Hofmeister's Curse: Impact of Cations on Diffusion, Hydrogen Bonding, and Clustering of Water

Water plays a role in the stability, reactivity, and dynamics of the solutes that it contains. The presence of ions alters this capacity by changing the dynamics and structure of water. However, our understanding of how and to what extent this occurs is still incomplete. Here, a study of the low-frequency Raman spectra of aqueous solutions of various cations by using optical Kerr-effect spectroscopy is presented. This technique allows for the measurement of the changes that ions cause in both the diffusive dynamics and the vibrations of the hydrogen-bond structure of water. It is found that when salts are added, some of the water molecules become part of the ion solvation layers, while the rest retain the same diffusional properties as those of pure water. The slowing of the dynamics of the water molecules in the solvation shell of each ion was found to depend on its charge density at infinite dilution conditions and on its position in the Hofmeister series at higher concentrations. It is also observed that all cations weaken the hydrogen-bond structure of the solution and that this weakening depends only on the size of the cation. Finally, evidence is found that ions tend to form amorphous aggregates, even at very dilute concentrations. This work provides a novel approach to water dynamics that can be used to better study the mechanisms of solute nucleation and crystallization, the structural stability of biomolecules, and the dynamic properties of complex solutions, such as water-in-salt electrolytes.

Examining DNA breathing with pyDNA-EPBD

MOTIVATION

The two strands of the DNA double helix locally and spontaneously separate and recombine in living cells due to the inherent thermal DNA motion. This dynamics results in transient openings in the double helix and is referred to as "DNA breathing" or "DNA bubbles." The propensity to form local transient openings is important in a wide range of biological processes, such as transcription, replication, and transcription factors binding. However, the modeling and computer simulation of these phenomena, have remained a challenge due to the complex interplay of numerous factors, such as, temperature, salt content, DNA sequence, hydrogen bonding, base stacking, and others.

RESULTS

We present pyDNA-EPBD, a parallel software implementation of the Extended Peyrard-Bishop-Dauxois (EPBD) nonlinear DNA model that allows us to describe some features of DNA dynamics in detail. The pyDNA-EPBD generates genomic scale profiles of average base-pair openings, base flipping probability, DNA bubble probability, and calculations of the characteristically dynamic length indicating the number of base pairs statistically significantly affected by a single point mutation using the Markov Chain Monte Carlo algorithm.

AVAILABILITY AND IMPLEMENTATION

pyDNA-EPBD is supported across most operating systems and is freely available at https://github.com/lanl/pyDNA_EPBD. Extensive documentation can be found at https://lanl.github.io/pyDNA_EPBD/.

Examining DNA Breathing with pyDNA-EPBD

Motivation

The two strands of the DNA double helix locally and spontaneously separate and recombine in living cells due to the inherent thermal DNA motion.This dynamics results in transient openings in the double helix and is referred to as "DNA breathing" or "DNA bubbles." The propensity to form local transient openings is important in a wide range of biological processes, such as transcription, replication, and transcription factors binding. However, the modeling and computer simulation of these phenomena, have remained a challenge due to the complex interplay of numerous factors, such as, temperature, salt content, DNA sequence, hydrogen bonding, base stacking, and others.

Results

We present pyDNA-EPBD, a parallel software implementation of the Extended Peyrard-Bishop- Dauxois (EPBD) nonlinear DNA model that allows us to describe some features of DNA dynamics in detail. The pyDNA-EPBD generates genomic scale profiles of average base-pair openings, base flipping probability, DNA bubble probability, and calculations of the characteristically dynamic length indicating the number of base pairs statistically significantly affected by a single point mutation using the Markov Chain Monte Carlo (MCMC) algorithm.

Tuning symmetry-protected quasi bound state in the continuum using terahertz meta-atoms of rotational and reflectional symmetry

Conventionally, a symmetry-protected quasi bound state of the continuum (BIC) becomes achievable by breaking the C₂ symmetry of meta-atoms. Our work exhibits a novel approach to achieving dual band quasi-BIC by breaking the C_2v symmetry into C_s symmetry. Also, we show that a single band quasi-BIC can be achieved by breaking the C_2v symmetry into C₂ symmetry. Our metasurface of C_2v symmetry is composed of double gaps split ring resonator (DSRR), and it degrades to C₂ symmetry when the double gaps are displaced in opposite directions. One band quasi-BIC can be observed occurring at around 0.36 and 0.61 THz respectively with the metasurface excited by x- and y-polarized terahertz radiation, respectively. A couple of dark dipole oscillator dominates the quasi-BIC at 0.36 THz, while a quadruple-like oscillator dominates the quasi-BIC at 0.61 THz. The damping ratio and coupling coefficients of the above single quasi-BIC are close to the orthogonal polarization of the incident terahertz wave. However, the metasurface of the DSRR array degrades down to C_s symmetry when the double gaps are displaced in the same directions. A dual band quasi-BIC (0.23 THz and 0.62 THz) is found to be sensitive to the y-polarized terahertz radiation. It is found that the inductive-capacitive (LC) resonance results in quasi-BIC at 0.23 THz, while a quadrupole-like oscillation results in quasi-BIC at 0.62 THz. The quasi-BIC at 0.62 THz has a higher coupling coefficient and lower damping ratio than quasi-BIC at 0.23 THz in a metasurface of C_s symmetry. The realization of the above locally symmetric breaking on the quasi-BIC of terahertz metasurfaces is helpful for the innovation of multi-band terahertz biosensors.

Intermolecular Dynamics and Structure in Aqueous Lidocaine Hydrochloride Solutions

We investigated the intermolecular dynamics and static structure in the aqueous solutions of lidocaine hydrochloride (LDHCl) in the concentration range of [LDHCl] = 0-2.00 M using femtosecond Raman-induced Kerr effect spectroscopy (fs-RIKES), small- and wide-angle X-ray scattering (SWAXS), and dynamic light scattering (DLS). For the fs-RIKES experiments, the concentration dependence of the difference low-frequency spectra of the aqueous LDHCl solutions relative to the neat water, which was mainly due to the intermolecular vibrations, was characterized using an exponential function with a characteristic concentration of ∼1 M. For the SWAXS experiments, we observed a manifestation of an excess scattering component centered within a range of 8-10 nm^-1 in the aqueous LDHCl solutions. The results of Fourier inversion and further deconvolution analyses unambiguously demonstrated that lidocaines assemble into a nanometer-sized micelle-like structure with the innermost core (∼0.3 nm) and outer shell (∼0.5 nm), respectively. The DLS experiments also found nanometer-sized aggregates and further indicated evidence of the clusters of the aggregates. The results of viscosities, densities, and surface tensions of the solutions and the quantum chemistry calculations supported the unique features of the microscopic intermolecular interaction and the micelle-like aggregation.

Long-range DNA-water interactions

DNA functions only in aqueous environments and adopts different conformations depending on the hydration level. The dynamics of hydration water and hydrated DNA leads to rotating and oscillating dipoles that, in turn, give rise to a strong megahertz to terahertz absorption. Investigating the impact of hydration on DNA dynamics and the spectral features of water molecules influenced by DNA, however, is extremely challenging because of the strong absorption of water in the megahertz to terahertz frequency range. In response, we have employed a high-precision megahertz to terahertz dielectric spectrometer, assisted by molecular dynamics simulations, to investigate the dynamics of water molecules within the hydration shells of DNA as well as the collective vibrational motions of hydrated DNA, which are vital to DNA conformation and functionality. Our results reveal that the dynamics of water molecules in a DNA solution is heterogeneous, exhibiting a hierarchy of four distinct relaxation times ranging from ∼8 ps to 1 ns, and the hydration structure of a DNA chain can extend to as far as ∼18 Å from its surface. The low-frequency collective vibrational modes of hydrated DNA have been identified and found to be sensitive to environmental conditions including temperature and hydration level. The results reveal critical information on hydrated DNA dynamics and DNA-water interfaces, which impact the biochemical functions and reactivity of DNA.

Low-Frequency Spectra of 1-Methyl-3-octylimidazolium Tetrafluoroborate Mixtures with Poly(ethylene glycol) by Femtosecond Raman-Induced Kerr Effect Spectroscopy

This is the first report on low-frequency spectra of ionic liquid (IL)/polymer mixtures using femtosecond Raman-induced Kerr effect spectroscopy. We studied mixtures of 1-methyl-3-octylimidazolium tetrafluoroborate ([MOIm][BF₄]) and poly(ethylene glycol) (PEG) with M_n = 400 (PEG400) at various concentrations. To elucidate the unique features of the IL/polymer mixture system, mixtures of PEG400 with a molecular liquid, 1-octhylimidazole (OIm), which is a neutral analog of the cation, were also studied. In addition, mixtures of [MOIm][BF₄] with ethylene glycol (EG) and poly(ethylene glycol) with M_n = 4000 (PEG4000) were also investigated. The first moments of broad low-frequency spectra, mainly due to intermolecular vibrations for the [MOIm][BF₄]/PEG400 and OIm/PEG400, increased slightly with increasing concentration of PEG400, indicating that microscopic intermolecular interactions, in general, are slightly enhanced. We also compared the [MOIm][BF₄] mixtures with EG, PEG400, and PEG4000 at concentrations of 5 and 10 wt % PEG or EG. The low-frequency spectra of samples with the same concentrations were quite similar, but a comparison of the normalized spectra showed that the spectral intensity in the low-frequency region below ∼50 cm^-1 of the [MOIm][BF₄] mixtures with PEG400 and PEG4000 is somewhat lower than that of the [MOIm][BF₄] mixtures with EG. Although the effect of the polymer is small compared to other polymer solution systems, this feature is attributed to a suppression of translational motion in the mixtures of [MOIm][BF₄] with PEG compared to the mixtures of [MOIm][BF₄] with EG due to the greater mass of PEG than EG. Density, surface tension, viscosity, and electrical conductivity were also estimated. From Walden plots, it was found that the [MOIm][BF₄]/PEG4000 system showed more ideal electrical conductive behavior than the [MOIm][BF₄]/PEG400 and [MOIm][BF₄]/EG systems.

Protein-Ligand Binding Molecular Details Revealed by Terahertz Optical Kerr Spectroscopy: A Simulation Study

Picosecond fast motions and their involvement in the biochemical processes such as protein-ligand binding has engaged significant attention. Terahertz optical Kerr spectroscopy (OKE) has the superior potential to probe these fast motions directly. Application of OKE in protein-ligand binding study is, however, limited by the difficulty of quantitative atomistic interpretation, and the calculation of Kerr spectrum for entire solvated protein complex was considered not yet feasible, due to the lack of one consistent polarizable model for both configuration sampling and polarizability calculation. Here, we analyzed the biochemical relevance of OKE to the lysozyme-triacetylchitotriose binding based on the first OKE simulation using one consistent Drude polarizable model. An analytical multipole and induced dipole scheme was employed to calculate the off-diagonal Drude polarizability more efficiently and accurately. Further theoretical analysis revealed how the subtle twisting and stiffening of aromatic protein residues' spatial arrangement as well as the confinement of small water clusters between ligand and protein cavity due to the ligand binding can be examined using Kerr spectroscopy. Comparison between the signals of bound complex and that of uncorrelated protein/ligand demonstrated that binding action alone has reflection in the OKE spectrum. Our study indicated OKE as a powerful terahertz probe for protein-ligand binding chemistry and dynamics.

Bubble lifetimes in DNA gene promoters and their mutations affecting transcription

Relative lifetimes of inherent double stranded DNA openings with lengths up to ten base pairs are presented for different gene promoters and corresponding mutants that either increase or decrease transcriptional activity in the framework of the Peyrard-Bishop-Dauxois model. Extensive microcanonical simulations are used with energies corresponding to physiological temperature. The bubble lifetime profiles along the DNA sequences demonstrate a significant reduction of the average lifetime at the mutation sites when the mutated promoter decreases transcription, while a corresponding enhancement of the bubble lifetime is observed in the case of mutations leading to increased transcription. The relative difference in bubble lifetimes between the mutated and wild type promoters at the position of mutation varies from 20% to more than 30% as the bubble length decreases.

Low-frequency vibrational modes in G-quadruplexes reveal the mechanical properties of nucleic acids

Low-frequency vibrations play an essential role in biomolecular processes involving DNA such as gene expression, charge transfer, drug intercalation, and DNA–protein recognition. However, understanding the vibrational basis of these mechanisms relies on theoretical models due to the lack of experimental evidence. Here we present the low-frequency vibrational spectra of G-quadruplexes (structures formed by four strands of DNA) and B-DNA characterized using femtosecond optical Kerr-effect spectroscopy. Contrary to expectation, we found that G-quadruplexes show several strongly underdamped delocalized phonon-like modes that have the potential to contribute to the biology of the DNA at the atomic level. In addition, G-quadruplexes present modes at a higher frequency than B-DNA demonstrating that changes in the stiffness of the molecule alter its gigahertz to terahertz vibrational profile.

Low-frequency vibrations play an essential role in biomolecular processes involving DNA such as gene expression, charge transfer, drug intercalation, and DNA–protein recognition.

Ultrafast coherent motion and helix rearrangement of homodimeric hemoglobin visualized with femtosecond X-ray solution scattering

Ultrafast motion of molecules, particularly the coherent motion, has been intensively investigated as a key factor guiding the reaction pathways. Recently, X-ray free-electron lasers (XFELs) have been utilized to elucidate the ultrafast motion of molecules. However, the studies on proteins using XFELs have been typically limited to the crystalline phase, and proteins in solution have rarely been investigated. Here we applied femtosecond time-resolved X-ray solution scattering (fs-TRXSS) and a structure refinement method to visualize the ultrafast motion of a protein. We succeeded in revealing detailed ultrafast structural changes of homodimeric hemoglobin involving the coherent motion. In addition to the motion of the protein itself, the time-dependent change of electron density of the hydration shell was tracked. Besides, the analysis on the fs-TRXSS data of myoglobin allows for observing the effect of the oligomeric state on the ultrafast coherent motion.

Femtosecond time-resolved X-ray solution scattering (fs-TRXSS) measurements provide information on the structural dynamics of proteins in solution. Here, the authors present a structure refinement method for the analysis of fs-TRXSS data and use it to characterise the ultrafast structural changes of homodimeric haemoglobin.

Broadband terahertz signatures and vibrations of dopamine

Dopamine (DA) is an essential neurotransmitter and hormone of the nervous system, its structural and conformational properties play critical roles in biological functions and signal transmission processes. Although this neuroactive molecule has been studied extensively, the low-frequency vibration features that are closely related to the conformation and molecular interactions in the terahertz (THz) band still remain unclear. In this study, a broadband THz time-domain spectroscopy (THz-TDS) system in the frequency band of 0.5-18 THz was used to characterize the unique THz fingerprint of DA. In addition, density functional theory (DFT) calculations were performed to analyze the vibrational properties of DA. The results suggest that each THz resonant absorption peak of DA corresponds to specific vibrational modes, and the collective vibration also exists in the broadband THz range. Moreover, the interactions between the DA ligand and the D2 and D3 receptors were investigated by docking, and the simulated THz spectra were obtained. The results indicate the dominant role of hydrogen bonding interactions and the specificity of molecular conformation. This work may help to understand the resonance coupling between THz electromagnetic waves and neurotransmitters.

Nonthermal excitation effects mediated by sub-terahertz radiation on hydrogen exchange in ubiquitin

Water dynamics in the hydration layers of biomolecules play crucial roles in a wide range of biological functions. A hydrated protein contains multiple components of diffusional and vibrational dynamics of water and protein, which may be coupled at ∼0.1-THz frequency (10-ps timescale) at room temperature. However, the microscopic description of biomolecular functions based on various modes of protein-water-coupled motions remains elusive. A novel approach for perturbing the hydration dynamics in the subterahertz frequency range and probing them at the atomic level is therefore warranted. In this study, we investigated the effect of klystron-based, intense 0.1-THz excitation on the slow dynamics of ubiquitin using NMR-based measurements of hydrogen-deuterium exchange. We demonstrated that the subterahertz irradiation accelerated the hydrogen-deuterium exchange of the amides located in the interior of the protein and hydrophobic surfaces while decelerating this exchange in the amides located in the surface loop and short 3₁₀ helix regions. This subterahertz-radiation-induced effect was qualitatively contradictory to the increased-temperature-induced effect. Our results suggest that the heterogeneous water dynamics occurring at the protein-water interface include components that are nonthermally excited by the subterahertz radiation. Such subterahertz-excited components may be linked to the slow function-related dynamics of the protein.

Fluorescent bacterial biosensor E. coli/pTdcR-TurboYFP sensitive to terahertz radiation

A fluorescent biosensor E. coli/pTdcR-TurboYFP sensitive to terahertz (THz) radiation was developed via transformation of Escherichia coli (E. coli) cells with plasmid, in which the promotor of the tdcR gene controls the expression of yellow fluorescent protein TurboYFP. The biosensor was exposed to THz radiation in various vessels and nutrient media. The threshold and dynamics of fluorescence were found to depend on irradiation conditions. Heat shock or chemical stress yielded the absence of fluorescence induction. The biosensor is applicable to studying influence of THz radiation on the activity of tdcR promotor that is involved in the transport and metabolism of threonine and serine in E. coli.

Gre factors prevent thermal and mechanical stresses induced by terahertz irradiation during transcription

During transcription in cells, the transcription complex consisting of RNA polymerase, DNA and nascent RNA is exposed to fluctuating temperature and pressure. However, little is known about the mechanism of transcriptional homeostasis under fluctuating physical parameters. In this study, we generated these fluctuating parameters using pulsed local heating and acoustic waves in the reaction system of transcription by Escherichia coli RNA polymerase, using a terahertz free-electron laser. We demonstrated that transcription processes, including abortive initiation and elongation pausing, and the fidelity of elongation are significantly affected by the laser-based local perturbations. We also found that all these functional alternations in the transcription process are almost completely mitigated by the presence of Gre proteins. It is well known that Gre proteins enhance RNA cleavage of polymerase by binding to the pore structure termed secondary channel. Recently, the chaperone activities have also been proposed for Gre proteins, yet the details directly associated with transcription are largely unknown. Our finding indicates that Gre proteins are necessary for maintaining transcriptional homeostasis under thermal and mechanical stresses.

Base-Pairs' Correlated Oscillation Effects on the Charge Transfer in Double-Helix B-DNA Molecules

By introducing a suitable renormalization process, the charge carrier and phonon dynamics of a double-stranded helical DNA molecule are expressed in terms of an effective Hamiltonian describing a linear chain, where the renormalized transfer integrals explicitly depend on the relative orientations of the Watson-Crick base pairs, and the renormalized on-site energies are related to the electronic parameters of consecutive base pairs along the helix axis, as well as to the low-frequency phonons' dispersion relation. The existence of synchronized collective oscillations enhancing the π-π orbital overlapping among different base pairs is disclosed from the study of the obtained analytical dynamical equations. The role of these phonon-correlated, long-range oscillation effects on the charge transfer properties of double-stranded DNA homopolymers is discussed in terms of the resulting band structure.

Coherence preservation and electron-phonon interaction in electron transfer in DNA

We analyze the influence of electron-phonon (e-ph) interaction in a model for electron transfer (ET) processes in DNA in terms of the envelope function approach for spinless electrons. We are specifically concerned with the effect of e-ph interaction on the coherence of the ET process and how to model the interaction of DNA with phonon reservoirs of biological relevance. We assume that the electron bearing orbitals are half filled and derive the physics of e-ph coupling in the vicinity in reciprocal space. We find that at half filling, the acoustical modes are decoupled to ET at first order, while optical modes are predominant. The latter are associated with inter-strand vibrational modes in consistency with previous studies involving polaron models of ET. Coupling to acoustic modes depends on electron doping of DNA, while optical modes are always coupled within our model. Our results yield e-ph coupling consistent with estimates in the literature, and we conclude that large polarons are the main result of such e-ph interactions. This scenario will have strong consequences on decoherence of ET under physiological conditions due to relative isolation from thermal equilibration of the ET mechanism.

Origin of heat capacity increment in DNA folding: The hydration effect

BACKGROUND

Understanding DNA folding thermodynamics is crucial for prediction of DNA thermal stability. It is now well established that DNA folding is accompanied by a decrease of the heat capacity ∆c_{p, F}, however its molecular origin is not understood. In analogy to protein folding it has been assumed that this is due to dehydration of DNA constituents, however no evidence exists to support this conclusion.

METHODS

Here we analyze partial molar heat capacity of nucleic bases and nucleosides in aqueous solutions obtained from calorimetric experiments and calculate the hydration heat capacity contribution ∆c_p^hyd.

RESULTS

We present hydration heat capacity contributions of DNA constituents and show that they correlate with the solvent accessible surface area. The average contribution for nucleic base dehydration is +0.56 J mol^-1 K^-1 Å^-2 and can be used to estimate the ∆c_{p, F} contribution for DNA folding.

CONCLUSIONS

We show that dehydration is one of the major sources contributing to the observed ∆c_{p, F} increment in DNA folding. Other possible sources contributing to the overall ∆c_{p, F} should be significant but appear to compensate each other to high degree. The calculated ∆c_p^hyd for duplexes and noncanonical DNA structures agree excellently with the overall experimental ∆c_{p, F} values. By contrast, empirical parametrizations developed for proteins result in poor ∆c_p^hyd predictions and should not be applied to DNA folding.

GENERAL SIGNIFICANCE

Heat capacity is one of the main thermodynamic quantities that strongly affects thermal stability of macromolecules. At the molecular level the heat capacity in DNA folding stems from removal of water from nucleobases.

Low-Frequency (Gigahertz to Terahertz) Depolarized Raman Scattering Off n-Alkanes, Cycloalkanes, and Six-Membered Rings: A Physical Interpretation

Molecular liquids have long been known to undergo various distinct intermolecular motions, from fast librations and cage-rattling oscillations to slow orientational and translational diffusion. However, their resultant gigahertz to terahertz spectra are far from simple, appearing as broad shapeless bands that span many orders of magnitude of frequency, making meaningful interpretation troublesome. Ad hoc spectral line shape fitting has become a notoriously fine art in the field; a unified approach to handling such spectra is long overdue. Here we apply ultrafast optical Kerr-effect (OKE) spectroscopy to study the intermolecular dynamics of room-temperature n-alkanes, cycloalkanes, and six-carbon rings, as well as liquid methane and propane. This work provides stress tests and converges upon an experimentally robust model across simple molecular series and range of temperatures, providing a blueprint for the interpretation of the dynamics of van der Waals liquids. This will enable the interpretation of low-frequency spectra of more complex liquids.

Towards spectrally selective catastrophic response

We study the large-amplitude response of classical molecules to electromagnetic radiation, showing the universality of the transition from linear to nonlinear response and breakup at sufficiently large amplitudes. We demonstrate that a range of models, from the simple harmonic oscillator to the successful Peyrard-Bishop-Dauxois type models of DNA, which include realistic effects of the environment (including damping and dephasing due to thermal fluctuations), lead to characteristic universal behavior: formation of domains of dissociation in driving force amplitude-frequency space, characterized by the presence of local boundary minima. We demonstrate that by simply following the progression of the resonance maxima in this space, while gradually increasing intensity of the radiation, one must necessarily arrive at one of these minima, i.e., a point where the ultrahigh spectral selectivity is retained. We show that this universal property, applicable to other oscillatory systems, is a consequence of the fact that these models belong to the fold catastrophe universality class of Thom's catastrophe theory. This in turn implies that for most biostructures, including DNA, high spectral sensitivity near the onset of the denaturation processes can be expected. Such spectrally selective molecular denaturation could find important applications in biology and medicine.

Experimental observation of nanophase segregation in aqueous salt solutions around the predicted liquid-liquid transition in water

The liquid-liquid transition in supercooled liquid water, predicted to occur around 220 K, is controversial due to the difficulty of studying it caused by competition from ice crystallization (the so-called "no man's land"). In aqueous solutions, it has been predicted to give rise to phase separation on a nanometer scale between a solute-rich high-density phase and a water-rich low-density phase. Here we report direct experimental evidence for the formation of a nanosegregated phase in eutectic aqueous solutions of LiCl and LiSCN where the presence of crystalline water can be experimentally excluded. Femtosecond infrared and Raman spectroscopies are used to determine the temperature-dependent structuring of water, the solvation of the SCN^- anion, and the size of the phase segregated domains.

Low-Frequency Vibrational Motions of Polystyrene in Carbon Tetrachloride: Comparison with Model Monomer and Dependence on Concentration and Molecular Weight

In this study, the low-frequency vibrational dynamics of polystyrene (PS) in CCl₄ was investigated by femtosecond Raman-induced Kerr effect spectroscopy. Ethylbenzene (EBz) was also investigated as a model monomer of the polymer to elucidate the unique dynamical features of PS in solution. The broadened low-frequency spectrum of the PS/CCl₄ in the frequency region below 150 cm^-1 is significantly different from that of the EBz/CCl₄. Difference spectra between the PS or EBz solutions and neat CCl₄, normalized to an internal vibrational mode of CCl₄, clearly show a much lower spectral intensity for the PS/CCl₄ than the EBz/CCl₄ in the low-frequency region below ca. 20 cm^-1. This indicates that translational motions are suppressed in the PS/CCl₄ compared to the EBz/CCl₄. Moreover, the high-frequency motion at ca. 70 cm^-1, mainly due to phenyl ring librations, occurs at higher frequency in PS (78 cm^-1) than EBz (65 cm^-1). In addition, the results of concentration-dependent experiments show that the first moment (M₁) of the low-frequency difference spectra of both PS/CCl₄ and EBz/CCl₄ is almost independent of the concentration. The molecular weight dependence of the low-frequency spectrum in the PS/CCl₄ shows that the M₁ value of the low-frequency spectral band of PS shifts to higher frequencies when the molecular weight of PS increases up to M_w = ∼1000, which corresponds approximately to the decamer, and then remains constant upon further increasing the molecular weight.

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.

Anticancer drug impact on DNA – a study by neutron spectroscopy coupled with synchrotron-based FTIR and EXAFS

Complementary information on drug–DNA interplay has been achieved for Pt/Pd anticancer agents, by a combined QENS, SR-FTIR-ATR and EXAFS approach.

Terahertz: dictating the frequency of life. Do macromolecular vibrational modes impose thermal limitations on terrestrial life?

Conditions on exoplanets include elevated temperatures and pressures. The response of carbon-based biological macromolecules to such conditions is then relevant to the viability of life. The capacity of proteins and ribozymes to catalyse reactions or bind receptors, and nucleic acids to convey information, depends on them sampling different conformational states. These are determined by macromolecular vibrational states, or phonon modes, accessible using terahertz (THz: 10¹²Hz) absorption spectroscopy. THz spectra of biological macromolecules exhibit broad absorption at approximately 6 THz (equating to approx. 280 K) corresponding to dense transitions between phonon modes. There are also troughs at approximately 10 THz (approx. 500 K) implying diminishing numbers of available conformational states at higher temperatures; hence, fewer routes by which biochemical processes can be realized, as equilibrium is approached. Could this conformational bottleneck hinder the operation of biological macromolecules at higher temperatures? We suggest that the troughs at approximately 10 THz in absorbance spectra indicate that the hydrogen bonds, charge interactions and geometry of biological macromolecules associated with terrestrial life impose fundamental vibrational properties that could limit the upper temperature at which they may function.

Uncovering the Connection Between Low-Frequency Dynamics and Phase Transformation Phenomena in Molecular Solids

The low-frequency motions of molecules in the condensed phase have been shown to be vital to a large number of physical properties and processes. However, in the case of disordered systems, it is often difficult to elucidate the atomic-level details surrounding these phenomena. In this work, we have performed an extensive experimental and computational study on the molecular solid camphor, which exhibits a rich and complex structure-dynamics relationship, and undergoes an order-disorder transition near ambient conditions. The combination of x-ray diffraction, variable temperature and pressure terahertz time-domain spectroscopy, ab initio molecular dynamics, and periodic density functional theory calculations enables a complete picture of the phase transition to be obtained, inclusive of mechanistic, structural, and thermodynamic phenomena. Additionally, the low-frequency vibrations of a disordered solid are characterized for the first time with atomic-level precision, uncovering a clear link between such motions and the phase transformation. Overall, this combination of methods allows for significant details to be obtained for disordered solids and the associated transformations, providing a framework that can be directly applied for a wide range of similar systems.

Ultrafast bridge planarization in donor-π-acceptor copolymers drives intramolecular charge transfer

Donor-π-acceptor conjugated polymers form the material basis for high power conversion efficiencies in organic solar cells. Large dipole moment change upon photoexcitation via intramolecular charge transfer in donor-π-acceptor backbone is conjectured to facilitate efficient charge-carrier generation. However, the primary structural changes that drive ultrafast charge transfer step have remained elusive thereby limiting a rational structure-function correlation for such copolymers. Here we use structure-sensitive femtosecond stimulated Raman spectroscopy to demonstrate that π-bridge torsion forms the primary reaction coordinate for intramolecular charge transfer in donor-π-acceptor copolymers. Resonance-selective Raman snapshots of exciton relaxation reveal rich vibrational dynamics of the bridge modes associated with backbone planarization within 400 fs, leading to hot intramolecular charge transfer state formation while subsequent cooling dynamics of backbone-centric modes probe the charge transfer relaxation. Our work establishes a phenomenological gating role of bridge torsions in determining the fundamental timescale and energy of photogenerated carriers, and therefore opens up dynamics-based guidelines for fabricating energy-efficient organic photovoltaics.

Kinetic theory for DNA melting with vibrational entropy

By treating DNA as a vibrating nonlinear lattice, an activated kinetic theory for DNA melting is developed to capture the breakage of the hydrogen bonds and subsequent softening of torsional and bending vibration modes. With a coarse-grained lattice model, we identify a key bending mode with GHz frequency that replaces the hydrogen vibration modes as the dominant out-of-phase phonon vibration at the transition state. By associating its bending modulus to a universal in-phase bending vibration modulus at equilibrium, we can hence estimate the entropic change in the out-of-phase vibration from near-equilibrium all-atom simulations. This and estimates of torsional and bending entropy changes lead to the first predictive and sequence-dependent theory with good quantitative agreement with experimental data for the activation energy of melting of short DNA molecules without intermediate hairpin structures.

Evaluating the role of coherent delocalized phonon-like modes in DNA cyclization

The innate flexibility of a DNA sequence is quantified by the Jacobson-Stockmayer's J-factor, which measures the propensity for DNA loop formation. Recent studies of ultra-short DNA sequences revealed a discrepancy of up to six orders of magnitude between experimentally measured and theoretically predicted J-factors. These large differences suggest that, in addition to the elastic moduli of the double helix, other factors contribute to loop formation. Here, we develop a new theoretical model that explores how coherent delocalized phonon-like modes in DNA provide single-stranded "flexible hinges" to assist in loop formation. We combine the Czapla-Swigon-Olson structural model of DNA with our extended Peyrard-Bishop-Dauxois model and, without changing any of the parameters of the two models, apply this new computational framework to 86 experimentally characterized DNA sequences. Our results demonstrate that the new computational framework can predict J-factors within an order of magnitude of experimental measurements for most ultra-short DNA sequences, while continuing to accurately describe the J-factors of longer sequences. Further, we demonstrate that our computational framework can be used to describe the cyclization of DNA sequences that contain a base pair mismatch. Overall, our results support the conclusion that coherent delocalized phonon-like modes play an important role in DNA cyclization.

Water Dynamics in the Hydration Shells of Biomolecules

The structure and function of biomolecules are strongly influenced by their hydration shells. Structural fluctuations and molecular excitations of hydrating water molecules cover a broad range in space and time, from individual water molecules to larger pools and from femtosecond to microsecond time scales. Recent progress in theory and molecular dynamics simulations as well as in ultrafast vibrational spectroscopy has led to new and detailed insight into fluctuations of water structure, elementary water motions, electric fields at hydrated biointerfaces, and processes of vibrational relaxation and energy dissipation. Here, we review recent advances in both theory and experiment, focusing on hydrated DNA, proteins, and phospholipids, and compare dynamics in the hydration shells to bulk water.

Spectrum of Slow and Super-Slow (Picosecond to Nanosecond) Water Dynamics around Organic and Biological Solutes

Water dynamics in the solvation shell of solutes plays a very important role in the interaction of biomolecules and in chemical reaction dynamics. However, a selective spectroscopic study of the solvation shell is difficult because of the interference of the solute dynamics. Here we report on the observation of heavily slowed down water dynamics in the solvation shell of different solutes by measuring the low-frequency spectrum of solvation water, free from the contribution of the solute. A slowdown factor of ∼50 is observed even for relatively low concentrations of the solute. We go on to show that the effect can be generalized to different solutes including proteins.

Perspective: Structure and ultrafast dynamics of biomolecular hydration shells

The structure and function of biomolecules can be strongly influenced by their hydration shells. A key challenge is thus to determine the extent to which these shells differ from bulk water, since the structural fluctuations and molecular excitations of hydrating water molecules within these shells can cover a broad range in both space and time. Recent progress in theory, molecular dynamics simulations, and ultrafast vibrational spectroscopy has led to new and detailed insight into the fluctuations of water structure, elementary water motions, and electric fields at hydrated biointerfaces. Here, we discuss some central aspects of these advances, focusing on elementary molecular mechanisms and processes of hydration on a femto- to picosecond time scale, with some special attention given to several issues subject to debate.

Phonon-like Hydrogen-Bond Modes in Protic Ionic Liquids

Gigahertz- to terahertz-frequency infrared and Raman spectra contain a wealth of information concerning the structure, intermolecular forces, and dynamics of ionic liquids. However, these spectra generally have a large number of contributions ranging from slow diffusional modes to underdamped librations and intramolecular vibrational modes. This makes it difficult to isolate effects such as the role of Coulombic and hydrogen-bonding interactions. We have applied far-infrared and ultrafast optical Kerr effect spectroscopies on carefully selected ions with a greater or lesser degree of symmetry in order to isolate spectral signals of interest. This has allowed us to demonstrate the presence of longitudinal and transverse optical phonon modes and a great similarity of alkylammonium-based protic ionic liquids to liquid water. The data show that such phonon modes will be present in all ionic liquids, requiring a reinterpretation of their spectra.

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.

Far infrared spectroscopy of hydrogen bonding collective motions in complex molecular systems

Far infrared spectroscopy as a tool for the study of inter and intramolecular interactions in complex molecular structures.

The significance of the amorphous potential energy landscape for dictating glassy dynamics and driving solid-state crystallisation

We show clear evidence for a theory proposing that the shape and structure of the PES is the fundamental factor underlying the dynamics at temperatures below the glass transition.

Concomitant polymorphism and the martensitic-like transformation of an organic crystal

Terahertz vibrational spectroscopy and solid-state density functional theory together reveal the true nature of a pseudo-continuous crystalline polymorphic phase transition.