Email updates

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

This article is part of the supplement: Highlights from the Third International Society for Computational Biology (ISCB) Student Council Symposium at the Fifteenth Annual International Conference on Intelligent Systems for Molecular Biology (ISMB)

Open Access Oral presentation

Variation of geometrical and physicochemical properties in protein binding pockets and their ligands

Abdullah Kahraman1*, Richard J Morris2, Roman A Laskowski1 and Janet M Thornton1

Author Affiliations

1 European Bioinformatics Institute, Thornton Group, Wellcome Trust Genome Campus, Cambridge, CB10 1SD, UK

2 John Innes Centre, Computational Biology Group, Norwich Research Park, Colney Lane, Norwich, NR4 7UH, UK

For all author emails, please log on.

BMC Bioinformatics 2007, 8(Suppl 8):S1  doi:10.1186/1471-2105-8-S8-S1

The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-2105/8/S8/S1


Published:20 November 2007

© 2007 Kahraman et al; licensee BioMed Central Ltd.

Background

Physicochemical complementarity is commonly believed to be the driving force for molecular binding. The complementarity for example of electrostatic potentials is regarded as the force that draws the ligand from the solvent into the binding site [1]. If this hypothesis is true, the same ligand should encounter complementarity environmental properties in all proteins to which it binds. We have used our recently published ligand and binding pocket matching algorithm [2] to test this common assumption by searching for property distributions that are similar for the same ligand bound to different proteins.

Methods

The algorithm bases on real spherical harmonic functions, which are applicable to approximate any property function on a unit sphere. These property functions can either be of geometrical or physicochemical nature. For our current analysis we used the shape of binding pockets to test their geometrical similarity and mapped electrostatic, van der Waals and hydrophobicity potentials of the protein on the ligand surface to simulate the physicochemical forces that a ligand may feel in its binding site.

Results

It was discovered that, of these properties the two that vary least for a given ligand are the binding conformation of the ligand followed by the shape of the binding pocket. Conversely, the same ligand encountered very different electrostatic and van der Waals potential environments in the different proteins to which it is bound. These properties were often found not to be complementary to the ligand's properties, which is in conflict with the general assumption stated above. However, the hydrophobicity of the binding pocket did seem to correlate with the properties of the ligand bound to the protein. Hydrophobic parts of the ligand are often confronted with hydrophobic parts of the protein, giving rise to similar hydrophobicity distributions within different binding pockets binding the same ligand (see Figure 1).

thumbnailFigure 1. A set of Adenosine-mono-phosphate (AMP) ligands bound to non-homologous binding sites is shown. Each row displays different geometrical and physicochemical properties of the binding site and ligand respectively. From top to bottom are shown the variation of the ligand shape, the binding pocket shape, the hydrophobicity of the protein mapped on the ligand shape, the van der Waals potential and the electrostatic potential both again mapped on the ligand shape. The properties were ordered according to their average degree of similarity among the different binding pockets from highest to lowest from top to bottom. In addition the AMP binding pockets were ranked according to the similarity of their bound ligand to the AMP ligand of the Protein Data Bank [3] structure 1amu [4].

Conclusion

These results demonstrate that binding sites that bind the same ligand can exhibit a large variation of properties by facing different physicochemical forces within different binding sites. The results urge a re-evaluation of the total contribution of some physicochemical properties to molecular recognition and the factors that drive molecular binding.

Acknowledgements

The molecules in the figure were rendered using PyMOL (W.L. DeLano, http://pymol.sourceforge.net/).

References

  1. Livesay DR, Jambeck P, Rojnuckarin A, Subramaniam S: Conservation of electrostatic properties within enzyme families and superfamilies.

    Biochemistry 2003, 42(12):3464-3473. PubMed Abstract | Publisher Full Text OpenURL

  2. Kahraman A, Morris RJ, Laskowski RA, Thornton JM: Shape variation in protein binding pockets and their ligands.

    J Mol Biol 2007, 368:283-301. PubMed Abstract | Publisher Full Text OpenURL

  3. Berman HM, Westbrook J, Feng Z, Gilliland G, Bhat TN, Weissig H, Shindyalov IN, Bourne PE: The Protein Data Bank.

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

  4. Conti E, Stachelhaus T, Marahiel MA, Brick P: Structural basis for the activation of phenylalanine in the non-ribosomal biosynthesis of gramicidin S.

    Embo J 1997, 16:4174-4183. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL