Skip directly to search Skip directly to A to Z list Skip directly to page options Skip directly to site content

NIOSHTIC-2 Publications Search

Search Results

Efficient calculation of configurational entropy from molecular simulations by combining the mutual-information expansion and nearest-neighbor methods.

Authors
Hnizdo-V; Tan-J; Killian-BJ; Gilson-MK
Source
J Comput Chem 2007 Dec; 29(10):1605-1614
NIOSHTIC No.
20034096
Abstract
Changes in the configurational entropies of molecules make important contributions to the free energies of reaction for processes such as protein-folding, noncovalent association, and conformational change. However, obtaining entropy from molecular simulations represents a long-standing computational challenge. Here, two recently introduced approaches, the nearest-neighbor (NN) method and the mutual-information expansion (MIE), are combined to furnish an efficient and accurate method of extracting the configurational entropy from a molecular simulation to a given order of correlations among the internal degrees of freedom. The resulting method takes advantage of the strengths of each approach. The NN method is entirely nonparametric (i.e., it makes no assumptions about the underlying probability distribution), its estimates are asymptotically unbiased and consistent, and it makes optimum use of a limited number of available data samples. The MIE, a systematic expansion of entropy in mutual information terms of increasing order, provides a well-characterized approximation for lowering the dimensionality of the numerical problem of calculating the entropy of a high-dimensional system. The combination of these two methods enables obtaining well-converged estimations of the configurational entropy that capture many-body correlations of higher order than is possible with the simple histogramming that was used in the MIE method originally. The combined method is tested here on two simple systems: an idealized system represented by an analytical distribution of six circular variables, where the full joint entropy and all the MIE terms are exactly known, and the R,S stereoisomer of tartaric acid, a molecule with seven internal-rotation degrees of freedom for which the full entropy of internal rotation has been already estimated by the NN method. For these two systems, all the expansion terms of the full MIE of the entropy are estimated by the NN method and, for comparison, the MIE approximations up to third order are also estimated by simple histogramming. The results indicate that the truncation of the MIE at the two-body level can be an accurate, computationally nondemanding approximation to the configurational entropy of anharmonic internal degrees of freedom. If needed, higher-order correlations can be estimated reliably by the NN method without excessive demands on the molecular-simulation sample size and computing time.
Keywords
Risk-analysis; Risk-factors; Toxic-effects; Computer-models; Statistical-analysis; Mathematical-models
Contact
Vladimir Hnizdo, Health Effects Laboratory Division, National Institute for Occupational Safety and Health, Morgantown, West Virginia 26505
CODEN
JCCHDD
Publication Date
20071201
Document Type
Journal Article
Email Address
vhnizdo@cdc.gov
Funding Type
Contract
Fiscal Year
2008
NTIS Accession No.
NTIS Price
Identifying No.
Contract-212-2006-M-16936
Issue of Publication
10
ISSN
0192-8651
NIOSH Division
HELD
Source Name
Journal of Computational Chemistry
State
WV; MD
TOP