Re: [gmx-users] distance restraints in mdp
Please don't reply to the entire digest. Copy the relevant portions and use a proper subject line. On 3/8/16 8:19 PM, Sana Saeed wrote: dear Justin Lemkuli already run that without restraints, it has no problems. even i used the gro file from the simple 5ns MD but still have these same errors and warnings. when i set Lincs constraints to all-bonds , simulation finishes fine but by using alchemical_analysis.py for free energies it gives errors and warnings. (Warning: states 29 and 37 have the same energies on the dataset.) i used 42 lambdas and about 32 states are having same energies. and when i set constraints to h-bonds then LINCS warnings. i tried to use different options for pressure coupling also because i read somewhere it also affects simulation. i am quite confused where the problem is. is it possible to send you topology file and you have a look?Thanks in advance Now we're talking about half a million different things. When reporting a problem, please keep it simple. It's impossible for anyone on this mailing list to be inside your head and know what you've done. If you want to troubleshoot distance restraints, we can do that. But let's take things one at a time. If you want to share files, please upload a tarball of *all* input files: coordinates, topology(ies), any relevant .mdp file(s), and a run script that lists all of your commands, exactly in sequence as you typed them. Please do not send this information privately; share with the list so others can contribute (I'm frankly swamped lately and have little time to delve deeply into such issues but I'll try to make a little time). -Justin -- == Justin A. Lemkul, Ph.D. Ruth L. Kirschstein NRSA Postdoctoral Fellow Department of Pharmaceutical Sciences School of Pharmacy Health Sciences Facility II, Room 629 University of Maryland, Baltimore 20 Penn St. Baltimore, MD 21201 jalem...@outerbanks.umaryland.edu | (410) 706-7441 http://mackerell.umaryland.edu/~jalemkul == -- Gromacs Users mailing list * Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting! * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists * For (un)subscribe requests visit https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a mail to gmx-users-requ...@gromacs.org.
Re: [gmx-users] distance restraints in mdp
On 3/7/16 10:35 PM, Sana Saeed wrote: i am performing protein ligand binding free energy calculation MD. i have applied distance restraint to 1 pair of atoms. now when i am running md it gives me LINCS warnings. first i thought it is because of wrong geometry so i changed complex, now i am doing it for Lysozyme+JZ4. but still same warning. i am confused about the distance restraints to be introduced in mdp file. disre = simple/ensemble? i read in manual, but still confused. i tried both with "simple it gives me LINCS warning and with ensemble it says that this version of gromacs doesnt support it. i am using 5.1 gromacs. Then simplify further and make sure the system runs without restraints or anything else. Just normal MD. Start simple to avoid chasing the wrong problem. -Justin -- == Justin A. Lemkul, Ph.D. Ruth L. Kirschstein NRSA Postdoctoral Fellow Department of Pharmaceutical Sciences School of Pharmacy Health Sciences Facility II, Room 629 University of Maryland, Baltimore 20 Penn St. Baltimore, MD 21201 jalem...@outerbanks.umaryland.edu | (410) 706-7441 http://mackerell.umaryland.edu/~jalemkul == -- Gromacs Users mailing list * Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting! * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists * For (un)subscribe requests visit https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a mail to gmx-users-requ...@gromacs.org.
[gmx-users] distance restraints in mdp
i am performing protein ligand binding free energy calculation MD. i have applied distance restraint to 1 pair of atoms. now when i am running md it gives me LINCS warnings. first i thought it is because of wrong geometry so i changed complex, now i am doing it for Lysozyme+JZ4. but still same warning. i am confused about the distance restraints to be introduced in mdp file. disre = simple/ensemble? i read in manual, but still confused. i tried both with "simple it gives me LINCS warning and with ensemble it says that this version of gromacs doesnt support it. i am using 5.1 gromacs. Sana Saeed Khan, -- Gromacs Users mailing list * Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting! * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists * For (un)subscribe requests visit https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a mail to gmx-users-requ...@gromacs.org.