Membrane models of E. coli containing cyclic moieties in the aliphatic lipid chain

Discovery of Distinctin-Like-Peptide-PH (DLP-PH) From the Skin Secretion of Phyllomedusa hypochondrialis, a Prototype of a Novel Family of Antimicrobial Peptide

Amphibian skin secretions are an important treasure house of bioactive antimicrobial peptides (AMPs). Despite having been the focus of decades of research in this context, investigations of phyllomedusine frogs continue to identify new AMPs from their skin secretions. In this study, the prototype of a novel family of AMP distinctin-like-peptide-PH (DLP-PH) was identified from the skin secretion of the otherwise well-studied Tiger-Legged Tree Frog Phyllomedusa hypochondrialis through cloning of its precursor-encoding cDNA from a skin secretion-derived cDNA library by a 3′-rapid amplification of cDNA ends (RACE) strategy. Subsequently, the mature peptide was isolated and characterized using reverse-phase HPLC and MS/MS fragmentation sequencing. DLP-PH adopted an α-helical conformation in membrane mimetic solution and demonstrated unique structural features with two distinct domains that differed markedly in their physiochemical properties. Chemically synthesized replicates of DLP-PH showed antimicrobial activity against planktonic bacterial and yeast cells, but more potent against Escherichia coli at 32 μg/mL. However, DLP-PH showed much weaker inhibitory activity against the growth of sessile cells in biofilms. In addition, DLP-PH exhibited anti-proliferative activity against human cancer cell lines, H157, and PC3, but with no major toxicity against normal human cell, HMEC-1. These combined properties make DLP-PH deserving further study as an antimicrobial agent and further investigations of its structure-activity relationship could provide valuable new insights into drug lead candidates for antimicrobial and/or anti-cancer purposes.

Unravelling a Mechanism of Action for a Cecropin A-Melittin Hybrid Antimicrobial Peptide: The Induced Formation of Multilamellar Lipid Stacks

An understanding of the mechanism of action of antimicrobial peptides is fundamental to the development of new and more active antibiotics. In the present work, we use a wide range of techniques (SANS, SAXD, DSC, ITC, CD, and confocal and electron microscopy) in order to fully characterize the interaction of a cecropin A-melittin hybrid antimicrobial peptide, CA(1-7)M(2-9), of known antimicrobial activity, with a bacterial model membrane of POPE/POPG in an effort to unravel its mechanism of action. We found that CA(1-7)M(2-9) disrupts the vesicles, inducing membrane condensation and forming an onionlike structure of multilamellar stacks, held together by the intercalated peptides. SANS and SAXD revealed changes induced by the peptide in the lipid bilayer thickness and the bilayer stiffening in a tightly packed liquid-crystalline lamellar phase. The analysis of the observed abrupt changes in the repeat distance upon the phase transition to the gel state suggests the formation of an L_γ phase. To the extent of our knowledge, this is the first time that the L_γ phase is identified as part of the mechanism of action of antimicrobial peptides. The energetics of interaction depends on temperature, and ITC results indicate that CA(1-7)M(2-9) interacts with the outer leaflet. This further supports the idea of a surface interaction that leads to membrane condensation and not to pore formation. As a result, we propose that this peptide exerts its antimicrobial action against bacteria through extensive membrane disruption that leads to cell death.

Simulations of simple Bovine and Homo sapiens outer cortex ocular lens membrane models with a majority concentration of cholesterol

The lipid composition of bovine and human ocular lens membranes has been probed, and a variety of lipids have been found including phosphatidylcholine (PC), phosphatidylethanolamine (PE), sphingomyelin (SM), and cholesterol (CHOL) with cholesterol being present in particularly high concentrations. In this study, we use the all-atom CHARMM36 force field to simulate binary, ternary, and quaternary mixtures as models of the ocular lens. High concentration of cholesterol, in combination with different and varying diversity of phospholipids (PL) and sphingolipids (SL), affect the structure of the ocular lens lipid bilayer. The following analyses were done for each simulation: surface area per lipid, component surface area per lipid, deuterium order parameters (S_CD), electron density profiles (EDP), membrane thickness, hydrogen bonding, radial distribution functions, clustering, and sterol tilt angle distribution. The S_CD show significant bilayer alignment and packing in cholesterol-rich bilayers. The EDP show the transition from liquid crystalline to liquid ordered with the addition of cholesterol. Hydrogen bonds in our systems show the tendency for intramolecular interactions between cholesterol and fully saturated lipid tails for less complex bilayers. But with an increased number of components in the bilayer, the acyl chain of the lipids becomes a less important characteristic, and the headgroup of the lipid becomes more significant. Overall, cholesterol is the driving force of membrane structure of the ocular lens membrane where interactions between cholesterol, PL, and SL determine structure and function of the biomembrane. The goal of this work is to develop a baseline for further study of more physiologically realistic ocular lens lipid membranes. This article is part of a Special Issue entitled: Emergence of Complex Behavior in Biomembranes edited by Marjorie Longo.

Simulations of simple linoleic acid-containing lipid membranes and models for the soybean plasma membranes

The all-atom CHARMM36 lipid force field (C36FF) has been tested with saturated, monounsaturated, and polyunsaturated lipids; however, it has not been validated against the 18:2 linoleoyl lipids with an unsaturated sn-1 chain. The linoleoyl lipids are common in plants and the main component of the soybean membrane. The lipid composition of soybean plasma membranes has been thoroughly characterized with experimental studies. However, there is comparatively less work done with computational modeling. Our molecular dynamics (MD) simulation results show that the pure linoleoyl lipids, 1-stearoyl-2-linoleoyl-sn-glycero-3-phosphocholine (18:0/18:2) and 1,2-dilinoleoyl-sn-glycero-3-phosphocholine (di-18:2), agree very well with the experiments, which demonstrates the accuracy of the C36FF for the computational study of soybean membranes. Based on the experimental composition, the soybean hypocotyl and root plasma membrane models are developed with each containing seven or eight types of linoleoyl phospholipids and two types of sterols (sitosterol and stigmasterol). MD simulations are performed to characterize soybean membranes, and the hydrogen bonds and clustering results demonstrate that the lipids prefer to interact with the lipids of the same/similar tail unsaturation. All the results suggest that these two soybean membrane models can be used as a basis for further research in soybean and higher plant membranes involving membrane-associated proteins.

An extensive simulation study of lipid bilayer properties with different head groups, acyl chain lengths, and chain saturations

Previous MD simulations of six phosphocholine (PC) lipid bilayers demonstrated the accuracy of the CHARMM36 force field (C36FF) for PC bilayer simulation at varied temperatures (BBA-Biomembranes, 1838 (2014): 2520-2529). In this work, we further examine the accuracy of C36FF over a wide temperature range for a broader range of lipid types such as various head groups (phosphatidic acid (PA), PC, phosphoethanolamine (PE), phosphoglycerol (PG), and phosphoserine (PS)), and tails (saturated, mono-, mixed- and poly-unsaturated acyl chains with varied chain lengths). The structural properties (surface area per lipid (SA/lip), overall bilayer thickness, hydrophobic thickness, headgroup-to-headgroup thickness, deuterium order parameter (S_CD), and spin-lattice relaxation time (T₁)) obtained from simulations agree well with nearly all available experimental data. Our analyses indicate that PS lipids have the most inter-lipid hydrogen bonds, while PG lipids have the most intra-lipid hydrogen bonds, which play the main role in their low SA/lip in PS lipids and low thicknesses in PG lipids, respectively. PS, PE, and PA lipids have the largest contact clusters with on average 5-8 lipids per cluster, while PC and PG have clusters of 4 lipids based on a cutoff distance of 6.5Å. PS lipids have much slower lipid wobble (i.e., higher correlation time) than other head groups at a given temperature as the hydrogen bonded network significantly reduces a lipid's mobility, and the rate of lipid wobble increases dramatically as temperature increases. These in-depth analyses facilitate further understanding of lipid bilayers at the atomic level.

CHARMM-GUI Input Generator for NAMD, GROMACS, AMBER, OpenMM, and CHARMM/OpenMM Simulations Using the CHARMM36 Additive Force Field

Proper treatment of nonbonded interactions is essential for the accuracy of molecular dynamics (MD) simulations, especially in studies of lipid bilayers. The use of the CHARMM36 force field (C36 FF) in different MD simulation programs can result in disagreements with published simulations performed with CHARMM due to differences in the protocols used to treat the long-range and 1-4 nonbonded interactions. In this study, we systematically test the use of the C36 lipid FF in NAMD, GROMACS, AMBER, OpenMM, and CHARMM/OpenMM. A wide range of Lennard-Jones (LJ) cutoff schemes and integrator algorithms were tested to find the optimal simulation protocol to best match bilayer properties of six lipids with varying acyl chain saturation and head groups. MD simulations of a 1,2-dipalmitoyl-sn-phosphatidylcholine (DPPC) bilayer were used to obtain the optimal protocol for each program. MD simulations with all programs were found to reasonably match the DPPC bilayer properties (surface area per lipid, chain order parameters, and area compressibility modulus) obtained using the standard protocol used in CHARMM as well as from experiments. The optimal simulation protocol was then applied to the other five lipid simulations and resulted in excellent agreement between results from most simulation programs as well as with experimental data. AMBER compared least favorably with the expected membrane properties, which appears to be due to its use of the hard-truncation in the LJ potential versus a force-based switching function used to smooth the LJ potential as it approaches the cutoff distance. The optimal simulation protocol for each program has been implemented in CHARMM-GUI. This protocol is expected to be applicable to the remainder of the additive C36 FF including the proteins, nucleic acids, carbohydrates, and small molecules.

Anti-infectives in Drug Delivery-Overcoming the Gram-Negative Bacterial Cell Envelope

Infectious diseases are becoming a major menace to the state of health worldwide, with difficulties in effective treatment especially of nosocomial infections caused by Gram-negative bacteria being increasingly reported. Inadequate permeation of anti-infectives into or across the Gram-negative bacterial cell envelope, due to its intrinsic barrier function as well as barrier enhancement mediated by resistance mechanisms, can be identified as one of the major reasons for insufficient therapeutic effects. Several in vitro, in silico, and in cellulo models are currently employed to increase the knowledge of anti-infective transport processes into or across the bacterial cell envelope; however, all such models exhibit drawbacks or have limitations with respect to the information they are able to provide. Thus, new approaches which allow for more comprehensive characterization of anti-infective permeation processes (and as such, would be usable as screening methods in early drug discovery and development) are desperately needed. Furthermore, delivery methods or technologies capable of enhancing anti-infective permeation into or across the bacterial cell envelope are required. In this respect, particle-based carrier systems have already been shown to provide the opportunity to overcome compound-related difficulties and allow for targeted delivery. In addition, formulations combining efflux pump inhibitors or antimicrobial peptides with anti-infectives show promise in the restoration of antibiotic activity in resistant bacterial strains. Despite considerable progress in this field however, the design of carriers to specifically enhance transport across the bacterial envelope or to target difficult-to-treat (e.g., intracellular) infections remains an urgently needed area of improvement. What follows is a summary and evaluation of the state of the art of both bacterial permeation models and advanced anti-infective formulation strategies, together with an outlook for future directions in these fields.

Formation and Characterization of Supported Lipid Bilayers Composed of Hydrogenated and Deuterated Escherichia coli Lipids

Supported lipid bilayers are widely used for sensing and deciphering biomolecular interactions with model cell membranes. In this paper, we present a method to form supported lipid bilayers from total lipid extracts of Escherichia coli by vesicle fusion. We show the validity of this method for different types of extracts including those from deuterated biomass using a combination of complementary surface sensitive techniques; quartz crystal microbalance, neutron reflection and atomic force microscopy. We find that the head group composition of the deuterated and the hydrogenated lipid extracts is similar (approximately 75% phosphatidylethanolamine, 13% phosphatidylglycerol and 12% cardiolipin) and that both samples can be used to reconstitute high-coverage supported lipid bilayers with a total thickness of 41 ± 3 Å, common for fluid membranes. The formation of supported lipid bilayers composed of natural extracts of Escherichia coli allow for following biomolecular interactions, thus advancing the field towards bacterial-specific membrane biomimics.

LipidBuilder: A Framework To Build Realistic Models for Biological Membranes

The physical and chemical characterization of biological membranes is of fundamental importance for understanding the functional role of lipid bilayers in shaping cells and organelles, steering vesicle trafficking and promoting membrane-protein signaling. Molecular dynamics simulations stand as a powerful tool to probe the properties of membranes at atomistic level. However, the biological membrane is highly complex, and closely mimicking its physiological constitution in silico is not a straightforward task. Here, we present LipidBuilder, a framework for creating and storing models of biologically relevant phospholipid species with acyl tails of heterogeneous composition. LipidBuilder also enables the assembly of these database-stored lipids into realistic bilayers featuring asymmetric distribution on layer leaflets and concentration of given membrane constituents as defined, for example, by lipidomics experiments. The ability of LipidBuilder to assemble robust membrane models was validated by simulating membranes of homogeneous lipid composition for which experimental data are available. Furthermore, taking advantage of the extensive lipid headgroup repertoire, we assembled models of membranes of heterogeneous nature as naturally found in viral (phage PRD1), bacterial (Salmonella enterica, Laurinavicius , S. ; Kakela , R. ; Somerharju , P. ; Bamford , D. H. ; Virology 2004 , 322 , 328 - 336 ) and plant (Chlorella kessleri, Rezanka , T. ; Podojil , M. ; J. Chromatogr. 1989 , 463 , 397 - 408 ) organisms. These realistic membrane models were built using a near-exact lipid composition revealed from analytical chemistry experiments. We suggest LipidBuilder as a useful tool to model biological membranes of near-biological complexity, and as a robust complement to the current efforts to characterize the biophysical properties of biological membrane using molecular simulation.

CHARMM-GUI Input Generator for NAMD, GROMACS, AMBER, OpenMM, and CHARMM/OpenMM Simulations Using the CHARMM36 Additive Force Field

Proper treatment of nonbonded interactions is essential for the accuracy of molecular dynamics (MD) simulations, especially in studies of lipid bilayers. The use of the CHARMM36 force field (C36 FF) in different MD simulation programs can result in disagreements with published simulations performed with CHARMM due to differences in the protocols used to treat the long-range and 1-4 nonbonded interactions. In this study, we systematically test the use of the C36 lipid FF in NAMD, GROMACS, AMBER, OpenMM, and CHARMM/OpenMM. A wide range of Lennard-Jones (LJ) cutoff schemes and integrator algorithms were tested to find the optimal simulation protocol to best match bilayer properties of six lipids with varying acyl chain saturation and head groups. MD simulations of a 1,2-dipalmitoyl-sn-phosphatidylcholine (DPPC) bilayer were used to obtain the optimal protocol for each program. MD simulations with all programs were found to reasonably match the DPPC bilayer properties (surface area per lipid, chain order parameters, and area compressibility modulus) obtained using the standard protocol used in CHARMM as well as from experiments. The optimal simulation protocol was then applied to the other five lipid simulations and resulted in excellent agreement between results from most simulation programs as well as with experimental data. AMBER compared least favorably with the expected membrane properties, which appears to be due to its use of the hard-truncation in the LJ potential versus a force-based switching function used to smooth the LJ potential as it approaches the cutoff distance. The optimal simulation protocol for each program has been implemented in CHARMM-GUI. This protocol is expected to be applicable to the remainder of the additive C36 FF including the proteins, nucleic acids, carbohydrates, and small molecules.

Mechanical properties of lipid bilayers from molecular dynamics simulation

Lipid areas (Aℓ), bilayer area compressibilities (KA), bilayer bending constants (KC), and monolayer spontaneous curvatures (c0) from simulations using the CHARMM36 force field are reported for 12 representative homogenous lipid bilayers. Aℓ (or their surrogate, the average deuterium order parameter in the "plateau region" of the chain) agree very well with experiment, as do the KA. Simulated KC are in near quantitative agreement with vesicle flicker experiments, but are somewhat larger than KC from X-ray, pipette aspiration, and neutron spin echo for saturated lipids. Spontaneous curvatures of bilayer leaflets from the simulations are approximately 30% smaller than experimental values of monolayers in the inverse hexagonal phase.

Probing the importance of lipid diversity in cell membranes via molecular simulation

Lipid membranes in prokaryotes and eukaryotes have a wide array of lipids that are necessary for proper membrane structure and function. In this paper, an introduction to lipid diversity in biology and a mini-review on how molecular simulations have been used to model biological membranes (primarily limited to one to three lipid types in most simulation-based models) is provided, which motivates the use of all-atom molecular dynamics (MD) simulations to study the effect of lipid diversity on properties of realistic membrane models of prokaryotes and eukaryotes. As an example, cytoplasmic membrane models of Escherichia coli were developed at different stages of the colony growth cycle (early-log, mid-log, stationary and overnight). The main difference between lipid compositions at each stage was the concentration of a cyclopropane-containing moiety on the sn-2 lipid acyl chain (cyC17:0). Triplicate MD simulations for each stage were run for 300 ns to study the influence of lipid diversity on the surface area per lipid, area compressibility modulus, deuterium order parameters, and electron density profiles. The overnight stage (also known as the death stage) had the highest average surface area per lipid, highest rigidity, and lowest bilayer thickness compare to other stages of E. coli cytoplasmic membrane. Although bilayer thickness did depend on the growth stage, the changes between these were small suggesting that the hydrophobic core of transmembrane proteins fit well with the membrane in all growth stages. Although it is still common practise in MD simulations of membrane proteins to use simple one- or two-component membranes, it can be important to use diverse lipid model membranes when membrane protein structure and function are influenced by changes in lipid membrane composition.

Studies of the interactions of ursane-type bioactive terpenes with the model of Escherichia coli inner membrane-Langmuir monolayer approach

Pentacyclic triterpenes (PT), ursolic acid (Urs), and α-amyrin (AMalf) are natural products exhibiting broad spectrum of antibacterial activity. These compounds are membrane-active and can disorder bacterial membranes when incorporated; however, the exact mechanism of their membrane activity is unknown. In our studies, we applied Langmuir monolayer technique supported by Brewster angle microscopy to model the interactions of the selected PT with the lipid matrix of E. coli inner membrane. As the model membrane, we applied mixtures (75/25 mole/.mole %) of the representative Escherichia coli phosphatidylethanolamine (POPE), with the cardiolipin (ECCL) or phosphatidylglycerol (ECPG) extracted from the E. coli inner membrane. On the basis of the recorded isotherms, we performed thermodynamic analysis and calculated free energy of mixing ΔGexc. It turned out that the phospholipids forming the inner membrane of E. coli are ideally miscible, whereas in binary systems composed of PT and POPE, negative deviations from ideality indicating attractive interactions between the investigated PT and POPE molecules were observed. On the other hand, in ternary systems composed of PT, POPE and one of the E. coli anionic phospholipids large positive changes in ΔGexc were observed. Thus, both PT exhibit disorganizing effect on the model E. coli membrane. It was also proved that at low terpene proportion, AMalf can be more active than Urs. However, at higher proportion Urs incorporation can lead to the disintegration of cardiolipin-rich domains present in bacterial membrane.

Microsecond molecular dynamics simulations of lipid mixing

Molecular dynamics (MD) simulations of membranes are often hindered by the slow lateral diffusion of lipids and the limited time scale of MD. In order to study the dynamics of mixing and characterize the lateral distribution of lipids in converged mixtures, we report microsecond-long all-atom MD simulations performed on the special-purpose machine Anton. Two types of mixed bilayers, POPE:POPG (3:1) and POPC:cholesterol (2:1), as well as a pure POPC bilayer, were each simulated for up to 2 μs. These simulations show that POPE:POPG and POPC:cholesterol are each fully miscible at the simulated conditions, with the final states of the mixed bilayers similar to a random mixture. By simulating three POPE:POPG bilayers at different NaCl concentrations (0, 0.15, and 1 M), we also examined the effect of salt concentration on lipid mixing. While an increase in NaCl concentration is shown to affect the area per lipid, tail order, and lipid lateral diffusion, the final states of mixing remain unaltered, which is explained by the largely uniform increase in Na(+) ions around POPE and POPG. Direct measurement of water permeation reveals that the POPE:POPG bilayer with 1 M NaCl has reduced water permeability compared with those at zero or low salt concentration. Our calculations provide a benchmark to estimate the convergence time scale of all-atom MD simulations of lipid mixing. Additionally, equilibrated structures of POPE:POPG and POPC:cholesterol, which are frequently used to mimic bacterial and mammalian membranes, respectively, can be used as starting points of simulations involving these membranes.

CHARMM-GUI Membrane Builder toward realistic biological membrane simulations

CHARMM-GUI Membrane Builder, http://www.charmm-gui.org/input/membrane, is a web-based user interface designed to interactively build all-atom protein/membrane or membrane-only systems for molecular dynamics simulations through an automated optimized process. In this work, we describe the new features and major improvements in Membrane Builder that allow users to robustly build realistic biological membrane systems, including (1) addition of new lipid types, such as phosphoinositides, cardiolipin (CL), sphingolipids, bacterial lipids, and ergosterol, yielding more than 180 lipid types, (2) enhanced building procedure for lipid packing around protein, (3) reliable algorithm to detect lipid tail penetration to ring structures and protein surface, (4) distance-based algorithm for faster initial ion displacement, (5) CHARMM inputs for P21 image transformation, and (6) NAMD equilibration and production inputs. The robustness of these new features is illustrated by building and simulating a membrane model of the polar and septal regions of E. coli membrane, which contains five lipid types: CL lipids with two types of acyl chains and phosphatidylethanolamine lipids with three types of acyl chains. It is our hope that CHARMM-GUI Membrane Builder becomes a useful tool for simulation studies to better understand the structure and dynamics of proteins and lipids in realistic biological membrane environments.

ST-analyzer: a web-based user interface for simulation trajectory analysis

Molecular dynamics (MD) simulation has become one of the key tools to obtain deeper insights into biological systems using various levels of descriptions such as all-atom, united-atom, and coarse-grained models. Recent advances in computing resources and MD programs have significantly accelerated the simulation time and thus increased the amount of trajectory data. Although many laboratories routinely perform MD simulations, analyzing MD trajectories is still time consuming and often a difficult task. ST-analyzer, http://im.bioinformatics.ku.edu/st-analyzer, is a standalone graphical user interface (GUI) toolset to perform various trajectory analyses. ST-analyzer has several outstanding features compared to other existing analysis tools: (i) handling various formats of trajectory files from MD programs, such as CHARMM, NAMD, GROMACS, and Amber, (ii) intuitive web-based GUI environment--minimizing administrative load and reducing burdens on the user from adapting new software environments, (iii) platform independent design--working with any existing operating system, (iv) easy integration into job queuing systems--providing options of batch processing either on the cluster or in an interactive mode, and (v) providing independence between foreground GUI and background modules--making it easier to add personal modules or to recycle/integrate pre-existing scripts utilizing other analysis tools. The current ST-analyzer contains nine main analysis modules that together contain 18 options, including density profile, lipid deuterium order parameters, surface area per lipid, and membrane hydrophobic thickness. This article introduces ST-analyzer with its design, implementation, and features, and also illustrates practical analysis of lipid bilayer simulations.

The high resolution structure of tyrocidine A reveals an amphipathic dimer

Tyrocidine A, one of the first antibiotics ever to be discovered, is a cyclic decapeptide that binds to membranes of target bacteria, disrupting their integrity. It is active against a broad spectrum of Gram-positive organisms, and has recently engendered interest as a potential scaffold for the development of new drugs to combat antibiotic-resistant pathogens. We present here the X-ray crystal structure of tyrocidine A at a resolution of 0.95Å. The structure reveals that tyrocidine forms an intimate and highly amphipathic homodimer made up of four beta strands that associate into a single, highly curved antiparallel beta sheet. We used surface plasmon resonance and potassium efflux assays to demonstrate that tyrocidine binds tightly to mimetics of bacterial membranes with an apparent dissociation constant (K(D)) of 10 μM, and efficiently permeabilizes bacterial cells at concentrations equal to and below the K(D). Using variant forms of tyrocidine in which the fluorescent probe p-cyano-phenylalanine had been inserted on either the polar or apolar face of the molecule, we performed fluorescence quenching experiments, using both water-soluble and membrane-embedded quenchers. The quenching results, together with the structure, strongly support a membrane association model in which the convex, apolar face of tyrocidine's beta sheet is oriented toward the membrane interior, while the concave, polar face is presented to the aqueous phase.

Structure, antimicrobial activities and mode of interaction with membranes of novel [corrected] phylloseptins from the painted-belly leaf frog, Phyllomedusa sauvagii

Transcriptomic and peptidomic analysis of skin secretions from the Painted-belly leaf frog Phyllomedusa sauvagii led to the identification of 5 novel phylloseptins (PLS-S2 to -S6) and also of phylloseptin-1 (PSN-1, here renamed PLS-S1), the only member of this family previously isolated in this frog. Synthesis and characterization of these phylloseptins revealed differences in their antimicrobial activities. PLS-S1, -S2, and -S4 (79–95% amino acid sequence identity; net charge  = +2) were highly potent and cidal against Gram-positive bacteria, including multidrug resistant S. aureus strains, and killed the promastigote stage of Leishmania infantum, L. braziliensis and L. major. By contrast, PLS-S3 (95% amino acid identity with PLS-S2; net charge  = +1) and -S5 (net charge  = +2) were found to be almost inactive against bacteria and protozoa. PLS-S6 was not studied as this peptide was closely related to PLS-S1. Differential scanning calorimetry on anionic and zwitterionic multilamellar vesicles combined with circular dichroism spectroscopy and membrane permeabilization assays on bacterial cells indicated that PLS-S1, -S2, and -S4 are structured in an amphipathic α-helix that disrupts the acyl chain packing of anionic lipid bilayers. As a result, regions of two coexisting phases could be formed, one phase rich in peptide and the other lipid-rich. After reaching a threshold peptide concentration, the disruption of lipid packing within the bilayer may lead to local cracks and disintegration of the microbial membrane. Differences in the net charge, α-helical folding propensity, and/or degree of amphipathicity between PLS-S1, -S2 and -S4, and between PLS-S3 and -S5 appear to be responsible for their marked differences in their antimicrobial activities. In addition to the detailed characterization of novel phylloseptins from P. sauvagii, our study provides additional data on the previously isolated PLS-S1 and on the mechanism of action of phylloseptins.

Improving the CHARMM force field for polyunsaturated fatty acid chains

CHARMM36 (C36) is the most up-to-date pairwise additive all-atom lipid force field and is able to accurately represent bilayer properties of saturated and monounsaturated lipid molecules in the natural constant particle, pressure, and temperature (NPT) ensemble. However, molecular dynamics (MD) simulations on 1-stearoyl-2-docosahexaenoyl-sn-glycerco-3-phosphocholine (SDPC) bilayers of the polyunsaturated fatty acid (PUFA) chains result in inaccuracies of the surface area per lipid (SA), deuterium order parameters (S(CD)), and X-ray form factors. Therefore, in this study, high-level quantum mechanical calculations are used to improve the dihedral potential of neighboring double bonds, and the corresponding force field is referred to as C36p. The SA for SDPC at 303 K increases from 63.2 ± 0.2 (C36) to 70.8 ± 0.2 (C36p) Å(2) and agrees favorably with X-ray diffraction results at 297 K. The resulting S(CD) are in excellent agreement with experimental values of both the sn-1 and sn-2 chains. Calculated NMR (13)C relaxation times and X-ray form factors from MD simulations of SDPC bilayers also agree with experiments. MD simulations of 1,2-diarachidonyl-phosphatidylcholine (DAPC) bilayers are used to further validate our force field parameters on a lipid with both chains containing PUFAs. As expected, the thickness of DAPC bilayers is reduced, and the SA is increased compared to the SDPC bilayers. This update in the PUFA force field should allow for accurate MD simulations of PUFA-containing bilayers in the NPT ensemble.