00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031
00032
00033 #ifndef DH_MDPOCKET
00034 #define DH_MDPOCKET
00035
00036
00037
00038 #include <math.h>
00039 #include <time.h>
00040 #include <stdio.h>
00041 #include <stdlib.h>
00042 #include <string.h>
00043 #include <ctype.h>
00044
00045 #include "fpocket.h"
00046 #include "fpout.h"
00047 #include "writepocket.h"
00048 #include "tpocket.h"
00049 #include "dparams.h"
00050 #include "descriptors.h"
00051 #include "neighbor.h"
00052 #include "pocket.h"
00053 #include "cluster.h"
00054 #include "refine.h"
00055 #include "aa.h"
00056 #include "utils.h"
00057 #include "calc.h"
00058 #include "mdparams.h"
00059 #include "mdpbase.h"
00060 #include "memhandler.h"
00061 #include "mdpout.h"
00062
00063
00064 #define M_MDP_OUTP_HEADER "snapshot pock_volume pock_asa pock_pol_asa pock_apol_asa pock_asa22 pock_pol_asa22 pock_apol_asa22 nb_AS mean_as_ray mean_as_solv_acc apol_as_prop mean_loc_hyd_dens hydrophobicity_score volume_score polarity_score charge_score prop_polar_atm as_density as_max_dst"
00065 #define M_MDP_OUTP_FORMAT "%d %4.2f %4.2f %4.2f %4.2f %4.2f %4.2f %4.2f %d %4.2f %4.2f %4.2f %4.2f %4.2f %4.2f %d %d %4.2f %4.2f %4.2f"
00066 #define M_MDP_OUTP_VAR(i, d) i, d->volume, \
00067 d->surf_vdw14, \
00068 d->surf_pol_vdw14, \
00069 d->surf_apol_vdw14, \
00070 d->surf_vdw22, \
00071 d->surf_pol_vdw22, \
00072 d->surf_apol_vdw22, \
00073 d->nb_asph,\
00074 d->mean_asph_ray, \
00075 d->masph_sacc, \
00076 d->apolar_asphere_prop, \
00077 d->mean_loc_hyd_dens, \
00078 d->hydrophobicity_score, \
00079 d->volume_score, \
00080 d->polarity_score, \
00081 d->charge_score, \
00082 d->prop_polar_atm, \
00083 d->as_density, \
00084 d->as_max_dst
00085
00086
00087 void mdpocket_detect(s_mdparams *par) ;
00088 c_lst_pockets* mdprocess_pdb(s_pdb *pdb, s_mdparams *mdparams,FILE *pqrout, int snnumber);
00089 void mdpocket_characterize(s_mdparams *par) ;
00090 s_pocket* extract_wanted_vertices(c_lst_pockets *pockets,s_pdb *pdb);
00091 s_pdb *open_pdb_file(char *pdbname);
00092 void write_md_descriptors(FILE *f, s_pocket *p, int i);
00093 int *get_wanted_atom_ids(s_pdb *prot,s_pdb *pocket, int *n);
00094
00095 #endif