Antimicrobial Peptide Simulations and the Influence of Force Field on the Free Energy for Pore Formation in Lipid Bilayers

ANTIMICROBIAL PEPTIDES — A POTENTIAL REPLACEMENT FOR TRADITIONAL ANTIBIOTICS

Antimicrobial peptides are a heterogeneous group of molecules involved in the innate and acquired immune response of various organisms, ranging from prokaryotes to mammals, including humans. They consist of 12–50 amino acid residues; have different physico-chemical and biological properties. The most common feature is their ability to destroy the prokaryotic cell membrane, which causes cell death. In the action, the molecules of antimicrobial peptides are embedded in the target bacteriological cells and change their conformation, forming structures in some cases resembling channels. Some other molecules of antimicrobial peptides can cover the surface of a bacteriological cell and form a carpet, when they reach a critical mass they act like detergents. In addition, being positively charged molecules of such peptides, penetrating through the membranes of parasitic and bacteriological cells, bind to polyanionic RNA and DNA molecules. Among the benefits of antimicrobial peptides is their high metabolic activity, low probability of occurrence of addictions and side effects. In addition, bacteriological pathogens that previously did not have resistance to any antimicrobial peptide are difficult to develop a strategy to control them. In this connection, these peptides are the most promising moleculessubstitutes for traditional antibiotics. The article discusses the approaches and strategies of therapeutic use, the studies of antimicrobial peptides identified in recent years; The most frequently encountered mechanisms of interaction of antimicrobial peptides and a bacteriological membrane are described, the physicochemical properties of peptide molecules are described; the results of studies on the detection of resistance of some strains of bacteria to antimicrobial peptides and antibiotics in general are summarized.

Antibacterial activity and its mechanisms of a recombinant Funme peptide against Cronobacter sakazakii in powdered infant formula

Cronobacter sakazakii (Cs) is a typical foodborne bacterium that infect powdered infant formula (PIF) worldwide. In this study, a recombinant antimicrobial peptide, branded as Funme peptide (FP)was applied to protect PIF from Cs contamination. The result from the antimicrobial activity assay showed that the minimum inhibitory concentration (MIC) of BMAP-27 peptide, FP and Ampicillin against Cs were 250.0, 125.0 and 15.6 μg/mL, respectively, indicating FP possessed higher MIC than that of Ampicillin, and lower MIC than that of BMAP-27. The minimum biofilm eradication concentration (MBEC) assay showed that FP at 2 × MIC (250.0 μg/mL) could completely eradicated Cs biofilms. The antibacterial activity of FP might be due to the increasing permeability and the release of cytoplasmic β-galactosidase of Cs. The results acquired from transmission electron microscopy and scanning electron microscopy indicated that FP induced the disruption and dysfunction of cell walls and membranes. Moreover, safety assay showed that FP had low cytotoxicity to human erythrocytes. The present study investigated the antibacterial effects and mechanisms of FP against Cs, providing promising evidence to apply this novel antimicrobial agent against Cs contamination in foods and food processing facilities.

Microscopic view of lipids and their diverse biological functions

Biological membranes and their diverse lipid constituents play key roles in a broad spectrum of cellular and physiological processes. Characterization of membrane-associated phenomena at a microscopic level is therefore essential to our fundamental understanding of such processes. Due to the semi-fluid and dynamic nature of lipid bilayers, and their complex compositions, detailed characterization of biological membranes at an atomic scale has been refractory to experimental approaches. Computational modeling and simulation offer a highly complementary toolset with sufficient spatial and temporal resolutions to fill this gap. Here, we review recent molecular dynamics studies focusing on the diversity of lipid composition of biological membranes, or aiming at the characterization of lipid-protein interaction, with the overall goal of dissecting how lipids impact biological roles of the cellular membranes.

Insights into Membrane Translocation of Protegrin Antimicrobial Peptides by Multistep Molecular Dynamics Simulations

Protegrin-1 (PG-1) is a cationic arginine-rich antimicrobial peptide. It is widely accepted that PG-1 induces membrane disruption by forming pores that lead to cell death. However, the insertion mechanism for these highly cationic peptides into the hydrophobic membrane environment is still poorly understood at the molecular scale. It has previously been determined that the association of arginine guanidinium and lipid phosphate groups results in strong bidentate bonds that stabilize peptide-lipid complexes. It has also been suggested that arginine residues are able to drag phosphate groups as they insert inside the membrane to form a toroidal pore. However, whether bidentate bonds play a significant role in inducing a pore formation remains unclear. To investigate the role of bidentate complexes in PG-1 translocation, we conducted molecular dynamics simulations. Two computational electroporation methods were implemented to examine the translocation process. We found that PG-1 could insert into the membrane, provided the external electric potential is large enough to first induce a water column or a pore within the lipid bilayer membrane. We also found that the highly charged PG-1 is capable in itself of inducing molecular electroporation. Substitution of arginines with charge-equivalent lysines showed a markedly reduced tendency for insertion. This indicates that the guanidinium group likely facilitates PG-1 translocation. Potential of mean force calculations suggests that peptide insertion inside the hydrophobic environment of the membrane core is not favored. We found that formation of a water column or a pore might be a prerequisite for PG-1 translocation. We also found that PG-1 can stabilize the pore after insertion. We suggest that PG-1 could be a pore inducer and stabilizer. This work sheds some light on PG-1 translocation mechanisms at the molecular level. Methods presented in this study may be extended to other arginine-rich antimicrobial and cell-penetrating peptides.

Single-Molecule Resolution of Antimicrobial Peptide Interactions with Supported Lipid A Bilayers

The molecular interactions between antimicrobial peptides (AMPs) and lipid A-containing supported lipid bilayers were probed using single-molecule total internal reflection fluorescence microscopy. Hybrid supported lipid bilayers with lipid A outer leaflets and phospholipid (1,2-dioleoyl-sn-glycero-3-phosphoethanolamine (DOPE)) inner leaflets were prepared and characterized, and the spatiotemporal trajectories of individual fluorescently labeled LL37 and Melittin AMPs were determined as they interacted with the bilayer surfaces comprising either monophosphoryl or diphosphoryl lipid A (from Escherichia coli) to determine the impact of electrostatic interactions. Large numbers of trajectories were obtained and analyzed to obtain the distributions of surface residence times and the statistics of the spatial trajectories. Interestingly, the AMP species were sensitive to subtle differences in the charge of the lipid, with both peptides diffusing more slowly and residing longer on the diphosphoryl lipid A. Furthermore, the single-molecule dynamics indicated a qualitative difference between the behavior of AMPs on hybrid Lipid A bilayers and on those composed entirely of DOPE. Whereas AMPs interacting with a DOPE bilayer exhibited two-dimensional Brownian diffusion with a diffusion coefficient of ∼1.7 μm2/s, AMPs adsorbed to the lipid A surface exhibited much slower apparent diffusion (on the order of ∼0.1 μm2/s) and executed intermittent trajectories that alternated between two-dimensional Brownian diffusion and desorption-mediated three-dimensional flights. Overall, these findings suggested that bilayers with lipid A in the outer leaflet, as it is in bacterial outer membranes, are valuable model systems for the study of the initial stage of AMP-bacterium interactions. Furthermore, single-molecule dynamics was sensitive to subtle differences in electrostatic interactions between cationic AMPs and monovalent or divalent anionic lipid A moieties.

Link between Membrane Composition and Permeability to Drugs

Prediction of membrane permeability to small molecules represents an important aspect of drug discovery. First-principles calculations of this quantity require an accurate description of both the thermodynamics and kinetics that underlie translocation of the permeant across the lipid bilayer. In this contribution, the membrane permeability to three drugs, or drug-like molecules, namely, 9-anthroic acid (ANA), 2',3'-dideoxyadenosine (DDA), and hydrocortisone (HYL), are estimated in a pure 1-palmitoyl-2-oleoylphosphatidylcholine (POPC) and in a POPC:cholesterol (2:1) mixture. On the basis of independent 2-5-μs free-energy calculations combined with a time-fractional Smoluchowski determination of the diffusivity, the estimated membrane permeabilities to these chemically diverse permeants fall within an order of magnitude from the experimental values obtained in egg-lecithin bilayers, with the exception of HYL in pure POPC. This exception is particularly interesting because the calculated permeability of the sterol-rich bilayer to HYL, in close agreement with the experimental value, is about 600 times lower than that of the pure POPC bilayer to HYL. In contrast, the permeabilities to ANA and DDA differ by less than a factor of 10 between the pure POPC and POPC:cholesterol bilayers. The unusual behavior of HYL, a large, amphiphilic compound, may be linked with the longer range perturbation of the lipid bilayer it induces, compared to ANA and DDA, suggestive of a possibly different translocation mechanism. We find that the tendency of lower permeabilities of the POPC:cholesterol bilayer relative to those of the pure POPC one is a consequence of increased free-energy barriers. Beyond reporting accurate estimates of the membrane permeability, the present contribution also demonstrates that rigorous free-energy calculations and a fractional-diffusion model are key in revealing the molecular phenomena linking the composition of a membrane to its permeability to drugs.

Implicit-solvent dissipative particle dynamics force field based on a four-to-one coarse-grained mapping scheme

A new set of efficient solvent-free dissipative particle dynamics (DPD) force fields was developed for phospholipids and peptides. To enhance transferability, this model maps around four heavy atoms and their connected hydrogen atoms into a coarse-grained elementary bead based on functional group. The effective hybrid potential between any pair of beads is composed of a short-range repulsive soft-core potential that directly adopts the form of an explicit-solvent DPD model and a long-range attractive hydrophobic potential. The parameters of the attractive potentials for lipid molecules were obtained by fitting the explicit-solvent DPD simulation of one bead of any type in a water box, then finely tuning it until the bilayer membrane properties obtained in the explicit-solvent model were matched. These parameters were further extended to amino acids according to bead type. The structural and elastic properties of bilayer membranes, free energy profiles for a lipid flip-flop and amino acid analogues translocating across the membrane, and membrane pore formation induced by antimicrobial peptides obtained from this solvent-free DPD force field considerably agreed with the explicit-solvent DPD results. Importantly, the efficiency of this method is guaranteed to accelerate the assembly of vesicles composed of several thousand lipids by up to 50-fold, rendering the experimental liposome dynamics as well as membrane-peptide interactions feasible at accessible computational expense.

Cooperative Modes of Action of Antimicrobial Peptides Characterized with Atomistic Simulations: A Study on Cecropin B

Antimicrobial peptides (AMPs) are widely occurring host defense agents of interest as one route for addressing the growing problem of multidrug-resistant pathogens. Understanding the mechanisms behind their antipathogen activity is instrumental in designing new AMPs. Herein, we present an all-atom molecular dynamics and free energy study on cecropin B (CB) and its constituent domains. We find a cooperative mechanism in which CB inserts into an anionic model membrane with its amphipathic N-terminal segment, supported by the hydrophobic C-terminal segment of a second peptide. The two peptides interact via a Glu···Lys salt bridge and together sustain a pore in the membrane. Using a modified membrane composition, we demonstrate that when the lower leaflet is overall neutral, insertion of the cationic segment is retarded and thus this mode of action is membrane specific. The observed mode of action utilizes a flexible hinge, a common structural motif among AMPs, which allows CB to insert into the membrane using either or both termini. Data from both unbiased trajectories and enhanced sampling simulations indicate that a requirement for CB to be an effective AMP is the interaction of its hydrophobic C-terminal segment with the membrane. Simulations of these segments in isolation reveal their aggregation in the membrane and a different mechanism of supporting pore formation. Together, our results show the complex interaction of different structural motifs of AMPs and, in particular, a potential role for electronegative side chains in an overall cationic AMP.

Molecular Dynamics Simulations Are Redefining Our View of Peptides Interacting with Biological Membranes

Ever since the first molecular mechanics computer simulations of biological molecules became possible, there has been the dream to study all complex biological phenomena in silico, simply bypassing the enormous experimental challenges and their associated costs. For this, two inherent requirements need to be met: First, the time scales achievable in simulations must reach up to the millisecond range and even longer. Second, the computational model must accurately reproduce what is measured experimentally. Despite some recent successes, the general consensus in the field to date has been that neither of these conditions have yet been met and that the dream will be realized, if at all, only in the distant future. In this Account, we show that this view is wrong; instead, we are actually in the middle of the in silico molecular dynamics (MD) revolution, which is reshaping how we think about protein function. The example explored in this Account is a recent advance in the field of membrane-active peptides (MAPs). MD simulations have succeeded in accurately capturing the process of peptide binding, folding, and partitioning into lipid bilayers as well as revealing how channels form spontaneously from polypeptide fragments and conduct ionic and other cargo across membranes, all at atomic resolution. These game-changing advances have been made possible by a combination of steadily advancing computational power, more efficient algorithms and techniques, clever accelerated sampling schemes, and thorough experimental verifications. The great advantage of MD is the spatial and temporal resolution, directly providing a molecular movie of a protein undergoing folding and cycling through a functional process. This is especially important for proteins with transitory functional states, such as pore-forming MAPs. Recent successes are demonstrated here for the large class of antimicrobial peptides (AMPs). These short peptides are an essential part of the nonadaptive immune system for many organisms, ubiquitous in nature, and of particular interest to the pharmaceutical industry in the age of rising bacterial resistance to conventional antibiotic treatments. Unlike integral membrane proteins, AMPs are sufficiently small to allow converged sampling with the unbiased high-temperature sampling methodology outlined here and are relatively easy to handle experimentally. At the same time, AMPs exhibit a wealth of complex and poorly understood interactions with lipid bilayers, which allow not only tuning and validation of the simulation methodology but also advancement of our knowledge of protein-lipid interactions at a fundamental level. Space constraints limit our discussion to AMPs, but the MD methodologies outlined here can be applied to all phenomena involving peptides in membranes, including cell-penetrating peptides, signaling peptides, viral channel forming peptides, and fusion peptides, as well as ab initio membrane protein folding and assembly. For these systems, the promise of MD simulations to predict the structure of channels and to provide complete-atomic-detail trajectories of the mechanistic processes underlying their biological functions appears to rapidly become a reality. The current challenge is to design joint experimental and computational benchmarks to verify and tune MD force fields. With this, MD will finally fulfill its promise to become an inexpensive, powerful, and easy-to-use tool providing atomic-detail insights to researchers as part of their investigations into membrane biophysics and beyond.

Computational studies of peptide-induced membrane pore formation

A variety of peptides induce pores in biological membranes; the most common ones are naturally produced antimicrobial peptides (AMPs), which are small, usually cationic, and defend diverse organisms against biological threats. Because it is not possible to observe these pores directly on a molecular scale, the structure of AMP-induced pores and the exact sequence of steps leading to their formation remain uncertain. Hence, these questions have been investigated via molecular modelling. In this article, we review computational studies of AMP pore formation using all-atom, coarse-grained, and implicit solvent models; evaluate the results obtained and suggest future research directions to further elucidate the pore formation mechanism of AMPs.This article is part of the themed issue 'Membrane pores: from structure and assembly, to medicine and technology'.

Recent computational studies of membrane interaction and disruption of human islet amyloid polypeptide: Monomers, oligomers and protofibrils

The amyloid deposits of human islet amyloid polypeptide (hIAPP) are found in type 2 diabetes patients. hIAPP monomer is intrinsically disordered in solution, whereas it can form amyloid fibrils both in vivo and in vitro. Extensive evidence suggests that hIAPP causes the disruption of cellular membrane, and further induces cytotoxicity and the death of islet β-cells in pancreas. The presence of membrane also accelerates the hIAPP fibril formation. hIAPP oligomers and protofibrils in the early stage of aggregation were reported to be the most cytotoxic, disrupting the membrane integrity and giving rise to the pathological process. The detailed molecular mechanisms of hIAPP-membrane interactions and membrane disruption are complex and remain mostly unknown. Here in this review, we focus on recent computational studies that investigated the interactions of full length and fragmentary hIAPP monomers, oligomers and protofibrils with anionic, zwitterionic and mixed anionic-zwitterionic lipid bilayers. We mainly discuss the binding orientation of monomers at membrane surface, the conformational ensemble and the oligomerization of hIAPP inside membranes, the effect of lipid composition on hIAPP oligomers/protofibrils-membrane interactions, and the hIAPP-induced membrane perturbation. This review provides mechanistic insights into the interactions between hIAPP and lipid bilayers with different lipid composition at an atomistic level, which is helpful to understand the hIAPP cytotoxicity mediated by membrane. This article is part of a Special Issue entitled: Protein Aggregation and Misfolding at the Cell Membrane Interface edited by Ayyalusamy Ramamoorthy.

Free Energy Calculations of Membrane Permeation: Challenges Due to Strong Headgroup-Solute Interactions

Understanding how different classes of molecules move across biological membranes is a prerequisite to predicting a solute's permeation rate, which is a critical factor in the fields of drug design and pharmacology. We use biased molecular dynamics computer simulations to calculate and compare the free energy profiles of translocation of several small molecules across 1,2-dioleoyl- sn-glycero-3-phosphocholine (DOPC) lipid bilayers as a first step toward determining the most efficient method for free energy calculations. We study the translocation of arginine, a sodium ion, alanine, and a single water molecule using the metadynamics, umbrella sampling, and replica exchange umbrella sampling techniques. Within the fixed lengths of our simulations, we find that all methods produce similar results for charge-neutral permeants, but not for polar or positively charged molecules. We identify the long relaxation time scale of electrostatic interactions between lipid headgroups and the solute to be the principal cause of this difference and show that this slow process can lead to an erroneous dependence of computed free energy profiles on the initial system configuration. We demonstrate the use of committor analysis to validate the proper sampling of the presumed transition state, which in our simulations is achieved only in replica exchange calculations. On the basis of these results we provide some useful guidance to perform and evaluate free energy calculations of membrane permeation.

Charged Antimicrobial Peptides Can Translocate across Membranes without Forming Channel-like Pores

How can highly charged, cationic antimicrobial peptides (AMPs) translocate across hydrophobic lipid bilayers despite the prohibitive energetic penalty to do so? A common explanation has been the formation of peptide-lined channels. However, for most AMPs, no structures of membrane pores have been found despite clear evidence of membrane leakage and antimicrobial activity. The study here suggests an alternative and simple reason: for the AMP PGLa from Xenopus laevis (charge +5), such pores are not needed to explain both leakage and peptide translocation. Elevated-temperature multimicrosecond equilibrium simulations at all-atomistic level reveal that peptides spontaneously translocate across the membrane individually on a timescale of tens of microseconds, without forming pores. Both surface-bound peptides and lipids assist in the one-by-one translocation of the charged side chains. Single peptides can remain in a transmembrane orientation for many microseconds, snorkeling some charged residues to one interface and some to the opposite, but without inducing a water channel. Instead of stable pores, short-lived water bridges occur when two or three peptides connect at their termini, allowing both ion translocation and lipid flip-flop via a brushlike mechanism usually involving the C terminus of one peptide. The results here suggest that for some specific antimicrobial and other membrane active peptides, pore formation may not have to be invoked at all to explain peptide translocation and membrane permeabilization, which may explain why no channel structures for them have been determined experimentally.

Probing the disparate effects of arginine and lysine residues on antimicrobial peptide/bilayer association

Antimicrobial peptides (AMPs) are key components of the innate immune response and represent promising templates for the development of broad-spectrum alternatives to conventional antibiotics. Most AMPs are short, cationic peptides that interact more strongly with negatively charged prokaryotic membranes than net neutral eukaryotic ones. Both AMPs and synthetic analogues with arginine-like side chains are more active against bacteria than those with lysine-like amine groups, though the atomistic mechanism for this increase in potency remains unclear. To examine this, we conducted comparative molecular dynamics simulations of a model negatively-charged membrane system interacting with two mutants of the AMP KR-12: one with lysine residues mutated to arginines (R-KR12) and one with arginine residues mutated to lysine (K-KR12). Simulations show that both partition analogously to the bilayer and display similar preferences for hydrogen bonding with the anionic POPGs. However, R-KR12 binds stronger to the bilayer than K-KR12 and forms significantly more hydrogen bonds, leading to considerably longer interaction times. Additional simulations with methylated R-KR12 and charge-modified K-KR12 mutants show that the extensive interaction seen in the R-KR12 system is partly due to arginine's strong atomic charge distribution, rather than being purely an effect of the greater number of hydrogen bond donors. Finally, free energy simulations reveal that both peptides are disordered in solution but form an amphipathic α-helix when inserted into the bilayer headgroup region. Overall, these results highlight the role of charge and hydrogen bond strength in peptide bilayer insertion, and offer potential insights for designing more potent analogues in the future.

Membrane localization and dynamics of geranylgeranylated Rab5 hypervariable region

The small GTPase Rab5 is a key regulator of endosomal trafficking processes and a marker for the early endosome. The C-terminal hypervariable region (HVR) of Rab5 is post-translationally modified at residues Cys²¹² and Cys²¹³ to accommodate two geranylgeranyl anchors (C20 carbon chain length) in order to associate Rab5 with the membrane. The structural role of the HVR regarding protein-early endosome membrane recruitment is not resolved due to its high degree of flexibility and lack of crystallographic information. Here, full-atomistic and coarse-grained molecular dynamics simulations of the truncated Rab5 HVR^206-215 in three model membranes of increasing complexity (pure phospholipid bilayer, ternary membrane with cholesterol, six-component early endosome) were performed. Specific electrostatic interactions between the HVR^206-215 Arg²⁰⁹ residue and the phosphate group of the inositol ring of PI(3)P were detected. This shows that PI(3)P acts as a first contact site of protein recruitment to the early endosome. The free energy change of HVR^206-215 extraction from the bilayer was largest for the physiological negatively charged membrane. 5μs coarse-grained simulations revealed an active recruitment of PI(3)P to the HVR^206-215 supporting the formation of Rab5- and PI(3)P enriched signaling platforms.

Probing a Continuous Polar Defect: A Reaction Coordinate for Pore Formation in Lipid Membranes

Various biophysical processes involve the formation of aqueous pores over lipid membranes, including processes of membrane fusion, antimicrobial peptide activity, lipid flip-flop, and membrane permeation. Reliable and efficient free-energy calculations of pore formation using molecular dynamics simulations remained challenging due to the lack of good reaction coordinates (RCs) for pore formation. We present a new RC for pore formation that probes the formation and rupture of a continuous polar defect over the membrane. Potential of mean force (PMF) calculations along the new RC rapidly converge and exhibit no hysteresis between pore-opening and pore-closing pathways, in contrast to calculations based on previous RCs. We show that restraints along the new RC may restrain the system tightly to the transition state of pore formation, rationalizing the absence of hysteresis. We observe that the PMF of pore formation in a tension-free membrane of dimyristoylphosphatidylcholine (DMPC) reveals a free-energy barrier for pore nucleation, confirming a long-hypothesized metastable prepore state. We test the influence of the lipid force field, the cutoff distance used for Lennard-Jones interactions, and the lateral membrane size on the free energies of pore formation. In contrast to PMF calculations based on previous RCs, we find that such parameters have a relatively small influence on the free energies of pore nucleation. However, the metastability of the open pore in DMPC may depend on such parameters. The RC has been implemented into an extension of the GROMACS simulation software. The new RC allows for reliable and computationally efficient free-energy calculations of pore formation in lipid membranes.