GROMACS version:2021.7
GROMACS modification: No
Here post your question
Hello everyone,
I’ve been trying to perform MD simulation for a protein-ligand complex in non-aqueous solvent; Since this solvent has a different isothermal compressibility factor to (1.3e-4bar-1) to that of water, should I use this value in my .mdp file or is it OK to use the compressibility factor of water?
I appreciate any help in advence.
Hi Ali,
You can find the answer here.
When you change the compressibility, do you get different properties?
Yes, use the value for your solvent.
But using an incorrect compressibility with effectively only give a different coupling time than the tau you supply.
Thank you @hess
is there any reference to find the right value for tau? or is there any way to determine it myself? According to some research articles the experimentally measured value for acetonitrile at 45C is nearly 1.3e-4 (bar-1) but unfortunately I have not found anything related to tau and as far as I know almost all MD simulations related to protein-ligand complex have been biased in one way or another towards aqueous medium.
Thank you @yehon
Since I have not found any reference to compare my results with, I honestly don’t know; but after equilibration in NVT and NPT ensembles the density of system is nearly equal to that of acetonitrile ( this actually happens when I use 4.5e-5) but when I change the compressibility to 1.3e-4 I my system does not equilibriate properly resulting in an error in MD production run. I don’t know if the density of the system counts as a reliable way :)
You choose tau_p yourself. For the c-rescale barostat I recommend using 2 ps, for parinello-rahman at least 5 ps.
I really appreciate your help