Skip to content

Instantly share code, notes, and snippets.

@sqamara
Last active February 8, 2016 20:00
Show Gist options
  • Save sqamara/6ed6ddfc904a5224e558 to your computer and use it in GitHub Desktop.
Save sqamara/6ed6ddfc904a5224e558 to your computer and use it in GitHub Desktop.
HETATM 1 Br Br 1 0.000 0.000 5.000 1.00 0.00
from simtk.openmm import app
pdb_file="./br.pdb"
#HETATM 1 Br Br 1 0.000 0.000 5.000 1.00 0.00
pdb = app.PDBFile(pdb_file)
for chain in pdb.topology.chains():
for res in chain.residues():
for atom in res.atoms():
print(atom.element)
# prints <Element boron>
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment