A Model for Membrane Degradation Using a Gelatin Invadopodia Assay.

One of the most crucial and lethal characteristics of solid tumors is represented by the increased ability of cancer cells to migrate and invade other organs during the so-called metastatic spread. This is allowed thanks to the production of matrix metalloproteinases (MMPs), enzymes capable of degrading a type of collagen abundant in the basal membrane separating the epithelial tissue from the connective one. In this work, we employ a synergistic experimental and mathematical modelling approach to explore the invasion process of tumor cells. A mathematical model composed of reaction-diffusion equations describing the evolution of the tumor cells density on a gelatin substrate, MMPs enzymes concentration and the degradation of the gelatin is proposed. This is completed with a calibration strategy. We perform a sensitivity analysis and explore a parameter estimation technique both on synthetic and experimental data in order to find the optimal parameters that describe the in vitro experiments. A comparison between numerical and experimental solutions ends the work.

Modeling ATP-mediated endothelial cell elongation on line patterns.

Endothelial cell (EC) migration is crucial for a wide range of processes including vascular wound healing, tumor angiogenesis, and the development of viable endovascular implants. We have previously demonstrated that ECs cultured on 15-µm wide adhesive line patterns exhibit three distinct migration phenotypes: (a) "running" cells that are polarized and migrate continuously and persistently on the adhesive lines with possible spontaneous directional changes, (b) "undecided" cells that are highly elongated and exhibit periodic changes in the direction of their polarization while maintaining minimal net migration, and (c) "tumbling-like" cells that migrate persistently for a certain amount of time but then stop and round up for a few hours before spreading again and resuming migration. Importantly, the three migration patterns are associated with distinct profiles of cell length. Because of the impact of adenosine triphosphate (ATP) on cytoskeletal organization and cell polarization, we hypothesize that the observed differences in EC length among the three different migration phenotypes are driven by differences in intracellular ATP levels. In the present work, we develop a mathematical model that incorporates the interactions between cell length, cytoskeletal (F-actin) organization, and intracellular ATP concentration. An optimization procedure is used to obtain the model parameter values that best fit the experimental data on EC lengths. The results indicate that a minimalist model based on differences in intracellular ATP levels is capable of capturing the different cell length profiles observed experimentally.

Numerical optimization of plasmid DNA delivery combined with hyaluronidase injection for electroporation protocol.

BACKGROUND AND OBJECTIVE: The paper focuses on the numerical strategies to optimize a plasmid DNA delivery protocol, which combines hyaluronidase and electroporation. METHODS: A well-defined continuum mechanics model of muscle porosity and advanced numerical optimization strategies have been used, to propose a substantial improvement of a pre-existing experimental protocol of DNA transfer in mice. Our work suggests that a computational model might help in the definition of innovative therapeutic procedures, thanks to the fine tuning of all the involved experimental steps. This approach is particularly interesting in optimizing complex and costly protocols, to make in vivo DNA therapeutic protocols more effective. RESULTS: Our preliminary work suggests that computational model might help in the definition of innovative therapeutic protocol, thanks to the fine tuning of all the involved operations. CONCLUSIONS: This approach is particularly interesting in optimizing complex and costly protocols for which the number of degrees of freedom prevents a experimental test of the possible configuration.

A Continuum Mechanics Model of Enzyme-Based Tissue Degradation in Cancer Therapies.

We propose a mathematical model to describe enzyme-based tissue degradation in cancer therapies. The proposed model combines the poroelastic theory of mixtures with the transport of enzymes or drugs in the extracellular space. The effect of the matrix-degrading enzymes on the tissue composition and its mechanical response are accounted for. Numerical simulations in 1D, 2D and axisymmetric (3D) configurations show how an injection of matrix-degrading enzymes alters the porosity of a biological tissue. We eventually exhibit numerically the main consequences of a matrix-degrading enzyme pretreatment in the framework of chemotherapy: the removal of the diffusive hindrance to the penetration of therapeutic molecules in tumors and the reduction of interstitial fluid pressure which improves transcapillary transport. Both effects are consistent with previous biological observations.

A discrete in continuous mathematical model of cardiac progenitor cells formation and growth as spheroid clusters (Cardiospheres).

We propose a discrete in continuous mathematical model describing the in vitro growth process of biophsy-derived mammalian cardiac progenitor cells growing as clusters in the form of spheres (Cardiospheres). The approach is hybrid: discrete at cellular scale and continuous at molecular level. In the present model, cells are subject to the self-organizing collective dynamics mechanism and, additionally, they can proliferate and differentiate, also depending on stochastic processes. The two latter processes are triggered and regulated by chemical signals present in the environment. Numerical simulations show the structure and the development of the clustered progenitors and are in a good agreement with the results obtained from in vitro experiments.

A Macroscopic Mathematical Model for Cell Migration Assays Using a Real-Time Cell Analysis.

Experiments of cell migration and chemotaxis assays have been classically performed in the so-called Boyden Chambers. A recent technology, xCELLigence Real Time Cell Analysis, is now allowing to monitor the cell migration in real time. This technology measures impedance changes caused by the gradual increase of electrode surface occupation by cells during the course of time and provide a Cell Index which is proportional to cellular morphology, spreading, ruffling and adhesion quality as well as cell number. In this paper we propose a macroscopic mathematical model, based on advection-reaction-diffusion partial differential equations, describing the cell migration assay using the real-time technology. We carried out numerical simulations to compare simulated model dynamics with data of observed biological experiments on three different cell lines and in two experimental settings: absence of chemotactic signals (basal migration) and presence of a chemoattractant. Overall we conclude that our minimal mathematical model is able to describe the phenomenon in the real time scale and numerical results show a good agreement with the experimental evidences.

The dynamics of p53 in single cells: physiologically based ODE and reaction-diffusion PDE models.

The intracellular signalling network of the p53 protein plays important roles in genome protection and the control of cell cycle phase transitions. Recently observed oscillatory behaviour in single cells under stress conditions has inspired several research groups to simulate and study the dynamics of the protein with the aim of gaining a proper understanding of the physiological meanings of the oscillations. We propose compartmental ODE and PDE models of p53 activation and regulation in single cells following DNA damage and we show that the p53 oscillations can be retrieved by plainly involving p53-Mdm2 and ATM-p53-Wip1 negative feedbacks, which are sufficient for oscillations experimentally, with no further need to introduce any delays into the protein responses and without considering additional positive feedback.

Gene therapy: the role of cytoskeleton in gene transfer studies based on biology and mathematics.

Gene therapy is a promising approach for treating a wide range of human pathologies such as genetic disorders as well as diseases acquired over time. Viral and non-viral vectors are used to convey sequences of genes that can be expressed for therapeutic purposes. Plasmid DNA is receiving considerable attention for intramuscular gene transfer due to its safety, simplicity and low cost of production. Nevertheless, strategies to improve DNA uptake into the nucleus of cells for its expression are required. Cytoskeleton plays an important role in the intracellular trafficking. The mechanism regulating this process must be elucidated. Here, we propose a new methodological approach based on the coupling of biology assays and predictive mathematical models, in order to clarify the mechanism of the DNA uptake and its expression into the cells. Once these processes are better clarified, we will be able to propose more efficient therapeutic gene transfer protocols for the treatment of human patients.

The p53 protein and its molecular network: modelling a missing link between DNA damage and cell fate.

Various molecular pharmacokinetic-pharmacodynamic (PK-PD) models have been proposed in the last decades to represent and predict drug effects in anticancer chemotherapies. Most of these models are cell population based since clearly measurable effects of drugs can be seen much more easily on populations of cells, healthy and tumour, than in individual cells. The actual targets of drugs are, however, cells themselves. The drugs in use either disrupt genome integrity by causing DNA strand breaks, and consequently initiate programmed cell death, or block cell proliferation mainly by inhibiting factors that enable cells to proceed from one cell cycle phase to the next through checkpoints in the cell division cycle. DNA damage caused by cytotoxic drugs (and also cytostatic drugs at high concentrations) activates, among others, the p53 protein-modulated signalling pathways that directly or indirectly force the cell to make a decision between survival and death. The paper aims to become the first-step in a larger scale enterprise that should bridge the gap between intracellular and population PK-PD models, providing oncologists with a rationale to predict and optimise the effects of anticancer drugs in the clinic. So far, it only sticks at describing p53 activation and regulation in single cells following their exposure to DNA damaging stress agents. We show that p53 oscillations that have been observed in individual cells can be reconstructed and predicted by compartmentalising cellular events occurring after DNA damage, either in the nucleus or in the cytoplasm, and by describing network interactions, using ordinary differential equations (ODEs), between the ATM, p53, Mdm2 and Wip1 proteins, in each compartment, nucleus or cytoplasm, and between the two compartments. This article is part of a Special Issue entitled: Computational Proteomics, Systems Biology & Clinical Implications.

Design and optimization of reaction chamber and detection system in dynamic labs-on-chip for proteins detection.

In this paper, the lab-on-chip section for a protein assay is designed and optimized. To avoid severe reliability problems related to activated surface stability, a dynamic assay approach is adopted: protein-to-protein neutralization is performed while proteins diffuse freely in the reaction chamber. The related refraction index change is detected via an integrated interferometer. The structure is also design to provide a functional test of the reference protein solution, which is generally required for qualification for medical uses.

A spatial physiological model for p53 intracellular dynamics.

In this paper we design and analyse a physiologically based model representing the accumulation of protein p53 in the nucleus after triggering of ATM by DNA damage. The p53 protein is known to have a central role in the response of the cell to cytotoxic or radiotoxic insults resulting in DNA damage. A reasonable requirement for a model describing intracellular signalling pathways is taking into account the basic feature of eukaryotic cells: the distinction between nucleus and cytoplasm. Our aim is to show, on a simple reaction network describing p53 dynamics, how this basic distinction provides a framework which is able to yield expected oscillatory dynamics without introducing either positive feedbacks or delays in the reactions. Furthermore we prove that oscillations appear only if some spatial constraints are respected, e.g. if the diffusion coefficients correspond to known biological values. Finally we analyse how the spatial features of a cell influence the dynamic response of the p53 network to DNA damage, pointing out that the protein oscillatory dynamics is indeed a response that is robust towards changes with respect to cellular environments. Even if we change the cell shape or its volume or better its ribosomal distribution, we observe that DNA damage yields sustained oscillations of p53.

A pressure model of immune response to mycobacterium tuberculosis infection in several space dimensions.

Mycobacterium tuberculosis (Mtb) is a widely diffused infection. However, in general, the human immune system is able to contain it. In this work, we propose a mathematical model which describes the early immune response to the Mtb infection in the lungs, also including the possible evolution of the infection in the formation of a granuloma. The model is based on coupled reaction-diffusion-transport equations with chemotaxis, which take into account the interactions among bacteria, macrophages and chemoattractant. The novelty of this approach is in the modeling of the velocity field, proportional to the gradient of the pressure developed between the cells, which makes possible to deal with a full multidimensional description and efficient numerical simulations. We perform a linear stability analysis of the model and propose a robust implicit-explicit scheme to deal with long time simulations. Both in one and two-dimensions, we find that there are threshold values in the parameters space, between a contained infection and the uncontrolled bacteria growth, and the generation of granuloma-like patterns can be observed numerically.

A model of ischemia-induced neuroblast activation in the adult subventricular zone.

We have developed a rat brain organotypic culture model, in which tissue slices contain cortex-subventricular zone-striatum regions, to model neuroblast activity in response to in vitro ischemia. Neuroblast activation has been described in terms of two main parameters, proliferation and migration from the subventricular zone into the injured cortex. We observed distinct phases of neuroblast activation as is known to occur after in vivo ischemia. Thus, immediately after oxygen/glucose deprivation (6-24 hours), neuroblasts reduce their proliferative and migratory activity, whereas, at longer time points after the insult (2 to 5 days), they start to proliferate and migrate into the damaged cortex. Antagonism of ionotropic receptors for extracellular ATP during and after the insult unmasks an early activation of neuroblasts in the subventricular zone, which responded with a rapid and intense migration of neuroblasts into the damaged cortex (within 24 hours). The process is further enhanced by elevating the production of the chemoattractant SDf-1alpha and may also be boosted by blocking the activation of microglia. This organotypic model which we have developed is an excellent in vitro system to study neurogenesis after ischemia and other neurodegenerative diseases. Its application has revealed a SOS response to oxygen/glucose deprivation, which is inhibited by unfavorable conditions due to the ischemic environment. Finally, experimental quantifications have allowed us to elaborate a mathematical model to describe neuroblast activation and to develop a computer simulation which should have promising applications for the screening of drug candidates for novel therapies of ischemia-related pathologies.