Molecular dynamics simulations in pre-polymerization mixtures for peptide recognition

Context: Molecularly imprinted polymers (MIPs) have promising applications as synthetic antibodies for protein and peptide recognition. A critical aspect of MIP design is the selection of functional monomers and their adequate proportions to achieve materials with high recognition capacity toward their targets. To contribute to this goal, we calibrated a molecular dynamics protocol to reproduce the experimental trends in peptide recognition of 13 pre-polymerization mixtures reported in the literature for the peptide toxin melittin. Methods: Three simulation conditions were tested for each mixture by changing the box size and the number of monomers and cross-linkers surrounding the template in a solvent-explicit environment. Fully atomistic MD simulations of 350 ns were conducted with the AMBER20 software, with ff19SB parameters for the peptide, gaff2 parameters for the monomers and cross-linkers, and the OPC water model. Template-monomer interaction energies under the LIE approach showed significant differences between high-affinity and low-affinity mixtures. Simulation systems containing 100 monomers plus cross-linkers in a cubic box of 90 ų successfully ranked the mixtures according to their experimental performance. Systems with higher monomer densities resulted in non-specific intermolecular contacts that could not account for the experimental trends in melittin recognition. The mixture with the best recognition capacity showed preferential binding to the 13–26-α-helix, suggesting a relevant role for this segment in melittin imprinting and recognition. Our findings provide insightful information to assist the computational design of molecularly imprinted materials with a validated protocol that can be easily extended to other templates.