The distribution of residues in a polypeptide sequence is a determinant of aggregation optimized by evolution

Unveiling the genetic tapestry: Rare disease genomics of spinal muscular atrophy and phenylketonuria proteins

Rare diseases, defined by their low prevalence, present significant challenges, including delayed detection, expensive treatments, and limited research. This study delves into the genetic basis of two noteworthy rare diseases in Saudi Arabia: Phenylketonuria (PKU) and Spinal Muscular Atrophy (SMA). PKU, resulting from mutations in the phenylalanine hydroxylase (PAH) gene, exhibits geographical variability and impacts intellectual abilities. SMA, characterized by motor neuron loss, is linked to mutations in the survival of motor neuron 1 (SMN1) gene. Recognizing the importance of unveiling signature genomics in rare diseases, we conducted a quantitative study on PAH and SMN1 proteins of multiple organisms by employing various quantitative techniques to assess genetic variations. The derived signature-genomics contributes to a deeper understanding of these critical genes, paving the way for enhanced diagnostics for disorders associated with PAH and SMN1.

A3D Model Organism Database (A3D-MODB): a database for proteome aggregation predictions in model organisms

Protein aggregation has been associated with aging and different pathologies and represents a bottleneck in the industrial production of biotherapeutics. Numerous past studies performed in Escherichia coli and other model organisms have allowed to dissect the biophysical principles underlying this process. This knowledge fuelled the development of computational tools, such as Aggrescan 3D (A3D) to forecast and re-design protein aggregation. Here, we present the A3D Model Organism Database (A3D-MODB) http://biocomp.chem.uw.edu.pl/A3D2/MODB, a comprehensive resource for the study of structural protein aggregation in the proteomes of 12 key model species spanning distant biological clades. In addition to A3D predictions, this resource incorporates information useful for contextualizing protein aggregation, including membrane protein topology and structural model confidence, as an indirect reporter of protein disorder. The database is openly accessible without any need for registration. We foresee A3D-MOBD evolving into a central hub for conducting comprehensive, multi-species analyses of protein aggregation, fostering the development of protein-based solutions for medical, biotechnological, agricultural and industrial applications.

Sequence Versus Composition: What Prescribes IDP Biophysical Properties?

Intrinsically disordered proteins (IDPs) represent a distinct class of proteins and are distinguished from globular proteins by conformational plasticity, high evolvability and a broad functional repertoire. Some of their properties are reminiscent of early proteins, but their abundance in eukaryotes, functional properties and compositional bias suggest that IDPs appeared at later evolutionary stages. The spectrum of IDP properties and their determinants are still not well defined. This study compares rudimentary physicochemical properties of IDPs and globular proteins using bioinformatic analysis on the level of their native sequences and random sequence permutations, addressing the contributions of composition versus sequence as determinants of the properties. IDPs have, on average, lower predicted secondary structure contents and aggregation propensities and biased amino acid compositions. However, our study shows that IDPs exhibit a broad range of these properties. Induced fold IDPs exhibit very similar compositions and secondary structure/aggregation propensities to globular proteins, and can be distinguished from unfoldable IDPs based on analysis of these sequence properties. While amino acid composition seems to be a major determinant of aggregation and secondary structure propensities, sequence randomization does not result in dramatic changes to these properties, but for both IDPs and globular proteins seems to fine-tune the tradeoff between folding and aggregation.

A Shift in Aggregation Avoidance Strategy Marks a Long-Term Direction to Protein Evolution

To detect a direction to evolution, without the pitfalls of reconstructing ancestral states, we need to compare "more evolved" to "less evolved" entities. But because all extant species have the same common ancestor, none are chronologically more evolved than any other. However, different gene families were born at different times, allowing us to compare young protein-coding genes to those that are older and hence have been evolving for longer. To be retained during evolution, a protein must not only have a function, but must also avoid toxic dysfunction such as protein aggregation. There is conflict between the two requirements: hydrophobic amino acids form the cores of protein folds, but also promote aggregation. Young genes avoid strongly hydrophobic amino acids, which is presumably the simplest solution to the aggregation problem. Here we show that young genes' few hydrophobic residues are clustered near one another along the primary sequence, presumably to assist folding. The higher aggregation risk created by the higher hydrophobicity of older genes is counteracted by more subtle effects in the ordering of the amino acids, including a reduction in the clustering of hydrophobic residues until they eventually become more interspersed than if distributed randomly. This interspersion has previously been reported to be a general property of proteins, but here we find that it is restricted to old genes. Quantitatively, the index of dispersion delineates a gradual trend, i.e., a decrease in the clustering of hydrophobic amino acids over billions of years.

Loss of charge mutations in solvent exposed Lys residues of superoxide dismutase 1 do not induce inclusion formation in cultured cell models

Mutations in superoxide dismutase 1 (SOD1) associated with familial amyotrophic lateral sclerosis (fALS) induce the protein to misfold and aggregate. Missense mutations at more than 80 different amino acid positions have been associated with disease. How these mutations heighten the propensity of SOD1 to misfold and aggregate is unclear. With so many mutations, it is possible that more than one mechanism of aggregation may be involved. Of many possible mechanisms to explain heightened aggregation, one that has been suggested is that mutations that eliminate charged amino acids could diminish repulsive forces that would inhibit aberrant protein:protein interactions. Mutations at twenty-one charged residues in SOD1 have been associated with fALS, but of the 11 Lys residues in the protein, only 1 has been identified as mutated in ALS patients. Here, we examined whether loss of positively charged surface Lys residues in SOD1 would induce misfolding and formation of intracellular inclusions. We mutated four different Lys residues (K30, K36, K75, K91) in SOD1 that are not particularly well conserved, and expressed these variants as fusion proteins with yellow fluorescent protein (YFP) to assess inclusion formation. We also assessed whether these mutations induced binding to a conformation-restricted SOD1 antibody, designated C4F6, which recognizes non-natively folded protein. Although we observed some mutations to cause enhanced C4F6 binding, we did not observe that mutations that reduce charge at these positions caused the protein to form intracellular inclusions. Our findings may have implications for the low frequency of mutations at Lys residues SOD1 in ALS patients.

Theory of Sequence Effects in Amyloid Aggregation

We present a simple model for the effect of amino acid sequences on amyloid fibril formation. Using the HP model we find the binding lifetimes of four simple sequences by solving the first passage time for the intermolecular H-bond reaction coordinate. We find that sequences with identical binding energies have widely varying binding times depending on where the aggregation prone amino acids are located in the sequence. In general, longer binding times occur when the aggregation prone amino acids are clustered in a single "hot spot". Similarly, binding times are shortened by clustering weakly bound residues. Both of these effects are explained by an increase in the multiplicity of unbinding trajectories that comes from adding weak binding residues. Our model predicts a transition from ordered to disordered fibrils as the concentration of monomers increases. We apply our model to Aβ, IAPP, and apomyoglobin using binding energy estimates derived from bioinformatics. We find that these sequences are highly selective of the in-register state. This selectivity arises from the having strongly bound segments of varying length and separation.

Hidden Structural Codes in Protein Intrinsic Disorder

Intrinsic disorder is a major structural category in biology, accounting for more than 30% of coding regions across the domains of life, yet consists of conformational ensembles in equilibrium, a major challenge in protein chemistry. Anciently evolved papillomavirus genomes constitute an unparalleled case for sequence to structure-function correlation in cases in which there are no folded structures. E7, the major transforming oncoprotein of human papillomaviruses, is a paradigmatic example among the intrinsically disordered proteins. Analysis of a large number of sequences of the same viral protein allowed for the identification of a handful of residues with absolute conservation, scattered along the sequence of its N-terminal intrinsically disordered domain, which intriguingly are mostly leucine residues. Mutation of these led to a pronounced increase in both α-helix and β-sheet structural content, reflected by drastic effects on equilibrium propensities and oligomerization kinetics, and uncovers the existence of local structural elements that oppose canonical folding. These folding relays suggest the existence of yet undefined hidden structural codes behind intrinsic disorder in this model protein. Thus, evolution pinpoints conformational hot spots that could have not been identified by direct experimental methods for analyzing or perturbing the equilibrium of an intrinsically disordered protein ensemble.

Effect of position-specific single-point mutations and biophysical characterization of amyloidogenic peptide fragments identified from lattice corneal dystrophy patients

Corneal stromal dystrophies are a group of genetic disorders that may be caused by mutations in the transforming growth factor β-induced (TGFBI) gene which results in the aggregation and deposition of mutant proteins in various layers of the cornea. The type of amino acid substitution dictates the age of onset, anatomical location of the deposits, morphological features of deposits (amyloid, amorphous powder or a mixture of both forms) and the severity of disease presentation. It has been suggested that abnormal turnover and aberrant proteolytic processing of the mutant proteins result in the accumulation of insoluble protein deposits. Using mass spectrometry, we identified increased abundance of a 32 amino acid-long peptide in the 4th fasciclin-like domain-1 (FAS-1) domain of transforming growth factor β-induced protein (amino acid 611-642) in the amyloid deposits of the patients with lattice corneal dystrophies (LCD). In vitro studies demonstrated that the peptide readily formed amyloid fibrils under physiological conditions. Clinically relevant substitution (M619K, N622K, N622H, G623R and H626R) of the truncated peptide resulted in profound changes in the kinetics of amyloid formation, thermal stability of the amyloid fibrils and cytotoxicity of fibrillar aggregates, depending on the position and the type of the amino acid substitution. The results suggest that reduction in the overall net charge, nature and position of cationic residue substitution determines the amyloid aggregation propensity and thermal stability of amyloid fibrils.

From the Evolution of Protein Sequences Able to Resist Self-Assembly to the Prediction of Aggregation Propensity

Folding of polypeptide chains into biologically active entities is an astonishingly complex process, determined by the nature and the sequence of residues emerging from ribosomes. While it has been long believed that evolution has pressed genomes so that specific sequences could adopt unique, functional three-dimensional folds, it is now clear that complex protein machineries act as quality control system and supervise folding. Notwithstanding that, events such as erroneous folding, partial folding, or misfolding are frequent during the life of a cell or a whole organism, and they can escape controls. One of the possible outcomes of this misbehavior is cross-β aggregation, a super secondary structure which represents the hallmark of self-assembled, well organized, and extremely ordered structures termed amyloid fibrils. What if evolution would have not taken into account such possibilities? Twenty years of research point toward the idea that, in fact, evolution has constantly supervised the risk of errors and minimized their impact. In this review we tried to survey the major findings in the amyloid field, trying to describe what the real pitfalls of protein folding are-from an evolutionary perspective-and how sequence and structural features have evolved to balance the need for perfect, dynamic, functionally efficient structures, and the detrimental effects implicit in the dangerous process of folding. We will discuss how the knowledge obtained from these studies has been employed to produce computational methods able to assess, predict, and discriminate the aggregation properties of protein sequences.

Antibody Aggregation: Insights from Sequence and Structure

Monoclonal antibodies (mAbs) are the fastest-growing biological therapeutics with important applications ranging from cancers, autoimmunity diseases and metabolic disorders to emerging infectious diseases. Aggregation of mAbs continues to be a major problem in their developability. Antibody aggregation could be triggered by partial unfolding of its domains, leading to monomer-monomer association followed by nucleation and growth. Although the aggregation propensities of antibodies and antibody-based proteins can be affected by the external experimental conditions, they are strongly dependent on the intrinsic antibody properties as determined by their sequences and structures. In this review, we describe how the unfolding and aggregation susceptibilities of IgG could be related to their cognate sequences and structures. The impact of antibody domain structures on thermostability and aggregation propensities, and effective strategies to reduce aggregation are discussed. Finally, the aggregation of antibody-drug conjugates (ADCs) as related to their sequence/structure, linker payload, conjugation chemistry and drug-antibody ratio (DAR) is reviewed.

Comparison of the aggregation of homologous β2-microglobulin variants reveals protein solubility as a key determinant of amyloid formation

The mouse and human β2-microglobulin protein orthologs are 70% identical in sequence and share 88% sequence similarity. These proteins are predicted by various algorithms to have similar aggregation and amyloid propensities. However, whilst human β2m (hβ2m) forms amyloid-like fibrils in denaturing conditions (e.g. pH2.5) in the absence of NaCl, mouse β2m (mβ2m) requires the addition of 0.3M NaCl to cause fibrillation. Here, the factors which give rise to this difference in amyloid propensity are investigated. We utilise structural and mutational analyses, fibril growth kinetics and solubility measurements under a range of pH and salt conditions, to determine why these two proteins have different amyloid propensities. The results show that, although other factors influence the fibril growth kinetics, a striking difference in the solubility of the proteins is a key determinant of the different amyloidogenicity of hβ2m and mβ2m. The relationship between protein solubility and lag time of amyloid formation is not captured by current aggregation or amyloid prediction algorithms, indicating a need to better understand the role of solubility on the lag time of amyloid formation. The results demonstrate the key contribution of protein solubility in determining amyloid propensity and lag time of amyloid formation, highlighting how small differences in protein sequence can have dramatic effects on amyloid formation.

Computational analysis of candidate prion-like proteins in bacteria and their role

Prion proteins were initially associated with diseases such as Creutzfeldt Jakob and transmissible spongiform encephalopathies. However, deeper research revealed them as versatile tools, exploited by the cells to execute fascinating functions, acting as epigenetic elements or building membrane free compartments in eukaryotes. One of the most intriguing properties of prion proteins is their ability to propagate a conformational assembly, even across species. In this context, it has been observed that bacterial amyloids can trigger the formation of protein aggregates by interacting with host proteins. As our life is closely linked to bacteria, either through a parasitic or symbiotic relationship, prion-like proteins produced by bacterial cells might play a role in this association. Bioinformatics is helping us to understand the factors that determine conformational conversion and infectivity in prion-like proteins. We have used PrionScan to detect prion domains in 839 different bacteria proteomes, detecting 2200 putative prions in these organisms. We studied this set of proteins in order to try to understand their functional role and structural properties. Our results suggest that these bacterial polypeptides are associated to peripheral rearrangement, macromolecular assembly, cell adaptability, and invasion. Overall, these data could reveal new threats and therapeutic targets associated to infectious diseases.

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.

On the role of aggregation prone regions in protein evolution, stability, and enzymatic catalysis: insights from diverse analyses

The various roles that aggregation prone regions (APRs) are capable of playing in proteins are investigated here via comprehensive analyses of multiple non-redundant datasets containing randomly generated amino acid sequences, monomeric proteins, intrinsically disordered proteins (IDPs) and catalytic residues. Results from this study indicate that the aggregation propensities of monomeric protein sequences have been minimized compared to random sequences with uniform and natural amino acid compositions, as observed by a lower average aggregation propensity and fewer APRs that are shorter in length and more often punctuated by gate-keeper residues. However, evidence for evolutionary selective pressure to disrupt these sequence regions among homologous proteins is inconsistent. APRs are less conserved than average sequence identity among closely related homologues (≥80% sequence identity with a parent) but APRs are more conserved than average sequence identity among homologues that have at least 50% sequence identity with a parent. Structural analyses of APRs indicate that APRs are three times more likely to contain ordered versus disordered residues and that APRs frequently contribute more towards stabilizing proteins than equal length segments from the same protein. Catalytic residues and APRs were also found to be in structural contact significantly more often than expected by random chance. Our findings suggest that proteins have evolved by optimizing their risk of aggregation for cellular environments by both minimizing aggregation prone regions and by conserving those that are important for folding and function. In many cases, these sequence optimizations are insufficient to develop recombinant proteins into commercial products. Rational design strategies aimed at improving protein solubility for biotechnological purposes should carefully evaluate the contributions made by candidate APRs, targeted for disruption, towards protein structure and activity.

Biotechnology requires the large-scale expression, yield, and storage of recombinant proteins. Each step in protein production has the potential to cause aggregation as proteins, not evolved to exist outside the cell, endure the various steps involved in commercial manufacturing processes. Mechanistic studies into protein aggregation have revealed that certain sequence regions contribute more to the aggregation propensity of a protein than other sequence regions do. Efforts to disrupt these regions have thus far indicated that rational sequence engineering is a useful technique to reduce the aggregation of biotechnologically relevant proteins. To improve our ability to rationally engineer proteins with enhanced expression, solubility, and shelf-life we conducted extensive analyses of aggregation prone regions (APRs) within protein sequences to characterize the various roles these regions play in proteins. Findings from this work indicate that protein sequences have evolved by minimizing their aggregation propensities. However, we also found that many APRs are conserved in protein families and are essential to maintain protein stability and function. Therefore, the contributions that APRs, targeted for disruption, make towards protein stability and function should be carefully evaluated when improving protein solubility via rational design.

Misfolding and amyloid aggregation of apomyoglobin

Apomyoglobin is an excellent example of a monomeric all α-helical globular protein whose folding pathway has been extensively studied and well characterized. Structural perturbation induced by denaturants or high temperature as well as amino acid substitution have been described to induce misfolding and, in some cases, aggregation. In this article, we review the molecular mechanism of the aggregation process through which a misfolded form of a mutated apomyoglobin aggregates at physiological pH and room temperature forming an amyloid fibril. The results are compared with data showing that either amyloid or aggregate formation occurs under particular denaturing conditions or upon cleavage of the residues corresponding to the C-terminal helix of apomyoglobin. The results are discussed in terms of the sequence regions that are more important than others in determining the amyloid aggregation process.

Analyses of the general rule on residue pair frequencies in local amino acid sequences of soluble, ordered proteins

The amino acid sequences of soluble, ordered proteins with stable structures have evolved due to biological and physical requirements, thus distinguishing them from random sequences. Previous analyses have focused on extracting the features that frequently appear in protein substructures, such as α-helix and β-sheet, but the universal features of protein sequences have not been addressed. To clarify the differences between native protein sequences and random sequences, we analyzed 7368 soluble, ordered protein sequences, by inspecting the observed and expected occurrences of 400 amino acid pairs in local proximity, up to 10 residues along the sequence in comparison with their expected occurrence in random sequence. We found the trend that the hydrophobic residue pairs and the polar residue pairs are significantly decreased, whereas the pairs between a hydrophobic residue and a polar residue are increased. This trend was universally observed regardless of the secondary structure content but was not observed in protein sequences that include intrinsically disordered regions, indicating that it can be a general rule of protein foldability. The possible benefits of this rule are discussed from the viewpoints of protein aggregation and disorder, which are both caused by low-complexity regions of hydrophobic or polar residues.

Evolutionary selection for protein aggregation

Protein aggregation is being found to be associated with an increasing number of human diseases. Aggregation can lead to a loss of function (lack of active protein) or to a toxic gain of function (cytotoxicity associated with protein aggregates). Although potentially harmful, protein sequences predisposed to aggregation seem to be ubiquitous in all kingdoms of life, which suggests an evolutionary advantage to having such segments in polypeptide sequences. In fact, aggregation-prone segments are essential for protein folding and for mediating certain protein-protein interactions. Moreover, cells use protein aggregates for a wide range of functions. Against this background, life has adapted to tolerate the presence of potentially dangerous aggregation-prone sequences by constraining and counteracting the aggregation process. In the present review, we summarize the current knowledge of the advantages associated with aggregation-prone stretches in proteomes and the strategies that cellular systems have developed to control the aggregation process.

Resolution of the effects induced by W → F substitutions on the conformation and dynamics of the amyloid-forming apomyoglobin mutant W7FW14F

Myoglobin is an alpha-helical globular protein containing two highly conserved tryptophanyl residues at positions 7 and 14 in the N-terminal region. The simultaneous substitution of the two residues increases the susceptibility of the polypeptide chain to misfold, causing amyloid aggregation under physiological condition, i.e., neutral pH and room temperature. The role played by tryptophanyl residues in driving the folding process has been investigated by examining three mutated apomyoglobins, i.e., W7F, W14F, and the amyloid-forming mutant W7FW14F, by an integrated approach based on far-ultraviolet (UV) circular dichroism (CD) analysis, fluorescence spectroscopy, and complementary proteolysis. Particular attention has been devoted to examine the conformational and dynamic properties of the equilibrium intermediate formed at pH 4.0, since it represents the early organized structure from which the native fold originates. The results show that the W → F substitutions at position 7 and 14 differently affect the structural organization of the AGH subdomain of apomyoglobin. The combined effect of the two substitutions in the double mutant impairs the formation of native-like contacts and favors interchain interactions, leading to protein aggregation and amyloid formation.

Aggregation gatekeepers modulate protein homeostasis of aggregating sequences and affect bacterial fitness

The most common mechanism by which proteins aggregate consists in the assembly of short hydrophobic primary sequence segments into extended β-structured agglomerates. A significant enrichment of charged residues is observed at the flank of these aggregation-prone sequence segments, suggesting selective pressure against aggregation. These so-called aggregation gatekeepers act by increasing the intrinsic solubility of aggregating sequences in vitro, but it has been suggested that they could also facilitate chaperone interactions. Here, we address whether aggregation gatekeepers affect bacterial fitness. In Escherichia coli MC4100 we overexpressed GFP fusions with an aggregation-prone segment of σ32 (further termed σ32β) flanked by gatekeeper and non-gatekeeper residues and measured pairwise competitive growth. We found that the identity of flanking residues had significant effect on bacterial growth. Overexpression of σ32β flanked by its natural gatekeepers displayed the greatest competitive fitness, followed by other combinations of gatekeepers, while absence of gatekeepers strongly affects bacterial fitness. Further analysis showed the diversity of effects of gatekeepers on the proteostasis of σ32β including synthesis and degradation rates, in vivo aggregation propensity and chaperone response. Our results suggest that gatekeeper residues affect bacterial fitness not only by modulating the intrinsic aggregation propensity of proteins but also by the manner in which they affect the processing of σ32β-GFP by the protein quality control machinery of the cell. In view of these observations, we hypothesize that variation at gatekeeper positions offers a flexible selective strategy to modulate the proteostatic regulation of proteins to the match intrinsic aggregation propensities of proteins with required expression levels.

The role of tryptophan in protein fibrillogenesis: relevance of Trp7 and Trp14 to the amyloidogenic properties of myoglobin

In order to understand the role of tryptophan in the mechanisms of fibrils formation, the ability of a series of analogs of the residue 7-18 span of myoglobin to form amyloid-like fibrils was investigated. Alternatively one or both tryptophans were substituted with alanine and leucine, to determine the contribution of hydrophobicity and aromaticity. The scale of aggregation propensity of the peptides determined indicates that tryptophan is crucial for the amyloidogenic process. Since the rare tryptophan residue is generally engaged in structural roles in proteins, or when exposed serves as binding sites, we surmise that its exposure in the amyloidogenic fragments allows for intermolecular clustering with residues from other molecules leading to the formation of amyloid aggregates.

AGGRESCAN: method, application, and perspectives for drug design

Protein aggregation underlies the development of an increasing number of conformational human diseases of growing incidence, such as Alzheimer's and Parkinson's diseases. Furthermore, the accumulation of recombinant proteins as intracellular aggregates represents a critical obstacle for the biotechnological production of polypeptides. Also, ordered protein aggregates constitute novel and versatile nanobiomaterials. Consequently, there is an increasing interest in the development of methods able to forecast the aggregation properties of polypeptides in order to modulate their intrinsic solubility. In this context, we have developed AGGRESCAN, a simple and fast algorithm that predicts aggregation-prone segments in protein sequences, compares the aggregation properties of different proteins or protein sets and analyses the effect of mutations on protein aggregation propensities.

An evolutionary trade-off between protein turnover rate and protein aggregation favors a higher aggregation propensity in fast degrading proteins

We previously showed the existence of selective pressure against protein aggregation by the enrichment of aggregation-opposing ‘gatekeeper’ residues at strategic places along the sequence of proteins. Here we analyzed the relationship between protein lifetime and protein aggregation by combining experimentally determined turnover rates, expression data, structural data and chaperone interaction data on a set of more than 500 proteins. We find that selective pressure on protein sequences against aggregation is not homogeneous but that short-living proteins on average have a higher aggregation propensity and fewer chaperone interactions than long-living proteins. We also find that short-living proteins are more often associated to deposition diseases. These findings suggest that the efficient degradation of high-turnover proteins is sufficient to preclude aggregation, but also that factors that inhibit proteasomal activity, such as physiological ageing, will primarily affect the aggregation of short-living proteins.

In order to carry out their biological function, proteins need to fold into well-defined three-dimensional structures. Protein aggregation is a process whereby proteins misfold into inactive and often toxic higher order structures, which is implied in about 30 human diseases such as Alzheimer's disease, Parkinson's disease and systemic amyloidosis. In earlier work it has been shown that although protein aggregation is an intrinsic property of polypeptide chains that cannot be entirely avoided, evolution has optimized protein sequences to minimize the risk of aggregation in a proteome. Here we show that this pressure is not uniform, but that proteins with a short lifetime have on average a higher aggregation propensity than long-living proteins. In addition, we show that high turnover proteins also make fewer interactions with chaperones. Taken together, these observations suggest that under normal physiological conditions the aggregation propensity of short-lived proteins does not represent a significant treat for the biochemistry of the cell. Presumably the strong dependence of these proteins on proteasomal degradation is sufficient to preclude the accumulation of aggregates. As proteasomal activity declines with age this would also explain why we observe a higher association of high turnover proteins with age-dependent aggregation-related diseases.

Proline periodicity modulates the self-assembly properties of elastin-like polypeptides

Elastin is a self-assembling protein of the extracellular matrix that provides tissues with elastic extensibility and recoil. The monomeric precursor, tropoelastin, is highly hydrophobic yet remains substantially disordered and flexible in solution, due in large part to a high combined threshold of proline and glycine residues within hydrophobic sequences. In fact, proline-poor elastin-like sequences are known to form amyloid-like fibrils, rich in β-structure, from solution. On this basis, it is clear that hydrophobic elastin sequences are in general optimized to avoid an amyloid fate. However, a small number of hydrophobic domains near the C terminus of tropoelastin are substantially depleted of proline residues. Here we investigated the specific contribution of proline number and spacing to the structure and self-assembly propensities of elastin-like polypeptides. Increasing the spacing between proline residues significantly decreased the ability of polypeptides to reversibly self-associate. Real-time imaging of the assembly process revealed the presence of smaller colloidal droplets that displayed enhanced propensity to cluster into dense networks. Structural characterization showed that these aggregates were enriched in β-structure but unable to bind thioflavin-T. These data strongly support a model where proline-poor regions of the elastin monomer provide a unique contribution to assembly and suggest a role for localized β-sheet in mediating self-assembly interactions.

A computational approach for identifying the chemical factors involved in the glycosaminoglycans-mediated acceleration of amyloid fibril formation

Background

Amyloid fibril formation is the hallmark of many human diseases, including Alzheimer's disease, type II diabetes and amyloidosis. Amyloid fibrils deposit in the extracellular space and generally co-localize with the glycosaminoglycans (GAGs) of the basement membrane. GAGs have been shown to accelerate the formation of amyloid fibrils in vitro for a number of protein systems. The high number of data accumulated so far has created the grounds for the construction of a database on the effects of a number of GAGs on different proteins.

Methodology/Principal Findings

In this study, we have constructed such a database and have used a computational approach that uses a combination of single parameter and multivariate analyses to identify the main chemical factors that determine the GAG-induced acceleration of amyloid formation. We show that the GAG accelerating effect is mainly governed by three parameters that account for three-fourths of the observed experimental variability: the GAG sulfation state, the solute molarity, and the ratio of protein and GAG molar concentrations. We then combined these three parameters into a single equation that predicts, with reasonable accuracy, the acceleration provided by a given GAG in a given condition.

Conclusions/Significance

In addition to shedding light on the chemical determinants of the protein∶GAG interaction and to providing a novel mathematical predictive tool, our findings highlight the possibility that GAGs may not have such an accelerating effect on protein aggregation under the conditions existing in the basement membrane, given the values of salt molarity and protein∶GAG molar ratio existing under such conditions.

Potential aggregation-prone regions in complementarity-determining regions of antibodies and their contribution towards antigen recognition: a computational analysis

PURPOSE

To analyze contribution of short aggregation-prone regions (APRs), which may self-associate via cross-beta motif and were earlier identified in therapeutic mAbs, towards antigen recognition via structural analyses of antibody-antigen complexes.

METHODS

A dataset of 29 publically available high-resolution crystal structures of Fab-antigen complexes was collected. Contribution of APRs towards the surface areas of the Fabs buried by the cognate antigens was computed. Propensities of amino acids to occur in APRs and to be involved in antigen binding were compared. Coincidence between APRs and individual CDR loops was examined.

RESULTS

All Fabs in the dataset contain at least one APR in CDR loops and adjacent framework beta-strands. The average contribution of APRs towards buried surface area of Fabs is 16.0 +/- 10.7%. Aggregation and antigen recognition may be coupled via aromatic residues (Tyr, Trp), which occur with high propensities in both APRs and antigen binding sites. APRs are infrequent in the heavy chain CDR 3 (H3) loops (7%), but are frequent in H2 loops (45%).

CONCLUSIONS

Co-incidence of APRs with antigen recognition sites can potentially lead to the loss of function upon aggregation. Rational structure-based design or selection strategies are suggested for biotherapeutics with improved druggability while maintaining potency.

Protein sequences encode safeguards against aggregation

Functional requirements shaped proteins into globular structures. Under these structural constraints, which require both regular secondary structure and a hydrophobic core, protein aggregation is an unavoidable corollary to protein structure. However, as aggregation results in reduced fitness, natural selection will tend to eliminate strongly aggregating sequences. The analysis of distribution and variation of aggregation patterns in the human proteome using the TANGO algorithm confirms the findings of a previous study on several proteomes: the flanks of aggregation-prone regions are enriched with charged residues and proline, the so-called gatekeeper-residues. Moreover, in this study, we observed a widespread redundancy in gatekeeper usage. Interestingly, aggregating regions from key proteins such as p53 or huntingtin are among the most extensive "gatekept" sequences. As a consequence, mutations that remove gatekeepers could therefore result in a strong increase in disease-susceptibility. In a set of disease-associated mutations from the UniProt database, we find a strong enrichment of mutations that disrupt gatekeeper motifs. Closer inspection of a number of case studies indicates clearly that removing gatekeepers may play a determining role in widely varying disorders, such as van der Woude syndrome (VWS), X-linked Fabry disease (FD), and limb-girdle muscular dystrophy.

NetCSSP: web application for predicting chameleon sequences and amyloid fibril formation

The calculation of contact-dependent secondary structure propensity (CSSP) is a unique and sensitive method that detects non-native secondary structure propensities in protein sequences. This method has applications in predicting local conformational change, which typically is observed in core sequences of protein aggregation and amyloid fibril formation. NetCSSP implements the latest version of the CSSP algorithm and provides a Flash chart-based graphic interface that enables an interactive calculation of CSSP values for any user-selected regions in a given protein sequence. This feature also can quantitatively estimate the mutational effect on changes in native or non-native secondary structural propensities in local sequences. In addition, this web tool provides precalculated non-native secondary structure propensities for over 1 400 000 fragments that are seven-residues long, collected from PDB structures. They are searchable for chameleon subsequences that can serve as the core of amyloid fibril formation. The NetCSSP web tool is available at http://cssp2.sookmyung.ac.kr/.

Fibrils with parallel in-register structure constitute a major class of amyloid fibrils: molecular insights from electron paramagnetic resonance spectroscopy

The deposition of amyloid- and amyloid-like fibrils is the main pathological hallmark of numerous protein misfolding diseases including Alzheimer's disease, transmissible spongiform encephalopathy, and type 2 diabetes. Besides the well-established role in disease, recent work on a variety of organisms ranging from bacteria to humans suggests that amyloid fibrils can also convey biological functions. To better understand the molecular mechanisms by which amyloidogenic proteins misfold in disease or perform biological functions, structural information is essential. Although high-resolution structural analysis of amyloid fibrils has been challenging, a combination of biophysical approaches is beginning to unravel the various structural features of amyloid fibrils. Here we review these recent developments with particular emphasis on amyloid fibrils that have been studied using site-directed spin labeling and electron paramagnetic resonance spectroscopy. This approach has been used to define the precise location of fibril-forming core regions and identify local secondary structures within such core regions. Perhaps one of the most remarkable findings arrived at by site-directed spin labeling was that most fibrils that contain an extensive core region of 20 amino acids or more share a common parallel in-register arrangement of beta strands. The preference for this arrangement can be explained on topological grounds and may be rationalized by the maximization of hydrophobic contact surface.

Aggregation propensity of the human proteome

Formation of amyloid-like fibrils is involved in numerous human protein deposition diseases, but is also an intrinsic property of polypeptide chains in general. Progress achieved recently now allows the aggregation propensity of proteins to be analyzed over large scales. In this work we used a previously developed predictive algorithm to analyze the propensity of the 34,180 protein sequences of the human proteome to form amyloid-like fibrils. We show that long proteins have, on average, less intense aggregation peaks than short ones. Human proteins involved in protein deposition diseases do not differ extensively from the rest of the proteome, further demonstrating the generality of protein aggregation. We were also able to reproduce some of the results obtained with other algorithms, demonstrating that they do not depend on the type of computational tool employed. For example, proteins with different subcellular localizations were found to have different aggregation propensities, in relation to the various efficiencies of quality control mechanisms. Membrane proteins, intrinsically disordered proteins, and folded proteins were confirmed to have very different aggregation propensities, as a consequence of their different structures and cellular microenvironments. In addition, gatekeeper residues at strategic positions of the sequences were found to protect human proteins from aggregation. The results of these comparative analyses highlight the existence of intimate links between the propensity of proteins to form aggregates with β-structure and their biology. In particular, they emphasize the existence of a negative selection pressure that finely modulates protein sequences in order to adapt their aggregation propensity to their biological context.

Amyloid-like fibrils are insoluble proteinaceous fibrillar aggregates with a characteristic structure (the cross-β core) that form and deposit in more than 40 pathological conditions in humans. These include Alzheimer's disease, Parkinson's disease, type II diabetes, and the spongiform encephalopathies. A number of proteins not involved in any disease can also form amyloid-like fibrils in vitro, suggesting that amyloid fibril formation is an intrinsic property of proteins in general. Recent efforts in understanding the physico-chemical grounds of amyloid fibril formation has led to the development of several algorithms, capable of predicting a number of aggregation-related parameters of a protein directly from its amino acid sequence. In order to study the predicted aggregation behavior of the human proteome, we have run one of these algorithms on the 34,180 human protein sequences. Our results demonstrate that molecular evolution has acted on protein sequences to finely modulate their aggregation propensities, depending on different parameters related to their in vivo environment. Together with cellular control mechanisms, this natural selection protects proteins from aggregation during their lifetime.

Methionine oxidation inhibits assembly and promotes disassembly of apolipoprotein C-II amyloid fibrils

Methionine residues are linked to the pathogenicity of several amyloid diseases; however, the mechanism of this relationship is largely unknown. These diseases are characterized, in vivo, by the accumulation of insoluble proteinaceous plaques, of which the major constituents are amyloid fibrils. In vitro, methionine oxidation has been shown to modulate fibril assembly in several well-characterized amyloid systems. Human apolipoprotein (apo) C-II contains two methionine residues (Met-9 and Met-60) and readily self-assembles in vitro to form homogeneous amyloid fibrils, thus providing a convenient system to examine the effect of methionine oxidation on amyloid fibril formation and stability. Upon oxidation of the methionine residues of apoC-II with hydrogen peroxide, fibril formation was inhibited. Oxidized apoC-II molecules did not inhibit native apoC-II assembly, indicating that the oxidized molecules had a reduced ability to interact with the growing fibrils. Single Met-Val substitutions were performed and showed that oxidation of Met-60 had a more significant inhibitory effect than oxidation of Met-9. In addition, Met-Gln substitutions designed to mimic the effect of oxidation on side chain hydrophilicity showed that a change in hydrophobicity at position 60 within the core region of the fibril had a potent inhibitory effect. The oxidation of preformed apoC-II fibrils caused their dissociation; however, mutants in which the Met-60 was substituted with a valine were protected from this peroxide-induced dissociation. This work highlights an important role for methionine in the formation of amyloid fibril structure and gives new insight into how oxidation affects the stability of mature fibrils.

Natural selection against protein aggregation on self-interacting and essential proteins in yeast, fly, and worm

Protein aggregation is the phenomenon of protein self-association potentially leading to detrimental effects on physiology, which is closely related to numerous human diseases such as Alzheimer's and Parkinson's disease. Despite progress in understanding the mechanism of protein aggregation, how natural selection against protein aggregation acts on subunits of protein complexes and on proteins with different contributions to organism fitness remains largely unknown. Here, we perform a proteome-wide analysis by using an experimentally validated algorithm TANGO and utilizing sequence, interactomic and phenotype-based functional genomic data from yeast, fly, and nematode. We find that proteins that are capable of forming homooligomeric complex have lower aggregation propensity compared with proteins that do not function as homooligomer. Further, proteins that are essential to the fitness of an organism have lower aggregation propensity compared with nonessential ones. Our finding suggests that the selection force against protein aggregation acts across different hierarchies of biological system.

Prevention of amyloid-like aggregation as a driving force of protein evolution

Uncontrolled protein aggregation is a constant challenge in all compartments of living organisms. The failure of a peptide or protein to remain soluble often results in pathology. So far, more than 40 human diseases have been associated with the formation of extracellular fibrillar aggregates - known as amyloid fibrils - or structurally related intracellular deposits. It is well known that molecular chaperones and elaborate quality control mechanisms exist in the cell to counteract aggregation. However, an increasing number of reports during the past few years indicate that proteins have also evolved structural and sequence-based strategies to prevent aggregation. This review describes these strategies and the selection pressures that exist on protein sequences to combat their uncontrolled aggregation. We will describe the different types of mechanism evolved by proteins that adopt different conformational states including normally folded proteins, intrinsically disordered polypeptide chains, elastomeric systems and multimodular proteins.