Effect of proline position on the antimicrobial mechanism of buforin II

Buforin II (BF2) is a histone-derived antimicrobial peptide that causes cell death by translocating across membranes and interacting with nucleic acids. It contains one proline residue critical for its function. Previous research found that mutations replacing proline lead to decreased membrane translocation and antimicrobial activity as well as increased membrane permeabilization. This study further investigates the role of proline in BF2's antimicrobial mechanism by considering the effect of changing proline position on membrane translocation, membrane permeabilization, and antimicrobial activity. For this purpose, four mutants were made with proline substitution (P11A) or relocation (P11A/G7P, P11A/V12P, P11A/V15P). These mutations altered the amount of helical content. Although antimicrobial activity correlated with the α-helical content for the peptides containing proline, membrane translocation did not. This observation suggests that factors in BF2's bactericidal mechanism other than translocation must be altered by these mutations. To better explain these trends we also measured the nucleic acid binding and membrane permeabilization of the mutant peptides. A comparison of mutant and wild type BF2 activity revealed that BF2 relies principally on membrane translocation and nucleic acid binding for antimicrobial activity, although membrane permeabilization may play a secondary role for some BF2 variants. A better understanding of the role of proline in the BF2 antimicrobial mechanism will contribute to the further design and development of BF2 analogs. Moreover, since proline residues are prevalent among other antimicrobial peptides, this systematic characterization of BF2 provides general insights that can promote our understanding of other systems.     

Modular analysis of hipposin,a histone-derived antimicrobial peptide consisting of membrane translocating and membrane permeabilizing fragments

Antimicrobial peptides continue to garner attention as potential alternatives to conventional antibiotics. Hipposin is a histone-derived antimicrobial peptide (HDAP) previously isolated from Atlantic halibut. Though potent against bacteria, its antibacterial mechanism had not been characterized. The mechanism of this peptide is particularly interesting to consider since the full hipposin sequence contains the sequences of parasin and buforin II (BF2), two other known antimicrobial peptides that act via different antibacterial mechanisms. While parasin kills bacteria by inducing membrane permeabilization, buforin II enters cells without causing significant membrane disruption, harming bacteria through interactions with intracellular nucleic acids. In this study, we used a modular approach to characterize hipposin and determine the role of the parasin and buforin II fragments in the overall hipposin mechanism. Our results show that hipposin kills bacteria by inducing membrane permeabilization, and this membrane permeabilization is promoted by the presence of the N-terminal domain. Portions of hipposin lacking the N-terminal sequence do not cause membrane permeabilization and function more similarly to buforin II. We also determined that the C-terminal portion of hipposin, HipC, is a cell-penetrating peptide that readily enters bacterial cells but has no measurable antimicrobial activity. HipC is the first membrane active histone fragment identified that does not kill bacterial or eukaryotic cells. Together, these results characterize hipposin and provide a useful starting point for considering the activity of chimeric peptides made by combining peptides with different antimicrobial mechanisms. This article is part of a Special Issue entitled: Interfacially Active Peptides and Proteins. Guest Editors: William C. Wimley and Kalina Hristova.     

Effect of lipid composition on buforin II structure and membrane entry

Buforin II is a 21-amino acid polycationic antimicrobial peptide derived from a peptide originally isolated from the stomach tissue of the Asian toad Bufo bufo gargarizans. It is hypothesized to target a wide range of bacteria by translocating into cells without membrane permeabilization and binding to nucleic acids. Previous research found that the structure and membrane interactions of buforin II are related to lipid composition. In this study, we used molecular dynamics (MD) simulations along with lipid vesicle experiments to gain insight into how buforin II interacts differently with phosphatidylcholine (PC), phosphatidylglycerol (PG), and phosphatidylethanolamine (PE) lipids. Fluorescent spectroscopic measurements agreed with the previous assertion that buforin II does not interact with pure PC vesicles. Nonetheless, the reduced entry of the peptide into anionic PG membranes versus neutral PC membranes during simulations correlates with the experimentally observed reduction in BF2 translocation through pure PG membranes. Simulations showing membrane entry into PC also provide insight into how buforin II may initially penetrate cell membranes. Our MD simulations also allowed us to consider how neutral PE lipids affect the peptide differently than PC. In particular, the peptide had a more helical secondary structure in simulations with PE lipids. A change in structure was also apparent in circular dichroism measurements. PE also reduced membrane entry in simulations, which correlates with decreased translocation in the presence of PE observed in previous studies. Together, these results provide molecular-level insight into how lipid composition can affect buforin II structure and function and will be useful in efforts to design peptides with desired antimicrobial and cell-penetrating properties.     

Novel histone-derived antimicrobial peptides use different antimicrobial mechanisms

The increase in multidrug resistant bacteria has sparked an interest in the development of novel antibiotics. Antimicrobial peptides that operate by crossing the cell membrane may also have the potential to deliver drugs to intracellular targets. Buforin 2 (BF2) is an antimicrobial peptide that shares sequence identity with a fragment of histone subunit H2A and whose bactericidal mechanism depends on membrane translocation and DNA binding. Previously, novel histone-derived antimicrobial peptides (HDAPs) were designed based on properties of BF2, and DesHDAP1 and DesHDAP3 showed significant antibacterial activity. In this study, their DNA binding, permeabilization, and translocation abilities were assessed independently and compared to antibacterial activity to determine whether they share a mechanism with BF2. To investigate the importance of proline in determining the peptides' mechanisms of action, proline to alanine mutants of the novel peptides were generated. DesHDAP1, which shows significant similarities to BF2 in terms of secondary structure, translocates effectively across lipid vesicle and bacterial membranes, while the DesHDAP1 proline mutant shows reduced translocation abilities and antimicrobial potency. In contrast, both DesHDAP3 and its proline mutant translocate poorly, though the DesHDAP3 proline mutant is more potent. Our findings suggest that a proline hinge can promote membrane translocation in some peptides, but that the extent of its effect on permeabilization depends on the peptide's amphipathic properties. Our results also highlight the different antimicrobial mechanisms exhibited by histone-derived peptides and suggest that histones may serve as a source of novel antimicrobial peptides with varied properties.     

The intracellular mechanism of action on Escherichia coli of BF2-A/C, two analogues of the antimicrobial peptide Buforin 2

In the present study, the antimicrobial peptides BF2-A and BF2-C, two analogues of Buforin 2, were chemically synthesized and the activities were assayed. To elucidate the bactericidal mechanism of BF2-A/C and their different antimicrobial activities, the influence of peptides to E. coli cell membrane and targets of intracellular action were researched. Obviously, BF2-A and BF2-C did not induce the influx of PI into the E. coli cells, indicating nonmemebrane permeabilizing killing action. The FITC-labeled BF2-A/C could penetrate the E. coli cell membrane and BF2-C penetrated the cells more efficiently. Furthermore, BF2-A/C could bind to DNA and RNA respectively, and the affinity of BF2-C to DNA was powerful at least over 4 times than that of BF2-A. The present results implied that BF2-A and BF2-C inhibited the cellular functions by binding to DNA and RNA of cells after penetrating the cell membranes, resulting in the rapid cell death. The structure-activity relationship analysis of BF2-A/C revealed that the cell-penetrating efficiency and the affinity ability to DNA were critical factors for determining the antimicrobial potency of both peptides. The more efficient cell-penetrating and stronger affinity to DNA caused that BF2-C displayed more excellent antimicrobial activity and rapid killing kinetics than BF2-A.     

Insights into buforin II membrane translocation from molecular dynamics simulations

Buforin II is a histone-derived antimicrobial peptide that readily translocates across lipid membranes without causing significant membrane permeabilization. Previous studies showed that mutating the sole proline of buforin II dramatically decreases its translocation. As well, researchers have proposed that the peptide crosses membranes in a cooperative manner by forming transient toroidal pores. This paper reports molecular dynamics simulations designed to investigate the structure of buforin II upon membrane entry and evaluate whether the peptide is able to form toroidal pore structures. These simulations showed a relationship between protein–lipid interactions and increased structural deformations of the buforin N-terminal region promoted by proline. Moreover, simulations with multiple peptides show how buforin II can embed deeply into membranes and potentially form toroidal pores. Together, these simulations provide structural insight into the translocation process for buforin II in addition to providing more general insight into the role proline can play in antimicrobial peptides.     

抗菌肽BuforinⅡ衍生物与大肠杆菌基因组DNA的作用机制

【目的】研究抗菌肽BuforinⅡ的衍生物BF2-A/B与大肠杆菌基因组DNA的作用机制。【方法】琼脂糖电泳检测肽对DNA的断裂作用,凝胶阻滞实验研究肽与DNA的结合作用,圆二色谱考察结合肽后DNA结构的变化,荧光光谱分析肽与溴化乙锭竞争性嵌入DNA以及磷酸根对肽与DNA相互作用的影响。【结果】BF2-A/B不断裂基因组DNA而是结合DNA,使DNA双螺旋结构变得松散,削弱碱基对间的堆积作用,并取代EB,使EB-DNA复合体系荧光减弱。而PO43-的加入减弱了肽对DNA-EB荧光的淬灭作用。【结论】衍生肽与DNA的结合方式是先靠静电引力吸附到DNA磷酸基团上,随即插入双螺旋沟槽,嵌入碱基对间。BF2-B有更多的正电荷,更强的插入沟槽和嵌入碱基对的能力,使得其结合DNA的能力比BF2-A强。     

Membrane translocation mechanism of the antimicrobial peptide buforin 2

The antimicrobial peptide magainin 2 isolated from the skin of the African clawed frog Xenopus laevis crosses lipid bilayers by transiently forming a peptide-lipid supramolecular complex pore inducing membrane permeabilization and flip-flop of membrane lipids [Matsuzaki, K., Murase, O., Fujii, N., and Miyajima, K. (1996) Biochemistry 35, 11361-11368]. In contrast, the antimicrobial peptide buforin 2 discovered in the stomach tissue of the Asian toad Bufo bufo gargarizans efficiently crosses lipid bilayers without inducing severe membrane permeabilization or lipid flip-flop, and the Pro(11) residue plays a key role in this unique property [Kobayashi, S, Takeshima, K., Park, C. B., Kim, S. C., and Matsuzaki, K. (2000) Biochemistry 39, 8648-8654]. To elucidate the translocation mechanism, the secondary structure and the orientation of the peptide in lipid bilayers as well as the effects of the peptide concentration, the lipid composition, and the cis-trans isomerization of the Pro peptide bond on translocation efficiency were investigated. The translocation efficiencies of F10W-buforin 2 (BF2), P11A-BF2, and F5W-magainin 2 (MG2) across egg yolk L-alpha-phosphatidyl-DL-glycerol (EYPG)/egg yolk L-alpha-phosphatidylcholine (1/1) bilayers were dependent supralinearly on the peptide concentration, suggesting that the translocation mechanisms of these peptides are similar. The incorporation of the negative curvature-inducing lipid egg yolk L-alpha -phosphatidylethanolamine completely suppressed the translocation of BF2, indicating the induction of the positive curvature by BF2 on the membrane is related to the translocation process, similarly to MG2. In pure EYPG, where the repulsion between polycationic BF2 molecules is reduced, membrane permeabilization and coupling lipid flip-flop were clearly observed. Structural studies by use of Fourier transform infrared-polarized attenuated total reflection spectroscopy indicated that BF2 assumed distorted helical structures in EYPG/EYPC bilayers. A BF2 analogue with an alpha-methylproline, which fixed the peptide bond to the trans configuration, translocated similarly to the parent peptide, suggesting the cis-trans isomerization of the Pro peptide bond is not involved in the translocation process. These results indicate that BF2 crosses lipid bilayers via a mechanism similar to that of MG2. The presence of Pro(11) distorts the helix, concentrating basic amino acid residues in a limited amphipathic region, thus destabilizing the pore by enhanced electrostatic repulsion, enabling efficient translocation.     

Recognition of damaged regions in DNA by oligopeptides and proteins

The binding of various damaged DNAs to the single-strand binding protein coded for by gene 32 from bacteriophage T4, on the one hand, and of oligopeptides containing tryptophan and lysine residues, on the other hand, is described. These molecules exhibit a higher affinity for modified DNA than for native DNA in so far as modification results in a local destabilization of the double-stranded structure of the nucleic acid. Stacking interactions between aromatic amino acids and nucleic acid bases appear to play a crucial role in the recognition of destabilized regions induced by chemical agents (carcinogens and antitumor drugs). These interactions confer to the peptide lysyl-tryptophyl-lysine an endonucleolytic activity specific for apurinic sites. From results obtained with such oligopeptides a model for the active sites of Ap-endonucleases is proposed which could account for the strategy used by the denV endonuclease from phage T4 during the first step of excision repair of pyrimidine dimers in DNA. The effect of the overall conformation of modified DNA on repair efficiency is discussed.     

Quantitative approaches to monitor protein-nucleic acid interactions using fluorescent probes

Sequence-specific recognition of nucleic acids by proteins is required for nearly every aspect of gene expression. Quantitative binding experiments are a useful tool to measure the ability of a protein to distinguish between multiple sequences. Here, we describe the use of fluorophore-labeled oligonucleotide probes to quantitatively monitor protein/nucleic acid interactions. We review two complementary experimental methods, fluorescence polarization and fluorescence electrophoretic mobility shift assays, that enable the quantitative measurement of binding affinity. We also present two strategies for post-synthetic end-labeling of DNA or RNA oligonucleotides with fluorescent dyes. The approaches discussed here are efficient and sensitive, providing a safe and accessible alternative to the more commonly used radio-isotopic methods.     

RRM proteins interacting with the cap region of topoisomerase I

RNA recognition motif (RRM) domains bind both nucleic acids and proteins. Several proteins that contain two closely spaced RRM domains were previously found in protein complexes formed by the cap region of human topoisomerase I, a nuclear enzyme responsible for DNA relaxation or phosphorylation of SR splicing proteins. To obtain molecular insight into specific interactions between the RRM proteins and the cap region of topo I we examined their binary interactions using the yeast two-hybrid system. The interactions were established for hnRNP A1, p54(nrb) and SF2/ASF, but not for hnRNP L or HuR. To identify the amino acid pattern responsible for binding, experimental mutagenesis was employed and computational modelling of these processes was carried out. These studies revealed that two RRM domains and six residues of the consensus sequence are required for the binding to the cap region. On the basis of the above data, a structural model for the hnRNP A1-topoisomerase I complex was proposed. The main component of the hnRNP A1 binding site is a hydrophobic pocket on the beta-surface of the first RRM domain, similar to that described for Y14 protein interacting with Mago. We demonstrated that the interaction between RRM domains and the cap region was important for the kinase reaction catalyzed by topoisomerase I. Together with the previously described inhibitory effect of RRM domains of SF2/ASF on DNA cleavage, the above suggests that the binding of RRM proteins could regulate the activity of topoisomerase I.     

Solution NMR studies of amphibian antimicrobial peptides: Linking structure to function?

The high-resolution three-dimensional structure of an antimicrobial peptide has implications for the mechanism of its antimicrobial activity, as the conformation of the peptide provides insights into the intermolecular interactions that govern the binding to its biological target. For many cationic antimicrobial peptides the negatively charged membranes surrounding the bacterial cell appear to be a main target. In contrast to what has been found for other classes of antimicrobial peptides, solution NMR studies have revealed that in spite of the wide diversity in the amino acid sequences of amphibian antimicrobial peptides (AAMPs), they all adopt amphipathic α-helical structures in the presence of membrane-mimetic micelles, bicelles or organic solvent mixtures. In some cases the amphipathic AAMP structures are directly membrane-perturbing (e.g. magainin, aurein and the rana-box peptides), in other instances the peptide spontaneously passes through the membrane and acts on intracellular targets (e.g. buforin). Armed with a high-resolution structure, it is possible to relate the peptide structure to other relevant biophysical and biological data to elucidate a mechanism of action. While many linear AAMPs have significant antimicrobial activity of their own, mixtures of peptides sometimes have vastly improved antibiotic effects. Thus, synergy among antimicrobial peptides is an avenue of research that has recently attracted considerable attention. While synergistic relationships between AAMPs are well described, it is becoming increasingly evident that analyzing the intermolecular interactions between these peptides will be essential for understanding the increased antimicrobial effect. NMR structure determination of hybrid peptides composed of known antimicrobial peptides can shed light on these intricate synergistic relationships. In this work, we present the first NMR solution structure of a hybrid peptide composed of magainin 2 and PGLa bound to SDS and DPC micelles. The hybrid peptide adopts a largely helical conformation and some information regarding the inter-helix organization of this molecule is reported. The solution structure of the micelle associated MG2-PGLa hybrid peptide highlights the importance of examining structural contributions to the synergistic relationships but it also demonstrates the limitations in the resolution of the currently used solution NMR techniques for probing such interactions. Future studies of antimicrobial peptide synergy will likely require stable isotope-labeling strategies, similar to those used in NMR studies of proteins.     

Experimental and theoretical investigations of Dy(III) complex with 2,2′-bipyridine ligand: DNA and BSA interactions and antimicrobial activity study

Abstract

In this study, the interactions of a novel metal complex [Dy(bpy)₂Cl₃.OH₂] (bpy is 2,2'-bipyridine) with fish salmon DNA (FS-DNA) and bovine serum albumin (BSA) were investigated by experimental and theoretical methods. All results suggested significant binding between the Dy(III) complex with FS-DNA and BSA. The binding constants (K_b), Stern-Volmer quenching constants (K_SV) of Dy(III)-complex with FS-DNA and BSA at various temperatures as well as thermodynamic parameters using Van’t Hoff equation were obtained. The experimental results from absorption, ionic strength, iodide ion quenching, ethidium bromide (EtBr) quenching studies and positive ΔH? and ΔS? suggested that hydrophobic groove-binding mode played a predominant role in the binding of Dy(III)-complex with FS-DNA. Indeed, the molecular docking results for DNA-binding were in agreement with experimental data. Besides, the results found from experimental and molecular modeling indicated that the Dy(III)-complex bound to BSA via Van der Waals interactions. Moreover, the results of competitive tests by phenylbutazone, ibuprofen, and hemin (as a site-I, site-II and site-III markers, respectively) considered that the site-III of BSA is the most possible binding site for Dy(III)-complex. In addition, Dy(III) complex was concurrently screened for its antimicrobial activities. The presented data provide a promising platform for the development of novel metal complexes that target nucleic acids and proteins with antimicrobial activity.

Communicated by Ramaswamy H. Sarma     

Unwinding of unnatural substrates by a DNA helicase

Helicases separate double-stranded DNA into single-stranded DNA intermediates that are required during replication and recombination. These enzymes are believed to transduce free energy available from ATPase activity to unwind the duplex and translocate along the nucleic acid lattice. The nature of enzyme-substrate interactions between helicases and duplex DNA substrates has not been well-defined. Most helicases require a single-stranded DNA overhang adjacent to duplex DNA in order to initiate unwinding. The strand containing the overhang is referred to as the loading strand whereas the complementary strand is referred to as the displaced strand. We have investigated the interactions between a DNA helicase and the DNA substrate by replacing the displaced strand with a nucleic acid mimic, peptide nucleic acid (PNA). PNA is capable of forming duplex structures with DNA according to Watson-Crick base pairing rules, but contains a N-(2-aminoethyl)glycine backbone in place of the deoxyribose phosphates. The PNA-DNA hybrids had higher melting temperatures than their DNA-DNA counterparts. Dda helicase, from bacteriophage T4, was able to unwind the DNA-PNA substrates at similar rates as DNA-DNA substrates. The results indicate that the rate-limiting step for unwinding is relatively insensitive to the chemical nature of the displaced strand and the thermal stability of oligonucleotide substrates.     

Effect of the bases flanking an abasic site on the recognition of nucleobase by amiloride

Background

We explain here the various non-covalent interactions which are responsible for the different binding modes of a small ligand with DNA.

Methods

The combination of experimental and theoretical methods was used.

Results

The interaction of amiloride with thymine was found to depend on the bases flanking the AP site and different binding modes were observed for different flanking bases. Molecular modeling, absorption studies and binding constant measurements support for the different binding patterns. The flanking base dependent recognition of AP site phosphates was investigated by ³¹P NMR experiments. The thermodynamics of the ligand–nucleotide interaction was demonstrated by isothermal titration calorimetry. The emission behavior of amiloride was found to depend on the bases flanking the AP site. Amiloride photophysics in the context of AP-site containing DNA is investigated by time-dependent density functional theory.

Conclusions

Flanking bases affect the ground and excited electronic states of amiloride when binding to AP site, which causes flanking base-dependent fluorescence signaling.

General significance

The various noncovalent interactions have been well characterized for the determination of nucleic acid structure and dynamics, and protein–DNA interactions. However, these are not clear for the DNA–small molecule interactions and we believe that our studies will bring a new insight into such phenomena.     

Studies on cleavage of DNA by N-phosphoryl branched peptides

It was found that Nalpha,Nepsilon-di[N-(O,O-diisopropyl)phosphoryl-L-leucy]-L-lysyl-methyl ester (1) and Nalpha,Nepsilon-di[N-(O,O-diisopropyl)phosphoryl-L-phenylalanyl]-L-lysyl-methyl ester (2) could cleave supercoiled DNA such as PUC19 efficiently in 40 mM Britton-Robinson buffer. The cleavage activities for both were investigated by agarose gel electrophoresis. The T4 ligase experiments implied that the cleavage of DNA occurs via a hydrolytic path. The results showed that the cleavage reaction of DNA is dependent on the value of pH and ionic strength in the solution. DNA cleavage is more efficient by N-phosphoryl branched peptide 2 than by N-phosphoryl branched peptide 1. The experiments also show that hydrolysis of DNA by N-phosphoryl branched peptide 1 was accelerated in the presence of Mg2+ or Zn2+ ions. The interactions of DNA with N-phosphoryl branched peptides were also characterized by melting temperature measurements and circular dichroism (CD) techniques. On the basis of experimental data, the possible mechanism of interactions between DNA with N-phosphoryl branched peptides was discussed.     

Simulating the antimicrobial mechanism of human β-defensin-3 with coarse-grained molecular dynamics

Human β-defensin-3 (HβD-3) is an endogenous antimicrobial peptide with potent and broad killing activity against various microorganisms, and thus, it is an attractive candidate for the development of novel peptide antibiotics, but its antimicrobial mechanism remains elusive. To characterize the mechanism, we used multi-microsecond coarse-grained simulations with the MARTINI force field. These simulations show HβD-3 peptides can form oligomers on the surface of bacterial membrane and make anionic lipids (POPG) clustered. Furthermore, two kinds of regions (one is composed of pure POPG lipids, and the other is enriched in POPE lipids) are formed in the membrane; on the border of them, there are some obvious defects, which result in the membrane disruption. By contrast, the simulations also reveal that the contacts between the HβD-3 peptides and mammalian membrane are not stable. These results provide biophysical insights into HβD-3 selectivity and suggest a possible antimicrobial mechanism.     

An Antibody as Surrogate Receptor Reveals Determinants of Activity of an Innate Immune Peptide Antibiotic

Drug discovery initiatives often depend critically on knowledge of ligand-receptor interactions. However, the identity or structure of the target receptor may not be known in every instance. The concept of receptor surrogate, a molecular environment mimic of natural receptor, may prove beneficial under such circumstances. Here, we demonstrate the potential of monoclonal antibodies (mAbs) to act as surrogate receptors for a class of innate immune peptide antibiotics, a strategy that can help comprehend their action mechanism and identify chemical entities crucial for activity. A panel of antibody surrogates was raised against indolicidin, a tryptophan-rich cationic broad spectrum antimicrobial peptide of innate immune origin. Employing an elegant combination of thermodynamics, crystallography, and molecular modeling, interactions of the peptide with a high affinity anti-indolicidin monoclonal antibody were analyzed and were used to identify a motif that contained almost the entire antibiotic activity of native indolicidin. The analysis clarified the interaction of the peptide with previously proposed targets such as bacterial cell membrane and DNA and could further be correlated with antimicrobial compounds whose actions involve varied other mechanisms. These features suggest a multipronged assault pathway for indolicidin. Remarkably, the anti-indolicidin mAb surrogate was able to isolate additional independent bactericidal sequences from a random peptide library, providing compelling evidence as to the physiological relevance of surrogate receptor concept and suggesting applications in receptor-based pharmacophore research.     

Antimicrobial selectivity based on zwitterionic lipids and underlying balance of interactions

An important feature of antimicrobial peptides is their ability to distinguish pro- from eukaryotic membranes. In vitro experiments on the antimicrobial peptide NK-2 indicate that the discrimination between zwitterionic phosphatidylethanolamine lipids exposed by prokaryotes and phosphatidylcholine lipids exposed by eukaryotes plays an important role. The underlying mechanism is not understood. Here we present molecular dynamics simulations in conjunction with a coarse grained model and thermodynamic integration showing that NK-2 binds more strongly to palmitoyloleoylphosphatidylethanolamine (POPE) than to palmitoyloleoylphosphatidylcholine (POPC) bilayers. Finite size effects on the relative free energy have been corrected for with a method that may also be useful in future studies of the affinities of macromolecules for lipid membranes. Our results support the previous hypothesis that the stronger binding to PE compared to PC arises from a better accessibility of the phosphates of the lipids to the cationic peptide in a sense that a similar number of peptide-lipid salt bridges requires to break more favorable electrostatic headgroup-headgroup interactions for PC relative to PE. The transfer of NK-2 from POPC to POPE is found to lead to a decrease in electrostatic peptide-lipid but an increase in lipid-lipid and ion-lipid interactions, correlating with a dehydration of the lipids and the ions but an increased hydration of the peptide. The increase in affinity of NK-2 for POPE compared to POPC hence arises from a complex interplay of competing interactions. This work opens the perspective to study how the affinity of antimicrobial peptides changes with amino acid sequence and lipid composition.