- This topic has 7 replies, 4 voices, and was last updated 11 years, 5 months ago by Anonymous.
June 11, 2012 at 9:58 pm #1303Anonymous
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 #7231Anonymous
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 #7232Anonymous
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 #7277Anonymous
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.html
Does 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 #7234Anonymous
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 #7235Anonymous
Try using score_jd2 with -renumber_pdb first.
June 20, 2012 at 6:37 pm #7278Anonymous
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:
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
(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 #7295Anonymous
BOUNDED in particular requires rswitch and tag, in contradiction to earlier documentation. It will silently eat through the next line of constraint code otherwise.
- You must be logged in to reply to this topic.