Non-native interactions play an effective role in protein folding dynamics

Systematic Monte Carlo simulations of simple lattice models show that the final stage of protein folding is an ordered process where native contacts get locked (i.e., the residues come into contact and remain in contact for the duration of the folding process) in a well‐defined order. The detailed study of the folding dynamics of protein‐like sequences designed as to exhibit different contact energy distributions, as well as different degrees of sequence optimization (i.e., participation of non‐native interactions in the folding process), reveals significant differences in the corresponding locking scenarios—the collection of native contacts and their average locking times, which are largely ascribable to the dynamics of non‐native contacts. Furthermore, strong evidence for a positive role played by non‐native contacts at an early folding stage was also found. Interestingly, for topologically simple target structures, a positive interplay between native and non‐native contacts is observed also toward the end of the folding process, suggesting that non‐native contacts may indeed affect the overall folding process. For target models exhibiting clear two‐state kinetics, the relation between the nucleation mechanism of folding and the locking scenario is investigated. Our results suggest that the stabilization of the folding transition state can be achieved through the establishment of a very small network of native contacts that are the first to lock during the folding process.     

Insight into the correlation between lag time and aggregation rate in the kinetics of protein aggregation

Under favorable conditions, many proteins can assemble into macroscopically large aggregates such as the amyloid fibrils that are associated with Alzheimer's, Parkinson's, and other neurological and systemic diseases. The overall process of protein aggregation is characterized by initial lag time during which no detectable aggregation occurs in the solution and by maximal aggregation rate at which the dissolved protein converts into aggregates. In this study, the correlation between the lag time and the maximal rate of protein aggregation is analyzed. It is found that the product of these two quantities depends on a single numerical parameter, the kinetic index of the curve quantifying the time evolution of the fraction of protein aggregated. As this index depends relatively little on the conditions and/or system studied, our finding provides insight into why for many experiments the values of the product of the lag time and the maximal aggregation rate are often equal or quite close to each other. It is shown how the kinetic index is related to a basic kinetic parameter of a recently proposed theory of protein aggregation. Proteins 2010. © 2010 Wiley‐Liss, Inc.     

Weak protein interactions and pH- and temperature-dependent aggregation of human Fc1

The Fc (fragment crystallizable) is a common structural region in immunoglobulin gamma (IgG) proteins, IgG-based multi-specific platforms, and Fc-fusion platform technologies. Changes in conformational stability, protein-protein interactions, and aggregation of NS0-produced human Fc1 were quantified experimentally as a function of pH (4 to 6) and temperature (30 to 77°C), using a combination of differential scanning calorimetry, laser light scattering, size-exclusion chromatography, and capillary electrophoresis. The Fc1 was O-glycosylated at position 3 (threonine), and confirmed to correspond to the intact IgG1 by comparison with Fc1 produced by cleavage of the parent IgG1. Changing the pH caused large effects for thermal unfolding transitions, but it caused surprisingly smaller effects for electrostatic protein-protein interactions. The aggregation behavior was qualitatively similar across different solution conditions, with soluble dimers and larger oligomers formed in most cases. Aggregation rates spanned approximately 5 orders of magnitude and could be divided into 2 regimes: (i) Arrhenius, unfolding-limited aggregation at temperatures near or above the midpoint-unfolding temperature of the C_H2 domain; (ii) a non-Arrhenius regime at lower temperatures, presumably as a result of the temperature dependence of the unfolding enthalpy for the C_H2 domain. The non-Arrhenius regime was most pronounced for lower temperatures. Together with the weak protein-protein repulsions, these highlight challenges that are expected for maintaining long-term stability of biotechnology products that are based on human Fc constructs.     

Weak protein interactions and pH- and temperature-dependent aggregation of human Fc1

The Fc (fragment crystallizable) is a common structural region in immunoglobulin gamma (IgG) proteins, IgG-based multi-specific platforms, and Fc-fusion platform technologies. Changes in conformational stability, protein-protein interactions, and aggregation of NS0-produced human Fc1 were quantified experimentally as a function of pH (4 to 6) and temperature (30 to 77°C), using a combination of differential scanning calorimetry, laser light scattering, size-exclusion chromatography, and capillary electrophoresis. The Fc1 was O-glycosylated at position 3 (threonine), and confirmed to correspond to the intact IgG1 by comparison with Fc1 produced by cleavage of the parent IgG1. Changing the pH caused large effects for thermal unfolding transitions, but it caused surprisingly smaller effects for electrostatic protein-protein interactions. The aggregation behavior was qualitatively similar across different solution conditions, with soluble dimers and larger oligomers formed in most cases. Aggregation rates spanned approximately 5 orders of magnitude and could be divided into 2 regimes: (i) Arrhenius, unfolding-limited aggregation at temperatures near or above the midpoint-unfolding temperature of the C_H2 domain; (ii) a non-Arrhenius regime at lower temperatures, presumably as a result of the temperature dependence of the unfolding enthalpy for the C_H2 domain. The non-Arrhenius regime was most pronounced for lower temperatures. Together with the weak protein-protein repulsions, these highlight challenges that are expected for maintaining long-term stability of biotechnology products that are based on human Fc constructs.     

Kinetics of thermal aggregation of tobacco mosaic virus coat protein

The kinetics of thermal aggregation of coat protein (CP) of tobacco mosaic virus (TMV) have been studied at 42 and 52°C in a wide range of protein concentrations, [P]₀. The kinetics of aggregation were followed by monitoring the increase in the apparent absorbance (A) at 320 nm. At 52°C the kinetic curves may be approximated by the exponential law in the range of TMV CP concentrations from 0.02 to 0.30 mg/ml, the first order rate constant being linearly proportional to [P]₀ (50 mM phosphate buffer, pH 8.0). The analogous picture was observed at 42°C in the range of TMV CP concentrations from 0.01 to 0.04 mg/ml (100 mM phosphate buffer, pH 8.0). At higher TMV CP concentrations the time of half-conversion approaches a limiting value with increasing [P]₀ and at sufficiently high protein concentrations the kinetic curves fall on a common curve in the coordinates {A/A _lim; t} (t is time and A _lim is the limiting value of A at t ). According to a mechanism of aggregation of TMV CP proposed by the authors at rather low protein concentrations the rate of aggregation is limited by the stage of growth of aggregate, which proceeds as a reaction of the pseudo-first order, whereas at rather high protein concentrations the rate-limiting stage is the stage of protein molecule unfolding.     

Protein misfolding and aggregation research: Some thoughts on improving quality and utility

Once misfolded and aggregated proteins were as interesting as yesterday's trash, just a bothersome byproduct of productive activities. Today, they attract sustained interest from both basic researchers and practicing engineers. In the burgeoning biopharmaceutical industry, protein misfolding and aggregation pose significant challenges to the economic manufacture of safe and effective protein products. In the clinic, protein aggregates are believed to be pathological agents in a number of serious neurodegenerative disorders, such as Alzheimer's and Parkinson's. Over the past few years, the quantity of research into biotechnological aspects of protein misfolding and aggregation has skyrocketed. However, the quality of the published work is quite variable. In this brief opinion piece, we describe what we believe are some key features of high‐quality publications in protein aggregation. We focus on experimental studies that may also have a kinetic modeling component. © 2013 American Institute of Chemical Engineers Biotechnol. Prog., 29:1109–1115, 2013     

Change in kinetic regime of protein aggregation with temperature increase. Thermal aggregation of rabbit muscle creatine kinase

Creatine kinase thermal aggregation kinetics has been studied in 30 mM Hepes-NaOH buffer, pH 8.0, at two temperatures: 50.6 and 60°C. Aggregation kinetics was analyzed by measuring the growth of apparent absorption (A) at 400 nm. It was found that the limiting value of apparent absorption (A _lim) is proportional to protein concentration at both temperatures. The first order rate constant (k _I) does not depend on protein concentration in the range 0.05–0.2 mg/ml at temperature 50.6°C, but at temperature 60°C it increases with the growth of protein concentration in the range 0.1–0.4 mg/ml. Kinetic curves, shown in coordinates {A/A _lim; t}, in experiments at 50.6°C fuse to a common curve, which coincides with the theoretical curve of creatine kinase denaturation calculated using the denaturation rate constant determined from differential scanning calorimetry. At temperature 60°C, half-transformation time t _1/2 = ln2/k _I decreases when protein concentration grows. We conclude that when temperature increased from 50.6 to 60°C, change in the kinetic regime of thermal creatine kinase aggregation took place: at 50.6°C aggregation rate is limited by the stage of protein molecule denaturation, but at 60°C it is limited by the stage of protein aggregate growth, which proceeds as a reaction of pseudo-first order. Small heat shock protein Hsp 16.3 Mycobacterium tuberculosis suppresses the creatine kinase aggregation. Published in Russian in Biokhimiya, 2006, Vol. 71, No. 3, pp. 408–416.     

Interplay between protein homeostasis networks in protein aggregation and proteotoxicity

The misfolding and aggregation of disease proteins is characteristic of numerous neurodegenerative diseases. Particular neuronal populations are more vulnerable to proteotoxicity while others are more apt to tolerate the misfolding and aggregation of disease proteins. Thus, the cellular environment must play a significant role in determining whether disease proteins are converted into toxic or benign forms. The endomembrane network of eukaryotes divides the cell into different subcellular compartments that possess distinct sets of molecular chaperones and protein interaction networks. Chaperones act as agonists and antagonists of disease protein aggregation to prevent the accumulation of toxic intermediates in the aggregation pathway. Interacting partners can also modulate the conformation and localization of disease proteins and thereby influence proteotoxicity. Thus, interplay between these protein homeostasis network components can modulate the self‐association of disease proteins and determine whether they elicit a toxic or benign outcome. © 2009 Wiley Periodicals, Inc. Biopolymers 93: 229–236, 2010. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com     

Interpreting the aggregation kinetics of amyloid peptides

Amyloid fibrils are insoluble mainly beta-sheet aggregates of proteins or peptides. The multi-step process of amyloid aggregation is one of the major research topics in structural biology and biophysics because of its relevance in protein misfolding diseases like Alzheimer's, Parkinson's, Creutzfeld-Jacob's, and type II diabetes. Yet, the detailed mechanism of oligomer formation and the influence of protein stability on the aggregation kinetics are still matters of debate. Here a coarse-grained model of an amphipathic polypeptide, characterized by a free energy profile with distinct amyloid-competent (i.e. beta-prone) and amyloid-protected states, is used to investigate the kinetics of aggregation and the pathways of fibril formation. The simulation results suggest that by simply increasing the relative stability of the beta-prone state of the polypeptide, disordered aggregation changes into fibrillogenesis with the presence of oligomeric on-pathway intermediates, and finally without intermediates in the case of a very stable beta-prone state. The minimal-size aggregate able to form a fibril is generated by collisions of oligomers or monomers for polypeptides with unstable or stable beta-prone state, respectively. The simulation results provide a basis for understanding the wide range of amyloid-aggregation mechanisms observed in peptides and proteins. Moreover, they allow us to interpret at a molecular level the much faster kinetics of assembly of a recently discovered functional amyloid with respect to the very slow pathological aggregation.     

Time scale of protein aggregation dictated by liquid-liquid demixing

The growing impact of protein aggregation pathologies, together with the current high need for extensive information on protein structures are focusing much interest on the physics underlying the nucleation and growth of protein aggregates and crystals. Sickle Cell Hemoglobin (HbS), a point-mutant form of normal human Hemoglobin (HbA), is the first recognized and best-studied case of pathologically aggregating protein. Here we reanalyze kinetic data on nucleation of deoxy-HbS aggregates by referring them to the (concentration-dependent) temperature T(s) characterizing the occurrence of the phase transition of liquid-liquid demixing (LLD) of the solution. In this way, and by appropriate scaling of kinetic data at different concentrations, so as to normalize their spans, the apparently disparate sets of data are seen to fall on a master curve. Expressing the master curve vs. the parameter epsilon = (T - T(s)) / T(s), familiar from phase transition theory, allows eliciting the role of anomalously large concentration fluctuations associated with the LLD phase transition and also allows decoupling quantitatively the role of such fluctuations from that of microscopic, inter-protein interactions leading to nucleation. Referring to epsilon shows how in a narrow temperature span, that is at T - T(s), nucleation kinetics can undergo orders-of-magnitude changes, unexpected in terms of ordinary chemical kinetics. The same is true for similarly small changes of other parameters (pH, salts, precipitants), capable of altering T(s) and consequently epsilon. This offers the rationale for understanding how apparently minor changes of parameters can dramatically affect protein aggregation and related diseases.     

Effects of arginine on heat-induced aggregation of concentrated protein solutions

Arginine is one of the commonly used additives to enhance refolding yield of proteins, to suppress aggregation of proteins, and to increase solubility of proteins, and yet the molecular interactions that contribute to the role of arginine are unclear. Here, we present experiments, using bovine serum albumin (BSA), lysozyme (LYZ), and β-lactoglobulin (BLG) as model proteins, to show that arginine can enhance heat-induced aggregation of concentrated protein solutions, contrary to the conventional belief that arginine is a universal suppressor of aggregation. Results show that the enhancement in aggregation is caused only for BSA and BLG, but not for LYZ, indicating that arginine's preferential interactions with certain residues over others could determine the effect of the additive on aggregation. We use this previously unrecognized behavior of arginine, in combination with density functional theory calculations, to identify the molecular-level interactions of arginine with various residues that determine arginine's role as an enhancer or suppressor of aggregation of proteins. The experimental and computational results suggest that the guanidinium group of arginine promotes aggregation through the hydrogen-bond-based bridging interactions with the acidic residues of a protein, whereas the binding of the guanidinium group to aromatic residues (aggregation-prone) contributes to the stability and solubilization of the proteins. The approach, we describe here, can be used to select suitable additives to stabilize a protein solution at high concentrations based on an analysis of the amino acid content of the protein.     

Proline inhibits aggregation during protein refolding

The in vitro refolding of hen egg-white lysozyme is studied in the presence of various osmolytes. Proline is found to prevent aggregation during protein refolding. However, other osmolytes used in this study fail to exhibit a similar property. Experimental evidence suggests that proline inhibits protein aggregation by binding to folding intermediate(s) and trapping the folding intermediate(s) into enzymatically inactive, "aggregation-insensitive" state(s). However, elimination of proline from the refolded protein mixture results in significant recovery of the bacteriolytic activity. At higher concentrations (>1.5 M), proline is shown to form loose, higher-order molecular aggregate(s). The supramolecular assembly of proline is found to possess an amphipathic character. Formation of higher-order aggregates is believed to be crucial for proline to function as a protein folding aid. In addition to its role in osmoregulation under water stress conditions, the results of this study hint at the possibility of proline behaving as a protein folding chaperone.     

Computational approaches to understanding protein aggregation in neurodegeneration

The generation of toxic non-native protein conformers has emerged as a unifying thread among disorders such as Alzheimer＇s disease, Parkinson＇s disease, and amyotrophic lateral sclerosis. Atomic-level detail regarding dynamical changes that facilitate protein aggre- gation, as well as the structural features of large-scale ordered aggregates and soluble non-native oligomers, would contribute signifi- cantly to current understanding of these complex phenomena and offer potential strategies for inhibiting formation of cytotoxic species. However, experimental limitations often preclude the acquisition of high-resolution structural and mechanistic information for aggregating systems. Computational methods, particularly those combine both aU-atom and coarse-grained simulations to cover a wide range of time and length scales, have thus emerged as crucial tools for investigating protein aggregation. Here we review the current state of computational methodology for the study of protein self-assembly, with a focus on the application of these methods toward understanding of protein aggregates in human neurodegenerative disorders.     

A new kinetic scheme for lysozyme refolding and aggregation

The competing first- and third-order reaction scheme for lysozyme is shown to not predict fed-batch lysozyme refolding when the model is parameterized using independent batch experiments, even when variations in chemical composition during the fed-batch experiment are accounted for. A new kinetic scheme is proposed that involves rapid partitioning between the alternative fates of refolding and aggregation, and which allows for aggregation via a sequential mechanism. The model assumes that monomeric lysozyme in different states, including native, is able to aggregate with intermediates, accounting for recent experimental evidence that native protein can be incorporated into aggregates and explaining why native protein in the refolding buffer reduces yield. Stopped-flow light-scattering measurements were used to measure the association rate for the sequential aggregation mechanism, and refolding rate constants were determined in a series of batch experiments designed to be "snapshots" of the composition during a fed-batch experiment. The new kinetic scheme gave a good a priori prediction of fed-batch refolding performance.     

Aggregated proteins accelerate but do not increase the aggregation of D-glyceraldehyde-3-phosphate dehydrogenase. Specificity of protein aggregation

The effect of protein aggregates on the aggregation of d-glyceraldehyde-3-phosphate dehydrogenase (GAPDH) during unfolding and refolding has been studied. The aggregation of GAPDH follows a sigmoid course. The presence of protein aggregates increases the aggregation rate during unfolding and refolding of GAPDH but does not change the extent of aggregation and the final renaturation yield. It is suggested that protein aggregates function as seeds for aggregation via hydrophobic interaction with only GAPDH folding intermediates destined to aggregate and do not affect the distribution between pathways leading to correct folding and aggregation. Moreover, two different proteins do not interfere with each other during their simultaneous refolding together in a buffer. These findings provide insight into a mechanism by which cells prevent protein folding against the interference from aggregation of other proteins.     

Preventing serpin aggregation: the molecular mechanism of citrate action upon antitrypsin unfolding

The aggregation of antitrypsin into polymers is one of the causes of neonatal hepatitis, cirrhosis, and emphysema. A similar reaction resulting in disease can occur in other human serpins, and collectively they are known as the serpinopathies. One possible therapeutic strategy involves inhibiting the conformational changes involved in antitrypsin aggregation. The citrate ion has previously been shown to prevent antitrypsin aggregation and maintain the protein in an active conformation; its mechanism of action, however, is unknown. Here we demonstrate that the citrate ion prevents the initial misfolding of the native state to a polymerogenic intermediate in a concentration-dependent manner. Furthermore, we have solved the crystal structure of citrate bound to antitrypsin and show that a single citrate molecule binds in a pocket between the A and B beta-sheets, a region known to be important in maintaining antitrypsin stability.