Re: [gmx-users] special molecule

2006-06-13 Thread Tamas Horvath
Does it mean, that the itp file, generated by the PRODRG server should only used with  the gmx forcefield, or may it be combined with the other gromos frocefields, like the G43a2?On 6/13/06, 
Mark Abraham <[EMAIL PROTECTED]> wrote:
Dongsheng Zhang wrote:> Sear Mark,>> Since PRODRG server only generates a itp file by using gromos type force> field, do I have to use gromos force field for my whole protein-ligand> system if my ligand force field parameters are got from PRODRG? In
> another word, can I mix two different force field (one for protein, one> for ligand) in one system? Thanks!It's technically possible, if fiddly, but I'd strongly advise againstit. Force fields are mathematical constructs that are optimized to
produce reasonable model physics as an interacting unit. There is noreason to suppose that a bond stretching parameter is a measure of bondstrength, for example, because the other parameters in this force field
might have some quirk that this parameter is compensating for because ofthe optimization procedure. A force field is merely a point in ahigh-dimensional space which when used the right way has been shown to
produce a useful physical model. Thus there is no reason to suppose thata few of those high-dimensional space coordinates plus some others fromsome other force field will combine to produce a useful physical model.
Mark___gmx-users mailing listgmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-usersPlease don't post (un)subscribe requests to the list. Use thewww interface or send it to [EMAIL PROTECTED]
.Can't post? Read http://www.gromacs.org/mailing_lists/users.php
___
gmx-users mailing listgmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to [EMAIL PROTECTED]
Can't post? Read http://www.gromacs.org/mailing_lists/users.php

Re: [gmx-users] special molecule

2006-06-13 Thread Tamas Horvath
What I asked is, if I have to reorganise the pdb file, or renumber the residues, etc.On 6/13/06, Mark Abraham <
[EMAIL PROTECTED]> wrote:Tamas Horvath wrote:> So if I generate an .itp with PRODRG, and link it in the .top file, then
> grompp can use the original .pdb file? Or it's not that simple?The order of the [ molecules ] section will have to correspond to theorder in the structure file you supply with -c to grompp, but if you
read the grompp manpage it is obvious that the format of that structurefile is irrelevant.Mark___gmx-users mailing list
gmx-users@gromacs.orghttp://www.gromacs.org/mailman/listinfo/gmx-usersPlease don't post (un)subscribe requests to the list. Use thewww interface or send it to 
[EMAIL PROTECTED].Can't post? Read http://www.gromacs.org/mailing_lists/users.php

___
gmx-users mailing listgmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to [EMAIL PROTECTED]
Can't post? Read http://www.gromacs.org/mailing_lists/users.php

Re: [gmx-users] special molecule

2006-06-12 Thread Mark Abraham

Dongsheng Zhang wrote:

Sear Mark,

Since PRODRG server only generates a itp file by using gromos type force
field, do I have to use gromos force field for my whole protein-ligand
system if my ligand force field parameters are got from PRODRG? In
another word, can I mix two different force field (one for protein, one
for ligand) in one system? Thanks!


It's technically possible, if fiddly, but I'd strongly advise against 
it. Force fields are mathematical constructs that are optimized to 
produce reasonable model physics as an interacting unit. There is no 
reason to suppose that a bond stretching parameter is a measure of bond 
strength, for example, because the other parameters in this force field 
might have some quirk that this parameter is compensating for because of 
the optimization procedure. A force field is merely a point in a 
high-dimensional space which when used the right way has been shown to 
produce a useful physical model. Thus there is no reason to suppose that 
a few of those high-dimensional space coordinates plus some others from 
some other force field will combine to produce a useful physical model.


Mark
___
gmx-users mailing listgmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to [EMAIL PROTECTED]

Can't post? Read http://www.gromacs.org/mailing_lists/users.php


Re: [gmx-users] special molecule

2006-06-12 Thread Mark Abraham

Tamas Horvath wrote:
So if I generate an .itp with PRODRG, and link it in the .top file, then 
grompp can use the original .pdb file? Or it's not that simple?


The order of the [ molecules ] section will have to correspond to the 
order in the structure file you supply with -c to grompp, but if you 
read the grompp manpage it is obvious that the format of that structure 
file is irrelevant.


Mark
___
gmx-users mailing listgmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to [EMAIL PROTECTED]

Can't post? Read http://www.gromacs.org/mailing_lists/users.php


Re: [gmx-users] special molecule

2006-06-12 Thread Florian Haberl
Hi,

On Monday 12 June 2006 17:25, Dongsheng Zhang wrote:
> Sear Mark,
>
> Since PRODRG server only generates a itp file by using gromos type force
> field, do I have to use gromos force field for my whole protein-ligand
> system if my ligand force field parameters are got from PRODRG? In
> another word, can I mix two different force field (one for protein, one
> for ligand) in one system? Thanks!

Mixing ff types is a bad idea, i don`t know if PRODRG supports other ff like 
the old and depracted gmx format.

I hope that there will be soon a PRODRG version which will support other force 
field formats, also amber ff in gromacs without heavy modifing the pdb file.

As a hint: Take a look at 
http://www.gromacs.org/gromacs/documentation/tutorial.html

for builing a system of ligand and protein.

>
>
> Best Wishes!
>
> Dongsheng
>
> On Tue, 2006-06-13 at 01:16 +1000, Mark Abraham wrote:
> > Tamas Horvath wrote:
> > > As I understand, if there are "special" molecules in a pdb file,
> > > pdb2gmx cannot convert it. However, PRODRG can create an *.itp file for
> > > that molecule, so that I can include it in the generated topology file.
> > > But how can I place the molecule in it's original position? Or even
> > > better, how can I convert the original ("full") .pdb file for gromacs?
> >
> > Where's the problem? pdb2gmx makes a .top file for your non-special
> > system, you #include the .itp file for your special part. There's
> > nothing about atomic positions yet have a look at a .top file. Now
> > grompp takes the .top and some file with atomic positions and makes a
> > .tpr file. Same general procedure as normal.
> >
> > Mark
> > ___
> > gmx-users mailing listgmx-users@gromacs.org
> > http://www.gromacs.org/mailman/listinfo/gmx-users
> > Please don't post (un)subscribe requests to the list. Use the
> > www interface or send it to [EMAIL PROTECTED]
> > Can't post? Read http://www.gromacs.org/mailing_lists/users.php
>
> ___
> gmx-users mailing listgmx-users@gromacs.org
> http://www.gromacs.org/mailman/listinfo/gmx-users
> Please don't post (un)subscribe requests to the list. Use the
> www interface or send it to [EMAIL PROTECTED]
> Can't post? Read http://www.gromacs.org/mailing_lists/users.php

Greetings,

Florian

-- 
---
 Florian Haberl
 Computer-Chemie-Centrum   
 Universitaet Erlangen/ Nuernberg
 Naegelsbachstr 25
 D-91052 Erlangen
 Mailto: florian.haberl AT chemie.uni-erlangen.de
---
___
gmx-users mailing listgmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to [EMAIL PROTECTED]
Can't post? Read http://www.gromacs.org/mailing_lists/users.php


Re: [gmx-users] special molecule

2006-06-12 Thread Tamas Horvath
So if I generate an .itp with PRODRG, and link it in the .top file, then grompp can use the original .pdb file? Or it's not that simple?On 6/12/06, Alan Dodd
 <[EMAIL PROTECTED]> wrote:It's entirely possible to mix and match forcefield
files to create a hybrid forcefield - the lipidforcefield I use has OPLS parameters for theheadgroups, custom parameters for lipid tails, and GMXparameters for everything else.  As long as the hybridforcefield is internally consistent, it will work.
How accurate it is, is another matter...So yes, you can technically use PRODRG itps with otherforcefields, but bear in mind firstly, atomtypes etcin your itp file will all need checking to make surethey are accounted for in both forcefields (are
'default' values the same in one forcefield asanother, for instance?  Probably not - you may need toalter the itp to describe defaults explicitly), andsecondly that they are not designed to be used in the
same system.--- Dongsheng Zhang <[EMAIL PROTECTED]>wrote:> Sear Mark,>> Since PRODRG server only generates a itp file by
> using gromos type force> field, do I have to use gromos force field for my> whole protein-ligand> system if my ligand force field parameters are got> from PRODRG? In> another word, can I mix two different force field
> (one for protein, one> for ligand) in one system? Thanks!>>> Best Wishes!>> Dongsheng On Tue, 2006-06-13 at 01:16 +1000, Mark Abraham> wrote:
> > Tamas Horvath wrote:> > > As I understand, if there are "special"> molecules in a pdb file, pdb2gmx> > > cannot convert it. However, PRODRG can create an> *.itp file for that
> > > molecule, so that I can include it in the> generated topology file. But> > > how can I place the molecule in it's original> position? Or even better,> > > how can I convert the original ("full") .pdb
> file for gromacs?> >> > Where's the problem? pdb2gmx makes a .top file for> your non-special> > system, you #include the .itp file for your> special part. There's> > nothing about atomic positions yet have a look
> at a .top file. Now> > grompp takes the .top and some file with atomic> positions and makes a> > .tpr file. Same general procedure as normal.> >> > Mark> > ___
> > gmx-users mailing listgmx-users@gromacs.org> > http://www.gromacs.org/mailman/listinfo/gmx-users
> > Please don't post (un)subscribe requests to the> list. Use the> > www interface or send it to> [EMAIL PROTECTED].
> > Can't post? Read> http://www.gromacs.org/mailing_lists/users.php> ___> gmx-users mailing list
gmx-users@gromacs.org> http://www.gromacs.org/mailman/listinfo/gmx-users> Please don't post (un)subscribe requests to the
> list. Use the> www interface or send it to> [EMAIL PROTECTED].> Can't post? Read> 
http://www.gromacs.org/mailing_lists/users.php>__Do You Yahoo!?Tired of spam?  Yahoo! Mail has the best spam protection around
http://mail.yahoo.com___gmx-users mailing listgmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-usersPlease don't post (un)subscribe requests to the list. Use thewww interface or send it to [EMAIL PROTECTED]
.Can't post? Read http://www.gromacs.org/mailing_lists/users.php
___
gmx-users mailing listgmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to [EMAIL PROTECTED]
Can't post? Read http://www.gromacs.org/mailing_lists/users.php

Re: [gmx-users] special molecule

2006-06-12 Thread Alan Dodd
It's entirely possible to mix and match forcefield
files to create a hybrid forcefield - the lipid
forcefield I use has OPLS parameters for the
headgroups, custom parameters for lipid tails, and GMX
parameters for everything else.  As long as the hybrid
forcefield is internally consistent, it will work. 
How accurate it is, is another matter...
So yes, you can technically use PRODRG itps with other
forcefields, but bear in mind firstly, atomtypes etc
in your itp file will all need checking to make sure
they are accounted for in both forcefields (are
'default' values the same in one forcefield as
another, for instance?  Probably not - you may need to
alter the itp to describe defaults explicitly), and
secondly that they are not designed to be used in the
same system.

--- Dongsheng Zhang <[EMAIL PROTECTED]>
wrote:

> Sear Mark,
> 
> Since PRODRG server only generates a itp file by
> using gromos type force
> field, do I have to use gromos force field for my
> whole protein-ligand
> system if my ligand force field parameters are got
> from PRODRG? In
> another word, can I mix two different force field
> (one for protein, one
> for ligand) in one system? Thanks!
> 
> 
> Best Wishes!
> 
> Dongsheng
> 
> 
> 
> On Tue, 2006-06-13 at 01:16 +1000, Mark Abraham
> wrote:
> > Tamas Horvath wrote:
> > > As I understand, if there are "special"
> molecules in a pdb file, pdb2gmx 
> > > cannot convert it. However, PRODRG can create an
> *.itp file for that 
> > > molecule, so that I can include it in the
> generated topology file. But 
> > > how can I place the molecule in it's original
> position? Or even better, 
> > > how can I convert the original ("full") .pdb
> file for gromacs?
> > 
> > Where's the problem? pdb2gmx makes a .top file for
> your non-special 
> > system, you #include the .itp file for your
> special part. There's 
> > nothing about atomic positions yet have a look
> at a .top file. Now 
> > grompp takes the .top and some file with atomic
> positions and makes a 
> > .tpr file. Same general procedure as normal.
> > 
> > Mark
> > ___
> > gmx-users mailing listgmx-users@gromacs.org
> > http://www.gromacs.org/mailman/listinfo/gmx-users
> > Please don't post (un)subscribe requests to the
> list. Use the 
> > www interface or send it to
> [EMAIL PROTECTED]
> > Can't post? Read
> http://www.gromacs.org/mailing_lists/users.php
> ___
> gmx-users mailing listgmx-users@gromacs.org
> http://www.gromacs.org/mailman/listinfo/gmx-users
> Please don't post (un)subscribe requests to the
> list. Use the 
> www interface or send it to
> [EMAIL PROTECTED]
> Can't post? Read
> http://www.gromacs.org/mailing_lists/users.php
> 


__
Do You Yahoo!?
Tired of spam?  Yahoo! Mail has the best spam protection around 
http://mail.yahoo.com 
___
gmx-users mailing listgmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
Please don't post (un)subscribe requests to the list. Use the
www interface or send it to [EMAIL PROTECTED]
Can't post? Read http://www.gromacs.org/mailing_lists/users.php


Re: [gmx-users] special molecule

2006-06-12 Thread Dongsheng Zhang
Sear Mark,

Since PRODRG server only generates a itp file by using gromos type force
field, do I have to use gromos force field for my whole protein-ligand
system if my ligand force field parameters are got from PRODRG? In
another word, can I mix two different force field (one for protein, one
for ligand) in one system? Thanks!


Best Wishes!

Dongsheng



On Tue, 2006-06-13 at 01:16 +1000, Mark Abraham wrote:
> Tamas Horvath wrote:
> > As I understand, if there are "special" molecules in a pdb file, pdb2gmx 
> > cannot convert it. However, PRODRG can create an *.itp file for that 
> > molecule, so that I can include it in the generated topology file. But 
> > how can I place the molecule in it's original position? Or even better, 
> > how can I convert the original ("full") .pdb file for gromacs?
> 
> Where's the problem? pdb2gmx makes a .top file for your non-special 
> system, you #include the .itp file for your special part. There's 
> nothing about atomic positions yet have a look at a .top file. Now 
> grompp takes the .top and some file with atomic positions and makes a 
> .tpr file. Same general procedure as normal.
> 
> Mark
> ___
> gmx-users mailing listgmx-users@gromacs.org
> http://www.gromacs.org/mailman/listinfo/gmx-users
> Please don't post (un)subscribe requests to the list. Use the 
> www interface or send it to [EMAIL PROTECTED]
> Can't post? Read http://www.gromacs.org/mailing_lists/users.php
___
gmx-users mailing listgmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to [EMAIL PROTECTED]
Can't post? Read http://www.gromacs.org/mailing_lists/users.php


Re: [gmx-users] special molecule

2006-06-12 Thread Dongsheng Zhang
I think you can use editconf to convert your original pdb to a gro file.


On Mon, 2006-06-12 at 13:48 +, Tamas Horvath wrote:
> As I understand, if there are "special" molecules in a pdb file,
> pdb2gmx cannot convert it. However, PRODRG can create an *.itp file
> for that molecule, so that I can include it in the generated topology
> file. But how can I place the molecule in it's original position? Or
> even better, how can I convert the original ("full") .pdb file for
> gromacs? 
> ___
> gmx-users mailing listgmx-users@gromacs.org
> http://www.gromacs.org/mailman/listinfo/gmx-users
> Please don't post (un)subscribe requests to the list. Use the 
> www interface or send it to [EMAIL PROTECTED]
> Can't post? Read http://www.gromacs.org/mailing_lists/users.php
___
gmx-users mailing listgmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to [EMAIL PROTECTED]
Can't post? Read http://www.gromacs.org/mailing_lists/users.php


Re: [gmx-users] special molecule

2006-06-12 Thread Mark Abraham

Tamas Horvath wrote:
As I understand, if there are "special" molecules in a pdb file, pdb2gmx 
cannot convert it. However, PRODRG can create an *.itp file for that 
molecule, so that I can include it in the generated topology file. But 
how can I place the molecule in it's original position? Or even better, 
how can I convert the original ("full") .pdb file for gromacs?


Where's the problem? pdb2gmx makes a .top file for your non-special 
system, you #include the .itp file for your special part. There's 
nothing about atomic positions yet have a look at a .top file. Now 
grompp takes the .top and some file with atomic positions and makes a 
.tpr file. Same general procedure as normal.


Mark
___
gmx-users mailing listgmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to [EMAIL PROTECTED]

Can't post? Read http://www.gromacs.org/mailing_lists/users.php