In this article, we generalize the univariate Macdonald distribution to the matrix case and give its derivation using matrix variate gamma distribution. We study several properties such as cumulative distribution function, marginal distribution of submatrix, triangular factorization, moment generating function, and expected values of several functions of the Macdonald matrix. Some of these results are expressed in terms of special functions of matrix arguments and zonal polynomials.