Ultrarapid mixing experiments reveal that Im7 folds via an on-pathway intermediate

Differentiating between the sequence of structural events on alternative pathways of folding of a heterodimeric protein

Distinguishing between competing pathways of folding of a protein, on the basis of how they differ in their progress of structure acquisition, remains an important challenge in protein folding studies. A previous study had shown that the heterodimeric protein, double chain monellin (dcMN) switches between alternative folding pathways upon a change in guanidine hydrochloride (GdnHCl) concentration. In the current study, the folding of dcMN has been characterized by the pulsed hydrogen exchange (HX) labeling methodology used in conjunction with mass spectrometry. Quantification of the extent to which folding intermediates accumulate and then disappear with time of folding at both low and high GdnHCl concentrations, where the folding pathways are known to be different, shows that the folding mechanism is describable by a triangular three-state mechanism. Structural characterization of the productive folding intermediates populated on the alternative pathways has enabled the pathways to be differentiated on the basis of the progress of structure acquisition that occurs on them. The intermediates on the two pathways differ in the extent to which the α-helix and the rest of the β-sheet have acquired structure that is protective against HX. The major difference is, however, that β2 has not acquired any protective structure in the intermediate formed on one pathway, but it has acquired significant protective structure in the intermediate formed on the alternative pathway. Hence, the sequence of structural events is different on the two alternative pathways.

The Molten Globule State of a Globular Protein in a Cell Is More or Less Frequent Case Rather than an Exception

Quite a long time ago, Oleg B. Ptitsyn put forward a hypothesis about the possible functional significance of the molten globule (MG) state for the functioning of proteins. MG is an intermediate between the unfolded and the native state of a protein. Its experimental detection and investigation in a cell are extremely difficult. In the last decades, intensive studies have demonstrated that the MG-like state of some globular proteins arises from either their modifications or interactions with protein partners or other cell components. This review summarizes such reports. In many cases, MG was evidenced to be functionally important. Thus, the MG state is quite common for functional cellular proteins. This supports Ptitsyn’s hypothesis that some globular proteins may switch between two active states, rigid (N) and soft (MG), to work in solution or interact with partners.

Advances in Mixer Design and Detection Methods for Kinetics Studies of Macromolecular Folding and Binding on the Microsecond Time Scale

Many important biological processes such as protein folding and ligand binding are too fast to be fully resolved using conventional stopped-flow techniques. Although advances in mixer design and detection methods have provided access to the microsecond time regime, there is room for improvement in terms of temporal resolution and sensitivity. To address this need, we developed a continuous-flow mixing instrument with a dead time of 12 to 27 µs (depending on solution viscosity) and enhanced sensitivity, sufficient for monitoring tryptophan or tyrosine fluorescence changes at fluorophore concentrations as low as 1 µM. Relying on commercially available laser microfabrication services, we obtained an integrated mixer/flow-cell assembly on a quartz chip, based on a cross-channel configuration with channel dimensions and geometry designed to minimize backpressure. By gradually increasing the width of the observation channel downstream from the mixing region, we are able to monitor a reaction progress time window ranging from ~10 µs out to ~3 ms. By combining a solid-state UV laser with a Galvano-mirror scanning strategy, we achieved highly efficient and uniform fluorescence excitation along the flow channel. Examples of applications, including refolding of acid-denatured cytochrome c triggered by a pH jump and binding of a peptide ligand to a PDZ domain, demonstrate the capability of the technique to resolve fluorescence changes down to the 10 µs time regime on modest amounts of reagents.

Successes and challenges in simulating the folding of large proteins

Computational simulations of protein folding can be used to interpret experimental folding results, to design new folding experiments, and to test the effects of mutations and small molecules on folding. However, whereas major experimental and computational progress has been made in understanding how small proteins fold, research on larger, multidomain proteins, which comprise the majority of proteins, is less advanced. Specifically, large proteins often fold via long-lived partially folded intermediates, whose structures, potentially toxic oligomerization, and interactions with cellular chaperones remain poorly understood. Molecular dynamics based folding simulations that rely on knowledge of the native structure can provide critical, detailed information on folding free energy landscapes, intermediates, and pathways. Further, increases in computational power and methodological advances have made folding simulations of large proteins practical and valuable. Here, using serpins that inhibit proteases as an example, we review native-centric methods for simulating the folding of large proteins. These synergistic approaches range from Gō and related structure-based models that can predict the effects of the native structure on folding to all-atom-based methods that include side-chain chemistry and can predict how disease-associated mutations may impact folding. The application of these computational approaches to serpins and other large proteins highlights the successes and limitations of current computational methods and underscores how computational results can be used to inform experiments. These powerful simulation approaches in combination with experiments can provide unique insights into how large proteins fold and misfold, expanding our ability to predict and manipulate protein folding.

The Effect of Proline cis-trans Isomerization on the Folding of the C-Terminal SH2 Domain from p85

SH2 domains are protein domains that modulate protein–protein interactions through a specific interaction with sequences containing phosphorylated tyrosines. In this work, we analyze the folding pathway of the C-terminal SH2 domain of the p85 regulatory subunit of the protein PI3K, which presents a proline residue in a cis configuration in the loop between the βE and βF strands. By employing single and double jump folding and unfolding experiments, we demonstrate the presence of an on-pathway intermediate that transiently accumulates during (un)folding. By comparing the kinetics of folding of the wild-type protein to that of a site-directed variant of C-SH2 in which the proline was replaced with an alanine, we demonstrate that this intermediate is dictated by the peptidyl prolyl cis-trans isomerization. The results are discussed in the light of previous work on the effect of peptidyl prolyl cis-trans isomerization on folding events.

Functional protein dynamics on uncharted time scales detected by nanoparticle-assisted NMR spin relaxation

Protein function depends critically on intrinsic internal dynamics, which is manifested in distinct ways, such as loop motions that regulate protein recognition and catalysis. Under physiological conditions, dynamic processes occur on a wide range of time scales from subpicoseconds to seconds. Commonly used NMR spin relaxation in solution provides valuable information on very fast and slow motions but is insensitive to the intermediate nanosecond to microsecond range that exceeds the protein tumbling correlation time. Presently, very little is known about the nature and functional role of these motions. It is demonstrated here how transverse spin relaxation becomes exquisitely sensitive to these motions at atomic resolution when studying proteins in the presence of nanoparticles. Application of this novel cross-disciplinary approach reveals large-scale dynamics of loops involved in functionally critical protein-protein interactions and protein-calcium ion recognition that were previously unobservable.

Structural characterization of an on-pathway intermediate and transition state in the folding of the N-terminal SH2 domain from SHP2

Src Homology 2 (SH2) domains are a class of protein domains that present a conserved three-dimensional structure and possess a crucial role in mediating protein-protein interactions. Despite their importance and abundance in the proteome, knowledge about the folding properties of SH2 domain is limited. Here we present an extensive mutational analysis (Φ value analysis) of the folding pathway of the N-SH2 domain of the Src homology region 2 domain-containing phosphatase-2 (SHP2) protein, a 104 residues domain that presents the classical SH2 domain fold (two α-helices flanking a central β-sheet composed of 3-5 antiparallel β-strands), with a fundamental role in mediating the interaction of SHP2 with its substrates and triggering key metabolic pathways in the cell. By analysing folding kinetic data we demonstrated that the folding pathway of N-SH2 presents an obligatory on-pathway intermediate that accumulates during the folding reaction. The production of 24 conservative site-directed variants allowed us to perform a Φ value analysis, by which we could fully characterize the intermediate and the transition state native-like interactions, providing a detailed quantitative analysis of the folding pathway of N-SH2. Results highlight the presence of a hydrophobic nucleus that stabilizes the intermediate, leading to a higher degree of native-like interactions in the transition state. Data are discussed and compared with previous works on SH2 domains.

Complex Folding Landscape of Apomyoglobin at Acidic pH Revealed by Ultrafast Kinetic Analysis of Core Mutants

Under mildly acidic conditions (pH 4-4.5) apomyoglobin (apoMb) adopts a partially structured equilibrium state ( M-state) that structurally resembles a kinetic intermediate encountered at a late stage of folding to the native structure at neutral pH. We have previously reported that the M-state is formed rapidly (<1 ms) via a multistate process and thus offers a unique opportunity for exploring early stages of folding by both experimental and computational techniques. In order to gain structural insight into intermediates and barriers at the residue level, we studied the folding/unfolding kinetics of 12 apoMb mutants at pH 4.2 using fluorescence-detected ultrafast mixing techniques. Global analysis of the submillisecond folding/unfolding kinetics vs urea concentration for each variant, based on a sequential four-state mechanism ( U ⇔ I ⇔ L ⇔ M), allowed us to determine elementary rate constants and their dependence on urea concentration for most transitions. Comparison of the free energy diagrams constructed from the kinetic data of the mutants with that of wild-type apoMb yielded quantitative information on the effects of mutations on the free energy (ΔΔ G) of both intermediates and the first two kinetic barriers encountered during folding. Truncation of conserved aliphatic side chains on helices A, G, and H gives rise to a stepwise increase in ΔΔ G as the protein advances from U toward M, consistent with progressive stabilization of native-like contacts within the primary core of apoMb. Helix-helix contacts in the primary core contribute little to the first folding barrier ( U ⇔ I) and thus are not required for folding initiation but are critical for the stability of the late intermediate, L, and the M-state. Alanine substitution of hydrophobic residues at more peripheral helix-helix contact sites of the native structure, which are still absent or unstable in the M-state, shows both positive (destabilizing) and negative (stabilizing) ΔΔ G, indicating that non-native contacts are formed initially and weakened or lost as a result of subsequent structural rearrangement steps.

Molecular recognition and packing frustration in a helical protein

Biomolecular recognition entails attractive forces for the functional native states and discrimination against potential nonnative interactions that favor alternate stable configurations. The challenge posed by the competition of nonnative stabilization against native-centric forces is conceptualized as frustration. Experiment indicates that frustration is often minimal in evolved biological systems although nonnative possibilities are intuitively abundant. Much of the physical basis of minimal frustration in protein folding thus remains to be elucidated. Here we make progress by studying the colicin immunity protein Im9. To assess the energetic favorability of nonnative versus native interactions, we compute free energies of association of various combinations of the four helices in Im9 (referred to as H1, H2, H3, and H4) by extensive explicit-water molecular dynamics simulations (total simulated time > 300 μs), focusing primarily on the pairs with the largest native contact surfaces, H1-H2 and H1-H4. Frustration is detected in H1-H2 packing in that a nonnative packing orientation is significantly stabilized relative to native, whereas such a prominent nonnative effect is not observed for H1-H4 packing. However, in contrast to the favored nonnative H1-H2 packing in isolation, the native H1-H2 packing orientation is stabilized by H3 and loop residues surrounding H4. Taken together, these results showcase the contextual nature of molecular recognition, and suggest further that nonnative effects in H1-H2 packing may be largely avoided by the experimentally inferred Im9 folding transition state with native packing most developed at the H1-H4 rather than the H1-H2 interface.

Biomolecules need to recognize one another with high specificity: promoting “native” functional intermolecular binding events while avoiding detrimental “nonnative” bound configurations; i.e., “frustration”—the tendency for nonnative interactions—has to be minimized. Folding of globular proteins entails a similar discrimination. To gain physical insight, we computed the binding affinities of helical structures of the protein Im9 in various native or nonnative configurations by atomic simulations, discovering that partial packing of the Im9 core is frustrated. This frustration is overcome when the entire core of the protein is assembled, consistent with experiment indicating no significant kinetic trapping in Im9 folding. Our systematic analysis thus reveals a subtle, contextual aspect of biomolecular recognition and provides a general approach to characterize folding frustration.

Kinetic and thermodynamic studies reveal chemokine homologues CC11 and CC24 with an almost identical tertiary structure have different folding pathways

BACKGROUND

Proteins with low sequence identity but almost identical tertiary structure and function have been valuable to uncover the relationship between sequence, tertiary structure, folding mechanism and functions. Two homologous chemokines, CCL11 and CCL24, with low sequence identity but similar tertiary structure and function, provide an excellent model system for respective studies.

RESULTS

The kinetics and thermodynamics of the two homologous chemokines were systematically characterized. Despite their similar tertiary structures, CCL11 and CCL24 show different thermodynamic stability in guanidine hydrochloride titration, with D_50% = 2.20 M and 4.96 M, respectively. The kinetics curves clearly show two phases in the folding/unfolding processes of both CCL11 and CCL24, which suggests the existence of an intermediate state in their folding/unfolding processes. The folding pathway of both CCL11 and CCL24 could be well described using a folding model with an on-pathway folding intermediate. However, the folding kinetics and stability of the intermediate state of CCL11 and CCL24 are obviously different.

CONCLUSION

Our results suggest homologous proteins with low sequence identity can display almost identical tertiary structure, but very different folding mechanisms, which applies to homologues in the chemokine protein family, extending the general applicability of the above observation.

Unveiling the folding mechanism of the Bromodomains

Bromodomains (BRDs) are small protein domains often present in large multidomain proteins involved in transcriptional regulation in eukaryotic cells. They currently represent valuable targets for the development of inhibitors of aberrant transcriptional processes in a variety of human diseases. Here we report urea-induced equilibrium unfolding experiments monitored by circular dichroism (CD) and fluorescence on two structurally similar BRDs: BRD2(2) and BRD4(1), showing that BRD4(1) is more stable than BRD2(2). Moreover, we report a description of their kinetic folding mechanism, as obtained by careful analysis of stopped-flow and temperature-jump data. The presence of a high energy intermediate for both proteins, suggested by the non-linear dependence of the folding rate on denaturant concentration in the millisec time regime, has been experimentally observed by temperature-jump experiments. Quantitative global analysis of all the rate constants obtained over a wide range of urea concentrations, allowed us to propose a common, three-state, folding mechanism for these two BRDs. Interestingly, the intermediate of BRD4(1) appears to be more stable and structurally native-like than that populated by BRD2(2). Our results underscore the role played by structural topology and sequence in determining and tuning the folding mechanism.

•

A three-state mechanism for the folding of two representative Bromodomains is proposed.

•

Global analyses of BRD2(2) and BRD4(1) folding kinetics highlights the presence of an on-pathway, folding intermediate.

•

The folding intermediate of BRD4(1) is proposed to be more native-like than that of BRD2(2).

The Folding Pathway of the KIX Domain

The KIX domain is an 89-residues globular domain with an important role in mediating protein-protein interactions. The presence of two distinct binding sites in such a small domain makes KIX a suitable candidate to investigate the effect of the potentially divergent demands between folding and function. Here, we report an extensive mutational analysis of the folding pathway of the KIX domain, based on 30 site-directed mutants, which allow us to assess the structures of both the transition and denatured states. Data reveal that, while the transition state presents mostly native-like interactions, the denatured state is somewhat misfolded. We mapped some of the non-native contacts in the denatured state using a second round of mutagenesis, based on double mutant cycles on 15 double mutants. Interestingly, such a misfolding arises from non-native interactions involving the residues critical for the function of the protein. The results described in this work appear to highlight the diverging demands between folding and function that may lead to misfolding, which may be observed in the early stages of folding.

Probing conformational dynamics in biomolecules via chemical exchange saturation transfer: a primer

Although Chemical Exchange Saturation Transfer (CEST) type NMR experiments have been used to study chemical exchange processes in molecules since the early 1960s, there has been renewed interest in the past several years in using this approach to study biomolecular conformational dynamics. The methodology is particularly powerful for the study of sparsely populated, transiently formed conformers that are recalcitrant to investigation using traditional biophysical tools, and it is complementary to relaxation dispersion and magnetization transfer experiments that have traditionally been used to study chemical exchange processes. Here we discuss the concepts behind the CEST experiment, focusing on practical aspects as well, we review some of the pulse sequences that have been developed to characterize protein and RNA conformational dynamics, and we discuss a number of examples where the CEST methodology has provided important insights into the role of dynamics in biomolecular function.

Folding Mechanism of Proteins Im7 and Im9: Insight from All-Atom Simulations in Implicit and Explicit Solvent

Im7 and Im9 are evolutionary related proteins with almost identical native structures. In spite of their structural similarity, experiments show that Im7 folds through a long-lived on-pathway intermediate, while Im9 folds according to two-state kinetics. In this work, we use a recently developed enhanced path sampling method to generate many folding trajectories for these proteins, using realistic atomistic force fields, in both implicit and explicit solvent. Overall, our results are in good agreement with the experimental ϕ values and with the result of ϕ-value-restrained molecular dynamics (MD) simulations. However, our implicit solvent simulations fail to predict a qualitative difference in the folding pathways of Im7 and Im9. In contrast, our simulations in explicit solvent correctly reproduce the fact that only protein Im7 folds through a on-pathway intermediate. By analyzing our atomistic trajectories, we provide a physical picture which explains the observed difference in the folding kinetics of these chains.

Microsecond Rearrangements of Hydrophobic Clusters in an Initially Collapsed Globule Prime Structure Formation during the Folding of a Small Protein

Determining how polypeptide chain collapse initiates structure formation during protein folding is a long standing goal. It has been challenging to characterize experimentally the dynamics of the polypeptide chain, which lead to the formation of a compact kinetic molten globule (MG) in about a millisecond. In this study, the sub-millisecond events that occur early during the folding of monellin from the guanidine hydrochloride-unfolded state have been characterized using multiple fluorescence and fluorescence resonance energy transfer probes. The kinetic MG is shown to form in a noncooperative manner from the unfolded (U) state as a result of at least three different processes happening during the first millisecond of folding. Initial chain compaction completes within the first 37μs, and further compaction occurs only after structure formation commences at a few milliseconds of folding. The transient nonnative and native-like hydrophobic clusters with side chains of certain residues buried form during the initial chain collapse and the nonnative clusters quickly disassemble. Subsequently, partial chain desolvation occurs, leading to the formation of a kinetic MG. The initial chain compaction and subsequent chain rearrangement appear to be barrierless processes. The two structural rearrangements within the collapsed globule appear to prime the protein for the actual folding transition.

Forces Driving Chaperone Action

It is still unclear what molecular forces drive chaperone-mediated protein folding. Here, we obtain a detailed mechanistic understanding of the forces that dictate the four key steps of chaperone-client interaction: initial binding, complex stabilization, folding, and release. Contrary to the common belief that chaperones recognize unfolding intermediates by their hydrophobic nature, we discover that the model chaperone Spy uses long-range electrostatic interactions to rapidly bind to its unfolded client protein Im7. Short-range hydrophobic interactions follow, which serve to stabilize the complex. Hydrophobic collapse of the client protein then drives its folding. By burying hydrophobic residues in its core, the client's affinity to Spy decreases, which causes client release. By allowing the client to fold itself, Spy circumvents the need for client-specific folding instructions. This mechanism might help explain how chaperones can facilitate the folding of various unrelated proteins.

Identification and Structural Characterization of an Intermediate in the Folding of the Measles Virus X Domain

Although most proteins fold by populating intermediates, the transient nature of such states makes it difficult to characterize their structures. In this work we identified and characterized the structure of an intermediate of the X domain of phosphoprotein (P) of measles virus. We obtained this result by a combination of equilibrium and kinetic measurements and NMR chemical shifts used as structural restraints in replica-averaged metadynamics simulations. The structure of the intermediate was then validated by rationally designing four mutational variants predicted to affect the stability of this state. These results provide a detailed view of an intermediate state and illustrate the opportunities offered by a synergistic use of experimental and computational methods to describe non-native states at atomic resolution.

Substrate protein folds while it is bound to the ATP-independent chaperone Spy

Chaperones assist the folding of many proteins in the cell. While the most well studied chaperones use cycles of ATP binding and hydrolysis to assist protein folding, a number of chaperones have been identified that promote protein folding in the absence of high-energy cofactors. Precisely how ATP-independent chaperones accomplish this feat is unclear. Here we have characterized the kinetic mechanism of substrate folding by the small, ATP-independent chaperone, Spy. Spy rapidly associates with its substrate, Immunity protein 7 (Im7), eliminating its potential for aggregation. Remarkably, Spy then allows Im7 to fully fold into its native state while remaining bound to the surface of the chaperone. These results establish a potentially widespread mechanism whereby ATP-independent chaperones can assist in protein refolding. They also provide compelling evidence that substrate proteins can fold while continuously bound to a chaperone.

Mapping the Topography of a Protein Energy Landscape

Protein energy landscapes are highly complex, yet the vast majority of states within them tend to be invisible to experimentalists. Here, using site-directed mutagenesis and exploiting the simplicity of tandem-repeat protein structures, we delineate a network of these states and the routes between them. We show that our target, gankyrin, a 226-residue 7-ankyrin-repeat protein, can access two alternative (un)folding pathways. We resolve intermediates as well as transition states, constituting a comprehensive series of snapshots that map early and late stages of the two pathways and show both to be polarized such that the repeat array progressively unravels from one end of the molecule or the other. Strikingly, we find that the protein folds via one pathway but unfolds via a different one. The origins of this behavior can be rationalized using the numerical results of a simple statistical mechanics model that allows us to visualize the equilibrium behavior as well as single-molecule folding/unfolding trajectories, thereby filling in the gaps that are not accessible to direct experimental observation. Our study highlights the complexity of repeat-protein folding arising from their symmetrical structures; at the same time, however, this structural simplicity enables us to dissect the complexity and thereby map the precise topography of the energy landscape in full breadth and remarkable detail. That we can recapitulate the key features of the folding mechanism by computational analysis of the native structure alone will help toward the ultimate goal of designed amino-acid sequences with made-to-measure folding mechanisms-the Holy Grail of protein folding.

A Native-like Intermediate Serves as a Branching Point between the Folding and Aggregation Pathways of the Mouse Prion Protein

Transient folding intermediates and/or partially unfolded equilibrium states are thought to play a key role in the formation of protein aggregates. However, there is only indirect evidence linking accumulation of folding intermediates to aggregation, and the underlying mechanism remains to be elucidated. Here, we show that a partially unfolded state of the prion protein accumulates both as a stable equilibrium state at acidic pH (A-state) and as a late folding intermediate. With a time resolution of approximately 60 μs, we systematically studied the kinetics of folding and unfolding, starting from various initial conditions including the U-, N-, and A-states. Quantitative modeling showed that the observed kinetic data are completely consistent with a sequential four-state mechanism where the A-state is a late folding intermediate. Combined with previous evidence linking A-state accumulation to aggregation, the results indicate that this native-like state serves as a branching point between the folding and aggregation pathways.

Evidence for a Shared Mechanism in the Formation of Urea-Induced Kinetic and Equilibrium Intermediates of Horse Apomyoglobin from Ultrarapid Mixing Experiments

In this study, the equivalence of the kinetic mechanisms of the formation of urea-induced kinetic folding intermediates and non-native equilibrium states was investigated in apomyoglobin. Despite having similar structural properties, equilibrium and kinetic intermediates accumulate under different conditions and via different mechanisms, and it remains unknown whether their formation involves shared or distinct kinetic mechanisms. To investigate the potential mechanisms of formation, the refolding and unfolding kinetics of horse apomyoglobin were measured by continuous- and stopped-flow fluorescence over a time range from approximately 100 μs to 10 s, along with equilibrium unfolding transitions, as a function of urea concentration at pH 6.0 and 8°C. The formation of a kinetic intermediate was observed over a wider range of urea concentrations (0–2.2 M) than the formation of the native state (0–1.6 M). Additionally, the kinetic intermediate remained populated as the predominant equilibrium state under conditions where the native and unfolded states were unstable (at ~0.7–2 M urea). A continuous shift from the kinetic to the equilibrium intermediate was observed as urea concentrations increased from 0 M to ~2 M, which indicates that these states share a common kinetic folding mechanism. This finding supports the conclusion that these intermediates are equivalent. Our results in turn suggest that the regions of the protein that resist denaturant perturbations form during the earlier stages of folding, which further supports the structural equivalence of transient and equilibrium intermediates. An additional folding intermediate accumulated within ~140 μs of refolding and an unfolding intermediate accumulated in <1 ms of unfolding. Finally, by using quantitative modeling, we showed that a five-state sequential scheme appropriately describes the folding mechanism of horse apomyoglobin.

Native contact density and nonnative hydrophobic effects in the folding of bacterial immunity proteins

The bacterial colicin-immunity proteins Im7 and Im9 fold by different mechanisms. Experimentally, at pH 7.0 and 10°C, Im7 folds in a three-state manner via an intermediate but Im9 folding is two-state-like. Accordingly, Im7 exhibits a chevron rollover, whereas the chevron arm for Im9 folding is linear. Here we address the biophysical basis of their different behaviors by using native-centric models with and without additional transferrable, sequence-dependent energies. The Im7 chevron rollover is not captured by either a pure native-centric model or a model augmented by nonnative hydrophobic interactions with a uniform strength irrespective of residue type. By contrast, a more realistic nonnative interaction scheme that accounts for the difference in hydrophobicity among residues leads simultaneously to a chevron rollover for Im7 and an essentially linear folding chevron arm for Im9. Hydrophobic residues identified by published experiments to be involved in nonnative interactions during Im7 folding are found to participate in the strongest nonnative contacts in this model. Thus our observations support the experimental perspective that the Im7 folding intermediate is largely underpinned by nonnative interactions involving large hydrophobics. Our simulation suggests further that nonnative effects in Im7 are facilitated by a lower local native contact density relative to that of Im9. In a one-dimensional diffusion picture of Im7 folding with a coordinate- and stability-dependent diffusion coefficient, a significant chevron rollover is consistent with a diffusion coefficient that depends strongly on native stability at the conformational position of the folding intermediate.

In order to fold correctly, a globular protein must avoid being trapped in wrong, i.e., nonnative conformations. Thus a biophysical account of how attractive nonnative interactions are bypassed by some amino acid sequences but not others is key to deciphering protein structure and function. We examine two closely related bacterial immunity proteins, Im7 and Im9, that are experimentally known to fold very differently: Whereas Im9 folds directly, Im7 folds through a mispacked conformational intermediate. A simple model we developed accounts for their intriguingly different folding kinetics in terms of a balance between the density of native-promoting contacts and the hydrophobicity of local amino acid sequences. This emergent principle is extensible to other biomolecular recognition processes.

Nonspecific shielding of unfavorable electrostatic intramolecular interactions in the erythropoietin native-state increase conformational stability and limit non-native aggregation

Previous equilibrium and kinetic folding studies of the glycoprotein erythropoietin indicate that sodium chloride increases the conformational stability of this therapeutically important cytokine, ostensibly by stabilizing the native-state [Banks DD, (2011) The Effect of Glycosylation on the Folding Kinetics of Erythropoietin. J Mol Biol 412:536-550]. The focus of the current report is to determine the underlying cause of the salt dependent increase in erythropoietin conformational stability and to understand if it has any impact on aggregation, an instability that remains a challenge to the biotech industry in maintaining the efficacy and shelf-life of protein therapeutics. Isothermal urea denaturation experiments conducted at numerous temperatures in the absence and presence of sodium chloride indicated that salt stabilizes erythropoietin primarily by increasing the difference in enthalpy between the native and unfolded sates. This result, and the finding that the salt induced increases in erythropoietin melting temperatures were independent of the identity of the salt cation and anion, indicates that salt likely increases the conformational stability of erythropoietin at neutral pH by nonspecific shielding of unfavorable electrostatic interaction(s) in the native-state. The addition of salt (even low concentrations of the strong chaotrope salt guanidinium hydrochloride) also exponentially decreased the initial rate of soluble erythropoietin non-native aggregation at 37 °C storage.

The assembly dynamics of the cytolytic pore toxin ClyA

Pore-forming toxins are protein assemblies used by many organisms to disrupt the membranes of target cells. They are expressed as soluble monomers that assemble spontaneously into multimeric pores. However, owing to their complexity, the assembly processes have not been resolved in detail for any pore-forming toxin. To determine the assembly mechanism for the ring-shaped, homododecameric pore of the bacterial cytolytic toxin ClyA, we collected a diverse set of kinetic data using single-molecule spectroscopy and complementary techniques on timescales from milliseconds to hours, and from picomolar to micromolar ClyA concentrations. The entire range of experimental results can be explained quantitatively by a surprisingly simple mechanism. First, addition of the detergent n-dodecyl-β-D-maltopyranoside to the soluble monomers triggers the formation of assembly-competent toxin subunits, accompanied by the transient formation of a molten-globule-like intermediate. Then, all sterically compatible oligomers contribute to assembly, which greatly enhances the efficiency of pore formation compared with simple monomer addition.

Measurement of histidine pKa values and tautomer populations in invisible protein states

The histidine imidazole side chain plays a critical role in protein function and stability. Its importance for catalysis is underscored by the fact that histidines are localized to active sites in ∼ 50% of all enzymes. NMR spectroscopy has become an important tool for studies of histidine side chains through the measurement of site-specific pK(a)s and tautomer populations. To date, such studies have been confined to observable protein ground states; however, a complete understanding of the role of histidine electrostatics in protein function and stability requires that similar investigations be extended to rare, transiently formed conformers that populate the energy landscape, yet are often "invisible" in standard NMR spectra. Here we present NMR experiments and a simple strategy for studies of such conformationally excited states based on measurement of histidine (13)Cγ, (13)Cδ2 chemical shifts and (1)Hε-(13)Cε one-bond scalar couplings. The methodology is first validated and then used to obtain pKa values and tautomer distributions for histidine residues of an invisible on-pathway folding intermediate of the colicin E7 immunity protein. Our results imply that the side chains of H40 and H47 are exposed in the intermediate state and undergo significant conformational rearrangements during folding to the native structure. Further, the pKa values explain the pH-dependent stability differences between native and intermediate states over the pH range 5.5-6.5 and they suggest that imidazole deprotonation is not a barrier to the folding of this protein.

Conformational dynamics is more important than helical propensity for the folding of the all α-helical protein Im7

Im7 folds via an on-pathway intermediate that contains three of the four native α-helices. The missing helix, helix III, is the shortest and its failure to be formed until late in the pathway is related to frustration in the structure. Im7H3M3, a 94-residue variant of the 87-residue Im7 in which helix III is the longest of the four native helices, also folds via an intermediate. To investigate the structural basis for this we calculated the frustration in the structure of Im7H3M3 and used NMR to investigate its dynamics. We found that the native state of Im7H3M3 is highly frustrated and in equilibrium with an intermediate state that lacks helix III, similar to Im7. Model-free analysis identified residues with chemical exchange contributions to their relaxation that aligned with the residues predicted to have highly frustrated interactions, also like Im7. Finally, we determined properties of urea-denatured Im7H3M3 and identified four clusters of interacting residues that corresponded to the α-helices of the native protein. In Im7 the cluster sizes were related to the lengths of the α-helices with cluster III being the smallest but in Im7H3M3 cluster III was also the smallest, despite this region forming the longest helix in the native state. These results suggest that the conformational properties of the urea-denatured states promote formation of a three-helix intermediate in which the residues that form helix III remain non-helical. Thus it appears that features of the native structure are formed early in folding linked to collapse of the unfolded state.

Protein destabilisation in ionic liquids: the role of preferential interactions in denaturation

The preferential binding of anions and cations in aqueous solutions of the ionic liquids (ILs) 1-butyl-3-methylimidazolium ([C4mim](+)) and 1-ethyl-3-methylimidazolium ([C2mim](+)) chloride and dicyanamide (dca(-)) with the small alpha-helical protein Im7 was investigated using a combination of differential scanning calorimetry, NMR spectroscopy and molecular dynamics (MD) simulations. Our results show that direct ion interactions are crucial to understand the effects of ILs on the stability of proteins and that an anion effect is dominant. We show that the binding of weakly hydrated anions to positively charged or polar residues leads to the partial dehydration of the backbone groups, and is critical to control stability, explaining why dca(-) is more denaturing than Cl(-). Direct cation-protein interactions also mediate stability; cation size and hydrophobicity are relevant to account for destabilisation as shown by the effect of [C4mim](+) compared to [C2mim](+). The specificity in the interaction of IL ions with protein residues established by weak favourable interactions is confirmed by NMR chemical shift perturbation, amide hydrogen exchange data and MD simulations. Differences in specificity are due to the balance of interaction established between ion pairs and ion-solvent that determine the type of residues affected. When the interaction of both cation and anion with the protein is strong the net result is similar to a non-specific interaction, leading ultimately to unfolding. Since the nature of the ions is a determinant of the level of interaction with the protein towards denaturation or stabilisation, ILs offer a unique possibility to modulate protein stabilisation or even folding events.

Native contacts determine protein folding mechanisms in atomistic simulations

The recent availability of long equilibrium simulations of protein folding in atomistic detail for more than 10 proteins allows us to identify the key interactions driving folding. We find that the collective fraction of native amino acid contacts, Q, captures remarkably well the transition states for all the proteins with a folding free energy barrier. Going beyond this global picture, we devise two different measures to quantify the importance of individual interresidue contacts in the folding mechanism: (i) the log-ratio of lifetimes of contacts during folding transition paths and in the unfolded state and (ii) a Bayesian measure of how predictive the formation of each contact is for being on a transition path. Both of these measures indicate that native, or near-native, contacts are important for determining mechanism, as might be expected. More remarkably, however, we found that for almost all the proteins, with the designed protein α3D being a notable exception, nonnative contacts play no significant part in determining folding mechanisms.

Microsecond folding experiments and simulations: a match is made

For the past two decades, protein folding experiments have been speeding up from the second or millisecond time scale to the microsecond time scale, and full-atom simulations have been extended from the nanosecond to the microsecond and even millisecond time scale. Where the two meet, it is now possible to compare results directly, allowing force fields to be validated and refined, and allowing experimental data to be interpreted in atomistic detail. In this perspective we compare recent experiments and simulations on the microsecond time scale, pointing out the progress that has been made in determining native structures from physics-based simulations, refining experiments and simulations to provide more quantitative underlying mechanisms, and tackling the problems of multiple reaction coordinates, downhill folding, and complex underlying structure of unfolded or misfolded states.

Probing slowly exchanging protein systems via ¹³Cα-CEST: monitoring folding of the Im7 protein

A¹³C(α) chemical exchange saturation transfer based experiment is presented for the study of protein systems undergoing slow interconversion between an 'observable' ground state and one or more 'invisible' excited states. Here a labeling strategy whereby [2-(13)C]-glucose is the sole carbon source is exploited, producing proteins with ¹³C at the C(α) position, while the majority of residues remain unlabeled at CO or C(β). The new experiment is demonstrated with an application to the folding reaction of the Im7 protein that involves an on-pathway excited state. The obtained excited state (13)C(α) chemical shifts are cross validated by comparison to values extracted from analysis of CPMG relaxation dispersion profiles, establishing the utility of the methodology.

Reassessing the folding of the KIX domain: evidence for a two-state mechanism

The debate about the presence and role of intermediates in the folding of proteins has been a critical issue, especially for fast folders. One of the classical methodologies to identify such metastable species is the "burst-phase analysis," whereby the observed signal amplitude from stopped-flow traces is determined as a function of denaturant concentration. However, a complication may arise when folding is sufficiently fast to jeopardize the reliability of the stopped-flow technique. In this study, we reassessed the folding of the KIX domain from cAMP Response Element-Binding (CREB)-binding protein, which has been proposed to involve the formation of an intermediate that accumulates in the dead time of the stopped flow. By using an in-house-built capillary continuous flow with a 50-μs dead time, we demonstrate that this intermediate is not present; the problem arose because of the instrumental limitation of the standard stopped flow to assess very fast refolding rate constants (e.g., ≥ 500 s⁻¹).

Understanding how small helical proteins fold: conformational dynamics of Im proteins relevant to their folding landscapes

Understanding the mechanism of folding of small proteins requires characterization of their starting unfolded states and any partially unfolded states populated during folding. Here, we review what is known from NMR about these states of Im7, a 4-helix bundle protein that folds via an on-pathway intermediate, and show that there is an alignment of non-native structure in urea-unfolded Im7 with the helices of native Im7 that is a consequence of hydrophobic helix-promoting residues also promoting cluster-formation in the unfolded protein. We suggest that this kind of alignment is present in other proteins and is relevant to how native state topology determines folding rates.

Quantifying millisecond exchange dynamics in proteins by CPMG relaxation dispersion NMR using side-chain 1H probes

A Carr-Purcell-Meiboom-Gill relaxation dispersion experiment is presented for quantifying millisecond time-scale chemical exchange at side-chain (1)H positions in proteins. Such experiments are not possible in a fully protonated molecule because of magnetization evolution from homonuclear scalar couplings that interferes with the extraction of accurate transverse relaxation rates. It is shown, however, that by using a labeling strategy whereby proteins are produced using {(13)C,(1)H}-glucose and D(2)O a significant number of 'isolated' side-chain (1)H spins are generated, eliminating such effects. It thus becomes possible to record (1)H dispersion profiles at the β positions of Asx, Cys, Ser, His, Phe, Tyr, and Trp as well as the γ positions of Glx, in addition to the methyl side-chain moieties. This brings the total of amino acid side-chain positions that can be simultaneously probed using a single (1)H dispersion experiment to 16. The utility of the approach is demonstrated with an application to the four-helix bundle colicin E7 immunity protein, Im7, which folds via a partially structured low populated intermediate that interconverts with the folded, ground state on the millisecond time-scale. The extracted (1)H chemical shift differences at side-chain positions provide valuable restraints in structural studies of invisible, excited states, complementing backbone chemical shifts that are available from existing relaxation dispersion experiments.

Cross-validation of the structure of a transiently formed and low populated FF domain folding intermediate determined by relaxation dispersion NMR and CS-Rosetta

We have recently reported the atomic resolution structure of a low populated and transiently formed on-pathway folding intermediate of the FF domain from human HYPA/FBP11 [Korzhnev, D. M.; Religa, T. L.; Banachewicz, W.; Fersht, A. R.; Kay, L.E. Science 2011, 329, 1312-1316]. The structure was determined on the basis of backbone chemical shift and bond vector orientation restraints of the invisible intermediate state measured using relaxation dispersion nuclear magnetic resonance (NMR) spectroscopy that were subsequently input into the database structure determination program, CS-Rosetta. As a cross-validation of the structure so produced, we present here the solution structure of a mimic of the folding intermediate that is highly populated in solution, obtained from the wild-type domain by mutagenesis that destabilizes the native state. The relaxation dispersion/CS-Rosetta structures of the intermediate are within 2 Å of those of the mimic, with the nonnative interactions in the intermediate also observed in the mimic. This strongly confirms the structure of the FF domain folding intermediate, in particular, and validates the use of relaxation dispersion derived restraints in structural studies of invisible excited states, in general.

Microsecond unfolding kinetics of sheep prion protein reveals an intermediate that correlates with susceptibility to classical scrapie

The microsecond folding and unfolding kinetics of ovine prion proteins (ovPrP) were measured under various solution conditions. A fragment comprising residues 94-233 of the full-length ovPrP was studied for four variants with differing susceptibilities to classical scrapie in sheep. The observed biexponential unfolding kinetics of ovPrP provides evidence for an intermediate species. However, in contrast to previous results for human PrP, there is no evidence for an intermediate under refolding conditions. Global analysis of the kinetic data, based on a sequential three-state mechanism, quantitatively accounts for all folding and unfolding data as a function of denaturant concentration. The simulations predict that an intermediate accumulates under both folding and unfolding conditions, but is observable only in unfolding experiments because the intermediate is optically indistinguishable from the native state. The relative population of intermediates in two ovPrP variants, both transiently and under destabilizing equilibrium conditions, correlates with their propensities for classical scrapie. The variant susceptible to classical scrapie has a larger population of the intermediate state than the resistant variant. Thus, the susceptible variant should be favored to undergo the PrP(C) to PrP(Sc) conversion and oligomerization.

Slow amyloid nucleation via α-helix-rich oligomeric intermediates in short polyglutamine-containing huntingtin fragments

The 17-amino-acid N-terminal segment (htt(NT)) that leads into the polyglutamine (polyQ) segment in the Huntington's disease protein huntingtin (htt) dramatically increases aggregation rates and changes the aggregation mechanism, compared to a simple polyQ peptide of similar length. With polyQ segments near or above the pathological repeat length threshold of about 37, aggregation of htt N-terminal fragments is so rapid that it is difficult to tease out mechanistic details. We describe here the use of very short polyQ repeat lengths in htt N-terminal fragments to slow this disease-associated aggregation. Although all of these peptides, in addition to htt(NT) itself, form α-helix-rich oligomeric intermediates, only peptides with Q(N) of eight or longer mature into amyloid-like aggregates, doing so by a slow increase in β-structure. Concentration-dependent circular dichroism and analytical ultracentrifugation suggest that the htt(NT) sequence, with or without added glutamine residues, exists in solution as an equilibrium between disordered monomer and α-helical tetramer. Higher order, α-helix rich oligomers appear to be built up via these tetramers. However, only htt(NT)Q(N) peptides with N=8 or more undergo conversion into polyQ β-sheet aggregates. These final amyloid-like aggregates not only feature the expected high β-sheet content but also retain an element of solvent-exposed α-helix. The α-helix-rich oligomeric intermediates appear to be both on- and off-pathway, with some oligomers serving as the pool from within which nuclei emerge, while those that fail to undergo amyloid nucleation serve as a reservoir for release of monomers to support fibril elongation. Based on a regular pattern of multimers observed in analytical ultracentrifugation, and a concentration dependence of α-helix formation in CD spectroscopy, it is likely that these oligomers assemble via a four-helix assembly unit. PolyQ expansion in these peptides appears to enhance the rates of both oligomer formation and nucleation from within the oligomer population, by structural mechanisms that remain unclear.

Nuclease colicins and their immunity proteins

It is more than 80 years since Gratia first described 'a remarkable antagonism between two strains of Escherichia coli'. Shown subsequently to be due to the action of proteins (or peptides) produced by one bacterium to kill closely related species with which it might be cohabiting, such bacteriocins have since been shown to be commonplace in the internecine warfare between bacteria. Bacteriocins have been studied primarily from the twin perspectives of how they shape microbial communities and how they penetrate bacteria to kill them. Here, we review the modes of action of a family of bacteriocins that cleave nucleic acid substrates in E. coli, known collectively as nuclease colicins, and the specific immunity (inhibitor) proteins that colicin-producing organisms make in order to avoid committing suicide. In a process akin to targeting in mitochondria, nuclease colicins engage in a variety of cellular associations in order to translocate their cytotoxic domains through the cell envelope to the cytoplasm. As well as informing on the process itself, the study of nuclease colicin import has also illuminated functional aspects of the host proteins they parasitize. We also review recent studies where nuclease colicins and their immunity proteins have been used as model systems for addressing fundamental problems in protein folding and protein-protein interactions, areas of biophysics that are intimately linked to the role of colicins in bacterial competition and to the import process itself.

NMR characterisation of the relationship between frustration and the excited state of Im7

Previous work shows that Im9 folds in a two-state transition while its homologue Im7 folds in a three-state transition via an on-pathway kinetic intermediate state (KIS), with this difference being related to frustration in the structure of Im7. We have used NMR spectroscopy to study conformational dynamics connected to the frustration. A combination of equilibrium peptide N(1)H/N(2)H exchange, model-free analyses of backbone NH relaxation data and relaxation dispersion (RD)-NMR shows that the native state of Im7 is in equilibrium with an intermediate state that is lowly populated [equilibrium intermediate state (EIS)]. Comparison of kinetic and thermodynamic parameters describing the EIS native-state equilibrium obtained by RD-NMR with previously reported parameters describing the KIS native-state equilibrium obtained from stopped-flow fluorescence studies of refolding His-tagged Im7 shows that the KIS and the EIS are the same species. (15)N chemical shifts of the EIS obtained from the RD-NMR analysis show that residues forming helix III in the native state are unstructured in the EIS while other residues experiencing frustration in the native state are in structured regions of the EIS. We show that binding of Im7 and its L53A/I54A variant (which resembles the EIS as shown in previous work) to the cognate partner for Im7, the DNase domain of colicin E7, causes the dynamic processes associated with the frustration to be dampened.

NMR analysis of a kinetically trapped intermediate of a disulfide-deficient mutant of the starch-binding domain of glucoamylase

A thermally unfolded disulfide-deficient mutant of the starch-binding domain of glucoamylase refolds into a kinetically trapped metastable intermediate when subjected to a rapid lowering of temperature. We attempted to characterise this intermediate using multidimensional NMR spectroscopy. The (1)H-(15)N heteronuclear single quantum coherence spectrum after a rapid temperature decrease (the spectrum of the intermediate) showed good chemical shift dispersion but was significantly different from that of the native state, suggesting that the intermediate adopts a nonnative but well-structured conformation. Large chemical shift changes for the backbone amide protons between the native and the intermediate states were observed for residues in the β-sheet consisting of strands 2, 3, 5, 6, and 7 as well as in the C-terminal region. These residues were found to be in close proximity to aromatic residues, suggesting that the chemical shift changes are mainly due to ring current shifts caused by the aromatic residues. The two-dimensional nuclear Overhauser enhancement (NOE) spectroscopy experiments showed that the intermediate contained substantial, native-like NOE connectivities, although there were fewer cross peaks in the spectrum of the intermediate compared with that of the native state. It was also shown that there were native-like interresidue NOEs for residues buried in the protein, whereas many of the NOE cross peaks were lost for the residues involved in a surface-exposed aromatic cluster. These results suggest that, in the intermediate, the aromatic cluster at the surface is structurally less organised, whereas the interior of the protein has relatively rigid, native-like side-chain packing.

Quantifying millisecond time-scale exchange in proteins by CPMG relaxation dispersion NMR spectroscopy of side-chain carbonyl groups

A new pulse sequence is presented for the measurement of relaxation dispersion profiles quantifying millisecond time-scale exchange dynamics of side-chain carbonyl groups in uniformly (13)C labeled proteins. The methodology has been tested using the 87-residue colicin E7 immunity protein, Im7, which is known to fold via a partially structured low populated intermediate that interconverts with the folded, ground state on the millisecond time-scale. Comparison of exchange parameters extracted for this folding 'reaction' using the present methodology with those obtained from more 'traditional' (15)N and backbone carbonyl probes establishes the utility of the approach. The extracted excited state side-chain carbonyl chemical shifts indicate that the Asx/Glx side-chains are predominantly unstructured in the Im7 folding intermediate. However, several crucial salt-bridges that exist in the native structure appear to be already formed in the excited state, either in part or in full. This information, in concert with that obtained from existing backbone and side-chain methyl relaxation dispersion experiments, will ultimately facilitate a detailed description of the structure of the Im7 folding intermediate.

Probing the folding intermediate of Bacillus subtilis RNase P protein by nuclear magnetic resonance

Protein folding intermediates are often imperative for overall folding processes and consequent biological functions. However, the low population and transient nature of the intermediate states often hinder their biochemical and biophysical characterization. Previous studies have demonstrated that Bacillus subtilis ribonuclease P protein (P protein) is conformationally heterogeneous and folds with multiphasic kinetics, indicating the presence of an equilibrium and kinetic intermediate in its folding mechanism. In this study, nuclear magnetic resonance (NMR) spectroscopy was used to study the ensemble corresponding to this intermediate (I). The results indicate that the N-terminal and C-terminal helical regions are mostly unfolded in I. 1H−15N heteronuclear single-quantum coherence NMR spectra collected as a function of pH suggest that the protonation of His 22 may play a major role in the energetics of the equilibria among the unfolded, intermediate, and folded state ensembles of P protein. NMR paramagnetic relaxation enhancement experiments were also used to locate the small anion binding sites in both the intermediate and folded ensembles. The results for the folded protein are consistent with the previously modeled binding regions. These structural insights suggest a possible role for I in the RNase P holoenzyme assembly process.

Dissecting key residues in folding and stability of the bacterial immunity protein 7

The small four-helix immunity protein, Im7, has previously been shown to fold via a compact intermediate containing three of the four native helices. The short, six-residue helix III only docks onto the developing Im7 structure after the rate-limiting second transition state has been traversed. Previous work demonstrated that mutation of the helix III sequence can be used to trap the protein in the on-pathway intermediate ensemble at equilibrium. Here the role played by individual residues in the native helix III sequence in locking Im7 into a stable native structure is further examined. This work commenced with an Im7 sequence trapped in the partially folded state by substitution of the six residues in helix III with a polyglycine sequence. Biophysical analysis of variants in which individual residues from the native helix III sequence, and combinations of these residues, were introduced into this background demonstrated a critical requirement for three residues, Leu 53, Ile 54 and Tyr 55, to lock Im7 into its unique native structure. The results demonstrate a stringent constraint on the evolution of the Im7 helix III sequence rationalizing its high-sequence identity in the fold family. Thus, Leu 53 and Ile 54 provide crucial stabilizing interactions in the hydrophobic core of native Im7, while Tyr 55 is required for both stability and function. In contrast, Tyr 56 is critical for colicin binding and has no role in maintaining a stable native fold.

Osmolyte-induced folding of an intrinsically disordered protein: folding mechanism in the absence of ligand

Understanding the interconversion between thermodynamically distinguishable states present in a protein folding pathway provides not only the kinetics and energetics of protein folding but also insights into the functional roles of these states in biological systems. The protein component of the bacterial RNase P holoenzyme from Bacillus subtilis (P protein) was previously shown to be unfolded in the absence of its cognate RNA or other anionic ligands. P protein was used in this study as a model system to explore general features of intrinsically disordered protein (IDP) folding mechanisms. The use of trimethylamine N-oxide (TMAO), an osmolyte that stabilizes the unliganded folded form of the protein, enabled us to study the folding process of P protein in the absence of ligand. Transient stopped-flow kinetic traces at various final TMAO concentrations exhibited multiphasic kinetics. Equilibrium "cotitration" experiments were performed using both TMAO and urea during the titration to produce a urea-TMAO titration surface of P protein. Both kinetic and equilibrium studies show evidence of a previously undetected intermediate state in the P protein folding process. The intermediate state is significantly populated, and the folding rate constants are relatively slow compared to those of intrinsically folded proteins similar in size and topology. The experiments and analysis described serve as a useful example for mechanistic folding studies of other IDPs.

The plastic energy landscape of protein folding: a triangular folding mechanism with an equilibrium intermediate for a small protein domain

Protein domains usually fold without or with only transiently populated intermediates, possibly to avoid misfolding, which could result in amyloidogenic disease. Whether observed intermediates are productive and obligatory species on the folding reaction pathway or dispensable by-products is a matter of debate. Here, we solved the crystal structure of a small protein domain, SAP97 PDZ2 I342W C378A, and determined its folding pathway. The presence of a folding intermediate was demonstrated both by single and double-mixing kinetic experiments using urea-induced (un)folding as well as ligand-induced folding. This protein domain was found to fold via a triangular scheme, where the folding intermediate could be either on- or off-pathway, depending on the experimental conditions. Furthermore, we found that the intermediate was present at equilibrium, which is rarely seen in folding reactions of small protein domains. The folding mechanism observed here illustrates the roughness and plasticity of the protein folding energy landscape, where several routes may be employed to reach the native state. The results also reconcile the folding mechanisms of topological variants within the PDZ domain family.