Difference between Magainin-2 and Melittin Assemblies in Phosphatidylcholine Bilayers: Results from Coarse-Grained Simulations

Vesicle protrusion induced by antimicrobial peptides suggests common carpet mechanism for short antimicrobial peptides

Short-cationic alpha-helical antimicrobial peptides (SCHAMPs) are promising candidates to combat the growing global threat of antimicrobial resistance. They are short-sequenced, selective against bacteria, and have rapid action by destroying membranes. A full understanding of their mechanism of action will provide key information to design more potent and selective SCHAMPs. Molecular Dynamics (MD) simulations are invaluable tools that provide detailed insights into the peptide-membrane interaction at the atomic- and meso-scale level. We use atomistic and coarse-grained MD to look into the exact steps that four promising SCHAMPs-BP100, Decoralin, Neurokinin-1, and Temporin L-take when they interact with membranes. Following experimental set-ups, we explored the effects of SCHAMPs on anionic membranes and vesicles at multiple peptide concentrations. Our results showed all four peptides shared similar binding steps, initially binding to the membrane through electrostatic interactions and then flipping on their axes, dehydrating, and inserting their hydrophobic moieties into the membrane core. At higher concentrations, fully alpha-helical peptides induced membrane budding and protrusions. Our results suggest the carpet mode of action is fit for the description of SCHAMPs lysis activity and discuss the importance of large hydrophobic residues in SCHAMPs design and activity.

Anti-Herpes Simplex Virus and Anti-Inflammatory Activities of the Melittin Peptides Derived from Apis mellifera and Apis florea Venom

Herpes simplex virus (HSV) is known to cause cold sores and various diseases in humans. Importantly, HSV infection can develop latent and recurrent infections, and it is also known to cause inflammation. These infections are difficult to control, and effective treatment of the disease remains a challenge. Thus, the search for new antiviral and anti-inflammatory agents is a necessity. Melittin is a major peptide that is present in the venom of the honeybee. It possesses a number of pharmacological properties. In this study, the effects of the melittin peptides from A. mellifera (MEL-AM) and A. florea (MEL-AF) against HSV-1 and HSV-2 were evaluated at different stages during the viral multiplication cycle in an attempt to define the mode of antiviral action using plaque reduction and virucidal assays. The results revealed a new finding that melittin at 5 µg/mL demonstrated the highest inhibitory effect on HSV through the direct inactivation of viral particles, and MEL-AF displayed a greater virucidal activity. Moreover, melittin was also observed to interfere with the process of HSV attachment to the host cells. MEL-AM exhibited anti-HSV-1 and anti-HSV-2 effects with EC50 values of 4.90 ± 0.15 and 4.39 ± 0.20 µg/mL, while MEL-AF demonstrated EC50 values of 4.47 ± 0.21 and 3.95 ± 0.61 µg/mL against HSV-1 and HSV-2, respectively. However, non-cytotoxic concentrations of both types of melittin produced only slight degrees of HSV-1 and HSV-2 inhibition after viral attachment, but melittin at 5 µg/mL was able to reduce the plaque size of HSV-2 when compared to the untreated group. In addition, MEL-AM and MEL-AF also exhibited anti-inflammatory activity via the inhibition of nitric oxide production in LPS-stimulated RAW 264.7 macrophage cells, and they were also found to down-regulate the expressions of the iNOS, COX-2 and IL-6 genes. The highest inhibition of IL-6 mRNA expression was found after treatment with 10 µg/mL of MEL-AM and MEL-AF. Therefore, melittin peptides have displayed strong potential to be used as an alternative treatment for HSV infection and inflammatory diseases in the future.

pSPICA Force Field Extended for Proteins and Peptides

Many coarse-grained (CG) molecular dynamics (MD) studies have been performed to investigate biological processes involving proteins and lipids. CG force fields (FFs) in these MD studies often use implicit or nonpolar water models to reduce computational costs. CG-MD using water models cannot properly describe electrostatic screening effects owing to the hydration of ionic segments and thus cannot appropriately describe molecular events involving water channels and pores through lipid membranes. To overcome this issue, we developed a protein model in the pSPICA FF, in which a polar CG water model showing the proper dielectric response was adopted. The developed CG model greatly improved the transfer free energy profiles of charged side chain analogues across the lipid membrane. Application studies on melittin-induced membrane pores and mechanosensitive channels in lipid membranes demonstrated that CG-MDs using the pSPICA FF correctly reproduced the structure and stability of the pores and channels. Furthermore, the adsorption behavior of the highly charged nona-arginine peptides on lipid membranes changed with salt concentration, indicating the pSPICA FF is also useful for simulating protein adsorption on membrane surfaces.

Alteration of Average Thickness of Lipid Bilayer by Membrane-Deforming Inclusions

Thickness of lipid bilayer membranes is a key physical parameter determining membrane permeability and stability with respect to formation of through pores. Most membrane inclusions or impurities like amphipathic peptides, transmembrane peptides, lipid inclusions of a different molecular shape, lipid domains, and protein-lipid domains, locally deform the membrane. The detailed structure of the locally deformed region of the membrane is a kind of "fingerprint" for the inclusion type. However, most experimental methods allow determining only averaged parameters of membranes with incorporated inclusions, thus preventing the direct obtaining of the characteristics of the inclusion. Here we developed a model that allows the obtaining of characteristic parameters of three types of membrane inclusions (amphipathic peptides, transmembrane peptides, monolayer lipid patches) from experimentally observable dependencies of the average thickness of lipid bilayer on the surface concentration of the inclusions. In the case of amphipathic peptides, the model provided the peptide parameters that were in qualitative agreement with the available experimental data.

Bacterial Outer-Membrane-Mimicking Giant Unilamellar Vesicle Model for Detecting Antimicrobial Permeability

The construction of bacterial outer membrane models with native lipids like lipopolysaccharide (LPS) is a barrier to understanding antimicrobial permeability at the membrane interface. Here, we engineer bacterial outer membrane (OM)-mimicking giant unilamellar vesicles (GUVs) by constituting LPS under different pH conditions and assembled GUVs with controlled dimensions. We quantify the LPS reconstituted in GUV membranes and reveal their arrangement in the leaflets of the vesicles. Importantly, we demonstrate the applications of OM vesicles by exploring antimicrobial permeability activity across membranes. Model peptides, melittin and magainin-2, are examined where both peptides exhibit lower membrane activity in OM vesicles than vesicles devoid of LPS. Our findings reveal the mode of action of antimicrobial peptides in bacterial-membrane-mimicking models. Notably, the critical peptide concentration required to elicit activity on model membranes correlates with the cell inhibitory concentrations that revalidate our models closely mimic bacterial membranes. In conclusion, we provide an OM-mimicking model capable of quantifying antimicrobial permeability across membranes.

Organizations of melittin peptides after spontaneous penetration into cell membranes

The antimicrobial peptide, melittin, is a potential next-generation antibiotic because melittin can spontaneously form pores in bacterial cell membranes and cause cytoplasm leakage. However, the organizations of melittin peptides in cell membranes remain elusive, which impedes the understanding of the poration mechanism. In this work, we use coarse-grained and all-atom molecular dynamics (MD) simulations to investigate the organizations of melittin peptides during and after spontaneous penetration into DPPC/POPG lipid bilayers. We find that the peptides in lipid bilayers adopt either a transmembrane conformation or a U-shaped conformation, which are referred to as T- and U-peptides, respectively. Several U-peptides and/or T-peptides aggregate to form stable pores. We analyze a T-pore consisting of four T-peptides and a U-pore consisting of three U-peptides and one T-peptide. In both pores, peptides are organized in a manner such that polar residues face inward and hydrophobic residues face outward, which stabilizes the pores and produces water channels. Compared with the U-pore, the T-pore has lower energy, larger pore diameter, and higher permeability. However, the T-pore occurs less frequently than the U-pore in our simulations, probably because the formation of the T-pore is kinetically slower than the U-pore. The stability and permeability of both pores are confirmed by 300 ns all-atom MD simulations. The peptide organizations obtained in this work should deepen the understanding of the stability, poration mechanism, and permeability of melittin, and facilitate the optimization of melittin to enhance the antibacterial ability.

Cooperative antimicrobial action of melittin on lipid membranes: A coarse-grained molecular dynamics study

We conducted a series of coarse-grained molecular dynamics (CG-MD) simulations to investigate the complicated actions of melittin, which is an antimicrobial peptide (AMP) derived from honey bee venom, on a lipid membrane. To accurately simulate the AMP action, we developed and used a protein CG model as an extension of the pSPICA force field (FF), which was designed to reproduce several thermodynamic quantities and structural properties. At a low peptide-to-lipid (P/L) ratio (1/102), no defect was detected. At P/L = 1/51, toroidal pore formation was observed due to collective insertion of multiple melittin peptides from the N-termini. The pore formation was initiated by a local increase in membrane curvature in the vicinity of the peptide aggregate. At a higher P/L ratio (1/26), two more modes were detected, seemingly not controlled by the P/L ratio but by a local arrangement of melittin peptides: 1. Pore formation accompanied by lipid extraction by melittin peptides:a detergent-like mechanism. 2. A rapidly formed large pore in a significantly curved membrane: bursting. Thus, we observed three pore formation modes (toroidal pore formation, lipid extraction, and bursting) depending on the peptide concentration and local arrangement. These observations were consistent with experimental observations and hypothesized melittin modes. Through this study, we found that the local arrangements and population of melittin peptides and the area expansion rate by membrane deformation were key to the initiation of and competition among the multiple pore formation mechanisms.

Measuring Thousands of Single-Vesicle Leakage Events Reveals the Mode of Action of Antimicrobial Peptides

Host defense or antimicrobial peptides hold promise for providing new pipelines of effective antimicrobial agents. Their activity quantified against model phospholipid membranes is fundamental to a detailed understanding of their structure–activity relationships. However, classical characterization assays often lack the ability to achieve this insight. Leveraging a highly parallelized microfluidic platform for trapping and studying thousands of giant unilamellar vesicles, we conducted quantitative long-term microscopy studies to monitor the membrane-disruptive activity of archetypal antimicrobial peptides with a high spatiotemporal resolution. We described the modes of action of these peptides via measurements of the disruption of the vesicle population under the conditions of continuous peptide dosing using a range of concentrations and related the observed modes to the molecular activity mechanisms of these peptides. The study offers an effective approach for characterizing membrane-targeting antimicrobial agents in a standardized manner and for assigning specific modes of action to the corresponding antimicrobial mechanisms.

Experimental and simulation studies reveal mechanism of action of human defensin derivatives

Antimicrobial peptides (AMPs) are naturally occurring promising candidates which can be used as antibiotics against a wide variety of bacteria. The key component for using them as a potent antibiotic is that their mechanism of action is less prone to bacterial resistance. However, the molecular details of their mechanism of action is not yet fully understood. In this study, we try to shed light on the mode of action of AMPs, possible reason behind it, and their interaction with lipid bilayers through experimental as well as molecular dynamics (MD) simulation studies. The focal of our study was Human beta defensin 3 (hBD-3) which is a naturally occurring AMP. We chose three derivatives of hBD-3, namely CHRG01, KSR, and KLR for the detailed analysis presented in this study. These three peptides are evaluated for their antibacterial potency, secondary structure analysis and mechanism of action. The experimental results reveal that these peptides are active against gram positive as well as gram negative bacteria and kill bacteria by forming membrane pores. The MD simulation results correlate well with the antibacterial activity and shed light into the early membrane insertion dynamics. Moreover, the specific amino acids responsible for membrane disruptions are also identified from the MD simulations. Understanding the molecular level interaction of individual amino acids with the lipid bilayer will greatly help in the design of more efficient antimicrobial peptides.

Characteristic lengths of transmembrane peptides controlling their tilt and lateral distribution between membrane domains

Lipids and proteins of plasma membranes of eukaryotic cells are supposed to form protein-lipid domains, characterized by a different molecular order, bilayer thickness, and elastic parameters. Several mechanisms of preferable distribution of transmembrane proteins to the ordered or disordered membrane domains have been revealed. The mismatch between the length of the protein transmembrane domain and hydrophobic thickness of the lipid bilayer is considered to be an important driving force of protein lateral sorting. Utilizing the continuum theory of elasticity, we analyzed optimal configurations and preferable membrane domains for single-pass transmembrane peptides of various hydrophobic lengths and effective molecular shapes. We obtained that short transmembrane peptides stand perpendicularly to the membrane plane. The exceedance of a certain characteristic length leads to the tilt of the peptide. This length depends on the bilayer thickness. Thus, in the membrane with coexisting ordered (thicker) and disordered (thinner) phases tilting of the peptide in each phase is governed by its individual characteristic length. The lateral distribution of the peptides between ordered and disordered membrane domains is shown to be described by two additional characteristic lengths. The exceedance of the smaller one drives the peptide towards a more ordered and thicker membrane, while the exceedance of the larger characteristic length switches the preferable membrane domain from ordered and thicker to disordered and thinner. Thus, membrane proteins with long enough transmembrane domains are predicted to accumulate in the thinner disordered membrane as compared to the thicker ordered bilayer. For hourglass-like and barrel-like shaped transmembrane peptides the specific regime of sorting was obtained: the peptides distributed almost equally between the phases in a wide range of peptide lengths. This finding allowed explaining the experimental data on lateral distribution of transmembrane peptide tLAT.

Vectorial insertion of a β-helical peptide into membrane: a theoretical study on polytheonamide B

Spontaneous unidirectional, or vectorial, insertion of transmembrane peptides is a fundamental biophysical process for toxin and viral actions. Polytheonamide B (pTB) is a potent cytotoxic peptide with a β^6.3-helical structure. Previous experimental studies revealed that the pTB inserts into the membrane in a vectorial fashion and forms a channel with its single molecular length long enough to span the membrane. Also, molecular dynamics simulation studies demonstrated that the pTB is prefolded in aqueous solution. These are unique features of pTB because most of the peptide toxins form channels through oligomerization of transmembrane helices. Here, we performed all-atom molecular dynamics simulations to examine the dynamic mechanism of the vectorial insertion of pTB, providing underlying elementary processes of the membrane insertion of a prefolded single transmembrane peptide. We find that the insertion of pTB proceeds with only the local lateral compression of the membrane in three successive phases: "landing," "penetration," and "equilibration" phases. The free energy calculations using the replica-exchange umbrella sampling simulations present an energy cost of 4.3 kcal/mol at the membrane surface for the membrane insertion of pTB from bulk water. The trajectories of membrane insertion revealed that the insertion process can occur in two possible pathways, namely "trapped" and "untrapped" insertions; in some cases, pTB is trapped in the upper leaflet during the penetration phase. Our simulations demonstrated the importance of membrane anchoring by the hydrophobic N-terminal blocking group in the landing phase, leading to subsequent vectorial insertion.

Mechanistic Understanding from Molecular Dynamics in Pharmaceutical Research 2: Lipid Membrane in Drug Design

We review the use of molecular dynamics (MD) simulation as a drug design tool in the context of the role that the lipid membrane can play in drug action, i.e., the interaction between candidate drug molecules and lipid membranes. In the standard "lock and key" paradigm, only the interaction between the drug and a specific active site of a specific protein is considered; the environment in which the drug acts is, from a biophysical perspective, far more complex than this. The possible mechanisms though which a drug can be designed to tinker with physiological processes are significantly broader than merely fitting to a single active site of a single protein. In this paper, we focus on the role of the lipid membrane, arguably the most important element outside the proteins themselves, as a case study. We discuss work that has been carried out, using MD simulation, concerning the transfection of drugs through membranes that act as biological barriers in the path of the drugs, the behavior of drug molecules within membranes, how their collective behavior can affect the structure and properties of the membrane and, finally, the role lipid membranes, to which the vast majority of drug target proteins are associated, can play in mediating the interaction between drug and target protein. This review paper is the second in a two-part series covering MD simulation as a tool in pharmaceutical research; both are designed as pedagogical review papers aimed at both pharmaceutical scientists interested in exploring how the tool of MD simulation can be applied to their research and computational scientists interested in exploring the possibility of a pharmaceutical context for their research.

Ligand-decoration determines the translational and rotational dynamics of nanoparticles on a lipid bilayer membrane

Nanoparticles (NPs) promise a huge potential for clinical diagnostic and therapeutic applications. However, nano-bio (e.g., the NP-cell membrane) interactions and underlying mechanisms are still largely elusive. In this study, two types of congeneric peptides, namely PGLa and magainin 2 (MAG2), with similar membrane activities were employed as model ligands for NP decoration, and the diffusion behaviours (including both translation and rotation) of the ligand-decorated NPs on a lipid bilayer membrane were studied via molecular dynamics simulations. It was found that, although both PGLa- and MAG2-coated NPs showed alternatively "hopping" and "jiggling" diffusions, the PGLa-coated ones had an enhanced circling at the hopping stage, while a much confined circling at the jiggling stage. In contrast, the MAG2-coated NPs demonstrated constant circling tendencies throughout the diffusion process. Such differences in the coupling between translational and rotational dynamics of these two types of NPs are ascribed to the different ligand-lipid interactions of PGLa and MAG2, in which the PGLa ligands prefer to vertically insert into the membrane, while MAG2 tends to lie flat on the membrane surface. Our results are helpful for the understanding the underlying associations between the NP motions and their interfacial membrane interactions, and shed light on the possibility of regulating NP behaviours on a cellular surface for better biomedical uses.

Opposite Regulatory Effects of Immobilized Cations on the Folding Vs. Assembly of Melittin

Ions are crucial in modulating the protein structure. For the free ions in bulk solution, ammonium is kosmotropic (structure forming) and guanidinium is chaotropic (structure breaking) to the protein structure within the Hofmeister series. However, the effect of immobilized ions on a protein surface is less explored. Herein, we explored the influence of two immobilized cations (ammonium in the side chain of lysine and guanidinium in the side chain of arginine) on the folding and assembly of melittin. Melittin adopts an α-helix structure and is driven by hydrophobic interactions to associate into a helical bundle. To test the influence of immobilized cations on the peptide structure, we designed the homozygous mutants exclusively containing ammonium (melittin-K) or guanidinium (melittin-R) and compared the differences of melittin-K vs. melittin-R in their folding, assembly, and molecular functions. The side chains of lysine and arginine differ in their influences on the folding and assembly of melittin. Specifically, the side chain of R increases the α-helical propensity of melittin relative to that of K, following an inverse Hofmeister series. In contrast, the side chain of K favors the assembly of melittin relative to the side chain of R in line with a direct Hofmeister series. The opposite regulatory effects of immobilized cations on the folding and assembly of melittin highlight the complexity of the noncovalent interactions that govern protein intermolecular architecture.

Advances in Molecular Understanding of α-Helical Membrane-Active Peptides

Biological membranes separate the interior of cells or cellular compartments from their outer environments. This barrier function of membranes can be disrupted by membrane-active peptides, some of which can spontaneously penetrate through the membranes or open leaky transmembrane pores. However, the origin of their activity/toxicity is not sufficiently understood for the development of more potent peptides. To this day, there are no design rules that would be generally valid, and the role of individual amino acids tends to be sequence-specific.In this Account, we describe recent progress in understanding the design principles that govern the activity of membrane-active peptides. We focus on α-helical amphiphilic peptides and their ability to (1) translocate across phospholipid bilayers, (2) form transmembrane pores, or (3) act synergistically, i.e., to produce a significantly more potent effect in a mixture than the individual components.We refined the description of peptide translocation using computer simulations and demonstrated the effect of selected residues. Our simulations showed the necessity to explicitly include charged residues in the translocation description to correctly sample the membrane perturbations they can cause. Using this description, we calculated the translocation of helical peptides with and without the kink induced by the proline/glycine residue. The presence of the kink had no effect on the translocation barrier, but it decreased the peptide affinity to the membrane and reduced the peptide stability inside the membrane. Interestingly, the effects were mainly caused by the peptide's increased polarity, not the higher flexibility of the kink.Flexibility plays a crucial role in pore formation and affects distinct pore structures in different ways. The presence of a kink destabilizes barrel-stave pores, because the kink prevents the tight packing of peptides in the bundle, which is characteristic of the barrel-stave structure. In contrast, the kink facilitates the formation of toroidal pores, where the peptides are only loosely arranged and do not need to closely assemble. The exact position of the kink in the sequence further determines the preferred arrangement of peptides in the pore, i.e., an hourglass or U-shaped structure. In addition, we demonstrated that two self-associated (via termini) helical peptides could mimic the behavior of peptides with a helix-kink-helix motif.Finally, we review the recent findings on the peptide synergism of the archetypal mixture of Magainin 2 and PGLa peptides. We focused on a bacterial plasma membrane mimic that contains negatively charged lipids and lipids with negative intrinsic curvature. We showed that the synergistic action of peptides was highly dependent on the lipid composition. When the lipid composition and peptide/lipid ratios were changed, the systems exhibited more complex behavior than just the previously reported pore formation. We observed membrane adhesion, fusion, and even the formation of the sponge phase in this regime. Furthermore, enhanced adhesion/partitioning to the membrane was reported to be caused by lipid-induced peptide aggregation.In conclusion, the provided molecular insight into the complex behavior of membrane-active peptides provides clues for the design and modification of antimicrobial peptides or toxins.

Molecular Dynamics Simulation of the Interaction of Two Linear Battacin Analogs with Model Gram-Positive and Gram-Negative Bacterial Cell Membranes

Antimicrobial peptides (AMPs) are a potential solution to the increasing threat of antibiotic resistance, but successful design of active but nontoxic AMPs requires understanding their mechanism of action. Molecular dynamics (MD) simulations can provide atomic-level information regarding how AMPs interact with the cell membrane. Here, we have used MD simulations to study two linear analogs of battacin, a naturally occurring cyclic, lipidated, nonribosomal AMP. Like battacin, these analogs are active against Gram-negative multidrug resistant and Gram-positive bacteria, but they are less toxic than battacin. Our simulations show that this activity depends upon a combination of positively charged and hydrophobic moieties. Favorable interactions with negatively charged membrane lipid head groups drive association with the membrane and insertion of hydrophobic residues, and the N-terminal lipid anchors the peptides to the membrane surface. Both effects are required for stable membrane binding.

Differential Stability of Aurein 1.2 Pores in Model Membranes of Two Probiotic Strains

Aurein 1.2 is an antimicrobial peptide from the skin secretion of an Australian frog. In the previous experimental work, we reported a differential action of aurein 1.2 on two probiotic strains Lactobacillus delbrueckii subsp. bulgaricus (CIDCA 331) and Lactobacillus delbrueckii subsp. lactis (CIDCA 133). The differences found were attributed to the bilayer compositions. Cell cultures and CIDCA 331-derived liposomes showed higher susceptibility than the ones derived from the CIDCA 133 strain, leading to content leakage and structural disruption. Here, we used molecular dynamics simulations to explore these systems at the atomistic level. We hypothesize that if the antimicrobial peptides organized themselves to form a pore, it will be more stable in membranes that emulate the CIDCA 331 strain than in those of the CIDCA 133 strain. To test this hypothesis, we simulated preassembled aurein 1.2 pores embedded into bilayer models that emulate the two probiotic strains. It was found that the general behavior of the systems depends on the composition of the membrane rather than the preassemble system characteristics. Overall, it was observed that aurein 1.2 pores are more stable in the CIDCA 331 model membranes. This fact coincides with the high susceptibility of this strain against antimicrobial peptide. In contrast, in the case of the CIDCA 133 model membranes, peptides migrate to the water-lipid interphase, the pore shrinks, and the transport of water through the pore is reduced. The tendency of glycolipids to make hydrogen bonds with peptides destabilizes the pore structures. This feature is observed to a lesser extent in CIDCA 331 due to the presence of anionic lipids. Glycolipid transverse diffusion (flip-flop) between monolayers occurs in the pore surface region in all the cases considered. These findings expand our understanding of the antimicrobial peptide resistance properties of probiotic strains.

Individual Roles of Peptides PGLa and Magainin 2 in Synergistic Membrane Poration

Synergy between antimicrobial peptides PGLa and Magainin 2 (MAG2) provides an efficient way to enhance their antimicrobial ability. However, the underlying molecular mechanism of such synergy, especially the individual roles of each peptide, remains poorly understood. We combined a giant unilamellar vesicle leakage assay, in situ interfacial photovoltage testing, and molecular dynamics to investigate membrane poration under the action of PGLa, MAG2, or a PGLa/MAG2 mixture. Our results clearly show the different membrane action modes of the three systems and demonstrate the importance of forming PGLa-MAG2 heterodimers in the membrane poration process. PGLa inserted into and extracted from a membrane rapidly and continually with minimal aggregation and produced only transient, small pores. In contrast, MAG2 peptides tended to aggregate together on the membrane surface or only shallowly embed in the membrane. Additionally, the PGLa and MAG2 residues were well integrated into the membrane via the formation of PGLa-MAG2 heterodimers. The membrane defect produced by the rapid insertion of PGLa was stabilized by MAG2, which further recruited other peptides for the formation of PGLa-MAG2 heterodimers and even heterodimer clusters. Growth in pore size then occurred in a step-by-step process involving the formation and assembly of heterodimer clusters within the membrane. Our results provide insight into the complicated synergy that occurs between PGLa and MAG2 during membrane poration and will assist in the design of new antimicrobial peptides.

Correlation between Single-Molecule Dynamics and Biological Functions of Antimicrobial Peptide Melittin

Many fundamental biological processes occur on cell membranes, and a typical example is the membrane permeabilization by peptides for an antimicrobial purpose. Previous studies of the underlying mechanism mostly focus on structural changes of membranes and peptides during their interactions. Herein, from a new perspective of single-molecule dynamics, the real-time three-dimensional motions of individual phospholipid and peptide molecules were monitored, and specifically, their correlation with the membrane poration function of melittin, a most representative natural antimicrobial peptide, was studied. We found that the adsorption and accumulation of melittin on the membrane surface significantly sped up the lateral diffusion of lipids surrounding the peptides, which in turn facilitated the peptide insertion at such heterogeneous regions. A unique "U"-bending pathway of melittin during membrane insertion and the ultimate formation of toroidal pores with dynamical translocations of peptides and lipids with several metastable states between the two leaflets of bilayer were observed.

Continuum Models of Membrane Fusion: Evolution of the Theory

Starting from fertilization, through tissue growth, hormone secretion, synaptic transmission, and sometimes morbid events of carcinogenesis and viral infections, membrane fusion regulates the whole life of high organisms. Despite that, a lot of fusion processes still lack well-established models and even a list of main actors. A merger of membranes requires their topological rearrangements controlled by elastic properties of a lipid bilayer. That is why continuum models based on theories of membrane elasticity are actively applied for the construction of physical models of membrane fusion. Started from the view on the membrane as a structureless film with postulated geometry of fusion intermediates, they developed along with experimental and computational techniques to a powerful tool for prediction of the whole process with molecular accuracy. In the present review, focusing on fusion processes occurring in eukaryotic cells, we scrutinize the history of these models, their evolution and complication, as well as open questions and remaining theoretical problems. We show that modern approaches in this field allow continuum models of membrane fusion to stand shoulder to shoulder with molecular dynamics simulations, and provide the deepest understanding of this process in multiple biological systems.

Lipid-specific interactions determine the organization and dynamics of membrane-active peptide melittin

The cell membranes of different cells deviate significantly in lipid compositions and thus provide varying biological environments to modulate the diffusion, organization and the resultant function of biomacromolecules. However, the detailed modulation mechanism remains elusive especially in consideration of the current overuse of the simplified membrane models such as the pure phosphatidylcholine (PC) membrane. In this work, with the typical membrane-active peptide melittin, we demonstrated that a more complicated membrane environment, such as the bacterial (IME) or plasma membrane (PM), would significantly change the organization and dynamics of melittin, by using molecular dynamics simulations as a "computational microscope". It was found that in these membrane systems, adding melittin would cause a varying degree of reduction in the lateral diffusion of lipids due to the different assembly states of peptides. Melittin tended to aggregate to oligomers in the pure PC membrane, mostly as a tetramer or trimer, while in IME or PM, its degree of oligomerization was significantly reduced. More surprisingly, melittin displayed a strong affinity with ganglioside GM3 in PM, leading to the formation of melittin-GM3 nanoclusters, which hindered its diffusion and further oligomerization. Additionally, small changes in the residue sequence of melittin could modulate the degree or structure of the peptide oligomer. Our work provides a typical example of a study on the organization and dynamics of pore-forming peptides in specific membrane environments and has great significance on the optimization of peptide sequences and the design of helix bundles in the membrane for target biological function.

Effect of helical kink in antimicrobial peptides on membrane pore formation

Every cell is protected by a semipermeable membrane. Peptides with the right properties, for example Antimicrobial peptides (AMPs), can disrupt this protective barrier by formation of leaky pores. Unfortunately, matching peptide properties with their ability to selectively form pores in bacterial membranes remains elusive. In particular, the proline/glycine kink in helical peptides was reported to both increase and decrease antimicrobial activity. We used computer simulations and fluorescence experiments to show that a kink in helices affects the formation of membrane pores by stabilizing toroidal pores but disrupting barrel-stave pores. The position of the proline/glycine kink in the sequence further controls the specific structure of toroidal pore. Moreover, we demonstrate that two helical peptides can form a kink-like connection with similar behavior as one long helical peptide with a kink. The provided molecular-level insight can be utilized for design and modification of pore-forming antibacterial peptides or toxins.

Molecular dynamics study of membrane permeabilization by wild-type and mutant lytic peptides from the non-enveloped Flock House virus

Flock House virus (FHV) serves as a model system for understanding infection mechanisms utilized by non-enveloped viruses to transport across cellular membranes. During the infection cycle of FHV, a fundamental stage involves disruption of the endosomal membrane by membrane active peptides, following externalization of the peptides from the capsid interior. The FHV lytic agents are the 44 C-terminal amino acids residues of the capsid protein, which are auto-catalytically cleaved during the capsid maturation process. The cleaved peptides are termed γ peptides. In this study, we perform multi-scale molecular dynamics simulations including 40 μs all-atom molecular dynamics simulations to study the behavior of pre-inserted transmembrane lytic peptides at a high concentration in a neutral membrane. We study the dynamical organization among peptides to form oligomeric bundles in four systems including the wild-type γ peptide and three mutant forms; namely, a truncation mutant in which the 23 C-terminal residues are deleted (γ₁), a construct where the 8 C-terminal residues of γ are fused to γ₁ (Δ385-399 γ) and a single-point mutant (F402A γ), all of which have been experimentally shown to drastically affect infectivity and lytic activity compared to the wild-type γ. Our results shed light on the actions of varied forms of the FHV lytic peptide including membrane insertion, trans-membrane stability, peptide oligomerization, water permeation activity and dynamic pore formation. Findings from this study provide detailed structural information and rationale for the differences in lytic activity among variants of FHV γ.

How Melittin Inserts into Cell Membrane: Conformational Changes, Inter-Peptide Cooperation, and Disturbance on the Membrane

Antimicrobial peptides (AMPs), as a key component of the immune defense systems of organisms, are a promising solution to the serious threat of drug-resistant bacteria to public health. As one of the most representative and extensively studied AMPs, melittin has exceptional broad-spectrum activities against microorganisms, including both Gram-positive and Gram-negative bacteria. Unfortunately, the action mechanism of melittin with bacterial membranes, especially the underlying physics of peptide-induced membrane poration behaviors, is still poorly understood, which hampers efforts to develop melittin-based drugs or agents for clinical applications. In this mini-review, we focus on recent advances with respect to the membrane insertion behavior of melittin mostly from a computational aspect. Membrane insertion is a prerequisite and key step for forming transmembrane pores and bacterial killing by melittin, whose occurrence is based on overcoming a high free-energy barrier during the transition of melittin molecules from a membrane surface-binding state to a transmembrane-inserting state. Here, intriguing simulation results on such transition are highlighted from both kinetic and thermodynamic aspects. The conformational changes and inter-peptide cooperation of melittin molecules, as well as melittin-induced disturbances to membrane structure, such as deformation and lipid extraction, are regarded as key factors influencing the insertion of peptides into membranes. The associated intermediate states in peptide conformations, lipid arrangements, membrane structure, and mechanical properties during this process are specifically discussed. Finally, potential strategies for enhancing the poration ability and improving the antimicrobial performance of AMPs are included as well.

Free energy analysis of membrane pore formation process in the presence of multiple melittin peptides

Understanding the molecular mechanism underlying pore formation in lipid membranes by antimicrobial peptides is of great importance in biological sciences as well as in drug design applications. Melittin has been widely studied as a pore forming peptide, though the molecular mechanism for pore formation is still illusive. We examined the free energy barrier for the creation of a pore in lipid membranes with and without multiple melittin peptides. It was found that six melittin peptides significantly stabilized a pore, though a small barrier (a few k_BT) for the formation still existed. With five melittin peptides or fewer, the pore formation barrier was much higher, though the established pore was in a local energy minimum. Although seven melittins effectively reduced the free energy barrier, a single melittin peptide left the pore after a long time MD simulation probably because of the overcrowded environment around the bilayer pore. Thus, it is highly selective for the number of melittin peptides to stabilize the membrane pore, as was also suggested by the line tension evaluations. The free energy cost required to insert a single melittin into the membrane is too high to explain the one-by-one insertion mechanism for pore formation, which also supports the collective melittin mechanism for pore formation.

Mechanistic Landscape of Membrane-Permeabilizing Peptides

Membrane permeabilizing peptides (MPPs) are as ubiquitous as the lipid bilayer membranes they act upon. Produced by all forms of life, most membrane permeabilizing peptides are used offensively or defensively against the membranes of other organisms. Just as nature has found many uses for them, translational scientists have worked for decades to design or optimize membrane permeabilizing peptides for applications in the laboratory and in the clinic ranging from antibacterial and antiviral therapy and prophylaxis to anticancer therapeutics and drug delivery. Here, we review the field of membrane permeabilizing peptides. We discuss the diversity of their sources and structures, the systems and methods used to measure their activities, and the behaviors that are observed. We discuss the fact that "mechanism" is not a discrete or a static entity for an MPP but rather the result of a heterogeneous and dynamic ensemble of structural states that vary in response to many different experimental conditions. This has led to an almost complete lack of discrete three-dimensional active structures among the thousands of known MPPs and a lack of useful or predictive sequence-structure-function relationship rules. Ultimately, we discuss how it may be more useful to think of membrane permeabilizing peptides mechanisms as broad regions of a mechanistic landscape rather than discrete molecular processes.

Potential of mean force for insertion of antimicrobial peptide melittin into a pore in mixed DOPC/DOPG lipid bilayer by molecular dynamics simulation

Antimicrobial peptides (AMPs) inactivate microorganisms by forming transmembrane pores in a cell membrane through adsorption and aggregation. Energetics of addition of an AMP to a transmembrane pore is important for evaluation of its formation and growth. Such information is essential for the characterization of pore forming ability of peptides in cell membranes. This study quantifies the potential of mean force through molecular dynamics (MD) simulation for the addition of melittin, a naturally occurring AMP, into a DOPC/DOPG mixed bilayer, a mimic of bacterial membrane, for different extents of insertion into either a bilayer or a pore consisting of three to six transmembrane peptides. The energy barrier for insertion of a melittin molecule into the bilayer was highest in the absence of transmembrane peptides and decreased for the number of transmembrane peptides from three to six, eventually approaching zero. The decrease in free energy for complete insertion of peptide was found to be higher for larger pore size. Water channel formation occurred only for insertion into pores consisting of three or more transmembrane peptides with the radius of water channel being larger for a larger number of transmembrane peptides. The structure of the pore was found to be paraboloid. The estimated free energy barrier for insertion of melittin into an ideal paraboloid pore accounting for different intermolecular interactions was consistent with MD simulation results. The results reported in this manuscript will be useful for the development of a model for nucleation of pores and a rational methodology for selection of synthetic antimicrobial peptides.

Penetration of antimicrobial peptides in a lung surfactant model

Molecular dynamics simulations were successfully performed to understand the absorption mechanism of antimicrobial peptides LL-37, CATH-2, and SMAP-29 in a lung surfactant model. The antimicrobial peptides quickly penetrate in the lung surfactant model in dozens or hundreds nanoseconds, but they electrostatically interact with the lipid polar heads during the simulation time of 2 μs. This electrostatic interaction should be the explanation for the inactivation of the antimicrobial peptides when co-administrated with lung surfactant. As they strongly interact with the lipid polar heads of the lung surfactant, there is no positive charge available on the antimicrobial peptide to attack the negatively charged bacteria membrane. In order to avoid the interaction of peptides with the lipid polar heads, sodium cholate was used to form nanoparticles which act as an absorption enhancer of all antimicrobial peptides used in this investigation. The nanoparticles of 150 molecules of sodium cholate with one peptide were inserted on the top of the lung surfactant model. The nanoparticles penetrated into the lung surfactant model, spreading the sodium cholate molecules around the lipid polar heads. The sodium cholate molecules seem to protect the peptides from the interaction with the lipid polar heads, leaving them free to be delivered to the water phase. The penetration of peptides alone or even the peptide nanoparticles with sodium cholate do not collapse the lung surfactant model, indicating to be a promisor drug delivery system to the lung. The implications of this finding are that antimicrobial peptides may only be co-administered with an absorption enhancer such as sodium cholate into lung surfactant in order to avoid inactivation of their antimicrobial activity.

Antimicrobial action of the cationic peptide, chrysophsin-3: a coarse-grained molecular dynamics study

Antimicrobial peptides (AMPs) are small cationic proteins that are able to destabilize a lipid bilayer structure through one or more modes of action. In this study, we investigate the processes of peptide aggregation and pore formation in lipid bilayers and vesicles by the highly cationic AMP, Chrysophsin-3 (chrys-3), using coarse-grained molecular dynamics (CG-MD) simulations and potential of mean force calculations. We study long 50 μs simulations of chrys-3 at different concentrations, both at the surface of dipalmitoylphosphatidylcholine (DPPC) and palmitoyloleoylphosphatidylcholine (POPC) bilayers, and also interacting within the interior of the lipid membrane. We show that aggregation of peptides at the surface, leads to pronounced deformation of lipid bilayers, leading in turn to lipid protrusions for peptide : ligand ratios > 1 : 12. In addition, aggregation of chrys-3 peptides within the centre of a lipid bilayer leads to spontaneous formation of pores and aggregates. Both mechanisms of interaction are consistent with previously reported experimental data for chrys-3. Similar results are observed also in POPC vesicles and mixed lipid bilayers composed of the zwitterionic lipid palmitoyloleoylphosphatidylethanolamine (POPE) and the negatively charged lipid palmitoyloleoylphosphatidylglycerol (POPG). The latter are employed as models of the bacterial membrane of Escherichia coli.

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'.

Differential Interaction of Antimicrobial Peptides with Lipid Structures Studied by Coarse-Grained Molecular Dynamics Simulations

In this work; we investigated the differential interaction of amphiphilic antimicrobial peptides with 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) lipid structures by means of extensive molecular dynamics simulations. By using a coarse-grained (CG) model within the MARTINI force field; we simulated the peptide-lipid system from three different initial configurations: (a) peptides in water in the presence of a pre-equilibrated lipid bilayer; (b) peptides inside the hydrophobic core of the membrane; and (c) random configurations that allow self-assembled molecular structures. This last approach allowed us to sample the structural space of the systems and consider cooperative effects. The peptides used in our simulations are aurein 1.2 and maculatin 1.1; two well-known antimicrobial peptides from the Australian tree frogs; and molecules that present different membrane-perturbing behaviors. Our results showed differential behaviors for each type of peptide seen in a different organization that could guide a molecular interpretation of the experimental data. While both peptides are capable of forming membrane aggregates; the aurein 1.2 ones have a pore-like structure and exhibit a higher level of organization than those conformed by maculatin 1.1. Furthermore; maculatin 1.1 has a strong tendency to form clusters and induce curvature at low peptide-lipid ratios. The exploration of the possible lipid-peptide structures; as the one carried out here; could be a good tool for recognizing specific configurations that should be further studied with more sophisticated methodologies.

Effect of physicochemical properties of peptides from soy protein on their antimicrobial activity

Antimicrobial peptides (AMPs) kill microbial cells through insertion and damage/permeabilization of the cytoplasmic cell membranes and has applications in food safety and antibiotic replacement. Soy protein is an attractive, abundant natural source for commercial production of AMPs. In this research, explicit solvent molecular dynamics (MD) simulation was employed to investigate the effects of (i) number of total and net charges, (ii) hydrophobicity (iii) hydrophobic moment and (iv) helicity of peptides from soy protein on their ability to bind to lipid bilayer and their transmembrane aggregates to form pores. Interaction of possible AMP segments from soy protein with 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine/1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoglycerol (POPC/POPG) bilayers, a mimic of bacterial cell membrane, was investigated. Pore formation was insensitive to helicity and occurred for hydrophobicity threshold in the range of -0.3-0kcal/mol, hydrophobic moment threshold of 0.3kcal/mol, net charge threshold of 2. Though low hydrophobicity and high number of charges help in the formation of water channel for transmembrane aggregates, insertion of peptides with these properties requires overcome of energy barrier, as shown by potential of mean force calculations, thereby resulting in low antimicrobial activity. Experimental evaluation of antimicrobial activity of these peptides against Gram positive L. monocytogenes and Gram negative E. coli as obtained by spot-on-lawn assay was consistent with simulation results. These results should help in the development of guidelines for selection of peptides with antimicrobial activity based on their physicochemical properties.

Protein Aggregation and Molecular Crowding: Perspectives From Multiscale Simulations

Cells are extremely crowded environments, thus the use of diluted salted aqueous solutions containing a single protein is too simplistic to mimic the real situation. Macromolecular crowding might affect protein structure, folding, shape, conformational stability, binding of small molecules, enzymatic activity, interactions with cognate biomolecules, and pathological aggregation. The latter phenomenon typically leads to the formation of amyloid fibrils that are linked to several lethal neurodegenerative diseases, but that can also play a functional role in certain organisms. The majority of molecular simulations performed before the last few years were conducted in diluted solutions and were restricted both in the timescales and in the system dimensions by the available computational resources. In recent years, several computational solutions were developed to get close to physiological conditions. In this review we summarize the main computational techniques used to tackle the issue of protein aggregation both in a diluted and in a crowded environment.

Coarse-grained simulations of hemolytic peptide δ-lysin interacting with a POPC bilayer

δ-lysin, secreted by a Gram-positive bacterium Staphylococcus aureus, is a 26-residue membrane active peptide that shares many common features with antimicrobial peptides (AMPs). However, it possesses a few unique features that differentiate itself from typical AMPs. In particular, δ-lysin has zero net charge, even though it has many charged residues, and it preferentially lyses eukaryotic cells over bacterial cells. Here, we present the results of coarse-grained molecular dynamics simulations of δ-lysin interacting with a zwitterionic membrane over a wide range of peptide concentrations. When the peptides concentration is low, spontaneous dimerization of peptides is observed on the membrane surface, but deep insertion of peptides or pore formation was not observed. However, the calculated free energy of peptide insertion suggests that a small fraction of peptides is likely to be present inside the membrane at the peptide concentrations typically seen in dye efflux experiments. When the simulations with multiple peptides are carried out with a single pre-inserted transmembrane peptide, spontaneous pore formation occurs with a peptide-to-lipid ratio (P/L) as low as P/L=1:42. Inter-peptide salt bridges among the transmembrane peptides seem to play a role in creating compact pores with very low level of hydration. More importantly, the transmembrane peptides making up the pore are constantly pushed to the opposite side of the membrane when the mass imbalance between the two sides of membrane is significant. Thus, the pore is very dynamic, allowing multiple peptides to translocate across the membrane simultaneously.

Interaction of peptides with cell membranes: insights from molecular modeling

The investigation of the interaction of peptides with cell membranes is the focus of active research. It can enhance the understanding of basic membrane functions such as membrane transport, fusion, and signaling processes, and it may shed light on potential applications of peptides in biomedicine. In this review, we will present current advances in computational studies on the interaction of different types of peptides with the cell membrane. Depending on the properties of the peptide, membrane, and external environment, the peptide-membrane interaction shows a variety of different forms. Here, on the basis of recent computational progress, we will discuss how different peptides could initiate membrane pores, translocate across the membrane, induce membrane endocytosis, produce membrane curvature, form fibrils on the membrane surface, as well as interact with functional membrane proteins. Finally, we will present a conclusion summarizing recent progress and providing some specific insights into future developments in this field.

Molecular dynamics simulations of biological membranes and membrane proteins using enhanced conformational sampling algorithms

This paper reviews various enhanced conformational sampling methods and explicit/implicit solvent/membrane models, as well as their recent applications to the exploration of the structure and dynamics of membranes and membrane proteins. Molecular dynamics simulations have become an essential tool to investigate biological problems, and their success relies on proper molecular models together with efficient conformational sampling methods. The implicit representation of solvent/membrane environments is reasonable approximation to the explicit all-atom models, considering the balance between computational cost and simulation accuracy. Implicit models can be easily combined with replica-exchange molecular dynamics methods to explore a wider conformational space of a protein. Other molecular models and enhanced conformational sampling methods are also briefly discussed. As application examples, we introduce recent simulation studies of glycophorin A, phospholamban, amyloid precursor protein, and mixed lipid bilayers and discuss the accuracy and efficiency of each simulation model and method. This article is part of a Special Issue entitled: Membrane Proteins edited by J.C. Gumbart and Sergei Noskov.

Pore Structure and Synergy in Antimicrobial Peptides of the Magainin Family

Magainin 2 and PGLa are among the best-studied cationic antimicrobial peptides. They bind preferentially to negatively charged membranes and apparently cause their disruption by the formation of transmembrane pores, whose detailed structure is still unclear. Here we report the results of 5–9 μs all-atom molecular dynamics simulations starting from tetrameric transmembrane helical bundles of these two peptides, as well as their stoichiometric mixture, and the analog MG-H2 in DMPC or 3:1 DMPC/DMPG membranes. The simulations produce pore structures that appear converged, although some effect of the starting peptide arrangement (parallel vs. antiparallel) is still observed on this timescale. The peptides remain mostly helical and adopt tilted orientations. The calculated tilt angles for PGLa are in excellent agreement with recent solid state NMR experiments. The antiparallel dimer structure in the magainin 2 simulations resembles previously determined NMR and crystal structures. More transmembrane orientations and a larger and more ordered pore are seen in the 1:1 heterotetramer with an antiparallel helix arrangement. Insights into the mechanism of synergy between these two peptides are obtained via implicit solvent modeling of homo- and heterodimers and analysis of interactions in the atomistic simulations. This analysis suggests stronger pairwise interactions in the heterodimer than in the two homodimers.

The emergence of antibiotic resistance has created a compelling need for new potent antibiotics. Antimicrobial peptides, naturally produced by many organisms, have long been pursued as a means to fill this gap. However, they have not yet realized their practical potential, partly due to a lack of full understanding of their mechanism of action. The formation of pores by these peptides in bacterial membranes has been demonstrated by multiple biophysical approaches, but the detailed structure of these pores is still unclear. The magainin family of antimicrobial peptides have been thoroughly studied over the last three decades. We have generated atomic-resolution models of membrane pores generated by tetrameric assemblies of different magainin peptides by means of molecular dynamics simulations on the multimicrosecond timescale. The results complement previous experimental findings and set the stage for further development of theoretical and experimental approaches that may ultimately allow the rational design of improved antibiotics.

Opening of the blood-brain barrier tight junction due to shock wave induced bubble collapse: a molecular dynamics simulation study

Passage of a shock wave across living organisms may produce bubbles in the blood vessels and capillaries. It was suggested that collapse of these bubbles imposed by an impinging shock wave can be responsible for the damage or even destruction of the blood-brain barrier. To check this possibility, we performed molecular dynamics computer simulations on systems that contained a model of tight junction from the blood-brain barrier. In our model, we represent the tight junction by two pairs of interacting proteins, claudin-15. Some of the simulations were done in the absence of a nanobubble, some in its presence. Our simulations show that when no bubble is present in the system, no damage to tight junction is observed when the shock wave propagates across it. In the presence of a nanobubble, even when the impulse of the shock wave is relatively low, the implosion of the bubble causes serious damage to our model tight junction.

SLIM: an improved generalized Born implicit membrane model

In most implicit continuum models, membranes are represented as heterogeneous dielectric environments, but their treatment within computationally efficient generalized Born (GB) models is challenging. Despite several previous attempts, an adequate description of multiple dielectric regions in implicit GB-based membrane models that reproduce the qualitative and quantitative features of Poisson-Boltzmann (PB) electrostatics remains an unmet prerequisite of qualitatively correct implicit membrane models. A novel scheme (SLIM) to decompose one environment consisting of multiple dielectric regions into a sum of multiple environments consisting only of two dielectric regions each is proposed to solve this issue. These simpler environments can be treated with established GB methods. This approach captures qualitative features of PB electrostatic that are not present in previous models. Simulations of three membrane proteins demonstrate that this model correctly reproduces known properties of these proteins in agreement with experimental or other computational studies.

Amphipathic membrane-active peptides recognize and stabilize ruptured membrane pores: exploring cause and effect with coarse-grained simulations

Induction of membrane pores has been suggested as the common molecular action by which a variety of amphipathic membrane-active peptides cause damage to cells. In this study, we have performed coarse-grained molecular dynamics simulations to establish two clear molecular processes that seem critical for the activity of amphipathic peptides. They are (i) the recognition and (ii) the stabilization of ruptured membrane pores. By considering 12 structurally different peptide types, we reveal that peptide secondary structure content, hydrophobicity, and length are important physicochemical factors that allow amphipathic peptides to aggregate in and stabilize ruptured membrane pores. The simulated inner diameters of peptide-stabilized membrane pores are in good agreement with available experimental data. However, the orientations of α-helical peptides in the membrane pore were found to be quite dispersed. This supports recent challenges to the traditional depictions to peptide orientations in the classical toroidal and barrel-stave pore models.

Experimental and computational investigation of the effect of hydrophobicity on aggregation and osteoinductive potential of BMP-2-derived peptide in a hydrogel matrix

An attractive approach to reduce the undesired side effects of bone morphogenetic proteins (BMPs) in regenerative medicine is to use osteoinductive peptide sequences derived from BMPs. Although the structure and function of BMPs have been studied extensively, there is limited data on structure and activity of BMP-derived peptides immobilized in hydrogels. The objective of this work was to investigate the effect of concentration and hydrophobicity of the BMP-2 peptide, corresponding to residues 73-92 of the knuckle epitope of BMP-2 protein, on peptide aggregation and osteogenic differentiation of human mesenchymal stem cells encapsulated in a polyethylene glycol (PEG) hydrogel. The peptide hydrophobicity was varied by capping PEG chain ends with short lactide segments. The BMP-2 peptide with a positive index of hydrophobicity had a critical micelle concentration (CMC) and formed aggregates in aqueous solution. Based on simulation results, there was a slight increase in the concentration of free peptide in solution with 1000-fold increase in peptide concentration. The dose-osteogenic response curve of the BMP-2 peptide was in the 0.0005-0.005 mM range, and osteoinductive potential of the BMP-2 peptide was significantly less than that of BMP-2 protein even at 1000-fold higher concentrations, which was attributed to peptide aggregation. Further, the peptide or PEG-peptide aggregates had significantly higher interaction energy with the cell membrane compared with the free peptide, which led to a higher nonspecific interaction with the cell membrane and loss of osteoinductive potential. Conjugation of the BMP-2 peptide to PEG increased CMC and osteoinductive potential of the peptide whereas conjugation to lactide-capped PEG reduced CMC and osteoinductive potential of the peptide. Experimental and simulation results revealed that osteoinductive potential of the BMP-2 peptide is correlated with its CMC and the free peptide concentration in aqueous medium and not the total concentration.

Synergistic effects of magainin 2 and PGLa on their heterodimer formation, aggregation, and insertion into the bilayer

We performed coarse-grained molecular dynamics simulations of antimicrobial peptides PGLa and magainin 2 in lipid bilayers.