JPWO2010001584A1 - インスリン抵抗性評価支援システム、インスリン抵抗性評価支援方法、及びコンピュータプログラム - Google Patents

Abstract

【課題】 インスリン抵抗性に関する生体情報を推定するために必要となる被験者の負担を、従来に比して少なくすることが可能なインスリン抵抗性評価支援システム、インスリン抵抗性評価支援方法、及びコンピュータプログラムを提供する。【解決手段】インスリン抵抗性評価支援システム１は、被験者を測定することにより得られた血中のグルコース濃度、インスリン濃度、及び遊離脂肪酸濃度に関する情報の入力を受け付ける入出力インタフェース１１ｆと、入力を受け付けた血中のグルコース濃度、インスリン濃度、及び遊離脂肪酸濃度に関する情報に基づいて、被験者のグルコース取込速度を推定するＣＰＵ１１ａと、推定されたグルコース取込速度を出力する画像出力インタフェース１１ｇとを備える。【選択図】 図１

Description

本発明は、被験者のインスリン抵抗性の評価を支援するインスリン抵抗性評価支援システム、その方法、及びコンピュータをインスリン抵抗性評価支援システムとして機能させるコンピュータプログラムに関する。

インスリン抵抗性は、糖尿病の病態の一つであり、またメタボリックシンドロームの重要な背景因子として知られている。従来、インスリン抵抗性の評価には、グルコースクランプ試験が用いられている。このグルコースクランプ試験は、人工膵臓を用いてインスリンを被験者の静脈に注射し、血糖値を一定に保つようにグルコースの注入速度を調節するものであり、被験者への侵襲度が高く、被験者の負担が大きい検査方法である。

特許文献１には、「空腹時インスリン値」、「食後２時間血糖値」、「ＨＯＭＡ−ＩＲ」及び「インスリンＯＧＴＴ頂値」を入力値に含み、これらの入力値のそれぞれについて、予め与えられた判定基準値と比較してスコア値を求める末梢のインスリン抵抗性判定処理を実行することにより、糖尿病及びメタボリックシンドロームの疾患リスクを分析する診断支援システムが開示されている。

特開２００６−３０４８３３号公報

しかしながら、特許文献１に記載の診断支援システムでは、経口糖負荷試験（ＯＧＴＴ：Oral Glucose Tolerance Test）の検査結果を入力情報として必要とする。ＯＧＴＴ は、口からブドウ糖を摂取して、所定時間経過後に数回採血して血糖値と血中インスリン濃度を測定する試験であり、グルコースクランプに比べると、被験者への負担が少ないが、試験時間として数時間を必要とする。

本発明は斯かる事情に鑑みてなされたものであり、その主たる目的は、上記課題を解決することができるインスリン抵抗性評価支援システム、インスリン抵抗性評価支援方法及びコンピュータプログラムを提供することにある。

上述した課題を解決するために、本発明の一の態様のインスリン抵抗性評価支援システムは、被験者のインスリン抵抗性の評価を支援するためのインスリン抵抗性評価支援システムであって、被験者を測定することにより得られた血中のグルコース濃度、インスリン濃度、及び遊離脂肪酸濃度に関する情報の入力を受け付ける入力手段と、前記入力手段によって入力を受け付けた血中のグルコース濃度、インスリン濃度、及び遊離脂肪酸濃度に関する情報に基づいて、被験者のグルコース取込速度を推定する推定手段と、を備える。

この態様において、前記推定手段は、血中のインスリン濃度及び遊離脂肪酸濃度に関する情報に基づいて、グルコーストランスポータの発現量に応じた数値を推定する第１推定手段と、前記第１推定手段によって推定されたグルコーストランスポータの発現量に応じた数値と血中のグルコース濃度とに基づいて、グルコース取込速度を推定する第２推定手段と、を備えることが好ましい。

また、この場合においては、前記第１推定手段が、血中の遊離脂肪酸濃度に関する情報に基づいて、脂肪酸アシル−補酵素Ａ複合体の濃度を推定する第３推定手段と、血中のインスリン濃度及び前記第１推定手段によって推定された脂肪酸アシル−補酵素Ａ複合体の濃度に基づいて、グルコーストランスポータの発現量に応じた数値を推定する第４推定手段と、を備える構成とすることが好ましい。

また、この場合においては、前記推定手段が、前記第２推定手段によって推定されたグルコース取込速度に基づいて、細胞内のピルビン酸の濃度を推定する第５推定手段と、前記第３推定手段によって推定された脂肪酸アシル−補酵素Ａ複合体の濃度と、前記第５推定手段によって推定されたピルビン酸の濃度とに基づいて、細胞内のアセチル補酵素の濃度を推定する第６推定手段と、前記第６推定手段によって推定されたアセチル補酵素の濃度に基づいて、脂肪酸アシル−補酵素Ａ複合体から生成されるアセチル補酵素の生成速度を調整する調整手段と、を更に備え、前記第６推定手段は、前記調整手段による調整後のアセチル補酵素の生成速度に基づいて、アセチル補酵素の濃度を再度推定するように構成されていることが好ましい。

また、上記態様においては、前記入力手段が、被験者の筋肉量に関する情報の入力をさらに受け付けるように構成されており、前記第２推定手段が、前記第１推定手段によって推定されたグルコーストランスポータの発現量に応じた数値と血中のグルコース濃度とに基づいて、被験者の単位筋肉量当たりのグルコース取込速度を推定するように構成されており、前記推定手段が、前記第２推定手段によって推定された被験者の単位筋肉量当たりのグルコース取込速度と、前記入力手段によって入力を受け付けた被験者の筋肉量に関する情報とに基づいて、被験者の単位重量当たりのグルコース取込速度を推定するように構成されていることが好ましい。

また、上記態様においては、前記推定手段が、生体内におけるグルコースからピルビン酸への変化に関連する複数の物質の生成反応において、生成反応前の各物質の濃度に基づいて、各物質の一定時間あたりの生成量を取得し、それぞれの物質について、前記生成反応前の物質の濃度にその物質の前記生成量を反映させることにより、前記一定時間後の各物質の濃度を取得し、取得された各物質の濃度に基づいて、グルコース取込速度を推定する処理を、繰り返し実行するように構成されていることが好ましい。

また、この場合においては、前記推定手段が、グルコース取込速度が定常状態に至ったか否かを判定し、グルコース取込速度が定常状態に至ったと判定されるまで、前記各物質の濃度を取得する処理を繰り返し実行するように構成されていることが好ましい。

また、上記態様においては、前記入力手段が、被験者の空腹時における血中のグルコース濃度、インスリン濃度、及び遊離脂肪酸濃度に関する情報の入力を受け付けるように構成されており、前記推定手段が、前記入力手段によって入力を受け付けた血中のグルコース濃度、インスリン濃度、及び遊離脂肪酸濃度に関する情報に基づいて、空腹時における各物質の濃度を取得し、取得された各物質の濃度に基づいて、空腹時におけるグルコース取込速度を推定する第１グルコース取込速度推定手段と、前記第１グルコース取込速度推定手段によって取得された空腹時における各物質の濃度に基づいて、インスリン濃度が所定値の場合における各物質の濃度を取得し、取得された各物質の濃度に基づいて、インスリン濃度が所定値の場合におけるグルコース取込速度を推定する第２グルコース取込速度推定手段と、を備えることが好ましい。

また、上記態様においては、前記推定手段によって推定されたグルコース取込速度を出力する出力手段をさらに備える構成とするか、又は、前記推定手段によって推定されたグルコース取込速度に関する情報に基づいて、被験者のインスリン抵抗性を推定するインスリン抵抗性推定手段と、前記インスリン抵抗性推定手段による推定結果を出力する出力手段と、をさらに備える構成とすることが好ましい。

本発明の一の態様のインスリン抵抗性評価支援方法は、入力装置を備えるコンピュータにより、被験者のインスリン抵抗性の評価を支援するためのインスリン抵抗性評価支援方法であって、被験者を測定することにより得られた血中のグルコース濃度、インスリン濃度、及び遊離脂肪酸濃度に関する情報の入力を前記入力装置により受け付けるステップと、前記入力装置によって入力を受け付けた血中のグルコース濃度、インスリン濃度、及び遊離脂肪酸濃度に関する情報に基づいて、前記コンピュータにより被験者のグルコース取込速度を推定するステップと、を有する。

この場合においては、前記被験者のグルコース取込速度を推定するステップが、血中のインスリン濃度及び遊離脂肪酸濃度に関する情報に基づいて、前記演算部によりグルコーストランスポータの発現量に応じた数値を推定するステップと、推定されたグルコーストランスポータの発現量に応じた数値と血中のグルコース濃度とに基づいて、前記コンピュータによりグルコース取込速度を推定するステップと、を含むことが好ましい。

本発明の一の態様のコンピュータプログラムは、入力装置を備えるコンピュータに、被験者を測定することにより得られた血中のグルコース濃度、インスリン濃度、及び遊離脂肪酸濃度に関する情報の入力を前記入力装置により受け付けるステップと、前記入力手段によって入力を受け付けた血中のグルコース濃度、インスリン濃度、及び遊離脂肪酸濃度、並びに被験者の血中のグルコース濃度に基づいて、被験者のグルコース取込速度を推定するステップと、を実行させる。

この場合においては、前記被験者のグルコース取込速度を推定するステップにおいて、前記コンピュータに、血中のインスリン濃度及び遊離脂肪酸濃度に関する情報に基づいて、グルコーストランスポータの発現量に応じた数値を推定するステップと、推定されたグルコーストランスポータの発現量に応じた数値と血中のグルコース濃度とに基づいて、グルコース取込速度を推定するステップと、を実行させることが好ましい。

本発明に係るインスリン抵抗性評価支援システム、インスリン抵抗性評価支援方法、及びコンピュータプログラムによれば、インスリン抵抗性に関する生体情報を推定するために必要となる被験者の負担を、従来に比して少なくすることが可能となる。

実施の形態１に係るインスリン抵抗性評価支援システムの構成を示すブロック図。 実施の形態１に係るインスリン抵抗性評価支援システムの機能的な構成を示す概念図。 実施の形態１に係るインスリン抵抗性評価支援システムにおける仮想的な物質の反応の流れを示す概念図。 実施の形態１に係るインスリン抵抗性評価支援プログラムの処理の流れを示すフローチャート。 実施の形態１に係るインスリン抵抗性評価支援システムの入力画面の一例を示す模式図。 実施の形態１に係るインスリン抵抗性評価支援システムの第１グルコース取込速度推定処理の手順を示すフローチャート。 実施の形態１に係るインスリン抵抗性評価支援システムの第２グルコース取込速度推定処理の手順を示すフローチャート。 実施の形態１に係るインスリン抵抗性評価支援システムの出力画面の一例を示す模式図。 実施の形態２に係るインスリン抵抗性評価支援システムの構成を示すブロック図。 実施の形態２に係るインスリン抵抗性評価支援システムにおける仮想的な物質の反応の流れを示す概念図。 実施の形態２に係るインスリン抵抗性評価支援プログラムの処理の流れを示すフローチャート。 実施の形態２に係るインスリン抵抗性評価支援システムの入力画面の一例を示す模式図。 実施の形態２に係るインスリン抵抗性評価支援システムの第１グルコース取込速度推定処理の手順を示すフローチャート。 実施の形態２に係るインスリン抵抗性評価支援システムの第２グルコース取込速度推定処理の手順を示すフローチャート。 実施の形態２に係るインスリン抵抗性評価支援システムの評価実験におけるシミュレーション結果表示画面を示す図。 実施の形態２に係るインスリン抵抗性評価支援システムの評価実験における血中グルコース濃度のシミュレーション結果を示すグラフ。 実施の形態２に係るインスリン抵抗性評価支援システムの評価実験における血中インスリン濃度のシミュレーション結果を示すグラフ。 実施の形態２に係るインスリン抵抗性評価支援システムの評価実験における血中脂肪酸濃度のシミュレーション結果を示すグラフ。 実施の形態２に係るインスリン抵抗性評価支援システムの評価実験におけるグルコース取込速度のシミュレーション結果を示すグラフ。 文献中の実測値とシミュレーションの結果との比較結果を示すグラフ。

以下、本発明の好ましい実施の形態を、図面を参照しながら説明する。

（実施の形態１）
本実施の形態は、被験者を測定することにより得られた血中のグルコース濃度、インスリン濃度、遊離脂肪酸濃度に関する情報、並びに骨格筋の単位時間あたりの酸素消費量及び二酸化炭素産生量の入力を受け付け、入力を受け付けた血中のグルコース濃度、インスリン濃度、遊離脂肪酸濃度に関する情報、骨格筋の単位時間あたりの酸素消費量、骨格筋の単位時間あたりの二酸化炭素産生量に基づいて、被験者のグルコース取込速度を推定し、推定されたグルコース取込速度に基づいて、被験者のインスリン抵抗性を推定し、推定結果を出力するインスリン抵抗性評価支援システムである。

［インスリン抵抗性評価支援システムの構成］
図１は、本実施の形態に係るインスリン抵抗性評価支援システムの構成を示すブロック図である。本実施の形態に係るインスリン抵抗性評価支援システム１は、コンピュータ１ａによって実現される。図１に示すように、コンピュータ１ａは、本体１１と、画像表示部１２と、入力部１３とを備えている。本体１１は、ＣＰＵ１１ａと、ＲＯＭ１１ｂ、ＲＡＭ１１ｃ、ハードディスク１１ｄ、読出装置１１ｅ、入出力インタフェース１１ｆ、及び画像出力インタフェース１１ｇを備えており、ＣＰＵ１１ａ、ＲＯＭ１１ｂ、ＲＡＭ１１ｃ、ハードディスク１１ｄ、読出装置１１ｅ、入出力インタフェース１１ｆ、および画像出力インタフェース１１ｇは、バス１１ｉによって接続されている。

ＣＰＵ１１ａは、ＲＡＭ１１ｃにロードされたコンピュータプログラムを実行することが可能である。そして、後述するようなインスリン抵抗性評価支援プログラム１４ａを当該ＣＰＵ１１ａが実行することにより、コンピュータ１ａがインスリン抵抗性評価支援システム１として機能する。

ＲＯＭ１１ｂは、マスクＲＯＭ、ＰＲＯＭ、ＥＰＲＯＭ、又はＥＥＰＲＯＭ等によって構成されており、ＣＰＵ１１ａに実行されるコンピュータプログラムおよびこれに用いるデータ等が記録されている。

ＲＡＭ１１ｃは、ＳＲＡＭまたはＤＲＡＭ等によって構成されている。ＲＡＭ１１ｃは、ハードディスク１１ｄに記録されているインスリン抵抗性評価支援プログラム１４ａの読み出しに用いられる。また、ＣＰＵ１１ａがコンピュータプログラムを実行するときに、ＣＰＵ１１ａの作業領域として利用される。

ハードディスク１１ｄは、オペレーティングシステムおよびアプリケーションプログラム等、ＣＰＵ１１ａに実行させるための種々のコンピュータプログラムおよび当該コンピュータプログラムの実行に用いられるデータがインストールされている。後述するインスリン抵抗性評価支援プログラム１４ａも、このハードディスク１１ｄにインストールされている。

読出装置１１ｅは、フレキシブルディスクドライブ、ＣＤ−ＲＯＭドライブ、またはＤＶＤ−ＲＯＭドライブ等によって構成されており、可搬型記録媒体１４に記録されたコンピュータプログラムまたはデータを読み出すことができる。また、可搬型記録媒体１４には、コンピュータをインスリン抵抗性評価支援システムとして機能させるためのインスリン抵抗性評価支援プログラム１４ａが格納されており、コンピュータ１ａが当該可搬型記録媒体１４からインスリン抵抗性評価支援プログラム１４ａを読み出し、当該インスリン抵抗性評価支援プログラム１４ａをハードディスク１１ｄにインストールすることが可能である。

なお、前記インスリン抵抗性評価支援プログラム１４ａは、可搬型記録媒体１４によって提供されるのみならず、電気通信回線（有線、無線を問わない）によってコンピュータ１ａと通信可能に接続された外部の機器から前記電気通信回線を通じて提供することも可能である。例えば、前記インスリン抵抗性評価支援プログラム１４ａがインターネット上のサーバコンピュータのハードディスク内に格納されており、このサーバコンピュータにコンピュータ１ａがアクセスして、当該コンピュータプログラムをダウンロードし、これをハードディスク１１ｄにインストールすることも可能である。

また、ハードディスク１１ｄには、例えば米マイクロソフト社が製造販売するWindows（登録商標）等のマルチタスクオペレーティングシステムがインストールされている。以下の説明においては、本実施の形態に係るインスリン抵抗性評価支援プログラム１４ａは当該オペレーティングシステム上で動作するものとしている。なお、インスリン抵抗性評価支援プログラム１４ａの構成の詳細については後述する。

入出力インタフェース１１ｆは、例えばUSB，IEEE1394，又はRS-232C等のシリアルインタフェース、SCSI，IDE，又はIEEE1284等のパラレルインタフェース、およびＤ／Ａ変換器、Ａ／Ｄ変換器等からなるアナログインタフェース等から構成されている。入出力インタフェース１１ｆには、キーボードおよびマウスからなる入力部１３が接続されており、ユーザが当該入力部１３を使用することにより、コンピュータ１ａにデータを入力することが可能である。

画像出力インタフェース１１ｇは、ＬＣＤまたはＣＲＴ等で構成された画像表示部１２に接続されており、ＣＰＵ１１ａから与えられた画像データに応じた映像信号を画像表示部１２に出力するようになっている。画像表示部１２は、入力された映像信号にしたがって、画像（画面）を表示する。

［インスリン抵抗性評価支援システムの機能的な構成］
次に、インスリン抵抗性評価支援プログラム１４ａについて更に詳細に説明する。図２は、インスリン抵抗性評価支援プログラム１４ａにより実現されるインスリン抵抗性評価支援システムの機能的な構成を示す概念図であり、図３は、当該インスリン抵抗性評価支援システムにおける仮想的な物質の反応の流れを示す概念図である。インスリン抵抗性評価支援プログラム１４ａは、被験者の体重、骨格筋の単位時間あたりの酸素消費量、二酸化炭素生成量、骨格筋率、空腹時の血中グルコース濃度、血漿インスリン濃度、及び血中遊離脂肪酸濃度の各生体情報を入力とし、被験者の末梢組織（骨格筋）の糖取込速度（糖取込率）の推定値を出力とする。かかるインスリン抵抗性評価支援プログラム１４ａは、身体の糖取込に関連する生体器官（機能）を機能別に分けたそれぞれの機能を仮想的に再現する４つのブロック２１〜２４を含んでいる。それぞれのブロック２１〜２４は、複数のパラメータを含んでおり、グルコースの取込に関連する物質の生成の反応速度を演算するように構成されている。また、インスリン抵抗性評価支援プログラム１４ａは、上記の反応速度から、グルコースの取込に関連する物質の生成速度（一定時間あたりの生成量）を演算し、この演算によって得られた一定時間後の生成量を反映した、一定時間後の物質濃度を演算するように構成されている。

以下に、各ブロックの詳細な構成について説明するが、その説明に先立って、各ブロックにおける演算の基礎となる一般的な代謝反応速度の概念について説明する。生体器官における一般的に不可逆な酵素反応は次式のように表される。
Ｘ＋Ｙ＋Ｅ_１→Ｖ＋Ｗ＋Ｅ_２ …（１）
ここで、Ｘ，Ｙは代謝される基質濃度を、Ｖ，Ｗは生成される基質濃度を示しており、Ｅ_１ ，Ｅ_２ はＡＴＰ，ＡＤＰ若しくはＡＤＰ，ＡＴＰ、及び／又は、ＮＡＤＨ，ＮＡＤ又はＮＡＤ，ＮＡＤＨを示している。

式（１）によって表される酵素反応は、次の式（２）のように表すこともできる。

上記の酵素反応の反応速度ｆ_{Ｘ＋Ｙ→Ｖ＋Ｗ}は、以下の式（３）で求められる。

ただし、Ｖ_{Ｘ＋Ｙ→Ｖ＋Ｗ}は飽和最大速度を、Ｃ_Ｘ ，Ｃ_Ｙ は基質Ｘ，Ｙの濃度を、ＰＳ^＋はＡＴＰ／ＡＤＰを、ＰＳ⁻はＡＤＰ／ＡＴＰを、ＲＳ^＋はＮＡＤＨ／ＮＡＤを、ＲＳ⁻はＮＡＤ／ＮＡＤＨをそれぞれ示しており、Ｋ_Ｘ ，Ｋ_Ｙ ，μ^±，ν^±は反応過程におけるミカエリスメンテン定数又は代謝制御に関するモデルパラメータを示している。

なお、以下の各ブロック２１〜２４及び補足演算処理の説明における物質の反応速度、生成速度、及び濃度は、特に断りのない限り単位重量当たりの反応速度（ｍｍｏｌ／ｋｇ／ｍｉｎ）、生成速度（ｍｍｏｌ／ｋｇ／ｍｉｎ）、及び濃度（ｍＭ）とする。

＜脂肪酸代謝ブロック＞
脂肪酸代謝ブロック２１は、生体器官の脂肪酸代謝機能を仮想的に再現する機能ブロックである。身体の脂肪酸代謝機能は、血中の遊離脂肪酸（ＦＦＡ）を細胞内に取り込んでジアシルグリセロール（ＤＡＧ）を介して中性脂肪（ＴＧ）を生成し、これとともに、遊離脂肪酸から脂肪酸アシル−補酵素Ａ複合体（ＦＡＣ）を生成する。脂肪酸代謝ブロック２１は、かかる身体の脂肪酸代謝機能を模したものであり、当該脂肪酸代謝ブロック２１をＣＰＵ１１ａが実行することにより、細胞内遊離脂肪酸濃度に基づいて、遊離脂肪酸が中性脂肪に変換される反応速度ｆ_{ＦＦＡ→ＴＧ}、遊離脂肪酸から脂肪酸アシル−補酵素Ａ複合体が生成される反応速度ｆ_{ＦＦＡ→ＦＡＣ}、中性脂肪から遊離脂肪酸が生成される反応速度ｆ_{ＴＧ→ＦＦＡ}が演算され、また、一定時間経過後の細胞内の遊離脂肪酸濃度、中性脂肪濃度、及び脂肪酸アシル−補酵素Ａ複合体濃度が演算される。

以下に、生体器官の脂肪酸代謝機能に関連する化学反応と、それに基づいた脂肪酸代謝ブロック２１の具体的な演算処理について説明する。なお、以下の式における遊離脂肪酸、中性脂肪、ＡＴＰ、ＡＤＰ、脂肪酸アシル−補酵素Ａ複合体、及びリン酸の細胞内濃度は、それぞれ所定の初期値が与えられており、１度目の演算では初期値が用いられ、その後の演算では更新された値が用いられる。まず、血中の遊離脂肪酸は細胞内に取り込まれるが、このときの血液から組織（細胞）への受動的な遊離脂肪酸の流入速度ｆ_ＦＦＡは、以下の式（４）で表される。

ただし、ｆ_Ｏ２は組織内での酸素の消費速度を示しており、入力された骨格筋の単位時間あたりの酸素消費量と骨格筋率とによって求められる。また、ｆ_ＣＯ２は組織内での二酸化炭素の生成速度を示しており、入力された骨格筋の単位時間あたりの二酸化炭素産生量と骨格筋率とによって求められる。

上記の遊離脂肪酸の流入速度（組織への取込速度）ｆ_ＦＦＡは、次式（５）によっても表すことができる。

ただし、Ｃ_ＦＦＡｂは血中遊離脂肪酸濃度を、Ｃ_ＦＦＡは細胞内の遊離脂肪酸濃度を、λ_ＦＦＡは膜透過性を含むＦＦＡの膜輸送の係数を、σ_ＦＦＡはＦＦＡに関する分配係数をそれぞれ示す。ＣＰＵ１１ａにより、上記の式（４）で求めたｆ_ＦＦＡと、入力された血中遊離脂肪酸濃度とを式（５）に適用することで、細胞内の遊離脂肪酸濃度ＦＦＡを演算する。

次に、細胞内に取り込まれた遊離脂肪酸（ＦＦＡ）は、ジアシルグリセロールを介して中性脂肪（ＴＧ）に変換される。
３ＦＦＡ＋６ＡＴＰ→ＴＧ＋６ＡＤＰ …（６）
この反応の速度は、次式（７）で示される。

また、遊離脂肪酸（ＦＦＡ）と補酵素Ａ（ＣｏＡ）から、脂肪酸アシル−補酵素Ａ複合体（ＦＡＣ）が生成される。
ＦＦＡ＋ＣｏＡ＋２ＡＴＰ→ＦＡＣ＋２ＡＤＰ＋２Ｐｉ …（８）
また、この反応の速度は、次式（９）で示される。

更に、中性脂肪は、次式（１０）の通り、分解されて遊離脂肪酸となる。
ＴＧ→３ＦＦＡ …（１０）
この反応の速度は、次式（１１）で示される。

上記の式（８）の反応により生成された脂肪酸アシル−補酵素Ａ複合体は、ミトコンドリアに与えられる。また、脂肪酸アシル−補酵素Ａ複合体の濃度情報は、下記のように更新される。

脂肪酸代謝ブロック２１では、ＣＰＵ１１ａにより、上記の式（７），（９），（１１）で表されるｆ_{ＦＦＡ→ＴＧ}，ｆ_{ＦＦＡ→ＦＡＣ}，ｆ_{ＴＧ→ＦＦＡ}の各反応速度が演算される。

また、脂肪酸代謝ブロック２１では、ＣＰＵ１１ａにより、下式（１２）で表される遊離脂肪酸の生成速度、下式（１３）で表される中性脂肪の生成速度、下式（１４）で表される脂肪酸アシル−補酵素Ａ複合体の生成速度がそれぞれ演算される。

ただし、式（１２）において、反応速度ｆ_{ＦＡＣ→ＡＣｏＡ}は、後述する式（３０）で表され、ミトコンドリアブロック２４において算出される。ミトコンドリアブロック２４の演算が一度も行われていないときには、反応速度ｆ_{ＦＡＣ→ＡＣｏＡ}の初期値が用いられる。

さらに脂肪酸代謝ブロック２１では、ＣＰＵ１１ａにより、上記のようにして得られた遊離脂肪酸、中性脂肪、及び脂肪酸アシル−補酵素Ａ複合体の生成速度から、一定時間に生成される遊離脂肪酸、中性脂肪、及び脂肪酸アシル−補酵素Ａ複合体の量がそれぞれ演算され、その時点での遊離脂肪酸濃度、中性脂肪濃度、脂肪酸アシル−補酵素Ａ複合体濃度にこれらの量を反映させることで、前記一定時間経過後の細胞内の遊離脂肪酸濃度、中性脂肪酸濃度、脂肪酸アシル−補酵素Ａ複合体濃度が演算される。なお、本実施の形態ではこの一定時間は定数とされるが、ユーザが設定可能とすることもできる。

また、脂肪酸代謝ブロック２１は、後述するように、第１グルコース取込速度推定処理と、第２グルコース取込速度推定処理の両方において用いられる。第１グルコース取込速度推定処理では、脂肪酸代謝ブロック２１を用いることにより、空腹時における細胞内の遊離脂肪酸濃度、中性脂肪酸濃度、脂肪酸アシル−補酵素Ａ複合体濃度を演算する処理が行われる。また、第２グルコース取込速度演算処理では、脂肪酸代謝ブロック２１を用いることにより、高インスリン血症時における前記各物質の細胞内濃度を演算する処理が行われる。

＜インスリンシグナルブロック＞
インスリンシグナルブロック２２は、生体器官のグルコーストランスポータ（ＧＬＵＴ４）の発現量調整機能を仮想的に再現する機能ブロックである。生体のＧＬＵＴ４発現量調整機能は、インスリン受容体に結合したインスリン量と、脂肪酸代謝の代謝生成物である脂肪酸アシル−補酵素Ａ複合体の濃度及びジアシルグリセロールの濃度に応じてＧＬＵＴ４の発現量を調整する。本実施の形態に係るインスリンシグナルブロック２２は、かかる身体のＧＬＵＴ４発現量調整機能を模したものであり、当該インスリンシグナルブロック２２をＣＰＵ１１ａが実行することにより、血漿インスリン濃度（ＰＩ）とＦＡＣ濃度とに基づいて、ＧＬＵＴ４の発現量に応じた数値が演算される。

以下に、生体器官のＧＬＵＴ４発現量調整機能に基づいたインスリンシグナルブロック２２の具体的な演算処理について説明する。生体器官のＧＬＵＴ４発現量調整機能は、血漿インスリン濃度の上昇に従い、ＧＬＵＴ４の発現量が増加し、ＦＡＣ濃度及びＤＡＧ濃度の上昇に従い、ＧＬＵＴ４の発現量が抑制されるという特徴を有している。かかる特徴を考慮し、インスリンシグナルブロック２２においては、下式（１５）に従って、入力された血漿インスリン濃度と、脂肪酸代謝ブロック２１の演算で得られたＦＡＣ濃度とを用いて、ＧＬＵＴ４の発現量に応じた数値であるグルコースの取込率（ＧＬＵＴ）が演算される。なお、本実施の形態においては、ＤＡＧ濃度を考慮せずにグルコースの取込率を演算することとしている。

ただし、Ｖ_ｍａｘは所定の係数、ｎＰＩ及びｎＦＡＣは定数である。なお、グルコース取込率ＧＬＵＴは、ＧＬＵＴ４の発現量に比例する実数である。

また、インスリンシグナルブロック２２は、後述するように、第１グルコース取込速度推定処理と、第２グルコース取込速度推定処理の両方において用いられる。第１グルコース取込速度推定処理では、インスリンシグナルブロック２２を用いることにより、空腹時におけるグルコースの取込率ＧＬＵＴを演算する処理が行われる。また、第２グルコース取込速度演算処理では、インスリンシグナルブロック２２を用いることにより、高インスリン血症時におけるグルコースの取込率ＧＬＵＴを演算する処理が行われる。

＜解糖系ブロック＞
解糖系ブロック２３は、生体器官のグルコース分解機能を仮想的に再現する機能ブロックである。身体のグルコース分解機能は、ＧＬＵＴ４の発現量に応じてグルコースを細胞内に取り込み、細胞内のグルコースを分解して、Ｇ６Ｐ（グルコース６リン酸）、ＧＡ３Ｐ（グリセルアルデヒド３リン酸）を介してピルビン酸を産生する。解糖系ブロック２３は、かかる身体のグルコース分解機能を模したものであり、当該解糖系ブロック２３をＣＰＵ１１ａが実行することにより、インスリンシグナルブロック２２で得られたグルコース取込率ＧＬＵＴ、組織内での酸素消費速度及び二酸化炭素産生速度に基づいてグルコース取込速度が演算され、グルコースがＧ６Ｐに変換される反応速度ｆ_{ＧＬＵ→Ｇ６Ｐ}、Ｇ６ＰがＧＡ３Ｐに変換される反応速度ｆ_{Ｇ６Ｐ→ＧＡ３Ｐ}、ＧＡ３Ｐがピルビン酸に変換される反応速度ｆ_{ＧＡ３Ｐ→ＰＹＲ}がそれぞれ演算される。また、これらの反応速度に基づき、一定時間経過後のグルコース濃度、Ｇ６Ｐ濃度、ＧＡ３Ｐ濃度、及びピルビン酸濃度が演算される。

以下に、生体器官のグルコース分解機能に関連する化学反応と、それに基づいた解糖系ブロック２３の具体的な演算処理について説明する。なお、以下の式におけるグルコース、Ｇ６Ｐ、ＧＡ３Ｐ、ピルビン酸、ＮＡＤ、及びＮＡＤＨの細胞内濃度は、それぞれ所定の初期値が与えられており、１度目の演算では初期値が用いられ、その後の演算では更新された値が用いられる。
まず、血中のグルコースが、ＧＬＵＴ４の発現量（細胞表面における発現量）に応じて細胞内に取り込まれる。この取り込み速度ｆ_ＧＬＵは、以下の式（１６）で表される。

ただし、Ｃ_ＧＬＵｂは、入力された血中のグルコース濃度を、Ｃ_ＧＬＵは細胞内のグルコース濃度を、σ_ＧＬＵはグルコースに関する分配係数をそれぞれ示す。

細胞内へ取り込まれたグルコース（ＧＬＵ）はリン酸化され、Ｇ６Ｐに変換される。
ＧＬＵ＋ＡＴＰ→Ｇ６Ｐ＋ＡＤＰ …（１７）
この反応の速度は、次式（１８）で表される。

また、変換されて生じたＧ６ＰはＡＴＰと反応して、次式の通りＧＡ３ＰとＡＤＰとを生じる。
Ｇ６Ｐ＋ＡＴＰ→２ＧＡ３Ｐ＋ＡＤＰ …（１９）
この反応の速度は、次式（２０）で示される。

変換されて生じたＧＡ３Ｐは、次式に従ってピルビン酸（ＰＹＲ）に変換される。
ＧＡ３Ｐ＋Ｐｉ＋ＮＡＤ＋２ＡＤＰ→ＰＹＲ＋ＮＡＤＨ＋２ＡＴＰ …（２１）
このＰＹＲに変換される反応速度は、次式（２２）で示される。

解糖系ブロック２３では、ＣＰＵ１１ａにより、上記の式（１８），（２０），（２２）で表されるｆ_{ＧＬＵ→Ｇ６Ｐ}，ｆ_{Ｇ６Ｐ→ＧＡ３Ｐ}，ｆ_{ＧＡ３Ｐ→ＰＹＲ}の各反応速度が演算される。

また、解糖系ブロック２３では、ＣＰＵ１１ａにより、下式（２３）で表されるグルコースの生成速度、下式（２４）で表されるＧ６Ｐの生成速度、下式（２５）で表されるＧＡ３Ｐの生成速度、下式（２６）で表されるピルビン酸の生成速度がそれぞれ演算される。

ただし、式（２６）において、反応速度ｆ_{ＰＹＲ→ＡＣｏＡ}は、後述する式（２８）で表され、ミトコンドリアブロック２４において算出される。ミトコンドリアブロック２４の演算が一度も行われていないときには、反応速度ｆ_{ＰＹＲ→ＡＣｏＡ}の初期値が用いられる。

さらに解糖系ブロック２３では、ＣＰＵ１１ａにより、上記のようにして得られたグルコース、Ｇ６Ｐ、ＧＡ３Ｐ、及びピルビン酸の生成速度から、一定時間に生成されるグルコース、Ｇ６Ｐ、ＧＡ３Ｐ、及びピルビン酸の量がそれぞれ演算され、その時点での細胞内におけるグルコース濃度、Ｇ６Ｐ濃度、ＧＡ３Ｐ濃度、ピルビン酸濃度にこれらの量を反映させることで、前記一定時間経過後の細胞内のグルコース濃度、Ｇ６Ｐ濃度、ＧＡ３Ｐ濃度、及びピルビン酸濃度が演算される。

また、解糖系ブロック２３は、後述するように、第１グルコース取込速度推定処理と、第２グルコース取込速度推定処理の両方において用いられる。第１グルコース取込速度推定処理では、解糖系ブロック２３を用いることにより、空腹時における細胞内のグルコース濃度、Ｇ６Ｐ濃度、ＧＡ３Ｐ濃度、及びピルビン酸濃度を演算する処理が行われる。また、第２グルコース取込速度演算処理では、解糖系ブロック２３を用いることにより、高インスリン血症時における前記各物質の細胞内濃度を演算する処理が行われる。

＜ミトコンドリアブロック＞
ミトコンドリアブロック２４は、生体器官のミトコンドリアの機能を仮想的に再現する機能ブロックである。ミトコンドリアは、ピルビン酸及び脂肪酸アシル−補酵素Ａ複合体をそれぞれ酸化させることでアセチル補酵素Ａ（ＡＣｏＡ）を生成し、このアセチル補酵素ＡをＴＣＡ回路で消費する。また、アセチル補酵素Ａの濃度により、脂肪酸酸化が抑制される。ミトコンドリアブロック２４は、かかるミトコンドリアの機能を模したものであり、当該ミトコンドリアブロック２４をＣＰＵ１１ａが実行することにより、解糖系ブロック２３で得られたピルビン酸からアセチル補酵素Ａが生成される反応速度ｆ_{ＰＹＲ→ＡＣｏＡ}、脂肪酸代謝ブロック２１で得られた脂肪酸アシル−補酵素Ａ複合体からアセチル補酵素Ａが生成される反応速度ｆ_{ＦＡＣ→ＡＣｏＡ}、アセチル補酵素Ａから二酸化炭素が生成される反応速度ｆ_{ＡＣｏＡ→ＣＯ２}、酸素が消費され、水が生成される反応速度ｆ_{Ｏ２→Ｈ２Ｏ}がそれぞれ演算される。また、これらの反応速度に基づき、一定時間経過後のアセチル補酵素Ａ、補酵素Ａ、酸素、及び二酸化炭素の細胞内濃度が演算される。

以下に、ミトコンドリアに関連する化学反応と、それに基づいたミトコンドリアブロック２４の具体的な演算処理について説明する。なお、以下の式におけるアセチル補酵素Ａ、補酵素Ａ、酸素、及び二酸化炭素の細胞内濃度は、それぞれ所定の初期値が与えられており、１度目の演算では初期値が用いられ、その後の演算では更新された値が用いられる。
まず、細胞内のピルビン酸が酸化され、アセチル補酵素Ａが生成される。
ＰＹＲ＋ＣｏＡ＋ＮＡＤ→ＡＣｏＡ＋ＮＡＤＨ＋ＣＯ_２ …（２７）
この反応の速度は、次式（２８）で表される。

また、脂肪酸代謝機能で生成された脂肪酸アシル−補酵素Ａ複合体が酸化され、次式（２９）の通りアセチル補酵素Ａが生成される。このルートはβ酸化と呼ばれる。

この反応の速度は、次式（３０）で表される。

ただし、Ｃ_ＭＣｏＡは、マロニルＣｏＡの濃度を示しており、ＡＣｏＡの濃度Ｃ_ＡＣｏＡに所定の係数を乗ずることにより得られる。このように、β酸化においては、マロニルＣｏＡの濃度によりアセチル補酵素Ａの生成速度が抑制される。したがって、マロニルＣｏＡの濃度により、β酸化におけるアセチルＣｏＡの生成速度が調整される。ここで、式（３０）の最初の演算では、Ｃ_ＭＣｏＡの所定の初期値が用いられ、２度目以降の演算では、演算によって得られたＣ_ＭＣｏＡの値が用いられる。

また、脂肪酸アシル−補酵素Ａ複合体は、ＴＣＡ回路により代謝され、新たにＡＴＰ及びＮＡＤＨが生成される。
ＡＣｏＡ＋ＡＤＰ＋Ｐｉ＋４ＮＡＤ→２ＣＯ_２ ＋ＣｏＡ＋ＡＴＰ＋４ＮＡＤＨ …（３１）
この反応の速度は、次式（３２）で表される。

一方、ミトコンドリアにおけるＮＡＤＨ、酸素、及びＡＤＰの消費、並びにＡＴＰの合成の関係式は以下の式（３３）ように示される。
Ｏ_２ ＋６ＡＤＰ＋６Ｐｉ＋２ＮＡＤＨ→２Ｈ_２ Ｏ＋６ＡＴＰ＋２ＮＡＤ …（３３）
この反応の速度は、次式（３４）で表される。

ミトコンドリアブロック２４では、ＣＰＵ１１ａにより、上記の式（２８），（３０），（３２），（３４）で表されるｆ_{ＰＹＲ→ＡＣｏＡ}，ｆ_{ＦＡＣ→ＡＣｏＡ}，ｆ_{ＡＣｏＡ→ＣＯ２}，ｆ_{Ｏ２→Ｈ２Ｏ}の各反応速度が演算される。

細胞内における酸素の消費速度ｆ_Ｏ２ は、前述の通り、入力された骨格筋の単位時間あたりの酸素消費量と骨格筋率とによって定まるが、次式（３５）のようにも示すことができる。

ただし、Ｃ_Ｏ２ｂは血中酸素濃度（定数）を、Ｃ_Ｏ２ は細胞内の酸素濃度を、λ_Ｏ２ は膜透過性を含む酸素の膜輸送の係数を、σ_Ｏ２ は酸素に関する分配係数をそれぞれ示す。ＣＰＵ１１ａにより、上記の式（３５）から細胞内の酸素濃度が求められる。

同様に、細胞内における二酸化炭素の産生速度ｆ_ＣＯ２は、前述の通り、入力された骨格筋の単位時間あたりの二酸化炭素産生量と骨格筋率とによって定まるが、次式のようにも示すことができる。

ただし、Ｃ_ＣＯ２ｂは血中二酸化炭素濃度（定数）を、Ｃ_ＣＯ２は細胞内の二酸化炭素濃度を、λ_ＣＯ２は膜透過性を含む二酸化炭素の膜輸送の係数を、σ_ＣＯ２は二酸化炭素に関する分配係数をそれぞれ示す。ＣＰＵ１１ａにより、上記の式（３６）から細胞内の二酸化炭素濃度が求められる。

また、ミトコンドリアブロック２４では、ＣＰＵ１１ａにより、下式（３７）で表されるアセチル補酵素Ａの生成速度、及び下式（３８）で表される補酵素Ａの生成速度がそれぞれ演算される。

さらにミトコンドリアブロック２４では、ＣＰＵ１１ａにより、上記のようにして得られた酸素、二酸化炭素、アセチル補酵素Ａ、及び補酵素Ａの生成（消費）速度から、一定時間に生成（消費）される酸素、二酸化炭素、アセチル補酵素Ａ、及び補酵素Ａの量がそれぞれ演算され、その時点での遊離脂肪酸濃度、中性脂肪濃度、脂肪酸アシル−補酵素Ａ複合体の濃度にこれらの量を反映させることで、前記一定時間経過後の細胞内の酸素濃度、二酸化炭素濃度、アセチル補酵素Ａ濃度、及び補酵素Ａ濃度が演算される。

また、ミトコンドリアブロック２４は、後述するように、第１グルコース取込速度推定処理と、第２グルコース取込速度推定処理の両方において用いられる。第１グルコース取込速度推定処理では、ミトコンドリアブロック２４を用いることにより、空腹時における細胞内の酸素濃度、二酸化炭素濃度、アセチル補酵素Ａ濃度、及び補酵素Ａ濃度を演算する処理が行われる。また、第２グルコース取込速度演算処理では、ミトコンドリアブロック２４を用いることにより、高インスリン血症時における前記各物質の細胞内濃度を演算する処理が行われる。

＜補足演算処理＞
また、ＣＰＵ１１ａにより、以下のような補足演算処理が行われる。

ＡＤＰの濃度が高いときは、クレアチンリン酸（ＰＣｒ）とＡＤＰとが反応することにより、クレアチン（Ｃｒ）とＡＴＰとが生成される。
ＰＣｒ＋ＡＤＰ→Ｃｒ＋ＡＴＰ …（３９）
この反応の速度は、次式（４０）で表される。

一方、ＡＴＰの濃度が高くなると、クレアチンとＡＴＰとによって、クレアチンリン酸とＡＤＰとが生成される。
Ｃｒ＋ＡＴＰ→ＰＣｒ＋ＡＤＰ …（４１）
この反応の速度は、次式（４２）で表される。

また、ＡＴＰは加水分解され、次式に従ってＡＤＰに変換される。
ＡＴＰ→ＡＤＰ＋Ｐｉ …（４３）
この反応の速度は、次式（４４）で表される。

補足演算処理では、ＣＰＵ１１ａにより、上記の式（４０），（４２），（４４）で表されるｆ_{ＰＣｒ→Ｃｒ}，ｆ_{Ｃｒ→ＰＣｒ}，ｆ_{ＡＴＰ→ＡＤＰ}の各反応速度が演算される。

また、ＣＰＵ１１ａにより、以下の式（４５）〜（５１）に従って、ＮＡＤ、ＮＡＤＨ、ＡＴＰ、ＡＤＰ、Ｐｉ、ＰＣｒ、Ｃｒの各生成速度が演算される。

さらに、ＣＰＵ１１ａにより、上記のようにして得られたＮＡＤ、ＮＡＤＨ、ＡＴＰ、ＡＤＰ、Ｐｉ、ＰＣｒ、Ｃｒの各生成速度から、一定時間に生成されるＮＡＤ、ＮＡＤＨ、ＡＴＰ、ＡＤＰ、Ｐｉ、ＰＣｒ、Ｃｒの量がそれぞれ演算され、その時点での細胞内におけるＮＡＤ、ＮＡＤＨ、ＡＴＰ、ＡＤＰ、Ｐｉ、ＰＣｒ、Ｃｒの濃度にこれらの量を反映させることで、前記一定時間経過後の細胞内のＮＡＤ、ＮＡＤＨ、ＡＴＰ、ＡＤＰ、Ｐｉ、ＰＣｒ、Ｃｒの各濃度が演算される。

また、補足演算処理は、後述するように、第１グルコース取込速度推定処理と、第２グルコース取込速度推定処理の両方において用いられる。第１グルコース取込速度推定処理では、補足演算処理を実行することにより、空腹時における細胞内のＮＡＤ、ＮＡＤＨ、ＡＴＰ、ＡＤＰ、Ｐｉ、ＰＣｒ、Ｃｒの各濃度を演算する処理が行われる。また、第２グルコース取込速度演算処理では、補足演算処理を実行することにより、高インスリン血症時における前記各物質の細胞内濃度を演算する処理が行われる。

［インスリン抵抗性評価支援システムの動作］
次に、本実施の形態に係るインスリン抵抗性評価支援システム１の動作について説明する。図４は、本実施の形態に係るインスリン抵抗性評価支援プログラムの処理の流れを示すフローチャートである。コンピュータ１ａは、インスリン抵抗性評価支援プログラム１４ａを実行することにより、以下のように動作する。まず、被験者は、予め体重計で体重を計測し、体組成計又は呼気ガス分析により骨格筋の単位時間あたりの酸素消費量を測定し、体組成計により骨格筋率の測定、血液検査、及び呼気ガス検査を行うことにより、被験者の体重、骨格筋の単位時間あたりの酸素消費量及び二酸化炭素産生量（以下、「酸素消費量」及び「二酸化炭素産生量」という）、骨格筋率、血中グルコース濃度、血漿インスリン濃度及び血中遊離脂肪酸濃度の各生体情報をそれぞれ取得しておく。なお、血中グルコース濃度、血漿インスリン濃度、及び血中遊離脂肪酸濃度は、それぞれ空腹時における測定値が用いられる。これらの生体情報は、インスリン抵抗性評価支援システムを使用する医師、オペレータ等のユーザに予め連絡される。

インスリン抵抗性評価支援プログラム１４ａが起動された後、まず、ＣＰＵ１１ａは、被験者の体重、酸素消費量、二酸化炭素産生量、骨格筋率、血中グルコース濃度、血漿インスリン濃度及び血中遊離脂肪酸濃度の入力をユーザに促すための入力画面を表示する（ステップＳ１）。図５は、入力画面の一例を示す模式図である。図５に示すように、入力画面３には、被験者の体重、酸素消費量、二酸化炭素産生量、骨格筋率、血中グルコース濃度、血漿インスリン濃度及び血中遊離脂肪酸濃度をそれぞれ入力するための入力エリア３１〜３７と、入力エリア３１〜３７に各生体情報を入力した後に、インスリン抵抗性の推定の実行を指示するための実行ボタン３８とが設けられている。ユーザは、入力部１３を操作することによりこれらの入力エリア３１〜３７に被験者の体重、酸素消費量、二酸化炭素産生量、骨格筋率、血中グルコース濃度、血漿インスリン濃度及び血中遊離脂肪酸濃度の各生体情報を入力し、実行ボタン３８を選択（クリック）することによりインスリン抵抗性の推定の実行を指示する。ＣＰＵ１１ａは、かかるユーザからの生体情報及び実行指示の入力を受け付ける（ステップＳ２）。かかる生体情報の入力があった場合、ＣＰＵ１１ａに割り込みが発生し、ステップＳ３以下の処理が呼び出される。

ユーザから生体情報及び実行指示の入力を受け付けると、ＣＰＵ１１ａは、第１グルコース取込速度推定処理（ステップＳ３）を実行する。図６は、第１グルコース取込速度推定処理の手順を示すフローチャートである。この第１グルコース取込速度推定処理は、空腹時におけるグルコース取込速度を推定する処理であり、ステップＳ３０１〜Ｓ３０８の処理を順次的に繰り返し実行する処理となっている。まず、ＣＰＵ１１ａは、変数に初期値をセット（初期化）する（ステップＳ３０１）。この処理では、上述した各物質の空腹時における細胞内濃度及び各反応速度の初期値をＲＡＭ１１ｃに記憶し、また、グルコース取込速度の値として、初期値０をＲＡＭ１１ｃに記憶する。次に、ＣＰＵ１１ａは、その時点でのグルコース取込速度の値を、前回のグルコース取込速度の値としてＲＡＭ１１ｃに記憶する（ステップＳ３０２）。そして、ＣＰＵ１１ａは、脂肪酸代謝ブロック２１における演算処理である脂肪酸代謝ブロック演算処理（ステップＳ３０３）、インスリンシグナルブロック２２における演算処理であるインスリンシグナルブロック演算処理（ステップＳ３０４）、解糖系ブロック２３における演算処理である解糖系ブロック演算処理（ステップＳ３０５）、ミトコンドリアブロック２４における演算処理であるミトコンドリアブロック演算処理（ステップＳ３０６）、及び補足演算処理（ステップＳ３０７）を順次実行する。かかるステップＳ３０３〜３０７において得られる物質の濃度、反応速度、ＧＬＵＴ４発現量に応じた数値及びグルコース取込速度は、空腹時における数値となる。

次にＣＰＵ１１ａは、上記の処理によって得られた空腹時におけるグルコース取込速度が定常状態に到達したか否かを判定する（ステップＳ３０８）。この処理は、本実施の形態においては、今回の演算（ターン）によって求められたグルコース取込速度の値と、ＲＡＭ１１ｃに記憶された前回の演算（ターン）のグルコース取込速度の値との差を求め、その差がグルコース取込速度が定常状態に至ったか否かを判定するための第１基準値（例えば、０．１）未満であるか否かを判定することにより行われる。そして、ＣＰＵ１１ａは、グルコース取込速度の今回値と前回値との差が第１基準値未満であるときは（ステップＳ３０８においてＹＥＳ）、メインルーチンにおける第１グルコース取込速度推定処理の呼出アドレスへ処理をリターンし、差が第１基準値以上であるときは（ステップＳ３０８においてＮＯ）、再度ステップＳ３０２以下の処理を繰り返す。

次に、ＣＰＵ１１ａは、第２グルコース取込速度推定処理を実行する（ステップＳ４）。図７は、第２グルコース取込速度推定処理の手順を示すフローチャートである。この第２グルコース取込速度推定処理は、高インスリン血症時におけるグルコース取込速度を推定する処理であり、ステップＳ４０１〜Ｓ４０８の処理を順次的に繰り返し実行する処理となっている。まず、ＣＰＵ１１ａは、インスリン濃度に所定値（高インスリン血症時のインスリン濃度）をセットする（ステップＳ４０１）。なお、第２グルコース取込速度推定処理では、インスリン濃度以外の各物質の空腹時における細胞内濃度及び各反応速度並びにグルコース取込率（ＧＬＵＴ）の初期値として、第１グルコース取込速度推定処理において最終的に得られた数値、すなわち、グルコース取込速度が定常状態に達したときの各物質の濃度及び反応速度並びにグルコース取込率がそのまま用いられる。血中グルコース濃度については、第１グルコース取込速度推定処理で使用した空腹時のグルコース濃度がそのまま用いられる。また、グルコース取込速度の値として、初期値０がＲＡＭ１１ｃに記憶される。次に、ＣＰＵ１１ａは、その時点でのグルコース取込速度の値を、前回のグルコース取込速度の値としてＲＡＭ１１ｃに記憶する（ステップＳ４０２）。そして、ＣＰＵ１１ａは、脂肪酸代謝ブロック２１における演算処理である脂肪酸代謝ブロック演算処理（ステップＳ４０３）、インスリンシグナルブロック２２における演算処理であるインスリンシグナルブロック演算処理（ステップＳ４０４）、解糖系ブロック２３における演算処理である解糖系ブロック演算処理（ステップＳ４０５）、ミトコンドリアブロック２４における演算処理であるミトコンドリアブロック演算処理（ステップＳ４０６）、及び補足演算処理（ステップＳ４０７）を順次実行する。かかるステップＳ４０３〜４０７において得られる物質の濃度、反応速度、ＧＬＵＴ４発現量に応じた数値及びグルコース取込速度は、高インスリン血症時における数値となる。

次にＣＰＵ１１ａは、上記の処理によって得られた高インスリン血症時におけるグルコース取込速度が定常状態に到達したか否かを判定する（ステップＳ４０８）。この処理は、本実施の形態においては、今回の演算（ターン）によって求められたグルコース取込速度の値と、ＲＡＭ１１ｃに記憶された前回の演算（ターン）のグルコース取込速度の値との差を求め、その差がグルコース取込速度が定常状態に至ったか否かを判定するための第１基準値（例えば、０．１）未満であるか否かを判定することにより行われる。なお、第１グルコース取込速度推定処理と第２グルコース取込速度推定処理とにおいて、同じ第１の基準値を用いてグルコース取込速度が定常状態に至ったか否かを判定する構成としたが、これに限られず、それぞれで異なる基準値を用いてもよい。そして、ＣＰＵ１１ａは、グルコース取込速度の今回値と前回値との差が第１基準値未満であるときは（ステップＳ４０８においてＹＥＳ）、メインルーチンにおける第２グルコース取込速度推定処理の呼出アドレスへ処理をリターンし、差が第１基準値以上であるときは（ステップＳ４０８においてＮＯ）、再度ステップＳ４０２以下の処理を繰り返す。

次に、ＣＰＵ１１ａは、入力された被験者の体重及び骨格筋率から、被験者の筋肉量を算出し（ステップＳ５）、この筋肉量、グルコース分布容量（定数）、及び第２グルコース取込速度推定処理によって得られたグルコース濃度を、第２グルコース取込速度推定処理によって得られたグルコース取込率（ＧＬＵＴ）に乗じ、単位重量当たりの高インスリン血症時のグルコース取込速度を算出する（ステップＳ６）。次に、ＣＰＵ１１ａは、インスリン抵抗性の有無を推定する（ステップＳ７）。この処理では、ステップＳ６で得た単位重量当たりのグルコース取込速度の値（推定値）が、インスリン抵抗性の有無を推定するための第２基準値（例えば、１２ｍｇ／ｋｇ／ｍｉｎ）以上か否かを判定することにより行われる。これにより、グルコース取込速度の推定値が前記第２基準値以上であれば、インスリン抵抗性がない、即ち、インスリン感受性がある、と推定することができ、グルコース取込速度の推定値が前記第２基準値未満であれば、インスリン抵抗性がある、即ち、インスリン感受性がない、と推定することができる。このようにして、ＣＰＵ１１ａは、インスリン抵抗性の推定を行う。

次にＣＰＵ１１ａは、インスリン抵抗性の推定結果の出力画面を表示する（ステップＳ４）。図８は、出力画面の一例を示す模式図である。図８に示すように、出力画面４には、上記の処理によって得られたインスリン抵抗性の推定結果４１と、最終的に得られた単位重量当たりのグルコース取込速度の推定値４２とが含まれる。この出力画面４により、ユーザに対して、インスリン抵抗性の推定結果及びグルコース取込速度の推定値が通知される。このようにインスリン抵抗性の推定結果及びグルコース取込速度の推定値をユーザに提供することにより、ユーザはこれらの情報を活用してインスリン抵抗性の評価を行うことができる。

上記のような構成とすることにより、被験者の負担の大きいグルコースクランプ及び経口糖負荷試験の検査結果を必要とせず、被験者の負担が少なく、簡便な検査によって得ることができる「体重」、「骨格筋の単位時間あたりの酸素消費量」、「骨格筋の単位時間あたりの二酸化炭素生成量」、「骨格筋率」、「空腹時の血中グルコース濃度」、「空腹時の血漿インスリン濃度」、及び「空腹時の血中遊離脂肪酸濃度」の各生体情報によってインスリン抵抗性の有無を推定することができる。このような入力情報の中でも、空腹時の血中グルコース濃度及び血漿インスリン濃度については、上記のように血液検査によって簡便に得られる検査値を用いることができるが、グルコースクランプ又は経口糖負荷試験によって得られた検査値を入力情報として用いることもできる。しかし、グルコースクランプ及び経口糖負荷試験のような被験者にとって負担の大きい検査結果を常に必要とする訳ではなく、これらの検査結果でなくとも、簡便な血液検査の結果を利用することができるという点、及び、入力情報に用いる情報を、血液検査、グルコースクランプ、及び経口糖負荷試験のいずれからも選択することができるという点で、本システムは有用である。

また、マロニルＣｏＡの濃度によるβ酸化におけるアセチルＣｏＡの生成速度の抑制を考慮しており、アセチルＣｏＡの濃度に基づいて求められたマロニルＣｏＡにより、前記β酸化におけるアセチルＣｏＡの生成速度を調整するように構成し、実際の生体におけるグルコースの代謝機能を忠実に再現した。このことからも、実際のグルコース取込速度を精度よく反映した推定値が得られることが期待できる。

また、上記の生体情報に基づき、実際の生体器官におけるグルコースの代謝機能を仮想的に再現し、グルコース取込速度の推定値を求めるため、実際のグルコース取込速度を精度よく反映した推定値が得られることが期待できる。また、本実施の形態では、高インスリン血症時におけるグルコース取込速度を推定し、この推定値を出力する。グルコースクランプ試験では、高インスリン血症の状態が維持されるようにグルコース及びインスリンが被験者に注入される。また、文献等に開示されるグルコース取込速度も、高インスリン血症時における値が一般的である。したがって、上述のように高インスリン血症時におけるグルコース取込速度を推定し、その推定値を出力することで、ユーザは上記推定値をグルコースクランプ試験の結果又は文献等に報告されている数値と比較することができ、その被験者のインスリン抵抗性を簡便に評価することができる。

さらに、本実施の形態では、単位重量当たりのグルコース取込速度を推定し、これを出力する構成とした。グルコースクランプ試験では、単位重量当たりのグルコース取込速度が得られる。また、文献等に開示されるグルコース取込速度も、単位重量当たりのものが一般的である。したがって、上述のように単位重量当たりのグルコース取込速度を推定し、その推定値を出力することで、ユーザは上記推定値をグルコースクランプ試験の結果又は文献等に報告されている数値と比較することができ、その被験者のインスリン抵抗性を簡便に評価することができる。

（実施の形態２）
本実施の形態は、被験者を測定することにより得られた血中のグルコース濃度、インスリン濃度、及び遊離脂肪酸濃度に関する情報の入力を受け付け、入力を受け付けた血中のグルコース濃度、インスリン濃度、及び遊離脂肪酸濃度に関する情報に基づいて、被験者のグルコース取込速度を推定し、推定されたグルコース取込速度に基づいて、被験者のインスリン抵抗性を推定し、推定結果を出力するインスリン抵抗性評価支援システムである。

［インスリン抵抗性評価支援システムの構成］
図９は、本実施の形態に係るインスリン抵抗性評価支援システムの構成を示すブロック図である。本実施の形態に係るインスリン抵抗性評価支援システム２０１は、コンピュータ２０１ａによって実現される。図９に示すように、コンピュータ２０１ａは、本体２１１と、画像表示部２１２と、入力部２１３とを備えている。本体２１１は、ＣＰＵ２１１ａと、ＲＯＭ２１１ｂ、ＲＡＭ２１１ｃ、ハードディスク２１１ｄ、読出装置２１１ｅ、入出力インタフェース２１１ｆ、及び画像出力インタフェース２１１ｇを備えており、ＣＰＵ２１１ａ、ＲＯＭ２１１ｂ、ＲＡＭ２１１ｃ、ハードディスク２１１ｄ、読出装置２１１ｅ、入出力インタフェース２１１ｆ、および画像出力インタフェース２１１ｇは、バス２１１ｉによって接続されている。

ＣＰＵ２１１ａは、ＲＡＭ２１１ｃにロードされたコンピュータプログラムを実行することが可能である。そして、後述するようなインスリン抵抗性評価支援プログラム２１４ａを当該ＣＰＵ２１１ａが実行することにより、コンピュータ２０１ａがインスリン抵抗性評価支援システム２０１として機能する。

ＲＡＭ１１ｃは、ＳＲＡＭまたはＤＲＡＭ等によって構成されている。ＲＡＭ１１ｃは、ハードディスク２１１ｄに記録されているインスリン抵抗性評価支援プログラム２１４ａの読み出しに用いられる。また、ＣＰＵ２１１ａがコンピュータプログラムを実行するときに、ＣＰＵ２１１ａの作業領域として利用される。

ハードディスク２１１ｄは、オペレーティングシステムおよびアプリケーションプログラム等、ＣＰＵ２１１ａに実行させるための種々のコンピュータプログラムおよび当該コンピュータプログラムの実行に用いられるデータがインストールされている。後述するインスリン抵抗性評価支援プログラム２１４ａも、このハードディスク２１１ｄにインストールされている。

読出装置２１１ｅは、フレキシブルディスクドライブ、ＣＤ−ＲＯＭドライブ、またはＤＶＤ−ＲＯＭドライブ等によって構成されており、可搬型記録媒体２１４に記録されたコンピュータプログラムまたはデータを読み出すことができる。また、可搬型記録媒体２１４には、コンピュータをインスリン抵抗性評価支援システムとして機能させるためのインスリン抵抗性評価支援プログラム２１４ａが格納されており、コンピュータ２０１ａが当該可搬型記録媒体２１４からインスリン抵抗性評価支援プログラム２１４ａを読み出し、当該インスリン抵抗性評価支援プログラム２１４ａをハードディスク２１１ｄにインストールすることが可能である。

なお、前記インスリン抵抗性評価支援プログラム２１４ａは、可搬型記録媒体２１４によって提供されるのみならず、電気通信回線（有線、無線を問わない）によってコンピュータ２０１ａと通信可能に接続された外部の機器から前記電気通信回線を通じて提供することも可能である。例えば、前記インスリン抵抗性評価支援プログラム２１４ａがインターネット上のサーバコンピュータのハードディスク内に格納されており、このサーバコンピュータにコンピュータ２０１ａがアクセスして、当該コンピュータプログラムをダウンロードし、これをハードディスク２１１ｄにインストールすることも可能である。

また、ハードディスク２１１ｄには、例えば米マイクロソフト社が製造販売するWindows（登録商標）等のマルチタスクオペレーティングシステムがインストールされている。以下の説明においては、本実施の形態に係るインスリン抵抗性評価支援プログラム２１４ａは当該オペレーティングシステム上で動作するものとしている。なお、インスリン抵抗性評価支援プログラム２１４ａの構成の詳細については後述する。

なお、本実施の形態２に係るコンピュータ２０１ａのその他の構成は、実施の形態１に係るコンピュータ１ａの構成と同様であるので、その説明を省略する。

［インスリン抵抗性評価支援システムの機能的な構成］
次に、インスリン抵抗性評価支援プログラム２１４ａについて更に詳細に説明する。図１０は、本実施の形態に係るインスリン抵抗性評価支援システムにおける仮想的な物質の反応の流れを示す概念図である。インスリン抵抗性評価支援プログラム２１４ａは、被験者の体重、骨格筋率、空腹時の血中グルコース濃度、血漿インスリン濃度、及び血中遊離脂肪酸濃度の各生体情報を入力とし、被験者の末梢組織（骨格筋）の糖取込速度（糖取込率）の推定値を出力とする。かかるインスリン抵抗性評価支援プログラム２１４ａは、身体の糖取込に関連する生体器官（機能）を機能別に分けたそれぞれの機能を仮想的に再現する４つのブロック２２１〜２２４を含んでいる。それぞれのブロック２２１〜２２４は、複数のパラメータを含んでおり、グルコースの取込に関連する物質の生成の反応速度を演算するように構成されている。また、インスリン抵抗性評価支援プログラム２１４ａは、上記の反応速度から、グルコースの取込に関連する物質の生成速度（一定時間あたりの生成量）を演算し、この演算によって得られた一定時間後の生成量を反映した、一定時間後の物質濃度を演算するように構成されている。以下に、各ブロック２２１〜２２４の詳細な構成について説明する。

＜脂肪酸代謝ブロック＞
脂肪酸代謝ブロック２２１は、筋肉組織や脂肪組織などの生体器官の脂肪酸代謝機能を仮想的に再現する機能ブロックである。身体の脂肪酸代謝機能は、血中の遊離脂肪酸を細胞内に取り込んで遊離脂肪酸（ＦＦＡ）から脂肪酸アシル−補酵素Ａ複合体（ＦＡＣ）を生成し、これとともに、ジアシルグリセロール（ＤＧ）を介して中性脂肪（ＴＧ）を生成する。脂肪酸代謝ブロック２２１は、かかる身体の脂肪酸代謝機能を模したものであり、当該脂肪酸代謝ブロック２２１をＣＰＵ１１ａが実行することにより、ＦＡＣとＤ−グリセルアルデヒド３-リン酸(ＧＡ３Ｐ)からＤＧを生成する反応速度ｆ_{ＧＡ３Ｐ→ＤＧ}、ＤＧとＦＡＣからＴＧを生成する反応速度ｆ_{ＤＧ→ＴＧ}、ＴＧからＤＧとＦＦＡを生成する反応速度ｆ_{ＴＧ→ＤＧ}、ＤＧからＦＦＡとグリセロール（ＧＬＲ）を生成する反応速度ｆ_{ＤＧ→ＦＦＡ}、ＦＦＡからＦＡＣが生成される反応速度ｆ_{ＦＦＡ→ＦＡＣ}、が演算され、また、一定時間経過後の細胞内ＦＦＡ、ＧＬＲ, ＤＧ、ＴＧ、ＦＡＣ及びＦＡＣ濃度が演算される。

以下に、筋肉組織及び脂肪組織等の生体器官の脂肪酸代謝機能に関連する化学反応、並びにそれに基づいた脂肪酸代謝ブロック２２１の具体的な演算処理について説明する。なお、以下の式におけるＦＦＡ、ＦＡＣ、ＤＧ、ＴＧ、ＧＡ３Ｐ、ＧＬＲ,ＡＴＰ、ＡＤＰ、ＮＡＤ, ＮＡＤＨ、及び無機リン酸(Ｐｉ)の細胞内濃度は、それぞれ所定の初期値が与えられており、１度目の演算では初期値が用いられ、その後の演算では更新された値が用いられる。

まず、血中の脂肪酸（ＦＦＡ）が細胞内に取り込まれる。この取り込み速度ｆ_ＦＦＡは以下の式（５２）で表される。

ただし、Ｃ_{ＦＦＡ，ｂ}は、入力された血中の脂肪酸濃度を、Ｃ_ＦＦＡは細胞内の脂肪酸濃度を、σ_ＦＦＡは脂肪酸に関する分配係数を、Ｑは筋肉組織の血流速をそれぞれ示す。

次に、血中のグリセロール（ＧＬＲ）が細胞内に取り込まれる。この取り込み速度ｆ_ＧＬＣは以下の式（５３）で表される。

ただし、Ｃ_{ＧＬＲ，ｂ}は、入力された血中のグリセロール濃度を、Ｃ_ＧＬＲは細胞内のグリセロール濃度を、σ_ＧＬＣはグリセロールに関する分配係数をそれぞれ示す。

実施の形態１の式（８）にしたがって、細胞内に流入した脂肪酸（ＦＦＡ）とＣｏＡからＡＴＰを使いＦＡＣを合成する。また、この反応の速度は、次式（５４）で示される。

次に、解糖系ブロックのＧＡ３Ｐと生成されたＦＡＣから複数の反応を経てＤＧを合成する。
ＧＡ３Ｐ＋２ＦＡＣ＋ＮＡＤＨ→ＤＧ＋２ＣｏＡ＋Ｐｉ＋ＮＡＤ …（５５）
また、この反応の速度は、次式（５６）で示される。

次に、ＤＧとＦＡＣからＴＧが生成される。
ＤＧ＋ＦＡＣ→ＴＧ＋ＣｏＡ …（５７）
また、この反応の速度は、次式（５８）で示される。

中性脂肪（ＴＧ）は、次式の通り、ホルモンセンシティブリパーゼ（ＨＳＬ）などによってＤＧとＦＦＡに分解される。
ＴＧ→ＤＧ＋ＦＦＡ …（５９）
また、この反応の速度は、次式（６０）で示される。

ただし、ＨＳＬ（ｉｎｓｕｌｉｎ）はインスリンシグナリングブロック２２２にて演算される。

ＤＧはＨＳＬなどによってＦＦＡとＧＬＲに分解される。
ＤＧ→ＧＬＲ＋２ＦＦＡ …（６１）
また、この反応の速度は、次式（６２）で示される。

脂肪酸代謝ブロック２２１では、ＣＰＵ２１１ａにより、上記の式（５２）、（５３）、（５４），（５６），（５８）、（６０）、（６２）で表されるｆ_ＦＦＡ，ｆ_ＧＬＣ，ｆ_{ＦＦＡ→ＦＡＣ}，ｆ_{ＧＡ３Ｐ→ＤＧ}，ｆ_{ＤＧ→ＴＧ}，ｆ_{ＴＧ→ＤＧ}，ｆ_{ＤＧ→ＧＬＲ}の各反応速度が演算される。

また、脂肪酸代謝ブロック２２１では、ＣＰＵ２１１ａにより、下式（６３）で表されるＦＦＡの生成速度、下式（６４）で表されるＦＡＣの生成速度、下式（６５）で表されるＤＧの生成速度、下式（６６）で表されるＴＧの生成速度、下式（６７）で表されるＧＬＣの生成速度がそれぞれ演算される。

ただし、式（６４）において、反応速度ｆ_{ＦＡＣ→ＡＣｏＡ}は、後述する式（１５２）で表され、ミトコンドリアブロック２２４において算出される。ミトコンドリアブロック２２４の演算が一度も行われていないときには、反応速度ｆ_{ＦＡＣ→ＡＣｏＡ}の初期値が用いられる。

さらに脂肪酸代謝ブロック２２１では、ＣＰＵ２１１ａにより、上記のようにして得られたＦＦＡ、ＦＡＣ、ＤＧ、ＴＧ及びＧＬＲの生成速度から、一定時間に生成されるＦＦＡ、ＦＡＣ、ＤＧ、ＴＧ及びＧＬＲの量がそれぞれ演算され、その時点での細胞内におけるＦＦＡ濃度、ＦＡＣ濃度、ＤＧ濃度、ＴＧ濃度、ＧＬＲ濃度にこれらの量を反映させることで、前記一定時間経過後の細胞内のＦＦＡ濃度、ＦＡＣ濃度、ＤＧ濃度、ＴＧ濃度及びＧＬＲ濃度が演算される。

また、脂肪酸代謝ブロック２２１は、後述するように、第１グルコース取込速度推定処理と、第２グルコース取込速度推定処理の両方において用いられる。第１グルコース取込速度推定処理では、脂肪酸代謝ブロック２２１を用いることにより、空腹時における細胞内のグルコース濃度、ＦＦＡ濃度、ＦＡＣ濃度、ＤＧ濃度、ＴＧ濃度及びＧＬＣ酸濃度を演算する処理が行われる。また、第２グルコース取込速度演算処理では、脂肪酸代謝ブロック２２１を用いることにより、高インスリン血症時における前記各物質の細胞内濃度を演算する処理が行われる。

＜インスリンシグナルブロック＞
インスリンシグナルブロック２２２は、生体器官、特に筋肉組織、脂肪組織において発現が知られているグルコーストランスポータ（ＧＬＵＴ４）の発現量調整機能を仮想的に再現する機能ブロックである。生体のＧＬＵＴ４発現量調整機能は、インスリン受容体に結合したインスリン量と、脂肪酸代謝の代謝生成物であるジアシルグリセロール（ＤＧ）の濃度に応じてグルコース輸送担体（ＧＬＵＴ４）の細胞膜上に占める割合を調整する。本実施の形態に係るインスリンシグナルブロック２２２は、かかる身体のＧＬＵＴ４の細胞膜上に占める割合や、インスリン作用が知られているヘキソキナーゼ（ＨＫ）、グリコーゲンシンターゼ（ＧＳ）、ピルビン酸デヒドロゲナーゼ（ＰＤＨ）、ホルモン感受性リパーゼ（ＨＳＬ）の各酵素と、脂肪酸酸化の調整機能を模したものであり、当該インスリンシグナルブロック２２２をＣＰＵ２１１ａが実行することにより、血漿インスリン濃度（ＰＩ）とＤＧ濃度に基づいて、ＧＬＵＴ４の発現量ならびにＨＫ，ＧＳ，ＰＤＨ，ＨＳＬの各酵素と脂肪酸酸化の活性に応じた数値が演算される。

以下に、生体器官のＧＬＵＴ４の細胞膜上発現量調整機能に基づいたインスリンシグナルブロック２２２の具体的な演算処理について説明する。生体器官のＧＬＵＴ４発現量調整機能は、血漿インスリン濃度の上昇に従い、ＧＬＵＴ４の発現量が増加し、ＦＡＣ濃度及びＤＧ濃度の上昇に従い、ＧＬＵＴ４の発現量が抑制されるという特徴を有している。かかる特徴を考慮し、インスリンシグナルブロック２２２においては、下式（６８）〜（１３１）従って、入力された血漿インスリン濃度と、脂肪酸代謝ブロック２１の演算で得られたＤＧ濃度とを用いて、ＧＬＵＴ４の発現量に応じた数値であるグルコースの取込率（ＧＬＵＴ）が演算される。

細胞膜上のインスリン未結合のインスリン受容体から、細胞膜上の１つのインスリンが結合したインスリン受容体へ遷移する速度ｆ_{ｘ２→ｘ３}は、下式（６８）のように血漿インスリン濃度ｘ_１と細胞膜上のインスリン未結合インスリン受容体濃度ｘ２及び反応速度定数ｋ_１によって決まる。

細胞膜上の１つのインスリンが結合したインスリン受容体から細胞膜上のインスリン未結合のインスリン受容体へ遷移する速度ｆ_{ｘ３→ｘ２}は、下式（６９）のように細胞膜上の１つのインスリンが結合したインスリン受容体濃度ｘ_３及び反応速度定数ｋ_−１によって決まる。

細胞膜上の１つのインスリンが結合したインスリン受容体から細胞膜上の１つのインスリンが結合しリン酸化により活性化されたインスリン受容体へ遷移する速度ｆ_{ｘ３→ｘ５}は、下式（７０）のように細胞膜上の１つのインスリンが結合したインスリン受容体濃度ｘ_３及び反応速度定数ｋ_３によって決まる。

細胞膜上の１つのインスリンが結合しリン酸化により活性化されたインスリン受容体から細胞膜上のインスリン未結合のインスリン受容体へ遷移する速度ｆ_{ｘ５→ｘ２}は、下式（７１）のように細胞膜上の１つのインスリンが結合しリン酸化により活性化されたインスリン受容体濃度ｘ_５、プロテインチロシンフォスファターゼ(ＰＴＰ)及び反応速度定数ｋ_−３によって決まる。

細胞膜上の１つのインスリンが結合しリン酸化により活性化されたインスリン受容体から細胞膜上の２つのインスリンが結合しリン酸化により活性化されたインスリン受容体へ遷移する速度ｆ_{ｘ５→ｘ４}は、下式（７２）のように細胞膜上の１つのインスリンが結合しリン酸化により活性化されたインスリン受容体濃度ｘ_５及び反応速度定数ｋ_２によって決まる。

細胞膜上の２つのインスリンが結合しリン酸化により活性化されたインスリン受容体から細胞膜上の１つのインスリンが結合しリン酸化により活性化されたインスリン受容体へ遷移する速度ｆ_{ｘ４→ｘ５}は、下式（７３）のように細胞膜上の１つのインスリンが結合しリン酸化により活性化されたインスリン受容体濃度ｘ_４及び反応速度定数ｋ_−２によって決まる。

細胞膜上のインスリン未結合のインスリン受容体から、細胞内インスリン未結合のインスリン受容体へ遷移する速度ｆ_{ｘ２→ｘ６}は、下式（７４）のように細胞膜上のインスリン未結合インスリン受容体濃度ｘ_２及び反応速度定数ｋ_４によって決まる。

細胞内インスリン未結合のインスリン受容体から細胞膜上の１つのインスリンが結合したインスリン受容体へ遷移する速度ｆ_{ｘ６→ｘ２}は、下式（７５）のように細胞内インスリン未結合のインスリン受容体濃度ｘ_６及び反応速度定数ｋ_−４によって決まる。

細胞膜上の１つのインスリンが結合しリン酸化により活性化されたインスリン受容体から細胞内の１つのインスリンが結合しリン酸化により活性化されたインスリン受容体へ遷移する速度ｆ_{ｘ５→ｘ８}は、下式（７６）のように細胞膜上の１つのインスリンが結合しリン酸化により活性化されたインスリン受容体濃度ｘ_５及び反応速度定数ｋ_４’によって決まる。

細胞内の１つのインスリンが結合しリン酸化により活性化されたインスリン受容体から細胞膜上の１つのインスリンが結合しリン酸化により活性化されたインスリン受容体へ遷移する速度ｆ_{ｘ８→ｘ５}は、下式（７７）のように細胞内の１つのインスリンが結合しリン酸化により活性化されたインスリン受容体濃度ｘ_８及び反応速度定数ｋ_―４’によって決まる。

細胞膜上の２つのインスリンが結合しリン酸化により活性化されたインスリン受容体から細胞内の２つのインスリンが結合しリン酸化により活性化されたインスリン受容体へ遷移する速度ｆ_{ｘ４→ｘ７}は、下式（７８）のように細胞膜上の２つのインスリンが結合しリン酸化により活性化されたインスリン受容体濃度ｘ_４及び反応速度定数ｋ_４’によって決まる。

細胞内の２つのインスリンが結合しリン酸化により活性化されたインスリン受容体から細胞膜上の２つのインスリンが結合しリン酸化により活性化されたインスリン受容体へ遷移する速度ｆ_{ｘ７→ｘ４}は、下式（７９）のように細胞内の２つのインスリンが結合しリン酸化により活性化されたインスリン受容体濃度ｘ_７及び反応速度定数ｋ_−４’によって決まる。

細胞内の２つのインスリンが結合しリン酸化により活性化されたインスリン受容体から細胞内のインスリン未結合のインスリン受容体へ遷移する速度ｆ_{ｘ７→ｘ６}は、下式（８０）のように細胞内の２つのインスリンが結合しリン酸化により活性化されたインスリン受容体濃度ｘ_７、プロテインチロシンフォスファターゼ(ＰＴＰ)及び反応速度定数ｋ_６によって決まる。

細胞内の１つのインスリンが結合しリン酸化により活性化されたインスリン受容体から細胞内のインスリン未結合のインスリン受容体へ遷移する速度ｆ_{ｘ８→ｘ６}は、下式（８１）のように細胞内の１つのインスリンが結合しリン酸化により活性化されたインスリン受容体濃度ｘ_８、プロテインチロシンフォスファターゼ(ＰＴＰ)及び反応速度定数ｋ_６によって決まる。

細胞内におけるインスリン受容体生成速度は、下式（８２）のように定数ｋ５で表される一定の値とする。

細胞内のインスリン未結合のインスリン受容体の分解速度ｆ_ｘ６→は、下式（８３）のように細胞内のインスリン未結合のインスリン受容体濃度ｘ_６及び反応速度定数ｋ_−５によって決まる。

未リン酸化インスリン受容体基質１（ＩＲＳ１）からチロシンリン酸化インスリン受容体基質１へ遷移する速度ｆ_{ｘ９→ｘ１０}は、下式（８４）のように未リン酸化インスリン受容体基質１濃度ｘ_９、胞膜上の２つのインスリンが結合しリン酸化により活性化されたインスリン受容体ｘ_４、胞膜上の１つのインスリンが結合しリン酸化により活性化されたインスリン受容体ｘ_５及びチロシンリン酸化反応速度定数ｋ_７によって決まる。

ただし、ＩＲ_ｐは最大のインスリン刺激下における細胞膜上のリン酸化されたインスリン受容体濃度をあらわす。

チロシンリン酸化インスリン受容体基質１から未リン酸化インスリン受容体基質１（ＩＲＳ１）へ遷移する速度ｆ_{ｘ１０→ｘ９}は、下式（８５）のようにチロシンリン酸化インスリン受容体基質１ｘ_１０、プロテインチロシンフォスファターゼ(ＰＴＰ)及び脱リン酸化反応速度定数ｋ_−７によって決まる。

未リン酸化インスリン受容体基質１（ＩＲＳ１）からセリンリン酸化インスリン受容体基質１へ遷移する速度ｆ_{ｘ９→ｘ１０ａ}は、下式（８６）のように未リン酸化インスリン受容体基質１濃度ｘ_９、プロテインキナーゼＣ（ＰＫＣ）及びセリンリン酸化反応速度定数ｋ_７’によって決まる。

セリンリン酸化インスリン受容体基質１から未リン酸化インスリン受容体基質１（ＩＲＳ１）へ遷移する速度ｆ_{ｘ１０ａ→ｘ９}は、下式（８７）のようにセリンリン酸化インスリン受容体基質１ｘ_１０ａ及び脱リン酸化反応速度定数ｋ_−７’によって決まる。

セリンリン酸化インスリン受容体基質１とフォスファチジルイノシトール３キナーゼ（ＰＩ３Ｋ）からセリンリン酸化インスリン受容体基質１とフォスファチジルイノシトール３キナーゼ（ＰＩ３Ｋ）の複合体が生成される速度ｆ_{ｘ１０→ｘ１２}は、下式（８８）のようにセリンリン酸化インスリン受容体基質１濃度ｘ_１０、フォスファチジルイノシトール３キナーゼ（ＰＩ３Ｋ）濃度ｘ_１１及び速度定数ｋ_８によって決まる。

セリンリン酸化インスリン受容体基質１とフォスファチジルイノシトール３キナーゼ（ＰＩ３Ｋ）の複合体からセリンリン酸化インスリン受容体基質１とフォスファチジルイノシトール３キナーゼ（ＰＩ３Ｋ）のそれぞれが生成される速度ｆ_{ｘ１２→ｘ１０}は、下式（８９）のようにセリンリン酸化インスリン受容体基質１とフォスファチジルイノシトール３キナーゼ（ＰＩ３Ｋ）の複合体濃度ｘ_１２及び速度定数ｋ_−８によって決まる。

フォスファチジルイノシトール（３，４，５）三リン酸が、Src homology 2 containinginositol 5' phosphatase(ＳＨＩＰ)などの５’脂質フォスファターゼによって、フォスファチジルイノシトール（３，４）二リン酸へ脱リン酸化される速度ｆ_{ｘ１３→ｘ１５}は、下式（９０）のようにフォスファチジルイノシトール（３，４，５）三リン酸の割合ｘ_１３、ＳＨＩＰ濃度、及び速度定数ｋ_−１０によって決まる。

フォスファチジルイノシトール（３，４）二リン酸がフォスファチジルイノシトール（３，４，５）三リン酸へリン酸化される速度ｆ_{ｘ１５→ｘ１３}は、下式（９１）のようにフォスファチジルイノシトール（３，４）二リン酸の割合ｘ_１５及び速度定数ｋ_１０によって決まる。

フォスファチジルイノシトール（３，４，５）三リン酸が、Phosphatase and Tensin Homolog Deleted from Chromosome 10 (ＰＴＥＮ)などの３’脂質フォスファターゼによって、フォスファチジルイノシトール（４，５）二リン酸へ脱リン酸化される速度ｆ_{ｘ１３→ｘ１４}は、下式（９２）のようにフォスファチジルイノシトール（３，４，５）三リン酸の割合ｘ_１３、ＰＴＥＮ濃度、及び速度定数ｋ_−９によって決まる。

なお、ｋ_９（ｓｔｉｍｕｌａｔｅｄ）とは、フォスファチジルイノシトール（４，５）二リン酸がフォスファチジルイノシトール（３，４，５）三リン酸へリン酸化される速度ｆ_{ｘ１４→ｘ１３}に関する定数であって、当該速度におけるセリンリン酸化インスリン受容体基質１とフォスファチジルイノシトール３キナーゼ（ＰＩ３Ｋ）との複合体に対して依存している成分に関するものである。

フォスファチジルイノシトール（４，５）二リン酸がフォスファチジルイノシトール（３，４，５）三リン酸へリン酸化される速度ｆ_{ｘ１４→ｘ１３}は、下式（９３）のようにフォスファチジルイノシトール（４，５）二リン酸の割合ｘ_１４、及び速度定数ｋ_９によって決まる。

なお、ｋ_９（ｂａｓａｌ）とは、フォスファチジルイノシトール（４，５）二リン酸がフォスファチジルイノシトール（３，４，５）三リン酸へリン酸化される速度ｆ_{ｘ１４→ｘ１３}に関する定数であって、当該速度におけるセリンリン酸化インスリン受容体基質１とフォスファチジルイノシトール３キナーゼ（ＰＩ３Ｋ）との複合体に対して依存していない成分に関するものである。

未リン酸化ＡＫＴがリン酸化ＡＫＴへリン酸化される速度ｆ_{ｘ１６→ｘ１７}は、下式（９４）のように未リン酸化ＡＫＴの割合ｘ_１６及び速度定数ｋ_１１によって決まる。

リン酸化ＡＫＴが未リン酸化ＡＫＴへ脱リン酸化される速度ｆ_{ｘ１７→ｘ１６}は、下式（９５）のようにリン酸化ＡＫＴの割合ｘ_１７及び速度定数ｋ_−１１によって決まる。

未リン酸化ＰＫＣがリン酸化ＰＫＣへリン酸化される速度ｆ_{ｘ１８→ｘ１９}は、下式（９６）のように未リン酸化ＰＫＣの割合ｘ_１８及び速度定数ｋ_１２によって決まる。

リン酸化ＰＫＣから未リン酸化ＰＫＣへ脱リン酸化される速度ｆ_{ｘ１９→ｘ１８}は、下式（９７）のようにリン酸化ＰＫＣの割合ｘ_１９及び速度定数ｋ_−１２によって決まる。

細胞内ＧＬＵＴ４が生成される速度ｆ_→ｘ２０は、下式（９８）のように生成速度ｋ_１４によって決まる。

細胞内ＧＬＵＴ４が分解される速度ｆ_ｘ２０→は、下式（９９）のように細胞内ＧＬＵＴ４の割合ｘ_２０及び速度定数ｋ_−１４によって決まる。

細胞内グルコース輸送担体ＧＬＵＴ４から細胞膜上のＧＬＵＴ４への移動速度ｆ_{ｘ２０→ｘ２１}は、下式（１００）のように細胞内ＧＬＵＴ４の割合ｘ_２０及び速度定数ｋ_１３、ｋ_１３’によって決まる。

細胞膜上グルコース輸送担体ＧＬＵＴ４から細胞内ＧＬＵＴ４への移動速度ｆ_{ｘ２１→ｘ２０}は、下式（１０１）のように細胞膜上ＧＬＵＴ４の割合ｘ_２１及び速度定数ｋ_−１３によって決まる。

非活性酵素から活性化された酵素ヘの活性化速度ｆ_{ｘ２２→ｘ２３}は、下式（１０２）のように酵素活性補助変数ｘ_２２及び速度定数ｋ_１５によって決まる。

活性化された酵素から非活性酵素ヘの非活性化速度ｆ_{ｘ２３→ｘ２２}は、下式（１０３）のように酵素活性補助変数ｘ_２３及び速度定数ｋ_−１５によって決まる。

ヘキソキナーゼ活性
活性化されたヘキソキナーゼ（ＨＫ）の値は、下式（１０４）によって演算される。

グリコーゲンシンターゼ活性
活性化されたグリコーゲンシンターゼ（ＧＳ）の値は、下式（１０５）によって演算される。

ピルビン酸デヒドロゲナーゼ活性
活性化されたピルビン酸デヒドロゲナーゼ（ＰＤＨ）の値は、下式（１０６）によって演算される。

ホルモンセンシティブリパーゼ活性
活性化されたホルモンセンシティブリパーゼ（ＨＳＬ）の値は、下式（１０７）によって演算される。

脂肪酸酸化（ベータ酸化）の活性
活性化されたヒドロキシアシルデヒドロゲナーゼ（ＨＡＤ）の値は、下式（１０８）によって演算される。

またインスリンシグナリングブロック２２２では、下式（１０９）で表される細胞膜上のインスリン未結合のインスリン受容体生成速度、下式（１１０）で表される細胞膜上の１つのインスリンが結合したインスリン受容体生成速度、下式（１１１）で表される細胞膜上の１つのインスリンが結合しリン酸化により活性化されたインスリン受容体生成速度、下式（１１２）で表される細胞膜上の２つのインスリンが結合しリン酸化により活性化されたインスリン受容体生成速度、下式（１１３）で表される胞膜内インスリン未結合のインスリン受容体生成速度、下式（１１４）で表される細胞内の１つのインスリンが結合しリン酸化により活性化されたインスリン受容体生成速度、下式（１１５）で表される細胞内の２つのインスリンが結合しリン酸化により活性化されたインスリン受容体生成速度、下式（１１６）で表される未リン酸化インスリン受容体基質１（ＩＲＳ１）生成速度、下式（１１７）で表されるチロシンリン酸化インスリン受容体基質１生成速度、下式（１１８）で表されるセリンリン酸化インスリン受容体基質１生成速度、下式（１１９）で表されるフォスファチジルイノシトール３キナーゼ（ＰＩ３Ｋ）生成速度、下式（１２０）で表されるセリンリン酸化インスリン受容体基質１とフォスファチジルイノシトール３キナーゼ（ＰＩ３Ｋ）の複合体生成速度、下式（１２１）で表されるフォスファチジルイノシトール（３，４，５）三リン酸生成速度、下式（１２２）で表されるフォスファチジルイノシトール（３，４）二リン酸生成速度、下式（１２３）で表されるフォスファチジルイノシトール（４，５）二リン酸生成速度、下式（１２４）で表される未リン酸ＡＫＴ生成速度、下式（１２５）で表されるリン酸化ＡＫＴ生成速度、下式（１２６）で表される未リン酸ＰＫＣ生成速度、下式（１２７）で表されるリン酸化ＰＫＣ生成速度、下式（１２８）で表される細胞内グルコース輸送担体（ＧＬＵＴ４）生成速度、下式（１２９）で表される細胞膜上グルコース輸送担体（ＧＬＵＴ４）生成速度、下式（１３０）で表される酵素活性補助変数ｘ_２２生成速度及び下式（１３１）で表される酵素活性補助変数ｘ_２３生成速度がそれぞれ演算される。

また、インスリンシグナルブロック２２２は、後述するように、第１グルコース取込速度推定処理と、第２グルコース取込速度推定処理の両方において用いられる。第１グルコース取込速度推定処理では、インスリンシグナルブロック２２２を用いることにより、空腹時におけるグルコースの取込率ＧＬＵＴを演算する処理が行われる。また、第２グルコース取込速度演算処理では、インスリンシグナルブロック２２２を用いることにより、高インスリン血症時におけるグルコースの取込率ＧＬＵＴを演算する処理が行われる。

＜解糖系ブロック＞
解糖系ブロック２２３は、生体器官のグルコース分解機能を仮想的に再現する機能ブロックである。身体のグルコース分解機能は、グルコース輸送担体（ＧＬＵＴ４）の細胞膜上に占める割合に応じてグルコースを細胞内に取り込み、細胞内のグルコースを分解して、グルコース６−リン酸（Ｇ６Ｐ）、Ｄ−グリセルアルデヒド３−リン酸（ＧＡ３Ｐ）、１，３−ビスホスホグリセリン酸（ＢＰＧ）を介してピルビン酸（ＰＹＲ）を産生する。またピルビン酸（ＰＹＲ）より乳酸（ＬＡＣ）を生成する。解糖系ブロック２２３は、かかる身体のグルコース分解機能を模したものであり、当該解糖系ブロック２２３をＣＰＵ２１１ａが実行することにより、インスリンシグナルブロック２２で得られたグルコース取込率ＧＬＵＴ、組織内での酸素消費速度及び二酸化炭素産生速度に基づいてグルコース取込速度が演算され、グルコースがＧ６Ｐに変換される反応速度ｆ_{ＧＬＵ→Ｇ６Ｐ}、Ｇ６ＰがＧＡ３Ｐに変換される反応速度ｆ_{Ｇ６Ｐ→ＧＡ３Ｐ}、ＧＡ３ＰがＢＰＧに変換される反応速度ｆ_{ＧＡ３Ｐ→ＢＰＧ}、ＢＰＧがＰＹＲに変換される反応速度ｆ_{ＢＰＧ→ＰＹＲ}、ＰＹＲがＬＡＣに変換される反応速度ｆ_{ＰＹＲ→ＬＡＣ}、ＬＡＣがＰＹＲに変換される反応速度ｆ_{ＬＡＣ→ＰＹＲ}がそれぞれ演算される。また、これらの反応速度に基づき、一定時間経過後のグルコース（ＧＬＵ）濃度、Ｇ６Ｐ濃度、ＧＡ３Ｐ濃度、ＢＰＧ濃度、ＰＹＲ濃度及びＬＡＣ濃度が演算される。

以下に、生体器官のグルコース分解機能に関連する化学反応と、それに基づいた解糖系ブロック２２３の具体的な演算処理について説明する。なお、以下の式におけるグルコース、Ｇ６Ｐ、ＧＡ３Ｐ、ＢＰＧ、ピルビン酸、ＬＡＣ、ＮＡＤ、及びＮＡＤＨの細胞内濃度は、それぞれ所定の初期値が与えられており、１度目の演算では初期値が用いられ、その後の演算では更新された値が用いられる。

まず、血中のグルコースが、ＧＬＵＴ４の発現量（細胞表面における発現量）に応じて細胞内に取り込まれる。この取り込み速度ｆ_ＧＬＵは、以下の式（１３２）で表される。

ただし、Ｃ_{ＧＬＵ，ｂ}は、入力された血中のグルコース濃度を、Ｃ_ＧＬＵは細胞内のグルコース濃度を、σ_ＧＬＵはグルコースに関する分配係数をそれぞれ示す。

まず、血中の乳酸（ＬＡＣ）が細胞内に取り込まれる。この取り込み速度ｆ_ＬＡＣは、以下の式（１３３）で表される。

ただし、Ｃ_{ＬＡＣ，ｂ}は、入力された血中の乳酸濃度を、Ｃ_ＬＡＣは細胞内の乳酸濃度を、σ_ＬＡＣは乳酸に関する分配係数をそれぞれ示す。

まず、血中のピルビン酸（ＰＹＲ）が細胞内に取り込まれる。この取り込み速度ｆ_ＰＹＲは、以下の式（１３４）で表される。

ただし、Ｃ_{ＰＹＲ，ｂ}は、入力された血中のピルビン酸濃度を、Ｃ_ＰＹＲは細胞内のピルビン酸濃度を、σ_ＰＹＲはピルビン酸に関する分配係数をそれぞれ示す。

細胞内へ取り込まれたグルコース（ＧＬＵ）はヘキソキナーゼによりリン酸化され、Ｇ６Ｐに変換される（式（１７））。この反応の速度は、次式（１３５）で表される。

ただし、ＨＫ（ｉｎｓｕｌｉｎ）はインスリンシグナリングブロック２２２において演算される。

また、変換されて生じたＧ６ＰはＡＴＰと反応して、式（１９）の通りＧＡ３ＰとＡＤＰとを生じる。この反応の速度は、次式（１３６）で示される。

変換されて生じたＧＡ３Ｐは、次式（１３７）に従って１，３−ビスホスホグリセリン酸（ＢＰＧ）に変換される。
ＧＡ３Ｐ＋Ｐｉ＋ＮＡＤ→ＢＰＧ＋ＮＡＤＨ …（１３７）
このＢＰＧに変換される反応速度は、次式（１３８）で示される。

変換されて生じたＢＰＧは次式（１３９）に従ってピルビン酸（ＰＹＲ）に変換される。
ＢＰＧ＋２ＡＤＰ→ＰＹＲ＋２ＡＴＰ …（１３９）
このＰＹＲに変換される反応速度は、次式（１４０）で示される。

生成されたピルビン酸（ＰＹＲ）は、次式（１４１）に従って乳酸（ＬＡＣ）に変換される。
ＰＹＲ＋ＮＡＤＨ→ＬＡＣ＋ＮＡＤ …（１４１）
この反応速度は、次式（１４２）で示される。

乳酸（ＬＡＣ）は、次式（１４３）に従ってピルビン酸（ＰＹＲ）に変換される。
ＬＡＣ＋ＮＡＤ→ＰＹＲ＋ＮＡＤＨ …（１４３）
この反応速度は、次式（１４４）で示される。

解糖系ブロック２２３では、ＣＰＵ２１１ａにより、上記の式（１３２），（１３３），（１３４），（１３５），（１３６），（１３８），（１４０），（１４２），（１４４）で表されるｆ_ＧＬＵ，ｆ_ＬＡＣ，ｆ_ＰＹＲ，ｆ_{ＧＬＵ→Ｇ６Ｐ}，ｆ_{Ｇ６Ｐ→ＧＡ３Ｐ}，ｆ_{ＧＡ３Ｐ→ＢＰＧ}，ｆ_{ＢＰＧ→ＰＹＲ}，ｆ_{ＰＹＲ→ＬＡＣ}，ｆ_{ＬＡＣ→ＰＹＲ}の各反応速度が演算される。

また、解糖系ブロック２２３では、ＣＰＵ２１１ａにより、下式（１４５）で表されるグルコースの生成速度、下式（１４６）で表されるＧ６Ｐの生成速度、下式（１４７）で表されるＧＡ３Ｐの生成速度、下式（１４８）で表されるＢＰＧの生成速度、下式（１４９）で表されるピルビン酸の生成速度、下式（１５０）で表される乳酸生成速度がそれぞれ演算される。

ただし、式（１４９）において、反応速度ｆ_{ＰＹＲ→ＡＣｏＡ}は、後述する式（１５１）で表され、ミトコンドリアブロック２２４において算出され、反応速度ｆ_{ＰＹＲ→ＡＬＡ}は、補足演算処理で算出される。ミトコンドリアブロック２２４と補足演算処理の演算が一度も行われていないときには、反応速度ｆ_{ＰＹＲ→ＡＣｏＡ}、ｆ_{ＰＹＲ→ＡＬＡ}の初期値が用いられる。

さらに解糖系ブロック２２３では、ＣＰＵ２１１ａにより、上記のようにして得られたグルコース、Ｇ６Ｐ、ＧＡ３Ｐ、ＢＰＧ、ピルビン酸及び乳酸の生成速度から、一定時間に生成されるグルコース、Ｇ６Ｐ、ＧＡ３Ｐ、ＢＰＧ、ピルビン酸及び乳酸の量がそれぞれ演算され、その時点での細胞内におけるグルコース濃度、Ｇ６Ｐ濃度、ＧＡ３Ｐ濃度、ＢＰＧ濃度、ピルビン酸濃度にこれらの量を反映させることで、前記一定時間経過後の細胞内のグルコース濃度、Ｇ６Ｐ濃度、ＧＡ３Ｐ濃度、ＢＰＧ濃度及びピルビン酸濃度が演算される。

また、解糖系ブロック２２３は、後述するように、第１グルコース取込速度推定処理と、第２グルコース取込速度推定処理の両方において用いられる。第１グルコース取込速度推定処理では、解糖系ブロック２２３を用いることにより、空腹時における細胞内のグルコース濃度、Ｇ６Ｐ濃度、ＧＡ３Ｐ濃度、ＢＰＧ濃度、ピルビン酸濃度、及び乳酸濃度を演算する処理が行われる。また、第２グルコース取込速度演算処理では、解糖系ブロック２２３を用いることにより、高インスリン血症時における前記各物質の細胞内濃度を演算する処理が行われる。

＜ミトコンドリアブロック＞
ミトコンドリアブロック２２４は、生体器官のミトコンドリアの機能を仮想的に再現する機能ブロックである。ミトコンドリアは、ピルビン酸（ＰＹＲ）及び脂肪酸アシル−補酵素Ａ複合体（ＦＡＣ）をそれぞれ酸化させることでアセチル補酵素Ａ（ＡＣｏＡ）を生成し、このアセチル補酵素Ａ（ＡＣｏＡ）をＴＣＡ回路で代謝し、Ｈ_２ＯとＣＯ_２に変換する。また、アセチル補酵素Ａ（ＡＣｏＡ）の濃度により、脂肪酸酸化が抑制される。ミトコンドリアブロック２２４は、かかるミトコンドリアの機能を模したものであり、当該ミトコンドリアブロック２２４をＣＰＵ２１１ａが実行することにより、解糖系ブロック２２３で得られたピルビン酸からアセチル補酵素Ａ（ＡＣｏＡ）が生成される反応速度ｆ_{ＰＹＲ→ＡＣｏＡ}、脂肪酸代謝ブロック２１で得られた脂肪酸アシル−補酵素Ａ複合体（ＦＡＣ）からアセチル補酵素Ａ（ＡＣｏＡ）が生成される反応速度ｆ_{ＦＡＣ→ＡＣｏＡ}、アセチル補酵素Ａ（ＡＣｏＡ）とオキサロ酢酸（ＯＸＡ）からクエン酸（ＣＩＴ）が生成され、その後、αケトグルタル酸（ａＫＧ）、スクシニルＣｏＡ（ＳＣｏＡ）、コハク酸（ＳＵＣ）、リンゴ酸（ＭＡＬ）、オキサロ酢酸（ＯＸＡ）の順に代謝・生成される各反応速度ｆ_{ＡＣｏＡ→ＣＩＴ}、ｆ_{ＣＩＴ→ａＫＧ}、ｆ_{ａＫＧ→ＳＣｏＡ}、ｆ_{ＳＣｏＡ→ＳＵＣ}、ｆ_{ＳＵＣ→ＭＡＬ}、ｆ_{ＭＡＬ→ＯＸＡ}、及び酸素が消費され水が生成される反応速度ｆ_{Ｏ２→Ｈ２Ｏ}がそれぞれ演算される。また、これらの反応速度に基づき、一定時間経過後のアセチル補酵素Ａ（ＡＣｏＡ）、クエン酸（ＣＩＴ）、αケトグルタル酸（ａＫＧ）、スクシニルＣｏＡ（ＳＣｏＡ）、コハク酸（ＳＵＣ）、リンゴ酸（ＭＡＬ）及びオキサロ酢酸（ＯＸＡ）、酸素（Ｏ_２）及び二酸化炭素（ＣＯ_２）の細胞内濃度が演算される。

以下に、ミトコンドリアに関連する生化学反応と、それに基づいたミトコンドリアブロック２２４の具体的な演算処理について説明する。なお、以下の式におけるアセチル補酵素Ａ（ＡＣｏＡ）、補酵素Ａ（ＣｏＡ）、酸素（Ｏ_２）、クエン酸（ＣＩＴ）、αケトグルタル酸（ａＫＧ）、スクシニルＣｏＡ（ＳＣｏＡ）、コハク酸（ＳＵＣ）、リンゴ酸（ＭＡＬ）及びオキサロ酢酸（ＯＸＡ）及び二酸化炭素（ＣＯ_２）の細胞内濃度は、それぞれ所定の初期値が与えられており、１度目の演算では初期値が用いられ、その後の演算では更新された値が用いられる。

まず、細胞内のピルビン酸（ＰＹＲ）が酸化され、アセチル補酵素Ａ（ＡＣｏＡ）が生成される（式（２７））。この反応の速度は、次式（１５１）で表される。

ただしＰＤＨ（ｉｎｓｕｌｉｎ）は、インスリンシグナリングブロック２２２において演算される。

また、脂肪酸代謝ブロック２２１で生成された脂肪酸アシル−補酵素Ａ複合体（ＦＡＣ）が酸化され、式（２９）の通りアセチル補酵素Ａ（ＡＣｏＡ）が生成される（β酸化）。この反応の速度は、次式（１５２）で表される。

ただしＨＡＤ（ｉｎｓｕｌｉｎ）はインスリンシグナリングブロック２２２にて演算される。

ピルビン酸（ＰＹＲ）と脂肪酸アシル−補酵素Ａ複合体（ＦＡＣ）より生成されたアセチル−補酵素Ａ複合体（ＡＣｏＡ）は、以下に示す一連の式によって表されるＴＣＡ回路により代謝され、新たにＡＴＰ及びＮＡＤＨが生成される。ＴＣＡ回路内では、オキサロ酢酸（ＯＸＡ）とアセチルＣｏＡからクエン酸（ＣＩＴ）が生成され、その後、αケトグルタル酸（ａＫＧ）、スクシニルＣｏＡ（ＳＣｏＡ）、コハク酸（ＳＵＣ）、リンゴ酸（ＭＡＬ）、オキサロ酢酸（ＯＸＡ）の順に生成され代謝される。

まずＡＣｏＡとＯＸＡからクエン酸シンターゼによりＣＩＴとＣｏＡに変換される（式（１５３））。
ＡＣｏＡ＋ＯＸＡ→ＣＩＴ＋ＣｏＡ …（１５３）
この反応の速度は、次式（１５４）で表される。

クエン酸（ＣＩＴ）が、ｃｉｓ−アコニット酸、イソクエン酸、オキサロコハク酸を経由してαケトグルタル酸（ａＫＧ）が生成される反応過程を簡略化すると次の式で表される。
ＣＩＴ＋ＮＡＤ→ａＫＧ＋ＮＡＤＨ＋ＣＯ_２ …（１５５）
この反応の速度は、次式（１５６）で表される。

αケトグルタル酸（ａＫＧ）は酸化されスクシニルＣｏＡ（ＳＣｏＡ）と二酸化炭素（ＣＯ_２）が生成される（式（１５７））。
ａＫＧ＋ＣｏＡ＋ＮＡＤ→ＳＣｏＡ＋ＮＡＤＨ＋ＣＯ_２ …（１５７）
この反応の速度は、次式（１５８）で表される。

スクシニルＣｏＡ（ＳＣｏＡ）は、スクシニルＣｏＡシンテターゼによってコハク酸（ＳＵＣ）が生成される（式（１５９））。
ＳＣｏＡ＋ＡＤＰ＋Ｐｉ→ＳＵＣ＋ＣｏＡ＋ＡＴＰ …（１５９）
この反応の速度は、次式（１６０）で表される。

コハク酸（ＳＵＣ）からフマル酸を介してリンゴ（ＭＡＬ）が生成される（式（１６１））。

この反応の速度は、次式（１６２）で表される。

リンゴ酸（ＭＡＬ）はリンゴ酸デヒドロゲナーゼによって酸化されオキサロ酢酸（ＯＸＡ）が生成される（式（１６３））。
ＭＡＬ＋ＮＡＤ→ＯＸＡ＋ＮＡＤＨ …（１６３）
この反応の速度は、次式（１６４）で表される。

一方、ミトコンドリアにおけるＮＡＤＨ、酸素、及びＡＤＰの消費、並びにＡＴＰの合成の関係式は以下の式（１６５）ように示される。

この反応の速度は、次式（１６６）で表される。

細胞内への酸素の流入速度ｆ_Ｏ２ は、次式（１６７）のようにも示すことができる。

ただし、Ｃ_Ｏ２，ｂは血中酸素濃度（定数）を、Ｃ_Ｏ２ は細胞内の酸素濃度を、σ_Ｏ２ は酸素に関する分配係数をそれぞれ示す。

同様に、細胞内における二酸化炭素の産生速度ｆ_ＣＯ２は、次式（１６８）のようにも示すことができる。

ただし、Ｃ_{ＣＯ２，ｂ}は血中二酸化炭素濃度（定数）を、Ｃ_ＣＯ２は細胞内の二酸化炭素濃度を、σ_ＣＯ２は二酸化炭素に関する分配係数をそれぞれ示す。

ミトコンドリアブロック２２４では、ＣＰＵ２１１ａにより、上記の式（１５１），（１５２），（１５４），（１５６），（１５８），（１６０），（１６２），（１６４），（１６６），（１６７），（１６８）で表されるｆ_{ＰＹＲ→ＡＣｏＡ}，ｆ_{ＦＡＣ→ＡＣｏＡ}，ｆ_{ＡＣｏＡ→ＣＩＴ}、ｆ_{ＣＩＴ→ａＫＧ}、ｆ_{ａＫＧ→ＳＣｏＡ}、ｆ_{ＳＣｏＡ→ＳＵＣ}、ｆ_{ＳＵＣ→ＭＡＬ}、ｆ_{ＭＡＬ→ＯＸＡ}、ｆ_{Ｏ２→Ｈ２Ｏ}、ｆ_Ｏ２ 、ｆ_ＣＯ２の各反応速度が演算される。

また、ミトコンドリアブロック２２４では、ＣＰＵ２１１ａにより、下式（１６９）で表されるアセチル補酵素Ａ（ＡＣｏＡ）の生成速度、下式（１７０）で表されるクエン酸（ＣＩＴ）の生成速度、下式（１７１）で表されるαケトグルタル酸（ａＫＧ）の生成速度、下式（１７２）で表されるスクシニルＣｏＡ（ＳＣｏＡ）の生成速度、下式（１７３）で表されるコハク酸（ＳＵＣ）の生成速度、下式（１７４）で表されるリンゴ酸（ＭＡＬ）の生成速度、下式（１７５）で表されるオキサロ酢酸（ＯＸＡ）の生成速度、下式（１７６）で表される酸素（Ｏ_２）の生成速度及び下式（１７７）で表される二酸化炭素（ＣＯ_２）の生成速度がそれぞれ演算される。

さらにミトコンドリアブロック２２４では、ＣＰＵ２１１ａにより、上記のようにして得られたアセチル補酵素Ａ（ＡＣｏＡ）、クエン酸（ＣＩＴ）、αケトグルタル酸（ａＫＧ）、スクシニルＣｏＡ（ＳＣｏＡ）、コハク酸（ＳＵＣ）、リンゴ酸（ＭＡＬ）及びオキサロ酢酸（ＯＸＡ）、酸素（Ｏ_２）及び二酸化炭素（ＣＯ_２）の生成速度から、一定時間に生成（消費）されるアセチル補酵素Ａ（ＡＣｏＡ）、クエン酸（ＣＩＴ）、αケトグルタル酸（ａＫＧ）、スクシニルＣｏＡ（ＳＣｏＡ）、コハク酸（ＳＵＣ）、リンゴ酸（ＭＡＬ）及びオキサロ酢酸（ＯＸＡ）、酸素（Ｏ_２）及び二酸化炭素（ＣＯ_２）の量がそれぞれ演算され、前記一定時間経過後の細胞内のアセチル補酵素Ａ（ＡＣｏＡ）、クエン酸（ＣＩＴ）、αケトグルタル酸（ａＫＧ）、スクシニルＣｏＡ（ＳＣｏＡ）、コハク酸（ＳＵＣ）、リンゴ酸（ＭＡＬ）及びオキサロ酢酸（ＯＸＡ）、酸素（Ｏ_２）及び二酸化炭素（ＣＯ_２）が演算される。

また、ミトコンドリアブロック２２４は、後述するように、第１グルコース取込速度推定処理と、第２グルコース取込速度推定処理の両方において用いられる。第１グルコース取込速度推定処理では、ミトコンドリアブロック２２４を用いることにより、空腹時における細胞内のアセチル補酵素Ａ（ＡＣｏＡ）、クエン酸（ＣＩＴ）、αケトグルタル酸（ａＫＧ）、スクシニルＣｏＡ（ＳＣｏＡ）、コハク酸（ＳＵＣ）、リンゴ酸（ＭＡＬ）及びオキサロ酢酸（ＯＸＡ）、酸素（Ｏ_２）及び二酸化炭素（ＣＯ_２）を演算する処理が行われる。また、第２グルコース取込速度演算処理では、ミトコンドリアブロック２２４を用いることにより、高インスリン血症時における前記各物質の細胞内濃度を演算する処理が行われる。

＜補足演算処理＞
また、ＣＰＵ２１１ａにより、以下のような補足演算処理が行われる。

ＡＤＰの濃度が高いときは、クレアチンリン酸（ＰＣｒ）からリン酸基が外れクレアチン（Ｃｒ）が生成されるとともに、ＡＤＰからＡＴＰが生成される（式（３９））。この反応の速度は、次式（１７８）で表される。

一方、ＡＴＰの濃度が高くなると、クレアチンとＡＴＰとによって、クレアチンリン酸とＡＤＰとが生成される（式（４１））。この反応の速度は、次式（１７９）で表される。

また、ＡＴＰは加水分解され、式（４３）に従ってＡＤＰに変換される。この反応の速度は、次式（１８０）で表される。

また、ＡＭＰとＡＴＰとはアデニル酸キナーゼによってＡＤＰが産生される（式（１８１））。
ＡＭＰ＋ＡＴＰ→２ＡＤＰ …（１８１）
この反応の速度は、次式（１８２）で表される。

２つのＡＤＰから、アデニル酸キナーゼによりＡＭＰとＡＴＰが生成される（式（１８３））。
２ＡＤＰ→ＡＭＰ＋ＡＴＰ …（１８３）
この反応の速度は、次式（１８４）で表される。

解糖系ブロック２２３のＧ６ＰとＡＴＰから複数の反応を経由し、グリコーゲン（ＧＬＹ）が生成される反応を簡略化すると反応式は次式のようになる。
Ｇ６Ｐ＋ＡＴＰ→ＧＬＹ＋ＡＤＰ＋２Ｐｉ …（１８５）
この反応の速度は、次式（１８６）で表される。

ただし、ＧＳ（ｉｎｓｕｌｉｎ）はインスリンシグナリングブロック２２２にて演算される。

グリコーゲン（ＧＬＹ）が複数の反応を経由して分解され解糖系ブロック２２３のＧ６Ｐが生成される反応を簡略化すると反応式は次式のようになる。
ＧＬＹ→Ｇ６Ｐ＋Ｐｉ …（１８７）
この反応の速度は、次式（１８８）で表される。

解糖系ブロック２２３で生成されたピルビン酸（ＰＹＲ）とグルタミン酸からアラニン（ＡＬＡ）が生成される反応を簡略化すると反応式は次式のようになる。
ＰＹＲ→ＡＬＡ …（１８９）
この反応の速度は、次式（１９０）で表される。

血中のアラニン（ＡＬＡ）が細胞内に取り込まれる。この取り込み速度ｆ_ＡＬＡは以下の式（１９１）で表される。

ただし、Ｃ_{ＡＬＡ，ｂ}は入力された血中の脂肪酸濃度を、Ｃ_ＡＬＡは細胞内の脂肪酸濃度を、σ_ＡＬＡは脂肪酸に関する分配係数をそれぞれ示す。

補足演算処理では、ＣＰＵ２１１ａにより、上記の式（１７８），（１７９），（１８０），（１８２），（１８４），（１８６），（１８８），（１９０），（１９１）で表されるｆ_{ＰＣｒ→Ｃｒ}，ｆ_{Ｃｒ→ＰＣｒ}，ｆ_{ＡＴＰ→ＡＤＰ}，ｆ_{ＡＭＰ→ＡＤＰ}，ｆ_{ＡＤＰ→ＡＭＰ}，ｆ_{Ｇ６Ｐ→ＧＬＹ}，ｆ_{ＧＬＹ→Ｇ６Ｐ}，ｆ_{ＰＹＲ→ＡＬＡ}，ｆ_ＡＬＡの各反応速度が演算される。

また、ＣＰＵ２１１ａにより、以下の式（１９２）〜（２０２）に従って、ＧＬＹ、ＡＬＡ、ＮＡＤ、ＮＡＤＨ、ＡＴＰ、ＡＤＰ、Ｐｉ、ＰＣｒ、Ｃｒ、ＣｏＡの各物質の生成速度が演算される。

さらに、ＣＰＵ２１１ａにより、上記のようにして得られたＧＬＹ、ＡＬＡ、ＮＡＤ、ＮＡＤＨ、ＡＴＰ、ＡＤＰ、ＡＭＰ、Ｐｉ、ＰＣｒ、Ｃｒ、ＣｏＡの各生成速度から、一定時間に生成されるＧＬＹ、ＡＬＡ、ＮＡＤ、ＮＡＤＨ、ＡＴＰ、ＡＤＰ、ＡＭＰ、Ｐｉ、ＰＣｒ、Ｃｒ、ＣｏＡの量がそれぞれ演算され、その時点での細胞内におけるＧＬＹ、ＡＬＡ、ＮＡＤ、ＮＡＤＨ、ＡＴＰ、ＡＤＰ、ＡＭＰ、Ｐｉ、ＰＣｒ、Ｃｒ、ＣｏＡの濃度にこれらの量を反映させることで、前記一定時間経過後の細胞内のＧＬＹ、ＡＬＡ、ＮＡＤ、ＮＡＤＨ、ＡＴＰ、ＡＤＰ、ＡＭＰ、Ｐｉ、ＰＣｒ、Ｃｒ、ＣｏＡの各濃度が演算される。

また、補足演算処理は、後述するように、第１グルコース取込速度推定処理と、第２グルコース取込速度推定処理の両方において用いられる。第１グルコース取込速度推定処理では、補足演算処理を実行することにより、空腹時における細胞内のＧＬＹ、ＡＬＡ、ＮＡＤ、ＮＡＤＨ、ＡＴＰ、ＡＤＰ、ＡＭＰ、Ｐｉ、ＰＣｒ、Ｃｒ、ＣｏＡの各濃度を演算する処理が行われる。また、第２グルコース取込速度演算処理では、補足演算処理を実行することにより、高インスリン血症時における前記各物質の細胞内濃度を演算する処理が行われる。

上記の本実施の形態に係る脂肪酸代謝ブロック２２１、インスリンシグナリングブロック２２２、解糖系ブロック２２３、及びミトコンドリアブロック２２４については、以下の文献を参考にして作成されたものである。
1. A computational model of skeletal muscle metabolism linking cellular adaptations induced by altered loading states to metabolic responses during exercise.
Dash RK, Dibella JA 2nd, Cabrera ME.
Biomed Eng Online. 2007 Apr 20;6:14. PMID: 17448235
2. Metabolic dynamics in skeletal muscle during acute reduction in blood flow and oxygen supply to mitochondria: in-silico studies using a multi-scale, top-down integrated model.
Dash RK, Li Y, Kim J, Beard DA, Saidel GM, Cabrera ME.
PLoS ONE. 2008 Sep 9;3(9):e3168. PMID: 18779864
3. A mathematical model of metabolic insulin signaling pathways.
Sedaghat AR, Sherman A, Quon MJ.
Am J Physiol Endocrinol Metab. 2002 Nov;283(5):E1084-101. PMID: 12376338

［インスリン抵抗性評価支援システムの動作］
次に、本実施の形態に係るインスリン抵抗性評価支援システム２０１の動作について説明する。図１１は、本実施の形態に係るインスリン抵抗性評価支援プログラムの処理の流れを示すフローチャートである。コンピュータ２０１ａは、インスリン抵抗性評価支援プログラム２１４ａを実行することにより、以下のように動作する。まず、被験者は、予め体重計で体重を計測し、体組成計により骨格筋率の測定、血液検査を行うことにより、被験者の体重、骨格筋率、血中グルコース濃度、血漿インスリン濃度及び血中遊離脂肪酸濃度の各生体情報をそれぞれ取得しておく。なお、血中グルコース濃度、血漿インスリン濃度、及び血中遊離脂肪酸濃度は、それぞれ空腹時における測定値が用いられる。これらの生体情報は、インスリン抵抗性評価支援システムを使用する医師、オペレータ等のユーザに予め連絡される。

インスリン抵抗性評価支援プログラム２１４ａが起動された後、まず、ＣＰＵ２１１ａは、被験者の体重、骨格筋率、血中グルコース濃度、血漿インスリン濃度及び血中遊離脂肪酸濃度の入力をユーザに促すための入力画面を表示する（ステップＳ２１）。図１２は、実施の形態２に係るインスリン抵抗性評価支援システム２０１の入力画面の一例を示す模式図である。図１２に示すように、入力画面２３０には、被験者の体重、骨格筋率、血中グルコース濃度、血漿インスリン濃度及び血中遊離脂肪酸濃度をそれぞれ入力するための入力エリア２３１，２３４〜２３７と、入力エリア２３１，２３４〜２３７に各生体情報を入力した後に、インスリン抵抗性の推定の実行を指示するための実行ボタン２３８とが設けられている。ユーザは、入力部２１３を操作することによりこれらの入力エリア２３１，２３４〜２３７に被験者の体重、骨格筋率、血中グルコース濃度、血漿インスリン濃度及び血中遊離脂肪酸濃度の各生体情報を入力し、実行ボタン２３８を選択（クリック）することによりインスリン抵抗性の推定の実行を指示する。ＣＰＵ２１１ａは、かかるユーザからの生体情報及び実行指示の入力を受け付ける（ステップＳ２２）。かかる生体情報の入力があった場合、ＣＰＵ２１１ａに割り込みが発生し、ステップＳ２３以下の処理が呼び出される。

ユーザから生体情報及び実行指示の入力を受け付けると、ＣＰＵ２１１ａは、第１グルコース取込速度推定処理（ステップＳ２３）を実行する。図１３は、実施の形態２に係る第１グルコース取込速度推定処理の手順を示すフローチャートである。この第１グルコース取込速度推定処理は、空腹時におけるグルコース取込速度を推定する処理であり、ステップＳ２３１〜Ｓ２３８の処理を順次的に繰り返し実行する処理となっている。まず、ＣＰＵ２１１ａは、変数に初期値をセット（初期化）する（ステップＳ２３１）。この処理では、上述した各物質の空腹時における細胞内濃度及び各反応速度の初期値をＲＡＭ２１１ｃに記憶し、また、グルコース取込速度の値として、初期値０をＲＡＭ２１１ｃに記憶する。次に、ＣＰＵ２１１ａは、その時点でのグルコース取込速度の値を、前回のグルコース取込速度の値としてＲＡＭ２１１ｃに記憶する（ステップＳ２３２）。そして、ＣＰＵ２１１ａは、脂肪酸代謝ブロック２２１における演算処理である脂肪酸代謝ブロック演算処理（ステップＳ２３３）、インスリンシグナルブロック２２２における演算処理であるインスリンシグナルブロック演算処理（ステップＳ２３４）、解糖系ブロック２２３における演算処理である解糖系ブロック演算処理（ステップＳ２３５）、ミトコンドリアブロック２２４における演算処理であるミトコンドリアブロック演算処理（ステップＳ２３６）、及び補足演算処理（ステップＳ２３７）を順次実行する。かかるステップＳ２３３〜２３７において得られる物質の濃度、反応速度、ＧＬＵＴ４発現量に応じた数値及びグルコース取込速度は、空腹時における数値となる。

次にＣＰＵ２１１ａは、上記の処理によって得られた空腹時におけるグルコース取込速度が定常状態に到達したか否かを判定する（ステップＳ２３８）。この処理は、本実施の形態においては、今回の演算（ターン）によって求められたグルコース取込速度の値と、ＲＡＭ２１１ｃに記憶された前回の演算（ターン）のグルコース取込速度の値との差を求め、その差がグルコース取込速度が定常状態に至ったか否かを判定するための第１基準値未満であるか否かを判定することにより行われる。そして、ＣＰＵ２１１ａは、グルコース取込速度の今回値と前回値との差が第１基準値未満であるときは（ステップＳ２３８においてＹＥＳ）、メインルーチンにおける第１グルコース取込速度推定処理の呼出アドレスへ処理をリターンし、差が第１基準値以上であるときは（ステップＳ２３８においてＮＯ）、再度ステップＳ２３２以下の処理を繰り返す。

次に、ＣＰＵ２１１ａは、第２グルコース取込速度推定処理を実行する（ステップＳ２４）。図１４は、実施の形態２に係る第２グルコース取込速度推定処理の手順を示すフローチャートである。この第２グルコース取込速度推定処理は、高インスリン血症時におけるグルコース取込速度を推定する処理であり、ステップＳ２４１〜Ｓ２４８の処理を順次的に繰り返し実行する処理となっている。まず、ＣＰＵ２１１ａは、インスリン濃度に所定値（高インスリン血症時のインスリン濃度）をセットする（ステップＳ２４１）。なお、第２グルコース取込速度推定処理では、インスリン濃度以外の各物質の空腹時における細胞内濃度及び各反応速度並びにグルコース取込率（ＧＬＵＴ）の初期値として、第１グルコース取込速度推定処理において最終的に得られた数値、すなわち、グルコース取込速度が定常状態に達したときの各物質の濃度及び反応速度並びにグルコース取込率がそのまま用いられる。血中グルコース濃度については、第１グルコース取込速度推定処理で使用した空腹時のグルコース濃度がそのまま用いられる。また、グルコース取込速度の値として、初期値０がＲＡＭ２１１ｃに記憶される。次に、ＣＰＵ２１１ａは、その時点でのグルコース取込速度の値を、前回のグルコース取込速度の値としてＲＡＭ２１１ｃに記憶する（ステップＳ２４２）。そして、ＣＰＵ２１１ａは、脂肪酸代謝ブロック２２１における演算処理である脂肪酸代謝ブロック演算処理（ステップＳ２４３）、インスリンシグナルブロック２２２における演算処理であるインスリンシグナルブロック演算処理（ステップＳ２４４）、解糖系ブロック２２３における演算処理である解糖系ブロック演算処理（ステップＳ２４５）、ミトコンドリアブロック２２４における演算処理であるミトコンドリアブロック演算処理（ステップＳ２４６）、及び補足演算処理（ステップＳ２４７）を順次実行する。かかるステップＳ２４３〜２４７において得られる物質の濃度、反応速度、ＧＬＵＴ４発現量に応じた数値及びグルコース取込速度は、高インスリン血症時における数値となる。

次にＣＰＵ２１１ａは、上記の処理によって得られた高インスリン血症時におけるグルコース取込速度が定常状態に到達したか否かを判定する（ステップＳ２４８）。この処理は、本実施の形態においては、今回の演算（ターン）によって求められたグルコース取込速度の値と、ＲＡＭ２１１ｃに記憶された前回の演算（ターン）のグルコース取込速度の値との差を求め、その差がグルコース取込速度が定常状態に至ったか否かを判定するための第１基準値未満であるか否かを判定することにより行われる。なお、第１グルコース取込速度推定処理と第２グルコース取込速度推定処理とにおいて、同じ第１の基準値を用いてグルコース取込速度が定常状態に至ったか否かを判定する構成としたが、これに限られず、それぞれで異なる基準値を用いてもよい。そして、ＣＰＵ２１１ａは、グルコース取込速度の今回値と前回値との差が第１基準値未満であるときは（ステップＳ２４８においてＹＥＳ）、メインルーチンにおける第２グルコース取込速度推定処理の呼出アドレスへ処理をリターンし、差が第１基準値以上であるときは（ステップＳ２４８においてＮＯ）、再度ステップＳ２４２以下の処理を繰り返す。

次に、ＣＰＵ２１１ａは、インスリン抵抗性の有無を推定する（ステップＳ２５）。この処理では、ステップＳ２４８で得た単位筋肉量当たりのグルコース取込速度の値（推定値）が、インスリン抵抗性の有無を推定するための第２基準値（人種によって異なるが、例えば、グルコースクランプ時のインスリン濃度が約３４８０ｐＭの場合に、１２．０ｍｇ／ｋｇ／ｍｉｎ）以上か否かを判定することにより行われる。これにより、グルコース取込速度の推定値が前記第２基準値以上であれば、インスリン抵抗性がない、即ち、インスリン感受性がある、と推定することができ、グルコース取込速度の推定値が前記第２基準値未満であれば、インスリン抵抗性がある、即ち、インスリン感受性がない、と推定することができる。このようにして、ＣＰＵ２１１ａは、インスリン抵抗性の推定を行う。なお、第２基準値を１２．０ｍｇ／ｋｇ／ｍｉｎとする根拠は、次の文献の記載による。
Critical evaluation of adult treatment panel III criteria in identifying insulin resistance with dyslipidemia.
Liao Y, Kwon S, Shaughnessy S, Wallace P, Hutto A, Jenkins AJ, Klein RL, Garvey WT.
Diabetes Care. 2004 Apr;27(4):978-83.PMID: 15047659

次にＣＰＵ１１ａは、インスリン抵抗性の推定結果の出力画面を表示する（ステップＳ２６）。かかる出力画面には、上記の処理によって得られたインスリン抵抗性の推定結果と、最終的に得られた単位筋肉量当たりのグルコース取込速度の推定値とが含まれる。この出力画面により、ユーザに対して、インスリン抵抗性の推定結果及びグルコース取込速度の推定値が通知される。このようにインスリン抵抗性の推定結果及びグルコース取込速度の推定値をユーザに提供することにより、ユーザはこれらの情報を活用してインスリン抵抗性の評価を行うことができる。また、この画面からは、後述するようなシミュレーション結果画面へと表示画面を遷移させることが可能である。シミュレーション結果画面の構成については後述する。

上記のような構成とすることにより、被験者の負担の大きいグルコースクランプ及び経口糖負荷試験の検査結果を必要とせず、被験者の負担が少なく、簡便な検査によって得ることができる「体重」、「骨格筋率」、「空腹時の血中グルコース濃度」、「空腹時の血漿インスリン濃度」、及び「空腹時の血中遊離脂肪酸濃度」の各生体情報によってインスリン抵抗性の有無を推定することができる。このような入力情報の中でも、空腹時の血中グルコース濃度及び血漿インスリン濃度については、上記のように血液検査によって簡便に得られる検査値を用いることができるが、グルコースクランプ又は経口糖負荷試験によって得られた検査値を入力情報として用いることもできる。しかし、グルコースクランプ及び経口糖負荷試験のような被験者にとって負担の大きい検査結果を常に必要とする訳ではなく、これらの検査結果でなくとも、簡便な血液検査の結果を利用することができるという点、及び、入力情報に用いる情報を、血液検査、グルコースクランプ、及び経口糖負荷試験のいずれからも選択することができるという点で、本システムは有用である。

なお、上記実施の形態１及び２においては、第２グルコース取込速度推定処理において、グルコースクランプ試験時の被験者の身体状態を再現するために、高インスリン血症時におけるインスリン濃度をセットし、グルコース濃度は第１グルコース取込速度推定処理から変更せず、空腹時におけるグルコース濃度をそのまま使用する構成としたが、これに限定されるものではない。インスリン抵抗性を評価するグルコースクランプ試験では、一般に正常血糖値かつ高インスリン状態で実施される。グルコースクランプ試験では、空腹時の血中グルコース濃度が維持されるようにインスリン及びグルコースが被験者に投与される場合、及び１００ｍｇ／ｄＬの血中グルコース濃度が維持されるようにインスリン及びグルコースが被験者に投与される場合等があり、後者を再現する場合には、第２グルコース取込速度推定処理において、高インスリン血症時におけるインスリン濃度をセットするだけでなく、血中グルコース濃度として１００ｍｇ／ｄＬをセットすればよい。

また、実施の形態２においては、単位筋肉量当たりのグルコース取込速度を推定し、これを出力する構成としたが、実施の形態１と同様に、単位重量当たりのグルコース取込速度を求め、これを出力する構成としてもよい。

また、実施の形態１及び２においては、グルコース取込速度が定常状態になるまで、第１代謝量演算処理〜補足演算処理を繰り返し実行し、各物質の反応速度、生成速度、及び細胞内濃度を繰り返し更新する構成とした。これにより、グルコース取込速度が定常状態に到達したときには、各物質の細胞内濃度も定常状態に到達していることが期待でき、この定常状態の各物質の濃度は、被験者の身体の状態を反映したものと考えることができる。従って、このような構成とすることにより、被験者の身体の状態を正確に反映したグルコースの取込速度を推定することができる。

［インスリン抵抗性評価支援システムの性能評価実験］
実施の形態２に係るインスリン抵抗性表化支援システム２０１の性能評価実験を行った。グルコースクランプ試験を実施した結果を報告している文献のひとつに、Basuら報告がある(Basu et al., Obesity and Type 2 Diabetes Impair Insulin-Induced Suppression of Glycogenolysis as well as gluconeogenesis, Diabetes 54:1942-1948, 2005)。Basuらの文献では、グルコースクランプ試験を用いて、非糖尿病痩せ型（Lean）１０名, 非糖尿病肥満型（Obese）１０名, ２型糖尿病患者（ＤＭ２）１１名を対象に比較し、主に肝臓インスリン抵抗性について主に論じている。本実験では、文献中のデータを使用し、ＧＩＲ値としてglucose appearance rateとendogeneous glucose productionの差を用いた。

以下の表１は、上記文献中に記載された、Ｌｅａｎ、Ｏｂｅｓｅ、ＤＭ２の測定データをまとめたものである。

本実験では、まず、上記表１に記載したＬｅａｎ、Ｏｂｅｓｅ、ＤＭ２の各測定データの空腹時の血糖値・インスリン値・遊離脂肪酸濃度をインスリン抵抗性評価支援システム２０１に入力し、第１グルコース取込速度推定処理の手順に従って定常状態の各物質濃度と反応速度を求めた。第１グルコース取り込み速度推定処理に用いた初期値と当該推定処理を行ったあとに定常状態となったときに得られた各変数（濃度、反応速度等の変数）の値（以下、「推定値」という）を表２〜４に示す。表２は、インスリン抵抗性評価支援システム２０１における物質の濃度の初期値及び推定値を示し、表３は、インスリンシグナリングブロック２２２におけるパラメータ（物質濃度及び反応速度以外の変数）の初期値及び推定値を示し、表４は、インスリン抵抗性評価支援システム２０１における物質の反応速度の初期値及び推定値を示している。

表２〜４で示された第１グルコース取込速度推定処理後の各変数の推定値を初期値として、第２グルコース取込速度推定手段によりグルコースクランプ試験シミュレーションを行った。このグルコースクランプ試験シミュレーションでは、クランプ時（インスリン注入前の６０分間）における血中グルコース、血中インスリン、及び血中脂肪酸の各濃度に表１の推定値を使用し、インスリン注入前６０分間、注入後２４０分間の計３００分間をシミュレートした。図１５は、本評価実験におけるシミュレーション結果表示画面を示す図である。図に示すように、シミュレーション結果表示画面４５０には、シミュレーションに用いられる各種パラメータの設定値が表示されるパラメータ表示領域４５１と、シミュレーション結果のグラフ４５２〜４５５が含まれる。グラフ４５２は、血中グルコース濃度の時間的変化を示し、グラフ４５３は、血中インスリン濃度の時間的変化を示し、グラフ４５４は、血中脂肪酸濃度の時間的変化を示す。グラフ４５５は、グルコース取込速度の時間的変化を示し、グラフ４５６は、血中脂肪酸の酸化速度の時間的変化を示している。

図１６〜図１９に、さらに詳しくシミュレーションの結果を示す。図１６は、血中グルコース濃度のシミュレーション結果を示すグラフであり、図１７は、血中インスリン濃度のシミュレーション結果を示すグラフであり、図１８は、血中脂肪酸濃度のシミュレーション結果を示すグラフであり、図１９は、グルコース取込速度のシミュレーション結果を示すグラフである。

図１６に示すように、Ｌｅａｎ及びＯｂｅｓｅの血中グルコース濃度は、シミュレーション期間全体を通じて約５ｍＭで概ね一定であった。一方、ＤＭ２の血中グルコース濃度は、インスリン注入前６０分間は約１０ｍＭで一定であり、インスリン注入後に減少した。図１７に示すように、インスリン注入前の６０分間において、Ｌｅａｎの血中インスリン濃度は約２０ｐＭで一定であり、同期間におけるＯｂｅｓｅ及びＤＭ２の血中インスリン濃度は約４５〜５０ｐＭで一定であった。また、インスリン注入後の期間における血中インスリン濃度は、Ｌｅａｎ、Ｏｂｅｓｅ、及びＤＭ２共に１４０ｐＭ前後の一定値を示した。図１８に示すように、インスリン注入前６０分間におけるＬｅａｎ、Ｏｂｅｓｅ、及びＤＭ２の血中脂肪酸濃度は、約０．３〜０．４ｍＭで一定値を示した。また、インスリン注入後におけるＬｅａｎ、Ｏｂｅｓｅ、及びＤＭ２は、いずれも減少し、時間経過と共にその減少率が低下した。図１９に示すように、Ｌｅａｎのグルコース取込速度は、インスリン注入前６０分間において約１８ｕｍｏｌ／ｋｇ／ｍｉｎの一定値を示し、インスリン注入後増加した。また、その増加率は時間経過と共に減少した。Ｏｂｅｓｅ及びＤＭ２のグルコース取込速度は、インスリン注入前６０分間において約２５ｕｍｏｌ／ｋｇ／ｍｉｎの一定値を示し、インスリン注入後増加した。上記のシミュレーションの結果、インスリン注入後２４０分後におけるＬｅａｎ、Ｏｂｅｓｅ、ＤＭ２のそれぞれのグルコース取込速度は、６３．３ｕｍｏｌ／ｋｇ／ｍｉｎ、５３．１ｕｍｏｌ／ｋｇ／ｍｉｎ、４９．８２ｕｍｏｌ／ｋｇ／ｍｉｎとなった。

図２０は、文献中の実測値とシミュレーションの結果との比較結果を示すグラフである。図２０において、縦軸は文献中に記載されたグルコース取込速度を示し、横軸はインスリン抵抗性評価支援システム２０１によって推定したグルコース取込速度を示している。図に示すように、文献に記載されたＬｅａｎのグルコース取込速度の実測値及びインスリン抵抗性評価支援システム２０１によるＬｅａｎのグルコース取込速度の推定値が対応する点２６１と、文献に記載されたＯｂｅｓｅのグルコース取込速度の実測値と、インスリン抵抗性評価支援システム２０１によるＯｂｅｓｅのグルコース取込速度の推定値とが対応する点２６２と、文献に記載されたＤＭ２のグルコース取込速度の実測値と、インスリン抵抗性評価支援システム２０１によるＤＭ２のグルコース取込速度の推定値とが対応する点２６３とは、直線上に並んでいる。このことから、Ｌｅａｎ、Ｏｂｅｓｅ、及びＤＭ２のそれぞれのグルコース取込速度の実測値の大小関係と推定値の大小関係は矛盾しておらず、十分に対応しており、インスリン抵抗性評価支援システム２０１が正確にグルコース取込速度を推定していることがわかる。また、より正確にグルコース取込速度を推定した場合には回帰直線が原点を通るはずである。これは上記の式中の係数等を調整することによって達成できる。

（その他の実施の形態）
上記の実施の形態１及び２においては、グルコース取込速度を推定し、その推定結果に基づいてインスリン抵抗性の有無の推定を行い、インスリン抵抗性の推定結果とグルコース取込速度の推定値とを出力画面に表示する構成としたが、これに限定されるものではない。他の実施の形態として、「インスリン抵抗性があると推定されます。」又は「インスリン抵抗性がないと推定されます。」のように、インスリン抵抗性の推定結果のみを示す画面を表示する構成としてもよい。また、インスリン抵抗性の有無を推定するステップＳ７の処理を行わず、高インスリン血症時におけるグルコース取込速度の演算値が定常状態に至ったときに、単位重量当たりのグルコース取込速度の推定値を示す出力画面を表示する構成としてもよい。この場合、「１２ｍｇ／ｋｇ／ｍｉｎ以上であれば、インスリン抵抗性がないと推定されます。」のような、インスリン抵抗性の評価を支援する情報を同時に表示してもよい。

また、上記の実施の形態１においては、測定された体重と骨格筋率とを入力情報に含め、これらの体重及び骨格筋率から筋肉量を求め、この筋肉量により単位重量当たりのグルコース取込速度を求める構成としたが、これに限定されるものではない。他の実施の形態として、測定された筋肉量を入力情報に含め、その筋肉量から単位重量当たりのグルコース取込速度を求める構成としてもよい。

また、上記の実施の形態１においては、マロニルＣｏＡの濃度を用いて、β酸化におけるアセチルＣｏＡの生成速度を調整する構成としたが、これに限定されるものではない。他の実施の形態として、マロニルＣｏＡによる抑制を考慮せずに、β酸化におけるアセチルＣｏＡの生成速度を求める構成としてもよい。

また、上記の実施の形態１及び２においては、遊離脂肪酸濃度を入力情報に含める構成としたが、これに限定されるものではない。他の実施の形態として、遊離脂肪酸濃度に代えて中性脂肪濃度を入力情報に含め、この中性脂肪濃度から遊離脂肪酸濃度を算出する構成としてもよい。

また、上記の実施の形態１においては、骨格筋の単位時間あたりの酸素消費量及び二酸化炭素産生量並びに骨格筋率を入力情報に含める構成としたが、これに限定されるものではない。他の実施の形態として、全身の単位時間あたりの酸素消費量及び二酸化炭素産生量並びに体脂肪率を入力情報としてもよい。この場合には、細胞内の酸素の消費速度を、全身の単位時間あたりの酸素消費量と体脂肪率とによって求め、細胞内の二酸化炭素の産生量を、全身の単位時間あたりの二酸化炭素産生量と体脂肪率とによって求める構成とする。

また、上記の実施の形態１及び２においては、ＡＴＰの加水分解反応速度がインスリン濃度には依存しないものとしているが、実際にはＡＴＰの加水分解反応を含むグリコーゲンの合成及びナトリウム−カリウムポンプがインスリン濃度に依存して活性化されることが知られている。従って、上記の構成に限るものではなく、ＡＴＰの加水分解反応速度がインスリン濃度によって変化するようにしてもよい。

また、上記の実施の形態１及び２においては、ジアシルグリセロールの濃度は考慮せず、インスリン量及び脂肪酸アシル−補酵素Ａ複合体の濃度を用いてＧＬＵＴ４の発現量に応じた数値を演算する構成としたが、これに限定されるものではない。ＧＬＵＴ４の発現量には、インスリン量及び脂肪酸アシル−補酵素Ａ複合体の濃度に加えて、ジアシルグリセロールの濃度が影響することが分かっている。そこで、他の実施の形態として、ＧＬＵＴ４の発現量に応じた数値ＧＬＵＴを、ジアシルグリセロールの濃度も利用して求める構成とすることもできる。

また、上述した実施の形態１及び２においては、１つのコンピュータ１ａ，２０１ａのＣＰＵ１１ａ，２１１ａにインスリン抵抗性評価支援プログラム１４ａ，２１４ａを実行させることにより、このコンピュータ１ａ，２０１ａをインスリン抵抗性評価支援システム１，２０１として機能させる構成について述べたが、これに限定されるものではなく、インスリン抵抗性評価支援プログラム１４ａ，２１４ａと実質的に同一の処理を実行するための専用のハードウェア回路によりインスリン抵抗性評価支援システムを構成することもできる。

また、上述した実施の形態１及び２においては、単一のコンピュータ１ａ，２０１ａによりインスリン抵抗性評価支援プログラムの全ての処理を実行する構成について述べたが、これに限定されるものではなく、上述したインスリン抵抗性評価支援プログラムと同様の処理を、複数の装置（コンピュータ）により分散して実行する分散システムとすることも可能である。

本発明のインスリン抵抗性評価支援システム、インスリン抵抗性評価支援方法、及びコンピュータプログラムは、被験者のインスリン抵抗性の評価を支援するためのインスリン抵抗性評価支援システム、インスリン抵抗性評価支援方法、及びコンピュータプログラムなどとして有用である。

１，２０１ インスリン抵抗性評価支援システム
１ａ，２０１ａ コンピュータ
１１，２１１ 本体
１２，２１２ 画像表示部
１３，２１３ 入力部
１１ｄ，２１１ｄ ハードディスク
１１ｇ，２１１ｇ 画像出力インタフェース
１１ｅ，２１１ｅ 読出装置
１１ｆ，２１１ｆ 入出力インタフェース
１４，２１４ 可搬型記録媒体
１４ａ，２１４ａ インスリン抵抗性評価支援プログラム
２１，２２１ 脂肪酸代謝ブロック
２２，２２２ インスリンシグナルブロック
２３，２２３ 解糖系ブロック
２４，２２４ ミトコンドリアブロック
３，２３０ 入力画面
３１〜３７，２３１，２３４〜２３７ 入力エリア
３８，２３８ 実行ボタン
４ 出力画面
４１ インスリン抵抗性の推定結果
４２ グルコース取込速度の推定値

Claims (14)

1. 被験者のインスリン抵抗性の評価を支援するためのインスリン抵抗性評価支援システムであって、
   被験者を測定することにより得られた血中のグルコース濃度、インスリン濃度、及び遊離脂肪酸濃度に関する情報の入力を受け付ける入力手段と、
   前記入力手段によって入力を受け付けた血中のグルコース濃度、インスリン濃度、及び遊離脂肪酸濃度に関する情報に基づいて、被験者のグルコース取込速度を推定する推定手段と、
   を備える、インスリン抵抗性評価支援システム。
2. 前記推定手段は、
   血中のインスリン濃度及び遊離脂肪酸濃度に関する情報に基づいて、グルコーストランスポータの発現量に応じた数値を推定する第１推定手段と、
   前記第１推定手段によって推定されたグルコーストランスポータの発現量に応じた数値と血中のグルコース濃度とに基づいて、グルコース取込速度を推定する第２推定手段と、
   を備える、請求項１に記載のインスリン抵抗性評価支援システム。
3. 前記第１推定手段は、
   血中の遊離脂肪酸濃度に関する情報に基づいて、脂肪酸アシル−補酵素Ａ複合体の濃度を推定する第３推定手段と、
   血中のインスリン濃度及び前記第３推定手段によって推定された脂肪酸アシル−補酵素Ａ複合体の濃度に基づいて、グルコーストランスポータの発現量に応じた数値を推定する第４推定手段と、
   を備える、請求項２に記載のインスリン抵抗性評価支援システム。
4. 前記推定手段は、
   前記第２推定手段によって推定されたグルコース取込速度に基づいて、細胞内のピルビン酸の濃度を推定する第５推定手段と、
   前記第３推定手段によって推定された脂肪酸アシル−補酵素Ａ複合体の濃度と、前記第５推定手段によって推定されたピルビン酸の濃度とに基づいて、細胞内のアセチル補酵素の濃度を推定する第６推定手段と、
   前記第６推定手段によって推定されたアセチル補酵素の濃度に基づいて、脂肪酸アシル−補酵素Ａ複合体から生成されるアセチル補酵素の生成速度を調整する調整手段と、
   を更に備え、
   前記第６推定手段は、前記調整手段による調整後のアセチル補酵素の生成速度に基づいて、アセチル補酵素の濃度を再度推定するように構成されている、請求項３に記載のインスリン抵抗性評価支援システム。
5. 前記入力手段は、被験者の筋肉量に関する情報の入力をさらに受け付けるように構成されており、
   前記第２推定手段は、前記第１推定手段によって推定されたグルコーストランスポータの発現量に応じた数値と血中のグルコース濃度とに基づいて、被験者の単位筋肉量当たりのグルコース取込速度を推定するように構成されており、
   前記推定手段は、前記第２推定手段によって推定された被験者の単位筋肉量当たりのグルコース取込速度と、前記入力手段によって入力を受け付けた被験者の筋肉量に関する情報とに基づいて、被験者の単位重量当たりのグルコース取込速度を推定するように構成されている、請求項２乃至４の何れかに記載のインスリン抵抗性評価支援システム。
6. 前記推定手段は、生体内におけるグルコースからピルビン酸への変化に関連する複数の物質の生成反応において、生成反応前の各物質の濃度に基づいて、各物質の一定時間あたりの生成量を取得し、それぞれの物質について、前記生成反応前の物質の濃度にその物質の前記生成量を反映させることにより、前記一定時間後の各物質の濃度を取得し、取得された各物質の濃度に基づいて、グルコース取込速度を推定する処理を、繰り返し実行するように構成されている、請求項１乃至５のいずれかに記載のインスリン抵抗性評価支援システム。
7. 前記推定手段は、グルコース取込速度が定常状態に至ったか否かを判定し、グルコース取込速度が定常状態に至ったと判定されるまで、前記各物質の濃度を取得する処理を繰り返し実行するように構成されている、請求項６に記載のインスリン抵抗性評価支援システム。
8. 前記入力手段は、被験者の空腹時における血中のグルコース濃度、インスリン濃度、及び遊離脂肪酸濃度に関する情報の入力を受け付けるように構成されており、
   前記推定手段は、
   前記入力手段によって入力を受け付けた血中のグルコース濃度、インスリン濃度、及び遊離脂肪酸濃度に関する情報に基づいて、空腹時における各物質の濃度を取得し、取得された各物質の濃度に基づいて、空腹時におけるグルコース取込速度を推定する第１グルコース取込速度推定手段と、
   前記第１グルコース取込速度推定手段によって取得された空腹時における各物質の濃度に基づいて、インスリン濃度が所定値の場合における各物質の濃度を取得し、取得された各物質の濃度に基づいて、インスリン濃度が所定値の場合におけるグルコース取込速度を推定する第２グルコース取込速度推定手段と、を備える、請求項６又は７のいずれかに記載のインスリン抵抗性評価支援システム。
9. 前記推定手段によって推定されたグルコース取込速度を出力する出力手段をさらに備える、請求項１乃至８のいずれかに記載のインスリン抵抗性評価支援システム。
10. 前記推定手段によって推定されたグルコース取込速度に関する情報に基づいて、被験者のインスリン抵抗性を推定するインスリン抵抗性推定手段と、
    前記インスリン抵抗性推定手段による推定結果を出力する出力手段と、
    をさらに備える、請求項１乃至８のいずれかに記載のインスリン抵抗性評価支援システム。
11. 入力装置を備えるコンピュータにより、被験者のインスリン抵抗性の評価を支援するためのインスリン抵抗性評価支援方法であって、
    被験者を測定することにより得られた血中のグルコース濃度、インスリン濃度、及び遊離脂肪酸濃度に関する情報の入力を前記入力装置により受け付けるステップと、
    前記入力装置によって入力を受け付けた血中のグルコース濃度、インスリン濃度、及び遊離脂肪酸濃度に関する情報に基づいて、前記コンピュータにより被験者のグルコース取込速度を推定するステップと、
    を有する、インスリン抵抗性評価支援方法。
12. 前記被験者のグルコース取込速度を推定するステップは、
    血中のインスリン濃度及び遊離脂肪酸濃度に関する情報に基づいて、前記演算部によりグルコーストランスポータの発現量に応じた数値を推定するステップと、
    推定されたグルコーストランスポータの発現量に応じた数値と血中のグルコース濃度とに基づいて、前記コンピュータによりグルコース取込速度を推定するステップと、
    を含む、請求項１１に記載のインスリン抵抗性評価支援方法。
13. 入力装置を備えるコンピュータに、
    被験者を測定することにより得られた血中のグルコース濃度、インスリン濃度、及び遊離脂肪酸濃度に関する情報の入力を前記入力装置により受け付けるステップと、
    前記入力手段によって入力を受け付けた血中のグルコース濃度、インスリン濃度、及び遊離脂肪酸濃度、並びに被験者の血中のグルコース濃度に基づいて、被験者のグルコース取込速度を推定するステップと、
    を実行させる、コンピュータプログラム。
14. 前記被験者のグルコース取込速度を推定するステップにおいて、前記コンピュータに、
    血中のインスリン濃度及び遊離脂肪酸濃度に関する情報に基づいて、グルコーストランスポータの発現量に応じた数値を推定するステップと、
    推定されたグルコーストランスポータの発現量に応じた数値と血中のグルコース濃度とに基づいて、グルコース取込速度を推定するステップと、
    を実行させる、請求項１３に記載のコンピュータプログラム。
    

Images