- This topic has 7 replies, 4 voices, and was last updated 10 years ago by Anonymous.
November 5, 2013 at 10:09 am #1745Anonymous
I am trying to dock a peptide to a receptor using FlexPepDock Refinement and generated 600 models,my top five models lie within sub-angstrom RMSDbb of the native structure.The lowest energy is -557.578.
Here is my flags:
Performing the command:
score.linuxgccrelease @flags > score.log
I got the total score of the native structure(-316.451).I’m confused that the total score of the native is greater than the models(-316.451>-557.578).Does it means there is something wrong with the scoring function(SCOREFUNCTION:standard SCOREFUNCTION PATCH:score12) or the default scoring function for FlexPepDock search is unfitted. My files have been attached.
Any replies is greatly appreciated.
November 5, 2013 at 10:33 am #9464Anonymous
This is a common observation: the crystal structure (native structure) is usually scored with less favorable, higher energy than any structure optimized using the Rosetta forcefield.
You should compare therefore your models to a near native structure obtained by local optimization, which would change the structure minimally, but reduce the energy significantly, by e.g. relieving clashes.
For FlexPepDock applications, you may use the flag -flexPepDockingMinimizeOnly to do so. More details can be found in the manual at https://www.rosettacommons.org/manuals/archive/rosetta3.5_user_guide/d7/d14/_flex_pep_dock.html
November 5, 2013 at 3:44 pm #9467Anonymous
Just to add to the discussion: Which version of Rosetta are you using? The default energy function and settings changed from 3.5 to the more recent weekly releases (talaris2013). If you are using a weekly, to get the same options for standard/score12 you will need to pass the option -restore_pre_talaris_2013_behavior
You may also want to look into relaxing your native structure first using the relax application with the option -constrain_relax_to_start_coords. This will use the more rigorous minimization of the relax algorithm, while not messing with your structure too much.
November 5, 2013 at 5:27 pm #9468Anonymous
could you please upload the the native or the optimized native and one of the decoy structures, so I can look into it?
November 7, 2013 at 8:59 am #9477Anonymous
The version of Rosetta i use is 3.4.Here is my PDB files and score files.
November 7, 2013 at 11:25 am #9478Anonymous
In order to compare structures that differ only at the interface, you need to use the same background (or look at different score terms, such as peptide score or interface score only). In your case it looks like you did not prepack the native when you did the optimization run. Therefore the receptor side chains might differ. I can see residues in the native structure that have unfavorable energy values resolved in the decoy (e.g. high repulsive value for His 263 that is flipped in the prepacking step) and are not at the interface.
November 5, 2013 at 2:07 pm #9465Anonymous
Thank you very much for your response Ora.
You’re right.I use the flag -flexPepDockingMinimizeOnly to perform a short minimization of the native structure and the energy of native structure reduces obviously(total_score:-494.152),but is still higher than the models.I wonder if i can modify the default scoring function and how to modify it.
November 5, 2013 at 3:03 pm #9466Anonymous
Before changing the scorefunction, I’d check to see what’s giving you the difference. If you have your models output as PDBs, there should be a table at the end of the PDB which shows the energies for each residue broken down by scoreterm. You can try comparing the two tables to see if it’s a particular residue/scoreterm which is giving you the difference. Once you know that, you can examine the native and the docked structures to see if you can rationalize the difference. This hopefully will give you insight as to what you should be doing differently.
- You must be logged in to reply to this topic.