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 inter-protein 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 allencompassing 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 semi-empirical 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 website (http://frequencymap.org) for vibrational spectroscopic maps is presented. We anticipate that a combination of the vibrational frequency map approach and state-of-theart 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.

Isotope-edited two-dimensional infrared spectroscopy has been used to characterize the conformational heterogeneity of the beta-hairpin peptide TrpZip2 (TZ2) across its thermal unfolding transition. Four isotopologues were synthesized to probe hydrogen bonding and solvent exposure of the beta-turn (K8), the N-terminus (S1), and the midstrand region (T10 and T3T10). Isotope-shifts, 2D lineshapes, and other spectral changes to the amide I 2D IR spectra of labeled TZ2 isotopologues were observed as a function of temperature. Data were interpreted on the basis of structure-based spectroscopic modeling of conformers obtained from extensive molecular dynamics simulations. The K8 spectra reveal two unique turn geometries, the type I' beta-turn observed in the NMR structure, and a less populated disordered or bulged loop. The data indicate that structures at low temperature resemble the folded NMR structure with typical cross-strand hydrogen bonds, although with a subpopulation of misformed turns. As the temperature is raised from 25 to 85 degrees C, the fraction of population with a type I' turn increases, but the termini also fray. Hydrogen bonding contacts in the midstrand region remain at all temperatures although with increasing thermal disorder. Our data show no evidence of an extended chain or random coil state for the TZ2 peptide at any temperature. The methods demonstrated here offer an approach to characterizing conformational variation within the folded or unfolded states of proteins and peptides.

The development of technologies for nuclear reactors based on molten salts has seen a big resurgence. The success of thermodynamic models for these hinges in part on our ability to predict at the atomistic level the behavior of pure salts and their mixtures under a range of conditions. In this letter, we present high-energy X-ray scattering experiments and molecular dynamics simulations that describe the molten structure of mixtures of MgCl2 and KCl. As one would expect, KCl is a prototypical salt in which structure is governed by simple charge alternation. In contrast, MgCl2 and its mixtures with KCl display more complex correlations including intermediate-range order and the formation of Cl–-decorated Mg2+ chains. A thorough computational analysis suggests that intermediate-range order beyond charge alternation may be traced to correlations between these chains. An analysis of the coordination structure for Mg2+ ions paints a more complex picture than previously understood, with multiple accessible states of distinct geometries.

This article addresses the non-Debye−Waller temperature behavior in the intermediate range order for molten MgCl 2 and its mixtures with KCl from a theory, Molecular Dynamics, and experimental X-ray scattering perspective and puts these findings in the context of discussions and controversies extending at least four decades. We find that these liquids are defined by two structural motifs. The first motif is associated with chains of positive−negative charge alternation; the second motif, which results in a prepeak in the structure function S(q), is associated with the interaction of Mg 2+ and Cl − ions that do not belong to the same charge alternation chain or aggregate. Our complementary X-ray scattering and computational results provide a quantitative explanation for the increase in intensity of the prepeak with temperature as opposed to the behavior of other peaks following normal Debye−Waller behavior. Temperature has opposite effects on the prevalence of each of the two structural motifs, and the enhancement of one pattern appears to be at the detriment of the other. Whereas the intensity in S(q) associated with the charge alternation motif is diminished at higher temperature, the opposite is true for the prepeak associated with intermediate range order due to the second structural motif.

Understanding the structure of water near cell membranes is crucial for characterizing water-mediated events such as molecular transport. To obtain structural information of water near a membrane, it is useful to have a surface-selective technique that can probe only interfacial water molecules. One such technique is vibrational sum-frequency generation (VSFG) spectroscopy. As model systems for studying membrane headgroup/water interactions, in this paper we consider lipid and surfactant monolayers on water. We adopt a theoretical approach combining molecular dynamics simulations and phase-sensitive VSFG to investigate water structure near these interfaces. Our simulated spectra are in qualitative agreement with experiments and reveal orientational ordering of interfacial water molecules near cationic, anionic, and zwitterionic interfaces. OH bonds of water molecules point toward an anionic interface leading to a positive VSFG peak, whereas the water hydrogen atoms point away from a cationic interface leading to a negative VSFG peak. Coexistence of these two interfacial water species is observed near interfaces between water and mixtures of cationic and anionic lipids, as indicated by the presence of both negative and positive peaks in their VSFG spectra. In the case of a zwitterionic interface, OH orientation is toward the interface on the average, resulting in a positive VSFG peak.

We present a mixed quantum-classical model for studying the amide I vibrational dynamics (predominantly CO stretching) in peptides and proteins containing proline. There are existing models developed for determining frequencies of and couplings between the secondary amide units. However, these are not applicable to proline because this amino acid has a tertiary amide unit. Therefore, a new parametrization is required for infrared-spectroscopic studies of proteins that contain proline, such as collagen, the most abundant protein in humans and animals. Here, we construct the electrostatic and dihedral maps accounting for solvent and conformation effects on frequency and coupling for the proline unit. We examine the quality and the applicability of these maps by carrying out spectral simulations of a number of peptides with proline in D 2 O and compare with experimental observations.

The peptide amide-I vibration of a proline turn encodes information on the turn structure. In this study, FTIR, two-dimensional IR spectroscopy and molecular dynamics simulations were employed to characterize the varying turn conformations that exist in the GVGX(L)PGVG family of disordered peptides. This analysis revealed that changing the size of the side chain at the X amino acid site from Gly to Ala to Val substantially alters the conformation of the peptide. To quantify this effect, proline peak shifts and intensity changes were compared to a structure-based spectroscopic model. These simulated spectra were used to assign the population of type-II β turns, bulged turns, and irregular β turns for each peptide. Of particular interest was the Val variant commonly found in the protein elastin, which contained a 25% population of irregular β turns containing two peptide hydrogen bonds to the proline C═O.

We present a theory for ion pair dissociation and association, motivated by the concepts of Marcus theory of electron transfer. Despite the extensive research on ion-pairing in many chemical and biological processes, much can be learned from the exploration of collective reaction coordinates. To this end, we explore two reaction coordinates, ion pair distance and coordination number. The study of the correlation between these reaction coordinates provides a new insight into the mechanism and kinetics of ion pair dissociation and association in water. The potential of mean force on these 2D surfaces computed from molecular dynamics simulations of different monovalent ion pairs reveal a Marcus-like mechanism for ion-pairing: Water molecules rearrange forming an activated coordination state prior to ion pair dissociation or association, followed by relaxation of the coordination state due to further water rearrangement. Like Marcus theory, we find the existence of an inverted region where the transition rates are slower with increasing exergonicity. This study provides a new perspective for the future investigations of ion-pairing and transport.