Hi Nitin,
It seems like the fraction of protonation from the concatenated files is very low. To confirm if there’s a trend of decreased protonation:
For the individual chunks of simulations, you can run cphstats for the original cpout files. If you know the length of each chunk, the total fraction of protonation can be calculated as a weighted average. No concatenation needed.
If it doesn’t look like an expected trend, maybe it’s due to the concatenation. You can compared with the weighted average and see if you’re concatenating in the way you actually want.
From: Nitin Kulhar <bo18resch11002.iith.ac.in>
Date: Tuesday, August 6, 2024 at 9:37 AM
To: He, Amy <he.1768.buckeyemail.osu.edu>, AMBER Mailing List <amber.ambermd.org>
Subject: Re: [AMBER] cphstats processing of pH-REMD 'chunks'
Hi Amy Thanks for the email. I apologise for any confusion caused by the term 'chunks'. To clarify, I was NOT referring to the various interval-based averages namely, cumulative, running, and chunk, that can be obtained from cphstats.
Hi Amy
Thanks for the email.
I apologise for any confusion
caused by the term 'chunks'.
To clarify, I was NOT referring to
the various interval-based
averages namely, cumulative,
running, and chunk, that can be
obtained from cphstats.
I meant 'chunk' in the context of
Jason Swail's tutorial on pH-REMD:
http://jswails.wikidot.com/ph-remd#toc10<
https://urldefense.com/v3/__http:/jswails.wikidot.com/ph-remd*toc10__;Iw!!KGKeukY!0QPqWCdXnVdLe8rKHuhDQJUoHGhMynPyu-b40QIv63g2Os1e-DpK5_jKEPGP4bL9_uzoA0pBmpXdBHrl_s3gAVsi_QK48jQAXHEKQHKF$>
To quote from the tutorial:
If you have run 10 ns of simulation
in 1 ns chunks, you will need to
run this program 10 times!
PFA two titration curves: one obtained from
by processing just the first segment, the other
obtained by processing concatenated cpouts
from all segments.
Thanks and regards
Nitin Kulhar
On Mon, Aug 5, 2024 at 4:35 AM He, Amy <he.1768.buckeyemail.osu.edu<mailto:he.1768.buckeyemail.osu.edu>> wrote:
Hi Nitin,
Why do you think the results are wrong? Did it cumulate, or not cumulate, or not space in the way you want?
The manual will give you a detailed difference with examples for the Cumulative, running, and “chunk” options of cphstats.
From: Nitin Kulhar via AMBER <amber.ambermd.org<mailto:amber.ambermd.org>>
Date: Friday, August 2, 2024 at 9:38 PM
To: AMBER Mailing List <amber.ambermd.org<mailto:amber.ambermd.org>>
Subject: [AMBER] cphstats processing of pH-REMD 'chunks'
Dear all
I performed pH-REMD of 8
REPLICAS. 25 ns of REMD
were run in 1 ns SEGMENTS
(or chunks) as delineated in
the tutorial by Jason Swails
available at link <
https://urldefense.com/v3/__http://jswails.wikidot.com/ph-remd__;!!KGKeukY!yNeLxboggV6vrz65BXBmJdJtWYdmVZcrdIOV1wBGc78GZvpfjF_-WCNsdRPGRjBJUhH7F7YBvR3VSKHDhd3UMV574E0$<
https://urldefense.com/v3/__http:/jswails.wikidot.com/ph-remd__;!!KGKeukY!yNeLxboggV6vrz65BXBmJdJtWYdmVZcrdIOV1wBGc78GZvpfjF_-WCNsdRPGRjBJUhH7F7YBvR3VSKHDhd3UMV574E0$> >.
*The problem:*
The pKa values obtained by
processing the cpout files from
all SEGMENTS (or chunks)
were questionable, nay, wrong.
However, the pKa values obtained
by processing the cpout file from
the first SEGMENT alone were
correct.
*Salient details of the REMD:*
The inputs for each SEGMENT
were generated by the
genremdinputs.py script.
The cpin and starting coordinates
for REPLICA 'i' of SEGMENT 'j'
were respectively the cprestart
and restart files generated by
REPLICA 'i' during SEGMENT 'j-1'.
While the seed (ig) value
differed among REPLICAS
of the same SEGMENT, it
was the same for a given
REPLICA in all SEGMENTS.
(Is that genremdinputs.py's doing?)
*Salient details of processing cpouts:*
The cpout files from consecutive
SEGMENTS were concatenated
with the cat command. This was
done for each REPLICA separately.
The concatenated cpout files of ALL
REPLICAS were sorted according to
pH values by using cphstats program
with the --fix-remd option.
The pH-sorted cpout files from above
were processed with cphstats to obtain
pKa values and population data.
Please help me spot the mistake
in my logic or execution of the steps.
Thanks and regards
Nitin Kulhar
--
Disclaimer:- This footer text is to convey that this email is sent by one
of the users of IITH. So, do not mark it as SPAM.
_______________________________________________
AMBER mailing list
AMBER.ambermd.org<mailto:AMBER.ambermd.org>
https://urldefense.com/v3/__http://lists.ambermd.org/mailman/listinfo/amber__;!!KGKeukY!yNeLxboggV6vrz65BXBmJdJtWYdmVZcrdIOV1wBGc78GZvpfjF_-WCNsdRPGRjBJUhH7F7YBvR3VSKHDhd3UQjpcOUQ$<https://urldefense.com/v3/__http:/lists.ambermd.org/mailman/listinfo/amber__;!!KGKeukY!yNeLxboggV6vrz65BXBmJdJtWYdmVZcrdIOV1wBGc78GZvpfjF_-WCNsdRPGRjBJUhH7F7YBvR3VSKHDhd3UQjpcOUQ$>
Disclaimer:- This footer text is to convey that this email is sent by one of the users of IITH. So, do not mark it as SPAM.
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Tue Aug 06 2024 - 13:00:04 PDT
