Matrix Mittag–Leffler distributions and modeling heavy-tailed risks