Re: [gmx-users] problem in energy minimization of the inflted system

2011-03-27 Thread Mark Abraham
On 28/03/11, Parul tew wrote: > Dear GMX users, > >  I am facing problem while doing energy minimization of the inflated system > of a protein with membrane. I inflated the system and updated the removed > lipids in the topology file. Now I am stuck during the energy minimization > state. B

[gmx-users] problem in energy minimization of the inflted system

2011-03-27 Thread Parul tew
Dear GMX users, I am facing problem while doing energy minimization of the inflated system of a protein with membrane. I inflated the system and updated the removed lipids in the topology file. Now I am stuck during the energy minimization state. Below is the .mdp file for minimization. -

Re: [gmx-users] g_clustsize and monomers

2011-03-27 Thread devicerandom
On 27/03/11 20:55, David van der Spoel wrote: On 2011-03-27 21.27, devicerandom wrote: On 27/03/11 16:29, David van der Spoel wrote: On 2011-03-27 16.42, devicerandom wrote: Hi, I am trying to use g_clustsize to follow aggregation of peptides. I am currently using the following line: g_clust

Re: [gmx-users] dihedral angles

2011-03-27 Thread Justin A. Lemkul
Sai Pooja wrote: Hi, I am running a simulation for capped alanine dipeptide and I now need to extract the dihedral angles. However, one of the nitrogen is named NT instead of NH due to which g_rama is not able to give the dihedral angles. Is there any other program I can use to get dihedra

[gmx-users] dihedral angles

2011-03-27 Thread Sai Pooja
Hi, I am running a simulation for capped alanine dipeptide and I now need to extract the dihedral angles. However, one of the nitrogen is named NT instead of NH due to which g_rama is not able to give the dihedral angles. Is there any other program I can use to get dihedral angles? I tried g_dih a

[gmx-users] Is there still interest in rigid-body simulation?

2011-03-27 Thread Adam Herbst
Hi all, I have seen a few posts on gmx-users indicating a desire to treat certain atom groups as rigid bodies in MD simulations. I just started implementing this, and so far I have it working for translational forces (not rotation, though this should be simple to add), even when the group is split

Re: [gmx-users] mathematical expression to obtain the C6 and C12 LJ terms for the GROMOSXX ff family

2011-03-27 Thread Justin A. Lemkul
sa wrote: Dear All, I would like to compute the C6 and C12 LJ terms for the GROMOS 53A6 ff for a new atom type and incorporate them in my ffnonbonded.itp. Does anybody know the mathematical expressions used to obtain these terms from a sigma (nm) and an epsilon (kj/mol) values manually.

Re: [gmx-users] g_clustsize and monomers

2011-03-27 Thread David van der Spoel
On 2011-03-27 21.27, devicerandom wrote: On 27/03/11 16:29, David van der Spoel wrote: On 2011-03-27 16.42, devicerandom wrote: Hi, I am trying to use g_clustsize to follow aggregation of peptides. I am currently using the following line: g_clustsize -f my.trr -s ../my.tpr -nc nclusters.xvg -

Re: [gmx-users] g_clustsize and monomers

2011-03-27 Thread devicerandom
On 27/03/11 16:29, David van der Spoel wrote: On 2011-03-27 16.42, devicerandom wrote: Hi, I am trying to use g_clustsize to follow aggregation of peptides. I am currently using the following line: g_clustsize -f my.trr -s ../my.tpr -nc nclusters.xvg -mol Problem is, it actually counts monome

Re: [gmx-users] heat capacity in 4.5.3 and 4.0.7

2011-03-27 Thread David van der Spoel
On 2011-03-27 19.23, Elisabeth wrote: Thanks David for replying me. I appreciate your attention to comment on the following short inquiries. 1- g_energy says: you will need to give the number of constraints per molecule –nconstr and for water this is 3

[gmx-users] GROMACS 4.5 and MARTINI Force Field Installation

2011-03-27 Thread Edroaldo Lummertz da Rocha
Dear GROMACS users, I am starting to setup a molecular dynamics simulation and I want to use the MARTINI force field. I was reading the tutorial in website: http://md.chem.rug.nl/cgmartini/index.php/tutorial/reverse-transformationand I have a beginner question for you: How can I use the MARTINI in

[gmx-users] Invitation to connect on LinkedIn

2011-03-27 Thread ram prasad via LinkedIn
LinkedIn ram prasad requested to add you as a connection on LinkedIn: -- Chinmay, I'd like to add you to my professional network on LinkedIn. - ram Accept invitation from ram prasad http://www.linkedin.com/e/-85v1n9-gls95twz-50/Xl9gjr6GAOlB4vI

[gmx-users] mathematical expression to obtain the C6 and C12 LJ terms for the GROMOSXX ff family

2011-03-27 Thread sa
Dear All, I would like to compute the C6 and C12 LJ terms for the GROMOS 53A6 ff for a new atom type and incorporate them in my ffnonbonded.itp. Does anybody know the mathematical expressions used to obtain these terms from a sigma (nm) and an epsilon (kj/mol) values manually. Thank you for your

Re: [gmx-users] heat capacity in 4.5.3 and 4.0.7

2011-03-27 Thread Elisabeth
> > Thanks David for replying me. I appreciate your attention to comment on >> the following short inquiries. >> >> 1- g_energy says: >> >> you will need to give the number of constraints per molecule –nconstr and > > for water this is 3. Can you elaborate on how to set this option. When constrai

Re: [gmx-users] g_clustsize and monomers

2011-03-27 Thread David van der Spoel
On 2011-03-27 16.42, devicerandom wrote: Hi, I am trying to use g_clustsize to follow aggregation of peptides. I am currently using the following line: g_clustsize -f my.trr -s ../my.tpr -nc nclusters.xvg -mol Problem is, it actually counts monomers as clusters. I have a system with 125 separa

[gmx-users] g_clustsize and monomers

2011-03-27 Thread devicerandom
Hi, I am trying to use g_clustsize to follow aggregation of peptides. I am currently using the following line: g_clustsize -f my.trr -s ../my.tpr -nc nclusters.xvg -mol Problem is, it actually counts monomers as clusters. I have a system with 125 separate molecules, and the number of cluster

Re: [gmx-users] RMSD Calculation

2011-03-27 Thread Tsjerk Wassenaar
Hey :) You probably want to fit on the protein and calculate the RMSD on the ligand. You may need to specify these groups in an index file. Hope it helps, Tsjerk On Mar 27, 2011 3:28 AM, "Justin A. Lemkul" wrote: Nancy wrote: > > Hi All, > > I need to determine the RMSD of a small molecule co

Re: [gmx-users] NADPH parameters

2011-03-27 Thread Mark Abraham
On 27/03/11, Kishore wrote: > Dear Mark/Justin, > You were asked yesterday to do things when you reply to digest emails. Do change the subject to something useful, and do not include the entire digest in your email. Not taking actions to play nicely with the other children will lead to the o

Re: [gmx-users] amide partial charges

2011-03-27 Thread Mark Abraham
On 27/03/11, XUEMING TANG wrote: > Hi > >   > > May I ask if anyone know a good reference for partial charges of > "-CH2-NH-CO-CH2-CH2"(amide)? I plan to use gromos 53A6 force field. > I searched the internet quite a while but not sure what is the best. I will > be appreciated if you coul

Re: [gmx-users] heat capacity in 4.5.3 and 4.0.7

2011-03-27 Thread David van der Spoel
On 2011-03-27 07.36, Elisabeth wrote: Dear all, I am trying to calculate Cv of a pure alkane 125 molecules with actual density, and here is what I have done: As indicated in an earlier post this is a non-solved problem, but we are working on it. The computation of heat capacity in g_energy w