pharmacophore alignment using RDKit

One of the my favorite paper is “SAR NEWS 19”.
link here
I wrote some script for p4core align and uploaded to my repository .
get_p4core.py makes pharmacophore object from input file and dump it.
phcoresearchmodule.py aligns input sdf to pharmacophore from get_p4core.py. But I must improve this script ASAP.

Basic usage…

To Get Pharmacophore.

iwatobipen$ python get_p4core.py  query_mol.sdf

You could make query_mol.sdf from your interested PDB file or any .
And after run the script, you could get pharmacophore object named query_mol_pcophore.pkl.

To Get Align molecules to the p4core.

iwatobipen$ python phcoresearchmodule.py -i mols.sdf -p query_mol_pcophore.pkl -o output.sdf

-o is option. default output file name is res.sdf.

At first this script reads input sdf and pharmacophore object.
Then it tries to generate 3D conf that satisfies the pharmacophore and optimize conf using universal force field.
Finally all generated conformers are aligned and write to sdf.
I tested these scripts, but does not work well. 😦
If anyone who has any advices, please post comment or advice for me.

広告

コメントを残す

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

WordPress.com ロゴ

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

Twitter 画像

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

Facebook の写真

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

Google+ フォト

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

%s と連携中