Direct and High Resolution Characterization of Cytochrome c Equilibrium Folding

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.

Genetically Encoded Fluorescent Probe for Detection of Heme-Induced Conformational Changes in Cytochrome c

Cytochrome c (Cytc) is a key redox protein for energy metabolism and apoptosis in cells. The activation of Cytc is composed of several steps, including its transfer to the mitochondrial membrane, binding to cytochrome c heme lyase (CCHL) and covalent attachment to heme. The spectroscopic methods are often applied to study the structural changes of Cytc. However, they require the isolation of Cytc from cells and have limited availability under physiological conditions. Despite recent studies to elucidate the tightly regulated folding mechanism of Cytc, the role of these events and their association with different conformational states remain elusive. Here, we provide a genetically encoded fluorescence method that allows monitoring of the conformational changes of Cytc upon binding to heme and CCHL. Cerulean and Venus fluorescent proteins attached at the N and C terminals of Cytc can be used to determine its unfolded, intermediate, and native states by measuring FRET amplitude. We found that the noncovalent interaction of heme in the absence of CCHL induced a shift in the FRET signal, indicating the formation of a partially folded state. The higher concentration of heme and coexpression of CCHL gave rise to the recovery of Cytc native structure. We also found that Cytc was weakly associated with CCHL in the absence of heme. As a result, a FRET-based fluorescence approach was demonstrated to elucidate the mechanism of heme-induced Cytc conformational changes with spatiotemporal resolution and can be applied to study its interaction with small molecules and other protein partners in living cells.

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.

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...

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.

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.

Thermally induced protein unfolding probed by isotope-edited IR spectroscopy

Infrared (IR) spectroscopy has been widely utilized for the study of protein folding, unfolding, and misfolding processes. We have previously developed a theoretical method for calculating IR spectra of proteins in the amide I region. In this work, we apply this method, in combination with replica-exchange molecular dynamics simulations, to study the equilibrium thermal unfolding transition of the villin headpiece subdomain (HP36). Temperature-dependent IR spectra and spectral densities are calculated. The spectral densities correctly reflect the unfolding conformational changes in the simulation. With the help of isotope labeling, we are able to capture the feature that helix 2 of HP36 loses its secondary structure before global unfolding occurs, in agreement with experiment.

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.

Native state dynamics drive the unfolding of the SH3 domain of PI3 kinase at high denaturant concentration

Little is known about the role of protein dynamics in directing protein unfolding along a specific pathway and about the role played by chemical denaturants in modulating the dynamics and the initiation of unfolding. In this study, deuterium-hydrogen exchange (HX) detected by electrospray ionization mass spectrometry (ESI-MS) was used to study the unfolding of the SH3 domain of the PI3 kinase. Unfolding on the principal unfolding pathway occurs in 2 steps, both in the absence and in the presence of 1.8 M guanidine hydrochloride (GdnHCl). In both cases, the first step leads to the formation of an intermediate, I(N), with 5 fewer protected amide hydrogen sites than in N. In the second step, I(N) loses the structure protecting the remaining 14 amide hydrogen sites from HX as it unfolds completely. ESI-MS analysis of fragments of the protein created by proteolytic digestion, after completion of the HX reaction, shows that I(N) has lost protection against HX in the same segments of native structure during unfolding in the absence and presence of 1.8 M GdnHCl. Hence, GdnHCl does not appear to play a direct active role in the initiation of unfolding. However, at higher GdnHCl concentrations, a second unfolding pathway is shown to compete effectively with the N <--> I(N) <--> U pathway. In this way, the denaturant modulates the energy landscape of unfolding.

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.

C-D modes of deuterated side chain of leucine as structural reporters via dual-frequency two-dimensional infrared spectroscopy

Perdeuteration of the side chains of amino acids such as leucine results in appearance of reasonably strong absorption peaks around 2050-2220 cm(-1) that belong to the CD stretching modes and exhibit extinction coefficients of up to 120 M(-1) cm(-1). The properties of the CD stretching transitions in leucine-d(10) as structural labels are studied via the methods of two-dimensional infrared (2DIR) spectroscopy. The cross peaks caused by interactions of the CD stretching modes with amide I (Am-I), CO, and amide II (Am-II) modes are obtained by the dual-frequency 2DIR method. The CD stretching peaks in leucine-d(10) are characterized using DFT computational modeling as well as relaxation-assisted 2DIR (RA 2DIR) measurements. The RA 2DIR measurements showed different enhancements and different energy transport times (arrival times) for different CD/Am-II and CD/CO cross peaks; a correlation between the intermode distance, the arrival time, and the amplification factor is reported. We demonstrated that the CD transitions of leucine-d(10) amino acid can serve as convenient structural reporters via the dual-frequency 2DIR methods and discussed the potential of leucine-d(10) and other amino acids with deuterium-labeled side chains as probes of protein structure and dynamics.

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.

Characterization of alkaline transitions in ferricytochrome c using carbon-deuterium infrared probes

The alkaline-induced structural transitions of ferricytochrome c have been studied intensively as a model for how changes in metal ligation contribute to protein function and folding. Previous studies have demonstrated that multiple non-native species accumulate with increasing pH. Here, we used a combination of experiments and simulations to provide a high-resolution view of the changes associated with increasing alkaline conditions. Alkaline-induced transitions were characterized under equilibrium conditions by following changes in the IR absorptions of carbon-deuterium chromophores incorporated at Leu68, Lys72, Lys73, Lys79, and Met80. The data suggest that at least four intermediates are formed as the pH is increased prior to complete unfolding of the protein. The first alkaline transition observed appears to be driven by a single deprotonation and occurs with a midpoint of pH 8.8, but surprisingly, the intermediate formed does not appear to be one of the well-characterized lysine misligates. At higher pH, second and third deprotonations, with a combined apparent midpoint pH of 10.2, induce transitions to Lys73- or Lys79-misligated species. Interestingly, the lysine misligates appear to undergo iron reduction by the coordinated amine. A transition from the lysine misligates to another intermediate, likely a hydroxide-misligated species, is associated with a fourth deprotonation and a midpoint of pH 10.7. Finally, the protein loses tertiary structure with a fifth deprotonation that occurs with a midpoint of pH 12.7. Native topology-based models with enforced misligation are employed to help understand the structures of the observed intermediates.

The conformational manifold of ferricytochrome c explored by visible and far-UV electronic circular dichroism spectroscopy

The oxidized state of cytochrome c is a subject of continuous interest, owing to the multitude of conformations which the protein can adopt in solution and on surfaces of artificial and cell membranes. The structural diversity corresponds to a variety of functions in electron transfer, peroxidase and apoptosis processes. In spite of numerous studies, a comprehensive analysis and comparison of native and non-native states of ferricytochrome c has thus far not been achieved. This results in part from the fact that the influence of solvent conditions (i.e., ionic strength, anion concentration, temperature dependence of pH values) on structure, function and equilibrium thermodynamics has not yet been thoroughly assessed. The current study is a first step in this direction, in that it provides the necessary experimental data to compare different non-native states adopted at high temperature and alkaline pH. To this end, we employed visible electronic circular dichroism (ECD) and absorption spectroscopy to probe structural changes of the heme environment in bovine and horse heart ferricytochrome c as a function of temperature between 278 and 363 K at different neutral and alkaline pH values. A careful selection of buffers enabled us to monitor the partial unfolding of the native state at room temperature while avoiding a change to an alkaline state at high temperatures. We found compelling evidence for the existence of a thermodynamic intermediate of the thermal unfolding/folding process, termed III h, which is structurally different from the alkaline states, IV 1 and IV 2, contrary to current belief. At neutral or slightly acidic pH, III h is populated in a temperature region between 320 and 345 K. The unfolded state of the protein becomes populated at higher temperatures. The ECD spectra of the B-bands of bovine and horse heart cytochrome c (pH 7.0) exhibit a pronounced couplet that is maintained below 343 K, before protein unfolding replaces it by a rather strong positive Cotton band. A preliminary vibronic analysis of the B-band profile reveals that the couplet reflects a B-band splitting of 350 cm (-1), which is mostly of electronic origin, due to the internal electric field in the heme cavity. Our results suggest that the conformational transition from the native state, III, into a thermally activated intermediate state, III h, does not substantially affect the internal electric field and causes only moderate rearrangements of the heme pocket, which involves changes, rather than a rupture, of the Fe (3+)-M80 linkage. In the unfolded state, as well as in the alkaline states IV and V, the band splitting is practically eliminated, but the positive Cotton effect observed for the B-band suggests that the proximal environment, encompassing H18 and the two cysteine residues 14 and 17, is most likely still intact and covalently bound to the heme chromophore. Both alkaline states IV and V were found to melt via intermediate states. Unfolded states probed at neutral and alkaline pH can be discriminated, owing to the different intensities of the Cotton bands of the respective B-band transitions. Differences between the ECD intensities of the B-bands of the different unfolded states and alkaline states most likely reflect different degrees of openness of the corresponding heme crevice.

Conformational substates of horse heart cytochrome c exhibit different thermal unfolding of the heme cavity

The charge transfer (CT) band at 695 nm in the spectrum of ferri-cytochrome c is highly asymmetric, indicating conformational heterogeneity due to the coexistence of different conformational substates. We have measured the respective band profile of horse heart ferri-cytochrome c as a function of temperature between 283 K (10 degrees C) and 333 K (60 degrees C) and found that the well-known decrease of the absorptivity is wavenumber-dependent and exhibits a biphasic behavior. This indicates that the underlying conformational substates differ in their thermodynamic stability with respect to the structural changes associated with the disappearance of the 695 nm band, which eventually (at high temperatures) involves the replacement of M80 by a nearby lysine residue. Our data further indicates that the thermal unfolding process involves two structurally different intermediate states.