TY - JOUR
T1 - CoCo-MD: A Simple and Effective Method for the Enhanced Sampling of Conformational Space
AU - Shkurti, Ardita
AU - Styliari, Ioanna Danai
AU - Balasubramanian, Vivek
AU - Bethune, Iain
AU - Pedebos, Conrado
AU - Jha, Shantenu
AU - Laughton, Charles A.
PY - 2019/1/8
Y1 - 2019/1/8
N2 - CoCo ("complementary coordinates") is a method for ensemble enrichment based on principal component analysis (PCA) that was developed originally for the investigation of NMR data. Here we investigate the potential of the CoCo method, in combination with molecular dynamics simulations (CoCo-MD), to be used more generally for the enhanced sampling of conformational space. Using the alanine penta-peptide as a model system, we find that an iterative workflow, interleaving short multiple-walker MD simulations with long-range jumps through conformational space informed by CoCo analysis, can increase the rate of sampling of conformational space up to 10 times for the same computational effort (total number of MD timesteps). Combined with the reservoir-REMD method, free energies can be readily calculated. An alternative, approximate but fast and practically useful, alternative approach to unbiasing CoCo-MD generated data is also described. Applied to cyclosporine A, we can achieve far greater conformational sampling than has been reported previously, using a fraction of the computational resource. Simulations of the maltose binding protein, begun from the "open" state, effectively sample the "closed" conformation associated with ligand binding. The PCA-based approach means that optimal collective variables to enhance sampling need not be defined in advance by the user but are identified automatically and are adaptive, responding to the characteristics of the developing ensemble. In addition, the approach does not require any adaptations to the associated MD code and is compatible with any conventional MD package.
AB - CoCo ("complementary coordinates") is a method for ensemble enrichment based on principal component analysis (PCA) that was developed originally for the investigation of NMR data. Here we investigate the potential of the CoCo method, in combination with molecular dynamics simulations (CoCo-MD), to be used more generally for the enhanced sampling of conformational space. Using the alanine penta-peptide as a model system, we find that an iterative workflow, interleaving short multiple-walker MD simulations with long-range jumps through conformational space informed by CoCo analysis, can increase the rate of sampling of conformational space up to 10 times for the same computational effort (total number of MD timesteps). Combined with the reservoir-REMD method, free energies can be readily calculated. An alternative, approximate but fast and practically useful, alternative approach to unbiasing CoCo-MD generated data is also described. Applied to cyclosporine A, we can achieve far greater conformational sampling than has been reported previously, using a fraction of the computational resource. Simulations of the maltose binding protein, begun from the "open" state, effectively sample the "closed" conformation associated with ligand binding. The PCA-based approach means that optimal collective variables to enhance sampling need not be defined in advance by the user but are identified automatically and are adaptive, responding to the characteristics of the developing ensemble. In addition, the approach does not require any adaptations to the associated MD code and is compatible with any conventional MD package.
UR - http://www.scopus.com/inward/record.url?scp=85062821676&partnerID=8YFLogxK
U2 - 10.1021/acs.jctc.8b00657
DO - 10.1021/acs.jctc.8b00657
M3 - Article
AN - SCOPUS:85062821676
SN - 1549-9618
VL - 15
SP - 2587
EP - 2596
JO - Journal of Chemical Theory and Computation
JF - Journal of Chemical Theory and Computation
IS - 4
ER -