The main objectives of bioprocesses are to reliably deliver drugs in a relatively short time frame with high quality within a tight regulatory framework. Bioprocesses are highly complex, the level of automation is moderate, and there is constant pressure to improve efficiency and costs. In addition, climate change and resource scarcity mandate a reduction in the environmental footprint of bioprocesses and production facilities. In the biopharmaceutical industry, two extreme production scenarios are applied: a fully disposable factory with the characteristics of full flexibility and speed, or a fixed large-scale plant with high capacity. Forward-looking solutions and ideas will be discussed how to combine new processes and environmental friendliness for the benefit of the patient, security of supply and profitability. The concept will be extended to large scale production of proteins for food and non-pharma applications, e.g., in material science and a roadmap towards a future plant will be laid out.

@en

Optimizing a biopharmaceutical chromatographic purification process is currently the greatest challenge during process development. A lack of process understanding calls for extensive experimental efforts in pursuit of an optimal process. In silico techniques, such as mechanistic or data driven modeling, enhance the understanding, allowing more cost-effective and time efficient process optimization. This work presents a modeling strategy integrating quantitative structure property relationship (QSPR) models and chromatographic mechanistic models (MM) to optimize a cation exchange (CEX) capture step, limiting experiments. In QSPR, structural characteristics obtained from the protein structure are used to describe physicochemical behavior. This QSPR information can be applied in MM to predict the chromatogram and optimize the entire process. To validate this approach, retention profiles of six proteins were determined experimentally from mixtures, at different pH (3.5, 4.3, 5.0, and 7.0). Four proteins at different pH's were used to train QSPR models predicting the retention volumes and characteristic charge, subsequently the equilibrium constant was determined. For an unseen protein knowing only the protein structure, the retention peak difference between the modeled and experimental peaks was 0.2% relative to the gradient length (60 column volume). Next, the CEX capture step was optimized, demonstrating a consistent result in both the experimental and QSPR-based methods. The impact of model parameter confidence on the final optimization revealed two viable process conditions, one of which is similar to the optimization achieved using experimentally obtained parameters. The multiscale modeling approach reduces the required experimental effort by identification of initial process conditions, which can be optimized.@en

Purification of recombinantly produced biopharmaceuticals involves removal of host cell material, such as host cell proteins (HCPs). For lysates of the common expression host Escherichia coli (E. coli) over 1500 unique proteins can be identified. Currently, understanding the behavior of individual HCPs for purification operations, such as preparative chromatography, is limited. Therefore, we aim to elucidate the elution behavior of individual HCPs from E. coli strain BLR(DE3) during chromatography. Understanding this complex mixture and knowing the chromatographic behavior of each individual HCP improves the ability for rational purification process design. Specifically, linear gradient experiments were performed using ion exchange (IEX) and hydrophobic interaction chromatography, coupled with mass spectrometry-based proteomics to map the retention of individual HCPs. We combined knowledge of protein location, function, and interaction available in literature to identify trends in elution behavior. Additionally, quantitative structure–property relationship models were trained relating the protein 3D structure to elution behavior during IEX. For the complete data set a model with a cross-validated R2 of 0.55 was constructed, that could be improved to a R2 of 0.70 by considering only monomeric proteins. Ultimately this study is a significant step toward greater process understanding.@en

Mechanistic models mostly focus on the target protein and some selected process- or product-related impurities. For a better process understanding, however, it is advantageous to describe also reoccurring host cell protein impurities. Within the purification of biopharmaceuticals, the binding of host cell proteins to a chromatographic resin is far from being described comprehensively. For a broader coverage of the binding characteristics, large-scale proteomic data and systems level knowledge on protein interactions are key. However, a method for determining binding parameters of the entire host cell proteome to selected chromatography resins is still lacking. In this work, we have developed a method to determine binding parameters of all detected individual host cell proteins in an Escherichia coli harvest sample from large-scale proteomics experiments. The developed method was demonstrated to model abundant and problematic proteins, which are crucial impurities to be removed. For these 15 proteins covering varying concentration ranges, the model predicts the independently measured retention time during the validation gradient well. Finally, we optimized the anion exchange chromatography capture step in silico using the determined isotherm parameters of the persistent host cell protein contaminants. From these results, strategies can be developed to separate abundant and problematic impurities from the target antigen.@en

The monoclonal antibody (mAb) industry is becoming increasingly digitalized. Digital twins are becoming increasingly important to test or validate processes before manufacturing. High-Throughput Process Development (HTPD) has been progressively used as a tool for process development and innovation. The combination of High-Throughput Screening with fast computational methods allows to study processes in-silico in a fast and efficient manner. This paper presents a hybrid approach for HTPD where equal importance is given to experimental, computational and decision-making stages. Equilibrium adsorption isotherms of 13 protein A and 16 Cation-Exchange resins were determined with pure mAb. The influence of other components in the clarified cell culture supernatant (harvest) has been under-investigated. This work contributes with a methodology for the study of equilibrium adsorption of mAb in harvest to different protein A resins and compares the adsorption behavior with the pure sample experiments. Column chromatography was modelled using a Lumped Kinetic Model, with an overall mass transfer coefficient parameter (k_ov). The screening results showed that the harvest solution had virtually no influence on the adsorption behavior of mAb to the different protein A resins tested. k_ov was found to have a linear correlation with the sample feed concentration, which is in line with mass transfer theory. The hybrid approach for HTPD presented highlights the roles of the computational, experimental, and decision-making stages in process development, and how it can be implemented to develop a chromatographic process. The proposed white-box digital twin helps to accelerate chromatographic process development.

@en

The challenging task of designing biopharmaceutical downstream processes is initially to select the type of unit operations, followed by optimizing their operating conditions. For complex flowsheet optimizations, the strategy becomes crucial in terms of duration and outcome. In this study, we compared three optimization strategies, namely, simultaneous, top-to-bottom, and superstructure decomposition. Moreover, all strategies were evaluated by either using chromatographic Mechanistic Models (MMs) or Artificial Neural Networks (ANNs). An overall evaluation of 39 flowsheets was performed, including a buffer-exchange step between the chromatography operations. All strategies identified orthogonal structures to be optimal, and the weighted overall performance values were generally consistent between the MMs and ANNs. In terms of time-efficiency, the decomposition method with MMs stands out when utilizing multiple cores on a multiprocessing system for simulations. This study analyses the influence of different optimization strategies on flowsheet optimization and advices on suitable strategies and modeling techniques for specific scenarios.@en

Protein-based biopharmaceuticals require high purity before final formulation to ensure product safety, making process development time consuming. Implementation of computational approaches at the initial stages of process development offers a significant reduction in development efforts. By preselecting process conditions, experimental screening can be limited to only a subset. One such computational selection approach is the application of Quantitative Structure Property Relationship (QSPR) models that describe the properties exploited during purification. This work presents a novel open-source Python tool capable of extracting a range of features from protein 3D models on a local computer allowing total transparency of the calculations. As open-source tool, it also impacts initial investments in constructing a QSPR workflow for protein property prediction for third parties, making it widely applicable within the field of bioprocess development. The focus of current calculated molecular features is projection onto the protein surface by constructing surface grid representations. Linear regression models were trained with the calculated features to predict chromatographic retention times/volumes. Model validation shows a high accuracy for anion and cation exchange chromatography data (cross-validated R² of 0.87 and 0.95). Hence, these models demonstrate the potential of the use of QSPR to accelerate process design.

@en

The hydrodynamics of the Expanded Bed Adsorption process is studied through simulations combining Computational Fluid Dynamics and the Discrete Element Method. A representative base case is defined, based on process design parameters commonly encountered in literature. Then, 19 other cases are defined, each representing a singular adjustment to the column design, material properties, or operating conditions. The parameters that are varied are the expansion factor, liquid viscosity, bed aspect ratio, mean particle density, width of the particle density distribution, width of the particle size distribution, column taper angle, and column alignment angle. The impact of each adjustment on the bed behaviour is discussed, using the local particle size distribution and solids dispersion coefficient as main indicators of bed stability. Optimal performance was found for an expansion factor of two to three, and the combination of particle size distribution and particle density distribution was found to greatly improve bed stability. The mixing process of the liquid and solid phases is concluded to be of highly complex nature, and cannot simply be predicted from the liquid flow velocity.

@en

An optimal purification process for biopharmaceutical products is important to meet strict safety regulations, and for economic benefits. To find the global optimum, it is desirable to screen the overall design space. Advanced model-based approaches enable to screen a broad range of the design-space, in contrast to traditional statistical or heuristic-based approaches. Though, chromatographic mechanistic modeling (MM), one of the advanced model-based approaches, can be speed-limiting for flowsheet optimization, which evaluates every purification possibility (e.g., type and order of purification techniques, and their operating conditions). Therefore, we propose to use artificial neural networks (ANNs) during global optimization to select the most optimal flowsheets. So, the number of flowsheets for final local optimization is reduced and consequently the overall optimization time. Employing ANNs during global optimization proved to reduce the number of flowsheets from 15 to only 3. From these three, one flowsheet was optimized locally and similar final results were found when using the global outcome of either the ANN or MM as starting condition. Moreover, the overall flowsheet optimization time was reduced by 50% when using ANNs during global optimization. This approach accelerates the early purification process design; moreover, it is generic, flexible, and regardless of sample material's type.

@en

The ultrafiltration/diafiltration (UF/DF) process is a crucial step in the canola protein isolation process from rapeseed meal. The process involves using a multi-stage membrane system to separate components of the mixture. As diafiltration dilutes the feed stream in the ultrafiltration system, a large amount of diafiltration water is required. Reducing the diafiltrate for sustainability reasons can be done by carefully selecting process variables or using recycle streams. However, finding the optimum process variables can be a meticulous process if performed experimentally or via trial and error. In this study, we performed an optimization using a mechanistic model integrated with a genetic algorithm to aid in finding an optimum combination of process variables. The mechanistic ultrafiltration model was derived by taking into account transport phenomena within the filtration system. Parameters were characterized experimentally in term of viscosity coefficient, membrane resistance, cake porosity, aggregate diameter equivalent, and material compressibility factors. Using the mechanistic model-based optimization in combination with actual experimental values, the performance of a four-stage UF/DF system could already be improved despite fixing the configuration, albeit at the cost of a reduction in purity. Further improvement could be achieved by using recycle streams. The optimized system achieved a diafiltrate reduction of up to 79%, an increase of purity of up to 31%, and an increase of dry matter content of up to 18%, while maintaining the product purity of the reference set-up.

@en

Mass-spectrometry-based proteomics is increasingly employed to monitor purification processes or to detect critical host cell proteins in the final drug substance. This approach is inherently unbiased and can be used to identify individual host cell proteins without prior knowledge. In process development for the purification of new biopharmaceuticals, such as protein subunit vaccines, a broader knowledge of the host cell proteome could promote a more rational process design. Proteomics can establish qualitative and quantitative information on the complete host cell proteome before purification (i.e., protein abundances and physicochemical properties). Such information allows for a more rational design of the purification strategy and accelerates purification process development. In this study, we present an extensive proteomic characterisation of two E. coli host cell strains widely employed in academia and industry to produce therapeutic proteins, BLR and HMS174. The established database contains the observed abundance of each identified protein, information relating to their hydrophobicity, the isoelectric point, molecular weight, and toxicity. These physicochemical properties were plotted on proteome property maps to showcase the selection of suitable purification strategies. Furthermore, sequence alignment allowed integration of subunit information and occurrences of post-translational modifications from the well-studied E. coli K12 strain.

@en

Parkinson's Disease (PD) is a common neurodegenerative disorder affecting millions of people worldwide for which there are only symptomatic therapies. Small molecules able to target key pathological processes in PD have emerged as interesting options for modifying disease progression. We have previously shown that a (poly)phenol-enriched fraction (PEF) of Corema album L. leaf extract modulates central events in PD pathogenesis, namely α-synuclein (αSyn) toxicity, aggregation and clearance. PEF was now subjected to a bio-guided fractionation with the aim of identifying the critical bioactive compound. We identified genipin, an iridoid, which relieves αSyn toxicity and aggregation. Furthermore, genipin promotes metabolic alterations and modulates lipid storage and endocytosis. Importantly, genipin was able to prevent the motor deficits caused by the overexpression of αSyn in a Drosophila melanogaster model of PD. These findings widens the possibility for the exploitation of genipin for PD therapeutics.

@en

The lack of process analytical technologies able to provide real-time information and process control over a biopharmaceutical process has long impaired the transition to continuous biomanufacturing. For the monoclonal antibody (mAb) production, aggregate formation is a major critical quality attribute (CQA) with several known process parameters (i.e., protein concentration and agitation) influencing this phenomenon. The development of a real-time tool to monitor aggregate formation is then crucial to gain control and achieve a continuous processing. Due to an inherent short operation time, miniaturized biosensors placed after each step can be a powerful solution. In this work, the development of a fluorescent dye-based microfluidic sensor for fast at-line PAT is described, using fluorescent dyes to examine possible mAb size differences. A zigzag microchannel, which provides 90% of mixing efficiency under 30 s, coupled to an UV–Vis detector, and using four FDs, was studied and validated. With different generated mAb aggregation samples, the FDs Bis-ANS and CCVJ were able to robustly detect from, at least, 2.5% to 10% of aggregation. The proposed FD-based micromixer is then ultimately implemented and validated in a lab-scale purification system, demonstrating the potential of a miniaturized biosensor to speed up CQAs measurement in a continuous process.

@en

Membrane technology is commonly used within food, bio- and pharmaceutical processes. Beside single-stage membranes, multi-stage membrane systems are become more popular to improve separation performance. In this review, we present a unified four-phase model-based optimization framework to optimize these systems, using mechanistic models, empirical models including machine learning models, or a combination of them. We begin by providing a general overview and outlining the steps to construct each phase in the framework. The importance of each stage and critical points to consider are discussed. We then provide detailed information for each phase, including the governing equations from known literature models. Finally, we explore the platform's potential applications and outlook. Despite the great potential of an integrated approach, studies thus far focus either on extensive membrane modeling with brute-force optimization via simple comparison or on meticulous optimization using an oversimplified membrane model. We believe that the integrated framework can bridge the well justified approaches in both filtration modeling and mathematical optimization and help in designing multi-unit processes.

@en

The implementation of continuous processing in the biopharmaceutical industry is hindered by the scarcity of process analytical technologies (PAT). To monitor and control a continuous process, PAT tools will be crucial to measure real-time product quality attributes such as protein aggregation. Miniaturizing these analytical techniques can increase measurement speed and enable faster decision-making. A fluorescent dye (FD)-based miniaturized sensor has previously been developed: a zigzag microchannel which mixes two streams under 30 s. Bis-ANS and CCVJ, two established FDs, were employed in this micromixer to detect aggregation of the biopharmaceutical monoclonal antibody (mAb). Both FDs were able to robustly detect aggregation levels starting at 2.5%. However, the real-time measurement provided by the microfluidic sensor still needs to be implemented and assessed in an integrated continuous downstream process. In this work, the micromixer is implemented in a lab-scale integrated system for the purification of mAbs, established in an ÄKTA™ unit. A viral inactivation and two polishing steps were reproduced, sending a sample of the product pool after each phase directly to the microfluidic sensor for aggregate detection. An additional UV sensor was connected after the micromixer and an increase in its signal would indicate that aggregates were present in the sample. The at-line miniaturized PAT tool provides a fast aggregation measurement, under 10 min, enabling better process understanding and control.

@en

Lubrication forces play a major role in the behaviour of fluid–solid systems, where they affect the collisions between particles. Current implementations of lubrication forces in unresolved simulations often suffer from shortcomings, such as neglecting parts of the physics or relying on arbitrarily defined parameters. In this short communication, we propose a novel implementation, rigorously defined based on physical and numerical factors. Both particle roughness and deformation are considered, and the model accuracy is demonstrated through comparison with experimental results.

@en

A major challenge in the transition to continuous biomanufacturing is the lack of process analytical technology (PAT) tools which are able to collect real-time information on the process and elicit a response to facilitate control. One of the critical quality attributes (CQAs) of interest during monoclonal antibodies production is aggregate formation. The development of a real-time PAT tool to monitor aggregate formation is then crucial to have immediate feedback and process control. Miniaturized sensors placed after each unit operation can be a powerful solution to speed up an analytical measurement due to their characteristic short reaction time. In this work, a micromixer structure capable of mixing two streams is presented, to be employed in the detection of mAb aggregates using fluorescent dyes. Computational fluid dynamics (CFD) simulations were used to compare the mixing performance of a series of the proposed designs. A final design of a zigzag microchannel with 45° angle was reached and this structure was subsequently fabricated and experimentally validated with colour dyes and, later, with a FITC-IgG molecule. The designed zigzag micromixer presents a mixing index of around 90%, obtained in less than 30 seconds. Therefore, a micromixer channel capable of a fast and efficient mixing is hereby demonstrated, to be used as a real-time PAT tool for a fluorescence based detection of protein aggregation.

@en

Biopharmaceuticals are becoming increasingly important in modern healthcare. Monoclonal antibodies (mAb) are one of the most widely used therapeutic proteins and are important for the treatment of cancer and autoimmune diseases, among others. After cell culture there are still large amounts of other impurities (e.g. host cell proteins) in solution. Chromatography is usually the first purification step, allowing to increase purity and reduce volume. This comes associated with high costs and chromatography accounts for a significant portion of total production costs for therapeutic proteins. Chromatographic process development may be time consuming and use large amounts of resins. Therefore, there is increased interest in finding cheaper techniques for chromatographic process development without compromising accuracy. This paper presents a highly sophisticated microfluidic chip approach for efficient adsorption isotherm determinations compared to current chromatographic process development. Implementation of an image analysis software ensures that chromatographic resin volume is accurately determined. The adsorption isotherm performance of microfluidics was compared to the robotic Liquid-handling Station (LHS) and labor intensive Eppendorf tubes. The microfluidic chip allows a 15-fold volume reduction and resin consumptions as low as 100/200 nl (200/100-fold reduction). The microfluidic chip performed comparably to the other miniaturized techniques, using less liquid and resin volume. For process development of expensive products (e.g. monoclonal antibodies), miniaturization (provided by the microfluidic chip) proved to be the most cost effective alternative whereas for less valuable products (e.g. lysozyme) automation (provided by the LHS) was the most cost effective alternative.

@en

The safety requirements for vaccines are extremely high since they are administered to healthy people. For that reason, vaccine development is time-consuming and very expensive. Reducing time-to-market is key for pharmaceutical companies, saving lives and money. Therefore the need is raised for systematic, general and efficient process development strategies to shorten development times and enhance process understanding. High throughput technologies tremendously increased the volume of process-related data available and, combined with statistical and mechanistic modeling, new high throughput process development (HTPD) approaches evolved. The introduction of model-based HTPD enabled faster and broader screening of conditions, and furthermore increased knowledge. Model-based HTPD has particularly been important for chromatography, which is a crucial separation technique to attain high purities. This review provides an overview of downstream process development strategies and tools used within the (bio)pharmaceutical industry, focusing attention on (protein subunit) vaccine purification processes. Subsequently high throughput process development and other combinatorial approaches are discussed and compared according to their experimental effort and understanding. Within a growing sea of information, novel modeling tools and artificial intelligence (AI) gain importance for finding patterns behind the data and thereby acquiring a deeper process understanding.

@en

Efficient recovery and utilization of valuable components from industrial food side streams is a main driver towards a circular economy. Among different available purification techniques, adsorption can effectively recover these components. However, the conventional batch mode of operation can limit its applicability in food processes due to limited efficiency. This work compares conventional batch packed bed adsorption with semi-continuous adsorption (so-called CaptureSMB) for the recovery of sinapic acid at industrial scale, using a food grade resin Amberlite^TM FPX66. A mathematical mechanistic model able to describe semi-continuous operation is successfully validated and used to identify optimum operating parameters to maximize productivity and resin capacity utilization in batch and semi-continuous operating modes. The results indicate that CaptureSMB outperforms batch operation, increasing productivity from 5.18 g/L/h to 10.3 g/L/h for a given yield (>97%). A resin capacity utilization (RU) of around 70% is observed in both operating modes when productivity is maximized. A 92% RU can be accomplished for a given yield using the CaptureSMB process at a productivity of 7.0 g/L/h, higher than for conventional batch operation. The use of semi-continuous adsorption operation in food industry contributes to more efficient processes at reduced purification costs.

@en