Temperature-Dependent Re-alignment of the Short Multifunctional Peptide BP100 in Membranes Revealed by Solid-State NMR Spectroscopy and Molecular Dynamics Simulations.

BP100 is a cationic undecamer peptide with antimicrobial and cell-penetrating activities. The orientation of this amphiphilic α-helix in lipid bilayers was examined under numerous conditions using solid-state 19 F, 15 N and 2 H NMR. At high temperatures in saturated phosphatidylcholine lipids, BP100 lies flat on the membrane surface, as expected. Upon lowering the temperature towards the lipid phase transition, the helix is found to flip into an upright transmembrane orientation. In thin bilayers, this inserted state was stable at low peptide concentration, but thicker membranes required higher peptide concentrations. In the presence of lysolipids, the inserted state prevailed even at high temperature. Molecular dynamics simulations suggest that BP100 monomer insertion can be stabilized by snorkeling lysine side chains. These results demonstrate that even a very short helix like BP100 can span (and thereby penetrate through) a cellular membrane under suitable conditions.

A least-squares-fitting procedure for an efficient preclinical ranking of passive transport across the blood-brain barrier endothelium.

The treatment of various disorders of the central nervous system (CNS) is often impeded by the limited brain exposure of drugs, which is regulated by the human blood-brain barrier (BBB). The screening of lead compounds for CNS penetration is challenging due to the biochemical complexity of the BBB, while experimental determination of permeability is not feasible for all types of compounds. Here we present a novel method for rapid preclinical screening of libraries of compounds by utilizing advancements in computing hardware, with its foundation in transition-based counting of the flux. This method has been experimentally validated for in vitro permeabilities and provides atomic-level insights into transport mechanisms. Our approach only requires a single high-temperature simulation to rank a compound relative to a library, with a typical simulation time converging within 24 to 72 h. The method offers unbiased thermodynamic and kinetic information to interpret the passive transport of small-molecule drugs across the BBB.

Development of membrane-active peptide therapeutics in oncology.

Membrane-active peptides play an essential role in many living organisms and their immune systems and counter many infectious diseases. Many have dual or multiple mechanisms and can synergize with other molecules, like peptides, proteins, and small molecules. Although membrane-active peptides have been intensively studied in the past decades and more than 3500 sequences have been identified, only a few received approvals from the US Food and Drug Administration. In this review, we investigated all the peptide therapeutics that have entered the market or were subjected to preclinical and clinical studies to understand how they succeeded. With technological advancement (e.g., chemical modifications and pharmaceutical formulations) and a better understanding of the mechanism of action and the potential targets, we found at least five membrane-active peptide drugs that have entered preclinical/clinical phases and show promising results for cancer treatment. We summarized our findings in this review and provided insights into membrane-active anticancer peptide therapeutics.

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

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

Simulation-Guided Rational de Novo Design of a Small Pore-Forming Antimicrobial Peptide.

In the age of failing small-molecule antibiotics, tapping the near-infinite structural and chemical repertoire of antimicrobial peptides (AMPs) offers one of the most promising routes toward developing next-generation antibacterial compounds. One of the key impediments en route is the lack of methodologies for systematic rational design and optimization of new AMPs. Here we present a new simulation-guided rational design approach and apply it to develop a potent new AMP. We show that unbiased atomic detail molecular dynamics (MD) simulations are able to predict structures formed by evolving peptide designs enabling structure-based rational fine-tuning of functional properties. Starting from a 14-residue poly leucine template we demonstrate the design of a minimalistic potent new AMP. Consisting of only four types of amino acids (LDKA), this peptide forms large pores in microbial membranes at very low peptide-to-lipid ratios (1:1000) and exhibits low micromolar activity against common Gram-positive and Gram-negative pathogenic bacteria. Remarkably, the four amino acids were sufficient to encode preferential poration of bacterial membranes with negligible damage to red blood cells at bactericidal concentrations. As the sequence is too short to span cellular membranes, pores are formed by stacking of channels in each bilayer leaflet.

Molecular Dynamics Simulations Are Redefining Our View of Peptides Interacting with Biological Membranes.

Ever since the first molecular mechanics computer simulations of biological molecules became possible, there has been the dream to study all complex biological phenomena in silico, simply bypassing the enormous experimental challenges and their associated costs. For this, two inherent requirements need to be met: First, the time scales achievable in simulations must reach up to the millisecond range and even longer. Second, the computational model must accurately reproduce what is measured experimentally. Despite some recent successes, the general consensus in the field to date has been that neither of these conditions have yet been met and that the dream will be realized, if at all, only in the distant future. In this Account, we show that this view is wrong; instead, we are actually in the middle of the in silico molecular dynamics (MD) revolution, which is reshaping how we think about protein function. The example explored in this Account is a recent advance in the field of membrane-active peptides (MAPs). MD simulations have succeeded in accurately capturing the process of peptide binding, folding, and partitioning into lipid bilayers as well as revealing how channels form spontaneously from polypeptide fragments and conduct ionic and other cargo across membranes, all at atomic resolution. These game-changing advances have been made possible by a combination of steadily advancing computational power, more efficient algorithms and techniques, clever accelerated sampling schemes, and thorough experimental verifications. The great advantage of MD is the spatial and temporal resolution, directly providing a molecular movie of a protein undergoing folding and cycling through a functional process. This is especially important for proteins with transitory functional states, such as pore-forming MAPs. Recent successes are demonstrated here for the large class of antimicrobial peptides (AMPs). These short peptides are an essential part of the nonadaptive immune system for many organisms, ubiquitous in nature, and of particular interest to the pharmaceutical industry in the age of rising bacterial resistance to conventional antibiotic treatments. Unlike integral membrane proteins, AMPs are sufficiently small to allow converged sampling with the unbiased high-temperature sampling methodology outlined here and are relatively easy to handle experimentally. At the same time, AMPs exhibit a wealth of complex and poorly understood interactions with lipid bilayers, which allow not only tuning and validation of the simulation methodology but also advancement of our knowledge of protein-lipid interactions at a fundamental level. Space constraints limit our discussion to AMPs, but the MD methodologies outlined here can be applied to all phenomena involving peptides in membranes, including cell-penetrating peptides, signaling peptides, viral channel forming peptides, and fusion peptides, as well as ab initio membrane protein folding and assembly. For these systems, the promise of MD simulations to predict the structure of channels and to provide complete-atomic-detail trajectories of the mechanistic processes underlying their biological functions appears to rapidly become a reality. The current challenge is to design joint experimental and computational benchmarks to verify and tune MD force fields. With this, MD will finally fulfill its promise to become an inexpensive, powerful, and easy-to-use tool providing atomic-detail insights to researchers as part of their investigations into membrane biophysics and beyond.

Role of the Interaction Motif in Maintaining the Open Gate of an Open Sodium Channel.

Voltage-gated sodium channels undergo transitions between open, closed, and inactivated states, enabling regulation of the translocation of sodium ions across membranes. A recently published crystal structure of the full-length prokaryotic NavMs crystal structure in the activated open conformation has revealed the presence of a novel motif consisting of an extensive network of salt bridges involving residues in the voltage sensor, S4-S5 linker, pore, and C-terminal domains. This motif has been proposed to be responsible for maintaining an open conformation that enables ion translocation through the channel. In this study, we have used long-time molecular dynamics calculations without artificial restraints to demonstrate that the interaction network of full-length NavMs indeed prevents a rapid collapse and closure of the gate, in marked difference to earlier studies of the pore-only construct in which the gate had to be restrained to remain open. Interestingly, a frequently discussed "hydrophobic gating" mechanism at nanoscopic level is also observed in our simulations, in which the discontinuous water wire close to the gate region leads to an energetic barrier for ion conduction. In addition, we demonstrate the effects of in silico mutations of several of the key residues in the motif on the open channel's stability and functioning, correlating them with existing functional studies on this channel and homologous disease-associated mutations in human sodium channels; we also examine the effects of truncating/removing the voltage sensor and C-terminal domains in maintaining an open gate.

Correction to: Computed Free Energies of Peptide Insertion into Bilayers are Independent of Computational Method.

The original version of the article unfortunately contained an error in NIH support grant number RO1-GM74639 in the Acknowledgements section. The correct grant number is RO1-GM74637. This has been corrected with this erratum.

Computed Free Energies of Peptide Insertion into Bilayers are Independent of Computational Method.

We show that the free energy of inserting hydrophobic peptides into lipid bilayer membranes from surface-aligned to transmembrane inserted states can be reliably calculated using atomistic models. We use two entirely different computational methods: high temperature spontaneous peptide insertion calculations as well as umbrella sampling potential-of-mean-force (PMF) calculations, both yielding the same energetic profiles. The insertion free energies were calculated using two different protein and lipid force fields (OPLS protein/united-atom lipids and CHARMM36 protein/all-atom lipids) and found to be independent of the simulation parameters. In addition, the free energy of insertion is found to be independent of temperature for both force fields. However, we find major difference in the partitioning kinetics between OPLS and CHARMM36, likely due to the difference in roughness of the underlying free energy surfaces. Our results demonstrate not only a reliable method to calculate insertion free energies for peptides, but also represent a rare case where equilibrium simulations and PMF calculations can be directly compared.

Charged Antimicrobial Peptides Can Translocate across Membranes without Forming Channel-like Pores.

How can highly charged, cationic antimicrobial peptides (AMPs) translocate across hydrophobic lipid bilayers despite the prohibitive energetic penalty to do so? A common explanation has been the formation of peptide-lined channels. However, for most AMPs, no structures of membrane pores have been found despite clear evidence of membrane leakage and antimicrobial activity. The study here suggests an alternative and simple reason: for the AMP PGLa from Xenopus laevis (charge +5), such pores are not needed to explain both leakage and peptide translocation. Elevated-temperature multimicrosecond equilibrium simulations at all-atomistic level reveal that peptides spontaneously translocate across the membrane individually on a timescale of tens of microseconds, without forming pores. Both surface-bound peptides and lipids assist in the one-by-one translocation of the charged side chains. Single peptides can remain in a transmembrane orientation for many microseconds, snorkeling some charged residues to one interface and some to the opposite, but without inducing a water channel. Instead of stable pores, short-lived water bridges occur when two or three peptides connect at their termini, allowing both ion translocation and lipid flip-flop via a brushlike mechanism usually involving the C terminus of one peptide. The results here suggest that for some specific antimicrobial and other membrane active peptides, pore formation may not have to be invoked at all to explain peptide translocation and membrane permeabilization, which may explain why no channel structures for them have been determined experimentally.

Transmembrane helices containing a charged arginine are thermodynamically stable.

Hydrophobic amino acids are abundant in transmembrane (TM) helices of membrane proteins. Charged residues are sparse, apparently due to the unfavorable energetic cost of partitioning charges into nonpolar phases. Nevertheless, conserved arginine residues within TM helices regulate vital functions, such as ion channel voltage gating and integrin receptor inactivation. The energetic cost of arginine in various positions along hydrophobic helices has been controversial. Potential of mean force (PMF) calculations from atomistic molecular dynamics simulations predict very large energetic penalties, while in vitro experiments with Sec61 translocons indicate much smaller penalties, even for arginine in the center of hydrophobic TM helices. Resolution of this conflict has proved difficult, because the in vitro assay utilizes the complex Sec61 translocon, while the PMF calculations rely on the choice of simulation system and reaction coordinate. Here we present the results of computational and experimental studies that permit direct comparison with the Sec61 translocon results. We find that the Sec61 translocon mediates less efficient membrane insertion of Arg-containing TM helices compared with our computational and experimental bilayer-insertion results. In the simulations, a combination of arginine snorkeling, bilayer deformation, and peptide tilting is sufficient to lower the penalty of Arg insertion to an extent such that a hydrophobic TM helix with a central Arg residue readily inserts into a model membrane. Less favorable insertion by the translocon may be due to the decreased fluidity of the endoplasmic reticulum (ER) membrane compared with pure palmitoyloleoyl-phosphocholine (POPC). Nevertheless, our results provide an explanation for the differences between PMF- and experiment-based penalties for Arg burial.

Molecular dynamics of ion transport through the open conformation of a bacterial voltage-gated sodium channel.

The crystal structure of the open conformation of a bacterial voltage-gated sodium channel pore from Magnetococcus sp. (NaVMs) has provided the basis for a molecular dynamics study defining the channel's full ion translocation pathway and conductance process, selectivity, electrophysiological characteristics, and ion-binding sites. Microsecond molecular dynamics simulations permitted a complete time-course characterization of the protein in a membrane system, capturing the plethora of conductance events and revealing a complex mixture of single and multi-ion phenomena with decoupled rapid bidirectional water transport. The simulations suggest specific localization sites for the sodium ions, which correspond with experimentally determined electron density found in the selectivity filter of the crystal structure. These studies have also allowed us to identify the ion conductance mechanism and its relation to water movement for the NavMs channel pore and to make realistic predictions of its conductance properties. The calculated single-channel conductance and selectivity ratio correspond closely with the electrophysiology measurements of the NavMs channel expressed in HEK 293 cells. The ion translocation process seen in this voltage-gated sodium channel is clearly different from that exhibited by members of the closely related family of voltage-gated potassium channels and also differs considerably from existing proposals for the conductance process in sodium channels. These studies simulate sodium channel conductance based on an experimentally determined structure of a sodium channel pore that has a completely open transmembrane pathway and activation gate.

How reliable are molecular dynamics simulations of membrane active antimicrobial peptides?

Membrane-active antimicrobial peptides (AMPs) are challenging to study experimentally, but relatively easy to investigate using molecular dynamics (MD) computer simulations. For this reason, a large number of MD studies of AMPs have been reported over recent years. Yet relatively little effort has focused on the validity of such simulations. Are these results reliable, and do they agree with what is known experimentally? And how much meaningful information can be obtained? To answer these questions, we demonstrate here some of the requirements and limitations of running MD simulations for several common AMPs: PGLa, melittin, maculatin and BP100. The two most important findings are: (a) simulation results depend strongly on force field parameters, making experimental verification of the simulations obligatory, and (b) slow orientational and conformational fluctuations mean that much longer sampling timescales (multi-µs) are needed if quantitative agreement between simulation averages and experimental data is to be achieved. This article is part of a Special Issue entitled: Interfacially Active Peptides and Proteins. Guest Editors: William C. Wimley and Kalina Hristova.

Insights from Micro-second Atomistic Simulations of Melittin in Thin Lipid Bilayers.

The membrane disruption and pore-forming mechanism of melittin has been widely explored by experiments and computational studies. However, the precise mechanism is still enigmatic, and further study is required to turn antimicrobial peptides into future promising drugs against microbes. In this study, unbiased microsecond (µs) time scale (total 17 µs) atomistic molecular dynamics simulation were performed on multiple melittin systems in 1,2-dimyristoyl-sn-glycero-3-phosphocholine membrane to capture the various events during the membrane disorder produced by melittin. We observed bent U-shaped conformations of melittin, penetrated deeply into the membrane in all simulations, and a special double U-shaped structure. However, no peptide transmembrane insertion, nor pore formation was seen, indicating that these processes occur on much longer timescales, and suggesting that many prior computational studies of melittin were not sufficiently unbiased.

Melittin can permeabilize membranes via large transient pores.

Membrane active peptides are known to porate lipid bilayers, but their exact permeabilization mechanism and the structure of the nanoaggregates they form in membranes have often been difficult to determine experimentally. For many sequences at lower peptide concentrations, transient leakage is observed in experiments, suggesting the existence of transient pores. For two well-know peptides, alamethicin and melittin, we show here that molecular mechanics simulations i) can directly distinguish equilibrium poration and non-equilibrium transient leakage processes, and ii) can be used to observe the detailed pore structures and mechanism of permeabilization in both cases. Our results are in very high agreement with numerous experimental evidence for these two peptides. This suggests that molecular simulations can capture key membrane poration phenomena directly and in the future may develop to be a useful tool that can assist experimental peptide design.

Conformational states of melittin at a bilayer interface.

The distribution of peptide conformations in the membrane interface is central to partitioning energetics. Molecular-dynamics simulations enable characterization of in-membrane structural dynamics. Here, we describe melittin partitioning into dioleoylphosphatidylcholine lipids using CHARMM and OPLS force fields. Although the OPLS simulation failed to reproduce experimental results, the CHARMM simulation reported was consistent with experiments. The CHARMM simulation showed melittin to be represented by a narrow distribution of folding states in the membrane interface.

Masking the transmembrane region of the amyloid ß precursor protein as a safe means to lower amyloid ß production.

Introduction: Reducing brain levels of both soluble and insoluble forms of amyloid beta (Aß) remains the primary goal of most therapies that target Alzheimer's disease (AD). However, no treatment has so far resulted in patient benefit, and clinical trials of the most promising drug candidates have generally failed due to significant adverse effects. This highlights the need for safer and more selective ways to target and modulate Aß biogenesis. Methods: Peptide technology has advanced to allow reliable synthesis, purification, and delivery of once-challenging hydrophobic sequences. This is opening up new routes to target membrane processes associated with disease. Here we deploy a combination of atomic detail molecular dynamics (MD) simulations, living-cell Förster resonance energy transfer (FRET), and in vitro assays to elucidate the atomic-detail dynamics, molecular mechanisms, and cellular activity and selectivity of a membrane-active peptide that targets the Aß precursor protein (APP). Results: We demonstrate that Aß biogenesis can be downregulated selectively using an APP occlusion peptide (APPOP). APPOP inhibits Aß production in a dose-dependent manner, with a mean inhibitory concentration (IC50) of 450 nM toward exogenous APP and 50 nM toward endogenous APP in primary rat cortical neuronal cultures. APPOP does not impact the γ-secretase cleavage of Notch-1, or exhibit toxicity toward cultured primary rat neurons, suggesting that it selectively shields APP from proteolysis. Discussion: Drugs targeting AD need to be given early and for very long periods to prevent the onset of clinical symptoms. This necessitates being able to target Aß production precisely and without affecting the activity of key cellular enzymes such as γ-secretase for other substrates. Peptides offer a powerful way for targeting key pathways precisely, thereby reducing the risk of adverse effects. Here we show that protecting APP from proteolytic processing offers a promising route to safely and specifically lower Aß burden. In particular, we show that the amyloid pathway can be targeted directly and specificically. This reduces the risk of off-target effects and paves the way for a safe prophylactic treatment.

Reorientation and dimerization of the membrane-bound antimicrobial peptide PGLa from microsecond all-atom MD simulations.

The membrane-active antimicrobial peptide PGLa from Xenopus laevis is known from solid-state (2)H-, (15)N-, and (19)F-NMR spectroscopy to occupy two distinct α-helical surface adsorbed states in membranes: a surface-bound S-state with a tilt angle of ~95° at low peptide/lipid molar ratio (P/L = 1:200), and an obliquely tilted T-state with a tilt angle of 127° at higher peptide concentration (P/L = 1:50). Using a rapid molecular-dynamics insertion protocol in combination with microsecond-scale simulation, we have characterized the structure of both states in detail. As expected, the amphiphilic peptide resides horizontally on the membrane surface in a monomeric form at a low P/L, whereas the T-state is seen in the simulations to be a symmetric antiparallel dimer, with close contacts between small glycine and alanine residues at the interface. The computed tilt angles and azimuthal rotations, as well as the quadrupolar splittings predicted from the simulations agree with the experimental NMR data. The simulations reveal many structural details previously inaccessible, such as the immersion depth of the peptide in the membrane and the packing of the dimerization interface. The study highlights the ability and limitations of current state-of-the-art multimicrosecond all-atom simulations of membrane-active peptides to complement experimental data from solid-state NMR.

Predicting Membrane-Active Peptide Dynamics in Fluidic Lipid Membranes.

Understanding the interactions between peptides and lipid membranes could not only accelerate the development of antimicrobial peptides as treatments for infections but also be applied to finding targeted therapies for cancer and other diseases. However, designing biophysical experiments to study molecular interactions between flexible peptides and fluidic lipid membranes has been an ongoing challenge. Recently, with hardware advances, algorithm improvements, and more accurate parameterizations (i.e., force fields), all-atom molecular dynamics (MD) simulations have been used as a "computational microscope" to investigate the molecular interactions and mechanisms of membrane-active peptides in cell membranes (Chen et al., Curr Opin Struct Biol 61:160-166, 2020; Ulmschneider and Ulmschneider, Acc Chem Res 51(5):1106-1116, 2018; Dror et al., Annu Rev Biophys 41:429-452, 2012). In this chapter, we describe how to utilize MD simulations to predict and study peptide dynamics and how to validate the simulations by circular dichroism, intrinsic fluorescent probe, membrane leakage assay, electrical impedance, and isothermal titration calorimetry. Experimentally validated MD simulations open a new route towards peptide design starting from sequence and structure and leading to desirable functions.

Integrated Design of a Membrane-Lytic Peptide-Based Intravenous Nanotherapeutic Suppresses Triple-Negative Breast Cancer.

Membrane-lytic peptides offer broad synthetic flexibilities and design potential to the arsenal of anticancer therapeutics, which can be limited by cytotoxicity to noncancerous cells and induction of drug resistance via stress-induced mutagenesis. Despite continued research efforts on membrane-perforating peptides for antimicrobial applications, success in anticancer peptide therapeutics remains elusive given the muted distinction between cancerous and normal cell membranes and the challenge of peptide degradation and neutralization upon intravenous delivery. Using triple-negative breast cancer as a model, the authors report the development of a new class of anticancer peptides. Through function-conserving mutations, the authors achieved cancer cell selective membrane perforation, with leads exhibiting a 200-fold selectivity over non-cancerogenic cells and superior cytotoxicity over doxorubicin against breast cancer tumorspheres. Upon continuous exposure to the anticancer peptides at growth-arresting concentrations, cancer cells do not exhibit resistance phenotype, frequently observed under chemotherapeutic treatment. The authors further demonstrate efficient encapsulation of the anticancer peptides in 20 nm polymeric nanocarriers, which possess high tolerability and lead to effective tumor growth inhibition in a mouse model of MDA-MB-231 triple-negative breast cancer. This work demonstrates a multidisciplinary approach for enabling translationally relevant membrane-lytic peptides in oncology, opening up a vast chemical repertoire to the arms race against cancer.