2026-08-12 | | Total: 18
Predicting how molecular changes affect large-scale brain activity is a difficult task because of the lack of appropriate methods to link scales. In this perspective, we review a class of mean-field models that can integrate biophysical details such as synaptic receptors or membrane ion channels. This leads to a multi-scale modeling approach that can be used to evaluate how microscopic changes can impact macroscopic brain activity. This approach is illustrated here for the case of anesthesia, where changes at the level of specific synaptic receptors can lead to a global change in brain activity and a disconnection from external inputs. This is only possible using mean-field models that can include enough detail about the microscopic biophysical properties. This biophysically-based mean-field approach could be generalized to study cellular or molecular origins of brain diseases, or to better understand how drugs acting at microscopic scales can influence global brain activity. Biophysical mean-field models also link different fields of neuroscience, from molecular studies to brain imaging.
Reproductive hormones regulate a wide range of physiological processes throughout the human lifespan. Estrogen, in particular, varies substantially across the menstrual cycle and is widely used in contraceptives and hormone replacement therapies. Despite its physiological importance, few experimental studies and even fewer mathematical models explicitly investigate how estrogen regulates smooth muscle function. As smooth muscle lines our blood vessels, airways, uterus, and several other organs, understanding how estrogen impacts smooth muscle is important to improving the understanding of sex differences in lifelong health. Here we extend an established mathematical model of smooth muscle cell calcium signalling to incorporate estrogen-dependent modulation of intracellular calcium transport pathways. Numerical simulations, global sensitivity analysis, and numerical bifurcation analysis are then used to quantify the influence of estrogen on intracellular calcium dynamics and the resulting steady-state and oscillatory behaviours. Our results demonstrate that physiologically relevant changes in estrogen shift intracellular calcium concentrations while leaving the underlying bifurcation structure and qualitative dynamics largely unchanged, suggesting that estrogen acts primarily as a quantitative modulator of smooth muscle calcium signalling. This work also serves to establish a foundation for future mechanistic models of hormone-dependent cell physiology.
Neural tracking - the time-locking of neural responses to continuous stimuli such as speech, music, and video - is widely used to study how the brain processes natural input. Tracking strength is typically quantified as the correlation between the recorded neural response and the stimulus, decoded and/or encoded through data-driven models, and this correlation is routinely used to compare stimulus features, models, or settings. However, its magnitude depends not only on how strongly the brain tracks the stimulus, but also on the statistical properties of the signals being correlated. For example, a smallband speech envelope carrying almost no information about speech content yields among the highest correlations, simply because it is easier to reconstruct. Meaningful interpretation therefore requires comparing each correlation to its null distribution: the correlations expected without any stimulus-response relationship. We show that the randomization procedures commonly used to construct this null distribution are not interchangeable: each implicitly encodes a different null hypothesis, and we motivate stimulus-response misalignment as the most practical and appropriate choice. Because reliable null distributions require many permutations, we introduce a semi-parametric model using the normal distribution after the Fisher transform that yields accurate significance levels from only 3-5 min of data and predicts them across analysis window lengths. Building on this, we propose the null-normalized tracking score, an interpretable measure placing features and models on a common scale, which relates directly to the widely used match-mismatch accuracy. Applied to EEG from 121 participants listening to continuous speech, the framework reverses conclusions drawn from raw correlations, providing an efficient and principled methodology for interpreting neural tracking correlations.
A key issue in using AI for pathology diagnosis is what image information should be given to the AI and how limited analysis resources should be used. This study compares two ways of processing different types of images under limited resources. The first is WSI-AI, in which AI automatically compresses information from the whole slide image (WSI). The second is Diagnostic Field Selection (DFS)-AI in which an expert first selects several regions, magnifications, and comparisons needed for diagnosis, and the AI then analyzes those selected fields. To compare these two types, we built three image models: a rarely localized lesion, lesion detection in a nonuniform background, and a spatially continuous lesion. As a result, WSI-AI was better when a coarse view did not give enough information about lesion location in advance, as may occur with very small tumor foci in lymph nodes. In contrast, when a low-cost coarse view provided useful location information, DFS-AI was better in an intermediate range of limited resources. When searching for a specific object such as infectious organism in a nonuniform background, the relative value of DFS-AI increased. For spatially continuous lesions, WSI-AI found the presence of a lesion more easily as the lesion became larger. However, for complete characterization of the entire lesion, DFS-AI could be better because it used the continuous structure to select fields for analysis. In conclusion, both relative performance changes in a systematic way with location information from the coarse view, background heterogeneity, local context, and the spatial structure of the lesion. A practical strategy under limited resources is for physicians to choose between WSI-AI and DIS-AI according to the information structure of the diagnostic task.
Protein backbone generation models are often credited with exploring novel fold space based solely on low full-chain similarity to known proteins, yet this cannot distinguish a genuinely new fold from a novel assembly of known structural units. We first ask whether this granularity mismatch alone explains the reported rates, and introduce the Domain Retrieval Rate (DRR), the fraction of generated backbones for which any constituent domain matches a known domain in CATH S40. Applied to eight backbone generation models spanning diffusion and flow-matching paradigms, DRR finds locally alignable known structure in most outputs, while the fraction containing a substantially covered complete domain is considerably smaller and depends on the scoring convention. To calibrate what retrieval alone can achieve, we propose RetFold, a zero-training baseline that constructs backbones by retrieving CATH domains and refining inter-domain connections through geometry-based helix-linker optimization, at two orders of magnitude lower cost on CPU alone.
Proteolysis-targeting chimeras (PROTACs) induce protein degradation by recruiting a target protein to an E3 ubiquitin ligase, making degradation a joint outcome of the degrader molecule and its biological context. Although public databases contain thousands of structured molecule-target-E3 records, degradation measurements are available for only a small fraction of them. Existing supervised approaches therefore leave most recorded chemical-biological relationships unused. We introduce DegradeQuery, a context-aware prediction framework that converts these label-missing records into a pretraining signal. Its counterfactual tuple pretraining objective contrasts recorded tuples with alternatives formed by replacing the target, the E3 ligase, or both, enabling the model to learn contextual associations without assigning activity pseudo-labels. The resulting representation is then fine-tuned to predict degradation from the complete molecule-target-E3 context. On the official PROTAC-8K benchmark, DegradeQuery achieves an area under the receiver operating characteristic curve of 0.9065 and an accuracy of 0.8500, outperforming the compared methods. Controlled analyses further show that the improvement is primarily attributable to tuple-level pretraining, can be recovered using only label-missing records, and remains complementary to protein language model representations. These findings demonstrate that incompletely labeled PROTAC databases contain useful relational supervision and provide a practical route for learning context-aware degradation predictors from scarce experimental labels.
How many distinguishable labels can a biological oscillator carry? Proposals invoking collective vibrational modes, endogenous electromagnetic fields, microtubule excitations and oscillatory phase codes are each debated on grounds particular to themselves, with no shared standard for comparison. We show that spectral distinguishability alone bounds the number of labels by the quality factor, M <= Q = 2 pi nu tau. This follows from the relation between linewidth and coherence time, so it is independent of substrate, of mechanism, and of any position on quantum effects in biology, and it can be evaluated from two published quantities. Applied to a recently proposed 30 GHz intracolumnar microwave field in cortex, it gives Q = 0.19: the linewidth exceeds the carrier five-fold. The obvious rescue, that a driven emitter can be spectrally narrower than its gain medium, requires a resonant cavity, and the model's own geometry forbids one. An independent bound on metabolic power is exceeded by five to nine orders of magnitude. Six further criteria follow from the same standpoint, including a two-sided persistence window requiring a label to be both readable and rewritable. Screening eleven carriers, only the low-frequency neural rhythms pass. High-frequency molecular carriers are eliminated by brevity, not by the fragility the debate has assumed.
Aquaporins (AQPs) are an ancient group of channels that arose early in phylogeny and radiated throughout the kingdoms of life. Classified as members of the broad Membrane Intrinsic Protein family, AQPs originally were envisioned as strictly water pores that were constitutively open. International research interest is driving new views of AQPs as multi-functional channels with complex control mechanisms that enable adaptive responses to physiological challenges. Expanding numbers of AQP classes are being discovered to allow permeation of diverse solutes in parallel to the archetypal substrate water, with channel activities that show subtype-specific regulation by intracellular and extracellular signals. To date, at least fifteen classes of AQPs from mammals, plants, insects and algae have been proposed to carry ions, enabling amino acid sequence comparisons that suggest ion channel functional domains arise from conserved transmembrane, loop and teminal regions of the tetrameric AQP protein, serving as selectivity filters, gates, and sites for modulation. More AQP ion channel classes likely await discovery, pending identification of relevant activators. Subtype-specific control of dual water-and-ion AQP classes orchestrate diverse roles for AQPs across the domains of life, promoting homeostasis, cell motility, nutrient acquisition, redox protection, sensory detection, and intercellular signaling in cells and tissues.
Characterizing the dissimilarity of neural representations between experimental conditions, and tracking it across time, is a central goal of multivariate pattern analysis. Guggenmos et al. (2018) assessed the reliability of many of the measures that can be used for that purpose on MEG data and recommended the use of either the cross-validated Euclidean distance or the within-class-corrected Pearson distance. In this commentary, we show that we can improve upon these distances. First, we show that the cross-validated Euclidean distance is equivalent to a sum of between-partition distances and that this equivalence can be leveraged to obtain a generalized variant, with increased reliability and accuracy. Second, we use the relationship between Euclidean distance and Pearson correlation to define a cross-validated correlation distance in a similar way. The resulting distance is more accurate and interpretable than a formulation proposed by Guggenmos and colleagues. Finally, we discuss the relationship between our generalized cross-validation and within-class correction, another strategy often used to increase reliability, and we show that generalized cross-validation results in higher accuracy for the correlation distance.
Computational head models are typically applied to isolated impacts, leaving repeated head loading largely unexplored. An Ogden-Roxburgh Mullins damage formulation was implemented in a high-fidelity finite element head model to represent loading-history-dependent softening during cyclic brain-tissue deformation. Repeated-loading histories derived from mixed martial arts head-impact data were applied and compared with damage-free hyperelastic (HE) and linear visco-hyperelastic (LVHE) model variants. Under five identical single-axis cycles, Mullins-type softening progressively increased strain and strain rate metrics relative to the HE model. Mullins-based injury probabilities progressively exceeded strain-based HE predictions and diverged from unchanged kinematics-based predictions, indicating that neglecting prior softening may underestimate injury risk. In a randomized twenty-cycle multiaxial sequence, cycles of similar kinematic intensity produced different deformation and injury-risk estimates depending on prior softening. HE and LVHE models predicted higher injury probabilities initially, whereas the Mullins-based model produced the largest later-cycle estimates and highest probability of at least one injury over the sequence. Regional amplification depended on loading direction and prior softening, with no direction-independent trend among brain substructures. Gyral elements exhibited higher cumulative maximum principal strain than sulcal elements, which showed greater amplification relative to initial responses. These findings demonstrate that short-term damage-mediated softening can substantially amplify tissue deformation and injury-risk estimates beyond damage-free head models under the same loading histories. Further experimental characterization of cyclic brain-tissue softening is needed to improve models of repeated head loading and traumatic brain injury.
Background: Brain information processing deteriorates as cortical neurons gradually transition from a rested state into fatigue. Subjectively, fatigue is experienced as a state of weariness, tiredness, or lack of energy that reduces the ability to work safely and effectively. Objective: In this theoretical paper, we pinpoint the physical origin of brain fatigue in the gradual deterioration of biochemical reaction quotients and transmembrane ion concentration gradients, which increase neuronal excitability and decrease the signal-to-noise ratio in the brain cortex. Methods: Brain performance in a rested state versus fatigue is examined by well-established, data-driven computer models for energy transport inside protein $α$-helices in the presence of thermal noise for different ATP energy states or for pyramidal neuron firing of action potentials under electric stimulation in a rested membrane state versus fatigue. Results: We found that reduced Gibbs free energy supply from ATP hydrolysis impairs the cooperative effect between amide I excitons propagating inside protein $α$-helices, with resulting decreased thermal stability of molecular solitons. Concurrent changes in Nernst reversal potentials for Na+ or K+ ions further led to neuronal hyperexcitability and a higher risk of neuronal depolarization block during mental fatigue. Conclusions: Detailed computational modeling showed that inefficient protein function due to diminished ATP energy status, alters the electrophysiological properties of individual neurons, thereby impairing their information processing capacity for the proper execution of cognitive tasks. Scheduling practices aimed at intermittent recovery of the rested brain state during intellectually challenging work could protect physical and mental wellbeing, prevent burnout, and enhance long-term productivity.
Many computational problems possess the following three features: (1) the problem could be solved efficiently if a linear operator could be diagonalized, (2) direct diagonalization is infeasible due to scale, but, (3) the operator is known to commute with a permutation since the problem is (a) symmetric with respect to a change of coordinates, or (b) the operator is block circulant. For example, high-resolution template matching demands repeated convolution of large, high resolution images with an exhaustive list of related templates. The associated linear operator could be compressed, via a low rank approximation, if diagonalized. The required decomposition expensive, but, the entire problem is symmetric to in-plane rotations. In these cases, the eigenfunctions are constrained by the symmetry. These constraints allow efficient decomposition. We illustrate a parallelized algorithm that allows fast diagonalization of any block-circulant matrix. When the index space can be partitioned into $l$ classes of $m$ interchangeable elements, the algorithm reduces storage costs by a factor of $m$ and, if $w$ workers are available, the algorithm demands $\mathcal{O}(l^2m\log(m)/w) + \mathcal{O}(ml^3/w)$ floating point operations per worker yielding a $m^2$ speedup. We use this procedure to decompose a high precision template matching matrix. We compare runtime on a reduced problem, where the fast procedure ran 205 times faster per recovered feature, for 22.5 times as many features. On a typical template matching matrix, the decomposition required 25 times less time than computing the matrix entries. We also demonstrate decomposition of a $\sim$30 times larger matrix that covers the complete set of possible projections represented in a cryo-EM image at 2~Å resolution in 14 minutes. This procedure is more stable and allows recovery of in-plane rotations to arbitrary precision.
Tumor progression is an inherently multiphysical phenomenon in which interstitial fluid dynamics, biochemical transport, and cellular mechanics interact across multiple spatiotemporal scales. Classical mesh-based solvers, although accurate, impose prohibitive computational costs for the repeated evaluations demanded by inverse parameter identification and future patient-specific predictive pipelines. In this work we introduce a Physics-Informed Neural Network (PINN) framework for a tractable chemo-fluidic continuum model of tumor growth that couples an advection-diffusion-reaction (ADR) equation for the tumor volume fraction with a quasi-static Darcy pressure equation for the interstitial fluid pressure. By intentionally decoupling the solid-mechanical equilibrium, we obtain a three-equation system whose gradient structure is stable under automatic differentiation, enabling robust deep-learning optimization. The network simultaneously learns both state variables from physics constraints alone (forward problem) and recovers hidden transport parameters from sparse, noisy synthetic measurements (Data-Assimilation PINN, DA-PINN, inverse problem). We verify the forward solver against a high-resolution finite-difference (FD) reference, achieving a mean absolute error below 0.002. For the inverse problem, starting from an initial permeability estimate of 0.08 (a factor of 4x above the true value of 0.02) with only 5% spatially sparse observations corrupted by 5% Gaussian noise, the DA-PINN recovers the permeability with a relative error below 5%. These results demonstrate that physics-informed deep learning constitutes a viable, computationally efficient route to multiphysics oncology modeling and lays the mathematical groundwork for future integration into clinical data assimilation pipelines.
Continuous hemodynamic monitoring guides treatment decisions in surgery and intensive care. However, gold-standard signals are only measured in severe cases due to risks associated with invasive measurement. In this work, we introduce HIPNO (Hemodynamic Inference via Physics-informed Neural Operators) to recover hemodynamic state from ubiquitous, non-invasive signals and expand access to advanced monitoring. HIPNO addresses a problem of scale symmetry in physics-informed hemodynamic inference, where different combinations of flow, resistance, and compliance can generate the same observed pressure. We identify the symmetry group of the observation model and parameterize the network in its quotient space. For the 3-element Windkessel model, the quotient coordinates are the compliance-normalized flow $U=Q/C$, the decay time constant $τ_{WK}=R_2 C$, and the characteristic-impedance coordinate $κ=R_1 C$. Across 945499 intraoperative windows from 2562 patients, HIPNO predicts $τ_{wave}$, a proxy for vascular decay derived from pressure, with 32% lower error on the log scale than a population baseline while preserving mean arterial pressure accuracy. Because vascular decay and flow drive occupy separate coordinates, counterfactual perturbations produce the expected directional responses in at least 90% of windows in almost all prespecified scenarios, a separation unavailable to pressure-only baselines. The coordinates are also used as inputs to a calibration model for monitored cardiac output. Finally, the formulation identifies the external compliance or flow reference required to recover absolute physical scale.
Yellowfin tuna (Thunnus albacares) is an economically important pelagic species in the Indian Ocean, but its monthly catch variability is difficult to predict because reported catch is influenced by ocean productivity, thermal habitat, monsoon seasonality, fishing-ground characteristics, and fleet behavior. This study develops a remote sensing-based statistical and machine learning stacked ensemble for modeling monthly yellowfin tuna catch across five fishing grounds in the western equatorial Indian Ocean during 2003-2024. Monthly fishery records were integrated with satellite-derived chlorophyll-a concentration (CHL) and sea surface temperature (SST). Linear mixed models (LMMs) and a generalized additive mixed model (GAMM) were used to represent interpretable current and lagged environmental effects, while Random Forest models captured nonlinear interactions among lagged, rolling, seasonal, anomaly, and spatial predictors. Expanding-window validation over the pre-COVID evaluation period from 2008 to 2019 was used to reduce temporal leakage. Among the individual models, the GAMM achieved the best log-scale performance, with an R-squared value of 0.137 and an RMSE of 1.284, while RF-E6 achieved the lowest RMSE on the kilogram scale. The nonnegative Ridge stacked ensemble improved the log-scale R-squared value to 0.177, with an RMSE of 1.254 and an MAE of 1.000. These results show that satellite-derived CHL and SST contain useful but incomplete predictive information for monthly catch modeling.
Active-quiescent switching occurs in biological populations in which growth is confined to a proliferative active state, while cells may reversibly enter a nonproliferative quiescent state. Experiments often observe only part of this process, through active-state markers, aggregate population measurements, or aggregate data supplemented by sparse active-state observations. Such measurements arise when marker panels are limited, active-state assays require fixation or endpoint sampling, or only total cell number, optical density, tumor burden, or aggregate fluorescence is reported. These observation choices complicate sparse regression methods such as sparse identification of nonlinear dynamics (SINDy), because the measured variable may satisfy a different equation from the underlying active-quiescent system. A library chosen for the wrong observable may therefore fit a trajectory without preserving mechanistic interpretation or transferability. We study this issue using a two-compartment ordinary differential equation model. We define observable reduction as the elimination of hidden states to obtain the differential equation satisfied by the measured variable. For several biologically relevant growth laws, we derive observation-specific reductions and use them to construct sparse-regression libraries. Using synthetic data, we compare these structured libraries with standard polynomial SINDy. Polynomial libraries can match training trajectories while failing coefficient-relation and transfer tests, whereas reduction-guided libraries recover interpretable coefficient maps when the observed regime is informative. These results show that interpretable equation learning in hidden-compartment systems requires matching both the regression target and candidate library to the observation process.
A subcritical age-structured branching population dies out almost surely. Conditioned on survival, it converges to its Yaglom limit, a quasi-stationary equilibrium that we take as the operating point for control. We model preventive removal (culling) as an age-dependent actuator that raises the mortality rate and leaves the offspring law untouched, and we show that its authority over this equilibrium is bounded for structural reasons. Two facts drive the result. First, the input is matched to the killing rate but unmatched with respect to the Foster--Lyapunov drift, so the transmission barrier $\Lbar$ set by reproduction alone is invariant under such actuation. Second, and this does not follow from invariance alone,} the supremum of the reachable decay rates is $\Lbar+ν^\star$, where $ν^\star\le0$ is the Malthusian parameter of the lineage conditioned never to die childless; the gap $|ν^\star|$ is given in closed form and vanishes exactly when no individual has two or more offspring. Consequently no removal law of this class reaches the barrier, and along the admissibility boundary the achievable decay rate is governed by the shape of the actuator rather than by its size. We illustrate these results on a model calibrated to the 2001 Cumbrian foot-and-mouth outbreak.
We present a model-agnostic approach to the validation of computational models, and a python library containing a suite of validations for mouse primary visual cortex models. By viewing the model not as a stand-in for the target system, but instead as a generator of predicted system properties, this approach allows validations to easily generalize to many different models, while also making them comparatively simple to write. Several limitations of this approach are outlined alongside possible means of addressing them. This presents an important step in addressing what we consider the most important problem in computational science: the standardization of model validation.