1
|
Frieboes HB, Smith BR, Wang Z, Kotsuma M, Ito K, Day A, Cahill B, Flinders C, Mumenthaler SM, Mallick P, Simbawa E, AL-Fhaid AS, Mahmoud SR, Gambhir SS, Cristini V. Predictive Modeling of Drug Response in Non-Hodgkin's Lymphoma. PLoS One 2015; 10:e0129433. [PMID: 26061425 PMCID: PMC4464754 DOI: 10.1371/journal.pone.0129433] [Citation(s) in RCA: 24] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/29/2015] [Accepted: 05/09/2015] [Indexed: 12/20/2022] Open
Abstract
We combine mathematical modeling with experiments in living mice to quantify the relative roles of intrinsic cellular vs. tissue-scale physiological contributors to chemotherapy drug resistance, which are difficult to understand solely through experimentation. Experiments in cell culture and in mice with drug-sensitive (Eµ-myc/Arf-/-) and drug-resistant (Eµ-myc/p53-/-) lymphoma cell lines were conducted to calibrate and validate a mechanistic mathematical model. Inputs to inform the model include tumor drug transport characteristics, such as blood volume fraction, average geometric mean blood vessel radius, drug diffusion penetration distance, and drug response in cell culture. Model results show that the drug response in mice, represented by the fraction of dead tumor volume, can be reliably predicted from these inputs. Hence, a proof-of-principle for predictive quantification of lymphoma drug therapy was established based on both cellular and tissue-scale physiological contributions. We further demonstrate that, if the in vitro cytotoxic response of a specific cancer cell line under chemotherapy is known, the model is then able to predict the treatment efficacy in vivo. Lastly, tissue blood volume fraction was determined to be the most sensitive model parameter and a primary contributor to drug resistance.
Collapse
Affiliation(s)
- Hermann B. Frieboes
- Department of Bioengineering, University of Louisville, Louisville, KY, 40202, United States of America
- James Graham Brown Cancer Center, University of Louisville, Louisville, KY, 40202, United States of America
- Department of Pathology, University of New Mexico, Albuquerque, NM, 87131, United States of America
- * E-mail: (HBF); (SGG); (VC)
| | - Bryan R. Smith
- Molecular Imaging Program at Stanford (MIPS), Department of Radiology, Stanford University, Stanford, CA, 94305, United States of America
| | - Zhihui Wang
- Department of Pathology, University of New Mexico, Albuquerque, NM, 87131, United States of America
| | - Masakatsu Kotsuma
- Molecular Imaging Program at Stanford (MIPS), Department of Radiology, Stanford University, Stanford, CA, 94305, United States of America
| | - Ken Ito
- Molecular Imaging Program at Stanford (MIPS), Department of Radiology, Stanford University, Stanford, CA, 94305, United States of America
| | - Armin Day
- Department of Pathology, University of New Mexico, Albuquerque, NM, 87131, United States of America
| | - Benjamin Cahill
- Department of Bioengineering, University of Louisville, Louisville, KY, 40202, United States of America
| | - Colin Flinders
- Department of Biological Chemistry, University of California at Los Angeles, Los Angeles, CA, 90095, United States of America
- Center for Applied Molecular Medicine, University of Southern California, Los Angeles, CA, 90033, United States of America
| | - Shannon M. Mumenthaler
- Center for Applied Molecular Medicine, University of Southern California, Los Angeles, CA, 90033, United States of America
| | - Parag Mallick
- Molecular Imaging Program at Stanford (MIPS), Department of Radiology, Stanford University, Stanford, CA, 94305, United States of America
| | - Eman Simbawa
- Department of Mathematics, Faculty of Science, King Abdulaziz University, Jeddah, 21589, Saudi Arabia
| | - A. S. AL-Fhaid
- Department of Mathematics, Faculty of Science, King Abdulaziz University, Jeddah, 21589, Saudi Arabia
| | - S. R. Mahmoud
- Department of Mathematics, Faculty of Science, King Abdulaziz University, Jeddah, 21589, Saudi Arabia
| | - Sanjiv S. Gambhir
- Molecular Imaging Program at Stanford (MIPS), Department of Radiology, Stanford University, Stanford, CA, 94305, United States of America
- Department of Bioengineering, Stanford University, Stanford, CA, 94305, United States of America
- Department of Materials Science & Engineering, and Bio-X, Stanford University, Stanford, CA, 94305, United States of America
- * E-mail: (HBF); (SGG); (VC)
| | - Vittorio Cristini
- Department of Pathology, University of New Mexico, Albuquerque, NM, 87131, United States of America
- Department of Mathematics, Faculty of Science, King Abdulaziz University, Jeddah, 21589, Saudi Arabia
- Department of Chemical Engineering and Center for Biomedical Engineering, University of New Mexico, Albuquerque, NM, 87131, United States of America
- * E-mail: (HBF); (SGG); (VC)
| |
Collapse
|
2
|
An integrated computational/experimental model of lymphoma growth. PLoS Comput Biol 2013; 9:e1003008. [PMID: 23555235 PMCID: PMC3610621 DOI: 10.1371/journal.pcbi.1003008] [Citation(s) in RCA: 29] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/05/2012] [Accepted: 02/13/2013] [Indexed: 12/27/2022] Open
Abstract
Non-Hodgkin's lymphoma is a disseminated, highly malignant cancer, with resistance to drug treatment based on molecular- and tissue-scale characteristics that are intricately linked. A critical element of molecular resistance has been traced to the loss of functionality in proteins such as the tumor suppressor p53. We investigate the tissue-scale physiologic effects of this loss by integrating in vivo and immunohistological data with computational modeling to study the spatiotemporal physical dynamics of lymphoma growth. We compare between drug-sensitive Eμ-myc Arf-/- and drug-resistant Eμ-myc p53-/- lymphoma cell tumors grown in live mice. Initial values for the model parameters are obtained in part by extracting values from the cellular-scale from whole-tumor histological staining of the tumor-infiltrated inguinal lymph node in vivo. We compare model-predicted tumor growth with that observed from intravital microscopy and macroscopic imaging in vivo, finding that the model is able to accurately predict lymphoma growth. A critical physical mechanism underlying drug-resistant phenotypes may be that the Eμ-myc p53-/- cells seem to pack more closely within the tumor than the Eμ-myc Arf-/- cells, thus possibly exacerbating diffusion gradients of oxygen, leading to cell quiescence and hence resistance to cell-cycle specific drugs. Tighter cell packing could also maintain steeper gradients of drug and lead to insufficient toxicity. The transport phenomena within the lymphoma may thus contribute in nontrivial, complex ways to the difference in drug sensitivity between Eμ-myc Arf-/- and Eμ-myc p53-/- tumors, beyond what might be solely expected from loss of functionality at the molecular scale. We conclude that computational modeling tightly integrated with experimental data gives insight into the dynamics of Non-Hodgkin's lymphoma and provides a platform to generate confirmable predictions of tumor growth. Non-Hodgkin's lymphoma is a cancer that develops from white blood cells called lymphocytes in the immune system, whose role is to fight disease throughout the body. This cancer can spread throughout the whole body and be very lethal – in the US, one third of patients will die from this disease within five years of diagnosis. Chemotherapy is a usual treatment for lymphoma, but the cancer can become highly resistant to it. One reason is that a critical gene called p53 can become mutated and help the cancer to survive. In this work we investigate how cells with this mutation affect the cancer growth by performing experiments in mice and using a computer model. By inputting the model parameters based on data from the experiments, we are able to accurately predict the growth of the tumor as compared to tumor measurements in living mice. We conclude that computational modeling integrated with experimental data gives insight into the dynamics of Non-Hodgkin's lymphoma, and provides a platform to generate confirmable predictions of tumor growth.
Collapse
|
3
|
Exploration of factors affecting the onset and maturation course of follicular lymphoma through simulations of the germinal center. Bull Math Biol 2009; 71:1432-62. [PMID: 19412639 DOI: 10.1007/s11538-009-9408-8] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/22/2007] [Accepted: 02/10/2009] [Indexed: 10/20/2022]
Abstract
Genetic mutations frequently observed in human follicular lymphoma (FL) B-cells result in aberrant expression of the anti-apoptotic protein bcl-2 and surface immunoglobulins (Igs) which display one or more novel variable (V) region N-glycosylation motifs. In the present study, we develop a simulation model of the germinal center (GC) to explore how these mutations might influence the emergence and clonal expansion of key mutants which provoke FL development. The simulations employ a stochastic method for calculating the cellular dynamics, which incorporates actual IgV region sequences and a simplified hypermutation scheme. We first bring our simulations into agreement with experimental data for well-characterized normal and bcl-2(+) anti-hapten GC responses in mice to provide a model for understanding how bcl-2 expression leads to permissive selection and memory cell differentiation of weakly competitive B-cells. However, as bcl-2 expression in the GC alone is thought to be insufficient for FL development, we next monitor simulated IgV region mutations to determine the emergence times of key mutants displaying aberrant N-glycosylation motifs recurrently observed in human FL IgV regions. Simulations of 26 germline V(H) gene segments indicate that particular IgV regions have a dynamical selective advantage by virtue of the speed with which one or more of their key sites can generate N-glycosylation motifs upon hypermutation. Separate calculations attribute the high occurrence frequency of such IgV regions in FL to an ability to produce key mutants at a fast enough rate to overcome stochastic processes in the GC that hinder clonal expansion. Altogether, these simulations characterize three pathways for FL maturation through positively selected N-glycosylations, namely, via one of two key sites within germline V(H) region gene segments, or via a site in the third heavy chain complementarity-determining region (CDR-H3) that is generated from VDJ recombination.
Collapse
|