1
|
Best J, Kim R, Reed M, Nijhout HF. A mathematical model of melatonin synthesis and interactions with the circadian clock. Math Biosci 2024; 377:109280. [PMID: 39243938 DOI: 10.1016/j.mbs.2024.109280] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/17/2024] [Accepted: 08/13/2024] [Indexed: 09/09/2024]
Abstract
A new mathematical model of melatonin synthesis in pineal cells is created and connected to a slightly modified previously created model of the circadian clock in the suprachiasmatic nucleus (SCN). The SCN influences the production of melatonin by upregulating two key enzymes in the pineal. The melatonin produced enters the blood and the cerebrospinal fluid and thus the SCN, influencing the circadian clock. We show that the model of melatonin synthesis corresponds well with extant experimental data and responds similarly to clinical experiments on bright light in the middle of the night. Melatonin is widely used to treat jet lag and sleep disorders. We show how the feedback from the pineal to the SCN causes phase resetting of the circadian clock. Melatonin doses early in the evening advance the clock and doses late at night delay the clock with a dead zone in between where the phase of the clock does not change.
Collapse
Affiliation(s)
- Janet Best
- Department of Mathematics, The Ohio State University, 231 W. 18th Ave., Columbus, 43210, OH, USA.
| | - Ruby Kim
- Department of Mathematics, University of Michigan, 2074 East Hall, 530 Church St., Ann Arbor, 48109, MI, USA
| | - Michael Reed
- Department of Mathematics, Duke University, 120 Science Drive, Campus box 90338, Durham, 27708, NC, USA
| | - H Frederik Nijhout
- Department of Biology, Duke University, Biological Sciences Building, Campus box 90320, Durham, 27708, NC, USA
| |
Collapse
|
2
|
Helenek C, Krzysztoń R, Petreczky J, Wan Y, Cabral M, Coraci D, Balázsi G. Synthetic gene circuit evolution: Insights and opportunities at the mid-scale. Cell Chem Biol 2024; 31:1447-1459. [PMID: 38925113 PMCID: PMC11330362 DOI: 10.1016/j.chembiol.2024.05.018] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/12/2024] [Revised: 05/07/2024] [Accepted: 05/30/2024] [Indexed: 06/28/2024]
Abstract
Directed evolution focuses on optimizing single genetic components for predefined engineering goals by artificial mutagenesis and selection. In contrast, experimental evolution studies the adaptation of entire genomes in serially propagated cell populations, to provide an experimental basis for evolutionary theory. There is a relatively unexplored gap at the middle ground between these two techniques, to evolve in vivo entire synthetic gene circuits with nontrivial dynamic function instead of single parts or whole genomes. We discuss the requirements for such mid-scale evolution, with hypothetical examples for evolving synthetic gene circuits by appropriate selection and targeted shuffling of a seed set of genetic components in vivo. Implementing similar methods should aid the rapid generation, functionalization, and optimization of synthetic gene circuits in various organisms and environments, accelerating both the development of biomedical and technological applications and the understanding of principles guiding regulatory network evolution.
Collapse
Affiliation(s)
- Christopher Helenek
- The Louis and Beatrice Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, NY 11794, USA; Department of Biomedical Engineering, Stony Brook University, Stony Brook, NY 11794, USA
| | - Rafał Krzysztoń
- The Louis and Beatrice Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, NY 11794, USA
| | - Julia Petreczky
- The Louis and Beatrice Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, NY 11794, USA; Department of Chemistry, Stony Brook University, Stony Brook, NY 11794, USA
| | - Yiming Wan
- The Louis and Beatrice Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, NY 11794, USA
| | - Mariana Cabral
- The Louis and Beatrice Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, NY 11794, USA; Department of Biomedical Engineering, Stony Brook University, Stony Brook, NY 11794, USA
| | - Damiano Coraci
- The Louis and Beatrice Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, NY 11794, USA
| | - Gábor Balázsi
- The Louis and Beatrice Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, NY 11794, USA; Department of Biomedical Engineering, Stony Brook University, Stony Brook, NY 11794, USA; Stony Brook Cancer Center, Stony Brook University, Stony Brook, NY 11794, USA.
| |
Collapse
|
3
|
Vaidehi Narayanan H, Xiang MY, Chen Y, Huang H, Roy S, Makkar H, Hoffmann A, Roy K. Direct observation correlates NFκB cRel in B cells with activating and terminating their proliferative program. Proc Natl Acad Sci U S A 2024; 121:e2309686121. [PMID: 39024115 PMCID: PMC11287273 DOI: 10.1073/pnas.2309686121] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/19/2023] [Accepted: 05/28/2024] [Indexed: 07/20/2024] Open
Abstract
Antibody responses require the proliferative expansion of B cells controlled by affinity-dependent signals. Yet, proliferative bursts are heterogeneous, varying between 0 and 8 divisions in response to the same stimulus. NFκB cRel is activated in response to immune stimulation in B cells and is genetically required for proliferation. Here, we asked whether proliferative heterogeneity is controlled by natural variations in cRel abundance. We developed a fluorescent reporter mTFP1-cRel for the direct observation of cRel in live proliferating B cells. We found that cRel is heterogeneously distributed among naïve B cells, which are enriched for high expressors in a heavy-tailed distribution. We found that high cRel expressors show faster activation of the proliferative program, but do not sustain it well, with population expansion decaying earlier. With a mathematical model of the molecular network, we showed that cRel heterogeneity arises from balancing positive feedback by autoregulation and negative feedback by its inhibitor IκBε, confirmed by mouse knockouts. Using live-cell fluorescence microscopy, we showed that increased cRel primes B cells for early proliferation via higher basal expression of the cell cycle driver cMyc. However, peak cMyc induction amplitude is constrained by incoherent feedforward regulation, decoding the fold change of cRel activity to terminate the proliferative burst. This results in a complex nonlinear, nonmonotonic relationship between cRel expression and the extent of proliferation. These findings emphasize the importance of direct observational studies to complement gene knockout results and to learn about quantitative relationships between biological processes and their key regulators in the context of natural variations.
Collapse
Affiliation(s)
- Haripriya Vaidehi Narayanan
- Signaling Systems Laboratory, Department of Microbiology, Immunology, and Molecular Genetics, University of California Los Angeles, Los Angeles, CA90095
- Institute for Quantitative and Computational Biosciences, University of California Los Angeles, Los Angeles, CA90095
| | - Mark Y. Xiang
- Signaling Systems Laboratory, Department of Microbiology, Immunology, and Molecular Genetics, University of California Los Angeles, Los Angeles, CA90095
- Institute for Quantitative and Computational Biosciences, University of California Los Angeles, Los Angeles, CA90095
| | - Yijia Chen
- Signaling Systems Laboratory, Department of Microbiology, Immunology, and Molecular Genetics, University of California Los Angeles, Los Angeles, CA90095
- Institute for Quantitative and Computational Biosciences, University of California Los Angeles, Los Angeles, CA90095
| | - Helen Huang
- Signaling Systems Laboratory, Department of Microbiology, Immunology, and Molecular Genetics, University of California Los Angeles, Los Angeles, CA90095
- Institute for Quantitative and Computational Biosciences, University of California Los Angeles, Los Angeles, CA90095
| | - Sukanya Roy
- Division of Microbiology and Immunology, Department of Pathology, University of Utah, Salt Lake City, UT84112
| | - Himani Makkar
- Division of Microbiology and Immunology, Department of Pathology, University of Utah, Salt Lake City, UT84112
| | - Alexander Hoffmann
- Signaling Systems Laboratory, Department of Microbiology, Immunology, and Molecular Genetics, University of California Los Angeles, Los Angeles, CA90095
- Institute for Quantitative and Computational Biosciences, University of California Los Angeles, Los Angeles, CA90095
| | - Koushik Roy
- Division of Microbiology and Immunology, Department of Pathology, University of Utah, Salt Lake City, UT84112
| |
Collapse
|
4
|
Brückner DB, Broedersz CP. Learning dynamical models of single and collective cell migration: a review. REPORTS ON PROGRESS IN PHYSICS. PHYSICAL SOCIETY (GREAT BRITAIN) 2024; 87:056601. [PMID: 38518358 DOI: 10.1088/1361-6633/ad36d2] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/07/2023] [Accepted: 03/22/2024] [Indexed: 03/24/2024]
Abstract
Single and collective cell migration are fundamental processes critical for physiological phenomena ranging from embryonic development and immune response to wound healing and cancer metastasis. To understand cell migration from a physical perspective, a broad variety of models for the underlying physical mechanisms that govern cell motility have been developed. A key challenge in the development of such models is how to connect them to experimental observations, which often exhibit complex stochastic behaviours. In this review, we discuss recent advances in data-driven theoretical approaches that directly connect with experimental data to infer dynamical models of stochastic cell migration. Leveraging advances in nanofabrication, image analysis, and tracking technology, experimental studies now provide unprecedented large datasets on cellular dynamics. In parallel, theoretical efforts have been directed towards integrating such datasets into physical models from the single cell to the tissue scale with the aim of conceptualising the emergent behaviour of cells. We first review how this inference problem has been addressed in both freely migrating and confined cells. Next, we discuss why these dynamics typically take the form of underdamped stochastic equations of motion, and how such equations can be inferred from data. We then review applications of data-driven inference and machine learning approaches to heterogeneity in cell behaviour, subcellular degrees of freedom, and to the collective dynamics of multicellular systems. Across these applications, we emphasise how data-driven methods can be integrated with physical active matter models of migrating cells, and help reveal how underlying molecular mechanisms control cell behaviour. Together, these data-driven approaches are a promising avenue for building physical models of cell migration directly from experimental data, and for providing conceptual links between different length-scales of description.
Collapse
Affiliation(s)
- David B Brückner
- Institute of Science and Technology Austria, Am Campus 1, 3400 Klosterneuburg, Austria
| | - Chase P Broedersz
- Department of Physics and Astronomy, Vrije Universiteit Amsterdam, 1081 HV Amsterdam, The Netherlands
- Arnold Sommerfeld Center for Theoretical Physics and Center for NanoScience, Department of Physics, Ludwig-Maximilian-University Munich, Theresienstr. 37, D-80333 Munich, Germany
| |
Collapse
|
5
|
Nano M, Montell DJ. Apoptotic signaling: Beyond cell death. Semin Cell Dev Biol 2024; 156:22-34. [PMID: 37988794 DOI: 10.1016/j.semcdb.2023.11.002] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/26/2023] [Revised: 11/02/2023] [Accepted: 11/04/2023] [Indexed: 11/23/2023]
Abstract
Apoptosis is the best described form of regulated cell death, and was, until relatively recently, considered irreversible once particular biochemical points-of-no-return were activated. In this manuscript, we examine the mechanisms cells use to escape from a self-amplifying death signaling module. We discuss the role of feedback, dynamics, propagation, and noise in apoptotic signaling. We conclude with a revised model for the role of apoptosis in animal development, homeostasis, and disease.
Collapse
Affiliation(s)
- Maddalena Nano
- Molecular, Cellular, and Developmental Biology Department, University of California, Santa Barbara, CA 93106, USA; Neuroscience Research Institute, University of California, Santa Barbara, CA 93106, USA.
| | - Denise J Montell
- Molecular, Cellular, and Developmental Biology Department, University of California, Santa Barbara, CA 93106, USA; Neuroscience Research Institute, University of California, Santa Barbara, CA 93106, USA.
| |
Collapse
|
6
|
Schindler-Johnson M, Petridou NI. Collective effects of cell cleavage dynamics. Front Cell Dev Biol 2024; 12:1358971. [PMID: 38559810 PMCID: PMC10978805 DOI: 10.3389/fcell.2024.1358971] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/20/2023] [Accepted: 03/05/2024] [Indexed: 04/04/2024] Open
Abstract
A conserved process of early embryonic development in metazoans is the reductive cell divisions following oocyte fertilization, termed cell cleavages. Cell cleavage cycles usually start synchronously, lengthen differentially between the embryonic cells becoming asynchronous, and cease before major morphogenetic events, such as germ layer formation and gastrulation. Despite exhibiting species-specific characteristics, the regulation of cell cleavage dynamics comes down to common controllers acting mostly at the single cell/nucleus level, such as nucleus-to-cytoplasmic ratio and zygotic genome activation. Remarkably, recent work has linked cell cleavage dynamics to the emergence of collective behavior during embryogenesis, including pattern formation and changes in embryo-scale mechanics, raising the question how single-cell controllers coordinate embryo-scale processes. In this review, we summarize studies across species where an association between cell cleavages and collective behavior was made, discuss the underlying mechanisms, and propose that cell-to-cell variability in cell cleavage dynamics can serve as a mechanism of long-range coordination in developing embryos.
Collapse
Affiliation(s)
- Magdalena Schindler-Johnson
- Developmental Biology Unit, European Molecular Biology Laboratory, Heidelberg, Germany
- Faculty of Biosciences, Heidelberg University, Heidelberg, Germany
| | - Nicoletta I. Petridou
- Developmental Biology Unit, European Molecular Biology Laboratory, Heidelberg, Germany
| |
Collapse
|
7
|
Fourneaux C, Racine L, Koering C, Dussurgey S, Vallin E, Moussy A, Parmentier R, Brunard F, Stockholm D, Modolo L, Picard F, Gandrillon O, Paldi A, Gonin-Giraud S. Differentiation is accompanied by a progressive loss in transcriptional memory. BMC Biol 2024; 22:58. [PMID: 38468285 DOI: 10.1186/s12915-024-01846-9] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/05/2023] [Accepted: 02/13/2024] [Indexed: 03/13/2024] Open
Abstract
BACKGROUND Cell differentiation requires the integration of two opposite processes, a stabilizing cellular memory, especially at the transcriptional scale, and a burst of gene expression variability which follows the differentiation induction. Therefore, the actual capacity of a cell to undergo phenotypic change during a differentiation process relies upon a modification in this balance which favors change-inducing gene expression variability. However, there are no experimental data providing insight on how fast the transcriptomes of identical cells would diverge on the scale of the very first two cell divisions during the differentiation process. RESULTS In order to quantitatively address this question, we developed different experimental methods to recover the transcriptomes of related cells, after one and two divisions, while preserving the information about their lineage at the scale of a single cell division. We analyzed the transcriptomes of related cells from two differentiation biological systems (human CD34+ cells and T2EC chicken primary erythrocytic progenitors) using two different single-cell transcriptomics technologies (scRT-qPCR and scRNA-seq). CONCLUSIONS We identified that the gene transcription profiles of differentiating sister cells are more similar to each other than to those of non-related cells of the same type, sharing the same environment and undergoing similar biological processes. More importantly, we observed greater discrepancies between differentiating sister cells than between self-renewing sister cells. Furthermore, a progressive increase in this divergence from first generation to second generation was observed when comparing differentiating cousin cells to self renewing cousin cells. Our results are in favor of a gradual erasure of transcriptional memory during the differentiation process.
Collapse
Affiliation(s)
- Camille Fourneaux
- Laboratoire de Biologie et Modélisation de la Cellule, Ecole Normale Supérieure de Lyon, CNRS, UMR5239, Université Claude Bernard Lyon 1, Lyon, France
| | - Laëtitia Racine
- Ecole Pratique des Hautes Etudes, PSL Research University, Sorbonne Université, INSERM, CRSA, Paris, 75012, France
| | - Catherine Koering
- Laboratoire de Biologie et Modélisation de la Cellule, Ecole Normale Supérieure de Lyon, CNRS, UMR5239, Université Claude Bernard Lyon 1, Lyon, France
| | - Sébastien Dussurgey
- Plateforme AniRA-Cytométrie, Université Claude Bernard Lyon 1, CNRS UAR3444, Inserm US8, ENS de Lyon, SFR Biosciences, Lyon, F-69007, France
| | - Elodie Vallin
- Laboratoire de Biologie et Modélisation de la Cellule, Ecole Normale Supérieure de Lyon, CNRS, UMR5239, Université Claude Bernard Lyon 1, Lyon, France
| | - Alice Moussy
- Ecole Pratique des Hautes Etudes, PSL Research University, Sorbonne Université, INSERM, CRSA, Paris, 75012, France
| | - Romuald Parmentier
- Ecole Pratique des Hautes Etudes, PSL Research University, Sorbonne Université, INSERM, CRSA, Paris, 75012, France
| | - Fanny Brunard
- Laboratoire de Biologie et Modélisation de la Cellule, Ecole Normale Supérieure de Lyon, CNRS, UMR5239, Université Claude Bernard Lyon 1, Lyon, France
| | - Daniel Stockholm
- Ecole Pratique des Hautes Etudes, PSL Research University, Sorbonne Université, INSERM, CRSA, Paris, 75012, France
| | - Laurent Modolo
- Laboratoire de Biologie et Modélisation de la Cellule, Ecole Normale Supérieure de Lyon, CNRS, UMR5239, Université Claude Bernard Lyon 1, Lyon, France
| | - Franck Picard
- Laboratoire de Biologie et Modélisation de la Cellule, Ecole Normale Supérieure de Lyon, CNRS, UMR5239, Université Claude Bernard Lyon 1, Lyon, France
| | - Olivier Gandrillon
- Laboratoire de Biologie et Modélisation de la Cellule, Ecole Normale Supérieure de Lyon, CNRS, UMR5239, Université Claude Bernard Lyon 1, Lyon, France
- Inria Center, Grenoble Rhone-Alpes, Equipe Dracula, Villeurbanne, F69100, France
| | - Andras Paldi
- Ecole Pratique des Hautes Etudes, PSL Research University, Sorbonne Université, INSERM, CRSA, Paris, 75012, France
| | - Sandrine Gonin-Giraud
- Laboratoire de Biologie et Modélisation de la Cellule, Ecole Normale Supérieure de Lyon, CNRS, UMR5239, Université Claude Bernard Lyon 1, Lyon, France.
| |
Collapse
|
8
|
Mold JE, Weissman MH, Ratz M, Hagemann-Jensen M, Hård J, Eriksson CJ, Toosi H, Berghenstråhle J, Ziegenhain C, von Berlin L, Martin M, Blom K, Lagergren J, Lundeberg J, Sandberg R, Michaëlsson J, Frisén J. Clonally heritable gene expression imparts a layer of diversity within cell types. Cell Syst 2024; 15:149-165.e10. [PMID: 38340731 DOI: 10.1016/j.cels.2024.01.004] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/24/2022] [Revised: 05/25/2023] [Accepted: 01/18/2024] [Indexed: 02/12/2024]
Abstract
Cell types can be classified according to shared patterns of transcription. Non-genetic variability among individual cells of the same type has been ascribed to stochastic transcriptional bursting and transient cell states. Using high-coverage single-cell RNA profiling, we asked whether long-term, heritable differences in gene expression can impart diversity within cells of the same type. Studying clonal human lymphocytes and mouse brain cells, we uncovered a vast diversity of heritable gene expression patterns among different clones of cells of the same type in vivo. We combined chromatin accessibility and RNA profiling on different lymphocyte clones to reveal thousands of regulatory regions exhibiting interclonal variation, which could be directly linked to interclonal variation in gene expression. Our findings identify a source of cellular diversity, which may have important implications for how cellular populations are shaped by selective processes in development, aging, and disease. A record of this paper's transparent peer review process is included in the supplemental information.
Collapse
Affiliation(s)
- Jeff E Mold
- Department of Cell and Molecular Biology, Karolinska Institute, Stockholm, Sweden
| | - Martin H Weissman
- Mathematics Department, University of California, Santa Cruz, CA, USA.
| | - Michael Ratz
- Department of Cell and Molecular Biology, Karolinska Institute, Stockholm, Sweden; SciLifeLab, Department of Gene Technology, KTH Royal Institute of Technology, Stockholm, Sweden
| | | | - Joanna Hård
- Department of Cell and Molecular Biology, Karolinska Institute, Stockholm, Sweden
| | - Carl-Johan Eriksson
- Department of Cell and Molecular Biology, Karolinska Institute, Stockholm, Sweden
| | - Hosein Toosi
- SciLifeLab, Computational Science and Technology Department, KTH Royal Institute of Technology, Stockholm, Sweden
| | - Joseph Berghenstråhle
- SciLifeLab, Department of Gene Technology, KTH Royal Institute of Technology, Stockholm, Sweden
| | - Christoph Ziegenhain
- Department of Cell and Molecular Biology, Karolinska Institute, Stockholm, Sweden
| | - Leonie von Berlin
- Department of Cell and Molecular Biology, Karolinska Institute, Stockholm, Sweden
| | - Marcel Martin
- Department of Biochemistry and Biophysics, National Bioinformatics Infrastructure Sweden, SciLifeLab, Stockholm University, Stockholm, Sweden
| | - Kim Blom
- Center for Infectious Medicine, Department of Medicine Huddinge, Karolinska Institute, Stockholm, Sweden
| | - Jens Lagergren
- SciLifeLab, Computational Science and Technology Department, KTH Royal Institute of Technology, Stockholm, Sweden
| | - Joakim Lundeberg
- SciLifeLab, Department of Gene Technology, KTH Royal Institute of Technology, Stockholm, Sweden
| | - Rickard Sandberg
- Department of Cell and Molecular Biology, Karolinska Institute, Stockholm, Sweden
| | - Jakob Michaëlsson
- Center for Infectious Medicine, Department of Medicine Huddinge, Karolinska Institute, Stockholm, Sweden.
| | - Jonas Frisén
- Department of Cell and Molecular Biology, Karolinska Institute, Stockholm, Sweden
| |
Collapse
|
9
|
Tserunyan V, Finley S. Information-Theoretic Analysis of a Model of CAR-4-1BB-Mediated NFκB Activation. Bull Math Biol 2023; 86:5. [PMID: 38038772 PMCID: PMC10691998 DOI: 10.1007/s11538-023-01232-6] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/09/2023] [Accepted: 11/01/2023] [Indexed: 12/02/2023]
Abstract
Systems biology utilizes computational approaches to examine an array of biological processes, such as cell signaling, metabolomics and pharmacology. This includes mathematical modeling of CAR T cells, a modality of cancer therapy by which genetically engineered immune cells recognize and combat a cancerous target. While successful against hematologic malignancies, CAR T cells have shown limited success against other cancer types. Thus, more research is needed to understand their mechanisms of action and leverage their full potential. In our work, we set out to apply information theory on a mathematical model of NFκB signaling initiated by the CAR following antigen encounter. First, we estimated channel capacity for CAR-4-1BB-mediated NFκB signal transduction. Next, we evaluated the pathway's ability to distinguish contrasting "low" and "high" antigen concentration levels, depending on the amount of variability in protein concentrations. Finally, we assessed the fidelity by which NFκB activation reflects the encountered antigen concentration, depending on the prevalence of antigen-positive targets in tumor population. We found that in most scenarios, fold change in the nuclear concentration of NFκB carries a higher channel capacity for the pathway than NFκB's absolute response. Additionally, we found that most errors in transducing the antigen signal through the pathway skew towards underestimating the concentration of encountered antigen. Finally, we found that disabling IKKβ deactivation could increase signaling fidelity against targets with antigen-negative cells. Our information-theoretic analysis of signal transduction can provide novel perspectives on biological signaling, as well as enable a more informed path to cell engineering.Kindly check and confirm whether the corresponding affiliation is correctly identified.this is correct.
Collapse
Affiliation(s)
- Vardges Tserunyan
- Department of Quantitative and Computational Biology, University of Southern California, Los Angeles, CA, USA
| | - Stacey Finley
- Department of Quantitative and Computational Biology, University of Southern California, Los Angeles, CA, USA.
- Alfred E. Mann Department of Biomedical Engineering, University of Southern California, Los Angeles, CA, USA.
- Mork Family Department of Chemical Engineering and Materials Science, University of Southern California, Los Angeles, CA, USA.
| |
Collapse
|
10
|
Rajakaruna H, Desai M, Das J. PASCAR: a multiscale framework to explore the design space of constitutive and inducible CAR T cells. Life Sci Alliance 2023; 6:e202302171. [PMID: 37507138 PMCID: PMC10387492 DOI: 10.26508/lsa.202302171] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/18/2023] [Revised: 07/08/2023] [Accepted: 07/12/2023] [Indexed: 07/30/2023] Open
Abstract
CAR T cells are engineered to bind and destroy tumor cells by targeting overexpressed surface antigens. However, healthy cells expressing lower abundances of these antigens can also be lysed by CAR T cells. Various CAR T cell designs increase tumor cell elimination, whereas reducing damage to healthy cells. However, these efforts are costly and labor-intensive, constraining systematic exploration of potential hypotheses. We develop a protein abundance structured population dynamic model for CAR T cells (PASCAR), a framework that combines multiscale population dynamic models and multi-objective optimization approaches with data from cytometry and cytotoxicity assays to systematically explore the design space of constitutive and tunable CAR T cells. PASCAR can quantitatively describe in vitro and in vivo results for constitutive and inducible CAR T cells and can successfully predict experiments outside the training data. Our exploration of the CAR design space reveals that optimal CAR affinities in the intermediate range of dissociation constants effectively reduce healthy cell lysis, whereas maintaining high tumor cell-killing rates. Furthermore, our modeling offers guidance for optimizing CAR expressions in synthetic notch CAR T cells. PASCAR can be extended to other CAR immune cells.
Collapse
Affiliation(s)
- Harshana Rajakaruna
- The Steve and Cindy Rasmussen Institute for Genomics, The Abigail Wexner Research Institute, Nationwide Children's Hospital, Columbus, OH, USA
| | - Milie Desai
- Department of Biology, Indian Institute of Science Education and Research, Pune, India
| | - Jayajit Das
- The Steve and Cindy Rasmussen Institute for Genomics, The Abigail Wexner Research Institute, Nationwide Children's Hospital, Columbus, OH, USA
- Department of Pediatrics and Pelotonia Institute for Immuno-Oncology, College of Medicine, Columbus, OH, USA
- Biophysics Program, The Ohio State University, Columbus, OH, USA
| |
Collapse
|
11
|
Zhang J, Han X, Ma L, Xu S, Lin Y. Deciphering a global source of non-genetic heterogeneity in cancer cells. Nucleic Acids Res 2023; 51:9019-9038. [PMID: 37587722 PMCID: PMC10516630 DOI: 10.1093/nar/gkad666] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/13/2022] [Revised: 07/09/2023] [Accepted: 08/09/2023] [Indexed: 08/18/2023] Open
Abstract
Cell-to-cell variability within a clonal population, also known as non-genetic heterogeneity, has created significant challenges for intervening with diseases such as cancer. While non-genetic heterogeneity can arise from the variability in the expression of specific genes, it remains largely unclear whether and how clonal cells could be heterogeneous in the expression of the entire transcriptome. Here, we showed that gene transcriptional activity is globally modulated in individual cancer cells, leading to non-genetic heterogeneity in the global transcription rate. Such heterogeneity contributes to cell-to-cell variability in transcriptome size and displays both dynamic and static characteristics, with the global transcription rate temporally modulated in a cell-cycle-coupled manner and the time-averaged rate being distinct between cells and heritable across generations. Additional evidence indicated the role of ATP metabolism in this heterogeneity, and suggested its implication in intrinsic cancer drug tolerance. Collectively, our work shed light on the mode, mechanism, and implication of a global but often hidden source of non-genetic heterogeneity.
Collapse
Affiliation(s)
- Jianhan Zhang
- Center for Quantitative Biology and Peking-Tsinghua Center for Life Sciences, Academy for Advanced Interdisciplinary Studies, Peking University, Beijing 100871, China
- The MOE Key Laboratory of Cell Proliferation and Differentiation, School of Life Sciences, Peking University, Beijing 100871, China
| | - Xu Han
- Center for Quantitative Biology and Peking-Tsinghua Center for Life Sciences, Academy for Advanced Interdisciplinary Studies, Peking University, Beijing 100871, China
- The MOE Key Laboratory of Cell Proliferation and Differentiation, School of Life Sciences, Peking University, Beijing 100871, China
| | - Liang Ma
- Center for Quantitative Biology and Peking-Tsinghua Center for Life Sciences, Academy for Advanced Interdisciplinary Studies, Peking University, Beijing 100871, China
- The MOE Key Laboratory of Cell Proliferation and Differentiation, School of Life Sciences, Peking University, Beijing 100871, China
| | - Shuhui Xu
- Center for Quantitative Biology and Peking-Tsinghua Center for Life Sciences, Academy for Advanced Interdisciplinary Studies, Peking University, Beijing 100871, China
- The MOE Key Laboratory of Cell Proliferation and Differentiation, School of Life Sciences, Peking University, Beijing 100871, China
| | - Yihan Lin
- Center for Quantitative Biology and Peking-Tsinghua Center for Life Sciences, Academy for Advanced Interdisciplinary Studies, Peking University, Beijing 100871, China
- The MOE Key Laboratory of Cell Proliferation and Differentiation, School of Life Sciences, Peking University, Beijing 100871, China
| |
Collapse
|
12
|
Grabowski F, Nałęcz-Jawecki P, Lipniacki T. Predictive power of non-identifiable models. Sci Rep 2023; 13:11143. [PMID: 37429934 DOI: 10.1038/s41598-023-37939-8] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/07/2023] [Accepted: 06/29/2023] [Indexed: 07/12/2023] Open
Abstract
Resolving practical non-identifiability of computational models typically requires either additional data or non-algorithmic model reduction, which frequently results in models containing parameters lacking direct interpretation. Here, instead of reducing models, we explore an alternative, Bayesian approach, and quantify the predictive power of non-identifiable models. We considered an example biochemical signalling cascade model as well as its mechanical analogue. For these models, we demonstrated that by measuring a single variable in response to a properly chosen stimulation protocol, the dimensionality of the parameter space is reduced, which allows for predicting the measured variable's trajectory in response to different stimulation protocols even if all model parameters remain unidentified. Moreover, one can predict how such a trajectory will transform in the case of a multiplicative change of an arbitrary model parameter. Successive measurements of remaining variables further reduce the dimensionality of the parameter space and enable new predictions. We analysed potential pitfalls of the proposed approach that can arise when the investigated model is oversimplified, incorrect, or when the training protocol is inadequate. The main advantage of the suggested iterative approach is that the predictive power of the model can be assessed and practically utilised at each step.
Collapse
Affiliation(s)
- Frederic Grabowski
- Institute of Fundamental Technological Research, Polish Academy of Sciences, Warsaw, Poland
| | - Paweł Nałęcz-Jawecki
- Institute of Fundamental Technological Research, Polish Academy of Sciences, Warsaw, Poland
| | - Tomasz Lipniacki
- Institute of Fundamental Technological Research, Polish Academy of Sciences, Warsaw, Poland.
| |
Collapse
|
13
|
Ramu A, Cohen BA. Transcription factor fluctuations underlie cell-to-cell variability in a signaling pathway response. Genetics 2023; 224:iyad094. [PMID: 37226217 PMCID: PMC10691749 DOI: 10.1093/genetics/iyad094] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/08/2023] [Revised: 04/25/2023] [Accepted: 05/05/2023] [Indexed: 05/26/2023] Open
Abstract
Stochastic differences among clonal cells can initiate cell fate decisions in development or cause cell-to-cell differences in the responses to drugs or extracellular ligands. One hypothesis is that some of this phenotypic variability is caused by stochastic fluctuations in the activities of transcription factors (TFs). We tested this hypothesis in NIH3T3-CG cells using the response to Hedgehog signaling as a model cellular response. Here, we present evidence for the existence of distinct fast- and slow-responding substates in NIH3T3-CG cells. These two substates have distinct expression profiles, and fluctuations in the Prrx1 TF underlie some of the differences in expression and responsiveness between fast and slow cells. Our results show that fluctuations in TFs can contribute to cell-to-cell differences in Hedgehog signaling.
Collapse
Affiliation(s)
- Avinash Ramu
- The Edison Family Center for Genome Sciences and Systems Biology, School of Medicine, Washington University in St. Louis, Saint Louis, MO 63110, USA
- Department of Genetics, School of Medicine, Washington University in St. Louis, Saint Louis, MO 63110, USA
| | - Barak A Cohen
- The Edison Family Center for Genome Sciences and Systems Biology, School of Medicine, Washington University in St. Louis, Saint Louis, MO 63110, USA
- Department of Genetics, School of Medicine, Washington University in St. Louis, Saint Louis, MO 63110, USA
| |
Collapse
|
14
|
Schneider M, Bird AD, Gidon A, Triesch J, Jedlicka P, Cuntz H. Biological complexity facilitates tuning of the neuronal parameter space. PLoS Comput Biol 2023; 19:e1011212. [PMID: 37399220 DOI: 10.1371/journal.pcbi.1011212] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/16/2023] [Accepted: 05/24/2023] [Indexed: 07/05/2023] Open
Abstract
The electrical and computational properties of neurons in our brains are determined by a rich repertoire of membrane-spanning ion channels and elaborate dendritic trees. However, the precise reason for this inherent complexity remains unknown, given that simpler models with fewer ion channels are also able to functionally reproduce the behaviour of some neurons. Here, we stochastically varied the ion channel densities of a biophysically detailed dentate gyrus granule cell model to produce a large population of putative granule cells, comparing those with all 15 original ion channels to their reduced but functional counterparts containing only 5 ion channels. Strikingly, valid parameter combinations in the full models were dramatically more frequent at -6% vs. -1% in the simpler model. The full models were also more stable in the face of perturbations to channel expression levels. Scaling up the numbers of ion channels artificially in the reduced models recovered these advantages confirming the key contribution of the actual number of ion channel types. We conclude that the diversity of ion channels gives a neuron greater flexibility and robustness to achieve a target excitability.
Collapse
Affiliation(s)
- Marius Schneider
- Frankfurt Institute for Advanced Studies, Frankfurt am Main, Germany
- Ernst Strüngmann Institute (ESI) for Neuroscience in cooperation with the Max Planck Society, Frankfurt am Main, Germany
- ICAR3R-Interdisciplinary Centre for 3Rs in Animal Research, Justus Liebig University Giessen, Giessen, Germany
- Faculty of Physics, Goethe University, Frankfurt/Main, Frankfurt am Main, Germany
| | - Alexander D Bird
- Frankfurt Institute for Advanced Studies, Frankfurt am Main, Germany
- Ernst Strüngmann Institute (ESI) for Neuroscience in cooperation with the Max Planck Society, Frankfurt am Main, Germany
- ICAR3R-Interdisciplinary Centre for 3Rs in Animal Research, Justus Liebig University Giessen, Giessen, Germany
| | - Albert Gidon
- Institute for Biology, Humboldt-Universität zu Berlin, Berlin, Germany
| | - Jochen Triesch
- Frankfurt Institute for Advanced Studies, Frankfurt am Main, Germany
- Faculty of Physics, Goethe University, Frankfurt/Main, Frankfurt am Main, Germany
- Faculty of Computer Science and Mathematics, Goethe University, Frankfurt am Main, Germany
| | - Peter Jedlicka
- Frankfurt Institute for Advanced Studies, Frankfurt am Main, Germany
- ICAR3R-Interdisciplinary Centre for 3Rs in Animal Research, Justus Liebig University Giessen, Giessen, Germany
- Institute of Clinical Neuroanatomy, Neuroscience Center, Goethe University, Frankfurt am Main, Germany
| | - Hermann Cuntz
- Frankfurt Institute for Advanced Studies, Frankfurt am Main, Germany
- Ernst Strüngmann Institute (ESI) for Neuroscience in cooperation with the Max Planck Society, Frankfurt am Main, Germany
- ICAR3R-Interdisciplinary Centre for 3Rs in Animal Research, Justus Liebig University Giessen, Giessen, Germany
| |
Collapse
|
15
|
Wan Y, Cohen J, Szenk M, Farquhar KS, Coraci D, Krzysztoń R, Azukas J, Van Nest N, Smashnov A, Chern YJ, De Martino D, Nguyen LC, Bien H, Bravo-Cordero JJ, Chan CH, Rosner MR, Balázsi G. Nonmonotone invasion landscape by noise-aware control of metastasis activator levels. Nat Chem Biol 2023; 19:887-899. [PMID: 37231268 PMCID: PMC10299915 DOI: 10.1038/s41589-023-01344-z] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/12/2022] [Accepted: 04/18/2023] [Indexed: 05/27/2023]
Abstract
A major pharmacological assumption is that lowering disease-promoting protein levels is generally beneficial. For example, inhibiting metastasis activator BACH1 is proposed to decrease cancer metastases. Testing such assumptions requires approaches to measure disease phenotypes while precisely adjusting disease-promoting protein levels. Here we developed a two-step strategy to integrate protein-level tuning, noise-aware synthetic gene circuits into a well-defined human genomic safe harbor locus. Unexpectedly, engineered MDA-MB-231 metastatic human breast cancer cells become more, then less and then more invasive as we tune BACH1 levels up, irrespective of the native BACH1. BACH1 expression shifts in invading cells, and expression of BACH1's transcriptional targets confirm BACH1's nonmonotone phenotypic and regulatory effects. Thus, chemical inhibition of BACH1 could have unwanted effects on invasion. Additionally, BACH1's expression variability aids invasion at high BACH1 expression. Overall, precisely engineered, noise-aware protein-level control is necessary and important to unravel disease effects of genes to improve clinical drug efficacy.
Collapse
Affiliation(s)
- Yiming Wan
- Department of Biomedical Engineering, Stony Brook University, Stony Brook, NY, USA
- Louis and Beatrice Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, NY, USA
| | - Joseph Cohen
- Department of Biomedical Engineering, Stony Brook University, Stony Brook, NY, USA
- Louis and Beatrice Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, NY, USA
| | - Mariola Szenk
- Department of Biomedical Engineering, Stony Brook University, Stony Brook, NY, USA
- Louis and Beatrice Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, NY, USA
| | - Kevin S Farquhar
- Department of Biomedical Engineering, Stony Brook University, Stony Brook, NY, USA
- Louis and Beatrice Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, NY, USA
- Genetics and Epigenetics Graduate Program, The University of Texas MD Anderson Cancer Center, UT Health Graduate School of Biomedical Sciences, Houston, TX, USA
| | - Damiano Coraci
- Louis and Beatrice Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, NY, USA
| | - Rafał Krzysztoń
- Department of Biomedical Engineering, Stony Brook University, Stony Brook, NY, USA
- Louis and Beatrice Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, NY, USA
| | - Joshua Azukas
- Department of Biomedical Engineering, Stony Brook University, Stony Brook, NY, USA
| | - Nicholas Van Nest
- Department of Biomedical Engineering, Stony Brook University, Stony Brook, NY, USA
- Louis and Beatrice Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, NY, USA
| | - Alex Smashnov
- Department of Biomedical Engineering, Stony Brook University, Stony Brook, NY, USA
- Louis and Beatrice Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, NY, USA
| | - Yi-Jye Chern
- Department of Pharmacological Sciences, Stony Brook University, Stony Brook, NY, USA
- Department of Molecular and Cellular Biology, Roswell Park Cancer Institute, Buffalo, NY, USA
| | - Daniela De Martino
- Department of Medicine, Division of Hematology and Oncology, Tisch Cancer Institute, Icahn School of Medicine at Mount Sinai, New York, NY, USA
| | - Long Chi Nguyen
- Ben May Department for Cancer Research, University of Chicago, Chicago, IL, USA
| | - Harold Bien
- Louis and Beatrice Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, NY, USA
| | - Jose Javier Bravo-Cordero
- Department of Medicine, Division of Hematology and Oncology, Tisch Cancer Institute, Icahn School of Medicine at Mount Sinai, New York, NY, USA
| | - Chia-Hsin Chan
- Department of Pharmacological Sciences, Stony Brook University, Stony Brook, NY, USA
- Department of Molecular and Cellular Biology, Roswell Park Cancer Institute, Buffalo, NY, USA
| | - Marsha Rich Rosner
- Ben May Department for Cancer Research, University of Chicago, Chicago, IL, USA
| | - Gábor Balázsi
- Department of Biomedical Engineering, Stony Brook University, Stony Brook, NY, USA.
- Louis and Beatrice Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, NY, USA.
- Stony Brook Cancer Center, Stony Brook University, Stony Brook, NY, USA.
| |
Collapse
|
16
|
Hu M, Zhang Y, Yuan Y, Ma W, Zheng Y, Gu Q, Xie XS. Correlated Protein Modules Revealing Functional Coordination of Interacting Proteins Are Detected by Single-Cell Proteomics. J Phys Chem B 2023. [PMID: 37368753 DOI: 10.1021/acs.jpcb.3c00014] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/29/2023]
Abstract
Single-cell proteomics has attracted a lot of attention in recent years because it offers more functional relevance than single-cell transcriptomics. However, most work to date has focused on cell typing, which has been widely accomplished by single-cell transcriptomics. Here we report the use of single-cell proteomics to measure the correlation between the translational levels of a pair of proteins in a single mammalian cell. In measuring pairwise correlations among ∼1000 proteins in a population of homogeneous K562 cells under a steady-state condition, we observed multiple correlated protein modules (CPMs), each containing a group of highly positively correlated proteins that are functionally interacting and collectively involved in certain biological functions, such as protein synthesis and oxidative phosphorylation. Some CPMs are shared across different cell types while others are cell-type specific. Widely studied in omics analyses, pairwise correlations are often measured by introducing perturbations into bulk samples. However, some correlations of gene or protein expression under the steady-state condition would be masked by perturbation. The single-cell correlations probed in our experiment reflect intrinsic steady-state fluctuations in the absence of perturbation. We note that observed correlations between proteins are experimentally more distinct and functionally more relevant than those between corresponding mRNAs measured in single-cell transcriptomics. By virtue of single-cell proteomics, functional coordination of proteins is manifested through CPMs.
Collapse
Affiliation(s)
- Mo Hu
- Beijing Advanced Innovation Center for Genomics, Peking University, Beijing 100871, China
- Biomedical Pioneering Innovation Center (BIOPIC), Peking University, Beijing 100871, China
- Changping Laboratory, Beijing 102206, China
| | - Yutong Zhang
- Beijing Advanced Innovation Center for Genomics, Peking University, Beijing 100871, China
- Biomedical Pioneering Innovation Center (BIOPIC), Peking University, Beijing 100871, China
- College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| | - Yuan Yuan
- Beijing Advanced Innovation Center for Genomics, Peking University, Beijing 100871, China
- Biomedical Pioneering Innovation Center (BIOPIC), Peking University, Beijing 100871, China
- College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| | - Wenping Ma
- Beijing Advanced Innovation Center for Genomics, Peking University, Beijing 100871, China
- Biomedical Pioneering Innovation Center (BIOPIC), Peking University, Beijing 100871, China
- Peking-Tsinghua Center for Life Sciences (CLS), Academy for Advanced Interdisciplinary Studies, Peking University, Beijing 100871, China
| | - Yinghui Zheng
- Beijing Advanced Innovation Center for Genomics, Peking University, Beijing 100871, China
- Biomedical Pioneering Innovation Center (BIOPIC), Peking University, Beijing 100871, China
| | | | - X Sunney Xie
- Beijing Advanced Innovation Center for Genomics, Peking University, Beijing 100871, China
- Biomedical Pioneering Innovation Center (BIOPIC), Peking University, Beijing 100871, China
- Changping Laboratory, Beijing 102206, China
| |
Collapse
|
17
|
Zhang Y, Sun H, Lian X, Tang J, Zhu F. ANPELA: Significantly Enhanced Quantification Tool for Cytometry-Based Single-Cell Proteomics. ADVANCED SCIENCE (WEINHEIM, BADEN-WURTTEMBERG, GERMANY) 2023; 10:e2207061. [PMID: 36950745 DOI: 10.1002/advs.202207061] [Citation(s) in RCA: 8] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/30/2022] [Revised: 02/13/2023] [Indexed: 05/27/2023]
Abstract
ANPELA is widely used for quantifying traditional bulk proteomic data. Recently, there is a clear shift from bulk proteomics to the single-cell ones (SCP), for which powerful cytometry techniques demonstrate the fantastic capacity of capturing cellular heterogeneity that is completely overlooked by traditional bulk profiling. However, the in-depth and high-quality quantification of SCP data is still challenging and severely affected by the large numbers of quantification workflows and extreme performance dependence on the studied datasets. In other words, the proper selection of well-performing workflow(s) for any studied dataset is elusory, and it is urgently needed to have a significantly enhanced and accelerated tool to address this issue. However, no such tool is developed yet. Herein, ANPELA is therefore updated to its 2.0 version (https://idrblab.org/anpela/), which is unique in providing the most comprehensive set of quantification alternatives (>1000 workflows) among all existing tools, enabling systematic performance evaluation from multiple perspectives based on machine learning, and identifying the optimal workflow(s) using overall performance ranking together with the parallel computation. Extensive validation on different benchmark datasets and representative application scenarios suggest the great application potential of ANPELA in current SCP research for gaining more accurate and reliable biological insights.
Collapse
Affiliation(s)
- Ying Zhang
- College of Pharmaceutical Sciences, The Second Affiliated Hospital, Zhejiang University School of Medicine, Zhejiang University, Hangzhou, 310058, China
| | - Huaicheng Sun
- College of Pharmaceutical Sciences, The Second Affiliated Hospital, Zhejiang University School of Medicine, Zhejiang University, Hangzhou, 310058, China
| | - Xichen Lian
- College of Pharmaceutical Sciences, The Second Affiliated Hospital, Zhejiang University School of Medicine, Zhejiang University, Hangzhou, 310058, China
| | - Jing Tang
- Department of Bioinformatics, Chongqing Medical University, Chongqing, 400016, China
| | - Feng Zhu
- College of Pharmaceutical Sciences, The Second Affiliated Hospital, Zhejiang University School of Medicine, Zhejiang University, Hangzhou, 310058, China
- Innovation Institute for Artificial Intelligence in Medicine of Zhejiang University, Alibaba-Zhejiang University Joint Research Center of Future Digital Healthcare, Hangzhou, 330110, China
| |
Collapse
|
18
|
Liu C, Kudo T, Ye X, Gascoigne K. Cell-to-cell variability in Myc dynamics drives transcriptional heterogeneity in cancer cells. Cell Rep 2023; 42:112401. [PMID: 37060565 DOI: 10.1016/j.celrep.2023.112401] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/14/2022] [Revised: 03/07/2023] [Accepted: 03/31/2023] [Indexed: 04/16/2023] Open
Abstract
Cell-to-cell heterogeneity is vital for tumor evolution and survival. How cancer cells achieve and exploit this heterogeneity remains an active area of research. Here, we identify c-Myc as a highly heterogeneously expressed transcription factor and an orchestrator of transcriptional and phenotypic diversity in cancer cells. By monitoring endogenous c-Myc protein in individual living cells, we report the surprising pulsatile nature of c-Myc expression and the extensive cell-to-cell variability in its dynamics. We further show that heterogeneity in c-Myc dynamics leads to variable target gene transcription and that timing of c-Myc expression predicts cell-cycle progression rates and drug sensitivities. Together, our data advocate for a model in which cancer cells increase the heterogeneity of functionally diverse transcription factors such as c-Myc to rapidly survey transcriptional landscapes and survive stress.
Collapse
Affiliation(s)
- Chad Liu
- Department of Discovery Oncology, Genentech, Inc., South San Francisco, CA 94080, USA
| | - Takamasa Kudo
- Department of Cellular and Tissue Genomics, Genentech, Inc., South San Francisco, CA 94080, USA
| | - Xin Ye
- Department of Discovery Oncology, Genentech, Inc., South San Francisco, CA 94080, USA
| | - Karen Gascoigne
- Department of Discovery Oncology, Genentech, Inc., South San Francisco, CA 94080, USA.
| |
Collapse
|
19
|
Barua A, Hatzikirou H. Cell Decision Making through the Lens of Bayesian Learning. ENTROPY (BASEL, SWITZERLAND) 2023; 25:e25040609. [PMID: 37190396 PMCID: PMC10137733 DOI: 10.3390/e25040609] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/06/2023] [Revised: 03/20/2023] [Accepted: 03/29/2023] [Indexed: 05/17/2023]
Abstract
Cell decision making refers to the process by which cells gather information from their local microenvironment and regulate their internal states to create appropriate responses. Microenvironmental cell sensing plays a key role in this process. Our hypothesis is that cell decision-making regulation is dictated by Bayesian learning. In this article, we explore the implications of this hypothesis for internal state temporal evolution. By using a timescale separation between internal and external variables on the mesoscopic scale, we derive a hierarchical Fokker-Planck equation for cell-microenvironment dynamics. By combining this with the Bayesian learning hypothesis, we find that changes in microenvironmental entropy dominate the cell state probability distribution. Finally, we use these ideas to understand how cell sensing impacts cell decision making. Notably, our formalism allows us to understand cell state dynamics even without exact biochemical information about cell sensing processes by considering a few key parameters.
Collapse
Affiliation(s)
- Arnab Barua
- Departement de Biochimie, Université de Montréal, Montréal, QC H3T 1C5, Canada
- Centre Robert-Cedergren en Bio-Informatique et Génomique, Université de Montréal, Montréal, QC H3C 3J7, Canada
| | - Haralampos Hatzikirou
- Center for Information Services and High Performance Computing, Technische Univesität Dresden, 01062 Dresden, Germany
- Mathematics Department, Khalifa University, Abu Dhabi P.O. Box 127788, United Arab Emirates
| |
Collapse
|
20
|
Hastings JF, Latham SL, Kamili A, Wheatley MS, Han JZ, Wong-Erasmus M, Phimmachanh M, Nobis M, Pantarelli C, Cadell AL, O’Donnell YE, Leong KH, Lynn S, Geng FS, Cui L, Yan S, Achinger-Kawecka J, Stirzaker C, Norris MD, Haber M, Trahair TN, Speleman F, De Preter K, Cowley MJ, Bogdanovic O, Timpson P, Cox TR, Kolch W, Fletcher JI, Fey D, Croucher DR. Memory of stochastic single-cell apoptotic signaling promotes chemoresistance in neuroblastoma. SCIENCE ADVANCES 2023; 9:eabp8314. [PMID: 36867694 PMCID: PMC9984174 DOI: 10.1126/sciadv.abp8314] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 02/28/2022] [Accepted: 01/27/2023] [Indexed: 06/18/2023]
Abstract
Gene expression noise is known to promote stochastic drug resistance through the elevated expression of individual genes in rare cancer cells. However, we now demonstrate that chemoresistant neuroblastoma cells emerge at a much higher frequency when the influence of noise is integrated across multiple components of an apoptotic signaling network. Using a JNK activity biosensor with longitudinal high-content and in vivo intravital imaging, we identify a population of stochastic, JNK-impaired, chemoresistant cells that exist because of noise within this signaling network. Furthermore, we reveal that the memory of this initially random state is retained following chemotherapy treatment across a series of in vitro, in vivo, and patient models. Using matched PDX models established at diagnosis and relapse from individual patients, we show that HDAC inhibitor priming cannot erase the memory of this resistant state within relapsed neuroblastomas but improves response in the first-line setting by restoring drug-induced JNK activity within the chemoresistant population of treatment-naïve tumors.
Collapse
Affiliation(s)
- Jordan F. Hastings
- Cancer Ecosystems Program, Garvan Institute of Medical Research, Sydney, NSW 2010, Australia
| | - Sharissa L. Latham
- Cancer Ecosystems Program, Garvan Institute of Medical Research, Sydney, NSW 2010, Australia
- School of Clinical Medicine, Faculty of Medicine and Health, UNSW Sydney, Sydney, NSW, Australia
| | - Alvin Kamili
- School of Clinical Medicine, Faculty of Medicine and Health, UNSW Sydney, Sydney, NSW, Australia
- Children’s Cancer Institute, Lowy Cancer Research Centre, UNSW Sydney, Sydney, NSW, Australia
| | - Madeleine S. Wheatley
- Children’s Cancer Institute, Lowy Cancer Research Centre, UNSW Sydney, Sydney, NSW, Australia
| | - Jeremy Z. R. Han
- Cancer Ecosystems Program, Garvan Institute of Medical Research, Sydney, NSW 2010, Australia
| | - Marie Wong-Erasmus
- Children’s Cancer Institute, Lowy Cancer Research Centre, UNSW Sydney, Sydney, NSW, Australia
| | - Monica Phimmachanh
- Cancer Ecosystems Program, Garvan Institute of Medical Research, Sydney, NSW 2010, Australia
| | - Max Nobis
- Cancer Ecosystems Program, Garvan Institute of Medical Research, Sydney, NSW 2010, Australia
- School of Clinical Medicine, Faculty of Medicine and Health, UNSW Sydney, Sydney, NSW, Australia
| | - Chiara Pantarelli
- Cancer Ecosystems Program, Garvan Institute of Medical Research, Sydney, NSW 2010, Australia
| | - Antonia L. Cadell
- Cancer Ecosystems Program, Garvan Institute of Medical Research, Sydney, NSW 2010, Australia
| | - Yolande E. I. O’Donnell
- Cancer Ecosystems Program, Garvan Institute of Medical Research, Sydney, NSW 2010, Australia
| | - King Ho Leong
- Cancer Ecosystems Program, Garvan Institute of Medical Research, Sydney, NSW 2010, Australia
| | - Sophie Lynn
- Cancer Ecosystems Program, Garvan Institute of Medical Research, Sydney, NSW 2010, Australia
| | - Fan-Suo Geng
- Cancer Ecosystems Program, Garvan Institute of Medical Research, Sydney, NSW 2010, Australia
| | - Lujing Cui
- Children’s Cancer Institute, Lowy Cancer Research Centre, UNSW Sydney, Sydney, NSW, Australia
| | - Sabrina Yan
- Children’s Cancer Institute, Lowy Cancer Research Centre, UNSW Sydney, Sydney, NSW, Australia
| | - Joanna Achinger-Kawecka
- Cancer Ecosystems Program, Garvan Institute of Medical Research, Sydney, NSW 2010, Australia
- School of Clinical Medicine, Faculty of Medicine and Health, UNSW Sydney, Sydney, NSW, Australia
| | - Clare Stirzaker
- Cancer Ecosystems Program, Garvan Institute of Medical Research, Sydney, NSW 2010, Australia
- School of Clinical Medicine, Faculty of Medicine and Health, UNSW Sydney, Sydney, NSW, Australia
| | - Murray D. Norris
- Children’s Cancer Institute, Lowy Cancer Research Centre, UNSW Sydney, Sydney, NSW, Australia
- University of New South Wales Centre for Childhood Cancer Research, UNSW Sydney, Sydney, NSW, Australia
| | - Michelle Haber
- School of Clinical Medicine, Faculty of Medicine and Health, UNSW Sydney, Sydney, NSW, Australia
- Children’s Cancer Institute, Lowy Cancer Research Centre, UNSW Sydney, Sydney, NSW, Australia
| | - Toby N. Trahair
- School of Clinical Medicine, Faculty of Medicine and Health, UNSW Sydney, Sydney, NSW, Australia
- Children’s Cancer Institute, Lowy Cancer Research Centre, UNSW Sydney, Sydney, NSW, Australia
- Kids Cancer Centre, Sydney Children’s Hospital, Randwick, NSW 2031, Australia
| | - Frank Speleman
- Center for Medical Genetics, Ghent University, Ghent, Belgium
- Cancer Research Institute Ghent, Ghent University, Ghent, Belgium
| | - Katleen De Preter
- Center for Medical Genetics, Ghent University, Ghent, Belgium
- Cancer Research Institute Ghent, Ghent University, Ghent, Belgium
| | - Mark J. Cowley
- School of Clinical Medicine, Faculty of Medicine and Health, UNSW Sydney, Sydney, NSW, Australia
- Children’s Cancer Institute, Lowy Cancer Research Centre, UNSW Sydney, Sydney, NSW, Australia
- University of New South Wales Centre for Childhood Cancer Research, UNSW Sydney, Sydney, NSW, Australia
| | - Ozren Bogdanovic
- Cancer Ecosystems Program, Garvan Institute of Medical Research, Sydney, NSW 2010, Australia
- School of Clinical Medicine, Faculty of Medicine and Health, UNSW Sydney, Sydney, NSW, Australia
| | - Paul Timpson
- Cancer Ecosystems Program, Garvan Institute of Medical Research, Sydney, NSW 2010, Australia
- School of Clinical Medicine, Faculty of Medicine and Health, UNSW Sydney, Sydney, NSW, Australia
| | - Thomas R. Cox
- Cancer Ecosystems Program, Garvan Institute of Medical Research, Sydney, NSW 2010, Australia
- School of Clinical Medicine, Faculty of Medicine and Health, UNSW Sydney, Sydney, NSW, Australia
| | - Walter Kolch
- Systems Biology Ireland, School of Medicine, University College Dublin, Belfield, Dublin 4, Ireland
- Conway Institute of Biomolecular and Biomedical Research, University College Dublin, Belfield, Dublin 4, Ireland
| | - Jamie I. Fletcher
- School of Clinical Medicine, Faculty of Medicine and Health, UNSW Sydney, Sydney, NSW, Australia
- Children’s Cancer Institute, Lowy Cancer Research Centre, UNSW Sydney, Sydney, NSW, Australia
- University of New South Wales Centre for Childhood Cancer Research, UNSW Sydney, Sydney, NSW, Australia
| | - Dirk Fey
- Systems Biology Ireland, School of Medicine, University College Dublin, Belfield, Dublin 4, Ireland
| | - David R. Croucher
- Cancer Ecosystems Program, Garvan Institute of Medical Research, Sydney, NSW 2010, Australia
- School of Clinical Medicine, Faculty of Medicine and Health, UNSW Sydney, Sydney, NSW, Australia
| |
Collapse
|
21
|
Shlyakhtina Y, Bloechl B, Portal MM. BdLT-Seq as a barcode decay-based method to unravel lineage-linked transcriptome plasticity. Nat Commun 2023; 14:1085. [PMID: 36841849 PMCID: PMC9968323 DOI: 10.1038/s41467-023-36744-1] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/26/2022] [Accepted: 02/14/2023] [Indexed: 02/26/2023] Open
Abstract
Cell plasticity is a core biological process underlying a myriad of molecular and cellular events taking place throughout organismal development and evolution. It has been postulated that cellular systems thrive to balance the organization of meta-stable states underlying this phenomenon, thereby maintaining a degree of populational homeostasis compatible with an ever-changing environment and, thus, life. Notably, albeit circumstantial evidence has been gathered in favour of the latter conceptual framework, a direct observation of meta-state dynamics and the biological consequences of such a process in generating non-genetic clonal diversity and divergent phenotypic output remains largely unexplored. To fill this void, here we develop a lineage-tracing technology termed Barcode decay Lineage Tracing-Seq. BdLT-Seq is based on episome-encoded molecular identifiers that, supported by the dynamic decay of the tracing information upon cell division, ascribe directionality to a cell lineage tree whilst directly coupling non-genetic molecular features to phenotypes in comparable genomic landscapes. We show that cell transcriptome states are both inherited, and dynamically reshaped following constrained rules encoded within the cell lineage in basal growth conditions, upon oncogene activation and throughout the process of reversible resistance to therapeutic cues thus adjusting phenotypic output leading to intra-clonal non-genetic diversity.
Collapse
Affiliation(s)
- Yelyzaveta Shlyakhtina
- Cell Plasticity & Epigenetics Lab, Cancer Research UK - Manchester Institute, The University of Manchester, SK10 4TG, Manchester, UK
| | - Bianca Bloechl
- Cell Plasticity & Epigenetics Lab, Cancer Research UK - Manchester Institute, The University of Manchester, SK10 4TG, Manchester, UK
| | - Maximiliano M Portal
- Cell Plasticity & Epigenetics Lab, Cancer Research UK - Manchester Institute, The University of Manchester, SK10 4TG, Manchester, UK.
| |
Collapse
|
22
|
Dichamp J, Cellière G, Ghallab A, Hassan R, Boissier N, Hofmann U, Reinders J, Sezgin S, Zühlke S, Hengstler JG, Drasdo D. In vitro to in vivo acetaminophen hepatotoxicity extrapolation using classical schemes, pharmacodynamic models and a multiscale spatial-temporal liver twin. Front Bioeng Biotechnol 2023; 11:1049564. [PMID: 36815881 PMCID: PMC9932319 DOI: 10.3389/fbioe.2023.1049564] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/20/2022] [Accepted: 01/10/2023] [Indexed: 02/05/2023] Open
Abstract
In vitro to in vivo extrapolation represents a critical challenge in toxicology. In this paper we explore extrapolation strategies for acetaminophen (APAP) based on mechanistic models, comparing classical (CL) homogeneous compartment pharmacodynamic (PD) models and a spatial-temporal (ST), multiscale digital twin model resolving liver microarchitecture at cellular resolution. The models integrate consensus detoxification reactions in each individual hepatocyte. We study the consequences of the two model types on the extrapolation and show in which cases these models perform better than the classical extrapolation strategy that is based either on the maximal drug concentration (Cmax) or the area under the pharmacokinetic curve (AUC) of the drug blood concentration. We find that an CL-model based on a well-mixed blood compartment is sufficient to correctly predict the in vivo toxicity from in vitro data. However, the ST-model that integrates more experimental information requires a change of at least one parameter to obtain the same prediction, indicating that spatial compartmentalization may indeed be an important factor.
Collapse
Affiliation(s)
- Jules Dichamp
- Group SIMBIOTX, INRIA Saclay-Île-de-France, Palaiseau, France,Leibniz Research Centre for Working Environment and Human Factors, Technical University Dortmund, Dortmund, Germany,Group MAMBA, INRIA Paris, Paris, France
| | | | - Ahmed Ghallab
- Leibniz Research Centre for Working Environment and Human Factors, Technical University Dortmund, Dortmund, Germany,Department of Forensic Medicine and Toxicology, Faculty of Veterinary Medicine, South Valley University, Qena, Egypt
| | - Reham Hassan
- Leibniz Research Centre for Working Environment and Human Factors, Technical University Dortmund, Dortmund, Germany,Department of Forensic Medicine and Toxicology, Faculty of Veterinary Medicine, South Valley University, Qena, Egypt
| | - Noemie Boissier
- Group SIMBIOTX, INRIA Saclay-Île-de-France, Palaiseau, France
| | - Ute Hofmann
- Dr. Margarete Fischer-Bosch Institute of Clinical Pharmacology and University of Tübingen, Stuttgart, Germany
| | - Joerg Reinders
- Leibniz Research Centre for Working Environment and Human Factors, Technical University Dortmund, Dortmund, Germany
| | - Selahaddin Sezgin
- Faculty of Chemistry and Chemical Biology, TU Dortmund, Dortmund, Germany
| | - Sebastian Zühlke
- Center for Mass Spectrometry (CMS), Faculty of Chemistry and Chemical Biology, TU Dortmund University, Dortmund, Germany
| | - Jan G. Hengstler
- Leibniz Research Centre for Working Environment and Human Factors, Technical University Dortmund, Dortmund, Germany
| | - Dirk Drasdo
- Group SIMBIOTX, INRIA Saclay-Île-de-France, Palaiseau, France,Leibniz Research Centre for Working Environment and Human Factors, Technical University Dortmund, Dortmund, Germany,Group MAMBA, INRIA Paris, Paris, France,*Correspondence: Dirk Drasdo,
| |
Collapse
|
23
|
Jain P, Corbo S, Mohammad K, Sahoo S, Ranganathan S, George JT, Levine H, Taube J, Toneff M, Jolly MK. Epigenetic memory acquired during long-term EMT induction governs the recovery to the epithelial state. J R Soc Interface 2023; 20:20220627. [PMID: 36628532 PMCID: PMC9832289 DOI: 10.1098/rsif.2022.0627] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/26/2022] [Accepted: 12/16/2022] [Indexed: 01/12/2023] Open
Abstract
Epithelial-mesenchymal transition (EMT) and its reverse mesenchymal-epithelial transition (MET) are critical during embryonic development, wound healing and cancer metastasis. While phenotypic changes during short-term EMT induction are reversible, long-term EMT induction has been often associated with irreversibility. Here, we show that phenotypic changes seen in MCF10A cells upon long-term EMT induction by TGFβ need not be irreversible, but have relatively longer time scales of reversibility than those seen in short-term induction. Next, using a phenomenological mathematical model to account for the chromatin-mediated epigenetic silencing of the miR-200 family by ZEB family, we highlight how the epigenetic memory gained during long-term EMT induction can slow the recovery to the epithelial state post-TGFβ withdrawal. Our results suggest that epigenetic modifiers can govern the extent and time scale of EMT reversibility and advise caution against labelling phenotypic changes seen in long-term EMT induction as 'irreversible'.
Collapse
Affiliation(s)
- Paras Jain
- Centre for BioSystems Science and Engineering, Indian Institute of Science, Bengaluru 560012, India
| | - Sophia Corbo
- Department of Biology, Widener University, Chester, PA 19013, USA
| | - Kulsoom Mohammad
- Department of Biology, Widener University, Chester, PA 19013, USA
| | - Sarthak Sahoo
- Centre for BioSystems Science and Engineering, Indian Institute of Science, Bengaluru 560012, India
| | | | - Jason T. George
- Department of Biomedical Engineering, Texas A&M University, College Station, TX 76798, USA
| | - Herbert Levine
- Center for Theoretical Biological Physics and Departments of Physics and Bioengineering, Northeastern University, Boston, MA 02115, USA
| | - Joseph Taube
- Department of Biology, Baylor University, Waco, TX 76706, USA
| | - Michael Toneff
- Department of Biology, Widener University, Chester, PA 19013, USA
| | - Mohit Kumar Jolly
- Centre for BioSystems Science and Engineering, Indian Institute of Science, Bengaluru 560012, India
| |
Collapse
|
24
|
Sarma U, Ripka L, Anyaegbunam UA, Legewie S. Modeling Cellular Signaling Variability Based on Single-Cell Data: The TGFβ-SMAD Signaling Pathway. Methods Mol Biol 2023; 2634:215-251. [PMID: 37074581 DOI: 10.1007/978-1-0716-3008-2_10] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/20/2023]
Abstract
Nongenetic heterogeneity is key to cellular decisions, as even genetically identical cells respond in very different ways to the same external stimulus, e.g., during cell differentiation or therapeutic treatment of disease. Strong heterogeneity is typically already observed at the level of signaling pathways that are the first sensors of external inputs and transmit information to the nucleus where decisions are made. Since heterogeneity arises from random fluctuations of cellular components, mathematical models are required to fully describe the phenomenon and to understand the dynamics of heterogeneous cell populations. Here, we review the experimental and theoretical literature on cellular signaling heterogeneity, with special focus on the TGFβ/SMAD signaling pathway.
Collapse
Affiliation(s)
- Uddipan Sarma
- Institute of Molecular Biology (IMB), Mainz, Germany
| | - Lorenz Ripka
- Institute of Molecular Biology (IMB), Mainz, Germany
- Department of Systems Biology, Institute for Biomedical Genetics, University of Stuttgart, Stuttgart, Germany
| | - Uchenna Alex Anyaegbunam
- Institute of Molecular Biology (IMB), Mainz, Germany
- Department of Systems Biology, Institute for Biomedical Genetics, University of Stuttgart, Stuttgart, Germany
| | - Stefan Legewie
- Institute of Molecular Biology (IMB), Mainz, Germany.
- Department of Systems Biology, Institute for Biomedical Genetics, University of Stuttgart, Stuttgart, Germany.
- Stuttgart Research Center for Systems Biology, University of Stuttgart, Stuttgart, Germany.
| |
Collapse
|
25
|
Durrieu L, Bush A, Grande A, Johansson R, Janzén D, Katz A, Cedersund G, Colman-Lerner A. Characterization of cell-to-cell variation in nuclear transport rates and identification of its sources. iScience 2022; 26:105906. [PMID: 36686393 PMCID: PMC9852351 DOI: 10.1016/j.isci.2022.105906] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/21/2022] [Revised: 11/10/2022] [Accepted: 12/25/2022] [Indexed: 12/30/2022] Open
Abstract
Nuclear transport is an essential part of eukaryotic cell function. Here, we present scFRAP, a model-assisted fluorescent recovery after photobleaching (FRAP)- based method to determine nuclear import and export rates independently in individual live cells. To overcome the inherent noise of single-cell measurements, we performed sequential FRAPs on the same cell. We found large cell-to-cell variation in transport rates within isogenic yeast populations. For passive transport, the variability in NPC number might explain most of the variability. Using this approach, we studied mother-daughter cell asymmetry in the active nuclear shuttling of the transcription factor Ace2, which is specifically concentrated in daughter cell nuclei in early G1. Rather than reduced export in the daughter cell, as previously hypothesized, we found that this asymmetry is mainly due to an increased import in daughters. These results shed light on cell-to-cell variation in cellular dynamics and its sources.
Collapse
Affiliation(s)
- Lucía Durrieu
- Department of Physiology, Molecular and Cellular Biology, School of Exact and Natural Sciences, University of Buenos Aires (UBA), C1428EGA, Argentina,Institute of Physiology, Molecular Biology and Neurosciences, National Council of Scientific and Technical Research (IFIBYNE-UBA-CONICET), C1428EGA, Argentina
| | - Alan Bush
- Department of Physiology, Molecular and Cellular Biology, School of Exact and Natural Sciences, University of Buenos Aires (UBA), C1428EGA, Argentina,Institute of Physiology, Molecular Biology and Neurosciences, National Council of Scientific and Technical Research (IFIBYNE-UBA-CONICET), C1428EGA, Argentina,Department of Biomedical Engineering, Linköping University, Linköping, Sweden
| | - Alicia Grande
- Department of Physiology, Molecular and Cellular Biology, School of Exact and Natural Sciences, University of Buenos Aires (UBA), C1428EGA, Argentina,Institute of Physiology, Molecular Biology and Neurosciences, National Council of Scientific and Technical Research (IFIBYNE-UBA-CONICET), C1428EGA, Argentina
| | - Rikard Johansson
- Department of Biomedical Engineering, Linköping University, Linköping, Sweden
| | - David Janzén
- Department of Biomedical Engineering, Linköping University, Linköping, Sweden
| | - Andrea Katz
- Department of Physiology, Molecular and Cellular Biology, School of Exact and Natural Sciences, University of Buenos Aires (UBA), C1428EGA, Argentina
| | - Gunnar Cedersund
- Department of Biomedical Engineering, Linköping University, Linköping, Sweden
| | - Alejandro Colman-Lerner
- Department of Physiology, Molecular and Cellular Biology, School of Exact and Natural Sciences, University of Buenos Aires (UBA), C1428EGA, Argentina,Institute of Physiology, Molecular Biology and Neurosciences, National Council of Scientific and Technical Research (IFIBYNE-UBA-CONICET), C1428EGA, Argentina,Corresponding author
| |
Collapse
|
26
|
Lee U, Mortola EN, Kim EJ, Long M. Evolution and maintenance of phenotypic plasticity. Biosystems 2022; 222:104791. [DOI: 10.1016/j.biosystems.2022.104791] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/07/2022] [Revised: 09/20/2022] [Accepted: 10/03/2022] [Indexed: 11/02/2022]
|
27
|
Govindaraj V, Sarma S, Karulkar A, Purwar R, Kar S. Transcriptional Fluctuations Govern the Serum-Dependent Cell Cycle Duration Heterogeneities in Mammalian Cells. ACS Synth Biol 2022; 11:3743-3758. [PMID: 36325971 DOI: 10.1021/acssynbio.2c00347] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Abstract
Mammalian cells exhibit a high degree of intercellular variability in cell cycle period and phase durations. However, the factors orchestrating the cell cycle duration heterogeneities remain unclear. Herein, by combining cell cycle network-based mathematical models with live single-cell imaging studies under varied serum conditions, we demonstrate that fluctuating transcription rates of cell cycle regulatory genes across cell lineages and during cell cycle progression in mammalian cells majorly govern the robust correlation patterns of cell cycle period and phase durations among sister, cousin, and mother-daughter lineage pairs. However, for the overall cellular population, alteration in the serum level modulates the fluctuation and correlation patterns of cell cycle period and phase durations in a correlated manner. These heterogeneities at the population level can be fine-tuned under limited serum conditions by perturbing the cell cycle network using a p38-signaling inhibitor without affecting the robust lineage-level correlations. Overall, our approach identifies transcriptional fluctuations as the key controlling factor for the cell cycle duration heterogeneities and predicts ways to reduce cell-to-cell variabilities by perturbing the cell cycle network regulations.
Collapse
Affiliation(s)
| | - Subrot Sarma
- Department of Chemistry, IIT Bombay, Powai, Mumbai 400076, India
| | - Atharva Karulkar
- Department of Biosciences and Bioengineering, IIT Bombay, Powai, Mumbai 400076, India
| | - Rahul Purwar
- Department of Biosciences and Bioengineering, IIT Bombay, Powai, Mumbai 400076, India
| | - Sandip Kar
- Department of Chemistry, IIT Bombay, Powai, Mumbai 400076, India
| |
Collapse
|
28
|
Mohammadi F, Visagan S, Gross SM, Karginov L, Lagarde JC, Heiser LM, Meyer AS. A lineage tree-based hidden Markov model quantifies cellular heterogeneity and plasticity. Commun Biol 2022; 5:1258. [PMID: 36396800 PMCID: PMC9671968 DOI: 10.1038/s42003-022-04208-9] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/09/2021] [Accepted: 11/01/2022] [Indexed: 11/18/2022] Open
Abstract
Individual cells can assume a variety of molecular and phenotypic states and recent studies indicate that cells can rapidly adapt in response to therapeutic stress. Such phenotypic plasticity may confer resistance, but also presents opportunities to identify molecular programs that could be targeted for therapeutic benefit. Approaches to quantify tumor-drug responses typically focus on snapshot, population-level measurements. While informative, these methods lack lineage and temporal information, which are particularly critical for understanding dynamic processes such as cell state switching. As new technologies have become available to measure lineage relationships, modeling approaches will be needed to identify the forms of cell-to-cell heterogeneity present in these data. Here we apply a lineage tree-based adaptation of a hidden Markov model that employs single cell lineages as input to learn the characteristic patterns of phenotypic heterogeneity and state transitions. In benchmarking studies, we demonstrated that the model successfully classifies cells within experimentally-tractable dataset sizes. As an application, we analyzed experimental measurements in cancer and non-cancer cell populations under various treatments. We find evidence of multiple phenotypically distinct states, with considerable heterogeneity and unique drug responses. In total, this framework allows for the flexible modeling of single cell heterogeneity across lineages to quantify, understand, and control cell state switching.
Collapse
Affiliation(s)
- Farnaz Mohammadi
- Department of Bioengineering, University of California, Los Angeles, CA, USA
| | - Shakthi Visagan
- Department of Bioengineering, University of California, Los Angeles, CA, USA
| | - Sean M Gross
- Department of Biomedical Engineering, Oregon Health & Science University, Portland, OR, USA
| | - Luka Karginov
- Department of Bioengineering, University of Illinois, Urbana Champaign, IL, USA
| | - J C Lagarde
- Department of Bioengineering, University of California, Los Angeles, CA, USA
| | - Laura M Heiser
- Department of Biomedical Engineering, Oregon Health & Science University, Portland, OR, USA
| | - Aaron S Meyer
- Department of Bioengineering, University of California, Los Angeles, CA, USA.
- Department of Bioinformatics, University of California, Los Angeles, CA, USA.
- Jonsson Comprehensive Cancer Center, University of California, Los Angeles, CA, USA.
- Eli and Edythe Broad Center of Regenerative Medicine and Stem Cell Research, University of California, Los Angeles, CA, USA.
| |
Collapse
|
29
|
Maier BD, Aguilera LU, Sahle S, Mutz P, Kalra P, Dächert C, Bartenschlager R, Binder M, Kummer U. Stochastic dynamics of Type-I interferon responses. PLoS Comput Biol 2022; 18:e1010623. [PMID: 36269758 PMCID: PMC9629604 DOI: 10.1371/journal.pcbi.1010623] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/13/2022] [Revised: 11/02/2022] [Accepted: 09/30/2022] [Indexed: 11/07/2022] Open
Abstract
Interferon (IFN) activates the transcription of several hundred of IFN stimulated genes (ISGs) that constitute a highly effective antiviral defense program. Cell-to-cell variability in the induction of ISGs is well documented, but its source and effects are not completely understood. The molecular mechanisms behind this heterogeneity have been related to randomness in molecular events taking place during the JAK-STAT signaling pathway. Here, we study the sources of variability in the induction of the IFN-alpha response by using MxA and IFIT1 activation as read-out. To this end, we integrate time-resolved flow cytometry data and stochastic modeling of the JAK-STAT signaling pathway. The complexity of the IFN response was matched by fitting probability distributions to time-course flow cytometry snapshots. Both, experimental data and simulations confirmed that the MxA and IFIT1 induction circuits generate graded responses rather than all-or-none responses. Subsequently, we quantify the size of the intrinsic variability at different steps in the pathway. We found that stochastic effects are transiently strong during the ligand-receptor activation steps and the formation of the ISGF3 complex, but negligible for the final induction of the studied ISGs. We conclude that the JAK-STAT signaling pathway is a robust biological circuit that efficiently transmits information under stochastic environments. We investigate the impact of intrinsic and extrinsic noise on the reliability of interferon signaling. Information must be transduced robustly despite existing biochemical variability and at the same time the system has to allow for cellular variability to tune it against changing environments. Getting insights into stochasticity in signaling networks is crucial to understand cellular dynamics and decision-making processes. To this end, we developed a detailed stochastic computational model based on single cell data. We are able to show that reliability is achieved despite high noise at the receptor level.
Collapse
Affiliation(s)
- Benjamin D. Maier
- Department of Modeling of Biological Processes, COS Heidelberg / Bioquant, Heidelberg University, Heidelberg, Germany
| | - Luis U. Aguilera
- Department of Modeling of Biological Processes, COS Heidelberg / Bioquant, Heidelberg University, Heidelberg, Germany
| | - Sven Sahle
- Department of Modeling of Biological Processes, COS Heidelberg / Bioquant, Heidelberg University, Heidelberg, Germany
| | - Pascal Mutz
- Division Virus-Associated Carcinogenesis, German Cancer Research Center (DKFZ), Heidelberg, Germany
- Department for Infectious Diseases, Molecular Virology, Medical Faculty, Heidelberg University, Heidelberg, Germany
| | - Priyata Kalra
- Department of Modeling of Biological Processes, COS Heidelberg / Bioquant, Heidelberg University, Heidelberg, Germany
| | - Christopher Dächert
- Research Group “Dynamics of early viral infection and the innate antiviral response”, German Cancer Research Center (DKFZ), Heidelberg, Germany
- Department for Infectious Diseases, Molecular Virology, Medical Faculty, Heidelberg University, Heidelberg, Germany
| | - Ralf Bartenschlager
- Division Virus-Associated Carcinogenesis, German Cancer Research Center (DKFZ), Heidelberg, Germany
- Department for Infectious Diseases, Molecular Virology, Medical Faculty, Heidelberg University, Heidelberg, Germany
| | - Marco Binder
- Research Group “Dynamics of early viral infection and the innate antiviral response”, German Cancer Research Center (DKFZ), Heidelberg, Germany
| | - Ursula Kummer
- Department of Modeling of Biological Processes, COS Heidelberg / Bioquant, Heidelberg University, Heidelberg, Germany
- * E-mail:
| |
Collapse
|
30
|
Analysis of the Human Protein Atlas Weakly Supervised Single-Cell Classification competition. Nat Methods 2022; 19:1221-1229. [PMID: 36175767 PMCID: PMC9550622 DOI: 10.1038/s41592-022-01606-z] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/14/2021] [Accepted: 08/10/2022] [Indexed: 11/09/2022]
Abstract
While spatial proteomics by fluorescence imaging has quickly become an essential discovery tool for researchers, fast and scalable methods to classify and embed single-cell protein distributions in such images are lacking. Here, we present the design and analysis of the results from the competition Human Protein Atlas - Single-Cell Classification hosted on the Kaggle platform. This represents a crowd-sourced competition to develop machine learning models trained on limited annotations to label single-cell protein patterns in fluorescent images. The particular challenges of this competition include class imbalance, weak labels and multi-label classification, prompting competitors to apply a wide range of approaches in their solutions. The winning models serve as the first subcellular omics tools that can annotate single-cell locations, extract single-cell features and capture cellular dynamics.
Collapse
|
31
|
Batra A, Banerjee SC, Sharma R. Persistent Correlation in Cellular Noise Determines Longevity of Viral Infections. J Phys Chem Lett 2022; 13:7252-7260. [PMID: 35913772 DOI: 10.1021/acs.jpclett.2c01875] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/15/2023]
Abstract
The slowly decaying viral dynamics, even after 2-3 weeks from diagnosis, is one of the characteristics of COVID-19 infection that is still unexplored in theoretical and experimental studies. This long-lived characteristic of viral infections in the framework of inherent variations or noise present at the cellular level is often overlooked. Therefore, in this work, we aim to understand the effect of these variations by proposing a stochastic non-Markovian model that not only captures the coupled dynamics between the immune cells and the virus but also enables the study of the effect of fluctuations. Numerical simulations of our model reveal that the long-range temporal correlations in fluctuations dictate the long-lived dynamics of a viral infection and, in turn, also affect the rates of immune response. Furthermore, predictions of our model system are in agreement with the experimental viral load data of COVID-19 patients from various countries.
Collapse
Affiliation(s)
- Abhilasha Batra
- Department of Chemistry, Indian Institute of Science Education and Research (IISER), Bhopal, Madhya Pradesh 462066, India
| | - Shoubhik Chandan Banerjee
- Department of Biological Sciences, Indian Institute of Science Education and Research (IISER), Bhopal, Madhya Pradesh 462066, India
| | - Rati Sharma
- Department of Chemistry, Indian Institute of Science Education and Research (IISER), Bhopal, Madhya Pradesh 462066, India
| |
Collapse
|
32
|
Que-Salinas U, Martinez-Peon D, Reyes-Figueroa AD, Ibarra I, Scheckhuber CQ. On the Prediction of In Vitro Arginine Glycation of Short Peptides Using Artificial Neural Networks. SENSORS (BASEL, SWITZERLAND) 2022; 22:5237. [PMID: 35890916 PMCID: PMC9324327 DOI: 10.3390/s22145237] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 06/22/2022] [Revised: 07/08/2022] [Accepted: 07/11/2022] [Indexed: 06/15/2023]
Abstract
One of the hallmarks of diabetes is an increased modification of cellular proteins. The most prominent type of modification stems from the reaction of methylglyoxal with arginine and lysine residues, leading to structural and functional impairments of target proteins. For lysine glycation, several algorithms allow a prediction of occurrence; thus, making it possible to pinpoint likely targets. However, according to our knowledge, no approaches have been published for predicting the likelihood of arginine glycation. There are indications that arginine and not lysine is the most prominent target for the toxic dialdehyde. One of the reasons why there is no arginine glycation predictor is the limited availability of quantitative data. Here, we used a recently published high-quality dataset of arginine modification probabilities to employ an artificial neural network strategy. Despite the limited data availability, our results achieve an accuracy of about 75% of correctly predicting the exact value of the glycation probability of an arginine-containing peptide without setting thresholds upon whether it is decided if a given arginine is modified or not. This contribution suggests a solution for predicting arginine glycation of short peptides.
Collapse
Affiliation(s)
- Ulices Que-Salinas
- Centro de Ciencias de la Tierra, Universidad Veracruzana, Xalapa 91090, VER, Mexico;
| | - Dulce Martinez-Peon
- Department of Electrical and Electronic Engineering, National Technological Institute of Mexico/IT, Monterrey 67170, NL, Mexico;
| | - Angel D. Reyes-Figueroa
- Consejo Nacional de Ciencia y Tecnología, Av. Insurgentes Sur 1582, Col. Crédito Constructor, Benito Juárez, Mexico City 03940, DF, Mexico;
- Centro de Investigación en Matemáticas Unidad Monterrey, Parque de Investigación e Innovación Tecnológica (PIIT), Av. Alianza Centro No. 502, Apodaca 66628, NL, Mexico
| | - Ivonne Ibarra
- Independent Researcher, Monterrey 66620, NL, Mexico;
| | - Christian Quintus Scheckhuber
- Departamento de Bioingeniería, Escuela de Ingeniería y Ciencias, Tecnologico de Monterrey, Ave. Eugenio Garza Sada 2501, Monterrey 64849, NL, Mexico
| |
Collapse
|
33
|
Seong J, Frias-Aldeguer J, Holzmann V, Kagawa H, Sestini G, Heidari Khoei H, Scholte Op Reimer Y, Kip M, Pradhan SJ, Verwegen L, Vivié J, Li L, Alemany A, Korving J, Darmis F, van Oudenaarden A, Ten Berge D, Geijsen N, Rivron NC. Epiblast inducers capture mouse trophectoderm stem cells in vitro and pattern blastoids for implantation in utero. Cell Stem Cell 2022; 29:1102-1118.e8. [PMID: 35803228 DOI: 10.1016/j.stem.2022.06.002] [Citation(s) in RCA: 17] [Impact Index Per Article: 8.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/12/2022] [Revised: 04/21/2022] [Accepted: 06/02/2022] [Indexed: 11/03/2022]
Abstract
The embryo instructs the allocation of cell states to spatially regulate functions. In the blastocyst, patterning of trophoblast (TR) cells ensures successful implantation and placental development. Here, we defined an optimal set of molecules secreted by the epiblast (inducers) that captures in vitro stable, highly self-renewing mouse trophectoderm stem cells (TESCs) resembling the blastocyst stage. When exposed to suboptimal inducers, these stem cells fluctuate to form interconvertible subpopulations with reduced self-renewal and facilitated differentiation, resembling peri-implantation cells, known as TR stem cells (TSCs). TESCs have enhanced capacity to form blastoids that implant more efficiently in utero due to inducers maintaining not only local TR proliferation and self-renewal, but also WNT6/7B secretion that stimulates uterine decidualization. Overall, the epiblast maintains sustained growth and decidualization potential of abutting TR cells, while, as known, distancing imposed by the blastocyst cavity differentiates TR cells for uterus adhesion, thus patterning the essential functions of implantation.
Collapse
Affiliation(s)
- Jinwoo Seong
- Institute of Molecular Biotechnology of the Austrian Academy of Sciences, Vienna Biocenter, Vienna, Austria
| | - Javier Frias-Aldeguer
- Hubrecht Institute for Developmental Biology and Stem Cell Research, Utrecht, the Netherlands; Maastricht University, Maastricht, the Netherlands
| | - Viktoria Holzmann
- Institute of Molecular Biotechnology of the Austrian Academy of Sciences, Vienna Biocenter, Vienna, Austria
| | - Harunobu Kagawa
- Institute of Molecular Biotechnology of the Austrian Academy of Sciences, Vienna Biocenter, Vienna, Austria
| | - Giovanni Sestini
- Institute of Molecular Biotechnology of the Austrian Academy of Sciences, Vienna Biocenter, Vienna, Austria
| | - Heidar Heidari Khoei
- Institute of Molecular Biotechnology of the Austrian Academy of Sciences, Vienna Biocenter, Vienna, Austria; Department of Stem Cells and Developmental Biology, Cell Science Research Center, Royan Institute for Stem Cell Biology and Technology, ACECR, Tehran, Iran
| | - Yvonne Scholte Op Reimer
- Institute of Molecular Biotechnology of the Austrian Academy of Sciences, Vienna Biocenter, Vienna, Austria
| | - Maarten Kip
- Hubrecht Institute for Developmental Biology and Stem Cell Research, Utrecht, the Netherlands
| | - Saurabh J Pradhan
- Institute of Molecular Biotechnology of the Austrian Academy of Sciences, Vienna Biocenter, Vienna, Austria
| | - Lucas Verwegen
- Department of Cell Biology, Erasmus MC, University Medical Center Rotterdam, Rotterdam, the Netherlands
| | - Judith Vivié
- Hubrecht Institute for Developmental Biology and Stem Cell Research, Utrecht, the Netherlands
| | - Linfeng Li
- Maastricht University, Maastricht, the Netherlands
| | - Anna Alemany
- Hubrecht Institute for Developmental Biology and Stem Cell Research, Utrecht, the Netherlands
| | - Jeroen Korving
- Hubrecht Institute for Developmental Biology and Stem Cell Research, Utrecht, the Netherlands
| | - Frank Darmis
- Hubrecht Institute for Developmental Biology and Stem Cell Research, Utrecht, the Netherlands
| | | | - Derk Ten Berge
- Department of Cell Biology, Erasmus MC, University Medical Center Rotterdam, Rotterdam, the Netherlands
| | - Niels Geijsen
- Hubrecht Institute for Developmental Biology and Stem Cell Research, Utrecht, the Netherlands; Department of Anatomy and Embryology, LUMC, Leiden University, Leiden, the Netherlands
| | - Nicolas C Rivron
- Institute of Molecular Biotechnology of the Austrian Academy of Sciences, Vienna Biocenter, Vienna, Austria; Hubrecht Institute for Developmental Biology and Stem Cell Research, Utrecht, the Netherlands; Maastricht University, Maastricht, the Netherlands.
| |
Collapse
|
34
|
Lannan R, Maity A, Wollman R. Epigenetic fluctuations underlie gene expression timescales and variability. Physiol Genomics 2022; 54:220-229. [PMID: 35476585 DOI: 10.1152/physiolgenomics.00051.2021] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022] Open
Abstract
Isogenic populations of mammalian cells exhibit significant gene expression variability. This variability can be separated into two sources, cis, or allele-specific sources, and trans and global processes. Furthermore, each source of variability has its own timescale. Fast timescales will result in rapid fluctuation of gene expression whereas slow timescales will result in longer persistence of gene expression levels over time. Here we investigated sources of gene expression that are intrinsic, i.e. coming from cis-regulatory factors and follow slow timescales. To do so, we developed a reporter system that isolates allele-specific variability and measures its persistence in imaging and long-term fluctuation analysis experiments. Our results identify a new source of gene expression variability that is allele-specific but that is fluctuating on timescales of days. We hypothesized that allele-specific fluctuations of epigenetic regulatory factors are responsible for the newly discovered allele-specific and slow source of gene expression variability. Using mathematical modeling we showed that the addition of this effect to the two-state model is sufficient to account for all empirical observation. Furthermore, using direct assays of chromatin markers we find fluctuation in H3K4me3 levels that match the observed changes in gene expression levels providing direct experimental support of our model. Collectively, our work shows that slow fluctuations of regulatory chromatin modifications contribute to the variability in gene expression.
Collapse
Affiliation(s)
- Ryan Lannan
- Department of Chemistry and Biochemistry, University of California, Los Angeles, California.,Department of Integrative Biology and Physiology, University of California, Los Angeles, California.,Institute of Quantitative Biosciences, University of California, Los Angeles, California
| | - Alok Maity
- Department of Chemistry and Biochemistry, University of California, Los Angeles, California.,Department of Integrative Biology and Physiology, University of California, Los Angeles, California.,Institute of Quantitative Biosciences, University of California, Los Angeles, California
| | - Roy Wollman
- Department of Chemistry and Biochemistry, University of California, Los Angeles, California.,Department of Integrative Biology and Physiology, University of California, Los Angeles, California.,Institute of Quantitative Biosciences, University of California, Los Angeles, California
| |
Collapse
|
35
|
Membrane marker selection for segmenting single cell spatial proteomics data. Nat Commun 2022; 13:1999. [PMID: 35422106 PMCID: PMC9010440 DOI: 10.1038/s41467-022-29667-w] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/04/2021] [Accepted: 03/25/2022] [Indexed: 12/21/2022] Open
Abstract
The ability to profile spatial proteomics at the single cell level enables the study of cell types, their spatial distribution, and interactions in several tissues and conditions. Current methods for cell segmentation in such studies rely on known membrane or cell boundary markers. However, for many tissues, an optimal set of markers is not known, and even within a tissue, different cell types may express different markers. Here we present RAMCES, a method that uses a convolutional neural network to learn the optimal markers for a new sample and outputs a weighted combination of the selected markers for segmentation. Testing RAMCES on several existing datasets indicates that it correctly identifies cell boundary markers, improving on methods that rely on a single marker or those that extend nuclei segmentations. Application to new spatial proteomics data demonstrates its usefulness for accurately assigning cell types based on the proteins expressed in segmented cells. Cell segmentation of single-cell spatial proteomics data remains a challenge and often relies on the selection of a membrane marker, which is not always known. Here, the authors introduce RAMCES, a method that selects the optimal membrane markers to use for more accurate cell segmentation.
Collapse
|
36
|
Cancer: More than a geneticist’s Pandora’s box. J Biosci 2022. [DOI: 10.1007/s12038-022-00254-x] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
|
37
|
Simoni A, Huber HA, Georgia SK, Finley SD. Phosphatases are predicted to govern prolactin-mediated JAK–STAT signaling in pancreatic beta cells. Integr Biol (Camb) 2022; 14:37-48. [DOI: 10.1093/intbio/zyac004] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/25/2021] [Revised: 02/18/2022] [Accepted: 02/21/2022] [Indexed: 11/12/2022]
Abstract
Abstract
Patients with diabetes are unable to produce a sufficient amount of insulin to properly regulate their blood glucose levels. One potential method of treating diabetes is to increase the number of insulin-secreting beta cells in the pancreas to enhance insulin secretion. It is known that during pregnancy, pancreatic beta cells proliferate in response to the pregnancy hormone, prolactin (PRL). Leveraging this proliferative response to PRL may be a strategy to restore endogenous insulin production for patients with diabetes. To investigate this potential treatment, we previously developed a computational model to represent the PRL-mediated JAK–STAT signaling pathway in pancreatic beta cells. Here, we applied the model to identify the importance of particular signaling proteins in shaping the response of a population of beta cells. We simulated a population of 10 000 heterogeneous cells with varying initial protein concentrations responding to PRL stimulation. We used partial least squares regression to analyze the significance and role of each of the varied protein concentrations in producing the response of the cell. Our regression models predict that the concentrations of the cytosolic and nuclear phosphatases strongly influence the response of the cell. The model also predicts that increasing PRL receptor strengthens negative feedback mediated by the inhibitor suppressor of cytokine signaling. These findings reveal biological targets that can potentially be used to modulate the proliferation of pancreatic beta cells to enhance insulin secretion and beta cell regeneration in the context of diabetes.
Collapse
Affiliation(s)
- Ariella Simoni
- Department of Biomedical Engineering, University of Southern California, Los Angeles, CA 90089, USA
| | - Holly A Huber
- Department of Biomedical Engineering, University of Southern California, Los Angeles, CA 90089, USA
| | - Senta K Georgia
- Departments of Pediatrics and Stem Cell Biology and Regenerative Medicine, Keck School of Medicine, University of Southern California, Los Angeles, CA 90089, USA
| | - Stacey D Finley
- Department of Biomedical Engineering, University of Southern California, Los Angeles, CA 90089, USA
- Department of Quantitative and Computational Biology, University of Southern California, Los Angeles, CA 90089, USA
- Mork Family Department of Chemical Engineering and Materials Science, University of Southern California, Los Angeles, CA 90089, USA
| |
Collapse
|
38
|
Ipiña EP, Camley BA. Collective gradient sensing with limited positional information. Phys Rev E 2022; 105:044410. [PMID: 35590664 DOI: 10.1103/physreve.105.044410] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/05/2021] [Accepted: 03/21/2022] [Indexed: 06/15/2023]
Abstract
Eukaryotic cells sense chemical gradients to decide where and when to move. Clusters of cells can sense gradients more accurately than individual cells by integrating measurements of the concentration made across the cluster. Is this gradient-sensing accuracy impeded when cells have limited knowledge of their position within the cluster, i.e., limited positional information? We apply maximum likelihood estimation to study gradient-sensing accuracy of a cluster of cells with finite positional information. If cells must estimate their location within the cluster, this lowers the accuracy of collective gradient sensing. We compare our results with a tug-of-war model where cells respond to the gradient by polarizing away from their neighbors without relying on their positional information. As the cell positional uncertainty increases, there is a trade-off where the tug-of-war model responds more accurately to the chemical gradient. However, for sufficiently large cell clusters or sufficiently shallow chemical gradients, the tug-of-war model will always be suboptimal to one that integrates information from all cells, even if positional uncertainty is high.
Collapse
Affiliation(s)
- Emiliano Perez Ipiña
- Department of Physics & Astronomy, Johns Hopkins University, Baltimore, Maryland 21218, USA
| | - Brian A Camley
- Department of Physics & Astronomy and Biophysics, Johns Hopkins University, Baltimore, Maryland 21218, USA
| |
Collapse
|
39
|
Wu G, Xiu H, Luo H, Ding Y, Li Y. A mathematical model for cell cycle control: graded response or quantized response. Cell Cycle 2022; 21:820-834. [PMID: 35107036 PMCID: PMC8973363 DOI: 10.1080/15384101.2022.2031770] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/17/2021] [Revised: 01/03/2022] [Accepted: 01/17/2022] [Indexed: 02/04/2023] Open
Abstract
Cell cycle is an important and complex biological system. A lot of efforts have been put in understanding cell cycle arrest for its vital role in clinical therapies. The cell-cycle-arrest outcomes upon stimulation are complicated. The response could be stringent or relaxed, and graded or quantized. A model fully addressing various cell-cycle-arrest outcomes is to be developed. Here, we developed a mathematical model of cell cycle control incorporating distinct characteristics of various cell-cycle-arrest outcomes. The model can simulate two typical properties of cell cycle arrest, quantized and graded. We also characterized the inheritable quiescence and refractory state, which were crucial in long-term response of the population. Then, we monitored cells respond to multiple stimulations, and the results indicated that cells responded to stimulations with small interval did not induce significantly sustained cell cycle arrest as the existence of refractory state. Our work will benefit fundamental research and make efforts to predicting outcomes of clinical therapeutics.
Collapse
Affiliation(s)
- Guoyu Wu
- School of Clinical Pharmacy, Guangdong Pharmaceutical University, Guangdong, China
- Key Specialty of Clinical Pharmacy, The First Affiliated Hospital of Guangdong Pharmaceutical University, Guangzhou, China
- NMPA Key Laboratory for Technology Research and Evaluation of Pharmacovigilance, Guangdong Pharmaceutical University, Guangzhou, China
- CONTACT Guoyu Wu
| | - Huiyu Xiu
- School of Clinical Pharmacy, Guangdong Pharmaceutical University, Guangdong, China
| | - Haiying Luo
- School of Clinical Pharmacy, Guangdong Pharmaceutical University, Guangdong, China
| | - Yu Ding
- School of Clinical Pharmacy, Guangdong Pharmaceutical University, Guangdong, China
| | - Yuchao Li
- MegaLab, MegaRobo Technologies Co., Ltd, Beijing, China
- Yuchao Li
| |
Collapse
|
40
|
Howard GR, Jost TA, Yankeelov TE, Brock A. Quantification of long-term doxorubicin response dynamics in breast cancer cell lines to direct treatment schedules. PLoS Comput Biol 2022; 18:e1009104. [PMID: 35358172 PMCID: PMC9004764 DOI: 10.1371/journal.pcbi.1009104] [Citation(s) in RCA: 13] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/16/2021] [Revised: 04/12/2022] [Accepted: 02/07/2022] [Indexed: 01/05/2023] Open
Abstract
While acquired chemoresistance is recognized as a key challenge to treating many types of cancer, the dynamics with which drug sensitivity changes after exposure are poorly characterized. Most chemotherapeutic regimens call for repeated dosing at regular intervals, and if drug sensitivity changes on a similar time scale then the treatment interval could be optimized to improve treatment performance. Theoretical work suggests that such optimal schedules exist, but experimental confirmation has been obstructed by the difficulty of deconvolving the simultaneous processes of death, adaptation, and regrowth taking place in cancer cell populations. Here we present a method of optimizing drug schedules in vitro through iterative application of experimentally calibrated models, and demonstrate its ability to characterize dynamic changes in sensitivity to the chemotherapeutic doxorubicin in three breast cancer cell lines subjected to treatment schedules varying in concentration, interval between pulse treatments, and number of sequential pulse treatments. Cell populations are monitored longitudinally through automated imaging for 600–800 hours, and this data is used to calibrate a family of cancer growth models, each consisting of a system of ordinary differential equations, derived from the bi-exponential model which characterizes resistant and sensitive subpopulations. We identify a model incorporating both a period of growth arrest in surviving cells and a delay in the death of chemosensitive cells which outperforms the original bi-exponential growth model in Akaike Information Criterion based model selection, and use the calibrated model to quantify the performance of each drug schedule. We find that the inter-treatment interval is a key variable in determining the performance of sequential dosing schedules and identify an optimal retreatment time for each cell line which extends regrowth time by 40%-239%, demonstrating that the time scale of changes in chemosensitivity following doxorubicin exposure allows optimization of drug scheduling by varying this inter-treatment interval. Acquired chemoresistance is a common cause of treatment failure in cancer. The scheduling of a multi-dose course of chemotherapeutic treatment may influence the dynamics of acquired chemoresistance, and drug schedule optimization may increase the duration of effectiveness of a particular chemotherapeutic agent for a particular patient. Here we present a method for experimentally optimizing an in vitro drug schedule through iterative rounds of experimentation and computational analysis, and demonstrate the method’s ability to improve the performance of doxorubicin treatment in three breast carcinoma cell lines. Specifically, we find that the interval between drug exposures can be optimized while holding drug concentration and number of treatments constant, suggesting that this may be a key variable to explore in future drug schedule optimization efforts. We further use this method’s model calibration and selection process to extract information about the underlying biology of the doxorubicin response, and find that the incorporation of delays on both cell death and regrowth are necessary for accurate parameterization of cell growth data.
Collapse
Affiliation(s)
- Grant R. Howard
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, Texas, United States of America
| | - Tyler A. Jost
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, Texas, United States of America
| | - Thomas E. Yankeelov
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, Texas, United States of America
- Department of Diagnostic Medicine, The University of Texas at Austin, Austin, Texas, United States of America
- Department of Oncology, The University of Texas at Austin, Austin, Texas, United States of America
- Oden Institute for Computational Engineering and Sciences, The University of Texas at Austin, Austin, Texas, United States of America
- Livestrong Cancer Institutes, The University of Texas at Austin, Austin, Texas, United States of America
- Department of Imaging Physics, The University of Texas MD Anderson Cancer Center, Houston, Texas, United States of America
| | - Amy Brock
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, Texas, United States of America
- Department of Cell and Molecular Biology, The University of Texas at Austin, Austin, Texas, United States of America
- * E-mail:
| |
Collapse
|
41
|
Jain P, Bhatia S, Thompson EW, Jolly MK. Population Dynamics of Epithelial-Mesenchymal Heterogeneity in Cancer Cells. Biomolecules 2022; 12:biom12030348. [PMID: 35327538 PMCID: PMC8945776 DOI: 10.3390/biom12030348] [Citation(s) in RCA: 10] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/21/2022] [Revised: 02/17/2022] [Accepted: 02/18/2022] [Indexed: 11/16/2022] Open
Abstract
Phenotypic heterogeneity is a hallmark of aggressive cancer behaviour and a clinical challenge. Despite much characterisation of this heterogeneity at a multi-omics level in many cancers, we have a limited understanding of how this heterogeneity emerges spontaneously in an isogenic cell population. Some longitudinal observations of dynamics in epithelial-mesenchymal heterogeneity, a canonical example of phenotypic heterogeneity, have offered us opportunities to quantify the rates of phenotypic switching that may drive such heterogeneity. Here, we offer a mathematical modeling framework that explains the salient features of population dynamics noted in PMC42-LA cells: (a) predominance of EpCAMhigh subpopulation, (b) re-establishment of parental distributions from the EpCAMhigh and EpCAMlow subpopulations, and (c) enhanced heterogeneity in clonal populations established from individual cells. Our framework proposes that fluctuations or noise in content duplication and partitioning of SNAIL—an EMT-inducing transcription factor—during cell division can explain spontaneous phenotypic switching and consequent dynamic heterogeneity in PMC42-LA cells observed experimentally at both single-cell and bulk level analysis. Together, we propose that asymmetric cell division can be a potential mechanism for phenotypic heterogeneity.
Collapse
Affiliation(s)
- Paras Jain
- Centre for BioSystems Science and Engineering, Indian Institute of Science, Bangalore 560012, India;
| | - Sugandha Bhatia
- School of Biomedical Sciences, Faculty of Health, Queensland University of Technology, Brisbane 4000, Australia;
- The University of Queensland Diamantina Institute, Faculty of Medicine, The University of Queensland, Woolloongabba 4102, Australia
- Translational Research Institute, Woolloongabba 4102, Australia
| | - Erik W. Thompson
- School of Biomedical Sciences, Faculty of Health, Queensland University of Technology, Brisbane 4000, Australia;
- Translational Research Institute, Woolloongabba 4102, Australia
- Correspondence: (E.W.T.); (M.K.J.)
| | - Mohit Kumar Jolly
- Centre for BioSystems Science and Engineering, Indian Institute of Science, Bangalore 560012, India;
- Correspondence: (E.W.T.); (M.K.J.)
| |
Collapse
|
42
|
Cele S, Jackson L, Khoury DS, Khan K, Moyo-Gwete T, Tegally H, San JE, Cromer D, Scheepers C, Amoako DG, Karim F, Bernstein M, Lustig G, Archary D, Smith M, Ganga Y, Jule Z, Reedoy K, Hwa SH, Giandhari J, Blackburn JM, Gosnell BI, Abdool Karim SS, Hanekom W, von Gottberg A, Bhiman JN, Lessells RJ, Moosa MYS, Davenport MP, de Oliveira T, Moore PL, Sigal A. Omicron extensively but incompletely escapes Pfizer BNT162b2 neutralization. Nature 2022; 602:654-656. [PMID: 35016196 PMCID: PMC8866126 DOI: 10.1038/s41586-021-04387-1] [Citation(s) in RCA: 781] [Impact Index Per Article: 390.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/14/2021] [Accepted: 12/23/2021] [Indexed: 11/09/2022]
Abstract
The emergence of the SARS-CoV-2 variant of concern Omicron (Pango lineage B.1.1.529), first identified in Botswana and South Africa, may compromise vaccine effectiveness and lead to re-infections1. Here we investigated Omicron escape from neutralization by antibodies from South African individuals vaccinated with Pfizer BNT162b2. We used blood samples taken soon after vaccination from individuals who were vaccinated and previously infected with SARS-CoV-2 or vaccinated with no evidence of previous infection. We isolated and sequence-confirmed live Omicron virus from an infected person and observed that Omicron requires the angiotensin-converting enzyme 2 (ACE2) receptor to infect cells. We compared plasma neutralization of Omicron relative to an ancestral SARS-CoV-2 strain and found that neutralization of ancestral virus was much higher in infected and vaccinated individuals compared with the vaccinated-only participants. However, both groups showed a 22-fold reduction in vaccine-elicited neutralization by the Omicron variant. Participants who were vaccinated and had previously been infected exhibited residual neutralization of Omicron similar to the level of neutralization of the ancestral virus observed in the vaccination-only group. These data support the notion that reasonable protection against Omicron may be maintained using vaccination approaches.
Collapse
Affiliation(s)
- Sandile Cele
- Africa Health Research Institute, Durban, South Africa
- School of Laboratory Medicine and Medical Sciences, University of KwaZulu-Natal, Durban, South Africa
| | | | - David S Khoury
- Kirby Institute, University of New South Wales, Sydney, New South Wales, Australia
| | - Khadija Khan
- Africa Health Research Institute, Durban, South Africa
- School of Laboratory Medicine and Medical Sciences, University of KwaZulu-Natal, Durban, South Africa
| | - Thandeka Moyo-Gwete
- National Institute for Communicable Diseases of the National Health Laboratory Service, Johannesburg, South Africa
- SA MRC Antibody Immunity Research Unit, School of Pathology, Faculty of Health Sciences, University of the Witwatersrand, Johannesburg, South Africa
| | - Houriiyah Tegally
- KwaZulu-Natal Research Innovation and Sequencing Platform, Durban, South Africa
- Centre for Epidemic Response and Innovation, School of Data Science and Computational Thinking, Stellenbosch University, Stellenbosch, South Africa
| | - James Emmanuel San
- KwaZulu-Natal Research Innovation and Sequencing Platform, Durban, South Africa
| | - Deborah Cromer
- Kirby Institute, University of New South Wales, Sydney, New South Wales, Australia
| | - Cathrine Scheepers
- National Institute for Communicable Diseases of the National Health Laboratory Service, Johannesburg, South Africa
- SA MRC Antibody Immunity Research Unit, School of Pathology, Faculty of Health Sciences, University of the Witwatersrand, Johannesburg, South Africa
| | - Daniel G Amoako
- School of Laboratory Medicine and Medical Sciences, University of KwaZulu-Natal, Durban, South Africa
- National Institute for Communicable Diseases of the National Health Laboratory Service, Johannesburg, South Africa
| | - Farina Karim
- Africa Health Research Institute, Durban, South Africa
- School of Laboratory Medicine and Medical Sciences, University of KwaZulu-Natal, Durban, South Africa
| | | | - Gila Lustig
- Centre for the AIDS Programme of Research in South Africa, Durban, South Africa
| | - Derseree Archary
- Centre for the AIDS Programme of Research in South Africa, Durban, South Africa
- Department of Medical Microbiology, University of KwaZulu-Natal, Durban, South Africa
| | - Muneerah Smith
- Department of Integrative Biomedical Sciences, Faculty of Health Sciences, University of Cape Town, Cape Town, South Africa
| | - Yashica Ganga
- Africa Health Research Institute, Durban, South Africa
| | - Zesuliwe Jule
- Africa Health Research Institute, Durban, South Africa
| | - Kajal Reedoy
- Africa Health Research Institute, Durban, South Africa
| | - Shi-Hsia Hwa
- Africa Health Research Institute, Durban, South Africa
- Division of Infection and Immunity, University College London, London, UK
| | - Jennifer Giandhari
- KwaZulu-Natal Research Innovation and Sequencing Platform, Durban, South Africa
| | - Jonathan M Blackburn
- Department of Integrative Biomedical Sciences, Faculty of Health Sciences, University of Cape Town, Cape Town, South Africa
- Institute of Infectious Disease and Molecular Medicine, University of Cape Town, Cape Town, South Africa
| | - Bernadett I Gosnell
- Department of Infectious Diseases, Nelson R. Mandela School of Clinical Medicine, University of KwaZulu-Natal, Durban, South Africa
| | - Salim S Abdool Karim
- Centre for the AIDS Programme of Research in South Africa, Durban, South Africa
- Department of Epidemiology, Mailman School of Public Health, Columbia University, New York, NY, USA
| | - Willem Hanekom
- Africa Health Research Institute, Durban, South Africa
- Division of Infection and Immunity, University College London, London, UK
| | - Anne von Gottberg
- National Institute for Communicable Diseases of the National Health Laboratory Service, Johannesburg, South Africa
- SA MRC Antibody Immunity Research Unit, School of Pathology, Faculty of Health Sciences, University of the Witwatersrand, Johannesburg, South Africa
| | - Jinal N Bhiman
- National Institute for Communicable Diseases of the National Health Laboratory Service, Johannesburg, South Africa
- SA MRC Antibody Immunity Research Unit, School of Pathology, Faculty of Health Sciences, University of the Witwatersrand, Johannesburg, South Africa
| | - Richard J Lessells
- KwaZulu-Natal Research Innovation and Sequencing Platform, Durban, South Africa
- Centre for the AIDS Programme of Research in South Africa, Durban, South Africa
| | - Mahomed-Yunus S Moosa
- Department of Infectious Diseases, Nelson R. Mandela School of Clinical Medicine, University of KwaZulu-Natal, Durban, South Africa
| | - Miles P Davenport
- Kirby Institute, University of New South Wales, Sydney, New South Wales, Australia
| | - Tulio de Oliveira
- KwaZulu-Natal Research Innovation and Sequencing Platform, Durban, South Africa
- Centre for Epidemic Response and Innovation, School of Data Science and Computational Thinking, Stellenbosch University, Stellenbosch, South Africa
- Centre for the AIDS Programme of Research in South Africa, Durban, South Africa
- Department of Global Health, University of Washington, Seattle, WA, USA
| | - Penny L Moore
- National Institute for Communicable Diseases of the National Health Laboratory Service, Johannesburg, South Africa
- SA MRC Antibody Immunity Research Unit, School of Pathology, Faculty of Health Sciences, University of the Witwatersrand, Johannesburg, South Africa
- Centre for the AIDS Programme of Research in South Africa, Durban, South Africa
- Institute of Infectious Disease and Molecular Medicine, University of Cape Town, Cape Town, South Africa
| | - Alex Sigal
- Africa Health Research Institute, Durban, South Africa.
- School of Laboratory Medicine and Medical Sciences, University of KwaZulu-Natal, Durban, South Africa.
- Max Planck Institute for Infection Biology, Berlin, Germany.
| |
Collapse
|
43
|
Wong DCS, Seinkmane E, Zeng A, Stangherlin A, Rzechorzek NM, Beale AD, Day J, Reed M, Peak‐Chew SY, Styles CT, Edgar RS, Putker M, O’Neill JS. CRYPTOCHROMES promote daily protein homeostasis. EMBO J 2022; 41:e108883. [PMID: 34842284 PMCID: PMC8724739 DOI: 10.15252/embj.2021108883] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/03/2021] [Revised: 11/07/2021] [Accepted: 11/09/2021] [Indexed: 11/29/2022] Open
Abstract
The daily organisation of most mammalian cellular functions is attributed to circadian regulation of clock-controlled protein expression, driven by daily cycles of CRYPTOCHROME-dependent transcriptional feedback repression. To test this, we used quantitative mass spectrometry to compare wild-type and CRY-deficient fibroblasts under constant conditions. In CRY-deficient cells, we found that temporal variation in protein, phosphopeptide, and K+ abundance was at least as great as wild-type controls. Most strikingly, the extent of temporal variation within either genotype was much smaller than overall differences in proteome composition between WT and CRY-deficient cells. This proteome imbalance in CRY-deficient cells and tissues was associated with increased susceptibility to proteotoxic stress, which impairs circadian robustness, and may contribute to the wide-ranging phenotypes of CRY-deficient mice. Rather than generating large-scale daily variation in proteome composition, we suggest it is plausible that the various transcriptional and post-translational functions of CRY proteins ultimately act to maintain protein and osmotic homeostasis against daily perturbation.
Collapse
Affiliation(s)
| | | | - Aiwei Zeng
- MRC Laboratory of Molecular BiologyCambridgeUK
| | | | | | | | - Jason Day
- Department of Earth SciencesUniversity of CambridgeCambridgeUK
| | - Martin Reed
- MRC Laboratory of Molecular BiologyCambridgeUK
| | | | | | - Rachel S Edgar
- Department of Infectious DiseasesImperial CollegeLondonUK
| | - Marrit Putker
- MRC Laboratory of Molecular BiologyCambridgeUK
- Present address:
Crown BioscienceUtrechtthe Netherlands
| | | |
Collapse
|
44
|
Wu G, Li Y. Distinct characteristics of correlation analysis at the single-cell and the population level. Stat Appl Genet Mol Biol 2022; 21:sagmb-2022-0015. [PMID: 35918809 DOI: 10.1515/sagmb-2022-0015] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/27/2022] [Accepted: 06/13/2022] [Indexed: 11/15/2022]
Abstract
Correlation analysis is widely used in biological studies to infer molecular relationships within biological networks. Recently, single-cell analysis has drawn tremendous interests, for its ability to obtain high-resolution molecular phenotypes. It turns out that there is little overlap of co-expressed genes identified in single-cell level investigations with that of population level investigations. However, the nature of the relationship of correlations between single-cell and population levels remains unclear. In this manuscript, we aimed to unveil the origin of the differences between the correlation coefficients at the single-cell level and that at the population level, and bridge the gap between them. Through developing formulations to link correlations at the single-cell and the population level, we illustrated that aggregated correlations could be stronger, weaker or equal to the corresponding individual correlations, depending on the variations and the correlations within the population. When the correlation within the population is weaker than the individual correlation, the aggregated correlation is stronger than the corresponding individual correlation. Besides, our data indicated that aggregated correlation is more likely to be stronger than the corresponding individual correlation, and it was rare to find gene-pairs exclusively strongly correlated at the single-cell level. Through a bottom-up approach to model interactions between molecules in a signaling cascade or a multi-regulator-controlled gene expression, we surprisingly found that the existence of interaction between two components could not be excluded simply based on their low correlation coefficients, suggesting a reconsideration of connectivity within biological networks which was derived solely from correlation analysis. We also investigated the impact of technical random measurement errors on the correlation coefficients for the single-cell level and the population level. The results indicate that the aggregated correlation is relatively robust and less affected. Because of the heterogeneity among single cells, correlation coefficients calculated based on data of the single-cell level might be different from that of the population level. Depending on the specific question we are asking, proper sampling and normalization procedure should be done before we draw any conclusions.
Collapse
Affiliation(s)
- Guoyu Wu
- School of Clinical Pharmacy, Guangdong Pharmaceutical University, Guangzhou, China
- Key Specialty of Clinical Pharmacy, The First Affiliated Hospital of Guangdong Pharmaceutical University, Guangzhou, China
- NMPA Key Laboratory for Technology Research and Evaluation of Pharmacovigilance, Guangdong Pharmaceutical University, Guangzhou, China
| | - Yuchao Li
- Max Planck Institute for Molecular Genetics, Berlin, Germany
- MegaLab, MegaRobo Technologies Co., Ltd, Beijing, China
| |
Collapse
|
45
|
Kim R, Witelski TP. Uncovering the dynamics of a circadian-dopamine model influenced by the light-dark cycle. Math Biosci 2021; 344:108764. [PMID: 34952036 DOI: 10.1016/j.mbs.2021.108764] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/28/2021] [Revised: 11/09/2021] [Accepted: 11/26/2021] [Indexed: 10/19/2022]
Abstract
The neurotransmitter dopamine (DA) is known to be influenced by the circadian timekeeping system in the mammalian brain. We have previously created a single-cell differential equations model to understand the mechanisms behind circadian rhythms of extracellular DA. In this paper, we investigate the dynamics in our model and study different behaviors such as entrainment to the 24-hour light-dark cycle and robust periodicity versus decoupling, quasiperiodicity, and chaos. Imbalances in DA are often accompanied by disrupted circadian rhythms, such as in Parkinson's disease, hyperactivity, and mood disorders. Our model provides new insights into the links between the circadian clock and DA. We show that the daily rhythmicity of DA can be disrupted by decoupling between interlocked loops of the clock circuitry or by quasiperiodic clock behaviors caused by misalignment with the light-dark cycle. The model can be used to further study how the circadian clock affects the dopaminergic system, and to help develop therapeutic strategies for disrupted DA rhythms. .
Collapse
Affiliation(s)
- Ruby Kim
- Department of Mathematics, Duke University, Durham, NC, USA.
| | | |
Collapse
|
46
|
Kim R, Nijhout HF, Reed MC. One-carbon metabolism during the menstrual cycle and pregnancy. PLoS Comput Biol 2021; 17:e1009708. [PMID: 34914693 PMCID: PMC8741061 DOI: 10.1371/journal.pcbi.1009708] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/22/2021] [Revised: 01/07/2022] [Accepted: 12/01/2021] [Indexed: 11/18/2022] Open
Abstract
Many enzymes in one-carbon metabolism (OCM) are up- or down-regulated by the sex hormones which vary diurnally and throughout the menstrual cycle. During pregnancy, estradiol and progesterone levels increase tremendously to modulate physiological changes in the reproductive system. In this work, we extend and improve an existing mathematical model of hepatic OCM to understand the dynamic metabolic changes that happen during the menstrual cycle and pregnancy due to estradiol variation. In particular, we add the polyamine drain on S-adenosyl methionine and the direct effects of estradiol on the enzymes cystathionine β-synthase (CBS), thymidylate synthase (TS), and dihydrofolate reductase (DHFR). We show that the homocysteine concentration varies inversely with estradiol concentration, discuss the fluctuations in 14 other one-carbon metabolites and velocities throughout the menstrual cycle, and draw comparisons with the literature. We then use the model to study the effects of vitamin B12, vitamin B6, and folate deficiencies and explain why homocysteine is not a good biomarker for vitamin deficiencies. Additionally, we compute homocysteine throughout pregnancy, and compare the results with experimental data. Our mathematical model explains how numerous homeostatic mechanisms in OCM function and provides new insights into how homocysteine and its deleterious effects are influenced by estradiol. The mathematical model can be used by others for further in silico experiments on changes in one-carbon metabolism during the menstrual cycle and pregnancy.
Collapse
Affiliation(s)
- Ruby Kim
- Department of Mathematics, Duke University, Durham, North Carolina, United States of America
| | - H. Frederik Nijhout
- Department of Biology, Duke University, Durham, North Carolina, United States of America
| | - Michael C. Reed
- Department of Mathematics, Duke University, Durham, North Carolina, United States of America
| |
Collapse
|
47
|
Peckys DB, Gaa D, de Jonge N. Quantification of EGFR-HER2 Heterodimers in HER2-Overexpressing Breast Cancer Cells Using Liquid-Phase Electron Microscopy. Cells 2021; 10:cells10113244. [PMID: 34831465 PMCID: PMC8623301 DOI: 10.3390/cells10113244] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/12/2021] [Revised: 11/10/2021] [Accepted: 11/15/2021] [Indexed: 12/25/2022] Open
Abstract
Currently, breast cancer patients are classified uniquely according to the expression level of hormone receptors, and human epidermal growth factor receptor 2 (HER2). This coarse classification is insufficient to capture the phenotypic complexity and heterogeneity of the disease. A methodology was developed for absolute quantification of receptor surface density ρR, and molecular interaction (dimerization), as well as the associated heterogeneities, of HER2 and its family member, the epidermal growth factor receptor (EGFR) in the plasma membrane of HER2 overexpressing breast cancer cells. Quantitative, correlative light microscopy (LM) and liquid-phase electron microscopy (LPEM) were combined with quantum dot (QD) labeling. Single-molecule position data of receptors were obtained from scanning transmission electron microscopy (STEM) images of intact cancer cells. Over 280,000 receptor positions were detected and statistically analyzed. An important finding was the subcellular heterogeneity in heterodimer shares with respect to plasma membrane regions with different dynamic properties. Deriving quantitative information about EGFR and HER2 ρR, as well as their dimer percentages, and the heterogeneities thereof, in single cancer cells, is potentially relevant for early identification of patients with HER2 overexpressing tumors comprising an enhanced share of EGFR dimers, likely increasing the risk for drug resistance, and thus requiring additional targeted therapeutic strategies.
Collapse
Affiliation(s)
- Diana B. Peckys
- Clinic of Operative Dentistry, Periodontology and Preventive Dentistry, University Hospital, Saarland University, 66421 Homburg, Germany;
| | - Daniel Gaa
- INM—Leibniz Institute for New Materials, 66123 Saarbrücken, Germany;
| | - Niels de Jonge
- INM—Leibniz Institute for New Materials, 66123 Saarbrücken, Germany;
- Department of Physics, Saarland University, 66123 Saarbrücken, Germany
- Correspondence:
| |
Collapse
|
48
|
Aguado-García A, Priego-Espinosa DA, Aldana A, Darszon A, Martínez-Mekler G. Mathematical model reveals that heterogeneity in the number of ion transporters regulates the fraction of mouse sperm capacitation. PLoS One 2021; 16:e0245816. [PMID: 34793454 PMCID: PMC8601445 DOI: 10.1371/journal.pone.0245816] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/26/2020] [Accepted: 10/20/2021] [Indexed: 12/03/2022] Open
Abstract
Capacitation is a complex maturation process mammalian sperm must undergo in the female genital tract to be able to fertilize an egg. This process involves, amongst others, physiological changes in flagellar beating pattern, membrane potential, intracellular ion concentrations and protein phosphorylation. Typically, in a capacitation medium, only a fraction of sperm achieve this state. The cause for this heterogeneous response is still not well understood and remains an open question. Here, one of our principal results is to develop a discrete regulatory network, with mostly deterministic dynamics in conjunction with some stochastic elements, for the main biochemical and biophysical processes involved in the early events of capacitation. The model criterion for capacitation requires the convergence of specific levels of a select set of nodes. Besides reproducing several experimental results and providing some insight on the network interrelations, the main contribution of the model is the suggestion that the degree of variability in the total amount and individual number of ion transporters among spermatozoa regulates the fraction of capacitated spermatozoa. This conclusion is consistent with recently reported experimental results. Based on this mathematical analysis, experimental clues are proposed for the control of capacitation levels. Furthermore, cooperative and interference traits that become apparent in the modelling among some components also call for future theoretical and experimental studies.
Collapse
Affiliation(s)
- Alejandro Aguado-García
- Instituto de Ciencias Básicas y Aplicadas, Universidad Autónoma del Estado de Morelos, Cuernavaca, Morelos, México
- Instituto de Ciencias Físicas, Universidad Nacional Autónoma de México, Cuernavaca, Morelos, México
| | | | - Andrés Aldana
- Instituto de Biotecnología, Universidad Nacional Autónoma de México, Cuernavaca, Morelos, México
- Centro de Ciencias de la Complejidad, Universidad Nacional Autónoma de México, CDMX, México
| | - Alberto Darszon
- Instituto de Biotecnología, Universidad Nacional Autónoma de México, Cuernavaca, Morelos, México
| | - Gustavo Martínez-Mekler
- Instituto de Ciencias Físicas, Universidad Nacional Autónoma de México, Cuernavaca, Morelos, México
- Centro de Ciencias de la Complejidad, Universidad Nacional Autónoma de México, CDMX, México
| |
Collapse
|
49
|
Chaves M, Gomes-Pereira LC, Roux J. Two-level modeling approach to identify the regulatory dynamics capturing drug response heterogeneity in single-cells. Sci Rep 2021; 11:20809. [PMID: 34675364 PMCID: PMC8531316 DOI: 10.1038/s41598-021-99943-0] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/26/2021] [Accepted: 09/27/2021] [Indexed: 11/09/2022] Open
Abstract
Single-cell multimodal technologies reveal the scales of cellular heterogeneity impairing cancer treatment, yet cell response dynamics remain largely underused to decipher the mechanisms of drug resistance they take part in. As the phenotypic heterogeneity of a clonal cell population informs on the capacity of each single-cell to recapitulate the whole range of observed behaviors, we developed a modeling approach utilizing single-cell response data to identify regulatory reactions driving population heterogeneity in drug response. Dynamic data of hundreds of HeLa cells treated with TNF-related apoptosis-inducing ligand (TRAIL) were used to characterize the fate-determining kinetic parameters of an apoptosis receptor reaction model. Selected reactions sets were augmented to incorporate a mechanism that leads to the separation of the opposing response phenotypes. Using a positive feedback loop motif to identify the reaction set, we show that caspase-8 is able to encapsulate high levels of heterogeneity by introducing a response delay and amplifying the initial differences arising from natural protein expression variability. Our approach enables the identification of fate-determining reactions that drive the population response heterogeneity, providing regulatory targets to curb the cell dynamics of drug resistance.
Collapse
Affiliation(s)
- Madalena Chaves
- Université Côte d'Azur, Inria, INRAE, CNRS, Sorbonne Université, Biocore Team, Sophia Antipolis, France
| | - Luis C Gomes-Pereira
- Université Côte d'Azur, Inria, INRAE, CNRS, Sorbonne Université, Biocore Team, Sophia Antipolis, France.,Université Côte d'Azur, CNRS UMR 7284, Inserm U 1081, Institut de Recherche sur le Cancer et le Vieillissement de Nice, Centre Antoine Lacassagne, 06107, Nice, France
| | - Jérémie Roux
- Université Côte d'Azur, CNRS UMR 7284, Inserm U 1081, Institut de Recherche sur le Cancer et le Vieillissement de Nice, Centre Antoine Lacassagne, 06107, Nice, France.
| |
Collapse
|
50
|
Adlung L, Stapor P, Tönsing C, Schmiester L, Schwarzmüller LE, Postawa L, Wang D, Timmer J, Klingmüller U, Hasenauer J, Schilling M. Cell-to-cell variability in JAK2/STAT5 pathway components and cytoplasmic volumes defines survival threshold in erythroid progenitor cells. Cell Rep 2021; 36:109507. [PMID: 34380040 DOI: 10.1016/j.celrep.2021.109507] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/06/2020] [Revised: 04/01/2021] [Accepted: 07/19/2021] [Indexed: 12/25/2022] Open
Abstract
Survival or apoptosis is a binary decision in individual cells. However, at the cell-population level, a graded increase in survival of colony-forming unit-erythroid (CFU-E) cells is observed upon stimulation with erythropoietin (Epo). To identify components of Janus kinase 2/signal transducer and activator of transcription 5 (JAK2/STAT5) signal transduction that contribute to the graded population response, we extended a cell-population-level model calibrated with experimental data to study the behavior in single cells. The single-cell model shows that the high cell-to-cell variability in nuclear phosphorylated STAT5 is caused by variability in the amount of Epo receptor (EpoR):JAK2 complexes and of SHP1, as well as the extent of nuclear import because of the large variance in the cytoplasmic volume of CFU-E cells. 24-118 pSTAT5 molecules in the nucleus for 120 min are sufficient to ensure cell survival. Thus, variability in membrane-associated processes is sufficient to convert a switch-like behavior at the single-cell level to a graded population-level response.
Collapse
Affiliation(s)
- Lorenz Adlung
- Division Systems Biology of Signal Transduction, German Cancer Research Center (DKFZ), 69120 Heidelberg, Germany; Department of Medicine, University Medical Center Hamburg-Eppendorf, 20246 Hamburg, Germany; Hamburg Center for Translational Immunology (HCTI), University Medical Center Hamburg-Eppendorf, 20246 Hamburg, Germany
| | - Paul Stapor
- Helmholtz Zentrum München - German Research Center for Environmental Health, Institute of Computational Biology, 85764 Neuherberg, Germany; Technische Universität München, Center for Mathematics, Chair of Mathematical Modeling of Biological Systems, 85748 Garching, Germany
| | - Christian Tönsing
- Institute of Physics and Freiburg Center for Data Analysis and Modelling (FDM), University of Freiburg, 79104 Freiburg, Germany; CIBSS-Centre for Integrative Biological Signalling Studies, University of Freiburg, 79104 Freiburg, Germany
| | - Leonard Schmiester
- Helmholtz Zentrum München - German Research Center for Environmental Health, Institute of Computational Biology, 85764 Neuherberg, Germany; Technische Universität München, Center for Mathematics, Chair of Mathematical Modeling of Biological Systems, 85748 Garching, Germany
| | - Luisa E Schwarzmüller
- Division Systems Biology of Signal Transduction, German Cancer Research Center (DKFZ), 69120 Heidelberg, Germany
| | - Lena Postawa
- Division Systems Biology of Signal Transduction, German Cancer Research Center (DKFZ), 69120 Heidelberg, Germany
| | - Dantong Wang
- Helmholtz Zentrum München - German Research Center for Environmental Health, Institute of Computational Biology, 85764 Neuherberg, Germany; Technische Universität München, Center for Mathematics, Chair of Mathematical Modeling of Biological Systems, 85748 Garching, Germany
| | - Jens Timmer
- Institute of Physics and Freiburg Center for Data Analysis and Modelling (FDM), University of Freiburg, 79104 Freiburg, Germany; CIBSS-Centre for Integrative Biological Signalling Studies, University of Freiburg, 79104 Freiburg, Germany.
| | - Ursula Klingmüller
- Division Systems Biology of Signal Transduction, German Cancer Research Center (DKFZ), 69120 Heidelberg, Germany; Translational Lung Research Center (TLRC), Member of the German Center for Lung Research (DZL), 69120 Heidelberg, Germany.
| | - Jan Hasenauer
- Helmholtz Zentrum München - German Research Center for Environmental Health, Institute of Computational Biology, 85764 Neuherberg, Germany; Technische Universität München, Center for Mathematics, Chair of Mathematical Modeling of Biological Systems, 85748 Garching, Germany; Faculty of Mathematics and Natural Sciences, University of Bonn, 53113 Bonn, Germany.
| | - Marcel Schilling
- Division Systems Biology of Signal Transduction, German Cancer Research Center (DKFZ), 69120 Heidelberg, Germany.
| |
Collapse
|