Rosetta
|
Functions | |
def | check_and_print_pdb (count, residue_buffer, residue_letter) |
Variables | |
string | local_pdb_database = "/lab/shared/scaffolds/" |
set your scaffold source here ################################# More... | |
bool | keep_ligand = True |
bool | shit_stat_insres = False |
bool | shit_stat_altpos = False |
bool | shit_stat_modres = False |
dont count MSEs as modiied residues (cos they're so common and get_pdb deal with them previosuly) More... | |
bool | shit_stat_misdns = False |
Print out the residue when the new line is for a new residue. More... | |
list | hetatm_list = [] |
bool | CA_MODEL = True |
bool | NO_OCCUPANCY = True |
string | fastaseq = "" |
string | pdbfile = "" |
pdbname = argv[1] | |
pdbcode = argv[1] | |
chainid = argv[2] | |
outfile = pdbname | |
int | nopdbout = 0 |
int | removechain = 0 |
int | ignorechain = 0 |
string | netpdbname = local_pdb_database + pdbname[1:3] + '/' + pdbname[0:4] + '/' + argv[1] + '_0.pdb' |
lines = popen( 'zcat '+netpdbname,'r').readlines() | |
string | oldresnum = ' ' |
residue_invalid = True More... | |
int | count = 1; |
string | modifiedres = '' |
list | residue_buffer = [] |
string | residue_letter = '' |
bool | residue_invalid = False |
line = lines[i] | |
line_edit = line | |
Is it a modified residue ? More... | |
bool | ok = False |
resnum = line_edit[22:27] | |
if line_edit[13:14]=='P': #Nucleic acid? Skip. More... | |
insres = line[26] | |
altpos = line[16] | |
longname = line_edit[17:20] | |
string | flag_altpos = "---" |
string | flag_insres = "---" |
string | flag_modres = "---" |
string | flag_misdns = "---" |
nres = len(fastaseq) | |
string | flag_successful = "OK" |
outid = open( outfile, 'w') | |
fastaid = stderr | |
def clean_pdb_keep_ligand.check_and_print_pdb | ( | count, | |
residue_buffer, | |||
residue_letter | |||
) |
References run_backbone.float, and ObjexxFCL.len().
clean_pdb_keep_ligand.altpos = line[16] |
bool clean_pdb_keep_ligand.CA_MODEL = True |
string clean_pdb_keep_ligand.chainid = argv[2] |
Referenced by MPInterfaceStatistics.get_size().
int clean_pdb_keep_ligand.count = 1; |
clean_pdb_keep_ligand.fastaid = stderr |
string clean_pdb_keep_ligand.fastaseq = "" |
Referenced by main().
clean_pdb_keep_ligand.flag_altpos = "---" |
clean_pdb_keep_ligand.flag_insres = "---" |
clean_pdb_keep_ligand.flag_misdns = "---" |
clean_pdb_keep_ligand.flag_modres = "---" |
string clean_pdb_keep_ligand.flag_successful = "OK" |
list clean_pdb_keep_ligand.hetatm_list = [] |
int clean_pdb_keep_ligand.ignorechain = 0 |
clean_pdb_keep_ligand.insres = line[26] |
bool clean_pdb_keep_ligand.keep_ligand = True |
string clean_pdb_keep_ligand.line_edit = line |
Is it a modified residue ?
What does this do ?
other substitution (of atoms mainly)
if so replace it with its canonical equivalent !
clean_pdb_keep_ligand.lines = popen( 'zcat '+netpdbname,'r').readlines() |
string clean_pdb_keep_ligand.local_pdb_database = "/lab/shared/scaffolds/" |
set your scaffold source here #################################
clean_pdb_keep_ligand.longname = line_edit[17:20] |
string clean_pdb_keep_ligand.modifiedres = '' |
clean_pdb_keep_ligand.netpdbname = local_pdb_database + pdbname[1:3] + '/' + pdbname[0:4] + '/' + argv[1] + '_0.pdb' |
bool clean_pdb_keep_ligand.NO_OCCUPANCY = True |
int clean_pdb_keep_ligand.nopdbout = 0 |
clean_pdb_keep_ligand.nres = len(fastaseq) |
bool clean_pdb_keep_ligand.ok = False |
clean_pdb_keep_ligand.oldresnum = ' ' |
residue_invalid = True
string clean_pdb_keep_ligand.outfile = pdbname |
string clean_pdb_keep_ligand.pdbfile = "" |
Referenced by densityTools(), main(), ReadCAsFromPDB(), relax_test(), run(), and spanfile_from_pdb().
clean_pdb_keep_ligand.pdbname = argv[1] |
int clean_pdb_keep_ligand.removechain = 0 |
list clean_pdb_keep_ligand.residue_buffer = [] |
bool clean_pdb_keep_ligand.residue_invalid = False |
string clean_pdb_keep_ligand.residue_letter = '' |
clean_pdb_keep_ligand.resnum = line_edit[22:27] |
if line_edit[13:14]=='P': #Nucleic acid? Skip.
resnum = line_edit[23:26] oldresnum = resnum while (resnum == oldresnum): print "HERE" i += 1 line = lines[i] resnum = line_edit[23:26]
clean_pdb_keep_ligand.shit_stat_altpos = False |
clean_pdb_keep_ligand.shit_stat_insres = False |
bool clean_pdb_keep_ligand.shit_stat_misdns = False |
Print out the residue when the new line is for a new residue.
if unsuccessful
is there a residue in the buffer ? If so print it out, then reinitialize
print "call from the loop:\n" ## new
if unsuccessful
bool clean_pdb_keep_ligand.shit_stat_modres = False |
dont count MSEs as modiied residues (cos they're so common and get_pdb deal with them previosuly)