[gmx-users] g_hbond output analysis

2008-08-26 Thread parthiban
Message: 1 Date: Tue, 26 Aug 2008 15:38:56 +0530 (IST) From: [EMAIL PROTECTED] Subject: [gmx-users] g_hbond output analysis To: gmx-users@gromacs.org Message-ID: [EMAIL PROTECTED] Content-Type: text/plain;charset=iso-8859-1 hi i used the -num option for g_hbond run. since i new

Re: [gmx-users] g_hbond distance distribution problem

2008-05-18 Thread Dechang Li
-- Message: 4 Date: Sun, 18 May 2008 00:50:30 +0200 From: Xavier Periole [EMAIL PROTECTED] Subject: Re: [gmx-users] g_hbond distance distribution problem To: Discussion list for GROMACS users gmx-users@gromacs.org Message-ID: [EMAIL PROTECTED] Content-Type: text/plain

Re: [gmx-users] g_hbond distance distribution problem

2008-05-18 Thread David van der Spoel
Dechang Li wrote: -- Message: 4 Date: Sun, 18 May 2008 00:50:30 +0200 From: Xavier Periole [EMAIL PROTECTED] Subject: Re: [gmx-users] g_hbond distance distribution problem To: Discussion list for GROMACS users gmx-users@gromacs.org Message-ID: [EMAIL PROTECTED

[gmx-users] g_hbond distance distribution problem

2008-05-17 Thread Dechang Li
Dear all, I used command g_hbond to calculated the distance distribution of hbonds in my system. There are two columns in the file hbdist.xvg, as showed follow: @title Hydrogen Bond Distribution @xaxis label Hydrogen - Acceptor Distance (nm) @yaxis label @TYPE xy 0.0025

Re: [gmx-users] g_hbond distance distribution problem

2008-05-17 Thread Xavier Periole
Dechang Li wrote: Dear all, I used command g_hbond to calculated the distance distribution of hbonds in my system. There are two columns in the file hbdist.xvg, as showed follow: @title Hydrogen Bond Distribution @xaxis label Hydrogen - Acceptor Distance (nm) @yaxis label

[gmx-users] g_hbond -ins

2008-04-29 Thread cfgauna
Hi all: I´m using gromacs version 3.3.2 and a box type truncated octahedron. I need to analyze hydrogen bonds with solvent (water)insertion but -ins option of g_hbond is broken... What should I do? I´ve allready tried to analyze my trayectory file (maked with 3.3.2) with gromacs version 3.1.4 but

Re: [gmx-users] g_hbond -contact -xpm does not work

2008-01-07 Thread David van der Spoel
Michel Cuendet wrote: Hi, I'm trying to look at contacts between two proteins at the individual atom level. I have a bunch of analysis scripts that use the .xpm file and work for h-bonds, so I wanted to use the output of g_hbond 3.3.2 again with the option -contact. But: g_hbond -hbn -hbm

[gmx-users] g_hbond -contact -xpm does not work

2008-01-06 Thread Michel Cuendet
Hi, I'm trying to look at contacts between two proteins at the individual atom level. I have a bunch of analysis scripts that use the .xpm file and work for h-bonds, so I wanted to use the output of g_hbond 3.3.2 again with the option -contact. But: g_hbond -hbn -hbm -num -contact -r 0.4

Re: [gmx-users] g_hbond -contact -xpm does not work

2008-01-06 Thread David van der Spoel
Michel Cuendet wrote: Hi, I'm trying to look at contacts between two proteins at the individual atom level. I have a bunch of analysis scripts that use the .xpm file and work for h-bonds, so I wanted to use the output of g_hbond 3.3.2 again with the option -contact. But: g_hbond -hbn -hbm

Re: [gmx-users] g_hbond -contact -xpm does not work

2008-01-06 Thread David van der Spoel
Michel Cuendet wrote: Hi, I'm trying to look at contacts between two proteins at the individual atom level. I have a bunch of analysis scripts that use the .xpm file and work for h-bonds, so I wanted to use the output of g_hbond 3.3.2 again with the option -contact. But: g_hbond -hbn -hbm

Re: [gmx-users] g_hbond -contact -xpm does not work

2008-01-06 Thread David van der Spoel
Michel Cuendet wrote: Hi, I'm trying to look at contacts between two proteins at the individual atom level. I have a bunch of analysis scripts that use the .xpm file and work for h-bonds, so I wanted to use the output of g_hbond 3.3.2 again with the option -contact. But: g_hbond -hbn -hbm

Re: [gmx-users] g_hbond -contact -xpm does not work

2008-01-06 Thread David van der Spoel
Michel Cuendet wrote: Hi, I'm trying to look at contacts between two proteins at the individual atom level. I have a bunch of analysis scripts that use the .xpm file and work for h-bonds, so I wanted to use the output of g_hbond 3.3.2 again with the option -contact. But: g_hbond -hbn -hbm

[gmx-users] g_hbond (Version 3.3 vs Version 3.14)

2007-12-20 Thread Anil Kumar
Dear David, Thanks for your concern. i am using g_hbond (After vacuum md, just for one step) Command with options is g_hbond -f *.pdb -s *.tpr -num *_num.xvg -hx *_hbhelix.xvg and when i compare the result of v 3.3.1 and v 3.14, i found there is difference in *_hbhelix.xvg. Even for standard

[gmx-users] g_hbond (Version 3.3 vs Version 3.14)

2007-12-19 Thread Anil Kumar
Dear All, I found there is some problem with version 3.3 to g_hbond command. Earlier i was using version 3.14, and when i was using g_hbond(v3.14) i found correct distribution of hydrogen bond pattern (i.e., n-n+1,n-n+2, n-n+3 so on). But when i use g_hbond (v3.3) then i found it only gives

[gmx-users] g_hbond (Version 3.3 vs Version 3.14)

2007-12-19 Thread Anil Kumar
Dear All, I found there is some problem with version 3.3 to g_hbond command. Earlier i was using version 3.14, and when i was using g_hbond(v3.14) i found correct distribution of hydrogen bond pattern (i.e., n-n+1,n-n+2, n-n+3 so on). But when i use g_hbond (v3.3) then i found it only gives

Re: [gmx-users] g_hbond (Version 3.3 vs Version 3.14)

2007-12-19 Thread David van der Spoel
Anil Kumar wrote: Dear All, I found there is some problem with version 3.3 to g_hbond command. Earlier i was using version 3.14, and when i was using g_hbond(v3.14) i found correct distribution of hydrogen bond pattern (i.e., n-n+1,n-n+2, n-n+3 so on). But when i use g_hbond (v3.3) then i

[gmx-users] g_hbond (v 3.14) vs g_hbond (3.3)

2007-12-13 Thread Anil Kumar
Dear Users, I have done clustering of equilibrium ensemble of octa-alanine model system, and when i was trying to find the no. of hydrogen bonds for each cluster, i found different result with different versions. Is there some bug to version 3.3 in comparison of version 3.14 or am i doing the

[gmx-users] g_hbond command problem

2007-12-03 Thread wangcl
Dear all users, I want to analyze Hydrogen-bond among water molecules using g_hbond command. We use periodic boundary conditions in the simulation box. The question is: I don't know whether this g_hbond tool considers the periodic boundary conditions, that is, whether this tool

Re: [gmx-users] g_hbond command problem

2007-12-03 Thread David van der Spoel
wangcl wrote: Dear all users, I want to analyze Hydrogen-bond among water molecules using g_hbond command. We use periodic boundary conditions in the simulation box. The question is: I don't know whether this g_hbond tool considers the periodic boundary conditions, that is, whether

[gmx-users] g_hbond

2007-10-17 Thread Dechang Li
Dear all, I have two moleculars in the .trr. I want to use g_hbond to analyze the hydrogen bonds between them. Can g_hbond analyze the hydrogen bonds between two moleculars associated by a H2O molecular? The hydrogen bond may be the following form: -N-HO-HO=C-

Re: [gmx-users] g_hbond

2007-10-17 Thread Mark Abraham
Dechang Li wrote: Dear all, I have two moleculars in the .trr. I want to use g_hbond to analyze the hydrogen bonds between them. Can g_hbond analyze the hydrogen bonds between two moleculars associated by a H2O molecular? The hydrogen bond may be the following form:

[gmx-users] G_hbond and periodicity

2007-09-28 Thread tangxuan
Hi all, Does periodic boundary affect the calculation of hydrogen bond if my protein is not in one box? What about other calculations? Thank you, Tang jiaowei ___ gmx-users mailing listgmx-users@gromacs.org

Re: [gmx-users] G_hbond and periodicity

2007-09-28 Thread Tsjerk Wassenaar
Hi Tang Jiaowei, g_hbond takes periodicity into account, like most tools do. Do note that this doesn't work if you've performed rotational transformations (lsq-fitting). Cheers, Tsjerk On 9/28/07, tangxuan [EMAIL PROTECTED] wrote: Hi all, Does periodic boundary affect the calculation of

[gmx-users] g_hbond (method)

2007-06-06 Thread Anil Kumar
Dear all, I would like to know the actual cut-off angle used by gromacs package to calculate Hydrogen bond. since in User manual of version 3.2 and earlier, it was written á equal to or lesser than áHB = 60 degree While in User manual of version 3.3, it is mentioned as á equal to or lesser

Re: [gmx-users] g_hbond (method)

2007-06-06 Thread David van der Spoel
Anil Kumar wrote: Dear all, I would like to know the actual cut-off angle used by gromacs package to calculate Hydrogen bond. since in User manual of version 3.2 and earlier, it was written á equal to or lesser than áHB = 60 degree While in User manual of version 3.3, it is mentioned as á

[gmx-users] g_hbond in/output

2006-12-20 Thread OZGE ENGIN
Hi, I have some questions about g_hbond option. 1) My system consists of a protein molecule (a capped 1 alanine residue) (CH3-CO-ALA-NH-CH3) and 250 water molecules. After selecting the g_hbond option, it says that protein molecule consists of 12 elements; however, it consists much more

Re: [gmx-users] g_hbond in/output

2006-12-20 Thread Erik Marklund
20 dec 2006 kl. 17.34 skrev OZGE ENGIN: Hi, I have some questions about g_hbond option. 1) My system consists of a protein molecule (a capped 1 alanine residue) (CH3-CO-ALA-NH-CH3) and 250 water molecules. After selecting the g_hbond option, it says that protein molecule consists of 12

Re: Re: [gmx-users] g_hbond

2006-12-19 Thread Erik Marklund
the analysis into parts: solvent-solute ... /Erik Thanks in advance Ozge -Original Message- From: Erik Marklund [EMAIL PROTECTED] To: Discussion list for GROMACS users gmx-users@gromacs.org Date: Mon, 18 Dec 2006 13:48:18 +0100 Subject: Re: [gmx-users] g_hbond On Mon, 2006-12

[gmx-users] g_hbond

2006-12-18 Thread OZGE ENGIN
Hi, I have a problem with g_hbond. I want to investigate the H bonds formed not only within the protein but also between water molecules and the protein.In this respect, I select the group 0 -- 'system' option for both selections, an error message like the following one appears: Segmentation

Re: [gmx-users] g_hbond

2006-12-18 Thread Erik Marklund
On Mon, 2006-12-18 at 13:34 +0200, OZGE ENGIN wrote: Hi, I have a problem with g_hbond. I want to investigate the H bonds formed not only within the protein but also between water molecules and the protein.In this respect, I select the group 0 -- 'system' option for both selections, an

Re: [gmx-users] g_hbond

2006-12-18 Thread Diego Enry
Maybe not enough memory. Try using -shell . On 12/18/06, Erik Marklund [EMAIL PROTECTED] wrote: On Mon, 2006-12-18 at 13:34 +0200, OZGE ENGIN wrote: Hi, I have a problem with g_hbond. I want to investigate the H bonds formed not only within the protein but also between water molecules and

Re: [gmx-users] g_hbond

2006-12-18 Thread Florian Haberl
hi, On Monday 18 December 2006 12:34, OZGE ENGIN wrote: Hi, I have a problem with g_hbond. I want to investigate the H bonds formed not only within the protein but also between water molecules and the protein.In this respect, I select the group 0 -- 'system' option for both selections, an

Re: Re: [gmx-users] g_hbond

2006-12-18 Thread OZGE ENGIN
Message- From: Erik Marklund [EMAIL PROTECTED] To: Discussion list for GROMACS users gmx-users@gromacs.org Date: Mon, 18 Dec 2006 13:48:18 +0100 Subject: Re: [gmx-users] g_hbond On Mon, 2006-12-18 at 13:34 +0200, OZGE ENGIN wrote: Hi, I have a problem with g_hbond. I want to investigate

Re: [gmx-users] g_hbond -contact option

2006-11-23 Thread Yang Ye
Besides used for h-bond, g_hbond can be used for determine whether the distance has fallen into a range (with -nocontact and -r and -r2). g_dist is not good for calculate for distance between two atoms but g_mindist can perform this job well. Actually, because of a compulsive input of tpr

Re: [gmx-users] g_hbond -contact option

2006-11-22 Thread Erik Marklund
carotti [EMAIL PROTECTED] To: gmx-users@gromacs.org Sent: Wednesday, November 22, 2006 10:46 AM Subject: [gmx-users] g_hbond -contact option Hi all, I have found this interesting option in g_hbond, and I would like to know something more about it. How does it work? Which is the difference

[gmx-users] g_hbond error Your computational box has shrunk too much

2006-10-18 Thread andrea carotti
Hi all, I'm beginning to analyze my trajectories, everything works fine with (g_rmsf, g_saltbr ...) but when i try to use g_hbond I receive this error: -- Program g_hbond, VERSION 3.3.1 Source code file: gmx_hbond.c, line: 631 Fatal error: Your computational box has shrunk too much.

[gmx-users] g_hbond

2006-10-11 Thread merc mertens
hello, after using g_hbond to generate an *.xpm file one sees in the *.xpm file the following: y-label: Hydrogen Bond Index and the y-axis ranges between 0 and 163. if i look in the additionally generated hbbond.ndx file, i see at the bottom of it the group [ hbonds_* ]. i assumed that the

Re: [gmx-users] g_hbond

2006-10-11 Thread Erik Marklund
[EMAIL PROTECTED] To: Discussion list for GROMACS users gmx-users@gromacs.org Sent: Wednesday, October 11, 2006 8:21 AM Subject: [gmx-users] g_hbond hello, after using g_hbond to generate an *.xpm file one sees in the *.xpm file the following: y-label: Hydrogen Bond Index and the y-axis ranges

Re: [gmx-users] g_hbond

2006-08-23 Thread Erik Marklund
Ok. I'm on it. A few things I can see right away, other things need a closer look. See comments below On Tue, 2006-08-22 at 16:02 -0500, Moore, Jonathan (J) wrote: I'm still having problems with g_hbond. I've submitted a bugzilla report. - The hbond existence matrix in the -hbm output does

RE: [gmx-users] g_hbond

2006-08-23 Thread Moore, Jonathan (J)
23, 2006 7:28 AM To: Discussion list for GROMACS users Subject: Re: [gmx-users] g_hbond Ok. I'm on it. A few things I can see right away, other things need a closer look. See comments below On Tue, 2006-08-22 at 16:02 -0500, Moore, Jonathan (J) wrote: I'm still having problems with g_hbond

[gmx-users] g_hbond

2006-08-22 Thread Moore, Jonathan (J)
I'm still having problems with g_hbond. I've submitted a bugzilla report. - The hbond existence matrix in the -hbm output does not match the list of hbonds in the ndx file from -hbn. - Also, it outputs No option -g when no -g command line option has been tried - Also, sometimes the number of

Re: [gmx-users] g_hbond

2006-08-10 Thread Erik Marklund
It was a hard nut to crack, but I think I've done it after some sifting the source code: Since the -merge option is switched on by default, the hydrogen indexes are more or less irrelevant for some parts of the code. What happens here is that g_hbond finds the hbond [10 12 22] (distance: 0.2971

Re: [gmx-users] g_hbond

2006-08-10 Thread Chinmay Das
Dear David, ok. Now I understand HDA criterion. Can you give me a reference where I can find the differences between HDA criteria and OHD criteria and why one of them is preferable than the other? vmd certainly uses the later criteria. Changing the angle criteria there 44 degrees in vmd picks

[gmx-users] g_hbond

2006-08-09 Thread Chinmay Das
Hi, I am having trouble understanding the criteria for hydrogen bond and how it is implemented in g_hbond. Defining hydrogen bond for O1--H..O2 by requiring distance(O1, O2) 0.35 nm and angle(O1-H, H-O2) 30 degrees, my own code consistently give lower number of hydrogen bonds than g_hbond. But

Re: [gmx-users] g_hbond

2006-08-09 Thread David van der Spoel
Chinmay Das wrote: Hi, I am having trouble understanding the criteria for hydrogen bond and how it is implemented in g_hbond. Defining hydrogen bond for O1--H..O2 by requiring distance(O1, O2) 0.35 nm and angle(O1-H, H-O2) 30 angle is defined as H-D-A in g_hbond. is that what you mean?

Re: [gmx-users] g_hbond

2006-08-09 Thread Chinmay Das
Dear David, Partly yes, I am not following H-D-A angle trying to visualize two water molecules I hope the mail program will not try to format the ascii cartoon below: H1--O1 \ /H3 H2 O2 \ H4 Which is the angle I should consider? Luzer-Chandler

Re: [gmx-users] g_hbond account dieletric medium ?

2006-08-03 Thread X.Periole
I have shown hydrogen bond formation between guanidino side chain of Arginine and carboxylate in sidechain of Aspartate out my simulation result using g_hbond. Now I see those two residues are well exposed in solvent. Since g_hbond compute hydrogen bonds based on only distance and angles,

Re: [gmx-users] g_hbond account dieletric medium ?

2006-08-03 Thread raja
Hi XAvier, Thanks again, I meant the resistance caused due to dielectric of water separating two charged residues. I doubt whether it will allow the transfer of charge between them in order for hydrogen bond to be effected. With regards, B.Nataraj On Thu, 03 Aug 2006 09:25:07 +0200, X.Periole

Re: [gmx-users] g_hbond account dieletric medium ?

2006-08-03 Thread David van der Spoel
raja wrote: Hi XAvier, Thanks again, I meant the resistance caused due to dielectric of water separating two charged residues. I doubt whether it will allow the transfer of charge between them in order for hydrogen bond to be effected. It's still unclear. The hydrogen bonds you find are there

Re: [gmx-users] g_hbond account dieletric medium ?

2006-08-03 Thread raja
Thanks spoel, The g_hbond merely report formation of hydrogen bonds because of two electronegative atoms come closure to certain distance and angle. The force field used to simulate is also not capable of modeling hydrogen bond (Afraid whether I am right). If the h_bond computes the hydrogen bonds

Re: [gmx-users] g_hbond account dieletric medium ?

2006-08-03 Thread David van der Spoel
raja wrote: Thanks spoel, The g_hbond merely report formation of hydrogen bonds because of two electronegative atoms come closure to certain distance and angle. The force field used to simulate is also not capable of modeling hydrogen bond (Afraid whether I am right). If the h_bond computes the

Re: [gmx-users] g_hbond account dieletric medium ?

2006-08-03 Thread raja
Dear Spoel, Thank for your elaborate answer and references. Some times this forum serves more than just for gromacs, yes it helps to understand some basics understanding of physical chemistry. I learned lot of fundamental concepts of protein modeling in this forum, which I could never understand

[gmx-users] g_hbond account dieletric medium ?

2006-08-02 Thread raja
Hi all, I have shown hydrogen bond formation between guanidino side chain of Arginine and carboxylate in sidechain of Aspartate out my simulation result using g_hbond. Now I see those two residues are well exposed in solvent. Since g_hbond compute hydrogen bonds based on only distance and angles,

[gmx-users] g_hbond -g

2006-06-29 Thread Diego Enry
Hi list, I couldn't find in the archives and 3.2 to 3.3 changes what happened to -g flag in g_hbond ? Is the old hbond.log information written somewhere else ? Diego Enry. // command line g_hbond_331i -g(intel compiled) // error message: Program g_hbond_331i, VERSION 3.3.1 Source code

RE: [gmx-users] g_hbond range checking error

2006-05-08 Thread Dallas B. Warren
Upgraded it to 3.3.1 now, still get the same error ... Using this command: g_hbond -f ../76.55-78.95ns.xtc -s ../md_000.tpr -n index.ndx -num hbnum_MGL-MGL I get this output: Reading file ../md_000.tpr,

Re: [gmx-users] g_hbond lifetime output

2006-04-24 Thread David van der Spoel
Michael Torrice wrote: I am confused about the output from the hblife.xvg file of g_hbond. What exactly is the p(t) output? Is it the percent distribution of a given lifetime? If so, I have received p(t) data that sums up to anything from 0.1 toi 2. Check my recent paper: David van der

RE: [gmx-users] g_hbond question

2006-03-22 Thread Moore, Jonathan (J)
-users] g_hbond question hi jonathan, the actually occuring hydrogen bonds are listed in the log-file. daniela On Tue, 2006-03-21 at 14:56 -0800, Moore, Jonathan (J) wrote: I created a hydrogen bond map using g_hbond of version 3.3. The manual indicates that for the hbmap the Ordering

[gmx-users] g_hbond question

2006-03-21 Thread Moore, Jonathan (J)
I created a hydrogen bond map using g_hbond of version 3.3. The manual indicates that for the hbmap the Ordering is identical to that in -hbn index file. However, in my case, the hbmap file contains 13 hydrogen bond indices, but the hbond.ndx file lists 20 sets of donors and acceptors. Why

Re: [gmx-users] g_hbond question

2006-03-21 Thread Daniela S. Mueller
hi jonathan, the actually occuring hydrogen bonds are listed in the log-file. daniela On Tue, 2006-03-21 at 14:56 -0800, Moore, Jonathan (J) wrote: I created a hydrogen bond map using g_hbond of version 3.3. The manual indicates that for the hbmap the Ordering is identical to that in

<    1   2   3