Decomposition of vibrational shifts of nitriles into electrostatic and hydrogen-bonding effects

Unnatural Amino Acids for Biological Spectroscopy and Microscopy

Due to advances in methods for site-specific incorporation of unnatural amino acids (UAAs) into proteins, a large number of UAAs with tailored chemical and/or physical properties have been developed and used in a wide array of biological applications. In particular, UAAs with specific spectroscopic characteristics can be used as external reporters to produce additional signals, hence increasing the information content obtainable in protein spectroscopic and/or imaging measurements. In this Review, we summarize the progress in the past two decades in the development of such UAAs and their applications in biological spectroscopy and microscopy, with a focus on UAAs that can be used as site-specific vibrational, fluorescence, electron paramagnetic resonance (EPR), or nuclear magnetic resonance (NMR) probes. Wherever applicable, we also discuss future directions.

Critical Evaluation of Polarizable and Nonpolarizable Force Fields for Proteins Using Experimentally Derived Nitrile Electric Fields

Molecular dynamics (MD) simulations are frequently carried out for proteins to investigate the role of electrostatics in their biological function. The choice of force field (FF) can significantly alter the MD results, as the simulated local electrostatic interactions lack benchmarking in the absence of appropriate experimental methods. We recently reported that the transition dipole moment (TDM) of the popular nitrile vibrational probe varies linearly with the environmental electric field, overcoming well-known hydrogen bonding (H-bonding) issues for the nitrile frequency and, thus, enabling the unambiguous measurement of electric fields in proteins (J. Am. Chem. Soc. 2022, 144 (17), 7562-7567). Herein, we utilize this new strategy to enable comparisons of experimental and simulated electric fields in protein environments. Specifically, previously determined TDM electric fields exerted onto nitrile-containing o-cyanophenylalanine residues in photoactive yellow protein are compared with MD electric fields from the fixed-charge AMBER FF and the polarizable AMOEBA FF. We observe that the electric field distributions for H-bonding nitriles are substantially affected by the choice of FF. As such, AMBER underestimates electric fields for nitriles experiencing moderate field strengths; in contrast, AMOEBA robustly recapitulates the TDM electric fields. The FF dependence of the electric fields can be partly explained by the presence of additional negative charge density along the nitrile bond axis in AMOEBA, which is due to the inclusion of higher-order multipole parameters; this, in turn, begets more head-on nitrile H-bonds. We conclude by discussing the implications of the FF dependence for the simulation of nitriles and proteins in general.

The hydrogen-bonding dynamics of water to a nitrile-functionalized electrode is modulated by voltage according to ultrafast 2D IR spectroscopy

We report the hydrogen-bonding dynamics of water to a nitrile-functionalized and plasmonic electrode surface as a function of applied voltage. The surface-enhanced two-dimensional infrared spectra exhibit hydrogen-bonded and non-hydrogen-bonded nitrile features in similar proportions, plus cross peaks between the two. Isotopic dilution experiments show that the cross peaks arise predominantly from chemical exchange between hydrogen-bonded and non-hydrogen-bonded nitriles. The chemical exchange rate depends upon voltage, with the hydrogen bond of the water to the nitriles breaking 2 to 3 times slower (>63 vs. 25 ps) under a positive as compared to a negative potential. Spectral diffusion created by hydrogen-bond fluctuations occurs on a ~1 ps timescale and is moderately potential-dependent. Timescales from molecular dynamics simulations agree qualitatively with the experiment and show that a negative voltage causes a small net displacement of water away from the surface. These results show that the voltage applied to an electrode can alter the timescales of solvent motion at its interface, which has implications for electrochemically driven reactions.

Good Vibrations Report on the DNA Quadruplex Binding of an Excited State Amplified Ruthenium Polypyridyl IR Probe

The nitrile containing Ru(II)polypyridyl complex [Ru(phen)2(11,12-dCN-dppz)]2+ (1) is shown to act as a sensitive infrared probe of G-quadruplex (G4) structures. UV-visible absorption spectroscopy reveals enantiomer sensitive binding for the hybrid htel(K) and antiparallel htel(Na) G4s formed by the human telomer sequence d[AG3(TTAG3)3]. Time-resolved infrared (TRIR) of 1 upon 400 nm excitation indicates dominant interactions with the guanine bases in the case of Λ-1/htel(K), Δ-1/htel(K), and Λ-1/htel(Na) binding, whereas Δ-1/htel(Na) binding is associated with interactions with thymine and adenine bases in the loop. The intense nitrile transient at 2232 cm-1 undergoes a linear shift to lower frequency as the solution hydrogen bonding environment decreases in DMSO/water mixtures. This shift is used as a sensitive reporter of the nitrile environment within the binding pocket. The lifetime of 1 in D2O (ca. 100 ps) is found to increase upon DNA binding, and monitoring of the nitrile and ligand transients as well as the diagnostic DNA bleach bands shows that this increase is related to greater protection from the solvent environment. Molecular dynamics simulations together with binding energy calculations identify the most favorable binding site for each system, which are in excellent agreement with the observed TRIR solution study. This study shows the power of combining the environmental sensitivity of an infrared (IR) probe in its excited state with the TRIR DNA "site effect" to gain important information about the binding site of photoactive agents and points to the potential of such amplified IR probes as sensitive reporters of biological environments.

Vibrational Spectroscopy of Phytochromes

Phytochromes are biological photoswitches that translate light into physiological functions. Spectroscopic techniques are essential tools for molecular research into these photoreceptors. This review is directed at summarizing how resonance Raman and IR spectroscopy contributed to an understanding of the structure, dynamics, and reaction mechanism of phytochromes, outlining the substantial experimental and theoretical challenges and describing the strategies to master them. It is shown that the potential of the various vibrational spectroscopic techniques can be most efficiently exploited using integral approaches via a combination of theoretical methods as well as other experimental techniques.

Frequency Changes in Terminal Alkynes Provide Strong, Sensitive, and Solvatochromic Raman Probes of Biochemical Environments

The C≡C stretching frequencies of terminal alkynes appear in the "clear" window of vibrational spectra, so they are attractive and increasingly popular as site-specific probes in complicated biological systems like proteins, cells, and tissues. In this work, we collected infrared (IR) absorption and Raman scattering spectra of model compounds, artificial amino acids, and model proteins that contain terminal alkyne groups, and we used our results to draw conclusions about the signal strength and sensitivity to the local environment of both aliphatic and aromatic terminal alkyne C≡C stretching bands. While the IR bands of alkynyl model compounds displayed surprisingly broad solvatochromism, their absorptions were weak enough that alkynes can be ruled out as effective IR probes. The same solvatochromism was observed in model compounds' Raman spectra, and comparisons to published empirical solvent scales (including a linear regression against four meta-aggregated solvent parameters) suggested that the alkyne C≡C stretching frequency mainly reports on local electronic interactions (i.e., short-range electron donor-acceptor interactions) with solvent molecules and neighboring functional groups. The strong solvatochromism observed here for alkyne stretching bands introduces an important consideration for Raman imaging studies based on these signals. Raman signals for alkynes (especially those that are π-conjugated) can be exceptionally strong and should permit alkynyl Raman signals to function as probes at very low concentrations, as compared to other widely used vibrational probe groups like azides and nitriles. We incorporated homopropargyl glycine into a transmembrane helical peptide via peptide synthesis, and we installed p-ethynylphenylalanine into the interior of the Escherichia coli fatty acid acyl carrier protein using a genetic code expansion technique. The Raman spectra from each of these test systems indicate that alkynyl C≡C bands can act as effective and unique probes of their local biomolecular environments. We provide guidance for the best possible future uses of alkynes as solvatochromic Raman probes, and while empirical explanations of the alkyne solvatochromism are offered, open questions about its physical basis are enunciated.

Potential Distribution across Model Membranes

Membrane models assembled on electrodes are widely used tools to study potential-dependent molecular processes at or in membranes. However, the relationship between the electrode potential and the potential across the membrane is not known. Here we studied lipid bilayers immobilized on mixed self-assembled monolayers (SAM) on Au electrodes. The mixed SAM was composed of thiol derivatives of different chain lengths such that between the islands of the short one, mercaptobenzonitrile (MBN), and the tethered lipid bilayer an aqueous compartment was formed. The nitrile function of MBN, which served as a reporter group for the vibrational Stark effect (VSE), was probed by surface-enhanced infrared absorption spectroscopy to determine the local electric field as a function of the electrode potential for pure MBN, mixed SAM, and the bilayer system. In parallel, we calculated electric fields at the VSE probe by molecular dynamics (MD) simulations for different charge densities on the metal, thereby mimicking electrode potential changes. The agreement with the experiments was very good for the calculations of the pure MBN SAM and only slightly worse for the mixed SAM. The comparison with the experiments also guided the design of the bilayer system in the MD setups, which were selected to calculate the electrode potential dependence of the transmembrane potential, a quantity that is not directly accessible by the experiments. The results agree very well with estimates in previous studies and thus demonstrate that the present combined experimental-theoretical approach is a promising tool for describing potential-dependent processes at biomimetic interfaces.

Origin of thiocyanate spectral shifts in water and organic solvents

Vibrational spectroscopy is a useful technique for probing chemical environments. The development of models that can reproduce the spectra of nitriles and azides is valuable because these probes are uniquely suited for investigating complex systems. Empirical vibrational spectroscopic maps are commonly employed to obtain the instantaneous vibrational frequencies during molecular dynamics simulations but often fail to adequately describe the behavior of these probes, especially in its transferability to a diverse range of environments. In this paper, we demonstrate several reasons for the difficulty in constructing a general-purpose vibrational map for methyl thiocyanate (MeSCN), a model for cyanylated biological probes. In particular, we found that electrostatics alone are not a sufficient metric to categorize the environments of different solvents, and the dominant features in intermolecular interactions in the energy landscape vary from solvent to solvent. Consequently, common vibrational mapping schemes do not cover all essential interaction terms adequately, especially in the treatment of van der Waals interactions. Quantum vibrational perturbation (QVP) theory, along with a combined quantum mechanical and molecular mechanical potential for solute-solvent interactions, is an alternative and efficient modeling technique, which is compared in this paper, to yield spectroscopic results in good agreement with experimental FTIR. QVP has been used to analyze the computational data, revealing the shortcomings of the vibrational maps for MeSCN in different solvents. The results indicate that insights from QVP analysis can be used to enhance the transferability of vibrational maps in future studies.

Characterization of Acetonitrile Isotopologues as Vibrational Probes of Electrolytes

Acetonitrile has emerged as a solvent candidate for novel electrolyte formulations in metal-ion batteries and supercapacitors. It features a bright local C≡N stretch vibrational mode whose infrared (IR) signature is sensitive to battery-relevant cations (Li⁺, Mg²⁺, Zn²⁺, Ca²⁺) both in pure form and in the presence of water admixture across a full possible range of concentrations from the dilute to the superconcentrated regime. Stationary and time-resolved IR spectroscopy thus emerges as a natural tool to study site-specific intermolecular interactions from the solvent perspective without introducing an extrinsic probe that perturbs solution morphology and may not represent the intrinsic dynamics in these electrolytes. The metal-coordinated acetonitrile, water-separated metal-acetonitrile pair, and free solvent each have a distinct vibrational signature that allows their unambiguous differentiation. The IR band frequency of the metal-coordinated acetonitrile depends on the ion charge density. To study the ion transport dynamics, it is necessary to differentiate energy-transfer processes from structural interconversions in these electrolytes. Isotope labeling the solvent is a necessary prerequisite to separate these processes. We discuss the design principles and choice of the CD₃¹³CN label and characterize its vibrational spectroscopy in these electrolytes. The Fermi resonance between ¹³C≡N and C-D stretches complicates the spectral response but does not prevent its effective utilization. Time-resolved two-dimensional (2D) IR spectroscopy can be performed on a mixture of acetonitrile isotopologues and much can be learned about the structural dynamics of various species in these formulations.

Association-Dissociation Dynamics of Ionic Electrolytes in Low Dielectric Medium

Ionic electrolytes are known to form various complexes which exist in dynamic equilibrium in a low dielectric medium. However, structural characterization of these complexes has always posed a great challenge to the scientific community. An additional challenge is the estimation of the dynamic association-dissociation time scales (lifetime of the complexes), which are key to the fundamental understanding of ion transport. In this work, we have used a combination of infrared absorption spectroscopy, two-dimensional infrared spectroscopy, molecular dynamics simulations, and density functional theory calculations to characterize the various ion complexes formed by the thiocyanate-based ionic electrolytes as a function of different cations in a low dielectric medium. Our results demonstrate that thiocyanate is an excellent vibrational reporter of the heterogeneous ion complexes undergoing association-dissociation dynamics. We find that the ionic electrolytes exist as contact ion pairs, dimers, and clusters in a low dielectric medium. The relative ratios of the various ion complexes are sensitive to the cations. In addition to the interactions between the thiocyanate anion and the countercation, the solute-solvent interactions drive the dynamic equilibrium. We have estimated the association-dissociation dynamics time scales from two-dimensional infrared spectroscopy. The exchange time scale involving the cluster is faster than that between a dimer and an ion pair. Moreover, we find that the dynamic equilibrium between the cluster and another ion complex is correlated to the solvent fluctuations.

Local Electric Field Changes during the Photoconversion of the Bathy Phytochrome Agp2

Phytochromes switch between a physiologically inactive and active state via a light-induced reaction cascade, which is initiated by isomerization of the tetrapyrrole chromophore and leads to the functionally relevant secondary structure transition of a protein segment (tongue). Although details of the underlying cause-effect relationships are not known, electrostatic fields are likely to play a crucial role in coupling chromophores and protein structural changes. Here, we studied local electric field changes during the photoconversion of the dark state Pfr to the photoactivated state Pr of the bathy phytochrome Agp2. Substituting Tyr165 and Phe192 in the chromophore pocket by para-cyanophenylalanine (pCNF), we monitored the respective nitrile stretching modes in the various states of photoconversion (vibrational Stark effect). Resonance Raman and IR spectroscopic analyses revealed that both pCNF-substituted variants undergo the same photoinduced structural changes as wild-type Agp2. Based on a structural model for the Pfr state of F192pCNF, a molecular mechanical-quantum mechanical approach was employed to calculate the electric field at the nitrile group and the respective stretching frequency, in excellent agreement with the experiment. These calculations serve as a reference for determining the electric field changes in the photoinduced states of F192pCNF. Unlike F192pCNF, the nitrile group in Y165pCNF is strongly hydrogen bonded such that the theoretical approach is not applicable. However, in both variants, the largest changes of the nitrile stretching modes occur in the last step of the photoconversion, supporting the view that the proton-coupled restructuring of the tongue is accompanied by a change of the electric field.

Mapping solvation heterogeneity in live cells by hyperspectral stimulated Raman scattering microscopy

Water provides a dynamic matrix in which all biochemical processes occur in living organisms. The structure and dynamics of intracellular water constitute the cornerstone for understanding all aspects of cellular function. Fundamentally, direct visualization of subcellular solvation heterogeneity is essential but remains challenging with commonly used nuclear magnetic resonance methods due to poor spatial resolution. To explore this question, we demonstrate a vibrational-shift imaging approach by combining the spectral-focusing hyperspectral stimulated Raman scattering technique with an environmentally sensitive nitrile probe. The sensing ability of a near-infrared nitrile-containing molecule is validated in the solution phase, microscopic droplets, and cellular environments. Finally, we quantitatively measure the subcellular solvation variance between the cytoplasm (29.5%, S.E. 1.8%) and the nucleus (57.3%, S.E. 1.0%), which is in good agreement with previous studies. This work sheds light on heterogeneous solvation in live systems using coherent Raman microscopy and opens up new avenues to explore environmental variance in complex systems with high spatiotemporal resolution.

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.

Electric Field Measurements Reveal the Pivotal Role of Cofactor-Substrate Interaction in Dihydrofolate Reductase Catalysis

The contribution of ligand–ligand electrostatic interaction to transition state formation during enzyme catalysis has remained unexplored, even though electrostatic forces are known to play a major role in protein functions and have been investigated by the vibrational Stark effect (VSE). To monitor electrostatic changes along important steps during catalysis, we used a nitrile probe (T46C-CN) inserted proximal to the reaction center of three dihydrofolate reductases (DHFRs) with different biophysical properties, Escherichia coli DHFR (EcDHFR), its conformationally impaired variant (EcDHFR-S148P), and Geobacillus stearothermophilus DHFR (BsDHFR). Our combined experimental and computational approach revealed that the electric field projected by the substrate toward the probe negates those exerted by the cofactor when both are bound within the enzymes. This indicates that compared to previous models that focus exclusively on subdomain reorganization and protein–ligand contacts, ligand–ligand interactions are the key driving force to generate electrostatic environments conducive for catalysis.

Biosynthetic Incorporation of Site-Specific Isotopes in β-Lactam Antibiotics Enables Biophysical Studies

A biophysical understanding of the mechanistic, chemical, and physical origins underlying antibiotic action and resistance is vital to the discovery of novel therapeutics and the development of strategies to combat the growing emergence of antibiotic resistance. The site-specific introduction of stable-isotope labels into chemically complex natural products is particularly important for techniques such as NMR, IR, mass spectrometry, imaging, and kinetic isotope effects. Toward this goal, we developed a biosynthetic strategy for the site-specific incorporation of ¹³C labels into the canonical β-lactam carbonyl of penicillin G and cefotaxime, the latter via cephalosporin C. This was achieved through sulfur-replacement with 1-¹³C-l-cysteine, resulting in high isotope incorporations and milligram-scale yields. Using ¹³C NMR and isotope-edited IR difference spectroscopy, we illustrate how these molecules can be used to interrogate interactions with their protein targets, e.g., TEM-1 β-lactamase. This method provides a feasible route to isotopically labeled penicillin and cephalosporin precursors for future biophysical studies.

Differential evolution algorithm approach for describing vibrational solvatochromism

We model the solvation-induced vibrational frequency shifts of the amide I and amide II modes of N-methylacetamide in water and the nitrile stretch mode of acetonitrile in water by expressing the frequency shift as a polynomial function expanded by the inverse power of interatomic distances. The coefficients of the polynomial are optimized to minimize the deviation between the predicted frequency shifts and those calculated with quantum chemistry methods. Here, we show that a differential evolution algorithm combined with singular value decomposition is useful to find the optimum set of coefficients of polynomial terms. The differential evolution optimization shows that only a few terms in the polynomial are dominant in the contribution to the vibrational frequency shifts. We anticipate that the present work paves the way for further developing different genetic algorithms and machine learning schemes for their applications to vibrational spectroscopic studies.

Heterogeneous and Highly Dynamic Interface in Plastocyanin-Cytochrome f Complex Revealed by Site-Specific 2D-IR Spectroscopy

Transient protein complexes are crucial for sustaining dynamic cellular processes. The complexes of electron-transfer proteins are a notable example, such as those formed by plastocyanin (Pc) and cytochrome f (cyt f) in the photosynthetic apparatus. The dynamic and heterogeneous nature of these complexes, however, makes their study challenging. To better elucidate the complex of Nostoc Pc and cyt f, 2D-IR spectroscopy coupled to site-specific labeling with cyanophenylalanine infrared (IR) probes was employed to characterize how the local environments at sites along the surface of Pc were impacted by cyt f binding. The results indicate that Pc most substantially engages with cyt f via the hydrophobic patch around the copper redox site. Complexation with cyt f led to an increase in inhomogeneous broadening of the probe absorptions, reflective of increased heterogeneity of interactions with their environment. Notably, most of the underlying states interconverted very rapidly (1 to 2 ps), suggesting a complex with a highly mobile interface. The data support a model of the complex consisting of a large population of an encounter complex. Additionally, the study demonstrates the application of 2D-IR spectroscopy with site-specifically introduced probes to reveal new quantitative insight about dynamic biochemical systems.

Site-specific 2D IR spectroscopy: a general approach for the characterization of protein dynamics with high spatial and temporal resolution

The conformational heterogeneity and dynamics of protein side chains contribute to function, but investigating exactly how is hindered by experimental challenges arising from the fast timescales involved and the spatial heterogeneity of protein structures. The potential of two-dimensional infrared (2D IR) spectroscopy for measuring conformational heterogeneity and dynamics with unprecedented spatial and temporal resolution has motivated extensive effort to develop amino acids with functional groups that have frequency-resolved absorptions to serve as probes of their protein microenvironments. We demonstrate the full advantage of the approach by selective incorporation of the probe p-cyanophenylalanine at six distinct sites in a Src homology 3 domain and the application of 2D IR spectroscopy to site-specifically characterize heterogeneity and dynamics and their contribution to cognate ligand binding. The approach revealed a wide range of microenvironments and distinct responses to ligand binding, including at the three adjacent, conserved aromatic residues that form the recognition surface of the protein. Molecular dynamics simulations performed for all the labeled proteins provide insight into the underlying heterogeneity and dynamics. Similar application of 2D IR spectroscopy and site-selective probe incorporation will allow for the characterization of heterogeneity and dynamics of other proteins, how heterogeneity and dynamics are affected by solvation and local structure, and how they might contribute to biological function.

New Insights into Quinine-DNA Binding Using Raman Spectroscopy and Molecular Dynamics Simulations

Quinine's ability to bind DNA and potentially inhibit transcription and translation has been examined as a mode of action for its antimalarial activity. UV absorption and fluorescence-based studies have lacked the chemical specificity to develop an unambiguous molecular-level picture of the binding interaction. To address this, we use Raman spectroscopy and molecular dynamics (MD) to investigate quinine-DNA interactions. We demonstrate that quinine's strongest Raman band in the fingerprint region, which derives from a symmetric stretching mode of the quinoline ring, is highly sensitive to the local chemical environment and pH. The frequency shifts observed for this mode in solvents of varying polarity can be explained in terms of the Stark effect using a simple Onsager solvation model, indicating that the vibration reports on the local electrostatic environment. However, specific chemical interactions between the quinoline ring and its environment, such as hydrogen bonding and π-stacking, perturb the frequency of this mode in a more complicated but predictable manner. We use this vibration as a spectroscopic probe to investigate the binding interaction between quinine and DNA. We find that, when the quinoline ring is protonated, quinine weakly intercalates into DNA by forming π-stacking interactions with the base pairs. The Raman spectra indicate that quinine can intercalate into DNA with a ratio reaching up to roughly one molecule per 25 base pairs. Our results are confirmed by MD simulations, which also show that the quinoline ring adopts a t-shaped π-stacking geometry with the DNA base pairs, whereas the quinuclidine head group weakly interacts with the phosphate backbone in the minor groove. We expect that the spectral correlations determined here will enable future studies to probe quinine's antimalarial activities, such as disrupting hemozoin biocrystallization, which is hypothesized to be, among other things, one of its primary modes of action against Plasmodium parasites.

Quantitative Measurement of Intrinsic GTP Hydrolysis for Carcinogenic Glutamine 61 Mutants in H-Ras

Mutations of human oncoprotein p²¹H-Ras (hereafter "Ras") at glutamine 61 are known to slow the rate of guanosine triphosphate (GTP) hydrolysis and transform healthy cells into malignant cells. It has been hypothesized that this glutamine plays a role in the intrinsic mechanism of GTP hydrolysis by interacting with an active site water molecule that stabilizes the formation of the charged transition state at the γ-phosphate during hydrolysis. However, there is no comprehensive data set of the effects of mutations to Q61 on the protein's intrinsic catalytic rate, structure, or interactions with water at the active site. Here, we present the first comprehensive and quantitative set of initial rates of intrinsic hydrolysis for all stable variants of RasQ61X. We further conducted enhanced molecular dynamics (MD) simulations of each construct to determine the solvent accessible surface area (SASA) of the side chain at position 61 and compared these results to previously measured changes in electric fields caused by RasQ61X mutations. For polar and negatively charged residues, we found that the rates are normally distributed about an optimal electrostatic contribution, close to that of the native Q61 residue, and the rates are strongly correlated to the number of waters in the active site. Together, these results support a mechanism of GTP hydrolysis in which Q61 stabilizes a transient hydronium ion, which then stabilizes the transition state while the γ-phosphate is undergoing nucleophilic attack by a second, catalytically active water molecule. We discuss the implications of such a mechanism on future strategies for combating Ras-based cancers.

Long-Range Modulations of Electric Fields in Proteins

Electrostatic interactions are essential for controlling the protein structure and function. Whereas so far experimental and theoretical efforts focused on the effect of local electrostatics, this work aims at elucidating the long-range modulation of electric fields in proteins upon binding to charged surfaces. The study is based on cytochrome c (Cytc) variants carrying nitrile reporters for the vibrational Stark effect that are incorporated into the protein via genetic engineering and chemical modification. The Cytc variants were thoroughly characterized with respect to possible structural perturbations due to labeling. For the proteins in solution, the relative hydrogen bond occupancy and the calculated electric fields, both obtained from molecular dynamics (MD) simulations, and the experimental nitrile stretching frequencies were used to develop a relationship for separating hydrogen-bonding and non-hydrogen-bonding electric field effects. This relationship provides an excellent description for the stable Cytc variants in solution. For the proteins bound to Au electrodes coated with charged self-assembled monolayers (SAMs), the underlying MD simulations can only account for the electric field changes Δ E_ads due to the formation of the electrostatic SAM-Cytc complexes but not for the additional contribution, Δ E_int, representing the consequences of the potential drops over the electrode/SAM/protein interfaces. Both Δ E_ads and Δ E_int, determined at distances between 20 and 30 Å with respect to the SAM surface, are comparable in magnitude to the non-hydrogen-bonding electric field in the unbound protein. This long-range modulation of the internal electric field may be of functional relevance for proteins in complexes with partner proteins (Δ E_ads) and attached to membranes (Δ E_ads + Δ E_int).

Synthesis of 5-Cyano-Tryptophan as a Two-Dimensional Infrared Spectroscopic Reporter of Structure

A concise synthesis of protected 5-cyano-l-tryptophan (Trp5CN ) has been developed for 2D IR spectroscopic investigations within either peptides or proteins. To assess the potential of differently substituted cyano-tryptophans, several model cyano-indole systems were characterized using IR spectroscopy. Upon assessment of their spectroscopic properties, Trp5CN was integrated into a model peptide sequence, Trp5CN -Gly-Phe4CN , to elucidate its structure. This peptide demonstrates the capability of this probe to capture structural information by 2D IR spectroscopy. The 2D IR spectrum of the peptide in water was simulated to reveal a unique spectral signature resulting from the presence of dipolar coupling. The coupling strength between cyano labels was determined to be 1.4 cm-1 by matching the slopes along the max contour for the simulated and experimental spectrum. Using transition dipole coupling, a distance between the two probes of 13 Å was calculated.

Cyanylated Cysteine Reports Site-Specific Changes at Protein-Protein-Binding Interfaces Without Perturbation

To investigate the cyanylated cysteine vibrational probe group’s ability to report on binding-induced changes along a protein–protein interface, the probe group was incorporated at several sites in a peptide of the calmodulin (CaM)-binding domain of skeletal muscle myosin light chain kinase. Isothermal titration calorimetry was used to determine the binding thermodynamics between calmodulin and each peptide. For all probe positions, the binding affinity was nearly identical to that of the unlabeled peptide. The CN stretching infrared band was collected for each peptide free in solution and bound to calmodulin. Binding-induced shifts in the IR spectral frequencies were correlated with estimated solvent accessibility based on molecular dynamics simulations. This work generally suggests (1) that site-specific incorporation of this vibrational probe group does not cause major perturbations to its local structural environment and (2) that this small probe group might be used quite broadly to map dynamic protein-binding interfaces. However, site-specific perturbations due to artificial labeling groups can be somewhat unpredictable and should be evaluated on a site-by-site basis through complementary measurements. A fully quantitative, simulation-based interpretation of the rich probe IR spectra is still needed but appears to be possible given recent advances in simulation techniques.

A Direct, Quantitative Connection between Molecular Dynamics Simulations and Vibrational Probe Line Shapes

A quantitative connection between molecular dynamics simulations and vibrational spectroscopy of probe-labeled systems would enable direct translation of experimental data into structural and dynamical information. To constitute this connection, all-atom molecular dynamics (MD) simulations were performed for two SCN probe sites (solvent-exposed and buried) in a calmodulin-target peptide complex. Two frequency calculation approaches with substantial nonelectrostatic components, a quantum mechanics/molecular mechanics (QM/MM)-based technique and a solvatochromic fragment potential (SolEFP) approach, were used to simulate the infrared probe line shapes. While QM/MM results disagreed with experiment, SolEFP results matched experimental frequencies and line shapes and revealed the physical and dynamic bases for the observed spectroscopic behavior. The main determinant of the CN probe frequency is the exchange repulsion between the probe and its local structural neighbors, and there is a clear dynamic explanation for the relatively broad probe line shape observed at the "buried" probe site. This methodology should be widely applicable to vibrational probes in many environments.

Exploring local solvation environments of a heme protein using the spectroscopic reporter 4-cyano-l-phenylalanine

The vibrational reporter unnatural amino acid (UAA) 4-cyano-l-phenylalanine (pCNF) was genetically incorporated individually at three sites (5, 36, and 78) in the heme protein Caldanaerobacter subterraneus H-NOX to probe local hydration environments. The UAA pCNF was incorporated site-specifically using an engineered, orthogonal tRNA synthetase in E. coli. The ability of all of the pCNF-containing H-NOX proteins to form the ferrous CO, NO, or O₂ ligated and unligated states was confirmed with UV-Vis spectroscopy. The solvation state at each site of the three sites of pCNF incorporation was assessed using temperature-dependent infrared spectroscopy. Specifically, the frequency-temperature line slope (FTLS) method was utilized to show that the nitrile group at site 36 was fully solvated and the nitrile group at site 78 was de-solvated (buried) in the heme pocket. The nitrile group at site 5 was found to be partially solvated suggesting that the nitrile group was involved in moderate strength hydrogen bonds. These results were confirmed by the determination of the X-ray crystal structure of the H-NOX protein construct containing pCNF at site 5.

Fermi resonance as a means to determine the hydrogen-bonding status of two infrared probes

The C[double bond, length as m-dash]O/C[triple bond, length as m-dash]N stretching vibration arising from a carbonyl/nitrile functional group in various molecular systems has been frequently used to assess, for example, local hydrogen-bonding interactions, among other applications. However, in practice it is not always easy to ascertain whether the carbonyl or nitrile group in question is engaged in such interactions. Herein, we use 4-cyanoindole and cyclopentanone as models to show that, when a fundamental C[double bond, length as m-dash]O or C[triple bond, length as m-dash]N stretching mode is involved in Fermi resonance, the underlying vibrational coupling constant (W) is a convenient reporter of the hydrogen-bonding status of the corresponding carbonyl or nitrile group. Specifically, we find that for both groups a W value of 7.7 cm^-1 or greater is indicative of their involvement in hydrogen-bonding interactions. Furthermore, we find that, as observed in similar studies, the Fermi resonance coupling leads to quantum beats in the two-dimensional infrared spectra of 4-cyanoindole in isopropanol, with a period of about 1.9 ps.

A structural model of the hierarchical assembly of an amyloid nanosheet by an infrared probe technique

A hierarchical structural model of an amyloid nanosheet by IR probe technique.

Strategy for Modeling the Electrostatic Responses of the Spectroscopic Properties of Proteins

For better understanding and more efficient use of the spectroscopic probes (vibrational and NMR) of the local electrostatic situations inside proteins, appropriate modeling of the properties of those probes is essential. The present study is devoted to examining the strategy for constructing such models. A more well-founded derivation than the ones in previous studies is given in constructing the models. Theoretical analyses are conducted on two representative example cases related to proteins, i.e., the peptide group of the main chains and the CO and NO ligands to the Fe²⁺ ion of heme, with careful treatment of the behavior of electrons in the electrostatic responses and with verification of consistency with observable quantities. It is shown that, for the stretching frequencies and NMR chemical shifts, it is possible to construct reasonable electrostatic interaction models that encompass the situations of hydration and uniform electric field environment and thus are applicable also to the cases of nonuniform electrostatic situations, which are highly expected for inside of proteins.

Electric Fields and Enzyme Catalysis

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.

Conformational Heterogeneity and the Affinity of Substrate Molecular Recognition by Cytochrome P450cam

The broad and variable substrate specificity of cytochrome P450 enzymes makes them a model system for studying the determinants of protein molecular recognition. The archetypal cytochrome P450cam (P450cam) is a relatively specific P450, a feature once attributed to the high rigidity of its active site. However, increasingly studies have provided evidence of the importance of conformational changes to P450cam activity. Here we used infrared (IR) spectroscopy to investigate the molecular recognition of P450cam. Toward this goal, and to assess the influence of a hydrogen bond (H-bond) between active site residue Y96 and substrates, two variants in which Y96 is replaced by a cyanophenyl (Y96CNF) or phenyl (Y96F) group were characterized in complexes with the substrates camphor, isoborneol, and camphane. These combinations allow for a comparison of complexes in which the moieties on both the protein and substrate can serve as a H-bond donor, acceptor, or neither. The IR spectra of heme-bound CO and the site-specifically incorporated CN of Y96CNF were analyzed to characterize the number and nature of environments in each protein, both in the free and bound states. Although the IR spectra do not support the idea that protein-substrate H-bonding is central to P450cam recognition, the data altogether suggest that the differing conformational heterogeneity in the active site of the P450cam variants and changes in heterogeneity upon binding of different substrates likely contribute to their variable affinities via a conformational selection mechanism. This study further extends our understanding of the molecular recognition of archetypal P450cam and demonstrates the application of IR spectroscopy combined with selective protein modification to delineate protein-ligand interactions.

Direct local solvent probing by transient infrared spectroscopy reveals the mechanism of hydrogen-bond induced nonradiative deactivation

The hydrogen-bond induced quenching of an excited chromophore is visualised by probing O–H vibrations of the interacting solvent molecules.

The fluorescence quenching of organic dyes via H-bonding interactions is a well-known phenomenon. However, the mechanism of this Hydrogen-Bond Induced Nonradiative Deactivation (HBIND) is not understood. Insight into this process is obtained by probing in the infrared the O–H stretching vibration of the solvent after electronic excitation of a dye with H-bond accepting cyano groups. The fluorescence lifetime of this dye was previously found to decrease from 1.5 ns to 110 ps when going from an aprotic solvent to the strongly protic hexafluoroisopropanol (HFP). Prompt strengthening of the H-bond with the dye was identified by the presence of a broad positive O–H band of HFP, located at lower frequency than the O–H band of the pure solvent. Further strengthening occurs within a few picoseconds before the excited H-bonded complex decays to the ground state in 110 ps. The latter process is accompanied by the dissipation of energy from the dye to the solvent and the rise of a characteristic hot solvent band in the transient spectrum. Polarization-resolved measurements evidence a collinear alignment of the nitrile and hydroxyl groups in the H-bonded complex, which persists during the whole excited-state lifetime. Measurements in other fluorinated alcohols and in chloroform/HFP mixtures reveal that the HBIND efficiency depends not only on the strength of the H-bond interactions between the dye and the solvent but also on the ability of the solvent to form an extended H-bond network. The HBIND process can be viewed as an enhanced internal conversion of an excited complex consisting of the dye molecule connected to a large H-bond network.

Extended timescale 2D IR probes of proteins: p-cyanoselenophenylalanine

The importance of dynamics to the function of proteins is well appreciated, but the difficulty in their measurement impedes investigation into their precise role(s). 2D IR spectroscopy is a developing approach for the study of dynamics and has motivated efforts to develop spectrally resolved IR probe groups that enable its application for measuring the dynamics at specific sites in a protein. A challenge with this approach is that the timescales accessible are limited by the vibrational lifetimes of the probes. Toward development of better probes for 2D IR spectroscopy of protein dynamics, we report the characterization of p-cyano-seleno-phenylalanine (CNSePhe), a derivative of the well established IR probe p-cyano-phenylalanine (CNPhe), by FT IR, pump-probe, and 2D IR spectroscopy. The incorporation of the heavy Se atom decouples the CN vibration from the rest in the molecule. Although this leads to a reduction of the transition dipole strength, and thus a reduction in signal intensity, it also dramatically increases the vibrational lifetime, enabling collection of 2D IR spectra for analysis of molecular dynamics on much longer timescales. Interestingly, we also find that the lifetime for CNSePhe shows increased sensitivity to the presence of hydrogen bonding interactions with the CN, suggesting that the probe should be useful for interpretation of CN spectra and possibly for the study of solvation.

Detection of electron tunneling across plasmonic nanoparticle-film junctions using nitrile vibrations

The significant electric field enhancements that occur in plasmonic nanogap junctions are instrumental in boosting the performance of spectroscopy, optoelectronics and catalysis. Electron tunneling, associated with quantum effects in small junctions, is reported to limit the electric field enhancement. However, observing and quantitatively determining how tunneling alters the electric fields within small gaps is challenging due to the nanoscale dimensions and heterogeneity present experimentally. Here, we report the use of a nitrile probe placed in the nanoparticle-film gap junctions to demonstrate that the change in the nitrile stretching band associated with the vibrational Stark effect can be directly correlated with the local electric field environment modulated by gap size variations. The emergence of Stark shifts correlates with plasmon resonance shifts associated with electron tunneling across the gap junction. Time dependent changes in the nitrile band with extended illumination further support a build up of charge associated with optical rectification in the coupled plasmon system. Computational models agree with our experimental observations that the frequency shifts arise from a vibrational Stark effect. Large local electric fields associated with the smallest gap junctions give rise to significant Stark shifts. These results indicate that nitrile Stark probes can measure the local field strengths in plasmonic junctions and monitor the subtle changes in the local electric fields resulting from electron tunneling.