skip to main content

Title: A generalized Rayleigh-Plesset equation for ions with solvent fluctuations
We introduce a mathematical modeling framework for the conformational dynamics of charged molecules (i.e., solutes) in an aqueous solvent (i.e., water or salted water). The solvent is treated as an incompressible fluid, and its fluctuating motion is described by the Stokes equation with the Landau–Lifschitz stochastic stress. The motion of the solute-solvent interface (i.e., the dielectric boundary) is determined by the fluid velocity together with the balance of the viscous force,hydrostatic pressure, surface tension, solute-solvent van der Waals interaction force, and electrostatic force. The electrostatic interactions are described by the dielectric Poisson–Boltzmann theory.Within such a framework, we derive a generalized Rayleigh–Plesset equation, a nonlinear stochastic ordinary differential equation (SODE), for the radius of a spherical charged molecule, such as anion. The spherical average of the stochastic stress leads to a multiplicative noise. We design and test numerical methods for solving the SODE and use the equation, together with explicit solvent molecular dynamics simulations, to study the effective radius of a single ion. Potentially, our general modeling framework can be used to efficiently determine the solute-solvent interfacial structures and predict the free energies of more complex molecular systems.
; ;
Award ID(s):
Publication Date:
Journal Name:
SIAM journal on applied mathematics
Page Range or eLocation-ID:
Sponsoring Org:
National Science Foundation
More Like this
  1. Mechanical properties of cellular structures, including the cell cytoskeleton, are increasingly used as biomarkers for disease diagnosis and fundamental studies in cell biology. Recent experiments suggest that the cell cytoskeleton and its permeating cytosol, can be described as a poroelastic (PE) material. Biot theory is the standard model used to describe PE materials. Yet, this theory does not account for the fluid viscous stress, which can lead to inaccurate predictions of the mechanics in the dilute filamentous network of the cytoskeleton. Here, we adopt a two-phase model that extends Biot theory by including the fluid viscous stresses in the fluid's momentum equation. We use generalized linear viscoelastic (VE) constitutive equations to describe the permeating fluid and the network stresses and assume a constant friction coefficient that couples the fluid and network displacement fields. As the first step in developing a computational framework for solving the resulting equations, we derive closed-form general solutions of the fluid and network displacement fields in spherical coordinates. To demonstrate the applicability of our results, we study the motion of a rigid sphere moving under a constant force inside a PE medium, composed of a linear elastic network and a Newtonian fluid. We find that themore »network compressibility introduces a slow relaxation of the sphere and a non-monotonic network displacements with time along the direction of the applied force. These novel features cannot be predicted if VE constitutive equation is used for the medium. We show that our results can be applied to particle-tracking microrheology to differentiate between PE and VE materials and to independently measure the permeability and VE properties of the fluid and the network phases.« less
  2. Classical molecular dynamics simulations of the hydration thermodynamics, structure, and dynamics of water in hydration shells of charged buckminsterfullerenes are presented in this study. Charging of fullerenes leads to a structural transition in the hydration shell, accompanied by creation of a significant population of dangling O–H bonds pointing toward the solute. In contrast to the well accepted structure–function paradigm, this interfacial structural transition causes nearly no effect on either the dynamics of hydration water or on the solvation thermodynamics. Linear response to the solute charge is maintained despite significant structural changes in the hydration shell, and solvation thermodynamic potentials are nearly insensitive to the altering structure. Only solvation heat capacities, which are higher thermodynamic derivatives of the solvation free energy, indicate some sensitivity to the local hydration structure. We have separated the solvation thermodynamic potentials into direct solute–solvent interactions and restructuring of the hydration shell and analyzed the relative contributions of electrostatic and nonpolar interactions to the solvation thermodynamics.
  3. To help better interpret experimental measurement of nanoparticle size, it is important to understand how their diffusion depends on the physical and chemical features of surface ligands. In this study, explicit solvent molecular dynamics simulations are used to probe the effect of ligand charge and flexibility on the diffusion of small gold nanoparticles. The results suggest that despite a high bare charge (+18 e), cationic nanoparticles studied here have reduced diffusion constants compared to a hydrophobic gold nanoparticle by merely a modest amount. Increasing the ligand length by 10 CH2 units also has a limited impact on the diffusion constant. For the three particles studied here, the difference between estimated hydrodynamic radius and radius of gyration is on the order of one solvent layer (3–5 Å), confirming that the significant discrepancies found in the size of similar nanoparticles by recent transmission electron microscopy and dynamic light scattering measurements were due to aggregation under solution conditions. The limited impact of electrostatic friction on the diffusion of highly charged nanoparticles is found to be due to the strong anticorrelation between electrostatic and van der Waals forces between nanoparticle and environment, supporting the generality of recent observation for proteins by Matyushov and co-workers.more »Including the first shell of solvent molecules as part of the diffusing particle has a minor impact on the total force autocorrelation function but reduces the disparity in relaxation time between the total force and its electrostatic and van der Waals components.« less
  4. The Langevin Dynamics (LD) method (also known in the literature as Brownian Dynamics) is routinely used to simulate aerosol particle trajectories for transport rate constant calculations as well as to understand aerosol particle transport in internal and external fluid flows. This tutorial intends to explain the methodological details of setting up a LD simulation of a population of aerosol particles and to deduce rate constants from an ensemble of classical trajectories. We discuss the applicability and limitations of the translational Langevin equation to model the combined stochastic and deterministic motion of particles in fields of force or fluid flow. The drag force and stochastic “diffusion” force terms that appear in the Langevin equation are discussed elaborately, along with a summary of common forces relevant to aerosol systems (electrostatic, gravity, van der Waals, …); a commonly used first order and a fourth order Runge-Kutta time stepping schemes for linear stochastic ordinary differential equations are presented. A MATLAB® implementation of a LD code for simulating particle settling under gravity using the first order scheme is included for illustration. Scaling analysis of aerosol transport processes and the selection of timestep and domain size for trajectory simulations are demonstrated through two specific aerosol processes:more »particle diffusion charging and coagulation. Fortran® implementations of the first order and fourth order time-stepping schemes are included for simulating the 3D motion of a particle in a periodic domain. Potential applications and caveats to the usage of LD are included as a summary.« less
  5. We develop a hybrid approach that combines the Monte Carlo (MC)method, a variational implicit-solvent model (VISM), and a binary level-set method forthe simulation of biomolecular binding in an aqueous solvent. The solvation free energy for the biomolecular complex is estimated by minimizing the VISM free-energy functional of all possible solute−solvent interfaces that are used as dielectric boundaries. This functional consists of the solute volumetric, solute−solvent interfacial, solute−solvent van der Waals interaction, and electrostatic free energy. A technique of shifting the dielectric boundary is used to accurately predict the electrostatic part of the solvation free energy.Minimizing such a functional in each MC move is made possible by our new and fast binary level-set method. This method is based on the approximation of surface area by the convolution of an indicator function with a compactly supported kernel and is implemented by simple flips of numerical grid cells locally around the solute−solvent interface. We apply our approach to the p53-MDM2 system for which the two molecules are approximated by rigid bodies. Our efficient approach captures some of the poses before the final bound state. All atom molecular dynamics simulations with most of such poses quickly reach the final bound state.Our work is amore »new step toward realistic simulations of biomolecular interactions. With further improvement of coarse graining and MC sampling, and combined with other models, our hybrid approach can be used to study the free-energy landscape and kinetic pathways of ligand binding to proteins.« less