Re: [gmx-users] Freezing a portion of a protein during simulation

2010-04-27 Thread Anirban Ghosh
Hello Justin, In my topology file I am declaring: --- ; Include Position restraint file #ifdef POSRES #include "posre.itp" #endif ; Strong position restraints on rest of B2AR #ifdef STRONG_

Re: [gmx-users] more than one peptide in one simulation box

2010-04-27 Thread Mark Abraham
On 28/04/10 15:13, shahid nayeem wrote: Hi Mark How one should be certain that this much trajectory is long enough to get coverged ensemble. When the observables of interest aren't changing... This is a "how long is a piece of string?"-type question. Read some literature about simulations of

Re: [gmx-users] more than one peptide in one simulation box

2010-04-27 Thread shahid nayeem
Hi Mark How one should be certain that this much trajectory is long enough to get coverged ensemble. Shahid On 4/27/10, Mark Abraham wrote: > > On 27/04/2010 8:58 PM, Justin A. Lemkul wrote: > >> >> >> shahid nayeem wrote: >> >>> Dear Mark >>> Following your advice I started using three peptid

Re: [gmx-users] more than one peptide in one simulation box

2010-04-27 Thread shahid nayeem
Hi Justin Should I try to do position restraint at 500k and then full MD simulation. shahid On 4/27/10, Justin A. Lemkul wrote: > > > > shahid nayeem wrote: > >> My peptide is 26 residue alpha helix obtained from crystal structure .pdb >> file. I am posting energy minimization, position restari

Re: [gmx-users] how many dihedral angles should be included for the naphthalene molecule if I use oplsaa

2010-04-27 Thread Justin A. Lemkul
Ming Han wrote: 1__2 6__ // 3\/ \\ |||4| \\___ /\___// 5 7 I want to know if 1-2-3-6 torsion should be included? And if both 2-3-4-7 and 6-3-4-5 should be included? Thanks. I would think that proper dihedrals would not even be used for such a molecule. The fuse

Re: [gmx-users] how many dihedral angles should be included for the naphthalene molecule if I use oplsaa

2010-04-27 Thread Ming Han
1__2 6__ // 3\/ \\ |||4| \\___ /\___// 5 7 I want to know if 1-2-3-6 torsion should be included? And if both 2-3-4-7 and 6-3-4-5 should be included? Thanks. Ming On Tue, Apr 27, 2010 at 6:43 PM, Justin A. Lemkul wrote: > > > Ming Han wrote: >> >> Dear all, >> >> I wa

Re: [gmx-users] how many dihedral angles should be included for the naphthalene molecule if I use oplsaa

2010-04-27 Thread Justin A. Lemkul
Ming Han wrote: Dear all, I want to test a MD simulation job for the naphthalene molecule by the oplsaa force field by. The x2top generated less number of dihedral angles than what I expected. Do you know how many dihedral angles should be included for the naphthalene molecule if I use oplsaa?

[gmx-users] how many dihedral angles should be included for the naphthalene molecule if I use oplsaa

2010-04-27 Thread Ming Han
Dear all, I want to test a MD simulation job for the naphthalene molecule by the oplsaa force field by. The x2top generated less number of dihedral angles than what I expected. Do you know how many dihedral angles should be included for the naphthalene molecule if I use oplsaa? Thanks a lot. Min

[gmx-users] RDF

2010-04-27 Thread Nilesh Dhumal
Hello, I am trying to plot radial distribution function between a atom and a center of two atoms. How can I calculate the centre of two atoms and further how can I use this center to plot radial distribution funciton? THanks Nilesh -- gmx-users mailing listgmx-users@gromacs.org http://li

Re: [gmx-users] h_bond error

2010-04-27 Thread Maurício Menegatti Rigo
Thanks J. My mistake. -- gmx-users mailing listgmx-users@gromacs.org http://lists.gromacs.org/mailman/listinfo/gmx-users Please search the archive at http://www.gromacs.org/search before posting! Please don't post (un)subscribe requests to the list. Use the www interface or send it to gmx-us

Re: [gmx-users] h_bond error

2010-04-27 Thread Justin A. Lemkul
Maurício Menegatti Rigo wrote: Hi, I'm trying to compute H bonds in a protein-ligand system. I got the following error: Program g_hbond, VERSION 4.0.7 Source code file: gmx_hbond.c, line: 565 Fatal error: Error in func_type Position Rest

[gmx-users] h_bond error

2010-04-27 Thread Maurício Menegatti Rigo
Hi, I'm trying to compute H bonds in a protein-ligand system. I got the following error: Program g_hbond, VERSION 4.0.7 Source code file: gmx_hbond.c, line: 565 Fatal error: Error in func_type Position Rest. ---

Re: [gmx-users] I noticed implicit_solvent in the mdout.mdp. Is it going to be implemented soon?

2010-04-27 Thread Justin A. Lemkul
Arthur Roberts wrote: Hi, all, I noticed that there was an implicit_solvent term in the mdout.mdp. However, it didn't seem operative. I am just curious. Is it going to be implemented soon? Implicit solvent is implemented in the development code and (as far as I have heard) will be part

[gmx-users] I noticed implicit_solvent in the mdout.mdp. Is it going to be implemented soon?

2010-04-27 Thread Arthur Roberts
Hi, all, I noticed that there was an implicit_solvent term in the mdout.mdp. However, it didn't seem operative. I am just curious. Is it going to be implemented soon? Sincerely, Art Art Roberts 3950 Mahaila Ave G18 San Diego, CA 92122 cell: 206-850-7468 email: aroberts99...@yahoo.com sk

[gmx-users] Solvent addition

2010-04-27 Thread teklebrh
Dear Gromacs users, I have a box 10 X 10 X 10 with 20 solute molecules in side it. I want to add a solvent of 25% X and 75% Y and so on to see the aggregation behavior of my molecules in different ratio of solvents X and Y. Can any body help me how to do this. The two solvent is mixed. e.

[gmx-users] Re: Solvent addition

2010-04-27 Thread Justin A. Lemkul
tekle...@ualberta.ca wrote: Dear Gromacs users, I have a box 10 X 10 X 10 with 20 solute molecules in side it. I want to add a solvent of 25% X and 75% Y and so on to see the aggregation behavior of my molecules in different ratio of solvents X and Y. Can any body help me how to do this. Th

[gmx-users] Ensemble averaging between two proteins in the same box.

2010-04-27 Thread Arthur Roberts
Hi, all, I am trying to do distance restraints under two different conditions of the protein. I have both of these proteins in the same box. Because the experimentally determined distance restraints are from the same signal, there are distance restraints between similar atoms in both st

Re: [gmx-users] help

2010-04-27 Thread Justin A. Lemkul
Lucio Ricardo Montero Valenzuela wrote: hello. I got the same error, because I was using an NPT ensemble to do a position-restrained simulation, and the error disappeared when I did the position-restrained simulation with NVT. I would be careful about ad hoc changes to the statistical mechani

Re: [gmx-users] help

2010-04-27 Thread Lucio Ricardo Montero Valenzuela
hello. I got the same error, because I was using an NPT ensemble to do a position-restrained simulation, and the error disappeared when I did the position-restrained simulation with NVT. Best regards. Lucio. El mar, 27-04-2010 a las 16:52 +0530, vani ms escribió: > hello > i got the warning as the

Re: [gmx-users] Constrained simulations crash when bilayer moves in the z direction in the box

2010-04-27 Thread XAvier Periole
you should remove the water and lipid bilayer COM separately. I am not sure what you should do with your small molecule though. Probably best to add it to the bilayer as you constrain their relative position! On Apr 27, 2010, at 4:32 PM, ERIKSSON, EMMA wrote: Hi again, Thanks Xavier for your r

Re: [gmx-users] Constrained simulations crash when bilayer moves in the z direction in the box

2010-04-27 Thread Justin A. Lemkul
ERIKSSON, EMMA wrote: Hi again, Thanks Xavier for your reply. I didn't know that this mdp option existed. However, I read the manual and also checked the mdout.mdp files for my previous simulations, and I understood it as if those are the default settings even if you don't specify any of the

RE: [gmx-users] Constrained simulations crash when bilayer moves in the z direction in the box

2010-04-27 Thread ERIKSSON, EMMA
Hi again, Thanks Xavier for your reply. I didn't know that this mdp option existed. However, I read the manual and also checked the mdout.mdp files for my previous simulations, and I understood it as if those are the default settings even if you don't specify any of them in the md.mdp file. The

Re: [gmx-users] mutation to/from proline

2010-04-27 Thread Mark Abraham
On 25/04/2010 10:18 PM, afsaneh maleki wrote: Hi, I want to calculate relative free energy associated to mutation of P1 (native protein) to p2 (mutated protein).In this mutation, Isolusine (in P1) is mutated to Proline(in P2). With using Thermodynamic cycle: Bilayer+P1= = => Bilayer-P1 delta G1

Re: [gmx-users] more than one peptide in one simulation box

2010-04-27 Thread Mark Abraham
On 27/04/2010 8:58 PM, Justin A. Lemkul wrote: shahid nayeem wrote: Dear Mark Following your advice I started using three peptide in one simulation box. Iwas able to add these with genconf as previously in ordered manner, generated .gro with genconf, solvated it and after energy minimization I

Re: [gmx-users] more than one peptide in one simulation box

2010-04-27 Thread Justin A. Lemkul
shahid nayeem wrote: My peptide is 26 residue alpha helix obtained from crystal structure .pdb file. I am posting energy minimization, position restarint and full MD simulation .mdp file ref_t = 300 300 Here, you're equilibrating at 300 K... ref_t = 500 500 and here, you're ru

Re: [gmx-users] more than one peptide in one simulation box

2010-04-27 Thread shahid nayeem
My peptide is 26 residue alpha helix obtained from crystal structure .pdb file. I am posting energy minimization, position restarint and full MD simulation .mdp file Energy minimization cpp = /usr/bin/cpp define = -DFLEX_SPC constraints = none integrator = steep nsteps = 3000 ; ; Energy mi

Re: [gmx-users] Normal Mode Analysis

2010-04-27 Thread Ran Friedman
Hi, NMA is not MD - for one thing you don't run an NMA simulation for a certain time. I suggest you read about NMA and make sure you understand what the method does and what it can achieve before continuing. There is some data on the manual, a lot of data on the web and even more in books. When yo

Re: [gmx-users] Constrained simulations crash when bilayer moves in the z direction in the box

2010-04-27 Thread XAvier Periole
The crash seems to be expected as by crossing the pbc the distance will change significantly and in way the algorithm can not handle. Note that the overall translational motion of your system should always be removed. The removal of the COM motion of your bilayer should be sufficient to

Re: [gmx-users] help

2010-04-27 Thread Justin A. Lemkul
vani ms wrote: hello i got the warning as the Warning: "1-4 interaction between 246 and 251 at distance 4.991 which is larger than the 1-4 table size 2.400 nm These are ignored for the rest of the simulation This usually means your system is exploding, if not, you should increase table-exte

[gmx-users] help

2010-04-27 Thread vani ms
hello i got the warning as the Warning: "1-4 interaction between 246 and 251 at distance 4.991 which is larger than the 1-4 table size 2.400 nm These are ignored for the rest of the simulation This usually means your system is exploding, if not, you should increase table-extension in your mdp file

[gmx-users] Normal Mode Analysis

2010-04-27 Thread Anirban Ghosh
Hi ALL, This may sound like a very basic question, but I am still pondering over it. I have simulated a membrane protein system for 30 ns after Steepest Descent minimization and now I want to perform NMA. From the help pages what I understand is that I need a very well minimized system (using l-bf

Re: [gmx-users] more than one peptide in one simulation box

2010-04-27 Thread Justin A. Lemkul
shahid nayeem wrote: Dear Mark Following your advice I started using three peptide in one simulation box. Iwas able to add these with genconf as previously in ordered manner, generated .gro with genconf, solvated it and after energy minimization I did MD run for 10ns. Everything ran well. In

[gmx-users] Constrained simulations crash when bilayer moves in the z direction in the box

2010-04-27 Thread ERIKSSON, EMMA
Hi all! I have a problem regarding lipid bilayer simulations in Gromacs 4. During some of my simulations the whole system is moving in the z direction in the box, meaning that after some time the lipids are going out in the bottom of the box and coming in in the top of the box, since I'm using

Re: [gmx-users] more than one peptide in one simulation box

2010-04-27 Thread shahid nayeem
Dear Mark Following your advice I started using three peptide in one simulation box. Iwas able to add these with genconf as previously in ordered manner, generated .gro with genconf, solvated it and after energy minimization I did MD run for 10ns. Everything ran well. In the end when I see the traj