In a 2009 paper Yu, Etheridge and Cuthbertson proposed a model that was intended to address two questions. The first was a question related to Muller's ratchet, "What ratio of mutations must be beneficial for the meanfitness of a population to increase in time?" The second question was related to the Hill-Robertson effect, "If many beneficial mutations are introduced into a population, how much will competition slow the rate of adaptation?" They introduced a model of an asexually reproducing population of fixed size N and mutation rate mu and conjectured that the rate of adaptation is O(logN/(log logN)^2) for large N so long as there is some positive ratio of beneficial mutations. I will present an outline of my proof of this conjecture.

Understanding the folding of RNA sequences into three-dimensional structures is one of the fundamental challenges in molecular biology. For example, the branching of an RNA secondary structure is an important molecular characteristic yet difficult to predict correctly, especially for sequences on the scale of viral genomes. However, results from enumerative, probabilistic, analytic, and geometric combinatorics yield insights into RNA structure formation, and suggest new directions in viral capsid assembly.

Viscoelastic interactions of Myxococcus xanthus cells in a low-density domain close to the edge of a swarm have been recently studied in [1] using a combination of a cell-based three-dimensional Subcellular Element (SCE) model [1,2] and cell-tracking experiments. The model takes into account the flexible nature of M. xanthus as well as the effects of adhesion between cells arising from the interaction of the capsular polysaccharide covering two cells in contact with each other. New image and dynamic cell curvature analysis algorithms were used to track and measure the change in cell shapes that occur as flexible cells undergo significant bending during collisions resulting in direct calibration of the model parameters. It will be shown in this talk that flexibility of cells and the adhesive cell–cell and cell–substrate interactions of M. xanthus together with cell to aspect-ratio and directional reversals [3], play an important role in smooth cell gliding and more efficient swarming. In the second part of the talk results of the analysis of the three dimensional structures of fibrin networks, with and without cells, reconstructed from two-dimensional z-stacks of confocal microscopy sections using novel image analysis algorithms, will be presented. These images were used to establish microstructure-based models for studying the relationship between the structural features and the mechanical properties of the fibrin networks in blood clots. The change in the fibrin network alignment under applied strain and the elastic modulus values will be shown to agree well with the experimental data [4]. 1. C.W. Harvey, F. Morcos, C.R. Sweet, D. Kaiser, S. Chatterjee, X. Lu, D. Chen and M. Alber [2011], Study of elastic collisions of M. xanthus in swarms, Physical Biology 8, 026016. 2. C.R. Sweet, S. Chatterjee, Z. Xu, K. Bisordi, E.D. Rosen and M. Alber [2011], Modeling Platelet-Blood Flow Interaction Using Subcellular Element Langevin Method, J R Soc Interface, 2011 May 18. [Epub ahead of print], doi: 10.1098/rsif.2011.0180. 3. Y. Wu, Y. Jiang, D. Kaiser and M. Alber [2009], Periodic reversal of direction allows Myxobacteria to swarm, Proc. Natl. Acad. Sci. USA 106 4 1222-1227. 4. E. Kim, O.V. Kim, K.R. Machlus, X. Liu, T. Kupaev, J. Lioi, A.S. Wolberg, D.Z. Chen, E.D. Rosen, Z. Xu and M. Alber [2011], Soft Matter 7, 4983-4992.

A conference in honor of Harold Layton's 60th Birthday (Morning Session)

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.

Inferential summaries of tree estimates are useful in the setting of evolutionary biology, where phylogenetic trees have been built from DNA data since the 1960's. In bioinformatics, psychometrics and data mining, hierarchical clustering techniques output the same mathematical objects, and practitioners have similar questions about the stability and `generalizability' of these summaries. I will present applications of the Billera, Holmes, Vogtman (2001) distance to inferential problems both in the frequentist (bootstrap) and Bayesian contexts. I will compare the tree of trees representation to the Euclidean approximations of treespace made available through Multidimensional Scaling of the matrix of distances between trees. We also provide applications of the distances between trees to hierarchical clustering trees constructed from microarrays and phylogenetic trees of metagenomic data of bacteria in the gut. This talk contains joint work with John Chakerian and Alfred Spormann.

In this talk, I will discuss a new result in fluid flows through channels with permeable membranes with simple pumping dynamics. Fluid will be exchanged and metabolized in a simple reservoir and I will demonstrate the existence of optimal reservoir properties that may either maximize or minimized the amount of fluid being extracted across the channel walls. The biological relevance of this work may be seen by noting that all living organisms of a sufficient size rely on complex systems of tubular networks to efficiently collect, transport and distribute nutrients or waste. These networks exchange material with the interstitium via embedded channels leading to effective permeabilities across the wall separating the channel interior from the interstitium. In many invertebrates, for example, respiratory systems are made of complex tracheal systems that branch out through the entire body allowing for passive exchange of oxygen and carbon dioxide. In many of these systems, certain animals utilize various pumping mechanisms that alter the flow of the air or fluid being transported. Although the net effect of pumping of the averaged rates of fluid flow through the channel is typically well understood, it is still a largely open problem to understand how, and in what circumstances, pumping enables and enhances the exchange of material across channel walls. It has been demonstrated experimentally, for example, that when certain insects flap their wings, compression of the trachea allow for more efficient oxygen extraction, however it is unclear if this pumping is optimized for flight, oxygen uptake or neither, and understanding this problem quantitatively will shed insight on this biological process. Many of these interesting scenarios occur at low Reynolds number and this regime will be the focus of the presentation.

Biological oscillators such as the cell cycle, circadian clocks, and metabolic rhythms are ubiquitous across the domains of life. These biochemical oscillators co-exist in the same cells, often sharing and competing for resources. Are there mechanisms and regulatory principles that ensure harmony between these oscillators? Recent studies have shown that in addition to the transcriptional circadian clock, many organisms (including Arabidopsis) have a circadian redox rhythm driven by the organism's metabolic activities. It has been hypothesized that the redox rhythm is linked to the circadian clock, but the mechanism and the biological significance of this link have only begun to be investigated. In the first half of my talk, I will describe our work (in collaboration with the Dong lab at Duke) on the coupling of redox rhythms and the plant circadian clock. In the second half of my talk, I will discuss our work on the coupling of yeast metabolic cycle and the cell division cycle.

We discuss the question: What properties of solutions to coupled cell network systems are invariant under changes of coordinates that preserve network structure? This question was motivated by trying to understand the biological phenomenon of homeostasis in a mathematically satisfactory way. In its simplest mathematical form homeostasis can be described as follows. Given a stable equilibrium $x(lambda)$ of a system that depends on an input parameter $lambda$: When is some coordinate (say $x_j(lambda)$) approximately constant? First, we translate approximately constant to derivative approximately 0. This allows us to search for regions of homeostasis in a model using bifurcation theory like formulas. Second, we claim that there is a sense in which homeostasis can be thought of as a network invariant. This is joint work with Ian Stewart.

Vaccines protect their recipients by inducing long-term structural changes in populations of immune cells. Part of that restructuring is exactly analogous to Darwinian Selection. New antibody molecules are created by somatic mutation of existing antibody genes. Subsequently, the immune cell populations that possess these mutated receptors overtake the "wild-type" immune cells due  to the selective advantage they have acquired. Thus the immune system is vastly better prepared to recognize and eliminate the eliciting pathogen the next time around.

New sequencing and biosynthesis technologies, together with mathematical and computational tools, now allow us to investigate this fascinating and important phenomenon more deeply than ever before. I will illustrate this development with examples from the immune response to HIV infection.

Understanding and predicting cancer progression requires detailed interacting models of tumor and stromal cells, all calibrated to experimental data. Work to date has been limited by a lack of standardization of data representations of multicellular systems, though this is now being addressed through MultiCellDS (MultiCellular Data Standard) and digital cell lines, which are standardized representations of microenvironment-dependent cell phenotypes. Computational cancer modelers require biologically and mathematically consistent initialization routines to seed simulations with cells defined in digital cell lines. In this talk, we will briefly introduce a 3-D agent-based model designed for use in integrative computational biology. We introduce a “snapshot generator” that can take a digital cancer cell line and produce for the agent-based model an initial cell arrangement and a phenotypic state based upon analyses of the digital cell line data elements. We demonstrate 2-D monolayer and 3-D hanging drop simulations up to 500k MCF7 cells, a common breast cancer cell line. We additionally demonstrate the production of digital snapshots, standardized simulation output that will facilitate computational model comparison with a common core of analytical tools. With an early version of these tools, we assess the match between simulations and in vitro experiments. In the future, this work will be used to create and simulate combinations of tumor and stromal cells from appropriate digital cell lines in realistic tissue environments in order to understand, predict, and eventually control cancer progression in individual patients.

Mathematical modeling has become a widely-used tool for integrating
biological data, designing experiments, and ultimately understanding
biological systems. In recent years two important challenges for the
successful use of mathematical models have become apparent. One is that
models contain parameters that determine the behavior of the model, and
the values of these parameters are often hard to determine from the
available biological data. The other challenge is that many biological
systems exhibit a great deal of heterogeneity in behavior, so even if the
model parameters could be perfectly calibrated by pooling cell behaviors
to produce an “average cell model”, this model may not provide a good
description of any single cell in the population. In this seminar I will
describe a technique that we are using to integrate mathematical modeling
into experimental studies in a way that addresses both of these challenges.
We study endocrine pituitary cells that release a variety of hormones into
the blood, and our aim is to develop an approach for modeling the
behaviors of these cells with enough accuracy so that we can use the
models to make and test predictions in real time.

Abstract: The construction of flagellar motors in motile bacteria such as Salmonella is a carefully regulated genetic process. Among the structures that are built are the hook and the filament. The length of the hook is tightly controlled while the length of filaments is less so. However, if a filament is broken off it will regrow, while a broken hook will not regrow. The question that will be addressed in this talk is how Salmonella detects and regulates the length of these structures. This is related to the more general question of how physical properties (such as size or length) can be detected by chemical signals and what those mechanisms are. In this talk, I will present mathematical models for the regulation of hook and filament length. The model for hook length regulation is based on the hypothesis that the hook length is determined by the rate of secretion of the length regulatory molecule FliK and a cleavage reaction with the gatekeeper molecule FlhB. A stochastic model for this interaction is built and analyzed, showing excellent agreement with hook length data. The model for filament length regulation is based on the hypothesis that the growth of filaments is diffusion limited and is measured by negative feedback involving the regulatory protein FlgM. Thus, the model includes diffusion on a one-dimensional domain with a moving boundary, coupled with a negative feedback chemical network. The model shows excellent qualitative agreement with data, although there are some interesting unresolved issues related to the quantitative results.

I will describe three different types of problems inspired by the dynamics of the cytoskeleton: (i) the structural and kinetic aspects of dynamic instability in macromolecular assemblies, (ii) the collective flagella-like dynamics of ordered assemblies of active particles and cells, and (iii) the onset of oscillations and the persistence of strain in disordered aggregates of motors and filaments. In all cases, I will show how simple aspects of geometry, chemical kinetics and statistical and continuum mechanics allow us to explain biological observations in a minimal setting.

Cartilage physiology is regulated by a single population of specialized cells called chondrocytes. The chondroyctes are sparsely distributed within the extracellular matrix (ECM) and maintain a state of homeostasis in healthy tissue. ECM degeneration due to osteoarthritis can lead to compete degradation of cartilage surfaces, necessitating total joint replacement. Chondrocytes can be utilized to regenerate cartilage via tissue engineering approaches in which these cells are seeded in biocompatible and degradable biopolymer or hydrogel scaffolds. In such systems, biosynthetic activity of the cells in response to their non-native environment results in regeneration and accumulation of ECM constituents concurrent with degradation of the surrounding scaffold material. In this talk, mixture models are presented for interactions between biosynthesis of ECM constituents and ECM linking in cell-seeded scaffolds. Both ODE-based (temporal) models for evolution of average apparent densities and PDE-based (spatio-temporal) models will be presented for variables including unlinked ECM, linked ECM and scaffold. Model extensions accounting for cell proliferation will also be discussed. Of particular interest are model predictions for the evolution of solid phase apparent density, which is correlated with the compressive elastic modulus of the tissue construct. These models provide a quantitative framework for assessing and optimizing the design of engineered cell-scaffold systems and guiding strategies for articular cartilage tissue engineering.

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.

Literature:

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

In recent years it has become increasingly clear that stochasticity plays
an important role in many biological processes. Examples include bistable
genetic switches, noise enhanced robustness of oscillations, and fluctuation
enhanced sensitivity or "stochastic focusing." Numerous cellular systems
rely on spatial stochastic noise for robust performance. We examine the
need for stochastic models, report on the state of the art of algorithms and
software for modeling and simulation of stochastic biochemical systems, and
identify some computational challenges.

The behavior of inextensible fibers immersed in a fluid is of interest in a variety of applications ranging from polymer suspensions to actin filament transport. In these cases, the dynamics of an immersed fiber can play a large role in the observed macroscale fluid dynamics. The method of regularized Stokeslets provides a way to calculate fluid velocities in the Stokes fluid flow regime due to a collection of regularized point-forces without computing fluid velocities on an underlying grid. In this discussion, the method of Regularized Stokeslets will be used to model the dynamics of an inextensible flexible fiber immersed in a two-dimensional cellular background flow in comparison with results found in the experimental and mathematical literature. Studying this scenario with regularized Stokeslets provides insight into the documented stretch-coil transition and macroscale random walk behavior supported by mathematical models and experimental results.

Vaccines protect their recipients by inducing long-term structural changes in populations of immune cells. Part of that restructuring is exactly analogous to Darwinian Selection. New antibody molecules are created by somatic mutation of existing antibody genes. Subsequently, the immune cell populations that possess these mutated receptors overtake the "wild-type" immune cells due to the selective advantage they have acquired. Thus the immune system is vastly better prepared to recognize and eliminate the eliciting pathogen the next time around. New sequencing and biosynthesis technologies, together with mathematical and computational tools, now allow us to investigate this fascinating and important phenomenon more deeply than ever before. I will illustrate this development with examples from the immune response to HIV infection. Second lecture will focus on specifically mathematical questions.

The job of embryogenesis is to diversify cells into the hundreds of specialized cellular functions required by an animal or plant and to place those cells in the correct location in the embryo. To accomplish that complex job transcriptional regulation provides control circuits directing each dichotomy, each cell movement, and the patterning of cellular assemblies so that the animal that emerges from embryogenesis can feed and perform the necessary functions for survival of the species. The questions being addressed in this talk are to ask how gene regulatory networks are assembled, how do they change with time, and how do they accomplish the enormous regulatory task needed for building an embryo.

Biochemical reaction networks can be represented as graphs, with chemical complexes serving as vertices and reaction channels serving as directed edges. To each such network there is an associated stochastic process that models the behavior of the underlying system of interest. In this talk I will present results that relate the structure of the reaction network to the possible dynamics of the model. I will pay particular attention to how the results pertaining to the stochastic models compare and contrast with similar results pertaining to deterministic models.

Cells are highly organized and complex structures, with the capacity to segregate specific factors to spatially disjoint regions in a process called polarization. Polarization, which specifies a spatial axis in the cell, is a highly conserved biological process and is required for proper embryonic development, wound healing, and many other normal and pathological biological functions. Despite the importance of polarization, we do not fully understand how this protein segregation is initiated and maintained. In this talk, I will show how we can use numerical and analytical approaches to investigate how symmetry breaking begins the process of polarization, and how the geometry of the cell may play a role in the establishment and maintenance of spatial patterns associated with polarization.

Cancers are caused by mutations that may be inherited, induced by environmental factors, or result from DNA replication errors (R). We studied the relationship between the number of normal stem cell divisions and the risk of 17 cancer types in 69 countries throughout the world. The data revealed a strong correlation (median = 0.80) between cancer incidence and normal stem cell divisions in all countries, regardless of their environment. The major role of R mutations in cancer etiology was supported by an independent approach, based solely on cancer genome sequencing and epidemiological data, which suggested that R mutations are responsible for two-thirds of the mutations in human cancers. All of these results are consistent with epidemiological estimates of the fraction of cancers that can be prevented by changes in the environment. Moreover, they accentuate the importance of early detection and intervention to reduce deaths from the many cancers arising from unavoidable R mutations.