Volume 1 Supplement 1

BioSysBio 2007: Systems Biology, Bioinformatics, Synthetic Biology

Open Access

Mathematical modelling of multi-site protein-ligand interactions

  • Najl V Valeyev1, 2Email author,
  • Nikolai V Kotov2,
  • Ian Postlethwaite1 and
  • Declan Bates1
BMC Systems Biology20071(Suppl 1):P65

DOI: 10.1186/1752-0509-1-S1-P65

Published: 8 May 2007


Protein molecules with multiple binding sites play a fundamental role in regulation of a large number of ferments, ion channels and other cellular components. Mathematical modelling of multisite protein-ligand interactions which should be based on realistic assumptions is very important for understanding complex biological networks.

Materials and methods

Frequently the experimental dose-response curves for multisite protein interactions with ligand are analysed using the Hill equation [1]. Although the Hill equation is widely used in biological pathways modelling as well as in pharmacodynamic studies [2], it is derived based on a physically unrealistic assumption according to which all ligand molecules simultaneously bind to a protein. In effect, this assumption is only valid in the case of extreme positive cooperativity.


We developed models for physically realistic consecutive and independent binding schemes and compared them with the Hill equation (Figure 1). It follows from this comparison that the "half-effect" concentration may significantly differ from the conventional Hill equation depending on the binding mechanism.
Figure 1

Comparison of multisite protein ligand binding models. The dose-response curves are shown for a protein with four binding sites using following kinetic schemes: a) the Hill equation; b) the consecutive binding; c) the independent binding. The same values have been used for all three models.

We used our models to rigorously investigate mechanisms for second messenger activation by multisite proteins. We show that differential and selective activation of different targets by the same protein-ligand complex (Ca2+ and calmodulin, for example) can be achieved by biologically active non-saturated intermediate multisite protein-ligand complexes.

Authors’ Affiliations

Department of Engineering, Control and Instrumentation Research Group,University of Leicester, University Road, Leicester
Biophysics & Bionics Lab, Department of Physics, Kazan State University


  1. Hill A: The combinations of haemoglobin with oxygen and with carbon monoxide. J Physiology. 1910, 40: 4-7.Google Scholar
  2. Aviles P, Aliouat EM, Martinez A, Dei-Cas E, Herreros E, Dujardin L, Gargallo-Viola D: In vitro harmacodynamic parameters of sordarin derivatives in comparison with those of marketed compounds against Pneumocystis carinii isolated from rats. Antimicrob Agents Chemother. 2000, 44: 1284-90. 10.1128/AAC.44.5.1284-1290.2000PubMedPubMed CentralView ArticleGoogle Scholar


© Valeyev et al; licensee BioMed Central Ltd. 2007

This article is published under license to BioMed Central Ltd.