Intensities of the Raman bands in the low-frequency spectra of DNA with light and heavy counterions

The approach for calculation of the mode intensities of DNA conformational vibrations in the Raman spectra is developed. It is based on the valence-optic theory and the model for description of conformational vibrations of DNA with counterions. The calculations for Na- and Cs-DNA low-frequency Raman spectra show that the vibrations of DNA backbone chains near 15 cm-1 have the greatest intensity. In the spectrum of Na-DNA at frequency range upper than 40 cm-1 the modes of H-bond stretching in base pairs have the greatest intensities, while the modes of ion-phosphate vibrations have the lowest intensity. In Cs-DNA spectra at this frequency range the mode of ion-phosphate vibrations is prominent. Its intensity is much higher than the intensities of Na-DNA modes of this spectra range. Other modes of Cs-DNA have much lower intensities than in the case of Na-DNA. The comparison of our calculations with the experimental data shows that developed approach gives the understanding of the sensitivity of DNA low-frequency Raman bands to the neutralization of the double helix by light and heavy counterions.

💡 Research Summary

The paper presents a theoretical framework for calculating Raman intensities of low‑frequency vibrational modes in DNA solutions containing either light (Na⁺) or heavy (Cs⁺) counter‑ions. Building on the authors’ previously developed four‑mass model of DNA conformational dynamics (phosphate, base, deoxyribose, and counter‑ion as discrete masses), the study incorporates the valence‑optic (or “bond‑polarizability”) theory to evaluate how the polarizability tensor of the DNA monomer changes with normal coordinates of vibration.

Key methodological steps:

1. Four‑mass model – Treats each nucleotide as a set of coupled pendulums (bases and sugars) attached to phosphate groups, with counter‑ions tethered to phosphates forming an ion‑phosphate lattice. The model yields seven normal modes: symmetric and antisymmetric base‑pendulum motions (B, B′), symmetric/antisymmetric ion‑phosphate vibrations (Ion, Ion′), hydrogen‑bond stretching modes (H, HS), and other low‑frequency deformations (S, S′).
2. Valence‑optic theory (zero‑order approximation) – Assumes that only the intrinsic polarizabilities of chemical bonds change with bond angles; contributions from neighboring bonds, phosphate motions, and ion motions are neglected. The polarizability tensor of each nucleoside is expressed as a sum of bond tensors, rotated into the molecular frame using small‑angle rotation matrices that depend on the generalized coordinates (θ₁, θ₂) of the model.
3. Derivation of intensity expression – Raman intensity for a mode m is given by the standard semi‑classical formula Iₘ ∝ (ν₀‑νₘ)⁴·|∂α/∂Qₘ|²·Qₘ². By relating the derivatives of the polarizability tensor to the small rotation angles and using the known relationship between generalized and normal coordinates, the authors obtain an explicit analytical expression (Eq. 19) that separates symmetric and antisymmetric contributions.

Using parameters (frequencies, amplitudes) previously calculated for Na‑DNA and Cs‑DNA, the authors compute Raman intensities for each of the seven modes. The main findings are:

• Backbone (B, B′) modes near 15 cm⁻¹ have the largest calculated intensities for both Na‑ and Cs‑DNA, confirming that backbone motions dominate the very low‑frequency region.
• For Na‑DNA, modes above 40 cm⁻¹ are dominated by hydrogen‑bond stretching (H, HS). The ion‑phosphate modes (Ion, Ion′) are weak (degenerate around 90 cm⁻¹) and contribute little to the observable spectrum.
• For Cs‑DNA, the ion‑phosphate modes become prominent, especially around 100 cm⁻¹, where their calculated intensity exceeds that of the H‑bond modes. This reflects the mass‑dependent softening of the ion‑phosphate lattice and the larger vibrational amplitudes associated with heavy counter‑ions. Other modes (S, S′) have markedly lower intensities in Cs‑DNA compared with Na‑DNA.

The theoretical intensity patterns are compared with experimental Raman spectra of aqueous DNA solutions. The experimental data show a strong band near 15 cm⁻¹ (backbone), a dominant H‑bond band above 40 cm⁻¹ for Na‑DNA, and a markedly enhanced band around 100 cm⁻¹ when Na⁺ is replaced by Cs⁺. These observations align well with the calculated intensities, supporting the validity of the model.

In conclusion, the study demonstrates that (i) a semi‑classical, bond‑polarizability approach can successfully predict Raman intensities of complex macromolecular vibrations, and (ii) the type of counter‑ion dramatically influences the low‑frequency Raman signature of DNA, primarily by modulating the strength of ion‑phosphate lattice vibrations. The work provides a quantitative tool for interpreting DNA Raman spectra in terms of counter‑ion composition and suggests avenues for extending the model to other metal ions, temperature effects, and solvent environments.