Numerical simulation of gel electrophoresis of DNA knots in weak and strong electric fields.

Gel electrophoresis allows one to separate knotted DNA (nicked circular) of equal length according to the knot type. At low electric fields, complex knots, being more compact, drift faster than simpler knots. Recent experiments have shown that the drift velocity dependence on the knot type is inverted when changing from low to high electric fields. We present a computer simulation on a lattice of a closed, knotted, charged DNA chain drifting in an external electric field in a topologically restricted medium. Using a Monte Carlo algorithm, the dependence of the electrophoretic migration of the DNA molecules on the knot type and on the electric field intensity is investigated. The results are in qualitative and quantitative agreement with electrophoretic experiments done under conditions of low and high electric fields.

Prediction of Infant Drug Exposure Through Breastfeeding: Population PK Modeling and Simulation of Fluoxetine Exposure.

The likelihood of significant exposure to drugs in infants through breast milk is poorly defined, given the difficulties of conducting pharmacokinetics (PK) studies. Using fluoxetine (FX) as an example, we conducted a proof-of-principle study applying population PK (popPK) modeling and simulation to estimate drug exposure in infants through breast milk. We simulated data for 1,000 mother-infant pairs, assuming conservatively that the FX clearance in an infant is 20% of the allometrically adjusted value in adults. The model-generated estimate of the milk-to-plasma ratio for FX (mean: 0.59) was consistent with those reported in other studies. The median infant-to-mother ratio of FX steady-state plasma concentrations predicted by the simulation was 8.5%. Although the disposition of the active metabolite, norfluoxetine, could not be modeled, popPK-informed simulation may be valid for other drugs, particularly those without active metabolites, thereby providing a practical alternative to conventional PK studies for exposure risk assessment in this population.

Shaping fission yeast cells by rerouting actin-based transport on microtubules.

Kinesins and myosins transport cargos to specific locations along microtubules and actin filaments, respectively. The relative contribution of the two transport systems for cell polarization varies extensively in different cell types, with some cells relying exclusively on actin-based transport while others mainly use microtubules. Using fission yeast, we asked whether one transport system can substitute for the other. In this organism, microtubules and actin cables both contribute to polarized growth by transporting cargos to cell poles, but with distinct roles: microtubules transport landmarks to label cell poles for growth and actin assembly but do not directly contribute to the growth process [1]. Actin cables serve as tracks for myosin V delivery of growth vesicles to cell poles [2-4]. We engineered a chimera between the motor domain of the kinesin 7 Tea2 and the globular tail of the myosin V Myo52, which we show transports Ypt3, a myosin cargo receptor, to cell poles along microtubules. Remarkably, this chimera restores polarized growth and viability to cells lacking actin cables. It also bypasses the normal microtubule-dependent marking of cell poles for polarized growth, but not for other functions. Thus, a synthetic motor protein successfully redirects cargos along a distinct cytoskeletal route.

Molecular players shaping the arbuscular mycorrhizal symbiosis of rice

SUMMARY : The arbuscular mycorrhizal (AM) symbiosis is an evolutionarily ancient association between most land plants and Glomeromycotan fungi that is based on the mutual exchange of nutrients between the two partners. Its structural and physiological establishment is a multi-step process involving a tightly regulated signal exchange leading to intracellular colonization of roots by the fungi. Most research on the molecular biology and genetics of symbiosis development has been performed in dicotyledonous model legumes. In these, a plant signaling pathway, the common SYM pathway, has been found to be required for accommodation of both root symbionts rhizobia and AM fungi. Rice, a monocotyledon model and the world's most important staple crop also forms AM symbioses, has been largely ignored for studies of the AM symbiosis. Therefore in this PhD work functional conservation of the common SYM pathway in rice was addressed and demonstrated. Mycorrhiza-specific marker genes were established that are expressed at different stages of AM development and therefore represent readouts for various AM-specific signaling events. These tools were successfully used to obtain evidence for a yet unknown signaling network comprising common SYM-dependent and -independent events. In legumes AM colonization induces common SYM signaling dependent changes in root system architecture. It was demonstrated that also in rice, root system architecture changes in response to AM colonization but these alterations occur independently of common SYM signaling. The rice root system is complex and contains three different root types. It was shown that root type identity influences the quantity of AM colonization, indicating root type specific symbiotic properties. Interestingly, the root types differed in their transcriptional responses to AM colonization and the less colonized root type responded more dramatically than the more strongly colonized root type. Finally, in an independent project a novel mutant, inhospitable (iho), was discovered. It is perturbed at the most early step of AM colonization, namely differentiation of the AM fungal hyphae into a hyphopodium at the root surface. As plant factors required for this early step are not known, identification of the IHO gene will greatly contribute to the advance of mycorrhiza RÉSUMÉ : La symbiose mycorhizienne arbusculaire (AM) est une association évolutionnairement ancienne entre la majorité des plantes terrestres et les champignons du type Glomeromycota, basée sur l'échange mutuel d'éléments nutritifs entre les deux partenaires. Son établissement structural et physiologique est un processus en plusieurs étapes, impliquant des échanges de signaux étroitement contrôlés, aboutissant à la colonisation intracellulaire des racines par le champignon. La plupart des recherches sur la biologie moléculaire et la génétique du développement de la symbiose ont été effectuées sur des légumineuses dicotylédones modèles. Dans ces dernières, une voie de signalisation, la voie SYM, s'est avérée nécessaire pour permettre la mise en place de la symbiose mycorhizienne. Chez les plantes monocotylédones, comme le riz, une des céréales les plus importantes, nourrissant la moitié de la population mondiale, peu de recherches ont été effectuées sur les bases de la cette symbiose. Dans ce travail de thèse, la conservation fonctionnelle de la voie commune SYM chez le riz a été étudiée et démontrée. De plus, des gènes marqueurs spécifiques des différentes étapes du développement de l'AM ont été identifiés, permettant ainsi d'avoir des traceurs de la colonisation. Ces outils ont été utilisés avec succès pour démontrer l'existence d'un nouveau réseau de signalisation, comprenant des éléments SYM dépendant et indépendant. Chez les légumineuses, la colonisation par les AM induit des changements dans l'architecture du système racinaire, via la signalisation SYM dépendantes. Cependant chez le riz, il a été démontré que l'architecture de système racinaire changeait suite à la colonisation de l'AM, mais ceux, de façon SYM indépendante. Le système racinaire du riz est complexe et contient trois types différents de racines. Il a été démontré que le type de racine pouvait influencer l'efficacité de la colonisation par l'AM, indiquant que les racines ont des propriétés symbiotiques spécifiques différentes. De façon surprenante, les divers types de racines répondent de différemment suite à colonisation par l'AM avec des changements de la expression des gènes. Le type de racine le moins colonisé, répondant le plus fortement a la colonisation, et inversement. En parallèle, dans un projet indépendant, un nouveau mutant, inhospitable (iho), a été identifié. Ce mutant est perturbé lors de l'étape la plus précoce de la colonisation par l'AM, à savoir la différentiation des hyphes fongiques de l'AM en hyphopodium, à la surface des racines. Les facteurs d'origine végétale requis pour cette étape étant encore inconnus, l'identification du gène IHO contribuera considérablement a accroître nos connaissance sur les bases de la mise en place de cette symbiose.

Accurate and efficient simulation of multiphase flow in a heterogeneous reservoir by using error estimate and control in the multiscale finite-volume framework

The multiscale finite-volume (MSFV) method is designed to reduce the computational cost of elliptic and parabolic problems with highly heterogeneous anisotropic coefficients. The reduction is achieved by splitting the original global problem into a set of local problems (with approximate local boundary conditions) coupled by a coarse global problem. It has been shown recently that the numerical errors in MSFV results can be reduced systematically with an iterative procedure that provides a conservative velocity field after any iteration step. The iterative MSFV (i-MSFV) method can be obtained with an improved (smoothed) multiscale solution to enhance the localization conditions, with a Krylov subspace method [e.g., the generalized-minimal-residual (GMRES) algorithm] preconditioned by the MSFV system, or with a combination of both. In a multiphase-flow system, a balance between accuracy and computational efficiency should be achieved by finding a minimum number of i-MSFV iterations (on pressure), which is necessary to achieve the desired accuracy in the saturation solution. In this work, we extend the i-MSFV method to sequential implicit simulation of time-dependent problems. To control the error of the coupled saturation/pressure system, we analyze the transport error caused by an approximate velocity field. We then propose an error-control strategy on the basis of the residual of the pressure equation. At the beginning of simulation, the pressure solution is iterated until a specified accuracy is achieved. To minimize the number of iterations in a multiphase-flow problem, the solution at the previous timestep is used to improve the localization assumption at the current timestep. Additional iterations are used only when the residual becomes larger than a specified threshold value. Numerical results show that only a few iterations on average are necessary to improve the MSFV results significantly, even for very challenging problems. Therefore, the proposed adaptive strategy yields efficient and accurate simulation of multiphase flow in heterogeneous porous media.

Monte Carlo simulation of a whole-body counter using IGOR phantoms.

Whole-body counting is a technique of choice for assessing the intake of gamma-emitting radionuclides. An appropriate calibration is necessary, which is done either by experimental measurement or by Monte Carlo (MC) calculation. The aim of this work was to validate a MC model for calibrating whole-body counters (WBCs) by comparing the results of computations with measurements performed on an anthropomorphic phantom and to investigate the effect of a change in phantom's position on the WBC counting sensitivity. GEANT MC code was used for the calculations, and an IGOR phantom loaded with several types of radionuclides was used for the experimental measurements. The results show a reasonable agreement between measurements and MC computation. A 1-cm error in phantom positioning changes the activity estimation by >2%. Considering that a 5-cm deviation of the positioning of the phantom may occur in a realistic counting scenario, this implies that the uncertainty of the activity measured by a WBC is ∼10-20%.

Simulation of interface dynamics across pore discontinuities

We study the dynamics of a water-oil meniscus moving from a smaller to a larger pore. The process is characterised by an abrupt change in the configuration, yielding a sudden energy release. A theoretic study for static conditions provides analytical solutions of the surface energy content of the system. Although the configuration after the sudden energy release is energetically more convenient, an energy barrier must be overcome before the process can happen spontaneously. The energy barrier depends on the system geometry and on the flow parameters. The analytical results are compared to numerical simulations that solve the full Navier-Stokes equation in the pore space and employ the Volume Of Fluid (VOF) method to track the evolution of the interface. First, the numerical simulations of a quasi-static process are validated by comparison with the analytical solutions for a static meniscus, then numerical simulations with varying injection velocity are used to investigate dynamic effects on the configuration change. During the sudden energy jump the system exhibits an oscillatory behaviour. Extension to more complex geometries might elucidate the mechanisms leading to a dynamic capillary pressure and to bifurcations in final distributions of fluid phases in porous

Statistical Modeling of the Eye for Multimodal Treatment Planning for External Beam Radiation Therapy of Intraocular Tumors.

PURPOSE: Ocular anatomy and radiation-associated toxicities provide unique challenges for external beam radiation therapy. For treatment planning, precise modeling of organs at risk and tumor volume are crucial. Development of a precise eye model and automatic adaptation of this model to patients' anatomy remain problematic because of organ shape variability. This work introduces the application of a 3-dimensional (3D) statistical shape model as a novel method for precise eye modeling for external beam radiation therapy of intraocular tumors. METHODS AND MATERIALS: Manual and automatic segmentations were compared for 17 patients, based on head computed tomography (CT) volume scans. A 3D statistical shape model of the cornea, lens, and sclera as well as of the optic disc position was developed. Furthermore, an active shape model was built to enable automatic fitting of the eye model to CT slice stacks. Cross-validation was performed based on leave-one-out tests for all training shapes by measuring dice coefficients and mean segmentation errors between automatic segmentation and manual segmentation by an expert. RESULTS: Cross-validation revealed a dice similarity of 95% ± 2% for the sclera and cornea and 91% ± 2% for the lens. Overall, mean segmentation error was found to be 0.3 ± 0.1 mm. Average segmentation time was 14 ± 2 s on a standard personal computer. CONCLUSIONS: Our results show that the solution presented outperforms state-of-the-art methods in terms of accuracy, reliability, and robustness. Moreover, the eye model shape as well as its variability is learned from a training set rather than by making shape assumptions (eg, as with the spherical or elliptical model). Therefore, the model appears to be capable of modeling nonspherically and nonelliptically shaped eyes.

Integration of local-scale hydrological and regional-scale geophysical based on a nonlinear Bayesian sequential simulation approach

Significant progress has been made with regard to the quantitative integration of geophysical and hydrological data at the local scale. However, extending the corresponding approaches to the regional scale represents a major, and as-of-yet largely unresolved, challenge. To address this problem, we have developed a downscaling procedure based on a non-linear Bayesian sequential simulation approach. The basic objective of this algorithm is to estimate the value of the sparsely sampled hydraulic conductivity at non-sampled locations based on its relation to the electrical conductivity, which is available throughout the model space. The in situ relationship between the hydraulic and electrical conductivities is described through a non-parametric multivariate kernel density function. This method is then applied to the stochastic integration of low-resolution, re- gional-scale electrical resistivity tomography (ERT) data in combination with high-resolution, local-scale downhole measurements of the hydraulic and electrical conductivities. Finally, the overall viability of this downscaling approach is tested and verified by performing and comparing flow and transport simulation through the original and the downscaled hydraulic conductivity fields. Our results indicate that the proposed procedure does indeed allow for obtaining remarkably faithful estimates of the regional-scale hydraulic conductivity structure and correspondingly reliable predictions of the transport characteristics over relatively long distances.

Whole anterior segment proton beam radiotherapy for diffuse iris melanoma.

AIM: To report the results of whole anterior segment proton beam radiotherapy (PBR) for diffuse iris melanoma. METHODS: Between 2000 and 2011, 12 patients with iris melanoma received PBR to the entire iris and ciliary body. RESULTS: Patients had a mean age of 57 years and a median follow-up of 3.5 years (range 1-11.6 years). Tumour iris involvement was 1-4 h in five patients, 5-8 h in four and 9-12 h in three. Angle involvement was 6-8 h in five patients and 9-12 h in seven. The visual acuity (VA) before treatment was 6/5-6/6 in six patients, 6/8-6/9 in three and 6/18-6/38 in three. No tumour recurrence occurred during the follow-up period. Glaucoma treatment was required in 11 of 12 patients. The visual acuity at the last follow-up was 6/5-6/9 in five patients, 6/18-6/24 in three, 6/60-1/60 in two and no light perception in two. Four patients developed varying non-severe degrees of limbal stem cell deficiency, which was treatable with conservative measures. CONCLUSIONS: Whole anterior segment PBR is a useful alternative to enucleation for diffuse iris melanoma. Most patients will need treatment for glaucoma and some may require treatment for tear-film instability and/or stem cell failure.

Influence of scatter reduction method and monochromatic beams on image quality and dose in mammography.

In mammography, the image contrast and dose delivered to the patient are determined by the x-ray spectrum and the scatter to primary ratio S/P. Thus the quality of the mammographic procedure is highly dependent on the choice of anode and filter material and on the method used to reduce the amount of scattered radiation reaching the detector. Synchrotron radiation is a useful tool to study the effect of beam energy on the optimization of the mammographic process because it delivers a high flux of monochromatic photons. Moreover, because the beam is naturally flat collimated in one direction, a slot can be used instead of a grid for scatter reduction. We have measured the ratio S/P and the transmission factors for grids and slots for monoenergetic synchrotron radiation. In this way the effect of beam energy and scatter rejection method were separated, and their respective importance for image quality and dose analyzed. Our results show that conventional mammographic spectra are not far from optimum and that the use of a slot instead of a grid has an important effect on the optimization of the mammographic process. We propose a simple numerical model to quantify this effect.

Dosage optimization of treatments using population pharmacokinetic modeling and simulation.

Pharmacokinetic variability in drug levels represent for some drugs a major determinant of treatment success, since sub-therapeutic concentrations might lead to toxic reactions, treatment discontinuation or inefficacy. This is true for most antiretroviral drugs, which exhibit high inter-patient variability in their pharmacokinetics that has been partially explained by some genetic and non-genetic factors. The population pharmacokinetic approach represents a very useful tool for the description of the dose-concentration relationship, the quantification of variability in the target population of patients and the identification of influencing factors. It can thus be used to make predictions and dosage adjustment optimization based on Bayesian therapeutic drug monitoring (TDM). This approach has been used to characterize the pharmacokinetics of nevirapine (NVP) in 137 HIV-positive patients followed within the frame of a TDM program. Among tested covariates, body weight, co-administration of a cytochrome (CYP) 3A4 inducer or boosted atazanavir as well as elevated aspartate transaminases showed an effect on NVP elimination. In addition, genetic polymorphism in the CYP2B6 was associated with reduced NVP clearance. Altogether, these factors could explain 26% in NVP variability. Model-based simulations were used to compare the adequacy of different dosage regimens in relation to the therapeutic target associated with treatment efficacy. In conclusion, the population approach is very useful to characterize the pharmacokinetic profile of drugs in a population of interest. The quantification and the identification of the sources of variability is a rational approach to making optimal dosage decision for certain drugs administered chronically.