Direct numerical simulations of interface dynamics to link capillary pressure and total surface energy

The flow of two immiscible fluids through a porous medium depends on the complex interplay between gravity, capillarity, and viscous forces. The interaction between these forces and the geometry of the medium gives rise to a variety of complex flow regimes that are difficult to describe using continuum models. Although a number of pore-scale models have been employed, a careful investigation of the macroscopic effects of pore-scale processes requires methods based on conservation principles in order to reduce the number of modeling assumptions. In this work we perform direct numerical simulations of drainage by solving Navier-Stokes equations in the pore space and employing the Volume Of Fluid (VOF) method to track the evolution of the fluid-fluid interface. After demonstrating that the method is able to deal with large viscosity contrasts and model the transition from stable flow to viscous fingering, we focus on the macroscopic capillary pressure and we compare different definitions of this quantity under quasi-static and dynamic conditions. We show that the difference between the intrinsic phase-average pressures, which is commonly used as definition of Darcy-scale capillary pressure, is subject to several limitations and it is not accurate in presence of viscous effects or trapping. In contrast, a definition based on the variation of the total surface energy provides an accurate estimate of the macroscopic capillary pressure. This definition, which links the capillary pressure to its physical origin, allows a better separation of viscous effects and does not depend on the presence of trapped fluid clusters.

Fate of the herbicide glyphosate and its metabolite AMPA in soils and their transfer to surface waters: A multi-scale approach in the Lavaux vineyards, western Switzerland

The use of herbicides in agriculture may lead to environmental problems, such as surface water pollution, with a potential risk for aquatic organisms. The herbicide glyphosate is the most used active ingredient in the world and in Switzerland. In the Lavaux vineyards it is nearly the only molecule applied. This work aimed at studying its fate in soils and its transfer to surface waters, using a multi-scale approach: from molecular (10-9 m) and microscopic scales (10-6 m), to macroscopic (m) and landscape ones (103 m). First of all, an analytical method was developed for the trace level quantification of this widely used herbicide and its main by-product, aminomethylphosphonic acid (AMPA). Due to their polar nature, their derivatization with 9-fluorenylmethyl chloroformate (FMOC-Cl) was done prior to their concentration and purification by solid phase extraction. They were then analyzed by ultra performance liquid chromatography coupled with tandem mass spectrometry (UPLC-MS/MS). The method was tested in different aqueous matrices with spiking tests and validated for the matrix effect correction in relevant environmental samples. Calibration curves established between 10 and 1000ng/l showed r2 values above 0.989, mean recoveries varied between 86 and 133% and limits of detection and quantification of the method were as low as 5 and 10ng/l respectively. At the parcel scale, two parcels of the Lavaux vineyard area, located near the Lutrive River at 6km to the east of Lausanne, were monitored to assess to which extent glyphosate and AMPA were retained in the soil or exported to surface waters. They were equipped at their bottom with porous ceramic cups and runoff collectors, which allowed retrieving water samples for the growing seasons 2010 and 2011. Results revealed that the mobility of glyphosate and AMPA in the unsaturated zone was likely driven by the precipitation regime and the soil characteristics, such as slope, porosity structure and layer permeability discrepancy. Elevated glyphosate and AMPA concentrations were measured at 60 and 80 cm depth at parcel bottoms, suggesting their infiltration in the upper parts of the parcels and the presence of preferential flow in the studied parcels. Indeed, the succession of rainy days induced the gradual saturation of the soil porosity, leading to rapid infiltration through macropores, as well as surface runoff formation. Furthermore, the presence of more impervious weathered marls at 100 cm depth induced throughflows, the importance of which for the lateral transport of the herbicide molecules was determined by the slope steepness. Important rainfall events (>10 mm/day) were clearly exporting molecules from the soil top layer, as indicated by important concentrations in runoff samples. A mass balance showed that total loss (10-20%) mainly occurred through surface runoff (96%) and, to a minor extent, by throughflows in soils (4%), with subsequent exfiltration to surface waters. Observations made in the Lutrive River revealed interesting details of glyphosate and AMPA dynamics in urbanized landscapes, such as the Lavaux vineyards. Indeed, besides their physical and chemical properties, herbicide dynamics at the catchment level strongly depend on application rates, precipitation regime, land use and also on the presence of drains or constructed channels. Elevated concentrations, up to 4970 ng/l, observed just after the application, confirmed the diffuse export of these compounds from the vineyard area by surface runoff during main rain events. From April to September 2011, a total load of 7.1 kg was calculated, with 85% coming from vineyards and minor urban sources and 15% from arable crops. Small vineyard surfaces could generate high concentrations of herbicides and contribute considerably to the total load calculated at the outlet, due to their steep slopes (~10%). The extrapolated total amount transferred yearly from the Lavaux vineyards to the Lake of Geneva was of 190kg. At the molecular scale, the possible involvement of dissolved organic matter (DOM) in glyphosate and copper transport was studied using UV/Vis fluorescence spectroscopy. Combined with parallel factor (PARAFAC) analysis, this technique allowed characterizing DOM of soil and surface water samples from the studied vineyard area. Glyphosate concentrations were linked to the fulvic-like spectroscopic signature of DOM in soil water samples, as well as to copper, suggesting the formation of ternary complexes. In surface water samples, its concentrations were also correlated to copper ones, but not in a significant way to the fulvic-like signature. Quenching experiments with standards confirmed field tendencies in the laboratory, with a stronger decrease in fluorescence intensity for fulvic-like fluorophore than for more aromatic ones. Lastly, based on maximum concentrations measured in the river, an environmental risk for these compounds was assessed, using laboratory tests and ecotoxicity data from the literature. In our case and with the methodology applied, the risk towards aquatic species was found negligible (RF<1).

Covalent cell surface functionalization of human fetal osteoblasts for tissue engineering.

The chemical functionalization of cell-surface proteins of human primary fetal bone cells with hydrophilic bioorthogonal intermediates was investigated. Toward this goal, chemical pathways were developed for click reaction-mediated coupling of alkyne derivatives with cellular azido-expressing proteins. The incorporation via a tetraethylene glycol linker of a dipeptide and a reporter biotin allowed the proof of concept for the introduction of cell-specific peptide ligands and to follow the reaction in living cells. Tuning the conditions of the click reaction resulted in chemical functionalization of living human fetal osteoblasts with excellent cell survival.

Assessment of diesel exhaust particulate exposure and surface characteristics in association with levels of oxidative stress biomarkers

Exposure to PM10 and PM2.5 (particulate matter with aerodynamic diameter smaller than 10 μm and 2.5 μm, respectively) is associated with a range of adverse health effects, including cancer, pulmonary and cardiovascular diseases. Surface characteristics (chemical reactivity, surface area) are considered of prime importance to understand the mechanisms which lead to harmful effects. A hypothetical mechanism to explain these adverse effects is the ability of components (organics, metal ions) adsorbed on these particles to generate Reactive Oxygen Species (ROS), and thereby to cause oxidative stress in biological systems (Donaldson et al., 2003). ROS can attack almost any cellular structure, like DNA or cellular membrane, leading to the formation of a wide variety of degradation products which can be used as a biomarker of oxidative stress. The aim of the present research project is to test whether there is a correlation between the exposure to Diesel Exhaust Particulate (DEP) and the oxidative stress status. For that purpose, a survey has been conducted in real occupational situations where workers were exposed to DEP (bus depots). Different exposure variables have been considered: - particulate number, size distribution and surface area (SMPS); - particulate mass - PM2.5 and PM4 (gravimetry); - elemental and organic carbon (coulometry); - total adsorbed heavy metals - iron, copper, manganese (atomic adsorption); - surface functional groups present on aerosols (Knudsen flow reactor). Several biomarkers of oxidative stress (8-hydroxy-2'-deoxyguanosine and several aldehydes) have been determined either in urine or serum of volunteers. Results obtained during the sampling campaign in several bus depots indicated that the occupational exposure to particulates in these places was rather low (40-50 μg/m3 for PM4). Bimodal size distributions were generally observed (5 μm and <1 μm). Surface characteristics of PM4 varied strongly, depending on the bus depot. They were usually characterized by high carbonyl and low acidic sites content. Among the different biomarkers which have been analyzed within the framework of this study, mean urinary levels of 8-hydroxy-2'-deoxyguanosine increased significantly (p<0.05) during two consecutive days of exposure for non-smoker workers. On the other hand, no statistically significant differences were observed for serum levels of hexanal, nonanal and 4- hydroxy-nonenal (p>0.05). Biomarkers levels will be compared to exposure variables to gain a better understanding of the relation between the particulate characteristics and the formation of ROS by-products. This project is financed by the Swiss State Secretariat for Education and Research. It is conducted within the framework of the COST Action 633 "Particulate Matter - Properties Related to Health Effects".

Surface modification of biomaterials for conjugation with human fetal osteoblasts.

Surface functionalization of hydroxyapatite (HA) and beta-tricalcium phosphate (TCP) bioceramics with chemical ligands containing a pyrrogallol moiety was developed to improve the adhesion of bone cell precursors to the biomaterials. Fast and biocompatible copper-free click reaction with azido-modified human fetal osteoblasts resulted in improved cell binding to both HA and TCP bioceramics, opening the way for using this methodology in the preparation of cell-engineered bone implants.

Simulation of surface waves in porous media

We present a novel numerical algorithm for the simulation of seismic wave propagation in porous media, which is particularly suitable for the accurate modelling of surface wave-type phenomena. The differential equations of motion are based on Biot's theory of poro-elasticity and solved with a pseudospectral approach using Fourier and Chebyshev methods to compute the spatial derivatives along the horizontal and vertical directions, respectively. The time solver is a splitting algorithm that accounts for the stiffness of the differential equations. Due to the Chebyshev operator the grid spacing in the vertical direction is non-uniform and characterized by a denser spatial sampling in the vicinity of interfaces, which allows for a numerically stable and accurate evaluation of higher order surface wave modes. We stretch the grid in the vertical direction to increase the minimum grid spacing and reduce the computational cost. The free-surface boundary conditions are implemented with a characteristics approach, where the characteristic variables are evaluated at zero viscosity. The same procedure is used to model seismic wave propagation at the interface between a fluid and porous medium. In this case, each medium is represented by a different grid and the two grids are combined through a domain-decomposition method. This wavefield decomposition method accounts for the discontinuity of variables and is crucial for an accurate interface treatment. We simulate seismic wave propagation with open-pore and sealed-pore boundary conditions and verify the validity and accuracy of the algorithm by comparing the numerical simulations to analytical solutions based on zero viscosity obtained with the Cagniard-de Hoop method. Finally, we illustrate the suitability of our algorithm for more complex models of porous media involving viscous pore fluids and strongly heterogeneous distributions of the elastic and hydraulic material properties.

Off-axis low-coherence interferometry for surface topology measurement.

In this communication we introduce a low or reduced coherence interferometry technique that can be used to retrieve surface topology on samples with high roughness. Moreover, we will show that the approach enables surface topology measurement also at the interface of so-called turbid media, where multiple scattering inside tissues can be a major issue, preventing accurate measurements.

Phorbol 12-myristate 13-acetate induces surface expression of T3 on human immature T cell lines with and without concomitant expression of the T cell antigen receptor complex.

The T3 complex is known to be expressed on the cell surface of mature T cells together with either the alpha-beta heterodimeric T cell receptor (TCR) or the TCR gamma protein. In a number of immature T cell malignancies, however, T3 has been described exclusively in the cytoplasm. We have investigated five such T cell lines with cytoplasmic T3 and could demonstrate by biosynthetic labeling the presence of the alpha and beta chains of the TCR in the cytoplasm of two of them, CEM and Ichikawa. No surface TCR alpha-beta protein could be detected by staining with the WT31 antibody. These observations, therefore, argue against the concept that expression of the TCR alpha chain controls the surface expression of the T3/TCR complex. Interestingly, phorbol 12-myristate 13-acetate (PMA) induced cell surface expression of T3 protein in these two cell lines only. Moreover, on surface-iodinated CEM cells no association of T3 and TCR molecules could be demonstrated after treatment with PMA, and expression of TCR alpha and beta chains was limited to the cytoplasm. In Ichikawa cells, however, PMA induced surface expression of a mature T3/TCR complex. Our findings indicate that separate regulatory mechanisms may exist for the surface expression of the T3 proteins and for the assembly of the T3/TCR complex.

Proteomic and transcriptomic profiling of Staphylococcus aureus surface LPXTG-proteins: correlation with agr genotypes and adherence phenotypes.

Staphylococcus aureus infections involve numerous adhesins and toxins, which expression depends on complex regulatory networks. Adhesins include a family of surface proteins covalently attached to the peptidoglycan via a conserved LPXTG motif. Here we determined the protein and mRNA expression of LPXTG-proteins of S. aureus Newman in time-course experiments, and their relation to fibrinogen adherence in vitro. Experiments were performed with mutants in the global accessory-gene regulator (agr), surface protein A (Spa), and fibrinogen-binding protein A (ClfA), as well as during growth in iron-rich or iron-poor media. Surface proteins were recovered by trypsin-shaving of live bacteria. Released peptides were analyzed by liquid chromatography coupled to tandem mass-spectrometry. To unambiguously identify peptides unique to LPXTG-proteins, the analytical conditions were refined using a reference library of S. aureus LPXTG-proteins heterogeneously expressed in surrogate Lactococcus lactis. Transcriptomes were determined by microarrays. Sixteen of the 18 LPXTG-proteins present in S. aureus Newman were detected by proteomics. Nine LPXTG-proteins showed a bell-shape agr-like expression that was abrogated in agr-negative mutants including Spa, fibronectin-binding protein A (FnBPA), ClfA, iron-binding IsdA, and IsdB, immunomodulator SasH, functionally uncharacterized SasD, biofilm-related SasG and methicillin resistance-related FmtB. However, only Spa and SasH modified their proteomic and mRNA profiles in parallel in the parent and its agr- mutant, whereas all other LPXTG-proteins modified their proteomic profiles independently of their mRNA. Moreover, ClfA became highly transcribed and active in fibrinogen-adherence tests during late growth (24 h), whereas it remained poorly detected by proteomics. On the other hand, iron-regulated IsdA-B-C increased their protein expression by >10-times in iron-poor conditions. Thus, proteomic, transcriptomic, and adherence-phenotype demonstrated differential profiles in S. aureus. Moreover, trypsin peptide signatures suggested differential protein domain exposures in various environments, which might be relevant for anti-adhesin vaccines. A comprehensive understanding of the S. aureus physiology should integrate all three approaches.

Gene transfer of the Ly-3 chain gene of the mouse CD8 molecular complex: co-transfer with the Ly-2 polypeptide gene results in detectable cell surface expression of the Ly-3 antigenic determinants.

The CD8 molecule is a glycoprotein expressed on a subset of mature T lymphocytes. It has been postulated to be a receptor for class I major histocompatibility complex molecules. In the mouse, CD8 is a heterodimer composed of Ly-2 and Ly-3 chains. We have isolated and analyzed cDNA and cosmid clones corresponding to the Ly-3 subunit. One of the isolated, cosmid clones was subsequently transfected, alone or in combination with the Ly-2 gene, into mouse Ltk- cells. Analysis of the Ly-2,3 molecules expressed at the surface of the double transfectants indicated that they are serologically and biochemically indistinguishable from their normal counterparts expressed on lymphoid cells. Ltk- cells transfected with the Ly-2 gene alone were shown to react with a subset of anti-CD8 monoclonal antibodies whereas Ly-3 transfectants did not stain with any of the anti-Ly-3 antibodies employed in this study. Since at least one of these antibodies (53-5.8) has been previously shown to recognize an epitope which is retained on the Ly-3 subunit after dissociation of the heterodimeric Ly-2,3 complex, these observations suggest that the expression of the Ly-2 polypeptide is required to permit the detectable cell surface expression of the antigenic determinants carried by the Ly-3 subunit.

Structure of the glycosyl-phosphatidylinositol membrane anchor of the Leishmania major promastigote surface protease.

In common with many other plasma membrane glycoproteins of eukaryotic origin, the promastigote surface protease (PSP) of the protozoan parasite Leishmania contains a glycosyl-phosphatidylinositol (GPI) membrane anchor. The GPI anchor of Leishmania major PSP was purified following proteolysis of the PSP and analyzed by two-dimensional 1H-1H NMR, compositional and methylation linkage analyses, chemical and enzymatic modifications, and amino acid sequencing. From these results, the structure of the GPI-containing peptide was found to be Asp-Gly-Gly-Asn-ethanolamine-PO4-6Man alpha 1-6Man alpha 1-4GlcN alpha 1-6myo-inositol-1-PO4-(1-alkyl-2-acyl-glycerol). The glycan structure is identical to the conserved glycan core regions of the GPI anchor of Trypanosoma brucei variant surface glycoprotein and rat brain Thy-1 antigen, supporting the notion that this portion of GPIs are highly conserved. The phosphatidylinositol moiety of the PSP anchor is unusual, containing a fully saturated, unbranched 1-O-alkyl chain (mainly C24:0) and a mixture of fully saturated unbranched 2-O-acyl chains (C12:0, C14:0, C16:0, and C18:0). This lipid composition differs significantly from those of the GPIs of T. brucei variant surface glycoprotein and mammalian erythrocyte acetylcholinesterase but is similar to that of a family of glycosylated phosphoinositides found uniquely in Leishmania.

β1- and β3- voltage-gated sodium channel subunits modulate cell surface expression and glycosylation of Nav1.7 in HEK293 cells.

Voltage-gated sodium channels (Navs) are glycoproteins composed of a pore-forming α-subunit and associated β-subunits that regulate Nav α-subunit plasma membrane density and biophysical properties. Glycosylation of the Nav α-subunit also directly affects Navs gating. β-subunits and glycosylation thus comodulate Nav α-subunit gating. We hypothesized that β-subunits could directly influence α-subunit glycosylation. Whole-cell patch clamp of HEK293 cells revealed that both β1- and β3-subunits coexpression shifted V ½ of steady-state activation and inactivation and increased Nav1.7-mediated I Na density. Biotinylation of cell surface proteins, combined with the use of deglycosydases, confirmed that Nav1.7 α-subunits exist in multiple glycosylated states. The α-subunit intracellular fraction was found in a core-glycosylated state, migrating at ~250 kDa. At the plasma membrane, in addition to the core-glycosylated form, a fully glycosylated form of Nav1.7 (~280 kDa) was observed. This higher band shifted to an intermediate band (~260 kDa) when β1-subunits were coexpressed, suggesting that the β1-subunit promotes an alternative glycosylated form of Nav1.7. Furthermore, the β1-subunit increased the expression of this alternative glycosylated form and the β3-subunit increased the expression of the core-glycosylated form of Nav1.7. This study describes a novel role for β1- and β3-subunits in the modulation of Nav1.7 α-subunit glycosylation and cell surface expression.

Brain surface non-rigid registration in auditory processing

Résumé: Le développement rapide de nouvelles technologies comme l'imagerie médicale a permis l'expansion des études sur les fonctions cérébrales. Le rôle principal des études fonctionnelles cérébrales est de comparer l'activation neuronale entre différents individus. Dans ce contexte, la variabilité anatomique de la taille et de la forme du cerveau pose un problème majeur. Les méthodes actuelles permettent les comparaisons interindividuelles par la normalisation des cerveaux en utilisant un cerveau standard. Les cerveaux standards les plus utilisés actuellement sont le cerveau de Talairach et le cerveau de l'Institut Neurologique de Montréal (MNI) (SPM99). Les méthodes de recalage qui utilisent le cerveau de Talairach, ou celui de MNI, ne sont pas suffisamment précises pour superposer les parties plus variables d'un cortex cérébral (p.ex., le néocortex ou la zone perisylvienne), ainsi que les régions qui ont une asymétrie très importante entre les deux hémisphères. Le but de ce projet est d'évaluer une nouvelle technique de traitement d'images basée sur le recalage non-rigide et utilisant les repères anatomiques. Tout d'abord, nous devons identifier et extraire les structures anatomiques (les repères anatomiques) dans le cerveau à déformer et celui de référence. La correspondance entre ces deux jeux de repères nous permet de déterminer en 3D la déformation appropriée. Pour les repères anatomiques, nous utilisons six points de contrôle qui sont situés : un sur le gyrus de Heschl, un sur la zone motrice de la main et le dernier sur la fissure sylvienne, bilatéralement. Evaluation de notre programme de recalage est accomplie sur les images d'IRM et d'IRMf de neuf sujets parmi dix-huit qui ont participés dans une étude précédente de Maeder et al. Le résultat sur les images anatomiques, IRM, montre le déplacement des repères anatomiques du cerveau à déformer à la position des repères anatomiques de cerveau de référence. La distance du cerveau à déformer par rapport au cerveau de référence diminue après le recalage. Le recalage des images fonctionnelles, IRMf, ne montre pas de variation significative. Le petit nombre de repères, six points de contrôle, n'est pas suffisant pour produire les modifications des cartes statistiques. Cette thèse ouvre la voie à une nouvelle technique de recalage du cortex cérébral dont la direction principale est le recalage de plusieurs points représentant un sillon cérébral. Abstract : The fast development of new technologies such as digital medical imaging brought to the expansion of brain functional studies. One of the methodolgical key issue in brain functional studies is to compare neuronal activation between individuals. In this context, the great variability of brain size and shape is a major problem. Current methods allow inter-individual comparisions by means of normalisation of subjects' brains in relation to a standard brain. A largerly used standard brains are the proportional grid of Talairach and Tournoux and the Montreal Neurological Insititute standard brain (SPM99). However, there is a lack of more precise methods for the superposition of more variable portions of the cerebral cortex (e.g, neocrotex and perisyvlian zone) and in brain regions highly asymmetric between the two cerebral hemipsheres (e.g. planum termporale). The aim of this thesis is to evaluate a new image processing technique based on non-linear model-based registration. Contrary to the intensity-based, model-based registration uses spatial and not intensitiy information to fit one image to another. We extract identifiable anatomical features (point landmarks) in both deforming and target images and by their correspondence we determine the appropriate deformation in 3D. As landmarks, we use six control points that are situated: one on the Heschl'y Gyrus, one on the motor hand area, and one on the sylvian fissure, bilaterally. The evaluation of this model-based approach is performed on MRI and fMRI images of nine of eighteen subjects participating in the Maeder et al. study. Results on anatomical, i.e. MRI, images, show the mouvement of the deforming brain control points to the location of the reference brain control points. The distance of the deforming brain to the reference brain is smallest after the registration compared to the distance before the registration. Registration of functional images, i.e fMRI, doesn't show a significant variation. The small number of registration landmarks, i.e. six, is obvious not sufficient to produce significant modification on the fMRI statistical maps. This thesis opens the way to a new computation technique for cortex registration in which the main directions will be improvement of the registation algorithm, using not only one point as landmark, but many points, representing one particular sulcus.

Surface functionalization of alumina ceramic foams with organic ligands.

Different anchoring groups have been studied with the aim of covalently binding organic linkers to the surface of alumina ceramic foams. The results suggested that a higher degree of functionalization was achieved with a pyrogallol derivative - as compared to its catechol analogue - based on the XPS analysis of the ceramic surface. The conjugation of organic ligands to the surface of these alumina materials was corroborated by DNP-MAS NMR measurements.

Quantification of the relation between surface morphodynamics and subsurface sedimentological product in sandy braided rivers

This study uses digital elevation models and ground-penetrating radar to quantify the relation between the surface morphodynamics and subsurface sedimentology in the sandy braided South Saskatchewan River, Canada. A unique aspect of the methodology is that both digital elevation model and ground-penetrating radar data were collected from the same locations in 2004, 2005, 2006 and 2007, thus enabling the surface morphodynamics to be tied explicitly to the associated evolving depositional product. The occurrence of a large flood in 2005 also allowed the influence of discharge to be assessed with respect to the processproduct relationship. The data demonstrate that the morphology of the study reach evolved even during modest discharges, but more extensive erosion was caused by the large flood. In addition, the study reach was dominated by compound bars before the flood, but switched to being dominated by unit bars during and after the flood. The extent to which the subsurface deposits (the product') were modified by the surface morphodynamics (the process') was quantified using the changes in radar-facies recorded in sequential ground-penetrating radar surveys. These surveys reveal that during the large flood there was an increase in the proportion of facies associated with bar margin accretion and larger dunes. In subsequent years, these facies became truncated and replaced with facies associated with smaller dune sets. This analysis shows that unit bars generally become truncated more laterally than vertically and, thus, they lose the high-angle bar margin deposits and smaller scale bar-top deposits. In general, the only fragments that remain of the unit bars are dune sets, thus making identification of the original unit barform problematic. This novel data set has implications for what may ultimately become preserved in the rock record.