Comparative molecular dynamics simulations of the antimicrobial peptide CM15 in model lipid bilayers

Structure of the Bacterial Cell Envelope and Interactions with Antimicrobials: Insights from Molecular Dynamics Simulations

Bacteria have evolved over 3 billion years, shaping our intrinsic and symbiotic coexistence with these single-celled organisms. With rising populations of drug-resistant strains, the search for novel antimicrobials is an ongoing area of research. Advances in high-performance computing platforms have led to a variety of molecular dynamics simulation strategies to study the interactions of antimicrobial molecules with different compartments of the bacterial cell envelope of both Gram-positive and Gram-negative species. In this review, we begin with a detailed description of the structural aspects of the bacterial cell envelope. Simulations concerned with the transport and associated free energy of small molecules and ions through the outer membrane, peptidoglycan, inner membrane and outer membrane porins are discussed. Since surfactants are widely used as antimicrobials, a section is devoted to the interactions of surfactants with the cell wall and inner membranes. The review ends with a discussion on antimicrobial peptides and the insights gained from the molecular simulations on the free energy of translocation. Challenges involved in developing accurate molecular models and coarse-grained strategies that provide a trade-off between atomic details with a gain in sampling time are highlighted. The need for efficient sampling strategies to obtain accurate free energies of translocation is also discussed. Molecular dynamics simulations have evolved as a powerful tool that can potentially be used to design and develop novel antimicrobials and strategies to effectively treat bacterial infections.

Peptide Power: Mechanistic Insights into the Effect of Mitochondria-Targeted Tetrapeptides on Membrane Electrostatics from Molecular Simulations

Mitochondrial dysfunction is implicated in nine of the ten leading causes of death in the US, yet there are no FDA-approved therapeutics to treat it. Synthetic mitochondria-targeted peptides (MTPs), including the lead compound SS-31, offer promise, as they have been shown to restore healthy mitochondrial function and treat a variety of common diseases. At the cellular level, research has shown that MTPs accumulate strongly at the inner mitochondrial membrane (IMM), slow energy sinks (e.g., proton leaks), and improve ATP production. Modulation of electrostatic fields around the IMM has been implicated as a key aspect in the mechanism of action (MoA) of these peptides; however, molecular and mechanistic details have remained elusive. In this study, we employed all-atom molecular dynamics simulations (MD) to investigate the interactions of four MTPs with lipid bilayers and calculate their effect on structural and electrostatic properties. In agreement with previous experimental findings, we observed the modulation of the membrane surface and dipole potentials by MTPs. The simulations reveal that the MTPs achieve a reduction in the dipole potential by acting to disorder both lipid head groups and water layers proximal to the bilayer surface. We also find that MTPs decrease the bilayer thickness and increase the membrane's capacitance. These changes suggest that MTPs may enhance how much potential energy can be stored across the IMM at a given transmembrane potential difference. The MTPs also displace cations away from the bilayer surface, modulating the surface potential and offering an alternative mechanism for how these MTPs reduce mitochondrial energy sinks like proton leaks and mitigate Ca2+ accumulation stress. In conclusion, this study highlights the therapeutic potential of MTPs and underlines how interactions of MTPs with lipid bilayers serve as a fundamental component of their MoA.

Effect of CM15 on Supported Lipid Bilayer Probed by Atomic Force Microscopy

Antimicrobial peptides are key components of the immune system. These peptides affect the membrane in various ways; some form nano-sized pores, while others only produce minor defects. Since these peptides are increasingly important in developing antimicrobial drugs, understanding the mechanism of their interactions with lipid bilayers is critical. Here, using atomic force microscopy (AFM), we investigated the effect of a synthetic hybrid peptide, CM15, on the membrane surface comprising E. coli polar lipid extract. Direct imaging of supported lipid bilayers exposed to various concentrations of the peptide revealed significant membrane remodeling. We found that CM15 interacts with supported lipid bilayers and forms membrane-spanning defects very quickly. It is found that CM15 is capable of remodeling both leaflets of the bilayer. For lower CM15 concentrations, punctate void-like defects were observed, some of which re-sealed themselves as a function of time. However, for CM15 concentrations higher than 5 µM, the defects on the bilayers became so widespread that they disrupted the membrane integrity completely. This work enhances the understanding of CM15 interactions with the bacterial lipid bilayer.

Adsorption/Desorption of Cationic-Hydrophobic Peptides on Zwitterionic Lipid Bilayer Is Associated with the Possibility of Proton Transfer

Cell-penetrating peptides (CPPs) are short peptides built up from dominantly cationic and hydrophobic amino acid residues with a distinguished ability to pass through the cell membrane. Due to the possibility of linking and delivering the appropriate cargo at the desired location, CPPs are considered an economic and less invasive alternative to antibiotics. Besides knowing that their membrane passage mechanism is a complex function of CPP chemical composition, the ionic strength of the solution, and the membrane composition, all other details on how they penetrate cell membranes are rather vague. The aim of this study is to elucidate the ad(de)sorption of arginine-/lysine- and phenylalanine-rich peptides on a lipid membrane composed of 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) lipids. DSC and temperature-dependent UV-Vis measurements confirmed the impact of the adsorbed peptides on thermotropic properties of DPPC, but in an inconclusive way. On the other hand, FTIR spectra acquired at 30 °C and 50 °C (when DPPC lipids are found in the gel and fluid phase, respectively) unambiguously confirmed the proton transfer between particular titratable functional groups of R5F2/K5F2 that highly depend on their immediate surroundings (DPPC or a phosphate buffer). Molecular dynamic simulations showed that both peptides may adsorb onto the bilayer, but K5F2 desorbs more easily and favors the solvent, while R5F2 remains attached. The results obtained in this work highlight the importance of proton transfer in the design of CPPs with their desired cargo, as its charge and composition dictates the possibility of entering the cell.

A dual functional polypeptide with antibacterial and anti-inflammatory properties for the treatment of periodontitis

Periodontitis has been reported as the sixth most prevalent disease in human beings. This destructive disease is closely related to systemic diseases. Existing local drug delivery systems for periodontitis suffer from poor antibacterial effect and drug resistance. Inspired by the pathogenesis of periodontitis, we implemented a strategy to construct a dual functional polypeptide LL37-C15, which exhibited remarkable antibacterial effect against P. gingivalis and A. actinomycetemcomitans. In addition, LL37-C15 inhibits the release of pro-inflammatory cytokines by controlling the inflammatory pathway and reversing macrophage M1. Furthermore, the anti-inflammatory effect of LL37-C15 was also verified in vivo in a periodontitis rat model through the morphometry and histological observations of alveolar bone, hematoxylin-eosin, and Trap staining in gingival tissue. The results of molecular dynamics simulations showed that LL37-C15 could selectively destroy the bacterial cell membrane and protect the animal cell membrane in a self-destructive manner. The results showed that the polypeptide LL37-C15, as a novel promising therapeutic agent, exhibited a great potential for the periodontitis management. What's more, this dual functional polypeptide provides a promising strategy for building a multifunctional therapeutic platform against the inflammation and other diseases.

Molecular Dynamics Simulation of Antimicrobial Peptide CM15 in Staphylococcus Aureus and Escherichia coli Model Bilayer Lipid

Background

In animals and plants, antimicrobial peptides (AMPs) are crucial components of defense mechanisms, as they play a crucial role in innate immunity, which protects hosts from pathogenic bacteria. The CM15 has attracted considerable interest as a novel antibiotic against gram-negative and positive pathogens.

Objective

The aim of this study was to investigate the permeation potential of the CM15 with membrane bilayers of Staphylococcus aureus and Escherichia coli.

Material and Methods

The bilayer membranes of Escherichia coli and Staphylococcus aureus were modelled with the resemblance in lipid composition to its biological sample. This study followed Protein-Membrane Interaction (PMI) through successive applications of molecular dynamics simulation by GROMACS and CHARMM36 force field for two sets of 120-ns simulations.

Results

Significant results were obtained from analyzing the trajectory of the unsuccessful insertion of CM15 during simulation. Our data suggested that Lysine residues in CM15 and Cardiolipins in membrane leaflets play a crucial role in stability and interaction terms.

Conclusion

The obtained results strengthen the insertion possibility through the toroidal model, which should consider for further studies on AMPs interaction.

Removal and identification of external protein corona members from RBC-derived extracellular vesicles by surface manipulating antimicrobial peptides

In the last years, extracellular vesicles (EVs), secreted by various cells and body fluids have shown extreme potential in biomedical applications. Increasing number of studies suggest that a protein corona could adhere to the surface of EVs which can have a fundamental effect on their function, targeting and therapeutical efficacy. However, removing and identifying these corona members is currently a challenging task to achieve. In this study we have employed red blood cell-derived extracellular vesicles (REVs) as a model system and three membrane active antimicrobial peptides (AMPs), LL-37, FK-16 and CM15, to test whether they can be used to remove protein corona members from the surface of vesicles. These AMPs were reported to preferentially exert their membrane-related activity via one of the common helical surface-covering models and do not significantly affect the interior of lipid bilayer bodies. The interaction between the peptides and the REVs was followed by biophysical techniques, such as flow-linear dichroism spectroscopy which provided the effective applicable peptide concentration for protein removal. REV samples were then subjected to subsequent size exclusion chromatography and to proteomics analysis. Based on the comparison of control REVs with the peptide treated samples, seventeen proteins were identified as external protein corona members. From the three investigated AMPs, FK-16 can be considered as the best candidate to further optimize EV-related applicability of AMPs. Our results on the REV model system envisage that membrane active peptides may become a useful set of tools in engineering and modifying surfaces of EVs and other lipid-based natural particles.

Application of a deep generative model produces novel and diverse functional peptides against microbial resistance

Antimicrobial resistance could threaten millions of lives in the immediate future. Antimicrobial peptides (AMPs) are an alternative to conventional antibiotics practice against infectious diseases. Despite the potential contribution of AMPs to the antibiotic's world, their development and optimization have encountered serious challenges. Cutting-edge methods with novel and improved selectivity toward resistant targets must be established to create AMPs-driven treatments. Here, we present AMPTrans-lstm, a deep generative network-based approach for the rational design of AMPs. The AMPTrans-lstm pipeline involves pre-training, transfer learning, and module identification. The AMPTrans-lstm model has two sub-models, namely, (long short-term memory) LSTM sampler and Transformer converter, which can be connected in series to make full use of the stability of LSTM and the novelty of Transformer model. These elements could generate AMPs candidates, which can then be tailored for specific applications. By analyzing the generated sequence and trained AMPs, we prove that AMPTrans-lstm can expand the design space of the trained AMPs and produce reasonable and brand-new AMPs sequences. AMPTrans-lstm can generate functional peptides for antimicrobial resistance with good novelty and diversity, so it is an efficient AMPs design tool.

A Molecular Dynamics Study of Antimicrobial Peptide Interactions with the Lipopolysaccharides of the Outer Bacterial Membrane

With rising bacterial resistance, antimicrobial peptides (AMPs) have been widely investigated as potential antibacterial molecules to replace conventional antibiotics. Our understanding of the molecular mechanisms for membrane disruption are largely based on AMP interactions with the inner phospholipid bilayers of both Gram-negative and Gram-positive bacteria. Mechanisms for AMP translocation across the outer membrane of Gram-negative bacteria composed of lipopolysaccharides and the asymmetric lipid bilayer are complicated by the secondary structure adopted by the peptide in the different membrane environments. We have employed atomistic molecular dynamics and umbrella-sampling simulations with an aggregate duration of [Formula: see text] 6 microseconds to obtain the free energy landscape of CM15 peptide translocating through the lipopolysaccharide region of Gram-negative bacteria, E. coli. The peptide has a favorable binding-free energy (- 130 kJ mol[Formula: see text]) in the O-antigen region with a large barrier (150 kJ mol[Formula: see text]) at the interface between the anionic core saccharides and upper bilayer leaflet made up of lipid-A molecules. Restraint-free molecular dynamics simulations show that the random coil structure is favored over the helix in both the extracellular aqueous region and the cation-rich core-saccharide regions of the outer membrane. The peptide and membrane properties are analyzed at each of the 100 ns duration of the umbrella-sampling windows to illustrate changes in peptide length, orientation, and hydration. Our study provides insights into the free energy landscape for the insertion of the AMP CM15 in the outer membrane of Gram-negative bacteria, and we discuss the implications of our findings with the broader question of how AMPs overcome this barrier during antimicrobial activity.

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.

Biomimetic antimicrobial polymers—Design, characterization, antimicrobial, and novel applications

Biomimetic antimicrobial polymers have been an area of great interest as the need for novel antimicrobial compounds grows due to the development of resistance. These polymers were designed and developed to mimic naturally occurring antimicrobial peptides in both physicochemical composition and mechanism of action. These antimicrobial peptide mimetic polymers have been extensively investigated using chemical, biophysical, microbiological, and computational approaches to gain a deeper understanding of the molecular interactions that drive function. These studies have helped inform SARs, mechanism of action, and general physicochemical factors that influence the activity and properties of antimicrobial polymers. However, there are still lingering questions in this field regarding 3D structural patterning, bioavailability, and applicability to alternative targets. In this review, we present a perspective on the development and characterization of several antimicrobial polymers and discuss novel applications of these molecules emerging in the field. This article is categorized under: Therapeutic Approaches and Drug Discovery > Emerging Technologies Therapeutic Approaches and Drug Discovery > Nanomedicine for Infectious Disease.

Identification of Bacterial Membrane Selectivity of Romo1-Derived Antimicrobial Peptide AMPR-22 via Molecular Dynamics

The abuse or misuse of antibiotics has caused the emergence of extensively drug-resistant (XDR) bacteria, rendering most antibiotics ineffective and increasing the mortality rate of patients with bacteremia or sepsis. Antimicrobial peptides (AMPs) are proposed to overcome this problem; however, many AMPs have attenuated antimicrobial activities with hemolytic toxicity in blood. Recently, AMPR-11 and its optimized derivative, AMPR-22, were reported to be potential candidates for the treatment of sepsis with a broad spectrum of antimicrobial activity and low hemolytic toxicity. Here, we performed molecular dynamics (MD) simulations to clarify the mechanism of lower hemolytic toxicity and higher efficacy of AMPR-22 at an atomic level. We found four polar residues in AMPR-11 bound to a model mimicking the bacterial inner/outer membranes preferentially over eukaryotic plasma membrane. AMPR-22 whose polar residues were replaced by lysine showed a 2-fold enhanced binding affinity to the bacterial membrane by interacting with bacterial specific lipids (lipid A or cardiolipin) via hydrogen bonds. The MD simulations were confirmed experimentally in models that partially mimic bacteremia conditions in vitro and ex vivo. The present study demonstrates why AMPR-22 showed low hemolytic toxicity and this approach using an MD simulation would be helpful in the development of AMPs.

Effect of biphosphate salt on dipalmitoylphosphatidylcholine bilayer deformation by Tat polypeptide

Translocation of positively charged cell penetrating peptides (CPP) through cell membrane is important in drug delivery. Here we report all-atom molecular dynamics simulations to investigate how a biphosphate salt in a solvent affects the interaction of a CPP, HIV-1 Tat peptide with model dipalmitoylphosphatidylcholine (DPPC) lipid bilayer. Tat peptide has a large number of basic arginines and a couple of polar glutamines. We observe that in absence of salt, the basic residues of the polypeptide get localized in the vicinity of the membrane without altering the bilayer properties much; polypeptide induce local thinning of the bilayer membrane at the area of localization. In presence of biphosphate salt, the basic residues, dressed by the biphosphate ions, are repelled by the phosphate head groups of the lipid molecules. However, polar glutamine prefers to stay in the vicinity of the bilayer. This leads to larger local bilayer thickness at the contact point by the polar residue and non-uniform bilayer thickness profile. The thickness deformation of bilayer structure disappears upon mutating the polar residue, suggesting importance of the polar residue in bilayer deformation. Our studies point to control bilayer deformation by appropriate peptide sequence and solvent conditions.

Helical structure motifs made searchable for functional peptide design

The systematic design of functional peptides has technological and therapeutic applications. However, there is a need for pattern-based search engines that help locate desired functional motifs in primary sequences regardless of their evolutionary conservation. Existing databases such as The Protein Secondary Structure database (PSS) no longer serves the community, while the Dictionary of Protein Secondary Structure (DSSP) annotates the secondary structures when tertiary structures of proteins are provided. Here, we extract 1.7 million helices from the PDB and compile them into a database (Therapeutic Peptide Design database; TP-DB) that allows queries of compounded patterns to facilitate the identification of sequence motifs of helical structures. We show how TP-DB helps us identify a known purification-tag-specific antibody that can be repurposed into a diagnostic kit for Helicobacter pylori. We also show how the database can be used to design a new antimicrobial peptide that shows better Candida albicans clearance and lower hemolysis than its template homologs. Finally, we demonstrate how TP-DB can suggest point mutations in helical peptide blockers to prevent a targeted tumorigenic protein-protein interaction. TP-DB is made available at http://dyn.life.nthu.edu.tw/design/ .

Recent advances in the design of antimicrobial peptide conjugates

Antimicrobial peptides (AMPs) are ubiquitous host defense peptides characterized by antibiotic activity and lower propensity for developing resistance compared to classic antibiotics. While several AMPs have shown activity against antibiotic-sensitive...

Controlling Peptide Function by Directed Assembly Formation: Mechanistic Insights Using Multiscale Modeling on an Antimicrobial Peptide-Drug-Membrane System

Owing to their potential applicability against multidrug-resistant bacteria, antimicrobial peptides (AMPs) or host defense peptides (HDPs) gain increased attention. Besides diverse immunomodulatory roles, their classical mechanism of action mostly involves membrane disruption of microbes. Notably, their unbalanced overexpression has also been associated with host cell cytotoxicity in various diseases. Relatedly, AMPs can be subject to aggregate formation, either via self-assembly or together with other compounds, which has demonstrated a modulation effect on their biological functions, thus highly relevant both for drug targeting projects and understanding their in vivo actions. However, the molecular aspects of the related assembly formation are not understood. Here, we focused in detail on an experimentally studied AMP-drug system, i.e., CM15-suramin, and performed all-atom and coarse-grain (CG) simulations. Results obtained for all systems were in close line with experimental observations and indicate that the CM15-suramin aggregation is an energetically favorable and dynamic process. In the presence of bilayers, the peptide-drug assembly formation was highly dependent on lipid composition, and peptide aggregates themselves were also capable of binding to the membranes. Interestingly, longer CG simulations with zwitterionic membranes indicated an intermediate state in the presence of both AMP-drug assemblies and monomeric peptides located on the membrane surface. In sharp contrast, larger AMP-drug aggregates could not be detected with a negatively charged membrane, rather the AMPs penetrated its surface in a monomeric form, in line with previous in vitro observations. Considering experimental and theoretical results, it is promoted that in biological systems, cationic AMPs may often form associates with anionic compounds in a reversible manner, resulting in lower bioactivity. This is only mildly affected by zwitterionic membranes; however, membranes with a negative charge strongly alter the energetic preference of AMP assemblies, resulting in the dissolution of the complexes into the membrane. The phenomenon observed here at a molecular level can be followed in several experimental systems studied recently, where peptides interact with food colors, drug molecules, or endogenous compounds, which strongly indicates that reversible associate formation is a general phenomenon for these complexes. These results are hoped to be exploited in novel therapeutic strategies aiming to use peptides as drug targets and control AMP bioactivity by directed assembly formation.

Nanocapsule designs for antimicrobial resistance

The pressing need of new antimicrobial products is growing stronger, particularly because of widespread antimicrobial resistance, endangering our ability to treat common infections. The recent coronavirus pandemic has dramatically highlighted the necessity of effective antibacterial and antiviral protection. This work explores at the molecular level the mechanism of action of antibacterial nanocapsules assembled in virus-like particles, their stability and their interaction with mammal and antimicrobial model membranes. We use Molecular Dynamics with force-fields of different granularity and protein design strategies to study the stability, self-assembly and membrane poration properties of these nanocapsules.

The effect of acidic pH on the adsorption and lytic activity of the peptides Polybia-MP1 and its histidine-containing analog in anionic lipid membrane: a biophysical study by molecular dynamics and spectroscopy

Antimicrobial peptides (AMPs) are part of the innate immune system of many species. AMPs are short sequences rich in charged and non-polar residues. They act on the lipid phase of the plasma membrane without requiring membrane receptors. Polybia-MP1 (MP1), extracted from a native wasp, is a broad-spectrum bactericide, an inhibitor of cancer cell proliferation being non-hemolytic and non-cytotoxic. MP1 mechanism of action and its adsorption mode is not yet completely known. Its adsorption to lipid bilayer and lytic activity is most likely dependent on the ionization state of its two acidic and three basic residues and consequently on the bulk pH. Here we investigated the effect of bulk acidic (pH 5.5) and neutral pH (7.4) solution on the adsorption, insertion, and lytic activity of MP1 and its analog H-MP1 to anionic (7POPC:3POPG) model membrane. H-MP1 is a synthetic analog of MP1 with lysines replaced by histidines. Bulk pH changes could modulate this peptide efficiency. The combination of different experimental techniques and molecular dynamics (MD) simulations showed that the adsorption, insertion, and lytic activity of H-MP1 are highly sensitive to bulk pH in opposition to MP1. The atomistic details, provided by MD simulations, showed peptides contact their N-termini to the bilayer before the insertion and then lay parallel to the bilayer. Their hydrophobic faces inserted into the acyl chain phase disturb the lipid-packing.

Molecular Simulations Guidelines for Biological Nanomaterials: From Peptides to Membranes

In studying biological processes and focusing on the molecular mechanisms at the basis of these, molecular dynamics (MD) simulations have demonstrated to be a very useful tool for the past 50 years. This suite of computational methods calculates the time-dependent evolution of a molecular system using physics-based first principles. In this chapter, we give a brief introduction to the theory and practical use of molecular dynamics simulations, highlighting the different models and algorithms that have been developed to tackle specific problems, with a special focus on classical force fields. Some examples of how simulations have been used in the past will help the reader in discerning their power, limitations, and significance.

Prediction of amphipathic helix-membrane interactions with Rosetta

Amphipathic helices have hydrophobic and hydrophilic/charged residues situated on opposite faces of the helix. They can anchor peripheral membrane proteins to the membrane, be attached to integral membrane proteins, or exist as independent peptides. Despite the widespread presence of membrane-interacting amphipathic helices, there is no computational tool within Rosetta to model their interactions with membranes. In order to address this need, we developed the AmphiScan protocol with PyRosetta, which runs a grid search to find the most favorable position of an amphipathic helix with respect to the membrane. The performance of the algorithm was tested in benchmarks with the RosettaMembrane, ref2015_memb, and franklin2019 score functions on six engineered and 44 naturally-occurring amphipathic helices using membrane coordinates from the OPM and PDBTM databases, OREMPRO server, and MD simulations for comparison. The AmphiScan protocol predicted the coordinates of amphipathic helices within less than 3Å of the reference structures and identified membrane-embedded residues with a Matthews Correlation Constant (MCC) of up to 0.57. Overall, AmphiScan stands as fast, accurate, and highly-customizable protocol that can be pipelined with other Rosetta and Python applications.

Amphipathic helices are important targets as antibacterial peptides and as domains of membrane proteins that play a role in sensing the membrane environment. Understanding how amphipathic helices interact with membrane enables us to design better peptides and understand how membrane proteins use them to interact with their environment. However, there is a limited number of tools available for the modeling of amphipathic helices in membranes. Implicit membrane models can be used for this purpose as simplistic representations of the membrane environment. In this work, we developed the AmphiScan protocol that can be used to predict membrane coordinates of amphipathic helices starting with a helix structure in an implicit membrane environment. We benchmarked the performance of AmphiScan on engineered LK peptides, naturally-occurring amphipathic helices, and hydrophobic and hydrophilic peptides. Our approach provides a reliable and customizable tool to model amphipathic helix–membrane interactions, and pose a platform for the screening of amphipathic helix properties in silico.

Membrane Binding of Antimicrobial Peptides Is Modulated by Lipid Charge Modification

Peptide interactions with lipid bilayers play a key role in a range of biological processes and depend on electrostatic interactions between charged amino acids and lipid headgroups. Antimicrobial peptides (AMPs) initiate the killing of bacteria by binding to and destabilizing their membranes. The multiple peptide resistance factor (MprF) provides a defense mechanism for bacteria against a broad range of AMPs. MprF reduces the negative charge of bacterial membranes through enzymatic conversion of the anionic lipid phosphatidyl glycerol (PG) to either zwitterionic alanyl-phosphatidyl glycerol (Ala-PG) or cationic lysyl-phosphatidyl glycerol (Lys-PG). The resulting change in the membrane charge is suggested to reduce the binding of AMPs to membranes, thus impeding downstream AMP activity. Using coarse-grained molecular dynamics to investigate the effects of these modified lipids on AMP binding to model membranes, we show that AMPs have substantially reduced affinity for model membranes containing Ala-PG or Lys-PG. More than 5000 simulations in total are used to define the relationship between lipid bilayer composition, peptide sequence (using five different membrane-active peptides), and peptide binding to membranes. The degree of interaction of a peptide with a membrane correlates with the membrane surface charge density. Free energy profile (potential of mean force) calculations reveal that the lipid modifications due to MprF alter the energy barrier to peptide helix penetration of the bilayer. These results will offer a guide to the design of novel peptides, which addresses the issue of resistance via MprF-mediated membrane modification.

Interaction between Antimicrobial Peptide CM15 and a Model Cell Membrane Affected by CM15 Terminal Amidation and the Membrane Phase State

Antimicrobial peptides (AMPs) have been proposed as an effective class of antimicrobial agents against microorganisms. In this work, the interaction between an antimicrobial peptide, CM15, and a negatively charged phospholipid bilayer, DPPG, was studied via sum frequency generation (SFG) vibrational spectroscopy. Two structurally correlated characteristic variables were introduced to reveal the interaction mechanism/efficiency, i.e. C-terminal amidation and temperature variation (∼20 °C, room temperature, and ∼35 °C, close to human body temperature). Experimental results indicated that owing to the increased positive charge, C-terminal amidation resulted in rapid adsorption onto the bilayer surface and efficient disruption of the outer layer, exhibiting less ordered insertion orientation. The elevated temperature (from ∼20 °C to ∼35 °C) promoted the penetration of both the outer and inner leaflets by the peptides and finally led to the disruption of the whole bilayer owing to the enhanced fluidity of the bilayer. From the perspective of the interaction mechanism, this experimental study provides two practical cues to understand the disruption process of the negatively charged model biomembranes, which can lay the structural foundation for designing and developing high-efficiency antimicrobial peptides.

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.

Design, Screening, and Testing of Non-Rational Peptide Libraries with Antimicrobial Activity: In Silico and Experimental Approaches

One of the challenges of modern biotechnology is to find new routes to mitigate the resistance to conventional antibiotics. Antimicrobial peptides (AMPs) are an alternative type of biomolecules, naturally present in a wide variety of organisms, with the capacity to overcome the current microorganism resistance threat. Here, we reviewed our recent efforts to develop a new library of non-rationally produced AMPs that relies on bacterial genome inherent diversity and compared it with rationally designed libraries. Our approach is based on a four-stage workflow process that incorporates the interplay of recent developments in four major emerging technologies: artificial intelligence, molecular dynamics, surface-display in microorganisms, and microfluidics. Implementing this framework is challenging because to obtain reliable results, the in silico algorithms to search for candidate AMPs need to overcome issues of the state-of-the-art approaches that limit the possibilities for multi-space data distribution analyses in extremely large databases. We expect to tackle this challenge by using a recently developed classification algorithm based on deep learning models that rely on convolutional layers and gated recurrent units. This will be complemented by carefully tailored molecular dynamics simulations to elucidate specific interactions with lipid bilayers. Candidate AMPs will be recombinantly-expressed on the surface of microorganisms for further screening via different droplet-based microfluidic-based strategies to identify AMPs with the desired lytic abilities. We believe that the proposed approach opens opportunities for searching and screening bioactive peptides for other applications.

Planarity and out-of-plane vibrational modes of tryptophan and tyrosine in biomolecular modeling

Tryptophan and tyrosine are amino acids that play significant roles in the folding processes of proteins at water-membrane interfaces because of their amphipathic heteroaromatic rings. Employing appropriate heteroaromatic molecular structures is essential for obtaining accurate dynamics and predictive capabilities in molecular simulations of these amino acids. In this study, molecular dynamics simulations that applied the most recent version of the CHARMM36 force field were conducted on aqueous solutions of tryptophan and of tyrosine. Geometric analysis and dynamics quantified how aromatic rings deviated from planar structures and exhibited out-of-plane fluctuations. Radial distribution functions showed possible biological significance because the extent of ring planarity slightly affected local water concentrations near aromatic rings. Instantaneous all-atom normal mode analysis (NMA) and Fourier transformation of time autocorrelation functions of out-of-plane displacements were applied to study out-of-plane vibrations of atoms in these rings. The NMA started with minimum energy configurations and then averaged over fluctuations in aqueous solution. The frequencies and frequency patterns that were obtained for tryptophan and tyrosine with CHARMM36 differed from literature reports of Raman spectra, infrared spectra, and frequencies calculated using quantum mechanics, with some out-of-plane modes found at higher frequencies. Effects of imposing improper torsion potentials and changing torsion angle force constants were investigated for all atoms in the rings of tryptophan and tyrosine. Results show that these coarse force field variations only affect planarity and out-of-plane vibrations of atoms within the rings, and not other vibrations. Although increasing improper torsion force constants reduced deviations from aromatic ring planarity significantly, it increased out-of-plane mode frequencies. Reducing torsion angle force constants (with and without improper torsions) shifted modes to lower frequencies. A combination of decreasing most torsion angle force constants for ring atoms in both amino acids and including improper torsion forces attained frequencies and frequency patterns for out-of-plane normal modes that were more similar to the literature spectra. These force field variations decreased the extents of out-of-plane vibrations within the heteroaromatic rings of tryptophan, especially around the nitrogen atom in the ring, but not within the heteroaromatic ring of tyrosine. Conclusions were unaffected by the peptide endgroup, water, or simulation ensemble.

Self-Association of Antimicrobial Peptides: A Molecular Dynamics Simulation Study on Bombinin

Antimicrobial peptides (AMPs) are a diverse group of membrane-active peptides which play a crucial role as mediators of the primary host defense against microbial invasion. Many AMPs are found to be fully or partially disordered in solution and to acquire secondary structure upon interaction with a lipid membrane. Here, we report molecular dynamics simulations studies on the solution behaviour of a specific AMP, bombinin H2. We show that in monomeric form in water solution the peptide is somewhat disordered and preferably adopts a helix-loop-helix conformation. However, when more than a single monomer is placed in the solution, the peptides self-associate in aggregates. Within the aggregate, the peptides provide each other with an amphipathic environment that mimics the water–membrane interface, which allows them to adopt a single-helix structure. We hypothesise that this is the mechanism by which bombinin H2 and, possibly, other small linear AMPs reach the target membrane in a functional folded state and are able to effectively exert their antimicrobial action on it.

Studies on the Interaction of Alyteserin 1c Peptide and Its Cationic Analogue with Model Membranes Imitating Mammalian and Bacterial Membranes

Antimicrobial peptides (AMPs) are effector molecules of the innate immune system and have been isolated from multiple organisms. Their antimicrobial properties are due to the fact that they interact mainly with the anionic membrane of the microorganisms, permeabilizing it and releasing the cytoplasmic content. Alyteserin 1c (+2), an AMP isolated from Alytes obstetricans and its more cationic and hydrophilic analogue (+5) were synthesized using the solid phase method, in order to study the interaction with model membranes by calorimetric and spectroscopic assays. Differential scanning calorimetry (DSC) showed that both peptides had a strong effect when the membrane contained phosphatidylcholine (PC) alone or was mixed with phosphatidylglycerol (PG), increasing membrane fluidization. Attenuated total reflectance Fourier transform infrared spectroscopy (ATR-FTIR) was used to study the secondary structure of the peptide. Peptide +2 exhibited a transition from β-sheet/turns to β-sheet/α-helix structures after binding with model membranes, whereas peptide +5 had a transition from aggregation/unordered to β-sheet/α-helix structures after binding with membrane-contained PC. Interestingly, the latter showed a β-sheet structure predominantly in the presence of PG lipids. Additionally, molecular dynamics (MD) results showed that the carboxy-terminal of the peptide +5 has the ability to insert into the surface of the PC/PG membranes, resulting in the increase of the membrane fluidity.

Lipid Head Group Parameterization for GROMOS 54A8: A Consistent Approach with Protein Force Field Description

Membranes are a crucial component of both bacterial and mammalian cells, being involved in signaling, transport, and compartmentalization. This versatility requires a variety of lipid species to tailor the membrane’s behavior as needed, increasing the complexity of the system. Molecular dynamics simulations have been successfully applied to study model membranes and their interactions with proteins, elucidating some crucial mechanisms at the atomistic detail and thus complementing experimental techniques. An accurate description of the functional interplay of the diverse membrane components crucially depends on the selected parameters that define the adopted force field. A coherent parameterization for lipids and proteins is therefore needed. In this work, we propose and validate new lipid head group parameters for the GROMOS 54A8 force field, making use of recently published parametrizations for key chemical moieties present in lipids. We make use additionally of a new canonical set of partial charges for lipids, chosen to be consistent with the parameterization of soluble molecules such as proteins. We test the derived parameters on five phosphocholine model bilayers, composed of lipid patches four times larger than the ones used in previous studies, and run 500 ns long simulations of each system. Reproduction of experimental data like area per lipid and deuterium order parameters is good and comparable with previous parameterizations, as well as the description of liquid crystal to gel-phase transition. On the other hand, the orientational behavior of the head groups is more realistic for this new parameter set, and this can be crucial in the description of interactions with other polar molecules. For that reason, we tested the interaction of the antimicrobial peptide lactoferricin with two model membranes showing that the new parameters lead to a weaker peptide–membrane binding and give a more realistic outcome in comparing binding to antimicrobial versus mammal membranes.

Lipid Bilayer Composition Influences the Activity of the Antimicrobial Peptide Dermcidin Channel

Antimicrobial peptides (AMPs) carry great potential as new antibiotics against "superbugs." Dermcidin (DCD), a broad-spectrum AMP in human sweat, has been recently crystallized in its oligomeric state and showed channel-like properties. In this work, we performed multiscale molecular dynamics simulations to study how the membrane composition influences the behavior of a transmembrane pore formed by the DCD oligomer in the hope of revealing the origin of the membrane selectivity of this AMP toward bacteria. Our results indicate that bilayers composed of various lipids (DMPC, DPPC, and DSPC) with different thicknesses result in different orientations of the DCD oligomer when embedded in lipid bilayers. The thicker the bilayer, the less tilted the channel. Cholesterol makes the bilayers more rigid and thicker, which also affects the orientation of the channel. Furthermore, we observed that the predicted conductance of the channel from computational electrophysiology simulations is related to its orientation in the lipid bilayer: the larger the tilt, the larger the conductance. Our results indicate that the membrane composition has a significant influence on the activity of the DCD channel, with thicker, cholesterol-rich membranes showing lower conductance than that of thinner membranes.

Molecular Dynamics Investigation into the Effect of Zinc(II) on the Structure and Membrane Interactions of the Antimicrobial Peptide Clavanin A

Clavanin A (ClavA) is an antimicrobial peptide (AMP) whose antimicrobial activity is enhanced in the presence of Zn(II) ions. The antimicrobial activity of ClavA has been shown to increase 16-fold in the presence of Zn(II) ions. In this study, we investigate the potential sources of this enhancement, namely, the effect of Zn(II) binding on the helical conformation of ClavA and on the ClavA interaction with a model for gram-negative bacterial membranes. In addition, we investigate the effect of Zn(II) on the membrane mechanical properties. We employed all-atom equilibrium molecular dynamics simulations initiated from both fully helical and random coil structures of ClavA. We observe that Zn(II) can stabilize an existing helical conformation in the Zn(II)-binding region, but we do not observe induction of helical conformations in systems initiated in random coil configurations. Zn(II) binding to ClavA provides more favorable electrostatics for membrane association in the C-terminal region. This is evidenced by longer and stronger C-terminal-lipid interactions. Zn(II) is also capable of modulating the membrane properties in a manner which favors ClavA insertion and the potential for enhanced translocation into the cell. This work provides insights into the role of divalent metal cations in the antimicrobial activity of ClavA. This information can be used for the development of synthetic AMPs containing motifs that can bind metals (metalloAMPs) for therapeutic and medical purposes.

Antimicrobial peptides: Promising alternatives in the post feeding antibiotic era

Antimicrobial peptides (AMPs), critical components of the innate immune system, are widely distributed throughout the animal and plant kingdoms. They can protect against a broad array of infection-causing agents, such as bacteria, fungi, parasites, viruses, and tumor cells, and also exhibit immunomodulatory activity. AMPs exert antimicrobial activities primarily through mechanisms involving membrane disruption, so they have a lower likelihood of inducing drug resistance. Extensive studies on the structure-activity relationship have revealed that net charge, hydrophobicity, and amphipathicity are the most important physicochemical and structural determinants endowing AMPs with antimicrobial potency and cell selectivity. This review summarizes the recent advances in AMPs development with respect to characteristics, structure-activity relationships, functions, antimicrobial mechanisms, expression regulation, and applications in food, medicine, and animals.

Designing effective anticancer-radiopeptides. A Molecular Dynamics study of their interaction with model tumor and healthy cell membranes

One of the greatest merit of the use of radiopeptides in oncology is their selectivity which, however, brings about the drawback that each radiopeptide is specific for a given tumor type. To overcome this problem the direction currently taken in drug design is that of radiolabelling peptide hormones (or their analogues), relying on their intrinsic ability to bind to specific receptors in precise areas of the human body, at the cost, however, of a poor selectivity against healthy cells. We present here an extensive Molecular Dynamics study of a promising alternative inspired by the mechanism through which antimicrobial peptides interact with the negatively charged bacterial membranes. Appropriately modifying the human antimicrobial peptide, LL-37, we designed a functionalized radionuclide carrier capable of binding more strongly to the negatively charged (model) tumor membranes than to the neutral healthy ones. The mechanism behind this behaviour relies on the fact that at the slight acidic pH surrounding tumor tissues the histidines belonging to the peptide get protonated thus making it positively charged. We have investigated by an extended numerical study the way in which this artificial peptide interacts with models of tumor and healthy cell membranes, proving by Potential Mean Force calculations that the affinity of the peptide to model tumor membranes is significantly larger than to healthy ones. These features (high affinity and generic tumor selectivity) recommend antimicrobial derived customized carriers as promising theranostic constructs in cancer diagnostic and therapy.

Molecular Dynamics Simulations of Human Antimicrobial Peptide LL-37 in Model POPC and POPG Lipid Bilayers

Cathelicidins are a large family of cationic antimicrobial peptides (AMPs) found in mammals with broad spectrum antimicrobial activity. LL-37 is the sole amphipathic α-helical AMP from human Cathelicidins family. In addition to its bactericidal capability, LL-37 has antiviral, anti-tumor, and immunoregulatory activity. Despite many experimental studies, its molecular mechanism of action is not yet fully understood. Here, we performed three independent molecular dynamics simulations (600 ns or more) of a LL-37 peptide in the presence of 256 lipid bilayers with 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoglycerol (POPG) mimicking bacterial and 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) mimicking mammalian membranes. We found that LL-37 can be quickly absorbed onto the POPG bilayer without loss of its helical conformation in the core region and with the helix lying in parallel to the bilayer. The POPG bilayer was deformed. In contrast, LL-37 is slower in reaching the POPC surface and loss much of its helical conformation during the interaction with the bilayer. LL-37 only partially entered the POPC bilayer without significant deformation of the membrane. The observed difference for different bilayers is largely due to the fact that LL-37 is positively charged, POPG is negatively charged, and POPC is neutral. Our simulation results demonstrated the initial stage of disruption of the bacterial membrane by LL-37 in atomic details. Comparison to experimental results on LL-37 and simulation studies in other systems was made.

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

Computational design of the helical hairpin structure of membrane-active antibacterial peptides based on RSV glycoprotein epitope scaffold

Peptides with helical hairpin conformation have been found to possess potent membrane activity and can be exploited as the structural scaffold of antibacterial peptides (ABPs). Here, we attempted to computationally design membrane-active ABPs based on the helical hairpin motif of respiratory syncytial virus (RSV) glycoprotein epitope. Dynamics simulations revealed that the epitope peptide Rfe (net charge = -1) cannot effectively interact with and permeabilize bacterial membrane due to the electrostatic repulsion between the negatively charged peptide and anionic membrane surface. The native Rfe can be modified to a cationic peptide Rfe-KKK (net charge = +6) by triple mutation of its positively charged residues Glu256, Asp263 and Asp269 to a basic lysine as well as by C-terminal amidation. As might be expected, the modified peptide was able to target membrane surface with a moderate antibacterial potency (MIC = 50-100 μg/ml). Next, a cyclized version of the linear Rfe-KKK was generated, termed as cycRfe-KKK, which was observed to have improved membrane activity and increased antibacterial potency (MIC < 50 μg/ml) by pre-stabilizing amphipathic hairpin conformation of the peptide.

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.

Antimicrobial Peptide Potency is Facilitated by Greater Conformational Flexibility when Binding to Gram-negative Bacterial Inner Membranes

The interaction of antimicrobial peptides (AMPs) with the inner membrane of Gram-negative bacteria is a key determinant of their abilities to exert diverse bactericidal effects. Here we present a molecular level understanding of the initial target membrane interaction for two cationic α-helical AMPs that share structural similarities but have a ten-fold difference in antibacterial potency towards Gram-negative bacteria. The binding and insertion from solution of pleurocidin or magainin 2 to membranes representing the inner membrane of Gram-negative bacteria, comprising a mixture of 128 anionic and 384 zwitterionic lipids, is monitored over 100 ns in all atom molecular dynamics simulations. The effects of the membrane interaction on both the peptide and lipid constituents are considered and compared with new and published experimental data obtained in the steady state. While both magainin 2 and pleurocidin are capable of disrupting bacterial membranes, the greater potency of pleurocidin is linked to its ability to penetrate within the bacterial cell. We show that pleurocidin displays much greater conformational flexibility when compared with magainin 2, resists self-association at the membrane surface and penetrates further into the hydrophobic core of the lipid bilayer. Conformational flexibility is therefore revealed as a key feature required of apparently α-helical cationic AMPs for enhanced antibacterial potency.

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.

Membrane Mediated Antimicrobial and Antitumor Activity of Cathelicidin 6: Structural Insights from Molecular Dynamics Simulation on Multi-Microsecond Scale

The cathelicidin derived bovine antimicrobial peptide BMAP27 exhibits an effective microbicidal activity and moderate cytotoxicity towards erythrocytes. Irrespective of its therapeutic and multidimensional potentiality, the structural studies are still elusive. Moreover, the mechanism of BMAP27 mediated pore formation in heterogeneous lipid membrane systems is poorly explored. Here, we studied the effect of BMAP27 in model cell-membrane systems such as zwitterionic, anionic, thymocytes-like (TLM) and leukemia-like membranes (LLM) by performing molecular dynamics (MD) simulation longer than 100 μs. All-atom MD studies revealed a stable helical conformation in the presence of anionic lipids, however, significant loss of helicity was identified in TLM and zwitterionic systems. A peptide tilt (~45˚) and central kink (at residue F10) was found in anionic and LLM models, respectively, with an average membrane penetration of < 0.5 nm. Coarse-grained (CG) MD analysis on a multi-μs scale shed light on the membrane-dependent peptide and lipid organization. Stable micelle and end-to-end like oligomers were formed in zwitterionic and TLM models, respectively. In contrast, unstable oligomer formation and monomeric BMAP27 penetration were observed in anionic and LLM systems with selective anionic lipid aggregation (in LLM). Peptide penetration up to ~1.5 nm was observed in CG-MD systems with the BMAP27 C-terminal oriented towards the bilayer core. Structural inspection suggested membrane penetration by micelle/end-to-end like peptide oligomers (carpet-model like) in the zwitterionic/TLM systems, and transmembrane-mode (toroidal-pore like) in the anionic/LLM systems, respectively. Structural insights and energetic interpretation in BMAP27 mutant highlighted the role of F10 and hydrophobic residues in mediating a membrane-specific peptide interaction. Free energy profiling showed a favorable (-4.58 kcal mol^-1 for LLM) and unfavorable (+0.17 kcal mol^-1 for TLM) peptide insertion in anionic and neutral systems, respectively. This determination can be exploited to regulate cell-specific BMAP27 cytotoxicity for the development of potential drugs and antibiotics.

In silico pharmacology: Drug membrane partitioning and crossing

Over the past decade, molecular dynamics (MD) simulations have become particularly powerful to rationalize drug insertion and partitioning in lipid bilayers. MD simulations efficiently support experimental evidences, with a comprehensive understanding of molecular interactions driving insertion and crossing. Prediction of drug partitioning is discussed with respect to drug families (anesthetics; β-blockers; non-steroidal anti-inflammatory drugs; antioxidants; antiviral drugs; antimicrobial peptides). To accurately evaluate passive permeation coefficients turned out to be a complex theoretical challenge; however the recent methodological developments based on biased MD simulations are particularly promising. Particular attention is paid to membrane composition (e.g., presence of cholesterol), which influences drug partitioning and permeation. Recent studies concerning in silico models of membrane proteins involved in drug transport (influx and efflux) are also reported here. These studies have allowed gaining insight in drug efflux by, e.g., ABC transporters at an atomic resolution, explicitly accounting for the mandatory forces induced by the surrounded lipid bilayer. Large-scale conformational changes were thoroughly analyzed.

Role of charged lipids in membrane structures - Insight given by simulations

Lipids and proteins are the main components of cell membranes. It is becoming increasingly clear that lipids, in addition to providing an environment for proteins to work in, are in many cases also able to modulate the structure and function of those proteins. Particularly charged lipids such as phosphatidylinositols and phosphatidylserines are involved in several examples of such effects. Molecular dynamics simulations have proved an invaluable tool in exploring these aspects. This so-called computational microscope can provide both complementing explanations for the experimental results and guide experiments to fruitful directions. In this paper, we review studies that have utilized molecular dynamics simulations to unravel the roles of charged lipids in membrane structures. We focus on lipids as active constituents of the membranes, affecting both general membrane properties as well as non-lipid membrane components, mainly proteins. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.

Computer modelling studies of the bilayer/water interface

This review summarises high resolution studies on the interface of lamellar lipid bilayers composed of the most typical lipid molecules which constitute the lipid matrix of biomembranes. The presented results were obtained predominantly by computer modelling methods. Whenever possible, the results were compared with experimental results obtained for similar systems. The first and main section of the review is concerned with the bilayer-water interface and is divided into four subsections. The first describes the simplest case, where the interface consists only of lipid head groups and water molecules and focuses on interactions between the lipid heads and water molecules; the second describes the interface containing also mono- and divalent ions and concentrates on lipid-ion interactions; the third describes direct inter-lipid interactions. These three subsections are followed by a discussion on the network of direct and indirect inter-lipid interactions at the bilayer interface. The second section summarises recent computer simulation studies on the interactions of antibacterial membrane active compounds with various models of the bacterial outer membrane. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.

Interactions of a class IIb bacteriocin with a model lipid bilayer, investigated through molecular dynamics simulations

The emergence of antibiotic resistant microorganisms poses an alarming threat to global health. Antimicrobial peptides (AMPs) are considered a possible effective alternative to conventional antibiotic therapies. An understanding of the mechanism of action of AMPs is needed in order to better control and optimize their bactericidal activity. Plantaricin EF is a heterodimeric AMP, consisting of two peptides Plantaricin E (PlnE) and Plantaricin F (PlnF). We studied the behavior of these peptides on the surface of a model lipid bilayer. We identified the residues that facilitate peptide-peptide interactions. We also identified residues that mediate interactions of the dimer with the membrane. PlnE interacts with the membrane through amino acids at both its termini, while only the N terminus of PlnF approaches the membrane. By comparing the activity of single-site mutants of the two-peptide bacteriocin and the simulations of the bacteriocin on the surface of a model lipid bilayer, structure activity relationships are proposed. These studies allow us to generate hypotheses that relate biophysical interactions observed in simulations with the experimentally measured activity. We find that single-site amino acid substitutions result in markedly stronger antimicrobial activity when they strengthen the interactions between the two peptides, while, concomitantly, they weaken peptide-membrane association. This effect is more pronounced in the case of the PlnE mutant (G20A), which interacts the strongest with PlnF and the weakest with the membrane while displaying the highest activity.