Qualified, stable, and efficient molecular surface meshing appears to be necessitated by recent developments for realistic mathematical modeling and numerical simulation of biomolecules, especially in implicit solvent modeling (e.g., see a review in B. Z. Lu et al. Commun. Comput. Phys. 2008, 3, 973-1009). In this paper, we present a new method: tracing molecular surface for meshing (TMSmesh) the Gaussian surface of biomolecules. The method computes the surface points by solving a nonlinear equation directly, polygonizes by connecting surface points through a trace technique, and finally outputs a triangulated mesh. TMSmesh has a linear complexity with respect to the number of atoms and is shown to be capable of handling molecules consisting of more than one million atoms, which is usually difficult for the existing methods for surface generation used in molecular visualization and geometry analysis. Moreover, the meshes generated by TMSmesh are successfully tested in boundary element solutions of the Poisson-Boltzmann equation, which directly gives rise to a route to simulate electrostatic solvation of large-scale molecular systems. The binary version of TMSmesh and a set of representative PQR benchmark molecules are downloadable at our Web page http://lsec.cc.ac.cn/∼lubz/Meshing.html .

A parallel finite element simulator, ichannel, is developed for ion transport through three-dimensional ion channel systems that consist of protein and membrane. The coordinates of heavy atoms of the protein are taken from the Protein Data Bank and the membrane is represented as a slab. The simulator contains two components: a parallel adaptive finite element solver for a set of Poisson-Nernst-Planck (PNP) equations that describe the electrodiffusion process of ion transport, and a mesh generation tool chain for ion channel systems, which is an essential component for the finite element computations. The finite element method has advantages in modeling irregular geometries and complex boundary conditions. We have built a tool chain to get the surface and volume mesh for ion channel systems, which consists of a set of mesh generation tools. The adaptive finite element solver in our simulator is implemented using the parallel adaptive finite element package Parallel Hierarchical Grid (PHG) developed by one of the authors, which provides the capability of doing large scale parallel computations with high parallel efficiency and the flexibility of choosing high order elements to achieve high order accuracy. The simulator is applied to a real transmembrane protein, the gramicidin A (gA) channel protein, to calculate the electrostatic potential, ion concentrations and I - V curve, with which both primitive and transformed PNP equations are studied and their numerical performances are compared. To further validate the method, we also apply the simulator to two other ion channel systems, the voltage dependent anion channel (VDAC) and α-Hemolysin (α-HL). The simulation results agree well with Brownian dynamics (BD) simulation results and experimental results. Moreover, because ionic finite size effects can be included in PNP model now, we also perform simulations using a size-modified PNP (SMPNP) model on VDAC and α-HL. It is shown that the size effects in SMPNP can effectively lead to reduced current in the channel, and the results are closer to BD simulation results.

We present a numerical method to approximate the memory functions in the generalized Langevin models for the collective dynamics of macromolecules. We first derive the exact expressions of the memory functions, obtained from projection to subspaces that correspond to the selection of coarse-grain variables. In particular, the memory functions are expressed in the forms of matrix functions, which will then be approximated by Krylov-subspace methods. It will also be demonstrated that the random noise can be approximated under the same framework, and the second fluctuation-dissipation theorem is automatically satisfied. The accuracy of the method is examined through several numerical examples.

No abstract