Terahertz (THz) spectroscopic investigations of condensed-phase biological samples are reviewed ranging from the simple crystalline forms of amino acids, carbohydrates and polypeptides to the more complex aqueous forms of small proteins, DNA and RNA. Vibrationally resolved studies of crystalline samples have revealed the exquisite sensitivity of THz modes to crystalline order, temperature, conformational form, peptide sequence and local solvate environment and have given unprecedented measures of the binding force constants and anharmonic character of the force fields, properties necessary to improve predictability but not readily obtainable using any other method. These studies have provided benchmark vibrational data on extended periodic structures for direct comparisons with classical (CHARMm) and quantum chemical (density functional theory) theories. For the larger amorphous and/or aqueous phase samples, the THz modes form a continuum-like absorption that arises because of the full accessibility to conformational space and/or the rapid time scale for inter-conversion in these environments. Despite severe absorption by liquid water, detailed investigations have uncovered the photo- and hydration-induced conformational flexibility of proteins, the solvent shell depth of the water/biomolecule boundary layers and the solvent reorientation dynamics occurring in these interfacial layers that occur on sub-picosecond time scales. As such, THz spectroscopy has enhanced and extended the accessibility to intermolecular forces, length- and timescales important in biological structure and activity.

High-resolution terahertz absorption spectra (0.06-3 THz) have been obtained at 4.2 K for three crystalline forms of trialanine [H2+-(Ala)3-O-]. The crystal structures differ in their beta-sheet forms (parallel vs antiparallel) and in their water composition (hydrated vs dehydrated antiparallel beta-sheet). The spectra are nearly vibrationally resolved, with little absorption below 1 THz. In sharp contrast to observations made in the mid-IR region, the spectral patterns of all three forms are qualitatively different, illustrating the extreme sensitivity to changes in the intermolecular hydrogen-bonding networks that stabilize peptide crystals. Predictions obtained from a classical force field model (CHARMM) and density functional theory (DFT/PW91) for periodic solids are compared with the X-ray structural data and the terahertz absorption spectra. In general, the results for the parallel beta-sheet are in better agreement with experiment than those for the antiparallel beta-sheet. For all three structures, however, most hydrogen bond distances are underestimated at both levels of theory, and the predicted absorption features are significantly red-shifted for the two antiparallel beta-sheet structures. Moreover, the nuclear motions predicted at the two levels of theory are qualitatively different. These results indicate that the PW91 functional is not sufficient to treat the weak intersheet hydrogen bonding present in the different beta-sheet forms and strongly suggest the need for improved force field models that include three-atom hydrogen-bonding terms for periodic solids.

Terahertz (THz) vibrational modes are characterized by nonlocal, collective molecular motions which are relevant to conformational changes and molecular functions in biological systems. We have investigated the THz spectra of a set of small bionanotubes which can serve as very simple models of membrane pores, and have examined the character of the THz modes which can impact transport processes. In this work, THz spectra of the crystalline VA class dipeptide nanotubes were calculated at both the harmonic and vibrational self-consistent field (VSCF) level using the CHARMM22 force field with periodic boundary conditions. Comparison of the calculated THz spectra against the experimental spectra revealed that the VSCF corrections generally improved the predictions in the low-frequency region. The improvements were especially manifested in the overall blue-shifts of the VSCF frequencies relative to the harmonic values, and blue shifts were attributed to the overall positive coupling strengths in all systems. Closer examination of the motions in the most significantly coupled normal mode pairs leads us to propose that, when two similar side-chain squeezing modes are coupled, the rapidly increased van der Waals interactions can lead to a stiffening of the effective potential, which in turn leads to the observed blue-shifts. However, we also noted that when the side-chain atoms become unphysically proximate and the van der Waals repulsion becomes too large, the VSCF calculations tend to deviate in the high frequency region and for the system of l-isoleucyl-l-valine. In addition, normal-mode analysis revealed a series of channel-breathing motions in all systems except l-valyl-l-alanine. We show that the inner products of the backbone vibrations between these channel-breathing motions divided the remaining VA class dipeptide systems into two subgroups. It is suggested that these modes may facilitate a pathway for the guest molecule absorption, substitution and removal in the VA class dipeptide nanotubes. Normal mode analysis also demonstrated that the THz motions may contribute to the pore permeability either directly by changing the pore size, or indirectly by affecting the solvent-host effective potentials.