Treatment scheduling effects on the evolution of drug resistance in heterogeneous cancer cell populations.
Ontology highlight
ABSTRACT:
SUBMITTER: Patwardhan GA
PROVIDER: S-EPMC8154902 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:Recent drug discovery and development efforts have created a large arsenal of targeted and chemotherapeutic drugs for precision medicine. However, drug resistance remains a major challenge as minor pre-existing resistant subpopulations are often found to be enriched at relapse. Current drug design has been heavily focused on initial efficacy, and we do not fully understand the effects of drug selective pressure on long-term drug resistance potential. Using a minimal two-population model, taking into account subpopulation proportions and growth/kill rates, we modeled long-term drug treatment and performed parameter sweeps to analyze the effects of each parameter on therapeutic efficacy. We found that drugs with the same overall initial kill may exert differential selective pressures, affecting long-term therapeutic outcome. We validated our conclusions experimentally using a preclinical model of Burkitt's lymphoma. Furthermore, we highlighted an intrinsic tradeoff between drug-imposed overall selective pressure and rate of adaptation. A principled approach in understanding the effects of distinct drug selective pressures on short-term and long-term tumor response enables better design of therapeutics that ultimately minimize relapse.
| S-EPMC5098152 | biostudies-literature
Project description:The evolution of multidrug resistance (MDR) is a pressing public health concern. Yet many aspects, such as the role played by population structure, remain poorly understood. Here, we argue that studying MDR evolution by focusing upon the dynamical equations for linkage disequilibrium (LD) can greatly simplify the calculations, generate more insight, and provide a unified framework for understanding the role of population structure. We demonstrate how a general epidemiological model of MDR evolution can be recast in terms of the LD equations. These equations reveal how the different forces generating and propagating LD operate in a dynamical setting at both the population and metapopulation levels. We then apply these insights to show how the LD perspective: (i) explains equilibrium patterns of MDR, (ii) provides a simple interpretative framework for transient evolutionary dynamics, and (iii) can be used to assess the consequences of different drug prescription strategies for MDR evolution.
| S-EPMC8208818 | biostudies-literature
Project description:Adaptation to inhibitory concentrations of the antifungal agent fluconazole was monitored in replicated experimental populations founded from a single, drug-sensitive cell of the yeast Candida albicans and reared over 330 generations. The concentration of fluconazole was maintained at twice the MIC in six populations; no fluconazole was added to another six populations. All six replicate populations grown with fluconazole adapted to the presence of drug as indicated by an increase in MIC; none of the six populations grown without fluconazole showed any change in MIC. In all populations evolved with drug, increased fluconazole resistance was accompanied by increased resistance to ketoconazole and itraconazole; these populations contained ergosterol in their cell membranes and were amphotericin sensitive. The increase in fluconazole MIC in the six populations evolved with drug followed different trajectories, and these populations achieved different levels of resistance, with distinct overexpression patterns of four genes involved in azole resistance: the ATP-binding cassette transporter genes, CDR1 and CDR2; the gene encoding the target enzyme of the azoles in the ergosterol biosynthetic pathway, ERG11; and the major facilitator gene, MDR1. Selective sweeps in these populations were accompanied by additional genomic changes with no known relationship to drug resistance: loss of heterozygosity in two of the five marker genes assayed and alterations in DNA fingerprints and electrophoretic karyotypes. These results show that chance, in the form of mutations that confer an adaptive advantage, is a determinant in the evolution of azole drug resistance in experimental populations of C. albicans.
| S-EPMC94447 | biostudies-literature
Project description:Single-cell measurement techniques can now probe gene expression in heterogeneous cell populations from the human body across a range of environmental and physiological conditions. However, new mathematical and computational methods are required to represent and analyze gene-expression changes that occur in complex mixtures of single cells as they respond to signals, drugs, or disease states. Here, we introduce a mathematical modeling platform, PopAlign, that automatically identifies subpopulations of cells within a heterogeneous mixture and tracks gene-expression and cell-abundance changes across subpopulations by constructing and comparing probabilistic models. Probabilistic models provide a low-error, compressed representation of single-cell data that enables efficient large-scale computations. We apply PopAlign to analyze the impact of 40 different immunomodulatory compounds on a heterogeneous population of donor-derived human immune cells as well as patient-specific disease signatures in multiple myeloma. PopAlign scales to comparisons involving tens to hundreds of samples, enabling large-scale studies of natural and engineered cell populations as they respond to drugs, signals, or physiological change.
| S-EPMC7682438 | biostudies-literature
Project description:Increasing body of experimental evidence suggests that anticancer and antimicrobial therapies may themselves promote the acquisition of drug resistance by increasing mutability. The successful control of evolving populations requires that such biological costs of control are identified, quantified and included to the evolutionarily informed treatment protocol. Here we identify, characterise and exploit a trade-off between decreasing the target population size and generating a surplus of treatment-induced rescue mutations. We show that the probability of cure is maximized at an intermediate dosage, below the drug concentration yielding maximal population decay, suggesting that treatment outcomes may in some cases be substantially improved by less aggressive treatment strategies. We also provide a general analytical relationship that implicitly links growth rate, pharmacodynamics and dose-dependent mutation rate to an optimal control law. Our results highlight the important, but often neglected, role of fundamental eco-evolutionary costs of control. These costs can often lead to situations, where decreasing the cumulative drug dosage may be preferable even when the objective of the treatment is elimination, and not containment. Taken together, our results thus add to the ongoing criticism of the standard practice of administering aggressive, high-dose therapies and motivate further experimental and clinical investigation of the mutagenicity and other hidden collateral costs of therapies.
| S-EPMC8491903 | biostudies-literature
Project description:Certain tumor phenomena, like metabolic heterogeneity and local stable regions of chronic hypoxia, signify a tumor's resistance to therapy. Although recent research has shed light on the intracellular mechanisms of cancer metabolic reprogramming, little is known about how tumors become metabolically heterogeneous or chronically hypoxic, namely the initial conditions and spatiotemporal dynamics that drive these cell population conditions. To study these aspects, we developed a minimal, spatially-resolved simulation framework for modeling tissue-scale mixed populations of cells based on diffusible particles the cells consume and release, the concentrations of which determine their behavior in arbitrarily complex ways, and on stochastic reproduction. We simulate cell populations that self-sort to facilitate metabolic symbiosis, that grow according to tumor-stroma signaling patterns, and that give rise to stable local regions of chronic hypoxia near blood vessels. We raise two novel questions in the context of these results: (1) How will two metabolically symbiotic cell subpopulations self-sort in the presence of glucose, oxygen, and lactate gradients? We observe a robust pattern of alternating striations. (2) What is the proper time scale to observe stable local regions of chronic hypoxia? We observe the stability is a function of the balance of three factors related to O2-diffusion rate, local vessel release rate, and viable and hypoxic tumor cell consumption rate. We anticipate our simulation framework will help researchers design better experiments and generate novel hypotheses to better understand dynamic, emergent whole-tumor behavior.
| S-EPMC5193460 | biostudies-literature
Project description:Emergence of tumor resistance to an anti-cancer therapy directed against a putative target raises several questions including: (1) do mutations in the target/pathway confer resistance? (2) Are these mutations pre-existing? (3) What is the relative fitness of cells with/without the mutation? We addressed these questions in patients with metastatic colorectal cancer (mCRC). We conducted an exhaustive review of published data to establish a median doubling time for CRCs and stained a cohort of CRCs to document mitotic indices. We analyzed published data and our own data to calculate rates of growth (g) and regression (d, decay) of tumors in patients with CRC correlating these results with the detection of circulating MT-KRAS DNA. Additionally we estimated mathematically the caloric burden of such tumors using data on mitotic and apoptotic indices. We conclude outgrowth of cells harboring intrinsic or acquired MT-KRAS cannot explain resistance to anti-EGFR (epidermal growth factor receptor) antibodies. Rates of tumor growth with panitumumab are unaffected by presence/absence of MT-KRAS. While MT-KRAS cells may be resistant to anti-EGFR antibodies, WT-KRAS cells also rapidly bypass this blockade suggesting inherent resistance mechanisms are responsible and a neutral evolution model is most appropriate. Using the above clinical data on tumor doubling times and mitotic and apoptotic indices we estimated the caloric intake required to support tumor growth and suggest it may explain in part cancer-associated cachexia.
| S-EPMC5628783 | biostudies-literature
Project description:Combinations of metabolic blockers (incl. fenofibrate) with chemotherapeutic drugs interfere with the drug-resistance of prostate cancer cells. However, their effect on cancer stem cells-dependent microevolution of prostate cancer malignancy remains unaddressed. Here, we hypothesize that the combined docetaxel/fenofibrate treatment prompts the selective expansion of cancer stem cells that affects the microevolution of their progenies. Accordingly, we adapted a combined in vitro/in vivo approach to identify biological and therapeutic consequences of this process. Minute subpopulations of docetaxel-resistant CD133high and/or CD44high cancer stem cell-like (SCL) cells were found in prostate cancer DU145 and PC3 cell populations. When pretreated with docetaxel, they readily differentiated into docetaxel-resistant CD44negative "bulk" cells, thus accounting for the microevolution of drug-resistant cell lineages. Combined docetaxel/fenofibrate treatment induced the generation of poly(morpho)nuclear giant cells and drug-resistant CD44high SCL cells. However, the CD44negative offspring of docetaxel- and docetaxel/fenofibrate-treated SCLs remained relatively sensitive to the combined treatment, while retaining enhanced resistance to docetaxel. Long-term propagation of drug-resistant SCL-derived lineages in the absence of docetaxel/fenofibrate resulted in their reverse microevolution toward the drug-sensitivity and invasive phenotype. Consequently, prostate tumors were able to recover from the combined docetaxel/fenofibrate stress after the initial arrest of their expansion in vivo. In conclusion, we have confirmed the potential of fenofibrate for the metronomic treatment of drug-resistant prostate tumors. However, docetaxel/fenofibrate-induced selective expansion of hyper-resistant CD44high SCL prostate cells and their "bulk" progenies prompts the microevolution of prostate tumor drug-resistance. This process can limit the implementation of metabolic chemotherapy in prostate cancer treatment.
| S-EPMC7756969 | biostudies-literature
Project description:We develop a combined experimental-mathematical framework to investigate heterogeneity in the context of breast cancer treated with doxorubicin. We engineer a cell line to over-express the multi-drug resistance 1 protein (MDR1), an ATP-dependent pump that effluxes intracellular drug. Co-culture experiments mixing the MDR1-overexpressing line with its parental line are evaluated via fluorescence microscopy. To quantify the impact of population heterogeneity on therapy response, these data are analyzed with a coupled pharmacokinetics/pharmacodynamics model. The proliferation and death rates of each line vary with co-culture condition (the relative fraction of each cell line at the time of seeding). For example, the death rate in the parental line under low-dose doxorubicin treatment is increased from 0.64 (± 0.22) × 10-2 to 1.46 (± 0.58) × 10-2 h-1 with increasing fractions of MDR1-overexpressing cells. The growth rate of the MDR1-overexpressing line increases 29% as its relative fraction is decreased. Simulations of the pharmacokinetics/pharmacodynamics model suggest increased efflux from MDR1-overexpressing cells contributes to the increased death rate in the parental cells. Experimentally, the death rate of parental cells is constant across co-culture conditions under co-treatment with an MDR1 inhibitor. These data indicate that intercellular pharmacokinetic variability should be considered in analyzing treatment response in heterogeneous populations.
| S-EPMC5935587 | biostudies-literature
Project description:In the analysis of survival times, the logrank test and the Cox model have been established as key tools, which do not require specific distributional assumptions. Under the assumption of proportional hazards, they are efficient and their results can be interpreted unambiguously. However, delayed treatment effects, disease progression, treatment switchers or the presence of subgroups with differential treatment effects may challenge the assumption of proportional hazards. In practice, weighted logrank tests emphasizing either early, intermediate or late event times via an appropriate weighting function may be used to accommodate for an expected pattern of non-proportionality. We model these sources of non-proportional hazards via a mixture of survival functions with piecewise constant hazard. The model is then applied to study the power of unweighted and weighted log-rank tests, as well as maximum tests allowing different time dependent weights. Simulation results suggest a robust performance of maximum tests across different scenarios, with little loss in power compared to the most powerful among the considered weighting schemes and huge power gain compared to unfavorable weights. The actual sources of non-proportional hazards are not obvious from resulting populationwise survival functions, highlighting the importance of detailed simulations in the planning phase of a trial when assuming non-proportional hazards.We provide the required tools in a software package, allowing to model data generating processes under complex non-proportional hazard scenarios, to simulate data from these models and to perform the weighted logrank tests.
| S-EPMC7818232 | biostudies-literature