Electronic Emulator of Biological Tissue as an Electrical Load during Electroporation

Electroporation is an emerging technology, with great potential in many different medical and biotechnological applications, food engineering and biomass processing. Large variations of biological load characteristics, however, represent a great challenge in electroporator design, which results in different solutions. Because a clinical electroporator is a medical device, it must comply with medical device regulative and standards. However, none of the existing standards directly address the operation or electroporator’s performance requirements. In order to evaluate clinical, laboratory and prototype electroporation devices during the development process, or to evaluate their final performance considering at least from the perspective of output pulse parameters, we present a case study on the design of an electronic emulator of biological tissue as an electrical load during electroporation. The proposed electronic load emulator is a proof of concept, which enables constant and sustainable testing and unbiased comparison of different electroporators’ operations. We developed an analog electrical circuit that has equivalent impedance to the beef liver tissue in combination with needle electrodes, during high voltage pulse delivery and/or electroporation. Current and voltage measurements during electroporation of beef liver tissue ex vivo, were analyzed and parametrized to define the analog circuit equation. An equivalent circuit was simulated, built and validated. The proposed concept of an electronic load emulator can be used for “classical” electroporator (i.e., not nanosecond) performance evaluation and comparison of their operation. Additionally, it facilitates standard implementation regarding the testing protocol and enables quality assurance.


Introduction
Electroporation is an emerging platform technology, with great potential [1][2][3]. One of the most successful electroporation-based applications in medicine is an anti-tumor therapy, called electrochemotherapy (ECT). Through the permeabilized cell membrane transport of chemotherapeutic drug, bleomycin or cisplatin is increased and considerably higher cytotoxicity is reached [4][5][6][7]. According to Standard Operating Procedure for ECT [8,9] in clinical practice, mainly eight square wave pulses often reported as amplitude to distance ratio in range from 1000 to 1500 V/cm, 100 µs pulse duration each and repetition frequency of 1 Hz or 5 kHz are applied to the tissue. A medical device that enables electroporation is called an electroporator; it is a high voltage pulse generator that generates high voltage pulses of a specific duration, pulse repetition frequency and shape [10,11]. From the point of view of an electroporator, the load is a tissue defined by the electrodes [12].
Electric properties of biological tissue vary considerably, tissue has an inhomogeneous composition [13][14][15] and different electrode type/geometry (different distances between the electrodes) are used. In addition due to electroporation [16][17][18] tissue conductivity changes during pulse delivery, which results in a non-linear current response [16,[19][20][21]. However the pre-set train of pulses should be generated and delivered to the load, for successful electroporation. This present a great challenge in electroporator design which results in different solutions. Operation errors mainly occur with low resistant loads, or when the pulse repetition rate is in the higher half of the device's operation range. Namely some pulse generators cannot deliver sufficiently high currents or energies on the output [11]. Electroporators are however not only pulse generators, but they also have built-in protection systems for current and energy limitation and systems for ensuring the effectiveness of the therapy (for protection of patients, operators, and electronics), the operation at boundary condition of safe operation area may deviate from the expectations. Implementation of protection systems is not necessarily always straight forward, because even at normal operation, high displacement current spikes occur at the beginning of the pulse (due to capacitive nature of biological tissue as the load) and ringing is sometimes present at pulse switching.
For more than a decade our group is evaluating electroporation devices [11,22,23] and currently, a rise in electroporation device providers, can be detected on the market. Unfortunately the transitions between devices from different manufacturers, or comparison of results obtained by different researchers using different electroporators, is becoming more and more difficult or even impossible, due to the pre-programmed electroporation procedures. The researchers do not even know and therefore also do not report basic pulse parameters (e.g., pulse shape, repetition rate, and voltage amplitude [24,25]).
The most stressful conditions for an electroporator are short-circuit or discharge between electrodes during the pulse delivery, which can be tested with short-circuiting or very closely placed electrodes. However until now, for the evaluation of an electroporator performance within the safe operation area mainly a resistor or more advanced, a resistor with a capacitance in parallel were used. But those evaluations are considered over-simplified, especially for electroporators operating near their current and energy limits, and for electroporators with over-current protection circuit. Some of the clinically approved electroporators significantly drop the voltage during the delivery of a pulse train due to increased demand for current and energy. And some electroporators even stop delivering electroporation pulses during the delivery of a pulse train, i.e., in the middle of the treatment to protect the device against over-currents. This is especially problematic in clinical treatments where in such a case the patient receives only part of the treatment. Simpler load models as resistive capacitive model cannot simulate increased demand for current and energy, and activation of protection against over-currents during the delivery of a train of pulses. Since in the case of ECT a typical load is tissue, real tissue or its model should be used [26]. Currently, however, there are no available models on the market, which incorporates inhomogeneous characteristics of tissue or are too expensive, or for single-use and are thus inadequate, therefore some researchers rather use for example, a potato [27,28]. This might be a better solution than a resistor, because at least a conductivity increase (due to electroporation) and heating are present. Only a load that raises the current from pulse to pulse represent real-life situation. We think that many problems concerning electroporation devices have evolved due to the lack of standards and evaluation protocols. These are the reasons why some times therapy may not be successful and why we think that an electronic emulator is needed to thoroughly test electroporation devices before usage.
In the scope of this paper we present an idea of an analog electrical circuit that emulates biological tissue, specifically beef liver tissue in combination with needle electrodes, undergoing high voltage pulse delivery and/or electroporation. Current and voltage measurements during electroporation of liver tissue ex vivo, were analyzed and parameterized in a way to fit the proposed analog circuit equation. An equivalent circuit was then simulated, built and validated as a proof of concept. The developed equivalent electronic emulator can be used for electroporators performance evaluation and comparison of the operation between different electroporators, or of an electroporator over its lifetime. The developed device emulates small volume beef liver tissue in combination with needle electrodes. The proposed process/concept can be however applied to any electroporation measurements and specific emulators can be developed.

Methods and Materials
For the development of an electronic emulation circuit, high-quality measurements of current and voltage, during pulse application to the biological tissue are needed. We used ex vivo measurements on beef liver tissue, for which accurate and high sample rate ex vivo measurements were available [19]. Langus et al. [19] introduced time-dependent effects into a finite element model developed specifically for electroporation, which was able to predict the time evolution of electric current, within a 5% error. Briefly ex vivo measurements on beef liver samples, were obtained, nine different sequences of electric pulses were delivered to the tissue, by commercially available needle electrodes (IGEA, Carpi, Italy) with a 10 mm of active part, a diameter of 1.2 mm and the distance between the centers of the electrodes being 10 mm. Applied voltage (with electroporator: i.e., ELECTRO cell S20, β-teh, Leroy, France) and current were monitored with an oscilloscope (WavePro 7300 A, LeCroy, Milpitas, CA, USA), by a current probe (AP015, LeCroy, USA) and a high-voltage differential voltage probe (ADP305, LeCroy, USA) and used for the development of the model. Current and voltage measurements of eight pulses in a sequence, 1000 V, 750 V, or 500 V voltage amplitude with a pulse duration of 100 µs, two different repetition frequencies 1 Hz and 4717 Hz, were re-processed and analyzed in order to model and emulate electroporation load with an analog electric circuit. Measurement data were imported into Matlab R2018b (MathWorks, Natick, MA, USA) computing environment, where values of optimal circuit elements were calculated and then the circuit was simulated in PSpice AD Lite (OrCAD, Cadence, CA, USA). When pulsed electric fields (PEF) are applied to a biological tissue, current response has a typical time course which consists of a fast, high capacitive current, followed by a slight amplitude rise due to tissue electroporation, or more or less square wave pulse, if tissue was not yet electroporated [17,29,30] (Figure 1a).

Mathematical Model
A typical current time course, during PEF application to tissue, is shown in the figure (Figure 1a), where dynamics associated with pulse rise and fall time have not been included (indicated with the red dashed line). The dynamic of measured current, during the pulse (pulse rise and fall time are not included) can be roughly described as a sum of a square wave pulse and two exponential functions, as presented with the equation (Equation (1)) ( Figure 1b), where t represents time and A, B, C, D, and E are abstract parameters that are in the circuit model substituted by applied voltage and load parameters. Such a form of the equation has been chosen, because it represents a resistive, capacitive response and increase of conductivity due to electroporation, but also because it can be easily converted to an equation, which defines an analog electric circuit: where T pulse is equal to the duration of the pulse. The first element of the equation, a constant A, represents a real (ohmic) part of the current, which is a result of resistive tissue component. But as current in tissue is also reactive, two additional elements are added. The second element B × e −t/C , an exponential decay function, represents a capacitive current spike and can be modeled as the capacitor discharge. The third element D × (1 − e −t/E ), an exponential rise function, illustrates a conductivity growth during pulse delivery, which is a result of tissue electroporation (pore formation) and tissue heating [16][17][18].  (1)).

Circuit Model
The mathematical model (Equation (1)) can be easily substituted by a circuit model (Scheme 1) with three parallel branches, namely, resistive, capacitive and inductive, which models increase of conductivity due to electroporation (Equation (2)). When the applied voltage is lower than the electroporation threshold voltage, the tissue is not electroporated and therefore the conductivity is roughly constant. In the equivalent circuit (Scheme 1) a Zener diode with breakdown voltage equal to the electroporation threshold voltage, enables or disables an inductive/electroporation branch. Additionally, at the end of each applied pulse, the stored energy in the capacitor and inductor will be released and will cause the transient exponential decay current on the trailing edge of the current waveform, resulting in a short negative value current spike. The negative spike is also present in the case of biological tissue, but we only modeled the dynamics of the pulse to facilitate calculations, therefore we neglected the rise, fall time and post pulse dynamics.

Single Monopolar Pulse Model
For a single positive monpolar pulse the following circuit is proposed (Scheme 1).

Scheme 1.
Single monopolar pulse circuit emulator: equivalent circuit model of tissue during electroporation. Where V ep represent applied electroporation pulse, with pulse voltage amplitude U amp and U th electroporation threshold voltage.
The values of the circuit elements are calculated using the following equations (Equations (2) and (3)): where U amp represents pulse voltage amplitude of applied electroporation pulse (V ep in Scheme 1) and U th electroporation threshold voltage. For the threshold voltage determination measurements at different applied voltages are required. While all the other values can be determined directly from the mathematical model: This solution is accurate only in case of one pulse application, or when pulse repetition frequency is low (pause between pulses is longer than 1 ms [31]) e.g., in our case 1 Hz, meaning all the applied pulses in one sequence have approximately the same current time course.
For numerical model, all measurement data (obtained and described earlier [19]), were imported into Matlab R2018b (MathWorks, Natick, MA, USA) computing environment, where values of optimal circuit elements were calculated. Due to a big number of undefined variables (R 1 , R 2 , R 3 , R 4 , C 1 , L 1 , U th ) a deterministic optimisation method was used. A flow chart ( Figure 2) visualises a calculation process. Five high-quality measurements of applied current and voltage for each preset voltage 500 V, 750 V and 1000 V (pulse duration: 100 µs and repetition frequency: 4.717 kHz), were averaged and down-sampled for further analyses. From average voltage and current signals, the active parts were cutout (values in the range from the maximal value of current to the FWHM minus 1 µs). Averaged current signals were split into two parts, a capacitive current spike part and electroporation part (the minimal value in the cutout signal was set as a split point). In order to facilitate calculations, signals were decimated. The capacitive current spike part was due to its short duration-dynamics decimated by 10 (sample rate was approximately 50 MHz) and electroporation part by 70 (sample rate was approximately 7 MHz). Those two values (decimation by 10 and 70) were chosen in a way that half of the values in the new decimated signal represented a capacitive current spike part and the other half an inductive/electroporation part of the signal. At first initial values for final least square error (LSE) fit, were determined. The dynamic conductivity was calculated from applied averaged and cut-out current and voltage signals, only for the first pulse in sequence. With the Matlab's curve fitting toolbox, optimal parameters (A, B, C, D, E) for proposed equation (Equation (1)) were defined, for each applied voltage separately. Than optimal electronic component values (R 1 , R 2 , R 3 , τ 1 , τ 2 ) were calculated for each applied voltage in accordance with the equation (Equation (4)), assuming the threshold voltage is equal to 200 V (assumption based on the slope of the current dynamic of electroporation part). Due to a large number of variables and three applied voltages, each variable range was predefined. For resistive components (R 1 , R 2 , R 3 ) and time constants (τ 1 , τ 2 ), fifteen linearly distributed values between a minimal and maximal proposed optimal values from individual fit were used. Then LSE regression was done for all three applied voltages together. While the threshold voltage was now limited between 200 V and 400 V, because in Langus et al. [19] optimal fitting maximum and minimum electric field magnitude limit were set to 200 V/cm and 400 V/cm. In this way all values of elements for single monopolar pulse circuit (Scheme 1) were determined.

Multiple Monopolar Pulse Model
With a higher repetition frequencies, in our case 4717 Hz, it was observed that the current reaches around 75% of the final amplitude at the first pulse and percentage than exponentially rises with pulses, while all pulses in sequence have the same shape and amplitude in case of 1 Hz repetition frequency [19]. Already Pavlin et al. showed that for pulse repetition frequencies above 1 kHz a significant change in pulse shape occurs, meaning the following pulses start at higher initial value [31]. Heat generation, diffusion of intracellular content and pore lifetime could explain these phenomena [31][32][33]. Therefore in parallel to inductor L 1 an additional resistor R 4 and diode in series are added (Scheme 2). A diode emulates an "electroporation" memory while resistor R 4 emulates the loss of memory. For the first pulse in a pulse sequence with higher repetition frequency, the situation is the same as for first/single pulse or pulses with lower repetition frequency (Scheme 1). The electroporation memory is not yet present. During the pulse, the inductor L 1 is magnetized and right after the pulse, during the pause between pulses, the L 1 current in the model regarded as I M0 starts to flow through the new resistor R 4 and the inductor is in a demagnetization phase until the next pulse arrival. I E is still present when the next pulse is generated because the pause between the pulses is shorter than inductor demagnetizing time. Therefore the I E increases, it is the sum of inductor magnetize current and I mem memory current. The inductor current for second and all other pulses is calculated in according to (Equations (5) and (6)): I mem (N) = I M0 (N) × e −T pause /τ 3 (6) in which, N is a pulse number, I M0 represents a current at the end of the previous pulse, T pause is equal to the duration of the pause between pulses and T per in pulse period (T per = T pulse + T pause ).
Scheme 2. Additional in series connected resistor R 4 and diode in parallel to the inductor L 1 , emulate an "electroporation" memory, which enables that the following pulse starts at the higher initial value.
All the calculated optimal component values calculated in previous section (R 1 , R 2 , R 3 , C 1 , L 1 ) are the same also in case of multiple monopolar pulse model, only threshold voltage was optimized again and R 4 which enables memory current was again defined by using LSE regression. In this case, the whole sequence, all eight pulses were, included in the calculation. Because only one variable R 4 , was left undefined, a range of potential R 4 values was extended, meaning a lot higher number of potential values were included in the regression. Final circuit emulator (Scheme 3) is based on equations (Equations (7)-(10)), which are only valid for the time of the pulse ((N − 1) × T per < t < (N − 1) × T per + T pulse ), and well describe the behaviour of our specific biological tissue.
Final load current for higher repetition frequencies is described with discrete equation (Equation (7)), where N is the pulse number, T per period (pulse duration + pause). When the applied voltage is lower than the electroporation threshold voltage, the load current consists only of the ohmic and capacitive components. However when the applied voltage is higher than the electroporation threshold voltage, additionally to an ohmic and capacitive component also the inductive/electroporation current is added, which emulates conductivity rise during electroporation (Equation (7)). In the second circuit (Scheme 3) an additional memory current, which is pulse number dependent, is added to the inductive component. At the first pulse, the memory current is equal to zero because the inductor L 1 is not yet magnetized. Then with each successive pulse the memory current is higher until it reaches the plateau (Equations (8) and (9)), due to inductor magnetization and demagnetization. From the time constant also the last missing components (C 1 , L 1 and R 4 ) are calculated and thus the circuit is fully defined (Equation (10)). The proposed solution can be also used for lower frequencies in our case 1 Hz repetition rate, where all the pulses in sequence have the same dynamic. Because the time constant (329 µs) of the circuit which adds the memory current is significantly smaller than pulse repetition rate (1 Hz − 1 s) the memory current can be in this case neglected.
Final electronic emulator of biological tissue as an electrical load during electroporation.

Electronic Load Emulator
Since electronic components with exactly the same values as defined in the numerical model are not available on the market, the prototype's components that have the closest values, were used. The circuit with commercially available components values was first simulated in PSpice AD Lite (OrCAD, Cadence, CA, USA). The load current was compared in Matlab to the measured current and an emulator was built and evaluated. The electric pulses, with the same pulse settings and with the same electroporation device (i.e., ELECTRO cell S20, β-teh, Leroy, France) as in Langus et al. [19], were applied to the electronic load emulator. For the measurements of prototype emulator, oscilloscope (HDO6104A-MS, LeCroy, USA), high voltage differential probe (HVD3206A, LeCroy, USA) and current Hall probe (CP031A, LeCroy, Chestnut Ridge, NY, USA) were used.

Numerical Calculations
For each applied voltage at 4717 Hz repetition rate individual LSE-fit of the mathematical model, gave a precise fit of the time evolution of electric current (Figure 3a), (Table 1). However the aim was to develop a prototype of electroporation load emulator, therefore the same circuit should be used for all applied voltages. Therefore a fit of a circuit model was done on all three voltages together, based on which a slightly higher deviation from the average measurement was obtained. For both 1000 V and 750 V voltage amplitude, the modeled current is however still within the standard deviation of measurements (Figure 3b), while for 500 V the optimal fit better matches the 1 Hz repetition rate measurement. Table 1. Optimal values of electronic components for the single positive monopolar pulse model, for each applied voltage, at assumption U th = 200 V (Figure 3a). Used as initial values for LSE regression (Figure 3b).

voltage [V]
R The first modeled pulse in a sequence with higher repetition frequency, has the same time course as in a single monopolar pulse model, while all consecutive pulses have higher initial values (Figure 4a, first row). Matching between measurement and modeled current, slightly decreases with pulse number, but it is still within the standard deviation (Figure 3, blue line) of measurements for 750 V and 1000 V applied voltage amplitude. Optimal values of electronic components for a single-first pulse circuit emulator and for a sequence of pulses with higher repetition frequency are presented in table (Table 2).

PSpice Model and Electronic Emulator
Based on the PSpice model, a prototype emulator was built ( Figure 5). Available component values (Table 3) were implemented into PSpice transient analysis. The analysis was run for 1.6 ms. Coil internal resistance was set to measured winded wire resistance 15 Ω. The input voltage, rise time was equal to the average rise time of all applied voltage pulses,i.e., 0.358 µs and the same was done also for fall time, which was considerably slower 2.44 µs. The pulse width was set to 99.5 µs and period to 212 µs in order to match the measurements. The PSpice simulation results of the prototype device are presented in the second row of Figure 4b. Zener diode was composed of two Zener diodes in series, one with breakdown voltage 200 V (1N5388B, MULTICOMP) and another with 160 V (1N5384B, MULTICOMP). In order to implement the most suitable inductor, that would optimally operate at the applied pulse frequency, the inductor was made especially for this application. An RM14-type ferrite core made from N41 material grade with 1.9 mm gap (EPCOS/N41/RM14, 160 nH, 42.2 × 34.8 × 30.2mm), was winded in a 0.25 mm thick copper wire (number of windings: 600). The inductor was analyzed with an LCR meter (E4980A, Keysight, Santa Rosa, CA, USA), at 5 kHz the serial inductance was 57.6 mH and serial resistance 15 Ω and at 10 kHz the serial inductance was 57.8 mH and serial resistance 17.6 Ω (The inductor maximal constant current is 150 mA). CREE Silicone Carbide diode C2D05120 was used as a memory diode. Final value list of components built-in a prototype, is given in the table (Table 3). Measurements of the prototype emulator are presented in the following figure (Figure 4c, third row). An additional current measurement was added at the applied voltage set to 200 V, in order to present a current response, when the applied voltage is lower than the electroporation threshold voltage. Figure 4. Average of five current measurements at 4717 Hz repetition frequency (green), (a) numerically modeled current by Matlab (red, first row), (b) modeled current by a PSpice (blue, second row) and (c) measured current trough prototype emulator (purple, third row). An additional current measurement in the case of the emulator was added at the applied voltage set to 200 V, in order to present a current response, when the applied voltage is lower than the threshold voltage. In the first column, the current at 500 V in the second 750 V and third 1000 V applied voltage amplitude is presented.

Discussion
Unfortunately, in vivo measurements are mainly collected by means of the ECT Cliniporator device, of which sampling frequency is too low to accurately measure current spikes. Additionally, due to muscle contractions (which are triggered with pulse delivery), some noise can be present in in vivo measurements. Therefore this study is based on measurements of already published research by Langus et al. [19]. In theory measurements of the first pulse, for both repetition frequencies, displayed in Figure 3 (green and yellow) should be the same. However, this reflects a reality of ex vivo and in vivo measurement and inherent variations (due to the tissue inhomogeneity, error in electrode placement, difference in temperature and water content of the tissue). The method should be robust and insensitive to such variations which we demonstrate in our method by using imperfect measurements. The proposed numerical model thus slightly deviates from the measured values obtained in ex vivo experiment. The best fit was reached for the current measurement at the input voltage 750 V. This may be explained by the fact that the electric field was due to needle electrode use, slightly heterogeneous, meaning that cells that are closer to the needles, electroporate at lower voltages. In the case of plate electrodes, we believe the fit would be even better, because of a more homogeneous electric field distribution. It is important to always consider the fact, that different electrode geometries will "define" different load. Also the negative spike at the trailing edge of the current waveform was not modeled. The biggest challenge in the post pulse calculations is the impedance of the pulse generator. After the pulse generation, the impedance is changed (due to the turn-off, of the switch in electroporator) and cannot be defined from the current and voltage measurements, due to high noise and low voltage.
PSpice model and the current measurement of the prototype emulator match well. A Zener diode that defines the electroporation threshold voltage behaves as expected and therefore enables simulation of electric current before and during electroporation. The proposed mathematical, numerical and circuit models seems thus appropriate for the development of an electronic emulator of biological tissue as an electrical load during electroporation. The performance summary table (Table 4) presents the root mean square error values, between the measurements on beef liver and mathematical model or prototype emulator or resistor and mean standard deviation of five measurements on beef liver. In case of 1 Hz repetition rate RMSE was calculated only for one pulse, as all eight pulses have the same dynamics. While in case of 4717 Hz RMSE was calculated for the whole sequence of eight pulses. For the comparison of operation, the optimal resistance value from E24 resistance value list was selected. The minimal square error criteria was used and the value of 560 Ω gave the best match. Additionally, to the three main features (capacitive spike, rise of conductivity during the pulse and electroporation memory) the developed emulator in comparison to the optimal resistor has significantly lower RMSE, which can be observed in (Table 4). The presented circuit model could be upgraded for electroporation pulses of both polarities, named also biphasic or bipolar pulses. The last branch in the circuit should be doubled and mirrored as shown in (Scheme 4). Current measurements of both polarities should be however analyzed in detail in order to determine the relation between both inductors. If there is a memory present also between a positive, and negative pulse inductors L 1 and L 2 should be coupled. The solution, presented in this paper, presents a new approach, which enables repeatable and unbiased testing and evaluation of electroporator performance. However, the proposed emulator has to be adjusted for specific biological load i.e., cell-suspension in cuvettes, plate electrodes, electrodes for a skin treatment or needle electrodes for deep-seated tumors. Therefore, for each specific biological load, a separate analysis should be done. At first, we have to collect accurate measurements, with a high sampling rate, and than parametrization in accordance with this paper should be done. Finally, the circuit components should be modified accordingly to the newly calculated values. From the prospect of standardization or evaluation and comparison of different electroporation devices, each of the components in the proposed circuit could also be replaced with an adjustable value component. A series of a different capacitor and Zener diode modulation would incorporate flexibility and provide versatile emulator. The maximal voltage and current of an emulator are defined with the maximal voltage and current of the selected components and the frequency is limited by the selected inductor. Characteristics of the ferrite core are frequency dependent and in case of saturation the inductor would not work as desired.
The study in this paper is focused only on "classical" electroporation pulses, which are used clinically for ECT and IRE. Recently, however a much broader range of pulse duration is investigated and therefore the frequency spectrum of induced transmembrane potential is getting more attention [34]. With lower frequencies, the proposed concept can be used, but in case of higher frequencies, actual physical components should be thoughtfully selected with minimal parasitic values, or if not possible, the parasitic contributions should be added to the PSpice model. Also, a switch from circuital to distributed models of pulse delivery may be needed.

Conclusions
The proposed concept and electronic emulator of biological tissue as an electrical load during electroporation accurately simulates the current during electroporation. It serves as a proof of concept with a great potential. For each specific biological load, first voltage and current measurements should be made and a specific model should be developed. But from the prospect of standardization or evaluation and comparison of different electroporation devices, each of the components in the proposed circuit could be replaced with an adjustable value component. A series of a different capacitor and Zener diode modulation would incorporate flexibility and provide a versatile emulator. One of the possible solution is also that for each electroporation application a typical or worst case standardized emulator is developed. However additional measurements on different biological loads and different electrodes should be done and analyzed for this purpose. The main innovation of the presented process/concept of an emulator is that enables sustainable, repeatable and unbiased testing, therefore, it facilitates standard implementation. Additionally, the emulator can be used as a load during the development of an electroporator, which greatly facilitates the development process. We believe this is the first step towards the development of an electroporator tester. Such testers are well known and indispensable in other fields e.g., the defibrillation [35]. In the case of a defibrillator, a load is well known and electrodes are well defined, while in case of electroporation the same device can be used with different electrodes and load can vary from the skin to liver.