Nov 28 2023 Feedback.
Contents
 
Introduction
Help Videos
Reference Guide
Getting Started
Protein Structure
Molecular Graphics
Slides & ActiveICM
Sequences & Alignments
Protein Modeling
Cheminformatics
Learn and Predict
Docking
Virtual Screening
Molecular Dynamics
MolScreen
3D Ligand Editor
Tables and Plots
Local Databases
ICM-Scarab
KNIME
Tutorials
FAQs
 
Index
PrevICM User's Guide
14 Molecular Dynamics
Next

Molecular Dynamics (MD) in ICM uses the functions in OpenMM - if you use this feature please cite OpenMM.

To run an MD simulations using Amber force field ff14SB:

  • Convert your protein to an ICM object.
  • MolMechanics/Run MD Simulation
  • Use the drop down box to select the protein.
  • Select whether or not you wish to solvate the protein.
  • Select the length of time for the simulation in ns.
  • Click OK.

To add restraints:

  • Make a copy of the protein you wish to simulate - right click on the name of the protein in the ICM workspace (left hand side) and choose the option 'clone'.
  • Select the region of the cloned molecule you wish to restrain and convert that selection to an 'orange' selection as described here.
  • Go to MolMechaics/Run MD Simulation as described above and choose the selection in the Restrain dialog box.
  • Click OK to run the simulation.

Results:

The conformation snapshots are stored at regular time intervals during the MD simulation and are embedded in the ICM object as a stack.

Results Analysis:

Right click on the stack of MD snapshots in the ICM workspace and choose Stack Calculations. In this dialog box you can analyze

  • Distance changes between snapshots.
  • RMSD to Reference - Read in a separate object as a reference and make the selection you want to analyze. You can choose to analyze the structures 'static' or superimpose. A table and plot of the RMSDs are displayed. This is a useful tool for analysis of a Molecular Dynamics run. Your x-axis will be the time in ns / number of conformations so if you chose 1ns and 50 conformations (snapshots) then each point is 1ns/50 = 0.02 ns
  • Torsion Changes - Select a pair of atoms on a rotatable bond to calculate the torsion angle.
  • Area Changes - use green and orange selections to calculate differences in area
  • RTCNN - use green and orange selections for receptor/ligand respectively to calculate RTCNN docking score.
  • ClusterTree - make a selection of atoms to compare and then a cluster tree will be displayed. You can then 'select-tree{select} a diverse set of conformations in the tree.


Prev
RIDE GPU Benchmark
Home
Up
Next
MolScreen

Copyright© 1989-2020, Molsoft,LLC - All Rights Reserved.
This document contains proprietary and confidential information of Molsoft, LLC.
The content of this document may not be disclosed to third parties, copied or duplicated in any form,
in whole or in part, without the prior written permission from Molsoft, LLC.