Title: On Toric Ideals of some Statistical Models.
Abstract: We introduce hierarchical models from statistics and their associated Markov bases. These bases are often large and difficult to compute. We introduce certain toric ideals and their algebraic properties as an alternative way of thinking about these objects. One challenge is to describe hierarchical models with infinitely many generators in a finite way. Using a symmetric group action, we describe certain classes of models including progress made for the non-reducible Models. This is joint work with Uwe Nagel.
Title: Efficient Methods for Enforcing Contiguity in Geographic Districting Problems
Abstract: Every ten years, United States Congressional Districts must be redesigned in response to a national census. While the size of practical political districting problems is typically too large for exact optimization approaches, heuristics such as local search can help stakeholders quickly identify good (but suboptimal) plans that suit their objectives. However, enforcing a district contiguity constraint during local search can require significant computation; tools that can reduce contiguity-based computations in large practical districting problems are needed. This talk introduces the geo-graph framework for modeling geographic districting as a graph partitioning problem, discusses two geo-graph contiguity algorithms, and applies these algorithms to the creation of United States Congressional Districts from census blocks in several states. The experimental results demonstrate that the geo-graph contiguity assessment algorithms reduce the average number of edges visited during contiguity assessments by at least three orders of magnitude in every problem instance when compared with simple graph search, suggesting that the geo-graph model and its associated contiguity algorithms provide a powerful constraint assessment tool to political districting stakeholders. Joint work with Douglas M. King and Edward C. Sewell
Abstract: What can mathematics tell us about the treatment of cancer? In this talk I will present some of work that I have done in the modeling of tumor growth and treatment over the last fifteen years. Cancer is a myriad of individual diseases, with the common feature that an individual's own cells have become malignant. Thus, the treatment of cancer poses great challenges, since an attack must be mounted against cells that are nearly identical to normal cells. Mathematical models that describe tumor growth in tissue, the immune response, and the administration of different therapies can suggest treatment strategies that optimize treatment efficacy and minimize negative side-effects. However, the inherent complexity of the immune system and the spatial heterogeneity of human tissue gives rise to mathematical models that pose unique challenges for the mathematician. In this talk I will give a few examples of how mathematicians can work with clinicians and immunologists to understand the development of the disease and to design effective treatments. I will use mathematical tools from dynamical systems, optimal control and network analysis. This talk is intended for a general math audience: no knowledge of biology will be assumed.
Title: A Mathematical Model for the Force and Energetics in Competitive Running
Abstract: Competitive running has been around for thousands of years and many people have wondered what the optimal form and strategy is for running a race. In his paper, Behncke develops a simple mathematical model that focuses on the relationships and dynamics between the forces and energetics at play in order to find an optimal strategy for racing various distances. In this talk, I will describe the biomechanics, energetics, and optimization of running in Behncke's model and present his findings. Note: you do not have to like running to come to this talk :)
Title: Mathematical deep learning for drug discovery
Abstract: Designing efficient drugs for curing diseases is of essential importance for the 21stcentury's life science. Computer-aided drug design and discovery has obtained a significant recognition recently. However, the geometric complexity of protein-drug complexes remains a grand challenge to conventional computational methods, including machine learning algorithms. We assume that the physics of interest of protein-drug complexes lies on low-dimensional manifolds or subspaces embedded in a high-dimensional data space. We devise topological abstraction, differential geometry reduction, graph simplification, and multiscale modeling to construct low-dimensional representations of biomolecules in massive and diverse datasets. These representations are integrated with various deep learning algorithms for the predictions of protein-ligand binding affinity, drug toxicity, drug solubility, drug partition coefficient and mutation induced protein stability change, and for the discrimination of active ligands from decoys. I will briefly discuss the working principle of various techniques and their performance in D3R Grand Challenges,a worldwide competition series in computer-aided drug design and discovery (http://users.math.msu.edu/users/wei/D3R_GC3.pdf).
Title: Mathematics for Breast Cancer Research: investigating the role of iron.
Abstract: Breast cancer cells are addicted to iron. The mechanisms by which malignant cells acquire and contain high levels of iron are not completely understood. Furthermore, other cell types in a tumor, such as immune cells, can either aid or inhibit cancer cells from acquiring high levels of iron. In order to shed light in the question of how iron affects breast cancer growth, we are applying mathematical tools including polynomial dynamical systems over finite fields and 3D multiscale mathematical modeling. In this talk we will survey how mathematics is aiding in understanding the mechanisms of this addictive iron behavior of malignant cells, and present some preliminary work.
Title: Preconditioning for Accurate Solutions of the Biharmonic Eigenvalue Problem
Abstract: Solving ill-conditioned systems poses two basic problems: convergence and accuracy. Preconditioning can overcome slow convergence, but this is only practical if the preconditioned system can be formed sufficiently accurately. In fact, for a fourth order operator, existing eigenvalue algorithms may compute smaller eigenvalues with little or no accuracy in standard double precision. In this talk, we combine standard matrix eigenvalue solvers with an accurate preconditioning scheme in order to compute the smallest eigenvalue of the biharmonic operator to machine precision in spite of ill-conditioning. The results on various domains are compared with the best known computations from the literature to demonstrate the accuracy and applicability of the method.
Abstract: Understanding why some human populations remain extremely poor despite current development trends around the world remains a mystery to the natural, social and mathematical sciences. The poor rely on their immediate natural environment for subsistence and suffer from high burdens of infectious diseases. We present a general framework for modeling the ecology of poverty and disease, focusing on infectious diseases and renewable resources. Interactions between these ecological drivers of poverty and economics create reinforcing feedbacks resulting in three possible development regimes: 1) globally stable wealthy/healthy development, 2) globally stable unwealthy/unhealthy development, and 3) bistability. We show that the proportion of parameters leading to poverty is larger than that resulting in healthy/wealthy development; bistability consistently emerges as a general property of generalized disease-economic systems and that the systems under consideration are most sensitive to human disease parameters. The framework highlights feedbacks, processes and parameters that are important to measure in future studies of development, to identify effective and sustainable pathways out of poverty.
Title: Simulating Within-Vector Generation of the Malaria Parasite Diversity
Abstract:Plasmodium falciparum, the malaria parasite causing the most severe disease in humans, undergoes an asexual stage within the human host, and a sexual stage within the vector host, Anopheles mosquitoes. Because mosquitoes may be superinfected with parasites of different genotypes, this sexual stage of the parasite life-cycle presents the opportunity to create genetically novel parasites. To investigate the role that mosquitoes’ biology plays on the generation of parasite diversity, which introduces bottlenecks in the parasites’ development, we first constructed a stochastic model of parasite development within-mosquito, generating a distribution of parasite densities at five parasite life-cycle stages: gamete, zygote, ookinete, oocyst, and sporozoite, over the lifespan of a mosquito. We then coupled a model of sequence diversity generation via recombination between genotypes to the stochastic parasite population model. Our model framework shows that bottlenecks entering the oocyst stage decrease diversity from the initial gametocyte population in a mosquito’s blood meal, but diversity increases with the possibility for recombination and proliferation in the formation of sporozoites. Furthermore, when we begin with only two distinct parasite genotypes in the initial gametocyte population, the probability of transmitting more than two unique genotypes from mosquito to human is over 50% for a wide range of initial gametocyte densities.
Title: Modeling RNA secondary structure with auxiliary information
Abstract: The secondary structure of an RNA sequence plays an important role in determining its function, but directly observing RNA secondary structure is costly and difficult. Therefore, researchers have developed computational tools to predict the secondary structure of RNAs. One of the most popular methods is the Nearest Neighbor Thermodynamic Model (NNTM). More recently, high-throughput data that correlates with the state of a nucleotide being paired or unpaired has been developed. This data, called SHAPE for `selective 2'-hydroxyl acylation analyzed by primer extension', has been incorporated as auxiliary information into the objective function of NNTM with the goal of improving the accuracy of the predictions. This type of prediction is referred to as SHAPE-directed RNA secondary structure modeling. The addition of auxiliary information usually improves the accuracy of the predictions of NNTM. This talk will discuss challenges in RNA secondary structure modeling using NNTM and will provide ideas for developing synthetic auxiliary information that can be incorporated into NNTM to improve the accuracy of the predictions.
