蛋白质的折叠

重点介绍了蛋白质折叠的热力学控制学说和动力学控制学说,简单介绍了几种蛋白质折叠模型并分析了多肽链在体内进行快速折叠的原因。     

Chaperonin-mediated Protein Folding

We have been studying chaperonins these past twenty years through an initial discovery of an action in protein folding, analysis of structure, and elucidation of mechanism. Some of the highlights of these studies were presented recently upon sharing the honor of the 2013 Herbert Tabor Award with my early collaborator, Ulrich Hartl, at the annual meeting of the American Society for Biochemistry and Molecular Biology in Boston. Here, some of the major findings are recounted, particularly recognizing my collaborators, describing how I met them and how our great times together propelled our thinking and experiments.     

Microfluidic Mixers for Studying Protein Folding

The process by which a protein folds into its native conformation is highly relevant to biology and human health yet still poorly understood. One reason for this is that folding takes place over a wide range of timescales, from nanoseconds to seconds or longer, depending on the protein¹. Conventional stopped-flow mixers have allowed measurement of folding kinetics starting at about 1 ms. We have recently developed a microfluidic mixer that dilutes denaturant ~100-fold in ~8 μs². Unlike a stopped-flow mixer, this mixer operates in the laminar flow regime in which turbulence does not occur. The absence of turbulence allows precise numeric simulation of all flows within the mixer with excellent agreement to experiment^3-4.Laminar flow is achieved for Reynolds numbers Re ≤100. For aqueous solutions, this requires micron scale geometries. We use a hard substrate, such as silicon or fused silica, to make channels 5-10 μm wide and 10 μm deep (See Figure 1). The smallest dimensions, at the entrance to the mixing region, are on the order of 1 μm in size. The chip is sealed with a thin glass or fused silica coverslip for optical access. Typical total linear flow rates are ~1 m/s, yielding Re~10, but the protein consumption is only ~0.5 nL/s or 1.8 μL/hr. Protein concentration depends on the detection method: For tryptophan fluorescence the typical concentration is 100 μM (for 1 Trp/protein) and for FRET the typical concentration is ~100 nM.The folding process is initiated by rapid dilution of denaturant from 6 M to 0.06 M guanidine hydrochloride. The protein in high denaturant flows down a central channel and is met on either side at the mixing region by buffer without denaturant moving ~100 times faster (see Figure 2). This geometry causes rapid constriction of the protein flow into a narrow jet ~100 nm wide. Diffusion of the light denaturant molecules is very rapid, while diffusion of the heavy protein molecules is much slower, diffusing less than 1 μm in 1 ms. The difference in diffusion constant of the denaturant and the protein results in rapid dilution of the denaturant from the protein stream, reducing the effective concentration of the denaturant around the protein. The protein jet flows at a constant rate down the observation channel and fluorescence of the protein during folding can be observed using a scanning confocal microscope⁵.     

First-Principles Protein Folding Simulations

Abstract

It is widely believed that the prediction of the three-dimensional structures of proteins from the first principles is impossible. This view is based on the fact that the number of possible structures for each protein is astronomically large. The question is then why a protein folds into its native structure with the proper biological functions in the time scale of milliseconds to minutes, and this is called Levinthal's paradox. In this article I will discuss our strategy for attacking the protein folding problem. Our approach consists of two elements: the inclusion of accurate solvent effects and the development of powerful simulation algorithms that can avoid getting trapped in states of energy local minima. For the former, we discuss several models varying in nature from crude (distance-dependent dielectric function) to rigorous (reference interaction site model). For the latter, we show the effectiveness of Monte Carlo simulated annealing and generalized-ensemble algorithms.     

蛋白质折叠和分子伴侣

一个有活性的蛋白质分子不但有特定的氨基酸序列,还处于特定的由氨基酸序列决定的三维空间结构。三维结构的完整性受到干扰,生物活性也会发生变化：有时即使只是轻微的破坏,都可能导致其生物活性全部丧失。所以蛋白质的生物功能是与其三维空间结构密切联系在一起的。     

Toward Correct Protein Folding Potentials

Empirical protein folding potentialfunctions should have a global minimum nearthe native conformationof globular proteins that fold stably, andthey should give the correct free energy offolding. We demonstrate that otherwise verysuccessful potentials fail to have even alocal minimumanywhere near the native conformation, anda seemingly well validated method ofestimatingthe thermodynamic stability of the nativestate is extremely sensitive to smallperturbations inatomic coordinates. These are bothindicative of fitting a great deal ofirrelevant detail. Here weshow how to devise a robust potentialfunction that succeeds very well at bothtasks, at least for alimited set of proteins, and this involvesdeveloping a novel representation of thedenatured state.Predicted free energies of unfolding for 25mutants of barnase are in close agreementwith theexperimental values, while for 17 mutantsthere are substantial discrepancies.     

Back to Units of Protein Folding

Abstract

In response to the criticism by A. Finkelstein (J. Biomol. Struct. Dyn. 20, 311–314, 2002) of our Communication (J. Biomol. Struct. Dyn. 20, 5–6, 2002) several issues are dealt with. Importance of the notion of elementary folding unit, its size and structure, and the necessity of further characterization of the units for the elucidation of the protein folding in vivo are discussed. The criticism (J. Biomol. Struct. Dyn. 20, 311–314, 2002) on the hierarchical protein folding is also briefly addressed.     

50+ Years of Protein Folding

The ability of proteins to spontaneously form their spatial structures is a long-standing puzzle in molecular biology. Experimentally measured rates of spontaneous folding of single-domain globular proteins range from microseconds to hours: the difference–10-11 orders of magnitude–is the same as between the lifespan of a mosquito and the age of the Universe. This review (based on the literature and some personal recollections) describes a winding road to understanding spontaneous folding of protein structure. The main attention is given to the free-energy landscape of conformations of a protein chain–especially to the barrier separating its unfolded (U) and the natively folded (N) states–and to physical the-ories of rates of crossing this barrier in both directions: from U to N, and from N to U. It is shown that theories of both these processes come to essentially the same result and outline the observed range of folding and unfolding rates for single-domain globular proteins. In addition, they predict the maximal size of protein domains that fold under solely thermodynamic (rather than kinetic) control, and explain the observed maximal size of “foldable” protein domains.     

Experimental and Theoretical Protein Folding

Abstract

Resonance Raman spectra with Q-band excitation are reported for microperoxidase-11, the cytochrome c analog. Spectra were acquired in the mid-frequency range for the oxidized, and reduced forms of the undecapeptide, as well as for the imidazole and carbonyl complexes. Oxidation and spin state marker bands of the undecapeptides are consistent with a six-coordinate, low spin iron in both oxidation states. Porphyrin core size correlations yield a porphyrin-centre to pyrrole-nitrogen distance of 2.00 Å for MP11, suggestive of a six-coordinate species in a distorted heme environment. Molecular dynamics results show that the non-planarity of the heme of the parent cytochrome is conserved in the microperoxidase and its carbonmonoxy analog.     

模拟蛋白质折叠过程的新算法研究

蛋白质折叠过程模拟是当前蛋白质研究领域的一个难点问题。针对这一问题,提出了一个描述蛋白质折叠过程的算法-拟蛇算法,并且从分子振荡和分子动力学理论两个方面来证明该算法的核心函数是可行和正确的。经过实验总结出所有蛋白质空间结构都可以通过两种类型函数构造出来,提出了描述蛋白质折叠过程模型。与其它蛋白质折叠过程模拟算法的实验结果比较表明,拟蛇算法所构造的空间结构能量值最小、相似度最好。进而说明拟蛇算法和蛋白质折叠过程模型在描述蛋白质折叠过程方面具有明显优势。     

The Role of Solvent in Protein Folding and in Aggregation

We discuss features of the effect of solvent on protein folding andaggregation, highlighting the physics related to the particulate nature and the peculiar structure of the aqueous solvent, and the biological significance of interactions between solvent and proteins. To this purpose we use a generalized energy landscape of extended dimensionality. A closer look at the properties of solvent induced interactions and forces proves useful for understanding the physical grounds of `ad hoc' interactions and for devising realistic ways of accounting for solvent effects. The solvent has long been known to be a crucially important part of biological systems, and times appear mature for it to be adequately accounted for in the protein folding problem. Use of the extended dimensionality energy landscape helpseliciting the possibility of coupling among conformational changes and aggregation, such as proved by experimental data in the literature.     

Ultrafast Protein Folding in Membrane-Mimetic Environments

Proteins fold on timescales from hours to microseconds. In addition to protein size, sequence, and topology, the environment represents an equally important factor in determining folding speed. This is particularly relevant for proteins that require a lipid membrane or a membrane mimic to fold. However, only little is known about how properties of such a hydrophilic/hydrophobic interface modulate the folding landscape of membrane-interacting proteins. Here, we studied the influence of different membrane-mimetic micellar environments on the folding and unfolding kinetics of the helical-bundle protein Mistic. Devising a single-molecule fluorescence spectroscopy approach, we extracted folding and unfolding rates under equilibrium conditions and dissected the contributions from different detergent moieties to the free-energy landscape. While both polar and nonpolar moieties contribute to stability, they exert differential effects on the free-energy barrier: Hydrophobic burial stabilizes the folded state but not the transition state in reference to a purely aqueous environment; by contrast, zwitterionic headgroup moieties stabilize the folded state and, additionally, lower the free-energy barrier to accelerate the folding of Mistic to achieve ultrafast folding times down to 35 μs.     

Going the Distance for Protein Function Prediction: A New Distance Metric for Protein Interaction Networks

In protein-protein interaction (PPI) networks, functional similarity is often inferred based on the function of directly interacting proteins, or more generally, some notion of interaction network proximity among proteins in a local neighborhood. Prior methods typically measure proximity as the shortest-path distance in the network, but this has only a limited ability to capture fine-grained neighborhood distinctions, because most proteins are close to each other, and there are many ties in proximity. We introduce diffusion state distance (DSD), a new metric based on a graph diffusion property, designed to capture finer-grained distinctions in proximity for transfer of functional annotation in PPI networks. We present a tool that, when input a PPI network, will output the DSD distances between every pair of proteins. We show that replacing the shortest-path metric by DSD improves the performance of classical function prediction methods across the board.     

Multiplexed-Replica Exchange Molecular Dynamics Method for Protein Folding Simulation

Simulating protein folding thermodynamics starting purely from a protein sequence is a grand challenge of computational biology. Here, we present an algorithm to calculate a canonical distribution from molecular dynamics simulation of protein folding. This algorithm is based on the replica exchange method where the kinetic trapping problem is overcome by exchanging noninteracting replicas simulated at different temperatures. Our algorithm uses multiplexed-replicas with a number of independent molecular dynamics runs at each temperature. Exchanges of configurations between these multiplexed-replicas are also tried, rendering the algorithm applicable to large-scale distributed computing (i.e., highly heterogeneous parallel computers with processors having different computational power). We demonstrate the enhanced sampling of this algorithm by simulating the folding thermodynamics of a 23 amino acid miniprotein. We show that better convergence is achieved compared to constant temperature molecular dynamics simulation, with an efficient scaling to large number of computer processors. Indeed, this enhanced sampling results in (to our knowledge) the first example of a replica exchange algorithm that samples a folded structure starting from a completely unfolded state.     

蛋白质折叠速率数据集的构建及分析

近年来,随着高精度的蛋白质折叠速率实验数据的不断积累,使得从蛋白质折叠速率角度研究蛋白质折叠机制的理论工作者,迎来了前所未有的机遇和挑战。然而,却有约100多个蛋白质的折叠速率实验数据散落在2个数据库和若干文献中。为了方便今后的理论工作分析,作者将这些散落数据汇集整理出来,构建了一个包含109个非冗余单体野生型蛋白质的折叠速率数据集,称为PFRD109(protein folding rate dataset 109)。PFRD109所包含的109个蛋白质中,有69个二态蛋白和40个多态蛋白,折叠速率从10-4到106s-1,跨度为10个数量级。链长最短的为16 aa,最长为390 aa,二态蛋白平均长度为78 aa,多态蛋白平均长度为137 aa。当前,生物信息学对蛋白质折叠速率的研究,主要集中于寻找与折叠速率和折叠动力学相关的各种生化参数或拓扑参数,进而实现对蛋白质折叠速率和蛋白质折叠动力学类型的预测。因此,本文还针对PFRD109数据集,就这两个方面进行了一些参数的统计分析。     

Stoichiometry-driven Protein Folding: A Comment

Abstract

We have found, with the aid of 2-D gel electrophoresis, that double-stranded human telomeric repeat, (T₂AG₃)₁₂·(C₃TA₂)₁₂, being cloned within a plasmid, forms a protonated superhelically-induced structure. Experiments on chemical and enzymatic probing also indicate that the human telomeric repeats adopt an unusual structure. We have proposed an eclectic model for this structure in which four different elements coexist: a non-orthodox intramolecular triplex stabilized by the canonical protonated C · G*C⁺ base-triads and highly enriched by non-canonical base-triads; the intramolecular quadruplex formed by a portion of the G-rich strand; the single-stranded region encompassing a portion of the G-rich strand and, probably, the (C,A)-hairpin formed by a portion of the C-rich strand.     

Stoichiometry versus Hydrophobicity in Protein Folding

Abstract

In a channel-forming bundle of five alpha-helices of poly-L-alanine, the replacement of all the alanyl side-chains lining the inner wall by serines is shown, by energy optimization, to produce only small modifications of the packing. The stability of the bundle is larger than that of the pure alanyl package, owing to hydrogen bonding between serine hydroxyls and carbonyl oxygens. The energy profile for sodium as well as the water-channel interactions are favored by the presence of the OH groups and by the lability of the seryl side chains. The possible general significance of the results is suggested.