
MBW:Modeling of Oxygen Transport Across Tumor Multicellular LayersFrom MathBioArticle review by Kumar Thurimella & Nikolas Kauffman ContentsArticleOur paper by Rod Braun and Alexis Beatty, discusses oxygen and its role in tumor cells ^{[1]}. Executive SummaryIn this article the authors were trying to understand how different levels of oxygen within tumors affect the response of those tumors with different treatments. They furthered their question by modeling the oxygen diffusion across a lab grown melanoma cell using Fick's diffusion equation. The authors had initially sought out to figure out any correlation between a cancer patients outcome with oxygen transport in tumor parenchyma. Each tumor was examined on a multicellular level and thus the authors tested different multicellular levels (MCL) of a certain tumor. At each MCL there was a measured oxygen tension and the authors fit those profiles to a ondimensional diffusion model. The model was justified with experimental data that was taken from the experiment to describe the oxygen tension at each of those multicellular levels. This model is able to analyze certain tumor cell lines in MCL that would not normally be analyzed in the more common, spheroid tumor. The MCL can be studied in this manner since, unlike spheroids, they grow on an underlying substrate. The authors discuss the importance of oxygen transport, specifically its impact on MCL and how it can affect the treatment of cancer. Context/Biological Events
HistoryThe notion behind studying oxygen levels and the correlation to tumor cells has been studied previously by researchers, W. MuellerKlieser, J. P. Freyer, and R. M. Sutherland. ^{[3]} This history extends to the notion to believe that tumor cells are dependent on oxygen and there is a correlation between a cell line survival with oxygen access. This topic has examined the idea that tumor spheroids (like those in human bodies) adapt their metabolism based on oxygen and glucose intake. At that time the results had shown to be inconclusive and the authors consensus was that: although the metabolism can vary with glucose and oxygen intake, there are too many other possibilities to simply state that growth is solely dependent on glucose and oxygen. When the researchers modeled their oxygen tension profiles in the tumor spheroids they were fitting experimental data based on oxygen diffusion and plotting against the diameter of the spheroid. There was no model presented formally and the ideas and notions of Fick's diffusion equation was an underlying part of the experiment. Fick's equation in our setting would relate the diffusive flux between the outside and the inside of the spheroid with the concentration of the oxygen in the environment. These concepts lay a foreground into the insight of this interesting research and is able to help start a connection between oxygen levels and tumor growth. Mathematical ModelParameters of the Model
To model the oxygen transport across OCM1 cultured MCL's we must use Fick's onedimensional oxygen diffusion equation: From Fick's diffusion equation and assumptions made for the method of experimentation that the change in oxygen tension per unit time (dP/dt) is equal to zero, i.e. the steadystate, we obtain a new equation that models oxygen tension with respect to distance x as follows: where is the oxygen consumption across the MCL. In order to adequately determine a bestfitted model for the oxygen diffusion across the MCL, the authors decided to fit four different diffusion models to the experimental data. These models only differed with respect to the number of boundary layers between the media layer and membrane layer of the MCL. By solving the secondorder inhomogenous ODE using the method of undetermined coefficients the following equation was obtained: The law of mass conservation tells us that the and oxygen consumption levels at the boundary of adjacent layers must be equal, yielding: The constants and are determined by the boundary conditions of the MCL which are defined as follows: To compare oxygen consumption rates across the MCL it was necessary to define these rates according to their changes across different boundary regimes. The rates were taken at the media interface, membrane interface, and across the MCL. These were calculated using the following relations: The latter two of the above three equations involve xvalues which are defined as the minimum values from the media and membrane interfaces where . The oxygen consumption rates for these two equations are thus calculated as the average rate from the boundary layer to these minimum xvalues. As stated earlier, four different diffusion models were used to determine the bestfit model across the MCL. The accuracy of the fit was determined using MATLAB written programs that evaluated the sumofsquares error, rootmeansquared error, and the coefficient of determination . Understanding the Mathematical ModelFrom Fick's diffusion model in the previous section the assumption that was made, thus, the model with respect to time is in steadystate. This is understood through the realization that upon experimentation the tumor MCL's have already assumed a uniform oxygen profile and therefore, oxygen levels across the cells only differ with respect to the x location within the MCL. Solving the steadystate then tells us that the secondderivative of the oxygen levels with respect to position equal the oxygen consumption rate across the entire MCL. The distance across the MCL was then broken into different layers as to obtain different profiles. This allowed for multiple profiles to be examined against the experimental data, which, allowed for the model to better represent the behavior of the test results. The secondorder ODE produced by finding the steadystate of the Fick's onedimensional diffusion model made it possible to then solve for the levels across each layer using the characteristic equation and then using the principle of superpostion to additively combine the homogeneous and particular solutions of the ODE. For a quick refresh, lets solve for the roots of the characteristic equation . This models C.E. parameters are a=1, b=0, c=0; yielding a repeated root of 0 by solving the quadratic equation. Since we have a repeated root solutions to the homogeneous ODE are not linearly independent and are made so by multiplying one of the terms in the linear combination by x. By setting r=0 and combining the particular solution (i.e. the oxygen consumption rate), we arrive at the solution: This defines the across each layer within the MCL. ResultsA. Measuring Profiles Across OCM1 MCL profiles were measured using polarographic microelectrodes in MCL's ranging from 400 to 1400 micrometers. These measurements demonstrated that thicker MCL's gave bigger regions of anoxia across the MCL. B. Diffusion Modeling of the Profiles Depending on the extent of anoxia across the MCL it was found that two of the four models used to describe the profiles most accurately displayed the behavior of the experimental data. For MCL's with insignificant anoxic regions the tworegion model was found to be the best representation for the data, whereas; thicker MCL's with longer anoxic regions were found to be best represented by the model consisting of one extra layer (the central layer) where the levels were set to be zero, i.e. the threeregion model. Upon reexamination of the mathematical model, the tworegion representation requires five unknown parameter fits for and the threeregion model requires seven unknown parameter fits for . Because the oneregion model does not accurately represent the data it can be said that the boundaries of the MCL at the membrane and media are metabolically different, i.e. the cells near these boundaries are consuming oxygen at different rates. C. Oxygen Consumption of OCM1 MCL It was found that the oxygen consumption rates of the cancerous cells at the media and membrane boundaries were different, indicating different metabolic rates at the two interfaces. The metabolic rate at the membrane surface was found to be higher than that of the free surface (lower values of at the MCL membrane boundary). Unfortunately, it is not known why this phenomena occurs, but the authors speculate that the metabolism levels of these cells close to the boundary could be affected by the collagen matrix interacting with the membrane interface. Another important characteristic of the oxygen consumption is that the average consumption decreased with increasing MCL thickness. Based upon correlations from earlier, it was seen that increasing MCL thickness indicated longer anoxic regions in the profiles. Trivial examination of this demonstrates that the inclusion of zero oxygen consumption for these anoxic regions decreased the average oxygen consumption, . There were also differences in the oxygen consumption parameters when the MCL was exposed to air saturated RPMI as compared to 5% saturated RPMI. In both instances, it was evident that more oxygen was being supplied at the free surface, but in the 5% saturated RPMI test procedure it was shown that all of the oxygen consumption parameters decreased in value. ConclusionsMeasurement of profiles The authors had measured various profiles between a thickness of 400 to 1400 and their finding in this area led them to look at anoxia at those various profiles. The authors generalized that the more thick the MCL the less oxygen (or the higher anoxia). This finding warrants itself to the idea that with very a thick cell layer there is not very much oxygen that can interact with the cells. However, that being said there is no quantitative element in describing at what specific thickness oxygen is completely cutoff or is very prevalent in. The differences in the profile within the same MCL region is thought to be cause by the various metabolic needs in several regions of the MCL. Figure one describes the general correlation between the thickness level of the MCL and the anoxia (lack of oxygen) there is around that region. Comparison of Oxygen Measurements in Spheroids vs. MCL The authors stand behind their claim that in order to properly study oxygen uptake in tumor cells using an MCL model is better than a spheroid model. For one spheroids are difficult to measure under. When measuring oxygen tension in a spheroid it becomes very difficult because everything has to be measured radially which is technically challenging and requires quite a bit of skilled manipulation. However, an MCL is grown perpendicularly with relation to a flat surface and there is much less positional error. There are some disadvantages to using an MCL model, however. Even within an MCL model not all cell lines can be grown which limits the ability of tumor cells to be studied. Since not all cell lines can be grown, oxygen consumption levels cannot even be determined which also limits sides of this experiments. The MCL model still opened up more possibility for study in the oxygen levels within tumor cells. Diffusion Modeling of Profiles Throughout this paper several region models were introduced and were tried in order to fit the experimental data. When there was no anoxia (presence of oxygen) the authors found that a 2 region model best described the data. The 1 region model was not able to fit data at different boundaries of the MCL but the 2 region model could account for the metabolic differences among the various boundaries. The authors found that cells growing near the membrane of the MCL had a higher metabolism even though the oxygen tension was lower at the free surface of the MCL which helps back the oddity found when using the 2 region model. When the authors studied the three region model for the MCL they found that it matched all of the oxygen tension profiles quite well. It matched the profiles so well that they didn't need to include particular parameters to model all of the data since there were 3 layers in the MCL. When previous research was done on spheroids there tended to be a zone where there was no oxygen diffusion. In the MCL there was no evidence of such a zone which may be because of the large step size of the electrode in the MCL. Everything was dependent on specific oxygen permeability and that was all calculated through thickness of the MCL. In fact there is a possibility that every diffusion difference between certain profiles can come from cell lines and other cell lines need to be tested to further resolve any issue presented. Oxygen Consumption of MCL One of the most interesting results of this paper was the idea that the metabolism of the cells on the membrane surface was larger than the rate of cells exposed to a free surface. This result could possibly be explained because the cells on the membrane surface was exposed to a collagen matrix, which was not the case for the free surface cells. On the whole the average oxygen consumption decreased with the increasing MCL thickness. This helped illustrate the correlation between anoxia and MCL thickness. The fact there is a thick layer of MCL, there is not enough oxygen just due to the sheer density around the region. This experiment cannot be fully evaluated in terms of oxygen consumption because only three MCL were exposed to the oxygen levels. The goal was not to seek the effects of oxygen level on consumption but to simply fit oxygen tension profiles under several different conditions. However, the overall consensus was that the consumption parameters decreased in the MCL in the 5% environment versus the 65% environment. Figure 4 shows the differences between oxygen consumption along the membrane surface and oxygen consumption along the free surface. Future Applications This experiment can be expanded out towards certain aspects of cancer research. Issues such as growth conditions, oxygen levels, glucose concentration and other factors can model tumor metabolism. Hopefully in the future there can be an impact towards understanding pharmaceutical drugs and their effect with respect to oxygenation. The MCL model brings a new aspect to the study of tumor oxygenation and can be very useful in the future. Our PresentationOur presentation (Thurimella, Kauffman) in .pdf format References
