Comparing methods for modelling spreading cell fronts

Spreading cell fronts play an essential role in many physiological processes. Classically, models of this process are based on the Fisher-Kolmogorov equation; however, such continuum representations are not always suitable as they do not explicitly represent behaviour at the level of individual cells. Additionally, many models examine only the large time asymptotic behaviour, where a travelling wave front with a constant speed has been established. Many experiments, such as a scratch assay, never display this asymptotic behaviour, and in these cases the transient behaviour must be taken into account. We examine the transient and asymptotic behaviour of moving cell fronts using techniques that go beyond the continuum approximation via a volume-excluding birth-migration process on a regular one-dimensional lattice. We approximate the averaged discrete results using three methods: (i) mean-field, (ii) pair-wise, and (iii) one-hole approximations. We discuss the performace of these methods, in comparison to the averaged discrete results, for a range of parameter space, examining both the transient and asymptotic behaviours. The one-hole approximation, based on techniques from statistical physics, is not capable of predicting transient behaviour but provides excellent agreement with the asymptotic behaviour of the averaged discrete results, provided that cells are proliferating fast enough relative to their rate of migration. The mean-field and pair-wise approximations give indistinguishable asymptotic results, which agree with the averaged discrete results when cells are migrating much more rapidly than they are proliferating. The pair-wise approximation performs better in the transient region than does the mean-field, despite having the same asymptotic behaviour. Our results show that each approximation only works in specific situations, thus we must be careful to use a suitable approximation for a given system, otherwise inaccurate predictions could be made.

Fragmentation pathways of 2,3-dimethyl-2,3-dinitrobutane cations in the gas phase

2,3-Dimethyl-2,3-dinitrobutane (DMNB) is an explosive taggant added to plastic explosives during manufacture making them more susceptible to vapour-phase detection systems. In this study, the formation and detection of gas-phase \[M+H](+), \[M+Li](+), \[M+NH(4)](+) and \[M+Na](+) adducts of DMNB was achieved using electrospray ionisation on a triple quadrupole mass spectrometer. The \[M+H](+) ion abundance was found to have a strong dependence on ion source temperature, decreasing markedly at source temperatures above 50 degrees C. In contrast, the \[M+Na](+) ion demonstrated increasing ion abundance at source temperatures up to 105 degrees C. The relative susceptibility of DMNB adduct ions toward dissociation was investigated by collision-induced dissociation. Probable structures of product ions and mechanisms for unimolecular dissociation have been inferred based on fragmentation patterns from tandem mass (MS/MS) spectra of source-formed ions of normal and isotopically labelled DMNB, and quantum chemical calculations. Both thermal and collisional activation studies suggest that the \[M+Na](+) adduct ions are significantly more stable toward dissociation than their protonated analogues and, as a consequence, the former provide attractive targets for detection by contemporary rapid screening methods such as desorption electrospray ionisation mass spectrometry. Copyright (C) 2009 Commonwealth of Australia. Published by John Wiley & Sons, Ltd.

Novel solutions for a model of wound healing angiogenesis

We prove the existence of novel, shock-fronted travelling wave solutions to a model of wound healing angiogenesis studied in Pettet et al (2000 IMA J. Math. App. Med. 17 395–413) assuming two conjectures hold. In the previous work, the authors showed that for certain parameter values, a heteroclinic orbit in the phase plane representing a smooth travelling wave solution exists. However, upon varying one of the parameters, the heteroclinic orbit was destroyed, or rather cut-off, by a wall of singularities in the phase plane. As a result, they concluded that under this parameter regime no travelling wave solutions existed. Using techniques from geometric singular perturbation theory and canard theory, we show that a travelling wave solution actually still exists for this parameter regime. We construct a heteroclinic orbit passing through the wall of singularities via a folded saddle canard point onto a repelling slow manifold. The orbit leaves this manifold via the fast dynamics and lands on the attracting slow manifold, finally connecting to its end state. This new travelling wave is no longer smooth but exhibits a sharp front or shock. Finally, we identify regions in parameter space where we expect that similar solutions exist. Moreover, we discuss the possibility of more exotic solutions.

Canards in advection-reaction-diffusion systems in one spatial dimension

This thesis contains a mathematical investigation of the existence of travelling wave solutions to singularly perturbed advection-reaction-diffusion models of biological processes. An enhanced mathematical understanding of these solutions and models is gained via the identification of canards (special solutions of fast/slow dynamical systems) and their role in the existence of the most biologically relevant, shock-like solutions. The analysis focuses on two existing models. A new proof of existence of a whole family of travelling waves is provided for a model describing malignant tumour invasion, while new solutions are identified for a model describing wound healing angiogenesis.

Influences of Allee effects in the spreading of malignant tumours

A recent study by Korolev et al. [Nat. Rev. Cancer, 14:371–379, 2014] evidences that the Allee effect—in its strong form, the requirement of a minimum density for cell growth—is important in the spreading of cancerous tumours. We present one of the first mathematical models of tumour invasion that incorporates the Allee effect. Based on analysis of the existence of travelling wave solutions to this model, we argue that it is an improvement on previous models of its kind. We show that, with the strong Allee effect, the model admits biologically relevant travelling wave solutions, with well-defined edges. Furthermore, we uncover an experimentally observed biphasic relationship between the invasion speed of the tumour and the background extracellular matrix density.

Role of electromechanical wave propagation in power systems

In a large interconnected power system, disturbances initiated by a fault or other events cause acceleration in the generator rotors with respect to their synchronous reference frame. This acceleration of rotors can be described by two different dynamic phenomena, as shown in existing literature. One of the phenomena is simultaneous acceleration and the other is electromechanical wave propagation, which is characterized by travelling waves in terms of a wave equation. This paper demonstrates that depending on the structure of the system, the exhibited dynamic response will be dominated by one phenomenon or the other or a mixture of both. Two system structures of choice are examined, with each structure exemplifying each phenomenon present to different degrees in their dynamic responses. Prediction of dominance of either dynamic phenomenon in a particular system can be determined by taking into account the relative sizes of the values of its reduced admittance matrix.

Ion-acoustic surface waves on a dielectric-dusty plasma interface

The theory of ion-acoustic surface wave propagation on the interface between a dusty plasma and a dielectric is presented. Both the constant and variable dust-charge cases are considered. It is found that massive negatively charged dust grains can significantly affect the propagation and damping of the surface waves. Application of the results to surface-wave generated plasmas is discussed. © 1998 IEEE.

Nonlinear ion-acoustic surface waves on a dusty plasma

A nonlinear theory for ion-acoustic surface waves propagating at the interface between a dusty plasma and a dielectric is presented. The nonlinear effects are associated with density modulations caused by surface-wave induced anomalous ionization. The negative charge of the massive dust grains is assumed to be constant. It is shown that the effect of the ionization nonlinearity arising from the ion-acoustic surface waves can result in the formation of surface envelope solitons. The wave phase shifts and the widths of the solitons are estimated for typical gas discharge plasmas.

Ion-acoustic waves in a complex plasma with negative ions

A wave propagation in a complex dusty plasma with negative ions was considered. The relevant processes such as ionization, electron attachment, diffusion, positive-negative ion recombination, plasma particle collisions, as well as elastic Coulomb and inelastic dust-charging collisions were taken self-consistently. It was found that the equilibrium of the plasma as well as the propagation of ion waves were modified to various degrees by these effects.

Characterization of acyl chain position in unsaturated phosphatidylcholines using differential mobility-mass spectrometry

Glycerophospholipids (GPs) that differ in the relative position of the two fatty acyl chains on the glycerol backbone (i.e., sn-positional isomers) can have distinct physicochemical properties. The unambiguous assignment of acyl chain position to an individual GP represents a significant analytical challenge. Here we describe a workflow where phosphatidylcholines (PCs) are subjected to ESI for characterization by a combination of differential mobility spectrometry and MS (DMS-MS). When infused as a mixture, ions formed from silver adduction of each phospholipid isomer {e.g., [PC (16:0/18:1) + Ag]+ and [PC (18:1/16:0) + Ag]+} are transmitted through the DMS device at discrete compensation voltages. Varying their relative amounts allows facile and unambiguous assignment of the sn-positions of the fatty acyl chains for each isomer. Integration of the well-resolved ion populations provides a rapid method (< 3 min) for relative quantification of these lipid isomers. The DMS-MS results show excellent agreement with established, but time-consuming, enzymatic approaches and also provide superior accuracy to methods that rely on MS alone. The advantages of this DMS-MS method in identification and quantification of GP isomer populations is demonstrated by direct analysis of complex biological extracts without any prior fractionation.

Electrochemistry of Heteroleptic Tris(phthalocyaninato) Rare Earth(III) Complexes

The electrochemical characteristics of a series of heteroleptic tris(phthalocyaninato) complexes with identical rare earths or mixed rare earths (Pc)M(OOPc)M(OOPc) [M = Eu...Lu, Y; H2Pc = unsubstituted phthalocyanine, H2(OOPc) = 3,4,12,13,21,22,30,31-octakis(octyloxy)phthalocyanine] and (Pc)Eu(OOPc)Er(OOPc) have been recorded and studied comparatively by cyclic voltammetry (CV) and differential pulse voltammetry (DPV) in CH2Cl2 containing 0.1 M tetrabutylammonium perchlorate (TBAP). Up to five quasi-reversible one-electron oxidations and four one-electron reductions have been revealed. The half-wave potentials of the first, second and fifth oxidations depend on the size of the metal center, but the fifth changes in the opposite direction to that of the first two. Moreover, the difference in redox potentials of the first oxidation and first reduction for (Pc)M(OOPc)M(OOPc), 0.85−0.98 V, also decreases linearly along with decreasing rare earth ion radius, clearly showing the rare earth ion size effect and indicating enhanced π−π interactions in the triple-deckers connected by smaller lanthanides. This order follows the red-shift seen in the lowest energy band of triple-decker compounds. The electronic differences between the lanthanides and yttrium are more apparent for triple-decker sandwich complexes than for the analogous double-deckers. By comparing triple-decker, double-decker and mononuclear [ZnII] complexes containing the OOPc ligand, the HOMO−LUMO gap has been shown to contract approximately linearly with the number of stacked phthalocyanine ligands.