diff PLIDflow/scripts/gpffilemaker.R @ 6:795e11fac81b draft default tip

Included new tools for standardization
author bitlab
date Wed, 22 Apr 2020 06:12:00 -0400
parents 6fcfa4756040
children
line wrap: on
line diff
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/PLIDflow/scripts/gpffilemaker.R	Wed Apr 22 06:12:00 2020 -0400
@@ -0,0 +1,116 @@
+#gpffilemaker makes the GPF file
+
+#!/usr/bin/env Rscript      
+args = commandArgs(trailingOnly=TRUE) 
+
+if(length(args) < 1){
+  stop("USE: Rscript gpffilemaker.R <receptor_name> <session_dir>")
+}
+
+#Read file containing PDB geometric center and number of points x, y ,z  
+#archivo_center_npts <- scan("/home/galaxy/galaxy/tools/proteindocking/scripts/pdbcenter_npts.txt", what = character(), quiet = TRUE)
+archivo_center_npts <- scan("pdbcenter_npts.txt", what = character(), quiet = TRUE)
+
+#Read file containing receptor atom types
+ver <- scan("receptor_atm_types.txt", what = character(), quiet = TRUE, na.strings = "NULLAJSKAJSFL")
+
+
+#Arguments definition 
+receptor_name <- args[1]
+session_dir <- args[2]
+aux_receptor_name <- basename(args[1])
+#aux_receptor_name <- paste("/home/galaxy/galaxy/tools/proteindocking/scripts/", basename(args[1]), sep="")
+npts_x <- archivo_center_npts[4]
+npts_y <- archivo_center_npts[5]
+npts_z <- archivo_center_npts[6] 
+gridcenter_x <- archivo_center_npts[1]
+gridcenter_y <- archivo_center_npts[2]
+gridcenter_z <- archivo_center_npts[3]
+
+#Creation of the GPF file
+#Row 1 
+fila_1_e1 <- paste("npts", npts_x, sep = " ")
+fila_1_e2 <- paste(fila_1_e1, npts_y, sep = " ")
+fila_1_e3 <- paste(fila_1_e2, npts_z, sep = " ")
+fila_1 <- paste(fila_1_e3, "# num.grid points in xyz", sep = " ")
+
+#Row 2 
+fila_2_e2 <- paste(aux_receptor_name, ".maps.fld", sep = "")
+fila_2_e1 <- paste("gridfld", fila_2_e2, sep = " " )
+fila_2 <- paste(fila_2_e1, "# grid_data_file", sep = " ")
+
+#Row 3 
+fila_3 <- paste("spacing 1.0", "# spacing(A)", sep = " ")
+
+#Row 4 
+#fila_4 <- paste("receptor_types A C HD N NA OA S SA", "# receptor atom types", sep = " ")
+cabecera_f4 <-  ver[1]
+vector_f4 <- c()
+for(f4 in 2:length(ver)){
+  cabecera_f4 <- paste(cabecera_f4, ver[f4], sep = " ")
+}
+
+fila_4 <- paste("receptor_types", cabecera_f4, "# receptor atom types", sep = " ")
+#print(fila_4)
+
+#Row 5  
+fila_5 <- paste("ligand_types C HD OA", "# ligand atom types", sep = " ")
+
+#Row 6 
+fila_6_e2 <- paste(aux_receptor_name, ".pdbqt", sep = "")
+fila_6_e1 <- paste("receptor", fila_6_e2, sep = " ")
+fila_6 <- paste(fila_6_e1, "# macromolecule", sep = " ")
+
+#Row 7 
+fila_7_e1 <- paste("gridcenter",gridcenter_x, sep = " ")
+fila_7_e2 <- paste(fila_7_e1, gridcenter_y, sep = " ")
+fila_7_e3 <- paste(fila_7_e2, gridcenter_z, sep = " ")
+fila_7 <- paste(fila_7_e3, "# xyz-coordinates or auto", sep = " ")
+
+#Row 8  
+fila_8 <- paste("smooth 0.5", "# store minimum energy w/in rad(A)", sep = " ")
+
+#Row 9 
+fila_9_e2 <- paste(aux_receptor_name, ".C.map", sep = "")
+fila_9_e1 <- paste("map", fila_9_e2, sep = " ")
+fila_9 <- paste(fila_9_e1, "# atom-specific affinity map", sep = " ")
+
+#Row 10 
+fila_10_e2 <- paste(aux_receptor_name, ".HD.map", sep = "")
+fila_10_e1 <- paste("map",fila_10_e2, sep = " ")
+fila_10 <- paste(fila_10_e1, "# atom-specific affinity map", sep = " ")
+
+#Row 11 
+fila_11_e2 <- paste(aux_receptor_name, ".OA.map", sep = "")
+fila_11_e1 <- paste("map", fila_11_e2, sep = " ")
+fila_11 <- paste(fila_11_e1, "# atom-specific affinity map", sep = " ")
+
+#Row 12 
+fila_12_e2 <- paste(aux_receptor_name, ".e.map", sep = "")
+fila_12_e1 <- paste("elecmap", fila_12_e2, sep = " ")
+fila_12 <- paste(fila_12_e1, "# electrostatic potential map", sep = " ")
+
+#Row 13
+fila_13_e2 <- paste(aux_receptor_name, ".d.map", sep = "")
+fila_13_e1 <- paste("dsolvmap", fila_13_e2, sep = " ")
+fila_13 <- paste(fila_13_e1, "# desolvation potential map", sep = " ")
+
+#Row 14
+fila_14 <- paste("dielectric -0.1465", "# <0, AD4 distance-dep.diel;>0, constant", sep = " ")
+
+#Paste all rows with a line jump for separation 
+filas <- c(fila_1, fila_2, fila_3, fila_4, fila_5, fila_6, fila_7, fila_8, fila_9, fila_10, fila_11, fila_12, fila_13, fila_14)
+
+datos<- c()
+for(i in 1:14){
+datos <- paste(c(datos, filas[i]))
+}
+
+#nombre_final <- paste("/home/galaxy/galaxy/tools/proteindocking/scripts/", basename(args[1]), ".gpf", sep="")
+nombre_final <- paste(basename(args[1]), ".gpf", sep="")
+
+write(datos, file=nombre_final, append = FALSE)
+
+#write(basename(args[1]), file= paste("/home/galaxy/galaxy/tools/proteindocking/scripts/", basename(args[1]), ".txt", sep="")) #parte añadida 18jun2018 
+write(basename(args[1]), file= paste(basename(args[1]), ".txt", sep="")) #parte añadida 18jun2018 
+