1
|
McMahon‑Cole H, Johnson A, Sadat Aghamiri S, Helikar T, Crawford LB. Modeling and Remodeling the Cell: How Digital Twins and HCMV Can Elucidate the Complex Interactions of Viral Latency, Epigenetic Regulation, and Immune Responses. CURRENT CLINICAL MICROBIOLOGY REPORTS 2023; 10:141-151. [PMID: 37901689 PMCID: PMC10601359 DOI: 10.1007/s40588-023-00201-w] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 07/12/2023] [Indexed: 10/31/2023]
Abstract
Purpose of Review Human cytomegalovirus (HCMV), while asymptomatic in most, causes significant complications during fetal development, following transplant or in immunosuppressed individuals. The host-virus interactions regulating viral latency and reactivation and viral control of the cellular environment (immune regulation, differentiation, epigenetics) are highly complex. Understanding these processes is essential to controlling infection and can be leveraged as a novel approach for understanding basic cell biology. Recent Findings Immune digital twins (IDTs) are digital simulations integrating knowledge of human immunology, physiology, and patient-specific clinical data to predict individualized immune responses and targeted treatments. Recent studies used IDTs to elucidate mechanisms of T cells, dendritic cells, and epigenetic control-all key to HCMV biology. Summary Here, we discuss how leveraging the unique biology of HCMV and IDTs will clarify immune response dynamics, host-virus interactions, and viral latency and reactivation and serve as a powerful IDT-validation platform for individualized and holistic health management.
Collapse
Affiliation(s)
- Hana McMahon‑Cole
- Department of Biochemistry, University of Nebraska-Lincoln, Lincoln, NE, USA
| | - Alicia Johnson
- Department of Biochemistry, University of Nebraska-Lincoln, Lincoln, NE, USA
| | - Sara Sadat Aghamiri
- Department of Biochemistry, University of Nebraska-Lincoln, Lincoln, NE, USA
| | - Tomáš Helikar
- Department of Biochemistry, University of Nebraska-Lincoln, Lincoln, NE, USA
| | - Lindsey B. Crawford
- Department of Biochemistry, University of Nebraska-Lincoln, Lincoln, NE, USA
- Nebraska Center for Virology, Lincoln, NE, USA
- Nebraska Center for Integrated Biomolecular Communication, Lincoln, NE, USA
| |
Collapse
|
2
|
Mazaya M, Kwon YK. In Silico Pleiotropy Analysis in KEGG Signaling Networks Using a Boolean Network Model. Biomolecules 2022; 12:biom12081139. [PMID: 36009032 PMCID: PMC9406064 DOI: 10.3390/biom12081139] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/20/2022] [Revised: 08/10/2022] [Accepted: 08/15/2022] [Indexed: 11/16/2022] Open
Abstract
Pleiotropy, which refers to the ability of different mutations on the same gene to cause different pathological effects in human genetic diseases, is important in understanding system-level biological diseases. Although some biological experiments have been proposed, still little is known about pleiotropy on gene–gene dynamics, since most previous studies have been based on correlation analysis. Therefore, a new perspective is needed to investigate pleiotropy in terms of gene–gene dynamical characteristics. To quantify pleiotropy in terms of network dynamics, we propose a measure called in silico Pleiotropic Scores (sPS), which represents how much a gene is affected against a pair of different types of mutations on a Boolean network model. We found that our model can identify more candidate pleiotropic genes that are not known to be pleiotropic than the experimental database. In addition, we found that many types of functionally important genes tend to have higher sPS values than other genes; in other words, they are more pleiotropic. We investigated the relations of sPS with the structural properties in the signaling network and found that there are highly positive relations to degree, feedback loops, and centrality measures. This implies that the structural characteristics are principles to identify new pleiotropic genes. Finally, we found some biological evidence showing that sPS analysis is relevant to the real pleiotropic data and can be considered a novel candidate for pleiotropic gene research. Taken together, our results can be used to understand the dynamics pleiotropic characteristics in complex biological systems in terms of gene–phenotype relations.
Collapse
Affiliation(s)
- Maulida Mazaya
- Research Center for Computing, National Research and Innovation Agency (BRIN), Cibinong Science Center, Jl. Raya Jakarta-Bogor KM 46, Cibinong 16911, West Java, Indonesia
| | - Yung-Keun Kwon
- School of IT Convergence, University of Ulsan, 93 Daehak-ro, Nam-gu, Ulsan 44610, Korea
- Correspondence:
| |
Collapse
|
3
|
Boolean analysis of lateral inhibition. J Math Biol 2020; 81:463-486. [PMID: 32728826 PMCID: PMC7427764 DOI: 10.1007/s00285-020-01515-9] [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/09/2019] [Revised: 12/09/2019] [Indexed: 11/03/2022]
Abstract
We study Boolean networks which are simple spatial models of the highly conserved Delta-Notch system. The models assume the inhibition of Delta in each cell by Notch in the same cell, and the activation of Notch in presence of Delta in surrounding cells. We consider fully asynchronous dynamics over undirected graphs representing the neighbour relation between cells. In this framework, one can show that all attractors are fixed points for the system, independently of the neighbour relation, for instance by using known properties of simplified versions of the models, where only one species per cell is defined. The fixed points correspond to the so-called fine-grained "patterns" that emerge in discrete and continuous modelling of lateral inhibition. We study the reachability of fixed points, giving a characterisation of the trap spaces and the basins of attraction for both the full and the simplified models. In addition, we use a characterisation of the trap spaces to investigate the robustness of patterns to perturbations. The results of this qualitative analysis can complement and guide simulation-based approaches, and serve as a basis for the investigation of more complex mechanisms.
Collapse
|
4
|
Enciso J, Pelayo R, Villarreal C. From Discrete to Continuous Modeling of Lymphocyte Development and Plasticity in Chronic Diseases. Front Immunol 2019; 10:1927. [PMID: 31481957 PMCID: PMC6710364 DOI: 10.3389/fimmu.2019.01927] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/12/2018] [Accepted: 07/30/2019] [Indexed: 12/12/2022] Open
Abstract
The molecular events leading to differentiation, development, and plasticity of lymphoid cells have been subject of intense research due to their key roles in multiple pathologies, such as lymphoproliferative disorders, tumor growth maintenance and chronic diseases. The emergent roles of lymphoid cells and the use of high-throughput technologies have led to an extensive accumulation of experimental data allowing the reconstruction of gene regulatory networks (GRN) by integrating biochemical signals provided by the microenvironment with transcriptional modules of lineage-specific genes. Computational modeling of GRN has been useful for the identification of molecular switches involved in lymphoid specification, prediction of microenvironment-dependent cell plasticity, and analyses of signaling events occurring downstream the activation of antigen recognition receptors. Among most common modeling strategies to analyze the dynamical behavior of GRN, discrete dynamic models are widely used for their capacity to capture molecular interactions when a limited knowledge of kinetic parameters is present. However, they are less powerful when modeling complex systems sensitive to biochemical gradients. To compensate it, discrete models may be transformed into regulatory networks that includes state variables and parameters varying within a continuous range. This approach is based on a system of differential equations dynamics with regulatory interactions described by fuzzy logic propositions. Here, we discuss the applicability of this method on modeling of development and plasticity processes of adaptive lymphocytes, and its potential implications in the study of pathological landscapes associated to chronic diseases.
Collapse
Affiliation(s)
- Jennifer Enciso
- Centro de Investigación Biomédica de Oriente, Instituto Mexicano del Seguro Social, Mexico City, Mexico
- Programa de Doctorado en Ciencias Biomédicas, Universidad Nacional Autónoma de México, Mexico City, Mexico
- Centro de Ciencias de la Complejidad, Universidad Nacional Autónoma de México, Mexico City, Mexico
| | - Rosana Pelayo
- Centro de Investigación Biomédica de Oriente, Instituto Mexicano del Seguro Social, Mexico City, Mexico
| | - Carlos Villarreal
- Centro de Ciencias de la Complejidad, Universidad Nacional Autónoma de México, Mexico City, Mexico
- Departamento de Física Cuántica y Fotónica, Instituto de Física, Universidad Nacional Autónoma de México, Mexico City, Mexico
| |
Collapse
|
5
|
Abstract
Being concerned by the understanding of the mechanism underlying chronic degenerative diseases , we presented in the previous chapter the medical systems biology conceptual framework that we present for that purpose in this volume. More specifically, we argued there the clear advantages offered by a state-space perspective when applied to the systems-level description of the biomolecular machinery that regulates complex degenerative diseases. We also discussed the importance of the dynamical interplay between the risk factors and the network of interdependencies that characterizes the biochemical, cellular, and tissue-level biomolecular reactions that underlie the physiological processes in health and disease. As we pointed out in the previous chapter, the understanding of this interplay (articulated around cellular phenotypic plasticity properties, regulated by specific kinds of gene regulatory networks) is necessary if prevention is chosen as the human-health improvement strategy (potentially involving the modulation of the patient's lifestyle). In this chapter we provide the medical systems biology mathematical and computational modeling tools required for this task.
Collapse
|
6
|
Varela PL, Ramos CV, Monteiro PT, Chaouiya C. EpiLog: A software for the logical modelling of epithelial dynamics. F1000Res 2018; 7:1145. [PMID: 30363398 PMCID: PMC6173114 DOI: 10.12688/f1000research.15613.2] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Accepted: 03/01/2019] [Indexed: 01/12/2023] Open
Abstract
Cellular responses are governed by regulatory networks subject to external signals from surrounding cells and to other micro-environmental cues. The logical (Boolean or multi-valued) framework proved well suited to study such processes at the cellular level, by specifying qualitative models of involved signalling pathways and gene regulatory networks. Here, we describe and illustrate the main features of EpiLog, a computational tool that implements an extension of the logical framework to the tissue level. EpiLog defines a collection of hexagonal cells over a 2D grid, which embodies a mono-layer epithelium. Basically, it defines a cellular automaton in which cell behaviours are driven by associated logical models subject to external signals. EpiLog is freely available on the web at http://epilog-tool.org. It is implemented in Java (version ≥1.7 required) and the source code is provided at https://github.com/epilog-tool/epilog under a GNU General Public License v3.0.
Collapse
Affiliation(s)
- Pedro L Varela
- Instituto Superior Técnico, Universidade de Lisboa, Lisbon, P-1049-001, Portugal.,INESC-ID, Lisbon, P-1000-029, Portugal.,Instituto Gulbenkian de Ciência, Oeiras, P-2780-156, Portugal
| | - Camila V Ramos
- INESC-ID, Lisbon, P-1000-029, Portugal.,Instituto Gulbenkian de Ciência, Oeiras, P-2780-156, Portugal
| | - Pedro T Monteiro
- Instituto Superior Técnico, Universidade de Lisboa, Lisbon, P-1049-001, Portugal.,INESC-ID, Lisbon, P-1000-029, Portugal
| | - Claudine Chaouiya
- Instituto Gulbenkian de Ciência, Oeiras, P-2780-156, Portugal.,CNRS, Centrale Marseille, l'Institut de Mathématiques de Marseille, Aix-Marseille University, Marseille, France
| |
Collapse
|
7
|
Deritei D, Aird WC, Ercsey-Ravasz M, Regan ER. Principles of dynamical modularity in biological regulatory networks. Sci Rep 2016; 6:21957. [PMID: 26979940 PMCID: PMC4793241 DOI: 10.1038/srep21957] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/06/2015] [Accepted: 02/02/2016] [Indexed: 01/02/2023] Open
Abstract
Intractable diseases such as cancer are associated with breakdown in multiple individual functions, which conspire to create unhealthy phenotype-combinations. An important challenge is to decipher how these functions are coordinated in health and disease. We approach this by drawing on dynamical systems theory. We posit that distinct phenotype-combinations are generated by interactions among robust regulatory switches, each in control of a discrete set of phenotypic outcomes. First, we demonstrate the advantage of characterizing multi-switch regulatory systems in terms of their constituent switches by building a multiswitch cell cycle model which points to novel, testable interactions critical for early G2/M commitment to division. Second, we define quantitative measures of dynamical modularity, namely that global cell states are discrete combinations of switch-level phenotypes. Finally, we formulate three general principles that govern the way coupled switches coordinate their function.
Collapse
Affiliation(s)
- Dávid Deritei
- Hungarian Physics Institute, Faculty of Physics, Babes¸-Bolyai University, Cluj-Napoca 400084, Romania.,Center for Network Science, Central European University, Budapest, 1051, Hungary
| | - William C Aird
- Center for Vascular Biology Research, Department of Medicine, Beth Israel Deaconess Medical Center, Harvard Medical School, Boston, MA 02215, USA
| | - Mária Ercsey-Ravasz
- Hungarian Physics Institute, Faculty of Physics, Babes¸-Bolyai University, Cluj-Napoca 400084, Romania
| | - Erzsébet Ravasz Regan
- Center for Vascular Biology Research, Department of Medicine, Beth Israel Deaconess Medical Center, Harvard Medical School, Boston, MA 02215, USA.,Biochemistry and Molecular Biology, The College of Wooster, Wooster, OH 44691, USA
| |
Collapse
|
8
|
Trinh HC, Kwon YK. Effective Boolean dynamics analysis to identify functionally important genes in large-scale signaling networks. Biosystems 2015; 137:64-72. [DOI: 10.1016/j.biosystems.2015.07.007] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/23/2015] [Revised: 07/13/2015] [Accepted: 07/16/2015] [Indexed: 01/18/2023]
|
9
|
Bérenguier D, Chaouiya C, Monteiro PT, Naldi A, Remy E, Thieffry D, Tichit L. Dynamical modeling and analysis of large cellular regulatory networks. CHAOS (WOODBURY, N.Y.) 2013; 23:025114. [PMID: 23822512 DOI: 10.1063/1.4809783] [Citation(s) in RCA: 36] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/02/2023]
Abstract
The dynamical analysis of large biological regulatory networks requires the development of scalable methods for mathematical modeling. Following the approach initially introduced by Thomas, we formalize the interactions between the components of a network in terms of discrete variables, functions, and parameters. Model simulations result in directed graphs, called state transition graphs. We are particularly interested in reachability properties and asymptotic behaviors, which correspond to terminal strongly connected components (or "attractors") in the state transition graph. A well-known problem is the exponential increase of the size of state transition graphs with the number of network components, in particular when using the biologically realistic asynchronous updating assumption. To address this problem, we have developed several complementary methods enabling the analysis of the behavior of large and complex logical models: (i) the definition of transition priority classes to simplify the dynamics; (ii) a model reduction method preserving essential dynamical properties, (iii) a novel algorithm to compact state transition graphs and directly generate compressed representations, emphasizing relevant transient and asymptotic dynamical properties. The power of an approach combining these different methods is demonstrated by applying them to a recent multilevel logical model for the network controlling CD4+ T helper cell response to antigen presentation and to a dozen cytokines. This model accounts for the differentiation of canonical Th1 and Th2 lymphocytes, as well as of inflammatory Th17 and regulatory T cells, along with many hybrid subtypes. All these methods have been implemented into the software GINsim, which enables the definition, the analysis, and the simulation of logical regulatory graphs.
Collapse
Affiliation(s)
- D Bérenguier
- Institut de Mathématiques de Luminy, Marseille, France
| | | | | | | | | | | | | |
Collapse
|