Folding Transition Research Articles

Variation of CI-2 Conformers upon Addition of Methanol to Water: An IMS-MS-Based Thermodynamic Analysis.

Chymotrypsin inhibitor 2 (CI-2) is a well-studied, textbook example of a cooperative, two-state, native ↔ denatured folding transition. A recent hybrid ion mobility spectrometry (IMS)/mass spectrometry (MS) thermal denaturation study of CI-2 (the well-studied truncated 64-residue model) in water reported evidence that this two-state transition involves numerous (∼41) unique native and non-native (denatured) solution conformations. The characterization of so many, often low-abundance, states is possible because of the very high dynamic range of IMS-MS measurements of ionic species that are produced upon electrospraying CI-2 solutions from a variable temperature electrospray ionization source. A thermodynamic analysis of these states revealed large changes in enthalpy (ΔH) and entropy (ΔS) at different temperatures, and it was suggested that such variation might arise because of temperature-dependent conformational changes of the protein in response to changes in the conformational entropy and the dielectric permeability of water, which drops from a value of ε ∼ 79 at 24 °C to ∼ 60 at 82 °C. Herein, we examine how adding methanol to water influences the distributions of CI-2 conformers and their ensuing stabilities. The dielectric constant of a 60:40 water:methanol (MeOH) drops from ε ∼ 60 at 24 °C to ∼ 51 at 64 °C. Although the same set of conformers observed in water appears to be present in 60:40 water:MeOH, the abundance of each is substantially altered by the presence of methanol. Relative free energy values (ΔG) and thermodynamic values [ΔH and ΔS and heat capacities (ΔCp)] are derived from a Gibbs-Helmholtz analysis. A comparison of these data from water and water:MeOH systems allows rare insight into how variations in solvation and temperature affect many-state protein equilibria. While these studies confirm that variations in solvent dielectric constant with temperature affect the distributions of conformers that are observed, our findings suggest that other solvent differences may also affect abundances.

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.

Cholesterol catalyzes unfolding in membrane-inserted motifs of the pore forming protein cytolysin A

Plasma membrane-induced protein folding and conformational transitions play a central role in cellular homeostasis. Several transmembrane proteins are folded in the complex lipid milieu to acquire a specific structure and function. Bacterial pore forming toxins (PFTs) are proteins expressed by a large class of pathogenic bacteria that exploit the plasma membrane environment to efficiently undergo secondary structure changes, oligomerize, and form transmembrane pores. Unregulated pore formation causes ion imbalance, leading to cell death and infection. Determining the free energy landscape of these membrane-driven-driven transitions remains a challenging problem. Although cholesterol recognition is required for lytic activity of several proteins in the PFT family of toxins, the regulatory role of cholesterol for the α-PFT, cytolysin A expressed by Escherichia coli remains unexplained. In a recent free energy computation, we showed that the β tongue, a critical membrane-inserted motif of the ClyA toxin, has an on-pathway partially unfolded intermediate that refolds into the helix-turn-helix motif of the pore state. To understand the molecular role played by cholesterol, we carry out string-method-based computations in membranes devoid of cholesterol, which reveals an increase of ∼30 times in the free energy barrier for the loss of β sheet secondary structure when compared with membranes containing cholesterol. Specifically, the tyrosine-cholesterol interaction was found to be critical to creating the unfolded intermediate. Cholesterol also increases the packing and hydrophobicity of the bilayer, resulting in enhanced interactions of the bound protein before complete membrane insertion. Our study illustrates that cholesterol is critical to catalyzing and stabilizing the membrane-inserted unfolded state of the β tongue motif of ClyA, opening up fresh insights into cholesterol-assisted unfolding of membrane proteins.

Detecting Molecular Folding from Noise Measurements

Detecting conformational transitions in molecular systems is key to understanding biological processes. Here, we investigate the force variance in single-molecule pulling experiments as an indicator of molecular folding transitions. We consider cases where Brownian force fluctuations are large, masking the force rips and jumps characteristics of conformational transitions. We compare unfolding and folding data for DNA hairpin systems of loop sizes 4, 8, and 20 and the 110-amino acid protein barnase, finding conditions that facilitate the detection of folding events at low forces where the signal-to-noise ratio is low. In particular, we discuss the role of temperature as a useful parameter to improve the detection of folding transitions in entropically driven processes where folding forces are temperature independent. The force variance approach might be extended to detect the elusive intermediate states in RNA and protein folding.

Features of Protein Unfolding Transitions and Their Relation to Domain Topology Probed by Single-Molecule FRET.

A protein fold is defined as a structural arrangement of a secondary structure in a three-dimensional space. It would be interesting to know whether a particular fold can be assigned to certain features of the corresponding folding/unfolding transitions. To understand the underlying principles of the manifold folding transitions in more detail, single-molecule FRET is the method of choice. Taking the two-domain protein phosphoglycerate kinase (PGK) as an example, we investigated denaturant-induced unfolded states of PGK using the above method. For this purpose, different intramolecular distances within the two domains were measured. In addition to the known two-state transition, a transition with a compact folding intermediate was also identified in each of the two domains. Based on the structural homology of the domains (characterized by a Rossmann fold) and the striking similarity in the features of the measured distance changes during unfolding, clear evidence emerged that the underlying domain topology plays an important role in determining the observed structural changes.

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.

Effects of sequence-dependent non-native interactions in equilibrium and kinetic folding properties of knotted proteins.

Determining the role of non-native interactions in folding dynamics, kinetics, and mechanisms is a classic problem in protein folding. More recently, this question has witnessed a renewed interest in light of the hypothesis that knotted proteins require the assistance of non-native interactions to fold efficiently. Here, we conduct extensive equilibrium and kinetic Monte Carlo simulations of a simple off-lattice C-alpha model to explore the role of non-native interactions in the thermodynamics and kinetics of three proteins embedding a trefoil knot in their native structure. We find that equilibrium knotted conformations are stabilized by non-native interactions that are non-local, and proximal to native ones, thus enhancing them. Additionally, non-native interactions increase the knotting frequency at high temperatures, and in partially folded conformations below the transition temperatures. Although non-native interactions clearly enhance the efficiency of transition from an unfolded conformation to a partially folded knotted one, they are not required to efficiently fold a knotted protein. Indeed, a native-centric interaction potential drives the most efficient folding transition, provided that the simulation temperature is well below the transition temperature of the considered model system.

Compression of a confined semiflexible polymer under direct and oscillating fields.

The folding transition of biopolymers from the coil to compact structures has attracted wide research interest in the past and is well studied in polymer physics. Recent seminal works on DNA in confined devices have shown that these long biopolymers tend to collapse under an external field, which is contrary to the previously reported stretching of the chain. In this work, we capture the compression of a confined semiflexible polymer under direct and oscillating fields using a coarse-grained computer simulation model in the presence of long-range hydrodynamics. In the case of a semiflexible polymer chain, the inhomogeneous hydrodynamic drag from the center to the periphery of the coil couples with the chain bending to cause a swirling movement of the chain segments, leading to structural intertwining and compaction. Contrarily, a flexible chain of the same length lacks such structural deformation and forms a well-established tadpole structure. While bending rigidity profoundly influences the chain's folding favorability, we also found that subject to the direct field, chains in stronger confinements exhibit substantial compaction, contrary to the one in moderate confinements or bulk where such compaction is absent. However, an alternating field within an optimum frequency can effectuate this compression even in moderate or no confinement. This field-induced collapse is a quintessential hydrodynamic phenomenon, resulting in intertwined knotted structures even for shorter chains, unlike other spontaneous knotting experiments where it happens exclusively for longer chains.

General trends in the effects of VX-661 and VX-445 on the plasma membrane expression of clinical CFTR variants

Cystic fibrosis (CF) is caused by mutations that compromise the expression and/or function of the cystic fibrosis transmembrane conductance regulator (CFTR) chloride channel. Most people with CF harbor a common misfolded variant (ΔF508) that can be partially rescued by therapeutic "correctors" that restore its expression. Nevertheless, many other CF variants are insensitive to correctors. Using deep mutational scanning, we quantitatively compare the effects of two correctors on the plasma membrane expression of 129 CF variants. Though structural calculations suggest corrector binding provides similar stabilization to most variants, it's those with intermediate expression and mutations near corrector binding pockets that exhibit the greatest response. Deviations in sensitivity appear to depend on the degree of variant destabilization and the timing of misassembly. Combining correctors appears to rescue more variants by doubling the binding energy and stabilizing distinct cotranslational folding transitions. These results provide an overview of rare CF variant expression and establish new tools for precision pharmacology.

Loss of Residues 119–136, Including the First β-strand of Human Prion Protein, Generates an Aggregation-competent Partially “Open” Form

In prion replication, the cellular form of prion protein (PrPC) must undergo a full conformational transition to its disease-associated fibrillar form. Transmembrane forms of PrP have been implicated in this structural conversion. The cooperative unfolding of a structural core in PrPC presents a substantial energy barrier to prion formation, with membrane insertion and detachment of parts of PrP presenting a plausible route to its reduction. Here, we examined the removal of residues 119–136 of PrP, a region which includes the first β-strand and a substantial portion of the conserved hydrophobic region of PrP, a region which associates with the ER membrane, on the structure, stability and self-association of the folded domain of PrPC. We see an “open” native-like conformer with increased solvent exposure which fibrilises more readily than the native state. These data suggest a stepwise folding transition, which is initiated by the conformational switch to this “open” form of PrPC.

Effect of environment and charge patterning on the dynamics of marginally disordered proteins.

Intrinsically disordered proteins (IDPs) and proteins with long intrinsically disordered regions (IDRs) compose a large fraction of the proteome in all domains of life. Despite an extensive body of literature on “ordinary” disordered proteins—ones with high charge to hydropathy ratio—disordered proteins with lower charge to hydropathy ratio remain underinvestigated. To what extent do such “borderline” proteins compactify or perhaps even fold under physiological conditions? In addition, how does the patterning of charged residues in the primary structure of these IDPs affect the three-dimensional structure? Here we aim to characterize three proteins with various degrees of disorder in vitro and in cell using a suite of biophysical techniques. Our results show that the chain dimensions of all three proteins are sensitive to temperature, ionic strength, and crowder concentration in vitro and inside cells, as predicted from analytical and computational models. Moreover, we observe cooperative collapse for two of the three proteins, reminiscent of folding transitions of globular proteins. We further show that charge decoration is the primary predictor of IDP chain dimensions and models that rely solely on the polymer net charge fail to capture the relative chain dimensions of the mutants of the IDPs in our study.

Flow-Matching: Efficient Coarse-Graining of Molecular Dynamics without Forces

Coarse-grained (CG) molecular simulations have become a standard tool to study molecular processes on time and length scales inaccessible to all-atom simulations. Parametrizing CG force fields to match all-atom simulations has mainly relied on force-matching or relative entropy minimization, which require many samples from costly simulations with all-atom or CG resolutions, respectively. Here we present flow-matching, a new training method for CG force fields that combines the advantages of both methods by leveraging normalizing flows, a generative deep learning method. Flow-matching first trains a normalizing flow to represent the CG probability density, which is equivalent to minimizing the relative entropy without requiring iterative CG simulations. Subsequently, the flow generates samples and forces according to the learned distribution in order to train the desired CG free energy model via force-matching. Even without requiring forces from the all-atom simulations, flow-matching outperforms classical force-matching by an order of magnitude in terms of data efficiency and produces CG models that can capture the folding and unfolding transitions of small proteins.

Nanopore single-molecule biosensor in protein denaturation analysis

Unclear issues in protein studies include but not limited to the stability and denaturation mechanism in the presence of denaturants. Herein, we report a dynamic monitoring approach based on nanopore single-molecule biosensor, which can detect the protein's folding and unfolding transitions by recording a nanopore ionic current. When gradually increasing the concentration of denaturant guanidine hydrochloride (GdmCl), sensitive responses were observed with lysozyme unfolding. The emergence of the featured biphasic-pulse demonstrated the existence of a stable intermediate. It was the first time to experimentally confirm the dynamic equilibrium between the intermediate and the native states at single molecule level, therefore consolidating the standpoint of lysozyme denaturation process following the three-state model. Additionally, we got more insights into the conformation about the intermediate as globular-like structure, larger gyration radius, and enhanced positive charge density. We considered that the manner of denaturant toward lysozyme adopts the “direct” model based on stronger electrostatic and van der Waals forces. Nanopore biosensor exhibited excellent sensitivity with a low detection concentration of 280 pM and reproducibility in analysing the folding intermediate of lysozyme.

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.

Topological origin of the protein folding transition.

In this paper, a geometrical and thermodynamical analysis of the global properties of the potential energy landscape of a minimalistic model of a polypeptide is presented. The global geometry of the potential energy landscape is supposed to contain relevant information about the properties of a given sequence of amino acids, that is, to discriminate between a random heteropolymer and a protein. By considering the SH3 and PYP protein-sequences and their randomized versions it turns out that, in addition to the standard signatures of the folding transition-discriminating between protein sequences of amino acids and random heteropolymer sequences-also peculiar geometric signatures of the equipotential hypersurfaces in configuration space can discriminate between proteins and random heteropolymers. Interestingly, these geometric signatures are the "shadows" of deeper topological changes that take place in correspondence with the protein folding transition. The protein folding transition takes place in systems with a small number of degrees of freedom (very far from the Avogadro number) and in the absence of a symmetry-breaking phenomenon. Nevertheless, seen from the deepest level of topology changes of equipotential submanifolds of phase space, the protein folding transition fully qualifies as a phase transition.

Analytical Approaches for Deriving Friction Coefficients for Selected α-Helical Peptides Based Entirely on Molecular Dynamics Simulations.

In this paper we derive analytically from molecular dynamics (MD) simulations the friction coefficients related to conformational transitions within several model peptides with α-helical structures. We study a series of alanine peptides with various length from ALA5 to ALA21 as well as their two derivatives, the (AAQAA)3 peptide and a 13-residue KR1 peptide that is a derivative of the (AAQAA)2 peptide with the formula GN(AAQAA)2G. We use two kinds of approaches to derive their friction coefficients. In the local approach, friction associated with fluctuations of single hydrogen bonds are studied. In the second approach, friction coefficients associated with a folding transitions within the studied peptides are obtained. In both cases, the respective friction coefficients differentiated very well the subtle structural changes between studied peptides and compared favorably to experimentally available data.

Structural and thermodynamic analyses of the β-to-α transformation in RfaH reveal principles of fold-switching proteins.

The two-domain protein RfaH, a paralog of the universally conserved NusG/Spt5 transcription factors, is regulated by autoinhibition coupled to the reversible conformational switch of its 60-residue C-terminal Kyrpides, Ouzounis, Woese (KOW) domain between an α-hairpin and a β-barrel. In contrast, NusG/Spt5-KOW domains only occur in the β-barrel state. To understand the principles underlying the drastic fold switch in RfaH, we elucidated the thermodynamic stability and the structural dynamics of two RfaH- and four NusG/Spt5-KOW domains by combining biophysical and structural biology methods. We find that the RfaH-KOW β-barrel is thermodynamically less stable than that of most NusG/Spt5-KOWs and we show that it is in equilibrium with a globally unfolded species, which, strikingly, contains two helical regions that prime the transition toward the α-hairpin. Our results suggest that transiently structured elements in the unfolded conformation might drive the global folding transition in metamorphic proteins in general.

The Effect of Time Resolution on Apparent Transition Path Times Observed in Single-Molecule Studies of Biomolecules.

Single-molecule experiments have now achieved a time resolution allowing observation of transition paths, the brief trajectory segments where the molecule undergoing an unfolding or folding transition enters the energetically or entropically unfavorable barrier region from the folded/unfolded side and exits to the unfolded/folded side, thereby completing the transition. This resolution, however, is yet insufficient to identify the precise entrance/exit events that mark the beginning and the end of a transition path: the nature of the diffusive dynamics is such that a molecular trajectory will recross the boundary between the barrier region and the folded/unfolded state, multiple times, at a time scale much shorter than that of the typical experimental resolution. Here we use theory and Brownian dynamics simulations to show that, as a result of such recrossings, the apparent transition path times are generally longer than the true ones. We quantify this effect using a simple model where the observed dynamics is a moving average of the true dynamics and discuss experimental implications of our results.

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.