-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathtest.py
30 lines (22 loc) · 809 Bytes
/
test.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
import os
# from LYRUS.lyrusClass import lyrusClass, lyrusPredict
from LYRUS.lyrusClass import lyrusClass, lyrusPredict
gene = 'PTEN'
uniprot = 'P60484'
currDir = os.getcwd()
outputDir = '{}/test1'.format(currDir)
try:
os.mkdir(outputDir)
except:
print('Output directory already exist')
#load model
lyrusModel = lyrusClass(gene, uniprot, outputDir, savFile='{}/{}/input.txt'.format(outputDir,gene))
#download orthologs
lyrusModel.getFasta()
#download PDB
lyrusModel.getPDB()
#calculate all the parameters except for fathmm
lyrusModel.getParameters(maestroDir='MAESTRO_OSX_x64',p2rankDir='p2rank_2.2', EVparam=5)
#run fathmm at http://fathmm.biocompute.org.uk/inherited.html and save results in fathmm file
# fathmmFile = 'test1/fathmm.txt'
# lyrusPredict(gene, fathmmFile, outputDir, uniprot)