Member Site › Forums › Rosetta 3 › Rosetta 3 – General › about Residue Number in Rosetta
- This topic has 7 replies, 4 voices, and was last updated 12 years, 2 months ago by
Anonymous.
-
AuthorPosts
-
-
June 11, 2012 at 9:58 pm #1303
Anonymous
Does Rosetta use the residue number written in the pdb files? Or it just take the first shown residue has the residue number 1?
I need to clarify this issue so that the constrains are applied correctly.Thank you very much.
-
June 12, 2012 at 12:23 am #7231
Anonymous
The answer to the general question is “it depends on where the numbering is being used”.
Generally, whenever a chain letter is specified along with the residue number, the numbering corresponds to the PDB numbering. Whenever a residue number is specified as just a single number (without an alphabetic chain designator) the numbering is “pose numbering” – that is, each residue numbered sequentially starting at 1. (Be careful, as occasionally residues in the input PDB aren’t loaded into the pose, and as such get skipped in pose numbering.)
For the standard constraint file format (http://www.rosettacommons.org/manuals/archive/rosetta3.4_user_guide/de/d50/constraint_file.html), the numbering is pose (sequential starting from 1) numbering, consistent with the fact that there isn’t an opportunity to specify the chain.
(Note that there’s been discussion about supporting PDB numbering in constraint files, but I don’t believe that’s been developed enough to be in/be usable with any of the released versions.)
-
June 12, 2012 at 1:00 am #7232
Anonymous
Generally speaking:
Resfiles: PDB numbering (includes chain; supports insertion codes)
Loops files: “pose numbering” as Rocco designated
Constraints files: pose numbering
Anything else: Pose numbering unless it also requests a chain, in which case it’s probably PDB numbering.
I have a Q&D app to strip PDB info and turn it into a Rosetta-fied PDB called “ForgetPDBInfo” if you need it. I think there may also be an option you can pass to the score or score_jd2 app that does the same thing.
Thanks for asking, by the way – we get a LOT of questions here that are caused by people using PDB numbering where they should have used pose numbering.
-
June 20, 2012 at 2:44 pm #7277
Anonymous
Hey there,
wondering about the same problem I found that at least for SiteConstraint a chain ID and a residue number are given in the Docking manual of 3.4 under Tips: http://www.rosettacommons.org/manuals/archive/rosetta3.4_user_guide/d0/de4/docking_protocol.htmlDoes that mean this residue number is the one in the pdb file then in the respective chain?
Also for the other constraint types residue numbers with chain identifiers are given in the RosettaCon2011Demo “dock_with_user_defined_constraints”…
Any clarification would be highly appreciated
Additionally, when trying to use AmbiguousConstraint Constraint1 … Constraint2 …
only Constraint1 is mentioned in the output/logfile. AmbiguousConstraint is not mentioned in the output/logfile.
Is that normal or indicative of a mistaken I may have made? -
June 12, 2012 at 3:23 pm #7234
Anonymous
Really appreciate your answers. I think I need this app, where can I find it?
Thank you very much!
-
June 12, 2012 at 3:45 pm #7235
Anonymous
Try using score_jd2 with -renumber_pdb first.
-
June 20, 2012 at 6:37 pm #7278
Anonymous
For constraints, you should be able to specify residues either way (at least for recent Rosetta releases). If you specify it with a postfixed chain letter, it uses PDB numbering (the number of the residue in the PDB file), whereas if you specify it as a bare number without a chain, you would use pose numbering (sequential numbers starting at 1). (Note that you currently can’t use PDB numbering with residues with insertion codes – you would have to use pose numbering for them.)
The constraint input code is sometimes rather brittle. I think the issue you’re running into with the AmbiguousConstraint is that the AtomPairConstraint (as well as SiteConstraint and probably a number of others) silently eats the rest of the line.
The “approved” way of dealing with AmbiguousConstraints and other types of multiconstraints is to put each subconstraint on its own line:AmbiguousConstraint
CoordinateConstraint OD1 6 CA 6 95.723 49.348 9.574 BOUNDED 0 0.1 0.5 0.5 tag
CoordinateConstraint ND2 6 CA 6 95.723 49.348 9.574 BOUNDED 0 0.1 0.5 0.5 tag
END_AMBIGUOUS(Note that as the residue number is not followed by a chain, it’s being given with pose numbering.)
-
June 24, 2012 at 8:53 pm #7295
Anonymous
BOUNDED in particular requires rswitch and tag, in contradiction to earlier documentation. It will silently eat through the next line of constraint code otherwise.
-
-
AuthorPosts
- You must be logged in to reply to this topic.