Coupled mathematical modeling of cisplatin electroporation.

The use of electrochemotherapy (ECT) is a well-established technique to increase the cellular uptake of cytotoxic agents within certain cancer treatment strategies. The study of the mechanisms that take part in this complex process is of high interest to gain a deeper knowledge of it, enabling the improvement of these strategies. In this work, we present a coupled multi-physics electroporation model based on a related previous one, to describe the effect of a set of electric pulses on cisplatin transport across the plasma membrane. The model applies a system of partial differential equations that includes Poisson's equation for the electric field, Nernst-Planck's equation for species transport, Maxwell's tensor and mechanical equilibrium equation for membrane deformation and Smoluchowski's equation for pore creation dynamics. Our numerical results were compared with previous numerical and experimental published data with good qualitative and quantitative agreement. These results indicate that pore aperture is favored at the cell poles by the electric field and mechanical stress forces, giving support to the dominant hypothesis of hydrophilic pore creation as the main mechanism of drug entry during an ECT treatment.

Cell membrane electroporation modeling: A multiphysics approach.

Electroporation-based techniques, i.e. techniques based on the perturbation of the cell membrane through the application of electric pulses, are widely used at present in medicine and biotechnology. However, the electric pulse - cell membrane interaction is not yet completely understood neither explicitly formalized. Here we introduce a Multiphysics (MP) model describing electric pulse - cell membrane interaction consisting on the Poisson equation for the electric field, the Nernst-Planck equations for ion transport (protons, hydroxides, sodium or calcium, and chloride), the Maxwell tensor and mechanical equilibrium equation for membrane deformations (with an explicit discretization of the cell membrane), and the Smoluchowski equation for membrane permeabilization. The MP model predicts that during the application of an electric pulse to a spherical cell an elastic deformation of its membrane takes place affecting the induced transmembrane potential, the pore creation dynamics and the ionic transport. Moreover, the coincidence among maximum membrane deformation, maximum pore aperture, and maximum ion uptake is predicted. Such behavior has been corroborated experimentally by previously published results in red blood and CHO cells as well as in supramolecular lipid vesicles.

Microenvironmental influence on microtumour infiltration patterns: 3D-mathematical modelling supported by in vitro studies.

Mathematical modelling approaches have become increasingly abundant in cancer research. Tumour infiltration extent and its spatial organization depend both on the tumour type and stage and on the bio-physicochemical characteristics of the microenvironment. This sets a complex scenario that often requires a multidisciplinary and individually adjusted approach. The ultimate goal of this work is to present an experimental/numerical combined method for the development of a three-dimensional mathematical model with the ability to reproduce the growth and infiltration patterns of a given avascular microtumour in response to different microenvironmental conditions. The model is based on a diffusion-convection reaction equation that considers logistic proliferation, volumetric growth, a rim of proliferative cells at the tumour surface, and invasion with diffusive and convective components. The parameter values of the model were fitted to experimental results while radial velocity and diffusion coefficients were made spatially variable in a case-specific way through the introduction of a shape function and a diffusion-limited-aggregation (DLA)-derived fractal matrix, respectively, according to the infiltration pattern observed. The in vitro model consists of multicellular tumour spheroids (MTSs) of an epithelial mammary tumour cell line (LM3) immersed in a collagen I gel matrix with a standard culture medium ("naive" matrix) or a conditioned medium from adipocytes or preadipocytes ("conditioned" matrix). It was experimentally determined that both adipocyte and preadipocyte conditioned media had the ability to change the MTS infiltration pattern from collective and laminar to an individual and atomized one. Numerical simulations were able to adequately reproduce qualitatively and quantitatively both kinds of infiltration patterns, which were determined by area quantification, analysis of fractal dimensions and lacunarity, and Bland-Altman analysis. These results suggest that the combined approach presented here could be established as a new framework with interesting potential applications at both the basic and clinical levels in the oncology area.

Adaptive numerical algorithms to simulate the dynamical Casimir effect in a closed cavity with different boundary conditions.

We present an alternative numerical approach to compute the number of particles created inside a cavity due to time-dependent boundary conditions. The physical model consists of a rectangular cavity, where a wall always remains still while the other wall of the cavity presents a smooth movement in one direction. The method relies on the setting of the boundary conditions (Dirichlet and Neumann) and the following resolution of the corresponding equations of modes. By a further comparison between the ground state before and after the movement of the cavity wall, we finally compute the number of particles created. To demonstrate the method, we investigate the creation of particle production in vibrating cavities, confirming previously known results in the appropriate limits. Within this approach, the dynamical Casimir effect can be investigated, making it possible to study a variety of scenarios where no analytical results are known. Of special interest is, of course, the realistic case of the electromagnetic field in a three-dimensional cavity, with transverse electric (TE)-mode and transverse magnetic (TM)-mode photon production. Furthermore, with our approach we are able to calculate numerically the particle creation in a tuneable resonant superconducting cavity by the use of the generalized Robin boundary condition. We compare the numerical results with analytical predictions as well as a different numerical approach. Its extension to three dimensions is also straightforward.

Mathematical modelling of microtumour infiltration based on in vitro experiments.

The present mathematical models of microtumours consider, in general, volumetric growth and spherical tumour invasion shapes. Nevertheless in many cases, such as in gliomas, a need for more accurate delineation of tumour infiltration areas in a patient-specific manner has arisen. The objective of this study was to build a mathematical model able to describe in a case-specific way as well as to predict in a probabilistic way the growth and the real invasion pattern of multicellular tumour spheroids (in vitro model of an avascular microtumour) immersed in a collagen matrix. The two-dimensional theoretical model was represented by a reaction-convection-diffusion equation that considers logistic proliferation, volumetric growth, a rim with proliferative cells at the tumour surface and invasion with diffusive and convective components. Population parameter values of the model were extracted from the experimental dataset and a shape function that describes the invasion area was derived from each experimental case by image processing. New possible and aleatory shape functions were generated by data mining and Monte Carlo tools by means of a satellite EGARCH model, which were fed with all the shape functions of the dataset. Then the main model is used in two different ways: to reproduce the growth and invasion of a given experimental tumour in a case-specific manner when fed with the corresponding shape function (descriptive simulations) or to generate new possible tumour cases that respond to the general population pattern when fed with an aleatory-generated shape function (predictive simulations). Both types of simulations are in good agreement with empirical data, as it was revealed by area quantification and Bland-Altman analysis. This kind of experimental-numerical interaction has wide application potential in designing new strategies able to predict as much as possible the invasive behaviour of a tumour based on its particular characteristics and microenvironment.

The role of additional pulses in electropermeabilization protocols.

Electropermeabilization (EP) based protocols such as those applied in medicine, food processing or environmental management, are well established and widely used. The applied voltage, as well as tissue electric conductivity, are of utmost importance for assessing final electropermeabilized area and thus EP effectiveness. Experimental results from literature report that, under certain EP protocols, consecutive pulses increase tissue electric conductivity and even the permeabilization amount. Here we introduce a theoretical model that takes into account this effect in the application of an EP-based protocol, and its validation with experimental measurements. The theoretical model describes the electric field distribution by a nonlinear Laplace equation with a variable conductivity coefficient depending on the electric field, the temperature and the quantity of pulses, and the Penne's Bioheat equation for temperature variations. In the experiments, a vegetable tissue model (potato slice) is used for measuring electric currents and tissue electropermeabilized area in different EP protocols. Experimental measurements show that, during sequential pulses and keeping constant the applied voltage, the electric current density and the blackened (electropermeabilized) area increase. This behavior can only be attributed to a rise in the electric conductivity due to a higher number of pulses. Accordingly, we present a theoretical modeling of an EP protocol that predicts correctly the increment in the electric current density observed experimentally during the addition of pulses. The model also demonstrates that the electric current increase is due to a rise in the electric conductivity, in turn induced by temperature and pulse number, with no significant changes in the electric field distribution. The EP model introduced, based on a novel formulation of the electric conductivity, leads to a more realistic description of the EP phenomenon, hopefully providing more accurate predictions of treatment outcomes.

Nature of inclined growth in thin-layer electrodeposition under uniform magnetic fields.

Electrochemical deposition (ECD) in thin cells in a vertical position relative to gravity, subject to an external uniform magnetic field, yields a growth pattern formation with dense branched morphology with branches tilted in the direction of the magnetic force. We study the nature of the inclined growth through experiments and theory. Experiments in ECD, in the absence of magnetic forces, reveal that a branch grows by allowing fluid to penetrate its tip and to be ejected from the sides through a pair of symmetric vortices attached to the tip. The upper vortices zone defines an arch separating an inner zone ion depleted and an outer zone in a funnel-like form with a concentrated solution through which metal ions are carried into the tip. When a magnetic field is turned on, vortex symmetry is broken, one vortex becoming weaker than the other, inducing an inclination of the funnel. Consequently, particles entering the funnel give rise to branch growth tilted in the same direction. Theory predicts, in the absence of a magnetic force, funnel symmetry induced through symmetric vortices driven by electric and gravitational forces; when the magnetic force is on, it is composed with the pair of clockwise and counterclockwise vortices, reducing or amplifying one or the other. In turn, funnel tilting modifies particle trajectories, thus, growth orientation.

Dynamics of a paramagnetic colloidal particle driven on a magnetic-bubble lattice.

We present a theoretical study of the recently observed dynamical regimes of paramagnetic colloidal particles externally driven above a regular lattice of magnetic bubbles [P. Tierno, T. H. Johansen, and T. M. Fischer, Phys. Rev. Lett. 99, 038303 (2007)]. An external precessing magnetic field alters the potential generated by the surface of the film in such a way to either drive the particle circularly around one bubble, ballistically through the array, or in triangular orbits on the interstitial regions between the bubbles. In the ballistic regime, we observe different trajectories performed by the particles phase locked with the external driving. Superdiffusive motion, which was experimentally found bridging the localized and delocalized dynamics, emerge only by introducing a certain degree of randomness into the bubbles size distribution.