分子構造の予測とその応用
国立医薬品食品衛生研究所 有機化学部
栗原 正明
分子
の構造
核磁気共鳴(NMR)
赤外分光法(IR)
円二色性(CD)
X線構造解析
分子の構造と機能
化学シミュレーションにより 分子の構造を予測する 分光学的手法で 分子の構造を解析する分子力学法(MM)
分子動力学法(MD)
分子軌道法(MO)
分子を設計できる
創薬
分子の機能
化学計算
分子力学法(分子力場法)
Molecular Mechanics (MM)
古典力学、コンピュータ負荷小、大きな分子
分子動力学
Molecular Dynamics (MD)
古典力学、動的挙動
分子軌道法
Molecular Orbital (MO)
量子力学、精密な計算、コンピュータ負荷大
半経験的分子軌道法
ab initio 分子軌道法
分子力学法(分子力場法)
Molecular Mechanics (MM)
分子の立体配座の安定性や配座間のエネルギー差 を原子間に働く力によるポテンシャルエネルギー の総和によって計算する手法のことである。
分子力場
静電相互作用エネルギー
van der Waalsエネルギー ねじれ角エネルギー 変角エネルギー 結合伸縮エネルギー
E
total= Σ E
bond+ Σ E
bending+ Σ E
torsio+ Σ E
el+ Σ E
van( )
(
)
2 0 2 1 l l k l Ebond = −( )
(
)
2 02
1
θ
θ
θ
= k
−
E
bending( )
ω
= V(
+(
nω
−γ
)
)
Etorsion n 1 cos 2 ij j i el r q q E 0 4πε = ⎥ ⎥ ⎦ ⎤ ⎢ ⎢ ⎣ ⎡ ⎟ ⎠ ⎞ ⎜ ⎝ ⎛ − ⎟ ⎠ ⎞ ⎜ ⎝ ⎛ = 6 12 4 r r Evan ε σ σ分子力場
MM2, MM3(Allinger)
AMBER (Kollmann)
MMFF (Merck社)
CHARM (Karplus)
OPLS
構造最適化
V(n) = f (a
1, a
2--- a
3)
分子のポテンシャルエネルギー
構造最適化:Vが極小値をとるような原子核配置を求める
δ a
1δ V
=
δ a
2δ V
=
=
δ a
nδ V
=
0
座標 ポテン シ ヤ ル
コンフォメーショナルサーチ
Global minimum Local minimum Local minimum 最安定構造配座探索
米国コロンビア大学のProf. W. Clark Still
MacroModel
ver.1 in 1988
1946 born (Augusta, Georgia)
1964 graduated fromWinter Haven HS (Polk County, Florida) 1969 B.S., Ph.D. (1972) at Emory U. (Adv. David Goldsmith)
1973 Postdoc at Princeton U. (Computer-related, adv. not known) 1974-5 Postdoc at Columbia U. (Adv. Gilbert Stork)
1975-6 Professor at Vanderbilt U. (Nashville, TN) 1977-98 Professor at Columbia U. (NY, NY)
OH O O O O HO O O O O O OH N H O O Paclitaxel (Taxol)
Conformational Search: MCMM
Global minimum
MolFeat Control v2.2
Presentation file: C:\Documents and Settings\kurihara Caption: On Medium
Comment: On Small
Mouse Operation: Rotation/Scale Auto Camera Rotation: On
O O OH O O O HO O O O O N H O OH O
タンパク質中のリガンドの構造
リガンド リガンド 溶媒 タンパク質 OH O O O O HO O O O O O OH N H O O OH O O O O HO O O O O O OH N H O OO O OH O O O HO O O O O N H O OH O O O OH O O O HO O O O O N H O OH O OH 代謝酵素で水酸化される パクリタキセル ドセタキセル (シトクロムP450 2C8) O O OH O O O HO O O O O N H O O OH O
どうしてドセタキセルは代謝されないか?
水酸化されないCYP2C8
Schoch, G. A., et al., J.Biol.Chem. 279, 9497 (2004)
OH O O O O HO O O O O O OH N H O O O Docetaxel OH O O O O HO O O O O O OH N H O O Paclitaxel (Taxol) hydroxylation by 2C8 水酸化されない パクリタキセル ドセタキセル 水酸化 CYP2C8
docetaxel OH O O O O HO O O O O O OH N H O O Paclitaxel (Taxol) hydroxylation by 2C8 OH O O O O HO O O O O O OH N H O O O Docetaxel
H N N H H N O O N H H N N H H N N H O O O O O O O MeO OMe OMe MeO
MeO MeO MeO
OMe OMe OMe OMe OMe OMe
MeO MeO MeO OMe O
どんな構造化か?
310-helix
α-helix
planar
Left handed helix Right handed helix
Conformational Search : MM & Monte Calro
AMBER*, MMFF by MacroModel
Method of Conformational Search
H N N H H N O O N H H N N H H N N H O O O O O O O MeO OMe OMe MeO
MeO MeO MeO
OMe OMe OMe OMe OMe OMe
MeO MeO MeO OMe O 4: Octapeptide
H N N H H N O O N H H N N H H N N H O O O O O O O MeO OMe OMe MeO
MeO MeO MeO
OMe OMe OMe OMe OMe OMe
MeO MeO MeO OMe O
どんな構造化か?
左巻き α-へリックスH N N H H N N H O O COOEt O CF3CO (s) H N N H H N O O N H H N N H H N N H O O O O O O O MeO OMe OMe MeO
MeO MeO MeO
OMe OMe OMe OMe OMe OMe
MeO MeO MeO OMe O H N N H H N N H O O O H N COOEt O CF3CO
3
10-helix
planar
α-helix
2
α-Substituted 1α, 25-Dihydroxyvitamin D
3Analogues
1α, 25-dihydroxyvitamin D3 (100) 2α-substituted analogues
1: R= Me (400) OH OH HO OH OH HO R 1 3 25 2 A 活性型ビタミンD3の2αに置換基を入れるとどうして活性が強くなるのか?
OH OH HO R A H HO H R OH H H H H OH R H OH H H H α-form R=H: (54 : 46) β-form R=CH3: (83 : 17)
Conformation of A-Ring
H11 H12 Ser-278 Tyr-143 Arg-274 Ser-237 His-305 His-397 活性型ビタミンD3とV D R の結合様式
N. Rochel et al., Mol. Cell, 5, 173(2000)
OH OH β-form
ビ タ ミ ン Dレ セプ タ ー( VDR) のX線構造解析
OH OH HO A C D 活性型ビタミンD3 1α,25-dihydroxyvitaminD3 AX-ray Structure Modeled Structure OH OH HO OH OH HO OH 100 300 J. Org. Chem., 66, 8760-8771 (2001)
2
α-(Hydroxypropyl)-1α, 25-dihydroxyvitamin D
3X-ray Structure Modeled Structure OH OH HO OH OH HO O HIS 397 TYR 236 SER 237 TYR 143 ARG 274 ASP 144 HO
2
α-(ω-Hydroxyalkoxy)-1α, 25-dihydroxyvitamin D
3 Org. Lett., 2, 2619-2622 (2000) 100 180核内レセプター変異疾患に対する薬物の分子設計と合成
“Pharmacological rescue” of genetically impaired protein
ビタミンDレセプター(VDR)をターゲットとする創薬
VDRアゴニスト
VDRアンタゴニスト
変異VDRアゴニスト
骨粗鬆症薬
抗がん薬
Paget’s骨病薬
遺伝的くる病薬
(ビタミンD受容体機構異常症 : 特定疾患)Hormone responsive gene Hormone responsive gene Cell membrance Nucleus VDR 正常細胞 ビタミンD受容機構異常症の細胞 1α,25(OH)2ビタミンD3 1α,25(OH)2ビタミンD3 救済型リガンド
核内レセプター変異疾患に対する薬物の分子設計と合成
R274
h yd ro g en b o nd
Hereditary Vitamin D-Resistant Rickets : 遺伝的くる病
OH OH HO Leu
水素結合できず結合力低下
骨の成長に障害
OH OH HO OH OH OH HO O HO Ser-278 Tyr-143 Arg-274 Ser-237 His-305 His-397 OH OH β-form HOH2CH2CH2C Tyr-143 Ser-278 Arg-274 Ser-237 His-397 His-305 Asp-144 OH OH β-form HOH2CH2CH2OC
Structure-Based Design of Vitamin D3 Analogues Accessible to
VitaminD Receptor Mutant Related to Hereditary Vitamin D-Resistant Rickets
OH OH HO R274 h yd ro g en b o nd OH OH HO O HO
a nchor side cha in R274L
D144
anchor side chain at C2α on t he A-ring ne w hy dro g en b on ds 1200 1000 800 600 400 200 0 0 20 40 60 80 100
Functional binding activity of R274A VDR (protease assay) nM % maxi mal bi ndi n g (WT) 1α,25(OH) 2D3 O1C3 O2C3 正常な受容体と活性型ビタミンD3 変異疾患の受容体とアンカー型ビタミンD3アナログ
O O OH HO OH * * OH HO OH H H 1 3 25 1α,25-(OH)2D3 LG190178
X-ray structure (1α,25-(OH)2D3) (PDB ID: 1DB1) (2000)
Calculated structure (LG190178)
M. F. Boehm, et. al., Chem. Biol., 6, 265-275 (1999).
Modeled Structures of Enantiomers of LG190178 in the VDR Ligand Binding Domain SR RR O O OH HO OH +3.56 kcal/mol * * O O OH HO OH O O OH HO OH R S R R His-305 His-397 Ser-237 Arg-274 Arg-274 Ser-237 His-397 His-305 SS +1.87 kcal/mol O O OH HO OH s s Arg-274 His-397 His-305 Ser-237 LG190178 RS +4.94 kcal/mol O O OH HO OH R S His-305 His-397 Ser-237 Arg-274 ΔE=0 kcal/mol
Transcriptional activity of isomers of LG190178 against VDR, VDRR274L Caco-2 / 5% FCS 3) HOS / 5% FCS 2) HOS/SF 1) 0.112 0.482 0.831 28.3 100 VDR affinity 0.69 131.4 83.2 15.6 YR304 (R,S) 0.17 22.1 16.8 4.68 YR303 (R,R) 1.0 21.1 9.9 1.66 YR302 (S,S) 0.22 1.80 0.79 0.0396 YR301 (S,R) 4.4 3.52 0.555 0.0106 1,25(OH)2D3 Transcription EC50 (nM) Mutant VDR Wild type VDR Compound YR303 (R, R) YR301(S, R) O O OH HO OH (S) (S) O O OH HO OH (S) (R) YR302(S, S) O O OH HO OH (R) (R) O O OH HO OH (R) (S) YR304 (R, S)
X-ray structure (PDB ID: 2ZFX) Calculated structure YR301 O O OH HO OH (S) (R)
OH HO OH H H 1 3 25 1α,25-(OH)2D3 YR301(S, R) O O OH HO OH (S) (R) TYR-143
無承認無許可医薬品を迅速に規制するには?
N HN N N O O S N N O OSildenafil (Viagra)
N HN N N O N O N HO O N HN N N O S N O O O N HN N N O N O N O N O N H H Cl O O O O N HN N N O N O N O N HN N N O O2N O N HN N N S S N O O O N N HN N N S S N O N O O N N N H O O H O ON HN N N O O S N N O O X線の構造はNature, 4 2 5 , 9 8 -1 0 2 ( 2 0 0 3 ) よ り
Structure of the catalytic domain of human phosphodiesterase 5 with bound drug molecules
Sildenafil N HN N N O O N N O Hongdenafil
N N N N O O S N N O O GLN817 GLN817と水素結 合を形成してい る NH O N H H O H H GLN817 Sildenafil Sildenafil の結合様式
N HN N N O O S N N O O N N HN N O O S N N O O N N O O H O O N H X線の構造はNature, 4 2 5 , 9 8 -1 0 2 ( 2 0 0 3 ) よ り
Structure of the catalytic domain of human phosphodiesterase 5 with bound drug molecules
N HN N N O S N O N O O N N N H O O H O O
Chem. Pharm. Bull., 57, 185-189(2009)
シクロペンチルナフィル