Hello Zimin,
Through the command "cpinutil.py --describe" you are able to access the full list of reference energies currently in the Amber library. There are residues in the library (like AS4, for example) with reference energies in the same range as yours.
finddgref.py is sort of a "fail safe" tool because its predictions come directly from the constant pH output(s) file(s) produced in the simulations. In the serial mode, for example, finddgref.py would only stop when the fraction of protonated species for a given residue is sufficiently close to 50% at the pH you set. Please bare in mind that your reference energy depends not only on the charge distributions you set for each protonation state, but also on your solvent model and igb choices. Other than that, as long as your initial structure is properly equilibrated and your simulation parameters are correct in your mdin file (make sure your MD simulations are long enough to ensure the fractions you are obtaining are converged), you should be fine.
I hope this helps,
Best,
Vinícius Wilian D. Cruzeiro
PhD Candidate
Department of Chemistry, Physical Chemistry Division
University of Florida, United States
Voice: +1(352)846-1633
________________________________
From: Feng, Zimin <Feng.Zimin.hydro.qc.ca>
Sent: Wednesday, February 27, 2019 3:25 PM
To: AMBER Mailing List
Subject: [AMBER] is this reference energy dG ~34 kcal/mol reasonable for acetic acid?
Dear AMBER community,
I have a question about the value of dG_ref for the constant pH MD.
Basically I was using finddgref.py to get the reference energy (deprotonated state was set to 0 energy so this will be the energy difference as well) for acetic acid CH3COOH. I got a seemingly big value as 34 kcal/mol or ~1.5 eV.
I'm new to AMBER and I wonder if this number is reasonable or not? During the calculation no error messages was shown but the result seems a lot bigger than what the other titratable residues have.
Any comments will be appreciated and I thank you in advance!
Zimin
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
https://urldefense.proofpoint.com/v2/url?u=http-3A__lists.ambermd.org_mailman_listinfo_amber&d=DwICAg&c=sJ6xIWYx-zLMB3EPkvcnVg&r=vg8iTdivJL1PEwCvcH8P2DFF-Rtc9lAQvIqnaSWm1Pc&m=JWDk2g-gT9mHEv-lJFhNOplySIl5OZHnnT6g0wZYOhw&s=3PwvN72kRtzl5f80MrkBoMW3eK7k088gneTxAMA9Tw0&e=
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Wed Feb 27 2019 - 15:00:03 PST
