Phenomenological Electropermeabilization Model

• Phenomenological electropermeabilization models are quantitative frameworks that predict electric field-induced pore creation, membrane resealing, and drug transport.
• They integrate nonlinear kinetics, electric field distributions, and thermal effects to optimize pulse protocols for safe and effective treatments.
• The models use sigmoid and exponential laws validated via finite-difference simulations to enable reproducible, patient-specific clinical optimization.

Phenomenological electropermeabilization models provide quantitative, mathematically explicit frameworks to predict the dynamics of membrane permeability, drug delivery, and tissue response under externally applied electric fields. At their core, these models represent electroporation-induced pore creation, resealing, and the consequent transient enhancement of mass transfer across the cell membrane. By coupling electric-field distribution, pore-density kinetics, drug transport through extracellular and intracellular compartments, and, where relevant, thermal effects due to Joule heating, such models enable parameter optimization for clinical applications including electrochemotherapy and gene transfer. They are distinguished from molecular-scale or atomistic approaches by their reliance on macroscopic, empirically parameterized laws—often sigmoid or exponential in structure—governing permeability states and transfer coefficients.

1. Mathematical Framework and Governing Equations

The foundational structure is an integro-differential system coupling Laplace-type electric field equations, nonlinear membrane permeabilization kinetics, mass-transfer models, and often thermal conduction. For a cubical tissue block between parallel electrodes, the potential $\phi(x,y)$ obeys

$$\nabla^2\phi(x,y)=0, \quad \phi(0,y)=\phi_0,\;\phi(L,y)=\phi_L,$$

with electric field $E=-\nabla\phi$ and transmembrane potential on a spherical cell given by

$V_m(\Psi) = 1.5\,E\,r_c\cos\Psi.$

Membrane pore population $N(t)$ evolves via a Krassowska–Filev-type ODE

$$\frac{dN}{dt} = \alpha\,A(V_m)\left[1-\frac{N}{N_0}A(V_m)^{-q}\right], \quad A(V_m)=\exp\left(\frac{V_m^2}{V_{ep}^2}\right), \quad N(0)=0,$$

with explicit spatial dependence for non-homogeneous tissue (Mondal et al., 2022).

The post-pulse resealing phase leads to time-dependent pore area

$A_P(t) = \pi\,R_p^2\,N(t_p)\,\exp(-t/\tau),$

and the crucial mass-transfer coefficient (MTC) for drug permeation,

$$k(t)=\mu(t) = \frac{4\pi^2\,R_p^2\,r_c^2}{V_0} P N(t_p) \exp(-t/\tau),$$

serves as the rate constant for ECS $\leftrightarrow$ ICS transport in the system of coupled reaction–diffusion PDEs: $$\frac{\partial C_E}{\partial t} = \nabla\cdot(D\nabla C_E) - \frac{1-\varepsilon}{\varepsilon} k(t) (C_E - C_{RE}),$$

$$\frac{\partial C_{RE}}{\partial t} = k(t) (C_E - C_{RE}).$$

Thermal evolution during electroporation pulses is captured by

$$\rho\,c\,\frac{\partial T}{\partial t} = k_T\,\nabla^2T + \sigma\,|\nabla\phi|^2,$$

with Joule heating suppressed during OFF intervals (Mondal et al., 2022).

2. Pore-Density Law, Sigmoid Thresholds, and Membrane State Variables

An essential phenomenological element is the empirically derived relationship between local electric field and pore formation/closure. Mondal & Dalal (Mondal et al., 2022) advocate a sigmoid pore-fraction law,

$$\varphi(E) = \left[1+\exp\left(\frac{E_f-E}{b_f}\right)\right]^{-1},$$

where $E_f$ is the field midpoint and $b_f$ the steepness parameter. The MTC is parametrized as

$k_m(t) = \frac{P\varphi}{r_c}\exp(-t/\tau),$

and pulse duration $t_p$ as well as resealing time $\tau$ critically affect uptake kinetics (Mondal et al., 2022).

State variables tracking membrane permeabilization ($X_1$, $X_2$) in multicell/aggregate models evolve under voltage-dependent ODEs of the form

$$\partial_t X_1 = [\exp(-V_{ep}^2/\Delta V_m^2) - X_1]/\tau_{ep},$$

$$\partial_t X_2 = \max\left\{[\exp(-X_{ep}^2/X_1^2) - X_2]/\tau_{perm}, [\cdots]/\tau_{res}\right\},$$

with stepwise transitions to high-conductance states as the field exceeds characteristic voltage thresholds (Mistani et al., 2018).

3. Parameter Dependence, Pulse Protocols, and Drug Uptake Dynamics

The quantitative performance of the model is governed by pulse number, pulse duration, inter-pulse interval, field strength, membrane permeability, extracellular space porosity, and drug diffusion coefficient. Key findings include:

• Increasing $P$ (membrane permeability) or $D$ (extracellular diffusion) accelerates intracellular drug loading, with diminishing returns if the source depletes.
• For long OFF intervals ($t_M$) and repeated low-voltage pulses ($E \le 28$ V/mm, $N \le 10$), drug uptake is maximized without risking thermal damage.
• MTC rises with $P$ and $t_p$; decays exponentially with $\tau$.
• An optimal protocol consists of $E \approx 20$–$28$ V/mm, $t_p \approx 50$–$80$ ms, $N \approx 5$–$10$, $t_M \approx 600$ s, keeping the temperature increase below $2$ K (Mondal et al., 2022).
• Multiple short pulses overcome the exponential decay of MTC and establish near-equilibrium between ECS and ICS after saturating the former (Mondal et al., 2022).

4. Numerical Solution Strategies and Boundary Effects

The models utilize finite-difference (FTCS) approaches on Cartesian grids (e.g., $101 \times 101$ nodes), with stability bounds enforced on the time step. Drug boundary loss is modeled via Neumann-type conditions,

$\frac{\partial C_E}{\partial n} = \beta C_E,$

where $\beta$ quantifies ECS leakage, demonstrating that higher boundary losses confine penetration to a thin layer near the injection face (Mondal et al., 2022). Explicit formulas for FTCS update coefficients as a function of $D$, $\varepsilon$, and grid parameters are stated, supporting reproducibility in simulation (Mondal et al., 2022).

5. Thermodynamic Safety and Thermal Modeling

Phenomenological models integrate thermal safety by evaluating the maximum tissue temperature during a pulse train. For typical clinical parameters (ten $80$ ms pulses at $28$ V/mm), simulated $T_{\max}\approx 314.3$ K never exceeds $315.15$ K ($42^\circ$C), thus avoiding heat-induced damage within reversible electroporation protocols (Mondal et al., 2022). Coupling mass transfer calculations with thermal profiles allows rigorous definition of boundaries between reversible, safe operation and damaging irreversible effects.

6. Parameter Optimization and Clinical Applicability

These models constitute tools for designing patient-specific dose regimens, enabling in silico optimization of pulse parameters prior to costly experimental validation. By tuning $E$, $P$, $D$, and tissue geometry, the framework can be adapted for personalized electrochemotherapy or gene transfer, predicting intracellular drug concentrations as functions of controllable variables (Mondal et al., 2022, Mondal et al., 2022). Extensions include relaxing spatial homogeneity, incorporating cell-size distributions, and coupling with irreversible electroporation models for tumor ablation scenarios.

7. Comparison to Molecular and Energy-Landscape Models

While sharing pore-kinetic elements with Smoluchowski or energy-barrier-based models (Gowrishankar et al., 2017, Yu et al., 2014), the phenomenological framework for tissue-level reversible electroporation omits atomistic pore size distributions and instead aggregates their effects into macroscopic coefficients or sigmoid state variables. This enables fast, scalable computation and direct parameter linkage to clinical pulse generators, tissue electrical properties, and observable uptake metrics.

---

The phenomenological electropermeabilization model systematically links electric field–induced pore formation, membrane resealing, drug mass transfer, and tissue thermal response. Explicit analytical and numerical treatment of each process enables comprehensive clinical protocol optimization and prediction of outcomes in reversible tissue electroporation-based drug delivery (Mondal et al., 2022, Mondal et al., 2022).