MantisBT - Rosetta |
| View Issue Details |
|
| ID | Project | Category | View Status | Date Submitted | Last Update |
| 0000268 | Rosetta | [All Projects] Crash | public | 2013-05-07 15:20 | 2013-05-08 17:59 |
|
| Reporter | momeara | |
| Assigned To | renfrew | |
| Priority | normal | Severity | crash | Reproducibility | always |
| Status | resolved | Resolution | fixed | |
| Platform | | OS | | OS Version | |
| Product Version | | |
| Fixed in Version | Trunk | |
| Application(s) Affected | all |
| Command Line Used | score_jd2.linuxgccrelease -s inputs/1jf8FH.pdb -score:weights mm_std.wts -ignore_unrecognized_res |
| Developer Options | Confirmed As Bug |
| Fixed in SVN Version | 1691dfed3fb982744aa0c7c71baf03fad73acccf |
|
| Summary | 0000268: K residues causes MMTwist to crash |
| Description | The K residue type (database/residue_type_sets/fa_standard/residue_types/metal_ions/K.params) has the potassium and two virtual atoms. This leads it to have 3 torsion angles, which are computed and stored in ResidueType::update_derived_data() line 1175. But when scoring with the MMTorsionEnergy score term, the MMTorsionLibrary::lookup() fails with
No parameters for X-X-X-X.
|
| Steps To Reproduce |
score_jd2.linuxgccrelease \
-s inputs/1jf8FH.pdb \
-score:weights mm_std.wts \
-ignore_unrecognized_res |
| Additional Information | |
| Tags | No tags attached. |
| Relationships | |
| Attached Files |  1jf8FH.pdb (204,777) 2013-05-07 15:20
https://carbon.structbio.vanderbilt.edu/mantisbt/file_download.php?file_id=30&type=bug
1h4xFH.pdb (361,614) 2013-05-07 15:45
https://carbon.structbio.vanderbilt.edu/mantisbt/file_download.php?file_id=31&type=bug |
|
| Issue History |
| Date Modified | Username | Field | Change |
| 2013-05-07 15:20 | momeara | New Issue | |
| 2013-05-07 15:20 | momeara | Status | new => assigned |
| 2013-05-07 15:20 | momeara | Assigned To | => renfrew |
| 2013-05-07 15:20 | momeara | File Added: 1jf8FH.pdb | |
| 2013-05-07 15:45 | momeara | File Added: 1h4xFH.pdb | |
| 2013-05-08 09:21 | renfrew | Note Added: 0000284 | |
| 2013-05-08 13:00 | renfrew | Note Added: 0000285 | |
| 2013-05-08 15:22 | renfrew | Note Added: 0000287 | |
| 2013-05-08 15:24 | momeara | Note Added: 0000288 | |
| 2013-05-08 15:35 | renfrew | Note Added: 0000289 | |
| 2013-05-08 16:21 | renfrew | Note Added: 0000290 | |
| 2013-05-08 17:59 | renfrew | Fixed in SVN Version | => 1691dfed3fb982744aa0c7c71baf03fad73acccf |
| 2013-05-08 17:59 | renfrew | Status | assigned => resolved |
| 2013-05-08 17:59 | renfrew | Fixed in Version | => Trunk |
| 2013-05-08 17:59 | renfrew | Resolution | open => fixed |
|
Notes |
|
|
|
This is actually a problem with the params file; the MMTorsion (mm_twist) code is fully compatible with virtual atoms. This has come up before with the other ions. Folks seem to think that the MMAtomType of "X" is equivalent to a (Rosetta)AtomType of "VIRT" when it is not. The MMTorsion code uses "X" internally as a wildcard value (e.g. X-CT2-CT2-X in the mm_torsion_params.txt would mean use these params for all torsions that have CT2 and CT2 as the central atoms).
The solution is to change the MMAtomType field of the params file for the virtual atoms to be "VIRT" just like it is for the rosetta AtomTypes. |
|
|
|
|
|
The problem of scoring potassium in 1jf8FH.pdb was fixed by adding potassium params to mm_atom_properties.txt and modifying the params files to reflect the additional atoms types and changing X to VIRT for the virtual atoms. |
|
|
|
|
The 1hxFH.pdb also fails in the MMTorsion lookup but for a different reason. It doesn't contain potassium. It does contain a phosphorlated serine labeled as SEP.
This was a bitch to track down but both...
database/chemical/residue_type_sets/fa_standard/patches/ser_phosphorylated.txt
database/chemical/residue_type_sets/fa_standard/residue_types/ptm-caa/SEP.params
...define a ResidueType with SEP as the three letter code. One has complete MMAtomTypes and the other does not. However even if the MMAtomTypes were fully defined. It probably would have crashed because we don't have MMTOrsion params for phosphate groups.
On a side note, the SEP side chain is missing a hydrogen. I assume that they are built with reduce? That looks like a bug. |
|
|
|
|
|
Should one of those be removed? |
|
|
|
|
So none the params in residue_types/ptm-caa are currently being read in by default. It was the however the first param file I found though and it was really confusing why there were still undefined MMAtomTypes because they are all set in there (they are not set correctly but they are present and not "X"). I later found the patch file that creates the same ResidueType and the MMAtomType were not complete there.
I am not sure who added either of them but one should probably go. I can try to add some torsion params for the phosphate group, pulling them from D/RNA params. To make it work.
When I first made all the NCAA I had params files for most of these phosphorylated and sulfonated residues but CHARMM did not have params for them so I never commited them. |
|
|
|
|
The problem with the phosphoserine in 1h4xFH.pdb was fixed by modifying the patch to add the correct MMAtomTypes. New MMTorsion params were added based on...
Feng, M.-H., Philippopoulos, M., MacKerell, Jr., A.D. and Lim, C. Structural Characterization of the Phosphotyrosine Binding Region of a High-Affinity aSH2 Domain-Phosphopeptide Complex by Molecular Dynamics Simulation and Chemical Shift Calculations. Journal of the American Chemical Society, 1996, 118: 11265-11277.
...and added to mm_torsions.txt. |
|