Reaction coordinates and transition pathways of rare events via forward flux sampling

A new approach is developed for identifying suitable reaction coordinates to describe the progression of rare events in complex systems. The method is based on the forward flux sampling (FFS) technique and standard least-square estimation (LSE) and it is denoted as FFS-LSE. The FFS algorithm generates trajectories for the transition between stable states as chains of partially connected paths, which can then be used to obtain "on-the-fly" estimates for the committor probability to the final region, p(B). These p(B) data are then used to screen a set of candidate collective properties for an optimal order parameter (i.e., reaction coordinate) that depends on a few relevant variables. LSE is used to find the coefficients of the proposed reaction coordinate model and an analysis of variance is used to determine the significant terms in the model. The method is demonstrated for several test systems, including the folding of a lattice protein. It is shown that a simple approximation to p(B) via a model linear on energy and number of native contacts is sufficient to describe the intrinsic dynamics of the protein system and to ensure an efficient sampling of pathways. In addition, since the p(B) surface found from the FFS-LSE approach leads to the identification of the transition state ensemble, mechanistic details of the dynamics of the system can be readily obtained during a single FFS-type simulation without the need to perform additional committor simulations.