The interface between protein solute and aqueous solvent exhibits complex geometries, and can undergo conformational changes by combined influences from electrostatic force, surface tension, and hydrodynamic force. To understand the role of solvent Stokes flow in this process, we develop a Dynamic Implicit Solvent Model (DISM). Based on this model, we first analytically study the linear stability of a cylindrical solute-solvent interface, where the asymptotic dispersion reveals a power law. Moreover, we develop a computational method to simulate the solvent Stokes flow and interface motion. The key components of our fluid solver are virtual node method, pressure Poisson equation, specially designed boundary condition, Schur complement and least square techniques. 3D numerical tests will demonstrate the convergence of our method, and show that this new approach can capture interesting dynamics of protein conformational change.

海报