Re: [gmx-users] g_lie and ligand only simulation

2013-11-06 Thread Kavyashree M
Thank you.. On Wed, Nov 6, 2013 at 7:39 PM, Justin Lemkul wrote: > > > On 11/6/13 5:47 AM, Kavyashree M wrote: > >> Dear users, >> >> Sorry for repeating the same question. I just wanted to know >> whether is it ok if I have rlist > rcoulomb in ligand-water and >> prot-lig-water rerun md (wi

Re: [gmx-users] g_lie and ligand only simulation

2013-11-06 Thread Justin Lemkul
On 11/6/13 5:47 AM, Kavyashree M wrote: Dear users, Sorry for repeating the same question. I just wanted to know whether is it ok if I have rlist > rcoulomb in ligand-water and prot-lig-water rerun md (with RF-0) while having rlist = rcoulomb in the original simulation using PME? The energi

Re: [gmx-users] g_lie and ligand only simulation

2013-11-06 Thread Kavyashree M
Dear users, Sorry for repeating the same question. I just wanted to know whether is it ok if I have rlist > rcoulomb in ligand-water and prot-lig-water rerun md (with RF-0) while having rlist = rcoulomb in the original simulation using PME? Thank you Regards Kavya On Wed, Nov 6, 2013 at 11

Re: [gmx-users] g_lie and ligand only simulation

2013-11-05 Thread Kavyashree M
Dear users, When the simulation was carried out with PME rcoulomb was set equal to rlist. But when I need to to ligand-water simulation without PME (with RF-0) then it requires rlist greater by 0.1-0.3 than rcoulomb. So if I rerun protein-ligand-water simulation there could be more differences in

Re: [gmx-users] g_lie useage

2013-11-04 Thread Justin Lemkul
On 11/4/13 1:52 AM, Saman Shahriyari wrote: dear gmx users i am trying to use g_lie tool regarding the fact i am a newbie to this. so i searched through the web for an appropriate protocol. although there are lots of staffs discussing the theory, but i couldn't find any thing describing the

[gmx-users] g_lie useage

2013-11-03 Thread Saman Shahriyari
dear gmx users i am trying to use g_lie tool regarding the fact i am a newbie to this. so i searched through the web for an appropriate protocol. although there are lots of staffs discussing the theory, but i couldn't find any thing describing the steps techniqally and with some details. genera

Re: [gmx-users] g_lie and ligand only simulation

2013-11-02 Thread Kavyashree M
Ok thank you. I thought it was for protein-ligand-water that needs to be rerun without PME. Thanks Regards Kavya On Sat, Nov 2, 2013 at 9:45 PM, Justin Lemkul wrote: > > > On 11/2/13 12:14 PM, Kavyashree M wrote: > >> Sir, >> >> Thank you. Should the ligand-water MD be done without PME? >> >>

Re: [gmx-users] g_lie and ligand only simulation

2013-11-02 Thread Justin Lemkul
On 11/2/13 12:14 PM, Kavyashree M wrote: Sir, Thank you. Should the ligand-water MD be done without PME? I already answered this. Please read my previous reply again. -Justin Thank you Regards Kavya On Sat, Nov 2, 2013 at 9:13 PM, Justin Lemkul wrote: On 11/2/13 1:22 AM, Kavyash

Re: [gmx-users] g_lie and ligand only simulation

2013-11-02 Thread Kavyashree M
Sir, Thank you. Should the ligand-water MD be done without PME? Thank you Regards Kavya On Sat, Nov 2, 2013 at 9:13 PM, Justin Lemkul wrote: > > > On 11/2/13 1:22 AM, Kavyashree M wrote: > >> Dear Users, >> >> Its mentioned in the list that it would be >> wrong to use g_lie on a simulation w

Re: [gmx-users] g_lie and ligand only simulation

2013-11-02 Thread Justin Lemkul
On 11/2/13 1:22 AM, Kavyashree M wrote: Dear Users, Its mentioned in the list that it would be wrong to use g_lie on a simulation which uses PME. So kindly suggest any other way available to get the free energy of ligand binding other using g_lie? The original simulation should be done wit

Re: [gmx-users] g_lie and ligand only simulation

2013-11-02 Thread Justin Lemkul
On 11/2/13 1:00 AM, Kavyashree M wrote: Dear Gromacs users, I have a protein-ligand in water simulation (Gmx 4.5.3), for calculating free energy of ligand binding, a separate simulation of ligand in water simulation is required (which I read from the list). The question is the protein-ligand i

[gmx-users] g_lie and ligand only simulation

2013-11-01 Thread Kavyashree M
Dear Users, Its mentioned in the list that it would be wrong to use g_lie on a simulation which uses PME. So kindly suggest any other way available to get the free energy of ligand binding other using g_lie? Thank you Regards kavya -- gmx-users mailing listgmx-users@gromacs.org http://lists

[gmx-users] g_lie and ligand only simulation

2013-11-01 Thread Kavyashree M
Dear Gromacs users, I have a protein-ligand in water simulation (Gmx 4.5.3), for calculating free energy of ligand binding, a separate simulation of ligand in water simulation is required (which I read from the list). The question is the protein-ligand is simulated as a dimeric system so is it nec

[gmx-users] g_lie

2013-07-07 Thread Mahboobeh Eslami
hi must all mdp files be similar for g_lie program? shouldn't i use of PME in all mdp files? what edr file must i use as input file for g_lie? many thanks / -- gmx-users mailing listgmx-users@gromacs.org http://lists.gromacs.org/mailman/listinfo/gmx-users * Please s

Re: [gmx-users] g_lie reproducibility (g_lie is lying)

2012-06-13 Thread Mark Abraham
On 14/06/2012 3:56 PM, Tom Dupree wrote: Greetings all, I can't manually reproduce g_lie results. After raging at excel for a while I think I have found a bug. Here is my first time point, Reported by g_lie to be 35.0073 > From energy file Lj_complex =-130.762 Coul_complex = -286.746 My constant

[gmx-users] g_lie reproducibility (g_lie is lying)

2012-06-13 Thread Tom Dupree
Greetings all, I can't manually reproduce g_lie results. After raging at excel for a while I think I have found a bug. Here is my first time point, Reported by g_lie to be 35.0073 >From energy file Lj_complex =-130.762 Coul_complex = -286.746 My constants specified to g_lie Clj = Alpha = lj_const

[gmx-users] g_lie question

2011-11-18 Thread EGY
I am trying to figure out how to use the g_lie tool in gromacs. I have the dGbind value for the ligand-solvent interactions and for the protein-ligand complex. However, I cannot find in the manual what units these values are. I get the following: DGbind = 1953867.080 (0.664) and DGbind = 185147

[gmx-users] g_lie (again..)

2011-09-12 Thread גדעון לפידות
Hi all, Could anyone please refer to me to a paper or other resource that explains in detail (specifically technical) on how to preform g_lie. I know the matter has been brought up here quit often lately but looking through the mailing list archives and reference papaers (Aquist et al. and similar

[gmx-users] g_lie

2011-09-11 Thread gal.fradin
Greetings to all! I have a question for the g_lie function users... Is there anyone who used this function successfully including the -Elj and -Eqq options??? When I'm trying to use this function with the above mentioned options I'm getting Zeros for both Elj and Eqq or any other value w

[gmx-users] g_lie

2011-08-29 Thread pawan raghav
g_lie calculates the deltaG= alpha(p - w) + beta(p - w) equation. But when calculating DGbind from g_lie it will ask for only one term at a time i.e. -Elj for either p or w and -Cqq for either p or w. My question is when using only p for -Elj and Vl-s el>p for -Cqq, how the LJ and EL terms calcul

Re: [gmx-users] g_lie query

2010-12-28 Thread Justin A. Lemkul
Anirban Ghosh wrote: Thanks Justin. Yes I did use PME during the simulations. But still when I process the .edr file with g_energy, it gives both the options...LJ-1-4 and LJ(SR). So if I am not wrong, I should be using the LJ(SR) value as the Elj in the LIE calculation. Right? I'll have

Re: [gmx-users] g_lie query

2010-12-28 Thread Anirban Ghosh
Thanks Justin. Yes I did use PME during the simulations. But still when I process the .edr file with g_energy, it gives both the options...LJ-1-4 and LJ(SR). So if I am not wrong, I should be using the LJ(SR) value as the Elj in the LIE calculation. Right? Thanks again. Anirban On Tue, Dec 28,

Re: [gmx-users] g_lie query

2010-12-28 Thread Justin A. Lemkul
Anirban Ghosh wrote: Thanks a lot Justin for the reply !!! While calculating the Elj and Eqq values should we consider the short-ranged LJ (LJ-SR) or LJ-14, the two components that are present in 1-4 interactions are intramolecular, thus should not be relevant to the LIE calculation. t

Re: [gmx-users] g_lie query

2010-12-28 Thread Anirban Ghosh
Thanks a lot Justin for the reply !!! While calculating the Elj and Eqq values should we consider the short-ranged LJ (LJ-SR) or LJ-14, the two components that are present in the .edr file? And for Coulomb also? Which one should we consider? Thanks again. Anirban On Tue, Dec 28, 2010 at 11:03

Re: [gmx-users] g_lie query

2010-12-28 Thread Anirban Ghosh
Thanks a lot Justin !!! --Anirban On Tue, Dec 28, 2010 at 11:03 PM, Justin A. Lemkul wrote: > > > Anirban Ghosh wrote: > >> Thanks a lot Justin for the reply. >> Yes I am going through all the relevant literature on LIE. >> Actually the lie.xvg file contains the same value of -25.4 for all the

Re: [gmx-users] g_lie query

2010-12-28 Thread Justin A. Lemkul
Anirban Ghosh wrote: Thanks a lot Justin for the reply. Yes I am going through all the relevant literature on LIE. Actually the lie.xvg file contains the same value of -25.4 for all the frames. So I am getting a straight line plot. Why is this happening? Am I missing out something? If the

Re: [gmx-users] g_lie query

2010-12-28 Thread Anirban Ghosh
Thanks a lot Justin for the reply. Yes I am going through all the relevant literature on LIE. Actually the lie.xvg file contains the same value of -25.4 for all the frames. So I am getting a straight line plot. Why is this happening? Am I missing out something? Thanks a lot again. Anirban On Tu

Re: [gmx-users] g_lie query

2010-12-28 Thread Justin A. Lemkul
Anirban Ghosh wrote: Thanks a lot Justin for the reply. So I ran a simulation with my ligand in water for 1 ns and using g_energy I calculated the LG-14 and Coulomb-14 values from the .edr file. I supplied the average of these two values as my Elj and Eqq to g_lie and I got the DGbind as -25

Re: [gmx-users] g_lie query

2010-12-28 Thread Anirban Ghosh
Thanks a lot Justin for the reply. So I ran a simulation with my ligand in water for 1 ns and using g_energy I calculated the LG-14 and Coulomb-14 values from the .edr file. I supplied the average of these two values as my Elj and Eqq to g_lie and I got the DGbind as -25.4. Is this the correct way

Re: [gmx-users] g_lie query

2010-12-27 Thread Justin A. Lemkul
Anirban Ghosh wrote: Thanks Justin for the reply. I have through the threads about g_lie, but cannot understand how to get the values for Elj and Eqq for a particular ligand. Like in my case for a system consisting of a beta2AR protein + dopamine (ligand) + POPC + water, what should be the v

Re: [gmx-users] g_lie query

2010-12-27 Thread Anirban Ghosh
Thanks Justin for the reply. I have through the threads about g_lie, but cannot understand how to get the values for Elj and Eqq for a particular ligand. Like in my case for a system consisting of a beta2AR protein + dopamine (ligand) + POPC + water, what should be the values for Elj and Eqq? Than

Re: [gmx-users] g_lie query

2010-07-17 Thread Justin A. Lemkul
Anirban Ghosh wrote: Hi ALL, I have run a protein + ligand (dopamine) simulation. Now I want to calculate the free energy of binding using g_lie. But g_lie asks for two values: Elj and Eqq. How or from where can I get these values for my ligand? Also, do I need to run a simulation with only

[gmx-users] g_lie query

2010-07-17 Thread Anirban Ghosh
Hi ALL, I have run a protein + ligand (dopamine) simulation. Now I want to calculate the free energy of binding using g_lie. But g_lie asks for two values: Elj and Eqq. How or from where can I get these values for my ligand? Also, do I need to run a simulation with only the ligand? And, is there a

[gmx-users] g_lie query

2010-07-16 Thread Anirban Ghosh
Hi ALL, I have run a protein + ligand (dopamine) simulation. Now I want to calculate the free energy of binding using g_lie. But g_lie asks for two values: Elj and Eqq. How or from where can I get these values for my ligand? Also, do I need to run a simulation with only the ligand? And, is there a

[gmx-users] g_lie query

2010-07-16 Thread Anirban Ghosh
Hi ALL, I have run a protein + ligand (dopamine) simulation. Now I want to calculate the free energy of binding using g_lie. But g_lie asks for two values: Elj and Eqq. How or from where can I get these values for my ligand? Also, do I need to run a simulation with only the ligand? And, is there a

Re: [gmx-users] g_lie + running MDS for a ligand only

2008-10-15 Thread Justin A. Lemkul
vivek sharma wrote: Hi All, I want to use g_lie function of gromacs to calculate the binding free energy og a ligand with the protein. For the same I need to run two MDS, one for ligand only and one for ligand+protein complex (please correct if I am wrong). What are the checks that I should

[gmx-users] g_lie + running MDS for a ligand only

2008-10-14 Thread vivek sharma
Hi All, I want to use g_lie function of gromacs to calculate the binding free energy og a ligand with the protein. For the same I need to run two MDS, one for ligand only and one for ligand+protein complex (please correct if I am wrong). What are the checks that I should take care in this case, e.

Re: [gmx-users] g_lie problem

2007-10-01 Thread Mark Abraham
Blaise Mathias-Costa wrote: hello all, In my ligand enzyme simulation, I have given *energygrps: protein non protein. *The reason is that when I gave energygrps* *give as : Protein DRG SOL CL, the mdrun stopped at 560ps saying that the system is physically unrealistic! The identity of energ

[gmx-users] g_lie problem

2007-10-01 Thread Blaise Mathias-Costa
hello all,In my ligand enzyme simulation, I have given energygrps: protein non protein.  The reason is that when I gave energygrps give as : Protein DRG SOL CL, the mdrun stopped at 560ps saying that the system is physically unrealistic! Now I want to analyze the g_lie for the ligand, is there any

[gmx-users] g_lie restraints, neutral ??

2006-09-13 Thread Mikko Hellgren
Hi MD people! Some questions about g_lie. My purpose is to score different ligands and get a rough estimation of the free energy of binding between free and bound state of the ligands and protein of interest. 1. When I run my ligands in water without the protein, should I add counterions (Cl