Update: I think I found the answer to the question just mentioned. From this manual section. https://amberhub.chpc.utah.edu/lie/
So looks like one could define a region A and a region B, and calculate the interaction energy between A and B atoms, using the lie program. Did not know this was technically an option. Good to know now!
> On Sep 23, 2021, at 10:05 PM, Liao <liaojunzhuo.aliyun.com> wrote:
>
> Thanks for the response Dave,
>
> Never performed LIE before, I wonder, very practically, where does LIE actually get the vdw and elec energy values? I’m influenced by the explicit solvent binding energy example, in which the energies are obtained from the md output file that writes the different types of energy of the whole system for each snapshot. From my understanding, LIE reads only the intermolecular interaction energy of the ligand with its environment, doesn’t it? Where exactly is this number calculated, such that the water energy or protein energy is totally absent?
> I think this is what’s puzzling me the most, not the model or the theory.
>
>>> On Sep 23, 2021, at 8:52 PM, David A Case <dacase1.gmail.com> wrote:
>>>
>>> On Fri, Sep 24, 2021, Liao wrote:
>>>
>>> ...it seems that LIE is actually direct calculation of binding energy in
>>> explicit solvent, is this correct?
>>
>> No, this is not correct. LIE does not attempt to compute free energy
>> differences, but rather to compute a quantity that is sometimes (often?)
>> *correlated with* binding free energies.
>>
>> Note that the weighting factors (alpha, beta, gamma in Fig.4 in the
>> tutorial) are far from one, and are only approximately universal. The LIE
>> estimate is intuitive, and fairly straightforward to compute. So it has its
>> pros and cons, which is why it shows up workflows like FEW.
>>
>> ...hope this helps....dac
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Thu Sep 23 2021 - 20:30:03 PDT
