Hostname: page-component-7c8c6479df-27gpq Total loading time: 0 Render date: 2024-03-26T14:26:57.545Z Has data issue: false hasContentIssue false

Automated Parallel and Body-Fitted Mesh Generation in Finite Element Simulation of Macromolecular Systems

Published online by Cambridge University Press:  16 March 2016

Yan Xie
Affiliation:
State Key Laboratory of Scientific and Engineering Computing, Institute of Computational Mathematics and Scientific/Engineering Computing, Academy of Mathematics and Systems Science, Chinese Academy of Sciences, Beijing 100190, China
Tiantian Liu
Affiliation:
State Key Laboratory of Scientific and Engineering Computing, Institute of Computational Mathematics and Scientific/Engineering Computing, Academy of Mathematics and Systems Science, Chinese Academy of Sciences, Beijing 100190, China
Bin Tu
Affiliation:
National Center for NanoScience and Technology, Chinese Academy of Sciences, Beijing 100190, China
Benzhuo Lu*
Affiliation:
State Key Laboratory of Scientific and Engineering Computing, Institute of Computational Mathematics and Scientific/Engineering Computing, Academy of Mathematics and Systems Science, Chinese Academy of Sciences, Beijing 100190, China
Linbo Zhang*
Affiliation:
State Key Laboratory of Scientific and Engineering Computing, Institute of Computational Mathematics and Scientific/Engineering Computing, Academy of Mathematics and Systems Science, Chinese Academy of Sciences, Beijing 100190, China
*
*Corresponding author. Email addresses:xieyan@lsec.cc.ac.cn (Y. Xie), liutt@lsec.cc.ac.cn (T. Liu), tubin@lsec.cc.ac.cn (B. Tu), bzlu@lsec.cc.ac.cn (B. Lu), zlb@lsec.cc.ac.cn (L. Zhang)
*Corresponding author. Email addresses:xieyan@lsec.cc.ac.cn (Y. Xie), liutt@lsec.cc.ac.cn (T. Liu), tubin@lsec.cc.ac.cn (B. Tu), bzlu@lsec.cc.ac.cn (B. Lu), zlb@lsec.cc.ac.cn (L. Zhang)
Get access

Abstract

Mesh generation is a bottleneck for finite element simulations of biomolecules. A robust and efficient approach, based on the immersed boundary method proposed in [8], has been developed and implemented to generate large-scale mesh body-fitted to molecular shape for general parallel finite element simulations. The molecular Gaussian surface is adopted to represent the molecular surface, and is finally approximated by piecewise planes via the tool phgSurfaceCut in PHG [43], which is improved and can reliably handle complicated molecular surfaces, through mesh refinement steps. A coarse background mesh is imported first and then is distributed into each process using a mesh partitioning algorithm such as space filling curve [5] or METIS [22]. A bisection method is used for the mesh refinements according to the molecular PDB or PQR file which describes the biomolecular region. After mesh refinements, the mesh is optionally repartitioned and redistributed for load balancing. For finite element simulations, the modification of region mark and boundary types is done in parallel. Our parallel mesh generation method has been successfully applied to a sphere cavity model, a DNA fragment, a gramicidin A channel and a huge Dengue virus system. The results of numerical experiments show good parallel efficiency. Computations of electrostatic potential and solvation energy also validate the method. Moreover, the meshing process and adaptive finite element computation can be integrated as one PHG project to avoid the mesh importing and exporting costs, and improve the convenience of application as well.

Type
Research Article
Copyright
Copyright © Global-Science Press 2016 

Access options

Get access to the full version of this content by using one of the access options below. (Log in options will check for institutional or personal access. Content may require purchase if you do not have access.)

References

[1]Andersen, O.S.. Gramicidin channels. Annu. Rev. Physiol., 46:531548, 1984.CrossRefGoogle ScholarPubMed
[2]Bashford, D.. Scientific computing in object-oriented parallel environments. Lecture notes in computer science, 1343:233240, 1997.CrossRefGoogle Scholar
[3]Bates, P.W., Wei, G.W., and Zhao, S.. Minimal molecular surfaces and their applications r applications. Journal of Computational Chemistry, 29(3):380391, 2008.CrossRefGoogle Scholar
[4]Boschitsch, A.H., Fenley, M.O., and Zhou, H.. Fast boundary element method for the linear poisson-boltzmann equation. J. Phys. Chem. B., 106(10):27412754, 2002.CrossRefGoogle Scholar
[5]Campbell, P.M., Devine, K.D., Flaherty, J.E., Gervasio, L.G., and J.D. Teresco, . Dynamic octree load balancing using space-filling curves. Technical Report CS-03-01, 2003.Google Scholar
[6]Chen, M. and Lu, B.. Tmsmesh: a robust method for molecular surface mesh generation using a trace technique. J. Chem. Theory Comput., 7(1):203212, 2010.CrossRefGoogle ScholarPubMed
[7]Chen, M., Tu, B., and Lu, B.. Surface triangular mesh and volume tetrahedral mesh generations for biomolecular modeling. In Image-Based Geometric Modeling and Mesh Generation, pages 85106. Springer, 2013.CrossRefGoogle Scholar
[8]Chen, Z., Xiao, Y., and Zhang, L.. The adaptive immersed interface finite element method for elliptic and maxwell interface problems. J. Comput. Phys., 228(14):50005019, 2009.CrossRefGoogle Scholar
[9]Cheng, H.L. and Shi, X.. Quality mesh generation for molecular skin surfaces using restricted union of balls. Computational Geometry, 42(3):196206, 2009.CrossRefGoogle Scholar
[10]Chern, I.L., Liu, J.G., and Wang, W.C.. Accurate evaluation of electrostatics for macro-molecules in solution. Methods Appl. Anal., 10(2):309328, 2003.CrossRefGoogle Scholar
[11]Chew, L.P., Chrisochoides, N., and Sukup, F.. Parallel constrained delaunay meshing. ASME APPLIED MECHANICS DIVISION-PUBLICATIONS-AMD, 220:8996, 1997.Google Scholar
[12]Connolly, M.L.. Molecular surfaces: A review. Network Science, 14, 1996.Google Scholar
[13]De Cougny, H.L. and Shephard, M.S.. Parallel refinement and coarsening of tetrahedral meshes. International Journal for Numerical Methods in Engineering, 46(7):11011125, 1999.3.0.CO;2-E>CrossRefGoogle Scholar
[14]De Cougny, H.L., Shephard, M.S., and Ozturan, C.. Parallel three-dimensional mesh generation. Computing Systems in Engineering, 5(4):311323, 1994.CrossRefGoogle Scholar
[15]Dolinsky, T.J., Nielsen, J.E., McCammon, J.A., and Baker, N.A.. PDB2PQR: an automated pipeline for the setup, execution, and analysis of poisson-boltzmann electrostatics calculations. Nucleic. Acids. Res., 32:W665W667, 2004.CrossRefGoogle Scholar
[16]Edelsbrunner, H.. Deformable smooth surface design. Discrete & Computational Geometry, 21(1):87115, 1999.CrossRefGoogle Scholar
[17]Fang, Q.. Iso2mesh: a 3D surface and volumetric mesh generator for matlab/octave, 2010.Google Scholar
[18]Frey, P.J.. Medit: An interactive mesh visualization software. Technical Report RT-0253, 2001.Google Scholar
[19]Gerstein, M.Richards, F.M., Chapman, M.S., and Connolly, M.L.. Protein surfaces and volumes: measurement and use. In International Tables for Crystallography Volume F: Crystallography of biological macromolecules, pages 531545. Springer, 2001.Google Scholar
[20]Holst, M.J.. The poisson-boltzmann equation: analysis and multilevel numerical solution. 1994.Google Scholar
[21]Ito, Y., Shih, A.M., Erukala, A.K., Soni, B.K., Chernikov, A., Chrisochoides, N.P., and Naka-hashi, K.. Parallel unstructured mesh generation by an advancing front method. Mathematics and Computers in Simulation, 75(5):200209, 2007.CrossRefGoogle Scholar
[22]Karypis, G. and Kumar, V.. Metis-unstructured graph partitioning and sparse matrix ordering system, version 2.0. Citeseer, 1995.Google Scholar
[23]Koeppe, R.E. and Anderson, O.S.. Engineering the gramicidin channel. Annu. Rev. Cell. Dev. Biol., 25(1):231258, 1996.Google ScholarPubMed
[24]Kruithofand, N.G.H. and Vegter, G.. Meshing skin surfaces with certified topology. In Computer Aided Design and Computer Graphics, 2005. Ninth International Conference on, pages 6-pp. IEEE, 2005.Google Scholar
[25]Kuo, S.S., Altman, M.D., Bardhan, J.P., Tidor, B., and White, J.K.. Fast methods for simulation of biomolecule electrostatics. In Proceedings of the 2002 IEEE/ACM international conference on Computer-aided design, pages 466473. ACM, 2002.Google Scholar
[26]Löhner, R.. A parallel advancing front grid generation scheme. International Journal for Numerical Methods in Engineering, 51(6):663678, 2001.CrossRefGoogle Scholar
[27]Lu, B., Zhang, D., and McCammon, J.A.. Computation of electrostatic forces between sol-vated molecules determined by the poisson-boltzmann equation using a boundary element method. J. Chem. Phys., 122:214102, 2005.CrossRefGoogle ScholarPubMed
[28]Lu, B., Zhou, Y., Holst, M.J., and McCammon, J.A.. Recent progress in numerical methods for the Poisson-Boltzmann equation in biophysical applications. Commun. Comput. Phys., 3(5):9731009, 2008.Google Scholar
[29]Lu, B., Zhou, Y., Huber, G.A., Bond, S.D., Holst, M.J., and McCammon, J.A.. Electrodiffusion: A continuum modeling framework for biomolecular systems with realistic spatiotemporal resolution. J. Chem. Phys., 127(13):135102, 2007.CrossRefGoogle ScholarPubMed
[30]Lu, B., Holst, M.J., McCammon, J.A., and Zhou, Y.C.. Poisson-Nernst-Planck equations for simulating biomolecular diffusion-reaction processes I: Finite element solutions. J. Comput. Phys., 229(19):69796994, 2010.CrossRefGoogle ScholarPubMed
[31]Mehrotra, P., Saltz, J., and Voigt, R.. Unstructured scientific computation on scalable multiprocessors. MIT Press, 1992.Google Scholar
[32]Owen, S.J.. A survey of unstructured mesh generation technology. In IMR, pages 239267, 1998.Google Scholar
[33]Rocchia, W., Alexov, E., and Honig, B.. Extending the applicability of the nonlinear poisson-boltzmann equation: Multiple dielectric constants and multivalent ions. J. Phys. Chem. B., 105(28):65076514, 2001.CrossRefGoogle Scholar
[34]Said, R, Weatherill, N.P., Morgan, K., and Verhoeven, N. A.. Distributed parallel delaunay mesh generation. Computer methods in applied mechanics and engineering, 177(1):109125, 1999.CrossRefGoogle Scholar
[35]Shestakov, A.I., Milovich, J.L., and Noy, A.. Solution of the nonlinear poisson-boltzmann equation using pseudo-transient continuation and the finite element method. J. Colloid. Inter f. Set, 247(1):6279, 2002.CrossRefGoogle ScholarPubMed
[36]TetGen, H.Si., a Delaunay-based quality tetrahedral mesh generator. ACM Transactions on Mathematical Software, 41(2):Article 11, 2015.Google Scholar
[37]Tu, B., Chen, M., Xie, Y., Zhang, L., Eisenberg, B., and Lu, B.. A parallel finite element simulator for ion transport through three-dimensional ion channel systems. Journal of Computational Chemistry, 34(24):20652078, 2013.CrossRefGoogle ScholarPubMed
[38]Weiser, J., Shenkin, P.S., and Still, W.C.. Optimization of gaussian surface calculations and extension to solvent-accessible surface areas. Journal of computational chemistry, 20(7):688703, 1999.3.0.CO;2-F>CrossRefGoogle ScholarPubMed
[39]Xie, D. and Zhou, S.. A new minimization protocol for solving nonlinear Poisson-Boltzmann mortar finite element equation. BIT Numerical Mathematics, 47(4):853871, 2007.CrossRefGoogle Scholar
[40]Xie, Y., Cheng, J., Lu, B., and Zhang, L.. Parallel adaptive finite element algorithms for solving the coupled electro-diffusion equations. Molecular Based Mathematical Biology, 1:90108, 2013.Google Scholar
[41]Zaharescu, A., Boyer, E., and Horaud, R.P.. Transformesh: a topology-adaptive mesh-based approach to surface evolution. In Proceedings of the Eighth Asian Conference on Computer Vision, II:166175, November 2007.Google Scholar
[42]Zhang, B., Peng, B., Huang, J., Pitsianis, N.P., Sun, X., and Lu, B.. Parallel AFMPB solver with automatic surface meshing for calculation of molecular solvation free energy. Computer Physics Communications, 190:173181, 2015.CrossRefGoogle Scholar
[43]Zhang, L.. A parallel algorithm for adaptive local refinement of tetrahedral meshes using bisection. Numer. Math. Theor. Meth. Appl., 2(1):6589, 2009.Google Scholar
[44]Zhang, Y., Xu, G., and Bajaj, C.. Quality meshing of implicit solvation models of biomolecular structures. Computer Aided Geometric Design, 23(6):510530, 2006.CrossRefGoogle ScholarPubMed