[gmx-users] Selecting force field for simulation in vacuum

Berk Hess gmx3 at hotmail.com
Thu Dec 4 13:36:57 CET 2008


Hi,

It depends on what you mean with vacuum simulation.
Do you mean really simulating how a protein would behave in vacuum?
Or do you want some type of implicit solvent?

OPLS has been parametrized on the properties of pure liquids of small compounds.
For Gromos at least the liquid properties of alkanes have been checked.
So both force fields would work well for real vacuum simulations.

Berk


> Date: Thu, 4 Dec 2008 07:17:12 -0500
> From: jalemkul at vt.edu
> To: nazoia at det.uminho.pt; gmx-users at gromacs.org
> Subject: Re: [gmx-users] Selecting force field for simulation in vacuum
> CC: 
> 
> GROMOS and OPLS were derived for condensed-phase properties, so none of those 
> will be useful.  I know nothing about ENCAD, however, but I do not often see it 
> in the literature.
> 
> -Justin
> 
> Nuno Azoia wrote:
> > Hello everyone!
> > 
> > I'm having some problems for selecting force field for a simulation in
> > vacuum. I'm using version 4.0.2.
> > According to pdb2gmx man page (pdb2gmx -h)
> > 
> >         DESCRIPTION
> >         -----------
> >         This program reads a pdb file, reads some database files, adds
> >         hydrogens to
> >         the molecules and generates coordinates in Gromacs (Gromos)
> >         format and a
> >         topology in Gromacs format. These files can subsequently be
> >         processed to
> >         generate a run input file. 
> >         
> >         The force fields in the distribution are currently:
> >         
> >         oplsaa OPLS-AA/L all-atom force field (2001 aminoacid dihedrals)
> >         G43b1  GROMOS96 43b1 Vacuum Forcefield 
> >         G43a1  GROMOS96 43a1 Forcefield 
> >         G43a2  GROMOS96 43a2 Forcefield (improved alkane dihedrals)
> >         G45a3  GROMOS96 45a3 Forcefield 
> >         G53a5  GROMOS96 53a5 Forcefield 
> >         G53a6  GROMOS96 53a6 Forcefield 
> >         gmx    Gromacs Forcefield (a modified GROMOS87, see manual)
> >         encads Encad all-atom force field, using scaled-down vacuum
> >         charges
> >         encadv Encad all-atom force field, using full solvent charges
> > 
> > 
> > Ok. I would like to use G43b1, but when I run
> > 
> > pdb2gmx -f protein.pdb -o protein.gro -p protein.top -n system.ndx -i
> > posre.itp  -inter -ignh
> > 
> > I just get this options, with no possibility for choosing G43b1:
> > 
> >         Select the Force Field:
> >          0: GROMOS96 43a1 force field 
> >          1: GROMOS96 43a2 force field (improved alkane dihedrals)
> >          2: GROMOS96 45a3 force field (Schuler JCC 2001 22 1205)
> >          3: GROMOS96 53a5 force field (JCC 2004 vol 25 pag 1656) 
> >          4: GROMOS96 53a6 force field (JCC 2004 vol 25 pag 1656) 
> >          5: OPLS-AA/L all-atom force field (2001 aminoacid dihedrals)
> >          6: [DEPRECATED] Gromacs force field (see manual)
> >          7: [DEPRECATED] Gromacs force field with hydrogens for NMR
> >          8: Encad all-atom force field, using scaled-down vacuum charges
> >          9: Encad all-atom force field, using full solvent charges    
> > 
> > 
> > I search on the mailing list for some solution and I get this message
> > from Berk Hess
> > (http://www.gromacs.org/component/option,com_wrapper/Itemid,165/)
> > 
> >         The GROMOS96 G43b1 "vacuum" force field has been removed.
> >         This confusing "vacuum" naming of this force field could make people
> >         think that this force field was optimized for real vacuum simulations,
> >         whereas it is actually G43a1 with neutralized charged residues.
> >         This force field has never been used seriously.
> > 
> > Now I'm confused! Ok, there is no more G43b1.
> > So I would like to ask witch force field I can use for simulations in vaccum? If I search on the mailing list I get the advice of using G43b1.
> > Can anyone help me on this.
> > 
> > And just another small thing. Maybe it would be good to remove the reference to G43b1 on the pdb2gmx man page.
> > 
> > Thank you in advance.
> > 
> > Nuno Azoia
> > 
> > --
> > 
> > Nuno Gonçalo Azoia Lopes
> > 
> > Laboratório de Investigação em Acabamento
> > Departamento de Engenharia Textil
> > Universidade do Minho
> > Campus de Azurém
> > 4800-058 Guimarães
> > Portugal
> > 
> > Tel: +351 253 510 280 - Ext: 517 289
> > Fax: +351 253 510 293
> > 
> > Mobile: +351 965 382 487
> > E-mail: nazoia at det.uminho.pt
> > 
> > _______________________________________________
> > gmx-users mailing list    gmx-users at gromacs.org
> > http://www.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-request at gromacs.org.
> > Can't post? Read http://www.gromacs.org/mailing_lists/users.php
> > 
> 
> -- 
> ========================================
> 
> Justin A. Lemkul
> Graduate Research Assistant
> 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 at gromacs.org
> http://www.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-request at gromacs.org.
> Can't post? Read http://www.gromacs.org/mailing_lists/users.php

_________________________________________________________________
Express yourself instantly with MSN Messenger! Download today it's FREE!
http://messenger.msn.click-url.com/go/onm00200471ave/direct/01/
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://maillist.sys.kth.se/pipermail/gromacs.org_gmx-users/attachments/20081204/7173ca88/attachment.html>


More information about the gromacs.org_gmx-users mailing list