- This topic has 1 reply, 2 voices, and was last updated 11 years, 1 month ago by Anonymous.
October 31, 2012 at 5:10 pm #1449Anonymous
I’m using Rosetta 3.4 and I want to make a homology model of a GPCR protein -> I know very well how the model should look like.
This is the input I’m using:
-loops:frag_files mc4_09 mc4_03 none
-loops:frag_sizes 9 3 1
The template pdb is s1p.txt and the model is mc4_S_001.txt.
What has gone wrong? Any ideas?
October 31, 2012 at 6:25 pm #8020Anonymous
“The” model? One of the common mistakes is not understanding that the abinito/homology modeling programs (well, all Rosetta programs, but abinitio/homology modeling in particular) are stochastic protocols. You can’t just run them for one output structure and expect the first output structure to be good. Typically you should make thousands of structures, and then do clustering and energy filtering to pick out the best structure. ( See http://www.meilerlab.org/index.php/jobs/resources -> Rosetta Tutorials for a introduction/overview of the process.)
The other issue you’re going to run into is the fact that it looks like you’re modeling a trans-membrane helical bundle. By default Rosetta (like nature) relies on hydrophobic packing to stabilize structures. With membrane proteins, you don’t have the apolar inside/polar outside to help organize things, so Rosetta is a bit behind the curve. The constraints to the template structure may correct for this, but if increasing the number of output structures doesn’t help, you might want to look into incorporating the membrane into the scoring. See http://www.rosettacommons.org/manuals/archive/rosetta3.4_user_guide/db/d38/membrane_abinitio.html for the description of membrane abinito.
That said, there’s some concern that the membrane scoring for Rosetta isn’t all that great. (For example, there was a membrane target in CASP10, and modeling it with membrane scoring didn’t give anything good, but modeling it with the as-yet-unpublished Baker lab CASP10 protocol did give good results.) I’d recommend using the regular homology modeling protocol to start with (but with a larger number of structures), and only try the membrane protocol adjustments if that isn’t working well.
- You must be logged in to reply to this topic.