未折叠蛋白质的普遍初始热力学亚稳态

探索和理解蛋白质折叠问题一直是分子生物学、结构生物学和生物物理学的终极挑战.未折叠的蛋白质应该存在一种普遍初始热力学亚稳态,否则无法解释蛋白质是如何在剧烈的热振动干扰下完成快速精确折叠的.本文通过分析水溶液环境和蛋白质折叠的相关性,揭示了一种由水分子屏蔽效应引起的未折叠蛋白质的普遍初始热力学亚稳态,该亚稳态的存在是水溶液环境中水分子的物理性质决定,并赋予未折叠蛋白质抵抗热扰动和避免错误折叠的能力.我们通过研究已发表的实验数据和建立分子模型,找到了该初始热力学亚稳态存在的相关证据,并推测了该亚稳态导致蛋白质精确折叠的相关物理学机制.     

RNA and protein folding: common themes and variations

Visualizing the navigation of an ensemble of unfolded molecules through the bumpy energy landscape in search of the native state gives a pictorial view of biomolecular folding. This picture, when combined with concepts in polymer theory, provides a unified theory of RNA and protein folding. Just as for proteins, the major folding free energy barrier for RNA scales sublinearly with the number of nucleotides, which allows us to extract the elusive prefactor for RNA folding. Several folding scenarios can be anticipated by considering variations in the energy landscape that depend on sequence, native topology, and external conditions. RNA and protein folding mechanism can be described by the kinetic partitioning mechanism (KPM) according to which a fraction (Phi) of molecules reaches the native state directly, whereas the remaining fraction gets kinetically trapped in metastable conformations. For two-state folders Phi approximately 1. Molecular chaperones are recruited to assist protein folding whenever Phi is small. We show that the iterative annealing mechanism, introduced to describe chaperonin-mediated folding, can be generalized to understand protein-assisted RNA folding. The major differences between the folding of proteins and RNA arise in the early stages of folding. For RNA, folding can only begin after the polyelectrolyte problem is solved, whereas protein collapse requires burial of hydrophobic residues. Cross-fertilization of ideas between the two fields should lead to an understanding of how RNA and proteins solve their folding problems.     

Observation of Continuous Contraction and a Metastable Misfolded State during the Collapse and Folding of a Small Protein

To obtain proper insight into how structure develops during a protein folding reaction, it is necessary to understand the nature and mechanism of the polypeptide chain collapse reaction, which marks the initiation of folding. Here, the time-resolved fluorescence resonance energy transfer technique, in which the decay of the fluorescence light intensity with time is used to determine the time evolution of the distribution of intra-molecular distances, has been utilized to study the folding of the small protein, monellin. It is seen that when folding begins, about one-third of the protein molecules collapse into a molten globule state (I_MG), from which they relax by continuous further contraction to transit to the native state. The larger fraction gets trapped into a metastable misfolded state. Exit from this metastable state occurs via collapse to the lower free energy I_MG state. This exit is slow, on a time scale of seconds, because of activation energy barriers. The trapped misfolded molecules as well as the I_MG molecules contract continuously and slowly as structure develops. A phenomenological model of Markovian evolution of the polymer chain undergoing folding, incorporating these features, has been developed, which fits well the experimentally observed time evolution of distance distributions. The observation that the “wrong turn” to the misfolded state has not been eliminated by evolution belies the common belief that the folding of functional protein sequences is very different from that of a random heteropolymer, and the former have been selected by evolution to fold quickly.     

Thermodynamics and kinetics of protein folding: an evolutionary perspective

This article appeals to an evolutionary model which postulates that primordial proteins were described by small polypeptide chains which (i) lack disulfide bridges, and (ii) display slow folding rates with multi-state kinetics, to determine relations between structural properties of proteins and their folding kinetics. We parameterize the energy landscape of proteins in terms of thermodynamic activation variables. The model studies evolutionary changes in these thermodynamic parameters, and we invoke relations between these activation variables and structural properties of the protein to predict the following correspondence between protein structure and folding kinetics. 1. Proteins with inter- and intra-chain disulfide bridges: large variability in both folding rates and stability of intermediates, multi-state kinetics. 2. Proteins which lack inter and intra-chain disulfide bridges. 2.1 Single-domain chains: fast folding rates; unstable intermediates; two-state kinetics. 2.2 Multi-domain monomers: intermediate rates; metastable intermediates; multi-state kinetics. 2.3 Multi-domain oligomers: slow rates; metastable intermediates; multi-state kinetics. The evolutionary model thus provides a kinetic characterization of one important subfamily of proteins which we describe by the following properties: Folding dynamics of single-domain proteins which lack disulfide bridges are described by two-state kinetics. Folding rate of this class of proteins is positively correlated with the thermodynamic stability of the folded state.     

How cooperative are protein folding and unfolding transitions?

A thermodynamically and kinetically simple picture of protein folding envisages only two states, native (N) and unfolded (U), separated by a single activation free energy barrier, and interconverting by cooperative two‐state transitions. The folding/unfolding transitions of many proteins occur, however, in multiple discrete steps associated with the formation of intermediates, which is indicative of reduced cooperativity. Furthermore, much advancement in experimental and computational approaches has demonstrated entirely non‐cooperative (gradual) transitions via a continuum of states and a multitude of small energetic barriers between the N and U states of some proteins. These findings have been instrumental towards providing a structural rationale for cooperative versus noncooperative transitions, based on the coupling between interaction networks in proteins. The cooperativity inherent in a folding/unfolding reaction appears to be context dependent, and can be tuned via experimental conditions which change the stabilities of N and U. The evolution of cooperativity in protein folding transitions is linked closely to the evolution of function as well as the aggregation propensity of the protein. A large activation energy barrier in a fully cooperative transition can provide the kinetic control required to prevent the accumulation of partially unfolded forms, which may promote aggregation. Nevertheless, increasing evidence for barrier‐less “downhill” folding, as well as for continuous “uphill” unfolding transitions, indicate that gradual non‐cooperative processes may be ubiquitous features on the free energy landscape of protein folding.     

Ion-exchange resins greatly facilitate refolding of like-charged proteins at high concentrations

Protein refolding is a crucial step for the production of therapeutic proteins expressed in bacteria as inclusion bodies. In vitro protein refolding is severely impeded by the aggregation of folding intermediates during the folding process, so inhibition of the aggregation is the most effective approach to high‐efficiency protein refolding. We have herein found that electrostatic repulsion between like‐charged protein and ion exchange gel beads can greatly suppress the aggregation of folding intermediates, leading to the significant increase of native protein recovery. This finding is extensively demonstrated with three different proteins and four kinds of ion‐exchange resins when the protein and ion‐exchange gel are either positively or negatively charged at the refolding conditions. It is remarkable that the enhancing effect is significant at very high protein concentrations, such as 4 mg/mL lysozyme (positively charged) and 2 mg/mL bovine serum albumin (negatively charged). Moreover, the folding kinetics is not compromised by the presence of the resins, so fast protein refolding is realized at high protein concentrations. It was not realistic by any other approaches. The working mechanism of the like‐charged resin is considered due to the charge repulsion that could induce oriented alignment of protein molecules near the charged surface, leading to the inhibition of protein aggregation. The molecular crowding effect induced by the charge repulsion may also contribute to accelerating protein folding. The refolding method with like‐charged ion exchangers is simple to perform, and the key material is easy to separate for recycling. Moreover, because ion exchangers can work as adsorbents of oppositely charged impurities, an operation of simultaneous protein refolding and purification is possible. All the characters are desirable for preparative refolding of therapeutic proteins expressed in bacteria as inclusion bodies. Bioeng. 2011; 108:1068–1077. © 2010 Wiley Periodicals, Inc.     

Folding nucleus structure persists in thermally‐aggregated FGF‐1

An efficient protein‐folding pathway leading to target structure, and the avoidance of aggregation, is essential to protein evolution and de novo design; however, design details to achieve efficient folding and avoid aggregation are poorly understood. We report characterization of the thermally‐induced aggregate of fibroblast growth factor‐1 (FGF‐1), a small globular protein, by solid‐state NMR. NMR spectra are consistent with residual structure in the aggregate and provide evidence of a structured region that corresponds to the region of the folding nucleus. NMR data on aggregated FGF‐1 also indicate the presence of unstructured regions that exhibit hydration‐dependent dynamics and suggest that unstructured regions of aggregated FGF‐1 lie outside the folding nucleus. Since it is known that regions outside the folding nucleus fold late in the folding pathway, we postulate that these regions unfold early in the unfolding pathway and that the partially folded state is more prone to intermolecular aggregation. This interpretation is further supported by comparison with a designed protein that shares the same FGF‐1 folding nucleus sequence, but has different 1° structure outside the folding nucleus, and does not thermally aggregate. The results suggest that design of an efficient folding nucleus, and the avoidance of aggregation in the folding pathway, are potentially separable design criteria – the latter of which could principally focus upon the physicochemical properties of 1° structure outside the folding nucleus.     

Iterative annealing mechanism explains the functions of the GroEL and RNA chaperones

Molecular chaperones are ATP‐consuming machines, which facilitate the folding of proteins and RNA molecules that are kinetically trapped in misfolded states. Unassisted folding occurs by the kinetic partitioning mechanism according to which folding to the native state, with low probability as well as misfolding to one of the many metastable states, with high probability, occur rapidly. GroEL is an all‐purpose stochastic machine that assists misfolded substrate proteins to fold. The RNA chaperones such as CYT‐19, which are ATP‐consuming enzymes, help the folding of ribozymes that get trapped in metastable states for long times. GroEL does not interact with the folded proteins but CYT‐19 disrupts both the folded and misfolded ribozymes. The structures of GroEL and RNA chaperones are strikingly different. Despite these differences, the iterative annealing mechanism (IAM) quantitatively explains all the available experimental data for assisted folding of proteins and ribozymes. Driven by ATP binding and hydrolysis and GroES binding, GroEL undergoes a catalytic cycle during which it samples three allosteric states, T (apo), R (ATP bound), and R^″ (ADP bound). Analyses of the experimental data show that the efficiency of the GroEL–GroES machinery and mutants is determined by the resetting rate k _{R ″ → T} , which is largest for the wild‐type (WT) GroEL. Generalized IAM accurately predicts the folding kinetics of Tetrahymena ribozyme and its variants. Chaperones maximize the product of the folding rate and the steady‐state native state fold by driving the substrates out of equilibrium. Neither the absolute yield nor the folding rate is optimized.     

The nascent polypeptide‐associated complex is a key regulator of proteostasis

The adaptation of protein synthesis to environmental and physiological challenges is essential for cell viability. Here, we show that translation is tightly linked to the protein‐folding environment of the cell through the functional properties of the ribosome bound chaperone NAC (nascent polypeptide‐associated complex). Under non‐stress conditions, NAC associates with ribosomes to promote translation and protein folding. When proteostasis is imbalanced, NAC relocalizes from a ribosome‐associated state to protein aggregates in its role as a chaperone. This results in a functional depletion of NAC from the ribosome that diminishes translational capacity and the flux of nascent proteins. Depletion of NAC from polysomes and re‐localisation to protein aggregates is observed during ageing, in response to heat shock and upon expression of the highly aggregation‐prone polyglutamine‐expansion proteins and Aβ‐peptide. These results demonstrate that NAC has a central role as a proteostasis sensor to provide the cell with a regulatory feedback mechanism in which translational activity is also controlled by the folding state of the cellular proteome and the cellular response to stress.     

How the folding rates of two‐ and multistate proteins depend on the amino acid properties

Proteins fold by either two‐state or multistate kinetic mechanism. We observe that amino acids play different roles in different mechanism. Many residues that are easy to form regular secondary structures (α helices, β sheets and turns) can promote the two‐state folding reactions of small proteins. Most of hydrophilic residues can speed up the multistate folding reactions of large proteins. Folding rates of large proteins are equally responsive to the flexibility of partial amino acids. Other properties of amino acids (including volume, polarity, accessible surface, exposure degree, isoelectric point, and phase transfer energy) have contributed little to folding kinetics of the proteins. Cysteine is a special residue, it triggers two‐state folding reaction and but inhibits multistate folding reaction. These findings not only provide a new insight into protein structure prediction, but also could be used to direct the point mutations that can change folding rate. Proteins 2014; 82:2375–2382. © 2014 Wiley Periodicals, Inc.     

Nucleation‐based prediction of the protein folding rate and its correlation with the folding nucleus size

The problem of protein self‐organization is in the focus of current molecular biology studies. Although the general principles are understood, many details remain unclear. Specifically, protein folding rates are of interest because they dictate the rate of protein aggregation which underlies many human diseases. Here we offer predictions of protein folding rates and their correlation with folding nucleus sizes. We calculated free energies of the transition state and sizes of folding nuclei for 84 proteins and peptides whose other parameters were measured at the point of thermodynamic equilibrium between their unfolded and native states. We used the dynamic programming method where each residue was considered to be either as folded as in its native state or completely disordered. The calculated and measured folding rates showed a good correlation at the temperature mid‐transition point (the correlation coefficient was 0.75). Also, we pioneered in demonstrating a moderate (‐0.57) correlation coefficient between the calculated sizes of folding nuclei and the folding rates. Predictions made by different methods were compared. The established good correlation between the estimated free energy barrier and the experimentally found folding rate of each studied protein/peptide indicates that our model gives reliable results for the considered data set. Proteins 2012; © 2012 Wiley Periodicals, Inc.     

Thermodynamic modulation of folding and aggregation energy landscape by DNA binding of functional domains of TDP-43

TDP-43 is a vital nucleic acid binding protein which forms stress-induced aberrant aggregates in around 97% cases of ALS, a fatal neurodegenerative disease. The functional tandem RRM domain of the protein (TDP-43^tRRM) has been shown to undergo amyloid-like aggregation under stress in a pH-dependent fashion. However, the underlying thermodynamic and molecular basis of aggregation and how the energy landscape of folding, stability, and aggregation are coupled and modulated by nucleic acid binding is poorly understood. Here, we show that the pH stress thermodynamically destabilizes the native protein and systematically populates the unfolded-like aggregation-prone molecules which leads to amyloid-like aggregation. We observed that specific DNA binding inhibits aggregation and populates native-like compact monomeric state even under low-pH stress as measured by circular dichroism, ANS binding, size exclusion chromatography, and transmission electron microscopy. We show that DNA-binding thermodynamically stabilizes and populates the native state even under stress and reduces the population of unfolded-like aggregation-prone molecules which leads to systematic aggregation inhibition. Our results suggest that thermodynamic modulation of the folding and aggregation energy landscape by nucleic-acid-like molecules could be a promising approach for effective therapeutic intervention in TDP-43-associated proteinopathies.     

Retarded protein folding of deficient human alpha 1-antitrypsin D256V and L41P variants

alpha(1)-Antitrypsin is the most abundant protease inhibitor in plasma and is the archetype of the serine protease inhibitor superfamily. Genetic variants of human alpha(1)-antitrypsin are associated with early-onset emphysema and liver cirrhosis. However, the detailed molecular mechanism for the pathogenicity of most variant alpha(1)-antitrypsin molecules is not known. Here we examined the structural basis of a dozen deficient alpha(1)-antitrypsin variants. Unlike most alpha(1)-antitrypsin variants, which were unstable, D256V and L41P variants exhibited extremely retarded protein folding as compared with the wild-type molecule. Once folded, however, the stability and inhibitory activity of these variant proteins were comparable to those of the wild-type molecule. Retarded protein folding may promote protein aggregation by allowing the accumulation of aggregation-prone folding intermediates. Repeated observations of retarded protein folding indicate that it is an important mechanism causing alpha(1)-antitrypsin deficiency by variant molecules, which have to fold into the metastable native form to be functional.     

Energy Landscapes of Protein Aggregation and Conformation Switching in Intrinsically Disordered Proteins

The protein folding problem was apparently solved recently by the advent of a deep learning method for protein structure prediction called AlphaFold. However, this program is not able to make predictions about the protein folding pathways. Moreover, it only treats about half of the human proteome, as the remaining proteins are intrinsically disordered or contain disordered regions. By definition these proteins differ from natively folded proteins and do not adopt a properly folded structure in solution. However these intrinsically disordered proteins (IDPs) also systematically differ in amino acid composition and uniquely often become folded upon binding to an interaction partner. These factors preclude solving IDP structures by current machine-learning methods like AlphaFold, which also cannot solve the protein aggregation problem, since this meta-folding process can give rise to different aggregate sizes and structures. An alternative computational method is provided by molecular dynamics simulations that already successfully explored the energy landscapes of IDP conformational switching and protein aggregation in multiple cases. These energy landscapes are very different from those of ‘simple’ protein folding, where one energy funnel leads to a unique protein structure. Instead, the energy landscapes of IDP conformational switching and protein aggregation feature a number of minima for different competing low-energy structures. In this review, I discuss the characteristics of these multifunneled energy landscapes in detail, illustrated by molecular dynamics simulations that elucidated the underlying conformational transitions and aggregation processes.     

A Fluorescence‐Based Sensor Assay that Monitors General Protein Aggregation in Human Cells

Protein conformational disorders are characterized by disruption of protein folding and toxic accumulation of protein aggregates. Here we describe a sensitive and simple method to follow and monitor general protein aggregation in human cells. Heat shock protein 27 (HSP27) is an oligomeric small heat shock protein that binds and keeps unfolded proteins in a folding competent state. This high specificity of HSP27 for aggregated proteins can be explored to monitor aggregation in living cells by fusing it to a fluorescent protein as Green Fluorescent Protein (GFP). We have constructed a HeLa stable cell line expressing a HSP27:GFP chimeric reporter protein and after validation, this stable cell line is exposed to different agents that interfere with proteostasis, namely Arsenite, MG132, and Aβ‐peptide. Exposure to proteome destabilizers lead to re‐localization of HSP27:GFP fluorescence to foci, confirming that our reporter system is functional and can be used to detect and follow protein aggregation in living cells. This reporter is a valuable tool to setup wide‐genetic screens to identify genes and pathways involved in protein misfolding and aggregation.     

Analyzing the effect of homogeneous frustration in protein folding

The energy landscape theory has been an invaluable theoretical framework in the understanding of biological processes such as protein folding, oligomerization, and functional transitions. According to the theory, the energy landscape of protein folding is funneled toward the native state, a conformational state that is consistent with the principle of minimal frustration. It has been accepted that real proteins are selected through natural evolution, satisfying the minimum frustration criterion. However, there is evidence that a low degree of frustration accelerates folding. We examined the interplay between topological and energetic protein frustration. We employed a C_α structure‐based model for simulations with a controlled nonspecific energetic frustration added to the potential energy function. Thermodynamics and kinetics of a group of 19 proteins are completely characterized as a function of increasing level of energetic frustration. We observed two well‐separated groups of proteins: one group where a little frustration enhances folding rates to an optimal value and another where any energetic frustration slows down folding. Protein energetic frustration regimes and their mechanisms are explained by the role of non‐native contact interactions in different folding scenarios. These findings strongly correlate with the protein free‐energy folding barrier and the absolute contact order parameters. These computational results are corroborated by principal component analysis and partial least square techniques. One simple theoretical model is proposed as a useful tool for experimentalists to predict the limits of improvements in real proteins.Proteins 2013; 81:1727–1737. © 2013 Wiley Periodicals, Inc.     

Existence of Molten Globule State in Homocysteine-Induced Protein Covalent Modifications

Homocysteine thiolactone is a toxic metabolite produced from homocysteine by amino-acyl t-RNA synthetase in error editing reaction. The basic cause of toxicity of homocysteine thiolactone is believed to be due to the adduct formation with lysine residues (known as protein N-homocysteinylation) leading to protein aggregation and loss of enzyme function. There was no data available until now that showed the effect of homocysteine thiolactone on the native state structural changes that led to aggregate formation. In the present study we have investigated the time dependent structural changes due to homocysteine thiolactone induced modifications on three different proteins having different physico-chemical properties (cytochrome-c, lysozyme and alpha lactalbumin). We discovered that N-homocysteinylation leads to the formation of molten globule state—an important protein folding intermediate in the protein folding pathway. We also found that the formation of the molten globule state might be responsible for the appearance of aggregate formation. The study indicates the importance of protein folding intermediate state in eliciting the homocysteine thiolactone toxicity.     

Mechanism for the alpha-helix to beta-hairpin transition

The aggregation of alpha-helix-rich proteins into beta-sheet-rich amyloid fibrils is associated with fatal diseases, such as Alzheimer's disease and prion disease. During an aggregation process, protein secondary structure elements-alpha-helices-undergo conformational changes to beta-sheets. The fact that proteins with different sequences and structures undergo a similar transition on aggregation suggests that the sequence nonspecific hydrogen bond interaction among protein backbones is an important factor. We perform molecular dynamics simulations of a polyalanine model, which is an alpha-helix in its native state and observe a metastable beta-hairpin intermediate. Although a beta-hairpin has larger potential energy than an alpha-helix, the entropy of a beta-hairpin is larger because of fewer constraints imposed by the hydrogen bonds. In the vicinity of the transition temperature, we observe the interconversion of the alpha-helix and beta-sheet states via a random coil state. We also study the effect of the environment by varying the relative strength of side-chain interactions for a designed peptide-an alpha-helix in its native state. For a certain range of side-chain interaction strengths, we find that the intermediate beta-hairpin state is destabilized and even disappears, suggesting an important role of the environment in the aggregation propensity of a peptide.