MATHEMATICAL MODELING OF MOLECULAR TRANSMEMBRANE TRANSPORT AND CHANGES OF TISSUES´ DIELECTRIC PROPERTIES DUE TO ELECTROPORATIONDERMOL-ČERNE, JANJA (Avtor) Miklavčič, Damijan (Mentor) Serša, Gregor (Komentor) Electroporationmultiscale modelingexcitable cellspredictive modelselectrochemotherapynumerical modelingirreversible electroporationElectroporation is a phenomenon, which occurs when short high voltage pulses are applied to cells and tissues resulting in a transient increase in membrane permeability or cell death, presumably due to pore formation. If cells recover after pulse application, this is reversible electroporation. If cells die, this is irreversible electroporation. Electroporation is used in biotechnology for biocompound extraction and cryopreservation, in food processing for sterilization and pasteurization of liquid food and in medicine for treating tumors by electrochemotherapy or irreversible electroporation as an ablation technique, for gene electrotransfer, transdermal drug delivery, DNA vaccination, and cell fusion. In electroporation-based medical treatments, we can treat tumors with predefined electrode geometry and parameters of electric pulses. When we treat larger tumors of irregular shape treatment plan of the position of the electrodes and parameters of the electric pulses has to be calculated before each treatment to assure coverage of the tumor with a sufficient electric field. In treatment plans, currently, 1) we assume that above an experimentally determined critical electric field all cells are affected and below not, although, in reality, the transition between non-electroporated and electroporated state is continuous. 2) We do not take into account the excitability of some tissues. 3) The increase in tissues’ conductivity is described phenomenologically and does not include mechanisms of electroporation. 4) Transport of chemotherapeutics into the tumor cells in electrochemotherapy treatments is not included in the treatment plan although it is vital for a successful treatment. We focused on the mathematical and numerical models of electroporation with the aim of including them in the treatment planning of electroporation-based medical treatments. We aimed to model processes happening during electroporation of tissues, relevant in the clinical procedures, by taking into account processes happening at the single cell level. First, we used mathematical models of cell membrane permeability and cell death which are phenomenological descriptions of experimental data. The models were chosen on the basis of the best fit with the experimental data. However, they did not include mechanisms of electroporation, and their transferability to tissues was questionable. We modeled time dynamics of dye uptake due to increased cell membrane permeability in several electroporation buffers with regard to the electrosensitization, i.e., delayed hypersensitivity to electric pulses caused by pretreating cells with electric pulses. We also modeled the strength-duration depolarization curve and cell membrane permeability curve of excitable and non-excitable cell lines which could be used to optimize pulse parameters to achieve maximal drug uptake at minimal tissue excitation. Second, we modeled change in dielectric properties of tissues during electroporation. Model of change in dielectric properties of tissues was built for skin and validated with current-voltage measurements. Dielectric properties of separate layers of skin before electroporation were determined by taking into account geometric and dielectric properties of single cells, i.e., keratinocytes, corneocytes. Dielectric properties of separate layers during electroporation were obtained from cell-level models of pore formation on single cells of lower skin layers (keratinocytes in epidermis and lipid spheres in papillary dermis) and local transport region formation in the stratum corneum. Current-voltage measurements of long low-voltage pulses were accurately described taking into account local transport region formation, pore formation in the cells of lower layers and electrode polarization. Voltage measurements of short high-voltage pulses were also accurately described in a similar way as with long low-voltage pulses; however, the model underestimated the current, probably due to electrochemical reactions taking place at the electrode-electrolyte interface. Third, we modeled the transport of chemotherapeutics during electrochemotherapy in vivo. In electrochemotherapy treatments, transport of chemotherapeutics in sufficient amounts into the cell is vital for a successful treatment. We performed experiments in vitro and measured the intracellular platinum mass as a function of pulse number and electric field by inductively coupled plasma – mass spectrometry. Using the dualporosity model, we calculated the in vitro permeability coefficient as a function of electric field and number of applied pulses. The in vitro determined permeability coefficient was then used in the numerical model of mouse melanoma tumor to describe the transport of cisplatin to the tumor cells. We took into account the differences in the transport of cisplatin in vitro and in vivo caused by the decreased mobility of molecules and decreased membrane area available for the uptake in vivo due to the high volume fraction of cells, the presence of cell matrix and close cell connections. Our model accurately described the experimental results obtained in electrochemotherapy of tumors and could be used to predict the efficiency of electrochemotherapy in vitro thus reducing the number of needed animal experiments. In the thesis, we connected the models at the cell level to the models at the tissue level with respect to cell membrane permeability and depolarization, cell death, change in dielectric properties and transport. Our models offer a step forward in modeling and understanding electroporation at the tissue level. In future, our models could be used to improve treatment planning of electroporation-based medical treatments.20182018-05-18 11:02:29Doktorsko delo/naloga101263VisID: 41881sl