Re: [gmx-users] Fe atom problems....in simulation

2012-04-16 Thread Kamalesh Roy
Can you send me the parameters I need to usefor FE On Mon, Apr 16, 2012 at 2:55 PM, francesco oteri wrote: > Hi, > You have to had LJ parameters in ffnonbonded.itp file in the subfolder > relative to your force-field > > Francesco > > Il giorno 16 aprile 2012 11:08, Kamalesh Roy ha > sc

Re: [gmx-users] Is GPU GTX 560 compatible with GROMACS?

2012-04-16 Thread Gaurav Goel
Dear Shouliang, Thanks. Maybe the developers can help undertand why GTX560 is not included on the GPU compatibility list for GROMACS? -G 2012/4/16 shouliang dong : > It is compatible with GROMACS. > > 在 2012年4月14日 上午11:32,Gaurav Goel 写道: >> >> Dear All, >> >> Can you please guide me on how to fi

[gmx-users] Re: temperature

2012-04-16 Thread Roy
Mark, According to the manual, the temperature of a group should always be calculated from the total KE of the group (see eq 3.13 in the manual), which presumably includes contributions from the translational motion of the center of mass of the group, plus the rotational motion around the center

Re: [gmx-users] What is epsilon_r?

2012-04-16 Thread Justin A. Lemkul
Andrew DeYoung wrote: Hi, What is the parameter epsilon_r mentioned in the manual (http://manual.gromacs.org/current/online/mdp_opt.html#el)? The manual says that it is the relative dielectric constant. My initial thought was that this would only be relevant for reaction-field electrostatics

[gmx-users] What is epsilon_r?

2012-04-16 Thread Andrew DeYoung
Hi, What is the parameter epsilon_r mentioned in the manual (http://manual.gromacs.org/current/online/mdp_opt.html#el)? The manual says that it is the relative dielectric constant. My initial thought was that this would only be relevant for reaction-field electrostatics, or somewhere where impli

Re: [gmx-users] ewald_rtol

2012-04-16 Thread Justin A. Lemkul
Andrew DeYoung wrote: Hi, Does the parameter ewald_rtol affect PME electrostatics (coulombtype = PME), or does ewald_rtol only affect Ewald electrostatics (coulombtype = Ewald)? In the manual description of the .mdp parameters (http://manual.gromacs.org/current/online/mdp_opt.html#el), it says

[gmx-users] ewald_rtol

2012-04-16 Thread Andrew DeYoung
Hi, Does the parameter ewald_rtol affect PME electrostatics (coulombtype = PME), or does ewald_rtol only affect Ewald electrostatics (coulombtype = Ewald)? In the manual description of the .mdp parameters (http://manual.gromacs.org/current/online/mdp_opt.html#el), it says in the Ewald section, "Th

Re: [gmx-users] Re: NPT simulation for mutation study (Justin A. Lemkul)

2012-04-16 Thread Justin A. Lemkul
Shyno Mathew wrote: Hey Justin, thanks for your reply. For gen_vel, the manual says "Generate velocities in grompp according to a Maxwell distribution at temperature gen_temp [K], with random seed gen_seed. This is only meaningful with integrator md" Since the integrator I am using is sd, I

Re: [gmx-users] GROMOS87 and CHARMM27

2012-04-16 Thread Justin A. Lemkul
Shima Arasteh wrote: So, I can not use the coordinates of the output files of gromos runs. Right? Not without significant modification of names, presence of H atoms, etc. You'll need to regenerate a suitable topology, as has been said, and then run thorough equilibration under the new forc

Re: [gmx-users] GROMOS87 and CHARMM27

2012-04-16 Thread Peter C. Lai
On 2012-04-16 10:14:01AM -0700, Shima Arasteh wrote: > So, I can not use the coordinates of the output files of gromos runs. Right? You can but you may need to rename the atoms for each residue for pdb2gmx to work. > > > > From: Peter C. Lai > To: Discussion l

[gmx-users] Re: NPT simulation for mutation study (Justin A. Lemkul)

2012-04-16 Thread Shyno Mathew
Hey Justin, thanks for your reply. For gen_vel, the manual says "Generate velocities in grompp according to a Maxwell distribution at temperature gen_temp [K], with random seed gen_seed. This is only meaningful with integrator md" Since the integrator I am using is sd, I put 'gen_vel no' Yes you a

Re: [gmx-users] GROMOS87 and CHARMM27

2012-04-16 Thread Shima Arasteh
So, I can not use the coordinates of the output files of gromos runs. Right? From: Peter C. Lai To: Discussion list for GROMACS users Sent: Monday, April 16, 2012 8:23 PM Subject: Re: [gmx-users] GROMOS87 and CHARMM27 On 2012-04-16 08:26:00AM -0700, Shima Ar

Re: [gmx-users] Re: Re: NPT simulation for mutation study

2012-04-16 Thread Justin A. Lemkul
Shyno Mathew wrote: Hello Prof. David, thanks for your reply. I will try using isotropic pressure scaling. But still I am not clear why the simulations ran fine with first .mdp file. As mentioned in the previous email only few parameters were different in my .mdp file compared to my advisor

[gmx-users] Re: Re: NPT simulation for mutation study

2012-04-16 Thread Shyno Mathew
Hello Prof. David, thanks for your reply. I will try using isotropic pressure scaling. But still I am not clear why the simulations ran fine with first .mdp file. As mentioned in the previous email only few parameters were different in my .mdp file compared to my advisor's. I am just copying those

Re: [gmx-users] Protein ligand molecular dynamics simulation

2012-04-16 Thread lina
On Mon, Apr 16, 2012 at 11:01 PM, sai nitin wrote: > Hi all, > > I have done complex (protein + ligand) complex from autodock software using > this complex im trying to follow > > http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin/gmx-tutorials/complex/01_pdb2gmx.html > tutorial > > But when

Re: [gmx-users] Tabulated potential segmentation fault

2012-04-16 Thread Mark Abraham
On 17/04/2012 1:50 AM, Laura Leay wrote: On Fri, 2012-04-13 at 11:01 +1000, Mark Abraham wrote: On 13/04/2012 2:48 AM, Laura Leay wrote: All, I'm trying to run a tabulated soft core potential with the form V = A + Br^2 + Cr^3 up to about r=0.1 A and the normal LJ 6-12 potential after this. I'

Re: [gmx-users] GROMOS87 and CHARMM27

2012-04-16 Thread Peter C. Lai
On 2012-04-16 08:26:00AM -0700, Shima Arasteh wrote: > Dear GROMACS users, > > I reproduced the results of a protein-membrane system by using force field > GROMOSE87. This protein forms ion channel in membrane. > Now if I wanna study the ion conduction through this channel using force > field CH

Re: [gmx-users] Tabulated potential segmentation fault

2012-04-16 Thread Laura Leay
On Fri, 2012-04-13 at 11:01 +1000, Mark Abraham wrote: > On 13/04/2012 2:48 AM, Laura Leay wrote: > > All, > > > > I'm trying to run a tabulated soft core potential with the form V = > > A + Br^2 + Cr^3 up to about r=0.1 A and the normal LJ 6-12 potential > > after this. > > > > I've chosen the

Re: [gmx-users] NPT simulation for mutation study

2012-04-16 Thread David van der Spoel
On 2012-04-16 17:23, Shyno Mathew wrote: Dear Gromcas users, I am doing some mutation study, NPT simulations. Initially I was using a generic .mdp file I got from my advisor and I was able to run the systems for ~50ns with out any issues. But then I spend time reading about gromacs and created a

[gmx-users] GROMOS87 and CHARMM27

2012-04-16 Thread Shima Arasteh
Dear GROMACS users, I reproduced the results of a protein-membrane system by using force field GROMOSE87. This protein forms ion channel in membrane. Now if I wanna study the ion conduction through this channel using force field CHARMM27 in umbrella sampling method, is it possible? Can I use the

[gmx-users] NPT simulation for mutation study

2012-04-16 Thread Shyno Mathew
Dear Gromcas users, I am doing some mutation study, NPT simulations. Initially I was using a generic .mdp file I got from my advisor and I was able to run the systems for ~50ns with out any issues. But then I spend time reading about gromacs and created a .mdp file. However, with my .mdp file the r

Re: [gmx-users] protein and DMPC in Charmm36 ff

2012-04-16 Thread Krzysztof Kuczera
An easy way to build a protein in a bilayer is through the charmm-gui website www.charmm-gui.org/ Partial drawbacks are that you need CHARMM to perform final equilibration and then may need to rename some atoms to work with GROMACS On plus side is that charmm36.ff retains the original CHARMM atom

Re: [gmx-users] Charmm27.ff with FEP

2012-04-16 Thread Mark Abraham
On 17/04/2012 1:14 AM, Sanku M wrote: Hi, Is there any known issue/problem in running FEP calculations with charm27.ff in gromacs4.5.4 ? I tried running an FEP calculation using charmm27.ff by interpolating A state and B state but it gives error that dihedral terms with multiple values can no

[gmx-users] Charmm27.ff with FEP

2012-04-16 Thread Sanku M
Hi,  Is there any known issue/problem in running FEP calculations with charm27.ff in gromacs4.5.4 ? I tried running an FEP calculation using charmm27.ff by interpolating A state and B state but it gives error that dihedral terms with multiple values can not be interpolated..One need to write all

Re: [gmx-users] Protein ligand molecular dynamics simulation

2012-04-16 Thread Justin A. Lemkul
sai nitin wrote: Hi all, I have done complex (protein + ligand) complex from autodock software using this complex im trying to follow http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin/gmx-tutorials/complex/01_pdb2gmx.html tutorial But when i take complex structure directly from a

[gmx-users] Protein ligand molecular dynamics simulation

2012-04-16 Thread sai nitin
Hi all, I have done complex (protein + ligand) complex from autodock software using this complex im trying to follow http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin/gmx-tutorials/complex/01_pdb2gmx.html tutorial But when i take complex structure directly from autodock result and run PDB

Re: [gmx-users] how to install DSSP in linux with Gromacs 4.0.7

2012-04-16 Thread Erik Marklund
16 apr 2012 kl. 13.32 skrev Justin A. Lemkul: > > > 李 麗花 wrote: >> Hi: >> I am getting upsep to install DSSP in linux with Gromacs 4.0.7 >> I do not know where is the problem >> which DSSP file I should download it >> and possible because gromacs 4.0.7 ? >> could someone can help me to figure o

Re: [gmx-users] Atom MNZ1 in residue LYS 107 was not found

2012-04-16 Thread Steven Neumann
Thankl you. Well... Indeed, just removed them manually and everything is ok. Steven On Mon, Apr 16, 2012 at 2:38 PM, Peter C. Lai wrote: > Sounds like a bad VMD selection that didn't remove MNZ1 when you did the > conversion. MNZ1 sounds like a virtual site for the extra hydrogen off NZ > in pr

Re: [gmx-users] Re: Simulation in the high temperature conditions

2012-04-16 Thread Justin A. Lemkul
James Starlight wrote: Justin, Thank you for explanation. Tomorrow I'll try to check results of simulation with the disres applied with its default values as well as with narrower -disre_dist values ( ignorring -disre_frac option at all ) and post here results of such simulations. 1) T

Re: [gmx-users] Re: Simulation in the high temperature conditions

2012-04-16 Thread James Starlight
Justin, Thank you for explanation. Tomorrow I'll try to check results of simulation with the disres applied with its default values as well as with narrower -disre_dist values ( ignorring -disre_frac option at all ) and post here results of such simulations. 1) The cut-off distance wich I've s

Re: [gmx-users] Atom MNZ1 in residue LYS 107 was not found

2012-04-16 Thread Peter C. Lai
Sounds like a bad VMD selection that didn't remove MNZ1 when you did the conversion. MNZ1 sounds like a virtual site for the extra hydrogen off NZ in protonated LYS. On 2012-04-16 02:32:09PM +0100, Steven Neumann wrote: > Dear Gmx Users, > > I run implicit simulation for 1 us with virtual sites o

Re: [gmx-users] Re: Simulation in the high temperature conditions

2012-04-16 Thread Justin A. Lemkul
James Starlight wrote: Justin, I've applied disres on each backbone atom of my potein within cutoff distance of 1nm ( Rc=1.0 nm). I've selected this value for cutoff to decrease overall ammount of the restains in my itp file. Also such value ( 1nm) was selected because of the relatively ti

[gmx-users] Re: pre-wham

2012-04-16 Thread lloyd riggs
Dear all, Im using wham on series of runs and have run into the following proble. The command line >g_wham_d -if pullf.dat -it tpr.dat -temp 300 -o whamf.xvg -hist outf.hist >-unit kCal will work in one run, but not the next one as I go through I found a couple that give the error Fatal er

[gmx-users] Atom MNZ1 in residue LYS 107 was not found

2012-04-16 Thread Steven Neumann
Dear Gmx Users, I run implicit simulation for 1 us with virtual sites on hydrogens, then using VMD extracted coordinates into the pdb file. As I want to run explicit solvent simulation now I removed hydrogens so that pdb2gmx will add them. Then I got an error while trying to pdb2gmx using Charmm27

Re: [gmx-users] Re: Simulation in the high temperature conditions

2012-04-16 Thread James Starlight
Justin, I've applied disres on each backbone atom of my potein within cutoff distance of 1nm ( Rc=1.0 nm). I've selected this value for cutoff to decrease overall ammount of the restains in my itp file. Also such value ( 1nm) was selected because of the relatively tight packing of the alpha helic

Re: [gmx-users] DD load balancing is limited by minimum cell size in dimension Z

2012-04-16 Thread Mark Abraham
On 16/04/2012 8:43 PM, Gavin Melaugh wrote: Hi all What does the following note mean in the log file DD load balancing is limited by minimum cell size in dimension Z Is is purely a performance related issue ? Yes. See 3.17.2 for description of DLB. Various algorithms constrain minimum cell

[gmx-users] Clustering

2012-04-16 Thread dina dusti
Hi Erik, Thank you very much from your response. Best Regards Dina From: Erik Marklund To: dina dusti ; Discussion list for GROMACS users Sent: Monday, April 16, 2012 12:07 PM Subject: Re: [gmx-users] Clustering Hi, It's been a while since I used g_clu

Re: [gmx-users] how to install DSSP in linux with Gromacs 4.0.7

2012-04-16 Thread Justin A. Lemkul
李 麗花 wrote: Hi: I am getting upsep to install DSSP in linux with Gromacs 4.0.7 I do not know where is the problem which DSSP file I should download it and possible because gromacs 4.0.7 ? could someone can help me to figure out it ? many thanks http://swift.cmbi.ru.nl/gv/dssp/ Under "Misce

[gmx-users] how to install DSSP in linux with Gromacs 4.0.7

2012-04-16 Thread 李 麗花
Hi:I am getting upsep to install DSSP in linux with Gromacs 4.0.7I do not know where is the problemwhich DSSP file I should download itand possible because gromacs 4.0.7 ?could someone can help me to figure out it ?many thanks Best Wishesli-hua -- gmx-us

Re: [gmx-users] Re: Simulation in the high temperature conditions

2012-04-16 Thread Justin A. Lemkul
James Starlight wrote: Dear Gromacs Users! By that moments I've completed 2 sets of simulation in high temperature 1- With applied posres on the backbone atoms ( fc= 200 ). The result was- that the posres prevented motion of the helixes as the rigid bodies so I've not noticed any conformati

[gmx-users] DD load balancing is limited by minimum cell size in dimension Z

2012-04-16 Thread Gavin Melaugh
Hi all What does the following note mean in the log file DD load balancing is limited by minimum cell size in dimension Z Is is purely a performance related issue ? Cheers Gavin -- gmx-users mailing listgmx-users@gromacs.org http://lists.gromacs.org/mailman/listinfo/gmx-users Please sear

Re: [gmx-users] Fe atom problems....in simulation

2012-04-16 Thread francesco oteri
Hi, You have to had LJ parameters in ffnonbonded.itp file in the subfolder relative to your force-field Francesco Il giorno 16 aprile 2012 11:08, Kamalesh Roy ha scritto: > CAn any one siggest me how can I run simulation a protein containing Fe > atom, > I have changed the iions.itpfile and inc

[gmx-users] Fe atom problems....in simulation

2012-04-16 Thread Kamalesh Roy
CAn any one siggest me how can I run simulation a protein containing Fe atom, I have changed the iions.itpfile and included there Fe in residue type.dat bu still it is returning an error. FE parameter not found. -- gmx-users mailing listgmx-users@gromacs.org http://lists.gromacs.org/mailman/l

[gmx-users] Re:Data Analysis

2012-04-16 Thread lloyd riggs
Dear All, Is there a way to extract the pull force (and not the total force or energy) for a specified group (indexed) or center of mass, rather than just the total pullf over a run? I wanted to break up the pullf into varied sub contributions, however as it is applied to the center of mass o

Re: [gmx-users] Clustering

2012-04-16 Thread Erik Marklund
Hi, It's been a while since I used g_clustsize, but if I'm not mistaken you get the number of molecules per cluster as output, from which you can estimate the physical size of the clusters with a few assumptions. If that's not good enough you also get the identities of the molecules in the clus

Re: [gmx-users] Re: Simulation in the high temperature conditions

2012-04-16 Thread James Starlight
Dear Gromacs Users! By that moments I've completed 2 sets of simulation in high temperature 1- With applied posres on the backbone atoms ( fc= 200 ). The result was- that the posres prevented motion of the helixes as the rigid bodies so I've not noticed any conformation sampling. Question : Cou