High-precision tests of fundamental symmetries are looking for the parity- ($\mathcal{P}$), time-reversal- ($\mathcal{T}$) violating electric dipole moment of the electron (eEDM) as proof of physics beyond the Standard Model. Particularly, in polyatomic molecules, the complex vibrational and rotational structure gives the possibility to reach high enhancement of the $\mathcal{P},\mathcal{T}$-odd effects in moderate electric fields, and with the possibility of increasing the statistical sensitivity by using laser cooling. In this work, we calculate the $\mathcal{P},\mathcal{T}$-odd molecular enhancement factor of the eEDM (${W}_{\mathrm{d}}$) and of the scalar-pseudoscalar interaction (${W}_{\mathrm{s}}$) necessary for the interpretation of future experiments on the promising candidates ${\mathrm{BaCH}}_{3}$ and ${\mathrm{YbCH}}_{3}$. We employ high-accuracy relativistic coupled cluster methods and systematically evaluate the uncertainties of our computational approach. Compared to other Ba- and Yb-containing molecules, ${\mathrm{BaCH}}_{3}$ and ${\mathrm{YbCH}}_{3}$ exhibit larger ${W}_{\mathrm{d}}$ and ${W}_{\mathrm{s}}$ associated to the increased covalent character of the $M$--C bond. The calculated values are $3.22\ifmmode\pm\else\textpm\fi{}0.12\ifmmode\times\else\texttimes\fi{}{10}^{24}\frac{h\text{Hz}}{e\text{cm}}$ and $13.80\ifmmode\pm\else\textpm\fi{}0.35\ifmmode\times\else\texttimes\fi{}{10}^{24}\frac{h\text{Hz}}{e\text{cm}}$ for ${W}_{\mathrm{d}}$, and $8.42\ifmmode\pm\else\textpm\fi{}0.29\phantom{\rule{4pt}{0ex}}h\mathrm{kHz}$ and $50.16\ifmmode\pm\else\textpm\fi{}1.27\phantom{\rule{4pt}{0ex}}h\mathrm{kHz}$ for ${W}_{\mathrm{s}}$, in ${\mathrm{BaCH}}_{3}$ and ${\mathrm{YbCH}}_{3}$, respectively. The robust, accurate, and cost-effective computational scheme reported in this work makes our results suitable for extracting the relevant fundamental properties from future measurements and also can be used to explore other polyatomic molecules sensitive to various violations of fundamental symmetries.