Fractional diffusion models of cardiac electrical propagation: Role of structural heterogeneity in dispersion of repolarization

Impulse propagation in biological tissues is known to be modulated by structural heterogeneity. In cardiac muscle, improved understanding on how this heterogeneity influences electrical spread is key to advancing our interpretation of dispersion of repolarization. We propose fractional diffusion models as a novel mathematical description of structurally heterogeneous excitable media, as a means of representing the modulation of the total electric field by the secondary electrical sources associated with tissue inhomogeneities. Our results, analysed against in vivo human recordings and experimental data of different animal species, indicate that structural heterogeneity underlies relevant characteristics of cardiac electrical propagation at tissue level. These include conduction effects on action potential (AP) morphology, the shortening of AP duration along the activation pathway and the progressive modulation by premature beats of spatial patterns of dispersion of repolarization. The proposed approach may also have important implications in other research fields involving excitable complex media.

Unimolecular reaction chemistry of a charge-tagged beta-hydroxyperoxyl radical

β-Hydroxyperoxyl radicals are formed during atmospheric oxidation of unsaturated volatile organic compounds such as isoprene. They are intermediates in the combustion of alcohols. In these environments the unimolecular isomerization and decomposition of β-hydroxyperoxyl radicals may be of importance, either through chemical or thermal activation. We have used ion-trap mass spectrometry to generate the distonic charge-tagged β-hydroxyalkyl radical anion, ˙CH2C(OH)(CH3)CH2C(O)O−, and investigated its subsequent reaction with O2 in the gas phase under conditions that are devoid of complicating radical–radical reactions. Quantum chemical calculations and master equation/RRKM theory modeling are used to rationalize the results and discern a reaction mechanism. Reaction is found to proceed via initial hydrogen abstraction from the γ-methylene group and from the β-hydroxyl group, with both reaction channels eventually forming isobaric product ions due to loss of either ˙OH + HCHO or ˙OH + CO2. Isotope labeling studies confirm that a 1,5-hydrogen shift from the β-hydroxyl functionality results in a hydroperoxyalkoxyl radical intermediate that can undergo further unimolecular dissociations. Furthermore, this study confirms that the facile decomposition of β-hydroxyperoxyl radicals can yield ˙OH in the gas phase.

A stochastic exponential Euler scheme for simulation of stiff biochemical reaction systems

In order to simulate stiff biochemical reaction systems, an explicit exponential Euler scheme is derived for multidimensional, non-commutative stochastic differential equations with a semilinear drift term. The scheme is of strong order one half and A-stable in mean square. The combination with this and the projection method shows good performance in numerical experiments dealing with an alternative formulation of the chemical Langevin equation for a human ether a-go-go related gene ion channel mode

Revisiting the CO oxidation reaction on various Au/TiO2catalysts: Roles of the surface OH groups and the reaction mechanism

This work aims to understand the influence of TiO2 surface structure in Au/TiO2 catalysts on CO oxidation. Au nanoparticles (3 wt%) in the range of 4 to 8 nm were loaded onto four kinds of TiO2 surfaces, which had different surface structures and were synthesized by calcining hydrogen titanate nanotubes at various temperatures and in different atmospheres. The Au catalyst supported on anatase nanorods exhibited the highest activity in CO oxidation at 30 °C among all the five Au/TiO2 catalysts including the reference catalyst of Au/TiO2-P25. X-ray photoelectron spectroscopy (XPS) and infrared emission spectra (IES) results indicate that the anatase nanorods have the most active surface on which water molecules can be strongly adsorbed and OH groups can be formed readily. Theoretical calculation indicates that the surface OH can facilitate the O2 adsorption on the anatase surface. Such active surface features are conducive to the O2 activation and CO oxidation

Numerical analysis of a new space-time variable fractional order advection-dispersion equation

Many physical processes appear to exhibit fractional order behavior that may vary with time and/or space. The continuum of order in the fractional calculus allows the order of the fractional operator to be considered as a variable. In this paper, we consider a new space–time variable fractional order advection–dispersion equation on a finite domain. The equation is obtained from the standard advection–dispersion equation by replacing the first-order time derivative by Coimbra’s variable fractional derivative of order α(x)∈(0,1]α(x)∈(0,1], and the first-order and second-order space derivatives by the Riemann–Liouville derivatives of order γ(x,t)∈(0,1]γ(x,t)∈(0,1] and β(x,t)∈(1,2]β(x,t)∈(1,2], respectively. We propose an implicit Euler approximation for the equation and investigate the stability and convergence of the approximation. Finally, numerical examples are provided to show that the implicit Euler approximation is computationally efficient.

Combining the diaza-Diels-Alder reaction and palladium-catalyzed aminations to prepare amino-substituted porphyrins

A strategy to tackle the synthesis of azoporphyrins with unsubstituted terminal meso positions was investigated. It comprised the combination of diaza-Diels–Alder (DADA) reaction of 1,3-dienes with dialkyl azodicarboxylates, decarboxylative hydrolysis of the bis(carbamates), palladium-catalyzed amination of bromoporphyrin precursors, and retro-DADA reactions to release the ultimate targets. The somewhat confused historical results on the DADA reactions of 1,3-cyclohexadiene were clarified, but the hydrolyses yielded extremely air-sensitive amines which decomposed completely in minutes via autooxidation and retro-DADA reaction. With anthracene or 2,3-dimethyl-1,3-butadiene as the diene, the synthesis of azoporphyrin was not achieved but three amino-substituted porphyrins were obtained in moderate yields under mild conditions. The X-ray crystal structures of several of the intermediates and the final aminoanthracene-porphyrin nickel(II) complex were determined.

Functional localization of the human color center by decreased water displacement using diffusion-weighted fMRI

Introduction Decreased water displacement following increased neural activity has been observed using diffusion-weighted functional MRI (DfMRI) at high b-values. The physiological mechanisms underlying the diffusion signal change may be unique from the standard blood oxygenation level-dependent (BOLD) contrast and closer to the source of neural activity. Whether DfMRI reflects neural activity more directly than BOLD outside the primary cerebral regions remains unclear. Methods Colored and achromatic Mondrian visual stimuli were statistically contrasted to functionally localize the human color center Area V4 in neurologically intact adults. Spatial and temporal properties of DfMRI and BOLD activation were examined across regions of the visual cortex. Results At the individual level, DfMRI activation patterns showed greater spatial specificity to V4 than BOLD. The BOLD activation patterns were more prominent in the primary visual cortex than DfMRI, where activation was localized to the ventral temporal lobe. Temporally, the diffusion signal change in V4 and V1 both preceded the corresponding hemodynamic response, however the early diffusion signal change was more evident in V1. Conclusions DfMRI may be of use in imaging applications implementing cognitive subtraction paradigms, and where highly precise individual functional localization is required.

In situ reaction furnace for real-time XRD studies

The new furnace at the Materials Characterization by X-ray Diffraction beamline at Elettra has been designed for powder diffraction measurements at high temperature (up to 1373 K at the present state). Around the measurement region the geometry of the radiative heating element assures a negligible temperature gradient along the capillary and can accommodate either powder samples in capillary or small flat samples. A double capillary holder allows flow-through of gas in the inner sample capillary while the outer one serves as the reaction chamber. The furnace is coupled to a translating curved imaging-plate detector, allowing the collection of diffraction patterns up to 2[theta] [asymptotically equal to] 130°.

Metal-free graphitic carbon nitride as mechano-catalyst for hydrogen evolution reaction

Graphitic carbon nitride (g-C3N4), as a promising metal-free catalyst for photo-catalytic and electrochemical water splitting, has recently attracted tremendous research interest. However, the underlying catalytic mechanism for the hydrogen evolution reaction (HER) is not fully understood. By using density functional theory calculations, here we have established that the binding free energy of hydrogen atom (ΔGH∗0) on g-C3N4 is very sensitive to mechanical strain, leading to substantial tuning of the HER performance of g-C3N4 at different coverages. The experimentally-observed high HER activity in N-doped graphene supported g-C3N4 (Zheng et al., 2014) is actually attributed to electron-transfer induced strain. A more practical strategy to induce mechanical strain in g-C3N4 is also proposed by doping a bridge carbon atom in g-C3N4 with an isoelectronic silicon atom. The calculated ΔGH∗0 on the Si-doped g-C3N4 is ideal for HER. Our results indicate that g-C3N4 would be an excellent metal-free mechano-catalyst for HER and this finding is expected to guide future experiments to efficiently split water into hydrogen based on the g-C3N4 materials.

Preliminary investigations into the use of a functionalised polymer to reduce diffusion in Fricke gel dosimeters

Purpose A modification of the existing PVA-​FX hydrogel has been made to investigate the use of a functionalised polymer in a Fricke gel dosimetry system to decrease Fe3+ diffusion. Methods The chelating agent, xylenol orange, was chem. bonded to the gelling agent, polyvinyl alc. (PVA) to create xylenol orange functionalised PVA (XO-​PVA)​. A gel was created from the XO-​PVA (20​% w​/v) with ferrous sulfate (0.4 mM) and sulfuric acid (50 mM)​. Results This resulted in an optical d. dose sensitivity of 0.014 Gy-​1, an auto-​oxidn. rate of 0.0005 h-​1, and a diffusion rate of 0.129 mm2 h-​1; an 8​% redn. compared to the original PVA-​FX gel, which in practical terms adds approx. 1 h to the time span between irradn. and accurate read-​out. Conclusions Because this initial method of chem. bonding xylenol orange to polyvinyl alc. has inherently low conversion, the improvement on existing gel systems is minimal when compared to the drawbacks. More efficient methods of functionalising polyvinyl alc. with xylenol orange must be developed for this system to gain clin. relevance.

Direct electrochemical formation of nanostructured amorphous Co(OH)2 on gold electrodes with enhanced activity for the oxygen evolution reaction.

The oxides of cobalt have recently been shown to be highly effective electrocatalysts for the oxygen evolution reaction (OER) under alkaline conditions. In general species such as Co3O4 and CoOOH have been investigated that often require an elevated temperature step during their synthesis to create crystalline materials. In this work we investigate the rapid and direct electrochemical formation of amorphous nanostructured Co(OH)2 on gold electrodes under room temperture conditions which is a highly active precursor for the OER. During the OER some conversion to crystalline Co3O4 occurs at the surface, but the bulk of the material remains amorphous. It is found that the underlying gold electrode is crucial to the materials enhanced performance and provides higher current density than can be achieved using carbon, palladium or copper support electrodes. This catalyst exhibits excellent activity with a current density of 10 mA cm-2 at an overpotential of 360 mV with a high turnover frequency of 2.1 s-1 in 1 M NaOH. A Tafel slope of 56 mV dec-1 at low overpotentials and a slope of 122 mV dec-1 at high overpotentials is consistent with the dual barrier model for the electrocatalytic evolution of oxygen. Significantly, the catalyst maintains excellent activity for up to 24 hr of continuous operation and this approach offers a facile way to create a highly effective and stable material.

Exact solutions of coupled multispecies linear reaction–diffusion equations on a uniformly growing domain

Embryonic development involves diffusion and proliferation of cells, as well as diffusion and reaction of molecules, within growing tissues. Mathematical models of these processes often involve reaction–diffusion equations on growing domains that have been primarily studied using approximate numerical solutions. Recently, we have shown how to obtain an exact solution to a single, uncoupled, linear reaction–diffusion equation on a growing domain, 0 < x < L(t), where L(t) is the domain length. The present work is an extension of our previous study, and we illustrate how to solve a system of coupled reaction–diffusion equations on a growing domain. This system of equations can be used to study the spatial and temporal distributions of different generations of cells within a population that diffuses and proliferates within a growing tissue. The exact solution is obtained by applying an uncoupling transformation, and the uncoupled equations are solved separately before applying the inverse uncoupling transformation to give the coupled solution. We present several example calculations to illustrate different types of behaviour. The first example calculation corresponds to a situation where the initially–confined population diffuses sufficiently slowly that it is unable to reach the moving boundary at x = L(t). In contrast, the second example calculation corresponds to a situation where the initially–confined population is able to overcome the domain growth and reach the moving boundary at x = L(t). In its basic format, the uncoupling transformation at first appears to be restricted to deal only with the case where each generation of cells has a distinct proliferation rate. However, we also demonstrate how the uncoupling transformation can be used when each generation has the same proliferation rate by evaluating the exact solutions as an appropriate limit.

Exact calculations of survival probability for diffusion on growing lines, disks, and spheres: The role of dimension

Unlike standard applications of transport theory, the transport of molecules and cells during embryonic development often takes place within growing multidimensional tissues. In this work, we consider a model of diffusion on uniformly growing lines, disks, and spheres. An exact solution of the partial differential equation governing the diffusion of a population of individuals on the growing domain is derived. Using this solution, we study the survival probability, S(t). For the standard nongrowing case with an absorbing boundary, we observe that S(t) decays to zero in the long time limit. In contrast, when the domain grows linearly or exponentially with time, we show that S(t) decays to a constant, positive value, indicating that a proportion of the diffusing substance remains on the growing domain indefinitely. Comparing S(t) for diffusion on lines, disks, and spheres indicates that there are minimal differences in S(t) in the limit of zero growth and minimal differences in S(t) in the limit of fast growth. In contrast, for intermediate growth rates, we observe modest differences in S(t) between different geometries. These differences can be quantified by evaluating the exact expressions derived and presented here.

Estimating cell diffusivity and cell proliferation rate by interpreting IncuCyte ZOOM™ assay data using the Fisher-Kolmogorov model

Background: Standard methods for quantifying IncuCyte ZOOM™ assays involve measurements that quantify how rapidly the initially-vacant area becomes re-colonised with cells as a function of time. Unfortunately, these measurements give no insight into the details of the cellular-level mechanisms acting to close the initially-vacant area. We provide an alternative method enabling us to quantify the role of cell motility and cell proliferation separately. To achieve this we calibrate standard data available from IncuCyte ZOOM™ images to the solution of the Fisher-Kolmogorov model. Results: The Fisher-Kolmogorov model is a reaction-diffusion equation that has been used to describe collective cell spreading driven by cell migration, characterised by a cell diffusivity, D, and carrying capacity limited proliferation with proliferation rate, λ, and carrying capacity density, K. By analysing temporal changes in cell density in several subregions located well-behind the initial position of the leading edge we estimate λ and K. Given these estimates, we then apply automatic leading edge detection algorithms to the images produced by the IncuCyte ZOOM™ assay and match this data with a numerical solution of the Fisher-Kolmogorov equation to provide an estimate of D. We demonstrate this method by applying it to interpret a suite of IncuCyte ZOOM™ assays using PC-3 prostate cancer cells and obtain estimates of D, λ and K. Comparing estimates of D, λ and K for a control assay with estimates of D, λ and K for assays where epidermal growth factor (EGF) is applied in varying concentrations confirms that EGF enhances the rate of scratch closure and that this stimulation is driven by an increase in D and λ, whereas K is relatively unaffected by EGF. Conclusions: Our approach for estimating D, λ and K from an IncuCyte ZOOM™ assay provides more detail about cellular-level behaviour than standard methods for analysing these assays. In particular, our approach can be used to quantify the balance of cell migration and cell proliferation and, as we demonstrate, allow us to quantify how the addition of growth factors affects these processes individually.

Survival probability for a diffusive process on a growing domain

We consider the motion of a diffusive population on a growing domain, 0 < x < L(t ), which is motivated by various applications in developmental biology. Individuals in the diffusing population, which could represent molecules or cells in a developmental scenario, undergo two different kinds of motion: (i) undirected movement, characterized by a diffusion coefficient, D, and (ii) directed movement, associated with the underlying domain growth. For a general class of problems with a reflecting boundary at x = 0, and an absorbing boundary at x = L(t ), we provide an exact solution to the partial differential equation describing the evolution of the population density function, C(x,t ). Using this solution, we derive an exact expression for the survival probability, S(t ), and an accurate approximation for the long-time limit, S = limt→∞ S(t ). Unlike traditional analyses on a nongrowing domain, where S ≡ 0, we show that domain growth leads to a very different situation where S can be positive. The theoretical tools developed and validated in this study allow us to distinguish between situations where the diffusive population reaches the moving boundary at x = L(t ) from other situations where the diffusive population never reaches the moving boundary at x = L(t ). Making this distinction is relevant to certain applications in developmental biology, such as the development of the enteric nervous system (ENS). All theoretical predictions are verified by implementing a discrete stochastic model.