SEARCH · Search PubMed

Results for “Dynamic network”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

SourceResearch connectionFrom yearThrough year

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

At least 55 records · Page 3Linked to original sources

PubMed2023

Profiling of drug resistance in Src kinase at scale uncovers a regulatory network coupling autoinhibition and catalytic domain dynamics.

Kinase inhibitors are effective cancer therapies, but resistance often limits clinical efficacy. Despite the cataloging of numerous resistance mutations, our understanding of kinase inhibitor resistance is still incomplete. Here, we comprehensively profiled the resistance of ∼3,500 Src tyrosine kinase mutants to four different ATP-competitive inhibitors. We found that ATP-competitive inhibitor resistance mutations are distributed throughout Src's catalytic domain. In addition to inhibitor contact residues, residues that participate in regulating Src's phosphotransferase activity were prone to the development of resistance. Unexpectedly, we found that a resistance-prone cluster of residues located on the top face of the N-terminal lobe of Src's catalytic domain contributes to autoinhibition by reducing catalytic domain dynamics, and mutations in this cluster led to resistance by lowering inhibitor affinity and promoting kinase hyperactivation. Together, our studies demonstrate how drug resistance profiling can be used to define potential resistance pathways and uncover new mechanisms of kinase regulation.

src-Family Kinases↗

PubMed2024

Regulating IL-2 Immune Signaling Function Via A Core Allosteric Structural Network.

Human interleukin-2 (IL-2) is a crucial cytokine for T cell regulation, with therapeutic potential in cancer and autoimmune diseases. However, IL-2's pleiotropic effects across different immune cell types often lead to toxicity and limited efficacy. Previous efforts to enhance IL-2's therapeutic profile have focused on modifying its receptor binding sites. Yet, the underlying dynamics and intramolecular networks contributing to IL-2 receptor recognition remain unexplored. This study presents a detailed characterization of IL-2 dynamics compared to two engineered IL-2 mutants, "superkines" S15 and S1, which exhibit biased signaling towards effector T cells. Using NMR spectroscopy and molecular dynamics simulations, we demonstrate significant variations in core dynamic pathways and conformational exchange rates across these three IL-2 variants. We identify distinct allosteric networks and minor state conformations in the superkines, despite their structural similarity to wild-type IL-2. Furthermore, we rationally design a mutation (L56A) in the S1 superkine's core network, which partially reverts its dynamics, receptor binding affinity, and T cell signaling behavior towards that of wild-type IL-2. Our results reveal that IL-2 superkine core dynamics play a critical role in their enhanced receptor binding and function, suggesting that modulating IL-2 dynamics and core allostery represents an untapped approach for designing immunotherapies with improved immune cell selectivity profiles.

Interleukin-2↗

PubMed2024

Seasonal hydrological dynamics affected the diversity and assembly process of the antibiotic resistome in a canal network.

The significant threat of antibiotic resistance genes (ARGs) to aquatic environments health has been widely acknowledged. To date, several studies have focused on the distribution and diversity of ARGs in a single river while their profiles in complex river networks are largely known. Here, the spatiotemporal dynamics of ARG profiles in a canal network were examined using high-throughput quantitative PCR, and the underlying assembly processes and its main environmental influencing factors were elucidated using multiple statistical analyses. The results demonstrated significant seasonal dynamics with greater richness and relative abundance of ARGs observed during the dry season compared to the wet season. ARG profiles exhibited a pronounced distance-decay pattern in the dry season, whereas no such pattern was evident in the wet season. Null model analysis indicated that deterministic processes, in contrast to stochastic processes, had a significant impact on shaping the ARG profiles. Furthermore, it was found that Firmicutes and pH emerged as the foremost factors influencing these profiles. This study enhanced our comprehension of the variations in ARG profiles within canal networks, which may contribute to the design of efficient management approaches aimed at restraining the propagation of ARGs.

Seasons↗

PubMed1979

Effect of boundaries on the response of a neural network.

The effect an abrupt boundary has upon the dynamical response of a neural network is investigated. The retina of the Limulus eye is used as a model system for studying this effect. A theoretical technique is presented for the quantitative prediction of the manner in which this neural network responds in the vicinity of its boundary. Corresponding experimental measurements of the response to moving stimuli by single optic neurons located near retinal boundaries are presented. Theory and experiment show detailed quantitative agreement.

Animals↗

PubMed2026

Differential regulation of CYP46A1 in ischemic core and peri-infarct regions of male mouse brain after permanent middle cerebral artery occlusion.

Cholesterol 24-hydroxylase (CYP46A1) regulates brain cholesterol homeostasis and synaptic plasticity, playing a crucial role in ischemic stroke. Although previous studies have reported post-ischemic CYP46A1 upregulation, its spatiotemporal dynamics remain poorly defined. To elucidate these dynamics, we investigated the expression of CYP46A1 and other essential cholesterol homeostasis-related genes from 6 h to 3 days after permanent middle cerebral artery occlusion (pMCAO) in CB-17 mice. We utilized single-cell and single-nucleus transcriptomics, regional quantitative PCR, and high-resolution immunohistochemistry. CYP46A1 is predominantly expressed in neurons. Following ischemia, the cholesterol network exhibited a dynamic spatiotemporal divergence. Acutely (6 h post-ischemia), surviving regions transiently upregulated cell-autonomous cholesterol synthesis genes and CYP46A1. Subacutely (3 days), this response shifted toward a widespread upregulation of glia-dependent cholesterol transport genes and general CYP46A1 downregulation. At 24 h, CYP46A1 protein was substantially reduced in the necrotic core and superficial layer II/III of the peri-infarct cortex, but upregulated in deeper layer V, hippocampus, and lateral striatum. Notably, this localized upregulation spatially coincided with reactive microglial hypertrophy. These findings indicate that CYP46A1 is dynamically modulated in viable tissues following ischemic stress. This spatial divergence likely reflects a synergistic interaction between inflammatory propagation and neural circuit-mediated oxidative stress. Resolving these spatiotemporal profiles provides a rigorous foundation for evaluating CYP46A1 functionality and developing stage-specific therapeutic interventions.

Cholesterol 24-hydroxylase↗

PubMed2026

Human Systems Immunology in the Omics Era: Challenges, Methods, and Emerging Directions.

The human immune system is a highly complex, dynamic, and heterogeneous network shaped by genetic, environmental, and temporal influences. Advances in high-throughput omics technologies have transformed our ability to study this complexity directly and comprehensively in human cohorts. These developments have positioned systems immunology as a powerful framework for investigating coordinated immune responses, identifying regulatory mechanisms, and linking molecular patterns to clinical phenotypes. However, the analytical challenges inherent to large-scale, multimodal datasets-including batch effects, small sample sizes, high dimensionality, and substantial interindividual heterogeneity-require rigorous study design, robust statistical modeling, and thoughtful data analysis strategies. In this review, we summarize key technological foundations enabling modern human systems immunology, outline common analytical pitfalls and effective mitigation approaches, discuss data integration concepts, and highlight emerging opportunities in the field. Together, these technological and analytical advances are redefining how immune function is measured and interpreted in real-world human biology and hold significant promise for enhancing mechanistic insight, biomarker discovery, and precision medicine across immunological diseases and interventions.

Humans↗

PubMed2026

Urban stormwater infrastructure as a microplastic superhighway: a critical review of transport dynamics, modelling, and mitigation across pavements and drainage networks.

This review examines the transport, fate, modelling, and mitigation of Microplastics (MPs) in urban stormwater infrastructure, with emphasis on pavements, runoff pathways, micro-drainage, and macro-drainage systems. Following a systematic review approach, more than 1000 records were screened and approximately 50 core studies were retained when they addressed urban stormwater or drainage-related MP transport with adequate methodological reporting; marine-only studies and biological-effect studies without direct relevance to transport processes were excluded. The evidence shows that stormwater systems function not merely as passive conduits but as dynamic reactive transport systems with temporary storage, where particle mobilisation, sedimentation, resuspension, and temporary retention regulate MP export. Road surfaces, especially high-traffic areas, are major reservoirs of tyre wear, road-marking, atmospheric, and litter-derived particles that are rapidly mobilised during rainfall. Conventional grab sampling may underestimate MP loads, which in some cases exceed treated wastewater effluent loads by up to six-fold. Drainage structures such as manholes can immobilise up to 17.3% of near-neutrally buoyant particles, while biofouling and aggregation may shift buoyant polymers from wash-load to bedload. Mitigation systems, including permeable pavements, bioretention, wetlands, and technical inserts, can achieve high removal of coarse MPs, but performance declines for fine particles below 100 µm. The review highlights the need for standardised flow-proportional sampling, physically informed modelling, and treatment-train strategies targeting both surface sources and in-network storage.

Microplastics↗

PubMed1979

Spatial firing patterns of auditory neuron network modelling by computer simulation.

This communication examines, in digital computer simulated network, input signals and response patterns established at excitatory neurons' level i.e. the membrane potential of neuron soma. It is restricted to spatial patterns of the auditory neuron networks and time factor for nervous conduction and transmission is neglected compared with long maintained membrane potentials of neuron somas. The model analyzes the change in the spatial patterns of the membrane potential in the two dimensional networks of the auditory system. In order to evaluate the contribution of the various parameters, it is started that the simplest model has only one parameter, lateral inhibition. The other parameters are then added, one at a time, to successive models. The lateral inhibition is a necessary condition in the auditory nervous system if any sharpening of the response areas in the single neurons is to occur. A necessary condition for the validity of the model is that is should be applicable to the other senses such as vision and chemical patterns, taste. The threshold feature of auditory neurons aids in producing a sharpening in the neuron of the auditory relay nuclei. It does this clipping the spatial response patterns in one dimensional arrays of excitatory neurons. Recurrent inhibition seems a necessary condition in the sensory nervous system that any kinds of input signals are to be preserved over a wide range of stimulus intensity. In other words, this network has a wide dynamic range against any kinds of input signals. A simple self-recurrent negative feedback does not contribute to the sharpening, but more complex socalled averaged type does. A neuron network is capable of responding stably to stimuli with a wide range of intensity and with any kind of spatial patterns if there is a simple negative feedback mechanism. When there is no negative feedback, input signals soon disappear or saturate in the neuron network. Therefore, recurrent inhibition is the most important mechanism. Spontaneous activity appears to aid in the sharpening by providing a kind of contrast, that is by reducting the amount of activity in neurons adjacent to the excitatory area. Moreover, the effect of spontaneous activity in the model seems to make repples around the excitatory area and suggests that an introduction of activity at any stage of the networks, from whatever source for example reticulum formation and thalamus, might appreciably alter the response patterns at subsequent neuron network. This suggests that the mechanism of the consciousness that might be controlled by the thalamus and or reticular formation. These two dimensional neuron networks may be expanded to three dimensional neuron networks. The former might simulate the auditory nervous system while the latter might simulate the visual system.

Animals↗

PubMed2025

MNMO: discover driver genes from a multi-omics data based-multi-layer network.

MOTIVATION: Cancer as a public health problem is driven by genomic variations in "cancer driver" genes. The identification of driver genes is critical for the discovery of key biomarkers and the development of personalized therapy. RESULTS: We propose a prediction method MNMO: a multi-layer network model based on multi-omics data. MNMO firstly constructs a dynamically adjusted four-layer network composed of miRNAs and three kinds of genes with different features. Then three kinds of scores, i.e. control capacity, mutation score, and network score, are devised and calculated by harmonic mean to produce the integrated gene score. Experiments were performed on three kinds of real cancer data to compare the identification performance of method MNMO with that of six state-of-the-art ones. The results indicate that method MNMO presents the best identification performance under most circumstances. The genes prioritized by method MNMO not only have a better match to the benchmark ones than those identified by the other methods, but also are all associated with the development and progression of cancers. In addition, some extended versions of method MNMO can further achieve better performance on most evaluation metrics for some specific datasets. They may be more conducive to identifying tissue-specific genes, which has been verified through a number of experiments. AVAILABILITY AND IMPLEMENTATION: The source code and the R package "MNMO" are available at https://github.com/Zheng-D/MNMO. The dataset and code are archived at https://doi.org/10.5281/zenodo.14969986.

Humans↗

PubMed2025

Cross-feeding percolation phase transitions of intercellular metabolic networks.

Intercellular cross-talk is essential for the adaptation capabilities of populations of cells. While direct diffusion-driven cell-to-cell exchanges are difficult to map, current nanotechnology enables one to probe single-cell exchanges with the medium. We introduce a mathematical method to reconstruct the dynamic unfolding of intercellular exchange networks from these data, applying it to an experimental coculture system. The exchange network, initially dense, progressively fragments into small disconnected clusters. To explain these dynamics, we develop a maximum-entropy multicellular metabolic model with diffusion-driven exchanges. The model predicts a transition from a dense network to a sparse one as nutrient consumption shifts. We characterize this crossover both numerically, revealing a power-law decay in the cluster-size distribution, and analytically, by connecting to percolation theory. Comparison with data suggests that populations evolve toward the sparse phase by remaining near the crossover. These findings offer insights into the collective organization driving the adaptive dynamics of cell populations.

Metabolic Networks and Pathways↗

PubMed2026

Resolving cellular signaling in space and time: From organelle proteomics to spatial phosphoproteomics.

Cellular signaling is inherently organized in space and time, requiring coordinated control of protein localization, molecular interactions, and enzymatic activity across subcellular compartments. Recent advances in chemical biology, protein engineering, and quantitative proteomics have made it possible to interrogate these dimensions in an integrated manner. Here, we highlight emerging strategies to resolve signaling organization across three interconnected dimensions: organelle-resolved proteome mapping to define spatial context, proximity labeling to capture local protein interaction networks, and spatially resolved phosphoproteomics to quantify signaling outputs. Developments in proximity labeling, including split, conditionally activated and light-gated enzymes, enable temporally controlled, context-dependent profiling of transient protein assemblies in living cells. Advances in high-throughput and low-input phosphoproteomics, together with improved computational frameworks for kinase activity inference and subcellular enrichment strategies, are enabling spatially resolved measurement of signaling activity. Together, these approaches are shifting the field from static localization maps toward dynamic models of signaling networks.

Proteomics↗

PubMed2026

Quantifying uncertainty of predictions from cancer progression models.

MOTIVATION: Cancer progresses through the accumulation of genomic events. Cancer progression models such as Mutual Hazard Networks (MHNs) describe this dynamic, enabling prediction of temporal event positions and patient-specific risks of acquiring mutations. However, current MHN analyses rely on single most likely models and do not quantify the uncertainty inherent to parameter estimation. Assessing forecast stability is essential before using them to anticipate treatment-relevant mutations, adapt targeted therapies, or prioritize monitoring of patients at elevated progression risk. RESULTS: We address a key prerequisite for the responsible clinical use of cancer progression models by making MHN-derived predictions uncertainty-aware. We present a Bayesian framework for MHN that uses Markov Chain Monte Carlo to sample from the posterior distributions of model parameters and derived predictions. For practical use we implemented the Random-Walk Metropolis, Metropolis-Adjusted Langevin Algorithm (MALA), and simplified manifold MALA samplers as part of the existing mhn Python package. Only MALA and smMALA were successful in sampling from MHN posteriors, with MALA performing best. While most MHN parameters and predictions showed low posterior variance, a small subset displayed greater variability across the posterior distribution. This differentiation cannot be obtained from a single most likely model, emphasizing the need for uncertainty quantification, especially in clinical contexts. As an illustrative example, posterior sampling identified a subgroup of STK11$-$, KRAS$+$ lung adenocarcinoma patients with a high predicted short-term risk-with low variance across posterior samples-to develop an STK11 mutation. This subgroup exhibited poorer survival under immunotherapy, resembling patterns observed in STK11+ patients. AVAILABILITY AND IMPLEMENTATION: Our implementation is part of version 1.2.0 of the mhn package (https://github.com/spang-lab/LearnMHN). All analyses including the code to produce all figures in this article can be found under https://github.com/huy29433/MCMC-sampling-for-MHN (https://doi.org/10.5281/zenodo.21160219).

Humans↗

PubMed2026

The Mechanism of Celosiae Semen in the Treatment of Diabetic Cataract: Based on Network Pharmacology.

INTRODUCTION: Diabetes mellitus can be complicated by a variety of ocular diseases, among which the postoperative complications of diabetic cataract (DC) are significantly higher than those of non-DC patients. Therefore, finding drugs with natural active ingredients is an urgent challenge in the prevention and treatment of DC. Discovering the potential molecular mechanism of celosiae semen (CS) for the treatment of DC and providing new ideas and programs for the treatment and prevention of DC. METHODS: In this study, network pharmacology, molecular docking, and molecular dynamics simulations were utilized to predict the binding and functional enrichment of the main active ingredients of CS with DC-related targets, and to explore the potential pathways and mechanisms of CS for the treatment of DC. RESULTS: Through database searching and screening, a total of 45 potential targets of CS for the treatment of DC were identified, functionally enriched, and a protein-protein interaction network was constructed, and the key target, SRC, was finally found. The results of molecular docking and molecular dynamics simulation showed that the main active ingredient of CS, stigmasterol, could bind stably to the key target SRC protein. DISCUSSION: This study not only elucidates the phyto-pharmacological basis of CS in DC management but also provides a framework for developing natural product-derived targeted therapies against diabetic ocular complications. The integration of modern genomics and computational chemistry to deconstruct the therapeutic effects of traditional Chinese herbal medicines has great clinical significance in expanding the scope of traditional Chinese medicines for the treatment of DC and promoting precision targeting. However, this requires verification through basic experiments. CONCLUSION: These computational findings suggest that CS may exert its anti-cataract effects through the multi-target modulation of diabetic metabolic pathways and SRC-mediated signaling cascades.

Humans↗

PubMed2026

Altered EEG microstate dynamics reflect depressive symptoms in temporal lobe epilepsy.

BACKGROUND: Depressive symptoms are a common and disabling comorbidity in temporal lobe epilepsy (TLE), yet the neural mechanisms linking seizure networks to affective symptoms remain unclear. Although limbic network dysfunction has been implicated in both epilepsy and depressive disorders, it is unknown whether the time-varying dynamics of large-scale electrophysiological brain states reflect depressive symptom severity in TLE. In this study, we examined whether EEG microstate dynamics capture network alterations associated with depressive symptoms in individuals with unilateral TLE. METHODS: We analyzed resting-state, visually normal scalp EEG from 26 individuals with unilateral TLE. EEG microstates were identified by clustering global field power peaks into four canonical classes, with electrode positions mirrored to align the ictal hemisphere across subjects. Microstate dwell time, fractional occupancy, global transition entropy, and Markov transition probabilities were quantified and related to Beck Depression Inventory-II (BDI) scores. RESULTS: Individuals with high depressive symptoms (BDI ≥ 13; N = 12) exhibited longer mean dwell time in the ictal hemisphere-aligned microstate compared with individuals with low depressive symptom burden (BDI < 13; N = 14). Across subjects, dwell time in this microstate correlated with depressive symptom severity (r = 0.57, p = 0.002). TLE individuals with higher depressive symptoms exhibited reduced global transition entropy (p = 0.02), which also correlated with depressive symptom severity (r = -0.54, p = 0.004), indicating decreased flexibility of microstate transitions. Despite similar fractional occupancy of this state between groups, individuals with higher depressive symptoms were less likely to transition into the ictal hemisphere-aligned microstate from non-ictal or posterior configurations. Once engaged, however, the ictal-aligned microstate showed increased persistence, indicating prolonged stabilization of this network configuration. CONCLUSION: Higher depressive symptom burden in unilateral TLE is associated with increased temporal rigidity of the ictal hemisphere-aligned brain microstate, reflecting impaired disengagement of epileptogenic network configurations. These findings suggest that depressive symptoms in TLE may be associated with epilepsy-related disruptions in large-scale neural dynamics.

Humans↗

PubMed2026

Exploring the mechanism of Shengmai San in treating lung adenocarcinoma based on bioinformatics and molecular dynamics simulation.

To investigate the mechanism of Shengmai San (SMS) in the treatment of lung adenocarcinoma (LUAD) based on an integrated strategy combining "network pharmacology, bioinformatics, molecular docking, and molecular dynamics simulation," aiming to provide a precise combination therapy strategy and identify potential bioactive compounds. Differentially expressed genes in LUAD were identified from the Gene Expression Omnibus database using R (originally developed at Bell Laboratories and currently managed by Lucent Technologies). SMS components (ginseng, Ophiopogon japonicus, and Schisandra chinensis) were retrieved from encyclopaedia of traditional Chinese medicine, with Lipinski-compliant compounds selected. Compound targets were predicted via SwissTargetPrediction and Similarity Ensemble Approach. Intersecting targets between differentially expressed genes and compound targets were identified for "herbs-compounds-targets-disease" network construction. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analyses were performed. Hub targets were identified by analyzing the protein-protein interaction network. High-prognostic relevance targets were screened from The Cancer Genome Atlas. Compounds targeting these were identified through the herbs-compounds-targets-disease network, and absorption, distribution, metabolism, excretion, and toxicity-compliant compounds were selected using SwissADME (a web-based tool provided by the Molecular Modeling Group of the Swiss Institute of Bioinformatics). Core regulatory targets were identified through molecular docking, with complex stability assessed by molecular dynamics simulations. The key bioactive compounds of SMS for treating LUAD were identified as 7-hydroxy-2,5-dimethyl-4H-1-benzopyran-4-one, N-trans-feruloyltyramine, paprazine, and (E)-N-[(2S)-2-hydroxy-2-(4-hydroxyphenyl)ethyl]-3-(4-hydroxyphenyl)prop-2-enamide. Hub targets included AURKA, CCNA2, CCNB1, CDK1, CHEK1, KIF11, NEK2, PLK1, TTK, and TYMS. Among these, CDK1, CHEK1, and PLK1 demonstrated both high-prognostic relevance and strong binding affinity with SMS, emerging as core regulatory targets for SMS in LUAD treatment. Mechanistically, SMS exerts its anticancer effects primarily by modulating the tumor necrosis factor, interleukin-17, cell cycle, and Lipid and atherosclerosis signaling pathways. The active components of SMS, such as paprazine, may exert antitumor effects partly through downregulating CDK1, CHEK1, and PLK1 expression. Although the present study did not examine drug-resistance models or combination regimens, our findings raise the possibility that, in patients with high expression of these genes, combining SMS with standard chemotherapy or targeted therapy could potentially enhance chemosensitivity and mitigate the development of resistance. This hypothesis, however, requires formal testing in appropriate preclinical models and functional validation studies.

Molecular Dynamics Simulation↗

PubMed2026

Dynamics of antibiotic resistance genes co-occurrence with pathogenic and non-pathogenic bacteria throughout wastewater treatment processes.

Wastewater treatment plants (WWTPs) are recognized hotspots for antibiotic resistance genes (ARGs) and pathogenic bacteria. Despite advancements in treatment technologies, the persistence of ARGs and pathogenic bacteria remains a concern. In this study, we analyzed the dynamic changes in ARGs and bacterial communities throughout the treatment processes within an anaerobic-anoxic-oxic (AAO) WWTP over one week by using HT-qPCR coupled with 16S rRNA gene amplicon sequencing. The connectedness index, based on network analysis, showed that the dynamics of ARGs and mobile genetic elements (MGEs) were more strongly associated with potentially pathogenic bacteria than with non-pathogenic bacteria, suggesting that ARG immigration and dissemination in the WWTP were likely driven by potentially pathogenic taxa. The AAO treatment significantly reduced ARGs in final effluent (EF) (∼64 %) and residual sludge (RS) (∼81 %); however, potential hosts of ARGs such as Comamonas testosteroni and Clostridioides difficile persisted with minimal changes in relative abundance and remained detectable in EF and RS. Notably, the abundance of ARGs was lower in RS than in EF, and source tracking analysis identified influent as the primary source of ARGs and potentially pathogenic taxa in EF, underscoring the greater health risks associated with effluent discharge.

Wastewater↗

PubMed2025

Genome-Wide Analysis of DtxR and HrrA Regulons Reveals Novel Targets and a High Level of Interconnectivity Between Iron and Heme Regulatory Networks in Corynebacterium glutamicum.

Iron is vital for most organisms, serving as a cofactor in enzymes, regulatory proteins, and respiratory cytochromes. In Corynebacterium glutamicum , iron and heme homeostasis are tightly interconnected and controlled by the global regulators DtxR and HrrA. While DtxR senses intracellular Fe2+, HrrSA is activated by heme. This study provides the first genome-wide analysis of DtxR and HrrA binding dynamics under varying iron and heme conditions using chromatin affinity purification and sequencing (ChAP-Seq). We revealed 25 novel DtxR targets and 210 previously unrecognized HrrA targets. Among these, metH, encoding homocysteine methyltransferase, and xerC, encoding a tyrosine recombinase, were bound by DtxR exclusively under heme conditions, underscoring condition-dependent variation. Activation of metH by DtxR links iron metabolism to methionine synthesis, potentially relevant for the mitigation of oxidative stress. Beyond novel targets, 16 shared targets between DtxR and HrrA, some with overlapping operator sequences, highlight their interconnected regulons. Strikingly, we demonstrate the significance of weak ChAP-Seq peaks that are often disregarded in global approaches, but feature an impact of the regulator on differential gene expression. These findings emphasize the importance of genome-wide profiling under different conditions to uncover novel targets and shed light on the complexity and dynamic nature of bacterial regulatory networks.

Corynebacterium glutamicum↗

PubMed2026

NLCD: A method to discover nonlinear causal relations among genes.

Distinguishing correlation from causation is a fundamental challenge in many scientific fields, including biology, especially when interventions like randomized controlled trials are infeasible and only observational data are available. Methods based on statistical tests of conditional independence within the Mendelian Randomization framework can detect causality between two observed variables that are each associated with a third instrumental variable. However, these methods for detecting causal relationships between traits (e.g., two gene expression or clinical traits associated with a genetic variant, all observed in the same population) often assume a linear relationship, thereby hindering the discovery of causal gene networks from genomics data. We have developed NLCD, a method for NonLinear Causal Discovery from genomics data based on nonlinear regression modeling and conditional feature importance scoring. NLCD uses these techniques to extend the statistical tests in an existing linear causal discovery method called the Causal Inference Test (CIT). We benchmarked NLCD against current state-of-the-art methods: CIT, Findr, and MRPC. On simulated datasets, NLCD performs comparably to most methods in detecting linear relations (Average AUPRC (Area Under the Precision-Recall Curve) of NLCD = 0.94, CIT = 0.94, Findr = 0.94, and MRPC = 0.99), and outperforms them in detecting nonlinear (sine and sawtooth type) relations between two genes (Average AUPRC of NLCD = 0.76, CIT = 0.60, Findr = 0.56, and MRPC = 0.73). When tested on a nonlinear subset of a yeast genomic dataset to recover known causal relations involving transcription factors, NLCD and CIT performed comparable to each other and slightly better than Findr and MRPC (Average AUPRC of NLCD = 0.82, CIT = 0.81, Findr = 0.71, and MRPC = 0.54). On application to a human genomic dataset, NLCD revealed active causal gene pairs (IRF1 → PSME1 and HLA-C → HLA-T) in the muscle tissue, and clarified the promises and challenges in discovering causal gene networks in tissues under in vivo human settings.

Humans↗