File: generate_amber.py

package info (click to toggle)
pymol 1.2r2-1.1
  • links: PTS, VCS
  • area: main
  • in suites: squeeze
  • size: 38,716 kB
  • ctags: 23,006
  • sloc: ansic: 480,417; python: 70,953; cpp: 12,928; sh: 10,266; makefile: 530; csh: 21
file content (19 lines) | stat: -rw-r--r-- 409 bytes parent folder | download | duplicates (8) 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
 
# python

from chempy import io
from chempy import protein
from chempy import protein_amber99

model= io.pdb.fromFile("../../test/dat/pept.pdb")

model= protein.generate(model,forcefield=protein_amber99)

sm = 0
for a in model.atom:
   sm = sm + a.partial_charge

print " prot: net partial charge on protein is %8.4f" % sm
print " prot: (this should be integral)!"

io.pkl.toFile(model,"generate_amber.pkl")