The growing adoption of generalized-ensemble algorithms for biomolecular simulation has resulted in a resurgence in the use of the weighted histogram analysis method (WHAM) to make use of all data generated by these simulations. Unfortunately, the original presentation of WHAM by Kumar et al. is not directly applicable to data generated by these methods. WHAM was originally formulated to combine data from independent samplings of the canonical ensemble, whereas many generalized-ensemble algorithms sample from mixtures of canonical ensembles at different temperatures. Sorting configurations generated from a parallel tempering simulation by temperature obscures the temporal correlation in the data and results in an improper treatment of the statistical uncertainties used in constructing the estimate of the density of states. Here we present variants of WHAM, STWHAM and PTWHAM, derived with the same set of assumptions, that can be directly applied to several generalized ensemble algorithms, including simulated tempering, parallel tempering (better known as replica-exchange among temperatures), and replica-exchange simulated tempering. We present methods that explicitly capture the considerable temporal correlation in sequentially generated configurations using autocorrelation analysis. This allows estimation of the statistical uncertainty in WHAM estimates of expectations for the canonical ensemble. We test the method with a one-dimensional model system and then apply it to the estimation of potentials of mean force from parallel tempering simulations of the alanine dipeptide in both implicit and explicit solvent.
Download full-text PDF |
Source |
|---|---|
| http://dx.doi.org/10.1021/ct0502864 | DOI Listing |
Publication Analysis
Top Keywords
Similar Publications
Strong electrostatic attraction drives milk heteroprotein complex coacervation.
Int J Biol Macromol
November 2024
Division of Computational Chemistry, Lund University, Naturvetarvägen 24, SE-223 62 Lund, Sweden; LINXS - Institute of advanced Neutron and X-ray Science, Lund University, Scheelevägen 19, 223 70 SE-Lund, Sweden. Electronic address:
Coacervates of oppositely charged milk proteins are used in functional food development, mainly to encapsulate bioactives. To uncover the driving forces behind coacervates formation, we study the association of lactoferrin and β-lactoglobulin at amino-acid level detail, using molecular simulations. Our findings show that inter-protein electrostatic interactions dominate and are, surprisingly, equally divided between an isotropic part, due to monopole-monopole attraction of the oppositely charged proteins, and an anisotropic part due to uneven surface charge distributions.
View Article and Find Full Text PDFAll-to-all reconfigurability with sparse and higher-order Ising machines.
Nat Commun
October 2024
Department of Electrical and Computer Engineering, University of California, Santa Barbara, Santa Barbara, CA, 93106, USA.
Domain-specific hardware to solve computationally hard optimization problems has generated tremendous excitement. Here, we evaluate probabilistic bit (p-bit) based Ising Machines (IM) on the 3-Regular 3-Exclusive OR Satisfiability (3R3X), as a representative hard optimization problem. We first introduce a multiplexed architecture that emulates all-to-all network functionality while maintaining highly parallelized chromatic Gibbs sampling.
View Article and Find Full Text PDFAdsorption of cytochrome c on different self-assembled monolayers: The role of surface chemistry and charge density.
Biointerphases
September 2024
School of Chemistry and Chemical Engineering, Guangdong Provincial Key Lab for Green Chemical Product Technology, South China University of Technology, Guangzhou 510640, People's Republic of China.
Article Synopsis
- The study investigates how cytochrome c (Cyt-c) adsorbs onto five types of self-assembled monolayers (SAMs) through simulations, focusing on factors like surface chemistry and charge density.
- Cyt-c displays strong binding to CH3-SAM due to hydrophobic interactions while encountering weak adsorption on OH-SAM; lower charge density NH2-SAM does not support stable binding under high ionic conditions.
- The structural integrity of Cyt-c is generally maintained on most SAMs, but its redox properties are compromised on highly negative OSO3--SAM surfaces due to changes in orientation and conformation.
Structural classification of Ag and Cu nanocrystals with machine learning.
Nanoscale
September 2024
Department of Chemical Engineering and Department of Physics, The Pennsylvania State University, University Park, PA 16802, USA.
We use machine learning (ML) to classify the structures of mono-metallic Cu and Ag nanoparticles. Our datasets comprise a broad range of structures - both crystalline and amorphous - derived from parallel-tempering molecular dynamics simulations of nanoparticles in the 100-200 atom size range. We construct nanoparticle features using common neighbor analysis (CNA) signatures, and we utilize principal component analysis to reduce the dimensionality of the CNA feature set.
View Article and Find Full Text PDFEnergetics and spectroscopic studies of clusters and the temperature dependencies of the isomers: An approach based on a combined recipe of parallel tempering and quantum chemical methods.
J Comput Chem
December 2024
Department of Chemistry, Adamas University, Kolkata, India.
A system associated with several number of weak interactions supports numerous number of stable structures within a narrow range of energy. Often, a deterministic search method fails to locate the global minimum geometry as well as important local minimum isomers for such systems. Therefore, in this work, the stochastic search technique, namely parallel tempering, has been executed on the quantum chemical surface of the system for -8 to generate global minimum as well as several number of local minimum isomers.
View Article and Find Full Text PDFWant AI Summaries of new PubMed Abstracts delivered to your In-box?
Enter search terms and have AI summaries delivered each week - change queries or unsubscribe any time!