Hello,
thank you for an advice. This might be relly the way, however, I do not
know how to perform this step. I dont know how to couple a system between
free and restrained state (how exatly to do it, as in the manual there is
not a word about this).
Please, can you give me any suggestions?
> I just realise that in both sander and pmemd, NMR restraints (nmropt=1)
> are scaled too as per the link in your original email. This may be the
> better option to pursue.
>
> ________________________________________
> From: Loeffler, Hannes (STFC,DL,SC)
> Sent: 25 March 2015 19:12
> To: AMBER Mailing List
> Subject: RE: [AMBER] Free energy calculation: Thermodynamic integration
> with use of restraints
>
> So, what you want to do is an "absolute" free energy calculation with the
> TI code. You also want to only couple a set of restraints to lambda. If
> I understand the Boresch' method correctly you would need to define six
> such restraints: one bond, two angles, three dihedrals. Right?
>
> Sorry, I'm not too well at the moment so I may get this wrong but wouldn't
> it be possible to define those restraints in terms of existing atoms and
> introduce additional bonds, angles and dihedrals between them. In the
> input file you would then define the TI region by including those atoms
> only while the softcore mask is empty. Maybe that's a plan but I would
> have to rethink it. The only other idea I can come up with is to create
> multiple topology files with those additional terms pre-scaled by the
> current lambda but then you would need to do a EXP/BAR estimation for the
> associated free energy as I can't see at the moment how this could be fit
> into the TI formalism.
>
> Cheers,
> Hannes.
>
> ________________________________________
> From: vladimir.palivec.marge.uochb.cas.cz
> [vladimir.palivec.marge.uochb.cas.cz]
> Sent: 25 March 2015 17:55
> To: amber.ambermd.org
> Subject: [AMBER] Free energy calculation: Thermodynamic integration with
> use of restraints
>
> Hello,
>
> this question is similar to this one:
>
> http://archive.ambermd.org/201205/0093.html
>
> however, I think I quite did not get it.
>
> I would like to do a thermodynamic integration simulations to calculate dG
> of binding of a ligand to a protein. After search through a literature I
> decided to do it in multiple steps:
>
> 1. constrain a ligand in a certain position using 6 restraints
> - here i need to calculate dG of proposing restraints
>
> 2. turn off el. i.
>
> 3. using soft core, turn off vww interactions
>
> 4. analytically remove restraints
>
>
> I am not sure how to perform the step 1. I thought it is possible to use
> pmemd with single topology and gradually propose restraints and thus get
> the dG (with use of nmropt=1 and DIANG file). However, I have no found any
> way to do it - either in manual or forum. Please, can somebody give me a
> hint how this calculation should be performed? How to set up the
> calculation and how to obtain energy of proposing restraints?
>
> Is it even possible to do it with one topology file?
>
> Thank you very much for any help!
>
> Vlada
>
>
>
>
>
>
>
> _______________________________________________
> AMBER mailing list
> AMBER.ambermd.org
> http://lists.ambermd.org/mailman/listinfo/amber
>
> _______________________________________________
> AMBER mailing list
> AMBER.ambermd.org
> http://lists.ambermd.org/mailman/listinfo/amber
>
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Thu Mar 26 2015 - 02:00:02 PDT
