Sergej Kakorin

Fundamentals of electroporative delivery of drugs and genes

Electrooptical and conductometrical relaxation methods have given a new insight in the molecular mechanisms of the electroporative delivery of drug-like dyes and genes (DNA) to cells and tissues. Key findings are: (1) Membrane electroporation (ME) and hence the electroporative transmembrane transport of macromolecules are facilitated by a higher curvature of the membrane as well as by a gradient of the ionic strength across charged membranes, affecting the spontaneous curvature. (2) The degree of pore formation as the primary field response increases continuously without a threshold field strength, whereas secondary phenomena, such as a dramatic increase in the membrane permeability to drug-like dyes and DNA (also called electropermeabilization), indicate threshold field strength ranges. (3) The transfer of DNA by ME requires surface adsorption and surface insertion of the permeant molecule or part of it. The diffusion coefficient for the translocation of DNA (M(r) approximately 3.5 x 10(6)) through the electroporated membrane is Dm = 6.7 x 10(-13) cm2 s-1 and Dm for the drug-like dye Serva Blue G (M(r) approximately 854) is Dm = 2.0 x 10(-12) cm2 s-1. The slow electroporative transport of both DNA and drugs across the electroporated membrane reflects highly interactive (electro-) diffusion, involving many small pores coalesced into large, but transiently occluded pores (DNA). The data on mouse B-cells and yeast cells provide directly the flow and permeability coefficients of Serva blue G and plasmid DNA at different electroporation protocols. The physico-chemical theory of ME and electroporative transport in terms of time-dependent flow coefficients has been developed to such a degree that analytical expressions are available to handle curvature and ionic strength effects on ME and transport. The theory presents further useful tools for the optimization of the ME techniques in biotechnology and medicine, in particular in the new field of electroporative delivery of drugs (electrochemotherapy) and of DNA transfer and gene therapy.

Kinetics, Statistics, and Energetics of Lipid Membrane Electroporation Studied by Molecular Dynamics Simulations

Membrane electroporation is the method to directly transfer bioactive substances such as drugs and genes into living cells, as well as preceding electrofusion. Although much information on the microscopic mechanism has been obtained both from experiment and simulation, the existence and nature of possible intermediates is still unclear. To elucidate intermediates of electropore formation by direct comparison with measured prepore formation kinetics, we have carried out 49 atomistic electroporation simulations on a palmitoyl-oleoyl-phosphatidylcholine bilayer for electric field strengths between 0.04 and 0.7 V/nm. A statistical theory is developed to facilitate direct comparison of experimental (macroscopic) prepore formation kinetics with the (single event) preporation times derived from the simulations, which also allows us to extract an effective number of lipids involved in each pore formation event. A linear dependency of the activation energy for prepore formation on the applied field is seen, with quantitative agreement between experiment and simulation. The distribution of preporation times suggests a four-state pore formation model. The model involves a first intermediate characterized by a differential tilt of the polar lipid headgroups on both leaflets, and a second intermediate (prepore), where a polar chain across the bilayer is formed by 3-4 lipid headgroups and several water molecules, thereby providing a microscopic explanation for the polarizable volume derived previously from the measured kinetics. An average pore radius of 0.47 +/- 0.15 nm is seen, in favorable agreement with conductance measurements and electrooptical experiments of lipid vesicles.

Mechanism of Electroporative Dye Uptake by Mouse B Cells

The color change of electroporated intact immunoglobulin G receptor (Fc gammaR-) mouse B cells (line IIA1.6) after direct electroporative transfer of the dye SERVA blue G (Mr 854) into the cell interior is shown to be dominantly due to diffusion of the dye after the electric field pulse. Hence the dye transport is described by Ficks first law, where, as a novelty, time-integrated flow coefficients are introduced. The chemical-kinetic analysis uses three different pore states (P) in the reaction cascade (C <==> P1 <==> P2 <==> P3), to model the sigmoid kinetics of pore formation as well as the biphasic pore resealing. The rate coefficient for pore formation k(p) is dependent on the external electric field strength E and pulse duration tE. At E = 2.1 kV cm(-1) and tE = 200 micros, k(p) = (2.4 +/- 0.2) x 10(3) s(-1) at T = 293 K; the respective (field-dependent) flow coefficient and permeability coefficient are k(f)0 = (1.0 +/- 0.1) x 10(-2) s(-1) and P0 = 2 cm s(-1), respectively. The maximum value of the fractional surface area of the dye-conductive pores is 0.035 +/- 0.003%, and the maximum pore number is Np = (1.5 +/- 0.1) x 10(5) per average cell. The diffusion coefficient for SERVA blue G, D = 10(-6) cm2 s(-1), is slightly smaller than that of free dye diffusion, indicating transient interaction of the dye with the pore lipids during translocation. The mean radii of the three pore states are r(P1) = 0.7 +/- 0.1 nm, r(P2) = 1.0 +/- 0.1 nm, and r(P3) = 1.2 +/- 0.1 nm, respectively. The resealing rate coefficients are k(-2) = (4.0 +/- 0.5) x 10(-2) s(-1) and k(-3) = (4.5 +/- 0.5) x 10)(-3) s(-1), independent of E. At zero field, the equilibrium constant of the pore states (P) relative to closed membrane states (C) is K(p)0 = [(P)]/[C] = 0.02 +/- 0.002, indicating 2.0 +/- 0.2% water associated with the lipid membrane. Finally, the results of SERVA blue G cell coloring and the new analytical framework may also serve as a guideline for the optimization of the electroporative delivery of drugs that are similar in structure to SERVA blue G, for instance, bleomycin, which has been used successfully in the new discipline of electrochemotherapy.

Thermodynamics of the dimer-decamer transition of reduced human and plant 2-cys peroxiredoxin.

Isothermal titration calorimetry (ITC) is a powerful technique for investigating self-association processes of protein complexes and was expected to reveal quantitative data on peroxiredoxin oligomerization by directly measuring the thermodynamic parameters of dimer-dimer interaction. Recombinant classical 2-cysteine peroxoredoxins from Homo sapiens, Arabidopsis thaliana, and Pisum sativum as well as a carboxy-terminally truncated variant were subjected to ITC analysis by stepwise injection into the reaction vessel under various redox conditions. The direct measurement of the decamer-dimer equilibrium of reduced peroxiredoxin revealed a critical concentration in the very low micromolar range. The data suggest a cooperative assembly above this critical transition concentration where a nucleus facilitates assembly. The rather abrupt transition indicates that assembly processes do not occur below the critical transition concentration while oligomerization is efficiently triggered above it. The magnitude of the measured enthalpy confirmed the endothermic nature of the peroxiredoxin oligomerization. Heterocomplexes between peroxiredoxin polypeptides from different species were not formed. We conclude that a functional constraint conserved the dimer-decamer transition with highly similar critical transition concentrations despite emerging sequence variation during evolution.

Mechanism for the conductivity changes caused by membrane electroporation of CHO cell-pellets

Electric field pulses, applied to densely packed pellets of Chinese hamster ovary (CHO) cells of mean radius āc = 7.5 ± 0.7 μm, cause major electric conductivity changes, described by three kinetic normal modes. The first mode reflects Wien effects of ionic atmosphere perturbations and ion pair dissociations (cell surfaces). Using Maxwell’s conductivity equation, the second and third mode are converted to the respective membrane conductivity modes. Electrothermodynamic analysis in terms of structural transitions from closed (C) to porated (P) membrane states of very different lifetimes, according to the scheme (C ⇌ C1) ⇌ (P2 ⇌ P3), yields the mean pore radii 2 = 1.00 ± 0.05 nm (P2) pores and 3 = 1.5 ± 0.1 nm (P3) at T = 293 K (20 °C). The relaxation time τ2 (P2-formation) reflects the rate limiting step (C ⇌ C1), associated with the activation dipole moment of Δm1 = 63 × 10−30 C m (or 19 Debye units), suggesting orientational changes of dipolar lipid head groups, in the solution membrane interfaces preceding the actual pore formations. Besides: the field-dependencies of the pore fractions f2 and f3 (order of 10−3), the field reduction factors fλ,i≤ 1 and the membrane voltage, we obtain the zero-field pore conductivities λ0p,2 = 1.7 × 10−2 mS cm−1 (P2) and λ0p,3 = 0.10 mS cm−1 (P3) and the membrane conductivity λm0 = 3.2 μS m−1. The post-field conductivity changes, due to the long-lived P3-pores, are analyzed in terms of time- (and field-) dependent efflux coefficients. The characteristic post-field pore resealing time τR = τ03 = 45 ± 3 s is independent of the field strength of the causative pulse and independent of the distance between the two electrodes. These results are an essential part for the optimization of the electrical pulse parameters, also for the clinical electrotransfer of bioactive substances into aggregated biological cells (tissue).

Electro-optics of membrane electroporation in diphenylhexatriene-doped lipid bilayer vesicles

The electric (linear) dichroisms observed in the membrane electroporation of salt-filled lipid bilayer vesicles (diameter O = 2 alpha = 0.32 micron; inside [NaCl] = 0.2 M) in isotonic aqueous 0.284 M sucrose-0.2 mM NaCl solution indicate orientation changes of the anisotropic light scattering centers (lipid head groups) and of the optical transition moments of the membrane-inserted probe 1,6-diphenyl-1,3,5-hexatriene (DPH). Both the turbidity dichroism and DPH absorbance dichroism show peculiar features: (1) at external electric fields E > or = Esat the time course of the dichroism shows a maximum value (reversal): Esat = 4.0 (+/- 0.2) MV m-1, T = 293 K (20 degrees C), (2) this reversal value is independent of the field strength for E > or = Esat, (3) the dichroism amplitudes exhibit a maximum value Emax = 3.0 (+/- 0.5) MV m-1, (4) for the pulse duration of 10 microseconds there is one dominant visible normal mode, the relaxation rate increases up to tau-1 approximately 0.6 x 10(6) s-1 at Esat and then decreases for E > Esat. The data can be described in terms of local lipid phase transitions involving clusters Ln of n lipids in the pore edges according to the three-state scheme C<-->HO<-->HI, C being the closed bilayer state, HO the hydrophobic pore state and HI the hydrophilic or inverted pore state with rotated lipid and DPH molecules. At E > or = Esat, further transitions HO<-->HO* and HI<-->HI* are rapidly coupled to the C<-->HO transition, which is rate-limiting. The vesicle geometry conditions a cos theta dependence of the local membrane field effects relative to the E direction and the data reflect cos theta averages. The stationary induced transmembrane voltage delta phi (theta, lambda m) = -1.5 aEf(lambda m) magnitude of cos theta does not exceed the limiting value delta phi sat = -0.53 V, corresponding to the field strength Em,sat = -delta phi sat/d = 100 MV m-1 (10(3) kV cm-1), due to increasing membrane conductivity lambda m. At E = Esat, f(lambda m) = 0.55, lambda m = 0.11 mS m-1. The lipid cluster phase transition model yields an average pore radius of rp = 0.35 (+/- 0.05) nm of the assumed cylindrical pore of thickness d = 5 nm, suggesting an average cluster size of = 12 (+/- 2) lipids per pore edge. For E > Esat, the total number of DPH molecules in pore states approaches a saturation value; the fraction of DPH molecules in HI pores is 12 (+/- 2)% and that in HO pores is 48 (+/- 2)%. The percentage of membrane area P approximately (lambda m/lambda i) x 100% of conductive openings filled with the intravesicular medium of conductance lambda i = 2.2 S m-1 linearly increases from P approximately 0% (E = 1.8 MV m-1) to P = 0.017% (E = 8.5 MV m-1). Analogous estimations made by Kinosita et al. (1993) on the basis of fluorescence imaging data for sea urchin eggs give the same order of magnitude for P (0.02-0.2%). The increase in P with the field strength is collinear with the increase in concentration of HI and HI* states with the field strength, whereas the HO and HO* states exhibit a sigmoid field dependence. Therefore our data suggest that it is only the HI and HI* pore states which are conductive. It is noted that the various peculiar features of the dichroism data cannot be described by simple whole particle deformation.

Electroporative deformation of salt filled lipid vesicles

Abstract Membrane electroporation, vesicle shape deformation and aggregation of small, NaCl-filled lipid vesicles (of radius a = 50 nm) in DC electric fields was characterized using conductometric and turbidimetrical data. At pulse durations tE≤ 55 ± 5 ms the increase in the conductivity of the vesicle suspension is due to the field-induced efflux of electrolyte through membrane electropores. Membrane electroporation and Maxwell stress on the vesicle membrane lead to vesicle elongation concomitant with small volume reduction (up to 0.6% in an electric field of E = 1 MV m–1). At tE > 55 ± 5 ms, further increases in the conductivity and the optical density suggest electroaggregation and electrofusion of vesicles. The conductivity changes after the electric pulse termination reflect salt ion efflux through slowly resealing electropores. The analysis of the volume reduction kinetics yields the bending rigidity κ = (4.1 ± 0.3) ⋅ 10–20 J of the vesicle membrane. If the flow of Na+ and Cl– ions from the vesicle interior is treated in terms of Hagen-Poiseuilles equation, the number of permeable electropores is N = 39 per vesicle with mean pore radius rp = 0.85 ± 0.05 nm at E = 1 MVm–1 and tE≤ 55 ± 5 ms. The turbidimetric and conductometric data suggest that small lipid vesicles (a ≤ 50 nm) are not associated with extensive membrane thermal undulations or superstructures. In particular with respect to membrane curvature, the vesicle results are suggestive for the design and optimization of electroporative delivery of drugs and genes to cell tissue at small field strengths (≤1 MVm–1) and large pulse durations (≤100 ms).

Electrooptics of membrane electroporation and vesicle shape deformation

Electrooptical and conductometric methods continue to reveal new and more detailed information on the dynamic properties of membranes in electric fields. In particular, the electric pore formation in lipid vesicles, doped with optical probes, has been successfully investigated with electrooptical techniques thus providing new insight into the lipid rearrangements underlying membrane electroporation (ME) and vesicle deformation. Progress in understanding the molecular mechanism of ME and related phenomena, such as electrofusion of cells or electroinsertion of foreign proteins into membranes, is crucially important for the numerous applications of ME, for example, direct electroporative gene transfer and drug delivery in the new medical discipline of electroporative chemotherapy.

Membrane electroporation and electromechanical deformation of vesicles and cells

Analysis of the reduced turbidity (delta T-/T0) and absorbance (delta A-/A0) relaxations of unilamellar lipid vesicles, doped with the diphenylhexatrienyl-phosphatidylcholine (beta-DPH pPC) lipids in high-voltage rectangular electrical field pulses, demonstrates that the major part of the turbidity and absorbance dichroism is caused by vesicle elongation under electric Maxwell stress. The kinetics of this electrochemomechanical shape deformation (time constants 0.1 < or = tau/microsecond < or = 3) is determined both by the entrance of water and ions into the bulk membrane phase to form local electropores, and by the faster processes of membrane stretching and smoothing of thermal undulations. Moreover, the absorbance dichroism indicates local displacements of the chromophore relative to the membrane normal in the field. The slightly slower relaxations of the chemical turbidity (delta T+/T0) and absorbance (delta A+/A0) modes are both associated with the entrance of solvent into the interface membrane/medium, caused by the alignment of the bipolar lipid head groups in one of the leaflets at the pole caps of the vesicle bilayer. In addition, (delta T+/T0) indicates changes in vesicle shape and volume. The results for lipid vesicles provide guidelines for the analysis of electroporative deformations of biological cells.

Ionic conductivity of electroporated lipid bilayer membranes

The ionic conductivity of lipid membrane pores has been theoretically analysed in terms of electrostatic interactions of the transported ions with the low-dielectric pore wall for a commonly encountered case of unequal concentrations of electrolyte on the two sides of curved lipid membranes. Theoretical analysis of the data on the conductivity of the electroporated membrane of lipid vesicles (Lecithin 20%) of radius a=90 nm yields the molar energy of interaction of a small monovalent ion with a pore wall w(0)=9+/-1 RT (or w(0)=22+/-kJ mol(-1)), corresponding to a mean pore radius of (-)r(p)=0.56+/-0.05 nm. The proposed theoretical approach provides a tool for the analysis and description of the nonlinear current-voltage dependencies in membrane pores and channels.