Member Site › Forums › PyRosetta › PyRosetta – General › Fragments file
- This topic has 4 replies, 3 voices, and was last updated 10 years, 9 months ago by
Anonymous.
-
AuthorPosts
-
-
November 26, 2013 at 10:40 am #1766
Anonymous
The fragments file:aat000_03_05.200_v1_3
position: 1 neighbors: 2001pg4 A 309 W L -108.068 162.569 -178.328 0.102 4.369 116.235 3 0.000 P 1 F 1
1pg4 is the name of protein.But what is the “A 309 W L -108.068 162.569 -178.328 0.102 4.369 116.235 3 0.000 P 1 F 1”?What’s the meaning of those?
“{series}{pdb}{chain}{size}_{strategy}.{depth}_{version}” is mentioned in Rosetta user guide.But it is different.When I used pose=pose_from_pdb(“xx.pdb”),there is a warning(core.io.pdb.file_data: [ WARNING ] discarding 2 atoms at position 3351 in file xx.pdb. Best match rsd_type: ASP).Why?
-
November 28, 2013 at 6:52 pm #9531
Anonymous
What you found in the guide refers just to the fragment file name. I am looking into the fragments as well now, here’s what I found in the code.
Method ConstantLengthFragSet::read_fragment_stream() from srccorefragmentConstantLengthFragSet.cc handles reading of the fragment file. The string it expects consists of fixed length fields:
{pdbid} {chain} {aa_index} {aa} {ss} {phi} {psi} {omega} {x} {y} {z}
where
{aa_index} – position of the residue in the original file.
{aa} – one-letter id of the residue
{ss} – secondary structure, H | L | E (Helix | Loop | Extended (beta-sheet))
{phi} {psi} {omega} – dihedral angles of the residue
{x} {y} {z} – coordinates of the Ca atoms (optional)The rest of the string seems to be never parsed.
Can’t say anything about the discarded atoms issue.
-
November 29, 2013 at 5:56 pm #9533
Anonymous
Regarding the discarded atoms issue, it’s hard to say without seeing the PDB file in question, but that message is usually because there’s some naming issue with the residue. Look at that amino acid – does it conform to the standard atom naming conventions? One issue may be if the structure doesn’t have a regular aspartate, but instead has some sort of posttranslational modification to the residue. If Rosetta isn’t appropriately told about the modification (through patch files), it will simply ignore the extra atoms, and do the best it can with what is left.
-
November 30, 2013 at 8:10 am #9534
Anonymous
Do you know the meaning of the rest string(3 0.000 P 1 F 1)?
-
December 2, 2013 at 3:28 pm #9536
Anonymous
Sorry, no. Looking at the code, it looks like those numbers aren’t even printed anymore, for more recent versions of the fragment picker. (Though I could be misreading things.)
From what little I can tell, it looks like the P X F Y designation is for “position X fragment Y”. The 3 and the 0.000 are still a mystery, though.
-
-
AuthorPosts
- You must be logged in to reply to this topic.