Conformational and functional analysis of molecular dynamics trajectories by Self-Organising Maps

Abstract

BackgroundMolecular dynamics (MD) simulations are powerful tools to investigate the conformational dynamics of proteins that is often a critical element of their function. Identification of functionally relevant conformations is generally done clustering the large ensemble of structures that are generated. Recently, Self-Organising Maps (SOMs) were reported performing more accurately and providing more consistent results than traditional clustering algorithms in various data mining problems. We present a novel strategy to analyse and compare conformational ensembles of protein domains using a two-level approach that combines SOMs and hierarchical clustering.ResultsThe conformational dynamics of the α-spectrin SH3 protein domain and six single mutants were analysed by MD simulations. The Cα's Cartesian coordinates of conformations sampled in the essential space were used as input data vectors for SOM training, then complete linkage clustering was performed on the SOM prototype vectors. A specific protocol to optimize a SOM for structural ensembles was proposed: the optimal SOM was selected by means of a Taguchi experimental design plan applied to different data sets, and the optimal sampling rate of the MD trajectory was selected. The proposed two-level approach was applied to single trajectories of the SH3 domain independently as well as to groups of them at the same time. The results demonstrated the potential of this approach in the analysis of large ensembles of molecular structures: the possibility of producing a topological mapping of the conformational space in a simple 2D visualisation, as well as of effectively highlighting differences in the conformational dynamics directly related to biological functions.ConclusionsThe use of a two-level approach combining SOMs and hierarchical clustering for conformational analysis of structural ensembles of proteins was proposed. It can easily be extended to other study cases and to conformational ensembles from other sources.

Structural modeling of the AhR:ARNT complex in the bHLH–PASA–PASB region elucidates the key determinants of dimerization

“…1000 sets of loop models were generated; then the ensemble of models of each loop was clustered on the basis of the backbone structural similarity by using the Self Organizing Map (SOM) approach previously described. 40 The models representing the cluster medoids were selected as representative of the conformational variability of that loop.…”

Use of restrained molecular dynamics to predict the conformations of phosphorylated receiver domains in two‐component signaling systems

“…When performing geometrical clustering on MD simulations, we assume that conformations grouped into the same cluster are structurally similar, and therefore lie within the same basin on the freeenergy landscape of the protein. 116 Clusters are typically different sizes, due to conformational sampling progressing according to a Boltzmann distribution during simulation. Sparsely populated clusters are more likely to represent transient and/or higher energy states.…”

Analysis of Compounds Isolated from Gnetum gnemon L. Seeds as Potential ACE Inhibitors through Molecular Docking and Molecular Dynamics Simulations

“…Analysis of RMSF values aims to observe the flexibility of proteins occurring in each amino acid residue during molecular dynamics simulations. 20 We focused on the amino acid residues at the active site of ACE, including Gln281, His353, Ala354, Glu384, Glu411, Asp415, Asp453, Lys454, Lys511, His513, Tyr520, and Tyr523. Based on the graphs in Figure 4, it demonstrated that at temperatures of 300 K and 310 K all ligands that bonded with amino acid residues at the active sites of ACE did not exhibit a high RMSF value, so ligands were able to form stable conformation at the active site of ACE.…”