Evaluation of Dimensionality-Reduction Methods from Peptide Folding-Unfolding Simulations

摘要

Dimensionality-reduction methods have been widely used to study the free energy landscapes and low-free-energy pathways of molecular systems. It was shown that the nonlinear dimensionality-reduction methods gave better embedding results than the linear methods, such as principal component analysis, in some simple systems. In this study, we have evaluated several nonlinear methods, locally linear embedding, Isomap, and diffusion maps, as well as principal component analysis from the equilibrium folding/unfolding trajectory of the second β-hairpin of the Bl domain of streptococcal protein G. The CHARMM parm19 polar hydrogen potential function was used. A series of criteria which reflects different aspects of the embedding qualities was employed in the evaluation. Our results show that principal component analysis is not worse than the nonlinear ones on this complex system. There is no clear winner in all aspects of the evaluation. Each dimensionality-reduction method has its limitations in a certain aspect. We emphasize that a fair, informative assessment of an embedding result requires a combination of multiple evaluation criteria rather than any single one. Caution should be used when dimensionality-reduction methods are employed, especially when only a few of the top embedding dimensions are used to describe the free energy landscape.