|
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)