# load forcefields
source leaprc.ff99SB
source leaprc.gaff

# load additional prepins, e.g. for ligands and ions 

# load structure
COM = loadpdb "1ahp.pdb"

# vacuum topologies and coordinates

# solvate
solvatebox COM TIP3PBOX  5 
savepdb COM prot_1.pdb
saveamberparm COM prot_1.top prot_1.crd

quit
