[gmx-users] TPI run

2012-10-10 Thread rikhiag
Dear users, I want to run TPI of water in a binary mixture of water DMSO system. My question is, how to add the extra water or extra DMSO molecule, or precisely, how to get the NEW set of coordinates for one water OR one DMSO molecule? I will appreciate your reply. Rikg -- View this message

[gmx-users] Secondary structure after InflateGRO step - protein in lipid membrane

2012-10-10 Thread Shima Arasteh
Hi, I'm doing the protein insertion into the lipid membrane. After shrinking for 30 times, I visualized the system, but the secondary structure of protein is still incomplete . How come? I finished the InflateGRO step. Thanks for your suggestions in advance. Sincerely, Shima -- gmx-users

Re: [gmx-users] rapid change in the angle of simulation cell - regd

2012-10-10 Thread ramesh cheerla
Dear Tsjerk, I am very thankful to you for your reply and sorry for delay as it has taken much time to understand the things properly. I have gone through the material that you have attached with this mail that is really useful. Now I have understood that large change in

Re: [gmx-users] Could not converge NPT constraints

2012-10-10 Thread juan-manuel.castillo
Thank you, At the end the problem was solved just increasing the time coupling constant of the barostat tau_p. The original value of 0.15 that I set at the beginning was too small, with 0.5 works fine. On Tue, 09 Oct 2012 14:05:43 -0400 Justin Lemkul jalem...@vt.edu wrote: On 10/9/12

Re: [gmx-users] rapid change in the angle of simulation cell - regd

2012-10-10 Thread Tsjerk Wassenaar
Hi Ramesh, You already got a good handle on it! You can subtract and add vectors without changing the lattice. Now, the box in gromacs will be something like ((ax, 0, 0), (bx,by,0),(cx,cy,cz)). If cx 0.5*ax, then you can change the vector to (cx-ax,cy-0,cz-0), and the new |cx'| will always be

[gmx-users] Re: Could not converge NPT constraints

2012-10-10 Thread escajarro
Thank you, At the end the problem was solved just increasing the time coupling constant of the barostat tau_p. The original value of 0.15 that I set at the beginning was too small, with 0.5 works fine. -- View this message in context:

Re: [gmx-users] Martini FF for Sec structure changes...

2012-10-10 Thread XAvier Periole
Martini FF cannot model changes in secondary structure ... other CG FF can. You'll find them easily in the literature. Notably the ones from Deserno or Derreumaux. On Oct 10, 2012, at 2:03 PM, rama david wrote: Hi friends, I planed to use the martini force-field for my simulation study

Re: [gmx-users] Martini FF for Sec structure changes...

2012-10-10 Thread rama david
Thank you for your reply, Are these Cg can be used in Gromacs. Thank you in advance. With best wishes and regards, Rama david On Wed, Oct 10, 2012 at 6:13 PM, XAvier Periole x.peri...@rug.nl wrote: Martini FF cannot model changes in secondary structure ... other CG FF can. You'll find

Re: [gmx-users] Martini FF for Sec structure changes...

2012-10-10 Thread XAvier Periole
Nope, but on other softwares. On Oct 10, 2012, at 2:50 PM, rama david wrote: Thank you for your reply, Are these Cg can be used in Gromacs. Thank you in advance. With best wishes and regards, Rama david On Wed, Oct 10, 2012 at 6:13 PM, XAvier Periole x.peri...@rug.nl wrote: Martini

Re: [gmx-users] Martini FF for Sec structure changes...

2012-10-10 Thread rama david
Hi thank you Please told me the name of Freely available software on which these FF can be used .. Thank you in advance With best wishes and regards, Rama david On Wed, Oct 10, 2012 at 6:26 PM, XAvier Periole x.peri...@rug.nl wrote: Nope, but on other softwares. On Oct 10, 2012, at

Re: [gmx-users] Martini FF for Sec structure changes...

2012-10-10 Thread francesco oteri
Hi rama, actually MARTINI has been further improved to allow secondary structure change. The title is: *Improving Internal Peptide Dynamics in the Coarse-Grained MARTINI Model: Toward Large-Scale Simulations of Amyloid- and Elastin-like Peptides* and here there is a link to the paper:

Re: [gmx-users] Martini FF for Sec structure changes...

2012-10-10 Thread XAvier Periole
No, it does NOT! HAve you red the paper?! This implementation is an adoc representation to mimic specific sequences of short peptides! They do never for any secondary structure! On Oct 10, 2012, at 3:05 PM, francesco oteri wrote: Hi rama, actually MARTINI has been further improved to

Re: [gmx-users] Secondary structure after InflateGRO step - protein in lipid membrane

2012-10-10 Thread Justin Lemkul
On 10/10/12 5:43 AM, Shima Arasteh wrote: Hi, I'm doing the protein insertion into the lipid membrane. After shrinking for 30 times, I visualized the system, but the secondary structure of protein is still incomplete . How come? I finished the InflateGRO step. Visualization oddities

[gmx-users] Secondary structure after InflateGRO step - protein in lipid membrane

2012-10-10 Thread Christopher Neale
What, precisely, do you mean the SS of the protein is still incomplete? When you use inflategro, the conformation of your protein is intended to remain as in your initial conformation. Did you forget to use position restraints on your protein during the contraction steps (

[gmx-users] pull=constraint gives zero forces

2012-10-10 Thread Christopher Neale
How can there be forces for holonomic constraints? Is this described by an equation in the manual? Just because there are values in the pullf.xvg file does not mean that these values are forces. If they are forces, what is the force constant and what is the equation that defines this force?

Re: [gmx-users] Secondary structure after InflateGRO step - protein in lipid membrane

2012-10-10 Thread Shima Arasteh
I exactly followed the steps in Justin's tutorial of KALP15-DPPC, so I applied position restrain to restrain the protein in the step before iteration. It doesn't seem to be artifact, because when I select only the protein in VMD, I see a broken secondary structure. Sometime ago, other gmx

[gmx-users] Secondary structure after InflateGRO step - protein in lipid membrane

2012-10-10 Thread Christopher Neale
Can you compute the Ca RMSD to the starting structure (using structural fitting) as a function of inflategro step and post it somewhere online and link it here? It is very difficult to know what is going on based simply on the fact that VMD doesn't render your protein as you expect it to.

[gmx-users] genbox Not enough memory

2012-10-10 Thread Christopher Neale
Sounds like you ran out of memory. Many clusters have a few large-memory nodes. Can you use one of those? It's failing on a call for 1.3 Gb of memory, which by itself isn't really a lot... Also, can you confirm 250 A box length, not 250 nm box length? Gromacs defines length in units of nm.

Re: [gmx-users] genbox Not enough memory

2012-10-10 Thread Juliette N.
Thanks..You are right...The last line of gro file says 250 so it is in nm!... On 10 October 2012 12:30, Christopher Neale chris.ne...@mail.utoronto.ca wrote: Sounds like you ran out of memory. Many clusters have a few large-memory nodes. Can you use one of those? It's failing on a call for

[gmx-users] genbox Not enough memory

2012-10-10 Thread Christopher Neale
So perhaps you don't need a box this big? a box length of 25 nm might be within your available memory. Otherwise, talk to your cluster's sysadmins. I am not sure that you're going to get any useful statistics on a box of length 250 nm in any event. A box that size will fit about 1 billion

[gmx-users] Re: Segmentation fault, mdrun_mpi

2012-10-10 Thread Ladasky
Update: Ladasky wrote Justin Lemkul wrote Random segmentation faults are really hard to debug. Can you resume the run using a checkpoint file? That would suggest maybe an MPI problem or something else external to Gromacs. Without a reproducible system and a debugging backtrace,

Re: [gmx-users] Re: Segmentation fault, mdrun_mpi

2012-10-10 Thread Justin Lemkul
On 10/10/12 1:33 PM, Ladasky wrote: Update: Ladasky wrote Justin Lemkul wrote Random segmentation faults are really hard to debug. Can you resume the run using a checkpoint file? That would suggest maybe an MPI problem or something else external to Gromacs. Without a reproducible

Re: [gmx-users] Re: genbox Not enough memory

2012-10-10 Thread Juliette N.
Hi again, The reason I have this big box is that I have fully extended chains of the length of ~ 250 nm,. In fact this 250 nm is the minimum size that I can fit the chain in the box; and I am going to fill this box with solvent and use NPT to increase the density. So I dont need to fill up the

Re: [gmx-users] Re: genbox Not enough memory

2012-10-10 Thread Justin Lemkul
On 10/10/12 4:24 PM, Juliette N. wrote: Hi again, The reason I have this big box is that I have fully extended chains of the length of ~ 250 nm,. In fact this 250 nm is the minimum size that I can fit the chain in the box; and I am going to fill this box with solvent and use NPT to increase

[gmx-users] Re: genbox Not enough memory

2012-10-10 Thread Dr. Vitaly Chaban
On Wed, Oct 10, 2012 at 10:24 PM, Juliette N. joojoojo...@gmail.com wrote: Hi again, The reason I have this big box is that I have fully extended chains of the length of ~ 250 nm,. In fact this 250 nm is the minimum size that I can fit the chain in the box; and I am going to fill this box

[gmx-users] Re: genbox Not enough memory

2012-10-10 Thread Dr. Vitaly Chaban
By the way, if you place your solvent molecule (-cs .gro) into a very large box, such as 10x10x10nm, genbox may perform well. I did NOT try, just an idea. For the overall system setup, I would rather follow the previous comments anyway. On Wed, Oct 10, 2012 at 10:49 PM, Dr. Vitaly Chaban

[gmx-users] (no subject)

2012-10-10 Thread Ananya Mondal
Hi friends, I'm interested in using Thole type model for water (polarizable) TTM2 J. Phys. Chem. A, 2006, 110 (11), pp 4100 or TTM3 ;J. Chem. Phys. 128, 074506 (2008) ) model. Can someone send me itp file for water.. Thanks Ananya -- gmx-users mailing listgmx-users@gromacs.org

[gmx-users] Re: genbox Not enough memory

2012-10-10 Thread Juliette N.
Hi, I will try to compress before doing genbox. As for the number of atoms, as I said I *dont *intend to fill up the whole box with solvent. So I put around 4000 solvent molecules to get the desired polymer Wt%, hence there would be a lot of free space for a total of 120 000 atoms approx..Also I

[gmx-users] .n2t file for the CHARMM forcefield

2012-10-10 Thread Elie M
Dear all, Is the form of the .n2t file (read by x2top to produce the top file) the same for all forcefields? I am using the CHARMM forcefield but I could not find an ffcharmm27.n2t file. Regards Elie -- gmx-users mailing listgmx-users@gromacs.org

Re: [gmx-users] .n2t file for the CHARMM forcefield

2012-10-10 Thread Justin Lemkul
On 10/10/12 10:35 PM, Elie M wrote: Dear all, Is the form of the .n2t file (read by x2top to produce the top file) the same for all forcefields? I am using the CHARMM forcefield but I could not find an ffcharmm27.n2t file. The .n2t format is universal. Most force fields do not have .n2t

RE: [gmx-users] .n2t file for the CHARMM forcefield

2012-10-10 Thread Elie M
Ok. This means that I have to form the .n2t file for the CHARMM field as it is required by x2top. For the OPLSAA field the format was: Copls_1570.000 12.011 4H 0.108 H 0.108 H 0.108 C 0.150C opls_1580.000 12.0114H 0.108 H 0.108 C 0.150 C

Re: [gmx-users] .n2t file for the CHARMM forcefield

2012-10-10 Thread Justin Lemkul
On 10/10/12 11:11 PM, Elie M wrote: Ok. This means that I have to form the .n2t file for the CHARMM field as it is required by x2top. For the OPLSAA field the format was: Copls_1570.00012.011 4H 0.108 H 0.108 H 0.108 C 0.150C opls_1580.000 12.0114H

RE: [gmx-users] .n2t file for the CHARMM forcefield

2012-10-10 Thread Elie M
Understandable. I meant the content of the .n2t file itself ; will it contain the same info as the OPLSAA with just changes in values according to the CHARMM field? Thank you Elie Date: Wed, 10 Oct 2012 23:13:34 -0400 From: jalem...@vt.edu To: gmx-users@gromacs.org Subject: Re: [gmx-users]