Search

Abstract Pangenomics is a relatively new scientific field which investigates the union of all the genomes of a clade. The word pan means everything in ancient Greek; the term pangenomics originally regarded genomes of bacteria and was later intended to refer to human genomes as well. Modern bioinformatics offers several tools to analyze pangenomics data, paving the way to an emerging field that we can call computational pangenomics. Current computational power available for the bioinformatics community has made computational pangenomic analyses easy to perform, but this higher accessibility to pangenomics analysis also increases the chances to make mistakes and to produce misleading or inflated results, especially by beginners. To handle this problem, we present here a few quick tips for efficient and correct computational pangenomic analyses with a focus on bacterial pangenomics, by describing common mistakes to avoid and experienced best practices to follow in this field. We believe our recommendations can help the readers perform more robust and sound pangenomic analyses and to generate more reliable results.

Background: Hypersensitivity pneumonitis (HP) is an inflammatory disorder affecting lung parenchyma and often evolves into fibrosis (fHP). The altered regulation of genes involved in the pathogenesis of the disease is not well comprehended, while the role of microRNAs in lung fibroblasts remains unexplored. Methods: We used integrated bulk RNA-Seq and enrichment pathway bioinformatic analyses to identify differentially expressed (DE)-miRNAs and genes (DEGs) associated with HP lungs. In vitro, we evaluated the expression and potential role of miR-155-5p in the phenotype of fHP lung fibroblasts. Loss and gain assays were used to demonstrate the impact of miR-155-5p on fibroblast functions. In addition, mir-155-5p and its target TP53INP1 were analyzed after treatment with TGF-β, IL-4, and IL-17A. Results: We found around 50 DEGs shared by several databases that differentiate HP from control and IPF lungs, constituting a unique HP lung transcriptional signature. Additionally, we reveal 18 DE-miRNAs that may regulate these DEGs. Among the candidates likely associated with HP pathogenesis was miR-155-5p. Our findings indicate that increased miR-155-5p in fHP fibroblasts coincides with reduced TP53INP1 expression, high proliferative capacity, and a lack of senescence markers compared to IPF fibroblasts. Induced overexpression of miR-155-5p in normal fibroblasts remarkably increases the proliferation rate and decreases TP53INP1 expression. Conversely, miR-155-5p inhibition reduces proliferation and increases senescence markers. TGF-β, IL-4, and IL-17A stimulated miR-155-5p overexpression in HP lung fibroblasts. Conclusion: Our findings suggest a distinctive signature of 53 DEGs in HP, including CLDN18, EEF2, CXCL9, PLA2G2D, and ZNF683, as potential targets for future studies. Likewise, 18 miRNAs, including miR-155-5p, could be helpful to establish differences between these two pathologies. The overexpression of miR-155-5p and downregulation of TP53INP1 in fHP lung fibroblasts may be involved in his proliferative and profibrotic phenotype. These findings may help differentiate and characterize their pathogenic features and understand their role in the disease.

Abstract Adeno-associated viruses 2 (AAV2) are minute viruses renowned for their capacity to infect human cells and akin organisms. They have recently emerged as prominent candidates in the field of gene therapy, primarily attributed to their inherent non-pathogenic nature in humans and the safety associated with their manipulation. The efficacy of AAV2 as gene therapy vectors hinges on their ability to infiltrate host cells, a phenomenon reliant on their competence to construct a capsid capable of breaching the nucleus of the target cell. To enhance their infection potential, researchers have extensively scrutinized various combinatorial libraries by introducing mutations into the capsid, aiming to boost their effectiveness. The emergence of high-throughput experimental techniques, like deep mutational scanning (DMS), has made it feasible to experimentally assess the fitness of these libraries for their intended purpose. Notably, machine learning is starting to demonstrate its potential in addressing predictions within the mutational landscape from sequence data. In this context, we introduce a biophysically-inspired model designed to predict the viability of genetic variants in DMS experiments. This model is tailored to a specific segment of the CAP region within AAV2’s capsid protein. To evaluate its effectiveness, we conduct model training with diverse datasets, each tailored to explore different aspects of the mutational landscape influenced by the selection process. Our assessment of the biophysical model centers on two primary objectives: (i) providing quantitative forecasts for the log-selectivity of variants and (ii) deploying it as a binary classifier to categorize sequences into viable and non-viable classes.

This literature review rigorously examines the growing scientific interest in computational methods for Tissue Engineering and Regenerative Medicine biofabrication, a leading-edge area in biomedical innovation, emphasizing the need for accurate, multi-stage, and multi-component biofabrication process models. The paper presents a comprehensive bibliometric and contextual analysis, followed by a literature review, to shed light on the vast potential of computational methods in this domain. It reveals that most existing methods focus on single biofabrication process stages and components, and there is a significant gap in approaches that utilize accurate models encompassing both biological and technological aspects. This analysis underscores the indispensable role of these methods in understanding and effectively manipulating complex biological systems and the necessity for developing computational methods that span multiple stages and components. The review concludes that such comprehensive computational methods are essential for developing innovative and efficient Tissue Engineering and Regenerative Medicine biofabrication solutions, driving forward advancements in this dynamic and evolving field.

Objective(s): MicroRNAs, which are micro-coordinators of gene expression, have been recently investigated as a potential treatment for cancer. The study used computational techniques to identify microRNAs that could target a set of genes simultaneously. Due to their multi-target-directed nature, microRNAs have the potential to impact multiple key pathways and their pathogenic cross-talk. Materials and Methods: We identified microRNAs that target a prostate cancer-associated gene set using integrated bioinformatics analyses and experimental validation. The candidate gene set included genes targeted by clinically approved prostate cancer medications. We used STRING, GO, and KEGG web tools to confirm gene-gene interactions and their clinical significance. Then, we employed integrated predicted and validated bioinformatics approaches to retrieve hsa-miR-124-3p, 16-5p, and 27a-3p as the top three relevant microRNAs. KEGG and DIANA-miRPath showed the related pathways for the candidate genes and microRNAsResults: The Real-time PCR results showed that miR-16-5p simultaneously down-regulated all genes significantly except for PIK3CA/CB in LNCaP; miR-27a-3p simultaneously down-regulated all genes significantly, excluding MET in LNCaP and PIK3CA in PC-3; and miR-124-3p could not down-regulate significantly PIK3CB, MET, and FGFR4 in LNCaP and FGFR4 in PC-3. Finally, we used a cell cycle assay to show significant G0/G1 arrest by transfecting miR-124-3p in LNCaP and miR-16-5p in both cell lines. Conclusion: Our findings suggest that this novel approach may have therapeutic benefits and these predicted microRNAs could effectively target the candidate genes.

As artificial intelligence (AI) integrates within the intersecting domains of healthcare and computational biology, developing interpretable models tailored to medical contexts is met with significant challenges. Explainable AI (XAI) is vital for fostering trust and enabling effective use of AI in healthcare, particularly in image-based specialties such as pathology and radiology where adjunctive AI solutions for diagnostic image analysis are increasingly utilized. Overcoming these challenges necessitates interdisciplinary collaboration, essential for advancing XAI to enhance patient care. This commentary underscores the critical role of interdisciplinary conferences in promoting the necessary cross-disciplinary exchange for XAI innovation. A literature review was conducted to identify key challenges, best practices, and case studies related to interdisciplinary collaboration for XAI in healthcare. The distinctive contributions of specialized conferences in fostering dialogue, driving innovation, and influencing research directions were scrutinized. Best practices and recommendations for fostering collaboration, organizing conferences, and achieving targeted XAI solutions were adapted from the literature. By enabling crucial collaborative junctures that drive XAI progress, interdisciplinary conferences integrate diverse insights to produce new ideas, identify knowledge gaps, crystallize solutions, and spur long-term partnerships that generate high-impact research. Thoughtful structuring of these events, such as including sessions focused on theoretical foundations, real-world applications, and standardized evaluation, along with ample networking opportunities, is key to directing varied expertise toward overcoming core challenges. Successful collaborations depend on building mutual understanding and respect, clear communication, defined roles, and a shared commitment to the ethical development of robust, interpretable models. Specialized conferences are essential to shape the future of explainable AI and computational biology, contributing to improved patient outcomes and healthcare innovations. Recognizing the catalytic power of this collaborative model is key to accelerating the innovation and implementation of interpretable AI in medicine.

Abstract Background Computational models in systems biology are becoming more important with the advancement of experimental techniques to query the mechanistic details responsible for leading to phenotypes of interest. In particular, Boolean models are well fit to describe the complexity of signaling networks while being simple enough to scale to a very large number of components. With the advance of Boolean model inference techniques, the field is transforming from an artisanal way of building models of moderate size to a more automatized one, leading to very large models. In this context, adapting the simulation software for such increases in complexity is crucial. Results We present two new developments in the continuous time Boolean simulators: MaBoSS.MPI, a parallel implementation of MaBoSS which can exploit the computational power of very large CPU clusters, and MaBoSS.GPU, which can use GPU accelerators to perform these simulations. Conclusion These implementations enable simulation and exploration of the behavior of very large models, thus becoming a valuable analysis tool for the systems biology community.

Qiang Wang,* Pinduan Bi,* Ding Luo, Pingli Cao, Weihong Chen, Bin Yang Department of Hepatobiliary Surgery, the First Affiliated Hospital of Kunming Medical University, Kunming, People’s Republic of China*These authors contributed equally to this workCorrespondence: Bin Yang, Department of Hepatobiliary Surgery, the First Affiliated Hospital of Kunming Medical University, No. 295, Xichang Road, Wuhua District, Kunming, People’s Republic of China, Tel +86 013987642890, Email yangbin@kmmu.edu.cnBackground: Gallstone disease (GS) is an important risk factor for Gallbladder cancer (GBC). However, the mechanisms of the progression of GS to GBC remain unclear. Long non-coding RNA (lncRNA), modulates DNA/RNA/proteins at epigenetic, pre-transcriptional, transcriptional and posttranscriptional levels, and plays a potential therapeutic role in various diseases. This study aims to identify lncRNAs that have a potential impact on GS-promoted GBC progression.Methods and Results: Six GBC patients without GS, six normal gallbladder tissues, nine gallstones and nine GBC patients with GS were admitted to our hospital. The next-generation RNA-sequencing was performed to analyze differentially expressed (DE) lncRNA and messenger RNA (mRNA) in four groups. Then overlapping and specific molecular signatures were analyzed. We identified 29 co-DEGs and 500 co-DElncRNAs related to gallstone or GBC. The intersection and concatenation of co-DEGs and co-DElncRNA functionally involved in focal adhesion, Transcriptional misregulation in cancers, Protein digestion and absorption, and ECM-receptor interaction signaling pathways may contribute to the development of gallbladder cancer. Further exploration is necessary for early diagnosis and the potential treatment of GBC. FXYD2, MPZL1 and PAH were observed in both co-DEGs and co-DElncRNA and validated by qRT-PCR.Conclusion: Our data identified a series of DEGs and DElncRNAs, which were involved in the progression of GBC and GS-related metabolism pathways. Compared to GBC, the GS profile was more similar to para-tumor tissues in transcriptome level and lower risk of cancer. Further exploration is necessary from GBC patients with different periods of follow-up gallstone.Keywords: gallstone, gallbladder cancer, long non-coding RNA, next-generation sequencing analysis, computational biology

Abstract Background The Biology System Description Language (BiSDL) is an accessible, easy-to-use computational language for multicellular synthetic biology. It allows synthetic biologists to represent spatiality and multi-level cellular dynamics inherent to multicellular designs, filling a gap in the state of the art. Developed for designing and simulating spatial, multicellular synthetic biological systems, BiSDL integrates high-level conceptual design with detailed low-level modeling, fostering collaboration in the Design-Build-Test-Learn cycle. BiSDL descriptions directly compile into Nets-Within-Nets (NWNs) models, offering a unique approach to spatial and hierarchical modeling in biological systems. Results BiSDL’s effectiveness is showcased through three case studies on complex multicellular systems: a bacterial consortium, a synthetic morphogen system and a conjugative plasmid transfer process. These studies highlight the BiSDL proficiency in representing spatial interactions and multi-level cellular dynamics. The language facilitates the compilation of conceptual designs into detailed, simulatable models, leveraging the NWNs formalism. This enables intuitive modeling of complex biological systems, making advanced computational tools more accessible to a broader range of researchers. Conclusions BiSDL represents a significant step forward in computational languages for synthetic biology, providing a sophisticated yet user-friendly tool for designing and simulating complex biological systems with an emphasis on spatiality and cellular dynamics. Its introduction has the potential to transform research and development in synthetic biology, allowing for deeper insights and novel applications in understanding and manipulating multicellular systems.

3D cell culture models replicate tissue complexity and aim to study cellular interactions and responses in a more physiologically relevant environment compared to traditional 2D cultures. However, the spherical structure of these models makes it difficult to extract meaningful data, necessitating advanced techniques for proper analysis. In silico simulations enhance research by predicting cellular behaviors and therapeutic responses, providing a powerful tool to complement experimental approaches. Despite their potential, these simulations often require advanced computational skills and significant resources, which creates a barrier for many researchers. To address these challenges, we developed an accessible pipeline using open-source software to facilitate virtual tissue simulations. Our approach employs the Cellular Potts Model, a versatile framework for simulating cellular behaviors in tissues. The simulations are constructed from real world 3D image stacks of cancer spheroids, ensuring that the virtual models are rooted in experimental data. By introducing a new metric for parameter optimization, we enable the creation of realistic simulations without requiring extensive computational expertise. This pipeline benefits researchers wanting to incorporate computational biology into their methods, even if they do not possess extensive expertise in this area. By reducing the technical barriers associated with advanced computational modeling, our pipeline enables more researchers to utilize these powerful tools. Our approach aims to foster a broader use of in silico methods in disease research, contributing to a deeper understanding of disease biology and the refinement of therapeutic interventions.

Summary: Background: Melanoma is a heterogeneous cancer influenced by the plasticity of melanoma cells and their dynamic adaptations to microenvironmental cues. Melanoma cells transition between well-defined transcriptional cell states that impact treatment response and resistance. Methods: In this study, we applied single-cell RNA sequencing to interrogate the molecular features of immunotherapy-naive and immunotherapy-resistant melanoma tumours in response to ex vivo BRAF/MEK inhibitor treatment. Findings: We confirm the presence of four distinct melanoma cell states - melanocytic, transitory, neural-crest like and undifferentiated, and identify enrichment of neural crest-like and undifferentiated melanoma cells in immunotherapy-resistant tumours. Furthermore, we introduce an integrated computational approach to identify subsets of responding and nonresponding melanoma cells within the transcriptional cell states. Interpretation: Nonresponding melanoma cells are identified in all transcriptional cell states and are predisposed to BRAF/MEK inhibitor resistance due to pro-inflammatory IL6 and TNFɑ signalling. Our study provides a framework to study treatment response within distinct melanoma cell states and indicate that tumour-intrinsic pro-inflammatory signalling contributes to BRAF/MEK inhibitor resistance. Funding: This work was supported by Macquarie University, Melanoma Institute Australia, and the National Health and Medical Research Council of Australia (NHMRC; grant 2012860, 2028055).

Abstract The emergence of resistance to current antimalarial drugs poses a significant challenge in the fight against malaria. This study aimed to investigate the in vivo antiplasmodial potential of the aqueous extract of fresh and dried leaves of A3 in Plasmodium berghei‐infected (P. berghei) mice. A 4‐day suppressive test was conducted, with infected BALB/c mice receiving artesunate and A3 extracts. The results showed that the tested doses of A3 attenuated the elevation of parasitemia induced by P. berghei, particularly at the dose of 400 mg/kg, and improved hematological indices. Computational techniques, including molecular docking, binding free energy calculations, and ADMET predictions, identified several bioactive compounds in A3 with promising inhibitory potential against lysyl‐tRNA synthetases and Dihydrofolate reductase (DHFR), the crucial enzymes targeted by antimalarial drugs. In this paper, Friedelin, Bauerenol, Epifriedelanol, Alpha‐Amyrenone, Stigmasterol, and beta‐Amyrin acetate were top‐ranked, having docking scores from −10.6 to −9.9 kcal/mol, compared with the −9.4 and −7.1 kcal/mol demonstrated by artesunate and chloroquine, respectively, as standard ligands. Also, it was shown that docking score from the Lysyl‐tRNA protein target (4YCV) ranged from −9.5 to −7.8 kcal/mol in comparison to artesunate (8.1 kcal/mol) and chloroquine (5.6 kcal/mol). The results suggest that the identified compounds in A3 could serve as potential candidates for the development of new anti‐malarial agents.

Abstract The categorization of human diseases is mainly based on the affected organ system and phenotypic characteristics. This is limiting the view to the pathological manifestations, while it neglects mechanistic relationships that are crucial to develop therapeutic strategies. This work aims to advance the understanding of diseases and their relatedness beyond traditional phenotypic views. Hence, the similarity among 502 diseases is mapped using six different data dimensions encompassing molecular, clinical, and pharmacological information retrieved from public sources. Multiple distance measures and multi‐view clustering are used to assess the patterns of disease relatedness. The integration of all six dimensions into a consensus map of disease relationships reveals a divergent disease view from the International Classification of Diseases (ICD), emphasizing novel insights offered by a multi‐view disease map. Disease features such as genes, pathways, and chemicals that are enriched in distinct disease groups are identified. Finally, an evaluation of the top similar diseases of three candidate diseases common in the Western population shows concordance with known epidemiological associations and reveals rare features shared between Type 2 diabetes (T2D) and Alzheimer's disease. A revision of disease relationships holds promise for facilitating the reconstruction of comorbidity patterns, repurposing drugs, and advancing drug discovery in the future.

The hypothalamus in the mammalian brain is responsible for regulating functions associated with survival and reproduction representing a complex set of highly interconnected, yet anatomically and functionally distinct, sub-regions. It remains unclear what factors drive the spatial organization of sub-regions within the hypothalamus. One potential factor may be structural connectivity of the network that promotes efficient function with well-connected sub-regions placed closer together geometrically, i.e., the strongest axonal signal transferred through the shortest geometrical distance. To empirically test for such efficiency, we use hypothalamic data derived from the Allen Mouse Brain Connectivity Atlas, which provides a structural connectivity map of mouse brain regions derived from a series of viral tracing experiments. Using both cost function minimization and comparison with a weighted, sphere-packing ensemble, we demonstrate that the sum of the distances between hypothalamic sub-regions are not close to the minimum possible distance, consistent with prior whole brain studies. However, if such distances are weighted by the inverse of the magnitude of the connectivity, their sum is among the lowest possible values. Specifically, the hypothalamus appears within the top 94th percentile of neural efficiencies of randomly packed configurations and within one standard deviation of the median efficiency when packings are optimized for maximal neural efficiency. Our results, therefore, indicate that a combination of geometrical and topological constraints help govern the structure of the hypothalamus.

Background: The role of Mast cells has not been thoroughly explored in the context of prostate cancer's (PCA) unpredictable prognosis and mixed immunotherapy outcomes. Our research aims to employs a comprehensive computational methodology to evaluate Mast cell marker gene signatures (MCMGS) derived from a global cohort of 1091 PCA patients. This approach is designed to identify a robust biomarker to assist in prognosis and predicting responses to immunotherapy. Methods: This study initially identified mast cell-associated biomarkers from prostate adenocarcinoma (PRAD) patients across six international cohorts. We employed a variety of machine learning techniques, including Random Forest, Support Vector Machine (SVM), Lasso regression, and the Cox Proportional Hazards Model, to develop an effective MCMGS from candidate genes. Subsequently, an immunological assessment of MCMGS was conducted to provide new insights into the evaluation of immunotherapy responses and prognostic assessments. Additionally, we utilized Gene Set Enrichment Analysis (GSEA) and pathway analysis to explore the biological pathways and mechanisms associated with MCMGS. Results: MCMGS incorporated 13 marker genes and was successful in segregating patients into distinct high- and low-risk categories. Prognostic efficacy was confirmed by survival analysis incorporating MCMGS scores, alongside clinical parameters such as age, T stage, and Gleason scores. High MCMGS scores were correlated with upregulated pathways in fatty acid metabolism and β-alanine metabolism, while low scores correlated with DNA repair mechanisms, homologous recombination, and cell cycle progression. Patients classified as low-risk displayed increased sensitivity to drugs, indicating the utility of MCMGS in forecasting responses to immune checkpoint inhibitors. Conclusion: The combination of MCMGS with a robust machine learning methodology demonstrates considerable promise in guiding personalized risk stratification and informing therapeutic decisions for patients with PCA.

An increase in acute myocardial infarction (AMI)-related deaths has been reported during the COVID-19 pandemic. Despite evidence suggesting the association between severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection and AMI, the underlying mechanisms remain unclear. Here, we integrated mRNA and microRNA expression profiles related to SARS-CoV-2 infection and AMI from public databases. We then performed transcriptomic analysis using bioinformatics and systems biology approaches to explore the potential molecular mechanisms of SARS-CoV-2 infection affects AMI. First, twenty-one common differentially expressed genes (DEGs) were identified from SARS-CoV-2 infection and AMI patients in endothelial cells datasets and then we performed functional analysis to predict the roles of these DEGs. The functional analysis emphasized that the endothelial cell response to cytokine stimulus due to excessive inflammation was essential in these two diseases. Importantly, the tumor necrosis factor and interleukin-17 signaling pathways appeared to be integral factors in this mechanism. Interestingly, most of these common genes were also upregulated in transcriptomic datasets of SARS-CoV-2-infected cardiomyocytes, suggesting that these genes may be shared in cardiac- and vascular-related injuries. We subsequently built a protein-protein interaction network and extracted hub genes and essential modules from this network. At the transcriptional and post-transcriptional levels, regulatory networks with common DEGs were also constructed, and some key regulator signatures were further identified and validated. In summary, our research revealed that a highly activated inflammatory response in patients with COVID-19 might be a crucial factor for susceptibility to AMI and we identified some candidate genes and regulators that could be used as biomarkers or potential therapeutic targets.

Background: Platelet function is driven by the expression of specialized surface markers. The concept of distinct circulating subpopulations of platelets has emerged in recent years, but their exact nature remains debatable. Objectives: To design a spectral flow cytometry–based phenotyping workflow to provide a more comprehensive characterization, at a global and individual level, of surface markers in resting and activated healthy platelets, and to apply this workflow to investigate how responses differ according to platelet age. Methods: A 14-marker flow cytometry panel was developed and applied to vehicle- or agonist-stimulated platelet-rich plasma and whole blood samples obtained from healthy volunteers, or to platelets sorted according to SYTO-13 (Thermo Fisher Scientific) staining intensity as an indicator of platelet age. Data were analyzed using both user-led and independent approaches incorporating novel machine learning–based algorithms. Results: The assay detected differences in marker expression in healthy platelets, at rest and on agonist activation, in both platelet-rich plasma and whole blood samples, that are consistent with the literature. Machine learning identified stimulated populations of platelets with high accuracy (>80%). Similarly, machine learning differentiation between young and old platelet populations achieved 76% accuracy, primarily weighted by forward scatter, cluster of differentiation (CD) 41, side scatter, glycoprotein VI, CD61, and CD42b expression patterns. Conclusion: Our approach provides a powerful phenotypic assay coupled with robust bioinformatic and machine learning workflows for deep analysis of platelet subpopulations. Cleavable receptors, glycoprotein VI and CD42b, contribute to defining shared and unique subpopulations. This adoptable, low-volume approach will be valuable in deep characterization of platelets in disease.

Background Global public health is seriously threatened by the escalating issue of antimicrobial resistance (AMR). Antimicrobial peptides (AMPs), pivotal components of the innate immune system, have emerged as a potent solution to AMR due to their therapeutic potential. Employing computational methodologies for the prompt recognition of these antimicrobial peptides indeed unlocks fresh perspectives, thereby potentially revolutionizing antimicrobial drug development. Methods In this study, we have developed a model named as deepAMPNet. This model, which leverages graph neural networks, excels at the swift identification of AMPs. It employs structures of antimicrobial peptides predicted by AlphaFold2, encodes residue-level features through a bi-directional long short-term memory (Bi-LSTM) protein language model, and constructs adjacency matrices anchored on amino acids’ contact maps. Results In a comparative study with other state-of-the-art AMP predictors on two external independent test datasets, deepAMPNet outperformed in accuracy. Furthermore, in terms of commonly accepted evaluation matrices such as AUC, Mcc, sensitivity, and specificity, deepAMPNet achieved the highest or highly comparable performances against other predictors. Conclusion deepAMPNet interweaves both structural and sequence information of AMPs, stands as a high-performance identification model that propels the evolution and design in antimicrobial peptide pharmaceuticals. The data and code utilized in this study can be accessed at https://github.com/Iseeu233/deepAMPNet.

Erzhi Wan (EZW), a classic Traditional Chinese Medicine formula, has shown promise as a potential therapeutic option for Alzheimer's disease (AD), yet its mechanism remains elusive. Herein, we employed an integrative in-silico approach to investigate the active components and their mechanisms against AD. We screened four active components with blood-brain barrier permeabilities from TCMSP, along with 307 corresponding targets predicted by SwissTargetPrediction, PharmMapper, and TCMbank websites. Then, we retrieved 2260 AD-related targets from Genecards, OMIM, and NCBI databases. Furthermore, we constructed the protein-protein interaction (PPI) network of the intersected targets via the STRING database and performed the GO and KEGG enrichment analyses using the “clusterProfiler” R package. The results showed that the intersected targets were intimately related to the p53/PI3K/Akt signaling pathway, serotonergic synapse, and response to oxygen level. Subsequently, 25 core targets were found differentially expressed in brain regions by bioinformatics analyses of GEO datasets of clinical samples from the Alzdata database. The binding sites and stabilities between the active components and the core targets were investigated by the molecular docking approach using Autodock 4.2.6 software, followed by pocket detection and druggability assessment via the DoGSiteScorer server. The results showed that acacetin, β-sitosterol, and 3-O-acetyldammarenediol-II strongly interacted with the druggable pockets of AR, CASP8, POLB, and PREP. Eventually, the docking results were further cross-referenced with the literature research and validated by 100 ns of molecular dynamics simulations using GROMACS software. Binding free energies were calculated via MM/PBSA strategy combined with interaction entropy. The simulation results indicated stable bindings between four docking pairs including acacetin-AR, acacetin-CASP8, β-sitosterol-POLB, and 3-O-acetyldammarenediol–II–PREP. Overall, our study demonstrated a theoretical basis for how three active components of EZW confer efficacy against AD. It provides a promising reference for subsequent research regarding drug discoveries and clinical applications.

Similar to plant pathogens, phloem-feeding insects such as aphids deliver effector proteins inside their hosts that act to promote host susceptibility and enable feeding and infestation. Despite exciting progress toward identifying and characterizing effector proteins from these insects, their functions remain largely unknown. The recent groundbreaking development in protein structure prediction algorithms, combined with the availability of proteomics and transcriptomic datasets for agriculturally important pests, provides new opportunities to explore the structural and functional diversity of effector repertoires. In this study, we sought to gain insight into the infection strategy used by the Myzus persicae (green peach aphid) by predicting and analyzing the structures of a set of 71 effector candidate proteins. We used two protein structure prediction methods, AlphaFold and OmegaFold, that produced mutually consistent results. We observed a wide continuous spectrum of structures among the effector candidates, from disordered proteins to globular enzymes. We made use of the structural information and state-of-the-art computational methods to predict M. persicae effector protein properties, including function and interaction with host plant proteins. Overall, our investigation provides novel insights into prediction of structure, function, and interaction of M. persicae effector proteins and will guide the necessary experimental characterization to address new hypotheses. [Graphic: see text] Copyright © 2024 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

RNA-sequencing (RNA-seq) is a technique used for providing global patterns of transcriptomes in samples. However, it can only provide the average gene expression across cells and does not address the heterogeneity within the samples. The advances in single-cell RNA sequencing (scRNA-seq) technology have revolutionized our understanding of heterogeneity and the dynamics of gene expression at the single-cell level. For example, scRNA-seq allows us to identify the cell types in complex tissues, which can provide information regarding the alteration of the cell population by perturbations, such as genetic modification. Since its initial introduction, scRNA-seq has rapidly become popular, leading to the development of a huge number of bioinformatic tools. However, the analysis of the big dataset generated from scRNA-seq requires a general understanding of the preprocessing of the dataset and a variety of analytical techniques. Here, we present an overview of the workflow involved in analyzing the scRNA-seq dataset. First, we describe the preprocessing of the dataset, including quality control, normalization, and dimensionality reduction. Then, we introduce the downstream analysis provided with the most commonly used computational packages. This review aims to provide a workflow guideline for new researchers interested in this field.

Molecular interaction maps (MIMs) are static graphical representations depicting complex biochemical networks that can be formalized using one of the Systems Biology Graphical Notation languages. Regardless of their extensive coverage of various biological processes, they are limited in terms of dynamic insights. However, MIMs can serve as templates for developing dynamic computational models. We present MetaLo, an open-source Python package that enables the coupling of Boolean models inferred from process description MIMs with generic core metabolic networks. MetaLo provides a framework to study the impact of signaling cascades, gene regulation processes, and metabolic flux distribution of central energy production pathways. MetaLo computes the Boolean model’s asynchronous asymptotic behavior, through the identification of trap-spaces, and extracts metabolic constraints to contextualize the generic metabolic network. MetaLo is able to handle large-scale Boolean models and genome-scale metabolic models without requiring kinetic information or manual tuning. The framework behind MetaLo enables in depth analysis of the regulatory model, and may allow tackling a lack of omics data in poorly addressed biological fields to contextualize generic metabolic networks along with improper automatic reconstructions of cell- and/or disease-specific metabolic networks. MetaLo is available at https://pypi.org/project/metalo/ under the terms of the GNU General Public License v3.

Abstract Background Alzheimer’s disease (AD) is a progressive neurodegenerative disease that can cause dementia. We aim to screen out the hub genes involved in AD based on microarray datasets. Methods Gene expression profiles GSE5281 and GSE28146 were retrieved from Gene Expression Omnibus database to acquire differentially expressed genes (DEGs). Gene Ontology and pathway enrichment were conducted using DAVID online tool. The STRING database and Cytoscape tools were employed to analyze protein-protein interactions and identify hub genes. The predictive value of hub genes was assessed by principal component analysis and receiver operating characteristic curves. AD mice model was constructed, and histology was then observed by hematoxylin-eosin staining. Gene expression levels were finally determined by real-time quantitative PCR. Results We obtained 197 overlapping DEGs from GSE5281 and GSE28146 datasets. After constructing protein-protein interaction network, three highly interconnected clusters were identified and 6 hub genes (RBL1, BUB1, HDAC7, KAT5, SIRT2, and ITGB1) were selected. The hub genes could be used as basis to predict AD. Histological abnormalities of brain were observed, suggesting successful AD model was constructed. Compared with the control group, the mRNA expression levels of RBL1, BUB1, HDAC7, KAT5 and SIRT2 were significantly increased, while the mRNA expression level of ITGB1 was significantly decreased in AD groups. Conclusion RBL1, BUB1, HDAC7, KAT5, SIRT2 and ITGB1 are promising gene signatures for diagnosis and therapy of AD.

Biomedical science, which initially required only conventional research in the laboratory, currently involves information technology and has created bioinformatics in its development. Bioinformatics, a branch of biology, quantitatively analyzes information within biological macromolecules using software. Contemporary applications of bioinformatics have advanced biotechnological, medical, and pharmaceutical practices. Among the established applications of bioinformatics is diagnosing lung diseases using the genome-wide association study (GWAS) technique. Owing to sequencing technology and rapid computational methods, this technique is applied to analyze the link between genes with essential traits in the population, thus mapping the target genes to diagnose and treat diseases. The lung diseases diagnosed using GWAS include the responsible locus in asthma, chronic obstructive pulmonary disease (COPD), and lung cancer. Moreover, it can identify the treatment for COPD and suggest a new locus in lung cancer. Advancing the current gene-mapping technology demands genotype and phenotype data to study disease-linked genomes. Currently, bioinformatics is barely known and receives little attention in Indonesia. However, it can grow rapidly through open-source basis data and cross-disciplinary collaboration.

Abstract The road from transcription to protein synthesis is paved with many obstacles, allowing for several modes of post-transcriptional regulation of gene expression. A fundamental player in mRNA biology is DDX3X, an RNA binding protein that canonically regulates mRNA translation. By monitoring dynamics of mRNA abundance and translation following DDX3X depletion, we observe stabilization of translationally suppressed mRNAs. We use interpretable statistical learning models to uncover GC content in the coding sequence as the major feature underlying RNA stabilization. This result corroborates GC content-related mRNA regulation detectable in other studies, including hundreds of ENCODE datasets and recent work focusing on mRNA dynamics in the cell cycle. We provide further evidence for mRNA stabilization by detailed analysis of RNA-seq profiles in hundreds of samples, including a Ddx3x conditional knockout mouse model exhibiting cell cycle and neurogenesis defects. Our study identifies a ubiquitous feature underlying mRNA regulation and highlights the importance of quantifying multiple steps of the gene expression cascade, where RNA abundance and protein production are often uncoupled.

Abstract Multi-omics approaches have been successfully applied to investigate pregnancy and health outcomes at a molecular and genetic level in several studies. As omics technologies advance, research areas are open to study further. Here we discuss overall trends and examples of successfully using omics technologies and techniques (e.g., genomics, proteomics, metabolomics, and metagenomics) to investigate the molecular epidemiology of pregnancy. In addition, we outline omics applications and study characteristics of pregnancy for understanding fundamental biology, causal health, and physiological relationships, risk and prediction modeling, diagnostics, and correlations.

Abstract Background Drug repurposing plays a significant role in providing effective treatments for certain diseases faster and more cost-effectively. Successful repurposing cases are mostly supported by a classical paradigm that stems from de novo drug development. This paradigm is based on the “one-drug-one-target-one-disease” idea. It consists of designing drugs specifically for a single disease and its drug’s gene target. In this article, we investigated the use of biological pathways as potential elements to achieve effective drug repurposing. Methods Considering a total of 4214 successful cases of drug repurposing, we identified cases in which biological pathways serve as the underlying basis for successful repurposing, referred to as DREBIOP. Once the repurposing cases based on pathways were identified, we studied their inherent patterns by considering the different biological elements associated with this dataset, as well as the pathways involved in these cases. Furthermore, we obtained gene-disease association values to demonstrate the diminished significance of the drug’s gene target in these repurposing cases. To achieve this, we compared the values obtained for the DREBIOP set with the overall association values found in DISNET, as well as with the drug’s target gene (DREGE) based repurposing cases using the Mann-Whitney U Test. Results A collection of drug repurposing cases, known as DREBIOP, was identified as a result. DREBIOP cases exhibit distinct characteristics compared with DREGE cases. Notably, DREBIOP cases are associated with a higher number of biological pathways, with Vitamin D Metabolism and ACE inhibitors being the most prominent pathways. Additionally, it was observed that the association values of GDAs in DREBIOP cases were significantly lower than those in DREGE cases (p-value

Understanding the molecular and physical complexity of the tissue microenvironment (TiME) in the context of its spatiotemporal organization has remained an enduring challenge. Recent advances in engineering and data science are now promising the ability to study the structure, functions, and dynamics of the TiME in unprecedented detail; however, many advances still occur in silos that rarely integrate information to study the TiME in its full detail. This review provides an integrative overview of the engineering principles underlying chemical, optical, electrical, mechanical, and computational science to probe, sense, model, and fabricate the TiME. In individual sections, we first summarize the underlying principles, capabilities, and scope of emerging technologies, the breakthrough discoveries enabled by each technology and recent, promising innovations. We provide perspectives on the potential of these advances in answering critical questions about the TiME and its role in various disease and developmental processes. Finally, we present an integrative view that appreciates the major scientific and educational aspects in the study of the TiME.

Abstract Recent studies have revealed that numerous lncRNAs can translate proteins under specific conditions, performing diverse biological functions, thus termed coding lncRNAs. Their comprehensive landscape, however, remains elusive due to this field's preliminary and dispersed nature. This study introduces codLncScape, a framework for coding lncRNA exploration consisting of codLncDB, codLncFlow, codLncWeb, and codLncNLP. Specifically, it contains a manually compiled knowledge base, codLncDB, encompassing 353 coding lncRNA entries validated by experiments. Building upon codLncDB, codLncFlow investigates the expression characteristics of these lncRNAs and their diagnostic potential in the pan‐cancer context, alongside their association with spermatogenesis. Furthermore, codLncWeb emerges as a platform for storing, browsing, and accessing knowledge concerning coding lncRNAs within various programming environments. Finally, codLncNLP serves as a knowledge‐mining tool to enhance the timely content inclusion and updates within codLncDB. In summary, this study offers a well‐functioning, content‐rich ecosystem for coding lncRNA research, aiming to accelerate systematic studies in this field.

The vast biological and biochemical diversity of the global ocean is the driver behind marine bioprospecting for novel bioproducts. As Marine Biotechnology is gaining momentum as one of the main pillars of the ‘Brue Growth’ revolution, the ability to screen for novel compounds of interest in species with little or no genomic resources is paramount. With this respect, proteins, which are easily metabolised, can be synthetised using convenient DNA recombinant methods and can easily be modified to better meet the needs of human society, making them prized targets. Evidently, proteins that hold natural bioactivity and specificity such as toxins and other venom components, have long captured the focus of biotechnologists, leading to the merger between environmental omics and toxinology termed as ‘venomics’. Indeed, bioactive proteins such as conopeptides, conotoxins, turripeptides and others are long deemed important subjects of research. Even though current mainstream paradigms set the focus on secondary metabolites from marine organisms, transcriptomics and proteomics approaches and their combination are rising strategies for screening for thousands of proteins and peptides in non-conventional biological models, emphasising, but not limited to, marine invertebrate animals due to their abundance, biodiversity and uncanny biochemical strategies to cope with selective pressure in literally every known marine habitat. Untargeted approaches, such as RNA-Seq – based transcriptomics and tandem mass spectrometry – based proteomics, can circumvent limitations related with absent or reduced genomic annotation. The present review will outline the main contributions of ‘omics’ and computational approaches for bioprospecting for proteinaceous marine bioactives. Despite the relatively low number of ‘omics’ studies with the main purpose of discover novel compounds, there is already important literature showcasing pipelines and approaches for revolutionising the exploration of the ocean.

Non-Newtonian fluids flow generated by “cilia” are critical in medicine and bio-medical engineering. Such investigations are created by the back-and-forth movement of a microscopic hair-like structure connected to the walls, which causes a metachronal wave to form and drive biological fluids. Motivated by a wide range of biological applications, this study aims to explain the incompressible flow of Ellis fluid caused by the propagation of an infinite metachronal wave train traveling along channel walls owing to constantly beating cilia. The problem is simplified by low Reynolds number and long wavelength assumption. The mathematical model is solved with the aid of symbolic computational software Mathematica 13 version. The consequences of emerging parameters are then shown in graphical form and discussed comprehensively. From the study, it is worth mentioning that velocity declines with increasing material constants. Temperature distribution is also improved in the core sector of the channel and reduced at the walls. It is predicted that this approach will make an essential contribution to the progress and enhancement of various types of drug delivery systems in the biomedical industry and biomechanics.

We develop a thermodynamic model describing the binding of RNA binding proteins (RBP) to oligomers in vitro. We apply expectation-maximization to infer the specificity of RBPs, represented as position-specific weight matrices (PWMs), by maximizing the likelihood of RNA Bind’n Seq data from the ENCODE project. Analyzing these public data we find sequence motifs that can partly explain the data for more than half of the studied 111 RBPs, and for 48 of the proteins these motifs are consistent with the known specificity. Our code is publicly available, facilitating analysis of RBP binding data.

Tasmanian devils are endangered by a transmissible cancer known as Tasmanian devil facial tumour 1 (DFT1). A 2020 study by Patton et al. (Science 370, eabb9772 (doi:10.1126/science.abb9772)) used genome data from DFT1 tumours to produce a dated phylogenetic tree for this transmissible cancer lineage, and thence, using phylodynamics models, to estimate its epidemiological parameters and predict its future trajectory. It concluded that the effective reproduction number for DFT1 had declined to a value of one, and that the disease had shifted from emergence to endemism. We show that the study is based on erroneous mutation calls and flawed methodology, and that its conclusions cannot be substantiated.

Advances in high-throughput technologies have enhanced our ability to describe microbial communities as they relate to human health and disease. Alongside the growth in sequencing data has come an influx of resources that synthesize knowledge surrounding microbial traits, functions, and metabolic potential with knowledge of how they may impact host pathways to influence disease phenotypes. These knowledge bases can enable the development of mechanistic explanations that may underlie correlations detected between microbial communities and disease. In this review, we survey existing resources and methodologies for the computational integration of broad classes of microbial and host knowledge. We evaluate these knowledge bases in their access methods, content, and source characteristics. We discuss challenges of the creation and utilization of knowledge bases including inconsistency of nomenclature assignment of taxa and metabolites across sources, whether the biological entities represented are rooted in ontologies or taxonomies, and how the structure and accessibility limit the diversity of applications and user types. We make this information available in a code and data repository at: https://github.com/lozuponelab/knowledge-source-mappings. Addressing these challenges will allow for the development of more effective tools for drawing from abundant knowledge to find new insights into microbial mechanisms in disease by fostering a systematic and unbiased exploration of existing information.

Objective To investigate the expression of COL6A1 in sarcoma tissue and its significance based on bioinformatics analysis. Methods The GEO database was used to download the GSE16088, GSE21122, and GSE49972 datasets to analyze the expression level of COL6A1-3. The Kaplan-Meier survival curves were used to analyze the impact of the expression level of COL6A1-3 on the prognosis of patients with sarcoma. A Spearman correlation analysis was performed using GSVA package to investigate the correlation between the expression level of COL6A1 and immune cell infiltration. The LinkedOmics database was used to analyze the genes co-expressed with COL6A1, and the ggstatsplot package of R software was used to investigate the correlation between COL6A1 and COL6A2. The ClusterProfiler package of R software was used to perform GO functional enrichment analysis and KEGG pathway analysis of the genes co-expressed with COL6A1. The GSEA analysis was used to investigate the important kinase, miRNA, and transcription factor targets of the genes co-expressed with COL6A1. Results Compared with normal tissue, there was a significant increase in the expression level of COL6A1 in osteosarcoma, soft tissue sarcoma, and clear cell sarcoma of the kidney. The survival analysis showed that the high expression of COL6A1 was associated with low overall survival rate and di-sease-specific survival rate (χ2=6.37,4.57,P

As next-generation sequencing has been developed and used widely, RNA-sequencing (RNA-seq) has rapidly emerged as the first choice of tools to validate global transcriptome profiling. With the significant advances in RNA-seq, various types of RNA-seq have evolved in conjunction with the progress in bioinformatic tools. On the other hand, it is difficult to interpret the complex data underlying the biological meaning without a general understanding of the types of RNA-seq and bioinformatic approaches. In this regard, this paper discusses the two main sections of RNA-seq. First, two major variants of RNA-seq are described and compared with the standard RNA-seq. This provides insights into which RNA-seq method is most appropriate for their research. Second, the most widely used RNA-seq data analyses are discussed: (1) exploratory data analysis and (2) pathway enrichment analysis. This paper introduces the most widely used exploratory data analysis for RNA-seq, such as principal component analysis, heatmap, and volcano plot, which can provide the overall trends in the dataset. The pathway enrichment analysis section introduces three generations of pathway enrichment analysis and how they generate enriched pathways with the RNA-seq dataset.

Objective To lay the groundwork for formulating novel immunotherapies through understanding the regulatory role of ferroptosis and its function in immune microenvironment in ANCA-associated vasculitis (AAV) with renal injury. Methods Kidney data sets on AAV were downloaded from the database of Gene Expression Omnibus (GEO). Based upon the results of differential expression, consensus clustering algorithm was employed for identifying the subtypes of iron mortality. Single sample gene set enrichment analysis (ssGSEA) and gene set variation analysis (GSVA) were performed for evaluating the immune invasion and signaling pathways of subtypes. Key ferroptosis-related genes (FRGs) were further screened by LASSO and SVM-RFE algorithms and verified by independent data sets. CIBERSORT was utilized for assessing the abundance of infiltrating immune cells for FRGs. Finally a clinical diagnostic model of FRGs was constructed by rsm package. And the correlation between FRGs and renal function was examined by Nephroseq V5 database. Results Twenty-four ferroptosis related genes were significantly differentially expressed in AAV with kidney injury. Unsupervised cluster analysis revealed the regulatory patterns of two different ferroptosis subtypes with significant differences in gene expression, immune infiltration and biological functional pathways. A set of FRGs composed of ALOX5, CD44, CRYAB, TP53, and FTH1 was selected by two machine learning algorithms. It could effectively distinguish AAV with kidney injury from normal control group. CCR, macrophage and HLA were correlated positively with FRGs. Area under ROC curve of clinical diagnostic model constructed by FRGs was 0. 961. FRGs were correlated negatively with glomerular filtration rate (GFR). However, except for CRYAB and FTH1, other FRGs were correlated positively with serum creatinine (Scr). It implicated their potential relevance to renal injury in ANCA-associated vasculitis (AAV) and their potential as prognostic biomarkers. Conclusions Ferroptosis in AAV with kidney injury is correlated closely with immune microenvironment. It provides rationales for devising new treatments for AAV with kidney injury through a better understanding the pattern of ferroptosis.

Carcinogenesis is closely related to the expression, maintenance, and stability of DNA. These processes are regulated by one-carbon metabolism (1CM), which involves several vitamins of the complex B (folate, B2, B6, and B12), whereas alcohol disrupts the cycle due to the inhibition of folate activity. The relationship between nutrients related to 1CM (all aforementioned vitamins and alcohol) in breast cancer has been reviewed. The interplay of genes related to 1CM was also analyzed. Single nucleotide polymorphisms located in those genes were selected by considering the minor allele frequency in the Caucasian population and the linkage disequilibrium. These genes were used to perform several in silico functional analyses (considering corrected p-values < 0.05 as statistically significant) using various tools (FUMA, ShinyGO, and REVIGO) and databases such as the Kyoto Encyclopedia of Genes and Genomes (KEGG) and GeneOntology (GO). The results of this study showed that intake of 1CM-related B-complex vitamins is key to preventing breast cancer development and survival. Also, the genes involved in 1CM are overexpressed in mammary breast tissue and participate in a wide variety of biological phenomena related to cancer. Moreover, these genes are involved in alterations that give rise to several types of neoplasms, including breast cancer. Thus, this study supports the role of one-carbon metabolism B-complex vitamins and genes in breast cancer; the interaction between both should be addressed in future studies.

Although the disease caused by chikungunya virus (CHIKV) is of great interest to public health organizations around the world, there are still no authorized antivirals for its treatment. Previously, dihalogenated anti-CHIKV compounds derived from L-tyrosine (dH-Y) were identified as being effective against in vitro infection by this virus, so the objective of this study was to determine the mechanisms of its antiviral action. Six dH-Y compounds (C1 to C6) dihalogenated with bromine or chlorine and modified in their amino groups were evaluated by different in vitro antiviral strategies and in silico tools. When the cells were exposed before infection, all compounds decreased the expression of viral proteins; only C4, C5 and C6 inhibited the genome; and C1, C2 and C3 inhibited infectious viral particles (IVPs). Furthermore, C1 and C3 reduce adhesion, while C2 and C3 reduce internalization, which could be related to the in silico interaction with the fusion peptide of the E1 viral protein. Only C3, C4, C5 and C6 inhibited IVPs when the cells were exposed after infection, and their effect occurred in late stages after viral translation and replication, such as assembly, and not during budding. In summary, the structural changes of these compounds determine their mechanism of action. Additionally, C3 was the only compound that inhibited CHIKV infection at different stages of the replicative cycle, making it a compound of interest for conversion as a potential drug.

We herein introduce voyAGEr, an online graphical interface to explore age-related gene expression alterations in 49 human tissues. voyAGEr offers a visualisation and statistical toolkit for the finding and functional exploration of sex- and tissue-specific transcriptomic changes with age. In its conception, we developed a novel bioinformatics pipeline leveraging RNA sequencing data, from the GTEx project, encompassing more than 900 individuals. voyAGEr reveals transcriptomic signatures of the known asynchronous ageing between tissues, allowing the observation of tissue-specific age periods of major transcriptional changes, associated with alterations in different biological pathways, cellular composition, and disease conditions. Notably, voyAGEr was created to assist researchers with no expertise in bioinformatics, providing a supportive framework for elaborating, testing and refining their hypotheses on the molecular nature of human ageing and its association with pathologies, thereby also aiding in the discovery of novel therapeutic targets. voyAGEr is freely available at https://compbio.imm.medicina.ulisboa.pt/app/voyAGEr.

Cryogenic electron microscopy (cryo-EM) is a key method in structural and cell biology. Analysis of cryo-EM images requires interpretation of noisy, low-resolution densities, which relies on identifying the most probable orientation of macromolecules in a target using template matching. Many method-specific template-matching software solutions exist for single-particle cryo-EM, cryo-electron tomography (cryo-ET), or fitting atomic structures into averaged 3D maps of macromolecules. Here, we report the Python Template Matching Engine (pyTME), a software engine that consolidates method-specific template matching problems. The underlying library provides abstract data structures for storing and manipulating input and output data. PyTME runs up to ten times faster without loss in accuracy compared to existing software with multiple CPUs and GPUs, enabling template matching of even unbinned cryo-ET data in hours, which was previously nearly impossible due to technical constraints. Any hardware-specific optimization needed for dealing with large data is automatically performed to increase ease of use and minimize user intervention. The efficiency and simplicity of pyTME will enable high throughput mining of a variety of cryo-EM and ET datasets in the future.

Identification of genes whose expression increases or decreases with age is central to understanding the mechanisms behind aging. Recent scRNA-seq studies have shown that changes in single-cell expression profiles with aging are complex and diverse. In this study, we introduce a novel workflow to detect changes in the distribution of arbitrary monotonic age-related changes in single-cell expression profiles. Since single-cell gene expression profiles can be analyzed as probability distributions, our approach uses information theory to quantify the differences between distributions and employs distance matrices for association analysis. We tested this technique on simulated data and confirmed that potential parameter changes could be detected in a set of probability distributions. Application of the technique to a public scRNA-seq dataset demonstrated its potential utility as a straightforward screening method for identifying aging-related cellular features.

Antimicrobial resistance genes (ARG) are commonly found on acquired mobile genetic elements (MGEs) such as plasmids or transposons. Understanding the spread of resistance genes associated with mobile elements (mARGs) across different hosts and environments requires linking ARGs to the existing mobile reservoir within bacterial communities. However, reconstructing mARGs in metagenomic data from diverse ecosystems poses computational challenges, including genome fragment reconstruction (assembly), high-throughput annotation of MGEs, and identification of their association with ARGs. Recently, several bioinformatics tools have been developed to identify assembled fragments of plasmids, phages, and insertion sequence (IS) elements in metagenomic data. These methods can help in understanding the dissemination of mARGs. To streamline the process of identifying mARGs in multiple samples, we combined these tools in an automated high-throughput open-source pipeline, MetaMobilePicker, that identifies ARGs associated with plasmids, IS elements and phages, starting from short metagenomic sequencing reads. This pipeline was used to identify these three elements on a simplified simulated metagenome dataset, comprising whole genome sequences from seven clinically relevant bacterial species containing 55 ARGs, nine plasmids and five phages. The results demonstrated moderate precision for the identification of plasmids (0.57) and phages (0.71), and moderate sensitivity of identification of IS elements (0.58) and ARGs (0.70). In this study, we aim to assess the main causes of this moderate performance of the MGE prediction tools in a comprehensive manner. We conducted a systematic benchmark, considering metagenomic read coverage, contig length cutoffs and investigating the performance of the classification algorithms. Our analysis revealed that the metagenomic assembly process is the primary bottleneck when linking ARGs to identified MGEs in short-read metagenomics sequencing experiments rather than ARGs and MGEs identification by the different tools.

Phylogenetic comparative methods comprise the general endeavor of using an estimated phylogenetic tree (or set of trees) to make secondary inferences: about trait evolution, diversification dynamics, biogeography, community ecology, and a wide range of other phenomena or processes. Over the past ten years or so, the phytools R package has grown to become an important research tool for phylogenetic comparative analysis. phytools is a diverse contributed R library now consisting of hundreds of different functions covering a variety of methods and purposes in phylogenetic biology. As of the time of writing, phytools included functionality for fitting models of trait evolution, for reconstructing ancestral states, for studying diversification on trees, and for visualizing phylogenies, comparative data, and fitted models, as well numerous other tasks related to phylogenetic biology. Here, I describe some significant features of and recent updates to phytools, while also illustrating several popular workflows of the phytools computational software.

Preserving data privacy is an important concern in the research use of patient data. The DataSHIELD suite enables privacy-aware advanced statistical analysis in a federated setting. Despite its many applications, it has a few open practical issues: the complexity of hosting a federated infrastructure, the performance penalty imposed by the privacy-preserving constraints, and the ease of use by non-technical users. In this work, we describe a case study in which we review different breast cancer classifiers and report our findings about the limits and advantages of such non-disclosive suite of tools in a realistic setting. Five independent gene expression datasets of breast cancer survival were downloaded from Gene Expression Omnibus (GEO) and pooled together through the federated infrastructure. Three previously published and two newly proposed 5-year cancer-free survival risk score classifiers were trained in a federated environment, and an additional reference classifier was trained with unconstrained data access. The performance of these six classifiers was systematically evaluated, and the results show that i) the published classifiers do not generalize well when applied to patient cohorts that differ from those used to develop them; ii) among the methods we tried, the classification using logistic regression worked better on average, closely followed by random forest; iii) the unconstrained version of the logistic regression classifier outperformed the federated version by 4% on average. Reproducibility of our experiments is ensured through the use of VisualSHIELD, an open-source tool that augments DataSHIELD with new functions, a standardized deployment procedure, and a simple graphical user interface.

Computational models of disease progression have been constructed for a myriad of pathologies. Typically, the conceptual implementation for pathology-related in silico intervention studies has been ad hoc and similar in design to experimental studies. We introduce a multi-scale interventional design (MID) framework toward two key goals: tracking of disease dynamics from within-body to patient to population scale; and tracking impact(s) of interventions across these same spatial scales. Our MID framework prioritizes investigation of impact on individual patients within virtual pre-clinical trials, instead of replicating the design of experimental studies. We apply a MID framework to develop, organize, and analyze a cohort of virtual patients for the study of tuberculosis (TB) as an example disease. For this study, we use HostSim: our next-generation whole patient-scale computational model of individuals infected with Mycobacterium tuberculosis. HostSim captures infection within lungs by tracking multiple granulomas, together with dynamics occurring with blood and lymph node compartments, the compartments involved during pulmonary TB. We extend HostSim to include a simple drug intervention as an example of our approach and use our MID framework to quantify the impact of treatment at cellular and tissue (granuloma), patient (lungs, lymph nodes and blood), and population scales. Sensitivity analyses allow us to determine which features of virtual patients are the strongest predictors of intervention efficacy across scales. These insights allow us to identify patient-heterogeneous mechanisms that drive outcomes across scales.

Background: The RNA-dependent RNA polymerase (RdRp) complex, essential in viral transcription and replication, is a key target for antiviral therapeutics. The core unit of RdRp comprises the nonstructural protein NSP12, with NSP7 and two copies of NSP8 (NSP81 and NSP82) binding to NSP12 to enhance its affinity for viral RNA and polymerase activity. Notably, the interfaces between these subunits are highly conserved, simplifying the design of molecules that can disrupt their interaction.Methods: We conducted a detailed quantum biochemical analysis to characterize the interactions within the NSP12-NSP7, NSP12-NSP81, and NSP12-NSP82 dimers. Our objective was to ascertain the contribution of individual amino acids to these protein-protein interactions, pinpointing hotspot regions crucial for complex stability.Results: The analysis revealed that the NSP12-NSP81 complex possessed the highest total interaction energy (TIE), with 14 pairs of residues demonstrating significant energetic contributions. In contrast, the NSP12-NSP7 complex exhibited substantial interactions in 8 residue pairs, while the NSP12-NSP82 complex had only one pair showing notable interaction. The study highlighted the importance of hydrogen bonds and π-alkyl interactions in maintaining these complexes. Intriguingly, introducing the RNA sequence with Remdesivir into the complex resulted in negligible alterations in both interaction energy and geometric configuration.Conclusion: Our comprehensive analysis of the RdRp complex at the protein-protein interface provides invaluable insights into interaction dynamics and energetics. These findings can guide the design of small molecules or peptide/peptidomimetic ligands to disrupt these critical interactions, offering a strategic pathway for developing effective antiviral drugs.