Oscillations in Neuronal Activity: A Neuron-Centered Spatiotemporal Model of the Unfolded Protein Response in Prion Diseases.

Many neurodegenerative diseases (NDs) are characterized by the slow spatial spread of toxic protein species in the brain. The toxic proteins can induce neuronal stress, triggering the Unfolded Protein Response (UPR), which slows or stops protein translation and can indirectly reduce the toxic load. However, the UPR may also trigger processes leading to apoptotic cell death and the UPR is implicated in the progression of several NDs. In this paper, we develop a novel mathematical model to describe the spatiotemporal dynamics of the UPR mechanism for prion diseases. Our model is centered around a single neuron, with representative proteins P (healthy) and S (toxic) interacting with heterodimer dynamics (S interacts with P to form two S's). The model takes the form of a coupled system of nonlinear reaction-diffusion equations with a delayed, nonlinear flux for P (delay from the UPR). Through the delay, we find parameter regimes that exhibit oscillations in the P- and S-protein levels. We find that oscillations are more pronounced when the S-clearance rate and S-diffusivity are small in comparison to the P-clearance rate and P-diffusivity, respectively. The oscillations become more pronounced as delays in initiating the UPR increase. We also consider quasi-realistic clinical parameters to understand how possible drug therapies can alter the course of a prion disease. We find that decreasing the production of P, decreasing the recruitment rate, increasing the diffusivity of S, increasing the UPR S-threshold, and increasing the S clearance rate appear to be the most powerful modifications to reduce the mean UPR intensity and potentially moderate the disease progression.

Networks of necessity: Simulating COVID-19 mitigation strategies for disabled people and their caregivers.

A major strategy to prevent the spread of COVID-19 is the limiting of in-person contacts. However, limiting contacts is impractical or impossible for the many disabled people who do not live in care facilities but still require caregivers to assist them with activities of daily living. We seek to determine which interventions can best prevent infections of disabled people and their caregivers. To accomplish this, we simulate COVID-19 transmission with a compartmental model that includes susceptible, exposed, asymptomatic, symptomatically ill, hospitalized, and removed/recovered individuals. The networks on which we simulate disease spread incorporate heterogeneity in the risk levels of different types of interactions, time-dependent lockdown and reopening measures, and interaction distributions for four different groups (caregivers, disabled people, essential workers, and the general population). Of these groups, we find that the probability of becoming infected is largest for caregivers and second largest for disabled people. Consistent with this finding, our analysis of network structure illustrates that caregivers have the largest modal eigenvector centrality of the four groups. We find that two interventions-contact-limiting by all groups and mask-wearing by disabled people and caregivers-most reduce the number of infections in disabled and caregiver populations. We also test which group of people spreads COVID-19 most readily by seeding infections in a subset of each group and comparing the total number of infections as the disease spreads. We find that caregivers are the most potent spreaders of COVID-19, particularly to other caregivers and to disabled people. We test where to use limited infection-blocking vaccine doses most effectively and find that (1) vaccinating caregivers better protects disabled people from infection than vaccinating the general population or essential workers and that (2) vaccinating caregivers protects disabled people from infection about as effectively as vaccinating disabled people themselves. Our results highlight the potential effectiveness of mask-wearing, contact-limiting throughout society, and strategic vaccination for limiting the exposure of disabled people and their caregivers to COVID-19.

From reaction kinetics to dementia: A simple dimer model of Alzheimer's disease etiology.

Oligomers of the amyloid ß-protein (Aß) have been implicated in the pathogenesis of Alzheimer's disease (AD) through their toxicity towards neurons. Understanding the process of oligomerization may contribute to the development of therapeutic agents, but this has been difficult due to the complexity of oligomerization and the metastability of the oligomers thus formed. To understand the kinetics of oligomer formation, and how that relates to the progression of AD, we developed models of the oligomerization process. Here, we use experimental data from cell viability assays and proxies for rate constants involved in monomer-dimer-trimer kinetics to develop a simple mathematical model linking Aß assembly to oligomer-induced neuronal degeneration. This model recapitulates the rapid growth of disease incidence with age. It does so through incorporation of age-dependent changes in rates of Aß monomer production and elimination. The model also describes clinical progression in genetic forms of AD (e.g., Down's syndrome), changes in hippocampal volume, AD risk after traumatic brain injury, and spatial spreading of the disease due to foci in which Aß production is elevated. Continued incorporation of clinical and basic science data into the current model will make it an increasingly relevant model system for doing theoretical calculations that are not feasible in biological systems. In addition, terms in the model that have particularly large effects are likely to be especially useful therapeutic targets.

Functional Kernel Density Estimation: Point and Fourier Approaches to Time Series Anomaly Detection.

We present an unsupervised method to detect anomalous time series among a collection of time series. To do so, we extend traditional Kernel Density Estimation for estimating probability distributions in Euclidean space to Hilbert spaces. The estimated probability densities we derive can be obtained formally through treating each series as a point in a Hilbert space, placing a kernel at those points, and summing the kernels (a "point approach"), or through using Kernel Density Estimation to approximate the distributions of Fourier mode coefficients to infer a probability density (a "Fourier approach"). We refer to these approaches as Functional Kernel Density Estimation for Anomaly Detection as they both yield functionals that can score a time series for how anomalous it is. Both methods naturally handle missing data and apply to a variety of settings, performing well when compared with an outlyingness score derived from a boxplot method for functional data, with a Principal Component Analysis approach for functional data, and with the Functional Isolation Forest method. We illustrate the use of the proposed methods with aviation safety report data from the International Air Transport Association (IATA).