Investigating protein aggregation in protein-carbohydrate interfaces using sequence and structural features.

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.

Yeast Proteins may Reversibly Aggregate like Amphiphilic Molecules

Colloidal stability is among the key challenges the pharmaceutical industry faces during the production and manufacturing of protein therapeutics. Self-association and aggregation processes can not only impair therapeutic efficacy but also induce immunogenic responses in patients. Aggregation-prone regions (APRs) consisting of hydrophobic patches are commonly identified as the source for colloidal instability, and rational strategies to mitigate aggregation propensity often require genetic engineering to eliminate hydrophobic amino acid residues. Here, we investigate cucurbit[7]uril (CB[7]), a water-soluble macrocycle able to form host-guest complexes with aromatic amino acid residues, as a potential excipient to mitigate protein aggregation propensity. Two monoclonal antibodies (mAbs), one harboring an APR and one lacking an APR, were first assessed for their colloidal stability (measured as the translational diffusion coefficient) in the presence and absence of CB[7] using dynamic light scattering. Due to the presence of a tryptophan residue within the APR, we were able to monitor changes in intrinsic fluorescence in response to increasing concentrations of CB[7]. Isothermal titration calorimetry and NMR spectroscopy were then used to characterize the putative host-guest interaction. Our results suggest a stabilizing effect of CB[7] on the aggregation-prone mAb, due to the specific interaction of CB[7] with aromatic amino acid residues located within the APR. This provides a starting point for exploring CB[7] as a candidate excipient for the formulation of aggregation-prone mAbs.

PurposeTo analyze contribution of short aggregation-prone regions (APRs), which may self-associate via cross-β motif and were earlier identified in therapeutic mAbs, towards antigen recognition via structural analyses of antibody-antigen complexes.MethodsA 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.ResultsAll Fabs in the dataset contain at least one APR in CDR loops and adjacent framework β-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%).ConclusionsCo-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.Electronic Supplementary MaterialThe online version of this article (doi:10.1007/s11095-010-0143-5) contains supplementary material, which is available to authorized users.

Most protein sequences contain one or several short aggregation prone regions (APR) that can nucleate protein aggregation. Under normal conditions these APRs are protected from aggregation by protein interactions or because they are buried in the hydrophobic core of native protein domains. However, mutation, physiological stress or age-related disregulation of protein homeostasis increases the probability that aggregation-nucleating regions become solvent exposed. Aggregation then results from the self-assembly of APRs into β-structured agglomerates that vary from small soluble oligomeric assemblies to large insoluble inclusions containing thousands of molecules. The functional effects of APR-driven aggregation are diverse and protein-specific leading to distinct disease phenotypes ranging from neurodegeneration to cancer. On a cellular and physiological level both wild type loss-of-function as well as aggregation-dependent gain-of-function effects have been shown to contribute to disease. Several molecular mechanism have been proposed to contribute to gain-of-function activity of protein aggregates including cellular membrane disregulation, saturation of the protein quality control machinery or the ability of aggregates to engage non-native interactions with proteins and nucleic acids. These different mechanisms will all, to some extent, contribute to gain-of-function as in essence they all contribute to the rewiring of the cellular interactome by aggregation-specific interactions, resulting for instance in the pronounced neurotoxicity of TDP43 aggregates by the sequestration of RNA molecules or the promotion of cell proliferation by the entrapment of homologous tumor suppressor proteins in p53 aggregates in cancer. In this review we discuss the mechanism of APR driven aggregation and how APRs contribute to modifying the cellular interactome by recruiting both misfolded as well as active proteins thereby inhibiting or activating specific cellular functions. Finally, we discuss the ubiquity of APRs in protein sequences and how selective pressure shaped protein sequences to minimize APR aggregation.

Aggregation of p53 mutants can result in loss-of-function, gain-of-function, and dominant-negative effects that contribute to tumor growth. Revealing the mechanisms underlying mutation-enhanced p53 aggregation and dissecting how small molecule inhibitors prevent the conversion of p53 into aggregation-primed conformations are fundamentally important for the development of novel therapeutics for p53 aggregation-associated cancers. A recent experimental study shows that resveratrol (RSV) has an inhibitory effect on the aggregation of hot-spot R248Q mutant of the p53 core domain (p53C), while pterostilbene (PT) exhibits a relatively poor inhibitory efficacy. However, the conformational properties of the R248Q mutant leading to its enhanced aggregation propensity and the inhibitory mechanism of RSV against p53C aggregation are not well understood. Herein, we performed extensive all-atom molecular dynamics simulations on R248Q p53C in the absence and presence of RSV/PT, as well as wild-type (WT) p53C. Our simulations reveal that loop L3, where the mutation resides, remains compact in WT p53C, while it becomes extended in the R248Q mutant. The extension of loop L3 weakens the interactions between loop L3 and two crucial aggregation-prone regions (APRs) of p53C, leading to impaired interactions within the APRs and their structural destabilization as well as p53C. The destabilized APRs in the R248Q mutant are more exposed than in WT p53C, which is conducive to p53C aggregation. RSV has a higher preference to bind to R248Q p53C than PT. This binding not only stabilizes loop L3 of R248Q mutant to its WT-like conformation, preventing L3-extension-caused APRs' destabilization but also reduces APRs' solvent exposure, thereby inhibiting R248Q p53C aggregation. However, PT exhibits a lower hydrogen-bonding capability and a higher self-association propensity, which would lead to a reduced p53C binding and a weakened inhibitory effect on R248Q mutant aggregation. Our study provides mechanistic insights into the R248Q mutation-enhanced aggregation propensity and RSV's potent inhibition against R248Q p53C aggregation.

The inner capsid protein of rotavirus, VP6, emerges as a promising candidate for next-generation vaccines against rotaviruses owing to its abundance in virion particles and high conservation. However, the formation of inclusion bodies during prokaryotic VP6 expression poses a significant hurdle to rotavirus research and applications. Here, we employed experimental and computational approaches to investigate inclusion body formation and aggregation-prone regions (APRs). Heterologous recombinant VP6 expression in Escherichia coli BL21(DE3) cells resulted in inclusion body formation, confirmed by transmission electron microscopy revealing amorphous aggregates. Thioflavin T assay demonstrated incubation temperature-dependent aggregation of VP6 inclusion bodies. Computational predictions of APRs in rotavirus A VP6 protein were performed using sequence-based tools (TANGO, AGGRESCAN, Zyggregator, Waltz, FoldAmyloid, ANuPP, Camsol intrinsic) and structure-based tools (SolubiS, CamSol structurally corrected, Aggrescan3D). A total of 24 consensus APRs were identified, with 21 of them being surface-exposed in VP6. All identified APRs display a predominance of hydrophobic amino acids, ranging from 33 to 100%. Computational identification of these APRs corroborates our experimental observation of VP6 inclusion body or aggregate formation. Characterization of VP6's aggregation propensity facilitates understanding of its behaviour during prokaryotic expression and opens avenues for protein engineering of soluble variants, advancing research on rotavirus VP6 in pathology, therapy, and diagnostics.

Aggregation is an ancient threat that must be overcome by proteins from all organisms to maintain their native functional states. This is essential for the maintenance of metabolic flux and viability of their cellular machineries. Here, we compare the aggregation-resistance strategies adapted by the thermophilic proteins and their mesophilic homologs using a dataset of 373 protein families. Like their mesophilic homologs, the thermophilic protein sequences also contain potential aggregation prone regions (APRs), capable of forming cross-β motif and amyloid-like fibrils. Tetrapeptide and hexapeptide amyloid-like fibril forming sequence patterns and experimentally proven amyloid-like fibril forming peptide sequences were also detected in the thermophilic proteins. Both the thermophilic and the mesophilic proteins use similar strategies to resist aggregation. However, the thermophilic proteins show superior utilization of these strategies. The thermophilic protein monomers show greater ability to "stow away" the APRs in the hydrophobic cores to protect them from solvent exposure. The thermophilic proteins are also better at gatekeeping the APRs by surrounding them with charged residues (Asp, Glu, Lys, and Arg) and Pro to a greater extent. While thermophilic and mesophilic proteins in our dataset are highly homologous and show strong overall sequence conservation, the APRs are not conserved between the homologs. These findings indicate that evolution is working to avoid amyloidogenic regions in proteins. Our results are also consistent with the observation that thermophilic cells often accumulate small molecule osmolytes capable of stabilizing their proteins and other macromolecules. This study has important implications for rational design and formulation of therapeutic proteins and antibodies.

The aggregation of therapeutic antibodies is a major issue for the pharmaceutical industry leading to loss of drug quality, increased dosage, and unwanted immune responses such as the production of anti-drug antibodies (ADA). As aggregation can occur at various stages of development and storage, much work has been performed to reduce or eliminate it. In this report we analyzed four antibodies available in the PDB (1IGT, 1IGY, 1HZH, and 5DK3) using the online software UCSF Chimera to study the structural features of the proteins and the associated N-linked glycans in the CH2 domains of the Fc region. To study antibody aggregation in silico we used the online software TANGO and AGGRESCAN to identify aggregation prone regions (APR) in the antibodies and the influence of the Fc glycans on hydrophobic and aromatic residues present in the APRs. In the 3D structures of 1IGT and 1IGY the glycan chains are in close enough proximity to influence and protect these hydrophobic regions. However, in the 3D structures of 1HZH and 5DK3 the glycans do not appear to influence the likely APRs of the antibodies. Therefore, in these structures we modified the Fc glycan regions by adjusting the glycosylated asparagine side chains and glycosidic bonds. We successfully adjusted the glycan chains of 1HZH and 5DK3 and reduced the distance between them and the APRs to show potential influence on aggregation. However, similar to 5DK3, the influence of glycosylation on the APRs of the antibody was limited due to the size of the glycans present in the 3D structure. This report is based on in silico studies to show how antibody glycans can influence aggregation.

Relationship Between Potential Aggregation-Prone Regions and HLA-DR-Binding T-Cell Immune Epitopes: Implications for Rational Design of Novel and Follow-on Therapeutic Antibodies

Recombinant proteins play pivotal roles in numerous applications including industrial biocatalysts or therapeutics. Despite the recent progress in computational protein structure prediction, protein solubility and reduced aggregation propensity remain challenging attributes to design. Identification of aggregation-prone regions is essential for understanding misfolding diseases or designing efficient protein-based technologies, and as such has a great socio-economic impact. Here, we introduce AggreProt, a user-friendly webserver that automatically exploits an ensemble of deep neural networks to predict aggregation-prone regions (APRs) in protein sequences. Trained on experimentally evaluated hexapeptides, AggreProt compares to or outperforms state-of-the-art algorithms on two independent benchmark datasets. The server provides per-residue aggregation profiles along with information on solvent accessibility and transmembrane propensity within an intuitive interface with interactive sequence and structure viewers for comprehensive analysis. We demonstrate AggreProt efficacy in predicting differential aggregation behaviours in proteins on several use cases, which emphasize its potential for guiding protein engineering strategies towards decreased aggregation propensity and improved solubility. The webserver is freely available and accessible at https://loschmidt.chemi.muni.cz/aggreprot/.

Aggregation of initially stably structured proteins is involved in more than 20 human amyloid diseases. Despite intense research, however, how this class of proteins assembles into amyloid fibrils remains poorly understood, principally because of the complex effects of amino acid substitutions on protein stability, solubility, and aggregation propensity. We address this question using β2-microglobulin (β2m) as a model system, focusing on D76N-β2m that is involved in hereditary amyloidosis. This amino acid substitution causes the aggregation-resilient wild-type protein to become highly aggregation prone in vitro, although the mechanism by which this occurs remained elusive. Here, we identify the residues key to protecting β2m from aggregation by coupling aggregation with antibiotic resistance in E. coli using a tripartite β-lactamase assay (TPBLA). By performing saturation mutagenesis at three different sites (D53X-, D76X-, and D98X-β2m) we show that residue 76 has a unique ability to drive β2m aggregation in vivo and in vitro. Using a randomly mutated D76N-β2m variant library, we show that all of the mutations found to improve protein behavior involve residues in a single aggregation-prone region (APR) (residues 60 to 66). Surprisingly, no correlation was found between protein stability and protein aggregation rate or yield, with several mutations in the APR decreasing aggregation without affecting stability. Together, the results demonstrate the power of the TPBLA to develop proteins that are resilient to aggregation and suggest a model for D76N-β2m aggregation involving the formation of long-range couplings between the APR and Asn76 in a nonnative state.

The recombinant human keratinocyte growth factor (rhKGF) is a highly aggregation-prone therapeutic protein. The high aggregation liability of rhKGF is manifested by loss of the monomeric state, and accumulation of the aggregated species even at moderate temperatures. Here, we analyzed the rhKGF for its vulnerability toward aggregation by detection of aggregation-prone regions (APRs) using several sequence-based computational tools including TANGO, ZipperDB, AGGRESCAN, Zyggregator, Camsol, PASTA, SALSA, WALTZ, SODA, Amylpred, AMYPDB, and structure-based tools including SolubiS, CamSol structurally corrected, Aggrescan3D and spatial aggregation propensity (SAP) algorithm. The sequence-based prediction of APRs in rhKGF indicated that they are mainly located at positions 10–30, 40–60, 61–66, 88–120, and 130–140. Mapping on the rhKGF structure revealed that most of these residues including F16-R25, I43, E45, R47-I56, F61, Y62, N66, L88-E91, E108-F110, A112, N114, T131, and H133-T140 are surface-exposed in the native state which can promote aggregation without major unfolding event, or the conformational change may occur in the oligomers. The other regions are buried in the native state and their contribution to non-native aggregation is mediated by a preceding unfolding event. The structure-based prediction of APRs using the SAP tool limited the number of identified APRs to the dynamically-exposed hydrophobic residues including V12, A50, V51, L88, I89, L90, I118, L135, and I139 mediating the native-state aggregation. Our analysis of APRs in rhKGF identified the regions determining the intrinsic aggregation propensity of the rhKGF which are the candidate positions for engineering the rhKGF to reduce its aggregation tendency. Communicated by Ramaswamy H. Sarma

Owing to its association with a diverse range of human diseases, the determinants of protein aggregation are studied intensively. It is generally accepted that the effective aggregation tendency of a protein depends on many factors such as folding efficiency towards the native state, thermodynamic stability of that conformation, intrinsic aggregation propensity of the polypeptide sequence and its ability to be recognized by the protein quality control system. The intrinsic aggregation propensity of a polypeptide sequence is related to the presence of short APRs (aggregation-prone regions) that self-associate to form intermolecular β-structured assemblies. These are typically short sequence segments (5-15 amino acids) that display high hydrophobicity, low net charge and a high tendency to form β-structures. As the presence of such APRs is a prerequisite for aggregation, a plethora of methods have been developed to identify APRs in amino acid sequences. In the present chapter, the methodological basis of these approaches is discussed, as well as some practical applications.

ANuPP: A Versatile Tool to Predict Aggregation Nucleating Regions in Peptides and Proteins