Design patterns for the construction of computational biological models.

Computational biological models have proven to be an invaluable tool for understanding and predicting the behaviour of many biological systems. While it may not be too challenging for experienced researchers to construct such models from scratch, it is not a straightforward task for early stage researchers. Design patterns are well-known techniques widely applied in software engineering as they provide a set of typical solutions to common problems in software design. In this paper, we collect and discuss common patterns that are usually used during the construction and execution of computational biological models. We adopt Petri nets as a modelling language to provide a visual illustration of each pattern; however, the ideas presented in this paper can also be implemented using other modelling formalisms. We provide two case studies for illustration purposes and show how these models can be built up from the presented smaller modules. We hope that the ideas discussed in this paper will help many researchers in building their own future models.

The ReCoDe addiction research consortium: Losing and regaining control over drug intake-Findings and future perspectives.

Substance use disorders (SUDs) are seen as a continuum ranging from goal-directed and hedonic drug use to loss of control over drug intake with aversive consequences for mental and physical health and social functioning. The main goals of our interdisciplinary German collaborative research centre on Losing and Regaining Control over Drug Intake (ReCoDe) are (i) to study triggers (drug cues, stressors, drug priming) and modifying factors (age, gender, physical activity, cognitive functions, childhood adversity, social factors, such as loneliness and social contact/interaction) that longitudinally modulate the trajectories of losing and regaining control over drug consumption under real-life conditions. (ii) To study underlying behavioural, cognitive and neurobiological mechanisms of disease trajectories and drug-related behaviours and (iii) to provide non-invasive mechanism-based interventions. These goals are achieved by: (A) using innovative mHealth (mobile health) tools to longitudinally monitor the effects of triggers and modifying factors on drug consumption patterns in real life in a cohort of 900 patients with alcohol use disorder. This approach will be complemented by animal models of addiction with 24/7 automated behavioural monitoring across an entire disease trajectory; i.e. from a naïve state to a drug-taking state to an addiction or resilience-like state. (B) The identification and, if applicable, computational modelling of key molecular, neurobiological and psychological mechanisms (e.g., reduced cognitive flexibility) mediating the effects of such triggers and modifying factors on disease trajectories. (C) Developing and testing non-invasive interventions (e.g., Just-In-Time-Adaptive-Interventions (JITAIs), various non-invasive brain stimulations (NIBS), individualized physical activity) that specifically target the underlying mechanisms for regaining control over drug intake. Here, we will report on the most important results of the first funding period and outline our future research strategy.

Adsorption in Pyrene-Based Metal-Organic Frameworks: The Role of Pore Structure and Topology.

Pore topology and chemistry play crucial roles in the adsorption characteristics of metal-organic frameworks (MOFs). To deepen our understanding of the interactions between MOFs and CO2 during this process, we systematically investigate the adsorption properties of a group of pyrene-based MOFs. These MOFs feature Zn(II) as the metal ion and employ a pyrene-based ligand, specifically 1,3,6,8-tetrakis(p-benzoic acid)pyrene (TBAPy). Including different additional ligands leads to frameworks with distinctive structural and chemical features. By comparing these structures, we could isolate the role that pore size, the presence of open-metal sites (OMS), metal-oxygen bridges, and framework charges play in the CO2 adsorption of these MOFs. Frameworks with constricted pore structures display a phenomenon known as the confinement effect, fostering stronger MOF-CO2 interactions and higher uptakes at low pressures. In contrast, entropic effects dominate at elevated pressures, and the MOF's pore volume becomes the driving factor. Through analysis of the CO2 uptakes of the benchmark materials ─some with narrower pores and others with larger pore volumes─it becomes evident that structures with narrower pores and high binding energies excel at low pressures. In contrast, those with larger volumes perform better at elevated pressures. Moreover, this research highlights that open-metal sites and inherent charges within the frameworks of ionic MOFs stand out as CO2-philic characteristics.

Machine Learning for Predicting Biologic Agent Efficacy in Ulcerative Colitis: An Analysis for Generalizability and Combination with Computational Models.

Machine learning (ML) has been applied to predict the efficacy of biologic agents in ulcerative colitis (UC). ML can offer precision, personalization, efficiency, and automation. Moreover, it can improve decision support in predicting clinical outcomes. However, it faces challenges related to data quality and quantity, overfitting, generalization, and interpretability. This paper comments on two recent ML models that predict the efficacy of vedolizumab and ustekinumab in UC. Models that consider multiple pathways, multiple ethnicities, and combinations of real-world and clinical trial data are required for optimal shared decision-making and precision medicine. This paper also highlights the potential of combining ML with computational models to enhance clinical outcomes and personalized healthcare. Key Insights: (1) ML offers precision, personalization, efficiency, and decision support for predicting the efficacy of biologic agents in UC. (2) Challenging aspects in ML prediction include data quality, overfitting, and interpretability. (3) Multiple pathways, multiple ethnicities, and combinations of real-world and clinical trial data should be considered in predictive models for optimal decision-making. (4) Combining ML with computational models may improve clinical outcomes and personalized healthcare.

Corrigendum: Biomarkers and computational models for predicting efficacy to tumor ICI immunotherapy.

[This corrects the article DOI: 10.3389/fimmu.2024.1368749.].

Neurocomputational Models of Interval Timing: Seeing the Forest for the Trees.

Extracting temporal regularities and relations from experience/observation is critical for organisms' adaptiveness (communication, foraging, predation, prediction) in their ecological niches. Therefore, it is not surprising that the internal clock that enables the perception of seconds-to-minutes-long intervals (interval timing) is evolutionarily well-preserved across many species of animals. This comparative claim is primarily supported by the fact that the timing behavior of many vertebrates exhibits common statistical signatures (e.g., on-average accuracy, scalar variability, positive skew). These ubiquitous statistical features of timing behaviors serve as empirical benchmarks for modelers in their efforts to unravel the processing dynamics of the internal clock (namely answering how internal clock "ticks"). In this chapter, we introduce prominent (neuro)computational approaches to modeling interval timing at a level that can be understood by general audience. These models include Treisman's pacemaker accumulator model, the information processing variant of scalar expectancy theory, the striatal beat frequency model, behavioral expectancy theory, the learning to time model, the time-adaptive opponent Poisson drift-diffusion model, time cell models, and neural trajectory models. Crucially, we discuss these models within an overarching conceptual framework that categorizes different models as threshold vs. clock-adaptive models and as dedicated clock/ramping vs. emergent time/population code models.

A Review of Personalised Cardiac Computational Modelling Using Electroanatomical Mapping Data.

Computational models of cardiac electrophysiology have gradually matured during the past few decades and are now being personalised to provide patient-specific therapy guidance for improving suboptimal treatment outcomes. The predictive features of these personalised electrophysiology models hold the promise of providing optimal treatment planning, which is currently limited in the clinic owing to reliance on a population-based or average patient approach. The generation of a personalised electrophysiology model entails a sequence of steps for which a range of activation mapping, calibration methods and therapy simulation pipelines have been suggested. However, the optimal methods that can potentially constitute a clinically relevant in silico treatment are still being investigated and face limitations, such as uncertainty of electroanatomical data recordings, generation and calibration of models within clinical timelines and requirements to validate or benchmark the recovered tissue parameters. This paper is aimed at reporting techniques on the personalisation of cardiac computational models, with a focus on calibrating cardiac tissue conductivity based on electroanatomical mapping data.

D2/D3 dopamine supports the precision of mental state inferences and self-relevance of joint social outcomes.

Striatal dopamine is important in paranoid attributions, although its computational role in social inference remains elusive. We employed a simple game-theoretic paradigm and computational model of intentional attributions to investigate the effects of dopamine D2/D3 antagonism on ongoing mental state inference following social outcomes. Haloperidol, compared with the placebo, enhanced the impact of partner behaviour on beliefs about the harmful intent of partners, and increased learning from recent encounters. These alterations caused substantial changes to model covariation and negative correlations between self-interest and harmful intent attributions. Our findings suggest that haloperidol improves belief flexibility about others and simultaneously reduces the self-relevance of social observations. Our results may reflect the role of D2/D3 dopamine in supporting self-relevant mentalising. Our data and model bridge theory between general and social accounts of value representation. We demonstrate initial evidence for the sensitivity of our model and short social paradigm to drug intervention and clinical dimensions, allowing distinctions between mechanisms that operate across traits and states.

Modulation of neuronal dynamics by sustained and activity-dependent continuous-wave near-infrared laser stimulation.

Significance: Near-infrared laser illumination is a non-invasive alternative/complement to classical stimulation methods in neuroscience but the mechanisms underlying its action on neuronal dynamics remain unclear. Most studies deal with high-frequency pulsed protocols and stationary characterizations disregarding the dynamic modulatory effect of sustained and activity-dependent stimulation. The understanding of such modulation and its widespread dissemination can help to develop specific interventions for research applications and treatments for neural disorders. Aim: We quantified the effect of continuous-wave near-infrared (CW-NIR) laser illumination on single neuron dynamics using sustained stimulation and an open-source activity-dependent protocol to identify the biophysical mechanisms underlying this modulation and its time course. Approach: We characterized the effect by simultaneously performing long intracellular recordings of membrane potential while delivering sustained and closed-loop CW-NIR laser stimulation. We used waveform metrics and conductance-based models to assess the role of specific biophysical candidates on the modulation. Results: We show that CW-NIR sustained illumination asymmetrically accelerates action potential dynamics and the spiking rate on single neurons, while closed-loop stimulation unveils its action at different phases of the neuron dynamics. Our model study points out the action of CW-NIR on specific ionic-channels and the key role of temperature on channel properties to explain the modulatory effect. Conclusions: Both sustained and activity-dependent CW-NIR stimulation effectively modulate neuronal dynamics by a combination of biophysical mechanisms. Our open-source protocols can help to disseminate this non-invasive optical stimulation in novel research and clinical applications.

From silos to synergy: Integrating approaches to investigate the role of prior knowledge and expectations on episodic memory.

Significant progress in the investigation of how prior knowledge influences episodic memory has been made using three sometimes isolated (but not mutually exclusive) approaches: strictly adult behavioral investigations, computational models, and investigations into the development of the system. Here we point out that these approaches are complementary, each approach informs and is informed by the other. Thus, a natural next step for research is to combine all three approaches to further our understanding of the role of prior knowledge in episodic memory. Here we use studies of memory for expectation-congruent and incongruent information from each of these often disparate approaches to illustrate how combining approaches can be used to test and revise theories from the other. This domain is particularly advantageous because it highlights important features of more general memory processes, further differentiates models of memory, and can shed light on developmental change in the memory system. We then present a case study to illustrate the progress that can be made from integrating all three approaches and highlight the need for more endeavors in this vein. As a first step, we also propose a new computational model of memory that takes into account behavioral and developmental factors that can influence prior knowledge and episodic memory interactions. This integrated approach has great potential for offering novel insights into the relationship between prior knowledge and episodic memory, and cognition more broadly.

Modelling the impact of wastewater flows and management practices on antimicrobial resistance in dairy farms.

Dairy slurry is a major source of environmental contamination with antimicrobial resistant genes and bacteria. We developed mathematical models and conducted on-farm research to explore the impact of wastewater flows and management practices on antimicrobial resistance (AMR) in slurry. Temporal fluctuations in cephalosporin-resistant Escherichia coli were observed and attributed to farm activities, specifically the disposal of spent copper and zinc footbath into the slurry system. Our model revealed that resistance should be more frequently observed with relevant determinants encoded chromosomally rather than on plasmids, which was supported by reanalysis of sequenced genomes from the farm. Additionally, lower resistance levels were predicted in conditions with lower growth and higher death rates. The use of muck heap effluent for washing dirty channels did not explain the fluctuations in cephalosporin resistance. These results highlight farm-specific opportunities to reduce AMR pollution, beyond antibiotic use reduction, including careful disposal or recycling of waste antimicrobial metals.

The neural representation of an auditory spatial cue in the primate cortex.

Humans make use of small differences in the timing of sounds at the two ears-interaural time differences (ITDs)-to locate their sources. Despite extensive investigation, however, the neural representation of ITDs in the human brain is contentious, particularly the range of ITDs explicitly represented by dedicated neural detectors. Here, using magneto- and electro-encephalography (MEG and EEG), we demonstrate evidence of a sparse neural representation of ITDs in the human cortex. The magnitude of cortical activity to sounds presented via insert earphones oscillated as a function of increasing ITD-within and beyond auditory cortical regions-and listeners rated the perceptual quality of these sounds according to the same oscillating pattern. This pattern was accurately described by a population of model neurons with preferred ITDs constrained to the narrow, sound-frequency-dependent range evident in other mammalian species. When scaled for head size, the distribution of ITD detectors in the human cortex is remarkably like that recorded in vivo from the cortex of rhesus monkeys, another large primate that uses ITDs for source localization. The data solve a long-standing issue concerning the neural representation of ITDs in humans and suggest a representation that scales for head size and sound frequency in an optimal manner.

A mechanobiological model of bone metastasis reveals that mechanical stimulation inhibits the pro-osteolytic effects of breast cancer cells.

Bone is highly susceptible to cancer metastasis, and both tumor and bone cells enable tumor invasion through a "vicious cycle" of biochemical signaling. Tumor metastasis into bone also alters biophysical cues to both tumor and bone cells, which are highly sensitive to their mechanical environment. However, the mechanobiological feedback between these cells that perpetuate this cycle has not been studied. Here, we develop highly advanced in vitro and computational models to provide an advanced understanding of how tumor growth is regulated by the synergistic influence of tumor-bone cell signaling and mechanobiological cues. In particular, we develop a multicellular healthy and metastatic bone model that can account for physiological mechanical signals within a custom bioreactor. These models successfully recapitulated mineralization, mechanobiological responses, osteolysis, and metastatic activity. Ultimately, we demonstrate that mechanical stimulus provided protective effects against tumor-induced osteolysis, confirming the importance of mechanobiological factors in bone metastasis development.

To what extent can mastication functionality be restored following mandibular reconstruction surgery? A computer modeling approach.

STATEMENT OF PROBLEM: Advanced cases of head and neck cancer involving the mandible often require surgical removal of diseased sections and subsequent replacement with donor bone. During the procedure, the surgeon must make decisions regarding which bones or tissues to resect. This requires balancing tradeoffs related to issues such as surgical access and post-operative function; however, the latter is often difficult to predict, especially given that long-term functionality also depends on the impact of post-operative rehabilitation programs. PURPOSE: To assist in surgical decision-making, we present an approach for estimating the effects of reconstruction on key aspects of post-operative mandible function. MATERIAL AND METHODS: We develop dynamic biomechanical models of the reconstructed mandible considering different defect types and validate them using literature data. We use these models to estimate the degree of functionality that might be achieved following post-operative rehabilitation. RESULTS: We find significant potential for restoring mandibular functionality, even in cases involving large defects. This entails an average trajectory error below 2 mm, bite force comparable to a healthy individual, improved condyle mobility, and a muscle activation change capped at a maximum of 20%. CONCLUSION: These results suggest significant potential for adaptability in the masticatory system and improved post-operative rehabilitation, leading to greater restoration of jaw function.

Lexical Feedback in the Time-Invariant String Kernel (TISK) Model of Spoken Word Recognition.

The Time-Invariant String Kernel (TISK) model of spoken word recognition (Hannagan, Magnuson & Grainger, 2013; You & Magnuson, 2018) is an interactive activation model with many similarities to TRACE (McClelland & Elman, 1986). However, by replacing most time-specific nodes in TRACE with time-invariant open-diphone nodes, TISK uses orders of magnitude fewer nodes and connections than TRACE. Although TISK performed remarkably similarly to TRACE in simulations reported by Hannagan et al., the original TISK implementation did not include lexical feedback, precluding simulation of top-down effects, and leaving open the possibility that adding feedback to TISK might fundamentally alter its performance. Here, we demonstrate that when lexical feedback is added to TISK, it gains the ability to simulate top-down effects without losing the ability to simulate the fundamental phenomena tested by Hannagan et al. Furthermore, with feedback, TISK demonstrates graceful degradation when noise is added to input, although parameters can be found that also promote (less) graceful degradation without feedback. We review arguments for and against feedback in cognitive architectures, and conclude that feedback provides a computationally efficient basis for robust constraint-based processing.

Joint EURADOS-EANM initiative for an advanced computational framework for the assessment of external dose rates from nuclear medicine patients.

BACKGROUND: In order to ensure adequate radiation protection of critical groups such as staff, caregivers and the general public coming into proximity of nuclear medicine (NM) patients, it is necessary to consider the impact of the radiation emitted by the patients during their stay at the hospital or after leaving the hospital. Current risk assessments are based on ambient dose rate measurements in a single position at a specified distance from the patient and carried out at several time points after administration of the radiopharmaceutical to estimate the whole-body retention. The limitations of such an approach are addressed in this study by developing and validating a more advanced computational dosimetry approach using Monte Carlo (MC) simulations in combination with flexible and realistic computational phantoms and time activity distribution curves from reference biokinetic models. RESULTS: Measurements of the ambient dose rate equivalent H*(10) at 1 m from the NM patient have been successfully compared against MC simulations with 5 different codes using the ICRP adult reference computational voxel phantoms, for typical clinical procedures with 99mTc-HDP/MDP, 18FDG and Na131I. All measurement data fall in the 95% confidence intervals, determined for the average simulated results. Moreover, the different MC codes (MCNP-X, PHITS, GATE, GEANT4, TRIPOLI-4®) have been compared for a more realistic scenario where the effective dose rate E of an exposed individual was determined in positions facing and aside the patient model at 30 cm, 50 cm and 100 cm. The variation between codes was lower than 8% for all the radiopharmaceuticals at 1 m, and varied from 5 to 16% for the face-to face and side-by-side configuration at 30 cm and 50 cm. A sensitivity study on the influence of patient model morphology demonstrated that the relative standard deviation of H*(10) at 1 m for the range of included patient models remained under 16% for time points up to 120 min post administration. CONCLUSIONS: The validated computational approach will be further used for the evaluation of effective dose rates per unit administered activity for a variety of close-contact configurations and a range of radiopharmaceuticals as part of risk assessment studies. Together with the choice of appropriate dose constraints this would facilitate the setting of release criteria and patient restrictions.

Biomarkers and computational models for predicting efficacy to tumor ICI immunotherapy.

Numerous studies have shown that immune checkpoint inhibitor (ICI) immunotherapy has great potential as a cancer treatment, leading to significant clinical improvements in numerous cases. However, it benefits a minority of patients, underscoring the importance of discovering reliable biomarkers that can be used to screen for potential beneficiaries and ultimately reduce the risk of overtreatment. Our comprehensive review focuses on the latest advancements in predictive biomarkers for ICI therapy, particularly emphasizing those that enhance the efficacy of programmed cell death protein 1 (PD-1)/programmed cell death-ligand 1 (PD-L1) inhibitors and cytotoxic T-lymphocyte antigen-4 (CTLA-4) inhibitors immunotherapies. We explore biomarkers derived from various sources, including tumor cells, the tumor immune microenvironment (TIME), body fluids, gut microbes, and metabolites. Among them, tumor cells-derived biomarkers include tumor mutational burden (TMB) biomarker, tumor neoantigen burden (TNB) biomarker, microsatellite instability (MSI) biomarker, PD-L1 expression biomarker, mutated gene biomarkers in pathways, and epigenetic biomarkers. TIME-derived biomarkers include immune landscape of TIME biomarkers, inhibitory checkpoints biomarkers, and immune repertoire biomarkers. We also discuss various techniques used to detect and assess these biomarkers, detailing their respective datasets, strengths, weaknesses, and evaluative metrics. Furthermore, we present a comprehensive review of computer models for predicting the response to ICI therapy. The computer models include knowledge-based mechanistic models and data-based machine learning (ML) models. Among the knowledge-based mechanistic models are pharmacokinetic/pharmacodynamic (PK/PD) models, partial differential equation (PDE) models, signal networks-based models, quantitative systems pharmacology (QSP) models, and agent-based models (ABMs). ML models include linear regression models, logistic regression models, support vector machine (SVM)/random forest/extra trees/k-nearest neighbors (KNN) models, artificial neural network (ANN) and deep learning models. Additionally, there are hybrid models of systems biology and ML. We summarized the details of these models, outlining the datasets they utilize, their evaluation methods/metrics, and their respective strengths and limitations. By summarizing the major advances in the research on predictive biomarkers and computer models for the therapeutic effect and clinical utility of tumor ICI, we aim to assist researchers in choosing appropriate biomarkers or computer models for research exploration and help clinicians conduct precision medicine by selecting the best biomarkers.

The Causes for Genomic Instability and How to Try and Reduce Them Through Rational Design of Synthetic DNA.

Genetic engineering has revolutionized our ability to manipulate DNA and engineer organisms for various applications. However, this approach can lead to genomic instability, which can result in unwanted effects such as toxicity, mutagenesis, and reduced productivity. To overcome these challenges, smart design of synthetic DNA has emerged as a promising solution. By taking into consideration the intricate relationships between gene expression and cellular metabolism, researchers can design synthetic constructs that minimize metabolic stress on the host cell, reduce mutagenesis, and increase protein yield. In this chapter, we summarize the main challenges of genomic instability in genetic engineering and address the dangers of unknowingly incorporating genomically unstable sequences in synthetic DNA. We also demonstrate the instability of those sequences by the fact that they are selected against conserved sequences in nature. We highlight the benefits of using ESO, a tool for the rational design of DNA for avoiding genetically unstable sequences, and also summarize the main principles and working parameters of the software that allow maximizing its benefits and impact.

Bet-hedging and variability in plant development: seed germination and beyond.

When future conditions are unpredictable, bet-hedging strategies can be advantageous. This can involve isogenic individuals producing different phenotypes, under the same environmental conditions. Ecological studies provide evidence that variability in seed germination time has been selected for as a bet-hedging strategy. We demonstrate how variability in germination time found in Arabidopsis could function as a bet-hedging strategy in the face of unpredictable lethal stresses. Despite a body of knowledge on how the degree of seed dormancy versus germination is controlled, relatively little is known about how differences between isogenic seeds in a batch are generated. We review proposed mechanisms for generating variability in germination time and the current limitations and new possibilities for testing the model predictions. We then look beyond germination to the role of variability in seedling and adult plant growth and review new technologies for quantification of noisy gene expression dynamics. We discuss evidence for phenotypic variability in plant traits beyond germination being under genetic control and propose that variability in stress response gene expression could function as a bet-hedging strategy. We discuss open questions about how noisy gene expression could lead to between-plant heterogeneity in gene expression and phenotypes. This article is part of a discussion meeting issue 'Causes and consequences of stochastic processes in development and disease'.

Personalized modeling of Alzheimer's disease progression estimates neurodegeneration severity from EEG recordings.

INTRODUCTION: Early identification of Alzheimer's disease (AD) is necessary for a timely onset of therapeutic care. However, cortical structural alterations associated with AD are difficult to discern. METHODS: We developed a cortical model of AD-related neurodegeneration accounting for slowing of local dynamics and global connectivity degradation. In a monocentric study we collected electroencephalography (EEG) recordings at rest from participants in healthy (HC, n = 17), subjective cognitive decline (SCD, n = 58), and mild cognitive impairment (MCI, n = 44) conditions. For each patient, we estimated neurodegeneration model parameters based on individual EEG recordings. RESULTS: Our model outperformed standard EEG analysis not only in discriminating between HC and MCI conditions (F1 score 0.95 vs 0.75) but also in identifying SCD patients with biological hallmarks of AD in the cerebrospinal fluid (recall 0.87 vs 0.50). DISCUSSION: Personalized models could (1) support classification of MCI, (2) assess the presence of AD pathology, and (3) estimate the risk of cognitive decline progression, based only on economical and non-invasive EEG recordings. Highlights: Personalized cortical model estimating structural alterations from EEG recordings.Discrimination of Mild Cognitive Impairment (MCI) and Healthy (HC) subjects (95%)Prediction of biological markers of Alzheimer's in Subjective Decline (SCD) Subjects (87%)Transition correctly predicted for 3/3 subjects that converted from SCD to MCI after 1y.