A steady-state modeling approach for simulation of antimicrobial peptide-cell membrane interaction.

Antimicrobial Peptides (AMPs) are host defense molecules that initiate microbial death by binding to the membrane. On membrane binding, AMPs undergo changes in conformation and aggregation state to enable killing action. Depending on the AMP and cell membrane characteristics, the nature of binding can be aggregating or non-aggregating, with high/low cooperativity, at single or multiple sites with high/low affinity leading to a unique killing action that needs to be studied individually. In the present study, a steady-state model that simulates AMP-membrane interaction was developed and was used to predict the mechanism of AMP binding. The predictions obtained from the model were validated with experimentally deciphered values available in literature. The model was further used to predict the mechanism for a set of designed AMPs with high sequence similarity to Myeloid Antimicrobial Peptide (MAP) family. Depending on the predicted mechanism, a unique half saturation constant and steepness of response (Hill coefficient) was obtained which was further validated with available data from literature. The model could reliably predict the mechanism, the half saturation constant and the Hill coefficient values. Further based on the analysis, it was observed that aggregation and oligomerization result in drastic killing action in a short range of peptide concentration owing to high Hill coefficient values. Mechanisms such as monomers binding at multiple sites with/without cooperativity result in antimicrobial activity at low half saturation constant though the killing action may not be steep. Thus, the methodology developed here can be used to develop hypothesis for studying AMP-membrane interaction mechanisms.