Cell damage and mitigation in Swiss albino mice: experiment and modelling
Cell damage and mitigation in Swiss albino mice: experiment and modelling
Dey, Urmimala;Ghosh, Archisman;Abbas, Syed;Taraphder, A;Roy, Madhumita
2019-04-05 00:00:00
Cell damage and mitigation in Swiss albino mice: experiment and modelling Urmimala Dey Centre for Theoretical Studies, Indian Institute of Technology Kharagpur, Kharagpur 721302, India Archisman Ghosh Department of Environmental Carcinogenesis and Toxicology, Chittaranjan National Cancer Institute, 37 S. P. Mukherjee Road, Kolkata 700026, India Syed Abbas School of Basic Sciences, Indian Institute of Technology, Mandi 751005, India A Taraphder Department of Physics, Indian Institute of Technology Kharagpur, Kharagpur 721302, India Centre for Theoretical Studies, Indian Institute of Technology Kharagpur, Kharagpur 721302, India School of Basic Sciences, Indian Institute of Technology, Mandi 751005, India Madhumita Roy Department of Environmental Carcinogenesis and Toxicology, Chittaranjan National Cancer Institute, 37 S. P. Mukherjee Road, Kolkata 700026, India E-mail: urmimaladey@iitkgp.ac.in (24 April 2019) Abstract. Chronic exposure to inorganic arsenic is a potential cause of carcinogenesis. It elicits its potential by generation of ROS, leading to DNA, protein and lipid damage. Therefore, the deleterious eect of arsenic can be mitigated by quenching ROS using antioxidants. There is a homology between the protein coding regions of mice and human. Eect of these alterations in human can be mimicked in mice. Therefore to understand the underlying mechanism of arsenic toxicity and its amelioration by black tea, studies have been conducted in mice model. Long term exposure to iAs leads to tumour growth, which has been found to be alleviated by black tea. Observations reveal that black tea has two salutary eects on the growth of tumour: the rate of growth of damaged cells was appreciably reduced and an early saturation of the level of damage is achieved. To take the experimental ndings further, the experimental data have been modelled with simple dynamical equations. The arXiv:1904.10394v1 [q-bio.TO] 5 Apr 2019 Cell damage and mitigation in Swiss albino mice: experiment and modelling 2 curves obtained from in vivo studies have been tted with the data obtained from the model. The corresponding steady states and their stabilities are analyzed. 1. Introduction Arsenic, the 20th most abundant element in the earth's crust [1] is found in nature in two oxidation states i.e. +3 (arsenite) and +5 (arsenate). It is a metalloid, found in both organic and inorganic forms. Arsenic is released into the ground water from rocks containing copper or lead, either naturally or through intensive anthropogenic activities like mining etc [2]. Humans may be exposed to As primarily from air, food and drinking water [3]. Chronic exposure to inorganic arsenic (iAs) has a strong correlation with many health hazards. More than 200 million people residing in 40 countries around the world are exposed to dangerous levels of inorganic As (iAs) through ground water, the level being much higher than the permissible limit of 10 g/l [4]. In the Indian subcontinent, at several places, arsenic levels of water range from 50-3200 mg/l [3]. Exposure to iAs leads to development of several diseases including cancer. Excessive generation of reactive oxygen species (ROS) by iAs is considered to be the main reason behind the imbalance and disturbed homeostasis of individuals [5]. ROS causes damage to DNA, protein and lipid and modulates expression [6] of several genes including those implicated in cancer. Chronic arsenic exposure has also been reported to decrease the activities of antioxidant enzymes which are required to counter the eect of free radicals due to ROS generation [5]. Carcinogenesis is a multistep process involving initiation, promotion, progression and nally metastasis. In order to check the deleterious eect of iAs, it is preferable to halt the phenomenon of carcinogenesis at the initial stages itself. Chemoprevention is the use of natural, synthetic, or biologic agents in order to reverse, suppress, or prevent the progression of cancer [8]. Chemopreventive strategies to mitigate the eect of iAs-induced toxicity using natural compounds is the need of the hour. Tea, the most preferred beverage is a good antioxidant and has anticancer properties [7]. Owing to its antioxidant properties, black tea can counter the deleterious eects of ROS generated due to iAs. To elucidate the role of black tea in the prevention of iAs-induced carcinogenesis, studies have been initiated in Swiss Albino mice. Black tea can eciently quench generation of ROS, diminish DNA damage and lipid peroxidation, eliciting its anti-carcinogenic potential. 2. Materials and methods 2.1. Chemicals Preparation of Tea Extracts 1.25 gm of Assam tea and 1.25 gm of Darjeeling tea was mixed together and brewed in 100 ml of freshly boiled water for 5 minutes. The liquor was then lyophilized using Cell damage and mitigation in Swiss albino mice: experiment and modelling 3 a SCANVAC lyophilizer and stored. The lyophilized powder was then weighed and reconstituted in water and administered to mice by gavage. The catechin and thea
avin contents of the particular brand of green tea and black tea has been analysed by HPLC. Dose was determined by epidemiological evidences and previous publications from our lab [5]. Accordingly 100 l of tea extract was administered to the mice. Figure 1. Experimental datasets for the proportion of cell damage due to iAs (upper curve). The lower curve is for the group that received both iAs and black tea. Note the amelioration due to black tea, in limiting the extent of damage and causing an early saturation. 2.2. Maintenance of animals Normal, male, 4-5 week old Swiss albino mice (Mus Musculus) were selected from an inbred colony maintained in the animal house of the Institute (IAEC 1774/MR- 3/2017/9). During the course of study the animals were fed on synthetic pellets speci ed for mice-feed and had access to water ad libitum. They were maintained under standard conditions with a 12 hour light and dark period and controlled room temperature (222 C). The animals were sacri ced after euthanasia with thiopentone sodium in overdose (100 mg/kg body wt.). 2.3. Treatment The mice were divided in four groups namely I, II, III and IV. Group I: Control group, where the mice were fed with normal food and tap water. Group II: Mice were fed with black tea by gavage. Group III: Mice were given arsenic water (500 g/l) for a period of 360 days. Group IV: Black tea (100 mg/ml, thrice daily) was administered to the mice Cell damage and mitigation in Swiss albino mice: experiment and modelling 4 by gavage along with arsenic water (500 g/ml) for the same period of duration. On sacri cing the mice at dierent time intervals, blood and tissue samples were collected and stored properly to perform the experiments. Two sets of data estimating the damaged cells with time (up to about 360 days) were generated : one is only with iAs and the other, with both iAs and black tea (for this group, the data were available for about 330 days). No dataset for Group II is presented because no toxicity is found due to the application of black tea. There were already some cell damage in the mice and therefore the damage cell count does not start from zero at zero time. Fig. 1 shows both these data: the application of black tea clearly brings down the cell damage and leads to early saturation. 2.4. Estimation of ROS generation Generation of Reactive Oxygen species was estimated in the blood of the mice according to Balasubramanyam et al. [9] with slight modi cations. Leukocytes were isolated from blood using solution A (NH Cl dissolved in TRIS at pH 7.2) and Solution B (meso- 0 0 inositol), then suspended in HBS (Hepes Buered Saline pH 7.4). 10 mM 2 ,7 dichloro- dihydro-
uorescein (DCFH-DA) was added and incubated for 45 minutes in the dark. The
uorescent intensity of the DCF was recorded at 530 nm . 2.5. Single Cell Gel Electrophoresis DNA damage was estimated by the Single Cell Gel Electrophoresis (SCGE) or Comet assay [10], using leukocytes isolated from blood. 20 l of leukocytes suspended in 0.6% Low melting agarose (LMA) and were layered over frosted slide pre-coated with 0.75% Normal melting agarose (NMA). Cells were subjected to lysis by immersing them in lysis buer (2.5M NaCl, 0.1M Na EDTA, 10mM TRIS, 0.3M NaOH, 10%DMSO and 1% TritonX; pH 10) for overnight at 4 C. Slides were then presoaked in electrophoresis buer (10mM NaOH, 0.2 m Na EDTA; pH>13) for 20 minutes and subsequently subjected to electrophoresis for 25 minutes at 15V, 220mA. After staining with ethidium bromide (EtBr), slides were viewed under a
uorescent microscope. 2.6. Lipid peroxidation assay Damage to the cell membrane due to ROS generation was measured using the lipid peroxidation assay following the principle of Okhawa et al. [11], with slight modi cations. The liver tissue of the sacri ced animals were collected and homogenized. 10% SDS, 20% acetic acid and 0.8% TBA was added to the tissue homogenate. It was then placed in boiling water bath for 1 hour and then immediately transferred to ice for 10 minutes. The samples were then centrifuged at 2500 rpm for 10 minutes. The supernatant was collected and absorbance was measured at 535 nm. Lipid peroxidation was calculated as the number of moles of MDA (Malondialdehyde) generated. Cell damage and mitigation in Swiss albino mice: experiment and modelling 5 3. Model for growth of cell damage and mitigation The dynamics of the eects of iAs and black tea on the colony of Swiss albino mice are modelled by a set of simple mathematical equations. As mentioned above, we have two datasets that represent the growth of damaged cells in the groups of mice studied, one with only iAs (at 500 g/l) and the other with tea (100 mg/ml, thrice daily), acting, as it does, as a mitigating agent and iAs with the same dose as above. The extent of cell damage is proportional to the generation of reactive oxygen species (ROS). First we ignore the eect of immunity and try to understand the eects of iAs and tea only. The eects of immune cell response will be brought on in the next section. Our aim in this and the following section is to try to model actual data that are obtained from our experiment using standard mathematical models. There are excellent review articles on such models and their merits and demerits (see, for example, Altrock, et al. [12] and Pillis, et al., [13]). We do not, therefore, digress on these details here. The dynamics of the damaged cells may be represented by the usual growth equation: dn n = rn(1 ) (1) dt K N in the above equation is a measure of the damaged cells as interpreted from ROS data. We de ne n = where N is the total number of cells at the beginning. Note that one does not need to know the value of either N or N , only the ratio (percentage) of damaged cells, represented by n, is needed for the model. There are two unknown constants above, r and K = B=N , where B is the well-known `carrying capacity'. Now iAs is administered at a constant rate daily from outside and is incorporated via the term An in the above equation while the mitigating eects of tea is introduced via the term Tn [14]. Then the overall equation for the cell growth is dn n = rn(1 ) + An Tn (2) dt K The constants and are unknown, to be tted from the data. The equations governing the concentration of iAs and tea, A and T respectively, are similar [14]. In principle, there could be additional eects, but in the minimal model we consider the eects of iAs and tea require the two terms shown below: dA = A
A (3) 0 A dt dT = T
T (4) 0 T dt Here A and T are the constant rate of administration of iAs and tea externally. Here 0 0 again, other eects of iAs and tea on the body could be added, which we will consider in the future. One eect that is usually considered important, but we neglect here (discussed in the next section), is the immune-surveillance. Usually, A and T , the 0 0 so-called dosage, are known in terms of mg (or ml)/day. However, to incorporate them Cell damage and mitigation in Swiss albino mice: experiment and modelling 6 in the Eq. (2) above, one needs to know the concentration thereof at the cellular level, which, at this moment, is not available to us. Therefore we leave them to tting. For the dataset with iAs only, the equation for tea is obviously redundant. The expulsion of iAs and tea from the body is also accounted for via decay rates (rough estimates of half-lives of elimination for these are available),
and
, from the body. A T Note that dierent ingredients decay dierently, depending also on the organ considered. So some approximate values are chosen from literature [15], though a considerable degree of latitude could not be ruled out. Solving Eq. (3) and Eq. (4) t
t A A A(t) = (1 e ) + A e (5) t
t T T T (t) = (1 e ) + T e (6) where, A = A(0) and T = T (0). Substituting the expressions of A(t) and T (t) in i i Eq. (2) : dn n t
t
t
t A A T T = rn(1 ) + mn(1 e ) + pne qn(1 e ) sne (7) dt K A T 0 0 where, m = , p = A , q = and s = T . It is safe to assume that there is no i i A T Figure 2. Curves tted to the given datasets. Here, red and green dots represent experimental data and the blue and orange solid lines represent the tted curves. prior presence of iAs and tea in the mice. Normal tap water which is given to the mice is tested for the presence of iAs (by atomic absorption spectroscopy (AAS)). The tap water was found to be free of the metalloid. Although the presence of iAs at a very small level in the body cannot be ruled out, but that is still negligible for the model. Cell damage and mitigation in Swiss albino mice: experiment and modelling 7 The initial values of iAs and tea at time t = 0 i.e., A , T = 0. This implies p; s = 0. In i i 1 1 our case,
= 0.04 day and
= 4 day . We reiterate that these values are merely A T approximations. Such rates depend on a large number of factors, primarily the organ being considered, the method of administration and the measurement protocol. The values used here are similar to Hughes et al. [15] Solving and tting Eq. (7) numerically to the given datasets, we obtain the values of the unknown coecients as : r = 0.00524, K = 0.14706, m = 0.01281 for the rst dataset and r = 0.01380, K = 0.21644, m = 0.00792, q = 0.00041 for the second dataset. Note that for the rst dataset, the eect of black tea is not included. So, q is not relevant for the rst dataset. Fig. 2 shows the curves tted to the given datasets, red and green dots represent experimental data. The tting appears to be excellent for both the curves. Note that although iAs has a much longer half-life compared to black tea, the competition between the damage by iAs and mitigation by tea is clear from the data and the model. The saturation value and time for the curve with tea and iAs is considerably lower than the upper curve with iAs alone. This is a redeeming feature from the experiment and the model that the eect of tea is quite appreciable. The fact that curves for iAs alone and iAs and tea together could both be tted with similar models gives con dence that the model could be used for further predictions. One such prediction is that the damage appears to saturate in both cases, earlier with tea. Such saturation would imply that the damage could be limited by the use of an external mitigant like tea. We have neglected the eects of immunogenicity above, which will abet early saturation. Although in the present paper we took the actual data and modelled them, which makes modelling much more relevant, it is important to underline the remit of data- tting by mathematical equations : such exercises should be carried out with due caution. In the absence of cellular level concentration of iAs and black tea, and the inherent uncertainty and
uctuation in their values at the site, the results of the model should be taken as indicative only. In addition, as we show below, the eect of immune-cell response dramatically alters the model and make it closer to reality [16]. 4. Growth model including the eect of immune cells Let I (t) denote the number of immune cells at any instance of time t with a constant in
ux rate of I and if the death rate of the immune cells be
, the overall dynamics of 0 I the system may be represented by the following equations : dn n = rn(1 ) + An nT nI (8) dt K dA = A
A (9) 0 A dt dT = T
T (10) 0 T dt dI In = I + c In
I (11) 0 1 I dt + n 1 Cell damage and mitigation in Swiss albino mice: experiment and modelling 8 The decay rates of iAs and black tea (
and
) are given in the previous section. We A T assume the constant in
ux rate of the immune cells I = 0.3 and the death rate of the immune cells
= 0.0208 day [17, 18]. Eqs. (9) and (10) can be solved analytically as given in Eqs. (5) and 6. Substituting the expressions of A(t) and T (t) in Eq. (8) as before, we obtain dn n t
t
t
t A A T T = rn(1 ) + mn(1 e ) + pne qn(1 e ) sne nI (12) dt K A T 0 0 where, the unknown coecients m (= ), p (= A ), q (= ), s (= T ), , , i i 1 A T and c are to be tted from the experimental data. We consider that initially at t = 0, no iAs, tea and immune cells are present. This implies A , T = 0. Therefore, the i i coecients p, s = 0. Numerically solving Eqs. (11) and (12) and after tting to the experimental datasets, we obtained the values of the unknown coecients as : r = 0.00173, K = 0.03166, m = 0.02802, = 0.00055, = 0.95103, = 1.30187, c = 0.63433 for the rst dataset 1 1 and r = 0.60516, K = 8.97523, m = 0.08006, q = 0.61760, = 0.00379, = 0.02733, = 2.15877, c = 0.02718 for the second dataset. For the rst dataset, the coecients 1 1 are calculated without considering the eect of black tea. The tted curves are given in Fig. 3 along with the experimental datasets, which are shown by the orange and blue dots. The existence and uniqueness of the solutions are discussed in the Appendix. Figure 3. Curves tted to the given datasets. Here, orange and blue dots represent experimental data and the green and red solid lines represent the tted curves. Cell damage and mitigation in Swiss albino mice: experiment and modelling 9 4.1. Equilibrium points 4.1.1. Tea-free equilibria In order to study the equilibrium solutions, we rst consider the case without tea. We look for the steady state solutions: n _ = 0 A = 0 I = 0 For experimental relevance, we look for the non-zero solution for n only. The three Eqs. above generate three null surfaces, the intersection of which give the equilibrium point (n , A , I ). The non-trivial solutions are : A I n = K 1 + (13) r r A = (14) I ( + n ) 0 1 I = (15) c n + c n +
+
n n 1 1 1 I 1 I Although the trivial equilibrium (no-tumour) n = 0 is studied for the sake of completeness by other authors [16], it does not correspond to any relevant experimental outcome. Therefore we neglect the trivial point. In Fig. 4, we plot the three null- surfaces given by Eqs. (13 15) using the values of the coecients, as obtained from the t. Solving Eqs. (13) and (15) numerically and using the obtained values of the coecients, we nd that equilibrium occurs at the point (A , I , n ) = (10, 4.010730, 0.502892), where we consider the constant in
ux rate of iAs is xed at A = 0.4. Once the non-trivial points of equilibrium are obtained, we study their linear stability. In order to ensure the positivity of equilibria, we need r + A > I and > c +
: 1 1 I The other condition we obtain by observing the roots of c n +(c +
)n +
= 1 1 1 I I 1 2 1 c n + bn +
= (n a )(n a ); where a = ( b b 4c
): 1 I 1 1 2 1;2 1 I 1 2c 4.1.2. Stability analysis We construct the Jacobian matrix: 0 1 2r A I + r n n n B C 0
0 J = @ A A I I n n c I 0 c n 2 1 1 I +n +n 1 ( +n ) 1 0 1 L n n B C = 0
0 @ A M 0 P I I n n 2r where, L = A I +r n , M = c I and P = c n 2 1 1 I K +n +n 1 ( +n ) 1 The three eigenvalues of J are given by: =
< 0 (16) 1 A Cell damage and mitigation in Swiss albino mice: experiment and modelling 10 = (L + P ) (L P ) 4n Q (17) For stability, eigenvalues must be negative. So, we have to nd the region of negative eigenvalues along the n I A axes. Substituting the values of n , A and I in Eq. (17), Figure 4. The three null-surfaces given by Eqs. (13 15) using the values of the coecients obtained from the t. Here, we have taken A = 0.4. The equilibrium points (n , A , T , I ) appear at the intersections of these surfaces. we get < 0. Since is already negative, it corresponds to stable equilibrium. 4.1.3. Equilibrium with Tea Following the same procedure, we solve the equations: n _ = 0 A = 0 T = 0 I = 0 The non-zero solutions are A T I n = K 1 + (18) r r r Cell damage and mitigation in Swiss albino mice: experiment and modelling 11 A = (19) T = (20) I ( + n ) 0 1 I = (21) c n + c n +
+
n n 1 1 1 I 1 I To ensure the positivity of equilibria, we need r + A > T + I and > c +
: 1 1 I The equilibrium points (T , I , n ) and (A , I , n ) are obtained from the upper and lower gures in Fig. 5. They occur at the intersections of these surfaces as before. 4.1.4. Stability analysis We write the 44 Jacobian Matrix now 0 1 L n n n B C 0
0 0 B A C J = B C @ 0 0
0 A M 0 0 P I I n 2r where, L = A T I + r n , M = c I and P = 2 1 K +n 1 ( +n ) c n
. The two eigenvalues of the Jacobian, =
and =
are 1 I 1 A 2 T +n both negative. The other two eigenvalues are = (L + P ) (L P ) 4n M (22) Again we have to look for the regions where < 0. In Fig. 5, we show the three null- surfaces given by Eqs. (18 21) using the values of the coecients, as obtained from the t; the upper one is plotted keeping A xed at A = 10. In the lower one, the null-surfaces are plotted keeping T xed at T = 0:1. The equilibrium points (T , I , n ) and (A , I , n ) are the intersections of these surfaces. Our calculations show that equilibrium occurs at the point (A , T , I , n ) = (10, 0.1, 11.48990, 0.35624), where the constant in
ux rates of iAs and black tea are kept xed at A , T = 0.4. When we 0 0 substitute these values of n , A , T and I in Eq. (22), we nd that < 0. Since the other two eigenvalues are already negative, this implies that the values of the coecients obtained from the t actually lead to a stable equilibrium. 4.2. Model with delay The immune response in any living organism is usually a delayed process and therefore it is natural to assume a delay in I (t) in the functional response term. In this situation, Cell damage and mitigation in Swiss albino mice: experiment and modelling 12 Figure 5. The three null-surfaces given by Eqs. (18 21) using the values of the coecients, as obtained from the t. Here, we have taken A = 0.4 and T = 0.4. 0 0 The upper one is plotted keeping A xed at A = A = 10 and in the lower one, the null-surfaces are plotted keeping T xed at T = T = 0:1. The equilibrium points (T , I , n ) and (A , I , n ) occur at the intersections of these surfaces. the overall dynamics of the system may be represented by the following equations: dn(t) n(t) = rn(t)(1 ) + A(t)n(t) n(t)T (t) n(t)I (t) dt K dA(t) = A
A(t) 0 A dt dT (t) = T
T (t) 0 T dt dI (t) I (t )n(t) = I + c I (t)n(t)
I (t); (23) 0 1 I dt + n(t) 1 Cell damage and mitigation in Swiss albino mice: experiment and modelling 13 with history condition n(0) = n ; A(0) = A ; T (0) = T ; I (t) = (t) for t 2 [ ; 0]: 0 0 0 4.2.1. Stability of Equilibrium without Tea Now we discuss the stability of non-trivial equilibrium solutions. It is to note that the set of equilibria are the same for this system as with Eqs. (8) (11). First, we analyze tea free system. We construct the Jacobian matrix : 0 1 2r (A I + r n ) n n B C 0
0 (; ) = det A @ A I I n n e c I 0 ( c n
) 2 1 1 I +n +n 1 ( +n ) 1 0 1 L n n B C = det 0
0 @ A M 0 P + P e 1 2 I I n 2r where, L = A I +r n , M = c I ; P = c n
and P 2 1 1 I 1 2 K +n ( +n ) 1 1 = : We use the following stability criteria for a nonlinear delay equation from [19]. +n Theorem 1 The necessary and sucient conditions for an equilibrium to be asymptotically stable for all 0 are the following: (1) The real part of the roots of (; 0) = 0 is negative. (2) For all real and all 0; (i; ) 6= 0. One eigenvalue is of course
and others are given by the root of the following equation (L )(P + P e ) + n M = 0: 1 2 Simplifying, we get + ( L P P e ) + L (P + P e ) + n M = 0: 1 2 1 2 We have already obtained condition for non delay case. To examine the other condition, let us substitute = i in the above equation, 2 i i + ( L P P e )i + L (P + P e ) + n M = 0: 1 2 1 2 Separating real and imaginary part, we obtain P sin( ) + P cos( ) + n M = 0 2 2 L P cos( ) P sin( ) = 0: (24) 2 2 Squaring and adding, we get 2 2 2 2 2 2 (n M ) + L = P + P ; 2 2 which further implies 2 2 2 2 2 4 2 2 ((P ) + 2n M ) + n M + L P = 0: 2 2 Cell damage and mitigation in Swiss albino mice: experiment and modelling 14 The roots are given by " # 2 2 2 2 2 2 2 2 2 = ((P ) + 2n M )) ((P ) + 2n M ) 4( n M + L P ) 2 2 2 So condition (2) of theorem violates implies that there exists a such that the above 2 2 2 2 relation holds. It holds if n M +L > P ; so the equilibrium is stable if in addition 2 2 2 2 to the condition for non-delay case (subsection 4.1.2), we have n M + L < P : Remark 1 Stability of Equilibrium with Tea : In this case, since the equation in I 2 2 2 2 does not depend on T; we get the same condition for stability n M + L < P in addition to the condition given in subsection (4.1.4). 5. Conclusion An attempt has been made to model the data obtained from the in vivo studies, using Swiss albino mice. Mice are maintained under the in
uence of iAs for a long period, with iAs-loaded water as the only source of drinking water. Black tea infusion is administered on a regular basis to a subset of the mice. The data obtained on the cell-damage in these two groups of mice were then modelled by the standard dynamical equations. It is redeeming to note that even simple models could t the data well and predict long-term behaviour of the cell damage by iAs and its alleviation by black tea. Further studies using experimentally obtained tting parameters would reveal more interesting underlying dynamics that are testable in real systems. The predictions thereof would be reliable and useful for the study of growth and limitation of tumours. Appendix We discuss the existence and uniqueness of the solution of the dynamical equations considered in the text above. Appendix A. Existence and uniqueness of the solution We consider the immune-free case (Eq. 7) rst. The case with immunity will follow similarly. The problem can be written as dy = f (t; y); y(0) = y ; t 2 [t ; t + T ]; 0 0 0 dt t
t
t
t A A T T where f (t; y) = ry(1 ) +my(1 e ) +pye qy(1 e ) sye . Since we are interested in only positive solutions, hence we assume the positivity of the solution when y is positive. By using bounds of exponentials, we obtain dy y ry(1 ) + my + py dt K (r + m + p)y y y((r + m + p) y): (A.1) K Cell damage and mitigation in Swiss albino mice: experiment and modelling 15 K(r+m+p) Hence, we obtain lim sup y(t) := M: Furthermore, we can write equation t!1 as [ + ]dy (r + m + p)dt; y r + m + p y which gives r + m + p y(t) : (r+m+p)t + ce Similarly, we get dy y ry(1 ) qy sy dt K (r q s)y y y((r q s) y): (A.2) K(r q s) Hence, we obtain lim inf y(t) := m: Since all coecients are positive, t!1 hence if y is positive, the solution remains positive in this case. So, we have achieved the bounds and positivity of the solution. In order to ensure the existence and uniqueness of solution, we need to check the Lipschitz condition. Suppose y ; y 2 R; we have 1 2 jf (t; y ) f (t; y )j 1 2 y y 1 2 rjy (1 ) y (1 )j + m(1 e )jy y j 1 2 1 2 K K t
t
t A T A + pe jy y j + q(1 e )jy y j + se jy y j 1 2 1 2 1 2 rjy y j + jy y j2M + (m + p + q + s)jy y j 1 2 1 2 1 2 2Mr (r + m + p + q + s + )jy y j := Ljy y j: (A.3) 1 2 1 2 Hence f is Lipschitz with respect to y: It is easy to see that f is continuous with respect to t: So, we can apply Picard theorem to ensure the existence and uniqueness of solution. Next, we discuss the existence and uniqueness of the solution for Eqs. (8) (11). One can convert this system of Eqs. in vector from by de ning X = (n; A; T; I ) : dX = F (X ) = AX + F (X ); where dt 2 3 r 0 0 0 6 7 0
0 0 6 A 7 A = 6 7 40 0
0 5 0 0 0 and 2 3 rn + An nT nI 6 7 6 0 7 F (X ) = : 6 7 4 T 5 In I + c In 0 1 +n It is straightforward to check that the function F is Lipschitz by assuming the boundedness, i.e. there exists a constant L such that kF (X ) F (X )k L kX X k 1 2 F 1 2 Cell damage and mitigation in Swiss albino mice: experiment and modelling 16 Also, continuity with respect to t is evident since the F does not depend on t: Hence, the existence and uniqueness of the solution is guaranteed by Picard-Lindelo theorem. Acknowledgments UD would like to acknowledge the Ministry of Human Resource Development (MHRD), India for research fellowship. References [1] I. Herath, M. Vithanage, J. Bundschuh, J. P. Maity, and P. Bhattacharya, Natural Arsenic in Global Groundwaters: Distribution and Geochemical Triggers for Mobilization, Curr Pollution Rep 2 : 68 89 (2016). [2] V. D. Martinez, E. A. Vucic, D. D. B.-Santos, L. Gil, and W. L. Lam, Arsenic Exposure and the Induction of Human Cancers, Journal of Toxicology 2011, 431287 (2011). [3] J. Biswas, S. Roy, S. Mukherjee, D. Sinha, and M. Roy, Indian spice curcumin may be an eective strategy to combat the genotoxicity of arsenic in Swiss albino mice, Asian Paci c J Cancer Prev 11(1) : 239 246 (2010). [4] I. Khairul, Q. Q. Wang, Y. H. Jiang, C. Wang, and H. Naranmandura, Metabolism, toxicity and anticancer activities of arsenic compounds, Oncotarget, 8(14) : 23905 23926 (2017). [5] D. Sinha, S. Roy, and M. Roy, Antioxidant potential of tea reduces arsenite induced oxidative stress in Swiss albino mice, Food and Chemical Toxicology 48(4) : 1032 1039 (2010). [6] R. R.-Ramos, L. L.-Carrillo, A. D. R.-Perez, A. D. V.-Ru z, M. E. Cebrian, Sodium arsenite induces ROS generation, DNA oxidative damage, HO-1 and c-Myc proteins, NF-B activation and cell proliferation in human breast cancer MCF-7 cells, Mutat Res. 674(1-2) : 109 115 (2009). [7] D. Sinha, M. Roy, M. Siddiqi, and R. K. Bhattacharya, Modulation of arsenic induced DNA damage by tea as assessed by single cell gel electrophoresis, International Journal of Cancer Prevention 2(2), 143 154 (2005). [8] A. S. Tsao, E. S. Kim and W. K. Hong, Chemoprevention of Cancer, CA A Cancer J Clin 54 : 150 180 (2004). [9] M. Balasubramanyam, A. A. Koteswari, R. S. Kumar, S. F. Monickaraj, J. U. Maheswari, and V. Mohan, Curcumin-induced inhibition of cellular reactive oxygen species generation: Novel therapeutic implications, J. Biosci. 28(6) : 715 721 (2003). [10] N. P. Singh, M. T. McCoy, R. R. Tice, E. L. Schneider, A Simple Technique for Quanti cation of Low Levels of DNA Damage in Individual Cells, Exp Cell Res. 175 : 184 191 (1988). [11] H. Ohkawa, N. Ohishi, and K. Yagi, Assay for Lipid Peroxides in Animal Tissues by Thiobarbituric Acid Reaction, Analytical Biochemistry 95 : 351 358 (1979). [12] P. M. Altrock, L. L. Liu, and F. Michor, The mathematics of cancer: integrating quantitative models, Nature Reviews Cancer15, 730 (2015). [13] L.G. de Pillis and A.E. Radunskaya, A Mathematical Tumor Model with Immune Resistance and Drug Therapy : an Optimal Control Approach, Journal of Theoretical Medicine 3 : 79 100 (2001). [14] P. Liu and X. Liu, Dynamics of a tumor-immune model considering targeted chemotherapy, Chaos, Solitons and Fractals 98, 7 13 (2017). [15] M. F. Hughes, E. M. Kenyon, B. C. Edwards, C. T. Mitchell, L. M. D. Razo, and D. J. Thomas, Accumulation and metabolism of arsenic in mice after repeated oral administration of arsenate, Toxicology and Applied Pharmacology 191(3) : 202210 (2003). [16] L. G. de Pillis and A. Radunskaya, The Dynamics of an Optimally Controlled Tumor Model: A Case Study, Mathematical and Computer Modelling 37, 1221 1244 (2003). Cell damage and mitigation in Swiss albino mice: experiment and modelling 17 [17] L. G. de Pillis, W. Gu, K. R. Fister, T. Head, K. Maples, A. Murugan, T. Neal, and K. Yoshida, Chemotherapy for tumors: An analysis of the dynamics and a study of quadratic and linear optimal controls, Mathematical Biosciences 209, 292315 (2007). [18] V. A. Kuznetsov, I. A. Makalkin, M. A. Taylor, and A. S. Perelson, Nonlinear dynamics of immunogenic tumors: parameter estimation and global bifurcation analysis, Bull. Math. Bio. 56(2): 295 321 (1994). [19] Theorem 3.7.1 in K. Gopalsamy, Stability and Oscillation in Delay Dierential Equations of Population Dynamics, Kluwer Academic, Dordrecht, (1992).
http://www.deepdyve.com/assets/images/DeepDyve-Logo-lg.pngQuantitative BiologyarXiv (Cornell University)http://www.deepdyve.com/lp/arxiv-cornell-university/cell-damage-and-mitigation-in-swiss-albino-mice-experiment-and-nujlPiKX3h
Cell damage and mitigation in Swiss albino mice: experiment and modelling