Analysis and evaluation of channel models: simulations of alamethicin

Lipid vesicle-based molecular robots

A molecular robot, which is a system comprised of one or more molecular machines and computers, can execute sophisticated tasks in many fields that span from nanomedicine to green nanotechnology. The core parts of molecular robots are fairly consistent from system to system and always include (i) a body to encapsulate molecular machines, (ii) sensors to capture signals, (iii) computers to make decisions, and (iv) actuators to perform tasks. This review aims to provide an overview of approaches and considerations to develop molecular robots. We first introduce the basic technologies required for constructing the core parts of molecular robots, describe the recent progress towards achieving higher functionality, and subsequently discuss the current challenges and outlook. We also highlight the applications of molecular robots in sensing biomarkers, signal communications with living cells, and conversion of energy. Although molecular robots are still in their infancy, they will unquestionably initiate massive change in biomedical and environmental technology in the not too distant future.

Effect of lipid saturation on the topology and oligomeric state of helical membrane polypeptides

Natural liquid crystalline membranes are made up of many different lipids carrying a mixture of saturated and unsaturated fatty acyl chains. Whereas in the past considerable attention has been paid to cholesterol content, the phospholipid head groups and the membrane surface charge the detailed fatty acyl composition was often considered less important. However, recent investigations indicate that the detailed fatty acyl chain composition has pronounced effects on the oligomerization of the transmembrane helical anchoring domains of the MHC II receptor or the membrane alignment of the cationic antimicrobial peptide PGLa. In contrast the antimicrobial peptides magainin 2 and alamethicin are less susceptible to lipid saturation. Using histidine-rich LAH4 designer peptides the high energetic contributions of lipid saturation in stabilizing transmembrane helical alignments are quantitatively evaluated. These observations can have important implications for the biological regulation of membrane proteins and should be taken into considerations during biophysical or structural experiments.

A Peptide from Budding Yeast GAPDH Serves as a Promising Antifungal against Cryptococcus neoformans

Infection of Cryptococcus neoformans is one of the leading causes of morbidity and mortality, particularly among immunocompromised patients. However, currently available drugs for the treatment of C. neoformans infection are minimal. Here, we report SP1, a peptide derived from glyceraldehyde-3-phosphate dehydrogenase (GAPDH) of Saccharomyces cerevisiae, efficiently kills C. neoformans and Cryptococcus gattii. SP1 causes damages to the capsule. Unlike many antimicrobial peptides, SP1 does not form pores on the cell membrane of C. neoformans. It interacts with membrane ergosterol and enters vacuole possibly through membrane trafficking. C. neoformans treated with SP1 show the apoptotic phenotypes such as imbalance of calcium ion homeostasis, reactive oxygen increment, phosphatidylserine exposure, and nuclear fragmentation. Our data imply that SP1 has the potential to be developed into a treatment option for cryptococcosis. IMPORTANCE Cryptococcus neoformans and Cryptococcus gattii can cause cryptococcosis, which has a high mortality rate. To treat the disease, amphotericin B and fluconazole are often used in clinic. However, amphotericin B has rather high renal toxicity, and tolerance to these drugs are quicky developed. The peptide SP1 derived from baker's yeast GAPDH shows antifungal function to kill Cryptococcus neoformans and Cryptococcus gattii efficiently with a high specificity, even for the drug-resistant strains. Our data demonstrate that SP1 induces the apoptosis-like death of Cryptococcus neoformans at low concentrations. The finding of this peptide may shed light on a new direction to treat cryptococcosis.

Antimicrobial peptides: mechanism of action and lipid-mediated synergistic interactions within membranes

Biophysical and structural studies of peptide-lipid interactions, peptide topology and dynamics have changed our view of how antimicrobial peptides insert and interact with membranes. Clearly, both peptides and lipids are highly dynamic, and change and mutually adapt their conformation, membrane penetration and detailed morphology on a local and a global level. As a consequence, peptides and lipids can form a wide variety of supramolecular assemblies in which the more hydrophobic sequences preferentially, but not exclusively, adopt transmembrane alignments and have the potential to form oligomeric structures similar to those suggested by the transmembrane helical bundle model. In contrast, charged amphipathic sequences tend to stay intercalated at the membrane interface. Although the membranes are soft and can adapt, at increasing peptide density they cause pronounced disruptions of the phospholipid fatty acyl packing. At even higher local or global concentrations the peptides cause transient membrane openings, rupture and ultimately lysis. Interestingly, mixtures of peptides such as magainin 2 and PGLa, which are stored and secreted naturally as a cocktail, exhibit considerably enhanced antimicrobial activities when investigated together in antimicrobial assays and also in pore forming experiments applied to biophysical model systems. Our most recent investigations reveal that these peptides do not form stable complexes but act by specific lipid-mediated interactions and the nanoscale properties of phospholipid bilayers.

Structure Changes of a Membrane Polypeptide under an Applied Voltage Observed with Surface-Enhanced 2D IR Spectroscopy

The structures of many membrane-bound proteins and polypeptides depend on the membrane potential. However, spectroscopically studying their structures under an applied field is challenging, because a potential is difficult to generate across more than a few bilayers. We study the voltage-dependent structures of the membrane-bound polypeptide, alamethicin, using a spectroelectrochemical cell coated with a rough, gold film to create surface plasmons. The plasmons sufficiently enhance the 2D IR signal to measure a single bilayer. The film is also thick enough to conduct current and thereby apply a potential. The 2D IR spectra resolve features from both 3₁₀- and α-helical structures and cross-peaks connecting the two. We observe changes in the peak intensity, not their frequencies, upon applying a voltage. A similar change occurs with pH, which is known to alter the angle of alamethicin relative to the surface normal. The spectra are modeled using a vibrational exciton Hamiltonian, and the voltage-dependent spectra are consistent with a change in angle of the 3₁₀- and α-helices in the membrane from 55 to 44°and from 31 to 60°, respectively. The 3₁₀- and α-helices are coupled by approximately 10 cm^-1. These experiments provide new structural information about alamethicin under a potential difference and demonstrate a technique that might be applied to voltage-gated membrane proteins and compared to molecular dynamics structures.

Influence of the dimensionality of the periodic boundary conditions on the transport of a drug-peptide complex across model cell membranes

Many research efforts are devoted to improving the efficiency of chemotherapy. One of the aspects is to facilitate the transport of drugs across the cell membranes by attaching the therapeutics to a carrier molecule. The current study focuses on computational investigation of such a system with doxorubicin as the model drug, which is covalently bound to a cell-penetrating peptide. The correct description of its membrane translocation at the molecular level requires proper choice of the model membrane and of the simulation parameters. For the purpose, two phospholipid bilayers are built, one containing solely DPPC and another with mixed lipid content mimicking the composition of a human erythrocyte membrane. Atomistic molecular dynamics simulations are carried out in two types of periodic boundary conditions (2D and 3D PBC), in order to assess the effect of the periodicity dimensionality on the intermolecular interactions. The evolution of some basic characteristics of the bilayers and of the drug-peptide complex is tracked: mass density profiles, electrostatic potentials, lateral diffusion coefficients and areas per lipid, lipid-complex radial distribution functions, secondary structure of the peptide and orientation of the drug relative to the membrane. Thus, the influence of the periodic boundary conditions is quantified and it shows that the mixed system in 3D PBC is the most suitable for analysis of the translocation of the transporting moiety across cell membranes.

The structure of the antimicrobial human cathelicidin LL-37 shows oligomerization and channel formation in the presence of membrane mimics

The human cathelicidin LL-37 serves a critical role in the innate immune system defending bacterial infections. LL-37 can interact with molecules of the cell wall and perforate cytoplasmic membranes resulting in bacterial cell death. To test the interactions of LL-37 and bacterial cell wall components we crystallized LL-37 in the presence of detergents and obtained the structure of a narrow tetrameric channel with a strongly charged core. The formation of a tetramer was further studied by cross-linking in the presence of detergents and lipids. Using planar lipid membranes a small but defined conductivity of this channel could be demonstrated. Molecular dynamic simulations underline the stability of this channel in membranes and demonstrate pathways for the passage of water molecules. Time lapse studies of E. coli cells treated with LL-37 show membrane discontinuities in the outer membrane followed by cell wall damage and cell death. Collectively, our results open a venue to the understanding of a novel AMP killing mechanism and allows the rational design of LL-37 derivatives with enhanced bactericidal activity.

The Mechanisms of Action of Cationic Antimicrobial Peptides Refined by Novel Concepts from Biophysical Investigations

Even 30 years after the discovery of magainins, biophysical and structural investigations on how these peptides interact with membranes can still bear surprises and add new interesting detail to how these peptides exert their antimicrobial action. Early on, using oriented solid-state NMR spectroscopy, it was found that the amphipathic helices formed by magainins are active when being oriented parallel to the membrane surface. More recent investigations indicate that this in-planar alignment is also found when PGLa and magainin in combination exert synergistic pore-forming activities, where studies on the mechanism of synergistic interaction are ongoing. In a related manner, the investigation of dimeric antimicrobial peptide sequences has become an interesting topic of research which bears promise to refine our views how antimicrobial action occurs. The molecular shape concept has been introduced to explain the effects of lipids and peptides on membrane morphology, locally and globally, and in particular of cationic amphipathic helices that partition into the membrane interface. This concept has been extended in this review to include more recent ideas on soft membranes that can adapt to external stimuli including membrane-disruptive molecules. In this manner, the lipids can change their shape in the presence of low peptide concentrations, thereby maintaining the bilayer properties. At higher peptide concentrations, phase transitions occur which lead to the formation of pores and membrane lytic processes. In the context of the molecular shape concept, the properties of lipopeptides, including surfactins, are shortly presented, and comparisons with the hydrophobic alamethicin sequence are made.

Pore Forming Properties of Alamethicin in Negatively Charged Floating Bilayer Lipid Membranes Supported on Gold Electrodes

Electrochemical impedance spectroscopy (EIS), atomic force microscopy (AFM), and photon polarization modulation infrared reflection absorption spectroscopy (PM-IRRAS) were employed to investigate the formation of alamethicin pores in negatively charged bilayers composed of a mixture of 1,2-dimyristoyl- sn-glycero-3-phosphocholine (DMPC) and egg-PG floating at gold (111) electrode surfaces modified by self-assembled monolayers of 1-thio-β-d-glucose (β-Tg). The EIS data showed that the presence of alamethicin decreases the membrane resistivity by about 1 order of magnitude. PM-IRRAS measurements provided information about the tilt angles of peptide helical axis with respect to the bilayer normal. The small tilt angles obtained for the peptide helical axis prove that the alamethicin molecules were inserted into the DMPC/egg-PG membranes. The tilt angles decreased when negative potentials were applied, which correlates with the observed decrease in membrane resistivity, indicating that ion pore formation is assisted by the transmembrane potential. Molecular resolution AFM images provided visual evidence that alamethicin molecules aggregate forming hexagonal porous 2D lattices with periodicities of 2.0 ± 0.2 nm. The pore formation by alamethicin in the negatively charged membrane was compared with the interaction of this peptide with a bilayer formed by zwitterionic lipids. The comparison of these results showed that alamethicin preferentially forms ion translocating pores in negatively charged phospholipid membranes.

Biophysical Investigations Elucidating the Mechanisms of Action of Antimicrobial Peptides and Their Synergism

Biophysical and structural investigations are presented with a focus on the membrane lipid interactions of cationic linear antibiotic peptides such as magainin, PGLa, LL37, and melittin. Observations made with these peptides are distinct as seen from data obtained with the hydrophobic peptide alamethicin. The cationic amphipathic peptides predominantly adopt membrane alignments parallel to the bilayer surface; thus the distribution of polar and non-polar side chains of the amphipathic helices mirror the environmental changes at the membrane interface. Such a membrane partitioning of an amphipathic helix has been shown to cause considerable disruptions in the lipid packing arrangements, transient openings at low peptide concentration, and membrane disintegration at higher peptide-to-lipid ratios. The manifold supramolecular arrangements adopted by lipids and peptides are represented by the 'soft membranes adapt and respond, also transiently' (SMART) model. Whereas molecular dynamics simulations provide atomistic views on lipid membranes in the presence of antimicrobial peptides, the biophysical investigations reveal interesting details on a molecular and supramolecular level, and recent microscopic imaging experiments delineate interesting sequences of events when bacterial cells are exposed to such peptides. Finally, biophysical studies that aim to reveal the mechanisms of synergistic interactions of magainin 2 and PGLa are presented, including unpublished isothermal titration calorimetry (ITC), circular dichroism (CD) and dynamic light scattering (DLS) measurements that suggest that the peptides are involved in liposome agglutination by mediating intermembrane interactions. A number of structural events are presented in schematic models that relate to the antimicrobial and synergistic mechanism of amphipathic peptides when they are aligned parallel to the membrane surface.

Channel current analysis estimates the pore-formation and the penetration of transmembrane peptides

We measured the current signal of the transmembrane model peptides using the barrel-stave, toroidal pore, and penetration models in order to establish a precise assignment of the channel signals.

Alamethicin self-assembling in lipid membranes: concentration dependence from pulsed EPR of spin labels

The antimicrobial action of the peptide antibiotic alamethicin (Alm) is commonly related to peptide self-assembling resulting in the formation of voltage-dependent channels in bacterial membranes, which induces ion permeation.

Structural Behavior of the Peptaibol Harzianin HK VI in a DMPC Bilayer: Insights from MD Simulations

Microsecond molecular dynamics simulations of harzianin HK VI (HZ) interacting with a dimyristoylphosphatidylcholine bilayer were performed at the condition of low peptide-to-lipid ratio. Two orientations of HZ molecule in the bilayer were found and characterized. In the orientation perpendicular to the bilayer surface, HZ induces a local thinning of the bilayer. When inserted into the bilayer parallel to its surface, HZ is located nearly completely within the hydrophobic region of the bilayer. A combination of solid-state NMR and circular dichroism experiments found the latter orientation to be dominant. An extended sampling simulation provided qualitative results and showed the same orientation to be a global minimum of free energy. The secondary structure of HZ was characterized, and it was found to be located in the 3₁₀-helical family. The specific challenges of computer simulation of nonpolar peptides are discussed briefly.

Alamethicin Supramolecular Organization in Lipid Membranes from 19F Solid-State NMR

Alamethicins (ALMs) are antimicrobial peptides of fungal origin. Their sequences are rich in hydrophobic amino acids and strongly interact with lipid membranes, where they cause a well-defined increase in conductivity. Therefore, the peptides are thought to form transmembrane helical bundles in which the more hydrophilic residues line a water-filled pore. Whereas the peptide has been well characterized in terms of secondary structure, membrane topology, and interactions, much fewer data are available regarding the quaternary arrangement of the helices within lipid bilayers. A new, to our knowledge, fluorine-labeled ALM derivative was prepared and characterized when reconstituted into phospholipid bilayers. As a part of these studies, C¹⁹F₃-labeled compounds were characterized and calibrated for the first time, to our knowledge, for ¹⁹F solid-state NMR distance and oligomerization measurements by centerband-only detection of exchange (CODEX) experiments, which opens up a large range of potential labeling schemes. The ¹⁹F-¹⁹F CODEX solid-state NMR experiments performed with ALM in POPC lipid bilayers and at peptide/lipid ratios of 1:13 are in excellent agreement with molecular-dynamics calculations of dynamic pentameric assemblies. When the peptide/lipid ratio was lowered to 1:30, ALM was found in the dimeric form, indicating that the supramolecular organization is tuned by equilibria that can be shifted by changes in environmental conditions.

Long-range electrostatic corrections in multipolar/polarizable QM/MM simulations

Taking long-range electrostatic effects into account in classical and hybrid quantum mechanics-molecular mechanics (QM/MM) simulations is necessary for an accurate description of the system under study. We have recently developed a method, termed long-range electrostatic corrections (LREC), for monopolar QM/MM calculations. Here, we present an extension of LREC for multipolar/polarizable QM/MM simulations within the LICHEM software package. Reaction barriers and QM-MM interaction energies converge with a LREC cutoff between 20 and 25 Å, in agreement with our previous results. Additionally, the LREC approach for the QM-MM interactions can be smoothly combined with standard shifting or Ewald summation methods in the MM calculations. We recommend the use of QM(LREC)/MM(PME), where the QM region is treated with LREC and the MM region is treated with particle mesh Ewald (PME) summation. This combination is an excellent compromise between simplicity, speed, and accuracy for large QM/MM simulations.

Interaction of peptides with cell membranes: insights from molecular modeling

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

Molecular dynamics simulation of water permeation through the alpha-hemolysin channel

The alpha-hemolysin (AHL) nanochannel is a non-selective channel that allows for uncontrolled transport of small molecules across membranes leading to cell death. Although it is a bacterial toxin, it has promising applications, ranging from drug delivery systems to nano-sensing devices. This study focuses on the transport of water molecules through an AHL nanochannel using molecular dynamics (MD) simulations. Our results show that AHL can quickly transport water across membranes. The first-passage time approach was used to estimate the diffusion coefficient and the mean exit time. To study the energetics of transport, the potential of mean force (PMF) of a water molecule along the AHL nanochannel was calculated. The results show that the energy barriers of water permeation across a nanopore are always positive along the channel and the values are close to thermal energy (kBT). These findings suggest that the observed quick permeation of water is due to small energy barriers and a hydrophobic inner channel surface resulting in smaller friction. We speculate that these physical mechanisms are important in how AHL causes cell death.

Multiscale modeling of droplet interface bilayer membrane networks

Droplet interface bilayer (DIB) networks are considered for the development of stimuli-responsive membrane-based materials inspired by cellular mechanics. These DIB networks are often modeled as combinations of electrical circuit analogues, creating complex networks of capacitors and resistors that mimic the biomolecular structures. These empirical models are capable of replicating data from electrophysiology experiments, but these models do not accurately capture the underlying physical phenomena and consequently do not allow for simulations of material functionalities beyond the voltage-clamp or current-clamp conditions. The work presented here provides a more robust description of DIB network behavior through the development of a hierarchical multiscale model, recognizing that the macroscopic network properties are functions of their underlying molecular structure. The result of this research is a modeling methodology based on controlled exchanges across the interfaces of neighboring droplets. This methodology is validated against experimental data, and an extension case is provided to demonstrate possible future applications of droplet interface bilayer networks.

A new smoothing function to introduce long-range electrostatic effects in QM/MM calculations

A new method to account for long range electrostatic contributions is proposed and implemented for quantum mechanics/molecular mechanics long range electrostatic correction (QM/MM-LREC) calculations. This method involves the use of the minimum image convention under periodic boundary conditions and a new smoothing function for energies and forces at the cutoff boundary for the Coulomb interactions. Compared to conventional QM/MM calculations without long-range electrostatic corrections, the new method effectively includes effects on the MM environment in the primary image from its replicas in the neighborhood. QM/MM-LREC offers three useful features including the avoidance of calculations in reciprocal space (k-space), with the concomitant avoidance of having to reproduce (analytically or approximately) the QM charge density in k-space, and the straightforward availability of analytical Hessians. The new method is tested and compared with results from smooth particle mesh Ewald (PME) for three systems including a box of neat water, a double proton transfer reaction, and the geometry optimization of the critical point structures for the rate limiting step of the DNA dealkylase AlkB. As with other smoothing or shifting functions, relatively large cutoffs are necessary to achieve comparable accuracy with PME. For the double-proton transfer reaction, the use of a 22 Å cutoff shows a close reaction energy profile and geometries of stationary structures with QM/MM-LREC compared to conventional QM/MM with no truncation. Geometry optimization of stationary structures for the hydrogen abstraction step by AlkB shows some differences between QM/MM-LREC and the conventional QM/MM. These differences underscore the necessity of the inclusion of the long-range electrostatic contribution.

Toward the correction of effective electrostatic forces in explicit-solvent molecular dynamics simulations: restraints on solvent-generated electrostatic potential and solvent polarization

Despite considerable advances in computing power, atomistic simulations under nonperiodic boundary conditions, with Coulombic electrostatic interactions and in systems large enough to reduce finite-size associated errors in thermodynamic quantities to within the thermal energy, are still not affordable. As a result, periodic boundary conditions, systems of microscopic size and effective electrostatic interaction functions are frequently resorted to. Ensuing artifacts in thermodynamic quantities are nowadays routinely corrected a posteriori, but the underlying configurational sampling still descends from spurious forces. The present study addresses this problem through the introduction of on-the-fly corrections to the physical forces during an atomistic molecular dynamics simulation. Two different approaches are suggested, where the force corrections are derived from special potential energy terms. In the first approach, the solvent-generated electrostatic potential sampled at a given atom site is restrained to a target value involving corrections for electrostatic artifacts. In the second approach, the long-range regime of the solvent polarization around a given atom site is restrained to the Born polarization, i.e., the solvent polarization corresponding to the ideal situation of a macroscopic system under nonperiodic boundary conditions and governed by Coulombic electrostatic interactions. The restraints are applied to the explicit-water simulation of a hydrated sodium ion, and the effect of the restraints on the structural and energetic properties of the solvent is illustrated. Furthermore, by means of the calculation of the charging free energy of a hydrated sodium ion, it is shown how the electrostatic potential restraint translates into the on-the-fly consideration of the corresponding free-energy correction terms. It is discussed how the restraints can be generalized to situations involving several solute particles. Although the present study considers a very simple system only, it is an important step toward the on-the-fly elimination of finite-size and approximate-electrostatic artifacts during atomistic molecular dynamics simulations.

In silico design of antimicrobial peptides

The rapid spread of drug-resistant pathogenic microbial strains has created an urgent need for the development of new anti-infective molecules, having different mechanism of action in comparison to existing drugs. Natural antimicrobial peptides (AMPs) represent a novel class of molecules with a broad spectrum of activity and a low rate in inducing bacterial resistance. In particular, linear alpha-helical cationic antimicrobial peptides are among the most widespread membrane-disruptive AMPs in nature, representing a particularly successful structural arrangement of the innate defense against microbes. However, until now, many AMPs have failed in clinical trials because of several drawbacks that strongly limit their applicability such as degradation, cytotoxicity, and high production cost. Thus, to overcome the limitations of native peptides, a rational in silico approach to AMPs design becomes a promising strategy that drastically reduce production costs and the time required for evaluation of activity and toxicity. This chapter focuses on the strategies and methods for de novo design of potentially active AMPs. In particular, statistical-based design strategies and MD methods for modelling AMPs are elucidated.

Long-timescale motions in glycerol-monopalmitate lipid bilayers investigated using molecular dynamics simulation

The occurrence of long-timescale motions in glycerol-1-monopalmitate (GMP) lipid bilayers is investigated based on previously reported 600 ns molecular dynamics simulations of a 2×8×8 GMP bilayer patch in the temperature range 302-338 K, performed at three different hydration levels, or in the presence of the cosolutes methanol or trehalose at three different concentrations. The types of long-timescale motions considered are: (i) the possible phase transitions; (ii) the precession of the relative collective tilt-angle of the two leaflets in the gel phase; (iii) the trans-gauche isomerization of the dihedral angles within the lipid aliphatic tails; and (iv) the flipping of single lipids across the two leaflets. The results provide a picture of GMP bilayers involving a rich spectrum of events occurring on a wide range of timescales, from the 100-ps range isomerization of single dihedral angles, via the 100-ns range of tilt precession motions, to the multi-μs range of phase transitions and lipid-flipping events.

Protein arcs may form stable pores in lipid membranes

Electron microscopy and atomic force microscopy images of cholesterol-dependent cytolysins and related proteins that form large pores in lipid membranes have revealed the presence of incomplete rings, or arcs. Some evidence indicates that these arcs are inserted into the membrane and induce membrane leakage, but other experiments seem to refute that. Could such pores, only partially lined by protein, be kinetically and thermodynamically stable? How would the lipids be structured in such a pore? Using the antimicrobial peptide protegrin-1 as a model, we test the stability of pores only partially lined by peptide using all-atom molecular dynamics simulations in POPC and POPE/POPG membranes. The data show that, whereas pure lipid pores close rapidly, pores partially lined by protegrin arcs are stable for at least 300 ns. Estimates of the thermodynamic stability of these arcs using line tension data and implicit solvent calculations show that these arcs can be marginally stable in both zwitterionic and anionic membranes. Arcs provide an explanation for the observed ion selectivity in protegrin electrophysiology experiments and could possibly be involved in other membrane permeabilization processes where lipids are thought to participate, such as those induced by antimicrobial peptides and colicins, as well as the Bax apoptotic pore.

Conformational and thermodynamic properties of non-canonical α,α-dialkyl glycines in the peptaibol Alamethicin: molecular dynamics studies

In this work, we investigate the structure, dynamic and thermodynamic properties of noncanonical disubstituted amino acids (α,α-dialkyl glycines), also known as non-natural amino acids, in the peptaibol Alamethicin. The amino acids under study are Aib (α-amino isobutyric acid or α-methyl alanine), Deg (α,α-diethyl glycine), Dpg (α,α-dipropyl glycine), Dibg (α,α-di-isobutyl glycine), Dhg (α,α-dihexyl glycine), DΦg (α,α-diphenyl glycine), Dbzg (α,α-dibenzyl glycine), Ac6c (α,α-cyclohexyl glycine), and Dmg (α,α-dihydroxymethyl glycine). It is hypothesized that these amino acids are able to induce well-defined secondary structure in peptidomimetics. To test this hypothesis, new peptidomimetics of Alamethicin were constructed by replacing the native Aib positions of Alamethicin by one or more new α,α-dialkyl glycines. Dhg and Ac6c demonstrated the capacity to induce well-defined α-helical structures. Dhg and Ac6c also promote the thermodynamic stabilization of these peptides in a POPC model membrane and are better alternatives to the Aib in Alamethicin. These noncanonical amino acids also improved secondary structure properties, revealing preorganization in water and maintenance of α helical structure in POPC. We show that it is possible to optimize the helicity and thermodynamic properties of native Alamethicin, and we suggest that these amino acids could be incorporated in other peptides with similar structural effect.

Treatment of microbial biofilms in the post-antibiotic era: prophylactic and therapeutic use of antimicrobial peptides and their design by bioinformatics tools

The treatment for biofilm infections is particularly challenging because bacteria in these conditions become refractory to antibiotic drugs. The reduced effectiveness of current therapies spurs research for the identification of novel molecules endowed with antimicrobial activities and new mechanisms of antibiofilm action. Antimicrobial peptides (AMPs) have been receiving increasing attention as potential therapeutic agents, because they represent a novel class of antibiotics with a wide spectrum of activity and a low rate in inducing bacterial resistance. Over the past decades, a large number of naturally occurring AMPs have been identified or predicted from various organisms as effector molecules of the innate immune system playing a crucial role in the first line of defense. Recent studies have shown the ability of some AMPs to act against microbial biofilms, in particular during early phases of biofilm development. Here, we provide a review of the antimicrobial peptides tested on biofilms, highlighting their advantages and disadvantages for prophylactic and therapeutic applications. In addition, we describe the strategies and methods for de novo design of potentially active AMPs and discuss how informatics and computational tools may be exploited to improve antibiofilm effectiveness.

Effects of PEGylation on the binding interaction of magainin 2 and tachyplesin I with lipid bilayer surface

Poly(ethylene glycol) (PEG)-grafted magainin 2 and tachyplesin I were simulated with lipid bilayers. In the simulations of PEGylated magainin 2 and tachyplesin I in water, both peptides are wrapped by PEG chains. The α-helical structure of PEGylated magainin 2 is broken in water, while the β-sheet of PEGylated tachyplesin I keeps stable, similar to the structural behavior of unPEGylated peptides, in agreement with experiments. Simulations of PEGylated peptides with lipid bilayers show that PEG chains block the electrostatic interaction between cationic residues of peptides and anionic phosphates of lipids, leading to the less binding of the peptide to the bilayer surface, which is observed more significantly for magainin 2 than for tachyplesin I. Since the random-coiled magainin 2 can be more completely covered by PEGs than does the β-sheet tachyplesin I, the PEGylation effect on the decreased binding is larger for magainin 2, showing the dependence of PEGylation on the peptide structure. These simulation findings qualitatively support the experimental observation of the different extents of the reduced membrane-permeabilizing activity for PEGylated magainin 2 and tachyplesin I.

A thermodynamic approach to alamethicin pore formation

The structure and energetics of alamethicin Rf30 monomer to nonamer in cylindrical pores of 5 to 11Å radius are investigated using molecular dynamics simulations in an implicit membrane model that includes the free energy cost of acyl chain hydrophobic area exposure. Stable, low energy pores are obtained for certain combinations of radius and oligomeric number. The trimer and the tetramer formed 6Å pores that appear closed while the larger oligomers formed open pores at their optimal radius. The hexamer in an 8Å pore and the octamer in an 11Å pore give the lowest effective energy per monomer. However, all oligomers beyond the pentamer have comparable energies, consistent with the observation of multiple conductance levels. The results are consistent with the widely accepted "barrel-stave" model. The N terminal portion of the molecule exhibits smaller tilt with respect to the membrane normal than the C terminal portion, resulting in a pore shape that is a hybrid between a funnel and an hourglass. Transmembrane voltage has little effect on the structure of the oligomers but enhances or decreases their stability depending on its orientation. Antiparallel bundles are lower in energy than the commonly accepted parallel ones and could be present under certain experimental conditions. Dry aggregates (without an aqueous pore) have lower average effective energy than the corresponding aggregates in a pore, suggesting that alamethicin pores may be excited states that are stabilized in part by voltage and in part by the ion flow itself.

Modeling peptide binding to anionic membrane pores

Peptide-induced pore formation in membranes can be dissected into two steps: pore formation and peptide binding to the pore. A computational method is proposed to study the second step in anionic membranes. The electrostatic potential is obtained from numerical solutions to the Poisson-Boltzmann equation and is then used in conjunction with IMM1 (implicit membrane model 1). A double charge layer model is used to incorporate the effects of the membrane dipole potential. Inhomogeneity of the charge density in the pore, characterized by explicit membrane simulations of toroidal pores, is included in the model. This approach was applied to two extensively studied peptides, magainin and melittin. In agreement with previous work, binding to toroidal pores is more favorable than binding to the flat membrane. The dependence of binding energy on anionic content exhibits different patterns for the two peptides, in correlation with the different lipid selectivity that has been observed experimentally.

Behavior of solvent-exposed hydrophobic groove in the anti-apoptotic Bcl-XL protein: clues for its ability to bind diverse BH3 ligands from MD simulations

Bcl-XL is a member of Bcl-2 family of proteins involved in the regulation of intrinsic pathway of apoptosis. Its overexpression in many human cancers makes it an important target for anti-cancer drugs. Bcl-XL interacts with the BH3 domain of several pro-apoptotic Bcl-2 partners. This helical bundle protein has a pronounced hydrophobic groove which acts as a binding region for the BH3 domains. Eight independent molecular dynamics simulations of the apo/holo forms of Bcl-XL were carried out to investigate the behavior of solvent-exposed hydrophobic groove. The simulations used either a twin-range cut-off or particle mesh Ewald (PME) scheme to treat long-range interactions. Destabilization of the BH3 domain-containing helix H2 was observed in all four twin-range cut-off simulations. Most of the other major helices remained stable. The unwinding of H2 can be related to the ability of Bcl-XL to bind diverse BH3 ligands. The loss of helical character can also be linked to the formation of homo- or hetero-dimers in Bcl-2 proteins. Several experimental studies have suggested that exposure of BH3 domain is a crucial event before they form dimers. Thus unwinding of H2 seems to be functionally very important. The four PME simulations, however, revealed a stable helix H2. It is possible that the H2 unfolding might occur in PME simulations at longer time scales. Hydrophobic residues in the hydrophobic groove are involved in stable interactions among themselves. The solvent accessible surface areas of bulky hydrophobic residues in the groove are significantly buried by the loop LB connecting the helix H2 and subsequent helix. These observations help to understand how the hydrophobic patch in Bcl-XL remains stable in the solvent-exposed state. We suggest that both the destabilization of helix H2 and the conformational heterogeneity of loop LB are important factors for binding of diverse ligands in the hydrophobic groove of Bcl-XL.

Direct visualization of the alamethicin pore formed in a planar phospholipid matrix

We present direct visualization of pores formed by alamethicin (Alm) in a matrix of phospholipids using electrochemical scanning tunneling microscopy (EC-STM). High-resolution EC-STM images show individual peptide molecules forming channels. The channels are not dispersed randomly in the monolayer but agglomerate forming 2D nanocrystals with a hexagonal lattice in which the average channel-channel distance is 1.90 ± 0.1 nm. The STM images suggest that each Alm is shared between the two adjacent channels. Every channel consists of six Alm molecules. Three or four of these molecules have the hydrophilic group oriented toward the center of the channel allowing for water column formation inside the channel. The dimensions of the central pore in the images are consistent with the dimension of the water column in a model of hexameric pore proposed in the literature. The images obtained in this work validate the barrel-stave model of the pore formed in phospholipid membranes by amphiphatic peptides. They also provide direct evidence for cluster formation by such pores.

The p7 protein of hepatitis C virus forms structurally plastic, minimalist ion channels

Hepatitis C virus (HCV) p7 is a membrane-associated oligomeric protein harboring ion channel activity. It is essential for effective assembly and release of infectious HCV particles and an attractive target for antiviral intervention. Yet, the self-assembly and molecular mechanism of p7 ion channelling are currently only partially understood. Using molecular dynamics simulations (aggregate time 1.2 µs), we show that p7 can form stable oligomers of four to seven subunits, with a bias towards six or seven subunits, and suggest that p7 self-assembles in a sequential manner, with tetrameric and pentameric complexes forming as intermediate states leading to the final hexameric or heptameric assembly. We describe a model of a hexameric p7 complex, which forms a transiently-open channel capable of conducting ions in simulation. We investigate the ability of the hexameric model to flexibly rearrange to adapt to the local lipid environment, and demonstrate how this model can be reconciled with low-resolution electron microscopy data. In the light of these results, a view of p7 oligomerization is proposed, wherein hexameric and heptameric complexes may coexist, forming minimalist, yet robust functional ion channels. In the absence of a high-resolution p7 structure, the models presented in this paper can prove valuable as a substitute structure in future studies of p7 function, or in the search for p7-inhibiting drugs.

Cyclodextrin-scaffolded alamethicin with remarkably efficient membrane permeabilizing properties and membrane current conductance

Bacterial resistance to classical antibiotics is a serious medical problem, which continues to grow. Small antimicrobial peptides represent a potential solution and are increasingly being developed as novel therapeutic agents. Many of these peptides owe their antibacterial activity to the formation of trans-membrane ion-channels resulting in cell lysis. However, to further develop the field of peptide antibiotics, a thorough understanding of their mechanism of action is needed. Alamethicin belongs to a class of peptides called peptaibols and represents one of these antimicrobial peptides. To examine the dynamics of assembly and to facilitate a thorough structural evaluation of the alamethicin ion-channels, we have applied click chemistry for the synthesis of templated alamethicin multimers covalently attached to cyclodextrin-scaffolds. Using oriented circular dichroism, calcein release assays, and single-channel current measurements, the α-helices of the templated multimers were demonstrated to insert into lipid bilayers forming highly efficient and remarkably stable ion-channels.

Charge distribution and imperfect amphipathicity affect pore formation by antimicrobial peptides

Antimicrobial peptides often permeabilize biological membranes via a pore mechanism. Two pore types have been proposed: toroidal, where the pore is partly lined by lipid, and barrel-stave, where a cylindrical pore is completely lined by peptides. What drives the preference of antimicrobial peptides for a certain pore type is not yet fully understood. According to neutron scattering and oriented circular dichroism, melittin and MG-H2 induce toroidal pores whereas alamethicin forms barrel-stave pores. In previous work we found that indeed melittin seems to favor toroidal pores whereas alamethicin favors cylindrical pores. Here we designed mutants of these two peptides and the magainin analog MG-H2, aimed to probe how the distribution of charges along the helix and its imperfectly amphipathic structure influence pore formation. Molecular dynamics (MD) simulations of the peptides in a pre-formed cylindrical pore have been performed. The duration of the simulations was 136ns to 216ns. We found that a melittin mutant with lysine 7 neutralized favors cylindrical pores whereas a MG-H2 mutant with lysines in the N-terminal half of these peptides neutralized and an alamethicin mutant with a positive charge at the position 7 form semitoroidal pores. These results suggest that charged residues within the N-terminal half are important for toroidal pore formation. Toroidal pores produced by MG-H2 are more disordered than the melittin pores, likely because of the charged residues located in the middle of the MG-H2 helix (K11 and K14). Imperfect amphipathicity of melittin seems to play a role in its preference for toroidal pores since the substitutions of charged residues located within the nonpolar face by hydrophobic residues suppress evolution of a toroidal pore. The mutations change the position of lysine 7 near the N-terminus, relative to the lower leaflet headgroups. The MD simulations also show that the melittin P14A mutant forms a toroidal pore, but its configuration diverges from that of melittin and it is probably metastable.

Molecular dynamics simulation of the antiamoebin ion channel: linking structure and conductance

Molecular-dynamics simulations were carried out to ascertain which of the potential multimeric forms of the transmembrane peptaibol channel, antiamoebin, is consistent with its measured conductance. Estimates of the conductance obtained through counting ions that cross the channel and by solving the Nernst-Planck equation yield consistent results, indicating that the motion of ions inside the channel can be satisfactorily described as diffusive. The calculated conductance of octameric channels is markedly higher than the conductance measured in single channel recordings, whereas the tetramer appears to be nonconducting. The conductance of the hexamer was estimated to be 115 ± 34 pS and 74 ± 20 pS, at 150 mV and 75 mV, respectively, in satisfactory agreement with the value of 90 pS measured at 75 mV. On this basis, we propose that the antiamoebin channel consists of six monomers. Its pore is large enough to accommodate K⁺ and Cl⁻ with their first solvation shells intact. The free energy barrier encountered by K⁺ is only 2.2 kcal/mol whereas Cl⁻ encounters a substantially higher barrier of nearly 5 kcal/mol. This difference makes the channel selective for cations. Ion crossing events are shown to be uncorrelated and follow Poisson statistics.