Cross-Interaction Chromatography-Based QSAR Model for Early-Stage Screening to Facilitate Enhanced Developability of Monoclonal Antibody Therapeutics

Current state of implementation of in silico tools in the biopharmaceutical industry-Proceedings of the 5th modeling workshop

The fifth modeling workshop (5MW) was held in June 2023 at Favrholm, Denmark and sponsored by Recovery of Biological Products Conference Series. The goal of the workshop was to assemble modeling practitioners to review and discuss the current state, progress since the last fourth mini modeling workshop (4MMW), gaps and opportunities for development, deployment and maintenance of models in bioprocess applications. Areas of focus were four categories: biophysics and molecular modeling, mechanistic modeling, computational fluid dynamics (CFD) and plant modeling. Highlights of the workshop included significant advancements in biophysical/molecular modeling to novel protein constructs, mechanistic models for filtration and initial forays into modeling of multiphase systems using CFD for a bioreactor and mapped strategically to cell line selection/facility fit. A significant impediment to more fully quantitative and calibrated models for biophysics is the lack of large, anonymized datasets. A potential solution would be the use of specific descriptors in a database that would allow for detailed analyzes without sharing proprietary information. Another gap identified was the lack of a consistent framework for use of models that are included or support a regulatory filing beyond the high-level guidance in ICH Q8-Q11. One perspective is that modeling can be viewed as a component or precursor of machine learning (ML) and artificial intelligence (AI). Another outcome was alignment on a key definition for "mechanistic modeling." Feedback from participants was that there was progression in all of the fields of modeling within scope of the conference. Some areas (e.g., biophysics and molecular modeling) have opportunities for significant research investment to realize full impact. However, the need for ongoing research and development for all model types does not preclude the application to support process development, manufacturing and use in regulatory filings. Analogous to ML and AI, given the current state of the four modeling types, a prospective investment in educating inter-disciplinary subject matter experts (e.g., data science, chromatography) is essential to advancing the modeling community.

A multiscale modeling method for therapeutic antibodies in ion exchange chromatography

The development of biopharmaceutical downstream processes relies on exhaustive experimental studies. The root cause is the poorly understood relationship between the protein structure of monoclonal antibodies (mAbs) and their macroscopic process behavior. Especially the development of preparative chromatography processes is challenged by the increasing structural complexity of novel antibody formats and accelerated development timelines. This study introduces a multiscale in silico model consisting of homology modeling, quantitative structure-property relationships (QSPR), and mechanistic chromatography modeling leading from the amino acid sequence of a mAb to the digital representation of its cation exchange chromatography (CEX) process. The model leverages the mAbs' structural characteristics and experimental data of a diverse set of 21 therapeutic antibodies to predict elution profiles of two mAbs that were removed from the training data set. QSPR modeling identified mAb-specific protein descriptors relevant for the prediction of the thermodynamic equilibrium and the stoichiometric coefficient of the adsorption reaction. The consideration of two discrete conformational states of IgG4 mAbs enabled prediction of split-peak elution profiles. Starting from the sequence, the presented multiscale model allows in silico development of chromatography processes before protein material is available for experimental studies.

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.

A Descriptor set for Quantitative Structure-Property Relationship Prediction in Biologics

In addition to attaining the desired binding to their targets, a crucial aspect in the development of biotherapeutics is 'developability', which includes several desirable properties such as high solubility, low viscosity and aggregation, physico-chemical stability and low immunogenicity. The lack of any of these properties can lead to significant obstacles in advancing them to clinic; thus in silico methods capable of raising warning flags in earlier stages of development are highly beneficial. We have developed a computational framework based on a large and diverse set of protein specific descriptors ideal for making liability predictions using a machine-learning approach. This set offers a high degree of feature diversity classifiable by sequence, structure and surface patches. We assess the sensitivity and applicability of these descriptors in four dedicated case studies that are believed to be representative of biophysical characterizations commonly employed during the development process. In addition to data sets obtained from public sources, we have validated the descriptors on novel experimental data sets in order to address antibody developability and to generate prospective predictions on Adnectins. The results demonstrate that the descriptors are well suited to assist in the improvement of properties of systems that exhibit poor solubility or aggregation.

Predicting Antibody Developability Profiles Through Early Stage Discovery Screening

Monoclonal antibodies play an increasingly important role for the development of new drugs across multiple therapy areas. The term 'developability' encompasses the feasibility of molecules to successfully progress from discovery to development via evaluation of their physicochemical properties. These properties include the tendency for self-interaction and aggregation, thermal stability, colloidal stability, and optimization of their properties through sequence engineering. Selection of the best antibody molecule based on biological function, efficacy, safety, and developability allows for a streamlined and successful CMC phase. An efficient and practical high-throughput developability workflow (100 s-1,000 s of molecules) implemented during early antibody generation and screening is crucial to select the best lead candidates. This involves careful assessment of critical developability parameters, combined with binding affinity and biological properties evaluation using small amounts of purified material (<1 mg), as well as an efficient data management and database system. Herein, a panel of 152 various human or humanized monoclonal antibodies was analyzed in biophysical property assays. Correlations between assays for different sets of properties were established. We demonstrated in two case studies that physicochemical properties and key assay endpoints correlate with key downstream process parameters. The workflow allows the elimination of antibodies with suboptimal properties and a rank ordering of molecules for further evaluation early in the candidate selection process. This enables any further engineering for problematic sequence attributes without affecting program timelines.

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.

QSAR Implementation for HIC Retention Time Prediction of mAbs Using Fab Structure: A Comparison between Structural Representations

Monoclonal antibodies (mAbs) constitute a rapidly growing biopharmaceutical sector. However, their growth is impeded by high failure rates originating from failed clinical trials and developability issues in process development. There is, therefore, a growing need for better in silico tools to aid in risk assessment of mAb candidates to promote early-stage screening of potentially problematic mAb candidates. In this study, a quantitative structure–activity relationship (QSAR) modelling workflow was designed for the prediction of hydrophobic interaction chromatography (HIC) retention times of mAbs. Three novel descriptor sets derived from primary sequence, homology modelling, and atomistic molecular dynamics (MD) simulations were developed and assessed to determine the necessary level of structural resolution needed to accurately capture the relationship between mAb structures and HIC retention times. The results showed that descriptors derived from 3D structures obtained after MD simulations were the most suitable for HIC retention time prediction with a R² = 0.63 in an external test set. It was found that when using homology modelling, the resulting 3D structures became biased towards the used structural template. Performing an MD simulation therefore proved to be a necessary post-processing step for the mAb structures in order to relax the structures and allow them to attain a more natural conformation. Based on the results, the proposed workflow in this paper could therefore potentially contribute to aid in risk assessment of mAb candidates in early development.

Charge and hydrophobicity are key features in sequence-trained machine learning models for predicting the biophysical properties of clinical-stage antibodies

Improved understanding of properties that mediate protein solubility and resistance to aggregation are important for developing biopharmaceuticals, and more generally in biotechnology and synthetic biology. Recent acquisition of large datasets for antibody biophysical properties enables the search for predictive models. In this report, machine learning methods are used to derive models for 12 biophysical properties. A physicochemical perspective is maintained in analysing the models, leading to the observation that models cluster largely according to charge (cross-interaction measurements) and hydrophobicity (self-interaction methods). These two properties also overlap in some cases, for example in a new interpretation of variation in hydrophobic interaction chromatography. Since the models are developed from differences of antibody variable loops, the next stage is to extend models to more diverse protein sets.

AVAILABILITY

The web application for the sequence-based algorithms are available on the protein-sol webserver, at https://protein-sol.manchester.ac.uk/abpred, with models and virtualisation software available at https://protein-sol.manchester.ac.uk/software.