I want to calculate the preferential interaction coefficient and subsequently the Gibbs free energy of solvation for protein in water and in water-co-solvent mixture.
According to the equation specified in the literature, if I can generate a list of no. of solvent molecules (for each local/bulk region) for each frame using g_select tool in GROMACS, it will give me nw, nx for local and bulk regions per frame.
With this information, how can I calculate the preferential interaction coefficient? For the entire average of the equation given, can I take these list of values and for each frame, compute the value of interaction coefficient. Thereby, if I get a list of all such values for every frame, I can average it to get to the value of preferential interaction coefficient over all frames.
Is this method correct?
Does anyone have a code /script to calculate this parameter?
Any other suggestions?