the interactions of biological macromolecules are critical to their physiological function and dependent on their molecular interaction fields . the electrostatic potential is a molecular interaction field of particular importance for determining the specificity and kinetics of molecular binding . pipsa ( protein interaction property similarity analysis ) ( 1 ) may be used to classify a large number of proteins according to the similarities or dissimilarities in their 3d molecular interaction fields , such as the electrostatic potential ( 2,3 ) . some of the applications of pipsa have been to ww domains ( 4 ) , electron transfer proteins ( 5 ) , ubiquitin conjugating enzymes ( 6 ) , complement control protein modules ( 7 ) and dihydrofolate reductases ( 8) . an extension of pipsa is qpipsa ( quantitative pipsa ) ( 9 ) . qpipsa enables the kinetic parameters of a set of enzymes sharing the same function to be related to the molecular interaction field , e.g. the electrostatic potential , at a functional region of the protein . such a comparison may enable estimation of unknown kinetic parameters for an enzymatic reaction and thereby assist in the modeling and simulation of biochemical pathways ( 9,10 ) . webpipsa allows the user to perform pipsa to compute and compare the electrostatic potentials of a set of structurally related proteins . other web servers such as pce ( http://bioserv.rpbs.jussieu.fr/pce ) ( 11 ) and pfplus ( http://pfp.technion.ac.il ) ( 12 ) also allow the calculation of protein electrostatic potentials . pfplus was designed to extract and display the largest positive electrostatic patch on a protein surface . these web servers , however , only allow the calculation of electrostatic potentials for single proteins . webpipsa on the other hand permits calculation of the electrostatic potentials of a large number of proteins and performs an all - versus - all pairwise comparison of the electrostatic potentials around the entire protein and , optionally , over a user - predefined region . the electrostatic potentials are then computed on a grid by solution of the linearized finite - difference poisson the similarity or dissimilarity of the electrostatic potential of each pair of proteins in the dataset is quantified for a user - defined region by means of similarity indices and distance measures . the results are displayed as an epogram ( tree - like diagram based on electrostatic differences ) and a colored matrix view . the results can be used for the classification of the proteins according to their interaction properties or for the correlation of the molecular interaction fields with functional properties of the proteins . the input required for using webpipsa is a set of protein coordinate files in pdb format . these can be either ( i ) user - supplied or ( ii ) specified by pdb identifier code in the rcsb . the user - supplied structures can be either experimentally determined or generated by comparative modeling techniques using , for example , modeller ( 15 ) or swissmodel ( 16 ) or structures from the respective databases , such as modbase ( 17 ) . an example dataset consisting of triosephosphate isomerases from various species ( 9 ) is provided on the web server and described in detail below . after the structures have been uploaded , the user can choose whether to use the uhbd ( 13 ) or the apbs ( 14 ) software to calculate the electrostatic potentials . example input files for uhbd and apbs are given in the online documentation for webpipsa . at this point , the user is asked to provide an email address to be notified when the pipsa calculation has been completed . the user receives an email with a link to the results page displaying the heat map and a clustered epogram . the email also provides a link to the results directory on the web server containing intermediate and additional results such as the potential grid files and the complete pipsa output . a description of how to visualize electrostatic potential grids from the results directory is given in the online documentation . triosephosphate isomerases from 12 different species are clustered according to the all pairwise distances between electrostatic potentials using a color code from red ( small distance ) to blue ( large distance ) . proteins with similar electrostatic properties are clustered on the top and left side of the graph . the proteins are assigned to a user - defined number of clusters according to their electrostatic potential similarities . the number of clusters can be specified by selecting from a drop down box at the bottom of the result page . here , triosephosphate isomerases from 12 species form four subclusters in a comparison of the electrostatic potentials around the whole protein surface . the type of distance measure [ hodgkin similarity index ( 18 ) , carbo similarity index or average potential differences ( 9 ) ] displayed in the heat map can also be selected from a dropdown box . the distance matrix may also be viewed with the proteins in their input order , without clustering . figure 1.screenshot of part of the results page for the example provided on the webpipsa web server for the comparison of the electrostatic potentials of triosephosphate isomerases from 12 different species . the colored matrix ( heat map ) is shown with color coding corresponding to the distances calculated from the hodgkin similarity indices for the electrostatic potentials . enzymes from species with highly similar electrostatic potentials , such as human , rabbit and chicken , are clustered together ( second subcluster from top / left ) . screenshot of part of the results page for the example provided on the webpipsa web server for the comparison of the electrostatic potentials of triosephosphate isomerases from 12 different species . the colored matrix ( heat map ) is shown with color coding corresponding to the distances calculated from the hodgkin similarity indices for the electrostatic potentials . enzymes from species with highly similar electrostatic potentials , such as human , rabbit and chicken , are clustered together ( second subcluster from top / left ) . after the upload of the protein structures is complete , the user can choose to superimpose the structures using one of two different methods . sup2pdb , starts with one structure ( the template ) and subsequently does a pairwise sequence alignment ( 19 ) between this template structure and the remaining coordinate files . the structures are then superimposed according to the respective alignments . since the outcome of this algorithm depends on the choice of protein template , in the second option , optimized sup2pdb , every coordinate file is considered as a template for superimposing the others . then the most successful superposition based on the maximum number of matched structures and minimal root mean square deviation ( rmsd ) of the superimposed structures is selected . for a large number of protein structures , however , this second option can be time consuming . the protein structures are a subset of the triosephosphate isomerases in the example analysis provided on the web server . the region for comparison of the electrostatic potential can be selected after upload of the coordinate files ( the substrate binding site in the example given ) . these distances are used to cluster the proteins according to the relations between their electrostatic potentials and the clustering is displayed in a tree - like diagram , an epogram . the protein structures are a subset of the triosephosphate isomerases in the example analysis provided on the web server . the region for comparison of the electrostatic potential can be selected after upload of the coordinate files ( the substrate binding site in the example given ) . these distances are used to cluster the proteins according to the relations between their electrostatic potentials and the clustering is displayed in a tree - like diagram , an epogram . the next step in the workflow is to add polar hydrogen atoms to the protein structures using whatif ( 20 ) . standard protonation states at ph 7 are assumed for all residues except for histidine which can be treated as singly or doubly protonated . for the superimposed set of coordinate files , electrostatic potentials are computed automatically using a set of default parameters assuming an ionic strength of 50 mm and a temperature of 300 k. uhbd or apbs can be used to solve the linearized poisson boltzmann equation ( lpbe ) treating the protein as a low dielectric with partial atomic charges embedded in a homogeneous high dielectric continuum representing the solvent . additionally , the user can specify spherical regions within which the electrostatic potentials in the skin are compared . the cartesian coordinates of the center of the sphere and its radius should be input by the user . this option can therefore only be used when the uploaded input structures have already been superimposed . the pipsa software is used to calculate hodgkin and carbo similarity indices of the protein electrostatic potentials as well as average electrostatic potential differences ( the difference in electrostatic potentials in kcal mol e of two proteins divided by the number of grid points in the comparison region where the two protein skins overlap ) described in ( 9 ) . the similarity indices range from 1 ( anti - correlated potential ) , through 0 ( uncorrelated ) to + 1 ( identical potentials ) . these values are converted into distances given by sqrt(2  2*si ) where si is the respective similarity index . these distances thus range from 0 ( identical ) to 2 ( anti - correlated potentials ) . the clustering analysis and generation of tree - like diagrams ( epograms ) on the results page are performed using the statistical program r ( 21 ) . for visualization purposes , both the heat map and epogram representations allow the fast identification of inter - protein relations and classifications for a large set of proteins . one tab is given for the analysis using the entire protein skin and additional tabs are used for each spherical region specified for analysis by the user . a further tab has a java applet with the astex viewer ( 22 ) to allow the user to visually check the superposition of the protein structures as well as the positioning of the spherical regions used for focused comparisons of the electrostatic potentials . as an example , the dataset studied in reference ( 9 ) is provided on the web server for download . these structures of triosephosphate isomerase ( tpi ) from 12 different species are already superimposed ( five of them are shown in figure 2 ) . the example shows how their electrostatic potentials in the vicinity of the active site can be analysed and related to enzymatic kinetic parameters . for example , they have been shown to affect the rate of ligand binding , the affinity of ligand binding and the stability of the transition state . whether a relation between electrostatic potential differences and kinetic parameters or binding affinities exists in a particular case under study depends on the relative contribution of electrostatic interactions to these quantities and the consistency of the structural and the kinetic or thermodynamic experimental data . after the upload of the protein coordinates using an applet , one can request comparison of the electrostatic potentials in specific regions . for the triosephosphate isomerase example , regions may be selected as follows ( with the cartesian coordinates of the center and the radius of a sphere in  given ) : 
 the substrate binding site ( region_km : 1.07 , 4.06 , 21.11 , 15 ) where the electrostatic potential correlates with substrate km values , andthe active site ( region_kcat_km : 8.15 , 3.54 , 34.83 , 15 ) where the catalytic turnover occurs and the electrostatic potential influences enzymatic kcat / km values ( 9 ) . 
 the electrostatic potentials are computed for the entire proteins using uhbd . for restricting the region of comparison in case ( i ) , a spherical region of radius 15  around w168 this residue is part of the flexible loop region of the tpis which closes over the active site when the substrate binds . a comparison of the similarity of the electrostatic potentials ( measured by the hodgkin si ) in this region for the five model proteins is shown on the bottom right corner of figure 2 . the five enzymes form two subclusters : tpis from chicken , human and rabbit form one subcluster , whereas yeast and escherichia coli form a distinct second subcluster . this is in agreement with the different substrate km values for the two subclusters : tpis from chicken , human and rabbit possess lower km values ( 0.47 mm , 0.49 mm and 0.43 mm , respectively ) and yeast and e. coli exhibit significantly higher km values ( 1.22 mm and 1.03 mm , respectively ) . for a thorough discussion , thus , the electrostatic potential differences in the tpi loop region can be used as a fingerprint for classifying tpis from various species . such a pipsa analysis can provide insights that are complementary to and not apparent from a protein sequence - based analysis . the substrate binding site ( region_km : 1.07 , 4.06 , 21.11 , 15 ) where the electrostatic potential correlates with substrate km values , and the active site ( region_kcat_km : 8.15 , 3.54 , 34.83 , 15 ) where the catalytic turnover occurs and the electrostatic potential influences enzymatic kcat / km values ( 9 ) . webpipsa is implemented using java servlets and java server pages ( jsp ) running on tomcat 5 ( http://tomcat.apache.org ) . the workflows with significant computational components are implemented as ant scripts ( http://ant.apache.org ) and are launched from a java messaging server ( jms , http://activemq.apache.org ) . this architecture allows the separation of the tomcat web server from the computationally demanding workflows . the user is informed about the current state of the workflow by messages sent from the ant script via the messaging server . all data are stored on the file system and may be removed after 2 weeks without further notice . currently , webpipsa provides a description and categorization of the electrostatic potential differences between the input protein structures . it does not include all the features of the downloadable pipsa software which can be used to analyze other types of molecular interaction field and to select conical as well as spherical regions . these features are planned to be added in future developments of webpipsa . in addition , it is planned to add tools to webpipsa to enable the user to study the relations between protein molecular interaction fields and enzymatic kinetic parameters , as in qpipsa ( 9 ) . on the other hand , webpipsa is much more user - friendly and , therefore , accessible to non - experts . it also has additional analysis features such as the simultaneous display of colored heat maps and epograms as well as protein structure visualization with the astex viewer .
protein molecular interaction fields are key determinants of protein functionality . 
 pipsa ( protein interaction property similarity analysis ) is a procedure to compare and analyze protein molecular interaction fields , such as the electrostatic potential . 
 pipsa may assist in protein functional assignment , classification of proteins , the comparison of binding properties and the estimation of enzyme kinetic parameters . 
 webpipsa is a web server that enables the use of pipsa to compare and analyze protein electrostatic potentials . 
 while pipsa can be run with downloadable software ( see http://projects.eml.org/mcm/software/pipsa ) , webpipsa extends and simplifies a pipsa run . 
 this allows non - expert users to perform pipsa for their protein datasets . with input protein 
 coordinates , the superposition of protein structures , as well as the computation and analysis of electrostatic potentials , is automated . 
 the results are provided as electrostatic similarity matrices from an all - pairwise comparison of the proteins which can be subjected to clustering and visualized as epograms ( tree - like diagrams showing electrostatic potential differences ) or heat maps . 
 webpipsa is freely available at : http://pipsa.eml.org .