Abstract
Molecular imprinting is a promising approach for developing polymeric materials as artificial receptors. However, only a few types of molecularly imprinted polymers (MIPs) are commercially available, and most research on MIP(S) is still in the experimental phase. The significant limitation has been a challenge for screening imprinting systems, particularly for weak functional target molecules. Herein, a combined method of quantum mechanics (QM) computations and molecular dynamics (MD) simulations was employed to screen an appropriate 2,4-dichlorophenoxyacetic acid (2,4-D) imprinting system. QM calculations were performed using the Gaussian 09 software. MD simulations were conducted using the Gromacs2018.8 software suite. The QM computation results were consistent with those of the MD simulations. In the MD simulations, a realistic model of the 'actual' pre-polymerisation mixture was obtained by introducing numerous components in the simulations to thoroughly investigate all non-covalent interactions during imprinting. This study systematically examined MIP systems using computer simulations and established a theoretical prediction model for the affinity and selectivity of MIPs. The combined method of QM computations and MD simulations provides a robust foundation for the rational design of MIPs.