Cerebroside sulfotransferase (CST) is a key enzyme in sulfatide biosynthesis and regulation of the myelin sheath in the nervous system. To counter sulfatide accumulation with the deficiency of aryl sulfatase A, CST is considered a target protein in substrate reduction therapy in metachromatic leukodystrophy. In this study, 461 phytoconstituents from four herbs of Medhya Rasayana were screened using multi-pronged virtual screening methods including molecular docking, molecular dynamics (MD) simulation, and reverse pharmacophore analysis. The initial screening of the top 15 hits was based on the binding affinity of the compounds toward the CST substrate-binding site using the lowest free energy of a binding score cutoff of ≤ -7.5 kcal/mol, with the number of conformations in the largest cluster more than 75. The absorption, distribution, metabolism, and excretion (ADME) and toxicity-based pharmacokinetic analysis delivered the top four hits: 18alpha-glycyrrhetinic acid, lupeol, alpha carotene, and beta-carotene, with high blood-brain barrier permeability and negligible toxicity. Furthermore, a 100-ns simulation of protein-ligand complexes with a trajectory analysis of structural deviation, compactness, intramolecular interactions, principal component analysis, free energy landscape, and dynamic cross-correlation analysis showed the binding potential and positioning of the four hits in the binding pocket. Thus, an in-depth analysis of protein-ligand interactions from pre- and post-molecular dynamics simulation, along with reverse pharmacophore mapping, suggests that 18alpha-glycyrrhetinic acid is the most potent and specific CST inhibitor, while beta-carotene could be considered the second most potent compound for CST inhibition as it also exhibited overall stability throughout the simulation. Therefore, the computational drug screening approach applied in this study may contribute to the development of oral drugs as a therapeutic option for metachromatic leukodystrophy.