(Note: These tutorials are meant to provide
illustrative examples of how to use the AMBER software suite to carry out
simulations that can be run on a simple workstation in a reasonable period of
time. They do not necessarily provide the optimal choice of parameters or
methods for the particular application area.)
Copyright Ross Walker 2006
TUTORIAL 5
Using Antechamber to Create Leap Input Files for Simulating
Sustiva
using the General Amber Force Field
By Ross Walker
In this tutorial we will make use of the Antechamber tools, which ship with AMBER 8, in order to create an input file that can be read by Leap so that we can create prmtop and inpcrd files for simulations of organic molecules.
Antechamber is designed to be used with the "general AMBER force field (GAFF)1". This force field has been specifically designed to cover most pharmaceutical molecules and is compatible with the traditional AMBER force fields in such a way that the two can be mixed during a simulation. Like the traditional AMBER force fields, GAFF uses a simple harmonic function form for bonds and angles but unlike the traditional protein and DNA orientated AMBER force fields the atom types used in GAFF are much more general such that they cover most of the organic chemical space. The current implementation of the GAFF force field consists of 33 basic atom types and 22 special atom types. The charge methods used can be HF/6-31G* RESP or AM1-BCC2.
By design, GAFF, is a complete force field (so that missing parameters rarely occur), it covers almost all the organic chemical space that is made up of C, N, O, S, P, H, F, Cl, Br and I. Transition metals such as gold (Au), silver (Ag), and tungsten (W) are not currently covered. Moreover, since GAFF is totally compatible with the AMBER macromolecular force fields it should prove to be an useful molecular mechanical tool for rational drug design. Especially in binding free energy calculations and molecular docking studies.
The Antechamber tool set is designed to allow the rapid generation of topology files for use with the AMBER simulation programs. Unlike the previous tutorial where we manually assigned our atom types and parameters here we will allow antechamber to do this automatically for us using GAFF. With Antechamber, one may solve the following problems:
Automatically identify bond and atom types
Judge atomic equivalence
Generate residue topology files
Find missing force field parameters and supply reasonable suggestions
You should note, however, that Antechamber is not a replacement for due dilligence. You should always closely examine the atom types that Antechamber assigns and verify to yourself that the choices are reasonable. You should never use scientific software in a "Black Box" approach!
CLICK HERE TO GO TO THE NEXT SECTION
(Note: These tutorials are meant to provide
illustrative examples of how to use the AMBER software suite to carry out
simulations that can be run on a simple workstation in a reasonable period of
time. They do not necessarily provide the optimal choice of parameters or
methods for the particular application area.)
Copyright Ross Walker 2006