Chromatography Analysis and Design Toolkit (CADET)

Abstract CADET is an open source modeling and simulation framework for column liquid chromatography. The software is freely distributed to both academia and industry under the GPL license ( http://github.com/modsim/cadet ). CADET is based on a core simulator that is written in object oriented C++ and applies modern mathematical algorithms for efficiently solving a variety of customary chromatography models. This simulation engine is interfaced to a suite of MATLAB tools for setting up and executing the most common scientific workflows, e.g., model calibration, process design, robustness analysis, statistical analysis, and experimental design. The model library and numerical methods are continuously extended and improved. For instance, binding models with multiple bound states, pH and/or temperature dependence of binding parameters, surface diffusion, and arbitrary spacing of the radial discretization have been recently added. Moreover, numerical accuracy and computational speed of the code are comprehensively benchmarked using high precision reference solutions and realistic model problems. Versatility of the CADET modeling platform is demonstrated with several examples that are also published as open source code and can be freely adapted to specific use cases. In one of several case studies, sequential and simultaneous optimization of elution gradient shape and cut times are compared for a three component separation. This process is designed to achieve Pareto optimal purity and yield of the central fraction. Moreover, the robustness of these designs with respect to typical process variations is systematically studied. The last case study illustrates the optimal design of experiments for estimating model parameters with maximal accuracy.

[1]  A Kremling,et al.  Exploiting the bootstrap method for quantifying parameter confidence intervals in dynamical systems. , 2006, Metabolic engineering.

[2]  J. Åkesson,et al.  Methods and Tools for Robust Optimal Control of Batch Chromatographic Separation Processes , 2015 .

[3]  S. Osher,et al.  Weighted essentially non-oscillatory schemes , 1994 .

[4]  Jürgen Hubbuch,et al.  Advances in inline quantification of co‐eluting proteins in chromatography: Process‐data‐based model calibration and application towards real‐life separation issues , 2015, Biotechnology and bioengineering.

[5]  New liquid chromatography method combining thermo-responsive material and inductive heating via alternating magnetic field. , 2008, Journal of chromatography. B, Analytical technologies in the biomedical and life sciences.

[6]  How Much Can We Learn from a Single Chromatographic Experiment? A Bayesian Perspective. , 2016, Analytical chemistry.

[7]  Karin Westerberg,et al.  Model‐Based Process Challenge of an Industrial Ion‐Exchange Chromatography Step , 2012 .

[8]  Jürgen Hubbuch,et al.  Detection, Quantification, and Propagation of Uncertainty in High‐Throughput Experimentation by Monte Carlo Methods , 2012 .

[9]  Niklas Andersson,et al.  Supporting Design and Control of a Reversed‐Phase Chromatography Step by Mechanistic Modeling , 2012 .

[10]  Massimo Morbidelli,et al.  Model based adaptive control of a continuous capture process for monoclonal antibodies production. , 2016, Journal of chromatography. A.

[11]  M. E. Pate,et al.  Principal Component Analysis of Nonlinear Chromatography , 2008, Biotechnology progress.

[12]  Jürgen Hubbuch,et al.  Model-based integrated optimization and evaluation of a multi-step ion exchange chromatography , 2014 .

[13]  Massimiliano Barolo,et al.  Online Model-Based Redesign of Experiments for Parameter Estimation in Dynamic Systems , 2009 .

[14]  Jürgen Hubbuch,et al.  Effects of ionic strength and mobile phase pH on the binding orientation of lysozyme on different ion-exchange adsorbents. , 2008, Journal of chromatography. A.

[15]  Antônio José da Silva Neto,et al.  Uncertainty Quantification in Chromatography Process Identification Based on Markov Chain Monte Carlo , 2016 .

[16]  P. Hewson Bayesian Data Analysis 3rd edn A. Gelman, J. B. Carlin, H. S. Stern, D. B. Dunson, A. Vehtari and D. B. Rubin, 2013 Boca Raton, Chapman and Hall–CRC 676 pp., £44.99 ISBN 1‐439‐84095‐4 , 2015 .

[17]  Marcel Ottens,et al.  Ion-exchange chromatographic protein refolding. , 2010, Journal of chromatography. A.

[18]  Marco Mazzotti,et al.  Local Equilibrium Theory for the Binary Chromatography of Species Subject to a Generalized Langmuir Isotherm , 2006 .

[19]  Eric von Lieres,et al.  Fast arbitrary order moments and arbitrary precision solution of the general rate model of column liquid chromatography with linear isotherm , 2016, Comput. Chem. Eng..

[20]  P. Marquis,et al.  Effect of pH on protein adsorption to hydroxyapatite and tricalcium phosphate ceramics. , 1997, Biomaterials.

[21]  Nien-Hwa Linda Wang,et al.  Pore and surface diffusion in multicomponent adsorption and liquid chromatography systems , 1996 .

[22]  G. Guiochon,et al.  On the relationship between band broadening and the particle-size distribution of the packing material in liquid chromatography: theory and practice. , 2011, Journal of chromatography. A.

[23]  Karin Westerberg,et al.  Modeling and robust pooling design of a preparative cation-exchange chromatography step for purification of monoclonal antibody monomer from aggregates. , 2014, Journal of chromatography. A.

[24]  I. Lundström Models of protein adsorption on solid surfaces , 1985 .

[25]  Eric von Lieres,et al.  Direct quantification of intraparticle protein diffusion in chromatographic media. , 2006, The journal of physical chemistry. B.

[26]  B. Nilsson,et al.  Model-Based Comparison of Antibody Dimerization in Continuous and Batch-Wise Downstream Processing , 2015 .

[27]  José P. S. Aniceto,et al.  Simulated Moving Bed Strategies and Designs: From Established Systems to the Latest Developments , 2015 .

[28]  Jörg Thömmes,et al.  Modeling the effects of column packing quality and residence time changes on protein monomer/aggregate separation. , 2009, Journal of chromatography. A.

[29]  Marcus Degerman,et al.  Pooling control in variable preparative chromatography processes , 2010, Bioprocess and biosystems engineering.

[30]  Ajoy Velayudhan,et al.  Modeling of Purification Operations in Biotechnology: Enabling Process Development, Optimization, and Scale‐Up , 2007, Biotechnology progress.

[31]  Jürgen Hubbuch,et al.  High‐throughput micro‐scale cultivations and chromatography modeling: Powerful tools for integrated process development , 2015, Biotechnology and bioengineering.

[32]  Suzanne S. Farid,et al.  Multiobjective evolutionary optimization in antibody purification process design , 2013 .

[33]  I. Langmuir THE CONSTITUTION AND FUNDAMENTAL PROPERTIES OF SOLIDS AND LIQUIDS , 1917 .

[34]  Ekta Mahajan,et al.  Improving affinity chromatography resin efficiency using semi-continuous chromatography. , 2012, Journal of chromatography. A.

[35]  Harvey Arellano-Garcia,et al.  Handling Uncertainty in Model-Based Optimal Experimental Design , 2010 .

[36]  Bill Wilson,et al.  In a relationship , 2013 .

[37]  Jörg Raisch,et al.  Optimization of startup and shutdown operation of simulated moving bed chromatographic processes. , 2011, Journal of chromatography. A.

[38]  Simon Kluters,et al.  Application of linear pH gradients for the modeling of ion exchange chromatography: Separation of monoclonal antibody monomer from aggregates. , 2016, Journal of separation science.

[39]  Giorgio Carta,et al.  pH Transients in hydroxyapatite chromatography columns-experimental evidence and phenomenological modeling. , 2010, Journal of chromatography. A.

[40]  Niklas Jakobsson,et al.  Designing robust preparative purification processes with high performance , 2008 .

[41]  J. Selker,et al.  Numerical estimation of multicomponent adsorption isotherms in preparative chromatography: implications of experimental error. , 2001, Journal of chromatography. A.

[42]  Melissa A. Holstein,et al.  Downstream bioprocessing: recent advances and future promise , 2011 .

[43]  L. Whitehead,et al.  A new model of protein adsorption kinetics derived from simultaneous measurement of mass loading and changes in surface energy. , 2007, Langmuir : the ACS journal of surfaces and colloids.

[44]  Marcel Ottens,et al.  pH-gradient ion-exchange chromatography: an analytical tool for design and optimization of protein separations. , 2007, Journal of chromatography. A.

[45]  Sandro Macchietto,et al.  Curvature-Based Methods for Designing Optimally Informative Experiments in Multiresponse Nonlinear Dynamic Situations , 2005 .

[46]  Joon-Ho Koh,et al.  Pore and Surface Diffusion and Bulk-Phase Mass Transfer in Packed and Fluidized Beds , 1998 .

[47]  Marcus Degerman,et al.  Determining Critical Process Parameters and Process Robustness in Preparative Chromatography - A Model-Based Approach , 2009 .

[48]  C. Landes,et al.  pH-dependence of single-protein adsorption and diffusion at a liquid chromatographic interface. , 2016, Journal of separation science.

[49]  Anita M. Katti,et al.  Fundamentals of Preparative and Nonlinear Chromatography , 1994 .

[50]  Massimo Morbidelli,et al.  Protein adsorption on ion exchange resins and monoclonal antibody charge variant modulation. , 2016, Journal of chromatography. A.

[51]  Achim Kienle,et al.  Dynamics and Control of Coupled Continuous Chromatography and Crystallization Processes for the Production of Pure Enantiomers , 2013 .

[52]  B J Geha Practice management: Integrated system , 2004, British Dental Journal.

[53]  Jacek Namieśnik,et al.  Estimating uncertainty in analytical procedures based on chromatographic techniques. , 2010, Journal of chromatography. A.

[54]  W. R. Jackson,et al.  Development of a temperature-responsive agarose-based ion-exchange chromatographic resin. , 2009, Journal of chromatography. A.

[55]  G. Fieg,et al.  Experiment based model development for the enzymatic hydrolysis in a packed-bed reactor with biphasic reactant flow , 2011 .

[56]  Andreas Griewank,et al.  Algorithm 755: ADOL-C: a package for the automatic differentiation of algorithms written in C/C++ , 1996, TOMS.

[57]  A. Ferrero,et al.  Measurement uncertainty , 2006, IEEE Instrumentation & Measurement Magazine.

[58]  George Leitmann,et al.  Dynamics and Control , 2020, Fundamentals of Robotics.

[59]  Massimo Morbidelli,et al.  A continuous multicolumn countercurrent solvent gradient purification (MCSGP) process , 2007, Biotechnology and bioengineering.

[60]  Eric von Lieres,et al.  Investigation of pore diffusion hindrance of monoclonal antibody in hydrophobic interaction chromatography using confocal laser scanning microscopy. , 2007, Journal of chromatography. A.

[61]  D. Antos,et al.  Effects of thermal heterogeneity in hydrophobic interaction chromatography. , 2009, Journal of chromatography. A.

[62]  Andreas Seidel-Morgenstern,et al.  Analytical solutions and moment analysis of general rate model for linear liquid chromatography , 2014 .

[63]  E. von Lieres,et al.  Model-based analysis and quantitative prediction of membrane chromatography: extreme scale-up from 0.08 ml to 1200 ml. , 2014, Journal of chromatography. A.

[64]  M. Morbidelli,et al.  Model simulation and experimental verification of a cation-exchange IgG capture step in batch and continuous chromatography. , 2011, Journal of chromatography. A.

[65]  Robert Tibshirani,et al.  Bootstrap Methods for Standard Errors, Confidence Intervals, and Other Measures of Statistical Accuracy , 1986 .

[66]  David B. Dunson,et al.  Bayesian Data Analysis , 2010 .

[67]  Irving Langmuir,et al.  The constitution and fundamental properties of solids and liquids. Part II.—Liquids , 1917 .

[68]  Jürgen Hubbuch,et al.  Model-integrated process development demonstrated on the optimization of a robotic cation exchange step , 2012 .

[69]  Steven M. Cramer,et al.  Steric mass‐action ion exchange: Displacement profiles and induced salt gradients , 1992 .

[70]  Steven M Cramer,et al.  Use of MiniColumns for linear isotherm parameter estimation and prediction of benchtop column performance. , 2015, Journal of chromatography. A.

[71]  A. Rathore,et al.  Enablers for QbD implementation: Mechanistic modeling for ion-exchange membrane chromatography , 2016 .

[72]  T. Fornstedt,et al.  Injection profiles in liquid chromatography II: predicting accurate injection-profiles for computer-assisted preparative optimizations. , 2011, Journal of chromatography. A.

[73]  Jürgen Hubbuch,et al.  Analytical characterization of complex, biotechnological feedstocks by pH gradient ion exchange chromatography for purification process development. , 2013, Journal of chromatography. A.

[74]  C. Horváth,et al.  Interplay of hydrophobic and electrostatic interactions in biopolymer chromatography. Effect of salts on the retention of proteins. , 1989, Journal of chromatography.

[75]  Dorota Antos,et al.  Altering efficiency of hydrophobic interaction chromatography by combined salt and temperature effects. , 2009, Journal of chromatography. A.

[76]  Torgny Fornstedt,et al.  Impact of an error in the column hold-up time for correct adsorption isotherm determination in chromatography I. Even a small error can lead to a misunderstanding of the retention mechanism. , 2008, Journal of chromatography. A.

[77]  Nicanor Quijano Control Systems Technology Lab Lecture Notes , 2005 .

[78]  Daniel G Bracewell,et al.  Integration of scale-down experimentation and general rate modelling to predict manufacturing scale chromatographic separations. , 2010, Journal of chromatography. A.

[79]  D. Antos,et al.  Multiple-injection technique for isolating a target protein from multicomponent mixtures. , 2011, Journal of chromatography. A.

[80]  Marcus Degerman,et al.  A Model‐Based Approach to Determine the Design Space of Preparative Chromatography , 2009 .

[81]  Dorota Antos,et al.  Protein separations with induced pH gradients using cation-exchange chromatographic columns containing weak acid groups. , 2008, Journal of chromatography. A.

[82]  Asif Ladiwala,et al.  Competitive adsorption of labeled and native protein in confocal laser scanning microscopy , 2006, Biotechnology and bioengineering.

[83]  Vincent Lefèvre,et al.  MPFR: A multiple-precision binary floating-point library with correct rounding , 2007, TOMS.

[84]  Jürgen Hubbuch,et al.  Downstream process development strategies for effective bioprocesses: Trends, progress, and combinatorial approaches , 2017, Engineering in life sciences.

[85]  Marcus Degerman,et al.  Combined effects of potassium chloride and ethanol as mobile phase modulators on hydrophobic interaction and reversed-phase chromatography of three insulin variants. , 2015, Journal of chromatography. A.

[86]  Richard D. Braatz,et al.  Control systems technology in the advanced manufacturing of biologic drugs , 2015, 2015 IEEE Conference on Control Applications (CCA).

[87]  G. Carta,et al.  Particle‐size distribution effects in batch adsorption , 2003 .

[88]  Marcel Ottens,et al.  Purifying biopharmaceuticals: knowledge-based chromatographic process development. , 2014, Trends in biotechnology.

[89]  Massimo Morbidelli,et al.  A semicontinuous 3-column countercurrent solvent gradient purification (MCSGP) process. , 2008, Biotechnology and bioengineering.

[90]  K. Miyabe,et al.  Moment equations for chromatography based on Langmuir type reaction kinetics. , 2014, Journal of chromatography. A.

[91]  Ales Podgornik,et al.  Characterization of ion exchange stationary phases via pH transition profiles. , 2008, Journal of chromatography. A.

[92]  Cristina Maria Quintella,et al.  Statistical designs and response surface techniques for the optimization of chromatographic systems. , 2007, Journal of chromatography. A.

[93]  Giorgio Carta,et al.  pH transitions in cation exchange chromatographic columns containing weak acid groups. , 2007, Journal of chromatography. A.

[94]  L. Schapira,et al.  Handling uncertainty , 2013, Supportive Care in Cancer.

[96]  Suresh Vunnum,et al.  Cation exchange surface-mediated denaturation of an aglycosylated immunoglobulin (IgG1). , 2012, Journal of chromatography. A.

[97]  Jon C. Helton,et al.  Survey of sampling-based methods for uncertainty and sensitivity analysis , 2006, Reliab. Eng. Syst. Saf..

[98]  Anders Axelsson,et al.  Model-based optimization of a preparative ion-exchange step for antibody purification. , 2004, Journal of chromatography. A.

[99]  M. Polakovič,et al.  Effect of pH on protein adsorption capacity of strong cation exchangers with grafted layer. , 2011, Journal of chromatography. A.

[100]  N. J. Titchener-Hooker,et al.  Ranking bioprocess variables using global sensitivity analysis: a case study in centrifugation , 2007, Bioprocess and biosystems engineering.

[101]  Ana Sofia da Costa Moreira,et al.  Downstream process development of novel valuable bio-therapeutics , 2016 .

[102]  S. Leweke,et al.  Utilizing algorithmic differentiation to efficiently compute chromatograms and parameter sensitivities , 2016 .

[103]  Eric von Lieres,et al.  Fast and accurate parameter sensitivities for the general rate model of column liquid chromatography , 2013, Comput. Chem. Eng..

[104]  Jürgen Hubbuch,et al.  High Throughput Screening for the Design and Optimization of Chromatographic Processes: Assessment of Model Parameter Determination from High Throughput Compatible Data , 2008 .

[105]  Bertrand Iooss,et al.  Global sensitivity analysis based on entropy , 2008 .

[106]  A. Atkinson,et al.  Optimal experimental design in chromatography. , 2008, Journal of chromatography. A.

[107]  Marcel Ottens,et al.  Model-based rational methodology for protein purification process synthesis , 2013 .

[108]  Sunil Chhatre,et al.  Global sensitivity analysis for the determination of parameter importance in the chromatographic purification of polyclonal antibodies , 2008 .

[109]  Qilong Ren,et al.  Simulation of Levulinic Acid Adsorption in Packed Beds using Parallel Pore/Surface Diffusion Model , 2010 .

[110]  Christian Langel,et al.  Extra-column dead volume in simulated moving bed separations: theory and experiments. , 2009, Journal of chromatography. A.

[111]  M. Hafner,et al.  Modeling of dual gradient elution in ion exchange and mixed-mode chromatography. , 2015, Journal of chromatography. A.

[112]  Beckley K. Nfor,et al.  High-throughput isotherm determination and thermodynamic modeling of protein adsorption on mixed mode adsorbents. , 2010, Journal of chromatography. A.

[113]  Eric von Lieres,et al.  A fast and accurate solver for the general rate model of column liquid chromatography , 2010, Comput. Chem. Eng..

[114]  M. Saska,et al.  Continuous Separation of Sugarcane Molasses with a Simulated Moving-Bed Adsorber. Adsorption Equilibria, Kinetics, and Application , 1992 .

[115]  M. Etzel,et al.  Evaluation of Three Kinetic Equations in Models of Protein Purification Using Ion-Exchange Membranes , 2003 .

[116]  A. Rathore,et al.  Mechanistic modeling of ion-exchange process chromatography of charge variants of monoclonal antibody products. , 2015, Journal of chromatography. A.

[117]  G. Guiochon,et al.  Influence of the errors made in the measurement of the extra-column volume on the accuracies of estimates of the column efficiency and the mass transfer kinetics parameters. , 2006, Journal of chromatography. A.

[118]  Bernt Nilsson,et al.  Effects of uncertainties in experimental conditions on the estimation of adsorption model parameters in preparative chromatography , 2013, Comput. Chem. Eng..

[119]  J. Smith,et al.  Chromatographic study of surface diffusion , 1968 .

[120]  K. Miyabe Surface diffusion in reversed-phase liquid chromatography using silica gel stationary phases of different C1 and C18 ligand densities. , 2007, Journal of chromatography. A.

[121]  Jürgen Hubbuch,et al.  Rational and systematic protein purification process development: the next generation. , 2009, Trends in biotechnology.

[122]  Marcel Ottens,et al.  Optimization of biopharmaceutical downstream processes supported by mechanistic models and artificial neural networks , 2017, Biotechnology progress.

[123]  Steven M Cramer,et al.  Methods development in multimodal chromatography with mobile phase modifiers using the steric mass action model. , 2013, Journal of chromatography. A.

[124]  Sandro Macchietto,et al.  Model-based design of experiments for parameter precision: State of the art , 2008 .

[125]  Frank Riske,et al.  Periodic counter-current chromatography -- design and operational considerations for integrated and continuous purification of proteins. , 2012, Biotechnology journal.

[126]  Steven M Cramer,et al.  Model‐based process development for the purification of a modified human growth hormone using multimodal chromatography , 2014, Biotechnology progress.

[127]  L. Petzold,et al.  Sensitivity analysis of differential-algebraic equations: A comparison of methods on a special problem ✩ , 2000 .

[128]  Marcus Degerman,et al.  Constrained optimization of a preparative ion-exchange step for antibody purification. , 2004, Journal of chromatography. A.

[129]  Annegret Wagler,et al.  Optimal cut-times finding strategies for collecting a target component from overloaded elution chromatograms , 2013, Comput. Chem. Eng..

[130]  Jørgen Mollerup,et al.  A review of the thermodynamics of protein association to ligands, protein adsorption, and adsorption isotherms , 2008 .

[131]  J. Hubbuch,et al.  Systematic generation of buffer systems for pH gradient ion exchange chromatography and their application. , 2013, Journal of chromatography. A.

[132]  Jürgen Hubbuch,et al.  Estimation of adsorption isotherm and mass transfer parameters in protein chromatography using artificial neural networks. , 2017, Journal of chromatography. A.

[133]  Eva Sorensen,et al.  Modelling of industrial biopharmaceutical multicomponent chromatography , 2014 .

[134]  E. von Lieres,et al.  Optimizing a chromatographic three component separation: a comparison of mechanistic and empiric modeling approaches. , 2012, Journal of chromatography. A.

[135]  G WERNER,et al.  The measurement of uncertainty , 1961, Clinical pharmacology and therapeutics.

[136]  Chi-Wang Shu,et al.  Efficient Implementation of Weighted ENO Schemes , 1995 .

[137]  F. Durbin,et al.  Numerical Inversion of Laplace Transforms: An Efficient Improvement to Dubner and Abate's Method , 1974, Comput. J..

[138]  Massimo Morbidelli,et al.  Improvement of an overloaded, multi-component, solvent gradient bioseparation through multiobjective optimization. , 2007, Journal of chromatography. A.

[139]  Carol S. Woodward,et al.  Enabling New Flexibility in the SUNDIALS Suite of Nonlinear and Differential/Algebraic Equation Solvers , 2020, ACM Trans. Math. Softw..

[140]  Andrew L. Zydney,et al.  Continuous downstream processing for high value biological products: A Review , 2016, Biotechnology and bioengineering.

[141]  C. Haynes,et al.  Improved isoelectric focusing chromatography on strong anion exchange media via a new model that custom designs mobile phases using simple buffers , 2014, Biotechnology and bioengineering.

[142]  Eva Sørensen,et al.  Optimal Economic Design and Operation of Single‐ and Multi‐column Chromatographic Processes , 2005, Biotechnology progress.

[143]  M. Kaspereit,et al.  Continuous production of single enantiomers at high yields by coupling single column chromatography, racemization, and nanofiltration , 2013 .

[144]  O. Thomas,et al.  Integrated system for temperature-controlled fast protein liquid chromatography comprising improved copolymer modified beaded agarose adsorbents and a travelling cooling zone reactor arrangement. , 2013, Journal of chromatography. A.

[145]  C. Haynes,et al.  A new mixed‐mode model for interpreting and predicting protein elution during isoelectric chromatofocusing , 2014, Biotechnology and bioengineering.

[146]  Marcel Ottens,et al.  Size-exclusion chromatographic protein refolding: fundamentals, modeling and operation. , 2010, Journal of chromatography. A.