Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Improve automatical setup for production protocols #898

Open
stefdoerr opened this issue Aug 12, 2019 · 0 comments
Open

Improve automatical setup for production protocols #898

stefdoerr opened this issue Aug 12, 2019 · 0 comments

Comments

@stefdoerr
Copy link
Contributor

stefdoerr commented Aug 12, 2019

Improve automatical setup for production (and maybe equilibration) protocols. Especially when starting production directly from i.e. pdb+prmtop you need to modify a large amount of parameters to get it to work.

We always assume that production will follow equlibration and will use the files produced by acemd but in many cases you might want to start from pre-equilibrated structures even made from other software.

Example:

from htmd.ui import *
from htmd.protocols.production_v6 import Production

pd = Production(_version=3)
pd.runtime = 100
pd.timeunits = 'ns'
pd.acemd.coordinates = 'structure.pdb'
pd.acemd.structure = 'structure.pdb'
pd.acemd.parameters = 'structure.prmtop'
pd.acemd.bincoordinates = None
pd.acemd.boxsize = '100.35 54.325 168.2823'
pd.acemd.extendedsystem = None
pd.acemd.parmfile = 'structure.prmtop'
pd.write('./equil_finished/', './prod/')

It should really be much simpler than that

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Projects
None yet
Development

No branches or pull requests

1 participant