使用的代码:
source leaprc.protein.ff14SB
pro = loadpdb protein_removeH.pdb
check pro
source leaprc.gaff
loadamberparams ligand.frcmod
lig = loadmol2 ligand_resp.mol2
check lig
com = combine {pro lig}
check com
source leaprc.water.tip3p
solvatebox com TIP3PBOX 10.0
charge com
addIons2 com Cl- 0
check com
saveamberparm com com.prmtop com.inpcrd