Hi gmx users,<br />
<br />
I used g_lie command to compute free energy.<br />
The command and the output is as follow<br />
<br />
[shikhar@venus 30078-dyn10ns]$ g_lie -f edr.edr -ligand UNK -o lie10ns.xvg<br />
<br />
-ligand string UNK Name of the ligand in the energy file<br />
<br />
Opened edr.edr as single precision energy file<br />
Using the following energy terms:<br />
LJ: LJ-SR:Protein-UNK LJ-LR:Protein-UNK LJ-14:Protein-UNK LJ-SR:SOL-UNK LJ-LR:SOL-UNK LJ-14:SOL-UNK LJ-SR:UNK-rest LJ-LR:UNK-rest LJ-14:UNK-rest<br />
Coul: Coul-SR:Protein-UNK Coul-14:Protein-UNK Coul-SR:SOL-UNK Coul-14:SOL-UNK Coul-SR:UNK-rest Coul-14:UNK-rest<br />
Last frame read 50000 time 10000.000<br />
DGbind = -88.393 (5.098)<br />
<br />
gcq#306: "Miggida-Miggida-Miggida-Mac" (Kriss Kross)<br />
<br />
The value is coming DGbind = -88.393 (5.098).<br />
Whether this is correct method to get the free energy of binding of a ligand.<br />
Any other command or method is there to get the free energy of binding of ligands or can we do the MM/PBSA and MM/GBSA analysis from the gromacs dynamics.<br />
<br><br>------------------<br />
Shikhar Gupta<br />
Senior Research Fellow<br />
Pharmacoinformatics Department<br />
Block- A (Room No.- 208)<br />
National Institute of Pharmaceutical Education & Research( NIPER )<br />
Sec- 67, S.A.S Nagar<br />
Mohali, Punjab (India)<br />
Web-Site: www.niper.ac.in<br />
PIN- 160062<br />
Email:shik_sun@rediffmail.com,shiksun@gmail.com<br><Table border=0 Width=644 Height=57 cellspacing=0 cellpadding=0 style="font-family:Verdana;font-size:11px;line-height:15px;"><TR><td><A HREF="http://sigads.rediff.com/RealMedia/ads/click_nx.ads/www.rediffmail.com/signatureline.htm@Middle?" target="_blank"><IMG SRC="http://sigads.rediff.com/RealMedia/ads/adstream_nx.ads/www.rediffmail.com/signatureline.htm@Middle"></A></td></TR></Table>