Copper and human health: biochemistry, genetics, and strategies for modeling dose-response relationships

Copper (Cu) and its alloys are used extensively in domestic and industrial applications. Cu is also an essential element in mammalian nutrition. Since both copper deficiency and copper excess produce adverse health effects, the dose-response curve is U-shaped, although the precise form has not yet been well characterized. Many animal and human studies were conducted on copper to provide a rich database from which data suitable for modeling the dose-response relationship for copper may be extracted. Possible dose-response modeling strategies are considered in this review, including those based on the benchmark dose and categorical regression. The usefulness of biologically based dose-response modeling techniques in understanding copper toxicity was difficult to assess at this time since the mechanisms underlying copper-induced toxicity have yet to be fully elucidated. A dose-response modeling strategy for copper toxicity was proposed associated with both deficiency and excess. This modeling strategy was applied to multiple studies of copper-induced toxicity, standardized with respect to severity of adverse health outcomes and selected on the basis of criteria reflecting the quality and relevance of individual studies. The use of a comprehensive database on copper-induced toxicity is essential for dose-response modeling since there is insufficient information in any single study to adequately characterize copper dose-response relationships. The dose-response modeling strategy envisioned here is designed to determine whether the existing toxicity data for copper excess or deficiency may be effectively utilized in defining the limits of the homeostatic range in humans and other species. By considering alternative techniques for determining a point of departure and low-dose extrapolation (including categorical regression, the benchmark dose, and identification of observed no-effect levels) this strategy will identify which techniques are most suitable for this purpose. This analysis also serves to identify areas in which additional data are needed to better define the characteristics of dose-response relationships for copper-induced toxicity in relation to excess or deficiency.

Combined bone and soft-tissue augmentation surgery in temporo-orbital contour reconstruction

Temporal hollowing due to temporal muscle atrophy after standard skull base surgery is common. Various techniques have been previously described to correct the disfiguring defect. Most often reconstruction is performed using freehand molded polymethylmethacrylate cement. This method and material are insufficient in terms of aesthetic results and implant characteristics. We herein propose reconstruction of such defects with a polyetheretherketone (PEEK)-based patient-specific implant (PSI) including soft-tissue augmentation to preserve normal facial topography. We describe a patient who presented with a large temporo-orbital hemangioma that had been repaired with polymethylmethacrylate 25 years earlier. Because of a toxic skin atrophy fistula, followed by infection and meningitis, this initial implant had to be removed. The large, disfiguring temporo-orbital defect was reconstructed with a PEEK-based PSI. The lateral orbital wall and the temporal muscle atrophy were augmented with computer-aided design and surface modeling techniques. The operative procedure to implant and adopt the reconstructed PEEK-based PSI was simple, and an excellent cosmetic outcome was achieved. The postoperative clinical course was uneventful over a 5-year follow-up period. Polyetheretherketone-based combined bony and soft contour remodeling is a feasible and effective method for cranioplasty including combined bone and soft-tissue reconstruction of temporo-orbital defects. Manual reconstruction of this cosmetically delicate area carries an exceptional risk of disfiguring results. Augmentation surgery in this anatomic location needs accurate PSIs to achieve satisfactory cosmetic results. The cosmetic outcome achieved in this case is superior compared with previously reported techniques.

A Communication Platform for Tele-monitoring and Tele-management of Type 1 Diabetes

Type 1 diabetes mellitus is a chronic disease characterized by blood glucose levels out of normal range due to inability of insulin production. This dysfunction leads to many short- and long-term complications. In this paper, a system for tele-monitoring and tele-management of Type 1 diabetes patients is proposed, aiming at reducing the risk of diabetes complications and improving quality of life. The system integrates Wireless Personal Area Networks (WPAN), mobile infrastructure, and Internet technology along with commercially available and novel glucose measurement devices, advanced modeling techniques, and tools for the intelligent processing of the available diabetes patients information. The integration of the above technologies enables intensive monitoring of blood glucose levels, treatment optimisation, continuous medical care, and improvement of quality of life for Type 1 diabetes patients, without restrictions in everyday life activities.

Evaluation of a new commercial enzyme-linked immunosorbent assay for the detection of porcine antibodies against Trichinella spp

Trichinellosis is a zoonotic disease that is caused by the nematode Trichinella spp. Both European Union regulations and guidelines from the World Organization for Animal Health foresee the possibility of conducting serological surveillance for Trichinella spp. A newly developed commercial enzyme-linked immunosorbent assay (ELISA) was evaluated against 2 existing diagnostic techniques: an in-house ELISA and an in-house Western blot. A total of 875 Trichinella larva-negative samples of pigs and 93 Trichinella larva-positive samples of both naturally and experimentally infected pigs were included in the study. Bayesian modeling techniques were used to correct for the absence of a perfect reference test. The sensitivity and specificity of the commercial ELISA was 97.1-97.8% and 99.5-99.8%, respectively. Sensitivity analysis demonstrated high stability in the models. In a serological surveillance system, ELISA-positive samples should be tested by a confirmatory test. The Western blot is a suitable test for this purpose. With the use of the results of the models, the sensitivity and specificity of a test protocol in both ELISA and Western blot were 95.9% and 99.9%, respectively. The high sensitivity and specificity were achieved with a lower limit of detection than that of the routine artificial digestion test, suggesting that serological surveillance is a valuable alternative in surveillance for Trichinella spp. in pig production.

Validation of a Western Blot for the detection of anti-Trichinella spp. antibodies in domestic pigs

Trichinellosis is a zoonotic disease in humans caused by Trichinella spp. According to international regulations and guidelines, serological surveillance can be used to demonstrate the absence of Trichinella spp. in a defined domestic pig population. Most enzyme-linked immunosorbent assay (ELISA) tests presently available do not yield 100% specificity, and therefore, a complementary test is needed to confirm the diagnosis of any initial ELISA seropositivity. The goal of the present study was to evaluate the sensitivity and specificity of a Western Blot assay based on somatic Trichinella spiralis muscle stage (L1) antigen using Bayesian modeling techniques. A total of 295 meat juice and serum samples from pigs negative for Trichinella larvae by artificial digestion, including 74 potentially cross-reactive sera of pigs with other nematode infections, and 93 meat juice samples from pigs infected with Trichinella larvae were included in the study. The diagnostic sensitivity and specificity of the Western Blot were ranged from 95.8% to 96.0% and from 99.5% to 99.6%, respectively. A sensitivity analysis showed that the model outcomes were hardly influenced by changes in the prior distributions, providing a high confidence in the outcomes of the models. This validation study demonstrated that the Western Blot is a suitable method to confirm samples that reacted positively in an initial ELISA.

Precise Orbit Determination of Low Earth Satellites at AIUB using GPS and SLR data

An ever increasing number of low Earth orbiting (LEO) satellites is, or will be, equipped with retro-reflectors for Satellite Laser Ranging (SLR) and on-board receivers to collect observations from Global Navigation Satellite Systems (GNSS) such as the Global Positioning Sys- tem (GPS) and the Russian GLONASS and the European Galileo systems in the future. At the Astronomical Insti- tute of the University of Bern (AIUB) LEO precise or- bit determination (POD) using either GPS or SLR data is performed for a wide range of applications for satellites at different altitudes. For this purpose the classical numeri- cal integration techniques, as also used for dynamic orbit determination of satellites at high altitudes, are extended by pseudo-stochastic orbit modeling techniques to effi- ciently cope with potential force model deficiencies for satellites at low altitudes. Accuracies of better than 2 cm may be achieved by pseudo-stochastic orbit modeling for satellites at very low altitudes such as for the GPS-based POD of the Gravity field and steady-state Ocean Circula- tion Explorer (GOCE).

Analyzing diffuse scattering with supercomputers

Two new approaches to quantitatively analyze diffuse diffraction intensities from faulted layer stacking are reported. The parameters of a probability-based growth model are determined with two iterative global optimization methods: a genetic algorithm (GA) and particle swarm optimization (PSO). The results are compared with those from a third global optimization method, a differential evolution (DE) algorithm [Storn & Price (1997). J. Global Optim. 11, 341–359]. The algorithm efficiencies in the early and late stages of iteration are compared. The accuracy of the optimized parameters improves with increasing size of the simulated crystal volume. The wall clock time for computing quite large crystal volumes can be kept within reasonable limits by the parallel calculation of many crystals (clones) generated for each model parameter set on a super- or grid computer. The faulted layer stacking in single crystals of trigonal three-pointedstar- shaped tris(bicylco[2.1.1]hexeno)benzene molecules serves as an example for the numerical computations. Based on numerical values of seven model parameters (reference parameters), nearly noise-free reference intensities of 14 diffuse streaks were simulated from 1280 clones, each consisting of 96 000 layers (reference crystal). The parameters derived from the reference intensities with GA, PSO and DE were compared with the original reference parameters as a function of the simulated total crystal volume. The statistical distribution of structural motifs in the simulated crystals is in good agreement with that in the reference crystal. The results found with the growth model for layer stacking disorder are applicable to other disorder types and modeling techniques, Monte Carlo in particular.

P-T-t estimation of syn-kinematic strain in low-grade quartz-feldspar bearing rocks using thermodynamic modeling and 40Ar/39Ar dating techniques: example of the Plan-de-Phasy shear zone unit (Briançonnais Zone, Western Alps)

Pressure–Temperature–time (P–T–t) estimates of the syn-kinematic strain at the peak-pressure conditions reached during shallow underthrusting of the Briançonnais Zone in the Alpine subduction zone was made by thermodynamic modelling and 40Ar/39Ar dating in the Plan-de-Phasy unit (SE of the Pelvoux Massif, Western Alps). The dated phengite minerals crystallized syn-kinematically in a shear zone indicating top-to-the-N motion. By combining X-ray mapping with multi-equilibrium calculations, we estimate the phengite crystallization conditions at 270 ± 50 °C and 8.1 ± 2 kbar at an age of 45.9 ± 1.1 Ma. Combining this P–T–t estimate with data from the literature allows us to constrain the timing and geometry of Alpine continental subduction. We propose that the Briançonnais units were scalped on top of the slab during ongoing continental subduction and exhumed continuously until collision.

A fixed-links modeling approach to assess individual differences in the attentional blink: Analysis of behavioral and psychophysiological data

The attentional blink (AB) is a fundamental limitation of the ability to select relevant information from irrelevant information. It can be observed with the detection rate in an AB task as well as with the corresponding P300 amplitude of the event-related potential. In previous research, however, correlations between these two levels of observation were weak and rather inconsistent. A possible explanation of this finding might be that multiple processes underlie the AB and, thus, obscure a possible relationship between AB-related detection rate and the corresponding P300 amplitude. The present study investigated this assumption by applying a fixed-links modeling approach to represent behavioral individual differences in the AB as a latent variable. Concurrently, this approach enabled us to control for additional sources of variance in AB performance by deriving two additional latent variables. The correlation between the latent variable reflecting behavioral individual differences in AB magnitude and a corresponding latent variable derived from the P300 amplitude was high (r=.70). Furthermore, this correlation was considerably stronger than the correlations of other behavioral measures of the AB magnitude with their psychophysiological counterparts (all rs<.40). Our findings clearly indicate that the systematic disentangling of various sources of variance by utilizing the fixed-links modeling approach is a promising tool to investigate behavioral individual differences in the AB and possible psychophysiological correlates of these individual differences.

MODELING THE DYNAMICS OF AIRWAY CONSTRICTION: EFFECTS OF AGONIST TRANSPORT AND BINDING

Recent advances have revealed that during exogenous airway challenge, airway diameters can not be adequately predicted by their initial diameters. Furthermore, airway diameters can also vary greatly in time on scales shorter than a breath. In order to better understand these phenomena, we developed a multiscale model which allows us to simulate aerosol challenge in the airways during ventilation. The model incorporates agonist-receptor binding kinetics to govern the temporal response of airway smooth muscle (ASM) contraction on individual airway segments, which together with airway wall mechanics, determines local airway caliber. Global agonist transport and deposition is coupled with pressure-driven flow, linking local airway constrictions with global flow dynamics. During the course of challenge, airway constriction alters the flow pattern, redistributing agonist to less constricted regions. This results in a negative feedback which may be a protective property of the normal lung. As a consequence, repetitive challenge can cause spatial constriction patterns to evolve in time, resulting in a loss of predictability of airway diameters. Additionally, the model offers new insight into several phenomena including the intra- and inter-breath dynamics of airway constriction throughout the tree structure.

Individual and population level effects of partner notification for Chlamydia trachomatis

Partner notification (PN or contact tracing) is an important aspect of treating bacterial sexually transmitted infections (STIs), such as Chlamydia trachomatis. It facilitates the identification of new infected cases that can be treated through individual case management. PN also acts indirectly by limiting onward transmission in the general population. However, the impact of PN, both at the level of individuals and the population, remains unclear. Since it is difficult to study the effects of PN empirically, mathematical and computational models are useful tools for investigating its potential as a public health intervention. To this end, we developed an individual-based modeling framework called Rstisim. It allows the implementation of different models of STI transmission with various levels of complexity and the reconstruction of the complete dynamic sexual partnership network over any time period. A key feature of this framework is that we can trace an individual's partnership history in detail and investigate the outcome of different PN strategies for C. trachomatis. For individual case management, the results suggest that notifying three or more partners from the preceding 18 months yields substantial numbers of new cases. In contrast, the successful treatment of current partners is most important for preventing re-infection of index cases and reducing further transmission of C. trachomatis at the population level. The findings of this study demonstrate the difference between individual and population level outcomes of public health interventions for STIs.

Labeling the human skeleton with 41Ca to assess changes in bone calcium metabolism

Bone research is limited by the methods available for detecting changes in bone metabolism. While dual X-ray absorptiometry is rather insensitive, biochemical markers are subject to significant intra-individual variation. In the study presented here, we evaluated the isotopic labeling of bone using 41Ca, a long-lived radiotracer, as an alternative approach. After successful labeling of the skeleton, changes in the systematics of urinary 41Ca excretion are expected to directly reflect changes in bone Ca metabolism. A minute amount of 41Ca (100 nCi) was administered orally to 22 postmenopausal women. Kinetics of tracer excretion were assessed by monitoring changes in urinary 41Ca/40Ca isotope ratios up to 700 days post-dosing using accelerator mass spectrometry and resonance ionization mass spectrometry. Isotopic labeling of the skeleton was evaluated by two different approaches: (i) urinary 41Ca data were fitted to an established function consisting of an exponential term and a power law term for each individual; (ii) 41Ca data were analyzed by population pharmacokinetic (NONMEM) analysis to identify a compartmental model that describes urinary 41Ca tracer kinetics. A linear three-compartment model with a central compartment and two sequential peripheral compartments was found to best fit the 41Ca data. Fits based on the use of the combined exponential/power law function describing urinary tracer excretion showed substantially higher deviations between predicted and measured values than fits based on the compartmental modeling approach. By establishing the urinary 41Ca excretion pattern using data points up to day 500 and extrapolating these curves up to day 700, it was found that the calculated 41Ca/40Ca isotope ratios in urine were significantly lower than the observed 41Ca/40Ca isotope ratios for both techniques. Compartmental analysis can overcome this limitation. By identifying relative changes in transfer rates between compartments in response to an intervention, inaccuracies in the underlying model cancel out. Changes in tracer distribution between compartments were modeled based on identified kinetic parameters. While changes in bone formation and resorption can, in principle, be assessed by monitoring urinary 41Ca excretion over the first few weeks post-dosing, assessment of an intervention effect is more reliable approximately 150 days post-dosing when excreted tracer originates mainly from bone.