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,-156808.57561135953,117.78484508289314
200,-152600.55253085264,181.60134556768423
300,-149526.1089161045,217.24953316215505
400,-147500.05970184834,243.66002130804375
500,-146258.66347540068,262.71544094964736
600,-144762.89568025118,270.98800722788894
700,-144050.29091723793,279.103269372535
800,-143119.49190526138,282.82598807907766
900,-142676.27153689557,285.6939732844154
1000,-141984.03938527062,290.3643719155838
