In silico approach has been carried out for the determination of drug candidates from anthocyanin derivative as inhibitors of the COVID-19 main protease. Geometry optimization has performed using the DFT/B3LYP/6-31G(d,p) method as an initial step to prepare candidate ligand. The results of molecular docking showed that candidates C5 and C6 had promising results with a grid score smaller than the ligand reference (X77) with a flexible conformation type. Studies on absorption, distribution, metabolism, excretion, and toxicity of C5 and C6 candidates were conducted to study the physicochemical properties of drug candidates and to show good predictive results as drugs. Molecular dynamic simulation uses the ffSB14 force field for 200 ns to study the interaction between ligand and receptor, the system stability, solvent accessibility, energy interactions, and hydrogen bonds. The results show good interaction stability on the C5 complex compared to the reference ligand which is characterized by the binding free energy value of C5 was -42.77 ± 0.37 kcal/mol and X77 was -42.37 ± 0.41 kcal/mol.