Lattice simulations of aggregation funnels for protein folding.

A computer model of protein aggregation competing with productive folding is proposed. Our model adapts techniques from lattice Monte Carlo studies of protein folding to the problem of aggregation. However, rather than starting with a single string of residues, we allow independently folding strings to undergo collisions and consider their interactions in different orientations. We first present some background into the nature and significance of protein aggregation and the use of lattice Monte Carlo simulations in understanding other aspects of protein folding. The results of a series of simulation experiments involving simple versions of the model illustrate the importance of considering aggregation in simulations of protein folding and provide some preliminary understanding of the characteristics of the model. Finally, we discuss the value of the model in general and of our particular design decisions and experiments. We conclude that computer simulation techniques developed to study protein folding can provide insights into protein aggregation, and that a better understanding of aggregation may in turn provide new insights into and constraints on the more general protein folding problem.     

Folding pathway of the b1 domain of protein G explored by multiscale modeling

The understanding of the folding mechanisms of single-domain proteins is an essential step in the understanding of protein folding in general. Recently, we developed a mesoscopic CA-CB side-chain protein model, which was successfully applied in protein structure prediction, studies of protein thermodynamics, and modeling of protein complexes. In this research, this model is employed in a detailed characterization of the folding process of a simple globular protein, the B1 domain of IgG-binding protein G (GB1). There is a vast body of experimental facts and theoretical findings for this protein. Performing unbiased, ab initio simulations, we demonstrated that the GB1 folding proceeds via the formation of an extended folding nucleus, followed by slow structure fine-tuning. Remarkably, a subset of native interactions drives the folding from the very beginning. The emerging comprehensive picture of GB1 folding perfectly matches and extends the previous experimental and theoretical studies.     

Matching simulation and experiment: a new simplified model for simulating protein folding.

Simulations of simplified protein folding models have provided much insight into solving the protein folding problem. We propose here a new off-lattice bead model, capable of simulating several different fold classes of small proteins. We present the sequence for an alpha/beta protein resembling the IgG-binding proteins L and G. The thermodynamics of the folding process for this model are characterized using the multiple multihistogram method combined with constant-temperature Langevin simulations. The folding is shown to be highly cooperative, with chain collapse nearly accompanying folding. Two parallel folding pathways are shown to exist on the folding free energy landscape. One pathway contains an intermediate--similar to experiments on protein G, and one pathway contains no intermediates-similar to experiments on protein L. The folding kinetics are characterized by tabulating mean-first passage times, and we show that the onset of glasslike kinetics occurs at much lower temperatures than the folding temperature. This model is expected to be useful in many future contexts: investigating questions of the role of local versus nonlocal interactions in various fold classes, addressing the effect of sequence mutations affecting secondary structure propensities, and providing a computationally feasible model for studying the role of solvation forces in protein folding.     

Folding-unfolding energy change of a simple sphere model protein and an energy landscape of the folding process

We have calculated the free energy of a spherical model of a protein or part of a protein generated in the way of protein folding. Two spherical models are examined; one is a homogeneous model consisting of only one residue type—hydrophobic. The other is a heterogeneous model consisting of two residue types—strong hydrophobic and weak hydrophobic. Both models show a folding transition state, and the latter model reproduces the trend of the experimental folded-unfolded energy change. The heterogeneous model suggests that in the folding process of a protein of more than 70 residues, a specific region of the protein folds first to form a stable region, then the other residues follow the folding process. The energy landscape of folding of a small protein is approximately a funnel model, whereas a flatter energy landscape is suggested for larger proteins of more than 55–70 residues. Proteins 33:408–416, 1998. © 1998 Wiley-Liss, Inc.     

Folding studies of two hydrostatic pressure sensitive proteins

High hydrostatic pressure combined with various spectroscopies is a powerful technique to study protein folding. An ideal model system for protein folding studies should have the following characteristics. (1) The protein should be sensitive to pressure, so that the protein can be unfolded under mild pressure. (2) The folding process of the protein should be easily modulated by several chemical or physical factors. (3) The folding process should be easily monitored by some spectroscopic parameters. Here, we summarized the pressure induced folding studies of two proteins isolated from spinach photosystem II, namely the 23-kDa and the 33-kDa protein. They have all the characteristics mention above and might be an ideal model protein system for pressure studies.     

Toward minimalist models of larger proteins: a ubiquitin-like protein

Our recently developed off-lattice bead model capable of simulating protein structures with mixed alpha/beta content has been extended to model the folding of a ubiquitin-like protein and provides a means for examining the more complex kinetics involved in the folding of larger proteins. Using trajectories generated from constant-temperature Langevin dynamics simulations and sampling with the multiple multi-histogram method over five-order parameters, we are able to characterize the free energy landscape for folding and find evidence for folding through compact intermediates. Our model reproduces the observation that the C-terminus loop structure in ubiquitin is the last to fold in the folding process and most likely plays a spectator role in the folding kinetics. The possibility of a productive metastable intermediate along the folding pathway consisting of collapsed states with no secondary structure, and of intermediates or transition structures involving secondary structural elements occurring early in the sequence, is also supported by our model. The kinetics of folding remain multi-exponential below the folding temperature, with glass-like kinetics appearing at T/T(f) approximately 0.86. This new physicochemical model, designed to be predictive, helps validate the value of modeling protein folding at this level of detail for genomic-scale studies, and motivates further studies of other protein topologies and the impact of more complex energy functions, such as the addition of solvation forces.     

基于HP模型的蛋白质折叠问题的研究

基于蛋白质二维HP模型提出改进的遗传算法对真实蛋白质进行计算机折叠模拟。结果显示疏水能量函数最小值的蛋白质构象对应含疏水核心的稳定结构,疏水作用在蛋白质折叠中起主要作用。研究表明二维HP模型在蛋白质折叠研究中是可行的和有效的并为进一步揭示蛋白质折叠机理提供重要参考信息。     

用计算机折叠蛋白质:我们走了多远?(英文)

蛋白质折叠是现代科学最具挑战性的难题之一。Anfinsen的先驱性工作已经过去数十年了,我们对蛋白质折叠的机理仍不甚了然。但值得庆幸的是,科学工作者们在解析几个模型蛋白质的折叠机理上已经取得了明显的进展。这篇综述将主要回顾在蛋白质折叠的计算研究方面的进展。受益于计算机技术的迅猛发展,在1998年首次实现了一个微秒的全原子水平的蛋白质折叠,从2000年开始,folding@home将全球分布式计算技术应用于蛋白质折叠。在经历了数十年艰苦的努力之后,天然蛋白质亚埃精度的折叠终于在2007年成功实现。近年来,一种全新的概念开始出现,人们越来越多地用网络来解释蛋白质折叠的机理。随着分子力场的不断完善,分子模拟将会在解析蛋白质折叠的机理上起到越来越重要的作用。     

Global and local prediction of protein folding rates based on sequence autocorrelation information

To understand the folding behavior of proteins is an important and challenging problem in modern molecular biology. In the present investigation, a large number of features representing protein sequences were developed based on sequence autocorrelation weighted by properties of amino acid residues. Genetic algorithm (GA) combined with multiple linear regression (MLR) was employed to select significant features related to protein folding rates, and to build global predictive model. Moreover, local lazy regression (LLR) method was also used to predict the protein folding rates. The obtained results indicated that LLR performed much better than the global MLR model. The important properties of amino acid residues affecting protein folding rates were also analyzed. The results of this study will be helpful to understand the mechanism of protein folding. Our results also demonstrate that the features of amino acid sequence autocorrelation is effective in representing the relationship between protein sequence and folding rates, and the local method is a powerful tool to predict the protein folding rates.     

Gatekeepers in the ribosomal protein s6: thermodynamics, kinetics, and folding pathways revealed by a minimalist protein model

We investigate the effect of structural gatekeepers on the folding of the ribosomal protein S6. Folding thermodynamics and early refolding kinetics are studied for this system utilizing computer simulations of a minimalist protein model. When gatekeepers are eliminated, the thermodynamic signature of a folding intermediate emerges, and a marked decrease in folding efficiency is observed. We explain the prerequisites that determine the "strength" of a given gatekeeper. The investigated gatekeepers are found to have distinct functions, and to guide the folding and time-dependent packing of non-overlapping secondary structure elements in the protein. Gatekeepers avoid kinetic traps during folding by favoring the formation of "productive topologies" on the way to the native state. The trends in folding rates in the presence/absence of gatekeepers observed for our minimalist model of S6 are in very good agreement with experimental data on this protein.     

Diffusion-collision of foldons elucidates the kinetic effects of point mutations and suggests control strategies of the folding process of helical proteins

In this article we use mutation studies as a benchmark for a minimal model of the folding process of helical proteins. The model ascribes a pivotal role to the collisional dynamics of a few crucial residues (foldons) and predicts the folding rates by exploiting information drawn from the protein sequence. We show that our model rationalizes the effects of point mutations on the kinetics of folding. The folding times of two proteins and their mutants are predicted. Stability and location of foldons have a critical role as the determinants of protein folding. This allows us to elucidate two main mechanisms for the kinetic effects of mutations. First, it turns out that the mutations eliciting the most notable effects alter protein stability through stabilization or destabilization of the foldons. Secondly, the folding rate is affected via a modification of the foldon topology by those mutations that lead to the birth or death of foldons. The few mispredicted folding rates of some mutants hint at the limits of the current version of the folding model proposed in the present article. The performance of our folding model declines in case the mutated residues are subject to strong long-range forces. That foldons are the critical targets of mutation studies has notable implications for design strategies and is of particular interest to address the issue of the kinetic regulation of single proteins in the general context of the overall dynamics of the interactome.     

Lysozyme folded in silico according to the limited conformational sub-space

The conformational sub-space oriented on early-stage protein folding is applied to lysozyme folding. The part of the Ramachandran map distinguished on the basis of a geometrical model of the polypeptide chain limited to the mutual orientation of the peptide bond planes is shown to deliver the initial structure of the polypeptide for the energy minimization procedure in the ab initio model of protein folding prediction. Two forms of energy minimization and molecular dynamics simulation procedures were applied to the assumed early-stage protein folding of lysozyme. One of them included the disulphide bond system and the other excluded it. The post-energy-minimization and post-dynamics structures were compared using RMS-D and non-bonding contact maps to estimate the degree of approach to the native, target structure of the protein molecule obtained using the limited conformational sub-space for the early stage of folding.     

Lysozyme Folded In Silico According to the Limited Conformational Sub-space

Abstract

The conformational sub-space oriented on early-stage protein folding is applied to lysozyme folding. The part of the Ramachandran map distinguished on the basis of a geometrical model of the polypeptide chain limited to the mutual orientation of the peptide bond planes is shown to deliver the initial structure of the polypeptide for the energy minimization procedure in the ab initio model of protein folding prediction. Two forms of energy minimization and molecular dynamics simulation procedures were applied to the assumed early-stage protein folding of lysozyme. One of them included the disulphide bond system and the other excluded it. The post-energy-minimization and post-dynamics structures were compared using RMS-D and non-bonding contact maps to estimate the degree of approach to the native, target structure of the protein molecule obtained using the limited conformational sub-space for the early stage of folding.     

Molecular basis of cooperativity in protein folding IV. Core: A general cooperative folding model

The cooperative nature of the protein folding process is independent of the characteristic fold and the specific secondary structure attributes of a globular protein. A general folding/unfolding model should, therefore, be based upon structural features that transcend the peculiarities of α-helices, β-sheets, and other structural motifs found in proteins. The studies presented in this paper suggest that a single structural characteristic common to all globular proteins is essential for cooperative folding. The formation of a partly folded state from the native state results in the exposure to solvent of two distinct regions: (1) the portions of the protein that are unfolded; and (2) the “complementary surfaces,” located in the regions of the protein that remain folded. The cooperative character of the folding/unfolding transition is determined largely by the energetics of exposing complementary surface regions to the solvent. By definition, complementary regions are present only in partly folded states; they are absent from the native and unfolded states. An unfavorable free energy lowers the probability of partly folded states and increases the cooperativity of the transition. In this paper we present a mathematical formulation of this behavior and develop a general cooperative folding/unfolding model, termed the “complementary region” (CORE) model. This model successfully reproduces the main properties of folding/unfolding transitions without limiting the number of partly folded states accessible to the protein, thereby permitting a systematic examination of the structural and solvent conditions under which intermediates become populated. It is shown that the CORE model predicts two-state folding/unfolding behavior, even though the two-state character is not assumed in the model. © 1993 Wiley-Liss, Inc.     

Kinetic model of BiP- and PDI-mediated protein folding and assembly

A mechanism for heavy chain binding protein (BiP)- and protein disulfide isomerase (PDI)- mediated protein folding and assembly has been proposed. It considers BiP chaperoning action and PDI catalytic activity. A kinetic model has been developed based on the proposed mechanism. The model was used for quantifying the influence of polypeptide concentration and ratio, and the effect of BiP and PDI concentration on the kinetics of folding and assembly. An optimum value for polypeptide concentration that minimizes assembly times was found, and different kinetic behaviors were identified for polypeptide concentrations higher or lower than the optimum. Pulse-chase experiments and the dependence of assembly time on unassembled polypeptides ratio predicted by the model are similar to those found during in vitro and in vivo folding and assembly of antibodies and human chorionic gonadotropin (hCG), as well as bovine pancreatic trypsin inhibitor (BPTI) folding. The model also explains the increase in folding and assembly rates during overexpression of BiP and PDI.     

蛋白质序列的混合特征值对折叠速率的影响

鉴于蛋白质折叠速率预测对研究其蛋白质功能的重要性,许多的科研工作者都开始对影响蛋白质折叠速率的因素进行研究。各种预测参数和方法被提出。利用蛋白质编码序列的不同特征参数,不同的二级结构及不同的折叠类的蛋白质对折叠速率的不同影响,我们选取蛋白质编码序列的新的特征值,即选取蛋白质序列的LZ复杂度,等电点等特征值。然后把这些特征值与20种氨基酸的属性αc、Cα、K0、Pβ、Ra、ΔASA、PI、ΔGhD、Nm、LZ、Mu、El融合,建立多元线性回归模型,并利用回归模型计算了13个全α类蛋白质、18个全β类蛋白质、13个混合类蛋白质和39个未分类蛋白质的ln（kf）与预测值之间的相关系数分别达到0.89、0.93、0.98、0.86。在Jack-knife方法的验证下发现在不同的结构中混合特征值与相应折叠速率有很好的相关性。结果表明,在蛋白质折叠过程中,蛋白质序列的LZ复杂度、等电点等特征值可能影响蛋白质的折叠速率及其结构。     

Protein modeling with reduced representation: statistical potentials and protein folding mechanism

A high resolution reduced model of proteins is used in Monte Carlo dynamics studies of the folding mechanism of a small globular protein, the B1 immunoglobulin-binding domain of streptococcal protein G. It is shown that in order to reproduce the physics of the folding transition, the united atom based model requires a set of knowledge-based potentials mimicking the short-range conformational propensities and protein-like chain stiffness, a model of directional and cooperative hydrogen bonds, and properly designed knowledge-based potentials of the long-range interactions between the side groups. The folding of the model protein is cooperative and very fast. In a single trajectory, a number of folding/unfolding cycles were observed. Typically, the folding process is initiated by assembly of a native-like structure of the C-terminal hairpin. In the next stage the rest of the four-ribbon beta-sheet folds. The slowest step of this pathway is the assembly of the central helix on the scaffold of the beta-sheet.     

Protein folding as a nonlinear nonequilibrium thermodynamic process

Biofurcational theory of protein folding is developed describing the process of formation of protein native structure as a sequence of non-equilibrium irreversible fluctuations, specific for a particular protein. The model gives explanation to all characteristic features of the folding process: stochastic mechanism, short time and precision of proteins self-assembling. A constructive role of entropy in the formation of a highly ordered structure out of disorder is discussed. A numerical method of a priori calculations of polypeptide structures basing on principles of bifractional folding model is presented.     

Minimalist protein model as a diagnostic tool for misfolding and aggregation

We propose a realistic coarse-grained protein model and a technique to "anchor" the model to available experimental data. We apply this procedure to characterize the effect of multiple mutations on the folding mechanism of protein S6. We show that the mutation of a few "gatekeeper" residues triggers significant changes on the folding landscape of S6. These results suggest that gatekeeper residues control the flexibility of critical regions of S6, that in turn regulates the delicate balance between folding and aggregation. Although obtained with a minimalist protein model, these results are fully consistent with experimental evidence and offer a clue to understand the interplay between folding and aggregation in protein S6.