WO2014115339A1 - 溶液中の生体高分子の構造揺らぎとダイナミクスを記述する方程式 - Google Patents
溶液中の生体高分子の構造揺らぎとダイナミクスを記述する方程式 Download PDFInfo
- Publication number
- WO2014115339A1 WO2014115339A1 PCT/JP2013/051940 JP2013051940W WO2014115339A1 WO 2014115339 A1 WO2014115339 A1 WO 2014115339A1 JP 2013051940 W JP2013051940 W JP 2013051940W WO 2014115339 A1 WO2014115339 A1 WO 2014115339A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- fluctuations
- dynamics
- biopolymer
- structural
- equation
- Prior art date
Links
Images
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C10/00—Computational theoretical chemistry, i.e. ICT specially adapted for theoretical aspects of quantum chemistry, molecular mechanics, molecular dynamics or the like
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
- G16B15/00—ICT specially adapted for analysing two-dimensional or three-dimensional molecular structures, e.g. structural or functional relations or structure alignment
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
- G16B15/00—ICT specially adapted for analysing two-dimensional or three-dimensional molecular structures, e.g. structural or functional relations or structure alignment
- G16B15/30—Drug targeting using structural data; Docking or binding prediction
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
- G16C20/50—Molecular design, e.g. of drugs
Definitions
- the present invention relates to an equation that describes structural fluctuation and dynamics of a biopolymer conjugated with density (concentration) fluctuation of a solution at an atomic level, and a calculation algorithm thereof.
- molecular simulation molecular dynamics method or Monte Carlo method
- molecular dynamics method Monte Carlo method
- Monte Carlo method is a very effective method for describing relatively small structural fluctuations in proteins, but it is a group conjugated with solvent fluctuations.
- a particularly fatal flaw in molecular simulation is that it is almost powerless for systems that contain ions other than water and other small molecules in the solvent. This is because the time required for the diffusion process of ions in the solution is much longer than the integration time step of the Newton equation, and the simulation does not converge in a finite time.
- the deficiencies in conventional theories and models describing the structural fluctuations and dynamics of biopolymers stem from the fact that the liquid structure and density fluctuations of the solvent are not properly taken into account.
- the problem to be solved by the present invention is exactly this point, and is a method of describing the structural fluctuation and dynamics of a biopolymer in consideration of the liquid structure and density fluctuation of the solvent at the atomic level.
- ⁇ R ⁇ (t) Displacement (fluctuation) of the atomic coordinates from the equilibrium position in the biopolymer
- M ⁇ Mass of atom in biopolymer
- ⁇ ⁇ Friction force acting on atom of biopolymer
- W ⁇ (t) Random force acting on atom of biopolymer
- Equation (2) is an equation describing the density fluctuation and dynamics of the solvent around the biopolymer.
- Equation (1) on the gist of the invention has the form of a generalized Langevin equation, but if we ignore the second term (friction term) and the third term (random force) on the right side of this equation, coupled vibration Has a similar form to the child. Then, the first term on the right-hand side has the physical interpretation of “restoring force” proportional to the displacement ( ⁇ R) from the equilibrium structure of the atoms of the biopolymer, and the proportional constant is “force constant” in the coupled oscillator.
- restoring force proportional to the displacement ( ⁇ R) from the equilibrium structure of the atoms of the biopolymer
- the dynamics of the biopolymer described by this equation are fluctuations around its equilibrium structure (the minimum point of the free energy surface), and the proportionality constant described above is the second derivative (at the moment) of the atomic displacement of this free energy surface ( In the case of a coupled oscillator, it has the physical meaning of (second derivative of potential energy between atoms).
- ⁇ R ⁇ which is an argument of F ( ⁇ R ⁇ ) in the equation (4) represents a deviation from the equilibrium structure of all atomic coordinates of the biopolymer.
- the free energy curved surface (F ( ⁇ R ⁇ )) of the biopolymer in the solution may be obtained.
- F ( ⁇ R ⁇ ) is an amount defined as the sum of the interaction energy between protein atoms and the solvent free energy, but the applicant has already proposed F ( ⁇ R ⁇ ) based on the 3D-RISM / RISM theory.
- a numerical algorithm for obtaining the dispersion / covariance matrix of the structural fluctuations of biopolymers is an amount defined as the sum of the interaction energy between protein atoms and the solvent free energy, but the applicant has already proposed F ( ⁇ R ⁇ ) based on the 3D-RISM / RISM theory.
- f ⁇ is a perturbation accompanying ligand binding
- ⁇ R ⁇ > 1 is a protein structural change (displacement of each atom) induced by the perturbation
- its response function is a dispersion / covariance matrix ⁇ R ⁇ ⁇ R ⁇ > 0 .
- this dispersion / covariance matrix is obtained by molecular dynamics simulation.
- the method cannot describe a large protein fluctuation (collective motion mode).
- a system containing ions or other small molecules in a solvent cannot be handled.
- the variance / covariance matrix can be obtained by performing second order differentiation on the free energy curved surface of the protein with respect to its atomic coordinates. Since the free energy curved surface of a protein can be easily obtained by the 3D-RISM / RISM method, which the applicant has an inherent copyright, the problem of inductive binding of a ligand to a protein can be solved.
- FIG. 1 shows a calculation flow chart regarding “induced binding” of a ligand by the above-described protein.
- FIG. 1 shows the flow of numerical calculation in the center column, and the flow of input / output data in the left and right columns.
- the arrows between the center and the left and right columns indicate the direction of data flow between the numerical calculation unit and the data storage unit.
- the correlation function ( ⁇ ⁇ (r, r ′)) of the density fluctuation in the solvent is obtained by solving the RISM equation with the interaction information between the solvent molecules and the geometric shape information of the solvent molecules as inputs.
- Step 2 By solving the 3D-RISM equation using the protein structure information (amino acid sequence, coordinate data, interaction parameters, etc.) and ⁇ ⁇ (r, r ′) obtained in (Step 1) as input, The solvent molecule distribution (h ⁇ (r)) and free energy (F ( ⁇ R ⁇ )) are determined. (Step 3) The second-order differentiation of the free energy (F ( ⁇ R ⁇ )) obtained in (Step 2) with respect to the position coordinates of the protein is performed to obtain the variance / covariance matrix ( ⁇ R ⁇ R>) of the structural fluctuation. (Step 4) Using the perturbation force f associated with drug binding and ⁇ R ⁇ R> obtained in (Step 3) as inputs, the structural change of the protein associated with the perturbation is obtained by linear response theory.
- Target Search The first step in developing a new drug begins with the search for proteins and DNA that are the targets of the drug. Currently, this process is based exclusively on physiological and biochemical experimental studies, and the eye of the analysis is to examine the function of biomolecules involved in disease in some way. However, all functions of biomolecules (catalyst, information transmission, etc.) involve molecular recognition processes (ligand binding to biomolecules) as elementary processes. Therefore, the most important task in drug target search is nothing other than analyzing the molecular recognition process, that is, the binding affinity between the target molecule and the substrate molecule (ligand).
- the 3D-RISM / RISM theory which is one of the basic methods of this application, provides a powerful weapon as a method for determining the binding affinity of a ligand for a target protein with a fixed structure. However, this is not applicable when structural fluctuations are essential.
- the invention of this application provides a new calculation method that can be applied to the process in which target protein structure fluctuation plays an essential role.
- Screening drug candidate compounds Another important process in the development of new drugs is the search for compounds that most efficiently inhibit or promote the function of target proteins and DNA. In this process as well, it is necessary to analyze the molecular recognition process of the candidate compound by the target protein.
- the method of the present invention is a powerful method for analyzing the binding affinity of a compound to a target for which structural fluctuation is essential.
Landscapes
- Spectroscopy & Molecular Physics (AREA)
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Theoretical Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Chemical & Material Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Bioinformatics & Computational Biology (AREA)
- Medical Informatics (AREA)
- Crystallography & Structural Chemistry (AREA)
- Biotechnology (AREA)
- Evolutionary Biology (AREA)
- Biophysics (AREA)
- Computing Systems (AREA)
- Pharmacology & Pharmacy (AREA)
- Medicinal Chemistry (AREA)
- Information Retrieval, Db Structures And Fs Structures Therefor (AREA)
- Organic Low-Molecular-Weight Compounds And Preparation Thereof (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
申請者は液体の平衡・非平衡統計力学に基づき、本申請の請求項１において示された２つの理論式（（１）生体高分子の揺らぎとダイナミクス、（２）溶媒の密度揺らぎおよびダイナミクス）を導出し、これらの方程式の解析から、溶液内の生体高分子の自由エネルギー曲面をその原子座標に関して２階微分を行うことで、生体高分子の構造揺らぎの分散・共分散行列を求める新しい理論的手法を考案した。 本発明は、イオンや他の小分子を含む複雑な水溶液系の液体構造や密度揺らぎを原子レベルで解析し、それらが生体高分子の構造揺らぎとダイナミクスに対して及ぼす影響を原子レベルで、高精度に予測、記述することを可能にする。 本発明の応用により、高親和性リガンドの探索、医薬品候補化合物デザインならびに新規機能性材料の設計、予測及び開発における知的分子設計手法として、医療およびナノテクノロジー分野の発展に大きく寄与しうる。
Description
本発明は、溶液の密度（濃度）揺らぎと共役した生体高分子の構造揺らぎとダイナミクスを原子レベルで記述する方程式、およびその計算アルゴリズムに関するものである。
溶液中の生体高分子の構造揺らぎやダイナミクスはその機能と密接に関わっており、従来から多くの理論的ないしは計算科学的方法論が提案されてきた。その理由は、生体高分子（特に、蛋白質）の構造揺らぎやダイナミクスが酵素活性などその機能を支配する要因となっているからである。また、現在、人類社会にとって焦眉の課題となっている新薬の創成にとっても、生体高分子の構造揺らぎやダイナミクスを記述する方法論の開発は実践的（実用的）な要請となっている。それは薬剤化合物が標的分子（蛋白質やＤＮＡ）に結合するプロセスにおいて、生体高分子の構造揺らぎやダイナミクスが本質的役割を演じるからである。
生体高分子の構造揺らぎやダイナミクスに関して最も初期に提案された方法は蛋白質の基準振動解析に基づく手法であるが、この手法は溶媒である水を含まない「真空中」の蛋白質に関するものであった。このため、その「構造揺らぎ」はいわば「非生命的」な揺らぎに過ぎず、生体内にある蛋白質の「生命的」な揺らぎを記述することはできない。一方、溶媒を現象論的に考慮したいくつかの手法（Ｐｏｉｓｓｏｎ−Ｂｏｌｔｚｍａｎｎ方程式、Ｌａｎｇｅｖｉｎ方程式など）が提案されているが、それらは分子レベルでの溶媒の構造を無視しているため、蛋白質の活性部位に対するリガンドの結合に伴う水和、脱水和など、蛋白質の構造揺らぎに本質的影響を与える因子を記述することができない。さらに、最近、普及しているいわゆる分子シミュレーション（分子動力学法やモンテカルロ法）は、蛋白質の比較的小さな構造揺らぎを記述する上で非常に有効な方法であるが、溶媒の揺らぎと共役した集団運動モードを求める上では、系のサイズと揺らぎの時間スケールの点で大きな技術的困難を免れない。分子シミュレーションの特に致命的な欠陥は、溶媒中に水以外のイオンや他の小分子を含む系に対して、ほとんど無力だということである。それは、溶液中のイオンなどの拡散過程に要する時間はニュートン方程式の積分時間ステップに比べて非常に長く、シミュレーションが有限の時間で収束しないからである。
生体高分子の構造揺らぎやダイナミクスに関して最も初期に提案された方法は蛋白質の基準振動解析に基づく手法であるが、この手法は溶媒である水を含まない「真空中」の蛋白質に関するものであった。このため、その「構造揺らぎ」はいわば「非生命的」な揺らぎに過ぎず、生体内にある蛋白質の「生命的」な揺らぎを記述することはできない。一方、溶媒を現象論的に考慮したいくつかの手法（Ｐｏｉｓｓｏｎ−Ｂｏｌｔｚｍａｎｎ方程式、Ｌａｎｇｅｖｉｎ方程式など）が提案されているが、それらは分子レベルでの溶媒の構造を無視しているため、蛋白質の活性部位に対するリガンドの結合に伴う水和、脱水和など、蛋白質の構造揺らぎに本質的影響を与える因子を記述することができない。さらに、最近、普及しているいわゆる分子シミュレーション（分子動力学法やモンテカルロ法）は、蛋白質の比較的小さな構造揺らぎを記述する上で非常に有効な方法であるが、溶媒の揺らぎと共役した集団運動モードを求める上では、系のサイズと揺らぎの時間スケールの点で大きな技術的困難を免れない。分子シミュレーションの特に致命的な欠陥は、溶媒中に水以外のイオンや他の小分子を含む系に対して、ほとんど無力だということである。それは、溶液中のイオンなどの拡散過程に要する時間はニュートン方程式の積分時間ステップに比べて非常に長く、シミュレーションが有限の時間で収束しないからである。
Ｈｉｒａｔａ Ｆ．Ｍｏｌｅｃｕｌａｒ Ｔｈｅｏｒｙ ｏｆ Ｓｏｌｖａｔｉｏｎ．（Ｄｏｒｄｒｅｃｈｔ：Ｋｌｕｗｅｒ Ａｃａｄｅｍｉｃ Ｐｕｂｌｉｓｈｅｒ），２００３
"新しい分子統計力学・統計的動力学記述法"、"巨大分子系の計算化学"第７章、８８−９９ページ、日本化学会、２０１２年３月３０日
"分子認識の統計力学と生体機能"、生物物理５１（５）、２２２−２２５ページ、日本生物物理学会、２０１１年１１月
生体高分子の構造揺らぎとダイナミクスを記述する従来の理論やモデルの欠陥は溶媒の液体構造や密度揺らぎを正しく考慮していないことに起因している。本発明が解決しようとする課題は、まさに、この点であり、溶媒の液体構造や密度揺らぎを原子レベルで考慮して生体高分子の構造揺らぎとダイナミクスを記述する方法である。
上記の課題を解決するための手段は液体の平衡統計力学理論（３Ｄ−ＲＩＳＭ／ＲＩＳＭ）および非平衡統計力学理論（一般化ランジェヴァン理論）であり、これらを結合することにより溶液中の生体高分子の揺らぎとダイナミクスを記述する新規の方程式を導出した。
上の（１）式は生体高分子の構造揺らぎとダイナミクスを記述する式である。式中の各記号は以下の物理的意味をもっている。
ΔＲα（ｔ）：生体高分子中の原子座標の平衡位置からの変位（揺らぎ）
Ｍα：生体高分子中の原子の質量
Γαβ：生体高分子の原子に働く摩擦力
Ｗα（ｔ）：生体高分子の原子に働くランダムな力
（１）式の右辺第一項は、変位（原子座標の揺らぎ）に比例する復元力を表しており、その係数に含まれる行列Ｌが次の式で表されることを発見した。
（３）式中のΔＲは平衡状態での生体高分子の原子座標の平衡位置からのずれ（揺らぎ）を表しており、したがって、行列Ｌは構造揺らぎの分散・共分散行列を表す。
（２）式は生体高分子の周りの溶媒の密度揺らぎとダイナミクスを記述する方程式である。式中の各記号は以下の物理的意味をもっている。
発明の要点
上の方程式（１）は一般化ランジェヴァン方程式の形式を持っているが、この方程式の右辺第２項（摩擦項）と第３項（ランダムな力）を無視すると、連成振動子と類似の形式をもつ。そうすると、右辺第１項は生体高分子の原子の平衡構造からの変位（ΔＲ）に比例する「復元力」という物理的解釈が成り立ち、その比例定数に関しては連成振動子における「力の定数」（Ｈｅｓｓｉａｎ）との類推が成立する。一方、この方程式で記述する生体高分子のダイナミクスはその平衡構造（自由エネルギー曲面の最小点）の周りの揺らぎであり、上に述べた比例定数はこの自由エネルギー曲面の原子変位に関する２次微分（連成振動子の場合は原子間のポテンシャルエネルギーの２次微分）という物理的意味をもつ。そこで、我々は下記の等式が成り立つことを発見した。
（４）式中のＦ（｛ΔＲ｝）の引数である｛ΔＲ｝は生体高分子のすべての原子座標の平衡構造からずれを表している。上式によれば、蛋白質の構造揺らぎの分散・共分散行列を得るためには、溶液内の生体高分子の自由エネルギー曲面（Ｆ（｛ΔＲ｝））を求めれば良いことになる。（Ｆ（｛ΔＲ｝））は蛋白質原子間の相互作用エネルギーと溶媒自由エネルギーの和として定義される量であるが、当申請者らはすでに３Ｄ−ＲＩＳＭ／ＲＩＳＭ理論に基づきＦ（｛ΔＲ｝）を求める方法を発表しており、生体高分子の構造揺らぎの分散・共分散行列を求める数値計算アルゴリズムも確立している。
ΔＲα（ｔ）：生体高分子中の原子座標の平衡位置からの変位（揺らぎ）
Ｍα：生体高分子中の原子の質量
Γαβ：生体高分子の原子に働く摩擦力
Ｗα（ｔ）：生体高分子の原子に働くランダムな力
（１）式の右辺第一項は、変位（原子座標の揺らぎ）に比例する復元力を表しており、その係数に含まれる行列Ｌが次の式で表されることを発見した。
（２）式は生体高分子の周りの溶媒の密度揺らぎとダイナミクスを記述する方程式である。式中の各記号は以下の物理的意味をもっている。
上の方程式（１）は一般化ランジェヴァン方程式の形式を持っているが、この方程式の右辺第２項（摩擦項）と第３項（ランダムな力）を無視すると、連成振動子と類似の形式をもつ。そうすると、右辺第１項は生体高分子の原子の平衡構造からの変位（ΔＲ）に比例する「復元力」という物理的解釈が成り立ち、その比例定数に関しては連成振動子における「力の定数」（Ｈｅｓｓｉａｎ）との類推が成立する。一方、この方程式で記述する生体高分子のダイナミクスはその平衡構造（自由エネルギー曲面の最小点）の周りの揺らぎであり、上に述べた比例定数はこの自由エネルギー曲面の原子変位に関する２次微分（連成振動子の場合は原子間のポテンシャルエネルギーの２次微分）という物理的意味をもつ。そこで、我々は下記の等式が成り立つことを発見した。
本申請の発明がその効果を発揮すると予想される問題のひとつは、蛋白質に対するリガンド（薬剤分子）の誘起結合（ｉｎｄｕｃｅｄ ｆｉｔｔｉｎｇ）である。蛋白質に対するリガンドの結合親和力は、第一義的には結合前後の自由エネルギー差によって決定される。しかしながら、多くの蛋白質で見られる現象として、蛋白質に対するリガンドの結合部位の開口部を構成するアミノ酸残基の揺らぎ（開口部の開閉運動）がリガンドの結合親和力に大きな影響を及ぼすことが知られている。このような問題を取り扱う上で有用な下記の式（線形応答理論）がすでに横浜市立大学の池口らによって提案されている。
この式において、ｆβはリガンド結合に伴う摂動、〈ΔＲα〉１は摂動によって誘起される蛋白質の構造変化（各原子の変位）であり、その応答関数が分散・共分散行列〈ΔＲαΔＲβ〉０である。従来の方法ではこの分散・共分散行列を分子動力学シミュレーションに求めているが、上に述べたように、その方法では蛋白質の大きな揺らぎ（集団運動モード）を記述することができない。また、イオンや他の小分子を溶媒中に含む系を取り扱うことができない。一方、本申請の方程式（４）によれば、蛋白質の自由エネルギー曲面をその原子座標に関して２階微分を行うことにより、分散・共分散行列を求めることができる。蛋白質の自由エネルギー曲面は本申請者が固有の著作権を有する３Ｄ−ＲＩＳＭ／ＲＩＳＭ法により容易に求めることができるため、蛋白質に対するリガンドの誘起結合の問題が解けることになる。
本発明を実施する形態は計算科学であり、本申請で提案した方程式をコンピュータによって数値的に解法する。一例として、上に述べた蛋白質によるリガンドの「誘起結合」に関する計算のフロー図を図１に示す。
図１は、中央の列に数値計算の流れ、また左右の列には入出力データの流れを示してある。中央と左右の列の間にある矢印は数値計算部とデータ記憶部の間のデータの流れの方向を表す。
（ステップ１）
溶媒分子間の相互作用情報と溶媒分子の幾何形状情報を入力として、ＲＩＳＭ方程式を解法することにより、溶媒中の密度揺らぎの相関関数（χαγ（ｒ，ｒ’））を求める。
（ステップ２）
蛋白質の構造情報（アミノ酸配列、座標データ、相互作用パラメタなど）、および（ステップ１）で求めたχαγ（ｒ，ｒ’）を入力として、３Ｄ−ＲＩＳＭ方程式を解法することにより、蛋白質周りの溶媒分子の分布（ｈγ（ｒ））および自由エネルギー（Ｆ（｛ΔＲ｝））を求める。
（ステップ３）
（ステップ２）で求めた自由エネルギー（Ｆ（｛ΔＲ｝））を蛋白質の位置座標に関して２階微分を行うことにより、構造揺らぎの分散・共分散行列（〈ΔＲΔＲ〉）を求める。
（ステップ４）
ドラッグの結合に伴う摂動力ｆ、および（ステップ３）で求めた〈ΔＲΔＲ〉を入力として、線形応答理論により、摂動に伴う蛋白質の構造変化を求める。
（ステップ１）
溶媒分子間の相互作用情報と溶媒分子の幾何形状情報を入力として、ＲＩＳＭ方程式を解法することにより、溶媒中の密度揺らぎの相関関数（χαγ（ｒ，ｒ’））を求める。
（ステップ２）
蛋白質の構造情報（アミノ酸配列、座標データ、相互作用パラメタなど）、および（ステップ１）で求めたχαγ（ｒ，ｒ’）を入力として、３Ｄ−ＲＩＳＭ方程式を解法することにより、蛋白質周りの溶媒分子の分布（ｈγ（ｒ））および自由エネルギー（Ｆ（｛ΔＲ｝））を求める。
（ステップ３）
（ステップ２）で求めた自由エネルギー（Ｆ（｛ΔＲ｝））を蛋白質の位置座標に関して２階微分を行うことにより、構造揺らぎの分散・共分散行列（〈ΔＲΔＲ〉）を求める。
（ステップ４）
ドラッグの結合に伴う摂動力ｆ、および（ステップ３）で求めた〈ΔＲΔＲ〉を入力として、線形応答理論により、摂動に伴う蛋白質の構造変化を求める。
本発明の産業上の利用可能性として、最も有望な対象は新薬開発における「標的探索」と「薬剤候補化合物のスクリーニング」である。
標的探索：新薬開発のまず第一歩は薬剤の標的となる蛋白質やＤＮＡの探索から始まる。現在、このプロセスはもっぱら生理学的および生化学的実験研究に基づいて行われており、その解析の眼目は何らかの形で疾病に関わる生体分子の機能を調べることである。ところが、生体分子のすべての機能（触媒、情報伝達、など）には、その素過程として、分子認識過程（生体分子に対するリガンドの結合）が関わっている。したがって、薬剤の標的探索における最も重要な課題は、この分子認識過程、すなわち、標的分子と基質分子（リガンド）の間の結合親和力を解析することに他ならない。本申請の基盤となる方法のひとつである３Ｄ−ＲＩＳＭ／ＲＩＳＭ理論は構造が固定されている標的蛋白質に対するリガンドの結合親和力を求める方法として強力な武器を提供しているが、誘起結合過程のように構造の揺らぎが本質的である場合には適用できない。本申請の発明は標的蛋白質構造揺らぎが本質的な役割を演じる過程に対して適用できる新しい計算手法を提供する。
薬剤候補化合物のスクリーニング：新薬の開発におけるもうひとつの重要なプロセスは標的となる蛋白質やＤＮＡの機能を最も効率よく阻害あるいは促進する化合物の探索である。このプロセスにおいても、標的蛋白質による候補化合物の分子認識過程を解析する必要がある。この分子認識過程を第一義的に決めるのは標的蛋白質が候補化合物を結合する際の自由エネルギー変化であるが、同時に、蛋白質の構造揺らぎが決定的役割を演じることも少なくない。本発明の方法は、構造揺らぎが本質的であるような標的に対する化合物の結合親和力を解析する有力な方法である。
標的探索：新薬開発のまず第一歩は薬剤の標的となる蛋白質やＤＮＡの探索から始まる。現在、このプロセスはもっぱら生理学的および生化学的実験研究に基づいて行われており、その解析の眼目は何らかの形で疾病に関わる生体分子の機能を調べることである。ところが、生体分子のすべての機能（触媒、情報伝達、など）には、その素過程として、分子認識過程（生体分子に対するリガンドの結合）が関わっている。したがって、薬剤の標的探索における最も重要な課題は、この分子認識過程、すなわち、標的分子と基質分子（リガンド）の間の結合親和力を解析することに他ならない。本申請の基盤となる方法のひとつである３Ｄ−ＲＩＳＭ／ＲＩＳＭ理論は構造が固定されている標的蛋白質に対するリガンドの結合親和力を求める方法として強力な武器を提供しているが、誘起結合過程のように構造の揺らぎが本質的である場合には適用できない。本申請の発明は標的蛋白質構造揺らぎが本質的な役割を演じる過程に対して適用できる新しい計算手法を提供する。
薬剤候補化合物のスクリーニング：新薬の開発におけるもうひとつの重要なプロセスは標的となる蛋白質やＤＮＡの機能を最も効率よく阻害あるいは促進する化合物の探索である。このプロセスにおいても、標的蛋白質による候補化合物の分子認識過程を解析する必要がある。この分子認識過程を第一義的に決めるのは標的蛋白質が候補化合物を結合する際の自由エネルギー変化であるが、同時に、蛋白質の構造揺らぎが決定的役割を演じることも少なくない。本発明の方法は、構造揺らぎが本質的であるような標的に対する化合物の結合親和力を解析する有力な方法である。
Claims (2)
- 液体の統計力学理論（３Ｄ−ＲＩＳＭ／ＲＩＳＭ）と非平衡統計力学（一般化ランジェヴァン方程式）を組み合わせることによって得られる溶液中の生体高分子の構造揺らぎ（およびダイナミクス）と溶液の密度揺らぎ（およびダイナミクス）を関連づける下記の三本の理論式。
（１）生体高分子の揺らぎとダイナミクス：
（２）溶液の密度揺らぎとダイナミクス：
（３）生体高分子の構造揺らぎの分散・共分散行列をその自由エネルギーから求める式
- 上記請求項１で示された３つの理論式を溶液中の生体高分子の揺らぎおよびダイナミクスに具体化する計算アルゴリズムおよび計算プログラム。
Priority Applications (4)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| PCT/JP2013/051940 WO2014115339A1 (ja) | 2013-01-23 | 2013-01-23 | 溶液中の生体高分子の構造揺らぎとダイナミクスを記述する方程式 |
| JP2014509536A JPWO2014115416A1 (ja) | 2013-01-23 | 2013-11-25 | 溶液中の生体高分子の構造揺らぎとダイナミクスを記述する方程式 |
| PCT/JP2013/081601 WO2014115416A1 (ja) | 2013-01-23 | 2013-11-25 | 溶液中の生体高分子の構造揺らぎとダイナミクスを記述する方程式 |
| US14/762,918 US20150356276A1 (en) | 2013-01-23 | 2013-11-25 | Equations describing the structural fluctuation and dynamics of biopolymers in solution |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| PCT/JP2013/051940 WO2014115339A1 (ja) | 2013-01-23 | 2013-01-23 | 溶液中の生体高分子の構造揺らぎとダイナミクスを記述する方程式 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| WO2014115339A1 true WO2014115339A1 (ja) | 2014-07-31 |
Family
ID=51227152
Family Applications (2)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| PCT/JP2013/051940 WO2014115339A1 (ja) | 2013-01-23 | 2013-01-23 | 溶液中の生体高分子の構造揺らぎとダイナミクスを記述する方程式 |
| PCT/JP2013/081601 WO2014115416A1 (ja) | 2013-01-23 | 2013-11-25 | 溶液中の生体高分子の構造揺らぎとダイナミクスを記述する方程式 |
Family Applications After (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| PCT/JP2013/081601 WO2014115416A1 (ja) | 2013-01-23 | 2013-11-25 | 溶液中の生体高分子の構造揺らぎとダイナミクスを記述する方程式 |
Country Status (3)
| Country | Link |
|---|---|
| US (1) | US20150356276A1 (ja) |
| JP (1) | JPWO2014115416A1 (ja) |
| WO (2) | WO2014115339A1 (ja) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110793982A (zh) * | 2019-11-21 | 2020-02-14 | 山东建筑大学 | 一种纳米晶化动力学过程的高能x射线表征方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20170212982A1 (en) * | 2014-06-25 | 2017-07-27 | Japan Science And Technology Agency | Thermostabilized mutant-predicting apparatus for membrane protein, a thermostabilized mutant-predicting method, and computer program product |
-
2013
- 2013-01-23 WO PCT/JP2013/051940 patent/WO2014115339A1/ja active Application Filing
- 2013-11-25 US US14/762,918 patent/US20150356276A1/en not_active Abandoned
- 2013-11-25 WO PCT/JP2013/081601 patent/WO2014115416A1/ja active Application Filing
- 2013-11-25 JP JP2014509536A patent/JPWO2014115416A1/ja active Pending
Non-Patent Citations (4)
| Title |
|---|
| B.KIM ET AL.: "An attempt toward the generalized Langevin dynamics simulation", CONDENSED MATTER PHYSICS, vol. 11, no. 1, 2008, pages 179 - 190 * |
| BONGSOO KIM ET AL.: "Structural Fluctuation of Protein in Water around Its Native State: A New Statistical Mechanics Formulation, cond-mat.stat-mech", 30 October 2012 (2012-10-30), pages 1 - 22, Retrieved from the Internet <URL:http://arxiv.org/pdf/1210.7892v1.pdf> [retrieved on 20130329] * |
| FUMIO HIRATA ET AL.: "Suichu no Shudan Reiki to Ion no Dynamics: Ekitai no Sogo Saiyoten Model ni Motozuku Kenkyu", ENSEMBLE BUNSHI SIMULATION KENKYUKAI KAISHI, vol. 8, no. 3, 31 July 2006 (2006-07-31), pages 27 - 36 * |
| HIROHITO TADANO ET AL.: "Application and Performance Evaluation of the Volumetric Parallel 3D-FFT to 3D-RISM on Massively Parallel Cluster", IPSJ SIG NOTES, HEISEI 21 NENDO ?4?, 15 December 2009 (2009-12-15) * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110793982A (zh) * | 2019-11-21 | 2020-02-14 | 山东建筑大学 | 一种纳米晶化动力学过程的高能x射线表征方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| US20150356276A1 (en) | 2015-12-10 |
| JPWO2014115416A1 (ja) | 2017-01-26 |
| WO2014115416A1 (ja) | 2014-07-31 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Bernetti et al. | Pressure control using stochastic cell rescaling | |
| Shukla et al. | Molecular dynamics simulation of protein and protein–ligand complexes | |
| Monticelli et al. | Force fields for classical molecular dynamics | |
| Guvench et al. | Comparison of protein force fields for molecular dynamics simulations | |
| Meller | Molecular dynamics | |
| Rudzinski et al. | The role of many-body correlations in determining potentials for coarse-grained models of equilibrium structure | |
| Lindert et al. | Accelerated molecular dynamics simulations with the AMOEBA polarizable force field on graphics processing units | |
| de Ruiter et al. | Efficient and accurate free energy calculations on trypsin inhibitors | |
| Ruiter et al. | Extended thermodynamic integration: Efficient prediction of lambda derivatives at nonsimulated points | |
| Krämer et al. | Interactions of water and alkanes: Modifying additive force fields to account for polarization effects | |
| Kondo et al. | Free-energy landscapes of protein domain movements upon ligand binding | |
| Rosenberger et al. | Transferability of local density-assisted implicit solvation models for homogeneous fluid mixtures | |
| Michel et al. | Evaluation of host–guest binding thermodynamics of model cavities with grid cell theory | |
| Michielssens et al. | Molecular dynamics in principal component space | |
| Panosetti et al. | Global materials structure search with chemically motivated coordinates | |
| Lundberg et al. | Dispersion corrections to the surface tension at planar surfaces | |
| Sedeh et al. | Computing nonequilibrium conformational dynamics of structured nucleic acid assemblies | |
| Thomas et al. | Direct Observation of Salt Effects on Molecular Interactions through Explicit-Solvent Molecular Dynamics Simulations: Differential Effects on Electrostatic and Hydrophobic Interactions and Comparisons to Poisson− Boltzmann Theory | |
| Anandakrishnan et al. | An n log n generalized Born approximation | |
| Kolev et al. | Ab initio molecular dynamics of Na+ and Mg2+ countercations at the backbone of RNA in water solution | |
| Redhu et al. | Molecular modelling: a new scaffold for drug design | |
| Mori et al. | Integrated Hamiltonian Sampling: a simple and versatile method for free energy simulations and conformational sampling | |
| Okamoto et al. | Prediction of ligand binding affinity by the combination of replica-exchange method and double-decoupling method | |
| Harada et al. | Self-Avoiding Conformational Sampling Based on Histories of Past Conformational Searches | |
| WO2014115339A1 (ja) | 溶液中の生体高分子の構造揺らぎとダイナミクスを記述する方程式 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| 121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 13872785 Country of ref document: EP Kind code of ref document: A1 |
|
| NENP | Non-entry into the national phase |
Ref country code: DE |
|
| 122 | Ep: pct application non-entry in european phase |
Ref document number: 13872785 Country of ref document: EP Kind code of ref document: A1 |
|
| NENP | Non-entry into the national phase |
Ref country code: JP |