Member Site › Forums › Rosetta 3 › Rosetta 3 – General › abinitiorelax: No funnel
- This topic has 6 replies, 3 voices, and was last updated 14 years, 1 month ago by
Anonymous.
-
AuthorPosts
-
-
July 27, 2010 at 1:14 pm #639
Anonymous
I’m attempting to do an ab initio prediction. I have created 360.000 structures using: AbinitioRelax.linuxgccrelease with the following flags:
-in:file:fasta inputs/mINT.fasta
-in:file:frag9 inputs/aat000_03_05.200_v1_3
-in:file:frag3 inputs/aat000_09_05.200_v1_3
-out:pdb
-out:file:silent OUT
-out:nstruct 1
-abinitio
-database /x/Rosetta3.1_Bundles/rosetta_databaseI assume that the rosetta energy score is the one in the first column of the score.fsc file just after “SCORE: “. All the scores I get are positive?!
I have identified the lowest scoring file (score = 1.7) and calculated RMSDs of all other structures compared to this. When making a plot of the (assumed) score as function of RMSD I get cloud but no funnel!?
I hope you can help!
-
July 27, 2010 at 3:40 pm #4513
Anonymous
Caveat: I don’t use ab initio, so I’m guessing.
A) I’m not sure what happens when you use both out::file::silent and out::pdb at the same time. Are you getting both PDBs and a silent file? and a score file too?
Are the columns not labeled in the score.fsc file…? What do the first few lines look like? I’m not familiar with .fsc, is that .fasc?
C) It’s perfectly possible you have a cloud with no funnel, it happens. What’s the spread in values? If the range is 1.2-2.3, then it’s probably not rosetta total score in that column. If it’s 1.7-500, then it probably is.
D) It’s possible to get all positive scores. It means Rosetta hates all the models it made.
-
July 28, 2010 at 9:02 am #4517
Anonymous
First line in scores.fsc is the sequence. Second line:
SCORE: score vdw cenpack pair env cbeta rg hs_pair ss_pair rsigma sheet clashes_total clashes_bb time descriptionThis there any documentation on the output to this file?
My scores range from 1.7-194.
It seems unlikely that running more predictions will fix the problem! Is the any way to get Rosetta to search other conformations or must I conclude that rosetta cannot to the prediction of this 160 aa protein?!
-
July 28, 2010 at 2:39 pm #4521
Anonymous
The first issue is that you have centroid poses, not fullatom poses. It is expected that centroid poses have positive scores. I thought you had fullatom poses. So, at least the only-positive scores aren’t a worry.
You might be able to try running abrelax with a pre-supplied secondary structure prediction, that might improve Rosetta’s results. You could also try running the post-abinitio relax phase (that’s what it’s abrelax) and/or getting fullatom results, but I don’t know if that will help any.
-
July 31, 2010 at 10:37 am #4535
Anonymous
I realized the missing relaxation and did fix this. This resulted in all negatives (range -330 to -180). But still no funnel
Is there any chance that changing from centroid poses to fullatom poses will fix my problem?! In any case: How is this done?
Thank you!
-
August 1, 2010 at 7:23 pm #4538
Anonymous
The manual doesn’t say. You can try the -help option; I believe the abrelax authors liked that particular documentation system (as opposed to the manual) and it may offer more advice.
-
August 5, 2010 at 8:28 pm #4553
Anonymous
>You might be able to try running abrelax with a pre-supplied secondary structure prediction, that might improve Rosetta’s results.
Hi,
how do you give a secondary structure prediction to AbinitioRelax.linuxgccrelease? it is not obvious from online manual nor “-help” command option.
Thanks!
-
-
AuthorPosts
- You must be logged in to reply to this topic.