Principal component analysis of lipid molecule conformational changes in molecular dynamics simulations

Abstract : Molecular dynamics simulations of lipid bilayers are ubiquitous nowadays. Usually, either global properties of the bilayer or some particular characteristics of each lipid molecule are evaluated in such simulations, but the structural properties of the molecules as a whole are rarely studied. Here, we show how a comprehensive quantitative description of conformational space and dynamics of a single lipid molecule can be achieved via the principal component analysis (PCA). We illustrate the approach by analyzing and comparing simulations of DOPC bilayers obtained using eight different force fields: all-atom generalized AMBER, CHARMM27, CHARMM36, Lipid14 and Slipids, and united-atom Berger, GROMOS43A1-S3 and GROMOS54A7. Similarly to proteins, most of the structural variance of a lipid molecule can be described by only a few principal components. These major components are similar in different simulations, although there are notable distinctions between the older and newer force fields, and between the all-atom and united-atom force fields. The DOPC molecules in the simulations generally equilibrate on the timescales of tens to hundreds of nanoseconds. The equilibration is the slowest in the GAFF simulation, and the fastest in the Slipids simulation. Somewhat unexpectedly, the equilibration in the united-atom force fields is generally slower than in the all-atom force fields. Overall, there is a clear separation between the more variable previous generation force fields, and significantly more similar new generation force fields (CHARMM36, Lipid14, Slipids). We expect that the presented approaches will be useful for quantitative analysis of conformations and dynamics of individual lipid molecules in other simulations of lipid bilayers.
Type de document :
Article dans une revue
Journal of Chemical Theory and Computation, American Chemical Society, 2016, 12 (3), pp.1019-1028. 〈10.1021/acs.jctc.5b01106〉
Liste complète des métadonnées

https://hal.inria.fr/hal-01258167
Contributeur : Nano-D Equipe <>
Soumis le : lundi 18 janvier 2016 - 17:03:52
Dernière modification le : mercredi 11 avril 2018 - 01:58:59

Identifiants

Collections

Citation

Pavel Buslaev, Valentin I. Gordeliy, Sergei Grudinin, Ivan Yu. Gushchin. Principal component analysis of lipid molecule conformational changes in molecular dynamics simulations. Journal of Chemical Theory and Computation, American Chemical Society, 2016, 12 (3), pp.1019-1028. 〈10.1021/acs.jctc.5b01106〉. 〈hal-01258167〉

Partager

Métriques

Consultations de la notice

516