RDKitでRグループを抜き出す。

化合物を色々と合成していると、中心構造が同じならそこにつく置換機の情報だけを抽出したい場合があります。
たとえば縦横マトリックスを書いて考えたいなんて場合とか。。。

from rdkit import Chem
mols=Chem.SDMolSupplier("hogehoge.sdf")
core = Chem.MolFromSmiles("c1ccccc1)
for m in mols:
    r=Chem.ReplaceCore(m, core)
    print Chem.MolToSmiles( r )

とすると
コア構造が抜かれて各側鎖がピリオドで区切られたsmilesが得られます。
split関数等で分割してあれこれ解析に使えると思います。

広告

コメントを残す

以下に詳細を記入するか、アイコンをクリックしてログインしてください。

WordPress.com ロゴ

WordPress.com アカウントを使ってコメントしています。 ログアウト / 変更 )

Twitter 画像

Twitter アカウントを使ってコメントしています。 ログアウト / 変更 )

Facebook の写真

Facebook アカウントを使ってコメントしています。 ログアウト / 変更 )

Google+ フォト

Google+ アカウントを使ってコメントしています。 ログアウト / 変更 )

%s と連携中