Age‐related central neurodegenerative diseases, such as Alzheimer's and Parkinson's disease, are a rising public health concern and have been plagued by repeated drug development failures. The complex nature and poor mechanistic understanding of the etiology of neurodegenerative diseases has hindered the discovery and development of effective disease‐modifying therapeutics. Quantitative systems pharmacology models of neurodegeneration diseases may be useful tools to enhance the understanding of pharmacological intervention strategies and to reduce drug attrition rates. Due to the similarities in pathophysiological mechanisms across neurodegenerative diseases, especially at the cellular and molecular levels, we envision the possibility of structural components that are conserved across models of neurodegenerative diseases. Conserved structural submodels can be viewed as building blocks that are pieced together alongside unique disease components to construct quantitative systems pharmacology (QSP) models of neurodegenerative diseases. Model parameterization would likely be different between the different types of neurodegenerative diseases as well as individual patients. Formulating our mechanistic understanding of neurodegenerative pathophysiology as a mathematical model could aid in the identification and prioritization of drug targets and combinatorial treatment strategies, evaluate the role of patient characteristics on disease progression and therapeutic response, and serve as a central repository of knowledge. Here, we provide a background on neurodegenerative diseases, highlight hallmarks of neurodegeneration, and summarize previous QSP models of neurodegenerative diseases.

This tutorial illustrates the use of data assimilation algorithms to estimate unobserved variables and unknown parameters of conductance-based neuronal models. Modern data assimilation (DA) techniques are widely used in climate science and weather prediction, but have only recently begun to be applied in neuroscience. The two main classes of DA techniques are sequential methods and variational methods. We provide computer code implementing basic versions of a method from each class, the Unscented Kalman Filter and 4D-Var, and demonstrate how to use these algorithms to infer several parameters of the Morris–Lecar model from a single voltage trace. Depending on parameters, the Morris–Lecar model exhibits qualitatively different types of neuronal excitability due to changes in the underlying bifurcation structure. We show that when presented with voltage traces from each of the various excitability regimes, the DA methods can identify parameter sets that produce the correct bifurcation structure even with initial parameter guesses that correspond to a different excitability regime. This demonstrates the ability of DA techniques to perform nonlinear state and parameter estimation and introduces the geometric structure of inferred models as a novel qualitative measure of estimation success. We conclude by discussing extensions of these DA algorithms that have appeared in the neuroscience literature.Electronic Supplementary MaterialThe online version of this article (10.1186/s13408-018-0066-8) contains supplementary material.

Circadian transcriptional timekeepers in pacemaker neurons drive profound daily rhythms in sleep and wake. Here we reveal a molecular pathway that links core transcriptional oscillators to neuronal and behavioral rhythms. Using two independent genetic screens, we identified mutants of Transport and Golgi organization 10 (Tango10) with poor behavioral rhythmicity. Tango10 expression in pacemaker neurons expressing the neuropeptide PIGMENT-DISPERSING FACTOR (PDF) is required for robust rhythms. Loss of Tango10 results in elevated PDF accumulation in nerve terminals even in mutants lacking a functional core clock. TANGO10 protein itself is rhythmically expressed in PDF terminals. Mass spectrometry of TANGO10 complexes reveals interactions with the E3 ubiquitin ligase CULLIN 3 (CUL3). CUL3 depletion phenocopies Tango10 mutant effects on PDF even in the absence of the core clock gene timeless. Patch clamp electrophysiology in Tango10 mutant neurons demonstrates elevated spontaneous firing potentially due to reduced voltage-gated Shaker-like potassium currents. We propose that Tango10/Cul3 transduces molecular oscillations from the core clock to neuropeptide release important for behavioral rhythms.

Circadian rhythms in mammals are orchestrated by a central clock within the suprachiasmatic nuclei (SCN). Our understanding of the electrophysiological basis of SCN activity comes overwhelmingly from a small number of nocturnal rodent species, and the extent to which these are retained in day-active animals remains unclear. Here, we recorded the spontaneous and evoked electrical activity of single SCN neurons in the diurnal rodent Rhabdomys pumilio, and developed cutting-edge data assimilation and mathematical modeling approaches to uncover the underlying ionic mechanisms. As in nocturnal rodents, R. pumilio SCN neurons were more excited during daytime hours. By contrast, the evoked activity of R. pumilio neurons included a prominent suppressive response that is not present in the SCN of nocturnal rodents. Our modeling revealed and subsequent experiments confirmed transient subthreshold A-type potassium channels as the primary determinant of this response, and suggest a key role for this ionic mechanism in optimizing SCN function to accommodate R. pumilio’s diurnal niche.

The standard of care treatment for neovascular age-related macular degeneration, delivered as ocular injection, is based on anti-vascular endothelial growth factor (anti-VEGF). The course of treatment may need to be modified quickly for certain patients based on their response. Models that track both the concentration and the response to an anti-VEGF treatment are presented. The specific focus is to assess the existence of analytical solutions for the different types of models. Both an ODE-based model and a map-based model illustrate the dependence of the solution on various biological parameters and allow the measurement of patient-specific parameters from experimental data. A PDE-based model incorporates diffusive effects. The results are consistent with observed values, and could provide a framework for practitioners to understand the effect of the therapy on the progression of the disease in both responsive and non-responsive patients.

No abstract