Carbon-deuterium bonds as probes of dihydrofolate reductase

Noncanonical Amino Acids in Biocatalysis

In recent years, powerful genetic code reprogramming methods have emerged that allow new functional components to be embedded into proteins as noncanonical amino acid (ncAA) side chains. In this review, we will illustrate how the availability of an expanded set of amino acid building blocks has opened a wealth of new opportunities in enzymology and biocatalysis research. Genetic code reprogramming has provided new insights into enzyme mechanisms by allowing introduction of new spectroscopic probes and the targeted replacement of individual atoms or functional groups. NcAAs have also been used to develop engineered biocatalysts with improved activity, selectivity, and stability, as well as enzymes with artificial regulatory elements that are responsive to external stimuli. Perhaps most ambitiously, the combination of genetic code reprogramming and laboratory evolution has given rise to new classes of enzymes that use ncAAs as key catalytic elements. With the framework for developing ncAA-containing biocatalysts now firmly established, we are optimistic that genetic code reprogramming will become a progressively more powerful tool in the armory of enzyme designers and engineers in the coming years.

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.

The structure of plastocyanin tunes the midpoint potential by restricting axial ligation of the reduced copper ion

Blue copper proteins are models for illustrating how proteins tune metal properties. Nevertheless, the mechanisms by which the protein controls the metal site remain to be fully elucidated. A hindrance is that the closed shell Cu(I) site is inaccessible to most spectroscopic analyses. Carbon deuterium (C-D) bonds used as vibrational probes afford nonperturbative, selective characterization of the key cysteine and methionine copper ligands in both redox states. The structural integrity of Nostoc plastocyanin was perturbed by disrupting potential hydrogen bonds between loops of the cupredoxin fold via mutagenesis (S9A, N33A, N34A), variably raising the midpoint potential. The C-D vibrations show little change to suggest substantial alteration to the Cu(II) coordination in the oxidized state or in the Cu(I) interaction with the cysteine ligand. They rather indicate, along with visible and NMR spectroscopy, that the methionine ligand distinctly interacts more strongly with the Cu(I) ion, in line with the increases in midpoint potential. Here we show that the protein structure determines the redox properties by restricting the interaction between the methionine ligand and Cu(I) in the reduced state.

Site-Selective Deuteration of Amino Acids through Dual-Protein Catalysis

Deuterated amino acids have been recognized for their utility in drug development, for facilitating nuclear magnetic resonance (NMR) analysis, and as probes for enzyme mechanism. Small molecule-based methods for the site-selective synthesis of deuterated amino acids typically involve de novo synthesis of the compound from deuterated precursors. In comparison, enzymatic methods for introducing deuterium offer improved efficiency, operating directly on free amino acids to achieve hydrogen-deuterium (H/D) exchange. However, site selectivity remains a significant challenge for enzyme-mediated deuteration, limiting access to desirable deuteration motifs. Here, we use enzyme-catalyzed deuteration, combined with steady-state kinetic analysis and ultraviolet (UV)-vis spectroscopy to probe the mechanism of a two-protein system responsible for the biosynthesis of l-allo-Ile. We show that an aminotransferase (DsaD) can pair with a small partner protein (DsaE) to catalyze Cα and Cβ H/D exchange of amino acids, while reactions without DsaE lead exclusively to Cα-deuteration. With conditions for improved catalysis, we evaluate the substrate scope for Cα/Cβ-deuteration and demonstrate the utility of this system for preparative-scale, selective labeling of amino acids.

Direct observation of peptide hydrogel self-assembly

The characterization of self-assembling molecules presents significant experimental challenges, especially when associated with phase separation or precipitation. Transparent window infrared (IR) spectroscopy leverages site-specific probes that absorb in the “transparent window” region of the biomolecular IR spectrum. Carbon–deuterium (C–D) bonds are especially compelling transparent window probes since they are non-perturbative, can be readily introduced site selectively into peptides and proteins, and their stretch frequencies are sensitive to changes in the local molecular environment. Importantly, IR spectroscopy can be applied to a wide range of molecular samples regardless of solubility or physical state, making it an ideal technique for addressing the solubility challenges presented by self-assembling molecules. Here, we present the first continuous observation of transparent window probes following stopped-flow initiation. To demonstrate utility in a self-assembling system, we selected the MAX1 peptide hydrogel, a biocompatible material that has significant promise for use in drug delivery and medical applications. C–D labeled valine was synthetically introduced into five distinct positions of the twenty-residue MAX1 β-hairpin peptide. Consistent with current structural models, steady-state IR absorption frequencies and linewidths of C–D bonds at all labeled positions indicate that these side chains occupy a hydrophobic region of the hydrogel and that the motion of side chains located in the middle of the hairpin is more restricted than those located on the hairpin ends. Following a rapid change in ionic strength to initiate self-assembly, the peptide absorption spectra were monitored as function of time, allowing determination of site-specific time constants. We find that within the experimental resolution, MAX1 self-assembly occurs as a cooperative process. These studies suggest that stopped-flow transparent window FTIR can be extended to other time-resolved applications, such as protein folding and enzyme kinetics.

To facilitate the characterization of phase-transitioning molecules, site-specific non-perturbative infrared probes are leveraged for continuous observation of the self-assembly of fibrils in a peptide hydrogel following stopped-flow initiation.

Altered coordination in a blue copper protein upon association with redox partner revealed by carbon-deuterium vibrational probes

Proteins tune the reactivity of metal sites; less understood is the impact of association with a redox partner. We demonstrate the utility of carbon-deuterium labels for selective analysis of delicate...

Infrared spectroscopy probes ion binding geometries

Infrared (IR) spectroscopy is a well-established technique for probing the structure, behavior, and surroundings of molecules in their native environments. Its characteristics-most specifically high structural sensitivity, ready applicability to aqueous samples, and broad availability-make it a valuable enzymological technique, particularly for the interrogation of ion binding sites. While IR spectroscopy of the "garden variety" (steady state at room temperature with wild-type proteins) is versatile and powerful in its own right, the combination of IR spectroscopy with specialized experimental schemes for leveraging ultrafast time resolution, protein labeling, and other enhancements further extends this utility. This book chapter provides the fundamental physical background and literature context essential for harnessing IR spectroscopy in the general context of enzymology with specific focus on interrogation of ion binding. Studies of lanthanide ions binding to calmodulin are highlighted as illustrative examples of this process. Appropriate sample preparation, data collection, and spectral interpretation are discussed from a detail-oriented and practical perspective with the goal of facilitating the reader's rapid progression from reading words in a book to collecting and analyzing their own data in the lab.

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.

Capturing Deuteration Effects in a Molecular Mechanics Force Field: Deuterated THF and the THF-Water Miscibility Gap

Deuteration is a common chemical modification used in conjunction with experiments such as neutron scattering, NMR, and Fourier-transform infrared for the study of molecular systems. Under the Born-Oppenheimer (BO) approximation, while the underlying potential energy surface remains unchanged by isotopic substitutions, isotopic substitution still alters intramolecular vibrations, which in turn may alter intermolecular interactions. Molecular mechanics (MM) force fields used in classical molecular dynamics simulations are assumed to represent local approximations of the BO potential energy surfaces, and hence, MD simulations using simple isotopic mass substitutions should capture BO-compatible isotope effects. However, standard MM force-field parameterizations do not directly fit to the local harmonic quantum mechanical (QM) Hessian that describes the BO surface, but rather to QM normal-modes and/or mass-dependent internal-coordinate derived distortion energies. Here, using tetrahydrofuran (THF)-water mixtures as our model system, we show that not only does a simple mass-substitution approach fail to capture an experimentally characterized deuteration effect (the loss of the closed-loop miscibility gap associated with the complete deuteration of THF) but also it is necessary to generate new MM force-field parameters that correctly describe isotopic dependent vibrations to capture the experimental deuteration effect. We show that the origin of this failure is a result of using mass-dependent features to fit the THF MM force field, which unintentionally biases the bonded terms of the force field to represent only the isotopologue used during the original force-field parameterization. In addition, we make use of our isotopologue-corrected force field for D₈THF to examine the molecular origins of the isotope-dependent loss of the THF-water miscibility gap.

Involvement of Local, Rapid Conformational Dynamics in Binding of Flexible Recognition Motifs

Flexible protein sequences populate ensembles of rapidly interconverting states differentiated by small-scale fluctuations; however, elucidating whether and how the ensembles determine function experimentally is challenged by the combined high spatial and temporal resolution needed to capture the states. We used carbon-deuterium (C-D) bond vibrations incorporated as infrared probes to characterize with residue-specific detail the heterogeneity of states adopted by proline-rich (PR) sequences and assess their involvement in recognition of Src homology 3 domains. The C-D absorption envelopes provided evidence for two or three sub-populations at all proline residues. The changes in the subpopulations induced by binding generally reflected recognition by conformational selection but depended on the residue and the state of the ligand to illuminate distinct mechanisms among the PR ligands. Notably, the spectral data indicate that greater adaptability among the states is associated with reduced recognition specificity and that perturbation to the ensemble populations contributes to differences in binding entropy. Broadly, the study quantifies rapidly interconverting ensembles with residue-specific detail and implicates them in function.

Site-Specific 1D and 2D IR Spectroscopy to Characterize the Conformations and Dynamics of Protein Molecular Recognition

Proteins exist as ensembles of interconverting states on a complex energy landscape. A complete, molecular-level understanding of their function requires knowledge of the populated states and thus the experimental tools to characterize them. Infrared (IR) spectroscopy has an inherently fast time scale that can capture all states and their dynamics with, in principle, bond-specific spatial resolution, and 2D IR methods that provide richer information are becoming more routine. Although application of IR spectroscopy for investigation of proteins is challenged by spectral congestion, the issue can be overcome by site-specific introduction of amino acid side chains that have IR probe groups with frequency-resolved absorptions, which furthermore enables selective characterization of different locations in proteins. Here, we briefly introduce the biophysical methods and summarize the current progress toward the study of proteins. We then describe our efforts to apply site-specific 1D and 2D IR spectroscopy toward elucidation of protein conformations and dynamics to investigate their involvement in protein molecular recognition, in particular mediated by dynamic complexes: plastocyanin and its binding partner cytochrome f, cytochrome P450s and substrates or redox partners, and Src homology 3 domains and proline-rich peptide motifs. We highlight the advantages of frequency-resolved probes to characterize specific, local sites in proteins and uncover variation among different locations, as well as the advantage of the fast time scale of IR spectroscopy to detect rapidly interconverting states. In addition, we illustrate the greater insight provided by 2D methods and discuss potential routes for further advancement of the field of biomolecular 2D IR spectroscopy.

pKa Determination of a Histidine Residue in a Short Peptide Using Raman Spectroscopy

Determining the pKa of key functional groups is critical to understanding the pH-dependent behavior of biological proteins and peptide-based biomaterials. Traditionally, ¹H NMR spectroscopy has been used to determine the pKa of amino acids; however, for larger molecules and aggregating systems, this method can be practically impossible. Previous studies concluded that the C-D stretches in Raman are a useful alternative for determining the pKa of histidine residues. In this study, we report on the Raman application of the C2-D probe on histidine’s imidazole side chain to determining the pKa of histidine in a short peptide sequence. The pKa of the tripeptide was found via difference Raman spectroscopy to be 6.82, and this value was independently confirmed via ¹H NMR spectroscopy on the same peptide. The C2-D probe was also compared to other Raman reporters of the protonation state of histidine and was determined to be more sensitive and reliable than other protonation-dependent signals. The C2-D Raman probe expands the tool box available to chemists interested in directly interrogating the pKa’s of histidine-containing peptide and protein systems.

Reduction of Folate by Dihydrofolate Reductase from Thermotoga maritima

Mammalian dihydrofolate reductases (DHFRs) catalyze the reduction of folate more efficiently than the equivalent bacterial enzymes do, despite typically having similar efficiencies for the reduction of their natural substrate, dihydrofolate. In contrast, we show here that DHFR from the hyperthermophilic bacterium Thermotoga maritima can catalyze reduction of folate to tetrahydrofolate with an efficiency similar to that of reduction of dihydrofolate under saturating conditions. Nuclear magnetic resonance and mass spectrometry experiments showed no evidence of the production of free dihydrofolate during either the EcDHFR- or TmDHFR-catalyzed reductions of folate, suggesting that both enzymes perform the two reduction steps without release of the partially reduced substrate. Our results imply that the reaction proceeds more efficiently in TmDHFR than in EcDHFR because the more open active site of TmDHFR facilitates protonation of folate. Because T. maritima lives under extreme conditions where tetrahydrofolate is particularly prone to oxidation, this ability to salvage folate may impart an advantage to the bacterium by minimizing the squandering of a valuable cofactor.

An expanded genetic code for probing the role of electrostatics in enzyme catalysis by vibrational Stark spectroscopy

BACKGROUND

To find experimental validation for electrostatic interactions essential for catalytic reactions represents a challenge due to practical limitations in assessing electric fields within protein structures.

SCOPE OF REVIEW

This review examines the applications of non-canonical amino acids (ncAAs) as genetically encoded probes for studying the role of electrostatic interactions in enzyme catalysis.

MAJOR CONCLUSIONS

ncAAs constitute sensitive spectroscopic probes to detect local electric fields by exploiting the vibrational Stark effect (VSE) and thus have the potential to map the protein electrostatics.

GENERAL SIGNIFICANCE

Mapping the electrostatics in proteins will improve our understanding of natural catalytic processes and, in beyond, will be helpful for biocatalyst engineering. This article is part of a Special Issue entitled "Biochemistry of Synthetic Biology - Recent Developments" Guest Editor: Dr. Ilka Heinemann and Dr. Patrick O'Donoghue.

Protein motions and dynamic effects in enzyme catalysis

The role of protein motions in promoting the chemical step of enzyme catalysed reactions remains a subject of considerable debate. Here, a unified view of the role of protein dynamics in dihydrofolate reductase catalysis is described. Recently the role of such motions has been investigated by characterising the biophysical properties of isotopically substituted enzymes through a combination of experimental and computational analyses. Together with previous work, these results suggest that dynamic coupling to the chemical coordinate is detrimental to catalysis and may have been selected against during DHFR evolution. The full catalytic power of Nature's catalysts appears to depend on finely tuning protein motions in each step of the catalytic cycle.

Methionine Ligand Interaction in a Blue Copper Protein Characterized by Site-Selective Infrared Spectroscopy

The reactivity of metal sites in proteins is tuned by protein-based ligands. For example, in blue copper proteins such as plastocyanin (Pc), the structure imparts a highly elongated bond between the Cu and a methionine (Met) axial ligand to modulate its redox properties. Despite extensive study, a complete understanding of the contribution of the protein to redox activity is challenged by experimentally accessing both redox states of metalloproteins. Using infrared (IR) spectroscopy in combination with site-selective labeling with carbon-deuterium (C-D) vibrational probes, we characterized the localized changes at the Cu ligand Met97 in the oxidized and reduced states, as well as the Zn(II) or Co(II)-substituted, the pH-induced low-coordinate, the apoprotein, and the unfolded states. The IR absorptions of (d3-methyl)Met97 are highly sensitive to interaction of the sulfur-based orbitals with the metal center and are demonstrated to be useful reporters of its modulation in the different states. Unrestricted Kohn-Sham density functional theory calculations performed on a model of the Cu site of Pc confirm the observed dependence. IR spectroscopy was then applied to characterize the impact of binding to the physiological redox partner cytochrome (cyt) f. The spectral changes suggest a slightly stronger Cu-S(Met97) interaction in the complex with cyt f that has potential to modulate the electron transfer properties. Besides providing direct, molecular-level comparison of the oxidized and reduced states of Pc from the perspective of the axial Met ligand and evidence for perturbation of the Cu site properties by redox partner binding, this study demonstrates the localized spatial information afforded by IR spectroscopy of selectively incorporated C-D probes.

Influence of Equatorial CH⋅⋅⋅O Interactions on Secondary Kinetic Isotope Effects for Methyl Transfer

DFT calculations for methyl cation complexed within a constrained cage of water molecules permit the controlled manipulation of the "axial" donor/acceptor distance and the "equatorial" distance to hydrogen-bond acceptors. The kinetic isotope effect k(CH3)/k(CT3) for methyl transfer within a cage with a short axial distance becomes less inverse for shorter equatorial C⋅⋅⋅O distances: a decrease of 0.5 Å results in a 3 % increase at 298 K. Kinetic isotope effects in AdoMet-dependent methyltransferases may be m∧odulated by CH⋅⋅⋅O hydrogen bonding, and factors other than axial compression may contribute, at least partially, to recently reported isotope-effect variations for catechol-O-methyltransferase and its mutant structures.

Carbon-deuterium bonds as non-perturbative infrared probes of protein dynamics, electrostatics, heterogeneity, and folding

Vibrational spectroscopy is uniquely able to characterize protein dynamics and microenvironmental heterogeneity because it possesses an inherently high temporal resolution and employs probes of ultimately high structural resolution-the bonds themselves. The use of carbon-deuterium (C-D) bonds as vibrational labels circumvents the spectral congestion that otherwise precludes the use of vibrational spectroscopy to proteins and makes the observation of single vibrations within a protein possible while being wholly non-perturbative. Thus, C-D probes can be used to site-specifically characterize conformational heterogeneity and thermodynamic stability. C-D probes are also uniquely useful in characterizing the electrostatic microenvironment experienced by a specific residue side chain or backbone due to its effect on the C-D absorption frequency. In this chapter we describe the experimental procedures required to use C-D bonds and FT IR spectroscopy to characterize protein dynamics, structural and electrostatic heterogeneity, ligand binding, and folding.

Experimental characterization of electrostatic and conformational heterogeneity in an SH3 domain

Electrostatic and conformational heterogeneity make central contributions to protein function, but their experimental characterization requires a combination of spatial and temporal resolution that is challenging to achieve. Src homology 3 (SH3) domains mediate protein-protein interactions, and NMR studies have demonstrated that most possess conformational heterogeneity, which could be critical for their function. Here, we use the IR absorptions of carbon-deuterium (C-D) bonds site-selectively incorporated throughout the N-terminal SH3 domain from the murine adapter protein Crk-II to characterize its different microenvironments with high spatial and temporal resolution. The C-D absorptions are only differentiated in the folded state of the protein where they show evidence of significant environmental heterogeneity. However, the spectra of the folded state are independent of temperature, and upon thermal denaturation the protein undergoes a single, global unfolding transition. While some evidence of conformational heterogeneity is found within the peptide backbone, the majority of the environmental heterogeneity appears to result from electrostatics.

Vibrational solvatochromism and electrochromism. II. Multipole analysis

Small infrared probe molecules have been widely used to study local electrostatic environment in solutions and proteins. Using a variety of time- and frequency-resolved vibrational spectroscopic methods, one can accurately measure the solvation-induced vibrational frequency shifts and the timescales and amplitudes of frequency fluctuations of such IR probes. Since the corresponding frequency shifts are directly related to the local electric field and its spatial derivatives of the surrounding solvent molecules or amino acids in proteins, one can extract information on local electric field around an IR probe directly from the vibrational spectroscopic results. Here, we show that, carrying out a multipole analysis of the solvatochromic frequency shift, the solvatochromic dipole contribution to the frequency shift is not always the dominant factor. In the cases of the nitrile-, thiocyanato-, and azido-derivatized molecules, the solvatochromic quadrupole contributions to the corresponding stretch mode frequency shifts are particularly large and often comparable to the solvatochromic dipole contributions. Noting that the higher multipole moment-solvent electric field interactions are short range effects in comparison to the dipole interaction, the H-bonding interaction-induced vibrational frequency shift can be caused by such short-range multipole-field interaction effects. We anticipate that the present multipole analysis method specifically developed to describe the solvatochromic vibrational frequency shifts will be useful to understand the intermolecular interaction-induced vibrational property changes and to find out a relationship between vibrational solvatochromism and electrochromism of IR probes in condensed phases.

Monitoring Intramolecular Proton Transfer with Two-Dimensional Infrared Spectroscopy: A Computational Prediction

Proton transfer processes are ubiquitous and play a vital role in a broad range of chemical and biochemical phenomena. The ability of two-dimensional infrared (2D IR) spectroscopy with a carbon-deuterium (C-D) reporter to monitor the kinetics of proton transfer in the model compound malonaldehyde was demonstrated computationally. One of the two carbonyl/enol carbon atoms in malonaldehyde was labeled with a C-D bond. The C-D stretch vibrational frequency provides ∼150 cm(-1) of sensitivity to the two tautomers of malonaldehyde. Mixed quantum mechanics/molecular mechanics simulations employing the self-consistent-charge density functional tight binding (SCC-DFTB) method were used to compute 2D IR line shapes for the C-D stretch of labeled malonaldehyde in aqueous solution. The 2D IR spectra reveal cross peaks from the chemical exchange of the proton. The kinetics for the growth of the cross-peaks (and the decay of the diagonal peaks) precisely match the proton transfer rate observed in the SCC-DFTB simulations.

Carbon-deuterium bonds as probes of protein thermal unfolding

We report a residue-specific characterization of the thermal unfolding mechanism of ferric horse heart cytochrome c using C-D bonds site-specifically incorporated at residues dispersed throughout three different structural elements within the protein. As the temperature increases, Met80 first dissociates from the heme center, and the protein populates a folding intermediate before transitioning to a solvent exposed state. With further increases in temperature, the C-terminal helix frays and then loses structure along with the core of the protein. Interestingly, the data also reveal that the state populated at high temperature retains some structure and possibly represents a molten globule. Elucidation of the detailed unfolding mechanism and the structure of the associated molten globule, both of which represent challenges to conventional techniques, highlights the utility of the C-D technique.

Electrostatic effects of mutations of Ras glutamine 61 measured using vibrational spectroscopy of a thiocyanate probe

Mutations of human oncoprotein p21(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 electrostatically stabilizes the formation of the charged transition state at the γ-phosphate during hydrolysis. We have tested the interactions between amino acids at this position and water by measuring changes in the electrostatic field experienced by a nitrile probe positioned near Ras Q61 using vibrational Stark effect (VSE) spectroscopy. We mutated this glutamine to every amino acid except cysteine and proline and then incubated these mutants with a Ral guanine nucleotide dissociation stimulator (Ral) containing the I18C mutation that was chemically labeled with a thiocyanate vibrational spectroscopic probe. The formation of the docked Ras Q61X-labeled Ral complex was confirmed by measurement of the dissociation constant of the interaction. We measured the absorption energy of this nitrile to determine any differences in electrostatic environment in the immediate vicinity of the thiocyanate probe between wild type and mutants of Ras. For each Ras Q61X mutant, we correlate the change in electrostatic field at position 61 with the solvent accessible surface area of polar components of the mutant side chain determined from a Boltzmann-weighted ensemble of structures, as well as the residue's hydration potential. These results support the hypothesis that the role of Ras Q61 is to stabilize water in or near the active site during GTP hydrolysis. The substantial effect that nonpolar side chains of Ras Q61X have on the absorption energy of the thiocyanate must be investigated with further experiments.

Site-Specific Spectroscopic Reporters of the Local Electric Field, Hydration, Structure, and Dynamics of Biomolecules

Elucidating the underlying molecular mechanisms of protein folding and function is a very exciting and active research area, but poses significant challenges. This is due in part to the fact that existing experimental techniques are incapable of capturing snapshots along the 'reaction coordinate' in question with both sufficient spatial and temporal resolutions. In this regard, recent years have seen increased interests and efforts in development and employment of site-specific probes to enhance the structural sensitivity of spectroscopic techniques in conformational and dynamical studies of biological molecules. In particular, the spectroscopic and chemical properties of nitriles, thiocyanates, and azides render these groups attractive for the interrogation of complex biochemical constructs and processes. Here, we review their signatures in vibrational, fluorescence and NMR spectra and their utility in the context of elucidating chemical structure and dynamics of protein and DNA molecules.

Two-dimensional infrared spectroscopy of azido-nicotinamide adenine dinucleotide in water

Mid-IR active analogs of enzyme cofactors have the potential to be important spectroscopic reporters of enzyme active site dynamics. Azido-nicotinamide adenine dinucleotide (NAD(+)), which has been recently synthesized in our laboratory, is a mid-IR active analog of NAD(+), a ubiquitous redox cofactor in biology. In this study, we measure the frequency-frequency time correlation function for the antisymmetric stretching vibration of the azido group of azido-NAD(+) in water. Our results are consistent with previous studies of pseudohalides in water. We conclude that azido-NAD(+) is sensitive to local environmental fluctuations, which, in water, are dominated by hydrogen-bond dynamics of the water molecules around the probe. Our results demonstrate the potential of azido-NAD(+) as a vibrational probe and illustrate the potential of substituted NAD(+)-analogs as reporters of local structural dynamics that could be used for studies of protein dynamics in NAD-dependent enzymes.

Effect of pH on hydride transfer by Escherichia coli dihydrofolate reductase

The kinetic isotope effect (KIE) on hydride transfer in the reaction catalysed by dihydrofolate reductase from Escherichia coli (EcDHFR) is known to be temperature dependent at pH 7, but essentially independent of temperature at elevated pH. Here, we show that the transition from the temperature-dependent regime to the temperature-independent regime occurs sharply between pH 7.5 and 8. The activation energy for hydride transfer is independent of pH. The mechanism leading to the change in behaviour of the KIEs is not clear, but probably involves a conformational change in the enzyme brought about by deprotonation of a key residue (or residues) at high pH. The KIE on hydride transfer at low pH suggests that the rate constant for the reaction is not limited by a conformational change to the enzyme under these conditions. The effect of pH on the temperature dependence of the rate constants and KIEs for hydride transfer catalysed by EcDHFR suggests that enzyme motions and conformational changes do not directly influence the chemistry, but that the reaction conditions affect the conformational ensemble of the enzyme prior to reaction and control the reaction though this route.

Nitrile bonds as infrared probes of electrostatics in ribonuclease S

Three different nitrile-containing amino acids, p-cyanophenylalanine, m-cyanophenylalanine, and S-cyanohomocysteine, have been introduced near the active site of the semisynthetic enzyme ribonuclease S (RNase S) to serve as probes of electrostatic fields. Vibrational Stark spectra, measured directly on the probe-modified proteins, confirm the predominance of the linear Stark tuning rate in describing the sensitivity of the nitrile stretch to external electric fields, a necessary property for interpreting observed frequency shifts as a quantitative measure of local electric fields that can be compared with simulations. The X-ray structures of these nitrile-modified RNase variants and enzymatic assays demonstrate minimal perturbation to the structure and function, respectively, by the probes and provide a context for understanding the influence of the environment on the nitrile stretching frequency. We examine the ability of simulation techniques to recapitulate the spectroscopic properties of these nitriles as a means to directly test a computational electrostatic model for proteins, specifically that in the ubiquitous Amber-99 force field. Although qualitative agreement between theory and experiment is observed for the largest shifts, substantial discrepancies are observed in some cases, highlighting the ongoing need for experimental metrics to inform the development of theoretical models of electrostatic fields in proteins.

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

Infrared (IR) band shifts of isolated vibrational transitions can serve as quantitative and directional probes of local electrostatic fields, due to the vibrational Stark effect. However, departures from the Stark model can arise when the probe participates in specific, chemical interactions, such as direct hydrogen bonding. We present a method to identify and correct for these departures based on comparison of (13)C NMR chemical shifts and IR frequencies each calibrated in turn by a solvatochromic model. We demonstrate how the tandem use of these experimental observables can be applied to a thiocyanate-modified protein, ketosteroid isomerase, and show, by comparison to structural models, that changes in electrostatic field can be measured within the complex protein environment even in the background of direct hydrogen bonding to the probe.

The effects of alpha-helical structure and cyanylated cysteine on each other

β-Thiocyanatoalanine, or cyanylated cysteine, is an artificial amino acid that can be introduced at solvent-exposed cysteine residues in proteins via chemical modification. Its facile post-translational synthesis means that it may find broad use in large protein systems as a probe of site-specific structure and dynamics. The C≡N stretching vibration of this artificial side chain provides an isolated infrared chromophore. To test both the perturbative effect of this side chain on local secondary structure and its sensitivity to structural changes, three variants of a model water-soluble alanine-repeat helix were synthesized containing cyanylated cysteine at different sites. The cyanylated cysteine side chain is shown to destabilize, but not completely disrupt, the helical structure of the folded peptide when substituted for alanine. In addition, the C≡N stretching bandwidth of the artificial side chain is sensitive to the helix−coil structural transition. These model system results indicate that cyanylated cysteine can be placed into protein sequences with a native helical propensity without destroying the helix, and further that the CN probe may be able to report local helix formation events even when it is water-exposed in both the ordered and disordered conformational states. These results indicate that cyanylated cysteine could be a widely useful probe of structure-forming events in proteins with large in vitro structural distributions.

Efforts toward developing probes of protein dynamics: vibrational dephasing and relaxation of carbon-deuterium stretching modes in deuterated leucine

The spectral position of C-D stretching absorptions in the so-called "transparent window" of protein absorption (1800-2300 cm(-1)) makes them well suited as probes of protein dynamics with high temporal and structural resolution. We have previously incorporated single deuterated amino acids into proteins to site-selectively follow protein folding and ligand binding by steady-state FT IR spectroscopy. Ultimately, our goal is to use C-D bonds as probes in time-resolved IR spectroscopy to study dynamics and intramolecular vibrational energy redistribution (IVR) in proteins. As a step toward this goal, we now present the first time-resolved experiments characterizing the population and dephasing dynamics of selectively excited C-D bonds in a deuterated amino acid. Three differently deuterated, Boc-protected leucines were selected to systematically alter the number of additional C-D bonds that may mediate IVR out of the initially populated bright C-D stretching mode. Three-pulse photon echo experiments show that the steady-state C-D absorption linewidths are broadened by both homogeneous and inhomogeneous effects, and transient grating experiments reveal that IVR occurs on a subpicosecond time scale and is nonstatistical. The results have important implications for the interpretation of steady-state C-D spectra and demonstrate the potential utility of C-D bonds as probes of dynamics and IVR within a protein.

Direct observation of structural heterogeneity in a beta-sheet

Structural heterogeneity is thought to be inherent in many proteins and may be important for their folding and/or function. However, it is difficult to detect by conventional methods. Carbon-deuterium bonds are environmentally sensitive, nonperturbative probes of protein environments whose observation and characterization are facilitated by their unique stretching absorption frequency in an otherwise unobscured region of the IR spectrum. We demonstrate that deuterium atoms incorporated at C(alpha) backbone positions (C(alpha)-D bonds) are sensitive to the local backbone structure and thus may be used not only to detect structural heterogeneity but also to help characterize it structurally. Density functional theory calculations are used to predict that C(alpha)-D bonds of glycine are sensitive to their local structure, with the absorptions red-shifted for an extended beta-sheet relative to gamma- and alpha-helix-like turns. These predictions are confirmed using the N-terminal Src homology 3 (nSH3) domain from the human CrkII adaptor protein, whose function as a signaling domain may require structural heterogeneity. Four nSH3 variants were synthesized in which individual glycine residues were site-specifically modified with C(alpha)D(2) glycine residues. Not only were the C(alpha)-D bonds incorporated within the beta-sheet of nSH3 more red-shifted than those incorporated within loops, but the data also reveal that nSH3 populates at least two discrete beta-sheet core structures. Moreover, the data suggest that the folded core of nSH3 may be less ordered than previously believed and also that the unfolded state may be more ordered than previously thought, and both of these factors may influence the folding and function of these important signaling domains. The C-D-based IR technique should be generally useful in the characterization of structure and heterogeneity of both folded and unfolded states.

Efforts toward the direct experimental characterization of enzyme microenvironments: tyrosine100 in dihydrofolate reductase

State secrets: Site-specific deuteration and FTIR studies reveal that Tyr100 in dihydrofolate reductase plays an important role in catalysis, with a strong electrostatic coupling occurring between Tyr100 and the charge that develops in the hydride-transfer transition state (see picture, NADP(+) purple, Tyr100 green). However, relaying correlated motions that facilitate catalysis from distal sites of the protein to the hydride donor may also be involved.

The determinants of stability and folding in evolutionarily diverged cytochromes c

Cytochrome c has served as a paradigm for the study of protein stability, folding, and molecular evolution, but it remains unclear how these aspects of the protein are related. For example, while the bovine and equine cytochromes c are known to have different stabilities, and possibly different folding mechanisms, it is not known how these differences arise from just three amino acid substitutions introduced during divergence. Using site-selectively incorporated carbon-deuterium bonds, we show that like the equine protein, bovine cytochrome c is induced to unfold by guanidine hydrochloride via a stepwise mechanism, but it does not populate an intermediate as is observed with the equine protein. The increased stability also results in more similar free energies of unfolding observed at different sites within the protein, giving the appearance of a more concerted mechanism. Furthermore, we show that the differences in stability and folding appear to result from a single amino acid substitution that stabilizes a helix by allowing for increased solvation of its N-terminus.