$ tleap
>source leaprc.protein.ff14SB
----- Source: /home/wangq/Programs/amber16/dat/leap/cmd/leaprc.protein.ff14SB ----- Source of /home/wangq/Programs/amber16/dat/leap/cmd/leaprc.protein.ff14SB done Log file: ./leap.log Loading parameters: /home/wangq/Programs/amber16/dat/leap/parm/parm10.dat Reading title: PARM99 + frcmod.ff99SB + frcmod.parmbsc0 + OL3 for RNA Loading parameters: /home/wangq/Programs/amber16/dat/leap/parm/frcmod.ff14SB Reading force field modification type file (frcmod) Reading title: ff14SB protein backbone and sidechain parameters Loading library: /home/wangq/Programs/amber16/dat/leap/lib/amino12.lib Loading library: /home/wangq/Programs/amber16/dat/leap/lib/aminoct12.lib Loading library: /home/wangq/Programs/amber16/dat/leap/lib/aminont12.lib
> loadoff mg.off
Loading library: ./mg.off
> mol = loadpdb wt1mg.pdb
Loading PDB file: ./wt1mg.pdb Added missing heavy atom: .R<CGLN 154>.A<OXT 18> total atoms in file: 1189 Leap added 1192 missing atoms according to residue templates: 1 Heavy 1191 H / lone pairs
> source leaprc.water.tip3p
----- Source: /home/wangq/Programs/amber16/dat/leap/cmd/leaprc.water.tip3p ----- Source of /home/wangq/Programs/amber16/dat/leap/cmd/leaprc.water.tip3p done Loading library: /home/wangq/Programs/amber16/dat/leap/lib/atomic_ions.lib Loading library: /home/wangq/Programs/amber16/dat/leap/lib/solvents.lib Loading parameters: /home/wangq/Programs/amber16/dat/leap/parm/frcmod.ionsjc_tip3p Reading force field modification type file (frcmod) Reading title: Monovalent ion parameters for Ewald and TIP3P water from Joung & Cheatham JPCB (2008) Loading parameters: /home/wangq/Programs/amber16/dat/leap/parm/frcmod.ions234lm_126_tip3p Reading force field modification type file (frcmod) Reading title: Li/Merz ion parameters of divalent to tetravalent ions for TIP3P water model (12-6 normal usage set)
> solvateBox mol TIP3P 10
Solute vdw bounding box: 49.995 53.684 37.063 Total bounding box for atom centers: 69.995 73.684 57.063 Solvent unit box: 18.774 18.774 18.774 Total vdw box size: 73.439 76.851 60.149 angstroms. Volume: 339472.593 A^3 Total mass 162173.797 amu, Density 0.793 g/cc Added 8064 residues.
> charge mol
Total unperturbed charge: 2.000000 Total perturbed charge: 2.000000
> addIons2 mol Cl- 0
2 Cl- ions required to neutralize. Adding 2 counter ions to "mol" using 1A grid Grid extends from solute vdw + 2.51 to 8.51 Resolution: 1.00 Angstrom. grid build: 1 sec Calculating grid charges charges: 149 sec Placed Cl- in mol at (0.54, -17.17, 29.94). Placed Cl- in mol at (34.54, -37.17, -19.06). Done adding ions.
addions: simply draw a grid around the solute and places ions at grid points where the energies are lowest. This approach will ignore water molecules in locating where to place the ion and if the chosen location overlaps a water molecule, the water is deleted and replaced with the ion.
addions2: it treats solvent molecules the same as solute.
扫描二维码关注公众号,回复:
1178162 查看本文章
> saveAmberParm mol wt1mg.parm7 wt1mg.crd
Checking Unit. Building topology. Building atom parameters. Building bond parameters. Building angle parameters. Building proper torsion parameters. Building improper torsion parameters. total 468 improper torsions applied Building H-Bond parameters. Incorporating Non-Bonded adjustments. Not Marking per-residue atom chain types. Marking per-residue atom chain types. (Residues lacking connect0/connect1 - these don't have chain types marked: res total affected CGLN 1 NCYS 1 WAT 8064 ) (no restraints)