Re: [gmx-users] How to set up constraint distcance pulling
On 6/7/12 7:38 AM, xiaowu759 wrote: Dear gmxers, I am performing COM pulling using GMX4.0.7, to obtain the PMF as a function of distance between the COMs of two molecules. However, quite strange resuluts are generated, i.e., the distance can not be fixed at the value defined by pull_init1. The pull opinions are given below. Could you please tell me how to deal with this? Thanks a lot for any reply. How far away from the target value are the results in terms of COM distance? Chaofu Wu ;Pull opinions pull = constraint pull_geometry = distance ;pull_dim = Y Y Y The comment here has no effect; pull_dim defaults to Y Y Y. pull_group0 = r1 pull_group1 = r2 pull_nstxout = 0 pull_nstfout = 5 pull_init1 = 0.80 How does this pull_init1 value compare to the COM separation in the coordinate file? If the distance is not very close, the constraint algorithm may fail. -Justin -- Justin A. Lemkul, Ph.D. Research Scientist Department of Biochemistry Virginia Tech Blacksburg, VA jalemkul[at]vt.edu | (540) 231-9080 http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin -- gmx-users mailing listgmx-users@gromacs.org http://lists.gromacs.org/mailman/listinfo/gmx-users Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/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/Support/Mailing_Lists
Re: [gmx-users] How to set up constraint distcance pulling
Date: Thu, 07 Jun 2012 07:59:27 -0400 From: Justin A. Lemkul jalem...@vt.edu Subject: Re: [gmx-users] How to set up constraint distcance pulling To: Discussion list for GROMACS users gmx-users@gromacs.org Message-ID: 4fd0979f.4020...@vt.edu Content-Type: text/plain; charset=ISO-8859-1; format=flowed Thank you very much! On 6/7/12 7:38 AM, xiaowu759 wrote: Dear gmxers, I am performing COM pulling using GMX4.0.7, to obtain the PMF as a function of distance between the COMs of two molecules. However, quite strange resuluts are generated, i.e., the distance can not be fixed at the value defined by pull_init1. The pull opinions are given below. Could you please tell me how to deal with this? Thanks a lot for any reply. How far away from the target value are the results in terms of COM distance? A: very far, several hundreds vs. 0.8 nm Chaofu Wu ;Pull opinions pull = constraint pull_geometry = distance ;pull_dim = Y Y Y The comment here has no effect; pull_dim defaults to Y Y Y. A: Yes pull_group0 = r1 pull_group1 = r2 pull_nstxout = 0 pull_nstfout = 5 pull_init1 = 0.80 How does this pull_init1 value compare to the COM separation in the coordinate file? If the distance is not very close, the constraint algorithm may fail. A: The COM separation in the coordinate file is 0.567 nm, which should be not too larger than 0.80 defined by pull_init1. -Justin -- Justin A. Lemkul, Ph.D. Research Scientist Department of Biochemistry Virginia Tech Blacksburg, VA jalemkul[at]vt.edu | (540) 231-9080 http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin -- Department of Chemistry and Materials Science Hunan University of Humanities, Science and Technology, Loudi 417000, the People's Republic of China (P.R. China)-- gmx-users mailing listgmx-users@gromacs.org http://lists.gromacs.org/mailman/listinfo/gmx-users Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/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/Support/Mailing_Lists
Re: [gmx-users] How to set up constraint distcance pulling
On 6/7/12 9:14 AM, xiaowu759 wrote: Date: Thu, 07 Jun 2012 07:59:27 -0400 From: Justin A. Lemkul jalem...@vt.edu mailto:jalem...@vt.edu Subject: Re: [gmx-users] How to set up constraint distcance pulling To: Discussion list for GROMACS users gmx-users@gromacs.org mailto:gmx-users@gromacs.org Message-ID: 4fd0979f.4020...@vt.edu mailto:4fd0979f.4020...@vt.edu Content-Type: text/plain; charset=ISO-8859-1; format=flowed Thank you very much! On 6/7/12 7:38 AM, xiaowu759 wrote: Dear gmxers, I am performing COM pulling using GMX4.0.7, to obtain the PMF as a function of distance between the COMs of two molecules. However, quite strange resuluts are generated, i.e., the distance can not be fixed at the value defined by pull_init1. The pull opinions are given below. Could you please tell me how to deal with this? Thanks a lot for any reply. How far away from the target value are the results in terms of COM distance? A: very far, several hundreds vs. 0.8 nm A distance of several hundred nm would indicate to me that the system blew up. Chaofu Wu ;Pull opinions pull = constraint pull_geometry = distance ;pull_dim = Y Y Y The comment here has no effect; pull_dim defaults to Y Y Y. A: Yes pull_group0 = r1 pull_group1 = r2 pull_nstxout = 0 pull_nstfout = 5 pull_init1 = 0.80 How does this pull_init1 value compare to the COM separation in the coordinate file? If the distance is not very close, the constraint algorithm may fail. A: The COM separation in the coordinate file is 0.567 nm, which should be not too larger than 0.80 defined by pull_init1. That's a 29% difference, which is quite large. If constraints are not easily satisfied, the simulation can blow up. Your statement above makes that seem very likely. Your restraint distance should match the actual coordinates very closely, otherwise you will likely need to equilibrate better under e.g. an umbrella potential and relatively weak force constant to obtain better starting structures if you're intent on using a rigid constraint for the PMF calculation. -Justin -- Justin A. Lemkul, Ph.D. Research Scientist Department of Biochemistry Virginia Tech Blacksburg, VA jalemkul[at]vt.edu | (540) 231-9080 http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin -- gmx-users mailing listgmx-users@gromacs.org http://lists.gromacs.org/mailman/listinfo/gmx-users Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/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/Support/Mailing_Lists