autodataming的个人博客分享 http://blog.sciencenet.cn/u/autodataming

博文

中药化合物库的下载

已有 6814 次阅读 2013-11-14 21:37 |个人分类:DrugDesign|系统分类:科研笔记| 中药, 下载, 化合物库, tcmsp

files.txt

背景:中药有效果是有目共睹的,中药中的化合物是繁多的,已经提示了我们可以多靶标治疗。

从中药中寻找活性化合物,是新药发现的一种手段。

tcmsp网站,http://tcmspnw.com/login_clearSession

作为一个开源网站,可以提供化合物的单个下载。

据说今后也会提供下载包,批量下载的。

现在我和以及以前的师姐等不及了,所以就自己写了脚本从网站上不这些分子下载下来。

方法:downtcm.pl

downtcm.pl

需要用到files.txt

++++++++++++++++++++++++++++++++

#!/usr/bin/perl -w
use strict;
use LWP;

use HTTP::Request::Common qw ( POST  GET);

#FUNCTION:  to download all the molecue in the tcmsp website
#author:  chenzhaoqiang 744891290@qq.com
#need input file:files.txt
#version1


my $browser=LWP::UserAgent->new(cookie_jar=>{});


$browser->requests_redirectable;

$browser->agent("MyApp/0.1");

my $username='tcmspuser';
my $password='tcmspuser';

my $req=POST('http://tcmspnw.com/login_login',['user.userComments'=>$username,'user.userPassword'=>$password]);

my $res=$browser->request($req);

$req=GET('http://tcmspnw.com/download?moleculeDownload=zhx_yzhh_Molecule219');
$res=$browser->request($req);
 if ($res->is_success) {
     print $res->content;
 }
 else {
     print $res->status_line, "\n";
 }
 
 #test success
 

#TO write a loop,to batch download chemical molecule
my $baseurl='http://tcmspnw.com/download?moleculeDownload=';
my @ids=(0..999);
my $filebasename;
my $filename;
open FH,"f:/chenzhaoqiang/files.txt";
while(<FH>)
{
chomp;
$filebasename=$_;
push @ids,"";
foreach my $id(@ids)
 {
#print "$id\n";
$filename=$filebasename.$id.'.mol2';
my $totalurl=$baseurl.$filebasename.'_Molecule'.$id;
print $totalurl,"\n";
$req=GET($totalurl);
$res=$browser->request($req);
    if ($res->content!~/Sorry\,File not found /)
    {
    open FF,">F:/downtcm/$filename";
     print FF $res->content;
     }
    else
    {
     print $res->status_line, "\n";
     }
 }
}
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
=pod
#http://tcmspnw.com/login_login
#open FH,"F:/files.txt";

my $baseurl='http://tcmspnw.com/download?moleculeDownload=';
my $fileclass='zhx_yzhh_Molecule';
my $fileid='219';

my @ids=(1..999);

foreach my $id(@ids)
{
#print "$id\n";
my $totalurl=$baseurl.$fileclass.$id;
print $totalurl,"\n";
#  my $content=getstore($totalurl,$id);
#  print $content,"n";
}

=cut

 



https://wap.sciencenet.cn/blog-950202-741780.html

上一篇:sybyl中小分子的能量优化
下一篇:如何下载cas已经知道的化合物的结构
收藏 IP: 210.73.50.*| 热度|

0

该博文允许注册用户评论 请点击登录 评论 (0 个评论)

数据加载中...

Archiver|手机版|科学网 ( 京ICP备07017567号-12 )

GMT+8, 2024-4-18 17:37

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部