Entropy of Two-Molecule Correlated Translational-Rotational Motions Using the kth Nearest Neighbor Method

Data Structures and Algorithms for k-th Nearest Neighbours Conformational Entropy Estimation

Entropy of multivariate distributions may be estimated based on the distances of nearest neighbours from each sample from a statistical ensemble. This technique has been applied on biomolecular systems for estimating both conformational and translational/rotational entropy. The degrees of freedom which mostly define conformational entropy are torsion angles with their periodicity. In this work, tree structures and algorithms to quickly generate lists of nearest neighbours for periodic and non-periodic data are reviewed and applied to biomolecular conformations as described by torsion angles. The effect of dimensionality, number of samples, and number of neighbours on the computational time is assessed. The main conclusion is that using proper data structures and algorithms can greatly reduce the complexity of nearest neighbours lists generation, which is the bottleneck step in nearest neighbours entropy estimation.

Entropic Path Sampling: Computational Protocol to Evaluate Entropic Profile along a Reaction Path

Fleeting intermediates constitute dynamically stepwise mechanisms. They have been characterized in molecular dynamics trajectories, but whether these intermediates form a free energy minimum to become entropic intermediates remains elusively defined. We developed a computational protocol known as entropic path sampling to evaluate the entropic variation of reacting species along a reaction path based on an ensemble of trajectories. Using cyclopentadiene dimerization as a model reaction, we observed an entropy maximum along the reaction path which originates from an enhanced conformational flexibility as the reacting species enter into a flat energy region. As the reacting species further approach product formation, unfavorable entropic restriction fails to offset the potential energy drop, resulting in no free energy minimum along the post-TS pathway. Our results show that cyclopentadiene dimerization involves an entropy maximum that leads to dynamic intermediates with elongated lifetimes, but the reaction does not involve entropic intermediates.