Two-state Folding Research Articles

The characteristics of molten globule states and folding pathways strongly depend on the sequence of a protein

ABSTRACTThe majority of proteins perform their cellular function after folding into a specific and stable native structure. Additionally, for many proteins less compact ‘molten globule’ states have been observed. Current experimental observations show that the molten globule state can show varying degrees of compactness and solvent accessibility; the underlying molecular cause for this variation is not well understood. While the specificity of protein folding can be studied using protein lattice models, current design procedures for these models tend to generate sequences without molten globule-like behaviour. Here we alter the design process so the distance between the molten globule ensemble and the native structure can be steered; this allows us to design protein sequences with a wide range of folding pathways, and sequences with well-defined heat-induced molten globules. Simulating these sequences we find that (1) molten globule states are compact, but have less specific configurations compared to the folded state, (2) the nature of the molten globule state is highly sequence dependent, (3) both two-state and multi-state folding proteins may show heat-induced molten globule states, as observed in heat capacity curves. The varying nature of the molten globules and typical heat capacity curves associated with the transitions closely resemble experimental observations.

Monitoring 15N Chemical Shifts During Protein Folding by Pressure-Jump NMR.

Pressure-jump hardware permits direct observation of protein NMR spectra during a cyclically repeated protein folding process. For a two-state folding protein, the change in resonance frequency will occur nearly instantaneously when the protein clears the transition state barrier, resulting in a monoexponential change of the ensemble-averaged chemical shift. However, protein folding pathways can be more complex and contain metastable intermediates. With a pseudo-3D NMR experiment that utilizes stroboscopic observation, we measure the ensemble-averaged chemical shifts, including those of exchange-broadened intermediates, during the folding process. Such measurements for a pressure-sensitized mutant of ubiquitin show an on-pathway kinetic intermediate whose 15N chemical shifts differ most from the natively folded protein for strands β5, its preceding turn, and the two strands that pair with β5 in the native structure.

Prediction of protein folding kinetics states using hybrid brainstorm optimization

ABSTRACT Protein folding procedure is tremendously vital in determining the molecular role. The protein folding kinetics states select if the molecule influences the natural structure of intermediates or not. The fold can be complete with stable intermediates (3-State/3S) or without stable intermediates (2-State/2S). Protein folding is regularly determined, using experiments and often takes time and dreary. Furthermore, there are significant numbers of unfolding mechanisms that are available in the PDB to originate unidentified. Henceforth, it made curiosity and we focused on classifying and envisaging the folding mechanism as three-state or two-state. The authors established the clustering models by modified brainstorm optimization (MBSO) using the Fuzzy C-Means (FCM) algorithm through known parameters (hydrophobicity, length, hydrophilicity, and secondary structural components) to predict the protein folding kinetic states. The prototypes implemented well enough for envisaging three-state and two-state folding using the recognized features. MBSO using FCM shows a better model when compared with MBSO. The MBSO using FCM algorithm produced good accuracy in Davies-Bouldin index and Dunn index for clustering. The result depicts the better performance of MBSO using FCM techniques with the best prediction when compared with MBSO. In future, we can further extend to predict protein fold recognition.

Transient intermediates are populated in the folding pathways of single-domain two-state folding protein L.

Small single-domain globular proteins, which are believed to be dominantly two-state folders, played an important role in elucidating various aspects of the protein folding mechanism. However, recent single molecule fluorescence resonance energy transfer experiments [H. Y. Aviram et al. J. Chem. Phys. 148, 123303 (2018)] on a single-domain two-state folding protein L showed evidence for the population of an intermediate state and it was suggested that in this state, a β-hairpin present near the C-terminal of the native protein state is unfolded. We performed molecular dynamics simulations using a coarse-grained self-organized-polymer model with side chains to study the folding pathways of protein L. In agreement with the experiments, an intermediate is populated in the simulation folding pathways where the C-terminal β-hairpin detaches from the rest of the protein structure. The lifetime of this intermediate structure increased with the decrease in temperature. In low temperature conditions, we also observed a second intermediate state, which is globular with a significant fraction of the native-like tertiary contacts satisfying the features of a dry molten globule.

Two-State Folding Energy Determination Based on Transition Points in Nonequilibrium Single-Molecule Experiments.

Many small protein domains or nucleic acid structures undergo two-state unfolding-refolding transitions during mechanical stretching using single-molecule techniques. Here, by applying the Jarzynski equality (JE), we analytically express the folding energy of a molecule as a function of the experimentally measured transition points ξ* obtained with two typical time-varying mechanical constraints: the force constraints F(t) and the position constraints R(t) of a Hookian spring attached to one end of the molecule. Compared to previous applications of JE based on the integration of accurately measured force-extension curves of a tether that typically contains the molecule of interest and handles, our approach just needs to accurately measure a single data point. In the case of the F(t) process, the calculation is handle-independent. The broad applications of the theory are demonstrated by measuring the folding energies of a DNA hairpin, a DNA G-quadruplex, and the titin I27 domain based on transition forces using magnetic tweezers.

Probing the Effect of the Ribosome on the Protein Folding Pathway using Single-Molecule Chemo-Mechanical Folding

In vivo, proteins function in a complex environment where they are subject to many stresses that can modulate the protein's energy landscape. Using a combination of force and chemical denaturant (chemo-mechanical unfolding), we have demonstrated that a simple, two-state folding protein can fold through multiple parallel pathways and that seemingly small changes in force, denaturant concentration and protein sequence can strongly modulate the flux between the pathways. This result suggests that in vivo, the heterogeneous crowded cellular environment could strongly influence the mechanisms of protein folding and unfolding. Another perturbant that could affect the folding pathway in vivo is the ribosome since many proteins fold while they are being translated. To explore how the ribosome affects the folding mechanism, we use chemo-mechanical unfolding to characterize folding pathways both on and off the ribosome.

Conservation of folding and association within a family of spidroin N-terminal domains

Web spiders synthesize silk fibres, nature’s toughest biomaterial, through the controlled assembly of fibroin proteins, so-called spidroins. The highly conserved spidroin N-terminal domain (NTD) is a pH-driven self-assembly device that connects spidroins to super-molecules in fibres. The degree to which forces of self-assembly is conserved across spider glands and species is currently unknown because quantitative measures are missing. Here, we report the comparative investigation of spidroin NTDs originating from the major ampullate glands of the spider species Euprosthenops australis, Nephila clavipes, Latrodectus hesperus, and Latrodectus geometricus. We characterized equilibrium thermodynamics and kinetics of folding and self-association using dynamic light scattering, stopped-flow fluorescence and circular dichroism spectroscopy in combination with thermal and chemical denaturation experiments. We found cooperative two-state folding on a sub-millisecond time scale through a late transition state of all four domains. Stability was compromised by repulsive electrostatic forces originating from clustering of point charges on the NTD surface required for function. pH-driven dimerization proceeded with characteristic fast kinetics yielding high affinities. Results showed that energetics and kinetics of NTD self-assembly are highly conserved across spider species despite the different silk mechanical properties and web geometries they produce.

Exploring the correlation between the folding rates of proteins and the entanglement of their native states

The folding of a protein towards its native state is a rather complicated process. However, there is empirical evidence that the folding time correlates with the contact order, a simple measure of the spatial organization of the native state of the protein. Contact order is related to the average length of the main chain loops formed by amino acids that are in contact. Here we argue that folding kinetics can also be influenced by the entanglement that loops may undergo within the overall three-dimensional protein structure. In order to explore such a possibility, we introduce a novel descriptor, which we call ‘maximum intrachain contact entanglement’. Specifically, we measure the maximum Gaussian entanglement between any looped portion of a protein and any other non-overlapping subchain of the same protein, which is easily computed by discretized line integrals on the coordinates of the atoms. By analyzing experimental data sets of two-state and multi-state folders, we show that the new index is also a good predictor of the folding rate. Moreover, being only partially correlated with previous methods, it can be integrated with them to yield more accurate predictions.

Two states or not two states: Single-molecule folding studies of protein L.

Experimental tools of increasing sophistication have been employed in recent years to study protein folding and misfolding. Folding is considered a complex process, and one way to address it is by studying small proteins, which seemingly possess a simple energy landscape with essentially only two stable states, either folded or unfolded. The B1-IgG binding domain of protein L (PL) is considered a model two-state folder, based on measurements using a wide range of experimental techniques. We applied single-molecule fluorescence resonance energy transfer (FRET) spectroscopy in conjunction with a hidden Markov model analysis to fully characterize the energy landscape of PL and to extract the kinetic properties of individual molecules of the protein. Surprisingly, our studies revealed the existence of a third state, hidden under the two-state behavior of PL due to its small population, ∼7%. We propose that this minority intermediate involves partial unfolding of the two C-terminal β strands of PL. Our work demonstrates that single-molecule FRET spectroscopy can be a powerful tool for a comprehensive description of the folding dynamics of proteins, capable of detecting and characterizing relatively rare metastable states that are difficult to observe in ensemble studies.

Important model protein found to fold in unexpected way by single-molecule FRET spectroscopy

FRET spectroscopy reveals protein L is not a two-state folder, identifying a third state for the first time.

SCooP: an accurate and fast predictor of protein stability curves as a function of temperature.

The molecular bases of protein stability remain far from elucidated even though substantial progress has been made through both computational and experimental investigations. One of the most challenging goals is the development of accurate prediction tools of the temperature dependence of the standard folding free energy ΔG(T). Such predictors have an enormous series of potential applications, which range from drug design in the biopharmaceutical sector to the optimization of enzyme activity for biofuel production. There is thus an important demand for novel, reliable and fast predictors. We present the SCooP algorithm, which is a significant step towards accurate temperature-dependent stability prediction. This automated tool uses the protein structure and the host organism as sole entries and predicts the full T-dependent stability curve of monomeric proteins assumed to follow a two-state folding transition. Equivalently, it predicts all the thermodynamic quantities associated to the folding transition, namely the melting temperature Tm, the standard folding enthalpy ΔHm measured at Tm, and the standard folding heat capacity ΔCp. The cross-validated performances are good, with correlation coefficients between predicted and experimental values equal to [0.80, 0.83, 0.72] for ΔHm, ΔCp and Tm, respectively, which increase up to [0.88, 0.90, 0.78] upon the removal of 10% outliers. Moreover, the stability curve prediction of a target protein is very fast: it takes less than a minute. SCooP can thus potentially be applied on a structurome scale. This opens new perspectives of large-scale analyses of protein stability, which is of considerable interest for protein engineering. The SCooP webserver is freely available at http://babylone.ulb.ac.be/SCooP. fapucci@ulb.ac.be, jkwasigr@ulb.ac.be or mrooman@ulb.ac.be. Supplementary data are available at Bioinformatics online.

MERS\u2010CoV papain-like protease (PLpro): expression, purification, and spectroscopic/thermodynamic characterization

Within a decade, MERS-CoV emerged with nearly four times higher case fatality rate than an earlier outbreak of SARS-CoV and spread out in 27 countries in short span of time. As an emerging virus, combating it requires an in-depth understanding of its molecular machinery. Therefore, conformational characterization studies of coronavirus proteins are necessary to advance our knowledge of the matter for the development of antiviral therapies. In this study, MERS-CoV papain-like protease (PLpro) was recombinantly expressed and purified. Thermal folding pathway and thermodynamic properties were characterized using dynamic multimode spectroscopy (DMS) and thermal shift assay. DMS study showed that the PLpro undergoes a single thermal transition and follows a pathway of two-state folding with Tm and van’t Hoff enthalpy values of 54.4 ± 0.1 °C and 317.1 ± 3.9 kJ/mol, respectively. An orthogonal technique based on intrinsic tryptophan fluorescence also showed that MERS-CoV PLpro undergoes a single thermal transition and unfolds via a pathway of two-state folding with a Tm value of 51.4 °C. Our findings provide significant understandings of the thermodynamic and structural properties of MERS-CoV PLpro.

Highly Heterogeneous Nature of the Native and Unfolded States of the B Domain of Protein A Revealed by Two-Dimensional Fluorescence Lifetime Correlation Spectroscopy

Elucidating the protein folding mechanism is crucial to understand how proteins acquire their unique structures to realize various biological functions. With this aim, the folding/unfolding of small globular proteins has been extensively studied. Interestingly, recent studies have revealed that even such small proteins represent considerably complex processes. In this study, we examined the folding/unfolding process of a small α-helical protein, the B domain of protein A (BdpA), at equilibrium using two-dimensional fluorescence lifetime correlation spectroscopy with 10 μs time resolution. The results showed that although the BdpA is a two-state folder, both the native and unfolded states are highly heterogeneous and the conformational conversion within each ensemble occurs within 10 μs. Furthermore, it was shown that the average structures of both ensembles gradually change and become more elongated as the denaturant concentration increases. The analysis on two mutants suggested that fraying of the N-terminal helix is the origin of the inhomogeneity of the native state. Because the direct observation of the ensemble nature of the native state at the single-molecule level has not been reported, the data obtained in this study give new insights into complex conformational properties of small proteins.

Simulated Force Quench Dynamics Shows GB1 Protein Is Not a Two State Folder.

Single molecule force spectroscopy is a useful technique for investigating mechanically induced protein unfolding and refolding under reduced forces by monitoring the end-to-end distance of the protein. The data is often interpreted via a "two-state" model based on the assumption that the end-to-end distance alone is a good reaction coordinate and the thermodynamic behavior is then ascribed to the free energy as a function of this one reaction coordinate. In this paper, we determined the free energy surface (PMF) of GB1 protein from atomistic simulations in explicit solvent under different applied forces as a function of two collective variables (the end-to-end-distance, and the fraction of native contacts ρ). The calculated 2-d free energy surfaces exhibited several distinct states, or basins, mostly visible along the ρ coordinate. Brownian dynamics (BD) simulations on the smoothed free energy surface show that the protein visits a metastable molten globule state and is thus a three state folder, not the two state folder inferred using the end-to-end distance as the sole reaction coordinate. This study lends support to recent experiments that suggest that GB1 is not a two-state folder.

Trifluoroethanol-induced conformational transition of the C-terminal sterile alpha motif (SAM) of human p73

The alpha splice variant of p73 (p73α), a homologue of the tumour suppressor p53, has at its C terminus a sterile alpha motif (SAM); this domain, SAMp73, is involved in lipid binding and it is thought to mediate in protein-protein interactions. As SAMp73 is a 68-residue-long helical bundle, it could be a good model to study the (2,2,2-trifluoroethanol) TFE-induced conformational transitions of α-helical proteins. Furthermore, as SAMp73 binds to lipids through a well-known polypeptide patch, we can test whether TFE is a good mimic of lipids and membranes. To address those questions, we used several biophysical probes, namely, fluorescence, circular dichroism, 1D, 2D and 3D-NMR spectroscopies, and dynamic light scattering. The TFE-induced conformational transition of SAMp73 was complex, involving several species as detected by the biophysical probes. The last TFE-induced transition occurred at a concentration of TFE of ∼20% (v/v), where the protein lost its compactness. None of those TFE-induced species accumulated during the two-state folding of SAMp73 in aqueous solution. The final state at 40% TFE was highly helical, but its structure was not rigid. For SAMp73, TFE did not properly mimic a membrane-like environment, since at very low TFE concentrations, other residues, together with those known to interact with lipids, were also affected by the co-solvent. Comparison with studies on isolated peptides, comprising the helical regions of SAMp73, suggests that peptides were good models of the intact protein in TFE.

Common folding processes of mini-proteins: Partial formations of secondary structures initiate the immediate protein folding.

The folding processes of mini-proteins (FSD-EY/FBPWW28 domain) were computationally investigated by an enhanced conformational sampling method. Through the analyses of trajectories, these mini-proteins had multiple folding pathways different from a simple two-state folding, and the multiple folding processes were initiated by partial formations of secondary structures. These findings can be used to understand the folding of large proteins, that is, which secondary structures are partially folded in the early process, and how the remaining parts are sequentially folded. It is found that FSD-EY (α/β topology) folds by a simple diffusion-collision mechanism, while the folding process of the FBPWW28 domain (all-β topology) requires a modification of the diffusion-collision theory to adequately treat the coil-sheet transition for the β sheet formation. © 2017 Wiley Periodicals, Inc.

Is Protein Single Molecule Dynamics under Force Described by Two or More States?

Single molecule force spectroscopy is a useful technique for investigating mechanically induced protein unfolding and refolding under reduced forces by monitoring the end-to-end distance of the protein. The data is often interpreted via a “two-state” model based on the assumption that the end-to-end distance alone is a good reaction coordinate and the thermodynamic behavior is then ascribed to the free energy as a function of this one reaction coordinate. In this paper, we determined the free energy surface (PMF) of GB1 protein from atomistic simulations in explicit solvent under different applied forces as a function of two collective variables (the end-to-end-distance, and the fraction of native contacts ρ). The calculated 2-d free energy surfaces exhibited several distinct states, or basins, mostly visible along the ρ coordinate. Brownian dynamics (BD) simulations on the smoothed free energy surface show that the protein visits a metastable molten globule state and is thus a three state folder, not the two state folder inferred using the end-to-end distance as the sole reaction coordinate. These BD simulations reproduce the unfolding and collapse-refolding patterns observed in the force-clamp experiments. This study lends support to recent experiments that suggest that GB1 is not a two-state folder.

Complex Dynamics in Single Molecule Force Spectroscopy from Simple Simulation Models

Single molecule force spectroscopy studies of protein folding dynamics have yielded unprecedented insight into the mechanisms of the folding reaction. Still, they usually allow for a very limited description of the process of interest (e.g. based on the reaction coordinate that is probed directly). Molecular simulations --even using very simplified models-- can give a much enriched view of the underlying molecular events. In a recent example, we explain the observations on a prototypical two-state folding protein, the cold shock protein, that in the AFM turns out to exhibit multiple intermediate states. Using a very simple topology based model for the protein we show that in fact this should not come as a surprise. Indeed, force can selectively stabilize different protein intermediates and act in a very different way compared to that of chemical denaturants.

D19S Mutation of the Cationic, Cysteine-Rich Protein PAF: Novel Insights into Its Structural Dynamics, Thermal Unfolding and Antifungal Function.

The cysteine-rich, cationic, antifungal protein PAF is abundantly secreted into the culture supernatant of the filamentous Ascomycete Penicillium chrysogenum. The five β-strands of PAF form a compact β-barrel that is stabilized by three disulphide bonds. The folding of PAF allows the formation of four surface-exposed loops and distinct charged motifs on the protein surface that might regulate the interaction of PAF with the sensitive target fungus. The growth inhibitory activity of this highly stable protein against opportunistic fungal pathogens provides great potential in antifungal drug research. To understand its mode of action, we started to investigate the surface-exposed loops of PAF and replaced one aspartic acid at position 19 in loop 2 that is potentially involved in PAF active or binding site, with a serine (Asp19 to Ser19). We analysed the overall effects, such as unfolding, electrostatic changes, sporadic conformers and antifungal activity when substituting this specific amino acid to the fairly indifferent amino acid serine. Structural analyses revealed that the overall 3D solution structure is virtually identical with that of PAF. However, PAFD19S showed slightly increased dynamics and significant differences in the surface charge distribution. Thermal unfolding identified PAFD19S to be rather a two-state folder in contrast to the three-state folder PAF. Functional comparison of PAFD19S and PAF revealed that the exchange at residue 19 caused a dramatic loss of antifungal activity: the binding and internalization of PAFD19S by target cells was reduced and the protein failed to trigger an intracellular Ca2+ response, all of which are closely linked to the antifungal toxicity of PAF. We conclude that the negatively charged residue Asp19 in loop 2 is essential for full function of the cationic protein PAF.

Synthesis and NMR Analysis of a Conformationally Controlled β-Turn Mimetic Torsion Balance.

The molecular torsion balance concept was applied to a new conformationally controlled scaffold and synthesized to accurately evaluate pairwise amino acid interactions in an antiparallel β-sheet motif. The scaffold's core design combines (ortho-tolyl)amide and o,o,o'-trisubstituted biphenyl structural units to provide a geometry better-suited for intramolecular hydrogen bonding. Like the dibenzodiazocine hinge of the traditional torsion balance, the (ortho-tolyl)amide unit offers restricted rotation around an N-aryl bond. The resulting two-state folding model is a powerful template for measuring hydrogen bond stability between two competing sequences. The aim of this study was to improve the alignment between the amino acid sequences attached to the upper and lower aromatic rings in order to promote hydrogen bond formation at the correct distance and antiparallel orientation. Bromine substituents were introduced ortho to the upper side chains and compared to a control to test our hypothesis. Hydrogen bond formation has been identified between the NH amide proton of the upper side chain (proton donor) and glycine acetamide of the lower side chain (proton acceptor).