Molecularly-imprinted polymers (MIPs) are a sample preparation technique with a functional cavity that can rebind the same or similar analyte to the template. It also can increase the selectivity of the analytical method. Computational approaches are capable of making predictions about the most optimal synthesis conditions, including the most appropriate monomers and solvents. This study used the DFT method, GFN2-xTB, and molecular dynamics simulation to predict MIP synthesis’s best functional monomer, solvent, and stoichiometric ratio. Remdesivir (RMD), used as a template, has also never been developed to be made as MIP. The monomers were acrylic acid (AA), 2-hydroxyethyl methacrylate (HEMA), and acrylamide (AM). Laboratory testing, association constant, and Job plot analysis are carried out to confirm the results of computational tests. The polymer was synthesized using precipitation and characterized using SEM, FTIR, and TGA. The adsorption capacity was evaluated in two different solvents, i.e. acetonitrile (ACN) and water: ACN=9: 1 (v/v). Computational study results identified acrylamide as the most interactive ligand with RMD as a template. The results of molecular dynamics simulations identify that the RMD: AM=1: 1 complex is the most stable in terms of the radial distribution function, hydrogen bond occupancy, and Gibbs bond free energy. The adsorption ability test results show that MIP 1 and MIP 2 have imprinting factor values of 1.36 and 1.15, respectively.