Quantitative analysis of molecular surface based on improved Marching Tetrahedra algorithm.

Quantitative analysis of molecular surface is a valuable technique for analyzing non-covalent interaction, studying molecular recognition mode, predicting reactive site and reactivity. An efficient way to realize the analysis was first proposed by Bulat et al. (J. Mol. Model., 16, 1679), in which Marching Tetrahedra (MT) approach commonly used in computer graphics is employed to generate vertices on molecular surface. However, it has been found that the computations of the electrostatic potential in the MT vertices are very expensive and some artificial surface extremes will be presented due to the uneven distribution of MT vertices. In this article, we propose a simple and reliable method to eliminate these unreasonably distributed surface vertices generated in the original MT. This treatment can save more than 60% of total analysis time of electrostatic potential, yet the loss in accuracy is almost negligible. The artificial surface extremes are also largely avoided as a byproduct of this algorithm. In addition, the bisection iteration procedure has been exploited to improve accuracy of linear interpolation in MT. The most appropriate grid spacing for surface analysis has also been investigated. 0.25 and 0.20 bohr are recommended to be used for surface analysis of electrostatic potential and average local ionization energy, respectively.

[1]  Peter Politzer,et al.  Average local ionization energies on the molecular surfaces of aromatic systems as guides to chemical reactivity , 1990 .

[2]  J. Murray,et al.  Electrostatic Potentials: Chemical Applications , 2002 .

[3]  Peter Politzer,et al.  σ‐Hole bonding and hydrogen bonding: Competitive interactions , 2007 .

[4]  William E. Lorensen,et al.  Decimation of triangle meshes , 1992, SIGGRAPH.

[5]  Enrico O. Purisima,et al.  Molecular surface generation using marching tetrahedra , 1998, J. Comput. Chem..

[6]  Peter Politzer,et al.  An overview of halogen bonding , 2007, Journal of molecular modeling.

[7]  J. Murray,et al.  Average local ionization energy: A review , 2010, Journal of molecular modeling.

[8]  Kok-Lim Low,et al.  Model simplification using vertex-clustering , 1997, SI3D.

[9]  Peter Politzer,et al.  The electrostatic potential: an overview , 2011 .

[10]  Enrico O. Purisima,et al.  Molecular surface generation using marching tetrahedra , 1998, J. Comput. Chem..

[11]  J. Murray,et al.  Statistically-based interaction indices derived from molecular surface electrostatic potentials: a general interaction properties function (GIPF) , 1994 .

[12]  J. Pople,et al.  Self—Consistent Molecular Orbital Methods. XII. Further Extensions of Gaussian—Type Basis Sets for Use in Molecular Orbital Studies of Organic Molecules , 1972 .

[13]  T. Todd Elvins,et al.  A survey of algorithms for volume visualization , 1992, COMG.

[14]  Hong Yi,et al.  A survey of the marching cubes algorithm , 2006, Comput. Graph..

[15]  Enrico O. Purisima,et al.  A new tetrahedral tesselation scheme for isosurface generation , 1998, Comput. Graph..

[16]  Thomas Lewiner,et al.  Efficient Implementation of Marching Cubes' Cases with Topological Guarantees , 2003, J. Graphics, GPU, & Game Tools.

[17]  Cheng Chang,et al.  Properties of atoms in molecules: atomic volumes , 1987 .

[18]  Kevin E. Riley,et al.  σ-Holes, π-holes and electrostatically-driven interactions , 2012, Journal of Molecular Modeling.

[19]  Andrew H. Gee,et al.  Regularised marching tetrahedra: improved iso-surface extraction , 1999, Comput. Graph..

[20]  P. C. Hariharan,et al.  The influence of polarization functions on molecular orbital hydrogenation energies , 1973 .

[21]  W E Bolch,et al.  Marching cube algorithm: review and trilinear interpolation adaptation for image-based dosimetric models. , 2003, Computerized medical imaging and graphics : the official journal of the Computerized Medical Imaging Society.

[22]  P. Schleyer.  Encyclopedia of computational chemistry , 1998 .

[23]  Peter Politzer,et al.  The fundamental nature and role of the electrostatic potential in atoms and molecules , 2002 .

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

[25]  Peter Politzer,et al.  Molecular Electrostatic Potentials and Chemical Reactivity , 2007 .

[26]  A. Becke.  Density-functional thermochemistry. III. The role of exact exchange , 1993 .

[27]  Peter Politzer,et al.  Quantitative analysis of molecular surfaces: areas, volumes, electrostatic potentials and average local ionization energies , 2010, Journal of molecular modeling.

[28]  Tian Lu,et al.  Multiwfn: A multifunctional wavefunction analyzer , 2012, J. Comput. Chem..

[29]  P. Politzer,et al.  Molecular electrostatic potentials: an effective tool for the elucidation of biochemical phenomena. , 1985, Environmental health perspectives.

[30]  Timothy Clark,et al.  Halogen bonding: the σ-hole , 2007 .

[31]  Paolo Cignoni,et al.  A comparison of mesh simplification algorithms , 1998, Comput. Graph..