特許 5671522 イオン移動度の予測で使用するための分子の断面積を推定する方法およびシステム

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】特許公報(B2)

(11)【特許番号】5671522

(24)【登録日】2014年12月26日

(45)【発行日】2015年2月18日

(54)【発明の名称】イオン移動度の予測で使用するための分子の断面積を推定する方法およびシステム

(51)【国際特許分類】

   G06F 19/00 20110101AFI20150129BHJP

   G01N 27/62 20060101ALI20150129BHJP

【ＦＩ】

   G06F19/00 110

   G01N27/62 D

【請求項の数】16

【全頁数】9

(21)【出願番号】特願2012-505241(P2012-505241)

(86)(22)【出願日】2010年4月15日

(65)【公表番号】特表2012-524316(P2012-524316A)

(43)【公表日】2012年10月11日

(86)【国際出願番号】GB2010050630

(87)【国際公開番号】WO2010119293

(87)【国際公開日】20101021

【審査請求日】2013年3月28日

(31)【優先権主張番号】0906467.6

(32)【優先日】2009年4月15日

(33)【優先権主張国】GB

(31)【優先権主張番号】0906466.8

(32)【優先日】2009年4月15日

(33)【優先権主張国】GB

(31)【優先権主張番号】61/185,337

(32)【優先日】2009年6月9日

(33)【優先権主張国】US

(73)【特許権者】

【識別番号】504142097

【氏名又は名称】マイクロマス ユーケー リミテッド

(74)【代理人】

【識別番号】110001173

【氏名又は名称】特許業務法人川口國際特許事務所

(72)【発明者】

【氏名】リチヤードソン，キース・ジヨージ

(72)【発明者】

【氏名】カンプサノ，イアン・デイビツド・グラント

【審査官】 宮地 匡人

(56)【参考文献】

【文献】 Lan Jin, et al.，Conformation of glycosaminoglycans by ion mobility mass spectrometry and molecular modelling，Physical chemistry chemical ohysics，２００７年 ７月，Vol.7 No.19，p.3464-3471

【文献】 Philip D. Mosier, et al.，Prediction of Peptide Ion Collision Cross Sections from Topological Molecular Structure and Amino Acid Parameters，Analytical Chemistry，２００２年 ３月１５日，Vol.74 No.6，p.1360-1370

(58)【調査した分野】（Int.Cl.，ＤＢ名）

Ｇ０６Ｆ １９／００

Ｇ０１Ｎ ２７／６２

ＰｕｂＭｅｄ

(57)【特許請求の範囲】

【請求項1】

分子の断面積を推定する推定方法であって、
Ａ）当該分子の幾何構造を予測するステップと、
Ｂ）分子内部の各原子の断面積に対して値を指定するステップと、
Ｃ）分子内部の各原子に対する断面積の値の総和Ａｍａｘを計算するステップと、
Ｄ）当該の前記分子の予測された幾何構造の配向をランダムに選択するステップと、
Ｅ）原子の断面積に比例する確率で複数Ｋ個の原子をランダムに選択するステップとを含み、ここで、Ｋは１以上であり、前記推定方法はさらに、
Ｆ）選択された原子のそれぞれの断面積内部の位置をランダムに選択するステップと、
Ｇ）各前記ランダムに選択された位置に関して、断面積が前記ランダムに選択された位置を含む原子の数ｍを識別するステップと、
Ｈ）各前記ランダムに選択された位置に関してＡｍａｘ／（ｍ×Ｋ）だけパラメータを増分することによって、前記ランダムに選択された配向での当該の前記分子の断面積に対する値を計算するステップと、
Ｉ）すべての前記ランダムに選択された配向全体の平均断面積を計算するステップと、
Ｊ）少なくともＮ回の反復の間、および／または平均断面積のＭ回の独立した計算がＸ回の反復の間、所定の範囲Ｒ以内に一致するまで、ステップＤ）からステップＩ）を反復するステップとを含む、推定方法。

【請求項2】

Ｎの値が１から５０までの間である、請求項１に記載の推定方法。

【請求項3】

Ｍの値が１から２０までの間である、請求項１に記載の推定方法。

【請求項4】

Ｒの値が１％から２０％までの間である、請求項１または２に記載の推定方法。

【請求項5】

Ｘの値が１から２０までの間である、請求項１に記載の推定方法。

【請求項6】

所定の分子の特性を決定する方法であって、一定範囲の分子構造を実験で決定するため、および分子構造の差をそれらの分子の選択された所定の作用の差と相関させるように請求項１から５のいずれか一項による推定方法を使用して得られる値とその範囲の値を比較するための複合イオン移動度−質量分析（ＩＭ−ＭＳ）技法の使用を含む、方法。

【請求項7】

分子構造の差を薬物の作用の差と相関させる、請求項６に記載の方法。

【請求項8】

前記薬物が有機金属ベースの薬物である、請求項７に記載の方法。

【請求項9】

前記有機金属ベースの薬物が抗癌剤である、請求項８に記載の方法。

【請求項10】

有機金属薬物が異性体ルテニウム（Ｒｕ）ベースである、請求項９に記載の方法。

【請求項11】

ＩＭ−ＭＳ技法が進行波（Ｔ波）移動度分離の使用を含む、請求項６から１０のいずれか一項に記載の方法。

【請求項12】

所定の分子の特性を決定するシステムであって、イオン移動度セルと、質量分析計と、一定範囲の分子構造を実験で決定するように、および分子構造の差をそれらの分子の選択された所定の作用の差と相関させるために請求項１から５のいずれか一項に記載の推定方法を使用して得られる値とその範囲の値を比較するようにプログラムされるまたは構成されるプロセッサとを含む、システム。

【請求項13】

請求項６から１１のいずれか一項に記載の方法を実行するようにさらに構成される、請求項１２に記載のシステム。

【請求項14】

請求項１から１１のいずれか一項に記載の方法を実現する手順をプロセッサに実行させるためのコンピュータ読出可能プログラムコード手段を含む、コンピュータプログラム要素。

【請求項15】

コンピュータ読出可能媒体上に具体化される、請求項１４に記載のコンピュータプログラム要素。

【請求項16】

記憶されたプログラムを有するコンピュータ読出可能媒体であって、プログラムが請求項１から１１のいずれか一項に記載の方法を実現する手順をコンピュータに実行させるためのものである、コンピュータ読出可能媒体。

【発明の詳細な説明】

【技術分野】

【0001】

本発明は、広範な分子種についてイオン移動度の予測で使用するために分子の断面積を推定する方法に関する。

【背景技術】

【0002】

イオン移動度は、気相衝突断面積（Ω）の差に基づき異性体種を分離する可能性を有し、理論モデルとの比較を通してイオン配座に関する価値のある情報を提供することができる。しかしながら、分子の衝突断面積から得られるイオン移動度を分子の構造に正確に関連付けるためには、構成原子の気相半径の予備知識を得る必要がある。

【0003】

分子断面は分子により投じられる「影」の回転平均面積として近似される。「影」という用語は、任意の平面上への分子、原子、または楕円体の投影を指すために使用される。利用される平面は計算の間中、固定されたままである。

【0004】

投影近似（ＰＡ）が、イオン移動度の予測で使用される分子の断面を計算する有用な方法であることが示された。投影近似断面を直接計算することはできないが、モンテカルロ（Ｍｏｎｔｅ Ｃａｒｌｏ）技法を使用して数値近似を得ることができる。モンテカルロ法は、計算アルゴリズムの結果を計算するために反復ランダムサンプリングに頼り、かつ決定論的アルゴリズムを使って正確な結果を計算することが実行できないまたは不可能であるときに使用されることが多い計算アルゴリズムの部類である。

【0005】

以下のステップ：
（１）Ａ＝Ｏ、ｎ＝Ｏに設定し、
（２）分子のランダムな配向を選択する
（３）分子の影の境界を定める長方形を選択し、長方形の面積＝ａｅｒａ Ａｂｏｘを計算する
（４）長方形内部の点をランダムに選択する
（５）点が分子の影の内部にある場合、Ａ＝Ａ＋Ａｂｏｘとする
（６）ｎ＝ｎ＋１とする
（７）ｎ＜ｎ_ｍａｘの場合、ＧＯＴＯ ２
（すなわち、次の反復で分子の別のランダムな配向を選択する）
（８）終了。結果はＡ／ｎである
を含むＭｏｂＣａｌ（ｈｔｔｐ：／／ｗｗｗ．ｉｎｄｉａｎａ．ｅｄｕ／〜ｎａｎｏ／Ｓｏｆｔｗａｒｅ．ｈｔｍｌ）として知られているアルゴリズムを用いて、分子の構成原子の気相半径の予備知識を得ることが知られている。

【発明の概要】

【発明が解決しようとする課題】

【0006】

この手法に伴い少なくとも２つの可能性のある（しかし関連する）問題がある。第一に、ｎ_ｍａｘの選択では、配向の数（ステップ２）が十分かを知ることが不可能でないにしても困難であり、第二に、単純な収束試験が誤りに導く可能性があるほどに収束が遅い状況がある。

【課題を解決するための手段】

【0007】

上記のアルゴリズムは、配向ごとに２つ以上の選択を可能にするように改善されることができ、改善されたアルゴリズムを与える。

【0008】

本発明の１つの態様または特徴が、たとえば所定のまたは試料の分子の特性を決定するための分子の断面積を推定する方法を提供し、方法は、以下のステップ：
（１）Ａ＝Ｏ、ｎ＝Ｏに設定し、
（２）分子のランダムな配向を選択する
（３）分子の影の境界を定める長方形を選択し、長方形の面積＝ａｅｒａ Ａｂｏｘを計算する
（４）長方形内部にＫ個の点をランダムに選択する
（５）分子の影の中にある各点に対して、Ａ＝Ａ＋Ａｂｏｘ／Ｋとする
（６）ｎ＝ｎ＋１とする
（７）ｎ＜ｎ_ｍａｘの場合、ＧＯＴＯ ２（次の反復）
（８）終了。結果はＡ／ｎである。

【0009】

上記の方法は、本明細書で以下に「長方形」法と呼ばれる。

【0010】

しかしながら、本発明はまた、新規のサンプリング技法を利用して、ユーザにより指定された精度まで得られる回転平均断面の収束に対して改善された信頼性を提供する投影近似の非常に効果的な実装である新しく開発された衝突断面積（ＣＣＳ）アルゴリズムの使用を含む。

【0011】

本発明の別の態様が、たとえば所定のまたは試料の分子の特性を決定するために、分子の断面積を推定する方法を提供し、方法は：
Ａ）当該分子の幾何構造を予測するステップと、
Ｂ）分子内部の各原子の断面積に対する値を指定するステップと、
Ｃ）分子内部の各原子に対する断面積の値の総和を計算するステップと、
Ｄ）当該の前記分子の予測された幾何構造の配向をランダムに選択するステップと、
Ｅ）原子の断面積に比例する確率で１つまたは複数の原子をランダムに選択するステップと、
Ｆ）選択された原子の断面積内部、または選択された原子のそれぞれの断面積内部の位置をランダムに選択するステップと、
Ｇ）断面積が前記ランダムに選択された位置を含む原子の数を識別するステップと、
Ｈ）前記ランダムに選択された配向での当該の前記分子の断面積に対する値を計算するステップと、
Ｉ）すべての前記ランダムに選択された配向全体の平均断面積を計算するステップと、
Ｊ）少なくともＮ回の反復の間、および／または平均断面積のＭ回の独立した計算値がＸ回の反復の間、所定の範囲Ｒ以内に一致するまで、ステップＤ）からステップＩ）を反復するステップとを含む。

【0012】

上記の方法は、本明細書で以下に「重要度サンプリング（ｉｍｐｏｒｔａｎｃｅ ｓａｍｐｌｉｎｇ）」法と呼ばれる。

【0013】

本発明の特徴によれば、値Ｎは０から５０までの間、たとえば０から３０までの間、または２０から５０までの間、または１０から４０までの間とすることができる。

【0014】

本発明の別の特徴によれば、値Ｍは１から２０までの間、たとえば１から１０までの間、または１から５までの間、または１０から２０までの間、または１５から２０までの間、または５から１５までの間とすることができる。

【0015】

本発明の別の特徴によれば、値Ｒは０．１％から２０％までの間、たとえば１％から２０％までの間、または０．１％から１％までの間とすることができる。［Ｒはパーセンテージまたは絶対値とすることができることに留意されたい］

【0016】

本発明の別の特徴によれば、値Ｘは１から２０までの間、たとえば１から１０までの間、または１から５までの間、または１０から２０までの間、または１５から２０までの間、または５から１５までの間とすることができる。

【0017】

本発明の別の態様が、一定範囲の分子構造を実験で決定し、分子構造の差をそれらの分子の選択された所定の作用の差と相関させるように上記のいずれかによる推定法により得られる値とその範囲の値を比較するための複合イオン移動度−質量分析（ＩＭ−ＭＳ）技法の使用を含む、所定の分子（たとえば小分子および小タンパク質、ペプチド、オリゴヌクレオチド、ならびにグルカン）の特性を決定する方法を提供する。

【0018】

本発明のこの態様の特徴によれば、方法は分子構造の差を所定の抗癌剤の抗癌作用の差と相関させることができる。好ましくは、抗癌剤は有機金属ベースの薬物である。好ましくは、有機金属薬物は異性体ルテニウム（Ｒｕ）ベースである。

【0019】

本発明のこの態様の別の特徴によれば、ＩＭ−ＭＳ技法は進行波（Ｔ波）移動度分離の使用を含むことができる。

【0020】

本発明の方法は、イオン移動度による構造的に関連する異性体の分離および特徴付けを提供するために、気相相互作用半径決定および断面アルゴリズム計算を与える。

【0021】

本発明の別の態様が、イオン移動度セルと、質量分析計と、一定範囲の分子構造を実験で決定するように、および分子構造の差をそれらの分子の選択された所定の作用の差と相関させるために上述された方法を評価することにより得られる値とその範囲の値を比較するようにプログラムされるまたは構成されるプロセッサとを含む、所定の分子の特性を決定するためのシステムを提供する。

【0022】

本発明のさらに別の態様が、上述された方法を実現するための手順をプロセッサに実行させるためのコンピュータ読出可能プログラムコード手段を含むコンピュータプログラム要素を提供する。

【0023】

コンピュータプログラム要素は、コンピュータ読出可能媒体上で具体化されることができる。

【0024】

本発明のさらに別の態様が、上述された方法を実現する手順をコンピュータに実行させるためのものである、コンピュータ読出可能媒体上に記憶されるプログラムを有するコンピュータ読出可能媒体を提供する。

【0025】

ここで、本発明の一実施形態が、添付の図面を参照して例示としてのみ説明される。

【図面の簡単な説明】

【0026】

【図1】本発明の一実施形態によるアルゴリズムを図示するフローチャートである。

【発明を実施するための形態】

【0027】

本発明の１つのアルゴリズムは、以下の通りである：
（１）Ａｍａｘ＝ｓｕｍ ｏｖｅｒ π^＊Ｒ_ｉ^＊Ｒ_ｉ、ここでＲ_ｉはｉ番目の原子の半径である
（２）Ａ＝０、ｎ＝０に設定し、
（３）次に、分子のランダムな配向を選択する
（４）原子の断面積に比例する確率でランダムにＫ個の原子を選択する
（５）次に、それぞれの選択された原子の円形の影の内部で一様な確率で点を選択する
（６）そのように選択された各点について、
Ａ＝Ａ＋Ａｍａｘ／（ｍ^＊Ｋ）
ここで、ｍは影が点に重なる原子の数である
（７）ｎ＝ｎ＋１
（８）ｎ＜ｎ_ｍａｘである場合、ＧＯＴＯ ３（次の反復）
（９）終了。結果はＡ／ｎである。

【0028】

多くの分子について、この重要度サンプリングアルゴリズムは、参照された従来技術のアルゴリズムの結果よりも速く所望の結果に収束する。

【0029】

上記のアルゴリズムと共に以下の終了手順が実装される：
（１）Ｎの投影近似計算値をインスタンス生成する（今後オブジェクトと呼ばれる）、
（２）これらのオブジェクトのそれぞれに対して順番に反復を実行し、各オブジェクトに対して面積推定を実行し続ける、
（３）すべてのオブジェクトに対する面積推定がユーザにより指定される（パーセンテージまたは絶対の）許容範囲と一致するとき終端し、これらの推定の全体平均を結果として報告する。

【0030】

終了した際に信頼性を達成するために、言及された従来技術アルゴリズムと本発明のアルゴリズムの両方を利用するオブジェクトの混合を使用することが有利となり得る。

【0031】

この新しく開発されたＣＣＳアルゴリズムは、新規のサンプリング技法を利用して、ユーザにより指定される精度まで得られる回転平均断面の収束に対する改善された信頼性を提供する投影近似の非常に効果的な実装である。

【0032】

質量分析計のＩＭＳ機器は、ヘモグロビンペプチド混合物、および知られている衝突断面積の別の化合物、たとえばポリグリシン（ｈｔｔｐ：／／ｗｗｗ．ｉｎｄｉａｎａ．ｅｄｕ／〜ｃｌｅｍｍｅｒ／Ｒｅｓｅａｒｃｈ／ｒｅｓｅａｒｃｈ．ｈｔｍ）を使用して較正され、次に、較正係数が得られた。各小分子のイオン移動度が測定され、したがって、実験的衝突断面積（ＣＣＳ）が決定された。次に、ＣＣＳアルゴリズムの出力が、分析された小分子のそれぞれの実験で得られたＣＣＳ値とぴったり一致するように、個々の小分子に特有の原子のそれぞれの半径がＣＣＳアルゴリズム内部で「調整」された。そうして、Ｃ、Ｈ、Ｏ、Ｓ、Ｃｌ、Ｆ、Ｎ、およびＲｕに対する最適相互作用半径値（Å^２）が得られた。

【0033】

ＣＣＳアルゴリズムを最適化するために使用された小分子は、Ｃ６０、Ｃ７０、ピレン、カンフェン、フェナントレン、トリフェニレン、ナフタレン、アンピシリン、スペルミン、ロラゼパム、カフェイン、レセルピン、ラフィノース、ニフェジピン、ニモジピン、デキサメタゾン、ジフェニルヒドラミン、メトプロロール、グルタチオン、システイン、アポモルヒネ、エリスロマイシン、オキシテトラサイクリン、ベラパミル、サルブタモール、アセチルサリチル酸、２０種類の天然の一般的なアミノ酸、アンジオテンシンＩＩ、オキシトシン、テストステロン、イブプロフェン、βシクロデキストリン、および化合物のｒｕ−ｏｒｔｏ（および関連するＨＣｌ中性損失種（ｎｅｕｔｒａｌ ｌｏｓｓ ｓｐｅｃｉｅｓ））、ｒｕ−ｍｅｔａ（および関連するＨＣｌ中性損失種）、ｒｕ−ｐａｒａ（および関連するＨＣｌ中性損失種）を含むルテニウムであった。

【0034】

低分子量異性体有機ルテニウム抗癌複合体の分離および配列分析のために分子モデリングと組み合わせられるイオン移動度−質量分析が、進行波技法に基づきイオン移動度を使用して分離され、実験で決定された衝突断面積が理論的計算値を比較された。

【0035】

ある種の有機金属薬物が癌との戦いで使用されている。たとえば、プラチナベースの薬物、シスプラチン［シス−ジアミンジクロリド白金（ＩＩ）］がこの範ちゅうに入り、癌との戦いで使用される主要な薬物の１つである。ＤＮＡは多くの金属ベースの抗癌剤にとって可能性のある標的であり、ＤＮＡ構造の歪みが抗癌作用と相関関係にあることが多い。別の金属ベースの抗癌剤、たとえばルテニウムベースの抗癌剤が、癌と闘うための代替治療薬として開発されてきている。具体的には、抗癌作用の範囲を広げること、望まれない副作用を低減すること、およびシスプラチンおよび関連する薬物との交差耐性を回避することが目的である。これらの新規のルテニウムベースの薬物の物理的なサイズおよび形状の洞察が、構造−作用関係の解明のため、および重要な相互作用たとえばＤＮＡの中への挿入を最適化するために重要である。３つの新規の異性体ルテニウムベース抗癌剤が、本発明による推定法と共に、複合イオン移動度および質量分析（ＩＭ−ＭＳ）手法を使用して説明された。

【0036】

独立の技法として、ＭＳは異性体種を分離できない、またはバルク構造配座情報を提供できない。しかしながら、ＩＭは、気相中の異性種の衝突断面積（ＣＣＳ；物理的なサイズおよび形状）の差に基づき異性体種を（ＭＳ取得の時間尺度で）迅速に分離する能力を有し、したがって、イオン構成に関する特有の情報を提供する。ＩＭとＭＳとの組合せが非常に強力な分析手段を提供する。

【0037】

ＩＭ−ＭＳ技法のあり得る利点を精査するために、「進行波」（Ｔ波）移動度分離に関して基づく器具が、３つの低分子量異性体ルテニウムテルフェニル抗癌複合体（^１０２Ｒｕに基づきｍ／ｚが４２７．１）の混合物を分析するために使用された。形状の個々の差が、形状の個々の差を抗癌作用の差と相関させる取り組みで研究された。分子モデリングも一定範囲の可能な構造を生成するために使用され、これらの構造に対する理論的ＣＣＳが、実験で得られたＴ波値と比較するために計算された。実験で得られたＣＣＳ測定値と理論的に得られたＣＣＳ測定値の間に良好な一致が観測された。

【0038】

オルト−、メタ−、またはパラ−テルフェニルアレーンリガンド（Ｍｗ４２７．１）およびその後の異性体ＨＣｌ中性損失種（Ｍｗ３９１．１）を含む３つの低分子量異性体有機ルテニウム抗癌複合体の分離および配座分析のために、新しいＣＣＳアルゴリズム、分子モデリング、およびイオン移動度質量分析と組み合わせたイオン移動度により得られる気相相互作用半径が本発明に従って使用された。６つの異性体複合体が、実験衝突断面積が計算され理論値と比較されることができるようにするよく規定された到来時間分布を示した。すべての実験で得られた結果とすべての理論的に得られた結果の間の良好な一致が観測された。より緻密な形状と比較した細長いテルフェニルアレーンの形状の差が、細長い構造の著しくより高い抗癌作用に対する重要な寄与を引き起こすようである。

【0039】

要約すると、本発明は、Ｈ、Ｃ、ならびにこれまで精査されなかったＯ、Ｓ、Ｃｌ、Ｆ、Ｎ、Ｆｅ、Ｐｔ、およびＲｕに対する推論される気相相互作用半径を報告している。さらに、本発明の方法は、ＣＣＳが９．０Å^２未満だけ異なる異性体ルテニウム（ＩＩ）テルフェニル抗癌複合体が分離され、進行波イオン移動度により特徴付けられることができることを初めて実証している。

【0040】

下記の表は、長方形法および重要度サンプリング法と比較された、Ｍｏｂｃａｌアルゴリズムを使用して得られた結果を示す。

【0041】

【表1】

【0042】

上述の特徴および／または添付図面に示される特徴の任意の数の組合せが、従来技術に対して明らかに有利な点を提供し、したがって、本明細書で説明される本発明の範囲に入ることが当業者により理解される。

【図1】