For internal use
For internal use
reduce sustiva.pdb > sustiva_h.pdb
python convert.py
convert.py:
#!/usr/bin/python
from pdblib.base import *
prot=Mol('./sustiva_h.pdb')
for res in getreses(prot):
res.name='SUS'
prot.renumber(1,1)
prot.write('sustiva_new.pdb')
antechamber -i sustiva_new.pdb -fi pdb -o sustiva.mol2 -fo mol2 -c bcc -s 2
parmchk -i sustiva.mol2 -f mol2 -o sustiva.frcmod
tleap -f leap.in
leap.in: source leaprc.ff99SB source leaprc.gaff SUS = loadmol2 sustiva.mol2 check SUS loadamberparams sustiva.frcmod saveoff SUS sus.lib saveamberparm SUS sustiva.prmtop sustiva.inpcrd quit
tleap -f leap2.in
leap2.in: source leaprc.ff99SB source leaprc.protein.ff14SB source leaprc.gaff loadamberparams /home/panying/work/MDdock_jiangp/analysis/1.2_IP4par_incr/test/sustiva.frcmod loadoff /home/panying/work/MDdock_jiangp/analysis/1.2_IP4par_incr/test/sus.lib complex = loadpdb 1FKO_trunc_sus.pdb saveamberparm complex 1FKO_sus.prmtop 1FKO_sus.inpcrd savepdb complex 1FKO_sus.pdb quit