IMP  2.1.0
The Integrative Modeling Platform
serialPeptideDocker.py
1 import IMP
2 import subprocess
3 import random
4 import IMP.domino
5 import IMP.core
6 import IMP.rmf
7 import RMF
8 import time
9 import IMP.algebra
10 import types
11 import re
12 import sys
13 import operator
14 import os
15 import peptideDocker
16 import atomicDomino
17 
18 parameterFileName = sys.argv[1]
19 
20 startTime = time.time()
21 
22 p = peptideDocker.PeptideDocker(parameterFileName)
23 
24 p.createModel()
25 
26 p.loadHelpers()
27 
28 p.initDof()
29 
30 p.addForceFieldRestraints()
31 p.setInitialPositions()
32 p.addClosePairNonBondedRestraints()
33 
34 d = atomicDomino.AtomicDomino(p.getModel(), p.getProtein(), parameterFileName)
35 
36 d.loadDominoHelpers()
37 
38 d.createSubsets()
39 d.writeVisualization()
40 p.logTime("Created subsets")
41 p.addCompleteNonBondedRestraints()
42 
43 p.logTime("Added complete nonbonded restraints")
44 
45 
46 
47 #p.addExcludedVolume()
48 
49 p.logTime("Setup")
50 
51 p.runMolecularDynamics()
52 
53 p.logTime("Run MD")
54 
55 p.runAllCg()
56 
57 p.logTime("Run CG")
58 
59 d.createGrid()
60 
61 d.discretizeNativeProtein()
62 
63 particleNameList = d.getDominoParticleNames()
64 
65 flexibleAtoms = p.getFlexibleAtoms()
66 
67 d.readMdTrajectory(particleNameList, flexibleAtoms)
68 
69 d.readCgTrajectories(particleNameList, flexibleAtoms)
70 
71 p.logTime("Read Trajectory")
72 
73 d.createParticleStatesTable()
74 
75 p.logTime("create particle states table")
76 
77 d.createAllSubsetAssignments()
78 
79 p.logTime("Create Leaf assignments")
80 
81 d.createSampler()
82 
83 d.runDomino()
84 
85 p.logTime("Ran Domino")
86 
87 p.writeOutput()
88 
89 p.outputTimes()
90 
91 d.writeOutput(flexibleAtoms, startTime)
See IMP.core for more information.
See IMP.algebra for more information.
See IMP.rmf for more information.
Definition: associations.h:20
See IMP.domino for more information.
Definition: analysis.h:15