Member Site › Forums › Rosetta 3 › Rosetta 3 – General › fixbb and HETATM
- This topic has 2 replies, 2 voices, and was last updated 11 years, 9 months ago by Anonymous.
-
AuthorPosts
-
-
January 14, 2013 at 3:50 pm #1496Anonymous
Hi,
I want to use fixbb on a pdb file that contains ions (the ions are classified as HETATM in the pdb file) in addition to the protein chain (whose atoms are classified as ATOM in the pdb file). My question is: If I use fixxb on the pdb file for repacking, will the algorithm consider the presence of the ions or will it ignore them?
If the fixbb algorithm ignores the ions as a default, how can I force it to include the ions in the calculation (I don’t know if it matters, but just so you’ll have the full picture, my protein is a transmembrane protein so I am using the -score:patch no_solvation flag with fixbb) ?Thank you,
Doran -
January 14, 2013 at 4:34 pm #8307Anonymous
Rosetta is somewhat capable of considering your ions in fixbb. To get the ions working, you need parameters files for them. Look in the folder rosetta_database/chemical/residue_type_sets/fa_standard/residue_types/metal_ions to see which metals Rosetta already has programmed. You may need to write your own parameters files if none of these are the right elements. Also check in rosetta_database/chemical/residue_type_sets/fa_standard/residue_types.txt to make sure your ions are turned on (not # in front of the line).
Secondarily, you need to make sure the ions in your input match the spelling and spacing the parameters file expects. For demonstration, I’ll substitute spaces in the PDB for underscores: if you have “MG_”, and it doesn’t work, try “_MG”. You may need to swap HETATM for ATOM.
Finally, the quality of Rosetta’s consideration is somewhat suspect. Rosetta doesn’t know about the orbitals and bonding patterns of ions, it will consider their van der Waals radius and that’s about it. You aren’t using solvation and the standard scorefunction’s pseudo-electrostatics (fa_pair) ignores ions. If you want beautifully tetrahedrally-coordinated zinc, you aren’t going to get it without some heavy-duty additions to the scorefunction. If you just want atoms taking up space to prevent some ion-dependent conformation from collapsing, it oughta work.
-
January 14, 2013 at 4:48 pm #8308Anonymous
Thank you very much,
Doran
-
-
AuthorPosts
- You must be logged in to reply to this topic.