Temperature Weighted Histogram Analysis Method, Replica Exchange, and Transition Paths
Citations Over TimeTop 10% of 2005 papers
Abstract
We analyzed the data from a replica exchange molecular dynamics simulation using the weighted histogram analysis method to combine data from all of the temperature replicas (T-WHAM) to obtain the room-temperature potential of mean force of the G-peptide (the C-terminal beta-hairpin of the B1 domain of protein G) in regions of conformational space not sampled at room temperature. We were able to determine the potential of mean force in the transition region between a minor alpha-helical population and the major beta-hairpin population and identify a possible transition path between them along which the peptide retains a significant amount of secondary structure. This observation provides new insights into a possible mechanism of formation of beta-sheet secondary structures in proteins. We developed a novel Bayesian statistical uncertainty estimation method for any quantity derived from WHAM and used it to validate the calculated potential of mean force. The feasibility of estimating regions of the potential of mean force with unfavorable free energy at room temperature by T-WHAM analysis of replica exchange simulations was further tested on a system that can be solved analytically and presented some of the same challenges found in more complex chemical systems.
Related Papers
- → An adaptive umbrella sampling procedure in conformational analysis using molecular dynamics and its application to glycol(1992)149 cited
- → Na+–Na+ and Cl−–Cl− ion pairs in water: Mean force potentials by constrained molecular dynamics(1991)167 cited
- → Evaluation of Constrained and Restrained Molecular Dynamics Simulation Methods for Predicting Skin Lipid Permeability(2021)13 cited
- → A canonical replica exchange molecular dynamics implementation with normal pressure in each replica(2016)6 cited
- → Performance analysis of parallel replica molecular dynamics(2018)