From purely academic standpoints as well as practical situations requiring the design of lightweight structures, Topology Optimization has become a powerful tool in automotive, aerospace, and other weight-sensitive industries. Topology optimization method automates the process of finding the optimal design of a product by distributing the material appropriately in a given domain to fulfill a user-defined objective. Finite Element analyses (FEA) usually underlie these optimization algorithms. However, discrete models have been shown to provide a much better understanding of some important material's behavior and processes. For instance, microstructure such as polycrystalline grains, cracks, phase boundaries, and twin boundaries in phase-transforming materials affects the macroscopic response of solids in a variety of ways. But, so far, discrete particle methods in the literature have been limited to particles connected with springs which cannot mimic the behavior of solids with any given constitutive relation. In this study, we use a discrete particle model with multi-body interactions [1] that can handle complicated constitutive relations and use it to optimize the topology of a domain using the Simplified Isotropic Material with Penalization (SIMP) algorithm [2]. Here, an objective function and related sensitivity for a typical problem of stiffness maximization under a discrete particle model have been defined. In this study, this method is employed for a simply supported beam and the predictions obtained have been validated with those obtained by the traditional topology optimization method [3]. An advantage of this method is that it can account for geometric nonlinearity without added complexity. The change from rollers to hinges for beam supports is an illustration of this, which produces a distinct optimal structure.