CONSPECTUS Over the past decade, we have developed a spectroscopic approach to measure electric fields inside matter with high spatial (<1 Å) and field (<1 MV/cm) resolution. The approach hinges on exploiting a physical phenomenon known as the vibrational Stark effect (VSE), which ultimately provides a direct mapping between observed vibrational frequencies and electric fields. Therefore, the frequency of a vibrational probe encodes information about the local electric field in the vicinity around the probe. The VSE method has enabled us to understand in great detail the underlying physical nature of several important biomolecular phenomena, such as drug–receptor selectivity in tyrosine kinases, catalysis by the enzyme ketosteroid isomerase, and unidirectional electron transfer in the photosynthetic reaction center. Beyond these specific examples, the VSE has provided a conceptual foundation for how to model intermolecular (noncovalent) interactions in a quantitative, consistent, and general manner. The starting point for research in this area is to choose (or design) a vibrational probe to interrogate the particular system of interest. Vibrational probes are sometimes intrinsic to the system in question, but we have also devised ways to build them into the system (extrinsic probes), often with minimal perturbation. With modern instruments, vibrational frequencies can increasingly be recorded with very high spatial, temporal, and frequency resolution, affording electric field maps correspondingly resolved in space, time, and field magnitude. In this Account, we set out to explain the VSE in broad strokes to make its relevance accessible to chemists of all specialties. Our intention is not to provide an encyclopedic review of published work but rather to motivate the underlying framework of the methodology and to describe how we make and interpret the measurements. Using certain vibrational probes, benchmarked against computer models, it is possible to use the VSE to measure absolute electric fields in arbitrary environments. The VSE approach provides an organizing framework for thinking generally about intermolecular interactions in a quantitative way and may serve as a useful conceptual tool for molecular design.

Enzymes use protein architecture to impose specific electrostatic fields onto their bound substrates, but the magnitude and catalytic effect of these electric fields have proven difficult to quantify with standard experimental approaches. Using vibrational Stark effect spectroscopy, we found that the active site of the enzyme ketosteroid isomerase (KSI) exerts an extremely large electric field onto the C=O chemical bond that undergoes a charge rearrangement in KSI’s rate-determining step. Moreover, we found that the magnitude of the electric field exerted by the active site strongly correlates with the enzyme’s catalytic rate enhancement, enabling us to quantify the fraction of the catalytic effect that is electrostatic in origin. The measurements described here may help explain the role of electrostatics in many other enzymes and biomolecular systems.

What happens inside an enzyme’s active site to allow slow and difficult chemical reactions to occur so rapidly? This question has occupied biochemists’ attention for a long time. Computer models of increasing sophistication have predicted an important role for electrostatic interactions in enzymatic reactions, yet this hypothesis has proved vexingly difficult to test experimentally. Recent experiments utilizing the vibrational Stark effect make it possible to measure the electric field a substrate molecule experiences when bound inside its enzyme’s active site. These experiments have provided compelling evidence supporting a major electrostatic contribution to enzymatic catalysis. Here, we review these results and develop a simple model for electrostatic catalysis that enables us to incorporate disparate concepts introduced by many investigators to describe how enzymes work into a more unified framework stressing the importance of electric fields at the active site.

Vibrational probes can provide a direct read-out of the local electrostatic field in complex molecular environments, such as protein binding sites and enzyme active sites. This information provides an experimental method to explore the underlying physical causes of important biomolecular processes such as binding and catalysis. However, specific chemical interactions such as hydrogen bonds can have complicated effects on vibrational probes and confound simple electrostatic interpretations of their frequency shifts. We employ vibrational Stark spectroscopy along with infrared spectroscopy of carbonyl probes in different solvent environments and in Ribonuclease S to understand the sensitivity of carbonyl frequencies to electrostatic fields, including those due to hydrogen bonds. Additionally, we carried out molecular dynamics simulations to calculate ensemble-averaged electric fields in solvents and in Ribonuclease S, and found excellent correlation between calculated fields and vibrational frequencies. These data enabled the construction of a robust field-frequency calibration curve for the C=O vibration. The present results suggest that carbonyl probes are capable of quantitatively assessing the electrostatics of hydrogen bonding, making them promising for future study of protein function.

Electrostatic interactions provide a primary connection between a protein’s three-dimensional structure and its function. Infrared (IR) probes are useful because vibrational frequencies of certain chemical groups, such as nitriles, are linearly sensitive to local electrostatic field, and can serve as a molecular electric field meter. IR spectroscopy has been used to study electrostatic changes or fluctuations in proteins, but measured peak frequencies have not been previously mapped to total electric fields, because of the absence of a field-frequency calibration and the complication of local chemical effects such as H-bonds. We report a solvatochromic model that provides a means to assess the H-bonding status of aromatic nitrile vibrational probes, and calibrates their vibrational frequencies to electrostatic field. The analysis involves correlations between the nitrile’s IR frequency and its 13C chemical shift, whose observation is facilitated by a robust method for introducing isotopes into aromatic nitriles. The method is tested on the model protein Ribonuclease S (RNase S) containing a labeled p-CN-Phe near the active site. Comparison of the measurements in RNase S against solvatochromic data gives an estimate of the average total electrostatic field at this location. The value determined agrees quantitatively with MD simulations, suggesting broader potential for the use of IR probes in the study of protein electrostatics.

The physical properties of solvents strongly affect the spectra of dissolved solutes, and this phenomenon can be exploited to gain insight into the solvent-solute interaction. The large solvatochromic shifts observed for many dye molecules in polar solvents are due to variations in the solvent reaction field, and these shifts are widely used to estimate the change in the dye’s dipole moment upon photoexcitation, which is typically on the order of ~1-10 Debye. In contrast, the change in dipole moment for vibrational transitions is approximately two orders of magnitude smaller. Nonetheless, vibrational chromophores display significant solvatochromism, and the relative contributions of specific chemical interactions and electrostatic interactions are debated, complicating the interpretation of vibrational frequency shifts in complex systems such as proteins. Here we present a series of substituted benzonitriles that display widely varying degrees of vibrational solvatochromism. In most cases, this variation can be quantitatively described by the experimentally determined Stark tuning rate, coupled with a simple Onsager-like model of solvation, reinforcing the view that vibrational frequency shifts are largely caused by electrostatic interactions. In addition, we discuss specific cases where continuum solvation models fail to predict solvatochromic shifts, revealing the necessity for more advanced theoretical models that capture local aspects of solute-solvent interactions.

The electric field created by a condensed phase environment is a powerful and convenient descriptor for intermolecular interactions. Not only does it provide a unifying language to compare many different types of interactions, but it also possesses clear connections to experimental observables, such as vibrational Stark effects. We calculate here the electric fields experienced by a vibrational chromophore (the carbonyl group of acetophenone) in an array of solvents of diverse polarities using molecular dynamics simulations with the AMOEBA polarizable force field. The mean and variance of the calculated electric fields correlate well with solvent-induced frequency shifts and band broadening, suggesting Stark effects as the underlying mechanism of these key solution phase spectral effects. Compared to fixed-charge and continuum models, AMOEBA was the only model examined that could describe non-polar, polar, and hydrogen bonding environments in a consistent fashion. Nevertheless, we found that fixed-charge force fields and continuum models were able to replicate some results of the polarizable simulations accurately, allowing us to clearly identify which properties and situations require explicit polarization and/or atomistic representations to be modeled properly, and for which properties and situations simpler models are sufficient. We also discuss the ramifications of these results for modeling electrostatics in complex environments, such as proteins.

Enzymes use protein architectures to create highly specialized structural motifs that can greatly enhance the rates of complex chemical transformations. Here, we use experiments, combined with ab initio simulations that exactly include nuclear quantum effects, to show that a triad of strongly hydrogen-bonded tyrosine residues within the active site of the enzyme ketosteroid isomerase (KSI) facilitates quantum proton delocalization. This delocalization dramatically stabilizes the deprotonation of an active-site tyrosine residue, resulting in a very large isotope effect on its acidity. When an intermediate analog is docked, it is incorporated into the hydrogen-bond network, giving rise to extended quantum proton delocalization in the active site. These results shed light on the role of nuclear quantum effects in the hydrogen-bond network that stabilizes the reactive intermediate of KSI, and the behavior of protons in biological systems containing strong hydrogen bonds.enzyme | hydrogen bonding | nuclear quantum effects | proton delocalization | ab initio path integral molecular dynamics A lthough many biological processes can be well-described with classical mechanics, there has been much interest and debate as to the role of quantum effects in biological systems ranging from photosynthetic energy transfer, to photoinduced isomerization in the vision cycle and avian magnetoreception (1). For example, nuclear quantum effects, such as tunneling and zero-point energy (ZPE), have been observed to lead to kinetic isotope effects of greater than 100 in biological proton and proton-coupled electron transfer processes (2, 3). However, the role of nuclear quantum effects in determining the groundstate thermodynamic properties of biological systems, which manifest as equilibrium isotope effects, has gained significantly less attention (4).Ketosteroid isomerase (KSI) possesses one of the highest enzyme unimolecular rate constants and thus, is considered a paradigm of proton transfer catalysis in enzymology (5-11). The remarkable rate of KSI is intimately connected to the formation of a hydrogen-bond network in its active site (Fig. 1A), which acts to stabilize a charged dienolate intermediate, lowering its free energy by ∼11 kcal/mol (1 kcal = 4.18 kJ) relative to solution (Fig. S1) (6). This extended hydrogen-bond network in the active site links the substrate to Asp103 and Tyr16, with the latter further hydrogen-bonded to Tyr57 and Tyr32, which is shown in Fig. 1A.The mutant KSI D40N preserves the structure of the wild-type (WT) enzyme while mimicking the protonation state of residue 40 in the intermediate complex (Fig. 1B), therefore permitting experimental investigation of an intermediate-like state of the enzyme (6,(12)(13)(14). Experiments have identified that, in the absence of an inhibitor, one of the residues in the active site of KSI D40N is deprotonated (12). Although one might expect the carboxylic acid of Asp103 to be deprotonated, the combination of recent 13 C NMR and ultraviolet visible spectroscopy (UVVis) e...