[gmx-users] To know

2012-05-02 Thread Anik Sen
Hi, Am Anik Sen, using gromacs 3.3.2 version. I have a very simple question to ask, In the mdp file, like the emmod.mdp or run.mdp needed for a MD run there are three factors namely, rlist, rcoulomb and rvdw. 1. Is there any minimum value for these three below which the results will be

Re: [gmx-users] To know

2012-05-02 Thread Mark Abraham
On 2/05/2012 5:12 PM, Anik Sen wrote: Hi, Am Anik Sen, using gromacs 3.3.2 version. I have a very simple question to ask, In the mdp file, like the emmod.mdp or run.mdp needed for a MD run there are three factors namely,rlist, rcoulomb and rvdw. 1. Is there any minimum value for these

Re: [gmx-users] To know

2012-05-02 Thread Mark Abraham
On 2/05/2012 5:37 PM, Mark Abraham wrote: On 2/05/2012 5:12 PM, Anik Sen wrote: Hi, Am Anik Sen, using gromacs 3.3.2 version. I have a very simple question to ask, In the mdp file, like the emmod.mdp or run.mdp needed for a MD run there are three factors namely,rlist, rcoulomb and

[gmx-users] Automation of selecting water around a molecule

2012-05-02 Thread Lara Bunte
Hi I need a little help in using Linux. I wrote this script for the bash, that I called g_select_script.sh: #!/bin/bash g_select -s mol_in_water.pdb -select 'Close to ISO resname SOL and within 0.4 of resname ISO' -on make_ndx -f mol_in_water.pdb -o index2.ndx cat index2.ndx index.ndx

[gmx-users] some hardware questions

2012-05-02 Thread melichercik
Hi people, we plan to buy some computer mostly for gromacs, so I want to ask you for some hints. I'd tried to do some research, but no success. :( Have someone ever compared AMD Bulldozer and Sandy/Ivy Bridge architectures? Are there planed any optimizations in gmx core in asm parts? (if

Re: [gmx-users] File editing - only one layer of water around a molecule

2012-05-02 Thread Teemu Murtola
On Tue, May 1, 2012 at 9:55 PM, Justin A. Lemkul jalem...@vt.edu wrote: On 5/1/12 2:40 PM, Lara Bunte wrote: g_select -s molecule.pdb -select 'Close to ISO resname SOL and within 1.2 of resname ISO' and trjconv -n index.ndx -f molecule.pdb -s molecule.pdb -o output.pdb I got only water in

Re: [gmx-users] some hardware questions

2012-05-02 Thread Thomas Evangelidis
For CPU performance comparison look at: www.spec.org/cpu2006/results for each CPU look at the statistics for floating point operations, namely for SPECfp_base2006 score for the programs of interest (GAMESS, GROMACS, NAMD, etc.). For graphics cards you can look at the same site for GPU

[gmx-users] Re: Re: g_wham problem with negative COM differences

2012-05-02 Thread Anni Kauko
Anni Kauko wrote: Hi! I try to perform pmf calculations for case where a peptide shifts through the membrane. My COM differences should vary from 2.3 to -2.5. My problem is that g_wham plots negative COM difference as

[gmx-users] How to remove H atom from residue in gro file?

2012-05-02 Thread Hagit G
Hi gmx users, Well, I saw this question but the answer was not understood. I'm trying to work with the file 1PPB.pdb. There are 2 chains connected with a disulfide bond. Gromacs automatically adds H atoms. Although the disulfide bond is there, Gromacs ignore it because *each cystein is on a

[gmx-users] Number of coordinate in topol.top and solv.gro not matching

2012-05-02 Thread Sangita Kachhap
Hello all I have to do MD simulation of ligand bound membrane protein in lipid bilayer. Thus I am doing tutorialfor Protein - Ligand tutorial: http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin/gmx-tutorials/complex/index.html I have used following command: pdb2gmx -f 3HTB_clean.pdb -o

[gmx-users] CHARMM36 format POPC

2012-05-02 Thread Anirban Ghosh
Hi ALL, I was looking for a CHARMM36 format (atom-types) equilibrated POPC bilayer (PDB/GRO) to use with the CHARMM36 FF in GROAMCS 4.5.5. I downloaded one from Dr. Klauda's site ( http://terpconnect.umd.edu/~jbklauda/research/download.html), but that popc.pdb (under CHARMM36 FF) seem to have

[gmx-users] Question about tabulated potentials

2012-05-02 Thread Marcelo Lopez
Hi all, I want to be sure to understand what I'm doing with tabulated potentials. I'm using tabulated bonds, angles, dihedrals and non-bnding potentials. 1) If I don't specify 1-4 interactions (nrexcl = 3 and gen-pairs = no), they're all zero for the 1-4 atoms involved in tabulated dihedrals?

Re: [gmx-users] GROMOS87 and CHARMM27

2012-05-02 Thread Justin A. Lemkul
On 5/2/12 12:21 AM, Shima Arasteh wrote: Dear Justin, I don't want to compare the forcefields results. Yes, Iwant to use some old information in conjunction with something new. In fact, I have to regenerate a structure, which has modulated using gmx.ff. After it, I am going to study it using

Re: [gmx-users] Re: Re: g_wham problem with negative COM differences

2012-05-02 Thread Justin A. Lemkul
On 5/2/12 6:49 AM, Anni Kauko wrote: Anni Kauko wrote: Hi! I try to perform pmf calculations for case where a peptide shifts through the membrane. My COM differences should vary from 2.3 to -2.5. My

Re: [gmx-users] How to remove H atom from residue in gro file?

2012-05-02 Thread Justin A. Lemkul
On 5/2/12 6:55 AM, Hagit G wrote: Hi gmx users, Well, I saw this question but the answer was not understood. I'm trying to work with the file 1PPB.pdb. There are 2 chains connected with a disulfide bond. Gromacs automatically adds H atoms. Although the disulfide bond is there, Gromacs ignore

Re: [gmx-users] How to remove H atom from residue in gro file?

2012-05-02 Thread Mark Abraham
On 2/05/2012 8:55 PM, Hagit G wrote: Hi gmx users, Well, I saw this question but the answer was not understood. I'm trying to work with the file 1PPB.pdb. There are 2 chains connected with a disulfide bond. Gromacs automatically adds H atoms. Although the disulfide bond is there, Gromacs

Re: [gmx-users] Number of coordinate in topol.top and solv.gro not matching

2012-05-02 Thread Justin A. Lemkul
On 5/2/12 8:14 AM, Sangita Kachhap wrote: Hello all I have to do MD simulation of ligand bound membrane protein in lipid bilayer. Thus I am doing tutorialfor Protein - Ligand tutorial: http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin/gmx-tutorials/complex/index.html I have used

Re: [gmx-users] Automation of selecting water around a molecule

2012-05-02 Thread Mark Abraham
On 2/05/2012 7:49 PM, Lara Bunte wrote: Hi I need a little help in using Linux. I wrote this script for the bash, that I called g_select_script.sh: #!/bin/bash g_select -s mol_in_water.pdb -select 'Close to ISO resname SOL and within 0.4 of resname ISO' -on make_ndx -f mol_in_water.pdb -o

[gmx-users] سلام استاد

2012-05-02 Thread Shima Arasteh
سلام یه خواهش دارم. من بعد قاب کردن نقاشی، ازش عکس نگرفتم. لطفا اگر تونستین برام چندتا عکس خوب بگیرین مرسی میشه لطفا؟-- gmx-users mailing listgmx-users@gromacs.org http://lists.gromacs.org/mailman/listinfo/gmx-users Please search the archive at

Re: [gmx-users] GROMOS87 and CHARMM27

2012-05-02 Thread Justin A. Lemkul
On 5/2/12 10:30 AM, Shima Arasteh wrote: Thanks for your reply. What do you mean by running for some additional time? You're introducing a new force field, which includes different representation of atoms and a different potential energy function. You can't just take a configuration

[gmx-users] The wrong email, Sorry

2012-05-02 Thread Shima Arasteh
Dear gmx users, I am so sorry to send a wrong email. My last mail sent to mailing list sent by mistake. Sorry Shima -- gmx-users mailing listgmx-users@gromacs.org http://lists.gromacs.org/mailman/listinfo/gmx-users Please search the archive at

Re: [gmx-users] GROMOS87 and CHARMM27

2012-05-02 Thread Shima Arasteh
Thanks for your kind reply . Sincerely, Shima From: Justin A. Lemkul jalem...@vt.edu To: Discussion list for GROMACS users gmx-users@gromacs.org Sent: Wednesday, May 2, 2012 7:02 PM Subject: Re: [gmx-users] GROMOS87 and CHARMM27 On 5/2/12 10:30 AM, Shima

Re: [gmx-users] سلام استاد

2012-05-02 Thread leila separdar
salam, man mishnasametun? mozoe naghashie chiye? 2012/5/2 Shima Arasteh shima_arasteh2...@yahoo.com سلام یه خواهش دارم. من بعد قاب کردن نقاشی، ازش عکس نگرفتم. لطفا اگر تونستین برام چندتا عکس خوب بگیرین مرسی میشه لطفا؟ -- gmx-users mailing listgmx-users@gromacs.org

Re: [gmx-users] some hardware questions

2012-05-02 Thread Peter C. Lai
When considering compute hardware today, you're gonna have to shop like a gamer. A rig that can get high Battlefield 3 or Arma 2/3 framerates should also be able to get you some nice ns/day. I think most people consider Bulldozer mostly a flop. AMD masks lower cycle speeds with more cores per

Re: [gmx-users] CHARMM36 format POPC

2012-05-02 Thread Peter C. Lai
http://uab.hyperfine.info/~pcl/files/popc36/ These were generated for the following work: (please reference): Lai, P.C. and Crasto, C.J. Beyond Modeling: All-Atom Olfactory Receptor Model Simulations Front. Gene. Volume 3 Year 2012 Number 61 DOI: 10.3389/fgene.2012.00061 On 2012-05-02

Re: [gmx-users] some hardware questions

2012-05-02 Thread Mirco Wahab
Hello Peter, Am 02.05.2012 17:44, schrieb Peter C. Lai: You can wait for ivy bridge, then stick some Kepler GPUs (nvidia gtx 680) in it. That should max the performance. asm is pretty much stagnant for general purpose procs since core2 came out. Is this true? To my knowledge, the Fermi GPU

Re: [gmx-users] some hardware questions

2012-05-02 Thread Peter C. Lai
On 2012-05-02 06:13:04PM +0200, Mirco Wahab wrote: Hello Peter, Am 02.05.2012 17:44, schrieb Peter C. Lai: You can wait for ivy bridge, then stick some Kepler GPUs (nvidia gtx 680) in it. That should max the performance. asm is pretty much stagnant for general purpose procs since core2

[gmx-users] Number of cooredinate in topol.top and solv.gro not matching

2012-05-02 Thread Sangita Kachhap
/Making_Disulfide_Bonds Mark -- next part -- An HTML attachment was scrubbed... URL: http://lists.gromacs.org/pipermail/gmx-users/attachments/20120502/9c7e1646/attachment-0001.html -- Message: 3 Date: Wed, 02 May 2012 09:12:13 -0400 From

Re: [gmx-users] doubt about OpenMM

2012-05-02 Thread Martin Hoefling
Am 30.04.2012 19:32, schrieb Thales Kronenberger: Hey Folks, SUDO_USER=laboratorio24 PWD=/home/laboratorio24/OpenMM4.0-Linux64/python File setup.py, line 165 print %s % outputString ^ SyntaxError: invalid syntax INSTALLATION FAILED which python version are you using?

Re: [gmx-users] doubt about OpenMM

2012-05-02 Thread Thales Kronenberger
I have Python 2.6.5 2012/5/2 Martin Hoefling martin.hoefl...@googlemail.com Am 30.04.2012 19:32, schrieb Thales Kronenberger: Hey Folks, SUDO_USER=laboratorio24 PWD=/home/laboratorio24/OpenMM4.0-Linux64/python File setup.py, line 165 print %s % outputString ^

[gmx-users] pullf.xvg

2012-05-02 Thread Krzysztof Kuczera
Could someone please tell me what exactly is printed in the pullf.xvg file from an umbrella sampling simulation? Is it the total force acting the COM of the pull group - i.e. force from regular potential plus the umbrella force or just the unbiased force? Krzysztof Kuczera -- Krzysztof

[gmx-users] Pull code for Constraint

2012-05-02 Thread xu zhijun
Hi, I try to pull the small molecule from the bulk solution to the surface by the Constraint function. However, it seems that it does not work well. Here my part of run parameters about the Constraint in my mdp file ## constraints = all-bonds; hbonds; all-bonds

Re: [gmx-users] Number of cooredinate in topol.top and solv.gro not matching

2012-05-02 Thread Justin A. Lemkul
On 5/2/12 1:09 PM, Sangita Kachhap wrote: I have to do MD simulation of ligand bound membrane protein in lipid bilayer. Thus I am doing tutorialfor Protein - Ligand tutorial: http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin/gmx-tutorials/complex/index.html I have used following

RE: [gmx-users] pullf.xvg

2012-05-02 Thread Marzinek, Jan
From: gmx-users-boun...@gromacs.org [gmx-users-boun...@gromacs.org] on behalf of Krzysztof Kuczera [kkucz...@ku.edu] Sent: Wednesday, May 02, 2012 6:38 PM To: Discussion list for GROMACS users Subject: [gmx-users] pullf.xvg Could someone please tell me

[gmx-users] Pull direction is always outward

2012-05-02 Thread cescoxonta
I am trying to pull a subdomain of protein towards the center of the channel it makes. To do this I need to pull some aminoacids in -1, 0, 0 direction. I tried with this pull code ; Pull code pull= umbrella pull_geometry = direction_periodic pull_start = yes pull_nstfout= 1

Re: [gmx-users] NBFIX to gromacs

2012-05-02 Thread Ricardo O. S. Soares
Dear Mr. Peter C. Lai, I appreciate your kind response. I apologize for the delay on my part. This NBFIX parameters that I'm working with came from the str file in http://terpconnect.umd.edu/~jbklauda/research/download/chol_new.str , which dates from 2012. I successfully converted this whole

[gmx-users] Re: Pull code for Constraint

2012-05-02 Thread cescoxonta
pull_k1 is missing in your input, so it is probably set to 0. I do not know whether it is used in constraint mode, though. -- View this message in context: http://gromacs.5086.n6.nabble.com/Pull-code-for-Constraint-tp4946900p4947436.html Sent from the GROMACS Users Forum mailing list archive at

Re: [gmx-users] NBFIX to gromacs

2012-05-02 Thread Peter C. Lai
The problem is that functype 1 for [ pairtypes ] when gromacs uses a CHARMM forcefield is that gromacs treats them as 1-4 interactions. You will have to ask/wait for someone familiar with the gromacs source to say if there is a different functype available to use to override non-1-4 interactions

[gmx-users] trjconv correct the pbc before analysis

2012-05-02 Thread mu xiaojia
Dear gmx users, I have a question about using the trjconv -pbc options before analyzing my trajectory. It's stated by Justin's tutorial that: use trjconv to account for any periodicity in the system. The protein will diffuse through the unit cell, and may appear to jump across to the other side

Re: [gmx-users] trjconv correct the pbc before analysis

2012-05-02 Thread Justin A. Lemkul
On 5/2/12 8:55 PM, mu xiaojia wrote: Dear gmx users, I have a question about using the trjconv -pbc options before analyzing my trajectory. It's stated by Justin's tutorial that: use trjconv to account for any periodicity in the system. The protein will diffuse through the unit cell, and may

Re: [gmx-users] trjconv correct the pbc before analysis

2012-05-02 Thread Mark Abraham
On 03/05/12, mu xiaojia muxiaojia2...@gmail.com wrote: Dear gmx users, I have a question about using the trjconv -pbc options before analyzing my trajectory. It's stated by Justin's tutorial that: use trjconv to account for any periodicity in the system. The protein will

[gmx-users] Nubmer of coordinate not matching in topol.top and sol.gro

2012-05-02 Thread Sangita Kachhap
the www interface or send it to gmx-users-requ...@gromacs.org. Can't post? Read http://www.gromacs.org/Support/Mailing_Lists -- next part -- An HTML attachment was scrubbed... URL: http://lists.gromacs.org/pipermail/gmx-users/attachments/20120502/676d2a14/attachment

Re: [gmx-users] CHARMM36 format POPC

2012-05-02 Thread Anirban Ghosh
Thanks a lot Peter. Will try out with this. -Anirban On Wed, May 2, 2012 at 9:42 PM, Peter C. Lai p...@uab.edu wrote: http://uab.hyperfine.info/~pcl/files/popc36/ These were generated for the following work: (please reference): Lai, P.C. and Crasto, C.J. Beyond Modeling: All-Atom Olfactory

[gmx-users] a question about ensemble

2012-05-02 Thread Albert
hello: I wondering are the three thermostat methods: Langevin, Berendsen and Nose-Hoover chain are all compatible with semi-isotropy coupling style? If I would like to use semi-isotropy coupling method, which one would be better? thank you very much best Albert -- gmx-users mailing list