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.

Toward in silico CMC: An industrial collaborative approach to model-based process development.

The Third Modeling Workshop focusing on bioprocess modeling was held in Kenilworth, NJ in May 2019. A summary of these Workshop proceedings is captured in this manuscript. Modeling is an active area of research within the biotechnology community, and there is a critical need to assess the current state and opportunities for continued investment to realize the full potential of models, including resource and time savings. Beyond individual presentations and topics of novel interest, a substantial portion of the Workshop was devoted toward group discussions of current states and future directions in modeling fields. All scales of modeling, from biophysical models at the molecular level and up through large scale facility and plant modeling, were considered in these discussions and are summarized in the manuscript. Model life cycle management from model development to implementation and sustainment are also considered for different stages of clinical development and commercial production. The manuscript provides a comprehensive overview of bioprocess modeling while suggesting an ideal future state with standardized approaches aligned across the industry.

Highland games: A benchmarking exercise in predicting biophysical and drug properties of monoclonal antibodies from amino acid sequences.

Biopharmaceutical product and process development do not yet take advantage of predictive computational modeling to nearly the degree seen in industries based on smaller molecules. To assess and advance progress in this area, spirited coopetition (mutually beneficial collaboration between competitors) was successfully used to motivate industrial scientists to develop, share, and compare data and methods which would normally have remained confidential. The first "Highland Games" competition was held in conjunction with the October 2018 Recovery of Biological Products Conference in Ashville, NC, with the goal of benchmarking and assessment of the ability to predict development-related properties of six antibodies from their amino acid sequences alone. Predictions included purification-influencing properties such as isoelectric point and protein A elution pH, and biophysical properties such as stability and viscosity at very high concentrations. Essential contributions were made by a large variety of individuals, including companies which consented to provide antibody amino acid sequences and test materials, volunteers who undertook the preparation and experimental characterization of these materials, and prediction teams who attempted to predict antibody properties from sequence alone. Best practices were identified and shared, and areas in which the community excels at making predictions were identified, as well as areas presenting opportunities for considerable improvement. Predictions of isoelectric point and protein A elution pH were especially good with all-prediction average errors of 0.2 and 1.6 pH unit, respectively, while predictions of some other properties were notably less good. This manuscript presents the events, methods, and results of the competition, and can serve as a tutorial and as a reference for in-house benchmarking by others. Organizations vary in their policies concerning disclosure of methods, but most managements were very cooperative with the Highland Games exercise, and considerable insight into common and best practices is available from the contributed methods. The accumulated data set will serve as a benchmarking tool for further development of in silico prediction tools.

Prediction uncertainty assessment of chromatography models using Bayesian inference.

Mechanistic modeling of chromatography has been around in academia for decades and has gained increased support in pharmaceutical companies in recent years. Despite the large number of published successful applications, process development in the pharmaceutical industry today still does not fully benefit from a systematic mechanistic model-based approach. The hesitation on the part of industry to systematically apply mechanistic models can often be attributed to the absence of a general approach for determining if a model is qualified to support decision making in process development. In this work a Bayesian framework for the calibration and quality assessment of mechanistic chromatography models is introduced. Bayesian Markov Chain Monte Carlo is used to assess parameter uncertainty by generating samples from the parameter posterior distribution. Once the parameter posterior distribution has been estimated, it can be used to propagate the parameter uncertainty to model predictions, allowing a prediction-based uncertainty assessment of the model. The benefit of this uncertainty assessment is demonstrated using the example of a mechanistic model describing the separation of an antibody from its impurities on a strong cation exchanger. The mechanistic model was calibrated at moderate column load density and used to make extrapolations at high load conditions. Using the Bayesian framework, it could be shown that despite significant parameter uncertainty, the model can extrapolate beyond observed process conditions with high accuracy and is qualified to support process development.

Modeling of complex antibody elution behavior under high protein load densities in ion exchange chromatography using an asymmetric activity coefficient.

A main requirement for the implementation of model-based process development in industry is the capability of the model to predict high protein load densities. The frequently used steric mass action isotherm assumes a thermodynamically ideal system and, hence constant activity coefficients. In this manuscript, an industrial antibody purification problem under high load conditions is considered where this assumption does not hold. The high protein load densities, as commonly applied in industrial downstream processing, may lead to complex elution peak shapes. Using Mollerup's generalized ion-exchange isotherm (GIEX), the observed elution peak shapes could be modeled. To this end, the GIEX isotherm introduced two additional parameters to approximate the asymmetric activity coefficient. The effects of these two parameters on the curvature of the adsorption isotherm and the resulting chromatogram are investigated. It could be shown that they can be determined by inverse peak fitting and conform with the mechanistic demands of model-based process development.

SAP couples Fyn to SLAM immune receptors.

SAP (SLAM-associated protein) is a small lymphocyte-specific signalling molecule that is defective or absent in patients with X-linked lymphoproliferative syndrome (XLP). Consistent with its single src homology 2 (SH2) domain architecture and unusually high affinity for SLAM (also called CD150), SAP has been suggested to function by blocking binding of SHP-2 or other SH2-containing signalling proteins to SLAM receptors. Additionally, SAP has recently been shown to be required for recruitment and activation of the Src-family kinase FynT after SLAM ligation. This signalling 'adaptor' function has been difficult to conceptualize, because unlike typical SH2-adaptor proteins, SAP contains only a single SH2 domain and lacks other recognized protein interaction domains or motifs. Here, we show that the SAP SH2 domain binds to the SH3 domain of FynT and directly couples FynT to SLAM. The crystal structure of a ternary SLAM-SAP-Fyn-SH3 complex reveals that SAP binds the FynT SH3 domain through a surface-surface interaction that does not involve canonical SH3 or SH2 binding interactions. The observed mode of binding to the Fyn-SH3 domain is expected to preclude the auto-inhibited conformation of Fyn, thereby promoting activation of the kinase after recruitment. These findings broaden our understanding of the functional repertoire of SH3 and SH2 domains.