Folding intermediates of wild-type and mutants of barnase. II. Correlation of changes in equilibrium amide exchange kinetics with the population of the folding intermediate

From covalent transition states in chemistry to noncovalent in biology: from β- to Φ-value analysis of protein folding

Solving the mechanism of a chemical reaction requires determining the structures of all the ground states on the pathway and the elusive transition states linking them. 2024 is the centenary of Brønsted's landmark paper that introduced the β-value and structure-activity studies as the only experimental means to infer the structures of transition states. It involves making systematic small changes in the covalent structure of the reactants and analysing changes in activation and equilibrium-free energies. Protein engineering was introduced for an analogous procedure, Φ-value analysis, to analyse the noncovalent interactions in proteins central to biological chemistry. The methodology was developed first by analysing noncovalent interactions in transition states in enzyme catalysis. The mature procedure was then applied to study transition states in the pathway of protein folding - 'part (b) of the protein folding problem'. This review describes the development of Φ-value analysis of transition states and compares and contrasts the interpretation of β- and Φ-values and their limitations. Φ-analysis afforded the first description of transition states in protein folding at the level of individual residues. It revealed the nucleation-condensation folding mechanism of protein domains with the transition state as an expanded, distorted native structure, containing little fully formed secondary structure but many weak tertiary interactions. A spectrum of transition states with various degrees of structural polarisation was then uncovered that spanned from nucleation-condensation to the framework mechanism of fully formed secondary structure. Φ-analysis revealed how movement of the expanded transition state on an energy landscape accommodates the transition from framework to nucleation-condensation mechanisms with a malleability of structure as a unifying feature of folding mechanisms. Such movement follows the rubric of analysis of classical covalent chemical mechanisms that began with Brønsted. Φ-values are used to benchmark computer simulation, and Φ and simulation combine to describe folding pathways at atomic resolution.

The Molten Globule, and Two-State vs. Non-Two-State Folding of Globular Proteins

From experimental studies of protein folding, it is now clear that there are two types of folding behavior, i.e., two-state folding and non-two-state folding, and understanding the relationships between these apparently different folding behaviors is essential for fully elucidating the molecular mechanisms of protein folding. This article describes how the presence of the two types of folding behavior has been confirmed experimentally, and discusses the relationships between the two-state and the non-two-state folding reactions, on the basis of available data on the correlations of the folding rate constant with various structure-based properties, which are determined primarily by the backbone topology of proteins. Finally, a two-stage hierarchical model is proposed as a general mechanism of protein folding. In this model, protein folding occurs in a hierarchical manner, reflecting the hierarchy of the native three-dimensional structure, as embodied in the case of non-two-state folding with an accumulation of the molten globule state as a folding intermediate. The two-state folding is thus merely a simplified version of the hierarchical folding caused either by an alteration in the rate-limiting step of folding or by destabilization of the intermediate.

PFDB: A standardized protein folding database with temperature correction

We constructed a standardized protein folding kinetics database (PFDB) in which the logarithmic rate constants of all listed proteins are calculated at the standard temperature (25 °C). A temperature correction based on the Eyring–Kramers equation was introduced for proteins whose folding kinetics were originally measured at temperatures other than 25 °C. We verified the temperature correction by comparing the logarithmic rate constants predicted and experimentally observed at 25 °C for 14 different proteins, and the results demonstrated improvement of the quality of the database. PFDB consists of 141 (89 two-state and 52 non-two-state) single-domain globular proteins, which has the largest number among the currently available databases of protein folding kinetics. PFDB is thus intended to be used as a standard for developing and testing future predictive and theoretical studies of protein folding. PFDB can be accessed from the following link: http://lee.kias.re.kr/~bala/PFDB.

Mechanical Folding and Unfolding of Protein Barnase at the Single-Molecule Level

The unfolding and folding of protein barnase has been extensively investigated in bulk conditions under the effect of denaturant and temperature. These experiments provided information about structural and kinetic features of both the native and the unfolded states of the protein, and debates about the possible existence of an intermediate state in the folding pathway have arisen. Here, we investigate the folding/unfolding reaction of protein barnase under the action of mechanical force at the single-molecule level using optical tweezers. We measure unfolding and folding force-dependent kinetic rates from pulling and passive experiments, respectively, and using Kramers-based theories (e.g., Bell-Evans and Dudko-Hummer-Szabo models), we extract the position of the transition state and the height of the kinetic barrier mediating unfolding and folding transitions, finding good agreement with previous bulk measurements. Measurements of the force-dependent kinetic barrier using the continuous effective barrier analysis show that protein barnase verifies the Leffler-Hammond postulate under applied force and allow us to extract its free energy of folding, ΔG0. The estimated value of ΔG0 is in agreement with our predictions obtained using fluctuation relations and previous bulk studies. To address the possible existence of an intermediate state on the folding pathway, we measure the power spectrum of force fluctuations at high temporal resolution (50 kHz) when the protein is either folded or unfolded and, additionally, we study the folding transition-path time at different forces. The finite bandwidth of our experimental setup sets the lifetime of potential intermediate states upon barnase folding/unfolding in the submillisecond timescale.

Partially Unfolded Forms of the Prion Protein Populated under Misfolding-promoting Conditions: CHARACTERIZATION BY HYDROGEN EXCHANGE MASS SPECTROMETRY AND NMR

The susceptibility of the cellular prion protein (PrP(C)) to convert to an alternative misfolded conformation (PrP(Sc)), which is the key event in the pathogenesis of prion diseases, is indicative of a conformationally flexible native (N) state. In the present study, hydrogen-deuterium exchange (HDX) in conjunction with mass spectrometry and nuclear magnetic resonance spectroscopy were used for the structural and energetic characterization of the N state of the full-length mouse prion protein, moPrP(23-231), under conditions that favor misfolding. The kinetics of HDX of 34 backbone amide hydrogens in the N state were determined at pH 4. In contrast to the results of previous HDX studies on the human and Syrian hamster prion proteins at a higher pH, various segments of moPrP were found to undergo different extents of subglobal unfolding events at pH 4, a pH at which the protein is known to be primed to misfold to a β-rich conformation. No residual structure around the disulfide bond was observed for the unfolded state at pH 4. The N state of the prion protein was observed to be at equilibrium with at least two partially unfolded forms (PUFs). These PUFs, which are accessed by stochastic fluctuations of the N state, have altered surface area exposure relative to the N state. One of these PUFs resembles a conformation previously implicated to be an initial intermediate in the conversion of monomeric protein into misfolded oligomer at pH 4.

The partially folded homodimeric intermediate of Escherichia coli aspartate aminotransferase contains a "molten interface" structure

The role of intersubunit side chain-side chain interactions in the stability of the Escherichia coli aspartate aminotransferase (eAATase) homodimer was investigated by directed mutagenesis at 10 different interface contacts. The urea-mediated unfolding pathway of this enzyme proceeds through the formation of a dimeric intermediate, D*, that retains only 40% of the native enzyme secondary structure as judged by circular dichroism. Disruption of any single intersubunit interaction results in a >2.6 kcal mol(-1) decrease in native state stability, independent of its location or nature. However, the stability of D* with respect to U, the unfolded monomer, is the same for all mutants. The stability of the eAATase interface cannot be ascribed to the contribution of a few hot spots, or to the accumulation of a large number of weak interactions, but only to the presence of multiple important and interconnected interactions. It is proposed that a "molten interface" structure, flexible enough to accommodate point mutations, accounts for the stability of D*. Nuclei of tertiary structure, which are not involved in native intersubunit contacts, likely provide a scaffold for the unstructured interface of D*. Such a scaffold would account for the cooperative unfolding of the intermediate.

Effects of co-operative ligand binding on protein amide NH hydrogen exchange

Amide protection factors have been determined from NMR measurements of hydrogen/deuterium amide NH exchange rates measured on assigned signals from Lactobacillus casei apo-DHFR and its binary and ternary complexes with trimethoprim (TMP), folinic acid and coenzymes (NADPH/NADP(+)). The substantial sizes of the residue-specific DeltaH and TDeltaS values for the opening/closing events in NH exchange for most of the measurable residues in apo-DHFR indicate that sub-global or global rather than local exchange mechanisms are usually involved. The amide groups of residues in helices and sheets are those most protected in apo-DHFR and its complexes, and the protection factors are generally related to the tightness of ligand binding. The effects of ligand binding that lead to changes in amide protection are not localised to specific binding sites but are spread throughout the structure via a network of intramolecular interactions. Although the increase in protein stability in the DHFR.TMP.NADPH complex involves increased ordering in the protein structure (requiring TDeltaS energy) this is recovered, to a large extent, by the stronger binding (enthalpic DeltaH) interactions made possible by the reduced motion in the protein. The ligand-induced protection effects in the ternary complexes DHFR.TMP.NADPH (large positive binding co-operativity) and DHFR.folinic acid.NADPH (large negative binding co-operativity) mirror the co-operative effects seen in the ligand binding. For the DHFR.TMP.NADPH complex, the ligand-induced protection factors result in DeltaDeltaG(o) values for many residues being larger than the DeltaDeltaG(o) values in the corresponding binary complexes. In contrast, for DHFR.folinic acid.NADPH, the DeltaDeltaG(o) values are generally smaller than many of those in the corresponding binary complexes. The results indicate that changes in protein conformational flexibility on formation of the ligand complex play an important role in determining the co-operativity in the ligand binding.

Stabilization of Native Protein Fold by Intein-Mediated Covalent Cyclization

A mutant version of the N-terminal domain of Escherichia coli DnaB helicase was used as a model system to assess the stabilization against unfolding gained by covalent cyclization. Cyclization was achieved in vivo by formation of an amide bond between the N and C termini with the help of a split mini-intein. Linear and circular proteins were constructed to be identical in amino acid sequence. Mutagenesis of Phe102 to Glu rendered the protein monomeric even at high concentration. A difference in free energy of unfolding, DeltaDeltaG, between circular and linear protein of 2.3(+/-0.5) kcal mol(-1) was measured at 10 degrees C by circular dichroism. A theoretical estimate of the difference in conformational entropy of linear and circular random chains in a three-dimensional cubic lattice model predicted DeltaDeltaG=2.3 kcal mol(-1), suggesting that stabilization by protein cyclization is driven by the reduced conformational entropy of the unfolded state. Amide-proton exchange rates measured by NMR spectroscopy and mass spectrometry showed a uniform, approximately tenfold decrease of the exchange rates of the most slowly exchanging amide protons, demonstrating that cyclization globally decreases the unfolding rate of the protein. The amide proton exchange was found to follow EX1 kinetics at near-neutral pH, in agreement with an unusually slow refolding rate of less than 4 min(-1) measured by stopped-flow circular dichroism. The linear and circular proteins differed more in their unfolding than in their folding rates. Global unfolding of the N-terminal domain of E.coli DnaB is thus promoted strongly by spatial separation of the N and C termini, whereas their proximity is much less important for folding.

Non-linear effects of temperature and urea on the thermodynamics and kinetics of folding and unfolding of hisactophilin

Extensive measurements and analysis of thermodynamic stability and kinetics of urea-induced unfolding and folding of hisactophilin are reported for 5-50 degrees C, at pH 6.7. Under these conditions hisactophilin has moderate thermodynamic stability, and equilibrium and kinetic data are well fit by a two-state transition between the native and the denatured states. Equilibrium and kinetic m values decrease with increasing temperature, and decrease with increasing denaturant concentration. The betaF values at different temperatures and urea concentrations are quite constant, however, at about 0.7. This suggests that the transition state for hisactophilin unfolding is native-like and changes little with changing solution conditions, consistent with a narrow free energy profile for the transition state. The activation enthalpy and entropy of unfolding are unusually low for hisactophilin, as is also the case for the corresponding equilibrium parameters. Conventional Arrhenius and Eyring plots for both folding and unfolding are markedly non-linear, but these plots become linear for constant DeltaG/T contours. The Gibbs free energy changes for structural changes in hisactophilin have a non-linear denaturant dependence that is comparable to non-linearities observed for many other proteins. These non-linearities can be fit for many proteins using a variation of the Tanford model, incorporating empirical quadratic denaturant dependencies for Gibbs free energies of transfer of amino acid constituents from water to urea, and changes in fractional solvent accessible surface area of protein constituents based on the known protein structures. Noteworthy exceptions that are not well fit include amyloidogenic proteins and large proteins, which may form intermediates. The model is easily implemented and should be widely applicable to analysis of urea-induced structural transitions in proteins.

Accelerated design of bioconversion processes using automated microscale processing techniques

Microscale processing techniques are rapidly emerging as a means to increase the speed of bioprocess design and reduce material requirements. Automation of these techniques can reduce labour intensity and enable a wider range of process variables to be examined. This article examines recent research on various individual microscale unit operations including microbial fermentation, bioconversion and product recovery techniques. It also explores the potential of automated whole process sequences operated in microwell formats. The power of the whole process approach is illustrated by reference to a particular bioconversion, namely the Baeyer-Villiger oxidation of bicyclo[3.2.0]hept-2-en-6-one for the production of optically pure lactones.

Directed evolution of barnase stability using proteolytic selection

We report the construction of a phage-displayed repertoire of mutants of the ribonuclease barnase from Bacillus amyloliquefaciens. The construction was guided by the natural variability between two closely related ribonucleases, barnase and binase from Bacillus intermedius. This repertoire was selected using a proteolytic selection method, allowing sorting of the library according to the resistance of the mutants toward proteolysis. Susceptibility toward proteolysis has been correlated with flexibility and unfolding, and is thus expected to yield mutants with increased thermal stability. Enrichment of barnase mutants with specific combinations of amino acid residues at four of the randomised positions was observed. Three of these enriched amino acid residues are present in neither barnase nor binase. For some of the mutations, the improvement in proteolytic stability does not lead to a pronounced improvement in thermodynamic stability, indicating that the factors governing the proteolytic stability in some cases may be different from those governing the thermodynamic stability, e.g. propensity to local unfolding.The results obtained add important knowledge to a novel use of phage display technology for selection of thermodynamically stable proteins. Only by carefully establishing the parameters that can be adjusted, and recognising the influence this will have on the outcome of selection, will it be possible to realise the powerful technique of proteolytic selection.

Lack of definable nucleation sites in the rate-limiting transition state of barnase under native conditions

It has been shown that the burst-phase (submillisecond) intermediate of barnase, if it exists, can be only marginally more stable than the fully unfolded state at pH 6.3 and 25 degrees C. In the study reported here, no stable burst-phase intermediate could be detected, even in the presence of stabilizing salt (0.4 M Na(2)SO(4)). These results suggest that a burst-phase intermediate with even marginal stability does not exist. The absence of such an intermediate in turn suggests the need for re-examination of the rate-limiting transition state (RLTS) under native conditions, which was previously characterized by using a three-state model with a stable intermediate and protein engineering. Surprisingly, mutations throughout the structure of barnase do not significantly affect the folding rate, suggesting a lack of specific favorable interactions among the side-chains in the RLTS. This RLTS is clearly different from that previously characterized under denaturing conditions, indicating that changes take place in the RLTS under native and denaturing conditions. The occurrence of such changes is further supported by the observation that the unfolding rate constants of barnase and its mutants were divergent or convergent as a function of denaturant concentrations. Consistent with changes in the RLTS, a re-analysis of data from native-state hydrogen exchange studies has shown that the logarithm of the unfolding rate constant inflects down under low concentrations of denaturant. Here, we discuss in detail the question of whether changes in the RLTS involve a kinetically silent intermediate that occurs after the initial RLTS.

Four-body potentials reveal protein-specific correlations to stability changes caused by hydrophobic core mutations

Mutational experiments show how changes in the hydrophobic cores of proteins affect their stabilities. Here, we estimate these effects computationally, using four-body likelihood potentials obtained by simplicial neighborhood analysis of protein packing (SNAPP). In this procedure, the volume of a known protein structure is tiled with tetrahedra having the center of mass of one amino acid side-chain at each vertex. Log-likelihoods are computed for the 8855 possible tetrahedra with equivalent compositions from structural databases and amino acid frequencies. The sum of these four-body potentials for tetrahedra present in a given protein yields the SNAPP score. Mutations change this sum by changing the compositions of tetrahedra containing the mutated residue and their related potentials. Linear correlation coefficients between experimental mutational stability changes, Delta(DeltaG(unfold)), and those based on SNAPP scoring range from 0.70 to 0.94 for hydrophobic core mutations in five different proteins. Accurate predictions for the effects of hydrophobic core mutations can therefore be obtained by virtual mutagenesis, based on changes to the total SNAPP likelihood potential. Significantly, slopes of the relation between Delta(DeltaG(unfold)) and DeltaSNAPP for different proteins are statistically distinct, and we show that these protein-specific effects can be estimated using the average SNAPP score per residue, which is readily derived from the analysis itself. This result enhances the predictive value of statistical potentials and supports previous suggestions that "comparable" mutations in different proteins may lead to different Delta(DeltaG(unfold)) values because of differences in their flexibility and/or conformational entropy.

Protein folding intermediates and pathways studied by hydrogen exchange

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

A kinetically significant intermediate in the folding of barnase

A series of studies on the small protein barnase in the 1990s established it as a paradigm for protein folding in which there is a kinetically important intermediate. But, a recent study in PNAS claims that there are no stable intermediates on the folding pathway. I summarize the evidence that proves that the folding kinetics of barnase is inconsistent with the absence of a folding intermediate. I reinterpret the major evidence presented against the intermediate (an inflection in the unfolding limb of a chevron plot) and show that the inflection is precisely what is predicted from the energy diagram for a three-state reaction with a kinetically significant on-pathway intermediate. The inflection is indicative of a change of rate determining step from the formation to breakdown of an intermediate on unfolding. Other evidence presented against the intermediate is, in fact, consistent with a kinetically important intermediate. I show how the complexities in the kinetics provide a means for measuring otherwise unobtainable rate constants and provide a strategy for mapping the structure of the early transition state in folding. Rather than refute multistate kinetics, the presence of the inflection in the unfolding plot constitutes a novel type of evidence for on-pathway folding intermediates.

Evidence for an unfolding and refolding pathway in cytochrome c

It has been suggested that three partially unfolded forms detected in a native state hydrogen exchange study of oxidized cytochrome c may represent sequential intermediates in an unfolding-refolding pathway. To better define the structure of each intermediate a 'stability labeling' method was used in which the stability of a given segment against unfolding was changed. The condition--folded or unfolded--of the stability-labeled segment in each intermediate could then be determined. The structures found are discrete and native-like and are just those necessary for a sequential pathway.

How do small single-domain proteins fold?

Many small, monomeric proteins fold with simple two-state kinetics and show wide variation in folding rates, from microseconds to seconds. Thus, stable intermediates are not a prerequisite for the fast, efficient folding of proteins and may in fact be kinetic traps and slow the folding process. Using recent studies, can we begin to search for trends which may lead to a better understanding of the protein folding process?

Hydrogen exchange and protein folding

Amide hydrogen-deuterium exchange is a sensitive probe of the structure, stability and dynamics of proteins. The significant increase in the number of small, model proteins that have been studied has allowed a better understanding of the structural fluctuations that lead to hydrogen exchange. Recent technical advances enable the methodology to be applied to the study of protein-protein interactions in much larger, more complex systems.