Populated Folding Intermediates Research Articles

Mapping the FF domain folding pathway via structures of transiently populated folding intermediates

Despite the tremendous accomplishments of AlpaFold2/3 in predicting biomolecular structure, the protein folding problem remains unsolved in the sense that accurate atomistic models of how protein molecules fold into their native conformations from an unfolded ensemble are still elusive. Here, using chemical exchange saturation transfer (CEST) NMR experiments and a comprehensive four-state kinetic model of the folding trajectory of a 71 residue four-helix bundle FF domain from human HYPA/FBP11 we present an atomic resolution structure of a transiently formed intermediate, I2, that along with the structure of a second intermediate, I1, provides a description of the FF domain folding trajectory. By recording CEST profiles as a function of urea concentration the extent of compaction along the folding pathway is evaluated. Our data establish that unlike the partially disordered I1 state, the I2 intermediate that is also formed before the rate-limiting folding barrier is well ordered and compact like the native conformer, while retaining nonnative interactions similar to those found in I1. The slow-interconversion from I2 to F, involving changes in secondary structure and the breaking of nonnative interactions, proceeds via a compact transition-state. Interestingly, the native state of the FF1 domain from human p190-A Rho GAP resembles the I2 conformation, suggesting that well-ordered folding intermediates can be repurposed by nature in structurally related proteins to assume functional roles. It is anticipated that the strategy for elucidation of sparsely populated and transiently formed structures of intermediates along kinetic pathways described here will be of use in other studies of protein dynamics.

Optimized fast mixing device for real-time NMR applications

We present an improved fast mixing device based on the rapid mixing of two solutions inside the NMR probe, as originally proposed by Hore and coworkers (J. Am. Chem. Soc. 125 (2003) 12484–12492). Such a device is important for off-equilibrium studies of molecular kinetics by multidimensional real-time NMR spectrsocopy. The novelty of this device is that it allows removing the injector from the NMR detection volume after mixing, and thus provides good magnetic field homogeneity independently of the initial sample volume placed in the NMR probe. The apparatus is simple to build, inexpensive, and can be used without any hardware modification on any type of liquid-state NMR spectrometer. We demonstrate the performance of our fast mixing device in terms of improved magnetic field homogeneity, and show an application to the study of protein folding and the structural characterization of transiently populated folding intermediates.

Folding analysis of the most complex Stevedore's protein knot.

DehI is a homodimeric haloacid dehalogenase from Pseudomonas putida that contains the most complex 61 Stevedore’s protein knot within its folding topology. To examine how DehI attains such an intricate knotted topology we combined far-UV circular dichroism (CD), intrinsic fluorescence spectroscopy and small angle X-ray scattering (SAXS) to investigate its folding mechanism. Equilibrium unfolding of DehI by chemical denaturation indicated the presence of two highly populated folding intermediates, I and I’. While the two intermediates vary in secondary structure contents and tertiary packing according to CD and intrinsic fluorescence, respectively, their overall dimension and compactness are similar according to SAXS. Three single-tryptophan variants (W34, W53, and W196) were generated to probe non-cooperative unfolding events localized around the three fluorophores. Kinetic fluorescence measurements indicated that the transition from the intermediate I’ to the unfolded state is rate limiting. Our multiparametric folding analyses suggest that DehI unfolds through a linear folding pathway with two distinct folding intermediates by initial hydrophobic collapse followed by nucleation condensation, and that knotting precedes the formation of secondary structures.

Native state dynamics affects the folding transition of porcine pancreatic phospholipase A2

Porcine pancreatic phospholipase A2, a small and disulfide rich protein, is extremely resistant against chemically or thermally induced unfolding. Despite this marked resistance, the protein displays broad unfolding transitions resulting in comparatively low apparent thermodynamic stability. Broad unfolding transitions may result from undetected folding intermediates, residual structures in the unfolded state or an inhomogeneity of the native state. Using circular dichroism, fluorescence, and NMR spectroscopy, we ruled out the existence of stably populated folding intermediates, whereas UV absorbance measurements hinted at stable residual structures in the unfolded state. These residual structures proved, however, to have no impact on the folding parameters. Studies by limited proteolysis, CD, and NMR spectroscopy under non-denaturing conditions suggested pronounced dynamics of the protein in the native state, which as long as unrestrained by acidic pH or bound Ca2+ ions exert considerable influence on the unfolding transition.

Under‐folded proteins: Conformational ensembles and their roles in protein folding, function, and pathogenesis

ABSTRACTFor decades, protein function was intimately linked to the presence of a unique, aperiodic crystal‐like structure in a functional protein. The two only places for conformational ensembles of under‐folded (or partially folded) protein forms in this picture were either the end points of the protein denaturation processes or transiently populated folding intermediates. Recent years witnessed dramatic change in this perception and conformational ensembles, which the under‐folded proteins are, have moved from the shadow. Accumulated to date data suggest that a protein can exist in at least three global forms–functional and folded, functional and intrinsically disordered (nonfolded), and nonfunctional and misfolded/aggregated. Under‐folded protein states are crucial for each of these forms, serving as important folding intermediates of ordered proteins, or as functional states of intrinsically disordered proteins (IDPs) and IDP regions (IDPRs), or as pathology triggers of misfolded proteins. Based on these observations, conformational ensembles of under‐folded proteins can be classified as transient (folding and misfolding intermediates) and permanent (IDPs and stable misfolded proteins). Permanently under‐folded proteins can further be split into intentionally designed (IDPs and IDPRs) and unintentionally designed (misfolded proteins). Although intrinsic flexibility, dynamics, and pliability are crucial for all under‐folded proteins, the different categories of under‐foldedness are differently encoded in protein amino acid sequences. © 2013 Wiley Periodicals, Inc. Biopolymers 99: 870–887, 2013.

Toward a Digital Gene Response: RNA G-Quadruplexes with Fewer Quartets Fold with Higher Cooperativity

Changes in RNA conformation can alter gene expression. The guanine quadruplex sequence (GQS) is an RNA motif that folds in the presence of K(+) ions. Changes in the conformation of this motif could be especially important in regulating gene expression in plants because intracellular K(+) concentrations often increase during drought stress. Little is known about the folding thermodynamics of RNA GQS. We show here that RNA GQS with tracts containing three G's [e.g., (GGGxx)(4)] have a modest dependence on the K(+) concentration, folding with no or even negative cooperativity (Hill coefficients ≤1), and are associated with populated folding intermediates. In contrast, GQS with tracts containing just two G's [e.g., (GGxx)(4)] have a steep dependence on the K(+) concentration and fold with positive cooperativity (Hill coefficients of 1.7-2.7) without significantly populating intermediate states. We postulate that in plants, the more stable G3 sequences are largely folded even under unstressed conditions, while the less stable G2 sequences fold only at the higher K(+) concentrations associated with cellular stress, wherein they respond sharply to changing K(+) concentrations. Given the binary nature of their folding, G2 sequences may find application in computation with DNA and in engineering of genetic circuits.

Sparsely populated folding intermediates of the Fyn SH3 domain: matching native-centric essential dynamics and experiment.

A complete description of how a protein folds requires the characterization of intermediate conformations traversed during the folding transition. We have calculated dynamics trajectories of a simplified model of the Fyn SH3 domain with a native-centric potential energy function. Analysis of the resulting site-resolved energy trajectory identifies an ensemble of intermediate conformations for folding and another for unfolding. The model's folding intermediate is structured in the three beta-strands that make up the protein's core and is strikingly similar to intermediates detected in a recent NMR study of Fyn SH3 folding and to folding transition states elucidated in mutagenesis studies of SH3 domains. The unfolding intermediate is formed by dissociation of the folded protein's two terminal beta-strands from its core. The presence of such an intermediate is consistent with the results of a protein-engineering study on the src SH3 domain showing that these strands separate before the rate-limiting step of unfolding. Despite the presence of these conformations intermediate between the native and fully unfolded states, the computed heat capacity vs. temperature profile of the model protein indicates that its thermodynamics satisfies the usual calorimetric criterion for two-state folding. This observation highlights the fact that, if not properly interpreted, methods such as calorimetry that do not probe multiple sites in a molecule can lead to an oversimplified view of folding. The close agreement between results from this simplified model and experimental work underscores the important contributions that computational methods can make in providing insights into protein folding.

An amino acid code for protein folding.

Direct structural information obtained for many proteins supports the following conclusions. The amino acid sequences of proteins can stabilize not only the final native state but also a small set of discrete partially folded native-like intermediates. Intermediates are formed in steps that use as units the cooperative secondary structural elements of the native protein. Earlier intermediates guide the addition of subsequent units in a process of sequential stabilization mediated by native-like tertiary interactions. The resulting stepwise self-assembly process automatically constructs a folding pathway, whether linear or branched. These conclusions are drawn mainly from hydrogen exchange-based methods, which can depict the structure of infinitesimally populated folding intermediates at equilibrium and kinetic intermediates with subsecond lifetimes. Other kinetic studies show that the polypeptide chain enters the folding pathway after an initial free-energy-uphill conformational search. The search culminates by finding a native-like topology that can support forward (native-like) folding in a free-energy-downhill manner. This condition automatically defines an initial transition state, the search for which sets the maximum possible (two-state) folding rate. It also extends the sequential stabilization strategy, which depends on a native-like context, to the first step in the folding process. Thus the native structure naturally generates its own folding pathway. The same amino acid code that translates into the final equilibrium native structure-by virtue of propensities, patterning, secondary structural cueing, and tertiary context-also produces its kinetic accessibility.

Protein Folding: Matching Theory and Experiment

The impact of folding funnels and folding simulations on the way experimentalists interpret results is examined. The image of the transition state has changed from a unique species that has a strained configuration, with a correspondingly high free energy, to a more ordinary folding intermediate, whose balance between limited conformational entropy and stabilizing contacts places it at the top of the free energy barrier. Evidence for a broad transition barrier comes from studies showing that mutations can change the position of the barrier. The main controversial issue now is whether populated folding intermediates are productive on-pathway intermediates or dead-end traps. Direct experimental evidence is needed. Theories suggesting that populated intermediates are trapped in a glasslike state are usually based on mechanisms which imply that trapping would only be extremely short-lived (e.g., nanoseconds) in water at 25°C. There seems to be little experimental evidence for long-lived trapping in monomers, if folding aggregates are excluded. On the other hand, there is good evidence for kinetic trapping in dimers. α-Helix formation is currently the fastest known process in protein folding, and incipient helices are present at the start of folding. Fast helix formation has the effect of narrowing drastically the choice of folding routes. Thus helix formation can direct folding. It changes the folding metaphor from pouring liquid down a folding funnel to a train leaving a switchyard with only a few choices of exit tracks.

The Disulfide Folding Pathway of Human Epidermal Growth Factor

Human epidermal growth factor (EGF) contains three disulfides and 53 amino acids. Reduced/denatured EGF refolds spontaneously in vitro to acquire its native structure. The mechanism of this folding process has been elucidated by structural analysis of both acid and iodoacetate trapped intermediates. The results reveal that the folding is accompanied by a sequential flow of unfolded EGF (0-disulfide) through three groups of folding intermediates, namely 1-disulfide, 2-disulfide, and 3-disulfide (scrambled) EGF isomers, to reach the native structure. Equilibrium occurs among isomers of each class of disulfide species, and the composition of intermediates appears to be highly heterogeneous. Together, at least 27 fractions of folding intermediates have been identified, but there exist only limited numbers of well populated species which constitute more than 80% of the total intermediates found during EGF folding. Six species of such well populated intermediates have been isolated, which included two 1-S-S, two 2-S-S, and two 3-S-S scrambled species. Their disulfide structures have been identified here. Both 1-S-S isomers are found to contain non-native disulfides. One of the 2-S-S species consists of two non-native disulfides and the other admits two native disulfides. Among the six disulfides of the two scrambled species, only one is native. Together, native disulfides constitute 25% of the total disulfides found in these six well populated intermediates. These results contrast sharply to those observed with bovine pancreatic trypsin inhibitor, which has shown that well populated folding intermediates consist of exclusively native disulfides (Weissman, J. S., and Kim, P. S. (1991) Science 253, 1386-1393). We propose that well populated folding intermediates, regardless of whether they contain native or non-native disulfides, do not necessarily represent the productive species and specify the folding pathway. Furthermore, conditions influencing the efficiency of EGF folding have been investigated. It is demonstrated here that under optimized compositions of redox agents, including the use of cysteine/cystine and protein disulfide isomerase, the in vitro folding of EGF could be achieved quantitatively within 1 min.

A reexamination of the conformational transitions of T4 thioredoxin

The unfolding and refolding of T4 thioredoxin was observed by equilibrium and kinetic size exclusion chromatographic measurements in guanidine hydrochloride at 4 degrees C and pH 7.0. All the observed chromatographic profiles can be simulated by a cubic mechanism using a consistent set of equilibrium and kinetic parameters describing each of the coupled transitions. The four components in the folded protein and in the unfolded protein are interrelated by configurational transitions having parameters characteristic for proline peptide isomerizations. Only two of the four folded conformations are significantly populated at equilibrium. Each of the four unfolded components can refold by a unique conformational transition. No transiently populated folding intermediates are detected having hydrodynamic volumes intermediate between those characteristics for the folded and unfolded protein.

Recombination of S-peptide with S-protein during folding of ribonuclease S: I. Folding pathways of the slow-folding and fast-folding classes of unfolded S-protein

The refolding kinetics of ribonuclease S have been measured by tyrosine absorbance, by tyrosine fluorescence emission, and by rapid binding of the specific inhibitor 2′CMP † † Abbreviations used: 2′CMP, cytidine 2′-phosphate; nuclease, staphylococcal nuclease; RNAase, bovine pancreatic ribonuclease, disulfide bonds intact; GuHCl, guanidinium hydrochloride τ, time constant of kinetic phase or reaction defined as the reciprocal of the apparent rate constant; tm, temperature midpoint of a thermal folding transition. to folded RNAase S. The S-protein is first unfolded at pH 1.7 and then either mixed with S-peptide as refolding is initiated by a stopped-flow pH jump to pH 6.8, or the same results are obtained if S-protein and S-peptide are present together before refolding is initiated. The refolding kinetics of RNAase S have been measured as a function of temperature (10 to 40 °C) and of protein concentration (10 to 120 μ m). The results are compared to the folding kinetics of S-protein alone and to earlier studies of RNAase A. A thermal folding transition of S-protein has been found below 30 °C at pH 1.7; its effects on the refolding kinetics are described in the following paper (Labhardt & Baldwin, 1979). In this paper we characterize the refolding kinetics of unfolded S-protein, as it is found above 30 °C at pH 1.7, together with the kinetics of combination between S-peptide and S-protein during folding at pH 6.8. Two classes of unfolded S-protein molecules are found, fast-folding and slow-folding molecules, in a 20: 80 ratio. This is the same result as that found earlier for RNAase A; it is expected if the slow-folding molecules are produced by the slow cis-trans isomerization of proline residues after unfolding, since S-protein contains all four proline residues of RNAase A. The refolding kinetics of the fast-folding molecules show clearly that combination between S-peptide and S-protein occurs before folding of S-protein is complete. If combination occurred only after complete folding, then the kinetics of formation of RNAase S should be rather slow (5 s and 100 s at 30 °C) and nearly independent of protein concentration, as shown by separate measurements of the folding kinetics of S-protein, and of the combination between S-peptide and folded S-protein. The observed folding kinetics are faster than predicted by this model and also the folding rate increases strongly with protein concentration (apparent 1.6 order kinetics). The fact that RNAase S is formed more rapidly than S-protein alone is sufficient by itself to show that combination with S-peptide precedes complete folding of S-protein. Computer simulation of a simple, parallel-pathway scheme is able to reproduce the folding kinetics of the fast-folding molecules. All three probes give the same folding kinetics. These results exclude the model for protein folding in which the rate-limiting step is an initial diffusion of the polypeptide chain into a restricted range of three-dimensional configurations (“nueleation”) followed by rapid folding (“propagation”). If this model were valid, one would expect comparable rates of folding for RNAase A and for S-protein and one would also expect to find no populated folding intermediates, so that combination between S-peptide and S-protein should occur after folding is complete. Instead, RNAase A folds 60 times more rapidly than S-protein and also combination with S-peptide occurs before folding of S-protein is complete. The results demonstrate that the folding rate of S-protein increases after the formation, or stabilization, of an intermediate which results from combination with S-peptide. They support a sequential model for protein folding in which the rates of successive steps in folding depend on the stabilities of preceding intermediates. The refolding kinetics of the slow-folding molecules are complex. Two results demonstrate the presence of folding intermediates: (1) the three probes show different kinetic progress curves, and (2) the folding kinetics are concentration-dependent, in contrast to the results expected if complete folding of S-protein precedes combination with S-peptide. A faster phase of the slow-refolding reaction is detected both by tyrosine absorbance and fluorescence emission but not by 2′CMP binding, indicating that native RNAase S is not formed in this phase. Comparison of the kinetic progress curves measured by different probes is made with the use of the kinetic ratio test, which is defined here.