[gmx-users] Re: no output in g_dist

2011-09-12 Thread aiswarya pawar
Hi everyone, when am doing this g_dist between the protein and water am getting the distance between water-water. why it is so?? and is there a way i can find distance in each time frame. Thanks On Thu, Sep 8, 2011 at 5:55 PM, aiswarya pawar aiswarya.pa...@gmail.comwrote: Hi users, To get

[gmx-users] Re: problem during Energy minimization

2011-09-12 Thread Parul tew
gmx-users@gromacs.org Dear Gmx users, I am working on a transmembrane protein and my system contains protein, DPPC bilayer, water (spc) and ions. After preparing my system for simulation I have successfully performed the energy minimization, I am facing a problem at the nvt equilibration phase

[gmx-users] g_dist error

2011-09-12 Thread aiswarya pawar
Hi Users, Am using g_dist to find the distance between water and protein. but my output has the values of SOL-water distance. t: 1 136 SOL 2336 OW 0.772373 (nm) Thanks -- gmx-users mailing listgmx-users@gromacs.org http://lists.gromacs.org/mailman/listinfo/gmx-users Please search the

[gmx-users] Convert an unnormal residues contained protein from all-atoms model to coarse grained model.

2011-09-12 Thread Du Jiangfeng (BIOCH)
Hi gmx-users, I am trying to a study of GLA domain simulation in membrane by using coarse gained models. There is a problem coming from converting all-atom model to coarse gained because GLA domain contains some unusual residues (10 glutamate residues were modified by carboxylation) and also

[gmx-users] RE:g_lie

2011-09-12 Thread lloyd riggs
Message: 1 Date: Sun, 11 Sep 2011 13:01:28 +0200 From: gal.fra...@live.biu.ac.il Subject: [gmx-users] g_lie To: gmx-users@gromacs.org Message-ID: blu160-w601d2ab913084c46676afee9...@phx.gbl Content-Type: text/plain; charset=windows-1255 Greetings to all! I have a question for the g_lie

RE: [gmx-users] RE:g_lie

2011-09-12 Thread gal.fradin
Lloyd thank you very much for your answer! I'm going to check it out... Best regards, Gal Date: Mon, 12 Sep 2011 10:34:59 +0200 From: lloyd.ri...@gmx.ch To: gmx-users@gromacs.org Subject: [gmx-users] RE:g_lie Message: 1 Date: Sun, 11 Sep 2011 13:01:28 +0200 From:

[gmx-users] problem of lipid molecules entering voids of solvent during equilibration

2011-09-12 Thread Parul tew
Dear Gmx users, I am working on a transmembrane protein and my system contains protein, DPPC bilayer, water (spc) and ions. After preparing my system for simulation I have successfully performed the energy minimization, I am facing a problem at the nvt equilibration phase of 100ps, the lipid

[gmx-users] Representative structure from MD trajectories

2011-09-12 Thread bipin singh
Sir, thanks for your reply Yes I should have limited number of structure for a given point(pc1,pc2) but in my case, I have taken pc1 and pc2 as a range of points where the minima in FEL lies rather than a single point.So I think it would be relevant to apply the clustering as per your

Re: [gmx-users] problem of lipid molecules entering voids of solvent during equilibration

2011-09-12 Thread Justin A. Lemkul
Parul tew wrote: Dear Gmx users, I am working on a transmembrane protein and my system contains protein, DPPC bilayer, water (spc) and ions. After preparing my system for simulation I have successfully performed the energy minimization, I am facing a problem at the nvt equilibration phase of

Re: [gmx-users] Re: no output in g_dist

2011-09-12 Thread Justin A. Lemkul
aiswarya pawar wrote: Hi everyone, when am doing this g_dist between the protein and water am getting the distance between water-water. why it is so?? and is there a way i can find distance in each time frame. The distance is calculated based on the groups you tell g_dist to analyze.

Re: [gmx-users] g_dist error

2011-09-12 Thread Justin A. Lemkul
aiswarya pawar wrote: Hi Users, Am using g_dist to find the distance between water and protein. but my output has the values of SOL-water distance. t: 1 136 SOL 2336 OW 0.772373 (nm) This is not a water-water distance, it is the output of the -dist option telling you that water

[gmx-users] g_dist output error

2011-09-12 Thread aiswarya pawar
Hi gromacs Users, Am using g_dist to find the distance between water and protein. for that i made a index file such as- a CA(protein atoms) a OW(water atoms) i then run this cmd- g_dist -f md.xtc -s md.tpr -n index.ndx -e 500 -dist 0.35 instead of getting the result for protein-SOL distance .

Re: [gmx-users] g_dist output error

2011-09-12 Thread Justin A. Lemkul
aiswarya pawar wrote: Hi gromacs Users, Am using g_dist to find the distance between water and protein. for that i made a index file such as- a CA(protein atoms) a OW(water atoms) i then run this cmd- g_dist -f md.xtc -s md.tpr -n index.ndx -e 500 -dist 0.35 instead of getting the

Re: [gmx-users] g_dist error

2011-09-12 Thread aiswarya pawar
Hi Justin, Am using g_dist to find the distance between water and protein. for that i made a index file such as- a CA(protein atoms) a OW(water atoms) i then run this cmd- g_dist -f md.xtc -s md.tpr -n index.ndx -e 500 -dist 0.35 instead of getting the result for protein-SOL distance . my

Re: [gmx-users] g_dist error

2011-09-12 Thread aiswarya pawar
hi Justin, As far i referred the OW,HW1 etc are water atoms so how can it be distance between the SOL protein atoms, instead it is SOL water atoms. Thanks On Mon, Sep 12, 2011 at 4:49 PM, Justin A. Lemkul jalem...@vt.edu wrote: aiswarya pawar wrote: Hi Users, Am using g_dist to find the

Re: [gmx-users] g_dist output error

2011-09-12 Thread aiswarya pawar
hi Justin, As far i referred the OW,HW1 etc are water atoms so how can it be distance between the SOL protein atoms, instead it is SOL water atoms. Thanks On Mon, Sep 12, 2011 at 4:53 PM, Justin A. Lemkul jalem...@vt.edu wrote: aiswarya pawar wrote: Hi gromacs Users, Am using g_dist to

Re: [gmx-users] g_dist error

2011-09-12 Thread Justin A. Lemkul
aiswarya pawar wrote: Hi Justin, Am using g_dist to find the distance between water and protein. for that i made a index file such as- a CA(protein atoms) a OW(water atoms) i then run this cmd- g_dist -f md.xtc -s md.tpr -n index.ndx -e 500 -dist 0.35 instead of getting the result for

Re: [gmx-users] g_dist error

2011-09-12 Thread Justin A. Lemkul
aiswarya pawar wrote: hi Justin, As far i referred the OW,HW1 etc are water atoms so how can it be distance between the SOL protein atoms, instead it is SOL water atoms. The printed distance indicates that there is a certain water molecule that is just over 2 hydrogen bonding lengths

[gmx-users] How do I examine the results didn't write issue

2011-09-12 Thread lina
Hi, I am confused, why the job is running, but they did not write the .edr, .log .xtc ... files. Everything showed so normal, just not writing, when I used thread, all is fine, the writing is normal. just use mpi, I tried different nodes, still not work, now two hours later, still no writing,

[gmx-users] Nucleic Acid in vacuo

2011-09-12 Thread Matthias Ernst
Dear Gromacs users, I am currently working on porting forcefields from Gromacs to another program with the aim of simulating DNA. Therefore, I wrote a little script which reads in a PDB file, evaluates the Lennard-Jones and the Coulomb energy in vacuo with the forcefield data from Gromacs. I

Re: [gmx-users] How do I examine the results didn't write issue

2011-09-12 Thread Justin A. Lemkul
lina wrote: Hi, I am confused, why the job is running, but they did not write the .edr, .log .xtc ... files. Everything showed so normal, just not writing, when I used thread, all is fine, the writing is normal. just use mpi, I tried different nodes, still not work, now two hours later,

Re: [gmx-users] How do I examine the results didn't write issue

2011-09-12 Thread Mark Abraham
On 12/09/2011 10:23 PM, lina wrote: Hi, I am confused, why the job is running, but they did not write the .edr, .log .xtc ... files. Everything showed so normal, just not writing, when I used thread, all is fine, the writing is normal. just use mpi, I tried different nodes, still not work,

Re: [gmx-users] Nucleic Acid in vacuo

2011-09-12 Thread Justin A. Lemkul
Matthias Ernst wrote: Dear Gromacs users, I am currently working on porting forcefields from Gromacs to another program with the aim of simulating DNA. Therefore, I wrote a little script which reads in a PDB file, evaluates the Lennard-Jones and the Coulomb energy in vacuo with the forcefield

Re: [gmx-users] How do I examine the results didn't write issue

2011-09-12 Thread lina
On Mon, Sep 12, 2011 at 8:33 PM, Justin A. Lemkul jalem...@vt.edu wrote: lina wrote: Hi, I am confused, why the job is running, but they did not write the .edr, .log .xtc ... files. Everything showed so normal, just not writing, when I used thread, all is fine, the writing is normal.

Re: [gmx-users] How do I examine the results didn't write issue

2011-09-12 Thread lina
On Mon, Sep 12, 2011 at 8:35 PM, Mark Abraham mark.abra...@anu.edu.auwrote: On 12/09/2011 10:23 PM, lina wrote: Hi, I am confused, why the job is running, but they did not write the .edr, .log .xtc ... files. Everything showed so normal, just not writing, when I used thread, all is

[gmx-users] adding Ca2+ ions to the solvent (either randomly or restrained)

2011-09-12 Thread m r
  Dear all,   In an attempt to simulate the effect of adding calcium ion to the solvent around my protein of interest, I wonder how to add Ca2+ ions in the following cases (in GROMACS 4.5.4):   1) adding Ca2+ ions randomly distributed in the solvent;   2) adding Ca2+ such that it is initially

[gmx-users] g_rms matrix between wt and mutant

2011-09-12 Thread Shay Teaching
Hi all, (Gromacs 4.0.7): I am trying to make rms matrix between one Wt trajectory and one mutant trajectory using the following command: g_rms -f wt.xtc -f2 mutant.xtc -s wt.tpr -m -fit rot+trans -n wt_backbone.ndx The file wt_backbone.ndx contains the backbone of the protein (Backbone indices

Re: [gmx-users] g_rms matrix between wt and mutant

2011-09-12 Thread Justin A. Lemkul
Shay Teaching wrote: Hi all, (Gromacs 4.0.7): I am trying to make rms matrix between one Wt trajectory and one mutant trajectory using the following command: g_rms -f wt.xtc -f2 mutant.xtc -s wt.tpr -m -fit rot+trans -n wt_backbone.ndx The file wt_backbone.ndx contains the backbone of

Re: [gmx-users] g_rms matrix between wt and mutant

2011-09-12 Thread Shay Teaching
When I try to work the command on a small portion of the backbone it seems to work just fine. But when I try the entire backbone (which is composed of several *separate* chains) I am getting segmentation fault. Any workaround for that, so I can use the entire backbone? Thanks again, -Shay On Mon,

Re: [gmx-users] g_rms matrix between wt and mutant

2011-09-12 Thread Justin A. Lemkul
Shay Teaching wrote: When I try to work the command on a small portion of the backbone it seems to work just fine. But when I try the entire backbone (which is composed of several _separate_ chains) I am getting segmentation fault. Any workaround for that, so I can use the entire backbone?

[gmx-users] g_lie (again..)

2011-09-12 Thread גדעון לפידות
was scrubbed... URL: http://lists.gromacs.org/pipermail/gmx-users/attachments/20110912/03f71465/attachment-0001.html -- Message: 2 Date: Mon, 12 Sep 2011 21:12:26 +0800 From: lina lina.lastn...@gmail.com Subject: Re: [gmx-users] How do I examine the results didn't

[gmx-users] Problem when adding new residue in itp file

2011-09-12 Thread Du Jiangfeng (BIOCH)
Hi Guys, I want to run a protein which contains some carboxylated residues into a membrane. I have had to add a special residue into itp file since there is no any description in normal itp for GLA residue. I admit some values I added for GLA residue were ambiguous because I don't know the

Re: [gmx-users] Problem when adding new residue in itp file

2011-09-12 Thread Justin A. Lemkul
Du Jiangfeng (BIOCH) wrote: Hi Guys, I want to run a protein which contains some carboxylated residues into a membrane. I have had to add a special residue into itp file since there is no any description in normal itp for GLA residue. I admit some values I added for GLA residue were

[gmx-users] Running two versions of GMX at a time

2011-09-12 Thread chandran karunakaran
Hi Gmx users,     Are we able to run both the lower and higher version at a time? If possible, please help me how to do it. with thanks ***+ Dr.Karunakaran Chandran+ Biophysics Department + Medical College of Wisconsin + Milwaukee,

Re: [gmx-users] Running two versions of GMX at a time

2011-09-12 Thread Dommert Florian
On Mon, 2011-09-12 at 21:57 +0530, chandran karunakaran wrote: Hi Gmx users, Are we able to run both the lower and higher version at a time? If possible, please help me how to do it. with thanks You can always run different instances of a program at the same time. As you

[gmx-users] Re: amb2gmx.pl to convert GLYCAM topology

2011-09-12 Thread Yun Shi
Hi all, I am not a CS person, but I did find something in acpype.py as . if phase in [0, 180]: properDihedralsGmx45.append([item[0].atoms, phaseRaw, kPhi, period]) if not self.gmx45: if kPhi 0: V[period]

Re: [gmx-users] g_rms matrix between wt and mutant

2011-09-12 Thread Shay Teaching
Thanks, I'll try that, and post again if it works. On Mon, Sep 12, 2011 at 6:20 PM, Justin A. Lemkul jalem...@vt.edu wrote: Shay Teaching wrote: When I try to work the command on a small portion of the backbone it seems to work just fine. But when I try the entire backbone (which is

[gmx-users] Twin-range cut-off

2011-09-12 Thread Rabab Toubar
Does anyone has a clear definition or a reference for the Twin-range cut-off? I searched the internet and looked at the manual yet I couldn't get the idea   Thanks Rabab Toubar -- gmx-users mailing listgmx-users@gromacs.org http://lists.gromacs.org/mailman/listinfo/gmx-users Please search

Re: [gmx-users] Re: amb2gmx.pl to convert GLYCAM topology

2011-09-12 Thread Mark Abraham
On 13/09/2011 6:21 AM, Yun Shi wrote: Hi all, I am not a CS person, but I did find something in acpype.py as . if phase in [0, 180]: properDihedralsGmx45.append([item[0].atoms, phaseRaw, kPhi, period]) if not self.gmx45:

Re: [gmx-users] Twin-range cut-off

2011-09-12 Thread Mark Abraham
On 13/09/2011 8:56 AM, Rabab Toubar wrote: Does anyone has a clear definition or a reference for the Twin-range cut-off? I searched the internet and looked at the manual yet I couldn't get the idea Manual 4.6.3 and 3.4.7, and the top hits from Googling twin-range cutoff look pretty

Re: [gmx-users] Twin-range cut-off

2011-09-12 Thread Justin A. Lemkul
Rabab Toubar wrote: Does anyone has a clear definition or a reference for the Twin-range cut-off? I searched the internet and looked at the manual yet I couldn't get the idea A twin-range cutoff just means that your short-range cutoffs aren't all the same value, such that they form two

Re: [gmx-users] adding Ca2+ ions to the solvent (either randomly or restrained)

2011-09-12 Thread Mark Abraham
On 13/09/2011 12:03 AM, m r wrote: Dear all, In an attempt to simulate the effect of adding calcium ion to the solvent around my protein of interest, I wonder how to add Ca2+ ions in the following cases (in GROMACS 4.5.4): 1) adding Ca2+ ions randomly distributed in the solvent; genion does

[gmx-users] V(R) around protein

2011-09-12 Thread E. Nihal Korkmaz
Dear all, Is there any function to get the Volume around a peptide/protein as a function of R? The problem is g_rdf -surf doesn't have a volume correction factor, so my aim is to apply a volume correction factor manually by having V(R) around the protein. Thanks, -- Elif Nihal Korkmaz

Re: [gmx-users] V(R) around protein

2011-09-12 Thread Mark Abraham
On 13/09/2011 9:29 AM, E. Nihal Korkmaz wrote: Dear all, Is there any function to get the Volume around a peptide/protein as a function of R? The problem is g_rdf -surf doesn't have a volume correction factor, so my aim is to apply a volume correction factor manually by having V(R) around

Re: [gmx-users] V(R) around protein

2011-09-12 Thread E. Nihal Korkmaz
R is any distance from the surface of the protein. On Mon, Sep 12, 2011 at 6:31 PM, Mark Abraham mark.abra...@anu.edu.auwrote: On 13/09/2011 9:29 AM, E. Nihal Korkmaz wrote: Dear all, Is there any function to get the Volume around a peptide/protein as a function of R? The problem is g_rdf

Re: [gmx-users] V(R) around protein

2011-09-12 Thread E. Nihal Korkmaz
And I don't think g_sas would do what i want. It calculates the volume of the protein itself. I am interested in the V of space at a distance R from the surface of the protein, excluding the volume of the protein. In other words, i need the surface accessible volume as a function of distance from

Re: [gmx-users] V(R) around protein

2011-09-12 Thread Mark Abraham
On 13/09/2011 9:47 AM, E. Nihal Korkmaz wrote: And I don't think g_sas would do what i want. It calculates the volume of the protein itself. I am interested in the V of space at a distance R from the surface of the protein, excluding the volume of the protein. In other words, i need the

[gmx-users] How to use comm_mode and accelerate together?

2011-09-12 Thread Wang Xiaofeng
Dear all, To make molecules move, accelerate should be used. nstcomm is the frequency for center of mass motion removal But if they are used together, the molecule will not move. I made a model system: a box of pure water. Adding accelerate to SOL in Z direction. Firstly, comm_mode was set

[gmx-users] Time-varying Ion Charges

2011-09-12 Thread Zach Levine
Hi everyone, I'm trying to implement time-varying charges for various ions located in ions.itp in GROMACS 4.5.4, but I'm not sure where to begin. I've already had some luck implementing AC external electric fields by modifying /src/mdlib/sim_util.c, but this required minimal work from me and a

[gmx-users] Problem in VMD

2011-09-12 Thread meisam valizadeh kiamahalleh
Dear Gmx users I have run MD on carbon nanotube with solvated drug molecules inside and outside the tubes. Kindly would you please let me know how I can visualize and calculate the amount of molecules only inside the carbon nanotube during MD process. May I know what keywords I should use in graph

Re: [gmx-users] g_dist error

2011-09-12 Thread aiswarya . pawar
Even if I specify an atom say 1277 atom number to find distance against the OW atoms. I get the same result of SOL-OW distance. Is it a bug cause even after specifying one atom from a protein why doesn't it give me the result for the SOL. Thanks Sent from my BlackBerry® on Reliance Mobile,

Re: [gmx-users] g_dist error

2011-09-12 Thread Justin A. Lemkul
aiswarya.pa...@gmail.com wrote: Even if I specify an atom say 1277 atom number to find distance against the OW atoms. I get the same result of SOL-OW distance. Is it a bug cause even after specifying one atom from a protein why doesn't it give me the result for the SOL. As was suggested

Re: [gmx-users] g_dist error

2011-09-12 Thread aiswarya . pawar
Iam -dist option because I need the distance between two groups, excluding -dist gives me X,Y,Z output which I don't want. And am not specifying an -o. Sent from my BlackBerry® on Reliance Mobile, India's No. 1 Network. Go for it! -Original Message- From: Justin A. Lemkul

[gmx-users] Regarding clustering of MD data

2011-09-12 Thread bipin singh
Sir, thanks for your reply Yes I should have limited number of structure for a given point(pc1,pc2) but in my case, I have taken pc1 and pc2 as a range of points where the minima in FEL lies rather than a single point.So I think it would be relevant to apply the clustering as per your

Re: [gmx-users] dialanine using charmm27 ff

2011-09-12 Thread Mark Abraham
On 10/09/2011 10:23 AM, Sandeep Somani wrote: Hi Mark It worked ! pdb2gmx is now properly processing the pdb, and moreover the potential energy from gmx is within 1.27 kJ/mol of that from Charmm. More on energy comparison shortly. However, further tinkering of the rtp and pdb file was

Re: [gmx-users] g_dist error

2011-09-12 Thread Mark Abraham
On 13/09/2011 2:27 PM, aiswarya.pa...@gmail.com wrote: Iam -dist option because I need the distance between two groups That is not what g_dist -dist does. Please read g_dist -h. excluding -dist gives me X,Y,Z output which I don't want. And other output which you do, but you have to use

[gmx-users] Strangeness in gro file

2011-09-12 Thread Sweta Iyer
Hi all, I am trying to use the g_membed tool to insert my protein into a DMPC membrane. My protein is a dimer and hence I separated the two monomers by a TER card in the original pdb file before pdb2gmx step for gromacs to identify it as teo separate entities.. I then merged my protein file

Re: [gmx-users] g_dist error

2011-09-12 Thread aiswarya pawar
Hi Mark, The -dist option says- print all the atoms in group 2 that are closer than a certain distance to the center of mass of group 1. That means it should give me the distance from OW to protein atom. And when am already specifying only one atom from protein ie say 1322. why do i get this

Re: [gmx-users] g_dist error

2011-09-12 Thread Mark Abraham
On 13/09/2011 3:40 PM, aiswarya pawar wrote: Hi Mark, The -dist option says- print all the atoms in group 2 that are closer than a certain distance to the center of mass of group 1. That means it should give me the distance from OW to protein atom. If you choose correct groups that are