Rapid and accurate in silico solubility screening of a monoclonal antibody library

Stability of Protein Pharmaceuticals: Recent Advances

There have been significant advances in the formulation and stabilization of proteins in the liquid state over the past years since our previous review. Our mechanistic understanding of protein-excipient interactions has increased, allowing one to develop formulations in a more rational fashion. The field has moved towards more complex and challenging formulations, such as high concentration formulations to allow for subcutaneous administration and co-formulation. While much of the published work has focused on mAbs, the principles appear to apply to any therapeutic protein, although mAbs clearly have some distinctive features. In this review, we first discuss chemical degradation reactions. This is followed by a section on physical instability issues. Then, more specific topics are addressed: instability induced by interactions with interfaces, predictive methods for physical stability and interplay between chemical and physical instability. The final parts are devoted to discussions how all the above impacts (co-)formulation strategies, in particular for high protein concentration solutions.'

Chagasin from Trypanosoma cruzi as a molecular scaffold to express epitopes of TSA-1 as soluble recombinant chimeras

Trypanosoma cruzi is the causative agent of Chagas disease, a global public health problem. New therapeutic drugs and biologics are needed. The TSA-1 recombinant protein of T. cruzi is one such promising antigen for developing a therapeutic vaccine. However, it is overexpressed in E. coli as inclusion bodies, requiring an additional refolding step. As an alternative, in this study, we propose the endogenous cysteine protease inhibitor chagasin as a molecular scaffold to generate chimeric proteins. These proteins will contain combinations of two of the five conserved epitopes (E1 to E5) of TSA-1 in the L4 and L6 chagasin loops. Twenty chimeras (Q1-Q20) were designed, and their solubility was predicted using bioinformatics tools. Nine chimeras with different degrees of solubility were selected and expressed in E. coli BL21 (DE3). Western blot assays with anti-6x-His and anti-chagasin antibodies confirmed the expression of soluble recombinant chimeras. Both theoretically and experimentally, the Q12 (E5-E3) chimera was the most soluble, and the Q20 (E4-E5) the most insoluble protein. Q4 (E5-E1) and Q8 (E5-E2) chimeras were classified as proteins with medium solubility that exhibited the highest yield in the soluble fraction. Notably, Q4 has a yield of 239 mg/L, well above the yield of recombinant chagasin (16.5 mg/L) expressed in a soluble form. The expression of the Q4 chimera was scaled up to a 7 L fermenter obtaining a yield of 490 mg/L. These data show that chagasin can serve as a molecular scaffold for the expression of TSA-1 epitopes in the form of soluble chimeras.

Automated optimization of the solubility of a hyper-stable α-amylase

Most successes in computational protein engineering to date have focused on enhancing one biophysical trait, while multi-trait optimization remains a challenge. Different biophysical properties are often conflicting, as mutations that improve one tend to worsen the others. In this study, we explored the potential of an automated computational design strategy, called CamSol Combination, to optimize solubility and stability of enzymes without affecting their activity. Specifically, we focus on Bacillus licheniformis α-amylase (BLA), a hyper-stable enzyme that finds diverse application in industry and biotechnology. We validate the computational predictions by producing 10 BLA variants, including the wild-type (WT) and three designed models harbouring between 6 and 8 mutations each. Our results show that all three models have substantially improved relative solubility over the WT, unaffected catalytic rate and retained hyper-stability, supporting the algorithm's capacity to optimize enzymes. High stability and solubility embody enzymes with superior resilience to chemical and physical stresses, enhance manufacturability and allow for high-concentration formulations characterized by extended shelf lives. This ability to readily optimize solubility and stability of enzymes will enable the rapid and reliable generation of highly robust and versatile reagents, poised to contribute to advancements in diverse scientific and industrial domains.

Investigating the effect of SUMO fusion on solubility and stability of amylase-catalytic domain from Pyrococcus abyssi

Alpha amylase belonging to starch hydrolyzing enzymes has significant contributions to different industrial processes. The enzyme production through recombinant DNA technology faces certain challenges related to their expression, solubility and purification, which can be overcome through fusion tags. This study explored the influence of SUMO, a protein tag reported to enhance the solubility and stability of target proteins when fused to the N-terminal of the catalytic domain of amylase from Pyrococcus abyssi (PaAD). The insoluble expression of PaAD in E. coli was overcome when the enzyme was expressed in a fusion state (S-PaAD) and culture was cultivated at 18 °C. Moreover, the activity of S-PaAD increased by 1.5-fold as compared to that of PaAD. The ligand binding and enzyme activity assays against different substrates demonstrated that it was more active against 1 % glycogen and amylopectin. The analysis of the hydrolysates through HPLC demonstrated that the enzyme activity is mainly amylolytic, producing longer oligosaccharides as the major end product. The secondary structure analyses by temperature ramping in CD spectroscopy and MD simulation demonstrated the enzymes in the free, as well as fusion state, were stable at 90 °C. The soluble production, thermostability and broad substrate specificity make this enzyme a promising choice for various foods, feed, textiles, detergents, pharmaceuticals, and many industrial applications.

Modulation of α-synuclein in vitro aggregation kinetics by its alternative splice isoforms

The misfolding and aggregation of α-synuclein is linked to a family of neurodegenerative disorders known as synucleinopathies, the most prominent of which is Parkinson's disease (PD). Understanding the aggregation process of α-synuclein from a mechanistic point of view is thus of key importance. SNCA, the gene encoding α-synuclein, comprises six exons and produces various isoforms through alternative splicing. The most abundant isoform is expressed as a 140-amino acid protein (αSyn-140), while three other isoforms, αSyn-126, αSyn-112, and αSyn-98, are generated by skipping exon 3, exon 5, or both exons, respectively. In this study, we performed a detailed biophysical characterization of the aggregation of these four isoforms. We found that αSyn-112 and αSyn-98 exhibit accelerated aggregation kinetics compared to αSyn-140 and form distinct aggregate morphologies, as observed by transmission electron microscopy. Moreover, we observed that the presence of relatively small amounts of αSyn-112 accelerates the aggregation of αSyn-140, significantly reducing the aggregation half-time. These results indicate a potential role of alternative splicing in the pathological aggregation of α-synuclein and provide insights into how this process could be associated with the development of synucleinopathies.

Targeting Efficient Features of Urate Oxidase to Increase Its Solubility

With the demand for mass production of protein drugs, solubility has become a serious issue. Extrinsic and intrinsic factors both affect this property. A homotetrameric cofactor-free urate oxidase (UOX) is not sufficiently soluble. To engineer UOX for optimum solubility, it is important to identify the most effective factor that influences solubility. The most effective feature to target for protein engineering was determined by measuring various solubility-related factors of UOX. A large library of homologous sequences was obtained from the databases. The data was reduced to six enzymes from different organisms. On the basis of various sequence- and structure-derived elements, the most and the least soluble enzymes were defined. To determine the best protein engineering target for modification, features of the most and least soluble enzymes were compared. Metabacillus fastidiosus UOX was the most soluble enzyme, while Agrobacterium globiformis UOX was the least soluble. According to the comparison-constant method, positive surface patches caused by arginine residue distribution are appropriate targets for modification. Two Arg to Ala mutations were introduced to the least soluble enzyme to test this hypothesis. These mutations significantly enhanced the mutant's solubility. While different algorithms produced conflicting results, it was difficult to determine which proteins were most and least soluble. Solubility prediction requires multiple algorithms based on these controversies. Protein surfaces should be investigated regionally rather than globally, and both sequence and structural data should be considered. Several other biotechnological products could be engineered using the data reduction and comparison-constant methods used in this study.

Production of Antibodies to Peptide Targets Using Hybridoma Technology

Hybridoma technology is a well-established and indispensable tool for generating high-quality monoclonal antibodies and has become one of the most common methods for monoclonal antibody production. In this process, antibody-producing B cells are isolated from mice following immunization of mice with a specific immunogen and fused with an immortal myeloma cell line to form antibody-producing hybridoma cell lines. Hybridoma-derived monoclonal antibodies not only serve as powerful research and diagnostic reagents but have also emerged as the most rapidly expanding class of therapeutic biologicals. In spite of the development of new high-throughput monoclonal antibody generation technologies, hybridoma technology still is applied for antibody production due to its ability to preserve innate functions of immune cells and to preserve natural cognate antibody paring information. In this chapter, an overview of hybridoma technology and the laboratory procedures used for hybridoma production and antibody screening of peptide-specific antibodies are presented.

Sequence-based prediction of the intrinsic solubility of peptides containing non-natural amino acids

Non-natural amino acids are increasingly used as building blocks in the development of peptide-based drugs as they expand the available chemical space to tailor function, half-life and other key properties. However, while the chemical space of modified amino acids (mAAs) such as residues containing post-translational modifications (PTMs) is potentially vast, experimental methods for measuring the developability properties of mAA-containing peptides are expensive and time consuming. To facilitate developability programs through computational methods, we present CamSol-PTM, a method that enables the fast and reliable sequence-based prediction of the intrinsic solubility of mAA-containing peptides in aqueous solution at room temperature. From a computational screening of 50,000 mAA-containing variants of three peptides, we selected five different small-size mAAs for a total number of 37 peptide variants for experimental validation. We demonstrate the accuracy of the predictions by comparing the calculated and experimental solubility values. Our results indicate that the computational screening of mAA-containing peptides can extend by over four orders of magnitude the ability to explore the solubility chemical space of peptides and confirm that our method can accurately assess the solubility of peptides containing mAAs. This method is available as a web server at https://www-cohsoftware.ch.cam.ac.uk/index.php/camsolptm .

Amyloid Aggregation Is Potently Slowed Down by Osmolytes Due to Compaction of Partially Folded State

Amyloid aggregation is a key process in amyloidoses and neurodegenerative diseases. Hydrophobicity is one of the major driving forces for this type of aggregation, as an increase in hydrophobicity generally correlates with aggregation susceptibility and rate. However, most experimental systems in vitro and prediction tools in silico neglect the contribution of protective osmolytes present in the cellular environment. Here, we assessed the role of hydrophobic mutations in amyloid aggregation in the presence of osmolytes. To achieve this goal, we used the model protein human muscle acylphosphatase (mAcP) and mutations to leucine that increased its hydrophobicity without affecting its thermodynamic stability. Osmolytes significantly slowed down the aggregation kinetics of the hydrophobic mutants, with an effect larger than that observed on the wild-type protein. The effect increased as the mutation site was closer to the middle of the protein sequence. We propose that the preferential exclusion of osmolytes from mutation-introduced hydrophobic side-chains quenches the aggregation potential of the ensemble of partially unfolded states of the protein by inducing its compaction and inhibiting its self-assembly with other proteins. Our results suggest that including the effect of the cellular environment in experimental setups and predictive softwares, for both mechanistic studies and drug design, is essential in order to obtain a more complete combination of the driving forces of amyloid aggregation.

Pangenomic and immunoinformatics based analysis of Nipah virus revealed CD4+ and CD8+ T-Cell epitopes as potential vaccine candidates

Introduction: Nipah (NiV) is the zoonotic deadly bat-borne virus that causes neurological and respiratory infections which ultimately lead to death. There are 706 infected cases reported up till now especially in Asia, out of which 409 patients died. There is no vaccine and effective treatment available for NiV infections and we have to timely design such strategies as world could not bear another pandemic situation. Methods: In this study, we screened viral proteins of NiV strains based on pangenomics analysis, antigenicity, molecular weight, and sub-cellular localization. The immunoproteomics based approach was used to predict T-cell epitopes of MHC class-I and II as potential vaccine candidates. These epitopes are capable to activate CD4+, CD8+, and T-cell dependent B-lymphocytes. Results: The two surface proteins including fusion glycoprotein (F) and attachment glycoprotein (G) are antigenic with molecular weights of 60 kDa and 67 kDa respectively. Three epitopes of F protein (VNYNSEGIA, PNFILVRNT, and IKMIPNVSN) were ranked and selected based on the binding affinity with MHC class-I, and 3 epitopes (VILNKRYYS, ILVRNTLIS, and VKLQETAEK) with MHC-II molecules. Similarly, for G protein, 3 epitopes each for MHC-I (GKYDKVMPY, ILKPKLISY, and KNKIWCISL) and MHC-II (LRNIEKGKY, FLIDRINWI, and FLLKNKIWC) with substantial binding energies were predicted. Based on the physicochemical properties, all these epitopes are non-toxic, hydrophilic, and stable. Conclusion: Our vaccinomics and system-level investigation could help to trigger the host immune system to prevent NiV infection.

Understanding and controlling the molecular mechanisms of protein aggregation in mAb therapeutics

In antibody development and manufacturing, protein aggregation is a common challenge that can lead to serious efficacy and safety issues. To mitigate this problem, it is important to investigate its molecular origins. This review discusses (1) our current molecular understanding and theoretical models of antibody aggregation, (2) how various stress conditions related to antibody upstream and downstream bioprocesses can trigger aggregation, and (3) current mitigation strategies employed towards inhibiting aggregation. We discuss the relevance of the aggregation phenomenon in the context of novel antibody modalities and highlight how in silico approaches can be exploited to mitigate it.

Surface patches induce nonspecific binding and phase separation of antibodies

Nonspecific interactions are a key challenge in the successful development of therapeutic antibodies. The tendency for nonspecific binding of antibodies is often difficult to reduce by rational design, and instead, it is necessary to rely on comprehensive screening campaigns. To address this issue, we performed a systematic analysis of the impact of surface patch properties on antibody nonspecificity using a designer antibody library as a model system and single-stranded DNA as a nonspecificity ligand. Using an in-solution microfluidic approach, we find that the antibodies tested bind to single-stranded DNA with affinities as high as KD = 1 µM. We show that DNA binding is driven primarily by a hydrophobic patch in the complementarity-determining regions. By quantifying the surface patches across the library, the nonspecific binding affinity is shown to correlate with a trade-off between the hydrophobic and total charged patch areas. Moreover, we show that a change in formulation conditions at low ionic strengths leads to DNA-induced antibody phase separation as a manifestation of nonspecific binding at low micromolar antibody concentrations. We highlight that phase separation is driven by a cooperative electrostatic network assembly mechanism of antibodies with DNA, which correlates with a balance between positive and negative charged patches. Importantly, our study demonstrates that both nonspecific binding and phase separation are controlled by the size of the surface patches. Taken together, these findings highlight the importance of surface patches and their role in conferring antibody nonspecificity and its macroscopic manifestation in phase separation.

Automated optimisation of solubility and conformational stability of antibodies and proteins

Biologics, such as antibodies and enzymes, are crucial in research, biotechnology, diagnostics, and therapeutics. Often, biologics with suitable functionality are discovered, but their development is impeded by developability issues. Stability and solubility are key biophysical traits underpinning developability potential, as they determine aggregation, correlate with production yield and poly-specificity, and are essential to access parenteral and oral delivery. While advances for the optimisation of individual traits have been made, the co-optimization of multiple traits remains highly problematic and time-consuming, as mutations that improve one property often negatively impact others. In this work, we introduce a fully automated computational strategy for the simultaneous optimisation of conformational stability and solubility, which we experimentally validate on six antibodies, including two approved therapeutics. Our results on 42 designs demonstrate that the computational procedure is highly effective at improving developability potential, while not affecting antigen-binding. We make the method available as a webserver at www-cohsoftware.ch.cam.ac.uk.

Sequence-based prediction of pH-dependent protein solubility using CamSol

Solubility is a property of central importance for the use of proteins in research in molecular and cell biology and in applications in biotechnology and medicine. Since experimental methods for measuring protein solubility are material intensive and time consuming, computational methods have recently emerged to enable the rapid and inexpensive screening of solubility for large libraries of proteins, as it is routinely required in development pipelines. Here, we describe the development of one such method to include in the predictions the effect of the pH on solubility. We illustrate the resulting pH-dependent predictions on a variety of antibodies and other proteins to demonstrate that these predictions achieve an accuracy comparable with that of experimental methods. We make this method publicly available at https://www-cohsoftware.ch.cam.ac.uk/index.php/camsolph, as the version 3.0 of CamSol.

Aggregation Limiting Cell-Penetrating Peptides Derived from Protein Signal Sequences

Alzheimer's disease (AD) is the most common neurodegenerative disease (ND) and the leading cause of dementia. It is characterized by non-linear, genetic-driven pathophysiological dynamics with high heterogeneity in the biological alterations and the causes of the disease. One of the hallmarks of the AD is the progression of plaques of aggregated amyloid-β (Aβ) or neurofibrillary tangles of Tau. Currently there is no efficient treatment for the AD. Nevertheless, several breakthroughs in revealing the mechanisms behind progression of the AD have led to the discovery of possible therapeutic targets. Some of these include the reduction in inflammation in the brain, and, although highly debated, limiting of the aggregation of the Aβ. In this work we show that similarly to the Neural cell adhesion molecule 1 (NCAM1) signal sequence, other Aβ interacting protein sequences, especially derived from Transthyretin, can be used successfully to reduce or target the amyloid aggregation/aggregates in vitro. The modified signal peptides with cell-penetrating properties reduce the Aβ aggregation and are predicted to have anti-inflammatory properties. Furthermore, we show that by expressing the Aβ-EGFP fusion protein, we can efficiently assess the potential for reduction in aggregation, and the CPP properties of peptides in mammalian cells.

Binding affinity prediction for antibody-protein antigen complexes: A machine learning analysis based on interface and surface areas

Specific antibodies can bind to protein antigens with high affinity and specificity, and this property makes them one of the best protein-based therapeutics. Accurate prediction of antibody‒protein antigen binding affinity is crucial for designing effective antibodies. The current predictive methods for protein‒protein binding affinity usually fail to predict the binding affinity of an antibody‒protein antigen complex with a comparable level of accuracy. Here, new models specific for antibody‒antigen binding affinity prediction are developed according to the different types of interface and surface areas present in antibody‒antigen complex. The contacts-based descriptors are also employed to construct or train different models specific for antibody‒protein antigen binding affinity prediction. The results of this study show that (i) the area-based descriptors are slightly better than the contacts-based descriptors in terms of the predictive power; (ii) the new models specific for antibody‒protein antigen binding affinity prediction are superior to the previously-used general models for predicting the protein‒protein binding affinities; (iii) the performances of the best area-based and contacts-based models developed in this work are better than the performances of a recently-developed graph-based model (i.e., CSM-AB) specific for antibody‒protein antigen binding affinity prediction. The new models developed in this work would not only help understand the mechanisms underlying antibody‒protein antigen interactions, but would also be of some applicable utility in the design and virtual screening of antibody-based therapeutics.

Non-specificity as the sticky problem in therapeutic antibody development

Antibodies are highly potent therapeutic scaffolds with more than a hundred different products approved on the market. Successful development of antibody-based drugs requires a trade-off between high target specificity and target binding affinity. In order to better understand this problem, we here review non-specific interactions and explore their fundamental physicochemical origins. We discuss the role of surface patches - clusters of surface-exposed amino acid residues with similar physicochemical properties - as inducers of non-specific interactions. These patches collectively drive interactions including dipole-dipole, π-stacking and hydrophobic interactions to complementary moieties. We elucidate links between these supramolecular assembly processes and macroscopic development issues, such as decreased physical stability and poor in vivo half-life. Finally, we highlight challenges and opportunities for optimizing protein binding specificity and minimizing non-specificity for future generations of therapeutics.

Hallucinating structure-conditioned antibody libraries for target-specific binders

Antibodies are widely developed and used as therapeutics to treat cancer, infectious disease, and inflammation. During development, initial leads routinely undergo additional engineering to increase their target affinity. Experimental methods for affinity maturation are expensive, laborious, and time-consuming and rarely allow the efficient exploration of the relevant design space. Deep learning (DL) models are transforming the field of protein engineering and design. While several DL-based protein design methods have shown promise, the antibody design problem is distinct, and specialized models for antibody design are desirable. Inspired by hallucination frameworks that leverage accurate structure prediction DL models, we propose the F_vHallucinator for designing antibody sequences, especially the CDR loops, conditioned on an antibody structure. Such a strategy generates targeted CDR libraries that retain the conformation of the binder and thereby the mode of binding to the epitope on the antigen. On a benchmark set of 60 antibodies, F_vHallucinator generates sequences resembling natural CDRs and recapitulates perplexity of canonical CDR clusters. Furthermore, the F_vHallucinator designs amino acid substitutions at the V_H-V_L interface that are enriched in human antibody repertoires and therapeutic antibodies. We propose a pipeline that screens F_vHallucinator designs to obtain a library enriched in binders for an antigen of interest. We apply this pipeline to the CDR H3 of the Trastuzumab-HER2 complex to generate in silico designs predicted to improve upon the binding affinity and interfacial properties of the original antibody. Thus, the F_vHallucinator pipeline enables generation of inexpensive, diverse, and targeted antibody libraries enriched in binders for antibody affinity maturation.

Computational Design and Experimental Validation of ACE2-Derived Peptides as SARS-CoV-2 Receptor Binding Domain Inhibitors

The COVID-19 pandemic has caused significant social and economic disruption across the globe. Cellular entry of SARS-CoV-2 into the human body is mediated via binding of the Receptor Binding Domain (RBD) on the viral Spike protein (SARS-CoV-2 RBD) to Angiotensin-Converting Enzyme 2 (ACE2) expressed on host cells. Molecules that can disrupt ACE2:RBD interactions are attractive therapeutic candidates to prevent virus entry into human cells. A computational strategy that combines our Peptide Binding Design (PepBD) algorithm with atomistic molecular dynamics simulations was used to design new inhibitory peptide candidates via sequence iteration starting with a 23-mer peptide, referred to as SBP1. SBP1 is derived from a region of the ACE2 Peptidase Domain α1 helix that binds to the SARS-CoV-2 RBD of the initial Wuhan-Hu-1 strain. Three peptides demonstrated a solution-phase RBD-binding dissociation constant in the micromolar range during tryptophan fluorescence quenching experiments, one peptide did not bind, and one was insoluble at micromolar concentrations. However, in competitive ELISA assays, none of these peptides could outcompete ACE2 binding to SARS-CoV-2-RBD up to concentrations of 50 μM, similar to the parent SBP1 peptide which also failed to outcompete ACE2:RBD binding. Molecular dynamics simulations suggest that P4 would have a good binding affinity for the RBD domain of Beta-B.1.351, Gamma-P.1, Kappa-B.1.617.1, Delta-B.1.617.2, and Omicron-B.1.1.529 variants, but not the Alpha variant. Consistent with this, P4 bound Kappa-B.1.617.1 and Delta-B.1.617.2 RBD with micromolar affinity in tryptophan fluorescence quenching experiments. Collectively, these data show that while relatively short unstructured peptides can bind to SARS-CoV-2 RBD with moderate affinity, they are incapable of outcompeting the strong interactions between RBD and ACE2.

Protein Design: From the Aspect of Water Solubility and Stability

Water solubility and structural stability are key merits for proteins defined by the primary sequence and 3D-conformation. Their manipulation represents important aspects of the protein design field that relies on the accurate placement of amino acids and molecular interactions, guided by underlying physiochemical principles. Emulated designer proteins with well-defined properties both fuel the knowledge-base for more precise computational design models and are used in various biomedical and nanotechnological applications. The continuous developments in protein science, increasing computing power, new algorithms, and characterization techniques provide sophisticated toolkits for solubility design beyond guess work. In this review, we summarize recent advances in the protein design field with respect to water solubility and structural stability. After introducing fundamental design rules, we discuss the transmembrane protein solubilization and de novo transmembrane protein design. Traditional strategies to enhance protein solubility and structural stability are introduced. The designs of stable protein complexes and high-order assemblies are covered. Computational methodologies behind these endeavors, including structure prediction programs, machine learning algorithms, and specialty software dedicated to the evaluation of protein solubility and aggregation, are discussed. The findings and opportunities for Cryo-EM are presented. This review provides an overview of significant progress and prospects in accurate protein design for solubility and stability.

A rational design of a multi-epitope vaccine against SARS-CoV-2 which accounts for the glycan shield of the spike glycoprotein

The ongoing global health crisis caused by Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2), the virus which leads to Coronavirus Disease 2019 (COVID-19) has impacted not only the health of people everywhere, but the economy in nations across the world. While vaccine candidates and therapeutics are currently undergoing clinical trials, there is a lack of proven effective treatments or cures for COVID-19. In this study, we have presented a synergistic computational platform, including molecular dynamics simulations and immunoinformatics techniques, to rationally design a multi-epitope vaccine candidate for COVID-19. This platform combines epitopes across Linear B Lymphocytes (LBL), Cytotoxic T Lymphocytes (CTL) and Helper T Lymphocytes (HTL) derived from both mutant and wild-type spike glycoproteins from SARS-CoV-2 with diverse protein conformations. In addition, this vaccine construct also takes the considerable glycan shield of the spike glycoprotein into account, which protects it from immune response. We have identified a vaccine candidate (a 35.9 kDa protein), named COVCCF, which is composed of 5 LBL, 6 HTL, and 6 CTL epitopes from the spike glycoprotein of SARS-CoV-2. Using multi-dose immune simulations, COVCCF induces elevated levels of immunoglobulin activity (IgM, IgG1, IgG2), and induces strong responses from B lymphocytes, CD4 T-helper lymphocytes, and CD8 T-cytotoxic lymphocytes. COVCCF induces cytokines important to innate immunity, including IFN-γ, IL4, and IL10. Additionally, COVCCF has ideal pharmacokinetic properties and low immune-related toxicities. In summary, this study provides a powerful, computational vaccine design platform for rapid development of vaccine candidates (including COVCCF) for effective prevention of COVID-19.Communicated by Ramaswamy H. Sarma.

Comparison of hydrophobicity scales for predicting biophysical properties of antibodies

While antibody-based therapeutics have grown to be one of the major classes of novel medicines, some antibody development candidates face significant challenges regarding expression levels, solubility, as well as stability and aggregation, under physiological and storage conditions. A major determinant of those properties is surface hydrophobicity, which promotes unspecific interactions and has repeatedly proven problematic in the development of novel antibody-based drugs. Multiple computational methods have been devised for in-silico prediction of antibody hydrophobicity, often using hydrophobicity scales to assign values to each amino acid. Those approaches are usually validated by their ability to rank potential therapeutic antibodies in terms of their experimental hydrophobicity. However, there is significant diversity both in the hydrophobicity scales and in the experimental methods, and consequently in the performance of in-silico methods to predict experimental results. In this work, we investigate hydrophobicity of monoclonal antibodies using hydrophobicity scales. We implement several scoring schemes based on the solvent-accessibility and the assigned hydrophobicity values, and compare the different scores and scales based on their ability to predict retention times from hydrophobic interaction chromatography. We provide an overview of the strengths and weaknesses of several commonly employed hydrophobicity scales, thereby improving the understanding of hydrophobicity in antibody development. Furthermore, we test several datasets, both publicly available and proprietary, and find that the diversity of the dataset affects the performance of hydrophobicity scores. We expect that this work will provide valuable guidelines for the optimization of biophysical properties in future drug discovery campaigns.

p53 TAD2 Domain (38-61) Forms Amyloid-like Aggregates in Isolation

A strong association between protein aggregation and human diseases (such as Alzheimer's, Parkinson's, and Huntington's disease) is well demonstrated. Misfolding and aggregation of p53, a central transcriptional mediator, has been revealed by various experimental evidence in different types of cancers. Aggregation studies focusing on different p53 domains, mostly, the central core domain and its mutants under the influence of various environmental conditions, and the p53 transactivation domain (TAD) (1-63) have been reported. However, the specific subdomains responsible for p53 aggregation are not known. p53 TADs interact with diverse cellular factors to modulate the function of p53 and elicit appropriate cellular responses under different stress conditions. In this study, the aggregation of the p53 TAD2 domain (38-61) has been studied in isolation. The aggregates were generated in vitro under acidic pH conditions after in silico scoring for amyloidogenic tendency and characterized using dye-based assays (ThT and bis-ANS fluorescence), CD spectroscopy, and microscopy (scanning electron microscoy, transmission electron microscopy, and atomic force microscopy). It was observed that p53 TAD2 forms characteristic β-sheet-rich amyloid-like fibrils. Via a reductionist approach, this study highlights the nature of p53 TAD2 domain (38-61) aggregation.

The effect of mutation on an aggregation-prone protein: An in vivo, in vitro, and in silico analysis

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.

Drosophila nicotinic acetylcholine receptor subunits and their native interactions with insecticidal peptide toxins

Drosophila nicotinic acetylcholine receptors (nAChRs) are ligand-gated ion channels that represent a target for insecticides. Peptide neurotoxins are known to block nAChRs by binding to their target subunits, however, a better understanding of this mechanism is needed for effective insecticide design. To facilitate the analysis of nAChRs we used a CRISPR/Cas9 strategy to generate null alleles for all ten nAChR subunit genes in a common genetic background. We studied interactions of nAChR subunits with peptide neurotoxins by larval injections and styrene maleic acid lipid particles (SMALPs) pull-down assays. For the null alleles, we determined the effects of α-Bungarotoxin (α-Btx) and ω-Hexatoxin-Hv1a (Hv1a) administration, identifying potential receptor subunits implicated in the binding of these toxins. We employed pull-down assays to confirm α-Btx interactions with the Drosophila α5 (Dα5), Dα6, Dα7 subunits. Finally, we report the localisation of fluorescent tagged endogenous Dα6 during Drosophila CNS development. Taken together, this study elucidates native Drosophila nAChR subunit interactions with insecticidal peptide toxins and provides a resource for the in vivo analysis of insect nAChRs.

Two physics‐based models for pH ‐dependent calculations of protein solubility

When engineering a protein for its biological function, many physicochemical properties are also optimized throughout the engineering process, and the protein's solubility is among the most important properties to consider. Here, we report two novel computational methods to calculate the pH-dependent protein solubility, and to rank the solubility of mutants. The first is an empirical method developed for fast ranking of the solubility of a large number of mutants of a protein. It takes into account electrostatic solvation energy term calculated using Generalized Born approximation, hydrophobic patches, protein charge, and charge asymmetry, as well as the changes of protein stability upon mutation. This method has been tested on over 100 mutations for 17 globular proteins, as well as on 44 variants of five different antibodies. The prediction rate is over 80%. The antibody tests showed a Pearson correlation coefficient, R, with experimental data from .83 to .91. The second method is based on a novel, completely force-field-based approach using CHARMm program modules to calculate the binding energy of the protein to a part of the crystal lattice, generated from X-ray structure. The method predicted with very high accuracy the solubility of Ribonuclease SA and its 3K and 5K mutants as a function of pH without any parameter adjustments of the existing BIOVIA Discovery Studio binding affinity model. Our methods can be used for rapid screening of large numbers of design candidates based on solubility, and to guide the design of solution conditions for antibody formulation.

Current advances in biopharmaceutical informatics: guidelines, impact and challenges in the computational developability assessment of antibody therapeutics

Therapeutic monoclonal antibodies and their derivatives are key components of clinical pipelines in the global biopharmaceutical industry. The availability of large datasets of antibody sequences, structures, and biophysical properties is increasingly enabling the development of predictive models and computational tools for the "developability assessment" of antibody drug candidates. Here, we provide an overview of the antibody informatics tools applicable to the prediction of developability issues such as stability, aggregation, immunogenicity, and chemical degradation. We further evaluate the opportunities and challenges of using biopharmaceutical informatics for drug discovery and optimization. Finally, we discuss the potential of developability guidelines based on in silico metrics that can be used for the assessment of antibody stability and manufacturability.

In Silico Identification and Experimental Validation of Peptide-Based Inhibitors Targeting Clostridium difficile Toxin A

Clostridium difficile infection is mediated by two major exotoxins: toxins A (TcdA) and B (TcdB). Inhibiting the biocatalytic activities of these toxins with targeted peptide-based drugs can reduce the risk of C. difficile infection. In this work, we used a computational strategy that integrates a peptide binding design (PepBD) algorithm and explicit-solvent atomistic molecular dynamics simulation to determine promising toxin A-targeting peptides that can recognize and bind to the catalytic site of the TcdA glucosyltransferase domain (GTD). Our simulation results revealed that two out of three in silico discovered peptides, viz. the neutralizing peptides A (NPA) and B (NPB), exhibit lower binding free energies when bound to the TcdA GTD than the phage-display discovered peptide, viz. the reference peptide (RP). These peptides may serve as potential inhibitors against C. difficile infection. The efficacy of the peptides RP, NPA, and NPB to neutralize the cytopathic effects of TcdA was tested in vitro in human jejunum cells. Both phage-display peptide RP and in silico peptide NPA were found to exhibit strong toxin-neutralizing properties, thereby preventing the TcdA toxicity. However, the in silico peptide NPB demonstrates a relatively low efficacy against TcdA.

Assessment of Therapeutic Antibody Developability by Combinations of In Vitro and In Silico Methods

Although antibodies have become the fastest-growing class of therapeutics on the market, it is still challenging to develop them for therapeutic applications, which often require these molecules to withstand stresses that are not present in vivo. We define developability as the likelihood of an antibody candidate with suitable functionality to be developed into a manufacturable, stable, safe, and effective drug that can be formulated to high concentrations while retaining a long shelf life. The implementation of reliable developability assessments from the early stages of antibody discovery enables flagging and deselection of potentially problematic candidates, while focussing available resources on the development of the most promising ones. Currently, however, thorough developability assessment requires multiple in vitro assays, which makes it labor intensive and time consuming to implement at early stages. Furthermore, accurate in vitro analysis at the early stage is compromised by the high number of potential candidates that are often prepared at low quantities and purity. Recent improvements in the performance of computational predictors of developability potential are beginning to change this scenario. Many computational methods only require the knowledge of the amino acid sequences and can be used to identify possible developability issues or to rank available candidates according to a range of biophysical properties. Here, we describe how the implementation of in silico tools into antibody discovery pipelines is increasingly offering time- and cost-effective alternatives to in vitro experimental screening, thus streamlining the drug development process. We discuss in particular the biophysical and biochemical properties that underpin developability potential and their trade-offs, review various in vitro assays to measure such properties or parameters that are predictive of developability, and give an overview of the growing number of in silico tools available to predict properties important for antibody development, including the CamSol method developed in our laboratory.

Progress and challenges for the machine learning-based design of fit-for-purpose monoclonal antibodies

Although the therapeutic efficacy and commercial success of monoclonal antibodies (mAbs) are tremendous, the design and discovery of new candidates remain a time and cost-intensive endeavor. In this regard, progress in the generation of data describing antigen binding and developability, computational methodology, and artificial intelligence may pave the way for a new era of in silico on-demand immunotherapeutics design and discovery. Here, we argue that the main necessary machine learning (ML) components for an in silico mAb sequence generator are: understanding of the rules of mAb-antigen binding, capacity to modularly combine mAb design parameters, and algorithms for unconstrained parameter-driven in silico mAb sequence synthesis. We review the current progress toward the realization of these necessary components and discuss the challenges that must be overcome to allow the on-demand ML-based discovery and design of fit-for-purpose mAb therapeutic candidates.

Development of QSAR models for in silico screening of antibody solubility

Although monoclonal antibodies (mAbs) have been shown to be extremely effective in treating a number of diseases, they often suffer from poor developability attributes, such as high viscosity and low solubility at elevated concentrations. Since experimental candidate screening is often materials and labor intensive, there is substantial interest in developing in silico tools for expediting mAb design. Here, we present a strategy using machine learning-based QSAR models for the a priori estimation of mAb solubility. The extrapolated protein solubilities of a set of 111 antibodies in a histidine buffer were determined using a high throughput PEG precipitation assay. 3D homology models of the antibodies were determined, and a large set of in house and commercially available molecular descriptors were then calculated. The resulting experimental and descriptor data were then used for the development of QSAR models of mAb solubilities. After feature selection and training with different machine learning algorithms, the models were evaluated with external test sets. The resulting regression models were able to estimate the solubility values of external test set data with R² of 0.81 and 0.85 for the two regression models developed. In addition, three class and binary classification models were developed and shown to be good estimators of mAb solubility behavior, with overall test set accuracies of 0.70 and 0.95, respectively. The analysis of the selected molecular descriptors in these models was also found to be informative and suggested that several charge-based descriptors and isotype may play important roles in mAb solubility. The combination of high throughput relative solubility experimental techniques in concert with efficient machine learning QSAR models offers an opportunity to rapidly screen potential mAb candidates and to design therapeutics with improved solubility characteristics.

A computational method for predicting the aggregation propensity of IgG1 and IgG4(P) mAbs in common storage buffers

The propensity for some monoclonal antibodies (mAbs) to aggregate at physiological and manufacturing pH values can prevent their use as therapeutic molecules or delay time to market. Consequently, developability assessments are essential to select optimum candidates, or inform on mitigation strategies to avoid potential late-stage failures. These studies are typically performed in a range of buffer solutions because factors such as pH can dramatically alter the aggregation propensity of the test mAbs (up to 100-fold in extreme cases). A computational method capable of robustly predicting the aggregation propensity at the pH values of common storage buffers would have substantial value. Here, we describe a mAb aggregation prediction tool (MAPT) that builds on our previously published isotype-dependent, charge-based model of aggregation. We show that the addition of a homology model-derived hydrophobicity descriptor to our electrostatic aggregation model enabled the generation of a robust mAb developability indicator. To contextualize our aggregation scoring system, we analyzed 97 clinical-stage therapeutic mAbs. To further validate our approach, we focused on six mAbs (infliximab, tocilizumab, rituximab, CNTO607, MEDI1912 and MEDI1912_STT) which have been reported to cover a large range of aggregation propensities. The different aggregation propensities of the case study molecules at neutral and slightly acidic pH were correctly predicted, verifying the utility of our computational method.

Design of millidevices to expedite apparent solubility measurements

A fast, automated and accurate millidevice for determination of the apparent solubility of proteins and impurities and different industrially relevant precipitating agents.

Explicit solvation thermodynamics in ionic solution: extending grid inhomogeneous solvation theory to solvation free energy of salt-water mixtures

Hydration thermodynamics play a fundamental role in fields ranging from the pharmaceutical industry to environmental research. Numerous methods exist to predict solvation thermodynamics of compounds ranging from small molecules to large biomolecules. Arguably the most precise methods are those based on molecular dynamics (MD) simulations in explicit solvent. One theory that has seen increased use is inhomogeneous solvation theory (IST). However, while many applications require accurate description of salt-water mixtures, no implementation of IST is currently able to estimate solvation properties involving more than one solvent species. Here, we present an extension to grid inhomogeneous solvation theory (GIST) that can take salt contributions into account. At the example of carbazole in 1 M NaCl solution, we compute the solvation energy as well as first and second order entropies. While the effect of the first order ion entropy is small, both the water-water and water-ion entropies contribute strongly. We show that the water-ion entropies are efficiently approximated using the Kirkwood superposition approximation. However, this approach cannot be applied to the water-water entropy. Furthermore, we test the quantitative validity of our method by computing salting-out coefficients and comparing them to experimental data. We find a good correlation to experimental salting-out constants, while the absolute values are overpredicted due to the approximate second order entropy. Since ions are frequently used in MD, either to neutralize the system or as a part of the investigated process, our method greatly extends the applicability of GIST. The use-cases range from biopharmaceuticals, where many assays require high salt concentrations, to environmental research, where solubility in sea water is important to model the fate of organic substances.

The Therapeutic Antibody Profiler for Computational Developability Assessment

The need to consider an antibody's "developability" (immunogenicity, solubility, specificity, stability, manufacturability, and storability) is now well understood in therapeutic antibody design. Predicting these properties rapidly and inexpensively is critical to industrial workflows, to avoid devoting resources to non-productive candidates. Here, we describe a high-throughput computational developability assessment tool, the Therapeutic Antibody Profiler (TAP), which assesses the physicochemical "druglikeness" of an antibody candidate. Input variable domain sequences are converted to three-dimensional structural models, and then five developability-linked molecular surface descriptors are calculated and compared to advanced-stage clinical therapeutics. Values at the extremes of/outside of the distributions seen in therapeutics imply an increased risk of developability issues. Therefore, TAP, starting only from sequence information, provides a route to rapidly identifying drug candidate antibodies that are likely to have poor developability. Our web application ( opig.stats.ox.ac.uk/webapps/tap ) profiles input antibody sequences against a continually updated reference set of clinical therapeutics.

An open-source automated PEG precipitation assay to measure the relative solubility of proteins with low material requirement

The solubility of proteins correlates with a variety of their properties, including function, production yield, pharmacokinetics, and formulation at high concentrations. High solubility is therefore a key requirement for the development of protein-based reagents for applications in life sciences, biotechnology, diagnostics, and therapeutics. Accurate solubility measurements, however, remain challenging and resource intensive, which limits their throughput and hence their applicability at the early stages of development pipelines, when long-lists of candidates are typically available in minute amounts. Here, we present an automated method based on the titration of a crowding agent (polyethylene glycol, PEG) to quantitatively assess relative solubility of proteins using about 200 µg of purified material. Our results demonstrate that this method is accurate and economical in material requirement and costs of reagents, which makes it suitable for high-throughput screening. This approach is freely-shared and based on a low cost, open-source liquid-handling robot. We anticipate that this method will facilitate the assessment of the developability of proteins and make it substantially more accessible.

Stress-dependent flexibility of a full-length human monoclonal antibody: Insights from molecular dynamics to support biopharmaceutical development

After several decades of advancements in drug discovery, product development of biopharmaceuticals remains a time- and resource-consuming endeavor. One of the main reasons is associated to the lack of fundamental understanding of conformational dynamics of such biologic entities, and how they respond to various stresses encountered during manufacturing. In this work, we have studied the conformational dynamics of human IgG₁κ b12 monoclonal antibody (mAb) using molecular dynamics simulations. The hundreds of nanoseconds long trajectories reveal that b12 mAb is highly flexible. Its variable domains show greater conformational fluctuations than the constant domains. Additionally, it collapses towards a more globular shape in response to thermal stress, leading to decrease in the total solvent exposed surface area and radius of gyration. This behavior is more pronounced for the deglycosylated b12 mAb, and it appears to correlate with increase in inter-domain contacts between specific regions of the antibody. Conformational fluctuations also cause temporary formation and disruption of hydrophobic and charged patches on the antibody surface, which is particularly important for the prediction of CMC properties during development phases of antibody-based biotherapeutics. The insights gained through these simulations may help the development of biologic drugs, especially with regards to manufacturing processes where antibodies may undergo significant thermal stress.

De novo design of peptides that coassemble into β sheet-based nanofibrils

Peptides’ hierarchical coassembly into nanostructures enables controllable fabrication of multicomponent biomaterials. In this work, we describe a computational and experimental approach to design pairs of charge-complementary peptides that selectively coassemble into β-sheet nanofibers when mixed together but remain unassembled when isolated separately. The key advance is a peptide coassembly design (PepCAD) algorithm that searches for pairs of coassembling peptides. Six peptide pairs are identified from a pool of ~10⁶ candidates via the PepCAD algorithm and then subjected to DMD/PRIME20 simulations to examine their co-/self-association kinetics. The five pairs that spontaneously aggregate in kinetic simulations selectively coassemble in biophysical experiments, with four forming β-sheet nanofibers and one forming a stable nonfibrillar aggregate. Solid-state NMR, which is applied to characterize the coassembling pairs, suggests that the in silico peptides exhibit a higher degree of structural order than the previously reported CATCH(+/−) peptides.

Machine Learning Attempts for Predicting Human Subcutaneous Bioavailability of Monoclonal Antibodies

PURPOSE

One knowledge gap related to subcutaneous (SC) delivery is unpredictable and variable bioavailability. This study was aimed to develop machine learning methods to predict whether mAb's bioavailability was ≥70% or below, without completely knowing the mechanism and causality between inputs and outputs.

METHODS

A database of mAb SC products was built. The model training and validation were accomplished based on this database and a set of the inputs (product properties) were mapped to the output (bioavailability) using different machine learning algorithms. Dimensionality reduction was undertaken using principal component analysis (PCA).

RESULTS

The bioavailability of the mAb products being investigated varied from 35% to 90%. The tree-based methods, including random forest (RF), Adaptive Boost (AdaBoost), and decision tree (DT) presented the best predictability and generalization power on bioavailability classification. The models based on Multi-layer perceptron (MLP), Gaussian Naïve Bayes (GaussianNB), and k nearest neighbor (kNN) algorithms also provided acceptable prediction accuracy.

CONCLUSION

Machine learning could be a potential tool to predict mAb's bioavailability. Since all input features were acquired using theoretical calculations and predictions rather than experiments, the models may be particularly applicable to some early-stage research activities such as mAb molecule triage, design/optimization, mutant screening, molecule selection, and formulation design.

Solubility-Weighted Index: fast and accurate prediction of protein solubility

Motivation

Recombinant protein production is a widely used technique in the biotechnology and biomedical industries, yet only a quarter of target proteins are soluble and can therefore be purified.

Results

We have discovered that global structural flexibility, which can be modeled by normalized B-factors, accurately predicts the solubility of 12 216 recombinant proteins expressed in Escherichia coli. We have optimized these B-factors, and derived a new set of values for solubility scoring that further improves prediction accuracy. We call this new predictor the ‘Solubility-Weighted Index’ (SWI). Importantly, SWI outperforms many existing protein solubility prediction tools. Furthermore, we have developed ‘SoDoPE’ (Soluble Domain for Protein Expression), a web interface that allows users to choose a protein region of interest for predicting and maximizing both protein expression and solubility.

Availability and implementation

The SoDoPE web server and source code are freely available at https://tisigner.com/sodope and https://github.com/Gardner-BinfLab/TISIGNER-ReactJS, respectively. The code and data for reproducing our analysis can be found at https://github.com/Gardner-BinfLab/SoDoPE_paper_2020.

Supplementary information

Supplementary data are available at Bioinformatics online.

Investigation of interactions between binding residues and solubility of grafted humanized anti-VEGF IgG antibodies expressed as full-length format in the cytoplasm of a novel engineered E. coli SHuffle strain

Monoclonal antibodies (mAbs) are one of the fastest-growing areas of biopharmaceutical industry and have been widely used for a broad spectrum of diseases. Meanwhile, the immunogenicity of non-human derived antibodies can generate side effects by inducing the human immune response to produce human anti-mouse-immunoglobulin antibody (HAMA). In this work, we aim to reduce the immunogenicity of muMAb A.4.6.1 by substitute human sequences for murine sequences. Humanized antibodies are constructed by grafting, specificity determining residues (SDR), complementary determining regions (CDR), and chimeric region of muMAb A.4.6.1, onto variable domain of Trastuzumab (Herceptin). The interactions between grafted antibodies and their target, Vascular endothelial growth factor (VEGF), were theoretically investigated by molecular dynamics simulation in order to evaluate the antibodies-antigen binding behavior. The obtained protein-protein interactions and calculated binding free energy suggested that the SDR-VEGF complex presented a significantly greater binding affinity, number of contact and total number of H-bonds compared to CDR and chimeric mAbs, significantly. Moreover, the Camsol program predicted that the solubility of SDR mAb exhibits the greatest solubility. This result was supported by performing a western blot analysis of the grafted mAbs with soluble and insoluble fractions in order to evaluate their solubility, in which SDR was found to have a much lower amount of insoluble proteins. Consequently, the enhanced binding affinity and solubility of the designed SDR was achieved by the single S106D mutation using computational methods. With the aim of low immunogenicity, high solubility, and high affinity, this SDR humanized antibody was expected to have greater efficacy than murine or chimeric antibodies for future use in humans.

Effect of Humanizing Mutations on the Stability of the Llama Single-Domain Variable Region

In vivo clinical applications of nanobodies (VHHs) require molecules that induce minimal immunoresponse and therefore possess sequences as similar as possible to the human VH domain. Although the relative sequence variability in llama nanobodies has been used to identify scaffolds with partially humanized signature, the transformation of the Camelidae hallmarks in the framework2 still represents a major problem. We assessed a set of mutants in silico and experimentally to elucidate what is the contribution of single residues to the VHH stability and how their combinations affect the mutant nanobody stability. We described at molecular level how the interaction among residues belonging to different structural elements enabled a model llama nanobody (C8WT, isolated from a naïve library) to be functional and maintain its stability, despite the analysis of its primary sequence would classify it as aggregation-prone. Five chimeras formed by grafting CDRs isolated from different nanobodies into C8WT scaffold were successfully expressed as soluble proteins and both tested clones preserved their antigen binding specificity. We identified a nanobody with human hallmarks that seems suitable for humanizing selected camelid VHHs by grafting heterologous CDRs in its scaffold and could serve for the preparation of a synthetic library of human-like single domains.

Discovery-stage identification of drug-like antibodies using emerging experimental and computational methods

There is intense and widespread interest in developing monoclonal antibodies as therapeutic agents to treat diverse human disorders. During early-stage antibody discovery, hundreds to thousands of lead candidates are identified, and those that lack optimal physical and chemical properties must be deselected as early as possible to avoid problems later in drug development. It is particularly challenging to characterize such properties for large numbers of candidates with the low antibody quantities, concentrations, and purities that are available at the discovery stage, and to predict concentrated antibody properties (e.g., solubility, viscosity) required for efficient formulation, delivery, and efficacy. Here we review key recent advances in developing and implementing high-throughput methods for identifying antibodies with desirable in vitro and in vivo properties, including favorable antibody stability, specificity, solubility, pharmacokinetics, and immunogenicity profiles, that together encompass overall drug developability. In particular, we highlight impressive recent progress in developing computational methods for improving rational antibody design and prediction of drug-like behaviors that hold great promise for reducing the amount of required experimentation. We also discuss outstanding challenges that will need to be addressed in the future to fully realize the great potential of using such analysis for minimizing development times and improving the success rate of antibody candidates in the clinic.

Machine Learning Feature Selection for Predicting High Concentration Therapeutic Antibody Aggregation

Protein aggregation can hinder the development, safety and efficacy of therapeutic antibody-based drugs. Developing a predictive model that evaluates aggregation behaviors during early stage development is therefore desirable. Machine learning is a widely used tool to train models that predict data with different attributes. However, most machine learning techniques require more data than is typically available in antibody development. In this work, we describe a rational feature selection framework to develop accurate models with a small number of features. We applied this framework to predict aggregation behaviors of 21 approved monospecific monoclonal antibodies at high concentration (150 mg/mL), yielding a correlation coefficient of 0.71 on validation tests with only two features using a linear model. The nearest neighbors and support vector regression models further improved the performance, which have correlation coefficients of 0.86 and 0.80, respectively. This framework can be extended to train other models that predict different physical properties.

From cell line development to the formulated drug product: The art of manufacturing therapeutic monoclonal antibodies

Therapeutic monoclonal antibodies and related products have steadily grown to become the dominant product class within the biopharmaceutical market. Production of antibodies requires special precautions to ensure safety and efficacy of the product. In particular, minimizing antibody product heterogeneity is crucial as drug substance variants may impair the activity, efficacy, safety, and pharmacokinetic properties of an antibody, consequently resulting in the failure of a product in pre-clinical and clinical development. This review will cover the manufacturing and formulation challenges and advances of therapeutic monoclonal antibodies, focusing on improved processes to minimize variants and ensure batch-to-batch consistency. Processes put in place by regulatory agencies, such as Quality-by-Design (QbD) and current Good Manufacturing Practices (cGMP), and how their implementation has aided drug development in pharmaceutical companies will be reviewed. Advances in formulation and considerations on the intended use of a therapeutic antibody, including the route of administration and patient compliance, will be discussed.

Cross-talk between individual phenol-soluble modulins in Staphylococcus aureus biofilm enables rapid and efficient amyloid formation

The infective ability of the opportunistic pathogen Staphylococcus aureus, recognized as the most frequent cause of biofilm-associated infections, is associated with biofilm-mediated resistance to host immune response. Phenol-soluble modulins (PSM) comprise the structural scaffold of S. aureus biofilms through self-assembly into functional amyloids, but the role of individual PSMs during biofilm formation remains poorly understood and the molecular pathways of PSM self-assembly are yet to be identified. Here we demonstrate high degree of cooperation between individual PSMs during functional amyloid formation. PSMα3 initiates the aggregation, forming unstable aggregates capable of seeding other PSMs resulting in stable amyloid structures. Using chemical kinetics we dissect the molecular mechanism of aggregation of individual PSMs showing that PSMα1, PSMα3 and PSMβ1 display secondary nucleation whereas PSMβ2 aggregates through primary nucleation and elongation. Our findings suggest that various PSMs have evolved to ensure fast and efficient biofilm formation through cooperation between individual peptides.

A rationally designed bicyclic peptide remodels Aβ42 aggregation in vitro and reduces its toxicity in a worm model of Alzheimer's disease

Bicyclic peptides have great therapeutic potential since they can bridge the gap between small molecules and antibodies by combining a low molecular weight of about 2 kDa with an antibody-like binding specificity. Here we apply a recently developed in silico rational design strategy to produce a bicyclic peptide to target the C-terminal region (residues 31–42) of the 42-residue form of the amyloid β peptide (Aβ42), a protein fragment whose aggregation into amyloid plaques is linked with Alzheimer’s disease. We show that this bicyclic peptide is able to remodel the aggregation process of Aβ42 in vitro and to reduce its associated toxicity in vivo in a C. elegans worm model expressing Aβ42. These results provide an initial example of a computational approach to design bicyclic peptides to target specific epitopes on disordered proteins.

Heterologous Expression of Immature Forms of Human Islet Amyloid Polypeptide in Yeast Triggers Intracellular Aggregation and Cytotoxicity

Diabetes is a major public health issue that has attained alarming levels worldwide. Pancreatic aggregates of human islet amyloid polypeptide (IAPP) represent a major histopathological hallmark of type 2 diabetes. IAPP is expressed in β-cells as pre-pro-IAPP (ppIAPP) that is first processed to pro-IAPP (pIAPP) and finally to its mature form (matIAPP), being released upon glucose stimulation together with insulin. Impairment and overload of the IAPP processing machinery seem to be associated with the accumulation of immature IAPP species and the formation of toxic intracellular oligomers, which have been associated with β-cell dyshomeostasis and apoptosis. Nevertheless, the pathological importance of these immature IAPP forms for the assembly and cytotoxicity of these oligomers is not completely understood. Here, we describe the generation and characterization of unprecedented Saccharomyces cerevisiae models recapitulating IAPP intracellular oligomerization. Expression of green fluorescent protein (GFP) fusions of human ppIAPP, pIAPP, and matIAPP proved to be toxic in yeast cells at different extents, with ppIAPP exerting the most deleterious effect on yeast growth and cell viability. Although expression of all IAPP constructs induced the formation of intracellular aggregates in yeast cells, our data point out the accumulation of insoluble oligomeric species enriched in immature ppIAPP as the trigger of the high toxicity mediated by this construct in cells expressing ppIAPP-GFP. In addition, MS/MS analysis indicated that oligomeric species found in the ppIAPP-GFP lysates contain the N-terminal sequence of the propeptide fused to GFP. These models represent powerful tools for future research focused on the relevance of immature forms in IAPP-induced toxicity. Furthermore, they are extremely useful in high-throughput screenings for genetic and chemical modulators of IAPP aggregation.