Dear Bioexcel users,
I use PMX to calculate when Ser mutates into Phe, when calculating ΔG in water, the repeatability is very poor, I use the version of gromacs2020.7.
In addition, I would like to ask if the trajectory file generated after using pmx can be used to analyze RMSD, RMSF and the number of hydrogen bonds? Because when I use pymol to view the mutation site, the key display of the mutation site is wrong, what should I do about this situation?
What exactly do you mean with “poor repeatability”? What is the standard error that you get from multiple repeats? Is each individual simulation sufficiently converged?
Hybrid structures contain residues that can represent WT and Mut states. For A-state, you can just hide the dummies. For B-state, you would need to map the atoms to those of an actual structure in the B-state. If you omit the mutated residue, the rest of the structure follows the standard atom naming
Thank you for your answer.
I repeated it a total of six times, with a standard error of 6.29, and each simulation was 40ns, checking that the RMSD already fluctuated very little.
From these 5 runs I see that the standard error of the mean is ~2.5 kJ/mol which seems reasonable. Only water3 appears to be somewhat of an outlier. You could check visually if you spot anything different happening in that trajectory. Also, how long are your transitions?
Thanks, I will use VMD to view the track file of this water3, in addition, I set the extraction 100 frames when unbalanced conversion, and the transitions time per frame is 100ps.
Finally, I would like to ask if PMX is compatible with gromacs 2020.7? I get an error when using the gmx rms command when parsing steps, but I get it right when using gromacs 5.1.4.
pmx is compatible with gmx2020
Set your GMXLIB to the path with the mutation force fields to resolve these errors
Ok, would need to add dummy masses to the ff files.
Would it be an option to follow the suggestion that you see in the error description, i.e. use a .tpr file?
