Artificial targeting of the NEIL1 DNA glycosylase to the mitochondria.

The human NEIL1 DNA glycosylase is one of 11 mammalian glycosylases that initiate base excision repair. While substrate preference, catalytic mechanism, and structural information of NEIL1's ordered residues are available, limited information on its subcellular localization, compounded by relatively low endogenous expression levels, have impeded our understanding of NEIL1. Here, we employed a previously developed computational framework to optimize the mitochondrial localization signal of NEIL1, enabling the visualization of its specific targeting to the mitochondrion via confocal microscopy. While we observed clear mitochondrial localization and increased glycosylase/lyase activity in mitochondrial extracts from low-moderate NEIL1 expression, high NEIL1 mitochondrial expression levels proved harmful, potentially leading to cell death.

ImmCellTyper facilitates systematic mass cytometry data analysis for deep immune profiling.

Mass cytometry is a cutting-edge high-dimensional technology for profiling marker expression at the single-cell level, advancing clinical research in immune monitoring. Nevertheless, the vast data generated by cytometry by time-of-flight (CyTOF) poses a significant analytical challenge. To address this, we describe ImmCellTyper (https://github.com/JingAnyaSun/ImmCellTyper), a novel toolkit for CyTOF data analysis. This framework incorporates BinaryClust, an in-house developed semi-supervised clustering tool that automatically identifies main cell types. BinaryClust outperforms existing clustering tools in accuracy and speed, as shown in benchmarks with two datasets of approximately 4 million cells, matching the precision of manual gating by human experts. Furthermore, ImmCellTyper offers various visualisation and analytical tools, spanning from quality control to differential analysis, tailored to users' specific needs for a comprehensive CyTOF data analysis solution. The workflow includes five key steps: (1) batch effect evaluation and correction, (2) data quality control and pre-processing, (3) main cell lineage characterisation and quantification, (4) in-depth investigation of specific cell types; and (5) differential analysis of cell abundance and functional marker expression across study groups. Overall, ImmCellTyper combines expert biological knowledge in a semi-supervised approach to accurately deconvolute well-defined main cell lineages, while maintaining the potential of unsupervised methods to discover novel cell subsets, thus facilitating high-dimensional immune profiling.

Identifying high-risk multiple myeloma patients: A novel approach using a clonal gene signature.

Multiple myeloma (MM) is a heterogeneous disease with a small subset of high-risk patients having poor prognoses. Identifying these patients is crucial for treatment management and strategic decisions. In this study, we developed a novel computational framework to define prognostic gene signatures by selecting genes with expression driven by clonal copy number alterations. We applied this framework to MM and developed a clonal gene signature (CGS) consisting of 22 genes and evaluated in five independent datasets. The CGS provided significant prognostic values after adjusting for well-established factors including cytogenetic abnormalities, International Staging System (ISS), and Revised ISS (R-ISS). Importantly, CGS demonstrated higher performance in identifying high-risk patients compared to the GEP70 and SKY92 signatures recommended for prognostic stratification of MM. CGS can further stratify patients into subgroups with significantly differential prognoses when applied to the high- and low-risk groups identified by GEP70 and SKY92. Additionally, CGS scores are significantly associated with patient response to dexamethasone, a commonly used treatment for MM. In summary, we proposed a computational framework that requires only gene expression data to identify CGSs for prognosis prediction. CGS provides a useful biomarker for improving prognostic stratification in MM, especially for identifying the highest-risk patients.

A Data-Driven Computational Framework for Assessing the Risk of Placental Exposure to Environmental Chemicals.

A computational framework based on placental gene networks was proposed in this work to improve the accuracy of the placental exposure risk assessment of environmental compounds. The framework quantitatively characterizes the ability of compounds to cross the placental barrier by systematically considering the interaction and pathway-level information on multiple placental transporters. As a result, probability scores were generated for 307 compounds crossing the placental barrier based on this framework. These scores were then used to categorize the compounds into different levels of transplacental transport range, creating a gradient partition. These probability scores not only facilitated a more intuitive understanding of a compound's ability to cross the placental barrier but also provided valuable information for predicting potential placental disruptors. Compounds with probability scores greater than 90% were considered to have significant transplacental transport potential, whereas those with probability scores less than 80% were classified as unlikely to cross the placental barrier. Furthermore, external validation set results showed that the probability score could accurately predict the compounds known to cross the placental barrier. In conclusion, the computational framework proposed in this study enhances the intuitive understanding of the ability of compounds to cross the placental barrier and opens up new avenues for assessing the placental exposure risk of compounds.

Computational identification of DNA damage-relevant lncRNAs for predicting therapeutic efficacy and clinical outcomes in cancer.

OBJECTIVES: The role of long non-coding RNAs (lncRNAs) in cancer treatment, particularly in modulating DNA repair programs, is an emerging field that warrants systematic exploration. This study aimed to systematically identify the lncRNA regulators that potentially regulate DNA damage response (DDR). METHODS: Using genome-wide mRNA and lncRNA expression profiles of the same tumor patients, we proposed a novel computational framework. This framework performed Gene Set Variation Analysis to calculate DDR pathway enrichment score, which relies on weighting by tumor purity to obtain DDR activity score for each patient. Then, an in-depth differential expression profiling was conducted to identify pathway activity lncRNAs between high- and low-activity groups, utilizing a bootstrap-based randomization method. RESULTS: We comprehensively charted the landscape of DDR-relevant lncRNAs across 23 epithelial-based cancer types. Its effectiveness was validated by assessing the consistency of these lncRNAs within various datasets of the same cancer type (hypergeometric test P < 0.001), examining the expression perturbation of these lncRNAs in response to treatment and demonstrating its application in prioritizing clinically-related lncRNAs. Furthermore, leveraging 820 epithelial ovarian cancer patients from four public datasets, we applied these lncRNAs identified by DDRLnc to establish and validate a risk stratification model to evaluate the benefits of platinum-based adjuvant chemotherapy for the improvement of clinical treatment outcomes. CONCLUSIONS: Comprehensive pan-cancer analysis illustrates the utility of computational framework in identifying potentially lncRNAs involved in DDR, thereby offering novel insights into the complex realm of cancer research, and it will become a valuable tool for identifying potential therapeutic targets for cancer.

A benchmarking framework for the accurate and cost-effective detection of clinically-relevant structural variants for cancer target identification and diagnosis.

BACKGROUND: Accurate clinical structural variant (SV) calling is essential for cancer target identification and diagnosis but has been historically challenging due to the lack of ground truth for clinical specimens. Meanwhile, reduced clinical-testing cost is the key to the widespread clinical utility. METHODS: We analyzed massive data from tumor samples of 476 patients and developed a computational framework for accurate and cost-effective detection of clinically-relevant SVs. In addition, standard materials and classical experiments including immunohistochemistry and/or fluorescence in situ hybridization were used to validate the developed computational framework. RESULTS: We systematically evaluated the common algorithms for SV detection and established an expert-reviewed SV call set of 1,303 tumor-specific SVs with high-evidence levels. Moreover, we developed a random-forest-based decision model to improve the true positive of SVs. To independently validate the tailored 'two-step' strategy, we utilized standard materials and classical experiments. The accuracy of the model was over 90% (92-99.78%) for all types of data. CONCLUSION: Our study provides a valuable resource and an actionable guide to improve cancer-specific SV detection accuracy and clinical applicability.

A computational framework for improving genetic variants identification from 5,061 sheep sequencing data.

BACKGROUND: Pan-genomics is a recently emerging strategy that can be utilized to provide a more comprehensive characterization of genetic variation. Joint calling is routinely used to combine identified variants across multiple related samples. However, the improvement of variants identification using the mutual support information from multiple samples remains quite limited for population-scale genotyping. RESULTS: In this study, we developed a computational framework for joint calling genetic variants from 5,061 sheep by incorporating the sequencing error and optimizing mutual support information from multiple samples' data. The variants were accurately identified from multiple samples by using four steps: (1) Probabilities of variants from two widely used algorithms, GATK and Freebayes, were calculated by Poisson model incorporating base sequencing error potential; (2) The variants with high mapping quality or consistently identified from at least two samples by GATK and Freebayes were used to construct the raw high-confidence identification (rHID) variants database; (3) The high confidence variants identified in single sample were ordered by probability value and controlled by false discovery rate (FDR) using rHID database; (4) To avoid the elimination of potentially true variants from rHID database, the variants that failed FDR were reexamined to rescued potential true variants and ensured high accurate identification variants. The results indicated that the percent of concordant SNPs and Indels from Freebayes and GATK after our new method were significantly improved 12%-32% compared with raw variants and advantageously found low frequency variants of individual sheep involved several traits including nipples number (GPC5), scrapie pathology (PAPSS2), seasonal reproduction and litter size (GRM1), coat color (RAB27A), and lentivirus susceptibility (TMEM154). CONCLUSION: The new method used the computational strategy to reduce the number of false positives, and simultaneously improve the identification of genetic variants. This strategy did not incur any extra cost by using any additional samples or sequencing data information and advantageously identified rare variants which can be important for practical applications of animal breeding.

COIMMR: a computational framework to reveal the contribution of herbal ingredients against human cancer via immune microenvironment and metabolic reprogramming.

Immune evasion and metabolism reprogramming have been regarded as two vital hallmarks of the mechanism of carcinogenesis. Thus, targeting the immune microenvironment and the reprogrammed metabolic processes will aid in developing novel anti-cancer drugs. In recent decades, herbal medicine has been widely utilized to treat cancer through the modulation of the immune microenvironment and reprogrammed metabolic processes. However, labor-based herbal ingredient screening is time consuming, laborious and costly. Luckily, some computational approaches have been proposed to screen candidates for drug discovery rapidly. Yet, it has been challenging to develop methods to screen drug candidates exclusively targeting specific pathways, especially for herbal ingredients which exert anti-cancer effects by multiple targets, multiple pathways and synergistic ways. Meanwhile, currently employed approaches cannot quantify the contribution of the specific pathway to the overall curative effect of herbal ingredients. Hence, to address this problem, this study proposes a new computational framework to infer the contribution of the immune microenvironment and metabolic reprogramming (COIMMR) in herbal ingredients against human cancer and specifically screen herbal ingredients targeting the immune microenvironment and metabolic reprogramming. Finally, COIMMR was applied to identify isoliquiritigenin that specifically regulates the T cells in stomach adenocarcinoma and cephaelin hydrochloride that specifically targets metabolic reprogramming in low-grade glioma. The in silico results were further verified using in vitro experiments. Taken together, our approach opens new possibilities for repositioning drugs targeting immune and metabolic dysfunction in human cancer and provides new insights for drug development in other diseases. COIMMR is available at https://github.com/LYN2323/COIMMR.

Structural insights, biocatalytic characteristics, and application prospects of lignin-modifying enzymes for sustainable biotechnology.

Lignin modifying enzymes (LMEs) have gained widespread recognition in depolymerization of lignin polymers by oxidative cleavage. LMEs are a robust class of biocatalysts that include lignin peroxidase (LiP), manganese peroxidase (MnP), versatile peroxidase (VP), laccase (LAC), and dye-decolorizing peroxidase (DyP). Members of the LMEs family act on phenolic, non-phenolic substrates and have been widely researched for valorization of lignin, oxidative cleavage of xenobiotics and phenolics. LMEs implementation in the biotechnological and industrial sectors has sparked significant attention, although its potential future applications remain underexploited. To understand the mechanism of LMEs in sustainable pollution mitigation, several studies have been undertaken to assess the feasibility of LMEs in correlating to diverse pollutants for binding and intermolecular interactions at the molecular level. However, further investigation is required to fully comprehend the underlying mechanism. In this review we presented the key structural and functional features of LMEs, including the computational aspects, as well as the advanced applications in biotechnology and industrial research. Furthermore, concluding remarks and a look ahead, the use of LMEs coupled with computational framework, built upon artificial intelligence (AI) and machine learning (ML), has been emphasized as a recent milestone in environmental research.

Optimal locations and computational frameworks of FSR and IMU sensors for measuring gait abnormalities.

Neuromuscular diseases cause abnormal joint movements and drastically alter gait patterns in patients. The analysis of abnormal gait patterns can provide clinicians with an in-depth insight into implementing appropriate rehabilitation therapies. Wearable sensors are used to measure the gait patterns of neuromuscular patients due to their non-invasive and cost-efficient characteristics. FSR and IMU sensors are the most popular and efficient options. When assessing abnormal gait patterns, it is important to determine the optimal locations of FSRs and IMUs on the human body, along with their computational framework. The gait abnormalities of different types and the gait analysis systems based on IMUs and FSRs have therefore been investigated. After studying a variety of research articles, the optimal locations of the FSR and IMU sensors were determined by analysing the main pressure points under the feet and prime anatomical locations on the human body. A total of seven locations (the big toe, heel, first, third, and fifth metatarsals, as well as two close to the medial arch) can be used to measure gate cycles for normal and flat feet. It has been found that IMU sensors can be placed in four standard anatomical locations (the feet, shank, thigh, and pelvis). A section on computational analysis is included to illustrate how data from the FSR and IMU sensors are processed. Sensor data is typically sampled at 100 Hz, and wireless systems use a range of microcontrollers to capture and transmit the signals. The findings reported in this article are expected to help develop efficient and cost-effective gait analysis systems by using an optimal number of FSRs and IMUs.

Risky cascading transitions in international relationships.

Changing attitudes in diplomatic relations is a common feature of international politics. However, such changes may trigger risky domino-like cascades of "friend-to-enemy" transitions among other counties and yielding catastrophic damage that could reshape the global network of international relationships. While previous attention has been focused on studying single pairs of international relationships, due to the lack of a systematic framework, it remains still unknown whether, and how, a single transition of attitude between two countries could trigger a cascade of attitude transitions among other countries. Here, we develop such a framework and construct a global evolving network of relations between country pairs based on 70,756,728 international events between 1,225 country pairs from January 1995 to March 2020. Our framework can identify and quantify the cascade of transitions following a given original transition. Surprisingly, weaker transitions are found to initiate most of the largest cascades. We also find that transitions are not only related to the balance of the local environment, but also global network properties such as betweenness centrality. Our results suggest that these transitions have a substantial impact on bilateral trade volumes and scientific collaborations. Our results reveal reaction chains of international relations, which could be helpful for designing early warning signals and mitigation methods for global international conflicts.

Visualizing Fluid Flows via Regularized Optimal Mass Transport with Applications to Neuroscience.

The regularized optimal mass transport (rOMT) problem adds a diffusion term to the continuity equation in the original dynamic formulation of the optimal mass transport (OMT) problem proposed by Benamou and Brenier. We show that the rOMT model serves as a powerful tool in computational fluid dynamics for visualizing fluid flows in the glymphatic system. In the present work, we describe how to modify the previous numerical method for efficient implementation, resulting in a significant reduction in computational runtime. Numerical results applied to synthetic and real-data are provided.

A multi-reference poly-conformational method for in silico design, optimization, and repositioning of pharmaceutical compounds illustrated for selected SARS-CoV-2 ligands.

Background: This work presents a novel computational multi-reference poly-conformational algorithm for design, optimization, and repositioning of pharmaceutical compounds. Methods: The algorithm searches for candidates by comparing similarities between conformers of the same compound and identifies target compounds, whose conformers are collectively close to the conformers of each compound in the reference set. Reference compounds may possess highly variable MoAs, which directly, and simultaneously, shape the properties of target candidate compounds. Results: The algorithm functionality has been case study validated in silico, by scoring ChEMBL drugs against FDA-approved reference compounds that either have the highest predicted binding affinity to our chosen SARS-CoV-2 targets or are confirmed to be inhibiting such targets in-vivo. All our top scoring ChEMBL compounds also turned out to be either high-affinity ligands to the chosen targets (as confirmed in separate studies) or show significant efficacy, in-vivo, against those selected targets. In addition to method case study validation, in silico search for new compounds within two virtual libraries from the Enamine database is presented. The library's virtual compounds have been compared to the same set of reference drugs that we used for case study validation: Olaparib, Tadalafil, Ergotamine and Remdesivir. The large reference set of four potential SARS-CoV-2 compounds has been selected, since no drug has been identified to be 100% effective against the virus so far, possibly because each candidate drug was targeting only one, particular MoA. The goal here was to introduce a new methodology for identifying potential candidate(s) that cover multiple MoA-s presented within a set of reference compounds.

Breakdown of utilitarian moral judgement after basolateral amygdala damage.

Most of us would regard killing another person as morally wrong, but when the death of one saves multiple others, it can be morally permitted. According to a prominent computational dual-systems framework, in these life-and-death dilemmas, deontological (nonsacrificial) moral judgments stem from a model-free algorithm that emphasizes the intrinsic value of the sacrificial action, while utilitarian (sacrificial) moral judgments are derived from a model-based algorithm that emphasizes the outcome of the sacrificial action. Rodent decision-making research suggests that the model-based algorithm depends on the basolateral amygdala (BLA), but these findings have not yet been translated to human moral decision-making. Here, in five humans with selective, bilateral BLA damage, we show a breakdown of utilitarian sacrificial moral judgments, pointing at deficient model-based moral decision-making. Across an established set of moral dilemmas, healthy controls frequently sacrifice one person to save numerous others, but BLA-damaged humans withhold such sacrificial judgments even at the cost of thousands of lives. Our translational research confirms a neurocomputational hypothesis drawn from rodent decision-making research by indicating that the model-based algorithm which underlies outcome-based, utilitarian moral judgements in humans critically depends on the BLA.

Combating Viral Diseases in the Era of Systems Medicine.

Viruses can cause many diseases resulting in disabilities and death. Fortunately, advances in systems medicine enable the development of effective therapies for treating viral diseases, of vaccines to prevent viral infections, as well as of diagnostic tools to mitigate the risk and reduce the death toll. Characterizing the SARS-CoV-2 gene sequence and the role of its spike protein in infection informs development of small molecule drugs, antibodies, and vaccines to combat infection and complication, as well as to end the pandemic. Drug repurposing of small molecule drugs is a viable strategy to combat viral diseases; the key concepts include (1) linking a drug candidate's pharmacological network to its pharmacodynamic response in patients; (2) linking a drug candidate's physicochemical properties to its pharmacokinetic characteristics; and (3) optimizing the safe and effective dosing regimen within its therapeutic window. Computational integration of drug-induced signaling pathways with clinical outcomes is useful to inform selection of potential drug candidates with respect to safety and effectiveness. Key pharmacokinetic and pharmacodynamic principles for computational optimization of drug development include a drug candidate's Cminss/IC95 ratio, pharmacokinetic characteristics, and systemic exposure-response relationship, where Cminss is the trough concentration following multiple dosing. In summary, systems medicine approaches play a vital role in global success in combating viral diseases, including global real-time information sharing, development of test kits, drug repurposing, discovery and development of safe, effective therapies, detection of highly transmissible and deadly variants, and development of vaccines.

Integrative Analysis of Rhythmicity: From Biology to Urban Environments and Sustainability.

From biological to socio-technical systems, rhythmic processes are pervasive in our environment. However, methods for their comprehensive analysis are prevalent only in specific fields that limit the transfer of knowledge across scientific disciplines. This hinders interdisciplinary research and integrative analyses of rhythms across different domains and datasets. In this paper, we review recent developments in cross-disciplinary rhythmicity research, with a focus on the importance of rhythmic analyses in urban planning and biomedical research. Furthermore, we describe the current state of the art of (integrative) computational methods for the investigation of rhythmic data. Finally, we discuss the further potential and propose necessary future developments for cross-disciplinary rhythmicity analysis to foster integration of heterogeneous datasets across different domains, as well as guide data-driven decision making beyond the boundaries of traditional intradisciplinary research, especially in the context of sustainable and healthy cities.

Analysis of crystallization during rewarming in suboptimal vitrification conditions: a semi-empirical approach.

Circumventing ice formation is critical to successful cryopreservation by vitrification of large organs. While ice formation during the cooling part of the cryogenic protocol is dictated by the evolving thermal conditions, ice formation during the rewarming part of the cryogenic protocol is also dependent on the history of cooling and storage conditions. Furthermore, while the exothermic effect of ice crystallization during cooling tends to adversely slow down the desired high cooling rates to ensure ice-free preservation, the same effect under some conditions tends to assist acceleration of rewarming during recovery of the specimen from cryogenic storage when limited crystallization does occur. The current study proposes a computational framework to study the thermal effects of crystallization during recovery from cryogenic storage, using a semi-empirical approach to account for the relationship between latent heat effects and the rewarming rate. This study adds another layer of computational capabilities to a recent study investigating similar effects during cooling. Results of this study demonstrate that the thermal effects of crystallization on the local cooling and rewarming rates cannot be neglected. It further explains how crystallization during rewarming helps in increasing the rewarming rate and, thereby, affects rewarming-phase crystallization. Counterintuitively, this study suggests that the fastest possible rewarming rate at the outer surface of the domain in an inwards rewarming problem is not always advantageous, while the proposed computational tool is essential to find an intermediate optimal rate.

Progress towards a cellularly resolved mouse mesoconnectome is empowered by data fusion and new neuroanatomy techniques.

Over the past decade there has been a rapid improvement in techniques for obtaining large-scale cellular level data related to the mouse brain connectome. However, a detailed mapping of cell-type-specific projection patterns is lacking, which would, for instance, allow us to study the role of circuit motifs in cognitive processes. In this work, we review advanced neuroanatomical and data fusion techniques within the context of a proposed Multimodal Connectomic Integration Framework for augmenting the cellularly resolved mouse mesoconnectome. First, we emphasize the importance of registering data modalities to a common reference atlas. We then review a number of novel experimental techniques that can provide data for characterizing cell-types in the mouse brain. Furthermore, we examine a number of data integration strategies, which involve fine-grained cell-type classification, spatial inference of cell densities, latent variable models for the mesoconnectome and multi-modal factorisation. Finally, we discuss a number of use cases which depend on connectome augmentation techniques, such as model simulations of functional connectivity and generating mechanistic hypotheses for animal disease models.

Strategies for Vaccine Prioritization and Mass Dispensing.

We propose a system that helps decision makers during a pandemic find, in real time, the mass vaccination strategies that best utilize limited medical resources to achieve fast containments and population protection. Our general-purpose framework integrates into a single computational platform a multi-purpose compartmental disease propagation model, a human behavior network, a resource logistics model, and a stochastic queueing model for vaccination operations. We apply the modeling framework to the current COVID-19 pandemic and derive an optimal trigger for switching from a prioritized vaccination strategy to a non-prioritized strategy so as to minimize the overall attack rate and mortality rate. When vaccine supply is limited, such a mixed vaccination strategy is broadly effective. Our analysis suggests that delays in vaccine supply and inefficiencies in vaccination delivery can substantially impede the containment effort. Employing an optimal mixed strategy can significantly reduce the attack and mortality rates. The more infectious the virus, the earlier it helps to open the vaccine to the public. As vaccine efficacy decreases, the attack and mortality rates rapidly increase by multiples; this highlights the importance of early vaccination to reduce spreading as quickly as possible to lower the chances for further mutations to evolve and to reduce the excessive healthcare burden. To maximize the protective effect of available vaccines, of equal importance are determining the optimal mixed strategy and implementing effective on-the-ground dispensing. The optimal mixed strategy is quite robust against variations in model parameters and can be implemented readily in practice. Studies with our holistic modeling framework strongly support the urgent need for early vaccination in combating the COVID-19 pandemic. Our framework permits rapid custom modeling in practice. Additionally, it is generalizable for different types of infectious disease outbreaks, whereby a user may determine for a given type the effects of different interventions including the optimal switch trigger.

Uncovering Statistical Links Between Gene Expression and Structural Connectivity Patterns in the Mouse Brain.

Finding links between genes and structural connectivity is of the utmost importance for unravelling the underlying mechanism of the brain connectome. In this study we identify links between the gene expression and the axonal projection density in the mouse brain, by applying a modified version of the Linked ICA method to volumetric data from the Allen Institute for Brain Science for identifying independent sources of information that link both modalities at the voxel level. We performed separate analyses on sets of projections from the visual cortex, the caudoputamen and the midbrain reticular nucleus, and we determined those brain areas, injections and genes that were most involved in independent components that link both gene expression and projection density data, while we validated their biological context through enrichment analysis. We identified representative and literature-validated cortico-midbrain and cortico-striatal projections, whose gene subsets were enriched with annotations for neuronal and synaptic function and related developmental and metabolic processes. The results were highly reproducible when including all available projections, as well as consistent with factorisations obtained using the Dictionary Learning and Sparse Coding technique. Hence, Linked ICA yielded reproducible independent components that were preserved under increasing data variance. Taken together, we have developed and validated a novel paradigm for linking gene expression and structural projection patterns in the mouse mesoconnectome, which can power future studies aiming to relate genes to brain function.