Orthotropic active strain models for the numerical simulation of cardiac biomechanics

An active strain formulation for orthotropic constitutive laws arising in cardiac mechanics modeling is introduced and studied. The passive mechanical properties of the tissue are described by the Holzapfel-Ogden relation. In the active strain formulation, the Euler-Lagrange equations for minimizing the total energy are written in terms of active and passive deformation factors, where the active part is assumed to depend, at the cell level, on the electrodynamics and on the specific orientation of the cardiac cells. The well-posedness of the linear system derived from a generic Newton iteration of the original problem is analyzed and different mechanical activation functions are considered. In addition, the active strain formulation is compared with the classical active stress formulation from both numerical and modeling perspectives. Taylor-Hood and MINI finite elements are employed to discretize the mechanical problem. The results of several numerical experiments show that the proposed formulation is mathematically consistent and is able to represent the main key features of the phenomenon, while allowing savings in computational costs.

Population PK-guided simulation study in children exposed to drugs in human milk

Risks of significant infant drug exposure through human milk arepoorly defined due to lack of large-scale PK data. We propose to useBayesian approach based on population PK (popPK)-guided modelingand simulation for risk prediction. As a proof-of-principle study, weexploited fluoxetine milk concentration data from 25 women. popPKparameters including milk-to-plasma ratio (MP ratio) were estimatedfrom the best model. The dose of fluoxetine the breastfed infant wouldreceive through mother's milk, and infant plasma concentrations wereestimated from 1000 simulated mother-infant pairs, using randomassignment of feeding times and milk volume. A conservative estimateof CYP2D6 activity of 20% of the allometrically-adjusted adult valuewas assumed. Derived model parameters, including MP ratio were consistentwith those reported in the literature. Visual predictive check andother model diagnostics showed no signs of model misspecifications.The model simulation predicted that infant exposure levels to fluoxetinevia mother's milk were below 10% of weight-adjusted maternal therapeuticdoses in >99% of simulated infants. Predicted median ratio ofinfant-mother serum levels at steady state was 0.093 (range 0.033-0.31),consistent with literature reported values (mean=0.07; range 0-0.59).Predicted incidence of relatively high infant-mother ratio (>0.2) ofsteady-state serum fluoxetine concentrations was <1.3%. Overall, ourpredictions are consistent with clinical observations. Our approach maybe valid for other drugs, allowing in silico prediction of infant drugexposure risks through human milk. We will discuss application of thisapproach to another drug used in lactating women.

Prediction of infant drug exposure through breastfeeding: population PK modeling and simulation of fluoxetine

BACKGROUND: Risks of significant infant drug exposurethrough breastmilk are poorly defined for many drugs, and largescalepopulation data are lacking. We used population pharmacokinetics(PK) modeling to predict fluoxetine exposure levels ofinfants via mother's milk in a simulated population of 1000 motherinfantpairs.METHODS: Using our original data on fluoxetine PK of 25breastfeeding women, a population PK model was developed withNONMEM and parameters, including milk concentrations, wereestimated. An exponential distribution model was used to account forindividual variation. Simulation random and distribution-constrainedassignment of doses, dosing time, feeding intervals and milk volumewas conducted to generate 1000 mother-infant pairs with characteristicssuch as the steady-state serum concentrations (Css) and infantdose relative to the maternal weight-adjusted dose (relative infantdose: RID). Full bioavailability and a conservative point estimate of1-month-old infant CYP2D6 activity to be 20% of the adult value(adjusted by weigth) according to a recent study, were assumed forinfant Css calculations.RESULTS: A linear 2-compartment model was selected as thebest model. Derived parameters, including milk-to-plasma ratios(mean: 0.66; SD: 0.34; range, 0 - 1.1) were consistent with the valuesreported in the literature. The estimated RID was below 10% in >95%of infants. The model predicted median infant-mother Css ratio was0.096 (range 0.035 - 0.25); literature reported mean was 0.07 (range0-0.59). Moreover, the predicted incidence of infant-mother Css ratioof >0.2 was less than 1%.CONCLUSION: Our in silico model prediction is consistent withclinical observations, suggesting that substantial systemic fluoxetineexposure in infants through human milk is rare, but further analysisshould include active metabolites. Our approach may be valid forother drugs. [supported by CIHR and Swiss National Science Foundation(SNSF)]

SIMULIT : un modèle de simulation pour l'analyse et la planification de l'activité hospitalière

Le modèle développé à l'Institut universitaire de médecine sociale et préventive de Lausanne utilise un programme informatique pour simuler les mouvements d'entrées et de sorties des hôpitaux de soins généraux. Cette simulation se fonde sur les données récoltées de routine dans les hôpitaux; elle tient notamment compte de certaines variations journalières et saisonnières, du nombre d'entrées, ainsi que du "Case-Mix" de l'hôpital, c'est-à-dire de la répartition des cas selon les groupes cliniques et l'âge des patients.

Intravenous streptomycin dosing regimen in a patient undergoing hemodialysis: Plasma level monitoring and pharmacokinetic simulation

Introduction: Streptomycin, as other aminoglycosides, exhibits concentration-dependent bacterial killing but has a narrow therapeutic window. It is primarily eliminated unchanged by the kidneys. Data and dosing information to achieve a safe regimen in patients with chronic renal failure undergoing hemodialysis (HD) are scarce. Although main adverse reactions are related to prolonged, elevated serum concentrations, literature recommendation is to administer streptomycin after each HD. Patients (or Materials) and Methods: We report the case of a patient with end-stage renal failure, undergoing HD, who was successfully treated with streptomycin for gentamicin-resistant Enterococcus faecalis bacteremia with prosthetic arteriovenous fistula infection. Streptomycin was administered intravenously 7.5 mg/kg, 3 hours before each dialysis (3 times a week) during 6 weeks in combination with amoxicillin. Streptomycin plasma levels were monitored with repeated blood sampling before, after, and between HD sessions. A 2-compartment model was used to reconstruct the concentration time profile over days on and off HD. Results: Streptomycin trough plasma-concentration was 2.8 mg/L. It peaked to 21.4 mg/L 30 minutes after intravenous administration, decreased to 18.2 mg/L immediately before HD, and dropped to 4.5 mg/L at the end of a 4-hour HD session. Plasma level increased again to 5.7 mg/L 2 hours after the end of HD and was 2.8 mg/L 48 hours later, before the next administration and HD. The pharmacokinetics of streptomycin was best described with a 2-compartment model. The computer simulation fitted fairly well to the observed concentrations during or between HD sessions. Redistribution between the 2 compartments after the end of HD reproduced the rebound of plasma concentrations after HD. No significant toxicity was observed during treatment. The outcome of the infection was favorable, and no sign of relapse was observed after a follow-up of 3 months. Conclusion: Streptomycin administration of 7.5 mg/kg 3 hours before HD sessions in a patient with end-stage renal failure resulted in an effective and safe dosing regimen. Monitoring plasma levels along with pharmacokinetic simulation document the suitability of this dosing scheme, which should replace current dosage recommendations for streptomycin in HD.

Under- and over-water halves of Gyrinidae beetle eyes harbor different corneal nanocoatings providing adaptation to the water and air environments.

Whirligig beetles (Gyrinidae) inhabit water surfaces and possess unique eyes which are split into the overwater and underwater parts. In this study we analyze the micro- and nanostructure of the split eyes of two Gyrinidae beetles genera, Gyrinus and Orectochilus. We find that corneae of the overwater ommatidia are covered with maze-like nanostructures, while the corneal surface of the underwater eyes is smooth. We further show that the overwater nanostructures possess no anti-wetting, but the anti-reflective properties with the spectral preference in the range of 450-600 nm. These findings illustrate the adaptation of the corneal nanocoating of the two halves of an insect's eye to two different environments. The novel natural anti-reflective nanocoating we describe may find future technological applications.

A pseudo-spectral method for the simulation of poro-elastic seismic wave propagation in 2D polar coordinates using domain decomposition

We present a novel numerical approach for the comprehensive, flexible, and accurate simulation of poro-elastic wave propagation in 2D polar coordinates. An important application of this method and its extensions will be the modeling of complex seismic wave phenomena in fluid-filled boreholes, which represents a major, and as of yet largely unresolved, computational problem in exploration geophysics. In view of this, we consider a numerical mesh, which can be arbitrarily heterogeneous, consisting of two or more concentric rings representing the fluid in the center and the surrounding porous medium. The spatial discretization is based on a Chebyshev expansion in the radial direction and a Fourier expansion in the azimuthal direction and a Runge-Kutta integration scheme for the time evolution. A domain decomposition method is used to match the fluid-solid boundary conditions based on the method of characteristics. This multi-domain approach allows for significant reductions of the number of grid points in the azimuthal direction for the inner grid domain and thus for corresponding increases of the time step and enhancements of computational efficiency. The viability and accuracy of the proposed method has been rigorously tested and verified through comparisons with analytical solutions as well as with the results obtained with a corresponding, previously published, and independently bench-marked solution for 2D Cartesian coordinates. Finally, the proposed numerical solution also satisfies the reciprocity theorem, which indicates that the inherent singularity associated with the origin of the polar coordinate system is adequately handled.

Genome-Wide Prediction Methods in Highly Diverse and Heterozygous Species: Proof-of-Concept through Simulation in Grapevine.

Nowadays, genome-wide association studies (GWAS) and genomic selection (GS) methods which use genome-wide marker data for phenotype prediction are of much potential interest in plant breeding. However, to our knowledge, no studies have been performed yet on the predictive ability of these methods for structured traits when using training populations with high levels of genetic diversity. Such an example of a highly heterozygous, perennial species is grapevine. The present study compares the accuracy of models based on GWAS or GS alone, or in combination, for predicting simple or complex traits, linked or not with population structure. In order to explore the relevance of these methods in this context, we performed simulations using approx 90,000 SNPs on a population of 3,000 individuals structured into three groups and corresponding to published diversity grapevine data. To estimate the parameters of the prediction models, we defined four training populations of 1,000 individuals, corresponding to these three groups and a core collection. Finally, to estimate the accuracy of the models, we also simulated four breeding populations of 200 individuals. Although prediction accuracy was low when breeding populations were too distant from the training populations, high accuracy levels were obtained using the sole core-collection as training population. The highest prediction accuracy was obtained (up to 0.9) using the combined GWAS-GS model. We thus recommend using the combined prediction model and a core-collection as training population for grapevine breeding or for other important economic crops with the same characteristics.

Tax simulation in the SHP

This paper describes a simulation package designed to estimate the annual income taxes paid by respondents of the Swiss Household Panel (SHP). In Switzerland, the 26 cantons have their own tax system. Additionally, tax levels vary between the over 2000 municipalities and over time. The simulation package takes account of this complexity by building on existing tables on tax levels which are provided by the Swiss Federal Tax Administration Office. Because these are limited to a few types of households and only 812 municipalities, they have to be extended to cover all households and municipalities. A further drawback of these tables is that they neglect several deductions. The tax simulation package fills this gap by taking additionally account of deductions for children, double-earner couples, third pillar and support for dependent persons according to cantonal legislation. The resulting variable on direct taxes not only serves to calculate household income net of taxes, but can also be a variable for analysis by its own account.

Biosphere/lithosphere interface : microbially-mediated CaCO3 precipitation in hypersaline and freshwater environments

Abstract: Microbial mats very efficiently cycle elements, such as C, 0, N, S and H, which makes them key players of redox processes at the biosphere-lithosphere interface. They are characterized by high metabolic activities and high turnover rates (production and consumption) of biomass, which mainly consists of cell material and of extracellular organic matter (EOM). The EOM forms a matrix, embedding the microbial cells and fulfilling various functions within the microbial mat, including: mat attachment to surfaces; creation of micro-domains within the mat; physical stabilization under hy- drodynamic stress and the protection of the cells in multiple other stress conditions. EOM mainly consists of polysaccharides, amino acids, and a variety of chemical func-tional groups {e.g., -C00H, - SH -OH). These groups strongly bind cations such as Ca2+ and Mg2+ and thus exert a strong control on carbonate mineral formation within the microbial mat. A feedback mechanism between community metabolisms, their prod¬ucts, and the surrounding physicochemical microenvironment thus influences the de¬gree of carbonate saturation favoring either carbonate precipitation or dissolution. We investigated the driving forces and mechanisms of microbialite formation in the Sari ne River, FR, Switzerland, the hypersaline lake, Big Pond, Bahamas and in labo¬ratory experiments. The two fundamentally different natural systems allowed us to compare the geochemical conditions and microbial metabolisms, necessary for car¬bonate formation in microbial mats. Although carbonates are oversaturated in both environments, precipitation does not occur on physicochemical substrates (i.e. out¬side the microbial mats). In the Sarine a high crystal nucleation threshold exceeds the carbonate saturation, despite the high carbonate alkalinity in the water column. Cyanobacterial photosynthesis strongly locally enhances the carbonate alkalinity, whereas the EOM attract and immobilize calcium, which increases the saturation state and finally leads to carbonate precipitation within the EOM (in this case the cyanobacterial sheath) as nucleation template. In Big Pond, the presence of calcium- chelating anions (i.e. sulfate) and EOM, as well as the presence of magnesium, lowers the calcium activity in the water column and mat, and thus inhibits carbonate pre¬cipitation. Coupled with other heterotrophic metabolisms, sulfate reduction uses the EOM as carbon source, degrading it. The resulting EOM consumption creates alkalin¬ity, releases calcium and consumes sulfate in mat-micro domains, which leads to the formation of carbonate layers at the top of the microbial mat. Résumé: Interface biosphère/lithosphère: médiation microbienne de la précipitation de CaC03 dans des environnements en eaux douces et hypersalines Les tapis microbiens engendrent une circulation très efficace des éléments, tels que C, 0, N, S et H, ce qui en fait des acteurs clé pour les processus d'oxydoréduction à l'inter¬face biosphère-lithosphère. Ils sont caractérisés par des taux élevés d'activité méta¬bolique, ainsi que par la production et la consommation de biomasse, principalement constituée de cellules microbiennes et de matière organique extracellulaire (MOE). Dans un tapis microbien, les cellules microbiennes sont enveloppées par une matrice de MOE qui a différentes fonctions dont l'attachement du tapis aux surfaces, la créa¬tion de micro-domaines dans le tapis, la stabilisation physique en situation de stress hydrodynamique, et la protection des cellules dans de multiples autres conditions de stress. La MOE se compose principalement de polysaccharides, d'acides aminés, et d'une variété de groupes fonctionnels chimiques (par exemple, COOH, -SH et -OH). Ces groupes se lient fortement aux cations, tels que Ca2+ et Mg2+, et exercent ainsi un contrôle fort sur la formation de CaC03 dans le tapis microbien. Un mécanisme de rétroaction, entre les métabolismes de la communauté microbienne, leurs produits, et le microenvironnement physico-chimique, influence le degré de saturation de car¬bonate, favorisant soit leur précipitation, soit leur dissolution. Nous avons étudié le moteur et les mécanismes de minéralisation dans des tapis de la Sarine, FR, Suisse et du lac hypersalin, Big Pond, aux Bahamas, ainsi que durant des expériences en laboratoire. Les deux systèmes naturels, fondamentalement dif¬férents, nous ont permis de comparer les conditions géochimiques et les métabolis¬mes nécessaires à la formation des carbonates dans des tapis microbiens. Bien que les carbonates soient sursaturés dans les deux environnements, la précipitation ne se produit pas sur des substrats physico-chimiques (en dehors du tapis microbien). Dans la Sarine, malgré un taux d'alcalinité élevé, les valeurs de seuil pour la nucléa- tion de carbonates sont plus hautes que la saturation du carbonate. La photosynthèse cyanobactérienne augmente localement l'alcalinité, alors que la MOE attire et immo¬bilise le calcium, ce qui augmente l'état de saturation et conduit finalement à la pré¬cipitation des carbonates, en utilisant la MOE comme substrat de nucléation. À Big Pond, la présence de chélateurs de calcium, notamment les anions (p.ex. le sulfate) et la MOE, ainsi que la présence de magnésium, réduit l'activité du calcium et inhibe en conséquence la précipitation des carbonates. Couplée avec d'autres métabolismes hétérotrophes, la réduction des sulfates utilise la MOE comme source de carbone, en la dégradant. Cette consommation de MOE crée l'alcalinité, consomme des sulfates et libère du calcium dans des micro-domaines, conduisant à la formation de couches de carbonates dans le haut du tapis microbien.

Direct simulation of interface dynamics: linking capillary pressure, interfacial area and surface energy

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 to model the transition from stable flow to viscous fingering, we focus on the definition of macroscopic capillary pressure. When the fluids are at rest, the difference between inlet and outlet pressures and the difference between the intrinsic phase average pressure coincide with the capillary pressure. However, when the fluids are in motion these quantities are dominated by viscous forces. In this case, only a definition based on the variation of the interfacial energy provides an accurate measure of the macroscopic capillary pressure and allows separating the viscous from the capillary pressure components.

Brownian dynamics simulation of DNA condensation.

DNA condensation observed in vitro with the addition of polyvalent counterions is due to intermolecular attractive forces. We introduce a quantitative model of these forces in a Brownian dynamics simulation in addition to a standard mean-field Poisson-Boltzmann repulsion. The comparison of a theoretical value of the effective diameter calculated from the second virial coefficient in cylindrical geometry with some experimental results allows a quantitative evaluation of the one-parameter attractive potential. We show afterward that with a sufficient concentration of divalent salt (typically approximately 20 mM MgCl(2)), supercoiled DNA adopts a collapsed form where opposing segments of interwound regions present zones of lateral contact. However, under the same conditions the same plasmid without torsional stress does not collapse. The condensed molecules present coexisting open and collapsed plectonemic regions. Furthermore, simulations show that circular DNA in 50% methanol solutions with 20 mM MgCl(2) aggregates without the requirement of torsional energy. This confirms known experimental results. Finally, a simulated DNA molecule confined in a box of variable size also presents some local collapsed zones in 20 mM MgCl(2) above a critical concentration of the DNA. Conformational entropy reduction obtained either by supercoiling or by confinement seems thus to play a crucial role in all forms of condensation of DNA.