Hello,
I would like to use TI to obtain pKa estimates of functional groups in RNA molecules. Water and ions will be modeled explicitly and the systems will typically consist of ~30,000 atoms. Because the functional group will be neutral at ë = 0 and charged at ë = 1, the overall system will also be neutral and charged, respectively. From what I have read, this should be okay as long as the system is sufficiently large, however, this is based on literature/discussion from 10+ years ago. Several suggestions were provided in an older post, found at
http://archive.ambermd.org/200708/0288.html. Is there now a preferred way of handling these types of simulations? Is it best to just let the charge of the total system change or should I implement something to keep the system neutral?
Thanks,
Drew Veenis
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Mon Oct 01 2018 - 12:30:02 PDT
