Rosetta crashed while designing with PDB.rotamer #141
Replies: 2 comments 1 reply
-
2000 rotamers is a huge rotamer library. It could easily result in you
running out of memory at any of the following steps:
- Read-in of the rotamer library, if the dimensionality is large.
- Setup of the polycubic interpolation used for scoring.
- Construction of the list of candidate rotamer instances for a particular
packing problem.
- Precomputation of the interaction graph, which needs to store pairwise
energies for all of the potentially interacting rotamers (so if you had
just two positions each with just 2000 rotamers, this would be 4 million
pairwise interactions).
- Precomputation of other quantities stored by design-centric guidance
terms (voids_penalty, hbnet, buried_unsatisfied_penalty).
Try running on a single core without the -mute all flag to see what Rosetta
was trying to do when it ran out of memory.
…On Tue, Sep 3, 2024, 7:08 AM Muhammad Atta Fouad ***@***.***> wrote:
Hi,
I tried out "Non-canonical amino acid (NCAA) parameterization using
Rosetta" from the Meiler Lab 2023 tutorial, I used the fake_rotlib script
to make libraries for some amino acids, each of 2000 rotamers (as the
structure is large), and when I moved to the next part of the tutorial
"Macrocyclic peptide design in Rosetta" which relies on research article of
@vmullig <https://github.com/vmullig> Mulligan et. al. Computationally
designed peptide macrocycle inhibitors of New Delhi metallo-beta-lactamase
1, the submitted job crashed after rising of the occupied size RAM (250 GB
of the HPC), without any details about why that happened.
Note, that using only one or 2 amino acids worked without issues, only
more numbers resulted in that.
—
Reply to this email directly, view it on GitHub
<#141>, or
unsubscribe
<https://github.com/notifications/unsubscribe-auth/ABACZUGEWQ3DKLXPMAMQDW3ZUU75DAVCNFSM6AAAAABNRKCO4SVHI2DSMVQWIX3LMV43ERDJONRXK43TNFXW4OZXGEZTCMJXHE>
.
You are receiving this because you were mentioned.Message ID:
***@***.***>
|
Beta Was this translation helpful? Give feedback.
1 reply
Answer selected by
MuhammadAttaElhamouly
-
For more verbose output, you can also add:
…-out:level 500
On Tue, Sep 3, 2024, 8:14 AM Muhammad Atta Fouad ***@***.***> wrote:
Thanks, I'll check it out again according to your recommendations.
—
Reply to this email directly, view it on GitHub
<#141 (reply in thread)>,
or unsubscribe
<https://github.com/notifications/unsubscribe-auth/ABACZUC6N4X7BRMU43Q5SZDZUVHUVAVCNFSM6AAAAABNRKCO4SVHI2DSMVQWIX3LMV43URDJONRXK43TNFXW4Q3PNVWWK3TUHMYTANJSG43DSMY>
.
You are receiving this because you were mentioned.Message ID:
***@***.***>
|
Beta Was this translation helpful? Give feedback.
0 replies
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
-
Hi,
I tried out "Non-canonical amino acid (NCAA) parameterization using Rosetta" from the Meiler Lab 2023 tutorial, I used the fake_rotlib script to make libraries for some amino acids, each of 2000 rotamers (as the structure is large), and when I moved to the next part of the tutorial "Macrocyclic peptide design in Rosetta" which relies on research article of @vmullig Mulligan et. al. Computationally designed peptide macrocycle inhibitors of New Delhi metallo-beta-lactamase 1, the submitted job crashed after rising of the occupied size RAM (250 GB of the HPC), without any details about why that happened.
Note, that using only one or 2 amino acids worked without issues, only more numbers resulted in that.
Beta Was this translation helpful? Give feedback.
All reactions