On-the-Fly Specifications of Reaction Coordinates in Parallel Cascade Selection Molecular Dynamics Accelerate Conformational Transitions of Proteins.

Abstract

Parallel Cascade Selection Molecular Dynamics (PaCS-MD) is an efficient conformational sampling method for generating a set of reactive trajectories that connect a given reactant and a product. In PaCS-MD, initial structures relevant to conformational transitions are reasonably selected by referring to a set of reaction coordinates (RCs), and short-time molecular dynamics (MD) simulations are independently launched from them. To efficiently perform PaCS-MD, specifications of RCs are essential, but specifying reasonable RCs is generally nontrivial. In the present study, we propose on-the-fly specifications of RCs as an extended PaCS-MD. In the present method, n types of RCs are provided as candidates a priori as follows: RC = (X1, X2, ..., X n), and one of the RCs is specified in a cycle-dependent manner, i.e. the reasonable RC is searched at every cycle by evaluating gradients of the RCs, i.e. RC with the steepest gradient for cycle is regarded as the reasonable RC, and conformational resampling proceeds along it, promoting conformational transition of a given protein. For a demonstration, the extended PaCS-MD was applied to reproduce the open-closed conformational transition of T4 lysozyme (T4L). As candidates of possible RCs, (1) root-mean square distance, (2) principal coordinates, (3) accessible surface area, (4) radius of gyration, and (5) end-to-end distance were adopted in the cycle-dependent specifications of RCs. Through the demonstration, the extended PaCS-MD successfully reproduced the conformational transition from the open to closed states of T4L. As a more complicated practice, a dimerization process of diubiquitin was efficiently reproduced with the extended PaCS-MD, showing the high conformational sampling efficiency of the present algorithm. In contrast, the conventional PaCS-MD with a fixed RC sometimes failed to generate a set of reactive trajectories when an unreasonable RC was specified, i.e. the conformational sampling efficiency of PaCS-MD might more or less depend on the specified RCs. Judging from the present demonstrations, on-the-fly specifications of RCs might be effective in reproducing/predicting essential transitions of a given protein.