SEARCH · PubMed Health

Results for “Computer Simulation”

Explore indexed PubMed citations for clinical trials, systematic reviews and public health research. Read source abstracts and follow each citation to its original PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 recordsLinked to original sources

PubMed2025

Inbreeding load in finite populations from dominant and overdominant mutations.

Inbreeding depression is a widespread phenomenon that reflects the burden of deleterious effects hidden in heterozygosis in non-inbred populations but exposed in homozygosis in inbred individuals, known as inbreeding load (B). This load can be due to partially or fully recessive deleterious mutations (dominance model) or to heterozygote advantage (overdominance model, where both homozygotes are deleterious relative to the heterozygote). There are many studies addressing the changes in inbreeding load in finite populations assuming the dominance model. However, the contribution of overdominance to inbreeding depression has been focused on infinite-size populations. We carried out computer simulations to investigate the joint impact of dominant and pure overdominant mutations on inbreeding load, both for self-fertilizing populations and for panmictic populations suffering from a drastic bottleneck. We found that the overdominant inbreeding load can be substantially reduced by drift even for symmetrical overdominance, at least when considering mutations of small effect. For panmictic bottlenecked populations, the reduction in inbreeding load under dominance and overdominance loci cannot be easily distinguished. However, while purging depletes inbreeding load from dominant loci, slowing inbreeding depression and leading to partial fitness recovery, for overdominant loci fitness declines monotonically.

Inbreeding↗

PubMed2025

Competition and cooperation: The plasticity of bacterial interactions across environments.

Bacteria live in diverse communities, forming complex networks of interacting species. A central question in bacterial ecology is whether species engage in cooperative or competitive interactions. But this question often neglects the role of the environment. Here, we use genome-scale metabolic networks from two different open-access collections (AGORA and CarveMe) to assess pairwise interactions of different microbes in varying environmental conditions (provision of different environmental compounds). By computationally simulating thousands of environments for 10,000 pairs of bacteria from each collection, we found that most pairs were able to both compete and cooperate depending on the availability of environmental resources. This modeling approach allowed us to determine commonalities between environments that could facilitate the potential for cooperation or competition between a pair of species. Namely, cooperative interactions, especially obligate, were most common in less diverse environments. Further, as compounds were removed from the environment, we found interactions tended to degrade towards obligacy. However, we also found that on average at least one compound could be removed from an environment to switch the interaction from competition to facultative cooperation or vice versa. Together our approach indicates a high degree of plasticity in microbial interactions in response to the availability of environmental resources.

Microbial Interactions↗

PubMed2026

Elucidating the Mechanism of Xiaoqinglong Decoction in Chronic Urticaria Treatment: An Integrated Approach of Network Pharmacology, Bioinformatics Analysis, Molecular Docking, and Molecular Dynamics Simulations.

INTRODUCTION: Xiaoqinglong Decoction (XQLD) is a traditional Chinese medicinal formula commonly used to treat chronic urticaria (CU). However, its underlying therapeutic mechanisms remain incompletely characterized. This study employed an integrated approach combining network pharmacology, bioinformatics, molecular docking, and molecular dynamics simulations to identify the active components, potential targets, and related signaling pathways involved in XQLD's therapeutic action against CU, thereby providing a mechanistic foundation for its clinical application. METHODS: The active components of XQLD and their corresponding targets were identified using the Traditional Chinese Medicine Systems Pharmacology (TCMSP) database. CU-related targets were retrieved from the OMIM and GeneCards databases. Subsequently, core components and targets were determined via protein-protein interaction (PPI) network analysis and component-target-pathway network construction. Topological analyses were performed using Cytoscape software to prioritize core nodes within these networks. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were conducted via the DAVID database to identify enriched biological processes and signaling pathways. Molecular docking was performed to evaluate binding interactions between key components and core targets, while molecular dynamics (MD) simulations were employed to assess the stability of the component-target complexes with the lowest binding energy. Finally, CU-related targets of XQLD were validated using datasets from the Gene Expression Omnibus (GEO) database. RESULTS: A total of 135 active components and 249 potential targets of XQLD were identified, alongside 1,711 CU-related targets. Core components, such as quercetin, kaempferol, beta-sitosterol, naringenin, stigmasterol, and luteolin, exhibited high degree values in the constructed networks. The core targets identified included AKT1, TNF, IL6, TP53, PTGS2, CASP3, BCL2, ESR1, PPARG, and MAPK3. GO and KEGG pathway enrichment analyses revealed the PI3K-Akt signaling pathway as a central regulatory mechanism. Molecular docking studies demonstrated strong binding affinities between active components and core targets, with the stigmasterol-AKT1 complex exhibiting the lowest binding energy (-11.4 kcal/mol) and high stability in MD simulations. Validation using GEO datasets identified 12 core genes shared between CU-related targets and XQLD-associated targets, including PTGS2 and IL6, which were also prioritized as core targets in the network pharmacology analyses. DISCUSSION: This study comprehensively integrates multidisciplinary approaches to clarify the potential molecular mechanisms of XQLD in treating CU, highlighting its multitarget and multipathway synergistic effects. Molecular docking and dynamics simulations confirm the stable interaction between stigmasterol and the core target AKT1. Additionally, GEO dataset analysis verifies the pathogenic relevance of targets such as PTGS2 and IL6, significantly enhancing the credibility of our findings. These results provide a modern scientific basis for the traditional therapeutic effects of XQLD on CU and have important implications for developing multitarget treatments for this condition. However, this study mainly relies on database mining and computational simulations. Further in vitro and in vivo experimental validations are needed to confirm the predicted component-target-pathway interactions. CONCLUSION: This study identifies the active components, potential targets, and pathways through which XQLD exerts therapeutic effects on CU. These findings provide a theoretical foundation for further mechanistic studies and support their clinical application in the treatment of CU.

Molecular Docking Simulation↗

PubMed2025

Design of nanobody targeting SARS-CoV-2 spike glycoprotein using CDR-grafting assisted by molecular simulation and machine learning.

The design of proteins capable effectively binding to specific protein targets is crucial for developing therapies, diagnostics, and vaccine candidates for viral infections. Here, we introduce a complementarity-determining region (CDR) grafting approach for designing nanobodies (Nbs) that target specific epitopes, with the aid of computer simulation and machine learning. As a proof-of-concept, we designed, evaluated, and characterized a high-affinity Nb against the spike protein of SARS-CoV-2, the causative agent of the COVID-19 pandemic. The designed Nb, referred to as Nb Ab.2, was synthesized and displayed high-affinity for both the purified receptor-binding domain protein and to the virus-like particle, demonstrating affinities of 9 nM and 60 nM, respectively, as measured with microscale thermophoresis. Circular dichroism showed the designed protein's structural integrity and its proper folding, whereas molecular dynamics simulations provided insights into the internal dynamics of Nb Ab.2. This study shows that our computational pipeline can be used to efficiently design high-affinity Nbs with diagnostic and prophylactic potential, which can be tailored to tackle different viral targets.

Spike Glycoprotein, Coronavirus↗

PubMed2025

RadiSeq: a single- and bulk-cell whole-genome DNA sequencing simulator for radiation-damaged cell models.

Objective.To build and validate a simulation framework to perform single-cell and bulk-cell whole genome sequencing simulation of radiation-exposed Monte Carlo (MC) cell models to assist radiation genomics studies.Approach.Sequencing the genomes of radiation-damaged cells can provide useful insight into radiation action for radiobiology research. However, carrying out post-irradiation sequencing experiments can often be challenging, expensive, and time-consuming. Although computational simulations have the potential to provide solutions to these experimental challenges, and aid in designing optimal experiments, the absence of tools currently limits such application. MC toolkits exist to simulate radiation exposures of cell models but there are no tools to simulate single- and bulk-cell sequencing of cell models containing radiation-damaged DNA. Therefore, we aimed to develop a MC simulation framework to address this gap by designing a tool capable of simulating sequencing processes for radiation-damaged cells. Main results.We developed RadiSeq-a multi-threaded whole-genome DNA sequencing simulator written in C++. RadiSeq can be used to simulate Illumina sequencing of radiation-damaged cell models produced by MC simulations. RadiSeq has been validated through comparative analysis, where simulated data were matched against experimentally obtained data, demonstrating reasonable agreement between the two. Additionally, it comes with numerous features designed to closely resemble actual whole-genome sequencing. RadiSeq is also highly customizable with a single input parameter file.Significance.RadiSeq enables the research community to perform complex simulations of radiation-exposed DNA sequencing, supporting the optimization, planning, and validation of costly and time-intensive radiation biology experiments. This framework provides a powerful tool for advancing radiation genomics research.

Monte Carlo Method↗

PubMed2025

Multi-criteria decision making and its application to in silico discovery of vaccine candidates for Toxoplasma gondii.

Vaccine discovery against eukaryotic parasites is not trivial and few exist. Reverse vaccinology is an in silico vaccine discovery approach, designed to identify vaccine candidates from the thousands of protein sequences encoded by a target genome. Previously, we produced the Vacceed bioinformatics pipeline for identification of parasite membrane and excreted/secreted proteins that were likely be exposed to the hosts immune system. More recently, we improved upon machine learning as the final decision-making process to identify parasite proteins that induce a protective response in an animal model. Subsequently, we combined Vacceed with metrics on B and T cell epitope types to produce a new in silico discovery workflow. In this study we extend this in silico workflow to the developability of proteins as vaccines by the incorporation of metrics on the physicochemical properties of proteins. To demonstrate this process, every Toxoplasma gondii protein was ranked in its capacity to provide exposure to the immune system (Vacceed exposure score), presence of epitopes and solubility characteristics by several multicriteria decision making (MCDM) tools (such as TOPSIS, VIKOR and MABAC). A consensus rank was subsequently generated from the results of these tools using a variety of aggregate ranking methods. Levels of uncertainty in the aggregate protein rankings was assessed by conformal interval prediction in association with a machine learning model. Several of the top ranked proteins identified by this approach were novel, uncharacterized membrane transporters or proteins associated with RNA metabolism. In conclusion, MCDM automated the decision making using well known algorithms while conformal prediction intervals varied significantly across the 8000+ proteins of T. gondii. Highly ranked proteins (e.g. the top 100) typically generated low prediction intervals, providing high levels of confidence in their ranks.

Toxoplasma↗

PubMed2026

In silico analysis and comparison of the metabolic capabilities of different organisms by reducing metabolic complexity.

BACKGROUND: Understanding how metabolic capabilities diverge across microbial species is essential for deciphering community function, ecological interactions, and the design of synthetic microbiomes. Despite shared core pathways, microbial phenotypes can differ markedly due to evolutionary adaptations and metabolic specialization. Genome-scale metabolic models (GEMs) provide a systems-level framework to explore these differences; however, their complexity hinders direct comparison. RESULTS: We introduce NIS (Neidhardt-Ingraham-Schaechter), a computational workflow that integrates the redGEM, lumpGEM, and redGEMX algorithms to systematically reduce genome-scale models into biologically interpretable modules. This approach enables direct, quantitative comparison of fueling pathways, biomass biosynthetic routes, and environmental exchange processes while retaining essential metabolic information. We first demonstrate the utility of NIS by analyzing Escherichia coli and Saccharomyces cerevisiae, which revealed both conserved and divergent strategies in central metabolism, biosynthetic cost, and substrate utilization. We then applied NIS to the core honeybee gut microbiome, uncovering distinct metabolic traits, functional redundancy, and complementarity that help explain auxotrophy, cross-feeding interactions, and microbial coexistence. CONCLUSIONS: NIS provides an automated, scalable, and reproducible framework for dissecting microbial metabolic networks beyond gene content or taxonomy. By linking metabolism to ecological function, NIS offers new opportunities to interpret microbial community dynamics and to support the rational design of microbiomes in health, agriculture, and environmental applications. Video Abstract.

Metabolic Networks and Pathways↗

PubMed2026

Noncanonical bactericidal activity of teleost type I interferon is conferred by a membrane-targeting C-terminal peptide.

Type I interferons (IFNs) are indispensable antiviral cytokines in nonspecific immunity, yet they play dual roles in bacterial infections in mammals. Recent studies have revealed a subset of strongly cationic type I IFNs possessing potent antimicrobial properties across nonmammalian vertebrates. In this study, we identified a type I IFN gene, CaIFNi, from Cromileptes altivelis that is characterized by a unique triple-disulfide bond architecture. In Vibrio harveyi-challenged models, overexpression of CaIFNi potentiated bacterial clearance capacity in tissues, whereas its knockdown exacerbated bacterial colonization, highlighting its ability to protect the host against bacterial infection in vivo. In vitro assays further confirmed that CaIFNi directly binds to and kills both gram-negative (G-) and gram-positive (G+) bacteria, which first revealed the antibacterial function of new subgroup IFNi within teleost type I IFNs. Furthermore, the α-helical peptide CaIFNi-18 derived from CaIFNi was identified as a novel antimicrobial peptide (AMP) that has broad-spectrum antibacterial efficacy against G- and G+ bacteria and membrane-targeting ability. Further mechanistic studies revealed that CaIFNi has bactericidal effects on both G- and G+ bacteria through membrane depolarization and disruption, alteration of the bacterial ultrastructure, and in vitro binding to genomic DNA. In addition, CaIFNi-18 also has significant in vivo therapeutic efficacy against bacterial infection, highlighting its great potential as an antibacterial agent. Encouragingly, the loss of antibacterial activity in the truncation mutant (rCaIFNiΔ148-165) lacking the CaIFNi-18 segment suggests that this region is essential for the bactericidal function of the full-length protein and likely acts as its core domain. Further computational simulations revealed that the deletion of the CaIFNi-18 region attenuated the interaction between the protein and the bacterial membrane. These findings not only expand the functional scope of type I IFNs beyond their canonical antiviral role but also identify their derivative CaIFNi-18 as both a promising antimicrobial candidate and the essential bactericidal domain of CaIFNi, thereby offering novel therapeutic strategies against bacterial infections in the aquaculture industry and beyond.

Animals↗

PubMed2025

Modeling airborne transmission of viral genome using computational fluid dynamics simulation: A case study for SARS-CoV-2 virus.

Predicting indoor air quality during infectious disease conditions relies on models simulating particle materials (PM)/bioaerosols distribution. Understanding the thermo-fluid properties of exhaled air is crucial for comprehending disease transmission dynamics. This study employs a computational fluid dynamics (CFD) model to simulate cough-induced particle dispersion in a closed space. Furthermore, the number of released particles and the presence of SARS-CoV-2 viral genomes by a cough were assessed (in eight COVID-19 patients). According to the CFD model, in the first 30 s of cough, the vertical height and lateral breadth of the particles' dispersion were up to 138cm and 92cm, respectively. As the distance from the patient's respiratory zone increased, the lateral distribution width of particles expanded, reaching 1.3 m at 2.4 m away. Larger droplets (> 62.5µ) were deposited at shorter distances, while smaller particles remained airborne longer. The comparison of experimental and simulated results focused on particle dispersion at specific distances from the patient, particularly in the 2.5µ range. The distribution pattern of PM2.5 and PM10 at a distance of 1 and 2 m for women, not men, is similar to the distribution pattern of PM in CFD modeling. Viral genome detection was more prevalent in particles near the left side of the body, especially within the first 20 min post-cough, exhibiting a correlation with CFD predictions.

Airborne transmission↗

PubMed2025

In silico, in vitro, and in vivo characterization of thiamin-binding proteins from plant seeds.

Thiamin, an essential micronutrient, is a cofactor for enzymes involved in the central carbon metabolism and amino acid pathways. Despite efforts to enhance thiamin content in rice by incorporating thiamin biosynthetic genes, increasing thiamin content in the endosperm remains challenging, possibly due to a lack of thiamin stability and/or a local sink. The introduction of storage proteins has been successful in several biofortification strategies, and similar efforts targeting thiamin have been performed, leading to a 3-4-fold increase in white rice. However, only one thiamin-binding protein (TBP) sequence has been described in plants, more specifically from sesame seeds. Therefore, we aimed to identify and characterize TBPs, as well as to evaluate the effect of their expression on thiamin concentration, using a comprehensive approach integrating in silico, in vitro, and in vivo methods. We identified the sequences of putative TBPs from Oryza sativa (Os, rice), Fagopyrum esculentum (Fe, buckwheat), and Zea mays (Zm, maize) and pinpointed the thiamin-binding pockets through molecular docking. FeTBP and OsTBP contained one pocket with binding affinities similar to the Escherichia coli TBP, a well-characterized TBP, supporting their function as TBPs. In vivo expression studies of TBPs in tobacco leaves and rice callus resulted in increased thiamin levels, with FeTBP and OsTBP showing the most pronounced effects. Additionally, thermal shift assays confirmed the thiamin-binding capabilities of FeTBP and OsTBP, as observed by the significant increases in melting temperatures upon thiamin binding, indicating protein stabilization. These findings offer new insights into the diversity and function of plant TBPs and highlight the potential of FeTBP and OsTBP to modulate thiamin levels in crop plants.

Thiamine↗

PubMed2024

An integrated in-silico approach for drug target identification in human pathogen Shigella dysenteriae.

Shigella dysenteriae, is a Gram-negative bacterium that emerged as the second most significant cause of bacillary dysentery. Antibiotic treatment is vital in lowering Shigella infection rates, yet the growing global resistance to broad-spectrum antibiotics poses a significant challenge. The persistent multidrug resistance of S. dysenteriae complicates its management and control. Hence, there is an urgent requirement to discover novel therapeutic targets and potent medications to prevent and treat this disease. Therefore, the integration of bioinformatics methods such as subtractive and comparative analysis provides a pathway to compute the pan-genome of S. dysenteriae. In our study, we analysed a dataset comprising 27 whole genomes. The S. dysenteriae strain SD197 was used as the reference for determining the core genome. Initially, our focus was directed towards the identification of the proteome of the core genome. Moreover, several filters were applied to the core genome, including assessments for non-host homology, protein essentiality, and virulence, in order to prioritize potential drug targets. Among these targets were Integration host factor subunit alpha and Tyrosine recombinase XerC. Furthermore, four drug-like compounds showing potential inhibitory effects against both target proteins were identified. Subsequently, molecular docking analysis was conducted involving these targets and the compounds. This initial study provides the list of novel targets against S. dysenteriae. Conclusively, future in vitro investigations could validate our in-silico findings and uncover potential therapeutic drugs for combating bacillary dysentery infection.

Shigella dysenteriae↗

PubMed2026

Comparing ARG Inference Methods Under Transmission of Reproductive Success: Tree Imbalance Matters.

Inferring coalescent trees from genomic data has become a major subject in population genetics, particularly with the recent advances in tree sequence reconstruction methods. However, it remains unclear how well these methods perform for imbalanced genealogies. Such imbalances can arise from processes such as cultural transmission of reproductive success (CTRS) or positive selection. Using simulated genomic data, we benchmarked three major software packages, SINGER, Relate, and tsinfer, by comparing the imbalance of reconstructed trees by these methods with that of the true simulated trees, for three indices that quantify this imbalance. The three methods performed well under scenarios yielding balanced trees. However, their accuracy declined as imbalance increased. Performances also varied with mutation rate, recombination rate, and sample size. This study opens possibilities for applying these methods to infer CTRS or positive selection in large-scale genomic datasets, using simulation-based inference such as approximate Bayesian computation.

Models, Genetic↗

PubMed2025

CINner: Modeling and simulation of chromosomal instability in cancer at single-cell resolution.

Cancer development is characterized by chromosomal instability, manifesting in frequent occurrences of different genomic alteration mechanisms ranging in extent and impact. Mathematical modeling can help evaluate the role of each mutational process during tumor progression, however existing frameworks can only capture certain aspects of chromosomal instability (CIN). We present CINner, a mathematical framework for modeling genomic diversity and selection during tumor evolution. The main advantage of CINner is its flexibility to incorporate many genomic events that directly impact cellular fitness, from driver gene mutations to copy number alterations (CNAs), including focal amplifications and deletions, missegregations and whole-genome duplication (WGD). We apply CINner to find chromosome-arm selection parameters that drive tumorigenesis in the absence of WGD in chromosomally stable cancer types from the Pan-Cancer Analysis of Whole Genomes (PCAWG, [Formula: see text]). We found that the selection parameters predict WGD prevalence among different chromosomally unstable tumors, hinting that the selective advantage of WGD cells hinges on their tolerance for aneuploidy and escape from nullisomy. Analysis of inference results using CINner across cancer types in The Cancer Genome Atlas ([Formula: see text]) further reveals that the inferred selection parameters reflect the bias between tumor suppressor genes and oncogenes on specific genomic regions. Direct application of CINner to model the WGD proportion and fraction of genome altered (FGA) in PCAWG uncovers the increase in CNA probabilities associated with WGD in each cancer type. CINner can also be utilized to study chromosomally stable cancer types, by applying a selection model based on driver gene mutations and focal amplifications or deletions (chronic lymphocytic leukemia in PCAWG, [Formula: see text]). Finally, we used CINner to analyze the impact of CNA probabilities, chromosome selection parameters, tumor growth dynamics and population size on cancer fitness and heterogeneity. We expect that CINner will provide a powerful modeling tool for the oncology community to quantify the impact of newly uncovered genomic alteration mechanisms on shaping tumor progression and adaptation.

Chromosomal Instability↗

PubMed2026

In silico analysis based on network pharmacology and biomolecular informatics to explore the mechanism of action of Erjing Pills (from Shengji Zonglu) in the treatment of leukotrichia.

This study aimed to explore the core active ingredients and potential molecular mechanisms of Erjing Pills, a prescription in the classic work of Traditional Chinese Medicine, "Shengji Zonglu," in the treatment of leukotrichia by utilizing network pharmacology and biomolecular docking techniques. The chemical components and potential targets of Chinese herbal medicines were analyzed through databases such as the Traditional Chinese Medicine Systems Pharmacology Database. The targets related to leukotrichia were collected using GeneCards. The intersection targets were obtained using RStudio. The protein-protein interaction (PPI) network map and the "drug-component-target-disease" visualization network were generated using Cytoscape and STRING to screen the core components and key targets. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analyses were carried out using the Database for Annotation, Visualization and Integrated Discovery and RStudio. Finally, molecular docking verification was performed by AutoDock and PyMOL (Schrödinger LLC). The key active ingredients of Erjing Pills in the treatment of leukotrichia are β-sitosterol, quercetin, baicalein, and stigmasterol. The top 5 PPI core target proteins, in order, are AKT serine/threonine kinase 1, interleukin 6, tumor protein p53, cysteine-aspartic acid protease 3, and interleukin 1 beta. The Gene Ontology enrichment analysis suggests that the biological processes mainly include responses to exogenous stimuli, membrane rafts, and DNA-binding transcription factor binding. The Kyoto Encyclopedia of Genes and Genomes pathways involve signal pathways such as lipid and atherosclerosis, hepatitis B, Kaposi sarcoma virus infection, chemical carcinogenesis, and human cytomegalovirus infection. The molecular docking results indicate that most of the main active ingredients in Erjing Pills have relatively stable binding activities with the key targets, such as AKT serine/threonine kinase 1, interleukin 6, tumor protein p53, cysteine-aspartic acid protease 3, and interleukin 1 beta, in the PPI network. The active ingredients of Erjing Pills may interfere with the pathological process of leukotrichia by regulating key targets and signal pathways. This study provides a theoretical basis for the clinical application of Erjing Pills and indicates the direction for subsequent experimental research.

Drugs, Chinese Herbal↗

PubMed2025

vcfgl: a flexible genotype likelihood simulator for VCF/BCF files.

MOTIVATION: Accurate quantification of genotype uncertainty is pivotal in ensuring the reliability of genetic inferences drawn from NGS data. Genotype uncertainty is typically modeled using Genotype Likelihoods (GLs), which can help propagate measures of statistical uncertainty in base calls to downstream analyses. However, the effects of errors and biases in the estimation of GLs, introduced by biases in the original base call quality scores or the discretization of quality scores, as well as the choice of the GL model, remain under-explored. RESULTS: We present vcfgl, a versatile tool for simulating genotype likelihoods associated with simulated read data. It offers a framework for researchers to simulate and investigate the uncertainties and biases associated with the quantification of uncertainty, thereby facilitating a deeper understanding of their impacts on downstream analytical methods. Through simulations, we demonstrate the utility of vcfgl in benchmarking GL-based methods. The program can calculate GLs using various widely used genotype likelihood models and can simulate the errors in quality scores using a Beta distribution. It is compatible with modern simulators such as msprime and SLiM, and can output data in pileup, Variant Call Format (VCF)/BCF, and genomic VCF file formats, supporting a wide range of applications. The vcfgl program is freely available as an efficient and user-friendly software written in C/C++. AVAILABILITY AND IMPLEMENTATION: vcfgl is freely available at https://github.com/isinaltinkaya/vcfgl.

Software↗

PubMed2025

In silico encounters: harnessing metabolic modelling to understand plant-microbe interactions.

Understanding plant-microbe interactions is vital for developing sustainable agricultural practices and mitigating the consequences of climate change on food security. Plant-microbe interactions can improve nutrient acquisition, reduce dependency on chemical fertilizers, affect plant health, growth, and yield, and impact plants' resistance to biotic and abiotic stresses. These interactions are largely driven by metabolic exchanges and can thus be understood through metabolic network modelling. Recent developments in genomics, metagenomics, phenotyping, and synthetic biology now enable researchers to harness the potential of metabolic modelling at the genome scale. Here, we review studies that utilize genome-scale metabolic modelling to study plant-microbe interactions in symbiotic, pathogenic, and microbial community systems. This review catalogues how metabolic modelling has advanced our understanding of the plant host and its associated microorganisms as a holobiont. We showcase how these models can contextualize heterogeneous datasets and serve as valuable tools to dissect and quantify underlying mechanisms. Finally, we consider studies that employ metabolic models as a testbed for in silico design of synthetic microbial communities with predefined traits. We conclude by discussing broader implications of the presented studies, future perspectives, and outstanding challenges.

Plants↗

PubMed2025

insilicoSV: a flexible grammar-based framework for structural variant simulation and placement.

SUMMARY: Structural variants (SVs) are key drivers of genetic variation and disease in the genome. Their discovery remains challenging, however, in large part due to the scarcity of validated SV callsets and comprehensive benchmarks, which are essential for method development and evaluation. The growing number of data-driven learning-based approaches for SV discovery, in particular, requires large, diverse, and well-balanced training datasets to achieve reliable performance. To address this need, SV simulation has served as a key tool for assessing method performance and training SV models. However, existing SV simulators only support a fixed and limited set of SV classes and do not provide fine-grained control over the placement of SVs within specific contexts of the genome. Here we present insilicoSV, a versatile framework for SV simulation, which models SVs using a simple and flexible grammar, allowing users to easily define standard and custom arbitrary genome rearrangements, as well as encode genome placement constraints. This design allows insilicoSV to naturally support new and bespoke SV types, such as the complex rearrangements of cancer genomes. In addition to grammar-based modeling, insilicoSV provides built-in support for 26 predefined SV types, placement of user-provided SVs, small variant simulation, streamlined workflows for the simulation of genome evolution and genome mixtures, read simulation, alignment, and visualization. These features enable the creation of comprehensive genomic datasets for a variety of downstream applications, such as in-depth benchmarking of alignment and variant calling methods, as well as training of data-driven learning-based approaches for SV detection. AVAILABILITY AND IMPLEMENTATION: insilicoSV is available under the MIT license at https://github.com/PopicLab/insilicoSV and https://doi.org/10.5281/zenodo.17402009.

Software↗

PubMed2025

Accounting for contact tracing in epidemiological birth-death models.

Phylodynamics bridges the gap between classical epidemiology and pathogen genome sequence data by estimating epidemiological parameters from time-scaled pathogen phylogenetic trees. The models used in phylodynamics typically assume that the sampling procedure is independent between infected individuals. However, this assumption does not hold for many epidemics, in particular for such sexually transmitted infections as HIV-1, for which contact tracing schemes are included in health policies of many countries. We extended phylodynamic multi-type birth-death (MTBD) models with contact tracing (CT), and developed a simulator to generate trees under MTBD and MTBD-CT models. We proposed a non-parametric test for detecting contact tracing in pathogen phylogenetic trees. Its application to simulated data showed that it is both highly specific and sensitive. For the simplest representative of the MTBD-CT family, the BD-CT(1) model, where only the last contact can be notified, we solved the differential equations and proposed a closed form solution for the likelihood function. We implemented a maximum-likelihood program, which estimates the BD-CT(1) model parameters and their confidence intervals from phylogenetic trees. It performed accurate parameter inference on BD and BD-CT(1) simulated data, and detected contact tracing in HIV-1 B epidemics in Zurich and the UK. Importantly, we showed that not accounting for contact tracing when it is present, leads to bias in parameter estimation with the BD model (overestimation of the becoming-non-infectious rate). This bias is also present, but greatly reduced, when the BD-CT(1) model is used on data where multiple contacts can be notified. Our CT test, MTBD-CT tree simulator and BD-CT(1) parameter estimator are freely available at GitHub (evolbioinfo/treesimulator and evolbioinfo/bdct).

Contact Tracing↗