2025 AIChE Annual Meeting

(261b) Hybrid Modeling Employing Neural Differential Algebraic System of Equations (DAEs) for Enhanced Monoclonal Antibody Production in Fed-Batch Reactors

Authors

Victor Alves - Presenter, West Virginia University
Yu Luo, GSK
Gregory Nierode, GlaxoSmithKline
Sameer Talwar, Duquesne University
Anne Robinson, Carnegie Mellon University
Monoclonal antibodies (mAb) are an essential asset for treating a wide range of diseases including certain types of cancers, autoimmune disorders, and infectious diseases, with an expected market of US$ 300 billion in 2025 [1] and to surpass US$ 700 billion by 2033 [2]. The manufacture of mAbs is mainly performed using Chinese hamster ovary (CHO) cell lines with fed-batch processes, corresponding to around 70% of all mAb production [3]. This process has complex dynamics associated with cell growth, metabolism, and production of mAbs. Nutrients (e.g., glucose, amino acids) are fed regularly in boluses, and undesired byproducts such as ammonia and lactate are also produced, which inhibit cell growth and decrease mAb productivity. Such complex dynamics make it challenging to model this process and require a differential algebraic system of equations (DAE) with nontrivial solution techniques. While modeling of CHO-based mAb production from a mechanistic perspective is present in the literature [4], [5], [6], there is a critical need to embed machine learning and model predictions beyond the capabilities of purely mechanistic approaches.

To achieve this goal, we developed a hybrid model that embeds both mechanistic model equations, which incorporates domain knowledge, with neural network components trained based on experimental data. We aim to generate a model that is both interpretable and able to capture the complexity of mAb production processes. Inspired by recent hybrid modeling approaches such as neural ordinary differential equations (ODEs) and universal differential equations (UDEs) [7], [8] we generated a novel system of so-called neural DAEs [9]. This hybrid formulation is built leveraging the Pyomo.DAE framework [10], which allows for the algebraic modeling and optimization of DAE systems. The training and simulation of the CHO-based mAb production hybrid model is performed using state-of-the art nonlinear programming algorithms, such as Ipopt [11] and dynamic optimization techniques, namely the simultaneous approach [12]. The Pyomo [13] framework can accurately model a non-smooth process that is interrupted by feeding and sampling events based on experimental data available. We train the neural DAE system using data from experiments at small and large scales and evaluate the hybrid model performance across scales.

The hybrid model offers enhanced prediction capabilities while still maintaining domain knowledge. This approach can mitigate extrapolation issues often observed in purely data-driven models and reduce data needed to train the embedded neural networks. In this presentation, we highlight the potential of hybrid modeling enabled by nonlinear programming as a new venue for enhancing biopharmaceutical manufacturing.

References
[1] R.-M. Lu et al., “Development of therapeutic antibodies for the treatment of diseases,” J. Biomed. Sci., vol. 27, no. 1, p. 1, Jan. 2020, doi: 10.1186/s12929-019-0592-z.
[2] “Monoclonal Antibodies Market Size To Worth USD 729.43 Bn By 2033.” Accessed: Apr. 03, 2025. [Online]. Available: https://www.precedenceresearch.com/monoclonal-antibodies-market
[3] K. Liang, H. Luo, and Q. Li, “Enhancing and stabilizing monoclonal antibody production by Chinese hamster ovary (CHO) cells with optimized perfusion culture strategies,” Front. Bioeng. Biotechnol., vol. 11, Jan. 2023, doi: 10.3389/fbioe.2023.1112349.
[4] Y. Luo, V. Kurian, L. Song, E. A. Wells, A. S. Robinson, and B. A. Ogunnaike, “Model-based control of titer and glycosylation in fed-batch mAb production: Modeling and control system development,” AIChE J., vol. 69, no. 4, p. e18075, 2023, doi: 10.1002/aic.18075.
[5] C. Kontoravdi, E. N. Pistikopoulos, and A. Mantalaris, “Systematic development of predictive mathematical models for animal cell cultures,” Comput. Chem. Eng., vol. 34, no. 8, pp. 1192–1198, Aug. 2010, doi: 10.1016/j.compchemeng.2010.03.012.
[6] C. Kontoravdi and I. Jimenez del Val, “Computational tools for predicting and controlling the glycosylation of biopharmaceuticals,” Curr. Opin. Chem. Eng., vol. 22, pp. 89–97, Dec. 2018, doi: 10.1016/j.coche.2018.08.007.
[7] R. T. Q. Chen, Y. Rubanova, J. Bettencourt, and D. Duvenaud, “Neural Ordinary Differential Equations,” Dec. 13, 2019, arXiv: arXiv:1806.07366. doi: 10.48550/arXiv.1806.07366.
[8] C. Rackauckas et al., “Universal Differential Equations for Scientific Machine Learning,” Nov. 02, 2021, arXiv: arXiv:2001.04385. doi: 10.48550/arXiv.2001.04385.
[9] L. R. Lueg, V. Alves, D. Schicksnus, J. R. Kitchin, C. D. Laird, and L. T. Biegler, “A Simultaneous Approach for Training Neural Differential-Algebraic Systems of Equations,” Apr. 07, 2025, arXiv: arXiv:2504.04665. doi: 10.48550/arXiv.2504.04665.
[10] B. Nicholson, J. D. Siirola, J.-P. Watson, V. M. Zavala, and L. T. Biegler, “pyomo.dae: a modeling and automatic discretization framework for optimization with differential and algebraic equations,” Math. Program. Comput., vol. 10, no. 2, pp. 187–223, Jun. 2018, doi: 10.1007/s12532-017-0127-0.
[11] A. Wächter and L. T. Biegler, “On the implementation of an interior-point filter line-search algorithm for large-scale nonlinear programming,” Math. Program., vol. 106, no. 1, pp. 25–57, Mar. 2006, doi: 10.1007/s10107-004-0559-y.
[12] L. T. Biegler, “An overview of simultaneous strategies for dynamic optimization,” Chem. Eng. Process. Process Intensif., vol. 46, no. 11, pp. 1043–1053, Nov. 2007, doi: 10.1016/j.cep.2006.06.021.
[13] M. L. Bynum et al., Pyomo - optimization modeling in Python, Third edition. in Springer Optimization and it’s applications, no. volume 67. Cham, Switzerland: Springer, 2021.