[gmx-users] ensemble distance restraints

2013-01-11 Thread Da-Wei Li
Dear All I tried to run MD with ensemble based distance restrain. I set disre=Ensemble and use -multi command argument in mdrun, according to the manu. However, gromacs4.6 give me an error "Sorry, distance restraints with ensemble averaging over multiple molecules in one system are not functional

Re: [gmx-users] Reordering of a large GRO file - how

2013-01-11 Thread Tsjerk Wassenaar
Hey, Residue numbers don't matter. You can reorder the .gro file to have all the solvent at the end: sed -e '/SOL/{H;d}' -e '${x;s/^\n//p;x}' file.gro Note that this doesn't work with BSD sed, like on a Mac :p Cheers, Tsjerk On Thu, Jan 10, 2013 at 1:34 PM, Justin Lemkul wrote: > > > On 1/1

Re: [gmx-users] Covariance analysis of X-ray ensemble

2013-01-11 Thread Tsjerk Wassenaar
Hi James, The default is non-mass weighted. If you want to relate the principal components to normal modes, and want to obtain mode frequencies, you (formally) need to include the masses. Cheers, Tsjerk On Fri, Jan 11, 2013 at 7:46 PM, James Starlight wrote: > Tsjerk, thanks for suggestion! >

Re: [gmx-users] Covariance analysis of X-ray ensemble

2013-01-11 Thread James Starlight
Tsjerk, thanks for suggestion! By the way I've found in the pca.log ( that time pca was done on trajectory as well as TRP ( not pdb) files as the inputs ): Read 30 frames from cam_xray_coors_bb.pdb (time 1 to 30 ps) Read reference structure for fit from x_ray.tpr Analysis group is 'System' (575

[gmx-users] Cutoff: PME vs Shift

2013-01-11 Thread Juliette N.
Hello all, I am trying to compare PME and shift electrostatics methods for my hydrocarbon system. In both cases rlist = 1.35 , but in case of PME rcoulomb= 1.35 while for shift ** rcoulomb= 1.1. I see no significant change (less than 1%) in total system pot

Re: [gmx-users] Re: How to choose appropriate size for a sliding average for smoothing distances between molecule types

2013-01-11 Thread Justin Lemkul
On 1/11/13 7:18 AM, anna.duncan wrote: Thank you for your reply. You're right that I want to measure the protein-lipid interaction by determining the length of time that the lipid is within a certain distance of the protein. However, I also want to get a measure of the length of time that th

Re: [gmx-users] mdp option for implicit solvent simulation

2013-01-11 Thread Justin Lemkul
On 1/11/13 9:34 AM, Sanku M wrote: Hi Erik, That statement was a typo. The simulation will be done in implicit solvent ( not in gas phase). Sanku From: Erik Marklund To: Sanku M ; Discussion list for GROMACS users Sent: Friday, January 11, 2013 5:04 A

Re: [gmx-users] about g_hbond

2013-01-11 Thread Justin Lemkul
On 1/11/13 1:10 PM, Miguel Ángel Mompeán García wrote: The donors are two H atoms from the amide group of Gln residues. Thus, I will try the -merge option. Thank for the comment. I wonder if the following issue is affecting gromacs 4.5.5: http://lists.gromacs.org/pipermail/gmx-users/2011-May/

Re: [gmx-users] about g_hbond

2013-01-11 Thread Miguel Ángel Mompeán García
The donors are two H atoms from the amide group of Gln residues. Thus, I will try the -merge option. Thank for the comment. I wonder if the following issue is affecting gromacs 4.5.5: http://lists.gromacs.org/pipermail/gmx-users/2011-May/061249.html 2013/1/11, Justin Lemkul : > > > On 1/11/13 6:

Re: [gmx-users] Covariance analysis of X-ray ensemble

2013-01-11 Thread Tsjerk Wassenaar
Hi James, You don't need a .tpr file. A .gro or .pdb file will do as well. Unless you insist on doing mass-weighted PCA. Cheers, Tsjerk On Fri, Jan 11, 2013 at 4:04 PM, James Starlight wrote: > Dear Gromacs users! > > > I want to perform Covariance analysis of my x-ray data sets via > > g_cova

[gmx-users] Covariance analysis of X-ray ensemble

2013-01-11 Thread James Starlight
Dear Gromacs users! I want to perform Covariance analysis of my x-ray data sets via g_covar -f trr.pdb -s x_ray.tpr -n index -o PCA_eigenval.xvg -v PCA_eigenvec.trr where trr.pdb is the trajectory made from the 30 x-ray structures (only c-alpha atoms were selected. How I could obtain x_ray.tpr

Re: [gmx-users] mdp option for implicit solvent simulation

2013-01-11 Thread Sanku M
Hi Erik,   That statement was a typo. The simulation will be done in implicit solvent ( not in gas phase). Sanku From: Erik Marklund To: Sanku M ; Discussion list for GROMACS users Sent: Friday, January 11, 2013 5:04 AM Subject: Re: [gmx-users] mdp option fo

[gmx-users] about g_hbond

2013-01-11 Thread Miguel Ángel Mompeán García
Hi all, When calculating the average number of hbonds during a simulation, does the program g_hbond contemplate the possibility of an atom involved in two hbond? If I visualize the structures with molden, setting the same cutoff criterion for hbonding than g_hbond utilizes, I found out that some O

Re: [gmx-users] about g_hbond

2013-01-11 Thread Justin Lemkul
On 1/11/13 6:23 AM, Miguel Ángel Mompeán García wrote: Hi all, When calculating the average number of hbonds during a simulation, does the program g_hbond contemplate the possibility of an atom involved in two hbond? If I visualize the structures with molden, setting the same cutoff criterion

Re: [gmx-users] The number of lipids

2013-01-11 Thread Kieu Thu Nguyen
Thank Justin ! :-) On Fri, Jan 11, 2013 at 5:26 PM, Justin Lemkul wrote: > > > On 1/10/13 11:37 PM, Kieu Thu Nguyen wrote: > >> Thank Justin and Peter so much ! >> Can i make a new box (it's size is not a multiples of the size of the >> previous box) from a smaller lipid bilayer box ? >> >> I d

Re: [gmx-users] cann't calculate a silver particle

2013-01-11 Thread Justin Lemkul
On 1/11/13 6:33 AM, Alexej Mazheika wrote: Dear Justin, The grompp doesn't show any fatal error. Using gmxdump I have checked the topol.tpr file, and I found that C6 and C12 parameters of LJ have there completely different values in contrast to those I wrote in silver.itp! Maybe it is the reas

Re: [gmx-users] cann't calculate a silver particle

2013-01-11 Thread Alexej Mazheika
Dear Justin, The grompp doesn't show any fatal error. Using gmxdump I have checked the topol.tpr file, and I found that C6 and C12 parameters of LJ have there completely different values in contrast to those I wrote in silver.itp! Maybe it is the reason. But I wrote parameters in *.itp in units of

[gmx-users] about g_hbond

2013-01-11 Thread Miguel Ángel Mompeán García
Hi all, When calculating the average number of hbonds during a simulation, does the program g_hbond contemplate the possibility of an atom involved in two hbond? If I visualize the structures with molden, setting the same cutoff criterion for hbonding than g_hbond utilizes, I found out that some O

Re: [gmx-users] cann't calculate a silver particle

2013-01-11 Thread Justin Lemkul
On 1/11/13 5:02 AM, Alexej Mazheika wrote: Hello I removed PBC and set space conditions as you recommended. I also modified a bit the grompp.mdp (new version is below). However, it did not help to prevent the scattering of silver atoms in all directions. The change of parameters in LJ doesn't

Re: [gmx-users] RE: missing atom

2013-01-11 Thread Justin Lemkul
On 1/11/13 4:55 AM, Subramaniam Boopathi wrote: Hello All I have created a nwe residue "NIC" in AMBER 03 Forcrfield and I have added corresponding data to the aminoacid.rtp In the aminoacid.rtp, I have added all atoms topology information then I type the command

Re: [gmx-users] The number of lipids

2013-01-11 Thread Justin Lemkul
On 1/10/13 11:37 PM, Kieu Thu Nguyen wrote: Thank Justin and Peter so much ! Can i make a new box (it's size is not a multiples of the size of the previous box) from a smaller lipid bilayer box ? I don't see genconf -nbox can do that. genbox -cs old_lipids.gro -box x y z -o new_lipids.gro

Re: [gmx-users] mdp option for implicit solvent simulation

2013-01-11 Thread Erik Marklund
Hi, If you're simulating in the gas-phase, how come you want implicit solvent? Erik On Jan 11, 2013, at 5:07 AM, Sanku M wrote: Hi, I am trying to do implicit solvent simulation for a protein in gas- phase. I have a few questions : 1) should I use sd integrator? 2) should I also use n

Re: [gmx-users] cann't calculate a silver particle

2013-01-11 Thread Alexej Mazheika
Hello I removed PBC and set space conditions as you recommended. I also modified a bit the grompp.mdp (new version is below). However, it did not help to prevent the scattering of silver atoms in all directions. The change of parameters in LJ doesn't influence on the final geometry. It seems, that