Mathematical modelling and computational study of two-dimensional and three-dimensional dynamics of receptor-ligand...

Cell signalling processes involve receptor trafficking through highly connected networks of interacting components. The binding of surface receptors to their specific ligands is a key factor for the control and triggering of signalling pathways. But the binding process still presents many enigmas and, by analogy with surface catalytic reactions, two different mechanisms can be conceived: the first mechanism is related to the Eley-Rideal (ER) mechanism, i.e. the bulk-dissolved ligand interacts directly by pure three-dimensional (3D) diffusion with the specific surface receptor; the second mechanism is similar to the Langmuir-Hinshelwood (LH) process, i.e. 3D diffusion of the ligand to the cell surface followed by reversible ligand adsorption and subsequent two-dimensional (2D) surface diffusion to the receptor. A situation where both mechanisms simultaneously contribute to the signalling process could also occur. The aim of this paper is to perform a computational study of the behavior of the signalling response when these different mechanisms for ligand-receptor interactions are integrated into a model for signal transduction and ligand transport. To this end, partial differential equations have been used to develop spatio-temporal models that show trafficking dynamics of ligands, cell surface components, and intracellular signalling molecules through the different domains of the system. The mathematical modeling developed for these mechanisms has been applied to the study of two situations frequently found in cell systems: (a) dependence of the signal response on cell density; and (b) enhancement of the signalling response in a synaptic environment.