Re: 答复: [gmx-users] topology file of ligand

2013-10-15 Thread Justin Lemkul



On 10/15/13 12:25 PM, sunyeping wrote:

Thank you Dr. Lemkul,

Could you recommand some primary literatures which are most usefule for me to 
understand the force field?



The ones cited in the Gromacs manual are a good start.  One should never attempt 
to use a force field without scrutinizing its history, assumptions, and 
parametrization methods...


-Justin

--
==

Justin A. Lemkul, Ph.D.
Postdoctoral Fellow

Department of Pharmaceutical Sciences
School of Pharmacy
Health Sciences Facility II, Room 601
University of Maryland, Baltimore
20 Penn St.
Baltimore, MD 21201

jalem...@outerbanks.umaryland.edu | (410) 706-7441

==
--
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] topology file of ligand

2013-10-15 Thread sunyeping
Thank you Dr. Lemkul,
 
Could you recommand some primary literatures which are most usefule for me to 
understand the force field? 

Yeping Sun
Institute of Microbiology, Chinese Academy of Sciences

--发件人:Justin 
Lemkul 发送时间:2013年10月16日(星期三) 00:05收件人:gromacs 
主 题:Re: [gmx-users] topology file of ligandOn 10/15/13 
11:24 AM, sunyeping wrote:> Dear gromacs users,>> I am study the tutorial for 
protein-ligand complex system (gromacs tutorial 5: T4 lysozyme in complex with 
JZ4). The tutorial use PRODRG to produce the topology file. The author state 
that he immediately notice three things that are wrong with this topology, 
which include charge group, H atom charge, etc.> And the author then 
reconstructed a topology file. I don't understand how the charge group and 
charge of individual atom are correctly assigned. Could you explain it to 
me?>The Gromos96 parameter sets assign charges and charge groups based on 
chemical functional groups that are transferable between molecules. Spotting 
inaccuracies and knowing how to fix them come with understanding how the force 
field was derived (read the primary literature and study existing 
molecules).-Justin-- ==Justin 
A. Lemkul, Ph.D.Postdoctoral FellowDepartment of Pharmaceutical SciencesSchool 
of PharmacyHealth Sciences Facility II, Room 601University of Maryland, 
Baltimore20 Penn St.Baltimore, MD 21201jalem...@outerbanks.umaryland.edu | 
(410) 706-7441==-- gmx-users 
mailing list 
gmx-users@gromacs.orghttp://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] topology file of ligand

2013-10-15 Thread Justin Lemkul



On 10/15/13 11:24 AM, sunyeping wrote:

Dear gromacs users,

I am study the tutorial for protein-ligand complex system (gromacs tutorial 5: 
T4 lysozyme in complex with JZ4). The tutorial use PRODRG to produce the 
topology file. The author state that he  immediately notice three things that 
are wrong with this topology, which include charge group, H atom charge, etc.
  And the author then reconstructed a topology file.  I don't understand 
how the charge group and charge of individual atom are correctly assigned. 
Could you explain it to me?



The Gromos96 parameter sets assign charges and charge groups based on chemical 
functional groups that are transferable between molecules.  Spotting 
inaccuracies and knowing how to fix them come with understanding how the force 
field was derived (read the primary literature and study existing molecules).


-Justin

--
==

Justin A. Lemkul, Ph.D.
Postdoctoral Fellow

Department of Pharmaceutical Sciences
School of Pharmacy
Health Sciences Facility II, Room 601
University of Maryland, Baltimore
20 Penn St.
Baltimore, MD 21201

jalem...@outerbanks.umaryland.edu | (410) 706-7441

==
--
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] topology file of ligand

2013-10-15 Thread sunyeping
Dear gromacs users,
 
I am study the tutorial for protein-ligand complex system (gromacs tutorial 5: 
T4 lysozyme in complex with JZ4). The tutorial use PRODRG to produce the 
topology file. The author state that he  immediately notice three things that 
are wrong with this topology, which include charge group, H atom charge, etc. 
     And the author then reconstructed a topology file.  I don't understand how 
the charge group and charge of individual atom are correctly assigned. Could 
you explain it to me? 
 
Thanks.
 
Yeping Sun

Institute of Microbiology, Chinese Academy of Sciences

--
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] topology and coordinate file not matching after grompp

2013-07-31 Thread Justin Lemkul



On 7/31/13 6:04 AM, chinnu657 wrote:

Hello,

I have a problem and although I've read almost everything the internet could
offer on how differently people have approached this problem, mine hasn't
been solved yet.

Basically, I am trying to mix 2 different proteins in one box and after
doing my grompp to neutralise the system, I get the error number of
coordinates in coordinate file (solvated.gro, 85059) does not match topology
(topol.top, 85079).

This is how my topology file looks:
;GROMACS toplogy

;
;Include the force field
#include "oplsaa.ff/forcefield.itp"


; Include chain topologies
#include "topol_1AKI.itp"
#include "topol_2CDS.itp"

;Include water topology
#include "spce.itp"

;Include generic ion topology
#include "ions.itp"

[ system ]
Two proteins in water

[ molecules ]
Protein_1AKI1
Protein_2CDS1
SOL 27053

How can I deduce from this data that the number of coordinates in topol.top
is 85079. And what exactly should I do to reduce the number of molecules
here so that the 2 files have the same number of coordinates.



Number of atoms in Protein_1AKI [moleculetype] + number of atoms in Protein_2CDS 
[moleculetype] + (3*27053) is the number of atoms that the topology specifies. 
Cross-check that against the coordinate file.  This error always arises for the 
same reason - incorrect bookkeeping of some sort.  You haven't said how you 
built the system or what commands you gave, for instance, to genbox, so it's 
hard to be more specific.


-Justin

--
==

Justin A. Lemkul, Ph.D.
Postdoctoral Fellow

Department of Pharmaceutical Sciences
School of Pharmacy
Health Sciences Facility II, Room 601
University of Maryland, Baltimore
20 Penn St.
Baltimore, MD 21201

jalem...@outerbanks.umaryland.edu | (410) 706-7441

==
--
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] topology and coordinate file not matching after grompp

2013-07-31 Thread chinnu657
Hello,

I have a problem and although I've read almost everything the internet could
offer on how differently people have approached this problem, mine hasn't
been solved yet.

Basically, I am trying to mix 2 different proteins in one box and after
doing my grompp to neutralise the system, I get the error number of
coordinates in coordinate file (solvated.gro, 85059) does not match topology
(topol.top, 85079). 

This is how my topology file looks:
;GROMACS toplogy

;
;Include the force field
#include "oplsaa.ff/forcefield.itp"


; Include chain topologies
#include "topol_1AKI.itp"
#include "topol_2CDS.itp"

;Include water topology
#include "spce.itp"

;Include generic ion topology
#include "ions.itp"

[ system ]
Two proteins in water

[ molecules ]
Protein_1AKI1
Protein_2CDS1
SOL 27053

How can I deduce from this data that the number of coordinates in topol.top
is 85079. And what exactly should I do to reduce the number of molecules
here so that the 2 files have the same number of coordinates. 

I will really appreciate the help. I am pretty new with MD simulations

Kind Regards,
Chinnu



--
View this message in context: 
http://gromacs.5086.x6.nabble.com/topology-and-coordinate-file-not-matching-after-grompp-tp5010221.html
Sent from the GROMACS Users Forum mailing list archive at Nabble.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


Re: [gmx-users] Topology file-residue unknown

2013-04-06 Thread Mark Abraham
On Sat, Apr 6, 2013 at 5:56 AM, Juliette N.  wrote:

> Thank you Justin. I have now the topology. I have a quick question
> regarding atom types. I used opls_143 and opls_144 for C and H of ethylene
> respec in rtp. However, in VMD I dont see double bonds C=C. I though maybe
> these are not the proper atomtypes for H2C=CH2.
>

http://www.gromacs.org/Downloads/Related_Software/Visualization_Software#Topology_bonds_vs_Rendered_bonds

Mark


>
> rtp entry:
>
> [ atoms ]
>C1opls_143-0.1201
>H11   opls_144 0.0601
>H12   opls_144 0.0601
>C2opls_143-0.1202
>H21   opls_144 0.0602
>H22   opls_144 0.0602
>
> >From atomtypes.atp
>
> opls_141   12.01100  ; alkene C (R2-C=)
>  opls_142   12.01100  ; alkene C (RH-C=)
> * opls_143   12.01100  ; alkene C (H2-C=)
>  opls_1441.00800  ; alkene H (H-C=)*
>
> Did I select the correct atomtypes?
> Thanks many times!
>
>
>
>
> On 4 April 2013 18:23, Justin Lemkul  wrote:
>
> >
> >
> > On 4/4/13 6:17 PM, Juliette N. wrote:
> >
> >> Hi Justin,
> >>
> >> Thanks a lot for your message. I am petrified why pdb2gmx is not
> >> recognizing the residue ETY. This the first residue added to
> ffoplsaa.rtp
> >> but as you may see below it is not read from rtp.
> >>
> >> I am sure this residue is added to ffoplsaa.rtp which is existing in the
> >> working directory:
> >>
> >> [ ETY ]
> >>   [ atoms ]
> >> C1opls_143-0.1201
> >> H11   opls_144 0.0601
> >> H12   opls_144 0.0601
> >> C2opls_143-0.1202
> >> H21   opls_144 0.0602
> >> H22   opls_144 0.0602
> >>
> >> [ bonds ]
> >> C1H11
> >> C1H12
> >> C1C2
> >> C2H21
> >> C2H22
> >>;C2C1
> >>
> >> Then I issue:
> >>
> >> pdb2gmx -f Ethylene.gro -o GRO.gro -p Ethylene.top -ff oplsaa
> >>   :-)  G  R  O  M  A  C  S  (-:
> >>
> >>GROwing Monsters And Cloning Shrimps
> >>
> >>  :-)  VERSION 4.5.4  (-:
> >>
> >> Using the Oplsaa force field in directory oplsaa.ff
> >>
> >> Opening force field file
> >> /usr/local/gromacs/share/**gromacs/top/oplsaa.ff/**watermodels.dat
> >>
> >> Select the Water Model:
> >>   1: TIP4P  TIP 4-point, recommended
> >>   2: TIP3P  TIP 3-point
> >>   3: TIP5P  TIP 5-point
> >>   4: SPCsimple point charge
> >>   5: SPC/E  extended simple point charge
> >>   6: None
> >> 6
> >> Opening force field file
> >> /usr/local/gromacs/share/**gromacs/top/oplsaa.ff/**aminoacids.r2b
> >> Reading Ethylene.gro...
> >> Read 'Giant Rising Ordinary Mutants for A Clerical Setup', 6 atoms
> >> Analyzing pdb file
> >> Splitting PDB chains based on TER records or changing chain id.
> >> There are 1 chains and 0 blocks of water and 1 residues with 6 atoms
> >>
> >>chain  #res #atoms
> >>1 ' ' 1  6
> >>
> >> No occupancies in Ethylene.gro
> >> Opening force field file
> >> /usr/local/gromacs/share/**gromacs/top/oplsaa.ff/**atomtypes.atp
> >> Atomtype 1
> >> Reading residue database... (oplsaa)
> >> Opening force field file
> >> /usr/local/gromacs/share/**gromacs/top/oplsaa.ff/**aminoacids.rtp
> >> Residue 56
> >> Sorting it all out...
> >> Opening force field file
> >> /usr/local/gromacs/share/**gromacs/top/oplsaa.ff/**aminoacids.hdb
> >> Opening force field file
> >> /usr/local/gromacs/share/**gromacs/top/oplsaa.ff/**aminoacids.n.tdb
> >> Opening force field file
> >> /usr/local/gromacs/share/**gromacs/top/oplsaa.ff/**aminoacids.c.tdb
> >>
> >> Back Off! I just backed up Ethylene.top to ./#Ethylene.top.1#
> >> Processing chain 1 (6 atoms, 1 residues)
> >> There are 0 donors and 0 acceptors
> >> There are 0 hydrogen bonds
> >> Warning: Starting residue ETY1 in chain not identified as
> Protein/RNA/DNA.
> >> Problem with chain definition, or missing terminal residues.
> >> This chain does not appear to contain a recognized chain molecule.
> >> If this is incorrect, you can edit residuetypes.dat to modify the
> >> behavior.
> >> 8 out of 8 lines of specbond.dat converted successfully
> >>
> >> --**-
> >> Program pdb2gmx, VERSION 4.5.4
> >> Source code file: resall.c, line: 581
> >>
> >> Fatal error:
> >> Residue 'ETY' not found in residue topology database
> >> For more information and tips for troubleshooting, please check the
> >> GROMACS
> >> website at http://www.gromacs.org/**Documentation/Errors<
> http://www.gromacs.org/Documentation/Errors>
> >>
> >>
> >> I am kind of baffled what could be wrong. I would really appreciate if
> you
> >> could comment on this. Is there a bug in pdb2gmx version 4.5.4 (which I
> >> doubt)?
> >> Do i have to edit esiduetypes.dat?
> >>
> >>
> > You've added the .rtp entry incorrectly.  pdb2gmx in version 4.5 and
> above
> > needs .rtp files to be organized in force field subdirectories (in your
> > case oplsaa.ff).  It will not recognize a file called "ffoplsa

Re: [gmx-users] Topology file-residue unknown

2013-04-05 Thread Juliette N.
Thank you Justin. I have now the topology. I have a quick question
regarding atom types. I used opls_143 and opls_144 for C and H of ethylene
respec in rtp. However, in VMD I dont see double bonds C=C. I though maybe
these are not the proper atomtypes for H2C=CH2.

rtp entry:

[ atoms ]
   C1opls_143-0.1201
   H11   opls_144 0.0601
   H12   opls_144 0.0601
   C2opls_143-0.1202
   H21   opls_144 0.0602
   H22   opls_144 0.0602

>From atomtypes.atp

opls_141   12.01100  ; alkene C (R2-C=)
 opls_142   12.01100  ; alkene C (RH-C=)
* opls_143   12.01100  ; alkene C (H2-C=)
 opls_1441.00800  ; alkene H (H-C=)*

Did I select the correct atomtypes?
Thanks many times!




On 4 April 2013 18:23, Justin Lemkul  wrote:

>
>
> On 4/4/13 6:17 PM, Juliette N. wrote:
>
>> Hi Justin,
>>
>> Thanks a lot for your message. I am petrified why pdb2gmx is not
>> recognizing the residue ETY. This the first residue added to ffoplsaa.rtp
>> but as you may see below it is not read from rtp.
>>
>> I am sure this residue is added to ffoplsaa.rtp which is existing in the
>> working directory:
>>
>> [ ETY ]
>>   [ atoms ]
>> C1opls_143-0.1201
>> H11   opls_144 0.0601
>> H12   opls_144 0.0601
>> C2opls_143-0.1202
>> H21   opls_144 0.0602
>> H22   opls_144 0.0602
>>
>> [ bonds ]
>> C1H11
>> C1H12
>> C1C2
>> C2H21
>> C2H22
>>;C2C1
>>
>> Then I issue:
>>
>> pdb2gmx -f Ethylene.gro -o GRO.gro -p Ethylene.top -ff oplsaa
>>   :-)  G  R  O  M  A  C  S  (-:
>>
>>GROwing Monsters And Cloning Shrimps
>>
>>  :-)  VERSION 4.5.4  (-:
>>
>> Using the Oplsaa force field in directory oplsaa.ff
>>
>> Opening force field file
>> /usr/local/gromacs/share/**gromacs/top/oplsaa.ff/**watermodels.dat
>>
>> Select the Water Model:
>>   1: TIP4P  TIP 4-point, recommended
>>   2: TIP3P  TIP 3-point
>>   3: TIP5P  TIP 5-point
>>   4: SPCsimple point charge
>>   5: SPC/E  extended simple point charge
>>   6: None
>> 6
>> Opening force field file
>> /usr/local/gromacs/share/**gromacs/top/oplsaa.ff/**aminoacids.r2b
>> Reading Ethylene.gro...
>> Read 'Giant Rising Ordinary Mutants for A Clerical Setup', 6 atoms
>> Analyzing pdb file
>> Splitting PDB chains based on TER records or changing chain id.
>> There are 1 chains and 0 blocks of water and 1 residues with 6 atoms
>>
>>chain  #res #atoms
>>1 ' ' 1  6
>>
>> No occupancies in Ethylene.gro
>> Opening force field file
>> /usr/local/gromacs/share/**gromacs/top/oplsaa.ff/**atomtypes.atp
>> Atomtype 1
>> Reading residue database... (oplsaa)
>> Opening force field file
>> /usr/local/gromacs/share/**gromacs/top/oplsaa.ff/**aminoacids.rtp
>> Residue 56
>> Sorting it all out...
>> Opening force field file
>> /usr/local/gromacs/share/**gromacs/top/oplsaa.ff/**aminoacids.hdb
>> Opening force field file
>> /usr/local/gromacs/share/**gromacs/top/oplsaa.ff/**aminoacids.n.tdb
>> Opening force field file
>> /usr/local/gromacs/share/**gromacs/top/oplsaa.ff/**aminoacids.c.tdb
>>
>> Back Off! I just backed up Ethylene.top to ./#Ethylene.top.1#
>> Processing chain 1 (6 atoms, 1 residues)
>> There are 0 donors and 0 acceptors
>> There are 0 hydrogen bonds
>> Warning: Starting residue ETY1 in chain not identified as Protein/RNA/DNA.
>> Problem with chain definition, or missing terminal residues.
>> This chain does not appear to contain a recognized chain molecule.
>> If this is incorrect, you can edit residuetypes.dat to modify the
>> behavior.
>> 8 out of 8 lines of specbond.dat converted successfully
>>
>> --**-
>> Program pdb2gmx, VERSION 4.5.4
>> Source code file: resall.c, line: 581
>>
>> Fatal error:
>> Residue 'ETY' not found in residue topology database
>> For more information and tips for troubleshooting, please check the
>> GROMACS
>> website at 
>> http://www.gromacs.org/**Documentation/Errors
>>
>>
>> I am kind of baffled what could be wrong. I would really appreciate if you
>> could comment on this. Is there a bug in pdb2gmx version 4.5.4 (which I
>> doubt)?
>> Do i have to edit esiduetypes.dat?
>>
>>
> You've added the .rtp entry incorrectly.  pdb2gmx in version 4.5 and above
> needs .rtp files to be organized in force field subdirectories (in your
> case oplsaa.ff).  It will not recognize a file called "ffoplsaa.rtp," which
> is the pre-4.5 syntax for force field files.
>
>
> -Justin
>
> --
> ==**==
>
> Justin A. Lemkul, Ph.D.
> Research Scientist
> Department of Biochemistry
> Virginia Tech
> Blacksburg, VA
> jalemkul[at]vt.edu | (540) 231-9080
> http://www.bevanlab.biochem.**vt.edu/Pages/Personal/justin
>
> ==**==
> --

Re: [gmx-users] Topology file-residue unknown

2013-04-04 Thread Justin Lemkul



On 4/4/13 6:17 PM, Juliette N. wrote:

Hi Justin,

Thanks a lot for your message. I am petrified why pdb2gmx is not
recognizing the residue ETY. This the first residue added to ffoplsaa.rtp
but as you may see below it is not read from rtp.

I am sure this residue is added to ffoplsaa.rtp which is existing in the
working directory:

[ ETY ]
  [ atoms ]
C1opls_143-0.1201
H11   opls_144 0.0601
H12   opls_144 0.0601
C2opls_143-0.1202
H21   opls_144 0.0602
H22   opls_144 0.0602

[ bonds ]
C1H11
C1H12
C1C2
C2H21
C2H22
   ;C2C1

Then I issue:

pdb2gmx -f Ethylene.gro -o GRO.gro -p Ethylene.top -ff oplsaa
  :-)  G  R  O  M  A  C  S  (-:

   GROwing Monsters And Cloning Shrimps

 :-)  VERSION 4.5.4  (-:

Using the Oplsaa force field in directory oplsaa.ff

Opening force field file
/usr/local/gromacs/share/gromacs/top/oplsaa.ff/watermodels.dat

Select the Water Model:
  1: TIP4P  TIP 4-point, recommended
  2: TIP3P  TIP 3-point
  3: TIP5P  TIP 5-point
  4: SPCsimple point charge
  5: SPC/E  extended simple point charge
  6: None
6
Opening force field file
/usr/local/gromacs/share/gromacs/top/oplsaa.ff/aminoacids.r2b
Reading Ethylene.gro...
Read 'Giant Rising Ordinary Mutants for A Clerical Setup', 6 atoms
Analyzing pdb file
Splitting PDB chains based on TER records or changing chain id.
There are 1 chains and 0 blocks of water and 1 residues with 6 atoms

   chain  #res #atoms
   1 ' ' 1  6

No occupancies in Ethylene.gro
Opening force field file
/usr/local/gromacs/share/gromacs/top/oplsaa.ff/atomtypes.atp
Atomtype 1
Reading residue database... (oplsaa)
Opening force field file
/usr/local/gromacs/share/gromacs/top/oplsaa.ff/aminoacids.rtp
Residue 56
Sorting it all out...
Opening force field file
/usr/local/gromacs/share/gromacs/top/oplsaa.ff/aminoacids.hdb
Opening force field file
/usr/local/gromacs/share/gromacs/top/oplsaa.ff/aminoacids.n.tdb
Opening force field file
/usr/local/gromacs/share/gromacs/top/oplsaa.ff/aminoacids.c.tdb

Back Off! I just backed up Ethylene.top to ./#Ethylene.top.1#
Processing chain 1 (6 atoms, 1 residues)
There are 0 donors and 0 acceptors
There are 0 hydrogen bonds
Warning: Starting residue ETY1 in chain not identified as Protein/RNA/DNA.
Problem with chain definition, or missing terminal residues.
This chain does not appear to contain a recognized chain molecule.
If this is incorrect, you can edit residuetypes.dat to modify the behavior.
8 out of 8 lines of specbond.dat converted successfully

---
Program pdb2gmx, VERSION 4.5.4
Source code file: resall.c, line: 581

Fatal error:
Residue 'ETY' not found in residue topology database
For more information and tips for troubleshooting, please check the GROMACS
website at http://www.gromacs.org/Documentation/Errors


I am kind of baffled what could be wrong. I would really appreciate if you
could comment on this. Is there a bug in pdb2gmx version 4.5.4 (which I
doubt)?
Do i have to edit esiduetypes.dat?



You've added the .rtp entry incorrectly.  pdb2gmx in version 4.5 and above needs 
.rtp files to be organized in force field subdirectories (in your case 
oplsaa.ff).  It will not recognize a file called "ffoplsaa.rtp," which is the 
pre-4.5 syntax for force field files.


-Justin

--


Justin A. Lemkul, Ph.D.
Research Scientist
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] Topology file-residue unknown

2013-04-04 Thread Juliette N.
Hi Justin,

Thanks a lot for your message. I am petrified why pdb2gmx is not
recognizing the residue ETY. This the first residue added to ffoplsaa.rtp
but as you may see below it is not read from rtp.

I am sure this residue is added to ffoplsaa.rtp which is existing in the
working directory:

[ ETY ]
 [ atoms ]
   C1opls_143-0.1201
   H11   opls_144 0.0601
   H12   opls_144 0.0601
   C2opls_143-0.1202
   H21   opls_144 0.0602
   H22   opls_144 0.0602

[ bonds ]
   C1H11
   C1H12
   C1C2
   C2H21
   C2H22
  ;C2C1

Then I issue:

pdb2gmx -f Ethylene.gro -o GRO.gro -p Ethylene.top -ff oplsaa
 :-)  G  R  O  M  A  C  S  (-:

  GROwing Monsters And Cloning Shrimps

:-)  VERSION 4.5.4  (-:

Using the Oplsaa force field in directory oplsaa.ff

Opening force field file
/usr/local/gromacs/share/gromacs/top/oplsaa.ff/watermodels.dat

Select the Water Model:
 1: TIP4P  TIP 4-point, recommended
 2: TIP3P  TIP 3-point
 3: TIP5P  TIP 5-point
 4: SPCsimple point charge
 5: SPC/E  extended simple point charge
 6: None
6
Opening force field file
/usr/local/gromacs/share/gromacs/top/oplsaa.ff/aminoacids.r2b
Reading Ethylene.gro...
Read 'Giant Rising Ordinary Mutants for A Clerical Setup', 6 atoms
Analyzing pdb file
Splitting PDB chains based on TER records or changing chain id.
There are 1 chains and 0 blocks of water and 1 residues with 6 atoms

  chain  #res #atoms
  1 ' ' 1  6

No occupancies in Ethylene.gro
Opening force field file
/usr/local/gromacs/share/gromacs/top/oplsaa.ff/atomtypes.atp
Atomtype 1
Reading residue database... (oplsaa)
Opening force field file
/usr/local/gromacs/share/gromacs/top/oplsaa.ff/aminoacids.rtp
Residue 56
Sorting it all out...
Opening force field file
/usr/local/gromacs/share/gromacs/top/oplsaa.ff/aminoacids.hdb
Opening force field file
/usr/local/gromacs/share/gromacs/top/oplsaa.ff/aminoacids.n.tdb
Opening force field file
/usr/local/gromacs/share/gromacs/top/oplsaa.ff/aminoacids.c.tdb

Back Off! I just backed up Ethylene.top to ./#Ethylene.top.1#
Processing chain 1 (6 atoms, 1 residues)
There are 0 donors and 0 acceptors
There are 0 hydrogen bonds
Warning: Starting residue ETY1 in chain not identified as Protein/RNA/DNA.
Problem with chain definition, or missing terminal residues.
This chain does not appear to contain a recognized chain molecule.
If this is incorrect, you can edit residuetypes.dat to modify the behavior.
8 out of 8 lines of specbond.dat converted successfully

---
Program pdb2gmx, VERSION 4.5.4
Source code file: resall.c, line: 581

Fatal error:
Residue 'ETY' not found in residue topology database
For more information and tips for troubleshooting, please check the GROMACS
website at http://www.gromacs.org/Documentation/Errors


I am kind of baffled what could be wrong. I would really appreciate if you
could comment on this. Is there a bug in pdb2gmx version 4.5.4 (which I
doubt)?
Do i have to edit esiduetypes.dat?

Thank you many times for your precious time,




On 4 April 2013 17:21, Justin Lemkul  wrote:

>
>
> On 4/4/13 11:38 AM, Juliette N. wrote:
>
>> Dear all,
>>
>> I am having trouble creating topology file for simple molecule ethylene.
>> Here is the steps. Below is the pdb file:
>>
>> Ethylene.pdb:
>>
>>
>> ATOM  1  C1  ETY 1  0.672  -0.000   0.000  1.00  0.00
>> C
>> ATOM  2  C2  ETY 1 -0.672  -0.000   0.000  1.00  0.00
>> C
>> ATOM  3  H11 ETY 1  1.238  -0.928   0.000  1.00  0.00
>> H
>> ATOM  4  H12 ETY 1  1.238   0.928   0.000  1.00  0.00
>> H
>> ATOM  5  H21 ETY 1 -1.238  -0.928   0.000  1.00  0.00
>> H
>> ATOM  6  H22 ETY 1 -1.238   0.928   0.000  1.00  0.00
>> H
>> END
>>
>> Then:
>> editconf -f Ethylene.pdb -o Ethylene.gro
>>
>> Ethylene.gro:
>>
>>  6
>>  1ETY C11   0.067  -0.000   0.000
>>  1ETY C22  -0.067  -0.000   0.000
>>  1ETYH113   0.124  -0.093   0.000
>>  1ETYH124   0.124   0.093   0.000
>>  1ETYH215  -0.124  -0.093   0.000
>>  1ETYH226  -0.124   0.093   0.000
>> 0.0   0.0   0.0
>>
>>
>> Then I add the following residue for ethylene to the rtp file
>> ffoplsaa.rtp:
>>
>> [ ETY ]
>>   [ atoms ]
>> C1opls_143-0.1201
>> H11   opls_144 0.0601
>> H12   opls_144 0.0601
>> C2opls_143-0.1202
>> H21   opls_144 0.0602
>> H22   opls_144 0.0602
>>
>> [ bonds ]
>> C1H11
>> C1H12
>> C1C2
>> C2H21
>> C2H22
>> C2C1
>>
>>
> "C1 C2" and "C2 C1" are redundant.  You don't need both.
>
>
>  and then
>>
>> ; pdb2gmx -f Ethylene.pdb -o GRO.gro -p Ethylene.top -ff oplsaa
>>
>> Q1_ i dont realize what is the difference between the already created
>> Ethylene.gro and the o

Re: [gmx-users] Topology file-residue unknown

2013-04-04 Thread Justin Lemkul



On 4/4/13 11:38 AM, Juliette N. wrote:

Dear all,

I am having trouble creating topology file for simple molecule ethylene.
Here is the steps. Below is the pdb file:

Ethylene.pdb:


ATOM  1  C1  ETY 1  0.672  -0.000   0.000  1.00  0.00
C
ATOM  2  C2  ETY 1 -0.672  -0.000   0.000  1.00  0.00
C
ATOM  3  H11 ETY 1  1.238  -0.928   0.000  1.00  0.00
H
ATOM  4  H12 ETY 1  1.238   0.928   0.000  1.00  0.00
H
ATOM  5  H21 ETY 1 -1.238  -0.928   0.000  1.00  0.00
H
ATOM  6  H22 ETY 1 -1.238   0.928   0.000  1.00  0.00
H
END

Then:
editconf -f Ethylene.pdb -o Ethylene.gro

Ethylene.gro:

 6
 1ETY C11   0.067  -0.000   0.000
 1ETY C22  -0.067  -0.000   0.000
 1ETYH113   0.124  -0.093   0.000
 1ETYH124   0.124   0.093   0.000
 1ETYH215  -0.124  -0.093   0.000
 1ETYH226  -0.124   0.093   0.000
0.0   0.0   0.0


Then I add the following residue for ethylene to the rtp file
ffoplsaa.rtp:

[ ETY ]
  [ atoms ]
C1opls_143-0.1201
H11   opls_144 0.0601
H12   opls_144 0.0601
C2opls_143-0.1202
H21   opls_144 0.0602
H22   opls_144 0.0602

[ bonds ]
C1H11
C1H12
C1C2
C2H21
C2H22
C2C1



"C1 C2" and "C2 C1" are redundant.  You don't need both.


and then

; pdb2gmx -f Ethylene.pdb -o GRO.gro -p Ethylene.top -ff oplsaa

Q1_ i dont realize what is the difference between the already created
Ethylene.gro and the one in the pdb2gmx(GRO.gro). I am not getting the
latter and was wondering what would be the diffrence?



pdb2gmx produces an output coordinate file as a byproduct of its massively 
complicated fix-build-process procedure.


If you're not getting the output file(s), then pdb2gmx is failing and should 
produce obvious error messages.



Q2) Below is the  top file which is not generated. Once I Issue pdb2gmx, I
enter "6" because pdb2gmx output says:

Opening force field file
/usr/local/gromacs/share/gromacs/top/oplsaa.ff/watermo$
  1: TIP4P  TIP 4-point, recommended
  2: TIP3P  TIP 3-point
  3: TIP5P  TIP 5-point
  4: SPCsimple point charge
  5: SPC/E  extended simple point charge
  6: None
Opening force field file
/usr/local/gromacs/share/gromacs/top/oplsaa.ff/aminoac$
No occupancies in Ethylene.gro
Opening force field file
/usr/local/gromacs/share/gromacs/top/oplsaa.ff/atomtyp$

And finally I get a blank top file::

;File 'Ethylene.top' was generated
;By user: onbekend (0)
;On host: onbekend
;At date:
;
;This is a standalone topology file
;
;It was generated using program:
;pdb2gmx - VERSION 4.5.4
;
;Command line was:
;pdb2gmx -f Ethylene.pdb -o GRO.gro -p Ethylene.top -ff oplsaa
;
;Force field was read from the standard Gromacs share directory.
;

; Include forcefield parameters
#include "oplsaa.ff/forcefield.itp"


-
output says:

Select the Water Model:
Reading Ethylene.gro...
Read 'GROningen Mixture of Alchemy and Childrens' Stories', 6 atoms
Analyzing pdb file
Splitting PDB chains based on TER records or changing chain id.
There are 1 chains and 0 blocks of water and 1 residues with 6 atoms

   chain  #res #atoms
   1 ' ' 1  6

Reading residue database... (oplsaa)
Processing chain 1 (6 atoms, 1 residues)
Warning: Starting residue ETY1 in chain not identified as Protein/RNA/DNA.
Problem with chain definition, or missing terminal residues.
This chain does not appear to contain a recognized chain molecule.
If this is incorrect, you can edit residuetypes.dat to modify the behavior.



There should be more output after this line.  Any reason for pdb2gmx to fail 
would be printed here.





Is there something wrong with my rtp? Could I trouble you to let me know
what wrong is?



I can see nothing wrong.  I was able to produce a sensible topology using all of 
your input files.


-Justin

--


Justin A. Lemkul, Ph.D.
Research Scientist
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] Topology file-residue unknown

2013-04-04 Thread Juliette N.
Dear all,

I am having trouble creating topology file for simple molecule ethylene.
Here is the steps. Below is the pdb file:

Ethylene.pdb:


ATOM  1  C1  ETY 1  0.672  -0.000   0.000  1.00  0.00
C
ATOM  2  C2  ETY 1 -0.672  -0.000   0.000  1.00  0.00
C
ATOM  3  H11 ETY 1  1.238  -0.928   0.000  1.00  0.00
H
ATOM  4  H12 ETY 1  1.238   0.928   0.000  1.00  0.00
H
ATOM  5  H21 ETY 1 -1.238  -0.928   0.000  1.00  0.00
H
ATOM  6  H22 ETY 1 -1.238   0.928   0.000  1.00  0.00
H
END

Then:
editconf -f Ethylene.pdb -o Ethylene.gro

Ethylene.gro:

6
1ETY C11   0.067  -0.000   0.000
1ETY C22  -0.067  -0.000   0.000
1ETYH113   0.124  -0.093   0.000
1ETYH124   0.124   0.093   0.000
1ETYH215  -0.124  -0.093   0.000
1ETYH226  -0.124   0.093   0.000
   0.0   0.0   0.0


Then I add the following residue for ethylene to the rtp file
ffoplsaa.rtp:

[ ETY ]
 [ atoms ]
   C1opls_143-0.1201
   H11   opls_144 0.0601
   H12   opls_144 0.0601
   C2opls_143-0.1202
   H21   opls_144 0.0602
   H22   opls_144 0.0602

[ bonds ]
   C1H11
   C1H12
   C1C2
   C2H21
   C2H22
   C2C1

and then

; pdb2gmx -f Ethylene.pdb -o GRO.gro -p Ethylene.top -ff oplsaa

Q1_ i dont realize what is the difference between the already created
Ethylene.gro and the one in the pdb2gmx(GRO.gro). I am not getting the
latter and was wondering what would be the diffrence?

Q2) Below is the  top file which is not generated. Once I Issue pdb2gmx, I
enter "6" because pdb2gmx output says:

Opening force field file
/usr/local/gromacs/share/gromacs/top/oplsaa.ff/watermo$
 1: TIP4P  TIP 4-point, recommended
 2: TIP3P  TIP 3-point
 3: TIP5P  TIP 5-point
 4: SPCsimple point charge
 5: SPC/E  extended simple point charge
 6: None
Opening force field file
/usr/local/gromacs/share/gromacs/top/oplsaa.ff/aminoac$
No occupancies in Ethylene.gro
Opening force field file
/usr/local/gromacs/share/gromacs/top/oplsaa.ff/atomtyp$

And finally I get a blank top file::

;File 'Ethylene.top' was generated
;By user: onbekend (0)
;On host: onbekend
;At date:
;
;This is a standalone topology file
;
;It was generated using program:
;pdb2gmx - VERSION 4.5.4
;
;Command line was:
;pdb2gmx -f Ethylene.pdb -o GRO.gro -p Ethylene.top -ff oplsaa
;
;Force field was read from the standard Gromacs share directory.
;

; Include forcefield parameters
#include "oplsaa.ff/forcefield.itp"


-
output says:

Select the Water Model:
Reading Ethylene.gro...
Read 'GROningen Mixture of Alchemy and Childrens' Stories', 6 atoms
Analyzing pdb file
Splitting PDB chains based on TER records or changing chain id.
There are 1 chains and 0 blocks of water and 1 residues with 6 atoms

  chain  #res #atoms
  1 ' ' 1  6

Reading residue database... (oplsaa)
Processing chain 1 (6 atoms, 1 residues)
Warning: Starting residue ETY1 in chain not identified as Protein/RNA/DNA.
Problem with chain definition, or missing terminal residues.
This chain does not appear to contain a recognized chain molecule.
If this is incorrect, you can edit residuetypes.dat to modify the behavior.



Is there something wrong with my rtp? Could I trouble you to let me know
what wrong is?

-- 
Thanks,
J. N.
-- 
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] topology

2013-03-18 Thread Justin Lemkul



On 3/18/13 10:49 PM, shahid nayeem wrote:

I am using Gromacs-4.5.4 and this does not have -cys option in pdb2gmx.
shahid



My apologies.  For some reason I always think there is a -cys option, when in 
fact it is -ss.  In any case, if you're trying to treat some strange form of 
cysteine, some force fields have special names for different forms (CYM for 
thiolate, etc) so naming and proper force field choice may be the only viable 
mechanism here.


-Justin


On Tue, Mar 19, 2013 at 2:06 AM, Justin Lemkul  wrote:



On 3/18/13 10:39 AM, shahid nayeem wrote:


Hi All
   How can I get a topology file using pdb2gmx from a single chain
polypeptide with one of the CYS, SH in the form of SG as if it is
involved in disulfide linkage, while the other chain with which I
expect it to form disulfide link is not in the input pdb file. or can
I use pdb2gmx command and get the separate .top and .gro file for each
chain using index file.



You may be able to use the -cys option to fool pdb2gmx, but I don't know if
that will cause problems.  Having half a disulfide present is odd, though
some force fields support thiolate or neutral (i.e. radical) forms of
cysteine.  You'll have to look into the force field files to see what's
available.

-Justin

--


Justin A. Lemkul, Ph.D.
Research Scientist
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


--


Justin A. Lemkul, Ph.D.
Research Scientist
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] topology

2013-03-18 Thread shahid nayeem
I am using Gromacs-4.5.4 and this does not have -cys option in pdb2gmx.
shahid

On Tue, Mar 19, 2013 at 2:06 AM, Justin Lemkul  wrote:
>
>
> On 3/18/13 10:39 AM, shahid nayeem wrote:
>>
>> Hi All
>>   How can I get a topology file using pdb2gmx from a single chain
>> polypeptide with one of the CYS, SH in the form of SG as if it is
>> involved in disulfide linkage, while the other chain with which I
>> expect it to form disulfide link is not in the input pdb file. or can
>> I use pdb2gmx command and get the separate .top and .gro file for each
>> chain using index file.
>>
>
> You may be able to use the -cys option to fool pdb2gmx, but I don't know if
> that will cause problems.  Having half a disulfide present is odd, though
> some force fields support thiolate or neutral (i.e. radical) forms of
> cysteine.  You'll have to look into the force field files to see what's
> available.
>
> -Justin
>
> --
> 
>
> Justin A. Lemkul, Ph.D.
> Research Scientist
> 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 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] topology

2013-03-18 Thread Justin Lemkul



On 3/18/13 10:39 AM, shahid nayeem wrote:

Hi All
  How can I get a topology file using pdb2gmx from a single chain
polypeptide with one of the CYS, SH in the form of SG as if it is
involved in disulfide linkage, while the other chain with which I
expect it to form disulfide link is not in the input pdb file. or can
I use pdb2gmx command and get the separate .top and .gro file for each
chain using index file.



You may be able to use the -cys option to fool pdb2gmx, but I don't know if that 
will cause problems.  Having half a disulfide present is odd, though some force 
fields support thiolate or neutral (i.e. radical) forms of cysteine.  You'll 
have to look into the force field files to see what's available.


-Justin

--


Justin A. Lemkul, Ph.D.
Research Scientist
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] topology

2013-03-18 Thread shahid nayeem
Hi All
 How can I get a topology file using pdb2gmx from a single chain
polypeptide with one of the CYS, SH in the form of SG as if it is
involved in disulfide linkage, while the other chain with which I
expect it to form disulfide link is not in the input pdb file. or can
I use pdb2gmx command and get the separate .top and .gro file for each
chain using index file.

shahid
-- 
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] Topology of Graphite Oxide

2013-03-07 Thread Justin Lemkul



On 3/7/13 3:36 AM, zhhxu wrote:

Hi  users:

I want to run a simulation of adsorption of protein on graphite oxide(GO) which 
is oxidized at random and
is near 80% coverage by hydroxyl and epoxyl group.
My first question is how to generate a topology filefor GO. Previous to search 
solution to this problem in
this forum, I had tried to make a topology for GO use g_x2top with command:
g_x2top -f input.pdb -nopbc -r outputname -o outputname.
But g_x2top did not work fine(fine for prestine graphite) and gave error messages like 
"Can not find forcefield for atom .. "



The .n2t files provided with the (few) force fields in Gromacs are rather 
inflexible, so you will probably need to modify them to suit your needs.



My second also is about definition of improper dihedral in opls-aa. Whether it 
is validate to assign a improper
dihedral for a central sp2 carbon atom with three other bonded atoms no matter 
of what is type of  other
three atom?  I have gotten (mentioned above) a topology file of pristine 
graphite for opls-aa with
  g_x2top(the command is similar to above),  but the topology file did not 
include any improper dihedral.



g_x2top doesn't write impropers, so you have to manually add them if needed. 
Assignment of proper parameters should be done based on analogous functional 
groups present in the force field or by a suitable derivation protocol.


-Justin

--


Justin A. Lemkul, Ph.D.
Research Scientist
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] Topology of Graphite Oxide

2013-03-07 Thread zhhxu
Hi  users:

I want to run a simulation of adsorption of protein on graphite oxide(GO) which 
is oxidized at random and
is near 80% coverage by hydroxyl and epoxyl group.
My first question is how to generate a topology filefor GO. Previous to search 
solution to this problem in
this forum, I had tried to make a topology for GO use g_x2top with command:
g_x2top -f input.pdb -nopbc -r outputname -o outputname.
But g_x2top did not work fine(fine for prestine graphite) and gave error 
messages like "Can not find forcefield for atom .. "

My second also is about definition of improper dihedral in opls-aa. Whether it 
is validate to assign a improper
dihedral for a central sp2 carbon atom with three other bonded atoms no matter 
of what is type of  other
three atom?  I have gotten (mentioned above) a topology file of pristine 
graphite for opls-aa with
 g_x2top(the command is similar to above),  but the topology file did not 
include any improper dihedral.

  I will appreciated you for giving any advice.

   
sincerely
zhhxu PhD student
Institute of Process Engineering, Chinese Academy of Sciences,
Department of Biochemical Engineering.






-- 
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] Topology and mdp file for coupling vdw and coulombic interactions for running FE calculations using BAR

2012-09-17 Thread Justin Lemkul



On 9/16/12 10:45 PM, Sonia Milena Aguilera Segura wrote:

Hi,

I'm preparing my mdp and topology files for running free energy calculations 
using BAR method. I´m using Justin Lemkul's tutorial as a reference (You can 
find it here 
http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin/gmx-tutorials/free_energy/08_advanced.html).
 According to this tutorial, the simultaneous coupling of both Coulombic and 
van der Waals terms leads to instability. So, it is usefull to prepare the mdp 
files as follows (off course including all other parameters):


van der Waals coupling:

  sc-alpha  = 0.5 ; use soft-core for LJ (de)coupling
  sc-sigma  = 0.3
  sc-power  = 1
  couple-moltype= LIG
  couple-intramol   = no
  couple-lambda0= none; non-interacting dummy in state A
  couple-lambda1= vdw ; only vdW terms on in state B


Coulombic coupling:

  sc-alpha  = 0   ; soft-core during (dis)charging can be 
unstable!
  sc-sigma  = 0
  couple-moltype= LIG
  couple-intramol   = no
  couple-lambda0= vdw ; only vdW terms in state A (the previous 
state B is now A)
  couple-lambda1= vdw-q   ; all nonbonded interactions are on in state B

However, I don't understand how can this leads to a fully interacting molecule 
if all the charges in the topology file have been set to zero. Does it mean 
that for the second calculation (coulombic coupling) I have to use the original 
topology file with all charges? Or, should I use the same topology with zero 
charges?



You do not need to make any changes to the topology.  The couple-lambda* 
settings dictate which parameters are involved in the calculations.  In the 
second examples, in the lambda=0 state, only van der Waals interactions are 
calculated, but when lambda=1, both van der Waals and Coulombic interactions 
(vdw-q) are calculated.


-Justin

--


Justin A. Lemkul, Ph.D.
Research Scientist
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] Topology and mdp file for coupling vdw and coulombic interactions for running FE calculations using BAR

2012-09-16 Thread Sonia Milena Aguilera Segura
Hi,

I'm preparing my mdp and topology files for running free energy calculations 
using BAR method. I´m using Justin Lemkul's tutorial as a reference (You can 
find it here 
http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin/gmx-tutorials/free_energy/08_advanced.html).
 According to this tutorial, the simultaneous coupling of both Coulombic and 
van der Waals terms leads to instability. So, it is usefull to prepare the mdp 
files as follows (off course including all other parameters):


van der Waals coupling:

 sc-alpha  = 0.5 ; use soft-core for LJ (de)coupling
 sc-sigma  = 0.3
 sc-power  = 1
 couple-moltype= LIG
 couple-intramol   = no
 couple-lambda0= none; non-interacting dummy in state A
 couple-lambda1= vdw ; only vdW terms on in state B


Coulombic coupling:

 sc-alpha  = 0   ; soft-core during (dis)charging can be 
unstable!
 sc-sigma  = 0
 couple-moltype= LIG
 couple-intramol   = no
 couple-lambda0= vdw ; only vdW terms in state A (the previous 
state B is now A)
 couple-lambda1= vdw-q   ; all nonbonded interactions are on in state B

However, I don't understand how can this leads to a fully interacting molecule 
if all the charges in the topology file have been set to zero. Does it mean 
that for the second calculation (coulombic coupling) I have to use the original 
topology file with all charges? Or, should I use the same topology with zero 
charges?

Thanks in advance,

Sonia Aguilera
Graduate student-Chemical Engineering Department
Universidad de los Andes
Colombia


--
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] Topology generated of a dimer protein

2012-09-09 Thread Justin Lemkul



On 9/9/12 8:42 AM, Shima Arasteh wrote:

Hi all,

I am simulating a system of protein and water. My protein is composed of 2 
monomers, each monomer has its own C- and N-terminus, so there are 2 
N-terminals and 2 C-terminals.
I have a pdb of two monomers as dimer.pdb.

I tried to generate topology of dimer.pdb. But 2 monomers are identified as one 
fragment which has one -N-terminus and one C-terminus! So the generated 
topology file and the total charge would be incorrect.
How can I solve it? Any suggestions please?



The chains need to be either separated by a TER entry or have different chain 
identifiers.  pdb2gmx uses the -chainsep option to control the criteria for a 
new chain, and by default it uses either of these.


-Justin

--


Justin A. Lemkul, Ph.D.
Research Scientist
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] Topology generated of a dimer protein

2012-09-09 Thread Shima Arasteh
Hi all,

I am simulating a system of protein and water. My protein is composed of 2 
monomers, each monomer has its own C- and N-terminus, so there are 2 
N-terminals and 2 C-terminals.
I have a pdb of two monomers as dimer.pdb. 

I tried to generate topology of dimer.pdb. But 2 monomers are identified as one 
fragment which has one -N-terminus and one C-terminus! So the generated 
topology file and the total charge would be incorrect.
How can I solve it? Any suggestions please?

Your suggestions would be appreciated.

Thanks in advance.



Sincerely,
Shima
-- 
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] Topology- make bonds with certain cut-off distance

2012-08-26 Thread Mark Abraham

On 25/08/2012 11:58 PM, mohan maruthi sena wrote:

Hi ,
  Thanks for your reply.  I actually want to build an elastic
network model (ENM) for Protein containing 372 residues. According to
literature, ENM model only considers C-alpha atoms of protein and
discards rest of the atoms.The potential only considers bonded term
only. It does not
contain angle,dihedral or non-bonded term.  To build an ENM model the
specified atom makes bond with all other atoms with certain cutoff
distance. I consider the distance between two c-alpha atoms is 0.38nm
and force constant of 69.1 kj/(nm)2.I have written a code which gives
the number of  atoms around a particular atom with in certain cut-off
distance, With this i will know what are the atoms that are needed to
connect a particular atom.I have generated .gro and .top files using
pdb2gmx command .  After generating topology file , I replace bonds
section with the  connectivity information i get from code. The next
steps are equilibration and production. I am not sure whether it takes
connectivity information from topology or not. If it is not then how
can i mention it forms bonds with in certain distance .


You can only do this by treating your bonds as what GROMACS calls 
non-bonded interactions, since only those are detected dynamically based 
on distance. So setting up a neighbour-search cut-off and tabulated 
interaction with the same cut-off could implement the above...



  Finally , if i
run the simulation with all the above mentioned steps  and load the
trajectory in vmd , the protien gets contracted after first step.


... but since this model has only an attractive inter-particle 
interaction and  no way to exclude atoms from each other's volume, it is 
garbage. Rather than inventing stuff randomly, find out what has already 
been done in this area and how this (very basic) issue is solved. As a 
general rule, one should build on existing ideas until one has the 
judgement and experience to try something radical and new.


Mark



Thanks,
Mohan

On Sat, Aug 25, 2012 at 3:43 PM, Erik Marklund  wrote:

Hi,

In case you want the "bonds" to be dynamic I think the only way is to use 
tabulated interactions. A bit tedious to set up, but doable. Note that this will NOT 
emulate angles or dihedrals, just the stretching term. I can't understand why you would 
want bonds between all atoms within a certain distance, however, and I hope you know what 
you're doing.

Erik

25 aug 2012 kl. 06.54 skrev mohan maruthi sena:


Hi all,
   I want build a toplogy for a protein , for which , each and
every atom  has to  make bonds with other atoms with in certain
specified cut-off distance. How can i do this ?

Please suggest me a way,


Thanks in advance,
Mohan
--
gmx-users mailing listgmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
* Only plain text messages are allowed!
* 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
Dept. of Cell and Molecular Biology, Uppsala University.
Husargatan 3, Box 596,75124 Uppsala, Sweden
phone:+46 18 471 6688fax: +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
* Only plain text messages are allowed!
* 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
* Only plain text messages are allowed!
* 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] Topology- make bonds with certain cut-off distance

2012-08-25 Thread mohan maruthi sena
Hi ,
 Thanks for your reply.  I actually want to build an elastic
network model (ENM) for Protein containing 372 residues. According to
literature, ENM model only considers C-alpha atoms of protein and
discards rest of the atoms.The potential only considers bonded term
only. It does not
contain angle,dihedral or non-bonded term.  To build an ENM model the
specified atom makes bond with all other atoms with certain cutoff
distance. I consider the distance between two c-alpha atoms is 0.38nm
and force constant of 69.1 kj/(nm)2.I have written a code which gives
the number of  atoms around a particular atom with in certain cut-off
distance, With this i will know what are the atoms that are needed to
connect a particular atom.I have generated .gro and .top files using
pdb2gmx command .  After generating topology file , I replace bonds
section with the  connectivity information i get from code. The next
steps are equilibration and production. I am not sure whether it takes
connectivity information from topology or not. If it is not then how
can i mention it forms bonds with in certain distance . Finally , if i
run the simulation with all the above mentioned steps  and load the
trajectory in vmd , the protien gets contracted after first step.

Thanks,
Mohan

On Sat, Aug 25, 2012 at 3:43 PM, Erik Marklund  wrote:
> Hi,
>
> In case you want the "bonds" to be dynamic I think the only way is to use 
> tabulated interactions. A bit tedious to set up, but doable. Note that this 
> will NOT emulate angles or dihedrals, just the stretching term. I can't 
> understand why you would want bonds between all atoms within a certain 
> distance, however, and I hope you know what you're doing.
>
> Erik
>
> 25 aug 2012 kl. 06.54 skrev mohan maruthi sena:
>
>> Hi all,
>>   I want build a toplogy for a protein , for which , each and
>> every atom  has to  make bonds with other atoms with in certain
>> specified cut-off distance. How can i do this ?
>>
>> Please suggest me a way,
>>
>>
>> Thanks in advance,
>> Mohan
>> --
>> gmx-users mailing listgmx-users@gromacs.org
>> http://lists.gromacs.org/mailman/listinfo/gmx-users
>> * Only plain text messages are allowed!
>> * 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
> Dept. of Cell and Molecular Biology, Uppsala University.
> Husargatan 3, Box 596,75124 Uppsala, Sweden
> phone:+46 18 471 6688fax: +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
> * Only plain text messages are allowed!
> * 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
* Only plain text messages are allowed!
* 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] Topology- make bonds with certain cut-off distance

2012-08-25 Thread Erik Marklund
Hi,

In case you want the "bonds" to be dynamic I think the only way is to use 
tabulated interactions. A bit tedious to set up, but doable. Note that this 
will NOT emulate angles or dihedrals, just the stretching term. I can't 
understand why you would want bonds between all atoms within a certain 
distance, however, and I hope you know what you're doing.

Erik
 
25 aug 2012 kl. 06.54 skrev mohan maruthi sena:

> Hi all,
>   I want build a toplogy for a protein , for which , each and
> every atom  has to  make bonds with other atoms with in certain
> specified cut-off distance. How can i do this ?
> 
> Please suggest me a way,
> 
> 
> Thanks in advance,
> Mohan
> -- 
> gmx-users mailing listgmx-users@gromacs.org
> http://lists.gromacs.org/mailman/listinfo/gmx-users
> * Only plain text messages are allowed!
> * 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
Dept. of Cell and Molecular Biology, Uppsala University.
Husargatan 3, Box 596,75124 Uppsala, Sweden
phone:+46 18 471 6688fax: +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
* Only plain text messages are allowed!
* 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] Topology- make bonds with certain cut-off distance

2012-08-25 Thread mohan maruthi sena
Hi all,
   I want build a toplogy for a protein , for which , each and
every atom  has to  make bonds with other atoms with in certain
specified cut-off distance. How can i do this ?

Please suggest me a way,


Thanks in advance,
Mohan
-- 
gmx-users mailing listgmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
* Only plain text messages are allowed!
* 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] Topology file

2012-08-15 Thread Ankita naithani
Hi,

I noticed that in my topology file, there is no inclusion of position
restraint file for my protein. For instance, my topology file looks
like this:

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

; Include chain topologies
#include "topol_Protein_chain_D.itp"
#include "topol_Protein_chain_E.itp"
#include "topol_Protein_chain_F.itp"
#include "topol_Protein_chain_G.itp"

;Include ligand topology
#include "FDP_Dp.itp"

; Include water topology
#include "gromos53a6.ff/spc.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 "gromos53a6.ff/ions.itp"


But, after running pdb2gmx, I do get 4 posre files for all the
individual chains. Do I need to add in them manually? Also, I am
unable to understand as to why did it not get included in the first
instance itself? I am using gromos53a6 force field.


Best Wishes,

-- 
Ankita Naithani
-- 
gmx-users mailing listgmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
* Only plain text messages are allowed!
* 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] Topology file for simulation of inorganic compounds

2012-08-13 Thread David van der Spoel

On 2012-08-13 21:39, Justin Lemkul wrote:



On 8/13/12 1:49 PM, Deepak Ojha wrote:

Dear All

I want to do simulation of some ions like azide ion,thiocynate ion and
other inorganic compounds in liquid water.
I am not sure how to generate topology file for the same.How can I go
about doing the same.



Likely you will have to do it manually, which means you will have to
understand Chapter 5 of the manual and the proper parameterization
protocol for your molecules under the chosen force field.

Some web servers exist for such purposes, depending on the force field
you want to use.  Otherwise consult the following page and be prepared
for a significant time investment:

http://www.gromacs.org/Documentation/How-tos/Parameterization

-Justin


Hm, that page is slightly outdated. You want to google for
antechamber / gaff /amber
respectively
charmm / cgenff

--
David van der Spoel, Ph.D., Professor of Biology
Dept. of Cell & Molec. Biol., Uppsala University.
Box 596, 75124 Uppsala, Sweden. Phone:  +46184714205.
sp...@xray.bmc.uu.sehttp://folding.bmc.uu.se
--
gmx-users mailing listgmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
* Only plain text messages are allowed!
* 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] Topology file for simulation of inorganic compounds

2012-08-13 Thread Justin Lemkul



On 8/13/12 1:49 PM, Deepak Ojha wrote:

Dear All

I want to do simulation of some ions like azide ion,thiocynate ion and
other inorganic compounds in liquid water.
I am not sure how to generate topology file for the same.How can I go
about doing the same.



Likely you will have to do it manually, which means you will have to understand 
Chapter 5 of the manual and the proper parameterization protocol for your 
molecules under the chosen force field.


Some web servers exist for such purposes, depending on the force field you want 
to use.  Otherwise consult the following page and be prepared for a significant 
time investment:


http://www.gromacs.org/Documentation/How-tos/Parameterization

-Justin

--


Justin A. Lemkul, Ph.D.
Research Scientist
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
* Only plain text messages are allowed!
* 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] Topology file for simulation of inorganic compounds

2012-08-13 Thread Deepak Ojha
Dear All

I want to do simulation of some ions like azide ion,thiocynate ion and
other inorganic compounds in liquid water.
I am not sure how to generate topology file for the same.How can I go
about doing the same.

--
DeepaK Ojha
School Of Chemistry

"Selfishness is not living as one wishes to live, it is asking others
to live as one wishes to live"
-- 
gmx-users mailing listgmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
* Only plain text messages are allowed!
* 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] Topology file

2012-07-26 Thread David van der Spoel

On 2012-07-26 22:09, Shima Arasteh wrote:

Hi all,

To check a topology, is it correct to compare topologies generated by 2 
different  force fields? For example with gmx and charmm36. I know the topology 
generated by gmx is correct, however the parameters might be different.

You can not compare directly in this way, since different force field 
may have different number of e.g. dihedral terms.  Never use the gmx 
force field by the way.



Thanks in advance .


Sincerely,
Shima




--
David van der Spoel, Ph.D., Professor of Biology
Dept. of Cell & Molec. Biol., Uppsala University.
Box 596, 75124 Uppsala, Sweden. Phone:  +46184714205.
sp...@xray.bmc.uu.sehttp://folding.bmc.uu.se
--
gmx-users mailing listgmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
* Only plain text messages are allowed!
* 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] Topology file

2012-07-26 Thread Shima Arasteh
Hi all,

To check a topology, is it correct to compare topologies generated by 2 
different  force fields? For example with gmx and charmm36. I know the topology 
generated by gmx is correct, however the parameters might be different. 

Thanks in advance .


Sincerely,
Shima
--
gmx-users mailing listgmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
* Only plain text messages are allowed!
* 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] topology file in DPPC

2012-05-15 Thread Justin A. Lemkul



On 5/15/12 11:11 AM, scapr...@uniroma3.it wrote:

Do you know where  I should look to find the solution to my problem?



The entire procedure for modifying the force field is different, as the 
directory structure and file names are very different.  The best solution is to 
upgrade your Gromacs version and follow the tutorial directly.


-Justin

--


Justin A. Lemkul, Ph.D.
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] topology file in DPPC

2012-05-15 Thread scaprari
Do you know where  I should look to find the solution to my problem?

Silvia



>
> On 5/15/12 10:47 AM, scapr...@uniroma3.it wrote:
>> Dear all,
>> I'm meticulously following the tutorial KALP-15 in DPPC in order to
>> carry
>> on the simulation of a protein of mine in a extended patch of DPPC. I
>> have
>> already modified the ffnonbonded.itp and ffbonded.itp. At this point I'm
>> reading that I should change my topology file but I'm not able to find
>> within my topology file (obtained by using pdb2gmx on my protein)
>> #include"gromos53a6.ff/forcefield.itp"statement so as to replace this
>> line
>> with #include "gromos53a6_lipid.ff/forcefield.itp", as reported in the
>> tutorial.
>> Indeed, at the top of my topology file I have got;
>>
>> ; Include forcefield parameters
>> #include "ffG53a6.itp"
>>
>> and, after the list of all atoms,  at the bottom of the file;
>>
>> ; Include Position restraint file
>> #ifdef POSRES
>> #include "posre.itp"
>> #endif
>>
>> ; Include water topology
>> #include "spc.itp"
>>
>> #ifdef POSRES_WATER
>> ; Position restraint for each water oxygen
>> [ position_restraints ]
>> ;  i funct   fcxfcyfcz
>> 11   1000   1000   1000
>> #endif
>>
>> ; Include generic topology for ions
>> #include "ions.itp"
>>
>> [ system ]
>> ; Name
>> PHOSPHOLIPASE A2, AMMODYTOXIN A
>>
>> [ molecules ]
>> ; Compound#mols
>> Protein_A   1
>>
>> I'm wondering if it may depend on the fact I have got the version
>> gromacs-4.0.7 and not the version 4.5.3 or newer or if I have to change
>> some options when I launch the command pdb2gmx in order to build my
>> topology file. I used the following command:
>>
>> pdb2gmx -f 3G8G.pdb -o 3G8G_membr.gro  -p 3G8G_membr.top  -i
>> 3G8G_membr.itp -ignh -water spc
>>
>> and I chose GROMOS96 53A6
>>
>> Please, let me know.
>>
>
> The problem is you're using an old version.  The tutorial states that you
> are
> expected to be using a version in the 4.5.x series.  If you don't, then
> the
> procedure will be somewhat different.
>
> -Justin
>
> --
> 
>
> Justin A. Lemkul, Ph.D.
> 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 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] topology file in DPPC

2012-05-15 Thread scaprari
Dear all,
I'm meticulously following the tutorial KALP-15 in DPPC in order to carry
on the simulation of a protein of mine in a extended patch of DPPC. I have
already modified the ffnonbonded.itp and ffbonded.itp. At this point I'm
reading that I should change my topology file but I'm not able to find
within my topology file (obtained by using pdb2gmx on my protein)
#include"gromos53a6.ff/forcefield.itp"statement so as to replace this line
with #include "gromos53a6_lipid.ff/forcefield.itp", as reported in the
tutorial.
Indeed, at the top of my topology file I have got;

; Include forcefield parameters
#include "ffG53a6.itp"

and, after the list of all atoms,  at the bottom of the file;

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

; Include water topology
#include "spc.itp"

#ifdef POSRES_WATER
; Position restraint for each water oxygen
[ position_restraints ]
;  i funct   fcxfcyfcz
   11   1000   1000   1000
#endif

; Include generic topology for ions
#include "ions.itp"

[ system ]
; Name
PHOSPHOLIPASE A2, AMMODYTOXIN A

[ molecules ]
; Compound#mols
Protein_A   1

I'm wondering if it may depend on the fact I have got the version  
gromacs-4.0.7 and not the version 4.5.3 or newer or if I have to change
some options when I launch the command pdb2gmx in order to build my
topology file. I used the following command:

pdb2gmx -f 3G8G.pdb -o 3G8G_membr.gro  -p 3G8G_membr.top  -i
3G8G_membr.itp -ignh -water spc

and I chose GROMOS96 53A6

Please, let me know.

Silvia

-- 
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] topology file in DPPC

2012-05-15 Thread Justin A. Lemkul



On 5/15/12 10:47 AM, scapr...@uniroma3.it wrote:

Dear all,
I'm meticulously following the tutorial KALP-15 in DPPC in order to carry
on the simulation of a protein of mine in a extended patch of DPPC. I have
already modified the ffnonbonded.itp and ffbonded.itp. At this point I'm
reading that I should change my topology file but I'm not able to find
within my topology file (obtained by using pdb2gmx on my protein)
#include"gromos53a6.ff/forcefield.itp"statement so as to replace this line
with #include "gromos53a6_lipid.ff/forcefield.itp", as reported in the
tutorial.
Indeed, at the top of my topology file I have got;

; Include forcefield parameters
#include "ffG53a6.itp"

and, after the list of all atoms,  at the bottom of the file;

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

; Include water topology
#include "spc.itp"

#ifdef POSRES_WATER
; Position restraint for each water oxygen
[ position_restraints ]
;  i funct   fcxfcyfcz
11   1000   1000   1000
#endif

; Include generic topology for ions
#include "ions.itp"

[ system ]
; Name
PHOSPHOLIPASE A2, AMMODYTOXIN A

[ molecules ]
; Compound#mols
Protein_A   1

I'm wondering if it may depend on the fact I have got the version
gromacs-4.0.7 and not the version 4.5.3 or newer or if I have to change
some options when I launch the command pdb2gmx in order to build my
topology file. I used the following command:

pdb2gmx -f 3G8G.pdb -o 3G8G_membr.gro  -p 3G8G_membr.top  -i
3G8G_membr.itp -ignh -water spc

and I chose GROMOS96 53A6

Please, let me know.



The problem is you're using an old version.  The tutorial states that you are 
expected to be using a version in the 4.5.x series.  If you don't, then the 
procedure will be somewhat different.


-Justin

--


Justin A. Lemkul, Ph.D.
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] Topology parameters

2012-05-06 Thread Shima Arasteh
All right. 

Thank you so much.
Shima




 From: Justin A. Lemkul 
To: Discussion list for GROMACS users  
Sent: Sunday, May 6, 2012 11:32 PM
Subject: Re: [gmx-users] Topology parameters
 


On 5/6/12 2:36 PM, Shima Arasteh wrote:
> IS the parameters for a residue in CHARMM36 and CHARMM27 the same?
>

I would suggest you do a bit of reading about the CHARMM force field hierarchy 
and why they are named as they are.  If you're parameterizing a new residue you 
can do it for either, or both, depending on the parameter set.  In your case, 
I'd suspect your formyl parameters are going to be universal between CHARMM27 
and CHARMM36, but don't let my suggestion be your justification.  If a reviewer 
asks, "some guy on the Internet said it" is not sufficient justification ;)

-Justin

> Thanks in advance,
> Shima
>
> 
> *From:* Justin A. Lemkul 
> *To:* Discussion list for GROMACS users 
> *Sent:* Saturday, May 5, 2012 9:16 PM
> *Subject:* Re: [gmx-users] Topology parameters
>
>
>
> On 5/5/12 11:39 AM, Shima Arasteh wrote:
>  > Can I use CHARMM-GUI to get the formyl parameters for CHARMM36 and 
>CHARMM27?
>  >
>
> Maybe, try it and see. Swiss-Param might also be an option.
>
> -Justin
>
> -- 
>
> Justin A. Lemkul
> Ph.D. Candidate
> ICTAS Doctoral Scholar
> MILES-IGERT Trainee
> Department of Biochemistry
> Virginia Tech
> Blacksburg, VA
> jalemkul[at]vt.edu <http://vt.edu> | (540) 231-9080
> http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin
>
> 
> -- 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
>
>

-- 


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] Topology parameters

2012-05-06 Thread Justin A. Lemkul



On 5/6/12 2:36 PM, Shima Arasteh wrote:

IS the parameters for a residue in CHARMM36 and CHARMM27 the same?



I would suggest you do a bit of reading about the CHARMM force field hierarchy 
and why they are named as they are.  If you're parameterizing a new residue you 
can do it for either, or both, depending on the parameter set.  In your case, 
I'd suspect your formyl parameters are going to be universal between CHARMM27 
and CHARMM36, but don't let my suggestion be your justification.  If a reviewer 
asks, "some guy on the Internet said it" is not sufficient justification ;)


-Justin


Thanks in advance,
Shima


*From:* Justin A. Lemkul 
*To:* Discussion list for GROMACS users 
*Sent:* Saturday, May 5, 2012 9:16 PM
*Subject:* Re: [gmx-users] Topology parameters



On 5/5/12 11:39 AM, Shima Arasteh wrote:
 > Can I use CHARMM-GUI to get the formyl parameters for CHARMM36 and CHARMM27?
 >

Maybe, try it and see. Swiss-Param might also be an option.

-Justin

-- 

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


-- 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




--


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] Topology parameters

2012-05-06 Thread Shima Arasteh
IS the parameters for a residue in CHARMM36 and CHARMM27 the same?

Thanks in advance,

Shima




 From: Justin A. Lemkul 
To: Discussion list for GROMACS users  
Sent: Saturday, May 5, 2012 9:16 PM
Subject: Re: [gmx-users] Topology parameters
 


On 5/5/12 11:39 AM, Shima Arasteh wrote:
> Can I use CHARMM-GUI to get the formyl parameters for CHARMM36 and CHARMM27?
> 

Maybe, try it and see.  Swiss-Param might also be an option.

-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 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] Topology parameters

2012-05-05 Thread Justin A. Lemkul



On 5/5/12 11:39 AM, Shima Arasteh wrote:

Can I use CHARMM-GUI to get the formyl parameters for CHARMM36 and CHARMM27?



Maybe, try it and see.  Swiss-Param might also be an option.

-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] Topology parameters

2012-05-05 Thread Shima Arasteh
Can I use CHARMM-GUI  to get the formyl parameters for CHARMM36 and CHARMM27?

Cheers,
Shima




 From: Justin A. Lemkul 
To: Shima Arasteh ; Discussion list for GROMACS 
users  
Sent: Saturday, May 5, 2012 5:58 PM
Subject: Re: [gmx-users] Topology parameters
 


On 5/5/12 9:21 AM, Shima Arasteh wrote:
> Dear gmx users,
> 
> I need to define a new residue topology ( Formyl ) in CHARMM27 force field. I
> have the formyl parameters for GROMOS. I'd like to know if the residue
> parameters are the same in different force fields?
> 

Almost certainly not.  Each force field has different parameterization 
procedures and makes use of different atom types.  A new force field requires a 
new topology.

-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] Topology parameters

2012-05-05 Thread Justin A. Lemkul



On 5/5/12 9:21 AM, Shima Arasteh wrote:

Dear gmx users,

I need to define a new residue topology ( Formyl ) in CHARMM27 force field. I
have the formyl parameters for GROMOS. I'd like to know if the residue
parameters are the same in different force fields?



Almost certainly not.  Each force field has different parameterization 
procedures and makes use of different atom types.  A new force field requires a 
new topology.


-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


[gmx-users] Topology parameters

2012-05-05 Thread Shima Arasteh
Dear gmx users,

I need to define a new residue topology ( Formyl ) in CHARMM27 force field. I 
have the formyl parameters for GROMOS. I'd like to know if the residue 
parameters are the same in different force fields?

Thanks,
Shima
-- 
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] Topology files

2012-02-28 Thread kirtana S
Thanks for the reply. As I am new to GROMACS I am not able to understand. I
will go through the manual and  will get back if I have further questions.


On Tue, Feb 28, 2012 at 6:00 PM, Justin A. Lemkul  wrote:

>
>
> kirtana S wrote:
>
>> Thanks for the documentation. Again I have a  ligand terminated silica
>> nanoparticle ,what can be the easiest way to approach  the topology file
>> for this.
>>
>>
>>
> I don't know what the structure is, but I'll assume the .rtp route is
> probably easier than writing a (presumably) repetitive topology for a
> relatively large particle.  You may or may not need to invoke specbond.dat
> for cross-links, as pdb2gmx can deal only with linear species via .rtp
> entries; other bonds are introduced with specbond.dat.
>
> -Justin
>
>
>> On Tue, Feb 28, 2012 at 5:46 PM, Justin A. Lemkul > jalem...@vt.edu>> wrote:
>>
>>
>>
>>kirtana S wrote:
>>
>>For example if I have a polystyrene polymer chains and monomers
>>of user defined solvent . I need to construct the .itp file by
>>hand .
>>
>>
>>Polystyrene is most easily generated using .rtp entries in pdb2gmx.
>> The solvent can be described by manually created .itp files.
>>
>>
>> http://www.gromacs.org/__**Documentation/How-tos/Polymers
>>
>> 
>> **>
>>
>>
>>-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
>>
>> 
>> >
>>
>>
>>
> --
> ==**==
>
> 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/Searchbefore
>  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] Topology files

2012-02-28 Thread Justin A. Lemkul



kirtana S wrote:
Thanks for the documentation. Again I have a  ligand terminated silica 
nanoparticle ,what can be the easiest way to approach  the topology file 
for this.





I don't know what the structure is, but I'll assume the .rtp route is probably 
easier than writing a (presumably) repetitive topology for a relatively large 
particle.  You may or may not need to invoke specbond.dat for cross-links, as 
pdb2gmx can deal only with linear species via .rtp entries; other bonds are 
introduced with specbond.dat.


-Justin



On Tue, Feb 28, 2012 at 5:46 PM, Justin A. Lemkul > wrote:




kirtana S wrote:

For example if I have a polystyrene polymer chains and monomers
of user defined solvent . I need to construct the .itp file by
hand .


Polystyrene is most easily generated using .rtp entries in pdb2gmx.
 The solvent can be described by manually created .itp files.

http://www.gromacs.org/__Documentation/How-tos/Polymers



-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





--


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] Topology files

2012-02-28 Thread kirtana S
Thanks for the documentation. Again I have a  ligand terminated silica
nanoparticle ,what can be the easiest way to approach  the topology file
for this.



On Tue, Feb 28, 2012 at 5:46 PM, Justin A. Lemkul  wrote:

>
>
> kirtana S wrote:
>
>> For example if I have a polystyrene polymer chains and monomers of user
>> defined solvent . I need to construct the .itp file by hand .
>>
>
> Polystyrene is most easily generated using .rtp entries in pdb2gmx.  The
> solvent can be described by manually created .itp files.
>
> http://www.gromacs.org/**Documentation/How-tos/Polymers
>
>
> -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/Searchbefore
>  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] Topology files

2012-02-28 Thread Justin A. Lemkul



kirtana S wrote:
For example if I have a polystyrene polymer chains and monomers of user 
defined solvent . I need to construct the .itp file by hand . 



Polystyrene is most easily generated using .rtp entries in pdb2gmx.  The solvent 
can be described by manually created .itp files.


http://www.gromacs.org/Documentation/How-tos/Polymers

-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] Topology files

2012-02-28 Thread kirtana S
For example if I have a polystyrene polymer chains and monomers of user
defined solvent . I need to construct the .itp file by hand .

Thanks for the reply
Kirtana

On Tue, Feb 28, 2012 at 5:31 PM, Justin A. Lemkul  wrote:

>
>
> kirtana S wrote:
>
>> Dear Users,
>>
>> In the gromacs manual on means to convert the pdb file to gromacs
>> topology file I need to use pdb2gmx. Here I have to specify the force field
>> .How to work with organic residues which do not exists . Since I had used
>> general amber force field, and could not find this in the /share/top/
>> directory . Do I need to manually construct the .itp file for each
>> residues. What will be the easiest way to construct
>> a solvated box with user defined residues in gromacs.
>>
>>
> If the residue of interest will be a component of a polymeric molecule
> like a protein, refer to http://www.gromacs.org/**
> Documentation/How-tos/Adding_**a_Residue_to_a_Force_Field.
>  If the molecule is a standalone species like a ligand, solvent, etc then
> all you have to do is #include an .itp file in the system topology.
>
> -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/Searchbefore
>  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] Topology files

2012-02-28 Thread Justin A. Lemkul



kirtana S wrote:

Dear Users,

In the gromacs manual on means to convert the pdb file to gromacs 
topology file I need to use pdb2gmx. 
Here I have to specify the force field .How to work with organic 
residues which do not exists . 
Since I had used general amber force field, and could not find this in 
the /share/top/ directory . 
Do I need to manually construct the .itp file for each residues. What 
will be the easiest way to construct

a solvated box with user defined residues in gromacs.



If the residue of interest will be a component of a polymeric molecule like a 
protein, refer to 
http://www.gromacs.org/Documentation/How-tos/Adding_a_Residue_to_a_Force_Field. 
 If the molecule is a standalone species like a ligand, solvent, etc then all 
you have to do is #include an .itp file in the system topology.


-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


[gmx-users] Topology files

2012-02-28 Thread kirtana S
Dear Users,

In the gromacs manual on means to convert the pdb file to gromacs topology
file I need to use pdb2gmx.
Here I have to specify the force field .How to work with organic residues
which do not exists .
Since I had used general amber force field, and could not find this in the
/share/top/ directory .
Do I need to manually construct the .itp file for each residues. What will
be the easiest way to construct
a solvated box with user defined residues in gromacs.

Thanks
Kirtana
-- 
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] topology missing

2012-01-25 Thread Mark Abraham

On 25/01/2012 9:39 PM, Anik Sen wrote:

Hi,
   I am using GROMACS 4.5.5. I made the pdb file of DNA and 
changed its residue files according to the given dna.rtp.
But still am getting an error 'OH' residue topology file missing where 
as the 'OH' is present in the atomtype as well as the residue topology 
file.


Whats the reason?


Copy and paste the actual error message, along with your command line, 
and any of the other output that you think might be relevant. It's a 
waste of time to give help based on information filtered through 
someone's head, because you don't know when their error might have 
occurred (if any!).


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] topology missing

2012-01-25 Thread Anik Sen
Hi,
   I am using GROMACS 4.5.5. I made the pdb file of DNA and changed its 
residue files according to the given dna.rtp.
But still am getting an error 'OH' residue topology file missing where as the  
'OH' is present in the atomtype as well as the residue topology file.

Whats the reason?

Thanking n advance
Anik

Anik Sen
Student
CSIR-Central Salt & Marine Chemicals Research Institute,
Gijubhai Badheka Marg.
Bhavnagar, Gujarat 364002
[www.csmcri.org]

-- 
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] topology and parameter for cobalt metal

2011-12-14 Thread Justin A. Lemkul



shahid nayeem wrote:

Dear all
I am interested in simulating a model protein with cobalt forming 
coordinate bond with His. I coud'nt find entry for cobalt in any 
forcefield .rtp file. I am using Gromacs 4.5.4 . If any one can help me 
in getting forcefield parameters for charmm27/ OPLSaa/Amber99 in gromacs 
format please respond. Please suggest where else I should search for these.


Standard, non-polarizable MM force fields are not particularly well-suited for 
dealing with transition metals.  Fixed point charges on such species, especially 
those coordinated to amino acids, are generally very incorrect, not to mention 
geometric coordination parameters that are not built-in.


http://www.gromacs.org/Documentation/How-tos/Parameterization#Exotic_Species

-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


[gmx-users] topology and parameter for cobalt metal

2011-12-14 Thread shahid nayeem
Dear all
I am interested in simulating a model protein with cobalt forming
coordinate bond with His. I coud'nt find entry for cobalt in any forcefield
.rtp file. I am using Gromacs 4.5.4 . If any one can help me in getting
forcefield parameters for charmm27/ OPLSaa/Amber99 in gromacs format please
respond. Please suggest where else I should search for these.
Thanking all
shahid Nayeem
-- 
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] Topology error

2011-11-21 Thread Mark Abraham

On 21/11/2011 9:18 PM, archana sonawani wrote:

Dear Sir/Madam,

For orienting the protein and membrane, I have used the .mdp and 
following topology file:

;
;   File 'topol_popc.top' was generated
;   By user: jalemkul (502)
;   On host:bevany.biochem.vt.edu  
;   At date: Fri Oct 20 13:26:53 2006
;
;   This is your topology file
;
; Include chain topologies
#include "gromos53a6_lipid.ff/forcefield.itp"
#include "popc.itp"

; Include water topology
#include "gromos53a6_lipid.ff/spc.itp"

; Include ion topologies
#include "gromos53a6_lipid.ff/ions.itp"

; System specifications
[ system ]
60-Lipid POPC Bilayer

[ molecules ]
; molecule name nr.
POPC 60
SOL 1536
After running grompp to generate .tpr file I get following error:

Program grompp, VERSION 4.5.4

Source code file: gmxcpp.c, line: 248

Fatal error:

Topology include file "gromos53a6_lipid.ff/forcefield.itp" not found


But in the above topology file, if I change the

#include "gromos53a6_lipid.ff/forcefield.itp" line to *;*include 
"gromos53a6_lipid.ff/forcefield.itp"




This file is not optional, so failing to include anything makes further 
problems. You need such a file. You are doing something non-standard, so 
your burden of following directions, and/or understanding how the 
include file and force field mechanism works is a bit higher than normal :-)


Mark



I get following error:


Program grompp, VERSION 4.5.4

Source code file: topio.c, line: 653

Fatal error:

Syntax error - File popc.itp, line 1

Last line read:

'[ moleculetype ]'

Invalid order for directive moleculetype


Can anyone help me out...Thanks in adv




-- 
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] Topology error

2011-11-21 Thread archana sonawani
Dear Sir/Madam,

For orienting the protein and membrane, I have used the .mdp and following
topology file:

;
;   File 'topol_popc.top' was generated
;   By user: jalemkul (502)
;   On host: bevany.biochem.vt.edu
;   At date: Fri Oct 20 13:26:53 2006
;
;   This is your topology file
;
; Include chain topologies
#include "gromos53a6_lipid.ff/forcefield.itp"
#include "popc.itp"

; Include water topology
#include "gromos53a6_lipid.ff/spc.itp"

; Include ion topologies
#include "gromos53a6_lipid.ff/ions.itp"

; System specifications
[ system ]
60-Lipid POPC Bilayer

[ molecules ]
; molecule name nr.
POPC 60
SOL 1536

After running grompp to generate .tpr file I get following error:

Program grompp, VERSION 4.5.4

Source code file: gmxcpp.c, line: 248



Fatal error:

Topology include file "gromos53a6_lipid.ff/forcefield.itp" not found


But in the above topology file, if I change the

 #include "gromos53a6_lipid.ff/forcefield.itp" line to *;*include
"gromos53a6_lipid.ff/forcefield.itp"


I get following error:


Program grompp, VERSION 4.5.4

Source code file: topio.c, line: 653



Fatal error:

Syntax error - File popc.itp, line 1

Last line read:

'[ moleculetype ]'

Invalid order for directive moleculetype


Can anyone help me out...Thanks in adv
-- 
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] Topology dihedral or improper dihedral

2011-11-13 Thread Justin A. Lemkul



XUEMING TANG wrote:
Thank you for your reply, Lina. I did try the prodrg server before. For 
A molecule, there is no improper dihedral. A just has one dihedral for 
the head group. For B molecule, it gives improper dihedral of the tetra 
nitrogen. I am wondering what is the criteria or at least a guidance of 
what could be the best choice for the topology file. Or any choices 
could be right and I have to try which one provide consistent results 
with experiment by trial?  



I see no reason why a quaternary nitrogen would need an improper dihedral 
assigned to it.  Some impropers in Gromos96 are applied to tetrahedral groups 
like alpha-carbons because there is no explicit H atom present, thus the 
stereochemistry is maintained with an improper so that the groups do not invert. 
 With four explicit substituents, an improper should not be necessary.


-Justin


Best!
Xueming

On Sun, Nov 13, 2011 at 9:57 AM, lina > wrote:


On Sun, Nov 13, 2011 at 10:26 AM, XUEMING TANG mailto:yujan2...@gmail.com>> wrote:
 > Hi there
 >
 > I am trying to build up a topology file for a tetra nitrogen
molecule using
 > Gromos UA and OPLS UA force field. Following are the molecules.
 >
 > 1. for  C
 >   |
 >  A.  C-C-C-C-C-C-C-N-C and
 >   |
 >  C
 >
 >   C
 >|
 >  B.   C-C-C-C-C-C-C-N-C-C-C-C-COOH
 >|
 >   C
 >
 > I have several questions:
 >
 > For A, the first molecule, should I use improper dihedral for the
N-(CH3)3
 > or proper dihedral for C-C-N-C with multi 3 for the head group?
Is there any
 > differences when using different force field (GROMOS and OPLS-UA)?
 > For B, the second molecule, should I use improper dihedral for
the N-(CH3)3
 > or proer dihedral for C-C-N-C with multi 2? Is there any
differences when
 > using different force field (GROMOS UA and OPLS UA)?
 >
 > The first expression about improper didedral is a planer
restriction. When
 > read the manual it also refers to the tetrahedral structure like
this tetra
 > nitrogen functional group. My question is when it has
 > really long hydrocarbon tail, I guess it will be hard to flip over to
 > the mirror images. Or it is not true in the simulation?

You may take one built by other server, such as
http://davapc1.bioch.dundee.ac.uk/prodrg/

as reference, and built yourself if necessary.
 >
 > Thank you very much!!!
 >
 > Best!
 > Xueming
 >
 >
 >
 >
 >
 > --
 > 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




--


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] Topology dihedral or improper dihedral

2011-11-13 Thread XUEMING TANG
Thank you for your reply, Lina. I did try the prodrg server before. For A
molecule, there is no improper dihedral. A just has one dihedral for the
head group. For B molecule, it gives improper dihedral of the tetra
nitrogen. I am wondering what is the criteria or at least a guidance of
what could be the best choice for the topology file. Or any choices could
be right and I have to try which one provide consistent results with
experiment by trial?

Best!
Xueming

On Sun, Nov 13, 2011 at 9:57 AM, lina  wrote:

> On Sun, Nov 13, 2011 at 10:26 AM, XUEMING TANG 
> wrote:
> > Hi there
> >
> > I am trying to build up a topology file for a tetra nitrogen
> molecule using
> > Gromos UA and OPLS UA force field. Following are the molecules.
> >
> > 1. for  C
> >   |
> >  A.  C-C-C-C-C-C-C-N-C and
> >   |
> >  C
> >
> >   C
> >|
> >  B.   C-C-C-C-C-C-C-N-C-C-C-C-COOH
> >|
> >   C
> >
> > I have several questions:
> >
> > For A, the first molecule, should I use improper dihedral for the
> N-(CH3)3
> > or proper dihedral for C-C-N-C with multi 3 for the head group? Is there
> any
> > differences when using different force field (GROMOS and OPLS-UA)?
> > For B, the second molecule, should I use improper dihedral for the
> N-(CH3)3
> > or proer dihedral for C-C-N-C with multi 2? Is there any differences when
> > using different force field (GROMOS UA and OPLS UA)?
> >
> > The first expression about improper didedral is a planer restriction.
> When
> > read the manual it also refers to the tetrahedral structure like this
> tetra
> > nitrogen functional group. My question is when it has
> > really long hydrocarbon tail, I guess it will be hard to flip over to
> > the mirror images. Or it is not true in the simulation?
>
> You may take one built by other server, such as
> http://davapc1.bioch.dundee.ac.uk/prodrg/
>
> as reference, and built yourself if necessary.
> >
> > Thank you very much!!!
> >
> > Best!
> > Xueming
> >
> >
> >
> >
> >
> > --
> > 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
>
-- 
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] Topology dihedral or improper dihedral

2011-11-13 Thread lina
On Sun, Nov 13, 2011 at 10:26 AM, XUEMING TANG  wrote:
> Hi there
>
> I am trying to build up a topology file for a tetra nitrogen molecule using
> Gromos UA and OPLS UA force field. Following are the molecules.
>
> 1. for     C
>   |
>  A.  C-C-C-C-C-C-C-N-C  and
>   |
>  C
>
>   C
>        |
>  B.   C-C-C-C-C-C-C-N-C-C-C-C-COOH
>    |
>   C
>
> I have several questions:
>
> For A, the first molecule, should I use improper dihedral for the N-(CH3)3
> or proper dihedral for C-C-N-C with multi 3 for the head group? Is there any
> differences when using different force field (GROMOS and OPLS-UA)?
> For B, the second molecule, should I use improper dihedral for the N-(CH3)3
> or proer dihedral for C-C-N-C with multi 2? Is there any differences when
> using different force field (GROMOS UA and OPLS UA)?
>
> The first expression about improper didedral is a planer restriction. When
> read the manual it also refers to the tetrahedral structure like this tetra
> nitrogen functional group. My question is when it has
> really long hydrocarbon tail, I guess it will be hard to flip over to
> the mirror images. Or it is not true in the simulation?

You may take one built by other server, such as
http://davapc1.bioch.dundee.ac.uk/prodrg/

as reference, and built yourself if necessary.
>
> Thank you very much!!!
>
> Best!
> Xueming
>
>
>
>
>
> --
> 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] Topology dihedral or improper dihedral

2011-11-12 Thread XUEMING TANG
Hi there

I am trying to build up a topology file for a tetra nitrogen molecule using
Gromos UA and OPLS UA force field. Following are the molecules.

1. for C
  |
 A.  C-C-C-C-C-C-C-N-C  and
  |
 C

  C
   |
 B.   C-C-C-C-C-C-C-N-C-C-C-C-COOH
   |
  C

I have several questions:

For A, the first molecule, should I use improper dihedral for the N-(CH3)3
or proper dihedral for C-C-N-C with multi 3 for the head group? Is there
any differences when using different force field (GROMOS and OPLS-UA)?
For B, the second molecule, should I use improper dihedral for the N-(CH3)3
or proer dihedral for C-C-N-C with multi 2? Is there any differences when
using different force field (GROMOS UA and OPLS UA)?

The first expression about improper didedral is a planer restriction. When
read the manual it also refers to the tetrahedral structure like this tetra
nitrogen functional group. My question is when it has
really long hydrocarbon tail, I guess it will be hard to flip over to
the mirror images. Or it is not true in the simulation?

Thank you very much!!!

Best!
Xueming
-- 
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] Topology for a new ligand

2011-08-11 Thread Mark Abraham

On 12/08/2011 1:48 AM, Kavyashree M wrote:

Dear users,

I was mentioning about OPLSAA force field for ATP
and other small molecule.  I just wanted to know the
procedure to be followed and some guidance from
people who have created topologies for such
molecules manually. I am going through chapter 5 of
the manual. But wanted some useful suggestions.


The link you mentioned, and links from it contain all the 
readily-available information. For ATP, you may well find a validated 
topology in the published literature.


Mark



Thank you
With Regards
M. Kavyashree


On Thu, Aug 11, 2011 at 5:03 PM, Justin A. Lemkul > wrote:




Kavyashree M wrote:

Dear gromacs users,

I wanted to know the steps to be followed
in order to generate a topology for a new
ligand. I went through the mailing list and
http://www.gromacs.org/Documentation/How-tos/Parameterization
but was not clear.


All force fields are different, and since you haven't said which
one you're trying to use there's nothing that anyone can tell you.
 Read the primary literature for the force field you want to use
and follow the procedure laid out therein.

-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 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] Topology for a new ligand

2011-08-11 Thread Justin A. Lemkul



Kavyashree M wrote:

Dear users,

I was mentioning about OPLSAA force field for ATP
and other small molecule.  I just wanted to know the
procedure to be followed and some guidance from
people who have created topologies for such
molecules manually. I am going through chapter 5 of
the manual. But wanted some useful suggestions.



Methods for parameterizing molecules for OPLS-AA are described in fairly 
extensive detail in the 2001 OPLS-AA paper cited in the manual.  The manual 
itself will only tell you how to implement parameters in Gromacs, not how to 
derive them.


-Justin


Thank you
With Regards
M. Kavyashree


On Thu, Aug 11, 2011 at 5:03 PM, Justin A. Lemkul > wrote:




Kavyashree M wrote:

Dear gromacs users,

I wanted to know the steps to be followed
in order to generate a topology for a new
ligand. I went through the mailing list and
http://www.gromacs.org/__Documentation/How-tos/__Parameterization

but was not clear.


All force fields are different, and since you haven't said which one
you're trying to use there's nothing that anyone can tell you.  Read
the primary literature for the force field you want to use and
follow the procedure laid out therein.

-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





--


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] Topology for a new ligand

2011-08-11 Thread Kavyashree M
Dear users,

I was mentioning about OPLSAA force field for ATP
and other small molecule.  I just wanted to know the
procedure to be followed and some guidance from
people who have created topologies for such
molecules manually. I am going through chapter 5 of
the manual. But wanted some useful suggestions.

Thank you
With Regards
M. Kavyashree


On Thu, Aug 11, 2011 at 5:03 PM, Justin A. Lemkul  wrote:

>
>
> Kavyashree M wrote:
>
>> Dear gromacs users,
>>
>> I wanted to know the steps to be followed
>> in order to generate a topology for a new
>> ligand. I went through the mailing list and
>> http://www.gromacs.org/**Documentation/How-tos/**Parameterization
>> but was not clear.
>>
>>
> All force fields are different, and since you haven't said which one you're
> trying to use there's nothing that anyone can tell you.  Read the primary
> literature for the force field you want to use and follow the procedure laid
> out therein.
>
> -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/Searchbefore
>  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] Topology for a new ligand

2011-08-11 Thread Justin A. Lemkul



Kavyashree M wrote:

Dear gromacs users,

I wanted to know the steps to be followed
in order to generate a topology for a new
ligand. I went through the mailing list and
http://www.gromacs.org/Documentation/How-tos/Parameterization
but was not clear.



All force fields are different, and since you haven't said which one you're 
trying to use there's nothing that anyone can tell you.  Read the primary 
literature for the force field you want to use and follow the procedure laid out 
therein.


-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


[gmx-users] Topology for a new ligand

2011-08-11 Thread Kavyashree M
Dear gromacs users,

I wanted to know the steps to be followed
in order to generate a topology for a new
ligand. I went through the mailing list and
http://www.gromacs.org/Documentation/How-tos/Parameterization
but was not clear.

Awaiting your suggestions

Thanking you
With regards
M. Kavyashree
-- 
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: Res: Res: [gmx-users] Topology for nBuOH

2011-07-13 Thread Justin A. Lemkul
   
 [ atoms ]  
;   nrtype   resnr  residuatomcgnrcharge  mass
 1  Hb  1   BuOH HO 1   0.38   1.0080   
 2  O   1   BuOH OH 1  -0.66  15.9994   
 3  C   1   BuOH CA 1   0.40  12.0110   
 4  C   1   BuOH CB 1  -0.22  12.0110   
 5  C   1   BuOH CC 1  -0.22  12.0110   
 6  C   1   BuOH CD 1  -0.33  12.0110   
 7  H   1   BuOH HA 1  -0.06   1.0080   
 8  H   1   BuOH Hb 1  -0.06   1.0080   
 9  H   1   BuOH HC 1   0.11   1.0080   
 10 H   1   BuOH HD 1   0.11   1.0080   
 11 H   1   BuOH HE 1   0.11   1.0080   
 12 H   1   BuOH HF 1   0.11   1.0080   
 13 H   1   BuOH HG 1   0.11   1.0080   
 14 H   1   BuOH HF 1   0.11   1.0080   
 15 H   1   BuOH HG 1   0.11   1.0080   
 
 [ dihedrals ] 
  1  2  3  4  30.41840   1.25520   0.0  -1.67360   0.0   0.0

  2  3  4  5  32.87441   0.58158   2.09200  -5.54799   0.0   0.0
  3  4  5  6  32.92880  -1.46440   0.20920  -1.67360   0.0   0.0
 
 [ system ]

 ; Name
 BuOH in water
 
 [ molecules ]

 ; Compound#mols
 BuOH   343
 SOL185


*De:* Justin A. Lemkul 
*Para:* Discussion list for GROMACS users 
*Enviadas:* Quarta-feira, 13 de Julho de 2011 12:10:03
*Assunto:* Re: Res: [gmx-users] Topology for nBuOH



Juliana Angeiras wrote:
 > Hi,
 >  really the problem was with the order, but now i have another 
problem (below). I have already read the chapter 5, but i didn't fix my 
problem. I would like to use an archive .itp for the butanol, although i 
use the tip3p for water in the system (mixture water-butanol), I want 
use different parameters for nonbonded interactions between 
water-butanol (because were parametrized specifically for these 
interactions). How can I do that? My atual archive .top is more below.
 > 

Again, your directives are out of order.  The [nonbond_params] are part 
of the top-level force field definition and thus must come before the 
molecule definitions.


-Justin

 > Archive .top:
 >  [ defaults ]
 > ; nbfunccomb-rule  gen-pairs  fudgeLJ fudgeQQ
 >13yes0.50.5
 >  [ atomtypes ]
 > ; name  at.num  mass  charge  ptype  V(c6)  W(c12)
 >Hb  11.00800.38  A  0.4149611E-07  0.9737356E-04
 >O  815.9994-0.66  A  0.9722461E-06  0.1247853E-02
 >C  612.01100.40  A  0.4895634E-05  0.2791672E-02
 >C  612.0110-0.22  A  0.4895634E-05  0.2791672E-02C  6  
  12.0110-0.33  A  0.4895634E-05  0.2791672E-02

 >H  11.0080-0.06  A  0.4415162E-06  0.5131711E-03
 >H  11.00800.11  A  0.4415162E-06  0.5131711E-03
 >  [ bondtype ]
  ;  ai  aj funct  c0  c1
Hb  O  1  0.10900  462750.4  
  O  C  1  0.14300  376560.0
C  C  1  0.15300  265265.6  
  H  C  1  0.10800  284512.0
 >  [ angletype ]
;  aiajak func  c0c1
COHb  1108.000460.240
  CCO  1108.000418.400  
CCC  1112.000488.273
  HCH  1109.470292.000
 >  [ moleculetype ]
; name  nrexcl  
BuOH3

[ atoms ]
  ;  nrtype  resnr  residuatomcgnrcharge  mass
 >  1  Hb  1  BuOHHO1  0.38  1.0080  
  2  O  1  BuOHOH1  -0.66  15.99943  
C  1  BuOHCA1  0.40  12.01104  
C  1  BuOHCB1  -0.22  12.01105  C  
1  BuOHCC1  -0.22  12.01106  C  1  
BuOH 

Re: Res: [gmx-users] Topology for nBuOH

2011-07-13 Thread Justin A. Lemkul
0.1247853E-02
  O   OW  1  0.1562282E-05  0.1785000E-02
  O   Hb  1  0.000E+00  0.000E+00
  OW  OW  1  0.2510400E-05  0.2552000E-02
  OW  Hb  1  0.000E+00  0.000E+00
  Hb  Hb  1  0.4149611E-07  0.9737356E-04
 
 [ system ]

 ; Name
 BuOH in water
 
 [ molecules ]

 ; Compound#mols
 BuOH   343
 SOL185
 
Error with grompp:
 
WARNING 1 [file BuOHaa.top, line 27]:

  Overriding atomtype C

WARNING 2 [file BuOHaa.top, line 28]:
  Overriding atomtype C

WARNING 3 [file BuOHaa.top, line 30]:
  Overriding atomtype H

ERROR 1 [file BuOHaa.top, line 32]:
  Invalid directive bondtype

WARNING 4 [file BuOHaa.top, line 34]:
  Too few parameters on line (source file toppush.c, line 246)

WARNING 5 [file BuOHaa.top, line 35]:
  Too few parameters on line (source file toppush.c, line 246)

WARNING 6 [file BuOHaa.top, line 36]:
  Too few parameters on line (source file toppush.c, line 246)

WARNING 7 [file BuOHaa.top, line 37]:
  Too few parameters on line (source file toppush.c, line 246)

ERROR 2 [file BuOHaa.top, line 39]:
  Invalid directive angletype

WARNING 8 [file BuOHaa.top, line 41]:
  Too few parameters on line (source file toppush.c, line 315)

WARNING 9 [file BuOHaa.top, line 42]:
  Too few parameters on line (source file toppush.c, line 315)

WARNING 10 [file BuOHaa.top, line 43]:
  Too few parameters on line (source file toppush.c, line 315)

WARNING 11 [file BuOHaa.top, line 44]:
  Too few parameters on line (source file toppush.c, line 315)
Generated 10 of the 10 non-bonded parameter combinations
Generating 1-4 interactions: fudge = 0.5
Generated 10 of the 10 1-4 parameter combinations
---
Program grompp, VERSION 4.5.4
Source code file: toppush.c, line: 1166
Fatal error:
Atomtype opls_111 not found
Thanks,
 
Juliana Angeiras B. da Silva.

Laboratório de Química Teórica e Computacional
Departamento de Química Fundamental
Universidade Federal de Pernambuco
Recife, PE - Brasil.


*De:* Justin A. Lemkul 
*Para:* Discussion list for GROMACS users 
*Enviadas:* Terça-feira, 12 de Julho de 2011 12:42:49
*Assunto:* Re: [gmx-users] Topology for nBuOH



Juliana Angeiras wrote:
 > Hi,
 >  I'm trying to include parameters of n-butanol in the topology file 
(below) (i'm not using opls or other force field for butanol), but when 
i run the grompp program, I obtain the following message:

 >  Fatal error:
 > Atomtype Hb not found
 >  How can I fix this problem?
 > 

Your directives are out of order.  If you define new atomtypes, etc. 
they have to be defined before they can be used.  See Chapter 5 of the 
manual.


-Justin

 > Thanks
 >  Topology file:
 >  ; Include forcefield parameters#include 
"oplsaa.ff/forcefield.itp"; Include forcefield parameters
#include "oplsaa.ff/tip5p.itp"  [ moleculetype ]
 > ; name  nrexcl BuOH 3  [ atoms ];  nrtype  resnr  
residuatomcgnrcharge  mass

 >  1  Hb  1  BuOHHO1  0.38  1.0080
 >  2  O  1  BuOHOH1  -0.66  15.9994
 >  3  C  1  BuOHCA1  0.40  12.0110
 >  4  C  1  BuOHCB1  -0.22  12.0110
 >  5  C  1  BuOHCC1  -0.22  12.0110
 >  6  C  1  BuOHCD1  -0.33  12.0110
 >  7  H  1  BuOHHA1  -0.06  1.0080
 >  8  H  1  BuOHHB1  -0.06  1.0080
 >  9  H  1  BuOHHC1  0.11  1.0080
 >  10H  1  BuOHHD1  0.11  1.0080
 >  11H  1  BuOHHE1  0.11  1.0080
 >  12H  1  BuOHHF1  0.11  1.0080
 >  13H  1  BuOHHG1  0.11  1.0080
 >  14H  1  BuOHHF1  0.11  1.0080
 >  15H  1  BuOHHG1  0.11  1.0080
 >  [ atomtypes ]
 > ; name  at.num  mass  charge  ptype  V(c6)  W(c12)
 >Hb  11.00800.38  A  0.4149611E-07  0.9737356E-04
 >O  815.9994-0.66  A  0.9722461E-06  0.1247853E-02
 >C  612.01100.40  A  0.4895634E-05  0.2791672E-02
 >C  612.0110-0.22  A  0.4895634E-05  0.2791672E-02C  6  
  12.0110-0.33  A  0.4895634E-05  0.2791672E-02

 >H  11.0080-0.06  A  0.4415162E-06  0.5131711E-03
 >H  11.00800.11  A  0.4415162E-06  0.5131711E-03
 >  [ bondtype ]  ;  ai  aj funct  c0  c1
HB  O  1  0.10900  462750.4  O  C  1  0.14300  
376560.0C  C  1  0.15300  265265.6  H  C  1  
0.10800  284512.0
 >  [ angletype ];  aiajak func  c0
c1COHB  1108.000460.240  CCO  1  
  108.000418

Res: [gmx-users] Topology for nBuOH

2011-07-13 Thread Juliana Angeiras
0.000E+00
  Hb  Hb  1  0.4149611E-07  0.9737356E-04

 [ system ]
 ; Name
 BuOH in water

 [ molecules ]
 ; Compound    #mols
 BuOH   343
 SOL    185


Error with grompp:

WARNING 1 [file BuOHaa.top, line 27]:
  Overriding atomtype C

WARNING 2 [file BuOHaa.top, line 28]:
  Overriding atomtype C

WARNING 3 [file BuOHaa.top, line 30]:
  Overriding atomtype H

ERROR 1 [file BuOHaa.top, line 32]:
  Invalid directive bondtype

WARNING 4 [file BuOHaa.top, line 34]:
  Too few parameters on line (source file toppush.c, line 246)

WARNING 5 [file BuOHaa.top, line 35]:
  Too few parameters on line (source file toppush.c, line 246)

WARNING 6 [file BuOHaa.top, line 36]:
  Too few parameters on line (source file toppush.c, line 246)

WARNING 7 [file BuOHaa.top, line 37]:
  Too few parameters on line (source file toppush.c, line 246)

ERROR 2 [file BuOHaa.top, line 39]:
  Invalid directive angletype

WARNING 8 [file BuOHaa.top, line 41]:
  Too few parameters on line (source file toppush.c, line 315)

WARNING 9 [file BuOHaa.top, line 42]:
  Too few parameters on line (source file toppush.c, line 315)

WARNING 10 [file BuOHaa.top, line 43]:
  Too few parameters on line (source file toppush.c, line 315)

WARNING 11 [file BuOHaa.top, line 44]:
  Too few parameters on line (source file toppush.c, line 315)
Generated 10 of the 10 non-bonded parameter combinations
Generating 1-4 interactions: fudge = 0.5
Generated 10 of the 10 1-4 parameter combinations
---
Program grompp, VERSION 4.5.4
Source code file: toppush.c, line: 1166
Fatal error:
Atomtype opls_111 not found

Thanks,

Juliana Angeiras B. da Silva.
Laboratório de Química Teórica e Computacional
Departamento de Química Fundamental
Universidade Federal de Pernambuco
Recife, PE - Brasil.





De: Justin A. Lemkul 
Para: Discussion list for GROMACS users 
Enviadas: Terça-feira, 12 de Julho de 2011 12:42:49
Assunto: Re: [gmx-users] Topology for nBuOH



Juliana Angeiras wrote:
> Hi,
>  I'm trying to include parameters of n-butanol in the topology file (below) 
>(i'm not using opls or other force field for butanol), but when i run the 
>grompp 
>program, I obtain the following message:
>  Fatal error:
> Atomtype Hb not found
>  How can I fix this problem?
>  

Your directives are out of order.  If you define new atomtypes, etc. they have 
to be defined before they can be used.  See Chapter 5 of the manual.

-Justin

> Thanks
>  Topology file:
>  ; Include forcefield parameters    #include "oplsaa.ff/forcefield.itp"       
> 
>; Include forcefield parameters    #include "oplsaa.ff/tip5p.itp"      [ 
>moleculetype ]
> ; name  nrexcl BuOH 3      [ atoms ]            ;  nr    type  resnr  residu  
>   
>atom    cgnr        charge          mass
>      1      Hb  1  BuOH    HO    1      0.38      1.0080
>      2      O  1  BuOH    OH    1      -0.66      15.9994
>      3      C  1  BuOH    CA    1      0.40      12.0110
>      4      C  1  BuOH    CB    1      -0.22      12.0110
>      5      C  1  BuOH    CC    1      -0.22      12.0110
>      6      C  1  BuOH    CD    1      -0.33      12.0110
>      7      H  1  BuOH    HA    1      -0.06      1.0080
>      8      H  1  BuOH    HB    1      -0.06      1.0080
>      9      H  1  BuOH    HC    1      0.11      1.0080
>      10    H  1  BuOH    HD    1      0.11      1.0080
>      11    H  1  BuOH    HE    1      0.11      1.0080
>      12    H  1  BuOH    HF    1      0.11      1.0080
>      13    H  1  BuOH    HG    1      0.11      1.0080
>      14    H  1  BuOH    HF    1      0.11      1.0080
>      15    H  1  BuOH    HG    1      0.11      1.0080
>  [ atomtypes ]
> ; name  at.num  mass  charge  ptype  V(c6)  W(c12)
>    Hb  1    1.0080    0.38  A  0.4149611E-07  0.9737356E-04
>    O  8    15.9994    -0.66  A  0.9722461E-06  0.1247853E-02
>    C  6    12.0110    0.40  A  0.4895634E-05  0.2791672E-02
>    C  6    12.0110    -0.22  A  0.4895634E-05  0.2791672E-02    C  6    
>12.0110    -0.33  A  0.4895634E-05  0.2791672E-02
>    H  1    1.0080    -0.06  A  0.4415162E-06  0.5131711E-03
>    H  1    1.0080    0.11  A  0.4415162E-06  0.5131711E-03
>  [ bondtype ]              ;  ai  aj funct          c0          c1            
>HB  O  1  0.10900  462750.4              O  C  1  0.14300  376560.0            
>C  C  1  0.15300  265265.6          H  C  1  0.10800  284512.0
>  [ angletype ]            ;  ai    aj    ak func          c0        c1        
>C    O    HB  1    108.000    460.240          C    C    O  1    108.000    
>418.400          C    C    C  1    112.000    488.273          H    C    H  1  
>  
>109.470    292.000
>  [ dihedrals ]            1  2  3  4  3    0.41840  1.25520  0.0  
>-1.67360  
>0.0  0.000

Re: [gmx-users] Topology for nBuOH

2011-07-12 Thread Justin A. Lemkul



Juliana Angeiras wrote:

Hi,
 
I'm trying to include parameters of n-butanol in the topology file 
(below) (i'm not using opls or other force field for butanol), but when 
i run the grompp program, I obtain the following message:
 
Fatal error:

Atomtype Hb not found
 
How can I fix this problem?
 


Your directives are out of order.  If you define new atomtypes, etc. they have 
to be defined before they can be used.  See Chapter 5 of the manual.


-Justin


Thanks
 
Topology file:
 
; Include forcefield parameters 
#include "oplsaa.ff/forcefield.itp" 

; Include forcefield parameters 
#include "oplsaa.ff/tip5p.itp"  
[ moleculetype ]
; name  nrexcl 
BuOH 3  
 [ atoms ] 
;   nrtype   resnr  residuatomcgnrcharge  mass

 1  Hb  1   BuOH HO 1   0.38   1.0080
 2  O   1   BuOH OH 1  -0.66  15.9994
 3  C   1   BuOH CA 1   0.40  12.0110
 4  C   1   BuOH CB 1  -0.22  12.0110
 5  C   1   BuOH CC 1  -0.22  12.0110
 6  C   1   BuOH CD 1  -0.33  12.0110
 7  H   1   BuOH HA 1  -0.06   1.0080
 8  H   1   BuOH HB 1  -0.06   1.0080
 9  H   1   BuOH HC 1   0.11   1.0080
 10 H   1   BuOH HD 1   0.11   1.0080
 11 H   1   BuOH HE 1   0.11   1.0080
 12 H   1   BuOH HF 1   0.11   1.0080
 13 H   1   BuOH HG 1   0.11   1.0080
 14 H   1   BuOH HF 1   0.11   1.0080
 15 H   1   BuOH HG 1   0.11   1.0080
 [ atomtypes ]
; name   at.num   mass   charge   ptype   V(c6)   W(c12)
   Hb  1 1.0080 0.38   A   0.4149611E-07  0.9737356E-04
   O   815.9994-0.66   A   0.9722461E-06  0.1247853E-02
   C   612.0110 0.40   A   0.4895634E-05  0.2791672E-02
   C   612.0110-0.22   A   0.4895634E-05  0.2791672E-02  
   C   612.0110-0.33   A   0.4895634E-05  0.2791672E-02

   H   1 1.0080-0.06   A   0.4415162E-06  0.5131711E-03
   H   1 1.0080 0.11   A   0.4415162E-06  0.5131711E-03
 [ bondtype ]   
;  ai   aj funct   c0   c1 
   HB  O   1   0.10900   462750.4
   O   C   1   0.14300   376560.0 
   C   C   1   0.15300   265265.6
   H   C   1   0.10800   284512.0
 [ angletype ] 
;  aiajak func   c0c1   
  COHB  1108.000460.240
  CCO   1108.000418.400
  CCC   1112.000488.273
  HCH   1109.470292.000
 [ dihedrals ]   
  1  2  3  4  30.41840   1.25520   0.0  -1.67360   0.0   0.0

  2  3  4  5  32.87441   0.58158   2.09200  -5.54799   0.0   0.0
  3  4  5  6  32.92880  -1.46440   0.20920  -1.67360   0.0   0.0
 [ nonbond_params ]
; i   j   func   V(c6)   W(c12)
  H   H   1  0.4149611E-07  0.9737356E-04
  H   HW  1  0.4149611E-07  0.9737356E-04
  H   C   1  0.4415162E-06  0.5131711E-03
  H   O   1  0.1927397E-06  0.3416812E-03
  H   OW  1  0.1927397E-06  0.3416812E-03
  H   Hb  1  0.4149611E-07  0.9737356E-04
  HW  HW  1  0.000E+00  0.000E+00
  HW  C   1  0.4415162E-06  0.5131711E-03
  HW  O   1  0.000E+00  0.000E+00
  HW  OW  1  0.000E+00  0.000E+00
  HW  Hb  1  0.4149611E-07  0.9737356E-04
  C   C   1  0.4895634E-05  0.2791672E-02
  C   O   1  0.2159724E-05  0.1866537E-02
  C   OW  1  0.2159724E-05  0.1866537E-02
  C   Hb  1  0.4415162E-06  0.5131711E-03
  O   O   1  0.9722461E-06  0.1247853E-02
  O   OW  1  0.1562282E-05  0.1785000E-02
  O   Hb  1  0.000E+00  0.000E+00
  OW  OW  1  0.2510400E-05  0.2552000E-02
  OW  Hb  1  0.000E+00  0.000E+00
  Hb  Hb  1  0.4149611E-07  0.9737356E-04
 [ system ]
 ; Name
 BuOH in water
 [ molecules ]
 ; Compound#mols
 BuOH   343
 SOL185
 
 
 



--


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] Topology for nBuOH

2011-07-12 Thread Juliana Angeiras
Hi,

I'm trying to include parameters of n-butanol in the topology file (below) (i'm 
not using opls or other force field for butanol), but when i run the grompp 
program, I obtain the following message:

Fatal error:
Atomtype Hb not found

How can I fix this problem?

Thanks

Topology file:

; Include forcefield parameters  
#include "oplsaa.ff/forcefield.itp"  
 
; Include forcefield parameters  
#include "oplsaa.ff/tip5p.itp"   
[ moleculetype ]
; name  nrexcl 
BuOH 3   
 [ atoms ]      
;   nr    type   resnr  residu    atom    cgnr    charge  mass
 1  Hb  1   BuOH HO 1   0.38   1.0080
 2  O   1   BuOH OH 1  -0.66  15.9994
 3  C   1   BuOH CA 1   0.40  12.0110
 4  C   1   BuOH CB 1  -0.22  12.0110
 5  C   1   BuOH CC 1  -0.22  12.0110
 6  C   1   BuOH CD 1  -0.33  12.0110
 7  H   1   BuOH HA 1  -0.06   1.0080
 8  H   1   BuOH HB 1  -0.06   1.0080
 9  H   1   BuOH HC 1   0.11   1.0080
 10 H   1   BuOH HD 1   0.11   1.0080
 11 H   1   BuOH HE 1   0.11   1.0080
 12 H   1   BuOH HF 1   0.11   1.0080
 13 H   1   BuOH HG 1   0.11   1.0080
 14 H   1   BuOH HF 1   0.11   1.0080
 15 H   1   BuOH HG 1   0.11   1.0080
 [ atomtypes ]
; name   at.num   mass   charge   ptype   V(c6)   W(c12)
   Hb  1 1.0080 0.38   A   0.4149611E-07  0.9737356E-04
   O   8    15.9994    -0.66   A   0.9722461E-06  0.1247853E-02
   C   6    12.0110 0.40   A   0.4895634E-05  0.2791672E-02
   C   6    12.0110    -0.22   A   0.4895634E-05  0.2791672E-02   
   C   6    12.0110    -0.33   A   0.4895634E-05  0.2791672E-02
   H   1 1.0080    -0.06   A   0.4415162E-06  0.5131711E-03
   H   1 1.0080 0.11   A   0.4415162E-06  0.5131711E-03
 [ bondtype ]    
;  ai   aj funct   c0   c1      
   HB  O   1   0.10900   462750.4     
   O   C   1   0.14300   376560.0  
   C   C   1   0.15300   265265.6 
   H   C   1   0.10800   284512.0
 [ angletype ]      
;  ai    aj    ak func   c0    c1    
  C    O    HB  1    108.000    460.240 
  C    C    O   1    108.000    418.400     
  C    C    C   1    112.000    488.273 
  H    C    H   1    109.470    292.000
 [ dihedrals ]    
  1  2  3  4  3    0.41840   1.25520   0.0  -1.67360   0.0   0.0
  2  3  4  5  3    2.87441   0.58158   2.09200  -5.54799   0.0   0.0
  3  4  5  6  3    2.92880  -1.46440   0.20920  -1.67360   0.0   0.0
 [ nonbond_params ]
; i   j   func   V(c6)   W(c12)
  H   H   1  0.4149611E-07  0.9737356E-04
  H   HW  1  0.4149611E-07  0.9737356E-04
  H   C   1  0.4415162E-06  0.5131711E-03
  H   O   1  0.1927397E-06  0.3416812E-03
  H   OW  1  0.1927397E-06  0.3416812E-03
  H   Hb  1  0.4149611E-07  0.9737356E-04
  HW  HW  1  0.000E+00  0.000E+00
  HW  C   1  0.4415162E-06  0.5131711E-03
  HW  O   1  0.000E+00  0.000E+00
  HW  OW  1  0.000E+00  0.000E+00
  HW  Hb  1  0.4149611E-07  0.9737356E-04
  C   C   1  0.4895634E-05  0.2791672E-02
  C   O   1  0.2159724E-05  0.1866537E-02
  C   OW  1  0.2159724E-05  0.1866537E-02
  C   Hb  1  0.4415162E-06  0.5131711E-03
  O   O   1  0.9722461E-06  0.1247853E-02
  O   OW  1  0.1562282E-05  0.1785000E-02
  O   Hb  1  0.000E+00  0.000E+00
  OW  OW  1  0.2510400E-05  0.2552000E-02
  OW  Hb  1  0.000E+00  0.000E+00
  Hb  Hb  1  0.4149611E-07  0.9737356E-04
 [ system ]
 ; Name
 BuOH in water
 [ molecules ]
 ; Compound    #mols
 BuOH   343
 SOL    185-- 
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] Topology

2011-03-11 Thread Justin A. Lemkul



chitrala kumaraswamynaidu wrote:
Iam working with a NMR protein containing Zn ions. After generating 
topology file using pdb2gmx, it is generating topology files containing 
Protein and ions separately. How can i solve this.




pdb2gmx -chainsep (-merge in old versions).

-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


[gmx-users] Topology

2011-03-11 Thread chitrala kumaraswamynaidu
Iam working with a NMR protein containing Zn ions. After generating topology
file using pdb2gmx, it is generating topology files containing Protein and
ions separately. How can i solve this.
-- 
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] topology files for ligands in MD Simulation under OPLS AA force field

2011-02-24 Thread Michael Brunsteiner

mircial at sjtu.edu.cn mircial at sjtu.edu.cn wrote

> I am using GROMACS package to do molecular dynamics
> simulations under OPLS_AA force field. I encounter
> some problems when preparing the topology files
> of small molecules (ligands).My questions are as follows: 
> 1, how to chose the atom type of each atom from the ligands? 
amber + GAFF would be easier, here the ligand parameterization
is pretty straight forward (there's probably a HOWTO)

with OPLS you want to work with homology, and if you don't need
any high throughput you can do it manually ... look up 
groups in amino acids that are similar to your the components
of your ligands, and choose the atomtypes accordingly.

> 2, how to define the charges of each atom? I know that, when
> the atom type is defined, there will be a corresponding charge
> to this atom type. Does it safe to use this charge in the
> simulations on the ligands (since these charges are designed
> for amino acids)?

yes, OPLS was mostly used for peptides, but AFAIK, the
original parameterization in Jorgensens group started with generic
organic molecules, there's certainly papers from this group that
present parameters for non-peptide molecules, e.g., various
hetero cycles ... (I am not sure whether the gmx topology folder
only contains the part of OPLS relevant for peptides, or all of it,
if not ask someone from the Jorgenson group, they might give you
the original parameter files)
However, if you have some generic organic molecules and assign
OPLS atom types there is no guarantee that the charges add up to
zero (for a neutral molecule) ... what i did before, with
reasonable success, is combine EPSD (SCF or AM1BCC) charges with
OPLS parameters...

In any case, as has been pointed out before many times in this list
mixing different FFs, or using ad-hoc values for some parameters, is 
dangerous, and you really want to test your molecules before using 
them in any extensive calculations. 
A good start might be to compare structures that were minimized
with your FF to the same structure minimized at the DFT or MP2 level
(HF-SCF might do as well depending on the molecule) ... but that's
probably not enough ...

> 3,When the atom type and the charge of atom are defined,
> how to prepare the file in the GROMACS format? Does there easier
> method to prepare such files than manually? 
As for the bonded interactions, in some cases pdb2gmx and
grompp can figure them out by themselves. You'd assign atom- 
(thereby bond-) types, and the bonds, then make a new entry
in the rtp file, including only the atoms and bonds sections. you 
choose a new 3 letter residue name, and starting from a PDB
file of the ligand (with this residue name, and the same atom order
and the same atom names as in the rtp entry) use pdb2gmx/grompp.
However often enough, especially for proper and improper dihedrals
you have to define them manually (again by anlogy) 
for anything but the simplest molecules its a painstaking task ...
so you either take the time, or you use amber (or you talk to 
Schrodinger, they will have some script for assigning OPLS parameters
but that probably wont work with Gromacs)

good luck!
mic


  
-- 
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] topology files for ligands in MD Simulation under OPLS AA force field

2011-02-24 Thread Justin A. Lemkul



mirc...@sjtu.edu.cn wrote:

Dear All:

I am using GROMACS package to do molecular dynamics simulations under 
OPLS_AA force field. I encounter some problems when preparing the 
topology files of small molecules (ligands).My questions are as follows:


1, how to chose the atom type of each atom from the ligands?

2, how to define the charges of each atom? I know that, when the atom 
type is defined, there will be a corresponding charge to this atom type. 
Does it safe to use this charge in the simulations on the ligands (since 
these charges are designed for amino acids)?




Proper parameterization is a difficult task.  See, for instance:

http://www.gromacs.org/Documentation/How-tos/Parameterization

3,When the atom type and the charge of atom are defined, how to prepare 
the file in the GROMACS format? Does there easier method to prepare such 
files than manually?




There is at least one script (topolgen.pl) in the User Contributions section of 
the Gromacs site, but it is not terribly intelligent and the resulting topology 
should not be viewed as something that you should use without proper validation.


http://www.gromacs.org/Downloads/User_contributions/Other_software

Otherwise, break out Chapter 5 of the manual (something you should probably do 
anyway) and your favorite text editor.


-Justin

Thank you very much for your time and your kindness. I really appreciate 
if any one can share their experience in preparing the topology files 
for ligands used under OPLS_AA force field and GROMACS package.


Sincerely Yours

Ruo-Xu Gu




--


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] topology files for ligands in MD Simulation under OPLS AA force field

2011-02-24 Thread mircial

Dear All:

I am using GROMACS package to do molecular dynamics simulations under  
OPLS_AA force field. I encounter some problems when preparing the  
topology files of small molecules (ligands).My questions are as follows:


1, how to chose the atom type of each atom from the ligands?

2, how to define the charges of each atom? I know that, when the atom  
type is defined, there will be a corresponding charge to this atom  
type. Does it safe to use this charge in the simulations on the  
ligands (since these charges are designed for amino acids)?


3,When the atom type and the charge of atom are defined, how to  
prepare the file in the GROMACS format? Does there easier method to  
prepare such files than manually?


Thank you very much for your time and your kindness. I really  
appreciate if any one can share their experience in preparing the  
topology files for ligands used under OPLS_AA force field and GROMACS  
package.


Sincerely Yours

Ruo-Xu Gu


--
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] topology with all hydrogens

2011-02-06 Thread Tsjerk Wassenaar
Hi Mohsen,

Now I wouldn't want to get in the way of excellent ideas :)
Check the PRODRG FAQ. It is explained there (somewhere) how to force
inclusion of hydrogens.

Cheers,

Tsjerk

On Feb 6, 2011 11:39 AM, "mohsen ramezanpour" 
wrote:

Dear Dr.Tsjerk

I know it.I know PRODRG is not very valid for evaluating free
energies(According to Dr.Justin article) too.
But the problem is here:
If I know all of hydrigenes are polar,how can I have a topology with all of
them?
Actually there are so many polar H in my small molecules which are not
included in topology file!!
These can change the number of Hydrogen bonds and as a resut my binding free
energy!!!


I want to do a new work about these,I have something excellent in my mind.
I only need a topology file with all hydrogenes in 43A1 force field.
Is there any way(except manually adding hydrogenes)?
Please let me know.
I will present  a new and excellent method instead of using PRODRG
Thanks in advance
Mohsen

On Sun, Feb 6, 2011 at 1:23 PM, Tsjerk Wassenaar  wrote:
> > Hi Mohsen, > > So ...

--
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] topology with all hydrogens

2011-02-06 Thread mohsen ramezanpour
Dear Dr.Tsjerk

I know it.I know PRODRG is not very valid for evaluating free
energies(According to Dr.Justin article) too.
But the problem is here:
If I know all of hydrigenes are polar,how can I have a topology with all of
them?
Actually there are so many polar H in my small molecules which are not
included in topology file!!
These can change the number of Hydrogen bonds and as a resut my binding free
energy!!!


I want to do a new work about these,I have something excellent in my mind.
I only need a topology file with all hydrogenes in 43A1 force field.
Is there any way(except manually adding hydrogenes)?
Please let me know.
I will present  a new and excellent method instead of using PRODRG
Thanks in advance
Mohsen

On Sun, Feb 6, 2011 at 1:23 PM, Tsjerk Wassenaar  wrote:

> Hi Mohsen,
>
> So you want to have a topology with all atoms for use in a united atom
> force field? Seems a bit awkward... PRODRG will only include polar
> hydrogens, in accordance with the GROMOS force field. It won't be a
> particular trustworthy topology, though, but you've probably caught the
> words of caution expressed regularly on this list.
>
> Cheers,
>
> Tsjerk
>
> On Feb 6, 2011 9:42 AM, "mohsen ramezanpour" 
> wrote:
>
> Dear All
>
> I want to generate a topology file in 43A1 force field for a small
> molecule.
> Of course I want one which contains all hydrogens in that.
> PRODRG doesn't generate like this.do you know another server with these
> conditions?
>
> Thanks in advance
> Mohsen
>
> --
> 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
>
-- 
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] topology with all hydrogens

2011-02-06 Thread Tsjerk Wassenaar
Hi Mohsen,

So you want to have a topology with all atoms for use in a united atom force
field? Seems a bit awkward... PRODRG will only include polar hydrogens, in
accordance with the GROMOS force field. It won't be a particular trustworthy
topology, though, but you've probably caught the words of caution expressed
regularly on this list.

Cheers,

Tsjerk

On Feb 6, 2011 9:42 AM, "mohsen ramezanpour" 
wrote:

Dear All

I want to generate a topology file in 43A1 force field for a small molecule.
Of course I want one which contains all hydrogens in that.
PRODRG doesn't generate like this.do you know another server with these
conditions?

Thanks in advance
Mohsen

--
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

[gmx-users] topology with all hydrogens

2011-02-06 Thread mohsen ramezanpour
Dear All

I want to generate a topology file in 43A1 force field for a small molecule.
Of course I want one which contains all hydrogens in that.
PRODRG doesn't generate like this.do you know another server with these
conditions?

Thanks in advance
Mohsen
-- 
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] topology file for two molecules

2010-09-28 Thread Gavin Melaugh
Justin

Cheers, I had just posted a more general query prior to you getting
back? It is perhaps not as vague.

Gavin

Justin A. Lemkul wrote:
>
>
> Gavin Melaugh wrote:
>> Hi Justin
>>
>> I Did the following and it worked. You were right about using one atom
>> types directive. If you don't mind I have another query. In all my
>> topology files to date (one molecule type), I have always listed the
>> atomtypes in the [atomtypes] directive, and used the set up in
>> [defaults] below. Is this all that is required for the calculation of
>> the LJ energies. I have never used a pair list, a pair types directive
>> or a nonbond_params directive. From the output I still get LJ energies
>> but I am slightly confused by all the different options.
>>
>
> I can't give a better explanation than what is already in the manual,
> specifically section 5.3.3 ([nonbond_params] and whether or not you
> need it). Whether or not you need [pairs] depends on the underlying
> mechanics of the force field you're using.  See section 5.3.5.
>
> -Justin
>
>> parameter level
>> [defaults]
>> ; nbfunccomb-rule gen-pairsfudgeLJ fudgeQQ
>>  1 3  yes0.5 0.5
>>
>> [atomtypes]
>> ;type mass   charge  ptype sigma(nm)   
>> epsilon(kjmol-1)
>>CB 12.011000  0.00   A  0.355000  0.292880
>>CA 12.011000 -0.115000   A  0.355000  0.292880
>>HC  1.008000  0.115000   A  0.242000  0.125520
>>CU 13.019000  0.265000   A  0.35  0.334720
>>NU 14.007000 -0.597000   A  0.325000  0.711280
>>CH 13.019000  0.332000   A  0.385000  0.334720
>>C3 15.035000  0.00   A  0.390500  0.732200
>>CT 12.011000 -0.18   A  0.35  0.276144
>>HH  1.008000  0.06   A  0.25  0.125520
>>
>> #include "tertpentane.itp"
>> #include "methylcage.itp"
>>
>> ;System level
>> [system]
>> ;the name of this system
>> cage and pentane
>>
>> [molecules]
>> ; moleculename number
>>  tertpentane  1
>>  methylcage   1
>>
>>
>> Justin A. Lemkul wrote:
>>>
>>> Gavin Melaugh wrote:
>>>> Hi Justin
>>>>
>>>> Cheers. I have attached my pentane itp file. The cage itp is formatted
>>>> the same. In the examples there never seems to be [atom types]
>>>> included.
>>>> I have included it. Does it seem O.K?
>>>>
>>> You only need to use an [atomtypes] directive if you're introducing
>>> new atom types into the force field.  I don't know if it's
>>> necessary. I also don't know if these sections can be discontinuous
>>> (i.e., one
>>> section present in cage.itp, another in tertpentane.itp).  Try the
>>> approach I sent you and see.  If you get fatal errors, then you may
>>> have to assemble a single [atomtypes] directive and put it either in
>>> the first .itp file you #include, or before any of the #includes in
>>> the .top (since #include really just means "cut-and-paste the stuff
>>> here," it works the same).
>>>
>>> -Justin
>>>
>>>> Gavin
>>>>
>>>> Justin A. Lemkul wrote:
>>>>> Gavin Melaugh wrote:
>>>>>> Hi Mark.
>>>>>>
>>>>>> I have looked this up, but I am a little uncertain. Does one of the
>>>>>> topologies have to be included within the other as an .itp file.
>>>>>> If so
>>>>> No.  You #include .itp files to simply add molecules within a
>>>>> .top; it
>>>>> would be very rare to #include an .itp within an .itp (although it
>>>>> can
>>>>> be done).  For your case, your .top is simple:
>>>>>
>>>>> #include (force field)
>>>>>
>>>>> #include "cage.itp"
>>>>> #include "pentane.itp"
>>>>>
>>>>> [ system ]
>>>>> My system
>>>>>
>>>>> [ molecules ]
>>>>> Cage1
>>>>> Pentane 1
>>>>>
>>>>>> do I only define the [defaults directive] in the .top file i.e. I
>>>>>> take
>>>>>> it I don't include this directi

Re: [gmx-users] topology file for two molecules

2010-09-28 Thread Justin A. Lemkul



Gavin Melaugh wrote:

Hi Justin

I Did the following and it worked. You were right about using one atom
types directive. If you don't mind I have another query. In all my
topology files to date (one molecule type), I have always listed the
atomtypes in the [atomtypes] directive, and used the set up in
[defaults] below. Is this all that is required for the calculation of
the LJ energies. I have never used a pair list, a pair types directive
or a nonbond_params directive. From the output I still get LJ energies
but I am slightly confused by all the different options.



I can't give a better explanation than what is already in the manual, 
specifically section 5.3.3 ([nonbond_params] and whether or not you need it). 
Whether or not you need [pairs] depends on the underlying mechanics of the force 
field you're using.  See section 5.3.5.


-Justin


parameter level
[defaults]
; nbfunccomb-rule gen-pairsfudgeLJ fudgeQQ
 1 3  yes0.5 0.5

[atomtypes]
;type mass   charge  ptype sigma(nm)
epsilon(kjmol-1)

   CB 12.011000  0.00   A  0.355000  0.292880
   CA 12.011000 -0.115000   A  0.355000  0.292880
   HC  1.008000  0.115000   A  0.242000  0.125520
   CU 13.019000  0.265000   A  0.35  0.334720
   NU 14.007000 -0.597000   A  0.325000  0.711280
   CH 13.019000  0.332000   A  0.385000  0.334720
   C3 15.035000  0.00   A  0.390500  0.732200
   CT 12.011000 -0.18   A  0.35  0.276144
   HH  1.008000  0.06   A  0.25  0.125520

#include "tertpentane.itp"
#include "methylcage.itp"

;System level
[system]
;the name of this system
cage and pentane

[molecules]
; moleculename number
 tertpentane  1
 methylcage   1


Justin A. Lemkul wrote:


Gavin Melaugh wrote:

Hi Justin

Cheers. I have attached my pentane itp file. The cage itp is formatted
the same. In the examples there never seems to be [atom types] included.
I have included it. Does it seem O.K?


You only need to use an [atomtypes] directive if you're introducing
new atom types into the force field.  I don't know if it's necessary. 
I also don't know if these sections can be discontinuous (i.e., one

section present in cage.itp, another in tertpentane.itp).  Try the
approach I sent you and see.  If you get fatal errors, then you may
have to assemble a single [atomtypes] directive and put it either in
the first .itp file you #include, or before any of the #includes in
the .top (since #include really just means "cut-and-paste the stuff
here," it works the same).

-Justin


Gavin

Justin A. Lemkul wrote:

Gavin Melaugh wrote:

Hi Mark.

I have looked this up, but I am a little uncertain. Does one of the
topologies have to be included within the other as an .itp file. If so

No.  You #include .itp files to simply add molecules within a .top; it
would be very rare to #include an .itp within an .itp (although it can
be done).  For your case, your .top is simple:

#include (force field)

#include "cage.itp"
#include "pentane.itp"

[ system ]
My system

[ molecules ]
Cage1
Pentane 1


do I only define the [defaults directive] in the .top file i.e. I take
it I don't include this directive in the .itp file as well.


The [defaults] directive is provided by whatever force field you
#include.  If you're using some custom parameters, then it needs to be
defined once (and only once) at the top of your .top file.  Otherwise,
you'll get fatal errors from grompp.

-Justin


Cheers

Gavin

Mark Abraham wrote:

Look up the include file mechanism on the webpage.

Mark

----- Original Message -
From: Gavin Melaugh 
Date: Tuesday, September 28, 2010 19:04
Subject: [gmx-users] topology file for two molecules
To: Discussion list for GROMACS users 


Hi all

I want to run a simulation of one tert-pentane molecule inside one
organic cage molecule. I have the topology files for both
molecules but
obviously I only need one .top file as the input for grompp.
What is the
most convenient way of having both topologies in one file?

Many Thanks

Gavin
--
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 





--


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

Re: [gmx-users] topology file for two molecules

2010-09-28 Thread Gavin Melaugh
Hi Justin

I Did the following and it worked. You were right about using one atom
types directive. If you don't mind I have another query. In all my
topology files to date (one molecule type), I have always listed the
atomtypes in the [atomtypes] directive, and used the set up in
[defaults] below. Is this all that is required for the calculation of
the LJ energies. I have never used a pair list, a pair types directive
or a nonbond_params directive. From the output I still get LJ energies
but I am slightly confused by all the different options.

parameter level
[defaults]
; nbfunccomb-rule gen-pairsfudgeLJ fudgeQQ
 1 3  yes0.5 0.5

[atomtypes]
;type mass   charge  ptype sigma(nm)
epsilon(kjmol-1)
   CB 12.011000  0.00   A  0.355000  0.292880
   CA 12.011000 -0.115000   A  0.355000  0.292880
   HC  1.008000  0.115000   A  0.242000  0.125520
   CU 13.019000  0.265000   A  0.35  0.334720
   NU 14.007000 -0.597000   A  0.325000  0.711280
   CH 13.019000  0.332000   A  0.385000  0.334720
   C3 15.035000  0.00   A  0.390500  0.732200
   CT 12.011000 -0.18   A  0.35  0.276144
   HH  1.008000  0.06   A  0.25  0.125520

#include "tertpentane.itp"
#include "methylcage.itp"

;System level
[system]
;the name of this system
cage and pentane

[molecules]
; moleculename number
 tertpentane  1
 methylcage   1


Justin A. Lemkul wrote:
>
>
> Gavin Melaugh wrote:
>> Hi Justin
>>
>> Cheers. I have attached my pentane itp file. The cage itp is formatted
>> the same. In the examples there never seems to be [atom types] included.
>> I have included it. Does it seem O.K?
>>
>
> You only need to use an [atomtypes] directive if you're introducing
> new atom types into the force field.  I don't know if it's necessary. 
> I also don't know if these sections can be discontinuous (i.e., one
> section present in cage.itp, another in tertpentane.itp).  Try the
> approach I sent you and see.  If you get fatal errors, then you may
> have to assemble a single [atomtypes] directive and put it either in
> the first .itp file you #include, or before any of the #includes in
> the .top (since #include really just means "cut-and-paste the stuff
> here," it works the same).
>
> -Justin
>
>> Gavin
>>
>> Justin A. Lemkul wrote:
>>>
>>> Gavin Melaugh wrote:
>>>> Hi Mark.
>>>>
>>>> I have looked this up, but I am a little uncertain. Does one of the
>>>> topologies have to be included within the other as an .itp file. If so
>>> No.  You #include .itp files to simply add molecules within a .top; it
>>> would be very rare to #include an .itp within an .itp (although it can
>>> be done).  For your case, your .top is simple:
>>>
>>> #include (force field)
>>>
>>> #include "cage.itp"
>>> #include "pentane.itp"
>>>
>>> [ system ]
>>> My system
>>>
>>> [ molecules ]
>>> Cage1
>>> Pentane 1
>>>
>>>> do I only define the [defaults directive] in the .top file i.e. I take
>>>> it I don't include this directive in the .itp file as well.
>>>>
>>> The [defaults] directive is provided by whatever force field you
>>> #include.  If you're using some custom parameters, then it needs to be
>>> defined once (and only once) at the top of your .top file.  Otherwise,
>>> you'll get fatal errors from grompp.
>>>
>>> -Justin
>>>
>>>> Cheers
>>>>
>>>> Gavin
>>>>
>>>> Mark Abraham wrote:
>>>>> Look up the include file mechanism on the webpage.
>>>>>
>>>>> Mark
>>>>>
>>>>> - Original Message -
>>>>> From: Gavin Melaugh 
>>>>> Date: Tuesday, September 28, 2010 19:04
>>>>> Subject: [gmx-users] topology file for two molecules
>>>>> To: Discussion list for GROMACS users 
>>>>>
>>>>>> Hi all
>>>>>>
>>>>>> I want to run a simulation of one tert-pentane molecule inside one
>>>>>> organic cage molecule. I have the topology files for both
>>>>>> molecules but
>>>>>> obviously I only need one .top file as the input for grompp.
>>>>>> What is the

Re: [gmx-users] topology file for two molecules

2010-09-28 Thread Gavin Melaugh
O.K

Cheers

Gavin

Justin A. Lemkul wrote:
>
>
> Gavin Melaugh wrote:
>> Hi Justin
>>
>> Cheers. I have attached my pentane itp file. The cage itp is formatted
>> the same. In the examples there never seems to be [atom types] included.
>> I have included it. Does it seem O.K?
>>
>
> You only need to use an [atomtypes] directive if you're introducing
> new atom types into the force field.  I don't know if it's necessary. 
> I also don't know if these sections can be discontinuous (i.e., one
> section present in cage.itp, another in tertpentane.itp).  Try the
> approach I sent you and see.  If you get fatal errors, then you may
> have to assemble a single [atomtypes] directive and put it either in
> the first .itp file you #include, or before any of the #includes in
> the .top (since #include really just means "cut-and-paste the stuff
> here," it works the same).
>
> -Justin
>
>> Gavin
>>
>> Justin A. Lemkul wrote:
>>>
>>> Gavin Melaugh wrote:
>>>> Hi Mark.
>>>>
>>>> I have looked this up, but I am a little uncertain. Does one of the
>>>> topologies have to be included within the other as an .itp file. If so
>>> No.  You #include .itp files to simply add molecules within a .top; it
>>> would be very rare to #include an .itp within an .itp (although it can
>>> be done).  For your case, your .top is simple:
>>>
>>> #include (force field)
>>>
>>> #include "cage.itp"
>>> #include "pentane.itp"
>>>
>>> [ system ]
>>> My system
>>>
>>> [ molecules ]
>>> Cage1
>>> Pentane 1
>>>
>>>> do I only define the [defaults directive] in the .top file i.e. I take
>>>> it I don't include this directive in the .itp file as well.
>>>>
>>> The [defaults] directive is provided by whatever force field you
>>> #include.  If you're using some custom parameters, then it needs to be
>>> defined once (and only once) at the top of your .top file.  Otherwise,
>>> you'll get fatal errors from grompp.
>>>
>>> -Justin
>>>
>>>> Cheers
>>>>
>>>> Gavin
>>>>
>>>> Mark Abraham wrote:
>>>>> Look up the include file mechanism on the webpage.
>>>>>
>>>>> Mark
>>>>>
>>>>> - Original Message -
>>>>> From: Gavin Melaugh 
>>>>> Date: Tuesday, September 28, 2010 19:04
>>>>> Subject: [gmx-users] topology file for two molecules
>>>>> To: Discussion list for GROMACS users 
>>>>>
>>>>>> Hi all
>>>>>>
>>>>>> I want to run a simulation of one tert-pentane molecule inside one
>>>>>> organic cage molecule. I have the topology files for both
>>>>>> molecules but
>>>>>> obviously I only need one .top file as the input for grompp.
>>>>>> What is the
>>>>>> most convenient way of having both topologies in one file?
>>>>>>
>>>>>> Many Thanks
>>>>>>
>>>>>> Gavin
>>>>>> -- 
>>>>>> 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] topology file for two molecules

2010-09-28 Thread Justin A. Lemkul



Gavin Melaugh wrote:

Hi Justin

Cheers. I have attached my pentane itp file. The cage itp is formatted
the same. In the examples there never seems to be [atom types] included.
I have included it. Does it seem O.K?



You only need to use an [atomtypes] directive if you're introducing new atom 
types into the force field.  I don't know if it's necessary.  I also don't know 
if these sections can be discontinuous (i.e., one section present in cage.itp, 
another in tertpentane.itp).  Try the approach I sent you and see.  If you get 
fatal errors, then you may have to assemble a single [atomtypes] directive and 
put it either in the first .itp file you #include, or before any of the 
#includes in the .top (since #include really just means "cut-and-paste the stuff 
here," it works the same).


-Justin


Gavin

Justin A. Lemkul wrote:


Gavin Melaugh wrote:

Hi Mark.

I have looked this up, but I am a little uncertain. Does one of the
topologies have to be included within the other as an .itp file. If so

No.  You #include .itp files to simply add molecules within a .top; it
would be very rare to #include an .itp within an .itp (although it can
be done).  For your case, your .top is simple:

#include (force field)

#include "cage.itp"
#include "pentane.itp"

[ system ]
My system

[ molecules ]
Cage1
Pentane 1


do I only define the [defaults directive] in the .top file i.e. I take
it I don't include this directive in the .itp file as well.


The [defaults] directive is provided by whatever force field you
#include.  If you're using some custom parameters, then it needs to be
defined once (and only once) at the top of your .top file.  Otherwise,
you'll get fatal errors from grompp.

-Justin


Cheers

Gavin

Mark Abraham wrote:

Look up the include file mechanism on the webpage.

Mark

- Original Message -
From: Gavin Melaugh 
Date: Tuesday, September 28, 2010 19:04
Subject: [gmx-users] topology file for two molecules
To: Discussion list for GROMACS users 


Hi all

I want to run a simulation of one tert-pentane molecule inside one
organic cage molecule. I have the topology files for both
molecules but
obviously I only need one .top file as the input for grompp.
What is the
most convenient way of having both topologies in one file?

Many Thanks

Gavin
--
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 




--


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] topology file for two molecules

2010-09-28 Thread Gavin Melaugh
Hi Justin

Cheers. I have attached my pentane itp file. The cage itp is formatted
the same. In the examples there never seems to be [atom types] included.
I have included it. Does it seem O.K?

Gavin

Justin A. Lemkul wrote:
>
>
> Gavin Melaugh wrote:
>> Hi Mark.
>>
>> I have looked this up, but I am a little uncertain. Does one of the
>> topologies have to be included within the other as an .itp file. If so
>
> No.  You #include .itp files to simply add molecules within a .top; it
> would be very rare to #include an .itp within an .itp (although it can
> be done).  For your case, your .top is simple:
>
> #include (force field)
>
> #include "cage.itp"
> #include "pentane.itp"
>
> [ system ]
> My system
>
> [ molecules ]
> Cage1
> Pentane 1
>
>> do I only define the [defaults directive] in the .top file i.e. I take
>> it I don't include this directive in the .itp file as well.
>>
>
> The [defaults] directive is provided by whatever force field you
> #include.  If you're using some custom parameters, then it needs to be
> defined once (and only once) at the top of your .top file.  Otherwise,
> you'll get fatal errors from grompp.
>
> -Justin
>
>> Cheers
>>
>> Gavin
>>
>> Mark Abraham wrote:
>>> Look up the include file mechanism on the webpage.
>>>
>>> Mark
>>>
>>> - Original Message -
>>> From: Gavin Melaugh 
>>> Date: Tuesday, September 28, 2010 19:04
>>> Subject: [gmx-users] topology file for two molecules
>>> To: Discussion list for GROMACS users 
>>>
>>>> Hi all
>>>>
>>>> I want to run a simulation of one tert-pentane molecule inside one
>>>> organic cage molecule. I have the topology files for both
>>>> molecules but
>>>> obviously I only need one .top file as the input for grompp.
>>>> What is the
>>>> most convenient way of having both topologies in one file?
>>>>
>>>> Many Thanks
>>>>
>>>> Gavin
>>>> -- 
>>>> 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 
>>
>

[atomtypes]
;type mass   charge  ptype sigma(nm) epsilon(kjmol-1)
   CT 12.011000 -0.18   A  0.35  0.276144
   HC  1.008000  0.06   A  0.25  0.125520

;Molecular level
[moleculetype]
;   name nrexcl
tertpentane  3

[atoms]
; atomnr  type   resnr  residue   namecgnr charge mass
1  CT   1PENT  CT   1  0. 12.0110
2  CT   1PENT  CT   1 -0.1800 12.0110
3  CT   1PENT  CT   1 -0.1800 12.0110
4  CT   1PENT  CT   1 -0.1800 12.0110
5  CT   1PENT  CT   1 -0.1800 12.0110
6  HC   1PENT  HC   1  0.0600  1.0080
7  HC   1PENT  HC   1  0.0600  1.0080
8  HC   1PENT  HC   1  0.0600  1.0080
9  HC   1PENT  HC   1  0.0600  1.0080
   10  HC   1PENT  HC   1  0.0600  1.0080
   11  HC   1PENT  HC   1  0.0600  1.0080
   12  HC   1PENT  HC   1  0.0600  1.0080
   13  HC   1PENT  HC   1  0.0600  1.0080
   14  HC   1PENT  HC   1  0.0600  1.0080
   15  HC   1PENT  HC   1  0.0600  1.0080
   16  HC   1PENT  HC   1  0.0600  1.0080
   17  HC   1PENT  HC   1  0.0600  1.0080

[bonds]
;  aiaj   func   b0(nm)   kb(kjmol-1 nm-2)
1 2 1 0.1529  224262.4
1 3 1 0.1529  224262.4
1 4 1 0.1529  224262.4
1 5 1 0.1529  224262.4
2 6 1 0.1090  284512.0
2 7 1 0.1090  284512.0
2 8 1 0.1090  284512.0
3 9 1 0.1090  284512.0
310 1 0.1090  284512.0
311 1 0.1090  284512.0
412 1 0.1090  284512.0
413 1 0.1090  284512.0
414 

Re: [gmx-users] topology file for two molecules

2010-09-28 Thread Justin A. Lemkul



Gavin Melaugh wrote:

Hi Mark.

I have looked this up, but I am a little uncertain. Does one of the
topologies have to be included within the other as an .itp file. If so


No.  You #include .itp files to simply add molecules within a .top; it would be 
very rare to #include an .itp within an .itp (although it can be done).  For 
your case, your .top is simple:


#include (force field)

#include "cage.itp"
#include "pentane.itp"

[ system ]
My system

[ molecules ]
Cage1
Pentane 1


do I only define the [defaults directive] in the .top file i.e. I take
it I don't include this directive in the .itp file as well.



The [defaults] directive is provided by whatever force field you #include.  If 
you're using some custom parameters, then it needs to be defined once (and only 
once) at the top of your .top file.  Otherwise, you'll get fatal errors from grompp.


-Justin


Cheers

Gavin

Mark Abraham wrote:

Look up the include file mechanism on the webpage.

Mark

- Original Message -
From: Gavin Melaugh 
Date: Tuesday, September 28, 2010 19:04
Subject: [gmx-users] topology file for two molecules
To: Discussion list for GROMACS users 


Hi all

I want to run a simulation of one tert-pentane molecule inside one
organic cage molecule. I have the topology files for both
molecules but
obviously I only need one .top file as the input for grompp.
What is the
most convenient way of having both topologies in one file?

Many Thanks

Gavin
--
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 




--


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] topology file for two molecules

2010-09-28 Thread Gavin Melaugh
Hi Mark.

I have looked this up, but I am a little uncertain. Does one of the
topologies have to be included within the other as an .itp file. If so
do I only define the [defaults directive] in the .top file i.e. I take
it I don't include this directive in the .itp file as well.

Cheers

Gavin

Mark Abraham wrote:
> Look up the include file mechanism on the webpage.
>
> Mark
>
> - Original Message -
> From: Gavin Melaugh 
> Date: Tuesday, September 28, 2010 19:04
> Subject: [gmx-users] topology file for two molecules
> To: Discussion list for GROMACS users 
>
> > Hi all
> >
> > I want to run a simulation of one tert-pentane molecule inside one
> > organic cage molecule. I have the topology files for both
> > molecules but
> > obviously I only need one .top file as the input for grompp.
> > What is the
> > most convenient way of having both topologies in one file?
> >
> > Many Thanks
> >
> > Gavin
> > --
> > 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] topology file for two molecules

2010-09-28 Thread Mark Abraham
Look up the include file mechanism on the webpage.

Mark

- Original Message -
From: Gavin Melaugh 
Date: Tuesday, September 28, 2010 19:04
Subject: [gmx-users] topology file for two molecules
To: Discussion list for GROMACS users 

> Hi all
> 
> I want to run a simulation of one tert-pentane molecule inside one
> organic cage molecule. I have the topology files for both 
> molecules but
> obviously I only need one .top file as the input for grompp. 
> What is the
> most convenient way of having both topologies in one file?
> 
> Many Thanks
> 
> Gavin
> -- 
> 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

[gmx-users] topology file for two molecules

2010-09-28 Thread Gavin Melaugh
Hi all

I want to run a simulation of one tert-pentane molecule inside one
organic cage molecule. I have the topology files for both molecules but
obviously I only need one .top file as the input for grompp. What is the
most convenient way of having both topologies in one file?

Many Thanks

Gavin
-- 
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] topology for a cyclic peptide

2010-07-25 Thread Justin A. Lemkul



Rebeca García Fandiño wrote:

Hi,
I am trying to get  the topology of a cyclic peptide, but when I try to 
do pdb2pgx I get automatically the correction to a terminal one.
I have looked in the Gromacs list and I only have found a entry about 
it,  http://osdir.com/ml/science.biology.gromacs.user/2006-08/msg00297.html

but I cannot find a solution from there.
Does anyone have any idea about it, please?


You'll have to do a bit of manual work yourself.  I don't entirely understand 
the approach quoted above, but running pdb2gxm -ter (choosing "None" for both 
termini), then adding all necessary bonded information (bonds, angles, 
dihedrals, and impropers) manually should do the trick.


-Justin


Thanks a lot for your help.
Best wishes,

Rebeca García
Santiago de Compostela University
Spain



Hotmail prepara novedades y sorpresas en breve, ¡estate atento! 





--


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/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/mailing_lists/users.php


[gmx-users] Topology file for phosphoserine (SEP) for OPLS force field

2010-07-22 Thread sumanth kumar
Dear users,

Does anyone have an idea about where I could get the *phosphoserine (SEP)
itp file???* for *opls aa forcefield?*

 Can some one give me information to be added in *ffoplsaa.rtp*???

Thank you very much in advance
sumanth

-- 
*SUMANTH KUMAR.M*
IV B.Tech Bioinformatics,
Centre for Plant Molecular Biology,
Tamilnadu Agricultural University,
Coimbatore-641003
-- 
gmx-users mailing listgmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/search before posting!
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/mailing_lists/users.php

[gmx-users] topology builiding

2010-07-20 Thread pankhuri arora
Hi There,
I am using x2top to build topology for a molecule using oplsaa forcefield
from a pdb file. I have added a few entries in the n2t file but i have got a
wrong topology as some additional bonds are made in the molecule. As far as
I know x2top ignores the connect information of the atoms. Is it possible to
modify the x2top code so that it reads the connect information and then
build bonds . Also, if I make some changes in the x2top code, do I need to
compile it the whole gromacs package again or whether I can compile it
separately?

Thanks,
Pankhuri
-- 
gmx-users mailing listgmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/search before posting!
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/mailing_lists/users.php

[gmx-users] topology for a cyclic peptide

2010-07-19 Thread Rebeca García Fandiño

Hi,

I am trying to get  the topology of a cyclic peptide, but when I try to 
do pdb2pgx I get automatically the correction to a terminal one.

I have looked in the Gromacs list and I only have found a entry about 
it,  
http://osdir.com/ml/science.biology.gromacs.user/2006-08/msg00297.html

but I cannot find a solution from there.

Does anyone have any idea about it, please?

Thanks a lot for your help.

Best wishes,



Rebeca García 

Santiago de Compostela University

Spain

  
_
No has visto nada como el nuevo Messenger, ¡te sorprenderá!
http://explore.live.com/windows-live-messenger-- 
gmx-users mailing listgmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/search before posting!
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/mailing_lists/users.php

[gmx-users] Topology Generation

2010-06-04 Thread sonali dhindwal
Hello All,

I have to do simulation on one of the protein which contains one dioxy group, 
as gromacs do not have information to generate toplogy for it, I tried it with 
PRODRG, but is giving error:
ERRDRG> PRODRG does not deal with mono/di-atomic molecules.
PRODRG> Program terminated unsuccessfully, sorry!So if there is any server by 
which I could do the same ? please help

Thanks and regards.

--
Sonali Dhindwal

-- 
gmx-users mailing listgmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/search before posting!
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/mailing_lists/users.php

Re: [gmx-users] Topology files for OPLSAA and RESP charges

2010-05-26 Thread Nicolas Sapay

Hello,

I just read your message... As far as I know, the RESP charge method is 
mainly employed forthe AMBER force field and associated (GLYCAM, GAFF, 
etc). I would check first that RESP is compatible with the method 
employed to develop OPLS-AA. Note that the AmberTools provide some 
scripts to calculate RESP charges. I have also used the RED server 
(http://q4md-forcefieldtools.org/REDS/), which is quite convenient when 
you want to take into account different conformations of a same 
molecule. In all cases, you will need to provide to those tools the 
molecular electrostatic potential of your molecule. None of this tools 
can generate a Gromacs topology file, however.


Regards,
Nicolas

Tanos C. C. Franca a écrit :

   Dear GROMACS users,
   I'm trying to start running jobs using RESP charges and OPLSAA 
force field but I am facing problems to generate the RESP charges and 
the topology files for the OPLSAA force field. Does someone knows a 
software to calculate the RESP charges for the ligands and/or, also, 
to generate the .itp files for OPLSAA ?

   With the best regards,
   Tanos C. C. Franca.
<>-- 
gmx-users mailing listgmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/search before posting!
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/mailing_lists/users.php

Re: [gmx-users] Topology files for OPLSAA and RESP charges

2010-05-23 Thread Mark Abraham
- Original Message -
From: "Tanos C. C. Franca" 
Date: Monday, May 24, 2010 1:03
Subject: [gmx-users] Topology files for OPLSAA and RESP charges
To: gmx-users@gromacs.org

>    Dear GROMACS users,
>    I'm trying to start running jobs using RESP charges 
> and OPLSAA force field but I am facing problems to generate the 
> RESP charges and the topology files for the OPLSAA force field. 
> Does someone knows a software to calculate the RESP charges for 
> the ligands and/or, also, to generate the .itp files for OPLSAA ?

The general advice is here: 
http://www.gromacs.org/Documentation/How-tos/Parametrization. I am not aware of 
an OPLSAA-specific helper tool. Be wary of presupposing that use of RESP is 
correct for OPLSAA.

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/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/mailing_lists/users.php


[gmx-users] Topology files for OPLSAA and RESP charges

2010-05-23 Thread Tanos C. C. Franca

   Dear GROMACS users,
   I'm trying to start running jobs using RESP charges and OPLSAA force 
field but I am facing problems to generate the RESP charges and the 
topology files for the OPLSAA force field. Does someone knows a software 
to calculate the RESP charges for the ligands and/or, also, to generate 
the .itp files for OPLSAA ?

   With the best regards,
   Tanos C. C. Franca.
--
gmx-users mailing listgmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/search before posting!
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.

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


Re: [gmx-users] Topology and force fields

2010-04-05 Thread Justin A. Lemkul



Gareth Tribello wrote:

Hello

I think what you have been told is wrong.  You need to include both 
"ffcharmm27.itp" and "tip3p.itp" as the atomtypes in your tip3p.itp will 


That is exactly what was in the original topology, no problem there.  There are 
no [atomtypes] in tip3p.itp; these are defined by ff*nb.itp.


otherwise not be defined.  There is an easy way to check though - take 
out the include ffcharmm27.itp and try to run grompp on your new 
topology.  My bet would be that it won't work and, if you are using 
rigid molecules, it will tell you that it finds a settles section in the 
wrong place.




More likely grompp will raise a fatal error that the atomtypes are not found, 
before anything else happens :)


You could, in theory, run a simulation of just water by specifying a all the 
directives listed at the parameter level (Table 5.3 and section 5.7.1 in the 
manual) in the .top file, then #include "tip3p.itp" - that's all that's really 
necessary.  Using the #include mechanism with an existing force field solves 
this problem for you.  The combination rules of the Gromos, OPLS, and CHARMM 
force fields are all different (manual section 5.3.3, equations 5.3 and 5.4), 
but I don't know offhand what effect that might have on your results.


TIP3P itself is not really a force field, it is a set of parameters that model a 
single molecule.  The interaction function (the other part of a force field) is 
defined by the functional forms defined by ff*.itp in the parameter-level 
directives.  Otherwise, the atom types and nonbonded parameters for a water 
model should be uniform, that is, independent of the force field chosen (with 
the only modifications being conversion between sigma/epsilon and C6/C12, per 
force field requirements).  Your treatment of other .mdp parameters will affect 
your results, as well, but many of these should be prescribed in the derivation 
of the parameters or in subsequent modifications and extensions.


-Justin


ciao
Gareth

On Mon, Apr 5, 2010 at 7:26 PM, Lum Nforbi > wrote:


Hello all,

Below are the contents of a water1.top file that I used to
simulate 2000 molecules of TIP3P water (NPT).  I made this topology
file according to the pattern on page 114 of the manual (topology
for Urea in water). My diffusion coefficient, viscosity and radial
distribution functions matched those for tip3p water in literature.
Recently I was told that I did not need to use another force
field in the topology like ffcharmm27.itp since tip3p.itp is already
a forcefield on its own. Could this be true? I just wanted to verify
the authenticity of my simulation. If the structure of the topology
below is not correct, please let me know.

; The force field files to be included
#include "ffcharmm27.itp"

; Include TIP3P water topology
#include "tip3p.itp"

[ system ]
Pure Water

[ molecules ]
;molecule name number
SOL  1936
SOL64

I appreciate your answers.

Lum

--
gmx-users mailing listgmx-users@gromacs.org

http://lists.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/search before
posting!
Please don't post (un)subscribe requests to the list. Use the
www interface or send it to gmx-users-requ...@gromacs.org
.
Can't post? Read http://www.gromacs.org/mailing_lists/users.php




--


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/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/mailing_lists/users.php


Re: [gmx-users] Topology and force fields

2010-04-05 Thread Gareth Tribello
Hello

I think what you have been told is wrong.  You need to include both
"ffcharmm27.itp" and "tip3p.itp" as the atomtypes in your tip3p.itp will
otherwise not be defined.  There is an easy way to check though - take out
the include ffcharmm27.itp and try to run grompp on your new topology.  My
bet would be that it won't work and, if you are using rigid molecules, it
will tell you that it finds a settles section in the wrong place.

ciao
Gareth

On Mon, Apr 5, 2010 at 7:26 PM, Lum Nforbi  wrote:

> Hello all,
>
> Below are the contents of a water1.top file that I used to simulate
> 2000 molecules of TIP3P water (NPT).  I made this topology file according to
> the pattern on page 114 of the manual (topology for Urea in water). My
> diffusion coefficient, viscosity and radial distribution functions matched
> those for tip3p water in literature.
> Recently I was told that I did not need to use another force field
> in the topology like ffcharmm27.itp since tip3p.itp is already a forcefield
> on its own. Could this be true? I just wanted to verify the authenticity of
> my simulation. If the structure of the topology below is not correct, please
> let me know.
>
> ; The force field files to be included
> #include "ffcharmm27.itp"
>
> ; Include TIP3P water topology
> #include "tip3p.itp"
>
> [ system ]
> Pure Water
>
> [ molecules ]
> ;molecule name number
> SOL  1936
> SOL64
>
> I appreciate your answers.
>
> Lum
>
> --
> gmx-users mailing listgmx-users@gromacs.org
> http://lists.gromacs.org/mailman/listinfo/gmx-users
> Please search the archive at http://www.gromacs.org/search before posting!
> Please don't post (un)subscribe requests to the list. Use the
> www interface or send it to gmx-users-requ...@gromacs.org.
> Can't post? Read http://www.gromacs.org/mailing_lists/users.php
>
-- 
gmx-users mailing listgmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/search before posting!
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/mailing_lists/users.php

[gmx-users] Topology and force fields

2010-04-05 Thread Lum Nforbi
Hello all,

Below are the contents of a water1.top file that I used to simulate
2000 molecules of TIP3P water (NPT).  I made this topology file according to
the pattern on page 114 of the manual (topology for Urea in water). My
diffusion coefficient, viscosity and radial distribution functions matched
those for tip3p water in literature.
Recently I was told that I did not need to use another force field
in the topology like ffcharmm27.itp since tip3p.itp is already a forcefield
on its own. Could this be true? I just wanted to verify the authenticity of
my simulation. If the structure of the topology below is not correct, please
let me know.

; The force field files to be included
#include "ffcharmm27.itp"

; Include TIP3P water topology
#include "tip3p.itp"

[ system ]
Pure Water

[ molecules ]
;molecule name number
SOL  1936
SOL64

I appreciate your answers.

Lum
-- 
gmx-users mailing listgmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/search before posting!
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/mailing_lists/users.php

  1   2   >