------------------
Shikhar Gupta
Senior Research Fellow
Pharmacoinformatics Department
Block- A (Room No.- 208)
National Institute of Pharmaceutical Education & Research( NIPER )
Sec- 67, S.A.S Nagar
Mohali, Punjab (India)
Web-Site: www.niper.ac.in
PIN- 160062
Email:shik_...@rediffmail.com,shik...@gmail.com
--- Begin Message ---
Hi gmx users,
I used g_lie command to compute free energy.
The command and the output is as follow
[shik...@venus 30078-dyn10ns]$ g_lie -f edr.edr -ligand UNK -o lie10ns.xvg
-ligand string UNK Name of the ligand in the energy file
Opened edr.edr as single precision energy file
Using the following energy terms:
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
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
Last frame read 50000 time 10000.000
DGbind = -88.393 (5.098)
gcq#306: "Miggida-Miggida-Miggida-Mac" (Kriss Kross)
The value is coming DGbind = -88.393 (5.098).
Whether this is correct method to get the free energy of binding of a ligand.
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.
------------------
Shikhar Gupta
Senior Research Fellow
Pharmacoinformatics Department
Block- A (Room No.- 208)
National Institute of Pharmaceutical Education & Research( NIPER )
Sec- 67, S.A.S Nagar
Mohali, Punjab (India)
Web-Site: www.niper.ac.in
PIN- 160062
Email:shik_...@rediffmail.com,shik...@gmail.com
--- End Message ---
--
gmx-users mailing list gmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/search before posting!
Please don't post (un)subscribe requests to the list. Use the
www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/mailing_lists/users.php
