Abstract: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.

Efficient Sampling of Cavity Hydration in Proteins with Nonequilibrium Grand Canonical Monte Carlo and Polarizable Force Fields

GPU-specific algorithms for improved solute sampling in grand canonical Monte Carlo simulations

Monte Carlo Sampling of Protein Folding by Combining an All-Atom Physics-Based Model with a Native State Bias

Sampling of Organic Solutes in Aqueous and Heterogeneous Environments Using Oscillating Excess Chemical Potentials in Grand Canonical-like Monte Carlo-Molecular Dynamics Simulations

Quantum Chemistry for Solvated Molecules on Graphical Processing Units (GPUs)using Polarizable Continuum Models

A quantum mechanical computational method for modeling electrostatic and solvation effects of protein

Quantum Fragment Based ab Initio Molecular Dynamics for Proteins.

A combined steepest descent and genetic algorithm (SD/GA) approach for the optimization of solvation parameters

A library-based Monte Carlo technique enables rapid equilibrium sampling of a protein model with atomistic components

Development of an Effective Polarizable Bond Method for Biomolecular Simulation.

A Force Balanced Fragmentation Method for Ab Initio Molecular Dynamic Simulation of Protein

Statistical Analysis, Investigation, and Prediction of the Water Positions in the Binding Sites of Proteins

Iterative derivation of effective potentials to sample the conformational space of proteins at atomistic scale

Constrained Hybrid Monte Carlo Sampling Made Simple for Chemical Reaction Simulations

Efficient Sampling over Rough Energy Landscapes with High Barriers: A Combination of Metadynamics with Integrated Tempering Sampling.

Benchmarking polarizable and non-polarizable force fields for Ca2+-peptides against a comprehensive QM dataset

Monte Carlo simulation of electrostatic interactions in inhomogeneous dielectric media: Correct sampling for the local lattice simulation algorithm

On delivering polar solvation free energy of proteins from energy minimized structures using a regularized super‐Gaussian Poisson–Boltzmann model

Grand-canonical molecular dynamics simulations powered by a hybrid 4D nonequilibrium MD/MC method: Implementation in LAMMPS and applications to electrolyte solutions

CONDUCTOR, a fast, efficient and robust tool for continuum solvation computations for macromolecules

Open Boundary Simulations of Proteins and Their Hydration Shells by Hamiltonian Adaptive Resolution Scheme