Biomolecular simulations under realistic macroscopic salt conditions Journal Article

Authors: Ross, G. A.; Rustenburg, A. S.; Grinaway, P. B.; Fass, J.; Chodera, J. D.
Article Title: Biomolecular simulations under realistic macroscopic salt conditions
Abstract: Biomolecular simulations are typically performed in an aqueous environment where the number of ions remains fixed for the duration of the simulation, generally with either a minimally neutralizing ion environment or a number of salt pairs intended to match the macroscopic salt concentration. In contrast, real biomolecules experience local ion environments where the salt concentration is dynamic and may differ from bulk. The degree of salt concentration variability and average deviation from the macroscopic concentration remains, as yet, unknown. Here, we describe the theory and implementation of a Monte Carlo osmostat that can be added to explicit solvent molecular dynamics or Monte Carlo simulations to sample from a semigrand canonical ensemble in which the number of salt pairs fluctuates dynamically during the simulation. The osmostat reproduces the correct equilibrium statistics for a simulation volume that can exchange ions with a large reservoir at a defined macroscopic salt concentration. To achieve useful Monte Carlo acceptance rates, the method makes use of nonequilibrium candidate Monte Carlo (NCMC) moves in which monovalent ions and water molecules are alchemically transmuted using short nonequilibrium trajectories, with a modified Metropolis-Hastings criterion ensuring correct equilibrium statistics for an (Delta mu, N, p, T) ensemble to achieve a similar to 10(46)x boost in acceptance rates. We demonstrate how typical protein (DHFR and the tyrosine kinase Src) and nucleic acid (Drew-Dickerson B-DNA dodecamer) systems exhibit salt concentration distributions that significantly differ from fixed-salt bulk simulations and display fluctuations that are on the same order of magnitude as the average.
Keywords: algorithm; thermodynamics; water; hydration; systems; efficient; dna-binding; monte-carlo simulations; molecular-dynamics simulations; free-energy
Journal Title: Journal of Physical Chemistry B
Volume: 122
Issue: 21
ISSN: 1520-6106
Publisher: American Chemical Society  
Date Published: 2018-05-31
Start Page: 5466
End Page: 5486
Language: English
ACCESSION: WOS:000434236900023
DOI: 10.1021/acs.jpcb.7b11734
PUBMED: 29649876
PMCID: PMC6078207
Notes: Article -- Source: Wos
Citation Impact
MSK Authors
  1. John Damon Chodera
    72 Chodera
  2. Joshua Hanna Fass
    7 Fass
  3. Gregory Antonio Ross
    5 Ross