Hello everyone! I'm new to the Rosetta software and I apologize in advance if someone has already asked this question before (I looked through the forum but couldn't find such topic).
I'm trying to refine a membrane protein structure. I obtained the crude structure from homology modeling with Modeller and then ran the mp_relax protocol of Rosetta Membrane (I used the available mp_relax files from the protocol capture folder of Rosetta 3.8, a distribution that was installed some time last August). Now I wan to refine further the loop portions of my protein. For that I tried to used the files from the loop_modelling tutorial that came with this distribution, but I encountered a potential problem:
If I refine a single loop, the program seems to be running fine and produces an otput where the geometry of the loop has changed with respect to the initial structure. However, if I ask for more than one loop to be refined at the same time, the code seems to start working on one of the loops (see below the fiirst loop that gets optimized)
kinematic initial perturb with start_res: 299 middle res: 313 end_res: 328
and I start getting lists of warnings of the kind:
core.kinematics.AtomTree: No proper DoF can be found for these four atoms: 44-1, 44-2, 44-3, 45-1!
where the four atoms belong to loops that are NOT undergoing the optimization at the moment. Once the code finishes with the first loop, it moves to the second one and then I get new warning messages for missing DoFs. However, this time, I don't have DoF messages for atoms from the very first loop AND from the second loop that is being currently optimized. This goes on and on, until I no longer get DoF messages during optimization of the last loop of the set. Then I get another iteration of loop refinement, as ROSETTA is building a second model (I've asked for nstruct 5) starting from the very first loop, and the DoF messages appear again.
So, my question is: Can I trust the structures obtained through this multiple loop refinement modeling or are the DoF warnings an indicator that something went very wrong during the simulation? I couldn't find much information about this on the internet - I guess it implies a problem with the Atom Tree, but I have no idea how it can be fixed.
When I looked at the 5 generated structures I discovered that all of the loops have underwent change with respect to the initial structure. However, initially I had some highly structured loops (from the .pdb files on which I based my un-refined homology model), whose structure has been destroyed during this refining protocol and they have become highly disordered. I don't think this is realistic (these structured loops appear in the crystal structures in a number of bacterial homologues of this protein, so they are obviously stable enough to get resolved by X-Ray diffraction). I would love to hear some opinions on this from users who are more experienced with ROSETTA's code and performance. Are there some options that I can explore in this case that would preserve the structure of the highly structured loops? Any keywords that can be of help?
Please let me know if there is any other information that you may need! Below are the contents of my rosetta flag file and the .loops file which has the list of loops that need to be optimized. I've used the loopmodel.linuxgccrelease executable on my office machine and I've also ran a serial ROSETTA calculation (with nstruct 2000) with the loopmodel.mpi.linuxiccrelease (of ROSETTA 3.8) on one of the clusters my group is using. The serial calculation behaved in the same way.
My flag file looks like this:
And my NIS_refine_loop.loops file looks like this:
#subtract 10 since first 10 residues are missing
LOOP 24 45 0 0 0
LOOP 102 116 0 0 0
LOOP 150 153 0 0 0
LOOP 168 177 0 0 0
LOOP 299 328 0 0 0
LOOP 362 369 0 0 0
LOOP 392 405 0 0 0
LOOP 427 431 0 0 0