Hi all, I'm trying to run relax on a pdb which has ACE and NME terminal caps. Is there a way I can preserve those in the output? I'm not an expert Rosetta user, but here are some more specifics of what I've tried:
- If I run relax using simple settings, it looks like Rosetta sees my ACE and NME and doesn't complain. In the output pdb's score section, the residues were renamed ILE:AcetylatedNtermProteinFull_1 and LEU:MethylatedCtermProteinFull_134. The problem is that the ACE and NME atoms were squashed into their neighboring ILE and LEU residues meaning that those residues are now incomprehensible to other non-Rosetta software packages I'd like to use. I'd like to get legitimately standard ILE and LEU back out again, ideally with those cap atoms placed back into separate ACE and NME residues. Is there a way to do that within Rosetta?
Based on advice in the structure preparation section, I thought I'd try the -chemical:exclude_patches option. It seems to work in that those "merged" residues aren't appearing in the trace output, however, I'm now getting an error:
ERROR: Error! Could not open file ncaa_rotamer_libraries/n_methyl_amino_acid/LEU_Nmethylated.rotlib for read.
which I can't get rid of even if I try to use a resfile to not repack that LEU.
I looked back at the other relax advised options and found that if I add the -correct flag, then I don't get the error and it appears like relax proceeds as expected. According to the documentation, the -correct flag shouldn't be used with the newer score functions, so I don't know why this would make it not need those rotamers? Regardless, the output pdb appears to drop the ACE and NME completely, which is not ideal, but preferred over the non-standard ILE and LEU residues.
I'm not sure if I'm hitting a bug or not... but any advice is welcome!
Here's the version of Rosetta I'm using:
core.init: Rosetta version: rosetta.binary.linux.release-171 r171 2018.09+release.333d996 333d99699777cd6a50d1c1736bacd689f8f1d1df https://www.rosettacommons.org 2018-02-27T09:57:00.705373