RAPID APPROXIMATION TO MOLECULAR-SURFACE AREA VIA THE USE OF BOOLEAN LOGIC AND LOOK-UP TABLES

We report the development of a new approximate method of calculating molecular surface areas. Our technique is based upon the method of Shrake and Rupley but incorporates several major advances. First, we represent the state of surface points as bits in a bit string so we can utilize Boolean operations to simultaneously turn off multiple test points in one Boolean AND operation. Second, we use a series of Boolean mask look-up tables to reduce the time complexity of the calculation of molecular surface area down to the same magnitude as doing a potential energy evaluation. When we use a 256 surface point sphere for all of the atoms in BPTI, a 454 nonhydrogen atom protein, and a 1.4-angstrom solvent probe, we in general underestimate the total solvent-accessible surface area (SASA) by approximately 1.25% with a correlation coefficient of 0.9990 over a wide range of conformations. The average CPU time required to calculate the SASA of a BPTI conformer is 0.58 s on an SGI 4D/220 workstation. We also describe a method by which we can calculate an approximate finite difference SASA gradient for BPTI in 0.79 of CPU time.