We present terahertz (THz) measurements of salt solutions that shed new light on the controversy over whether salts act as kosmotropes (structure makers) or chaotropes (structure breakers), which enhance or reduce the solvent order, respectively. We have carried out precise measurements of the concentration-dependent THz absorption coefficient of 15 solvated alkali halide salts around 85 cm(-1) (2.5 THz). In addition, we recorded overview spectra between 30 and 300 cm(-1) using a THz Fourier transform spectrometer for six alkali halides. For all solutions we found a linear increase of THz absorption compared to pure water (THz excess) with increasing solute concentration. These results suggest that the ions may be treated as simple defects in an H-bond network. They therefore cannot be characterized as either kosmotropes or chaotropes. Below 200 cm(-1), the observed THz excess of all salts can be described by a linear superposition of the water absorption and an additional absorption that is attributed to a rattling motion of the ions within the water network. By providing a comprehensive set of data for different salt solutions, we find that the solutions can all be very well described by a model that includes damped harmonic oscillations of the anions and cations within the water network. We find this model predicts the main features of THz spectra for a variety of salt solutions. The assumption of the existence of these ion rattling motions on sub-picosecond time scales is supported by THz Fourier transform spectroscopy of six alkali halides. Above 200 cm(-1) the excess is interpreted in terms of a change in the wing of the water network librational mode. Accompanying molecular dynamics simulations using the TIP3P water model support our conclusion and show that the fast sub-picosecond motions of the ions and their surroundings are almost decoupled. These findings provide a complete description of the solute-induced changes in the THz solvation dynamics for the investigated salts. Our results show that THz spectroscopy is a powerful experimental tool to establish a new view on the contributions of anions and cations to the structuring of water.

Frequency-resolved communication maps provide a coarse-grained picture of energy transport in nanoscale systems. We calculate communication maps for homodimeric hemoglobin from Scapharca inaequivalvis and sample them to elucidate energy transfer pathways between the binding sites and other parts of the protein with focus on the role of the cluster of water molecules at the interface between the globules. We complement analysis of communication maps with molecular simulations of energy flow. Both approaches reveal that excess energy in one heme flows mainly to regions of the interface where early hydrogen bond rearrangements occur in the allosteric transition. In particular, energy is carried disproportionately by the water molecules, consistent with the larger thermal conductivity of water compared to proteins.

Water confined in proteins exhibits dynamics distinct from the dynamics of water in the bulk or near the surface of a biomolecule. We examine the water dynamics at the interface of the two globules of the homodimeric hemoglobin from Scapharca inaequivalvis (HbI) by molecular dynamics (MD) simulations, with focus on water-protein hydrogen bond lifetimes and rotational anisotropy of the interfacial waters. We find that relaxation of the waters at the interface of both deoxy- and oxy-HbI, which contain a cluster of 17 and 11 interfacial waters, respectively, is well described by stretched exponentials with exponents from 0.1 to 0.6 and relaxation times of tens to thousands of picoseconds. The interfacial water molecules of oxy-HbI exhibit slower rotational relaxation and hydrogen bond rearrangement than those of deoxy-HbI, consistent with an allosteric transition from unliganded to liganded conformers involving the expulsion of several water molecules from the interface. Though the interfacial waters are translationally and rotationally static on the picosecond time scale, they contribute to fast communication between the globules via vibrations. We find that the interfacial waters enhance vibrational energy transport across the interface by ≈10%.

Accurate prediction of protein-ligand binding affinities is essential for hit-to-lead optimization and virtual screening. The reliability of scoring functions can be improved by including quantum effects. Here, we demonstrate the ranking power of the semiempirical quantum mechanics (SQM)/implicit solvent (COSMO) scoring function by using a challenging set of 10 inhibitors binding to carbonic anhydrase II through Zn in the active site. This new dataset consists of the high-resolution (1.1-1.4 Å) crystal structures and experimentally determined inhibitory constant (K ) values. It allows for evaluation of the common approximations, such as representing the solvent implicitly or by using a single target conformation combined with a set of ligand docking poses. SQM/COSMO attained a good correlation of R of 0.56-0.77 with the experimental inhibitory activities, benefiting from careful handling of both noncovalent interactions (e.g. charge transfer) and solvation. This proof-of-concept study of SQM/COSMO ranking for metalloprotein-ligand systems demonstrates its potential for hit-to-lead applications.