Insights into the dynamics of ligand-induced dimerisation via mathematical modelling and analysis

The vascular endothelial growth factor (VEGF) receptor (VEGFR) system plays a role in cancer and many other diseases. It is widely accepted that VEGFR receptors dimerise in response to VEGF binding. However, analysis of these mechanisms and their implications for drug development still requires further exploration. In this paper, we present a mathematical model representing the binding of VEGF to VEGFR and the subsequent ligand-induced dimerisation. A key factor in this work is the qualitative and quantitative effect of binding cooperativity, which describes the effect that the binding of a ligand to a receptor has on the binding of that ligand to a second receptor, and the dimerisation of these receptors. We analyse the ordinary differential equation system at equilibrium, giving analytical solutions for the total amount of ligand bound. For time-course dynamics, we use numerical methods to explore possible behaviours under various parameter regimes, while perturbation analysis is used to understand the intricacies of these behaviours. Our simulation results show an excellent fit to experimental data, towards validating the model.