Protein folding: Hypotheses and experiments

The current state of the problem of protein folding is reviewed with special attention to the novel molten globule state of the protein molecule, intermediate between the native and unfolded states. Experimental evidence on the existence of this state and its role in protein folding are compared with the sequential model of protein folding proposed by the author in 1972–1973.     

Protein folding following synthesis in vitro and in vivo: association of newly synthesized protein with 50S subunit of E. coli ribosome

In the accompanying paper, it was shown that a protein, while reverting to native form from the unfolded state in vitro with the help of bacterial 70S ribosome, split the latter into its subunits (50S and 30S) and remains associated with the 50S subunit. Here, we follow the fate of nascent proteins both in case of in vivo and in vitro translation system. The newly synthesised protein was found to associate with the 50S subunit in both the cases.     

How strong are side chain interactions in the folding intermediate?

Influence of 12 nonpolar amino acids residues from the hydrophobic core of apomyoglobin on stability of its native state and folding intermediate was studied. Six of the selected residues are from the A, G and H helices; these are conserved in structure of the globin family, although nonfunctional, that is, not involved in heme binding. The rest are nonconserved hydrophobic residues that belong to the B, C, D, and E helices. Each residue was substituted by alanine, and equilibrium pH‐induced transitions in apomyoglobin and its mutants were studied by circular dichroism and fluorescent spectroscopy. The obtained results allowed estimating changes in their free energy during formation of the intermediate state. It was first shown that the strength of side chain interactions in the apomyoglobin intermediate state amounts to 15–50% of that in its native state for conserved residues, and practically to 0% for nonconserved residues. These results allow a better understanding of interactions occurring in the intermediate state and shed light on involvement of certain residues in protein folding at different stages.     

Denatured states of yeast cytochrome c induced by heat and guanidinium chloride are structurally and thermodynamically different

A sequence alignment of mammalian cytochromes c with yeast iso-1-cytochrome c (y-cyt-c) shows that the yeast protein contains five extra N-terminal residues. We have been interested in understanding the question: What is the role of these five extra N-terminal residues in folding and stability of the protein? To answer this question we have prepared five deletants of y-cyt-c by sequentially removing these extra residues. During our studies on the wild type (WT) protein and its deletants, we observed that the amount of secondary structure in the guanidinium chloride (GdmCl)-induced denatured (D) state of each protein is different from that of the heat-induced denatured (H) state. This finding is confirmed by the observation of an additional cooperative transition curve of optical properties between H and D states on the addition of different concentrations of GdmCl to the already heat denatured WT y-cyt-c and its deletants at pH 6.0 and 68°C. For each protein, analysis of transition curves representing processes, native (N) state ? D state, N state ? H state, and H state ? D state, was done to obtain Gibbs free energy changes associated with all the three processes. This analysis showed that, for each protein, thermodynamic cycle accommodates Gibbs free energies associated with transitions between N and D states, N and H states, and H and D states, the characteristics required for a thermodynamic function. All these experimental observations have been supported by our molecular dynamics simulation studies.     

From Kuru to Alzheimer: A personal outlook

Seventy years ago, we learned from Chris Anfinsen that the stereochemical code necessary to fold a protein is embedded into its amino acid sequence. In water, protein morphogenesis is a spontaneous reversible process leading from an ensemble of disordered structures to the ordered functionally competent protein; conforming to Aristotle''s definition of substance, the synolon of matter and form. The overall process of folding is generally consistent with a two state transition between the native and the denatured protein: not only the denatured state is an ensemble of several structures, but also the native protein populates distinct functionally relevant conformational (sub)states. This two‐state view should be revised, given that any globular protein can populate a peculiar third state called amyloid, characterized by an overall architecture that at variance with the native state, is by‐and‐large independent of the primary structure. In a nut shell, we should accept that beside the folded and unfolded states, any protein can populate a third state called amyloid which gained center stage being the hallmark of incurable neurodegenerative disorders, such as Alzheimer''s and Parkinson''s diseases as well as others. These fatal diseases are characterized by clear‐cut clinical differences, yet display some commonalities such as the presence in the brain of amyloid deposits constituted by one misfolded protein specific for each disease. Some aspects of this complex problem are summarized here as an excursus from the prion''s fibrils observed in the brain of aborigines who died of Kuru to the amyloid detectable in the cortex of Alzheimer''s patients.     

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

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

Do proteins with similar folds have similar transition state structures? A diffuse transition state of the 169 residue apoflavodoxin

Apoflavodoxin from Anabaena PCC 7119 is a 169 residue globular protein of known structure and energetics. Here, we present a comprehensive Phi-value analysis to characterize the structure of its transition state. A total of 34 non-disruptive mutations are made throughout the structure and a range of Phi-values from zero to one are observed. In addition, a small set of eight aliphatic small-to-large mutations have been introduced in the hydrophobic core of the protein and they have been analyzed to investigate the feasibility of stabilizing the unfolding transition state by creating new non-native interactions. We find that the transition state of apoflavodoxin (so far the largest protein subjected to Phi-analysis) is diffuse and that it can be stabilized by unspecific hydrophobic interactions that can speed up the folding reaction. The data gathered on the apoflavodoxin transition state are compared with results from experimental studies in other proteins to revisit the relationship between the native state topology and transition state structure.     

Molecular dynamics simulations of photoactive yellow protein (PYP) in three states of its photocycle: a comparison with X-ray and NMR data and analysis of the effects of Glu46 deprotonation and mutation

Photoactive yellow protein (PYP) is a prototype of the PAS domain superfamily of signaling proteins. The signaling process is coupled to a three-state photocycle. After the photoinduced trans-cis isomerization of the chromophore, 4-hydroxycinnamic acid (pCA), an early intermediate (pR) is formed, which proceeds to a second intermediate state (pB) on a sub-millisecond time scale. The signaling process is thought to be connected to the conformational changes upon the formation of pB and its recovery to the ground state (pG), but the exact signaling mechanism is not known. Experimental studies of PYP by solution NMR and X-ray crystallography suggest a very flexible protein backbone in the ground as well as in the signaling state. The relaxation from the pR to the pB state is accompanied by the protonation of the chromophore's phenoxyl group. This was found to be of crucial importance for the relaxation process. With the goal of gaining a better understanding of these experimental observations on an atomistic level, we performed five MD simulations on the three different states of PYP: a 1 ns simulation of PYP in its ground state [pG(MD)], a 1 ns simulation of the pR state [pR(MD)], a 2 ns simulation of the pR state with the chromophore protonated (pRprot), a 2 ns simulation of the pR state with Glu46 exchanged by Gln (pRGln) and a 2 ns simulation of PYP in its signaling state [pB(MD)]. Comparison of the pG simulation results with X-ray and NMR data, and with the results obtained for the pB simulation, confirmed the experimental observations of a rather flexible protein backbone and conformational changes during the recovery of the pG from the pB state. The conformational changes in the region around the chromophore pocket in the pR state were found to be crucially dependent on the strength of the Glu46-pCA hydrogen bond, which restricts the mobility of the chromophore in its unprotonated form considerably. Both the mutation of Glu46 with Gln and the protonation of the chromophore weaken this hydrogen bond, leading to an increased mobility of pCA and large structural changes in its surroundings. These changes, however, differ considerably during the pRGln and pRprot simulations, providing an atomistic explanation for the enhancement of the rate constant in the Gln46 mutant. Electronic supplementary material to this article is available at and is accessible for athorized users. Electronic Publication     

Protein translocation through a tunnel induces changes in folding kinetics: a lattice model study

Compaction of a nascent polypeptide chain inside the ribosomal exit tunnel, before it leaves the ribosome, has been proposed to accelerate the folding of newly synthesized proteins following their release from the ribosome. Thus, we used Kinetic Monte Carlo simulations of a minimalist on-lattice model to explore the effect that polypeptide translocation through a variety of channels has on protein folding kinetics. Our results demonstrate that tunnel confinement promotes faster folding of a well-designed protein relative to its folding in free space by displacing the unfolded state towards more compact structures that are closer to the transition state. Since the tunnel only forbids rarely visited, extended configurations, it has little effect on a "poorly designed" protein whose unfolded state is largely composed of low-energy, compact, misfolded configurations. The beneficial effect of the tunnel depends on its width; for example, a too-narrow tunnel enforces unfolded states with limited or no access to the transition state, while a too-wide tunnel has no effect on the unfolded state entropy. We speculate that such effects are likely to play an important role in the folding of some proteins or protein domains in the cellular environment and might dictate whether a protein folds co-translationally or post-translationally.     

Equilibrium and kinetics of the unfolding and refolding of Escherichia coli Malate Synthase G monitored by circular dichroism and fluorescence spectroscopy

The equilibrium and kinetics studies of an 82 kDa large monomeric Escherichia coli protein Malate Synthase G (MSG) was investigated by far and near-UV CD, intrinsic tryptophan fluorescence and extrinsic fluorescence spectroscopy. We find that despite of its large size, folding is reversible, in vitro. Equilibrium unfolding process of MSG exhibited three-state transition thus, indicating the presence of at least a stable equilibrium intermediate. Thermodynamic parameters suggest this intermediate resembles the unfolded state. However, the equilibrium intermediate exhibits pronounced secondary structure as measured by far-UV CD, partial tertiary structure as delineated by near-UV CD, compactness (m value) and exposed hydrophobic surface area as assessed by ANS binding, typically depicting a molten globule state. The stopped-flow kinetic data provide clear evidence for the presence of a burst phase during the refolding pathway due to the formation of an early Intermediate, within the dead time of the instrument. Refolding from 4 M to various lower concentrations until 0.4 M of GdnHCl follow biphasic kinetics at lower concentrations of GdnHCl (<0.8 M), whereas monophasic kinetics at concentrations above 1.5 M. Also, rollover in the refolding and unfolding limbs of chevron plot verifies the presence of a fast kinetic intermediate at lower concentration of GdnHCl. Based upon the above observations we hereby propose the folding pathway of a large multi-domain protein Malate Synthase G.     

Molecular dynamics simulations reveal that apo‐HisJ can sample a closed conformation

The Escherichia coli histidine binding protein HisJ is a type II periplasmic binding protein (PBP) that preferentially binds histidine and interacts with its cytoplasmic membrane ABC transporter, HisQMP₂, to initiate histidine transport. HisJ is a bilobal protein where the larger Domain 1 is connected to the smaller Domain 2 via two linking strands. Type II PBPs are thought to undergo “Venus flytrap” movements where the protein is able to reversibly transition from an open to a closed conformation. To explore the accessibility of the closed conformation to the apo state of the protein, we performed a set of all‐atom molecular dynamics simulations of HisJ starting from four different conformations: apo‐open, apo‐closed, apo‐semiopen, and holo‐closed. The simulations reveal that the closed conformation is less dynamic than the open one. HisJ experienced closing motions and explored semiopen conformations that reverted to closed forms resembling those found in the holo‐closed state. Essential dynamics analysis of the simulations identified domain closing/opening and twisting as main motions. The formation of specific inter‐hinge strand and interdomain polar interactions contributed to the adoption of the closed apo‐conformations although they are up to 2.5‐fold less prevalent compared with the holo‐closed simulations. The overall sampling of the closed form by apo‐HisJ provides a rationale for the binding of unliganded PBPs with their cytoplasmic membrane ABC transporters. Proteins 2014; 82:386–398. © 2013 Wiley Periodicals, Inc.     

An NMR view of the unfolding process of rusticyanin: Structural elements that maintain the architecture of a beta-barrel metalloprotein

The unfolding process of the blue copper protein rusticyanin (Rc) as well as its dynamic and D(2)O/H(2)O exchange properties in an incipient unfolded state have been studied by heteronuclear NMR spectroscopy. Titrations of apo, Cu(I), and Cu(II)Rc with guanidinium chloride (GdmCl) show that the copper ion stabilizes the folded species and remains bound in the completely unfolded state. The oxidized state of the copper ion is more efficient than the reduced form in this respect. The long loop of Rc (where the first ligand of the copper ion is located) is one of the most mobile domains of the protein. This region has no defined secondary structure elements and is prone to exchange its amide protons. In contrast, the last loop (including a short alpha-helix) and the last beta-strand (where the other three ligands of the metal ion are located) form the most rigid domain of the protein. The results taken as a whole suggest that the first ligand detaches from the metal ion when the protein unfolds, while the other three ligands remain bound to it. The implications of these findings for the biological folding process of Rc are also discussed.     

Capturing the photo‐signaling state of a photoreceptor in a steady‐state fashion by binding a transition metal complex

Binding a small molecule to proteins causes conformational changes, but often to a limited extent. Here, we demonstrate that the interaction of a CO‐releasing molecule (CORM3) with a photoreceptor photoactive yellow protein (PYP) drives large structural changes in the latter. The interaction of CORM3 and a mutant of PYP, Met100Ala, not only trigger the isomerization of its chromophore, p‐coumaric acid, from its anionic trans configuration to a protonated cis configuration, but also increases the content of β‐sheet at the cost of α‐helix and random coil in the secondary structure of the protein. The CORM3 derived Met100Ala is found to highly resemble the signaling state, which is one of the key photo‐intermediates of this photoactive protein, in both protein local conformation and chromophore configuration. The organometallic reagents hold promise as protein engineering tools. This work highlights a novel approach to structurally accessing short lived intermediates of proteins in a steady‐state fashion.     

Easy method to predict solvent accessibility from multiple protein sequence alignments

An easy and uncomplicated method to predict the solvent accessibility state of a site in a multiple protein sequence alignment is described. The approach is based on amino acid exchange and compositional preference matrices for each of three accessibility states: buried, exposed, and intermediate. Calculations utilized a modified version of the 3D―ali databank, a collection of multiple sequence alignments anchored through protein tertiary structural superpositions. The technique achieves the same accuracy as much more complex methods and thus provides such advantages as computational affordability, facile updating, and easily understood residue substitution patterns useful to biochemists involved in protein engineering, design, and structural prediction. The program is available from the authors; and, due to its simplicity, the algorithm can be readily implemented on any system. For a given alignment site, a hand calculation can yield a comparative prediction. Proteins 32:190–199, 1998. © 1998 Wiley-Liss, Inc.     

Actin-stimulated myosin Mg2+-ATPase inhibition by brain protein

A low-molecular-weight protein, isolated from bovine brain, inhibits the actin-stimulated Mg-ATPase activity of myosin from striated muscle. This inhibition is probably related to its ability to cause actin to revert from its polymerized to its depolymerized state and to prevent the polymerization of actin. Its effect on the polymeric state of the actin has been demonstrated by viscosity studies. DNase inhibition assay, and electron microscopy. Heavy meromyosin can overcome the effect of the brain protein and stimulate the polymerization of actin. The inhibition of ATPase activity is in part dependent upon the relative amounts of brain protein, actin, and myosin. The apparent molecular weight of the brain protein is approximately 20,000 daltons. It appears to be a heat-labile glycoprotein containing 5-6% carbohydrate.     

核磁共振波谱法在蛋白质三维结构解析中的应用

蛋白质特定的三维结构与其生物功能密切相关,因此,研究蛋白质的三维结构有助于揭示其生物功能机制。将核磁共振（NMR）波谱法应用于研究溶液状态下蛋白质的三维结构,能够更加准确地揭示蛋白质结构与生物功能之间的关系。本文综述了NMR解析蛋白质三维结构的理论和技术方法,以及NMR结合其他生物物理手段,并辅以分子建模计算法研究蛋白质三维结构的研究进展和最新方法,为精准解析蛋白质的三维结构提供思路及策略。     

The Oligomerization of Phage T4 DNA-(Adenine-N6)-Methyltransferase and Its Effect on the Catalytic Characteristics of the Enzyme

The structural and catalytic properties of the phage T4 DNA-(adenine-N ⁶)-methyltransferase (EC 2.1.1.72) were studied at different enzyme–substrate concentration ratios by chemical crosslinking of the protein subunits and by measuring the presteady state kinetics of the reactions. Various structural states of the methyltransferase were correlated with its catalytic activity, and it was shown that the oligomeric forms of the enzyme are catalytically active but are characterized by the reaction parameters different from those of the monomer.     

Crystal structures reveal transient PERK luminal domain tetramerization in endoplasmic reticulum stress signaling

Stress caused by accumulation of misfolded proteins within the endoplasmic reticulum (ER) elicits a cellular unfolded protein response (UPR) aimed at maintaining protein‐folding capacity. PERK, a key upstream component, recognizes ER stress via its luminal sensor/transducer domain, but the molecular events that lead to UPR activation remain unclear. Here, we describe the crystal structures of mammalian PERK luminal domains captured in dimeric state as well as in a novel tetrameric state. Small angle X‐ray scattering analysis (SAXS) supports the existence of both crystal structures also in solution. The salient feature of the tetramer interface, a helix swapped between dimers, implies transient association. Moreover, interface mutations that disrupt tetramer formation in vitro reduce phosphorylation of PERK and its target eIF2α in cells. These results suggest that transient conversion from dimeric to tetrameric state may be a key regulatory step in UPR activation.     

Intrinsically disordered proteins: modes of binding with emphasis on disordered domains

Our notions of protein function have long been determined by the protein structure–function paradigm. However, the idea that protein function is dictated by a prerequisite complementarity of shapes at the binding interface is becoming increasingly challenged. Interactions involving intrinsically disordered proteins (IDPs) have indicated a significant degree of disorder present in the bound state, ranging from static disorder to complete disorder, termed ‘random fuzziness’. This review assesses the anatomy of an IDP and relates how its intrinsic properties permit promiscuity and allow for the various modes of interaction. Furthermore, a mechanistic overview of the types of disordered domains is detailed, while also relating to a recent example and the kinetic and thermodynamic principles governing its formation.     

Prediction of titration properties of structures of a protein derived from molecular dynamics trajectories.

This paper explores the dependence of the molecular dynamics (MD) trajectory of a protein molecule on the titration state assigned to the molecule. Four 100-ps MD trajectories of bovine pancreatic trypsin inhibitor (BPTI) were generated, starting from two different structures, each of which was held in two different charge states. The two starting structures were the X-ray crystal structure and one of the solution structures determined by NMR, and the charge states differed only in the ionization state of N terminus. Although it is evident that the MD simulations were too short to sample fully the equilibrium distribution of structures in each case, standard Poisson-Boltzmann titration state analysis of the resulting configurations shows general agreement between the overall titration behavior of the protein and the charge state assumed during MD simulation: at pH 7, the total net charge of the protein resulting from the titration analysis is consistently lower for the protein with the N terminus assumed to be neutral than for the protein with the N terminus assumed to be charged. For most of the ionizable residues, the differences in the calculated pKaS among the four trajectories are statistically negligible and remain in good agreement with the data obtained by crystal structure titration and by experiment. The exceptions include the N terminus, which responds directly to the change of its imposed charge; the C terminus, which in the NMR structure interacts strongly with the former; and a few other residues (Arg 1, Glu 7, Tyr 35, and Arg 42) whose pKaS reflect the initial structure and the limited trajectory lengths. This study illustrates the importance of the careful assignment of protonation states at the start of MD simulations and points to the need for simulation methods that allow for the variation of the protonation state in the calculation of equilibrium properties.