(57)【要約】【課題】生きている細胞内のカーゴにかかる力を非侵襲的な方法で測定する技術を提供すること。【解決手段】本発明の一態様は、非侵襲的な手段によって取得したカーゴの移動の軌跡を示すタイムコースデータに含まれるデータのうち、カーゴの単位時間の変位が略一定である期間におけるカーゴの位置の時間変化を示す第一のデータを取得する第一取得ステップと、前記第一のデータに基づく確率分布であって、単位時間における前記カーゴの所定の位置の変化量が所定の量である確率の分布を示す第一の確率分布に基づいて、前記カーゴを輸送するモータータンパク質が前記カーゴにくわえる力に比例するゆらぎ値を算出するゆらぎ値算出ステップと、算出した前記ゆらぎ値に基づいて、前記力を算出する力算出ステップと、を有する力測定方法である。【選択図】図２
2019082432.tif 000002
【特許請求の範囲】
【請求項１】
非侵襲的な手段によって取得したカーゴの移動の軌跡を示すタイムコースデータに含まれるデータのうち、カーゴの単位時間の変位が略一定である期間におけるカーゴの位置の時間変化を示す第一のデータを取得する第一取得ステップと、
前記第一のデータに基づく確率分布であって、単位時間における前記カーゴの所定の位置の変化量が所定の量である確率の分布を示す第一の確率分布に基づいて、前記カーゴを輸送するモータータンパク質が前記カーゴにくわえる力に比例するゆらぎ値を算出するゆらぎ値算出ステップと、
算出した前記ゆらぎ値に基づいて、前記力を算出する力算出ステップと、
を有する力測定方法。
【請求項２】
前記力算出ステップは、前記ゆらぎ値に対してクラスタリングを行うことで前記力を算出する、
請求項１に記載の力測定方法。
【請求項３】
非侵襲的な手段によって取得したカーゴの移動の軌跡を示すタイムコースデータのうち、前記第一のデータと、前記カーゴの位置が時間によらず略同一の位置である期間の前記カーゴの位置の時間変化を示す第二のデータとを含むタイムコースデータを取得する第二取得ステップと、
前記力と前記ゆらぎ値との換算係数であるゆらぎ係数を前記第一及び第二のデータに基づいて算出するゆらぎ係数算出ステップと、
を有する請求項１又は２に記載の力測定方法。
【請求項４】
前記ゆらぎ係数算出ステップは、前記第二のデータに基づいて、前記カーゴを囲む媒質の粘性係数を算出する、
請求項３に記載の力測定方法。
【請求項５】
前記ゆらぎ係数算出ステップは、前記第二取得ステップにおいて取得された前記第一のデータに基づいて、前記カーゴが移動する速さを算出する請求項４に記載の力測定方法。
【請求項６】
前記ゆらぎ係数算出ステップは、前記粘性係数と前記速さとに基づいて前記ゆらぎ係数を算出する、
請求項５に記載の力測定方法。
【請求項７】
前記ゆらぎ値算出ステップは、前記位置の変化量をΔＸとし、前記単位時間におけるΔＸの値がΔＸ１である確率を第一の確率Ｐ１（ΔＸ１）とし、前記ゆらぎ値をχとした場合に、χ＝ｌｎ［Ｐ１（ΔＸ１）/Ｐ１（−ΔＸ１）］／ΔＸ１の式によってχを算出する、
請求項１に記載の力測定方法。
【請求項８】
非侵襲的な手段によって取得したカーゴの移動の軌跡を示すタイムコースデータに含まれるデータのうち、カーゴの単位時間の変位が略一定である期間におけるカーゴの位置の時間変化を示す第一のデータを取得する第一取得部と、
前記第一のデータに基づく確率分布であって、単位時間における前記カーゴの所定の位置の変化量が所定の量である確率の分布を示す第一の確率分布に基づいて、前記カーゴを輸送するモータータンパク質が前記カーゴにくわえる力に比例するゆらぎ値を算出するゆらぎ値算出部と、
算出した前記ゆらぎ値に基づいて、前記力を算出する力算出部と、
を備える力測定装置。
【請求項９】
非侵襲的な手段によって取得したカーゴの移動の軌跡を示すタイムコースデータに含まれるデータのうち、カーゴの単位時間の変位が略一定である期間におけるカーゴの位置の時間変化を示す第一のデータを取得する第一取得部と、
前記第一のデータに基づく確率分布であって、単位時間における前記カーゴの所定の位置の変化量が所定の量である確率の分布を示す第一の確率分布に基づいて、前記カーゴを輸送するモータータンパク質が前記カーゴにくわえる力に比例するゆらぎ値を算出するゆらぎ値算出部と、
算出した前記ゆらぎ値に基づいて、前記力を算出する力算出部と、
を備える力測定システム。
【請求項１０】
コンピュータに、非侵襲的な手段によって取得したカーゴの移動の軌跡を示すタイムコースデータに含まれるデータのうち、カーゴの単位時間の変位が略一定である期間におけるカーゴの位置の時間変化を示す第一のデータを取得させ、
前記コンピュータに、前記第一のデータに基づく確率分布であって、単位時間における前記カーゴの所定の位置の変化量が所定の量である確率の分布を示す第一の確率分布に基づいて、前記カーゴを輸送するモータータンパク質が前記カーゴにくわえる力に比例するゆらぎ値を算出させ、
前記コンピュータに、算出した前記ゆらぎ値に基づいて、前記力を算出するための力測定プログラム。
【請求項１１】
コンピュータに、非侵襲的な手段によって取得したカーゴの移動の軌跡を示すタイムコースデータに含まれるデータのうち、カーゴの単位時間の変位が略一定である期間におけるカーゴの位置の時間変化を示す第一のデータを取得させ、
前記コンピュータに、前記第一のデータに基づく確率分布であって、単位時間における前記カーゴの所定の位置の変化量が所定の量である確率の分布を示す第一の確率分布に基づいて、前記カーゴを輸送するモータータンパク質が前記カーゴにくわえる力に比例するゆらぎ値を算出させ、
前記コンピュータに、算出した前記ゆらぎ値に基づいて、前記力を算出するためのプログラムを記憶する
記録媒体。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、力測定方法、力測定装置、力測定システム、力測定プログラム及び記録媒体に関する。
【背景技術】
【０００２】
多くの真核細胞は発達した細胞内輸送システムを備えている。従来、神経細胞が細胞内輸送システムの研究対象として用いられてきた。神経細胞は、細胞体と細胞体から延びる長い突起（軸索）とを有する。神経細胞の形態は細胞骨格を構成するタンパク質群によって形成されており、軸索は主に微小管を構成するチューブリンによって形成されている。微小管は、細胞体と軸索の先端にあるシナプスとの間の物流システム（軸索輸送）の幹線として機能する。軸索輸送で運ばれるカーゴ（Cargo）の例として、脂質二重膜によって覆われた細胞内小器官（オルガネラ）や小胞（ベシクル）が知られている。細胞小器官の例としては、例えば、ミトコンドリアやエンドソームが知られている。小胞内には、軸索やシナプスの形成材料又は情報伝達物質としてのタンパク質、核酸高分子、低分子等が含まれている。カーゴの輸送には、細胞体から軸索先端側へ向かう順行性輸送と、その逆方向の逆行性輸送とがある。順行性輸送では、キネシンと呼ばれるモータータンパク質が、カーゴを背負って微小管の上を辿りながら移動することが知られている。同様に、逆行性輸送ではダイニンと呼ばれるモータータンパク質が、カーゴの輸送を担っている。モータータンパク質は、ＡＴＰの加水分解によって得たエネルギーを用いてカーゴを輸送する。軸索輸送の乱れは、アルツハイマー病やパーキンソン病などの神経性疾患の一因とも考えられており、軸索輸送の解析が神経性疾患の理解や治療に貢献すると期待されている。
ここまで、軸索輸送を具体例としてモーター蛋白質によるカーゴの輸送を説明したが、モータータンパク質によるカーゴの輸送は、軸索輸送に限られない。神経細胞においても、軸索のみならず、樹状突起などの突起内あるいは細胞体におけるカーゴの輸送が行われる。その際、カーゴは必ずしも脂質二重膜に覆われているとは限らず、タンパク質複合体あるいはタンパク質・核酸複合体などの場合がある。また、モータータンパク質は、必ずしも微小管の上を辿って移動するキネシンあるいはダイニンだけが用いられるわけではなく、アクチン繊維上を移動するミオシンなど様々なモーター分子が用いられる場合がある。また、細胞内輸送は、神経細胞に限らず、全ての真核細胞で同様の機構により行われている。例えば皮膚のメラニン細胞(メラノサイト)では、メラニン色素を含む細胞内顆粒がカーゴとして輸送され、膵臓のインスリン分泌細胞では、インスリンを含む細胞内小胞がカーゴとして輸送される。すなわち、ここまで、神経細胞を具体例としてモータータンパク質によるカーゴの輸送を説明したが、モータータンパク質によるカーゴの輸送は、必ずしも神経細胞だけで行われるわけではなく、神経細胞以外の細胞であっても行われる。例えば、神経細胞以外の細胞では、モータータンパク質によるカーゴの輸送は、軸索の他、微小管の上をモータータンパク質が移動することで行われる。
【０００３】
従来、このようなカーゴの輸送を解析するため、生きている細胞内のカーゴにかかる力を測定することが行われていた。生きている細胞内のカーゴにかかる力を測定する方法としては、例えば、光ピンセットを用いた方法があった。光ピンセットを用いた方法は、光ピンセットによって観察対象のカーゴに力をくわえて、カーゴを引っ張ることでカーゴにかかる力を測定する方法であった。なお、カーゴにかかる力とは、例えば、モータータンパク質がカーゴにかける力である。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２００９−７７６３５号公報
【非特許文献】
【０００５】
【非特許文献１】K. Hayashi, C.G. Pack, M.K. Sato, K. Mouri, K. Kaizu, K. Takahashi, and Y. Okada、"Viscosity and drag force involved in organelle transport: Investigation of the fluctuation dissipation theorem"、The European Physical Jounal E, 36: 136 (2013)
【発明の概要】
【発明が解決しようとする課題】
【０００６】
しかしながら、光ピンセットを用いた方法は、カーゴとその周囲の細胞質の屈折率差及びカーゴの大きさとレーザー光の強度が発生力を決めるため、油滴のような特定のカーゴ以外への適用が難しい。その結果、光ピンセットを用いた方法を油滴等以外のカーゴに適用しようとすると過大な強度のレーザー光が必要であり、細胞侵害・熱的ダメージを生じる場合があった。
このように、光ピンセットを用いた方法は、レーザー光を用いる侵襲的な方法であるため観察対象の細胞を傷つける可能性があるという問題があった。
【０００７】
上記事情に鑑み、本発明は、生きている細胞内のカーゴにかかる力を非侵襲的な方法で測定する技術を提供することを目的としている。
【課題を解決するための手段】
【０００８】
本発明の一態様は、非侵襲的な手段によって取得したカーゴの移動の軌跡を示すタイムコースデータに含まれるデータのうち、カーゴの単位時間の変位が略一定である期間におけるカーゴの位置の時間変化を示す第一のデータを取得する第一取得ステップと、前記第一のデータに基づく確率分布であって、単位時間における前記カーゴの所定の位置の変化量が所定の量である確率の分布を示す第一の確率分布に基づいて、前記カーゴを輸送するモータータンパク質が前記カーゴにくわえる力に比例するゆらぎ値を算出するゆらぎ値算出ステップと、算出した前記ゆらぎ値に基づいて、前記力を算出する力算出ステップと、を有する力測定方法である。
【０００９】
本発明の一態様は、上記の力測定方法であって、前記力算出ステップは、前記ゆらぎ値に対してクラスタリングを行うことで前記力を算出する。
【００１０】
本発明の一態様は、上記の力測定方法であって、非侵襲的な手段によって取得したカーゴの移動の軌跡を示すタイムコースデータのうち、前記第一のデータと、前記カーゴの位置が時間によらず略同一の位置である期間の前記カーゴの位置の時間変化を示す第二のデータとを含むタイムコースデータを取得する第二取得ステップと、
前記力と前記ゆらぎ値との換算係数であるゆらぎ係数を前記第一及び第二のデータに基づいて算出するゆらぎ係数算出ステップとを有する。
【００１１】
本発明の一態様は、上記の力測定方法であって、前記ゆらぎ係数算出ステップは、前記第二のデータに基づいて、前記カーゴを囲む媒質の粘性係数を算出する。
【００１２】
本発明の一態様は、上記の力測定方法であって、前記ゆらぎ係数算出ステップは、前記第二取得ステップにおいて取得された前記第一のデータに基づいて、前記カーゴが移動する速さを算出する。
【００１３】
本発明の一態様は、上記の力測定方法であって、前記ゆらぎ係数算出ステップは、前記粘性係数と前記速さとに基づいて前記ゆらぎ係数を算出する。
【００１４】
本発明の一態様は、上記の力測定方法であって、前記ゆらぎ値算出ステップは、前記位置の変化量をΔＸとし、前記単位時間におけるΔＸの値がΔＸ１である確率を第一の確率Ｐ１（ΔＸ１）とし、前記ゆらぎ値をχとした場合に、χ＝ｌｎ［Ｐ１（ΔＸ１）/Ｐ１（−ΔＸ１）］／ΔＸ１の式によってχを算出する。
【００１５】
本発明の一態様は、非侵襲的な手段によって取得したカーゴの移動の軌跡を示すタイムコースデータに含まれるデータのうち、カーゴの単位時間の変位が略一定である期間におけるカーゴの位置の時間変化を示す第一のデータを取得する第一取得部と、前記第一のデータに基づく確率分布であって、単位時間における前記カーゴの所定の位置の変化量が所定の量である確率の分布を示す第一の確率分布に基づいて、前記カーゴを輸送するモータータンパク質が前記カーゴにくわえる力に比例するゆらぎ値を算出するゆらぎ値算出部と、算出した前記ゆらぎ値に基づいて、前記力を算出する力算出部と、を備える力測定装置である。
【００１６】
本発明の一態様は、非侵襲的な手段によって取得したカーゴの移動の軌跡を示すタイムコースデータに含まれるデータのうち、カーゴの単位時間の変位が略一定である期間におけるカーゴの位置の時間変化を示す第一のデータを取得する第一取得部と、前記第一のデータに基づく確率分布であって、単位時間における前記カーゴの所定の位置の変化量が所定の量である確率の分布を示す第一の確率分布に基づいて、前記カーゴを輸送するモータータンパク質が前記カーゴにくわえる力に比例するゆらぎ値を算出するゆらぎ値算出部と、算出した前記ゆらぎ値に基づいて、前記力を算出する力算出部と、を備える力測定システムである。
【００１７】
本発明の一態様は、コンピュータに、非侵襲的な手段によって取得したカーゴの移動の軌跡を示すタイムコースデータに含まれるデータのうち、カーゴの単位時間の変位が略一定である期間におけるカーゴの位置の時間変化を示す第一のデータを取得させ、前記コンピュータに、前記タイムコースデータに基づく確率分布であって、単位時間中における前記カーゴの所定の位置の変化量が所定の量である確率の分布を示す第一の確率分布に基づいて、前記カーゴを輸送するモータータンパク質が前記カーゴにくわえる力に比例するゆらぎ値を算出させ、前記コンピュータに、算出した前記ゆらぎ値に基づいて、前記力を算出させるための力測定プログラムである。
【００１８】
本発明の一態様は、コンピュータに、非侵襲的な手段によって取得したカーゴの移動の軌跡を示すタイムコースデータに含まれるデータのうち、カーゴの単位時間の変位が略一定である期間におけるカーゴの位置の時間変化を示す第一のデータを取得させ、前記コンピュータに、前記タイムコースデータに基づく確率分布であって、単位時間中における前記カーゴの所定の位置の変化量が所定の量である確率の分布を示す第一の確率分布に基づいて、前記カーゴを輸送するモータータンパク質が前記カーゴにくわえる力に比例するゆらぎ値を算出させ、前記コンピュータに、算出した前記ゆらぎ値に基づいて、前記力を算出させるためのプログラムを記憶する記録媒体である。
【発明の効果】
【００１９】
本発明により、生きている細胞内のカーゴにかかる力を非侵襲的な方法で測定することが可能となる。
【図面の簡単な説明】
【００２０】
【図１】ゆらぎの定理が力Ｆｃを反映することを説明する説明図。
【図２】実施形態の力測定装置１を用いて構成された力測定システム１００の具体例を示す図。
【図３】実施形態の観察データの具体例を示す図。
【図４】実施形態の複数の観察データが観察対象の移動を示すことを説明する説明図。
【図５】実施形態の力測定装置１の機能構成の具体例を示す図。
【図６】実施形態のタイムコースデータ、定速移動データ及び定速期間の具体例を示す図。
【図７】実施形態の第一確率分布算出処理における変位量算出処理の具体的な手順を説明する説明図。
【図８】実施形態の第一確率分布の具体例な算出結果を示す図。
【図９】実施形態のゆらぎ値χの具体的な算出結果を示す図。
【図１０】実施形態のクラスタリング部１４によってクラスタリングされたゆらぎ値χを示す図。
【図１１】実施形態のゆらぎ係数算出装置２の機能構成の具体例を示す図。
【図１２】実施形態の静止データ及び定速移動データを含むタイムコースデータの具体例を示す図。
【図１３】実施形態の力測定装置１がモータータンパク質にかかる力を算出する具体的な処理の流れを示すフローチャート。
【図１４】実施形態の単位時間繰り返し制御（ステップＳ１０３）の具体的な処理の流れを示すフローチャート。
【図１５】Ciliobrevinの細胞への添加によってゆらぎ値χの度数分布が変化することが、実施形態の力測定装置１を用いたχの測定によって明らかになったことを示す具体的な結果を示す図。
【図１６】エンドソームが細胞中心と細胞膜との間をキネシン又はダイニンによって輸送されることを説明する説明図。
【図１７】エンドソームがキネシンによって輸送場合とダイニンによって輸送される場合とのゆらぎ値χの度数分布の違いが実施形態の力測定装置１を用いた測定で明らかになったことを示す図。
【発明を実施するための形態】
【００２１】
（原理）
本発明は、カーゴにかかる力に対してゆらぎの定理によって関係づけられた物理量を測定することでカーゴにかかる力を測定する。そこで、本発明の説明にあたって、まずはゆらぎの定理を説明した後、本発明によってカーゴにかかる力が測定される原理を説明する。
【００２２】
ゆらぎの定理はD.J.Evansらによって１９９３年に発表された、物理学の１分野である非平衡統計力学の定理である。ゆらぎの定理は、熱力学的に非平衡な状態にある系における所定の物理量と、その所定の物理量を変化させる物理量とを関係づける定理である。熱力学的に非平衡な状態にある系とは、例えば、カーゴが移動する系であり、その場合、ゆらぎの定理は、その系におけるカーゴの位置と、カーゴにかかる力とを関係づける。その場合、移動するカーゴの位置とその移動するカーゴにかかる力とを関係づけるゆらぎの定理は、次の式で表わされることが知られている。
【００２３】
【数１】
2019082432.tif 000003
【００２４】
ΔＸは、単位時間中のカーゴの変位量を表し、Ｐ（ΔＸ）は、カーゴの変位量がΔＸである確率を表す。なお、ΔＸは、以下の式（２）で示される。
【００２５】
【数２】
2019082432.tif 000004
【００２６】
式（２）において、Ｘ（ｔ）は、時刻ｔにおけるカーゴの位置である。式（２）において、Ｘ（ｔ＋Δｔ）は、時刻ｔ＋Δｔにおけるカーゴの位置である。
ΔＸは、力Ｆｃ（式（１）ではＦと表記）は、カーゴにかかる力の大きさを表す。Ａは、式（１）の左辺が表す物理量と力Ｆｃとの換算係数である。以下、式（１）の左辺が表す物理量をゆらぎ値χといい、Ａをゆらぎ係数という。ゆらぎ係数Ａは、カーゴとモータータンパク質との組み合わせに固有の値である。
【００２７】
カーゴの変位量が測定されれば、ゆらぎ値χが算出される。そのため、ゆらぎ係数Ａが与えられれば、式（１）によりカーゴにかかる力Ｆｃが算出され得る。
なお、カーゴの変位量は、カーゴの位置の変化量を表す。また、カーゴの位置とは、カーゴの部位の位置であればどのような位置であってもよく、例えば、カーゴの重心位置である。
【００２８】
図１は、ゆらぎの定理が力Ｆｃを反映することを説明する説明図である。図１において、図１（Ａ）及び（Ｂ）は、どちらもカーゴに力がかかる場合におけるΔＸの確率分布を表す。ただし、図１（Ａ）は、カーゴにかかるノイズが大きい場合のΔＸの確率分布を表す。
図１（Ａ）及び（Ｂ）における確率分布は、例えば、どちらもガウス分布に従う確率分布である。図１（Ａ）と（Ｂ）とを比較すると、図１（Ａ）においては、カーゴにかかるノイズが小さいため、ΔＸがとり得る値の範囲が図１（Ｂ）におけるΔＸがとり得る値の範囲よりも狭い。
そのため、図１（Ａ）における確率分布のピークの高さは、図１（Ａ）における確率分布のピークの高さよりも低い。しかしながらP(ΔX)とP(-ΔX)の比はどちらも等しく、P(ΔX)そのものでなくP(ΔX)とP(-ΔX)の比が力Ｆｃを反映していることを表す。
そのため、ゆらぎの定理は力Ｆｃを反映する定理であると言える。
【００２９】
カーゴの変位が測定された場合には、式（１）の左辺の値だけでなくゆらぎ係数Ａも算出され得る。このことを以下で説明する。粘性のある媒質中を速さＶで定速移動するカーゴにおいては、モータータンパク質によってカーゴにくわえられる力と、媒質による粘性抵抗の力とが釣り合う。そのため、媒質の粘性係数Γがわかれば、モータータンパク質によってカーゴにくわえられる力Ｆｃが算出され得る。具体的には、力Ｆｃと粘性係数Γとは以下の式（３）の関係を満たす。
【００３０】
【数３】
2019082432.tif 000005
【００３１】
そのため、ゆらぎ値χと、粘性係数Γに基づいて算出された力Ｆｃとに基づいて、ゆらぎ係数Ａが算出される。なお、ゆらぎ係数Ａは、ボルツマン定数ｋｂと、非平衡状態における観察対象のカーゴを含む系の温度Ｔｅｆｆ（すなわち、一般に非平衡温度又は有効温度といわれる温度）との積である。非平衡状態は、カーゴが輸送されている状態である。モータータンパク質によるカーゴの輸送においては、Ｔｅｆｆは一定値である。このことは実験的に検証済みである。
【００３２】
（実施形態）
図２は、実施形態の力測定装置１を用いて構成された力測定システム１００の具体例を示す図である。力測定システム１００は、非侵襲的な手段によって細胞中のカーゴを検出し、そのカーゴの動きを記録することで、そのカーゴにかかる力を測定するシステムである。力測定システム１００は、力測定装置１、ゆらぎ係数算出装置２及び撮像装置９を備える。
【００３３】
なお、本発明が観察対象とする細胞は、カーゴ及びモータータンパク質を有し、細胞内でモータータンパク質がカーゴを輸送する細胞であれば、どのような細胞であってもよい。一般的に知られる真核細胞の多くはこのようなカーゴ及びモータータンパク質を有する。細胞内におけるカーゴの輸送を観察し易いことから、観察対象は神経細胞であることが好ましい。神経細胞以外で観察対象として好適な細胞としては、例えばグリア細胞、血管内皮細胞、皮膚などの繊維芽細胞、メラニン細胞、一般的に用いられる株化された動物細胞、卵細胞等が挙げられる。
【００３４】
本発明が観察対象とする神経細胞の種類は特に限定されず、例えば、神経細胞内の動画撮影が容易である観点からin vitroで培養可能な神経細胞が好ましい。また、神経細胞を補助する細胞あるいはその他の細胞を併存させた細胞であってもよい。このような細胞として、例えば、アストロサイト(星状膠細胞)、オリゴデンドロサイト、シュワン細胞などのグリア細胞、血管上皮細胞などの血管系細胞など、繊維芽細胞などの間葉系細胞などが挙げられる。また、このような神経細胞の例として、アミノ酸系神経伝達物質を分泌する神経細胞、ペプチド系神経伝達物質を分泌する神経細胞、モノアミン類やアセチルコリンを分泌する神経細胞等が挙げられる。神経細胞を培養する際には、生体内環境に近づける目的で、神経細胞を補助する細胞を併存させてもよい。このような補助細胞として、例えば、アストロサイト(星状膠細胞)、オリゴデンドロサイト、シュワン細胞等のグリア細胞等が挙げられる。
【００３５】
また、本発明が観察対象とする細胞内で輸送されるカーゴの種類は特に限定されず、光学機器によってゆらぎを測定可能なカーゴであればよい。ゆらぎを測定する際のカーゴの蛍光標識の有無は限定されない。カーゴが脂質二重膜で包まれた小胞又は細胞小器官である場合、例えば、ＤｉＤ、ＤｉＩ、ＤｉＯ、ＤｉＲ等のカルボシアニン色素によって蛍光標識することができる。あるいは、カーゴに含まれるタンパク質をGFPなどの蛍光タンパク質によって蛍光標識することができる。その他、蛍光色素を用いた染色、ビオチン・アビジン、抗体などを用いた公知の標識方法も適用可能である。観察するカーゴの輸送経路は特に限定されず、軸索輸送でもよいし、その他の細胞内輸送であってもよく、解析の目的に応じて適宜選択することができる。
【００３６】
力測定装置１は、観察対象のカーゴの動きを示す画像のデータ群に基づいて、観察対象のカーゴにかかる力Ｆｃを算出する。以下、観察対象のカーゴの動きを示す画像のデータ群を、観察データという。観察対象のカーゴの動きを示す画像は、例えば、二次元画像であってもよい。また、観察データは、例えば動画であってもよい。
【００３７】
ゆらぎ係数算出装置２は、観察データに基づいて、観察対象のカーゴのゆらぎ係数Ａを算出する。ゆらぎ係数算出装置２は、算出したゆらぎ係数Ａを力測定装置１に出力する。
【００３８】
撮像装置９は、非侵襲的な手段によって、観察対象のカーゴを撮影し、そのカーゴの動きを記録した観察データを生成する。非侵襲的な手段とは、例えば、可視光で観察対象に関する情報を取得する手段であってもよい。観察対象に関する情報とは、例えば、観察対象の重心位置を示す信号である。また、非侵襲的な手段とは、例えば、観察対象を蛍光標識しその蛍光によって観察対象に関する情報を取得する手段であってもよい。これらの非侵襲的な手段は、蛍光顕微鏡や、微分干渉顕微鏡等の顕微鏡を用いて実現されてもよい。また、非侵襲的な手段とは、観察対象への影響が少ない強度の磁場によって観察対象に関する情報を取得する手段であってもよい。なお、侵襲的な手段とは、例えば、光ピンセットのような、観察対象へレーザーを照射することで観察対象に関する情報を取得する手段である。
【００３９】
なお、カーゴの動きを記録した観察データを生成する方法は特に限定されず、例えば、市販の細胞動態観察用デジタルビデオカメラを備えた共焦点レーザー顕微鏡などを用いる公知方法が挙げられる。撮影のフレームレートとしては、カーゴの軸索輸送の速度等の観点から、１０〜１００フレーム／秒程度またはそれ以上が好ましく、１００フレーム／秒が特に好ましい。
また、撮影した画像フレーム中の軸索領域を抽出する工程又は手段では、種々の方法を適用することが可能であるが、例えば、目視によって軸索領域を判断する方法が行われる。
【００４０】
図３は、実施形態の観察データの具体例を示す図である。
図３は、ある時刻ｔにおける二次元の静止画像である。図３の静止画像における点線で囲まれた領域は軸索を表す。図３において、線で囲まれた領域１０１内の発光する点は、観察対象のカーゴである。
【００４１】
図４は、実施形態の複数の観察データが観察対象の移動を示すことを説明する説明図である。
図４において、画像１０２は、時刻ｔ＝０ｓにおける観察対象のカーゴの位置を示す。なお、ここで観察対象のカーゴは、図３における領域１０１内の観察対象のカーゴである。図４において、画像１０３は、０ｓより後であって１．８ｓより前の時刻における観察対象のカーゴの位置を示す。図４において、画像１０４は、時刻ｔ＝１．８ｓにおける観察対象のカーゴの位置を示す。
【００４２】
観察データは、図３に示される二次元の静止画像を、時刻を変えて複数備えたデータ群である。そのため、観察データは、自身が備える異なる時刻における静止画像のデータによって、図４に示すような観察対象のカーゴの移動を示す。 
【００４３】
図５は、実施形態の力測定装置１の機能構成の具体例を示す図である。力測定装置１は、第一軌跡データ取得部１１、ゆらぎ値算出部１２、繰り返し制御部１３、クラスタリング部１４、ゆらぎ係数記憶部１５及び力算出部１６を備える。
【００４４】
第一軌跡データ取得部１１は、観察データに基づいて生成された定速移動データを含むタイムコースデータを取得する。タイムコースデータは、観察データが示す観察対象の軌跡の所定の一軸上の射影を示すデータであって、時刻ｔと時刻ｔにおけるカーゴの位置Ｘとの二つの値を有するデータの集合である。所定の一軸とは、どのような軸であってもよく、例えば、カーゴが輸送される方向に平行な軸であってもよい。また、定速移動データは、タイムコースデータに含まれるデータであって、カーゴの単位時間の変位が略一定である期間（以下「定速期間」という。）におけるカーゴの位置の時間変化を示すデータである。すなわち、定速移動データは、カーゴが概ね一定の速度で移動する期間のカーゴの位置の時間変化を示すデータである。定速移動データもまた、時刻ｔとカーゴの位置Ｘとの二つの値を有するデータの集合である。
【００４５】
図６は、実施形態のタイムコースデータ、定速移動データ及び定速期間の具体例を示す図である。
図６のグラフは、時刻ｔ＝０ｓからｔ＝１０ｓまでの時間におけるカーゴの位置Ｘを示すタイムコースデータである。図６において、横軸は時刻を表し、縦軸はカーゴの位置Ｘを表す。図６において縦軸のＸは、位置Ｘが時刻ｔに依存することを明記するためにＸ（ｔ）と表されている。図６において、時刻０から時刻ｔｓ１までの期間のデータは、観察対象のカーゴが概ね一定の速度で移動することを示すデータであり、定速移動データの具体例である。また、時刻０から時刻ｔｓ１までの期間は、定速期間の具体例である。図６には、定速移動データを示すグラフの一部分の拡大図も示されている。拡大図の直線は、定速移動データに対する回帰直線である。回帰直線は、Ｘの値が大きくなる方向に一定の速さでカーゴが移動することを示す。拡大図は、５０ｍｓ以下の短時間ではカーゴの移動の速さは一定ではなく、５０ｍｓで２００ｎｍ以下の位置のゆらぎを有する速さで移動することを示す。
【００４６】
図５の説明に戻る。ゆらぎ値算出部１２は、定速移動データに基づいてゆらぎ値χを算出する。具体的にはゆらぎ値算出部１２は、第一確率分布算出処理を実行した後にゆらぎ値算出処理を実行することで、ゆらぎ値を算出する。以下、第一確率分布算出処理とゆらぎ値算出処理の具体的な処理内容を説明する。まず、第一確率分布算出処理を説明する。
【００４７】
第一確率分布算出処理は、ひとつの定速移動データに基づいて定速期間中の単位時間Δｔにおけるカーゴの位置の変位量ΔＸの確率分布を算出する処理である。具体的には、第一確率分布算出処理は、まず、定速移動データにおける時間原点をｔ０とし、時間原点ｔ０におけるカーゴの位置をＸ０とし、時刻ｔ０＋ｎ×Δｔ（ｎは整数）におけるカーゴの位置をＸｎとした場合の、カーゴの位置の変位量ΔＸｎ＝Ｘｎ−Ｘ（ｎ−１）をさまざまなｎに対して取得する。以下、カーゴの位置の変位量ΔＸｎをさまざまなｎに対して取得する処理を変位量算出処理という。次に第一確率分布算出処理は、カーゴの位置の変位量ΔＸの値がΔＸｎである確率分布（以下「第一確率分布」という。）を算出する。
【００４８】
図７は、実施形態の第一確率分布算出処理における変位量算出処理の具体的な手順を説明する説明図である。
図７のグラフは、図６におけるタイムコースデータの一部分の拡大図におけるグラフである。図７において、ゆらぎ値算出部１２は、まず、変位量算出処理を実行することで、所定の単位時間Δｔごとに時刻ｔ０から時刻ｔｍまでのカーゴの位置（すなわち、Ｘ０、Ｘ１、Ｘ２、・・・Ｘｍ）を取得する。次に、ゆらぎ値算出部１２は、変位量算出処理を実行することで、カーゴの位置の変位量を算出する。すなわち、ΔＸ１＝Ｘ１−Ｘ０、ΔＸ２＝Ｘ２−Ｘ１、ΔＸ３＝Ｘ３−Ｘ２、・・・ΔＸｍ＝Ｘｍ−Ｘｍ−１を算出する。このようにして、ゆらぎ値算出部１２は、所定の単位時間Δｔごとに時刻ｔ０から時刻ｔｍまでのカーゴの位置を取得する。
【００４９】
図５の説明に戻る。ゆらぎ値算出部１２が実行するゆらぎ値算出処理を説明する。ゆらぎ値算出部１２は、ゆらぎ値算出処理によって、第一確率分布と式（１）とに基づいてゆらぎ値χを算出する。ゆらぎ値χは上述したように、式（１）の左辺が表す物理量である。ゆらぎ値算出処理において、式（１）の左辺におけるＰ（ΔＸ）に代入される値は、第一確率分布が表す確率の値であって、カーゴの位置の変位量がΔＸである確率の値である。
【００５０】
繰り返し制御部１３は、複数のゆらぎ値χを算出するようにゆらぎ値算出部１２を制御する。具体的には、繰り返し制御部１３は、まず、ひとつの定速移動データ対して、最短単位時間から最長単位時間までの所定の複数の単位時間Δｔにおけるゆらぎ値χを算出するようにゆらぎ値算出部１２を制御する。以下、この制御を単位時間変更繰り返し制御という。最長単位時間は、ユーザによって予め定められた所定の時間であって、単位時間Δｔがとり得る時間のうち、最短の時間である。最長単位時間は、ユーザによって予め定められた所定の時間であって、単位時間Δｔがとり得る時間のうち、最長の時間である。最短単位時間は例えば、１０ｍｓであってもよい。最長単位時間は例えば、１００ｍｓであってもよい。また、単位時間は例えば、５０ｍｓであってもよい。さらに、繰り返し制御部１３は、単位時間変更繰り返し制御を、第一軌跡データ取得部１１が取得した定速移動データを含むタイムコースデータのそれぞれに対して行う。
なお、単位時間変更繰り返し制御は、例えば、最短単位時間から最長単位時間まで、時間の短い方の単位時間Δｔから順番に各Δｔにおける、ゆらぎ値χを算出するようにゆらぎ値算出部１２を制御する制御であってもよい。
【００５１】
図８及び図９によって、ゆらぎ値算出部１２及び繰り返し制御部１３によって算出された第一確率分布及びゆらぎ値χの具体例を示す。
【００５２】
図８は、実施形態の第一確率分布の具体的な算出結果を示す図である。図８は、単位時間Δｔを１０ｍｓ、５０ｍｓ及び１００ｍｓとした場合の第一確率分布を表す。第一確率分布は、単位時間Δｔが長くなるほど、分布のピークにおけるΔＸの値が大きくなる方向にピークが移動する。この結果は、カーゴが一定の速さで動いているために生じた結果であり、時間が長くなればカーゴが移動する距離が長くなるために生じた結果である。また、第一確率分布は、単位時間Δｔが長くなるほど分散が大きくなる。この結果は、時間が長くなるほどカーゴが移動する距離が長くなり、単位時間におけるカーゴの変位量のばらつきが大きくなるために生じた結果である。
【００５３】
図９は、実施形態のゆらぎ値χの具体的な算出結果を示す図である。より具体的には、図９は、ひとつのタイムコースデータに対する複数の単位時間Δｔごとのχの値を示すグラフである。図９は、単位時間Δｔが大きくなるにつれてゆらぎ値χが一定値に収束することを示す。Δｔが大きくなるにつれてゆらぎ値χが一定の値に収束する結果は、系が緩和することが反映された結果である。系が緩和するとは、観察対象のカーゴを含む系が平衡状態に緩和することを意味する。緩和は、系に含まれる物質のミクロな運動や、その系内のモータータンパク質のＡＴＰ加水分解等によって引き起こされる。
【００５４】
図５の説明に戻る。クラスタリング部１４は、ゆらぎ値算出部１２が算出した複数のゆらぎ値χを対象としてクラスタリングを実行し、測定に関するばらつきを減らしたχの値を出力する。なお、測定に関するばらつきとは、例えば、測定時の観察対象を含む系の温度や湿度等の変化や測定装置の振動等によって引き起こされる取得データのばらつきである。具体的には、クラスタリング部１４はまず、繰り返し制御部１３の制御によってゆらぎ値算出部１２が算出した複数のゆらぎ値χに対してクラスタリングを実行しゆらぎ値χを分類する。クラスタリングは、クラスタ数を自動推定するための基準値を有するクラスタリング方法であればどのような方法で行われてもよい。例えば、基準値であるＡＩＣ（Akaike's Information Criterion)の算出を伴ったk-means methodであってもよい。クラスタ数を自動推定するための基準値はＡＩＣに限らず、例えばＢＩＣ（Bayesian Information Criterion）や、ＭＤＬ（Minimum Description Length)等の基準値であってもよい。クラスタリングする方法はk-means methodに限らず、例えば、k-means++ methodであってもよい。
次にクラスタリング部１４は、クラスタリングした結果の代表値を、測定に関するばらつきを減らしたχとして出力する。なお、代表値はクラスタリングの方法ごとにクラスタを代表する値であればどのような値であってもよく、例えば、平均値であってもよいし、中央値であってもよい。
【００５５】
図１０は、実施形態のクラスタリング部１４によってクラスタリングされたゆらぎ値χを示す図である。図１０は、代表値が離散化された３つのクラスタ（すなわち、第１〜３のクラスタ）を表す。第１のクラスタに属するゆらぎ値χの大きさは、Δｔ＝１００ｍｓにおいて、０以上０．２以下の値である。第２のクラスタに属するゆらぎ値χの大きさは、Δｔ＝１００ｍｓにおいて、０．１以上０．２以下の値である。第３のクラスタに属するゆらぎ値χの大きさは、Δｔ＝１００ｍｓにおいて、０．３以上０．６以下の値である。なお、図１０における代表値は平均値である。第１〜３のクラスタの代表値の値は、小さな順に、第１のクラスタの代表値、第２のクラスタの代表値、第３のクラスタの代表値である。また、第２のクラスタの代表値及び第３のクラスタの代表値は、それぞれ第１のクラスタの代表値の整数倍である。
【００５６】
ゆらぎ係数記憶部１５は、ゆらぎ係数算出装置２が算出したゆらぎ係数Ａを記憶する。ゆらぎ係数記憶部１５は、磁気ハードディスク装置や半導体記憶装置などの記憶装置を用いて構成される。
【００５７】
力算出部１６は、ゆらぎ値χとゆらぎ係数Ａとに基づいて、観察対象のカーゴにかかる力Ｆｃを算出する。具体的には、式（１）と式（１）の右辺がχであることを用いて、Ｆｃを算出する。
【００５８】
図１１は、実施形態のゆらぎ係数算出装置２の機能構成の具体例を示す図である。ゆらぎ係数算出装置２は、第二軌跡データ取得部２１、粘性係数算出部２２及びゆらぎ係数算出部２３を備える。
【００５９】
第二軌跡データ取得部２１は、観察データに基づいて生成された定速移動データ及び静止データを含むタイムコースデータを取得する。静止データは、タイムコースに含まれるデータであって、カーゴの位置が時間によらず略同一の位置である期間（以下「静止期間」という。）のカーゴの位置の時間変化を示すデータである。すなわち、静止データは、カーゴが時刻によらず同じ位置にとどまっている期間のカーゴの位置の時間変化を示すデータである。静止データもまたタイムコースデータと同様に、時刻ｔとカーゴの位置Ｘとの二つの値を有するデータの集合である。なお、定速移動データを含むタイムコースデータのうち、定速移動データだけではなく、さらに静止データをも含むタイムコースデータの割合は、１割程度である。
【００６０】
図１２は、実施形態の静止データ及び定速移動データを含むタイムコースデータの具体例を示す図である。
図１２のグラフは、図６のグラフと同様のグラフであり、時刻ｔ＝０ｓからｔ＝１０ｓまでの期間におけるカーゴの位置Ｘを示すタイムコースデータを示すグラフである。図１２において、横軸は時刻を表し、縦軸はカーゴの位置Ｘを表す。図１２においても、図６と同様に縦軸のＸは、位置Ｘが時刻ｔに依存することを明記するためにＸ（ｔ）と表されている。図１２において、時刻ｔｓ１から時刻ｔｓ2までの期間のデータは、観察対象のカーゴが時刻によらず同じ位置にとどまっていることを示すデータである。すなわち、図１２における時刻ｔｓ１から時刻ｔｓ2までの期間のデータは、静止データの具体例である。また、時刻ｔｓ１から時刻ｔｓ2までの期間は、静止期間の具体例である。なお、時刻０から時刻ｔｓ1までの期間のデータは、図５と同様の定速移動データの具体例である。
【００６１】
図１１の説明に戻る。粘性係数算出部２２は、静止データに基づいて、粘性係数を算出する。具体的には、粘性係数算出部２２は、第二確率分布算出処理、分布係数算出処理、パワースペクトル算出処理、パワー係数算出処理及び粘性係数算出処理を順次実行することで、粘性係数を算出する。
【００６２】
第二確率分布算出処理は、ひとつの静止データに基づいて静止期間中の時刻ｔにおけるカーゴの位置Ｘの確率分布（以下「第二確率分布」という。）を算出する処理である。第二確率分布は、ボルツマン分布にしたがう分布である。
【００６３】
分布係数算出処理は、第二確率分布算出処理によって算出された第二確率分布に対して、フィッティングパラメータとしてＫを含む関数であるボルツマン関数ＰＢ（Ｘ）によるフィッティングを行い、フィッティングパラメータＫの値を算出する。ボルツマン関数ＰＢ（Ｘ）の具体的な関数形は、式（４）で表される関数形である。以下、フィッティングパラメータＫを分布係数という。なお、フィッティングは、どのようなフィッティング方法で行われてもよく、例えば、最小二乗法によって行われてもよい。
【００６４】
【数４】
2019082432.tif 000006
【００６５】
ＰＢ（Ｘ）は、カーゴの位置がＸである確率を表す。Ｚ及びＴはそれぞれ、分配関数及び環境の温度を表す。環境の温度とは、観察対象のカーゴを取り囲む媒質の温度である。
【００６６】
パワースペクトル算出処理は、第二確率分布を周波数νによってフーリエ変換し、パワースペクトルを算出する処理である。
【００６７】
パワー係数算出処理は、フィッティングパラメータとしてＡｐ及びＢｐを含むローレンチアン型の関数ｇ（ν）によってパワースペクトルをフィッティングし、フィッティングパラメータＢｐを算出処理である。関数ｇ（ν）の具体的な関数形は式（５）で表される関数形である。以下、フィッティングパラメータＢｐをパワー係数という。フィッティングは、どのようなフィッティング方法で行われてもよく、例えば、最小二乗法によって行われてもよい。
【００６８】
【数５】
2019082432.tif 000007
【００６９】
粘性係数算出処理は、分布係数Ｋとパワー係数Ｂｐとを用いて、以下の式（６）によって表される粘性係数Γを算出する処理である。
【００７０】
【数６】
2019082432.tif 000008
【００７１】
ゆらぎ係数算出部２３は、粘性係数Γとゆらぎ値χとに基づいて、ゆらぎ係数Ａを算出する。粘性係数算出部は、速さ算出処理及びゆらぎ係数算出処理を順次実行することで、ゆらぎ係数Ａを算出する。
【００７２】
速さ算出処理は、静止データを含むタイムコースデータの定速移動データを用いて、定速移動期間におけるカーゴの時間平均の速さＶを算出する処理である。具体的には、速さ算出処理は、横軸を時刻とし縦軸をカーゴの位置とするグラフにおける定速移動データの回帰直線の傾きを算出する処理である。
【００７３】
ゆらぎ係数算出処理は、速さＶ、粘性係数Γ及びゆらぎ値χによって、ゆらぎ係数Ａを算出する処理である。ゆらぎ係数算出処理において、ゆらぎ係数算出部２３はまず、式（３）によって算出される力Ｆｃとゆらぎ値算出部１２によって算出されたゆらぎ値χとの二つの値を有するデータ（以下「力ゆらぎ値対応データ」という。）を、静止データと定速移動データとを含む複数のタイムコースデータについて取得する。次に、ゆらぎ係数算出処理においてゆらぎ係数算出部２３は、取得した複数の力ゆらぎ値対応データを、フィッティングパラメータＡ１を含む以下の式（７）に示す一次関数によってフィッティングする。
【００７４】
【数７】
2019082432.tif 000009
【００７５】
Ａ１はフィッティングパラメータであり、フィッティングによって得られたＡ１の値はゆらぎ係数Ａである。フィッティングによって得られたＡ１の値がゆらぎ係数Ａであることは式（１）及び式（６）によって明らかである。
【００７６】
図１３は、実施形態の力測定装置１がモータータンパク質にかかる力を算出する具体的な処理の流れを示すフローチャートである。
撮像装置９は、観察対象のカーゴの観察データを生成する（ステップＳ１０１）。第一軌跡データ取得部１１は、撮像装置９が取得した観察データに基づいて生成された定速移動データを含むタイムコースデータを取得する（ステップＳ１０２）。繰り返し制御部１３は、定速移動データを含むタイムコースデータのうちのひとつに対して（以下「算出対象データ」という。）、所定の複数の単位時間Δｔにおける単位時間繰り返し制御を実行することで、複数の単位時間Δｔにおけるゆらぎ値χを算出する（ステップＳ１０３）。繰り返し制御部１３は、第一軌跡データ取得部１１が取得した定速移動データを含むタイムコースデータのうち、第一確率分布算出処理が行われていないタイムコースデータがあるか否かを判定する（ステップＳ１０４）。第一確率分布算出処理が行われていないタイムコースデータがない場合（ステップＳ１０４：Ｎｏ）、クラスタリング部１４は、ゆらぎ値χに対してクラスタリングを実行し、測定に関するばらつきを減らしたゆらぎ値χを算出する（ステップＳ１０５）。
【００７７】
クラスタリング部１４が測定に関するばらつきを減らしたゆらぎ値χを算出した後、第二軌跡データ取得部２１は、静止データ及び定速移動データを含むタイムコースデータを取得する（ステップＳ１０６）。第二軌跡データ取得部２１が静止データ及び定速移動データを含むタイムコースデータを取得した後、粘性係数算出部２２は、第二確率分布算出処理を実行することで、第二軌跡データ取得部２１が取得した静止データに基づいて第二確率分布を算出する（ステップＳ１０７）。粘性係数算出部２２は、第二確率分布を算出した後、第二確率分布に対して分布係数算出処理を実行し、式（４）で表されるボルツマン関数ＰＢ（Ｘ）における分布係数Ｋの値を算出する（ステップＳ１０８）。分布係数Ｋを算出した、粘性係数算出部２２は、第二確率分布に対してパワースペクトル算出処理を実行し、パワースペクトルを算出する（ステップＳ１０９）。パワースペクトルを算出した粘性係数算出部２２は、粘性係数算出処理を実行し粘性係数Γを算出する（ステップＳ１１０）。
【００７８】
粘性係数算出部２２が粘性係数Γを算出した後、ゆらぎ係数算出部２３は、ゆらぎ係数算出処理を実行し、ゆらぎ係数Ａを算出する（ステップＳ１１１）。ゆらぎ係数Ａが算出された後、力算出部１６は、ステップＳ１０５において算出されたχとステップＳ１１１において算出されたゆらぎ係数Ａとに基づいて、カーゴにかかる力Ｆｃを算出する（ステップＳ１１２）。
【００７９】
一方、ステップＳ１０４において、第一確率分布算出処理が行われていないタイムコースデータがある場合（ステップＳ１０４：Ｙｅｓ）、繰り返し制御部１３は、そのタイムコースデータを新たな算出対象データとする（ステップＳ１１３）。繰り返し制御部１３は、新たな算出対象データに対して、単位時間繰り返し制御を実行する（ステップＳ１０３）。
【００８０】
なお、ステップＳ１０６〜１１０は必ずしも、ステップ１０７の後に実行される必要はなく、ステップＳ１０３〜１０５と並列して実行されてもよいし、ステップＳ１０３〜１０５よりも先に実行されてもよい。
【００８１】
図１４は、実施形態の単位時間繰り返し制御（ステップＳ１０３）の具体的な処理の流れを示すフローチャートである。ゆらぎ値算出部１２は、第一軌跡データ取得部１１が取得したタイムコースデータを取得し、第一確率分布算出処理を実行する（ステップＳ２０１）。ゆらぎ値算出部１２は、第一確率分布算出処理を実行することで、第一確率分布を算出する。第一確率分布を算出したゆらぎ値算出部１２は、ゆらぎ値算出処理を実行することによって、ゆらぎ値χを算出する（ステップＳ２０２）。
【００８２】
このように構成された実施形態の力測定装置１によれば、カーゴの動きを示す二次元画像のデータ群に基づいて観察対象のカーゴにかかる力Ｆｃに比例する値χを算出するゆらぎ値算出部１２を備えるため、ユーザは生きている細胞内のカーゴにかかる力を非侵襲的な方法で測定することができる。
【００８３】
（神経細胞の培養方法（例））
３週齢ＩＣＲマウス（オス）から上頚神経節を採取し、０．５％トリプシンと０．５％コラゲナーゼなどの酵素で処理した。個々に分かれた神経細胞を１０％非働化牛血清入りＤＭＥＭ／Ｆ１２培地でリンスした後、マトリゲルコートされたガラスボトムディッシュに撒いた。この神経細胞を３７℃、ＣＯ２ ５％のインキュベータ内で１０％非働化牛血清と２００ｎｇ／ｍｌ ２．５Ｓ神経成長因子入りＤＭＥＭ／Ｆ１２培地にて２〜４日培養し、力測定装置１を用いた実験に用いた。
【００８４】
（細胞小器官（例：エンドソーム）の染色方法）
培養中の神経細胞の培地に最終濃度１００ｎＭとなるように蛍光色素ＤｉＩを添加し、１０分間染色した後、蛍光観察に用いた。
【００８５】
（蛍光観察の方法（例））
染色後のディッシュをヒーティングプレート（CU-201, Live Cell Instrument）上に設置し、３７℃で保温しながら、メラニン細胞内の染色されたメラニン色素小胞を暗視野顕微鏡（IX71, Olympus）で観察した。100倍の対物レンズ（UPlanFL 100x/1.3, Plympus）で98 frames per secondでEMCCDカメラ（LucaS (Andor)）を用いて動画を撮影した。動画からImageJ(Rasband, 1997)を用いて、メラニン色素小胞の重心位置のタイムコースデータを計算した。
【００８６】
図１５は、モータータンパク質の機能を阻害する試薬であるCiliobrevinのメラニン細胞への添加によってゆらぎ値χの度数分布が変化することが、実施形態の力測定装置１を用いたχの測定によって明らかになったことを示す具体的な結果を示す図である。
【００８７】
図１５（Ａ）はCiliobrevinが細胞に添加されていない場合に力測定装置１によって測定されたゆらぎ値χの度数分布を示す図である。図１５（Ａ）において、χは、データの分布が４つのクラスタ（クラスタ１１４〜１１７）に分類される。
図１５（Ｂ）は細胞培養培地中のCiliobrevin濃度が２０μＭである場合に力測定装置１によって測定されたゆらぎ値χの度数分布を示す図である。図１５（Ｂ）において、χは、データの分布が２つのクラスタ（クラスタ１１８及び１１９）に分類される。
図１５（Ｃ）は細胞培養培地中のCiliobrevin濃度が４０μＭである場合に力測定装置１によって測定されたゆらぎ値χの度数分布を示す図である。図１５（Ｃ）において、χは、データの分布が一つのクラスタ１２０に分類される。
【００８８】
図１５（Ａ）〜（Ｃ）の比較によって、Ciliobrevinの添加量が増加するとクラスタの数が減少することが示される。
【００８９】
図１６は、エンドソームが細胞中心と細胞膜との間をキネシン又はダイニンによって輸送されることを説明する説明図である。
図１６（Ａ）は、エンドソーム２０１が自身に結合するキネシン２０２によって、細胞中心から細胞膜まで微小管にそって輸送されることを示す。
図１６（Ｂ）は、エンドソーム２０１が自身に結合するダイニン２０３によって、細胞膜から細胞中心まで微小管にそって輸送されることを示す。
【００９０】
図１７は、エンドソームがキネシンによって輸送場合とダイニンによって輸送される場合とのゆらぎ値χの度数分布の違いが実施形態の力測定装置１を用いた測定で明らかになったことを示す図である。
【００９１】
図１７（Ａ）は、図１６（Ａ）に示すキネシン２０２によってエンドソーム２０１が輸送された場合のχの分布を示す図である。図１７（Ａ）において、ゆらぎ値χは、データの分布が４つのクラスタ（クラスタ１２１〜１２４）に分類される。クラスタ化はクラスタリング解析の結果である。
図１７（Ｂ）は、図１６（Ｂ）に示すダイニン２０３によってエンドソーム２０１が輸送された場合のχの分布を示す図である。図１７（Ｂ）において、ゆらぎ値χは、データの分布が３つのクラスタ（クラスタ１２５〜１２７）に分類される。
図１７（Ａ）と（Ｂ）との比較によって、エンドソーム２０１の輸送は、細胞中心から細胞膜に向かう場合と、細胞膜から細胞中心に向かう場合とでゆらぎ値χの度数分布が異なることが示される。
【００９２】
（変形例）
力測定装置１及びゆらぎ係数算出装置２の機能部は必ずしも異なる筐体によって実装されなくてもよく、ひとつの筐体に実装されてもよい。
【００９３】
なお、第二軌跡データ取得部２１は、必ずしも、観察データに基づいて、定速移動データ及び静止データを含むタイムコースデータを取得する必要はなく、第一軌跡データ取得部１１が取得したタイムコースデータを取得し、取得したタイムコースデータのうち定速移動データ及び静止データを含むタイムコースデータを取得してもよい。
【００９４】
なお、第一軌跡データ取得部１１が取得する定速移動データを含むタイムコースデータは、観察データに基づいてユーザによって生成されてもよいし、装置によって生成されてもよい。また、第二軌跡データ取得部２１が取得する定速移動データ及び静止データを含むタイムコースデータは、観察データに基づいてユーザによって生成されてもよいし、装置によって生成されてもよい。ゆらぎ値算出部１２がゆらぎ値の算出に使用する定速移動データは、ユーザによって手動で選定されてもよいし、装置によって選定されてもよい。ゆらぎ係数算出装置２によって使用される定速移動データ及び静止データもまた、ユーザによって手動で選定されてもよいし、装置によって選定されてもよい。また、観察対象のカーゴの選定もユーザによって手動で選定されてもよいし、装置によって選定されてもよい。
【００９５】
なお、画像のデータ群と、観察データとは必ずしも二次元のデータでなくてもよく、三次元又は一次元のデータであってもよい。
【００９６】
上述した実施形態における力測定装置１及びゆらぎ係数算出装置２をコンピュータで実現するようにしてもよい。その場合、この機能を実現するためのプログラムをコンピュータ読み取り可能な記録媒体に記録して、この記録媒体に記録されたプログラムをコンピュータシステムに読み込ませ、実行することによって実現してもよい。なお、ここでいう「コンピュータシステム」とは、ＯＳや周辺機器等のハードウェアを含むものとする。また、「コンピュータ読み取り可能な記録媒体」とは、フレキシブルディスク、光磁気ディスク、ＲＯＭ、ＣＤ−ＲＯＭ等の可搬媒体、コンピュータシステムに内蔵されるハードディスク等の記憶装置のことをいう。さらに「コンピュータ読み取り可能な記録媒体」とは、インターネット等のネットワークや電話回線等の通信回線を介してプログラムを送信する場合の通信線のように、短時間の間、動的にプログラムを保持するもの、その場合のサーバやクライアントとなるコンピュータシステム内部の揮発性メモリのように、一定時間プログラムを保持しているものも含んでもよい。また上記プログラムは、前述した機能の一部を実現するためのものであってもよく、さらに前述した機能をコンピュータシステムにすでに記録されているプログラムとの組み合わせで実現できるものであってもよく、ＦＰＧＡ（Field Programmable Gate Array）等のプログラマブルロジックデバイスを用いて実現されるものであってもよい。
【００９７】
上述した各実施形態では、本発明を力測定装置として表現したが、本発明は、力測定方法として表現することもできる。この場合、実施形態における力測定装置に対応した本発明による力測定方法は、非侵襲的な手段によって取得したカーゴの移動の軌跡を示すタイムコースデータに含まれるデータのうち、カーゴの単位時間の変位が略一定である期間におけるカーゴの位置の時間変化を示す第一のデータを取得する第一取得ステップと、前記第一のデータに基づく確率分布であって、単位時間における前記カーゴの所定の位置の変化量が所定の量である確率の分布を示す第一の確率分布に基づいて、前記カーゴを輸送するモータータンパク質が前記カーゴにくわえる力に比例するゆらぎ値を算出するゆらぎ値算出ステップと、算出した前記ゆらぎ値に基づいて、前記力を算出する力算出ステップと、を有する力測定方法として表現することができる。
【００９８】
上記の実施形態における力測定装置に対応した力測定方法において、第一取得ステップは、実施形態における第一軌跡データ取得部１１の動作または処理に相当し、ゆらぎ値算出ステップは、実施形態におけるゆらぎ値算出部１２及び繰り返し制御部１３の動作または処理に相当する。また、力算出ステップは、実施形態における力算出部１６の動作または処理に相当する。
【００９９】
なお、実施形態における第二軌跡データ取得部２１の動作または処理は、第二取得ステップの一例であり、実施形態におけるゆらぎ係数算出部２３の動作または処理はゆらぎ係数算出ステップの一例である。
なお、第一軌跡データ取得部１１は、第一取得部の一例である。なお、第二軌跡データ取得部２１は、第二取得部の一例である。なお、定速移動データは第一のデータの一例である。なお、静止データは第二のデータの一例である。
【０１００】
以上、この発明の実施形態について図面を参照して詳述してきたが、具体的な構成はこの実施形態に限られるものではなく、この発明の要旨を逸脱しない範囲の設計等も含まれる。
【符号の説明】
【０１０１】
１…力測定装置、  ２…ゆらぎ係数算出装置、  １１…第一軌跡データ取得部、  １２…ゆらぎ値算出部、  １３…繰り返し制御部、  １４…クラスタリング部、  １５…ゆらぎ係数記憶部、  １６…力算出部、  ２１…第二軌跡データ取得部、  ２２…粘性係数算出部、  ２３…ゆらぎ係数算出部、  １００…力測定システム
【図１】
2019082432.tif 000010
【図２】
2019082432.tif 000011
【図３】
2019082432.tif 000012
【図４】
2019082432.tif 000013
【図５】
2019082432.tif 000014
【図６】
2019082432.tif 000015
【図７】
2019082432.tif 000016
【図８】
2019082432.tif 000017
【図９】
2019082432.tif 000018
【図１０】
2019082432.tif 000019
【図１１】
2019082432.tif 000020
【図１２】
2019082432.tif 000021
【図１３】
2019082432.tif 000022
【図１４】
2019082432.tif 000023
【図１５】
2019082432.tif 000024
【図１６】
2019082432.tif 000025
【図１７】
2019082432.tif 000026