Stochastic modelling of T cell homeostasis for two competing clonotypes via the master equation

Stochastic models for competing clonotypes of T cells by multivariate, continuous-time, discrete state, Markov processes have been proposed in the literature by Stirk, Molina-París and van den Berg (2008). A stochastic modelling framework is important because of rare events associated with small populations of some critical cell types. Usually, computational methods for these problems employ a trajectory-based approach, based on Monte Carlo simulation. This is partly because the complementary, probability density function (PDF) approaches can be expensive but here we describe some efficient PDF approaches by directly solving the governing equations, known as the Master Equation. These computations are made very efficient through an approximation of the state space by the Finite State Projection and through the use of Krylov subspace methods when evolving the matrix exponential. These computational methods allow us to explore the evolution of the PDFs associated with these stochastic models, and bimodal distributions arise in some parameter regimes. Time-dependent propensities naturally arise in immunological processes due to, for example, age-dependent effects. Incorporating time-dependent propensities into the framework of the Master Equation significantly complicates the corresponding computational methods but here we describe an efficient approach via Magnus formulas. Although this contribution focuses on the example of competing clonotypes, the general principles are relevant to multivariate Markov processes and provide fundamental techniques for computational immunology.

Analytical solutions of the multi-term modified power law wave equations in a finite domain

Fractional partial differential equations with more than one fractional derivative term in time, such as the Szabo wave equation, or the power law wave equation, describe important physical phenomena. However, studies of these multi-term time-space or time fractional wave equations are still under development. In this paper, multi-term modified power law wave equations in a finite domain are considered. The multi-term time fractional derivatives are defined in the Caputo sense, whose orders belong to the intervals (1, 2], [2, 3), [2, 4) or (0, n) (n > 2), respectively. Analytical solutions of the multi-term modified power law wave equations are derived. These new techniques are based on Luchko’s Theorem, a spectral representation of the Laplacian operator, a method of separating variables and fractional derivative techniques. Then these general methods are applied to the special cases of the Szabo wave equation and the power law wave equation. These methods and techniques can also be extended to other kinds of the multi term time-space fractional models including fractional Laplacian.

The time course of in vivo recovery of transverse strain in high-stress tendons following exercise

Objective To evaluate the time course of the recovery of transverse strain in the Achilles and patellar tendon following a bout of resistance exercise. Methods Seventeen healthy adults underwent sonographic examination of the right patellar (n=9) and Achilles (n=8) tendons immediately prior to and following 90 repetitions of weight-bearing quadriceps and gastrocnemius-resistance exercise performed against an effective resistance of 175% and 250% body weight, respectively. Sagittal tendon thickness was determined 20 mm from the enthesis and transverse strain, as defined by the stretch ratio, was repeatedly monitored over a 24 h recovery period. Results Resistance exercise resulted in an immediate decrease in Achilles (t7=10.6, p<0.01) and patellar (t8=8.9, p<0.01) tendon thickness, resulting in an average transverse stretch ratio of 0.86±0.04 and 0.82±0.05, which was not significantly different between tendons. The magnitude of the immediate transverse strain response, however, was reduced with advancing age (r=0.63, p<0.01). Recovery in transverse strain was prolonged compared with the duration of loading and exponential in nature. The average primary recovery time was not significantly different between the Achilles (6.5±3.2 h) and patellar (7.1±3.2 h) tendons. Body weight accounted for 62% and 64% of the variation in recovery time, respectively. Conclusions Despite structural and biochemical differences between the Achilles and patellar tendon, the mechanisms underlying transverse creep recovery in vivo appear similar and are highly time dependent. These novel findings have important implications concerning the time required for the mechanical recovery of high-stress tendons following an acute bout of exercise.

A simulation-based large deflection and inelastic analysis of steel frames under fire

This paper presents an accurate and robust geometric and material nonlinear formulation to predict structural behaviour of unprotected steel members at elevated temperatures. A fire analysis including large displacement effects for frame structures is presented. This finite element formulation of beam-column elements is based on the plastic hinge approach to model the elasto-plastic strain-hardening material behaviour. The Newton-Raphson method allowing for the thermal-time dependent effect was employed for the solution of the non-linear governing equations for large deflection in thermal history. A combined incremental and total formulation for determining member resistance is employed in this nonlinear solution procedure for the efficient modeling of nonlinear effects. Degradation of material strength with increasing temperature is simulated by a set of temperature-stress-strain curves according to both ECCS and BS5950 Part 8, which implicitly allows for creep deformation. The effects of uniform or non-uniform temperature distribution over the section of the structural steel member are also considered. Several numerical and experimental verifications are presented.

A linear-elastic model of anisotropic tumour growth

This paper examines the effect of anisotropic growth on the evolution of mechanical stresses in a linear-elastic model of a growing, avascular tumour. This represents an important improvement on previous linear-elastic models of tissue growth since it has been shown recently that spatially-varying isotropic growth of linear-elastic tissues does not afford the necessary stress-relaxation for a steady-state stress distribution upon reaching a nutrient-regulated equilibrium size. Time-dependent numerical solutions are developed using a Lax-Wendroff scheme, which show the evolution of the tissue stress distributions over a period of growth until a steady-state is reached. These results are compared with the steady-state solutions predicted by the model equations, and key parameters influencing these steady-state distributions are identified. Recommendations for further extensions and applications of this model are proposed.

Identification of eusynstyelamide B as a potent cell cycle inhibitor following the generation and screening of an ascidian-derived extract library using a real time cell analyzer

Ascidians are marine invertebrates that have been a source of numerous cytotoxic compounds. Of the first six marine-derived drugs that made anticancer clinical trials, three originated from ascidian specimens. In order to identify new anti-neoplastic compounds, an ascidian extract library (143 samples) was generated and screened in MDA-MB-231 breast cancer cells using a real-time cell analyzer (RTCA). This resulted in 143 time-dependent cell response profiles (TCRP), which are read-outs of changes to the growth rate, morphology, and adhesive characteristics of the cell culture. Twenty-one extracts affected the TCRP of MDA-MB-231 cells and were further investigated regarding toxicity and specificity, as well as their effects on cell morphology and cell cycle. The results of these studies were used to prioritize extracts for bioassay-guided fractionation, which led to the isolation of the previously identified marine natural product, eusynstyelamide B (1). This bis-indole alkaloid was shown to display an IC50 of 5 μM in MDA-MB-231 cells. Moreover, 1 caused a strong cell cycle arrest in G2/M and induced apoptosis after 72 h treatment, making this molecule an attractive candidate for further mechanism of action studies.

Time-independent charge carrier mobility in a model polymer:fullerene organic solar cell

The technique of photo-CELIV (charge extraction by linearly increasing voltage) is one of the more straightforward and popular approaches to measure the faster carrier mobility in measurement geometries that are relevant for operational solar cells and other optoelectronic devices. It has been used to demonstrate a time-dependent photocarrier mobility in pristine polymers, attributed to energetic relaxation within the density of states. Conversely, in solar cell blends, the presence or absence of such energetic relaxation on transport timescales remains under debate. We developed a complete numerical model and performed photo-CELIV experiments on the model high efficiency organic solar cell blend poly[3,6-dithiophene-2-yl-2,5-di(2-octyldodecyl)-pyrrolo[3,4-c]pyrrole-1,4-dione-alt-naphthalene] (PDPP-TNT):[6,6]-phenyl-C71-butyric-acid-methyl-ester (PC70BM). In the studied solar cells a constant, time-independent mobility on the scale relevant to charge extraction was observed, where thermalisation of photocarriers occurs on time scales much shorter than the transit time. Therefore, photocarrier relaxation effects are insignificant for charge transport in these efficient photovoltaic devices.

Unsteady flow and heat transfer between rotating coaxial disks

A study is made on the flow and heat transfer of a viscous fluid confined between two parallel disks. The disks are allowed to rotate with different time dependent angular velocities, and the upper disk is made to approach the lower one with a constant speed. Numerical solutions of the governing parabolic partial differential equations are obtained through a fourth-order accurate compact finite difference scheme. The normal forces and torques that the fluid exerts on the rotating surfaces are obtained at different nondimensional times for different values of the rate of squeezing and disk angular velocities. The temperature distribution and heat transfer are also investigated in the present analysis.

Tyrosine kinases and calcium dependent activation of endothelial cell phospholipase D by diperoxovanadate

Reactive oxygen species (ROS) mediated modulation of signal transduction pathways represent an important mechanism of cell injury and barrier dysfunction leading to the development of vascular disorders. Towards understanding the role of ROS in vascular dysfunction, we investigated the effect of diperoxovanadate (DPV), derived from mixing hydrogen peroxide and vanadate, on the activation of phospholipase D (PLD) in bovine pulmonary artery endothelial cells (BPAECs). Addition of DPV to BPAECs in the presence of .05% butanol resulted in an accumulation of [P-32] phosphatidylbutanol (PBt) in a dose- and time-dependent manner. DPV also caused an increase in tyrosine phosphorylation of several protein bands (Mr 20-200 kD), as determined by Western blot analysis with antiphosphotyrosine antibodies. The DPV-induced [P-32] PBt-accumulation was inhibited by putative tyrosine kinase inhibitors such as genistein, herbimycin, tyrphostin and by chelation of Ca2+ with either EGTA or BAPTA, however, pretreatment of BPAECs with the inhibitor PKC bisindolylmaleimide showed minimal inhibition. Also down-regulation of PKC alpha and epsilon, the major isotypes of PKC in BPAECs, by TPA (100 nM, 18 h) did not attenuate the DPV-induced PLD activation. The effects of putative tyrosine kinase and PKC inhibitors were specific as determined by comparing [P-32] PBt formation between DPV and TPA. In addition to tyrosine kinase inhibitors, antioxidants such as N-acetylcysteine and pyrrolidine dithiocarbamate also attenuated DPV-induced protein tyrosine phosphorylation and PLD stimulation. These results suggest that oxidation, prevented by reduction with thiol compounds, is involved in DPV-dependent protein tyrosine phosphorylation and PLD activation.

Unsteady flow over a stretching surface with a magnetic field in a rotating fluid

The effect of the magnetic field on the unsteady flow over a stretching surface in a rotating fluid has been studied. The unsteadiness in the flow field is due to the time-dependent variation of the velocity of the stretching surface and the angular velocity of the rotating fluid. The Navier-Stokes equations and the energy equation governing the flow and the heat transfer admit a self-similar solution if the velocity of the stretching surface and the angular velocity of the rotating fluid vary inversely as a linear function of time. The resulting system of ordinary differential equations is solved numerically using a shooting method. The rotation parameter causes flow reversal in the component of the velocity parallel to the strerching surface and the magnetic field tends to prevent or delay the flow reversal. The surface shear stresses dong the stretching surface and in the rotating direction increase with the rotation parameter, but the surface heat transfer decreases. On the other hand, the magnetic field increases the surface shear stress along the stretching surface, but reduces the surface shear stress in the rotating direction and the surface heat transfer. The effect of the unsteady parameter is more pronounced on the velocity profiles in the rotating direction and temperature profiles.

Unsteady free convection flow in the stagnation-point region of a rotating sphere

The unsteady free convection boundary-layer flow in the forward stagnation-point region of a sphere, which is rotating with time-dependent angular velocity in an ambient fluid, has been studied. Both constant wall temperature and constant hear flux conditions have been considered. The non-linear coupled parabolic partial differential equations governing the flow have been solved numerically using an implicit finite-difference scheme. The skin friction and the heat transfer are enhanced by the buoyancy force. The effect of the buoyancy force is found to be more pronounced for smaller Prandtl numbers than for larger Prandtl numbers. For a given buoyancy force, the heat transfer increases with an increase in Prandtl number, but the skin friction decreases.

Inertial mass of a vortex in cuprate superconductors

We present here a calculation of the inertial mass of a moving vortex in cuprate superconductors. This is a poorly known basic quantity of obvious interest in vortex dynamics. The motion of a vortex causes a dipolar density distortion and an associated electric field which is screened. The energy cost of the density distortion as well as the related screened electric field contributes to the vortex mass, which is small because of efficient screening. As a preliminary, we present a discussion and calculation of the vortex mass using a microscopically derivable phase-only action functional for the far region which shows that the contribution from the far region is negligible and that most of it arises from the (small) core region of the vortex. A calculation based on a phenomenological Ginzburg-Landau functional is performed in the core region. Unfortunately such a calculation is unreliable; the reasons for it are discussed. A credible calculation of the vortex mass thus requires a fully microscopic non-coarse-grained theory. This is developed, and results are presented for an s-wave BCS-like gap, with parameters appropriate to the cuprates. The mass, about 0.5m(e) per layer, for a magnetic field along the c axis arises from deformation of quasiparticle states bound in the core and screening effects mentioned above. We discuss earlier results, possible extensions to d-wave symmetry, and observability of effects dependent on the inertial mass. [S0163-1829(97)05534-3].

Simulation of the interaction between blood flow and atherosclerotic plaque

It has been well accepted that over 50% of cerebral ischemic events are the result of rupture of vulnerable carotid atheroma and subsequent thrombosis. Such strokes are potentially preventable by carotid interventions. Selection of patients for intervention is currently based on the severity of carotid luminal stenosis. It has been, however, widely accepted that luminal stenosis alone may not be an adequate predictor of risk. To evaluate the effects of degree of luminal stenosis and plaque morphology on plaque stability, we used a coupled nonlinear time-dependent model with flow-plaque interaction simulation to perform flow and stress/strain analysis for stenotic artery with a plaque. The Navier-Stokes equations in the Arbitrary Lagrangian-Eulerian (ALE) formulation were used as the governing equations for the fluid. The Ogden strain energy function was used for both the fibrous cap and the lipid pool. The plaque Principal stresses and flow conditions were calculated for every case when varying the fibrous cap thickness from 0.1 to 2mm and the degree of luminal stenosis from 10% to 90%. Severe stenosis led to high flow velocities and high shear stresses, but a low or even negative pressure at the throat of the stenosis. Higher degree of stenosis and thinner fibrous cap led to larger plaque stresses, and a 50% decrease of fibrous cap thickness resulted in a 200% increase of maximum stress. This model suggests that fibrous cap thickness is critically related to plaque vulnerability and that, even within presence of moderate stenosis, may play an important role in the future risk stratification of those patients when identified in vivo using high resolution MR imaging.

Unsteady incompressible boundary layer flow of a micropolar fluid at a stagnation point

The flow, heat and mass transfer on the unsteady laminar incompressible boundary layer in micropolar fluid at the stagnation point of a 2-dimensional and an axisymmetric body have been studied when the free stream velocity and the wall temperature vary arbitrarily with time. The partial defferential equations governing the flow have been solved numerically using a quasilinear finite-difference scheme. The skin friction, microrotation gradient and heat transfer parameters are found to be strongly dependent on the coupling parameter, mass transfer and time, whereas the effect of the microrotation parameter on the skin friction and heat transfer is rather weak, but microrotation gradient is strongly affected by it. The Prandtl number and the variation of the wall temperature with time affect the heat-transfer very significantly but the skin friction and micrortation gradient are unaffected by them.

A revisit to Pope's cohort analysis

Gulland's [Gulland, J.A., 1965. Estimation of mortality rates. Annex to Arctic Fisheries Working Group Report (meeting in Hamburg, January 1965). ICES. C.M. 1965, Doc. No. 3 (mimeographed)] virtual population analysis (VPA) is commonly used for studying the dynamics of harvested fish populations. However, it necessitates the solving of a nonlinear equation for the instantaneous rate of fishing mortality of the fish in a population. Pope [Pope, J.G., 1972. An investigation of the accuracy of Virtual Population Analysis using cohort analysis. ICNAF Res. Bull. 9, 65-74. Also available in D.H. Cushing (ed.) (1983), Key Papers on Fish Populations, p. 291-301, IRL Press, Oxford, 405 p.] eliminated this necessity in his cohort analysis by approximating its underlying age- and time-dependent population model. His approximation has since become one of the most commonly used age- and time-dependent fish population models in fisheries science. However, some of its properties are not well understood. For example, many assert that it describes the dynamics of a fish population, from which the catch of fish is taken instantaneously in the middle of the year. Such an assertion has never been proven, nor has its implied instantaneous rate of fishing mortality of the fish of a particular age at a particular time been examined, nor has its implied catch equation been derived from a general catch equation. In this paper, we prove this assertion, examine its implied instantaneous rate of fishing mortality of the fish of a particular age at a particular time, derive its implied catch equation from a general catch equation, and comment on how to structure an age- and time-dependent population model to ensure its internal consistency. This work shows that Gulland's (1965) virtual population analysis and Pope's (1972) cohort analysis lie at the opposite end of a continuous spectrum as a general model for a seasonally occurring fishery; Pope's (1972) approximation implies an infinitely large instantaneous rate of fishing mortality of the fish of a particular age at a particular time in a fishing season of zero length; and its implied catch equation has an undefined instantaneous rate of fishing mortality of the fish in a population, but a well-defined cumulative instantaneous rate of fishing mortality of the fish in the population. This work also highlights a need for a more careful treatment of the times of start and end of a fishing season in fish population models.