Member Site › Forums › Rosetta 3 › Rosetta 3 – General › How does one prevent Rosetta from connecting two separate chains?
- This topic has 7 replies, 3 voices, and was last updated 12 years, 3 months ago by
Church.
-
AuthorPosts
-
-
December 15, 2011 at 10:08 pm #1119
Anonymous
I have a homodimer of two separate chains, each one 31 residues long. I am attempting to run a comparative modeling protocol. Every pdb output has the two chains connected, which screws up subsequent relaxation.
How do I change the input fasta file, align file, and/or flags file to indicate that there are two separate chains?
(The two separate chains are: KMQMLKLDKENALDRAEQAEADKKAAEERSK and KMQMLKLDKENALDRAEQAEADKKAAEERSK
Current files:
fasta file:
> Tm7_37
KMQMLKLDKENALDRAEQAEADKKAAEERSKKMQMLKLDKENALDRAEQAEADKKAAEERSKalignment file:
Tm7_37 1 KMQMLKLDKENALDRAEQAEADKKAAEERSKKMQMLKLDKENALDRAEQAEADKKAAEERSK
2zta 1 RMKQLEDKVEELLSKNYHLENEVARLKKLVGRMKQLEDKVEELLSKNYHLENEVARLKKLVGcommand file:
/sware/rosetta/linux/rosetta_source/bin/minirosetta.linuxgccrelease @Tm7_37.flags -database /sware/rosetta/linux/rosetta_database -run:constant_seed -nodelay -out:overwrite true > Tm7_37RosCM.logflags file:
-run:protocol cm-in:file:fasta Inputs/Tm7_37.fasta
-in:file:alignment Inputs/Tm7_37.2zta.aln
-in:file:template_pdb Inputs/2ZTAnoHET.pdb
-in:file:extended_pose 1
-out:nstruct 1-in:file:fullatom
-out:file:fullatom
-loops:frag_sizes 9 3 1
-loops:frag_files Inputs/aat000_09_05.200_v1_3 Inputs/aat000_03_05.200_v1_3 none
-loops:remodel quick_ccd
-loops:relax fastrelax-cm:min_loop_size 4
-cm:loop_close_level 0-out:file:silent_struct_type pdb
-unmute protocols:looprelax
-run:test_cycles
…
Thanks,
-Jerry -
December 16, 2011 at 6:03 pm #6426
Anonymous
By the way, adding “-read_all_chain” to either the command file or flag file reulted in an ERROR
ERROR: Option matching -read_all_chains not found in command line top-level context
-
December 19, 2011 at 12:08 am #6441
Anonymous
Is it expected to be a symmetric homodimer? Can you predict on one chain and dock? What about fold_and_dock? I don’t think the basic ab initio app can do this.
-
December 20, 2011 at 3:59 pm #6452
Anonymous
It is expected be (part of ) a parallel alpha-helical homodimer, being nearly in-register but having a small (~1 Angstrom) axial offset from perfect two-fold symmetry.The two alpha-helical chains are intimately interacting with one another throughout their lengths in a coiled coil.
I never even attempted ab initio modelling for a number of reasons, foremost is that there is not enough information in the sequence of the peptide to indicate that it is a near in-register parallel coiled coil. (In vivo, this peptide is a fragment of a much longer sequence, and the information for the near in-register and parallel nature of the coiled coil resides elsewhere, beyond the fragment.)
Hence, I had thought to do a comparative modeling, based on a known structure that is a parallel in-register homodimer, and then, ?in the cm protocol? relax/refine the structure to see what Rosetta predicts, i.e. would Rosetta introduce the expected 1.0 Angstrom axial shift.
It seems a bit unusual that I can’t figure out how to tell Rosetta’s comparative modelling that the two chains are not covalently attached to each other, but in any case, perhaps I could try docking one chain to the other. Do you think that docking one alpha-helix to another to form an initmate coiled coil is appropriate for docking? Also, do you have a Rosetta3.3 compatible docking protocol that takes as input a modeled structure that is close to correct?
Thanks so much, -Jerry
-
December 22, 2011 at 8:47 pm #6463
Anonymous
Thanks. I had thought of putting in a long fake linker, but was hoping to avoid it. But perhaps it is worth a try.
-
June 1, 2012 at 7:40 am #7170
Church
Thanks smlewis, adding a faux GSGSGS linker and then just deleting it afterwards has worked out! I was trying to follow the initial guidelines provided by the topic starter but failed. I was in a hurry to fix the issue so my work really resembled typing speed test
Now it works well, thanks a bunch!
-
December 22, 2011 at 1:09 am #6461
Anonymous
You’re right that for a tightly coiled-coil, docking is unlikely to work well.
Ab initio may well work fine on this problem – the issue is that the basic abinitio_relax executable isn’t used much by the developers anymore, so nobody knows its quirks, and the newer, fancier abinitio protocols are either unreleased or undocumented.
I’ve sent this along to one of the comparative modeling people who may know a way to do it…
-
December 22, 2011 at 1:58 am #6462
Anonymous
The answer seems to be, “you can’t get there from here”. Ab initio is meant for single-domain single-chain predictions. I guess you can try something like adding a faux GSGSGS linker and then just deleting it afterwards…?
-
-
AuthorPosts
- You must be logged in to reply to this topic.