Free energy difference in indolicidin attraction to eukaryotic and prokaryotic model cell membranes

Accelerating Lipid Flip-Flop at Low Concentrations: A General Mechanism for Membrane Binding Peptides

We report a physicochemical investigation of the lipid transport properties of model lipid membranes in the presence of the antimicrobial peptide indolicidin through comparisons of experimental SANS/SAXS scattering techniques to fully atomistic molecular dynamics simulations. In agreement with the experiment, we show that upon peripheral binding of the peptides, even at low concentrations, lipid flip-flop dynamics is greatly accelerated. Computer modeling elucidates the interplay between structural changes and lipid dynamics induced by peptides and proposes a mechanism for the mode of action of antimicrobial peptides, assessing the major role of entropy for the catalysis of the flipping events. The mechanism introduced here is universal for all peptides with preferential peripheral binding to the membrane as it does not depend on the specific amino acid sequence.

Interaction of Uperin Peptides with Model Membranes: Molecular Dynamics Study

The interaction of antimicrobial and amyloid peptides with cell membranes is a critical step in their activities. Peptides of the uperin family obtained from the skin secretion of Australian amphibians demonstrate antimicrobial and amyloidogenic properties. All-atomic molecular dynamics and an umbrella sampling approach were used to study the interaction of uperins with model bacterial membrane. Two stable configurations of peptides were found. In the bound state, the peptides in helical form were located right under the head group region in parallel orientation with respect to the bilayer surface. Stable transmembrane configuration was observed for wild-type uperin and its alanine mutant in both alpha-helical and extended unstructured forms. The potential of mean force characterized the process of peptide binding from water to the lipid bilayer and its insertion into the membrane, and revealed that the transition of uperins from the bound state to the transmembrane position was accompanied by the rotation of peptides and passes through the energy barrier of 4-5 kcal/mol. Uperins have a weak effect on membrane properties.

Design of therapeutically improved analogue of the antimicrobial peptide, indolicidin, using a glycosylation strategy

Indolicidin is a member of cathelicidin family which displays broad spectrum antimicrobial activity. Severe toxicity and aggregation propensity associated with indolicidin pose a huge limitation to its probable therapeutic application. We are reporting the use of glycosylation strategy to design an analogue of indolicidin and subsequently explore structural and functional effects of sugar on it. Our study led to the design of a potent antibacterial glycosylated peptide, [βGlc-T9,K7]indolicidin, which showed decreased toxicity against erythrocytes and macrophage cells and thus a higher therapeutic selectivity. The incorporation of sugar also increased the solubility of the peptide. The mode of bacterial killing, functional stability, LPS binding, and cytokine inhibitory potential of the peptide, however, seemed unaffected upon glycosylation. Absence of significant changes in structure upon glycosylation accounts for the possibly retained functions and mode of action of the peptide. Our report thus presents the designing of an indolicidin analogue with improved therapeutic potential by substituting aromatic amino acid with glycosylated amino acid as a promising strategy for the first time.

Molecular Dynamics Simulation Study of Adsorption of Bioinspired Oligomers on Alumina Surfaces

The adsorption of small oligomers on a model metal oxide surface was studied with atomistically detailed molecular dynamics simulations. The oligomers consisted of two different repeat units: a maleimide, which contains a catechol functional group as in the dopamine residue found in marine adhesive proteins, and a methyl acrylate. A hydroxylated alumina surface was used as the model metal oxide surface. Adsorption interactions were investigated in aqueous as well as anhydrous conditions. In anhydrous conditions, the model oligomers displayed strong adsorption interactions with the surface. However, in aqueous conditions, the adsorption interactions were significantly weakened because of the competition with the water molecules for adsorption sites near the surface. Catechol functional groups in the model oligomers were found to play an important role in adsorption interactions with the alumina surface via hydrogen bonds. However, diverse adsorption properties were observed depending on compositions and sequences of two different repeat units and self-aggregations, indicating that the hydrogen bonding capability of catechol groups is not the sole factor determining adsorption properties.

Assessing the interaction between surfactant-like peptides and lipid membranes

Positively charged A_nK peptides have a higher preference for interacting with the DPPG membrane than remaining hydrated in bulk water.

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

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

Sampling errors in free energy simulations of small molecules in lipid bilayers

Free energy simulations are a powerful tool for evaluating the interactions of molecular solutes with lipid bilayers as mimetics of cellular membranes. However, these simulations are frequently hindered by systematic sampling errors. This review highlights recent progress in computing free energy profiles for inserting molecular solutes into lipid bilayers. Particular emphasis is placed on a systematic analysis of the free energy profiles, identifying the sources of sampling errors that reduce computational efficiency, and highlighting methodological advances that may alleviate sampling deficiencies. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.

Free Energy Calculations for the Peripheral Binding of Proteins/Peptides to an Anionic Membrane. 1. Implicit Membrane Models

The binding of peptides and proteins to the surface of complex lipid membranes is important in many biological processes such as cell signaling and membrane remodeling. Computational studies can aid experiments by identifying physical interactions and structural motifs that determine the binding affinity and specificity. However, previous studies focused on either qualitative behaviors of protein/membrane interactions or the binding affinity of small peptides. Motivated by this observation, we set out to develop computational protocols for bimolecular binding to charged membranes that are applicable to both peptides and large proteins. In this work, we explore a method based on an implicit membrane/solvent model (generalized Born with a simple switching in combination with the Gouy-Chapman-Stern model for a charged interface), which we expect to lead to useful results when the binding does not implicate significant membrane deformation and local demixing of lipids. We show that the binding free energy can be efficiently computed following a thermodynamic cycle similar to protein-ligand binding calculations, especially when a Bennett acceptance ratio based protocol is used to consider both the membrane bound and solution conformational ensembles. Test calculations on a series of peptides show that our computational approach leads to binding affinities in encouraging agreement with experimental data, including for the challenging example of the bringing of flexible MARCKS-ED peptides to membranes. The calculations highlight that for a membrane with a significant fraction of anionic lipids, it is essential to include the effect of ion adsorption using the Stern model, which significantly modifies the effective surface charge. This implicit membrane model based computational protocol helps lay the groundwork for more systematic analysis of protein/peptide binding to membranes of complex shape and composition.

Indolicidin binding induces thinning of a lipid bilayer

We use all-atom molecular dynamics simulations on a massive scale to compute the standard binding free energy of the 13-residue antimicrobial peptide indolicidin to a lipid bilayer. The analysis of statistical convergence reveals systematic sampling errors that correlate with reorganization of the bilayer on the microsecond timescale and persist throughout a total of 1.4 ms of sampling. Consistent with experimental observations, indolicidin induces membrane thinning, although the simulations significantly overestimate the lipophilicity of the peptide.

ChemNetworks: a complex network analysis tool for chemical systems

Many intermolecular chemical interactions persist across length and timescales and can be considered to form a "network" or "graph." Obvious examples include the hydrogen bond networks formed by polar solvents such as water or alcohols. In fact, there are many similarities between intermolecular chemical networks like those formed by hydrogen bonding and the complex and distributed networks found in computer science. Contemporary network analyses are able to dissect the complex local and global changes that occur within the network over multiple time and length scales. This work discusses the ChemNetworks software, whose purpose is to process Cartesian coordinates of chemical systems into a network/graph formalism and apply topological network analyses that include network neighborhood, the determination of geodesic paths, the degree census, direct structural searches, and the distribution of defect states of network. These properties can help to understand the network patterns and organization that may influence physical properties and chemical reactivity. The focus of ChemNetworks is to quantitatively describe intermolecular chemical networks of entire systems at both the local and global levels and as a function of time. The code is highly general, capable of converting a wide variety of systems into a chemical network formalism, including complex solutions, liquid interfaces, or even self-assemblies.

Biological activity of antibacterial peptides matches synergism between electrostatic and non electrostatic forces

Substitution of Ala 108 and Ala 111 in the 107-115 human lysozyme (hLz) fragment results in a 20-fold increased anti-staphylococcal activity while its hemolytic activity becomes significant (30%) at very high concentrations. This analog displays an additional positive charge near the N-terminus (108) and an extra Trp residue at the center of the molecule (111), indicating that this particular amino acid sequence improves its interaction with the bacterial plasma membrane. In order to understand the role of this arrangement in the membrane interaction, studies with model lipid membranes were carried out. The interactions of peptides, 107-115 hLz and the novel analog ([K(108)W(111)]107-115 hLz) with liposomes and lipid monolayers were evaluated by monitoring the changes in the fluorescence of the Trp residues and the variation of the monolayers surface pressure, respectively. Results obtained with both techniques revealed a significant affinity increase of [K(108)W(111)]107-115 hLz for lipids, especially when the membranes containing negatively charged lipids, such as phosphatidylglycerol. However, there is also a significant interaction with zwitterionic lipids, suggesting that other forces in addition to electrostatic interactions are involved in the binding. The analysis of adsorption isotherms and the insertion kinetics suggest that relaxation processes of the membrane structure are involved in the insertion process of novel peptide [K(108)W(111)]107-115 hLz but not in 107-115 hLz, probably by imposing a reorganization of water at the interphases. In this regard, the enhanced activity of peptide [K(108)W(111)]107-115 hLz may be explained by a synergistic effect between the increased electrostatic forces as well as the increased hydrophobic interactions.