Premium
Protein conformational landscapes: Energy minimization and clustering of a long molecular dynamics trajectory
Author(s) -
Troyer John M.,
Cohen Fred E.
Publication year - 1995
Publication title -
proteins: structure, function, and bioinformatics
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.699
H-Index - 191
eISSN - 1097-0134
pISSN - 0887-3585
DOI - 10.1002/prot.340230111
Subject(s) - dihedral angle , maxima and minima , energy minimization , molecular dynamics , trajectory , chemistry , root mean square , energy landscape , cluster (spacecraft) , nanosecond , kinetic energy , minification , cluster analysis , crystallography , biological system , physics , computational chemistry , mathematics , molecule , computer science , mathematical analysis , statistics , classical mechanics , mathematical optimization , hydrogen bond , laser , optics , biology , biochemistry , quantum mechanics , programming language , astronomy , organic chemistry
Abstract Using energy minimization and cluster analysis, we have analyzed a 1020 ps molecular dynamics trajectory of solvated bovine pancreatic trypsin inhibitor. Elucidation of conformational sub states in this way both illustrates the degree of conformational convergence in the simulation and reduces the structural data to a tractable subset. The relative movement of structures upon energy minimization was used to estimate the sizes of features on the protein potential energy surface. The structures were analyzed using their pairwise root‐mean‐square C α deviations, which gave a global measure of conformational changes that would not be apparent by monitoring single degrees of freedom. At time scales of 0.1 ps, energy minimization detected sharp transitions between energy minima separated by 0.1 Å rms deviation. Larger conformational clusters containing these smaller minima and separated by 0.25 Å were seen at 1 ps time scales. Both of these small features of the conformational landscape were characterized by movements in loop regions associated with small, correlated backbone dihedral angle shifts. On a nanosecond time scale, the main features of the protein energy landscape were clusters separated by over 0.7 Å rms deviation, with only seven of these sub states visited over the 1 ns trajectory. These substates, discernible both before and after energy minimization, differ mainly in a monotonic pivot of the loop residues 11–18 over the course of the simulation. This loop contains lysine 17, which specifically binds to trypsin in the active site. The trajectory did not return to previously visited clusters, indicating that this trajectory has not been shown to have completely sampled the conformational substates available to it. Because the apparent convergence to a single region of conformation space depends on both the time scale of observation and the size of the conformational features examined, convergence must be operationally defined within the context of the simulation. © 1995 Wiley‐Liss, Inc.