IFAC-PapersOnLine

The dynamics of gene regulatory networks are governed by intrinsic noise, stemming from the random nature of biochemical reactions, and by extrinsic noise, arising from fluctuations in cellular components and environmental conditions. Together, these sources can compromise the reliability of predictive computational models if not properly accounted for, and capturing both effects within a single framework remains a non-trivial task in computational biology. In this work, we propose an uncertainty quantification framework that addresses these two contributions jointly: intrinsic stochasticity is described through a partial integro-differential equation (PIDE) for the protein probability density function, whereas extrinsic noise is represented as parametric uncertainty in the kinetic parameters. The propagation of the uncertainty is carried out via an intrusive polynomial chaos expansion (PCE), in which the PCE coefficients are obtained from a stochastic Galerkin projection of the PIDE, yielding a coupled deterministic system that is solved with standard numerical methods. We illustrate the approach on a positive autoregulatory gene network with one and two uncertain kinetic parameters. The proposed approach accurately reproduces the mean, variance, and full protein probability density function, including the bimodal distributions, at a substantially lower computational cost.

Planning husbandry tasks and experiments with Drosophila melanogaster requires converting a target date into development times that depend on the rearing temperature. This calculation needs to be done for each cross, genotype, and temperature, and the risk of error grows quickly. Available laboratory management tools let users register stocks, crosses, and track them, but they do not create schedules based on a clear, adjustable thermal model. To fill that gap, we developed DrosoTracker, a self-contained web application that works offline and predicts Drosophila development with a thermal summation model recalibrated through regression on data from Powsner (1935) (T0 = 11.78 {degrees}C, DD = 116.38 {degrees}C{middle dot}days, R{superscript 2} = 0.997). The model offers an optional two-level calibration driven by user observations. A wild-type strain first adjusts the model to the laboratorys own conditions. Then each genotype is calibrated against that reference using a random-effects shrinkage estimator that accounts for measurement error and between-batch variability. The model creates schedules for husbandry tasks, evaluates adult cohort survival with the Kaplan-Meier estimator and the log-rank test, and calculates sample size for lifespan studies using Schoenfelds formula. The quantitative components were checked against independent references, including Rs survival package and manual calculations. Ongoing work is focused on validating the calibrated model using cohorts specifically bred for this purpose. DrosoTracker runs entirely in the browser, stores data locally, and is available in English and Spanish.

Microbial fuel cells (MFCs) have been widely investigated as decentralized bioelectrochemical systems capable of converting organic substrates into electricity. However, their long-term autonomous operation is constrained by substrate depletion in the anode compartment, leading to metabolic starvation of electroactive biofilms and a decline in power output. Conventional MFC design treats substrate crossover through the membrane separator as a parasitic loss that reduces coulombic efficiency. In this work, we propose a conceptual inversion of this paradigm by considering controlled cathodic-to-anodic substrate crossover as a passive mechanism to sustain basal microbial metabolism during periods of substrate scarcity. A transport-reaction framework is developed to quantify the balance between membrane-mediated substrate flux and microbial maintenance demand within the anode biofilm. Based on this balance, a dimensionless maintenance crossover Damkohler number (Dam) is introduced to define three operational regimes: starvation-dominated (Dam >> 1), balanced autonomous (Dam {approx} 1), and crossover-dominated (Dam << 1). The framework integrates membrane transport theory with biofilm kinetics to evaluate the effects of separator properties, substrate gradients, and current-dependent electro-osmotic transport on system stability. Order-of-magnitude analysis indicates that achievable crossover fluxes span several orders of magnitude depending on separator characteristics, suggesting that membrane properties critically influence system behavior. This perspective reframes substrate crossover from a loss mechanism to a potential design variable, offering a conceptual tool for enhancing resilience and guiding separator selection in MFCs intended for long-duration, and low-maintenance operation. HighlightsO_LIControlled crossover can sustain microbial metabolism in MFCs C_LIO_LIIntroduces maintenance crossover Damkohler number (Dam) C_LIO_LIIdentifies regimes for autonomous and starvation operation C_LIO_LILinks membrane properties to long-term system stability C_LIO_LIReframes crossover as a design variable, not only a loss C_LI

High-dimensional population balance equations (PBEs) provide a natural framework for modeling heterogeneous cell populations, but their direct numerical solution becomes computationally prohibitive when the internal state space contains many molecular variables. We propose a hybrid mechanistic-machine learning framework for reducing and simulating PBEs defined over high-dimensional intracellular coordinates. The cell population is described by a number density n(x, t), where x [isin] [R]N represents gene and protein states associated with macrophage activation. A dynamics-preserving autoencoder maps this state space to a low-dimensional latent coordinate z [isin] [R]d, with d << N, while retaining key qualitative features of the underlying gene regulatory network, including attractor structure and multistability. Mechanistic information from the original regulatory dynamics is used to construct interpretable drift and diffusion terms for the reduced latent-space PBE. The reduced PBE is solved using a stochastic Lagrangian particle representation, in which particles evolve according to stochastic differential equations (SDEs) corresponding to the latent drift and diffusion fields. The resulting latent-space solution is subsequently decoded and propagated back into the original state space to recover physically interpretable cellular dynamics. We demonstrate the framework on macrophage polarization under cytokine-dependent regulation, including gene knockout perturbations. Overall, the proposed framework provides a computationally tractable and mechanistically interpretable route for integrating single-cell genomic data with population balance models of cell-state dynamics.

Tumors consist of heterogeneous phenotypic cells, such as normoxic cells, which are highly proliferative, and hypoxic cells, which are less proliferative. Their phenotypic switching depends on tumor microenvironmental factors, such as oxygen and nutrient concentrations supplied by local blood vessels. However, during ongoing angiogenesis, the process of sprouting new blood vessels at the tumor site from pre-existing blood vessels, and how this phenotypic switching affects and impacts tumor growth, remains poorly understood. In this article, we formulate a mathematical model to elucidate the crosstalk between vasculature and tumor cellular heterogeneity during tumor progression. The model results show a strong agreement with the experimental data. Our simulation results demonstrate that ongoing angiogenesis increases tumor growth rate. In addition, we observe that the influence of hypoxic cells on phenotypic switching from normoxic to hypoxic is more pronounced than their influence on the transition from hypoxic to normoxic. Furthermore, we perform a global sensitivity analysis using the Sobol's method to assess the importance of the model's parameters. It highlights that the volume at which blood vessels attain half-maximal rate has the maximum effect on the model.

Multi-stage continuous bioprocessing can increase volumetric productivity, operational consistency, and process throughput, but its design is complicated by coupling among dilution rate, reactor volume, feed allocation, and cellular physiology. Here, we present ContiDesigner, available at https://chemnettools.anc.univie.ac.at/ContiDesigner/, a mechanistic steady-state framework and interactive web tool for the system-level design of continuous fermentation cascades. Comparing one- and two-stage configurations at equal total reactor volume and outlet flow, ContiDesigner reveals how internal flow and reactor volume allocation shape space-time yield and identifies productivity-maximizing operating conditions. Compared with one-stage processes, two-stage cascades favor lower over-all dilution rates, thereby preserving residence time in the production stage. The first-stage dilution rate approaches the corresponding one-stage productivity optimum, but the cascade optimum occurs earlier, reflecting a system-level compromise between biomass generation and production-stage residence time. However, two-stage operation outperforms optimized one-stage operation only when non-growth-associated production in the second stage is sufficiently strong, whereas increasing growth coupling favors one-stage operation. Two case studies demonstrate both the potential and limits of process intensification. An optimized two-stage design is predicted to achieve a more than 1.5 fold increase in space-time yield for poly-R-3-hydroxybutyrate (PHB) production compared with a published experimental five-stage cascade, whereas the lactic acid case study identifies conditions under which staging offers no advantage. ContiDesigner translates these design principles into an accessible workflow to explore feasible operating regions and prioritize cascade designs for experimental evaluation. Graphical Abstract O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=130 SRC="FIGDIR/small/743657v1_ufig1.gif" ALT="Figure 1"> View larger version (37K): org.highwire.dtl.DTLVardef@ef58faorg.highwire.dtl.DTLVardef@1ba88a4org.highwire.dtl.DTLVardef@160edd3org.highwire.dtl.DTLVardef@9dda34_HPS_FORMAT_FIGEXP M_FIG C_FIG O_LIContiDesigner enables system-level design of continuous fermentation cascades C_LIO_LIHigh stage-one dilution supports biomass generation C_LIO_LILow stage-two dilution preserves productive residence time C_LIO_LIYet two-stage cascades favor lower overall dilution than one-stage systems C_LIO_LITwo-stage advantage requires strong non-growth-associated production in stage two C_LI

Maintaining a prescribed composition in engineered microbial consortia is difficult because small fitness differences can drive competitive exclusion. We study a two-strain consortium in continuous culture and develop a feedback architecture that regulates composition by selectively slowing the fast strain as a function of the population ratio. At the population level, we derive an idealized ratio-feedback law with a tunable positive coexistence equilibrium. We then propose a biomolecular realization using orthogonal quorum sensing, an sRNA-based ratiometric controller, and a ppGpp-mediated growth actuator. Exploiting the separation between slow population growth and faster intracellular controller dynamics, we use singular perturbation theory to show that, for sufficiently fast controller dynamics, the full implementation model inherits the coexistence equilibrium and its local stability properties from the reduced model. Numerical simulations validate the reduction and show how weaker timescale separation or loss of the assumed molecular regime degrades performance.

Selecting stochastic gene-expression models from single-cell counts requires accurate parameter inference and efficient model selection. Likelihood methods in count space can be costly when full stationary count distributions are unavailable, whereas approximate methods may lose accuracy. Probability generating functions (PGFs) offer a compact analytical alternative, but existing PGF workflows are generally not likelihood based and therefore rely on computationally intensive cross-validation. We develop a likelihood-based PGF framework for both tasks. Correlated empirical PGF values are used to construct a Gaussian quasi-likelihood for parameter inference and PGF-based Bayesian information criterion (BIC) for model selection. We show that the empirical PGF is exactly unbiased and that the parameter estimator is consistent, converges at the inverse-square-root sample-size rate, and is first-order asymptotically unbiased. For large samples and a uniquely preferred model, PGF-BIC selects the same model as leave-one-out cross-validation in PGF space.

Physical human-device interaction is central to many emerging technologies in neurorehabilitation and assistive robotics, but simulation-based research in this area remains fragmented across musculoskeletal models, assistive-device representations, task definitions, and controller-development workflows. This fragmentation limits the accessibility, reproducibility, and extensibility of studies on prostheses, exoskeletons, wearable rehabilitation devices, and related human-device systems. Here we introduce MyoAssist 1.0, an open-source framework for neuromechanical simulation of physical human-device interaction built within the MyoSuite ecosystem. MyoAssist organizes each simulation environment as a composed human-device-task system that combines compatible musculoskeletal, assistive-device, and task-scenario components through a shared composition pipeline. The current release includes 15 assistive-device models spanning gait assistance, upper-body support, manipulation, and seated mobility and supports compatible musculoskeletal models ranging from reduced lower-limb models to a 416-muscle full-body model. These human-device systems can be simulated within the broad task scenarios provided by MyoSuite, while MyoAssist adds locomotion-specific task scenarios with configurable terrain and target-velocity conditions for gait-assistive studies. MyoAssist also provides two complementary controller-development frameworks: a reinforcement-learning framework for training adaptive policies and a controller-optimization framework for tuning structured, interpretable human and device controllers. Both frameworks operate on the same simulation environments and provide standardized evaluation outputs for inspecting, comparing, reusing, and extending learned and structured control strategies. By integrating modular human models, assistive-device models, task scenarios, and training workflows under a shared open-source interface, MyoAssist aims to lower the barrier to reproducible simulation-based research and to support collaborative development of assistive technologies for neurorehabilitation and physical human-device interaction.

Deep learning tools are increasingly used today, particularly in medical segmentation. A gap nonetheless remains in automating segmentation for insects. This work addresses the following question: can a generalist segmentation model, trained on several phylogenetically related orthopteran species, reliably automate head tissue segmentation from micro-CT images? To answer this, we used nnU-Net, a self-configuring 3D deep learning segmentation framework originally developed for medical imaging, whose core function, learning to recognize tissues of interest, applies directly to this context. Six anatomical classes were automated, comparing two training strategies: sequential fine-tuning, which adds species one at a time under the assumption that progressive learning would strengthen predictive power, and from-scratch training, in which the model learns the entire dataset simultaneously. The fine-tuning model (ModelB) reached a Dice coefficient (a measure of overlap between automated segmentation and manual ground truth, ranging from 0 to 1) of 0.7715, compared to 0.7664 for the from-scratch model (ModelC). Although both models produced accurate automated segmentations, no significant difference was found between the two training strategies (paired Wilcoxon test, n = 24, p = 0.243). Despite a dataset limited to 20 individuals and the absence of one method clearly outperforming the other, the models remain usable across the three species studied (Gryllus bimaculatus, Loxoblemmus equestris, L. doenitzi), including in the presence of pronounced sexual dimorphism. It reduces a 20 hour segmentation task to under a minute.

Isotopic Metabolic Flux Analysis (I-MFA) is a standard approach for estimating intracellular metabolic fluxes. I-MFA infers fluxes by comparing simulated and measured metabolite isotopologue distributions (MIDs) of metabolites from isotope labeling experiments. MIDs represent fractional abundances that strictly sum to one for any given metabolite, thus they are inherently compositional data. However, state-of-the-art estimation approaches rely on calculating standard Euclidean distances between MIDs in a non-compositional paradigm, introducing a systemic bias. To resolve this, our study proposes compositional I-MFA. We demonstrate how to construct a meaningful orthonormal basis for MIDs via ordered sequential binary partitioning, which can be used to perform isometric log-ratio (ILR) transformation. As a minimal change to existing I-MFA workflows, we suggest estimating fluxes by minimizing Euclidean distances between ILR-transformed MIDs. We validated this framework against traditional methods using both a toy model and a biologically realistic model, evaluating point estimates, sensitivity across varied true fluxes, and confidence intervals. In the two examples, compositional I-MFA consistently outperformed traditional approaches, reducing mean squared error of flux point estimates by an average of 42.6% and substantially narrowing confidence intervals. We conclude that compositional data analysis significantly improves I-MFA and can be implemented as a simple drop-in replacement for current pipelines. Graphical Abstract O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=156 SRC="FIGDIR/small/742769v1_ufig1.gif" ALT="Figure 1"> View larger version (34K): org.highwire.dtl.DTLVardef@1a53aa4org.highwire.dtl.DTLVardef@ad225aorg.highwire.dtl.DTLVardef@aa430eorg.highwire.dtl.DTLVardef@1880ca_HPS_FORMAT_FIGEXP M_FIG C_FIG HighlightsO_LINew compositional data approach improves metabolic flux estimation. C_LIO_LIThis data transformation requires minimal changes to existing workflows. C_LIO_LIThe new method reduced MSE of flux estimates by 42.6% in two examples tested. C_LIO_LIThe confidence intervals of the estimated fluxes were substantially narrowed. C_LIO_LIEstimation accuracy remained robust across a wide range of metabolic fluxes. C_LI

Electrical stimulation (ES) and transepithelial/transendothelial electrical resistance (TEER) measurements are essential techniques in cell biology and tissue engineering, yet commercial devices for these applications cost between USD 2,500-9,000 and typically offer only one functionality. We present LATEER (Low-cost Arduino-based TEER and Electrical stimulation device), an open-source hardware platform that combines both ES and TEER measurement capabilities at a total cost below USD 100. The device features four independent channels, configurable pulsatile signals (amplitude up to 8.2 V, frequency 0.1-500 Hz, pulse width [&ge;]0.1 ms), and a resistance measurement range of 300 {Omega} to 1 M{Omega}, with <5% error for R {gtrsim} 4.7 k{Omega}. LATEER uses commercially available graphite pencil leads as electrodes ([~]USD 2 vs. USD 350 for commercial Ag/AgCl electrodes), which demonstrated excellent biocompatibility in cell culture. The system includes 3D-printed electrode holders compatible with standard 12-well and 24-well plates, allowing microscope visualization without electrode removal, and a Python-based graphical user interface for parameter configuration and real-time data acquisition. Because the electrodes remain fixed in the plate lid and only a single cable enters the incubator, both stimulation and resistance measurement can run continuously under standard culture conditions (37 {degrees}C, 5% CO2) without removing the plate or repositioning the electrodes, avoiding the temperature excursions and placement variability inherent to manual chopstick measurements. Validation with human pluripotent stem cell-derived cardiomyocytes demonstrated reliable frequency capture (electrical pacing) of the contracting monolayer, with a capture threshold between 250 and 400 mV/mm and controlled pacing across the 0.5-5 Hz range. TEER functionality was verified with mesenchymal stem cells, where the device resolved cell-density-dependent differences in electrical resistance in real time. All design files, firmware, and software are freely available under the CERN-OHL-S v2 license, enabling replication and customization by research laboratories worldwide. HighlightsO_LIAn open-source device combines electrical stimulation and TEER measurement under $100 C_LIO_LIGraphite electrodes offer biocompatibility at 0.6% cost of commercial alternatives C_LIO_LIFour independent channels with configurable parameters and real-time data logging. C_LIO_LIContinuous run setup in-incubator; no electrode repositioning needed C_LIO_LIValidated with stem cell-derived cardiomyocytes, achieving frequency capture (threshold 250-400 mV/mm) C_LIO_LI3D-printed holders enable microscope visualization without electrode removal C_LI Graphical abstract O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=78 SRC="FIGDIR/small/743263v1_ufig1.gif" ALT="Figure 1"> View larger version (28K): org.highwire.dtl.DTLVardef@98f9deorg.highwire.dtl.DTLVardef@13c73aborg.highwire.dtl.DTLVardef@1cdf099org.highwire.dtl.DTLVardef@16ed0cf_HPS_FORMAT_FIGEXP M_FIG C_FIG Specifications Table O_TBL View this table: org.highwire.dtl.DTLVardef@4ef802org.highwire.dtl.DTLVardef@7c651borg.highwire.dtl.DTLVardef@d20013org.highwire.dtl.DTLVardef@102fc89org.highwire.dtl.DTLVardef@111b927_HPS_FORMAT_FIGEXP M_TBL C_TBL

Neural temporal point processes (NTPPs) are powerful tools for modeling sequences of timestamped events with statistical temporal structure. Density-based NTPPs, in particular, are an interesting opportunity to merge the universal function approximation capability of neural networks with a defined statistical model in a way that has many potential applications. We demonstrate one such application to heartbeat dynamics, a physiologic point process. We specifically apply a lognormal mixture NTPP to compute instantaneous estimates of the mean and standard deviation of beat-to-beat intervals. We compare our results to the state of art (Barbieri et al.) point process model for heartbeat dynamics, which uses a more physiologically rigorous inverse Gaussian model. We find that the NTPP model maintains reasonable accuracy while improving upon robustness to noise.

Pathological tremor is a neurological condition that impairs fine motor tasks, affecting 1% of the general population and 4% of the elderly. Tremors arise when muscles micro-oscillations synchronize and phase lock, typically within a 4-12 Hz frequency range. Administering beta-blockers can reduce tremor severity, but doses are hard to personalize, with heavy doses of propranolol correlating with low blood pressure, dizziness, and nausea. In this project, we aimed to model tremor and create a closed-loop control framework to suppress tremor amplitude while minimizing pharmacological dependence. Because side effects constrain the use of pharmacological suppression alone, we investigated noninvasive neuromodulation. We used vibrotactile stimulation (VTS) to disrupt pathological tremor synchronization and reduce oscillatory amplitude. We hypothesized that tremor suppression involving VTS followed a nonmonotonic relationship, tested by determining whether maximum relief requires an adaptable framework. The procedure consisted of constructing a propranolol-reduction simulation by implementing a Hill curve, where we calculated and utilized tremor reduction, heart rate (HR) drop, and blood pressure (BP) drop. We then built a device to capture tremor-related data and create vibration using two linear resonant actuator (LRA) coin motors. We connected it to a microcontroller, where we determined optimal vibration frequencies through a feedback loop. Across 50 trials, VTS alone reduced tremor amplitude by an average of 37.3%, reducing the propranolol dose needed to reach 50% total tremor reduction by 71.9%, lowering the modeled blood pressure drop from 38.1 to 18.9 mmHg. This device demonstrates proof-of-concept for a nonmonotonic tremor-vibration relationship to reduce dependency on propranolol in the treatment of pathological tremor. These propranolol dose-reduction estimates are derived from computational simulation and have not been clinically validated; they are not intended as a recommendation to alter prescribed medication.

Accurate tree volume estimation is central to forest management and carbon accounting. Allometric equations are widely used but limited in transferability across species, regions, and environmental conditions. Mobile Laser Scanning (MLS) offers a promising alternative through direct measurement of tree geometry; however, the influence of tree shape on MLS accuracy remains poorly understood. This study evaluated MLS-derived estimates of stem diameters, total tree height, and merchantable stem volume against destructive reference measurements from 176 trees spanning eight species (four hardwood, four softwood) in Wallonia, Belgium. A Zeb Horizon RT scanner was used; tree architectural descriptors extracted from the point cloud were tested for associations with measurement error. Across 7,824 stem diameter measurements, MLS achieved a mean error of 0.46 cm, with precision declining above 15 m. MLS-derived height outperformed Vertex IV clinometer measurements for hardwood species (RMSE% = 6.88 vs. 8.78) but performed slightly less well for softwoods (RMSE% = 7.36 vs. 6.14). QSM-based volume estimates systematically underestimated reference values, while taper-based reconstruction produced nearly unbiased estimates with an RMSE of 15.72%. Correlation analyses and PCA showed that tree architectural variables explained only a small fraction of MLS error variability. Diameter and height errors were largely independent of structural attributes, while volume errors showed moderate associations with tree size and crown density. These findings indicate that tree architecture is not a primary source of MLS measurement uncertainty. Future MLS-based forest inventory efforts should prioritize acquisition and processing optimization, as scanning conditions and forest structure appear more influential than tree shape.

Modern medicine implicitly assumes that physiological responses to intervention are predictably determined by administered treatments. However, physiological systems containing intrinsic delays between the detection of a stimulus and the biological response may violate this assumption. We investigate the human glucose-insulin system as described by the Ultradian model and mathematically demonstrate that clinically relevant forcing protocols-such as pulsatile insulin delivery and step-wise glucose infusion, both commonly used in intensive care units (ICUs)-can induce sustained temporal chaos that may hamper accurate prediction of the physiological response. If not accounted for, these chaotic dynamics could create difficulties in achieving optimal dosing and timing when administering glucose and insulin in clinical or home care settings. This phenomenon, termed delay-induced uncertainty (DIU), arises from the interaction between physiological delay, intrinsic shear near a limit cycle, and external forcing. Using the Ultradian glucose-insulin model, we compute top Lyapunov exponents to quantify predictability. Across a range of pulsatile and step-wise forcing regimes, including stochastic amplitudes drawn from Markov processes, we observe positive Lyapunov exponents, indicating sustained chaos. Our results suggest that delayed endocrine regulation may fundamentally limit the predictive value of the models used to develop glycemic management strategies, with implications for clinical protocols in the ICU.

Microorganisms interact through the exchange of metabolites and competition for shared substrates, and this metabolic coupling shapes the composition and function of microbial communities. Community flux balance analysis (cFBA) can predict such behaviour - the maximum community growth rate, the metabolic fluxes and the relative abundances of the species - from stoichiometric models of their metabolism, but existing formulations are either complex and hard to scale as communities grow or cannot predict optimal growth rates. Here we present a physiology-based formulation of cFBA in which each species' metabolism is reduced to a few macrochemical equations, one for each 'metabolic mode' the species can use, and the whole community is then solved as a single linear program. From this, the method predicts the optimal composition of the community, its maximum growth rate, the metabolites exchanged between the species, and the net conversion the community carries out as a whole; its ecological service. This reduction makes it far simpler to build and solve models of larger communities. We illustrate the approach on a two-species synergistic community that can be verified by hand, apply it to a five-member anaerobic digestion community, and use it to predict the metabolic interactions of a genome-scale syngas-fermenting coculture. Characterising these communities at their optimal steady states, we show that each species is driven to a distinct metabolic strategy. We discuss the method both as a practical tool for larger microbial communities and as a means of uncovering the ecological principles that govern them.

BackgroundContinuous glucose monitoring (CGM) systems are used not only for retrospective assessment of the glycemic profile but also for real-time decision-making, including automated insulin delivery. Accordingly, CGM performance characterization must capture not only the agreement of individual paired values but also the systems ability to reproduce the direction, rate, amplitude, and shape of glucose concentration change. Summary metrics, most notably MARD, cannot establish whether an observed deviation reflects an error in the formation of the test profile itself, a constant sensor offset, amplitude compression, a change in response rate, temporal misalignment, or hysteresis. ObjectiveTo adapt a programmable flow-based in vitro platform for the separate assessment of the experimentally delivered glucose profile and the dynamic response of CGM systems, and to propose a set of metrics that decomposes dynamic error into its components. MethodsGLU profiles were generated by programmable mixing of solutions at a constant total flow rate of 2 mL/min. Actual GLU concentration was independently measured with a SUPER GL2 glucose analyzer. Four static levels, three repeats of a 5.5[-&gt;]12.0[-&gt;]5.5 mmol/L profile, three repeats of a 6.0[-&gt;]3.0[-&gt;]6.0 mmol/L hypoglycemic profile, three 5.0[-&gt;]15.0[-&gt;]5.0 mmol/L profiles at different rates, one complex 4[-&gt;]18[-&gt;]3[-&gt;]12[-&gt;]5.5 mmol/L profile, and two proof-of-concept sensor experiments at 100- and 200-min transitions were investigated. Dynamic response was characterized by bias, MAE, RMSE, MARD, amplitude transfer coefficient K_A, rate transfer coefficients K_up and K_down, normalized shape RMSE, residual shift, and hysteresis loop area. ResultsAt the static levels, measured GLU exceeded the programmed value by 0.234-0.780 mmol/L. In the repeated 5.5[-&gt;]12.0[-&gt;]5.5 profiles, the ratio of actual to programmed rate was 0.978-1.083 on the rising phase and 0.987-1.157 on the falling phase, while the amplitude transfer coefficient was 0.967-1.066. In the hypoglycemic profile, minimum GLU was 2.55- 2.96 mmol/L, and time below 3.0 mmol/L was 15.2-72.6 min. The measured rates of 0.0519, 0.1045, and 0.2027 mmol/L/min preserved the intended ratio of approximately 1:2:4. In the complex profile, the programmed plateau of 18 mmol/L was not reached: mean measured GLU was 16.20 mmol/L. For CGM-A, K_A was 0.682 and 0.650, and K_up/K_down were 0.666/0.730 and 0.634/0.626; the corresponding values for CGM-B were 1.228 and 1.128, and 1.564/1.328 and 1.276/1.145. Hysteresis loop area differed 5- to 10-fold between the two sensor responses, exceeding an order of magnitude at the 100-min transition. ConclusionThe programmed concentration should be treated as a control setpoint, rather than as a reference measurement. The "programmed trajectory -- measured glucose -- CGM output" cascade first allows quantitative assessment of the agreement between the programmed and actually realized profile and only then separate characterization of sensor response. Decomposition of dynamic error into amplitude, rate, shape, and hysteresis components reveals differences that a single MARD value or correlation coefficient cannot capture.

Continuous wave diffuse correlation spectroscopy (cw-DCS) is a noninvasive optical technique to monitor the tissues blood flow changes. This technique measures the tissue blood flow index (BFI) by evaluating the decay rate of the autocorrelation function. The derived BFI is proportional to mean squared displacements of the red blood cells considered as the fast-dynamic scatterer component of tissue in time. However, biological tissue contains static scatterer component and slow-dynamic scatterer component which affect the decay rate of autocorrelation function and as a result the derived BFI. In this study, we assessed the fractional contribution of static, slow-dynamic and fast-dynamic scatterer components of a medium in the flow index derived by cw-DCS. The measurements performed on Agar-based phantom with tube showed that presence of static scatterer component and slow-dynamic scatterer component led to substantial underestimation ({approx} 123%) of the flow index derived by Siegert relation, compared to effective diffusion coefficient of fast-dynamic scatterers components derived by modified Siegert relation and bi-exponential model. The less underestimation was observed for the corresponding parameters obtained from the liquid phantom measurements ({approx} 25%) as well as during the forearm occlusion test and respiratory challenges ({approx} 16% - 26%).

Artificial intelligence (AI) is increasingly permeating the drug development pipeline. Numerous algorithms for accelerating this multi-stage and multi-task process have been constructed, which depends heavily on expert design and labor-intensive task-specific optimization. Given that AI-driven acceleration of drug development is recognized as a cumulative, often synergistic, effect across multiple stages, the autonomous evolution of existing algorithms across the entire pipeline is demanded to achieve a holistic advancement. Here, we present DrugEvolve, a multi-role large language model system for systematic and autonomous algorithm evolution in drug development. DrugEvolve realizes a closed-loop evolution process by incorporating Researcher, Engineer, and Analyst domains, and enables an iterative design, implementation, evaluation, and refinement of algorithm by leveraging scientific knowledge and accumulated evolutionary experience. Across eleven representative tasks spanning target identification, drug discovery, preclinical study, and clinical trial, DrugEvolve autonomously evolved the corresponding task-specific algorithms and achieved substantial performance enhancement on 120 benchmark test sets. Moreover, it showed robust generalizabilities across heterogeneous data modalities (ranging from biological sequence and graph to molecular topology and textual language), and realized gains in both predictive and generative tasks. Collectively, this AI system can serve not only as an algorithmic infrastructure for drug development, but also as a transferable paradigm for broader scientific domains.