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.

Advertisements

Leave a Reply

Fill in your details below or click an icon to log in:

WordPress.com Logo

You are commenting using your WordPress.com account. Log Out / Change )

Twitter picture

You are commenting using your Twitter account. Log Out / Change )

Facebook photo

You are commenting using your Facebook account. Log Out / Change )

Google+ photo

You are commenting using your Google+ account. Log Out / Change )

Connecting to %s