Geometric and potential driving formation and evolution of biomolecular surfaces

This paper presents new geometrical flow equations for the theoretical modeling of biomolecular surfaces in the context of multiscale implicit solvent models. To account for the local variations near the biomolecular surfaces due to interactions between solvent molecules, and between solvent and solute molecules, we propose potential driven geometric flows, which balance the intrinsic geometric forces that would occur for a surface separating two homogeneous materials with the potential forces induced by the atomic interactions. Stochastic geometric flows are introduced to account for the random fluctuation and dissipation in density and pressure near the solvent–solute interface. Physical properties, such as free energy minimization (area decreasing) and incompressibility (volume preserving), are realized by some of our geometric flow equations. The proposed approach for geometric and potential forces driving the formation and evolution of biological surfaces is illustrated by extensive numerical experiments and compared with established minimal molecular surfaces and molecular surfaces. Local modification of biomolecular surfaces is demonstrated with potential driven geometric flows. High order geometric flows are also considered and tested in the present work for surface generation. Biomolecular surfaces generated by these approaches are typically free of geometric singularities. As the speed of surface generation is crucial to implicit solvent model based molecular dynamics, four numerical algorithms, a semi-implicit scheme, a Crank–Nicolson scheme, and two alternating direction implicit (ADI) schemes, are constructed and tested. Being either stable or conditionally stable but admitting a large critical time step size, these schemes overcome the stability constraint of the earlier forward Euler scheme. Aided with the Thomas algorithm, one of the ADI schemes is found to be very efficient as it balances the speed and accuracy.

[1]  G. Huisken.  The volume preserving mean curvature flow. , 1987 .

[2]  A. Bertozzi,et al.  Low‐curvature image simplifiers: Global regularity of smooth solutions and Laplacian limiting schemes , 2004 .

[3]  Q. Du,et al.  Modelling and simulations of multi-component lipid membranes and open membranes via diffuse interface approaches , 2006, Journal of Mathematical Biology.

[4]  J. Andrew Grant,et al.  A smooth permittivity function for Poisson–Boltzmann solvation methods , 2001, J. Comput. Chem..

[5]  D. Swigon,et al.  Mathematics of DNA structure, function and interactions , 2009 .

[6]  Shan Zhao,et al.  The minimal molecular surface , 2006, q-bio/0610038.

[7]  Richard A. Friesner,et al.  Numerical solution of the Poisson-Boltzmann equation using tetrahedral finite-element meshes , 1997, J. Comput. Chem..

[8]  J. Sethian,et al.  Subjective surfaces: a geometric model for boundary completion , 2000 .

[9]  G. W. Wei,et al.  Generalized Perona-Malik equation for image restoration , 1999, IEEE Signal Processing Letters.

[10]  Martin Rumpf,et al.  A finite element method for surface restoration with smooth boundary conditions , 2004, Comput. Aided Geom. Des..

[11]  Shin Yoshizawa,et al.  Fair triangle mesh generation with discrete elastica , 2002, Geometric Modeling and Processing. Theory and Applications. GMP 2002. Proceedings.

[12]  Guo-Wei Wei,et al.  Synchronization-based image edge detection , 2002 .

[13]  Michael L. Connolly,et al.  Molecular surface Triangulation , 1985 .

[14]  Leif Kobbelt,et al.  Generating fair meshes with G/sup 1/ boundary conditions , 2000, Proceedings Geometric Modeling and Processing 2000. Theory and Applications.

[15]  R. Levy,et al.  Computer simulations with explicit solvent: recent progress in the thermodynamic decomposition of free energies and in modeling electrostatic effects. , 1998, Annual review of physical chemistry.

[16]  L. Rudin,et al.  Feature-oriented image enhancement using shock filters , 1990 .

[17]  Arvid Lundervold,et al.  Noise removal using fourth-order partial differential equation with applications to medical magnetic resonance images in space and time , 2003, IEEE Trans. Image Process..

[18]  Peter Smereka,et al.  Semi-Implicit Level Set Methods for Curvature and Surface Diffusion Motion , 2003, J. Sci. Comput..

[19]  Kim A. Sharp,et al.  Electrostatic interactions in macromolecules , 1994 .

[20]  Ron Kimmel,et al.  A general framework for low level vision , 1998, IEEE Trans. Image Process..

[21]  M. Fixman,et al.  The Poisson–Boltzmann equation and its application to polyelectrolytes , 1979 .

[22]  B. Ninham,et al.  Specific ion effects in solutions of globular proteins: comparison between analytical models and simulation. , 2005, Journal of Physical Chemistry B.

[23]  D. M,et al.  A level set formulation for Willmore flow , 2004 .

[24]  K. Mikula,et al.  A direct method for solving an anisotropic mean curvature flow of plane curves with an external force , 2004 .

[25]  R. Zauhar,et al.  A new method for computing the macromolecular electric potential. , 1985, Journal of molecular biology.

[26]  Yuying Li,et al.  A computational algorithm for minimizing total variation in image restoration , 1996, IEEE Trans. Image Process..

[27]  S. Osher.  A level set formulation for the solution of the Dirichlet problem for Hamilton-Jacobi equations , 1993 .

[28]  J. Escher,et al.  The surface diffusion flow for immersed hypersurfaces , 1998 .

[29]  N. Sochen,et al.  Image Sharpening by Flows Based on Triple Well Potentials , 2004 .

[30]  Ray Luo,et al.  Implicit nonpolar solvent models. , 2007, The journal of physical chemistry. B.

[31]  Shulin Zhou,et al.  Existence and uniqueness of weak solutions for a fourth-order nonlinear parabolic equation , 2007 .

[32]  Harold A. Scheraga,et al.  A fast adaptive multigrid boundary element method for macromolecular electrostatic computations in a solvent , 1997, J. Comput. Chem..

[33]  Lin He,et al.  Incorporating topological derivatives into shape derivatives based level set methods , 2007, J. Comput. Phys..

[34]  Guillermo Sapiro,et al.  From active contours to anisotropic diffusion: connections between basic PDE's in image processing , 1996, Proceedings of 3rd IEEE International Conference on Image Processing.

[35]  J A McCammon,et al.  Coupling nonpolar and polar solvation free energies in implicit solvent models. , 2006, The Journal of chemical physics.

[36]  W. Im,et al.  Continuum solvation model: Computation of electrostatic forces from numerical solutions to the Poisson-Boltzmann equation , 1998 .

[37]  Chandrajit L. Bajaj,et al.  Quality meshing of implicit solvation models of biomolecular structures , 2006, Comput. Aided Geom. Des..

[38]  Douglas A. Lauffenburger,et al.  NUMERICAL SOLUTION OF THE NONLINEAR POISSON-BOLTZMANN EQUATION FOR A MEMBRANE-ELECTROLYTE SYSTEM , 1994 .

[39]  D. Mumford,et al.  Optimal approximations by piecewise smooth functions and associated variational problems , 1989 .

[40]  Jianhong Shen,et al.  A Stochastic-Variational Model for Soft Mumford-Shah Segmentation , 2005, Int. J. Biomed. Imaging.

[41]  F M Richards,et al.  Areas, volumes, packing and protein structure. , 1977, Annual review of biophysics and bioengineering.

[42]  Guillermo Sapiro,et al.  Geodesic Active Contours , 1995, International Journal of Computer Vision.

[43]  Ray Luo,et al.  A Poisson–Boltzmann dynamics method with nonperiodic boundary condition , 2003 .

[44]  H. Ashbaugh.  Convergence of Molecular and Macroscopic Continuum Descriptions of Ion Hydration , 2000 .

[45]  T. Chan,et al.  A Variational Level Set Approach to Multiphase Motion , 1996 .

[46]  Andreas Prohl,et al.  Analysis of a fully discrete finite element method for the phase field model and approximation of its sharp interface limits , 2003, Math. Comput..

[47]  A. Bertozzi,et al.  $H^1$ Solutions of a class of fourth order nonlinear equations for image processing , 2003 .

[48]  Li-Tien Cheng,et al.  Application of the level-set method to the implicit solvation of nonpolar molecules. , 2007, The Journal of chemical physics.

[49]  I. Holopainen.  Riemannian Geometry , 1927, Nature.

[50]  Alexander D. MacKerell,et al.  All-atom empirical potential for molecular modeling and dynamics studies of proteins. , 1998, The journal of physical chemistry. B.

[51]  Olivier D. Faugeras,et al.  Using the Vector Distance Functions to Evolve Manifolds of Arbitrary Codimension , 2001, Scale-Space.

[52]  Tony F. Chan,et al.  High-Order Total Variation-Based Image Restoration , 2000, SIAM J. Sci. Comput..

[53]  C. Brooks,et al.  Recent advances in the development and application of implicit solvent models in biomolecule simulations. , 2004, Current opinion in structural biology.

[54]  Yang Wang,et al.  Total Variation Wavelet-Based Medical Image Denoising , 2006, Int. J. Biomed. Imaging.

[55]  M. Droske,et al.  A level set formulation for Willmore flow , 2004 .

[56]  J. Warwicker,et al.  Calculation of the electric potential in the active site cleft due to alpha-helix dipoles. , 1982, Journal of molecular biology.

[57]  Richard H. Henchman,et al.  Revisiting free energy calculations: a theoretical connection to MM/PBSA and direct calculation of the association free energy. , 2004, Biophysical journal.

[58]  J. Sethian.  Evolution, implementation, and application of level set and fast marching methods for advancing fronts , 2001 .

[59]  Frank Eisenhaber,et al.  Improved strategy in analytic surface calculation for molecular systems: Handling of singularities and computational efficiency , 1993, J. Comput. Chem..

[60]  H. Lawson.  Lectures On Minimal Submanifolds , 1980 .

[61]  Thomas C. Cecil,et al.  A numerical method for computing minimal surfaces in arbitrary dimension , 2005 .

[62]  G. Wei,et al.  Molecular multiresolution surfaces , 2005, math-ph/0511001.

[63]  Charles L. Brooks,et al.  Performance comparison of generalized born and Poisson methods in the calculation of electrostatic solvation energies for protein structures , 2004, J. Comput. Chem..

[64]  L. Rudin,et al.  Nonlinear total variation based noise removal algorithms , 1992 .

[65]  Nathan A. Baker,et al.  Solvent reaction field potential inside an uncharged globular protein: a bridge between implicit and explicit solvent models? , 2007, The Journal of chemical physics.

[66]  S. Osher,et al.  Level set methods: an overview and some recent results , 2001 .

[67]  M. Gage,et al.  The heat equation shrinking convex plane curves , 1986 .

[68]  Kris Sikorski.  Handbook of Numerical Analysis. Volume 1: Finite Difference Methods (Part 1) and Solution of Equations in Rn (Part 1) (P. G. Ciarlet and J L. Lions, eds.) , 1992, SIAM Rev..

[69]  Matthew P Jacobson,et al.  Surfaces affect ion pairing. , 2005, The journal of physical chemistry. B.

[70]  M. Sanner,et al.  Reduced surface: an efficient way to compute molecular surfaces. , 1996, Biopolymers.

[71]  Kenneth A. Brakke,et al.  The motion of a surface by its mean curvature , 2015 .

[72]  B. Honig,et al.  Classical electrostatics in biology and chemistry. , 1995, Science.

[73]  Gieri Simonett,et al.  The Willmore flow near spheres , 2001, Differential and Integral Equations.

[74]  Ge Wang,et al.  Evolution-Operator-Based Single-Step Method for Image Processing , 2006, Int. J. Biomed. Imaging.

[75]  A Bultheel,et al.  Book review 'Handbook of numerical analysis, volume I. (Finite difference methods - part 1; Solution of equations in R^n - part 1) P.G. Ciarlet and J.L. Lions (eds.)' , 1990 .

[76]  Andres Fco. Solé,et al.  3D Curves Reconstruction Based on Deformable Models , 2003, Journal of Mathematical Imaging and Vision.

[77]  Guillermo Sapiro,et al.  Anisotropic diffusion of multivalued images with applications to color filtering , 1996, IEEE Trans. Image Process..

[78]  Jan Forsman,et al.  A Simple Correlation-Corrected Poisson-Boltzmann Theory , 2004 .

[79]  Weihua Geng,et al.  Treatment of geometric singularities in implicit solvent models. , 2007, The Journal of chemical physics.

[80]  Nathan A. Baker,et al.  Improving implicit solvent simulations: a Poisson-centric view. , 2005, Current opinion in structural biology.

[81]  J. Sethian,et al.  Fronts propagating with curvature-dependent speed: algorithms based on Hamilton-Jacobi formulations , 1988 .

[82]  Leif Kobbelt,et al.  Geometric fairing of irregular meshes for free-form surface design , 2001, Comput. Aided Geom. Des..

[83]  B. Lee,et al.  The interpretation of protein structures: estimation of static accessibility. , 1971, Journal of molecular biology.

[84]  S. Osher,et al.  Algorithms Based on Hamilton-Jacobi Formulations , 1988 .

[85]  Mostafa Kaveh,et al.  Fourth-order partial differential equations for noise removal , 2000, IEEE Trans. Image Process..

[86]  Stanley Osher,et al.  Shocks and other nonlinear filtering applied to image processing , 1991, Optics & Photonics.

[87]  Shan Zhao,et al.  Minimal molecular surfaces and their applications , 2008, J. Comput. Chem..

[88]  Jitendra Malik,et al.  Scale-Space and Edge Detection Using Anisotropic Diffusion , 1990, IEEE Trans. Pattern Anal. Mach. Intell..

[89]  D. Chopp.  Computing Minimal Surfaces via Level Set Curvature Flow , 1993 .

[90]  B. Roux,et al.  Implicit solvent models. , 1999, Biophysical chemistry.

[91]  L. Kobbelt,et al.  Generating Fair Meshes with $G^1$ Boundary Conditions , 2000 .

[92]  Nathan A. Baker,et al.  Solvation forces on biomolecular structures: A comparison of explicit solvent and Poisson–Boltzmann models , 2004, J. Comput. Chem..

[93]  Tony F. Chan,et al.  Color TV: total variation methods for restoration of vector-valued images , 1998, IEEE Trans. Image Process..

[94]  Valentin Gogonea,et al.  Implementation of the solvent effect in molecular mechanics. 1. Model development and analytical algorithm for the solvent-accessible surface area , 1994 .