Search

The human gut microbiota is a complex ecosystem that affects a range of human physiology. In order to explore the dynamics of the human gut microbiota, we used a system of ordinary differential equations to model mathematically the biomass of three microorganism populations: $\textit{Bacteroides thetaiotaomicron}$, $\textit{Eubacterium rectale}$, and $\textit{Methanobrevibacter smithii}$. Additionally, we modeled the concentrations of relevant nutrients necessary to sustain these populations over time. Our model highlights the interactions and the competition among these three species. These three microorganisms were specifically chosen due to the system's end product, butyrate, which is a short chain fatty acid that aids in developing and maintaining the intestinal barrier in the human gut. The basis of our mathematical model assumes the gut is structured such that bacteria and nutrients exit the gut at a rate proportional to its volume, the rate of volumetric flow, and the biomass or concentration of the particular population or nutrient. We performed global sensitivity analyses using Sobol' sensitivities to estimate the relative importance of model parameters on simulation results.

Mounting evidence for the role of oxidative stress in the degeneration of articular cartilage after an injurious impact requires our modeling & simulation efforts to temporarily shift from just describing the effect of mechanical stress and inflammation on osteoarthritis (OA). The hypothesis that the injurious impact causes irreversible damage to chondrocyte mitochondria, which in turn increase their production of free radicals, affecting their energy production and their ability to rebuild the extracellular matrix, has to be modeled and the processes quantified in order to further the understanding of OA, its causes, and viable treatment options. The current article presents a calibrated model that captures the damage oxidative stress incurs on the cell viability, ATP production, and cartilage stability in a cartilage explant after a drop-tower impact. The model validates the biological hypothesis and will be used in further efforts to identify possibilities for treatment and be a part of a bigger modeling & simulation framework for the development of OA.

Multiple myeloma (MM), a plasma cell cancer, is associated with many health challenges, including damage to the kidney by tubulointerstitial fibrosis. We develop a mathematical model which captures the qualitative behavior of the cell and protein populations involved. Specifically, we model the interaction between cells in the proximal tubule of the kidney, free light chains, renal fibroblasts, and myeloma cells. We analyze the model for steady-state solutions to find a mathematically and biologically relevant stable steady-state solution. This foundational model provides a representation of dynamics between key populations in tubulointerstitial fibrosis that demonstrates how these populations interact to affect patient prognosis in patients with MM and renal impairment., Comment: Included version of model without tumor with steady-state analysis, corrected equations for free light chains and renal fibroblasts in model with tumor to reflect steady-state analysis, updated abstract, updated and added references

A severe application of stress on articular cartilage can initiate a cascade of biochemical reactions that can lead to the development of osteoarthritis. We constructed a multiscale mathematical model of the process with three components: cellular, chemical, and mechanical. The cellular component describes the different chondrocyte states according to the chemicals these cells release. The chemical component models the change in concentrations of those chemicals. The mechanical component contains a simulation of pressure application onto a cartilage explant and the resulting strains that initiate the biochemical processes. The model creates a framework for incorporating explicit mechanics, simulated by finite element analysis, into a theoretical biology framework.

Multiple myeloma is a plasma cell cancer that leads to a dysregulated bone remodeling process. We present a partial differential equation model describing the dynamics of bone remodeling with the presence of myeloma tumor cells. The model explicitly takes into account the roles of osteoclasts, osteoblasts, precursor cells, stromal cells, osteocytes, and tumor cells. Previous models based on ordinary differential equations make the simplifying assumption that the bone and tumor cells are adjacent to each other. However, in actuality, these cell populations are separated by the bone marrow. Our model takes this separation into account by including the diffusion of chemical factors across the marrow, which can be viewed as communication between the tumor and bone. Additionally, this model incorporates the growth of the tumor and the diminishing bone mass by utilizing a ``moving boundary.'' We present numerical simulations that qualitatively validate our model's description of the cell population dynamics.

We present a model of articular cartilage lesion formation to simulate the effects of cyclic loading. This model extends and modifies the reaction-diffusion-delay model by Graham et al. 2012 for the spread of a lesion formed though a single traumatic event. Our model represents "implicitly" the effects of loading, meaning through a cyclic sink term in the equations for live cells. Our model forms the basis for in silico studies of cartilage damage relevant to questions in osteoarthritis, for example, that may not be easily answered through in vivo or in vitro studies. Computational results are presented that indicate the impact of differing levels of EPO on articular cartilage lesion abatement.

Until recently many studies of bone remodeling at the cellular level have focused on the behavior of mature osteoblasts and osteoclasts, and their respective precursor cells, with the role of osteocytes and bone lining cells left largely unexplored. This is particularly true with respect to the mathematical modeling of bone remodeling. However, there is increasing evidence that osteocytes play important roles in the cycle of targeted bone remodeling, in serving as a significant source of RANKL to support osteoclastogenesis, and in secreting the bone formation inhibitor sclerostin. Moreover, there is also increasing interest in sclerostin, an osteocyte-secreted bone formation inhibitor, and its role in regulating local response to changes in the bone microenvironment. Here we develop a cell population model of bone remodeling that includes the role of osteocytes, sclerostin, and allows for the possibility of RANKL expression by osteocyte cell populations. This model extends and complements many of the existing mathematical models for bone remodeling but can be used to explore aspects of the process of bone remodeling that were previously beyond the scope of prior modeling work. Through numerical simulations we demonstrate that our model can be used to theoretically explore many of the most recent experimental results for bone remodeling, and can be utilized to assess the effects of novel bone-targeting agents on the bone remodeling process.

In this paper we develop mathematical models for collective cell motility. Initially we develop a model using a linear diffusion-advection type equation and fit the parameters to data from cell motility assays. This approach is helpful in classifying the results of cell motility assay experiments. In particular, this model can determine degrees of directed versus undirected collective cell motility. Next we develop a model using a nonlinear diffusion term that is able capture in a unified way directed and undirected collective cell motility. Finally we apply the nonlinear diffusion approach to a problem in tumor cell invasion, noting that neither chemotaxis or haptotaxis are present in the system under consideration in this article.

Injuries to articular cartilage result in the development of lesions that form on the surface of the cartilage. Such lesions are associated with articular cartilage degeneration and osteoarthritis. The typical injury response often causes collateral damage, primarily an effect of inflammation, which results in the spread of lesions beyond the region where the initial injury occurs. We present a minimal mathematical model based on known mechanisms to investigate the spread and abatement of such lesions. In particular we represent the "balancing act" between pro-inflammatory and anti-inflammatory cytokines that is hypothesized to be a principal mechanism in the expansion properties of cartilage damage during the typical injury response. We present preliminary results of in vitro studies that confirm the anti-inflammatory activities of the cytokine erythropoietin (EPO). We assume that the diffusion of cytokines determine the spatial behaviour of injury response and lesion expansion so that a reaction-diffusion system involving chemical species and chondrocyte cell state population densities is a natural way to represent cartilage injury response. We present computational results using the mathematical model showing that our representation is successful in capturing much of the interesting spatial behaviour of injury associated lesion development and abatement in articular cartilage. Further, we discuss the use of this model to study the possibility of using EPO as a therapy for reducing the amount of inflammation induced collateral damage to cartilage during the typical injury response. The mathematical model presented herein suggests that not only are anti-inflammatory cytokines, such as EPO necessary to prevent chondrocytes signaled by pro-inflammatory cytokines from entering apoptosis, they may also influence how chondrocytes respond to signaling by pro-inflammatory cytokines., Comment: 14 pages, 6 figures

Irregular bone remodeling is associated with a number of bone diseases such as osteoporosis and multiple myeloma. Computational and mathematical modeling can aid in therapy and treatment as well as understanding fundamental biology. Different approaches to modeling give insight into different aspects of a phenomena so it is useful to have an arsenal of various computational and mathematical models. Here we develop a mathematical representation of bone remodeling that can effectively describe many aspects of the complicated geometries and spatial behavior observed. There is a sharp interface between bone and marrow regions. Also the surface of bone moves in and out, i.e. in the normal direction, due to remodeling. Based on these observations we employ the use of a level-set function to represent the spatial behavior of remodeling. We elaborate on a temporal model for osteoclast and osteoblast population dynamics to determine the change in bone mass which influences how the interface between bone and marrow changes. We exhibit simulations based on our computational model that show the motion of the interface between bone and marrow as a consequence of bone remodeling. The simulations show that it is possible to capture spatial behavior of bone remodeling in complicated geometries as they occur \emph{in vitro} and \emph{in vivo}. By employing the level set approach it is possible to develop computational and mathematical representations of the spatial behavior of bone remodeling. By including in this formalism further details, such as more complex cytokine interactions and accurate parameter values, it is possible to obtain simulations of phenomena related to bone remodeling with spatial behavior much as \emph{in vitro} and \emph{in vivo}. This makes it possible to perform \emph{in silica} experiments more closely resembling experimental observations., Comment: Math. Biosci. Eng., 9(2), 2012

We present models of dormancy in a planktonic culture and in biofilm, and examine the relative advantage of short dormancy versus long dormancy times in each case. Simulations and analyses indicate that in planktonic batch cultures and in chemostats, live biomass is maximized by the fastest possible exit from dormancy. The lower limit of time to reawakening is thus perhaps governed by physiological, biochemical or other constraints within the cells. In biofilm we see that the slower waker has a defensive advantage over the fast waker due to a larger amount of dormant biomass, without an appreciable difference in total live biomass. Thus it would seem that typical laboratory culture conditions can be unrepresentative of the natural state. We discuss the computational methods developed for this work.

In this paper we consider chemotherapy in a spatial model of tumor growth. The model, which is of reaction-diffusion type, takes into account the complex interactions between the tumor and surrounding stromal cells by including densities of endothelial cells and the extra-cellular matrix. When no treatment is applied the model reproduces the typical dynamics of early tumor growth. The initially avascular tumor reaches a diffusion limited size of the order of millimeters and initiates angiogenesis through the release of vascular endothelial growth factor (VEGF) secreted by hypoxic cells in the core of the tumor. This stimulates endothelial cells to migrate towards the tumor and establishes a nutrient supply sufficient for sustained invasion. To this model we apply cytostatic treatment in the form of a VEGF-inhibitor, which reduces the proliferation and chemotaxis of endothelial cells. This treatment has the capability to reduce tumor mass, but more importantly, we were able to determine that inhibition of endothelial cell proliferation is the more important of the two cellular functions targeted by the drug. Further, we considered the application of a cytotoxic drug that targets proliferating tumor cells. The drug was treated as a diffusible substance entering the tissue from the blood vessels. Our results show that depending on the characteristics of the drug it can either reduce the tumor mass significantly or in fact accelerate the growth rate of the tumor. This result seems to be due to complicated interplay between the stromal and tumor cell types and highlights the importance of considering chemotherapy in a spatial context., Comment: revised version, 25 pages, 9 figures, minor misprints corrected

We derive a physiologically structured multiscale model for biofilm development. The model has components on two spatial scales, which induce different time scales into the problem. The macroscopic behavior of the system is modeled using growth-induced flow in a domain with a moving boundary. Cell-level processes are incorporated into the model using a so-called physiologically structured variable to represent cell senescence, which in turn affects cell division and mortality. We present computational results for our models which shed light on modeling the combined role senescence and the biofilm state play in the defense strategy of bacteria.

We present models and computational results which indicate that the spatial and temporal regularity seen in Proteus mirabilis swarm-colony development is largely an expression of a sharp age of dedifferentiation in the cell cycle from motile swarmer cells to immotile dividing cells (also called swimmer or vegetative cells.) This contrasts strongly with reaction-diffusion models of Proteus behavior that ignore or average out the age structure of the cell population and instead use only density-dependent mechanisms. We argue the necessity of retaining the explicit age structure, and suggest experiments that may help determine the underlying mechanisms empirically. Consequently, we advocate Proteus as a model organism for a multiscale understanding of how and to what extent the life cycle of individual cells affects the macroscopic behavior of a biological system.

We present multiscale models of cancer tumor invasion with components at the molecular, cellular, and tissue levels. We provide biological justifications for the model components, present computational results from the model, and discuss the scientific-computing methodology used to solve the model equations. The models and methodology presented in this paper form the basis for developing and treating increasingly complex, mechanistic models of tumor invasion that will be more predictive and less phenomenological. Because many of the features of the cancer models, such as taxis, aging and growth, are seen in other biological systems, the models and methods discussed here also provide a template for handling a broader range of biological problems.

In this paper we present continuous age- and space-structured models and numerical computations of Proteus mirabilis swarm-colony development. We base the mathematical representation of the cell-cycle dynamics of Proteus mirabilis on those developed by Esipov and Shapiro, which are the best understood aspects of the system, and we make minimum assumptions about less-understood mechanisms, such as precise forms of the spatial diffusion. The models in this paper have explicit age-structure and, when solved numerically, display both the temporal and spatial regularity seen in experiments, whereas the Esipov and Shapiro model, when solved accurately, shows only the temporal regularity. The composite hyperbolic-parabolic partial differential equations used to model Proteus mirabilis swarm-colony development are relevant to other biological systems where the spatial dynamics depend on local physiological structure. We use computational methods designed for such systems, with known convergence properties, to obtain the numerical results presented in this paper.