ALT Protein–protein interactions (PPIs) regulate signalling pathways and cell phenotypes, and the visualization of spatially resolved dynamics of PPIs would thus shed light on the activation and crosstalk of signalling networks. Here we report a method that leverages a sequential proximity ligation assay for the multiplexed profiling of PPIs with up to 47 proteins involved in multisignalling crosstalk pathways. We applied the method, followed by conventional immunofluorescence, to cell cultures and tissues of non-small-cell lung cancers with a mutated epidermal growth-factor receptor to determine the co-localization of PPIs in subcellular volumes and to reconstruct changes in the subcellular distributions of PPIs in response to perturbations by the tyrosine kinase inhibitor osimertinib. We also show that a graph convolutional network encoding spatially resolved PPIs can accurately predict the cell-treatment status of single cells. Multiplexed proximity ligation assays aided by graph-based