Hi Josiah,
I've only done a few ABFE runs with Amber on neutral molecules, so I can't claim to be expert - interested to hear input from others on the list. Here is my protocol; values for T4 Lysozyme were sensible.
Complex+Ligand
- restraint leg: transform one identical copy of your molecule without restraints into a second copy that does have Boresch restraints. So you need to set up a prmtop with two identical copies of your ligand. Then collect DV/DL values as you normally would. Here is the TI portion of my input (ligands are residue 1 and 2 respectively, with restraints on residue 2 only):
nmropt=1,
gti_add_sc = 1,
gti_lam_sch = 1,
tishake=1, noshakemask=':1,2', ntc=2, ntf=1,
icfe = 1, clambda = 0.00922, scalpha = 0.2, scbeta = 50.0,
logdvdl = 1,
ifmbar=1,
mbar_states = 12,
mbar_lambda = 0.00922 0.04794 0.11505 0.20634 0.31608 0.43738 0.56262 0.68392 0.79366 0.88495 0.95206 0.99078 ,
ifsc = 0,
crgmask = ''
timask1 = ':1', timask2 = ':2',
scmask1 = '', scmask2 = ''
/
&ewald
/
&wt
type='END',
&end
DISANG=disang.RST
- decouple leg: I used the one-step protocol to decouple electrostatics and vdW in a single step, with the same Boresch restraints applied to ligand (residue 1), now you have a prmtop with only the single ligand of interest. Note, I have to modify disang.RST file since prmtops are different between restraint and decouple step. You could split this up to separately remove electrostatics first followed by vdW.
nmropt=1,
gti_add_sc = 1,
gti_lam_sch = 1,
tishake=1, noshakemask=':1', ntc=2, ntf=1,
icfe = 1, clambda = 0.00922, scalpha = 0.2, scbeta = 50.0,
logdvdl = 1,
ifmbar=1,
mbar_states = 12,
mbar_lambda = 0.00922 0.04794 0.11505 0.20634 0.31608 0.43738 0.56262 0.68392 0.79366 0.88495 0.95206 0.99078 ,
ifsc = 1,
crgmask = ''
timask1 = ':1', timask2 = '',
scmask1 = ':1', scmask2 = ''
/
&ewald
/
&wt
type='END',
&end
DISANG=disang.RST
Solvent+Ligand
-decouple leg: settings as with the complex decouple leg but no restraints applied
-restraint term: calculate analytically
All the best,
Callum
-----Original Message-----
From: Josiah Bones <josiah.bones.anu.edu.au>
Sent: Friday, June 18, 2021 2:30 AM
To: amber.ambermd.org
Subject: [AMBER] Absolute Binding Free Energy of Protein+Ligand Complex in AMBER20 - Boresch restraints, Smoothstep
This Message is from an External Sender. Do not click links or open attachments unless you trust the sender.
Hi there,
I asked a question about some parts of this elsewhere, but I realised I needed to simplify my previous question into something a bit more basic and address some issues in a single location that I can't find mentioned in any tutorials.
I am wondering how to successfully complete what I believe is a fairly standard ABFE simulation in the context of AMBER20 on GPU. The goal is to capture the absolute binding of a ligand to a protein from a high-quality crystal structure.
I believe I need 3x simulation phases +1 post simulation calculation to correctly capture the dG.
1. Complex+Ligand: Simulate turning on Boresch restraints.
* I have defined these restraints in a DISANG file based on a high quality crystal structure using general recommendations mentioned elsewhere, but can't find any clue as to the correct way to simulate this leg anywhere, and how the .mdin file/s should look.
* Is this done without softcore potentials?
* What do the lambda windows look like, and how should this change with reference to the restraint file?
* Does the restraint file change between windows?
* What is V0 and V1 in this context?
* What other key .mdin toggles need attention beyond the defaults? Is dvdl_norest significant here, despite being deprecated?
* I would be absolutely thrilled to get commentary in particularly for this step, as it is completely elusive to me, and I have no idea how to solve it. Comments in a recent 2020 paper of TI changes in Amber made me assume that this was now essentially automated in Step 2 below, but I am not clear on this.
2. Complex+Ligand: Simulate decoupling of the restrained ligand in the protein complex.
* This I intend to do with SC potentials, in particular, the new Smoothstep default settings as published recently.
* Is there special attention I need to give to the lambda scheduling given the context of my calculation, or are default settings likely reasonable?
* Given I am decoupling a whole ligand with no shared common atoms, am I safe to use SHAKE in this context? Or is SHAKE a big no-no the moment SC is involved, requiring me to use the noshakemask setting for the ligand?
* Are there key default .mdin toggles that need my attention in this situation?
1. Ligand+Water: Essentially the same as Step 2, although in just water+ions and without any restraints.
* Same generic questions as Step 2.
2. Solve the equation relating to the Boresch restraints as they apply to my specific system.
In the context of a charged ligand, I believe I could condense steps 2 and 3 into a single step where the ligand is "recoupled" in bulk water in a single leg of the cycle, and is held a decent distance from the protein with COM restraints. Would this require any further attention in capturing the correct dG (despite the likely requirement for increased sampling given the larger overall perturbation)?
Any feedback on this process in the context of AMBER20 would be most appreciated.
Thank you kindly,
Josiah
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
https://urldefense.proofpoint.com/v2/url?u=http-3A__lists.ambermd.org_mailman_listinfo_amber&d=DwICAg&c=ZbgFmJjg4pdtrnL2HUJUDw&r=HwsFjSfOtLupDR-NuCP430rdz1DD2LkJxNM3BsKSjrw&m=WwGH4_77-S6VPUjHRAHLjOgYytcMJb0gaGkVXFOXBQI&s=CMfUrFd4iGUAR0JDDc68h3KqRbR-IzfFPvHdeq1K45o&e=
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Fri Jun 18 2021 - 08:00:03 PDT