WO2023077683A1 - Procédé d'estimation en ligne d'état de culture cellulaire et de régulation et de commande d'optimisation de réapprovisionnement - Google Patents
Procédé d'estimation en ligne d'état de culture cellulaire et de régulation et de commande d'optimisation de réapprovisionnement Download PDFInfo
- Publication number
- WO2023077683A1 WO2023077683A1 PCT/CN2022/073489 CN2022073489W WO2023077683A1 WO 2023077683 A1 WO2023077683 A1 WO 2023077683A1 CN 2022073489 W CN2022073489 W CN 2022073489W WO 2023077683 A1 WO2023077683 A1 WO 2023077683A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- cell
- state
- cell culture
- estimation
- optimal
- Prior art date
Links
- 238000000034 method Methods 0.000 title claims abstract description 116
- 238000004113 cell culture Methods 0.000 title claims abstract description 93
- 238000005457 optimization Methods 0.000 title claims abstract description 64
- 230000033228 biological regulation Effects 0.000 title claims abstract description 15
- 238000004519 manufacturing process Methods 0.000 claims abstract description 62
- 230000008901 benefit Effects 0.000 claims abstract description 61
- 230000008569 process Effects 0.000 claims abstract description 31
- 230000001580 bacterial effect Effects 0.000 claims abstract description 13
- 238000005096 rolling process Methods 0.000 claims abstract description 10
- 238000005259 measurement Methods 0.000 claims description 37
- 239000000758 substrate Substances 0.000 claims description 32
- 238000000691 measurement method Methods 0.000 claims description 17
- 239000011159 matrix material Substances 0.000 claims description 15
- 230000006870 function Effects 0.000 claims description 14
- 238000005070 sampling Methods 0.000 claims description 10
- 238000001914 filtration Methods 0.000 claims description 7
- 238000006243 chemical reaction Methods 0.000 claims description 6
- 239000002245 particle Substances 0.000 claims description 6
- 230000003595 spectral effect Effects 0.000 claims description 6
- 238000001228 spectrum Methods 0.000 claims description 4
- OILXMJHPFNGGTO-UHFFFAOYSA-N (22E)-(24xi)-24-methylcholesta-5,22-dien-3beta-ol Natural products C1C=C2CC(O)CCC2(C)C2C1C1CCC(C(C)C=CC(C)C(C)C)C1(C)CC2 OILXMJHPFNGGTO-UHFFFAOYSA-N 0.000 claims description 3
- RQOCXCFLRBRBCS-UHFFFAOYSA-N (22E)-cholesta-5,7,22-trien-3beta-ol Natural products C1C(O)CCC2(C)C(CCC3(C(C(C)C=CCC(C)C)CCC33)C)C3=CC=C21 RQOCXCFLRBRBCS-UHFFFAOYSA-N 0.000 claims description 3
- MSWZFWKMSRAUBD-IVMDWMLBSA-N 2-amino-2-deoxy-D-glucopyranose Chemical compound N[C@H]1C(O)O[C@H](CO)[C@@H](O)[C@@H]1O MSWZFWKMSRAUBD-IVMDWMLBSA-N 0.000 claims description 3
- OQMZNAMGEHIHNN-UHFFFAOYSA-N 7-Dehydrostigmasterol Natural products C1C(O)CCC2(C)C(CCC3(C(C(C)C=CC(CC)C(C)C)CCC33)C)C3=CC=C21 OQMZNAMGEHIHNN-UHFFFAOYSA-N 0.000 claims description 3
- DNVPQKQSNYMLRS-NXVQYWJNSA-N Ergosterol Natural products CC(C)[C@@H](C)C=C[C@H](C)[C@H]1CC[C@H]2C3=CC=C4C[C@@H](O)CC[C@]4(C)[C@@H]3CC[C@]12C DNVPQKQSNYMLRS-NXVQYWJNSA-N 0.000 claims description 3
- 101000802640 Homo sapiens Lactosylceramide 4-alpha-galactosyltransferase Proteins 0.000 claims description 3
- 102100035838 Lactosylceramide 4-alpha-galactosyltransferase Human genes 0.000 claims description 3
- 238000010923 batch production Methods 0.000 claims description 3
- MSWZFWKMSRAUBD-UHFFFAOYSA-N beta-D-galactosamine Natural products NC1C(O)OC(CO)C(O)C1O MSWZFWKMSRAUBD-UHFFFAOYSA-N 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 claims description 3
- DNVPQKQSNYMLRS-SOWFXMKYSA-N ergosterol Chemical compound C1[C@@H](O)CC[C@]2(C)[C@H](CC[C@]3([C@H]([C@H](C)/C=C/[C@@H](C)C(C)C)CC[C@H]33)C)C3=CC=C21 DNVPQKQSNYMLRS-SOWFXMKYSA-N 0.000 claims description 3
- 238000002474 experimental method Methods 0.000 claims description 3
- 230000002068 genetic effect Effects 0.000 claims description 3
- 229960002442 glucosamine Drugs 0.000 claims description 3
- 239000001963 growth medium Substances 0.000 claims description 3
- 238000013507 mapping Methods 0.000 claims description 3
- 230000007246 mechanism Effects 0.000 claims description 3
- 238000001216 nucleic acid method Methods 0.000 claims description 3
- 230000004044 response Effects 0.000 claims description 3
- 230000007704 transition Effects 0.000 claims description 3
- 238000009826 distribution Methods 0.000 claims description 2
- 239000002609 medium Substances 0.000 claims description 2
- 230000002194 synthesizing effect Effects 0.000 claims 1
- 230000010261 cell growth Effects 0.000 abstract description 11
- 238000012258 culturing Methods 0.000 abstract 1
- 239000000047 product Substances 0.000 description 28
- 230000012010 growth Effects 0.000 description 12
- 238000013406 biomanufacturing process Methods 0.000 description 11
- 230000015572 biosynthetic process Effects 0.000 description 9
- 230000006872 improvement Effects 0.000 description 9
- RFSUNEUAIZKAJO-ARQDHWQXSA-N Fructose Chemical compound OC[C@H]1O[C@](O)(CO)[C@@H](O)[C@@H]1O RFSUNEUAIZKAJO-ARQDHWQXSA-N 0.000 description 8
- 229930091371 Fructose Natural products 0.000 description 8
- 239000005715 Fructose Substances 0.000 description 8
- XSQUKJJJFZCRTK-UHFFFAOYSA-N Urea Chemical compound NC(N)=O XSQUKJJJFZCRTK-UHFFFAOYSA-N 0.000 description 8
- 239000004202 carbamide Substances 0.000 description 8
- 238000003786 synthesis reaction Methods 0.000 description 8
- 239000002253 acid Substances 0.000 description 7
- 238000013178 mathematical model Methods 0.000 description 7
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 description 6
- 241000894006 Bacteria Species 0.000 description 5
- 230000008859 change Effects 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 238000000855 fermentation Methods 0.000 description 4
- 230000004151 fermentation Effects 0.000 description 4
- 235000015097 nutrients Nutrition 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 3
- 229910052757 nitrogen Inorganic materials 0.000 description 3
- 238000004497 NIR spectroscopy Methods 0.000 description 2
- 238000001069 Raman spectroscopy Methods 0.000 description 2
- 230000001133 acceleration Effects 0.000 description 2
- 230000019522 cellular metabolic process Effects 0.000 description 2
- 238000012937 correction Methods 0.000 description 2
- 230000005764 inhibitory process Effects 0.000 description 2
- 230000002503 metabolic effect Effects 0.000 description 2
- 230000004060 metabolic process Effects 0.000 description 2
- 238000010238 partial least squares regression Methods 0.000 description 2
- 238000012628 principal component regression Methods 0.000 description 2
- 238000004886 process control Methods 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000006227 byproduct Substances 0.000 description 1
- 229910052799 carbon Inorganic materials 0.000 description 1
- -1 cell concentration Substances 0.000 description 1
- 230000009134 cell regulation Effects 0.000 description 1
- 238000004140 cleaning Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000013135 deep learning Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 238000004993 emission spectroscopy Methods 0.000 description 1
- 238000005265 energy consumption Methods 0.000 description 1
- 230000008713 feedback mechanism Effects 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 235000003642 hunger Nutrition 0.000 description 1
- 238000001427 incoherent neutron scattering Methods 0.000 description 1
- 238000011081 inoculation Methods 0.000 description 1
- 230000003834 intracellular effect Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 239000002207 metabolite Substances 0.000 description 1
- 244000005700 microbiome Species 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000037351 starvation Effects 0.000 description 1
- 230000001954 sterilising effect Effects 0.000 description 1
- 238000004659 sterilization and disinfection Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012706 support-vector machine Methods 0.000 description 1
- 238000000844 transformation Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/25—Design optimisation, verification or simulation using particle-based methods
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
- G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
- G01N21/35—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light
- G01N21/359—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light using near infrared light
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/62—Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light
- G01N21/63—Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light optically excited
- G01N21/65—Raman scattering
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/27—Design optimisation, verification or simulation using machine learning, e.g. artificial intelligence, neural networks, support vector machines [SVM] or training a model
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/12—Computing arrangements based on biological models using genetic models
- G06N3/126—Evolutionary algorithms, e.g. genetic algorithms or genetic programming
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
- G16B5/00—ICT specially adapted for modelling or simulations in systems biology, e.g. gene-regulatory networks, protein interaction networks or metabolic networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/04—Constraint-based CAD
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/06—Multi-objective optimisation, e.g. Pareto optimisation using simulated annealing [SA], ant colony algorithms or genetic algorithms [GA]
Definitions
- the invention relates to the technical field of biomanufacturing control, in particular to a method for online estimation of cell culture status and optimization of feeding regulation.
- the formation of products is closely related to the growth and metabolism of bacterial cells, and the supply of nutrient substrates directly affects cell growth and metabolism.
- the feeding of substrates can avoid substrate inhibition, cell starvation and metabolism. Repression, which optimizes the feed rate of substrates according to the different conditions of bacterial cell growth and product synthesis, provides an optimal growth environment for microorganisms, and is an effective means of regulation to improve the efficiency of biomanufacturing processes.
- the feeding method of the cell culture process initially relied on manual experience, and gradually studied according to the mathematical model of cell culture and production goals to solve the optimal feeding trajectory in advance, and then use the feeding trajectory as the process control curve in the actual production process.
- actual cell growth involves complex physical, chemical and biological reactions, and its metabolic flow and kinetic evolution are complex.
- the optimal process curve obtained based on mathematical models may not be able to achieve optimal results. Even seriously affect cell metabolism and product synthesis.
- the optimization control method of cell culture and product production process based on mathematical model is essentially an open-loop control, which fails to feed back the optimization effect and adjust the future feeding strategy.
- the purpose of the present invention is to provide a method for online estimation of cell culture status and optimization of feeding regulation, which integrates estimation of cell growth status, optimization of production benefits, online process control and real-time feedback of results, and can sense the growth status of bacterial cells in a timely manner.
- Rolling optimized feeding at different stages provides an optimal culture environment and maximizes economic benefits.
- the present invention provides a method for online estimation of cell culture status and optimization of feed regulation, comprising the following steps:
- establishing a cell culture state model specifically includes the following steps:
- x(k), u(k) are the cell state and feeding rate at time k, respectively, and f[x(k), u(k)] is the linear or nonlinear relationship between x(k), u(k) function, considering the disturbance noise of the cultivation process as w(k+1).
- step S2 specifically includes the following steps:
- x(T f ), u(T f ) are the cell state and feeding rate at the final moment, respectively, and m[x(T f ),u(T f )] is about x(T f ), u(T f ) is a linear or nonlinear function corresponding to the volume of the culture solution, n[x(k), u(k)] is a linear or nonlinear function about x(k), x(k) corresponding to the substrate concentration, u max and u min denote the upper and lower limits of the cell feeding rate u(k) at time k, respectively.
- the economic benefit optimization target includes:
- the optimization objective J is defined as the maximum output P(T f ) at the end point.
- the cell culture state x(k) in is related to the feeding rate u(k), that is, the benefit optimization objective is expressed as follows:
- L[x(k),u(k)] is the functional relationship between the terminal output, cell culture state and feeding rate
- the cost factor is introduced, and the economic benefit optimization goal selects the ratio of the net profit of a batch production to the production time, that is, the process benefit.
- the specific benefit optimization goal is expressed as follows:
- r is the sales price of the unit product
- c is the cost price of the unit feed
- T p is the time interval between two adjacent production batches of the same bioreactor
- the conversion rate from substrate to product is pursued, and the economic benefit optimization goal selects the ratio of the product amount of a batch to the total amount of feed, that is, the product yield.
- the specific benefit optimization goal is expressed as follows :
- an indirect measurement method is used to estimate the cell state, or a direct measurement method is used to obtain an estimated value of the cell state at the current moment.
- the indirect measurement method when the indirect measurement method is selected to estimate the cell state, it includes the following steps: based on biochemical mechanism or experiment, analyze and synthesize the relationship between the basic variable y(k) in the bacterial cell and the cell state x(k), Construct the measurement equation:
- g[x(k)] is the constructed measurement function, assuming that v(k) is the measurement noise, first measure y(k) by glucosamine method, ergosterol method or nucleic acid method, and then obtain indirectly by estimation method Estimated cell state at the current moment
- the estimation method selects Kalman filter method, extended Karl One of the methods of Mann filter, rolling time domain estimation, unscented Kalman filter, Bayesian estimation, particle filter or finite impulse response filter is used to estimate the cell state. Predict the cell state, update and correct the prediction to obtain the estimated value of the cell state at the current moment according to the current measured value of the basic variable y(k)
- the extended Kalman filter algorithm is used to estimate the cell culture state, including the following steps:
- P(k) is the estimated value of covariance at time k
- Q is the covariance matrix of process noise
- F(k) is the state transition matrix, if the cell culture model is nonlinear, then
- solving the optimal feed rate trajectory satisfying the constraint conditions under the optimized economic benefits specifically includes the following steps:
- the nonlinear programming optimization algorithm is used to solve the economic benefit J optimization problem satisfying the constraints, that is, to solve Get the optimal feeding trajectory U k , where the optimization algorithm includes interior point method, exterior point method, sequential quadratic programming method, genetic algorithm and particle swarm algorithm;
- the present invention takes the economic benefit of the biomanufacturing process as the optimization target, and is different from the commonly used optimal feeding method, and directly integrates production optimization, online control and feedback mechanism into one framework for implementation, at each sampling time , estimate the state of cell culture through the current production basic variables, use the mathematical model to predict the future state and economic benefits, use the feeding rate in the future production time domain as the decision variable to optimize the benefits, and implement the optimization results in time.
- the present invention is based on the measurable state of cell culture. Since most states of the cell growth process (such as intracellular metabolites, cell concentration, substrate concentration, etc.) The measurement of variables, and then use the state estimation method to estimate the cell state that is difficult to obtain online. On the other hand, advanced measuring instruments such as Raman spectroscopy, near-infrared spectroscopy, and emission spectroscopy can be used to directly measure the state of some cell cultures. , these estimation and measurement tools provide the support for rolling optimized feeding.
- Fig. 1 is a schematic diagram of the implementation process of the present invention
- Fig. 2 is to utilize fructose to produce polyhydroxybutyric acid process to optimize feeding rate locus figure;
- Fig. 3 is to utilize urea to produce polyhydroxybutyric acid process to optimize feeding rate locus figure
- Fig. 4 is the graph of the change of bacterial cell specific growth rate
- Fig. 5 is the state change diagram of substrate concentration in the production process of polyhydroxybutyric acid
- Fig. 6 is a diagram showing the state change of fructose concentration during the production of polyhydroxybutyric acid
- Fig. 7 is a diagram of the state change of urea concentration in the production process of polyhydroxybutyric acid
- Fig. 8 is a graph showing the state change of fermentation broth volume during the production of polyhydroxybutyric acid.
- the present invention provides a method for online estimation of cell culture status and optimization of feed regulation, comprising the following steps:
- the present invention establishes the cell culture state model (1), analyzes the production process requirements, and determines the constraint conditions (2)-(4) that the cell culture manufacturing process needs to meet; selects the economic benefit optimization target , if the feeding cost is low, use the maximum yield (5) as the economic benefit index; if the pursuit of high yield and take into account the substrate cost, use the maximum process benefit (6) as the index; if the conversion rate from substrate to product is emphasized, then Select the maximum product yield (7) as the benefit index. Benefit optimization depends on the state of cell culture.
- the first measurement method is to construct the measurement equation (8), by measuring the basic variable y(k) related to the cell state, and indirectly obtain the estimated value of the cell state by the estimation method
- the second measurement method is to use advanced spectral technology to directly measure the estimated value of the state
- the extended Kalman filter and other estimation methods are used to predict the current state from the previous cell state according to (9) (10), and then according to (11)-(13) from the current measurement of the basic variables value, update the revised predicted value to obtain the estimated value Estimated by cell state and the proposed feeding rate U k as input, use the computer numerical solution method to solve the cell state model, and obtain the cell state vector in the future; then use the nonlinear programming optimization method to solve A new optimal feeding trajectory U k is obtained.
- Implement u(k) at the current moment in the feeding trajectory U k in the cell culture production process, and repeat the above process until the production process ends.
- the whole process involves the feed regulation technology of the nutrient substrate in the biomanufacturing process.
- Optimizing the feed rate of the substrate according to the different conditions of bacterial cell growth and product synthesis is an effective regulation method to improve the efficiency of the biomanufacturing process.
- the present invention takes the economic benefits of the biomanufacturing process as the optimization goal, integrates the cell culture process requirements, growth state estimation, online rolling optimization, real-time implementation and other links in a framework, and at each sampling time , estimate the state of cell culture through the current production variables, use the mathematical model to predict the future state and economic benefits, use the feeding rate in the future production time domain as the decision variable to optimize the benefits, and implement the optimization results in time, so rolling operation, this method is suitable for Cell culture models and processes require known biomanufacturing processes.
- Step 1 Establish a mathematical model of the cell culture state:
- x(k) and u(k) are the cell state and feeding rate at time k, respectively, f[ ] is a linear or nonlinear function, and w(k+1) is considered to be the disturbance noise in the culture process.
- the Euler method is used to discretize it with T f sampling intervals.
- Step 2 Determine the constraints of the cultivation process:
- m[ ] and n[ ] are linear or nonlinear functions, u max and u min represent the upper limit and lower limit of the feeding rate respectively;
- formula (2) represents the terminal constraints of the culture process, such as the volume of culture medium in production Do not overflow the bioreactor at the end;
- formula (3) represents the restrictive conditions that need to be met throughout the manufacturing process, such as in order to avoid the metabolic inhibition of bacterial cells caused by excessive substrate concentration, the substrate concentration in the production process must be is lower than a certain upper limit;
- formula (4) is the feed rate constraint imposed by considering the operating capacity of the actual equipment and the impact on the growth of bacterial cells;
- Step 3 Select the economic benefit optimization target:
- the economic benefit objectives of the biomanufacturing process are diverse and should be set according to the actual needs of production. If the substrate cost is low, the economic benefit optimization objective is to maximize the amount of optional products, and the optimization objective J is defined as the maximum output P(T f ) at the end point, and the end point output is generally the product of the product concentration at the end time and the volume of the fermentation broth, and is related to The cell culture state x(k) in the production process is related to the feeding rate u(k), that is, the benefit index is expressed as follows:
- L[ ⁇ ] is to describe the functional relationship between the end point production, cell culture state and feed rate.
- the cost factor is introduced, and the economic benefit optimization goal selects the ratio of the net profit of a batch production to the production time, that is, the process benefit.
- the specific benefit indicators are expressed as follows:
- r is the sales price of the unit product
- c is the cost price of the unit feed
- Tp is the time interval between two adjacent production batches of the same bioreactor, that is, the time required for tank release, cleaning, sterilization, inoculation, etc. operating time.
- the conversion rate of substrate to product is pursued, and the economic benefit optimization goal is to choose the ratio of the product amount of a batch to the total amount of fed materials, that is, the product yield.
- the specific benefit indicators are expressed as follows:
- Cost-effective optimization goals also include best cell growth, lowest by-products, and lowest energy consumption.
- Step 4 Select the cell culture state measurement method:
- the end-point yield P(T f ), process efficiency, conversion rate and other indicators of manufacturing directly depend on the state of cell culture, and the detection of cell culture state is difficult.
- g[ ⁇ ] is the constructed measurement function, assuming that v(k) is the measurement noise.
- Step 1 Spectral preprocessing, the methods include smoothing, wavelet transform, multivariate scattering correction (MSC), standard normal variable transformation (SNV), orthogonal signal correction (OSC) , derivative algorithm (Der), etc.; step 2: selection of characteristic bands, the methods include continuous projection method (SPA), partial least squares method (PLS), uninformative variable elimination method (UVE), etc.; step 3: spectral data and
- the establishment of the cell culture state mapping model includes principal component regression (PCR), partial least squares regression (PLSR), support vector machine regression (SVMR), deep learning methods, etc.
- Step 5 Estimation of cell state based on state model and measurement equation:
- the cell state can be obtained indirectly by using the measurable basic variables and selecting estimation or filtering methods.
- estimation or filtering methods According to different situations such as the form of the state model, the process and the statistical distribution of the measurement noise, Kalman filter, extended Kalman filter, rolling time domain estimation, unscented Kalman filter, Bayesian estimation, particle filter, finite impulse response can be used State estimation methods such as filtering and corresponding extended forms.
- f[ ] and g[ ] are nonlinear in the biological manufacturing process.
- the extended Kalman filter algorithm is used to estimate the cell culture state as follows:
- Step 1 Prediction of cell state and its covariance:
- P(k) is the estimated value of covariance at time k
- Q is the covariance matrix of process noise
- F(k) is the state transition matrix, if the cell culture model f[ ⁇ ] is nonlinear, then
- Step 2 Update of cell state and its covariance:
- Step 6 Online solution and rolling implementation of optimal feeding rate trajectory:
- the benefit index J is related to the future state and the feeding rate trajectory, based on the estimated value of the state at the current moment and the required rate trajectory U k , use the production process state model to iteratively calculate the cell state in the future time domain, and then calculate the benefit index.
- Step 3 According to the future X k+1 , use nonlinear programming optimization method to solve Get the new feeding trajectory U k , the optimization algorithm includes interior point method, exterior point method, sequential quadratic programming method, genetic algorithm, particle swarm algorithm, etc.;
- Step 4 Implement u(k) at the current moment in the feeding trajectory U k in the cell culture production process
- Step 6 If you use the fourth step scheme 2 to directly measure and obtain the cell state value Then go to step 9 of the sixth step; otherwise, apply the solution 1 of the fourth step and go to the following step 7;
- Step 7 From step 1 of the fifth step, predict the state of the current moment from the state of the cell at the previous moment;
- Step 8 Use step 2 of the fifth step to update and correct the predicted value from the newly collected measured value y(k) to obtain the estimated value of the cell state
- Step 9 If the production process is not over, go to Step 2 of Step 6.
- the online estimation and optimization control method of cell culture state in biomanufacturing proposed by the present invention is applied to the process of producing polyhydroxybutyric acid (PHB) by fed-batch fermentation with fructose and urea, and its continuous kinetic state model is:
- x 1 is the concentration of non-PHB substances in the cell
- x 2 is the concentration of product PHB
- x 3 is the concentration of fructose
- x 4 is the concentration of urea
- x 5 is the volume of fermentation broth
- u 1 is the feed rate of fructose
- u 2 is the feed rate of urea.
- the input constraints are:
- the covariance matrix is set to
- the method of the present invention is used to realize the optimization and regulation of cell culture in the PHB manufacturing process, and its optimal feeding trajectory, specific growth rate and state changes are shown in Figure 2- Figure 8 respectively. It can be seen from Figures 2-8 that the entire manufacturing process can be roughly divided into two stages: cell growth and product synthesis.
- the first stage is the cell growth stage (0-25h).
- fructose and urea are supplemented at the same time as a carbon source and a nitrogen source, respectively, to maintain the growth of the cell at the maximum specific growth rate, wherein between 0-10h, Since the initial substrate concentration in the culture medium was sufficient for the growth of the bacteria, the feeding amount of fructose and urea was very small, and then the feeding was gradually increased with the consumption of nutrients.
- the second stage is the product synthesis stage (25-49h). Due to the reduction of the nitrogen source utilization capacity of the bacteria, the number of bacteria could not be obtained in the early stage.
- the online rolling optimization feeding method proposed by the present invention can timely regulate the flow acceleration rate of the nutrient substrate according to the actual situation of bacterial cell culture, thereby effectively improving the yield of biomanufacturing and increasing economic benefits.
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Biophysics (AREA)
- Evolutionary Computation (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Evolutionary Biology (AREA)
- Biochemistry (AREA)
- Chemical & Material Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- General Engineering & Computer Science (AREA)
- Analytical Chemistry (AREA)
- Software Systems (AREA)
- Molecular Biology (AREA)
- Physiology (AREA)
- Medical Informatics (AREA)
- Artificial Intelligence (AREA)
- Computer Hardware Design (AREA)
- Geometry (AREA)
- Biomedical Technology (AREA)
- Data Mining & Analysis (AREA)
- Computing Systems (AREA)
- Mathematical Physics (AREA)
- Computational Linguistics (AREA)
- Biotechnology (AREA)
- Genetics & Genomics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Apparatus Associated With Microorganisms And Enzymes (AREA)
- Micro-Organisms Or Cultivation Processes Thereof (AREA)
Abstract
Un procédé d'estimation en ligne d'état de culture cellulaire et de régulation et de commande d'optimisation de réapprovisionnement est divulgué. Le procédé comprend les étapes suivantes consistant : à établir un modèle d'état de culture cellulaire ; à déterminer une condition de contrainte qui doit être respectée pendant la culture cellulaire ; à sélectionner un objectif d'optimisation de bénéfices économiques ; à effectuer une estimation d'état cellulaire, et à obtenir une valeur d'estimation d'état cellulaire à un moment présent ; à résoudre le modèle d'état de culture cellulaire, et à obtenir un vecteur d'état cellulaire à un moment ultérieur ; à l'aide d'une programmation non linéaire, à obtenir une trajectoire de taux de réapprovisionnement optimale respectant la condition de contrainte ; et pendant la production de culture cellulaire, à mettre en œuvre le taux de réapprovisionnement au moment présent dans la trajectoire de taux de réapprovisionnement, et à optimiser et à résoudre de manière répétée la trajectoire de taux de réapprovisionnement jusqu'à ce que la production soit terminée. Dans la présente invention, une estimation d'état de croissance cellulaire, une optimisation de bénéfices de production, une commande en ligne de processus, et une rétroaction de résultat en temps réel sont intégrées dans un ensemble. Les conditions de croissance de cellules bactériennes peuvent être détectées rapidement, une optimisation de réapprovisionnement par laminage est effectuée pour différents stades, et l'environnement de culture optimal est fourni, maximisant ainsi les bénéfices économiques.
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111301428.6 | 2021-11-04 | ||
CN202111301428.6A CN114036810A (zh) | 2021-11-04 | 2021-11-04 | 一种细胞培养状态在线估计及优化补料调控方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2023077683A1 true WO2023077683A1 (fr) | 2023-05-11 |
Family
ID=80142883
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/CN2022/073489 WO2023077683A1 (fr) | 2021-11-04 | 2022-01-24 | Procédé d'estimation en ligne d'état de culture cellulaire et de régulation et de commande d'optimisation de réapprovisionnement |
Country Status (2)
Country | Link |
---|---|
CN (1) | CN114036810A (fr) |
WO (1) | WO2023077683A1 (fr) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117009831B (zh) * | 2023-10-07 | 2023-12-08 | 山东世纪阳光科技有限公司 | 一种精细化工事故风险预测评估方法 |
CN117422170A (zh) * | 2023-10-19 | 2024-01-19 | 江苏一家园健康科技有限公司 | 基于酸碱度控制灭菌工艺调优方法及系统 |
CN117807548A (zh) * | 2024-02-29 | 2024-04-02 | 江苏新希望生态科技有限公司 | 一种豆芽生长培育环境监测方法 |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114036810A (zh) * | 2021-11-04 | 2022-02-11 | 江南大学 | 一种细胞培养状态在线估计及优化补料调控方法 |
CN115090200B (zh) * | 2022-05-27 | 2023-04-07 | 福建省龙氟新材料有限公司 | 用于电子级氢氟酸制备的自动配料系统及其配料方法 |
CN116300479B (zh) * | 2023-05-22 | 2023-08-22 | 山东卫康生物医药科技有限公司 | 基于干细胞培养法的人参皂甙生产装置的控制方法及系统 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101484572A (zh) * | 2006-07-14 | 2009-07-15 | Abb研究有限公司 | 使产物收率最大化的补料-分批发酵设备的在线优化方法 |
CN101859106A (zh) * | 2010-06-23 | 2010-10-13 | 浙江大学 | 一种发酵生产过程控制方法及应用 |
CN106434414A (zh) * | 2016-07-08 | 2017-02-22 | 广州市微生物研究所 | 一种基于在线监控pH变化率的乳酸菌高细胞密度培养方法 |
CN111210867A (zh) * | 2020-02-10 | 2020-05-29 | 江南大学 | 基于代谢分析及酶调控的微生物生长状态估计方法 |
WO2020224779A1 (fr) * | 2019-05-08 | 2020-11-12 | Insilico Biotechnology Ag | Procédé et moyens d'optimisation d'une production biotechnologique |
CN114036810A (zh) * | 2021-11-04 | 2022-02-11 | 江南大学 | 一种细胞培养状态在线估计及优化补料调控方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP7092879B2 (ja) * | 2017-12-29 | 2022-06-28 | エフ.ホフマン-ラ ロシュ アーゲー | 細胞培養物の代謝状態の予測 |
US11542564B2 (en) * | 2020-02-20 | 2023-01-03 | Sartorius Stedim Data Analytics Ab | Computer-implemented method, computer program product and hybrid system for cell metabolism state observer |
CN113326618B (zh) * | 2021-06-02 | 2022-07-15 | 江南大学 | 连续发酵过程中估计培养基初始条件的方法 |
-
2021
- 2021-11-04 CN CN202111301428.6A patent/CN114036810A/zh active Pending
-
2022
- 2022-01-24 WO PCT/CN2022/073489 patent/WO2023077683A1/fr unknown
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101484572A (zh) * | 2006-07-14 | 2009-07-15 | Abb研究有限公司 | 使产物收率最大化的补料-分批发酵设备的在线优化方法 |
CN101859106A (zh) * | 2010-06-23 | 2010-10-13 | 浙江大学 | 一种发酵生产过程控制方法及应用 |
CN106434414A (zh) * | 2016-07-08 | 2017-02-22 | 广州市微生物研究所 | 一种基于在线监控pH变化率的乳酸菌高细胞密度培养方法 |
WO2020224779A1 (fr) * | 2019-05-08 | 2020-11-12 | Insilico Biotechnology Ag | Procédé et moyens d'optimisation d'une production biotechnologique |
CN111210867A (zh) * | 2020-02-10 | 2020-05-29 | 江南大学 | 基于代谢分析及酶调控的微生物生长状态估计方法 |
CN114036810A (zh) * | 2021-11-04 | 2022-02-11 | 江南大学 | 一种细胞培养状态在线估计及优化补料调控方法 |
Non-Patent Citations (2)
Title |
---|
"Master's Thesis", 1 June 2020, JIANGNAN UNIVERSITY, CN, article SHI, BOWEN: "Feeding Control of Fermentation Process Based on Control Vector Parameterization Approach", pages: 1 - 65, XP009545530, DOI: 10.27169/d.cnki.gwqgu.2020.000633 * |
WU, JIE ET AL.: "Optimal Feeding Strategy for Biomanufacturing Based on Economic Predictive Control", PROCEEDINGS OF 32ND CHINESE PROCESS CONTROL CONFERENCE (CPCC2021), 31 August 2021 (2021-08-31) * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117009831B (zh) * | 2023-10-07 | 2023-12-08 | 山东世纪阳光科技有限公司 | 一种精细化工事故风险预测评估方法 |
CN117422170A (zh) * | 2023-10-19 | 2024-01-19 | 江苏一家园健康科技有限公司 | 基于酸碱度控制灭菌工艺调优方法及系统 |
CN117422170B (zh) * | 2023-10-19 | 2024-04-09 | 江苏一家园健康科技有限公司 | 基于酸碱度控制灭菌工艺调优方法及系统 |
CN117807548A (zh) * | 2024-02-29 | 2024-04-02 | 江苏新希望生态科技有限公司 | 一种豆芽生长培育环境监测方法 |
CN117807548B (zh) * | 2024-02-29 | 2024-05-10 | 江苏新希望生态科技有限公司 | 一种豆芽生长培育环境监测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN114036810A (zh) | 2022-02-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
WO2023077683A1 (fr) | Procédé d'estimation en ligne d'état de culture cellulaire et de régulation et de commande d'optimisation de réapprovisionnement | |
Mears et al. | A review of control strategies for manipulating the feed rate in fed-batch fermentation processes | |
EP2147355B1 (fr) | Contrôle prédictif à modèle non linéaire d'un processus de fermentation de biocarburant | |
Komives et al. | Bioreactor state estimation and control | |
CN111210867B (zh) | 基于代谢分析及酶调控的微生物生长状态估计方法 | |
Yüzgeç et al. | On-line evolutionary optimization of an industrial fed-batch yeast fermentation process | |
US20220208297A1 (en) | Predicting cell culture performance in bioreactors | |
Ochoa et al. | Integrating real-time optimization and control for optimal operation: Application to the bio-ethanol process | |
CN103792845B (zh) | 一种金霉素发酵过程补糖速率优化控制的方法和系统 | |
Scaglia et al. | Linear algebra based controller design applied to a bench-scale oenological alcoholic fermentation | |
Yüzgeç | Performance comparison of differential evolution techniques on optimization of feeding profile for an industrial scale baker’s yeast fermentation process | |
Montague et al. | Fermentation monitoring and control: a perspective | |
Luna et al. | Iterative modeling and optimization of biomass production using experimental feedback | |
Luna et al. | Optimal design of dynamic experiments in the development of cybernetic models for bioreactors | |
Julien et al. | Bioreactor monitoring, modeling, and simulation | |
Rogers et al. | Investigating physics-informed neural networks for bioprocess hybrid model construction | |
Kiran et al. | Control of continuous fed-batch fermentation process using neural network based model predictive controller | |
Gunawan et al. | Washout and non-washout solutions of a system describing microbial fermentation process under the influence of growth inhibitions and maximal concentration of yeast cells | |
Yin et al. | Modeling and parameter identification for a nonlinear multi-stage system for dha regulon in batch culture | |
Assar et al. | Reconciling competing models: a case study of wine fermentation kinetics | |
CN116661517B (zh) | 基于物联网的复合微生物肥料发酵温度智能调控系统 | |
CN111308893B (zh) | 生物反应器流加调控的非均匀网格优化方法 | |
CN117813371A (zh) | 用于预测用于监测消化池的状态的指标的方法和设备 | |
CN116913391A (zh) | 一种生物制造过程的代谢通量优化求解方法和系统 | |
Petre et al. | Adaptive control of a fermentation bioprocess for lactic acid production |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 22888724 Country of ref document: EP Kind code of ref document: A1 |