JP7092879B2 - 細胞培養物の代謝状態の予測 - Google Patents
細胞培養物の代謝状態の予測 Download PDFInfo
- Publication number
- JP7092879B2 JP7092879B2 JP2020534849A JP2020534849A JP7092879B2 JP 7092879 B2 JP7092879 B2 JP 7092879B2 JP 2020534849 A JP2020534849 A JP 2020534849A JP 2020534849 A JP2020534849 A JP 2020534849A JP 7092879 B2 JP7092879 B2 JP 7092879B2
- Authority
- JP
- Japan
- Prior art keywords
- extracellular
- flux
- metabolites
- cell
- metabolite
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000004113 cell culture Methods 0.000 title claims description 133
- 235000020938 metabolic status Nutrition 0.000 title claims description 12
- 239000002207 metabolite Substances 0.000 claims description 547
- 230000004907 flux Effects 0.000 claims description 531
- 230000003834 intracellular effect Effects 0.000 claims description 228
- 230000002503 metabolic effect Effects 0.000 claims description 215
- 238000000034 method Methods 0.000 claims description 198
- 238000012549 training Methods 0.000 claims description 133
- 238000005259 measurement Methods 0.000 claims description 114
- 230000008569 process Effects 0.000 claims description 81
- 239000002609 medium Substances 0.000 claims description 77
- 238000006243 chemical reaction Methods 0.000 claims description 64
- 239000000126 substance Substances 0.000 claims description 52
- 238000013528 artificial neural network Methods 0.000 claims description 46
- 238000004458 analytical method Methods 0.000 claims description 40
- 230000006870 function Effects 0.000 claims description 39
- 238000004364 calculation method Methods 0.000 claims description 33
- 239000001963 growth medium Substances 0.000 claims description 30
- 102000003855 L-lactate dehydrogenase Human genes 0.000 claims description 28
- 108700023483 L-lactate dehydrogenases Proteins 0.000 claims description 28
- 238000001739 density measurement Methods 0.000 claims description 24
- 238000003860 storage Methods 0.000 claims description 23
- 150000001413 amino acids Chemical class 0.000 claims description 21
- 239000000203 mixture Substances 0.000 claims description 14
- 238000004519 manufacturing process Methods 0.000 claims description 13
- 230000036961 partial effect Effects 0.000 claims description 11
- 230000004663 cell proliferation Effects 0.000 claims description 10
- 230000001413 cellular effect Effects 0.000 claims description 10
- 230000000670 limiting effect Effects 0.000 claims description 10
- 238000010801 machine learning Methods 0.000 claims description 7
- 238000012544 monitoring process Methods 0.000 claims description 7
- 238000012258 culturing Methods 0.000 claims description 6
- 230000010412 perfusion Effects 0.000 claims description 6
- 230000000144 pharmacologic effect Effects 0.000 claims description 6
- 238000012886 linear function Methods 0.000 claims description 4
- 210000004027 cell Anatomy 0.000 description 342
- 238000000855 fermentation Methods 0.000 description 46
- 230000004151 fermentation Effects 0.000 description 45
- WQZGKKKJIJFFOK-GASJEMHNSA-N Glucose Natural products OC[C@H]1OC(O)[C@H](O)[C@@H](O)[C@@H]1O WQZGKKKJIJFFOK-GASJEMHNSA-N 0.000 description 34
- 239000008103 glucose Substances 0.000 description 33
- 239000000463 material Substances 0.000 description 33
- 238000012360 testing method Methods 0.000 description 33
- 239000011159 matrix material Substances 0.000 description 31
- 239000013598 vector Substances 0.000 description 31
- 230000004060 metabolic process Effects 0.000 description 29
- 239000000047 product Substances 0.000 description 29
- 230000008859 change Effects 0.000 description 27
- 108090000623 proteins and genes Proteins 0.000 description 26
- 230000000875 corresponding effect Effects 0.000 description 25
- JVTAAEKCZFNVCJ-UHFFFAOYSA-N lactic acid Chemical compound CC(O)C(O)=O JVTAAEKCZFNVCJ-UHFFFAOYSA-N 0.000 description 22
- JVTAAEKCZFNVCJ-UHFFFAOYSA-M Lactate Chemical compound CC(O)C([O-])=O JVTAAEKCZFNVCJ-UHFFFAOYSA-M 0.000 description 20
- 238000005070 sampling Methods 0.000 description 20
- 235000001014 amino acid Nutrition 0.000 description 19
- 230000019522 cellular metabolic process Effects 0.000 description 19
- 210000004978 chinese hamster ovary cell Anatomy 0.000 description 17
- 238000009826 distribution Methods 0.000 description 17
- 102000004169 proteins and genes Human genes 0.000 description 17
- 238000010923 batch production Methods 0.000 description 16
- 238000013459 approach Methods 0.000 description 15
- ZDXPYRJPNDTMRX-UHFFFAOYSA-N glutamine Natural products OC(=O)C(N)CCC(N)=O ZDXPYRJPNDTMRX-UHFFFAOYSA-N 0.000 description 15
- 235000004554 glutamine Nutrition 0.000 description 15
- 239000002028 Biomass Substances 0.000 description 12
- 230000015572 biosynthetic process Effects 0.000 description 12
- 230000014509 gene expression Effects 0.000 description 12
- 230000032258 transport Effects 0.000 description 12
- DHMQDGOQFOQNFH-UHFFFAOYSA-N Glycine Chemical compound NCC(O)=O DHMQDGOQFOQNFH-UHFFFAOYSA-N 0.000 description 11
- 238000007792 addition Methods 0.000 description 11
- 238000005206 flow analysis Methods 0.000 description 11
- 239000004310 lactic acid Substances 0.000 description 11
- 235000014655 lactic acid Nutrition 0.000 description 11
- 238000012545 processing Methods 0.000 description 10
- 239000012429 reaction media Substances 0.000 description 10
- 230000000694 effects Effects 0.000 description 9
- 238000005457 optimization Methods 0.000 description 9
- QGZKDVFQNNGYKY-UHFFFAOYSA-N Ammonia Chemical compound N QGZKDVFQNNGYKY-UHFFFAOYSA-N 0.000 description 8
- 235000015097 nutrients Nutrition 0.000 description 8
- 230000004102 tricarboxylic acid cycle Effects 0.000 description 8
- WHUUTDBJXJRKMK-UHFFFAOYSA-N Glutamic acid Natural products OC(=O)C(N)CCC(O)=O WHUUTDBJXJRKMK-UHFFFAOYSA-N 0.000 description 7
- 235000018102 proteins Nutrition 0.000 description 7
- 239000000758 substrate Substances 0.000 description 7
- 230000009897 systematic effect Effects 0.000 description 7
- QNAYBMKLOCPYGJ-REOHCLBHSA-N L-alanine Chemical compound C[C@H](N)C(O)=O QNAYBMKLOCPYGJ-REOHCLBHSA-N 0.000 description 6
- CKLJMWTZIZZHCS-REOHCLBHSA-N L-aspartic acid Chemical compound OC(=O)[C@@H](N)CC(O)=O CKLJMWTZIZZHCS-REOHCLBHSA-N 0.000 description 6
- 235000004279 alanine Nutrition 0.000 description 6
- 230000008901 benefit Effects 0.000 description 6
- 238000009472 formulation Methods 0.000 description 6
- 239000000543 intermediate Substances 0.000 description 6
- MTCFGRXMJLQNBG-UHFFFAOYSA-N serine Chemical compound OCC(N)C(O)=O MTCFGRXMJLQNBG-UHFFFAOYSA-N 0.000 description 6
- 239000004471 Glycine Substances 0.000 description 5
- 235000003704 aspartic acid Nutrition 0.000 description 5
- OQFSQFPPLPISGP-UHFFFAOYSA-N beta-carboxyaspartic acid Natural products OC(=O)C(N)C(C(O)=O)C(O)=O OQFSQFPPLPISGP-UHFFFAOYSA-N 0.000 description 5
- 230000001276 controlling effect Effects 0.000 description 5
- 239000000284 extract Substances 0.000 description 5
- 235000013922 glutamic acid Nutrition 0.000 description 5
- 239000004220 glutamic acid Substances 0.000 description 5
- 230000007246 mechanism Effects 0.000 description 5
- QGZKDVFQNNGYKY-UHFFFAOYSA-O Ammonium Chemical compound [NH4+] QGZKDVFQNNGYKY-UHFFFAOYSA-O 0.000 description 4
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical compound [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 description 4
- 108091028043 Nucleic acid sequence Proteins 0.000 description 4
- LCTONWCANYUPML-UHFFFAOYSA-M Pyruvate Chemical compound CC(=O)C([O-])=O LCTONWCANYUPML-UHFFFAOYSA-M 0.000 description 4
- LCTONWCANYUPML-UHFFFAOYSA-N Pyruvic acid Chemical compound CC(=O)C(O)=O LCTONWCANYUPML-UHFFFAOYSA-N 0.000 description 4
- 230000003698 anagen phase Effects 0.000 description 4
- 230000006399 behavior Effects 0.000 description 4
- 238000012937 correction Methods 0.000 description 4
- 230000001419 dependent effect Effects 0.000 description 4
- 230000016784 immunoglobulin production Effects 0.000 description 4
- 229940076788 pyruvate Drugs 0.000 description 4
- 230000002829 reductive effect Effects 0.000 description 4
- 235000004400 serine Nutrition 0.000 description 4
- 230000003068 static effect Effects 0.000 description 4
- 239000011573 trace mineral Substances 0.000 description 4
- 235000013619 trace mineral Nutrition 0.000 description 4
- 238000010200 validation analysis Methods 0.000 description 4
- 102000004190 Enzymes Human genes 0.000 description 3
- 108090000790 Enzymes Proteins 0.000 description 3
- ZDXPYRJPNDTMRX-VKHMYHEASA-N L-glutamine Chemical compound OC(=O)[C@@H](N)CCC(N)=O ZDXPYRJPNDTMRX-VKHMYHEASA-N 0.000 description 3
- 229910021529 ammonia Inorganic materials 0.000 description 3
- 238000005842 biochemical reaction Methods 0.000 description 3
- 238000004422 calculation algorithm Methods 0.000 description 3
- 239000006143 cell culture medium Substances 0.000 description 3
- 230000002596 correlated effect Effects 0.000 description 3
- 230000007423 decrease Effects 0.000 description 3
- 230000007812 deficiency Effects 0.000 description 3
- 210000003527 eukaryotic cell Anatomy 0.000 description 3
- 238000002474 experimental method Methods 0.000 description 3
- 102000037865 fusion proteins Human genes 0.000 description 3
- 108020001507 fusion proteins Proteins 0.000 description 3
- 230000002414 glycolytic effect Effects 0.000 description 3
- 239000007788 liquid Substances 0.000 description 3
- 238000012269 metabolic engineering Methods 0.000 description 3
- 230000037353 metabolic pathway Effects 0.000 description 3
- 210000002569 neuron Anatomy 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 3
- 238000004886 process control Methods 0.000 description 3
- 230000000717 retained effect Effects 0.000 description 3
- 241000894007 species Species 0.000 description 3
- 235000000346 sugar Nutrition 0.000 description 3
- 150000008163 sugars Chemical class 0.000 description 3
- 238000012706 support-vector machine Methods 0.000 description 3
- 238000001890 transfection Methods 0.000 description 3
- 238000012546 transfer Methods 0.000 description 3
- MTCFGRXMJLQNBG-REOHCLBHSA-N (2S)-2-Amino-3-hydroxypropansäure Chemical compound OC[C@H](N)C(O)=O MTCFGRXMJLQNBG-REOHCLBHSA-N 0.000 description 2
- BJEPYKJPYRNKOW-REOHCLBHSA-N (S)-malic acid Chemical compound OC(=O)[C@@H](O)CC(O)=O BJEPYKJPYRNKOW-REOHCLBHSA-N 0.000 description 2
- 241000588724 Escherichia coli Species 0.000 description 2
- WOZVBTAIEJXDNA-UHFFFAOYSA-N P(=O)(O)(O)O.C(O)(O)=O.C(O)(O)=O.C(O)(O)=O.C(O)(O)=O.C(O)(O)=O Chemical compound P(=O)(O)(O)O.C(O)(O)=O.C(O)(O)=O.C(O)(O)=O.C(O)(O)=O.C(O)(O)=O WOZVBTAIEJXDNA-UHFFFAOYSA-N 0.000 description 2
- 230000004913 activation Effects 0.000 description 2
- 230000006978 adaptation Effects 0.000 description 2
- BJEPYKJPYRNKOW-UHFFFAOYSA-N alpha-hydroxysuccinic acid Natural products OC(=O)C(O)CC(O)=O BJEPYKJPYRNKOW-UHFFFAOYSA-N 0.000 description 2
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 2
- 238000012512 characterization method Methods 0.000 description 2
- 238000000205 computational method Methods 0.000 description 2
- 238000002790 cross-validation Methods 0.000 description 2
- 210000004748 cultured cell Anatomy 0.000 description 2
- 230000034994 death Effects 0.000 description 2
- 238000009795 derivation Methods 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000007613 environmental effect Effects 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 230000002349 favourable effect Effects 0.000 description 2
- 239000007789 gas Substances 0.000 description 2
- 230000007967 glucose restriction Effects 0.000 description 2
- 230000012010 growth Effects 0.000 description 2
- 230000005764 inhibitory process Effects 0.000 description 2
- 239000002054 inoculum Substances 0.000 description 2
- 229910052742 iron Inorganic materials 0.000 description 2
- 230000002427 irreversible effect Effects 0.000 description 2
- 238000002372 labelling Methods 0.000 description 2
- 239000001630 malic acid Substances 0.000 description 2
- 235000011090 malic acid Nutrition 0.000 description 2
- 239000001301 oxygen Substances 0.000 description 2
- 229910052760 oxygen Inorganic materials 0.000 description 2
- 230000037361 pathway Effects 0.000 description 2
- 229940107700 pyruvic acid Drugs 0.000 description 2
- 230000001105 regulatory effect Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000002441 reversible effect Effects 0.000 description 2
- 230000019491 signal transduction Effects 0.000 description 2
- 238000000528 statistical test Methods 0.000 description 2
- 230000002123 temporal effect Effects 0.000 description 2
- 231100000331 toxic Toxicity 0.000 description 2
- 230000002588 toxic effect Effects 0.000 description 2
- 230000002103 transcriptional effect Effects 0.000 description 2
- 230000014616 translation Effects 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 241000699802 Cricetulus griseus Species 0.000 description 1
- 101710088194 Dehydrogenase Proteins 0.000 description 1
- 230000005526 G1 to G0 transition Effects 0.000 description 1
- WHUUTDBJXJRKMK-VKHMYHEASA-N L-glutamic acid Chemical compound OC(=O)[C@@H](N)CCC(O)=O WHUUTDBJXJRKMK-VKHMYHEASA-N 0.000 description 1
- 241001465754 Metazoa Species 0.000 description 1
- 238000003335 Production assurance Methods 0.000 description 1
- 241001522306 Serinus serinus Species 0.000 description 1
- 241000700605 Viruses Species 0.000 description 1
- 239000002253 acid Substances 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- VFRROHXSMXFLSN-KCDKBNATSA-N aldehydo-D-galactose 6-phosphate Chemical compound OP(=O)(O)OC[C@@H](O)[C@H](O)[C@H](O)[C@@H](O)C=O VFRROHXSMXFLSN-KCDKBNATSA-N 0.000 description 1
- 230000003190 augmentative effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- WQZGKKKJIJFFOK-VFUOTHLCSA-N beta-D-glucose Chemical compound OC[C@H]1O[C@@H](O)[C@H](O)[C@@H](O)[C@@H]1O WQZGKKKJIJFFOK-VFUOTHLCSA-N 0.000 description 1
- 238000010364 biochemical engineering Methods 0.000 description 1
- 230000000740 bleeding effect Effects 0.000 description 1
- 239000000872 buffer Substances 0.000 description 1
- 239000006227 byproduct Substances 0.000 description 1
- 230000006652 catabolic pathway Effects 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 230000030833 cell death Effects 0.000 description 1
- 230000003833 cell viability Effects 0.000 description 1
- 238000010367 cloning Methods 0.000 description 1
- 239000005515 coenzyme Substances 0.000 description 1
- 238000011109 contamination Methods 0.000 description 1
- 238000010219 correlation analysis Methods 0.000 description 1
- 230000009089 cytolysis Effects 0.000 description 1
- 210000000172 cytosol Anatomy 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 230000000593 degrading effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 230000003292 diminished effect Effects 0.000 description 1
- 230000029142 excretion Effects 0.000 description 1
- 239000003925 fat Substances 0.000 description 1
- 235000019197 fats Nutrition 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000001595 flow curve Methods 0.000 description 1
- 239000012737 fresh medium Substances 0.000 description 1
- 229930182830 galactose Natural products 0.000 description 1
- 230000002068 genetic effect Effects 0.000 description 1
- 150000002303 glucose derivatives Chemical class 0.000 description 1
- 125000000291 glutamic acid group Chemical group N[C@@H](CCC(O)=O)C(=O)* 0.000 description 1
- 238000011478 gradient descent method Methods 0.000 description 1
- 238000004128 high performance liquid chromatography Methods 0.000 description 1
- 238000000126 in silico method Methods 0.000 description 1
- 238000013101 initial test Methods 0.000 description 1
- 238000011081 inoculation Methods 0.000 description 1
- 238000003780 insertion Methods 0.000 description 1
- 230000037431 insertion Effects 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 230000006525 intracellular process Effects 0.000 description 1
- 230000010189 intracellular transport Effects 0.000 description 1
- 238000009533 lab test Methods 0.000 description 1
- 230000037447 lactate metabolism Effects 0.000 description 1
- 229920002521 macromolecule Polymers 0.000 description 1
- 210000004962 mammalian cell Anatomy 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 239000013028 medium composition Substances 0.000 description 1
- 238000006241 metabolic reaction Methods 0.000 description 1
- 230000002438 mitochondrial effect Effects 0.000 description 1
- 238000010295 mobile communication Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000001537 neural effect Effects 0.000 description 1
- 238000003062 neural network model Methods 0.000 description 1
- 229910000069 nitrogen hydride Inorganic materials 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 235000018343 nutrient deficiency Nutrition 0.000 description 1
- 238000011017 operating method Methods 0.000 description 1
- 210000001672 ovary Anatomy 0.000 description 1
- KHPXUQMNIQBQEV-UHFFFAOYSA-N oxaloacetic acid Chemical compound OC(=O)CC(=O)C(O)=O KHPXUQMNIQBQEV-UHFFFAOYSA-N 0.000 description 1
- 230000001590 oxidative effect Effects 0.000 description 1
- -1 pH regulators Substances 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 230000035790 physiological processes and functions Effects 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000000513 principal component analysis Methods 0.000 description 1
- 102000004196 processed proteins & peptides Human genes 0.000 description 1
- 108090000765 processed proteins & peptides Proteins 0.000 description 1
- 238000001243 protein synthesis Methods 0.000 description 1
- 238000000275 quality assurance Methods 0.000 description 1
- 238000011002 quantification Methods 0.000 description 1
- 230000006798 recombination Effects 0.000 description 1
- 238000005215 recombination Methods 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 230000008844 regulatory mechanism Effects 0.000 description 1
- 238000012827 research and development Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 150000003839 salts Chemical class 0.000 description 1
- 230000001932 seasonal effect Effects 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000007447 staining method Methods 0.000 description 1
- 230000002194 synthesizing effect Effects 0.000 description 1
- 230000008685 targeting Effects 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
- 210000001519 tissue Anatomy 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 230000009261 transgenic effect Effects 0.000 description 1
- 230000001960 triggered effect Effects 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
- 235000013343 vitamin Nutrition 0.000 description 1
- 239000011782 vitamin Substances 0.000 description 1
- 229940088594 vitamin Drugs 0.000 description 1
- 229930003231 vitamin Natural products 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
Classifications
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12M—APPARATUS FOR ENZYMOLOGY OR MICROBIOLOGY; APPARATUS FOR CULTURING MICROORGANISMS FOR PRODUCING BIOMASS, FOR GROWING CELLS OR FOR OBTAINING FERMENTATION OR METABOLIC PRODUCTS, i.e. BIOREACTORS OR FERMENTERS
- C12M41/00—Means for regulation, monitoring, measurement or control, e.g. flow regulation
- C12M41/46—Means for regulation, monitoring, measurement or control, e.g. flow regulation of cellular or enzymatic activity or functionality, e.g. cell viability
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/48—Biological material, e.g. blood, urine; Haemocytometers
- G01N33/50—Chemical analysis of biological material, e.g. blood, urine; Testing involving biospecific ligand binding methods; Immunological testing
- G01N33/5005—Chemical analysis of biological material, e.g. blood, urine; Testing involving biospecific ligand binding methods; Immunological testing involving human or animal cells
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12M—APPARATUS FOR ENZYMOLOGY OR MICROBIOLOGY; APPARATUS FOR CULTURING MICROORGANISMS FOR PRODUCING BIOMASS, FOR GROWING CELLS OR FOR OBTAINING FERMENTATION OR METABOLIC PRODUCTS, i.e. BIOREACTORS OR FERMENTERS
- C12M41/00—Means for regulation, monitoring, measurement or control, e.g. flow regulation
- C12M41/30—Means for regulation, monitoring, measurement or control, e.g. flow regulation of concentration
- C12M41/38—Means for regulation, monitoring, measurement or control, e.g. flow regulation of concentration of metabolites or enzymes in the cells
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12N—MICROORGANISMS OR ENZYMES; COMPOSITIONS THEREOF; PROPAGATING, PRESERVING, OR MAINTAINING MICROORGANISMS; MUTATION OR GENETIC ENGINEERING; CULTURE MEDIA
- C12N5/00—Undifferentiated human, animal or plant cells, e.g. cell lines; Tissues; Cultivation or maintenance thereof; Culture media therefor
- C12N5/0018—Culture media for cell or tissue culture
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N20/00—Machine learning
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
-
- 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
- G16B40/00—ICT specially adapted for biostatistics; ICT specially adapted for bioinformatics-related machine learning or data mining, e.g. knowledge discovery or pattern finding
- G16B40/20—Supervised data analysis
-
- 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
- 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
- G16B5/30—Dynamic-time models
Landscapes
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Chemical & Material Sciences (AREA)
- Physics & Mathematics (AREA)
- Biotechnology (AREA)
- Biomedical Technology (AREA)
- General Health & Medical Sciences (AREA)
- Theoretical Computer Science (AREA)
- Zoology (AREA)
- Wood Science & Technology (AREA)
- Organic Chemistry (AREA)
- General Engineering & Computer Science (AREA)
- Molecular Biology (AREA)
- Medical Informatics (AREA)
- Genetics & Genomics (AREA)
- Biochemistry (AREA)
- Cell Biology (AREA)
- Microbiology (AREA)
- Software Systems (AREA)
- Data Mining & Analysis (AREA)
- Biophysics (AREA)
- Analytical Chemistry (AREA)
- Artificial Intelligence (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Evolutionary Biology (AREA)
- Bioinformatics & Computational Biology (AREA)
- Immunology (AREA)
- Computing Systems (AREA)
- Mathematical Physics (AREA)
- Sustainable Development (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physiology (AREA)
- Hematology (AREA)
- Urology & Nephrology (AREA)
- Computational Linguistics (AREA)
- Medicinal Chemistry (AREA)
Description
これに関連して、上記の欠点が少なくとも部分的に回避される程度まで、細胞の代謝状態を予測するための改善された方法と、それに対応して改善されたシステムが必要とされている。
- その時点で測定された複数の測定値を受け取る工程であって、測定値が、細胞培養物の培養培地における代謝モデルの複数種の細胞外代謝産物の濃度と、細胞培養物における細胞の細胞密度測定値を含む、工程;
- 受け取った測定値を、入力パラメータ値として、トレーニング済み機械学習プログラムロジック(MLP)に入力する工程;
- 受け取った測定値を用いて、MLPによって、将来の時点での細胞外代謝産物の細胞外フラックスを予測する工程であって、該将来の時点が、測定値を受け取った時点より後の時点であり、該細胞外フラックスが、細胞内への細胞外代謝産物の取り込み速度および/または細胞から培地への細胞外代謝産物の放出速度である、工程;
- 細胞外代謝産物の細胞外フラックス予測値と代謝モデルの化学量論式とを用いて、将来の時点での細胞内フラックスを計算するために代謝フラックス解析を行う工程。
- 前記トレーニング時点で測定された複数の測定値を受け取る工程であって、前記測定値が、前記少なくとも1つのトレーニング細胞培養物の培養培地における代謝モデルの複数種の細胞外代謝産物の濃度と、前記少なくとも1つのトレーニング細胞培養物における細胞の細胞密度測定値を含む、工程;
- 現在のトレーニング時点の時間表示を受け取る工程;ならびに
- その時点で受け取った測定値と、それぞれの前の時点で受け取った測定値との関数として細胞外代謝産物の細胞外フラックスを計算する工程であって、細胞外フラックスが、細胞内への細胞外代謝産物の取り込み速度および/または培地への細胞外代謝産物の放出速度である、工程;
- MLPをトレーニングすることであって、
○トレーニング時点のそれぞれにおいて受け取った測定値を入力パラメータ値としてMLPに入力する工程、および、そのトレーニング時点の後続の各時点において該後続の時点について計算した細胞外フラックスを、その入力パラメータ値に関連する出力パラメータ値としてMLPに入力する工程;ならびに
○MLPが、入力パラメータ値に基づいて、それぞれの関連する出力パラメータを予測するように学習するように、MLPによる学習プロセスを実施する工程
を含む、MLPをトレーニングすること;
- トレーニングされた該MLPを揮発性または非揮発性の記憶媒体に記憶すること。
(a)細胞外代謝産物の第1のセットを定める工程であって、第1のセットが全ての候補入力代謝産物を含む、工程;
(b)第1のセットにおける各細胞外代謝産物の第1の関連性スコアを、その代謝産物の濃度測定値の関数として計算する工程であって、第1の関連性スコアが、該単一の細胞外代謝産物の細胞外フラックスに関する、それぞれの細胞外代謝産物の濃度の予測能力を示す、工程;
(c)細胞外代謝産物のうち最も高い第1の関連性スコアを有するものだけを、第1のセットから、まだ空である細胞外代謝産物の第2のセットに移し、この代謝産物を第1のセットから除去する工程;
(d)第1のセットにおける各細胞外代謝産物のさらなる関連性スコアを、その代謝産物濃度の測定値と、第2のセットに含まれる全ての細胞外代謝産物濃度の測定値との関数として計算する工程であって、さらなる関連性スコアが、既に第2のセット含まれている代謝産物を考慮した、該単一の細胞外代謝産物の細胞外フラックスに関する、それぞれの候補入力代謝産物の濃度の予測能力を示す、工程;
(e)第1のセットの細胞外代謝産物のうち最も高いさらなる関連性スコアを有するものだけを第2のセットに移し、この代謝産物を第1のセットから除去する工程であって、この代謝産物を含めることによって、該単一の細胞外代謝産物の細胞外フラックスの予測に関する、第2のセットに含まれる代謝産物の情報冗長性の最大限度を第2のセットが超えない場合にだけ、この代謝産物が移される、工程;
(f)第2のセットが情報冗長性の最大限度を超えることなしでは、これ以上の代謝産物を第1のセットから第2のセットに移すことができなくなるまで、工程d)およびe)を繰り返す工程;ならびに
(g)該単一の細胞外代謝産物の細胞外フラックスを予測するために、第2のセットに移された代謝産物だけを選択的に入力パラメータ値として使用する工程
を含む。
- いくつかのバイオリアクターの中で異なる細胞培養物を培養する工程であって、異なる細胞培養物が、特定の細胞タイプの細胞の異なるクローンを含む、工程;
- 本明細書に記載の細胞の代謝状態を予測するための方法の態様および実施例に従って、細胞クローンのそれぞれについて別々にいくつかの時点で、1種または複数種の細胞内代謝産物の細胞内フラックスを計算する工程;
- 代謝的に最も好ましい、1種または複数種の細胞内代謝産物の細胞内フラックス計算値を有する、細胞クローンのうちの1つを特定する工程
を含む。
- 現在の時点および以前の時点で測定された細胞外代謝産物濃度から、細胞外代謝産物のうちの1つまたは複数の現在の細胞外フラックスを計算する工程;
- さらなる代謝フラックス解析を行って、細胞外代謝産物の現在の細胞外フラックス計算値と代謝モデルの化学量論式を用いて現在の時点での現在の細胞内フラックスを計算する工程;および
- 細胞培養物の細胞の現在の代謝状態を表すものとして現在の細胞内フラックス計算値を使用する工程
を含む。
- 第1のインターフェイスを介して、その時点で測定された複数の測定値を受け取る工程であって、前記測定値が、細胞培養物の培養培地における代謝モデルの複数種の細胞外代謝産物の濃度と、細胞培養物における細胞の細胞密度測定値を含む、工程;
- 受け取った測定値を入力パラメータ値としてMLPに入力する工程;
- 前記受け取った測定値を用いて、前記MLPによって将来の時点での前記細胞外代謝産物の細胞外フラックスを予測する工程であって、前記将来の時点が、前記測定値を受け取った時点より後の時点であり、前記細胞外フラックスが、細胞への前記細胞外代謝産物の取り込み速度および/または細胞から前記培地への前記細胞外代謝産物の放出速度である、工程;
- 細胞外フラックス予測値と代謝モデルの化学量論式とを用いて、将来の時点での細胞内フラックスを計算するために代謝フラックス解析を行う工程
を含む。
- 細胞内フラックス予測値を、それぞれの細胞内代謝産物の許容可能な細胞内フラックスの参照値もしくは参照値範囲と比較する;
- 細胞内代謝産物の少なくとも1つの細胞内フラックス計算値の、そのそれぞれの参照値もしくは参照値範囲からの逸脱が限界値を超えたら、ユーザーインターフェイスを介して警告を発する;および/または
- 第2のインターフェイスを介して制御コマンドをバイオリアクターに送る
ように訓練されており、制御コマンドは、逸脱が小さくなるように、バイオリアクターの状態もしくはその中に含まれる培地を変えるように訓練されている。
[本発明1001]
以下の工程を含む、特定の細胞タイプの細胞の細胞培養物の代謝状態を予測するための方法:
- 該特定の細胞タイプの細胞の代謝モデル(402)を準備する工程(102)であって、該代謝モデルが、複数種の細胞内代謝産物および細胞外代謝産物ならびに複数種の細胞内フラックス(410)および細胞外フラックス(408)を含み、該代謝モデルが、該細胞内代謝産物(406)のうちの1つと該細胞外代謝産物(404)のうちの1つとの間の少なくとも1つの化学量論的関係を指定する化学量論式を含む、工程;
- 該細胞培養物の培養中の複数の時点のそれぞれにおける、
○該時点で測定された複数の測定値を受け取る工程(106)であって、該測定値が、該細胞培養物の培養培地における該代謝モデルの複数種の細胞外代謝産物の濃度と、該細胞培養物における該細胞の細胞密度測定値とを含む、工程;
○受け取った該測定値を、入力パラメータ値として、トレーニング済みの機械学習プログラムロジック(MLP)(218)に入力する工程(108);
○受け取った該測定値を用いて、該MLPによって、将来の時点での細胞外代謝産物の細胞外フラックス(408)を予測する工程(110)であって、該将来の時点が、該測定値を受け取った時点より後の時点であり、該細胞外フラックスが、細胞内への細胞外代謝産物の取り込み速度および/または細胞から培地への細胞外代謝産物の放出速度である、工程;
○該細胞外代謝産物の該細胞外フラックスの予測値と該代謝モデルの該化学量論式とを用いて、該将来の時点での細胞内フラックスを計算するために代謝フラックス解析を行う工程(112)。
[本発明1002]
前記方法が、機械学習によりMLPを作成する工程をさらに含み、該作成する工程が、
- トレーニングデータセットを作成することであって、
○前記特定の細胞タイプの細胞の少なくとも1つのトレーニング細胞培養物の培養中の複数のトレーニング時点のそれぞれにおいて、
■該トレーニング時点で測定された複数の測定値を受け取る工程であって、該測定値が、該少なくとも1つのトレーニング細胞培養物の培養培地における代謝モデルの複数種の細胞外代謝産物の濃度と、該少なくとも1つのトレーニング細胞培養物における細胞の細胞密度測定値を含む、工程、
■現在のトレーニング時点の時間表示を受け取る工程、および
■その時点で受け取った測定値と、それぞれの前の時点で受け取った測定値との関数として、該細胞外代謝産物の細胞外フラックスを計算する工程であって、該細胞外フラックスが、細胞内への該細胞外代謝産物の取り込み速度および/または培地への該細胞外代謝産物の放出速度である、工程
を含む、トレーニングデータセットを作成すること、
- 該MLPをトレーニングすることであって、
○該トレーニング時点のそれぞれにおいて受け取った測定値を、入力パラメータ値として該MLPに入力し、そのトレーニング時点の後続の各時点において該後続の時点について計算された細胞外フラックスを、該入力パラメータ値に関連する出力パラメータ値として該MLPに入力する工程;ならびに
○該MLPが、該入力パラメータ値に基づいて、それぞれの該関連する出力パラメータ値を予測するように学習するように、該MLPによる学習プロセスを実施する工程
を含む、該MLPをトレーニングすること、
- トレーニングされた該MLPを揮発性または非揮発性の記憶媒体に記憶すること
を含む、本発明1001の方法。
[本発明1003]
前記特定の細胞タイプの細胞の複数のトレーニング細胞培養物の培養中の複数のトレーニング時点のそれぞれにおいて、前記測定値および時間表示が受け取られ、かつ前記細胞外代謝産物の細胞外フラックスが計算されるように、前記トレーニングデータセットが作成され、ここで、該細胞培養物が、異なるタイプのバイオリアクターの中で培養され、該バイオリアクターのタイプが、フェドバッチバイオリアクター、バッチバイオリアクター、灌流リアクター、ケモスタット、およびスプリットバッチバイオリアクターを含むバイオリアクタータイプのうちの少なくとも2つの異なるバイオリアクタータイプを含む、本発明1002の方法。
[本発明1004]
前記MLPが、1つのニューラルネットワークであるか、またはいくつかのニューラルネットワークからなるシステムである、前記本発明のいずれかの方法。
[本発明1005]
前記複数種の細胞外代謝産物の濃度測定値がそれぞれ、
- 該代謝産物の体積当たりの含有量、特に、質量濃度もしくは物質量濃度、または
- 該体積当たりの含有量と直線的もしくは少なくともほぼ直線的に相関する値、特に、補正もしくは標準化された、該細胞外代謝産物の代謝産物濃度測定値もしくはフラックス測定値
である、前記本発明のいずれかの方法。
[本発明1006]
前記MLPが、いくつかのサブMLPからなるシステムであり、該システムに含まれる個々のサブMLPがそれぞれ、単一の細胞外代謝産物の細胞外フラックスを予測するようにトレーニングされており、かつ、該単一の細胞外代謝産物の将来の時点での細胞外フラックスを予測するのに選択的に用いられる、前記本発明のいずれかの方法。
[本発明1007]
前記MLPが、前記細胞外代謝産物のうちの1つの細胞外フラックスを予測するために、複数種の細胞外代謝産物の濃度測定値を入力パラメータ値として使用し、ここで、細胞外フラックスを求めようとする該細胞外代謝産物のうちの少なくとも2つについて、入力パラメータ値として使用する該複数種の細胞外代謝産物が異なる、前記本発明のいずれかの方法。
[本発明1008]
- 全ての入力候補代謝産物の濃度を複数の時点にわたって測定する工程であって、該入力候補代謝産物が、前記特定のタイプの細胞培養物を有する参照バイオリアクター内で測定可能な程度に利用可能な全ての細胞外代謝産物を含むか、または前記代謝モデルの全ての細胞外代謝産物を含む、工程;
- 細胞外フラックスを予測しようとする単一の細胞外代謝産物のそれぞれについて、該単一の代謝産物の細胞外フラックスを予測するために入力パラメータ値として使用される複数種の細胞外代謝産物を特定するための選択プロセスを行う工程であって、該選択プロセスが、該単一の代謝産物それぞれに関して、
(a)細胞外代謝産物の第1のセットを定める工程であって、該第1のセットが全ての候補入力代謝産物を含む、工程;
(b)該第1のセットにおける各細胞外代謝産物の第1の関連性スコアを、その代謝産物の濃度測定値の関数として計算する工程であって、該第1の関連性スコアが、該単一の細胞外代謝産物の細胞外フラックスに関する、それぞれの細胞外代謝産物の濃度の予測能力を示す、工程;
(c)該細胞外代謝産物のうち最も高い第1の関連性スコアを有するものだけを、該第1のセットから、まだ空である細胞外代謝産物の第2のセットに移し、この代謝産物を該第1のセットから除去する工程;
(d)該第1のセットにおける各細胞外代謝産物のさらなる関連性スコアを、その代謝産物の濃度測定値と、該第2のセットに含まれる全ての細胞外代謝産物の濃度測定値との関数として計算する工程であって、該さらなる関連性スコアが、既に該第2のセット含まれている代謝産物を考慮した、該単一の細胞外代謝産物の細胞外フラックスに関する、それぞれの候補入力代謝産物の濃度の予測能力を示す、工程;
(e)該第1のセットの細胞外代謝産物のうち最も高いさらなる関連性スコアを有するものだけを該第2のセットに移し、この代謝産物を該第1のセットから除去する工程であって、この代謝産物を含めることによって、該単一の細胞外代謝産物の細胞外フラックスの予測に関する、第2のセットに含まれる代謝産物の情報冗長性の最大限度を第2のセットが超えない場合にだけ、この代謝産物が移される、工程;
(f)該第2のセットが該情報冗長性の最大限度を超えることなしでは、これ以上の代謝産物を該第1のセットから該第2のセットに移すことができなくなるまで、工程(d)および(e)を繰り返す工程;ならびに
(g)該単一の細胞外代謝産物の細胞外フラックスを予測するために、該第2のセットに移された代謝産物だけを選択的に入力パラメータ値として使用する工程
を含む、選択プロセスを行う工程
をさらに含む、本発明1007の方法。
[本発明1009]
- 前記第1のセットの代謝産物と、細胞外フラックスを予測しようとする前記単一の代謝産物との間の部分相互情報量スコア(PMIスコア)として、前記第1の関連性スコアが計算され;
- 既に前記第2のセットに含まれている全ての代謝産物を考慮して、該第1のセットの代謝産物と、細胞外フラックスを予測しようとする該単一の代謝産物との間のPMIスコアとして、第2の関連性スコアが計算される、
本発明1008の方法。
[本発明1010]
前記MLPが、単一の代謝産物のそれぞれの細胞外フラックスを予測するために、複数種の細胞外代謝産物の濃度測定値を入力パラメータ値として使用し、該複数種の細胞外代謝産物が、少なくとも1種のアミノ酸、好ましくは少なくとも2種のアミノ酸を含む、前記本発明のいずれかの方法。
[本発明1011]
前記複数の時点が、10分~48時間、好ましくは1~24時間の時間間隔だけ隔てられている、前記本発明のいずれかの方法。
[本発明1012]
前記細胞タイプが、生体分子を得る目的でバイオリアクター内で維持されかつ/または増殖する遺伝子組換えされた細胞タイプである、前記本発明のいずれかの方法。
[本発明1013]
前記将来の時点のそれぞれで、前記代謝モデルの細胞内フラックスのうちのいくつかまたは全ての計算が行われる、前記本発明のいずれかの方法。
[本発明1014]
- それぞれの参照値または参照値範囲から閾値を超えて逸脱している細胞内フラックス(410)計算値を有する、前記代謝モデルの1種または複数種の細胞内代謝産物を特定する工程;
- 細胞増殖または所望の生体分子の産生の制限要因として作用する細胞内フラックスを自動的に特定する工程
をさらに含む、前記本発明のいずれかの方法。
[本発明1015]
以下の工程を含む、特定の細胞タイプの細胞の細胞培養物を含むバイオリアクターをモニタリングおよび/または制御するための方法:
- 本発明1001~1014のいずれかの方法に従って将来の時点での該細胞の細胞内フラックスを計算する工程;
- 該細胞内フラックスの予測値を、それぞれの1種もしくは複数種の細胞内代謝産物の許容可能な細胞内フラックスの参照値もしくは参照値範囲と比較する工程;
- 少なくとも1つの細胞内代謝産物の細胞内フラックスの計算値の、そのそれぞれの参照値もしくは参照値範囲からの逸脱が限界値を超えた場合に、警告を発する工程;ならびに/または
- 該逸脱を小さくするように該バイオリアクターの状態もしくはその中に含まれる培地を変える工程を自動的に開始するための制御コマンドを該バイオリアクターに送る工程であって、該自動工程が、特に培養培地の量および/もしくは組成を変えることを含む、工程。
[本発明1016]
- 本発明1014に従って、細胞増殖もしくは所望の生体分子の産生の制限要因として作用する、細胞の代謝モデル内の反応を特定する工程;および
- 細胞増殖もしくは該生体分子の産生または該生体分子の質が促進されるように、制限要因として作用する細胞内フラックスを改変する物質を選択的に自動的に追加するか、あるいは、このような追加の要求をユーザーインターフェイスを介して発行する、工程
をさらに含む、本発明1015の方法。
[本発明1017]
以下の工程を含む、特定の細胞タイプの細胞の代謝的に好ましいクローンを特定するための方法:
- 複数のバイオリアクターの中で異なる細胞培養物を培養する工程であって、該異なる細胞培養物が、該特定の細胞タイプの細胞の異なるクローンを含む、工程;
- 各細胞クローンについて別々に複数の時点で、本発明1001~1016のいずれかに従って、1種または複数種の細胞内代謝産物の細胞内フラックスを計算する工程;
- 代謝的に最も好ましい、該1種または複数種の細胞内代謝産物の細胞内フラックス計算値を有する、該細胞クローンのうちの1つを特定する工程。
[本発明1018]
- 1種または複数種の細胞外代謝産物の現在の細胞外フラックスを、現在の時点および以前の時点で測定された該細胞外代謝産物の濃度から計算する工程;
- 該細胞外代謝産物の該現在の細胞外フラックスの計算値と代謝モデルの化学量論式とを用いて、さらなる代謝フラックス解析を行って、現在の時点での現在の細胞内フラックス(410)を計算する工程;ならびに
- 該現在の細胞内フラックスの計算値を、細胞培養物の細胞の現在の代謝状態の特徴付けとして使用する工程
をさらに含む、前記本発明のいずれかの方法。
[本発明1019]
- 前記細胞培養物の培養中の前記時点のそれぞれにおいて、該細胞培養物の培地中の乳酸デヒドロゲナーゼ(LDH)濃度の測定値がさらに受け取られ;
- 前記将来の時点のそれぞれにおける前記細胞外代謝産物の細胞外フラックスの予測が、前記細胞密度測定値ではなく補正された細胞密度を用いて、前記MLPによって行われ、該補正された細胞密度の計算が、該時点のそれぞれにおける、
○該LDH濃度測定値の関数として、該細胞培養物の培地中の溶解細胞の密度を計算する工程であって、該関数が、該培地中のLDH濃度と前記特定の細胞タイプの溶解細胞の数との関係を表す、経験的に求められたヒューリスティックな一次関数である、工程;
○該培地中の該細胞密度測定値と該溶解細胞の密度計算値との和として、該補正された細胞密度を計算する工程
を含む、
前記本発明のいずれかの方法。
[本発明1020]
以下を備える、特定の細胞タイプの細胞の細胞培養物の代謝状態を予測するためのシステム(200):
- 1つまたは複数のプロセッサ(202)、
- 該細胞培養物を含むバイオリアクター(204、206、208)からの測定値を受け取るための第1のインターフェイス(210)、ならびに
- 以下を備える、揮発性または非揮発性の記憶媒体(212):
○該特定の細胞タイプの細胞の代謝モデル(402)であって、複数種の細胞内代謝産物および細胞外代謝産物ならびに複数種の細胞内フラックス(410)および細胞外フラックス(408)を含み、該細胞内代謝産物(406)のうちの1つと該細胞外代謝産物(404)のうちの1つとの間の少なくとも1つの化学量論的関係を指定する化学量論式を含む、代謝モデル(402);
○トレーニング済みの機械学習プログラムロジック(MLP)(218);
○該細胞培養物の培養中の複数の時点のそれぞれにおいて以下の工程を含む方法を行うように訓練された、プログラムロジック(220):
■第1のインターフェイスを介して該時点で測定された複数の測定値を受け取る工程(106)であって、該測定値が、該細胞培養物の培養培地における該代謝モデルの複数種の細胞外代謝産物の濃度と、該細胞培養物における該細胞の細胞密度測定値とを含む、工程;
■受け取った該測定値を入力パラメータ値として該MLPに入力する工程(108);
■受け取った該測定値を用いて、該MLPによって、将来の時点での該細胞外代謝産物の細胞外フラックス(408)を予測する工程(110)であって、該将来の時点が、該測定値を受け取った時点より後の時点であり、該細胞外フラックスが、細胞内への該細胞外代謝産物の取り込み速度および/または細胞から培地への該細胞外代謝産物の放出速度である、工程;
■該細胞外フラックス予測値と該代謝モデルの該化学量論式とを用いて、該将来の時点での細胞内フラックスを計算するために代謝フラックス解析を行う工程(112)。
バイオリアクター内のプロセスは数式で表すことができる。最初に、反応培地中の関連性のある物質濃度または量の時間変化(例えば、基質量、産物量、細胞密度の経過)の写像を考慮しなければならない。式は、物質収支に基づいており、物質の反応を表す項と、リアクター内外への任意の物質フローを含む対流項からなる。これは、一般式(付属書の[51]、セクション4.2を参照されたい):
または式:
で成り立つ。
式中、mは反応培地中の物質の量、tはプロセス時間
または
はそれぞれ内外の流量体積、CzuまたはCabはそれぞれ内外の対応する物質の濃度、ならびにQは時間および体積あたりに変換される物質の量である。細胞外代謝産物の量が考慮される場合、反応項は主として細胞による物質の取り込みまたは放出を含む。ファーメンター内の細胞密度が記述される場合、これは細胞の形成と死を含む。
で記録することができ、行は様々な代謝産物に対応し、列は反応に対応する。細胞外代謝産物は省略されている。具体例を付属書のセクション2に示した。
に簡約する。
と、未知のフラックスのサブベクトル
に分け、それに対応して、行列Aを部分行列
および
に分けると、Av=Amvm+Auvuになり、式(2.3)が
になる。
は通常、成り立たない。なぜなら、通常、行列Auは可逆でないからである。このような場合、解空間は無限になることがあるか、または解が存在しないように矛盾が生じることがある。van der Heijden et al.(1994)によって以下の項が導入され、これは、可解性および一貫性の基準に従って前記の連立方程式またはまたはマテリアルフローを分類し[17,39,40]、本発明の態様に従ってモデルを作成するのにも使用される。
Rang(Au)<k-kmが成り立つ場合、前記の式(2.4)は過小決定されている。この場合、全ての未知のフラックスが明確に計算され得るわけではない。なぜなら、代謝ネットワークは、あまりに少ない制限を含むからである。Rang(Au)=k-kmの場合、前記の式は最大で1つの解を有し、決定されたと呼ばれる。
Rang(Au)<lの場合、前記の式は冗長である。これは、線形従属線がAuに存在することを意味する。vmの決定に測定誤差があるために、または代謝ネットワークモデルに不正確さがあるために、これは、通常、解が存在しない一貫性のない式につながり(
が成り立つ)、後の方の項は拡大係数行列を表す。Rang(Au)=lであれば、前記の式は冗長でなく、従って常に一貫している。
フラックスvuは、式(2.4)を用いて明確に計算されれば計算可能とされ、そうでなければ計算可能でない。一貫性のある式が本明細書において想定される。前記の式が過小決定されたら、計算可能でない少なくとも1つのフラックスがある。
フラックスvmは、その値が前記の式の一貫性に影響を及ぼすのであれば均衡可能だと呼ばれ、そうでなければ均衡可能でない。均衡可能なフローは冗長な式でしか生じない。
を用いて解を式:
で表すことができる。これは全てのAuについて定義される。
Klamtらにより発表された、計算可能なフラックスを特定するための方法は[17]に示されている。
が、Auのコアの基底を形成する列を有する行列である場合、
AuT=0
が成り立つ。
は、対応する列ベクトルRがヌルベクトルであれば、フラックス測定値が式(2.4)の可解性に影響を及ぼさないことを示している。上記の式において、rjはRのj番目の列ベクトルを示し、νm,jはj番目のフラックス測定値を示す。
式が過小決定されている場合、すなわち、無限の数の解を有する場合、ユニークな解に到達するために問題を変更するいくつかの手法がある。1つの選択肢は、可能であれば、結果として生じた制限を、代謝ネットワークモデルを拡大することによって、従って、決められた式を作り出すことによって厳しくすることであろう。さらに、さらなる実験定量化によって既知のフラックスの数を増やすことが試みられるかもしれない。細胞内フラックスの決定は13C標識実験によって成し遂げられ得る(付属書の文献[43]を参照されたい)。
一般式:
が得られる。
virrev≧0
によってさらに制限することができる。式中、virrevは全ての不可逆的フラックスのベクトルである。
冗長な連立方程式の場合、通常、均衡化することができるフローに矛盾があるために、式(2.4)を解くフラックス分布は無い。完全に正確に定義されたネットワークモデルでも、vmの決定における測定誤差によって引き起こされる不一致が通常、生じる。実際に測定されたフロー値は、真の測定可能な値とはっきりと区別するために、以下において
と指定すべきである。文献には、冗長なMFA問題の近似解
を計算するいくつかの方法がある。
δは正規分布していると想定される。期待ベクトル0および共分散行列
は、式
を有する(付属書の参考文献[23]を参照されたい)。Iは単位行列である。R'はRの誘導形であり、線形従属線を除去することによって作成され、従って、最大ランクを有する(これはユニークでない)。これは、対応する線形変換:
R'=ΓR (2.11)
を行う、非正方行列Γとの掛け算によって作成され得る。
が対角であり)、
の標準偏差が測定値
の大きさと比例し、一定の比例定数b>0であると仮定して、最適化問題(2.9)は(2.8)に等しい。以下:
が成り立つ。
推論された生化学的ネットワークモデルの質は以前の論評において扱われたことがない。しかしながら、その定式化における質が十分でないことでMFA結果に重大な欠陥が生じ得ることは明らかである。高い有意性と、できるだけ大きな単純化の中を取るものとして、意味のあるモデルを作り出すために検証方法が必要とされる。1994からのvan der Heijdenらの刊行物において、統計学的に動機づけられた検定が示されており、可能性のある系統的な誤差発生源の検出が説明されている(付属書の文献目録[40]を参照されたい)。これらの研究は、均衡化することができるフローの解析と、モデルにおける不一致への影響に基づいている。従って、これらの研究は、冗長な式が存在する場合にしか成り立たない。
ε≠0
が成り立つ。εの共分散行列Cεは、
によって計算され得る。従って、これはフラックス測定値の共分散行列に依存し、測定値の不確実性を考慮に入れている。検定のために、検定統計量Hεが使用され、その観測値は、
によって示される。
検定
H0:考慮された代謝モデルの不一致は有意でない
対
H1:考慮されている代謝モデルの不一致は有意である
は、(1-α)・χ2分布の100%分位数を示す。Rang(Cε)は自由度である。
前もって定義された検定が不一致を示すのであれば、これは、行列
に反映され、従って、検定の結果に影響を及ぼす測定ノイズの過小評価によるものかもしれない。3つのさらなる、可能性がある誤差発生源が付属書の文献[40]において議論されている。
一態様では、試料は、細胞培養物のバイオリアクター208の中で、細胞培養物の培養中のいくつかの時点で採取され、図2に示したように1つまたは複数の分析装置250に自動的に、または手作業で移される。分析装置は、1つまたは複数の分析器のシステム、例えば、Thomasチャンバーまたは細胞密度を求める光学的計数ステーションでもよい。例えば、高速液体クロマトグラフまたは当技術分野において公知の他の適切な方法が個々のアミノ酸の濃度を求めるのに使用され得る。試料は、例えば、24時間間隔で採取されてもよい。このように得られた測定データはデータ処理システム252に伝送される。データ処理システム252は、例えば、制御装置として1つまたは複数のバイオリアクターをモニタリングおよび/または制御するコンピュータでもよい。
データ処理システム252は、決められた時点で測定された入力パラメータ値(特に、細胞外代謝産物の濃度および細胞密度)に基づいて代謝モデル254の1種または複数種の細胞外フラックスを予測または推定するためにトレーニングされたことのあるMLP、例えば、ニューラルネットワーク(NN)またはいくつかのニューラルネットワークの協働システムを含む。例えば、データ処理システム252は、入力として、ある時点で得られた測定データを、将来の時点(例えば、翌日)での代謝状態を予測しようとする細胞培養物の細胞と同じタイプの細胞の細胞培養物から得られた試験データセットに対してトレーニングされたMLPに自動的に渡すすプログラムロジックを含んでもよい。
ニューラルネットワークを使用して、時点t(n)で現在測定されている細胞外代謝産物濃度
に基づいて、一段階予測において、自由に選択された将来の時点の時点t(n+1)での細胞外フラックスが予測または評価される。従って、この工程では、測定された入力パラメータ値がMLPに入力されたことに応じて、MLPは将来の時点の1種または複数種の細胞外フラックスを計算し、戻す。
ニューラルネットワークを介して推定された細胞外フラックスはまた、本発明の態様に従って、その時点t(n)での現在の濃度
に基づいて、任意に選択された将来の時点t(n+1)での代謝産物濃度を一段階予測するのにも使用される。
次に、将来の時点についてMLPによって予測された細胞外代謝産物の取り込みと放出速度(細胞外フラックス)に基づいて、代謝フラックス解析が本発明の態様に従って行われる。これも、代謝モデルにおいて定式化されたように細胞内フラックスと化学量論式を組み込んでいる。この代謝モデルでは細胞外フラックスと細胞内フラックスは1種または複数種の細胞内代謝産物を介して互いにつなげられているので、細胞外フラックス予測値に基づいて将来の時点での細胞内フラックスを予測することも数学的に可能である。代謝フラックス解析を行うための、対応するプログラムルーチンはMatlabおよび市販されている他のソフトウェアソリューションにおいて実行することができる。
は、測定値の質と、ニューラルネットワークによる測定値の推定の質が共分散行列に組み込まれていることを示している。しかしながら、リワーディングが推定を容易にしない可能性がある。なぜなら、差
は右辺の真ん中の項において独立確率変数と期待値を表さず、従って、右辺の真ん中の項にある期待値は未知だからである。従って、以下で説明される、付属書の式4.5に従う、この論文において使用した式は実際の共分散行列の概算にすぎない。
将来のマテリアルフローを予測するための反応速度モデルの仕様と比較して、トレーニング済みMLPの使用には、その作成が、通常、半自動方法では簡単であり、かつ速いという利点がある。トレーニングによってNNの形でMLPがどのように作成され得るかという例の1つを下記で説明する。
MLPトレーニングのために、できるだけ広いデータベースを得るために、いくつかのトレーニング培養物を、好ましくは、いくつかのバイオリアクターの中で培養する。好ましくは、これらのバイオリアクターは、1つまたは複数のフェドバッチリアクターと、他のリアクタータイプからの1つまたは複数の追加のバイオリアクターを備える。
● 第2のバイオリアクターでも最初にバッチ方法を使用した。指数増殖期の終わりに近づくにつれて部分回収を行い、リアクターを新鮮な培地で満たした。その結果、(体積および接種物質濃度の点で)発酵の始まりに似た条件を実現した。その後、バッチ方法を継続した(いわゆるスプリットバッチ方法).
● 残りの6つのリアクター内での発酵を初期バッチ期ではフェドバッチとして行った。発酵の終わりに近づくにつれて連続供給と、パルス様の栄養分添加を両方とも行った。フェドバッチ発酵は培地中のグルコース濃度の点で異なり、グルコース濃度を異なる供給戦略によって調整した。以下の説明では、フェドバッチアプローチに番号が付いていることが多い。この番号に従って、最初の2つのアプローチでは、プロセスの後半に短期間の完全グルコース制限を行った後に、大量投与添加を行った。第3および第4のアプローチには同じ制限があったが、その後の大量投与は、さらに高いグルコース濃度を設定した。第5および第6のアプローチでは、プラスの最小グルコース濃度が常にあった。
バイオリアクター内の現在の状態から、現在のサンプリング時間と次のサンプリング時間との間の細胞外代謝産物の平均フラックスを推定するために、ニューラルネットワークを作成した(いわゆる一段階予測)。一部の態様によれば、出力変数として、それぞれの細胞外マテリアルフローを得るために別々のネットワークをトレーニングした。現在主流となっている、いろいろな細胞外代謝産物濃度が入力変数として役立った。
本発明の好ましい態様によれば、入力パラメータ値を選択するために、代謝モデルにおいて述べた、いくつかのまたは全ての細胞外代謝産物は、それぞれのフラックスの予測のために関連性に従ってソートされる。好ましくは、冗長な情報内容を有する細胞外代謝産物は考慮されない。
イテレーションおよび隠れニューロンの数を決定した後に、それぞれのネットにおける初回学習工程のために出力層に値η=0.1、隠れ層にη=0.02を選択する。総イテレーション数の1/10になるたびに、ステップサイズを初期値の1/10にする。重みの初期値を
均等に分布する乱数にした。調整された値が経験的平均0と経験的標準偏差0.5を有するように、入力トレーニングデータと出力トレーニングデータを別々に代謝産物によって標準化した。トレーニングデータセットの平均値および標準偏差を用いて試験データを同じように変換した。
推定は、トレーニングデータセットおよび試験データセットへのデータのグループ化の点で異なる3つの異なるニューラルネットワークによって行った。
1. ネットワーク1:第1のネットワークのトレーニングデータセットは、フェドバッチ発酵およびバッチ発酵のうちの3つからのデータからなり、試験データセットは、他の3つのフェドバッチ発酵およびスプリットバッチ発酵からのデータからなった。
2. ネットワーク2:第2のネットワークは、トレーニングデータセットにフェドバッチ発酵のうちの3つからのデータと、試験データセットにフェドバッチ発酵、バッチ発酵、およびスプリットバッチ発酵のうちの残りの3つからのデータを有した。
3. ネットワーク3:第3のネットワークのトレーニングデータセットは、フェドバッチ発酵のうちの3つからのデータを含み、試験データセットは残りの3つのフェドバッチ発酵からのデータを含んだ。
この目的は、細胞代謝の最も正確な予測を可能にするトレーニング済みMLPを作成することである。
に基づいており、以下で説明される。
が成り立つ。時間tでの、ある成分の細胞外フラックスvは、付属書のセクション2.2.2において既に述べられているように、時間あたりの細胞によって吸収または放出された物質の量である。従って、これは、
によって示される。式中、VC(t)は、時間tでの反応培地中の生細胞の数である。従って、付属書の式2.1によれば、2つの不連続点間の領域において、
が成り立つ。
によって推定し得る。
に基づいて、任意の選択された将来の時点t(n+1)、例えば、翌日の細胞外代謝産物濃度を一段階予測するために使用され得る。
は、この
を用いた密度差によって補足される。式中、
n番目の時点でのLDH濃度が指定され、n0が発酵を開始した時点であれば、
が成り立つように、a'0が選択される。
を計算するために使用される。これには、MFAの結果としての、一致しないフラックス分布の割合が大幅に低下するという利点がある。
を定式化するために、好ましくは、モデル検証の結果を含める。これにより、どの標準偏差が、おおむね一貫したモデルにつながるか、場合によっては、どの測定値が信頼できないと分類すべきかが分かった。
1. セットχに、可能性がある全ての入力変数、ここでは、全ての現在の細胞外代謝産物濃度をまとめる。このセットuは、既に選択されている全ての入力変数を含んでいる(最初は空である)。Yは、現在のサンプリング時間と将来のサンプリング時間との間の出力変数、すなわち、j番目の細胞外代謝産物フラックスである。
2. 所定の標準化されたデータセットに基づいて、uにおける要素を考慮して、χとYにある可能性のある、それぞれの入力変数間の部分トランスインフォメーションのために近似値を計算する。例えば、標準化により、全ての入力変数と出力変数の平均値は0となり、標準偏差は0.5となり、従って、異なる桁のせいで変数の関連性をゆがめる効果が取り除かれる。好ましくは、この標準化は、トレーニングデータセットの測定値に関してニューラルネットワークをトレーニングする前に、ならびに現在得られた測定値に関してトレーニング済みMLPを用いて現在の測定値を検定手順に入力する時に行われる。
3. 最も高いPMI値を有する、χにある変数を記録する。これをuに追加し、χから除去する。
4. χが空になるまで工程2および3を繰り返す。
残差
および
式中、gは周辺分布または共通分布の密度関数である。残差は、Uにまだ含まれていないXとYの情報しか含んでいない。I'の値が大きければ大きいほど、依存性は強くなる。
(x(n),y(n))、n=1,...,N、XとYの標本対
(x(n)',y(n)')は、関連付けられたgx',y'分布の残差対である。そして次に、普通知られていない密度関数はコア密度推定量によって近似され得る。これもまたN個の標本からの情報を使用する。これらの推定量から、簡単に言うと、標本のヒストグラムに似た連続密度関数が得られる。これはN個のコア関数の重みつき重ね合わせによって生じる。そして次に、これらは、ベル形であり、かつ標本値のうちの1つについて対称的な密度関数である。
は、おおまかに推定されるはずである。この目的にナダラヤ・ワトソン推定量が使用されることがある。これは、以前に適用された原理に基づいており、以下のように導出することができる。
最初に、全ての測定可能で利用可能な入力代謝産物または少なくとも、細胞外代謝産物として代謝モデルにおいて生じた全ての代謝産物を「第1のセット」に移し、関連性に従ってソートする。すなわち、PMI基準を使用することによって、どの代謝産物が、出力のための最大の表現/予測能力(細胞外フラックスを予測しようとする細胞外代謝産物-「出力代謝産物」の速度)を持つかを確かめる。この代謝産物を「第2のセット」に移す。次に、「第2のセット」は実際の入力パラメータ値を提供する。この第1のソーティング工程において確かめた「関連性」または予測能力は第2のセットにある代謝産物に依存しない。
好ましくは、MLPトレーニングまたはMLP適用のための入力として全ての入力代謝産物を使用するわけではない(これは、オーバーフィッティングにより非常に不十分な予測能力をもたらす)。従って、最も関連性のある入力代謝産物のうちどれだけ使用すべきかを確かめる。これは、試験データセットから「x」個の最も関連性のある入力代謝産物の予測能力を確かめ、xを変え、次に、予測能力が最も良い数xを選択することによって行われる。
102~112 工程
200 システム
202 プロセッサ
204 バイオリアクター
206 バイオリアクター
208 バイオリアクター
210 測定値インターフェイス
212 記憶媒体
214 代謝モデル
216 参照値
218 機械学習ロジック
220 プログラムロジック
222 制御インターフェイス
224 ユーザーインターフェイス
250 代謝産物の濃度を求めるための装置
252 計算および予測のためのコンピュータシステム
254 フラックスのある代謝モデル
256 細胞外代謝産物の濃度測定値
402 代謝モデル
404 細胞外代謝産物
406 細胞内代謝産物
408 細胞外フラックス
410 細胞内フラックス
502 細胞外代謝産物の濃度測定値の経過
504 細胞外フラックス予測値の経過
506 MFAによって計算した細胞内フラックスの経過
802 フェドバッチバイオリアクター1のグラフ
804 フェドバッチバイオリアクター5のグラフ
806 バッチバイオリアクターのグラフ
808 スプリットバッチバイオリアクターのグラフ
1502 出力パラメータ値
1504 入力パラメータ値
Claims (20)
- 以下の工程を含む、特定の細胞タイプの細胞の細胞培養物の代謝状態を予測するための方法:
- 該特定の細胞タイプの細胞の代謝モデル(402)を準備する工程(102)であって、該代謝モデルが、複数種の細胞内代謝産物および細胞外代謝産物ならびに複数種の細胞内フラックス(410)および細胞外フラックス(408)を含み、該代謝モデルが、該細胞内代謝産物(406)のうちの1つと該細胞外代謝産物(404)のうちの1つとの間の少なくとも1つの化学量論的関係を指定する化学量論式を含む、工程;
- 該細胞培養物の培養中の複数の時点のそれぞれにおける、
○該時点で測定された複数の測定値を受け取る工程(106)であって、該測定値が、該細胞培養物の培養培地における該代謝モデルの複数種の細胞外代謝産物の濃度と、該細胞培養物における該細胞の細胞密度測定値とを含む、工程;
○受け取った該測定値を、入力パラメータ値として、トレーニング済みの機械学習プログラムロジック(MLP)(218)に入力する工程(108);
○受け取った該測定値を用いて、該MLPによって、将来の時点での細胞外代謝産物の細胞外フラックス(408)を予測する工程(110)であって、該将来の時点が、該測定値を受け取った時点より後の時点であり、該細胞外フラックスが、細胞内への細胞外代謝産物の取り込み速度および/または細胞から培地への細胞外代謝産物の放出速度である、工程;
○該細胞外代謝産物の該細胞外フラックスの予測値と該代謝モデルの該化学量論式とを用いて、該将来の時点での細胞内フラックスを計算するために代謝フラックス解析を行う工程(112)。 - 前記方法が、機械学習によりMLPを作成する工程をさらに含み、該作成する工程が、
- トレーニングデータセットを作成することであって、
○前記特定の細胞タイプの細胞の少なくとも1つのトレーニング細胞培養物であって、該トレーニングデータセットを作成するためにデータが収集される少なくとも1つのトレーニング細胞培養物の培養中の複数のトレーニング時点であって、トレーニングデータが収集される時点である複数のトレーニング時点のそれぞれにおいて、
■該トレーニング時点で測定された複数の測定値を受け取る工程であって、該測定値が、該少なくとも1つのトレーニング細胞培養物の培養培地における代謝モデルの複数種の細胞外代謝産物の濃度と、該少なくとも1つのトレーニング細胞培養物における細胞の細胞密度測定値を含む、工程、
■現在のトレーニング時点の時間表示を受け取る工程、および
■その時点で受け取った測定値と、それぞれの前の時点で受け取った測定値との関数として、該細胞外代謝産物の細胞外フラックスを計算する工程であって、該細胞外フラックスが、細胞内への該細胞外代謝産物の取り込み速度および/または培地への該細胞外代謝産物の放出速度である、工程
を含む、トレーニングデータセットを作成すること、
- 該MLPをトレーニングすることであって、
○該トレーニング時点のそれぞれにおいて受け取った測定値を、入力パラメータ値として該MLPに入力し、そのトレーニング時点の後続の各時点において該後続の時点について計算された細胞外フラックスを、該入力パラメータ値に関連する出力パラメータ値として該MLPに入力する工程;ならびに
○該MLPが、該入力パラメータ値に基づいて、それぞれの該関連する出力パラメータ値を予測するように学習するように、該MLPによる学習プロセスを実施する工程
を含む、該MLPをトレーニングすること、
- トレーニングされた該MLPを揮発性または非揮発性の記憶媒体に記憶すること
を含む、請求項1記載の方法。 - 前記特定の細胞タイプの細胞の複数のトレーニング細胞培養物の培養中の前記複数のトレーニング時点のそれぞれにおいて、前記測定値および時間表示が受け取られ、かつ前記細胞外代謝産物の細胞外フラックスが計算されるように、前記トレーニングデータセットが作成され、ここで、該細胞培養物が、異なるタイプのバイオリアクターの中で培養され、該バイオリアクターのタイプが、フェドバッチバイオリアクター、バッチバイオリアクター、灌流リアクター、ケモスタット、およびスプリットバッチバイオリアクターを含むバイオリアクタータイプのうちの少なくとも2つの異なるバイオリアクタータイプを含む、請求項2記載の方法。
- 前記MLPが、1つのニューラルネットワークであるか、またはいくつかのニューラルネットワークからなるシステムである、請求項1~3のいずれか一項記載の方法。
- 前記複数種の細胞外代謝産物の濃度測定値がそれぞれ、
- 該代謝産物の体積当たりの含有量、特に、質量濃度もしくは物質量濃度、または
- 該体積当たりの含有量と直線的もしくは問題となっている濃度範囲の少なくとも90%であるようにほぼ直線的に相関する値、特に、補正もしくは標準化された、該細胞外代謝産物の代謝産物濃度測定値もしくはフラックス測定値
である、請求項1~4のいずれか一項記載の方法。 - 前記MLPが、いくつかのサブMLPからなるシステムであり、該システムに含まれる個々のサブMLPがそれぞれ、単一の細胞外代謝産物の細胞外フラックスを予測するようにトレーニングされており、かつ、該単一の細胞外代謝産物の将来の時点での細胞外フラックスを予測するのに選択的に用いられる、請求項1~5のいずれか一項記載の方法。
- 前記MLPが、前記細胞外代謝産物のうちの1つの細胞外フラックスを予測するために、複数種の細胞外代謝産物の濃度測定値を入力パラメータ値として使用し、ここで、細胞外フラックスを求めようとする該細胞外代謝産物のうちの少なくとも2つについて、入力パラメータ値として使用する該複数種の細胞外代謝産物が異なる、請求項1~6のいずれか一項記載の方法。
- - 全ての入力候補代謝産物の濃度を複数の時点にわたって測定する工程であって、該入力候補代謝産物が、前記特定のタイプの細胞培養物を有する参照バイオリアクター内で測定可能な程度に利用可能な全ての細胞外代謝産物を含むか、または前記代謝モデルの全ての細胞外代謝産物を含む、工程;
- 細胞外フラックスを予測しようとする単一の細胞外代謝産物のそれぞれについて、該単一の代謝産物の細胞外フラックスを予測するために入力パラメータ値として使用される複数種の細胞外代謝産物を特定するための選択プロセスを行う工程であって、該選択プロセスが、該単一の代謝産物それぞれに関して、
(a)細胞外代謝産物の第1のセットを定める工程であって、該第1のセットが全ての候補入力代謝産物を含む、工程;
(b)該第1のセットにおける各細胞外代謝産物の第1の関連性スコアを、その代謝産物の濃度測定値の関数として計算する工程であって、該第1の関連性スコアが、該単一の細胞外代謝産物の細胞外フラックスに関する、それぞれの細胞外代謝産物の濃度の予測能力を示す、工程;
(c)該細胞外代謝産物のうち最も高い第1の関連性スコアを有するものだけを、該第1のセットから、まだ空である細胞外代謝産物の第2のセットに移し、この代謝産物を該第1のセットから除去する工程;
(d)該第1のセットにおける各細胞外代謝産物のさらなる関連性スコアを、その代謝産物の濃度測定値と、該第2のセットに含まれる全ての細胞外代謝産物の濃度測定値との関数として計算する工程であって、該さらなる関連性スコアが、既に該第2のセット含まれている代謝産物を考慮した、該単一の細胞外代謝産物の細胞外フラックスに関する、それぞれの候補入力代謝産物の濃度の予測能力を示す、工程;
(e)該第1のセットの細胞外代謝産物のうち最も高いさらなる関連性スコアを有するものだけを該第2のセットに移し、この代謝産物を該第1のセットから除去する工程であって、この代謝産物を含めることによって、該単一の細胞外代謝産物の細胞外フラックスの予測に関する、第2のセットに含まれる代謝産物の情報冗長性の最大限度を第2のセットが超えない場合にだけ、この代謝産物が移される、工程;
(f)該第2のセットが該情報冗長性の最大限度を超えることなしでは、これ以上の代謝産物を該第1のセットから該第2のセットに移すことができなくなるまで、工程(d)および(e)を繰り返す工程;ならびに
(g)該単一の細胞外代謝産物の細胞外フラックスを予測するために、該第2のセットに移された代謝産物だけを選択的に入力パラメータ値として使用する工程
を含む、選択プロセスを行う工程
をさらに含む、請求項7記載の方法。 - - 前記第1のセットの代謝産物と、細胞外フラックスを予測しようとする前記単一の代謝産物との間の部分相互情報量スコア(PMIスコア)として、前記第1の関連性スコアが計算され;
- 既に前記第2のセットに含まれている全ての代謝産物を考慮して、該第1のセットの代謝産物と、細胞外フラックスを予測しようとする該単一の代謝産物との間のPMIスコアとして、第2の関連性スコアであって、該第1のセットにおける各細胞外代謝産物の前記さらなる関連性スコアである第2の関連性スコアが計算される、
請求項8記載の方法。 - 前記MLPが、単一の代謝産物のそれぞれの細胞外フラックスを予測するために、複数種の細胞外代謝産物の濃度測定値を入力パラメータ値として使用し、該複数種の細胞外代謝産物が、少なくとも1種のアミノ酸、好ましくは少なくとも2種のアミノ酸を含む、請求項1~9のいずれか一項記載の方法。
- 前記複数の時点が、10分~48時間、好ましくは1~24時間の時間間隔だけ隔てられている、請求項1~10のいずれか一項記載の方法。
- 前記細胞タイプが、生体分子を得る目的でバイオリアクター内で維持されかつ/または増殖する遺伝子組換えされた細胞タイプである、請求項1~11のいずれか一項記載の方法。
- 前記将来の時点のそれぞれにおいて、前記代謝モデルの細胞内フラックスのうちのいくつかまたは全ての計算が行われる、請求項1~12のいずれか一項記載の方法。
- - それぞれの参照値または参照値範囲から閾値を超えて逸脱している細胞内フラックス(410)計算値を有する、前記代謝モデルの1種または複数種の細胞内代謝産物を特定する工程;
- 細胞増殖または所望の生体分子の産生の制限要因として作用する細胞内フラックスを自動的に特定する工程
をさらに含む、請求項1~13のいずれか一項記載の方法。 - 以下の工程を含む、特定の細胞タイプの細胞の細胞培養物を含むバイオリアクターをモニタリングおよび/または制御するための方法:
- 請求項1~14のいずれか一項記載の方法に従って将来の時点での該細胞の細胞内フラックスを計算する工程;
- 該細胞内フラックスの予測値を、それぞれの1種もしくは複数種の細胞内代謝産物の許容可能な細胞内フラックスの参照値もしくは参照値範囲と比較する工程;
- 少なくとも1つの細胞内代謝産物の細胞内フラックスの計算値の、そのそれぞれの参照値もしくは参照値範囲からの逸脱が限界値を超えた場合に、警告を発する工程;ならびに/または
- 該逸脱を小さくするように該バイオリアクターの状態もしくはその中に含まれる培地を変える工程を自動的に開始するための制御コマンドを該バイオリアクターに送る工程であって、該自動工程が、特に培養培地の量および/もしくは組成を変えることを含む、工程。 - - 請求項14に従って、細胞増殖もしくは所望の生体分子の産生の制限要因として作用する、細胞の代謝モデル内の反応を特定する工程;および
- 細胞増殖もしくは該生体分子の産生または該生体分子の質が促進されるように、制限要因として作用する細胞内フラックスを改変する物質を選択的に自動的に追加するか、あるいは、このような追加の要求をユーザーインターフェイスを介して発行する、工程
をさらに含む、請求項15記載の方法。 - 以下の工程を含む、特定の細胞タイプの細胞の代謝的に好ましいクローンを特定するための方法:
- 複数のバイオリアクターの中で異なる細胞培養物を培養する工程であって、該異なる細胞培養物が、該特定の細胞タイプの細胞の異なるクローンを含む、工程;
- 各細胞クローンについて別々に複数の時点で、請求項1~16のいずれか一項に従って、1種または複数種の細胞内代謝産物の細胞内フラックスを計算する工程;
- 代謝的に最も好ましい、該1種または複数種の細胞内代謝産物の細胞内フラックス計算値を有する、該細胞クローンのうちの1つを特定する工程。 - - 1種または複数種の細胞外代謝産物の現在の細胞外フラックスを、現在の時点および以前の時点で測定された該細胞外代謝産物の濃度から計算する工程;
- 該細胞外代謝産物の該現在の細胞外フラックスの計算値と代謝モデルの化学量論式とを用いて、さらなる代謝フラックス解析を行って、現在の時点での現在の細胞内フラックス(410)を計算する工程;ならびに
- 該現在の細胞内フラックスの計算値を、細胞培養物の細胞の現在の代謝状態の特徴付けとして使用する工程
をさらに含む、請求項1~17のいずれか一項記載の方法。 - - 前記細胞培養物の培養中の前記時点のそれぞれにおいて、該細胞培養物の培地中の乳酸デヒドロゲナーゼ(LDH)濃度の測定値がさらに受け取られ;
- 前記将来の時点のそれぞれにおける前記細胞外代謝産物の細胞外フラックスの予測が、前記細胞密度測定値ではなく補正された細胞密度を用いて、前記MLPによって行われ、該補正された細胞密度の計算が、該時点のそれぞれにおける、
○該LDH濃度測定値の関数として、該細胞培養物の培地中の溶解細胞の密度を計算する工程であって、該関数が、該培地中のLDH濃度と前記特定の細胞タイプの溶解細胞の数との関係を表す、経験的に求められたヒューリスティックな一次関数である、工程;
○該培地中の該細胞密度測定値と該溶解細胞の密度計算値との和として、該補正された細胞密度を計算する工程
を含む、
請求項1~18のいずれか一項記載の方法。 - 以下を備える、特定の細胞タイプの細胞の細胞培養物の代謝状態を予測するためのシステム(200):
- 1つまたは複数のプロセッサ(202)、
- 該細胞培養物を含むバイオリアクター(204、206、208)からの測定値を受け取るための第1のインターフェイス(210)、ならびに
- 以下を備える、揮発性または非揮発性の記憶媒体(212):
○該特定の細胞タイプの細胞の代謝モデル(402)であって、複数種の細胞内代謝産物および細胞外代謝産物ならびに複数種の細胞内フラックス(410)および細胞外フラックス(408)を含み、該細胞内代謝産物(406)のうちの1つと該細胞外代謝産物(404)のうちの1つとの間の少なくとも1つの化学量論的関係を指定する化学量論式を含む、代謝モデル(402);
○トレーニング済みの機械学習プログラムロジック(MLP)(218);
○該細胞培養物の培養中の複数の時点のそれぞれにおいて以下の工程を含む方法を行うように構成された、プログラムロジック(220):
■第1のインターフェイスを介して該時点で測定された複数の測定値を受け取る工程(106)であって、該測定値が、該細胞培養物の培養培地における該代謝モデルの複数種の細胞外代謝産物の濃度と、該細胞培養物における該細胞の細胞密度測定値とを含む、工程;
■受け取った該測定値を入力パラメータ値として該MLPに入力する工程(108);
■受け取った該測定値を用いて、該MLPによって、将来の時点での該細胞外代謝産物の細胞外フラックス(408)を予測する工程(110)であって、該将来の時点が、該測定値を受け取った時点より後の時点であり、該細胞外フラックスが、細胞内への該細胞外代謝産物の取り込み速度および/または細胞から培地への該細胞外代謝産物の放出速度である、工程;
■該細胞外フラックス予測値と該代謝モデルの該化学量論式とを用いて、該将来の時点での細胞内フラックスを計算するために代謝フラックス解析を行う工程(112)。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
EP17211217.9 | 2017-12-29 | ||
EP17211217 | 2017-12-29 | ||
PCT/EP2019/050006 WO2019129891A1 (de) | 2017-12-29 | 2019-01-02 | Vorhersage des metabolischen zustands einer zellkultur |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2021508872A JP2021508872A (ja) | 2021-03-11 |
JP7092879B2 true JP7092879B2 (ja) | 2022-06-28 |
Family
ID=61017738
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2020534849A Active JP7092879B2 (ja) | 2017-12-29 | 2019-01-02 | 細胞培養物の代謝状態の予測 |
Country Status (5)
Country | Link |
---|---|
US (2) | US20200377844A1 (ja) |
EP (1) | EP3732485A1 (ja) |
JP (1) | JP7092879B2 (ja) |
CN (1) | CN112119306B (ja) |
WO (1) | WO2019129891A1 (ja) |
Families Citing this family (26)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
MX2021011623A (es) * | 2019-03-29 | 2021-10-13 | Amgen Inc | Prediccion del rendimiento de cultivos celulares en biorreactores. |
CN113994281A (zh) * | 2019-07-03 | 2022-01-28 | 富士胶片株式会社 | 优化支援装置、方法及程序 |
JPWO2021059577A1 (ja) * | 2019-09-24 | 2021-04-01 | ||
EP4015615A4 (en) * | 2019-09-24 | 2022-10-12 | FUJIFILM Corporation | INFORMATION PROCESSING DEVICE, INFORMATION PROCESSING METHOD AND INFORMATION PROCESSING PROGRAM |
WO2021076855A1 (en) * | 2019-10-18 | 2021-04-22 | Janssen Biotech, Inc. | Dynamic monosaccharide control processes |
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 |
AU2021355381A1 (en) * | 2020-10-01 | 2023-05-04 | Amgen Inc. | Predictive modeling and control of cell culture |
EP3979010B1 (en) * | 2020-10-02 | 2023-07-12 | Sartorius Stedim Data Analytics AB | Monitoring and control of bioprocesses |
CN112767995B (zh) * | 2020-12-16 | 2023-05-23 | 广东省微生物研究所(广东省微生物分析检测中心) | 一种微生物次级代谢产物制备过程中产物产率评估方法 |
US20220237471A1 (en) * | 2021-01-22 | 2022-07-28 | International Business Machines Corporation | Cell state transition features from single cell data |
US20220282199A1 (en) * | 2021-03-03 | 2022-09-08 | Applied Materials, Inc. | Multi-level machine learning for predictive and prescriptive applications |
EP4339958A1 (en) * | 2021-05-13 | 2024-03-20 | Shimadzu Corporation | Antibody production assistance method and antibody production assistance program |
US20240242781A1 (en) * | 2021-05-13 | 2024-07-18 | Shimadzu Corporation | Antibody Production Assisting Method and Antibody Production Assisting Program |
EP4116403A1 (en) * | 2021-07-07 | 2023-01-11 | Sartorius Stedim Data Analytics AB | Monitoring, simulation and control of bioprocesses |
CN113838530B (zh) * | 2021-08-18 | 2022-07-15 | 江南大学 | 利用细胞代谢网络监测生物制造过程的方法 |
KR20230060189A (ko) * | 2021-10-27 | 2023-05-04 | 프레스티지바이오로직스 주식회사 | 인공지능을 이용하여 세포 배양조건을 결정하기 위한 장치 및 장치의 동작 방법 |
CN114036810A (zh) * | 2021-11-04 | 2022-02-11 | 江南大学 | 一种细胞培养状态在线估计及优化补料调控方法 |
CN114121163B (zh) * | 2021-11-30 | 2023-10-27 | 深圳太力生物技术有限责任公司 | 基于整合学习的培养基预测系统、训练及培养基预测方法 |
US20230281362A1 (en) * | 2022-01-31 | 2023-09-07 | X Development Llc | Parameter and state initialization for model training |
WO2023167802A1 (en) * | 2022-03-01 | 2023-09-07 | Amgen Inc. | Hybrid predictive modeling for control of cell culture |
WO2023190137A1 (ja) * | 2022-03-30 | 2023-10-05 | 富士フイルム株式会社 | 培地組成の探索及び細胞特性の推定を行う方法、装置、及びプログラム |
WO2024062609A1 (ja) * | 2022-09-22 | 2024-03-28 | 日本電信電話株式会社 | 異常検知装置、方法およびプログラム |
WO2024064890A1 (en) * | 2022-09-23 | 2024-03-28 | Metalytics, Inc. | Using the concepts of metabolic flux rate calculations and limited data to direct cell culture. media optimization and enable the creation of digital twin software platforms |
WO2024095564A1 (ja) * | 2022-10-31 | 2024-05-10 | 株式会社島津製作所 | 代謝フラックス分布の推定方法、解析装置、プログラム |
CN117235673B (zh) * | 2023-11-15 | 2024-01-30 | 中南大学 | 一种细胞培养预测方法及装置、电子设备、存储介质 |
CN117598391B (zh) * | 2024-01-18 | 2024-04-16 | 华农恒青科技股份有限公司 | 发酵工艺条件的确定方法以及饲料发酵过程监测方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2005034149A (ja) | 2003-06-30 | 2005-02-10 | Ajinomoto Co Inc | 同位体で標識された基質を用いた細胞内代謝フラックスの解析方法 |
JP2005293394A (ja) | 2004-04-02 | 2005-10-20 | Ajinomoto Co Inc | 代謝フラックスの決定方法 |
JP2007047994A (ja) | 2005-08-09 | 2007-02-22 | Ajinomoto Co Inc | 物質の生産プロセスのシミュレーション方法 |
JP2013085516A (ja) | 2011-10-18 | 2013-05-13 | Hitachi Plant Technologies Ltd | 細胞培養制御方法、細胞培養制御装置及びこれを備える細胞培養装置 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4612864B2 (ja) * | 2005-05-19 | 2011-01-12 | キヤノン株式会社 | 通信装置、当該装置における受信方法、コーデック、デコーダ、通信モジュール、通信部及びデコード方法 |
WO2010104963A1 (en) * | 2009-03-10 | 2010-09-16 | Bayer Healthcare Llc | Method for monitoring cell culture |
CN102175809A (zh) * | 2011-01-26 | 2011-09-07 | 中国药科大学 | 细胞代谢物相对含量作为细胞数量指标进行数据校正的新方法 |
JP6581975B2 (ja) * | 2013-06-26 | 2019-09-25 | バリタセル リミテッド | 細胞の特徴を判定または予測する方法 |
-
2019
- 2019-01-02 EP EP19700004.5A patent/EP3732485A1/de active Pending
- 2019-01-02 CN CN201980015112.5A patent/CN112119306B/zh active Active
- 2019-01-02 WO PCT/EP2019/050006 patent/WO2019129891A1/de unknown
- 2019-01-02 JP JP2020534849A patent/JP7092879B2/ja active Active
-
2020
- 2020-06-22 US US16/907,786 patent/US20200377844A1/en active Pending
-
2023
- 2023-06-08 US US18/331,411 patent/US20230313113A1/en active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2005034149A (ja) | 2003-06-30 | 2005-02-10 | Ajinomoto Co Inc | 同位体で標識された基質を用いた細胞内代謝フラックスの解析方法 |
JP2005293394A (ja) | 2004-04-02 | 2005-10-20 | Ajinomoto Co Inc | 代謝フラックスの決定方法 |
JP2007047994A (ja) | 2005-08-09 | 2007-02-22 | Ajinomoto Co Inc | 物質の生産プロセスのシミュレーション方法 |
JP2013085516A (ja) | 2011-10-18 | 2013-05-13 | Hitachi Plant Technologies Ltd | 細胞培養制御方法、細胞培養制御装置及びこれを備える細胞培養装置 |
Non-Patent Citations (3)
Title |
---|
Nolan R.P. et al.,Dynamic model of CHO cell metabolism,Metabolic Engineering,2010年10月07日,検索日 : 令和3年10月20日, URL : https://doi.org/10.1016/j.ymben.2010.09.003,DOI : 10.1016/j.ymben.2010.09.003 |
Panagiotou, G. et al.,Monitoring novel metabolic pathways using metabolomics and machine learning: induction of the phosphoketolase pathway in Aspergillus nidulans cultivations,Metabolomics,2007年12月,Vol. 3, No. 4,p.503-516,検索日 : 令和3年10月20日, URL : https://doi.org/10.1007/s11306-007-0061-7,DOI : 10.1007/s11306-007-0061-7 |
Wu S.G. et al.,Rapid Prediction of Bacterial Heterotrophic Fluxomics Using Machine Learning and Constraint Programming,PLOS Computational Biology,2016年04月19日,p.1-22,検索日 : 令和3年10月20日, URL : https://doi.org/10.1371/journal.pcbi.1004838,DOI : 10.1371/journal.pcbi.1004838 |
Also Published As
Publication number | Publication date |
---|---|
CN112119306B (zh) | 2024-09-03 |
JP2021508872A (ja) | 2021-03-11 |
WO2019129891A1 (de) | 2019-07-04 |
EP3732485A1 (de) | 2020-11-04 |
CN112119306A (zh) | 2020-12-22 |
US20230313113A1 (en) | 2023-10-05 |
US20200377844A1 (en) | 2020-12-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP7092879B2 (ja) | 細胞培養物の代謝状態の予測 | |
Narayanan et al. | A new generation of predictive models: The added value of hybrid models for manufacturing processes of therapeutic proteins | |
JP7554774B2 (ja) | 生物工学的生産を最適化するための方法および手段 | |
US11636917B2 (en) | Simulating the metabolic pathway dynamics of an organism | |
US20190219992A1 (en) | Controlling and monitoring a process to produce a chemical, pharmaceutical or biotechnological product | |
US20200202051A1 (en) | Method for Predicting Outcome of an Modelling of a Process in a Bioreactor | |
CN115151869A (zh) | 细胞代谢状态观测的计算机实现方法、程序和混合系统 | |
US20230323275A1 (en) | Monitoring and control of bioprocesses | |
Ohadi et al. | Development of a soft‐sensor based on multi‐wavelength fluorescence spectroscopy and a dynamic metabolic model for monitoring mammalian cell cultures | |
KR102690117B1 (ko) | 세포 배양 프로세스에서 프로세스 변수를 결정하기 위한 방법 | |
Hernández Rodríguez et al. | Predicting industrial‐scale cell culture seed trains–A Bayesian framework for model fitting and parameter estimation, dealing with uncertainty in measurements and model parameters, applied to a nonlinear kinetic cell culture model, using an MCMC method | |
JP2023544029A (ja) | 細胞培養の予測モデル化及び制御 | |
TW201643744A (zh) | 用於產生發酵模式之電腦實施方法 | |
US20210257045A1 (en) | Method for verifying cultivation device performance | |
WO2022168774A1 (ja) | 推定装置、学習装置、最適化装置、推定方法、学習方法、および最適化方法 | |
EP4296350A1 (en) | A concept for training and using at least one machine-learning model for modelling kinetic aspects of a biological organism | |
TW202346567A (zh) | 用於控制細胞培養之混合預測建模 | |
Yatipanthalawa et al. | Predictive models for upstream mammalian cell culture development-A review | |
Maton et al. | Selection of a minimal suboptimal set of EFMs for dynamic metabolic modelling | |
Gao et al. | DMFA-based operation model for fermentation processes | |
González-Vargas et al. | Validation methods for population models of gene expression dynamics | |
Kontoravdi | Systematic methodology for the development of mathematical models for biological processes | |
US20240170098A1 (en) | Method for predicting cell culture result, culture result prediction program, and culture result prediction apparatus | |
Vásquez | Strategic Reparametrization as a Model Enhancer: Application to Wine Fermenter Digital Twins | |
Torrealba Vásquez | Strategic reparametrization as a model enhancer: application to wine fermenter digital twins |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20200820 |
|
RD04 | Notification of resignation of power of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7424 Effective date: 20210225 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20211025 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20220124 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20220530 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20220616 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 7092879 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |