Conformational populations of ligand‐sized molecules by replica exchange molecular dynamics and temperature reweighting
- 30 October 2009
- journal article
- research article
- Published by Wiley in Journal of Computational Chemistry
- Vol. 31 (7), 1357-1367
- https://doi.org/10.1002/jcc.21419
Abstract
The use of the replica exchange (RE) molecular dynamics (MD) method for the efficient estimation of conformational populations of ligand-sized molecules in solution is investigated. We compare the computational efficiency of the traditional constant temperature MD technique with that of the parallel RE molecular dynamics method for a series of alkanes and rilpivirine (TMC278), an inhibitor against HIV-1 reverse transcriptase, with implicit solvation. We show that conformational populations are accurately estimated by both methods; however, replica exchange estimates converge at a faster rate, especially for rilpivirine, which is characterized by multiple stable states separated by high-free energy barriers. Furthermore, convergence is enhanced when the weighted histogram analysis method (WHAM) is used to estimate populations from the data collected from multiple RE temperature replicas. For small drug-like molecules with energetic barriers separating the stable states, the use of RE with WHAM is an efficient computational approach for estimating the contribution of ligand conformational reorganization to binding affinities. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2010Keywords
This publication has 39 references indexed in Scilit:
- High-resolution structures of HIV-1 reverse transcriptase/TMC278 complexes: Strategic flexibility explains potency against resistance mutationsProceedings of the National Academy of Sciences of the United States of America, 2008
- Integrated Modeling Program, Applied Chemical Theory (IMPACT)Journal of Computational Chemistry, 2005
- Generalized ensemble simulations for complex systemsComputer Physics Communications, 2002
- Generalized-ensemble algorithms for molecular simulations of biopolymersPeptide Science, 2001
- Replica-exchange molecular dynamics method for protein foldingChemical Physics Letters, 1999
- Multicanonical Ensemble Generated by Molecular Dynamics Simulation for Enhanced Conformational Sampling of PeptidesThe Journal of Physical Chemistry B, 1997
- Molecular dynamics, Langevin and hydrid Monte Carlo simulations in a multicanonical ensembleChemical Physics Letters, 1996
- Exchange Monte Carlo Method and Application to Spin Glass SimulationsJournal of the Physics Society Japan, 1996
- Multicanonical ensemble: A new approach to simulate first-order phase transitionsPhysical Review Letters, 1992
- Multicanonical algorithms for first order phase transitionsPhysics Letters B, 1991