On the nature of conformational openings: native and unfolded-state hydrogen and thiol-disulfide exchange studies of ferric aquomyoglobin

Mass Spectrometry-Based Protein Footprinting for Higher-Order Structure Analysis: Fundamentals and Applications

Proteins adopt different higher-order structures (HOS) to enable their unique biological functions. Understanding the complexities of protein higher-order structures and dynamics requires integrated approaches, where mass spectrometry (MS) is now positioned to play a key role. One of those approaches is protein footprinting. Although the initial demonstration of footprinting was for the HOS determination of protein/nucleic acid binding, the concept was later adapted to MS-based protein HOS analysis, through which different covalent labeling approaches "mark" the solvent accessible surface area (SASA) of proteins to reflect protein HOS. Hydrogen-deuterium exchange (HDX), where deuterium in D2O replaces hydrogen of the backbone amides, is the most common example of footprinting. Its advantage is that the footprint reflects SASA and hydrogen bonding, whereas one drawback is the labeling is reversible. Another example of footprinting is slow irreversible labeling of functional groups on amino acid side chains by targeted reagents with high specificity, probing structural changes at selected sites. A third footprinting approach is by reactions with fast, irreversible labeling species that are highly reactive and footprint broadly several amino acid residue side chains on the time scale of submilliseconds. All of these covalent labeling approaches combine to constitute a problem-solving toolbox that enables mass spectrometry as a valuable tool for HOS elucidation. As there has been a growing need for MS-based protein footprinting in both academia and industry owing to its high throughput capability, prompt availability, and high spatial resolution, we present a summary of the history, descriptions, principles, mechanisms, and applications of these covalent labeling approaches. Moreover, their applications are highlighted according to the biological questions they can answer. This review is intended as a tutorial for MS-based protein HOS elucidation and as a reference for investigators seeking a MS-based tool to address structural questions in protein science.

Global Implications of Local Unfolding Phenomena, Probed by Cysteine Reactivity in Human Frataxin

Local events that affect specific regions of proteins are of utmost relevance for stability and function. The aim of this study is to quantitatively assess the importance of locally-focused dynamics by means of a simple chemical modification procedure. Taking human Frataxin as a working model, we investigated local fluctuations of the C-terminal region (the last 16 residues of the protein) by means of three L → C replacement mutants: L98C, L200C and L203C. The conformation and thermodynamic stability of each variant was assessed. All the variants exhibited native features and high stabilities: 9.1 (wild type), 8.1 (L198C), 7.0 (L200C) and 10.0 kcal mol^-1 (L203C). In addition, kinetic rates of Cys chemical modification by DTNB and DTDPy were measured, conformational dynamics data were extracted and free energy for the local unfolding of the C-terminal region was estimated. The analysis of these results indicates that the conformation of the C-terminal region fluctuates with partial independence from global unfolding events. Additionally, numerical fittings of the kinetic model of the process suggest that the local transition occurs in the seconds to minutes timescale. In fact, standard free energy differences for local unfolding were found to be significantly lower than those of the global unfolding reaction, showing that chemical modification results may not be explained in terms of the global unfolding reaction alone. These results provide unequivocal experimental evidence of local phenomena with global effects and contribute to understanding how global and local stability are linked to protein dynamics.

Simulated and experimental force spectroscopy of lysozyme on silica

The force spectra of proteins detaching from oxide surfaces measured by atomic force microscopy (AFM) often present complex patterns of peaks, which are difficult to correlate with individual bond-breaking events at the atomic scale. In this work we rationalize experimental AFM force spectra of hen-egg-white lysozyme detaching from silica by means of all-atom steered molecular dynamics (SMD) simulations. In particular, we demonstrate that the native tertiary structure of lysozyme is preserved if, and only if, its four intramolecular disulfide bridges are intact. Otherwise, the protein pulled off the surface undergoes severe unfolding, which is well captured by SMD simulations in explicit solvent. Implicit solvent simulations, on the contrary, wrongly predict protein unfolding even in the presence of S-S bridges, due to the lack of additional structural stabilization provided by the water's hydrogen-bond network within and surrounding the protein. On the basis of our combined experimental and theoretical findings, we infer that the rugged force spectra characteristic of lysozyme/silica interfaces are not due to the successive breaking of internal disulfide bonds leading to partial unfolding events. Rather, they reflect the detachment of several molecules bound to the same AFM tip, each anchored to the surface via multiple hydrogen and ionic bonds.

Spontaneous Unfolding-Refolding of Fibronectin Type III Domains Assayed by Thiol Exchange: THERMODYNAMIC STABILITY CORRELATES WITH RATES OF UNFOLDING RATHER THAN FOLDING

Globular proteins are not permanently folded but spontaneously unfold and refold on time scales that can span orders of magnitude for different proteins. A longstanding debate in the protein-folding field is whether unfolding rates or folding rates correlate to the stability of a protein. In the present study, we have determined the unfolding and folding kinetics of 10 FNIII domains. FNIII domains are one of the most common protein folds and are present in 2% of animal proteins. FNIII domains are ideal for this study because they have an identical seven-strand β-sandwich structure, but they vary widely in sequence and thermodynamic stability. We assayed thermodynamic stability of each domain by equilibrium denaturation in urea. We then assayed the kinetics of domain opening and closing by a technique known as thiol exchange. For this we introduced a buried Cys at the identical location in each FNIII domain and measured the kinetics of labeling with DTNB over a range of urea concentrations. A global fit of the kinetics data gave the kinetics of spontaneous unfolding and refolding in zero urea. We found that the folding rates were relatively similar, ∼0.1-1 s^-1, for the different domains. The unfolding rates varied widely and correlated with thermodynamic stability. Our study is the first to address this question using a set of domains that are structurally homologous but evolved with widely varying sequence identity and thermodynamic stability. These data add new evidence that thermodynamic stability correlates primarily with unfolding rate rather than folding rate. The study also has implications for the question of whether opening of FNIII domains contributes to the stretching of fibronectin matrix fibrils.

Structural and kinetic mapping of side-chain exposure onto the protein energy landscape

Identification and characterization of structural fluctuations that occur under native conditions is crucial for understanding protein folding and function, but such fluctuations are often rare and transient, making them difficult to study. Native-state hydrogen exchange (NSHX) has been a powerful tool for identifying such rarely populated conformations, but it generally reveals no information about the placement of these species along the folding reaction coordinate or the barriers separating them from the folded state and provides little insight into side-chain packing. To complement such studies, we have performed native-state alkyl-proton exchange, a method analogous to NSHX that monitors cysteine modification rather than backbone amide exchange, to examine the folding landscape of Escherichia coli ribonuclease H, a protein well characterized by hydrogen exchange. We have chosen experimental conditions such that the rate-limiting barrier acts as a kinetic partition: residues that become exposed only upon crossing the unfolding barrier are modified in the EX1 regime (alkylation rates report on the rate of unfolding), while those exposed on the native side of the barrier are modified predominantly in the EX2 regime (alkylation rates report on equilibrium populations). This kinetic partitioning allows for identification and placement of partially unfolded forms along the reaction coordinate. Using this approach we detect previously unidentified, rarely populated conformations residing on the native side of the barrier and identify side chains that are modified only upon crossing the unfolding barrier. Thus, in a single experiment under native conditions, both sides of the rate-limiting barrier are investigated.

A miniaturized technique for assessing protein thermodynamics and function using fast determination of quantitative cysteine reactivity

Protein thermodynamic stability is a fundamental physical characteristic that determines biological function. Furthermore, alteration of thermodynamic stability by macromolecular interactions or biochemical modifications is a powerful tool for assessing the relationship between protein structure, stability, and biological function. High-throughput approaches for quantifying protein stability are beginning to emerge that enable thermodynamic measurements on small amounts of material, in short periods of time, and using readily accessible instrumentation. Here we present such a method, fast quantitative cysteine reactivity, which exploits the linkage between protein stability, sidechain protection by protein structure, and structural dynamics to characterize the thermodynamic and kinetic properties of proteins. In this approach, the reaction of a protected cysteine and thiol-reactive fluorogenic indicator is monitored over a gradient of temperatures after a short incubation time. These labeling data can be used to determine the midpoint of thermal unfolding, measure the temperature dependence of protein stability, quantify ligand-binding affinity, and, under certain conditions, estimate folding rate constants. Here, we demonstrate the fQCR method by characterizing these thermodynamic and kinetic properties for variants of Staphylococcal nuclease and E. coli ribose-binding protein engineered to contain single, protected cysteines. These straightforward, information-rich experiments are likely to find applications in protein engineering and functional genomics.

Probing local structural fluctuations in myoglobin by size-dependent thiol-disulfide exchange

All proteins undergo local structural fluctuations (LSFs) or breathing motions. These motions are likely to be important for function but are poorly understood. LSFs were initially defined by amide hydrogen exchange (HX) experiments as opening events, which expose a small number of backbone amides to (1)H/(2)H exchange, but whose exchange rates are independent of denaturant concentration. Here, we use size-dependent thiol-disulfide exchange (SX) to characterize LSFs in single cysteine-containing variants of myoglobin (Mb). SX complements HX by providing information on motions that disrupt side chain packing interactions. Most importantly, probe reagents of different sizes and chemical properties can be used to characterize the size of structural opening events and the properties of the open state. We use thiosulfonate reagents (126-274 Da) to survey access to Cys residues, which are buried at specific helical packing interfaces in Mb. In each case, the free energy of opening increases linearly with the radius of gyration of the probe reagent. The slope and the intercept are interpreted to yield information on the size of the opening events that expose the buried thiol groups. The slope parameter varies by over 10-fold among Cys positions tested, suggesting that the sizes of breathing motions vary substantially throughout the protein. Our results provide insight to the longstanding question: how rigid or flexible are proteins in their native states?

Picomole-scale characterization of protein stability and function by quantitative cysteine reactivity

The Gibbs free energy difference between native and unfolded states ("stability") is one of the fundamental characteristics of a protein. By exploiting the thermodynamic linkage between ligand binding and stability, interactions of a protein with small molecules, nucleic acids, or other proteins can be detected and quantified. Determination of protein stability can therefore provide a universal monitor of biochemical function. Yet, the use of stability measurements as a functional probe is underutilized, because such experiments traditionally require large amounts of protein and special instrumentation. Here we present the quantitative cysteine reactivity (QCR) technique to determine protein stabilities rapidly and accurately using only picomole quantities of material and readily accessible laboratory equipment. We demonstrate that QCR-derived stabilities can be used to measure ligand binding over a wide range of ligand concentrations and affinities. We anticipate that this technique will have broad applications in high-throughput protein engineering experiments and functional genomics.

Electron-electron distances in spin-labeled low-spin metmyoglobin variants by relaxation enhancement

Thirteen single-cysteine variants of myoglobin were prepared by overexpression of apoprotein, spin labeling, and reconstitution with hemin. This procedure resulted in a protein with fewer hemichrome impurities than was obtained by an overexpression of holo-protein followed by spin labeling. Coordination of cyanide to the met heme formed low-spin complexes. Iron-nitroxyl interspin distances in the range of 17-30 A were determined by saturation recovery measurements of the enhancement of the nitroxyl spin lattice relaxation rates between approximately 30-140 K, and by spin-echo measurements of the enhancement of spin-spin relaxation rates at 10-30 K. Interspin distances were also calculated, using the molecular modeling program Insight II (Accelrys, San Diego, CA). For most variants, distances determined from the temperature dependence of spin-echo intensities at a pulse spacing of 200 ns agree with distances measured by saturation recovery and calculated with Insight II within about an angstrom, which is within experimental uncertainties. Measurements of interspin distances via spin-spin relaxation enhancement have the advantages that maximum effects are observed for slower metal relaxation rates than are required for spin-lattice relaxation enhancement, and the impact diminishes as r(-3) instead of r(-6), as with spin-lattice relaxation enhancement, which permits measurements at longer distances.

Expanding the repertoire of an ERK2 recruitment site: cysteine footprinting identifies the D-recruitment site as a mediator of Ets-1 binding

Many substrates of ERK2 contain a D-site, a sequence recognized by ERK2 that is used to promote catalysis. Despite lacking a canonical D-site, the substrate Ets-1 is displaced from ERK2 by peptides containing one. This suggests that Ets-1 may contain a novel or cryptic D-site. To investigate this possibility a protein footprinting strategy was developed to elucidate ERK2-ligand interactions. Using this approach, single cysteine reporters were placed in the D-recruitment site (DRS) of ERK2 and the resulting ERK2 proteins subjected to alkylation by iodoacetamide. The ability of residues 1-138 of Ets-1 to protect the cysteines from alkylation was determined. The pattern of protection observed is consistent with Ets-1 occupying a hydrophobic binding site within the DRS of ERK2. Significantly, a peptide derived from the D-site of Elk-1, which is known to bind the DRS, exhibits a similar pattern of cysteine protection. This analysis expands the repertoire of the DRS on ERK2 and suggests that other targeting sequences remain to be identified. Furthermore, cysteine-footprinting is presented as a useful way to interrogate protein-ligand interactions at the resolution of a single amino acid.

A Thiol Labelling Competition Experiment as a Probe for Sidechain Packing in the Kinetic Folding Intermediate of N-PGK

Protein folding is directed by the sequence of sidechains along the polypeptide backbone, but despite this the developement of sidechain interactions during folding is not well understood. Here, the thiol-active reagent, dithio-nitrobenzoic acid (DTNB), is used to probe the exposure of the cysteine sidechain thiols in the kinetic folding intermediates of the N-terminal domain of phosphoglycerate kinase (N-PGK) and a number of conservative (I-, L-, or V-to-C) single cysteine variants. Rapid dilution of chemically denatured protein into folding conditions in the presence of DTNB allowed the degree of sidechain protection in any rapidly formed intermediate to be determined through the analysis of the kinetics of labelling. The protection factors derived for the intermediate(s) were generally small (<25), indicating only partial burial of the sidechains. The distribution of protection parallels the previously reported backbone amide protection for the folding intermediate of N-PGK. These observations are consistent with the hypothesis that such intermediates resemble molten globule states; i.e. with native-like backbone hydrogen bonding and overall tertiary structure, but with the sidechains that make up the hydrophobic protein core dynamic and intermittently solvent exposed. The success of the competition technique in characterizing this kinetic intermediate invites application to other model systems.

How Cytochrome c Folds, and Why: Submolecular Foldon Units and their Stepwise Sequential Stabilization

Native state hydrogen exchange experiments have shown that the cytochrome c (Cyt c) protein consists of five cooperative folding-unfolding units, called foldons. These are named, in the order of increasing unfolding free energy, the nested-Yellow, Red, Yellow, Green, and Blue foldons. Previous results suggest that these units unfold in a stepwise sequential way so that each higher energy partially unfolded form includes all of the previously unfolded lower free energy units. If this is so, then selectively destabilizing any given foldon should equally destabilize each subsequent unfolding step above it in the unfolding ladder but leave the lower ones before it unaffected. To perform this test, we introduced the mutation Glu62Gly, which deletes a salt link in the Yellow unit and destabilizes the protein by 0.8 kcal/mol. Native state hydrogen exchange and other experiments show that the stability of the Yellow unit and the states above it in the free energy ladder are destabilized by about the same amount while the lower lying states are unaffected. These results help to confirm the sequential stepwise nature of the Cyt c unfolding pathway and therefore a similar refolding pathway. The steps in the pathway are dictated by the concerted folding-unfolding property of the individual unit foldons; the order of steps is determined by the sequential stabilization of progressively added foldons in the native context. Much related information for Cyt c strongly conforms with this mechanism. Its generality is supported by available information for other proteins.

Absence of kinetic thermal stabilization in a hyperthermophile rubredoxin indicated by 40 microsecond folding in the presence of irreversible denaturation

The striking kinetic stability of many proteins derived from hyperthermophilic organisms has led to the proposal that such stability may result from a heightened activation barrier for unfolding independent of a corresponding increase in the thermodynamic stability. This in turn implies a corresponding retardation of the folding reaction. A commonly cited model for kinetic thermal stabilization is the rubredoxin from Pyrococcus furiosus (Pf), which exhibits an irreversible denaturation lifetime at 100 degrees C of nearly a week. Utilizing protein resonances shifted well outside of the random coil chemical shift envelope, nuclear magnetic resonance (NMR) chemical exchange measurements on Pf rubredoxin as well as on the mesophile Clostridium pasteurianum (Cp) rubredoxin demonstrate reversible thermal transition temperatures of 144 degrees C (137 degrees C for the N-terminal modified A2K variant) and 104 degrees C, respectively, with similar (un)folding rates of approximately 25,000 s(-1), only modestly slower than the diffusion controlled rate. The absence of a substantial activation barrier to rubredoxin folding as well as the similar folding kinetics of the mesophile protein indicate that kinetic stabilization has not been utilized by the hyperthermophile rubredoxin in achieving its extreme thermal stability. The two-state folding kinetics observed for Pf rubredoxin contradict a previous assertion of multiphasic folding based on hydrogen exchange data extrapolated to an estimated midpoint of transition temperature (T(m)) of nearly 200 degrees C. This discrepancy is resolved by the observation that the base-catalyzed hydrogen exchange of the model dipeptide (N-acetyl-L-cysteine-N-methylamide)4-Cd2+ is 23-fold slower than that of the free cysteine model dipeptide used to normalize the Pf rubredoxin hydrogen exchange data.

The kinetics of side chain stabilization during protein folding

The rate of stabilization of side chains during protein folding has never been carefully studied. Recent developments in labeling proteins with (19)F-labeled amino acids coupled with real-time NMR measurements have allowed such measurements to be made. This paper describes the application of this method to the study of several proteins using 6-(19)F-tryptophan as the reporting group. It is found that these side chains adopt their final stable state at the last stages of the folding process and that the stabilization of side chains into their final conformation is a highly cooperative process. It is also possible to show the presence of intermediates in which the side chains are not correctly packed. The technique should be applicable to many systems.

Conformational dynamics of partially denatured myoglobin studied by time-resolved electrospray mass spectrometry with online hydrogen-deuterium exchange

This study demonstrates the use of electrospray mass spectrometry in conjunction with rapid online mixing ("time-resolved" ESI-MS) for monitoring protein conformational dynamics under equilibrium conditions. The hydrogen/deuterium exchange (HDX) kinetics of mildly denatured myoglobin (Mb) at pD 9.3, in the presence of 27% acetonitrile, were studied with millisecond time resolution. Analytical ultracentrifugation indicates that the average protein compactness under these solvent conditions is similar to that of native holomyoglobin (hMb). The mass spectrum shows protein ions in a wide array of charge and heme binding states, indicating the presence of multiple coexisting conformations. The experimental approach used allows the HDX kinetics of all of these species to be monitored separately. A combination of EX1 and EX2 behavior was observed for hMb ions in charge states 7+ to 9+, which predominantly represent nativelike hMb in solution. The EX1 kinetics are biphasic, indicating the presence of two protein populations that undergo conformational opening events with different rate constants. The EX2 kinetics observed for nativelike hMb are biphasic as well. All other charge and heme binding states represent non-native protein conformations that are involved in rapid interconversion processes, thus leading to monoexponential EX2 kinetics with a common rate constant. Burst phase labeling for these non-native proteins occurs at 125 sites. In contrast, the nativelike protein conformation shows burst phase labeling only for 88 sites. A kinetic model is developed which is based on the assumption of three distinct (un)folding units in Mb. The model implies that the free energy landscape of the protein exhibits a major barrier. The crossing of this barrier is most likely associated with slow, cooperative opening/closing events of the heme binding pocket. Rapid conformational fluctuations on either side of the barrier give rise to the observed EX2 kinetics. Simulated HDX kinetics based on this model are in excellent agreement with the experimental data.

Protein-folding kinetics and mechanisms studied by pulse-labeling and mass spectrometry

The "protein-folding problem" refers to the question of how and why a denatured polypeptide chain can spontaneously fold into a compact and highly ordered conformation. The classical description of this process in terms of reaction pathways has been complemented by models that describe folding as a biased conformational diffusion on a multidimensional energy landscape. The identification and characterization of short-lived intermediates provide important insights into the mechanism of folding. Pulsed hydrogen/deuterium exchange (HDX) methods are among the most powerful tools for studying the properties of kinetic intermediates. Analysis of pulse-labeled proteins by mass spectrometry (MS) provides information that is complementary to that obtained in nuclear magnetic resonance (NMR) studies; NMR data represent an average of entire protein ensembles, whereas MS can detect co-existing protein species. MS-based pulse-labeling experiments can distinguish between folding scenarios that involve parallel pathways, and those where folding is channeled through obligatory intermediates. The proteolytic digestion/MS technique provides spatially resolved information on the HDX pattern of folding intermediates. This method is especially important for proteins that are too large to be studied by NMR. Although traditional pulsed HDX protocols are based on quench-flow techniques, it is also possible to use electrospray (ESI) MS to analyze the reaction mixture on-line and "quasi-instantaneously" after labeling. This approach allows short-lived protein conformations to be studied by their HDX level, their ESI charge-state distribution, and their ligand-binding state. Covalent labeling of free cysteinyl residues provides an alternative approach to pulsed HDX experiments. Another promising development is the use of synchrotron X-rays to induce oxidation at specific sites within a protein for studying their solvent accessibility during folding.

The equilibrium unfolding pathway of a (beta/alpha)8 barrel

The (beta/alpha)(8) barrel is the most commonly occurring fold among enzymes. A key step towards rationally engineering (beta/alpha)(8) barrel proteins is to understand their underlying structural organization and folding energetics. Using misincorporation proton-alkyl exchange (MPAX), a new tool for solution structural studies of large proteins, we have performed a native-state exchange analysis of the prototypical (beta/alpha)(8) barrel triosephosphate isomerase. Three cooperatively unfolding subdomains within the structure are identified, as well as two partially unfolded forms of the protein. The C-terminal domain coincides with domains reported to exist in four other (beta/alpha)(8) barrels, but the two N-terminal domains have not been observed previously. These partially unfolded forms may represent sequential intermediates on the folding pathway of triosephosphate isomerase. The methods reported here should be applicable to a variety of other biological problems involving protein conformational changes.

Structural characterization of non-native states of sperm whale myoglobin in aqueous ethanol or 2,2,2-trifluoroethanol media

The effects of aqueous ethanol or 2,2,2-trifluoroethanol media on the structure of sperm whale myoglobin have been investigated by absorption, CD, and NMR spectra. The structural properties of myoglobin such as heme environments, helix contents, protein folding, and interactions between heme and the protein moiety have been sharply manifested in these spectra. The characterization demonstrated that alcohol-induced conformational change of myoglobin depends on the nature of alcohol and its concentration. It was shown for the first time that, upon the alcohol-induced denaturation of myoglobin, heme is released from partially denatured protein of which helix contents is altered by only about 20% relative to that of native state. Myoglobin has shown to unfold and refold reversibly by controlling the alcohol concentration. Novel methods for the preparation of apomyoglobin and in situ reconstitution of apomyoglobin with heme, based on the alcohol-induced denaturation of the protein, were presented.