From 6e45140c59253fb2f286bfd73f763c18e7ef104b Mon Sep 17 00:00:00 2001
From: Navan Chauhan <navanchauhan@gmail.com>
Date: Mon, 6 Jul 2020 13:03:43 +0530
Subject: added docking score feeature to main.sh and report generation

---
 scripts/get_dock_score.py | 46 ++++++++++++++++++++++++++++++++++++++++++++++
 1 file changed, 46 insertions(+)
 create mode 100644 scripts/get_dock_score.py

(limited to 'scripts/get_dock_score.py')

diff --git a/scripts/get_dock_score.py b/scripts/get_dock_score.py
new file mode 100644
index 0000000..a47312b
--- /dev/null
+++ b/scripts/get_dock_score.py
@@ -0,0 +1,46 @@
+#!/usr/bin/python3
+
+import argparse
+
+parser = argparse.ArgumentParser(description=desc_text)
+parser.add_argument("-p","--protein",help="Path to protein file")
+parser.add_argument("-l","--ligand",help="Path to ligand_out file")
+
+args = parser.parse_args()
+
+if args.protein == None:
+    print("Error: Please specify protein file")
+    exit(1)
+if args.ligand == None:
+    print("Error: Please specify ligand file")
+    exit(1)
+
+
+protein = args.protein
+ligand = args.ligand
+
+from os.path import basename
+
+print("# " + str(basename(protein)).replace(".pdbqt","") + "-" + str(basename(ligand)).replace("_out.pdbqt",""), end="\n\n")
+
+from tabulate import tabulate
+
+file = open(ligand, "r")
+lines = file.readlines()
+results = []
+i = 1
+for line in lines:
+    ta = []
+    if line.find('REMARK VINA') == 0 and line.split()[3] != "":
+        l = line.split()
+        ta.append(i)
+        ta.append(l[3])
+        ta.append(l[4])
+        ta.append(l[5])
+        i += 1
+    if ta != []:
+        results.append(ta)
+
+print("## Docking Scores",end="\n\n")
+print(tabulate(results,headers=["No.","Affinity","rmsd l.b","rmsd u.b"]))
+print("",end="\n\n")
-- 
cgit v1.2.3