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.

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.

Late Proterozoic thrust tectonics, high-pressure metamorphism and Uranium mineralization in the domes area, Lufilian arc, northwestern Zambia

The following main lithostratigraphic units have been distinguished in the Domes Area. The Kibaran basement complex composed of gneisses, migmatites with amphibolite bands and metagranites is exposed in dome structures; metamorphic features of Kibaran age have been almost completely obliterated by extensive Lufilian reactivation. The post-Kibaran cover sequence is subdivided into the Lower Roan Group consisting of well-preserved quartzites with high Mg content, talc-bearing, extremely foliated schists intercalated with pseudo-conglomerates of tectonic origin and the Upper Roan Group including dolomitic marbles with rare stromatolites, metapelites and a sequence of detrital metasediments, with local volcano-sedimentary components and interlayered banded ironstones. The sediments of the Lower Roan Group are interpreted as continental to lagoonal-evaporitic deposits partly converted into the talc-kyanite + garnet assemblage characteristic of ``white schists''. The dolomites and metapelites of the Upper Roan Group are attributed to a carbonate platform sequence progressively subsiding under terrigenous deposits, whilst the detrital metasediments and BIF may be interpreted as a basinal sequence, probably deposited on oceanic crust grading laterally into marbles. Metagabbros and metabasalts are considered as remnants of an ocean-floor-type crustal unit probably related to small basins. Alkaline stocks of Silurian age intruded the post-Kibaran cover. Significant ancestral tectonic discontinuities promoted the development of a nappe pile that underwent high-pressure metamorphism during the Lufilian orogeny and all lithostratigraphic units. Rb-Sr and K-Ar and U-Pb data indicate an age of 700 Ma for the highest grade metamorphism and 500 Ma for blocking of the K-Ar and Rb-Sr system in micas, corresponding to the time when the temperature dropped below 350-degrees-400-degrees-C and to an age of about 400 Ma for the emplacement of hypabyssal syenitic bodies. A first phase of crustal shortening by decoupling of basement and cover slices along shallow shear zones has been recognized. Fluid-rich tectonic slabs of cover sediments were thus able to transport fluids into the anhydrous metamorphic basement or mafic units. During the subsequent metamorphic re-equilibration stage of high pressure, pre-existing thrusts horizons were converted into recrystallized mylonites. Due to uplift, rocks were re-equilibrated into assemblages compatible with lower pressures and slightly lower temperatures. This stage occurs under a decompressional (nearly adiabatic) regime, with P(fluid) almost-equal-to P(lithostatic). It is accompanied by metasomatic development of minerals, activated by injection of hot fluids. New or reactivated shear zones and mylonitic belts were the preferred conduits of fluids. The most evident regional-scale effect of these processes is the intense metasomatic scapolitization of formerly plagioclase-rich lithologies. Uraninite mineralization can probably be assigned to the beginning of the decompressional stage. A third regional deformation phase characterized by open folds and local foliation is not accompanied by significant growth of new minerals. However, pitchblende mineralization can be ascribed to this phase as late-stage, short-range remobilization of previously existing deposits. Finally, shallow alkaline massifs were emplaced when the level of the Domes Area now exposed was already subjected to exchange with meteoric circuits, activated by residual geothermal gradients generally related to intrusions or rifting. Most of the superficial U-showings with U-oxidation products were probably generated during this relatively recent phase.

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.

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.

The yield of high-density electric source imaging in focal epilepsy

Background: EEG is the cornerstone of epilepsy diagnostics and mandatory to determine the underlying epilepsy syndrome (e.g. focal vs idiopathic generalized). However, its potential as imaging tool is still underrecognized. In the present study, we aim to determine the prerequisites of maximal benefit of electric source imaging (ESI) to localize the irritative zone in patients with focal epilepsy. Methods: 150 patients suffering from focal epilepsy and with minimum 1 year post-operative follow-up were studied prospectively by reviewers blinded to the underlying diagnosis and outcome. We evaluated the influence of two important factors on sensitivity and specificity of ESI: the number of electrodes (low resolution, LR-ESI: \30 vs. high resolution, HR-ESI: 128-256 electrodes), and the use of individual MRI (i-MRI) vs. template MRI (t-MRI) as head model.Results: ESI had a sensitivity of 85% and a specificity of 87% when HR-ESI with i-MRI was used. Using LR-ESI, sensitivity decreased to 68%, or even 57% when only t-MRI was available. The sensitivity of HR-ESI/i-MRI compared favorably with those of MRI (76%), PET (69%) and ictal/interictal SPECT (64%).Interpretation: This study on a large patient group shows excellent sensitivity and specificity of ESI if 128 EEG channels or more are used for ESI and if the results are co-registered to the patient's individual MRI. Localization precision is as high as or even higher than established brain imaging techniques, providing excellent costeffectiveness in epilepsy evaluation. HR-ESI appears to be a valuable additional imaging tool, given that larger electrode arrays are easily and rapidly applied with modern EEG equipment and that structural MRI is nearly always available for these patients.

The contribution of the young Cretaceous Caribbean Oceanic Plateau to the genesis of late Cretaceous arc magmatism in the Cordillera Occidental of Ecuador

The eastern part of the Cordillera Occidental of Ecuador comprises thick buoyant oceanic plateaus associated with island-arc tholeiites and subduction-related calc-alkaline series, accreted to the Ecuadorian Continental Margin from Late Cretaceous to Eocene times. One of these plateau sequences, the Guaranda Oceanic Plateau is considered as remnant of the Caribbean-Colombian Oceanic Province (CCOP) accreted to the Ecuadorian Margin in the Maastrichtien. Samples studied in this paper were taken from four cross-sections through two arc-sequences in the northern part of the Cordillera Occidental of Ecuador, dated as (Rio Cala) or ascribed to (Macuchi) the Late Cretaceous and one arc-like sequence in the Chogon-Colonche Cordillera (Las Orquideas). These three island-arcs can clearly be identified and rest conformably on the CCOP. In all four localities, basalts with abundant large clinopyroxene phenocrysts can be found, mimicking a picritic or ankaramitic facies. This mineralogical particularity, although not uncommon in island arc lavas, hints at a contribution of the CCOP in the genesis of these island arc rocks. The complete petrological and geochemical study of these rocks reveals that some have a primitive island-arc nature (MgO values range from 6 to 11 wt.%). Studied samples display marked Nb, Ta and Ti negative anomalies relative to the adjacent elements in the spidergrams characteristic of subductionrelated magmatism. These rocks are LREE-enriched and their clinopyroxenes show a tholeiitic affinity (FeO(1)-TiO(2) enrichment and CaO depletion from core to rim within a single crystal). The four sampled cross-sections through the island-arc sequences display homogeneous initial Nd, and Pb isotope ratios that suggest a unique mantellic source for these rocks resulting from the mixing of three components: an East-Pacific MORB end-member, an enriched pelagic sediment component, and a HIMU component carried by the CCOP. Indeed, the ankaramite and Mg-basalt sequences that form part of the Caribbean-Colombian Oceanic Plateau are radiogenically enriched in (206)Pb/(204)Pb and (207)Pb/(204)Pb and contain a HIMU component similar to that observed in the Gorgona basalts and Galapagos lavas. The subduction zone that generated the Late Cretaceous arcs occurred far from the continental margin, in an oceanic environment. This implies that no terrigenous detrital sediments interacted with the source at this period. Thus, the enriched component can only result from the melting of subducted pelagic sediments. We have thus defined the East-Pacific MORB, enriched (cherts, pelagic sediments) and HIMU components in an attempt to constrain and model the genesis of the studied island-arc magmatism, using a compilation of carefully selected isotopic data from literature according to rock age and paleogeographic location at the time of arc edification. Tripolar mixing models reveal that proportions of 12-15 wt.% of the HIMU component, 7-15 wt.% of the pelagic sediment end-member and 70-75 wt.% of an East-pacific MORB end-member are needed to explain the measured isotope ratios. These surprisingly high proportions of the HIMU/CCOP component could be explained by the young age of the oceanic plateau (5-15 Ma) during the Late Cretaceous arc emplacement. The CCOP, basement of these arc sequences, was probably still hot and easily assimilated at the island-arc lava source. (C) 2008 Elsevier Ltd. All rights reserved,

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.

Numerical simulation of ambient seismic wavefield modification caused by pore-fluid effects in an oil reservoir

We have modeled numerically the seismic response of a poroelastic inclusion with properties applicable to an oil reservoir that interacts with an ambient wavefield. The model includes wave-induced fluid flow caused by pressure differences between mesoscopic-scale (i.e., in the order of centimeters to meters) heterogeneities. We used a viscoelastic approximation on the macroscopic scale to implement the attenuation and dispersion resulting from this mesoscopic-scale theory in numerical simulations of wave propagation on the kilometer scale. This upscaling method includes finite-element modeling of wave-induced fluid flow to determine effective seismic properties of the poroelastic media, such as attenuation of P- and S-waves. The fitted, equivalent, viscoelastic behavior is implemented in finite-difference wave propagation simulations. With this two-stage process, we model numerically the quasi-poroelastic wave-propagation on the kilometer scale and study the impact of fluid properties and fluid saturation on the modeled seismic amplitudes. In particular, we addressed the question of whether poroelastic effects within an oil reservoir may be a plausible explanation for low-frequency ambient wavefield modifications observed at oil fields in recent years. Our results indicate that ambient wavefield modification is expected to occur for oil reservoirs exhibiting high attenuation. Whether or not such modifications can be detected in surface recordings, however, will depend on acquisition design and noise mitigation processing as well as site-specific conditions, such as the geologic complexity of the subsurface, the nature of the ambient wavefield, and the amount of surface noise.