A relatively simple protein solvent-accessible surface triangulation method for continuum electrostatics applications employing the boundary element method is presented. First, the protein is placed onto a three-dimensional lattice with a specified lattice spacing. To each lattice point, a box is assigned. Boxes located in the solvent region and in the interior of the protein are removed from the set. Improper connections between boxes and the possible existence of cavities in the interior of the protein which would destroy the proper connectivity of the triangulated surface are taken care of. The remaining set of boxes define the outer contour of the protein. Each free face exposed to the solvent of the remaining set of boxes is triangulated after the surface defined by the free faces has been smoothed to follow the shape of the macromolecule more accurately. The final step consists of a mapping of triangle vertices onto a set of surface points which define the solvent-accessible surface. Normal vectors at triangle vertices are obtained also from the free faces which define the orientation of the surface. The algorithm was tested for six molecules including four proteins; a dipeptide, a helical peptide consisting of 25 residues, calbindin, lysozyme, calmodulin and cutinase. For each molecule, total areas have been calculated and compared with the result computed from a dotted solvent-accessible surface. Since the boundary element method requires a low number of vertices and triangles to reduce the number of unknowns for reasons of efficiency, the number of triangles should not be too high. Nevertheless, credible results are obtained for the total area with relative errors not exceeding 12% for a large lattice spacing (0.30 nm) while close to zero for a smaller lattice spacing (down to 0.16 nm). The output of the triangulation computer program (written in C++) is rather simple so that it can be easily converted to any format acceptable for any molecular graphics programs.
Read full abstract