Tutorial 7.3: Simulated Annealing with Monte Carlo
This tutorial is nearly identical to Tutorial 7.2 (Simulated Annealing with Molecular Dynamics), except that here we use the Monte Carlo algorithm (MC), rather than MD. Both MC and MD generate a trajectory, but the two algorithms use different rules for generating successive conformations. MD takes small steps in time, numerically integrating Newton's equations of motion. It generates successive conformations that follow the evolution of the system in time. In contrast, MC takes small random jumps in conformational space, with the probability of selecting new conformations dictated by the Boltzmann distribution.
SET UP UNIVERSE:
Create a new folder for your work on this tutorial.
Open Chimera, Launch Sophia, and initiate a new simulation with:
PDB File: LJ10-2.pdb
Force Field: Lennard-Jones
FF Mod File: LJ_E0_0.1.frcmod
Box length: 2 nm
Non-bond cutoff: 1 nm
(Click on "Load Universe")
In the UCSF Chimera window, click Actions / Atoms/Bonds / sphere to generate a space-filling model.
FIRST ANNEALING RUN:
In the Sophia Recipe window:
Click on the "MC" tab (1) to be sure the active recipe is MD.
Click "New" (2). The default MD recipe will be loaded.
Click on the recipe (3). The blue band means that the recipe is active and can be edited.
NOTE THAT WE DON'T HAVE A 2D IMPLEMENTATION OF MC, AND I'VE NOT EDITED THIS PAGE BELOW THIS POINT TO MAKE IT CORRECT FOR MC.... I COPIED TUTORIAL 7.2 TO USE AS A TEMPLATE HERE; THE FIGURES ARE FROM THAT PAGE AND NEED REPLACING....
Edit the ingredient parameters (4). If necessary, drag the slider down (5) to expose the lower part of the window so you can edit the values there (6). Set the parameters as follows:
• Number of Steps: 100,000
• Output Interval: 100 steps
• Time Step: 10 fs
• Initialize Velocities: xy (keeps atoms in xy plane)
• Thermostat: Rescale
• Initial Temperature: 20K
• Final Temperature: 1K
• Rescale Interval: 50 steps
(Ignore the next four parameters)
• Remove Translation Interval: 0
• Remove Rotation Interval: 0
Be sure the recipe is still highlighted in blue (3).
Click on "Run Next Ingredient" (7).
You should get results something like this:
(You can plot temperature by checking the appropriate box in the Viewer Controls window.)
ANALYZE FIRST ANNEALING RUN:
The trajectory you have just completed has the information you need to select appropriate values for T1 and T2, the upper and lower temperatures for simulated annealing on this system.
T1 needs to be high enough that the initial structure is completely broken up, so the system freely explores conformational space.
T2 needs to be low enough that the number and pattern of interatomic contacts does not change, so the system is at a minimum energy conformation. This will allow us to produce a collection of minimum energy conformations in a series of simulated annealing runs. If our annealing protocol is well designed, we should have a reasonable chance of finding the global minimum.
Scroll through the trajectory using the controls in the Trajectory Viewer window. Watch the atomic motions in the Chimera window. Note how the pattern of motions changes as the system is cooled, and how those changes are correlated with changes in the slope of the plot of the potential energy.
Write a short summary of the behavior of the system at different temperatures. You should include snapshots like the one above, to clarify the correlations between structure and temperature.
Use this description to select appropriate values of T1 and T2.
SET UP AND RUN ANNEALING PROTOCOL:
In the Sophia Recipe window, be sure that Recipe #1 is still highlighted in blue (3, in the figure above).
Click the "Clone" button to make a copy of the original recipe.
Click the box "Make all ingredients deletable".
Be sure that Recipe #1 is still highlighted in blue.
Click on the "Delete" button. Recipe #1 should disappear, and the new recipe will change from #2 to #1.
Click on the new recipe #1 so you can edit it.
Edit this recipe by changing the values of the initial and final temperatures to the values you selected for T1 and T2. Set the number of steps to 10,000 and the output interval to 10 steps, so you will again generate a trajectory with 1000 frames.
While recipe #1 is highlighted in blue, click on the "Clone" button nine times, creating recipes #2 – #10.
Click on "Run All Ingredients".
The energy plot should look something like the frame at the left. You can expand the view as shown at the right by using the "zoom rectangle" option.
This particular simulation reached the global minimum energy structure three times (red circles), so the protocol had a success rate of 30%.
What was the success rate for your protocol?
VARY THE LENGTH OF THE ANNEALING PROTOCOL:
In the Sophia Recipe window, delete recipes #2-10.
Clone recipe #1 (whose status is "Done") to create recipe #2 ("Initialized").
Delete recipe #1.
Following the procedure above, edit the recipe to change the number of steps from 10,000 to 100,000 and the output interval from 10 to 100, so you still get 1000 frames per trajectory.
Make nine more clones of this recipe and run all ten of them.
What is the success rate now?
Create ten recipes with 1,000,000 steps each and an output interval of 1000.
What is the success rate now?