Aggregation-Prone Motifs in Human Immunoglobulin G

Therapeutic antibodies of many different IgG subclasses (IgG1, IgG2 and IgG4) are used in the treatment of various cancers, rheumatoid arthritis and other inflammatory and infectious diseases. These antibodies are stored for long durations under high concentrations as required in the disease treatment. Unfortunately, these antibodies aggregate under these storage conditions, leading to a decrease in antibody activity and raising concerns about causing an immunological response. Thus, there is a tremendous need to identify the aggregation-prone regions in different classes of antibodies. We use the SAP (spatial-aggregation-propensity) technology based on molecular simulations to determine the aggregation-prone motifs in the constant regions of IgG1 classes of antibodies. Mutations engineered on these aggregation-prone motif regions led to antibodies of enhanced stability. Fourteen aggregation-prone motifs are identified, with each motif containing one to seven residues. While some of these motifs contain residues that are neighbors in primary sequence, others contain residues that are far apart in primary sequence but are close together in the tertiary structure. Comparison of the IgG1 sequence with those of other subclasses (IgG2, IgG3 and IgG4) showed that these aggregation-prone motifs are largely preserved among all IgG subclasses. Other broader classes of antibodies (IgA1, IgD, IgE and IgM), however, differed in these motif regions. The aggregation-prone motifs identified were therefore common to all IgG subclasses, but differ from those of non-IgG classes. Moreover, since the motifs identified are in the constant regions, they are applicable for all antibodies within the IgG class irrespective of the variable region. Thus, the motif regions identified could be modified on all IgGs to yield antibodies of enhanced stability.     

Prediction of protein binding regions

Identifying protein binding sites provides important clues to the function of a protein. Experimental methods to identify the binding sites such as determining the crystal structures of protein complexes are extremely laborious and expensive. Here, we present a computational technique called spatial aggregation propensity (SAP) based on molecular simulations to predict protein binding sites. We apply this technique to two model proteins, an IgG1 antibody and epidermal growth factor receptor (EGFR) and demonstrate that SAP predicts protein binding regions with very good accuracy. In the case of the IgG1 antibody, SAP accurately predicts binding regions with the Fc-receptor, protein-A, and protein-G. For EGFR, SAP accurately predicts binding regions with EGF, TGFα, and with another EGFR. The resolution of SAP is varied to obtain a detailed picture of these binding sites. We also show that some of these binding sites overlap with protein self-aggregation prone regions. We demonstrate how SAP analysis can be used to engineer the protein to remove unfavorable aggregation prone regions without disturbing protein binding regions. The SAP technique could be also used to predict the yet unknown binding sites of numerous proteins, thereby providing clues to their function.     

Potential aggregation prone regions in biotherapeutics

Aggregation of a biotherapeutic is of significant concern and judicious process and formulation development is required to minimize aggregate levels in the final product. Aggregation of a protein in solution is driven by intrinsic and extrinsic factors. In this work we have focused on aggregation as an intrinsic property of the molecule. We have studied the sequences and Fab structures of commercial and non-commercial antibody sequences for their vulnerability towards aggregation by using sequence based computational tools to identify potential aggregation-prone motifs or regions. The mAbs in our dataset contain 2 to 8 aggregation-prone motifs per heavy and light chain pair. Some of these motifs are located in variable domains, primarily in CDRs. Most aggregation-prone motifs are rich in β branched aliphatic and aromatic residues. Hydroxyl-containing Ser/Thr residues are also found in several aggregation-prone motifs while charged residues are rare. The motifs found in light chain CDR3 are glutamine (Q)/asparagine (N) rich. These motifs are similar to the reported aggregation promoting regions found in prion and amyloidogenic proteins that are also rich in Q/N, aliphatic and aromatic residues. The implication is that one possible mechanism for aggregation of mAbs may be through formation of cross-β structures and fibrils. Mapping on the available Fab – receptor/antigen complex structures reveals that these motifs in CDRs might also contribute significantly towards receptor/antigen binding. Our analysis identifies the opportunity and tools for simultaneous optimization of the therapeutic protein sequence for potency and specificity while reducing vulnerability towards aggregation.     

Potential aggregation prone regions in biotherapeutics: A survey of commercial monoclonal antibodies

Aggregation of a biotherapeutic is of significant concern and judicious process and formulation development is required to minimize aggregate levels in the final product. Aggregation of a protein in solution is driven by intrinsic and extrinsic factors. In this work we have focused on aggregation as an intrinsic property of the molecule. We have studied the sequences and Fab structures of commercial and non-commercial antibody sequences for their vulnerability towards aggregation by using sequence based computational tools to identify potential aggregation-prone motifs or regions. The mAbs in our dataset contain 2 to 8 aggregation-prone motifs per heavy and light chain pair. Some of these motifs are located in variable domains, primarily in CDRs. Most aggregation-prone motifs are rich in β branched aliphatic and aromatic residues. Hydroxyl-containing Ser/Thr residues are also found in several aggregation-prone motifs while charged residues are rare. The motifs found in light chain CDR3 are glutamine (Q)/asparagine (N) rich. These motifs are similar to the reported aggregation promoting regions found in prion and amyloidogenic proteins that are also rich in Q/N, aliphatic and aromatic residues. The implication is that one possible mechanism for aggregation of mAbs may be through formation of cross-β structures and fibrils. Mapping on the available Fab—receptor/antigen complex structures reveals that these motifs in CDRs might also contribute significantly towards receptor/antigen binding. Our analysis identifies the opportunity and tools for simultaneous optimization of the therapeutic protein sequence for potency and specificity while reducing vulnerability towards aggregation.Key words: monoclonal antibody, aggregation, antibody sequence, aggregation-prone region, aggregation prediction     

Intrinsically Disordered and Aggregation Prone Regions Underlie β-Aggregation in S100 Proteins

S100 proteins are small dimeric calcium-binding proteins which control cell cycle, growth and differentiation via interactions with different target proteins. Intrinsic disorder is a hallmark among many signaling proteins and S100 proteins have been proposed to contain disorder-prone regions. Interestingly, some S100 proteins also form amyloids: S100A8/A9 forms fibrils in prostatic inclusions and S100A6 fibrillates in vitro and seeds SOD1 aggregation. Here we report a study designed to investigate whether β-aggregation is a feature extensive to more members of S100 family. In silico analysis of seven human S100 proteins revealed a direct correlation between aggregation and intrinsic disorder propensity scores, suggesting a relationship between these two independent properties. Averaged position-specific analysis and structural mapping showed that disorder-prone segments are contiguous to aggregation-prone regions and that whereas disorder is prominent on the hinge and target protein-interaction regions, segments with high aggregation propensity are found in ordered regions within the dimer interface. Acidic conditions likely destabilize the seven S100 studied by decreasing the shielding of aggregation-prone regions afforded by the quaternary structure. In agreement with the in silico analysis, hydrophobic moieties become accessible as indicated by strong ANS fluorescence. ATR-FTIR spectra support a structural inter-conversion from α-helices to intermolecular β-sheets, and prompt ThT-binding takes place with no noticeable lag phase. Dot blot analysis using amyloid conformational antibodies denotes a high diversity of conformers; subsequent analysis by TEM shows fibrils as dominant species. Altogether, our data suggests that β-aggregation and disorder-propensity are related properties in S100 proteins, and that the onset of aggregation is likely triggered by loss of protective tertiary and quaternary interactions.     

Amyloidogenic Regions and Interaction Surfaces Overlap in Globular Proteins Related to Conformational Diseases

Protein aggregation underlies a wide range of human disorders. The polypeptides involved in these pathologies might be intrinsically unstructured or display a defined 3D-structure. Little is known about how globular proteins aggregate into toxic assemblies under physiological conditions, where they display an initially folded conformation. Protein aggregation is, however, always initiated by the establishment of anomalous protein-protein interactions. Therefore, in the present work, we have explored the extent to which protein interaction surfaces and aggregation-prone regions overlap in globular proteins associated with conformational diseases. Computational analysis of the native complexes formed by these proteins shows that aggregation-prone regions do frequently overlap with protein interfaces. The spatial coincidence of interaction sites and aggregating regions suggests that the formation of functional complexes and the aggregation of their individual subunits might compete in the cell. Accordingly, single mutations affecting complex interface or stability usually result in the formation of toxic aggregates. It is suggested that the stabilization of existing interfaces in multimeric proteins or the formation of new complexes in monomeric polypeptides might become effective strategies to prevent disease-linked aggregation of globular proteins.     

Charge-Rich Regions Modulate the Anti-Aggregation Activity of Hsp90

Protein aggregation can have dramatic effects on cellular function and plays a causative role in many human diseases. In all cells, molecular chaperones bind to aggregation-prone proteins and hinder aggregation. The ability of a protein to resist aggregation and remain soluble in aqueous solution is linked to the physical properties of the protein. Numerous physical studies demonstrate that charged atoms favor solubility. We note that many molecular chaperones possess a substantial negative charge that may allow them to impart solubility on aggregation-prone proteins. Hsp90 is one such negatively charged molecular chaperone. The charge on Hsp90 is largely concentrated in two highly acidic regions. To investigate the relationship between chaperone charge and protein solubility, we deleted these charge-rich regions and analyzed the resulting Hsp90 constructs for anti-aggregation activity. We found that deletion of both charge-rich regions dramatically impaired Hsp90 anti-aggregation activity. The anti-aggregation role of the deleted charge-rich regions could be due to net charge or sequence-specific features. To distinguish these possibilities, we attached an acid-rich region with a distinct amino acid sequence to our double-deleted Hsp90 construct. This charge rescue construct displayed effective anti-aggregation activity indicating that the net charge of Hsp90 contributes to its anti-aggregation activity.     

Yeast Proteins may Reversibly Aggregate like Amphiphilic Molecules

More than a hundred proteins in yeast reversibly aggregate and phase-separate in response to various stressors, such as nutrient depletion and heat shock. We know little about the protein sequence and structural features behind this ability, which has not been characterized on a proteome-wide level. To identify the distinctive features of aggregation-prone protein regions, we apply machine learning algorithms to genome-scale limited proteolysis-mass spectrometry (LiP-MS) data from yeast proteins. LiP-MS data reveals that 96 proteins show significant structural changes upon heat shock. We find that in these proteins the propensity to phase separate cannot be solely driven by disordered regions, because their aggregation-prone regions (APRs) are not significantly disordered. Instead, the phase separation of these proteins requires contributions from both disordered and structured regions. APRs are significantly enriched in aliphatic residues and depleted in positively charged amino acids. Aggregator proteins with longer APRs show a greater propensity to aggregate, a relationship that can be explained by equilibrium statistical thermodynamics. Altogether, our observations suggest that proteome-wide reversible protein aggregation is mediated by sequence-encoded properties. We propose that aggregating proteins resemble supra-molecular amphiphiles, where APRs are the hydrophobic parts, and non-APRs are the hydrophilic parts.     

Protein aggregation and amyloidosis: confusion of the kinds?

Recent years have witnessed major advances in our understanding of the structural basis of protein aggregation on several fronts. Firstly, high-resolution structural information that remained elusive for many years was provided by a series of studies of amyloid fibers using NMR, X-ray crystallography and electron microscopy, thereby confirming earlier models based on lower resolution observations. Secondly, studies of the sequence determinants of protein aggregation culminated in the development of computer algorithms that predict aggregation-prone sequences with good accuracy, allowing the design of mutations that reduce aggregation. Thirdly, based on the first results from such predictions and on statistical analysis of naturally occurring aggregating sequences, a picture is emerging in which aggregation-prone sequences are capped by gatekeeper residues that oppose aggregation. In addition to their aggregation-opposing function, it seems that gatekeeper residues are also important in determining chaperone selectivity for strongly aggregating regions. Finally, recent computational and experimental work shows that preventing aggregation does not necessarily mean that amyloid formation is prevented and vice versa. Thus, although aggregation and amyloidosis correlate to a certain extent, they are different processes and should be treated as such.     

Folding on the chaperone: yield enhancement through loose binding

A variety of small cageless chaperones have been discovered that can assist protein folding without the consumption of ATP. These include mini-chaperones (catalytically active fragments of larger chaperones), as well as small proteins such as alpha-casein and detergents acting as "artificial chaperones." These chaperones all possess exposed hydrophobic patches on their surface that act as recognition sites for misfolded proteins. They lack the complexity of chaperonins (that encapsulate proteins in their inner rings) and their study can offer insight into the minimal requirements for chaperone function. We use molecular dynamics simulations to investigate how a cageless chaperone, modeled as a sphere of tunable hydrophobicity, can assist folding of a substrate protein. We find that under steady-state (non-stress) conditions, cageless chaperones that bind to a single substrate protein increase folding yields by reducing the time the substrate spends in an aggregation-prone state in a dual manner: (a) by competing for aggregation-prone hydrophobic sites on the surface of a protein, hence reducing the time the protein spends unprotected in the bulk and (b) by accelerating folding rates of the protein. In both cases, the chaperone must bind to and hold the protein loosely enough to allow the protein to change its conformation and fold while bound. Loose binding may enable small cageless chaperones to help proteins fold and avoid aggregation under steady-state conditions, even at low concentrations, without the consumption of ATP.     

Aggregation-prone Regions in HYPK Help It to Form Sequestration Complex for Toxic Protein Aggregates

Protein aggregates result from altered structural conformations and they can perturb cellular homeostasis. Prevention mechanisms, which function against protein aggregation by modulatory processes, are diverse and redundant. In this study, we have characterized Huntingtin interacting protein K (HYPK) as a global aggregation-regulatory protein. We report the mechanistic details of how HYPK's aggregation-prone regions allow it to sense and prevent other toxic protein's aggregation by forming unique annular-shaped sequestration complexes. Screenings for interacting partners of different aggregation-prone proteins identify HYPK as a global interacting partner/regulator of Huntingtin97Qexon1, α-Synuclein-A53T and Superoxide dismutase1-G93A. C-terminal hydrophobic region in HYPK makes direct contacts with aggregates to initiate the formation of sequestration complexes. HYPK acts as aggregate sensor by existing in a seeded amyloid-like state which also favors its own concentration-dependent self-oligomerization. Oligomerization of HYPK leads to annular and non-fibrillar/amorphous aggregates. Two hydrophobic segments in the C-terminus of HYPK are responsible for its own aggregations. Self-association of HYPK follows seed nucleation, in which oligomeric HYPK seeds nucleate to annular structures. Annular oligomers of HYPK fuse with each other to form amorphous aggregates. HYPK shows differential interactions with aggregation-prone and non-aggregating proteins, as it preferentially binds to aggregation-prone proteins with higher affinity than native/non-aggregating proteins. This favors the formation of HYPK's sequestration complexes both in cytosol and in ribosome. Besides having aggregation-preventive property, HYPK also reduces the cellular level of toxic proteins. In vivo, HYPK sequestration complexes prevent the formation of toxic protein aggregates to physiologically show positive impact on cell survival and restoration of normal cell physiology.     

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.     

Polyglutamine expansion diseases: More than simple repeats

Polyglutamine (polyQ) repeat-containing proteins are widespread in the human proteome but only nine of them are associated with highly incapacitating neurodegenerative disorders. The genetic expansion of the polyQ tract in disease-related proteins triggers a series of events resulting in neurodegeneration. The polyQ tract plays the leading role in the aggregation mechanism, but other elements modulate the aggregation propensity in the context of the full-length proteins, as implied by variations in the length of the polyQ tract required to trigger the onset of a given polyQ disease. Intrinsic features such as the presence of aggregation-prone regions (APRs) outside the polyQ segments and polyQ-flanking sequences, which synergistically participate in the aggregation process, are emerging for several disease-related proteins. The inherent polymorphic structure of polyQ stretches places the polyQ proteins in a central position in protein–protein interaction networks, where interacting partners may additionally shield APRs or reshape the aggregation course. Expansion of the polyQ tract perturbs the cellular homeostasis and contributes to neuronal failure by modulating protein–protein interactions and enhancing toxic oligomerization. Post-translational modifications further regulate self-assembly either by directly altering the intrinsic aggregation propensity of polyQ proteins, by modulating their interaction with different macromolecules or by modifying their withdrawal by the cell quality control machinery. Here we review the recent data on the multifaceted aggregation pathways of disease-related polyQ proteins, focusing on ataxin-3, the protein mutated in Machado-Joseph disease. Further mechanistic understanding of this network of events is crucial for the development of effective therapies for polyQ diseases.     

Sequestration of the Aβ Peptide Prevents Toxicity and Promotes Degradation In Vivo

Protein aggregation, arising from the failure of the cell to regulate the synthesis or degradation of aggregation-prone proteins, underlies many neurodegenerative disorders. However, the balance between the synthesis, clearance, and assembly of misfolded proteins into neurotoxic aggregates remains poorly understood. Here we study the effects of modulating this balance for the amyloid-beta (Aβ) peptide by using a small engineered binding protein (Z_Aβ3) that binds with nanomolar affinity to Aβ, completely sequestering the aggregation-prone regions of the peptide and preventing its aggregation. Co-expression of Z_Aβ3 in the brains of Drosophila melanogaster expressing either Aβ₄₂ or the aggressive familial associated E22G variant of Aβ₄₂ abolishes their neurotoxic effects. Biochemical analysis indicates that monomer Aβ binding results in degradation of the peptide in vivo. Complementary biophysical studies emphasize the dynamic nature of Aβ aggregation and reveal that Z_Aβ3 not only inhibits the initial association of Aβ monomers into oligomers or fibrils, but also dissociates pre-formed oligomeric aggregates and, although very slowly, amyloid fibrils. Toxic effects of peptide aggregation in vivo can therefore be eliminated by sequestration of hydrophobic regions in monomeric peptides, even when these are extremely aggregation prone. Our studies also underline how a combination of in vivo and in vitro experiments provide mechanistic insight with regard to the relationship between protein aggregation and clearance and show that engineered binding proteins may provide powerful tools with which to address the physiological and pathological consequences of protein aggregation.     

Rational design of therapeutic mAbs against aggregation through protein engineering and incorporation of glycosylation motifs applied to bevacizumab

The aggregation of biotherapeutics is a major hindrance to the development of successful drug candidates; however, the propensity to aggregate is often identified too late in the development phase to permit modification to the protein's sequence. Incorporating rational design for the stability of proteins in early discovery has numerous benefits. We engineered out aggregation-prone regions on the Fab domain of a therapeutic monoclonal antibody, bevacizumab, to rationally design a biobetter drug candidate. With the purpose of stabilizing bevacizumab with respect to aggregation, 2 strategies were undertaken: single point mutations of aggregation-prone residues and engineering a glycosylation site near aggregation-prone residues to mask these residues with a carbohydrate moiety. Both of these approaches lead to comparable decreases in aggregation, with an up to 4-fold reduction in monomer loss. These single mutations and the new glycosylation pattern of the Fab domain do not modify binding to the target. Biobetters with increased stability against aggregation can therefore be generated in a rational manner, by either removing or masking the aggregation-prone region or crowding out protein-protein interactions.     

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

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

Native structure protects SUMO proteins from aggregation into amyloid fibrils

SUMO proteins belong to the Ubiquitin-like protein family, all sharing a common fold and a similar mechanism of conjugation to target polypeptides. SUMO is ubiquitous in all eukaryotes and participates in many crucial pathways. Native SUMO proteins are highly soluble, a property that is exploited in biotechnology. Moreover, SUMO regulates the solubility of aggregation-prone proteins in neurodegenerative disorders. Despite these properties, we show here that human SUMO1, SUMO2, and SUMO3 proteins are at risk of aggregation into amyloid structures if their native conformation is perturbed. Aggregation is mediated by specific regions, which overlap with SUMO functional interfaces, illustrating a competition between function and aggregation. Aggregation of SUMOs might have important physiological implications because disruption of the SUMO pathway is lethal in different organisms. It appears that functional constraints make it difficult to avoid the competition between productive folding and deleterious aggregation in globular proteins, even for essential polypeptides.     

Sequence and structural determinants of Cu, Zn superoxide dismutase aggregation

Diverse point mutations in the enzyme Cu, Zn superoxide dismutase (SOD1) are linked to its aggregation in the familial form of the disease amyotrophic lateral sclerosis. The disease-associated mutations are known to destabilize the protein, but the structural basis of the aggregation of the destabilized protein and the structure of aggregates are not well understood. Here, we investigate in silico the sequence and structural determinants of SOD1 aggregation: (1) We identify sequence fragments in SOD1 that have a high aggregation propensity, using only the sequence of SOD1, and (2) we perform molecular dynamics simulations of the SOD1 dimer folding and misfolding. In both cases, we identify identical regions of the protein as having high propensity to form intermolecular interactions. These regions correspond to the N- and C-termini, and two crossover loops and two beta-strands in the Greek-key native fold of SOD1. Our results suggest that the high aggregation propensity of mutant SOD1 may result from a synergy of two factors: the presence of highly amyloidogenic sequence fragments ("hot spots"), and the presence of these fragments in regions of the protein that are structurally most likely to form intermolecular contacts under destabilizing conditions. Therefore, we postulate that the balance between the self-association of aggregation-prone sequences and the specific structural context of these sequences in the native state determines the aggregation propensity of proteins.