Search

Digital twins, driven by data and mathematical modelling, have emerged as powerful tools for simulating complex biological systems. In this work, we focus on modelling the clearance on a liver-on-chip as a digital twin that closely mimics the clearance functionality of the human liver. Our approach involves the creation of a compartmental physiological model of the liver using ordinary differential equations (ODEs) to estimate pharmacokinetic (PK) parameters related to on-chip liver clearance. The objectives of this study were twofold: first, to predict human clearance values, and second, to propose a framework for bridging the gap between in vitro findings and their clinical relevance. The methodology integrated quantitative Organ-on-Chip (OoC) and cell-based assay analyses of drug depletion kinetics and is further enhanced by incorporating an OoC-digital twin model to simulate drug depletion kinetics in humans. The in vitro liver clearance for 32 drugs was predicted using a digital-twin model of the liver-on-chip and in vitro to in vivo extrapolation (IVIVE) was assessed using time series PK data. Three ODEs in the model define the drug concentrations in media, interstitium and intracellular compartments based on biological, hardware, and physicochemical information. A key issue in determining liver clearance appears to be the insufficient drug concentration within the intracellular compartment. The digital twin establishes a connection between the hardware chip structure and an advanced mapping of the underlying biology, specifically focusing on the intracellular compartment. Our modelling offers the following benefits: i) better prediction of intrinsic liver clearance of drugs compared to the conventional model and ii)explainability of behaviour based on physiological parameters. Finally, we illustrate the clinical significance of this approach by applying the findings to humans, utilising propranolol as a proof-of-concept example. This study stands out as the biggest cross-organ-on-chip platform investigation to date, systematically analysing and predicting human clearance values using data obtained from various in vitro liver-on-chip systems. Accurate prediction of in vivo clearance from in vitro data is important as inadequate understanding of the clearance of a compound can lead to unexpected and undesirable outcomes in clinical trials, ranging from underdosing to toxicity. Physiologically based pharmacokinetic (PBPK) model estimation of liver clearance is explored. The aim is to develop digital twins capable of determining better predictions of clinical outcomes, ultimately reducing the time, cost, and patient burden associated with drug development. Various hepatic in vitro systems are compared and their effectiveness for predicting human clearance is investigated. The developed tool, DigiLoCs, focuses explicitly on accurately describing complex biological processes within liver-chip systems. ODE-constrained optimisation is applied to estimate the clearance of compounds. DigiLoCs enable differentiation between active biological processes (metabolism) and passive processes (permeability and partitioning) by incorporating detailed information on compound-specific characteristics and hardware-specific data. These findings signify a significant stride towards more accurate and efficient drug development methodologies., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2025 Aravindakshan et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

Prenatal alcohol exposure (PAE) can result in mild to severe consequences for children throughout their lives, with this range of symptoms referred to as Fetal Alcohol Spectrum Disorders (FASD). These consequences are thought to be linked to changes in gene expression and transcriptional programming in the brain, but the identity of those changes, and how they persist into adolescence are unclear. In this study, we isolated RNA from the hippocampus of adolescent rats exposed to ethanol during prenatal development and compared gene expression to controls. Briefly, dams were either given free access to standard chow ad libitum (AD), pair-fed a liquid diet (PF) or were given a liquid diet with ethanol (6.7% ethanol, ET) throughout gestation (gestational day (GD) 0-20). All dams were given control diet ad libitum beginning on GD 20 and throughout parturition and lactation. Hippocampal tissue was collected from adolescent male and female offspring (postnatal day (PD) 35-36). Exposure to ethanol caused widespread downregulation of many genes as compared to control rats. Gene ontology analysis demonstrated that affected pathways included cell adhesion, toxin metabolism, and immune responses. Interestingly, these differences were not strongly affected by sex. Furthermore, these changes were consistent when comparing ethanol-exposed rats to pair-fed controls provided with a liquid diet and those fed ad libitum on a standard chow diet. We conclude from this study that changes in genetic architecture and the resulting neuronal connectivity after prenatal exposure to alcohol continue through adolescent development. Further research into the consequences of specific gene expression changes on neural and behavioral changes will be vital to our understanding of the FASD spectrum of diseases., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2024 Khalifa et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

Kaposi's sarcoma herpesvirus (KSHV) is a leading cause of malignancy in AIDS and current therapies are limited. Like all herpesviruses, KSHV infection can be latent or lytic. KSHV latency-associated nuclear antigen (LANA) is essential for viral genome persistence during latent infection. LANA also maintains latency by antagonizing expression and function of the KSHV lytic switch protein, RTA. Here, we find LANA null KSHV is not capable of lytic replication, indicating a requirement for LANA. While LANA promoted both lytic and latent gene expression in cells partially permissive for lytic infection, it repressed expression in non-permissive cells. Importantly, forced RTA expression in non-permissive cells led to induction of lytic infection and LANA switched to promote, rather than repress, most lytic viral gene expression. When basal viral gene expression levels were high, LANA promoted expression, but repressed expression at low basal levels unless RTA expression was forcibly induced. LANA's effects were broad, but virus gene specific, extending to an engineered, recombinant viral GFP under control of host EF1α promoter, but not to host EF1α. Together, these results demonstrate that, in addition to its essential role in genome maintenance, LANA broadly regulates viral gene expression, and is required for high levels of lytic gene expression during lytic infection. Strategies that target LANA are expected to abolish KSHV infection., Competing Interests: The authors have declared that they have no competing interests., (Copyright: © 2024 Li et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

A new computational framework for FLow cytometric Analysis of Rare Events (FLARE) has been developed specifically for fast and automatic identification of rare cell populations in very large samples generated by platforms like multi-parametric flow cytometry. Using a hierarchical Bayesian model and information-sharing via parallel computation, FLARE rapidly explores the high-dimensional marker-space to detect highly rare populations that are consistent across multiple samples. Further it can focus within specified regions of interest in marker-space to detect subpopulations with desired precision., Competing Interests: The authors have declared that no competing interests exist.

Long noncoding RNAs (lncRNAs) localize in the cell nucleus and influence gene expression through a variety of molecular mechanisms. Chromatin-enriched RNAs (cheRNAs) are a unique class of lncRNAs that are tightly bound to chromatin and putatively function to locally cis-activate gene transcription. CheRNAs can be identified by biochemical fractionation of nuclear RNA followed by RNA sequencing, but until now, a rigorous analytic pipeline for nuclear RNA-seq has been lacking. In this study, we survey four computational strategies for nuclear RNA-seq data analysis and develop a new pipeline, Tuxedo-ch, which outperforms other approaches. Tuxedo-ch assembles a more complete transcriptome and identifies cheRNA with higher accuracy than other approaches. We used Tuxedo-ch to analyze benchmark datasets of K562 cells and further characterize the genomic features of intergenic cheRNA (icheRNA) and their similarity to enhancer RNAs (eRNAs). We quantify the transcriptional correlation of icheRNA and adjacent genes and show that icheRNA is more positively associated with neighboring gene expression than eRNA or cap analysis of gene expression (CAGE) signals. We also explore two novel genomic associations of cheRNA, which indicate that cheRNAs may function to promote or repress gene expression in a context-dependent manner. IcheRNA loci with significant levels of H3K9me3 modifications are associated with active enhancers, consistent with the hypothesis that enhancers are derived from ancient mobile elements. In contrast, antisense cheRNA (as-cheRNA) may play a role in local gene repression, possibly through local RNA:DNA:DNA triple-helix formation., Competing Interests: The authors declare no competing interests.

Proteins are involved in almost all functions in a living cell, and functions of proteins are realized by their tertiary structures. Obtaining a global perspective of the variety and distribution of protein structures lays a foundation for our understanding of the building principle of protein structures. In light of the rapid accumulation of low-resolution structure data from electron tomography and cryo-electron microscopy, here we map and classify three-dimensional (3D) surface shapes of proteins into a similarity space. Surface shapes of proteins were represented with 3D Zernike descriptors, mathematical moment-based invariants, which have previously been demonstrated effective for biomolecular structure similarity search. In addition to single chains of proteins, we have also analyzed the shape space occupied by protein complexes. From the mapping, we have obtained various new insights into the relationship between shapes, main-chain folds, and complex formation. The unique view obtained from shape mapping opens up new ways to understand design principles, functions, and evolution of proteins., Competing Interests: The authors have declared that no competing interests exist.

Motivation: Systems biology faces two key challenges when dealing with large amounts of disparate data produced by different experiments: the integration of results across different experiments, and the extraction of meaningful information from the data produced by these experiments. An ongoing challenge is to provide better tools that can mine data patterns that could not have been discovered through simple visualization. Such mining capabilities also need to be coupled with intuitive visualization to portray those findings. We introduce a software toolbox entitled BioNetApp to mine these patterns and visualize them across all experiments., Results: BioNetApp is an interactive visual data mining software for analyzing high-volume molecular expression data obtained from multiple 'omics experiments. By integrating visualization, statistical methods, and data mining techniques, BioNetApp can perform interactive correlative and comparative analysis along time-course studies of molecular expression data. Correlation analysis provides several visualization features such as Kamada-Kawai, Fruchterman-Reingold Spring embedding network layouts, in addition to single circle, multiple circle and heatmap layouts, whereas comparative analysis presents expression-data distributions across samples, groups, and time points with boxplot display, outlier detection, and data curve fitting. BioNetApp also provides data clustering based on molecular concentrations using Self Organizing Maps (SOM), K-Means, K-Medoids, and Farthest First algorithms., Conclusion: BioNetApp has been utilized in a metabolomics study to investigate the metabolite abundance changes in alcohol induced fatty liver, where pair-wise analyses of metabolome concentration revealed correlation networks and interesting patterns in the metabolomics dataset. This study case demonstrates the effectiveness of the BioNetApp software as an interactive visual analysis tool for molecular expression data in systems biology. The BioNetApp software is freely available under GNU GPL license and can be downloaded (including the case-study data and user-manual) at: https://doi.org/10.5281/zenodo.2563129., Competing Interests: The authors have declared that no competing interests exist.

Protein-protein interactions are the cornerstone of numerous biological processes. Although an increasing number of protein complex structures have been determined using experimental methods, relatively fewer studies have been performed to determine the assembly order of complexes. In addition to the insights into the molecular mechanisms of biological function provided by the structure of a complex, knowing the assembly order is important for understanding the process of complex formation. Assembly order is also practically useful for constructing subcomplexes as a step toward solving the entire complex experimentally, designing artificial protein complexes, and developing drugs that interrupt a critical step in the complex assembly. There are several experimental methods for determining the assembly order of complexes; however, these techniques are resource-intensive. Here, we present a computational method that predicts the assembly order of protein complexes by building the complex structure. The method, named Path-LzerD, uses a multimeric protein docking algorithm that assembles a protein complex structure from individual subunit structures and predicts assembly order by observing the simulated assembly process of the complex. Benchmarked on a dataset of complexes with experimental evidence of assembly order, Path-LZerD was successful in predicting the assembly pathway for the majority of the cases. Moreover, when compared with a simple approach that infers the assembly path from the buried surface area of subunits in the native complex, Path-LZerD has the strong advantage that it can be used for cases where the complex structure is not known. The path prediction accuracy decreased when starting from unbound monomers, particularly for larger complexes of five or more subunits, for which only a part of the assembly path was correctly identified. As the first method of its kind, Path-LZerD opens a new area of computational protein structure modeling and will be an indispensable approach for studying protein complexes.

Disordered protein-protein interactions (PPIs), those involving a folded protein and an intrinsically disordered protein (IDP), are prevalent in the cell, including important signaling and regulatory pathways. IDPs do not adopt a single dominant structure in isolation but often become ordered upon binding. To aid understanding of the molecular mechanisms of disordered PPIs, it is crucial to obtain the tertiary structure of the PPIs. However, experimental methods have difficulty in solving disordered PPIs and existing protein-protein and protein-peptide docking methods are not able to model them. Here we present a novel computational method, IDP-LZerD, which models the conformation of a disordered PPI by considering the biophysical binding mechanism of an IDP to a structured protein, whereby a local segment of the IDP initiates the interaction and subsequently the remaining IDP regions explore and coalesce around the initial binding site. On a dataset of 22 disordered PPIs with IDPs up to 69 amino acids, successful predictions were made for 21 bound and 18 unbound receptors. The successful modeling provides additional support for biophysical principles. Moreover, the new technique significantly expands the capability of protein structure modeling and provides crucial insights into the molecular mechanisms of disordered PPIs.

While RNAs are well known to possess complex structures, functionally similar RNAs often have little sequence similarity. While the exact size and spacing of base-paired regions vary, functionally similar RNAs have pronounced similarity in the arrangement, or topology, of base-paired stems. Furthermore, predicted RNA structures often lack pseudoknots (a crucial aspect of biological activity), and are only partially correct, or incomplete. A topological approach addresses all of these difficulties. In this work we describe each RNA structure as a graph that can be converted to a topological spectrum (RNA fingerprint). The set of subgraphs in an RNA structure, its RNA fingerprint, can be compared with the fingerprints of other RNA structures to identify and correctly classify functionally related RNAs. Topologically similar RNAs can be identified even when a large fraction, up to 30%, of the stems are omitted, indicating that highly accurate structures are not necessary. We investigate the performance of the RNA fingerprint approach on a set of eight highly curated RNA families, with diverse sizes and functions, containing pseudoknots, and with little sequence similarity-an especially difficult test set. In spite of the difficult test set, the RNA fingerprint approach is very successful (ROC AUC > 0.95). Due to the inclusion of pseudoknots, the RNA fingerprint approach both covers a wider range of possible structures than methods based only on secondary structure, and its tolerance for incomplete structures suggests that it can be applied even to predicted structures. Source code is freely available at https://github.rcac.purdue.edu/mgribsko/XIOS_RNA_fingerprint., Competing Interests: The first author, Jiajie Huang, is currently employed by Thermo Fisher Scientific. The second author, Kejie Li, is currently employed by Biogen Idec. These commercial affiliations do not alter the authors' adherence to all PLOS ONE policies on sharing data and materials.