Minimal Mechanisms of Microtubule Length Regulation in Living Cells.

The microtubule cytoskeleton is responsible for sustained, long-range intracellular transport of mRNAs, proteins, and organelles in neurons. Neuronal microtubules must be stable enough to ensure reliable transport, but they also undergo dynamic instability, as their plus and minus ends continuously switch between growth and shrinking. This process allows for continuous rebuilding of the cytoskeleton and for flexibility in injury settings. Motivated by in vivo experimental data on microtubule behavior in Drosophila neurons, we propose a mathematical model of dendritic microtubule dynamics, with a focus on understanding microtubule length, velocity, and state-duration distributions. We find that limitations on microtubule growth phases are needed for realistic dynamics, but the type of limiting mechanism leads to qualitatively different responses to plausible experimental perturbations. We therefore propose and investigate two minimally-complex length-limiting factors: limitation due to resource (tubulin) constraints and limitation due to catastrophe of large-length microtubules. We combine simulations of a detailed stochastic model with steady-state analysis of a mean-field ordinary differential equations model to map out qualitatively distinct parameter regimes. This provides a basis for predicting changes in microtubule dynamics, tubulin allocation, and the turnover rate of tubulin within microtubules in different experimental environments.

Parameter Identifiability in PDE Models of Fluorescence Recovery After Photobleaching.

Identifying unique parameters for mathematical models describing biological data can be challenging and often impossible. Parameter identifiability for partial differential equations models in cell biology is especially difficult given that many established in vivo measurements of protein dynamics average out the spatial dimensions. Here, we are motivated by recent experiments on the binding dynamics of the RNA-binding protein PTBP3 in RNP granules of frog oocytes based on fluorescence recovery after photobleaching (FRAP) measurements. FRAP is a widely-used experimental technique for probing protein dynamics in living cells, and is often modeled using simple reaction-diffusion models of the protein dynamics. We show that current methods of structural and practical parameter identifiability provide limited insights into identifiability of kinetic parameters for these PDE models and spatially-averaged FRAP data. We thus propose a pipeline for assessing parameter identifiability and for learning parameter combinations based on re-parametrization and profile likelihoods analysis. We show that this method is able to recover parameter combinations for synthetic FRAP datasets and investigate its application to real experimental data.

Simulated actin reorganization mediated by motor proteins.

Cortical actin networks are highly dynamic and play critical roles in shaping the mechanical properties of cells. The actin cytoskeleton undergoes significant reorganization in many different contexts, including during directed cell migration and over the course of the cell cycle, when cortical actin can transition between different configurations such as open patched meshworks, homogeneous distributions, and aligned bundles. Several types of myosin motor proteins, characterized by different kinetic parameters, have been involved in this reorganization of actin filaments. Given the limitations in studying the interactions of actin with myosin in vivo, we propose stochastic agent-based models and develop a set of data analysis measures to assess how myosin motor proteins mediate various actin organizations. In particular, we identify individual motor parameters, such as motor binding rate and step size, that generate actin networks with different levels of contractility and different patterns of myosin motor localization, which have previously been observed experimentally. In simulations where two motor populations with distinct kinetic parameters interact with the same actin network, we find that motors may act in a complementary way, by tuning the actin network organization, or in an antagonistic way, where one motor emerges as dominant. This modeling and data analysis framework also uncovers parameter regimes where spatial segregation between motor populations is achieved. By allowing for changes in kinetic rates during the actin-myosin dynamic simulations, our work suggests that certain actin-myosin organizations may require additional regulation beyond mediation by motor proteins in order to reconfigure the cytoskeleton network on experimentally-observed timescales.

Identification of approximate symmetries in biological development.

Virtually all forms of life, from single-cell eukaryotes to complex, highly differentiated multicellular organisms, exhibit a property referred to as symmetry. However, precise measures of symmetry are often difficult to formulate and apply in a meaningful way to biological systems, where symmetries and asymmetries can be dynamic and transient, or be visually apparent but not reliably quantifiable using standard measures from mathematics and physics. Here, we present and illustrate a novel measure that draws on concepts from information theory to quantify the degree of symmetry, enabling the identification of approximate symmetries that may be present in a pattern or a biological image. We apply the measure to rotation, reflection and translation symmetries in patterns produced by a Turing model, as well as natural objects (algae, flowers and leaves). This method of symmetry quantification is unbiased and rigorous, and requires minimal manual processing compared to alternative measures. The proposed method is therefore a useful tool for comparison and identification of symmetries in biological systems, with potential future applications to symmetries that arise during development, as observed in vivo or as produced by mathematical models. This article is part of the theme issue 'Recent progress and open frontiers in Turing's theory of morphogenesis'.

Topological Data Analysis Approaches to Uncovering the Timing of Ring Structure Onset in Filamentous Networks.

In developmental biology as well as in other biological systems, emerging structure and organization can be captured using time-series data of protein locations. In analyzing this time-dependent data, it is a common challenge not only to determine whether topological features emerge, but also to identify the timing of their formation. For instance, in most cells, actin filaments interact with myosin motor proteins and organize into polymer networks and higher-order structures. Ring channels are examples of such structures that maintain constant diameters over time and play key roles in processes such as cell division, development, and wound healing. Given the limitations in studying interactions of actin with myosin in vivo, we generate time-series data of protein polymer interactions in cells using complex agent-based models. Since the data has a filamentous structure, we propose sampling along the actin filaments and analyzing the topological structure of the resulting point cloud at each time. Building on existing tools from persistent homology, we develop a topological data analysis (TDA) method that assesses effective ring generation in this dynamic data. This method connects topological features through time in a path that corresponds to emergence of organization in the data. In this work, we also propose methods for assessing whether the topological features of interest are significant and thus whether they contribute to the formation of an emerging hole (ring channel) in the simulated protein interactions. In particular, we use the MEDYAN simulation platform to show that this technique can distinguish between the actin cytoskeleton organization resulting from distinct motor protein binding parameters.

Renewal Reward Perspective on Linear Switching Diffusion Systems in Models of Intracellular Transport.

In many biological systems, the movement of individual agents is characterized having multiple qualitatively distinct behaviors that arise from a variety of biophysical states. For example, in cells the movement of vesicles, organelles, and other intracellular cargo is affected by their binding to and unbinding from cytoskeletal filaments such as microtubules through molecular motor proteins. A typical goal of theoretical or numerical analysis of models of such systems is to investigate effective transport properties and their dependence on model parameters. While the effective velocity of particles undergoing switching diffusion dynamics is often easily characterized in terms of the long-time fraction of time that particles spend in each state, the calculation of the effective diffusivity is more complicated because it cannot be expressed simply in terms of a statistical average of the particle transport state at one moment of time. However, it is common that these systems are regenerative, in the sense that they can be decomposed into independent cycles marked by returns to a base state. Using decompositions of this kind, we calculate effective transport properties by computing the moments of the dynamics within each cycle and then applying renewal reward theory. This method provides a useful alternative large-time analysis to direct homogenization for linear advection-reaction-diffusion partial differential equation models. Moreover, it applies to a general class of semi-Markov processes and certain stochastic differential equations that arise in models of intracellular transport. Applications of the proposed renewal reward framework are illustrated for several case studies such as mRNA transport in developing oocytes and processive cargo movement by teams of molecular motor proteins.

Cardiovascular regulation in response to multiple hemorrhages: analysis and parameter estimation.

Mathematical models can provide useful insights explaining behavior observed in experimental data; however, rigorous analysis is needed to select a subset of model parameters that can be informed by available data. Here we present a method to estimate an identifiable set of parameters based on baseline left ventricular pressure and volume time series data. From this identifiable subset, we then select, based on current understanding of cardiovascular control, parameters that vary in time in response to blood withdrawal, and estimate these parameters over a series of blood withdrawals. These time-varying parameters are first estimated using piecewise linear splines minimizing the mean squared error between measured and computed left ventricular pressure and volume data over four consecutive blood withdrawals. As a final step, the trends in these splines are fit with empirical functional expressions selected to describe cardiovascular regulation during blood withdrawal. Our analysis at baseline found parameters representing timing of cardiac contraction, systemic vascular resistance, and cardiac contractility to be identifiable. Of these parameters, vascular resistance and cardiac contractility were varied in time. Data used for this study were measured in a control Sprague-Dawley rat. To our knowledge, this is the first study to analyze the response to multiple blood withdrawals both experimentally and theoretically, as most previous studies focus on analyzing the response to one large blood withdrawal. Results show that during each blood withdrawal both systemic vascular resistance and contractility decrease acutely and partially recover, and they decrease chronically across the series of blood withdrawals.

Model reconstruction from temporal data for coupled oscillator networks.

In a complex system, the interactions between individual agents often lead to emergent collective behavior such as spontaneous synchronization, swarming, and pattern formation. Beyond the intrinsic properties of the agents, the topology of the network of interactions can have a dramatic influence over the dynamics. In many studies, researchers start with a specific model for both the intrinsic dynamics of each agent and the interaction network and attempt to learn about the dynamics of the model. Here, we consider the inverse problem: given data from a system, can one learn about the model and the underlying network? We investigate arbitrary networks of coupled phase oscillators that can exhibit both synchronous and asynchronous dynamics. We demonstrate that, given sufficient observational data on the transient evolution of each oscillator, machine learning can reconstruct the interaction network and identify the intrinsic dynamics.

Modeling microtubule-based transport and anchoring of mRNA.

Localization of messenger RNA (mRNA) at the vegetal cortex plays an important role in the early development of Xenopus laevis oocytes. While it is known that molecular motors are responsible for the transport of mRNA cargo along microtubules to the cortex, the mechanisms of localization remain unclear. We model cargo transport along microtubules using partial differential equations with spatially-dependent rates. A theoretical analysis of reduced versions of our model predicts effective velocity and diffusion rates for the cargo and shows that randomness of microtubule networks enhances effective transport. A more complex model using parameters estimated from fluorescence microscopy data reproduces the spatial and timescales of mRNA localization observed in Xenopus oocytes, corroborates experimental hypotheses that anchoring may be necessary to achieve complete localization, and shows that anchoring of mRNA complexes actively transported to the cortex is most effective in achieving robust accumulation at the cortex.

Analysis of Active Transport by Fluorescence Recovery after Photobleaching.

Fluorescence recovery after photobleaching (FRAP) is a well-established experimental technique to study binding and diffusion of molecules in cells. Although a large number of analytical and numerical models have been developed to extract binding and diffusion rates from FRAP recovery curves, active transport of molecules is typically not included in the existing models that are used to estimate these rates. Here we present a validated numerical method for estimating diffusion, binding/unbinding rates, and active transport velocities using FRAP data that captures intracellular dynamics through partial differential equation models. We apply these methods to transport and localization of mRNA molecules in Xenopus laevis oocytes, where active transport processes are essential to generate developmental polarity. By providing estimates of the effective velocities and diffusion, as well as expected run times and lengths, this approach can help quantify dynamical properties of localizing and nonlocalizing RNA. Our results confirm the distinct transport dynamics in different regions of the cytoplasm, and suggest that RNA movement in both the animal and vegetal directions may influence the timescale of RNA localization in Xenopus oocytes. We also show that model initial conditions extracted from FRAP postbleach intensities prevent underestimation of diffusion, which can arise from the instantaneous bleaching assumption. The numerical and modeling approach presented here to estimate parameters using FRAP recovery data is a broadly applicable tool for systems where intracellular transport is a key molecular mechanism.

Minimal Mechanisms of Microtubule Length Regulation in Living Cells.

The microtubule cytoskeleton is responsible for sustained, long-range intracellular transport of mRNAs, proteins, and organelles in neurons. Neuronal microtubules must be stable enough to ensure reliable transport, but they also undergo dynamic instability, as their plus and minus ends continuously switch between growth and shrinking. This process allows for continuous rebuilding of the cytoskeleton and for flexibility in injury settings. Motivated by in vivo experimental data on microtubule behavior in Drosophila neurons, we propose a mathematical model of dendritic microtubule dynamics, with a focus on understanding microtubule length, velocity, and state-duration distributions. We find that limitations on microtubule growth phases are needed for realistic dynamics, but the type of limiting mechanism leads to qualitatively different responses to plausible experimental perturbations. We therefore propose and investigate two minimally-complex length-limiting factors: limitation due to resource (tubulin) constraints and limitation due to catastrophe of large-length microtubules. We combine simulations of a detailed stochastic model with steady-state analysis of a mean-field ordinary differential equations model to map out qualitatively distinct parameter regimes. This provides a basis for predicting changes in microtubule dynamics, tubulin allocation, and the turnover rate of tubulin within microtubules in different experimental environments.

Characterizing emerging features in cell dynamics using topological data analysis methods.

Filament-motor interactions inside cells play essential roles in many developmental as well as other biological processes. For instance, actin-myosin interactions drive the emergence or closure of ring channel structures during wound healing or dorsal closure. These dynamic protein interactions and the resulting protein organization lead to rich time-series data generated by using fluorescence imaging experiments or by simulating realistic stochastic models. We propose methods based on topological data analysis to track topological features through time in cell biology data consisting of point clouds or binary images. The framework proposed here is based on computing the persistent homology of the data at each time point and on connecting topological features through time using established distance metrics between topological summaries. The methods retain aspects of monomer identity when analyzing significant features in filamentous structure data, and capture the overall closure dynamics when assessing the organization of multiple ring structures through time. Using applications of these techniques to experimental data, we show that the proposed methods can describe features of the emergent dynamics and quantitatively distinguish between control and perturbation experiments.

How do classroom-turnover times depend on lecture-hall size?

Academic spaces in colleges and universities span classrooms for 10 students to lecture halls that hold over 600 people. During the break between consecutive classes, students from the first class must leave and the new class must find their desks, regardless of whether the room holds 10 or 600 people. Here we address the question of how the size of large lecture halls affects classroom-turnover times, focusing on non-emergency settings. By adapting the established social-force model, we treat students as individuals who interact and move through classrooms to reach their destinations. We find that social interactions and the separation time between consecutive classes strongly influence how long it takes entering students to reach their desks, and that these effects are more pronounced in larger lecture halls. While the median time that individual students must travel increases with decreased separation time, we find that shorter separation times lead to shorter classroom-turnover times overall. This suggests that the effects of scheduling gaps and lecture-hall size on classroom dynamics depends on the perspective-individual student or whole class-that one chooses to take.

Influenza spread on context-specific networks lifted from interaction-based diary data.

Studying the spread of infections is an important tool in limiting or preventing future outbreaks. A first step in understanding disease dynamics is constructing networks that reproduce features of real-world interactions. In this paper, we generate networks that maintain some features of the partial interaction networks that were recorded in an existing diary-based survey at the University of Warwick. To preserve realistic structure in our artificial networks, we use a context-specific approach. In particular, we propose different algorithms for producing larger home, work and social networks. Our networks are able to maintain much of the interaction structure in the original diary-based survey and provide a means of accounting for the interactions of survey participants with non-participants. Simulating a discrete susceptible-infected-recovered model on the full network produces epidemic behaviour which shares characteristics with previous influenza seasons. Our approach allows us to explore how disease transmission and dynamic responses to infection differ depending on interaction context. We find that, while social interactions may be the first to be reduced after influenza infection, limiting work and school encounters may be significantly more effective in controlling the overall severity of the epidemic.

A mechanism for neurofilament transport acceleration through nodes of Ranvier.

Neurofilaments are abundant space-filling cytoskeletal polymers in axons that are transported along microtubule tracks. Neurofilament transport is accelerated at nodes of Ranvier, where axons are locally constricted. Strikingly, these constrictions are accompanied by sharp decreases in neurofilament number, no decreases in microtubule number, and increases in the packing density of these polymers, which collectively bring nodal neurofilaments closer to their microtubule tracks. We hypothesize that this leads to an increase in the proportion of time that the filaments spend moving and that this can explain the local acceleration. To test this, we developed a stochastic model of neurofilament transport that tracks their number, kinetic state, and proximity to nearby microtubules in space and time. The model assumes that the probability of a neurofilament moving is dependent on its distance from the nearest available microtubule track. Taking into account experimentally reported numbers and densities for neurofilaments and microtubules in nodes and internodes, we show that the model is sufficient to explain the local acceleration of neurofilaments within nodes of Ranvier. This suggests that proximity to microtubule tracks may be a key regulator of neurofilament transport in axons, which has implications for the mechanism of neurofilament accumulation in development and disease.

JUSTFAIR: Judicial System Transparency through Federal Archive Inferred Records.

In the United States, the public has a constitutional right to access criminal trial proceedings. In practice, it can be difficult or impossible for the public to exercise this right. We present JUSTFAIR: Judicial System Transparency through Federal Archive Inferred Records, a database of criminal sentencing decisions made in federal district courts. We have compiled this data set from public sources including the United States Sentencing Commission, the Federal Judicial Center, the Public Access to Court Electronic Records system, and Wikipedia. With nearly 600,000 records from the years 2001-2018, JUSTFAIR is the first large scale, free, public database that links information about defendants and their demographic characteristics with information about their federal crimes, their sentences, and, crucially, the identity of the sentencing judge.