Comparative modeling of breakthrough curves of bovine serum albumin in anion-exchange chromatography

On the apparent dispersion coefficient of the equilibrium dispersion model: An asymptotic analysis

To model chromatography, researchers have developed several approaches. These cover a broad range of applications and, depending on the assumptions adopted, have different levels of accuracy. In general, the most suitable modelling approach is the simplest that can describe a process with the desired accuracy. A model that often meets this criterion is the equilibrium dispersion model (EDM). This features one mass balance equation per analyte, including an axial dispersion term, and assumes the analyte concentrations in the mobile and stationary phases to be in local equilibrium. To account for the finite mass transfer rate between the phases, the model employs an apparent dispersion coefficient. Two expressions are available for this coefficient, one being used much more frequently than the other. In this paper, we aimed to clarify which one should be favoured. A desirable feature of simple models is that they can be derived from more general ones with appropriate physical assumptions and rigorous mathematical methods. Thus, to answer our research question, we derived the EDM from the more general pore diffusion model (POR), using an asymptotic method. The expression obtained for the apparent dispersion coefficient does agree with one of the two reported in the literature - the less frequently used. To test the validity of this expression, we simulated elution profiles using the two versions of the EDM and compared the results against those from the POR model. The simulations were conducted in the range where the POR and EDM models should be essentially equivalent, their results confirming the outcome of the asymptotic analysis. This work offers a solid theoretical grounding for the EDM, clarifies which formulation of the model is correct, and provides usable applicability conditions for the model.

A comprehensive method for modeling and simulating ion exchange chromatography of complex mixtures

In recent years, many biological-based products have been developed, representing a significant fraction of income in the pharmaceutical market. Ion exchange chromatography is an important downstream step for the purification of target recombinant proteins present in clarified cell extracts, together with many other unknown impurities. This work develops a robust approach to model and simulate the purification of untagged heterologous proteins, so that the improved conditions to carry out an ion exchange chromatography are identified in a rational basis prior to the real purification run itself. Purification of the pneumococcal surface protein A (PspA4Pro) was used as a case study. This protein is produced by recombinant Escherichia coli and is a candidate for the manufacture of improved pneumococcal vaccines. The developed method combined experimental and computational procedures. Different anion exchange operating conditions were mapped in order to gather a broad range of representative experimental data. The equilibrium dispersive and the steric mass action equations were used to model and simulate the process. A training strategy to fit the model and separately describe the elution profiles of PspA4Pro and other proteins of the cell extract was applied. Based on the simulation results, a reduced ionic strength was applied for PspA4Pro elution, leading to increases of 14.9% and 11.5% for PspA4Pro recovery and purity, respectively, compared to the original elution profile. These results showed the potential of this method, which could be further applied to improve the performance of ion exchange chromatography in the purification of other target proteins under real process conditions.

Numerical Approximation of the Nonequilibrium Model of Gradient Elution Chromatography Considering Linear and Nonlinear Solvent Strength Models

In both linear and nonlinear chromatography, the lumped kinetic model is a suitable model for predicting elution bands when appropriate equilibrium functions and mass transfer coefficients are accessible. This model also works well in the case of gradient elution chromatography if variations in the equilibrium functions due to changes in the mobile phase composition are known. The rational selection of an optimum gradient is explored in this study from three different perspectives using the lumped kinetic model. Elution profiles generated by using (a) linear solvent strength, (b) quadratic solvent strength, and (c) power law are investigated. The effectiveness and reliability of the suggested numerical approach, utilizing the flux-limiting finite volume method, are demonstrated through numerical simulations. The impacts of axial dispersion, nonlinearity coefficient, Henry's constant, mass transfer coefficient, and gradient parameters are studied on single and two-component elution profiles.

An Experimental and Modeling Combined Approach in Preparative Hydrophobic Interaction Chromatography

Chromatography is a technique widely used in the purification of biopharmaceuticals, and generally consists of several chromatographic steps. In this work, Hydrophobic Interaction Chromatography (HIC) is investigated as a polishing step for the purification of therapeutic proteins. Adsorption mechanisms in hydrophobic interaction chromatography are still not completely clear and a limited amount of published data is available. In addition to new data on adsorption isotherms for some proteins (obtained both by high-throughput and frontal analysis method), and a comparison of different models proposed in the literature, two different approaches are compared in this work to investigate HIC. The predictive approach exploits an in-house code that simulates the behavior of the component in the column using the model parameters found from the fitting of experimental data. The estimation approach, on the other hand, exploits commercial software in which the model parameters are found by the fitting of a few experimental chromatograms. The two approaches are validated on some bind-elute runs: the predictive approach is very informative, but the experimental effort needed is high; the estimation approach is more effective, but the knowledge gained is lower. The second approach is also applied to an in-development industrial purification process and successfully resulted in predicting the behavior of the system, allowing for optimization with a reduction in the time and amount of sample needed.

Optimal performance comparison of the simulated moving bed process variants based on the modulation of the length of zones and the feed concentration

The VariCol and ModiCon processes are two variants of the simulated moving bed (SMB) process, characterized by the modulation of the length of zones of the chromatographic column train and the feed concentration. These features give more flexibility than the conventional operation, leading to essential improvements in the separation and purification of mixtures. The optimal performance comparison of these two variants, the hybrid formed by their combination, and the conventional SMB process are scarce in the literature. This comparison helps discover new characteristics of each single and combined operation mode and creates guidelines to select the appropriate operation mode for possible real applications. In this work, the performance comparison of the ModiCon, VariCol, ModiCon+VariCol, and SMB processes is carried out in terms of maximal throughput for specific product purity values. Particular emphasis is placed on both the ModiCon and the hybrid ModiCon+VariCol processes characteristics. A strategy for combining and optimizing the ModiCon and the VariCol processes was determined. As a case study, the enantioseparation of guaifenesin was considered. In the ModiCon process, more than two modulation subintervals did not improve the performance in the separation. The optimal pattern, based on two subintervals, has zero feed concentration in the first subinterval and the maximal concentration in the second one. The best result for the hybrid operation (ModiCon+VariCol) was reached when the feed port moves simultaneously as the SMB process switching period. The optimal throughput of the ModiCon and the ModiCon+VariCol processes was almost doubled than that of the SMB process. These performances were based on larger zones I and II and not in zones II and III as occur with the SMB and VariCol process. The throughput in the hybrid operation increases more significantly than the ModiCon process when 5 columns were considered instead of 6. The hybrid operation could be more attractive for a system with a few numbers of columns.

Chromatography bioseparation technologies and in-silico modelings for continuous production of biotherapeutics

The potential of continuous bioprocessing is hindered by the bottlenecks of chromatography processing, which continues to be executed in batch mode. Highlighting the critical drawbacks of batch chromatography, this review underscores the transition that the industry has made by implementing continuous upstream process without devising a working model for downstream chromatography operations. Even though multitude of process development initiatives have commenced, the review emphasizes the first principle models of chromatography on which these initiatives are built. Various models of continuous chromatography, which are essential, but not limited to multi-column systems, employed to congeal a unified process are reviewed. Advancements made by several mechanistic models and simulations to maximize productivity and performance are described, in an attempt to provide the integral tools. The modeling tools can be used for development of a strong model based control strategy and can be embedded into the continuous chromatography framework. The review addresses the limitations and challenges of the current modeling methods for development of robust mechanistic modeling and efficient unit operation platform in continuous chromatography.

Optimizing Column Length and Particle Size in Preparative Batch Chromatography Using Enantiomeric Separations of Omeprazole and Etiracetam as Models: Feasibility of Taguchi Empirical Optimization

The overreaching purpose of this study is to evaluate new approaches for determining the optimal operational and column conditions in chromatography laboratories, i.e., how best to select a packing material of proper particle size and how to determine the proper length of the column bed after selecting particle size. As model compounds, we chose two chiral drugs for preparative separation: omeprazole and etiracetam. In each case, two maximum allowed pressure drops were assumed: 80 and 200 bar. The processes were numerically optimized (mechanistic modeling) with a general rate model using a global optimization method. The numerical predictions were experimentally verified at both analytical and pilot scales. The lower allowed pressure drop represents the use of standard equipment, while the higher allowed drop represents more modern equipment. For both compounds, maximum productivity was achieved using short columns packed with small-particle size packing materials. Increasing the allowed backpressure in the separation leads to an increased productivity and reduced solvent consumption. As advanced numerical calculations might not be available in the laboratory, we also investigated a statistically based approach, i.e., the Taguchi method (empirical modeling), for finding the optimal decision variables and compared it with advanced mechanistic modeling. The Taguchi method predicted that shorter columns packed with smaller particles would be preferred over longer columns packed with larger particles. We conclude that the simpler optimization tool, i.e., the Taguchi method, can be used to obtain "good enough" preparative separations, though for accurate processes, optimization, and to determine optimal operational conditions, classical numerical optimization is still necessary.

Accounting for host cell protein behavior in anion-exchange chromatography

Host cell proteins (HCP) are a problematic set of impurities in downstream processing (DSP) as they behave most similarly to the target protein during separation. Approaching DSP with the knowledge of HCP separation behavior would be beneficial for the production of high purity recombinant biologics. Therefore, this work was aimed at characterizing the separation behavior of complex mixtures of HCP during a commonly used method: anion-exchange chromatography (AEX). An additional goal was to evaluate the performance of a statistical methodology, based on the characterization data, as a tool for predicting protein separation behavior. Aqueous two-phase partitioning followed by two-dimensional electrophoresis provided data on the three physicochemical properties most commonly exploited during DSP for each HCP: pI (isoelectric point), molecular weight, and surface hydrophobicity. The protein separation behaviors of two alternative expression host extracts (corn germ and E. coli) were characterized. A multivariate random forest (MVRF) statistical methodology was then applied to the database of characterized proteins creating a tool for predicting the AEX behavior of a mixture of proteins. The accuracy of the MVRF method was determined by calculating a root mean squared error value for each database. This measure never exceeded a value of 0.045 (fraction of protein populating each of the multiple separation fractions) for AEX. © 2016 American Institute of Chemical Engineers Biotechnol. Prog., 32:1453-1463, 2016.

Purification of Monoclonal Antibodies Using a Fiber Based Cation-Exchange Stationary Phase: Parameter Determination and Modeling

Monoclonal antibodies (mAb) currently dominate the market for protein therapeutics. Because chromatography unit operations are critical for the purification of therapeutic proteins, the process integration of novel chromatographic stationary phases, driven by the demand for more economic process schemes, is a field of ongoing research. Within this study it was demonstrated that the description and prediction of mAb purification on a novel fiber based cation-exchange stationary phase can be achieved using a physico-chemical model. All relevant mass-transport phenomena during a bind and elute chromatographic cycle, namely convection, axial dispersion, boundary layer mass-transfer, and the salt dependent binding behavior in the fiber bed were described. This work highlights the combination of model adaption, simulation, and experimental parameter determination through separate measurements, correlations, or geometric considerations, independent from the chromatographic cycle. The salt dependent binding behavior of a purified mAb was determined by the measurement of adsorption isotherms using batch adsorption experiments. Utilizing a combination of size exclusion and protein A chromatography as analytic techniques, this approach can be extended to a cell culture broth, describing the salt dependent binding behavior of multiple components. Model testing and validation was performed with experimental bind and elute cycles using purified mAb as well as a clarified cell culture broth. A comparison between model calculations and experimental data showed a good agreement. The influence of the model parameters is discussed in detail.

Enantioseparation of racemic aminoglutethimide using asynchronous simulated moving bed chromatography

The separation of aminoglutethimide enantiomers by the continuous multicolumn chromatographic processes were investigated experimentally and theoretically, where the columns were packed with cellulose tris 3,5-dimethylphenyl-carbamate stationary phase (brand name Chiralcel OD) and mobile phase was a mixture of n-hexane and ethanol with monoethanolamine additive. The continuous enantioseparation processes included a synchronous shifting process (SMB) and an asynchronous shifting process (VARICOL), which allowed reducing the column number (here from six-column SMB to five-column VARICOL process). Transport-dispersive model with the consideration of both intraparticle mass transfer resistance and axial dispersion was adopted to design and optimize the operation conditions for the separation of aminoglutethimide enantiomers by SMB process and VARICOL process. According to the optimized operation conditions, experiments were carried out on VARICOL-Micro unit using five-column VARICOL process with 1/1.5/1.5/1 configuration and six-column SMB process with 1/2/2/1 configuration. Products of R-aminoglutethimide (R-AG) enantiomer and S-aminoglutethimide (S-AG) enantiomer with more than 99.0% purity were obtained continuously from extract stream and raffinate stream, respectively. Furthermore, the experiemntal data obtained from five-column VARICOL process were compared with that from six-column SMB process, the feasibility and efficiency for the separation of guaifenesin enantiomers by VARICOL processes were evaluated.

Combining Synthetic and Analytical Strategies for Preparative HPLC Enantioseparation of Monastrol Racemic Mixture

Large-scale resolution of racemic monastrol has been carried out by normal-phase mode HPLC on an amylose-based chiral stationary phase. Because monastrol solubility, in media of proper compositions for normal-mode HPLC separations (in terms of retention factors and selectivity), was significantly low and impractical for preparative scopes, racemic monastrol was transformed into the corresponding O-tert-butyldimethylsilyl derivative. The tert-butyldimethylsilyl group was chosen as a suitable derivatizing agent because it induced approximatively a six-times higher solubility and allowed for an almost quantitative recovery of pure monastrol from the derivatization-deprotection sequence. The competitive isotherms of the O-tert-butyldimethylsilyl compounds, measured through frontal analysis, were fitted to competitive Langmuir and four-parameter bi-Langmuir models. The equilibrium dispersive model of chromatography was used for modeling the nonlinear separation of the racemate and to optimize the experimental conditions for collection of highly concentrated fractions of pure (R,S)-O-tert-butyldimethylsilyl compounds, from which significant amounts of the corresponding enantiomers of monastrol (about 100 mg of each enantiomer with 30 runs on an analytical-scale column) were obtained by quantitative back-derivatization.