Skip to content

Instantly share code, notes, and snippets.

@stefdoerr
Last active October 4, 2016 14:24
Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save stefdoerr/28b7a1cced1ad2dade640f6101eef258 to your computer and use it in GitHub Desktop.
Save stefdoerr/28b7a1cced1ad2dade640f6101eef258 to your computer and use it in GitHub Desktop.
building a system
In [40]: from htmd import *
In [41]: mol = Molecule('3PTB')
In [42]: mol.filter('protein')
In [43]: smol = solvate(mol)
In [44]: smol.chain
Out[44]: array(['A', 'A', 'A', ..., 'W', 'W', 'W'], dtype=object)
In [45]: bmol = charmm.build(smol, saltconc=0.15, outdir='/tmp/test/')
In [46]: bmol.chain
Out[46]: array(['A', 'A', 'A', ..., 'I', 'I', 'I'], dtype=object)
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment