Dataset Information

Determination of Ionic Hydration Free Energies with Grand Canonical Monte Carlo/Molecular Dynamics Simulations in Explicit Water.

ABSTRACT: Grand canonical Monte Carlo (GCMC) simulations of ionic solutions with explicit solvent models are known to be challenging. One challenge arises from the treatment of long-range electrostatics and finite-box size in Monte Carlo simulations when periodic boundary condition and Ewald summation methods are used. Another challenge is that constant excess chemical potential GCMC simulations for charged solutes suffer from inadequate insertion and deletion acceptance ratios. In this work, we address those problems by implementing an oscillating excess chemical potential GCMC algorithm with smooth particle mesh Ewald and finite-box-size corrections to treat the long-range electrostatics. The developed GCMC simulation program was combined with GROMACS to perform GCMC/MD simulations of ionic solutions individually containing Li+, Na+, K+, Rb+, Cs+, F-, Cl-, Br-, I-, Ca2+, and Mg2+, respectively. Our simulation results show that the combined GCMC/MD approach can approximate the ionic hydration free energies with proper treatment of long-range electrostatics. Our developed simulation approach can open up new avenues for simulating complex chemical and biomolecular systems and for drug discovery.

SUBMITTER: Sun D

PROVIDER: S-EPMC6195813 | biostudies-literature | 2018 Oct

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Determination of Ionic Hydration Free Energies with Grand Canonical Monte Carlo/Molecular Dynamics Simulations in Explicit Water.

Sun Delin D Lakkaraju Sirish Kaushik SK Jo Sunhwan S MacKerell Alexander D AD

Journal of chemical theory and computation 20180914 10

Grand canonical Monte Carlo (GCMC) simulations of ionic solutions with explicit solvent models are known to be challenging. One challenge arises from the treatment of long-range electrostatics and finite-box size in Monte Carlo simulations when periodic boundary condition and Ewald summation methods are used. Another challenge is that constant excess chemical potential GCMC simulations for charged solutes suffer from inadequate insertion and deletion acceptance ratios. In this work, we address t ...[more]

PMID: 30183291

Similar Datasets

Project description:The Grand Canonical Monte Carlo (GCMC) ensemble defined by the excess chemical potential, μex , volume, and temperature, in the context of molecular simulations allows for variations in the number of particles in the system. In practice, GCMC simulations have been widely applied for the sampling of rare gasses and water, but limited in the context of larger molecules. To overcome this limitation, the oscillating μex GCMC method was introduced and shown to be of utility for sampling small solutes, such as formamide, propane, and benzene, as well as for ionic species such as monocations, acetate, and methylammonium. However, the acceptance of GCMC insertions is low, and the method is computationally demanding. In the present study, we improved the sampling efficiency of the GCMC method using known cavity-bias and configurational-bias algorithms in the context of GPU architecture. Specifically, for GCMC simulations of aqueous solution systems, the configurational-bias algorithm was extended by applying system partitioning in conjunction with a random interval extraction algorithm, thereby improving the efficiency in a highly parallel computing environment. The method is parallelized on the GPU using CUDA and OpenCL, allowing for the code to run on both Nvidia and AMD GPUs, respectively. Notably, the method is particularly well suited for GPU computing as the large number of threads allows for simultaneous sampling of a large number of configurations during insertion attempts without additional computational overhead. In addition, the partitioning scheme allows for simultaneous insertion attempts for large systems, offering considerable efficiency. Calculations on the BK Channel, a transporter, including a lipid bilayer with over 760,000 atoms, show a speed up of ~53-fold through the use of system partitioning. The improved algorithm is then combined with an enhanced μex oscillation protocol and shown to be of utility in the context of the site-identification by ligand competitive saturation (SILCS) co-solvent sampling approach as illustrated through application to the protein CDK2.

Dataset Information

Determination of Ionic Hydration Free Energies with Grand Canonical Monte Carlo/Molecular Dynamics Simulations in Explicit Water.

Publications

Determination of Ionic Hydration Free Energies with Grand Canonical Monte Carlo/Molecular Dynamics Simulations in Explicit Water.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets