A multimodal zebrafish developmental atlas reveals the state-transition dynamics of late-vertebrate pluripotent axial progenitors.

Elucidating organismal developmental processes requires a comprehensive understanding of cellular lineages in the spatial, temporal, and molecular domains. In this study, we introduce Zebrahub, a dynamic atlas of zebrafish embryonic development that integrates single-cell sequencing time course data with lineage reconstructions facilitated by light-sheet microscopy. This atlas offers high-resolution and in-depth molecular insights into zebrafish development, achieved through the sequencing of individual embryos across ten developmental stages, complemented by reconstructions of cellular trajectories. Zebrahub also incorporates an interactive tool to navigate the complex cellular flows and lineages derived from light-sheet microscopy data, enabling in silico fate-mapping experiments. To demonstrate the versatility of our multimodal resource, we utilize Zebrahub to provide fresh insights into the pluripotency of neuro-mesodermal progenitors (NMPs) and the origins of a joint kidney-hemangioblast progenitor population.

Mapping transcriptomic vector fields of single cells.

Single-cell (sc)RNA-seq, together with RNA velocity and metabolic labeling, reveals cellular states and transitions at unprecedented resolution. Fully exploiting these data, however, requires kinetic models capable of unveiling governing regulatory functions. Here, we introduce an analytical framework dynamo (https://github.com/aristoteleo/dynamo-release), which infers absolute RNA velocity, reconstructs continuous vector fields that predict cell fates, employs differential geometry to extract underlying regulations, and ultimately predicts optimal reprogramming paths and perturbation outcomes. We highlight dynamo's power to overcome fundamental limitations of conventional splicing-based RNA velocity analyses to enable accurate velocity estimations on a metabolically labeled human hematopoiesis scRNA-seq dataset. Furthermore, differential geometry analyses reveal mechanisms driving early megakaryocyte appearance and elucidate asymmetrical regulation within the PU.1-GATA1 circuit. Leveraging the least-action-path method, dynamo accurately predicts drivers of numerous hematopoietic transitions. Finally, in silico perturbations predict cell-fate diversions induced by gene perturbations. Dynamo, thus, represents an important step in advancing quantitative and predictive theories of cell-state transitions.

Cannabinoid receptor 1 antagonist genistein attenuates marijuana-induced vascular inflammation.

Epidemiological studies reveal that marijuana increases the risk of cardiovascular disease (CVD); however, little is known about the mechanism. Δ9-tetrahydrocannabinol (Δ9-THC), the psychoactive component of marijuana, binds to cannabinoid receptor 1 (CB1/CNR1) in the vasculature and is implicated in CVD. A UK Biobank analysis found that cannabis was an risk factor for CVD. We found that marijuana smoking activated inflammatory cytokines implicated in CVD. In silico virtual screening identified genistein, a soybean isoflavone, as a putative CB1 antagonist. Human-induced pluripotent stem cell-derived endothelial cells were used to model Δ9-THC-induced inflammation and oxidative stress via NF-κB signaling. Knockdown of the CB1 receptor with siRNA, CRISPR interference, and genistein attenuated the effects of Δ9-THC. In mice, genistein blocked Δ9-THC-induced endothelial dysfunction in wire myograph, reduced atherosclerotic plaque, and had minimal penetration of the central nervous system. Genistein is a CB1 antagonist that attenuates Δ9-THC-induced atherosclerosis.

Meiotic recombination mirrors patterns of germline replication in mice and humans.

Genetic recombination generates novel trait combinations, and understanding how recombination is distributed across the genome is key to modern genetics. The PRDM9 protein defines recombination hotspots; however, megabase-scale recombination patterning is independent of PRDM9. The single round of DNA replication, which precedes recombination in meiosis, may establish these patterns; therefore, we devised an approach to study meiotic replication that includes robust and sensitive mapping of replication origins. We find that meiotic DNA replication is distinct; reduced origin firing slows replication in meiosis, and a distinctive replication pattern in human males underlies the subtelomeric increase in recombination. We detected a robust correlation between replication and both contemporary and historical recombination and found that replication origin density coupled with chromosome size determines the recombination potential of individual chromosomes. Our findings and methods have implications for understanding the mechanisms underlying DNA replication, genetic recombination, and the landscape of mammalian germline variation.

Metabolic modeling of single Th17 cells reveals regulators of autoimmunity.

Metabolism is a major regulator of immune cell function, but it remains difficult to study the metabolic status of individual cells. Here, we present Compass, an algorithm to characterize cellular metabolic states based on single-cell RNA sequencing and flux balance analysis. We applied Compass to associate metabolic states with T helper 17 (Th17) functional variability (pathogenic potential) and recovered a metabolic switch between glycolysis and fatty acid oxidation, akin to known Th17/regulatory T cell (Treg) differences, which we validated by metabolic assays. Compass also predicted that Th17 pathogenicity was associated with arginine and downstream polyamine metabolism. Indeed, polyamine-related enzyme expression was enhanced in pathogenic Th17 and suppressed in Treg cells. Chemical and genetic perturbation of polyamine metabolism inhibited Th17 cytokines, promoted Foxp3 expression, and remodeled the transcriptome and epigenome of Th17 cells toward a Treg-like state. In vivo perturbations of the polyamine pathway altered the phenotype of encephalitogenic T cells and attenuated tissue inflammation in CNS autoimmunity.

Systematic Discovery of Short Linear Motifs Decodes Calcineurin Phosphatase Signaling.

Short linear motifs (SLiMs) drive dynamic protein-protein interactions essential for signaling, but sequence degeneracy and low binding affinities make them difficult to identify. We harnessed unbiased systematic approaches for SLiM discovery to elucidate the regulatory network of calcineurin (CN)/PP2B, the Ca2+-activated phosphatase that recognizes LxVP and PxIxIT motifs. In vitro proteome-wide detection of CN-binding peptides, in vivo SLiM-dependent proximity labeling, and in silico modeling of motif determinants uncovered unanticipated CN interactors, including NOTCH1, which we establish as a CN substrate. Unexpectedly, CN shows SLiM-dependent proximity to centrosomal and nuclear pore complex (NPC) proteins-structures where Ca2+ signaling is largely uncharacterized. CN dephosphorylates human and yeast NPC proteins and promotes accumulation of a nuclear transport reporter, suggesting conserved NPC regulation by CN. The CN network assembled here provides a resource to investigate Ca2+ and CN signaling and demonstrates synergy between experimental and computational methods, establishing a blueprint for examining SLiM-based networks.

Viral T-cell epitopes - Identification, characterization and clinical application.

T-cell immunity, mediated by CD4+ and CD8+ T cells, represents a cornerstone in the control of viral infections. Virus-derived T-cell epitopes are represented by human leukocyte antigen (HLA)-presented viral peptides on the surface of virus-infected cells. They are the prerequisite for the recognition of infected cells by T cells. Knowledge of viral T-cell epitopes provides on the one hand a diagnostic tool to decipher protective T-cell immune responses in the human population and on the other hand various prophylactic and therapeutic options including vaccination approaches and the transfer of virus-specific T cells. Such approaches have already been proven to be effective against various viral infections, particularly in immunocompromised patients lacking sufficient humoral, antibody-based immune response. This review provides an overview on the state of the art as well as current studies regarding the identification and characterization of viral T-cell epitopes and approaches of clinical application. In the first chapter in silico prediction tools and direct, mass spectrometry-based identification of viral T-cell epitopes is compared. The second chapter provides an overview of commonly used assays for further characterization of T-cell responses and phenotypes. The final chapter presents an overview of clinical application of viral T-cell epitopes with a focus on human immunodeficiency virus (HIV), human cytomegalovirus (HCMV) and severe acute respiratory syndrome coronavirus type 2 (SARS-CoV-2), being representatives of relevant viruses.

Advances in phage-host interaction prediction: in silico method enhances the development of phage therapies.

Phages can specifically recognize and kill bacteria, which lead to important application value of bacteriophage in bacterial identification and typing, livestock aquaculture and treatment of human bacterial infection. Considering the variety of human-infected bacteria and the continuous discovery of numerous pathogenic bacteria, screening suitable therapeutic phages that are capable of infecting pathogens from massive phage databases has been a principal step in phage therapy design. Experimental methods to identify phage-host interaction (PHI) are time-consuming and expensive; high-throughput computational method to predict PHI is therefore a potential substitute. Here, we systemically review bioinformatic methods for predicting PHI, introduce reference databases and in silico models applied in these methods and highlight the strengths and challenges of current tools. Finally, we discuss the application scope and future research direction of computational prediction methods, which contribute to the performance improvement of prediction models and the development of personalized phage therapy.

Synthesis of geometrically realistic and watertight neuronal ultrastructure manifolds for in silico modeling.

Understanding the intracellular dynamics of brain cells entails performing three-dimensional molecular simulations incorporating ultrastructural models that can capture cellular membrane geometries at nanometer scales. While there is an abundance of neuronal morphologies available online, e.g. from NeuroMorpho.Org, converting those fairly abstract point-and-diameter representations into geometrically realistic and simulation-ready, i.e. watertight, manifolds is challenging. Many neuronal mesh reconstruction methods have been proposed; however, their resulting meshes are either biologically unplausible or non-watertight. We present an effective and unconditionally robust method capable of generating geometrically realistic and watertight surface manifolds of spiny cortical neurons from their morphological descriptions. The robustness of our method is assessed based on a mixed dataset of cortical neurons with a wide variety of morphological classes. The implementation is seamlessly extended and applied to synthetic astrocytic morphologies that are also plausibly biological in detail. Resulting meshes are ultimately used to create volumetric meshes with tetrahedral domains to perform scalable in silico reaction-diffusion simulations for revealing cellular structure-function relationships. Availability and implementation: Our method is implemented in NeuroMorphoVis, a neuroscience-specific open source Blender add-on, making it freely accessible for neuroscience researchers.

Atomistic simulations reveal impacts of missense mutations on the structure and function of SynGAP1.

De novo mutations in the synaptic GTPase activating protein (SynGAP) are associated with neurological disorders like intellectual disability, epilepsy, and autism. SynGAP is also implicated in Alzheimer's disease and cancer. Although pathogenic variants are highly penetrant in neurodevelopmental conditions, a substantial number of them are caused by missense mutations that are difficult to diagnose. Hence, in silico mutagenesis was performed for probing the missense effects within the N-terminal region of SynGAP structure. Through extensive molecular dynamics simulations, encompassing three 150-ns replicates for 211 variants, the impact of missense mutations on the protein fold was assessed. The effect of the mutations on the folding stability was also quantitatively assessed using free energy calculations. The mutations were categorized as potentially pathogenic or benign based on their structural impacts. Finally, the study introduces wild-type-SynGAP in complex with RasGTPase at the inner membrane, while considering the potential effects of mutations on these key interactions. This study provides structural perspective to the clinical assessment of SynGAP missense variants and lays the foundation for future structure-based drug discovery.

Integrating network pharmacology and in silico analysis deciphers Withaferin-A's anti-breast cancer potential via hedgehog pathway and target network interplay.

This study examines the remarkable effectiveness of Withaferin-A (WA), a withanolide obtained from Withania somnifera (Ashwagandha), in encountering the mortiferous breast malignancy, a global peril. The predominant objective is to investigate WA's intrinsic target proteins and hedgehog (Hh) pathway proteins in breast cancer targeting through the application of in silico computational techniques and network pharmacology predictions. The databases and webtools like Swiss target prediction, GeneCards, DisGeNet and Online Mendelian Inheritance in Man were exploited to identify the common target proteins. The culmination of the WA network and protein-protein interaction network were devised using Stitch and String web tools, through which the drug-target network of 30 common proteins was constructed employing Cytoscape-version 3.9. Enrichment analysis was performed by incorporating Gprofiler, Metascape and Cytoscape plugins. David compounded the Gene Ontology and Kyoto Encyclopedia of Genes and Genomes, and enrichment was computed through bioinformatics tools. The 20 pivotal proteins were docked harnessing Glide, Schrodinger Suite 2023-2. The investigation was governed by docking scores and affinity. The shared target proteins underscored the precise Hh and WA network roles with the affirmation enrichment P-value of <0.025. The implications for hedgehog and cancer pathways were profound with enrichment (P < 0.01). Further, the ADMET and drug-likeness assessments assisted the claim. Robust interactions were noticed with docking studies, authenticated through molecular dynamics, molecular mechanics generalized born surface area scores and bonds. The computational investigation emphasized WA's credible anti-breast activity, specifically with Hh proteins, implying stem-cell-level checkpoint restraints. Rigorous testament is imperative through in vitro and in vivo studies.

BEERS2: RNA-Seq simulation through high fidelity in silico modeling.

Simulation of RNA-seq reads is critical in the assessment, comparison, benchmarking and development of bioinformatics tools. Yet the field of RNA-seq simulators has progressed little in the last decade. To address this need we have developed BEERS2, which combines a flexible and highly configurable design with detailed simulation of the entire library preparation and sequencing pipeline. BEERS2 takes input transcripts (typically fully length messenger RNA transcripts with polyA tails) from either customizable input or from CAMPAREE simulated RNA samples. It produces realistic reads of these transcripts as FASTQ, SAM or BAM formats with the SAM or BAM formats containing the true alignment to the reference genome. It also produces true transcript-level quantification values. BEERS2 combines a flexible and highly configurable design with detailed simulation of the entire library preparation and sequencing pipeline and is designed to include the effects of polyA selection and RiboZero for ribosomal depletion, hexamer priming sequence biases, GC-content biases in polymerase chain reaction (PCR) amplification, barcode read errors and errors during PCR amplification. These characteristics combine to make BEERS2 the most complete simulation of RNA-seq to date. Finally, we demonstrate the use of BEERS2 by measuring the effect of several settings on the popular Salmon pseudoalignment algorithm.

The emerging challenge of Enterococcus faecalis endocarditis after transcatheter aortic valve implantation: time for innovative treatment approaches.

SUMMARYInfective endocarditis (IE) is a life-threatening infection that has nearly doubled in prevalence over the last two decades due to the increase in implantable cardiac devices. Transcatheter aortic valve implantation (TAVI) is currently one of the most common cardiac procedures. TAVI usage continues to exponentially rise, inevitability increasing TAVI-IE. Patients with TAVI are frequently nonsurgical candidates, and TAVI-IE 1-year mortality rates can be as high as 74% without valve or bacterial biofilm removal. Enterococcus faecalis, a historically less common IE pathogen, is the primary cause of TAVI-IE. Treatment options are limited due to enterococcal intrinsic resistance and biofilm formation. Novel approaches are warranted to tackle current therapeutic gaps. We describe the existing challenges in treating TAVI-IE and how available treatment discovery approaches can be combined with an in silico "Living Heart" model to create solutions for the future.

Cortical cell size regulates root metabolic cost.

It has been hypothesized that vacuolar occupancy in mature root cortical parenchyma cells regulates root metabolic cost and thereby plant fitness under conditions of drought, suboptimal nutrient availability, and increased soil mechanical impedance. However, the mechanistic role of vacuoles in reducing root metabolic cost was unproven. Here we provide evidence to support this hypothesis. We first show that root cortical cell size is determined by both cortical cell diameter and cell length. Significant genotypic variation for both cortical cell diameter (~1.1- to 1.5-fold) and cortical cell length (~ 1.3- to 7-fold) was observed in maize and wheat. GWAS and QTL analyses indicate cortical cell diameter and length are heritable and under independent genetic control. We identify candidate genes for both phenes. Empirical results from isophenic lines contrasting for cortical cell diameter and length show that increased cell size, due to either diameter or length, is associated with reduced root respiration, nitrogen content, and phosphorus content. RootSlice, a functional-structural model of root anatomy, predicts that an increased vacuolar: cytoplasmic ratio per unit cortical volume causes reduced root respiration and nutrient content. Ultrastructural imaging of cortical parenchyma cells with varying cortical diameter and cortical cell length confirms the in silico predictions and shows that an increase in cell size is correlated with increased vacuolar volume and reduced cytoplasmic volume. Vacuolar occupancy and its relationship with cell size merits further investigation as a phene for improving crop adaptation to edaphic stress.

Calibration of computational tools for missense variant pathogenicity classification and ClinGen recommendations for PP3/BP4 criteria.

Recommendations from the American College of Medical Genetics and Genomics and the Association for Molecular Pathology (ACMG/AMP) for interpreting sequence variants specify the use of computational predictors as "supporting" level of evidence for pathogenicity or benignity using criteria PP3 and BP4, respectively. However, score intervals defined by tool developers, and ACMG/AMP recommendations that require the consensus of multiple predictors, lack quantitative support. Previously, we described a probabilistic framework that quantified the strengths of evidence (supporting, moderate, strong, very strong) within ACMG/AMP recommendations. We have extended this framework to computational predictors and introduce a new standard that converts a tool's scores to PP3 and BP4 evidence strengths. Our approach is based on estimating the local positive predictive value and can calibrate any computational tool or other continuous-scale evidence on any variant type. We estimate thresholds (score intervals) corresponding to each strength of evidence for pathogenicity and benignity for thirteen missense variant interpretation tools, using carefully assembled independent data sets. Most tools achieved supporting evidence level for both pathogenic and benign classification using newly established thresholds. Multiple tools reached score thresholds justifying moderate and several reached strong evidence levels. One tool reached very strong evidence level for benign classification on some variants. Based on these findings, we provide recommendations for evidence-based revisions of the PP3 and BP4 ACMG/AMP criteria using individual tools and future assessment of computational methods for clinical interpretation.

Comprehensive assessment of nine target prediction web services: which should we choose for target fishing?

Identification of potential targets for known bioactive compounds and novel synthetic analogs is of considerable significance. In silico target fishing (TF) has become an alternative strategy because of the expensive and laborious wet-lab experiments, explosive growth of bioactivity data and rapid development of high-throughput technologies. However, these TF methods are based on different algorithms, molecular representations and training datasets, which may lead to different results when predicting the same query molecules. This can be confusing for practitioners in practical applications. Therefore, this study systematically evaluated nine popular ligand-based TF methods based on target and ligand-target pair statistical strategies, which will help practitioners make choices among multiple TF methods. The evaluation results showed that SwissTargetPrediction was the best method to produce the most reliable predictions while enriching more targets. High-recall similarity ensemble approach (SEA) was able to find real targets for more compounds compared with other TF methods. Therefore, SwissTargetPrediction and SEA can be considered as primary selection methods in future studies. In addition, the results showed that k = 5 was the optimal number of experimental candidate targets. Finally, a novel ensemble TF method based on consensus voting is proposed to improve the prediction performance. The precision of the ensemble TF method outperforms the individual TF method, indicating that the ensemble TF method can more effectively identify real targets within a given top-k threshold. The results of this study can be used as a reference to guide practitioners in selecting the most effective methods in computational drug discovery.

Bioinformatics toolbox for exploring target mutation-induced drug resistance.

Drug resistance is increasingly among the main issues affecting human health and threatening agriculture and food security. In particular, developing approaches to overcome target mutation-induced drug resistance has long been an essential part of biological research. During the past decade, many bioinformatics tools have been developed to explore this type of drug resistance, and they have become popular for elucidating drug resistance mechanisms in a low cost, fast and effective way. However, these resources are scattered and underutilized, and their strengths and limitations have not been systematically analyzed and compared. Here, we systematically surveyed 59 freely available bioinformatics tools for exploring target mutation-induced drug resistance. We analyzed and summarized these resources based on their functionality, data volume, data source, operating principle, performance, etc. And we concisely discussed the strengths, limitations and application examples of these tools. Specifically, we tested some predictive tools and offered some thoughts from the clinician's perspective. Hopefully, this work will provide a useful toolbox for researchers working in the biomedical, pesticide, bioinformatics and pharmaceutical engineering fields, and a good platform for non-specialists to quickly understand drug resistance prediction.

Ultraliser: a framework for creating multiscale, high-fidelity and geometrically realistic 3D models for in silico neuroscience.

Ultraliser is a neuroscience-specific software framework capable of creating accurate and biologically realistic 3D models of complex neuroscientific structures at intracellular (e.g. mitochondria and endoplasmic reticula), cellular (e.g. neurons and glia) and even multicellular scales of resolution (e.g. cerebral vasculature and minicolumns). Resulting models are exported as triangulated surface meshes and annotated volumes for multiple applications in in silico neuroscience, allowing scalable supercomputer simulations that can unravel intricate cellular structure-function relationships. Ultraliser implements a high-performance and unconditionally robust voxelization engine adapted to create optimized watertight surface meshes and annotated voxel grids from arbitrary non-watertight triangular soups, digitized morphological skeletons or binary volumetric masks. The framework represents a major leap forward in simulation-based neuroscience, making it possible to employ high-resolution 3D structural models for quantification of surface areas and volumes, which are of the utmost importance for cellular and system simulations. The power of Ultraliser is demonstrated with several use cases in which hundreds of models are created for potential application in diverse types of simulations. Ultraliser is publicly released under the GNU GPL3 license on GitHub (BlueBrain/Ultraliser). SIGNIFICANCE: There is crystal clear evidence on the impact of cell shape on its signaling mechanisms. Structural models can therefore be insightful to realize the function; the more realistic the structure can be, the further we get insights into the function. Creating realistic structural models from existing ones is challenging, particularly when needed for detailed subcellular simulations. We present Ultraliser, a neuroscience-dedicated framework capable of building these structural models with realistic and detailed cellular geometries that can be used for simulations.

Mechanoresponsive regulation of tissue regeneration during distraction osteogenesis.

Distraction osteogenesis is widely used for bone tissue engineering. Mechanical stimulation plays a central role in the massive tissue regeneration observed during distraction osteogenesis. Although distraction osteogenesis has been a boon for patients with bone defects, we still have limited knowledge about the intrinsic mechanotransduction that converts physical forces into biochemical signals capable of inducing cell behavior changes and new tissue formation. In this review, we summarize the findings for mechanoresponsive factors, including cells, genes, and signaling pathways, during the distraction osteogenesis different phases. These elements function for coupling of osteogenesis and angiogenesis via the Integrin-FAK, TGF-ß/BMP, Wnt/ß-catenin, Hippo, MAPK, PI3K/Akt, and HIF-1α signaling pathways in a mechanoresponsive niche. The available evidence further suggests the existence of a balance between the epithelial-mesenchymal transition and mesenchymal-epithelial transition under hypoxic stress. We also briefly summarize the current in silico simulation algorithms and propose several future research directions that may advance understanding of distraction osteogenesis in the era of bioinformation, particularly the integration of artificial intelligence models with reliable single-cell RNA sequencing datasets. The objective of this review is to utilize established knowledge to further optimize existing distraction protocols and to identify potential therapeutic targets.

Fast, accurate ranking of engineered proteins by target-binding propensity using structure modeling.

Deep-learning-based methods for protein structure prediction have achieved unprecedented accuracy, yet their utility in the engineering of protein-based binders remains constrained due to a gap between the ability to predict the structures of candidate proteins and the ability toprioritize proteins by their potential to bind to a target. To bridge this gap, we introduce Automated Pairwise Peptide-Receptor Analysis for Screening Engineered proteins (APPRAISE), a method for predicting the target-binding propensity of engineered proteins. After generating structural models of engineered proteins competing for binding to a target using an established structure prediction tool such as AlphaFold-Multimer or ESMFold, APPRAISE performs a rapid (under 1 CPU second per model) scoring analysis that takes into account biophysical and geometrical constraints. As proof-of-concept cases, we demonstrate that APPRAISE can accurately classify receptor-dependent vs. receptor-independent adeno-associated viral vectors and diverse classes of engineered proteins such as miniproteins targeting the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike, nanobodies targeting a G-protein-coupled receptor, and peptides that specifically bind to transferrin receptor or programmed death-ligand 1 (PD-L1). APPRAISE is accessible through a web-based notebook interface using Google Colaboratory (https://tiny.cc/APPRAISE). With its accuracy, interpretability, and generalizability, APPRAISE promises to expand the utility of protein structure prediction and accelerate protein engineering for biomedical applications.