Reversible electroporation (RE) is a well-established methodology for delivering a variety of normally membrane-impermeable compounds and molecules into cells1,2,3. Here, pulsed electric fields (PEFs) are used to achieve a critical electric field intensity that induces transient, localized nanopores along portions of the cell membrane. RE has been used to deliver a variety of molecules and compounds ranging from dyes and nucleic acids to chemotherapeutics and other drugs4,5,6,7,8,9,10. However, if the electric field intensity is increased sufficiently, it can meet a higher critical threshold, which results in a rapid decline in cell viability1,2,11,12. This process is known as irreversible electroporation (IRE) and has been studied for soft-tissue ablation, including as a treatment for otherwise inoperable tumors and cardiac arrhythmias11,12,13,14,15,16. To understand treatment outcomes and inform future treatments, quantification of these RE and IRE thresholds is common practice in electroporation research12,17,18,19,20. However, RE and IRE thresholds vary depending on many factors, including cell type, applied voltage, pulse duration, number of pulses, and pulse delivery rate1,17,18,21,22.
Electroporation protocols have a number of parameters that result in varied RE and IRE thresholds1,17,18,21,22. Optimal treatment voltages and waveforms have long been investigated. Most commercially available electroporation technologies utilize protocols with monopolar PEFs on the order of 100 µs-10 ms23. While these protocols achieve clinically relevant outcomes, they induce intense muscle contractions in vivo24,25,26. Because of this, shorter, bipolar microsecond pulses, which alleviate this stimulation, have become a recent topic of interest1,11,24. However, these pulses introduce new factors which affect RE and IRE thresholds, including pulse symmetry, burst duration, and temperature-dependence1,17,22,27. Additionally, cell type and tissue properties have also been shown to significantly affect both RE and IRE outcomes19,28,29. As the field continues to evolve, there is a need to further characterize how these factors and their interactions with each other and other possible environmental factors affect RE and IRE thresholds, making it increasingly necessary to improve experimental efficiency to enable higher-throughput experimentation.
Typical in vitro electroporation experiments test protocols by treating cells in a cuvette with parallel conductive plates, which results in a uniform electric field distribution within the cuvette during treatment19,30,31,32,33,34. As aforementioned, it is necessary to identify a protocol's RE and/or IRE thresholds to identify an optimal applied voltage for a given treatment. For this to be achieved in cuvette studies, multiple cuvettes must be treated, each with discrete applied voltages. This requires a significant number of samples that must be prepared, treated, and analyzed separately, limiting experimental throughput. To improve this process, both 2D35,36 and 3D13,20,27,37,38,39 culture models have been developed to evaluate a continuum of field strengths in a single sample by using electrodes that create non-uniform field distributions13,20,27,37,38,39. RE and IRE thresholds are then derived from computational models of the electric field distribution13,20,27,37,38. This greatly reduces the number of samples needed, enabling a more high-throughput approach to threshold quantification.
This article aims to demonstrate this approach for RE and IRE field quantification in a 3D in vitro model to demonstrate its potential as a means for efficient evaluation of electroporation protocols across several environments (e.g., cell types, treatment parameters, molecules to be delivered).