Stochastic models and simulation of ion channel dynamics

The behaviour of ion channels within cardiac and neuronal cells is intrinsically stochastic in nature. When the number of channels is small this stochastic noise is large and can have an impact on the dynamics of the system which is potentially an issue when modelling small neurons and drug block in cardiac cells. While exact methods correctly capture the stochastic dynamics of a system they are computationally expensive, restricting their inclusion into tissue level models and so approximations to exact methods are often used instead. The other issue in modelling ion channel dynamics is that the transition rates are voltage dependent, adding a level of complexity as the channel dynamics are coupled to the membrane potential. By assuming that such transition rates are constant over each time step, it is possible to derive a stochastic differential equation (SDE), in the same manner as for biochemical reaction networks, that describes the stochastic dynamics of ion channels. While such a model is more computationally efficient than exact methods we show that there are analytical problems with the resulting SDE as well as issues in using current numerical schemes to solve such an equation. We therefore make two contributions: develop a different model to describe the stochastic ion channel dynamics that analytically behaves in the correct manner and also discuss numerical methods that preserve the analytical properties of the model.

Hypoxia-mimicking mesoporous bioactive glass scaffolds with controllable cobalt ion release for bone tissue engineering

Low oxygen pressure (hypoxia) plays an important role in stimulating angiogenesis; there are, however, few studies to prepare hypoxia-mimicking tissue engineering scaffolds. Mesoporous bioactive glass (MBG) has been developed as scaffolds with excellent osteogenic properties for bone regeneration. Ionic cobalt (Co) is established as a chemical inducer of hypoxia-inducible factor (HIF)-1α, which induces hypoxia-like response. The aim of this study was to develop hypoxia-mimicking MBG scaffolds by incorporating ionic Co2+ into MBG scaffolds and investigate if the addition of Co2+ ions would induce a cellular hypoxic response in such a tissue engineering scaffold system. The composition, microstructure and mesopore properties (specific surface area, nano-pore volume and nano-pore distribution) of Co-containing MBG (Co-MBG) scaffolds were characterized and the cellular effects of Co on the proliferation, differentiation, vascular endothelial growth factor (VEGF) secretion, HIF-1α expression and bone-related gene expression of human bone marrow stromal cells (BMSCs) in MBG scaffolds were systematically investigated. The results showed that low amounts of Co (< 5%) incorporated into MBG scaffolds had no significant cytotoxicity and that their incorporation significantly enhanced VEGF protein secretion, HIF-1α expression, and bone-related gene expression in BMSCs, and also that the Co-MBG scaffolds support BMSC attachment and proliferation. The scaffolds maintain a well-ordered mesopore channel structure and high specific surface area and have the capacity to efficiently deliver antibiotics drugs; in fact, the sustained released of ampicillin by Co-MBG scaffolds gives them excellent anti-bacterial properties. Our results indicate that incorporating cobalt ions into MBG scaffolds is a viable option for preparing hypoxia-mimicking tissue engineering scaffolds and significantly enhanced hypoxia function. The hypoxia-mimicking MBG scaffolds have great potential for bone tissue engineering applications by combining enhanced angiogenesis with already existing osteogenic properties.

The undesirable acetylation of cellulose by the acetate ion of 1-ethyl-3-methylimidazolium acetate

The ionic liquid (IL) 1-ethyl-3-methylimidazolium acetate ([C2mim]OAc) is considered to be an inert solvent of cellulose and lignocellulosic biomass. Acetylation (1.7 % mol, or DS 0.017) of cellulose after dissolution in [C2mim]OAc (150 °C for 20 min), is demonstrated by compositional analysis, FTIR analysis and 13C NMR spectroscopy (in [C2min]OAc with 13C enriched acetate). This acetylation, in the absence of added acylating agents, has not been reported before and may limit [C2mim]OAc utility in industrial scale biomass processing, even at this low extent. For example, cellulose acetylation may contribute to IL loss in processes where the IL is recovered and reused and inhibit enzyme saccharification of cellulose in lignocellulosic biofuel production processes based on saccharification and fermentation.

Modeling ion channel dynamics through reflected stochastic differential equations

Ion channels are membrane proteins that open and close at random and play a vital role in the electrical dynamics of excitable cells. The stochastic nature of the conformational changes these proteins undergo can be significant, however current stochastic modeling methodologies limit the ability to study such systems. Discrete-state Markov chain models are seen as the "gold standard," but are computationally intensive, restricting investigation of stochastic effects to the single-cell level. Continuous stochastic methods that use stochastic differential equations (SDEs) to model the system are more efficient but can lead to simulations that have no biological meaning. In this paper we show that modeling the behavior of ion channel dynamics by a reflected SDE ensures biologically realistic simulations, and we argue that this model follows from the continuous approximation of the discrete-state Markov chain model. Open channel and action potential statistics from simulations of ion channel dynamics using the reflected SDE are compared with those of a discrete-state Markov chain method. Results show that the reflected SDE simulations are in good agreement with the discrete-state approach. The reflected SDE model therefore provides a computationally efficient method to simulate ion channel dynamics while preserving the distributional properties of the discrete-state Markov chain model and also ensuring biologically realistic solutions. This framework could easily be extended to other biochemical reaction networks. © 2012 American Physical Society.

Microcosm study of ion mobilization and greenhouse gas evolution in soils of a plantation-forested subtropical coastal catchment

This study examined the potential for Fe mobilization and greenhouse gas (GHG, e.g. CO2, and CH4) evolution in SEQ soils associated with a range of plantation forestry practices and water-logged conditions. Intact, 30-cm-deep soil cores collected from representative sites were saturated and incubated for 35 days in the laboratory, with leachate and headspace gas samples periodically collected. Minimal Fe dissolution was observed in well-drained sand soils associated with mature, first-rotation Pinus and organic Fe complexation, whereas progressive Fe dissolution occurred over 14 days in clear-felled and replanted Pinus soils with low organic matter and non-crystalline Fe fractions. Both CO2 and CH4 effluxes were relatively lower in clear-felled and replanted soils compared with mature, first-rotation Pinus soils, despite the lack of statistically significant variations in total GHG effluxes associated with different forestry practices. Fe dissolution and GHG evolution in low-lying, water-logged soils adjacent to riparian and estuarine, native-vegetation buffer zones were impacted by mineral and physical soil properties. Highest levels of dissolved Fe and GHG effluxes resulted from saturation of riparian loam soils with high Fe and clay content, as well as abundant organic material and Fe-metabolizing bacteria. Results indicate Pinus forestry practices such as clear-felling and replanting may elevate Fe mobilization while decreasing CO2 and CH4 emissions from well-drained, SEQ plantation soils upon heavy flooding. Prolonged water-logging accelerates bacterially mediated Fe cycling in low-lying, clay-rich soils, leading to substantial Fe dissolution, organic matter mineralization, and CH4 production in riparian native-vegetation buffer zones.

Comparison of the binding stoichiometries of positively charged DNA-binding drugs using positive and negative ion electrospray ionization mass spectrometry

Positive and negative ion electrospray ionization (ESI) mass spectra of complexes of positively charged small molecules (distamycin, Hoechst 33258, [Ru(phen)2dpq]Cl2 and [Ru(phen)2dpqC]Cl2) have been compared. [Ru(phen)2dpq]Cl2 and [Ru(phen)2dpqC]Cl2 bind to DNA by intercalation. Negative ion ESI mass spectra of mixtures of [Ru(phen)2dpq]Cl2 or [Ru(phen)2dpqC]Cl2 with DNA showed ions from DNA-ligand complexes consistent with solution studies. In contrast, only ions from freeDNAwere present in positive ion ESI mass spectra of mixtures of [Ru(phen)2dpq]Cl2 or [Ru(phen)2dpqC]Cl2 with DNA, highlighting the need for obtaining ESI mass spectra of non-covalent complexes under a range of experimental conditions. Negative ion spectra of mixtures of the minor groove binder Hoechst 33258 with DNA containing a known minor groove binding sequence were dominated by ions from a 1:1 complex. In contrast, in positive ion spectra there were also ions present from a 2:1 (Hoechst 33258: DNA) complex, suggesting an alternative binding mode was possible either in solution or in the gas phase. When Hoechst 33258 was mixed with a DNA sequence lacking a high affinity minor groove binding site, the negative ion ESI mass spectra showed that 1:1 and 2:1 complexes were formed, consistent with existence of binding modes other than minor groove binding. The data presented suggest that comparison of positive and negative ion ESI-MS spectra might provide an insight into various binding modes in both solution and the gas phase.

IL-18 inhibits growth of murine orthotopic prostate carcinomas via both adaptive and innate immune mechanisms

Interleukin(IL)-18 is a pleiotrophic cytokine with functions in immune modulation, angiogenesis and bone metabolism. In this study, the potential of IL-18 as an immunotherapy for prostate cancer (PCa) was examined using the murine model of prostate carcinoma, RM1 and a bone metastatic variant RM1(BM)/B4H7-luc. RM1 and RM1(BM)/B4H7-luc cells were stably transfected to express bioactive IL-18. These cells were implanted into syngeneic immunocompetent mice, with or without an IL-18-neutralising antibody (αIL-18, SK113AE4). IL-18 significantly inhibited the growth of both subcutaneous and orthotopic RM1 tumors and the IL-18 neutralizing antibody abrogated the tumor growth-inhibition. In vivo neutralization of interferon-gamma (IFN-γ) completely eliminated the anti-tumor effects of IL-18 confirming an essential role of IFN-γ as a down-stream mediator of the anti-tumor activity of IL-18. Tumors from mice in which IL-18 and/or IFN-γ was neutralized contained significantly fewer CD4+ and CD8+ T cells than those with functional IL-18. The essential role of adaptive immunity was demonstrated as tumors grew more rapidly in RAG1−/− mice or in mice depleted of CD4+ and/or CD8+ cells than in normal mice. The tumors in RAG1−/− mice were also significantly smaller when IL-18 was present, indicating that innate immune mechanisms are involved. IL-18 also induced an increase in tumor infiltration of macrophages and neutrophils but not NK cells. In other experiments, direct injection of recombinant IL-18 into established tumors also inhibited tumor growth, which was associated with an increase in intratumoral macrophages, but not T cells. These results suggest that local IL-18 in the tumor environment can significantly potentiate anti-tumor immunity in the prostate and clearly demonstrate that this effect is mediated by innate and adaptive immune mechanisms.

A two-metal ion mechanism operates in the hammerhead ribozyme-mediated cleavage of an RNA substrate

Evidence for a two-metal ion mechanism for cleavage of the HH16 hammerhead ribozyme is provided by monitoring the rate of cleavage of the RNA substrate as a function of La3+ concentration in the presence of a constant concentration of Mg2+. We show that a bell-shaped curve of cleavage activation is obtained as La3+ is added in micromolar concentrations in the presence of 8 mM Mg2+, with a maximal rate of cleavage being attained in the presence of 3 microM La3+. These results show that two-metal ion binding sites on the ribozyme regulate the rate of the cleavage reaction and, on the basis of earlier estimates of the Kd values for Mg2+ of 3.5 mM and > 50 mM, that these sites bind La3+ with estimated Kd values of 0.9 and > 37.5 microM, respectively. Furthermore, given the very different effects of these metal ions at the two binding sites, with displacement of Mg2+ by La3+ at the stronger (relative to Mg2+) binding site activating catalysis and displacement of Mg2+ by La3+ at the weaker (relative to Mg2+) (relative to Mg2+) binding site inhibiting catalysis, we show that the metal ions at these two sites play very different roles. We argue that the metal ion at binding site 1 coordinates the attacking 2'-oxygen species in the reaction and lowers the pKa of the attached proton, thereby increasing the concentration of the attacking alkoxide nucleophile in an equilibrium process. In contrast, the role of the metal ion at binding site 2 is to catalyze the reaction by absorbing the negative charge that accumulates at the leaving 5'-oxygen in the transition state. We suggest structural reasons why the Mg(2+)-La3+ ion combination is particularly suited to demonstrating these different roles of the two-metal ions in the ribozyme cleavage reaction.

Observations on catalysis by hammerhead ribozymes are consistent with a two-divalent-metal-ion mechanism

Significant cleavage by hammerhead ribozymes requires activation by divalent metal ions. Several models have been proposed to account for the influence of metal ions on hammerhead activity. A number of recent papers have presented data that have been interpreted as supporting a one-metal-hydroxide-ion mechanism. In addition, a solvent deuterium isotope effect has been taken as evidence against a proton transfer in the rate-limiting step of the cleavage reaction. We propose that these data are more easily explained by a two-metal-ion mechanism that does not involve a metal hydroxide, but does involve a proton transfer in the rate-limiting step.

Monitoring charged particle and ion emissions from a laser printer

While the emission rate of ultrafine particles has been measured and quantified, there is very little information on the emission rates of ions and charged particles from laser printers. This paper describes a methodology that can be adopted for measuring the surface charge density on printed paper and the ion and charged particle emissions during operation of a high-emitting laser printer and shows how emission rates of ultrafine particles, ions and charged particles may be quantified using a controlled experiment within a closed chamber.