Folding Transition State Research Articles

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.

Hydrogen bonding heterogeneity correlates with protein folding transition state passage time as revealed by data sonification

Protein-protein and protein-water hydrogen bonding interactions play essential roles in the way a protein passes through the transition state during folding or unfolding, but the large number of these interactions in molecular dynamics (MD) simulations makes them difficult to analyze. Here, we introduce a state space representation and associated "rarity" measure to identify and quantify transition state passage (transit) events. Applying this representation to a long MD simulation trajectory that captured multiple folding and unfolding events of the GTT WW domain, a small protein often used as a model for the folding process, we identified three transition categories: Highway (faster), Meander (slower), and Ambiguous (intermediate). We developed data sonification and visualization tools to analyze hydrogen bond dynamics before, during, and after these transition events. By means of these tools, we were able to identify characteristic hydrogen bonding patterns associated with "Highway" versus "Meander" versus "Ambiguous" transitions and to design algorithms that can identify these same folding pathways and critical protein-water interactions directly from the data. Highly cooperative hydrogen bonding can either slow down or speed up transit. Furthermore, an analysis of protein-water hydrogen bond dynamics at the surface of WW domain shows an increase in hydrogen bond lifetime from folded to unfolded conformations with Ambiguous transitions as an outlier. In summary, hydrogen bond dynamics provide a direct window into the heterogeneity of transits, which can vary widely in duration (by a factor of 10) due to a complex energy landscape.

The Mechanism of Folding of Human Frataxin in Comparison to the Yeast Homologue – Broad Energy Barriers and the General Properties of the Transition State

The funneled energy landscape theory suggests that the folding pathway of homologous proteins should converge at the late stages of folding. In this respect, proteins displaying a broad energy landscape for folding are particularly instructive, allowing inferring both the early, intermediate and late stages of folding. In this paper we explore the folding mechanisms of human frataxin, an essential mitochondrial protein linked to the neurodegenerative disorder Friedreich's ataxia. Building upon previous studies on the yeast homologue, the folding pathway of human frataxin is thoroughly examined, revealing a mechanism implying the presence of a broad energy barrier, reminiscent of the yeast counterpart. Through an extensive site-directed mutagenesis, we employed a Φ -value analysis to map native-like contacts in the folding transition state. The presence of a broad energy barrier facilitated the exploration of such contacts in both early and late folding events. We compared results from yeast and human frataxin providing insights into the impact of native topology on the folding mechanism and elucidating the properties of the underlying free energy landscape. The findings are discussed in the context of the funneled energy landscape theory of protein folding.

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.

Topological Reaction Coordinate Captures the Folding Transition State Ensemble in a Pierced Lasso Protein.

Proteins with a pierced lasso topology (PLT) have a covalent loop created by a disulfide bond, and the backbone circles back to thread the loop. This threaded topology has unique features compared to knotted topologies; notably, the topology is controlled by the chemical environment and the covalent loop remains intact even when denatured. In this work, we use the hormone leptin as our model PLT system and study its folding using molecular dynamics simulations that employ a structure-based (Go̅-like) model. We find that the reduced protein has a two-state folding mechanism with a transition state ensemble (TSE) that can be characterized by the reaction coordinate Q, the fraction of native contacts formed. In contrast, the oxidized protein, which must thread part of the polypeptide chain through a covalent loop, has a folding process that is poorly characterized by Q. Instead, we find that a topological coordinate that monitors the residue crossing the loop can identify the TSE of oxidized leptin. By precisely identifying the predicted TSE, one may now reliably calculate theoretical phi-values for the PLT protein, thereby enabling a comparison with experimental measurements. We find the loop-threading constraint leads to noncanonical phi-values that are uniformly small because this PLT protein has a flat energy landscape through the TSE.

Protein Folding Stability and Kinetics in Alginate Hydrogels.

Proteins are commonly encapsulated in alginate gels for drug delivery and tissue-engineering applications. However, there is limited knowledge of how encapsulation impacts intrinsic protein properties such as folding stability or unfolding kinetics. Here, we use fast relaxation imaging (FReI) to image protein unfolding in situ in alginate hydrogels after applying a temperature jump. Based on changes in the Förster resonance energy transfer (FRET) response of FRET-labeled phosphoglycerate kinase (PGK), we report the quantitative impact of multiple alginate hydrogel concentrations on protein stability and folding dynamics. The gels stabilize PGK by increasing its melting temperature up to 18.4 °C, and the stabilization follows a nonmonotonic dependence on the alginate density. In situ kinetic measurements also reveal that PGK deviates more from two-state folding behavior in denser gels and that the gel decreases the unfolding rate and accelerates the folding rate of PGK, compared to buffer. Phi-value analysis suggests that the folding transition state of an encapsulated protein is structurally similar to that of folded protein. This work reveals both beneficial and negative impacts of gel encapsulation on protein folding, as well as potential mechanisms contributing to altered stability.

Characterizing the Folding Transition-State Ensembles in the Energy Landscape of an RNA Tetraloop.

Molecular dynamics (MD) simulations have become increasingly powerful and can now describe the folding/unfolding of small biomolecules in atomic detail. However, a major challenge in MD simulations is to represent the complex energy landscape of biomolecules using a small number of reaction coordinates. In this study, we investigate the folding pathways of an RNA tetraloop, gcGCAAgc, using five classical MD simulations with a combined simulation time of approximately 120 μs. Our approach involves analyzing the tetraloop dynamics, including the folding transition state ensembles, using the energy landscape visualization method (ELViM). The ELViM is an approach that uses internal distances to compare any two conformations, allowing for a detailed description of the folding process without requiring root mean square alignment of structures. This method has previously been applied to describe the energy landscape of disordered β-amyloid peptides and other proteins. The ELViM results in a non-linear projection of the multidimensional space, providing a comprehensive representation of the tetraloop's energy landscape. Our results reveal four distinct transition-state regions and establish the paths that lead to the folded tetraloop structure. This detailed analysis of the tetraloop's folding process has important implications for understanding RNA folding, and the ELViM approach can be used to study other biomolecules.

Evolutionary conservation of amino acids contributing to the protein folding transition state.

The question of whether amino acids critical to protein folding kinetics are evolutionarily conserved has been investigated intensively in the past, but no consensus has yet been reached. Recently, we have demonstrated that the transition state, dictating folding kinetics, is characterized as the state of maximum dynamic cooperativity, i.e., the state of maximum correlations between amino acid contact formations. Here, we investigate the evolutionary conservation of those amino acids contributing significantly to the dynamic cooperativity. We find a strong indication of a new kind of relationship-necessary but not sufficient causality-between the evolutionary conservation and the dynamic cooperativity: larger contributions to the dynamic cooperativity arise from more conserved residues, but not vice versa. This holds for all the protein systems for which long folding simulation trajectories are available. To our knowledge, this is the first systematic demonstration of any kind of evolutionary conservation of amino acids relevant to folding kinetics.

Ligand-Dependent Volumetric Characterization of Manganese Riboswitch Folding: A High-Pressure Single-Molecule Kinetic Study.

Nanoscopic differences in free volume result in pressure-dependent changes in free energies which can therefore impact folding/unfolding stability of biomolecules. Although such effects are typically insignificant under ambient pressure conditions, they are crucially important for deep ocean marine life, where the hydraulic pressure can be on the kilobar scale. In this work, single molecule FRET spectroscopy is used to study the effects of pressure on both the kinetics and overall thermodynamics for folding/unfolding of the manganese riboswitch. Detailed pressure-dependent analysis of the conformational kinetics allows one to extract precision changes (σ ≲ 4-8 Å3) in free volumes not only between the fully folded/unfolded conformations but also with respect to the folding transition state of the manganese riboswitch. This permits first extraction of a novel "reversible work" free energy (PΔV) landscape, which reveals a monotonic increase in manganese riboswitch volume along the folding coordinate. Furthermore, such a tool permits exploration of pressure-dependent effects on both Mn2+ binding and riboswitch folding, which demonstrate that ligand attachment stabilizes the riboswitch under pressure by decreasing the volume increase upon folding (ΔΔV < 0). Such competition between ligand binding and pressure-induced denaturation dynamics could be of significant evolutionary advantage, compensating for a weakening in riboswitch tertiary structure with pressure-mediated ligand binding and promotion of folding response.

A short structural extension dictates the early stages of folding of a PDZ domain

PDZ domains are highly abundant protein-protein interaction modules in human. One of the most extensively characterized PDZ domain, the third PDZ domain from PSD-95 (PDZ3), contains an α-helical C-terminal extension that has a key role in the function of the domain. Here we compared the folding of PDZ3 with a truncated variant (PDZ3Δα3), lacking the additional helix, by means of the so-called Φ-value analysis, an experimental technique that allows inferring the structure of folding transition states. Experiments reveal subtle but detectable differences in the folding of PDZ3Δα3 versus PDZ3, as probed by structural characterization of the folding transition states. These differences appear more remarkable in the early stages of folding, with a detectable shift of the folding nucleus. The presented results allow demonstrating that the native state exerts a weak bias at the early stages of folding, which appear to be characterized by alternative pathways.

Mutational Analysis of Protein Folding Transition States: Phi Values.

The analysis of protein folding reactions by monitoring the kinetic effects of specifically designed single-point mutations, the so-termed phi-value analysis, has been a favorite technique to experimentally probe the mechanisms of protein folding. The idea behind phi-value analysis is that the effects that mutations have on the folding and unfolding rate constants report on the energetic/structural features of the folding transition state ensemble (TSE), which is the highest point in the free energy surface connecting the native and unfolded states, and thus the rate limiting step that ultimately defines the folding mechanism. For single-domain, two-state folding proteins, the general procedure to perform the phi-value analysis of protein folding is relatively simple to implement in the lab. Once the mutations have been produced and purified, the researcher needs to follow a few specific guidelines to perform the experiments and to analyze the data so produced. In this chapter, a step-by-step description of how to measure and interpret the effects induced by site-directed mutations on the folding and unfolding rate constants of a protein of interest is provided. Some possible solutions to the most typical problems that arise when performing phi-value analysis in the lab are also provided.

Engineered Metal-Binding Sites to Probe Protein Folding Transition States: Psi Analysis.

The formation of the transition state ensemble (TSE) represents the rate-limiting step in protein folding. The TSE is the least populated state on the pathway, and its characterization remains a challenge. Properties of the TSE can be inferred from the effects on folding and unfolding rates for various perturbations. A difficulty remains on how to translate these kinetic effects to structural properties of the TSE. Several factors can obscure the translation of point mutations in the frequently used method, "mutational Phi analysis." We take a complementary approach in "Psi analysis," employing rationally inserted metal binding sites designed to probe pairwise contacts in the TSE. These contacts can be confidently identified and used to construct structural models of the TSE. The method has been applied to multiple proteins and consistently produces a considerably more structured and native-like TSE than Phi analysis. This difference has significant implications to our understanding of protein folding mechanisms. Here we describe the application of the method and discuss how it can be used to study other conformational transitions such as binding.

Measuring Excess Heat Capacities of Deoxyribonucleic Acid (DNA) Folding at the Single-Molecule Level.

Measurements of the thermodynamic properties of biomolecular folding (ΔG°, ΔH°, ΔS°, etc.) provide a wealth of information on the folding process and have long played a central role in biophysical investigation. In particular, the excess heat capacity of folding (ΔCP) is crucial, as typically measured in bulk ensemble studies by differential scanning calorimetry (DSC) and isothermal titration calorimetry (ITC). Here, we report the first measurements of ΔCP at the single-molecule level using the single-molecule fluorescence resonance energy transfer (smFRET) as well as the very first measurements of the heat capacity change associated with achieving the transition state (ΔC‡P) for nucleic acid folding. The deoxyribonucleic acid (DNA) hairpin used in these studies exhibits an excess heat capacity for hybridization (ΔCP = -340 ± 60 J/mol/K per base pair) consistent with the range of literature expectations (ΔCP = -100 to -420 J/mol/K per base pair). Furthermore, the measured activation heat capacities (ΔC‡P) for such hairpin unfolding are consistent with a folding transition state containing few fully formed base pairs, in agreement with prevailing models of DNA hybridization.

The Zero-Order Loop in Apoazurin Modulates Folding Mechanism In Silico.

Pseudomonas aeruginosa apoazurin (apo, without the copper cofactor) has a single disulfide bond between residues 3 and 26 and unfolds in a two-state reaction in vitro. The disulfide bond covalently connects the N-termini of β-strands 1 and 3; thereby, it creates a zero-order loop or a "cinch" that restricts conformational space. Covalent loops and threaded topologies are emerging as important structural elements in folded proteins and may be important for function. In order to understand the role of a zero-order loop in the folding process of a protein, here we used coarse-grained molecular dynamics (CGMD) simulations in silico to compare two variants of apoazurin: one named "loop" which contained the disulfide, and another named "open" in which the disulfide bond between residues 3 and 26 was removed. CGMD simulations were performed to probe the stability and unfolding pathway of the two apoazurin variants at different urea concentrations and temperatures. Our results show that the covalent loop plays a prominent role in the unfolding mechanism of apoazurin; its removal alters both the folding-transition state and the unfolded-state ensemble of conformations. We propose that modulation of azurin's folding landscape by the disulfide bridge may be related to both copper capturing and redox sensing.

Time-dependent communication between multiple amino acids during protein folding

Cooperativity is considered to be a key organizing principle behind biomolecular assembly, recognition and folding. However, it has remained very challenging to quantitatively characterize how cooperative processes occur on a concerted, multiple-interaction basis. Here, we address how and when the folding process is cooperative on a molecular scale. To this end, we analyze multipoint time-correlation functions probing time-dependent communication between multiple amino acids, which were computed from long folding simulation trajectories. We find that the simultaneous multiple amino-acid contact formation, which is absent in the unfolded state, starts to develop only upon entering the folding transition path. Interestingly, the transition state, whose presence is connected to the macrostate cooperative behavior known as the two-state folding, can be identified as the state in which the amino-acid cooperativity is maximal. Thus, our work not only provides a new mechanistic view on how protein folding proceeds on a multiple-interaction basis, but also offers a conceptually novel characterization of the folding transition state and the molecular origin of the phenomenological cooperative folding behavior. Moreover, the multipoint correlation function approach adopted here is general and can be used to expand the understanding of cooperative processes in complex chemical and biomolecular systems.

Friction-Limited Folding of Disulfide-Reduced Monomeric SOD1

The folding reaction of a stable monomeric variant of Cu/Zn superoxide dismutase (mSOD1), an enzyme responsible for the conversion of superoxide free radicals into hydrogen peroxide and oxygen, is known to be among the slowest folding processes that adhere to two-state behavior. The long lifetime, ∼10 s, of the unfolded state presents ample opportunities for the polypeptide chain to transiently sample nonnative structures before the formation of the productive folding transition state. We recently observed the formation of a nonnative structure in a peptide model of the C-terminus of SOD1, a sequence that might serve as a potential source of internal chain friction-limited folding. To test for friction-limited folding, we performed a comprehensive thermodynamic and kinetic analysis of the folding mechanism of mSOD1 in the presence of the viscogens glycerol and glucose. Using a, to our knowledge, novel analysis of the folding reactions, we found the disulfide-reduced form of the protein that exposes the C-terminal sequence, but not its disulfide-oxidized counterpart that protects it, experiences internal chain friction during folding. The sensitivity of the internal friction to the disulfide bond status suggests that one or both of the cross-linked regions play a critical role in driving the friction-limited folding. We speculate that the molecular mechanisms giving rise to the internal friction of disulfide-reduced mSOD1 might play a role in the amyotrophic lateral sclerosis-linked aggregation of SOD1.

Role of Zero-Order Loop in Protein Unfolding Case Study with Apoazurin

Pseudomonas aeruginosa apoazurin (apo, without the copper cofactor) is a two-state folder that has a single disulfide bond between residues 3 and 26. This bond covalently connects the N-termini of beta-strands 1 and 3; thereby it creates a zero-order loop. The loop restricts the conformational space for the apoazurin polypeptide. In order to understand the role played by the zero-order loop, we used molecular dynamics (MD) simulations to compare two variants of apoazurin; one variant called “loop” which contained the disulfide and another called “open” in which the disulfide bond between residues 3 and 26 was removed. MD simulations were performed to probe the unfolding pathway and stability of the two apoazurin variants at different urea concentrations and temperatures. Our results show that the folded structure apoazurin is somewhat more stable due to the presence of the disulfide bond. However, the disulfide bond plays a prominent role in the apoazurin unfolding mechanism: we find that it changes both the folding-transition state and the unfolded-state ensemble of conformations.

Folding of the Ig-Like Domain of the Dengue Virus Envelope Protein Analyzed by High-Hydrostatic-Pressure NMR at a Residue-Level Resolution.

Dengue fever is a mosquito-borne endemic disease in tropical and subtropical regions, causing a significant public health problem in Southeast Asia. Domain III (ED3) of the viral envelope protein contains the two dominant putative epitopes and part of the heparin sulfate receptor binding region that drives the dengue virus (DENV)’s fusion with the host cell. Here, we used high-hydrostatic-pressure nuclear magnetic resonance (HHP-NMR) to obtain residue-specific information on the folding process of domain III from serotype 4 dengue virus (DEN4-ED3), which adopts the classical three-dimensional (3D) ß-sandwich structure known as the Ig-like fold. Interestingly, the folding pathway of DEN4-ED3 shares similarities with that of the Titin I27 module, which also adopts an Ig-like fold, but is functionally unrelated to ED3. For both proteins, the unfolding process starts by the disruption of the N- and C-terminal strands on one edge of the ß-sandwich, yielding a folding intermediate stable over a substantial pressure range (from 600 to 1000 bar). In contrast to this similarity, pressure-jump kinetics indicated that the folding transition state is considerably more hydrated in DEN4-ED3 than in Titin I27.

Sampling the Folding Transition State Ensemble in a Tube-like Model of Protein

We used the tube model with Go-like potential for native contacts to study the folding transition of a designed three-helix bundle and a designed protein G-like structure. It is shown that both proteins in this model are two-state folders with a cooperative folding transition coincided with the collapse transition. We defined the transition states as protein conformations in a small region around the saddle point on a free energy surface with the energy and the conformationalroot mean square deviation (rmsd) from the native state as the coordinates. The transition state region on the free energy surface then was sampled by using umbrella sampling technique. We show that the transition state ensemble is broad consisting of different conformations that have different folded and unfolded elements.

Helix-Coil Transition Courses Through Multiple Pathways and Intermediates: Fast Kinetic Measurements and Dimensionality Reduction.

Nanosecond laser temperature jumps with tryptophan fluorescence detection and molecular dynamics simulation with kinetic dimensionality reduction were used to study the helix-coil transition in a 21-residue α-helical heteropeptide. Analysis of the temperature- dependent relaxation dynamics of this heteropeptide identified a distinct faster component of 20-35 ns, besides a slower component of 300-400 ns at temperatures between 296 and 280 K. To understand the mechanism of progression from a non-structured coil state to a structured helical state, we carried out a 12 μs molecular dynamics simulation of this peptide system. Clustering and optimal dimensionality reduction were applied to the molecular dynamics trajectory to generate low-dimensional coarse-grained models of the underlying kinetic network in terms of 2-5 metastable states. In accord with the generally accepted understanding of the multiple conformations and high entropy of the unfolded ensemble of states, the "coil" metastable set contains the largest number of structures. Interestingly, the helix metastable state was also found to be structurally heterogeneous, consisting of the completely helical form and several partly folded conformers that interconvert at a time scale faster that global folding. The intermediate states contain the fewest structures, have lowest populations, and have the shortest lifetimes. As the number of considered metastable states increases, more intermediates and more folding paths appear in the coarse-grained models. One of these intermediates corresponds to the transition state for folding, which involves an "off-center" helical region over residues 11-16. The kinetic network model is consistent with a statistical picture of folding following a simple reaction coordinate counting the helical population of individual residues. On the basis of simulations, we propose that the fast relaxation time should be assigned to cooperative folding/unfolding of segments of 1-4 neighboring residues.