Detection and characterization of an early folding intermediate of T4 lysozyme using pulsed hydrogen exchange and two-dimensional NMR

DMSO-Quenched H/D-Exchange 2D NMR Spectroscopy and Its Applications in Protein Science

Hydrogen/deuterium (H/D) exchange combined with two-dimensional (2D) NMR spectroscopy has been widely used for studying the structure, stability, and dynamics of proteins. When we apply the H/D-exchange method to investigate non-native states of proteins such as equilibrium and kinetic folding intermediates, H/D-exchange quenching techniques are indispensable, because the exchange reaction is usually too fast to follow by 2D NMR. In this article, we will describe the dimethylsulfoxide (DMSO)-quenched H/D-exchange method and its applications in protein science. In this method, the H/D-exchange buffer is replaced by an aprotic DMSO solution, which quenches the exchange reaction. We have improved the DMSO-quenched method by using spin desalting columns, which are used for medium exchange from the H/D-exchange buffer to the DMSO solution. This improvement has allowed us to monitor the H/D exchange of proteins at a high concentration of salts or denaturants. We describe methodological details of the improved DMSO-quenched method and present a case study using the improved method on the H/D-exchange behavior of unfolded human ubiquitin in 6 M guanidinium chloride.

OUP accepted manuscript

Summary

Motivation. Predicting the native state of a protein has long been considered a gateway problem for understanding protein folding. Recent advances in structural modeling driven by deep learning have achieved unprecedented success at predicting a protein’s crystal structure, but it is not clear if these models are learning the physics of how proteins dynamically fold into their equilibrium structure or are just accurate knowledge-based predictors of the final state.

Results. In this work, we compare the pathways generated by state-of-the-art protein structure prediction methods to experimental data about protein folding pathways. The methods considered were AlphaFold 2, RoseTTAFold, trRosetta, RaptorX, DMPfold, EVfold, SAINT2 and Rosetta. We find evidence that their simulated dynamics capture some information about the folding pathway, but their predictive ability is worse than a trivial classifier using sequence-agnostic features like chain length. The folding trajectories produced are also uncorrelated with experimental observables such as intermediate structures and the folding rate constant. These results suggest that recent advances in structure prediction do not yet provide an enhanced understanding of protein folding.

Availability. The data underlying this article are available in GitHub at https://github.com/oxpig/structure-vs-folding/

Supplementary information

Supplementary data are available at Bioinformatics online.

The Effects of p-Azidophenylalanine Incorporation on Protein Structure and Stability

Functionalization is often needed to harness the power of proteins for beneficial use but can cause losses to stability and/or activity. State of the art methods to limit these deleterious effects accomplish this by substituting an amino acid in the wild-type molecule into an unnatural amino acid, such as p-azidophenylalanine (pAz), but selecting the residue for substitution a priori remains an elusive goal of protein engineering. The results of this work indicate that all-atom molecular dynamics simulation can be used to determine whether substituting pAz for a natural amino acid will be detrimental to experimentally determined protein stability. These results offer significant hope that local deviations from wild-type structure caused by pAz incorporation observed in simulations can be a predictive metric used to reduce the number of costly experiments that must be done to find active proteins upon substitution with pAz and subsequent functionalization.

A Double-Resonance CEST Experiment To Study Multistate Protein Conformational Exchange: An Application to Protein Folding

Despite the importance of protein dynamics to function, studying exchange between multiple conformational states remains a challenge because sparsely populated states are invisible to conventional techniques. CEST NMR experiments can detect minor states with lifetimes between 5 and 200 ms populated to a level of just ∼1%. However, CEST often cannot provide the exchange mechanism for processes involving three or more states, leaving the role of the detected minor states unknown. Here a double-resonance CEST experiment to determine the kinetics of multistate exchange is presented. The approach that involves irradiating resonances from two minor states simultaneously is used to study the exchange of T4 lysozyme (T4L) between the dominant native state and two minor states, the unfolded state and a second minor state (B), each populated to only ∼4%. Regular CEST does not provide the folding mechanism, but double-resonance CEST clearly shows that T4L can fold directly without going through B.

Probing the folded state and mechanical unfolding pathways of T4 lysozyme using all-atom and coarse-grained molecular simulation

The Bacteriophage T4 Lysozyme (T4L) is a prototype modular protein comprised of an N-terminal and a C-domain domain, which was extensively studied to understand the folding/unfolding mechanism of modular proteins. To offer detailed structural and dynamic insights to the folded-state stability and the mechanical unfolding behaviors of T4L, we have performed extensive equilibrium and steered molecular dynamics simulations of both the wild-type (WT) and a circular permutation (CP) variant of T4L using all-atom and coarse-grained force fields. Our all-atom and coarse-grained simulations of the folded state have consistently found greater stability of the C-domain than the N-domain in isolation, which is in agreement with past thermostatic studies of T4L. While the all-atom simulation cannot fully explain the mechanical unfolding behaviors of the WT and the CP variant observed in an optical tweezers study, the coarse-grained simulations based on the Go model or a modified elastic network model (mENM) are in qualitative agreement with the experimental finding of greater unfolding cooperativity in the WT than the CP variant. Interestingly, the two coarse-grained models predict different structural mechanisms for the observed change in cooperativity between the WT and the CP variant--while the Go model predicts minor modification of the unfolding pathways by circular permutation (i.e., preserving the general order that the N-domain unfolds before the C-domain), the mENM predicts a dramatic change in unfolding pathways (e.g., different order of N/C-domain unfolding in the WT and the CP variant). Based on our simulations, we have analyzed the limitations of and the key differences between these models and offered testable predictions for future experiments to resolve the structural mechanism for cooperative folding/unfolding of T4L.

Environmental polarity induces conformational transitions in a helical peptide sequence from bacteriophage T4 lysozyme and its tandem duplicate: a molecular dynamics simulation study

Our recent molecular dynamics (MD) simulation of an insertion/duplication mutant 'L20' of bacteriophage T4 lysozyme demonstrated a solvent induced α→β transition in a loosely held duplicate helical region, while α-helical conformation in the parent region was relatively stabilized by its tertiary interactions with the neighboring residues. The solution NMR of the parent helical sequence, sans its protein context, showed no inherent tendency to adopt a particular secondary structure in pure water but showed α-helical propensity in TFE/water and SDS micelles. In this study we investigate the conformational preference of the 'parent' and 'duplicate' sequences, sans the protein context, in pure water and an apolar TFE/water solution. Apolar TFE/water solution is a model for non-polar protein context. We performed MD simulations of the two peptides, in explicit water and 80% (v/v) TFE/water, using GROMOS 53a6 force field, at 300 K and 1 bar (under NPT conditions). We show that in TFE/water mixture, salt bridges are stabilized by apolar TFE molecules and main chain-main chain hydrogen bonds promote the α-helical conformation, particularly in the duplicate peptide. Solvent exposure, in pure water, resulted in an α→β transition to form a triple stranded β-sheet structure in the 'duplicate' sequence, with a rare psi-loop topology, while a mixture of turn/bend conformations were adopted by the 'parent' sequence. Thus the differences in conformational preference of the parent and duplicate sequence sans protein context, in pure water and TFE/water, implicate the importance of the environment polarity in dictating the peptide conformation. Mechanism of folding of the observed psi-loop in the duplicate sequence gives insights into folding of this rare β-sheet topology.

Modulation of a protein free-energy landscape by circular permutation

Circular permutations usually retain the native structure and function of a protein while inevitably perturbing its folding dynamics. By using simulations with a structure-based model and a rigorous methodology to determine free-energy surfaces from trajectories, we evaluate the effect of a circular permutation on the free-energy landscape of the protein T4 lysozyme. We observe changes which, although subtle, largely affect the cooperativity between the two subdomains. Such a change in cooperativity has been previously experimentally observed and recently also characterized using single molecule optical tweezers and the Crooks relation. The free-energy landscapes show that both the wild type and circular permutant have an on-pathway intermediate, previously experimentally characterized, in which one of the subdomains is completely formed. The landscapes, however, differ in the position of the rate-limiting step for folding, which occurs before the intermediate in the wild type and after in the circular permutant. This shift of transition state explains the observed change in the cooperativity. The underlying free-energy landscape thus provides a microscopic description of the folding dynamics and the connection between circular permutation and the loss of cooperativity experimentally observed.

Probing protein interactions with hydrogen/deuterium exchange and mass spectrometry-a review

Assessing the functional outcome of protein interactions in structural terms is a goal of structural biology, however most techniques have a limited capacity for making structure-function determinations with both high resolution and high throughput. Mass spectrometry can be applied as a reader of protein chemistries in order to fill this void, and enable methodologies whereby protein structure-function determinations may be made on a proteome-wide level. Protein hydrogen/deuterium exchange (H/DX) offers a chemical labeling strategy suitable for tracking changes in "dynamic topography" and thus represents a powerful means of monitoring protein structure-function relationships. This review presents the exchange method in the context of interaction analysis. Applications involving interface detection, quantitation of binding, and conformational responses to ligation are discussed, and commentary on recent analytical developments is provided.

Folding of an all-helical Greek-key protein monitored by quenched-flow hydrogen-deuterium exchange and NMR spectroscopy

To advance our understanding of the protein folding process, we use stopped-flow far-ultraviolet (far-UV) circular dichroism and quenched-flow hydrogen-deuterium exchange coupled with nuclear magnetic resonance (NMR) spectroscopy to monitor the formation of hydrogen-bonded secondary structure in the C-terminal domain of the Fas-associated death domain (Fadd-DD). The death domain superfamily fold consists of six α-helices arranged in a Greek-key topology, which is shared by the all-β-sheet immunoglobulin and mixed α/β-plait superfamilies. Fadd-DD is selected as our model death domain protein system because the structure of this protein has been solved by NMR spectroscopy, and both thermodynamic and kinetic analysis indicate it to be a stable, monomeric protein with a rapidly formed hydrophobic core. Stopped-flow far-UV circular dichroism spectroscopy revealed that the folding process was monophasic and the rate is 23.4 s(-1). Twenty-two amide hydrogens in the backbone of the helices and two in the backbone of the loops were monitored, and the folding of all six helices was determined to be monophasic with rate constants between 19 and 22 s(-1). These results indicate that the formation of secondary structure is largely cooperative and concomitant with the hydrophobic collapse. This study also provides unprecedented insight into the formation of secondary structure within the highly populated Greek-key fold more generally. Additional insights are gained by calculating the exchange rates of 23 residues from equilibrium hydrogen-deuterium exchange experiments. The majority of protected amide protons are found on helices 2, 4, and 5, which make up core structural elements of the Greek-key topology.

High-pressure EPR reveals conformational equilibria and volumetric properties of spin-labeled proteins

Identifying equilibrium conformational exchange and characterizing conformational substates is essential for elucidating mechanisms of function in proteins. Site-directed spin labeling has previously been employed to detect conformational changes triggered by some event, but verifying conformational exchange at equilibrium is more challenging. Conformational exchange (microsecond-millisecond) is slow on the EPR time scale, and this proves to be an advantage in directly revealing the presence of multiple substates as distinguishable components in the EPR spectrum, allowing the direct determination of equilibrium constants and free energy differences. However, rotameric exchange of the spin label side chain can also give rise to multiple components in the EPR spectrum. Using spin-labeled mutants of T4 lysozyme, it is shown that high-pressure EPR can be used to: (i) demonstrate equilibrium between spectrally resolved states, (ii) aid in distinguishing conformational from rotameric exchange as the origin of the resolved states, and (iii) determine the relative partial molar volume (ΔV°) and isothermal compressibility (Δβ(T)) of conformational substates in two-component equilibria from the pressure dependence of the equilibrium constant. These volumetric properties provide insight into the structure of the substates. Finally, the pressure dependence of internal side-chain motion is interpreted in terms of volume fluctuations on the nanosecond time scale, the magnitude of which may reflect local backbone flexibility.

Rapid collapse into a molten globule is followed by simple two-state kinetics in the folding of lysozyme from bacteriophage λ

Stopped-flow fluorescence and circular dichroism spectroscopy have been used in combination with quenched-flow hydrogen exchange labeling, monitored by two-dimensional NMR and electrospray ionization mass spectrometry, to investigate the folding kinetics of lysozyme from bacteriophage λ (λ lysozyme) at pH 5.6, 20 °C. The first step in the folding of λ lysozyme occurs very rapidly (τ < 1 ms) after refolding is initiated and involves both hydrophobic collapse and formation of a high content of secondary structure but only weak protection from (1)H/(2)H exchange and no fixed tertiary structure organization. This early folding step is reflected in the dead-time events observed in the far-UV CD and ANS fluorescence experiments. Following accumulation of this kinetic molten globule species, the secondary structural elements are stabilized and the majority (ca. 88%) of refolding molecules acquire native-like properties in a highly cooperative two-state process, with τ = 0.15 ± 0.03 s. This is accompanied by the acquisition of substantial native-like protection from hydrogen exchange. A double-mixing experiment and the absence of a denaturant effect reveal that slow (τ = 5 ± 1 s) folding of the remaining (ca. 12%) molecules is rate limited by the cis/trans isomerization of prolines that are trans in the folded enzyme. In addition, native state hydrogen exchange and classical denaturant unfolding experiments have been used to characterize the thermodynamic properties of the enzyme. In good agreement with previous crystallographic evidence, our results show that λ lysozyme is a highly dynamic protein, with relatively low conformational stability (ΔG°(N-U) = 25 ± 2 kJ·mol(-1)).

Lessons from the lysozyme of phage T4

An overview is presented of some of the major insights that have come from studies of the structure, stability, and folding of T4 phage lysozyme. A major purpose of this review is to provide the reader with a complete tabulation of all of the variants that have been characterized, including melting temperatures, crystallographic data, Protein Data Bank access codes, and references to the original literature. The greatest increase in melting temperature (T(m)) for any point mutant is 5.1 degrees C for the mutant Ser 117 --> Val. This is achieved in part not only by hydrophobic stabilization but also by eliminating an unusually short hydrogen bond of 2.48 A that apparently has an unfavorable van der Waals contact. Increases in T(m) of more than 3-4 degrees C for point mutants are rare, whereas several different types of destabilizing substitutions decrease T(m) by 20 degrees C or thereabouts. The energetic cost of cavity creation and its relation to the hydrophobic effect, derived from early studies of "large-to-small" mutants in the core of T4 lysozyme, has recently been strongly supported by related studies of the intrinsic membrane protein bacteriorhodopsin. The L99A cavity in the C-terminal domain of the protein, which readily binds benzene and many other ligands, has been the subject of extensive study. Crystallographic evidence, together with recent NMR analysis, suggest that these ligands are admitted by a conformational change involving Helix F and its neighbors. A total of 43 nonisomorphous crystal forms of different monomeric lysozyme mutants were obtained plus three more for synthetically-engineered dimers. Among the 43 space groups, P2(1)2(1)2(1) and P2(1) were observed most frequently, consistent with the prediction of Wukovitz and Yeates.

Directionality in protein fold prediction

BACKGROUND

Ever since the ground-breaking work of Anfinsen et al. in which a denatured protein was found to refold to its native state, it has been frequently stated by the protein fold prediction community that all the information required for protein folding lies in the amino acid sequence. Recent in vitro experiments and in silico computational studies, however, have shown that cotranslation may affect the folding pathway of some proteins, especially those of ancient folds. In this paper aspects of cotranslational folding have been incorporated into a protein structure prediction algorithm by adapting the Rosetta program to fold proteins as the nascent chain elongates. This makes it possible to conduct a pairwise comparison of folding accuracy, by comparing folds created sequentially from each end of the protein.

RESULTS

A single main result emerged: in 94% of proteins analyzed, following the sense of translation, from N-terminus to C-terminus, produced better predictions than following the reverse sense of translation, from the C-terminus to N-terminus. Two secondary results emerged. First, this superiority of N-terminus to C-terminus folding was more marked for proteins showing stronger evidence of cotranslation and second, an algorithm following the sense of translation produced predictions comparable to, and occasionally better than, Rosetta.

CONCLUSIONS

There is a directionality effect in protein fold prediction. At present, prediction methods appear to be too noisy to take advantage of this effect; as techniques refine, it may be possible to draw benefit from a sequential approach to protein fold prediction.

A general and efficient method for the site-specific dual-labeling of proteins for single molecule fluorescence resonance energy transfer

A general strategy for the site-specific dual-labeling of proteins for single-molecule fluorescence resonance energy transfer is presented. A genetically encoded unnatural ketone amino acid was labeled with a hydroxylamine-containing fluorophore with high yield (>95%) and specificity. This methodology was used to construct dual-labeled T4 lysozyme variants, allowing the study of T4 lysozyme folding at single-molecule resolution. The presented strategy is anticipated to expand the scope of single-molecule protein structure and function studies.

Subdomain competition, cooperativity, and topological frustration in the folding of CheY

The folding of multidomain proteins often proceeds in a hierarchical fashion with individual domains folding independent of one another. A large single-domain protein, however, can consist of multiple modules whose folding may be autonomous or interdependent in ways that are unclear. We used coarse-grained simulations to explore the folding landscape of the two-subdomain bacterial response regulator CheY. Thermodynamic and kinetic characterization shows the landscape to be highly analogous to the four-state landscape reported for another two-subdomain protein, T4 lysozyme. An on-pathway intermediate structured in the more stable nucleating subdomain was observed, as were transient states frustrated in off-pathway contacts prematurely structured in the weaker subdomain. Local unfolding, or backtracking, was observed in the frustrated state before the native conformation could be reached. Nonproductive frustration was attributable to competition for van der Waals contacts between the two subdomains. In an accompanying article, stopped-flow kinetic measurements support an off-pathway burst-phase intermediate, seemingly consistent with our prediction of early frustration in the folding landscape of CheY. Comparison of the folding mechanisms for CheY, T4 lysozyme, and interleukin-1 beta leads us to postulate that subdomain competition is a general feature of large single-domain proteins with multiple folding modules.

Atomic force microscopy reveals parallel mechanical unfolding pathways of T4 lysozyme: evidence for a kinetic partitioning mechanism

Kinetic partitioning is predicted to be a general mechanism for proteins to fold into their well defined native three-dimensional structure from unfolded states following multiple folding pathways. However, experimental evidence supporting this mechanism is still limited. By using single-molecule atomic force microscopy, here we report experimental evidence supporting the kinetic partitioning mechanism for mechanical unfolding of T4 lysozyme, a small protein composed of two subdomains. We observed that on stretching from its N and C termini, T4 lysozyme unfolds by multiple distinct unfolding pathways: the majority of T4 lysozymes unfold in an all-or-none fashion by overcoming a dominant unfolding kinetic barrier; and a small fraction of T4 lysozymes unfold in three-state fashion involving unfolding intermediate states. The three-state unfolding pathways do not follow well defined routes, instead they display variability and diversity in individual unfolding pathways. The unfolding intermediate states are local energy minima along the mechanical unfolding pathways and are likely to result from the residual structures present in the two subdomains after crossing the main unfolding barrier. These results provide direct evidence for the kinetic partitioning of the mechanical unfolding pathways of T4 lysozyme, and the complex unfolding behaviors reflect the stochastic nature of kinetic barrier rupture in mechanical unfolding processes. Our results demonstrate that single-molecule atomic force microscopy is an ideal tool to investigate the folding/unfolding dynamics of complex multimodule proteins that are otherwise difficult to study using traditional methods.

Dissecting the Non-specific and Specific Components of the Initial Folding Reaction of Barstar by Multi-site FRET Measurements

Initial polypeptide chain collapse plays a major role in the development of subsequent structure during protein folding, but it has been difficult to elucidate the coupling between its cooperativity and specificity. To better understand this important aspect of protein folding, nine different intramolecular distances in the protein have been measured by fluorescence resonance energy transfer (FRET) in the product(s) of the initial, sub-millisecond collapse reaction during the folding of barstar, under different folding conditions. All nine distances contract in these initial folding products, when the denaturant concentration is reduced. Two of these distances were also measured in peptides corresponding to sequence segments 38-55 and 51-69 of the protein. Surprisingly, both distances do not contract in the peptides which remain fully unfolded when the denaturant concentration is reduced. This suggests that the contraction of at least some segments of the polypeptide chain may be facilitated only by contraction of other segments. In the case of the initial product of folding of the protein, the dependence on denaturant concentration of the relative change in each distance suggests that there are two components to the initial folding reaction. One is a nonspecific component, which appears to be driven by the change in denaturant concentration that is used to initiate refolding. This component corresponds to the collapse of completely unfolded protein (U) to unfolded protein in refolding conditions (U(C)). The extent of nonspecific collapse can be predicted by the response of completely unfolded protein to a change in denaturant concentration. All distances undergo such solvent-induced contraction, but each distance contracts to a different extent. There is also a specific component to initial sub-millisecond folding, in which some distances (but not all) contract more than that predicted by solvent-induced contraction. The observation that only some of the distances undergo contraction over and above solvent-induced contraction, suggest that this specific component is associated with the formation of a specific intermediate (I(E)). FRET efficiency and distance change differently for the different donor-acceptor pairs, with a change in denaturant concentration, indicating that the formation or dissolution of structure in U(C) and I(E) does not happen in a synchronized manner across different regions of the protein molecule. Also, all nine FRET efficiencies and intramolecular distances in the product(s) of sub-ms folding, change continuously with a change in denaturant concentration. Hence, it appears that the transitions from U to U(C) and to I(E) are gradual transformations, and not all-or-none structural transitions. Nevertheless, the product of these gradual transitions, I(E), possesses specific structure.

Exploring subdomain cooperativity in T4 lysozyme I: structural and energetic studies of a circular permutant and protein fragment

Small proteins are generally observed to fold in an apparent two-state manner. Recently, however, more sensitive techniques have demonstrated that even seemingly single-domain proteins are actually made up of smaller subdomains. T4 lysozyme is one such protein. We explored the relative autonomy of its two individual subdomains and their contribution to the overall stability of T4 lysozyme by examining a circular permutation (CP13*) that relocates the N-terminal A-helix, creating subdomains that are contiguous in sequence. By determining the high-resolution structure of CP13* and characterizing its energy landscape using native state hydrogen exchange (NSHX), we show that connectivity between the subdomains is an important determinant of the energetic cooperativity but not structural integrity of the protein. The circular permutation results in a protein more easily able to populate a partially unfolded form in which the C-terminal subdomain is folded and the N-terminal subdomain is unfolded. We also created a fragment model of this intermediate and demonstrate using X-ray crystallography that its structure is identical to the corresponding residues in the full-length protein with the exception of a small network of hydrophobic interactions. In sum, we conclude that the C-terminal subdomain dominates the energetics of T4 lysozyme folding, and the A-helix serves an important role in coupling the two subdomains.

Exploring subdomain cooperativity in T4 lysozyme II: uncovering the C-terminal subdomain as a hidden intermediate in the kinetic folding pathway

Intermediates along a protein's folding pathway can play an important role in its biology. Previous kinetics studies have revealed an early folding intermediate for T4 lysozyme, a small, well-characterized protein composed of an N-terminal and a C-terminal subdomain. Pulse-labeling hydrogen exchange studies suggest that residues from both subdomains contribute to the structure of this intermediate. On the other hand, equilibrium native state hydrogen experiments have revealed a high-energy, partially unfolded form of the protein that has an unstructured N-terminal subdomain and a structured C-terminal subdomain. To resolve this discrepancy between kinetics and equilibrium data, we performed detailed kinetics analyses of the folding and unfolding pathways of T4 lysozyme, as well as several point mutants and large-scale variants. The data support the argument for the presence of two distinct intermediates, one present on each side of the rate-limiting transition state barrier. The effects of circular permutation and site-specific mutations in the wild-type and circular permutant background, as well as a fragment containing just the C-terminal subdomain, support a model for the unfolding intermediate with an unfolded N-terminal and a folded C-terminal subdomain. Our results suggest that the partially unfolded form identified by native state hydrogen exchange resides on the folded side of the rate-limiting transition state and is, therefore, under most conditions, a "hidden" intermediate.

Cores and pH-dependent dynamics of ferredoxin-NADP+ reductase revealed by hydrogen/deuterium exchange

NMR-detected hydrogen/deuterium (H/D) exchange of amide protons is a powerful way for investigating the residue-based conformational stability and dynamics of proteins in solution. Maize ferredoxin-NADP(+) reductase (FNR) is a relatively large protein with 314 amino acid residues, consisting of flavin adenine dinucleotide (FAD) and nicotinamide adenine dinucleotide phosphate (NADP(+))-binding domains. To address the structural stability and dynamics of FNR, H/D exchange of amide protons was performed using heteronuclear NMR at pD(r) values 8.0 and 6.0, physiologically relevant conditions mimicking inside of chloroplasts. At both pD(r) values, the exchange rate varied widely depending on the residues. The profiles of protected residues revealed that the highly protected regions matched well with the hydrophobic cores suggested from the crystal structure, and that the NADP(+)-binding domain can be divided into two subdomains. The global stability of FNR obtained by H/D exchange with NMR was higher than that by chemical denaturation, indicating that H/D exchange is especially useful for analyzing the residue-based conformational stability of large proteins, for which global unfolding is mostly irreversible. Interestingly, more dynamic conformation of the C-terminal subdomain of the NADP(+)-binding domain at pD(r) 8.0, the daytime pH in chloroplasts, than at pD(r) 6.0 is likely to be involved in the increased binding of NADP(+) for elevating the activity of FNR. In light of photosynthesis, the present study provides the first structure-based relationship of dynamics with function for the FNR-type family in solution.

The folding pathway of T4 lysozyme: the high-resolution structure and folding of a hidden intermediate

Folding intermediates have been detected and characterized for many proteins. However, their structures at atomic resolution have only been determined for two small single domain proteins: Rd-apocytochrome b(562) and engrailed homeo domain. T4 lysozyme has two easily distinguishable but energetically coupled domains: the N and C-terminal domains. An early native-state hydrogen exchange experiment identified an intermediate with the C-terminal domain folded and the N-terminal domain unfolded. We have used a native-state hydrogen exchange-directed protein engineering approach to populate this intermediate and demonstrated that it is on the folding pathway and exists after the rate-limiting step. Here, we determined its high-resolution structure and the backbone dynamics by multi-dimensional NMR methods. We also characterized the folding behavior of the intermediate using stopped-flow fluorescence, protein engineering, and native-state hydrogen exchange. Unlike the folding intermediates of the two single-domain proteins, which have many non-native side-chain interactions, the structure of the hidden folding intermediate of T4 lysozyme is largely native-like. It folds like many small single domain proteins. These results have implications for understanding the folding mechanism and evolution of multi-domain proteins.

The folding pathway of T4 lysozyme: an on-pathway hidden folding intermediate

T4 lysozyme has two easily distinguishable but energetically coupled domains: the N and C-terminal domains. In earlier studies, an amide hydrogen/deuterium exchange pulse-labeling experiment detected a stable submillisecond intermediate that accumulates before the rate-limiting transition state. It involves the formation of structures in both the N and C-terminal regions. However, a native-state hydrogen exchange experiment subsequently detected an equilibrium intermediate that only involves the formation of the C-terminal domain. Here, using stopped-flow circular dichroism and fluorescence, amide hydrogen exchange-folding competition, and protein engineering methods, we re-examined the folding pathway of T4-lysozyme. We found no evidence for the existence of a stable folding intermediate before the rate-limiting transition state at neutral pH. In addition, using native-state hydrogen exchange-directed protein engineering, we created a mimic of the equilibrium intermediate. We found that the intermediate mimic folds with the same rate as the wild-type protein, suggesting that the equilibrium intermediate is an on-pathway intermediate that exists after the rate-limiting transition state.

Sequential reorganization of beta-sheet topology by insertion of a single strand

Insertions, duplications, and deletions of sequence segments are thought to be major evolutionary mechanisms that increase the structural and functional diversity of proteins. Alternative splicing, for example, is an intracellular editing mechanism that is thought to generate isoforms for 30%-50% of all human genes. Whereas the inserted sequences usually display only minor structural rearrangements at the insertion site, recent observations indicate that they may also cause more dramatic structural displacements of adjacent structures. In the present study we test how artificially inserted sequences change the structure of the beta-sheet region in T4 lysozyme. Copies of two different beta-strands were inserted into two different loops of the beta-sheet, and the structures were determined. Not surprisingly, one insert "loops out" at its insertion site and forms a new small beta-hairpin structure. Unexpectedly, however, the second insertion leads to displacement of adjacent strands and a sequential reorganization of the beta-sheet topology. Even though the insertions were performed at two different sites, looping out occurred at the C-terminal end of the same beta-strand. Reasons as to why a non-native sequence would be recruited to replace that which occurs in the native protein are discussed. Our results illustrate how sequence insertions can facilitate protein evolution through both local and nonlocal changes in structure.

Multiple methionine substitutions are tolerated in T4 lysozyme and have coupled effects on folding and stability

In order to further explore the tolerance of proteins to amino acid substitutions within the interior, a series of core residues was replaced by methionine within the C-terminal domain of T4 lysozyme. By replacing leucine, isoleucine, valine and phenylalanine residues a total of 10 methionines could be introduced, which corresponds to a third of the residues that are buried in this domain. As more methionines are incorporated the protein gradually loses stability. This is attributed in part to a reduction in hydrophobic stabilization, in part to the increased entropic cost of localizing the long, flexible methionine sidechains, and in part to steric clashes. The changes in structure of the mutants relative to the wildtype protein are modest but tend to increase in an additive fashion as more methionines are included. In the most extreme case, namely the 10-methionine mutant, much of the C-terminal domain remains quite similar to wildtype (root-mean-square backbone shifts of 0.56 A), while the F and G helices undergo rotations of approximately 20 degrees and center-of-mass shifts of approximately 1.4 A. For up to six methionine substitutions the changes in stability are additive. Beyond this point, however, the multiple mutants are somewhat more stable than suggested from the sum of their constituents, especially for those including the replacement Val111-->Met. This is interpreted in terms of the larger structural changes associated with this substitution. The substituted sidechains in the mutant structures have somewhat higher crystallographic thermal factors than their counterparts in WT*. Nevertheless, the interiors of the mutant proteins retain a well-defined structure with little suggestion of molten-globule characteristics. Lysozymes in which selenomethionine has been incorporated rather than methionine tend to have increased stability. At the same time they also fold faster. This provides further evidence that, at the rate-limiting step in folding, the structure of the C-terminal domain of T4 lysozyme is similar to that of the fully folded protein.

A first-order-like state transition for recombinant protein folding

Normally, proteins will aggregate and precipitate by direct folding processes. In this study, we report that quasi-static processes can restore both the structure and bio-function of two kinds of fish recombinant growth hormones (Plecoglossus altivelis and Epinephelus awoara). The conformational changes and the particle-size-distribution (PSD) of each refolding intermediate can be monitored by circular dichroism spectroscopy (CD) and dynamic light scattering (DLS), respectively. Conformation analysis of the CD spectra of the refolding intermediates indicated that the secondary structures were restored in the initial refolding intermediate. However, the tertiary interactions of the proteins were restored during the last two refolding stages, as elucidated by thermal stability tests. This is consistent with a sequential model. DLS analysis suggested that the average hydrodynamic radii of the refolding intermediates shrank to their native-like sizes after the first refolding stage. This is consistent with a collapse model. After comparison with the data on the direct folding process, it is concluded that the denaturant-containing protein folding reaction is a first-order-like state transition process.

Protein folding by a quasi-static-like process: a first-order state transition

In this paper we report that quasi-static-like processes, in which stable intermediates were introduced carefully and deliberately, may be used to reversibly unfold and refold purified native porcine growth hormone. Through circular dichroism (CD) and dynamic light scattering (DLS), we were able to study the secondary structure conformational changes, tertiary structure thermal stabilities, and the particle size distributions of both the intermediates and the final folded product. The CD data showed that the secondary structure was restored in the initial folding stage, whereas the tertiary structure within the protein was restored one step before the last folding stage, as elucidated by thermal stability experiments. DLS analysis suggested that the average hydrodynamic radii of the folding intermediates shrunk to nativelike size immediately after the first folding stage. Our data suggested that the denaturant-containing protein folding reaction is a first-order-like state transition process. This quasi-static-like process may be useful in the prevention of aggregate formation in protein purification and thus can be used in protein engineering to improve the overall yield from harvesting proteins.

Crystal structures of a T4-lysozyme duplication-extension mutant demonstrate that the highly conserved beta-sheet region has low intrinsic folding propensity

Residues 24 to 35 of T4 lysozyme correspond to the second and third strands of a region of beta-sheet that is highly conserved in all known lysozyme and chitinase structures. To evaluate the intrinsic propensity of these amino acid residues to form a defined structure they were added at the C terminus of the native protein, together with a dipeptide linker. Two crystal structures of this active, mutant protein were obtained, to 1.9A and 2.3A resolution, respectively. Even though the crystal conditions are similar, the appended sequence adopts very different secondary structures. In one case it is weakly structured and appears to extend through the active-site cleft, perhaps in part adding an extra strand to the original beta-sheet. In the other crystal form the extension is largely alpha-helical. The formation of these alternative structures shows that the sequence does not have a strong intrinsic propensity to form a unique fold (either beta-sheet or otherwise). The results also suggest that structural conservation during evolution does not necessarily depend on sequence conservation or the conservation of folding propensity.

Comparison of the folding processes of T. thermophilus and E. coli ribonucleases H

In order to examine how the stabilization of thermophilic proteins affects their folding, we have characterized the folding process of Thermus thermophilus ribonuclease H using circular dichroism, fluorescence, and pulse-labeling hydrogen exchange. Like its homolog from Escherichia coli, this thermophilic protein populates a partially folded kinetic intermediate within the first few milliseconds of folding. The structure of this intermediate is similar to that of E.coli RNase H and corresponds remarkably well to a partially folded form that is populated at low levels in the native state of the protein. Proline isomerization appears to partly limit the folding of the thermophilic but not the mesophilic protein. Lastly, unlike other thermophilic proteins, which unfold much more slowly than their mesophilic counterparts, T.thermophilus RNase H folds and unfolds with overall rates similar to those of E.coli RNase H.

Highly divergent dihydrofolate reductases conserve complex folding mechanisms

To test the hypothesis that protein folding mechanisms are better conserved than amino acid sequences, the mechanisms for dihydrofolate reductases (DHFR) from human (hs), Escherichia coli (ec) and Lactobacillus casei (lc) were elucidated and compared using intrinsic Trp fluorescence and fluorescence-detected 8-anilino-1-naphthalenesulfonate (ANS) binding. The development of the native state was monitored using either methotrexate (absorbance at 380 nm) or NADPH (extrinsic fluorescence) binding. All three homologs displayed complex unfolding and refolding kinetic mechanisms that involved partially folded states and multiple energy barriers. Although the pairwise sequence identities are less than 30 %, folding to the native state occurs via parallel folding channels and involves two types of on-pathway kinetic intermediates for all three homologs. The first ensemble of kinetic intermediates, detected within a few milliseconds, has significant secondary structure and exposed hydrophobic cores. The second ensemble is obligatory and has native-like side-chain packing in a hydrophobic core; however, these intermediates are unable to bind active-site ligands. The formation of the ensemble of native states occurs via three channels for hsDHFR, and four channels for lcDHFR and ecDHFR. The binding of active-site ligands (methotrexate and NADPH) accompanies the rate-limiting formation of the native ensemble. The conservation of the fast, intermediate and slow-folding events for this complex alpha/beta motif provides convincing evidence for the hypothesis that evolutionarily related proteins achieve the same fold via similar pathways.

Structural characterization of protein kinase A as a function of nucleotide binding. Hydrogen-deuterium exchange studies using matrix-assisted laser desorption ionization-time of flight mass spectrometry detection

Transient state kinetic studies indicate that substrate phosphorylation in protein kinase A is partially rate-limited by conformational changes, some of which may be associated with nucleotide binding (Shaffer, J., and Adams, J. A. (1999) Biochemistry 38, 12072-12079). To assess whether specific structural changes are associated with the binding of nucleotides, hydrogen-deuterium exchange experiments were performed on the enzyme in the absence and presence of ADP. Four regions of the protein are protected from exchange in the presence of ADP. Two regions encompass the catalytic and glycine-rich loops and are integral parts of the active site. Conversely, protection of probes in the C terminus is consistent with nucleotide-induced domain closure. One protected probe encompasses a portion of helix C, a secondary structural element that does not make any direct contacts with the nucleotide but has been reported to undergo segmental motion upon the activation of some protein kinases. The combined data suggest that binding of the nucleotide has distal structural effects that may include stabilizing the closed state of the enzyme and altering the position of a critical helix outside the active site. The latter represents the first evidence that the nucleotide alone can induce changes in helix C in solution.

Protein folding intermediates and pathways studied by hydrogen exchange

In order to solve the immensely difficult protein-folding problem, it will be necessary to characterize the barriers that slow folding and the intermediate structures that promote it. Although protein-folding intermediates are not accessible to the usual structural studies, hydrogen exchange (HX) methods have been able to detect and characterize intermediates in both kinetic and equilibrium modes--as transient kinetic folding intermediates on a subsecond time scale, as labile equilibrium molten globule intermediates under destabilizing conditions, and as infinitesimally populated intermediates in the high free-energy folding landscape under native conditions. Available results consistently indicate that protein-folding landscapes are dominated by a small number of discrete, metastable, native-like partially unfolded forms (PUFs). The PUFs appear to be produced, one from another, by the unfolding and refolding of the protein's intrinsically cooperative secondary structural elements, which can spontaneously create stepwise unfolding and refolding pathways. Kinetic experiments identify three kinds of barrier processes: (a) an initial intrinsic search-nucleation-collapse process that prepares the chain for intermediate formation by pinning it into a condensed coarsely native-like topology; (b) smaller search-dependent barriers that put the secondary structural units into place; and (c) optional error-dependent misfold-reorganization barriers that can cause slow folding, intermediate accumulation, and folding heterogeneity. These conclusions provide a coherent explanation for the grossly disparate folding behavior of different globular proteins in terms of distinct folding pathways.

Structural events during the refolding of an all beta-sheet protein

The refolding kinetics of the 140-residue, all beta-sheet, human fibroblast growth factor (hFGF-1) is studied using a variety of biophysical techniques such as stopped-flow fluorescence, stopped-flow circular dichroism, and quenched-flow hydrogen exchange in conjunction with multidimensional NMR spectroscopy. Urea-induced unfolding of hFGF-1 under equilibrium conditions reveals that the protein folds via a two-state (native <--> unfolded) mechanism without the accumulation of stable intermediates. However, measurement of the unfolding and refolding rates in various concentrations of urea shows that the refolding of hFGF-1 proceeds through accumulation of kinetic intermediates. Results of the quenched-flow hydrogen exchange experiments reveal that the hydrogen bonds linking the N- and C-terminal ends are the first to form during the refolding of hFGF-1. The basic beta-trefoil framework is provided by the simultaneous formation of beta-strands I, IV, IX, and X. The other beta-strands comprising the beta-barrel structure of hFGF-1 are formed relatively slowly with time constants ranging from 4 to 13 s.

A kinetic folding intermediate probed by native state hydrogen exchange

Stopped-flow fluorescence studies on the N-terminal domain of rat CD2 (CD2.d1) have demonstrated that folding from the fully denatured state (U) proceeds via the transient accumulation of an apparent intermediate (I) in a so-called burst phase that precedes the rate-limiting transition leading to the native state (N). A previous pH-dependent equilibrium hydrogen exchange (HX) study identified a subset of amides in CD2.d1 which, under EX2 conditions, exchange from N with free energies greater than or equal to the free energy difference between the N and I states calculated from the stopped-flow data. Under EX1 conditions the rates of HX for these amides tend towards an asymptote that matches the global unfolding rate calculated from the stopped-flow data, suggesting that exchange for these amides requires traversing the N-to-I transition state barrier. Exchange for these amides presumably occurs from exchange-competent forms comprising the kinetic burst phase therefore. To explore this idea further, native state HX (NHX) data have been collected for CD2.d1 under EX2 conditions using denaturant concentrations which span either side of the denaturant concentration where, according to the stopped-flow data, the apparent U and I states are iso-energetic. The data fit to a two-component, sub-global (sg)/global (g) NHX mechanism, yielding Delta G and m value parameters (where the m value is a measure of hydrocarbon solvation). Regression analysis demonstrates that the (m(sg), Delta G(sg)) and (m(g), Delta G(g)) values calculated for this subset of amides correspond with those describing the kinetic burst phase transition. This result confirms the ability of the NHX technique to explore the structural and energetic properties of kinetic folding intermediates.

Anatomy of protein structures: visualizing how a one-dimensional protein chain folds into a three-dimensional shape

Here, we depict the anatomy of protein structures in terms of the protein folding process. Via an iterative, top-down dissecting procedure, tertiary structures are spliced down to reveal their anatomy: first, to produce domains (defined by visual three-dimensional inspection criteria); then, hydrophobic folding units (HFU); and, at the end of a multilevel process, a set of building blocks. The resulting anatomy tree organization not only clearly depicts the organization of a one-dimensional polypeptide chain in three-dimensional space but also straightforwardly describes the most likely folding pathway(s). Comparison of the tree with the formation of the hydrophobic folding units through combinatorial assembly of the building blocks illustrates how the chain folds in a sequential or a complex folding pathway. Further, the tree points to the kinetics of the folding, whether the chain is a fast or a slow folder, and the probability of misfolding. Our ability to successfully dissect the protein into an anatomy tree illustrates that protein folding is a hierarchical process and further validates a building blocks protein folding model.

A statistical appraisal of native state hydrogen exchange data: evidence for a burst phase continuum?

For a number of proteins, folding occurs via the rapid accumulation of secondary and tertiary structural features in a so-called burst phase, preceding the relatively slow, highly activated transition leading to the native state. A fundamental question is: do these burst phase reactions comprise two phase-separated thermodynamic states or a continuum of states? Ribonuclease HI (RNase H) from Escherichia coli and phage T4 lysozyme (T4L) both exhibit such a phenomenon. Native-state hydrogen exchange (NHX) data have been collected for these proteins, providing residue-specific free energies and m-values (a measure of hydrocarbon solvation) for the manifold of partially unfolded, exchange-competent forms that are accessible from the native state (DeltaG(sg) and m(sg), where the sg subscript denotes sub-global). There is good evidence that these parameters pertain to exchange-competent species comprising the burst phase observed in the global folding kinetics. We combine the results from the global folding kinetics of these proteins with a statistical analysis of their NHX parameters to determine if the distribution of experimental (m(sg), DeltaG(sg)) values derive from a mechanism where the burst phase is two-state. For RNase H, this analysis demonstrates that the burst phase of this protein is not two-state; the results imply a distribution of states, m and DeltaG exhibiting a linear functional relationship consistent with the global folding parameters. For T4L, it is difficult to distinguish the observed distribution of m(sg), DeltaG(sg) values from that expected for a mechanism where the burst phase is two-state. The results for RNase H* lend support for the idea that the burst phase reaction of this protein comprises a continuum of states. This has important implications for how we model the process of structural acquisition in protein folding reactions.

A physical basis for protein secondary structure

A physical theory of protein secondary structure is proposed and tested by performing exceedingly simple Monte Carlo simulations. In essence, secondary structure propensities are predominantly a consequence of two competing local effects, one favoring hydrogen bond formation in helices and turns, the other opposing the attendant reduction in sidechain conformational entropy on helix and turn formation. These sequence specific biases are densely dispersed throughout the unfolded polypeptide chain, where they serve to preorganize the folding process and largely, but imperfectly, anticipate the native secondary structure.

Structure-based analysis of protein dynamics: comparison of theoretical results for hen lysozyme with X-ray diffraction and NMR relaxation data

An analytical approach based on Gaussian network model (GNM) is proposed for predicting the rotational dynamics of proteins. The method, previously shown to successfully reproduce X-ray crystallographic temperature factors for a series of proteins is extended here to predict bond torsional mobilities and reorientation of main chain amide groups probed by 15N-H nuclear magnetic resonance (NMR) relaxation. The dynamics of hen egg-white lysozyme (HEWL) in the folded state is investigated using the proposed approach. Excellent agreement is observed between theoretical results and experimental (X-ray diffraction and NMR relaxation) data. The analysis reveals the important role of coupled rotations, or cross-correlations between dihedral angle librations, in defining the relaxation mechanism on a local scale. The crystal and solution structures exhibit some differences in their local motions, but their global motions are identical. Hinge residues mediating the cooperative movements of the alpha- and beta-domains are identified, which comprise residues in helix C, Glu35 and Ser36 on the loop succeeding helix B, Ile55 and Leu56 at the turn between strands II and III. The central part of the beta-domain long loop and the turn between strands I and II display an enhanced mobility. Finally, kinetically hot residues and key interactions are identified, which point at helix B and beta-strand III as the structural elements underlying the stability of the tertiary structure.

The cooperativity of burst phase reactions explored

The denaturant-dependence of the major, observable relaxation rates for folding (kobs) of ribonuclease HI from Escherichia coli (RNase H) and phage T4 lysozyme (T4L) reveal that, for both proteins, folding begins with the rapid and transient accumulation of intermediate species in a "burst phase" which precedes the rate-limiting formation of the native state; this is evidenced by a "rollover" in the folding limb of the rate profiles (kobs versus denaturant, or chevron plot). These rate profiles are most simply described by a three-state mechanism (unfolded-to-intermediate-to-native), which implies that the burst phase represents a transition between two distinct thermodynamic states. It is shown here that the equilibrium properties of these burst phase reactions can be equally well modeled by a mechanism involving a continuum of states where the free energy of each state is linearly related to its m-value (the parameter describing the linear relationship between free energy and denaturant). A numerical model is also developed to describe the time evolution of such a system, which exhibits nearly perfect exponential behavior. Both models emphasize how a continuum of states operating under a linear free energy relationship may behave like a two state system. Such a scheme finds experimental justification from an interpretation of recent native state hydrogen exchange data. The analytical model described for a continuum can account for the observed kinetic profiles of several other model proteins. The results, however, appear context specific, suggesting that burst phase reactions are not entirely random and non-specific. The results reported in this study have important implications for the concept of cooperativity in protein folding reactions.

Identification of side-chain clusters in protein structures by a graph spectral method

This paper presents a novel method to detect side-chain clusters in protein three-dimensional structures using a graph spectral approach. Protein side-chain interactions are represented by a labeled graph in which the nodes of the graph represent the Cbeta atoms and the edges represent the distance between the Cbeta atoms. The distance information and the non-bonded connectivity of the residues are represented in the form of a matrix called the Laplacian matrix. The constructed matrix is diagonalized and clustering information is obtained from the vector components associated with the second lowest eigenvalue and cluster centers are obtained from the vector components associated with the top eigenvalues. The method uses global information for clustering and a single numeric computation is required to detect clusters of interest. The approach has been adopted here to detect a variety of side-chain clusters and identify the residue which makes the largest number of interactions among the residues forming the cluster (cluster centers). Detecting such clusters and cluster centers are important from a protein structure and folding point of view. The crucial residues which are important in the folding pathway as determined by PhiF values (which is a measure of the effect of a mutation on the stability of the transition state of folding) as obtained from protein engineering methods, can be identified from the vector components corresponding to the top eigenvalues. Expanded clusters are detected near the active and binding site of the protein, supporting the nucleation condensation hypothesis for folding. The method is also shown to detect domains in protein structures and conserved side-chain clusters in topologically similar proteins.

The hydrogen exchange core and protein folding

A database of hydrogen-deuterium exchange results has been compiled for proteins for which there are published rates of out-exchange in the native state, protection against exchange during folding, and out-exchange in partially folded forms. The question of whether the slow exchange core is the folding core (Woodward C, 1993, Trends Biochem Sci 18:359-360) is reexamined in a detailed comparison of the specific amide protons (NHs) and the elements of secondary structure on which they are located. For each pulsed exchange or competition experiment, probe NHs are shown explicitly; the large number and broad distribution of probe NHs support the validity of comparing out-exchange with pulsed-exchange/competition experiments. There is a strong tendency for the same elements of secondary structure to carry NHs most protected in the native state, NHs first protected during folding, and NHs most protected in partially folded species. There is not a one-to-one correspondence of individual NHs. Proteins for which there are published data for native state out-exchange and theta values are also reviewed. The elements of secondary structure containing the slowest exchanging NHs in native proteins tend to contain side chains with high theta values or be connected to a turn/loop with high theta values. A definition for a protein core is proposed, and the implications for protein folding are discussed. Apparently, during folding and in the native state, nonlocal interactions between core sequences are favored more than other possible nonlocal interactions. Other studies of partially folded bovine pancreatic trypsin inhibitor (Barbar E, Barany G, Woodward C, 1995, Biochemistry 34:11423-11434; Barber E, Hare M, Daragan V, Barany G, Woodward C, 1998, Biochemistry 37:7822-7833), suggest that developing cores have site-specific energy barriers between microstates, one disordered, and the other(s) more ordered.

Distinguishing between sequential and nonsequentially folded proteins: implications for folding and misfolding

We describe here an algorithm for distinguishing sequential from nonsequentially folding proteins. Several experiments have recently suggested that most of the proteins that are synthesized in the eukaryotic cell may fold sequentially. This proposed folding mechanism in vivo is particularly advantageous to the organism. In the absence of chaperones, the probability that a sequentially folding protein will misfold is reduced significantly. The problem we address here is devising a procedure that would differentiate between the two types of folding patterns. Footprints of sequential folding may be found in structures where consecutive fragments of the chain interact with each other. In such cases, the folding complexity may be viewed as being lower. On the other hand, higher folding complexity suggests that at least a portion of the polypeptide backbone folds back upon itself to form three-dimensional (3D) interactions with noncontiguous portion(s) of the chain. Hence, we look at the mechanism of folding of the molecule via analysis of its complexity, that is, through the 3D interactions formed by contiguous segments on the polypeptide chain. To computationally splice the structure into consecutively interacting fragments, we either cut it into compact hydrophobic folding units or into a set of hypothetical, transient, highly populated, contiguous fragments ("building blocks" of the structure). In sequential folding, successive building blocks interact with each other from the amino to the carboxy terminus of the polypeptide chain. Consequently, the results of the parsing differentiate between sequentially vs. nonsequentially folded chains. The automated assessment of the folding complexity provides insight into both the likelihood of misfolding and the kinetic folding rate of the given protein. In terms of the funnel free energy landscape theory, a protein that truly follows the mechanism of sequential folding, in principle, encounters smoother free energy barriers. A simple sequentially folded protein should, therefore, be less error prone and fold faster than a protein with a complex folding pattern.

Hydrogen exchange kinetics of proteins in denaturants: a generalized two-process model

The recent progress in measurements on the amide hydrogen exchange (HX) in proteins under varying denaturing conditions, both at equilibrium and in transient relaxation, necessitates the development of a unifying theory which quantitatively relates the HX rates to the conformational energetics of the proteins. We present here a comprehensive kinetic model for the site-specific HX of proteins under varying solvent denaturing conditions based on the two-state protein folding model. The generalized two-process model considers both conformational fluctuations and residual protections, respectively, within the folded and unfolded states of a protein, as well as a global kinetic folding-unfolding transition between the two states. The global transition can be either rapid or slow, depending on the solvent condition for the protein. This novel model is applicable to the traditional equilibrium HX measurements in both EX2 and EX1 regimes, and also the recently introduced transient pulse-labeling HX experiments. A set of simple analytical equations is provided for quantitative interpretation of experimental data. The model emphasizes the use of full time-course of bi-exponential HX kinetics, rather than fitting time-course data to single rate constants, to obtain quantitative information about fluctuating conformers within the folded and unfolded states of proteins. This HX kinetic model naturally unfolds into a simple two-state and two-stage kinetic interpretation for protein folding. It suggests that the various observed intermediates of a protein can be interpreted as dominant isomers of either the folded or the unfolded state under different solvent conditions. This simple, minimalist's view of protein folding is consistent with various recent experimental observations on folding kinetics by HX.

Structure, thermostability, and conformational flexibility of hen egg-white lysozyme dissolved in glycerol

Hen egg-white lysozyme dissolved in glycerol containing 1% water was studied by using CD and amide proton exchange monitored by two-dimensional 1H NMR. The far- and near-UV CD spectra of the protein showed that the secondary and tertiary structures of lysozyme in glycerol were similar to those in water. Thermal melting of lysozyme in glycerol followed by CD spectral changes indicated unfolding of the tertiary structure with a Tm of 76.0 +/- 0.2 degreesC and no appreciable loss of the secondary structure up to 85 degreesC. This is in contrast to the coincident denaturation of both tertiary and secondary structures with Tm values of 74.8 +/- 0.4 degreesC and 74.3 +/- 0.7 degreesC, respectively, under analogous conditions in water. Quenched amide proton exchange experiments revealed a greater structural protection of amide protons in glycerol than in water for a majority of the slowly exchanging protons. The results point to a highly ordered, native-like structure of lysozyme in glycerol, with the stability exceeding that in water.

Amide protection in an early folding intermediate of cytochrome c

BACKGROUND

For many proteins, compact states appear long before the rate-limiting step in the formation of the native structure. A key issue is whether the initial collapse of the chain is driven by random or more specific hydrophobic interactions.

RESULTS

Hydrogen-exchange labeling coupled with NMR was used to monitor the formation of stable hydrogen-bonded and solvent-excluded structure in horse cytochrome c (cyt c). Protection was measured using a hydrogen exchange/folding competition protocol at variable pH and short competition time (2 ms). Protection factors of threefold to eightfold were observed in all three alpha helices of cyt c, whereas other regions showed no significant protection. This suggests that the compact states that are present contain segments of marginally stable hydrogen-bonded structure. When the intermediate(s) are destabilized, only amide protons from Cys14, Ala15 and His18 show significant protection, indicating a region of persistent residual structure near the covalently bound heme group in the unfolded protein. Fluorescence-detected stopped-flow studies showed that the maximum protection factor in the early intermediate is consistent with its unfolding equilibrium constant.

CONCLUSIONS

Together with previous fluorescence and CD results, the observed pattern of amide protection is consistent with the early formation of an alpha-helical core domain in an ensemble of compact states, indicating that efficient folding is facilitated by stepwise acquisition of native structural elements. These specific early interactions are established on the sub-millisecond time scale, prior to the rate-limiting step for folding.

Generation of ligand binding sites in T4 lysozyme by deficiency-creating substitutions

Several variants of T4 lysozyme have been identified that sequester small organic ligands in cavities or clefts. To evaluate potential binding sites for non-polar molecules, we screened a number of hydrophobic large-to-small mutants for stabilization in the presence of benzene. In addition to Leu99-->Ala, binding was indicated for at least five other mutants. Variants Met102-->Ala and Leu133-->Gly, and a crevice mutant, Phe104-->Ala, were further characterized using X-ray crystallography and thermal denaturation. As predicted from the shape of the cavity in the benzene complex, mutant Leu133-->Gly also bound p-xylene. We attempted to enlarge the cavity of the Met102-->Ala mutant into a deep crevice through an additional substitution, but the double mutant failed to bind ligands because an adjacent helix rearranged into a non-helical structure, apparently due to the loss of packing interactions. In general, the protein structure contracted slightly to reduce the volume of the void created by truncating substitutions and expanded upon binding the non-polar ligand, with shifts similar to those resulting from the mutations.A polar molecule binding site was also created by truncating Arg95 to alanine. This creates a highly complementary buried polar environment that can be utilized as a specific "receptor" for a guanidinium ion. Our results suggest that creating a deficiency through truncating mutations of buried residues generates "binding potential" for ligands with characteristics similar to the deleted side-chain. Analysis of complex and apo crystal structures of binding and non-binding mutants suggests that ligand size and shape as well as protein flexibility and complementarity are all determinants of binding. Binding at non-polar sites is governed by hydrophobicity and steric interactions and is relatively permissive. Binding at a polar site is more restrictive and requires extensive complementarity between the ligand and the site.

Hydrogen exchange studies of protein structure

Hydrogen exchange techniques, with their residue-level specificity, exquisite sensitivity, and adaptability to many solution conditions, are becoming essential to the study of protein stability, folding and dynamics. Recent studies have elucidated the structures of intermediates formed transiently during protein folding and rare partially folded ensembles present at equilibrium. Analysis of hydrogen exchange mechanisms has revealed protein stability and kinetics at the level of individual residues.

Subdomain interactions as a determinant in the folding and stability of T4 lysozyme

The folding of large, multidomain proteins involves the hierarchical assembly of individual domains. It remains unclear whether the stability and folding of small, single-domain proteins occurs through a comparable assembly of small, autonomous folding units. We have investigated the relationship between two subdomains of the protein T4 lysozyme. Thermodynamically, T4 lysozyme behaves as a cooperative unit and the unfolding transition fits a two-state model. The structure of the protein, however, resembles a dumbbell with two potential subdomains: an N-terminal subdomain (residues 13-75), and a C-terminal subdomain (residues 76-164 and 1-12). To investigate the effect of uncoupling these two subdomains within the context of the native protein, we created two circular permutations, both at the subdomain interface (residues 13 and 75). Both variants adopt an active wild-type T4 lysozyme fold. The protein starting with residue 13 is 3 kcal/mol less stable than wild type, whereas the protein beginning at residue 75 is 9 kcal/mol less stable, suggesting that the placement of the termini has a major effect on protein stability while minimally affecting the fold. When isolated as protein fragments, the C-terminal subdomain folds into a marginally stable helical structure, whereas the N-terminal subdomain is predominantly unfolded. ANS fluorescence studies indicate that, at low pH, the C-terminal subdomain adopts a loosely packed acid state. An acid state intermediate is also seen for all of the full-length variants. We propose that this acid state is comprised of an unfolded N-terminal subdomain and a loosely folded C-terminal subdomain.