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,-153994.86970545992,148.07974778375555
200,-150505.82946791913,199.45699371655266
300,-148033.13439363384,231.95264167248695
400,-146654.86137275695,256.7863474543448
500,-145128.7831932183,269.54475258276545
600,-143973.43240186953,281.78452974574395
700,-143177.10946037897,285.580674199192
800,-143181.7195059284,297.4403384309955
900,-142380.79630130428,293.58805205463676
1000,-141898.89963412852,294.12376694534345
