Entering edit mode
4.0 years ago
tpm
▴
30
Hello. I ran this code below
from Bio.PDB import *
parser = PDBParser()
io = PDBIO()
structure = parser.get_structure('3F9R', '3F9R.pdb')
model = structure[0]
dssp = DSSP(model, '3F9R.pdb')
residues = list(dssp)
residues
How do I select columns for calculated RSA and residue from this dictionary?
If I want to see the contents in this for residues, I have an output file that looks in short like this.
[(1,
'G',
'-',
0.7857142857142857,
360.0,
-167.6,
0,
0.0,
2,
-0.1,
0,
0.0,
444,
-0.1),
(2,
'M',
'-',
0.14893617021276595,
-67.2,
140.1,
1,
-0.1,
36,
-0.2,
2,
-0.0,
35,
-0.1),
(3,
'K',
'-',
0.23902439024390243,
-78.0,
-14.6,
34,
-2.9,
35,
-0.1,
2,
-0.2,
197,
-0.1),
Have you checked the manual? https://biopython.org/docs/1.75/api/Bio.PDB.DSSP.html
To be honest Fatima, I did check, I genuinely think that I don't know how I can extract RSA and its corresponding residue. The actual code is what I am kindly asking in case someone may have passed through this phase. Maybe its because I am not a good programmer, and it can be quite confusing especially when I can`t bring it to perspective.
Maybe you can try a simpler database, for example:
https://swift.cmbi.umcn.nl/gv/dssp/DSSP_3.html
If you prefer to use
biopython
, then try this to get the relative ASA:You may need to read more about python dictionaries. Also, I think it's a good idea to take a look at similar questions here. You can see a list of these questions on the right panel.
Thank you, so this code computes relative ASA directly. I will read the link you provided through.