Re: [gmx-users] Position restraints

2011-06-16 Thread Mark Abraham

On 16/06/2011 11:54 AM, Tom Dupree wrote:

Greetings all,

I have run into a little bit of a problem.

I am trying to simulate a hetro-dimer. Through previous work we have identified 
302  C-alphas (of 960 odd residues) that don't move much. Previously we 
position restrained these atoms in our simulations using charmm which gave us a 
significant speed up in simulation time. My supervisor is keen for me to do a 
comparison of that method with a completely unrestrained system using the 53A6 
FF. I managed to generate a position restraint file using the .gro file and 
genrestr.
However when I get to the simulation step I can't get grompp to work.


As genrestr -h will warn you, its output is only useful for the first 
molecule.


GROMACS position restraints will not give you any implementation-related 
speed-up, only one from reducing the sampling volume. Freeze groups will 
give you some speed-up (at the cost of making your results depend on 
your frozen configuration), but you will usually not be able to use both 
constraints and NPT as well.



I ran

genrestr -f em.gro -n subset.ndx -o subset.itp -disre -constr


Such a constraint or distance restraint matrix can only work 
intra-[moleculetype], and so is futile for your case.



As I understand my situation grompp is looking for atom numbers in subset.itp 
and matching them to the atom numbers in topol_Protein_chain_x.itp and since 
the atom numbers in subset.itp are based on the em.gro file it fails/goes out 
of range.


Yep


Thanks to some chain breaks in my structure I have 6 chains instead of 2, is 
there a way convert the .gro numbers into the .itp numbers? and hence generate 
my restraint files based on the topol_Protein_chain_x.itp files?


Sure, use editconf with some appropriate index groups (and maybe -resnr) 
to create subset structures that will then produce output from genrestr 
suitable for use with grompp.



   Furthermore can I actually do a constraint matrix over multiple molecules?


Not without uniting the [moleculetypes]. Do note that the various kinds 
of restraints, constraints and freeze groups are all distinct 
approaches, and attempting to mix them can be asking for trouble. Be 
sure to use the right terminology when asking for help :-)


Mark


Or is there a better way of achieving my desired result? (302 atoms across 2 
(6) chains constrained in their relative positions)


All the best,

Tom


Fatal error:
[ file subset.itp, line 145 ]:
Atom index (5222) in constraints out of bounds (1-5149).
This probably means that you have inserted topology section constraints
in a part belonging to a different molecule than you intended to.
In that case move the constraints section to the right molecule.
For more information and tips for troubleshooting, please check the GROMACS
website at http://www.gromacs.org/Documentation/Errors
---

My topology file
;
;   File 'topol.top' was generated
;   By user: onbekend (0)
;   On host: onbekend
;   At date: Tue May 31 11:44:23 2011
;
;   This is a standalone topology file
;
;   It was generated using program:
;   pdb2gmx - VERSION 4.5.3
;
;   Command line was:
;   pdb2gmx -f 3LP2.pdb -o 3lp2 -ignh -v
;
;   Force field was read from the standard Gromacs share directory.
;

; Include forcefield parameters
#include gromos53a6.ff/forcefield.itp

; Include chain topologies
#include topol_Protein_chain_A.itp
#ifdef POSRES
#include posre_Protein_chain_A.itp
#endif
#include topol_Protein_chain_A2.itp
#ifdef POSRES
#include posre_Protein_chain_A2.itp
#endif
#ifdef REST
#include subset.itp
#endif
#include topol_Protein_chain_B.itp
#ifdef POSRES
#include posre_Protein_chain_B.itp
#endif
#ifdef REST
#include subset.itp
#endif
#include topol_Protein_chain_B2.itp
#ifdef POSRES
#include posre_Protein_chain_B2.itp
#endif
#ifdef REST
#include subset.itp
#endif
#include topol_Protein_chain_B3.itp
#ifdef POSRES
#include posre_Protein_chain_B3.itp
#endif
#ifdef REST
#include subset.itp
#endif
#include topol_Protein_chain_B4.itp
#ifdef POSRES
#include posre_Protein_chain_B4.itp
#endif
#ifdef REST
#include subset.itp
#endif
___




--
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 post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.

Can't post? Read http://www.gromacs.org/Support/Mailing_Lists


[gmx-users] RE: Including quadrupole-charge interaction in GROMACS (WU Yanbin)

2011-06-16 Thread Ran Friedman
Hi,
Not implemented in Gromacs AFAIK.
Note that (1) even the definition of a quadruple isn't straightforwards and 
depends on the orientation and (2) common FF parameters were not developed with 
higher ordered multiple taken into account.

@Article{Plattner2009,
author = Plattner, N and Meuwly, M,
title = {Higher order multipole moments for molecular dynamics simulations},
journal = J Mol Model,
year = 2009,
volume = 15,
pages = 687-694
}

Ran

Ran Friedman
Biträdande Lektor (Assistant Professor)

Linnaeus University
School of Natural Sciences
391 82 Kalmar, Sweden

Norrgård, room 328d
+46 480 446 290 Telephone
+46 76 207 8763 Mobile
ran.fried...@lnu.se
http://lnu.se/ccbg



From: gmx-users-boun...@gromacs.org [gmx-users-boun...@gromacs.org] On Behalf 
Of gmx-users-requ...@gromacs.org [gmx-users-requ...@gromacs.org]
Sent: 16 June 2011 09:00
To: gmx-users@gromacs.org
Subject: gmx-users Digest, Vol 86, Issue 96

Send gmx-users mailing list submissions to
gmx-users@gromacs.org

To subscribe or unsubscribe via the World Wide Web, visit
http://lists.gromacs.org/mailman/listinfo/gmx-users
or, via email, send a message with subject or body 'help' to
gmx-users-requ...@gromacs.org

You can reach the person managing the list at
gmx-users-ow...@gromacs.org

When replying, please edit your Subject line so it is more specific
than Re: Contents of gmx-users digest...


Today's Topics:

   1. Re: CHARMM forcefield (Justin A. Lemkul)
   2. Re: local pressure calcuation for Gromacs-4.5 (Justin A. Lemkul)
   3. Including quadrupole-charge interaction in GROMACS (WU Yanbin)
   4. g_density (chris.ne...@utoronto.ca)
   5. Position restraints (Tom Dupree)
   6. Re: Position restraints (Mark Abraham)


--

Message: 1
Date: Wed, 15 Jun 2011 17:11:08 -0400
From: Justin A. Lemkul jalem...@vt.edu
Subject: Re: [gmx-users] CHARMM forcefield
To: Discussion list for GROMACS users gmx-users@gromacs.org
Message-ID: 4df91fec.8080...@vt.edu
Content-Type: text/plain; charset=ISO-8859-1; format=flowed



simon sham wrote:
 Hi,
 I have two questions about using the charmm force field.
 1. Do we still need to use the two perl scripts, convert_charmm_
 to_gromacs.pl and fix_top_for_charmm.pl?

No.

 2. I got the following note when I tried to do energy minim. with grompp:

 NOTE 1 [file topol.top]:
   The largest charge group contains 12 atoms.
   Since atoms only see each other when the centers of geometry of the charge
   groups they belong to are within the cut-off distance, too large charge
   groups can lead to serious cut-off artifacts.
   For efficiency and accuracy, charge group should consist of a few atoms.
   For all-atom force fields use: CH3, CH2, CH, NH2, NH, OH, CO2, CO, etc.

 Is this a problem?


Yes.  CHARMM topologies should not use charge groups.  You should invoke pdb2gmx
with the -nochargegrp option to create a proper topology.

-Justin

 Thanks for your insight!

 Simon


--


Justin A. Lemkul
Ph.D. Candidate
ICTAS Doctoral Scholar
MILES-IGERT Trainee
Department of Biochemistry
Virginia Tech
Blacksburg, VA
jalemkul[at]vt.edu | (540) 231-9080
http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin




--

Message: 2
Date: Wed, 15 Jun 2011 17:17:53 -0400
From: Justin A. Lemkul jalem...@vt.edu
Subject: Re: [gmx-users] local pressure calcuation for Gromacs-4.5
To: Discussion list for GROMACS users gmx-users@gromacs.org
Message-ID: 4df92181.2090...@vt.edu
Content-Type: text/plain; charset=ISO-8859-1; format=flowed



Amit Choubey wrote:
 Dear all,

 Could anyone direct me to the manual for local pressure calculation or a
 place where everything is mentioned about it ? I have been only able to
 collect bits and pieces from the mailing lists.


That's probably all there is to be found.  There's no manual linked from the git
web interface or the ftp site.

-Justin

 Thank you
 Amit

 On Wed, Jun 15, 2011 at 5:35 AM, Mark Abraham mark.abra...@anu.edu.au
 mailto:mark.abra...@anu.edu.au wrote:

 On 15/06/2011 9:09 PM, Jianguo Li wrote:
 Dear all,

 I have made a test calculation of local pressure using version 4.5
 for my membrane simulation using CHARMM FF. When rerun the
 simulation, mdrun gives the localpressure data. Howeve, instead of
 giving an anveraged data of the local pressure, mdrun gives a
 separate file for each frame, so I got many files:
 localpressure.dat0, localpressure.dat1, localpressure.dat2,
 localpressure.dat3 ..
 Then I need to calculate the pressure tensor for each frame and
 make average. but these localpressure.dat files are very big (each
 file is about 30 Mb), occupying large space of the hard disk. Can
 

[gmx-users] Implicit solvent calculatiom

2011-06-16 Thread Netaly Khazanov
Hi,
My system is protein +DNA.
I am trying to perform implicit solvent calculation using Amber 99.
In the begging I wanted to minimize the system and got this error message.


GB parameter(s) missing or negative for atom type 'OS'

GB parameter(s) missing or negative for atom type 'H2'

GB parameter(s) missing or negative for atom type 'N*'

GB parameter(s) missing or negative for atom type 'CM'

GB parameter(s) missing or negative for atom type 'P'

---
Program grompp, VERSION 4.5.4
Source code file: grompp.c, line: 1123

Fatal error:
Can't do GB electrostatics; the implicit_genborn_params section of the
forcefield is missing parameters for 5 atomtypes or they might be negative.
For more information and tips for troubleshooting, please check the GROMACS
website at http://www.gromacs.org/Documentation/Errors

All the atoms, in which parametres are missing  belong to DNA.
Should I add them manually to gbsa.itp? ( but i don't know the right
parameters?).

Thanks in advance.
Netaly Khazanov.
Post-doctoral student.
-- 
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 post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/Support/Mailing_Lists

[gmx-users] Position restraints with multiple identical molecules

2011-06-16 Thread Craig Kitchen
Dear All, 

My system consists of 32 identical molecules (with 91 atoms per molecule). 
I would like to apply a position restraint to the 91st atom of each 
molecules. Since all my molecules are identical I only have one 
[moleculetype] which contains 91 atoms. I have defined -DPOSRES in the .mdp 
file and have placed:


#ifdef POSRES
#include posre.itp
#endif

in (hopefully) the correct part of the topology file (immediately below the 
only [moleculetype]). There are no solvent molecules or ions.


Applying a position restraint to atom 91 works fine with posre.itp 
containing:


[ position_restraints ]
; ai funct fc (in x, y and z)
91 1 2000 2000 2000

This results in just one restraint to atom 91. But if I want to restrain 
the 91st atom in the second molecule (atom 182 in the .gro file), how do I 
do it?


[ position_restraints ]
; ai funct fc (in x, y and z)
91 1 2000 2000 2000
182 1 2000 2000 2000

Does not work since 182 is outside the atom numbers in the [moleculetype] 
(1 to 91). Similarly, having 32 lines in posre.itp reading 91 1 2000 2000 
2000 gives 32 position restraints to atom 91.


I expect I can get around this problem by duplicating the current topology 
32 times and creating a huge topol.top file with all 2912 atoms defined 
explicitly. However, it seems there should be a more elegant way to achieve 
this. Any help would be greatly appreciated!


I have included below a simplified example with all the other atoms 
removed. In this case, how do you restrain the position of atom 2,3,4 and 
5?


Craig

Craig Kitchen
Department of Chemistry
University of Cambridge
Lensfield Road
Cambridge
CB2 1EW
UK

topol.top:

[ moleculetype ]
; Atom
; Name nrexcl
UNK3

[ atoms ]
;   nr  type  resnr resid  atom  cgnr   charge mass
1   amber99_41 1  UNK OAI 1   -0.50762  15.9994

; Include Position restraint file
#ifdef POSRES
#include posre.itp
#endif

[ system ]
; Name
Multiple Atoms

[ molecules ]
; Compound#mols
UNK   5

Starting conf.gro:

Multiple Atoms
5
   1UNKOAI1   0.767   1.527   0.970
   2UNKOAI2   0.988   0.716   1.774
   3UNKOAI3   0.556   0.587   0.434
   4UNKOAI4   0.335   1.656   2.436
   5UNKOAI5   0.766   1.527   2.373
  2.64620   3.76260   2.80720

File used for restraint coordinates (grompp -r) :

Restraint coordinates
   5
   1UNKOAI1   0.000   1.500   1.000
   2UNKOAI2   1.100   0.621   2.000
   3UNKOAI3   0.399   0.621   0.349
   4UNKOAI4   0.221   1.562  -0.343
   5UNKOAI5   0.923   1.561   2.458
  2.64620   3.76260   2.80720

--
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 post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.

Can't post? Read http://www.gromacs.org/Support/Mailing_Lists


Re: [gmx-users] Position restraints with multiple identical molecules

2011-06-16 Thread Justin A. Lemkul



Craig Kitchen wrote:

Dear All,
My system consists of 32 identical molecules (with 91 atoms per 
molecule). I would like to apply a position restraint to the 91st atom 
of each molecules. Since all my molecules are identical I only have one 
[moleculetype] which contains 91 atoms. I have defined -DPOSRES in the 
.mdp file and have placed:


#ifdef POSRES
#include posre.itp
#endif

in (hopefully) the correct part of the topology file (immediately below 
the only [moleculetype]). There are no solvent molecules or ions.


Applying a position restraint to atom 91 works fine with posre.itp 
containing:


[ position_restraints ]
; ai funct fc (in x, y and z)
91 1 2000 2000 2000

This results in just one restraint to atom 91. But if I want to restrain 
the 91st atom in the second molecule (atom 182 in the .gro file), how do 
I do it?




You already are.  The [position_restraints] directive is applied to the atom(s) 
of the [moleculetype] that contains it, so if you only have one [moleculetype] 
then any instance of it will have the restraint applied to whatever atoms are 
indicated.


-Justin


[ position_restraints ]
; ai funct fc (in x, y and z)
91 1 2000 2000 2000
182 1 2000 2000 2000

Does not work since 182 is outside the atom numbers in the 
[moleculetype] (1 to 91). Similarly, having 32 lines in posre.itp 
reading 91 1 2000 2000 2000 gives 32 position restraints to atom 91.


I expect I can get around this problem by duplicating the current 
topology 32 times and creating a huge topol.top file with all 2912 atoms 
defined explicitly. However, it seems there should be a more elegant way 
to achieve this. Any help would be greatly appreciated!


I have included below a simplified example with all the other atoms 
removed. In this case, how do you restrain the position of atom 2,3,4 
and 5?


Craig

Craig Kitchen
Department of Chemistry
University of Cambridge
Lensfield Road
Cambridge
CB2 1EW
UK

topol.top:

[ moleculetype ]
; Atom
; Name nrexcl
UNK3

[ atoms ]
;   nr  type  resnr resid  atom  cgnr   charge mass
1   amber99_41 1  UNK OAI 1   -0.50762  15.9994

; Include Position restraint file
#ifdef POSRES
#include posre.itp
#endif

[ system ]
; Name
Multiple Atoms

[ molecules ]
; Compound#mols
UNK   5

Starting conf.gro:

Multiple Atoms
5
   1UNKOAI1   0.767   1.527   0.970
   2UNKOAI2   0.988   0.716   1.774
   3UNKOAI3   0.556   0.587   0.434
   4UNKOAI4   0.335   1.656   2.436
   5UNKOAI5   0.766   1.527   2.373
  2.64620   3.76260   2.80720

File used for restraint coordinates (grompp -r) :

Restraint coordinates
   5
   1UNKOAI1   0.000   1.500   1.000
   2UNKOAI2   1.100   0.621   2.000
   3UNKOAI3   0.399   0.621   0.349
   4UNKOAI4   0.221   1.562  -0.343
   5UNKOAI5   0.923   1.561   2.458
  2.64620   3.76260   2.80720



--


Justin A. Lemkul
Ph.D. Candidate
ICTAS Doctoral Scholar
MILES-IGERT Trainee
Department of Biochemistry
Virginia Tech
Blacksburg, VA
jalemkul[at]vt.edu | (540) 231-9080
http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin


--
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 post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.

Can't post? Read http://www.gromacs.org/Support/Mailing_Lists


[gmx-users] Re: gmx-users Digest, Vol 86, Issue 97

2011-06-16 Thread ITHAYARAJA
Dear Sir,

I found an error which is given below as i run grompp for position-restraint


Fatal error:
Group DRG not found in indexfile.
Maybe you have non-default goups in your mdp file, while not using the '-n'
option of grompp.
In that case use the '-n' option.

The error was attempted with -n but the following error also found,

File input/output error:
index.ndx

I can do nothing by this simple statement, please help me, Sir

Regards with,

Ithayaraja M..
-- 
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 post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/Support/Mailing_Lists

Re: [gmx-users] Group DRG not found in indexfile

2011-06-16 Thread Justin A. Lemkul


Please choose an informative subject line, especially when initiating a new 
thread.

ITHAYARAJA wrote:

Dear Sir,

I found an error which is given below as i run grompp for 
position-restraint


Fatal error:
Group DRG not found in indexfile.
Maybe you have non-default goups in your mdp file, while not using the 
'-n' option of grompp.

In that case use the '-n' option.

The error was attempted with -n but the following error also found,

File input/output error:
index.ndx

I can do nothing by this simple statement, please help me, Sir



You've specified a group somewhere in the .mdp that doesn't exist.  Groups can 
be assigned by the name given in the [moleculetype] directive or by a custom 
name given in an index file.  Apparently you've invoked grompp -n but you don't 
have an index file.  You may or may not need one.  If you've used DRG but no 
[moleculetype] entries are named as such, use a proper name.  If you need to use 
some custom subset of atoms, make an index group.


-Justin

--


Justin A. Lemkul
Ph.D. Candidate
ICTAS Doctoral Scholar
MILES-IGERT Trainee
Department of Biochemistry
Virginia Tech
Blacksburg, VA
jalemkul[at]vt.edu | (540) 231-9080
http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin


--
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 post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.

Can't post? Read http://www.gromacs.org/Support/Mailing_Lists


Re: [gmx-users] Position restraints with multiple identical molecules

2011-06-16 Thread Craig Kitchen
You already are. The [position_restraints] directive is applied to the 
atom(s) of the [moleculetype] that contains it, so if you only have one 
[moleculetype] then any instance of it will have the restraint applied 
to whatever atoms are indicated.


-Justin


Hi Justin, thanks for your reply. 

If this is the case then should I not expect to see all the position 
restraints in the .tpr file?


Using gmxdump -s md1.tpr to inspect the file I can only see

Position Rest.:
nr: 2
iatoms:
   0 type=1 (POSRES) 0

(this is on the simplified 5 atom case). Which implies there is one 
position restraint on atom 1 (index 0 I guess).


Craig

Craig Kitchen wrote: Dear All, My system consists of 32 identical molecules 
(with 91 atoms per molecule). I would like to apply a position restraint to 
the 91st atom of each molecules. Since all my molecules are identical I 
only have one [moleculetype] which contains 91 atoms. I have defined 
-DPOSRES in the .mdp file and have placed: #ifdef POSRES #include 
posre.itp #endif in (hopefully) the correct part of the topology file 
(immediately below the only [moleculetype]). There are no solvent molecules 
or ions. Applying a position restraint to atom 91 works fine with posre.itp 
containing: [ position_restraints ] ; ai funct fc (in x, y and z) 91 1 2000 
2000 2000 This results in just one restraint to atom 91. But if I want to 
restrain the 91st atom in the second molecule (atom 182 in the .gro file), 
how do I do it?


You already are. The [position_restraints] directive is applied to the 
atom(s) of the [moleculetype] that contains it, so if you only have one 
[moleculetype] then any instance of it will have the restraint applied to 
whatever atoms are indicated.


-Justin

[ position_restraints ] ; ai funct fc (in x, y and z) 91 1 2000 2000 2000 
182 1 2000 2000 2000 Does not work since 182 is outside the atom numbers in 
the [moleculetype] (1 to 91). Similarly, having 32 lines in posre.itp 
reading 91 1 2000 2000 2000 gives 32 position restraints to atom 91. I 
expect I can get around this problem by duplicating the current topology 32 
times and creating a huge topol.top file with all 2912 atoms defined 
explicitly. However, it seems there should be a more elegant way to achieve 
this. Any help would be greatly appreciated! I have included below a 
simplified example with all the other atoms removed. In this case, how do 
you restrain the position of atom 2,3,4 and 5? Craig Craig Kitchen 
Department of Chemistry University of Cambridge Lensfield Road Cambridge 
CB2 1EW UK topol.top: [ moleculetype ] ; Atom ; Name nrexcl UNK 3 [ atoms ] 
; nr type resnr resid atom cgnr charge mass

  1   amber99_41 1  UNK OAI 1   -0.50762  15.9994
; Include Position restraint file
#ifdef POSRES
#include posre.itp
#endif
[ system ]
; Name
Multiple Atoms
[ molecules ]
; Compound#mols
UNK   5
Starting conf.gro:
Multiple Atoms
5
 1UNKOAI1   0.767   1.527   0.970
 2UNKOAI2   0.988   0.716   1.774
 3UNKOAI3   0.556   0.587   0.434
 4UNKOAI4   0.335   1.656   2.436
 5UNKOAI5   0.766   1.527   2.373
2.64620   3.76260   2.80720
File used for restraint coordinates (grompp -r) :
Restraint coordinates
 5
 1UNKOAI1   0.000   1.500   1.000
 2UNKOAI2   1.100   0.621   2.000
 3UNKOAI3   0.399   0.621   0.349
 4UNKOAI4   0.221   1.562  -0.343
 5UNKOAI5   0.923   1.561   2.458
2.64620   3.76260   2.80720

--


Justin A. Lemkul
Ph.D. Candidate
ICTAS Doctoral Scholar
MILES-IGERT Trainee
Department of Biochemistry
Virginia Tech
Blacksburg, VA
jalemkul[at]vt.edu | (540) 231-9080
http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin


-- gmx-users mailing list gmx-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 post (un)subscribe requests to the list. Use the www 
interface or send it to gmx-users-requ...@gromacs.org. Can't post? Read 
http://www.gromacs.org/Support/Mailing_Lists


--
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 post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.

Can't post? Read http://www.gromacs.org/Support/Mailing_Lists


[gmx-users] an error for .tpr files

2011-06-16 Thread Maria Hamilton
Hi all

I want to run a program in a Quad Core PC and in cluster. the program was
run succesfuly in my PC but when I run it in cluster I have the following
error. Would you please help me?


Can not open file:
topol.tpr

the .tpr file filename.tpr that I want run it is in my folder.


Thanks alot

Regards

Elena
-- 
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 post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/Support/Mailing_Lists

[gmx-users] residue HISB not found 4.5.3 versus 4.0.7 OPLSAA

2011-06-16 Thread maria goranovic
Hi

I had a simulation run with 4.0.7, using opls-aa for a protein in water.

I now use part of the protein coordinates output from the simulation, and
use pdb2gmx as:

pdb2gmx -f conf.gro -ff oplsaa -ignh

pdb2gmx complains that it cannot find HISB. I checked that the topology
files in 4.5.3 do not contain HISB, while those in 4.0.7 do.

It does not suffice to create another copy of the HISE topology in the .rtp
file and call it HISB, because a number of hydrogen atoms also have been
renamed?

What can be done? I want to use 4.5.3 because I would prefer to retain
residue numbering

Maria


-- 
Maria G.
Technical University of Denmark
Copenhagen
-- 
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 post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/Support/Mailing_Lists

Re: [gmx-users] Position restraints with multiple identical molecules

2011-06-16 Thread Craig Kitchen
Hi Justin, 

Ok - the .tpr file only lists the atoms from the molecule type. I get it 
now. Thanks again!


Craig



On 16 Jun 2011, at 12:58, Justin A. Lemkul wrote:



Craig Kitchen wrote: Dear All, My system consists of 32 identical molecules 
(with 91 atoms per molecule). I would like to apply a position restraint to 
the 91st atom of each molecules. Since all my molecules are identical I 
only have one [moleculetype] which contains 91 atoms. I have defined 
-DPOSRES in the .mdp file and have placed: #ifdef POSRES #include 
posre.itp #endif in (hopefully) the correct part of the topology file 
(immediately below the only [moleculetype]). There are no solvent molecules 
or ions. Applying a position restraint to atom 91 works fine with posre.itp 
containing: [ position_restraints ] ; ai funct fc (in x, y and z) 91 1 2000 
2000 2000 This results in just one restraint to atom 91. But if I want to 
restrain the 91st atom in the second molecule (atom 182 in the .gro file), 
how do I do it?


You already are. The [position_restraints] directive is applied to the 
atom(s) of the [moleculetype] that contains it, so if you only have one 
[moleculetype] then any instance of it will have the restraint applied to 
whatever atoms are indicated.


-Justin

[ position_restraints ] ; ai funct fc (in x, y and z) 91 1 2000 2000 2000 
182 1 2000 2000 2000 Does not work since 182 is outside the atom numbers in 
the [moleculetype] (1 to 91). Similarly, having 32 lines in posre.itp 
reading 91 1 2000 2000 2000 gives 32 position restraints to atom 91. I 
expect I can get around this problem by duplicating the current topology 32 
times and creating a huge topol.top file with all 2912 atoms defined 
explicitly. However, it seems there should be a more elegant way to achieve 
this. Any help would be greatly appreciated! I have included below a 
simplified example with all the other atoms removed. In this case, how do 
you restrain the position of atom 2,3,4 and 5? Craig Craig Kitchen 
Department of Chemistry University of Cambridge Lensfield Road Cambridge 
CB2 1EW UK topol.top: [ moleculetype ] ; Atom ; Name nrexcl UNK 3 [ atoms ] 
; nr type resnr resid atom cgnr charge mass

  1   amber99_41 1  UNK OAI 1   -0.50762  15.9994
; Include Position restraint file
#ifdef POSRES
#include posre.itp
#endif
[ system ]
; Name
Multiple Atoms
[ molecules ]
; Compound#mols
UNK   5
Starting conf.gro:
Multiple Atoms
5
 1UNKOAI1   0.767   1.527   0.970
 2UNKOAI2   0.988   0.716   1.774
 3UNKOAI3   0.556   0.587   0.434
 4UNKOAI4   0.335   1.656   2.436
 5UNKOAI5   0.766   1.527   2.373
2.64620   3.76260   2.80720
File used for restraint coordinates (grompp -r) :
Restraint coordinates
 5
 1UNKOAI1   0.000   1.500   1.000
 2UNKOAI2   1.100   0.621   2.000
 3UNKOAI3   0.399   0.621   0.349
 4UNKOAI4   0.221   1.562  -0.343
 5UNKOAI5   0.923   1.561   2.458
2.64620   3.76260   2.80720

--


Justin A. Lemkul
Ph.D. Candidate
ICTAS Doctoral Scholar
MILES-IGERT Trainee
Department of Biochemistry
Virginia Tech
Blacksburg, VA
jalemkul[at]vt.edu | (540) 231-9080
http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin


-- gmx-users mailing list gmx-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 post (un)subscribe requests to the list. Use the www 
interface or send it to gmx-users-requ...@gromacs.org. Can't post? Read 
http://www.gromacs.org/Support/Mailing_Lists


--
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 post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.

Can't post? Read http://www.gromacs.org/Support/Mailing_Lists


Re: [gmx-users] residue HISB not found 4.5.3 versus 4.0.7 OPLSAA

2011-06-16 Thread Mark Abraham

On 16/06/2011 11:01 PM, maria goranovic wrote:

Hi

I had a simulation run with 4.0.7, using opls-aa for a protein in water.

I now use part of the protein coordinates output from the simulation, 
and use pdb2gmx as:


pdb2gmx -f conf.gro -ff oplsaa -ignh

pdb2gmx complains that it cannot find HISB. I checked that the 
topology files in 4.5.3 do not contain HISB, while those in 4.0.7 do.


It does not suffice to create another copy of the HISE topology in the 
.rtp file and call it HISB, because a number of hydrogen atoms also 
have been renamed?


What can be done? I want to use 4.5.3 because I would prefer to retain 
residue numbering


Rename all the HISB residues in your input coordinate file to whatever 
OPLS/AA wants.


Mark
--
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 post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.

Can't post? Read http://www.gromacs.org/Support/Mailing_Lists


Re: [gmx-users] an error for .tpr files

2011-06-16 Thread Mark Abraham

On 16/06/2011 10:35 PM, Maria Hamilton wrote:

Hi all
I want to run a program in a Quad Core PC and in cluster. the program 
was run succesfuly in my PC but when I run it in cluster I have the 
following error. Would you please help me?

Can not open file:
topol.tpr
the .tpr file filename.tpr that I want run it is in my folder.


You're doing it wrong, but you haven't told us how you're doing 
anything. This error means you haven't specified -s filename.tpr, one 
way or another.


Mark
--
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 post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.

Can't post? Read http://www.gromacs.org/Support/Mailing_Lists


Re: [gmx-users] residue HISB not found 4.5.3 versus 4.0.7 OPLSAA

2011-06-16 Thread maria goranovic
thought of that, but the trajectory also has HISB data .. so I guess I will
have to alter it frame by frame ?

On Thu, Jun 16, 2011 at 4:32 PM, Mark Abraham mark.abra...@anu.edu.auwrote:

 On 16/06/2011 11:01 PM, maria goranovic wrote:

 Hi

 I had a simulation run with 4.0.7, using opls-aa for a protein in water.

 I now use part of the protein coordinates output from the simulation, and
 use pdb2gmx as:

 pdb2gmx -f conf.gro -ff oplsaa -ignh

 pdb2gmx complains that it cannot find HISB. I checked that the topology
 files in 4.5.3 do not contain HISB, while those in 4.0.7 do.

 It does not suffice to create another copy of the HISE topology in the
 .rtp file and call it HISB, because a number of hydrogen atoms also have
 been renamed?

 What can be done? I want to use 4.5.3 because I would prefer to retain
 residue numbering


 Rename all the HISB residues in your input coordinate file to whatever
 OPLS/AA wants.

 Mark
 --
 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 post (un)subscribe requests to the list. Use the www interface
 or send it to gmx-users-requ...@gromacs.org.
 Can't post? Read http://www.gromacs.org/Support/Mailing_Lists




-- 
Maria G.
Technical University of Denmark
Copenhagen
-- 
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 post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/Support/Mailing_Lists

Re: [gmx-users] residue HISB not found 4.5.3 versus 4.0.7 OPLSAA

2011-06-16 Thread Mark Abraham

On 17/06/2011 12:33 AM, maria goranovic wrote:
thought of that, but the trajectory also has HISB data .. 


.xtc and .trr trajectories have no atom identification. That's why lots 
of the tools require a .tpr or coordinate file, together with a 
trajectory file in order to function. Using bash, look at gmxdump -f 
traj.xtc 21|less to learn what is there. You merely need to be 
confident the atom ordering hasn't changed from whatever generated the 
old trajectory.



so I guess I will have to alter it frame by frame ?


Even if your trajectory is some pseudo-format like .pdb or .gro, you can 
convert it back to one of the above easily.


Mark

On Thu, Jun 16, 2011 at 4:32 PM, Mark Abraham mark.abra...@anu.edu.au 
mailto:mark.abra...@anu.edu.au wrote:


On 16/06/2011 11:01 PM, maria goranovic wrote:

Hi

I had a simulation run with 4.0.7, using opls-aa for a protein
in water.

I now use part of the protein coordinates output from the
simulation, and use pdb2gmx as:

pdb2gmx -f conf.gro -ff oplsaa -ignh

pdb2gmx complains that it cannot find HISB. I checked that the
topology files in 4.5.3 do not contain HISB, while those in
4.0.7 do.

It does not suffice to create another copy of the HISE
topology in the .rtp file and call it HISB, because a number
of hydrogen atoms also have been renamed?

What can be done? I want to use 4.5.3 because I would prefer
to retain residue numbering


Rename all the HISB residues in your input coordinate file to
whatever OPLS/AA wants.

Mark
-- 
gmx-users mailing list gmx-users@gromacs.org

mailto:gmx-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 post (un)subscribe requests to the list. Use the www
interface or send it to gmx-users-requ...@gromacs.org
mailto:gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/Support/Mailing_Lists




--
Maria G.
Technical University of Denmark
Copenhagen


-- 
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 post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/Support/Mailing_Lists

Re: [gmx-users] Including quadrupole-charge interaction in GROMACS

2011-06-16 Thread WU Yanbin
Hi, Erik,

Thanks for the information.

Regarding using more than two partial charges to mimic quadrupole
explicitly, is there a one-to-one relation between the quadrupole moments
and position and magnitude of the partial charges? For example, if I would
like to include quadrupole moment for carbon in graphene, which has only
non-zero Q20 quadrupole component, how should I place partial charges to
mimic that quadrupole moment?

Any hint or direction for literature is appreciated.

Best,
Yanbin

On Thu, Jun 16, 2011 at 3:24 AM, Erik Marklund er...@xray.bmc.uu.se wrote:

 Quadrupole interactions are implicit whenever you have a molecule withe
 more than two partial charges. You don't have to include it explicitly.
 There are, however, multipole expansion techniques to speed up
 electrostatics calculations, but that's another story.

 Erik

 15 jun 2011 kl. 23.45 skrev WU Yanbin:

  Dear GMXers,
 
  Is there a way in GROMACS to include quadrupole-charge interaction?
  Or is there a standard to way to mimic quadrupole moment using partial
 charge?
 
  Thank you.
 
  Best,
  Yanbin
  --
  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 post (un)subscribe requests to the list. Use the
  www interface or send it to gmx-users-requ...@gromacs.org.
  Can't post? Read http://www.gromacs.org/Support/Mailing_Lists

 ---
 Erik Marklund, PhD student
 Dept. of Cell and Molecular Biology, Uppsala University.
 Husargatan 3, Box 596,75124 Uppsala, Sweden
 phone:+46 18 471 4537fax: +46 18 511 755
 er...@xray.bmc.uu.se
 http://www2.icm.uu.se/molbio/elflab/index.html

-- 
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 post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/Support/Mailing_Lists

[gmx-users] g_dist

2011-06-16 Thread Nilesh Dhumal
Hello,

I have a system with 128 emi (cations) and 128 Cl (anions).

I want to calculate how many CL atoms are in cutoff distance relative to
hydorgen aotm of cation.  I considered all CL atoms are distinguishable.

Basically I want to calcualte the distance between each CL atom and
correponing hydorgen and registered those CL atoms which are whthin
cutoff.

How can I do ?

I am using Gromacs 4.0.7 version.

Thanks

Nilesh








-- 
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 post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/Support/Mailing_Lists


Re: [gmx-users] Including quadrupole-charge interaction in GROMACS

2011-06-16 Thread Erik Marklund
Hi,

For neutral species, I think it's straightforward (unless the paper Ran 
Friedman mentioned says otherwise). Just use the definition of the quadrupole 
moment to figure out where to put the charges and what magnitude to assign to 
them. The electric field from a quadrupole declines like r^-3, however, so its 
importance for graphene interactions seem very small to me.

Erik


16 jun 2011 kl. 17.10 skrev WU Yanbin:

 Hi, Erik,
 
 Thanks for the information.
 
 Regarding using more than two partial charges to mimic quadrupole explicitly, 
 is there a one-to-one relation between the quadrupole moments and position 
 and magnitude of the partial charges? For example, if I would like to 
 include quadrupole moment for carbon in graphene, which has only non-zero Q20 
 quadrupole component, how should I place partial charges to mimic that 
 quadrupole moment?
 
 Any hint or direction for literature is appreciated.
 
 Best,
 Yanbin
 
 On Thu, Jun 16, 2011 at 3:24 AM, Erik Marklund er...@xray.bmc.uu.se wrote:
 Quadrupole interactions are implicit whenever you have a molecule withe more 
 than two partial charges. You don't have to include it explicitly. There are, 
 however, multipole expansion techniques to speed up electrostatics 
 calculations, but that's another story.
 
 Erik
 
 15 jun 2011 kl. 23.45 skrev WU Yanbin:
 
  Dear GMXers,
 
  Is there a way in GROMACS to include quadrupole-charge interaction?
  Or is there a standard to way to mimic quadrupole moment using partial 
  charge?
 
  Thank you.
 
  Best,
  Yanbin
  --
  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 post (un)subscribe requests to the list. Use the
  www interface or send it to gmx-users-requ...@gromacs.org.
  Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
 
 ---
 Erik Marklund, PhD student
 Dept. of Cell and Molecular Biology, Uppsala University.
 Husargatan 3, Box 596,75124 Uppsala, Sweden
 phone:+46 18 471 4537fax: +46 18 511 755
 er...@xray.bmc.uu.se
 http://www2.icm.uu.se/molbio/elflab/index.html
 -- 
 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 post (un)subscribe requests to the list. Use the 
 www interface or send it to gmx-users-requ...@gromacs.org.
 Can't post? Read http://www.gromacs.org/Support/Mailing_Lists

---
Erik Marklund, PhD student
Dept. of Cell and Molecular Biology, Uppsala University.
Husargatan 3, Box 596,75124 Uppsala, Sweden
phone:+46 18 471 4537fax: +46 18 511 755
er...@xray.bmc.uu.se
http://www2.icm.uu.se/molbio/elflab/index.html

-- 
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 post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/Support/Mailing_Lists

[gmx-users] g_hbond

2011-06-16 Thread simon sham
Hi,
I am little confused about setting up the two groups in g_hbond.
1. Are the two groups corresponding to the donor (N) and the acceptor (O) atoms?

I want to measure the lifetime, angle and distance of the triplet N-HO 
hydrogen bond. Is this the triplet that the manual talks about?

Thanks for your insight.

Simon
-- 
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 post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/Support/Mailing_Lists

[gmx-users] g_dist

2011-06-16 Thread chris . neale

Dear Nilesh:

You don't seem to have made any progress since Mark gave you some  
reading hints. Your best chance to get a useful response follows the  
general form:


1. I want to do A.
2. I tried method B, here are the exact commands that I used (copy and paste)
3. With that method, I obtain C
4. But I am confused about why C does not look like D.
5. What is the problem?

Basically, it helps to show that you are doing your own work too.

For example, try a command on a single frame and then look at that  
frame with VMD and count for yourself and verify the answer that way.


Chris.

-- original message --

I have a system with 128 emi (cations) and 128 Cl (anions).

I want to calculate how many CL atoms are in cutoff distance relative to
hydorgen aotm of cation.  I considered all CL atoms are distinguishable.

Basically I want to calcualte the distance between each CL atom and
correponing hydorgen and registered those CL atoms which are whthin
cutoff.

How can I do ?

I am using Gromacs 4.0.7 version.

Thanks

Nilesh










--
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 post (un)subscribe requests to the list. Use the
www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/Support/Mailing_Lists


[gmx-users] Can't unfold the protein

2011-06-16 Thread Hsin-Lin Chiang
Hi,

I have question about unfolding.
I use three different ways respectively but all failed.
The three way is, 
1. 600K in 20ns, then the system explode.
2. 400K in 10ns, the protein is very stable and the value is almost the same in 
radius of gyration.
3. heating up from 300K to 400K in 2ns and the other 2ns for 400K only, then 
the protein is still stable and the value is almost the same in radius of 
gyration.
Below is the mdp file of the heating.
What's wrong with my system?

regards,
Hsin-Lin
-
title#160;#160;#160;#160;#160;#160;#160; = ttt
cpp#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;
 =#160; /lib/cpp
constraints#160;#160;#160;#160;#160;#160;#160;#160; =#160; hbonds
;define#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;
 =#160; -DFLEX_SPC
integrator#160;#160;#160;#160;#160;#160;#160;#160;#160; =#160; md
emtol#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;
 =#160; 100.0
emstep#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;
 =#160; 0.005
dt#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;
 =#160; 0.002#160;#160;#160; ; ps !
nsteps#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;
 =#160; 200#160; ; total 4 ns
nstcomm#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;#160; 
=#160; 500
nstxout#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;#160; 
=#160; 500
nstvout#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;#160; 
=#160; 500
nstfout#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;#160; 
=#160; 500
nstlog#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;
 =#160; 500
nstenergy#160;#160;#160;#160;#160;#160;#160;#160;#160;#160; =#160; 
500
nstlist#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;#160; 
=#160; 5
ns_type#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;#160; 
=#160; grid
rlist#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;
 =#160; 1.
rcoulomb#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;#160; 
=#160; 1.
rvdw#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;
 =#160; 1.
coulombtype#160;#160;#160;#160;#160;#160;#160;#160; =#160; PME
fourierspacing#160;#160;#160;#160;#160; =#160; 0.12
pme_order#160;#160;#160;#160;#160;#160;#160;#160;#160;#160; =#160; 4
optimize_fft#160;#160;#160;#160;#160;#160;#160; =#160; yes
Tcoupl#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;
 =#160; v-rescale
tc-grps#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;#160; 
=#160; Protein Non-Protein
;tau_t#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;
 =#160; 0.1#160; 0.1
tau_t#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;
 =#160; 0.2 0.2
ref_t#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;
 =#160; 300 300
energygrps#160;#160;#160;#160;#160;#160;#160;#160;#160; =#160; 
A-chain B-chain SOL#160; NA+
Pcoupl#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;
 =#160; berendsen
Pcoupltype#160;#160;#160;#160;#160;#160;#160;#160;#160; =#160; 
isotropic
;tau_p#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;
 =#160; 0.1
tau_p#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;
 =#160; 0.25
compressibility#160;#160;#160;#160; =#160; 5.4e-5
ref_p#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;
 =#160; 1.0
gen_vel#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;#160; 
=#160; yes
gen_temp#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;#160; 
=#160; 300
gen_seed#160;#160;#160;#160;#160;#160;#160;#160;#160;#160;#160; 
=#160; 173529
; Annealing
annealing#160;#160;#160;#160;#160;#160;#160;#160;#160;#160; =#160; 
single single
annealing_npoints#160;#160; =#160; 21 21
annealing_time#160;#160;#160;#160;#160; =#160; 0 100 200 300 400 500 600 
700 800 900 1000 1100 1200 1300 1400 1500 1600 1700 1800 1900 2000 0 100 200 
300 400 500 600 700 800 900 1000 1100 1200 1300 1400 1500 1600 1700 1800 1900 
2000
annealing_temp#160;#160;#160;#160;#160; =#160; 300 305 310 315 320 325 
330 335 340 345 350 355 360 365 370 375 380 385 390 395 400 300 305 310 315 320 
325 330 335 340 345 350 355 360 365 370 375 380 385 390 395 400
-

 
-- 
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 post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/Support/Mailing_Lists

[gmx-users] Can't unfold the protein

2011-06-16 Thread chris . neale

Dear Hsin-Lin:

The first thing that comes to mind is using very high temperatures with NVT.
Obtaining the desired denatured state is a complex challenge, but you  
might try this paper:


Phys. Rev. Lett. 93, 238105 (2004)
Reversible Temperature and Pressure Denaturation of a Protein  
Fragment: A Replica Exchange Molecular Dynamics Simulation Study


Still, your post seems to be about simply being unable to unfold at  
all. I have unfolded a very stable beta-barrel by either simulating at  
3000 K or using the pull code to pull the termini apart. This was not  
a real study -- I only used this to make a reverse folding movie for  
non-scientists -- but I can verify that it is possible to unfold even  
very stable proteins by these methods.


Chris.

-- original message --

Hi,

I have question about unfolding.
I use three different ways respectively but all failed.
The three way is,
1. 600K in 20ns, then the system explode.
2. 400K in 10ns, the protein is very stable and the value is almost  
the same in radius of gyration.
3. heating up from 300K to 400K in 2ns and the other 2ns for 400K  
only, then the protein is still stable and the value is almost the  
same in radius of gyration.

Below is the mdp file of the heating.
What's wrong with my system?



--
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 post (un)subscribe requests to the list. Use the
www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/Support/Mailing_Lists


Re: [gmx-users] Implicit solvent calculatiom

2011-06-16 Thread Mark Abraham

On 16/06/2011 7:59 PM, Netaly Khazanov wrote:

Hi,
My system is protein +DNA.
I am trying to perform implicit solvent calculation using Amber 99.
In the begging I wanted to minimize the system and got this error message.


GB parameter(s) missing or negative for atom type 'OS'

GB parameter(s) missing or negative for atom type 'H2'

GB parameter(s) missing or negative for atom type 'N*'

GB parameter(s) missing or negative for atom type 'CM'

GB parameter(s) missing or negative for atom type 'P'

---
Program grompp, VERSION 4.5.4
Source code file: grompp.c, line: 1123

Fatal error:
Can't do GB electrostatics; the implicit_genborn_params section of the 
forcefield is missing parameters for 5 atomtypes or they might be 
negative.
For more information and tips for troubleshooting, please check the 
GROMACS

website at http://www.gromacs.org/Documentation/Errors

All the atoms, in which parametres are missing  belong to DNA.
Should I add them manually to gbsa.itp? ( but i don't know the right 
parameters?).


I expect you're going to have to check in the literature whether such 
parameters exist. Unfortunately, not every combination of everything can 
be available in a nice shrink-wrapped package :) Once you've found them, 
then you'll have to fill in the columns of the above-mentioned section 
accordingly.


Mark
--
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 post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.

Can't post? Read http://www.gromacs.org/Support/Mailing_Lists


[gmx-users] Can't unfold the protein

2011-06-16 Thread Hsin-Lin Chiang
Dear Chris,

Thank you for your reply. 
My protein is very stable.

I simulated it in 300ns in 300K before but there was almost no change. 
That's why I want to do denaturation now.

I want to start a new simulation on the protein which is unfolded. 
And I'll read the paper you told me. 
Thank you.

But I don't understand what you recommend me to do. 
You use 3000K on your protein and the protein unfolded. 
I use 600K on my protein but the system explode(box become very big and protein 
locate at corner).

I saw the second way on this web, 
http://manual.gromacs.org/online/protunf.html 
But it also useless to me.

And in my mdp I use thermostat and barostat, which means my system is NPT. 
Does anything wrong in my methods?

Sincerely yours,

Hsin-Lin 
 
-- 
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 post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/Support/Mailing_Lists

[gmx-users] Gromacs to Amber trajectory

2011-06-16 Thread Raja Pandian
Dear All,

I have created Amber input file (.prmtop and .rst7) using Xleap. And then I
have converted this input file into Gromacs input file format (*.gro and
*.top) with amb2gmx.pl script. I did simulation using Gromacs Package (Amber
Force Field). Now I would like to convert Gromacs trajectory file to Amber
trajectory file format. So that I can easily analyses this trajectory. I
tried to convert Gromacs trajectory file to Amber trajectory file using VMD.
This VMD converted trajectory is not working fine. Pls help me in this
regards

Eagerly waiting for your suggestion.



Regards

Raja
-- 
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 post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/Support/Mailing_Lists

[gmx-users] Can't unfold the protein

2011-06-16 Thread chris . neale

Dear Hsin-Lin:

I am no expert in this area. I am just saying that if you get  
densities that are way too low in NPT, then you might alleviate this  
problem with NVT.


In fact, I would personally do this in vacuum without pbc and use the  
sd integrator. Then you can sample extended conformations. I'd  
consider this a method to generate unfolded conformations that  
probably have no relation to the true temperature denatured state. But  
let's be honest, you're never going to come close to Boltzmann  
sampling of an unfolded state in 300 ns so why bother with the water?


Your link ( http://manual.gromacs.org/online/protunf.html ) is only  
concerned with analysis, so it is irrelevant.


If you are determined to use NPT then I suppose that you could use the  
Berendsen barostat (more stable in simulations but you get the wrong  
ensemble) with a very low compressibility. That's about all I can say,  
perhaps somebody else will comment.


Good luck,

Chris.

-- original message --

Dear Chris,

Thank you for your reply.
My protein is very stable.

I simulated it in 300ns in 300K before but there was almost no change.
That's why I want to do denaturation now.

I want to start a new simulation on the protein which is unfolded.
And I'll read the paper you told me.
Thank you.

But I don't understand what you recommend me to do.
You use 3000K on your protein and the protein unfolded.
I use 600K on my protein but the system explode(box become very big  
and protein locate at corner).


I saw the second way on this web,
http://manual.gromacs.org/online/protunf.html
But it also useless to me.

And in my mdp I use thermostat and barostat, which means my system is NPT.
Does anything wrong in my methods?

Sincerely yours,

Hsin-Lin



--
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 post (un)subscribe requests to the list. Use the
www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/Support/Mailing_Lists


[gmx-users] gmx4.5.4 genion problem: No line with moleculetype 'SOL' found

2011-06-16 Thread Ye Yang
Hi, everyone:
  I am a new user of Gromacs, and I am running through the tutorial.
When I am trying to run the ligand-receptor binding tutorial from

http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin/gmx-tutorials/complex/02_topology.html
 I met some trouble in adding ion.
 Each time when I use genion, it shows:

Will try to add 0 NA ions and 6 CL ions.
Select a continuous group of solvent molecules
Group 0 ( System) has 33046 elements
Group 1 (Protein) has  1693 elements
Group 2 (  Protein-H) has  1301 elements
Group 3 (C-alpha) has   163 elements
Group 4 (   Backbone) has   489 elements
Group 5 (  MainChain) has   653 elements
Group 6 (   MainChain+Cb) has   805 elements
Group 7 (MainChain+H) has   815 elements
Group 8 (  SideChain) has   878 elements
Group 9 (SideChain-H) has   648 elements
Group10 (Prot-Masses) has  1693 elements
Group11 (non-Protein) has 31353 elements
Group12 (  Other) has15 elements
Group13 (JZ4) has15 elements
Group14 (  Water) has 31338 elements
Group15 (SOL) has 31338 elements
Group16 (  non-Water) has  1708 elements
Select a group: 15
Selected 15: 'SOL'
Number of (3-atomic) solvent molecules: 10446

Processing topology

Back Off! I just backed up temp.top to ./#temp.top.3#

---
Program genion_d, VERSION 4.5.4
Source code file: gmx_genion.c, line: 285

Fatal error:
No line with moleculetype 'SOL' found the [ molecules ] section of file
'topol.top'

I checked my topology file and it looks fine:
; Include Position restraint file
#ifdef POSRES
#include posre.itp
#endif

; Include water topology
#include gromos43a1.ff/spc.itp

; Include ligand topoloty
#include drg.itp
#ifdef POSRES_WATER
; Position restraint for each water oxygen
[ position_restraints ]
;  i funct   fcxfcyfcz
   11   1000   1000   1000
#endif

; Include topology for ions
#include gromos43a1.ff/ions.itp

[ system ]
; Name

Protein
[ molecules ]
; Compound#mols
Protein_chain_A 1
JZ4 1
SOL 10446

One thing that might happen is in the genion source file,  I read through
it, and the problem either happens in loading the line to buf2, or in the
gmx_strcmp(buf2,SOL), since clearly the SOL_line is -1 aftermath.

The other thing I tried is to remove the ligand JZ4 part in both topology
and coordinates file, and in this case, it works perfectly for adding ion.
But in this way, I do not know how to insert my ligand into the system since
it might collide with solvent.

Can someone help me with this problem?

Thank you all very much.

Ye Yang
-- 
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 post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/Support/Mailing_Lists

Re: [gmx-users] Can't unfold the protein

2011-06-16 Thread felmer...@uchile.cl

Hey,

 

NPT is not the appropriate way to do this kind of simulations. I am not sure 
whether or not the water models available for classic MD simulations are able 
to reproduce the phase behavior. Indeed what you see when your system explodes 
and gets huge is the water evaporating. What is generally done is to allow the 
system to reach the density of liquid water at the temperature that you desire 
(see the Dagget's papers of that) and then use NVT for the rest of the 
simulation. 

Regarding the kinetics of the process it is up to the protein. There are some 
proteins that unfold really fast and others reallly really slow. For instance, 
on one of my systems i have to simulate 50 ns at 600 K to only see a 40% 
decrease in structure. 

As you may already imagine, the pressure gets quite high in these kind of 
simulations. Extreme pressure can also induce protein unfolding, but the 
depence here is again protein dependent since moderately high pressure tends to 
stabilize the structure. Of course what is moderate and what is high is case 
dependent.

Please consider also that biomolecular forcefields were not designed with these 
extreme temperatures in mind, so while doing high temperature induced unfolding 
you are taking the FF to the limit and you should be extremely careful of what 
kind of conclusion you obtain from your simulation.

 

Regards

 

Felipe
Mensaje originalDe: chris.neale@utoronto.caFecha: 16-jun-2011 
12:58Para: gmx-users@gromacs.orgAsunto: [gmx-users] Canamp;#39;t unfold the 
proteinDear Hsin-Lin:I am no expert in this area. I am just saying that if you 
get  densities that are way too low in NPT, then you might alleviate this  
problem with NVT.In fact, I would personally do this in vacuum without pbc and 
use the  sd integrator. Then you can sample extended conformations. I'd  
consider this a method to generate unfolded conformations that  probably have 
no relation to the true temperature denatured state. But  let's be honest, 
you're never going to come close to Boltzmann  sampling of an unfolded state in 
300 ns so why bother with the water?Your link ( 
http://manual.gromacs.org/online/protunf.html ) is only  concerned with 
analysis, so it is irrelevant.If you are determined to use NPT then I suppose 
that you could use the  Berendsen barostat (more stable in simulations but you 
get the wrong  ensemble) with a very low compressibility. That's about all I 
can say,  perhaps somebody else will comment.Good luck,Chris.-- original 
message --Dear Chris,Thank you for your reply.My protein is very stable.I 
simulated it in 300ns in 300K before but there was almost no change.That's why 
I want to do denaturation now.I want to start a new simulation on the protein 
which is unfolded.And I'll read the paper you told me.Thank you.But I don't 
understand what you recommend me to do.You use 3000K on your protein and the 
protein unfolded.I use 600K on my protein but the system explode(box become 
very big  and protein locate at corner).I saw the second way on this 
web,http://manual.gromacs.org/online/protunf.htmlBut it also useless to me.And 
in my mdp I use thermostat and barostat, which means my system is NPT.Does 
anything wrong in my methods?Sincerely yours,Hsin-Lin--gmx-users mailing list   
 gmx-users@gromacs.orghttp://lists.gromacs.org/mailman/listinfo/gmx-usersPlease 
search the archive at http://www.gromacs.org/Support/Mailing_Lists/Search 
before posting!Please don't post (un)subscribe requests to the list. Use thewww 
interface or send it to gmx-users-requ...@gromacs.org.can't post? Read 
http://www.gromacs.org/Support/Mailing_Lists

-- 
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 post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/Support/Mailing_Lists

[gmx-users] Heat Capacity Question

2011-06-16 Thread Fabian Casteblanco
Hello,

I am trying to compare the heat capacity (at NPT) of 1000 molecules of
methanol.  I ran it all to equilibrium and used g_energy and -nmol 1000 to
calculate the heat capacity.  The value I achieved, ~140 J/mol*K is far from
what I see is 79.5 J/mol.  I have read on some past posts that there was
some bug in calculating heat capacity but I'm not sure if that applies.  I
tried plotting Enthalpy vs Temp inorder to get dH/dT since that should be
equal to heat capacity but I'm getting ~140 J/mol*K.  I'm not sure if I need
to include the number of constraints and how would I do that.  I used
constr=all-bonds.  Is there something I'm doing wrong?  when it says 140
J/mol*K, it is per mole of what?  I had calculated heat of vaporization
using Epotg - Epotliq+RT=Hvap and I got it close to the experimental
value (34.09 compared to 35.2 kJ/mol).

Thanks for the help.


-- 
*Best regards,*
**
*Fabian F. Casteblanco*
*Rutgers University -- *
*Chemical Engineering PhD Student*
*C: +908 917 0723*
*E:  **fabian.castebla...@gmail.com* fabian.castebla...@gmail.com
-- 
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 post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/Support/Mailing_Lists

[gmx-users] Including a solvation model after converting from AMBER to GROMACS

2011-06-16 Thread Massih Khorvash
Hello,

I first used Antechamber and Leap to generate the topology and coordinate
files for benzamide.pdb.

I then converted them to gromacs formats and added the water molecules using
the following commands:

editconf -f benz.gro -bt dodecahedron -d 1.2 -o benz_box.gro

genbox -cp benz_box.gro -cs tip4p.gro (there was no tip3p.gro) -p benz.top
-o benz_solvated.gro


*
; benz.top created by rdparm2gmx.pl Wed Jun  8 14:58:07 MDT 2011

[ defaults ]
; nbfunccomb-rule   gen-pairs   fudgeLJ fudgeQQ
1   2   yes 0.5 0.8333

[ atomtypes ]
;name  bond_typemasscharge   ptype  sigma  epsilon
n  n  0.  0.  A   3.25000e-01  7.11280e-01
c  c  0.  0.  A   3.39967e-01  3.59824e-01
haha  0.  0.  A   2.59964e-01  6.27600e-02
hnhn  0.  0.  A   1.06908e-01  6.56888e-02
caca  0.  0.  A   3.39967e-01  3.59824e-01
o  o  0.  0.  A   2.95992e-01  8.78640e-01

[ moleculetype ]
; Namenrexcl
solute 3

[ atoms ]
;   nr   type  resnr residue  atom   cgnr charge   mass
typeBchargeB
 1 ca  1MOL  C  1   -0.07700  12.00
 2 ca  1MOL C1  2   -0.13800  12.00
 3 ca  1MOL C2  3   -0.10900  12.00
 4 ca  1MOL C3  4   -0.13900  12.00
 5 ca  1MOL C4  5   -0.10500  12.00
 6 ca  1MOL C5  6   -0.14160  12.00
 7  c  1MOL C6  70.67070  12.00
 8  o  1MOL  O  8   -0.61010  16.00
 9  n  1MOL  N  9   -0.67400  14.00
10 ha  1MOL  H 100.15600   1.00
11 ha  1MOL H1 110.13800   1.00
12 ha  1MOL H2 120.13500   1.00
13 ha  1MOL H3 130.13600   1.00
14 ha  1MOL H4 140.13300   1.00
15 hn  1MOL H5 150.31750   1.00
16 hn  1MOL H6 160.30750   1.00

[ bonds ]
;  aiaj funct  r  k
110 1  1.0870e-01  2.8811e+05
211 1  1.0870e-01  2.8811e+05
312 1  1.0870e-01  2.8811e+05
413 1  1.0870e-01  2.8811e+05
514 1  1.0870e-01  2.8811e+05
915 1  1.0090e-01  3.4326e+05
916 1  1.0090e-01  3.4326e+05
1 2 1  1.3870e-01  4.0033e+05
1 6 1  1.3870e-01  4.0033e+05
2 3 1  1.3870e-01  4.0033e+05
3 4 1  1.3870e-01  4.0033e+05
4 5 1  1.3870e-01  4.0033e+05
5 6 1  1.3870e-01  4.0033e+05
6 7 1  1.4870e-01  2.9263e+05
7 8 1  1.2140e-01  5.4225e+05
7 9 1  1.3450e-01  4.0016e+05

[ pairs ]
;  aiaj funct
 1 12  1
 1 14  1
 2 13  1
10  3  1
 3 14  1
 4 11  1
 5 12  1
10  5  1
 6 11  1
 6 13  1
 6 15  1
 6 16  1
10  7  1
 7 14  1
 8 15  1
 8 16  1
10 11  1
11 12  1
12 13  1
13 14  1
 1  4  1
 1  8  1
 1  9  1
 2  5  1
 2  7  1
 6  3  1
 4  7  1
 5  8  1
 5  9  1

[ angles ]
;  aiajak funct  theta   cth
1 211 1  1.2001e+02  4.0551e+02
2 110 1  1.2001e+02  4.0551e+02
2 312 1  1.2001e+02  4.0551e+02
3 211 1  1.2001e+02  4.0551e+02
3 413 1  1.2001e+02  4.0551e+02
4 312 1  1.2001e+02  4.0551e+02
4 514 1  1.2001e+02  4.0551e+02
5 413 1  1.2001e+02  4.0551e+02
6 110 1  1.2001e+02  4.0551e+02
6 514 1  1.2001e+02  4.0551e+02
7 915 1  1.1846e+02  4.1179e+02
7 916 1  1.1846e+02  4.1179e+02
   15 916 1  1.1785e+02  3.3246e+02
1 2 3 1  1.1997e+02  5.6216e+02
1 6 5 1  1.1997e+02  5.6216e+02
1 6 7 1  1.2014e+02  5.4091e+02
2 1 6 1  1.1997e+02  5.6216e+02
2 3 4 1  1.1997e+02  5.6216e+02
3 4 5 1  1.1997e+02  5.6216e+02
4 5 6 1  1.1997e+02  5.6216e+02
5 6 7 1  1.2014e+02  5.4091e+02
6 7 8 1  1.2344e+02  5.7463e+02
6 7 9 1  1.1514e+02  5.7296e+02
8 7 9 1  1.2203e+02  6.3455e+02

[ dihedrals ]
;i  j   k  l funcC0  ...  C5
12312 

[gmx-users]gmx4.5.4 genion problem: No line with moleculetype 'SOL' found

2011-06-16 Thread Ye Yang
 Hi, everyone:
   I am a new user of Gromacs, and I am running through the tutorial.
 When I am trying to run the ligand-receptor binding tutorial from

 http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin/gmx-tutorials/complex/02_topology.html
  I met some trouble in adding ion.
  Each time when I use genion, it shows:

 Will try to add 0 NA ions and 6 CL ions.
 Select a continuous group of solvent molecules
 Group 0 ( System) has 33046 elements
 Group 1 (Protein) has  1693 elements
 Group 2 (  Protein-H) has  1301 elements
 Group 3 (C-alpha) has   163 elements
 Group 4 (   Backbone) has   489 elements
 Group 5 (  MainChain) has   653 elements
 Group 6 (   MainChain+Cb) has   805 elements
 Group 7 (MainChain+H) has   815 elements
 Group 8 (  SideChain) has   878 elements
 Group 9 (SideChain-H) has   648 elements
 Group10 (Prot-Masses) has  1693 elements
 Group11 (non-Protein) has 31353 elements
 Group12 (  Other) has15 elements
 Group13 (JZ4) has15 elements
 Group14 (  Water) has 31338 elements
 Group15 (SOL) has 31338 elements
 Group16 (  non-Water) has  1708 elements
 Select a group: 15
 Selected 15: 'SOL'
 Number of (3-atomic) solvent molecules: 10446

 Processing topology

 Back Off! I just backed up temp.top to ./#temp.top.3#

 ---
 Program genion_d, VERSION 4.5.4
 Source code file: gmx_genion.c, line: 285

 Fatal error:
 No line with moleculetype 'SOL' found the [ molecules ] section of file
 'topol.top'

 I checked my topology file and it looks fine:
 ; Include Position restraint file
 #ifdef POSRES
 #include posre.itp
 #endif

 ; Include water topology
 #include gromos43a1.ff/spc.itp

 ; Include ligand topoloty
 #include drg.itp
 #ifdef POSRES_WATER
 ; Position restraint for each water oxygen
 [ position_restraints ]
 ;  i funct   fcxfcyfcz
11   1000   1000   1000
 #endif

 ; Include topology for ions
 #include gromos43a1.ff/ions.itp

 [ system ]
 ; Name

 Protein
 [ molecules ]
 ; Compound#mols
 Protein_chain_A 1
 JZ4 1
 SOL 10446

 One thing that might happen is in the genion source file,  I read through
 it, and the problem either happens in loading the line to buf2, or in the
 gmx_strcmp(buf2,SOL), since clearly the SOL_line is -1 aftermath.

 The other thing I tried is to remove the ligand JZ4 part in both topology
 and coordinates file, and in this case, it works perfectly for adding ion.
 But in this way, I do not know how to insert my ligand into the system since
 it might collide with solvent.

 Can someone help me with this problem?

 Thank you all very much.

 Ye Yang



-- 
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 post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/Support/Mailing_Lists

[gmx-users] g_hbond

2011-06-16 Thread simon sham
Hi,
To help clarify my previous g_hbond question, I have these two groups in my 
index file
[ a_100_a_101]  # 100 is N, #101 is H
   100 101
[ a_200 ]    # 200 is O
   200 

when I ran g_hbond, I got Found 1 donors and 2 acceptors. The part that I 
don't quite understand is why it says 2 acceptors. Did I set up the groups 
correctly?

Again, hope this clarifies my question.
Thanks for your insight.

Simon


-- 
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 post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/Support/Mailing_Lists

[gmx-users] g_hbond

2011-06-16 Thread chris . neale

From g_hbond -h

OH and NH groups are regarded as donors, O is an acceptor always, N  
is an acceptor by default, but this can be switched using -nitacc.  
Dummy hydrogen atoms are assumed to be connected to the first  
preceding non-hydrogen atom.


The two groups are to find hbonds between group A and group B. If you  
want to find hbonds within group A, then you specify an index group  
that contains group A two times.


I think that you want:

[ group ]
100 101 200

Then give group 2x as the identical index group.

PS, there were problems with g_hbond a few months ago reported on  
list, I suggest that you check to ensure that these problems were fixed.


Chris.

-- original message --

Hi,
To help clarify my previous g_hbond question, I have these two groups  
in my index file

[ a_100_a_101]  # 100 is N, #101 is H
   100 101
[ a_200 ]# 200 is O
   200

when I ran g_hbond, I got Found 1 donors and 2 acceptors. The part  
that I don't quite understand is why it says 2 acceptors. Did I set up  
the groups correctly?


Again, hope this clarifies my question.
Thanks for your insight.

Simon



--
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 post (un)subscribe requests to the list. Use the
www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/Support/Mailing_Lists


[gmx-users] Re: Gromacs error

2011-06-16 Thread bharat gupta
Hi,

I have installed gromacs 4.5 on fedora core 15 and whenever I try to run any
command like pdb2gmx ...  I am getting the following error :

error while loading shared libraries: libgmx.so.6: cannot enable executable
stack as shared object requires: Permission denied

Pls help??


-- 
Bharat
-- 
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 post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/Support/Mailing_Lists

RE: [gmx-users] Re: Gromacs error

2011-06-16 Thread Dallas Warren
Someone with some more in depth knowledge will hopefully pipe up, but until 
then there are some things you can read up on and see if it will solve your 
issue.

Google is your friend.

Put error while loading shared libraries: libgmx.so.6: Permission denied into 
it, and  you get a number of topics on forums discussing the issue you have 
had.  Not directly with libgmx.so.6 but that issue appears to be the same.  And 
they all relate to Fedora from what I can see.

This one appears to solve the issue - 
http://forums.fedoraforum.org/showthread.php?t=253870

With some more reading you should be able to sort out what the issue actually 
is, seems to be something to do with how files are handled with Fedora, and 
then the correct course of action.

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 9903 9304
-
When the only tool you own is a hammer, every problem begins to resemble a nail.

From: gmx-users-boun...@gromacs.org [mailto:gmx-users-boun...@gromacs.org] On 
Behalf Of bharat gupta
Sent: Friday, 17 June 2011 10:40 AM
To: Discussion list for GROMACS users
Subject: [gmx-users] Re: Gromacs error

Hi,

I have installed gromacs 4.5 on fedora core 15 and whenever I try to run any 
command like pdb2gmx ...  I am getting the following error :

error while loading shared libraries: libgmx.so.6: cannot enable executable 
stack as shared object requires: Permission denied

Pls help??


--
Bharat
-- 
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 post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/Support/Mailing_Lists

Re: [gmx-users] Gromacs to Amber trajectory

2011-06-16 Thread Mark Abraham

On 17/06/2011 2:57 AM, Raja Pandian wrote:



Dear All,

I have created Amber input file (.prmtop and .rst7) using Xleap. And 
then I have converted this input file into Gromacs input file format 
(*.gro and *.top) with amb2gmx.pl http://amb2gmx.pl script. I did 
simulation using Gromacs Package (Amber Force Field). Now I would like 
to convert Gromacs trajectory file to Amber trajectory file format. So 
that I can easily analyses this trajectory. I tried to convert Gromacs 
trajectory file to Amber trajectory file using VMD. This VMD converted 
trajectory is not working fine. Pls help me in this regards




You really haven't said what you've done at the critical point, nor why 
you think it is not working. So we are unable to help at this time.


Maek
-- 
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 post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/Support/Mailing_Lists

Re: [gmx-users] gmx4.5.4 genion problem: No line with moleculetype 'SOL' found

2011-06-16 Thread Justin A. Lemkul



Ye Yang wrote:

Hi, everyone:
  I am a new user of Gromacs, and I am running through the tutorial. 
When I am trying to run the ligand-receptor binding tutorial from
  
http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin/gmx-tutorials/complex/02_topology.html

 I met some trouble in adding ion.
 Each time when I use genion, it shows:

Will try to add 0 NA ions and 6 CL ions.
Select a continuous group of solvent molecules
Group 0 ( System) has 33046 elements
Group 1 (Protein) has  1693 elements
Group 2 (  Protein-H) has  1301 elements
Group 3 (C-alpha) has   163 elements
Group 4 (   Backbone) has   489 elements
Group 5 (  MainChain) has   653 elements
Group 6 (   MainChain+Cb) has   805 elements
Group 7 (MainChain+H) has   815 elements
Group 8 (  SideChain) has   878 elements
Group 9 (SideChain-H) has   648 elements
Group10 (Prot-Masses) has  1693 elements
Group11 (non-Protein) has 31353 elements
Group12 (  Other) has15 elements
Group13 (JZ4) has15 elements
Group14 (  Water) has 31338 elements
Group15 (SOL) has 31338 elements
Group16 (  non-Water) has  1708 elements
Select a group: 15
Selected 15: 'SOL'
Number of (3-atomic) solvent molecules: 10446

Processing topology

Back Off! I just backed up temp.top to ./#temp.top.3#

---
Program genion_d, VERSION 4.5.4
Source code file: gmx_genion.c, line: 285

Fatal error:
No line with moleculetype 'SOL' found the [ molecules ] section of file 
'topol.top'




Something doesn't match up here - the command backs up temp.top, but then 
complains it can't find SOL in topol.top.  Either genion is looking at the wrong 
file or your command is somehow wrong.



I checked my topology file and it looks fine:
; Include Position restraint file
#ifdef POSRES
#include posre.itp
#endif

; Include water topology
#include gromos43a1.ff/spc.itp

; Include ligand topoloty
#include drg.itp
#ifdef POSRES_WATER
; Position restraint for each water oxygen
[ position_restraints ]
;  i funct   fcxfcyfcz
   11   1000   1000   1000
#endif

; Include topology for ions
#include gromos43a1.ff/ions.itp

[ system ]
; Name

Protein
[ molecules ]
; Compound#mols
Protein_chain_A 1
JZ4 1 
SOL 10446


One thing that might happen is in the genion source file,  I read 
through it, and the problem either happens in loading the line to buf2, 
or in the gmx_strcmp(buf2,SOL), since clearly the SOL_line is -1 
aftermath.


The other thing I tried is to remove the ligand JZ4 part in both 
topology and coordinates file, and in this case, it works perfectly for 
adding ion. But in this way, I do not know how to insert my ligand into 
the system since it might collide with solvent.




If removal of the JZ4 line relieves the problem, perhaps there's a problem with 
that line, i.e. the line ending?  What type of system are you using?  Sometimes 
Windows line endings cause problems.  Always use a plain text editor and use 
dos2unix to process text files if you're on Windows.


Otherwise, the work-around is to not have genion work on the topology.  Make the 
corrections yourself.  Simple addition of ions and subtraction of water 
molecules is trivial.


-Justin


Can someone help me with this problem?

Thank you all very much.

Ye Yang




--


Justin A. Lemkul
Ph.D. Candidate
ICTAS Doctoral Scholar
MILES-IGERT Trainee
Department of Biochemistry
Virginia Tech
Blacksburg, VA
jalemkul[at]vt.edu | (540) 231-9080
http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin


--
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 post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.

Can't post? Read http://www.gromacs.org/Support/Mailing_Lists


Re: [gmx-users] Re: Gromacs error

2011-06-16 Thread bharat gupta
thanks..  I fixed the problem by using the command execstack -c filename


but I have another issue .. I am am preparing the structure for simulation
which is a docked complex containing Phospharylated tyrosine. I am using
Amber 99 force field and update the residuetypes.dat, aminoacid.rtp,
ffbonded.itp files for pTYR but still the error that Residue 'PTYR' not
found in residue topology database ... but all the changes I have made then
where else could the problem be ??


On Fri, Jun 17, 2011 at 10:04 AM, Dallas Warren dallas.war...@monash.eduwrote:

  Someone with some more in depth knowledge will hopefully pipe up, but
 until then there are some things you can read up on and see if it will solve
 your issue.

 ** **

 Google is your friend.

 ** **

 Put “error while loading shared libraries: libgmx.so.6: Permission denied”
 into it, and  you get a number of topics on forums discussing the issue you
 have had.  Not directly with libgmx.so.6 but that issue appears to be the
 same.  And they all relate to Fedora from what I can see.

 ** **

 This one appears to solve the issue -
 http://forums.fedoraforum.org/showthread.php?t=253870

 ** **

 With some more reading you should be able to sort out what the issue
 actually is, seems to be something to do with how files are handled with
 Fedora, and then the correct course of action.

 ** **

 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 9903 9304
 -
 When the only tool you own is a hammer, every problem begins to resemble a
 nail. 

 ** **

 *From:* gmx-users-boun...@gromacs.org [mailto:
 gmx-users-boun...@gromacs.org] *On Behalf Of *bharat gupta
 *Sent:* Friday, 17 June 2011 10:40 AM
 *To:* Discussion list for GROMACS users
 *Subject:* [gmx-users] Re: Gromacs error

 ** **

 Hi,

 I have installed gromacs 4.5 on fedora core 15 and whenever I try to run
 any command like pdb2gmx ...  I am getting the following error :

 error while loading shared libraries: libgmx.so.6: cannot enable executable
 stack as shared object requires: Permission denied

 Pls help??


 --
 Bharat

 --
 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 post (un)subscribe requests to the list. Use the
 www interface or send it to gmx-users-requ...@gromacs.org.
 Can't post? Read http://www.gromacs.org/Support/Mailing_Lists




-- 
Bharat
Ph.D. Candidate
Room No. : 7202A, 2nd Floor
Biomolecular Engineering Laboratory
Division of Chemical Engineering and Polymer Science
Pusan National University
Busan -609735
South Korea
Lab phone no. - +82-51-510-3680, +82-51-583-8343
Mobile no. - 010-5818-3680
E-mail : monu46...@yahoo.com
-- 
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 post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/Support/Mailing_Lists

[gmx-users] load imbalance

2011-06-16 Thread E. Nihal Korkmaz
Hi all,

I am trying to run GB model simulation of a small protein. I keep getting
these errors for every step printed to the log file.

DD  load balancing is limited by minimum cell size in dimension X
DD  step 35999  vol min/aver 0.799! load imb.: force 154.9%

Below i copied my mdp parameters. I'd appreciate any suggestion to help me
fix this.

Thanks,
Nihal


integrator   = sd
tinit= 0
dt   = 0.002
nsteps   = 500
simulation_part  = 1
init_step= 1 %start from 5ns


nstxout  = 5000
nstvout  = 5000
nstenergy= 500
nstxtcout= 500
nstlog   = 500

xtc_grps = System
energygrps   = System
comm_mode= Linear
; neighbor searching and vdw/pme setting up
nstlist  = 10
ns_type  = grid
pbc  = xyz
rlist= 2.0

coulombtype  = pme
fourierspacing   = 0.1
pme_order= 6
rcoulomb = 2.0

vdwtype  = Cut-off
rvdw_switch  = 1.0
rvdw = 2.0

; cpt control
tcoupl   = Berendsen
tc-grps  = System
tau_t= 0.1
ref_t= 300.0
Pcoupl   = Berendsen
pcoupltype   = isotropic
tau_p= 1.0
compressibility  = 4.5e-5
ref_p= 1.0

; velocity  temperature control
gen_vel  = yes
gen_temp = 300.0
annealing= no
constraints  = hbonds
constraint_algorithm = lincs
morse= no


-- 
Elif Nihal Korkmaz

Research Assistant
University of Wisconsin - Biophysics
Member of Qiang Cui  Thomas Record Labs
1101 University Ave, Rm. 8359
Madison, WI 53706
Phone:  608-265-3644
Email:   kork...@wisc.edu
-- 
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 post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/Support/Mailing_Lists

Re: [gmx-users] load imbalance

2011-06-16 Thread Mark Abraham

On 17/06/2011 12:44 PM, E. Nihal Korkmaz wrote:

Hi all,

I am trying to run GB model simulation of a small protein. I keep 
getting these errors for every step printed to the log file.


DD  load balancing is limited by minimum cell size in dimension X
DD  step 35999  vol min/aver 0.799! load imb.: force 154.9%


Not all systems can efficiently parallelize on arbitrary numbers of 
processors for a given implementation. There's an analysis at the top of 
the .log file that describes the issues leading to the minimum cell 
size. Possibly there's an issue there, but more likely you've got not 
enough work for your processors.


Mark

Below i copied my mdp parameters. I'd appreciate any suggestion to 
help me fix this.


Thanks,
Nihal


integrator   = sd
tinit= 0
dt   = 0.002
nsteps   = 500
simulation_part  = 1
init_step= 1 %start from 5ns


nstxout  = 5000
nstvout  = 5000
nstenergy= 500
nstxtcout= 500
nstlog   = 500

xtc_grps = System
energygrps   = System
comm_mode= Linear
; neighbor searching and vdw/pme setting up
nstlist  = 10
ns_type  = grid
pbc  = xyz
rlist= 2.0

coulombtype  = pme
fourierspacing   = 0.1
pme_order= 6
rcoulomb = 2.0

vdwtype  = Cut-off
rvdw_switch  = 1.0
rvdw = 2.0

; cpt control
tcoupl   = Berendsen
tc-grps  = System
tau_t= 0.1
ref_t= 300.0
Pcoupl   = Berendsen
pcoupltype   = isotropic
tau_p= 1.0
compressibility  = 4.5e-5
ref_p= 1.0

; velocity  temperature control
gen_vel  = yes
gen_temp = 300.0
annealing= no
constraints  = hbonds
constraint_algorithm = lincs
morse= no


--
Elif Nihal Korkmaz

Research Assistant
University of Wisconsin - Biophysics
Member of Qiang Cui  Thomas Record Labs
1101 University Ave, Rm. 8359
Madison, WI 53706
Phone:  608-265-3644
Email: kork...@wisc.edu mailto:kork...@wisc.edu




-- 
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 post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/Support/Mailing_Lists

[gmx-users] monitoring the water molecules

2011-06-16 Thread Ramachandran G
Hi gmx-users.
I am doing MD simulation with the green fluorescent protein, where
confined water molecules exist inside the protein barrel.
During the coarse of the simulation i noticed these confined water molecule
get exchange with outside bulk water.

Is it possible to monitor or mark the water molecules during the
simulations? For my analysis it would be very much useful.
Can anyone guide me? Thanks

Rama
-- 
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 post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/Support/Mailing_Lists

Re: [gmx-users] Gromacs to Amber trajectory

2011-06-16 Thread Raja Pandian
Dear All,

I have created Amber input file (.prmtop and .rst7) using Xleap. And then I
have converted this input file into Gromacs input file format (*.gro and
*.top) with amb2gmx.pl script. I did simulation using Gromacs Package (Amber
Force Field). Now I would like to convert Gromacs trajectory file to Amber
trajectory file format. So that I can easily analyses this trajectory. I
tried to convert Gromacs (.xtc) trajectory file to Amber (.mdcrd) trajectory
file using VMD. This VMD converted trajectory is not working fine for DCCM
(Dynamic cross correlation Matrix) analysis.

Now i'm trying to convert Gromacs trajectory into dcd format, and then use
ptraj to read the dcd file and write the mdcrd  (or AMBER netcdf)
trajectory.



Example:

---

 [1] Gromacs20ns.xtc à nojump.xtc (using trjcat/trjconv)

[2] nojump.xtc à amber20ns.dcd (using VMD)

[3] Ptraj file conversion

trajin amber20ns.dcd

trajout amber20ns.nc netcdf

(or)

Is it necessary to do image for this file.



trajin test.dcd

center origin

image origin center

trajout amber20ns.nc netcdf



Is this correct procedure?



I want to do DCCM (Dynamic cross correlation analysis).

DCCM Ptraj file



trajin test.rst7

trajin amber20ns.nc

strip :WAT

strip :Cl-

center :1-85 origin

image origin center

matrix correl name corr @CA out amber20ns.dccm byres



If you have any ptraj file format for this.

 Actually i did simulation based on NPT ensemble. Is it necessary to mention
the BOX dimensions in any were in the dcd file creation?

Already i tried to convert Gromacs trajectory into PDB (total trajectory)
but it is not working well for the analysis purpose of AMBER.


Eagerly waiting for your suggestion

Regards
Raja



On Fri, Jun 17, 2011 at 6:46 AM, Mark Abraham mark.abra...@anu.edu.auwrote:

 **
 On 17/06/2011 2:57 AM, Raja Pandian wrote:



 Dear All,

 I have created Amber input file (.prmtop and .rst7) using Xleap. And then I
 have converted this input file into Gromacs input file format (*.gro and
 *.top) with amb2gmx.pl script. I did simulation using Gromacs Package
 (Amber Force Field). Now I would like to convert Gromacs trajectory file to
 Amber trajectory file format. So that I can easily analyses this trajectory.
 I tried to convert Gromacs trajectory file to Amber trajectory file using
 VMD. This VMD converted trajectory is not working fine. Pls help me in this
 regards


 You really haven't said what you've done at the critical point, nor why you
 think it is not working. So we are unable to help at this time.

 Maek

 --
 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 post (un)subscribe requests to the list. Use the
 www interface or send it to gmx-users-requ...@gromacs.org.
 Can't post? Read http://www.gromacs.org/Support/Mailing_Lists

-- 
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 post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/Support/Mailing_Lists

Re: [gmx-users] monitoring the water molecules

2011-06-16 Thread Mark Abraham

On 17/06/2011 12:55 PM, Ramachandran G wrote:

Hi gmx-users.
I am doing MD simulation with the green fluorescent protein, where 
confined water molecules exist inside the protein barrel.
During the coarse of the simulation i noticed these confined water 
molecule get exchange with outside bulk water.


Is it possible to monitor or mark the water molecules during the 
simulations? For my analysis it would be very much useful.

Can anyone guide me? Thanks


They have the same atom number throughout. Whatever you want to do can 
be done using them, possibly via g_select creating dynamic selection groups.


Mark
--
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 post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.

Can't post? Read http://www.gromacs.org/Support/Mailing_Lists


[gmx-users] Re: gmx-users Digest, Vol 86, Issue 97

2011-06-16 Thread ITHAYARAJA
Dear Sir,
Above error was solved but the system has raised a following new error,
i.e.,

Fatal error:
52 atoms are not part of any of the T-Coupling groups

How do I find which atom file has problem?
help me to solve this problem...



-- 
**
Ithayaraja M,
Research Scholar,
Department of Bionformatics,
Bharathiar University,
Coimbatore 641 046,
Tamil Nadu
India
-- 
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 post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/Support/Mailing_Lists

[gmx-users] cross correlations

2011-06-16 Thread E. Nihal Korkmaz
Dear all,

Is there any built in function that gives me the pairwise correlation of the
fluctuation (unit vector between two coordinates of a residue) of residues
(averaged over the input trajectory). I tried g_covar but that's not what
i'm looking for. The result I want should be an NxN matrix with values
ranging from -1 to 1.

Thanks in advance,
Nihal

-- 
Elif Nihal Korkmaz

Research Assistant
University of Wisconsin - Biophysics
Member of Qiang Cui  Thomas Record Labs
1101 University Ave, Rm. 8359
Madison, WI 53706
Phone:  608-265-3644
Email:   kork...@wisc.edu
-- 
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 post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/Support/Mailing_Lists