Amber

Download and general information: http://ambermd.org

The Amber 10 molecular dynamics package has been installed on the SciNet GPC cluster, together with the amber tools package. In order to have access to this package, you need to "module load intel intelmpi fftw amber10" (best included in your .bashrc). This installation includes both serial and parallel (MPI) versions of amber. It was compiled with the Intel compilers, and uses the Intel MPI libraries. The pmemd program was compiled with fftw3, Intel compilers, and Intel MPI.

A typical workflow with amber, as exemplified in the Amber 10 Users Manual, is:

Step 1. Generate some starting coordinates. The first step is to obtain starting coordinates. We begin with the bovine pancreatic trypsin inhibitor, and consider the file 6pti.pdb, exactly as distributed by the Protein Data Bank. This file (as with most PDB files) needs some editing before it can be used by Amber. First, alternate conformations are provided for residues 39 and 50, so we need to figure out which one we want. For this example, we choose the "A" conformation, and manually edit the file to remove the alternate conformers. Second, coordinates are provided for a phosphate group and a variety of water molecules. These are not needed for the calculation we are pursuing here, so we also edit the file to remove these. Third, the cysteine residues are involved in disulfide bonds, and need to have their residue names changed in an editor from CYS to CYX to reflect this. Finally, since we removed the phosphate groups, some of the CONECT records now refer to non-existent atoms; if you are not sure that the CONECT records are all correct then it may be safest to remove all of them, as we do for this example. Let’s call this modified file 6pti.mod.pdb. Although Amber tries hard to understand pdb-format files, it is typical to have to do some manual editing before proceeding. A general prescription is: "keep running the loadPdb step in LEaP (see step 2, below), and editing the pdb file, until there are no error messages."

Step 2. Run LEaP to generate the parameter and topology file. This is a fairly straightforward exercise in loading in the pdb file, adding the disulfide cross links, and saving the resulting files. Typing the following commands should work in either tleap or xleap: source leaprc.ff03 bpti = loadPdb 6pti.mod.pdb bond bpti.5.SG bpti.55.SG bond bpti.14.SG bpti.38.SG bond bpti.30.SG bpti.51.SG saveAmberParm bpti prmtop prmcrd quit

Step 3. Perform some minimization. Use this script (I call it "bpti.pbs"), appropriate for the SciNet systems, to run parallel sander (sander.MPI) in the batch queue of the GPC machine:

<source lang="bash">

1. !/bin/bash
2. PBS -l nodes=1:ppn=8,walltime=1:30:00,os=centos53computeA
3. PBS -N sander

1. the following line should be in the user's .bashrc
2. module add intel openmpi fftw amber10

1. change to the submission directory

cd $PBS_O_WORKDIR

1. Running minimization for BPTI

1. Create the input control file

cat << eof > min.in

1. 200 steps of minimization, generalized Born solvent model

&cntrl maxcyc=200, imin=1, cut=12.0, igb=1, ntb=0, ntpr=10, / eof

1. Run the parallel version of sander, using all 8 cores in the node

mpirun -r ssh -np 8 $AMBERHOME/bin/sander.MPI -i min.in -o 6pti.min1.out -c prmcrd -r 6pti.min1.xyz

</source>