Member Site › Forums › Rosetta 3 › Rosetta 3 – Applications › Loop modeling: KIC with fragments producing unusual models
- This topic has 1 reply, 2 voices, and was last updated 9 years, 2 months ago by
Anonymous.
-
AuthorPosts
-
-
June 1, 2015 at 7:06 pm #2214
Anonymous
Dear Friends,
I have used KIC with fragment protocol for loop modeling as follows
-s test.pdb
-in:file:fullatom
-loops:loop_file test.loops
-loops:frag_files frags/test.200.9mers frags/test.200.3mers none
-loops:frag_sizes 9 3 1
-loops:remodel perturb_kic_with_fragments
-loops:refine refine_kic_with_fragments
-out:nstruct 50
-out:file:silent test.out/test.silent
-out:file:fullatom
-out:file:silent_struct_type binary
-ex1
-ex2
-overwrite
-mute core.io.database
-mute protocols.looprelax.FragmentPerturber
-mute core.fragments.ConstantLengthFragSet
-mute protocols.loop_modeling.loggers.ProgressBarI am using rosetta_2014wk35_bundle.
The loop file without and with extension designated as follows
1. LOOP 147 161 160 0 0
2. LOOP 147 161 160 0 1
My output have been low energy structures with unusual secondary structure and very long bond lengths for the loop modeled region (pdb001_without extension) with a loop file without extension. While with loop extension the secondary structure looks normal but atomic detail has C-O distance of ~66 Å for each residue in the loop model region see pdb002_with_extension.
Could someone help me to figure out what is causing this. I have used the CCD and Next generation KIC protocol ( with -loops:frag_sizes and -loops:frag_files) call removed) and I do get properly modeled loops see pdb003_KICPS: Example pdb files below need to be renamed with a .gz extension and then uncompressed as they have been compressed using gzip
Thank for your time
Kamau
-
June 19, 2015 at 9:11 pm #11058
Anonymous
Looking at the code, it looks like at one point there is a specific check for extended loops, and it deliberately skips rebuilding loops that are not extended. I’m guessing in the models that look strange you’re running up against this, resulting in loops which are not fully remodeled.
The fix is simply to skip doing the non-extended runs for the protocols which result in badly modeled loops under those conditions. (Fragment insertion should be thorough, anyway, so there wouldn’t be an appreciable difference between doing with/without loop extension anyway.)
-
-
AuthorPosts
- You must be logged in to reply to this topic.