Vibrational Spectroscopic Map, Vibrational Spectroscopy, and Intermolecular Interaction.

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 interprotein 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 all-encompassing 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 semiempirical 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 Web site (http://frequencymap.org) for vibrational spectroscopic maps is presented. We anticipate that a combination of the vibrational frequency map approach and state-of-the-art 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.

Coupled local mode method for simulating vibrational spectroscopy.

Experimental and theoretical studies have highlighted protonated water clusters (PWCs) as important models of the excess proton in aqueous systems. A significant focus has been characterizing the spectral signatures associated with different excess proton solvation motifs. Accurate vibrational frequency calculations are crucial for connecting the measured spectra to the structure of PWCs. In this paper, we extend and characterize a coupled local mode (CLM) approach for calculating the infrared spectra of PWCs using the H+(H2O)4 cluster as a benchmark system. The CLM method is relatively low cost and incorporates the anharmonicity and coupling of OH vibrations. Here, we demonstrate the accuracy of the technique compared to experiments. We also illustrate the dependence of calculated spectral features on the underlying electronic structure theory and basis sets used in the local mode frequency and coupling calculations.

Nonlinear measurements of kinetics and generalized dynamical modes. II. Application to a simulation of solvation dynamics in an ionic liquid.

Solvation dynamics in ionic liquids show features that are often associated with supercooled liquids, including "stretched" nonexponential relaxation. To better understand the mechanism behind the stretching, the nonlinear mode-correlation methods proposed in Paper I [S. R. Hodge and M. A. Berg, J. Chem. Phys. 155, 024122 (2021)] are applied to a simulation of a prototypical ionic liquid. A full Green's function is recovered. In addition, specific tests for non-Gaussian dynamics are made. No deviations from Gaussian dynamics are found. This finding is incompatible with rate heterogeneity as a cause of the nonexponential relaxation and appears to be in conflict with an earlier multidimensional analysis of the same data. Although this conflict is not resolved here, this work does demonstrate the practicality of mode-correlation analysis in the face of finite datasets and calculations.

Self-assembly of hydrogen-bonded two-dimensional quasicrystals.

The process of molecular self-assembly on solid surfaces is essentially one of crystallization in two dimensions, and the structures that result depend on the interplay between intermolecular forces and the interaction between adsorbates and the underlying substrate. Because a single hydrogen bond typically has an energy between 15 and 35 kilojoules per mole, hydrogen bonding can be a strong driver of molecular assembly; this is apparent from the dominant role of hydrogen bonding in nucleic-acid base pairing, as well as in the secondary structure of proteins. Carboxylic acid functional groups, which provide two hydrogen bonds, are particularly promising and reliable in creating and maintaining surface order, and self-assembled monolayers of benzoic acids produce structure that depends on the number and relative placement of carboxylic acid groups. Here we use scanning tunnelling microscopy to study self-assembled monolayers of ferrocenecarboxylic acid (FcCOOH), and find that, rather than producing dimeric or linear structures typical of carboxylic acids, FcCOOH forms highly unusual cyclic hydrogen-bonded pentamers, which combine with simultaneously formed FcCOOH dimers to form two-dimensional quasicrystallites that exhibit local five-fold symmetry and maintain translational and rotational order (without periodicity) for distances of more than 400 ångströms.

Counter Cations Affect Transport in Aqueous Hydroxide Solutions with Ion Specificity.

The anomalously high mobility of hydroxide and hydronium ions in aqueous solutions is related to proton transfer and structural diffusion. The role of counterions in these solutions, however, is often considered to be negligible. Herein, we explore the impact of alkali metal counter cations on hydroxide solvation and mobility. Impedance measurements demonstrate that hydroxide mobility is attenuated by lithium relative to sodium and potassium. These results are explained by ab initio molecular dynamics simulations and experimental vibrational hydration shell spectroscopy, which reveal substantially stronger ion pairing between OH- and Li+ than with other cations. Hydration shell spectra and theoretical vibrational frequency calculations together imply that lithium and sodium cations have different effects on the delocalization of water protons donating a hydrogen bond to hydroxide. Specifically, lithium leads to enhanced proton delocalization compared with sodium. However, proton delocalization and the overall diffusion process are not necessarily correlated.

Structural Polymorphism as the Result of Kinetically Controlled Self-Assembly.

Traditionally, the goal of self-assembly and supramolecular chemistry is to engineer an equilibrium structure with a desired geometry and functionality; this is achieved through careful choice of molecular monomers, growth conditions, and substrate. Supramolecular assemblies produced under nonequilibrium conditions, in contrast, can form metastable structures with conformations quite different from those accessible in equilibrium self-assembly. The study of nonequilibrium growth of clusters potentially impacts the study of nucleation in atmospheric aerosols, nucleation in organic crystallization, and mesoscale organization for systems ranging from biological molecules to molecular electronics. In our experiments, we prepare surface monolayers of small organic and organometallic molecules through direct injection of a solution onto a substrate in high vacuum. During this process, the rapid evaporation of small solution droplets in high vacuum can lead to nonequilibrium growth conditions. The resulting structures are then characterized by scanning tunneling microscopy. Among the features observed in these experiments are cyclic, hydrogen-bonded pentamers. For carboxylic acids, the two-molecule ring dimer is the common binding motif. Large, cyclic hydrogen-bonded systems are uncommon, especially so for rings with five members. Despite this, pentagonal clusters appear to be a general phenomenon for systems containing adjacent strong and weak hydrogen-bonding elements on five-member aromatic rings. Regular pentamers have been observed as metastable structures for ferrocenecarboxylic acid, indole-2-carboxylic acid, and isatin (1-H-indole-2,3-dione). Electronic structure calculations confirm the relative stability of these structures with respect to the dimer or catemer conformations which are observed in the solid-state crystal structures. For ferrocenecarboxylic acid, cyclic pentamers undergo further self-assembly, resulting in long-range order in conjunction with local 5-fold rotational symmetry. This system is the first reported self-assembled molecular quasicrystal, and it remains the only example of a hydrogen-bonded quasicrystal. This supramolecular structure forms as a result of the cocrystallization of hydrogen-bonded cyclic pentamers with intercalated molecular dimers. The shared bonding to a single dimer is responsible for locking the adjacent pentamers in specific distances and orientations, which produces the quasicrystal. Careful analysis of experimental data provides evidence that, in some cases, metastable clusters are formed in solution and then subsequently adsorb on the surface. This is a unusual mechanism for supramolecular assembly, and it has important implications for understanding questions in crystal growth, namely: what the initial stages of crystal growth are as molecules are first precipitating from solution; what role the solvent plays in determining crystal structure; and whether solvent-mediated clustering is important in the broader phenomenon of solid-state polymorphism.

Effect of pH and Salt on Surface pKa of Phosphatidic Acid Monolayers.

The pH-induced surface speciation of organic surfactants such as fatty acids and phospholipids in monolayers and coatings is considered to be an important factor controlling their interfacial organization and properties. Yet, correctly predicting the surface speciation requires the determination of the surface dissociation constants (surface pKa) of the protic functional group(s) present. Here, we use three independent methods-compression isotherms, surface tension pH titration, and infrared reflection-absorption spectroscopy (IRRAS)-to study the protonation state of dipalmitoylphosphatidic acid (DPPA) monolayers on water and NaCl solutions. By examining the molecular area expansion at basic pH, the pKa to remove the second proton of DPPA (surface pKa2) at the aqueous interface is estimated. In addition, utilizing IRRAS combined with density functional theory calculations, the vibrational modes of the phosphate headgroup were directly probed and assigned to understand DPPA charge speciation with increasing pH. We find that all three experimental techniques give consistent surface pKa2 values in good agreement with each other. Results show that a condensed DPPA monolayer has a surface pKa2 of 11.5, a value higher than previously reported (∼7.9-8.5). This surface pKa2 was further altered by the presence of Na+ cations in the aqueous subphase, which reduced the surface pKa2 from 11.5 to 10.5. It was also found that the surface pKa2 value of DPPA is modulated by the packing density (i.e., the surface charge density) of the monolayer, with a surface pKa2 as low as 9.2 for DPPA monolayers in the two-dimensional gaseous phase over NaCl solutions. The experimentally determined surface pKa2 values are also found to be in agreement with those predicted by Gouy-Chapman theory, validating these methods and proving that surface charge density is the driving factor behind changes to the surface pKa2.

Modeling Sequence-Dependent Peptide Fluctuations in Immunologic Recognition.

In cellular immunity, T cells recognize peptide antigens bound and presented by major histocompatibility complex (MHC) proteins. The motions of peptides bound to MHC proteins play a significant role in determining immunogenicity. However, existing approaches for investigating peptide/MHC motional dynamics are challenging or of low throughput, hindering the development of algorithms for predicting immunogenicity from large databases, such as those of tumor or genetically unstable viral genomes. We addressed this by performing extensive molecular dynamics simulations on a large structural database of peptides bound to the most commonly expressed human class-I MHC protein, HLA-A*0201. The simulations reproduced experimental indicators of motion and were used to generate simple models for predicting site-specific, rapid motions of bound peptides through differences in their sequence and chemical composition alone. The models can easily be applied on their own or incorporated into immunogenicity prediction algorithms. Beyond their predictive power, the models provide insight into how amino acid substitutions can influence peptide and protein motions and how dynamic information is communicated across peptides. They also indicate a link between peptide rigidity and hydrophobicity, two features known to be important in influencing cellular immune responses.

Electric-field-driven electron-transfer in mixed-valence molecules.

Molecular quantum-dot cellular automata is a computing paradigm in which digital information is encoded by the charge configuration of a mixed-valence molecule. General-purpose computing can be achieved by arranging these compounds on a substrate and exploiting intermolecular Coulombic coupling. The operation of such a device relies on nonequilibrium electron transfer (ET), whereby the time-varying electric field of one molecule induces an ET event in a neighboring molecule. The magnitude of the electric fields can be quite large because of close spatial proximity, and the induced ET rate is a measure of the nonequilibrium response of the molecule. We calculate the electric-field-driven ET rate for a model mixed-valence compound. The mixed-valence molecule is regarded as a two-state electronic system coupled to a molecular vibrational mode, which is, in turn, coupled to a thermal environment. Both the electronic and vibrational degrees-of-freedom are treated quantum mechanically, and the dissipative vibrational-bath interaction is modeled with the Lindblad equation. This approach captures both tunneling and nonadiabatic dynamics. Relationships between microscopic molecular properties and the driven ET rate are explored for two time-dependent applied fields: an abruptly switched field and a linearly ramped field. In both cases, the driven ET rate is only weakly temperature dependent. When the model is applied using parameters appropriate to a specific mixed-valence molecule, diferrocenylacetylene, terahertz-range ET transfer rates are predicted.

Hydrogen-bonded clusters of 1, 1'-ferrocenedicarboxylic acid on Au(111) are initially formed in solution.

Low-temperature scanning tunneling microscopy is used to observe self-assembled structures of ferrocenedicarboxylic acid (Fc(COOH)2) on the Au(111) surface. The surface is prepared by pulse-deposition of Fc(COOH)2 dissolved in methanol, and the solvent is evaporated before imaging. While the rows of hydrogen-bonded dimers that are common for carboxylic acid species are observed, the majority of adsorbed Fc(COOH)2 is instead found in six-molecule clusters with a well-defined and chiral geometry. The coverage and distribution of these clusters are consistent with a random sequential adsorption model, showing that solution-phase species are determinative of adsorbate distribution for this system under these reaction conditions.

Synthesis of a Neutral Mixed-Valence Diferrocenyl Carborane for Molecular Quantum-Dot Cellular Automata Applications.

The preparation of 7-Fc(+) -8-Fc-7,8-nido-[C2 B9 H10 ](-) (Fc(+) FcC2 B9 (-) ) demonstrates the successful incorporation of a carborane cage as an internal counteranion bridging between ferrocene and ferrocenium units. This neutral mixed-valence Fe(II) /Fe(III) complex overcomes the proximal electronic bias imposed by external counterions, a practical limitation in the use of molecular switches. A combination of UV/Vis-NIR spectroscopic and TD-DFT computational studies indicate that electron transfer within Fc(+) FcC2 B9 (-) is achieved through a bridge-mediated mechanism. This electronic framework therefore provides the possibility of an all-neutral null state, a key requirement for the implementation of quantum-dot cellular automata (QCA) molecular computing. The adhesion, ordering, and characterization of Fc(+) FcC2 B9 (-) on Au(111) has been observed by scanning tunneling microscopy.

Evidence of an unusual N-H···N hydrogen bond in proteins.

Many residues within proteins adopt conformations that appear to be stabilized by interactions between an amide N-H and the amide N of the previous residue. To explore whether these interactions constitute hydrogen bonds, we characterized the IR stretching frequencies of deuterated variants of proline and the corresponding carbamate, as well as the four proline residues of an Src homology 3 domain protein. The CδD2 stretching frequencies are shifted to lower energies due to hyperconjugation with Ni electron density, and engaging this density via protonation or the formation of the Ni+1-H···Ni interaction ablates this hyperconjugation and thus induces an otherwise difficult to explain blue shift in the C-D absorptions. Along with density functional theory calculations, the data reveal that the Ni+1-H···Ni interactions constitute H-bonds and suggest that they may play an important and previously underappreciated role in protein folding, structure, and function.

Dynamic allostery in the peptide/MHC complex enables TCR neoantigen selectivity.

The inherent cross-reactivity of the T cell receptor (TCR) is balanced by high specificity, which often manifests in confounding ways not easily interpretable from static structures. We show here that TCR discrimination between an HLA-A*03:01 (HLA-A3)-restricted public neoantigen derived from mutant PIK3CA and its wild-type (WT) counterpart emerges from motions within the HLA binding groove that vary with the identity of the peptide's first primary anchor. The motions form a dynamic gate that in the complex with the WT peptide impedes a large conformational change required for TCR binding. The more rigid neoantigen is insusceptible to this limiting dynamic, and with the gate open, is able to transit its central tryptophan residue underneath the peptide backbone to the contralateral side of the HLA-A3 peptide binding groove, facilitating TCR binding. Our findings reveal a novel mechanism driving TCR specificity for a cancer neoantigen that is rooted in the dynamic and allosteric nature of peptide/MHC-I complexes, with implications for resolving long-standing and often confounding questions about the determinants of T cell specificity.

Adsorption of diferrocenylacetylene on Au(111) studied by scanning tunneling microscopy.

Scanning tunneling microscopy images of diferrocenylacetylene (DFA) coadsorbed with benzene on Au(111) show individual and close-packed DFA molecules, either adsorbed alongside benzene or on top of a benzene monolayer. Images acquired over a range of positive and negative tip-sample bias voltages show a shift in contrast, with the acetylene linker appearing brighter than the ferrocenes at positive sample bias (where unoccupied states primarily contribute) and the reverse contrast at negative bias. Density functional theory was used to calculate the electronic structure of the gas-phase DFA molecule, and simulated images produced through two-dimensional projections of these calculations approximate the experimental images. The symmetry of both experimental and calculated molecular features for DFA rules out a cis adsorption geometry, and comparison of experiment to simulation indicates torsion around the inter-ferrocene axis between 90° and 180° (trans); the cyclopentadienyl rings are thus angled with respect to the surface.

Limitations of time-resolved fluorescence suggested by molecular simulations: assessing the dynamics of T cell receptor binding loops.

Time-resolved fluorescence anisotropy (TRFA) has a rich history in evaluating protein dynamics. Yet as often employed, TRFA assumes that the motional properties of a covalently tethered fluorescent probe accurately portray the motional properties of the protein backbone at the probe attachment site. In an extensive survey using TRFA to study the dynamics of the binding loops of a αß T cell receptor, we observed multiple discrepancies between the TRFA data and previously published results that led us to question this assumption. We thus simulated several of the experimentally probed systems using a protocol that permitted accurate determination of probe and protein time correlation functions. We found excellent agreement in the decays of the experimental and simulated correlation functions. However, the motional properties of the probe were poorly correlated with those of the backbone of both the labeled and unlabeled protein. Our results warrant caution in the interpretation of TRFA data and suggest further studies to ascertain the extent to which probe dynamics reflect those of the protein backbone. Meanwhile, the agreement between experiment and computation validates the use of molecular dynamics simulations as an accurate tool for exploring the molecular motion of T cell receptors and their binding loops.

Dynamical signature of abasic damage in DNA.

Time-dependent Stokes shift (TDSS) responses in proteins and DNA exhibit a broad range of long time scales (>10 ps) that are not present in bulk aqueous solution. The physical interpretation of the long TDSS time scales in biomolecular systems is a matter of considerable debate because of the many different components present in the sample (water, biomolecule, counterions), which have highly correlated motions and intrinsically different abilities to adapt to local perturbations. Here we use molecular dynamics (MD) simulations to show that the surprisingly slow (∼10 ns) TDSS response of coumarin 102 (C102), a base pair replacement, reflects a distinct dynamical signature for DNA damage. When the C102 molecule is covalently incorporated into DNA, an abasic site is created on the strand opposite the C102 probe. The abasic sugar exhibits a reversible interchange between intra- and extrahelical conformations that are kinetically stable on a nanosecond time scale. This conformational change, only possible in damaged DNA, was found to be responsible for the long time scales in the measured TDSS response. For the first time, a TDSS measurement has been attributed to a specific biomolecular motion. This finding directly contradicts the prevailing notion that the TDSS response in biomolecular contexts is dominated by hydration dynamics. It also suggests that TDSS experiments can be used to study ultrafast biomolecular dynamics that are inaccessible to other techniques.

Phosphate vibrations probe local electric fields and hydration in biomolecules.

The role of electric fields in important biological processes such as binding and catalysis has been studied almost exclusively by computational methods. Experimental measurements of the local electric field in macromolecules are possible using suitably calibrated vibrational probes. Here we demonstrate that the vibrational transitions of phosphate groups are highly sensitive to an electric field and show how that sensitivity can be quantified, allowing electric field measurements to be made in phosphate-containing biological systems without chemical modification.

Coupled Local-Mode Approach for the Calculation of Vibrational Spectra: Application to Protonated Water Clusters.

Spectroscopic studies of protonated water clusters (PWCs) have yielded enormous insights into the fundamental nature of the hydrated proton. Here, we introduce a new coupled local-mode (CLM) approach to calculate PWC OH stretch vibrational spectra. The CLM method combines a sampling of representative configurations from density functional theory (DFT)-based ab initio molecular dynamics (AIMD) simulations with DFT calculations of local-mode vibrational frequencies and couplings. Calculations of inhomogeneous OH stretch vibrational spectra for H+(H2O)4 and H+(H2O)21 agree well with experiment and higher-level calculations, and decompositions of the calculated spectra in terms of the coupled modes aids in the interpretation of the spectra. This observation is consistent with the idea that capturing anharmonicity and coupling is as important to accuracy as the underlying level of electronic structure theory. The CLM calculations can easily discern the configuration that dominates the experimental measurement for H+(H2O)5, which can adopt several low-energy conformations.

Local Electric Fields in Aqueous Electrolytes.

Vibrational Stark shifts were explored in aqueous solutions of organic molecules with carbonyl- and nitrile-containing constituents. In many cases, the vibrational resonances from these moieties shifted toward lower frequency as salt was introduced into solution. This is in contrast to the blue-shift that would be expected based upon Onsager's reaction field theory. Salts containing well-hydrated cations like Mg2+ or Li+ led to the most pronounced Stark shift for the carbonyl group, while poorly hydrated cations like Cs+ had the greatest impact on nitriles. Moreover, salts containing I- gave rise to larger Stark shifts than those containing Cl-. Molecular dynamics simulations indicated that cations and anions both accumulate around the probe in an ion- and probe-dependent manner. An electric field was generated by the ion pair, which pointed from the cation to the anion through the vibrational chromophore. This resulted from solvent-shared binding of the ions to the probes, consistent with their positions in the Hofmeister series. The "anti-Onsager" Stark shifts occur in both vibrational spectroscopy and fluorescence measurements.