| |
Metabolic remodeling of microbes by obligate intracellular parasites alters species evenness in mutualistic communities
Ave Bisesi and William Harcombe
University of Minnesota, Minneapolis, Minnesota, United States
Bacteria carry many types of obligate intracellular parasites, including plasmids and bacteriophage (viruses that infect bacteria). During infection, these parasites redirect intracellular resources away from bacterial processes toward parasite production. Because parasite-induced metabolic changes influence host phenotypes, parasitic infection is predicted to alter how microbes contribute to important community and ecosystem functions. Yet little is known about how infection shapes interactions between host and non-host species. Here we integrate a genome-scale metabolic modeling approach with an in vitro obligate cross-feeding system to investigate the metabolic consequences of two intracellular parasites of Escherichia coli – the F-plasmid and the filamentous phage M13 – on interactions between bacteria in a multispecies community composed of E. coli, Salmonella enterica, and Methylorubrum extorquens. Modeling predicts that metabolic conflict due to parasitic infection of E. coli changes the types and quantities of compounds secreted by infected hosts, increasing species evenness when bacterial species are engaged in obligate mutualism. These theoretical results are corroborated by in vitro experiments. Our work emphasizes that microbes infected by diverse obligate intracellular parasites have fundamentally different metabolisms from uninfected cells and demonstrates that these metabolic shifts can have significant consequences for microbial community structure and function.
|
| |
Mathematical Modeling of Acute Inflammatory Response to determine onset of Sepsis and Septic shock in ICU patients
Ayush Ranjan and Venkatesh Kareenhalli
University of Pittsburgh, Pittsburgh, Pennsylvania, United States
The immune system works as a well-oiled machine during a pathogen attack. The dysfunction of this well-oiled machine can be catastrophic, such as sepsis, with death being one of the ultimate consequences. Although several studies have shown the progression of sepsis and the complex dynamic interaction between the immune and metabolism playing a role, the prediction of the onset of sepsis and sepsis shock is still challenging. The intellectual construct of sepsis is based on this immune dysfunction either in its failure to clear out the foreign invaders attacking the body or the breakdown of the counter-attacking immune system due to internal failures of their mechanisms. The current work here proposes an ODE-based model that accounts for immune, metabolic, and physiological variables to address this issue of the onset of sepsis. The model developed was used to analyze system-level dynamics and predict the healthy and sepsis response. The interplay between pro-inflammatory and anti-inflammatory cytokines and Nitric oxide and lactate levels form the basis of sensitivity toward sepsis. Sensitivity Analysis gives an understanding of the different stages of the infection and the various parameter regimes that progress to sepsis and shock. The model makes use of states projected along with parameter regimes for the ICU patients with the model tweaked for the ICU patients to match the onset of sepsis and make predictions for its progression thus creating a model-based control system for the prediction of the onset of sepsis and septic shock in ICU patients. |
| |
MetNetComp: Database for minimal and maximal gene-deletion strategies for growth-coupled production of genome-scale metabolic networks
Takeyuki Tamura
Kyoto University, Kyōto, Japan
Cancer metastasis accounts for nearly 90% of cancer-related deaths. Epithelial to Mesenchymal Transition (EMT) of cancer cells is key to cancer metastasis. Recent preclinical and clinical data suggest that hybrid E/M cell states harbor the highest metastatic fitness instead of the fully mesenchymal cells. However, the mechanistic details of their survival strategies during metastasis remain unclear. We model the dynamics of a minimalistic gene regulatory network (GRN) between regulators of EMT and PD-L1, a known immune suppressor. We show that hybrid E/M states are highly likely to exhibit high PD-L1 levels, like those seen in mesenchymal cells, thus obviating the need to undergo a complete EMT to cause immunosuppression. We show that the switch from an epithelial/low-PDL1 state to a hybrid/mesenchymal state with high-PDL1 is reversible. We demonstrate that acquiring phenotypic resistance to targeted therapy can co-occur with high PD-L1 levels, enabling cross-resistance and enhancing breast cancer cell fitness during metastasis. We validate our model predictions by extensive analysis of transcriptomic datasets across multiple cancers at bulk and single-cell levels. Our results highlight how the emergent dynamics of interconnected GRNs can coordinate various axes of cellular fitness during metastasis, thus laying the foundation for the rational design of clinical therapies. |
| |
Unraveling the Dual Inhibitory Mechanism of Compound 22ac: A Molecular Dynamics Investigation into ERK1 and ERK5 In-hibition in Cancer
Elliasu Salifu, Pritika Ramharack, Mbuso A. Faya and James Abugri
South African Medical Research Council, Tygerberg, South Africa
Cancer remains a major challenge in the field of medicine, necessitating innovative therapeutic strategies. Mitogen-activated protein kinase (MAPK) signaling pathways, particularly Extracellular Signal-Regulated Kinase 1 and 2 (ERK1/2), play pivotal roles in cancer pathogenesis. Recently, ERK5 (also known as MAPK7) has emerged as an attractive target due to its compensatory role in cancer progression upon termination of ERK1 signaling. This study explores the potential of Compound 22ac, a novel small molecule inhibitor, to simultaneously target both ERK1 and ERK5 in cancer cells. Using molecular dynamics simulations, we investigate the binding affinity, conformational dynamics, and stability of Compound 22ac when interacting with ERK1 and ERK5. Our results indicate that compound 22ac forms strong interactions with key residues in the ATP-binding pocket of both ERK1 and ERK5, effectively inhibiting their catalytic activity. Further, the simulations reveal subtle differences in the binding modes of Compound 22ac within the two kinases, shedding light on the dual inhibitory mechanism. This research not only elucidates a structural mechanism of action of Compound 22ac, but also highlights its potential as a promising therapeutic agent for cancer treatment. The dual inhibition of ERK1 and ERK5 by Compound 22ac offers a novel approach to disrupting the MAPK signaling cascade, thereby hindering cancer progression. These findings may contribute to the development of targeted therapies that could improve the prognosis for cancer patients.
|
| |
Atomic elementary flux modes explain the steady state flow of metabolites in genome-scale flux networks
Justin G. Chitpin and Theodore J. Perkins
University of Ottawa, Ottawa, Canada
Elementary flux modes (EFMs) are minimal sets of reactions that can sustain steady state flux in a metabolic flux network. First proposed by Schuster and Hilgetag, any set of steady state fluxes in a metabolic network can be explained as a positive, linear combination of their EFM weights. This decomposition, however, is generally not unique, as there are typically exponentially many EFMs. We recently addressed this flux decomposition problem for the special case of unimolecular reaction networks involving single substrates and products. By imposing a Markovian constraint on EFMs, we modelled steady state particle fluxes as a cycle-history Markov chain (CHMC) to uniquely identify EFM weights that reconstruct any network fluxes. Here, we generalize our CHMC analysis to networks containing multispecies reactions. We propose to concept of atomic EFMs which are steady state pathways that trace the flow of individual atoms through a network. Using a state-of-the-art atom mapping algorithm, we enumerate carbon and nitrogen EFMs from four genome-scale networks and compute atomic EFM weights from a fifth HepG2 dataset. We find that enumerating atomic EFMs is computationally tractable compared to standard EFMs, and show how atomic EFMs characterize nutrient source metabolism through the network. Lastly, our analysis of the HepG2 atomic EFM weights reveals that the majority of steady state fluxes are explained by a small fraction of atomic EFMs. Our atomic CHMC method is a fast and novel tool to quantify metabolic remodelling as technologies to generate large-scale metabolic flux networks advance.
|
| |
Mitovolve Models the Evolution of the Prevalence of Mitochondrial Mutations in Pediatric Leukemia
Yonghui Ni, Melissa Franco, Kelly McCastlain, Ti-Cheng Chang, Gang Wu, Catherine Welsh, Knostantin Khrapko, Mondira Kundu and Stanley Pounds
St. Jude Children’s Research Hospital, Memphis, Tennessee, United States
We propose Mitovolve as a statistical model that utilizes single-cell sequencing data to infer the evolutionary history of somatic mitochondrial DNA (mtDNA) mutations. Mitovolve models the evolution of a mtDNA mutation from one cell of origin with a specified number of mutant and wild-type mitochondrial genomes to its observed abundance across the numerous cells represented in single-cell sequencing datasets. The distribution of the number of mutant mitochondrial genomes in each daughter cell is derived by hypergeometric sampling from a pool of duplicated mitochondrial genomes in the parent cell, a process that is reiterated for each successive generation of new cells. Given the initial abundance of mutant mitochondrial genomes, the model calculates an exact probability distribution for the number of mutated mitochondrial genomes per cell after a specified number of replicative generations in the presence or absence of selective pressure parameters represented as a cubic polynomial function of the mutant allele frequency (MAF). Maximum likelihood estimation is used to fit many models (varying the initial mutant allele fraction and number of replicative generations) with and without selection to the observed data. A likelihood ratio test is then used to compute confidence sets and p-values for the presence and nature of selective pressure operating on the mtDNA mutation. As proof of concept, we applied Mitovolve to single cell sequencing data from a primary pediatric leukemia sample and found statistically significant evidence in favor of selective pressure operating on a tumor-enriched mtDNA mutation. Mitovolve will soon be available as an R package.
|
| |
Python Hurdle Model for Differential Expression in Single-Cell RNA-seq
Negin Rahimzadeh and Vivek Swarup
University of California, California, United States
Background and Motivation: Single-cell RNA sequencing (scRNA-seq) has significantly enhanced the resolution and scale at which transcriptomics can be studied. Despite these advancements, it has also led to an expansion in the volume of data generated by single-cell studies, further complicated by the bimodal distribution of gene expression and ‘dropout’ events. To address these intricacies, MAST, a versatile statistical framework in R, employs a hurdle model—a dual-component generalized linear model adept at characterizing the discrete occurrence of gene expression from the continuous spectrum of expression levels. While R continues to be an invaluable tool for detailed statistical analysis, Python’s strengths in handling large datasets, scalability, and data integration techniques make it a suitable platform for projects involving large-scale genomic data.
Method and Results: We propose developing a Python-based hurdle model that integrates logistic regression for zero inflation and a negative binomial model for expression quantification. Adjustments include a refined approach to managing zeros and leveraging Python’s extensive libraries for optimization and model fitting. To demonstrate the efficacy of our framework, we validated it by analyzing the Peripheral Blood Mononuclear Cells dataset. Preliminary comparisons of our results with those obtained from the MAST indicate a strong positive Spearman’s correlation (ρ = 0.731) with a highly significant p-value (1.43e-257), showcasing the potential of our approach to offer comparable analytical depth. Further optimizations are underway to enhance model accuracy and data handling.
Conclusion: Our Python-based model offers a powerful approach for identifying differentially expressed genes across conditions in large scRNA-seq datasets. |
| |
Using constraint-based modeling to identify metabolic signatures associated with recurrence in ductal carcinoma in situ
Nirvana Nursimulu and Sushant Kumar
Princess Margaret Cancer Centre, Toronto, Canada
At least 50,000 women in the United States are expected to be diagnosed this year with ductal carcinoma in situ (DCIS), a type of breast precancer. To decrease the risk of DCIS progression to invasive tumours, breast-conserving lumpectomies or mastectomies are typically performed. However, patient treatment and follow-up strategies remain to be improved given significant DCIS recurrence rate despite surgery as well as growing concern that such aggressive interventions may be unnecessary in some cases. To complement the search for genetic and epigenetic markers of recurrence, we employ constraint-based metabolic modeling to identify metabolic signatures of DCIS recurrence as a means to improve patient outcome. We leverage bulk RNA-sequencing data from the Translational Breast Cancer Research Consortium to identify expression patterns at initial diagnosis distinguishing patients who subsequently experienced recurrence within the next 5 years (n=63) from those with no recurrence (n=95). We mapped these expression differences to the most recent human genome-scale metabolic reconstruction (Human1) and identified cellular metabolic tasks that were significantly up and downregulated in patient cohorts (for instance, fatty acid synthesis upregulation and phospholipid synthesis downregulation). We further applied deconvolution techniques to identify metabolic rewiring at cell type resolution. Through our cell type-specific map of metabolism in DCIS, we provide a more comprehensive view of the interaction between precancerous cells and the tumour microenvironment which may, in turn, provide a unique perspective into recurrence and the transition from DCIS to tumour invasion. |
| |
Deciphering the Pathogenic Mechanisms of LPMOs using Systems Biology
Pavinap Priyaa, Mezya Sezen and Ragothaman Yennamalli
SASTRA Deemed to be University, India
Lytic polysaccharide monooxygenases (LPMOs) are recently discovered oxidoreductases that can oxidatively cleave the glycosidic bonds found in carbohydrate polymers enabling the depolymerisation of recalcitrant biomasses. While most have explored the role of LPMOs in a (plant) biomass conversion/biofuel context, alternative roles have emerged such as in pathogenesis of communicable diseases. Various experimental studies on LPMOs from Vibrio cholerae, Paenibacillus larvae, Listeria monocytogenes, Pseudomonas aeruginosa and Lysobacter enzymogenes have shown involvement in pathogenesis and virulence roles. However, the molecular mechanism of pathogenesis remains unclear. Here, through systems biology and network analysis, we aimed to decipher the interplay between LPMOs and other bacterial proteins, utilizing data from the STRING database to construct interaction maps. Using Gene Ontology (GO) enrichment analysis, protein-protein docking, network perturbation analysis, we have attempted to decipher the pathogenic role of LPMOs. We focused on V. cholerae LPMOs, namely gbpA and VC_A0140. GO enrichment analysis of gbpA revealed extracellular localization, metalloendopeptidase activity, metal ion binding as molecular functions , and proteolysis as biological processes. The protein-protein docking was done using HADDOCK for two pairs, gbpA with hap and gbpA with VC_1952 and analyzed using PDBsum, where the presence of salt bridges and hydrogen-bonded interactions indicate strong interactions. Similarly, VC_A0140 showed that it is extracellularly located, ATP binding, Metal ion binding, and metalloendopeptidase for molecular function, and proteolysis for biological processes. This study can help understand LPMOs’ role in pathogenesis, for future computational studies, and identify targets for designing inhibitors to prevent bacterial pathogens from invading host cells.
|
| |
Fusion of Spatiotemporal and Network Models to Quantify Variations in ScRNA-seq
Osafu Augustine Egbon and Benedict Anchang
National Institute of Environmental Health Sciences, Durham, North Carolina, United States
A major challenge with analyzing spatially resolved single-cell data collected at multiple time points is temporal misalignment. That is, spatial coordinates captured at different time snapshots may represent different tissue regions. Despite technological efforts to address this issue, the inherent variability among different subjects hampers complete success. Analyzing such data using existing pipelines may lead to problems, inefficiencies, and some degree of inaccuracy. To mitigate these challenges, we propose SpatialDNet, which fuses network and spatio-temporal models to quantify cellular behavior within a cellular microenvironment. SpatialDNet takes both spatial location information of cells or spots and genomic data as input. It constructs a rectangular mesh over the spatial region, creates network nodes from these rectangles, and utilizes multidimensional gene expression data to connect the nodes, forming a comprehensive network. We validated the analysis framework using a CODEX dataset obtained from the study of advanced melanoma tumors to estimate spatio-temporal variable genes, determine gene regulatory networks, and annotate the trajectory of the tumor-immune cell network. In the experimental protocol, CD8+ T cells activated ex-vivo with gp100 antigen and IL-2 were transferred into mice with B16-F10 tumors. Tumors were harvested and imaged at day 0, 1, 3, 5, and 12 post-treatment. The results of the analysis revealed distinct dynamic tumor-immune relationships following T-cell therapy. Specifically, there was a significant increase in tumor and immune cell proliferation, which continued until day 5 but declined by day 12. These findings could provide valuable insights into the efficacy of T-cell therapy among patients.
|
| |
A computational platform to simulate cell state transitions in single-cell RNA-seq data
Dan Peng and Patrick Cahan
Johns Hopkins University, Baltimore, Maryland, United States
Computational modelling of cell state transitions has been a great interest of many researchers in the field of developmental biology, cancer biology and cell fate engineering. It allows us to study the biological systems and perform in-silico perturbations that would otherwise be expensive to do experimentally. Recent advancement in single-cell RNA sequencing (scRNA-seq) technology allows us to capture high-resolution snapshots of cells transitioning states across trajectories. Using these high-throughput datasets, we can train computational models that generate in-silico synthetic cells to faithfully mimic the developmental processes. Here we present OneSC, a platform that can simulate synthetic cells across developmental trajectories using systems of stochastic differential equations govern by a core transcription factors (TFs) regulatory network. Different from the current network inference methods, OneSC prioritizes on generating Boolean network that produces faithful cell state transitions and steady cell states that mimic real biological systems. Applying OneSC to real data, we inferred a core TF network using a mouse myeloid progenitor scRNA-seq dataset and showed that the dynamical simulations of that network generate synthetic single-cell expression profiles that faithfully recapitulate the four myeloid differentiation trajectories going into differentiated cell states (erythrocytes, megakaryocytes, granulocytes and monocytes). Finally, through the in-silico perturbations of the mouse myeloid progenitor core network, we showed that OneSC can accurately predict cell fate decision biases of TF perturbations that closely match with previous experimental observations. In summary, OneSC is a computational platform that simulates synthetic cells to mimic cell state transitions observed in real single-cell expression data.
|
| |
Dynamic modeling of hepatitis B virus infected hepatocytes using peripheral viral markers and validation with liver biopsy long-read DNA sequencing
Cole Campton, David Pan, Mario Cortese, Jeff Wallin, Xiaobai Sun, Chi-yi Chen, Tai-Chung Tseng, Scott Balsitis, Yao-Chun Hsu and Bryan Downie
Duke University, Durham, North Carolina, United States
Background: Functional cure in chronic hepatitis B (CHB) is defined as persistent loss of detectable serum hepatitis B surface antigen (HBsAg) and HBV DNA. Integrated HBV DNA (intDNA) in hepatocytes is a major source of HBsAg, but current methods to quantify intDNA require invasive liver biopsies.
Methods: We developed a mathematical model to estimate the relative abundance of hepatocytes with intDNA or covalently closed circular DNA (cccDNA) using serum viral measurements. This non-linear mixed-effects model was fit using stochastic approximate expectation maximum (SAEM) to data from CHB patients (n=119) in a randomized, double-blind, placebo-controlled study of Tenofovir Disoproxil Fumarate
(TDF) 300 mg/day or placebo over 3 years (NCT01522625). Paired core liver biopsies (n=42) collected per protocol at baseline and year 3 were used to determine observational relative cccDNA and intDNA content and for a novel comparison to modeled cccDNA and intDNA.
Results: Fitted model results provide individual-specific estimates of viral protein production rates and infection cell quantities at each collection time. There are significant correlations between modeled and observed cccDNA cell population (Pearson’s r=.40, p<.001) as well as intDNA cell population (Pearson’s r=.44, p<.001). Resulting viral protein rates indicate most sub-viral particles are produced by intDNA, in-line with biological expectation for this HBeAg- population.
Conclusions: Model results reflect biological expectation and can reproduce liver biopsy-derived results from only serum measurement. These results provide evidence for the utility of modeling of HBV infection dynamics in understanding cccDNA and intDNA disease burden through non-invasive means for HBeAg- CHB patients.
|
| |
Evaluating Transcriptomic Integration for Cyanobacterial Metabolic Modelling
Thomas Pugsley, Christopher Duffy and Guy Hanke
Queen Mary University of London, London, United Kingdom
Transcriptomic integration with genome-scale metabolic models is crucial for deepening our understanding of complex systems in biology. While existing validation studies for transcriptomic integration offer some insights, their scope is limited, particularly for organisms like cyanobacteria, for which little metabolic flux data are available. The impact of pre-processing decisions on integration has scarcely been assessed beyond human models, with no thorough exploration of parameter choices in valve-based integration methods. This study evaluates these methodological decisions using the genome-scale model for Synechocystis sp. PCC 6803 (iSynCJ816 [1]) with existing transcriptomic data to explore their structural impact on the metabolic network and their effect in biomass-optimised scenarios. Preliminary analyses indicate that selecting an appropriate integration method requires consideration for the intended application, highlighting the need for tailored strategies that capture the desired biological insight. Upon completion of this study, we expect to establish foundational guidelines for accurately modelling the cellular phenotype, suitable for downstream analysis and broadly applicable across systems biology.
[1] Joshi, C.J., Peebles, C.A., & Prasad, A. (2017). Modeling and analysis of flux distribution and bioproduct formation in Synechocystis sp. PCC 6803 using a new genome-scale metabolic reconstruction. Algal Research-Biomass Biofuels and Bioproducts, 27, 295-310.
|
| |
Docking simulation studies of potent BCR-ABL inhibitors containing the T315I gatekeeper resistance mutation
Jung Woo Park and Junehawk Lee
Korea Institute of Science and Technology Information, Korea
Modeling the interaction between the target protein and the inhibitor and analyzing the binding energy to understand the structure-activity relationship of the inhibitor is very helpful in new drug design research. In particular, recent improvements in computing power and advances in analysis equipment technology such as X-ray crystallography have made accurate structure prediction and rapid calculation possible. Accordingly, new drug development using computers is receiving attention.
Design of kinase inhibitors targeting the BCR-ABL oncoprotein is an effective strategy for the treatment of chronic myeloid leukemia (CML). In this study, we conducted a study to develop a pan-BCR-ABL inhibitor containing the imatinib resistance mutation T315I. To explain the experimentally observed Abl kinase inhibition results, molecular docking of indazole-based derivatives against the catalytic kinase domain of AblWT (PDB code: 3OXZ) and AblT315I (PDB code: 3OY3) was performed from a three-dimensional structural perspective. In particular, structural differences in the candidate substances that affect the Bcr-Abl kinase inhibitory effect were revealed through analysis of binding method, binding energy, and ligand-protein binding.
|
| |
DCS Clustering and the Dopamine-Angiotensin KEGG Contingency
Brian Westwood, Liliya Yamaleyeva and Rong Chen
Wake Forest School of Medicine, North Carolina, United States
Motivation: Swiss-Prot keyword searches provide a wealth of biocurated proteins with which to gauge a topic. We provide a paradigm for directed investigation of KEGG pathway analysis, beginning with discrete correlate summation (DCS) clustering of Swiss-Prot sets. Dopamine/dopaminergic and angiotensin/renin related pole proteins assemble in pathway analyses upon a field of proteins also associated with 16 other designated keywords to provide a juxtaposed trial of their biological indices. |
| |
Lightweight k-mer model generation for nanopore sequencing
Hiruna Samarakoon, Yuk Kei Wan, Sri Parameswaran, Jonathan Göke, Hasindu Gamaarachchi and Ira W. Deveson
University of New South Wales, Sydney, Australia
Our study focuses on generating a lightweight k-mer model tailored for RNA004 chemistry in nanopore sequencing, addressing challenges in accurate alignment and signal simulation, particularly in the absence of optimized official models. Leveraging the move table and k-mer sequences, our approach ensures precise event alignment from raw signals to basecalled sequences. We showcase the effectiveness of our custom k-mer model through high alignment rates (97.48%) and improved signal simulation accuracy compared to default models. Additionally, our 5-mer model exhibits similar performance as the default 9-mer models in methylation detection in m6A site detection. This work underscores the significance of lightweight k-mer models for computational efficiency and reliable signal simulation in nanopore sequencing applications.
|
| |
Fibonacci’s Blueprint in Lifespan: Integrating Mathematical Predictions with RNA Dynamics and DNA Heritability in Systems Biology
Rob Sacco
Fibonacci, LifeChart
Lifespan, characterized by significant interspecies variation, is shaped by an intricate interplay of genetic, environmental, and physiological factors. This study, anchored in Systems Biology, explores the influence of entropy, fractals, and Fibonacci energy optimization on lifespan evolution. Using the AnAge database, we elucidated the patterns linking the complex mechanics of aging with the Fibonacci sequence. Notably, our findings reveal a significant association between maximum lifespan and Fibonacci age, fostering an intriguing discourse on their connectedness. This connection sparks additional inquiries into the mechanisms and evolutionary significance of Fibonacci energy optimization in lifespan regulation. By weaving RNA structomics and DNA heritability considerations into a Systems Biology framework, we aim to unravel the sophisticated interactions bridging genotypes and phenotypes, offering fresh insights into the principles governing lifespan variation. Our study, synergizing Systems Biology principles with empirical examination, provides valuable insights into the multifaceted nature of lifespan and sets the foundation for future investigations in this fascinating field.
|
| |
Integrative Machine Learning Frameworks for Clinical Outcomes from multimodal data
Sonika Tyagi
School of Computing Technologies, RMIT University Melbourne Australia
The intersection of genomic information and healthcare data presents a promising avenue for personalized medicine and improved patient outcomes. However, the integration of these diverse data sources poses significant challenges. We explore the imperative need for integrative frameworks that effectively harmonize genomic data with clinical records and other healthcare information. By bridging the gap between genomics and healthcare, such frameworks hold the potential to unlock insights into disease mechanisms, optimize treatment strategies, and enhance preventive care. Key considerations in developing these frameworks include interoperability, data privacy, standardization, and scalability. In this work we will present and discuss emerging methodologies and technologies aimed at addressing these challenges, including data integration techniques, machine learning algorithms, and secure data sharing platforms.
|
| |
Cell-cycle dependent DNA repair and replication unifies patterns of chromosome instability
Bingxin Lu, Samuel Winnall, William Cross and Chris Barnes
University of Surrey, United Kingdom
|
| |
BMDx 2: A Tool for Integrating Toxicogenomics-Based Dose-Dependency Analysis and AOP-Based Mechanistic Insights
Angela Serra, Michele Fratello, Giorgia Migliaccio, Laura Saarimäki, Giusy del Giudice, Alisa Pavel, Jack Morikka and Dario Greco
Tampere University, Tampere, Finnland
Abstract. Transcriptomics-based benchmark dose (BMD) analysis is a promising approach for chemical risk assessment and toxicity testing. It provides mechanistic insights into chemical exposures and associated alerts. To enhance the implementation of toxicogenomics-based BMD analysis in the regulatory framework, the use of mechanistic evidence, such as the adverse outcome pathway (AOP) concept, is promoted.
In this work, we present BMDx 2, a new version of the BMDx tool for dose-dependent analysis of toxicogenomic data, which is currently under development. BMDx 2 allows the user to perform an AOP-based characterization of the mechanism of action of the tested chemicals by linking the dose-dependent genes with possible molecular initiating events and adverse outcomes.
The software efficiently handles complex BMD analysis of toxicogenomic data, including microarray and RNA-Seq data. BMDx 2 includes an expanded model list and enables consensus averaging of multiple models for each gene. Moreover, it implements a re-engineered code for easier extension to other models. With a comprehensive pipeline from data filtering to visualizing transcriptomic data and AOP-level BMD results, BMDx 2 enables efficient and robust analysis. Additionally, the tool offers visualization of AOPs categorized by various human health safety indicators, further enhancing its utility in hazard assessment.
|
| |
Cell-cell interaction analysis using a novel method for estimating ligand diffusion distances based on spatial omics data
Haruka Hirose, Yasuhiro Kojima and Teppei Shimamura
Department of Computational and Systems Biology, Medical Research Institute, Tokyo Medical and Dental University
Spatial omics techniques have made significant advancements in recent years. These techniques have enhanced our comprehension of the spatial localization of cells and cell-cell interactions. Cell-cell interactions are essential for maintaining tissue homeostasis. Dysregulation of these interactions can cause a variety of diseases, underscoring the importance of understanding their mechanisms. Ligand proteins, one of the signaling molecules responsible for cell-cell interactions, bind to cell surface receptors and trigger intracellular signaling pathways. The effective range of ligand signaling may play an important role in determining the occurrence of cell-cell interactions, but has not been systematically studied. In this presentation, we developed a new model of ligand spatial diffusion from spatial transcriptome data to estimate the spatial extent of ligand action. Specifically, using Visium data, we learned the diffusion distance as a parameter by regressing the expression of downstream target genes against the diffused ligand expression. This model succeeded in distinguishing between ligands with long and short effective distances, allowing analysis of cell-cell interactions taking into account the effective distance of the ligand. In this study, we present the results of our analysis using colorectal cancer specimens.
|
| |
Biocompilation assembles languages to express biological connectivity
Adriana Climescu-Haulica
Avantiv Limited Cambridge, United Kingdom
Next generation computational biology is needed to assist the emergent field of molecular system medicine. The development speed of omics tools fuels the increase of multi-omics knowledge databases, moving the medicine towards molecular grounds. The question is the integration of omics with physiology, for which new conceptual methods are in demand. Our work comes from the following observation: the bio-molecular relationships are extremely complex while the investigating methods, such as bio-networks and features extraction from bio-data sets, are too restrictive to model living organism complexity. Important dynamical information is lost between. The knowledge acquired until now is appealing to infer, using advanced mathematical objects, more biological insights, so better medical logic. With this vision we introduce a pioneering approach which address a central concept: the biological connectivity. Indeed, the network model, an engineering object, is too schematic to
describe connection in living system sense. We create a bio-compiler, as a comprehensive and flexible model, which main characteristic is ‘modus
communicatio’ as ‘modus operandi.’ Based on regular grammars, defined by means of bio-algebras having theorems/rules biological inputs, the
architecture of this compiler is created so that there is a unique minimal automaton accepting it. This construction is not only rigorous but it
produces a decidable state complexity. As an example of biocompilation architecture, we show how to asses breakpoints of the physiological
dynamic normalization between endocrino-immuno-lymphatic system, using communication between Endogrammar, Immunogrammar and
Lymphogrammar, associated language entities. Biocompilation opens rigorous and comprehensive perspectives on healthcare future, featuring
bio-connectivity as language communication.
|
| |
Multiscale modeling of oligomer assembly
Federico Fontana, Mustafa Ozen and Carlos Lopez
Altos Labs, Inc.
The self-assembly of oligomers is essential for cellular differentiation and response to various stimuli. These processes are not random but result from intricate interplay within extensive molecular interaction networks. On the nanoscale, techniques like molecular dynamics (MD) simulations, including steered MD simulations and umbrella sampling approaches, can provide valuable insights into oligomer formation mechanisms. These techniques offer profound insights into the thermodynamics and kinetics of oligomer formation, providing a detailed view of the molecular forces and structural changes involved. On the microscale, studying the dynamics of Chemical Reaction Networks (CRNs) can help understand and predict oligomers’ holistic role in the behavior of more complex chemical systems such as large-scale protein-protein interaction networks. Furthermore, despite the challenges, studying the causal effects of these multiscale processes in harmony can offer a more comprehensive insight into protein assembly and interaction, contributing significantly to our understanding of cellular behavior and its regulation by molecular interactions. In this research, we present a multiscale modeling framework allowing the integrated study of nano and microscale biology of oligomers, exemplified using the PP1-GADD34-eIF2B complex that is critical in the termination of cellular stress response. We exploited Umbrella Sampling approaches to estimate the dissociation constants for each monomer and get structural insights into the PP1 holoenzyme complex assembly. Then, we leveraged CRN dynamics to assess the combinatorial complexity leading to the multimer organization and relating it to different cellular states.
|
| |
BioRECIPE (Biological system Representation for Evaluation, Curation, Interoperability, Preserving, and Execution)
Gaoxiang Zhou and Natasa Miskov-Zivanov
University of Pittsburgh, Pennsylvania, United States
Many representation formats for computational modeling, such as SBML and CellML, rely on markup languages and prioritize computer-friendliness, but often lack accessibility for human modelers in tasks like model curation, creation, verification, evaluation and expansion. There is a growing demand for a format that enables human users to preview and modify various aspects of models, including individual elements, interactions between elements, and the attributes associated with these interactions. In this work, we introduce the BioRECIPE (Biological system Representation for Evaluation, Curation, Interoperability, Preserving, and Execution) format to address this demand. BioRECIPE is a tabular format ideal for models structured as directed graphs, maintaining machine-readability and compatibility with various model development and analysis tools. BioRECIPE supports two distinct formats: event-based lists of interactions and element-based executive models. In the event-based format, each row corresponds to a single interaction, with column headers aligning with interaction attribute names, categorizing attributes into elements, interactions, context, and provenance. On the other hand, the element-based format assigns each row to a model element and supports a model’s static graph structure along with attributes necessary for dynamic analysis. This format also accommodates various schemes, ranging from less to more detailed, including static and dynamic parameters such as rate and timing. The BioRECIPE format is compatible with established formats like SBML and JSON, facilitating seamless translation between various output formats obtained from INDRA. It also serves as a standardized format within DySE framework and is integrated into many analysis tools, including FLUTE, VIOLIN, CLARINET, ACCORDION, FIDDLE, and PIANO.
|
| |
Biophysical modeling of anisotropic brain tumor growth
Mutaz Dwairy and Junuthula Reddy
Department of Civil Engineering, Yarmouk University, Irbid, Jordan
Solid tumors have high interstitial fluid pressure (IFP), high mechanical stress, and low oxygen levels. Solid stresses may induce apoptosis, stimulate the invasiveness and metastasis of cancer cells and lower their proliferation rate, while oxygen concentration may affect the response of cancer cells to treatment. Although tumors grow in a nonhomogeneous environment, many existing theoretical models assume homogeneous growth and tissue has uniform mechanical properties. For example, the brain consists of three primary materials: white matter, gray matter, and Cerebrospinal fluid (CSF). Therefore, tissue inhomogeneity should be considered in the analysis.
This study established a physical model based on convection-diffusion equations and continuum mechanics principles. The model considers the geometrical inhomogeneity of the brain by including the three different matters in the analysis: white matter, gray matter, and CSF. The model also considers fluid-solid interaction and explicitly describes the effect of mechanical factors, e.g., solid stresses and IFP, chemical factors, e.g., oxygen concentration, and biological factors, e.g., cancer cell concentration, on growing tumors.
In this article, we applied the model on a brain tumor positioned within the white matter, considering the brain inhomogeneity to estimate solid stresses, IFP, the cancer cell concentration, oxygen concentration, and the deformation of the tissues within the neoplasm and the surrounding. Tumor size was estimated at different time points. This model might be clinically crucial for cancer detection and treatment planning by measuring mechanical stresses, IFP, and oxygen levels in the tissue.
|
in Montreal, Quebec, Canada. The meeting will feature two keynote talks and contributed presentations.
