We here present the first simulation of a complete molecularly imprinted polymer prepolymerization system. Molecular dynamics studies were performed for a system comprising a total of 1199 discrete molecules, replicating the components and concentrations employed in the corresponding polymer synthesis. The observed interactions correlate well with results obtained from (1)H NMR spectroscopic studies. Comparison with simulations performed in the absence of cross-linking agent (ethylene dimethacrylate) demonstrated its significance in the formation of ligand recognition sites. Moreover, the influence of events such as template-template (bupivacaine) and monomer-monomer (methacrylic acid) self-association, porogen-template interactions, and template conformational variability was revealed. The template recognition capacity of the modeled polymer system was verified by synthesis of imprinted and reference polymers and subsequent radioligand binding analysis. Collectively, through a series of statistical analyses of molecular trajectories in conjunction with spectroscopic data it was demonstrated that an ensemble of complex structures is present in the prepolymerization mixture and that this diversity is the basis for the binding site heterogeneity observed in molecularly imprinted polymers (MIPs) prepared using the noncovalent strategy.

No abstract

A combination of theoretical and experimental studies has provided correlations between molecularly imprinted polymer composition, morphology, and recognition behavior obtained using a series of bupivacaine-imprinted methacrylic acid (MAA)–ethylene glycol dimethacrylate copolymers differing in molar ratios of the respective monomers. Results extracted from analysis of molecular dynamics (MD) trajectory data demonstrated that stability and frequency of interactions between bupivacaine and the monomers in the prepolymerization phase were strongly affected by minor changes in polymer composition, which in turn affected binding site affinity and heterogeneity of the imprinted polymers. Moreover, through the characterization of polymer morphology, we show that higher molar fractions of MAA resulted in polymeric materials with increased pore size, a feature that enhanced the binding capacity of the polymers. Furthermore, the results presented point at the strength of MD for predicting MIP-template binding capacity and affinity.

The photophysical behavior of the isomers of the anticoagulant drug warfarin in various solvents and solvent mixtures was investigated using absorption, 1H NMR, and steady-state and time-resolved fluorescence spectroscopies in conjunction with B3LYP-based theoretical treatments. Complex absorption patterns were observed, indicative of the presence of different isomers of warfarin in the various solvents studied. In alkaline aqueous solution, the deprotonated open side form of warfarin is highly dominant and only one S0-->S1 singlet transition could be observed in the absorption spectrum centered at 320 nm. These observations were supported by theoretical density functional calculations (B3LYP) in which the geometries of nine isomers of warfarin were optimized and their respective eight lowest singlet and three lowest triplet excitation energy levels were predicted. Examination of the fluorescence excitation and emission spectra of the isomers in nonpolar and polar organic solvents showed the presence of the deprotonated open side chain form of warfarin in 2-propanol, ethanol, and acetonitrile. Time-resolved fluorescence experiments revealed a short decay time constant, tau1, in all solvents studied while in more polar environments a second longer one, tau2, was evident varying between 0.5 and 1.6 ns depending on solvent polarity. The variation of number and length of fluorescence lifetimes as a function of solvent environment has provided a tool for examining warfarin protein binding. Studies on the binding of warfarin to human serum albumin (HSA) have been undertaken, and different modes of binding were observed which are indicative of binding to the anion-selective Sudlow I and, second, a lower affinity mode of interaction.

The correlation of the recognition properties of a molecularly imprinted polymer (MIP) with the recognition events in pre-polymerisation mixtures is of central importance to our understanding of the molecular imprinting technique. Using the NSAID naproxen as a model template, we have applied parallel theoretical (molecular dynamics) and practical ((1)H-NMR, X-ray crystallography, HPLC, radioligand binding) methods to examine the nature of template-functional monomer complexation. An effective imprint is achieved, despite the presence of only one site on the template which provides for the formation of effective electrostatic interactions with the functional monomer used, 4-vinylpyridine. This is attributed to the creation of a well-defined receptor site for the acidic terminus of the molecule and complementary van der Waals interactions, as described in preliminary simulations of the pre-polymerisation system, and as confirmed for the resultant MIP by HPLC data. Qualitative agreement is also observed between simulation and proton NMR data examining monomer self-association in the presence and absence of the template. On the basis of the data obtained, the role of a cross-linker appears to be more significant for this system than previously anticipated.

A series of molecular dynamics simulations of prepolymerization mixtures for phenylalanine anilide imprinted co-(ethylene glycol dimethacrylate-methacrylic acid) molecularly imprinted polymers have been employed to investigate the mechanistic basis for template selective recognition in these systems. This has provided new insights on the mechanisms underlying template recognition, in particular the significant role played by the crosslinking agent. Importantly, the study supports the occurrence of template self-association events that allows us to resolve debate between the two previously proposed models used to explain this system's underlying recognition mechanisms. Moreover, the complexity of the molecular level events underlying template complexation is highlighted by this study, a factor that should be considered in rational molecularly imprinted polymer design, especially with respect to recognition site heterogeneity.

In principle, molecularly imprinted polymer science and technology provides a means for ready access to nano-structured polymeric materials of predetermined selectivity. The versatility of the technique has brought it to the attention of many working with the development of nanomaterials with biological or biomimetic properties for use as therapeutics or in medical devices. Nonetheless, the further evolution of the field necessitates the development of robust predictive tools capable of handling the complexity of molecular imprinting systems. The rapid growth in computer power and software over the past decade has opened new possibilities for simulating aspects of the complex molecular imprinting process. We present here a survey of the current status of the use of in silico-based approaches to aspects of molecular imprinting. Finally, we highlight areas where ongoing and future efforts should yield information critical to our understanding of the underlying mechanisms sufficient to permit the rational design of molecularly imprinted polymers.

A model for the molecular basis for ligand recognition in bupivacaine imprinted methacrylic acid-ethylene glycol dimethacrylate co-polymers has been developed based upon a series of (1)H-NMR studies in conjunction with HPLC and radioligand binding analyses. (1)H-NMR studies indicated that functional monomer-template complexes survive the polymerisation process, at least up until the stage of gelation. Polymers were synthesised and characterised by surface area analysis (BET), FT-IR and SEM. A combination of zonal and frontal chromatographic studies in aqueous and non-polar media indicate that selectivity arises from a combination of hydrophobic and electrostatic interactions. However, in the concentration regime employed for LC-based studies, ligand recognition in aqueous media was shown to be predominantly non-specific and hydrophobic in character. Radioligand binding studies, in lower ligand binding concentration regimes, permitted closer examination of the higher affinity binding sites. It was shown that the presence of a polar modifier in a non-polar solvent, or an organic modifier in water, produced enhanced selectivity. Variable temperature studies showed that the temperature of binding influences selectivity as well as the apparent number of sites available and that this effect is different in organic and aqueous environments. This indicates that the system studied is more complex in character than is generally appreciated. A comparison of the techniques employed here indicates that although chromatographic studies provide a valuable first-round screen for polymer-ligand selectivities, the level of detail obtainable using radioligand binding studies (lower concentrations and true equilibrium binding) makes them superior for detailed evaluations of molecularly imprinted polymers.