You are here

Is it possible to do Abinitio folding from a linear structure.?

12 posts / 0 new
Last post
Is it possible to do Abinitio folding from a linear structure.?
#1

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.

Post Situation: 
Sat, 2011-04-16 09:13
monos_morpheus

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.?

Sat, 2011-04-16 10:03
monos_morpheus

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"?

Sat, 2011-04-16 13:12
smlewis

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.

Sun, 2011-04-17 06:11
monos_morpheus

Were you using regular docking or the flexpepdock app? Flexpepdock does flexible peptide docking.

Sun, 2011-04-17 06:37
smlewis

I was using the flexpepdock module.

Sun, 2011-04-17 07:24
monos_morpheus

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.

Sun, 2011-04-17 07:30
smlewis

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.

Sun, 2011-04-17 08:50
monos_morpheus

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.

Wed, 2011-04-20 03:00
londonir

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

Wed, 2011-04-20 22:17
monos_morpheus

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.

Thu, 2011-04-21 08:30
londonir

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.

Sun, 2011-04-24 23:43
monos_morpheus