Understanding the electrostatic forces and features within highly heterogeneous, anisotropic, and chemically complex enzyme active sites and their connection to biological catalysis remains a longstanding challenge, in part due to the paucity of incisive experimental probes of electrostatic properties within proteins. To quantitatively assess the landscape of electrostatic fields at discrete locations and orientations within an enzyme active site, we have incorporated site-specific thiocyanate vibrational probes into multiple positions within bacterial ketosteroid isomerase. A battery of X-ray crystallographic, vibrational Stark spectroscopy, and NMR studies revealed electrostatic field heterogeneity of 8 MV∕cm between active site probe locations and widely differing sensitivities of discrete probes to common electrostatic perturbations from mutation, ligand binding, and pH changes. Electrostatic calculations based on active site ionization states assigned by literature precedent and computational pK a prediction were unable to quantitatively account for the observed vibrational band shifts. However, electrostatic models of the D40N mutant gave qualitative agreement with the observed vibrational effects when an unusual ionization of an active site tyrosine with a pK a near 7 was included. UV-absorbance and 13 C NMR experiments confirmed the presence of a tyrosinate in the active site, in agreement with electrostatic models. This work provides the most direct measure of the heterogeneous and anisotropic nature of the electrostatic environment within an enzyme active site, and these measurements provide incisive benchmarks for further developing accurate computational models and a foundation for future tests of electrostatics in enzymatic catalysis. E xtensive structural studies of enzymes have revealed that biological catalysis occurs within sequestered active site crevices that solvate reacting substrates with an anisotropic and chemically complex constellation of charged, polar, and hydrophobic groups. But beyond visualization of the chemical composition and architecture of active sites permitted by the rich library of available protein structures, our understanding of the electrostatic nature and properties of this highly heterogeneous environment and its role in molecular recognition and catalysis is largely based on simulations. Computations using three-dimensional structures routinely provide electrostatic potentials, as in Fig. 1A. These potentials have been used to identify and characterize protein-protein and protein-ligand binding sites, and interaction energies derived from these calculations have suggested key contributions from the electrostatic environment to enzymatic rate enhancement and specificity (1-5). However, even the most sophisticated computational methods have multiple approximations and limitations (6-8) and, most importantly, have not been rigorously evaluated through cycles of nontrivial computational predictions and experimental tests. Incisive and quantitative experimental measures o...