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.

The flammability of metallic materials in normal and reduced gravity

Metallic materials exposed to oxygen-enriched atmospheres – as commonly used in the medical, aerospace, aviation and numerous chemical processing industries – represent a significant fire hazard which must be addressed during design, maintenance and operation. Hence, accurate knowledge of metallic materials flammability is required. Reduced gravity (i.e. space-based) operations present additional unique concerns, where the absence of gravity must also be taken into account. The flammability of metallic materials has historically been quantified using three standardised test methods developed by NASA, ASTM and ISO. These tests typically involve the forceful (promoted) ignition of a test sample (typically a 3.2 mm diameter cylindrical rod) in pressurised oxygen. A test sample is defined as flammable when it undergoes burning that is independent of the ignition process utilised. In the standardised tests, this is indicated by the propagation of burning further than a defined amount, or „burn criterion.. The burn criterion in use at the onset of this project was arbitrarily selected, and did not accurately reflect the length a sample must burn in order to be burning independent of the ignition event and, in some cases, required complete consumption of the test sample for a metallic material to be considered flammable. It has been demonstrated that a) a metallic material.s propensity to support burning is altered by any increase in test sample temperature greater than ~250-300 oC and b) promoted ignition causes an increase in temperature of the test sample in the region closest to the igniter, a region referred to as the Heat Affected Zone (HAZ). If a test sample continues to burn past the HAZ (where the HAZ is defined as the region of the test sample above the igniter that undergoes an increase in temperature of greater than or equal to 250 oC by the end of the ignition event), it is burning independent of the igniter, and should be considered flammable. The extent of the HAZ, therefore, can be used to justify the selection of the burn criterion. A two dimensional mathematical model was developed in order to predict the extent of the HAZ created in a standard test sample by a typical igniter. The model was validated against previous theoretical and experimental work performed in collaboration with NASA, and then used to predict the extent of the HAZ for different metallic materials in several configurations. The extent of HAZ predicted varied significantly, ranging from ~2-27 mm depending on the test sample thermal properties and test conditions (i.e. pressure). The magnitude of the HAZ was found to increase with increasing thermal diffusivity, and decreasing pressure (due to slower ignition times). Based upon the findings of this work, a new burn criterion requiring 30 mm of the test sample to be consumed (from the top of the ignition promoter) was recommended and validated. This new burn criterion was subsequently included in the latest revision of the ASTM G124 and NASA 6001B international test standards that are used to evaluate metallic material flammability in oxygen. These revisions also have the added benefit of enabling the conduct of reduced gravity metallic material flammability testing in strict accordance with the ASTM G124 standard, allowing measurement and comparison of the relative flammability (i.e. Lowest Burn Pressure (LBP), Highest No-Burn Pressure (HNBP) and average Regression Rate of the Melting Interface(RRMI)) of metallic materials in normal and reduced gravity, as well as determination of the applicability of normal gravity test results to reduced gravity use environments. This is important, as currently most space-based applications will typically use normal gravity information in order to qualify systems and/or components for reduced gravity use. This is shown here to be non-conservative for metallic materials which are more flammable in reduced gravity. The flammability of two metallic materials, Inconel® 718 and 316 stainless steel (both commonly used to manufacture components for oxygen service in both terrestrial and space-based systems) was evaluated in normal and reduced gravity using the new ASTM G124-10 test standard. This allowed direct comparison of the flammability of the two metallic materials in normal gravity and reduced gravity respectively. The results of this work clearly show, for the first time, that metallic materials are more flammable in reduced gravity than in normal gravity when testing is conducted as described in the ASTM G124-10 test standard. This was shown to be the case in terms of both higher regression rates (i.e. faster consumption of the test sample – fuel), and burning at lower pressures in reduced gravity. Specifically, it was found that the LBP for 3.2 mm diameter Inconel® 718 and 316 stainless steel test samples decreased by 50% from 3.45 MPa (500 psia) in normal gravity to 1.72 MPa (250 psia) in reduced gravity for the Inconel® 718, and 25% from 3.45 MPa (500 psia) in normal gravity to 2.76 MPa (400 psia) in reduced gravity for the 316 stainless steel. The average RRMI increased by factors of 2.2 (27.2 mm/s in 2.24 MPa (325 psia) oxygen in reduced gravity compared to 12.8 mm/s in 4.48 MPa (650 psia) oxygen in normal gravity) for the Inconel® 718 and 1.6 (15.0 mm/s in 2.76 MPa (400 psia) oxygen in reduced gravity compared to 9.5 mm/s in 5.17 MPa (750 psia) oxygen in normal gravity) for the 316 stainless steel. Reasons for the increased flammability of metallic materials in reduced gravity compared to normal gravity are discussed, based upon the observations made during reduced gravity testing and previous work. Finally, the implications (for fire safety and engineering applications) of these results are presented and discussed, in particular, examining methods for mitigating the risk of a fire in reduced gravity.

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.

Controlled growth of carbon nanotubes for electronic and photovoltaic applications

Carbon nanotubes (CNTs), experimentally observed for the first time twenty years ago, have triggered an unprecedented research effort, on the account of their astonishing structural, mechanical and electronic properties. Unfortunately, the current inability in predicting the CNTs’ properties and the difficulty in controlling their position on a substrate are often limiting factors for the application of this material in actual devices. This research aims at the creation of specific methodologies for controlled synthesis of CNTs, leading to effectively employ them in various fields of electronics, e.g. photovoltaics. Focused Ion Beam (FIB) patterning of Si surfaces is here proposed as a means for ordering the assembly of vertical-aligned CNTs. With this technique, substrates with specific nano-structured morphologies have been prepared, enabling a high degree of control over CNTs’ position and size. On these nano-structured substrates, the growth of CNTs has been realized by chemical vapor deposition (CVD), i.e. thermal decomposition of hydrocarbon gases over a heated catalyst. The most common materials used as catalysts in CVD are transition metals like Fe and Ni; however, their presence in the CNT products often results in shortcomings for electronic applications, especially for those based on silicon, being the metallic impurities incompatible with very-large-scale integration (VLSI) technology. In the present work the role of Ge dots as an alternative catalysts for CNTs synthesis on Si substrates has been thoroughly assessed, finding a close connection between the catalytic activity of such material and the CVD conditions, which can affect both size and morphology of the dots. Successful CNT growths from Ge dots have been obtained by CVD at temperatures ranging from 750 to 1000°C, with mixtures of acetylene and hydrogen in an argon carrier gas. The morphology of the Si surface is observed to play a crucial role for the outcome of the CNT synthesis: natural (i.e. chemical etching) and artificial (i.e. FIB patterning, nanoindentation) means of altering this morphology in a controlled way have been then explored to optimize the CNTs yield. All the knowledge acquired in this study has been finally applied to synthesize CNTs on transparent conductive electrodes (indium-tin oxide, ITO, coated glasses), for the creation of a new class of anodes for organic photovoltaics. An accurate procedure has been established which guarantees a controlled inclusion of CNTs on ITO films, preserving their optical and electrical properties. By using this set of conditions, a CNTenhanced electrode has been built, contributing to improve the power conversion efficiency of polymeric solar cells.

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.

Linear and nonlinear propagation of localised plasmon in metallic nanostructures

A major challenge in modern photonics and nano-optics is the diffraction limit of light which does not allow field localisation into regions with dimensions smaller than half the wavelength. Localisation of light into nanoscale regions (beyond its diffraction limit) has applications ranging from the design of optical sensors and measurement techniques with resolutions as high as a few nanometres, to the effective delivery of optical energy into targeted nanoscale regions such as quantum dots, nano-electronic and nano-optical devices. This field has become a major research direction over the last decade. The use of strongly localised surface plasmons in metallic nanostructures is one of the most promising approaches to overcome this problem. Therefore, the aim of this thesis is to investigate the linear and non-linear propagation of surface plasmons in metallic nanostructures. This thesis will focus on two main areas of plasmonic research –– plasmon nanofocusing and plasmon nanoguiding. Plasmon nanofocusing – The main aim of plasmon nanofocusing research is to focus plasmon energy into nanoscale regions using metallic nanostructures and at the same time achieve strong local field enhancement. Various structures for nanofocusing purposes have been proposed and analysed such as sharp metal wedges, tapered metal films on dielectric substrates, tapered metal rods, and dielectric V-grooves in metals. However, a number of important practical issues related to nanofocusing in these structures still remain unclear. Therefore, one of the main aims of this thesis is to address two of the most important of issues which are the coupling efficiency and heating effects of surface plasmons in metallic nanostructures. The method of analysis developed throughout this thesis is a general treatment that can be applied to a diversity of nanofocusing structures, with results shown here for the specific case of sharp metal wedges. Based on the geometrical optics approximation, it is demonstrated that the coupling efficiency from plasmons generated with a metal grating into the nanofocused symmetric or quasi-symmetric modes may vary between ~50% to ~100% depending on the structural parameters. Optimal conditions for nanofocusing with the view to minimise coupling and dissipative losses are also determined and discussed. It is shown that the temperature near the tip of a metal wedge heated by nanosecond plasmonic pulses can increase by several hundred degrees Celsius. This temperature increase is expected to lead to nonlinear effects, self-influence of the focused plasmon, and ultimately self-destruction of the metal tip. This thesis also investigates a different type of nanofocusing structure which consists of a tapered high-index dielectric layer resting on a metal surface. It is shown that the nanofocusing mechanism that occurs in this structure is somewhat different from other structures that have been considered thus far. For example, the surface plasmon experiences significant backreflection and mode transformation at a cut-off thickness. In addition, the reflected plasmon shows negative refraction properties that have not been observed in other nanofocusing structures considered to date. Plasmon nanoguiding – Guiding surface plasmons using metallic nanostructures is important for the development of highly integrated optical components and circuits which are expected to have a superior performance compared to their electronicbased counterparts. A number of different plasmonic waveguides have been considered over the last decade including the recently considered gap and trench plasmon waveguides. The gap and trench plasmon waveguides have proven to be difficult to fabricate. Therefore, this thesis will propose and analyse four different modified gap and trench plasmon waveguides that are expected to be easier to fabricate, and at the same time acquire improved propagation characteristics of the guided mode. In particular, it is demonstrated that the guided modes are significantly screened by the extended metal at the bottom of the structure. This is important for the design of highly integrated optics as it provides the opportunity to place two waveguides close together without significant cross-talk. This thesis also investigates the use of plasmonic nanowires to construct a Fabry-Pérot resonator/interferometer. It is shown that the resonance effect can be achieved with the appropriate resonator length and gap width. Typical quality factors of the Fabry- Pérot cavity are determined and explained in terms of radiative and dissipative losses. The possibility of using a nanowire resonator for the design of plasmonic filters with close to ~100% transmission is also demonstrated. It is expected that the results obtained in this thesis will play a vital role in the development of high resolution near field microscopy and spectroscopy, new measurement techniques and devices for single molecule detection, highly integrated optical devices, and nanobiotechnology devices for diagnostics of living cells.

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.