Re: [AMBER] Perturbing charges

From: Delwakkada Liyanage, Senal Dinuka <sdd313.msstate.edu>
Date: Mon, 21 Sep 2020 18:34:56 +0000

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
Custom Search