First SimulationΒΆ
This is a good starting point for learning how to use OpenMM.
It loads a PDB file villin.pdb, which is the villin headpiece protein in a box of water. In then parameterizes it using the Amber14 forcefield and TIP3P-FB water model, energy minimizes it, and simulates it for 1000 steps with a langevin intergator.
[1]:
from openmm.app import *
from openmm import *
from openmm.unit import *
from sys import stdout
# Load in the PDB strucure
pdb = PDBFile('villin.pdb')
# Specifiy the forcefield
forcefield = ForceField('amber14-all.xml', 'amber14/tip3pfb.xml')
# Combine the molecular topology and the forcefield
system = forcefield.createSystem(pdb.topology, nonbondedMethod=PME,
nonbondedCutoff=1*nanometer, constraints=HBonds)
# Create the integrator to use for advacing the equations of motion.
# It specifies a Langevin integrator.
# The paramters set are temperature, friction coefficient, and timestep.
integrator = LangevinMiddleIntegrator(300*kelvin, 1/picosecond, 0.004*picoseconds)
# Combines the molecular topology, system, and integrator
# to begin a new simulation.
simulation = Simulation(pdb.topology, system, integrator)
simulation.context.setPositions(pdb.positions)
# Perform local energy minimization
print("Minimizing energy...")
simulation.minimizeEnergy(maxIterations=100)
# Write the trajectory to a file called "output.pdb"
simulation.reporters.append(PDBReporter('output.pdb', 1000))
# Report infomation to the screen as the simulation runs
simulation.reporters.append(StateDataReporter(stdout, 100, step=True,
potentialEnergy=True, temperature=True))
#Run the simulation for 1000 timsteps
print("Running simulation...")
simulation.step(1000)
Minimizing energy...
Running simulation...
#"Step","Potential Energy (kJ/mole)","Temperature (K)"
100,-156931.0346892771,120.06221506140332
200,-152813.02042911295,180.5264760128298
300,-149855.82419335766,219.61778066500156
400,-147462.5606267873,243.8422724307354
500,-145754.30790993015,260.08388038934595
600,-144368.09850578316,272.1621746708566
700,-143762.1432821709,285.84618341733983
800,-143233.06469224993,285.857237324155
900,-142705.47430668795,287.63641339336044
1000,-141817.4966432689,289.89614604970603
