Ctrl+K 
from htmd.ui import *
config(viewer='ngl')
os.chdir('/webdata/nc983hu3brda/') # Don't use this command

Using docking to generate starting poses for simulations#

Download the files for this tutorial from this link

Dock the protein with the ligand#

prot = Molecule('bentryp/trypsin.pdb')
prot.center()
lig = Molecule('bentryp/benzamidine.pdb')
poses, scores = dock(prot, lig)

Visualize the docked poses#

mol = Molecule()
mol.append(prot)
for i, p in enumerate(poses):
    mol.append(p)
mol.view(sel='protein', style='NewCartoon', hold=True)
mol.view(sel='resname MOL', style='Licorice', color=1)

Build systems from docked poses#

molbuilt = []
for i, p in enumerate(poses):
    prot = Molecule('bentryp/trypsin.pdb')
    prot.filter('chain A and (protein or water or resname CA)')
    prot.set('segid', 'P', sel='protein and noh')
    prot.set('segid', 'W', sel='water')
    prot.set('segid', 'CA', sel='resname CA')
    prot.center()
    from moleculekit.util import maxDistance
    D = maxDistance(prot, 'all')

    ligand = p
    ligand.set('segid','L')
    ligand.set('resname','MOL')

    mol = Molecule(name='combo')
    mol.append(prot)
    mol.append(ligand)

    D = D + 15
    smol = solvate(mol, minmax=[[-D, -D, -D], [D, D, D]])
    topos  = ['top/top_all22star_prot.rtf', 'top/top_water_ions.rtf', './bentryp/benzamidine.rtf']
    params = ['par/par_all22star_prot.prm', 'par/par_water_ions.prm', './bentryp/benzamidine.prm']

    molbuilt.append(charmm.build(smol, topo=topos, param=params, outdir='./docked/build/{}/'.format(i+1), saltconc=0.15))
    if i==1: # For time purposes lets only build the two first
        break

Equilibrate the build systems#

from htmd.protocols.equilibration_v3 import Equilibration
md = Equilibration()
md.numsteps = 1000
md.temperature = 298

builds = natsort(glob('docked/build/*/'))
for i, b in enumerate(builds):
    md.write(b, 'docked/equil/{}/'.format(i+1))

mdx = LocalGPUQueue()
mdx.submit(glob('./docked/equil/*/'))

mdx.wait()

Create the production folder#

from htmd.protocols.production_v6 import Production
md = Production()
md.runtime = 50
md.timeunits = 'ns'
md.temperature = 298

equils = natsort(glob('docked/equil/*/'))
for i, b in enumerate(equils):
    md.write(b, 'docked/generators/{}/'.format(i+1))

mdx = LocalGPUQueue()
mdx.submit(glob('./docked/generators/*/'))

mdx.wait()
On this page