We address the problem of gradient estimation with respect to four characterizing parameters of the Meixner distribution and Lévy process. With the help of the explicit marginal probability density function, the likelihood ratio method is directly applicable, while unbiased estimators may contain infinite random series in their score function. We quantify the estimator bias arising when the infinite series is truncated to finite term. We further propose a substantially simple exact simulation method for the Meixner distribution, based on acceptance-rejection sampling and the Esscher density transform. Numerical results are presented in the context of financial Greeks to illustrate the effectiveness of our formulas along with bias estimates.