PLoS Computational Biology最新文献

Slow brain rhythms, for example during slow-wave sleep or pathological conditions like seizures and spreading depolarization, can be accompanied by oscillations in extracellular potassium concentration. Such slow brain rhythms typically have a lower frequency than tonic action-potential firing. They are assumed to arise from network-level mechanisms, involving synaptic interactions and delays, or from intrinsically bursting neurons. Neuronal burst generation is commonly attributed to ion channels with slow kinetics. Here, we explore an alternative mechanism generically available to all neurons with class I excitability. It is based on the interplay of fast-spiking voltage dynamics with a one-dimensional slow dynamics of the extracellular potassium concentration, mediated by the activity of the Na+/K+-ATPase. We use bifurcation analysis of the complete system as well as the slow-fast method to reveal that this coupling suffices to generate a hysteresis loop organized around a bistable region that emerges from a saddle-node loop bifurcation-a common feature of class I excitable neurons. Depending on the strength of the Na+/K+-ATPase, bursts are generated from pump-induced shearing the bifurcation structure, spiking is tonic, or cells are silenced via depolarization block. We suggest that transitions between these dynamics can result from disturbances in extracellular potassium regulation, such as glial malfunction or hypoxia affecting the Na+/K+-ATPase activity. The identified minimal mechanistic model outlining the sodium-potassium pump's generic contribution to burst dynamics can, therefore, contribute to a better mechanistic understanding of pathologies such as epilepsy syndromes and, potentially, inform therapeutic strategies.

Most essential cellular functions are performed by proteins assembled into larger complexes. Fluorescence Polarization Microscopy (FPM) is a powerful technique that goes beyond traditional imaging methods by allowing researchers to measure not only the localization of proteins within cells, but also their orientation or alignment within complexes or cellular structures. FPM can be easily integrated into standard widefield microscopes with the addition of a polarization modulator. However, the extensive image processing and analysis required to interpret the data have limited its widespread adoption. To overcome these challenges and enhance accessibility, we introduce OOPS (Object-Oriented Polarization Software), a MATLAB package for object-based analysis of FPM data. By combining flexible image segmentation and novel object-based analyses with a high-throughput FPM processing pipeline, OOPS empowers researchers to simultaneously study molecular order and orientation in individual biological structures; conduct population assessments based on morphological features, intensity statistics, and FPM measurements; and create publication-quality visualizations, all within a user-friendly graphical interface. Here, we demonstrate the power and versatility of our approach by applying OOPS to punctate and filamentous structures.

Cancer is a significant global health issue, with treatment challenges arising from intratumor heterogeneity. This heterogeneity stems mainly from somatic evolution, causing genetic diversity within the tumor, and phenotypic plasticity of tumor cells leading to reversible phenotypic changes. However, the interplay of both factors has not been rigorously investigated. Here, we examine the complex relationship between somatic evolution and phenotypic plasticity, explicitly focusing on the interplay between cell migration and proliferation. This type of phenotypic plasticity is essential in glioblastoma, the most aggressive form of brain tumor. We propose that somatic evolution alters the regulation of phenotypic plasticity in tumor cells, specifically the reaction to changes in the microenvironment. We study this hypothesis using a novel, spatially explicit model that tracks individual cells' phenotypic and genetic states. We assume cells change between migratory and proliferative states controlled by inherited and mutation-driven genotypes and the cells' microenvironment. We observe that cells at the tumor edge evolve to favor migration over proliferation and vice versa in the tumor bulk. Notably, different genetic configurations can result in this pattern of phenotypic heterogeneity. We analytically predict the outcome of the evolutionary process, showing that it depends on the tumor microenvironment. Synthetic tumors display varying levels of genetic and phenotypic heterogeneity, which we show are predictors of tumor recurrence time after treatment. Interestingly, higher phenotypic heterogeneity predicts poor treatment outcomes, unlike genetic heterogeneity. Our research offers a novel explanation for heterogeneous patterns of tumor recurrence in glioblastoma patients.

Human listeners have the ability to direct their attention to a single speaker in a multi-talker environment. The neural correlates of selective attention can be decoded from a single trial of electroencephalography (EEG) data. In this study, leveraging the source-reconstructed and anatomically-resolved EEG data as inputs, we sought to employ CNN as an interpretable model to uncover task-specific interactions between brain regions, rather than simply to utilize it as a black box decoder. To this end, our CNN model was specifically designed to learn pairwise interaction representations for 10 cortical regions from five-second inputs. By exclusively utilizing these features for decoding, our model was able to attain a median accuracy of 77.56% for within-participant and 65.14% for cross-participant classification. Through ablation analysis together with dissecting the features of the models and applying cluster analysis, we were able to discern the presence of alpha-band-dominated inter-hemisphere interactions, as well as alpha- and beta-band dominant interactions that were either hemisphere-specific or were characterized by a contrasting pattern between the right and left hemispheres. These interactions were more pronounced in parietal and central regions for within-participant decoding, but in parietal, central, and partly frontal regions for cross-participant decoding. These findings demonstrate that our CNN model can effectively utilize features known to be important in auditory attention tasks and suggest that the application of domain knowledge inspired CNNs on source-reconstructed EEG data can offer a novel computational framework for studying task-relevant brain interactions.

Although depolymerization of complex carbohydrates is a growth-limiting bottleneck for microbial decomposers, we still lack understanding about how the production of different types of extracellular enzymes affect individual microbes and in turn the performance of whole decomposer communities. In this work we use a theoretical model to evaluate the potential trade-offs faced by microorganisms in biopolymer decomposition which arise due to the varied biochemistry of different depolymerizing enzyme classes. We specifically consider two broad classes of depolymerizing extracellular enzymes, which are widespread across microbial taxa: exo-enzymes that cleave small units from the ends of polymer chains and endo-enzymes that act at random positions generating degradation products of varied sizes. Our results demonstrate a fundamental trade-off in the production of these enzymes, which is independent of system's complexity and which appears solely from the intrinsically different temporal depolymerization dynamics. As a consequence, specialists that produce either exo- or only endo-enzymes limit their growth to high or low substrate conditions, respectively. Conversely, generalists that produce both enzymes in an optimal ratio expand their niche and benefit from the synergy between the two enzymes. Finally, our results show that, in spatially-explicit environments, consortia composed of endo- and exo-specialists can only exist under oligotrophic conditions. In summary, our analysis demonstrates that the (evolutionary or ecological) selection of a depolymerization pathway will affect microbial fitness under low or high substrate conditions, with impacts on the ecological dynamics of microbial communities. It provides a possible explanation why many polysaccharide degraders in nature show the genetic potential to produce both of these enzyme classes.

In silico clinical trials (ISCTs) are an emerging method in modeling and simulation where medical interventions are evaluated using computational models of patients. ISCTs have the potential to provide cost-effective, time-efficient, and ethically favorable alternatives for evaluating the safety and effectiveness of medical devices. However, ensuring the credibility of ISCT results is a significant challenge. This paper aims to identify unique considerations for assessing the credibility of ISCTs and proposes an ISCT credibility assessment workflow based on recently published model assessment frameworks. First, we review various ISCTs described in the literature, carefully selected to showcase the range of methodological options available. These studies cover a wide variety of devices, reasons for conducting ISCTs, patient model generation approaches including subject-specific versus 'synthetic' virtual patients, complexity levels of devices and patient models, incorporation of clinician or clinical outcome models, and methods for integrating ISCT results with real-world clinical trials. We next discuss how verification, validation, and uncertainty quantification apply to ISCTs, considering the range of ISCT approaches identified. Based on our analysis, we then present a hierarchical workflow for assessing ISCT credibility, using a general credibility assessment framework recently published by the FDA's Center for Devices and Radiological Health. Overall, this work aims to promote standardization in ISCTs and contribute to the wider adoption and acceptance of ISCTs as a reliable tool for evaluating medical devices.

The rapid development of vaccines against SARS-CoV-2 altered the course of the COVID-19 pandemic. In most countries, vaccinations were initially targeted at high-risk populations, including older individuals and healthcare workers. Now, despite substantial infection- and vaccine-induced immunity in host populations worldwide, waning immunity and the emergence of novel variants continue to cause significant waves of infection and disease. Policy makers must determine how to deploy booster vaccinations, particularly when constraints in vaccine supply, delivery and cost mean that booster vaccines cannot be administered to everyone. A key question is therefore whether older individuals should again be prioritised for vaccination, or whether alternative strategies (e.g. offering booster vaccines to the individuals who have most contacts with others and therefore drive infection) can instead offer indirect protection to older individuals. Here, we use mathematical modelling to address this question, considering SARS-CoV-2 transmission in a range of countries with different socio-economic backgrounds. We show that the population structures of different countries can have a pronounced effect on the impact of booster vaccination, even when identical booster vaccination targeting strategies are adopted. However, under the assumed transmission model, prioritising older individuals for booster vaccination consistently leads to the most favourable public health outcomes in every setting considered. This remains true for a range of assumptions about booster vaccine supply and timing, and for different assumed policy objectives of booster vaccination.

Human behaviors have non-negligible impacts on spread of contagious disease. For instance, large-scale gathering and high mobility of population could lead to accelerated disease transmission, while public behavioral changes in response to pandemics may effectively reduce contacts and suppress the peak of the outbreak. In order to understand how spatial characteristics like population mobility and clustering interplay with epidemic outbreaks, we formulate a stochastic-statistical environment-epidemic dynamic system (SEEDS) via an agent-based biased random walk model on a two-dimensional lattice. The "popularity" and "awareness" variables are taken into consideration to capture human natural and preventive behavioral factors, which are assumed to guide and bias agent movement in a combined way. It is found that the presence of the spatial heterogeneity, like social influence locality and spatial clustering induced by self-aggregation, potentially suppresses the contacts between agents and consequently flats the epidemic curve. Surprisedly, disease responses might not necessarily reduce the susceptibility of informed individuals and even aggravate disease outbreak if each individual responds independently upon their awareness. The disease control is achieved effectively only if there are coordinated public-health interventions and public compliance to these measures. Therefore, our model may be useful for quantitative evaluations of a variety of public-health policies.

Single-cell RNA sequencing (scRNA-seq) has emerged as a powerful tool in genomics research, enabling the analysis of gene expression at the individual cell level. However, scRNA-seq data often suffer from a high rate of dropouts, where certain genes fail to be detected in specific cells due to technical limitations. This missing data can introduce biases and hinder downstream analysis. To overcome this challenge, the development of effective imputation methods has become crucial in the field of scRNA-seq data analysis. Here, we propose an imputation method based on robust and non-negative matrix factorization (scRNMF). Instead of other matrix factorization algorithms, scRNMF integrates two loss functions: L2 loss and C-loss. The L2 loss function is highly sensitive to outliers, which can introduce substantial errors. We utilize the C-loss function when dealing with zero values in the raw data. The primary advantage of the C-loss function is that it imposes a smaller punishment for larger errors, which results in more robust factorization when handling outliers. Various datasets of different sizes and zero rates are used to evaluate the performance of scRNMF against other state-of-the-art methods. Our method demonstrates its power and stability as a tool for imputation of scRNA-seq data.