Hi Noureen,
In the tutorial all replicas started from the exact same structure and protonation states. But it is also possible to let them start from individually equilibrated structures at different pH, which may have different states. You just have to modify the inputs individually in the group file. Both ways make sense to me. The first is more convenient (needs less equilibrium) but some replicas may experience a drastic conformation change at the beginning, especially if their pHs are further away from the pH of choice in the previous equilibrium calculation.
To restart, it makes more sense to proceed from the respective protonation states at the end of the previous trajectories, for a smooth transition from the previous and the new simulations.
Best regards,
Amy
From: Abdelrahman, Noureen via AMBER <amber.ambermd.org>
Date: Saturday, June 8, 2024 at 6:03 PM
To: AMBER Mailing List <amber.ambermd.org>
Subject: [AMBER] Restarting pH REMD
Hello all,
In Jason Swail's pH replica exchange tutorial (linked here: https://urldefense.com/v3/__http://jswails.wikidot.com/ph-remd__;!!KGKeukY!wOc8Qyq17CCFYVTtfvM8yYGX6ITnIl6cKzmqY7oj_hcLYHXA1X5SXqshWeglCcB4ZUitmBmu8zMmOrbo3Dr_ZtWrePQ$<https://urldefense.com/v3/__http:/jswails.wikidot.com/ph-remd__;!!KGKeukY!wOc8Qyq17CCFYVTtfvM8yYGX6ITnIl6cKzmqY7oj_hcLYHXA1X5SXqshWeglCcB4ZUitmBmu8zMmOrbo3Dr_ZtWrePQ$> ), we generate one cpin file to include in the groupfile for all replicas using the line below:
bash$ cpinutil.py -p 1AKI.prmtop -resname AS4,GL4,HIP -igb 2 -o 1AKI.cpin
However, in the restart, the cprestrt file from each replica becomes the cpin file for that same replica in the restart. Can someone please explain why we initially use one cpin file for all replicas but then in the restart, we use a different cpin file for each?
Thanks,
Noureen
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
https://urldefense.com/v3/__http://lists.ambermd.org/mailman/listinfo/amber__;!!KGKeukY!wOc8Qyq17CCFYVTtfvM8yYGX6ITnIl6cKzmqY7oj_hcLYHXA1X5SXqshWeglCcB4ZUitmBmu8zMmOrbo3Dr_pWZfF4c$<
https://urldefense.com/v3/__http:/lists.ambermd.org/mailman/listinfo/amber__;!!KGKeukY!wOc8Qyq17CCFYVTtfvM8yYGX6ITnIl6cKzmqY7oj_hcLYHXA1X5SXqshWeglCcB4ZUitmBmu8zMmOrbo3Dr_pWZfF4c$>
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Sat Jun 08 2024 - 19:00:02 PDT