Local dynamics and stability of apocytochrome b562 examined by hydrogen exchange

Exploring Residue-Level Interactions between the Biofilm-Driving R2ab Protein and Polystyrene Nanoparticles

In biological systems, proteins can bind to nanoparticles to form a "corona" of adsorbed molecules. The nanoparticle corona is of significant interest because it impacts an organism's response to a nanomaterial. Understanding the corona requires knowledge of protein structure, orientation, and dynamics at the surface. A residue-level mapping of protein behavior on nanoparticle surfaces is needed, but this mapping is difficult to obtain with traditional approaches. Here, we have investigated the interaction between R2ab and polystyrene nanoparticles (PSNPs) at the level of individual residues. R2ab is a bacterial surface protein from Staphylococcus epidermidis and is known to interact strongly with polystyrene, leading to biofilm formation. We have used mass spectrometry after lysine methylation and hydrogen-deuterium exchange (HDX) NMR spectroscopy to understand how the R2ab protein interacts with PSNPs of different sizes. Lysine methylation experiments reveal subtle but statistically significant changes in methylation patterns in the presence of PSNPs, indicating altered protein surface accessibility. HDX rates become slower overall in the presence of PSNPs. However, some regions of the R2ab protein exhibit faster than average exchange rates in the presence of PSNPs, while others are slower than the average behavior, suggesting conformational changes upon binding. HDX rates and methylation ratios support a recently proposed "adsorbotope" model for PSNPs, wherein adsorbed proteins consist of unfolded anchor points interspersed with partially structured regions. Our data also highlight the challenges of characterizing complex protein-nanoparticle interactions using these techniques, such as fast exchange rates. While providing insights into how R2ab adsorbs onto PSNP surfaces, this research emphasizes the need for advanced methods to comprehend residue-level interactions in the nanoparticle corona.

Exploring the Residue-Level Interactions between the R2ab Protein and Polystyrene Nanoparticles

In biological systems, proteins can bind to nanoparticles to form a "corona" of adsorbed molecules. The nanoparticle corona is of high interest because it impacts the organism's response to the nanomaterial. Understanding the corona requires knowledge of protein structure, orientation, and dynamics at the surface. Ultimately, a residue-level mapping of protein behavior on nanoparticle surfaces is needed, but this mapping is difficult to obtain with traditional approaches. Here, we have investigated the interaction between R2ab and polystyrene nanoparticles (PSNPs) at the level of individual residues. R2ab is a bacterial surface protein from Staphylococcus epidermidis and is known to interact strongly with polystyrene, leading to biofilm formation. We have used mass spectrometry after lysine methylation and hydrogen-deuterium exchange (HDX) NMR spectroscopy to understand how the R2ab protein interacts with PSNPs of different sizes. Through lysine methylation, we observe subtle but statistically significant changes in methylation patterns in the presence of PSNPs, indicating altered protein surface accessibility. HDX measurements reveal that certain regions of the R2ab protein undergo faster exchange rates in the presence of PSNPs, suggesting conformational changes upon binding. Both results support a recently proposed "adsorbotope" model, wherein adsorbed proteins consist of unfolded anchor points interspersed with regions of partial structure. Our data also highlight the challenges of characterizing complex protein-nanoparticle interactions using these techniques, such as fast exchange rates. While providing insights into how proteins respond to nanoparticle surfaces, this research emphasizes the need for advanced methods to comprehend these intricate interactions fully at the residue level.

Cryptic intermediates and metastable states of proteins as predicted by OneG computational method

Understanding structural excursions of proteins under folding conditions is crucial to map energy landscapes of proteins. In the present study, OneG computational tool has been used for analyzing possible existence of cryptic intermediates and metastable states of 26 proteins for which three prerequisite inputs of the OneG such as atomic coordinates of proteins, free energy of unfolding (ΔG_U) and free energy of exchange (ΔG_HX) determined in the absence of denaturant were available during the course of the study. The veraciousness of the tool on predicting the partially folded states of the proteins has been comprehensively described using experimental data available for 15 of the 26 proteins. Meanwhile, possible existence of partially structured states in the folding pathways of 11 other proteins has also been delineated as predicted by the OneG. In addition to mapping the folding pathways of proteins, the salient merits of the tool on systematically addressing the discrepancy between the ΔG_U and the ΔG_HX of the proteins have also been dealt.Communicated by Ramaswamy H. Sarma.

Wilson D. HDX-MS: An Analytical Tool to Capture Protein Motion in Action

Virtually all protein functions in the cell, including pathogenic processes, require coordinated motion of atoms or domains, i.e., conformational dynamics. Understanding protein dynamics is therefore critical both for drug development and to learn about the underlying molecular causes of many diseases. Hydrogen–Deuterium Exchange Mass Spectrometry (HDX-MS) provides valuable information about protein dynamics, which is highly complementary to the static picture provided by conventional high-resolution structural tools (i.e., X-ray crystallography and structural NMR). The amount of protein required to carry out HDX-MS experiments is a fraction of the amount required by alternative biophysical techniques, which are also usually lower resolution. Use of HDX-MS is growing quickly both in industry and academia, and it has been successfully used in numerous drug and vaccine development efforts, with important roles in understanding allosteric effects and mapping binding sites.

Cytochrome c folds through foldon-dependent native-like intermediates in an ordered pathway

Previous hydrogen exchange (HX) studies of the spontaneous reversible unfolding of Cytochrome c (Cyt c) under native conditions have led to the following conclusions. Native Cyt c (104 residues) is composed of five cooperative folding units, called foldons. The high-energy landscape is dominated by an energy ladder of partially folded forms that differ from each other by one cooperative foldon unit. The reversible equilibrium unfolding of native Cyt c steps up through these intermediate forms to the unfolded state in an energy-ordered sequence, one foldon unit at a time. To more directly study Cyt c intermediates and pathways during normal energetically downhill kinetic folding, the present work used HX pulse labeling analyzed by a fragment separation-mass spectrometry method. The results show that 95% or more of the Cyt c population folds by stepping down through the same set of foldon-dependent pathway intermediates as in energetically uphill equilibrium unfolding. These results add to growing evidence that proteins fold through a classical pathway sequence of native-like intermediates rather than through a vast number of undefinable intermediates and pathways. The present results also emphasize the condition-dependent nature of kinetic barriers, which, with less informative experimental methods (fluorescence, etc.), are often confused with variability in intermediates and pathways.

Domain-swapped cytochrome cb562 dimer and its nanocage encapsulating a Zn-SO4 cluster in the internal cavity

Three domain-swapped cytochrome cb₅₆₂ dimers formed a unique cage structure with a Zn–SO₄ cluster inside the cavity.

Protein nanostructures have been gaining in interest, along with developments in new methods for construction of novel nanostructures. We have previously shown that c-type cytochromes and myoglobin form oligomers by domain swapping. Herein, we show that a four-helix bundle protein cyt cb₅₆₂, with the cyt b₅₆₂ heme attached to the protein moiety by two Cys residues insertion, forms a domain-swapped dimer. Dimeric cyt cb₅₆₂ did not dissociate to monomers at 4 °C, whereas dimeric cyt b₅₆₂ dissociated under the same conditions, showing that heme attachment to the protein moiety stabilizes the domain-swapped structure. According to X-ray crystallographic analysis of dimeric cyt cb₅₆₂, the two helices in the N-terminal region of one protomer interacted with the other two helices in the C-terminal region of the other protomer, where Lys51–Asp54 served as a hinge loop. The heme coordination structure of the dimer was similar to that of the monomer. In the crystal, three domain-swapped cyt cb₅₆₂ dimers formed a unique cage structure with a Zn–SO₄ cluster inside the cavity. The Zn–SO₄ cluster consisted of fifteen Zn²⁺ and seven SO₄^2– ions, whereas six additional Zn²⁺ ions were detected inside the cavity. The cage structure was stabilized by coordination of the amino acid side chains of the dimers to the Zn²⁺ ions and connection of two four-helix bundle units through the conformation-adjustable hinge loop. These results show that domain swapping can be applied in the construction of unique protein nanostructures.

Visualizing transient dark states by NMR spectroscopy

Myriad biological processes proceed through states that defy characterization by conventional atomic-resolution structural biological methods. The invisibility of these 'dark' states can arise from their transient nature, low equilibrium population, large molecular weight, and/or heterogeneity. Although they are invisible, these dark states underlie a range of processes, acting as encounter complexes between proteins and as intermediates in protein folding and aggregation. New methods have made these states accessible to high-resolution analysis by nuclear magnetic resonance (NMR) spectroscopy, as long as the dark state is in dynamic equilibrium with an NMR-visible species. These methods - paramagnetic NMR, relaxation dispersion, saturation transfer, lifetime line broadening, and hydrogen exchange - allow the exploration of otherwise invisible states in exchange with a visible species over a range of timescales, each taking advantage of some unique property of the dark state to amplify its effect on a particular NMR observable. In this review, we introduce these methods and explore two specific techniques - paramagnetic relaxation enhancement and dark state exchange saturation transfer - in greater detail.

OneG-Vali: a computational tool for detecting, estimating and validating cryptic intermediates of proteins under native conditions

Unfolding pathway of T4 lysozyme under native conditions as predicted by the OneG-Vali has been illustrated. Also, structural contexts of various states (native (N), cryptic intermediates (CIs) and unfolded (U) conformations) of the protein and the population of three CIs are depicted.

Using the COREX/BEST server to model the native-state ensemble

Protein structures under normal conditions exist as ensembles of interconverting, transient microstates. A computer algorithm known as COREX/BEST (Biology using Ensemble-based Structural Thermodynamics) was developed to model microstate structures and describe the native ensembles of proteins in statistical thermodynamic terms. This algorithm has been tested extensively and validated through experimental comparisons examining a range of biophysical and functional phenomena, such as structural cooperativity, pH-dependent stability, and cold denaturation. Here, we describe a Web-based implementation of the COREX/BEST algorithm, called the COREX/BEST Server, and demonstrate how to use this online resource to characterize the structural and thermodynamic properties of the native protein ensemble.

Folding of a large protein at high structural resolution

Kinetic folding of the large two-domain maltose binding protein (MBP; 370 residues) was studied at high structural resolution by an advanced hydrogen-exchange pulse-labeling mass-spectrometry method (HX MS). Dilution into folding conditions initiates a fast molecular collapse into a polyglobular conformation (<20 ms), determined by various methods including small angle X-ray scattering. The compaction produces a structurally heterogeneous state with widespread low-level HX protection and spectroscopic signals that match the equilibrium melting posttransition-state baseline. In a much slower step (7-s time constant), all of the MBP molecules, although initially heterogeneously structured, form the same distinct helix plus sheet folding intermediate with the same time constant. The intermediate is composed of segments that are distant in the MBP sequence but adjacent in the native protein where they close the longest residue-to-residue contact. Segments that are most HX protected in the early molecular collapse do not contribute to the initial intermediate, whereas the segments that do participate are among the less protected. The 7-s intermediate persists through the rest of the folding process. It contains the sites of three previously reported destabilizing mutations that greatly slow folding. These results indicate that the intermediate is an obligatory step on the MBP folding pathway. MBP then folds to the native state on a longer time scale (~100 s), suggestively in more than one step, the first of which forms structure adjacent to the 7-s intermediate. These results add a large protein to the list of proteins known to fold through distinct native-like intermediates in distinct pathways.

Stepwise protein folding at near amino acid resolution by hydrogen exchange and mass spectrometry

The kinetic folding of ribonuclease H was studied by hydrogen exchange (HX) pulse labeling with analysis by an advanced fragment separation mass spectrometry technology. The results show that folding proceeds through distinct intermediates in a stepwise pathway that sequentially incorporates cooperative native-like structural elements to build the native protein. Each step is seen as a concerted transition of one or more segments from an HX-unprotected to an HX-protected state. Deconvolution of the data to near amino acid resolution shows that each step corresponds to the folding of a secondary structural element of the native protein, termed a "foldon." Each folded segment is retained through subsequent steps of foldon addition, revealing a stepwise buildup of the native structure via a single dominant pathway. Analysis of the pertinent literature suggests that this model is consistent with experimental results for many proteins and some current theoretical results. Two biophysical principles appear to dictate this behavior. The principle of cooperativity determines the central role of native-like foldon units. An interaction principle termed "sequential stabilization" based on native-like interfoldon interactions orders the pathway.

Protein dynamics viewed by hydrogen exchange

To examine the relationship between protein structural dynamics and measurable hydrogen exchange (HX) data, the detailed exchange behavior of most of the backbone amide hydrogens of Staphylococcal nuclease was compared with that of their neighbors, with their structural environment, and with other information. Results show that H-bonded hydrogens are protected from exchange, with HX rate effectively zero, even when they are directly adjacent to solvent. The transition to exchange competence requires a dynamic structural excursion that removes H-bond protection and allows exposure to solvent HX catalyst. The detailed data often make clear the nature of the dynamic excursion required. These range from whole molecule unfolding, through smaller cooperative unfolding reactions of secondary structural elements, and down to local fluctuations that involve as little as a single peptide group or side chain or water molecule. The particular motion that dominates the exchange of any hydrogen is the one that allows the fastest HX rate. The motion and the rate it produces are determined by surrounding structure and not by nearness to solvent or the strength of the protecting H-bond itself or its acceptor type (main chain, side chain, structurally bound water). Many of these motions occur over time scales that are appropriate for biochemical function.

OneG: a computational tool for predicting cryptic intermediates in the unfolding kinetics of proteins under native conditions

Understanding the relationships between conformations of proteins and their stabilities is one key to address the protein folding paradigm. The free energy change (ΔG) of unfolding reactions of proteins is measured by traditional denaturation methods and native hydrogen-deuterium (H/D) exchange methods. However, the free energy of unfolding (ΔG(U)) and the free energy of exchange (ΔG(HX)) of proteins are not in good agreement, though the experimental conditions of both methods are well matching to each other. The anomaly is due to any one or combinations of the following reasons: (i) effects of cis-trans proline isomerisation under equilibrium unfolding reactions of proteins (ii) inappropriateness in accounting the baselines of melting curves (iii) presence of cryptic intermediates, which may elude the melting curve analysis and (iv) existence of higher energy metastable states in the H/D exchange reactions of proteins. Herein, we have developed a novel computational tool, OneG, which accounts the discrepancy between ΔG(U) and ΔG(HX) of proteins by systematically accounting all the four factors mentioned above. The program is fully automated and requires four inputs: three-dimensional structures of proteins, ΔG(U), ΔG(U)(*) and residue-specific ΔG(HX) determined under EX2-exchange conditions in the absence of denaturants. The robustness of the program has been validated using experimental data available for proteins such as cytochrome c and apocytochrome b(562) and the data analyses revealed that cryptic intermediates of the proteins detected by the experimental methods and the cryptic intermediates predicted by the OneG for those proteins were in good agreement. Furthermore, using OneG, we have shown possible existence of cryptic intermediates and metastable states in the unfolding pathways of cardiotoxin III and cobrotoxin, respectively, which are homologous proteins. The unique application of the program to map the unfolding pathways of proteins under native conditions have been brought into fore and the program is publicly available at http://sblab.sastra.edu/oneg.html.

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.

What lessons can be learned from studying the folding of homologous proteins?

The studies of the folding of structurally related proteins have proved to be a very important tool for investigating protein folding. Here we review some of the insights that have been gained from such studies. Our highlighted studies show just how such an investigation should be designed and emphasise the importance of the synergy between experiment and theory. We also stress the importance of choosing the right system carefully, exploiting the excellent structural and sequence databases at our disposal.

AMORE-HX: a multidimensional optimization of radial enhanced NMR-sampled hydrogen exchange

The Cartesian sampled three-dimensional HNCO experiment is inherently limited in time resolution and sensitivity for the real time measurement of protein hydrogen exchange. This is largely overcome by use of the radial HNCO experiment that employs the use of optimized sampling angles. The significant practical limitation presented by use of three-dimensional data is the large data storage and processing requirements necessary and is largely overcome by taking advantage of the inherent capabilities of the 2D-FT to process selective frequency space without artifact or limitation. Decomposition of angle spectra into positive and negative ridge components provides increased resolution and allows statistical averaging of intensity and therefore increased precision. Strategies for averaging ridge cross sections within and between angle spectra are developed to allow further statistical approaches for increasing the precision of measured hydrogen occupancy. Intensity artifacts potentially introduced by over-pulsing are effectively eliminated by use of the BEST approach.

Compressing the free energy range of substructure stabilities in iso-1-cytochrome c

Evolutionary conservation of substructure architecture between yeast iso-1-cytochrome c and the well-characterized horse cytochrome c is studied with limited proteolysis, the alkaline conformational transition and global unfolding with guanidine-HCl. Mass spectral analysis of limited proteolysis cleavage products for iso-1-cytochrome c show that its least stable substructure is the same as horse cytochrome c. The limited proteolysis data yield a free energy of 3.8 +/- 0.4 kcal mol(-1) to unfold the least stable substructure compared with 5.05 +/- 0.30 kcal mol(-1) for global unfolding of iso-1-cytochrome c. Thus, substructure stabilities of iso-1-cytochrome c span only approximately 1.2 kcal mol(-1) compared with approximately 8 kcal mol(-1) for horse cytochrome c. Consistent with the less cooperative folding thus expected for the horse protein, the guanidine-HCl m-values are approximately 3 kcal mol(-1)M(-1) versus approximately 4.5 kcal mol(-1)M(-1) for horse versus yeast cytochrome c. The tight free energy spacing of the yeast cytochrome c substructures suggests that its folding has more branch points than for horse cytochrome c. Studies on a variant of iso-1-cytochrome c with an H26N mutation indicate that the least and most stable substructures unfold sequentially and the two least stable substructures unfold independently as for horse cytochrome c. Thus, important aspects of the substructure architecture of horse cytochrome c, albeit compressed energetically, are preserved evolutionally in yeast iso-1-cytochrome c.

Reverse micelles in integral membrane protein structural biology by solution NMR spectroscopy

Integral membrane proteins remain a significant challenge to structural studies by solution NMR spectroscopy. This is due not only to spectral complexity, but also because the effects of slow molecular reorientation are exacerbated by the need to solubilize the protein in aqueous detergent micelles. These assemblies can be quite large and require deuteration for optimal use of the TROSY effect. In principle, another approach is to employ reverse micelle encapsulation to solubilize the protein in a low-viscosity solvent in which the rapid tumbling of the resulting particle allows for use of standard triple-resonance methods. The preparation of such samples of membrane proteins is difficult. Using a 54 kDa construct of the homotetrameric potassium channel KcsA, we demonstrate a strategy that employs a hybrid surfactant to transfer the protein to the reverse micelle system.

The foldon substructure of staphylococcal nuclease

To search for submolecular foldon units, the spontaneous reversible unfolding and refolding of staphylococcal nuclease under native conditions was studied by a kinetic native-state hydrogen exchange (HX) method. As for other proteins, it appears that staphylococcal nuclease is designed as an assembly of well-integrated foldon units that may define steps in its folding pathway and may regulate some other functional properties. The HX results identify 34 amide hydrogens that exchange with solvent hydrogens under native conditions by way of large transient unfolding reactions. The HX data for each hydrogen measure the equilibrium stability (Delta G(HX)) and the kinetic unfolding and refolding rates (k(op) and k(cl)) of the unfolding reaction that exposes it to exchange. These parameters separate the 34 identified residues into three distinct HX groupings. Two correspond to clearly defined structural units in the native protein, termed the blue and red foldons. The remaining HX grouping contains residues, not well separated by their HX parameters alone, that represent two other distinct structural units in the native protein, termed the green and yellow foldons. Among these four sets, a last unfolding foldon (blue) unfolds with a rate constant of 6 x 10(-6) s(-1) and free energy equal to the protein's global stability (10.0 kcal/mol). It represents part of the beta-barrel, including mutually H-bonding residues in the beta 4 and beta 5 strands, a part of the beta 3 strand that H-bonds to beta 5, and residues at the N-terminus of the alpha2 helix that is capped by beta 5. A second foldon (green), which unfolds and refolds more rapidly and at slightly lower free energy, includes residues that define the rest of the native alpha2 helix and its C-terminal cap. A third foldon (yellow) defines the mutually H-bonded beta1-beta2-beta 3 meander, completing the native beta-barrel, plus an adjacent part of the alpha1 helix. A final foldon (red) includes residues on remaining segments that are distant in sequence but nearly adjacent in the native protein. Although the structure of the partially unfolded forms closely mimics the native organization, four residues indicate the presence of some nonnative misfolding interactions. Because the unfolding parameters of many other residues are not determined, it seems likely that the concerted foldon units are more extensive than is shown by the 34 residues actually observed.

Protein folding and misfolding: mechanism and principles

Two fundamentally different views of how proteins fold are now being debated. Do proteins fold through multiple unpredictable routes directed only by the energetically downhill nature of the folding landscape or do they fold through specific intermediates in a defined pathway that systematically puts predetermined pieces of the target native protein into place? It has now become possible to determine the structure of protein folding intermediates, evaluate their equilibrium and kinetic parameters, and establish their pathway relationships. Results obtained for many proteins have serendipitously revealed a new dimension of protein structure. Cooperative structural units of the native protein, called foldons, unfold and refold repeatedly even under native conditions. Much evidence obtained by hydrogen exchange and other methods now indicates that cooperative foldon units and not individual amino acids account for the unit steps in protein folding pathways. The formation of foldons and their ordered pathway assembly systematically puts native-like foldon building blocks into place, guided by a sequential stabilization mechanism in which prior native-like structure templates the formation of incoming foldons with complementary structure. Thus the same propensities and interactions that specify the final native state, encoded in the amino-acid sequence of every protein, determine the pathway for getting there. Experimental observations that have been interpreted differently, in terms of multiple independent pathways, appear to be due to chance misfolding errors that cause different population fractions to block at different pathway points, populate different pathway intermediates, and fold at different rates. This paper summarizes the experimental basis for these three determining principles and their consequences. Cooperative native-like foldon units and the sequential stabilization process together generate predetermined stepwise pathways. Optional misfolding errors are responsible for 3-state and heterogeneous kinetic folding.

Branching in the sequential folding pathway of cytochrome c

Previous results indicate that the folding pathways of cytochrome c and other proteins progressively build the target native protein in a predetermined stepwise manner by the sequential formation and association of native-like foldon units. The present work used native state hydrogen exchange methods to investigate a structural anomaly in cytochrome c results that suggested the concerted folding of two segments that have little structural relationship in the native protein. The results show that the two segments, an 18-residue omega loop and a 10-residue helix, are able to unfold and refold independently, which allows a branch point in the folding pathway. The pathway that emerges assembles native-like foldon units in a linear sequential manner when prior native-like structure can template a single subsequent foldon, and optional pathway branching is seen when prior structure is able to support the folding of two different foldons.

A unified mechanism for protein folding: predetermined pathways with optional errors

There is a fundamental conflict between two different views of how proteins fold. Kinetic experiments and theoretical calculations are often interpreted in terms of different population fractions folding through different intermediates in independent unrelated pathways (IUP model). However, detailed structural information indicates that all of the protein population folds through a sequence of intermediates predetermined by the foldon substructure of the target protein and a sequential stabilization principle. These contrary views can be resolved by a predetermined pathway--optional error (PPOE) hypothesis. The hypothesis is that any pathway intermediate can incorporate a chance misfolding error that blocks folding and must be reversed for productive folding to continue. Different fractions of the protein population will then block at different steps, populate different intermediates, and fold at different rates, giving the appearance of multiple unrelated pathways. A test of the hypothesis matches the two models against extensive kinetic folding results for hen lysozyme which have been widely cited in support of independent parallel pathways. The PPOE model succeeds with fewer fitting constants. The fitted PPOE reaction scheme leads to known folding behavior, whereas the IUP properties are contradicted by experiment. The appearance of a conflict with multipath theoretical models seems to be due to their different focus, namely on multitrack microscopic behavior versus cooperative macroscopic behavior. The integration of three well-documented principles in the PPOE model (cooperative foldons, sequential stabilization, optional errors) provides a unifying explanation for how proteins fold and why they fold in that way.

Native state energetics of the Src SH2 domain: evidence for a partially structured state in the denatured ensemble

We have defined the free-energy profile of the Src SH2 domain using a variety of biophysical techniques. Equilibrium and kinetic experiments monitored by tryptophan fluorescence show that Src SH2 is quite stable and folds rapidly by a two-state mechanism, without populating any intermediates. Native state hydrogen-deuterium exchange confirms this two-state behavior; we detect no cooperative partially unfolded forms in equilibrium with the native conformation under any conditions. Interestingly, the apparent stability of the protein from hydrogen exchange is 2 kcal/mol greater than the stability determined by both equilibrium and kinetic studies followed by fluorescence. Native-state proteolysis demonstrates that this "super protection" does not result from a deviation from the linear extrapolation model used to fit the fluorescence data. Instead, it likely arises from a notable compaction in the unfolded state under native conditions, resulting in an ensemble of conformations with substantial solvent exposure of side chains and flexible regions sensitive to proteolysis, but backbone amides that exchange with solvent approximately 30-fold slower than would be expected for a random coil. The apparently simple behavior of Src SH2 in traditional unfolding studies masks the significant complexity present in the denatured-state ensemble.

Order of steps in the cytochrome C folding pathway: evidence for a sequential stabilization mechanism

Previous work used hydrogen exchange (HX) experiments in kinetic and equilibrium modes to study the reversible unfolding and refolding of cytochrome c (Cyt c) under native conditions. Accumulated results now show that Cyt c is composed of five individually cooperative folding units, called foldons, which unfold and refold as concerted units in a stepwise pathway sequence. The first three steps of the folding pathway are linear and sequential. The ordering of the last two steps has been unclear because the fast HX of the amino acid residues in these foldons has made measurement difficult. New HX experiments done under slower exchange conditions show that the final two foldons do not unfold and refold in an obligatory sequence. They unfold separately and neither unfolding obligately contains the other, as indicated by their similar unfolding surface exposure and the specific effects of destabilizing and stabilizing mutations, pH change, and oxidation state. These results taken together support a sequential stabilization mechanism in which folding occurs in the native context with prior native-like structure serving to template the stepwise formation of subsequent native-like foldon units. Where the native structure of Cyt c requires sequential folding, in the first three steps, this is found. Where structural determination is ambiguous, in the final two steps, alternative parallel folding is found.

GuHCl and NaCl-dependent hydrogen exchange in MerP reveals a well-defined core with an unusual exchange pattern

We have analysed hydrogen exchange at amide groups to characterise the energy landscape of the 72 amino acid residue protein MerP. From the guanidine hydrochloride (GuHCl) dependence of exchange in the pre-transitional region we have determined free energy values of exchange (DeltaG(HX)) and corresponding m-values for individual amide protons. Detailed analysis of the exchange patterns indicates that for one set of amide protons there is a weak dependence on denaturant, indicating that the exchange is dominated by local fluctuations. For another set of amide protons a linear, but much stronger, denaturant dependence is observed. Notably, the plots of free energy of exchange versus [GuHCl] for 16 amide protons show pronounced upward curvature, and a close inspection of the structure shows that these residues form a well-defined core in the protein. The hydrogen exchange that was measured at various concentrations of NaCl shows an apparent selective stabilisation of this core. Detailed analysis of this exchange pattern indicates that it may originate from selective destabilisation of the unfolded state by guanidinium ions and/or selective stabilisation of the core in the native state by chloride ions.

Functional role of a protein foldon-An Ω-loop foldon controls the alkaline transition in ferricytochrome c

Hydrogen exchange results for cytochrome c and several other proteins show that they are composed of a number of foldon units which continually unfold and refold and account for some functional properties. Previous work showed that one Omega-loop foldon controls the rate of the structural switching and ligand exchange behavior of cytochrome c known as the alkaline transition. The present work tests the role of foldons in the alkaline transition equilibrium. We measured the effects of denaturant and 14 destabilizing mutations. The results show that the ligand exchange equilibrium is controlled by the stability of the same foldon unit implicated before. In addition, the results obtained confirm the epsilon-amino group of Lys79 and Lys73 as the alkaline replacement ligands and bear on the search for a triggering group.

Local conformational fluctuations can modulate the coupling between proton binding and global structural transitions in proteins

Local conformational fluctuations in proteins can affect the coupling between ligand binding and global structural transitions. This finding was established by monitoring quantitatively how the population distribution in the ensemble of microstates of staphylococcal nuclease was affected by proton binding. Analysis of acid unfolding and proton-binding data with an ensemble-based model suggests that local fluctuations: (i) can be effective modulators of ligand-binding affinities, (ii) are important determinants of the cooperativity of ligand-driven global structural transitions, and (iii) are well represented thermodynamically as local unfolding processes. These studies illustrate how an ensemble-based description of proteins can be used to describe quantitatively the interdependence of local conformational fluctuations, ligand-binding processes, and global structural transitions. This level of understanding of the relationship between conformation, energy, and dynamics is required for a detailed mechanistic understanding of allostery, cooperativity, and other complex functional and regulatory properties of macromolecules.

Protein folding: the stepwise assembly of foldon units

Equilibrium and kinetic hydrogen exchange experiments show that cytochrome c is composed of five foldon units that continually unfold and refold even under native conditions. Folding proceeds by the stepwise assembly of the foldon units rather than one amino acid at a time. The folding pathway is determined by a sequential stabilization process; previously formed foldons guide and stabilize subsequent foldons to progressively build the native protein. Four other proteins have been found to show similar behavior. These results support stepwise protein folding pathways through discrete intermediates.

The N-terminal to C-terminal motif in protein folding and function

Essentially all proteins known to fold kinetically in a two-state manner have their N- and C-terminal secondary structural elements in contact, and the terminal elements often dock as part of the experimentally measurable initial folding step. Conversely, all N-C no-contact proteins studied so far fold by non-two-state kinetics. By comparison, about half of the single domain proteins in the Protein Data Bank have their N- and C-terminal elements in contact, more than expected on a random probability basis but not nearly enough to account for the bias in protein folding. Possible reasons for this bias relate to the mechanisms for initial protein folding, native state stability, and final turnover.

Protein misfolding: optional barriers, misfolded intermediates, and pathway heterogeneity

To investigate the character and role of misfolded intermediates in protein folding, a recombinant cytochrome c without the normally blocking histidine to heme misligation was studied. Folding remains heterogeneous as in the wild-type protein. Half of the population folds relatively rapidly to the native state in a two-state manner. The other half collapses (fluorescence quenching) and forms a full complement of helix (CD) with the same rate and denaturant dependence as the fast folding fraction but then is blocked and reaches the native structure (695nm absorbance) much more slowly. The factors that transiently block folding are not intrinsic to the folding process but depend on ambient conditions, including protein aggregation (f(concentration)), N terminus to heme misligation (f(pH)), and proline mis-isomerization (f(U state equilibration time)). The misfolded intermediate populated by the slowly folding fraction was characterized by hydrogen exchange pulse labeling. It is very advanced with all of the native-like elements fairly stably formed but not the final Met80-S to heme iron ligation, similar to a previously studied molten globule form induced by low pH. To complete final native state acquisition, some small back unfolding is required (error repair) but the misfolded intermediate does not revisit the U state before proceeding to N. These properties show that the intermediate is a normal on-pathway form that contains, in addition, adventitious misfolding errors that transiently block its forward progress. Related observations for other proteins (partially misfolded intermediates, pathway heterogeneity) might be similarly explained in terms of the optional insertion of error-dependent barriers into a classical folding pathway.

New Insights into FAK Signaling and Localization Based on Detection of a FAT Domain Folding Intermediate

Mounting evidence suggests that the focal adhesion targeting (FAT) domain, an antiparallel four-helix bundle, exists in alternative conformations that may modulate phosphorylation, ligand binding, and the subcellular localization of focal adhesion kinase (FAK). In order to characterize the conformational dynamics of the FAT domain, we have developed a novel method for reconstructing the folding pathway of the FAT domain by using discrete molecular dynamics (DMD) simulations, with free energy constraints derived from NMR hydrogen exchange data. The DMD simulations detect a folding intermediate, in which a cooperative unfolding event causes helix 1 to lose helical character while separating from the helix bundle. The conformational dynamic features of helix 1 in the intermediate state of the FAT domain are likely to facilitate Y926 phosphorylation, yet interfere with paxillin binding. The presence of this intermediate state in vivo may promote FAK signaling via the ERK/MAPK pathway and by release of FAK from focal adhesions.

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.

Proton dependence of tobacco mosaic virus dissociation by pressure

Tobacco mosaic virus (TMV) is an intensely studied model of viruses. This paper reports an investigation into the dissociation of TMV by pH and pressure up to 220 MPa. The viral solution (0.25 mg/ml) incubated at 277 K showed a significant decrease in light scattering with increasing pH, suggesting dissociation. This observation was confirmed by HPLC gel filtration and electron microscopy. The calculated volume change of dissociation (DeltaV) decreased (absolute value) from -49.7 ml/mol of subunit at pH 3.8 to -21.7 ml/mol of subunit at pH 9.0. The decrease from pH 9.0 to 3.8 caused a stabilization of 14.1 kJ/mol of TMV subunit. The estimated proton release calculated from pressure-induced dissociation curves was 0.584 mol H(+)/mol of TMV subunit. These results suggest that the degree of virus inactivation by pressure and the immunogenicity of the inactivated structures can be optimized by modulating the surrounding pH.

Beyond the EX1 limit: probing the structure of high-energy states in protein unfolding

Hydrogen exchange kinetics in native solvent conditions have been used to explore the conformational fluctuations of an immunoglobulin domain (CD2.domain1). The global folding/unfolding kinetics of the protein are unaltered between pH 4.5 and pH 9.5, allowing us to use the pH-dependence of amide hydrogen/deuterium exchange to characterise conformational states with energies up to 7.2kcal/mol higher than the folded ground state. The study was intended to search for discreet unfolding intermediates in this region of the energy spectrum, their presence being revealed by the concerted exchange behaviour of subsets of amide groups that become accessible at a given free energy, i.e. the spectrum would contain discreet groupings. Protection factors for 58 amide groups were measured across the pH range and the hydrogen-exchange energy profile is described. More interestingly, exchange behaviour could be grouped into three categories; the first two unremarkable, the third unexpected. (1) In 33 cases, amide exchange was dominated by rapid fluctuation, i.e. the free energy difference between the ground state and the rapidly accessed open state is sufficiently low that the contribution from crossing the unfolding barrier is negligible. (2) In 18 cases exchange is dominated by the global folding transition barrier across the whole pH range measured. The relationship between hydroxyl ion concentration and observed exchange rate is hyperbolic, with the limiting rate being that for global unfolding; the so-called EX1 limit. For these, the free energy difference between the folded ground state and any rapidly-accessed open state is too great for the proton to be exchanged through such fluctuations, even at the highest pH employed in this study. (3) For the third group, comprising five cases, we observe a behaviour that has not been described. In this group, as in category 2, the rate of exchange reaches a plateau; the EX1 limit. However, as the intrinsic exchange rate (k(int)) is increased, this limit is breached and the rate begins to rise again. This unintuitive behaviour does not result from pH instability, rather it is a consequence of amide groups experiencing two processes; rapid fluctuation of structure and crossing the global barrier for unfolding. The boundary at which the EX1 limit is overcome is determined by the equilibrium distribution of the fluctuating open and closed states (K(O/C)) and the rate constant for unfolding (k(u)). This critical boundary is reached when k(int)K(O/C)=k(u). Given that, in a simple transition state formalism: k(u)=K(#)k' (where K(#) describes the equilibrium distribution between the transition and ground state and k' describes the rate of a barrierless rearrangement), it follows that if the pH is raised to a level where k(int)=k', then the entire free energy spectrum from ground state to transition state could be sampled.

Many residues in cytochromec populate alternative states under equilibrium conditions

A curved temperature dependence of an amide proton NMR chemical shift indicates that it explores discrete alternative conformations at least 1% of the time; that is, it accesses conformations that lie within 5 kcal/mol(-1) of the ground state. The simulations presented show how curvature varies with the nature of the alternative state, and are compared to experimental results. From studies in different denaturant concentrations, it is concluded that at least 25% of residues in reduced horse cytochrome c, covering most of the protein, with the exception of the center of the N- and C-terminal helices, visit alternative states under equilibrium conditions. The conformational ensemble of the protein therefore has high structural entropy. The density of alternative states is particularly high near the heme ligand Met80, which is of interest because both redox change and the first identified stage in unfolding are associated with change in Met80 ligation. By combining theoretical and experimental approaches, it is concluded that the alternative states each comprise approximately five residues, have in general less structure than the native state, and are accessed independently. They are therefore locally unfolded structures. The locations of the alternative states match what is known of the global unfolding pathway of cytochrome c, suggesting that they may determine the pathway.

Hidden intermediates and levinthal paradox in the folding of small proteins

It has long been suggested that existence of partially folded intermediates may be essential for proteins to fold in a biologically meaningful time scale. Although partially folded intermediates have been commonly observed in larger proteins, they are generally not detectable in the kinetic folding of smaller proteins (approximately 100 amino acids or less). Recent native-state hydrogen exchange studies suggest that partially folded intermediates may exist behind the rate-limiting transition state in small proteins and evade detection by conventional kinetic methods.

Hydrogen-exchange stability analysis of Bergerac-Src homology 3 variants allows the characterization of a folding intermediate in equilibrium

Amide hydrogendeuterium exchange rates have been determined for two mutants of alpha-spectrin Src homology 3 domain (WT), containing an elongated stable (SHH) and unstable (SHA) distal loop. SHA, similarly to WT, follows a two-state transition, whereas SHH apparently folds via a three-state mechanism. Native-state amide hydrogen exchange is effective in ascribing energetic readjustments observed in kinetic experiments to species stabilized within the denatured base and distinguishing those from high-energy barrier crossings. Comparison of DeltaG(ex) and m(ex) parameters for amide protons of these mutants demonstrates the existence of an intermediate and allows the identification of protons protected in this state. The consolidation of a form containing a prefolded long beta-hairpin induces the switch to a three-state mechanism in an otherwise two-state folder. It can be inferred that the unbalanced high stability of individual elements of secondary structure in a polypeptide could ultimately complicate its folding mechanism.

Partial NMR assignments and secondary structure mapping of the isolated alpha subunit of Escherichia coli tryptophan synthase, a 29-kD TIM barrel protein

The alpha subunit of tryptophan synthase (alphaTS) from S. typhimurium belongs to the triosephosphate isomerase (TIM) or the (beta/alpha)(8) barrel fold, one of the most common structures in biology. To test the conservation of the global fold in the isolated Escherichia coli homolog, we have obtained a majority of the backbone assignments for the 29-kD alphaTS by using standard heteronuclear multidimensional NMR methods on uniformly (15)N- and (15)N/(13)C-labeled protein and on protein selectively (15)N-labeled at key hydrophobic residues. The secondary structure mapped by chemical shift index, nuclear Overhauser enhancements (NOEs), and hydrogen-deuterium (H-D) exchange, and several abnormal chemical shifts are consistent with the conservation of the global TIM barrel fold of the isolated E. coli alphaTS. Because most of the amide protons that are slow to exchange with solvent correspond to the beta-sheet residues, the beta-barrel is likely to play an important role in stabilizing the previously detected folding intermediates for E. coli alphaTS. A similar combination of uniform and selective labeling can be extended to other TIM barrel proteins to obtain insight into the role of the motif in stabilizing what appear to be common partially folded forms.

Identification of rare partially unfolded states in equilibrium with the native conformation in an all beta-barrel protein

Human acidic fibroblast growth factor 1 (hFGF-1) is an all beta-barrel protein, and the secondary structural elements in the protein include 12 antiparallel beta-strands arranged into a beta-trefoil fold. In the present study, we investigate the stability of hFGF-1 by hydrogen-deuterium exchange as a function of urea concentration. Urea-induced equilibrium unfolding of hFGF-1 monitored by fluorescence and CD spectroscopy suggests that the protein unfolds by a two-state (native to denatured) mechanism. Hydrogen exchange in hFGF-1, under the experimental conditions used, occurs by the EX2 mechanism. In contrast to the equilibrium unfolding events monitored by optical probes, native state hydrogen exchange data show that the beta-trefoil architecture of hFGF-1 does not behave as a single cooperative unit. There are at least two structurally independent units with differing stabilities in hFGF-1. Beta-strands I, II, III, VI, VII, X, XI, and XII fit into the global unfolding isotherm. By contrast, residues in beta-strands IV, V, VIII, and IX exchange by the subfolding isotherm and could be responsible for the occurrence of high-energy partially unfolded state(s) in hFGF-1. There appears to be a broad continuum of stabilities among the four beta-strands (beta-strands IV, V, VIII, and IX) constituting the subglobal folding unit. The slow exchanging residues in hFGF-1 do not represent the folding nucleus of the protein.

The temperature dependence of the hydrogen exchange in the SH3 domain of alpha-spectrin

The amide hydrogen-deuterium exchange (HX) in the Src homology region 3 (SH3) domain of alpha-spectrin has been measured by nuclear magnetic resonance as a function of temperature between 8 and 46 degrees C. The analysis of the temperature dependence of HX from a statistical thermodynamic point of view has allowed us to estimate the enthalpies and entropies of the conformational processes leading to HX. The results indicate that under native conditions the domain undergoes a wide variety of conformational fluctuations, ranging from local motions, mainly located in loops, turns and chain ends and involving only low enthalpy and entropy, to extensive structural disruptions affecting its core and involving enthalpies and entropies that come fairly close to those observed during global unfolding.

The linkage between protein folding and functional cooperativity: two sides of the same coin?

During the course of their biological function, proteins undergo different types of structural rearrangements ranging from local to large-scale conformational changes. These changes are usually triggered by their interactions with small-molecular-weight ligands or other macromolecules. Because binding interactions occur at specific sites and involve only a small number of residues, a chain of cooperative interactions is necessary for the propagation of binding signals to distal locations within the protein structure. This process requires an uneven structural distribution of protein stability and cooperativity as revealed by NMR-detected hydrogen/deuterium exchange experiments under native conditions. The distribution of stabilizing interactions does not only provide the architectural foundation to the three-dimensional structure of a protein, but it also provides the required framework for functional cooperativity. In this review, the statistical thermodynamic linkage between protein stability, functional cooperativity, and ligand binding is discussed.

pH dependence of the hydrogen exchange in the SH3 domain of alpha-spectrin

Using nuclear magnetic resonance we have measured the hydrogen exchange (HX) in the Src homology region 3 (SH3) domain of alpha-spectrin as a function of pH*. At very acidic pH* values the exchange of most residues appears to occur via global unfolding, although several residues show abnormally large Gibbs energies of exchange, suggesting the presence of some residual structure in the unfolded state. At higher pH* HX occurs mainly via local or partial unfoldings. We have been able to characterize the coupling between the electrostatic interactions in this domain and the conformational fluctuations occurring under native conditions by analyzing the dependence upon pH* of the Gibbs energy of exchange. The SH3 domain seems to be composed of a central core, which requires large structural disruptions to become exposed to the solvent, surrounded by smaller subdomains, which fluctuate independently.

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

Native-state amide hydrogen exchange (HX) of proteins in the presence of denaturant has provided valuable details on the structures of equilibrium folding intermediates. Here, we extend HX theory to model thiol group exchange (SX) in single cysteine-containing variants of sperm whale ferric aquomyoglobin. SX is complementary to HX in that it monitors conformational opening events that expose side-chains, rather than the main chain, to solvent. A simple two-process model, consisting of EX2-limited local structural fluctuations and EX1-limited global unfolding, adequately accounts for all HX data. SX is described by the same model except at very low denaturant concentrations and when the bulky labeling reagent 5,5'-dithiobis (2-nitrobenzoic acid) (DTNB) is used. Under these conditions SX can occur by a novel denaturant-dependent process. This anomalous behavior is not observed when the smaller labeling reagent methyl methanethiosulfonate is employed, suggesting that it reflects a denaturant-induced increase in the amplitudes of local structural fluctuations. It also is not seen in heme-free apomyoglobin, which may indicate that local openings are sufficiently large in the absence of denaturant to allow DTNB unhindered access. Differences in SX kinetics obtained using the two labeling reagents provide estimates of the sizes of local opening reactions at different sites in the protein. At all sequence positions examined except for position 73, the same opening event appears to facilitate exchange of both backbone amide and side-chain thiol groups. The C73 thiol group is exposed by a low-energy fluctuation that does not expose its amide group to exchange.