Hyperspectral images (HSIs) are widely used in various tasks such as earth observation and target detection. However, during the imaging process, HSIs are often corrupted by various noises. In this article, we firstly investigate the advantages of traditional physical restoration models and the denoising convolutional neural networks (CNN) for HSIs denoising tasks. The sparse based low-rank representation can explore the global correlations in both the spatial and spectral domains, and the CNN-based denoiser can represent the deep prior which cannot be designed by traditional restoration models. Then, we propose a HSI denoising model with low-rank representation and CNN denoiser prior in the flexible and extensible plug-and-play framework by combining the advantages of the two methods. The proposed model is user-friendly, requiring no retraining. Simulated data experiments show that, compared with competitive methods, the proposed one achieves better denoising results for both additive Gaussian noise and Poissonian noise in various quantitative evaluation indicators. Real data experiments show that the proposed model yields the best performance.