Skip to main content
  • Oral presentation
  • Open access
  • Published:

Analysis and visual summarization of molecular dynamics simulation

Molecular dynamics (MD) simulation, a standard technique used to study the dynamical properties of biomolecules, is very useful in collecting the trajectories, a series of snapshots – the coordinates of the system - of larger systems for longer simulation times. These MD generated trajectories are huge in size (many gigabytes) and the data analysis may take much longer time than the data generation. Managing the large amount of data and presenting them in a flexible and comprehensible manner are the major challenges. Analyzing these trajectories with standard parameter like root-mean square deviation (RMSD) may not reveal the most interesting properties of the dynamics.

To overcome these challenges, Cα torsion angles [1] - torsion angles build by four consecutive Cα atoms - are highly valuable as similarity measure on a substructure scale and to find major events - the information on the time at which a transition occurs (temporal domain) and the local structural changes (spatial domain) of it is combinedly called as “event”- occurring in the course of the MD simulation. By calculating the time series of the Cα torsion angles and their clustering it is possible to determine the mechanistic details on a residual length scale and find major events occurring in the simulation of large proteins or protein complexes. The main advantage of the Cα torsion angle criterion is that it does not depend on a previous alignment of the structures, and that the direction of the change is also defined. Heat maps of Cα torsion angle give nice graphical representations of processes described by the MD simulations. Clustering of snapshots according to the specific Cα torsion angles is used to automatically find the spatial domains of the structural changes. If all the snapshots are assigned to a single cluster, then those residues are considered as rigid core and the remaining residues are considered as flexible parts. The temporal domain can be characterized in more detail by finding continuous time intervals assigned to a single cluster as (meta) stable structures and time intervals where the assignment jumps between two clusters as transitional periods. Since the outliers can be removed from the fuzzy clusters, starts and ends of time patches now qualify as important events for the underlying substructure and structural changes of larger regions are caused by an accumulation of such substructure events.

DNA polymerase I – the open ternary complex of the large fragment of Thermus aquaticus DNA polymerase I (Klentaq1), which is used here as a practical application for Cα torsion angle based analysis, shows a hand-like arrangement, including a thumb, a palm and a finger domain [2]. The catalytic cycle leading to nucleotide insertion comprises several steps including a large structural rearrangement in the form of a movement of the finger domain towards the thumb domain, i.e. the transition from the open to the closed form. Molecular dynamics simulation were performed using the AMBER 10 suite of programs [3]. To get the visual picture of the ongoing processes, the Cα torsion angles with the differences to the crystal structure of the open form were plotted as heat map. The rigid and the flexible parts were clearly seen with no or a large number of significant changes, respectively, from the heat map. Once the Cα torsion angles corresponding to the rigid parts are removed, the remaining regions change only in a specific time interval of the simulation. The spatial and temporal domains of the structural changes were identified automatically by clustering of snapshots (using KNIME [4]) and finding the continuous time intervals, respectively.

References

  1. Flocco MM, Mowbray SL: C alpha-based torsion angles: a simple tool to analyze protein conformational changes. Protein Sci. 1995, 4 (10): 2118-2122. 10.1002/pro.5560041017.

    Article  CAS  Google Scholar 

  2. Li Y, Korolev S, Waksman G: Crystal structures of open and closed forms of binary and ternary complexes of the large fragment of Thermus aquaticus DNA polymerase I: structural basis for nucleotide incorporation. EMBO J. 1998, 17: 7514-7525. 10.1093/emboj/17.24.7514.

    Article  CAS  Google Scholar 

  3. AMBER: Assisted Model Building with Energy Refinement. [http://ambermd.org]

  4. KNIME: The Konstanz Information Miner. [http://www.knime.org/]

Download references

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Fredrick Robin Devadoss.

Rights and permissions

Open Access  This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.

The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.

To view a copy of this licence, visit https://creativecommons.org/licenses/by/4.0/.

Reprints and permissions

About this article

Cite this article

Devadoss, F.R., Paul Raj, V. Analysis and visual summarization of molecular dynamics simulation. J Cheminform 6 (Suppl 1), O16 (2014). https://doi.org/10.1186/1758-2946-6-S1-O16

Download citation

  • Published:

  • DOI: https://doi.org/10.1186/1758-2946-6-S1-O16

Keywords