Reduced Levels of Misfolded and Aggregated Mutant p53 by Proteostatic Activation.

In malignant cancer, excessive amounts of mutant p53 often lead to its aggregation, a feature that was recently identified as druggable. Here, we describe that induction of a heat shock-related stress response mediated by Foldlin, a small-molecule tool compound, reduces the protein levels of misfolded/aggregated mutant p53, while contact mutants or wild-type p53 remain largely unaffected. Foldlin also prevented the formation of stress-induced p53 nuclear inclusion bodies. Despite our inability to identify a specific molecular target, Foldlin also reduced protein levels of aggregating SOD1 variants. Finally, by screening a library of 778 FDA-approved compounds for their ability to reduce misfolded mutant p53, we identified the proteasome inhibitor Bortezomib with similar cellular effects as Foldlin. Overall, the induction of a cellular heat shock response seems to be an effective strategy to deal with pathological protein aggregation. It remains to be seen however, how this strategy can be translated to a clinical setting.

Thermodynamic and Evolutionary Coupling between the Native and Amyloid State of Globular Proteins.

The amyloid-like aggregation propensity present in most globular proteins is generally considered to be a secondary side effect resulting from the requirements of protein stability. Here, we demonstrate, however, that mutations in the globular and amyloid state are thermodynamically correlated rather than simply associated. In addition, we show that the standard genetic code couples this structural correlation into a tight evolutionary relationship. We illustrate the extent of this evolutionary entanglement of amyloid propensity and globular protein stability. Suppressing a 600-Ma-conserved amyloidogenic segment in the p53 core domain fold is structurally feasible but requires 7-bp substitutions to concomitantly introduce two aggregation-suppressing and three stabilizing amino acid mutations. We speculate that, rather than being a corollary of protein evolution, it is equally plausible that positive selection for amyloid structure could have been a driver for the emergence of globular protein structure.

Assessing computational predictions of the phenotypic effect of cystathionine-beta-synthase variants.

Accurate prediction of the impact of genomic variation on phenotype is a major goal of computational biology and an important contributor to personalized medicine. Computational predictions can lead to a better understanding of the mechanisms underlying genetic diseases, including cancer, but their adoption requires thorough and unbiased assessment. Cystathionine-beta-synthase (CBS) is an enzyme that catalyzes the first step of the transsulfuration pathway, from homocysteine to cystathionine, and in which variations are associated with human hyperhomocysteinemia and homocystinuria. We have created a computational challenge under the CAGI framework to evaluate how well different methods can predict the phenotypic effect(s) of CBS single amino acid substitutions using a blinded experimental data set. CAGI participants were asked to predict yeast growth based on the identity of the mutations. The performance of the methods was evaluated using several metrics. The CBS challenge highlighted the difficulty of predicting the phenotype of an ex vivo system in a model organism when classification models were trained on human disease data. We also discuss the variations in difficulty of prediction for known benign and deleterious variants, as well as identify methodological and experimental constraints with lessons to be learned for future challenges.

Exposure of a cryptic Hsp70 binding site determines the cytotoxicity of the ALS-associated SOD1-mutant A4V.

The accumulation of toxic protein aggregates is thought to play a key role in a range of degenerative pathologies, but it remains unclear why aggregation of polypeptides into non-native assemblies is toxic and why cellular clearance pathways offer ineffective protection. We here study the A4V mutant of SOD1, which forms toxic aggregates in motor neurons of patients with familial amyotrophic lateral sclerosis (ALS). A comparison of the location of aggregation prone regions (APRs) and Hsp70 binding sites in the denatured state of SOD1 reveals that ALS-associated mutations promote exposure of the APRs more than the strongest Hsc/Hsp70 binding site that we could detect. Mutations designed to increase the exposure of this Hsp70 interaction site in the denatured state promote aggregation but also display an increased interaction with Hsp70 chaperones. Depending on the cell type, in vitro this resulted in cellular inclusion body formation or increased clearance, accompanied with a suppression of cytotoxicity. The latter was also observed in a zebrafish model in vivo. Our results suggest that the uncontrolled accumulation of toxic SOD1A4V aggregates results from insufficient detection by the cellular surveillance network.

SolubiS: Optimizing Protein Solubility by Minimal Point Mutations.

Protein solubility is adapted to endogeneous protein abundance in the cell where protein folding is also assisted by multiple chaperones. During recombinant protein production, purification and storage proteins are frequently handled at concentrations that are several orders of magnitude above their physiological concentration, often resulting in protein aggregation. Here we describe SolubiS, a method allowing for (1) detection of aggregation prone linear segments within a protein sequence and (2) identification of mutations that abolish the aggregation propensity of these segments without affecting the thermodynamic stability of the protein. Provided the availability of structural information this method is applicable to all globular proteins including antibodies, resulting both in increased in vitro protein solubility and in better protein production yields.

Prediction and Reduction of the Aggregation of Monoclonal Antibodies.

Protein aggregation remains a major area of focus in the production of monoclonal antibodies. Improving the intrinsic properties of antibodies can improve manufacturability, attrition rates, safety, formulation, titers, immunogenicity, and solubility. Here, we explore the potential of predicting and reducing the aggregation propensity of monoclonal antibodies, based on the identification of aggregation-prone regions and their contribution to the thermodynamic stability of the protein. Although aggregation-prone regions are thought to occur in the antigen binding region to drive hydrophobic binding with antigen, we were able to rationally design variants that display a marked decrease in aggregation propensity while retaining antigen binding through the introduction of artificial aggregation gatekeeper residues. The reduction in aggregation propensity was accompanied by an increase in expression titer, showing that reducing protein aggregation is beneficial throughout the development process. The data presented show that this approach can significantly reduce liabilities in novel therapeutic antibodies and proteins, leading to a more efficient path to clinical studies.

Solubis: a webserver to reduce protein aggregation through mutation.

Protein aggregation is a major factor limiting the biotechnological and therapeutic application of many proteins, including enzymes and monoclonal antibodies. The molecular principles underlying aggregation are by now sufficiently understood to allow rational redesign of natural polypeptide sequences for decreased aggregation tendency, and hence potentially increased expression and solubility. Given that aggregation-prone regions (APRs) tend to contribute to the stability of the hydrophobic core or to functional sites of the protein, mutations in these regions have to be carefully selected in order not to disrupt protein structure or function. Therefore, we here provide access to an automated pipeline to identify mutations that reduce protein aggregation by reducing the intrinsic aggregation propensity of the sequence (using the TANGO algorithm), while taking care not to disrupt the thermodynamic stability of the native structure (using the empirical force-field FoldX). Moreover, by providing a plot of the intrinsic aggregation propensity score of APRs corrected by the local stability of that region in the folded structure, we allow users to prioritize those regions in the protein that are most in need of improvement through protein engineering. The method can be accessed at http://solubis.switchlab.org/.

Structural hot spots for the solubility of globular proteins.

Natural selection shapes protein solubility to physiological requirements and recombinant applications that require higher protein concentrations are often problematic. This raises the question whether the solubility of natural protein sequences can be improved. We here show an anti-correlation between the number of aggregation prone regions (APRs) in a protein sequence and its solubility, suggesting that mutational suppression of APRs provides a simple strategy to increase protein solubility. We show that mutations at specific positions within a protein structure can act as APR suppressors without affecting protein stability. These hot spots for protein solubility are both structure and sequence dependent but can be computationally predicted. We demonstrate this by reducing the aggregation of human α-galactosidase and protective antigen of Bacillus anthracis through mutation. Our results indicate that many proteins possess hot spots allowing to adapt protein solubility independently of structure and function.

Solubis: optimize your protein.

MOTIVATION: Protein aggregation is associated with a number of protein misfolding diseases and is a major concern for therapeutic proteins. Aggregation is caused by the presence of aggregation-prone regions (APRs) in the amino acid sequence of the protein. The lower the aggregation propensity of APRs and the better they are protected by native interactions within the folded structure of the protein, the more aggregation is prevented. Therefore, both the local thermodynamic stability of APRs in the native structure and their intrinsic aggregation propensity are a key parameter that needs to be optimized to prevent protein aggregation. RESULTS: The Solubis method presented here automates the process of carefully selecting point mutations that minimize the intrinsic aggregation propensity while improving local protein stability.

WALTZ-DB: a benchmark database of amyloidogenic hexapeptides.

Accurate prediction of amyloid-forming amino acid sequences remains an important challenge. We here present an online database that provides open access to the largest set of experimentally characterized amyloid forming hexapeptides. To this end, we expanded our previous set of 280 hexapeptides used to develop the Waltz algorithm with 89 peptides from literature review and by systematic experimental characterisation of the aggregation of 720 hexapeptides by transmission electron microscopy, dye binding and Fourier transform infrared spectroscopy. This brings the total number of experimentally characterized hexapeptides in the WALTZ-DB database to 1089, of which 244 are annotated as positive for amyloid formation.

Selectivity of aggregation-determining interactions.

Protein aggregation is sequence specific, favoring self-assembly over cross-seeding with non-homologous sequences. Still, as the majority of proteins in a proteome are aggregation prone, the high level of homogeneity of protein inclusions in vivo both during recombinant overexpression and in disease remains surprising. To investigate the selectivity of protein aggregation in a proteomic context, we here compared the selectivity of aggregation-determined interactions with antibody binding. To that purpose, we synthesized biotin-labeled peptides, corresponding to aggregation-determining sequences of the bacterial protein ß-galactosidase and two human disease biomarkers: C-reactive protein and prostate-specific antigen. We analyzed the selectivity of their interactions in Escherichia coli lysate, human serum and human seminal plasma, respectively, using a Western blot-like approach in which the aggregating peptides replace the conventional antibody. We observed specific peptide accumulation in the same bands detected by antibody staining. Combined spectroscopic and mutagenic studies confirmed accumulation resulted from binding of the peptide on the identical sequence of the immobilized target protein. Further, we analyzed the sequence redundancy of aggregating sequences and found that about 90% of them are unique within their proteome. As a result, the combined specificity and low sequence redundancy of aggregating sequences therefore contribute to the observed homogeneity of protein aggregation in vivo. This suggests that these intrinsic proteomic properties naturally compartmentalize aggregation events in sequence space. In the event of physiological stress, this might benefit the ability of cells to respond to proteostatic stress by allowing chaperones to focus on specific aggregation events rather than having to face systemic proteostatic failure.

A genome-wide sequence-structure analysis suggests aggregation gatekeepers constitute an evolutionary constrained functional class.

Protein aggregation is geared by aggregation-prone regions that self-associate by ß-strand interactions. Charged residues and prolines are enriched at the flanks of aggregation-prone regions resulting in decreased aggregation. It is still unclear what drives the overrepresentation of these "aggregation gatekeepers", that is, whether their presence results from structural constraints determining protein stability or whether they constitute a bona fide functional class selectively maintained to control protein aggregation. As functional residues are typically conserved regardless of their cost to protein stability, we compared sequence conservation and thermodynamic cost of these residues in 2659 protein families in Escherichia coli. Across protein families, we find gatekeepers to be under strong selective conservation while at the same time representing a significant thermodynamic cost to protein structure. This finding supports the notion that aggregation gatekeepers are not structurally determined but evolutionary selected to control protein aggregation.

α-Galactosidase aggregation is a determinant of pharmacological chaperone efficacy on Fabry disease mutants.

Fabry disease is a lysosomal storage disorder caused by loss of α-galactosidase function. More than 500 Fabry disease mutants have been identified, the majority of which are structurally destabilized. A therapeutic strategy under development for lysosomal storage diseases consists of using pharmacological chaperones to stabilize the structure of the mutant protein, thereby promoting lysosomal delivery over retrograde degradation. The substrate analog 1-deoxygalactonojirimycin (DGJ) has been shown to restore activity of mutant α-galactosidase and is currently in clinical trial for treatment of Fabry disease. However, only ∼65% of tested mutants respond to treatment in cultured patient fibroblasts, and the structural underpinnings of DGJ response remain poorly explained. Using computational modeling and cell culture experiments, we show that the DGJ response is negatively affected by protein aggregation of α-galactosidase mutants, revealing a qualitative difference between misfolding-associated and aggregation-associated loss of function. A scoring function combining predicted thermodynamic stability and intrinsic aggregation propensity of mutants captures well their aggregation behavior under overexpression in HeLa cells. Interestingly, the same classifier performs well on DGJ response data of patient-derived cultured lymphoblasts, showing that protein aggregation is an important determinant of chemical chaperone efficiency under endogenous expression levels as well. Our observations reinforce the idea that treatment of aggregation-associated loss of function observed for the more severe α-galactosidase mutants could be enhanced by combining pharmacological chaperone treatment with the suppression of mutant aggregation, e.g. via proteostatic regulator compounds that increase cellular chaperone expression.

SNPeffect 4.0: on-line prediction of molecular and structural effects of protein-coding variants.

Single nucleotide variants (SNVs) are, together with copy number variation, the primary source of variation in the human genome and are associated with phenotypic variation such as altered response to drug treatment and susceptibility to disease. Linking structural effects of non-synonymous SNVs to functional outcomes is a major issue in structural bioinformatics. The SNPeffect database (http://snpeffect.switchlab.org) uses sequence- and structure-based bioinformatics tools to predict the effect of protein-coding SNVs on the structural phenotype of proteins. It integrates aggregation prediction (TANGO), amyloid prediction (WALTZ), chaperone-binding prediction (LIMBO) and protein stability analysis (FoldX) for structural phenotyping. Additionally, SNPeffect holds information on affected catalytic sites and a number of post-translational modifications. The database contains all known human protein variants from UniProt, but users can now also submit custom protein variants for a SNPeffect analysis, including automated structure modeling. The new meta-analysis application allows plotting correlations between phenotypic features for a user-selected set of variants.

A graphical interface for the FoldX forcefield.

SUMMARY: A graphical user interface for the FoldX protein design program has been developed as a plugin for the YASARA molecular graphics suite. The most prominent FoldX commands such as free energy difference upon mutagenesis and interaction energy calculations can now be run entirely via a windowed menu system and the results are immediately shown on screen. AVAILABILITY AND IMPLEMENTATION: The plugin is written in Python and is freely available for download at http://foldxyasara.switchlab.org/ and supported on Linux, MacOSX and MS Windows.

Redox proteomics of protein-bound methionine oxidation.

We here present a new method to measure the degree of protein-bound methionine sulfoxide formation at a proteome-wide scale. In human Jurkat cells that were stressed with hydrogen peroxide, over 2000 oxidation-sensitive methionines in more than 1600 different proteins were mapped and their extent of oxidation was quantified. Meta-analysis of the sequences surrounding the oxidized methionine residues revealed a high preference for neighboring polar residues. Using synthetic methionine sulfoxide containing peptides designed according to the observed sequence preferences in the oxidized Jurkat proteome, we discovered that the substrate specificity of the cellular methionine sulfoxide reductases is a major determinant for the steady-state of methionine oxidation. This was supported by a structural modeling of the MsrA catalytic center. Finally, we applied our method onto a serum proteome from a mouse sepsis model and identified 35 in vivo methionine oxidation events in 27 different proteins.

Increased monomerization of mutant HSPB1 leads to protein hyperactivity in Charcot-Marie-Tooth neuropathy.

Small heat shock proteins are molecular chaperones capable of maintaining denatured proteins in a folding-competent state. We have previously shown that missense mutations in the small heat shock protein HSPB1 (HSP27) cause distal hereditary motor neuropathy and axonal Charcot-Marie-Tooth disease. Here we investigated the biochemical consequences of HSPB1 mutations that are known to cause peripheral neuropathy. In contrast to other chaperonopathies, our results revealed that particular HSPB1 mutations presented higher chaperone activity compared with wild type. Hyperactivation of HSPB1 was accompanied by a change from its wild-type dimeric state to a monomer without dissociation of the 24-meric state. Purification of protein complexes from wild-type and HSPB1 mutants showed that the hyperactive isoforms also presented enhanced binding to client proteins. Furthermore, we show that the wild-type HSPB1 protein undergoes monomerization during heat-shock activation, strongly suggesting that the monomer is the active form of the HSPB1 protein.

Proteome-wide substrate analysis indicates substrate exclusion as a mechanism to generate caspase-7 versus caspase-3 specificity.

Caspase-3 and -7 are considered functionally redundant proteases with similar proteolytic specificities. We performed a proteome-wide screen on a mouse macrophage lysate using the N-terminal combined fractional diagonal chromatography technology and identified 46 shared, three caspase-3-specific, and six caspase-7-specific cleavage sites. Further analysis of these cleavage sites and substitution mutation experiments revealed that for certain cleavage sites a lysine at the P5 position contributes to the discrimination between caspase-7 and -3 specificity. One of the caspase-7-specific substrates, the 40 S ribosomal protein S18, was studied in detail. The RPS18-derived P6-P5' undecapeptide retained complete specificity for caspase-7. The corresponding P6-P1 hexapeptide still displayed caspase-7 preference but lost strict specificity, suggesting that P' residues are additionally required for caspase-7-specific cleavage. Analysis of truncated peptide mutants revealed that in the case of RPS18 the P4-P1 residues constitute the core cleavage site but that P6, P5, P2', and P3' residues critically contribute to caspase-7 specificity. Interestingly, specific cleavage by caspase-7 relies on excluding recognition by caspase-3 and not on increasing binding for caspase-7.

Accurate prediction of DnaK-peptide binding via homology modelling and experimental data.

Molecular chaperones are essential elements of the protein quality control machinery that governs translocation and folding of nascent polypeptides, refolding and degradation of misfolded proteins, and activation of a wide range of client proteins. The prokaryotic heat-shock protein DnaK is the E. coli representative of the ubiquitous Hsp70 family, which specializes in the binding of exposed hydrophobic regions in unfolded polypeptides. Accurate prediction of DnaK binding sites in E. coli proteins is an essential prerequisite to understand the precise function of this chaperone and the properties of its substrate proteins. In order to map DnaK binding sites in protein sequences, we have developed an algorithm that combines sequence information from peptide binding experiments and structural parameters from homology modelling. We show that this combination significantly outperforms either single approach. The final predictor had a Matthews correlation coefficient (MCC) of 0.819 when assessed over the 144 tested peptide sequences to detect true positives and true negatives. To test the robustness of the learning set, we have conducted a simulated cross-validation, where we omit sequences from the learning sets and calculate the rate of repredicting them. This resulted in a surprisingly good MCC of 0.703. The algorithm was also able to perform equally well on a blind test set of binders and non-binders, of which there was no prior knowledge in the learning sets. The algorithm is freely available at http://limbo.vib.be.

Analysis of protein processing by N-terminal proteomics reveals novel species-specific substrate determinants of granzyme B orthologs.

Using a targeted peptide-centric proteomics approach, we performed in vitro protease substrate profiling of the apoptotic serine protease granzyme B resulting in the delineation of more than 800 cleavage sites in 322 human and 282 mouse substrates, encompassing the known substrates Bid, caspase-7, lupus La protein, and fibrillarin. Triple SILAC (stable isotope labeling by amino acids in cell culture) further permitted intra-experimental evaluation of species-specific variations in substrate selection by the mouse or human granzyme B ortholog. For the first time granzyme B substrate specificities were directly mapped on a proteomic scale and revealed unknown cleavage specificities, uncharacterized extended specificity profiles, and macromolecular determinants in substrate selection that were confirmed by molecular modeling. We further tackled a substrate hunt in an in vivo setup of natural killer cell-mediated cell death confirming in vitro characterized granzyme B cleavages next to several other unique and hitherto unreported proteolytic events in target cells.