No Arabic abstract
In this paper, we present a spatio-temporal mathematical model for simulating the formation and growth of a thrombus. Blood is treated as a multi-constituent mixture comprised of a linear fluid phase and a thrombus (solid) phase. The transport and reactions of 10 chemical and biological species are incorporated using a system of coupled convection-reaction-diffusion (CRD) equations to represent three processes in thrombus formation: initiation, propagation and stabilization. Computational fluid dynamic (CFD) simulations using the libraries of OpenFOAM were performed for two illustrative benchmark problems: in vivo thrombus growth in an injured blood vessel and in vitro thrombus deposition in micro-channels (1.5mm x 1.6mm x 0.1mm) with small crevices (125{mu}m x 75{mu}m and 125{mu}m x 137{mu}m). For both problems, the simulated thrombus deposition agreed very well with experimental observations, both spatially and temporally. Based on the success with these two benchmark problems, which have very different flow conditions and biological environments, we believe that the current model will provide useful insight into the genesis of thrombosis in blood-wetted devices, and provide a tool for the design of less thrombogenic devices.
Osteocytes and their cell processes reside in a large, interconnected network of voids pervading the mineralized bone matrix of most vertebrates. This osteocyte lacuno-canalicular network (OLCN) is believed to play important roles in mechanosensing, mineral homeostasis, and for the mechanical properties of bone. While the extracellular matrix structure of bone is extensively studied on ultrastructural and macroscopic scales, there is a lack of quantitative knowledge on how the cellular network is organized. Using a recently introduced imaging and quantification approach, we analyze the OLCN in different bone types from mouse and sheep that exhibit different degrees of structural organization not only of the cell network but also of the fibrous matrix deposited by the cells. We define a number of robust, quantitative measures that are derived from the theory of complex networks. These measures enable us to gain insights into how efficient the network is organized with regard to intercellular transport and communication. Our analysis shows that the cell network in regularly organized, slow-growing bone tissue from sheep is less connected, but more efficiently organized compared to irregular and fast-growing bone tissue from mice. On the level of statistical topological properties (edges per node, edge length and degree distribution), both network types are indistinguishable, highlighting that despite pronounced differences at the tissue level, the topological architecture of the osteocyte canalicular network at the subcellular level may be independent of species and bone type. Our results suggest a universal mechanism underlying the self-organization of individual cells into a large, interconnected network during bone formation and mineralization.
Bubbles introduced to the arterial circulation during invasive medical procedures can have devastating consequences for brain function but their effects are currently difficult to quantify. Here we present a Monte-Carlo simulation investigating the impact of gas bubbles on cerebral blood flow. For the first time, this model includes realistic adhesion forces, bubble deformation, fluid dynamical considerations, and bubble dissolution. This allows investigation of the effects of buoyancy, solubility, and blood pressure on embolus clearance. Our results illustrate that blockages depend on several factors, including the number and size distribution of incident emboli, dissolution time and blood pressure. We found it essential to model the deformation of bubbles to avoid overestimation of arterial obstruction. Incorporation of buoyancy effects within our model slightly reduced the overall level of obstruction but did not decrease embolus clearance times. We found that higher blood pressures generate lower levels of obstruction and improve embolus clearance. Finally, we demonstrate the effects of gas solubility and discuss potential clinical applications of the model.
Radiotherapy can effectively kill malignant cells, but the doses required to cure cancer patients may inflict severe collateral damage to adjacent healthy tissues. Hyperthermia (HT) is a promising option to improve the outcome of radiation treatment (RT) and is increasingly applied in hospital. However, the synergistic effect of simultaneous thermoradiotherapy is not well understood yet, while its mathematical modelling is essential for treatment planning. To better understand this synergy, we propose a theoretical model in which the thermal enhancement ratio (TER) is explained by the fraction of cells being radiosensitised by the infliction of sublethal damage through mild HT. Further damage finally kills the cell or inhibits its proliferation in a non-reversible process. We suggest the TER to be proportional to the energy invested in the sensitisation, which is modelled as a simple rate process. Assuming protein denaturation as the main driver of HT-induced sublethal damage and considering the temperature dependence of the heat capacity of cellular proteins, the sensitisation rates were found to depend exponentially on temperature; in agreement with previous empirical observations. Our predictions well reproduce experimental data from in-vitro and in-vivo studies, explaining the thermal modulation of cellular radioresponse for simultaneous thermoradiotherapy.
Severe Acute Respiratory Syndrome-CoronaVirus 2 (SARS-CoV2) caused the ongoing pandemic. This pandemic devastated the world by killing more than a million people, as of October 2020. It is imperative to understand the transmission dynamics of SARS-CoV2 so that novel and interdisciplinary prevention, diagnostic, and therapeutic techniques could be developed. In this work, we model and analyze the transmission of SARS-CoV2 through the human respiratory tract from a molecular communication perspective. We consider that virus diffusion occurs in the mucus layer so that the shape of the tract does not have a significant effect on the transmission. Hence, this model reduces the inherent complexity of the human respiratory system. We further provide the impulse response of SARS-CoV2-ACE2 receptor binding event to determine the proportion of the virus population reaching different regions of the respiratory tract. Our findings confirm the results in the experimental literature on higher mucus flow rate causing virus migration to the lower respiratory tract. These results are especially important to understand the effect of SARS-CoV2 on the different human populations at different ages who have different mucus flow rates and ACE2 receptor concentrations in the different regions of the respiratory tract.
Non-extensive statistical physics has allowed to generalize mathematical functions such as exponential and logarithms. The same framework is used to generalize sum and product so that the operations allow a more fluid way to work with mathematical expressions emerging from non-additive formulation of statistical physics. In this work we employ the generalization of the exponential, logarithm and product to obtain a formula for the survival fraction corresponding to the application of several radiation doses on a living tissue. Also we provide experimental recommendations to determine the universal characteristics of living tissues in interaction with radiation. These results have a potential application in radiobiology and radiation oncology.