Aggregation risk prediction for antibodies and its application to biotherapeutic development

Aggregation is a common problem affecting biopharmaceutical development that can have a significant effect on the quality of the product, as well as the safety to patients, particularly because of the increased risk of immune reactions. Here, we describe a new high-throughput screening algorithm developed to classify antibody molecules based on their propensity to aggregate. The tool, constructed and validated on experimental aggregation data for over 500 antibodies, is able to discern molecules with a high aggregation propensity as defined by experimental criteria relevant to bioprocessing and manufacturing of these molecules. Furthermore, we show how this tool can be combined with other computational approaches during early drug development to select molecules with reduced risk of aggregation and optimal developability properties.     

The competition between protein folding and aggregation: off-lattice minimalist model studies

Protein aggregation has been associated with a number of human diseases, and is a serious problem in the manufacture of recombinant proteins. Of particular interest to the biotechnology industry is deleterious aggregation that occurs during the refolding of proteins from inclusion bodies. As a complement to experimental efforts, computer simulations of multi-chain systems have emerged as a powerful tool to investigate the competition between folding and aggregation. Here we report results from Langevin dynamics simulations of minimalist model proteins. Order parameters are developed to follow both folding and aggregation. By mapping natural units to real units, the simulations are shown to be carried out under experimentally relevant conditions. Data pertaining to the contacts formed during the association process show that multiple mechanisms for aggregation exist, but certain pathways are statistically preferred. Kinetic data show that there are multiple time scales for aggregation, although most association events take place at times much shorter than those required for folding. Last, we discuss results presented here as a basis for future work aimed at rational design of mutations to reduce aggregation propensity, as well as for development of small-molecular weight refolding enhancers.     

Manipulating the aggregation activity of human prion-like proteins

Considerable advances in understanding the protein features favoring prion formation in yeast have facilitated the development of effective yeast prion prediction algorithms. Here we discuss a recent study in which we systematically explored the utility of the yeast prion prediction algorithm PAPA for designing mutations to modulate the aggregation activity of the human prion-like protein hnRNPA2B1. Mutations in hnRNPA2B1 cause multisystem proteinopathy in humans, and accelerate aggregation of the protein in vitro. Additionally, mutant hnRNPA2B1 forms cytoplasmic inclusions when expressed in Drosophila, and the mutant prion-like domain can substitute for a portion of a yeast prion domain in supporting prion activity in yeast. PAPA was quite successful at predicting the effects of PrLD mutations on prion activity in yeast and on in vitro aggregation propensity. Additionally, PAPA successfully predicted the effects of most, but not all, mutations in the PrLD of the hnRNPA2B1 protein when expressed in Drosophila. These results suggest that PAPA is quite effective at predicting the effects of mutations on intrinsic aggregation propensity, but that intracellular factors can influence aggregation and prion-like activity in vivo. A more complete understanding of these intracellular factors may inform the next generation of prion prediction algorithms.     

A molecular dynamics study of acylphosphatase in aggregation-promoting conditions: the influence of trifluoroethanol/water solvent

The 98-residue protein acylphosphatase exhibits a high propensity for aggregation under certain conditions. Aggregates formed from wild-type acylphosphatase in the presence of 2,2,2-trifluoroethanol and from highly destabilized mutants are essentially identical in structure. Furthermore, it has been shown by mutational studies that different regions of the protein are important for aggregation and folding. In the present molecular dynamics study, we compare the behavior of the protein in aqueous solution and in a 25% (v/v) 2,2,2-trifluoroethanol/water environment mimicking the experimental conditions. The 2,2,2-trifluoroethanol surrounding affects the structure of the protein mostly in the regions important for aggregation, in good agreement with experimental data. This suggests that the early step of (partly) unfolding, which precedes the aggregation process, has been observed.     

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.     

Sequence Determinants of the Aggregation of Proteins Within Condensates Generated by Liquid-liquid Phase Separation

The transition between the native and amyloid states of proteins can proceed via a deposition pathway via oligomeric intermediates or via a condensation pathway involving liquid droplet intermediates generated through liquid-liquid phase separation. While several computational methods are available to perform sequence-based predictions of the propensity of proteins to aggregate via the deposition pathway, much less is known about the physico-chemical principles that underlie aggregation within condensates. Here we investigate the sequence determinants of aggregation via the condensation pathway, and identify three relevant features: droplet-promoting propensity, aggregation-promoting propensity and multimodal interactions quantified by the binding mode entropy. By using this approach, we show that it is possible to predict aggregation-promoting mutations in droplet-forming proteins associated with amyotrophic lateral sclerosis (ALS). This analysis provides insights into the amino acid code for the conversion of proteins between liquid-like and solid-like condensates.     

SynAggreg: A Multifunctional High-Throughput Technology for Precision Study of Amyloid Aggregation and Systematic Discovery of Synergistic Inhibitor Compounds

Numerous proteins can coalesce into amyloid self-assemblies, which are responsible for a class of diseases called amyloidoses, but which can also fulfill important biological functions and are of great interest for biotechnology. Amyloid aggregation is a complex multi-step process, poorly prone to detailed structural studies. Therefore, small molecules interacting with amyloids are often used as tools to probe the amyloid aggregation pathway and in some cases to treat amyloidoses as they prevent pathogenic protein aggregation. Here, we report on SynAggreg, an in vitro high-throughput (HT) platform dedicated to the precision study of amyloid aggregation and the effect of modulator compounds. SynAggreg relies on an accurate bi-fluorescent amyloid-tracer readout that overcomes some limitations of existing HT methods. It allows addressing diverse aspects of aggregation modulation that are critical for pathomechanistic studies, such as the specificity of compounds toward various amyloids and their effects on aggregation kinetics, as well as the co-assembly propensity of distinct amyloids and the influence of prion-like seeding on self-assembly. Furthermore, SynAggreg is the first HT technology that integrates tailored methodology to systematically identify synergistic compound combinations—an emerging strategy to improve fatal amyloidoses by targeting multiple steps of the aggregation pathway. To this end, we apply analytical combinatorial scores to rank the inhibition efficiency of couples of compounds and to readily detect synergism. Finally, the SynAggreg platform should be suited for the characterization of a broad class of amyloids, whether of interest for drug development purposes, for fundamental research on amyloid functions, or for biotechnological applications.     

A comparative study of the relationship between protein structure and beta-aggregation in globular and intrinsically disordered proteins

A growing number of proteins are being identified that are biologically active though intrinsically disordered, in sharp contrast with the classic notion that proteins require a well-defined globular structure in order to be functional. At the same time recent work showed that aggregation and amyloidosis are initiated in amino acid sequences that have specific physico-chemical properties in terms of secondary structure propensities, hydrophobicity and charge. In intrinsically disordered proteins (IDPs) such sequences would be almost exclusively solvent-exposed and therefore cause serious solubility problems. Further, some IDPs such as the human prion protein, synuclein and Tau protein are related to major protein conformational diseases. However, this scenario contrasts with the large number of unstructured proteins identified, especially in higher eukaryotes, and the fact that the solubility of these proteins is often particularly good. We have used the algorithm TANGO to compare the beta aggregation tendency of a set of globular proteins derived from SCOP and a set of 296 experimentally verified, non-redundant IDPs but also with a set of IDPs predicted by the algorithms DisEMBL and GlobPlot. Our analysis shows that the beta-aggregation propensity of all-alpha, all-beta and mixed alpha/beta globular proteins as well as membrane-associated proteins is fairly similar. This illustrates firstly that globular structures possess an appreciable amount of structural frustration and secondly that beta-aggregation is not determined by hydrophobicity and beta-sheet propensity alone. We also show that globular proteins contain almost three times as much aggregation nucleating regions as IDPs and that the formation of highly structured globular proteins comes at the cost of a higher beta-aggregation propensity because both structure and aggregation obey very similar physico-chemical constraints. Finally, we discuss the fact that although IDPs have a much lower aggregation propensity than globular proteins, this does not necessarily mean that they have a lower potential for amyloidosis.     

A screening tool for therapeutic monoclonal antibodies: Identifying the most stable protein and its best formulation based on thioflavin T binding

The lack of a fast selection method to identify the most stable protein is one of the major challenges for developing successful therapeutic protein formulations more rapidly. The swift and accurate detection of small amounts of aggregates is another problem since aggregates may trigger an immunological response and the aggregation decreases the biological activity of the antibody. Here we present an alternative method for initial screening of the aggregation propensity of proteins, using monoclonal antibodies (mAb) as an example and thioflavin T (ThT) binding. The major advantage of ThT binding is the short duration of testing compared with size-exclusion chromatography (SEC) measurements that can take 6 months or more even under accelerated conditions. The tendency to aggregate of each therapeutic human mAb probed with the ThT assay, together with SEC, is employed to formulate the ranking of mAb aggregation. ThT binding can determine the propensity of proteins to aggregate in a few days, illustrating that ThT binding would be a valuable screening tool.     

Pathways and intermediates of amyloid fibril formation

The lack of understanding of amyloid fibril formation at the molecular level is a major obstacle in devising strategies to interfere with the pathologies linked to peptide or protein aggregation. In particular, little is known on the role of intermediates and fibril elongation pathways as well as their dependence on the intrinsic tendency of a polypeptide chain to self-assembly by β-sheet formation (β-aggregation propensity). Here, coarse-grained simulations of an amphipathic polypeptide show that a decrease in the β-aggregation propensity results in a larger heterogeneity of elongation pathways, despite the essentially identical structure of the final fibril. Protofibrillar intermediates that are thinner, shorter and less structured than the final fibril accumulate along some of these pathways. Moreover, the templated formation of an additional protofilament on the lateral surface of a protofibril is sometimes observed as a collective transition. Conversely, for a polypeptide model with a high β-aggregation propensity, elongation proceeds without protofibrillar intermediates. Therefore, changes in intrinsic β-aggregation propensity modulate the relative accessibility of parallel routes of aggregation.     

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.     

Physical degradation of proteins in well-defined fluid flows studied within a four-roll apparatus

In most applications of biotechnology and downstream processing proteins are exposed to fluid stresses in various flow configurations which often lead to the formation of unwanted protein aggregates. In this paper we present physical degradation experiments for proteins under well-defined flow conditions in a four-roll apparatus. The flow field was characterized numerically by computational fluid dynamics (CFD) and experimentally by particle image velocimetry (PIV). The local shear strain rate as well as the local shear and elongation rate was used to characterize the hydrodynamic stress environment acting on the proteins. Lysozyme was used as a model protein and subjected to well-defined fluid stresses in high and low stress environment. By using in situ turbidity measurements during stressing the aggregate formation was monitored directly in the fluid flow. An increase in absorbance at 350 nm was attributed to a higher content of visible particles (>1 μm). In addition to lysozyme, the formation of aggregates was confirmed for two larger proteins (bovine serum albumin and alcohol dehydrogenase). Thus, the presented experimental setup is a helpful tool to monitor flow-induced protein aggregation with high reproducibility. For instance, screening experiments for formulation development of biopharmaceuticals for fill and finish operations can be performed in the lab-scale in a short time-period if the stress distributions in the application are transferred and applied in the four-roll mill.     

Mammalian antimicrobial peptide protegrin‐4 self assembles and forms amyloid‐like aggregates: Assessment of its functional relevance

Protegrin‐4 (PG‐4) is a member of the porcine leukocyte protegrins family of cysteine‐rich antimicrobial peptides (AMPs) isolated from Sus scrofa. It consists of 18 amino acid residues and works as a part of innate immune system. In this study, we examined the intrinsic aggregation propensity of this AMP using multiple computational algorithms, namely, TANGO, AGGRESCAN, FOLDAMYLOID, AMYLPRED, and ZYGGREGATOR, and found that the peptide is predicted to have a high propensity for the β sheet formation that disposes this peptide to be amyloidogenic. Under in vitro conditions, PG‐4 formed visible aggregates and displayed the hallmark properties of typical amyloids such as enhanced binding of Congo red, increased fluorescence with Thioflavin‐T, and fibrillar morphology under transmission electron microscopy. Then we examined its antimicrobial activity against Bacillus subtilis and found that the aggregated peptide retained its antimicrobial activity. Additionally, the aggregates remain non‐toxic to the HEK293 and Caco2 cells. Our study suggests that the inherent aggregation properties of AMP can rationally be explored as a potential source of peptide‐based antimicrobials with enhanced stability.     

Computational investigation on the effects of H50Q and G51D mutations on the α-Synuclein aggregation propensity

The aggregation of α-synuclein is linked directly to the histopathology of Parkinson’s disease (PD). However, several missense mutations present in the α-synuclein gene (SNCA) have been known to be associated with PD. Several studies have highlighted the effect of SNCA mutations on the α-synuclein aggregation, but their pathological roles are not completely established. In this study, we have focused on the effects of the recently discovered α-synuclein missense mutants (H50Q and G51D) on the aggregation using computational approaches. We performed all atom molecular dynamics (MD) simulation on these mutants and compared their conformational dynamics with Wild-Type (WT) α-synuclein. We noticed the solvent accessible surface area (SASA), radius of gyration, atomic fluctuations, and beta strand content to be higher in H50Q than G51D and WT. Using PDBSum online server; we analyzed the inter-molecular interactions that drive the association of monomeric units of H50Q, WT, and G51D in forming the respective homo-dimer. We noticed the interface area, number of interacting residues and binding free energy to be higher for H50Q homo-dimer than the WT and G51D homo-dimers. Our findings in this study suggest that in comparison to WT and G51D, H50Q mutation to have a positive effect on increasing the α-synuclein aggregation propensity. Hence, we see that H50Q and G51D mutation show conflicting effect on the aggregation propensity of α-synuclein.     

Mechanism for the alpha-helix to beta-hairpin transition

The aggregation of alpha-helix-rich proteins into beta-sheet-rich amyloid fibrils is associated with fatal diseases, such as Alzheimer's disease and prion disease. During an aggregation process, protein secondary structure elements-alpha-helices-undergo conformational changes to beta-sheets. The fact that proteins with different sequences and structures undergo a similar transition on aggregation suggests that the sequence nonspecific hydrogen bond interaction among protein backbones is an important factor. We perform molecular dynamics simulations of a polyalanine model, which is an alpha-helix in its native state and observe a metastable beta-hairpin intermediate. Although a beta-hairpin has larger potential energy than an alpha-helix, the entropy of a beta-hairpin is larger because of fewer constraints imposed by the hydrogen bonds. In the vicinity of the transition temperature, we observe the interconversion of the alpha-helix and beta-sheet states via a random coil state. We also study the effect of the environment by varying the relative strength of side-chain interactions for a designed peptide-an alpha-helix in its native state. For a certain range of side-chain interaction strengths, we find that the intermediate beta-hairpin state is destabilized and even disappears, suggesting an important role of the environment in the aggregation propensity of a peptide.     

Protein Aggregation Profile of the Bacterial Cytosol

Background

Protein misfolding is usually deleterious for the cell, either as a consequence of the loss of protein function or the buildup of insoluble and toxic aggregates. The aggregation behavior of a given polypeptide is strongly influenced by the intrinsic properties encoded in its sequence. This has allowed the development of effective computational methods to predict protein aggregation propensity.

Methodology/Principal Findings

Here, we use the AGGRESCAN algorithm to approximate the aggregation profile of an experimental cytosolic Escherichia coli proteome. The analysis indicates that the aggregation propensity of bacterial proteins is associated with their length, conformation, location, function, and abundance. The data are consistent with the predictions of other algorithms on different theoretical proteomes.

Conclusions/Significance

Overall, the study suggests that the avoidance of protein aggregation in functional environments acts as a strong evolutionary constraint on polypeptide sequences in both prokaryotic and eukaryotic organisms.     

Recruitment into stress granules prevents irreversible aggregation of FUS protein mislocalized to the cytoplasm

Fused in sarcoma (FUS) belongs to the group of RNA-binding proteins implicated as underlying factors in amyotrophic lateral sclerosis (ALS) and certain other neurodegenerative diseases. Multiple FUS gene mutations have been linked to hereditary forms, and aggregation of FUS protein is believed to play an important role in pathogenesis of these diseases. In cultured cells, FUS variants with disease-associated amino acid substitutions or short deletions affecting nuclear localization signal (NLS) and causing cytoplasmic mislocalization can be sequestered into stress granules (SGs). We demonstrated that disruption of motifs responsible for RNA recognition and binding not only prevents SG recruitment, but also dramatically increases the protein propensity to aggregate in the cell cytoplasm with formation of juxtanuclear structures displaying typical features of aggresomes. Functional RNA-binding domains from TAR DNA-binding protein of 43 kDa (TDP-43) fused to highly aggregation-prone C-terminally truncated FUS protein restored the ability to enter SGs and prevented aggregation of the chimeric protein. Truncated FUS was also able to trap endogenous FUS molecules in the cytoplasmic aggregates. Our data indicate that RNA binding and recruitment to SGs protect cytoplasmic FUS from aggregation, and loss of this protection may trigger its pathological aggregation in vivo.     

Using extensional flow to reveal diverse aggregation landscapes for three IgG1 molecules

Monoclonal antibodies (mAbs) currently dominate the biopharmaceutical sector due to their potency and efficacy against a range of disease targets. These proteinaceous therapeutics are, however, susceptible to unfolding, mis‐folding, and aggregation by environmental perturbations. Aggregation thus poses an enormous challenge to biopharmaceutical development, production, formulation, and storage. Hydrodynamic forces have also been linked to aggregation, but the ability of different flow fields (e.g., shear and extensional flow) to trigger aggregation has remained unclear. To address this question, we previously developed a device that allows the degree of extensional flow to be controlled. Using this device we demonstrated that mAbs are particularly sensitive to the force exerted as a result of this flow‐field. Here, to investigate the utility of this device to bio‐process/biopharmaceutical development, we quantify the effects of the flow field and protein concentration on the aggregation of three mAbs. We show that the response surface of mAbs is distinct from that of bovine serum albumin (BSA) and also that mAbs of similar sequence display diverse sensitivity to hydrodynamic flow. Finally, we show that flow‐induced aggregation of each mAb is ameliorated by different buffers, opening up the possibility of using the device as a formulation tool. Perturbation of the native state by extensional flow may thus allow identification of aggregation‐resistant mAb candidates, their bio‐process parameters and formulation to be optimized earlier in the drug‐discovery pipeline using sub‐milligram quantities of material.     

Pseudophosphorylation of tau protein alters its ability for self-aggregation

Filamentous tau protein deposits are a pathological hallmark of a group of neurodegenerative disorders (tauopathies). Tau protein in these aggregates is highly phosphorylated at different phosphorylation sites. Although tau filaments can be formed by heparin-induced aggregation of unphosphorylated recombinant tau, it is not known how tau phosphorylation modulates aggregation behaviour. Analysis of the effect of tau phosphorylation at defined single or multiple sites is hampered by the low specificity of protein kinases and the highly dynamic turnover of phosphorylation in vivo. To overcome this problem we employed site-directed mutagenesis to convert serine and threonine to aspartic acid or glutamic acid, which introduce a negative charge and conformational change that mimic phosphorylation. We tested 14 different mutated tau proteins for their propensity for self-aggregation and formation of tau filaments. Tau aggregation was monitored with thioflavin S fluorescence in the presence of different inducers such as heparin, Al3+, Fe2+ and Fe3+. We found that mutations in the N-terminal portion up to amino acid 208 mainly suppress tau aggregation, whereas mutations in the C-terminal region mainly lead to an enhanced aggregation. Mutations in the middle portion of tau showed a mixed picture of suppression and enhancement of aggregation. A single amino acid change Ser422Glu has aggregation-favouring properties with all four inducers.     

Testosterone disrupts human collaboration by increasing egocentric choices

Collaboration can provide benefits to the individual and the group across a variety of contexts. Even in simple perceptual tasks, the aggregation of individuals' personal information can enable enhanced group decision-making. However, in certain circumstances such collaboration can worsen performance, or even expose an individual to exploitation in economic tasks, and therefore a balance needs to be struck between a collaborative and a more egocentric disposition. Neurohumoral agents such as oxytocin are known to promote collaborative behaviours in economic tasks, but whether there are opponent agents, and whether these might even affect information aggregation without an economic component, is unknown. Here, we show that an androgen hormone, testosterone, acts as such an agent. Testosterone causally disrupted collaborative decision-making in a perceptual decision task, markedly reducing performance benefit individuals accrued from collaboration while leaving individual decision-making ability unaffected. This effect emerged because testosterone engendered more egocentric choices, manifest in an overweighting of one's own relative to others' judgements during joint decision-making. Our findings show that the biological control of social behaviour is dynamically regulated not only by modulators promoting, but also by those diminishing a propensity to collaborate.