  A goal in the kinetic characterization of a macromolecular system is the
description of its slow relaxation processes, involving (i) identification of
the structural changes involved in these processes, and (ii) estimation of the
rates or timescales at which these slow processes occur. Most of the approaches
to this task, including Markov models, Master-equation models, and kinetic
network models, start by discretizing the high-dimensional state space and then
characterize relaxation processes in terms of the eigenvectors and eigenvalues
of a discrete transition matrix. The practical success of such an approach
depends very much on the ability to finely discretize the slow order
parameters. How can this task be achieved in a high-dimensional configuration
space without relying on subjective guesses of the slow order parameters? In
this paper, we use the variational principle of conformation dynamics to derive
an optimal way of identifying the "slow subspace" of a large set of prior order
parameters - either generic internal coordinates (distances and dihedral
angles), or a user-defined set of parameters. It is shown that a method to
identify this slow subspace exists in statistics: the time-lagged independent
component analysis (TICA). Furthermore, optimal indicators-order parameters
indicating the progress of the slow transitions and thus may serve as reaction
coordinates-are readily identified. We demonstrate that the slow subspace is
well suited to construct accurate kinetic models of two sets of molecular
dynamics simulations, the 6-residue fluorescent peptide MR121-GSGSW and the
30-residue natively disordered peptide KID. The identified optimal indicators
reveal the structural changes associated with the slow processes of the
molecular system under analysis.
