Abstract: A computational approach to simulate the formation of possible imprinted polymers in acetonitrile solution for theophylline (THO) is proposed, using combined molecular dynamics (MD), molecular mechanics (MM), docking and site mapping computational techniques. Methacrylic acid (MAA) and methylmethacrylate (MMA) monomers are used to simulate possible homo and copolymer structures. The model is able predict binding affinity and selectivity when considering THO analogues, such as caffeine, theobromine, xanthine and 3-methylxanthine. Comparison with available experimental data is proposed Template and target information:theophylline, THO Author keywords:molecular imprinting, Molecular dynamics, molecular mechanics