[gmx-users] li lin wants to share a link | Gromacs

2012-02-24 Thread ros...@kth.se
li lin wants to a share a link on the Gromacs wiki: http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List/Search li lin says: Dear all, Recently I made some tests, and I want to get same potential with the following method. (gmx 4.0.5.) Firstly, I run a complete simulation as the

Re: [gmx-users] li lin wants to share a link | Gromacs

2012-02-24 Thread Mark Abraham
On 24/02/2012 7:05 PM, ros...@kth.se wrote: li lin wants to a share a link on the Gromacs wiki: http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List/Search li lin says: Dear all, Recently I made some tests, and I want to get same potential with the following method. (gmx 4.0.5.)

Re: [gmx-users] Options for the buckingham potential

2012-02-24 Thread ramesh cheerla
Thank you for your valuble suggestion. On Fri, Feb 24, 2012 at 12:02 PM, Mark Abraham mark.abra...@anu.edu.auwrote: On 24/02/2012 5:15 PM, Mark Abraham wrote: On 24/02/2012 4:06 PM, ramesh cheerla wrote: Dear Gromacs users, i am planning to use buckingham

Re: [gmx-users] Umbrella_pull_simulation

2012-02-24 Thread shahid nayeem
My system is a protein-protein complex. After pulling I selected windows at 0.1nm from an initial COM distance of 3.65nm to 5.7nm and thereafter windows at 0.2nm spacing was selected upto 7.5nm. In each window 10ns MD was done and g_wham command was run neglecting the fist 1ns run for

[gmx-users] g_analyze -ee

2012-02-24 Thread dina dusti
Dear Gromacs Specialists, I have one problem about g_analyze -ee, May I ask you to help me, Please? When I do this program as g_analyze -f .xvg -av -ee, it is as followed: Read 4 sets of 83001 points, dt = 6   std. dev.    relative deviation of    

Re: [gmx-users] Adding new residue to the the force field

2012-02-24 Thread ramesh cheerla
Dear Mark, Thanks for your reply. I am not using both LJ and Buckingham in the same simulation, but trying to use only Buckingham . As my system is a polymer which is made of the monomer unit that's what i am adding to the force field. Our polymer system uses the parameters

[gmx-users] about protein in periodic boundary condition..

2012-02-24 Thread rama david
Hi GROMACS user, With the help of command genbox -ci molecule name -nmol molecule no required -o out put file -box dimension I get the desired no of added protein to cell. But I found that some part of protein is outside the cell .. Is such system is good for

Re: [gmx-users] Umbrella_pull_simulation

2012-02-24 Thread Justin A. Lemkul
shahid nayeem wrote: My system is a protein-protein complex. After pulling I selected windows at 0.1nm from an initial COM distance of 3.65nm to 5.7nm and thereafter windows at 0.2nm spacing was selected upto 7.5nm. In each window 10ns MD was done and g_wham command was run neglecting the

Re: [gmx-users] about protein in periodic boundary condition..

2012-02-24 Thread Justin A. Lemkul
rama david wrote: Hi GROMACS user, With the help of command genbox -ci molecule name -nmol molecule no required -o out put file -box dimension I get the desired no of added protein to cell. But I found that some part of protein is outside the cell .. Is such

Re: [gmx-users] Umbrella_pull_simulation

2012-02-24 Thread shahid nayeem
I thought this time to be sufficient without any reasonable basis. shahid Nayeem On Fri, Feb 24, 2012 at 7:40 PM, Justin A. Lemkul jalem...@vt.edu wrote: shahid nayeem wrote: My system is a protein-protein complex. After pulling I selected windows at 0.1nm from an initial COM distance of

Re: [gmx-users] Umbrella_pull_simulation

2012-02-24 Thread Justin A. Lemkul
shahid nayeem wrote: I thought this time to be sufficient without any reasonable basis. If you pick an arbitrary time frame, you're going to get results with arbitrary accuracy. Your time frame should be decided based on a whole host of factors, not the least of which include an

[gmx-users] Calculate Dihedrals

2012-02-24 Thread francesca vitalini
Hi GROMACS community, I have one problem: I need to calculate the dihedrals for a big system (so not possible to do it by hand using vmd) starting from a gro file. Do you know if it exists an efficient way of doing it? I've seen that the g_chi option exists, however, it doesn't suit for me as it

Re: [gmx-users] Calculate Dihedrals

2012-02-24 Thread Justin A. Lemkul
francesca vitalini wrote: Hi GROMACS community, I have one problem: I need to calculate the dihedrals for a big system (so not possible to do it by hand using vmd) starting from a gro file. Do you know if it exists an efficient way of doing it? I've seen that the g_chi option exists, however,

Re: [gmx-users] Umbrella_pull_simulation

2012-02-24 Thread shahid nayeem
My protein complex interface has a hydrophobic core and on each side of this core at the edge are two hydrogen bonds. The Hydrogen bond on one side is between Arg and Asp and another side it is between Arg and Glu. Its experimental Kd is in nanomolar regions. How should I decide the length of the

[gmx-users] gen_temp, gen_seed

2012-02-24 Thread Steven Neumann
Dear Gmx Users, Dear Justin I am wondering which approach for generating velocities in protein or protein-ligand system is a proper one: a) NVT gen_vel = yes gen_temp=298 gen_seed=-1 NPT gen_vel=no MDgen_vel=no b) NVT gen_vel = yes gen_temp=298

Re: [gmx-users] Positive Coul. recip. term

2012-02-24 Thread Denny Frost
Thank you both for your replies. I currently have another ionic liquid running just fine on the same gromacs build (compiled the tpr file yesterday), so I am reluctant to conclude that the problem is with the linking. Please let me know if you disagree. The force field I am using was published

[gmx-users] Domain decomposition problems

2012-02-24 Thread Olivia Waring
Dear Gromacs users, Could someone please shed light on the following problem? My system consists of an alkanethiol SAM (in the xy plane) with a layer of water on top. The mdrun command works when the box vectors are 1.1, 1.3, and 2.5, but the SAM flies apart and the simulation crashes. I'm trying

[gmx-users] Re: gen_temp, gen_seed

2012-02-24 Thread Justin A. Lemkul
Steven Neumann wrote: Dear Gmx Users, Dear Justin I am wondering which approach for generating velocities in protein or protein-ligand system is a proper one: a) NVT gen_vel = yes gen_temp=298 gen_seed=-1 NPT gen_vel=no MDgen_vel=no b) NVT gen_vel = yes

Re: [gmx-users] Umbrella_pull_simulation

2012-02-24 Thread Justin A. Lemkul
shahid nayeem wrote: My protein complex interface has a hydrophobic core and on each side of this core at the edge are two hydrogen bonds. The Hydrogen bond on one side is between Arg and Asp and another side it is between Arg and Glu. Its experimental Kd is in nanomolar regions. How should

Re: [gmx-users] Re: Internal water in the membrane receptor

2012-02-24 Thread James Starlight
Mark, So as I understood the ussage of X-ray structure as the starting model where the internal water is already present might be good in case to avoid those sterric issues doesn't it ? What are additional options should I use for preparation of such system with pdb2gmx ? Should I use posres on

Re: [gmx-users] Umbrella sampling and WHAM along a curved pathway?

2012-02-24 Thread HAO JIANG
Hi, Jochen, Thanks a lot for your reply. I looked at the example header provided by g_wham -h. I just wonder if you can answer some more questions for me. The example header is # UMBRELLA 3.0 # Component selection: 0 0 1 # nSkip 1 # Ref. Group 'TestAtom' # Nr. of pull groups 2 # Group 1

[gmx-users] The solvation free energy of small ligand in water

2012-02-24 Thread Tanping Li
Dear gmx users, I know this is a old topic. I searched the mailing list however haven't find a answer. I am calculating the solvation free energy of a small ligand in water, using FEP theory. Two steps are applied: firstly decrease the atomic charge to zero; secondly decrease the VDW to

Re: [gmx-users] rc or rlist in minimum image convention?

2012-02-24 Thread Juliette N.
On 23 February 2012 21:18, Mark Abraham mark.abra...@anu.edu.au wrote: On 24/02/12, *Juliette N. *joojoojo...@gmail.com wrote: On 23 February 2012 20:07, Mark Abraham mark.abra...@anu.edu.au wrote: On 24/02/2012 10:55 AM, Juliette N. wrote: Hi all, My average size is 2.9 nm obtained

[gmx-users] HBonds with VMD

2012-02-24 Thread Hovakim Grabski
Dear GROMACS users, I'm a novice user, I've been trying to find a way to visualize hydrogen bonds with VMD,but  I haven't been successful. After running a simulation  of 13 Hypericin molecules solved in Water, I loaded the.gro file and then the trajectory file .xtc, after that in the Create

Re: [gmx-users] HBonds with VMD

2012-02-24 Thread Justin A. Lemkul
Hovakim Grabski wrote: Dear GROMACS users, I'm a novice user, I've been trying to find a way to visualize hydrogen bonds with VMD,but I haven't been successful. After running a simulation of 13 Hypericin molecules solved in Water, I loaded the.gro file and then the trajectory file .xtc,

Re: [gmx-users] rc or rlist in minimum image convention?

2012-02-24 Thread Justin A. Lemkul
Juliette N. wrote: On 23 February 2012 21:18, Mark Abraham mark.abra...@anu.edu.au mailto:mark.abra...@anu.edu.au wrote: On 24/02/12, *Juliette N. *joojoojo...@gmail.com mailto:joojoojo...@gmail.com wrote: On 23 February 2012 20:07, Mark Abraham mark.abra...@anu.edu.au

RE: [gmx-users] HBonds with VMD

2012-02-24 Thread Marzinek, Jan
Try the representation: all not water Then set it to the HBONDS And refer to the VMD mailing list next time! :) Jan From: gmx-users-boun...@gromacs.org [gmx-users-boun...@gromacs.org] on behalf of Justin A. Lemkul [jalem...@vt.edu] Sent: Friday, February

RE: [gmx-users] HBonds with VMD

2012-02-24 Thread Marzinek, Jan
Dear GROMACS users, I'm a novice user, I've been trying to find a way to visualize hydrogen bonds with VMD,but I haven't been successful. After running a simulation of 13 Hypericin molecules solved in Water, I loaded the.gro file and then the trajectory file .xtc, after that in the Create

Re: [gmx-users] Domain decomposition problems

2012-02-24 Thread Mark Abraham
On 25/02/2012 4:52 AM, Olivia Waring wrote: Dear Gromacs users, Could someone please shed light on the following problem? My system consists of an alkanethiol SAM (in the xy plane) with a layer of water on top. The mdrun command works when the box vectors are 1.1, 1.3, and 2.5, but the SAM

Re: [gmx-users] Adding new residue to the the force field

2012-02-24 Thread Mark Abraham
On 24/02/2012 11:56 PM, ramesh cheerla wrote: Dear Mark, Thanks for your reply. I am not using both LJ and Buckingham in the same simulation, but trying to use only Buckingham . As my system is a polymer which is made of the monomer unit that's what i am adding to the

Re: [gmx-users] Re: Internal water in the membrane receptor

2012-02-24 Thread Mark Abraham
On 25/02/2012 5:37 AM, James Starlight wrote: Mark, So as I understood the ussage of X-ray structure as the starting model where the internal water is already present might be good in case to avoid those sterric issues doesn't it ? Yes. This was one of the options I suggested earlier.

[gmx-users] Acpype Parameter Error?

2012-02-24 Thread Kyle Greenway
Hello, This email is directed mainly to Alan, who created Acpype. I've noticed that Acpype has assigned dihedral constants as 0.65084 for many dihedrals of the form X -c3-n4-X, X -c3-c3-X, and others, in my generated GROMACS .itp files. These dihedrals have values of 1.400 in the amber

Re: [gmx-users] rc or rlist in minimum image convention?

2012-02-24 Thread Juliette N.
On 24 February 2012 16:52, Justin A. Lemkul jalem...@vt.edu wrote: Juliette N. wrote: On 23 February 2012 21:18, Mark Abraham mark.abra...@anu.edu.au mailto: mark.abra...@anu.edu.**au mark.abra...@anu.edu.au wrote: On 24/02/12, *Juliette N. *joojoojo...@gmail.com

[gmx-users] extending a simulation after it get stopped because of queue time limit

2012-02-24 Thread priya thiyagarajan
hello sir, i was performing simulation for 30ns. due to queue time limit my mdrun stopped at 11.6ns.. then i extended my simulation using these two commands *tpbconv -s md.tpr -extend 2 -o newmd.tpr mdrun -s newmd.tpr -o md.trr -c md.gro -e md.edr -g md.log -cpi state.cpt -x traj.xtc

Re: [gmx-users] extending a simulation after it get stopped because of queue time limit

2012-02-24 Thread Chandan Choudhury
On Sat, Feb 25, 2012 at 12:09 PM, priya thiyagarajan priya.thiyagaraja...@gmail.com wrote: hello sir, i was performing simulation for 30ns. due to queue time limit my mdrun stopped at 11.6ns.. then i extended my simulation using these two commands *tpbconv -s md.tpr -extend 2 -o

[gmx-users] gromacs: pdb2gmx error

2012-02-24 Thread venkatesh s
Hi i got that error in the step of pdb2gmx what can i do i read also http://www.gromacs.org/Documentation/Errors#Residue_%27XXX%27_not_found_in_residue_topology_database can any one tell exact reason and how i handle further thanking you - All occupancies are one Opening