A most wanted list of conserved microbial protein families with no known domains

The number and proportion of genes with no known function are growing rapidly. To quantify this phenomenon and provide criteria for prioritizing genes for functional characterization, we developed a bioinformatics pipeline that identifies robustly defined protein families with no annotated domains, ranks these with respect to phylogenetic breadth, and identifies them in metagenomics data. We applied this approach to 271 965 protein families from the SFams database and discovered many with no functional annotation, including >118 000 families lacking any known protein domain. From these, we prioritized 6 668 conserved protein families with at least three sequences from organisms in at least two distinct classes. These Function Unknown Families (FUnkFams) are present in Tara Oceans Expedition and Human Microbiome Project metagenomes, with distributions associated with sampling environment. Our findings highlight the extent of functional novelty in sequence databases and establish an approach for creating a "most wanted" list of genes to prioritize for further characterization.

Genomic region detection via Spatial Convex Clustering

Several modern genomic technologies, such as DNA-Methylation arrays, measure spatially registered probes that number in the hundreds of thousands across multiple chromosomes. The measured probes are by themselves less interesting scientifically; instead scientists seek to discover biologically interpretable genomic regions comprised of contiguous groups of probes which may act as biomarkers of disease or serve as a dimension-reducing pre-processing step for downstream analyses. In this paper, we introduce an unsupervised feature learning technique which maps technological units (probes) to biological units (genomic regions) that are common across all subjects. We use ideas from fusion penalties and convex clustering to introduce a method for Spatial Convex Clustering, or SpaCC. Our method is specifically tailored to detecting multi-subject regions of methylation, but we also test our approach on the well-studied problem of detecting segments of copy number variation. We formulate our method as a convex optimization problem, develop a massively parallelizable algorithm to find its solution, and introduce automated approaches for handling missing values and determining tuning parameters. Through simulation studies based on real methylation and copy number variation data, we show that SpaCC exhibits significant performance gains relative to existing methods. Finally, we illustrate SpaCC's advantages as a pre-processing technique that reduces large-scale genomics data into a smaller number of genomic regions through several cancer epigenetics case studies on subtype discovery, network estimation, and epigenetic-wide association.

Functional connectomics from neural dynamics: probabilistic graphical models for neuronal network of Caenorhabditis elegans

We propose an approach to represent neuronal network dynamics as a probabilistic graphical model (PGM). To construct the PGM, we collect time series of neuronal responses produced by the neuronal network and use singular value decomposition to obtain a low-dimensional projection of the time-series data. We then extract dominant patterns from the projections to get pairwise dependency information and create a graphical model for the full network. The outcome model is a functional connectome that captures how stimuli propagate through the network and thus represents causal dependencies between neurons and stimuli. We apply our methodology to a model of the Caenorhabditis elegans somatic nervous system to validate and show an example of our approach. The structure and dynamics of the C. elegans nervous system are well studied and a model that generates neuronal responses is available. The resulting PGM enables us to obtain and verify underlying neuronal pathways for known behavioural scenarios and detect possible pathways for novel scenarios.This article is part of a discussion meeting issue 'Connectome to behaviour: modelling C. elegans at cellular resolution'.

Data Assimilation Methods for Neuronal State and Parameter Estimation

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.

The role of intracellular signaling in the stripe formation in engineered Escherichia coli populations

Recent experiments showed that engineered Escherichia coli colonies grow and self-organize into periodic stripes with high and low cell densities in semi-solid agar. The stripes develop sequentially behind a radially propagating colony front, similar to the formation of many other periodic patterns in nature. These bacteria were created by genetically coupling the intracellular chemotaxis pathway of wild-type cells with a quorum sensing module through the protein CheZ. In this paper, we develop multiscale models to investigate how this intracellular pathway affects stripe formation. We first develop a detailed hybrid model that treats each cell as an individual particle and incorporates intracellular signaling via an internal ODE system. To overcome the computational cost of the hybrid model caused by the large number of cells involved, we next derive a mean-field PDE model from the hybrid model using asymptotic analysis. We show that this analysis is justified by the tight agreement between the PDE model and the hybrid model in 1D simulations. Numerical simulations of the PDE model in 2D with radial symmetry agree with experimental data semi-quantitatively. Finally, we use the PDE model to make a number of testable predictions on how the stripe patterns depend on cell-level parameters, including cell speed, cell doubling time and the turnover rate of intracellular CheZ.

Unintended consequences and the paradox of control: Management of emerging pathogens with age-specific virulence

We project forward total Zika virus disease (ZVD) under varying hazards of infection and consider how the age distribution of disease burden varies between these scenarios. Pathogens with age structured disease outcomes, such as rubella and Zika virus, require that management decisions consider their impact not only on total disease incidence but also on distribution of disease burden within a population. Some situations exhibit a "paradox of control" in which reductions of overall transmission decrease the total incidence but increase the incidence of severe disease. This happens because of corresponding increases in the average age of infection. Beginning with the current population structure and demographic rates of Brazil, we project forward total ZVD burden as measured by cases occurring in pregnant women and document the scenarios under which a paradox of control for ZVD management emerges. We conclude that while a paradox of control can occur for ZVD, the higher total costs from increasing the average age of infection will only be realized after several decades and vanish under conservative discounting of future costs. This indicates that managers faced with an emerging pathogen are justified to prioritize current disease incidence over potential increases in severe disease outcomes in the endemic state.

Receptor recharge time drastically reduces the number of captured particles

Many diverse biological systems are described by randomly moving particles that can be captured by traps in their environment. Examples include neurotransmitters diffusing in the synaptic cleft before binding to receptors and prey roaming an environment before capture by predators. In most cases, the traps cannot capture particles continuously. Rather, each trap must wait a transitory “recharge” time after capturing a particle before additional captures. This recharge time is often overlooked. In the case of instant recharge, the average number of particles captured before they escape grows linearly in the total number of particles. In stark contrast, we prove that for any nonzero recharge time, the average number of captured particles grows at most logarithmically in the total particle number. This is a fundamental effect of recharge, as it holds under very general assumptions on particle motion and spatial domain. Furthermore, we characterize the parameter regime in which a given recharge time will dramatically affect a system, allowing researchers to easily verify if they need to account for recharge in their specific system. Finally, we consider a few examples, including a neural system in which recharge reduces neurotransmitter bindings by several orders of magnitude.

Consider particles that are released into an environment (think diffusing molecules or plankton), and suppose that there are traps in the environment. How many particles will be captured by the traps before they escape? In a standard model, the number of captured particles is proportional to the initial number released. In this paper, we show that for a more realistic model of a trap (one in which traps must recharge after every capture), the number of captures is proportional to the logarithm of the initial number released. That means that if 10⁶ particles are released, only about 6 will be captured. We prove this result mathematically, and then consider a number of applications, including neuronal synapses and ambush predators.

Computable performance guarantees for compressed sensing matrices

The null space condition for ℓ1 minimization in compressed sensing is a necessary and sufficient condition on the sensing matrices under which a sparse signal can be uniquely recovered from the observation data via ℓ1 minimization. However, verifying the null space condition is known to be computationally challenging. Most of the existing methods can provide only upper and lower bounds on the proportion parameter that characterizes the null space condition. In this paper, we propose new polynomial-time algorithms to establish upper bounds of the proportion parameter. We leverage on these techniques to find upper bounds and further develop a new procedure-tree search algorithm-that is able to precisely and quickly verify the null space condition. Numerical experiments show that the execution speed and accuracy of the results obtained from our methods far exceed those of the previous methods which rely on linear programming (LP) relaxation and semidefinite programming (SDP).

Perturbing the Hypothalamic-Pituitary-Adrenal Axis: A Mathematical Model for Interpreting PTSD Assessment Tests

We use a dynamical systems model of the hypothalamic-pituitary-adrenal (HPA) axis to understand the mechanisms underlying clinical protocols used to probe patient stress response. Specifically, we address dexamethasone (DEX) and ACTH challenge tests, which probe pituitary and adrenal gland responses, respectively. We show that some previously observed features and experimental responses can arise from a bistable mathematical model containing two steady-states, rather than relying on specific and permanent parameter changes due to physiological disruption. Moreover, we show that the timing of a perturbation relative to the intrinsic oscillation of the HPA axis can affect challenge test responses. Conventional mechanistic hypotheses supported and refuted by the challenge tests are reexamined by varying parameters in our mathematical model associated with these hypotheses. We show that (a) adrenal hyposensitivity can give rise to the responses seen in ACTH challenge tests and (b) enhanced cortisol-mediated suppression of the pituitary in subjects with PTSD is not necessary to explain the responses observed in DEX stress tests. We propose a new two-stage DEX/external stressor protocol to more clearly distinguish between the conventional hypothesis of enhanced suppression of the pituitary and bistable dynamics hypothesized in our model.

Exploiting the geometry of the solution space to reduce sensitivity to neuromotor noise

Throwing is a uniquely human skill that requires a high degree of coordination to successfully hit a target. Timing of ball release appears crucial as previous studies report required timing accuracies as short as 1-2ms, which however appear physiologically challenging. This study mathematically and experimentally demonstrates that humans can overcome these seemingly stringent timing requirements by shaping their hand trajectories to create extended timing windows, where ball releases achieve target hits despite temporal imprecision. Subjects practiced four task variations in a virtual environment, each with a distinct geometry of the solution space and different demands for timing. Model-based analyses of arm trajectories revealed that subjects first decreased timing error, followed by lengthening timing windows in their hand trajectories. This pattern was invariant across solution spaces, except for a control case. Hence, the exquisite skill that humans evolved for throwing is achieved by developing strategies that are less sensitive to temporal variability arising from neuromotor noise. This analysis also provides an explanation why coaches emphasize the "follow-through" in many ball sports.

Robust Exponential Memory in Hopfield Networks

The Hopfield recurrent neural network is a classical auto-associative model of memory, in which collections of symmetrically coupled McCulloch-Pitts binary neurons interact to perform emergent computation. Although previous researchers have explored the potential of this network to solve combinatorial optimization problems or store reoccurring activity patterns as attractors of its deterministic dynamics, a basic open problem is to design a family of Hopfield networks with a number of noise-tolerant memories that grows exponentially with neural population size. Here, we discover such networks by minimizing probability flow, a recently proposed objective for estimating parameters in discrete maximum entropy models. By descending the gradient of the convex probability flow, our networks adapt synaptic weights to achieve robust exponential storage, even when presented with vanishingly small numbers of training patterns. In addition to providing a new set of low-density error-correcting codes that achieve Shannon's noisy channel bound, these networks also efficiently solve a variant of the hidden clique problem in computer science, opening new avenues for real-world applications of computational models originating from biology.

Modeling the assembly order of multimeric heteroprotein complexes

Protein-protein interactions are the cornerstone of numerous biological processes. Although an increasing number of protein complex structures have been determined using experimental methods, relatively fewer studies have been performed to determine the assembly order of complexes. In addition to the insights into the molecular mechanisms of biological function provided by the structure of a complex, knowing the assembly order is important for understanding the process of complex formation. Assembly order is also practically useful for constructing subcomplexes as a step toward solving the entire complex experimentally, designing artificial protein complexes, and developing drugs that interrupt a critical step in the complex assembly. There are several experimental methods for determining the assembly order of complexes; however, these techniques are resource-intensive. Here, we present a computational method that predicts the assembly order of protein complexes by building the complex structure. The method, named Path-LzerD, uses a multimeric protein docking algorithm that assembles a protein complex structure from individual subunit structures and predicts assembly order by observing the simulated assembly process of the complex. Benchmarked on a dataset of complexes with experimental evidence of assembly order, Path-LZerD was successful in predicting the assembly pathway for the majority of the cases. Moreover, when compared with a simple approach that infers the assembly path from the buried surface area of subunits in the native complex, Path-LZerD has the strong advantage that it can be used for cases where the complex structure is not known. The path prediction accuracy decreased when starting from unbound monomers, particularly for larger complexes of five or more subunits, for which only a part of the assembly path was correctly identified. As the first method of its kind, Path-LZerD opens a new area of computational protein structure modeling and will be an indispensable approach for studying protein complexes.

Ebola could be eradicated through voluntary vaccination

Ebola virus disease (EVD) is a severe infection with an extremely high fatality rate spread through direct contact with body fluids. A promising Ebola vaccine (rVSV-ZEBOV) may soon become universally available. We constructed a game-theoretic model of Ebola incorporating individual decisions to vaccinate. We found that if a population adopts selfishly optimal vaccination strategies, then the population vaccination coverage falls negligibly short of the herd immunity level. We concluded that eradication of Ebola is feasible if voluntary vaccination programmes are coupled with focused public education efforts. We conducted uncertainty and sensitivity analysis to demonstrate that our findings do not depend on the choice of the epidemiological model parameters.

The role of backward cell migration in two-hit mutants' production in the stem cell niche

It has been discovered that there are two stem cell groups in the intestinal crypts: central stem cells (CeSCs), which are at the very bottom of the crypt, and border stem cells (BSCs), which are located between CeSCs and transit amplifying cells (TAs). Moreover, backward cell migration from BSCs to CeSCs has been observed. Recently, a bi-compartmental stochastic model, which includes CeSCs and BSCs, has been developed to investigate the probability of two-hit mutant production in the stem cell niche. In this project, we improve this stochastic model by adding the probability of backward cell migration to the model. The model suggests that the probability of two-hit mutant production increases when the frequency of backward cell migration increases. Furthermore, a small non-zero probability of backward cell migration leads to the largest range of optimal values for the frequency of symmetric divisions and the portion of divisions at each stem cell compartment in terms of delaying 2-hit mutant production. Moreover, the probability of two-hit mutant production is more sensitive to the probability of symmetric divisions than to the rate of backward cell migrations. The highest probability of two-hit mutant production corresponds to the case when all stem cell’s divisions are asymmetric.

A geometric method for eigenvalue problems with low-rank perturbations

We consider the problem of finding the spectrum of an operator taking the form of a low-rank (rank one or two) non-normal perturbation of a well-understood operator, motivated by a number of problems of applied interest which take this form. We use the fact that the system is a low-rank perturbation of a solved problem, together with a simple idea of classical differential geometry (the envelope of a family of curves) to completely analyse the spectrum. We use these techniques to analyse three problems of this form: a model of the oculomotor integrator due to Anastasio & Gad (2007 J. Comput. Neurosci.22, 239-254. (doi:10.1007/s10827-006-0010-x)), a continuum integrator model, and a non-local model of phase separation due to Rubinstein & Sternberg (1992 IMA J. Appl. Math.48, 249-264. (doi:10.1093/imamat/48.3.249)).

VirFinder: a novel k-mer based tool for identifying viral sequences from assembled metagenomic data

BACKGROUND

Identifying viral sequences in mixed metagenomes containing both viral and host contigs is a critical first step in analyzing the viral component of samples. Current tools for distinguishing prokaryotic virus and host contigs primarily use gene-based similarity approaches. Such approaches can significantly limit results especially for short contigs that have few predicted proteins or lack proteins with similarity to previously known viruses.

METHODS

We have developed VirFinder, the first k-mer frequency based, machine learning method for virus contig identification that entirely avoids gene-based similarity searches. VirFinder instead identifies viral sequences based on our empirical observation that viruses and hosts have discernibly different k-mer signatures. VirFinder's performance in correctly identifying viral sequences was tested by training its machine learning model on sequences from host and viral genomes sequenced before 1 January 2014 and evaluating on sequences obtained after 1 January 2014.

RESULTS

VirFinder had significantly better rates of identifying true viral contigs (true positive rates (TPRs)) than VirSorter, the current state-of-the-art gene-based virus classification tool, when evaluated with either contigs subsampled from complete genomes or assembled from a simulated human gut metagenome. For example, for contigs subsampled from complete genomes, VirFinder had 78-, 2.4-, and 1.8-fold higher TPRs than VirSorter for 1, 3, and 5 kb contigs, respectively, at the same false positive rates as VirSorter (0, 0.003, and 0.006, respectively), thus VirFinder works considerably better for small contigs than VirSorter. VirFinder furthermore identified several recently sequenced virus genomes (after 1 January 2014) that VirSorter did not and that have no nucleotide similarity to previously sequenced viruses, demonstrating VirFinder's potential advantage in identifying novel viral sequences. Application of VirFinder to a set of human gut metagenomes from healthy and liver cirrhosis patients reveals higher viral diversity in healthy individuals than cirrhosis patients. We also identified contig bins containing crAssphage-like contigs with higher abundance in healthy patients and a putative Veillonella genus prophage associated with cirrhosis patients.

CONCLUSIONS

This innovative k-mer based tool complements gene-based approaches and will significantly improve prokaryotic viral sequence identification, especially for metagenomic-based studies of viral ecology.

Modeling the dynamics of mouse iron body distribution: hepcidin is necessary but not sufficient

BACKGROUND

Iron is an essential element of most living organisms but is a dangerous substance when poorly liganded in solution. The hormone hepcidin regulates the export of iron from tissues to the plasma contributing to iron homeostasis and also restricting its availability to infectious agents. Disruption of iron regulation in mammals leads to disorders such as anemia and hemochromatosis, and contributes to the etiology of several other diseases such as cancer and neurodegenerative diseases. Here we test the hypothesis that hepcidin alone is able to regulate iron distribution in different dietary regimes in the mouse using a computational model of iron distribution calibrated with radioiron tracer data.

RESULTS

A model was developed and calibrated to the data from adequate iron diet, which was able to simulate the iron distribution under a low iron diet. However simulation of high iron diet shows considerable deviations from the experimental data. Namely the model predicts more iron in red blood cells and less iron in the liver than what was observed in experiments.

CONCLUSIONS

These results suggest that hepcidin alone is not sufficient to regulate iron homeostasis in high iron conditions and that other factors are important. The model was able to simulate anemia when hepcidin was increased but was unable to simulate hemochromatosis when hepcidin was suppressed, suggesting that in high iron conditions additional regulatory interactions are important.

Using electronic health records and Internet search information for accurate influenza forecasting

BACKGROUND

Accurate influenza activity forecasting helps public health officials prepare and allocate resources for unusual influenza activity. Traditional flu surveillance systems, such as the Centers for Disease Control and Prevention's (CDC) influenza-like illnesses reports, lag behind real-time by one to 2 weeks, whereas information contained in cloud-based electronic health records (EHR) and in Internet users' search activity is typically available in near real-time. We present a method that combines the information from these two data sources with historical flu activity to produce national flu forecasts for the United States up to 4 weeks ahead of the publication of CDC's flu reports.

METHODS

We extend a method originally designed to track flu using Google searches, named ARGO, to combine information from EHR and Internet searches with historical flu activities. Our regularized multivariate regression model dynamically selects the most appropriate variables for flu prediction every week. The model is assessed for the flu seasons within the time period 2013-2016 using multiple metrics including root mean squared error (RMSE).

RESULTS

Our method reduces the RMSE of the publicly available alternative (Healthmap flutrends) method by 33, 20, 17 and 21%, for the four time horizons: real-time, one, two, and 3 weeks ahead, respectively. Such accuracy improvements are statistically significant at the 5% level. Our real-time estimates correctly identified the peak timing and magnitude of the studied flu seasons.

CONCLUSIONS

Our method significantly reduces the prediction error when compared to historical publicly available Internet-based prediction systems, demonstrating that: (1) the method to combine data sources is as important as data quality; (2) effectively extracting information from a cloud-based EHR and Internet search activity leads to accurate forecast of flu.

A mathematical model coupling polarity signaling to cell adhesion explains diverse cell migration patterns

Protrusion and retraction of lamellipodia are common features of eukaryotic cell motility. As a cell migrates through its extracellular matrix (ECM), lamellipod growth increases cell-ECM contact area and enhances engagement of integrin receptors, locally amplifying ECM input to internal signaling cascades. In contrast, contraction of lamellipodia results in reduced integrin engagement that dampens the level of ECM-induced signaling. These changes in cell shape are both influenced by, and feed back onto ECM signaling. Motivated by experimental observations on melanoma cells lines (1205Lu and SBcl2) migrating on fibronectin (FN) coated topographic substrates (anisotropic post-density arrays), we probe this interplay between intracellular and ECM signaling. Experimentally, cells exhibited one of three lamellipodial dynamics: persistently polarized, random, or oscillatory, with competing lamellipodia oscillating out of phase (Park et al., 2017). Pharmacological treatments, changes in FN density, and substrate topography all affected the fraction of cells exhibiting these behaviours. We use these observations as constraints to test a sequence of hypotheses for how intracellular (GTPase) and ECM signaling jointly regulate lamellipodial dynamics. The models encoding these hypotheses are predicated on mutually antagonistic Rac-Rho signaling, Rac-mediated protrusion (via activation of Arp2/3 actin nucleation) and Rho-mediated contraction (via ROCK phosphorylation of myosin light chain), which are coupled to ECM signaling that is modulated by protrusion/contraction. By testing each model against experimental observations, we identify how the signaling layers interact to generate the diverse range of cell behaviors, and how various molecular perturbations and changes in ECM signaling modulate the fraction of cells exhibiting each. We identify several factors that play distinct but critical roles in generating the observed dynamic: (1) competition between lamellipodia for shared pools of Rac and Rho, (2) activation of RhoA by ECM signaling, and (3) feedback from lamellipodial growth or contraction to cell-ECM contact area and therefore to the ECM signaling level.

Forecasting influenza in Hong Kong with Google search queries and statistical model fusion

BACKGROUND

The objective of this study is to investigate predictive utility of online social media and web search queries, particularly, Google search data, to forecast new cases of influenza-like-illness (ILI) in general outpatient clinics (GOPC) in Hong Kong. To mitigate the impact of sensitivity to self-excitement (i.e., fickle media interest) and other artifacts of online social media data, in our approach we fuse multiple offline and online data sources.

METHODS

Four individual models: generalized linear model (GLM), least absolute shrinkage and selection operator (LASSO), autoregressive integrated moving average (ARIMA), and deep learning (DL) with Feedforward Neural Networks (FNN) are employed to forecast ILI-GOPC both one week and two weeks in advance. The covariates include Google search queries, meteorological data, and previously recorded offline ILI. To our knowledge, this is the first study that introduces deep learning methodology into surveillance of infectious diseases and investigates its predictive utility. Furthermore, to exploit the strength from each individual forecasting models, we use statistical model fusion, using Bayesian model averaging (BMA), which allows a systematic integration of multiple forecast scenarios. For each model, an adaptive approach is used to capture the recent relationship between ILI and covariates.

RESULTS

DL with FNN appears to deliver the most competitive predictive performance among the four considered individual models. Combing all four models in a comprehensive BMA framework allows to further improve such predictive evaluation metrics as root mean squared error (RMSE) and mean absolute predictive error (MAPE). Nevertheless, DL with FNN remains the preferred method for predicting locations of influenza peaks.

CONCLUSIONS

The proposed approach can be viewed a feasible alternative to forecast ILI in Hong Kong or other countries where ILI has no constant seasonal trend and influenza data resources are limited. The proposed methodology is easily tractable and computationally efficient.

Multi-scale computational study of the mechanical regulation of cell mitotic rounding in epithelia

Mitotic rounding during cell division is critical for preventing daughter cells from inheriting an abnormal number of chromosomes, a condition that occurs frequently in cancer cells. Cells must significantly expand their apical area and transition from a polygonal to circular apical shape to achieve robust mitotic rounding in epithelial tissues, which is where most cancers initiate. However, how cells mechanically regulate robust mitotic rounding within packed tissues is unknown. Here, we analyze mitotic rounding using a newly developed multi-scale subcellular element computational model that is calibrated using experimental data. Novel biologically relevant features of the model include separate representations of the sub-cellular components including the apical membrane and cytoplasm of the cell at the tissue scale level as well as detailed description of cell properties during mitotic rounding. Regression analysis of predictive model simulation results reveals the relative contributions of osmotic pressure, cell-cell adhesion and cortical stiffness to mitotic rounding. Mitotic area expansion is largely driven by regulation of cytoplasmic pressure. Surprisingly, mitotic shape roundness within physiological ranges is most sensitive to variation in cell-cell adhesivity and stiffness. An understanding of how perturbed mechanical properties impact mitotic rounding has important potential implications on, amongst others, how tumors progressively become more genetically unstable due to increased chromosomal aneuploidy and more aggressive.

Causal dissipation for the relativistic dynamics of ideal gases

We derive a general class of relativistic dissipation tensors by requiring that, combined with the relativistic Euler equations, they form a second-order system of partial differential equations which is symmetric hyperbolic in a second-order sense when written in the natural Godunov variables that make the Euler equations symmetric hyperbolic in the first-order sense. We show that this class contains a unique element representing a causal formulation of relativistic dissipative fluid dynamics which (i) is equivalent to the classical descriptions by Eckart and Landau to first order in the coefficients of viscosity and heat conduction and (ii) has its signal speeds bounded sharply by the speed of light. Based on these properties, we propose this system as a natural candidate for the relativistic counterpart of the classical Navier-Stokes equations.

Biophysical assay for tethered signaling reactions reveals tether-controlled activity for the phosphatase SHP-1

Tethered enzymatic reactions are ubiquitous in signaling networks but are poorly understood. A previously unreported mathematical analysis is established for tethered signaling reactions in surface plasmon resonance (SPR). Applying the method to the phosphatase SHP-1 interacting with a phosphorylated tether corresponding to an immune receptor cytoplasmic tail provides five biophysical/biochemical constants from a single SPR experiment: two binding rates, two catalytic rates, and a reach parameter. Tether binding increases the activity of SHP-1 by 900-fold through a binding-induced allosteric activation (20-fold) and a more significant increase in local substrate concentration (45-fold). The reach parameter indicates that this local substrate concentration is exquisitely sensitive to receptor clustering. We further show that truncation of the tether leads not only to a lower reach but also to lower binding and catalysis. This work establishes a new framework for studying tethered signaling processes and highlights the tether as a control parameter in clustered receptor signaling.

Stochastic Simulation Service: Bridging the Gap between the Computational Expert and the Biologist

We present StochSS: Stochastic Simulation as a Service, an integrated development environment for modeling and simulation of both deterministic and discrete stochastic biochemical systems in up to three dimensions. An easy to use graphical user interface enables researchers to quickly develop and simulate a biological model on a desktop or laptop, which can then be expanded to incorporate increasing levels of complexity. StochSS features state-of-the-art simulation engines. As the demand for computational power increases, StochSS can seamlessly scale computing resources in the cloud. In addition, StochSS can be deployed as a multi-user software environment where collaborators share computational resources and exchange models via a public model repository. We demonstrate the capabilities and ease of use of StochSS with an example of model development and simulation at increasing levels of complexity.

Scaling Properties of Dimensionality Reduction for Neural Populations and Network Models

Recent studies have applied dimensionality reduction methods to understand how the multi-dimensional structure of neural population activity gives rise to brain function. It is unclear, however, how the results obtained from dimensionality reduction generalize to recordings with larger numbers of neurons and trials or how these results relate to the underlying network structure. We address these questions by applying factor analysis to recordings in the visual cortex of non-human primates and to spiking network models that self-generate irregular activity through a balance of excitation and inhibition. We compared the scaling trends of two key outputs of dimensionality reduction—shared dimensionality and percent shared variance—with neuron and trial count. We found that the scaling properties of networks with non-clustered and clustered connectivity differed, and that the in vivo recordings were more consistent with the clustered network. Furthermore, recordings from tens of neurons were sufficient to identify the dominant modes of shared variability that generalize to larger portions of the network. These findings can help guide the interpretation of dimensionality reduction outputs in regimes of limited neuron and trial sampling and help relate these outputs to the underlying network structure.

We seek to understand how billions of neurons in the brain work together to give rise to everyday brain function. In most current experimental settings, we can only record from tens of neurons for a few hours at a time. A major question in systems neuroscience is whether our interpretation of how neurons interact would change if we monitor orders of magnitude more neurons and for substantially more time. In this study, we use realistic networks of model neurons, which allow us to analyze the activity from as many model neurons as we want for as long as we want. For these models, we found that we can identify the salient interactions among neurons and interpret their activity meaningfully within the range of neurons and recording time available in current experiments. Furthermore, we studied how the neural activity from the models reflects how the neurons are connected. These results help to guide the interpretation of analyses using populations of neurons in the context of the larger network to understand brain function.

A tale of two antennules: the performance of crab odour-capture organs in air and water

Odour capture is an important part of olfaction, where dissolved chemical cues (odours) are brought into contact with chemosensory structures. Antennule flicking by marine crabs is an example of discrete odour capture (sniffing) where an array of chemosensory hairs is waved through the water to create a flow-no flow pattern based on a narrow range of speeds, diameters of and spacings between hairs. Changing the speed of movement and spacing of hairs at this scale to manipulate flow represents a complicated fluid dynamics problem. In this study, we use numerical simulation of the advection and diffusion of a chemical gradient to reveal how morphological differences of the hair arrays affect odour capture. Specifically, we simulate odour capture by a marine crab (Callinectes sapidus) and a terrestrial crab (Coenobita rugosus) in both air and water to compare performance. We find that the antennule morphologies of each species are adaptions to capturing odours in their native habitats. Sniffing is an important part of odour capture for marine crabs in water where the diffusivity of odorant molecules is low and flow through the array is necessary. On the other hand, flow within the hair array diminishes odour-capture performance in air where diffusivities are high. This study highlights some of the adaptations necessary to transition from water to air.

Stochastic Network Models in Neuroscience: A Festschrift for Jack Cowan. Introduction to the Special Issue

Jack Cowan's remarkable career has spanned, and molded, the development of neuroscience as a quantitative and mathematical discipline combining deep theoretical contributions, rigorous mathematical work and groundbreaking biological insights. The Banff International Research Station hosted a workshop in his honor, on Stochastic Network Models of Neocortex, July 17-24, 2014. This accompanying Festschrift celebrates Cowan's contributions by assembling current research in stochastic phenomena in neural networks. It combines historical perspectives with new results including applications to epilepsy, path-integral methods, stochastic synchronization, higher-order correlation analysis, and pattern formation in visual cortex.

Reconfigurable origami-inspired acoustic waveguides

We combine numerical simulations and experiments to design a new class of reconfigurable waveguides based on three-dimensional origami-inspired metamaterials. Our strategy builds on the fact that the rigid plates and hinges forming these structures define networks of tubes that can be easily reconfigured. As such, they provide an ideal platform to actively control and redirect the propagation of sound. We design reconfigurable systems that, depending on the externally applied deformation, can act as networks of waveguides oriented along one, two, or three preferential directions. Moreover, we demonstrate that the capability of the structure to guide and radiate acoustic energy along predefined directions can be easily switched on and off, as the networks of tubes are reversibly formed and disrupted. The proposed designs expand the ability of existing acoustic metamaterials and exploit complex waveguiding to enhance control over propagation and radiation of acoustic energy, opening avenues for the design of a new class of tunable acoustic functional systems.

A mathematical model of GTPase pattern formation during single-cell wound repair

Rho GTPases are regulatory proteins whose patterns on the surface of a cell affect cell polarization, cell motility and repair of single-cell wounds. The stereotypical patterns formed by two such proteins, Rho and Cdc42, around laser-injured frog oocytes permit experimental analysis of GTPase activation, inactivation, segregation and crosstalk. Here, we review the development and analysis of a spatial model of GTPase dynamics that describe the formation of concentric zones of Rho and Cdc42 activity around wounds, and describe how this model has provided insights into the roles of the GTPase effector molecules protein kinase C (PKCβ and PKCη) and guanosine nucleotide dissociation inhibitor (GDI) in the wound response. We further demonstrate how the use of a 'sharp switch' model approximation in combination with bifurcation analysis can aid mapping the model behaviour in parameter space (approximate results confirmed with numerical simulation methods). Using these methods in combination with experimental manipulation of PKC activity (PKC overexpression (OE) and dominant negative conditions), we have shown that: (i) PKCβ most probably acts by enhancing existing positive feedbacks (from Rho to itself via the guanosine nucleotide exchange factor domain of Abr, and from Cdc42 to itself), (ii) PKCη most probably increases basal rates of inactivation (or possibly decreases basal rates of activation) of Rho and Cdc42, and (iii) the graded distribution of PKCη and its effect on initial Rho activity accounts for inversion of zones in a fraction (20%) of PKCη OE cells. Finally, we speculate that GDIs (which sequester GTPases) may have a critical role in defining the spatial domain, where the wound response may occur. This paper provides a more thorough exposition of the methods of analysis used in the investigation, whereas previous work on this topic was addressed to biologists and abbreviated such discussion.

Modeling and Simulation of Mucus Flow in Human Bronchial Epithelial Cell Cultures - Part I: Idealized Axisymmetric Swirling Flow

A multi-mode nonlinear constitutive model for mucus is constructed directly from micro- and macro-rheology experimental data on cell culture mucus, and a numerical algorithm is developed for the culture geometry and idealized cilia driving conditions. This study investigates the roles that mucus rheology, wall effects, and HBE culture geometry play in the development of flow profiles and the shape of the air-mucus interface. Simulations show that viscoelasticity captures normal stress generation in shear leading to a peak in the air-mucus interface at the middle of the culture and a depression at the walls. Linear and nonlinear viscoelastic regimes can be observed in cultures by varying the hurricane radius and mean rotational velocity. The advection-diffusion of a drug concentration dropped at the surface of the mucus flow is simulated as a function of Peclet number.

In the lungs, the airway surface liquid protects the airway epithelium from inhaled pathogens and particulates. It is well known that failure to properly clear mucus from the airways leads to chronic, even fatal, lung infections. To date, there is no validated constitutive model capable of recapitulating mucus rheology under diverse, physiological stress and deformation conditions. This gap has hindered studies into the causal relationship between mucus rheology and mucociliary clearance. Our modeling-experimental approach fulfills several purposes: to implement linear and nonlinear constitutive modeling of mucus from micro- and macro-rheology, to test constitutive modeling in an independent experimental system, to build a coarse-grained model of the PCL-mucus boundary condition, to measure and understand modifications in mucociliary transport during and after deposition of a controlled drug concentration, to measure and simulate both the flow and stress fields throughout the mucus layer, and to measure and simulate how the advection profiles in the culture couple with diffusion of particulates landing on the air-mucus interface. These results lay the groundwork for extension of the code to three dimensions and more realistic metachronal wave boundary conditions, both in cell cultures and in airways.

Do terrestrial hermit crabs sniff? Air flow and odorant capture by flicking antennules

Capture of odorant molecules by olfactory organs from the surrounding fluid is the first step of smelling. Sniffing intermittently moves fluid across sensory surfaces, increasing delivery rates of molecules to chemosensory receptors and providing discrete odour samples. Aquatic malacostracan crustaceans sniff by flicking olfactory antennules bearing arrays of chemosensory hairs (aesthetascs), capturing water in the arrays during downstroke and holding the sample during return stroke. Terrestrial malacostracans also flick antennules, but how their flicking affects odour capture from air is not understood. The terrestrial hermit crab, Coenobita rugosus, uses antennules bearing shingle-shaped aesthetascs to capture odours. We used particle image velocimetry to measure fine-scale fluid flow relative to a dynamically scaled physical model of a flicking antennule, and computational simulations to calculate diffusion to aesthetascs by odorant molecules carried in that flow. Air does not flow into the aesthetasc array during flick downstrokes or recovery strokes. Odorants are captured from air flowing around the outside of the array during flick downstrokes, when aesthetascs face upstream and molecule capture rates are 21% higher than for stationary antennules. Bursts of flicking followed by pauses deliver discrete odour samples to olfactory sensors, causing intermittency in odour capture by a different mechanism than aquatic crustaceans use.

Modelling stripe formation in zebrafish: an agent-based approach

Zebrafish have distinctive black stripes and yellow interstripes that form owing to the interaction of different pigment cells. We present a two-population agent-based model for the development and regeneration of these stripes and interstripes informed by recent experimental results. Our model describes stripe pattern formation, laser ablation and mutations. We find that fish growth shortens the necessary scale for long-range interactions and that iridophores, a third type of pigment cell, help align stripes and interstripes.