Re: [AMBER-Developers] force-field naming (redux)

From: Ross Walker <ross.rosswalker.co.uk>
Date: Sun, 8 Mar 2015 12:40:47 -0700

I would like there to be a 'default' force field that is either:

1) Loaded by leap by default.

2) Loaded with 'source leaprc.default'

This would make things a lot simpler for beginning users. Option 1 might prove troublesome since I don't know how well one can unload parameter sets in leap if you wanted to switch things out so option 2 with documentation might be optimal.

I would suggest it be:

source leaprc.ff14SB
source leaprc.gaff
source leaprc.GLYCAM_06j-1
source leaprc.lipid14
loadamberparams frcmod.ionsjc_tip3p

This would be good since it will also give us early warning in the future of atom type clashes, incompatibilities etc.

Does it make sense to make TIP3P with ionsjc the default water model and ion parameters? - These are probably the most commonly used.

What do people think?

Then we just have a 'quickstart' section as it were in the manual (and we update the tutorials to use this default) and then in the future we can update this with whatever is the latest set of force fields we recommend.

All the best
Ross

> On Mar 8, 2015, at 6:31 AM, David A Case <case.biomaps.rutgers.edu> wrote:
>
> We have not yet addressed the problem of force field naming that was brought
> up at the Amber developers' meeting in February. Some notes:
>
> 1. Section 3.2.1 of the manual still reads as though ff14SB is only for
> proteins. There is one tiny sentence, many paragraphs later, that indicates
> it could also be used for RNA; nothing about whether it might be good for a
> combination of proteins and nucleic acids.
>
> 2. Table 3.1 is not as helpful to users as it should be. It lists of bunch
> of modifications, but gives no indication of *how* to load any of the force
> fields that are listed there.
>
> 3. We talked about having separate leaprc files for proteins, DNA and RNA,
> but Jason has noted potential problems with this. I don't have a suggested
> solution (and I've not really heard one): but we need very soon to come to
> some agreement, and to have appropriate leaprc files for various recommended
> combintations of parameters.
>
> 4. On a related topic: we need a better way for a user to specify water and
> ion parameters. Specifying TIP3P water should automatically load the
> associated ion paramters; same for other supported water models. This should
> probably be implemented by having a leaprc.tip3p that performs these actions.
>
> 5. (I can do this one): the beginning of section 3.1 should probably suggest
> to people that they load a "leap.in" file, which in turn would load *several*
> leaprc files for the components they need. That is, the "standard" would be
> to expect to load at least a protein and a water/ion file; if there are
> ligands, one would usually load leaprc.gaff; if one has lipids
> or nucleic acids or carbohydrates, one would load additional files.
>
> Suggestions are both welcome and needed....thx...dac
>
>
> _______________________________________________
> AMBER-Developers mailing list
> AMBER-Developers.ambermd.org
> http://lists.ambermd.org/mailman/listinfo/amber-developers


_______________________________________________
AMBER-Developers mailing list
AMBER-Developers.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber-developers
Received on Sun Mar 08 2015 - 13:00:02 PDT
Custom Search