Dimension Estimate of Uniform Attractor for a Model of High Intensity Focussed Ultrasound-Induced Thermotherapy.

High intensity focussed ultrasound (HIFU) has emerged as a novel therapeutic modality, for the treatment of various cancers, that is gaining significant traction in clinical oncology. It is a cancer therapy that avoids many of the associated negative side effects of other more well-established therapies (such as surgery, chemotherapy and radiotherapy) and does not lead to the longer recuperation times necessary in these cases. The increasing interest in HIFU from biomedical researchers and clinicians has led to the development of a number of mathematical models to capture the effects of HIFU energy deposition in biological tissue. In this paper, we study the simplest such model that has been utilized by researchers to study temperature evolution under HIFU therapy. Although the model poses significant theoretical challenges, in earlier work, we were able to establish existence and uniqueness of solutions to this system of PDEs (see Efendiev et al. Adv Appl Math Sci 29(1):231-246, 2020). In the current work, we take the next natural step of studying the long-time dynamics of solutions to this model, in the case where the external forcing is quasi-periodic. In this case, we are able to prove the existence of uniform attractors to the corresponding evolutionary processes generated by our model and to estimate the Hausdorff dimension of the attractors, in terms of the physical parameters of the system.

Spatial invasion dynamics on random and unstructured meshes: implications for heterogeneous tumor populations.

In this work we discuss a spatial evolutionary model for a heterogeneous cancer cell population. We consider the gain-of-function mutations that not only change the fitness potential of the mutant phenotypes against normal background cells but may also increase the relative motility of the mutant cells. The spatial modeling is implemented as a stochastic evolutionary system on a structured grid (a lattice, with random neighborhoods, which is not necessarily bi-directional) or on a two-dimensional unstructured mesh, i.e. a bi-directional graph with random numbers of neighbors. We present a computational approach to investigate the fixation probability of mutants in these spatial models. Additionally, we examine the effect of the migration potential on the spatial dynamics of mutants on unstructured meshes. Our results suggest that the probability of fixation is negatively correlated with the width of the distribution of the neighborhood size. Also, the fixation probability increases given a migration potential for mutants. We find that the fixation probability (of advantaged, disadvantaged and neutral mutants) on unstructured meshes is relatively smaller than the corresponding results on regular grids. More importantly, in the case of neutral mutants the introduction of a migration potential has a critical effect on the fixation probability and increases this by orders of magnitude. Further, we examine the effect of boundaries and as intuitively expected, the fixation probability is smaller on the boundary of regular grids when compared to its value in the bulk. Based on these computational results, we speculate on possible better therapeutic strategies that may delay tumor progression to some extent.

Poroelastic theory of transcapillary flow: effects of endothelial glycocalyx deterioration.

The luminal surface of endothelial cells is lined with a carbohydrate-rich layer known as the endothelial glycocalyx. Identification of the structural properties of the glycocalyx has led investigators to examine its various functions and it has since been recognized as playing a role in many physiological processes, one of which is the regulation of fluid and protein exchange across the capillary wall. Experimental observations in which the glycocalyx was degraded in rat myocardial capillaries showed fluid accumulation in the tissue, suggesting that the glycocalyx acts as a protective barrier against edema. In this work we seek to quantify the observed edema formation by using our earlier poroelastic model to examine the consequences of glycocalyx deterioration on transcapillary filtration. Upon enzymatic treatment the properties of the glycocalyx, such as its thickness and permeability, will be modified, and our purpose here is to investigate quantitatively how changes in these parameters affect the magnitude of the fluid filtration through the capillary wall. We compare our results with both experimental data as well as other theoretical models where applicable, discussing the implications of the models as well as the limitations of comparison. This work provides the basis for further experiments that may better characterise many of the parameters involved in this process.

Dynamics of tumor growth and combination of anti-angiogenic and cytotoxic therapies.

Tumors cannot grow beyond a certain size (about 1-2 mm in diameter) through simple diffusion of oxygen and other essential nutrients into the tumor. Angiogenesis, the formation of blood vessels from pre-existing vessels, is a crucial and observed step, through which a tumor obtains its own blood supply. Thus, strategies that interfere with the development of this tumor vasculature, known as anti-angiogenic therapy, represent a novel approach to controlling tumor growth. Several pre-clinical studies have suggested that currently available angiogenesis inhibitors are unlikely to yield significant sustained improvements in tumor control on their own, but rather will need to be used in combination with conventional treatments to achieve maximal benefit. Optimal sequencing of anti-angiogenic treatment and radiotherapy or chemotherapy is essential to the success of these combined treatment strategies. Hence, a major challenge to mathematical modeling and computer simulations is to find appropriate dosages, schedules and sequencing of combination therapies to control or eliminate tumor growth. Here, we present a mathematical model that incorporates tumor cells and the vascular network, as well as their interplay. We can then include the effects of two different treatments, conventional cytotoxic therapy and anti-angiogenic therapy. The results are compared with available experimental and clinical data.

Mathematical modeling of brain tumors: effects of radiotherapy and chemotherapy.

Gliomas, the most common primary brain tumors, are diffusive and highly invasive. The standard treatment for brain tumors consists of a combination of surgery, radiation therapy and chemotherapy. Over the past few years, mathematical models have been applied to study untreated and treated brain tumors. In an effort to improve treatment strategies, we consider a simple spatio-temporal mathematical model, based on proliferation and diffusion, that incorporates the effects of radiotherapeutic and chemotherapeutic treatments. We study the effects of different schedules of radiation therapy, including fractionated and hyperfractionated external beam radiotherapy, using a generalized linear quadratic (LQ) model. The results are compared with published clinical data. We also discuss the results for combination therapy (radiotherapy plus temozolomide, a new chemotherapy agent), as proposed in recent clinical trials. We use the model to predict optimal sequencing of the postoperative (combination of radiotherapy and adjuvant, neo-adjuvant or concurrent chemotherapy) treatments for brain tumors.

The constitutive properties of the brain paraenchyma Part 2. Fractional derivative approach.

Fractional models have proven to be very useful for studying viscoelastic materials. We consider the fractional Zener model (also called four-parameter model) to study both the relaxation function and creep compliance. The analytical results are compared with the known experimental results of the human brain tissue to obtain the best fit and brain mechanical parameters. The results are also compared to the non-fractional Zener model and four-parameter Burgers model, indicating that the four-parameter fractional model gives a substantially better fit for the all experimental data.

The constitutive properties of the brain parenchyma Part 1. Strain energy approach.

In this paper we study several constitutive equations for the brain based on the strain energy density function. We use the polynomial function and hyper-elastic Ogden model for the strain energy and include the energy dissipation by a Prony series expansion. The models are compared with known unconfined compression experimental results of the human brain tissue to obtain the best fitted model and brain mechanical parameters. Finite element simulations are also performed using the given constitutive equations, and numerical solutions match the analytical results very closely. The results are compared with other analytical and numerical calculations.

Frequency dependence of complex moduli of brain tissue using a fractional Zener model.

Brain tissue exhibits viscoelastic behaviour. If loading times are substantially short, static tests are not sufficient to determine the complete viscoelastic behaviour of the material, and dynamic test methods are more appropriate. The concept of complex modulus of elasticity is a powerful tool for characterizing the frequency domain behaviour of viscoelastic materials. On the other hand, it is well known that classical viscoelastic models can be generalized by means of fractional calculus to describe more complex viscoelastic behaviour of materials. In this paper, the fractional Zener model is investigated in order to describe the dynamic behaviour of brain tissue. The model is fitted to experimental data of oscillatory shear tests of bovine brain tissue to verify its behaviour and to obtain the material parameters.

Modeling Invasion Dynamics with Spatial Random-Fitness Due to Micro-Environment.

Numerous experimental studies have demonstrated that the microenvironment is a key regulator influencing the proliferative and migrative potentials of species. Spatial and temporal disturbances lead to adverse and hazardous microenvironments for cellular systems that is reflected in the phenotypic heterogeneity within the system. In this paper, we study the effect of microenvironment on the invasive capability of species, or mutants, on structured grids (in particular, square lattices) under the influence of site-dependent random proliferation in addition to a migration potential. We discuss both continuous and discrete fitness distributions. Our results suggest that the invasion probability is negatively correlated with the variance of fitness distribution of mutants (for both advantageous and neutral mutants) in the absence of migration of both types of cells. A similar behaviour is observed even in the presence of a random fitness distribution of host cells in the system with neutral fitness rate. In the case of a bimodal distribution, we observe zero invasion probability until the system reaches a (specific) proportion of advantageous phenotypes. Also, we find that the migrative potential amplifies the invasion probability as the variance of fitness of mutants increases in the system, which is the exact opposite in the absence of migration. Our computational framework captures the harsh microenvironmental conditions through quenched random fitness distributions and migration of cells, and our analysis shows that they play an important role in the invasion dynamics of several biological systems such as bacterial micro-habitats, epithelial dysplasia, and metastasis. We believe that our results may lead to more experimental studies, which can in turn provide further insights into the role and impact of heterogeneous environments on invasion dynamics.

Brain biomechanics: mathematical modeling of hydrocephalus.

The considerable amount of literature on mathematical models of hydrocephalus and other brain abnormalities is critically reviewed. These models have various degrees of mathematical sophistication, and have influenced not only the diagnosis of hydrocephalus, but also its treatment with CSF shunts. The mathematical models are classified into two classes, pressure-volume models, and consolidation models. Advantages and disadvantages of both types are pointed out with a view to removing the confusion frequently generated by the technical aspects of the subject. The conclusion is reached that, while none of the current models are good enough to be of immediate use to the neurosurgeon, mathematical models are likely in the future to be a powerful tool for the understanding and the treatment of hydrocephalus, as well as other conditions related to brain biomechanics. The amount of mathematics has been kept to the absolute minimum, but it is cited and appended for those who would like to dig further into this fascinating area of research.

Medical evidence for the defence in homicide trials.

The legal system for the investigation of major crime appears to be rather loose, in that there are no specific statutory provisions available to a defendant to present his or her case properly and adequately in court. There are no guidelines issued to the coroners nor are there any legal powers available to them to facilitate and authorize post-mortem examinations on behalf of the defence. This paper focuses attention on this legal dilemma, and attempts to suggest possible remedies.

Fabry's disease--a rare cause of sudden death.

A case of Fabry's disease, an X-linked in-born error of metabolism, presenting as a sudden death in an elderly female, is reported, along with a review of the literature. The pathological findings of the cellular vacuolation and the characteristic inclusion bodies in the tissues are described.

Traumatic sub-arachnoid haemorrhage as part of the NAI syndrome.

Case of fatal traumatic subarachnoid haemorrhage, as part of non-accidental injury in children is reported. The possible mechanisms are discussed.

Decapitation--a rare complication in hanging.

Decapitation from hanging is rare and there appears to be a lack of reports in scientific literature. Two cases of this phenomenon are reported.

The Bradford Fire Disaster. Part 1. The initial investigations: who died, where and how?

In May 1985 a football ground fire resulted in 50 deaths. All the bodies were severely burned. This paper provides an account of the planning and methods used in the initial investigation and assesses their value and the problems encountered.

The Bradford Fire Disaster. Part 2. Accident reconstruction: who died, when and why.

The circumstances of the Bradford Football Stadium Fire have been described in Part 1. Part 2 relates the positions in which the bodies were found, their carbon monoxide concentrations, and the physical configuration of the stand to the seat of the fire. It also examines the causes of the high death toll and makes recommendations for the future planning and management of sports stadiums.

Death from diesel fumes.

A rare fatality from inhalation of diesel motor exhaust fumes is reported and the toxicity of diesel fumes is discussed briefly.

Efficacy of dose escalation on TCP, recurrence and second cancer risks: a mathematical study.

OBJECTIVE: We investigated the effects of conventional and hypofractionation protocols by modelling tumour control probability (TCP) and tumour recurrence time, and examined their impact on second cancer risks. The main objectives of this study include the following: (a) incorporate tumour recurrence time and second cancer risks into the TCP framework and analyse the effects of variable doses and (b) investigate an efficient protocol to reduce the risk of a secondary malignancy while maximizing disease-free survival and tumour control. METHODS: A generalized mathematical formalism was developed that incorporated recurrence and second cancer risk models into the TCP dynamics. RESULTS: Our results suggest that TCP and relapse time are almost identical for conventional and hypofractionated regimens; however, second cancer risks resulting from hypofractionation were reduced by 22% when compared with the second cancer risk associated with a conventional protocol. The hypofractionated regimen appears to be sensitive to dose escalation and the corresponding impact on tumour recurrence time and reduction in second cancer risks. The reduction in second cancer risks is approximately 20% when the dose is increased from 60 to 72 Gy in a hypofractionated protocol. CONCLUSION: Our results suggest that hypofractionation may be a more efficient regimen in the context of TCP, relapse time and second cancer risks. Overall, our study demonstrates the importance of including a second cancer risk model in designing an efficient radiation regimen. ADVANCES IN KNOWLEDGE: The impact of various fractionation protocols on TCP and relapse in conjunction with second cancer risks is an important clinical question that is as yet unexplored.

Treatment of distal radial fractures with the DVR-A plate--the early Bristol experience.

AIM: The aim of the study was to assess how adequately distal radial fracture reduction was reproduced and maintained with the distal volar radius anatomic DVR-A (Biomet, Inc) locking plate. METHODS: We looked at a consecutive series of 111 patients treated with the DVR-A plate at our institution from 2007-2010. The preoperative, intra-operative, and postoperative films were reviewed. The AO and Frykman classification was recorded. The sagittal tilt, radial inclination, and radial length were measured on intra-operative X-rays and compared with final follow-up X-rays. The last recorded range of motion at follow-up and a functional assessment using the Disabilities of the Arm, Shoulder, and Hand (DASH) score was recorded. RESULTS: At final follow-up a mean radial inclination of 22.22 degrees (11 to 38 degrees), radial height of 11.85 mm (6 to 18 mm), and sagittal tilt of 6.71 degrees volar (-9 to 19 degrees) was achieved. From initial post-operative films, to final follow-up X-rays there was a mean increase of 0.17 degrees in radial inclination, a mean loss of 0.36 mm radial height, and a mean loss of 2.17 degrees volar tilt. The mean extension was 46.8 degrees, flexion 48.3 degrees, pronation 77.4 degrees, supination 74.8 degrees, radial deviation 15.3 degrees and the ulnar deviation 19 degrees. The mean DASH score was 12.8 (0-68). CONCLUSION: The DVR-A plate achieved a highly satisfactory reduction of radial length and radial inclination, with a small loss of volar sagittal tilt at final follow-up. A good functional outcome was reported, with a satisfactory range of motion achieved. The DVR-A plate is a safe and effective treatment for unstable and intra-articular displaced distal radius fractures, particularly in younger patients, in the short term.

Aging impact on brain biomechanics with applications to hydrocephalus.

Hydrocephalus is a neurological disorder whose clinical symptoms and treatment outcome are correlated with patient age. In Wilkie et al. (2010, A theoretical study of the effect of intraventricular pulsations on the pathogenesis of hydrocephalus. Appl. Math. Comput., 215, 3181-3191), the fractional Zener model was used to investigate the role of cerebrospinal fluid pressure pulsations in the development of hydrocephalus in infants and adults. In this paper, we determine the mechanical parameters of the fractional Zener model for the infant and adult brains using age-dependent shear complex modulus data (Thibault, K. L. & Margulies, S. S. (1998) Age-dependent material properties of the porcine cerebrum: effect on pediatric inertial head injury criteria. J. Biomech., 31, 1119-1126). The displacement of brain tissue under conditions representing the onset of hydrocephalus are then calculated. The infant brain was found to produce tissue displacements that are unphysical for our model geometry and a new boundary condition is proposed to replace the stress-free outer boundary condition used in Wilkie et al. (2010). The steadystate elastic modulus is identified as the parameter of interest in the development of hydrocephalus: it is found to increase from the infant value of 621 Pa to the young adult value of 955 Pa and we hypothesize that it then decreases with age. The low steady-state elastic modulus of the infant brain (and possibly the aged brain) increases the tissue's susceptibility to large deformations and thus to the ventricular expansion characteristic of hydrocephalus.