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

Investigating the Influence of Photoswitchable Lipids on the Structure and Dynamics of Lipid Membranes: Fundamentals and Potential Applications

In this work, we delve into the impact of photoisomerization of photoswitchable lipids (PSLs) on the membrane structure and dynamics at a molecular level. Through all-atom molecular dynamics simulations, we explore how UV irradiation-induced trans-to-cis isomerization of these lipids, particularly the azobenzene-derivatized phosphatidylcholine (AzoPC) lipid, influences the structure and dynamics of a simplified lipid membrane, mimicking those of E. coli bacteria across different temperatures. Our findings align with previous experimental observations regarding membrane properties and offer insights into localized effects and microscopic heterogeneity. Additionally, we estimate the relaxation time scale of the lipid membrane following AzoPC photoisomerization. Moreover, we demonstrate the feasibility of photoactivated drug release, exemplified by the controlled liberation of doxorubicin, an anticancer agent, through the membrane, suggesting the potential of PSLs in engineering photoactivated liposomes, coined as photoazosomes, for precise targeted drug delivery applications.

Processes and mechanisms underlying burst of giant unilamellar vesicles induced by antimicrobial peptides and compounds

To clarify the damage of lipid bilayer region in bacterial cell membrane caused by antimicrobial peptides (AMPs) and antimicrobial compounds (AMCs), their interactions with giant unilamellar vesicles (GUVs) of various lipid compositions have been examined. The findings revealed two main causes for the leakage: nanopore formation in the membrane and burst of GUVs. Although GUV burst has been explained previously based on the carpet model, the supporting evidence is limited. In this review, to better clarify the mechanism of GUV burst by AMPs, AMCs, and other membrane-active peptides, we described current knowledge of the conditions, characteristics, and detailed processes of GUV burst and the changes in the shape of the GUVs during burst. We identified several physical factors that affect GUV burst, such as membrane tension, electrostatic interaction, structural changes of GUV membrane such as membrane folding, and oil in the membrane. We also clarified one of the physical mechanisms underlying the instability of lipid bilayers that are associated with leakage in the carpet model. Based on these results, we propose a mechanism underlying some types of GUV burst induced by these substances: the growth of a nanopore to a micropore, resulting in GUV burst.

Discovery of potential anti- Staphylococcus aureus natural products and their mechanistic studies using machine learning and molecular dynamic simulations

The structure-activity analysis (SAR) and machine learning were used to investigate potential anti-S. aureus agents in a faster method. In this study, 24 oxygenated benzene ring components with S. aureus inhibition capacity were confirmed by literature exploring and in-house experiments, and the SAR analysis suggested that the hydroxyl group position may affect the anti-S. aureus activity. The 2D-MLR-QSAR model with 9 descriptors was further evaluated as the best model among the 21 models. After that, hesperetic acid and 2-HTPA were further explored and evaluated as the potential anti-S. aureus agents screening in the natural product clustering library through the best QSAR model calculation. The antibacterial capacities of hesperetic acid and 2-HTPA had been investigated and proved the similar predictive pMIC value resulting from the QSAR model. Besides, the two novel components were able to inhibit the growth of S. aureus by disrupting the cell membrane through the molecular dynamics simulation (MD), which further evidenced by scanning electron microscopy (SEM) test and PI dye results. Overall, these results are highly suggested that QSAR can be used to predict the antibacterial agents targeting S. aureus, which provides a new paradigm to research the molecular structure-antibacterial capacity relationship.

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.

Structural characterization of the antimicrobial peptides myxinidin and WMR in bacterial membrane mimetic micelles and bicelles

Antimicrobial peptides are a promising class of potential antibiotics that interact selectively with negatively charged lipid bilayers. This paper presents the structural characterization of the antimicrobial peptides myxinidin and WMR associated with bacterial membrane mimetic micelles and bicelles by NMR, CD spectroscopy, and molecular dynamics simulations. Both peptides adopt a different conformation in the lipidic environment than in aqueous solution. The location of the peptides in micelles and bicelles has been studied by paramagnetic relaxation enhancement experiments with paramagnetic tagged 5- and 16-doxyl stearic acid (5-/16-SASL). Molecular dynamics simulations of multiple copies of the peptides were used to obtain an atomic level of detail on membrane-peptide and peptide-peptide interactions. Our results highlight an essential role of the negatively charged membrane mimetic in the structural stability of both myxinidin and WMR. The peptides localize predominantly in the membrane's headgroup region and have a noticeable membrane thinning effect on the overall bilayer structure. Myxinidin and WMR show a different tendency to self-aggregate, which is also influenced by the membrane composition (DOPE/DOPG versus DOPE/DOPG/CL) and can be related to the previously observed difference in the ability of the peptides to disrupt different types of model membranes.

Structure, Function, and Physicochemical Properties of Pore-forming Antimicrobial Peptides

Antimicrobial peptides (AMPs), a class of antimicrobial agents, possess considerable potential to treat various microbial ailments. The broad range of activity and rare complete bacterial resistance to AMPs make them ideal candidates for commercial development. These peptides with widely varying compositions and sources share recurrent structural and functional features in mechanisms of action. Studying the mechanisms of AMP activity against bacteria may lead to the development of new antimicrobial agents that are more potent. Generally, AMPs are effective against bacteria by forming pores or disrupting membrane barriers. The important structural aspects of cytoplasmic membranes of pathogens and host cells will also be outlined to understand the selective antimicrobial actions. The antimicrobial activities of AMPs are related to multiple physicochemical properties, such as length, sequence, helicity, charge, hydrophobicity, amphipathicity, polar angle, and also self-association. These parameters are interrelated and need to be considered in combination. So, gathering the most relevant available information will help to design and choose the most effective AMPs.

How Cell-Penetrating Peptides Behave Differently from Pore-Forming Peptides: Structure and Stability of Induced Transmembrane Pores

Peptide-induced transmembrane pore formation is commonplace in biology. Examples of transmembrane pores include pores formed by antimicrobial peptides (AMPs) and cell-penetrating peptides (CPPs) in bacterial membranes and eukaryotic membranes, respectively. In general, however, transmembrane pore formation depends on peptide sequences, lipid compositions, and intensive thermodynamic variables and is difficult to observe directly under realistic solution conditions, with structures that are challenging to measure directly. In contrast, the structure and phase behavior of peptide-lipid systems are relatively straightforward to map out experimentally for a broad range of conditions. Cubic phases are often observed in systems involving pore-forming peptides; however, it is not clear how the structural tendency to induce negative Gaussian curvature (NGC) in such phases is quantitatively related to the geometry of biological pores. Here, we leverage the theory of anisotropic inclusions and devise a facile method to estimate transmembrane pore sizes from geometric parameters of cubic phases measured from small-angle X-ray scattering (SAXS) and show that such estimates compare well with known pore sizes. Moreover, our model suggests that although AMPs can induce stable transmembrane pores for membranes with a broad range of conditions, pores formed by CPPs are highly labile, consistent with atomistic simulations.

Molecular Dynamics Simulations of Drug-Conjugated Cell-Penetrating Peptides

Cell-penetrating peptides (CPPs) are small peptides capable of translocating through biological membranes carrying various attached cargo into cells and even into the nucleus. They may also participate in transcellular transport. Our in silico study intends to model several peptides and their conjugates. We have selected three CPPs with a linear backbone, including penetratin, a naturally occurring oligopeptide; two of its modified sequence analogues (6,14-Phe-penetratin and dodeca-penetratin); and three natural CPPs with a cyclic backbone: Kalata B1, the Sunflower trypsin inhibitor 1 (SFT1), and Momordica cochinchinensis trypsin inhibitor II (MCoTI-II). We have also built conjugates with the small-molecule drug compounds doxorubicin, zidovudine, and rasagiline for each peptide. Molecular dynamics (MD) simulations were carried out with explicit membrane models. The analysis of the trajectories showed that the interaction of penetratin with the membrane led to spectacular rearrangements in the secondary structure of the peptide, while cyclic peptides remained unchanged due to their high conformational stability. Membrane-peptide and membrane-conjugate interactions have been identified and compared. Taking into account well-known examples from the literature, our simulations demonstrated the utility of computational methods for CPP complexes, and they may contribute to a better understanding of the mechanism of penetration, which could serve as the basis for delivering conjugated drug molecules to their intracellular targets.

How cell penetrating peptides behave differently from pore forming peptides: structure and stability of induced transmembrane pores

Peptide induced trans-membrane pore formation is commonplace in biology. Examples of transmembrane pores include pores formed by antimicrobial peptides (AMPs) and cell penetrating peptides (CPPs) in bacterial membranes and eukaryotic membranes, respectively. In general, however, transmembrane pore formation depends on peptide sequences, lipid compositions and intensive thermodynamic variables and is difficult to observe directly under realistic solution conditions, with structures that are challenging to measure directly. In contrast, the structure and phase behavior of peptide-lipid systems are relatively straightforward to map out experimentally for a broad range of conditions. Cubic phases are often observed in systems involving pore forming peptides; however, it is not clear how the structural tendency to induce negative Gaussian curvature (NGC) in such phases is quantitatively related to the geometry of biological pores. Here, we leverage the theory of anisotropic inclusions and devise a facile method to estimate transmembrane pore sizes from geometric parameters of cubic phases measured from small angle X-ray scattering (SAXS) and show that such estimates compare well with known pore sizes. Moreover, our model suggests that whereas AMPs can induce stable transmembrane pores for membranes with a broad range of conditions, pores formed by CPPs are highly labile, consistent with atomistic simulations.

A top-down and bottom-up combined strategy for parameterization of coarse-grained force fields for phospholipids

Coarse-graining (CG) molecular dynamics (MD) simulations are widely used in interpreting experimental observations and predicting assembly morphology as well as collective behaviour but also face the problem of poor accuracy. A main issue is that cross-termed interactions between different CG beads are inadequately parameterized. This work proposes a novel top-down and bottom-up combined strategy to parameterize both self- and cross-termed interactions of zwitterionic phospholipids in water solution based on a piecewise Morse potential describing nonbonded van der Waals interactions. The self-interacting force parameters were optimized by matching experimental density, heat vapourization, and surface tension in a top-down manner, while the cross-termed interactions were optimized by fitting pseudo properties obtained from atomistic simulations in a bottom-up way, including mixing density, intermolecular energy, and radial mixing coefficient. The transferability of the CG force field (FF) was confirmed by reproducing a variety of structural and thermodynamic properties of lipid membranes in both liquid and gel phases. This FF can well depict vesicle self-assembly and vesicle fusion processes. Matching pseudo properties opens a new way to develop CG FF with increased accuracy and transferability.

Molecular dynamics simulations reveal the importance of amyloid-beta oligomer β-sheet edge conformations in membrane permeabilization

Oligomeric aggregates of the amyloid-beta peptide(1-42) (Aβ42) are regarded as a primary cause of cytotoxicity related to membrane damage in Alzheimer's disease. However, a dynamical and structural characterization of pore-forming Aβ42 oligomers at atomic detail has not been feasible. Here, we used Aβ42 oligomer structures previously determined in a membrane-mimicking environment as putative model systems to study the pore formation process in phospholipid bilayers with all-atom molecular dynamics simulations. Multiple Aβ42 oligomer sizes, conformations, and N-terminally truncated isoforms were investigated on the multi-μs time scale. We found that pore formation and ion permeation occur via edge conductivity and exclusively for β-sandwich structures that feature exposed side-by-side β-strand pairs formed by residues 9 to 21 of Aβ42. The extent of pore formation and ion permeation depends on the insertion depth of hydrophilic residues 13 to 16 (HHQK domain) and thus on subtle differences in the overall stability, orientation, and conformation of the aggregates in the membrane. Additionally, we determined that backbone carbonyl and polar side-chain atoms from the edge strands directly contribute to the coordination sphere of the permeating ions. Furthermore, point mutations that alter the number of favorable side-chain contacts correlate with the ability of the Aβ42 oligomer models to facilitate ion permeation in the bilayer center. Our findings suggest that membrane-inserted, layered β-sheet edges are a key structural motif in pore-forming Aβ42 oligomers independent of their size and play a pivotal role in aggregate-induced membrane permeabilization.

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.

Simulation Study of the Effect of Antimicrobial Peptide Associations on the Mechanism of Action with Bacterial and Eukaryotic Membranes

Antimicrobial peptides (AMPs) can be directed to specific membranes based on differences in lipid composition. In this study, we performed atomistic and coarse-grained simulations of different numbers of the designed AMP adepantin-1 with a eukaryotic membrane, cytoplasmic Gram-positive and Gram-negative membranes, and an outer Gram-negative membrane. At the core of adepantin-1's behavior is its amphipathic α-helical structure, which was implemented in its design. The amphipathic structure promotes rapid self-association of peptide in water or upon binding to bacterial membranes. Aggregates initially make contact with the membrane via positively charged residues, but with insertion, the hydrophobic residues are exposed to the membrane's hydrophobic core. This adaptation alters the aggregate's stability, causing the peptides to diffuse in the polar region of the membrane, mostly remaining as a single peptide or pairing up to form an antiparallel dimer. Thus, the aggregate's proposed role is to aid in positioning the peptide into a favorable conformation for insertion. Simulations revealed the molecular basics of adepantin-1 binding to various membranes, and highlighted peptide aggregation as an important factor. These findings contribute to the development of novel anti-infective agents to combat the rapidly growing problem of bacterial resistance to antibiotics.

Bacterial Membranes Are More Perturbed by the Asymmetric Versus Symmetric Loading of Amphiphilic Molecules

Characterizing the biophysical properties of bacterial membranes is critical for understanding the protective nature of the microbial envelope, interaction of biological membranes with exogenous materials, and designing new antibacterial agents. Presented here are molecular dynamics simulations for two cationic quaternary ammonium compounds, and the anionic and nonionic form of a fatty acid molecule interacting with a Staphylococcus aureus bacterial inner membrane. The effect of the tested materials on the properties of the model membranes are evaluated with respect to various structural properties such as the lateral pressure profile, lipid tail order parameter, and the bilayer’s electrostatic potential. Conducting asymmetric loading of molecules in only one leaflet, it was observed that anionic and cationic amphiphiles have a large impact on the Staphylococcus aureus membrane’s electrostatic potential and lateral pressure profile as compared to a symmetric distribution. Nonintuitively, we find that the cationic and anionic molecules induce a similar change in the electrostatic potential, which points to the complexity of membrane interfaces, and how asymmetry can induce biophysical consequences. Finally, we link changes in membrane structure to the rate of electroporation for the membranes, and again find a crucial impact of introducing asymmetry to the system. Understanding these physical mechanisms provides critical insights and viable pathways for the rational design of membrane-active molecules, where controlling the localization is key.

Small Amphiphilic Peptides: Activity Against a Broad Range of Drug-Resistant Bacteria and Structural Insight into Membranolytic Properties

We report the synthesis and antibacterial activities of a series of amphiphilic membrane-active peptides composed, in part, of various nongenetically coded hydrophobic amino acids. The lead cyclic peptides, 8C and 9C, showed broad-spectrum activity against drug-resistant Gram-positive (minimum inhibitory concentration (MIC) = 1.5-6.2 μg/mL) and Gram-negative (MIC = 12.5-25 μg/mL) bacteria. The cytotoxicity study showed the predominant lethal action of the peptides against bacteria as compared with mammalian cells. A plasma stability study revealed approximately 2-fold higher stability of lead cyclic peptides as compared to their linear counterparts after 24 h of incubation. A calcein dye leakage experiment revealed the membranolytic effect of the cyclic peptides. Nuclear magnetic resonance spectroscopy and molecular dynamics simulation studies of the interaction of the peptides with the phospholipid bilayer provided a solid structural basis to explain the membranolytic action of the peptides with atomistic details. These results highlight the potential of newly designed amphiphilic peptides as the next generation of peptide-based antibiotics.

Molecular Dynamics Studies on the Bacterial Membrane Pore Formation by Small Molecule Antimicrobial Agents

Antimicrobial peptides (AMPs) act on the membrane bilayer of pathogens, causing leakage in the membrane and cell death. Amphiphilic kaempferol derivatives possessing basic functional groups show excellent antibacterial activities, which has been proven through experimental techniques. These compounds are known to target negatively charged bacterial membranes. However, the detailed mechanism of action and their structure-activity relationship are not clear. In this work, we reported theoretical investigation on the mechanism of action of two previously reported kaempferol derivatives on a DMPC/DMPG mixed bilayer. Despite the rigid structure of the compounds when compared to AMPs, spontaneous pore formation in the membrane was not observed in 400 ns molecular dynamics (MD) simulations. MD simulations with biasing forces resulted in the formation of pores in the bilayer for the derivatives and not for kaempferol. The stability of the pores was assessed by pore closure timescales in unbiased MD simulations, which was found to be 5.3 and 17.0 ns for 2 and 3, respectively. Free energy change for the permeation into the bilayer for kaempferol (1), tertiary amine derivative (2), and arginine derivative (3) was calculated to be -1.5, -48.2, and -100.3 kJ/mol, respectively, which correlate with their antibacterial activity. Furthermore, our results indicate that compound 3 forms a stable toroidal pore in the membrane when multiple molecules are oriented in a transmembrane configuration. Our work sheds light on the mechanism of action of small molecule antimicrobial agents, which can be exploited for the rational design of drug candidates.

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.

Rational design of hyperstable antibacterial peptides for food preservation

We describe the design of peptides with properties like thermostability, pH stability, and antibacterial activity against a few bacterial food pathogens. Insights obtained from classical structure-function analysis of natural peptides and their mutants through antimicrobial and enzymatic assays are used to rationally develop a set of peptides. pH and thermostability assays were performed to demonstrate robust antimicrobial activity post-treatment with high temperatures and at wide pH ranges. We have also investigated the mode of action of these hyperstable peptides using membrane permeability assays, electron microscopy, and molecular dynamics simulations. Notably, through mutational studies, we show that these peptides elicit their antibacterial action via both membrane destabilization and inhibition of intracellular trypsin-the two functions attributable to separate peptide segments. Finally, toxicity studies and food preservation assays demonstrate the safety and efficacy of the designed peptides for food preservation. Overall, the study provides a general 'blueprint' for the development of stable antimicrobial peptides (AMPs). Insights obtained from this work may also be combined with combinatorial methods in high-throughput studies for future development of antimicrobials for various applications.

Computational Methods and Tools in Antimicrobial Peptide Research

The evolution of antibiotic-resistant bacteria is an ongoing and troubling development that has increased the number of diseases and infections that risk going untreated. There is an urgent need to develop alternative strategies and treatments to address this issue. One class of molecules that is attracting significant interest is that of antimicrobial peptides (AMPs). Their design and development has been aided considerably by the applications of molecular models, and we review these here. These methods include the use of tools to explore the relationships between their structures, dynamics, and functions and the increasing application of machine learning and molecular dynamics simulations. This review compiles resources such as AMP databases, AMP-related web servers, and commonly used techniques, together aimed at aiding researchers in the area toward complementing experimental studies with computational approaches.

Molecular dynamics simulations informed by membrane lipidomics reveal the structure-interaction relationship of polymyxins with the lipid A-based outer membrane of Acinetobacter baumannii

BACKGROUND

MDR bacteria represent an urgent threat to human health globally. Polymyxins are a last-line therapy against life-threatening Gram-negative 'superbugs', including Acinetobacter baumannii. Polymyxins exert antimicrobial activity primarily via permeabilizing the bacterial outer membrane (OM); however, the mechanism of interaction between polymyxins and the OM remains unclear at the atomic level.

METHODS

We constructed a lipid A-based OM model of A. baumannii using quantitative membrane lipidomics data and employed all-atom molecular dynamics simulations with umbrella sampling techniques to elucidate the structure-interaction relationship and thermodynamics governing the penetration of polymyxins [B1 and E1 (i.e. colistin A) representing the two clinically used polymyxins] into the OM.

RESULTS

Polymyxin B1 and colistin A bound to the A. baumannii OM by the initial electrostatic interactions between the Dab residues of polymyxins and the phosphates of lipid A, competitively displacing the cations from the headgroup region of the OM. Both polymyxin B1 and colistin A formed a unique folded conformation upon approaching the hydrophobic centre of the OM, consistent with previous experimental observations. Polymyxin penetration induced reorientation of the headgroups of the OM lipids near the penetration site and caused local membrane disorganization, thereby significantly increasing membrane permeability and promoting the subsequent penetration of polymyxin molecules into the OM and periplasmic space.

CONCLUSIONS

The thermodynamics governing the penetration of polymyxins through the outer leaflet of the A. baumannii OM were examined and novel structure-interaction relationship information was obtained at the atomic and membrane level. Our findings will facilitate the discovery of novel polymyxins against MDR Gram-negative pathogens.

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.

Translocation of the nonlabeled antimicrobial peptide PGLa across lipid bilayers and its entry into vesicle lumens without pore formation

Fluorescent-probe-labeled peptides are used to study the interactions of peptides with cells and lipid vesicles but labeling peptides with fluorescent probes can significantly change these interactions. We recently developed a new method to detect the entry of nonlabeled peptides into the lumen of single giant unilamellar vesicles (GUVs). Here we applied this method to examine the interaction of the antimicrobial peptide PGLa with single GUVs to elucidate whether PGLa can enter the GUV lumen without pore formation. First, we examined the interaction of nonlabeled PGLa with single GUVs comprising dioleoylphosphatidylglycerol (DOPG) and dioleoylphosphatidylcholine (DOPC) (4/6) whose lumens contain the fluorescent probe AF647 and DOPG/DOPC (8/2)-large unilamellar vesicles encapsulating a high concentration of calcein. After a large lag period from starting the interaction with PGLa, the fluorescence intensity of the GUV lumen due to calcein (I_calcein) increased gradually without leakage of AF647, indicating that PGLa enters the GUV lumen without pore formation in the GUV membrane. The fraction of entry of PGLa increased with increasing PGLa concentration. Simultaneous measurement of the fractional area change of the GUV membrane (δ) and PGLa-induced increase in I_calcein showed that the entry of PGLa occurs only during the second increase in δ, indicating that PGLa enters the lumen during its translocation from the outer leaflet to the inner leaflet. The fraction of entry of PGLa without pore formation increased with increasing membrane tension. Based on these results, we discuss the elementary processes and the mechanism of the entry of PGLa into the GUV lumen.

Molecular Dynamics Simulations of Human Beta-Defensin Type 3 Crossing Different Lipid Bilayers

Human β defensin type 3 (hBD-3) is a small cationic cysteine-rich peptide. It has a broad spectrum of antimicrobial activities. However, at high concentrations, it also shows hemolytic activity by interrupting red blood cells. To understand the selectivity of hBD-3 disrupting cell membranes, investigating the capability of hBD-3 translocating through different membranes is important. Since hBD-3 in the analogue form in which all three pairs of disulfide bonds are broken has similar antibacterial activities to the wild-type, this project investigates the structure and dynamics of an hBD-3 analogue in monomer, dimer, and tetramer forms through both zwitterionic and negatively charged lipid bilayers using molecular dynamics (MD) simulations. One tetramer structure of hBD-3 was predicted by running all-atom MD simulations on hBD-3 in water at a high concentration, which was found to be stable in water during 400 ns all-atom simulations based on root-mean-squared deviation, root-mean-squared fluctuation, buried surface area, and binding interaction energy calculations. After that, hBD-3 in different forms was placed inside different membranes, and then steered MD simulation was conducted to pull the hBD-3 out of the membrane along the z-direction to generate different configurational windows to set up umbrella-sampling (US) simulations. Because extensive sampling is important to obtain accurate free energy barriers, coarse-grained US MD simulations were performed in each window. Based on the long-term simulation result, membrane thinning was found near hBD-3 in different lipid bilayers and in different hBD-3 oligomer systems. By calculating the root-mean-squared deviation of the z-coordinate of hBD-3 molecules, rotation of the oligomer inside the bilayer and stretching of the oligomer structure along the z-direction were observed. Although reorientation of lipid heads toward the hBD-3 tetramer was observed based on the density profile calculation, the order parameter calculation shows that hBD-3 disrupts 1-palmitoyl-2-oleoyl-sn-glycero-3-phospho-l-serine (POPS) lipids more significantly and makes it less ordered than on 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) lipids. Calculating the free energy of hBD-3 through different lipid bilayers, it was found that generally hBD-3 encounters a lower energy barrier through negatively charged lipid membranes than the zwitterionic membrane. hBD-3 in different forms needs to overcome a lower energy barrier crossing the combined POPC+POPS bilayer through the POPS leaflet than through the POPC leaflet. Besides that, the potential of mean force result suggests that hBD-3 forms an oligomer translocating negatively charged lipid membranes at a low concentration. This study supplied new insight into the antibacterial mechanism of hBD-3 through different membranes.

Rational design of antimicrobial peptides targeting Gram-negative bacteria

Membrane-targeting host antimicrobial peptides (AMPs) can kill or inhibit the growth of Gram-negative bacteria. However, the evolution of resistance among microbes poses a substantial barrier to the long-term utility of the host AMPs. Combining experiment and molecular dynamics simulations, we show that terminal carboxyl capping enhances both membrane insertion and antibacterial activity of an AMP called P1. Furthermore, we show that a bacterial strain with evolved resistance to this peptide becomes susceptible to P1 variants with either backbone capping or lysine-to-arginine substitutions. Our results suggest that cocktails of closely related AMPs may be useful in overcoming evolved resistance.

Molecular dynamics simulations of bacterial outer membrane lipid extraction: Adequate sampling?

The outer membrane of Gram-negative bacteria is almost exclusively composed of lipopolysaccharide in its outer leaflet, whereas the inner leaflet contains a mixture of phospholipids. Lipopolysaccharide diffuses at least an order of magnitude slower than phospholipids, which can cause issues for molecular dynamics simulations in terms of adequate sampling. Here, we test a number of simulation protocols for their ability to achieve convergence with reasonable computational effort using the MARTINI coarse-grained force-field. This is tested in the context both of potential of mean force (PMF) calculations for lipid extraction from membranes and of lateral mixing within the membrane phase. We find that decoupling the cations that cross-link the lipopolysaccharide headgroups from the extracted lipid during PMF calculations is the best approach to achieve convergence comparable to that for phospholipid extraction. We also show that lateral lipopolysaccharide mixing/sorting is very slow and not readily addressable even with Hamiltonian replica exchange. We discuss why more sorting may be unrealistic for the short (microseconds) timescales we simulate and provide an outlook for future studies of lipopolysaccharide-containing membranes.

Tuning of a Membrane-Perforating Antimicrobial Peptide to Selectively Target Membranes of Different Lipid Composition

The use of designed antimicrobial peptides as drugs has been impeded by the absence of simple sequence-structure-function relationships and design rules. The likely cause is that many of these peptides permeabilize membranes via highly disordered, heterogeneous mechanisms, forming aggregates without well-defined tertiary or secondary structure. We suggest that the combination of high-throughput library screening with atomistic computer simulations can successfully address this challenge by tuning a previously developed general pore-forming peptide into a selective pore-former for different lipid types. A library of 2916 peptides was designed based on the LDKA template. The library peptides were synthesized and screened using a high-throughput orthogonal vesicle leakage assay. Dyes of different sizes were entrapped inside vesicles with varying lipid composition to simultaneously screen for both pore size and affinity for negatively charged and neutral lipid membranes. From this screen, nine different LDKA variants that have unique activity were selected, sequenced, synthesized, and characterized. Despite the minor sequence changes, each of these peptides has unique functional properties, forming either small or large pores and being selective for either neutral or anionic lipid bilayers. Long-scale, unbiased atomistic molecular dynamics (MD) simulations directly reveal that rather than rigid, well-defined pores, these peptides can form a large repertoire of functional dynamic and heterogeneous aggregates, strongly affected by single mutations. Predicting the propensity to aggregate and assemble in a given environment from sequence alone holds the key to functional prediction of membrane permeabilization.

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.

Free energy of hydrophilic and hydrophobic pores in lipid bilayers by free energy perturbation of a restraint

The free energy of pore formation in lipid bilayers has been previously calculated using a variety of reaction coordinates. Here, we use free energy perturbation of a cylindrical lipid exclusion restraint to compute the free energy profile as a function of pore radius in dimyristoylphosphatidylcholine (DMPC) and dioleoylphosphatidylcholine (DOPC) bilayers. Additionally restraining the headgroups to lie on the membrane surface allows us to also calculate the free energy profile of hydrophobic pores, i.e., cylindrical pores lined by acyl chains. For certain pore radii, the free energy of wetting of hydrophobic pores is calculated using the density bias method. It is found that wetting of hydrophobic pores becomes thermodynamically favorable at 5.0 Å for DMPC and 6.5 Å for DOPC, although significant barriers prevent spontaneous wetting of the latter on a nanosecond time scale. The free energy of transformation of hydrophilic pores to hydrophobic ones is also calculated using free energy perturbation of headgroup restraints along the bilayer normal. This quantity, along with wetting and pore growth free energies, provides complete free energy profiles as a function of radius. Pore line tension values for the hydrophilic pores obtained from the slope of the free energy profiles are 37.6 pN for DMPC and 53.7 pN for DOPC. The free energy profiles for the hydrophobic pores are analyzed in terms of elementary interfacial tensions. It is found that a positive three-phase line tension is required to explain the results. The estimated value for this three-phase line tension (51.2 pN) lies within the expected range.

Mode-of-Action of Antimicrobial Peptides: Membrane Disruption vs. Intracellular Mechanisms

Antimicrobial peptides are an attractive alternative to traditional antibiotics, due to their physicochemical properties, activity toward a broad spectrum of bacteria, and mode-of-actions distinct from those used by current antibiotics. In general, antimicrobial peptides kill bacteria by either disrupting their membrane, or by entering inside bacterial cells to interact with intracellular components. Characterization of their mode-of-action is essential to improve their activity, avoid resistance in bacterial pathogens, and accelerate their use as therapeutics. Here we review experimental biophysical tools that can be employed with model membranes and bacterial cells to characterize the mode-of-action of antimicrobial peptides.

Twisting and tilting of a mechanosensitive molecular probe detects order in membranes

Lateral forces in biological membranes affect a variety of dynamic cellular processes. Recent synthetic efforts have introduced fluorescent "flippers" as environment-sensitive planarizable push-pull probes that can detect lipid packing and membrane tension, and respond to lipid-induced mechanical forces by a shift in their spectroscopic properties. Herein, we investigate the molecular origin of the mechanosensitivity of the best known flipper, Flipper-TR, by an extended set of molecular dynamics (MD) simulations in membranes of increasing complexity and under different physicochemical conditions, revealing unprecedented details of the sensing process. Simulations enabled by accurate refinement of Flipper-TR force field using quantum mechanical calculations allowed us to unambiguously correlate the planarization of the two fluorescent flippers to spectroscopic response. In particular, Flipper-TR conformation exhibits bimodal distribution in disordered membranes and a unimodal distribution in highly ordered membranes. Such dramatic change was associated with a shift in Flipper-TR excitation spectra, as supported both by our simulated and experimentally-measured spectra. Flipper-TR sensitivity to phase-transition is confirmed by a temperature-jump protocol that alters the lipid phase of an ordered membrane, triggering an instantaneous mechanical twisting of the probe. Simulations show that the probe is also sensitive to surface tension, since even in a naturally disordered membrane, the unimodal distribution of coplanar flippers can be achieved if a sufficiently negative surface tension is applied to the membrane. MD simulations in ternary mixtures containing raft-like nanodomains show that the probe can discriminate lipid domains in phase-separated complex bilayers. A histogram-based approach, called DOB-phase classification, is introduced that can differentiate regions of disordered and ordered lipid phases by comparing dihedral distributions of Flipper-TR. Moreover, a new sensing mechanism involving the orientation of Flipper-TR is elucidated, corroborating experimental evidence that the probe tilt angle is strongly dependent on lipid ordering. The obtained atomic-resolution description of Flipper-TR mechanosensitivity is key to the interpretation of experimental data and to the design of novel mechanosensors with improved spectroscopic properties.

What Makes a Good Pore Former: A Study of Synthetic Melittin Derivatives

Pore formation by membrane-active peptides, naturally encountered in innate immunity and infection, could have important medical and technological applications. Recently, the well-studied lytic peptide melittin has formed the basis for the development of combinatorial libraries from which potent pore-forming peptides have been derived, optimized to work under different conditions. We investigate three such peptides, macrolittin70, which is most active at neutral pH; pHD15, which is active only at low pH; and MelP5_Δ6, which was rationally designed to be active at low pH but formed only small pores. There are three, six, and six acidic residues in macrolittin70, pHD15, and MelP5_Δ6, respectively. We perform multi-microsecond simulations in 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) of hexamers of these peptides starting from transmembrane orientations at neutral pH (all residues at standard protonation), low pH (acidic residues and His protonated), and highly acidic environments in which C-termini are also protonated. Previous simulations of the parent peptides melittin and MelP5 in 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) are repeated in POPC. We find that the most potent pore-forming peptides exhibit strong interpeptide interactions, including salt bridges, H-bonds, and polar interactions. Protonation of the C-terminus promotes helicity and pore size. The proximity of the peptides allows fewer lipid headgroups to line the pores than in previous simulations, making the pores intermediate between barrel stave and toroidal. Based on these structures and geometrical arguments, we attempt to rationalize the factors that under different conditions can increase or decrease pore stability and propose mutations that could be tested experimentally.

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

Antimicrobial peptide ROAD-1 triggers phase change in local membrane environment to execute its activity

Emergence of antibiotic-resistant pathogens has paved way for development of newer class of drugs that would not be susceptible to resistance. Antimicrobial peptides such as defensins that target the microbial membrane are promising candidates. ROAD-1 is an alpha-defensin present in the oral cavity of rhesus macaque and shares very high sequence similarity to human enteric defensin 5. In this study we have performed microsecond long all atom molecular dynamic simulations to understand the mechanism of action of ROAD-1. We find that ROAD-1 is able to adopt an energetically stable conformation predominantly stabilized by electrostatic interactions only in presence of bacterial membranes. In mammalian membrane even though it gets absorbed onto the bilayer, it is unable to adopt an equilibrium conformation. Binding of ROAD-1 to bilayer induces clustering of POPG molecules up to 15 Å around the peptide. POPG molecules show higher order parameters than the neighboring POPE implying coexistence of different phases. Analysis of binding free energy of ROAD-1-membrane complex indicates Arg1, Arg2, Arg7, and Arg25 to play key role in its antimicrobial activity. Unlike its homolog HD5, ROAD-1 is not observed to form a dimer. Our study gives insight into the membrane-bound conformation of ROAD-1 and its mechanism of action that can aid in designing defensin-based therapeutics. Graphical abstract Antimicrobial peptide ROAD-1 adopts a different membrane-bound conformation as compared with HD5 even though they belong to the same family implying a different mechanism of action.

Accelerating Membrane Simulations with Hydrogen Mass Repartitioning

The time step of atomistic molecular dynamics (MD) simulations is determined by the fastest motions in the system and is typically limited to 2 fs. An increasingly popular approach is to increase the mass of the hydrogen atoms to ∼3 amu and decrease the mass of the parent atom by an equivalent amount. This approach, known as hydrogen-mass repartitioning (HMR), permits time steps up to 4 fs with reasonable simulation stability. While HMR has been applied in many published studies to date, it has not been extensively tested for membrane-containing systems. Here, we compare the results of simulations of a variety of membranes and membrane-protein systems run using a 2 fs time step and a 4 fs time step with HMR. For pure membrane systems, we find almost no difference in structural properties, such as area-per-lipid, electron density profiles, and order parameters, although there are differences in kinetic properties such as the diffusion constant. Conductance through a porin in an applied field, partitioning of a small peptide, hydrogen-bond dynamics, and membrane mixing show very little dependence on HMR and the time step. We also tested a 9 Å cutoff as compared to the standard CHARMM cutoff of 12 Å, finding significant deviations in many properties tested. We conclude that HMR is a valid approach for membrane systems, but a 9 Å cutoff is not.

Folding profiles of antimicrobial scorpion venom-derived peptides on hydrophobic surfaces: a molecular dynamics study

Most helical antimicrobial peptides (AMPs) are usually unfolded in aqueous solution; however they acquire their secondary structure in the presence of a hydrophobic environment such as lipid membranes. Being the biological membranes the main target of many AMPs it is necessary to understand their way of action. Pandinin 2 (Pin2) is an alpha-helical AMP isolated from the venom of the African scorpion Pandinus imperator which shows high antimicrobial activity against Gram-positive bacteria and it is less active against Gram-negative bacteria, nevertheless, it has strong hemolytic activity. Its chemically synthesized Pin2GVG analog has low hemolytic activity while keeping its antimicrobial activity. With the aim of exploring the partition and subsequent folding of these peptides, in this work we report the results of extensive molecular dynamics simulations of Pin2 and Pin2GVG peptides in the presence of 2 hydrophobic environments such as dodecyl-phosphocholine (DPC) micelle and 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocoline (POPC) membrane. Our results indicate that Pin2 folds in DPC with a 79% of alpha-helical content, which is in agreement with the experimental results, while in POPC it has 62.5% of alpha-helical content. On the other hand, Pin2GVG presents a higher percentage of alpha-helical structure in POPC and a smaller content in DPC when compared with Pin2. These results can help to better choose the starting structures in future molecular dynamics simulations of AMPs, because these peptides can adopt slightly different conformations depending on the hydrophobic environment.Communicated by Ramaswamy H. Sarma.

Substitution of lysine for isoleucine at the center of the nonpolar face of the antimicrobial peptide, piscidin-1, leads to an increase in the rapidity of bactericidal activity and a reduction in toxicity

Purpose: Piscidin-1 is an effective antimicrobial peptide (AMP) against a variety of microbes. However, its toxicity has been reported as a limitation for its potential therapeutic applications. The toxicity of piscidin-1 may be related to the long nonpolar face of this AMP. Here, we investigated different piscidin-1 analogs to reach a peptide with the reduced toxicity.

Material and methods: In vitro and in vivo antibacterial activity and toxicity of piscidin-1 analogs generated by replacement of isoleucine at the border (I9) or the center (I16) of the nonpolar face of piscidin-1 by alanine or lysine were investigated.

Results: The results indicated that among all peptides, piscidin-1 with the highest HPLC retention time (RT) and I16K-piscidin-1 with the lowest RT had the highest and lowest cytotoxicity, respectively. Although I16K-piscidin-1 possessed the same MIC value as the parent peptide (piscidin-1) and other analogs, I16K-piscidin-1 exhibited a higher rapidity of bactericidal action at 5×MIC. The β-galactosidase leakage and propidium iodide staining assays indicated a higher pore-forming capacity of I16K-piscidin-1 relative to the parent peptide (piscidin-1). Taken together, RT is suggested to have a direct association with the toxicity and an inverse association with the rapidity of bactericidal action and pore-forming capacity. After infection of mice with clinical colistin-resistant Acinetobacter baumannii or clinical methicillin-resistant Staphylococcus aureus strains, treatment with I16K-piscidin-1, but not piscidin-1 and other analogs, resulted in a significantly stronger bactericidal potency. Furthermore, I16K-piscidin-1 exhibited the lowest in vivo toxicity. 

Conclusion: Overall, in vitro and in vivo comparison of piscidin-1 and its analogs together documented that replacement of isoleucine at the center of the nonpolar face of piscidin-1(I16) by lysine leads to not only a decrease in toxicity potential but also an increase in bactericidal potential.

Activity and characterization of a pH-sensitive antimicrobial peptide

Antimicrobial peptides (AMPs) have been an area of great interest, due to the high selectivity of these molecules toward bacterial targets over host cells and the limited development of bacterial resistance to these molecules throughout evolution. Previous work showed that when Histidine was incorporated into the peptide C18G it lost antimicrobial activity. The role of pH on activity and biophysical properties of the peptide was investigated to explain this phenomenon. Minimal inhibitory concentration (MIC) results demonstrated that decreased media pH increased antimicrobial activity. Trichloroethanol (TCE) quenching and red-edge excitation spectroscopy (REES) showed a clear pH dependence on peptide aggregation in solution. Trp fluorescence was used to monitor binding to lipid vesicles and demonstrated the peptide binds to anionic bilayers at all pH values tested, however, binding to zwitterionic bilayers was enhanced at pH 7 and 8 (above the His pKa). Dual Quencher Analysis (DQA) confirmed the peptide inserted more deeply in PC:PG and PE:PG membranes, but could insert into PC bilayers at pH conditions above the His pKa. Bacterial membrane permeabilization assays which showed enhanced membrane permeabilization at pH 5 and 6 but vesicle leakage assays indicate enhanced permeabilization of PC and PC:PG bilayers at neutral pH. The results indicate the ionization of the His side chain affects the aggregation state of the peptide in solution and the conformation the peptide adopts when bound to bilayers, but there are likely more subtle influences of lipid composition and properties that impact the ability of the peptide to form pores in membranes.

Characterization of Lipid-Protein Interactions and Lipid-Mediated Modulation of Membrane Protein Function through Molecular Simulation

The cellular membrane constitutes one of the most fundamental compartments of a living cell, where key processes such as selective transport of material and exchange of information between the cell and its environment are mediated by proteins that are closely associated with the membrane. The heterogeneity of lipid composition of biological membranes and the effect of lipid molecules on the structure, dynamics, and function of membrane proteins are now widely recognized. Characterization of these functionally important lipid-protein interactions with experimental techniques is however still prohibitively challenging. Molecular dynamics (MD) simulations offer a powerful complementary approach with sufficient temporal and spatial resolutions to gain atomic-level structural information and energetics on lipid-protein interactions. In this review, we aim to provide a broad survey of MD simulations focusing on exploring lipid-protein interactions and characterizing lipid-modulated protein structure and dynamics that have been successful in providing novel insight into the mechanism of membrane protein function.

Understanding the Role of Lipids in Signaling Through Atomistic and Multiscale Simulations of Cell Membranes

Cell signaling controls essentially all cellular processes. While it is often assumed that proteins are the key architects coordinating cell signaling, recent studies have shown more and more clearly that lipids are also involved in signaling processes in a number of ways. Lipids do, for instance, act as messengers, modulate membrane receptor conformation and dynamics, and control membrane receptor partitioning. Further, through structural modifications such as oxidation, the functions of lipids as part of signaling processes can be modified. In this context, in this article we discuss the understanding recently revealed by atomistic and coarse-grained computer simulations of nanoscale processes and underlying physicochemical principles related to lipids’ functions in cellular signaling.

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.

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.

Different effects of cholesterol on membrane permeation of arginine and tryptophan revealed by bias-exchange metadynamics simulations

Experiments have shown that cholesterol influences the membrane permeability of small molecules, amino acids, and cell-penetrating peptides. However, their exact translocation mechanisms under the influence of cholesterol remain poorly understood. Given the practical importance of cell-penetrating peptides and the existence of varied cholesterol contents in different cell types, it is necessary to examine the permeation of amino acids in cholesterol-containing membranes at atomic level of details. Here, bias-exchange metadynamics simulations were employed to investigate the molecular mechanism of the membrane permeation of two amino acids Arg and Trp important for cell-penetrating peptides in the presence of different concentrations of cholesterol. We found that the free energy barrier of Arg+ (the protonated form) permeation increased linearly as the cholesterol concentration increased, whereas the barrier of Trp permeation had a rapid increase from 0 mol. % to 20 mol. % cholesterol-containing membranes and nearly unchanged from 20 mol. % to 40 mol. % cholesterol-containing membranes. Arg0 becomes slightly more stable than Arg+ at the center of the dipalmitoylphosphatidylcholine (DPPC) membrane with 40 mol. % cholesterol concentrations. As a result, Arg+ has a similar permeability as Trp at 0 mol. % and 20 mol. % cholesterol, but a significantly lower permeability than Trp at 40 mol. % cholesterol. This difference is caused by the gradual reduction of water defects for Arg+ as the cholesterol concentration increases but lack of water defects for Trp in cholesterol-containing membranes. Strong but different orientation dependence between Arg+ and Trp permeations is observed. These results provide an improved microscopic understanding of amino-acid permeation through cholesterol-containing DPPC membrane systems.

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.

Sugar-based bactericides targeting phosphatidylethanolamine-enriched membranes

Anthrax is an infectious disease caused by Bacillus anthracis, a bioterrorism agent that develops resistance to clinically used antibiotics. Therefore, alternative mechanisms of action remain a challenge. Herein, we disclose deoxy glycosides responsible for specific carbohydrate-phospholipid interactions, causing phosphatidylethanolamine lamellar-to-inverted hexagonal phase transition and acting over B. anthracis and Bacillus cereus as potent and selective bactericides. Biological studies of the synthesized compound series differing in the anomeric atom, glycone configuration and deoxygenation pattern show that the latter is indeed a key modulator of efficacy and selectivity. Biomolecular simulations show no tendency to pore formation, whereas differential metabolomics and genomics rule out proteins as targets. Complete bacteria cell death in 10 min and cellular envelope disruption corroborate an effect over lipid polymorphism. Biophysical approaches show monolayer and bilayer reorganization with fast and high permeabilizing activity toward phosphatidylethanolamine membranes. Absence of bacterial resistance further supports this mechanism, triggering innovation on membrane-targeting antimicrobials.

Bacillus anthracis causes the infectious disease anthrax. Here, the authors synthesized deoxy glycosides that are effective against B. anthracis and related bacteria and found that these amphiphilic compounds kill bacteria via an unusual mechanism of action.

Theoretical Insights into the Interactions between Star-Shaped Antimicrobial Polypeptides and Bacterial Membranes

A structurally nanoengineered antimicrobial polypeptide consisting of lysine and valine residues is a new class of antimicrobial agent with superior antibacterial activity against multidrug-resistant bacteria and low toxicity toward mammalian cells. Utilizing coarse-grained models, we studied the interactions of microbial cytoplasmic membranes with polypeptides of either (K₂V₁)₅ (star-KV) or CM15 (star-CM15). Our computational results verify the low toxicity of polypeptides of (K₂V₁)₅ toward the dipalmitoyl phosphatidylcholine bilayer. This low toxicity is demonstrated to originate from weakened hydrophobicity combined with its random coil conformation for (K₂V₁)₅ because of the highly abundant valine residues, compared with the typical antimicrobial peptides, such as CM15. In the interactions with a palmitoyl-oleoyl-phosphatidylethanolamine/palmitoyl-oleoyl-phosphatidylglycerol bilayer, star-KV has greater ability in phase separation and generation of phase boundary defects not only in lipid redistribution but also in lateral dynamic movements, although both star-KV and star-CM15 can extract the phosphatidylglycerol lipids and purify the phosphatidylethanolamine lipids into continuum domains. We suggest that the polypeptide of (K₂V₁)₅ can nondisruptively kill bacteria by hampering bacterial metabolism through reorganizing lipid domain distribution and simultaneously "freezing" lipid movement.

Characterizing the Molecular Mechanisms for Flipping Charged Peptide Flanking Loops across a Lipid Bilayer

The cell membrane largely prevents the passive diffusion of charged molecules due to the large free energy barrier associated with translocating charged groups across the hydrophobic lipid bilayer core. Despite this barrier, some peptides can interconvert between transmembrane and surface-adsorbed states by "flipping" charged flanking loops across the bilayer on a surprisingly rapid second-minute time scale. The transmembrane helices of some multispanning membrane proteins undergo similar reorientation processes, suggesting that loop-flipping may be a mechanism for regulating membrane protein topology; however, the molecular mechanisms underlying this behavior remain unknown. In this work, we study the loop-flipping behavior exhibited by a peptide with a hydrophobic transmembrane helix, charged flanking loops, and a central, membrane-exposed aspartate residue of varying protonation state. We utilize all-atom temperature accelerated molecular dynamics simulations to predict the likelihood of loop-flipping without predefining specific loop-flipping pathways. We demonstrate that this approach can identify multiple possible flipping pathways, with the prevalence of each pathway depending on the protonation state of the central residue. In particular, we find that a charged central residue facilitates loop-flipping by stabilizing membrane water defects, enabling the "self-catalysis" of charge translocation. These findings provide detailed molecular-level insights into charged loop-flipping pathways that may generalize to other charge translocation processes, such as lipid flip-flop or the large-scale conformational rearrangements of multispanning membrane proteins.