A protein folding pathway with multiple folding intermediates at atomic resolution

Enzymatic and synthetic regulation of polypeptide folding

Proper folding is essential for the biological functions of all proteins. The folding process is intrinsically error-prone, and the misfolding of a polypeptide chain can cause the formation of toxic aggregates related to pathological outcomes such as neurodegenerative disease and diabetes. Chaperones and some enzymes are involved in the cellular proteostasis systems that assist polypeptide folding to diminish the risk of aggregation. Elucidating the molecular mechanisms of chaperones and related enzymes is important for understanding proteostasis systems and protein misfolding- and aggregation-related pathophysiology. Furthermore, mechanistic studies of chaperones and related enzymes provide important clues to designing chemical mimics, or chemical chaperones, that are potentially useful for recovering proteostasis activities as therapeutic approaches for treating and preventing protein misfolding-related diseases. In this Perspective, we provide a comprehensive overview of the latest understanding of the folding-promotion mechanisms by chaperones and oxidoreductases and recent progress in the development of chemical mimics that possess activities comparable to enzymes, followed by a discussion of future directions.

Pathfinder: Protein folding pathway prediction based on conformational sampling

The study of protein folding mechanism is a challenge in molecular biology, which is of great significance for revealing the movement rules of biological macromolecules, understanding the pathogenic mechanism of folding diseases, and designing protein engineering materials. Based on the hypothesis that the conformational sampling trajectory contain the information of folding pathway, we propose a protein folding pathway prediction algorithm named Pathfinder. Firstly, Pathfinder performs large-scale sampling of the conformational space and clusters the decoys obtained in the sampling. The heterogeneous conformations obtained by clustering are named seed states. Then, a resampling algorithm that is not constrained by the local energy basin is designed to obtain the transition probabilities of seed states. Finally, protein folding pathways are inferred from the maximum transition probabilities of seed states. The proposed Pathfinder is tested on our developed test set (34 proteins). For 11 widely studied proteins, we correctly predicted their folding pathways and specifically analyzed 5 of them. For 13 proteins, we predicted their folding pathways to be further verified by biological experiments. For 6 proteins, we analyzed the reasons for the low prediction accuracy. For the other 4 proteins without biological experiment results, potential folding pathways were predicted to provide new insights into protein folding mechanism. The results reveal that structural analogs may have different folding pathways to express different biological functions, homologous proteins may contain common folding pathways, and α-helices may be more prone to early protein folding than β-strands.

Synergetic Effects of Alanine and Glycine in Blob-Based Methods for Predicting Protein Folding Times

The polypeptide PGlyAlaGlu was prepared with 20 mol % glycine (Gly), 36 mol % d,l-alanine (Ala), and 44 mol % d,l-glutamic acid (Glu) and labeled with the dye 1-pyrenemethylamine to yield a series of Py-PGlyAlaGlu samples. The fluorescence decays of the Py-PGlyAlaGlu samples were analyzed according to the fluorescence blob model (FBM) to obtain the number N_blob^exp of amino acids (aa's) encompassed inside the subvolume V_blob of the polypeptide probed by an excited pyrene. An N_blob^exp value of 29 (±2) was retrieved for Py-PGlyAlaGlu, which was much larger than for any of the copolypeptide PGlyGlu or PAlaGlu prepared with either Gly and Glu or Ala and Glu, respectively. The continuous increase in N_blob^exp with decreasing side chain size (SCS) from 10 aa's for PGlu to 16 aa's for PAlaGlu and 22 aa's for PGlyGlu was used earlier to define the reach of an aa and determine the groups of aa's that could interact with each other along a polypeptide backbone according to their SCS. These groups of aa's, referred to as blobs, led to the implementation of blob-based models (BBM) to predict the folding time τ_F^theo,BBM of 145 proteins, which was found to match their experimental folding time τ_F^exp with a relatively high 0.71 correlation coefficient. Nevertheless, the much higher N_blob^exp value found for Py-PGlyAlaGlu compared to all other pyrene-labeled polypeptides studied to date indicates that the reach of aa's along a polypeptide sequence is affected not only by SCS but also by synergetic effects between different aa's. Following this new insight, a revised BBM was implemented to predict τ_F^theo,BBM for 195 proteins assuming the existence or absence of synergies to control the interactions between aa's along a polypeptide sequence. Similarly good correlation coefficients of 0.71 and 0.74 were obtained for a direct 1:1 comparison of τ_F^exp and τ_F^theo,BBM for the 195 proteins without and with synergies, respectively. This result suggests that synergetic effects between different aa's have little effect on τ_F^theo,BBM predicted from BBM underlying the robustness of this methodology.

Polyol and sugar osmolytes stabilize the molten globule state of α-lactalbumin and inhibit amyloid fibril formation

Protein misfolding and aggregation are associated with several human diseases such as Alzheimer's, Parkinson's, prion related disorders, type-II diabetes, etc. Different strategies using molecular chaperones, synthetic and naturally occurring small molecules, osmolytes, etc. have been used to prevent protein aggregation and amyloid fibril formation. In this study, we have used bovine α-lactalbumin at pH 1.6, 37 °C, and shaking conditions to promote amyloid fibril formation. Polyol and sugar osmolytes like glycerol, sorbitol, and trehalose have been used to inhibit the fibrillation of a number of proteins. In the present work, amyloid fibril formation of α-lactalbumin has been shown by ThT assay and AFM, while changes in the secondary structure during fibrillation has been followed by circular dichroism spectroscopy. Our results show that glycerol, sorbitol, and trehalose affect amyloid fibril formation of α-lactalbumin in a concentration-dependent manner. There is a delay in the lag phase of amyloid fibril formation in sorbitol and trehalose and complete inhibition in 6 M glycerol. Our results indicate that delay in the lag phase and inhibition of amyloid fibril formation are due to the stabilization of molten globule state by these osmolytes. At pH 1.6, the molten globule as well as the amyloid fibrils bind to ANS. However, when pH was shifted from 1.6 to 7, only the oligomeric and the fibrillar species bind to ANS due to refolding of molten globule state. The outcome of this study might be useful in designing small molecules which may stabilize the intermediate states, thus preventing amyloid fibril formation.

Effects of the Temperature and Salt Concentration on the Structural Characteristics of the Protein (PDB Code 1BBL)

The effect of the temperature and salt solution on the structural characteristics of the protein 1BBL was investigated by molecular dynamics simulations. The paper presents simulation results regarding the non-bonded energy and the structural stability of the protein immersed in salt solutions with different concentrations and temperatures. Our work demonstrates that the electrostatic potential energy and van der Waals energy of the system show the opposite changes with the influence of the external environment. Since the electrostatic potential energy changes more obviously, it is dominated in the non-bonding interactions. The structural parameters, such as the root mean square deviation and the radius of gyration, increased initially and decreased afterward with the increase of the salt concentration. The protein presented the loose structure with a relative low stability when it was immersed in a monovalent solution with a salt concentration of 0.8 mol/L. The salt concentration corresponding to the maximum value of structural parameters in the monovalent salt solution was double that in the divalent salt solution. It was also concluded that the protein presented a compact and stable structure when immersed in salt solutions with a high concentration of 2.3 mol/L. The analysis of the root mean square deviation and root mean square fluctuation of the protein sample also exhibited that the structural stability and chain flexibility are strongly guided by the effect of the temperature. These conclusions help us to understand the structural characteristics of the protein immersed in the salt solutions with different concentrations and temperatures.

Protein folding in vitro and in the cell: From a solitary journey to a team effort

Correct protein folding is essential for the health and function of living organisms. Yet, it is not well understood how unfolded proteins reach their native state and avoid aggregation, especially within the cellular milieu. Some proteins, especially small, single-domain and apparent two-state folders, successfully attain their native state upon dilution from denaturant. Yet, many more proteins undergo misfolding and aggregation during this process, in a concentration-dependent fashion. Once formed, native and aggregated states are often kinetically trapped relative to each other. Hence, the early stages of protein life are absolutely critical for proper kinetic channeling to the folded state and for long-term solubility and function. This review summarizes current knowledge on protein folding/aggregation mechanisms in buffered solution and within the bacterial cell, highlighting early stages. Remarkably, teamwork between nascent chain, ribosome, trigger factor and Hsp70 molecular chaperones enables all proteins to overcome aggregation propensities and reach a long-lived bioactive state.

Pyrene Excimer Formation (PEF) and Its Application to the Study of Polypeptide Dynamics

This Perspective describes how the fluorescence blob model (FBM) has been developed and applied over the past 30 years to characterize the long-range backbone dynamics (LRBD) of polymers in solution. In these experiments, the polymers are randomly labeled with the dye pyrene, which forms an excimer upon the encounter between an excited and a ground-state pyrenyl label inside a finite subvolume of the polymer coil referred to as a blob representing the volume probed by the excited pyrene. By compartmentalizing the polymer coil into a cluster of identical blobs, FBM analysis of the fluorescence decays acquired with the polymers yields the number N_blob of structural units inside a blob. Since a flexible or rigid backbone will result in an N_blob that is either large or small, N_blob can be used as a measure of the flexibility of a given polymer. After having established that these experiments based on pyrene excimer formation (PEF) yielded quantitative information about the LRBD of a variety of polymers in solution, control experiments were carried out to characterize the effects that different molecular variables, such as the side-chain size (SCS) of a structural unit or the length of the linker connecting pyrene to the polymeric backbone, had on the parameters retrieved with the FBM. At this point, the FBM was applied to study the LRBD of polypeptides prepared from racemic mixtures of amino acids (aa's). These studies led to the establishment of simple rules that could be developed into mathematical equations to describe the LRBD of polypeptides. The N_blob values retrieved from the FBM analysis of the fluorescence decays acquired with the pyrene-labeled polypeptides could then be employed to predict the total conformational search time (τ_tcs) of any polypeptide based on their sequence. Strong correlations were found between the predicted τ_tcs and the experimental folding times of 145 proteins. The good quality of these correlations suggests that the blob-based approach described in this report might represent an interesting mathematical means for studying protein folding.

Cysteine-based protein folding modulators for trapping intermediates and misfolded forms

Folding is a key process to form functional conformations of proteins. Folding via on-pathway intermediates leads to the formation of native structures, while folding through off-pathways affords non-native and disease-causing forms. Trapping folding intermediates and misfolded forms is important for investigating folding mechanisms and disease-related biological properties of the misfolded proteins. We developed cysteine-containing dipeptides conjugated with amino acids possessing mono- and diamino-groups. In oxidative protein folding involving disulfide-bond formation, the addition of cysteine and oxidized glutathione readily promoted the folding to afford native forms. In contrast, despite the acceleration of disulfide-bond formation, non-native isomers formed in significantly increased yields upon the addition of the dipeptides. This study provides a molecular design of cysteine-based protein-folding modulators that afford proteins adopting non-native conformations through intermolecular disulfide-bond formation. Because of the intrinsic reversibility of the disulfide bonds upon redox reactions, the disulfide bond-based approach demonstrated here is expected to lead to the development of reversible methodologies for trapping transient and misfolded forms by intermolecular disulfide bond formation and restarting the folding processes of the trapped forms by subsequent cleavage of the intermolecular disulfide bonds.

In this study, cysteine-containing dipeptides conjugated with amino acids possessing mono- and diamino-groups were developed as protein-folding modulators affording non-native forms through intermolecular disulfide-bond formation.

Effects of an Electric Field on the Conformational Transition of the Protein: Pulsed and Oscillating Electric Fields with Different Frequencies

The effect of pulsed and oscillating electric fields with different frequencies on the conformational properties of all-α proteins was investigated by molecular dynamics simulations. The root mean square deviation, the root mean square fluctuation, the dipole moment distribution, and the secondary structure analysis were used to assess the protein samples’ structural characteristics. In the simulation, we found that the higher frequency of the electric field influences the rapid response to the secondary structural transitions. However, the conformational changes measured by RMSD are diminished by applying the electrical field with a higher frequency. During the dipole moment analysis, we found that the magnitude and frequency of the dipole moment was directly related to the strength and frequency of the external electric field. In terms of the type of electric fields, we found that the average values of RMSD and RMSF of whole molecular protein are larger when the protein is exposed in the pulsed electric field. Concerning the typical sample 1BBL, the secondary structure analysis showed that two alpha-helix segments both transit to turns or random coils almost simultaneously when it is exposed in a pulsed electric field. Meanwhile, two segments present the different characteristic times when the transition occurs in the condition of an oscillating electric field. This study also demonstrated that the protein with fewer charged residues or more residues in forming α-helical structures display the higher conformational stability. These conclusions, achieved using MD simulations, provide a theoretical understanding of the effect of the frequency and expression form of external electric fields on the conformational changes of the all-α proteins with charged residues and the guidance for anticipative applications.

Conjecture on the Design of Helical Proteins

In an important advance in our understanding of protein folding, Wolynes and Onuchic found that the frustration ratio, T_f/T_s, for funneled energy landscapes is T_f/T_s ∼1.6. In our recent work on four heme proteins, we showed that when a protein unfolds from the native state to an early unfolded state, the degree of departure is characterized by a ratio f ∼1.6, where f is a measure of the elongation of n-residue segments of the polypeptide chain. Our analysis, which accounts for this apparent similarity in calculated signatures, is based on a logistic-map model of unfolding. We offer an important take home for the de novo protein synthesis community: in order to increase the probability of obtaining good quality crystals, nearest-neighbor repulsive interactions between adjacent residues (or sequences of residues) in the polypeptide chain must be propagated correctly.

Unfolding cytochromes c-b562 and Rd apo b562

We have analyzed the early stages of unfolding of cytochromes c-b₅₆₂ (PDB ID: 2BC5) and Rd apo b₅₆₂ (PDB ID: 1YYJ). Our geometrical approach proceeds from an analysis of the crystal structure reported for each protein. We quantify, residue-by-residue and region-by-region, the spatial and angular changes in the structure as the protein denatures, and quantify differences that result from the seven residues that differ in the two proteins. Using two independent analyses, one based on spatial metrics and the second on angular metrics, we establish the order of unfolding of the five helices in cyt c-b₅₆₂ and the four helices in the apo protein. For the two helices nearest the N-terminal end of both proteins, the ones in the apo protein unfold first. For the two helices nearest the C-terminal end, the interior helix of the apo protein unfolds first, whereas the terminal helix of the holo protein unfolds first. Excluded-volume effects (repulsive interactions) are minimized in turning regions; the overall range in Δ values is Δ = 36.3 Å³ for cyt c-b₅₆₂ and Δ = 36.6 Å³ for the apo protein, whereas the span for all 20 amino acids is Δ = 167.7 Å³. As our work indicates that the interior helix of cytochrome c-b₅₆₂ is the first to fold, we suggest that this helix protects the heme from misligation, consistent with ultrafast folding over a minimally frustrated funneled landscape.

Funneled angle landscapes for helical proteins

We use crystallographic data for four helical iron proteins (cytochrome c-b₅₆₂, cytochrome c', sperm whale myoglobin, human cytoglobin) to calculate radial and angular signatures as each unfolds from the native state stepwise though four unfolded states. From these data we construct an angle phase diagram to display the evolution of each protein from its native state; and, in turn, the phase diagram is used to construct a funneled angle landscape for comparison with the topography of its folding energy landscape. We quantify the departure of individual helical and turning regions from the areal, angular profile of corresponding regions of the native state. This procedure allows us to identify the similarities and differences among individual helical and turning regions in the early stages of unfolding of the four helical heme proteins.

Effects of an Electric Field on the Conformational Transition of the Protein: A Molecular Dynamics Simulation Study

The effect of the electric field on the conformational properties of the protein 1BBL was investigated by molecular dynamics simulations. Our simulation results clearly capture the structural transitions of the protein sample from helix to turn or random coil conformation induced by the increasing strength of the electric field. During our analysis, we found that the conformational stability is weakened, and the protein sample is stretched as an unfolded structure when it was exposed in a sufficiently high electric field. The characteristic time when the jump occurs in the time evolution curves of root mean square deviation (RMSD) and radius of gyration R_g decreases with increasing electric strength, which demonstrates the rapidly conformational transition that occurs. The number of intra-protein hydrogen bonds, which is the key factor for stabilizing the protein structure, is related to the overall size of the protein. The value of the dipole moment and characteristic time are both influenced by the strength, but are independent of the direction of the external field. The protein sample becomes rotated with the electric field direction. These conclusions provide a theoretical realization of understanding the protein conformational transition in an electric field and the guidance for anticipative applications.

The case for defined protein folding pathways

We consider the differences between the many-pathway protein folding model derived from theoretical energy landscape considerations and the defined-pathway model derived from experiment. A basic tenet of the energy landscape model is that proteins fold through many heterogeneous pathways by way of amino acid-level dynamics biased toward selecting native-like interactions. The many pathways imagined in the model are not observed in the structure-formation stage of folding by experiments that would have found them, but they have now been detected and characterized for one protein in the initial prenucleation stage. Analysis presented here shows that these many microscopic trajectories are not distinct in any functionally significant way, and they have neither the structural information nor the biased energetics needed to select native vs. nonnative interactions during folding. The opposed defined-pathway model stems from experimental results that show that proteins are assemblies of small cooperative units called foldons and that a number of proteins fold in a reproducible pathway one foldon unit at a time. Thus, the same foldon interactions that encode the native structure of any given protein also naturally encode its particular foldon-based folding pathway, and they collectively sum to produce the energy bias toward native interactions that is necessary for efficient folding. Available information suggests that quantized native structure and stepwise folding coevolved in ancient repeat proteins and were retained as a functional pair due to their utility for solving the difficult protein folding problem.

Protein folding by NMR

Protein folding is a highly complex process proceeding through a number of disordered and partially folded nonnative states with various degrees of structural organization. These transiently and sparsely populated species on the protein folding energy landscape play crucial roles in driving folding toward the native conformation, yet some of these nonnative states may also serve as precursors for protein misfolding and aggregation associated with a range of devastating diseases, including neuro-degeneration, diabetes and cancer. Therefore, in vivo protein folding is often reshaped co- and post-translationally through interactions with the ribosome, molecular chaperones and/or other cellular components. Owing to developments in instrumentation and methodology, solution NMR spectroscopy has emerged as the central experimental approach for the detailed characterization of the complex protein folding processes in vitro and in vivo. NMR relaxation dispersion and saturation transfer methods provide the means for a detailed characterization of protein folding kinetics and thermodynamics under native-like conditions, as well as modeling high-resolution structures of weakly populated short-lived conformational states on the protein folding energy landscape. Continuing development of isotope labeling strategies and NMR methods to probe high molecular weight protein assemblies, along with advances of in-cell NMR, have recently allowed protein folding to be studied in the context of ribosome-nascent chain complexes and molecular chaperones, and even inside living cells. Here we review solution NMR approaches to investigate the protein folding energy landscape, and discuss selected applications of NMR methodology to studying protein folding in vitro and in vivo. Together, these examples highlight a vast potential of solution NMR in providing atomistic insights into molecular mechanisms of protein folding and homeostasis in health and disease.

Protein Folding-How and Why: By Hydrogen Exchange, Fragment Separation, and Mass Spectrometry

Advanced hydrogen exchange (HX) methodology can now determine the structure of protein folding intermediates and their progression in folding pathways. Key developments over time include the HX pulse labeling method with nuclear magnetic resonance analysis, the fragment separation method, the addition to it of mass spectrometric (MS) analysis, and recent improvements in the HX MS technique and data analysis. Also, the discovery of protein foldons and their role supplies an essential interpretive link. Recent work using HX pulse labeling with MS analysis finds that a number of proteins fold by stepping through a reproducible sequence of native-like intermediates in an ordered pathway. The stepwise nature of the pathway is dictated by the cooperative foldon unit construction of the protein. The pathway order is determined by a sequential stabilization principle; prior native-like structure guides the formation of adjacent native-like structure. This view does not match the funneled energy landscape paradigm of a very large number of folding tracks, which was framed before foldons were known and is more appropriate for the unguided residue-level search to surmount an initial kinetic barrier rather than for the overall unfolded-state to native-state folding pathway.

Network representation of conformational transitions between hidden intermediates of Rd-apocytochrome b562

The folding kinetics of Rd-apocytochrome b562 is two-state, but native-state hydrogen exchange experiments show that there are discrete partially unfolded (PUF) structures in equilibrium with the native state. These PUF structures are called hidden intermediates because they are not detected in kinetic experiments and they exist after the rate-limiting step. Structures of the mimics of hidden intermediates of Rd-apocytochrome b562 are resolved by NMR. Based upon their relative stability and structural features, the folding mechanism was proposed to follow a specific pathway (unfolded → rate-limiting transition state → PUF1 → PUF2 → native). Investigating the roles of equilibrium PUF structures in folding kinetics and their interrelationship not only deepens our understanding of the details of folding mechanism but also provides guides in protein design and prevention of misfolding. We performed molecular dynamics simulations starting from a hidden intermediate and the native state of Rd-apocytochrome b562 in explicit solvent, for a total of 37.18 μs mainly with Anton. We validated our simulations by detailed comparison with experimental data and other computations. We have verified that we sampled the post rate-limiting transition state region only. Markov state model was used to analyze the simulation results. We replace the specific pathway model with a network model. Transition-path theory was employed to calculate the net effective flux from the most unfolded state towards the most folded state in the network. The proposed sequential folding pathway via PUF1 then more stable, more native-like PUF2 is one of the routes in our network, but it is not dominant. The dominant path visits PUF2 without going through PUF1. There is also a route from PUF1 directly to the most folded state in the network without visiting PUF2. Our results indicate that the PUF states are not necessarily sequential in the folding. The major routes predicted in our network are testable by future experiments such as single molecule experiment.

The Knotted Protein UCH-L1 Exhibits Partially Unfolded Forms under Native Conditions that Share Common Structural Features with Its Kinetic Folding Intermediates

The human ubiquitin C-terminal hydrolase, UCH-L1, is an abundant neuronal deubiquitinase that is associated with Parkinson's disease. It contains a complex Gordian knot topology formed by the polypeptide chain alone. Using a combination of fluorescence-based kinetic measurements, we show that UCH-L1 has two distinct kinetic folding intermediates that are transiently populated on parallel pathways between the denatured and native states. NMR hydrogen-deuterium exchange (HDX) experiments indicate the presence of partially unfolded forms (PUFs) of UCH-L1 under native conditions. HDX measurements as a function of urea concentration were used to establish the structure of the PUFs and pulse-labelled HDX NMR was used to show that the PUFs and the folding intermediates are likely the same species. In both cases, a similar stable core encompassing most of the central β-sheet is highly structured and α-helix 3, which is partially formed, packs against it. In contrast to the stable β-sheet core, the peripheral α-helices display significant local fluctuations leading to rapid exchange. The results also suggest that the main difference between the two kinetic intermediates is structure and packing of α-helices 3 and 7 and the degree of structure in β-strand 5. Together, the fluorescence and NMR results establish that UCH-L1 neither folds through a continuum of pathways nor by a single discrete pathway. Its folding is complex, the β-sheet core forms early and is present in both intermediate states, and the rate-limiting step which is likely to involve the threading of the chain to form the 52-knot occurs late on the folding pathway.

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

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

Cytochrome unfolding pathways from computational analysis of crystal structures

We have developed a model to study the role of geometrical factors in influencing the early stages of unfolding in three cytochromes: cyt c', cyt c-b562 and cyt c. Each stage in unfolding is quantified by the spatial extension λ̂i of n-residue segments, and by their angular extension 〈βn〉. Similarities and differences between and among the three cytochromes in the unfolding of helical and non-helical regions can be determined by analyzing the data for each signature separately. Definite conclusions can be drawn when spatial and angular changes are considered in tandem. To facilitate comparisons, we present graphical portraits of the three cytochromes at the same stage of unfolding, and in relation to their native state structures. We also display specific segments at different stages of unfolding to illustrate differences in stability of defined domains thereby allowing us to make specific predictions on the unfolding of corresponding internal and terminal helices in cyt c' and cyt c-b562. Our work accords with an earlier experimental report on the presence and persistence of a hydrophobic core in cyt c.

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

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

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

Statistical dictionaries for hypothetical in silico model of the early-stage intermediate in protein folding

The polypeptide chain folding process appears to be a multi-stage phenomenon. The scientific community has recently devoted much attention to early stages of this process, with numerous attempts at simulating them—either experimentally or in silico. This paper presents a comparative analysis of the predicted and observed results of folding simulations. The proposed technique, based on statistical dictionaries, yields a global accuracy of 57 %—a marked improvement over older approaches (with an accuracy of approximately 46 %).

Some (dis)assembly required: partial unfolding in the Par-6 allosteric switch

Allostery is commonly described as a functional connection between two distant sites in a protein, where a binding event at one site alters affinity at the other. Here we review the conformational dynamics that encode an allosteric switch in the PDZ domain of Par-6. Par-6 is a scaffold protein that organizes other proteins into a complex required to initiate and maintain cell polarity. NMR measurements revealed that the PDZ domain samples an evolutionarily conserved unfolding intermediate allowing rearrangement of two adjacent loop residues that control ligand binding affinity. Cdc42 binding to Par-6 creates a novel interface between the PDZ domain and the adjoining CRIB motif that stabilizes the high-affinity PDZ conformation. Thermodynamic and kinetic studies suggest that partial PDZ unfolding is an integral part of the Par-6 switching mechanism. The Par-6 CRIB-PDZ module illustrates two important structural aspects of protein evolution: the interface between adjacent domains in the same protein can give rise to allosteric regulation, and thermodynamic stability may be sacrificed to increase the sampling frequency of an unfolding intermediate required for conformational switching.

Reversibility and two state behaviour in the thermal unfolding of oligomeric TIM barrel proteins

The reversible thermal unfolding of oligomeric TIM barrels results from a delicate balance of physicochemical properties related to the sequence, the native and unfolded states and the transition between them.

The nature of protein folding pathways

How do proteins fold, and why do they fold in that way? This Perspective integrates earlier and more recent advances over the 50-y history of the protein folding problem, emphasizing unambiguously clear structural information. Experimental results show that, contrary to prior belief, proteins are multistate rather than two-state objects. They are composed of separately cooperative foldon building blocks that can be seen to repeatedly unfold and refold as units even under native conditions. Similarly, foldons are lost as units when proteins are destabilized to produce partially unfolded equilibrium molten globules. In kinetic folding, the inherently cooperative nature of foldons predisposes the thermally driven amino acid-level search to form an initial foldon and subsequent foldons in later assisted searches. The small size of foldon units, ∼ 20 residues, resolves the Levinthal time-scale search problem. These microscopic-level search processes can be identified with the disordered multitrack search envisioned in the "new view" model for protein folding. Emergent macroscopic foldon-foldon interactions then collectively provide the structural guidance and free energy bias for the ordered addition of foldons in a stepwise pathway that sequentially builds the native protein. These conclusions reconcile the seemingly opposed new view and defined pathway models; the two models account for different stages of the protein folding process. Additionally, these observations answer the "how" and the "why" questions. The protein folding pathway depends on the same foldon units and foldon-foldon interactions that construct the native structure.

Non-native structure appears in microseconds during the folding of E. coli RNase H

The folding pathway of Escherichia coli RNase H is one of the best experimentally characterized for any protein. In spite of this, spectroscopic studies have never captured the earliest events. Using continuous-flow microfluidic mixing, we have now observed the first several milliseconds of folding by monitoring the tryptophan fluorescence lifetime (60 μs dead time). Two folding intermediates are observed, the second of which is the previously characterized I(core) millisecond intermediate. The new earlier intermediate is likely on-pathway and appears to have long-range non-native structure, providing a rare example of such non-native structure formation in a folding pathway. The tryptophan fluorescence lifetimes also suggest a deviation from native packing in the second intermediate, I(core). Similar results from a fragment of RNase H demonstrate that only half of the protein is significantly involved in this early structure formation. These studies give us a view of the formation of tertiary structure on the folding pathway, which complements previous hydrogen-exchange studies that monitored only secondary structure and observed sequential native structure formation. Our results provide detailed folding information on both a timescale and a size-scale accessible to all-atom molecular dynamics simulations of protein folding.

Evidence for close side-chain packing in an early protein folding intermediate previously assumed to be a molten globule

The molten globule, a conformational ensemble with significant secondary structure but only loosely packed tertiary structure, has been suggested to be a ubiquitous intermediate in protein folding. However, it is difficult to assess the tertiary packing of transiently populated species to evaluate this hypothesis. Escherichia coli RNase H is known to populate an intermediate before the rate-limiting barrier to folding that has long been thought to be a molten globule. We investigated this hypothesis by making mimics of the intermediate that are the ground-state conformation at equilibrium, using two approaches: a truncation to generate a fragment mimic of the intermediate, and selective destabilization of the native state using point mutations. Spectroscopic characterization and the response of the mimics to further mutation are consistent with studies on the transient kinetic intermediate, indicating that they model the early intermediate. Both mimics fold cooperatively and exhibit NMR spectra indicative of a closely packed conformation, in contrast to the hypothesis of molten tertiary packing. This result is important for understanding the nature of the subsequent rate-limiting barrier to folding and has implications for the assumption that many other proteins populate molten globule folding intermediates.

Methionine mutations of outer membrane protein X influence structural stability and beta-barrel unfolding

We report the biochemical and biophysical characterization of outer membrane protein X (OmpX), an eight-stranded transmembrane β-barrel from E. coli, and compare the barrel behavior with a mutant devoid of methionine residues. Transmembrane outer membrane proteins of bacterial origin are known to display high tolerance to sequence rearrangements and mutations. Our studies with the triple mutant of OmpX that is devoid of all internal methionine residues (M18L; M21L; M118L) indicate that Met replacement has no influence on the refolding efficiency and structural characteristics of the protein. Surprisingly, the conserved substitution of Met→Leu leads to barrel destabilization and causes a lowering of the unfolding free energy by a factor of ∼8.5 kJ/mol, despite the mutations occurring at the loop regions. We report that the barrel destabilization is accompanied by a loss in cooperativity of unfolding in the presence of chemical denaturants. Furthermore, we are able to detect an unfolding intermediate in the Met-less barrel, whereas the parent protein exhibits a classic two-state unfolding. Thermal denaturation measurements also suggest a greater susceptibility of the OmpX barrel to heat, in the Met-less construct. Our studies reveal that even subtle variations in the extra-membrane region of rigid barrel structures such as OmpX, may bear severe implications on barrel stability. We propose that methionines contribute to efficient barrel structuring and protein-lipid interactions, and are therefore important elements of OmpX stability.

Visualizing side chains of invisible protein conformers by solution NMR

Sparsely populated and transiently formed protein conformers can play key roles in many biochemical processes. Understanding the structure function paradigm requires, therefore, an atomic-resolution description of these rare states. However, they are difficult to study because they cannot be observed using standard biophysical techniques. In the past decade, NMR methods have been developed for structural studies of these elusive conformers, focusing primarily on backbone (1)H, (15)N and (13)C nuclei. Here we extend the methodology to include side chains by developing a (13)C-based chemical exchange saturation transfer experiment for the assignment of side-chain aliphatic (13)C chemical shifts in uniformly (13)C labeled proteins. A pair of applications is provided, involving the folding of β-sheet Fyn SH3 and α-helical FF domains. Over 96% and 89% of the side-chain (13)C chemical shifts for excited states corresponding to the unfolded conformation of the Fyn SH3 domain and a folding intermediate of the FF domain, respectively, have been obtained, providing insight into side-chain packing and dynamics.

Folding of a large protein at high structural resolution

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

Hypothetical in silico model of the early-stage intermediate in protein folding

This paper presents a method for determining the structure of the early stage (ES) intermediate in the multistage protein folding process. ES structure is modeled on the basis of a limited conformational subspace of the Ramachandran plot. The model distinguishes seven structural motifs corresponding to seven local probability maxima within the limited conformational subspace. Three of these are assigned to well-defined secondary structures, while the remaining four are found to represent various types of random coils. The presented heuristic approach also provides insight into the reasons behind incorrect predictions occurring when the folding process depends on external factors (e.g., ligands, ions or other proteins) rather than on the characteristics of the backbone itself. The accuracy of the presented method is estimated at around 48 %.

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

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

A Euclidean Perspective on the Unfolding of Azurin: Spatial Correlations

We investigate the stability to structural perturbation of Pseudomonas aeruginosa azurin using a previously developed geometric model. Our analysis considers Ru(2,2',6',2″-terpyridine)(1,10-phenanthroline)(His83)-labeled wild-type azurin and five variants with mutations to Cu-ligating residues. We find that in the early stages of unfolding, the β-strands exhibit the most structural stability. The conserved residues comprising the hydrophobic core are dislocated only after nearly complete unfolding of the β-barrel. Attachment of the Ru-complex at His83 does not destabilize the protein fold, despite causing some degree of structural rearrangement. Notably, replacing the Cys112 and/or Met121 Cu ligands does not affect the conformational integrity of the protein. Notably, these results are in accord with experimental evidence, as well as molecular dynamics simulations of the denaturation of azurin.

Cold denaturation of a protein dimer monitored at atomic resolution

Protein folding and unfolding are crucial for a range of biological phenomena and human diseases. Defining the structural properties of the involved transient species is therefore of prime interest. Using a combination of cold denaturation with NMR spectroscopy, we reveal detailed insight into the unfolding of the homodimeric repressor protein CylR2. Seven three-dimensional structures of CylR2 at temperatures from 25 °C to -16 °C reveal a progressive dissociation of the dimeric protein into a native-like monomeric intermediate followed by transition into a highly dynamic, partially folded state. The core of the partially folded state seems critical for biological function and misfolding.

Structural analysis and dimerization profile of the SCAN domain of the pluripotency factor Zfp206

Zfp206 (also named as Zscan10) belongs to the subfamily of C₂H₂ zinc finger transcription factors, which is characterized by the N-terminal SCAN domain. The SCAN domain mediates self-association and association between the members of SCAN family transcription factors, but the structural basis and selectivity determinants for complex formation is unknown. Zfp206 is important for maintaining the pluripotency of embryonic stem cells presumably by combinatorial assembly of itself or other SCAN family members on enhancer regions. To gain insights into the folding topology and selectivity determinants for SCAN dimerization, we solved the 1.85 Å crystal structure of the SCAN domain of Zfp206. In vitro binding studies using a panel of 20 SCAN proteins indicate that the SCAN domain Zfp206 can selectively associate with other members of SCAN family transcription factors. Deletion mutations showed that the N-terminal helix 1 is critical for heterodimerization. Double mutations and multiple mutations based on the Zfp206SCAN–Zfp110SCAN model suggested that domain swapped topology is a possible preference for Zfp206SCAN–Zfp110SCAN heterodimer. Together, we demonstrate that the Zfp206SCAN constitutes a protein module that enables C₂H₂ transcription factor dimerization in a highly selective manner using a domain-swapped interface architecture and identify novel partners for Zfp206 during embryonal development.

Transiently populated intermediate functions as a branching point of the FF domain folding pathway

Studies of protein folding and the intermediates that are formed along the folding pathway provide valuable insights into the process by which an unfolded ensemble forms a functional native conformation. However, because intermediates on folding pathways can serve as initiation points of aggregation (implicated in a number of diseases), their characterization assumes an even greater importance. Establishing the role of such intermediates in folding, misfolding, and aggregation remains a major challenge due to their often low populations and short lifetimes. We recently used NMR relaxation dispersion methods and computational techniques to determine an atomic resolution structure of the folding intermediate of a small protein module--the FF domain--with an equilibrium population of 2-3% and a millisecond lifetime, 25 °C. Based on this structure a variant FF domain has been designed in which the native state is selectively destabilized by removing the carboxyl-terminal helix in the native structure to produce a highly populated structural mimic of the intermediate state. Here, we show via solution NMR studies of the designed mimic that the mimic forms distinct conformers corresponding to monomeric and dimeric (K(d) = 0.2 mM) forms of the protein. The conformers exchange on the seconds timescale with a monomer association rate of 1.1 · 10(4) M(-1) s(-1) and with a region responsible for dimerization localized to the amino-terminal residues of the FF domain. This study establishes the FF domain intermediate as a central player in both folding and misfolding pathways and illustrates how incomplete folding can lead to the formation of higher-order structures.

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.

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.

Optimal use of data in parallel tempering simulations for the construction of discrete-state Markov models of biomolecular dynamics

Parallel tempering (PT) molecular dynamics simulations have been extensively investigated as a means of efficient sampling of the configurations of biomolecular systems. Recent work has demonstrated how the short physical trajectories generated in PT simulations of biomolecules can be used to construct the Markov models describing biomolecular dynamics at each simulated temperature. While this approach describes the temperature-dependent kinetics, it does not make optimal use of all available PT data, instead estimating the rates at a given temperature using only data from that temperature. This can be problematic, as some relevant transitions or states may not be sufficiently sampled at the temperature of interest, but might be readily sampled at nearby temperatures. Further, the comparison of temperature-dependent properties can suffer from the false assumption that data collected from different temperatures are uncorrelated. We propose here a strategy in which, by a simple modification of the PT protocol, the harvested trajectories can be reweighted, permitting data from all temperatures to contribute to the estimated kinetic model. The method reduces the statistical uncertainty in the kinetic model relative to the single temperature approach and provides estimates of transition probabilities even for transitions not observed at the temperature of interest. Further, the method allows the kinetics to be estimated at temperatures other than those at which simulations were run. We illustrate this method by applying it to the generation of a Markov model of the conformational dynamics of the solvated terminally blocked alanine peptide.

Probing molecular kinetics with Markov models: metastable states, transition pathways and spectroscopic observables

Markov (state) models (MSMs) have attracted a lot of interest recently as they (1) can probe long-term molecular kinetics based on short-time simulations, (2) offer a way to analyze great amounts of simulation data with relatively little subjectivity of the analyst, (3) provide insight into microscopic quantities such as the ensemble of transition pathways, and (4) allow simulation data to be reconciled with measurement data in a rigorous and explicit way. Here we sketch our current perspective of Markov models and explain in short their theoretical basis and assumptions. We describe transition path theory which allows the entire ensemble of protein folding pathways to be investigated and that combines naturally with Markov models. Experimental observations can be naturally linked to Markov models with the dynamical fingerprint theory, by which experimentally observable timescales can be equipped with an understanding of the structural rearrangement processes that take place at these timescales. The concepts of this paper are illustrated by a simple kinetic model of protein folding.

Intermediates in the protein folding process: a computational model

The paper presents a model for simulating the protein folding process in silico. The two-step model (which consists of the early stage—ES and the late stage—LS) is verified using two proteins, one of which is treated (according to experimental observations) as the early stage and the second as an example of the LS step. The early stage is based solely on backbone structural preferences, while the LS model takes into account the water environment, treated as an external hydrophobic force field and represented by a 3D Gauss function. The characteristics of 1ZTR (the ES intermediate, as compared with 1ENH, which is the LS intermediate) confirm the link between the gradual disappearance of ES characteristics in LS structural forms and the simultaneous emergence of LS properties in the 1ENH protein. Positive verification of ES and LS characteristics in these two proteins (1ZTR and 1ENH respectively) suggest potential applicability of the presented model to in silico protein folding simulations.

Structural gymnastics of multifunctional metamorphic proteins

The classic structure-function paradigm holds that a protein exhibits a single well-defined native state that gives rise to its biological function. Nonetheless, over the past few decades, numerous examples of proteins exhibiting biological function arising from multiple structural states of varying disorder have been identified. Most recently, several examples of 'metamorphic proteins', able to interconvert between vastly different native-like topologies under physiological conditions, have been characterised with multiple functions. In this review, we look at the concept of protein metamorphosis in relation to the current understanding of the protein structure-function landscape. Although structural dynamism observed for metamorphic proteins provides a novel source of functional versatility, the dynamic nature of the metamorphic proteins generally makes them difficult to identify and probe using conventional protein structure determination methods. However, as the existence of metamorphic proteins has now been established and techniques enabling the analysis of multiple protein conformers are improving, it is likely that this class will continue to grow in number.

Nonnative interactions in the FF domain folding pathway from an atomic resolution structure of a sparsely populated intermediate: an NMR relaxation dispersion study

Several all-helical single-domain proteins have been shown to fold rapidly (microsecond time scale) to a compact intermediate state and subsequently rearrange more slowly to the native conformation. An understanding of this process has been hindered by difficulties in experimental studies of intermediates in cases where they are both low-populated and only transiently formed. One such example is provided by the on-pathway folding intermediate of the small four-helix bundle FF domain from HYPA/FBP11 that is populated at several percent with a millisecond lifetime at room temperature. Here we have studied the L24A mutant that has been shown previously to form nonnative interactions in the folding transition state. A suite of Carr-Purcell-Meiboom-Gill relaxation dispersion NMR experiments have been used to measure backbone chemical shifts and amide bond vector orientations of the invisible folding intermediate that form the input restraints in calculations of atomic resolution models of its structure. Despite the fact that the intermediate structure has many features that are similar to that of the native state, a set of nonnative contacts is observed that is even more extensive than noted previously for the wild-type (WT) folding intermediate. Such nonnative interactions, which must be broken prior to adoption of the native conformation, explain why the transition from the intermediate state to the native conformer (millisecond time scale) is significantly slower than from the unfolded ensemble to the intermediate and why the L24A mutant folds more slowly than the WT.