Hi Satyajit,
I wanna share my suggestions for WHAM and if it’s not correct, anyone in the mailing list can correct me.
First the unit for WHAM should be kcal/(mol*A^2). Second, the tolerance value I’d use is 0.0001.
Best,
Lin
> On Sep 18, 2019, at 4:06 AM, SATYAJIT KHATUA <satyajitkhatua09.gmail.com> wrote:
>
> Dear AMBER users,
> Recently, I have done umbrella sampling for
> PMF calculation using PLUMED plugin coupled with AMBER for a ligand-bound
> protein. The conventional unit of spring constant is [energy units of the
> program] / [units of the CV]^2. So after completing my umbrella sampling
> with let's say 34 bins, we go for WHAM to generate PMF. For PLUMED, units
> of the spring constant are kJ/mol/[nm]^2. Now for the calculation in WHAM,
> I have to transfer units to kcal/mol/[nm]^2. For my case, I have taken the
> spring constant of 300 kJ/mol/nm^2 and after conversion, the value is like
> 71.702 kcal/mol/nm^2. Although the WHAM program uses 0.5k(x-x0)^2, so we
> have to multiply the original force constant by 2. So, the final value of
> the spring constant is 143.404 kcal/mol/nm^2. After running WHAM, I am
> getting the PMF values which are continuously increasing. So, am I doing
> anything wrong with the unit conversion step? or there is something wrong
> about tolerance value from my side (taken as 0.01) in WHAM?? Please help me
> out.
>
> thanks in advance,
> Satyajit Khatua,
> Dept. of Chemistry,
> IIT Kharagpur,
> Kharagpur-721302.
> _______________________________________________
> AMBER mailing list
> AMBER.ambermd.org
> https://urldefense.proofpoint.com/v2/url?u=http-3A__lists.ambermd.org_mailman_listinfo_amber&d=DwICAg&c=nE__W8dFE-shTxStwXtp0A&r=3FEhv0ODXmqNZDIrYntxr0EDkmWDtT3vpmO7cMu9ps0&m=1A138UfkIICh6kBWaTlngtHznEScs-3d-zIycUDbL5I&s=2g1oDDpxsu3RAU3FFWtweBCHnCRUIFeyxsZ6zj0Plo8&e=
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Wed Sep 18 2019 - 09:30:01 PDT
