Kinetic role of early intermediates in protein folding

Burst Phase Analysis of the Aggregation Prone α-synuclein Amyloid Protein

While some studies inferred that valid information can be retrieved for the refolding of proteins and consequent identification of folding intermediates in the stopped-flow spectrometry collapse phase, other studies report that these burst phase folding intermediates can be questioned, implying a solvent-dependent modification of the still unfolded polypeptide chain. We therefore decided to investigate the burst phase occurring for the α-synuclein (Syn) amyloid protein by stopped-flow spectrometry. Solvent-dependent modification effects indeed occurred for the Nα-acetyl-L-tyrosinamide (NAYA) parent small compound and for the folded monomeric ubiquitin protein. More complex was the burst phase analysis of the disordered Syn amyloid protein. While this amyloid protein was determined to be aggregated at pH 7 and pH 2, in particular, this protein at pH 3 appears to be in a monomeric state in the burst phase analysis performed. In addition, the protein at pH 3 appears to suffer a hydrophobic collapse with the formation of a possible folded intermediate. This folded intermediate seems to result from a fast contraction of the disordered amyloid polypeptide chain, which is proceeded by an expansion of the protein, due to the occurrence of solvent-dependent modification effects in a milliseconds time scale of the burst phase. Generally, it can be argued that both literature criteria of solvent-dependent modifications of the disordered Syn amyloid protein and of the formation of its possible folded intermediate are very likely to occur in the burst phase.

Biosynthetic Protein Folding and Molecular Chaperons

The problem of linear polypeptide chain folding into a unique tertiary structure is one of the fundamental scientific challenges. The process of folding cannot be fully understood without its biological context, especially for big multidomain and multisubunit proteins. The principal features of biosynthetic folding are co-translational folding of growing nascent polypeptide chains and involvement of molecular chaperones in the process. The review summarizes available data on the early events of nascent chain folding, as well as on later advanced steps, including formation of elements of native structure. The relationship between the non-uniformity of translation rate and folding of the growing polypeptide is discussed. The results of studies on the effect of biosynthetic folding features on the parameters of folding as a physical process, its kinetics and mechanisms, are presented. Current understanding and hypotheses on the relationship of biosynthetic folding with the fundamental physical parameters and current views on polypeptide folding in the context of energy landscapes are discussed.

The change of conditions does not affect Ros87 downhill folding mechanism

Downhill folding has been defined as a unique thermodynamic process involving a conformations ensemble that progressively loses structure with the decrease of protein stability. Downhill folders are estimated to be rather rare in nature as they miss an energetically substantial folding barrier that can protect against aggregation and proteolysis. We have previously demonstrated that the prokaryotic zinc finger protein Ros87 shows a bipartite folding/unfolding process in which a metal binding intermediate converts to the native structure through a delicate barrier-less downhill transition. Significant variation in folding scenarios can be detected within protein families with high sequence identity and very similar folds and for the same sequence by varying conditions. For this reason, we here show, by means of DSC, CD and NMR, that also in different pH and ionic strength conditions Ros87 retains its partly downhill folding scenario demonstrating that, at least in metallo-proteins, the downhill mechanism can be found under a much wider range of conditions and coupled to other different transitions. We also show that mutations of Ros87 zinc coordination sphere produces a different folding scenario demonstrating that the organization of the metal ion core is determinant in the folding process of this family of proteins.

Folding of copper proteins: role of the metal?

Copper is a redox-active transition metal ion required for the function of many essential human proteins. For biosynthesis of proteins coordinating copper, the metal may bind before, during or after folding of the polypeptide. If the metal binds to unfolded or partially folded structures of the protein, such coordination may modulate the folding reaction. The molecular understanding of how copper is incorporated into proteins requires descriptions of chemical, thermodynamic, kinetic and structural parameters involved in the formation of protein-metal complexes. Because free copper ions are toxic, living systems have elaborate copper-transport systems that include particular proteins that facilitate efficient and specific delivery of copper ions to target proteins. Therefore, these pathways become an integral part of copper protein folding in vivo. This review summarizes biophysical-molecular in vitro work assessing the role of copper in folding and stability of copper-binding proteins as well as protein-protein copper exchange reactions between human copper transport proteins. We also describe some recent findings about the participation of copper ions and copper proteins in protein misfolding and aggregation reactions in vitro.

Ultrafast folding kinetics of WW domains reveal how the amino acid sequence determines the speed limit to protein folding

Protein (un)folding rates depend on the free-energy barrier separating the native and unfolded states and a prefactor term, which sets the timescale for crossing such barrier or folding speed limit. Because extricating these two factors is usually unfeasible, it has been common to assume a constant prefactor and assign all rate variability to the barrier. However, theory and simulations postulate a protein-specific prefactor that contains key mechanistic information. Here, we exploit the special properties of fast-folding proteins to experimentally resolve the folding rate prefactor and investigate how much it varies among structural homologs. We measure the ultrafast (un)folding kinetics of five natural WW domains using nanosecond laser-induced temperature jumps. All five WW domains fold in microseconds, but with a 10-fold difference between fastest and slowest. Interestingly, they all produce biphasic kinetics in which the slower phase corresponds to reequilibration over the small barrier (<3 RT) and the faster phase to the downhill relaxation of the minor population residing at the barrier top [transition state ensemble (TSE)]. The fast rate recapitulates the 10-fold range, demonstrating that the folding speed limit of even the simplest all-β fold strongly depends on the amino acid sequence. Given this fold's simplicity, the most plausible source for such prefactor differences is the presence of nonnative interactions that stabilize the TSE but need to break up before folding resumes. Our results confirm long-standing theoretical predictions and bring into focus the rate prefactor as an essential element for understanding the mechanisms of folding.

Ephemeral states in protein folding under force captured with a magnetic tweezers design

Magnetic tape heads are ubiquitously used to read and record on magnetic tapes in technologies as diverse as old VHS tapes, modern hard-drive disks, or magnetic bands on credit cards. Their design highlights the ability to convert electric signals into fluctuations of the magnetic field at very high frequencies, which is essential for the high-density storage demanded nowadays. Here, we twist this conventional use of tape heads to implement one in a magnetic tweezers design, which offers the unique capability of changing the force with a bandwidth of ∼10 kHz. We calibrate our instrument by developing an analytical expression that predicts the magnetic force acting on a superparamagnetic bead based on the Karlqvist approximation of the magnetic field created by a tape head. This theory is validated by measuring the force dependence of protein L unfolding/folding step sizes and the folding properties of the R3 talin domain. We demonstrate the potential of our instrument by carrying out millisecond-long quenches to capture the formation of the ephemeral molten globule state in protein L, which has never been observed before. Our instrument provides the capability of interrogating individual molecules under fast-changing forces with a control and resolution below a fraction of a piconewton, opening a range of force spectroscopy protocols to study protein dynamics under force.

Global computational mutagenesis of domain structures associated with inherited eye disease

Multidomain proteins account for 70% of the eukaryotic proteome. In genetic disease, multidomain proteins are often affected by numerous mutations, but the effects of these mutations on protein stability and their roles in genetic disease are not well understood. Here, we analyzed protein globular domains to understand how genetic mutations affect the stability of multidomain proteins in inherited disease. In total, 291 domain atomic structures from nine multidomain proteins were modeled by homology, equilibrated using molecular dynamics in water, and subjected to global computational mutagenesis. The domains were separated into 7 groups based on protein fold homology. Mutation propensities within each group of domains were then averaged to select residues critical for domain fold stability. The consensus derived from the sequence alignment shows that the critical residues determined by global mutagenesis are conserved within each group. From this analysis, we concluded that 80% of known disease-related genetic variants are associated with critical residues and are expected to have significant destabilizing effects on domain structure. Our work provides an in silico quantification of protein stability and could help to analyze the complex relationship among missense mutations, multidomain protein stability, and disease phenotypes in inherited eye disease.

Dihedral Angle Calculations To Elucidate the Folding of Peptides through Its Main Mechanical Forces

This study reports a general method to calculate dihedral angles (φ and ψ) of a given amino acid sequence, focusing on potential energy and torque moment concepts. By defining these physical measures in relation to the chemical interactions that occur on each single amino acid residue within a peptide, we analyze the folding process as the result of main mechanical forces (MMFs) exerted in the specific amino acid chain of interest. As a proof of concept, Leu-enkephalin was initially used as a model peptide to carry out the theoretical study. Our data show agreement between calculated Leu-enkephalin backbone dihedral angles and the corresponding experimentally determined X-ray values. Hence, we used calcitonin to validate our MMF-based method on a larger peptide, i.e., 32 amino acid residues forming an α-helix. Through a similar approach (although simplified with regard to electrostatic interactions), the calculations for calcitonin also demonstrate a good agreement with experimental values. This study offers new opportunities to analyze peptides' amino acid sequences and to help in the prediction of how they must fold, assisting in the development of new computational techniques in the field.

Initial Protein Unfolding Events in Ubiquitin, Cytochrome c and Myoglobin Are Revealed with the Use of 213 nm UVPD Coupled to IM-MS

The initial stages of protein unfolding may reflect the stability of the entire fold and can also reveal which parts of a protein can be perturbed, without restructuring the rest. In this work, we couple UVPD with activated ion mobility mass spectrometry to measure how three model proteins start to unfold. Ubiquitin, cytochrome c and myoglobin ions produced via nESI from salty solutions are subjected to UV irradiation pre-mobility separation; experiments are conducted with a range of source conditions which alter the conformation of the precursor ion as shown by the drift time profiles. For all three proteins, the compact structures result in less fragmentation than more extended structures which emerge following progressive in-source activation. Cleavage sites are found to differ between conformational ensembles, for example, for the dominant charge state of cytochrome c [M + 7H]⁷⁺, cleavage at Phe10, Thr19 and Val20 was only observed in activating conditions whilst cleavage at Ala43 is dramatically enhanced. Mapping the photo-cleaved fragments onto crystallographic structures provides insight into the local structural changes that occur as protein unfolding progresses, which is coupled to global restructuring observed in the drift time profiles. Graphical Abstract.

Complex Folding Landscape of Apomyoglobin at Acidic pH Revealed by Ultrafast Kinetic Analysis of Core Mutants

Under mildly acidic conditions (pH 4-4.5) apomyoglobin (apoMb) adopts a partially structured equilibrium state ( M-state) that structurally resembles a kinetic intermediate encountered at a late stage of folding to the native structure at neutral pH. We have previously reported that the M-state is formed rapidly (<1 ms) via a multistate process and thus offers a unique opportunity for exploring early stages of folding by both experimental and computational techniques. In order to gain structural insight into intermediates and barriers at the residue level, we studied the folding/unfolding kinetics of 12 apoMb mutants at pH 4.2 using fluorescence-detected ultrafast mixing techniques. Global analysis of the submillisecond folding/unfolding kinetics vs urea concentration for each variant, based on a sequential four-state mechanism ( U ⇔ I ⇔ L ⇔ M), allowed us to determine elementary rate constants and their dependence on urea concentration for most transitions. Comparison of the free energy diagrams constructed from the kinetic data of the mutants with that of wild-type apoMb yielded quantitative information on the effects of mutations on the free energy (ΔΔ G) of both intermediates and the first two kinetic barriers encountered during folding. Truncation of conserved aliphatic side chains on helices A, G, and H gives rise to a stepwise increase in ΔΔ G as the protein advances from U toward M, consistent with progressive stabilization of native-like contacts within the primary core of apoMb. Helix-helix contacts in the primary core contribute little to the first folding barrier ( U ⇔ I) and thus are not required for folding initiation but are critical for the stability of the late intermediate, L, and the M-state. Alanine substitution of hydrophobic residues at more peripheral helix-helix contact sites of the native structure, which are still absent or unstable in the M-state, shows both positive (destabilizing) and negative (stabilizing) ΔΔ G, indicating that non-native contacts are formed initially and weakened or lost as a result of subsequent structural rearrangement steps.

Co-translational protein targeting facilitates centrosomal recruitment of PCNT during centrosome maturation in vertebrates

As microtubule-organizing centers of animal cells, centrosomes guide the formation of the bipolar spindle that segregates chromosomes during mitosis. At mitosis onset, centrosomes maximize microtubule-organizing activity by rapidly expanding the pericentriolar material (PCM). This process is in part driven by the large PCM protein pericentrin (PCNT), as its level increases at the PCM and helps recruit additional PCM components. However, the mechanism underlying the timely centrosomal enrichment of PCNT remains unclear. Here, we show that PCNT is delivered co-translationally to centrosomes during early mitosis by cytoplasmic dynein, as evidenced by centrosomal enrichment of PCNT mRNA, its translation near centrosomes, and requirement of intact polysomes for PCNT mRNA localization. Additionally, the microtubule minus-end regulator, ASPM, is also targeted co-translationally to mitotic spindle poles. Together, these findings suggest that co-translational targeting of cytoplasmic proteins to specific subcellular destinations may be a generalized protein targeting mechanism.

Before a cell divides, it creates a copy of its genetic material (DNA) and evenly distributes it between the new ‘daughter’ cells with the help of a complex called the mitotic spindle. This complex is made of long cable-like protein chains called microtubules.

To ensure that each daughter cell receives an equal amount of DNA, structures known as centrosomes organize the microtubules during the division process. Centrosomes have two rigid cores, called centrioles, which are surrounded by a matrix of proteins called the pericentriolar material. It is from this material that the microtubules are organized.

The pericentriolar material is a dynamic structure and changes its size by assembling and disassembling its protein components. The larger the pericentriolar material, the more microtubules can form. Before a cell divides, it rapidly expands in a process called centrosome maturation. A protein called pericentrin initiates the maturation by helping to recruit other proteins to the centrosome. Pericentrin molecules are large, and it takes the cell between 10 and 20 minutes to make each one. Nevertheless, the cell can produce and deliver large quantities of pericentrin to the centrosome in a matter of minutes. We do not yet know how this happens.

To investigate this further, Sepulveda, Antkowiak, Brust-Mascher et al. used advanced microscopy to study zebrafish embryos and human cells grown in the laboratory. The results showed that cells build and transport pericentrin at the same time. Cells use messenger RNA molecules as templates to build proteins. These feed into protein factories called ribosomes, which assemble the building blocks in the correct order. Rather than waiting for the pericentrin production to finish, the cell moves the active factories to the centrosome with the help of a molecular motor called dynein. By the time the pericentrin molecules are completely made by ribosomes, they are already at the centrosome, ready to help with the recruitment of other proteins during centrosome maturation.

These findings improve our understanding of centrosome maturation. The next step is to find out how the cell coordinates this process with the recruitment of other proteins to the centrosome. It is also possible that the cell uses similar processes to deliver other proteins to different parts of the cell.

Study of protein folding under native conditions by rapidly switching the hydrostatic pressure inside an NMR sample cell

Development of specialized instrumentation enables rapid switching of the hydrostatic pressure inside an operating NMR spectrometer. This technology allows observation of protein signals during the repeated folding process. Applied to ubiquitin, a previously extensively studied model of protein folding, the methodology reveals an initially highly dynamic state that deviates relatively little from random coil behavior but also provides evidence for numerous repeatedly failed folding events, previously only observed in computer simulations. Above room temperature, direct NMR evidence shows a ∼50% fraction of proteins folding through an on-pathway kinetic intermediate, thereby revealing two equally efficient parallel folding pathways.

In general, small proteins rapidly fold on the timescale of milliseconds or less. For proteins with a substantial volume difference between the folded and unfolded states, their thermodynamic equilibrium can be altered by varying the hydrostatic pressure. Using a pressure-sensitized mutant of ubiquitin, we demonstrate that rapidly switching the pressure within an NMR sample cell enables study of the unfolded protein under native conditions and, vice versa, study of the native protein under denaturing conditions. This approach makes it possible to record 2D and 3D NMR spectra of the unfolded protein at atmospheric pressure, providing residue-specific information on the folding process. ¹⁵N and ¹³C chemical shifts measured immediately after dropping the pressure from 2.5 kbar (favoring unfolding) to 1 bar (native) are close to the random-coil chemical shifts observed for a large, disordered peptide fragment of the protein. However, ¹⁵N relaxation data show evidence for rapid exchange, on a ∼100-μs timescale, between the unfolded state and unstable, structured states that can be considered as failed folding events. The NMR data also provide direct evidence for parallel folding pathways, with approximately one-half of the protein molecules efficiently folding through an on-pathway kinetic intermediate, whereas the other half fold in a single step. At protein concentrations above ∼300 μM, oligomeric off-pathway intermediates compete with folding of the native state.

Protein folding and unfolding under force

The recent revolution in optics and instrumentation has enabled the study of protein folding using extremely low mechanical forces as the denaturant. This exciting development has led to the observation of the protein folding process at single molecule resolution and its response to mechanical force. Here, we describe the principles and experimental details of force spectroscopy on proteins, with a focus on the optical tweezers instrument. Several recent results will be discussed to highlight the importance of this technique in addressing a variety of questions in the protein folding field.

The effects of organic solvents on the folding pathway and associated thermodynamics of proteins: a microscopic view

Protein folding is subject to the effects of solvation environment. A variety of organic solvents are used as additives for in vitro refolding of denatured proteins. Examination of the solvent effects on protein folding could be of fundamental importance to understand the molecular interactions in determining protein structure. This article investigated the folding of α-helix and β-hairpin structures in water and the solutions of two representative refolding additives (methanol (MeOH) and 1-Ethyl-3-methylimidazolium chloride (EMIM-Cl) ionic liquid) using REMD simulations. For both α-helix and β-hairpin in MeOH/water solution or α-helix in EMIM-Cl/water solution, the transient structures along the folding pathway are consistent with the counterparts in water but the relative statistical weights are changed, leading to the decrease in the overall folding free energy barrier. Accordingly, MeOH promotes the folding of both α-helix and β-hairpin but EMIM-Cl ionic liquid only promotes the folding of α-helix, consistent with experimental observations. The present study reveals for the first time the trivial effects on folding route but significant effects on folding thermodynamics from MeOH and EMIM-Cl, explaining the function of protein refolding additives and testifying the validity of the folding mechanism revealed by in vitro protein folding study using refolding additives.

Accurate label-free reaction kinetics determination using initial rate heat measurements

Accurate label-free methods or assays to obtain the initial reaction rates have significant importance in fundamental studies of enzymes and in application-oriented high throughput screening of enzyme activity. Here we introduce a label-free approach for obtaining initial rates of enzyme activity from heat measurements, which we name initial rate calorimetry (IrCal). This approach is based on our new finding that the data recorded by isothermal titration calorimetry for the early stages of a reaction, which have been widely ignored, are correlated to the initial rates. Application of the IrCal approach to various enzymes led to accurate enzyme kinetics parameters as compared to spectroscopic methods and enabled enzyme kinetic studies with natural substrate, e.g. proteases with protein substrates. Because heat is a label-free property of almost all reactions, the IrCal approach holds promise in fundamental studies of various enzymes and in use of calorimetry for high throughput screening of enzyme activity.

Nonspecific shielding of unfavorable electrostatic intramolecular interactions in the erythropoietin native-state increase conformational stability and limit non-native aggregation

Previous equilibrium and kinetic folding studies of the glycoprotein erythropoietin indicate that sodium chloride increases the conformational stability of this therapeutically important cytokine, ostensibly by stabilizing the native-state [Banks DD, (2011) The Effect of Glycosylation on the Folding Kinetics of Erythropoietin. J Mol Biol 412:536-550]. The focus of the current report is to determine the underlying cause of the salt dependent increase in erythropoietin conformational stability and to understand if it has any impact on aggregation, an instability that remains a challenge to the biotech industry in maintaining the efficacy and shelf-life of protein therapeutics. Isothermal urea denaturation experiments conducted at numerous temperatures in the absence and presence of sodium chloride indicated that salt stabilizes erythropoietin primarily by increasing the difference in enthalpy between the native and unfolded sates. This result, and the finding that the salt induced increases in erythropoietin melting temperatures were independent of the identity of the salt cation and anion, indicates that salt likely increases the conformational stability of erythropoietin at neutral pH by nonspecific shielding of unfavorable electrostatic interaction(s) in the native-state. The addition of salt (even low concentrations of the strong chaotrope salt guanidinium hydrochloride) also exponentially decreased the initial rate of soluble erythropoietin non-native aggregation at 37 °C storage.

Single-molecule spectroscopy reveals chaperone-mediated expansion of substrate protein

Molecular chaperones are an essential part of the machinery that avoids protein aggregation and misfolding in vivo. However, understanding the molecular basis of how chaperones prevent such undesirable interactions requires the conformational changes within substrate proteins to be probed during chaperone action. Here we use single-molecule fluorescence spectroscopy to investigate how the DnaJ-DnaK chaperone system alters the conformational distribution of the denatured substrate protein rhodanese. We find that in a first step the ATP-independent binding of DnaJ to denatured rhodanese results in a compact denatured ensemble of the substrate protein. The following ATP-dependent binding of multiple DnaK molecules, however, leads to a surprisingly large expansion of denatured rhodanese. Molecular simulations indicate that hard-core repulsion between the multiple DnaK molecules provides the underlying mechanism for disrupting even strong interactions within the substrate protein and preparing it for processing by downstream chaperone systems.

Misfolding of luciferase at the single-molecule level

The folding of complex proteins can be dramatically affected by misfolding transitions. Directly observing misfolding and distinguishing it from aggregation is challenging. Experiments with optical tweezers revealed transitions between the folded states of a single protein in the absence of mechanical tension. Nonfolded chains of the multidomain protein luciferase folded within seconds to different partially folded states, one of which was stable over several minutes and was more resistant to forced unfolding than other partially folded states. Luciferase monomers can thus adopt a stable misfolded state and can do so without interacting with aggregation partners. This result supports the notion that luciferase misfolding is the cause of the low refolding yields and aggregation observed with this protein. This approach could be used to study misfolding transitions in other large proteins, as well as the factors that affect misfolding.

The role of key residues in structure, function, and stability of cytochrome-c

Cytochrome-c (cyt-c), a multi-functional protein, plays a significant role in the electron transport chain, and thus is indispensable in the energy-production process. Besides being an important component in apoptosis, it detoxifies reactive oxygen species. Two hundred and eighty-five complete amino acid sequences of cyt-c from different species are known. Sequence analysis suggests that the number of amino acid residues in most mitochondrial cyts-c is in the range 104 ± 10, and amino acid residues at only few positions are highly conserved throughout evolution. These highly conserved residues are Cys14, Cys17, His18, Gly29, Pro30, Gly41, Asn52, Trp59, Tyr67, Leu68, Pro71, Pro76, Thr78, Met80, and Phe82. These are also known as "key residues", which contribute significantly to the structure, function, folding, and stability of cyt-c. The three-dimensional structure of cyt-c from ten eukaryotic species have been determined using X-ray diffraction studies. Structure analysis suggests that the tertiary structure of cyt-c is almost preserved along the evolutionary scale. Furthermore, residues of N/C-terminal helices Gly6, Phe10, Leu94, and Tyr97 interact with each other in a specific manner, forming an evolutionary conserved interface. To understand the role of evolutionary conserved residues on structure, stability, and function, numerous studies have been performed in which these residues were substituted with different amino acids. In these studies, structure deals with the effect of mutation on secondary and tertiary structure measured by spectroscopic techniques; stability deals with the effect of mutation on T m (midpoint of heat denaturation), ∆G D (Gibbs free energy change on denaturation) and folding; and function deals with the effect of mutation on electron transport, apoptosis, cell growth, and protein expression. In this review, we have compiled all these studies at one place. This compilation will be useful to biochemists and biophysicists interested in understanding the importance of conservation of certain residues throughout the evolution in preserving the structure, function, and stability in proteins.

Folding dynamics of the Trp-cage miniprotein: evidence for a native-like intermediate from combined time-resolved vibrational spectroscopy and molecular dynamics simulations

Trp-cage is a synthetic 20-residue miniprotein which folds rapidly and spontaneously to a well-defined globular structure more typical of larger proteins. Due to its small size and fast folding, it is an ideal model system for experimental and theoretical investigations of protein folding mechanisms. However, Trp-cage's exact folding mechanism is still a matter of debate. Here we investigate Trp-cage's relaxation dynamics in the amide I' spectral region (1530-1700 cm(-1)) using time-resolved infrared spectroscopy. Residue-specific information was obtained by incorporating an isotopic label ((13)C═(18)O) into the amide carbonyl group of residue Gly11, thereby spectrally isolating an individual 310-helical residue. The folding-unfolding equilibrium is perturbed using a nanosecond temperature-jump (T-jump), and the subsequent re-equilibration is probed by observing the time-dependent vibrational response in the amide I' region. We observe bimodal relaxation kinetics with time constants of 100 ± 10 and 770 ± 40 ns at 322 K, suggesting that the folding involves an intermediate state, the character of which can be determined from the time- and frequency-resolved data. We find that the relaxation dynamics close to the melting temperature involve fast fluctuations in the polyproline II region, whereas the slower process can be attributed to conformational rearrangements due to the global (un)folding transition of the protein. Combined analysis of our T-jump data and molecular dynamics simulations indicates that the formation of a well-defined α-helix precedes the rapid formation of the hydrophobic cage structure, implying a native-like folding intermediate, that mainly differs from the folded conformation in the orientation of the C-terminal polyproline II helix relative to the N-terminal part of the backbone. We find that the main free-energy barrier is positioned between the folding intermediate and the unfolded state ensemble, and that it involves the formation of the α-helix, the 310-helix, and the Asp9-Arg16 salt bridge. Our results suggest that at low temperature (T ≪ Tm) a folding path via formation of α-helical contacts followed by hydrophobic clustering becomes more important.

Microsecond folding dynamics of apomyoglobin at acidic pH

Apomyolgobin (apoMb) is an important model for understanding the folding mechanism of helical proteins. This study focuses on a partially structured state of sperm whale apoMb populated at pH 4.2 (M-state), which structurally resembles a late kinetic intermediate in the formation of the native state (N) at higher pH. The thermodynamics and cooperativity of apoMb folding at pH 4.2 and 6.2 were studied by global analysis of the urea-induced unfolding transitions monitored by tryptophan fluorescence and circular dichroism. The kinetics of folding and unfolding of apoMb at pH 4.2 was measured over a time window from 40 to 850 μs, using fluorescence-detected continuous-flow measurements. Our observation of biphasic kinetics provides clear evidence for rapid (<100 μs) accumulation of previously unresolved intermediate states in both refolding and unfolding experiments. Quantitative kinetic modeling of the results, using a four-state mechanism with two intermediates on a direct route between the unfolded and folded states (U↔I↔L↔M), gave new insight into the conformational states and barriers that precede the rate-limiting step in the formation of the N-state of apoMb.

Overview of protein folding mechanisms: experimental and theoretical approaches to probing energy landscapes

We present an overview of the current experimental and theoretical approaches to studying protein folding mechanisms, set against current models of the folding energy landscape. We describe how stability and folding kinetics can be determined experimentally and how this data can be interpreted in terms of the characteristic features of various models from the simplest two-state pathway to a multi-state mechanism. We summarize the pros and cons of a range of spectroscopic methods for measuring folding rates and present a theoretical framework, coupled with protein engineering approaches, for elucidating folding mechanisms and structural features of folding transition states. A series of case studies are used to show how experimental kinetic data can be interpreted in the context of non-native interactions, populated intermediates, parallel folding pathways, and sequential transition states. We also show how computational methods now allow transient species of high energy, such as folding transition states, to be modeled on the basis of experimental Φ-value analysis derived from the effects of point mutations on folding kinetics.

Role of metal in folding and stability of copper proteins in vitro

Metal coordination is required for function of many proteins. For biosynthesis of proteins coordinating a metal, the question arises if the metal binds before, during or after folding of the polypeptide. Moreover, when the metal is bound to the protein, how does its coordination affect biophysical properties such as stability and dynamics? Understanding how metals are utilized by proteins in cells on a molecular level requires accurate descriptions of the thermodynamic and kinetic parameters involved in protein-metal complexes. Copper is one of the essential transition metals found in the active sites of many key proteins. To avoid toxicity of free copper ions, living systems have developed elaborate copper-transport systems that involve dedicated proteins that facilitate efficient and specific delivery of copper to target proteins. This review describes in vitro and in silico biophysical work assessing the role of copper in folding and stability of copper-binding proteins. Examples of proteins discussed are: a blue-copper protein (Pseudomonas aeruginosa azurin), members of copper-transport systems (bacterial CopZ, human Atox1 and ATP7B domains) and multi-copper ferroxidases (yeast Fet3p and human ceruloplasmin). The consequences of interactions between copper proteins and platinum-complexes are also discussed. This article is part of a Special Issue entitled: Cell Biology of Metals.

Microsecond unfolding kinetics of sheep prion protein reveals an intermediate that correlates with susceptibility to classical scrapie

The microsecond folding and unfolding kinetics of ovine prion proteins (ovPrP) were measured under various solution conditions. A fragment comprising residues 94-233 of the full-length ovPrP was studied for four variants with differing susceptibilities to classical scrapie in sheep. The observed biexponential unfolding kinetics of ovPrP provides evidence for an intermediate species. However, in contrast to previous results for human PrP, there is no evidence for an intermediate under refolding conditions. Global analysis of the kinetic data, based on a sequential three-state mechanism, quantitatively accounts for all folding and unfolding data as a function of denaturant concentration. The simulations predict that an intermediate accumulates under both folding and unfolding conditions, but is observable only in unfolding experiments because the intermediate is optically indistinguishable from the native state. The relative population of intermediates in two ovPrP variants, both transiently and under destabilizing equilibrium conditions, correlates with their propensities for classical scrapie. The variant susceptible to classical scrapie has a larger population of the intermediate state than the resistant variant. Thus, the susceptible variant should be favored to undergo the PrP(C) to PrP(Sc) conversion and oligomerization.

NMR analysis of a kinetically trapped intermediate of a disulfide-deficient mutant of the starch-binding domain of glucoamylase

A thermally unfolded disulfide-deficient mutant of the starch-binding domain of glucoamylase refolds into a kinetically trapped metastable intermediate when subjected to a rapid lowering of temperature. We attempted to characterise this intermediate using multidimensional NMR spectroscopy. The (1)H-(15)N heteronuclear single quantum coherence spectrum after a rapid temperature decrease (the spectrum of the intermediate) showed good chemical shift dispersion but was significantly different from that of the native state, suggesting that the intermediate adopts a nonnative but well-structured conformation. Large chemical shift changes for the backbone amide protons between the native and the intermediate states were observed for residues in the β-sheet consisting of strands 2, 3, 5, 6, and 7 as well as in the C-terminal region. These residues were found to be in close proximity to aromatic residues, suggesting that the chemical shift changes are mainly due to ring current shifts caused by the aromatic residues. The two-dimensional nuclear Overhauser enhancement (NOE) spectroscopy experiments showed that the intermediate contained substantial, native-like NOE connectivities, although there were fewer cross peaks in the spectrum of the intermediate compared with that of the native state. It was also shown that there were native-like interresidue NOEs for residues buried in the protein, whereas many of the NOE cross peaks were lost for the residues involved in a surface-exposed aromatic cluster. These results suggest that, in the intermediate, the aromatic cluster at the surface is structurally less organised, whereas the interior of the protein has relatively rigid, native-like side-chain packing.

The effect of glycosylation on the folding kinetics of erythropoietin

Glycosylation is a common posttranslational modification that generally increases protein solubility and thermodynamic stability. Less is known about how this modification influences protein folding, particularly folding processes involving intermediate species. In the present report, folding comparisons of a nonglycosylated erythropoietin (EPO) mutant are made with the fully glycosylated EPO, which was recently shown to fold by a three-state on-pathway mechanism. The absence of glycosylation did not alter the folding mechanism of EPO but did greatly decrease the stability of the intermediate species, change the rate-limiting step of the folding reaction, and accelerate the folding kinetics to both the intermediate state and the native state. Surprisingly, glycosylation stabilized the intermediate species to a greater extent than it increased the EPO equilibrium stability. These results suggest that glycosylation impedes the latter EPO folding steps rather than accelerating them by biasing particular folding pathways, as previously proposed for folding reactions initiated from unfolded ensembles with minimal residual structure. Due to the specific biological processes modulated by EPO glycosylation, however, there may be little evolutionary pressure to fold on a faster, more direct pathway at the expense of biological function, particularly given the protective role glycosylation has at preventing EPO aggregation. Lastly, evidence that is consistent with glycosylation destabilizing the unfolded state to some degree and contributing to the greater equilibrium stability of the glycosylated EPO is presented.

Probing the folding intermediate of Bacillus subtilis RNase P protein by nuclear magnetic resonance

Protein folding intermediates are often imperative for overall folding processes and consequent biological functions. However, the low population and transient nature of the intermediate states often hinder their biochemical and biophysical characterization. Previous studies have demonstrated that Bacillus subtilis ribonuclease P protein (P protein) is conformationally heterogeneous and folds with multiphasic kinetics, indicating the presence of an equilibrium and kinetic intermediate in its folding mechanism. In this study, nuclear magnetic resonance (NMR) spectroscopy was used to study the ensemble corresponding to this intermediate (I). The results indicate that the N-terminal and C-terminal helical regions are mostly unfolded in I. 1H−15N heteronuclear single-quantum coherence NMR spectra collected as a function of pH suggest that the protonation of His 22 may play a major role in the energetics of the equilibria among the unfolded, intermediate, and folded state ensembles of P protein. NMR paramagnetic relaxation enhancement experiments were also used to locate the small anion binding sites in both the intermediate and folded ensembles. The results for the folded protein are consistent with the previously modeled binding regions. These structural insights suggest a possible role for I in the RNase P holoenzyme assembly process.

Inhibition of unfolding and aggregation of lens protein human gamma D crystallin by sodium citrate

Cataract affects 1 in 6 Americans over the age of 40, and represents a global health problem. Mature onset cataract is associated with the aggregation of partially unfolded or damaged proteins in the lens, which accumulate as an individual ages. Currently, surgery is the primary effective treatment for cataract. As an alternative preventive approach, small molecules have been suggested as potential therapeutic agents. In this work, we study the effect of sodium citrate on the stability of Human γD Crystallin (HγD-Crys), a structural protein of the eye lens, and two cataract-related mutants, L5S HγD-Crys and I90F HγD-Crys. In equilibrium unfolding-refolding studies, the presence of 250 mM sodium citrate increased the transition midpoint of the N-terminal domain (N-td) of WT HγD-Crys and L5S HγD-Crys by 0.3 M GuHCl, the C-terminal domain (C-td) by 0.6 M GuHCl, and the single transition of I90F HγD-Crys by 0.4 M GuHCl. In kinetic unfolding reactions, sodium citrate stabilization effect was observed only for the mutant I90F HγD-Crys. In the presence of citrate, a kinetic unfolding intermediate of I90F HγD-Crys was observed, which was not populated in the absence of citrate. The rates of aggregation were measured using solution turbidity. Sodium citrate demonstrated negligible effect on rate of aggregation of WT HγD-Crys, but considerably slowed the rate of aggregation of both L5S HγD-Crys and I90F HγD-Crys. The presence of sodium citrate dramatically slowed refolding of WT HγD-Crys and I90F HγD-Crys, but had a significantly smaller effect on the refolding of L5S HγD-Crys. The differential stabilizing effect of sodium citrate suggests that the ion is binding to a partially unfolded conformation of the C-td, but a solution-based Hofmeister effect cannot be eliminated as a possible explanation for the effects observed. These results indicate that assessment of potential anti-cataract agents needs to include effects on the unfolding and aggregation pathways, as well as the native state.

Site-specific hydration dynamics in the nonpolar core of a molten globule by dynamic nuclear polarization of water

Water-protein interactions play a direct role in protein folding. The chain collapse that accompanies protein folding involves extrusion of water from the nonpolar core. For many proteins, including apomyoglobin (apoMb), hydrophobic interactions drive an initial collapse to an intermediate state before folding to the final structure. However, the debate continues as to whether the core of the collapsed intermediate state is hydrated and, if so, what the dynamic nature of this water is. A key challenge is that protein hydration dynamics is significantly heterogeneous, yet suitable experimental techniques for measuring hydration dynamics with site-specificity are lacking. Here, we introduce Overhauser dynamic nuclear polarization at 0.35 T via site-specific nitroxide spin labels as a unique tool to probe internal and surface protein hydration dynamics with site-specific resolution in the molten globular, native, and unfolded protein states. The (1)H NMR signal enhancement of water carries information about the local dynamics of the solvent within ∼10 Å of a spin label. EPR is used synergistically to gain insights on local polarity and mobility of the spin-labeled protein. Several buried and solvent-exposed sites of apoMb are examined, each bearing a covalently bound nitroxide spin label. We find that the nonpoloar core of the apoMb molten globule is hydrated with water bearing significant translational dynamics, only 4-6-fold slower than that of bulk water. The hydration dynamics of the native state is heterogeneous, while the acid-unfolded state bears fast-diffusing hydration water. This study provides a high-resolution glimpse at the folding-dependent nature of protein hydration dynamics.

The effects of pK(a) tuning on the thermodynamics and kinetics of folding: design of a solvent-shielded carboxylate pair at the a-position of a coiled-coil

The tuning of the pK(a) of ionizable residues plays a critical role in various protein functions, such as ligand-binding, catalysis, and allostery. Proteins harness the free energy of folding to position ionizable groups in highly specific environments that strongly affect their pK(a) values. To investigate the interplay among protein folding kinetics, thermodynamics, and pK(a) modulation, we introduced a pair of Asp residues at neighboring interior positions of a coiled-coil. A single Asp residue was replaced for an Asn side chain at the a-position of the coiled-coil from GCN4, which was also crosslinked at the C-terminus via a flexible disulfide bond. The thermodynamic and kinetic stability of the system was measured by circular dichroism and stopped-flow fluorescence as a function of pH and concentration of guanidine HCl. Both sets of data are consistent with a two-state equilibrium between fully folded and unfolded forms. Distinct pK(a) values of 6.3 and 5.35 are assigned to the first and second protonation of the Asp pair; together they represent an energetic difference of 5 kcal/mol relative to the protonation of two Asp residues with unperturbed pK(a) values. Analysis of the rate data as a function of pH and denaturant concentration allowed calculation of the kinetic constants for the conformational transitions of the peptide with the Asp residues in the doubly protonated, singly protonated, and unprotonated forms. The doubly and singly protonated forms fold rapidly, and a ϕ-value analysis shows that their contribution to folding occurs subsequent to the transition state ensemble for folding. By contrast, the doubly charged state shows a reduced rate of folding and a ϕ-value near 0.5 indicative of a repulsive interaction, and possibly also heterogeneity in the transition state ensemble.

Dynamics of protein folding: probing the kinetic network of folding-unfolding transitions with experiment and theory

The problem of spontaneous folding of amino acid chains into highly organized, biologically functional three-dimensional protein structures continues to challenge the modern science. Understanding how proteins fold requires characterization of the underlying energy landscapes as well as the dynamics of the polypeptide chains in all stages of the folding process. In recent years, important advances toward these goals have been achieved owing to the rapidly growing interdisciplinary interest and significant progress in both experimental techniques and theoretical methods. Improvements in the experimental time resolution led to determination of the timescales of the important elementary events in folding, such as formation of secondary structure and tertiary contacts. Sensitive single molecule methods made possible probing the distributions of the unfolded and folded states and following the folding reaction of individual protein molecules. Discovery of proteins that fold in microseconds opened the possibility of atomic-level theoretical simulations of folding and their direct comparisons with experimental data, as well as of direct experimental observation of the barrier-less folding transition. The ultra-fast folding also brought new questions, concerning the intrinsic limits of the folding rates and experimental signatures of barrier-less "downhill" folding. These problems will require novel approaches for even more detailed experimental investigations of the folding dynamics as well as for the analysis of the folding kinetic data. For theoretical simulations of folding, a main challenge is how to extract the relevant information from overwhelmingly detailed atomistic trajectories. New theoretical methods have been devised to allow a systematic approach towards a quantitative analysis of the kinetic network of folding-unfolding transitions between various configuration states of a protein, revealing the transition states and the associated folding pathways at multiple levels, from atomistic to coarse-grained representations. This article is part of a Special Issue entitled: Protein Dynamics: Experimental and Computational Approaches.

Theoretical perspectives on protein folding

Understanding how monomeric proteins fold under in vitro conditions is crucial to describing their functions in the cellular context. Significant advances in theory and experiments have resulted in a conceptual framework for describing the folding mechanisms of globular proteins. The sizes of proteins in the denatured and folded states, cooperativity of the folding transition, dispersions in the melting temperatures at the residue level, and timescales of folding are, to a large extent, determined by N, the number of residues. The intricate details of folding as a function of denaturant concentration can be predicted by using a novel coarse-grained molecular transfer model. By watching one molecule fold at a time, using single-molecule methods, investigators have established the validity of the theoretically anticipated heterogeneity in the folding routes and the N-dependent timescales for the three stages in the approach to the native state. Despite the successes of theory, of which only a few examples are documented here, we conclude that much remains to be done to solve the protein folding problem in the broadest sense.

Potentially amyloidogenic conformational intermediates populate the unfolding landscape of transthyretin: insights from molecular dynamics simulations

Protein aggregation into insoluble fibrillar structures known as amyloid characterizes several neurodegenerative diseases, including Alzheimer's, Huntington's and Creutzfeldt-Jakob. Transthyretin (TTR), a homotetrameric plasma protein, is known to be the causative agent of amyloid pathologies such as FAP (familial amyloid polyneuropathy), FAC (familial amyloid cardiomiopathy) and SSA (senile systemic amyloidosis). It is generally accepted that TTR tetramer dissociation and monomer partial unfolding precedes amyloid fibril formation. To explore the TTR unfolding landscape and to identify potential intermediate conformations with high tendency for amyloid formation, we have performed molecular dynamics unfolding simulations of WT-TTR and L55P-TTR, a highly amyloidogenic TTR variant. Our simulations in explicit water allow the identification of events that clearly discriminate the unfolding behavior of WT and L55P-TTR. Analysis of the simulation trajectories show that (i) the L55P monomers unfold earlier and to a larger extent than the WT; (ii) the single alpha-helix in the TTR monomer completely unfolds in most of the L55P simulations while remain folded in WT simulations; (iii) L55P forms, early in the simulations, aggregation-prone conformations characterized by full displacement of strands C and D from the main beta-sandwich core of the monomer; (iv) L55P shows, late in the simulations, severe loss of the H-bond network and consequent destabilization of the CBEF beta-sheet of the beta-sandwich; (v) WT forms aggregation-compatible conformations only late in the simulations and upon extensive unfolding of the monomer. These results clearly show that, in comparison with WT, L55P-TTR does present a much higher probability of forming transient conformations compatible with aggregation and amyloid formation.

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.

Measuring cotranslational folding of nascent polypeptide chains on ribosomes

Protein folding has been studied extensively in vitro, but much less is known about how folding proceeds in vivo. A particular distinction of folding in vivo is that folding begins while the nascent polypeptide chain is still undergoing synthesis by the ribosome. Studies of cotranslational protein folding are inherently much more complex than classical in vitro protein folding studies, and historically there have been few methods available to produce the quantities of pure material required for biophysical studies of the nascent chain, or assays to specifically interrogate its conformation. However, the past few years have produced dramatic methodological advances, which now place cotranslational folding studies within reach of more biochemists, enabling a detailed comparison of the earliest stages of protein folding on the ribosome to the wealth of information available for the refolding of full-length polypeptide chains in vitro.