Hello Matthew,
Let’s imagine your Lys-8 pKa is ~ 12. This means you need to go at least 2 pH units up to get a better fit. Therefore, 2 to 14 is likely a safe range for you to consider.
I hope this help,
All the best,
Vinícius Wilian D. Cruzeiro, PhD
Postdoctoral Researcher
Department of Chemistry
Stanford University
Twitter: @vwcruzeiro
On Feb 8, 2023, at 10:36 PM, Matthew Guberman-Pfeffer <matthew.guberman-pfeffer.uconn.edu> wrote:
Hi Vinicius,
Thanks for your quick reply. Your suggestion was very informative! The output showed me that the protonated fraction for Lys-8 and AS4-114 is always >0.99 and <0.005, respectively. So, the pKas are strange because only extremely rare transitions are being sampled.
What are the lowest and highest extremes in pH you’d recommend testing? Is it 2 to 12 (the pKas range, I think can be measured in aqueous solutions)?
Best,
Matthew
From: Vinicius Wilian Cruzeiro <vwcruz.stanford.edu>
Date: Thursday, February 9, 2023 at 12:58 AM
To: amber.ambermd.org <amber.ambermd.org>, Matthew Guberman-Pfeffer <matthew.guberman-pfeffer.uconn.edu>
Subject: Re: Suspicious pKas from CpHMD
*Message sent from a system outside of UConn.*
Hello Matthew,
You should take a look at the fractions of protonated species as a function of pH for these residues. This will give you an idea on why you are getting these fitted values.
One possibility is that you are having a conformational change and then your fractions of protonated species vs pH will not have a nice sigmoidal shape, what could explain the fitted values you got.
Also, I would not exclude the possibility that you need to add more pH values, as you pointed out. Notice that the plots I am suggesting will tell you if this is what you need to do.
I hope this helps,
Vinícius Wilian D. Cruzeiro, PhD
Postdoctoral Researcher
Department of Chemistry
Stanford University
Twitter: @vwcruzeiro
________________________________
From: Matthew Guberman-Pfeffer via AMBER <amber.ambermd.org>
Sent: Wednesday, February 8, 2023 9:08 PM
To: amber.ambermd.org <amber.ambermd.org>
Subject: [AMBER] Suspicious pKas from CpHMD
Dear AMBER Community,
I’m getting some really strange pKas for CpHMD simulations performed at pHs from 5 to 10. For example,
LYS 8: pKa = 564.07106 (+- 9292.89711 ) ; Hill coefficient = 0.004 (+- 0.060 )
AS4 114: pKa = -52.65171 (+- 144.03991 ) ; Hill coefficient = 0.046 (+- 0.112 )
Each simulation was propagated for 40 ns. My CpH-relevant settings were icnstph=2, ntcnstph=100, ntrelax=100, and solvph = 5, 6, 7, 8, 9, or 10.
Is there anything besides running the simulations for longer that I should try? Should I expand the range of considered pHs? What else can I try to get more reasonable pKas for these residues? Or, what should I inspect to see if something is going wrong?
Best,
Matthew
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber<
https://nam10.safelinks.protection.outlook.com/?url=http%3A%2F%2Flists.ambermd.org%2Fmailman%2Flistinfo%2Famber&data=05%7C01%7Cmatthew.guberman-pfeffer%40uconn.edu%7Cd0f7dd35d34843a3f8a208db0a629785%7C17f1a87e2a254eaab9df9d439034b080%7C0%7C0%7C638115190802844122%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C&sdata=yunrTinsYZZqT9Shpad0gmCSHNA7WgvJ0UZURaFhHos%3D&reserved=0>
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Wed Feb 08 2023 - 23:30:02 PST
