Sl. Kazmirski et al., Analysis methods for comparison of multiple molecular dynamics trajectories: Applications to protein unfolding pathways and denatured ensembles, J MOL BIOL, 290(1), 1999, pp. 283-304
In molecular dynamics simulations of protein unfolding, the pathway of one
protein molecule is studied at a time. In contrast, experimental denaturati
on studies sample from large ensembles of molecules passing from the native
to unfolded state. If reasonable comparisons with experiment are to be mad
e, then the generality of the simulations needs to be con firmed by perform
ing multiple unfolding simulations. Given that protein unfolding trajectori
es are very complicated functions of the proteins and the environment, comp
aring different trajectories, even under the same conditions, is not straig
htforward. Several methods are presented here that attempt to accomplish th
is task at different levels of complexity. The simpler methods are geometry
based and make use: of the root-mean-squared deviations between structures
, while the more complicated methods are based on the time variation of the
various properties of the system during the unfolding process. These metho
ds are applied to multiple simulations of three different proteins, bovine
pancreatic trypsin inhibitor, chymotrypsin inhibitor 2, and barnase. In gen
eral, for these three proteins protein unfolding proceeded via expansion of
the core and fraying of secondary structure to yield the major transition
state. Once past the transition state, the trajectories for a given protein
diverged as the protein lost further secondary and tertiary structure by a
variety of mechanisms. Although the unfolding pathways diverged, similar c
onformations were populated in the denatured state even when the unfolding
occurred via different pathways. The multitude of different pathways leadin
g to the denatured state agrees with the funnel description of protein fold
ing. Although the pathways differed in conformational space, the physical p
roperties of the conformations were often similar, highlighting the danger
of assuming that similar observed properties imply similar conformations. I
n fact, there may be many different "conformational pathways" of unfolding
that fit within a preferred "property space pathway". (C) 1999 Academic Pre
ss.