Thank you Dr. DasGupta for reaching out. I went through your suggestion and as for the clarification purpose please let me know if I'm correct.
For the removal of electrostatic step, I will be needing a prmtop which includes my ligand and its copy (under 2 residues names). Then I can use 'crgmask' and set charges to 0 in one of them. For this set of calculation with sufficient lambda windows 'scmask' will be set to blank (softcore potential is not used).
For the removal of vdW step, the same prmtop will be used while turning on 'scmask1' and 'scmask2' for the two residues respectively.
Looking forward to hearing from you soon.
Senal Dinuka
________________________________
From: Debarati DasGupta <debarati_dasgupta.hotmail.com>
Sent: Sunday, September 20, 2020 3:00 PM
To: AMBER Mailing List <amber.ambermd.org>
Subject: Re: [AMBER] Perturbing charges
Hello Senal Dinuka
For absolute free energy calculations I do not think you need to use the tiMerge option. The timasks are important and the scmasks.
I have not used tiMerge in disappearing a molecule type TI runs. I may be wrong and hopefully experts on the forum can correct me on this.
Debarati
Sent from Mail<
https://secure-web.cisco.com/1-Hoi5DsWDkDZyoA6IeGIULH4Lyisdc_ZHJtK917zpb2T_WIgAcek8QSd6n2uMi_O6wEJAxzCGH5I-lL6Y6DrGqLVFYWfqG1tU0_ESt6HB_xXBcOfI9oM3pnqtftCjfvGv9slPF_7SDHjGM1HAo7GdrPyWIZbDNnWlrKI_LQsz7tHOg4EY2IcqCCTLox2DzTLv9fairnqmSAihzi7RB1bL3IGSRakrDOykYQUI9nmZyyHMVUpNwMwd4tVcJeM4j7C_CkSh1zc1Drg5YZpYBA1FKEQaJUJpL8J1R1QmhspabcXrB1WiRHxE1EFbvs53UoHpm-XlvGRXJb2QqVYhZYRgw/https%3A%2F%2Fgo.microsoft.com%2Ffwlink%2F%3FLinkId%3D550986> for Windows 10
From: David A Case<mailto:david.case.rutgers.edu>
Sent: 20 September 2020 09:02
To: AMBER Mailing List<mailto:amber.ambermd.org>
Subject: Re: [AMBER] Perturbing charges
On Fri, Sep 18, 2020, Delwakkada Liyanage, Senal Dinuka wrote:
>Thank you very much for your suggestion, Dr. Case. I'm struggling with
>how to implement this approach with pmemd. I followed the tutorial
>and they have considered two different molecules combine in a single
>prmtop. But in my case since I'm trying to disappear the ligand and I'm
>lost as to how to form a combined prmtop with the part which have vdW
>while the other doesn't.
It's unclear what tutorial you are following here. The pmemd approach
is described in Section 23.1.8 of the manual, and uses the tiMerge
facility in parmed.
I'm not sure if/how crgmask works in pmemd. Others can chime in here.
As I noted before, the "change" action in parmed, which I think would be
used after tiMege, might help. But there are plenty of people on the
list who have done this more than I have, so I hope they will contribute
to the discussion.
....dac
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://secure-web.cisco.com/160ZjY-CF6_VveEt5WRJy0RPrfd-tG6sAYv2B36CvARSBhZc_T7IHE8jY15-3822VZOjHDIwwPhCrWKP31Fb4TxJOE3FIaMcwG55vCueShSUR80RER-ncuVOm65bMDewGi3QdwQLk2HxtiB6CMsFYvmOST_Bj9VU8GohDVpJA6am0LoEM6vU_x-OcTT7z_Mme8qAdPpx5WpWIqow7DsQZWJs-ogEQAsi4i-fKSV9iHfyK31KtGzvT-oKOK-eOzjLu6lLj_JEu20jgYVpsMfEvGb3rqNvMXfi8nDoL5yezU_2iEoqMz2DyqIdOk6t5jwydbpnGSDwj6WZrkYsjs2w2vA/http%3A%2F%2Flists.ambermd.org%2Fmailman%2Flistinfo%2Famber
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://secure-web.cisco.com/160ZjY-CF6_VveEt5WRJy0RPrfd-tG6sAYv2B36CvARSBhZc_T7IHE8jY15-3822VZOjHDIwwPhCrWKP31Fb4TxJOE3FIaMcwG55vCueShSUR80RER-ncuVOm65bMDewGi3QdwQLk2HxtiB6CMsFYvmOST_Bj9VU8GohDVpJA6am0LoEM6vU_x-OcTT7z_Mme8qAdPpx5WpWIqow7DsQZWJs-ogEQAsi4i-fKSV9iHfyK31KtGzvT-oKOK-eOzjLu6lLj_JEu20jgYVpsMfEvGb3rqNvMXfi8nDoL5yezU_2iEoqMz2DyqIdOk6t5jwydbpnGSDwj6WZrkYsjs2w2vA/http%3A%2F%2Flists.ambermd.org%2Fmailman%2Flistinfo%2Famber
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Mon Sep 21 2020 - 12:00:02 PDT
