In silico partitioning and transmembrane insertion of hydrophobic peptides under equilibrium conditions

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.

The Pathway of a Transmembrane Helix Insertion into the Membrane Assisted by Sec61α Channel

The significant inconsistency between the experimental and simulation results of the free energy for the translocon-assisted insertion of the transmembrane helix (TMH) has not been reasonably explained. Understanding the mechanism of TMH insertion through the translocon is the key to solving this problem. In this study, we performed a series of coarse-grained molecular dynamics simulations and calculated the potential mean forces (PMFs) for three insertion processes of a hydrophobic TMH. The simulations reveal the pathway of the TMH insertion assisted by a translocon. The results indicate that the TMH contacts the top of the lateral gate first and then inserts down the lateral gate, which agrees with the sliding model. The TMH begins to transfer laterally to the bilayer when it is blocked by the plug and reaches the exit of the lateral gate, where there is a free energy minimum point. We also found that the connecting section between TM2 and TM3 of Sec61α prevented TMH from leaving the lateral gate and directly transitioning to the surface-bound state. These findings provide insight into the mechanism of the insertion of TMH through the translocon.

Spontaneous transfer of small peripheral peptides between supported lipid bilayer and giant unilamellar vesicles

Vesicular trafficking facilitates material transport between membrane-bound organelles. Membrane protein cargos are trafficked for relocation, recycling, and degradation during various physiological processes. In vitro fusion studies utilized synthetic lipid membranes to study the molecular mechanisms of vesicular trafficking and to develop synthetic materials mimicking the biological membrane trafficking. Various fusogenic conditions which can induce vesicular fusion have been used to establish synthetic systems that can mimic biological systems. Despite these efforts, the mechanisms underlying vesicular trafficking of membrane proteins remain limited and robust in vitro methods that can construct synthetic trafficking systems for membrane proteins between large membranes (>1 μm2) are unavailable. Here, we provide data to show the spontaneous transfer of small membrane-bound peptides (∼4 kD) between a supported lipid bilayer (SLB) and giant unilamellar vesicles (GUVs). We found that the contact between the SLB and GUVs led to the occasional but notable transfer of membrane-bound peptides in a physiological saline buffer condition (pH 7.4, 150 mM NaCl). Quantitative and dynamic time-lapse analyses suggested that the observed exchange occurred through the formation of hemi-fusion stalks between the SLB and GUVs. Larger protein cargos with a size of ∼77 kD could not be transferred between the SLB and GUVs, suggesting that the larger-sized cargos limited diffusion across the hemi-fusion stalk, which was predicted to have a highly curved structure. Compositional study showed Ni-chelated lipid head group was the essential component catalyzing the process. Our system serves as an example synthetic platform that enables the investigation of small-peptide trafficking between synthetic membranes and reveals hemi-fused lipid bridge formation as a mechanism of peptide transfer.

Antibacterial and Antiviral Properties of Chenopodin-Derived Synthetic Peptides

Antimicrobial peptides have been developed based on plant-derived molecular scaffolds for the treatment of infectious diseases. Chenopodin is an abundant seed storage protein in quinoa, an Andean plant with high nutritional and therapeutic properties. Here, we used computer- and physicochemical-based strategies and designed four peptides derived from the primary structure of Chenopodin. Two peptides reproduce natural fragments of 14 amino acids from Chenopodin, named Chen1 and Chen2, and two engineered peptides of the same length were designed based on the Chen1 sequence. The two amino acids of Chen1 containing amide side chains were replaced by arginine (ChenR) or tryptophan (ChenW) to generate engineered cationic and hydrophobic peptides. The evaluation of these 14-mer peptides on Staphylococcus aureus and Escherichia coli showed that Chen1 does not have antibacterial activity up to 512 µM against these strains, while other peptides exhibited antibacterial effects at lower concentrations. The chemical substitutions of glutamine and asparagine by amino acids with cationic or aromatic side chains significantly favoured their antibacterial effects. These peptides did not show significant hemolytic activity. The fluorescence microscopy analysis highlighted the membranolytic nature of Chenopodin-derived peptides. Using molecular dynamic simulations, we found that a pore is formed when multiple peptides are assembled in the membrane. Whereas, some of them form secondary structures when interacting with the membrane, allowing water translocations during the simulations. Finally, Chen2 and ChenR significantly reduced SARS-CoV-2 infection. These findings demonstrate that Chenopodin is a highly useful template for the design, engineering, and manufacturing of non-toxic, antibacterial, and antiviral peptides.

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.

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.

In Search of a Molecular View of Peptide-Lipid Interactions in Membranes

Lipid bilayer membranes are often represented as a continuous nonpolar slab with a certain thickness bounded by two more polar interfaces. Phenomena such as peptide binding to the membrane surface, folding, insertion, translocation, and diffusion are typically interpreted on the basis of this view. In this Perspective, I argue that this membrane representation as a hydrophobic continuum solvent is not adequate to understand peptide-lipid interactions. Lipids are not small compared to membrane-active peptides: their sizes are similar. Therefore, peptide diffusion needs to be understood in terms of free volume, not classical continuum mechanics; peptide solubility or partitioning in membranes cannot be interpreted in terms of hydrophobic mismatch between membrane thickness and peptide length; peptide folding and translocation, often involving cationic peptides, can only be understood if realizing that lipids adapt to the presence of peptides and the membrane may undergo considerable lipid redistribution in the process. In all of those instances, the detailed molecular interactions between the peptide residues and the lipid components are essential to understand the mechanisms involved.

Fifty Years of Biophysics at the Membrane Frontier

The author first describes his childhood in the South and the ways in which it fostered the values he has espoused throughout his life, his development of a keen fascination with science, and the influences that supported his progress toward higher education. His experiences in ROTC as a student, followed by two years in the US Army during the Vietnam War, honed his leadership skills. The bulk of the autobiography is a chronological journey through his scientific career, beginning with arrival at the University of California, Irvine in 1972, with an emphasis on the postdoctoral students and colleagues who have contributed substantially to each phase of his lab's progress. White's fundamental findings played a key role in the development of membrane biophysics, helping establish it as fertile ground for research. A story gradually unfolds that reveals the deeply collaborative and painstakingly executed work necessary for a successful career in science.

Three-stage model of helical membrane protein folding: Role of membrane-water interface as the intermediate stage vestibule for TM helices during their in membrano assembly

Integral membrane proteins (MPs) are dominated by transmembrane α-helical (TMH) proteins playing critical roles in cellular signaling processes. These proteins display a wide range of sizes from one TMH domain to at least 26 TMH domains and diverse structural folds. A common feature of most of these folds is the TM orientation of the helical domains and the approximately parallel packing of these domains into helical bundles of varying stability, however, it has been challenging to study the folding of these proteins experimentally. The contribution of helix stabilization in membrane and interface to the folding energy landscape are investigated here for the full range of TMH protein sizes containing 1 TM domain (1-TMH protein) to 24 TM domains (24-TMH protein) for all TMH proteins with available structures using structural bioinformatics based hydropathy analysis. The TM helix insertion stabilization energies from Water to membrane-water Interface (WAT→INT energies) are on average half of those insertion energies from water to transmembrane orientation (WAT→TM energies) for the whole polytopic helical membrane proteome (1-TMH to 24-TMH proteins). This suggests a potentially dominant role of the membrane-water interface as a viable holding vestibule for the TM helices during their release from the translocon. This provides proteome-level evidence for the broadly applicable four-step thermodynamic framework by White and co-workers as well as a natural extension of Popot and Engelman's original two-stage model of helical MP folding to a three-stage model, where, in the new intermediate stage, the membrane-water interface acts as a holding vestibule for the translated TM helices, reconciling the interface's critical role in MP folding seen in many previous studies. Support for this model is provided by showing the stability of hydrophobic TM helices at the membrane-water interface through several microsecond long molecular dynamics simulations of five hydrophobic helical domains and a helical hairpin pre-folded from the ribosomal exit vestibule.

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.

Conformational Change from U- to I-Shape of Ion Transporters Facilitates K+ Transport across Lipid Bilayers

We have investigated, at the atomic level, the ion-fishing mechanism underlying the ion transport across membranes mediated by an artificial ion transporter composed of a hydroxyl-rich cholesterol group, a flexible alkyl chain, and a crown ether. Our results show that the transporter can spontaneously insert into the membrane and switch between the folded (U-shaped) and extended (I-shaped) conformations. The free-energy profile associated with the conformational transition indicates that compared with the U-shaped conformation of the transporter, the I-shaped one is thermodynamically more favorable. Furthermore, the free-energy profiles describing the ion translocation reveal that the transporter capturing the ion in U-shape on one side of the membrane and releasing it in I-shape on the other side constitutes a key way for the highly efficient transport of K⁺ ions. We present herewith a rigorous and rational framework to decipher the detailed ion-fishing mechanism of transmembrane ion transport with exceptionally high activity.

Untangling the complexity of membrane protein folding

Delineating the folding steps of helical-bundle membrane proteins has been a challenging task. Many questions remain unanswered, including the conformation and stability of the states populated during folding, the shape of the energy barriers between the states, and the role of lipids as a solvent in mediating the folding. Recently, theoretical frames have matured to a point that permits detailed dissection of the folding steps, and advances in experimental techniques at both single-molecule and ensemble levels enable selective modulation of specific steps for quantitative determination of the folding energy landscapes. We also discuss how lipid molecules would play an active role in shaping the folding energy landscape of membrane proteins, and how folding of multi-domain membrane proteins can be understood based on our current knowledge. We conclude this review by offering an outlook for emerging questions in the study of membrane protein folding.

Vectorial insertion of a β-helical peptide into membrane: a theoretical study on polytheonamide B

Spontaneous unidirectional, or vectorial, insertion of transmembrane peptides is a fundamental biophysical process for toxin and viral actions. Polytheonamide B (pTB) is a potent cytotoxic peptide with a β^6.3-helical structure. Previous experimental studies revealed that the pTB inserts into the membrane in a vectorial fashion and forms a channel with its single molecular length long enough to span the membrane. Also, molecular dynamics simulation studies demonstrated that the pTB is prefolded in aqueous solution. These are unique features of pTB because most of the peptide toxins form channels through oligomerization of transmembrane helices. Here, we performed all-atom molecular dynamics simulations to examine the dynamic mechanism of the vectorial insertion of pTB, providing underlying elementary processes of the membrane insertion of a prefolded single transmembrane peptide. We find that the insertion of pTB proceeds with only the local lateral compression of the membrane in three successive phases: "landing," "penetration," and "equilibration" phases. The free energy calculations using the replica-exchange umbrella sampling simulations present an energy cost of 4.3 kcal/mol at the membrane surface for the membrane insertion of pTB from bulk water. The trajectories of membrane insertion revealed that the insertion process can occur in two possible pathways, namely "trapped" and "untrapped" insertions; in some cases, pTB is trapped in the upper leaflet during the penetration phase. Our simulations demonstrated the importance of membrane anchoring by the hydrophobic N-terminal blocking group in the landing phase, leading to subsequent vectorial insertion.

Cotranslational recruitment of ribosomes in protocells recreates a translocon-independent mechanism of proteorhodopsin biogenesis

The emergence of lipid membranes and embedded proteins was essential for the evolution of cells. Translocon complexes mediate cotranslational recruitment and membrane insertion of nascent proteins, but they already contain membrane-integral proteins. Therefore, a simpler mechanism must exist, enabling spontaneous membrane integration while preventing aggregation of unchaperoned protein in the aqueous phase. Here, we used giant unilamellar vesicles encapsulating minimal translation components to systematically interrogate the requirements for insertion of the model protein proteorhodopsin (PR) – a structurally ubiquitous membrane protein. We show that the N-terminal hydrophobic domain of PR is both necessary and sufficient for cotranslational recruitment of ribosomes to the membrane and subsequent membrane insertion of PR. Insertion of N-terminally truncated PR was restored by artificially attaching ribosomes to the membrane. Our findings offer a self-sufficient protein-inherent mechanism as a possible explanation for effective membrane protein biogenesis in a “pretranslocon” era, and they offer new opportunities for generating artificial cells.

•

Generated a simple artificial cell model for membrane protein insertion

•

We identified protein-inherent control of translational targeting without chaperones

•

Ribosomes, artificially tethered to GUVs increased membrane protein insertion

Antidepressant drugs act by directly binding to TRKB neurotrophin receptors

It is unclear how binding of antidepressant drugs to their targets gives rise to the clinical antidepressant effect. We discovered that the transmembrane domain of tyrosine kinase receptor 2 (TRKB), the brain-derived neurotrophic factor (BDNF) receptor that promotes neuronal plasticity and antidepressant responses, has a cholesterol-sensing function that mediates synaptic effects of cholesterol. We then found that both typical and fast-acting antidepressants directly bind to TRKB, thereby facilitating synaptic localization of TRKB and its activation by BDNF. Extensive computational approaches including atomistic molecular dynamics simulations revealed a binding site at the transmembrane region of TRKB dimers. Mutation of the TRKB antidepressant-binding motif impaired cellular, behavioral, and plasticity-promoting responses to antidepressants in vitro and in vivo. We suggest that binding to TRKB and allosteric facilitation of BDNF signaling is the common mechanism for antidepressant action, which may explain why typical antidepressants act slowly and how molecular effects of antidepressants are translated into clinical mood recovery.

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.

In Silico Prediction of the Binding, Folding, Insertion, and Overall Stability of Membrane-Active Peptides

Membrane-active peptides (MAPs) are short-length peptides used for potential biomedical applications in diagnostic imaging of tissues, targeted drug delivery, gene delivery, and antimicrobials and antibiotics. The broad appeal of MAPs is that they are infinitely variable, relatively low cost, and biocompatible. However, experimentally characterizing the specific properties of a MAP or its many variants is a low-resolution and potentially time-consuming endeavor; molecular dynamics (MD) simulations have emerged as an invaluable tool in identifying the biophysical interactions that are fundamental to the function of MAPs. In this chapter, a step-by-step approach to discreetly model the binding, folding, and insertion of a membrane-active peptide to a model lipid bilayer using MD simulations is described. Detailed discussion is devoted to the critical aspects of running these types of simulations: prior knowledge of the system, understanding the strengths and weaknesses of molecular mechanics force fields, proper construction and equilibration of the system, realistically estimating both experimental and computational timescales, and leveraging analysis to make direct comparisons to experimental results as often as possible.

Towards a Quantitative Understanding of Protein-Lipid Bilayer Interactions at the Single Molecule Level: Opportunities and Challenges

Protein-lipid interfaces are among the most fundamental in biology. Yet applying conventional techniques to study the biophysical attributes of these systems is challenging and has left many unknowns. For example, what is the kinetic pathway and energy landscape experienced by a polypeptide chain when in close proximity to a fluid lipid bilayer? Here we review the experimental and theoretical progress we have made in addressing this question from a single molecule perspective. Some remaining impediments are also discussed.

Mechanisms of a Small Membrane-Active Antimicrobial Peptide from Hyla punctata

Thousands of antimicrobial peptides have been observed and studied in the past decades; however, their membrane-active mechanisms are ambiguous due to their dynamic structure in the cell membrane. Here, we applied both molecular dynamics (MD) simulations and biophysical experiments to study the small membrane-active antimicrobial peptide Hylaseptin P1 (HSP1), which has significant selectivity towards anionic 1-palmitoyl-2-oleoyl-sn-glycero-3-phospho-(1′-rac-glycerol) (POPG) and bacterial model membranes. HSP1 does not bind and fold onto human red blood cell model membranes, and it only binds, but does not fold, in zwitterionic 1-palmitoyl-2-oleoyl-glycero-3-phosphocholine (POPC) membranes. This suggests that the lipid chemistry and membrane rigidity are key to prevent HSP1 binding onto membranes, and the lipid headgroup charge may further promote peptide folding in the membrane. Our experiment-validated MD simulations suggest a carpet-like model mechanism for HSP1 through peptide binding, folding, aggregation, and assembly. HSP1 is shorter than the membrane thickness; therefore, the folded peptides aggregate on the surface, cross the membrane, and the oligomeric structure is supported by several surface-bound peptides in both bilayer leaflets.

Mechanistic Landscape of Membrane-Permeabilizing Peptides

Membrane permeabilizing peptides (MPPs) are as ubiquitous as the lipid bilayer membranes they act upon. Produced by all forms of life, most membrane permeabilizing peptides are used offensively or defensively against the membranes of other organisms. Just as nature has found many uses for them, translational scientists have worked for decades to design or optimize membrane permeabilizing peptides for applications in the laboratory and in the clinic ranging from antibacterial and antiviral therapy and prophylaxis to anticancer therapeutics and drug delivery. Here, we review the field of membrane permeabilizing peptides. We discuss the diversity of their sources and structures, the systems and methods used to measure their activities, and the behaviors that are observed. We discuss the fact that "mechanism" is not a discrete or a static entity for an MPP but rather the result of a heterogeneous and dynamic ensemble of structural states that vary in response to many different experimental conditions. This has led to an almost complete lack of discrete three-dimensional active structures among the thousands of known MPPs and a lack of useful or predictive sequence-structure-function relationship rules. Ultimately, we discuss how it may be more useful to think of membrane permeabilizing peptides mechanisms as broad regions of a mechanistic landscape rather than discrete molecular processes.

Chia Seed (Salvia hispanica L.) Pepsin Hydrolysates Inhibit Angiotensin-Converting Enzyme by Interacting with its Catalytic Site

High blood pressure can lead to cardiovascular diseases. The objective of this work was to obtain protein hydrolysates with antihypertensive potential from chia oil industry meal byproduct. Chia seed protein isolates (CPIs) were obtained from chia seed meal byproduct. CPI was hydrolyzed using different proteases (alcalase, pepsin, trypsin, and α-chymotrypsin) and their biological potential was evaluated using in vitro and in silico approaches. Chia seed pepsin protein hydrolysate showed the highest angiotensin-converting enzyme inhibition potential IC₅₀ of 0.128 mg/mL (P < 0.05) compared to the rest of hydrolysates. Peptide sequence LIVSPLAGRL presented the lowest predicted binding energy and highest inhibition potential (-9.5 kcal/mol) compared to other sequenced peptides and positive controls (captopril and lisinopril). Chia peptides showed potential to block angiotensin-converting enzyme by interacting with its catalytic site. Chia seed oil industry meal byproduct could be used as an inexpensive source of protein and bioactive peptides with antihypertensive potential. PRACTICAL APPLICATION: This research shows an upcycling alternative for chia oil industry byproduct. Chia meal is a rich source of protein and can be used to generate bioactive peptides with antihypertensive potential. Chia protein isolate was obtained from chia meal and hydrolyzed using different enzymes, pepsin showed the highest antihypertensive potential. Chia meal waste could be a low-cost source of protein and protein hydrolysates that could be used as a food ingredient with antihypertensive potential.

An experimentally validated approach to calculate the blood-brain barrier permeability of small molecules

Drug development for the treatment of central nervous system (CNS) diseases is extremely challenging, in large part due to the difficulty in crossing the blood-brain barrier (BBB). Here we develop and experimentally validate a new in silico method to predict quantitatively the BBB permeability for small-molecule drugs. We show accurate prediction of solute permeabilities at physiological temperature using high-temperature unbiased atomic detail molecular dynamics simulations of spontaneous drug diffusion across BBB bilayers. These simulations provide atomic detail insights into the transport mechanisms, as well as converged kinetics and thermodynamics. The method is validated computationally against physiological temperature simulations for fast-diffusing compounds, as well as experimentally by direct determination of the compound permeabilities using a transwell assay as an in vitro BBB model. The overall agreement of the predicted values with both direct simulations at physiological temperatures and experimental data is excellent. This new tool has the potential to replace current semi-empirical in silico screening and in vitro permeability measurements in CNS drug discovery.

Asymmetric Spontaneous Intercalation of Lutein into a Phospholipid Bilayer, a Computational Study

Lutein, a hydroxylated carotenoid, is a pigment synthesised by plants and bacteria. Animals are unable to synthesise lutein, nevertheless, it is present in animal tissues, where its only source is dietary intake. Both in plants and animals, carotenoids are associated mainly with membranes where they carry out important physiological functions. Their trafficking to and insertion into membranes are not well recognised due to experimental difficulties. In this paper, a computational approach is used to elucidate details of the dynamics and energetics of lutein intercalation from the water to the phospholipid bilayer phase. The dynamics is studied using molecular dynamics simulation, and the energetics using umbrella sampling. Lutein spontaneous insertion into the bilayer and translocation across it proceed via formation of hydrogen bonds between its hydroxyl groups and water and/or phospholipid oxygen atoms as well as desolvation of its polyene chain. As lutein molecule is asymmetric, its bilayer intercalation is also asymmetric. The course of events and timescale of the intercalation are different from those of helical peptides. The time of full lutein intercalation ranges from 20 to 100 ns and its final orientation is predominately vertical. Nevertheless, some lutein molecules are in the final horizontal position and some aggregate in the water phase and remain there for the whole simulation time. The highest energy barrier for the intercalation process is ~2.2 kcal/mol and the energy gain is ~18 kcal/mol. The results obtained for lutein can be applied to other xanthophylls and molecules of a similar structure.

•

Lutein as an amphiphilic molecule readily intercalates into a lipid bilayer.

•

MD simulations with dense sampling revealed details of the intercalation process.

•

Time of full lutein membrane intercalation ranges from 20 to 100 ns.

•

Lutein intercalation from the ε ring end is less probable than from the β ring end.

•

Horizontal position of lutein in the bilayer is less probable than vertical.

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.

Cooperative Nonbonded Forces Control Membrane Binding of the pH-Low Insertion Peptide pHLIP

Peptides with the ability to bind and insert into the cell membrane have immense potential in biomedical applications. pH (low) insertion peptide (pHLIP), a water-soluble polypeptide derived from helix C of bacteriorhodopsin, can insert into a membrane at acidic pH to form a stable transmembrane α-helix. The insertion process takes place in three stages: pHLIP is unstructured and soluble in water at neutral pH (state I), unstructured and bound to the surface of a membrane at neutral pH (state II), and inserted into the membrane as an α-helix at low pH (state III). Using molecular dynamics simulations, we have modeled state II of pHLIP and a fast-folding variant of pHLIP, in which each peptide is bound to a 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine bilayer surface. Our results provide strong support for recently published spectroscopic studies, namely that pHLIP preferentially binds to the bilayer surface as a function of location of anionic amino acids and that backbone dehydration occurs upon binding. Unexpectedly, we also observed several instances of segments of pHLIP folding into a stable helical turn. Our results provide a molecular level of detail that is essential to providing new insights into pHLIP function and to facilitate design of variants with improved membrane-active capabilities.

The importance of the membrane interface as the reference state for membrane protein stability

The insertion of nascent polypeptide chains into lipid bilayer membranes and the stability of membrane proteins crucially depend on the equilibrium partitioning of polypeptides. For this, the transfer of full sequences of amino-acid residues into the bilayer, rather than individual amino acids, must be understood. Earlier studies have revealed that the most likely reference state for partitioning very hydrophobic sequences is the membrane interface. We have used μs-scale simulations to calculate the interface-to-transmembrane partitioning free energies ΔG_S→TM for two hydrophobic carrier sequences in order to estimate the insertion free energy for all 20 amino acid residues when bonded to the center of a partitioning hydrophobic peptide. Our results show that prior single-residue scales likely overestimate the partitioning free energies of polypeptides. The correlation of ΔG_S→TM with experimental full-peptide translocon insertion data is high, suggesting an important role for the membrane interface in translocon-based insertion. The choice of carrier sequence greatly modulates the contribution of each single-residue mutation to the overall partitioning free energy. Our results demonstrate the importance of quantifying the observed full-peptide partitioning equilibrium, which is between membrane interface and transmembrane inserted, rather than combining individual water-to-membrane amino acid transfer free energies.

Membrane Active Peptides and Their Biophysical Characterization

In the last 20 years, an increasing number of studies have been reported on membrane active peptides. These peptides exert their biological activity by interacting with the cell membrane, either to disrupt it and lead to cell lysis or to translocate through it to deliver cargos into the cell and reach their target. Membrane active peptides are attractive alternatives to currently used pharmaceuticals and the number of antimicrobial peptides (AMPs) and peptides designed for drug and gene delivery in the drug pipeline is increasing. Here, we focus on two most prominent classes of membrane active peptides; AMPs and cell-penetrating peptides (CPPs). Antimicrobial peptides are a group of membrane active peptides that disrupt the membrane integrity or inhibit the cellular functions of bacteria, virus, and fungi. Cell penetrating peptides are another group of membrane active peptides that mainly function as cargo-carriers even though they may also show antimicrobial activity. Biophysical techniques shed light on peptide–membrane interactions at higher resolution due to the advances in optics, image processing, and computational resources. Structural investigation of membrane active peptides in the presence of the membrane provides important clues on the effect of the membrane environment on peptide conformations. Live imaging techniques allow examination of peptide action at a single cell or single molecule level. In addition to these experimental biophysical techniques, molecular dynamics simulations provide clues on the peptide–lipid interactions and dynamics of the cell entry process at atomic detail. In this review, we summarize the recent advances in experimental and computational investigation of membrane active peptides with particular emphasis on two amphipathic membrane active peptides, the AMP melittin and the CPP pVEC.

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.

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.

Modeling of the Passive Permeation of Mercury and Methylmercury Complexes Through a Bacterial Cytoplasmic Membrane

Cellular uptake and export are important steps in the biotransformation of mercury (Hg) by microorganisms. However, the mechanisms of transport across biological membranes remain unclear. Membrane-bound transporters are known to be relevant, but passive permeation may also be involved. Inorganic Hg^II and methylmercury ([CH₃Hg^II]⁺) are commonly complexed with thiolate ligands. Here, we have performed extensive molecular dynamics simulations of the passive permeation of Hg^II and [CH₃Hg^II]⁺ complexes with thiolate ligands through a model bacterial cytoplasmic membrane. We find that the differences in free energy between the individual complexes in bulk water and at their most favorable position within the membrane are ∼2 kcal mol^-1. We provide a detailed description of the molecular interactions that drive the membrane crossing process. Favorable interactions with carbonyl and tail groups of phospholipids stabilize Hg-containing solutes in the tail-head interface region of the membrane. The calculated permeability coefficients for the neutral compounds CH₃S-Hg^II-SCH₃ and CH₃Hg^II-SCH₃ are on the order of 10^-5 cm s^-1. We conclude that small, nonionized Hg-containing species can permeate readily through cytoplasmic membranes.

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.

Simulations of Membrane-Disrupting Peptides I: Alamethicin Pore Stability and Spontaneous Insertion

An all-atom molecular dynamics simulation of the archetype barrel-stave alamethicin (alm) pore in a 1,2-dioleoyl-sn-glycero-3-phosphocholine bilayer at 313 K indicates that ∼7 μs is required for equilibration of a preformed 6-peptide pore; the pore remains stable for the duration of the remaining 7 μs of the trajectory, and the structure factors agree well with experiment. A 5 μs simulation of 10 surface-bound alm peptides shows significant peptide unfolding and some unbinding, but no insertion. Simulations at 363 and 413 K with a -0.2 V electric field yield peptide insertion in 1 μs. Insertion is initiated by the folding of residues 3-11 into an α-helix, and mediated by membrane water or by previously inserted peptides. The stability of five alm pore peptides at 413 K with a -0.2 V electric field demonstrates a significant preference for a transmembrane orientation. Hence, and in contrast to the cationic antimicrobial peptide described in the following article, alm shows a strong preference for the inserted over the surface-bound state.

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 Modelling, Dynamics, and Docking of Membrane Proteins

Computational tools and techniques are now most popular and promising to progress the research at rapid rate. Molecular modelling studies contribute their maximum role in wide variety of disciplines especially in proteomics and drug discovery strategies. Molecular dynamics and molecular docking algorithms are now became an essential part in daily research activities of every laboratory throughout the world. These strategies are now well established and standardised to study any specific protein of interest and drug molecule. But still there exist considerable drawbacks in a special concern with membrane proteins as the presently available tools and methods cannot be applied directly to them. Modelling, dynamics and docking studies of membrane proteins need a special care and attention as several challenges are to be crossed with an intensive care to produce a reliable result. This chapter is aimed to discuss such challenges and solutions to handle membrane proteins.

Spontaneous formation of structurally diverse membrane channel architectures from a single antimicrobial peptide

Many antimicrobial peptides (AMPs) selectively target and form pores in microbial membranes. However, the mechanisms of membrane targeting, pore formation and function remain elusive. Here we report an experimentally guided unbiased simulation methodology that yields the mechanism of spontaneous pore assembly for the AMP maculatin at atomic resolution. Rather than a single pore, maculatin forms an ensemble of structurally diverse temporarily functional low-oligomeric pores, which mimic integral membrane protein channels in structure. These pores continuously form and dissociate in the membrane. Membrane permeabilization is dominated by hexa-, hepta- and octamers, which conduct water, ions and small dyes. Pores form by consecutive addition of individual helices to a transmembrane helix or helix bundle, in contrast to current poration models. The diversity of the pore architectures—formed by a single sequence—may be a key feature in preventing bacterial resistance and could explain why sequence–function relationships in AMPs remain elusive.

Antimicrobial peptides (AMPs) selectively form pores in microbial membranes in process not fully understood. Here the authors use experimentally guided molecular dynamics to study maculatin pore formation, showing how this AMP assembles into transient and structurally diverse oligomeric pores in cell membranes.

Coarse-Grained Simulations of Membrane Insertion and Folding of Small Helical Proteins Using the CABS Model

The CABS coarse-grained model is a well-established tool for modeling globular proteins (predicting their structure, dynamics, and interactions). Here we introduce an extension of the CABS representation and force field (CABS-membrane) to the modeling of the effect of the biological membrane environment on the structure of membrane proteins. We validate the CABS-membrane model in folding simulations of 10 short helical membrane proteins not using any knowledge about their structure. The simulations start from random protein conformations placed outside the membrane environment and allow for full flexibility of the modeled proteins during their spontaneous insertion into the membrane. In the resulting trajectories, we have found models close to the experimental membrane structures. We also attempted to select the correctly folded models using simple filtering followed by structural clustering combined with reconstruction to the all-atom representation and all-atom scoring. The CABS-membrane model is a promising approach for further development toward modeling of large protein-membrane systems.

Submillisecond Dynamics of Mastoparan X Insertion into Lipid Membranes

The mechanism of protein insertion into a lipid bilayer is poorly understood because the kinetics of this process is difficult to measure. We developed a new approach to study insertion of the antimicrobial peptide Mastoparan X into zwitterionic lipid vesicles, using a laser-induced temperature-jump to initiate insertion on the microsecond time scale and infrared and fluorescence spectroscopies to follow the kinetics. Infrared probes the desolvation of the peptide backbone and yields biphasic kinetics with relaxation lifetimes of 12 and 117 μs, whereas fluorescence probes the intrinsic tryptophan residue located near the N-terminus and yields a single exponential phase with a lifetime of 440 μs. Arrhenius analysis of the temperature-dependent rates yields an activation energy for insertion of 96 kJ/mol. These results demonstrate the complexity of the insertion process and provide mechanistic insight into the interplay between peptides and the lipid bilayer required for peptide transport across cellular membranes.

Multiscale Simulations Suggest a Mechanism for the Association of the Dok7 PH Domain with PIP-Containing Membranes

Dok7 is a peripheral membrane protein that is associated with the MuSK receptor tyrosine kinase. Formation of the Dok7/MuSK/membrane complex is required for the activation of MuSK. This is a key step in the complex exchange of signals between neuron and muscle, which lead to neuromuscular junction formation, dysfunction of which is associated with congenital myasthenic syndromes. The Dok7 structure consists of a Pleckstrin Homology (PH) domain and a Phosphotyrosine Binding (PTB) domain. The mechanism of the Dok7 association with the membrane remains largely unknown. Using multi-scale molecular dynamics simulations we have explored the formation of the Dok7 PH/membrane complex. Our simulations indicate that the PH domain of Dok7 associates with membranes containing phosphatidylinositol phosphates (PIPs) via interactions of the β1/β2, β3/β4, and β5/β6 loops, which together form a positively charged surface on the PH domain and interact with the negatively charged headgroups of PIP molecules. The initial encounter of the Dok7 PH domain is followed by formation of additional interactions with the lipid bilayer, and especially with PIP molecules, which stabilizes the Dok7 PH/membrane complex. We have quantified the binding of the PH domain to the model bilayers by calculating a density landscape for protein/membrane interactions. Detailed analysis of the PH/PIP interactions reveal both a canonical and an atypical site to be occupied by the anionic lipid. PH domain binding leads to local clustering of PIP molecules in the bilayer. Association of the Dok7 PH domain with PIP lipids is therefore seen as a key step in localization of Dok7 to the membrane and formation of a complex with MuSK.

Neuromuscular junction formation and maintenance is an essential biological process, the disruption of which leads to congenital myasthenic syndromes and premature death. Dok7 is a key member in formation, maintenance and signaling in neuromuscular junctions. Dok7 is a peripheral membrane protein that is necessary for full activation of the receptor tyrosine kinase MuSK, a receptor tyrosine kinase residing in the postsynaptic membrane. The structure of Dok7 consists of both a PH domain and a PTB domain. The interaction of Dok7 with cell membranes is not well understood. Here, we use molecular simulations to show that the Dok7 PH domain preferentially binds to PIP lipid molecules when associating with a membrane. Dok7 interacts with the bilayer in both a canonical binding mode and an alternative binding mode. Our simulations also demonstrate the presence of both a canonical and an atypical binding site for PIPs in agreement with recent crystallographic studies of the ASAP1 PH domain. Analysis of density landscapes for the interaction of the Dok7 PH domain with PIP-containing lipid bilayers is also able to identify both canonical and alternative binding modes. Our improved understanding of how Dok7 interacts with a membrane is key to examining Dok7/MuSK signaling.

Interactions of Pleckstrin Homology Domains with Membranes: Adding Back the Bilayer via High-Throughput Molecular Dynamics

A molecular simulation pipeline for determining the mode of interaction of pleckstrin homology (PH) domains with phosphatidylinositol phosphate (PIP)-containing lipid bilayers is presented. We evaluate our methodology for the GRP1 PH domain via comparison with structural and biophysical data. Coarse-grained simulations yield a 2D density landscape for PH/membrane interactions alongside residue contact profiles. Predictions of the membrane localization and interactions of 13 PH domains reveal canonical, non-canonical, and dual PIP-binding sites on the proteins. Thus, the PH domains associate with the PIP molecules in the membrane via a highly positively charged loop. Some PH domains exhibit modes of interaction with PIP-containing membranes additional to this canonical binding mode. All 13 PH domains cause a degree of local clustering of PIP molecules upon binding to the membrane. This provides a global picture of PH domain interactions with membranes. The high-throughput approach could be extended to other families of peripheral membrane proteins.

Mechanisms of Membrane Pore Formation by Amyloidogenic Peptides in Amyotrophic Lateral Sclerosis

Using unbiased atomic-detailed molecular dynamics simulations, the C-terminal fragments of TDP-43 are observed to aggregate and form disordered-toroidal pores in a lipid bilayer. Cytotoxicity of TDP-43 may be inferred from the observation that the membrane pores catalyze lipid flip-flop between bilayer leaflets and conduct water at high rates.