In this study, Density-functional theory/Time-dependent density-functional theory (DFT/TDDFT) and Molecular docking method was used to investigate the effect of methyl acetate, tetrahydrofuran and cyanobenzylidene substituents on the electronic structure and antiviral activity of favipiravir for treating COVID-19. The DFT and TDDFT computations were employed using the Gaussian 09 software package. The values were calculated using the 6-311++G(d, p) basis set and the hybrid B3LYP functional method. Autodock vina software was used for simulations to better predictions and to validate the modified compounds' binding affinities and poses. Results of the study indicate that compounds 1 to 6 all displayed a planar structure, where the pyrazine ring, carboxamide, hydroxyl groups, and other substituents are all situated within the same plane. In addition, the energy gaps (Egap) of these six compounds (Cpd 1, 2, 3, 4, 5, and 6) were compared. The significant dipole moment and binding affinity achieved implies a particular orientation for binding within the target protein, signaling the anticipated strength of the binding interaction. In all six compounds, the electrophilic domain is situated in the vicinity of the amine functional group within the carboxamide compound, whereas the nucleophilic domain encompasses both the carbonyl and hydroxyl groups. The most negatively charged sites are susceptible to electrophilic interactions. In conclusion, compounds 5 and 6 exhibit a high binding affinity of the target protein, while compound 6 has a high energy gap, which could enhance its antiviral activity against the COVID-19 virus.