AnteChamber PYthon Parser interfacE
A tool based in Python to use Antechamber to generate topologies for chemical compounds and to interface with others python applications like CCPN or ARIA.
acpype is pronounced as "ace + pipe"
Topologies files to be generated so far: CNS/XPLOR, GROMACS, CHARMM and AMBER.
NB: Topologies generated by acpype/Antechamber are based on General Amber Force Field (GAFF) and should be used only with compatible forcefields like AMBER and its variant.
Several flavours of AMBER FF are ported already for GROMACS (see ffamber: http://ffamber.cnsm.csulb.edu/) as well as to XPLOR/CNS (see xplor-nih: http://ambermd.org/xplor-nih.html) and CHARMM.
This code is released under GNU General Public License V3.
It was inspired by:
-
amb2gmx.pl (Eric Sorin, David Mobley and John Chodera) and depends on Antechamber and Openbabel
-
YASARA Autosmiles: http://www.yasara.org/autosmiles.htm (Elmar Krieger)
-
topolbuild (Bruce Ray)
-
xplo2d (G.J. Kleywegt)
For Non-uniform 1-4 scale factor conversion (e.g. if using GLYCAM06), please cite:
Austen Bernardi, Roland Faller, Dirk Reith, Karl N. Kirschner, "Conversion of GLYCAM06 Force-Field Parameters from AMBER to GROMACS Topologies" submitted for review, April 2018.
For Antechamber, please cite:
- Wang, J., Wang, W., Kollman P. A.; Case, D. A. "Automatic atom type and bond type perception in molecular mechanical calculations". Journal of Molecular Graphics and Modelling , 25, 2006, 247260.
- Wang, J., Wolf, R. M.; Caldwell, J. W.; Kollman, P. A.; Case, D. A. "Development and testing of a general AMBER force field". Journal of Computational Chemistry, 25, 2004, 1157-1174.
If you use this code, I am glad if you cite:
SOUSA DA SILVA, A. W. & VRANKEN, W. F. ACPYPE - AnteChamber PYthon Parser interfacE. BMC Research Notes 2012, 5:367 doi:10.1186/1756-0500-5-367 http://www.biomedcentral.com/1756-0500/5/367
BATISTA, P. R.; WILTER, A.; DURHAM, E. H. A. B. & PASCUTTI, P. G. Molecular Dynamics Simulations Applied to the Study of Subtypes of HIV-1 Protease. Cell Biochemistry and Biophysics, 44, 395-404, 2006.
Alan Wilter Sousa da Silva, D.Sc. Bioinformatician, UniProt, EMBL-EBI Hinxton, Cambridge CB10 1SD, UK. http://www.ebi.ac.uk/~awilter
alanwilter at gmail dot com
To run acpype with its all functionalities, you need ANTECHAMBER from package AmberTools and Open Babel if your input files are of PDB format.
However, if one wants acpype just to emulate amb2gmx.pl, one needs nothing at all but Python.
At the moment, is possible to download acpype using git:
git clone https://github.com/alanwilter/acpype.git
At folder acpype/test, type:
../acpype.py -i FFF.pdb
It'll create a folder called FFF.acpype, and inside it one may find topology files for GROMACS and CNS/XPLOR.
To get help and more information, type:
../acpype.py -h
At folder acpype, type:
ln -s $PWD/acpype.py /usr/local/bin/acpype
And re-login or start another shell session.
Yet, it is possible to install using Anaconda package Python 3.6:
conda install -c acpype acpype
GROMACS < v.5.0
cd FFF.acpype/
grompp -c FFF_GMX.gro -p FFF_GMX.top -f em.mdp -o em.tpr
mdrun -v -deffnm em
# And if you have VMD
vmd em.gro em.trr
GROMACS > v.5.0
cd FFF.acpype/
gmx grompp -c FFF_GMX.gro -p FFF_GMX.top -f em.mdp -o em.tpr
gmx mdrun -v -deffnm em
# And if you have VMD
vmd em.gro em.trr
GROMACS < v.5.0
grompp -c em.gro -p FFF_GMX.top -f md.mdp -o md.tpr
mdrun -v -deffnm md
vmd md.gro md.trr
GROMACS > v.5.0
gmx grompp -c em.gro -p FFF_GMX.top -f md.mdp -o md.tpr
gmx mdrun -v -deffnm md
vmd md.gro md.trr
GROMACS < v.5.0
grompp -c FFF_GMX.gro -p FFF_GMX.top -f em.mdp -o em.tpr
om-mpirun -n 2 mdrun_mpi -v -deffnm em
grompp -c em.gro -p FFF_GMX.top -f md.mdp -o md.tpr
om-mpirun -n 2 mdrun_mpi -v -deffnm md
vmd md.gro md.trr
GROMACS > v.5.0
gmx grompp -c FFF_GMX.gro -p FFF_GMX.top -f em.mdp -o em.tpr
gmx mdrun -ntmpi 2 -v -deffnm em
For any given prmtop and inpcrd files (outputs from AMBER LEaP), type:
acpype -p FFF_AC.prmtop -x FFF_AC.inpcrd
The output files FFF_GMX.gro
and FFF_GMX.top
will be generated at the same
folder of the input files.
At folder FFF.acpype, type:
cns < FFF_CNS.inp
- see [TutorialNAMD]