Antechamber — Force Field Parameterization

When to Use

• User needs GAFF/GAFF2 atom types and parameters for a small molecule
• User wants AM1-BCC or RESP partial atomic charges
• User is preparing a molecule for classical MD with AMBER or LAMMPS
• User needs to generate topology (prmtop) and coordinate (inpcrd) files

Prerequisites

1. AmberTools installed (antechamber --help)
2. Input molecule structure (mol2, pdb, or Gaussian output for RESP)
3. For RESP charges: Gaussian output with pop=mk iop(6/33=2) ESP data

Workflow Steps

1. Generate GAFF parameters with AM1-BCC charges

catgo_workflow_engine(action="add_task", params={
  "workflow_id": "wf_xxx",
  "task_type": "shell",
  "name": "antechamber_param",
  "command": "bash parameterize.sh",
  "input_files": {
    "parameterize.sh": "<script content>",
    "molecule.pdb": "<molecule structure>"
  },
  "system_name": "ligand_param"
})

Parameterization Script — AM1-BCC Charges

#!/bin/bash
set -e

# Step 1: Assign GAFF2 atom types and AM1-BCC charges
antechamber -i molecule.pdb -fi pdb \
            -o molecule.mol2 -fo mol2 \
            -c bcc -at gaff2 \
            -nc 0 -m 1 \
            -rn LIG

# Step 2: Check for missing parameters
parmchk2 -i molecule.mol2 -f mol2 \
         -o molecule.frcmod -a Y

# Step 3: Build topology with tleap
cat > tleap.in << 'EOF'
source leaprc.gaff2
LIG = loadmol2 molecule.mol2
loadamberparams molecule.frcmod
check LIG
saveamberparm LIG molecule.prmtop molecule.inpcrd
savemol2 LIG molecule_leap.mol2 1
quit
EOF

tleap -f tleap.in > tleap.log 2>&1
echo "Generated: molecule.prmtop molecule.inpcrd"

Parameterization Script — RESP Charges

For higher-quality charges, use RESP from Gaussian ESP:

#!/bin/bash
set -e

# Step 1: Generate Gaussian input for ESP
antechamber -i molecule.pdb -fi pdb \
            -o molecule.gjf -fo gcrt \
            -gm "%mem=4GB" -gn "%nproc=8" \
            -ge "molecule.gesp" \
            -gk "#HF/6-31G* opt pop=mk iop(6/33=2)"

# Step 2: Run Gaussian (separate task)
# g16 molecule.gjf

# Step 3: Extract RESP charges from Gaussian output
antechamber -i molecule.log -fi gout \
            -o molecule.mol2 -fo mol2 \
            -c resp -at gaff2 \
            -nc 0 -m 1 -rn LIG

# Step 4: Check and build topology (same as AM1-BCC)
parmchk2 -i molecule.mol2 -f mol2 -o molecule.frcmod -a Y

cat > tleap.in << 'EOF'
source leaprc.gaff2
LIG = loadmol2 molecule.mol2
loadamberparams molecule.frcmod
check LIG
saveamberparm LIG molecule.prmtop molecule.inpcrd
quit
EOF

tleap -f tleap.in > tleap.log 2>&1

Solvation with tleap

After parameterization, solvate the molecule:

source leaprc.water.tip3p
LIG = loadmol2 molecule.mol2
loadamberparams molecule.frcmod
solvateBox LIG TIP3PBOX 12.0
addIonsRand LIG Na+ 0 Cl- 0   # Neutralize
saveamberparm LIG solvated.prmtop solvated.inpcrd
quit

Parameter Guidance

Flag	Purpose
-c bcc	AM1-BCC charges (fast, good for most organic molecules)
-c resp	RESP charges (requires Gaussian ESP, higher quality)
-at gaff2	GAFF2 atom types (preferred over gaff)
-nc 0	Net charge of molecule
-m 1	Spin multiplicity
-rn LIG	Residue name (3 chars max)
-fi / -fo	Input/output format (pdb, mol2, gout, gcrt)

Common Pitfalls

1. Wrong net charge — -nc must match the actual molecular charge. Wrong charge gives wrong AM1-BCC charges.
2. Missing parameters — parmchk2 with -a Y fills gaps by analogy but warns. Check molecule.frcmod for ATTN lines.
3. Atom name conflicts — tleap requires unique atom names within a residue. antechamber usually handles this.
4. RESP without ESP — RESP charges require a Gaussian calculation with pop=mk iop(6/33=2). AM1-BCC does not.
5. GAFF vs GAFF2 — use GAFF2 (-at gaff2, leaprc.gaff2) for improved parameters. Do not mix GAFF and GAFF2.
6. Radical species — antechamber struggles with radicals and metals. Consider manual parameterization.
7. Large molecules — AM1 optimization in antechamber can be slow for >100 atoms. Pre-optimize geometry.