Journal of The Royal Society Interface最新文献

Selective scleral crosslinking has been proposed as a novel treatment to increase scleral stiffness to counteract biomechanical changes associated with glaucoma and high myopia. Scleral stiffening has been shown by transpupillary peripapillary scleral photocrosslinking in rats, where the photosensitizer, methylene blue (MB), was injected retrobulbarly and red light initiated crosslinking reactions with collagen. Here, we adapted a computational model previously developed to model this treatment in rat eyes to additionally model MB photocrosslinking in minipigs and humans. Increased tissue length and subsequent diffusion and light penetration limitations were found to be barriers to achieving the same extent of crosslinking as in rats. Per cent inspired O₂, injected MB concentration and laser fluence were simultaneously varied to overcome these limitations and used to determine optimal combinations of treatment parameters in rats, minipigs and humans. Increasing these three treatment parameters simultaneously resulted in maximum crosslinking, except in rats, where the highest MB concentrations decreased crosslinking. Additionally, the kinetics and diffusion of photocrosslinking reaction intermediates and unproductive side products were modelled across space and time. The model provides a mechanistic understanding of MB photocrosslinking in scleral tissue and a basis for adapting and screening treatment parameters in larger animal models and, eventually, human eyes.

The speed of evolution on structured populations is crucial for biological and social systems. The likelihood of invasion is key for evolutionary stability. But it makes little sense if it takes long. It is far from known what population structure slows down evolution. We investigate the absorption time of a single neutral mutant for all the 112 non-isomorphic undirected graphs of size 6. We find that about three-quarters of the graphs have an absorption time close to that of the complete graph, less than one-third are accelerators, and more than two-thirds are decelerators. Surprisingly, determining whether a graph has a long absorption time is too complicated to be captured by the joint degree distribution. Via the largest sojourn time, we find that echo-chamber-like graphs, which consist of two homogeneous graphs connected by few sparse links, are likely to slow down absorption. These results are robust for large graphs, mutation patterns as well as evolutionary processes. This work serves as a benchmark for timing evolution with complex interactions, and fosters the understanding of polarization in opinion formation.

The health and well-being of a host are deeply influenced by the interactions with its gut microbiota. Contrasted environmental conditions, such as diseases or dietary habits, play a pivotal role in modulating these interactions, impacting microbiota composition and functionality. Such conditions can also lead to transitions from beneficial to detrimental symbiosis, viewed as alternative stable states of the host-microbiota dialogue. This article introduces a novel mathematical model exploring host-microbiota interactions, integrating dynamics of the colonic epithelial crypt, microbial metabolic functions, inflammation sensitivity and colon flows in a transverse section. The model considers metabolic shifts in epithelial cells based on butyrate and hydrogen sulfide concentrations, innate immune pattern recognition receptor activation, microbial oxygen tolerance and the impact of antimicrobial peptides on the microbiota. Using the model, we demonstrated that a high-protein, low-fibre diet exacerbates detrimental interactions and compromises beneficial symbiotic resilience, underscoring a destabilizing effect towards an unhealthy state. Moreover, the proposed model provides essential insights into oxygen levels, fibre and protein breakdown, and basic mechanisms of innate immunity in the colon and offers a crucial understanding of factors influencing the colon environment.

Many animals employ a second frequency filter beyond the initial filtering of the eardrum (or tympanal membrane). In the field cricket ear, both the filtering mechanism and the transmission path from the posterior tympanal membrane (PTM) have remained unclear. A mismatch between PTM vibrations and sensilla tuning has prompted speculations of a second filter. PTM coupling to the tracheal branches is suggested to support a transmission pathway. Here, we present three independent lines of evidence converging on the same conclusion: the existence of a series of linked membranes with distinct resonant frequencies serving both filtering and transmission functions. Micro-computed tomography (µ-CT) highlighted the 'dividing membrane (DivM)', separating the tracheal branches and connected to the PTM via the dorsal membrane of the posterior tracheal branch (DM-PTB). Thickness analysis showed the DivM to share significant thinness similarity with the PTM. Laser Doppler vibrometry indicated the first of two PTM vibrational peaks, at 6 and 14 kHz, originates not from the PTM but from the coupled DM-PTB. This result was corroborated by µ-CT-based finite element analysis. These findings clarify further the biophysical source of neuroethological pathways in what is an important model of behavioural neuroscience. Tuned microscale coupled membranes may also hold biomimetic relevance.

The component Allee effect (AE) is the positive correlation between an organism's fitness component and population density. Depending on the population spatial structure, which determines the interactions between organisms, a component AE might lead to positive density dependence in the population per-capita growth rate and establish a demographic AE. However, existing spatial models impose a fixed population spatial structure, which limits the understanding of how a component AE and spatial dynamics jointly determine the existence of demographic AEs. We introduce a spatially explicit theoretical framework where spatial structure and population dynamics are emergent properties of the individual-level demographic and movement rates. This framework predicts various spatial patterns depending on its specific parametrization, including evenly spaced aggregates of organisms, which determine the demographic-level by-products of the component AE. We find that aggregation increases population abundance and allows population survival in harsher environments and at lower global population densities when compared with uniformly distributed organisms. Moreover, aggregation can prevent the component AE from manifesting at the population level or restrict it to the level of each independent aggregate. These results provide a mechanistic understanding of how component AEs might operate for different spatial structures and manifest at larger scales.

Cell polarity is important for controlling cell shape, motility and cell division processes. Vimentin intermediate filaments are important for cell migration and cell polarization in mesenchymal cells and assembly of vimentin and microtubule networks is dynamically coordinated, but the precise details of how vimentin mediates cell polarity remain unclear. Here, we characterize the effects of vimentin on the structure and function of the centrosome and the stability of microtubule filaments in wild-type and vimentin-null mouse embryonic fibroblasts. We find that vimentin mediates the structure of the pericentriolar material, promotes centrosome-mediated microtubule regrowth and increases the level of stable acetylated microtubules in the cell. Loss of vimentin also impairs centrosome repositioning during cell polarization and migration processes that occur during wound closure. Our results suggest that vimentin modulates centrosome structure and function as well as microtubule network stability, which has important implications for how cells establish proper cell polarization and persistent migration.

The Earth's magnetic field can provide reliable directional information, allowing migrating animals to orient themselves using a magnetic compass or estimate their position relative to a target using map-based orientation. Here we show for the first time that young, inexperienced herring (Clupea harengus, Ch) have a magnetic compass when they migrate hundreds of kilometres to their feeding grounds. In birds, such as the European robin (Erithacus rubecula), radical pair-based magnetoreception involving cryptochrome 4 (ErCRY4) was demonstrated; the molecular basis of magnetoreception in fish is still elusive. We show that cry4 expression in the eye of herring is upregulated during the migratory season, but not before, indicating a possible use for migration. The amino acid structure of herring ChCRY4 shows four tryptophans and a flavin adenine dinucleotide-binding site, a prerequisite for a magnetic receptor. Using homology modelling, we successfully reconstructed ChCRY4 of herring, DrCRY4 of zebrafish (Danio rerio) and StCRY4 of brown trout (Salmo trutta) and showed that ChCRY4, DrCRY4 and ErCRY4a, but not StCRY4, exhibit very comparable dynamic behaviour. The electron transfer could take place in ChCRY4 in a similar way to ErCRY4a. The combined behavioural, transcriptomic and simulation experiments provide evidence that CRY4 could act as a magnetoreceptor in Atlantic herring.

Here, employing computer simulation tools, we present a study on the development of a bacterial biofilm from a single starter cell on a flat inert surface overlaid by an aqueous solution containing nutrients. In our simulations, surface colonization involves an initial stage of two-dimensional cell proliferation to eventually transition to three-dimensional growth leading to the formation of biofilm colonies with characteristic three-dimensional semi-ellipsoids shapes. Thus, we have introduced the influence of the nutrient concentration on bacterial growth, and calculated the cell growth rate as a function of nutrient uptake, which in turn depends on local nutrient concentration in the vicinity of each bacterial cell. Our results show that the combination of cell growth and nutrient uptake and diffusion leads to the formation of stratified colonies containing an inner core in which nutrients are depleted and cells cannot grow or divide, surrounded by an outer, shallow crust in which cells have access to nutrients from the bulk medium and continue growing. This phenomenon is more apparent at high uptake rates that enable fast nutrient depletion. Our simulations also predict that the shape and internal structure of the biofilm are largely conditioned by the balance between nutrient diffusion and uptake.

Virtual balancing tasks facilitate the study of human motion control: human reaction to the change of artificially introduced parameters can be studied in a computer environment. In this article, the dynamics of human stick balancing are generalized using fractional-order derivatives. Reaction delay sets a strong limitation on the length of the shortest stick that human subjects can balance. Human processing of visual input also exhibits a memory effect, which can be modelled by fractional-order derivatives. Therefore, we hypothesize a delayed fractional-order PD control of the unstable fractional-order process. The resulting equation of motion is investigated in a dimensionless framework, and stabilizability limits are determined as a function of the dynamics's order. These theoretical limits are then compared with the results of a systematic series of virtual balancing tests performed by 18 subjects. The comparison shows that the theoretical stabilizability limits for controllers with fixed fractional order correspond to the measured data points. The best fit is obtained if the fractional order of the underlying control law is 0.475.

In recent years, blending mechanistic knowledge with machine learning has had a major impact in digital healthcare. In this work, we introduce a computational pipeline to build certified digital replicas of cardiac electrophysiology in paediatric patients with congenital heart disease. We construct the patient-specific geometry by means of semi-automatic segmentation and meshing tools. We generate a dataset of electrophysiology simulations covering cell-to-organ level model parameters and using rigorous mathematical models based on differential equations. We previously proposed Branched Latent Neural Maps (BLNMs) as an accurate and efficient means to recapitulate complex physical processes in a neural network. Here, we employ BLNMs to encode the parametrized temporal dynamics of in silico 12-lead electrocardiograms (ECGs). BLNMs act as a geometry-specific surrogate model of cardiac function for fast and robust parameter estimation to match clinical ECGs in paediatric patients. Identifiability and trustworthiness of calibrated model parameters are assessed by sensitivity analysis and uncertainty quantification.