Identification of electroporation sites in the complex lipid organization of the plasma membrane

  1. Lea Rems  Is a corresponding author
  2. Xinru Tang
  3. Fangwei Zhao
  4. Sergio Pérez-Conesa
  5. Ilaria Testa
  6. Lucie Delemotte  Is a corresponding author
  1. University of Lubljana, Slovenia
  2. KTH Royal Institute of Technology, Sweden

Abstract

The plasma membrane of a biological cell is a complex assembly of lipids and membrane proteins, which tightly regulate transmembrane transport. When a cell is exposed to strong electric field, the membrane integrity becomes transiently disrupted by formation of transmembrane pores. This phenomenon termed electroporation is already utilized in many rapidly developing applications in medicine including gene therapy, cancer treatment, and treatment of cardiac arrythmias. However, the molecular mechanisms of electroporation are not yet sufficiently well understood; in particular, it is unclear where exactly pores form in the complex organization of the plasma membrane. In this study we combine coarse-grained molecular dynamics simulations, machine learning methods, and Bayesian survival analysis to identify how formation of pores depends on the local lipid organization. We show that pores do not form homogeneously across the membrane, but colocalize with domains that have specific features, the most important being high density of polyunsaturated lipids. We further show that knowing the lipid organization is sufficient to reliably predict poration sites with machine learning. Additionally, by analysing poration kinetics with Bayesian survival analysis we show that poration does not depend solely on local lipid arrangement, but also on membrane mechanical properties and the polarity of the electric field. Finally, we discuss how the combination of atomistic and coarse-grained molecular dynamics simulations, machine learning methods, and Bayesian survival analysis can guide the design of future experiments and help us to develop an accurate description of plasma membrane electroporation on the whole-cell level. Achieving this will allow us to shift the optimization of electroporation applications from blind trial-and-error approaches to mechanistic-driven design.

Data availability

The current manuscript is a computational study, so no data have been generated for this manuscript. Analysis codes are available at https://github.com/learems/Electroporation-CGmem-MemSurfer, https://github.com/learems/Electroporation-CGmem-MachineLearning, https://docs.pymc.io/en/v3/pymc-examples/examples/survival_analysis/survival_analysis.html and https://github.com/sperezconesa/electroporation_modeling

The following data sets were generated

Article and author information

Author details

  1. Lea Rems

    University of Lubljana, Ljubljana, Slovenia
    For correspondence
    lea.rems@fe.uni-lj.si
    Competing interests
    No competing interests declared.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0001-7470-4367
  2. Xinru Tang

    Department of Applied Physics, KTH Royal Institute of Technology, Stockholm, Sweden
    Competing interests
    No competing interests declared.
  3. Fangwei Zhao

    Department of Applied Physics, KTH Royal Institute of Technology, Stockholm, Sweden
    Competing interests
    No competing interests declared.
  4. Sergio Pérez-Conesa

    Department of Applied Physics, KTH Royal Institute of Technology, Stockholm, Sweden
    Competing interests
    No competing interests declared.
  5. Ilaria Testa

    Department of Applied Physics, KTH Royal Institute of Technology, Stockholm, Sweden
    Competing interests
    Ilaria Testa, Reviewing editor, eLife.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0003-4005-4997
  6. Lucie Delemotte

    Department of Applied Physics, KTH Royal Institute of Technology, Stockholm, Sweden
    For correspondence
    lucie.delemotte@scilifelab.se
    Competing interests
    Lucie Delemotte, Reviewing editor, eLife.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-0828-3899

Funding

Science for Life Laboratory

  • Ilaria Testa
  • Lucie Delemotte

Vetenskapsrådet (2018-04905)

  • Lucie Delemotte

Gustafssons Stiftelse

  • Lucie Delemotte

Horizon 2020 Marie Skłodowska-Curie (893077)

  • Lea Rems

Slovenian Research Agency (J2-2503)

  • Lea Rems

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Copyright

© 2022, Rems et al.

This article is distributed under the terms of the Creative Commons Attribution License permitting unrestricted use and redistribution provided that the original author and source are credited.

Metrics

  • 2,221
    views
  • 330
    downloads
  • 16
    citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Lea Rems
  2. Xinru Tang
  3. Fangwei Zhao
  4. Sergio Pérez-Conesa
  5. Ilaria Testa
  6. Lucie Delemotte
(2022)
Identification of electroporation sites in the complex lipid organization of the plasma membrane
eLife 11:e74773.
https://doi.org/10.7554/eLife.74773

Share this article

https://doi.org/10.7554/eLife.74773

Further reading

    1. Biochemistry and Chemical Biology
    2. Structural Biology and Molecular Biophysics
    Cristina Paissoni, Sarita Puri ... Carlo Camilloni
    Research Article

    Both immunoglobulin light-chain (LC) amyloidosis (AL) and multiple myeloma (MM) share the overproduction of a clonal LC. However, while LCs in MM remain soluble in circulation, AL LCs misfold into toxic-soluble species and amyloid fibrils that accumulate in organs, leading to distinct clinical manifestations. The significant sequence variability of LCs has hindered the understanding of the mechanisms driving LC aggregation. Nevertheless, emerging biochemical properties, including dimer stability, conformational dynamics, and proteolysis susceptibility, distinguish AL LCs from those in MM under native conditions. This study aimed to identify a2 conformational fingerprint distinguishing AL from MM LCs. Using small-angle X-ray scattering (SAXS) under native conditions, we analyzed four AL and two MM LCs. We observed that AL LCs exhibited a slightly larger radius of gyration and greater deviations from X-ray crystallography-determined or predicted structures, reflecting enhanced conformational dynamics. SAXS data, integrated with molecular dynamics simulations, revealed a conformational ensemble where LCs adopt multiple states, with variable and constant domains either bent or straight. AL LCs displayed a distinct, low-populated, straight conformation (termed H state), which maximized solvent accessibility at the interface between constant and variable domains. Hydrogen-deuterium exchange mass spectrometry experimentally validated this H state. These findings reconcile diverse experimental observations and provide a precise structural target for future drug design efforts.

    1. Structural Biology and Molecular Biophysics
    Kingsley Y Wu, Ta I Hung, Chia-en A Chang
    Research Article

    PROteolysis TArgeting Chimeras (PROTACs) are small molecules that induce target protein degradation via the ubiquitin-proteasome system. PROTACs recruit the target protein and E3 ligase; a critical first step is forming a ternary complex. However, while the formation of a ternary complex is crucial, it may not always guarantee successful protein degradation. The dynamics of the PROTAC-induced degradation complex play a key role in ubiquitination and subsequent degradation. In this study, we computationally modelled protein complex structures and dynamics associated with a series of PROTACs featuring different linkers to investigate why these PROTACs, all of which formed ternary complexes with Cereblon (CRBN) E3 ligase and the target protein bromodomain-containing protein 4 (BRD4BD1), exhibited varying degrees of degradation potency. We constructed the degradation machinery complexes with Culling-Ring Ligase 4A (CRL4A) E3 ligase scaffolds. Through atomistic molecular dynamics simulations, we illustrated how PROTAC-dependent protein dynamics facilitating the arrangement of surface lysine residues of BRD4BD1 into the catalytic pocket of E2/ubiquitin cascade for ubiquitination. Despite featuring identical warheads in this PROTAC series, the linkers were found to affect the residue-interaction networks, and thus governing the essential motions of the entire degradation machine for ubiquitination. These findings offer a structural dynamic perspective on ligand-induced protein degradation, providing insights to guide future PROTAC design endeavors.