Search tips
Search criteria 


Logo of nihpaAbout Author manuscriptsSubmit a manuscriptHHS Public Access; Author Manuscript; Accepted for publication in peer reviewed journal;
Biochim Biophys Acta. Author manuscript; available in PMC 2013 June 27.
Published in final edited form as:
PMCID: PMC3694338

Characterization of a potent antimicrobial lipopeptide via coarse-grained molecular dynamics


The prevalence of antibiotic-resistant pathogens is a major medical concern, prompting increased interest in the development of novel antimicrobial compounds. One such set of naturally occurring compounds, known as antimicrobial peptides (AMPs), have broad-spectrum activity, but come with many limitations for clinical use. Recent work has resulted in a set of antimicrobial lipopeptides (AMLPs) with micromolar minimum inhibitory concentrations and excellent selectivity for bacterial membranes. To characterize a potent, synthetic lipopeptide, C16-KGGK, we used multi-microsecond coarse-grained simulations with the MARTINI forcefield, with a total simulation time of nearly 46 microseconds. These simulations show rapid binding of C16-KGGK, which forms micelles in solution, to model bacterial lipid bilayers. Furthermore, upon binding to the surface of the bilayer, these lipopeptides alter the local lipid organization by recruiting negatively charged POPG lipids to the site of binding. It is likely that this drastic reorganization of the bilayer has major effects on bilayer dynamics and cellular processes that depend on specific bilayer compositions. By contrast, the simulations revealed no association between the lipopeptides and model mammalian bilayers. These simulations provide biophysical insights into lipopeptide selectivity and suggest a possible mechanism for antimicrobial action.

Keywords: antimicrobial peptides, lipopeptides, molecular dynamics

1. Introduction

Antimicrobial peptides (AMPs) are naturally occurring compounds found in virtually all multicellular organisms, and serve as critical components of the innate immune system(1). First isolated from insects in the early 1980s(2) and from tree-frogs in the late 1980s(3), these compounds were quickly noted for their potency against bacte-ria(4,5). To date, over 1,700 antimicrobial peptides have been catalogued from a variety of species(6,7). These compounds tend to share a common set of characteristics, specifically a positive charge and an amphipathic structure(8).

AMPs have garnered increased interest over the last two decades as potential new drug candidates(9), as many inhibit bacterial growth at micromolar or even nanomolar concentrations(10). This potential lies in the fact that AMPs are less likely to induce evolved resistance in their targets, as they permeabilize the lipid bilayer by targeting the specific composition of the bacterial membrane and directly affecting membrane lipids(11). Their cationic nature provides selectivity to bacterial membranes, which are generally anionic due to large concentrations of negatively charged lipid species. Their amphipathic structure allows for binding and interaction with the lipid bilayer. Beyond this general hypothesis, the specific mechanism of membrane perturbation is unknown. A number of models have been suggested, including pore formation(12), detergent-like permeabilization of the bilayer(13), and membrane destabilization after AMPs coat the bilayer surface(14). Moreover, it seems likely that there is no single mechanism to explain AMP action; rather, different AMPs may be described by one or more of the above models.

Despite the optimism, most AMPs are not ideal drug candidates. Because of their size, they tend to be prohibitively expensive to produce in quantities large enough to be effective. Also, bioavailability is a concern, as peptidases would degrade free peptides in the body(15). With this in mind, the Shai group examined the effects of lipidating small peptides to construct easily-synthesized molecules with the same basic properties of AMPs. They demonstrated that conjugating fatty acids to short peptides that were in themselves membrane-inert could bestow membrane-active antimicrobial properties(16,17). In fairly recent work, the Shai lab developed a set of antimicrobial lipopep-tides (AMLPs) with broad-spectrum antibacterial and antifungal activity. They have a common architecture, including a fatty acid chain conjugated to a 4-residue peptide, of which 2 are lysines, yielding a net +2 charge. These molecules also include 1 D-amino acid to prevent peptidase activity (D-enantiomer denoted with bold lettering)(18). The most potent of these molecules, C16-KGGK (where C16 indicates a 16-carbon saturated fatty acid chain attached to the N-terminus), has a minimum inhibitory concentration (MIC) in the micromolar range for common strains of various pathogens. Further work showed C16-KGGK has a similar MIC against plant-pathogenic bacteria(19).

Probing the mechanism of action behind these specific AMLPs provides an approachable problem for molecular dynamics simulation. Great strides have already been made using MD to address problems in membrane biophysics(2022), including antimicrobial peptides(23). For many AMPs, sampling conformational space can prove challenging, especially if they are relatively unstructured in solution. C16-KGGK, however, is small, making conformational sampling computationally tractable. The limiting factor then becomes sampling the slow-timescale adjustments of the bilayer due to the presence of the lipopeptides.

To address this issue, coarse-graining (CG) approaches are frequently employed. In this kind of model, individual atoms are merged into larger particles that retain the underlying properties of the abstracted atoms. This dramatically reduces computational cost by reducing the number of degrees of freedom in the system, allowing for the use of significantly larger timesteps for the integration of the equations of motion. Recent work has shown success in simulating AMPs using the MARTINI coarse-grained force field, which boasts speeds two orders of magnitude greater than conventional all-atom models(24,25). This allows investigators to simulate many more systems, as well as larger systems, than are usually feasible in these kinds of studies.

In this study, we explore the selectivity and mechanism of activity of C16-KGGK using coarse-grained molecular dynamics simulation. Multiple multi-microsecond simulations employing the MARTINI force field were performed with two systems designed to serve as model replicas of “bacteria-like” and “mammal-like” bilayers. We demonstrate selectivity of C16-KGGK for bacterial bilayers and hypothesize that the effect is electrostatically driven. Further analysis of the bacterial systems shows some striking bilayer effects and implies a possible mechanism for the antibacterial activity of this AMLP.

2. Methods

2.1. System construction

The MARTINI coarse-grained force field was used to model our system of inter-est(26,27). Construction of C16-KGGK was done by merging the MARTINI models for lipid tails and peptides. Initial configurations for CG C16-KGGK were constructed by overlaying the CG model on an all-atom representation (see Figure 1). To model the peptide’s innate flexibility, no secondary structure restraints were applied.

Figure 1
Coarse-grained representation of the lipopeptide C16-KGGK (MARTINI model) overlaying a ball-and-stick all-atom representation. The non-polar saturated fatty acid tail is represented by gray spheres. Individual peptide backbone segments are represented ...

We simulated two unique membrane systems. The first represents a “gram-negative bacteria-like” bilayer, with a 2:1 phosphatidylethanolamine (POPE):phosphatidylglycerol (POPG) bilayer. This simplified representation and specific composition was chosen for its similarity to model membranes commonly used in experiment. The second was a “mammal-like” bilayer, with pure phosphatidylcholine (POPC). Each system totaled 480 lipids, split and distributed evenly between the two leaflets. The 2:1 ratio of POPE:POPG was maintained for each leaflet in the bacterial system. These bilayers were equilibrated for several hundred nanoseconds each before adding lipopeptides.

Before creating membrane simulations, we first built large water boxes with randomly distributed C16-KGGK lipopeptides. This system was run for several hundred nanoseconds, during which time the lipopeptides rapidly aggregated into micelles. A single large micelle was chosen arbitrarily and extracted. We removed enough lipopep-tides to reduce the number to 48, which allowed us to easily make systems with a 10:1 lipid to lipopeptide ratio.

For each bilayer, we created 4 unique systems by manually placing the 48 lipopep-tide micelle using VMD(28) such that the centroid of the micelle was between 60 and 80 Å from the centroid of the bilayer. We then added neutralizing sodium and chloride ions, with an excess to bring the concentration of free salt ions to approximately 100 mM.

To probe electrostatic selectivity, we created two more systems with the positive charges on the lysines removed, effectively rendering the peptides polar but chargeless. These systems are identical in composition to the POPE/POPG systems, with the exception of a minor change in the number of ions to maintain a system with a net zero charge.

All systems contained about 37,000 particles, including about 24,000 water beads, equivalent to about 300,000 total true atoms.

2.2. Simulation protocol

Our simulations were run using versions 4.0.5 and 4.5.4 of the GROMACS pack-age(29,30). We employed a time step of 10 fs as suggested by Winger et al, for accurate integration(31,32). Furthermore, the neighborlist was updated every 5 steps. We held the temperature at 300 K using Nose-Hoover temperature coupling(33,34). The pressure was treated semi-isotropically using the Parrinello-Rahman barostat(35), with a reference of 1 bar. Electrostatics were accounted for using a shift function with a coulomb cutoff of 12 Å. Shift was used for Van der Waals as well, with a switch distance of 9 Å and a cutoff of 12 Å.

We ran each simulation to an actual simulation time of at least 3 µs, though some trajectories were as long as 5 µs. Our total simulation time for this work, over all 10 systems, is about 46 µs. This can be thought of as an effective time of about 184 µs when we consider the suggested increase in kinetics due to the reduction of “friction” in our systems(27). All the times we report here are actual simulation times. Snapshots were saved every 100 ps and it is at this resolution that all analysis was performed.

2.3. Simulation analysis

All analysis was done using tools developed with the Lightweight Object Oriented Structure library (LOOS)(36). LOOS is an object-oriented structure library implemented in C++ and Boost, which provides a powerful library for creating new tools for the analysis of molecular dynamics simulations. LOOS is available for download at

2.3.1. Fractional contacts analysis

To assess the chemical composition of the environment immediately surrounding the lipopeptides, we used fractional contact analysis. Specifically, we counted the number of beads of lipid species and water within an 8 Å radius. We then report a time series of the fraction of those particles that were water, lipids and other lipopeptides. We ignored ions in this calculation, as they make up very few of the total number of atoms in the system.

2.3.2. Lipid density map

To measure the local enrichment of particular lipid species in the plane of the membrane, we generated a lipid density map. Specifically, for each frame we determined the centroid of each POPG molecule and created a 2D histogram using 1 Å2 bins, normalizing to produce density in lipids/Å2.

2.3.3. Radial distribution function

We computed the 2D radial distribution function (RDF) of various lipid species relative to lipopeptides in the membrane plane using the LOOS xy rdf tool(37). Each molecule was treated as a single point, located at its centroid. We also tracked the evolution of the RDF over time, breaking the trajectory into 10 ns windows.

2.3.4. Angle distribution

To assess the lateral ordering of lipid headgroups, we looked at the distribution of angles between sets of three POPG lipid headgroups (atom GL0, type P4 in the MARTINI model). For each frame of the simulation, the angle was computed for each combination of 3 atoms where the distance between the end atoms and the center atom was shorter than 8 Å.

3. Results

3.1. Selectivity for POPE/POPG bilayers

Choosing to simulate models for both bacterial and mammalian bilayers allowed us to investigate the selectivity of C16-KGGK for bacterial systems. Figure 2 shows the distance between the centroids of the AMLP micelle and the lipid bilayer, projected along the membrane normal. Part A shows the time series for each of the POPE/POPG (“bacterial”) simulations. In each case, the micelle binds to the surface of the membrane (distance of 35 Å) within the first few hundred nanoseconds. At this point, the micelle remains intact, with a well-formed hydrophobic core; we will henceforth refer to this as the “bound” state.

Figure 2
Distance, in Å, between the centroids of the lipopeptide micelle and the lipid bilayer along the membrane normal as a function of time for both (A) POPE/POPG and (B) POPC. The apparent plateau around 110 Å for the POPC systems reflects ...

In one of the four POPE/POPG simulations, the micelle actually inserts into the bilayer; this is seen in Figure 2A as the curve that drops below the bound state, with the final distance of 20 Å. After insertion, the micelle rapidly dissolves and the AMLPs disperse throughout the membrane, likely driven by AMLP-AMLP electrostatic repulsion. Figure 3 shows images of this process.

Figure 3
Snapshots from the POPE/POPG binding simulation. Highlighted are the (A) initial structure; (B) bilayer binding at 120 ns; (C) insertion at 220 ns; (D) and lateral dispersion at 490 ns. The micelle is colored green and red, for the peptide and fatty acid ...

By contrast, Figure 2B shows that in the POPC (“mammalian”) systems, the micelle does not approach—let alone bind—the membrane surface. The two simulations with a POPE/POPG bilayer and neutral lysines displayed a pattern identical to that for the POPC systems (data not shown). This makes intuitive sense, confirming that electrostatics are the primary driving force for AMLP-lipid binding. On the timescales simulated, the POPC bilayer is essentially invisible to the cationic micelle. For the POPE/POPG system, the net negative charge for the bilayer draws the micelle to the bound state only when the lipopeptides are cationic.

3.2. Preferential binding to POPG lipids

Considering the intermolecular contacts during the process of C16-KGGK micelle binding can provide insight into a mechanism for interaction and selectivity. Figure 4 highlights the evolution of these intermolecular contacts as a function of time as the AMLPs bind and insert into the membrane. The micelle begins completely surrounded by water; inter-AMLP contacts are not shown, but are at their greatest early in the trajectory when the lipopeptides are in the micelle. As the micelle approaches the membrane, water contacts are replaced by contacts to lipids. A two-step process of lipid contact is clearly visible, with the first step, around 100 ns, showing association with the bilayer and progression to the bound state. In the second step, at ~200 ns, complete insertion occurs and there is a drastic increase in contact with lipids in the bilayer.

Figure 4
Time series for the intermolecular contacts with C16-KGGK. Shows the POPE/POPG trajectory with complete lipopeptide insertion into the bilayer. The initial binding event occurs at about 120 ns, and insertion occurs at 200 ns. The missing remainder in ...

The significantly higher proportion of lipopeptides contacting POPG lipids, rather than POPE lipids, is particularly noteworthy. The difference is striking: in later portions of the trajectory, on average 40% of the contacts are to POPG, while about 25% are to POPE. This happens despite the fact that the relative proportion of POPE lipids in the bilayer is twice that of POPG. This demonstrates C16-KGGK’s strong affinity for POPG lipids.

3.3. Micelle binding locally enriches POPG

Figure 5 shows POPG density for one trajectory, after the micelle has bound (but not inserted). In this bound state, POPG lipids are recruited to the space directly beneath the micelle, forming a small patch of nearly pure POPG. When we compare the two leaflets, reorganization is clearly limited to the leaflet to which the lipopeptides bound (Part A) and not the distal leaflet (Part B).

Figure 5
Density of POPG lipids under the bound micelle, centered at the origin, for both the bound leaflet (A) and the unbound leaflet (B).

This effect can be further quantified through the use of a radial distribution function (RDF) in the plane of the bilayer, as seen in Figure 6. Here we show the probability density of POPG lipids, POPE lipids and the complete bilayer as a function of distance from lipopeptides. Again, we only consider time in the trajectory after binding. The curve for all lipids is nearly a straight line at 1, indicating that the headgroups are not induced to pack more tightly near the micelle binding site. However, the POPG curve indicates an enrichment of POPG lipids at distances less than 40 Å, while the POPE density is diminished; this is precisely what would be expected if the POPG were preferentially attracted to the AMLPs, locally displacing the POPE. The tail end of each curve shows the opposite trend, with an enrichment of POPE lipids at long distances. We do not believe this is due to repulsion between the POPE and the AMLPs, but rather simply due to the limited box size, as the POPE lipids displaced at short distances are pushed to the outer edges of the box. In ideal conditions, with a significantly larger box, both the POPE and POPG curves would decay to 1 to represent bulk concentrations.

Figure 6
Atomic density of POPE, POPG and all lipids as a function of distance from lipopeptides. Shown is the average of the bound state portions of the three trajectories that never inserted (with standard deviations for the three simulations).

The lateral reorganization of the bilayer becomes more evident when we consider the time dependence of the RDF between POPG and the lipopeptides. Figure 7 shows how the short-ranged enrichment of POPG evolves over time, both in the simulations where the micelle inserts (Part A) and in another where it merely remains bound to the membrane surface (Part B). Part A shows that the bilayer rapidly reorganizes as the micelle binds to the bilayer surface (the period roughly 100–200 ns into the trajectory), as POPG is recruited to the region under the micelle. Shortly after 200 ns insertion occurs, the POPG RDF drops dramatically, likely due to the dissipation of the lipopeptide micelle and the lateral dispersion of lipopeptides across the bilayer surface.

Figure 7
Lateral radial distribution function for POPG lipids as a function of distance from lipopeptides plotted over time with 10 ns windows. Shown are (A) a simulation with lipopeptide insertion and (B) a simulation that remains in the bound state throughout ...

For the bound simulation, on the other hand, the recruitment of POPG lipids to the micelle continues throughout the trajectory. After about 600 ns, the maximum limit of aggregation appears to have been reached, and changes in the RDF are likely due to local fluctuations.

4. Discussion and Conclusions

4.1. Mechanism of AMLP selectivity

In this work, the aim was to take advantage of the drastic speed increases available through coarse-grained molecular dynamics to characterize a potent antimicrobial lipopeptide, C16-KGGK. We ran a number of simulations to demonstrate the selectivity of C16-KGGK for bacteria-like model bilayers to begin to characterize the specific mechanism by which the molecule may inhibit bacterial growth.

One intuitive hypothesis, based on the simple physical properties of the molecule, is that electrostatics are the driving force for selectivity for the bilayer. Our simulations are in agreement with this hypothesis, as micelles of C16-KGGK showed rapid association with the bilayer in the systems with model 2:1 POPE:POPG bilayers. However, if the electrostatic attraction is removed, either by using neutral POPC lipids or artificially neutralizing the lysines, binding is abolished. We built these systems to be as similar in composition as possible, with the exception of the lipid composition and minor differences in the number of ions used for system charge neutralization. With all of this in mind, electrostatics are the most likely cause of selectivity, drawing the lipopeptides to the bacterial bilayer prior to antibacterial activity.

4.2. Bilayer effects

While understanding selectivity may be fairly straightforward, probing the mechanism of action for C16-KGGK is not trivial. The processes of binding to the bilayer and eventual insertion are complicated, with many degrees of freedom and multiple factors. Our analysis of simulations with the bacteria-like bilayers has elucidated some key features in the binding process that may provide clues to a possible mechanism of action.

One of the most striking features of the simulations are the distinction between the bound and inserted states. Of the four binding systems, only one actually fully inserted. After complete micelle insertion, the lipopeptides, likely due to electrostatic repulsion, spread laterally across the bilayer in the leaflet on which binding occurred. At high concentrations of lipopeptides, this is probably the method of antibacterial activity, as interactions between the lipopeptides would drive insertion and lateral dispersion, leading to detergent-like permeabilization of the bilayer.

However, at lower concentrations, such as the micromolar concentrations at which these lipopeptides are able to inhibit microbial growth, this bilayer permeabilization is unlikely to be the primary mode of action. Our simulations suggest an alternative mechanism of action, based on the relatively “stable” bound state seen in the other three simulations. In this state, a raft-like region of highly ordered POPG lipids beneath the micelle provides stability and affords the micelle an extended period of time to further rearrange the bilayer (see Section 4.3 below). Based on a small number of trajectories, it appears that insertion either occurs very quickly (soon after binding, before the full formation of the pure POPG phase under the micelle), or very slowly (longer than the microsecond timescale of the simulations).

4.3. Bilayer crystallization

As discussed above, the presence of the lipopeptide micelle in the bound state causes drastic reorganization of the lipid bilayer. This reorganization leads to a high local concentration of POPG lipids in the region of the bilayer directly below the micelle. Upon visual inspection of this nearly pure patch of POPG bilayer, it appears that the bilayer may be crystallizing, with the formation of a small raft with high order and packing. Figure 8 highlights this effect.

Figure 8
Crystallization effects as a result of micelle binding. (A) Micelle bound to the bilayer (Peptides in green, AMLP tails in red, POPG headgroups in pink, POPE headgroups in orange, and POPE/POPG phosphate groups in blue). (B) Visualization of headgroup ...

After roughly 500 ns in the bound state, the patch of bilayer beneath the micelle is composed entirely of POPG lipids (Part A). These lipids show ordered packing, with POPG headgroups aligned neatly into rows (Part B). This packing appears to be mediated by a similarly ordered interaction with the peptide portion of the lipopeptides (Part C); the peptides curl around the phosphate, allowing both charged lysine chains to rest between a pair of neighboring lipid phosphates.

We demonstrate the shift in structural organization as a result of micelle binding by showing the distribution of angles formed by sets of 3 headgroups from POPG lipids for the bound state in one trajectory (Part D). The curves corresponding to the neat bilayer (free of lipopeptides) and the distal leaflet in bound simulations are nearly identical. Most sets of POPG headgroups form angles within two main regions, around 60 degrees and 90 degrees. The distributions are broad and diffuse, indicating variability in positioning between headgroup beads.

The third curve, corresponding to the leaflet in the bound simulation to which the lipopeptide micelle binds, shows a very different distribution. The primary modes are now positioned around 90 degrees and nearly 180 degrees, as expected in a rectangular lattice. The peaks are very narrow compared to the other curves, demonstrating that a large number of triplet sets are forming these 90 and 180 degree angles. This confirms our visual hypothesis that the POPG lipids form an ordered rectangular lattice when the lipopeptide micelle is present.

Although we expected to see some ordering due to favorable peptide-lipid interactions, the degree of crystallization seems problematic. We suspect that this reflects a flaw in the force field. In contrast to an all-atom force field, where each particle has a partial charge, the only charged moieties in this systems are the lipid phosphates, POPE’s NH4, the lysine amines, and the free salt ions; water, which would ordinarily make strong polar interactions with these moieties, has zero charge in the MARTINI model. This, combined with the fact that MARTINI water is roughly four times larger than a water molecule, enhances the favorability of desolvating the lipid headgroups to create a lattice of oppositely charged phosphates and amines. Artifacts have been seen previously in simulations of charged peptides binding to bilayers using this model(38), suggesting that this may be a universal problem. For this reason, we plan to explore using the recently developed polarizable MARTINI water model in future work(39). Coarse-grained models with higher-resolution electrostatics models are also very in-triguing(40,41).

4.4. A possible mechanism of action

Recent work by many groups has described changes in bilayer reorganization and the localization of anionic lipids in response to the binding of cationic antimicrobial agents(42,43). These results, combined with the present simulations, suggest a number of possible modes of action for AMPs and AMLPs. First, the presence of POPE and POPG rich domains could destabilize the bilayer, either by rapidly creating boundary defects between domains or by causing local changes in bilayer tension and curvature due to the tendency of lipids with PE headgroups to promote high negative curva-ture(44). One could imagine anionic clusters also having an effect on membrane polarization, possibly resulting in slow cell leakage(45,46). Alternatively, the presence of bound AMLPs could alter the behavior of previously formed functional lipid domains. Membrane proteins that require certain local distributions of lipids may be impacted, as would proteins that are positioned along domain boundaries, which would be drastically redrawn(47). Also, cellular processes that require large-scale reorganization of specific lipids, such as cell division, would be impacted.

A scaffold for successful bilayer phase separation by an antimicrobial agent has been suggested by Epand and Epand(48). In their model, the molecule requires multiple cationic residues, significant hydrophobicity, and conformational flexibility. For C16-KGGK, these requirements are met, but micellization is likely an important part of the mechanism of action.

Given its structure, C16-KGGK is unlikely to exist as a monomer in solution except at very low concentrations; forming a micelle allows it to shield the hydrophobic acyl chain from water. Micelles also provide some level of stability necessary for C16-KGGK to bind to the bilayer surface and reorganize the bilayer without being immediately driven to insertion. Lastly, it allows C16-KGGK to cluster, creating a large structure of positive charge that, when bound to the surface, can create a large region of demixed POPG. Small clusters of the AMLP, or even individual molecules, if they were able to bind without inserting, would not be able to cause demixing at this level. Experimental work has shown that similar peptides are capable of forming macrostructures, including tubules and micelles(49).

4.5. Future work and conclusions

The present results suggest a number of avenues for future exploration. For example, it would be interesting to explore the differences in binding affinity and bound conformation between a lipopeptide such as C16-KGGK, which is potent but not very selective, other related lipopeptides that are selective (e.g. C16-KAAK), and bacteria-inert lipopeptides (e.g. C16-KLLK) . Given the structural similarity of these lipopeptides—in the MARTINI model, they are nearly identical—it seems likely that a higher-resolution approach would be required, such as all-atom modeling. It would also be interesting to try to elucidate the transition state of micelle binding and insertion using a combination of umbrella sampling and commitor analysis(50).

The simplicity of our chosen systems implies another avenue for future work. For instance, major constituents to the systems of interest, such as cholesterol in mam- malian bilayers, were omitted to maintain similarity to bilayers commonly used by experimental labs(18). Exploring the effects of other components of the bilayer and other ratios of POPE to POPG in bacterial systems to reflect the known diversity in bilayers would provide significant insights into the relative selectivity and potency of C16-KGGK.

We are optimistic that understanding the biophysical properties of these novel antimicrobial lipopeptides and their mechanism of action at micromolar concentrations is critical to developing more potent and selective AMP and AMLP based antibiotics.


  • C16-KGGK is a potent antimicrobial lipopeptide.
  • Coarse-grained molecular dynamics simulations yield detailed insights into peptide-membrane interactions.
  • Selectivity for bacterial membranes is due to favorable electrostatic interactions.
  • Lipopeptide binding reorganizes the membrane by recruiting anionic lipids to the binding site.


Thanks goes to the Center for Research Computing at the University of Rochester for providing the computational resources necessary for this research. This work and publication was made possible by grant number GM068411 from the Institutional Ruth L. Kirschstein National Research Service Award. Thanks also goes to Nick Leioatts for critical comments on this manuscript.


Publisher's Disclaimer: This is a PDF file of an unedited manuscript that has been accepted for publication. As a service to our customers we are providing this early version of the manuscript. The manuscript will undergo copyediting, typesetting, and review of the resulting proof before it is published in its final form. Please note that during the production process errors may be discovered which could affect the content, and all legal disclaimers that apply to the journal pertain.


1. Boman HG. Peptide antibiotics and their role in innate immunity. Annu Rev Immunol. 1995;13:61–92. [PubMed]
2. Boman HG, Steiner H. Humoral immunity in cecropia pupae. Curr Top Micro-biol Immunol. 1981;94–95:75–91. [PubMed]
3. Zasloff M. Magainins a class of antimicrobial peptides from xenopus skin: isolation characterization of two active forms and partial cDNA sequence of a precursor. Proc Natl Acad Sci U S A. 1987;84:5449–5453. [PubMed]
4. Westerhoff HV, Jureti D, Hendler RW, Zasloff M. Magainins and the disruption of membrane-linked free-energy transduction. Proc Natl Acad Sci U S A. 1989;86:6597–6601. [PubMed]
5. Mygind PH, Fischer RL, Schnorr KM, Hansen MT, Snksen CP, Lud-vigsen S, Ravents D, Buskov S, Christensen B, Maria LD, Taboureau O, Yaver D, Elvig-Jrgensen SG, Srensen MV, Christensen BE, Kjaerulff S, Frimodt-Moller N, Lehrer RI, Zasloff M, Kristensen H-H. Plectasin is a pep-tide antibiotic with therapeutic potential from a saprophytic fungus. Nature. 2005;437:975–980. [PubMed]
6. Wang Z, Wang G. APD: the antimicrobial peptide database. Nucleic Acids Res. 2004;32:D590–D592. [PMC free article] [PubMed]
7. Wang G, Li X, Wang Z. APD2: the updated antimicrobial peptide database and its application in peptide design. Nucleic Acids Res. 2009;37:D933–D937. [PMC free article] [PubMed]
8. Epand RM, Vogel HJ. Diversity of antimicrobial peptides and their mechanism of action. Biochim. Biophys Acta. 1999;1462:11–28. [PubMed]
9. Koczulla AR, Bals R. Antimicrobial peptides: current status and therapeutic potential. Drugs. 2003;63:389–406. [PubMed]
10. Zasloff M. Antimicrobial peptides of multicellular organisms. Nature. 2002;415:389–395. [PubMed]
11. Jenssen H, Hamill P, Hancock REW. Peptide antimicrobial agents. Clin Microbiol Rev. 2006;19:491–511. [PMC free article] [PubMed]
12. Brogden KA. Antimicrobial peptides: pore formers or metabolic inhibitors in bacteria? Nat Rev Micro. 2005;3:238–250. [PubMed]
13. Bechinger B, Lohner K. Detergent-like actions of linear amphipathic cationic antimicrobial peptides. Biochimica et Biophysica Acta (BBA) - Biomembranes. 2006;1758:1529–1539. [PubMed]
14. Shai Y, Oren Z. From “carpet” mechanism to de-novo designed diastereomeric cell-selective antimicrobial peptides. Peptides. 2001;22:1629–1641. [PubMed]
15. Hancock REW, Sahl H-G. Antimicrobial and host-defense peptides as new anti-infective therapeutic strategies. Nat Biotechnol. 2006;24:1551–1557. [PubMed]
16. Avrahami D, Shai Y. Bestowing antifungal and antibacterial activities by lipophilic acid conjugation to D,L-amino acid-containing antimicrobial peptides: a plausible mode of action. Biochemistry. 2003;42:14946–14956. [PubMed]
17. Avrahami D, Shai Y. A new group of antifungal and antibacterial lipopeptides derived from non-membrane active peptides conjugated to palmitic acid. J. Biol. Chem. 2004;279:12277–12285. [PubMed]
18. Makovitzki A, Avrahami D, Shai Y. Ultrashort antibacterial and antifungal lipopeptides. Proc Natl Acad Sci U S A. 2006;103:15997–16002. [PubMed]
19. Makovitzki A, Viterbo A, Brotman Y, Chet I, Shai Y. Inhibition of fungal and bacterial plant pathogens in vitro and in planta with ultrashort cationic lipopep-tides. Appl Environ Microbiol. 2007;73:6629–6636. [PMC free article] [PubMed]
20. Heller H, Schaefer M, Schulten K. Molecular dynamics simulation of a bilayer of 200 lipids in the gel and in the liquid-crystal phases. J. Phys. Chem. 1993;97:8343–8360.
21. Pastor RW, Feller SE. Time scales of lipid dynamics and molecular dynamics. In: Merz KM, Roux B Jr, editors. Biological Membranes: molecular A Computation perspective from Experiment. Birkhauser; 1996. pp. 3–30.
22. Lindahl E, Sansom MS. Membrane proteins: molecular dynamics simulations. Curr Opin Struct Biol. 2008;18:425–431. [PubMed]
23. Mtyus E, Kandt C, Tieleman DP. Computer simulation of antimicrobial pep-tides. Curr Med Chem. 2007;14:2789–2798. [PubMed]
24. Rzepiela AJ, Sengupta D, Goga N, Marrink SJ. Membrane poration by antimicrobial peptides combining atomistic and coarse-grained descriptions. Faraday Discuss. 2010;144:431–43. discussion 445-81. [PubMed]
25. Louhivuori M, Risselada HJ, van der Giessen E, Marrink SJ. Release of content through mechano-sensitive gates in pressurized liposomes. Proc Natl Acad Sci U S A. 2010;107:19856–19860. [PubMed]
26. Risselada JH, Mark AE, Marrink SJ. Application of mean field boundary potentials in simulations of lipid vesicles. J. Phys. Chem. B. 2008 [PubMed]
27. Marrink SJ, Risselada HJ, Yefimov S, Tieleman DP, de Vries AH. The MARTINI force field: coarse grained model for biomolecular simulations. J Phys Chem B. 2007;111:7812–7824. [PubMed]
28. Humphrey W, Dalke A, Schulten K. VMD: visual molecular dynamics. J Mol Graph. 1996;14:33–38. 27-8. [PubMed]
29. van der Spoel D, Lindahl E, Hess B, Groenhof G, Mark AE, Berend-sen HJC. GROMACS: fast flexible and free. J Comput Chem. 2005;26:1701–1718. [PubMed]
30. Hess B, Kutzner C, van der Spoel D, Lindahl E. GROMACS 4: Algorithms for highly efficient, load-balanced, and scalable molecular simulation. Journal of Chemical Theory and Computation. 2008;4:435–447.
31. Winger M, Trzesniak D, Baron R, van Gunsteren WF. On using a too large integration time step in molecular dynamics simulations of coarse-grained molecular models. Phys Chem Chem Phys. 2009;11:1934–1941. [PubMed]
32. Marrink SJ, Periole X, Tieleman DP, de Vries AH. Comment on “on using a too large integration time step in molecular dynamics simulations of coarsegrained molecular models” In: Winger M, Trzesniak D, Baron R, van Gunsteren WF, editors. Phys Chem. Chem. Phys. Vol. 11. 2009. p. 1934. Phys Chem Chem Phys 12 (2010) 2254-6, author reply 2257-8. [PubMed]
33. Nose S, Klein ML. Constant pressure molecular dynamics for molecular systems. Mol. Phys. 1983;50:1055–1076.
34. Hoover WG. Canonical dynamics: Equilibrium phase-space distributions. Phys. Rev. A. 1985;31:1695–1697. [PubMed]
35. Parrinello M, Rahman A. Polymorphic transitions in single crystals: A new molecular dynamics method. Journal of Applied Physics. 1981;52:7182–7190.
36. Romo TD, Grossfield A. LOOS: an extensible platform for the structural analysis of simulations. Conf Proc IEEE Eng Med Biol Soc. 2009;2009:2332–2335. [PubMed]
37. Romo TD, Grossfield A. LOOS: A lightweight object-oriented software library. 2010
38. Vorobyov I, Li L, Allen TW. Assessing atomistic and coarse-grained force fields for protein-lipid interactions: the formidable challenge of an ionizable side chain in a membrane. J Phys Chem B. 2008;112:9588–9602. [PubMed]
39. Yesylevskyy SO, Schfer LV, Sengupta D, Marrink SJ. Polarizable water model for the coarse-grained MARTINI force field. PLoS Comput Biol. 2010;6:e1000810. [PMC free article] [PubMed]
40. Golubkov PA, Ren P. Generalized coarse-grained model based on point multi-pole and Gay-Berne potentials. J Chem Phys. 2006;125:64103. [PubMed]
41. Golubkov PA, Wu JC, Ren P. A transferable coarse-grained model for hydrogen-bonding liquids. Phys Chem Chem Phys. 2008;10:2050–2057. [PMC free article] [PubMed]
42. Jean-Franois F, Castano S, Desbat B, Odaert B, Roux M, Metz-Boutigue M-H, Dufourc EJ. Aggregation of cateslytin beta-sheets on negatively charged lipids promotes rigid membrane domains. a new mode of action for antimicrobial pep-tides? Biochemistry. 2008;47:6394–6402. [PubMed]
43. Epand RF, Maloy L, Ramamoorthy A, Epand RM. Amphipathic helical cationic antimicrobial peptides promote rapid formation of crystalline states in the presence of phosphatidylglycerol: lipid clustering in anionic membranes. Biophys J. 2010;98:2564–2573. [PubMed]
44. Arouri A, Dathe M, Blume A. Peptide induced demixing in PG/PE lipid mixtures: a mechanism for the specificity of antimicrobial peptides towards bacterial membranes? Biochim Biophys Acta. 2009;1788:650–659. [PubMed]
45. Radzishevsky IS, Rotem S, Bourdetsky D, Navon-Venezia S, Carmeli Y, Mor A. Improved antimicrobial peptides based on acyl-lysine oligomers. Nat Biotechnol. 2007;25:657–659. [PubMed]
46. Rotem S, Radzishevsky IS, Bourdetsky D, Navon-Venezia S, Carmeli Y, Mor A. Analogous oligo-acyl-lysines with distinct antibacterial mechanisms. FASEB J. 2008;22:2652–2661. [PubMed]
47. Epand RM, Epand RF. Bacterial membrane lipids in the action of antimicrobial agents. J Pept Sci. 2011;17:298–305. [PubMed]
48. Epand RM, Epand RF. Lipid domains in bacterial membranes and the action of antimicrobial agents. Biochim Biophys Acta. 2009;1788:289–294. [PubMed]
49. Makovitzki A, Baram J, Shai Y. Antimicrobial lipopolypeptides composed of palmitoyl di- and tricationic peptides: In vitro and in vivo activities, self-assembly to nanostructures, and a plausible mode of action. Biochemistry. 2008 [PubMed]
50. Kasson PM, Lindahl E, Pande VS. Atomic-resolution simulations predict a transition state for vesicle fusion defined by contact of a few lipid tails. PLoS Comput Biol. 2010;6:e1000829. [PMC free article] [PubMed]