  Vascular endothelial growth factor receptors (VEGFRs) are receptor tyrosine
kinases (RTKs) that regulate proliferation, migration, angiogenesis and
vascular permeability of endothelial cells. VEGFR1 and VEGFR2 bind vascular
endothelial growth factors (VEGFs), inducing receptor dimerisation and
activation, characterised by phosphorylation of tyrosine residues in their
cytoplasmic domain. Although experimental evidence suggests that RTK signalling
occurs both on the plasma membrane and intra-cellularly, and reveals the role
of endocytosis in RTK signal transduction, we still lack knowledge of VEGFR
phosphorylation-site use and of the spatiotemporal regulation of VEGFR
signalling. In this paper, we introduce four stochastic mathematical models to
study the binding kinetics of vascular endothelial growth factor VEGF-A to
VEGFR1 and VEGFR2, and phosphorylation. The formation of phosphorylated dimers
on the cell surface is a two-step process: diffusive transport and binding. The
first two of our models only consider VEGFR2, which allows us to introduce new
stochastic descriptors making use of a matrix-analytic approach. The two
remaining models describe the competition of VEGFR1 and VEGFR2 for ligand
availability, and are analysed making use of Gillespie simulations and the van
Kampen approximation. Under the hypothesis that bound phosphorylated receptor
dimers are the signalling units, we study the time to reach a threshold number
of such complexes. Our results indicate that the presence of VEGFR1 does not
only affect the timescale to reach a given signalling threshold, but it also
affects the maximum attainable threshold. This result is consistent with the
conjectured role of VEGFR1 as a decoy receptor, that prevents VEGF-A binding to
VEGFR2, and thus, VEGFR2 attaining suitable phosphorylation levels. We identify
an optimum range of ligand concentration for sustained dimer phosphorylation.
