That's great, @mensur. Thank you. I have a run that operates on 1ubq_Repair.pdb and produces a 'Dif_1ubq_Repair.fxout' now, even with numberOfRuns=3:
FoldX 5.1 (2011)
by the FoldX Consortium
Jesper Borg, Frederic Rousseau, Joost Schymkowitz,
Luis Serrano and Francois Stricher
------------------------------------------------------
PDB file analysed: batch
Output type: BuildModel
Pdb total energy Backbone Hbond Sidechain Hbond Van der Waals Electrostatics Solvation Polar
1ubq_Repair_1.pdb 0.998679 0.00173548 0 0.867633 0 -0.917758 1.44122 ->
1ubq_Repair_1_0.pdb 1.1988 0.018838 0 0.526343 0 -0.354662 1.02289 0 ->
1ubq_Repair_1_1.pdb 1.14399 0.0682475 -5.41737e-09 0.35607 -0.000133908 -0.198468 0.673921 >
1ubq_Repair_1_2.pdb 1.16184 0.0781384 0 0.347325 0 -0.118711 0.719832
However, it's not clear that the mutation has been done correctly.
The original 1ubq.pdb at position 8:
$ awk '$6 == 8 {print $0}' /tmp/foldx_temp/1ubq_Repair.pdb
ATOM 59 N LEU A 8 30.132 40.069 18.642 1.00 9.84
ATOM 60 CA LEU A 8 29.607 41.180 19.467 1.00 14.15
ATOM 61 C LEU A 8 30.075 42.538 18.984 1.00 17.37
ATOM 62 O LEU A 8 29.586 43.570 19.483 1.00 17.01
ATOM 63 CB LEU A 8 29.919 40.890 20.938 1.00 16.63
ATOM 64 CG LEU A 8 28.876 40.104 21.721 1.00 18.88
ATOM 65 CD1 LEU A 8 28.679 38.712 21.146 1.00 19.31
ATOM 66 CD2 LEU A 8 29.258 39.960 23.190 1.00 18.59
I understand the mutation should be in the following files (based on the .fxout file): 1ubq_Repair_1.pdb, 1ubq_Repair_1_0.pdb, 1ubq_Repair_1_1.pdb, 1ubq_Repair_1_2.pdb. Please help my misunderstanding -
$ awk '$6 == 8 {print FILENAME,"|", $0}' /tmp/foldx_temp/1ubq_Repair_1*.pdb
/tmp/foldx_temp/1ubq_Repair_1_0.pdb | ATOM 59 N VAL A 8 30.132 40.069 18.642 1.00 9.84
/tmp/foldx_temp/1ubq_Repair_1_0.pdb | ATOM 60 CA VAL A 8 29.607 41.180 19.467 1.00 14.15
/tmp/foldx_temp/1ubq_Repair_1_0.pdb | ATOM 61 C VAL A 8 30.075 42.538 18.984 1.00 17.37
/tmp/foldx_temp/1ubq_Repair_1_0.pdb | ATOM 62 O VAL A 8 29.586 43.570 19.483 1.00 17.01
/tmp/foldx_temp/1ubq_Repair_1_0.pdb | ATOM 63 CB VAL A 8 30.013 41.002 20.941 1.00 16.47
/tmp/foldx_temp/1ubq_Repair_1_0.pdb | ATOM 64 CG1 VAL A 8 29.191 39.901 21.592 1.00 16.47
/tmp/foldx_temp/1ubq_Repair_1_0.pdb | ATOM 65 CG2 VAL A 8 31.484 40.628 21.047 1.00 16.47
/tmp/foldx_temp/1ubq_Repair_1_1.pdb | ATOM 59 N VAL A 8 30.132 40.069 18.642 1.00 9.84
/tmp/foldx_temp/1ubq_Repair_1_1.pdb | ATOM 60 CA VAL A 8 29.607 41.180 19.467 1.00 14.15
/tmp/foldx_temp/1ubq_Repair_1_1.pdb | ATOM 61 C VAL A 8 30.075 42.538 18.984 1.00 17.37
/tmp/foldx_temp/1ubq_Repair_1_1.pdb | ATOM 62 O VAL A 8 29.586 43.570 19.483 1.00 17.01
/tmp/foldx_temp/1ubq_Repair_1_1.pdb | ATOM 63 CB VAL A 8 30.013 41.002 20.941 1.00 16.47
/tmp/foldx_temp/1ubq_Repair_1_1.pdb | ATOM 64 CG1 VAL A 8 29.193 39.898 21.591 1.00 16.47
/tmp/foldx_temp/1ubq_Repair_1_1.pdb | ATOM 65 CG2 VAL A 8 31.485 40.631 21.047 1.00 16.47
/tmp/foldx_temp/1ubq_Repair_1_2.pdb | ATOM 59 N VAL A 8 30.132 40.069 18.642 1.00 9.84
/tmp/foldx_temp/1ubq_Repair_1_2.pdb | ATOM 60 CA VAL A 8 29.607 41.180 19.467 1.00 14.15
/tmp/foldx_temp/1ubq_Repair_1_2.pdb | ATOM 61 C VAL A 8 30.075 42.538 18.984 1.00 17.37
/tmp/foldx_temp/1ubq_Repair_1_2.pdb | ATOM 62 O VAL A 8 29.586 43.570 19.483 1.00 17.01
/tmp/foldx_temp/1ubq_Repair_1_2.pdb | ATOM 63 CB VAL A 8 30.013 41.002 20.941 1.00 16.47
/tmp/foldx_temp/1ubq_Repair_1_2.pdb | ATOM 64 CG1 VAL A 8 29.289 39.814 21.555 1.00 16.47
/tmp/foldx_temp/1ubq_Repair_1_2.pdb | ATOM 65 CG2 VAL A 8 31.511 40.761 21.056 1.00 16.47
/tmp/foldx_temp/1ubq_Repair_1.pdb | ATOM 59 N ALA A 8 30.132 40.069 18.642 1.00 9.84
/tmp/foldx_temp/1ubq_Repair_1.pdb | ATOM 60 CA ALA A 8 29.607 41.180 19.467 1.00 14.15
/tmp/foldx_temp/1ubq_Repair_1.pdb | ATOM 61 C ALA A 8 30.075 42.538 18.984 1.00 17.37
/tmp/foldx_temp/1ubq_Repair_1.pdb | ATOM 62 O ALA A 8 29.586 43.570 19.483 1.00 17.01
/tmp/foldx_temp/1ubq_Repair_1.pdb | ATOM 63 CB ALA A 8 29.919 40.890 20.938 1.00 16.47
It looks like the permutation is done correctly (L --> V) on 1ubq_Repair_1_0.pdb, 1ubq_Repair_1_1.pdb, 1ubq_Repair_1_2.pdb but not on 1ubq_Repair_1.pdb (L --> A).
I would like to accept the answer that I shouldn't worry about 1ubq_Repair_1 and just take my ddg values from 1ubq_Repair_1_*. So, I tested again with numberOfRuns=1, in which case I only get the 1ubq_Repair_1.pdb file. The .fxout file:
FoldX 5.1 (2011)
by the FoldX Consortium
Jesper Borg, Frederic Rousseau, Joost Schymkowitz,
Luis Serrano and Francois Stricher
------------------------------------------------------
PDB file analysed: batch
Output type: BuildModel
Pdb total energy ...
1ubq_Repair_1.pdb 1.1988 ...
and the mutation position:
$ awk '$6 == 8 {print FILENAME,"|", $0}' /tmp/foldx_temp/1ubq_Repair_1.pdb
/tmp/foldx_temp/1ubq_Repair_1.pdb | ATOM 59 N VAL A 8 30.132 40.069 18.642 1.00 9.84
/tmp/foldx_temp/1ubq_Repair_1.pdb | ATOM 60 CA VAL A 8 29.607 41.180 19.467 1.00 14.15
/tmp/foldx_temp/1ubq_Repair_1.pdb | ATOM 61 C VAL A 8 30.075 42.538 18.984 1.00 17.37
/tmp/foldx_temp/1ubq_Repair_1.pdb | ATOM 62 O VAL A 8 29.586 43.570 19.483 1.00 17.01
/tmp/foldx_temp/1ubq_Repair_1.pdb | ATOM 63 CB VAL A 8 30.013 41.002 20.941 1.00 16.47
/tmp/foldx_temp/1ubq_Repair_1.pdb | ATOM 64 CG1 VAL A 8 29.191 39.901 21.592 1.00 16.47
/tmp/foldx_temp/1ubq_Repair_1.pdb | ATOM 65 CG2 VAL A 8 31.484 40.628 21.047 1.00 16.47
I am encouraged because:
- L --> V as it should be
- total energy is 1.1988, which is the same as the total energy of 1ubq_Repair_1_0.pdb when numberOfRuns = 3.
So, it seems that numberOfRuns=1 produces a good result but numberOfRuns>1 is only valid for the *_Repair_<integer>.pdb files. Is that a valid conclusion that aligns with your experience, or am I still doing something wrong because numberOfRuns gives unexpected behavior?
The program is definitely working. It doesn't matter that it didn't produce
TERandENDlines at the end. Those are optional in PDB files.For the rest, it is almost impossible to read this as the whole thing is improperly formatted. If you put three backtick signs (like this `, on the same key in the upper left corner as ~ or tilde) as a single line at the start of the output you want to format, and the same three backtick signs as a single line at the end, it will be formatted like this:
Someone might be able to help from a properly formatted post.
By the way, a simple Google search will give you a clear clue how to fix this:
https://www.google.com/search?q=Sorry+you+specified+an+option_FOLDXDEFINED+for+the+order+of+mutations+that+does+not+exist%2C+should+be+_USERDEFINED+or+_FOLDXDEFINED&sca_esv=1e4238bc8dc008dd&rlz=1C1UEAD_enUS1016US1016&ei=Fv7laPewF6WQm9cPtaOg6QE&ved=0ahUKEwj3k8TC9pOQAxUlyOYEHbURKB0Q4dUDCBA&uact=5&oq=Sorry+you+specified+an+option_FOLDXDEFINED+for+the+order+of+mutations+that+does+not+exist%2C+should+be+_USERDEFINED+or+_FOLDXDEFINED
Responding to your answer.
Request: