The wall shear stress (WSS) spatial and temporal gradients are two hemodynamics parameters correlated with endothelial damage. Those two gradients become well pronounced in a bypass graft anastomosis geometry where the blood flow patterns are quite disturbed. The WSS gradient minimisation on the host artery floor can be achieved by optimising the anastomosis shape and hence may lead to an improved long-term post-surgical performance of the graft. The anastomosis shape optimisation can be executed via an integrated computational tool comprised of a meshless computational fluid dynamics (CFD) solver and a genetic algorithm (GA) shape optimiser. The meshless CFD solver serves to evaluate the WSS gradients and the GA optimiser serves to search for the end-to-side distal anastomosis (ETSDA) optimal shape that best minimises those gradients. We utilise a meshless CFD method to resolve hemodynamics and a GA for the purpose of optimisation. We consider three different anastomotic models: the conventional ETSDA, the Miller Cuff ETSDA and the hood ETSDA. The results reported herein demonstrate that the graft calibre should always be maximised whether a conventional or Miller Cuff ETSDA model is utilised. Also, it was noted that the Miller Cuff height should be minimised. The choice of an optimal anastomotic angle should be optimised to achieve a compromise between the concurrent minimisations of both the spatial WSS gradient and the temporal WSS gradient.