Hello Reza
I had a similar issue before and the order in which you source files is indeed relevant. In order to get the OPC water parameters it needs to be sourced last to be prioritised. The order of priority was "bottom up”.
Best regards
// Gustaf
> On 14 Apr 2021, at 19:24, Reza Esmaeeli <goldeagle.93.gmail.com> wrote:
>
> We have been trying to use OPC water model in our simulations. It looks like the order of loading water and macromolecule forcefields is important. Since ff19SB loads parm10, it will overwrite parameters for OW from tip3p, but not EP. This leads to abnormal behavior of water molecules which act as if they are positionally restrained. We have tested this with OL3 and OL15 as well. The pictures below show coordinates of two water molecules across a 50-ns trajectory with the aforementioned order of loading forcefields:
>
> source leaprc.water.opc
> source leaprc.protein.ff19SB
>
>
>
> source leaprc.protein.ff19SB
> source leaprc.water.opc
>
>
> I wonder in case of other 3-point water models, tip3p could be overwriting SPC parameters and while it might diffuse properly, the actual values may not be correct.
> Any ideas if this should be included as a warning in tleap or some other way of preventing it?
> _______________________________________________
> AMBER mailing list
> AMBER.ambermd.org
> http://lists.ambermd.org/mailman/listinfo/amber
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Thu Apr 15 2021 - 00:30:02 PDT
