(57)【要約】【課題】比較的少数個の内部状態変数を定め、これらの内部状態変数で細胞の内部を記述する数学的モデルを組み立てることにより、観測可能な外部変数とあわせて、培養槽内の状態をより精密に制御することが可能な培養槽の運転制御装置を提供する。【解決手段】動物細胞又は微生物を培養する培地が封入された培養槽１と、培地中の栄養成分濃度，生産物濃度，細胞数密度を測定する計測装置と、培地中に栄養成と酸素を補給する供給装置と、計測装置の計測値を入力し供給装置を制御する演算処理装置７とを備え、演算処理装置７は、栄養成分濃度，生産物濃度，細胞数密度の時間変化率を入力データとして、記憶装置８に記憶された細胞内反応速度の回路網から導かれる方程式を解いて細胞内反応速度を計算し、計算結果により供給装置を制御して培地中の成分濃度を制御する。【選択図】図１
2007202500.tif 000002
【特許請求の範囲】
【請求項１】
動物細胞又は微生物を培養する培地が封入された培養槽と、前記培地中の栄養成分濃度，生産物濃度，細胞数密度を測定する計測装置と、前記培地中に栄養成分と酸素を補給する供給装置と、前記計測装置の計測値を入力し前記供給装置を制御する演算処理装置とを備え、該演算処理装置は、栄養成分濃度，生産物濃度，細胞数密度の時間変化率を入力データとして、記憶装置に記憶された細胞内反応速度の回路網から導かれる方程式を解いて細胞内反応速度を計算し、該計算結果により前記供給装置を制御して前記培地中の成分濃度を制御する培養槽の運転制御装置。
【請求項２】
動物細胞又は微生物を培養する培地が封入された培養槽と、前記培地中の栄養成分濃度，生産物濃度，細胞数密度を測定する計測装置と、前記培地中に栄養成分と酸素を補給する供給装置と、前記計測装置の計測値を入力し前記供給装置を制御する演算処理装置とを備え、該演算処理装置は、生成反応速度などの目標条件を定め、栄養成分濃度，生産物濃度，細胞数密度の時間変化率を入力データとして、記憶装置に記憶された細胞内反応速度の回路網から導かれる方程式を解いて前記目標条件を満足する細胞内反応速度の組を計算し、該計算結果により前記供給装置を制御して前記培地中の成分濃度を制御する培養槽の運転制御装置。
【請求項３】
前記細胞内反応速度の回路網は、細胞内の反応速度を細胞の内部状態変数として細胞内反応速度を記述する数学的モデルを組み立てたものである請求項１又は２に記載の培養槽の運転制御装置。
【請求項４】
前記記憶装置には、回路網の各々の物質名と化学反応を表す矢印の繋がりが記憶されている請求項３に記載の培養槽の運転制御装置。
【請求項５】
出力表示装置を具備し、前記栄養成分濃度，生産物濃度，細胞数密度の時間変化率と、細胞内反応速度の組からなる数値ベクトルを画面に表示する請求項１又は２に記載の培養槽の運転制御装置。
【請求項６】
前記演算処理装置は、栄養成分濃度，生産物濃度，細胞数密度の関数として記述するモデル推定式を記憶するものであって、栄養成分濃度，生産物濃度，細胞数密度の時間変化率を、前記モデル推定式中のモデルパラメータを計測値からフィッティングして推定する請求項１又は２に記載の培養槽の運転制御装置。
【請求項７】
栄養成分濃度，生産物濃度，細胞数密度の時間変化率を与えるモデル推定式を時間積分して、時間的未来の栄養成分濃度，生産物濃度，細胞数密度を予測し、該予測値と実測値を出力表示装置の画面に表示する請求項６に記載の培養槽の運転制御装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、動物細胞あるいは微生物などを培養して有用生産物を収穫する培養槽の運転制御装置に関する。
【背景技術】
【０００２】
動物細胞あるいは微生物などの培養は、工業的には、培地中に糖やアミノ酸などの栄養成分および酸素などを供給しつつ培養液の攪拌を行って細胞を増殖し、細胞が生産する有用物質を収穫している。
【０００３】
例えば、（特許文献１）には、動物細胞の培養を行うに際し、第一段階として、培養中の細胞密度、及び細胞活性度の測定値から培養層に供給する全ガス量をファジイ推論により決定し、第二段階として、培養液の溶存酸素密度、及びｐＨの測定値からガス組成をファジイ推論により決定し、第一段階と第二段階で求めた値から各成分のガス量を決定して培養槽に供給する動物細胞の培養制御方法が記載されている。
【０００４】
（特許文献２）には、基質と微生物を含む水系内における微生物個数濃度，基質濃度及び溶存酸素濃度に関する非定常非線形現象を微分方程式で記述し、この微分方程式を差分方程式に変換し、この差分方程式で用いられるパラメータを変数とする微分方程式の数値解を求め、基質物質に対する微生物の分解能力を測定し、差分方程式の数値解と基質物質に対する微生物の分解能力の測定値を比較し、測定値に最も近い数値解を水槽内における微生物個数濃度，基質濃度及び溶存酸素濃度の挙動として測定する浄化反応の予知方法が記載されている。
【０００５】
（特許文献３）には、酵母の製造において、菌体性能の定量的数値の目標値について予め決定し、それに基づいて酵母の成長速度を律する栄養源の培養への供給及び温度調整等の各種操作を数学的モデルを使用して制御する酵母の培養方法が記載されている。
【０００６】
【特許文献１】特開平５−１０３６６３号公報
【特許文献２】特開２００２−９５４９６号公報
【特許文献３】特開平９−６５８７３号公報
【発明の開示】
【発明が解決しようとする課題】
【０００７】
培養には細胞の増殖と生産物生成に好適な槽内条件を整える必要があり、(特許文献１)，（特許文献２），（特許文献３）に開示の技術は、培地中の栄養成分濃度，酸素濃度，攪拌回転数，細胞密度，微生物個数濃度，溶存酸素濃度等の計測できるものを計測して、制御している。しかし、本来の物質生産プロセスは細胞内の代謝反応によっており、細胞内部の内部状態は直接に計測制御できないブラックボックスとなっている。
【０００８】
このため、栄養成分濃度，生産物濃度，細胞数密度などの細胞外部の計測情報だけでは、しばしば、好適な制御を行うことができない状況が生じる。このように、従来の技術では、経験に頼った運転、あるいはこれまでの経験をルールとして計算機に覚えこませ、観測データに応じて事象推定と制御ルールの適用を行うルールベース制御などが用いられているものであった。
【０００９】
本発明の第１の目的は、比較的少数個の内部状態変数を定め、これらの内部状態変数で細胞の内部を記述する数学的モデルを組み立てることにより、観測可能な外部変数とあわせて、培養槽内の状態をより精密に制御することが可能な培養槽の運転制御装置を提供することにある。
【００１０】
本発明の第２の目的は、栄養成分濃度，生産物濃度，細胞数密度などの外部観測変数とこれらの内部状態変数の組からなる数値ベクトルを培養槽内の状態変数として用いて、好適な培養槽の運転条件を見出す計算アルゴリズムを備える培養槽の運転制御装置を提供することにある。
【課題を解決するための手段】
【００１１】
上記の目的を達成するために、本発明の培養槽の運転制御装置は、動物細胞又は微生物を培養する培地が封入された培養槽と、培地中の栄養成分濃度，生産物濃度，細胞数密度を測定する計測装置と、培地中に栄養成分と酸素を補給する供給装置と、計測装置の計測値を入力し前記供給装置を制御する演算処理装置とを備え、演算処理装置は、栄養成分濃度，生産物濃度，細胞数密度の時間変化率を入力データとして、記憶装置に記憶された細胞内反応速度の回路網から導かれる方程式を解いて細胞内反応速度を計算し、該計算結果により前記供給装置を制御して前記培地中の成分濃度を制御するものである。
【００１２】
又、動物細胞又は微生物を培養する培地が封入された培養槽と、培地中の栄養成分濃度，生産物濃度，細胞数密度を測定する計測装置と、培地中に栄養成分と酸素を補給する供給装置と、計測装置の計測値を入力し前記供給装置を制御する演算処理装置とを備え、演算処理装置は、生成反応速度などの目標条件を定め、栄養成分濃度，生産物濃度，細胞数密度の時間変化率を入力データとして、記憶装置に記憶された細胞内反応速度の回路網から導かれる方程式を解いて前記目標条件を満足する細胞内反応速度の組を計算し、該計算結果により前記供給装置を制御して前記培地中の成分濃度を制御するものである。
【発明の効果】
【００１３】
本発明により、栄養成分濃度，生産物濃度，細胞数密度などの培地中観測変数と細胞内反応速度を培養槽の状態変数として用いて、好適な培養槽の運転制御ができる。
【発明を実施するための最良の形態】
【００１４】
本発明の一実施例を図１から図４により説明する。図１は、本実施例である培養槽の運転制御装置の構成図である。
【００１５】
培養槽１内には、上下に翼を有する攪拌翼３が設置され、攪拌翼３は培養槽１の上部に設けられた駆動装置１０により回転駆動される。攪拌槽１の上部にはバルブ１１ａを有する通気管２ａが、下部にはバルブ１１ｂを有する通気管２ｂが設けられている。通気管
２ａからは液面に向けて酸素ガスが通気され、通気管２ｂからは液中に酸素ガスが通気されるようになっている。培養槽１の上部には、排気管１３が設けられ、酸素ガスや二酸化炭素ガスが排気される。培養槽１内には、細胞を培養するための培地と動物細胞あるいは微生物が封入されており、培養槽１には、培地に糖やアミノ酸などの栄養分を供給するための栄養成分供給装置５が設けられている。
【００１６】
排気管１３にはガスクロマトグラフ１４が接続されており、排気管１３から排気される酸素ガス，二酸化炭素ガス等の分析を行うようになっている。培養槽１には培地の状態を検出するためのサンプラー１５が接続され、培地の一部を採取し、サンプラー１５に接続された測定装置４により、培地の状態等を検出するようになっている。測定装置４により、グルコース濃度，グルタミン濃度など糖類あるいはアミノ酸類の栄養成分濃度，乳酸濃度，アンモニア濃度，たんぱく質濃度などの生産物濃度、および酸素濃度，細胞数密度などが計測される。
【００１７】
ガスクロマトグラフ１４及び測定装置４は、演算処理装置７に接続され、演算処理装置７は記憶装置８，出力表示装置９に接続されている。制御装置６には測定装置４が接続され、測定値がフィードバックされるようになっている。演算処理装置７で演算された制御目標値は、制御装置６に入力され、フィードバックされた測定値に基づいてＰＩＤ制御信号を駆動装置１０、バルブ１１ａ，１１ｂ、栄養成分供給装置５等に出力する。出力表示装置９では、後述するように、栄養成分濃度，生産物濃度，細胞数密度の時間変化率と、細胞内反応速度の組からなる数値ベクトルを、例えば図２に示すように画面２０で表示し、また、栄養成分濃度，生産物濃度，細胞数密度の予測値と実測値を、例えば図３に示すように画面２０で表示する。
【００１８】
記憶装置８には、演算処理装置７で演算処理を行うためのプログラムが記憶されており、演算処理装置７では、プログラムに従って、次のような処理を行う。
【００１９】
測定装置４によって、グルコース濃度，グルタミン濃度など糖類あるいはアミノ酸類の栄養成分濃度，乳酸濃度，アンモニア濃度，たんぱく質濃度などの生産物濃度、および酸素濃度，細胞数密度などが計測される。第１段階では、２つの時刻間の測定濃度をその時間間隔Δｔと細胞数密度Ｘa で割って、単位時間当たり単位細胞数当りの消費速度、あるいは生成速度を計算する。
【００２０】
例えば、グルタミン濃度Ｇｌｎについては、数（１）で計算される。ここで、添え字
ｏｂｓは観測値を表す。又、Ｇｌｎはグルタミン濃度で単位が（mg／Ｌ）であり、Ｘa は細胞数密度で単位は（cells／ｍＬ）である。
【００２１】
【数１】
2007202500.tif 000003
【００２２】
又、細胞数密度については、数（２）で計算され、μは比増殖速度と呼ばれ、単位は
（１／ｓ）である。
【００２３】
【数２】
2007202500.tif 000004
【００２４】
細胞の内部状態を完全に同定することは困難であるが、細胞内部で起こっていることを化学反応のレベルで捉えると、基質成分を取り込み、多数の中間代謝物を経て生産物を外部へ分泌する、物質変換を伴う複数の化学反応の集まりである。細胞内で起こるすべての反応を拾い上げることは困難であるが、比較的少数の主要な反応を集めて細胞内化学反応系の数学モデルを作ることは可能である。細胞内の反応速度を細胞の内部状態を記述する内部状態変数とし、細胞内反応系を記述する数学的モデルを組み立て、これを解くことにより細胞内反応速度を求める。
【００２５】
記憶装置８には、細胞内反応速度の回路網が記憶されており、図２に示される回路網の各々の物質名とＦ１からＦ３０の矢印の繋がりが記憶されている。図２において、
Gultamine，Lactateなどは、細胞外の培地中の栄養成分，生産物濃度であり、Pyruvate，Oxaloacetateなどは細胞内の中間代謝物質濃度である。矢印Ｆ１からＦ３０は、それぞれ１つの化学反応を表わしている。例えば、Ｆ１はGlucoseからPyruvate に変換される化学反応であることを示し、その反応速度がＦ１で与えられることを示す。
【００２６】
次に、細胞内反応速度の回路網から導かれる方程式について説明する。細胞を記述する内部状態変数としてこれらの反応の反応速度を一列に並べたベクトルを考えることができる。とくに培養環境が劇的に変化しない条件のもとでは、細胞内の代謝反応は概ね定常状態にあるので、おのおのの中間代謝物の濃度は一定であり、おのおのの中間代謝物の生成と消費に係わる反応速度の総和は０となる。この条件と基質の流入速度、および生産物の流出速度を用いると細胞内化学反応系の数学モデルに含まれるすべての反応速度を求めることができる。
【００２７】
例えば、図２に示す中間代謝物質Pryuvate濃度は、定常状態の仮定のもとでは一定であるから、Pyruvateに入り、あるいはPryuvateから出ていく反応速度の総和は０であり、図２から分かるように数（３）が成り立つ。
【００２８】
【数３】
2007202500.tif 000005
【００２９】
このようにして、図２に示す回路網に含まれる中間代謝物質の個数分Ｍ個の方程式が導かれる。又、例えば、図２においてGlutamineからGlutamateへ変換される反応Ｆ６の反応速度は、Glutamine濃度の計測値を用いて数（１）で計算されるから、数（４）が成り立つ。
【００３０】
【数４】
2007202500.tif 000006
【００３１】
このように、培地中の成分濃度の時間変化率は方程式系の入力データとして扱われ、この関係式は、測定されている培地中の成分濃度の個数分Ｎ個与えられる。以上のようにして導かれた方程式は、総数がＭ＋Ｎ個の連立１次方程式である。
【００３２】
連立１次方程式の総数Ｍ＋Ｎが、モデルで考慮している反応の数、例えば図２に示す例では３０個と等しければ、連立１次方程式は解くことができ、すべての反応速度が求まるので、図２の対応する矢印上に値を表示する。又、求められた反応速度を図３に示すように時間軸上に表示する。このように、すべての反応速度が求まるので、どのように制御するとどのような結果が得られるかを計算することができ、目標を選定して細胞内反応速度も考慮した制御が行える。
【００３３】
連立１次方程式の総数Ｍ＋Ｎが、モデルで考慮している反応速度の数を超えていれば、測定された培地中の成分濃度の時間変化率のいくつかは入力データとして用いる必要がなく、検証データとして用いることができる。
【００３４】
又、連立１次方程式の総数Ｍ＋Ｎが、モデルで考慮している反応速度の数よりも少ないときは、一意に解を求めることができないが、この場合はモデルをさらに簡略化して考慮する反応速度の数を減らすか、あるいは、例えばアンモニアの生成反応速度Ｆ３０を最小化するといった目的関数を定めて、線形計画問題に変更することにより、目標条件を満足する反応速度の組を求めることができる。
【００３５】
次に、上述したようにして求められる細胞内反応速度を培地中の栄養成分濃度，生産物濃度の関数として表現する。例えば、GlucoseからPyruvateへの変換を与える反応速度
Ｆ１を培地中のグルコース濃度(Ｇｌｃ)，グルタミン濃度(ｇｌｎ)，乳酸濃度(Ｌａｃ)，アンモニア濃度（Ａｍｍ）の関数として、Ｆ１＝Ｆ１(Ｇｌｃ，Ｇｌｎ，Ｌａｃ，Ａｍｍ)と表す。この関数が陽に求まれば、制御的に望ましい細胞内反応速度をＧｌｃ，Ｇｌｎ，Ｌａｃ，Ａｍｍなどの培地中の成分濃度を調節することにより制御することができる。
【００３６】
上述したように、細胞内反応速度は、栄養成分濃度，生産物濃度，細胞数密度の時間変化率を入力データとして回路網の方程式により求めることができるから、栄養成分濃度，生産物濃度，細胞数密度の時間変化率が、栄養成分濃度，生産物濃度の関数として与えられれば関数Ｆ１が陽に求まることになる。そこで、これらの関数を測定装置から得られる実測値を用いて推定する。
【００３７】
すなわち、数（５）から数（１５）に示すような、栄養成分濃度，生産物濃度，細胞数密度の時間変化率のモデル推定式を定める。
【００３８】
【数５】
2007202500.tif 000007
【００３９】
【数６】
2007202500.tif 000008
【００４０】
【数７】
2007202500.tif 000009
【００４１】
【数８】
2007202500.tif 000010
【００４２】
【数９】
2007202500.tif 000011
【００４３】
【数１０】
2007202500.tif 000012
【００４４】
【数１１】
2007202500.tif 000013
【００４５】
【数１２】
2007202500.tif 000014
【００４６】
【数１３】
2007202500.tif 000015
【００４７】
【数１４】
2007202500.tif 000016
【００４８】
【数１５】
2007202500.tif 000017
【００４９】
ここで、Ｋd は、死滅速度で単位は（１／ｓ）、Ａｍｍは、アンモニア濃度で単位は
（mg／Ｌ）、Ｌａｃは、乳酸濃度で単位は（mg／Ｌ）、Ｇｌｃは、グルコース濃度で単位は（mg／Ｌ）である。又、μmax，κd0，α，Ｋi，Ｋji，Ｙi，ｍi，ｑEij，κj(ｉ，ｊ＝Ｇｌｃ，Ｇｌｎ，Ｌａｃ，Ａｍｍ）はモデル推定式中のモデルパラメータである。栄養成分濃度，生産物濃度，細胞数密度の実測値と上述したアルゴリズムで計算される時間変化率を用いて、モデルパラメータを最小２乗法を用いてフィッティングして最良推定式を定める。
【００５０】
この最良推定式を入力データとして回路網方程式を解くことにより、任意の細胞内反応速度ＦｉをＦｉ＝Ｆｉ（Ｇｌｃ，Ｇｌｎ，Ｌａｃ，Ａｍｍ）のように栄養成分濃度，生産物濃度の関数として求めることができる。
【００５１】
求めたモデル推定式により培地中成分濃度の予測ができる。この結果、栄養成分濃度，生産物濃度，細胞数密度の時間変化率のモデル推定式を図４に示すようなアルゴリズムにより数値積分することにより、時間的未来の栄養成分濃度，生産物濃度，細胞数密度を計算でき、図３に示すように表示することができる。
【００５２】
数値積分では、図４に示すように、ステップ３１で、積分する関数の初期値を設定し、ステップ３２でタイムステップを１増加して更新し、ステップ３３で関数の差分を計算する。ステップ３４で、タイムステップ更新した回数ｎがｎｍａｘになったか否かを判断し、ｎｍａｘになっていなければ、ステップ３２に戻り、ステップ３３，３４を繰り返す。タイムステップ更新回数がｎｍａｘとなれば、ステップ３５で出力する。
【００５３】
このような計算結果を用いて、培養槽の予測制御を行う。例えば、ある生産物濃度を最大にするように制御を行う場合には、その生産物濃度を最大化するという目的関数を定めて、目標条件を満足する反応速度の組を求めることができる。ここで、細胞数密度を目的関数とすることもできる。上述した最良推定式を入力データとして回路網方程式を解くことにより制御すべき栄養成分濃度，生産物濃度等の動的な制御目標値を求めることができるので、サンプラー１５により計測された培地の状態、測定装置４で測定された生産物濃度，酸素濃度，細胞数密度などの測定値に基づいて、制御装置６により栄養成分供給装置５，バルブ１１，駆動装置１０を制御して、糖やアミノ酸などの栄養分の供給制御、酸素ガスや二酸化炭素濃度，攪拌制御が行われる。これらの制御結果は、図２に示す画面、図３に示す画面で表示され、視覚的に制御が正常に行われているかを通知する。
【００５４】
本実施例によれば、栄養成分濃度，生産物濃度，細胞数密度などの培地中観測変数と細胞内反応速度を培養槽の状態変数として用いて、好適な培養槽の運転制御ができる。
【図面の簡単な説明】
【００５５】
【図１】本発明の一実施例である培養槽の運転制御装置の構成図である。
【図２】本実施例の出力表示装置の表示画面の一例を示す平面図である。
【図３】本実施例の出力表示装置の表示画面の一例を示す平面図である。
【図４】モデル推定式の時間積分アルゴリズムを示す図である。
【符号の説明】
【００５６】
１…培養槽、２…通気装置、３…攪拌翼、４…測定装置、５…栄養成分供給装置、６…制御装置、７…演算処理装置、８…記憶装置、９…出力表示装置。
【図１】
2007202500.tif 000018
【図２】
2007202500.tif 000019
【図３】
2007202500.tif 000020
【図４】
2007202500.tif 000021