- This topic has 2 replies, 2 voices, and was last updated 9 years, 1 month ago by Anonymous.
October 12, 2014 at 4:00 pm #2029Anonymous
I am trying to use the loopmodel ccd to build a loop region with missing residue positions in the original PDB file.
The whole sequence is of 228 amino acid and the missing region is 223-228.
Can I ask
1) As the loop to be built is at the end of the sequence (i.e. termini), how can I specify the loop file?
If I use “LOOP 222 228 0 0 0”, I was told (ccd_1.log attached):
“core.kinematics.FoldTree: FoldTree::reorder( 1 ) failed, new/old edge_list_ size mismatch”
If I use “LOOP 222 228 228 0 0”, I was told the same thing (ccd_2.log attached).
It is exciplitly indicated that “quick_ccd can also remodel termini”
So how can I really do this?
2) I am building homology model. I would like to minimize the relax in the rest of the structure (i.e. before residue 223).
I am going to use the options below following “bcorreia” said in
Is this suitable?
# NOT use -loops:extended
-loops:frag_sizes 9 3 1
3) Does “-relax::fastrelax_repeats” exist in the “loopmodel.linuxgccrelease”?
I was told
ERROR: ERROR: Option matching -relax:fastrelax_repeats not found in command line top-level context
when using it.
Thank you very much.
October 13, 2014 at 3:34 pm #10419Anonymous
There is at least one demo on loop length changing. The one I know of is at Rosetta/demos/protocol_capture/2010/AnchoredDesign/loop_length_changing. This demo assumes you have an internal loop, but the technique will probably work for termini.
The key is that Rosetta is not doing the extension. Give Rosetta a PDB with the correct number of residues. Generate that by adding residues in a text editor. Garbage coordinates are ok so long as the backbone heavyatoms (C, CA, N, O) are present. Rosetta has occasional bugs with too many (0,0,0) atoms, so you may want to slightly randomize the coordinates, or just copy coordinates from elsewhere in the protein. Make sure you are using the “extended loop” option in the loop file and it will probably work.
The name of the fastrelax option appears to have been changed to relax:default_repeats.
October 13, 2014 at 5:10 pm #10420Anonymous
Thank you for your help. I will try your methods later. In the meantime, I have been trying to use “floppy tail” as I think the scenario is exactly what I want.
I will ask the question in another post.
- You must be logged in to reply to this topic.