Modelling and analysis of the complement system signalling pathways: roles of C3, C5a and pro-inflammatory cytokines in SARS-CoV-2 infection.

The complement system is an essential part of innate immunity. It is activated by invading pathogens causing inflammation, opsonization, and lysis via complement anaphylatoxins, complement opsonin's and membrane attack complex (MAC), respectively. However, in SARS-CoV-2 infection overactivation of complement system is causing cytokine storm leading to multiple organs damage. In this study, the René Thomas kinetic logic approach was used for the development of biological regulatory network (BRN) to model SARS-CoV-2 mediated complement system signalling pathways. Betweenness centrality analysis in cytoscape was adopted for the selection of the most biologically plausible states in state graph. Among the model results, in strongly connected components (SCCs) pro-inflammatory cytokines (PICyts) oscillatory behaviour between recurrent generation and downregulation was found as the main feature of SARS-CoV-2 infection. Diversion of trajectories from the SCCs leading toward hyper-inflammatory response was found in agreement with in vivo studies that overactive innate immunity response caused PICyts storm during SARS-CoV-2 infection. The complex of negative regulators FI, CR1 and DAF in the inhibition of complement peptide (C5a) and PICyts was found desirable to increase immune responses. In modelling role of MAC and PICyts in lowering of SARS-CoV-2 titre was found coherent with experimental studies. Intervention in upregulation of C5a and PICyts by C3 was found helpful in back-and-forth variation of signalling pattern linked with the levels of PICyts. Moreover, intervention in upregulation of PICyts by C5a was found productive in downregulation of all activating factors in the normal SCCs. However, the computational model predictions require experimental studies to be validated by exploring the activation role of C3 and C5a which could change levels of PICyts at various phases of SARS-CoV-2 infection.

Decoding the Role of Epigenetics in Breast Cancer Using Formal Modeling and Machine-Learning Methods.

Breast carcinogenesis is known to be instigated by genetic and epigenetic modifications impacting multiple cellular signaling cascades, thus making its prevention and treatments a challenging endeavor. However, epigenetic modification, particularly DNA methylation-mediated silencing of key TSGs, is a hallmark of cancer progression. One such tumor suppressor gene (TSG) RUNX3 (Runt-related transcription factor 3) has been a new insight in breast cancer known to be suppressed due to local promoter hypermethylation mediated by DNA methyltransferase 1 (DNMT1). However, the precise mechanism of epigenetic-influenced silencing of the RUNX3 signaling resulting in cancer invasion and metastasis remains inadequately characterized. In this study, a biological regulatory network (BRN) has been designed to model the dynamics of the DNMT1-RUNX3 network augmented by other regulators such as p21, c-myc, and p53. For this purpose, the René Thomas qualitative modeling was applied to compute the unknown parameters and the subsequent trajectories signified important behaviors of the DNMT1-RUNX3 network (i.e., recovery cycle, homeostasis, and bifurcation state). As a result, the biological system was observed to invade cancer metastasis due to persistent activation of oncogene c-myc accompanied by consistent downregulation of TSG RUNX3. Conversely, homeostasis was achieved in the absence of c-myc and activated TSG RUNX3. Furthermore, DNMT1 was endorsed as a potential epigenetic drug target to be subjected to the implementation of machine-learning techniques for the classification of the active and inactive DNMT1 modulators. The best-performing ML model successfully classified the active and least-active DNMT1 inhibitors exhibiting 97% classification accuracy. Collectively, this study reveals the underlined epigenetic events responsible for RUNX3-implicated breast cancer metastasis along with the classification of DNMT1 modulators that can potentially drive the perception of epigenetic-based tumor therapy.

Parameter estimation of qualitative biological regulatory networks on high performance computing hardware.

BACKGROUND: Biological Regulatory Networks (BRNs) are responsible for developmental and maintenance related functions in organisms. These functions are implemented by the dynamics of BRNs and are sensitive to regulations enforced by specific activators and inhibitors. The logical modeling formalism by René Thomas incorporates this sensitivity with a set of logical parameters modulated by available regulators, varying with time. With the increase in complexity of BRNs in terms of number of entities and their interactions, the task of parameters estimation becomes computationally expensive with existing sequential SMBioNET tool. We extend the existing sequential implementation of SMBioNET by using a data decomposition approach using a Java messaging library called MPJ Express. The approach divides the parameters space into different regions and each region is then explored in parallel on High Performance Computing (HPC) hardware. RESULTS: The performance of the parallel approach is evaluated on BRNs of different sizes, and experimental results on multicore and cluster computers showed almost linear speed-up. This parallel code can be executed on a wide range of concurrent hardware including laptops equipped with multicore processors, and specialized distributed memory computer systems. To demonstrate the application of parallel implementation, we selected a case study of Hexosamine Biosynthetic Pathway (HBP) in cancer progression to identify potential therapeutic targets against cancer. A set of logical parameters were computed for HBP model that directs the biological system to a state of recovery. Furthermore, the parameters also suggest a potential therapeutic intervention that restores homeostasis. Additionally, the performance of parallel application was also evaluated on a network (comprising of 23 entities) of Fibroblast Growth Factor Signalling in Drosophila melanogaster. CONCLUSIONS: Qualitative modeling framework is widely used for investigating dynamics of biological regulatory networks. However, computation of model parameters in qualitative modeling is computationally intensive. In this work, we presented results of our Java based parallel implementation that provides almost linear speed-up on both multicore and cluster platforms. The parallel implementation is available at https://psmbionet.github.io .

Formal Modeling of mTOR Associated Biological Regulatory Network Reveals Novel Therapeutic Strategy for the Treatment of Cancer.

Cellular homeostasis is a continuous phenomenon that if compromised can lead to several disorders including cancer. There is a need to understand the dynamics of cellular proliferation to get deeper insights into the prevalence of cancer. Mechanistic Target of Rapamycin (mTOR) is implicated as the central regulator of the metabolic pathway involved in growth whereas its two distinct complexes mTORC1 and mTORC2 perform particular functions in cellular propagation. To date, mTORC1 is a well defined therapeutic target to inhibit uncontrolled cell division, while the role of mTORC2 is not well characterized. Therefore, the current study is designed to understand the signaling dynamics of mTOR and its partner proteins such as PI3K, PTEN, mTORC2, PKB (Akt), mTORC1, and FOXO. For this purpose, a qualitative model of mTOR-associated Biological Regulatory Network (BRN) is constructed to predict its regulatory behaviors which may not be predictable otherwise. The depleted expression of PTEN and FOXO along with the overexpression of PI3K, mTORC2, mTORC1 and Akt is predicted as a stable steady state which is in accordance with their observed expression levels in the progression of various cancers. The qualitative model also predicts the homeostasis of all the entities in the form of qualitative cycles. The significant qualitative (discrete) cycle is identified by analyzing betweenness centralities of the qualitative (discrete) states. This cycle is further refined as a linear hybrid automaton model with the production (activation) and degradation (inhibition) time delays in order to analyze the real-time constraints for its existence. The analysis of the hybrid model provides a formal proof that during homeostasis the inhibition time delay of Akt is less than the inhibition time delay of mTORC2. In conclusion, our observations characterize that in homeostasis Akt is degraded with a faster rate than mTORC2 which suggests that the inhibition of Akt along with the activation of mTORC2 may be a better therapeutic strategy for the treatment of cancer.

Formal modeling and analysis of ER-α associated Biological Regulatory Network in breast cancer.

BACKGROUND: Breast cancer (BC) is one of the leading cause of death among females worldwide. The increasing incidence of BC is due to various genetic and environmental changes which lead to the disruption of cellular signaling network(s). It is a complex disease in which several interlinking signaling cascades play a crucial role in establishing a complex regulatory network. The logical modeling approach of René Thomas has been applied to analyze the behavior of estrogen receptor-alpha (ER-α) associated Biological Regulatory Network (BRN) for a small part of complex events that leads to BC metastasis. METHODS: A discrete model was constructed using the kinetic logic formalism and its set of logical parameters were obtained using the model checking technique implemented in the SMBioNet software which is consistent with biological observations. The discrete model was further enriched with continuous dynamics by converting it into an equivalent Petri Net (PN) to analyze the logical parameters of the involved entities. RESULTS: In-silico based discrete and continuous modeling of ER-α associated signaling network involved in BC provides information about behaviors and gene-gene interaction in detail. The dynamics of discrete model revealed, imperative behaviors represented as cyclic paths and trajectories leading to pathogenic states such as metastasis. Results suggest that the increased expressions of receptors ER-α, IGF-1R and EGFR slow down the activity of tumor suppressor genes (TSGs) such as BRCA1, p53 and Mdm2 which can lead to metastasis. Therefore, IGF-1R and EGFR are considered as important inhibitory targets to control the metastasis in BC. CONCLUSION: The in-silico approaches allow us to increase our understanding of the functional properties of living organisms. It opens new avenues of investigations of multiple inhibitory targets (ER-α, IGF-1R and EGFR) for wet lab experiments as well as provided valuable insights in the treatment of cancers such as BC.

Formal modeling and analysis of the hexosamine biosynthetic pathway: role of O-linked N-acetylglucosamine transferase in oncogenesis and cancer progression.

The alteration of glucose metabolism, through increased uptake of glucose and glutamine addiction, is essential to cancer cell growth and invasion. Increased flux of glucose through the Hexosamine Biosynthetic Pathway (HBP) drives increased cellular O-GlcNAcylation (hyper-O-GlcNAcylation) and contributes to cancer progression by regulating key oncogenes. However, the association between hyper-O-GlcNAcylation and activation of these oncogenes remains poorly characterized. Here, we implement a qualitative modeling framework to analyze the role of the Biological Regulatory Network in HBP activation and its potential effects on key oncogenes. Experimental observations are encoded in a temporal language format and model checking is applied to infer the model parameters and qualitative model construction. Using this model, we discover step-wise genetic alterations that promote cancer development and invasion due to an increase in glycolytic flux, and reveal critical trajectories involved in cancer progression. We compute delay constraints to reveal important associations between the production and degradation rates of proteins. O-linked N-acetylglucosamine transferase (OGT), an enzyme used for addition of O-GlcNAc during O-GlcNAcylation, is identified as a key regulator to promote oncogenesis in a feedback mechanism through the stabilization of c-Myc. Silencing of the OGT and c-Myc loop decreases glycolytic flux and leads to programmed cell death. Results of network analyses also identify a significant cycle that highlights the role of p53-Mdm2 circuit oscillations in cancer recovery and homeostasis. Together, our findings suggest that the OGT and c-Myc feedback loop is critical in tumor progression, and targeting these mediators may provide a mechanism-based therapeutic approach to regulate hyper-O-GlcNAcylation in human cancer.