public 01:14:53

Stephan Huckemann : Statistical challenges in shape prediction of biomolecules

  -   Uploaded by schrett ( 45 Views )

The three-dimensional higher-order structure of biomolecules determines their functionality. While assessing primary structure is fairly easily accessible, reconstruction of higher order structure is costly. It often requires elaborate correction of atomic clashes, frequently not fully successful. Using RNA data, we describe a purely statistical method, learning error correction, drawing power from a two-scale approach. Our microscopic scale describes single suites by dihedral angles of individual atom bonds; here, addressing the challenge of torus principal component analysis (PCA) leads to a fundamentally new approach to PCA building on principal nested spheres by Jung et al. (2012). Based on an observed relationship with a mesoscopic scale, landmarks describing several suites, we use Fréchet means for angular shape and size-and-shape, correcting within-suite-backbone-to-backbone clashes. We validate this method by comparison to reconstructions obtained from simulations approximating biophysical chemistry and illustrate its power by the RNA example of SARS-CoV-2.

This is joint work with Benjamin Eltzner, Kanti V. Mardia and Henrik Wiechers.


Eltzner, B., Huckemann, S. F., Mardia, K. V. (2018): Torus principal component analysis with applications to RNA structure. Ann. Appl. Statist. 12(2), 1332?1359.

Jung, S., Dryden, I. L., Marron, J. S. (2012): Analysis of principal nested spheres. Biometrika, 99 (3), 551-568

Mardia, K. V., Wiechers, H., Eltzner, B., Huckemann, S. F. (2022). Principal component analysis and clustering on manifolds. Journal of Multivariate Analysis, 188, 104862, https://www.sciencedirect.com/science/article/pii/S0047259X21001408

Wiechers, H., Eltzner, B., Mardia, K. V., Huckemann, S. F. (2021). Learning torus PCA based classification for multiscale RNA backbone structure correction with application to SARS-CoV-2. To appear in the Journal of the Royal Statistical Society, Series C, bioRxiv https://doi.org/10.1101/2021.08.06.455406

public 01:34:56

Casey Diekman : Data Assimilation and Dynamical Systems Analysis of Circadian Rhythmicity and Entrainment

  -   Uploaded by schrett ( 50 Views )

Circadian rhythms are biological oscillations that align our physiology and behavior with the 24-hour environmental cycles conferred by the Earth’s rotation. In this talk, I will discuss two projects that focus on circadian clock cells in the brain and the entrainment of circadian rhythms to the light-dark cycle. Most of what we know about the electrical activity of circadian clock neurons comes from studies of nocturnal (night-active) rodents, hindering the translation of this knowledge to diurnal (day-active) humans. In the first part of the talk, we use data assimilation and patch-clamp recordings from the diurnal rodent Rhabdomys pumilio to build the first mathematical models of the electrophysiology of circadian neurons in a day-active species. We find that the electrical activity of circadian neurons is similar overall between nocturnal and diurnal rodents but that there are some interesting differences in their responses to inhibition. In the second part of the talk, we use tools from dynamical systems theory to study the reentrainment of a model of the human circadian pacemaker following perturbations that simulate jet lag. We show that the reentrainment dynamics are organized by invariant manifolds of fixed points of a 24-hour stroboscopic map and use these manifolds to explain a rapid reentrainment phenomenon that occurs under certain jet lag scenarios.

public 01:14:42

Rick Durrett : Overview of the semester

  -   Uploaded by root ( 82 Views )

public 01:14:46

David Basanta : Computational modeling of bone metastatic prostate cancer

  -   Uploaded by root ( 123 Views )

public 01:29:48

Jake Taylor-King : Generalized Jump Processes and Osteocyte Network Formation

  -   Uploaded by root ( 102 Views )

My talk will have two parts. PART I, From Birds to Bacteria: Generalised Velocity Jump Processes. There are various cases of animal movement where behaviour broadly switches between two modes of operation, corresponding to a long distance movement state and a resting or local movement state. In this talk, I will give a mathematical description of this process, adapted from Friedrich et. al. (2006). The approach allows the specification any running or waiting time distribution along with any angular and speed distributions. The resulting system of partial integro-differential equations are tumultuous and therefore it is necessary to both simplify and derive summary statistics. We derive an expression for the mean squared displacement, which shows good agreement with experimental data from the bacterium Escherichia coli and the gull Larus fuscus. Finally a large time diffusive approximation is considered via a Cattaneo approximation (Hillen, 2004). This leads to the novel result that the effective diffusion constant is dependent on the mean and variance of the running time distribution but only on the mean of the waiting time distribution. We also consider the Levy regime where the variance of the running distribution tends to infinity. This leads to a fractional diffusion equation for superdiffusive Levy walks and can be solved analytically. Our theory opens up new perspectives both for the systematic derivation of such equations, and for experimental data analysis of intermittent motion. I will also briefly discuss recent developments (by other researchers) within the field of velocity jump processes. PART II: Modelling Osteocyte Network Formation: Healthy and Cancerous Environments. Advanced prostate, breast, and lung cancer can metastasize to bone. In pathological bone, the highly regulated bone remodeling signaling pathway is disrupted. Within bone dendritic osteocytes form a spatial network allowing communication between osteocytes and the osteoblasts located on the bone surface. This communication network facilitates coordinated bone formation. In the presence of a cancerous microenvironment, the morphology of this network changes. Commonly osteocytes appear to be either overdifferentiated (i.e., there are more dendrites) or underdeveloped (i.e., dendrites do not fully form). In addition to structural changes, preliminary studies measuring the number of osteocytes per unit area using pathology slides show that the number density of osteocytes change from healthy to metastatic prostate and breast cancer xenografted mice. We present a stochastic agent-based model for bone formation incorporating osteoblasts and osteocytes that allows us to probe both network structure and number density of osteocytes in bone. Our model both allows for the simulation of our spatial network model and analysis of mean-field equations in the form of integro-partial differential equations. We consider variations of our model to test specific physiological hypotheses related to osteoblast differentiation; for example we can predict how changing measurable biological parameters, such as rates of bone secretion, rates of dendrite growth and rates of osteoblast differentiation can allow for qualitatively different network morphologies, and vice versa. We thenuse our model to hypothesize reasons for the limited efficacy of zoledronate therapy on metastatic breast cancer.

public 01:14:51

None : NO TALK

  -   Uploaded by root ( 84 Views )