Abstract:A multi-scale framework was recently proposed for more realistic molecular dynamics simulations in continuum solvent models by coupling a molecular mechanics treatment of solute with a fluid mechanics treatment of solvent, where we formulated the physical model and developed a numerical fluid dynamics integrator. In this study, we incorporated the fluid dynamics integrator with the Amber simulation engine to conduct atomistic simulations of biomolecules. At this stage of the development, only nonelectrostatic interactions, i.e., van del Waals and hydrophobic interactions are included in the multi-scale model. Nevertheless numerical challenges exist in accurately interpolating the highly nonlinear van del Waals term when solving the finite-difference fluid dynamics equations. We were able to bypass the challenge rigorously by merging the van del Waals potential and pressure together when solving the fluid dynamics equations and by considering its contribution in the free-boundary condition analytically. The multi-scale simulation engine was first validated by reproducing the solute-solvent interface of a single atom with analytical solution. Next, we performed the relaxation simulation of a restrained symmetrical monomer and observed a symmetrical solvent interface at equilibrium with detailed surface features resembling those found on the solvent excluded surface. Four typical small molecular complexes were then tested, both volume and force balancing analysis showing that these simple complexes can reach equilibrium within the simulation time window. Finally, we studied the quality of the multi-scale solute-solvent interfaces for the four tested dimer complexes and found they agree well with the boundaries as sampled in the explicit water simulations.

A three-dimensional atomistic kinetic Monte Carlo study of dynamic solute-interface interaction

Stochastic model and kinetic Monte Carlo simulation of solute interactions with stationary and moving grain boundaries. I. Model formulation and application to one-dimensional systems

A Moving Front Kinetic Monte Carlo Algorithm for Moving Interface Systems

A Front-Tracking Method for Simulating Interfacial Flows with Particles and Soluble Surfactants

Kinetic Monte Carlo methods for three-dimensional diffusive capture problems in exterior domains

A kinetic Monte Carlo approach to study fluid transport in pore networks.

Multiscale modeling of solute diffusion in triblock copolymer membranes

Stochastic model and kinetic Monte Carlo simulation of solute interactions with stationary and moving grain boundaries. II. Application to two-dimensional systems

Convection rolls and three-dimensional particle dynamics in merging solute streams

Bridging the gap between mesoscopic and molecular models of solid/liquid interfaces out-of-equilibrium

A Semi-Empirical Force Balance-Based Model to Capture Sessile Droplet Spread on Smooth Surfaces: A Moving Front Kinetic Monte Carlo Study

A Unified, One Fluid Model for the Drag of Fluid and Solid Dispersals by Permeate Flux towards a Membrane Surface

Mesoscale modeling of near-contact interactions for complex flowing interfaces

The surface diffusivity of nanoparticles physically adsorbed at a solid-liquid interface

Modeling Mass Transfer and Reaction of Dilute Solutes in A Ternary Phase System by the Lattice Boltzmann Method

Exploring A Multi-Scale Method for Molecular Simulations in Continuum Solvent Model: Explicit Simulation of Continuum Solvent As An Incompressible Fluid

Mass transfer through vapour–liquid interfaces: a molecular dynamics simulation study

Bridging the scales in capillary rise dynamics with complexity-reduced models

Simulations of Interfacial Tension of Liquid-Liquid Ternary Mixtures Using Optimized Parametrization for Coarse-Grained Models