skip to main content

Explore Research Products in the NSF-PAR

  • NSF Public Access
  • Search Results
  • Targeting the protein folding transition state by mutation: Large scale (un)folding rate accelerations without altering native stability

This content will become publicly available on July 1, 2025

Title: Targeting the protein folding transition state by mutation: Large scale (un)folding rate accelerations without altering native stability
Proteins are constantly undergoing folding and unfolding transitions, with rates that determine their homeostasis in vivo and modulate their biological function. The ability to optimize these rates without affecting overall native stability is hence highly desirable for protein engineering and design. The great challenge is, however, that mutations generally affect folding and unfolding rates with inversely complementary fractions of the net free energy change they inflict on the native state. Here we address this challenge by targeting the folding transition state (FTS) of chymotrypsin inhibitor 2 (CI2), a very slow and stable two‐state folding protein with an FTS known to be refractory to change by mutation. We first discovered that the CI2's FTS is energetically taxed by the desolvation of several, highly conserved, charges that form a buried salt bridge network in the native structure. Based on these findings, we designed a CI2 variant that bears just four mutations and aims to selectively stabilize the FTS. This variant has >250‐fold faster rates in both directions and hence identical native stability, demonstrating the success of our FTS‐centric design strategy. With an optimized FTS, CI2 also becomes 250‐fold more sensitive to proteolytic degradation by its natural substrate chymotrypsin, and completely loses its activity as inhibitor. These results indicate that CI2 has been selected through evolution to have a very unstable FTS in order to attain the kinetic stability needed to effectively function as protease inhibitor. Moreover, the CI2 case showcases that protein (un)folding rates can critically pivot around a few key residues‐interactions, which can strongly modify the general effects of known structural factors such as domain size and fold topology. From a practical standpoint, our results suggest that future efforts should perhaps focus on identifying such critical residues‐interactions in proteins as best strategy to significantly improve our ability to predict and engineer protein (un)folding rates.

  more » « less
Award ID(s):
2112710 2112675
NSF-PAR ID:
10523534
Author(s) / Creator(s):
;
Publisher / Repository:
Wiley
Date Published:
Journal Name:
Protein Science
Volume:
33
Issue:
7
ISSN:
0961-8368
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; ( , PLOS Computational Biology)
    Kasson, Peter M. (Ed.)
    Atomistic simulations can provide valuable, experimentally-verifiable insights into protein folding mechanisms, but existing ab initio simulation methods are restricted to only the smallest proteins due to severe computational speed limits. The folding of larger proteins has been studied using native-centric potential functions, but such models omit the potentially crucial role of non-native interactions. Here, we present an algorithm, entitled DBFOLD, which can predict folding pathways for a wide range of proteins while accounting for the effects of non-native contacts. In addition, DBFOLD can predict the relative rates of different transitions within a protein’s folding pathway. To accomplish this, rather than directly simulating folding, our method combines equilibrium Monte-Carlo simulations, which deploy enhanced sampling, with unfolding simulations at high temperatures. We show that under certain conditions, trajectories from these two types of simulations can be jointly analyzed to compute unknown folding rates from detailed balance. This requires inferring free energies from the equilibrium simulations, and extrapolating transition rates from the unfolding simulations to lower, physiologically-reasonable temperatures at which the native state is marginally stable. As a proof of principle, we show that our method can accurately predict folding pathways and Monte-Carlo rates for the well-characterized Streptococcal protein G. We then show that our method significantly reduces the amount of computation time required to compute the folding pathways of large, misfolding-prone proteins that lie beyond the reach of existing direct simulation. Our algorithm, which is available online , can generate detailed atomistic models of protein folding mechanisms while shedding light on the role of non-native intermediates which may crucially affect organismal fitness and are frequently implicated in disease. 
    more » « less
  2. ; ; ; ; ; ; ; ( , Proceedings of the National Academy of Sciences)

    Protein (un)folding rates depend on the free-energy barrier separating the native and unfolded states and a prefactor term, which sets the timescale for crossing such barrier or folding speed limit. Because extricating these two factors is usually unfeasible, it has been common to assume a constant prefactor and assign all rate variability to the barrier. However, theory and simulations postulate a protein-specific prefactor that contains key mechanistic information. Here, we exploit the special properties of fast-folding proteins to experimentally resolve the folding rate prefactor and investigate how much it varies among structural homologs. We measure the ultrafast (un)folding kinetics of five natural WW domains using nanosecond laser-induced temperature jumps. All five WW domains fold in microseconds, but with a 10-fold difference between fastest and slowest. Interestingly, they all produce biphasic kinetics in which the slower phase corresponds to reequilibration over the small barrier (<3RT) and the faster phase to the downhill relaxation of the minor population residing at the barrier top [transition state ensemble (TSE)]. The fast rate recapitulates the 10-fold range, demonstrating that the folding speed limit of even the simplest all-β fold strongly depends on the amino acid sequence. Given this fold’s simplicity, the most plausible source for such prefactor differences is the presence of nonnative interactions that stabilize the TSE but need to break up before folding resumes. Our results confirm long-standing theoretical predictions and bring into focus the rate prefactor as an essential element for understanding the mechanisms of folding.

     
    more » « less
  3. ( , ChemRxiv)
    Hydrostatic pressure together with the temperature is an important environmental variable that plays an essential role in biological adaptation of extremophilic organisms. In particular, the effects of hy-drostatic pressure on the rates of the protein folding/unfolding reaction are determined by the magni-tude and sign of the activation volume changes. Here we provide computational description of the ac-tivation volume changes for folding/unfolding reaction, and compare them with the experimental data for six different globular proteins. We find that the volume of the transition state ensemble is always in-between the folded and unfolded states. Based on this, we conclude that hydrostatic pressure will invariably slow down protein folding and accelerate protein unfolding. 
    more » « less
  4. ; ; ; ( , Proceedings of the National Academy of Sciences)

    Many large proteins suffer from slow or inefficient folding in vitro. It has long been known that this problem can be alleviated in vivo if proteins start folding cotranslationally. However, the molecular mechanisms underlying this improvement have not been well established. To address this question, we use an all-atom simulation-based algorithm to compute the folding properties of various large protein domains as a function of nascent chain length. We find that for certain proteins, there exists a narrow window of lengths that confers both thermodynamic stability and fast folding kinetics. Beyond these lengths, folding is drastically slowed by nonnative interactions involving C-terminal residues. Thus, cotranslational folding is predicted to be beneficial because it allows proteins to take advantage of this optimal window of lengths and thus avoid kinetic traps. Interestingly, many of these proteins’ sequences contain conserved rare codons that may slow down synthesis at this optimal window, suggesting that synthesis rates may be evolutionarily tuned to optimize folding. Using kinetic modeling, we show that under certain conditions, such a slowdown indeed improves cotranslational folding efficiency by giving these nascent chains more time to fold. In contrast, other proteins are predicted not to benefit from cotranslational folding due to a lack of significant nonnative interactions, and indeed these proteins’ sequences lack conserved C-terminal rare codons. Together, these results shed light on the factors that promote proper protein folding in the cell and how biomolecular self-assembly may be optimized evolutionarily.

     
    more » « less
  5. ; ( , Proceedings of the National Academy of Sciences)

    Single amino acid mutations provide quantitative insight into the energetics that underlie the dynamics and folding of membrane proteins. Chemical denaturation is the most widely used assay and yields the change in unfolding free energy (ΔΔG). It has been applied to >80 different residues of bacteriorhodopsin (bR), a model membrane protein. However, such experiments have several key limitations: 1) a nonnative lipid environment, 2) a denatured state with significant secondary structure, 3) error introduced by extrapolation to zero denaturant, and 4) the requirement of globally reversible refolding. We overcame these limitations by reversibly unfolding local regions of an individual protein with mechanical force using an atomic-force-microscope assay optimized for 2 μs time resolution and 1 pN force stability. In this assay, bR was unfolded from its native bilayer into a well-defined, stretched state. To measure ΔΔG, we introduced two alanine point mutations into an 8-amino-acid region at the C-terminal end of bR’s G helix. For each, we reversibly unfolded and refolded this region hundreds of times while the rest of the protein remained folded. Our single-molecule–derived ΔΔGfor mutant L223A (−2.3 ± 0.6 kcal/mol) quantitatively agreed with past chemical denaturation results while our ΔΔGfor mutant V217A was 2.2-fold larger (−2.4 ± 0.6 kcal/mol). We attribute the latter result, in part, to contact between Val217and a natively bound squalene lipid, highlighting the contribution of membrane protein–lipid contacts not present in chemical denaturation assays. More generally, we established a platform for determining ΔΔGfor a fully folded membrane protein embedded in its native bilayer.

     
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email