Commit 82c1835f by Chloe Dequeker

modif bug analyse complex

parent 7b5159b1
.*.swo .*.swo
.*.swp .*.swp
.*.swn
src/.*.swo src/.*.swo
src/.*.swp src/.*.swp
src/.*.swn
INTBuilder
...@@ -511,7 +511,9 @@ int main(int argc, char** argv){ ...@@ -511,7 +511,9 @@ int main(int argc, char** argv){
fflush(stdout); fflush(stdout);
/* Get all the docking informations regarding our complex */ /* Get all the docking informations regarding our complex */
dock_res = getDataForComplex(); if(!complexPDB){
dock_res = getDataForComplex();
}
/* Get the PDB structure for the receptor and the ligand */ /* Get the PDB structure for the receptor and the ligand */
pdbR = readPDB(receptor); pdbR = readPDB(receptor);
...@@ -522,8 +524,19 @@ int main(int argc, char** argv){ ...@@ -522,8 +524,19 @@ int main(int argc, char** argv){
*/ */
getAnglesFromCA1andCA5_sophie(pdbL,&alpha,&beta,&gamma); getAnglesFromCA1andCA5_sophie(pdbL,&alpha,&beta,&gamma);
if(complexPDB){ /* If we are building for a given complex, no rotation */
/* Write the number of the conformation and get the interface */
fprintf(outputFile_lig,"%d ",i);
fprintf(outputFile_rec,"%d ",i);
getInterface(pdbR,pdbL);
fprintf(outputFile_lig,"\n");
fprintf(outputFile_rec,"\n");
if(constructPDB){
writePDB(newPDB, ligand,target_conf);
}
if(target_conf > 0){ /* If we are building a specific conformation */ }else if(target_conf > 0){ /* If we are building a specific conformation */
dalpha = dock_res->alpha[target_conf] - alpha; dalpha = dock_res->alpha[target_conf] - alpha;
dbeta = dock_res->beta[target_conf] - beta; dbeta = dock_res->beta[target_conf] - beta;
dgamma = dock_res->gamma[target_conf] - gamma; dgamma = dock_res->gamma[target_conf] - gamma;
...@@ -544,18 +557,6 @@ int main(int argc, char** argv){ ...@@ -544,18 +557,6 @@ int main(int argc, char** argv){
} }
freePDB(newPDB); freePDB(newPDB);
newPDB = NULL; newPDB = NULL;
}else if(complexPDB){ /* If we are building for a given complex, no rotation */
/* Write the number of the conformation and get the interface */
fprintf(outputFile_lig,"%d ",i);
fprintf(outputFile_rec,"%d ",i);
getInterface(pdbR,pdbL);
fprintf(outputFile_lig,"\n");
fprintf(outputFile_rec,"\n");
if(constructPDB){
writePDB(newPDB, ligand,target_conf);
}
}else{ /* For all conformations */ }else{ /* For all conformations */
for(i=1;i<dock_res->nbConf;i++){ for(i=1;i<dock_res->nbConf;i++){
...@@ -592,7 +593,8 @@ int main(int argc, char** argv){ ...@@ -592,7 +593,8 @@ int main(int argc, char** argv){
freePDB(pdbL); freePDB(pdbL);
freePDB(pdbR); freePDB(pdbR);
freeDock(dock_res); if(!complexPDB)
freeDock(dock_res);
free_argLine(); free_argLine();
printf(" done\n"); printf(" done\n");
......
...@@ -142,13 +142,16 @@ void free_argLine(){ ...@@ -142,13 +142,16 @@ void free_argLine(){
* start of the program * start of the program
*/ */
free(outputDir); free(outputDir);
free(dockingFile);
free(pdbDir); free(pdbDir);
free(receptor); free(receptor);
free(ligand); free(ligand);
fclose(outputFile_rec); fclose(outputFile_rec);
fclose(outputFile_lig); fclose(outputFile_lig);
if(!complexPDB){
free(dockingFile);
}
if(constructPDB){ if(constructPDB){
free(outputPDB); free(outputPDB);
} }
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment