Scoring

Hello,

I'd like to decompose the energy in one body energy terms (for each residue) and pairwise energy terms for residue-residue interactions.
I tried to check if I could decompose these terms and sum them up to retrieve the score returned by the scoring function.
The problem is that I try 3 different ways of computing the energy, and I get 3 different numbers.
Could someone please tell me what's wrong ?

Here is my code:

Is it possible to approxmate the CHARMM energy (function) by the rosetta scoring function?

As a starting point I could work with fixed angles and bond lengths, so I guess setting the weghts for fa_atr, fa_rep and fa_elec to 1 and all the rest to zero is a good start. But is it possible to somehow load the CHARMM diheral parameters?

Thank you for your help,
Ajasja Ljubetič

I'm using MPI with the PeptideDeriver pilot app (which I reworked to support JD2). I am running into some problems:

1. nodes that have finished jobs stall with 100% cpu usage while the entire MPI run isn't complete.
2. some tracer.log files don't contain all of the output -- for instance, the master node might have gotten a success message from a node, but not all the output of that job appear in the tracer.log file (the last bit is missing)
3. some jobs aren't finishing.

We are interested in creating a protein that will bind to a hormone at physiological pH, but unbind to the hormone at acidic pH, encountered when the complex is cycled through the cell. To do this, we would like to use the pH_mode functionality developed by Kilambi and Gray (Biophys Journal 103:587-595, 2012).