Rapid adaptation to CDK2 inhibition exposes intrinsic cell-cycle plasticity.

CDK2 is a core cell-cycle kinase that phosphorylates many substrates to drive progression through the cell cycle. CDK2 is hyperactivated in multiple cancers and is therefore an attractive therapeutic target. Here, we use several CDK2 inhibitors in clinical development to interrogate CDK2 substrate phosphorylation, cell-cycle progression, and drug adaptation in preclinical models. Whereas CDK1 is known to compensate for loss of CDK2 in Cdk2-/- mice, this is not true of acute inhibition of CDK2. Upon CDK2 inhibition, cells exhibit a rapid loss of substrate phosphorylation that rebounds within several hours. CDK4/6 activity backstops inhibition of CDK2 and sustains the proliferative program by maintaining Rb1 hyperphosphorylation, active E2F transcription, and cyclin A2 expression, enabling re-activation of CDK2 in the presence of drug. Our results augment our understanding of CDK plasticity and indicate that co-inhibition of CDK2 and CDK4/6 may be required to suppress adaptation to CDK2 inhibitors currently under clinical assessment.

Ultrasensitivity dynamics of diverse aryl hydrocarbon receptor modulators in a hepatoma cell line.

The aryl hydrocarbon receptor (AhR) is a nuclear receptor that facilitates a wide transcriptional response and causes a variety of adaptive and maladaptive physiological functions. Such functions are entirely dependent on the type of ligand activating it, and therefore, the nuances in the activation of this receptor at the single-cell level have become a research interest for different pharmacological and toxicological applications. Here, we investigate the activation of the AhR by diverse classes of compounds in a Hepa1c1c7-based murine hepatoma cell line. The exogenous compounds analyzed produced different levels of ultrasensitivity in AhR activation as measured by XRE-coupled EGFP production and analyzed by both flow cytometric and computational simulation techniques. Interestingly, simulation experiments reported herein were able to reproduce and quantitate the natural single-cell stochasticity inherent to mammalian cell lines as well as the ligand-specific differences in ultrasensitivity. Classical AhR modulators 2,3,7,8-tetrachlorodibenzodioxin (10- 1-105 pM), PCB-126 (10- 1-107 pM), and benzo[a]pyrene (10- 1-107 pM) produced the greatest levels of single-cell ultrasensitivity and most maximal responses, while consumption-based ligands indole-3-carbinol (103-109 pM), 3,3'-diindolylmethane (103-108 pM), and cannabidiol (103-108 pM) caused low-level AhR activation in more purely graded single-cell fashions. All compounds were tested and analyzed over a 24 h period for consistency. The comparative quantitative results for each compound are presented within. This study aids in defining the disparity between different types of AhR modulators that produce distinctly different physiological outcomes. In addition, the simulation tool developed for this study can be used in future studies to predict the quantitative effects of diverse types of AhR ligands in the context of pharmacological therapies or toxicological concerns.

Reversible Intracellular Gelation of MCF10A Cells Enables Programmable Control Over 3D Spheroid Growth.

In nature, some organisms survive extreme environments by inducing a biostatic state wherein cellular contents are effectively vitrified. Recently, a synthetic biostatic state in mammalian cells is achieved via intracellular network formation using bio-orthogonal strain-promoted azide-alkyne cycloaddition (SPAAC) reactions between functionalized poly(ethylene glycol) (PEG) macromers. In this work, the effects of intracellular network formation on a 3D epithelial MCF10A spheroid model are explored. Macromer-transfected cells are encapsulated in Matrigel, and spheroid area is reduced by ≈50% compared to controls. The intracellular hydrogel network increases the quiescent cell population, as indicated by increased p21 expression. Additionally, bioenergetics (ATP/ADP ratio) and functional metabolic rates are reduced. To enable reversibility of the biostasis effect, a photosensitive nitrobenzyl-containing macromer is incorporated into the PEG network, allowing for light-induced degradation. Following light exposure, cell state, and proliferation return to control levels, while SPAAC-treated spheroids without light exposure (i.e., containing intact intracellular networks) remain smaller and less proliferative through this same period. These results demonstrate that photodegradable intracellular hydrogels can induce a reversible slow-growing state in 3D spheroid culture.

Multiple cancer types rapidly escape from multiple MAPK inhibitors to generate mutagenesis-prone subpopulations.

Many cancers harbor pro-proliferative mutations of the mitogen-activated protein kinase (MAPK) pathway and many targeted inhibitors now exist for clinical use, but drug resistance remains a major issue. We recently showed that BRAF-driven melanoma cells treated with BRAF inhibitors can non-genetically adapt to drug within 3-4 days to escape quiescence and resume slow proliferation. Here we show that this phenomenon is not unique to melanomas treated with BRAF inhibitors but rather is widespread across many clinical MAPK inhibitors and cancer types driven by EGFR, KRAS, and BRAF mutations. In all treatment contexts examined, a subset of cells can escape drug-induced quiescence within four days to resume proliferation. These escapee cells broadly experience aberrant DNA replication, accumulate DNA lesions, spend longer in G2-M cell cycle phases, and mount an ATR-dependent stress response. We further identify the Fanconi anemia (FA) DNA repair pathway as critical for successful mitotic completion in escapees. Long-term cultures, patient samples, and clinical data demonstrate a broad dependency on ATR- and FA-mediated stress tolerance. Together, these results highlight the pervasiveness with which MAPK-mutant cancers are able to rapidly escape drug and the importance of suppressing early stress tolerance pathways to potentially achieve more durable clinical responses to targeted MAPK pathway inhibitors.

Multiple cancers escape from multiple MAPK pathway inhibitors and use DNA replication stress signaling to tolerate aberrant cell cycles.

Many cancers harbor pro-proliferative mutations of the mitogen-activated protein kinase (MAPK) pathway. In BRAF-driven melanoma cells treated with BRAF inhibitors, subpopulations of cells escape drug-induced quiescence through a nongenetic manner of adaptation and resume slow proliferation. Here, we found that this phenomenon is common to many cancer types driven by EGFR, KRAS, or BRAF mutations in response to multiple, clinically approved MAPK pathway inhibitors. In 2D cultures and 3D spheroid models of various cancer cell lines, a subset of cells escaped drug-induced quiescence within 4 days to resume proliferation. These "escapee" cells exhibited DNA replication deficits, accumulated DNA lesions, and mounted a stress response that depended on the ataxia telangiectasia and RAD3-related (ATR) kinase. We further identified that components of the Fanconi anemia (FA) DNA repair pathway are recruited to sites of mitotic DNA synthesis (MiDAS) in escapee cells, enabling successful completion of cell division. Analysis of patient tumor samples and clinical data correlated disease progression with an increase in DNA replication stress response factors. Our findings suggest that many MAPK pathway-mutant cancers rapidly escape drug action and that suppressing early stress tolerance pathways may achieve more durable clinical responses to MAPK pathway inhibitors.

Intracellular Crowding by Bio-Orthogonal Hydrogel Formation Induces Reversible Molecular Stasis.

To survive extreme conditions, certain animals enter a reversible protective stasis through vitrification of the cytosol by polymeric molecules such as proteins and polysaccharides. In this work, synthetic gelation of the cytosol in living cells is used to induce reversible molecular stasis. Through the sequential lipofectamine-mediated transfection of complementary poly(ethylene glycol) macromers into mammalian cells, intracellular crosslinking occurs through bio-orthogonal strain-promoted azide-alkyne cycloaddition click reactions. This achieves efficient polymer uptake with minimal cell death (99% viable). Intracellular crosslinking decreases DNA replication and protein synthesis, and increases the quiescent population by 2.5-fold. Real-time tracking of single cells containing intracellular crosslinked polymers identifies increases in intermitotic time (15 h vs 19 h) and decreases in motility (30 µm h-1  vs 15 µm h-1 ). The cytosol viscosity increases threefold after intracellular crosslinking and results in disordered cytoskeletal structure in addition to the disruption of cellular coordination in a scratch assay. By incorporating photodegradable nitrobenzyl moieties into the polymer backbone, the effects of intracellular crosslinking are reversed upon exposure to light, thereby restoring proliferation (80% phospho-Rb+ cells), protein translation, and migration. Reversible intracellular crosslinking provides a novel method for dynamic manipulation of intracellular mechanics, altering essential processes that determine cellular function.

Melanoma subpopulations that rapidly escape MAPK pathway inhibition incur DNA damage and rely on stress signalling.

Despite the increasing number of effective anti-cancer therapies, successful treatment is limited by the development of drug resistance. While the contribution of genetic factors to drug resistance is undeniable, little is known about how drug-sensitive cells first evade drug action to proliferate in drug. Here we track the responses of thousands of single melanoma cells to BRAF inhibitors and show that a subset of cells escapes drug via non-genetic mechanisms within the first three days of treatment. Cells that escape drug rely on ATF4 stress signalling to cycle periodically in drug, experience DNA replication defects leading to DNA damage, and yet out-proliferate other cells over extended treatment. Together, our work reveals just how rapidly melanoma cells can adapt to drug treatment, generating a mutagenesis-prone subpopulation that expands over time.

Network Simulations Reveal Molecular Signatures of Vulnerability to Age-Dependent Stress and Tau Accumulation.

Alzheimer's disease (AD) is the leading cause of dementia and one of the most common causes of death worldwide. As an age-dependent multifactorial disease, the causative triggers of AD are rooted in spontaneous declines in cellular function and metabolic capacity with increases in protein stressors such as the tau protein. This multitude of age-related processes that cause neurons to change from healthy states to ones vulnerable to the damage seen in AD are difficult to simultaneously investigate and even more difficult to quantify. Here we aimed to diminish these gaps in our understanding of neuronal vulnerability in AD development by using simulation methods to theoretically quantify an array of cellular stress responses and signaling molecules. This temporally-descriptive molecular signature was produced using a novel multimethod simulation approach pioneered by our laboratory for biological research; this methodology combines hierarchical agent-based processes and continuous equation-based modeling in the same interface, all while maintaining intrinsic distributions that emulate natural biological stochasticity. The molecular signature was validated for a normal organismal aging trajectory using experimental longitudinal data from Caenorhabditis elegans and rodent studies. In addition, we have further predicted this aging molecular signature for cells impacted by the pathogenic tau protein, giving rise to distinct stress response conditions needed for cytoprotective aging. Interestingly, our simulation experiments showed that oxidative stress signaling (via daf-16 and skn-1 activities) does not substantially protect cells from all the early stressors of aging, but that it is essential in preventing a late-life degenerative cellular phenotype. Together, our simulation experiments aid in elucidating neurodegenerative triggers in the onset of AD for different genetic conditions. The long-term goal of this work is to provide more detailed diagnostic and prognostic tools for AD development and progression, and to provide more comprehensive preventative measures for this disease.

A multimethod computational simulation approach for investigating mitochondrial dynamics and dysfunction in degenerative aging.

Research in biogerontology has largely focused on the complex relationship between mitochondrial dysfunction and biological aging. In particular, the mitochondrial free radical theory of aging (MFRTA) has been well accepted. However, this theory has been challenged by recent studies showing minimal increases in reactive oxygen species (ROS) as not entirely deleterious in nature, and even beneficial under the appropriate cellular circumstances. To assess these significant and nonintuitive observations in the context of a functional system, we have taken an in silico approach to expand the focus of the MFRTA by including other key mitochondrial stress response pathways, as they have been observed in the nematode Caenorhabditis elegans. These include the mitochondrial unfolded protein response (UPRmt ), mitochondrial biogenesis and autophagy dynamics, the relevant DAF-16 and SKN-1 axes, and NAD+ -dependent deacetylase activities. To integrate these pathways, we have developed a multilevel hybrid-modeling paradigm, containing agent-based elements among stochastic system-dynamics environments of logically derived ordinary differential equations, to simulate aging mitochondrial phenotypes within a population of energetically demanding cells. The simulation experiments resulted in accurate predictions of physiological parameters over time that accompany normal aging, such as the declines in both NAD+ and ATP and an increase in ROS. Additionally, the in silico system was virtually perturbed using a variety of pharmacological (e.g., rapamycin, pterostilbene, paraquat) and genetic (e.g., skn-1, daf-16, sod-2) schemes to quantitate the temporal alterations of specific mechanistic targets, supporting insights into molecular determinants of aging as well as cytoprotective agents that may improve neurological or muscular healthspan.