In this paper we study a methodology for the numerical simulation of stable structures of fluid membranes and vesicles in biological organisms. In particular, we discuss the effects of spontaneous curvature on vesicle cell membranes under the bending energy for given volume and surface area. The geometric modeling of the vesicle shapes are undertaken by means of surfaces generated as Partial Differential Equations (PDEs). We combine PDE based geometric modeling with numerical optimization in order to study the stable shapes adopted by the vesicle membranes. Thus, through the PDE method we generate a generic template of a vesicle membrane which is then efficiently parameterized. The parameterization is taken as a basis to set up a numerical optimization procedure which enables us to predict a series of vesicle shapes subject to given surface area and volume.