A general and Robust Ray-Casting-Based Algorithm for Triangulating Surfaces at the Nanoscale

We present a general, robust, and efficient ray-casting-based approach to triangulating complex manifold surfaces arising in the nano-bioscience field. This feature is inserted in a more extended framework that: i) builds the molecular surface of nanometric systems according to several existing definitions, ii) can import external meshes, iii) performs accurate surface area estimation, iv) performs volume estimation, cavity detection, and conditional volume filling, and v) can color the points of a grid according to their locations with respect to the given surface. We implemented our methods in the publicly available NanoShaper software suite (www.electrostaticszone.eu). Robustness is achieved using the CGAL library and an ad hoc ray-casting technique. Our approach can deal with any manifold surface (including nonmolecular ones). Those explicitly treated here are the Connolly-Richards (SES), the Skin, and the Gaussian surfaces. Test results indicate that it is robust to rotation, scale, and atom displacement. This last aspect is evidenced by cavity detection of the highly symmetric structure of fullerene, which fails when attempted by MSMS and has problems in EDTSurf. In terms of timings, NanoShaper builds the Skin surface three times faster than the single threaded version in Lindow et al. on a 100,000 atoms protein and triangulates it at least ten times more rapidly than the Kruithof algorithm. NanoShaper was integrated with the DelPhi Poisson-Boltzmann equation solver. Its SES grid coloring outperformed the DelPhi counterpart. To test the viability of our method on large systems, we chose one of the biggest molecular structures in the Protein Data Bank, namely the 1VSZ entry, which corresponds to the human adenovirus (180,000 atoms after Hydrogen addition). We were able to triangulate the corresponding SES and Skin surfaces (6.2 and 7.0 million triangles, respectively, at a scale of 2 grids per Å) on a middle-range workstation.

[1]  Thomas Ertl,et al.  Parallel Contour-Buildup algorithm for the molecular surface , 2011, 2011 IEEE Symposium on Biological Data Visualization (BioVis)..

[2]  Amitabh Varshney,et al.  Parallel, stochastic measurement of molecular surface area. , 2008, Journal of molecular graphics & modelling.

[3]  Daniel Baum,et al.  Accelerated Visualization of Dynamic Molecular Surfaces , 2010, Comput. Graph. Forum.

[4]  Klaus Schulten,et al.  Fast Visualization of Gaussian Density Surfaces for Molecular Dynamics and Particle System Trajectories , 2012, EuroVis.

[5]  Yuan-Fang Wang,et al.  Efficient molecular surface generation using level-set methods. , 2006, Journal of molecular graphics & modelling.

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

[7]  Philipp Slusallek,et al.  Measuring properties of molecular surfaces using ray casting , 2010, 2010 IEEE International Symposium on Parallel & Distributed Processing, Workshops and Phd Forum (IPDPSW).

[8]  Emil Alexov,et al.  Rapid grid‐based construction of the molecular surface and the use of induced surface charge to calculate reaction field energies: Applications to the molecular systems and geometric objects , 2002, J. Comput. Chem..

[9]  Andrea Clematis,et al.  CUDA Accelerated Blobby Molecular Surface Generation , 2011, PPAM.

[10]  Yang Zhang,et al.  Generating Triangulated Macromolecular Surfaces by Euclidean Distance Transform , 2009, PloS one.

[11]  James F. Blinn,et al.  A Generalization of Algebraic Surface Drawing , 1982, TOGS.

[12]  Herbert Edelsbrunner,et al.  Triangulating the Surface of a Molecule , 1996, Discret. Appl. Math..

[13]  Tomas Akenine-Möller,et al.  Fast, Minimum Storage Ray-Triangle Intersection , 1997, J. Graphics, GPU, & Game Tools.

[14]  William E. Lorensen,et al.  Marching cubes: A high resolution 3D surface construction algorithm , 1987, SIGGRAPH.

[15]  Tomas Akenine-Möller,et al.  Fast, Minimum Storage Ray-Triangle Intersection , 1997, J. Graphics, GPU, & Game Tools.

[16]  Martin Falk,et al.  Interactive Exploration of Protein Cavities , 2011, Comput. Graph. Forum.

[17]  S. Sridharan,et al.  A rapid method for calculating derivatives of solvent accessible surface areas of molecules , 1995, J. Comput. Chem..

[18]  Mariette Yvinec,et al.  Triangulations in CGAL , 2002, Comput. Geom..

[19]  T. N. Bhat,et al.  The Protein Data Bank , 2000, Nucleic Acids Res..

[20]  Jianmin Zheng,et al.  Kernel modeling for molecular surfaces using a uniform solution , 2010, Comput. Aided Des..

[21]  Herbert Edelsbrunner,et al.  Deformable Smooth Surface Design , 1999, Discret. Comput. Geom..

[22]  J. A. Grant,et al.  A Gaussian Description of Molecular Shape , 1995 .

[23]  Ho-Lun Cheng,et al.  Dynamic Skin Triangulation , 2001, SODA '01.

[24]  M. Spagnuolo,et al.  Between algorithm and model: different Molecular Surface definitions for the Poisson-Boltzmann based electrostatic characterization of biomolecules in solution. , 2013, Communications in computational physics.

[25]  Matthieu Chavent,et al.  MetaMol: high-quality visualization of molecular skin surface. , 2008, Journal of molecular graphics & modelling.

[26]  S Subramaniam,et al.  Analytical shape computation of macromolecules: I. molecular area and volume through alpha shape , 1998, Proteins.

[27]  Byungjoo Kim,et al.  GPU Accelerated Molecular Surface Computing , 2012 .

[28]  Barry Honig,et al.  Extending the Applicability of the Nonlinear Poisson−Boltzmann Equation: Multiple Dielectric Constants and Multivalent Ions† , 2001 .

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

[30]  Thomas Ertl,et al.  Interactive Visualization of Molecular Surface Dynamics , 2009, IEEE Transactions on Visualization and Computer Graphics.

[31]  P. McAree,et al.  Using Sturm sequences to bracket real roots of polynomial equations , 1990 .

[32]  M. L. Connolly.  Analytical molecular surface calculation , 1983 .

[33]  Ho-Lun Cheng,et al.  Quality mesh generation for molecular skin surfaces using restricted union of balls , 2005, VIS 05. IEEE Visualization, 2005..

[34]  Tony You,et al.  An analytical algorithm for the rapid determination of the solvent accessibility of points in a three‐dimensional lattice around a solute molecule , 1995, J. Comput. Chem..

[35]  R. Mann,et al.  The role of DNA shape in protein-DNA recognition , 2009, Nature.

[36]  Gert Vegter,et al.  Meshing skin surfaces with certified topology , 2005, Ninth International Conference on Computer Aided Design and Computer Graphics (CAD-CG'05).

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

[38]  J Andrew McCammon,et al.  Feature-preserving adaptive mesh generation for molecular shape modeling and simulation. , 2008, Journal of molecular graphics & modelling.

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