[gmx-users] Re: A common error Atom O11 in residue bDM 1 was

2010-09-27 Thread ABEL Stephane 175950
OK you thank Justin for your explanation. Accordingly, I have changed the C* and O* atom names in the pdb and now it works. A bientôt Stefane ABEL Stephane 175950 wrote: > You are right Justin, > > The 1O1 atoms and (others 1O2, etc.) are changed to O11, O21, etc. (?). Why > this problem h

Re: [gmx-users] Query regarding protonation and deprotonation of some residues

2010-09-27 Thread Tsjerk Wassenaar
Hi Sonali, First of all, you'll have to find a force field that supports tyrosinate and serinate. These aren't generally considered titratable and are consequently not in the list for interactive selections with pdb2gmx. Once you've found a suitable force field, you'll have to set the protonation

Re: [gmx-users] Reg: Adding multiple molecules

2010-09-27 Thread Justin A. Lemkul
vinothkumar mohanakrishnan wrote: Hi all I have a diffculty in adding multiple molecules of hexane to my box. my box size is 4.72*2.36*2.36 (nm) according to my number density calculations the box should fit 124 molecules of hexane whereas it adds only 76 molecules to the box below is my co

Re: [gmx-users] How to make a lipid bilayer with specific dimensions?

2010-09-27 Thread Justin A. Lemkul
NG HUI WEN wrote: Dear gmxusers, I am trying to make a lipid bilayer with specific dimensions using gromacs. So far, I have got up to: 1) Download a lipid POPC128a.pdb from Peter Tieleman’s website 2) Use genconf –f popc128a.pdb –o popcx2.pdb –nbox 2 2 1 to multiply t

[gmx-users] Re: Ask for a question

2010-09-27 Thread Justin A. Lemkul
Please keep all Gromacs-related correspondence on the gmx-users list. I am not a private tutor. Most of your questions can be answered by reading molecular modeling textbooks or review articles about pulling. Some of this information is also in the manual, so you've certainly missed it.

Re: [gmx-users] Reg: Adding multiple molecules

2010-09-27 Thread vinothkumar mohanakrishnan
Hi justin if genbox is not able to find sufficient space then why should one calculate number density for the given box dimension?.similarly the the number of molecules are not getting updated in the molecules section of the topology file when inserting multiple molecules using the genbox command

Re: [gmx-users] Reg: Adding multiple molecules

2010-09-27 Thread Justin A. Lemkul
vinothkumar mohanakrishnan wrote: Hi justin if genbox is not able to find sufficient space then why should one calculate number density for the given box dimension?.similarly the the Just because you can calculate a certain number density does not necessarily mean that the algorithm genb

[gmx-users] Pulling Simulation Query

2010-09-27 Thread Natalie Stephenson
Hi ... Just wondering if there's a way, during pulling simulations, to immobilise say the C-terminus of the protein. I've been running through the 'Umbrella Sampling' tutorial with my own protein. Within this tutorial the ChainB is used as the immobile reference - so it is this that is posit

Re: [gmx-users] Reg: Adding multiple molecules

2010-09-27 Thread vinothkumar mohanakrishnan
Hi Justin Thank you for your suggestion. Regards Vinoth On Mon, Sep 27, 2010 at 4:52 PM, Justin A. Lemkul wrote: > > > vinothkumar mohanakrishnan wrote: > >> Hi justin >> >> if genbox is not able to find sufficient space then why should one >> calculate number density for the given box dimensi

Re: [gmx-users] Pulling Simulation Query

2010-09-27 Thread Justin A. Lemkul
Natalie Stephenson wrote: Hi ... Just wondering if there's a way, during pulling simulations, to immobilise say the C-terminus of the protein. I've been running through the 'Umbrella Sampling' tutorial with my own protein. Within this tutorial the ChainB is used as the immobile referenc

[gmx-users] How to make a lipid bilayer with specific dimensions?

2010-09-27 Thread chris . neale
I don't think that this is currently causing anybody any problems, but note that genbox is going to cut any lipids that cross out of the central unit box (either because genbox is unaware of PBC or because these lipids now clash across PBC). Therefore: genbox -cs popc128a.pdb -o popc_new.p

[gmx-users] Salt concentration

2010-09-27 Thread Anthony Cruz Balberdi
Dear Users: I am interested in perform a simulation of a protein in different salt concentrations. How I can calculate how much Na+ and Cl- ions to include in the simulation to achieve certain concentration? Could someone point me in the right direction? Is there a software to do this type of calcu

Re: [gmx-users] Salt concentration

2010-09-27 Thread Gonçalo C . Justino
Hi Anthony, You can use the '-conc' tag from genion during file preparation, as per http://manual.gromacs.org/current/online/genion.html. Gonçalo On 27 September 2010 15:32, Anthony Cruz Balberdi wrote: > Dear Users: > I am interested in perform a simulation of a protein in different salt > co

[gmx-users] LJ potential

2010-09-27 Thread nishap . patel
Hello, Is there a way to turn off the 6-term in LJ 6-12 potential? -Nisha P. -- gmx-users mailing listgmx-users@gromacs.org http://lists.gromacs.org/mailman/listinfo/gmx-users Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/Search before posting! Please don't po

Re: [gmx-users] LJ potential

2010-09-27 Thread David van der Spoel
On 2010-09-27 18.37, nishap.pa...@utoronto.ca wrote: Hello, Is there a way to turn off the 6-term in LJ 6-12 potential? -Nisha P. Set the parameters to zero? -- David van der Spoel, Ph.D., Professor of Biology Dept. of Cell & Molec. Biol., Uppsala University. Box 596, 75124 Uppsala, Sweden.

[gmx-users] (no subject)

2010-09-27 Thread Miroslawa aleksandrowa
Dear users of GROMACS, I have a problem with force field . Can you tell me what other ff except OPLS-AA I can use for Al3+ or other trivalent ions ? Thank you Mira -- gmx-users mailing listgmx-users@gromacs.org http://lists.gromacs.org/mailman/listinfo/gmx-users Please search the

Re: [gmx-users] LJ potential

2010-09-27 Thread nishap . patel
How can I do that using OPLSAA, because in the nb.itp file the values are for sigma/epsilon? Can I change C6=0 in the code file? Quoting David van der Spoel : On 2010-09-27 18.37, nishap.pa...@utoronto.ca wrote: Hello, Is there a way to turn off the 6-term in LJ 6-12 potential? -Nisha P.

[gmx-users] MPI and dual-core laptop

2010-09-27 Thread simon sham
Hi, I wanted to test the GROMACS MPI version in my dual-processors laptop. I have installed openmpi 1.4.2 version. However, when I tried to configure GROMACS 4.5.1 with --enable-mpi option, I got the following configuration problem: "checking whether the MPI cc command works... configure: error:

Re: [gmx-users] MPI and dual-core laptop

2010-09-27 Thread Gonçalo C . Justino
Hi, My 2 cents on your problem: I've been running gromacs on nearly everything, from one to six cores. It works. Besides mpi being in your /usr/local/bin, is it on the path of your system? You can do, in the terminal export PATH=$PATH:/usr/local/bin Or, for a more permanent solution, go to your

[gmx-users] g_hbond segmentation fault

2010-09-27 Thread Yao Yao
-- Found 1048 different hydrogen bonds in trajectory Found 2292 different atom-pairs within hydrogen bonding distance Merging hbonds with Acceptor and Donor swapped 2/10915Segmentation fault __ Hi guys, I got the msg like above when

Re: [gmx-users] Tables with forcefield

2010-09-27 Thread Sai Pooja
Thanks M! I am using the standard 6-12 tables available with the gromacs package. For -table and -tablep options to start with. I want to understand 1 thing. The forcefield files and the topology file specifies sigma and epsilon parameters. If I change the combination rule to 1 in the forcefield.

Re: [gmx-users] g_hbond segmentation fault

2010-09-27 Thread Erik Marklund
Yao Yao skrev 2010-09-27 21.22: -- Found 1048 different hydrogen bonds in trajectory Found 2292 different atom-pairs within hydrogen bonding distance Merging hbonds with Acceptor and Donor swapped 2/10915Segmentation fault __ Hi gu

[gmx-users] Hardware-specific crash with 4.5.1

2010-09-27 Thread Justin A. Lemkul
Hi All, I'm hoping I might get some tips in tracking down the source of an issue that appears to be hardware-specific, leading to crashes in my system. The failures are occurring on our supercomputer (Mac OSX 10.3, PowerPC). Running the same .tpr file on my laptop (Mac OSX 10.5.8, Intel Cor

RE: [gmx-users] Salt concentration

2010-09-27 Thread Dallas Warren
Calculator or spreadsheet will do the job fine. Catch ya, Dr. Dallas Warren Medicinal Chemistry and Drug Action Monash Institute of Pharmaceutical Sciences, Monash University 381 Royal Parade, Parkville VIC 3010 dallas.war...@monash.edu +61 3 9909 9304 - When the o

[gmx-users] Gromacs 4.0.7 failing in parallel mode

2010-09-27 Thread Jyoti Mahalik
Hi, I am simulating two protein units in water and it works fine in the single processor but when I use the parallel version of the GROMACS, it fails. My script: grompp_mpi_d -f md.mdp -c hba_b4md.gro -r hba_b4md.gro -p hba.top -o hba_md.tpr mpirun -np $NPROCS mdrun_mpi_d -np $NPROCS -ddo

Re: [gmx-users] Gromacs 4.0.7 failing in parallel mode

2010-09-27 Thread Mark Abraham
- Original Message - From: Jyoti Mahalik Date: Tuesday, September 28, 2010 9:19 Subject: [gmx-users] Gromacs 4.0.7 failing in parallel mode To: gmx-users@gromacs.org > Hi, > I am simulating two protein units in water and it works > fine in the single processor but when I use the paral

Re: [gmx-users] Hardware-specific crash with 4.5.1

2010-09-27 Thread Roland Schulz
Justin, I think the interaction kernel is not OK on your PowerPC machine. I assume that from: 1) The force seems to be zero (minimization output). 2) When you use the all-to-all kernel which is not available for the powerpc kernel, it automatically falls back to the C kernel and then it works. Wh

[gmx-users] Psuedo ONIOM gromacs MD run

2010-09-27 Thread TJ Mustard
  Hello all,   Is it possible to do a pseudo ONIOM MD run where you hold all atoms rigid except for a small sphere of some radius. I was wondering if I could do this with posres.itp files or a define argument. It would be beneficial if the atoms held rigid had no energi

Re: [gmx-users] Psuedo ONIOM gromacs MD run

2010-09-27 Thread Mark Abraham
Freeze groups do something along these lines. See manual. Mark - Original Message - From: TJ Mustard Date: Tuesday, September 28, 2010 10:46 Subject: [gmx-users] Psuedo ONIOM gromacs MD run To: "gmx-users@gromacs.org" -

Re: [gmx-users] Hardware-specific crash with 4.5.1

2010-09-27 Thread Justin A. Lemkul
Roland Schulz wrote: Justin, I think the interaction kernel is not OK on your PowerPC machine. I assume that from: 1) The force seems to be zero (minimization output). 2) When you use the all-to-all kernel which is not available for the powerpc kernel, it automatically falls back to the C k

Re: [gmx-users] Psuedo ONIOM gromacs MD run

2010-09-27 Thread TJ Mustard
Thanks Mark I will. I have found that the more I learn about gromacs the less I know about grimacs. On September 28, 2010 at 2:52 AM Mark Abraham wrote: Freeze groups do something along these lines. See manual. Mark

[gmx-users] Another error that keeps comming up.

2010-09-27 Thread TJ Mustard
Hey all,   Here is another error that I keep getting. I am trying to "speed up" my md runs with -heavyh and longer time steps. I don't get LINCS errors but I do get this...       Back Off! I just backed up prlog.log to ./#prlog.log.1# Getting Loaded..

Re: [gmx-users] Hardware-specific crash with 4.5.1

2010-09-27 Thread Mark Abraham
- Original Message - From: "Justin A. Lemkul" Date: Tuesday, September 28, 2010 11:11 Subject: Re: [gmx-users] Hardware-specific crash with 4.5.1 To: Gromacs Users' List > > > Roland Schulz wrote: > >Justin, > > > >I think the interaction kernel is not OK on your PowerPC > machine.

Re: [gmx-users] Hardware-specific crash with 4.5.1

2010-09-27 Thread Justin A. Lemkul
Mark Abraham wrote: - Original Message - From: "Justin A. Lemkul" Date: Tuesday, September 28, 2010 11:11 Subject: Re: [gmx-users] Hardware-specific crash with 4.5.1 To: Gromacs Users' List > > > Roland Schulz wrote: > >Justin, > > > >I think the interaction kernel is not OK

Re: [gmx-users] Hardware-specific crash with 4.5.1

2010-09-27 Thread Mark Abraham
- Original Message - From: "Justin A. Lemkul" Date: Tuesday, September 28, 2010 11:39 Subject: Re: [gmx-users] Hardware-specific crash with 4.5.1 To: Discussion list for GROMACS users > > > Mark Abraham wrote: > > > > > >- Original Message - > >From: "Justin A. Lemkul" > >Da

RE: [gmx-users] How to make a lipid bilayer with specific dimensions?

2010-09-27 Thread NG HUI WEN
Thanks for that guys! I will try them out. Just a quick question here with regards to my lipids being possibly too small, does it have something to do with the minimum image criterion for PBC? I have ensured that the distance of edge of the box and the atoms at the side of the protein to be great

Re: [gmx-users] How to make a lipid bilayer with specific dimensions?

2010-09-27 Thread Justin A. Lemkul
NG HUI WEN wrote: Thanks for that guys! I will try them out. Just a quick question here with regards to my lipids being possibly too small, does it have something to do with the minimum image criterion for PBC? I have ensured that the distance of edge of the box and the atoms at the side of t

RE: [gmx-users] How to make a lipid bilayer with specific dimensions?

2010-09-27 Thread NG HUI WEN
Points noted, thanks Justin! -Original Message- From: gmx-users-boun...@gromacs.org [mailto:gmx-users-boun...@gromacs.org] On Behalf Of Justin A. Lemkul Sent: Tuesday, September 28, 2010 10:06 AM To: Discussion list for GROMACS users Subject: Re: [gmx-users] How to make a lipid bilayer wit

[gmx-users] Compilation error - GROMACS 4.5.1

2010-09-27 Thread Shachi Katira
Hi all, I get the following errors when I try to compile the double precision version of GROMACS 4.5.1 using: ./configure --prefix=/home/shachi/gromacs --enable-double --with-fft=fft3 --without-x Errors: In file included from nb_kernel400_x86_64_sse2.c:24: ../../../../include/gmx_sse2_double.h:

Re: [gmx-users] Compilation error - GROMACS 4.5.1

2010-09-27 Thread Mark Abraham
- Original Message - From: Shachi Katira Date: Tuesday, September 28, 2010 12:56 Subject: [gmx-users] Compilation error - GROMACS 4.5.1 To: gmx-users@gromacs.org > Hi all, > > I get the following errors when I try to compile the double precision > version of GROMACS 4.5.1 using: > >

Re: [gmx-users] Hardware-specific crash with 4.5.1

2010-09-27 Thread Roland Schulz
On Mon, Sep 27, 2010 at 9:58 PM, Mark Abraham wrote: > > > - Original Message - > From: "Justin A. Lemkul" > Date: Tuesday, September 28, 2010 11:39 > Subject: Re: [gmx-users] Hardware-specific crash with 4.5.1 > To: Discussion list for GROMACS users > > > > > > > Mark Abraham wrote: > >