RDKitでQED

QEDはNature Chemistry に報告されている 指標です。
“Quantitative Estimate of Drug-Likeness” (QED)
QEDとケミストの得点を比較したデータ等、興味深いデータが乗っています。
SIにpipeline pilotのプロトコールがあるのでppがあればさくっとできますが、同様のことをRDKitでも実現できます。
silicos-it のHPに紹介されているBiscu-itを入れるとよいです。インストールはHPにある通りで

$ cd /usr/local/src
$ mv ~/Downloads/biscu-it-*.tar.gz .
$ sudo tar -xvf biscu-it-*.tar.gz
$ cd biscu-it-*



$ python setup.py install

以上

from silicos_it.descriptors import qed
from rdkit import Chem
mols = [m for m in Chem.SDMolSupplier("./hoge.sdf")]
for m in mols:
    print qed.default(m)

とすると読み込んだSDFのQEDが算出されます。

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