Email updates

Keep up to date with the latest news and content from BMC Structural Biology and BioMed Central.

Open Access Software

ProFace: a server for the analysis of the physicochemical features of protein-protein interfaces

Rudra P Saha, Ranjit P Bahadur, Arumay Pal, Saptarshi Mandal and Pinak Chakrabarti*

Author Affiliations

Department of Biochemistry, Bose Institute, P-1/12 CIT Scheme VIIM, Calcutta 700 054, India

For all author emails, please log on.

BMC Structural Biology 2006, 6:11  doi:10.1186/1472-6807-6-11


The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1472-6807/6/11


Received:20 February 2006
Accepted:7 June 2006
Published:7 June 2006

© 2006 Saha et al; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

Molecular recognition is all pervasive in biology. Protein molecules are involved in enzyme regulation, immune response, signal transduction, oligomer assembly, etc. Delineation of physical and chemical features of the interface formed by protein-protein association would allow us to better understand protein interaction networks on one hand, and to design molecules that can engage a given interface and thereby control protein function on the other hand.

Results

ProFace is a suite of programs that uses a file, containing atomic coordinates of a multi-chain molecule, as input and analyzes the interface between any two or more subunits. The interface residues are shown segregated into spatial patches (if such a clustering is possible based on an input threshold distance) and/or core and rim regions. A number of physicochemical parameters defining the interface is tabulated. Among the different output files, one contains the list of interacting residues across the interface. Results can be used to infer if a particular interface belongs to a homodimeric molecule.

Conclusion

A web-server, ProFace (available at http://www.boseinst.ernet.in/resources/bioinfo/stag.html webcite) has been developed for dissecting protein-protein interfaces and deriving various physicochemical parameters.

Background

Most proteins function by interacting with other molecules; the binding sites have evolved for achieving specific interactions and avoiding undesirable associations that would be deleterious to the normal functioning of the cell. Thus the interfaces between two protein subunits provide context for understanding the principles of molecular recognition. A large volume of structural data on protein interactions, either complexes between independent polypeptide chains, or oligomeric assembly of subunits, is available in the Protein Data Bank (PDB) [1], which has been used to generate diverse datasets of protein-protein interfaces [2]. The physical and chemical features of the interfaces have been analyzed [3-8] and softwares/websites, such as Protein-Protein Interaction Server [6], MolSurfer [9], SPIN-PP [10], etc. are available for their calculations. Nevertheless, our understanding of the biomolecular interactions is not adequate enough, for example, to infer unambiguously the arrangement of the subunits in an oligomeric protein from crystallographic studies [11], or to ascertain a high success rate for the prediction of models of protein-protein complexes through docking methods [12].

Recently, protein-protein interfaces have been dissected from new perspectives [13,14]. It has been shown that many large interfaces are not contiguous, but built of spatially demarcated surface patches. Such segregation into patches is also indicative of the location and distribution of water molecules held in the interface [15]. Additionally, one can also divide the interface into core and rim regions using the difference of solvent accessibilities of residues and the chemical properties of each region are quite distinct. Interestingly, this division also mirrors the degree of conservation of interface residues in a family of homologous proteins [16], and this represents an important signature of protein interaction sites. Various other physicochemical parameters have also been developed [17,18], which in combination, can distinguish the true oligomeric state (dimer, in particular) from the lattice contacts observed in protein crystals. In this article we describe a web-server, ProFace that dissects a given

    pro
tein-protein inter
    face
and obtains various parameters to characterize it.

Implementation and results

Input file and parameters

All the protein chains should be contained in the input file in the PDB format and the user must indicate which chains (a maximum of three allowed) constitute each of the two components forming the interface between them. Also, one has to specify the way to display the dissected interface, i.e., to show the residues belonging to core and rim and/or in spatial patches. For clustering into patches the threshold distance has to be supplied. This distance should typically be half the maximum distance between any two interface atoms on a given protein chain – the latter distance is listed along with the other parameters in the output. Ideally, the number of patches should be the same on both the components and if this is not the case the threshold value may have to be slightly changed (increase to reduce the number of patches and vice-versa) to achieve this. The suggested values are 15 Å for protein-protein complexes [13] and 22 Å for homodimers [14], as these gave patches that were visually meaningful in the vast majority of the cases.

Output files and parameters

There are five types of output: a) plot of interface residues with secondary structural information; b) statistics of interface parameters; c) coordinates of interface atoms and the PDB files in which the interface residues are tagged; d) list of residue contacts across interface; and e) the view of the interface atoms.

Plot of interface residues with secondary structural information

The secondary structural elements (α-helix and β-strand) are computed using the program DSSP [19] and shown below the residue names (one-letter code) along the sequence for the individual chains. The sequence information is based on residues for which coordinates are available (and not on the basis of SEQRES records). There are three options to show the interface residues: (i) to simply show the interface residues (in red color); (ii) to show them dissected into core/rim regions (red/blue color); and to show them dissected in two different ways – spatial patches (in different colors) and core/rim regions (upper/lower case). An example of option (iii) is displayed in Figure 1.

thumbnailFigure 1. The interface residues shown against the sequence of c-AMP dependent protein kinase in complex with H7 protein kinase inhibitor 1-(5-isoquinolinesulfonyl)-2-methylpiperazine (PDB file, 1ydr) [24]. There are two patches and the residues belonging to them are shown in orange and magenta (in decreasing patch size). Core and rim residues are distinguished by upper and lower-case letters, respectively. An α-helix is represented by red undulation and a β-strand by blue arrow.

Statistics of interface parameters

A typical example of output parameters is shown in Table 1. The interface area is the sum of the solvent accessible surface areas (ASA) of the two components less that of the pair. ASA is calculated using program NACCESS [20]. All protein atoms or residues contributing more than 0.1 Å2 to the interface area are counted as interface atoms or residues, whose numbers are tabulated. Non-polar interface area is the area contributed by non-polar interface atoms (i.e., all atoms excluding O, N and S). Interface area/surface area is the ratio of the interface area to the rest area of the protein surface in the two components. Fraction of non-polar atoms is based not on the area contributed, but on the number of atoms. Fraction of fully buried atoms is the ratio of interface atoms that are completely buried in the complex (with ASA = 0) to the total number of interface atoms (which also include atoms that do not have zero ASA in the complex). Residue propensity score and local density are defined in Bahadur et al. [17]. Residues with at least one fully buried interface atom are designated as core residues, while rim residues do not contain any interface atom that is fully-buried. Once a residue is identified as core, all its constituent atoms are assumed to be in core (irrespective of the atom being fully or partially buried) and the interface area contributed by the atoms of the residue is part of the core region. Statistics also include atoms/residues/areas divided into core and rim regions (Table 2). Also the number of patches in individual chains and their respective sizes are tabulated (Table 3).

Table 1. Interface parameters of c-AMP-dependent protein kinase complex (PDB code, 1ydr) [24]

Table 2. Statistics on the core and rim regions of the interface in the file, 1ydr

Table 3. Areas of individual patches in the interface of the two components in 1ydr

Output files

The 4-digit code used to name the output files are randomly generated and does not have any correspondence to the input file name. The coordinates are stored in two types of files (with extensions .pdb and .int) and there are two files (corresponding to individual components) of each type. In the .pdb file the interface residues are distinguished from the remaining atoms in the structure on the basis of the content in the two columns – occupancy factor and B-factor. The non-interface residues have a value of 0.00 in these columns. For the interface residues, a) the occupancy is replaced by -5.00 (if it is a core residue) or 5.00 (if it is a rim residue); b) the B-factor column is replaced by a value 1.00 through 9.00, depending on the patch to which the residue belongs. In the .int file, only the interface atoms are kept, with the occupancy and the B-factor column modified as above (and an additional information on patches is also provided by appending labels a, b, c, ... to the keyword ATOM to correspond to patch numbers 1, 2, 3,...). Moreover, there are two additional columns, in which the ASAs of the constituent atoms in the individual component and in the complex are provided. One can use this information to calculate the interface area contributed by individual residues and, for example, correlate with the thermodynamic data on the free energy of binding [16].

Another output file (with extension .cont) provides the list of residue contacts across the interface. For an interface residue in the first component the list shows the interface residues from the other component which are within a distance of 4.5 Å. If a pair of residues in contact have the same residue name and number, this is indicated by the symbol '<< ---' at the end of the line. This interaction has been designated as self-contact and indicates that the interface may have been formed by components/chains related by a 2-fold symmetry [18]. An example of the presence of self-contacting residues in a homodimer structure is presented in Figure 2. Some of the parameters in Table 1, along with the information on self-contacting residues may be used to ascertain if a 2-fold related contact observed in a crystal structure truly represents a biological homodimer.

thumbnailFigure 2. Self-contacting residues in the dimeric structure of wheat germ agglutinin (9wga) [25]. Residues in the two subunits are in two different colors (and those of one chain labeled), with the 2-fold axis running vertically.

View of the interface atoms

This can be done using either RasMol [21] or CHIME [22], depending on whichever program has been configured by the user on the machine. Clicking on the RasMol link will first enable the user to download the PDB file (with interface atoms), which can then be viewed by either program. Clicking on the CHIME link loads the PDB file directly in CHIME. As the B-factor column of the PDB file has been replaced by number codes indicating the patch to which the atoms belong, the interface atoms can be colored on the basis of patches using RasMol. Also, the PDB file generated by the program can be used in GRASP [23] to color the molecular surface according to the criterion of patch or core/rim region.

Conclusion

ProFace can be used to dissect a protein-protein interface, deriving physicochemical parameters. The output can be used to display the interface with standard softwares and understand the biological significance of the interaction.

Availability and requirements

• Project name: ProFace

• Project home page: http://www.boseinst.ernet.in/resources/bioinfo/stag.html webcite

• Operating system(s): Platform independent

• Programming language: Java, C++

• Other requirements: JRE 1.4.2.04 or higher, Chime plug-in 2.6 or higher; all of them are available for download at the above web address

• License: Free

• Any restrictions to use by non-academics: None

Authors' contributions

RPS and RPB wrote the source codes, participated in developing the server. AP and SM participated in developing the server. PC conceived the study, and participated in its design, analysis, and coordination. RPS, RPB, AP, SM and PC all contributed to writing the final manuscript and interpretation of data.

Acknowledgements

We are grateful to Prof. J. Janin and Dr. F. Rodier, with whose collaboration many of the parameters discussed here were developed. The computational facility has been provided by the Department of Biotechnology, India, and the Council of Scientific and Industrial Research provided fellowships to RPS, RPB, AP and SM.

References

  1. Berman HM, Westbrook J, Feng Z, Gilliland G, Bhat TN, Weissig H, Shindyalov IN, Bourne PE: The Protein Data Bank. [http://www.rcsb.org/] webcite

    Nucleic Acids Res 2000, 28:235-242. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  2. Keskin O, Tsai CJ, Wolfson H, Nussinov R: A new, structurally nonredundant, diverse data set of protein-protein interfaces and its implications.

    Protein Sci 2004, 13:1043-1055. PubMed Abstract | Publisher Full Text OpenURL

  3. Chothia C, Janin J: Principles of protein-protein recognition.

    Nature 1975, 256:705-708. PubMed Abstract | Publisher Full Text OpenURL

  4. Argos P: An investigation of protein subunits and domain interfaces.

    Protein Eng 1988, 2:101-113. PubMed Abstract OpenURL

  5. Lawrence MC, Colman PM: Shape complementarity at protein/protein interfaces.

    J Mol Biol 1993, 234:946-50. PubMed Abstract | Publisher Full Text OpenURL

  6. Jones S, Thornton JM: Principles of protein-protein interactions. [http://www.biochem.ucl.ac.uk/bsm/PP/server/] webcite

    Proc Natl Acad Sci USA 1996, 93:13-20. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  7. Tsai CJ, Nussinov R: Hydrophobic folding units at protein-protein interfaces: implications to protein folding and to protein-protein association.

    Protein Sci 1997, 6:1426-1437. PubMed Abstract | Publisher Full Text OpenURL

  8. Lo Conte L, Chothia C, Janin J: The atomic structure of protein-protein recognition sites.

    J Mol Biol 1999, 285:2177-2198. PubMed Abstract | Publisher Full Text OpenURL

  9. Gabdoulline RR, Wade RC, Walther D: MolSurfer: a macromolecular interface navigator. [http://projects.villa-bosch.de/mcm/software/molsurfer] webcite

    Nucleic Acids Res 2003, 31:3349-3351. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  10. SPIN-PP [http://honiglab.cpmc.columbia.edu/SPIN/main.html] webcite

  11. Ponstingl H, Kabir T, Thornton JM: Automatic inference of protein quaternary structure from crystals.

    J Appl Crystallogr 2003, 36:1116-1122. Publisher Full Text OpenURL

  12. Janin J: Assessing predictions of protein-protein interaction: the CAPRI experiment.

    Protein Sci 2005, 14:278-283. PubMed Abstract | Publisher Full Text OpenURL

  13. Chakrabarti P, Janin J: Dissecting protein-protein recognition sites.

    Proteins 2002, 47:334-343. PubMed Abstract | Publisher Full Text OpenURL

  14. Bahadur RP, Chakrabarti P, Rodier F, Janin J: Dissecting subunit interfaces in homodimeric proteins.

    Proteins 2003, 53:708-719. PubMed Abstract | Publisher Full Text OpenURL

  15. Rodier F, Bahadur RP, Chakrabarti P, Janin J: Hydration of protein-protein interfaces.

    Proteins 2005, 60:36-45. PubMed Abstract | Publisher Full Text OpenURL

  16. Guharoy M, Chakrabarti P: Conservation and relative importance of residues across protein-protein interfaces.

    Proc Natl Acad Sci USA 2005, 102:15447-15452. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  17. Bahadur RP, Chakrabarti P, Rodier F, Janin J: A dissection of specific and non-specific protein-protein interfaces.

    J Mol Biol 2004, 336:943-955. PubMed Abstract | Publisher Full Text OpenURL

  18. Saha RP, Bahadur RP, Chakrabarti P: Interresidue contacts in proteins and protein-protein interfaces and their use in characterizing the interface.

    J Proteome Res 2005, 4:1600-1609. PubMed Abstract | Publisher Full Text OpenURL

  19. Kabsch W, Sander C: Dictionary of protein secondary structure. Pattern recognition of hydrogen-bonded and geometrical features. [http://www.cmbi.kun.nl/gv/dssp] webcite

    Biopolymers 1983, 22:2577-2637. PubMed Abstract | Publisher Full Text OpenURL

  20. Hubbard SJ: NACCESS: A program for calculating accessibilities. [http://wolf.bms.umist.ac.uk/naccess/] webcite

    Department of Biochemistry and Molecular Biology University College of London; 1992. OpenURL

  21. Sayle RA, Milner-White EJ: RasMol: biomolecular graphics for all. [http://www.bernstein-plus-sons.com/software/rasmol/] webcite

    Trends Biochem Sci 1995, 20:374-375. PubMed Abstract | Publisher Full Text OpenURL

  22. CHIME [http://www.mdlchime.com] webcite

  23. Nicholls A, Sharp K, Honig B: Protein folding and association: insights from the interfacial and thermodynamic properties of hydrocarbons. [http://trantor.bioc.columbia.edu/grasp/] webcite

    Proteins 1991, 11:281-296. PubMed Abstract | Publisher Full Text OpenURL

  24. Engh RA, Girod A, Kinzel V, Huber R, Bossemeyer D: Crystal structures of catalytic subunit of cAMP-dependent protein kinase in complex with isoquinolinesulfonyl protein kinase inhibitors H7, H8, and H89. Structural implications for selectivity.

    J Biol Chem 1996, 271:26157-26164. PubMed Abstract | Publisher Full Text OpenURL

  25. Wright CS: 2.2 Å resolution structure analysis of two refined N-acetylneuraminyl-lactose-wheat germ agglutinin isolectin complexes.

    J Mol Biol 1990, 215:635-651. PubMed Abstract OpenURL