The aggregation or clustering of proteins and other macromolecules plays an important role in the formation of large-scale molecular assemblies within cell membranes. Examples of such assemblies include lipid rafts, and postsynaptic domains (PSDs) at excitatory and inhibitory synapses in neurons. PSDs are rich in scaffolding proteins that can transiently trap transmembrane neurotransmitter receptors, thus localizing them at specific spatial positions. Hence, PSDs play a key role in determining the strength of synaptic connections and their regulation during learning and memory. Recently, a two-dimensional (2D) diffusion-mediated aggregation model of PSD formation has been developed in which the spatial locations of the clusters are determined by a set of fixed anchoring sites. The system is kept out of equilibrium by the recycling of particles between the cell membrane and interior. This results in a stationary distribution consisting of multiple clusters, whose average size can be determined using an effective mean-field description of the particle concentration around each anchored cluster. In this paper, we derive corrections to the mean-field approximation by applying the theory of diffusion in singularly perturbed domains. The latter is a powerful analytical method for solving two-dimensional (2D) and three-dimensional (3D) diffusion problems in domains where small holes or perforations have been removed from the interior. Applications range from modeling intracellular diffusion, where interior holes could represent subcellular structures such as organelles or biological condensates, to tracking the spread of chemical pollutants or heat from localized sources. In this paper, we take the bounded domain to be the cell membrane and the holes to represent anchored clusters. The analysis proceeds by partitioning the membrane into a set of inner regions around each cluster, and an outer region where mean-field interactions occur. Asymptotically matching the inner and outer stationary solutions generates an asymptotic expansion of the particle concentration, which includes higher-order corrections to mean-field theory that depend on the positions of the clusters and the boundary of the domain. Motivated by a recent study of light-activated protein oligomerization in cells, we also develop the analogous theory for cluster formation in a three-dimensional (3D) domain. The details of the asymptotic analysis differ from the 2D case due to the contrasting singularity structure of 2D and 3D Green’s functions.Graphical abstract2D model of diffusion-based protein cluster formation in the presence of anchoring cites and particle recycling. a A set of N anchoring sites at positions xj\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$${\ extbf{x}}_j$$\\end{document}, j=1,…,N\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$j=1,\\ldots ,N$$\\end{document}, in a bounded domain Ω\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\Omega $$\\end{document}. b Diffusing particles accumulate at the anchoring sites resulting in the formation of particle aggregates or clusters Uj\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$${{\\mathcal {U}}}_j$$\\end{document}. c The clusters are dynamically maintained by a combination of lateral diffusion outside the clusters and particle recycling