A review on polycyclic aromatic hydrocarbons distribution in freshwater ecosystems and their toxicity to benthic fauna

Polycyclic aromatic hydrocarbons (PAHs) are a group of organic compounds, found ubiquitously in all environmental compartments. PAHs are considered hazardous pollutants, being of concern to both the environmental and human health. In the aquatic environment, PAHs tend to accumulate in the sediment due to their high hydrophobicity, and thus sediments can be considered their ultimate sink. Concurrently, sediments comprise important habitats for benthic species. This raises concern over the toxic effects of PAHs to benthic communities. Despite PAHs have been the subject of several reviews, their toxicity to freshwater benthic species has not been comprehensively discussed. This review aimed to provide an overview on PAHs distribution in freshwater environments and on their toxicity to benthic fauna species. The distribution of PAHs between sediments and the overlying water column, given by the sediment-water partition coefficient, revealed that PAHs concentrations were 2 to 4 orders of magnitude higher in sediments than in water. The sediment-water partition coefficient was positively correlated to PAHs hydrophobicity. Toxicity of PAHs to benthic fauna was addressed through Species Sensitivity Distributions. The derived hazardous concentration for 5% of the species (HC₅) decreased as follows: NAP (376 μg L^-1) > PHE > PYR > FLT > ANT (0.854 μg L^-1), varying by 3 orders of magnitude. The hazardous concentrations (HC₅) to benthic species were inversely correlated to the hydrophobicity of the individual PAHs. These findings are pertinent for environmental risk assessment of these compounds. This review also identified future challenges regarding the environmental toxicity of PAHs to freshwater benthic communities, namely the need for updating the PAHs priority list and the importance of comprehensively and more realistically assess the toxicity of PAHs in combination with other stressors, both chemical and climate-related.

Finding the match between healthcare worker and expert for optimal audit and feedback on antimicrobial resistance prevention measures

Background

The potentials of audit and feedback (AF) to improve healthcare are currently not exploited. To unlock the potentials of AF, this study focused on the process of making sense of audit data and translating data into actionable feedback by studying a specific AF-case: limiting antimicrobial resistance (AMR). This was done via audit and feedback of AMR prevention measures (APM) that are executed by healthcare workers (HCW) in their day-to-day contact with patients. This study’s aim was to counterbalance the current predominantly top-down, expert-driven audit and feedback approach for APM, with needs and expectations of HCW.

Methods

Qualitative semi-structured interviews were held with sixteen HCW (i.e. physicians, residents and nurses) from high-risk AMR departments at a regional hospital in The Netherlands. Deductive coding was succeeded by open and axial coding to establish main codes, subcodes and variations within codes.

Results

HCW demand insights from audits into all facets of APM in their working routines (i.e. diagnostics, treatment and infection control), preferably in the form of simple and actionable feedback that invites interdisciplinary discussions, so that substantiated actions for improvement can be implemented. AF should not be seen as an isolated ad-hoc intervention, but as a recurrent, long-term, and organic improvement strategy that balances the primary aims of HCW (i.e. improving quality and safety of care for individual patients and HCW) and AMR-experts (i.e. reducing the burden of AMR).

Conclusions

To unlock the learning and improvement potentials of audit and feedback, HCW’ and AMR-experts’ perspectives should be balanced throughout the whole AF-loop (incl. data collection, analysis, visualization, feedback and planning, implementing and monitoring actions). APM-AF should be flexible, so that both audit (incl. collecting and combining the right data in an efficient and transparent manner) and feedback (incl. persuasive and actionable feedback) can be tailored to the needs of various target groups. To balance HCW’ and AMR-experts’ perspectives a participatory holistic AF development approach is advocated.

Applying Infrared Thermography to Soil Surface Temperature Monitoring: Case Study of a High-Resolution 48 h Survey in a Vineyard (Anadia, Portugal)

The soil surface albedo decreases with an increasing biochar application rate as a power decay function, but the net impact of biochar application on soil temperature dynamics remains to be clarified. The objective of this study was to assess the potential of infrared thermography (IRT) sensing by monitoring soil surface temperature (SST) with a high spatiotemporal and thermal resolution in a scalable agricultural application. We monitored soil surface temperature (SST) variations over a 48 h period for three treatments in a vineyard: bare soil (plot S), 100% biochar cover (plot B), and biochar-amended topsoil (plot SB). The SST of all plots was monitored at 30 min intervals with a tripod-mounted IR thermal camera. The soil temperature at 10 cm depth in the S and SB plots was monitored continuously with a 5 min resolution probe. Plot B had greater daily SST variations, reached a higher daily temperature peak relative to the other plots, and showed a faster rate of T increase during the day. However, on both days, the SST of plot B dipped below that of the control treatment (plot S) and biochar-amended soil (plot SB) from about 18:00 onward and throughout the night. The diurnal patterns/variations in the IRT-measured SSTs were closely related to those in the soil temperature at a 10 cm depth, confirming that biochar-amended soils showed lower thermal inertia than the unamended soil. The experiment provided interesting insights into SST variations at a local scale. The case study may be further developed using fully automated SST monitoring protocols at a larger scale for a range of environmental and agricultural applications.

Helping stakeholders select and apply appraisal tools to mitigate soil threats: Researchers' experiences from across Europe

Soil improvement measures need to be ecologically credible, socially acceptable and economically affordable if they are to enter widespread use. However, in real world decision contexts not all measures can sufficiently meet these criteria. As such, developing, selecting and using appropriate tools to support more systematic appraisal of soil improvement measures in different decision-making contexts represents an important challenge. Tools differ in their aims, ranging from those focused on appraising issues of cost-effectiveness, wider ecosystem services impacts and adoption barriers/opportunities, to those seeking to foster participatory engagement and social learning. Despite the growing complexity of the decision-support tool landscape, comprehensive guidance for selecting tools that are best suited to appraise soil improvement measures, as well as those well-adapted to enable participatory deployment, has generally been lacking. We address this gap using the experience and survey data from an EU-funded project (RECARE: Preventing and REmediating degradation of soils in Europe through land CARE). RECARE applied different socio-cultural, biophysical and monetary appraisal tools to assess the costs, benefits and adoption of soil improvement measures across Europe. We focused on these appraisal tools and evaluated their performance against three broad attributes that gauge their differences and suitability for widespread deployment to aid stakeholder decision making in soil management. Data were collected using an online questionnaire administered to RECARE researchers. Although some tools worked better than others across case studies, the information collated was used to provide guiding strategies for choosing appropriate tools, considering resources and data availability, characterisation of uncertainty, and the purpose for which a specific soil improvement measure is being developed or promoted. This paper provides insights to others working in practical soil improvement contexts as to why getting the tools right matters. It demonstrates how use of the right tools can add value to decision-making in ameliorating soil threats, supporting the sustainable management of the services that our soil ecosystems provide.

Cross-border comparison of antimicrobial resistance (AMR) and AMR prevention measures: the healthcare workers' perspective

Background

Cross-border healthcare may promote the spread of multidrug-resistant microorganisms (MDRO) and is challenging due to heterogeneous antimicrobial resistance (AMR) prevention measures (APM). The aim of this article is to compare healthcare workers (HCW) from Germany (DE) and The Netherlands (NL) on how they perceive and experience AMR and APM, which is important for safe patient exchange and effective cross-border APM cooperation.

Methods

A survey was conducted amongst HCW (n = 574) in hospitals in DE (n = 305) and NL (n = 269), using an online self-administered survey between June 2017 and July 2018. Mann-Whitney U tests were used to analyse differences between answers of German and Dutch physicians (n = 177) and German and Dutch nurses (n = 397) on 5-point Likert Items and Scales.

Results

Similarities between DE and NL were a high awareness about the AMR problem and the perception that the possibility to cope with AMR is limited (30% respondents perceive their contribution to limit AMR as insufficient). Especially Dutch nurses scored significantly lower than German nurses on their contribution to limit AMR (means 2.6 vs. 3.1, p ≤ 0.001). German HCW were more optimistic about their potential role in coping with AMR (p ≤ 0.001), and scored higher on feeling sufficiently equipped to perform APM (p ≤ 0.003), although the mean scores did not differ much between German and Dutch respondents.

Conclusions

Although both German and Dutch HCW are aware of the AMR problem, they should be more empowered to contribute to limiting AMR through APM (i.e. screening diagnostics, infection diagnosis, treatment and infection control) in their daily working routines. The observed differences reflect differences in local, national and cross-border structures, and differences in needs of HCW, that need to be considered for safe patient exchange and effective cross-border APM.

Antegrade intramedullary Kirschner-wire fixation of displaced metacarpal shaft fractures

PURPOSE

The objective of this study was to analyze complications and patient-related functional outcome after antegrade intramedullary Kirschner-wire fixation of metacarpal shaft fractures.

METHODS

All consecutive patients treated from January 2010 until December 2015 were retrospectively analyzed using patient logs and radiographic images. Indications for operative fixation were angulation > 40°, shortening > 2 mm, or rotational deficit. Complications were registered from the patient logs. Functional outcome was assessed with the Patient-rated wrist/hand evaluation (PRWHE) and Disabilities of the Arm, Shoulder, and Hand score (DASH) questionnaire both ranging from 1 to 100 after a minimum follow-up of 6 months.

RESULTS

During the study period, 34 fractures of 27 patients could be included. Mean outpatient follow-up was 11 weeks (range 4-24 weeks). The mean interval for functional assessment was 30 months (range 8-62 months) and 19 patients (70%) responded to the questionnaires. During outpatient follow-up, all fractures proceeded to union with no signs of secondary fracture dislocation or implant migration. One re-fracture after a new adequate trauma was seen and one patient underwent tenolysis due to persistent pain and impaired function. In 26 cases (81%), the K-wires were removed of which 23 (68%) were planned removals. Functional outcome was excellent with mean PRWHE and DASH scores of 7 and 5 points, respectively.

CONCLUSIONS

If surgical treatment for metacarpal shaft fractures is considered, we recommend antegrade intramedullary K-wire fixation. This technique results in low complication rates and excellent functional outcome.

THERMOREVERSIBLE CROSS-LINKING OF RUBBER COMPOUNDS: FROM PROOF-OF-CONCEPT TOWARD AN INDUSTRIAL PROCESS

It is demonstrated that the concept of thermoreversible cross-linking of functionalized maleic anhydride grafted ethylene–propylene (EPM-g-MA) rubber using Diels–Alder chemistry is limited neither to laboratory scale using a solvent route nor to gum rubber. The use of an internal mixer is the first step toward an industrial process, since it greatly reduces the processing time and allows for a solventless process for the furan-functionalization and subsequent bismaleimide cross-linking of EPM rubber. Practical rubber compounds were prepared by mixing thermoreversibly cross-linked EPM with carbon black and mineral oil in the same batch mixer. This resulted in reinforcement of the rubber without affecting the thermoreversible character of the cross-linking. The pendant furan groups of the (non)cross-linked EPM-g-furan interact with the carbon black filler. Finally, crystalline EPM rubber compounds were prepared, which show excellent material properties and property retention over multiple reprocessing cycles.

Ulnar styloid process nonunion and outcome in patients with a distal radius fracture: a meta-analysis of comparative clinical trials

PURPOSE

There is no consensus on the relation between ulnar styloid process nonunion and outcome in patients with distal radius fractures. The aim of this study was to analyze whether patient-reported outcome is influenced by the nonunion of the accompanying ulnar styloid fracture in distal radius fracture patients.

METHODS

A meta-analysis of published studies comparing outcomes after distal radius fractures with a united versus a non-united ulnar styloid process was performed. In addition, if provided by the authors, the raw data of these studies were pooled and analysed as one study. The outcome measures of the analyses included patient-reported outcome, functional outcome, grip-strength, pain, and distal radioulnar joint (DRUJ) instability.

RESULTS

Data from six comparative studies were included, concerning 365 patients with a distal radius fracture. One hundred and thirty-five patients with an ulnar styloid union were compared with 230 patients with a nonunion of the ulnar styloid. No significant differences were found between groups regarding any outcome measure.

CONCLUSION

Based on this meta-analysis, there is no relation between the nonunion of the ulnar styloid process and function in patients with a distal radius fracture.

Qualms regarding the range of validity of the glansdorff-prigogine criterion for stability of non-equilibrium States

Doubt is raised concerning the range of validity of a stability criterion for non-equilibrium states which has been proposed by Glansdorff and Prigogine. In the case of a particular autocatalytic reaction, the stability analysis presented by Glansdorff and Prigogine, and by Eigen and by Katchalsky in their reviews of this problem, does not agree with our analysis, which is based upon exact solution of the relevant rate equations. We also find disagreement between the analysis based upon the Glansdorff-Prigogine criterion and our analysis of a second example which involves non-equilibrium steady states. The situation is quite delicate because seemingly innocent approximations (e.g., the use of specialized conditions in the autocatalytic reaction X + Y right arrow over left arrow 2X discussed in the sequel) may lead to the impression that the scope of validity of the criterion is wider than it actually is. By considering the stability of the equilibrium state, we conclude that the second differential of the entropy, which is at the heart of the Glansdorff-Prigogine criterion, is likely to be relevant for stability questions close to equilibrium only.

Maxwell-type constructions for multiple nonequilibrium steady states

A stochastic analysis is used to obtain Maxwell-type constructions for multiple nonequilibrium steady states in homogeneous systems. The stochastic theory shows that homogeneous molecular fluctuations can nucleate the Maxwell construction only on a time scale that is the order of a Poincaré recurrence time and indicates that hysteresis should be observed. The nonrelevance of homogeneous fluctuations for systems that undergo phase separation or in which spatial structures occur is discussed.

On the mechanism of oscillations in the "beating mercury heart"

The discovery of a class of electrochemical-mechanical oscillators-similar to the beating mercury heart -are reported and characterized. Chemical and physical effects on the voltage-time oscillations in these systems have been measured and a mechanism consistent with these effects is given. Interesting cross-catalytic effects involving the electrocapillarity of mercury are proposed to be the origin of the oscillations in these systems.

Modeling temporal behavior of postnatal cat retinal ganglion cells

During development, mammalian retinal ganglion cells (RGCs) go through marked ontogenetic changes with respect to their excitable membrane properties. Voltage-clamp studies conducted in our laboratory have shown that the amplitude, voltage-dependence and kinetics of activation and inactivation (where present) of Na(+), K(+) and Ca(2+) conductances all exhibit developmental changes during a time when the firing patterns of mammalian ganglion cells shift from being transient to being predominantly sustained in nature. In order to better understand the contribution of each conductance to the generation of spikes and spiking patterns, we have developed a model based on our experimental data. For simplicity, we have initially used experimental data obtained from postnatal ganglion cells. At this age the ontogenetic changes observed in the characteristics of the various ionic currents are complete. Utilizing the methods adopted by Hodgkin and Huxley for the giant squid axon, we have determined rate equations for the activation and inactivation properties of the I(A), I(K dr), I(Na), I(Ca L), I(Ca N), and I(leak) currents in postnatal cat RGCs. Combining these with a simplified model of the calcium-activated potassium current (I(KCa)), we have solved and analysed the resulting differential equations. While spikes and spiking patterns resembling experimental data could be obtained from a model in which [Ca(2+)i] was averaged across the whole cell, more accurate simulations were obtained when the diffusion of intracellular Ca(2+) was modeled spatially. The resulting spatial calcium gradients were more effective in gating I(KCa), and our simulations more accurately matched the recorded amplitude and shape of individual spikes as well as the frequency of maintained discharges observed in mammalian postnatal RGCs.

Fire-diffuse-fire model of dynamics of intracellular calcium waves

When Ca2+ is released from internal stores in living cells, the resulting wave of increased concentration can travel without deformation (continuous propagation) or with burst-like behavior (saltatory propagation). We analyze the "fire-diffuse-fire" model in order to illuminate the differences between these two modes of propagation. We show that the Ca2+ release wave in immature Xenopus oocytes and cardiac myocytes is saltatory, whereas the fertilization wave in the mature oocyte is continuous.

Simulation of the fertilization Ca2+ wave in Xenopus laevis eggs

In the preceding paper Fontanilla and Nuccitelli (Biophysical Journal 75:2079-2087 (1998)) present detailed measurements of the shape and speed of the fertilization Ca2+ wave in Xenopus laevis eggs. In order to help interpret their results, we develop here a computational technique based on the finite element method that allows us to carry out realistic simulations of the fertilization wave. Our simulations support the hypothesis that the physiological state of the mature egg is bistable, i.e., that its cytoplasm can accommodate two alternative physiological Ca2+ concentrations: a low concentration characteristic of the prefertilization state and a greatly elevated concentration characteristic of the state following the passage of the wave. We explore this hypothesis by assuming that the bistability is due to the release and re-uptake properties of the endoplasmic reticulum (ER) as determined by inositol trisphosphate (IP3) receptor/Ca2+ channels and sarcoendoplasmic reticulum calcium ATPase (SERCA) pumps. When combined with buffered diffusion of Ca2+ in the cytoplasm, our simulations show that inhomogeneities in the Ca2+ release properties near the plasma membrane are required to explain the temporal and spatial dependences of the shape and speed of these waves. Our results are consistent with an elevated IP3 concentration near the plasma membrane in the unfertilized egg that is augmented significantly near the site of fertilization. These gradients are essential in determining the concave shape of the Ca2+ fertilization wave front.

Saltatory propagation of Ca2+ waves by Ca2+ sparks

Punctate releases of Ca2+, called Ca2+ sparks, originate at the regular array of t-tubules in cardiac myocytes and skeletal muscle. During Ca2+ overload sparks serve as sites for the initiation and propagation of Ca2+ waves in myocytes. Computer simulations of spark-mediated waves are performed with model release sites that reproduce the adaptive Ca2+ release observed for the ryanodine receptor. The speed of these waves is proportional to the diffusion constant of Ca2+, D, rather than D, as is true for reaction-diffusion equations in a continuous excitable medium. A simplified "fire-diffuse-fire" model that mimics the properties of Ca2+-induced Ca2+ release (CICR) from isolated sites is used to explain this saltatory mode of wave propagation. Saltatory and continuous wave propagation can be differentiated by the temperature and Ca2+ buffer dependence of wave speed.

Spark-to-wave transition: saltatory transmission of calcium waves in cardiac myocytes

Using a modular approach, in which kinetic models of various mechanisms of calcium handling in cells are fine-tuned to in vivo and in vitro measurements before combining them into whole-cell models, three distinct modes of transmission of calcium waves in mature and immature frog eggs have been defined. Two modes of transmission are found in immature eggs, where the inositol 1,4,5-trisphosphate receptor (IP3R) controls release of calcium from the endoplasmic reticulum (ER). The first mode corresponds to an excitable physiological state of the cytoplasm and results in solitary waves that can appear as circular or spiral waves in two dimensions with the wave speed proportional to the square root of the diffusion constant of calcium. A second mode occurs when the state of the cytoplasm is oscillatory and because of the small size of the buffered diffusion constant for calcium, the wave speed can appear to be weakly dependent on diffusion. In the mature frog egg, where the sperm-induced Ca2+ fertilization wave is a propagating front, the cytoplasm appears to be bistable and in this mode the wave speed is also proportional to the square root of the diffusion constant. Here we investigate a fourth mode of propagation for cardiac myocytes, in which calcium release from the sarcoplasmic reticulum (SR) is dominated by clusters of ryanodine receptors spaced at regular intervals. In myocytes a stochastically excitable myoplasm leads to the spontaneous production of calcium 'sparks' that under certain conditions can merge into saltatory waves with a speed proportional to the diffusion constant.

Initial high anti-emetic efficacy of granisetron with dexamethasone is not maintained over repeated cycles

We have reported previously that the anti-emetic efficacy of single agent 5HT3 antagonists is not maintained when analysed with the measurement of cumulative probabilities. Presently, the most effective anti-emetic regimen is a combination of a 5HT3 antagonist plus dexamethasone. We, therefore, assessed the sustainment of efficacy of such a combination in 125 patients, scheduled to receive cisplatin > or = 70 mg m(-2) either alone or in combination with other cytotoxic drugs. Anti-emetic therapy was initiated with 10 mg of dexamethasone and 3 mg of granisetron intravenously, before cisplatin. On days 1-6, patients received 8 mg of dexamethasone and 1 mg of granisetron twice daily by oral administration. Protection was assessed during all cycles and calculated based on cumulative probability analyses using the method of Kaplan-Meier and a model for transitional probabilities. Irrespective of the type of analysis used, the anti-emetic efficacy of granisetron/dexamethasone decreased over cycles. The initial complete acute emesis protection rate of 66% decreased to 30% according to the method of Kaplan-Meier and to 39% using the model for transitional probabilities. For delayed emesis, the initial complete protection rate of 52% decreased to 21% (Kaplan-Meier) and to 43% (transitional probabilities). In addition, we observed that protection failure in the delayed emesis period adversely influenced the acute emesis protection in the next cycle. We conclude that the anti-emetic efficacy of a 5HT3 antagonist plus dexamethasone is not maintained over multiple cycles of highly emetogenic chemotherapy, and that the acute emesis protection is adversely influenced by protection failure in the delayed emesis phase.

Model of beta-cell mitochondrial calcium handling and electrical activity. II. Mitochondrial variables

In the preceding article [Am. J. Physiol. 274 (Cell Physiol. 43): C1158-C1173, 1998], we describe the development of a kinetic model for the interaction of mitochondrial Ca2+ handling and electrical activity in the pancreatic beta-cell. Here we describe further results of those simulations, focusing on mitochondrial variables, the rate of respiration, and fluxes of metabolic intermediates as a function of D-glucose concentration. Our simulations predict relatively smooth increases of O2 consumption, adenine nucleotide transport, oxidative phosphorylation, and ATP production by the tricarboxylic acid cycle as D-glucose concentrations are increased from basal to 20 mM. On the other hand, we find that the active fraction of pyruvate dehydrogenase saturates, due to increases in matrix Ca2+, near the onset of bursting electrical activity and that the NADH/NAD+ ratio in the mitochondria increases by roughly an order of magnitude as glucose concentrations are increased. The mitochondrial ATP/ADP ratio increases by factor of < 2 between the D-glucose threshold for bursting and continuous spiking. According to our simulations, relatively small changes in mitochondrial membrane potential (approximately 1 mV) caused by uptake of Ca2+ are sufficient to alter the cytoplasmic ATP/ADP ratio and influence ATP-sensitive K+ channels in the plasma membrane. In the simulations, these cyclic changes in the mitochondrial membrane potential are due to synchronization of futile cycle of Ca2+ from the cytoplasm through mitochondria via Ca2+ uniporters and Na+/Ca2+ exchange. Our simulations predict steady mitochondrial Ca2+ concentrations on the order of 0.1 microM at low glucose concentrations that become oscillatory with an amplitude on the order of 0.5 microM during bursting. Abrupt increases in mitochondrial Ca2+ concentration > 5 microM may occur during continuous electrical activity.

Model of beta-cell mitochondrial calcium handling and electrical activity. I. Cytoplasmic variables

We continue our development of a kinetic model of bursting electrical activity in the pancreatic beta-cell (J. Keizer and G. Magnus. Biophys. J. 56: 229-242, 1989), including the influence of Ca2+ handling by the mitochondria. Our minimal model of mitochondrial Ca2+ handling [G. Magnus and J. Keizer. Am. J. Physiol. 273 (Cell Physiol. 42): C717-C733, 1997] is expanded to include the D-glucose dependence of the rate of production of mitochondrial reducing equivalents. The Ca2+ dependence of the mitochondrial dehydrogenases, which is also included in the model, plays only a small role in the simulations, since the dehydrogenases appear to be maximally activated when D-glucose concentrations are sufficient to produce bursting. A previous model of ionic currents in the plasma membrane is updated using a recent experimental characterization of the dependence of the conductance of the ATP-sensitive K+ (KATP) current on adenine nucleotides. The resulting whole cell model is complex, involving 12 dynamic variables that couple Ca2+ handling in the cytoplasm and the mitochondria with electrical activity in the plasma and inner mitochondrial membranes. Simulations with the whole cell model give rise to bursting electrical activity similar to that seen in pancreatic islets and clusters of pancreatic beta-cells. The full D-glucose dose response of electrical activity is obtained if the cytosolic rate of ATP hydrolysis is a sigmoidal function of glucose. The simulations give the correct shape, period, and phase of the associated oscillations in cytosolic Ca2+, predict that the conductance of the KATP current oscillates out of phase with electrical activity [as recently observed in ob/ob mice (O. Larsson, H. Kindmark, R. Bränstrom, B. Fredholm, and P.-O. Berggren. Proc. Natl. Acad. Sci. USA 93: 5161-5165, 1996)], and make other novel predictions. In this model, bursting results because Ca2+ uptake into mitochondria during the active phase reduces the mitochondrial inner membrane potential, reducing the rate of production of ATP, which in turn activates the KATP current and repolarizes the plasma membrane.

A fluorometric method for estimating the calcium content of internal stores

The concentration of Ca2+ in intracellular stores is an important factor in many aspects of Ca2+ signaling, including the generation of Ca2+ spikes, oscillations and waves, control of mitochondrial respiration, and activation of store-operated Ca2+ channels. Here we describe a consistent method for estimating the content of stores, based on the release of stored Ca2+ by thapsigargin (TG) or ionomycin (IO). Once released from stores, Ca2+ elevates [Ca2+]i transiently before it is pumped across the plasma membrane. If the dependence of the pump rate on [Ca2+]i is known, then the kinetics and amplitude of the Ca2+ transient allows the total amount of releasable Ca2+ to be estimated. We develop this quantitative approach and validate its use in human T cells, in which the Ca2+ clearance rate is an approximately linear function of [Ca2+]i. Our results support the assumption that the ER Ca2+ leak in resting T cells is unregulated, i.e. its rate is proportional to luminal [Ca2+]. The characteristic time constant for basal Ca2+ release is 110-140 s, comparable to that for activation of Ca2+ release-activated Ca2+ (CRAC) channels by TG and consistent with the dependence of ICRAC on store depletion. This method for estimating store content may be useful for quantifying the overlap between functionally distinct stores and for defining the relation between store content and cellular responses.

Agonist-induced calcium waves in oscillatory cells: a biological example of Burgers' equation

Oscillations in cytosolic Ca2+ concentrations in living cells are often a manifestation of propagating waves of Ca2+. Numerical simulations with a realistic model of inositol 1,4,5-trisphosphate (IP3)-induced Ca2+ wave trains lead to wave speeds that increase linearly at long times when (a) IP3 levels are in the range for Ca2+ oscillations, (b) a gradient of phase is established by either an initial ramp or pulse of IP3, and (c) IP3 concentrations asymptotically become uniform. We explore this phenomenon with analytical and numerical methods using a simple two-variable reduction of the De Young-Keizer model of the IP3 receptor that includes the influence of Ca2+ buffers. For concentrations of IP3 in the oscillatory regime, numerical solution of the resulting reaction diffusion equations produces nonlinear wave trains that shows the same asymptotic growth of wave speed. Due to buffering, diffusion of Ca2+ is quite slow and, as previously noted, these waves occur without appreciable bulk movement of Ca2+. Thus, following Neu and Murray, we explore the behavior of these waves using an asymptotic expansion based on the small size of the buffered diffusion constant for Ca2+. We find that the gradient in phase of the wave obeys Burgers' equation asymptotically in time. This result is used to explain the linear increase of the wave speed observed in the simulations.

Minimal model of beta-cell mitochondrial Ca2+ handling

We develop a simplified, but useful, mathematical model to describe Ca2+ handling by mitochondria in the pancreatic beta-cell. The model includes the following six transport mechanisms in the inner mitochondrial membrane: proton pumping via respiration and proton uptake by way of the F1Fzero-ATPase (adapted from D. Pietrobon and S. Caplan. Biochemistry 24: 5764-5778, 1985), a proton leak, adenine nucleotide exchange, the Ca2+ uniporter, and Na+/Ca2+ exchange. Each mechanism is developed separately into a kinetic model for the rate of transport, with parameters taken from experiments on isolated mitochondrial preparations. These mechanisms are combined in a modular fashion first to describe state 4 (nonphosphorylating) and state 3 (phosphorylating) mitochondria with mitochondrial NADH and Ca2+ concentrations as fixed parameters and then to describe Ca2+ handling with variable mitochondrial Ca2+ concentration. Simulations are compared to experimental measurements and agree well with the threshold for Ca2+ uptake, measured mitochondrial Ca2+ levels, and the influence of Ca2+ on oxygen uptake. In the absence of Ca2+ activation of mitochondrial dehydrogenases, the simulations predict a significant reduction in the rate of production of ATP that involves a "short circuit" via Ca2+ uptake through the uniporter. This effect suggests a potential role for mitochondrial Ca2+ handling in determining the ATP-ADP ratio in the pancreatic beta-cell.

Sensing and refilling calcium stores in an excitable cell

Inositol 1,4,5-trisphosphate (IP3)-induced Ca2+ mobilization leads to depletion of the endoplasmic reticulum (ER) and an increase in Ca2+ entry. We show here for the gonadotroph, an excitable endocrine cell, that sensing of ER Ca2+ content can occur without the Ca2+ release-activated Ca2+ current (Icrac), but rather through the coupling of IP3-induced Ca2+ oscillations to plasma membrane voltage spikes that gate Ca2+ entry. Thus we demonstrate that capacitative Ca2+ entry is accomplished through Ca(2+)-controlled Ca2+ entry. We develop a comprehensive model, with parameter values constrained by available experimental data, to simulate the spatiotemporal behavior of agonist-induced Ca2+ signals in both the cytosol and ER lumen of gonadotrophs. The model combines two previously developed models, one for ER-mediated Ca2+ oscillations and another for plasma membrane potential-driven Ca2+ oscillations. Simulations show agreement with existing experimental records of store content, cytosolic Ca2+ concentration ([Ca2+]i), and electrical activity, and make a variety of new, experimentally testable predictions. In particular, computations with the model suggest that [Ca2+]i in the vicinity of the plasma membrane acts as a messenger for ER content via Ca(2+)-activated K+ channels and Ca2+ pumps in the plasma membrane. We conclude that, in excitable cells that do not express Icrac, [Ca2+]i profiles provide a sensitive mechanism for regulating net calcium flux through the plasma membrane during both store depletion and refilling.

Ryanodine receptor adaptation and Ca2+(-)induced Ca2+ release-dependent Ca2+ oscillations

A simplified mechanism that mimics "adaptation" of the ryanodine receptor (RyR) has been developed and its significance for Ca2+(-)induced Ca2+ release and Ca2+ oscillations investigated. For parameters that reproduce experimental data for the RyR from cardiac cells, adaptation of the RyR in combination with sarco/endoplasmic reticulum Ca2+ ATPase Ca2+ pumps in the internal stores can give rise to either low [Cai2+] steady states or Ca2+ oscillations coexisting with unphysiologically high [Cai2+] steady states. In this closed-cell-type model rapid, adaptation-dependent Ca2+ oscillations occur only in limited ranges of parameters. In the presence of Ca2+ influx and efflux from outside the cell (open-cell model) Ca2+ oscillations occur for a wide range of physiological parameter values and have a period that is determined by the rate of Ca2+ refilling of the stores. Although the rate of adaptation of the RyR has a role in determining the shape and the period of the Ca2+ spike, it is not essential for their existence. This is in marked contrast with what is observed for the inositol 1,4,5-trisphosphate receptor for which the biphasic activation and inhibition of its activity by Ca2+ are sufficient to produce oscillations. Results for this model are compared with those based on Ca2+(-)induced Ca2+ release alone in the bullfrog sympathetic neuron. This kinetic model should be suitable for analyzing phenomena associated with "Ca2+ sparks," including their merger into Ca2+ waves in cardiac myocytes.

Modeling plasma membrane and endoplasmic reticulum excitability in pituitary cells

The response of gonadotrophs to secretagogues involves dose-dependent, complex dynamic patterns of electrical activity and inositol 1,4,5-trisphosphate (InsP(3))-induced Ca(2+) mobilization, including pulsatility and oscillations on multiple time scales from milliseconds to minutes. Detailed in vitro experiments have enabled the identification of key mechanisms that underlie the plasma membrane (PM) electrical excitability and endoplasmic reticulum (ER) calcium excitability. We summarize these findings and review computer simulations of a biophysical model that resynthesizes and couples these components and that reproduces quantitatively the observed time courses and dose-response characteristics, as well as effects of various pharamacological manipulations. The theory suggests that cytosolic calcium is the primary messenger in coordinating the PM and ER regenerative behaviors during ER depletion and refilling.

Validity of the rapid buffering approximation near a point source of calcium ions

In the presence of rapid buffers the full reaction-diffusion equations describing Ca2+ transport can be reduced using the rapid buffering approximation to a single transport equation for [Ca2+]. Here we simulate the full and reduced equations, exploring the conditions necessary for the validity of the rapid buffering approximation for an isolated Ca2+ channel or a cluster of channels. Using a point source and performing numerical simulations of different durations, we quantify the error of the rapid buffering approximation as a function of buffer and source parameters as well as the time and spatial scale set by the resolution of confocal microscopic measurements. We carry out simulations of Ca2+ "sparks" and "puffs," both with and without the indicator dye Ca2+ Green-1, and find that the rapid buffering approximation is excellent. These calculations also show that the traditional calculation of [Ca2+] from a fluorescence signal may grossly underestimate the true value of [Ca2+] near a source. Finally, we use the full model to simulate the transient Ca2+ domain near the pore of an open Ca2+ channel in a cell dialyzed with millimolar concentrations of 1,2-bis(2-aminophenoxy)ethane-N,N,N,N-tetraacetic acid or EGTA. In this regime, where the rapid buffering approximation is poor. Neher's equation for the steady-state Ca2+ profile is shown to be a reliable approximation adjacent to the pore.

Effect of Ca2+ influx on intracellular free Ca2+ responses in antigen-stimulated RBL-2H3 cells

We undertake a quantitative investigation of changes in intracellular free Ca2+ concentration ([Ca2+]i) in antigen-stimulated rat basophilic leukemia (RBL-2H3) cells, which include contributions of both Ca2+ store release and Ca2+ influx from the medium. Following Keizer and De Young (J. Keizer and G. De Young. Biophys. J. 61: 649-660, 1992), we develop a highly constrained mathematical model for [Ca2+]i oscillations in RBL-2H3 cells, which includes activation of the inositol trisphosphate receptor (IP3R) by inositol 1,4,5-trisphospate, indirect Ca2+ activation of the IP3R via Ca2+ -dependent activity of phospholipase C-gamma, slow inhibition of the IP3R by cytosolic Ca2+, refilling of Ca2+ stores by a Ca2+ -ATPase (SERCA)-type pump, and a simple representation of the dependence of plasma membrane (PM) fluxes on experimental conditions. Using this full (open cell) model, we simulate [Ca2+]i responses for protocols in which antigen concentration and external Ca2+ are manipulated and compare out calculations with experimental data. In protocol A, cells are stimulated in the presence of external Ca2+, in protocols B and C, cells are stimulated in the absence of external Ca2+, with external Ca2+ later reapplied in protocol C. We are able to reproduce quantitatively the important features of all three protocols, including the dose response of protocol B, the [Ca2+]i response to thapsigargin, and lag time results, and we provide qualitative explanations for the responses derived from our calculations. We also develop a simplified (closed cell) version of the model in which PM fluxes are neglected and total free Ca2+ concentration ([Ca2+]T) is a slowly varying parameter. This permits us to explain in a simple graphical fashion how PM fluxes may influence [Ca2+]i responses in RBH-2H3 cells through modulation of [Ca2+]T.

Mechanisms of calcium oscillations and waves: a quantitative analysis

Oscillations and waves of increased intracellular free calcium concentration ([Ca2+]i) are observed in a wide range of cell types. Because of their inherent nonlinear nature and the consequent unreliability of intuitive approaches, mathematical modeling has an important role to play in the study of these phenomena. One important class of oscillations and waves is dependent on the presence of inositol (1,4,5)-trisphosphate (IP3), which releases Ca2+ from internal stores via the IP3 receptor/Ca2+ channel. With the minimum possible mathematical formalism, we review mechanistic models for IP3-dependent Ca2+ oscillations and waves. These models are based on the regulation of the IP3 receptor by both IP3 and Ca2+, and incorporate experimental data on the steady-state and kinetic properties of the receptor. The extension of the models to describe intracellular and intercellular Ca2+ waves is considered.

Ca2+ excitability of the ER membrane: an explanation for IP3-induced Ca2+ oscillations

Recent research dealing with experiments and theoretical models of Ca2+ excitability of the endoplasmic reticulum (ER) membrane induced by inositol 1,4,5-trisphosphate (IP3) is reviewed. Ca2+ excitability refers to the ability of a small increment of cytoplasmic Ca2+ concentration ([Ca2+]i) to trigger a large [Ca2+]i pulse or oscillations. Such nonlinear regenerative behavior is conferred by the existence of IP3 channels and Ca(2+)-ATPase transporters on the ER membrane, which extends throughout the cytoplasm. Ca2+ excitability resembles the plasma membrane electrical excitability of neurons and other cells: it is driven by the ionic concentration gradient across the ER membrane (higher Ca2+ concentration inside the ER); each [Ca2+]i spike partially consumes the prestored energy that is reestablished through ATP-dependent active transport; and [Ca2+]i, the excitation variable, controls the nonlinear dynamic release rate of ER Ca2+. This review focuses on the kinetic models based on these features and on experiments dealing with the kinetic properties of [Ca2+]i-dependent gating of the IP3 receptor channel. We summarize evidence in favor of two roles for [Ca2+]i in gating the channel's opening: activation at a rapid time scale and inactivation on a slower time scale. Exploiting an analogy to the well-known Hodgkin-Huxley model for neuronal electrical excitability, we show how Ca2+ excitability of the ER membrane can be explained by these gating properties combined with the ER Ca2+ pump activity. The theory's ability to predict is illustrated by comparing calculated with experimental [Ca2+]i responses for pituitary gonadotrophs under various stimulus conditions.

On the roles of Ca2+ diffusion, Ca2+ buffers, and the endoplasmic reticulum in IP3-induced Ca2+ waves

We have investigated the effects of Ca2+ diffusion, mobile and stationary Ca2+ buffers in the cytosol, and Ca2+ handling by the endoplasmic reticulum on inositol 1,4,5-trisphosphate-induced Ca2+ wave propagation. Rapid equilibration of free and bound Ca2+ is used to describe Ca2+ sequestration by buffers in both the cytosol and endoplasmic reticulum (ER) lumen. Cytosolic Ca2+ regulation is based on a kinetic model of the inositol 1,4,5-trisphosphate (IP3) receptor of De Young and Keizer that includes activation and inhibition of the IP3 receptor Ca2+ channel in the ER membrane and SERCA Ca2+ pumps in the ER. Diffusion of Ca2+ in the cytosol and the ER and the breakdown and diffusion of IP3 are also included in our calculations. Although Ca2+ diffusion is severely limited because of buffering, when conditions are chosen just below the threshold for Ca2+ oscillations, a pulse of IP3 or Ca2+ results in a solitary trigger wave that requires diffusion of Ca2+ for its propagation. In the oscillatory regime repetitive wave trains are observed, but for this type of wave neither the wave shape nor the speed is strongly dependent on the diffusion of Ca2+. Local phase differences lead to waves that are predominately kinematic in nature, so that the wave speed (c) is related to the wavelength (lambda) and the period of the oscillations (tau) approximately by the formula c = lambda/tau. The period is determined by features that control the oscillations, including [IP3] and pump activity, which are related to recent experiments. Both solitary waves and wave trains are accompanied by a Ca2+ depletion wave in the ER lumen, similar to that observed in cortical preparations from sea urchin eggs. We explore the effect of endogenous and exogenous Ca2+ buffers on wave speed and wave shape, which can be explained in terms of three distinct effects of buffering, and show that exogenous buffers or Ca2+ dyes can have considerable influence on the amplitude and width of the waves.

Enzymological differences of AChE and diazinon hepatic metabolism: correlation of in vitro data with the selective toxicity of diazinon to fish species

The in vitro hepatic metabolism of diazinon, as well as the sensitivity of the brain acetylcholine esterase, to diazoxon inhibitory action have been studied in order to explain the different toxicity of diazinon to Oncorhynchus mykiss (rainbow trout), Poecilia reticulata (guppy), Brachydanio rerio (zebra fish) and Cyprinus carpio (carp). In spite of a very sensitive acetylcholine esterase the carp is very resistant to diazinon toxicity because of its very low rate of bioactivation and relatively high activity of detoxicating enzymes. The trout is very sensitive towards diazinon in spite of its low activity of bioactivation, because of its lack of detoxicating enzymes and a very sensitive acetylcholine esterase. Diazinon is very toxic for the guppy, because this fish combines a relatively sensitive acetylcholine esterase with a high rate of bioactivation. The zebra fish has the most insensitive acetylcholine esterase, associated with a limited activation rate, thus resulting a rather resistant species. The results obtained indicate that diazinon toxicity differences among the fish species studied can largely be explained in relation to metabolic balances in the liver and with the features of the target enzyme.

InsP3-induced Ca2+ excitability of the endoplasmic reticulum

Oscillations in intracellular Ca2+ can be induced by a variety of cellular signalling processes (Woods et al., 1986; Berridge 1988; Jacob et al., 1988) and appear to play a role in secretion (Stojilković et al., 1994), fertilization (Miyazaki et al., 1993), and smooth muscle contraction (Iino and Tsukioka, 1994). Recently, great progress has been made in understanding the mechanisms involved in a particular class of Ca2+ oscillation, associated with the second messenger inositol 1,4,5-trisphosphate (InsP3) (Berridge, 1993). Working in concert with intracellular Ca2+, InsP3 controls Ca2+ release via the InsP3 receptor in the endoplasmic reticulum (ER) (Berridge and Irvine, 1989). The IP3 receptor is regulated by its coagonists InsP3 and Ca2+, which both activate and inhibit Ca2+ release (Finch et al., 1991; Bezprozvanny et al., 1991; De Young and Keizer, 1992). These processes, together with the periodic activation of Ca2+ uptake into the ER, have been identified as key features in the mechanism of InsP3-induced Ca2+ oscillations in pituitary gonadotrophs (Li et al., 1994), Xenopus laevis oocytes (Lechleiter and Clapham, 1992; Atri et al., 1993), and other cell types (Keizer and De Young, 1993). Earlier discussions and models of InsP3-induced Ca2+ oscillations focused on the nature and number of internal releasable pools of Ca2+ (Goldbeter et al., 1990; Swillens and Mercan, 1990; Somogyi and Stucki, 1991), the importance of oscillations in InsP3 (Meyer and Stryer, 1988), and other issues not based on detailed experimental findings in specific cells types.

Mathematical analysis of a proposed mechanism for oscillatory insulin secretion in perifused HIT-15 cells

Oscillatory secretion of insulin has been observed in many different experimental preparations ranging from pancreatic islets to the whole pancreas. Here we examine the mathematical features underlying a possible model for oscillatory secretion from the perifused, insulin-secreting cell line, HIT-15. The model includes the kinetics of uptake of glucose by GLUT transporters, the rate of glucose metabolism within the cell, and the effect of glucose on the rate of insulin secretion. Putative feedback by insulin on the rate of glucose transport into the cells is treated phenomenologically and leads to insulin oscillations similar to those observed experimentally in HIT cells. The resulting set of ordinary differential equations is simplified by time-scale analysis to a two-variable set of ordinary differential equations. Because of this simplification we can explore, in great detail, the characteristics of the oscillations and their sensitivity to parameter variation using phase plane analysis.

Analysis of possible mechanisms for in vitro oscillations of insulin secretion

We explore possible kinetic mechanisms responsible for the oscillatory (pulsatile) secretion of insulin observed in vitro when pancreatic islets or islet-derived cells are perifused with glucose. Three primary processes are included: 1) glucose stimulation of insulin secretion, controlled by glucokinase; 2) uptake of glucose through GLUT transporters; and 3) glucose metabolism. Perifusion is approximated as a limiting case of a chemical flow reactor. Using experimentally determined rate laws for the three processes, we examine the effects of exogenous insulin as a phenomenological activator and inhibitor of secretion. The resulting differential equations support oscillations using either direct or indirect inhibition by insulin. The oscillations have many of the characteristics observed in vitro, although the indirect model is in better overall agreement with experiment. We conclude that the mechanisms explored here may help explain insulin oscillations for HIT cells, but not for islets, and predict that oscillations can be modulated by including insulin in the perifusion medium.

Diffusion of inositol 1,4,5-trisphosphate but not Ca2+ is necessary for a class of inositol 1,4,5-trisphosphate-induced Ca2+ waves

Combining a realistic model of inositol 1,4,5-trisphosphate (IP3)-induced Ca2+ oscillations with the diffusion of IP3 and buffered diffusion of Ca2+, we have found that diffusion of Ca2+ plays only a minor role in a class of agonist-induced Ca2+ wave trains. These waves are primarily kinematic in nature, with variable wavelengths and speeds that depend primarily on the phase differences between oscillators at different spatial points. The period is set by the steady-state value of IP3, while the wave speed approximately equals the wavelength/period. Ca2+ diffusion, which is much slower than that of IP3 because of endogenous buffers, is shown to have only a small effect on the wave trains and not to be necessary for the apparent wave propagation. Diffusion of IP3 sets the phase gradient responsible for these wave trains, which consist primarily of localized cycles of Ca2+ uptake and release. Our results imply a possible previously undisclosed role for IP3 in cell signaling.

Effects of rapid buffers on Ca2+ diffusion and Ca2+ oscillations

Based on realistic mechanisms of Ca2+ buffering that include both stationary and mobile buffers, we derive and investigate models of Ca2+ diffusion in the presence of rapid buffers. We obtain a single transport equation for Ca2+ that contains the effects caused by both stationary and mobile buffers. For stationary buffers alone, we obtain an expression for the effective diffusion constant of Ca2+ that depends on local Ca2+ concentrations. Mobile buffers, such as fura-2, BAPTA, or small endogenous proteins, give rise to a transport equation that is no longer strictly diffusive. Calculations are presented to show that these effects can modify greatly the manner and rate at which Ca2+ diffuses in cells, and we compare these results with recent measurements by Allbritton et al. (1992). As a prelude to work on Ca2+ waves, we use a simplified version of our model of the activation and inhibition of the IP3 receptor Ca2+ channel in the ER membrane to illustrate the way in which Ca2+ buffering can affect both the amplitude and existence of Ca2+ oscillations.

Calcium oscillations in pituitary gonadotrophs: comparison of experiment and theory

We have developed a mathematical model that describes several aspects of agonist-induced Ca2+ signaling in single pituitary gonadotrophs. Our model is based on fast activation of the inositol 1,4,5-trisphosphate (InsP3) receptor Ca2+ channels at low free cytosolic Ca2+ concentration ([Ca2+]i) and slow inactivation at high [Ca2+]i. Previous work has shown that these gating properties, when combined with a Ca(2+)-ATPase, are sufficient to generate simulated Ca2+ oscillations. The Hodgkin-Huxley-like description we formulate here incorporates these different gating properties explicitly and renders their effects transparent and easy to modulate. We introduce regulatory mechanisms of channel opening which enable the model, both in the absence and in the presence of Ca2+ entry, to give responses to a wide range of agonist doses that are in good agreement with experimental findings, including subthreshold responses, superthreshold oscillations with frequency determined by [InsP3], and nonoscillatory "biphasic" responses followed occasionally by small-amplitude oscillations. A particular added feature of our model, enhanced channel opening by reduced concentration of Ca2+ in the lumen of the endoplasmic reticulum, allows oscillations to continue during pool depletion. The model predicts that ionomycin and thapsigargin can induce oscillations with basal [InsP3] and zero Ca2+ entry, while Ca2+ injection cannot. Responses to specific pairings of sub- or superthreshold stimuli of agonist, ionomycin, and thapsigargin are also correctly predicted. Since this model encompasses a wide range of observed dynamic behaviors within a single framework, based on well-established mechanisms, its relevance should not be restricted to gonadotrophs.

Effect of voltage-gated plasma membrane Ca2+ fluxes on IP3-linked Ca2+ oscillations

We present a theoretical investigation of the effect of Ca2+ influx through ion channels in the plasma membrane on Ca2+ oscillations induced by agonist stimulation of Ca2+ release from internal stores. We expand a recent model of internal Ca2+ oscillations based on activation and inhibition of the IP3-receptor in the endoplasmic reticulum by introducing plasma membrane voltage-gated Ca2+ and K+ channels based on patch-clamp experiments on mouse pancreatic beta cells. Simulations of voltage- and current-clamped experiments are carried out. The amplitude and frequency of the oscillations as well as the sensitivity to agonist are strongly affected by altering the value of the voltage-clamp. We show that the fundamental quantity governing voltage clamp measurements is the inward flux of Ca2+, which acts as a control parameter to alter characteristics of the oscillations. Under current clamp conditions (I(applied) = 0) that lead to continuous spiking electrical activity, we examine the coupling between electrical activity and Ca2+ by introducing a small whole cell Ca(2+)-activated K+ conductance. Results of calculations are similar to those observed in agonist-stimulated beta cells and gonadotrophs. We conclude, however, that glucose-induced bursting and agonist induced bursting in beta cells are distinct phenomena and suggest a mechanism by which agonist-induced Ca2+ oscillations might potentiate insulin secretion. This leads to a general principle for selective signal transduction by Ca2+ oscillations.

The regioselective binding of CHCl3 reactive intermediates to microsomal phospholipids

Microsomal phospholipids (PL) are a good target for the reactive intermediates produced by either the oxidative or the reductive biotransformation of CHCl3 (Testai et al. (1990), Toxicol. Appl. Pharmacol. 104, 496-503). In order to preliminarily characterize the different PL with CHCl3 reactive intermediates, two common methods of PL breakdown have been exploited: the acid-catalyzed transmethylation and the enzymatic hydrolysis with phospholipase C. The results indicated that radioactivity derived from the adducts of PL with the oxidation metabolite, phosgene, partitioned preferentially in the aqueous phase (the ratio of aqueous to organic phase radioactivity contents was about 10); the opposite occurred (ratio about 0.1) when the PL adducts were produced by the reductive process metabolites (dichloromethyl radicals). Therefore, the two methods of PL adduct breakdown can be used to detect and quantitate selectively the two reactive intermediates of CHCl3 biotransformation. The use of phospholipase C, which specifically cleaves the bond between the glyceryl-oxygen and the phosphor atom of PL also gave some structural information. Indeed, the radioactivity partitioning in the aqueous phase after enzymatic hydrolysis of CHCl3 oxidation-associated PL adducts, indicated the selective covalent binding of phosgene residues with the PL polar heads. The clear-cut different partition of radioactivity observed after hydrolysis of PL adducts with CHCl3 reduction intermediates, analogously indicated that dichloromethyl radicals were selectively bound to the PL fatty acyl chains. Using this method we could confirm that in in vitro experimental conditions resembling the physiological status of the liver, both metabolic pathways were concurrently active in hepatic microsomes of B6C3F1 mice. Extents of reactive metabolites similar to those found in B6C3F1 mouse liver microsomes, could be measured in Sprague-Dawley rat liver microsomes only after pretreatment of the animals with PB and incubation with higher CHCl3 concentrations. The toxicological implications of these findings are discussed.

A single-pool inositol 1,4,5-trisphosphate-receptor-based model for agonist-stimulated oscillations in Ca2+ concentration

Relying on quantitative measurements of Ca2+ activation and inhibition of the inositol 1,4,5-trisphosphate (IP3) receptor in the endoplasmic reticulum, we construct a simplified kinetic model to describe the properties of this channel. Selecting rate constants to fit key kinetic and equilibrium data, we find that the model reproduces a variety of in vivo and in vitro experiments. In combination with Ca(2+)-ATPase activity for Ca2+ uptake into the endoplasmic reticulum, the model leads to cytoplasmic oscillations in Ca2+ concentration at fixed IP3 concentration and only a single pool of releasable Ca2+, the endoplasmic reticulum. Incorporation of a positive-feedback mechanism of Ca2+ on IP3 production by phospholipase C enriches the properties of the oscillations and leads to oscillations in Ca2+ concentration accompanied by oscillations in IP3 concentration. We discuss the possible significance of these results for the interpretation of experiments.

Conditional probability analysis for a domain model of Ca(2+)-inactivation of Ca2+ channels

The domain model of Ca2+ inactivation of Ca2+ channels, which has been used to explain rapid inactivation of whole cell Ca2+ currents in pancreatic beta cells, is applied to single-time and conditional open probability measurements on guinea pig ventricular myocyte Ca2+ channels. These two measurements greatly constrain the choice of kinetic constants in the model. Calculations with the model provide a simple quantitative explanation of recent experimental results, including a slow increase in the inactivation rate.

Slow voltage inactivation of Ca2+ currents and bursting mechanisms for the mouse pancreatic beta-cell

Recent whole-cell electrophysiological data concerning the properties of the Ca2+ currents in mouse beta-cells are fitted by a two-current model of Ca2+ channel kinetics. When the beta-cell K+ currents are added to this model, only large modifications of the measured Ca2+ currents will reproduce the bursting pattern normally observed in mouse islets. However, when the measured Ca2+ currents are modified only slightly and used in conjunction with a K+ conductance that can be modulated dynamically by ATP concentration, reasonable bursting is obtained. Under these conditions it is the K-ATP conductance, rather than the slow voltage inactivation of the Ca2+ current, that determines the interburst interval. We find that this latter model can be reconciled with experiments that limit the possible periodic variation of the K-ATP conductance and with recent observations of intracellular Ca2+ bursting in islets.