Re: [gmx-users] Ewald summation: n=0

2009-04-23 Thread Mark Abraham
Steve Fiedler wrote: Gromacs Users: While porting field parameters to Gromacs, I ran across unexpected electrostatic energies returned by the ewald summation method. To simplify the problem, I reduced the system to two partially charged particles. For the simplest case of only one box, n =

Re: [gmx-users] Identical energies generated in a rerun calculation ... but ...

2009-04-23 Thread Mark Abraham
Mark Abraham wrote: Luca Mollica wrote: Dear users, I am experiencing some troubles in using GMX 4.0.4 trying to rerun trajectories and obtaining energies for single residues and/or portion of proteins. According to the mdrun command help and to my experience with the previous versions of GR

[gmx-users] Re: Re:g_cluster output

2009-04-23 Thread sarbani chattopadhyay
Hi, Actually I am not confused regarding the clustering method that is being used in default. I know that the method that is used is "Single linkage " method. I got to know the options that I had used in my "g_cluster" command after going through g_cluster -h. Regarding the command "g_cluster

[gmx-users] topology for oxygen molecule

2009-04-23 Thread Sunil Thapa
Respectable Justin/Mark For the study of diffusion of oxygen in water, I tried to search the itp file for oxygen molecule but could not find. Therefore, I have tried to make it myself.   Please tell me whether it is fine.   [ moleculetype ] ; molname nrexcl OMOL 1 ; as there are two atoms in an o

[gmx-users] about salt bridge

2009-04-23 Thread mary
To whom it may concern, when I used g_saltbr to compute some salt bridges in my system, I couldn't understand the option "-t". The tutorial tell me that it's trunc distance and the default is 1000. Then I want to find the definition of the salt bridge in Gromacs. However, I couldn't fi

[gmx-users] Re: gmx-users Digest, Vol 60, Issue 116

2009-04-23 Thread nitu sharma
ilayer > > [ molecules ] > ; Compound#mols > Protein_A 1 > Protein_B 1 > DMPC 128 > > but when I gave this topology file to grompp it shws error- > no. of coordinates in coordinate file doesn,t mach the topology file . > > Can yo

[gmx-users] g_sas

2009-04-23 Thread Cheong Wee Loong, Daniel
Dear all, I am interested to calculate the hydrophobic and hydrophilic area of the surface of the protein layer I am simulating. It looked like g_sas would be able to give me what I was looking for. But I was wondering what the difference is between the calculation group and the output group.

[gmx-users] Inclusion of position_restraints

2009-04-23 Thread Tsjerk Wassenaar
Hi Nitu, >> When replying, please edit your Subject line so it is more specific >> than "Re: Contents of gmx-users digest..." Please!... And don't include the whole digest. > > Dear Mark > Thanks for your reply. as u ask- how many atoms is--Protein A- 5244 > > protein B- 4658 >        

Re: [gmx-users] g_sas

2009-04-23 Thread Justin A. Lemkul
Cheong Wee Loong, Daniel wrote: Dear all, I am interested to calculate the hydrophobic and hydrophilic area of the surface of the protein layer I am simulating. It looked like g_sas would be able to give me what I was looking for. But I was wondering what the difference is between the

[gmx-users] problem in making topology file for protein inlipid bilayer

2009-04-23 Thread nitu sharma
Dear Mark Thanks for your consideration. I am replying you again I think it will be more easy for u in understanding my problem. as u ask- how many atoms is--Protein A- 5244 Protein B- 4658 And in one dmpc molecule there is 46 atoms so

Re: [gmx-users] topology for oxygen molecule

2009-04-23 Thread Justin A. Lemkul
Sunil Thapa wrote: Respectable Justin/Mark For the study of diffusion of oxygen in water, I tried to search the itp file for oxygen molecule but could not find. Therefore, I have tried to make it myself. Please tell me whether it is fine. [ moleculetype ] ; molname nrexcl OMOL 1 ; as the

Re: [gmx-users] about salt bridge

2009-04-23 Thread Mark Abraham
mary wrote: > > > To whom it may concern, > > when I used g_saltbr to compute some salt bridges in my system, I > couldn't understand the option "-t". The tutorial tell me that it's > trunc distance and the default is 1000. Then I want to find the > definition of the salt bridge in Gr

Re: [gmx-users] problem in making topology file for protein inlipid bilayer

2009-04-23 Thread Mark Abraham
nitu sharma wrote: Dear Mark Thanks for your consideration. I am replying you again I think it will be more easy for u in understanding my problem. as u ask- how many atoms is--Protein A- 5244 Protein B- 4658 And in one dmpc molecul

Re: [gmx-users] Re: Re:g_cluster output

2009-04-23 Thread Mark Abraham
sarbani chattopadhyay wrote: Hi, Actually I am not confused regarding the clustering method that is being used in default. I know that the method that is used is "Single linkage " method. I got to know the options that I had used in my "g_cluster" command after going through g_cluster -h. R

[gmx-users] problem with genbox -ci option

2009-04-23 Thread Ken Rotondi
Hello All, I have managed to get urea into the G53a6 ffield with suitable adjustments to the appropriate .pdb, .rtp, and .hdb files. I am trying to use the resulting urea.gro file (I have also tried with a urea.pdb file generated from the urea.gro file using pdb2gmx) in genbox using the -

Re: [gmx-users] problem with genbox -ci option

2009-04-23 Thread Justin A. Lemkul
Which version of Gromacs? There was a bug in genbox that was fixed for version 4.0.4. -Justin Ken Rotondi wrote: Hello All, I have managed to get urea into the G53a6 ffield with suitable adjustments to the appropriate .pdb, .rtp, and .hdb files. I am trying to use the resulting urea.gro

Re: [gmx-users] problem with genbox -ci option

2009-04-23 Thread Ken Rotondi
Sorry, vers 4.0.4 here on an intel macbook running OSX.5.5 On Apr 23, 2009, at 11:06 AM, Justin A. Lemkul wrote: Which version of Gromacs? There was a bug in genbox that was fixed for version 4.0.4. -Justin Ken Rotondi wrote: Hello All, I have managed to get urea into the G53a6 ffield w

Re: [gmx-users] problem with genbox -ci option

2009-04-23 Thread Justin A. Lemkul
Ken Rotondi wrote: Sorry, vers 4.0.4 here on an intel macbook running OSX.5.5 Does it work if you use -ci to insert your molecule, then in a separate step, use -cs to solvate? -Justin On Apr 23, 2009, at 11:06 AM, Justin A. Lemkul wrote: Which version of Gromacs? There was a bug in

Re: [gmx-users] problem with genbox -ci option

2009-04-23 Thread Ken Rotondi
Dear Justin, All, Yes, running as two separate steps works, is this a known genbox bug? Ken -works- pdb2gmx -f diala.pdb -o diala.gro -p diala.top -ff G53a6 -ignh editconf -f diala.gro -o dialabox.gro -bt dodecahedron -box 3.1 genbox -cp dialabox.gro -ci urea.gro -nmol 1 -o dialaurea.gro -p

Re: [gmx-users] problem with genbox -ci option

2009-04-23 Thread Justin A. Lemkul
Ken Rotondi wrote: Dear Justin, All, Yes, running as two separate steps works, is this a known genbox bug? I don't know that it's a bug per se, more of a limitation. The program is trying to solvate the structure, as well as cram in X molecules of something else all at once and probably

Re: [gmx-users] problem with genbox -ci option

2009-04-23 Thread Ken Rotondi
Indeed this seems likely. I was able to add 125 urea and solvate w/out incident doing it in two steps. Thanks again, Ken On Apr 23, 2009, at 1:31 PM, Justin A. Lemkul wrote: Ken Rotondi wrote: Dear Justin, All, Yes, running as two separate steps works, is this a known genbox bug? I don

[gmx-users] Coarse-graining in gromacs

2009-04-23 Thread ??????? ?????
Dear gmx-user I am a very new user in gromacs. I encounter the problem with coarse-grained option in gromacs. I have installed the option from the MARTINI tutorials which is “gromacs_reverse.tar.gz”. It is ok with following the tutorial to do the coarse-graining but when I want to do all atomi

[gmx-users] problem with the calculation of a propene molecule

2009-04-23 Thread Saskia Frenzel
Dear Ladies and Gentlemen, I from the University of Applied Sciences in Mittweida currently deals with the calculation of a propene molecule using the program Gromacs. For testing proposer I tried to generate my own force field. However, when creating the top-file with this command "grompp_d -f pr

Re: [gmx-users] Coarse-graining in gromacs

2009-04-23 Thread Mark Abraham
??? ? wrote: Dear gmx-user I am a very new user in gromacs. I encounter the problem with coarse-grained option in gromacs. I have installed the option from the MARTINI tutorials which is “gromacs_reverse.tar.gz”. It is ok with following the tutorial to do the coarse-graining but

Re: [gmx-users] problem with the calculation of a propene molecule

2009-04-23 Thread Mark Abraham
Saskia Frenzel wrote: Dear Ladies and Gentlemen, I from the University of Applied Sciences in Mittweida currently deals with the calculation of a propene molecule using the program Gromacs. For testing proposer I tried to generate my own force field. However, when creating the top-file with this

RE: [gmx-users] g_sas

2009-04-23 Thread Cheong Wee Loong, Daniel
Thanks Justin for the reply. I did read the help and manual and understand that the output file can be the whole or part of the calculation group. What I don't quite understand is what the output group represents. I am assuming that the calculation group would be the group of atoms that will b

Re: [gmx-users] g_sas

2009-04-23 Thread Justin A. Lemkul
Cheong Wee Loong, Daniel wrote: Thanks Justin for the reply. I did read the help and manual and understand that the output file can be the whole or part of the calculation group. What I don't quite understand is what the output group represents. I am assuming that the calculation group would

RE: [gmx-users] g_sas

2009-04-23 Thread Cheong Wee Loong, Daniel
Thanks for the explanation. It is much clearer now. Although I still don't quite understand why we can't just use Protein A as the calculation group AND output group to find the SASA of protein A. I know it states that the calculation group should be all non-solvent atoms, but I guess I am ju

Re: [gmx-users] g_sas

2009-04-23 Thread Justin A. Lemkul
Cheong Wee Loong, Daniel wrote: Thanks for the explanation. It is much clearer now. Although I still don't quite understand why we can't just use Protein A as the calculation group AND output group to find the SASA of protein A. I know it states that the calculation group should be all non-s

[gmx-users] about salt bridge

2009-04-23 Thread xiongxuqiong
To whom it may concern, when I used g_saltbr to compute some salt bridges in my system, I couldn't understand the option "-t". The tutorial tell me that it's trunc distance and the default is 1000. Then I want to find the definition of the salt bridge in Gromacs. However, I couldn't f

Re: [gmx-users] about salt bridge

2009-04-23 Thread Justin A. Lemkul
xiongxuqiong wrote: To whom it may concern, when I used g_saltbr to compute some salt bridges in my system, I couldn't understand the option "-t". The tutorial tell me that it's trunc distance and the default is 1000. Then I want to find the definition of the salt bridge in Gromacs

Re: [gmx-users] about salt bridge (documentation suggestion)

2009-04-23 Thread Justin A. Lemkul
Justin A. Lemkul wrote: From g_saltbr -h: "A minimum distance can be given, (eg. the cut-off), then groups that are never closer than that distance will not be plotted." Basically, if a salt bridge occurs within the truncation distance specified by -t, it will be plotted. Using the defau

RE: [gmx-users] g_sas

2009-04-23 Thread Cheong Wee Loong, Daniel
Ah ok. Fair enough. Thanks! -Original Message- From: gmx-users-boun...@gromacs.org [mailto:gmx-users-boun...@gromacs.org] On Behalf Of Justin A. Lemkul Sent: Friday, April 24, 2009 10:48 AM To: Gromacs Users' List Subject: Re: [gmx-users] g_sas Cheong Wee Loong, Daniel wrote: > Than

Re: [gmx-users] g_sas

2009-04-23 Thread Mark Abraham
Cheong Wee Loong, Daniel wrote: Thanks for the explanation. It is much clearer now. Although I still don't quite understand why we can't just use Protein A as the calculation group AND output group to find the SASA of protein A. If A is complexed to B in solvent you might want the true SASA o

[gmx-users] grompp error

2009-04-23 Thread nitu sharma
Dear MARK Now the problem of matching of topology file and coordinate file solved but the other two errors come during running grompp in the second step of inflategro . my command line is like this- grompp -f em.mdp -p topolinflate.top -c inflated_system.gro -o inflated-em.tp

RE: [gmx-users] grompp error

2009-04-23 Thread Dallas B. Warren
Both of those errors are telling you want you need to do. In your topol_A.itp file, there are no parameters for a particular proper dihedral. So, define one. Same with the next message. And it even tells you which line the issue is on in the .itp file, so very easy to track down which bonds are

Re: [gmx-users] grompp error

2009-04-23 Thread Mark Abraham
nitu sharma wrote: Dear MARK Now the problem of matching of topology file and coordinate file solved but the other two errors come during running grompp in the second step of inflategro . my command line is like this- grompp -f em.mdp -p topolinflate.top -c inflated_syst

Re: [gmx-users] g_sas

2009-04-23 Thread Tsjerk Wassenaar
Hi, Of course you can. But if part of Protein A is buried in an interface, doing the SAS calculation over A only will also include the surface interface; it will give the total surface of A. That may actually be handy if you want to determine what the buried surface is. You calculate the total sur

RE: [gmx-users] g_sas

2009-04-23 Thread Cheong Wee Loong, Daniel
Hi all. Thanks for all your replies. I think I understand it now and can definitely see how having the calculation and output group would add to the flexibility and utilty of g_sas. Thanks. Daniel -Original Message- From: gmx-users-boun...@gromacs.org [mailto:gmx-users-boun...@grom