Determining the orientation of protegrin-1 in DLPC bilayers using an implicit solvent-membrane model

Improving the antimicrobial potential of the peptide CIDEM-501 through acylation: A computational approach

Acylation is a common method used to modify antimicrobial peptides to enhance their effectiveness. It increases the interactions between the peptide and the bacterial cell membranes. However, acylation can also reduce the selectivity of the peptides by making them more active on eukaryotic membranes, which can lead to unintended toxicity. This study examines the potential of using in silico tools to evaluate the interaction and selectivity of the antimicrobial peptide CIDEM-501 when acylated with decanoic acid at the N-terminus, compared to the non-acylated counterpart. Circular dichroism, microdilution, and hemolysis assays were used to determine the peptide's secondary structure, antimicrobial activity, and selectivity to validate the theoretical predictions. The acylated peptide showed a more stable interaction with the bacterial membrane by inserting the acyl chain into the membrane's hydrophobic core, which led to tighter adsorption and a greater buried surface area. Additionally, it significantly altered membrane order more than the non-acylated counterpart, suggesting superior antimicrobial potential. Finally, in vitro activity assays confirmed theoretical predictions, showing that the acylated peptide had lower Minimum Inhibitory Concentration (MIC) values than the non-acylated peptide. Neither peptide showed significant hemolytic activity at their MIC. The computational techniques used in this study displayed strong predictive capability and helped to elucidate the interaction between the peptide and the membranes.

Novel antimicrobial peptides based on Protegrin-1: In silico and in vitro assessments

The development of antibiotic resistance has caused significant health problems. Antimicrobial peptides (AMPs) are considered next-generation antibiotics. Protegrin-1 (PG-1) is a β-hairpin AMP with a membrane-binding capacity. This study used twelve PG-1 analogs with different amino acid substitutions. Coarse-grained molecular dynamics (MD) simulations were used to assess these analogs, and their physicochemical properties were computed using the Antimicrobial Peptide Database. Three AMPs, PEP-D, PEP-C, and PEP-H, were chosen and synthesized for antibacterial testing. The microbroth dilution technique and hemolytic assays evaluated the antimicrobial efficacy and cellular toxicity. The checkerboard method was used to test the combined activity of AMP and standard antibiotics. Cell membrane permeability and electron microscopy were used to evaluate the mode of action. The chemical stability of the selective AMP, PEP-D, was assessed by a validated HPLC method. PEP-D consists of 16-18 amino acid residues and has a charge of +7 and a hydrophobicity of 44 %, similar to PG-1. It can efficiently inactivate bacteria by disrupting cell membranes and significantly reducing hemolytic activity. Chemical stability studies indicated that AMP was stable at 40 °C for six months under autoclave conditions. This study could introduce the potential therapeutic application of selective AMP as an anti-infective agent.

Genomewide Analysis of the Antimicrobial Peptides in Python bivittatus and Characterization of Cathelicidins with Potent Antimicrobial Activity and Low Cytotoxicity

In this study, we sought to identify novel antimicrobial peptides (AMPs) in Python bivittatus through bioinformatic analyses of publicly available genome information and experimental validation. In our analysis of the python genome, we identified 29 AMP-related candidate sequences. Of these, we selected five cathelicidin-like sequences and subjected them to further in silico analyses. The results showed that these sequences likely have antimicrobial activity. The sequences were named Pb-CATH1 to Pb-CATH5 according to their sequence similarity to previously reported snake cathelicidins. We predicted their molecular structure and then chemically synthesized the mature peptide for three putative cathelicidins and subjected them to biological activity tests. Interestingly, all three peptides showed potent antimicrobial effects against Gram-negative bacteria but very weak activity against Gram-positive bacteria. Remarkably, ΔPb-CATH4 showed potent activity against antibiotic-resistant clinical isolates and also was observed to possess very low hemolytic activity and cytotoxicity. ΔPb-CATH4 also showed considerable serum stability. Electron microscopic analysis indicated that ΔPb-CATH4 exerts its effects via toroidal pore preformation. Structural comparison of the cathelicidins identified in this study to previously reported ones revealed that these Pb-CATHs are representatives of a new group of reptilian cathelicidins lacking the acidic connecting domain. Furthermore, Pb-CATH4 possesses a completely different mature peptide sequence from those of previously described reptilian cathelicidins. These new AMPs may be candidates for the development of alternatives to or complements of antibiotics to control multidrug-resistant pathogens.

Mapping activity elements of protegrin antimicrobial peptides by HomoSAR

HomoSAR has been able to shed light on the relationship between sequences of protegrin peptides and their activity on six specific micro-organisms.

Cobra cytotoxins: structural organization and antibacterial activity

Cardiotoxins (cytotoxins, CT) are β-structured proteins isolated from the venom of cobra. They consist of 59-61 amino acid residues, whose antiparallel chains form three 'fingers'. In contrast to neurotoxins with an overall similar fold, CTs are amphiphilic. The amphiphilicity is caused by positively charged lysine and arginine residues flanking the tips of the loops that consist primarily of hydrophobic amino acids. A similar distribution of amino acid residues is typical for linear (without disulfide bonds) cationic cytolytic peptides from the venoms of other snakes and insects. Many of them are now considered to be lead compounds in combatting bacterial infections and cancer. In the present review, we summarize the data on the antibacterial activity of CTs and compare it to the activity of linear peptides.

Membrane interactions and pore formation by the antimicrobial peptide protegrin

Protegrin is an antimicrobial peptide with a β-hairpin structure stabilized by a pair of disulfide bonds. It has been extensively studied by solid-state NMR and computational methods. Here we use implicit membrane models to examine the binding of monomers on the surface and in the interior of the membrane, the energetics of dimerization, the binding to membrane pores, and the stability of different membrane barrel structures in pores. Our results challenge a number of conclusions based on previous experimental and theoretical work. The burial of monomers into the membrane interior is found to be unfavorable for any membrane thickness. Because of its imperfect amphipathicity, protegrin binds weakly, at most, on the surface of zwitterionic membranes. However, it binds more favorably onto toroidal pores. Anionic charge on the membrane facilitates the binding due to electrostatic interactions. Solid-state NMR results have suggested a parallel NCCN association of monomers in dimers and association of dimers to form octameric or decameric β-barrels. We find that this structure is not energetically plausible for binding to bilayers, because in this configuration the hydrophobic sides of two monomers point in opposite directions. In contrast, the antiparallel NCCN and especially the parallel NCNC octamers are stable and exhibit a favorable binding energy to the pore. The results of 100-ns simulations in explicit bilayers corroborate the higher stability of the parallel NCNC barrel compared with the parallel NCCN barrel. The ability to form pores in zwitterionic membranes provides a rationalization for the peptide's cytotoxicity. The discrepancies between our results and experiment are discussed, and new experiments are proposed to resolve them and to test the validity of the models.

Multiscale models of the antimicrobial peptide protegrin-1 on gram-negative bacteria membranes

Antimicrobial peptides (AMPs) are naturally-occurring molecules that exhibit strong antibiotic properties against numerous infectious bacterial strains. Because of their unique mechanism of action, they have been touted as a potential source for novel antibiotic drugs. We present a summary of computational investigations in our lab aimed at understanding this unique mechanism of action, in particular the development of models that provide a quantitative connection between molecular-level biophysical phenomena and relevant biological effects. Our work is focused on protegrins, a potent class of AMPs that attack bacteria by associating with the bacterial membrane and forming transmembrane pores that facilitate the unrestricted transport of ions. Using fully atomistic molecular dynamics simulations, we have computed the thermodynamics of peptide-membrane association and insertion, as well as peptide aggregation. We also present a multi-scale analysis of the ion transport properties of protegrin pores, ranging from atomistic molecular dynamics simulations to mesoscale continuum models of single-pore electrodiffusion to models of transient ion transport from bacterial cells. Overall, this work provides a quantitative mechanistic description of the mechanism of action of protegrin antimicrobial peptides across multiple length and time scales.

Thermodynamic analysis of protegrin-1 insertion and permeation through a lipid bilayer

Molecular dynamics (MD) simulations are used to study the pathway for the insertion of the cationic antimicrobial peptide protegrin-1 (PG1) into mixed anionic lipid bilayers composed of palmitoyl-oleoyl-phosphatidylglycerol (POPG) and palmitoyl-oleoyl-phosphatidylethanolamine (POPE) in a 1:3 ratio (POPG/POPE). We calculate the potential of mean force (PMF) during the transfer of the peptide from the bulk aqueous phase to the transmembrane (TM) configuration using the adaptive biasing force (ABF) method. We find that the PMF has two energy minima separated by an energy barrier. One minimum corresponds to the fully transmembrane inserted state, with a free energy of -20.1 kcal/mol. The second PMF minimum, which corresponds to adsorption to the membrane surface, has a value of -2.5 kcal/mol. The PMF also shows the existence of a free energy barrier of +6.3 kcal/mol for the insertion process. Using the Kramers theory Langevin equation and the Grote-Hynes theory generalized Langevin equation, we calculated the transmission coefficient for PG1 diffusion through the potential barrier. We focus on the use of the PMF and the time correlation function of the fluctuation of the instantaneous force to calculate the rate constants for insertion/deinsertion of PG1 from the mixed POPG/POPE membrane. The influence of the activation free energy barrier on the dynamics of the insertion and permeation of peptides through the membrane are discussed.

Computational studies of protegrin antimicrobial peptides: a review

Antimicrobial peptides (AMPs) are small, naturally occurring peptides that exhibit strong antibacterial properties generally believed to be a result of selective bacterial membrane disruption. As a result, there has been significant interest in the development of therapeutic antibiotics based on AMPs; however, the poor understanding of the fundamental mechanism of action of these peptides has largely hampered such efforts. We present a summary of computational and theoretical investigations of protegrin, a particularly potent peptide that is both an excellent model for the mechanism of action of AMPs and a promising therapeutic candidate. Experimental investigations have shed light on many of the key steps in the action of protegrin: protegrin monomers are known to dimerize in various lipid environments; protegrin peptides interact strongly with lipid bilayer membranes, particularly anionic lipids; protegrins have been shown to form pores in lipid bilayers, which results in uncontrolled ion transport and may be a key factor in bacterial death. In this work, we present a comprehensive review of the computational and theoretical studies that have complemented and extended the information obtained from experimental work with protegrins, as well as a brief survey of the experimental biophysical studies that are most pertinent to such computational work. We show that a consistent, mechanistic description of the bactericidal mechanism of action of protegrins is emerging, and briefly outline areas where the current understanding is deficient. We hope that the research reviewed herein offers compelling evidence of the benefits of computational investigations of protegrins and other AMPs, as well as providing a useful guide to future work in this area.

Dimerization of protegrin-1 in different environments

The dimerization of the cationic β-hairpin antimicrobial peptide protegrin-1 (PG1) is investigated in three different environments: water, the surface of a lipid bilayer membrane, and the core of the membrane. PG1 is known to kill bacteria by forming oligomeric membrane pores, which permeabilize the cells. PG1 dimers are found in two distinct, parallel and antiparallel, conformations, known as important intermediate structural units of the active pore oligomers. What is not clear is the sequence of events from PG1 monomers in solution to pores inside membranes. The step we focus on in this work is the dimerization of PG1. In particular, we are interested in determining where PG1 dimerization is most favorable. We use extensive molecular dynamics simulations to determine the potential of mean force as a function of distance between two PG1 monomers in the aqueous subphase, the surface of model lipid bilayers and the interior of these bilayers. We investigate the two known distinct modes of dimerization that result in either a parallel or an antiparallel β-sheet orientation. The model bilayer membranes are composed of anionic palmitoyl-oleoyl-phosphatidylglycerol (POPG) and palmitoyl-oleoyl-phosphatidylethanolamine (POPE) in a 1:3 ratio (POPG:POPE). We find the parallel PG1 dimer association to be more favorable than the antiparallel one in water and inside the membrane. However, we observe that the antiparallel PG1 β-sheet dimer conformation is somewhat more stable than the parallel dimer association at the surface of the membrane. We explore the role of hydrogen bonds and ionic bridges in peptide dimerization in the three environments. Detailed knowledge of how networks of ionic bridges and hydrogen bonds contribute to peptide stability is essential for the purpose of understanding the mechanism of action for membrane-active peptides as well as for designing peptides which can modulate membrane properties. The findings are suggestive of the dominant pathways leading from individual PG1 molecules in solution to functional pores in bacterial membranes.

Antimicrobial protegrin-1 forms ion channels: molecular dynamic simulation, atomic force microscopy, and electrical conductance studies

Antimicrobial peptides (AMPs) are an emerging class of antibiotics for controlling health effects of antibiotic-resistant microbial strains. Protegrin-1 (PG-1) is a model antibiotic among beta-sheet AMPs. Antibiotic activity of AMPs involves cell membrane damage, yet their membrane interactions, their 3D membrane-associated structures and the mechanism underlying their ability to disrupt cell membrane are poorly understood. Using complementary approaches, including molecular dynamics simulations, atomic force microscopy (AFM) imaging, and planar lipid bilayer reconstitution, we provide computational and experimental evidence that PG-1, a beta-hairpin peptide, forms ion channels. Simulations indicate that PG-1 forms channel-like structures with loosely attached subunits when reconstituted in anionic lipid bilayers. AFM images show the presence of channel-like structures when PG-1 is reconstituted in dioleoylphosphatidylserine/palmitoyloleoyl phosphatidylethanolamine bilayers or added to preformed bilayers. Planar lipid bilayer electrical recordings show multiple single channel conductances that are consistent with the heterogeneous oligomeric channel structures seen in AFM images. PG-1 channel formation seems to be lipid-dependent: PG-1 does not easily show ion channel electrical activity in phosphatidylcholine membranes, but readily shows channel activity in membranes rich in phosphatidylethanolamine or phosphatidylserine. The combined results support a model wherein the beta-hairpin PG-1 peptide acts as an antibiotic by altering cell ionic homeostasis through ion channel formation in cell membranes.

Relative free energy of binding between antimicrobial peptides and SDS or DPC micelles

We present relative binding free energy calculations for six antimicrobial peptide-micelle systems, three peptides interacting with two types of micelles. The peptides are the scorpion derived antimicrobial peptide (AMP), IsCT and two of its analogues. The micelles are dodecylphosphatidylcholine (DPC) and sodium dodecylsulphate (SDS) micelles. The interfacial electrostatic properties of DPC and SDS micelles are assumed to be similar to those of zwitterionic mammalian and anionic bacterial membrane interfaces, respectively. We test the hypothesis that the binding strength between peptides and the anionic micelle SDS can provide information on peptide antimicrobial activity, since it is widely accepted that AMPs function by binding to and disrupting the predominantly anionic lipid bilayer of the bacterial cytoplasmic membrane. We also test the hypothesis that the binding strength between peptides and the zwitterionic micelle DPC can provide information on peptide haemolytic activities, since it is accepted that they also bind to and disrupt the zwitterionic membrane of mammalian cells. Equilibrium structures of the peptides, micelles and peptide-micelle complexes are obtained from more than 300 ns of molecular dynamics simulations. A thermodynamic cycle is introduced to compute the binding free energy from electrostatic, non-electrostatic and entropic contributions. We find relative binding free energy strengths between peptides and SDS to correlate with the experimentally measured rankings for peptide antimicrobial activities, and relative free energy binding strengths between peptides and DPC to correlate with the observed rankings for peptide haemolytic toxicities. These findings point to the importance of peptide-membrane binding strength for antimicrobial activity and haemolytic activity.

Coupling molecular dynamics simulations with experiments for the rational design of indolicidin-analogous antimicrobial peptides

Antimicrobial peptides (AMPs) have attracted much interest in recent years because of their potential use as new-generation antibiotics. Indolicidin (IL) is a 13-residue cationic AMP that is effective against a broad spectrum of bacteria, fungi, and even viruses. Unfortunately, its high hemolytic activity retards its clinical applications. In this study, we adopted molecular dynamics (MD) simulations as an aid toward the rational design of IL analogues exhibiting high antimicrobial activity but low hemolysis. We employed long-timescale, multi-trajectory all-atom MD simulations to investigate the interactions of the peptide IL with model membranes. The lipid bilayer formed by the zwitterionic 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) was chosen as the model erythrocyte membrane; lipid bilayers formed from a mixture of POPC and the negatively charged 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoglycerol were chosen to model bacterial membranes. MD simulations with a total simulation time of up to 4 micros revealed the mechanisms of the processes of IL adsorption onto and insertion into the membranes. The packing order of these lipid bilayers presumably correlated to the membrane stability upon IL adsorption and insertion. We used the degree of local membrane thinning and the reduction in the order parameter of the acyl chains of the lipids to characterize the membrane stability. The order of the mixed 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoglycerol/POPC lipid bilayer reduced significantly upon the adsorption of IL. On the other hand, although the order of the pure-POPC lipid bilayer was perturbed slightly during the adsorption stage, the value was reduced more dramatically upon the insertion of IL into the membrane's hydrophobic region. The results imply that enhancing IL adsorption on the microbial membrane may amplify its antimicrobial activity, while the degree of hemolysis may be reduced through inhibition of IL insertion into the hydrophobic region of the erythrocyte membrane. In addition, through simulations, we identified the amino acids that are most responsible for the adsorption onto or insertion into the two model membranes. Positive charges are critical to the peptide's adsorption, whereas the presence of hydrophobic Trp8 and Trp9 leads to its deeper insertion. Combining the hypothetical relationships between the membrane disordering and the antimicrobial and hemolytical activities with the simulated results, we designed three new IL-analogous peptides: IL-K7 (Pro7-->Lys), IL-F89 (Trp8 and Trp9-->Phe), and IL-K7F89 (Pro7-->Lys; Trp8 and Trp9-->Phe). The hemolytic activity of IL-F89 is considerably lower than that of IL, whereas the antimicrobial activity of IL-K7 is greatly enhanced. In particular, the de novo peptide IL-K7F89 exhibits higher antimicrobial activity against Escherichia coli; its hemolytic activity decreased to only 10% of that of IL. Our simulated and experimental results correlated well. This approach-coupling MD simulations with experimental design-is a useful strategy toward the rational design of AMPs for potential therapeutic use.