A new algorithm meant for biomedical image retrieval application is presented in this paper. The local region of image is represented by peak valley edge patterns (PVEP), which are calculated by the first-order derivatives in 0º, 45º, 90º and 135º directions. The PVEP differs from the existing local binary pattern (LBP) in a manner that it extracts the directional edge information based on first-order derivative in an image. Further, the effectiveness of our algorithm is confirmed by combining it with Gabor transform. The performance of the proposed method is tested on VIA/I-ELCAP database which includes region of interest computer tomography (ROI-CT) images. Performance analysis shows that the proposed method improves retrieval results from 79.21% to 86.13% and 51.91% to 55.06% as compared to LBP in terms of average precision when number of top matches considered is 10 and 100 respectively.