構造活性相関部会・ニュースレター <1 April, 2024>

SAR News No.46「目次」///// Perspective/Retrospective /////
ヒストンのアセチル化を介したクロマチンの制御研究
･･･

1

･･･

9

･･･

16

広川 貴次

･･･

29

加藤 幸一郎

･･･

31

･･･

32

･･･

34

理化学研究所梅原 崇史///// Cutting Edge /////NSD2 によるヌクレオソーム上ヒストン H3K36 メチル化の制御と破綻
横浜市立大学仙石 徹///// SAR Presentation Award /////
受賞コメントおよび講演要旨///// Activities /////＜報告＞
第 51 回構造活性相関シンポジウム 開催報告
筑波大学
＜会告＞
構造活性フォーラム 2024 会告九州大学
第 52 回構造活性相関シンポジウム会告
国立医薬品食品衛生研究所

編集後記

古濱 彩子

SAR News No.46 (Apr. 2024)///// Perspective/Retrospective /////

ヒストンのアセチル化を介したクロマチンの制御研究
理化学研究所 生命機能科学研究センター

梅原 崇史

1. はじめに真核細胞の核内に存在するエピゲノムでは、ゲノム DNA を収納するヒストンタンパク質に様々な種類の翻訳後修飾が見られる。この翻訳後修飾の付加、消去とその認識は、細胞の恒常性の維持や特異性の発揮に関わり、その制御機構に異常が生じると様々な疾患につながることが知られている。ヒストンに対する様々な翻訳後修飾の中で、リジン残基のアセチル化は遺伝子転写の活性化とよく相関する。そのため、近年のケミカルバイオロジーではヒストンのアセチル化制御に関わる化合物が数多く開発され、それらの化合物を用いてエピゲノムの制御分子機構の理解やその人工的な調節を指向した研究が進んできている。この総説では、ヒストンのアセチル化を介したクロマチンの制御研究について、前半ではその分子機構に関わるタンパク質とそれらのタンパク質に対する代表的な制御分子を紹介する。後半では、そのようなヒストンアセチル化の制御分子を活用してクロマチン制御の分子機構解析に応用した研究を紹介する。

2. ヒストンのアセチル化の制御分子基盤
2.1 概論真核生物のエピゲノムはその多くの部分において、ヒストンタンパク質と DNA からなるヌク
レオソームと呼ばれる構造を凝縮単位とするクロマチン構造を形成している。ヌクレオソームの形成とその制御に関わるヒストンタンパク質は、コアヒストン（H2A, H2B, H3, H4）とリンカ
ーヒストン（H1 など）に分けられる。エピゲノム制御に関わる化学修飾の主要標的はコアヒストンである。コアヒストンは特定の古細菌においては 1 種類のタンパク質のみが存在し、酵母やヒトなどの真核生物においてはそのタンパク質を進化的な起源とすると考えられる H2A, H2B,
H3, H4 の 4 種類のヒストンが存在する。
真核生物のコアヒストンの C 末端側はいずれも折り畳まれたコア領域を形成し、H2A は H2Bと、H3 は H4 とそれぞれ握手しているようなヘテロ二量体を形成する。通常、H3 と H4 のヘテロ二量体はさらに二量体を形成して[H3-H4]2 のヘテロ四量体を形成し､その左右に H2A-H2B のヘテロ二量体が一つずつ結合してヒストン八量体が形成される。このヒストン八量体の外側に145〜147 塩基対の DNA が巻きつくことでヌクレオソームコア粒子と呼ばれる凝縮構造が形成される。ヌクレオソームコア粒子の DNA の両側にはヌクレオソーム構造を形成しない 20〜30
塩基対のリンカーDNA が通常存在し、この位置にリンカーヒストンが結合するとヌクレオソー
ムから DNA がほどけにくいさらに安定した凝縮構造（クロマトソーム）が形成される。
4 種類のコアヒストンの N 末端側はいずれも特定の構造を取らないテイル領域を持ち、ヌク
レオソームを巻いている DNA の外側に出ることができる。そのため、ヒストンの N 末端テイルは外部から様々な翻訳後修飾（アセチル化やメチル化など）を受けることで情報をやり取りでき、細胞核内においてエピゲノムの位置ごとに可逆的な情報を発信できる拠点となっている（図 1）。
このヒストンの N 末端テイルに存在する複数のリジン残基がアセチル化修飾を受けた場合、通常、そのヌクレオソーム上または近傍に位置する DNA は RNA に活発に転写される。このヒストンのリジン残基のアセチル化と RNA の転写活性の相関は 1964 年に報告された [1]。以来 60年にわたってヒストンのアセチル化が生物の機能制御にもたらす研究が続けられてきている。
現在、ヒストンのアセチル化は DNA 配列の変化に依存せずに特定の遺伝子発現の状態を確立す
るエピジェネティックな細胞記憶の一端を担うと考えられている。
ヒストンのリジン残基のアセチル化を制御するタンパク質の作用は主に、(1)アセチル化の付加（ヒストンアセチル化酵素）、(2)アセチル化の消去（ヒストン脱アセチル化酵素）、(3)アセチ
ル化の認識（アセチル化ヒストン結合因子）の 3 種類が知られている。これらのタンパク質はい-1-

SAR News No.46 (Apr. 2024)

ずれもエピゲノムからの遺伝子転写の制御に深く関わっている。以下の節では、この 3 種類のカ
テゴリーごとに代表的なタンパク質とその制御分子を紹介する。

図 1. ヌクレオソームを構成するヒストンと代表的な翻訳後修飾
2.2 ヒストンアセチル化酵素と代表的な制御分子
ヒストンのアセチル化は、細胞核内のクロマチンのヒストンに対して起こる反応と、細胞質においてクロマチンに取り込まれる前の新生ヒストンに対して起こる反応がある。真核細胞においてそれぞれのアセチル化反応は、前者を担う A 型と後者を担う B 型のヒストンアセチル化酵素によって触媒される。ヒストンのアセチル化を介したエピジェネティックな細胞記憶は主として A 型のヒストンアセチル化酵素によって担われると考えられている。A 型の代表的なヒス
トンアセチル化酵素としては p300/CBP や GCN5 関連ファミリー、MYST ファミリーなどが挙げられる。また、B 型のヒストンアセチル化酵素としては Hat1 が挙げられる。
4 種類のコアヒストンのどの N 末端テイルをアセチル化するかは酵素ごとに異なる。細胞の恒常性を制御する p300/CBP は多数の転写因子と相互作用して特定の遺伝子転写に関わるが、遺伝子の急速な転写活性化には p300/CBP のアセチル化酵素活性が必要である。p300/CBP は 4 種類のコアヒストンすべての N 末端テイルをアセチル化できる点が他のアセチル化酵素に見られない特徴である。一般的に GCN5 関連ファミリーは H3 に対するアセチル化活性が高く、MYST
ファミリーは H4 に対するアセチル化活性が高い。また、Hat1 は細胞質において主に新生ヒスト
ン H4 の N 末端から 5 番目と 12 番目のリジン残基をアセチル化する。このクロマチンに取り込まれる前のヒストンのアセチル化もエピジェネティクスの制御に貢献していると考えられる。
これらのヒストンアセチル化酵素のいくつかは、分子内にヒストンのアセチル化を認識するドメインを持っている。例えば p300/CBP はヒストンアセチル化の酵素活性ドメイン以外にヒストンのアセチル化状態を認識するブロモドメインを持つ。そのため、p300/CBP はクロマチンにおけるヒストンのアセチル化領域を認識し、アセチル化を「読み書き」する活性を介してクロマチンの特定のヒストンアセチル化状態の維持やその娘細胞への継承に働くと考えられている。
p300/CBP のアセチル化酵素ドメインに対する阻害剤としては、急性骨髄性白血病細胞などの増殖を抑制する A-485 が挙げられる [2]
（図 2）
。また、
ヒストンのアセチル化を認識する p300/CBPのブロモドメインに対しては、CBP30（図 2）などの阻害剤が開発されている [3]。

図 2. ヒストンアセチル化酵素の制御分子の例-2-

SAR News No.46 (Apr. 2024)

2.3 ヒストン脱アセチル化酵素と代表的な制御分子
細胞内でアセチル化されたヒストンに対しては、能動的に脱修飾（脱アセチル化）する酵素群が知られている。これらのヒストン脱アセチル化酵素の存在は、ヒストンアセチル化酵素と同様に、真核生物においてよく保存されている。ヒトでは 18 種類のヒストン脱アセチル化酵素が知られており、配列の相同性などに基づいて 4 種類のクラスに分類されている。これらの酵素は、酵素活性中心に亜鉛イオンを含むクラス（HDAC ファミリー, クラス I, II, IV）と、ニコチンア
ミドアデニンジヌクレオチドを補酵素とするクラス（SIRT ファミリー, クラス III）で構成されている [4]。このうち、HDAC ファミリーのヒストン脱アセチル化酵素に対しては、がん細胞の増殖制御などに関わるトリコスタチン A などの天然物が低濃度で作用する阻害剤であることが比較的初期から知られている。トリコスタチン A の類縁体の SAHA（図 3）は皮膚 T 細胞性リンパ腫に対する治療薬として開発され、ヒストンの翻訳後修飾を制御標的とする医薬品開発の端緒と位置づけられる。HDAC ファミリーのヒストン脱アセチル化酵素に対しては、その他にも
皮膚 T 細胞性リンパ腫や末梢性 T 細胞性リンパ腫の治療薬として用いられる FK228（図 3）などの制御分子が開発されており、エピゲノム創薬をリードする重要な研究が数多く存在する。
また、SIRT ファミリーのヒストン脱アセチル化酵素に対しても、SIRT1 に対する阻害剤の EX527 や活性化剤の resveratrol などが知られている。

図 3. ヒストン脱アセチル化酵素の制御分子の例
2.4 アセチル化ヒストン結合因子と代表的な制御分子
アセチル化ヒストンを認識するドメインも真核生物において酵母からヒトまでよく保存されている。アセチル化ヒストン結合ドメインとしては、ブロモドメインと YEATS ドメインが代表例として挙げられる。ヒトではブロモドメインを持つタンパク質は 46 種類あり、YEATS ドメイ
ンを持つタンパク質が 4 種類存在する。これらのドメインを持つタンパク質やその複合体によるアセチル化ヒストンへの結合は通常、クロマチン上の特定の遺伝子の転写を活性化し、その結合の異常が疾患の発症と相関するケースが知られている。
ブロモドメインは複数のサブファミリーに分類され、その一つは p300/CBP や GCN5 関連ファミリーなどのヒストンアセチル化酵素が持つブロモドメインである。他に注目されているブロ
モドメインとしては、BET（Bromodomain and Extra-Terminal domain）ファミリーのタンパク質に含まれる BET ブロモドメインが挙げられる。ヒトでは BRD2, BRD3, BRD4, BRDT と呼ばれる 4
種類のタンパク質が BET ファミリーに属している。BET ブロモドメインに対しては JQ1（図 4）が代表的な阻害剤として挙げられる [5]。BET ファミリーのうち、BRD4 は遺伝子転写が亢進しているクロマチン領域に局所的に集合する。このようなクロマチン領域はしばしば、エンハンサーと呼ばれる遺伝子発現を強めるために転写因子が結合する DNA 配列を複数持つ 10-kb 以上の長い DNA 領域で構成される。これらはスーパーエンハンサーと呼ばれ、通常の 1-kb 以下の短い
DNA 配列のエンハンサーと区別されるようになってきている [6]。BRD4 が関わるスーパーエン
ハンサーの形成や機能発揮は、細胞種に依存しない遺伝子転写の活性化機構であり、この仕組み
によってがん細胞の増殖機構の一端が包括的に説明できるという考え方がある。そのため、多種類のがんに対して、BET ブロモドメイン阻害剤を用いた臨床試験が現在進められている。

-3-

SAR News No.46 (Apr. 2024)

図 4. アセチル化ヒストン結合因子の制御分子の例
BET ブロモドメインの阻害剤は、クロマチンにおけるスーパーエンハンサーの位置を同定したり、その機能を解析したりする上で重要な役割を果たしている。また、ケミカルバイオロジー
で特定の機能を持つ制御分子を開発・実証するモデル系としても利用されている。例えば dBET1（図 4）は、JQ1 による BET ブロモドメインとの結合能を利用して BET ファミリータンパク質を細胞内で選択的に分解する PROTACs の化学ツールである [7]。次節では、BET ファミリーや
p300/CBP のブロモドメインに対する阻害剤を活用した最近のクロマチン制御研究を紹介する。

3. ブロモドメイン阻害剤を活用したクロマチンの制御研究
3.1 化合物によるアセチル化ヒストン結合因子の選択的分解
細胞内に存在する個々のタンパク質の機能を理解するために、特定のタンパク質を細胞から任意のタイミングで分解して表現型などの変化を調べる方法がある。この方法には遺伝子レベ
ルまたはタンパク質レベルでアプローチする方法があり、どちらの方法でも特定の化合物の投与に応じて特定の遺伝子産物やタンパク質を選択的に分解させる方法が開発されている。前者の方法としては主にオーキシンデグロン系が挙げられ、後者では主に PROTACs や SNIPERs と呼ばれる方法が挙げられる [8]。PROTACs や SNIPERs では、標的タンパク質に結合する低分子化合物をリンカーを介して E3 ユビキチンリガーゼのリガンドと連結した化合物を用いる。この化合物を細胞や組織に投与すると、タンパク質レベルで標的タンパク質をプロテアソームに誘導し、短時間で選択的に分解することができる。この方法は、細胞内で特定のタンパク質の機能解析に有用なだけでなく、疾患治療への応用も期待されている。
先に述べたように、BET ファミリータンパク質はアセチル化ヒストンへの結合を介してスーパーエンハンサーで機能することにより、遺伝子転写の活性化やがん細胞の増殖に重要である。
BET ブロモドメイン阻害剤による BET ファミリータンパク質の機能阻害はエピゲノム創薬の重要な制御標的と位置づけられているため、BET ブロモドメイン阻害剤の類縁体を用いて BET ファミリータンパク質を選択的に分解できれば、がん細胞などに対するさらなる制御効果が期待できる。そこで PROTACs や SNIPERs を BET ブロモドメイン阻害剤に適用することで BET ファ
ミリータンパク質を選択的に分解させる研究が進められてきている（図 5）[7,9,10]。

図 5. 化合物による BET ファミリータンパク質の選択的な分解-4-

SAR News No.46 (Apr. 2024)

dBET1（図 4）は、BET ファミリータンパク質を選択的に分解するために開発された化合物である。この化合物は、BET ブロモドメイン阻害剤の JQ1 の構造骨格にサリドマイド系の化合物
がリンカーを介して繋げられている。サリドマイド系化合物は、Cullin・RING ユビキチンリガーゼ複合体のサブユニットであるセレブロン E3 リガーゼのリガンドとして働く。そのため、
dBET1 の JQ1 部分に結合した BET ファミリータンパク質は、dBET1 のサリドマイド部分がセレ
ブロンに結合することでプロテアソームに運ばれて選択的に分解される（図 5）。dBET1 は細胞膜を透過できるため、細胞に投与すると体細胞性の BET ファミリータンパク質である BRD2、
BRD3 および BRD4 の量が選択的に減少する。BRD4 は急性骨髄性白血病などのがん細胞の増殖制御に関わることが知られているが、dBET1 を急性骨髄性白血病細胞のマウスモデルに投与す
ると急性骨髄性白血病の病変が抑えられることがわかった [7]。
PROTACs は現在、アンドロゲン受容体やエストロゲン受容体の選択的分解を制御標的としてそれぞれ前立腺がん患者や乳がん患者の治療を目的とする臨床研究が進められている [11]。ブロモドメイン含有タンパク質に対しては、現在、BRD9 の PROTACs 戦略によって滑膜肉腫の治療を目的とする臨床研究が進められている。
3.2 化合物によるアセチル化ヒストン結合因子のクロマチン結合部位の同定
PROTACs や SNIPERs では、標的タンパク質に結合する分子（JQ1 など）と E3 リガーゼに結合する分子（サリドマイド系化合物など）をリンカーで連結し、細胞内で標的タンパク質をタン
パク質分解系に誘導して選択的に分解させる点に特徴がある。このアイデアを抽象化すると「機能が異なる化合物を連結することで細胞内の標的タンパク質に何らかの機能を付与できる」こ
とが考えられる。付与できる機能は連結する化合物の性質に依存するが、工夫によっては標的タ
ンパク質に分解以外の制御的な機能を付与したり、細胞内での標的タンパク質の局在を検出したりできる。ここでは細胞内での標的タンパク質の検出について、クロマチンに対する結合部位
や組織内での分布を解析した研究を紹介する（図 6）。
クロマチンに対する結合部位の同定方法としては、Chem-seq 法 [12]や Chem-map 法 [13]が挙げられる。これらの方法では、BET ブロモドメイン阻害剤の JQ1 をリンカーを介してビオチンと連結した化合物を用いる。Chem-seq 法では、細胞やその抽出液のクロマチン DNA とタンパク質を架橋剤で架橋し、ビオチン化 JQ1 を添加してクロマチン上の BET タンパク質に化合物を結合させ、クロマチンを剪断する。その後、ビオチンとアビジンの親和性を利用して BET タンパ
ク質を含むクロマチン断片を取得し、次世代シーケンサーによる DNA 配列の網羅的解析によってビオチン化 JQ1 がクロマチンに結合する部位をマッピングできる [12]。JQ1 は試験管内では体細胞性の BET ファミリータンパク質である BRD2、BRD3 および BRD4 のブロモドメインとほぼ同程度の親和性で結合するため、ビオチン化 JQ1 のクロマチン結合部位は原理的にはこれらの BET ファミリータンパク質のクロマチン結合部位を総合的に反映すると考えられる。しかし多発性骨髄腫細胞での検討では、ビオチン化 JQ1 のクロマチン局在は BRD4 の局在と最もよく相関し、次いで BRD2 と相関が高く、BRD3 との相関は最も低かった。この結果は、解析に用いた細胞種での各タンパク質の機能的差異を反映している可能性が考えられる。

図 6. 化合物による BET ファミリータンパク質のクロマチン結合部位の同定-5-

SAR News No.46 (Apr. 2024)

一方、Chem-seq 法の改良版である Chem-map 法では、クロマチンを剪断する代わりに Tn5 ト
ランスポゼースを用いた DNA のタグメンテーションを行う。ビオチン化 JQ1 と結合した BET
ファミリータンパク質を含む細胞核内のクロマチンを抗ビオチン抗体で検出し、その抗体に対する二次抗体に Tn5 を融合しておくことで、BET タンパク質を含むクロマチンの局所的な DNA
断片をタグ付加して取得できる。この方法により、Chem-map 法では Chem-seq 法よりもシグナ
ル品質に優れたビオチン化 JQ1 のクロマチン結合部位の in situ マッピングを実現している [13]。
ここまでに紹介した PROTACs や Chem-seq 法、Chem-map 法では 2 種類の化合物をリンカーで連結した化合物を用いる。この連結によって化合物の分子量が大きくなると、細胞内への膜透過性や細胞内における化合物・標的タンパク質間の相互作用を損なう可能性があり、標的タンパ
ク質による作用機構の解明が難しくなる場合がある。そこで、BET 阻害剤に、クリック反応による直交ライゲーションが可能な官能基を導入した化合物が開発された [14]。このクリック可能な化合物を用いることにより、ビオチンなどのアフィニティータグだけでなく、蛍光プローブなどを細胞内や個体内で連結することができ、標的タンパク質を制御・検出するための方法が多様化された。
この方法を用いて、急性骨髄性白血病のマウスモデルでクリック型の JQ1 プローブの細胞内での蓄積を検討した結果、急性骨髄性白血病細胞の 2 種類の組織源である脾臓と骨髄では化合物の蓄積が異なり、脾臓に含まれる急性骨髄性白血病細胞では骨髄よりも化合物濃度が有意に高かった [14]。これらの結果は、BET 阻害剤は末梢血液を循環している急性骨髄性白血病細胞や脾臓内に局在する急性骨髄性白血病細胞を速やかに除去できるが、骨髄内の急性骨髄性白血
病細胞を除去できないという報告 [15]と整合性がある。
3.3 ヒストンのアセチル化が伝播する分子機構の解析
真核細胞においてヒストンのアセチル化は遺伝子転写の活性化に重要であるが、その直接的な読み書きの分子機構についてはこれまでよく分かっていなかった。最近、極低温電子顕微鏡構
造解析によって、p300/CBP がヌクレオソーム内でどのようにヒストンのアセチル化を読み書きするかを示唆する研究が報告された [16]。この構造では、p300/CBP はアセチル化の読み取りを担うブロモドメインのポケット内部でヒストン H4 の N 末端テイルのアセチル化を認識し、ヌクレオソーム上で多方向に回転することで同じヌクレオソームのヒストン H4 以外のヒストンの N 末端テイルをアセチル化できることが示唆された（図 7 左）
。また、p300/CBP のブロモド
メインに対する阻害剤の CBP30（図 2）を用いて、p300 がヌクレオソーム内のヒストンテイルのアセチル化をどのような方向性で伝播するかが調べられた。

図 7. 構造解析と化合物の利用によるヒストンアセチル化の伝播機構の解析
p300/CBP は 4 種類のコアヒストンをいずれもアセチル化できることから、ヌクレオソームのヒストンテイルを無差別にアセチル化すると従来考えられてきている。p300/CBP がヌクレオソーム内のヒストンアセチル化を伝播する際に特定の方向性があるかどうかは、あらかじめヌク-6-SAR News No.46 (Apr. 2024)

レオソームのヒストンの特定のリジン残基をアセチル化した場合に p300 がどのヒストンテイル
をアセチル化するか、また CBP30 の有無によってそのアセチル化の度合いがどの程度影響を受けるかを調べればわかる。この解析の結果、p300 はヒストン H3-H4 四量体内ではアセチル化を双方向に伝播する一方、ヒストン H2B テイルに対してはヒストン H3 テイルまたはヒストン H4
テイルにあらかじめアセチル化がある場合に依存して一方向的にアセチル化を伝播した（図 7右）。この結果から、p300 がヌクレオソーム内のヒストンをアセチル化する方向性には「双方向性」と「一方向性」の 2 種類があることがわかった。
文献[16]の構造解析では、p300/CBP はヒストン H4 テイルのアセチル化を認識し、ヌクレオソーム上で回転して多様な結合様式でヌクレオソームと結合することで、認識した H4 テイルの近傍に位置する H4 以外の複数のヒストンテイルをアセチル化しうることが示された。なかでも、
p300/CBP はヒストン H4 テイルのアセチル化を認識した場合にヒストン H2B の N 末端テイルを主要なターゲットとしてアセチル化することがわかった。さらに H2B の N 末端テイルのアセチ
ル化は、ヌクレオソームからの H2A-H2B 二量体の解離を選択的に促進することが示された。これらの結果を総合して、p300/CBP が H3-H4 四量体内でヒストン N 末端テイルのアセチル化を「複製」してエピジェネティックな情報（ヒストンのアセチル化）を継承し、それを H3-H4 四量体から H2B-H2A 二量体に「転写」して局所的なヌクレオソームの不安定化を介して特定の遺伝子転写を活性化するモデル（エピセントラル・モデル）が提唱された（図 8）。
このモデルは、H3-H4 四量体と H2B-H2A 二量体のアセチル化が異なる役割を果たす点に特徴がある。すなわち、ヒストンのアセチル化の伝播に関して、H3-H4 四量体の役割は DNA の役
割と類似しており、p300/CBP が H3-H4 四量体内で「複製」したアセチル化情報を細胞周期を越えて自己永続させるための情報媒体と位置づけられる。また、H2B-H2A 二量体の役割は自己永続せずに用いられる RNA の役割と類似しており、p300/CBP が H3-H4 四量体から H2B-H2A 二量体にアセチル化を「転写」することでヌクレオソームから H2B-H2A 二量体の脱離が促され、その位置の DNA を RNA ポリメラーゼが生産的に転写できるようになる。これらの結果をまとめると、真核生物の遺伝子転写制御の本質の一端は、クロマチンの個々のヌクレオソームが「エピジェネティック情報を継承する H3-H4 四量体」と「エピジェネティック情報を発現する H2BH2A 二量体」を併せ持つ点にあると考えられた。

図 8. ヒストンアセチル化の情報伝達を介した遺伝子転写制御モデル

4. 終わりに
本節では、ヒストンアセチル化酵素、ヒストン脱アセチル化酵素、アセチル化ヒストン結合因子とそれらのタンパク質に対する代表的な制御分子を紹介し、アセチル化ヒストン結合因子に対する阻害剤やその類縁化合物を用いた最近のクロマチン制御研究を紹介した。エピゲノムの制御は今回述べたヒストンのアセチル化を介した制御だけでなく、ヒストンのメチル化、リン酸
化、ユビキチン化などの制御やゲノム DNA のメチル化の制御など、多岐にわたる経路によって担われている。そのため、2 種類の異なる化学修飾の制御経路を同時に阻害する制御分子も開発されている。また、近年では個体の世代を超えて継承されるエピジェネティック情報にも注目が集まっている。今後は、食事や化合物、ストレスなどの環境条件に応じたエピゲノムの変化とその制御機構の解析にも注目が集まることが期待される。

-7-

SAR News No.46 (Apr. 2024)

参考文献[1] Allfrey, V.G., Faulkner, R., Mirsky, A.E., Acetylation and methylation of histones and their possible
role in the regulation of RNA synthesis, Proc. Natl. Acad. Sci. USA, 51, 786-794 (1964).
[2] Lasko, L.M., Jakob, C.G., Edalji, R.P., et al. Discovery of a selective catalytic p300/CBP inhibitor thattargets lineage-specific tumours, Nature, 550, 128-132 (2017).
[3] Hammitzsch, A., Tallant, C., Fedorov, O., et al. CBP30, a selective CBP/p300 bromodomain inhibitor,
suppresses human Th17 responses, Proc. Natl. Acad. Sci. USA, 112, 10768-10773 (2015).
[4] Seto, E., Yoshida, M., Erasers of histone acetylation: the histone deacetylase enzymes, Cold Spring
Harb. Perspect. Biol., 6, a018713 (2014).
[5] Filippakopoulos, P., Qi, J., Picaud, S., et al. Selective inhibition of BET bromodomains, Nature, 468,
1067-1073 (2010).
[6] Lovén, J., Hoke, H.A., Lin, C.Y., et al. Selective inhibition of tumor oncogenes by disruption of superenhancers, Cell, 153, 320-334 (2013).
[7] Winter, G.E. Buckley, D.L., Paulk, J., et al. DRUG DEVELOPMENT. Phthalimide conjugation as a
strategy for in vivo target protein degradation, Science, 348, 1376-1381 (2015).
[8] Naito, M., Ohoka, N., Shibata, N., et al. Targeted Protein Degradation by Chimeric Small Molecules,
PROTACs and SNIPERs, Front. Chem., 7, 849 (2019).
[9] Zengerle, M., Chan, K.-H., Ciulli, A., Selective small molecule induced degradation of the BET
bromodomain protein BRD4, ACS Chem. Biol., 10, 1770-1777 (2015).
[10] Ohoka, N., Okuhira, K., Ito, M., et al. In vivo knockdown of pathogenic proteins via specific andnongenetic inhibitor of apoptosis protein (IAP)-dependent protein erasers (SNIPERs), J. Biol. Chem.,
292, 4556-4570 (2017).
[11] Békés, M., Langley, D.R., Crews, C.M. PROTAC targeted protein degraders: the past is prologue, Nat.
Rev. Drug. Discov., 21, 181-200 (2022).
[12] Anders, L., Guenther, M.G., Qi, J., et al. Genome-wide localization of small molecules, Nat.
Biotechnol., 32, 92-96 (2014).
[13] Yu, Z., Spiegel, J., Melidis, L., et al. Chem-map profiles drug binding to chromatin in cells, Nat.
Biotechnol., 41, 1265-1271 (2023).
[14] Tyler, D.S., Vappiani, J., Cañeque, T., et al. Click chemistry enables preclinical evaluation of targeted
epigenetic therapies, Science, 356, 1397-1401 (2017).
[15] Dawson, M.A., Prinjha, R.K., Dittmann, A., et al. Inhibition of BET recruitment to chromatin as aneffective treatment for MLL-fusion leukaemia, Nature, 478, 529-533 (2011).
[16] Kikuchi, M., Morita, S., Wakamori, M., et al. Epigenetic mechanisms to propagate histone acetylation
by p300/CBP, Nat. Commun., 14, 4103 (2023).

-8-

SAR News No.46 (Apr. 2024)///// Cutting Edge /////NSD2 によるヌクレオソーム上ヒストン H3K36 メチル化の制御と破綻
横浜市立大学 医学部

仙石 徹

1. はじめにヒストンは様々な翻訳後修飾を受け、それにより転写・複製・修復などのゲノム機能が制御されている。一般に、ヒストンの修飾状態はそれを導入する酵素である「writer」と除去する酵素である「eraser」の拮抗する作用のせめぎあいによって確立され、それを特異的に認識するタンパク質である「reader」によって下流の生物学的応答が引き起こされる（図 1）[1]。ヒストン修飾の破綻はがんなどの様々な疾患の原因となっており、それを正常化する薬剤として writer・eraser・reader を標的とした阻害剤の開発が盛んに行われている（表 1）[2]。

図 1. ヒストン修飾に関係するタンパク質の機能
表 1. ヒストン修飾に関係する創薬標的の例（なお、しばしば単一のタンパク質が reader
ドメインと writer/eraser ドメインの両方を持つ）
タンパク質やドメイン
修飾
機能
Eraser
リジンアセチル化
HDAC
Reader
リジンアセチル化
ブロモドメインWriter
リジンアセチル化
p300
Eraser
H3K4 メチル化
LSD1Writer
H3K9 メチル化
G9aWriter
H3K27 メチル化
EZH2Writer
H3K36 メチル化
NSD2Writer
H3K79 メチル化
Dot1L
ヒストンのリジン残基は主に H3 の K4・K9・K27・K36・K79 と H4 の K20 でメチル化を受け、その残基位置とメチル化の程度（モノ・ジ・トリメチル化）によって異なる生物学的応答が引き起こされる [3]。例えば、ヒストン H3K36 のトリメチル化（H3K36me3）は活発に転写される領域に、ジメチル化（H3K36me2）は遺伝子間領域に多く存在し、どちらも転写活性化に関与する。一方で、H3K27 のトリまたはジメチル化（H3K27me3/2）は条件的ヘテロクロマチン形成を伴う転写不活性化に働き、細胞種特異的な発現パターンを示す遺伝子の異所的な発現を抑制
している。

-9-

SAR News No.46 (Apr. 2024)

本研究で対象とした NSD2（別名 MMSET・WHSC1）はヒストン H3K36 ジメチル化活性を持
つ NSD ファミリー（NSD1・2・3）[4]の一員であり、NSD ファミリーによる H3K36me2 レベルの厳密な制御が正常な発生に必要である。例えば、NSD2 のハプロ不全は成長障害や精神遅滞を伴う Wolf-Hirschhorn 症候群を引き起こす [5]。
一方で、NSD2 活性の亢進による細胞内 H3K36me2
レベルの異常な上昇は様々ながん（主に血液がん）を引き起こす。多発性骨髄腫の約 20%で染色体転座による NSD2 の過剰発現が見られ [6]、また一部のリンパ性白血病において NSD2 のメチル化活性を異常亢進する点変異である E1099K や T1150A が見られる [7,8]。さらに、リンカー
ヒストン H1 は NSD2 によるヌクレオソーム上 H3K36 メチル化を阻害することが知られており、一部のリンパ腫では H1 遺伝子の機能喪失変異により H3K36me2 レベルの異常な上昇が観察される [9,10]。これらの事実から NSD2 阻害剤は有望な抗がん剤候補化合物と考えられ [11]、その特 異 的 阻 害 剤 KTX-1001 の 多 発 性 骨 髄 腫 に 対 す る 治 験 が 2023 年 に 開 始 さ れ て い る（NCT05651932）。
ヌクレオソーム中ではヒストン 8 量体の球状ドメインが密接な相互作用によりコア構造を形成し、DNA はその外周に巻き付いている。ヒストンの N 末端領域に存在する「ヒストンテール」
は特定の立体構造を形成せずコア構造から外側に突き出ており、H3K36 はコア構造とヒストン
テールの境界である DNA 近傍に存在している（図 2 左）。そのため、H3K36 の周囲は立体的に込み入っており、そのままでは酵素がアクセスすることが難しい状態にある。NSD2 がヌクレオソームをどのように認識して込み入った位置の H3K36 をメチル化するかは未解明であった。
我々は、NSD2 がヌクレオソーム上の H3K36 にアクセスする機構、リンパ球の正常な増殖に必要な H1 による NSD2 活性の抑制機構、そして発がん性変異 E1099K や T1150A によるメチル化活性の異常亢進機構について知見を得るために、NSD2 とヌクレオソームとの複合体の立体構
造解析を行った [12]。

図 2. 左: 標準的なヌクレオソームの構造（PDB 1KX5）。中: NSD2-ヌクレオソーム複合体の構造（PDB 7E8D）。右: H1-ヌクレオソーム複合体の構造（PDB 5NL0）。

2. 立体構造
クライオ電子顕微鏡を用いた単粒子解析を行い、NSD2 の触媒ドメインがヌクレオソームに結合した構造を決定した（全体分解能約 2.8 Å）。
2.1 構造決定
構造解析に用いたサンプルの調製法について記す。予備的実験により長いリンカーDNA を含
むヌクレオソームが NSD2 の良い基質になることが分かったので、両端に約 20 bp のリンカー領域を持つ DNA を用いた。NSD2 には発がん性変異 E1099K を導入し、また複合体の安定化のた
めに NSD2 の触媒ドメインと H4 との融合タンパク質を発現させ他のコアヒストンと DNA と混合してヌクレオソームを再構成した。さらに、標的リジン残基がメチオニンに置換された H3 変

- 10 -

SAR News No.46 (Apr. 2024)

異体（H3K36M）を用いて、H3 と NSD2 との相互作用のさらなる安定化を行った [13]。メチル基供与体 SAM のアナログとして sinefungin をサンプルに添加した。
2.2 NSD2 は DNA が一部剥がれたヌクレオソームを結合する全体構造を図 2 中に示す。NSD2 は DNA が一部剝がれたヌクレオソームに結合し、それにより H3K36 へのアクセスを可能にしていた。溶液中でヌクレオソーム DNA の外部領域は結合と乖離を繰り返しており（breathing）、DNA が乖離した状態においてのみ NSD2 が結合すると考えられる。
この構造はリンカーヒストン H1 による NSD2 阻害機構を説明できる。図 2 右に示す通り、H1
はヌクレオソームから伸びるリンカーDNA を束ねるような部位に結合し、結果として DNA と
ヒストンとの相互作用を安定化する [14]。これにより、DNA の breathing は起こりにくくなり、
NSD2 による H3K36 へのアクセスが妨げられると考えられる。
2.3 ヒストンと DNA との相互作用
NSD2 はヒストン H3、H2A、そして DNA と直接相互作用を形成していた（図 3）。特に、
H2AK119 と NSD2 の E1031 の間で塩橋の形成が見られた。H2AK119 はモノユビキチン化を受けて条件的ヘテロクロマチン形成を介した転写抑制に働くことが知られている。一方で、H2AK119
がモノユビキチン化されると、E1031 との塩橋は形成できなくなる。これらの観察から、H2AK119のモノユビキチン化は H3K36 メチル化の阻害を介して転写抑制に関与するという新たな機構の存在が示唆された。
通常のヌクレオソーム構造において H3R49 と H3R52 は DNA リン酸基と塩橋を形成している。
今回の NSD2 複合体構造においては、NSD2 の主鎖あるいは側鎖の酸素原子が H3R49、H3R52 の側鎖と水素結合を形成していた。これらの相互作用は DNA が剥がれることによるエネルギー損
失を補償しているのであろう。

図 3. NSD2 と H2A・H3・DNA との相互作用。紫：NSD2、黄：H2A、青：H3、AtomColor：DNA。
2.4 自己阻害ループの構造変化先行研究で H3 フリーの NSD2 触媒ドメインの X 線結晶構造が解析されていた（図 4 左）[15]。
この構造において、NSD2 の「自己阻害ループ」と呼ばれる領域（オレンジ色で示す）は閉じた構造をとり、本来 H3 を結合すべき溝を埋めるように位置していた。自己阻害ループ上の D1182
は K1152 と分子内で塩橋を形成しつつ T1150 とは水を介した水素結合を形成し、一方で L1184
は M1119、T1121、I1127 などと疎水性相互作用を形成していた。これらの相互作用は自己阻害
ループの閉じた構造の安定化に寄与している可能性がある。
H3 フリーの構造と比較すると、ヌクレオソーム結合 NSD2 構造においては自己阻害ループが大きく動いて開いた構造をとり、それによって H3 の N 末端領域（H3A31-H3Y41）の結合が可能になっていた（図 4 右）。メチオニン置換が導入された H3K36M の側鎖は活性中心に伸びてsinefungin の近傍に位置しており、本構造はメチル化反応における複合体構造をよく表現した構- 11 -

SAR News No.46 (Apr. 2024)

造になっていると考えられる。
H3V35 と H3P38 は NSD2 の疎水性残基によって認識されており、
これが H3K36（または H3K36M）側鎖の活性中心への適切な配向に重要と考えられる。
2 か所の発がん性変異残基である E1099 と T1150 は共に活性中心と自己阻害ループの近傍に位置しており、H3 フリー構造において E1099 は K1124 と塩橋を形成していた。本研究のヌクレ
オソーム複合体構造において E1099K の側鎖は DNA リン酸基の近傍に位置しているが、その密度は不明瞭であり、DNA リン酸基との塩橋形成は確認できなかった。

図 4. NSD2 の H3 フリー構造（左）とヌクレオソーム結合構造（右）の比較

3. 変異体の生化学的解析
変異が導入された NSD2 触媒ドメインを発現・精製し、そのヌクレオソーム上 H3K36 メチル
化活性を測定した（図 5）。既報の通り、発がん性変異である E1099K と T1150A は共にメチル化
活性を大きく向上させた。興味深いことに、E1099 を Arg・Glu・Ala に変異させた変異体におい- 12 -

SAR News No.46 (Apr. 2024)ても程度の差こそあれ同様にメチル化活性の上昇が認められた。この事実は、E1099K の側鎖とDNA リン酸基の間の塩橋形成による親和性向上では説明がつかない活性亢進機構が働いている
ことを示す。
活性亢進機構について更なる知見を得るため、野生型・E1099K・T1150A 変異体について kinetics
解析とヌクレオソーム結合アッセイを行った（表 2）。その結果、両変異体においてヌクレオソーム親和性はあまり変化しておらず、一方で kcat が大きく上昇していることが明らかになった。

図 5. ヌクレオソーム上 H3K36 メチル化アッセイ

タンパク質
野生型
E1099K
T1150A

表 2. 変異体の kinetics 解析と結合アッセイ
kcatappkcatapp/KM
ヌクレオソーム-1(min )(µM-1min-1)
KM (µM)
0.14 ± 0.01
0.242 ± 0.043
0.57
2.16 ± 0.07
1.22 ± 0.03

0.232 ± 0.027
0.167 ± 0.018

9.33
7.35

ヌクレオソーム
KD (µM)
0.088 ± 0.031
0.202 ± 0.050
0.070 ± 0.028

4. MD シミュレーション
表 3. MD トラジェクトリ中で自己阻害ループが開く割合
システム
野生型
E1099K
T1150AE1099K-T1150A

開いた状態 (%)
0.86
9.7
6.3
18.8

図 6. 自己阻害ループの状態判定に用いた距離
構造解析と変異体解析の結果から、「E1099K と T1150A が自己阻害ループの構造ダイナミクスに影響を与えることで kcat を上昇させている」という仮説を我々は考えた、その仮説を検証する- 13 -

SAR News No.46 (Apr. 2024)

ため、NSD2 触媒ドメインの野生型・E1099K・T1150A・E1099K-T1150A（二重変異体）の 4 種類について MD シミュレーション解析を行った。すべてのシステムについて、MD トラジェクト
リ中で自己阻害ループはフレキシブルに動いていたが、その動き方は大きく異なっていた。解析のために、図 6 のように自己阻害ループ状の 2 つの疎水性残基（L1181 と L1184）と別の残基について 4 種類の距離（D1~D4）を計測し、それらが全て設定された閾値より大きい場合は自己阻害ループが「開いた」状態、そうでない場合は「閉じた」状態であると定義した（図 6）。この定義を用いて判定したところ、野生型 NSD2 の自己阻害ループは開いた状態をトラジェクトリ中でほとんどとらなかった（約 1％未満）のに対し、3 種類の変異体すべてで開いた状態をとる頻度が大幅に増加していた（表 3）。

5. Discussion
立体構造解析・変異体解析・MD シミュレーションの結果から、発がん性変異による活性亢進機構として我々は以下を提唱する。野生型 NSD2 においては、自己阻害ループはフレキシブルな構造をとりつつ H3 結合部位を覆い、そのことで過剰な H3K36 メチル化を防いでいる。一方で、
E1099K あるいは T1150A 変異は自己阻害ループの運動性に影響を与え、それが完全に開いた状
態を取りやすくさせる。結果としてヌクレオソームと相互作用した NSD2 が H3K36 を活性中心に結合してメチル化反応を起こす効率は上昇し（kcat の上昇に反映される）、メチル化活性の亢進が引き起こされる。
発がん性変異が自己阻害ループの運動性に影響を及ぼす機構の詳細な理解は今後の課題であるが、トラジェクトリの解析から、活性中心と自己阻害ループの近傍で複雑な塩橋と水素結合の動的なネットワークが形成されており、E1099K や T1150A においてはそれらの形成頻度が変化
していることが分かった。また、がん患者由来の点変異に加えて、新たな活性亢進変異として
我々は T1121A を同定した（図 5）。H3 フリー構造において T1121 は自己阻害ループの近傍で水素結合を形成しており、T1121A は発がん性変異と同様に水素結合ネットワークを変化させることで自己阻害ループの運動性に影響しているのかもしれない。今後、大規模な変異体解析に加えてヒストンペプチドやヌクレオソームを含む巨大な系を用いたシミュレーション計算を行うことで自己阻害ループによる NSD2 活性の制御とその破綻が詳細に理解されると期待される。

謝辞
立体構造解析と生化学的解析は筆者と佐藤光博士（横浜市立大学）が中心となって行い、東京大学の濡木理教授と西澤知宏博士（現・横浜市立大学教授）に解析上の助言を受け、クライオ電顕装置として東京大学医学部の Titan Krios を用いた。MD シミュレーションは量子科学技術研究
開発機構の河野秀俊チームリーダーと Amarjeet Kumar 博士らが行った。本研究は JSPS 科研費の助成（JP20K08717、JP18H05534、JP20H05394、JP21H05161）と AMED BINDS の支援（P20am0101115、
JP20am0101106）を受けた。

参考文献[1] Millán-Zambrano, G., Burton, A., Bannister, A. J., et al., Histone post-translational modifications —
cause

and

consequence of genome function, Nat Rev Genet. 23, 563–580 (2022).

https://doi.org/10.1038/s41576-022-00468-7[2] Feehley, T., O’Donnell, C. W., Mendlein, J., et al., Drugging the epigenome in the age of precision
medicine, Clin Epigenet. 15, 6 (2023). https://doi.org/10.1186/s13148-022-01419-z[3] Jambhekar, A., Dhall, A., Shi, Y., Roles and regulation of histone methylation in animal development,
Nat Rev Mol Cell Biol. 20, 625–641 (2019). https://doi.org/10.1038/s41580-019-0151-1[4] Bennett, R. L., Swaroop, A., Troche, C., et al., The role of nuclear receptor–binding set domain family- 14 -

SAR News No.46 (Apr. 2024)

histone lysine methyltransferases in cancer, Cold Spring Harb Perspect Med. 7, a026708 (2017).
https://doi.org/10.1101/cshperspect.a026708[5] Bergemann, A. D., Cole, F., Hirschhorn, K., The etiology of Wolf–Hirschhorn syndrome, Trends Genet.
21, 188–195 (2005). https://doi.org/10.1016/j.tig.2005.01.008[6] Keats, J. J., Maxwell, C. A., Taylor, B. J., et al., Overexpression of transcripts originating from the
MMSET locus characterizes all t(4;14)(p16;q32)-positive multiple myeloma patients, Blood. 105,
4060–4069 (2005). https://doi.org/10.1182/blood-2004-09-3704[7] Oyer, J. A., Huang, X., Zheng, Y., et al., Point mutation E1099K in MMSET/NSD2 enhances its
methyltranferase activity and leads to altered global chromatin methylation in lymphoid malignancies,
Leukemia. 28, 198–201 (2014). https://doi.org/10.1038/leu.2013.204[8] Bea, S., Valdes-Mas, R., Navarro, A., et al., Landscape of somatic mutations and clonal evolution in
mantle cell lymphoma, Proceedings of the National Academy of Sciences. 110, 18250–18255 (2013).
https://doi.org/10.1073/pnas.1314608110[9] Willcockson, M. A., Healton, S. E., Weiss, C. N., et al., H1 histones control the epigenetic landscapeby local chromatin compaction, Nature. (2020). https://doi.org/10.1038/s41586-020-3032-z[10] Yusufova, N., Kloetgen, A., Teater, M., et al., Histone h1 loss drives lymphoma by disrupting 3d
chromatin architecture, Nature. (2020). https://doi.org/10.1038/s41586-020-3017-y[11] He, L., Cao, Y., Sun, L., NSD family proteins: Rising stars as therapeutic targets, Cell Insight. 3,
100151 (2024). https://doi.org/10.1016/j.cellin.2024.100151[12] Sato, K., Kumar, A., Hamada, K., et al., Structural basis of the regulation of the normal and oncogenic
methylation of nucleosomal histone H3 Lys36 by NSD2, Nat Commun. 12, 6605 (2021).
https://doi.org/10.1038/s41467-021-26913-5[13] Herz, H.-M., Morgan, M., Gao, X., et al., Histone H3 lysine-to-methionine mutants as a paradigm tostudy chromatin signaling, Science. 345, 1065–1070 (2014). https://doi.org/10.1126/science.1255104[14] Bednar, J., Garcia-Saez, I., Boopathi, R., et al., Structure and dynamics of a 197 bp nucleosome incomplexwithlinkerhistone

H1,

Molecular

Cell.

66,

384-397.e8(2017).

https://doi.org/10.1016/j.molcel.2017.04.012[15] Tisi, D., Chiarparin, E., Tamanini, E., et al., Structure of the epigenetic oncogene MMSET and
inhibition by N-alkyl sinefungin derivatives, ACS Chem. Biol. 11, 3093–3105 (2016).
https://doi.org/10.1021/acschembio.6b00308- 15 -

SAR News No.46 (Apr. 2024)///// SAR Presentation Award /////SAR Presentation Award について
「SAR Presentation Award」は、構造活性相関シンポジウムにおける若手研究者の発表を奨励し、構造活性相関研究の発展を促進するため、2010 年度に創設された。当初は応募制として審査対象講演の募集を行った。2012 年度からは、正式名称を「構造活性相関シンポジウ
ム優秀発表賞」（英語表記 SAR Presentation Award）と定めた。
2023 年度 SAR Presentation Award について
2023 年度は、第 51 回構造活性相関シンポジウムにおける 40 歳以下の発表者（日本薬学会会員または受賞後に日本薬学会に入会いただける方）による一般講演（口頭発表・ポスタ
ー発表）を選考対象とすることとした。
2023 年度 SAR Presentation Award 受賞者（演題番号順）
口頭発表

： 吉澤

竜哉

（横浜市立大学生命医科学研究科）

ポスター発表： 立石優輔（熊本大学大学院自然科学教育部）

ポスター発表： 鍋谷朋哉（横浜市立大学理学部理学科）

ポスター発表： 平尾

巧（筑波大学大学院フロンティア医科学学位プログラム）

受賞者の選考について
2023 年 11 月 10‐11 日に各審査員からオンラインにて提出いただいた審査票を集計し、口頭発表 1 名、ポスター発表 3 名を受賞者として選出した。口頭発表の審査は点数方式、ポスター発表の審査は 3 演題選出するという方式で行った。審査員が共同演者となっている
演題は選出不可としている。後日受賞者には、賞状と副賞を贈呈した。なお、審査にあたっての審査項目は下記の通りである。
審査項目a) 講演要旨: 講演要旨は発表内容を反映して適切に作成されているか。
b) 講演資料: スライドは専門領域の異なる参加者にもわかりやすく、見やすく、かつ発
表時間に見合って適切に作成されているか。
c) プレゼンテーション: 発表時に参加者にわかりやすく説明しているか。
d) 研究の目的: 研究の背景と目的、先行研究との関係、研究の新規性あるいは有用性が明確になっているか。
e) 研究成果: 価値のある成果が得られているか。
f) 質疑応答: 質問等に対し、的確な応答・議論がなされたか。活発な討論がなされたか。
g) 将来性: 研究内容について、将来の発展が期待できるか。

- 16 -

SAR News No.46 (Apr. 2024)

審査員
第 51 回構造活性相関シンポジウムに参加した 2023 年度常任幹事および幹事- 17 -

SAR News No.46 (Apr. 2024)＜受賞者コメント＞
KO-03
氏名
所属
演題

吉澤 竜哉（よしざわ たつや）
横浜市立大学 生命医科学研究科
複数の特性予測モデルの信頼性を考慮した分子の多目的最適化

この度は名誉ある賞を頂き、大変光栄に存じます。評価をしていただいた先生方、日本薬学会
構造活性相関部会の関係者の皆様、シンポジウム実行委員の皆様に厚く御礼申し上げます。
本研究では、分子の特性を予測する予測モデルと分子生成 AI を組み合わせた多目的最適化における、特性予測の信頼度の問題に取り組みました。予測の信頼度は、予測モデルの適用範囲を定義することで考慮できるものの、高い信頼度で定義された複数の適用範囲が重複するとは限りません。さらに、適用範囲が重複するかどうかを、分子生成を実施する前に判別することが困難です。そこで本研究では、特性予測モデルの信頼度を探索的に調整することで、予測の信頼度を保ちつつ多目的最適化を実施する手法を開発しました。本手法を用いて Epidermal growth factor
receptor (EGFR) 阻害剤の設計を試みたところ、上市薬の一つを再設計することに成功しました。
本研究が、実践的な分子設計における分子生成 AI の活用の促進に貢献できることを期待しています。
今回のシンポジウムでは、オンサイトでの開催が再開されたことで、私が過去に参加した際には実現しなかった、創薬がご専門の先生方と直接議論をする機会を得ることができました。また、その議論を通じて自身の研究の価値や課題を改めて認識することができました。この貴重な経験を糧に、今後とも研究に精進します。最後に、本研究を遂行するにあたりご指導賜りました、
横浜市立大学 寺山慧准教授、石田祥一特任助教、共同研究者の先生方に、この場を借りて御礼申し上げます。

KP-06
氏名
所属
演題

立石 優輔（たていし ゆうすけ）
熊本大学大学院自然科学教育部
電子状態インフォマティクスによる-グルコシダーゼ阻害剤の探索

この度は、第 51 回薬学会構造活性相関シンポジウム SAR Presentation Award（ポスター賞）という身に余る賞を頂きましたこと大変光栄に存じます。当研究では、糖尿病治療薬の一つである
-glucosidase 阻害剤の新規探索を目的として、天然物からの in silico スクリーニングを行いました。探索に使用した電子状態インフォマティクス（ESI）手法は、電子状態的な記述子を用いた手法であり、電子状態的な観点で物性や活性などを予測、解析したり、化合物探索を行うアプローチです。ESI 手法を応用し、阻害剤の探索を行った結果、天然物から有望な-glucosidase 阻害剤候補を発見することが出来ました。加えて、ESI 手法により広範な化学空間から化合物を探索可能であることや、それにより従来の構造依存的な探索手法では、手の届かなかった未知の阻害- 18 -

SAR News No.46 (Apr. 2024)

剤候補が発見される可能性も示唆されました。これらの成果は天然物の新たな価値を見出す上で非常に有意義であり、また化学構造に制約されない新たな薬剤探索手法として、さらなる研究開発への道を拓くものと考えております。
このような研究成果を挙げることが出来たのは、杉本学准教授をはじめとする研究室メンバ
ー、そして日々の研究活動におけるすべての方々の指導・助言の賜物です。今後もこの受賞を励
みに、さらなる研究成果の創出に努めて参ります。最後に、研究発表の機会を提供してください
ました実行委員会の先生方に心から感謝申し上げます。

KP-13
氏名
所属
演題

鍋谷 朋哉（なべたに ともや）
横浜市立大学 理学部理学科
スーパーコンピュータ富岳を用いた MD ベースのハイスループットスクリーニング

この度は第 51 回構造活性相関シンポジウム SAR Presentation Award (ポスター)を賜りましたこと、大変光栄に思います。ご評価くださいました審査員の先生方、並びに日本薬学会構造活性相
関部会の先生方に心より感謝申し上げます。本シンポジウムでは多くの専門家の方々から研究
内容についてご意見やご質問を頂き、様々な知見や指針を得ることができました。
本研究では、ドッキングポーズにおけるリガンドの歪みによるフィルタリングや、スーパーコ
ンピュータ「富岳」を用いた 3000 回の MD シミュレーションによるドッキングポーズ安定性に基づくフィルタリング、さらには Absolute Free Energy Perturbation（ABFEP）法を用いた結合自由エネルギー計算を実施し、MD シミュレーションを用いた効率的なスクリーニングプロトコルの開発を行いました。結果、本研究で用いたデータセットの範囲において、開発した MD シミュ
レーションベースのプロトコルでは既存の静的構造ベースのプロトコルに比べて、ヒット率が18 倍程度向上することが明らかになりました。今後は、本手法を他のデータセットにも適用し、様々な標的タンパク質にも適用可能なインシリコスクリーニングプロトコルとして確立したいと考えています。
最後に、本研究を進めるにあたりご指導賜りました、池口満徳教授をはじめ、横浜市立大学生
命情報科学研究室の皆様に心より感謝申し上げます。

KP-21
氏名
所属
演題

平尾 巧（ひらお たくみ）
筑波大学大学院人間総合科学研究群フロンティア医科学学位プログラム
バーチャルスクリーニングモデル作成のための Parallel Cascade Selection Molecular
Dynamics の応用

この度は、第 51 回構造活性相関シンポジウム SAR Presentation Award（ポスター）という名- 19 -

SAR News No.46 (Apr. 2024)

誉ある賞を頂戴し、大変光栄に存じます。ご評価いただいた先生方、日本薬学会構造活性相関部
会の関係者の方々、シンポジウム実行委員の方々に厚く御礼申し上げます。
Structure-Based Virtual Screening (SBVS, VS) は、標的タンパク質の構造情報に基づいて低分子化合物を評価し、標的タンパク質に結合する可能性の高い化合物を選択する手法です。
VS には、
実験構造やホモロジーモデルによる予測構造が用いられ、VS の精度はこれらの構造に大きく依存することが知られています。近年では、AlphaFold2（AF2）によって、VS に利用可能な 3 次元タンパク質構造の潜在的な数は急速に増加していますが、AF2 によって予測された構造は、ドッ
キングシミュレーションに不適切であることが、いくつかの先行研究によって報告されています。本研究では、構造サンプリング手法である Parallel Cascade Selection Molecular Dynamics (PaCSMD)を用いて、標的タンパク質の AF2 構造から VS に適した構造に精密化する手法を提案しました。本手法を複数のターゲットに適用したところ、ベンチマークスクリーニングの精度向上が見られました。今後の展望として、さらなるデータセットでの検証やタンパク質-リガンド間相互作用評価の導入を検討しています。さらに、本研究をもとに創薬研究における VS の精度向上に貢献していければと考えています。
最後に、本研究を遂行するにあたりご指導賜りました筑波大学医学医療系の広川貴次教授をはじめ、吉野龍ノ介助教、筑波大学計算科学研究センターの原田隆平准教授、筑波大学物理学専
攻の工藤玄己さんに深く感謝申し上げます。

- 20 -

KO04

複数の特性予測モデルの信頼性を考慮した分子の多目的最適化○吉澤竜哉 1、石田祥一 1、佐藤朋広 2、大田雅照 3、
本間光貴 2、寺山慧 1,4,5（1 横浜市大•生命医、2 理研•BDR、3 理研•R-CCS、
4
理研•AIP、5 東工大•MDXES）
E-mail: w225435a@yokohama-cu.ac.jp

1. 背景と目的

2.2. 信頼性を 調整

医薬品候補の探索においては、疾患に対する薬
効や、体内動態等を含む多数の特性が同時に最適化された分子が要求される[1]。この探索は多目的最適化問題としてみなすことができ、近年では分子構造生成 AI を用いた分子特性の多目的最適化が試みられている[2-5]。
分子構造生成では、生成された分子構造の各特性は、教師あり学習に基づく予測モデルで評価されることが多い[6]。ただし、予測モデルの適用範
囲は限定されているため、生成構造に対して常に信頼性の高い予測を行えるとは限らない。そこで予測の信頼性を担保するために、予測モデルの適用範囲内で構造生成を実施する手法が提案されている[3,7,8]。しかし、これらの手法では単一の予測モデルの適用範囲のみに焦点が当てられており、多目的最適化において使用される複数の予測モデルの適用範囲は考慮されていなかった。予測モデルの適用範囲は、複数の予測モデル間で重複するとは限らないため、複数の予測モデルの信頼性を同時に向上させることは困難である。そこで本研究では、複数の予測モデルの信頼性それぞれ調整し、可能な限り信頼性の高い予測結果に基づき多目的最適化を行う手法を開発した。

2. 方法
2.1. ベイズ最適化による信頼性の調整
本研究では、予測モデルごとの信頼性の調整に、ベイズ最適化[9]を活用した。ベイズ最適化に基づ
き、信頼性の調整と構造生成を繰り返し実施することで、モデルごとに適切な信頼性を設定することを試みた。ベイズ最適化の計算は、予測モデルの信頼性の設定、構造生成、目的関数の計算の三段階で構成した(図１)。それぞれ詳細を 2.2–2.4 に記す。
2.2 予測モデルの信頼性の設定
本研究では、予測の信頼性として、予測モデルごとに適用範囲を定義し、適用範囲の厳しさを変2.3. 適用範囲内で構造を 生成

2.4.
目的関数
を 計算

図 1. 提案手法の計算の流れ 3 つのステップを繰り返すことで予測モデルの信頼性を調整しながら分子構造を生成する。
動させることにより調整した。適用範囲の定義は、
既存研究[3]にならい、訓練データとの最近接タニ
モト類似度が与えられた閾値を上回るか否かとした。適用範囲の厳しさは、訓練データとの最近
接タニモト類似度の閾値の大きさとした。
2.3. 適用範囲を考慮した構造生成
構造生成 には、 深層学習 手法の 一種である
recurrent neural network と、強化学習に基づく探索
アルゴリズムである Monte-Carlo tree search で構成される ChemTSv2[10]を使用した。ChemTSv2 は、予め定義した報酬関数を最大化するように、分子構造を生成する。報酬関数は、予測モデルの適用範囲内で特性が最適化されると値が大きくなる
ように設計した。また、適用範囲内で特性を最適化するために、生成構造が前項で設定された適用範囲を満たすか否かで場合分けをした。ここでは、生成構造と各予測モデルの訓練データの最近接
タニモト類似度が前項で設定された閾値を全て上回った場合、生成構造が適用範囲内である判定
し、報酬として各特性値の幾何平均を計算した。
各特性値は、予測値の望ましさに応じて 0 から 1に正規化し、報酬計算に用いた。最近接タニモト
類似度が閾値を満たさなかった場合には、適用範囲外と判定し、報酬は 0 とした。このような報酬関数を設計することで、適用範囲内での分子特性の多目的最適化を目指した。

- 21 -

2.4 ベイズ最適化の目的関数の計算
目的関数は、構造生成の結果と、入力したタニ
モト類似度の閾値の望ましさの二つの要素で構成した。構造生成の結果は、生成された構造の報
酬の上位一割の平均値で評価した。入力した閾値
の望ましさは、閾値が高いほど 1 に近付くように
0 から 1 に正規化し、正規化後の値で評価した。
それらの重み付き幾何平均をベイズ最適化の目的関数とした。

A-1 (Gefitinib)

A-2

B-1

B-2

図 2. 生成構造例 (A)EGFR に対する阻害活性を有する既知阻害剤と同一な構造(B)報酬関数の値が上位の構造
表 1. 生成構造例の各特性の予測値及び訓練データとの最近接タニモト類似度
ID は図 2 と対応している。EGFR に対する阻害活性、代謝安定性、膜透過性の予測値の単位は、それぞれ pIC50、1 hour remaining(%)、log10(μcm/s)である。

2.5. ベイズ最適化及び構造生成の計算条件
ベイズ最適化の探索回数は 40 回とした。またベイズ最適化の効率を評価するために、ランダム探索との比較を行った。構造生成において最適化する特性としては、epidermal growth factor receptor(EGFR)に対する阻害活性と、ヒト肝ミクロソームにおける代謝安定性、Caco-2 細胞における膜透過性の 3 項目を標的とした。1 回の探索当たりの生成構造数は 10,000 個とした。

EGFR活性

代謝安定性

特性

類似度

特性

類似度

特性

類似度

A-1

7.9

1.0

55

1.0

1.1

0.43

3. 結果とまとめ

A-2

7.7

1.0

60

0.73

1.3

0.73

ランダム探索とベイズ最適化における 40 回の探索で得られた目的関数の最大値を比較すると、ランダム探索では 0.20 だったものの、ベイズ最適化では 0.39 と、2 倍程度の高い値を得た。ベイズ
最適化により、目的関数の値を高くする閾値を、効率的に探索することに成功した。
ベイズ最適化の目的関数が最大の時の構造生
成では、生成の過程を通して、生成構造と訓練データとの最近接タニモト類似度が全特性において閾値を超えており、適用範囲の内側で分子構造が生成されたこと明らかになった。さらに生成構造の特性値は、全ての特性で正規化後の値が 0.7以上となっており、概ね最適化されていた。これ
らの結果から、予測モデルの適用範囲内で多目的最適化が達成されたことが示唆された。
図 2 に一部の生成構造、表 1 にそれらの特性の予測値、及び訓練データとの最近接タニモト類似度を示す。本手法により、上市薬の一種である
Gefitinib を含む既知 EGFR 阻害剤と同一の構造が生成された(図 2, 表 1, A)。また、それ以外にも望ましい特性を持つと予測された構造を複数生成することに成功した(図 2, 表 1, B)。
本研究では、複数の予測モデルの信頼性を同時に考慮して、分子特性の多目的最適化を実施する
ための手法を開発した。各予測モデルの信頼性を
ベイズ最適化に基づいて調整し、可能な限り信頼性が高い特性予測に基づいた多目的最適化を実施した。その結果、適用範囲を考慮した多目的最適化を達成し、EGFR 阻害剤候補の生成に成功した。本手法は、医薬品開発におけるリード最適化のプロセスへの応用が期待される。

B-1

8.1

0.73

68

0.58

1.1

0.49

B-2

7.4

0.81

71

0.58

1.4

0.58

ID

膜透過性

4. 謝辞
本研究は、AMED〈産学連携による次世代創薬
AI 開発（DAIIA）〉の課題番号 JP22nk0101111、
AMED BINDS の課題番号 JP22ama121023、「富岳」
で目指すシミュレーション・AI 駆動型次世代医療・創薬(JPMXP1020230120)の支援を受けて実施されたものである。

5. 参考文献[1] Nicolaou, C. A., et al. Drug Discov. Today
Technol. 10, e427–e435 (2013).
[2] Winter, R., et al. Chem. Sci. 10, 8016–8024 (2019).
[3] Perron, Q., et al. J. Comput. Chem. 43, 692–703(2022).
[4] Bung, N., et al. J. Chem. Inf. Model. 62, 2685–
2695 (2022).
[5] Yoshizawa, T., et al. J. Chem. Inf. Model. 62,
5351–5360 (2022).
[6] Luukkonen, S., et al. Curr. Opin. Struct. Biol. 79,
102537 (2023).
[7] Kaneko, H., et al. Mol. Inf. 36, 1600032 (2017).
[8] Langevin, M., et al. ACS Omega 8, 23148‒23167(2023).
[9] Shahriari, B., et al. Proc. IEEE 104, 148–175(2016).
[10] Ishida, S., et al. WIREs Comput. Mol. Sci., e1680(2023).

- 22 -

KP06

電子状態インフォマティクスによるa-グルコシ
ダーゼ阻害剤の探索○立石優輔 1、杉本 学 1-2（1 熊本大院、2 沼津高専グリーンアンモニア研）
E-mail: 235d8761@st.kumamoto-u.ac.jp

1. 背景と目的
糖尿病治療であるa-glucosidase 阻害剤の新規探索を目的とし、当研究室で開発している電子状態
インフォマティクス（ESI）手法 [1]を用いて天然物から、in silico スクリーニングを行った。
医薬品探索において、天然物の影響力は大きく、これまでに複数の承認薬を生み出してきたと共に、多様な化学構造を持つ薬剤分子の発見に貢献している [2]。ESI 法は、分子の構造的特徴に直接参照することなく、分子の電子的特徴に注目する
ため、化合物の探索に有用であると期待されている [3]。すなわち、既知の骨格とは類似しない多様な足場の発見 (Scaffold Hopping) に繋がる可能性がある。

2. 計算方法
ChEMBL データベース[4]からa-glucosidase 阻
害剤の構造と活性値(IC50)を得た。この構造を初期構造として配座探索、電子状態計算を行い、最終的に 433 分子の ESI 記述子を得た。
また従来の ESI
記述子に加えて、スピン 1 重項励起状態への励起
エネルギーΔEn (n = 1-30)の逆数を加えた。分子物性の記述に必要な 1 次の摂動波動関数が次式(1)で表されることから、予測精度と解釈性の向上を期待し、導入した。
(#)(#)&𝜓! '𝑋)' 𝜓# *(#)(#)
𝜓 = 𝜓! + %
𝜓#
(1)
∆𝐸!

高活性と予測された天然物の、a-glucosidase に対する結合能力を評価するために、Autodock vina
[7]を用いて、ドッキングシミュレーションを実施した。
3. 結果と考察得られた回帰モデルは、それぞれ IC50 の実験
値を良好に再現していた（表１）。ESI 記述子によるモデルは ECFP4 や RDKit 記述子による予測モ
表 1. ESI 記述子、ECFP4 および RDKit 記述子を用いた予測モデルの精度比較
ESI desc.
ECFP4
RDKit desc.
R2(train)
0.99380.98230.9924R2(test)0.80670.85390.8665デルと比較してやや劣るが、電子状態由来の ESIモデルが構造由来の予測モデルと同程度の予測精度を持つことがわかった。新規阻害剤を探索するために、ESI モデルを用いて天然物に対し insilico スクリーニングを行った。その結果、高活性と予測された上位の化合物は、紅茶などの含むフラボノイド類の Theasinensin A [8]や、タンニンのActive!%#予測モデルの構築には、教師あり機械学習手法の一つである XGBoost[5]を用いた。全データ(433 分子)を、訓練用(346 分子)、テスト用(87 分子)に分割し、訓練用データでモデル構築を行い、モデルの精度検証にはテスト用データを用いた。比較のために、ECFP4 や RDKit 記述子による予測モデルも作成した。
天然物（2649 分子）は富山大学が公開するデータベース(KampoDB)[6]から収集し、これらの天然物に対して in silico スクリーニングを行った。

ActiveTheasinensin A-8.21 kcal/mol(6.34)

Chebulagic acid-9.28 kcal/mol(6.32)

CID: 71448959-11.3 kcal/mol(6.42)

Procyanidin B-5 3,3'-di-O-gallate-11.4 kcal/mol(5.89)

図 1. スクリーニングにより得られた高活性化合物。
化合物の下部にはa-glucosidase への結合親
和性、括弧内は pIC50 の予測値である。

- 23 -

Almost all
steroidal compounds

Highly active
molecule groups
predicted by ESI,
mostly flavonoids

Most active
reported compounds

Small molecules

図 2. ECFP4 により描画したケミカルスペース。黄色は教師データのa-glucosidase 阻害剤、緑は天然物を表す。また ESI 記述子、ECFP4、および RDKit 記述子を用いた予測モデルにより高活性と予測された化合物の上位 50 化合物を、それぞれ赤、橙、青で示した。

一種である Chebulagic acid [9]などの既知の aglucosidase 阻害剤を含んでいた（図１）。また多くのフラボノイド類がa-glucosidase 阻害能を持つことが一般的に知られているが [10]、ESI モデルにより発見されたれた幾つかの天然物がフラボノイド類であることは非常に興味深い結果である。
スクリーニング結果の妥当性を検証するため、高活性と予測された上位の天然物に対し、 aglucosidase への結合親和性を評価したところ、多くの化合物が高い結合親和性を持つことがわかった。中でも CID: 71448959 や Procyanidin B-5 3,3'di-O-gallate（図１）はそれぞれ結合親和性が-11.3
kcal/mol、-11.4 kcal/mol と、上市薬の Acarbose の結合親和性(-8.0 kcal/mol)と比較して非常に高かった。以上の結果から、Theasinensin A などのように
文 献 に 報 告 の 無 か っ た 天 然 物 も 有 望 な aglucosidase 阻害剤の候補であること期待される。
また、ESI 記述子、ECFP4、RDKit 記述子による予測モデルでスクリーニングした結果、各予測モ
デルがそれぞれ個性的な探索範囲を持ち、中でもESI モデルはより広い化学空間を探索可能である
ことがわかった（図 2）。これは、ESI 手法によって構造情報に制約されない自由な薬剤探索が可能であることを支持する結果である。
電 子 状 態 を 考 慮 す る ESI 手 法 は 、 Scaffold
Hopping のための新たなアプローチとなり得る。

4. 参考文献
1.

T. Ideo, M. Sugimoto et al., Chem. Lett., 50, 849852(2021).
2. D. J. Newman, G. M. Cragg, J. Nat. Prod. 83, 770803(2020).
3. Y. Tateishi, M. Sugimoto, J. Comput. Chem.
Japan, 2023. (in press)
4. D. Mendez, et al, Nucleic Acids Res., 47, 930940(2019).
5. T. Chen, C. Guestrin, KDD’16: Proceedings of the
22nd ACM SIGKDD International Conference on
Knowledge Discovery and Data Mining, 785794(2016).
6. R. Sawada, et al, Sci. Rep., 8, 11216(2018).
7. J. Eberhardt, et al, J. Chem. Inf. Model., 61, 38913898(2021).
8. S. Tao et al., Food Funct., 11, 3527-3538(2020).
9. H. Gao et al., Biosci. Biotechnol. Biochem., 72,
601-603(2008).
10. Y. Jia et. al., J. Agric. Food Chem., 67, 1052110533(2019).

- 24 -

KP13

スーパーコンピュータ富岳を用いた
MD ベースのハイスループットスクリーニング○鍋谷 朋哉 1、浴本 亨 1,2、山根 努 3、池口満徳 1,2,3（1 横浜市大・理、2 横浜市大・生命医科学、3 理研・計算科学研究セ）
E-mail: s202076a@yokohama-cu.ac.jp

1. 背景と目的
ヒット化合物探索過程では、化合物ライブラリ
内の多数の化合物を標的タンパク質立体構造へ
ド ッ キ ン グ さ せ 、 経 験 的 な Docking score や
MMGBSA などの物理化学的スコアに基づいて候補化合物をランク付けし、アッセイ実験を実施する。これらの静的構造ベースのスクリーニングによるヒット率は一般に 1%程度[1]とされている。
生体内では、タンパク質は動的に揺らいでおり、静的な構造から得た結合ポーズが必ずしも正確である保証はなく、動的構造を考慮した化合物探索法が求められている。リード最適化過程においては、FEP+等の分子動力学（MD）シミュレーシ
ョン導入による数十化合物程度の検証が普及しているものの、数千～数万規模の化合物を扱うハイスループットスクリーニング(HTS)に MD を適用するには計算コストが高い。近年、富岳のような超大型スーパーコンピュータの登場により、数千、数万個の MD 実行が現実的な時間で可能になってきたため、MD ベースのインシリコスクリー
ニングが現実的になりつつある。しかし、MD を用いた HTS 手法はまだ確立されていない。そこ
で、本研究ではスーパーコンピュータ富岳を使った数千回の MD を実行し、MD ベースのフィルタリングや結合自由エネルギー計算を実施・既存指標との性能比較によって、有効な MD ベース HTS
手法を開発することを目的とする。

2. 方法
2.1.データセット
我 々 の 先 行 研 究 で 標 的 と し た UHRF1 (PDB
ID:5YY9)のインシリコスクリーニングで選択し、アッセイ実験を実施した 130 個の化合物群を用いた。2 つのヒット化合物[2]と、別グループから報告された、同ポケットに結合する別の化合物[3]も含め、合計でヒット化合物３つを含む、131 化合物を用いた。計算では、構造異性体を考慮したた
め、合計で 600 構造中 3 個の活性リガンドを含む
データセットとなった(図１)。

5A-DMP

2.4-lutidine

Compound2

図 1 UHRF1 及び活性リガンド
2.2.ドッキング及び MD シミュレーション
作成したデータセットの化合物全てについて
Glide[4]を使ってドッキングシミュレーションを行った。次に生成した 600 個のドッキングポーズ
について GROMACS(2021.4)を用いて、初速度の異なる 10 ns の MD を 5 回(計 3000 回)行った。

2.2. RMSD filter
ドッキングポーズの安定性を、短い MD によって評価する方法がヒット率向上に有効であると報告されている[2][5]。本研究では MD によって得られたトラジェクトリから、リガンド RMSD を計算し、5 回の MD のうち、1 回でもリガンド
RMSD が 2Åを下回るものがあれば、そのリガンドを候補化合物として選抜した。RMSD でフィルタリングした後、Docking score と MMGBSA でランキング化し、Enrichment が RMSD filter でどの程度向上するか検討した。
2.3. Strain energy filter
無理なコンフォメーションでドッキングされている結合ポーズが、ドッキングポーズにおけるリガンドの歪みの大きさで評価できると示唆されている[6]。具体的には、リガンドの Strain energy
で化合物をフィルタリングすれば、無理な結合ポーズの化合物が除去でき、ヒット率が向上すると報告されている[6]。本研究においても、Strainenergy を計算し、Total strain energy < 7 TEUs でフィルタリングした後、Docking score や MMGBSA
でランキング化して Enrichment がどの程度向上するか検討した。

- 25 -

2.4 リガンドの結合自由エネルギー評価
RMSD filter を 突 破 し た 化 合 物 に 関 し て 、GROMACS(2021.4)を用いて FEP 計算(Free Energy
Perturbation)を実施し、結合自由エネルギー(ΔG bind)を求めた。摂動計算は、複合体系、孤立系でそれぞれ 1ns ずつ行い、中間状態は複合体系
で 42λ、孤立系で 31λとした。また複合体系ではリガンドとタンパク質間の相対配置を拘束した[7][8]。ΔG bind 値の順にリガンドをランキング化
し、Docking score や MMGBSA でのランキングと性能比較した。

3. 結果とまとめ

本研究では生成したドッキングポーズに構造ベース及び結合自由エネルギーベースで評価を行
い 、 そ れ ぞ れ Docking score や MMGBSA とEnrichment を比較検討した。その結果、MD 中のドッキングポーズの安定性が低いリガンドを
RMSD filter を用いてフィルタリングすること、またドッキングポーズの歪みが大きいリガンドをStrain filter を用いて、フィルタリングすること、さらに FEP を用いて化合物をランキング化する
ことが、ヒット率の向上に有効であることが分かった。これらの手法を組み合わせてスクリーニン
グすると、既存の静的構造ベースの指標に比べ約10 倍程度、Enrichment が良化した。 (図 3)

RMSD filter 及び Strain filter で Filtering した後に、Docking score や MMGBSA でランキング化した時の Enrichment factor の推移を示す（表１，２）。
表 1 Docking score でランキングした時のEnrichment factor の推移
Enrichment factor

5%

10%

20%

Docking score

6.67

3.33

3.33

Docking score(RMSD filter)

6.67

6.67

5.00

Docking score(Strain filter)

13.33

10.00

5.00

表 2

MMGBSA score でランキングした際のEnrichment factor の推移

Enrichment factor

5%

10%

20%

MMGBSA

0.00

0.00

1.67

MMGBSA(RMSD filter)

0.00

3.33

5.00

MMGBSA(Strain filter)

6.67

10.00

5.00

RMSD filter の 導 入 に よ り 、 Docking score や
MMGBSA 単体でランキング化した場合に比べ、Enrichment factor が良化した(表１)。また Strain
filter の導入においても同様に、Enrichment factorが良化した(表２)。
次に RMSD filter を突破した 144 化合物について FEP 計算を行い、その後ΔG bind 値で化合物をランキング化した（図 2）。なお Compound 2 に関しては Thermal shift assay における変性温度が僅かで、Isothermal Titration Calorimetry 測定を実施しておらず、実験値が不明なため FEP 計算の対象にしなかった。FEP によるΔG bind 値のランキング
で、Enrichment が良化した。(図 2)

図 2 評価指標別 Enrichment curve

本研究においては結合実験データのある UHRF1をデータセットの対象としていたが、今後は他のデータセットについても今回の手法を適用し、様々な標的タンパク質にも適用可能な MD ベースのハイスループットスクリーニング手法の性能評価をしていきたいと考えている。

4. 謝辞
本研究の一部は、HPCI システム利用研究課題(課題番号：hp220236,hp230385(浴本))を通じて、理化学研究所が 提供するスーパーコンピュータ「富岳」の計算資源の提供を受け、実施しました。

5. 参考文献
1. Zhu T,et al.,J Med Chem.12;56(17):6560-72
2. Kori, S. et al., Bioorg. Med. Chem. 52, 116500(2021).
3. Chang,L.et al., Sci Rep 11, 1121 (2021)
4. Friesner, .et al., J. Med. Chem., 2006, 49, 6177–6196
5.Liu, K. et al., J. Comput. Aided Mol.Des 2017
Feb;31(2):201-211.
6Shuo Gu .et al., J. Chem. Inf. Model. 2021, 61, 9,
4331–4341
7. Aldeghi et al. 2016 Chem. Sci., 2016,7, 207-218
8. Alibay, I.et at al, Commun Chem 5, 105 (2022).

図 2 指標別 Enrichment curve (RMSD filter 後)- 26 -

KP21

バーチャルスクリーニングモデル作成のためのParallel Cascade Selection Molecular Dynamics の応用○平尾巧 1 工藤玄己 2 吉野龍ノ介 3,4 広川貴次 3,4（1 筑波大・医科学、2 筑波大・物理、3 筑波大・医、
4
筑波大・TMRC）
E-mail: s2221383@u.tsukuba.ac.jp

1. 背景と目的
Structure-Based Virtual Screening (SBVS)は、標的タンパク質の構造情報に基づいて低分子化合物
ライブラリーを評価し、標的タンパク質に結合する可能性の高い化合物を選択する手法である[1]。
VS には、実験構造やホモロジーモデルなどの標的タンパク質の 3 次元構造が用いられ、VS の精度はこれらの構造に大きく依存する。
近年、AlphaFold2（AF2）[2]の登場により、VSに利用可能な 3 次元タンパク質構造の潜在的な数は急速に増加している。しかし、これまでの研究
では、AF2 によって予測された構造は、しばしばドッキングシミュレーションに適さないことが報告されている [3]。また、分子動力学シミュレー
ションを用いた先進的なドッキング手法として、
ア ン サ ン ブ ル ド ッ キ ン グ [4] や Induced Fit
Docking[5]など、受容体の柔軟性を考慮してタンパク質構造を改良する手法が提案されているが、これらの手法では広い構造空間から最適な VS モ
デルを自動的に得ることには限界がある。
本研究では、複数の MD を繰り返し、適切な配置を選択するコンフォメーションサンプリング
法 で あ る Parallel Cascade Selection Molecular
Dynamics (PaCS-MD)[6]を用いて、様々な初期構造から標的タンパク質を VS に適した構造に精密化する手法を提案した。

アポ・ホロプロトコルのワークフローを図 1 に示す。評価用データセットを用いたベンチマーク
スクリーニングの結果を PaCS-MD のモデル選択に使用する。ベンチマークスクリーニングは、
ROC および AUC を用いて評価され、AUC の高い構造が次のサイクルへ採用される。ホロプロトコ
ルの場合、ベンチマークスクリーニングにおいて、ドッキングスコアがトップの化合物とターゲッ
トタンパク質の複合体として MD 計算を行う。
PaCS-MD のサイクルを両手法で 40 回行い、評価
用データセットを用いたベンチマークスクリー
ニングの AUC が最も高い構造を最終的なモデルとした。
Pre-MD for AF2 model(APO)

Model-1

Model-2•••

Model-10

Short MD (100 ps*2)
Traj-1Traj-2•••Traj-10Evaluation Screening (ranked by ROC/AUC)

cyc i*3

Select Models and Make Topology

Rank-1

Rank-2•••

Rank-102. 方法

図 1. 2 手法のフローチャート
ホロプロトコルの場合、リガンドを含んだ系を生成し、タンパク質-リガンドの複合体構造でMD を行う。

本研究では、PaCS-MD 法を用いたスクリーニン
グモデル探索法であるアポプロトコル、ホロプロ
トコルの 2 手法を提案する。はじめに、DUD-E[7]および先行研究[8]より、ターゲットタンパク質および活性化合物、デコイ化合物の情報を得た。ターゲットごとに、10 個の活性化合物、90 個のデコイ化合物を MOE の Diverse Dataset を用いて選択し評価用データセットとした。さらに、残りの活性化合物とデコイ化合物の比率が１対９となる
ように同様に選択し、最終的なベンチマーク用のテスト用データセットとした。

本プロトコルで得られたモデルを評価するた
めに、テスト用データセットを用いて、ベンチマ
ークスクリーニングを行なった。同様に、PDB に登録されている結晶構造、AF2 により生成されたモデル、conventional MD (cMD)により平衡化された構造に対してベンチマークスクリーニングを行なった。
結合ポーズの解析として、タンパク質-リガンド
間相互作用フィンガープリント(PLIF)から、類似度を評価する Tc50[9]を計算し、モデルの評価を行- 27 -

なった。

3. 結果とまとめ
各ターゲットにおけるテスト用データセット
を用いた、ベンチマークスクリーニングの AUC を
表 1 に示す。結晶構造が既知である４つのターゲットにおいて、ホロプロトコルから得られたモデル(ホロモデル)は、PDB、AF2、c MD 構造より高い AUC 値を記録した。また、結晶構造未知のVAChT において、アポ・ホロモデルでは AUC が向上し、スクリーニングに適した構造に最適化された。
表 1. ベンチマークスクリーニング結果の AUC
AMPC AKT1 CXCR4 HIVPR VAChTPDB 0.612 0.814 0.694 0.833
AF2
0.568 0.800 0.691 0.816 0.256
cMD 0.638 0.680 0.656 0.808 0.729
アポ 0.719 0.860 0.669 0.737 0.840
ホロ 0.713 0.855 0.853 0.840 0.832
各ベンチマークスクリーニング結果の結合ポーズおよび PLIF から算出した Tc50 の値を表 2 に示す。PDB、AF2、c MD 構造と比較して、アポ・ホロモデルは Tc50 が全般的に高い値を示す。
VAChT ではホロモデルの Tc50 が最も高く、ドッキングにより予測された各化合物の結合ポーズが類似していることが示唆される。実際のポーズにおいても、アポモデルと比較してスクリーニン
グの結合ポーズと MD 中のポーズが似た配向を示
している。
（図２）
表 2. ベンチマークスクリーニング結果の Tc50
AMPC AKT1 CXCR4 HIVPR VAChT
PDB
0.09
0.22
0.11
0.36
AF2
0.09
0.17
0.10
0.27
0.10
cMD
0.20
0.18
0.14
0.36
0.14
0.22
0.31
0.16
0.13
0.37
アポ
0.26
0.23
0.13
0.26
0.44
ホロ
Test screening

LigandHOLO

F335
Test screening
APO

W331

図 2. VAChT における結合ポーズ
F335/W331 はトランスポーターの輸送に関与する残基。

以上より、本プロトコルを用いて、AF2 モデルからスクリーニングに適したモデルを自動的にサンプリングすることができた。また、構造未知のターゲットに対しても有効であり、他のモデルよりもファーマコフォアを考慮した構造を得られることが示唆された。

4. 謝辞
本研究は、以下の支援のもと実施された。
• Program for Promoting Researches on the
Supercomputer Fugaku (JPMXP1020230120)• Research Support Project for Life Science and Drug
Discovery (BINDS, P23ama121029j0002)• KAKENHI (23K16987)• The Cygnus computational resources provided by
the Multidisciplinary Cooperative Research
Program at the Center for Computational Science(Project Code : CADD) at the University of Tsukuba.

5. 参考文献[1] Zhu H, et al. A Comprehensive Survey ofProspective Structure-Based Virtual Screening for
Early Drug Discovery in the Past Fifteen Years. Int J
Mol Sci. 23(24):15961, (2022).
[2] Jumper J, et al. Highly accurate protein structure
prediction with AlphaFold. Nature. 596(7873):583-589,(2021).
[3] Scardino V, et al. How good are AlphaFold models
for docking-based virtual screening?. iScience.
26(1):105920, (2022).
[4] Amaro RE, et al. Ensemble Docking in Drug
Discovery. Biophys J. 114(10):2271-2278, (2018).
[5] Miller EB, et al. Reliable and Accurate Solution to
the Induced Fit Docking Problem for Protein-Ligand
Binding. J Chem Theory Comput. 17(4):2630-2639,(2021).
[6] Harada R, et al. Parallel Cascade Selection
Molecular Dynamics (PaCS-MD) to generate
conformational transition pathway. J Chem Phys.
139(3):035103, (2013).
[7] Mysinger MM, et al. Directory of useful decoys,enhanced (DUD-E): better ligands and decoys forbetter benchmarking. J Med Chem. 55(14):6582-94,(2012).
[8] Kersten C, et al. Hic Sunt Dracones: MolecularDocking in Uncharted Territories with Structures from
AlphaFold2 and RoseTTAfold. J Chem Inf Model.
63(7):2218-2225, (2023).
[9] Sato M, et al. Extended template-based modelingand evaluation method using consensus of binding
mode of GPCRs for virtual screening. J Chem Inf
Model. 54(11):3153-61, (2014).

- 28 -

SAR News No.46 (Apr. 2024)///// Activities /////
第 51 回構造活性相関シンポジウム開催報告
日時： 令和 5 年 11 月 20 日（月）・21 日（火）
会場： 日本薬学会長井記念ホール（オンラインとのハイブリット開催）
主催： 日本薬学会構造活性相関部会
協賛： 情報計算化学生物学会(CBI 学会)、 日本農薬学会、日本バイオインフォマティクス学会、日本薬学会医薬化学部会
2023 年 11 月 20 日（月）～21 日（火）、日本薬学会長井記念ホールにおいて、第 51 回構造活性相関シンポジウムを開催いたしました。新型コロナ感染状況も落ち着き、4 年ぶりに対面をメインとした開催（オンラインとのハイブリット）となりました。招待者と実行委員などを含めて
約 138 名の参加をいただきました．その内、オンライン参加が 25 名でしたので、100 名以上の方に会場にお越しいただけたことになり、参加者の皆様が対面開催を望んでいたことが伺えました。初日のセッションの終了後には、日本薬学会長井記念ホールの受付前スペースを利用して懇親会（ケータリングスタイル）も開催し、久しぶりの再会等もあって大いに歓談を楽しんでいた
だきました。
演題数は、特別講演が 1 件、招待講演が 3 件、一般講演として口頭が 11 件、ポスターが 40件となりました。ポスター発表（フラッシュトークは配信）以外は、全てオンライン配信をいたしました。特別講演では、理化学研究所の大田雅照先生に、「変化、進化、深化に適応し、促進するための技術としてのインシリコ創薬」という演題でご講演をいただきました。大田先生は、企業とアカデミアの双方において創薬の実績をお持ちであり、これまでのインシリコ創薬の経緯、
AI 創薬が主流となっている現状、今後の方向性について、先生ご自身の研究成果を交えた示唆に富んだ内容でした。招待講演では、東京大学大学院の森脇由隆先生、名古屋大学大学院の山西芳裕先生、横浜市立大学の仙石徹先生にご講演いただきました。森脇先生には、「AlphaFold に
よってさらに発展するタンパク質の計算科学」という演題でご講演をいただきました。森脇先生
は AlphaFold(Colab Fold)の開発の共著者として唯一の日本人であり、AlphaFold を深く理解し活用されているトップランナーとして、AlphaFold をカスタマイズする技術や問題点など参加者の多くの方にとって関心の高い実務レベルの内容のご講演をいただきました。山西先生には、「AI・ビッグデータ時代の創薬の可能性」という演題でご講演いただきました。ご講演内容は、疾患に関するマルチオミクス情報・臨床情報、化合物に関する化学構造・表現形・標的分子情報などの多様なデータを融合解析し、創薬における様々な課題を解決すための機械学習アルゴリズ
ムの開発と、その技術を疾患・生体分子・制御化合物ネットワークを大規模に予測する問題へ適用した事例など、非常に多岐にわたる内容で、参加者の関心を集めていました。仙石先生には、「転写制御と環状ペプチドの構造生物学」という演題でご講演いただきました。現在、クライオ
電子顕微鏡技術の進展により、構造創薬の重要度が高まっています。仙石先生には、転写因子
Nrf2 による標的 DNA の配列特異的認識の構造基盤、ヒストン H3K36 メチル化酵素 NSD2 のヌ- 29 -

SAR News No.46 (Apr. 2024)

クレオソーム認識と発がん性点変異による活性亢進の構造基盤、アゾール・ピリジン閉環骨格を持つ de novo チオペプチドによる TNIK キナーゼ阻害の構造基盤、など多くの構造生物学研究の最新成果についてご講演をいただきました。
一般講演でも AlphaFold 予測構造や、化合物生成 AI を活用した演題や、PROTAC など新しい
モダリティを標的とした演題が多く、創薬プロセスとしてパラダイムシフトが起きている印象を持ちました。ポスター発表は、会議室 A・B を拡張して行いました。40 演題と議論する参加者のスペースでちょうど収まる広さでした（今後、長井記念ホールを利用する方のために記録として残しておきます）。ポスター発表時間の会場は、どのポスターでも質問や議論をする参加者で賑わっており、対面開催をしてよかったと実感するものでした。また、シンポジウム開催前に、長井記念ホールの Wifi 環境も新しくなったこともあり、参加者にとって通信環境も最適でありました。

今回のシンポジウムでは、対面形式や懇親会開催の決断についてかなり悩みましたが、開催期
間中、深刻な感染状況の再発もなく、また会場での参加者間の対面での議論や歓談の様子を拝見し、対面開催（オンラインハイブリット開催）にして本当に良かったと胸を撫でおろしています。
開催にあたり、実行委員の先生方には多大なるご尽力いただきした。アルバイトの皆様も当日、的確にサポートいただきました。また小田先生、高木先生、杉本先生には、これまでのご経験からのアドバイスを多くいただきました。この場をかりてお礼申し上げます。
SAR Award は、以下の４名の発表者に決定いたしました。受賞おめでとうございます。
・KO04 吉澤竜哉 （横浜市大・生命医科学）「複数の特性予測モデルの信頼性を考慮した分子の多目的最適化」（口頭）・KP06 立石優輔
（熊本大院）「電子状態インフォマティクスによるα-グルコシダーゼ阻害剤の探索」（ポスター）
・KP13 鍋谷朋哉
（横浜市大・理）「スーパーコンピュータ富岳を用いた MD ベースのハイスループットスクリーニング」（ポスター）
・KP21 平尾巧
（筑波大・医科学）「バ ーチャ ルス クリー ニン グモデ ル作 成のた めの Parallel Cascade Selection Molecular
Dynamics の応用」
（ポスター）

第51回構造活性相関シンポジウム
実行委員長 広川- 30 -

貴次

SAR News No.46 (Apr. 2024)///// Activities /////＜会告＞構造活性フォーラム 2024
「分子シミュレーション・実験・機械学習が導く創薬の未来」
主催:
日本薬学会構造活性相関部会
協賛・後援: 計算情報科学生物学会（CBI 学会）
会期:2024 年 6 月 21 日 (金)会場:
Zoom によるオンライン開催
フォーラムホームページ： https://sites.google.com/view/qsarforum2024（4 月公開予定）
開催趣旨：近年、分子シミュレーション、機械学習、実験がそれぞれ相互に連携しあい、高度な研究開発が進められている。創薬分野においても AI 創薬など非常に盛んな取り組みがなされているが、材料・化学・生物分野においても様々な研究が展開されている。そこで、本フォーラム
では創薬に限らず、様々な分野での先端的な取り組みをご紹介いただき、創薬分野への適用の可能性や今後の新たな切り口を検討・議論したい。
プログラム：
講演１． 「化学研究のための大規模言語モデルの活用と構築の試み」
畠山 歓（東京工業大学）
講演２． 「原子レベルの汎用 AI シミュレータ「MATLANTIS」と ENEOS における Materials
Informatics 事例のご紹介」
小島 隆嗣（ENEOS 株式会社）
講演３．「機械学習が導くタンパク質の加速進化」
梅津 光央（東北大学）
講演４．「脂質膜・脂質ナノ粒子の分子シミュレーション －細胞内動態の解析に向けて－」
篠田 渉（岡山大学）
講演５．「QM/MM 法の基礎と創薬応用への新展開」
八木 清・杉田 有治（理化学研究所）
参加登録および申込締切日： 6 月 3 日（月）までに、フォーラムホームページから事前参加登録をお願いいたします。参加人数が上限に達しましたら、参加登録を打ち切らせていただきます
ので、予めご了承ください。
参加費： 一般会員・関連学会員 1000 円、非会員 2000 円、学生 無料
問合先：
構造活性フォーラム 2024 実行委員会
九州大学 加藤 幸一郎（実行委員長）〒819-0395 福岡市西区元岡 744
Tel: 092-802-2922
E-mail: kato.koichiro.957@m.kyushu-u.ac.jp- 31 -

SAR News No.46 (Apr. 2024)///// Activities /////＜会告＞第 52 回構造活性相関シンポジウム
主催： 日本薬学会構造活性相関部会
協賛： 情報計算化学生物学会(CBI 学会)、日本農薬学会、日本薬学会医薬化学部会、日本薬学会レギュラトリーサイエンス部会、理論化学会、日本バイオインフォマティクス学会 (JSBi)
会期： 2024 年 12 月 12 日（木）～13 日（金）会場： 川崎市産業振興会館（神奈川県川崎市幸区堀川町 66 番地 20）1 階ホール・4 階展示場
発表形式：口頭発表・ポスター発表
特別講演（予定・仮題）：
天能精一郎（神戸大学）「特異な相互作用や機能に関わる複雑な電子状態と理論化学的アプローチ」
ミニシンポジウム「レギュラトリーサイエンスと QSAR」（予定・仮題）
基調講演：
本間正充（国立医薬品食品衛生研究所）「レギュラトリーサイエンスと国立医薬品食品衛生研究所」
招待講演：
小山直己（エーザイ株式会社）「ICH M7 ガイドラインに準ずる QSAR を活用した医薬品不純物の変異原性予測評価の基礎と実際」
伊丹悠人（国立環境研究所）「生態リスク評価における QSAR の活用」
他一名
参加登録費・懇親会参加費（予定）
：

参加登録費

懇親会参加費*

金額は不課税

薬学会会員

一般

8,000 円

学生

2,000 円

非会員

一般

10,000 円

学生

3,000 円

金額は税込額

一般

6,000 円

学生

4,000 円

金額は税込額*懇親会会場は川崎市産業振興会館 4 階展示場の予定- 32 -（適用対象外）SAR News No.46 (Apr. 2024)

その他、最新情報はホームページにてご確認ください。
HP：https://www.qsarj.org/52sympo/ （2024 年 4 月以降開設予定）

問い合わせ先：第 52 回構造活性相関シンポジウム実行委員会
国立医薬品食品衛生研究所 古濱彩子（実行委員長）
E-mail: sar2024@qsarj.org- 33 -

SAR News No.46 (Apr. 2024)

部会役員人事
2024 年度 常任世話人
部会長
副部会長
副部会長
会計幹事
庶務幹事
広報幹事
SAR News 編集長
ホームページ委員長

竹田–志鷹 真由子（北里大学 薬学部）
田上 宇乃（味の素（株））
前田 美紀（農業・食品産業技術総合研究機構）
川下 理日人（近畿大学 理工学部）
河合 健太郎（摂南大学 薬学部）
加藤 博明（広島商船高等専門学校）
合田 浩明（昭和大学 薬学部）
高木 達也（大阪大学大学院 薬学研究科）

構造活性相関部会の沿革と趣旨
1970年代の前半、医農薬を含む生理活性物質の活性発現の分子機構、立体構造・電子構造の計算や活性データ処理に対するコンピュータの活用など、関連分野のめざましい発展にともなって、構造活性相関と分子設計に対する新しい方法論が世界的に台頭してきた。このような情勢に呼応するとともに、研究者の交流と情報交換、研究発表と方法論の普及の場を提供することを目的に設立されたのが本部会の前身の構造活性相関懇話会である。1975年5月京都において第1回の「懇話会」（シンポジウム）が旗揚げされ、1980
年からは年1回の「構造活性相関シンポジウム」が関係諸学会の共催の下で定期的に開催されるようになった。
1993年より同シンポジウムは日本薬学会医薬化学部会の主催の下、関係学会の共催を得て行なわれることとなった。構造活性相関懇話会は1995年にその名称を同研究会に改め、シンポジウム開催の実務担当グ
ループとしての役割を果すこととなった。2002年4月からは、日本薬学会の傘下組織の構造活性相関部会として再出発し、関連諸学会と密接な連携を保ちつつ、生理活性物質の構造活性相関に関する学術・研究の振興と推進に向けて活動している。現在それぞれ年1回のシンポジウムとフォーラムを開催するとともに、部会誌のSAR Newsを年2回発行し、関係領域の最新の情勢に関する啓蒙と広報活動を行っている。本部会の沿革と趣旨および最新の動向などの詳細に関してはホームページを参照頂きたい。
(https://sar.pharm.or.jp/)

編集後記
今号では「エピジェネティクス」をテーマとして、2 名の先生にご執筆をお願いしました。理化学研究所の梅原崇史先生には、ヒストンのアセチル化／脱アセチル化に関与する制御分子およびクロマチンの制御研究について、詳細にご紹介いただきました。横浜市立大学の仙石徹先生には、NSD2 をターゲットとした生化学実験、構造解析、MD シミュレーションを組み合わせた包括的な研究をご紹介いただきました。
エピジェネティクスの複雑さに対して、どのような手法やアイディアを元に研究を進められているかについて、本号がみなさまの理解の一端になれば幸いです。ご多忙の中、快くご執筆していただいた先生方に深く感謝申し上げます。
Activities として、昨年 11 月に開催された第 51 回構造活性相関シンポジウムの開催報告とともに、今年
6 月、12 月にそれぞれ開催予定の構造活性フォーラム 2024、第 52 回構造活性相関シンポジウムの会告がございます。みなさまのご参加をお待ち申し上げております。
（編集委員会）SAR News No.46 2024 年 4 月 1 日
発行:日本薬学会 構造活性相関部会長

竹田–志鷹 真由子
SAR News 編集委員会（委員長）合田 浩明、浴本 亨、遠藤 智史、仲西 功、原田 俊幸、幸 瞳*本誌の全ての記事、図表等の無断複写・転載を禁じます。

- 34 -