Computational approaches to study protein unfolding: Hen egg white lysozyme as a case study

Four methods are compared to drive the unfolding of a protein: (1) high temperature (T-run), (2) high pressure (P-run), (3) by imposing a gradual increase in the mean radius of the protein using a penalty function added to the physical interaction function (F-run, radial force driven unfolding), and (4) by weak coupling of the difference between the temperature of the radially outward moving atoms and the radially inward moving atoms to an external temperature bath (K-run, kinetic energy driven unfolding). The characteristic features of the four unfolding pathways are analyzed in order to detect distortions due to the size or the type of the applied perturbation, as well as the features that are common to all of them. Hen egg white lysozyme is used as a test system. The simulations are analyzed and compared to experimental data like ¹H-NMR amide proton exchange-folding competition, heat capacity, and compressibility measurements. © 1995 Wiley-Liss, Inc.     

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.     

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.     

Tryptophan to Glycine mutation in the position 116 leads to protein aggregation and decreases the stability of the LITAF protein

Mutations in the gene-encoding vesicle lipopolysaccharide-induced tumor necrosis factor (LITAF) protein cause Charcot–Marie–Tooth type 1C (CMT1C) disease, a neurological disorder. The LITAF gene is mapped to chromosome number 16 and can be found at cytogenetic location 16p13 of the chromosome. CMT1C-linked small integral membrane protein of lysosome/late endosome mutants are loss-of-function mutants that act in a dominant negative manner to impair endosomal trafficking, leading to prolonged extracellular signal-regulated kinases 1/2 signaling downstream of ErbB activation. Mutation W116G in the LITAF decreases the stability of the protein and also interrupts the functioning of gene. We have analyzed the single nucleotide polymorphism (SNP) results of 28 nsSNPs obtained from dbSNP. We also carried out multiple molecular dynamics simulations of 200 ns and obtained results of root-mean-square deviation, root-mean-square fluctuation, radius of gyration, solvent-accessible surface area, H-bond, and principal component analysis to check and prove the stability of both the wild type and the mutant. The protein was then checked for its aggregation and the results showed loss of helix. The loss of helix leads to the instability of the protein.     

The neurodegenerative-disease-related protein sacsin is a molecular chaperone

Various human neurodegenerative disorders are associated with processes that involve misfolding of polypeptide chains. These so-called protein misfolding disorders include Alzheimer's and Parkinson's diseases and an increasing number of inherited syndromes that affect neurons involved in motor control circuits throughout the central nervous system. The reasons behind the particular susceptibility of neurons to misfolded proteins are currently not known. The main function of a class of proteins known as molecular chaperones is to prevent protein misfolding and aggregation. Although neuronal cells contain the major known classes of molecular chaperones, central-nervous-system-specific chaperones that maintain the neuronal proteome free from misfolded proteins are not well defined. In this study, we assign a novel molecular chaperone activity to the protein sacsin responsible for autosomal recessive spastic ataxia of Charlevoix-Saguenay, a degenerative disorder of the cerebellum and spinal cord. Using purified components, we demonstrate that a region of sacsin that contains a segment with homology to the molecular chaperone Hsp90 is able to enhance the refolding efficiency of the model client protein firefly luciferase. We show that this region of sacsin is highly capable of maintaining client polypeptides in soluble folding-competent states. Furthermore, we demonstrate that sacsin can efficiently cooperate with members of the Hsp70 chaperone family to increase the yields of correctly folded client proteins. Thus, we have identified a novel chaperone directly involved in a human neurodegenerative disorder.     

Minimalist protein model as a diagnostic tool for misfolding and aggregation

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

Protein misfolding, aggregation, and degradation in disease

Pathologies associated with protein misfolding have been observed in neurodegenerative diseases such as Alzheimer’s disease, metabolic diseases like phenylketonuria, and diseases affecting structural proteins like collagen or keratin. Misfolding of mutant proteins in these and many other diseases may result in premature degradation, formation of toxic aggregates, or incorporation of toxic conformations into structures. We review common traits of these diverse diseases under the unifying view of protein misfolding. The molecular pathogenesis is discussed in the context of protein quality control systems consisting of molecular chaperones and intracellular proteases that assist the folding and supervise the maintenance of the folded structure. Furthermore, genetic and environmental factors that may modify the severity of these diseases are underscored. The present article represents a partly revised and updated version of chapter 1 published earlier in volume 232 of the series Methods in Molecular Biology (Walker, J. M., ed., Humana Press, Totowa, NJ), Protein Misfolding and Disease: Principles and Protocols (Bross, P. & Gregersen, N., eds.), pp. 3–16 (2003).     

Probing the energy landscape of protein folding/unfolding transition states

Previous molecular dynamics (MD) simulations of the thermal denaturation of chymotrypsin inhibitor 2 (CI2) have provided atomic-resolution models of the transition state ensemble that is well supported by experimental studies. Here, we use simulations to further investigate the energy landscape around the transition state region. Nine structures within approximately 35 ps and 3 A C(alpha) RMSD of the transition state ensemble identified in a previous 498 K thermal denaturation simulation were quenched under the quasi-native conditions of 335 K and neutral pH. All of the structures underwent hydrophobically driven collapse in response to the drop in temperature. Structures less denatured than the transition state became structurally more native-like, while structures that were more denatured than the transition state tended to show additional loss of native structure. The structures in the immediate region of the transition state fluctuated between becoming more and less native-like. All of the starting structures had the same native-like topology and were quite similar (within 3.5 A C(alpha) RMSD). That the structures all shared native-like topology, yet diverged into either more or less native-like structures depending on which side of the transition state they occupied on the unfolding trajectory, indicates that topology alone does not dictate protein folding. Instead, our results suggest that a detailed interplay of packing interactions and interactions with water determine whether a partially denatured protein will become more native-like under refolding conditions.     

Manipulation of unfolding-induced protein aggregation by peptides selected for aggregate-binding ability through phage display library screening

A phage-displayed library of peptides (12-mer) was screened for the ability to bind to thermally aggregated bovine carbonic anhydrase (BCA), with a view toward examining whether peptides possessing this ability might bind to partially structured intermediates on the protein's unfolding pathway and, therefore, constitute useful tools for manipulation of the kinetic partitioning of molecules between the unfolded and aggregated states. Two peptides [N-HPSTMGLRTMHP-C and N-TPSAWKTALVKA-C] were identified and tested. While neither showed thermal aggregation autonomously, both peptides individually elicited remarkable increases in the levels of thermal aggregation of BCA. A possible explanation is that both peptides bind to surfaces on molten BCA that are not directly involved in aggregation. Such binding could slow down interconversions between folded and unfolded states and stabilize aggregation-prone intermediate(s) to make them more prone to aggregation, while failing to achieve any steric prevention of aggregation. The approach has the potential of yielding useful aggregation-aiding/inhibiting agents, and may provide clues to whether amorphous aggregates are "immobilized" forms of folding intermediates.     

Heat shock protein 90: The cancer chaperone

Heat shock protein 90 (Hsp90) is a molecular chaperone required for the stability and function of a number of conditionally activated and/or expressed signalling proteins, as well as multiple mutated, chimeric, and/or over-expressed signalling proteins, that promote cancer cell growth and/or survival. Hsp90 inhibitors are unique in that, although they are directed towards a specific molecular target, they simultaneously inhibit multiple cellular signalling pathways. By inhibiting nodal points in multiple overlapping survival pathways utilized by cancer cells, combination of an Hsp90 inhibitor with standard chemotherapeutic agents may dramatically increase the in vivo efficacy of the standard agent. Hsp90 inhibitors may circumvent the characteristic genetic plasticity that has allowed cancer cells to eventually evade the toxic effects of most molecularly targeted agents. The mechanism-based use of Hsp90 inhibitors, both alone and in combination with other drugs, should be effective toward multiple forms of cancer. Further, because Hsp90 inhibitors also induce Hsf-1-dependent expression of Hsp70, and because certain mutated Hsp90 client proteins are neurotoxic, these drugs display ameliorative properties in several neurodegenerative disease models, suggesting a novel role for Hsp90 inhibitors in treating multiple pathologies involving neurodegeneration.     

Accelerated simulation of unfolding and refolding of a large single chain globular protein

We have developed novel strategies for contracting simulation times in protein dynamics that enable us to study a complex protein with molecular weight in excess of 34 kDa. Starting from a crystal structure, we produce unfolded and then refolded states for the protein. We then compare these quantitatively using both established and new metrics for protein structure and quality checking. These include use of the programs Concoord and Darvols. Simulation of protein-folded structure well beyond the molten globule state and then recovery back to the folded state is itself new, and our results throw new light on the protein-folding process. We accomplish this using a novel cooling protocol developed for this work.     

Thermodynamic analysis of the aggregation propensity of oxidized Alzheimer's beta-amyloid variants

We have determined the critical concentrations of a set of 18 variants of Alzheimer's Abeta(1-40) peptide, each carrying a different residue at position 18. We find that the critical concentrations depend on the hydrophobicity and beta-sheet propensity of residue 18, and therefore on properties that we identified previously to affect also the kinetics by which these peptides aggregate. Since the critical concentrations can be related to the Gibbs free energy of aggregation (DeltaG), these data imply a link between the thermodynamics and the kinetics of aggregation in that sequences that form very stable aggregates are also those that form such aggregates very rapidly.     

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.     

Molecular mechanism of β-sheet self-organization at water-hydrophobic interfaces

The capacity to form β‐sheet structure and to self‐organize into amyloid aggregates is a property shared by many proteins. Severe neurodegenerative pathologies such as Alzheimer's disease are thought to involve the interaction of amyloidogenic protein oligomers with neuronal membranes. To understand the experimentally observed catalysis of amyloid formation by lipid membranes and other water‐hydrophobic interfaces, we examine the physico‐chemical basis of peptide adsorption and aggregation in a model membrane using atomistic molecular simulations. Blocked octapeptides with simple, repetitive sequences, (Gly‐Ala)₄, and (Gly‐Val)₄, are used as models of β‐sheet‐forming polypeptide chains found in the core of amyloid fibrils. In the presence of an n‐octane phase mimicking the core of lipid membranes, the peptides spontaneously partition at the octane‐water interface. The adsorption of nonpolar sidechains displaces the peptides' conformational equilibrium from a heterogeneous ensemble characterized by a high degree of structural disorder toward a more ordered ensemble favoring β‐hairpins and elongated β‐strands. At the interface, peptides spontaneously aggregate and rapidly evolve β‐sheet structure on a 10 to 100 ns time scale, while aqueous aggregates remain amorphous. Catalysis of β‐sheet formation results from the combination of the hydrophobic effect and of reduced conformational entropy of the polypeptide chain. While the former drives interfacial partition and displaces the conformational equilibrium of monomeric peptides, the planar interface further facilitates β‐sheet organization by increasing peptide concentration and reducing the dimensionality of self‐assembly from three to two. These findings suggest a general mechanism for the formation of β‐sheets on the surface of globular proteins and for amyloid self‐organization at hydrophobic interfaces. Proteins 2010. © 2010 Wiley‐Liss, Inc.     

Simulation and experiment conspire to reveal cryptic intermediates and a slide from the nucleation-condensation to framework mechanism of folding

There is a change from three-state to two-state kinetics of folding across the homeodomain superfamily of proteins as the mechanism slides from framework to nucleation-condensation. The tendency for framework folding in this family correlates with inherent helical propensity. The cellular myeloblastis protein (c-Myb) falls in the mechanistic transition region. An earlier, preliminary report of protein engineering experiments and molecular dynamics simulations (MD) showed that the folding mechanism for this protein has aspects of both the nucleation-condensation and framework models. In the more in-depth analysis of the MD trajectories presented here, we find that folding may be attributed to both of these mechanisms in different regions of the protein. The folding of the loop, middle helix, and turn is best described by nucleation-condensation, whereas folding of the N and C-terminal helices may be described by the framework model. Experimentally, c-Myb folds by apparent two-state kinetics, but the MD simulations predict that the kinetics hide a high-energy intermediate. We stabilized this hypothetical folding intermediate by deleting a residue (P174) in the loop between its second and third helices, and the mutant intermediate is long-lived in the simulations. Equilibrium and kinetic experiments demonstrate that folding of the DeltaP174 mutant is indeed three-state. The presence and shape of the intermediate observed in the simulations were confirmed by small angle X-ray scattering experiments.     

Enhanced sampling of peptide and protein conformations using replica exchange simulations with a peptide backbone biasing-potential

During replica exchange molecular dynamics (RexMD) simulations, several replicas of a system are simulated at different temperatures in parallel allowing for exchange between replicas at frequent intervals. This technique allows significantly improved sampling of conformational space and is increasingly being used for structure prediction of peptides and proteins. A drawback of the standard temperature RexMD is the rapid increase of the replica number with increasing system size to cover a desired temperature range. In an effort to limit the number of replicas, a new Hamiltonian-RexMD method has been developed that is specifically designed to enhance the sampling of peptide and protein conformations by applying various levels of a backbone biasing potential for each replica run. The biasing potential lowers the barrier for backbone dihedral transitions and promotes enhanced peptide backbone transitions along the replica coordinate. The application on several peptide cases including in all cases explicit solvent indicates significantly improved conformational sampling when compared with standard MD simulations. This was achieved with a very modest number of 5-7 replicas for each simulation system making it ideally suited for peptide and protein folding simulations as well as refinement of protein model structures in the presence of explicit solvent.     

Understanding the folding and stability of a designed WW domain protein with replica exchange molecular dynamics simulations

WW domain proteins are usually regarded as simple models for understanding the folding mechanism of β-sheet. CC45 is an artificial protein that is capable of folding into the same structure as WW domain. In this article, the replica exchange molecular dynamics simulations are performed to investigate the folding mechanism of CC45. The analysis of thermal stability shows that β-hairpin 1 is more stable than β-hairpin 2 during the unfolding process. Free energy analysis shows that the unfolding of this protein substantially proceeds through solvating the smaller β-hairpin 2, followed by the unfolding of β-hairpin 1. We further propose the unfolding process of CC45 and the folding mechanism of two β-hairpins. These results are similar to the previous folding studies of formin binding protein 28 (FBP28). Compared with FBP28, it is found that CC45 has more aromatic residues in N-terminal loop, and these residues contact with C-terminal loop to form the outer hydrophobic core, which increases the stability of CC45. Knowledge about the stability and folding behaviour of CC45 may help in understanding the folding mechanisms of the β-sheet and in designing new WW domains.     

Interactions between l‐arginine/l‐arginine derivatives and lysozyme and implications to their inhibition effects on protein aggregation

l ‐Arginine (Arg), l ‐homoarginine (HArg), l ‐arginine ethylester (ArgEE), and l ‐arginine methylester (ArgME) were found effective in inhibiting protein aggregation, but the molecular mechanisms are not clear. Herein, stopped‐flow fluorescence spectroscopy, isothermal titration calorimetry, and mass spectroscopy were used to investigate the folding kinetics of lysozyme and the interactions of the additives with lysozyme. It was found that the interactions of ArgME and ArgEE with lysozyme were similar to that of guanidine hydrochloride and were much stronger than those of Arg and HArg. The binding forces were all mainly hydrogen bonding and cation‐π interaction from the guanidinium group, but their differences in molecular states led to the significantly different binding strengths. The additives formed molecular clusters in an increasing order of ArgEE, ArgME, HArg, and Arg. Arg and HArg mainly formed annular clusters with head‐to‐tail bonding, while ArgME and ArgEE formed linear clusters with guanidinium groups stacking. The interactions between the additives and lysozyme were positively related to the monomer contents. That is, the monomers were the primary species that participated in the direct interactions due to their intact guanidinium groups and small sizes, while the clusters performed as barriers to crowd out the protein–protein interactions for aggregation. Thus, it is concluded that the effects of Arg and its derivatives on protein aggregation stemmed from the direct interactions by the monomers and the crowding effects by the clusters. Interplay of the two effects led to the differences in their inhibition effects on protein aggregation. © 2013 American Institute of Chemical Engineers Biotechnol. Prog., 29:1316–1324, 2013