Adhesion and insertion of curvature-mediating proteins can induce dramatic structural changes in cell membranes, allowing them to participate in several key cellular tasks. The way proteins interact to generate curvature remains largely unclear, especially at early stages of membrane remodeling. Using a coarse-grained model of Bin/ amphiphysin/Rvs domain with an N-terminal helix (N-BAR) interacting with flat membranes and vesicles, we demonstrate that at low protein surface densities, binding of N-BAR domain proteins to the membrane is followed by a linear aggregation and the formation of meshes on the surface. In this process, the proteins assemble at the base of emerging membrane buds. Our work shows that beyond a more straightforward scaffolding mechanism at high bound densities, the interplay of anisotropic interactions and the local stress imposed by the N-BAR proteins results in deep invaginations and endocytic vesicular bud-like deformations, an order of magnitude larger than the size of the individual protein. Our results imply that by virtue of this mechanism, cell membranes may achieve rapid local increases in protein concentration.coarse-grained simulation | membrane stress | self-assembly | membrane curvature | Gaussian curvature L ipid membranes protect cells and their organelles and serve as mechanical support for proteins involved in signal transduction and cellular trafficking (1). Owing to the way lipids are assembled into bilayers of mesoscopic length and molecular width, membranes behave as elastic and highly dynamic molecular sheets. Their innate multiscale nature is further evident in the local interplay between proteins and lipids, which results in large-scale membrane remodeling characterized by an impressive array of morphologies (2). These properties allow biological membranes to take part in several important cellular processes, as their restructuring is key to enabling communication between cells, formation of organelles, trafficking, division, and cell migration (1).The interactions between individual proteins and lipid molecules alone are insufficient to remodel flat membrane sheets to an appreciable extent, implying that the experimentally observed membrane remodeling is a result of the concerted action of multiple proteins (3) in ways mechanistically differing from the mode of action of a single protein. Considering the multiple timescales required to capture protein dynamics as the membrane deformations are induced, as well as the difficulty in separating individual events of the remodeling process, this mechanism remains largely unclear. It previously was demonstrated that the binding of curvature-inducing nanoparticles may form tubular and vesicular structures, in which particles interact with one another solely via curvature (3-6). Furthermore, analytical theory predicts that anisotropic inclusions may induce budding on the membrane surface (7). However, the way proteins assemble on the membrane and how their oligomerization couples to both local and long-wavelength transformat...