Member Site › Forums › Rosetta 3 › Rosetta 3 – General › Is it possible to do Abinitio folding from a linear structure.?
- This topic has 11 replies, 3 voices, and was last updated 13 years, 8 months ago by
Anonymous.
-
AuthorPosts
-
-
April 16, 2011 at 4:13 pm #874
Anonymous
Hi, I was thinking is it possible I generate a linear structure with AMBER and then use it in Abinitio module for de novo folding or I can generate a linear structure in Rosetta and use it.? Because from what I see, Abinitio only needs the fasta format and an input structure. So if my input structure is from PDB, then it doesn’t make any sense to do de novo folding. In addition, how can I generate a fasta file with a certain protein seq i have in mind.? Thanks.
-
April 16, 2011 at 5:03 pm #5438
By the way, I also realised that Robetta is unable to generate fragments for protein sequences with less than 20 residues, so how can I generate fragments for a nonapeptide if I want to fold it.?
-
April 16, 2011 at 8:12 pm #5439
To your first question: ab initio does not require an input structure. It can take an input for RMSD benchmarking (where the correct answer is known, and you are testing if Rosetta gets it right). It only needs the fasta file.
To your second question: A nonapeptide is unlikely to have a rigid structure in solution, so folding it ab initio doesn’t make sense, and folding it with fragments drawn from _folded_ proteins is inappropriate. Why do you want a free nonapeptide “structure”?
-
April 17, 2011 at 1:11 pm #5441
Oh, because I’m attempting to dock a nonapeptide onto a protein receptor, and I used the linear structure of the nonapeptide and it just seems to stay linear which is definitely inaccurate. Therefore, I’m thinking of folding the nonapeptide a bit, before I dock it onto the protein.
-
April 17, 2011 at 1:37 pm #5442
Were you using regular docking or the flexpepdock app? Flexpepdock does flexible peptide docking.
-
April 17, 2011 at 2:24 pm #5443
I was using the flexpepdock module.
-
April 17, 2011 at 2:30 pm #5444
Is it staying linear, or rigid? If it’s rigid, there’s a problem. If it’s linear, that’s probably what Rosetta thinks is best. I’ve brought this to the attention of the flexpepdock folks.
-
April 17, 2011 at 3:50 pm #5445
It stays rigid for the majority of the nonapeptide, while a small portion bent. Because I have done an MD simulation on this protein-peptide complex, so that’s why I know that Rosetta didnt give me the results I wanted. Of course, it could be due to my flags. But I was thinking, even if Im using the basic flags, I would expect the linear structure of the peptide to change.
-
April 21, 2011 at 5:17 am #5474
Okie, my cmd-line is FlexPepDocking.linuxgccrelease @flags.basic > dock.log &. Anyway, I just read the manual again and I realise that the starting structure is to be an approximate of the final complex structure that I’m expecting. If that’s the case, then it defeats my purpose in docking a linear peptide to the protein and hoping to see what and where would my final conformation of my peptide be like. So is it true that way or did I misinterpret anything.? Btw, I have added the extend_peptide flag but it seems to have no changes in the results at all.
Filename: flags.basic
-s input/complex.pdb
-native input/complex.pdb
-rbMCM
-torsionsMCM
#-lowres_preoptimize
-flexpep_prepack
#-flexPepDockingMinimizeOnly
-ex1
-ex2aro
-use_input_sc
-nstruct 20
#-unboundrot input/unbound.pdb
-database /home/u1/Programs/rosetta3.2.1_bundles/rosetta_database
-rep_ramp_cycles 10
-mcm_cycles 8
-smove_angle_range 6.0
-extend_peptide true -
April 25, 2011 at 6:43 am #5484
Okie, if let’s say I have the native structure, but I wish to only change 1 or 2 residues, do you think it might still give me a relatively accurate result.? since I presume that there shouldn’t be major changes to the secondary structure of the peptide.
-
April 20, 2011 at 10:00 am #5459
Hi,
Could you provide the cmd-line you used? (Just to assess if the problem is not there…)
By the way, flexpepdock has an intrinsic option to let you start the simulation from an extended peptide (“-extend_peptide true” : start the protocol with the peptide in extended conformation (neglect original peptide conformation ; extend from the anchor residue)
Cheers,
Nir. -
April 21, 2011 at 3:30 pm #5480
You should remove the -flexpep_prepack flag, this flag should be used once to create ‘prepacked’ structure and then for the production run you should not use it.
Regarding your question, if you believe the peptide is ‘close to linear’ in its native conformation, FlexPepDock should still serve your needs, however, it wouldn’t fold it into a helix or a beta hair-pin. If this is the level of sampling you’re looking for, the next rosetta release will include the FlexPepDock ab-initio protocol that would have a much larger sampling range.
-
-
AuthorPosts
- You must be logged in to reply to this topic.