Transfer of ultrasmall iron oxide nanoparticles from human brain-derived endothelial cells to human glioblastoma cells.

Nanoparticles (NPs) are being used or explored for the development of biomedical applications in diagnosis and therapy, including imaging and drug delivery. Therefore, reliable tools are needed to study the behavior of NPs in biological environment, in particular the transport of NPs across biological barriers, including the blood-brain tumor barrier (BBTB), a challenging question. Previous studies have addressed the translocation of NPs of various compositions across cell layers, mostly using only one type of cells. Using a coculture model of the human BBTB, consisting in human cerebral endothelial cells preloaded with ultrasmall superparamagnetic iron oxide nanoparticles (USPIO NPs) and unloaded human glioblastoma cells grown on each side of newly developed ultrathin permeable silicon nitride supports as a model of the human BBTB, we demonstrate for the first time the transfer of USPIO NPs from human brain-derived endothelial cells to glioblastoma cells. The reduced thickness of the permeable mechanical support compares better than commercially available polymeric supports to the thickness of the basement membrane of the cerebral vascular system. These results are the first report supporting the possibility that USPIO NPs could be directly transferred from endothelial cells to glioblastoma cells across a BBTB. Thus, the use of such ultrathin porous supports provides a new in vitro approach to study the delivery of nanotherapeutics to brain cancers. Our results also suggest a novel possibility for nanoparticles to deliver therapeutics to the brain using endothelial to neural cells transfer.

Hybrid Multiscale Finite Volume method for two-phase flow in porous media

We present a novel hybrid (or multiphysics) algorithm, which couples pore-scale and Darcy descriptions of two-phase flow in porous media. The flow at the pore-scale is described by the Navier?Stokes equations, and the Volume of Fluid (VOF) method is used to model the evolution of the fluid?fluid interface. An extension of the Multiscale Finite Volume (MsFV) method is employed to construct the Darcy-scale problem. First, a set of local interpolators for pressure and velocity is constructed by solving the Navier?Stokes equations; then, a coarse mass-conservation problem is constructed by averaging the pore-scale velocity over the cells of a coarse grid, which act as control volumes; finally, a conservative pore-scale velocity field is reconstructed and used to advect the fluid?fluid interface. The method relies on the localization assumptions used to compute the interpolators (which are quite straightforward extensions of the standard MsFV) and on the postulate that the coarse-scale fluxes are proportional to the coarse-pressure differences. By numerical simulations of two-phase problems, we demonstrate that these assumptions provide hybrid solutions that are in good agreement with reference pore-scale solutions and are able to model the transition from stable to unstable flow regimes. Our hybrid method can naturally take advantage of several adaptive strategies and allows considering pore-scale fluxes only in some regions, while Darcy fluxes are used in the rest of the domain. Moreover, since the method relies on the assumption that the relationship between coarse-scale fluxes and pressure differences is local, it can be used as a numerical tool to investigate the limits of validity of Darcy's law and to understand the link between pore-scale quantities and their corresponding Darcy-scale variables.

Error estimate and control in the MSFV method for multiphase flow in porous media

n this paper the iterative MSFV method is extended to include the sequential implicit simulation of time dependent problems involving the solution of a system of pressure-saturation equations. To control numerical errors in simulation results, an error estimate, based on the residual of the MSFV approximate pressure field, is introduced. In the initial time steps in simulation iterations are employed until a specified accuracy in pressure is achieved. This initial solution is then used to improve the localization assumption at later time steps. Additional iterations in pressure solution are employed only when the pressure residual becomes larger than a specified threshold value. Efficiency of the strategy and the error control criteria are numerically investigated. This paper also shows that it is possible to derive an a-priori estimate and control based on the allowed pressure-equation residual to guarantee the desired accuracy in saturation calculation.

Use of quantum dots in aqueous solution to detect blood fingermarks on non-porous surfaces.

A new and original reagent based on the use of highly fluorescent cadmium telluride (CdTe) quantum dots (QDs) in aqueous solution is proposed to detect weak fingermarks in blood on non-porous surfaces. To assess the efficiency of this approach, comparisons were performed with one of the most efficient blood reagents on non-porous surfaces, Acid Yellow 7 (AY7). To this end, four non-porous surfaces were studied, i.e. glass, transparent polypropylene, black polyethylene, and aluminium foil. To evaluate the sensitivity of both reagents, sets of depleted fingermarks were prepared, using the same finger, initially soaked with blood, which was then successively applied on the same surface without recharging it with blood or latent secretions. The successive marks were then cut in halves and the halves treated separately with each reagent. The results showed that QDs were equally efficient to AY7 on glass, polyethylene and polypropylene surfaces, and were superior to AY7 on aluminium. The use of QDs in new, sensitive and highly efficient latent and blood mark detection techniques appears highly promising. Health and safety issues related to the use of cadmium are also discussed. It is suggested that applying QDs in aqueous solution (and not as a dry dusting powder) considerably lowers the toxicity risks.

Dual role of Bmi1 loss in the preservation of photoreceptor layers in the Rd1 mouse

Purpose: In the Rd1 and Rd10 mouse models of retinitis pigmentosa, a mutation in the Pde6ß gene leads to the rapid loss of photoreceptors. As in several neurodegenerative diseases, Rd1 and Rd10 photoreceptors re-express cell cycle proteins prior to death. Bmi1 regulates cell cycle progression through inhibition of CDK inhibitors, and its deletion efficiently rescues the Rd1 retinal degeneration. The present study evaluates the effects of Bmi1 loss in photoreceptors and Müller glia, since in lower vertebrates, these cells respond to retinal injury through dedifferentiation and regeneration of retinal cells. Methods: Cell death and Müller cell activation were analyzed by immunostaining of wild-type, Rd1 and Rd1;Bmi1-/- eye sections during retinal degeneration, between P10 and P20. Lineage tracing experiments use the GFAP-Cre mouse (JAX) to target Müller cells. Results: In Rd1 retinal explants, inhibition of CDKs reduces the amount of dying cells. In vivo, Bmi1 deletion reduces CDK4 expression and cell death in the P15 Rd1;Bmi1-/- retina, although cGMP accumulation and TUNEL staining are detected at the onset of retinal degeneration (P12). This suggests that another process acts in parallel to overcome the initial loss of Rd1;Bmi1-/- photoreceptors. We demonstrate here that Bmi1 loss in the Rd1 retina enhances the activation of Müller glia by downregulation of p27Kip1, that these cells migrate toward the ONL, and that some cells express the retinal progenitor marker Pax6 at the inner part of the ONL. These events are also observed, but to a lesser extent, in Rd1 and Rd10 retinas. At P12, EdU incorporation shows proliferating cells with atypical elongated nuclei at the inner border of the Rd1;Bmi1-/- ONL. Lineage tracing targeting Müller cells is in process and will determine the implication of this cell population in the maintenance of the Rd1;Bmi1-/- ONL thickness and whether downregulation of Bmi1 in Rd10 Müller cells equally stimulates their activation. Conclusions: Our results show a dual role of Bmi1 deletion in the rescue of photoreceptors in the Rd1;Bmi1-/- retina. Indeed, the loss of Bmi1 reduces Rd1 retinal degeneration, and as well, enhances the Müller glia activation. In addition, the emergence of cells expressing a retinal progenitor marker in the ONL suggests Bmi1 as a blockade to the regeneration of retinal cells in mammals.

Functionalization of microstructured open-porous bioceramic scaffolds with human fetal bone cells.

Bone substitute materials allowing trans-scaffold migration and in-scaffold survival of human bone-derived cells are mandatory for development of cell-engineered permanent implants to repair bone defects. In this study, we evaluated the influence on human bone-derived cells of the material composition and microstructure of foam scaffolds of calcium aluminate. The scaffolds were prepared using a direct foaming method allowing wide-range tailoring of the microstructure for pore size and pore openings. Human fetal osteoblasts (osteo-progenitors) attached to the scaffolds, migrated across the entire bioceramic depending on the scaffold pore size, colonized, and survived in the porous material for at least 6 weeks. The long-term biocompatibility of the scaffold material for human bone-derived cells was evidenced by in-scaffold determination of cell metabolic activity using a modified MTT assay, a repeated WST-1 assay, and scanning electron microscopy. Finally, we demonstrated that the osteo-progenitors can be covalently bound to the scaffolds using biocompatible click chemistry, thus enhancing the rapid adhesion of the cells to the scaffolds. Therefore, the different microstructures of the foams influenced the migratory potential of the cells, but not cell viability. Scaffolds allow covalent biocompatible chemical binding of the cells to the materials, either localized or widespread integration of the scaffolds for cell-engineered implants.

Seismic attenuation and velocity dispersion in heterogeneous partially saturated porous rocks

Using a numerical approach, we explore wave-induced fluid flow effects in partially saturated porous rocks in which the gas-water saturation patterns are governed by mesoscopic heterogeneities associated with the dry frame properties. The link between the dry frame properties and the gas saturation is defined by the assumption of capillary pressure equilibrium, which in the presence of heterogeneity implies that neighbouring regions can exhibit different levels of saturation. To determine the equivalent attenuation and phase velocity of the synthetic rock samples considered in this study, we apply a numerical upscaling procedure, which permits to take into account mesoscopic heterogeneities associated with the dry frame properties as well as spatially continuous variations of the pore fluid properties. The multiscale nature of the fluid saturation is taken into account by locally computing the physical properties of an effective fluid, which are then used for the larger-scale simulations. We consider two sets of numerical experiments to analyse such effects in heterogeneous partially saturated porous media, where the saturation field is determined by variations in porosity and clay content, respectively. In both cases we also evaluate the seismic responses of corresponding binary, patchy-type saturation patterns. Our results indicate that significant attenuation and modest velocity dispersion effects take place in this kind of media for both binary patchy-type and spatially continuous gas saturation patterns and in particular in the presence of relatively small amounts of gas. The numerical experiments also show that the nature of the gas distribution patterns is a critical parameter controlling the seismic responses of these environments, since attenuation and velocity dispersion effects are much more significant and occur over a broader saturation range for binary patchy-type gas-water distributions. This analysis therefore suggests that the physical mechanisms governing partial saturation should be accounted for when analysing seismic data in a poroelastic framework. In this context, heterogeneities associated with the dry frame properties, which do not play important roles in wave-induced fluid flow processes per se, should be taken into account since they may determine the kind of gas distribution pattern taking place in the porous rock.

Stoneley wave modeling in heterogeneous porous media with viscous pore fluids

We implemented Biot-type porous wave equations in a pseudo-spectral numerical modeling algorithm for the simulation of Stoneley waves in porous media. Fourier and Chebyshev methods are used to compute the spatial derivatives along the horizontal and vertical directions, respectively. To prevent from overly short time steps due to the small grid spacing at the top and bottom of the model as a consequence of the Chebyshev operator, the mesh is stretched in the vertical direction. As a large benefit, the Chebyshev operator allows for an explicit treatment of interfaces. Boundary conditions can be implemented with a characteristics approach. The characteristic variables are evaluated at zero viscosity. We use this approach to model seismic wave propagation at the interface between a fluid and a porous medium. Each medium is represented by a different mesh and the two meshes are connected through the above described characteristics domain-decomposition method. We show an experiment for sealed pore boundary conditions, where we first compare the numerical solution to an analytical solution. We then show the influence of heterogeneity and viscosity of the pore fluid on the propagation of the Stoneley wave and surface waves in general.

A comparison of frailty of primary neurons, embryonic, and aging mouse cortical layers.

Superficial layers I to III of the human cerebral cortex are more vulnerable toward Aβ peptides than deep layers V to VI in aging. Three models of layers were used to investigate this pattern of frailty. First, primary neurons from E14 and E17 embryonic murine cortices, corresponding respectively to future deep and superficial layers, were treated either with Aβ1-42, okadaic acid, or kainic acid. Second, whole E14 and E17 embryonic cortices, and third, in vitro separated deep and superficial layers of young and old C57BL/6J mice, were treated identically. We observed that E14 and E17 neurons in culture were prone to death after the Aβ and particularly the kainic acid treatment. This was also the case for the superficial layers of the aged cortex, but not for the embryonic, the young cortex, and the deep layers of the aged cortex. Thus, the aged superficial layers appeared to be preferentially vulnerable against Aβ and kainic acid. This pattern of vulnerability corresponds to enhanced accumulation of senile plaques in the superficial cortical layers with aging and Alzheimer's disease.

Frequency-dependent attenuation as a potential indicator of oil saturation

At seismic frequencies, wave-induced fluid flow is a major cause of P-wave attenuation in partially saturated porous rocks. Attenuation is of great importance for the oil industry in the interpretation of seismic field data. Here, the effects on P-wave attenuation resulting from changes in oil saturation are studied for media with coexisting water, oil, and gas. For that, creep experiments are numerically simulated by solving Biot's equations for consolidation of poroelastic media with the finite-element method. The experiments yield time-dependent stress?strain relations that are used to calculate the complex P-wave modulus from which frequency-dependent P-wave attenuation is determined. The models are layered media with periodically alternating triplets of layers. Models consisting of triplets of layers having randomly varying layer thicknesses are also considered. The layers in each triplet are fully saturated with water, oil, and gas. The layer saturated with water has lower porosity and permeability than the layers saturated with oil and gas. These models represent hydrocarbon reservoirs in which water is the wetting fluid preferentially saturating regions of lower porosity. The results from the numerical experiments showed that increasing oil saturation, connected to a decrease in gas saturation, resulted in a significant increase of attenuation at low frequencies (lower than 2 Hz). Furthermore, replacing the oil with water resulted in a distinguishable behavior of the frequency-dependent attenuation. These results imply that, according to the physical mechanism of wave-induced fluid flow, frequency-dependent attenuation in media saturated with water, oil, and gas is a potential indicator of oil saturation.

Multiscale descriptions of density-driven flow instabilities in porous media

Les instabilités engendrées par des gradients de densité interviennent dans une variété d'écoulements. Un exemple est celui de la séquestration géologique du dioxyde de carbone en milieux poreux. Ce gaz est injecté à haute pression dans des aquifères salines et profondes. La différence de densité entre la saumure saturée en CO2 dissous et la saumure environnante induit des courants favorables qui le transportent vers les couches géologiques profondes. Les gradients de densité peuvent aussi être la cause du transport indésirable de matières toxiques, ce qui peut éventuellement conduire à la pollution des sols et des eaux. La gamme d'échelles intervenant dans ce type de phénomènes est très large. Elle s'étend de l'échelle poreuse où les phénomènes de croissance des instabilités s'opèrent, jusqu'à l'échelle des aquifères à laquelle interviennent les phénomènes à temps long. Une reproduction fiable de la physique par la simulation numérique demeure donc un défi en raison du caractère multi-échelles aussi bien au niveau spatial et temporel de ces phénomènes. Il requiert donc le développement d'algorithmes performants et l'utilisation d'outils de calculs modernes. En conjugaison avec les méthodes de résolution itératives, les méthodes multi-échelles permettent de résoudre les grands systèmes d'équations algébriques de manière efficace. Ces méthodes ont été introduites comme méthodes d'upscaling et de downscaling pour la simulation d'écoulements en milieux poreux afin de traiter de fortes hétérogénéités du champ de perméabilité. Le principe repose sur l'utilisation parallèle de deux maillages, le premier est choisi en fonction de la résolution du champ de perméabilité (grille fine), alors que le second (grille grossière) est utilisé pour approximer le problème fin à moindre coût. La qualité de la solution multi-échelles peut être améliorée de manière itérative pour empêcher des erreurs trop importantes si le champ de perméabilité est complexe. Les méthodes adaptatives qui restreignent les procédures de mise à jour aux régions à forts gradients permettent de limiter les coûts de calculs additionnels. Dans le cas d'instabilités induites par des gradients de densité, l'échelle des phénomènes varie au cours du temps. En conséquence, des méthodes multi-échelles adaptatives sont requises pour tenir compte de cette dynamique. L'objectif de cette thèse est de développer des algorithmes multi-échelles adaptatifs et efficaces pour la simulation des instabilités induites par des gradients de densité. Pour cela, nous nous basons sur la méthode des volumes finis multi-échelles (MsFV) qui offre l'avantage de résoudre les phénomènes de transport tout en conservant la masse de manière exacte. Dans la première partie, nous pouvons démontrer que les approximations de la méthode MsFV engendrent des phénomènes de digitation non-physiques dont la suppression requiert des opérations de correction itératives. Les coûts de calculs additionnels de ces opérations peuvent toutefois être compensés par des méthodes adaptatives. Nous proposons aussi l'utilisation de la méthode MsFV comme méthode de downscaling: la grille grossière étant utilisée dans les zones où l'écoulement est relativement homogène alors que la grille plus fine est utilisée pour résoudre les forts gradients. Dans la seconde partie, la méthode multi-échelle est étendue à un nombre arbitraire de niveaux. Nous prouvons que la méthode généralisée est performante pour la résolution de grands systèmes d'équations algébriques. Dans la dernière partie, nous focalisons notre étude sur les échelles qui déterminent l'évolution des instabilités engendrées par des gradients de densité. L'identification de la structure locale ainsi que globale de l'écoulement permet de procéder à un upscaling des instabilités à temps long alors que les structures à petite échelle sont conservées lors du déclenchement de l'instabilité. Les résultats présentés dans ce travail permettent d'étendre les connaissances des méthodes MsFV et offrent des formulations multi-échelles efficaces pour la simulation des instabilités engendrées par des gradients de densité. - Density-driven instabilities in porous media are of interest for a wide range of applications, for instance, for geological sequestration of CO2, during which CO2 is injected at high pressure into deep saline aquifers. Due to the density difference between the C02-saturated brine and the surrounding brine, a downward migration of CO2 into deeper regions, where the risk of leakage is reduced, takes place. Similarly, undesired spontaneous mobilization of potentially hazardous substances that might endanger groundwater quality can be triggered by density differences. Over the last years, these effects have been investigated with the help of numerical groundwater models. Major challenges in simulating density-driven instabilities arise from the different scales of interest involved, i.e., the scale at which instabilities are triggered and the aquifer scale over which long-term processes take place. An accurate numerical reproduction is possible, only if the finest scale is captured. For large aquifers, this leads to problems with a large number of unknowns. Advanced numerical methods are required to efficiently solve these problems with today's available computational resources. Beside efficient iterative solvers, multiscale methods are available to solve large numerical systems. Originally, multiscale methods have been developed as upscaling-downscaling techniques to resolve strong permeability contrasts. In this case, two static grids are used: one is chosen with respect to the resolution of the permeability field (fine grid); the other (coarse grid) is used to approximate the fine-scale problem at low computational costs. The quality of the multiscale solution can be iteratively improved to avoid large errors in case of complex permeability structures. Adaptive formulations, which restrict the iterative update to domains with large gradients, enable limiting the additional computational costs of the iterations. In case of density-driven instabilities, additional spatial scales appear which change with time. Flexible adaptive methods are required to account for these emerging dynamic scales. The objective of this work is to develop an adaptive multiscale formulation for the efficient and accurate simulation of density-driven instabilities. We consider the Multiscale Finite-Volume (MsFV) method, which is well suited for simulations including the solution of transport problems as it guarantees a conservative velocity field. In the first part of this thesis, we investigate the applicability of the standard MsFV method to density- driven flow problems. We demonstrate that approximations in MsFV may trigger unphysical fingers and iterative corrections are necessary. Adaptive formulations (e.g., limiting a refined solution to domains with large concentration gradients where fingers form) can be used to balance the extra costs. We also propose to use the MsFV method as downscaling technique: the coarse discretization is used in areas without significant change in the flow field whereas the problem is refined in the zones of interest. This enables accounting for the dynamic change in scales of density-driven instabilities. In the second part of the thesis the MsFV algorithm, which originally employs one coarse level, is extended to an arbitrary number of coarse levels. We prove that this keeps the MsFV method efficient for problems with a large number of unknowns. In the last part of this thesis, we focus on the scales that control the evolution of density fingers. The identification of local and global flow patterns allows a coarse description at late times while conserving fine-scale details during onset stage. Results presented in this work advance the understanding of the Multiscale Finite-Volume method and offer efficient dynamic multiscale formulations to simulate density-driven instabilities. - Les nappes phréatiques caractérisées par des structures poreuses et des fractures très perméables représentent un intérêt particulier pour les hydrogéologues et ingénieurs environnementaux. Dans ces milieux, une large variété d'écoulements peut être observée. Les plus communs sont le transport de contaminants par les eaux souterraines, le transport réactif ou l'écoulement simultané de plusieurs phases non miscibles, comme le pétrole et l'eau. L'échelle qui caractérise ces écoulements est définie par l'interaction de l'hétérogénéité géologique et des processus physiques. Un fluide au repos dans l'espace interstitiel d'un milieu poreux peut être déstabilisé par des gradients de densité. Ils peuvent être induits par des changements locaux de température ou par dissolution d'un composé chimique. Les instabilités engendrées par des gradients de densité revêtent un intérêt particulier puisque qu'elles peuvent éventuellement compromettre la qualité des eaux. Un exemple frappant est la salinisation de l'eau douce dans les nappes phréatiques par pénétration d'eau salée plus dense dans les régions profondes. Dans le cas des écoulements gouvernés par les gradients de densité, les échelles caractéristiques de l'écoulement s'étendent de l'échelle poreuse où les phénomènes de croissance des instabilités s'opèrent, jusqu'à l'échelle des aquifères sur laquelle interviennent les phénomènes à temps long. Etant donné que les investigations in-situ sont pratiquement impossibles, les modèles numériques sont utilisés pour prédire et évaluer les risques liés aux instabilités engendrées par les gradients de densité. Une description correcte de ces phénomènes repose sur la description de toutes les échelles de l'écoulement dont la gamme peut s'étendre sur huit à dix ordres de grandeur dans le cas de grands aquifères. Il en résulte des problèmes numériques de grande taille qui sont très couteux à résoudre. Des schémas numériques sophistiqués sont donc nécessaires pour effectuer des simulations précises d'instabilités hydro-dynamiques à grande échelle. Dans ce travail, nous présentons différentes méthodes numériques qui permettent de simuler efficacement et avec précision les instabilités dues aux gradients de densité. Ces nouvelles méthodes sont basées sur les volumes finis multi-échelles. L'idée est de projeter le problème original à une échelle plus grande où il est moins coûteux à résoudre puis de relever la solution grossière vers l'échelle de départ. Cette technique est particulièrement adaptée pour résoudre des problèmes où une large gamme d'échelle intervient et évolue de manière spatio-temporelle. Ceci permet de réduire les coûts de calculs en limitant la description détaillée du problème aux régions qui contiennent un front de concentration mobile. Les aboutissements sont illustrés par la simulation de phénomènes tels que l'intrusion d'eau salée ou la séquestration de dioxyde de carbone.

Atomic Force Microscopy Stiffness Tomography on Living Arabidopsis thaliana Cells Reveals the Mechanical Properties of Surface and Deep Cell-Wall Layers during Growth.

Cell-wall mechanical properties play a key role in the growth and the protection of plants. However, little is known about genuine wall mechanical properties and their growth-related dynamics at subcellular resolution and in living cells. Here, we used atomic force microscopy (AFM) stiffness tomography to explore stiffness distribution in the cell wall of suspension-cultured Arabidopsis thaliana as a model of primary, growing cell wall. For the first time that we know of, this new imaging technique was performed on living single cells of a higher plant, permitting monitoring of the stiffness distribution in cell-wall layers as a function of the depth and its evolution during the different growth phases. The mechanical measurements were correlated with changes in the composition of the cell wall, which were revealed by Fourier-transform infrared (FTIR) spectroscopy. In the beginning and end of cell growth, the average stiffness of the cell wall was low and the wall was mechanically homogenous, whereas in the exponential growth phase, the average wall stiffness increased, with increasing heterogeneity. In this phase, the difference between the superficial and deep wall stiffness was highest. FTIR spectra revealed a relative increase in the polysaccharide/lignin content.