BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Michael W. Deem

Michael W. Deem appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
Populations and Evolutioncond-mat.stat-mechTissues and OrgansBiomoleculescond-mat.softMolecular NetworksNeurons and Cognitionphysics.comp-phphysics.soc-phq-bioq-fin.GNSocial and Information Networks

Catalog footprint

What is connected

28works
12topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

28 published item(s)

preprint2016arXiv

Influenza Evolution and H3N2 Vaccine Effectiveness, with Application to the 2014/2015 Season

Influenza A is a serious disease that causes significant morbidity and mortality, and vaccines against the seasonal influenza disease are of variable effectiveness. In this paper, we discuss use of the $p_{\rm epitope}$ method to predict the dominant influenza strain and the expected vaccine effectiveness in the coming flu season. We illustrate how the effectiveness of the 2014/2015 A/Texas/50/2012 [clade 3C.1] vaccine against the A/California/02/2014 [clade 3C.3a] strain that emerged in the population can be estimated via pepitope. In addition, we show by a multidimensional scaling analysis of data collected through 2014, the emergence of a new A/New Mexico/11/2014-like cluster [clade 3C.2a] that is immunologically distinct from the A/California/02/2014-like strains.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Modular knowledge systems accelerate human migration in asymmetric random environments

Migration is a key mechanism for expansion of communities. In spatially heterogeneous environments, rapidly gaining knowledge about the local environment is key to the evolutionary success of a migrating population. For historical human migration, environmental heterogeneity was naturally asymmetric in the north-south (NS) and east-west (EW) directions. We here consider the human migration process in the Americas, modeled as random, asymmetric, modularly correlated environments. Knowledge about the environments determines the fitness of each individual. We present a phase diagram for asymmetry of migration as a function of carrying capacity and fitness threshold. We find that the speed of migration is proportional to the inverse complement of the spatial environmental gradient, and in particular we find that north-south migration rates are lower than east-west migration rates when the environmental gradient is higher in the north-south direction. Communication of knowledge between individuals can help to spread beneficial knowledge within the population. The speed of migration increases when communication transmits pieces of knowledge that contribute in a modular way to the fitness of individuals. The results for the dependence of migration rate on asymmetry and modularity are consistent with existing archaeological observations. The results for asymmetry of genetic divergence are consistent with patterns of human gene flow.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Prediction of Influenza B Vaccine Effectiveness from Sequence Data

Influenza is a contagious respiratory illness that causes significant human morbidity and mortality, affecting 5-15% of the population in a typical epidemic season. Human influenza epidemics are caused by types A and B, with roughly 25% of human cases due to influenza B. Influenza B is a single-stranded RNA virus with a high mutation rate, and both prior immune history and vaccination put significant pressure on the virus to evolve. Due to the high rate of viral evolution, the influenza B vaccine component of the annual influenza vaccine is updated, roughly every other year in recent years. To predict when an update to the vaccine is needed, an estimate of expected vaccine effectiveness against a range of viral strains is required. We here introduce a method to measure antigenic distance between the influenza B vaccine and circulating viral strains. The measure correlates well with effectiveness of the influenza B component of the annual vaccine in humans between 1979 and 2014. We discuss how this measure of antigenic distance may be used in the context of annual influenza vaccine design and prediction of vaccine effectiveness.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Development of modularity in the neural activity of children's brains

We study how modularity of the human brain changes as children develop into adults. Theory suggests that modularity can enhance the response function of a networked system subject to changing external stimuli. Thus, greater cognitive performance might be achieved for more modular neural activity, and modularity might likely increase as children develop. The value of modularity calculated from fMRI data is observed to increase during childhood development and peak in young adulthood. Head motion is deconvolved from the fMRI data, and it is shown that the dependence of modularity on age is independent of the magnitude of head motion. A model is presented to illustrate how modularity can provide greater cognitive performance at short times, i.e.\ task switching. A fitness function is extracted from the model. Quasispecies theory is used to predict how the average modularity evolves with age, illustrating the increase of modularity during development from children to adults that arises from selection for rapid cognitive function in young adults. Experiments exploring the effect of modularity on cognitive performance are suggested. Modularity may be a potential biomarker for injury, rehabilitation, or disease.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Hierarchy in Gene Expression is Predictive for Adult Acute Myeloid Leukemia

Cancer progresses with a change in the structure of the gene network in normal cells. We define a measure of organizational hierarchy in gene networks of affected cells in adult acute myeloid leukemia (AML) patients. With a retrospective cohort analysis based on the gene expression profiles of 116 acute myeloid leukemia patients, we find that the likelihood of future cancer relapse and the level of clinical risk are directly correlated with the level of organization in the cancer related gene network. We also explore the variation of the level of organization in the gene network with cancer progression. We find that this variation is non-monotonic, which implies the fitness landscape in the evolution of AML cancer cells is nontrivial. We further find that the hierarchy in gene expression at the time of diagnosis may be a useful biomarker in AML prognosis.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Modularity Enhances the Rate of Evolution in a Rugged Fitness Landscape

Biological systems are modular, and this modularity affects the evolution of biological systems over time and in different environments. We here develop a theory for the dynamics of evolution in a rugged, modular fitness landscape. We show analytically how horizontal gene transfer couples to the modularity in the system and leads to more rapid rates of evolution at short times. The model, in general, analytically demonstrates a selective pressure for the prevalence of modularity in biology. We use this model to show how the evolution of the influenza virus is affected by the modularity of the proteins that are recognized by the human immune system. Approximately 25\% of the observed rate of fitness increase of the virus could be ascribed to a modular viral landscape.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Quasispecies Theory for Evolution of Modularity

Biological systems are modular, and this modularity evolves over time and in different environments. A number of observations have been made of increased modularity in biological systems under increased environmental pressure. We here develop a quasispecies theory for the dynamics of modularity in populations of these systems. We show how the steady-state fitness in a randomly changing environment can be computed. We derive a fluctuation dissipation relation for the rate of change of modularity and use it to derive a relationship between rate of environmental changes and rate of growth of modularity. We also find a principle of least action for the evolved modularity at steady state. Finally, we compare our predictions to simulations of protein evolution and find them to be consistent.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Evolution: Life has Evolved to Evolve

Jim Shapiro synthesizes a great many observations about the mechanisms of evolution to reach the remarkable conclusion that large-scale modification, exchange, and rearrangement of the genome are common and should be viewed as fundamental features of life. In other words, the genome should be viewed not as mostly read-only with a few rare mutations, but rather as a fully-fledged read-write library of genetic functions under continuous revision. Revision of the genome occurs during cellular replication, during multicellular development, and during evolution of a population of individuals. DNA formatting controls the timing and location of genetic rearrangements, gene expression, and genetic repair. Each of these events is under the control of precise cellular circuits. Shapiro reviews the toolbox of natural genetic engineering that provides the functionalities necessary for efficient long-term genome restructuring.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Hierarchy of Gene Expression Data is Predictive of Future Breast Cancer Outcome

We calculate measures of hierarchy in gene and tissue networks of breast cancer patients. We find that the likelihood of metastasis in the future is correlated with increased values of network hierarchy for expression networks of cancer-associated genes, due to correlated expression of cancer-specific pathways. Conversely, future metastasis and quick relapse times are negatively correlated with values of network hierarchy in the expression network of all genes, due to dedifferentiation of gene pathways and circuits. These results suggest that hierarchy of gene expression may be useful as an additional biomarker for breast cancer prognosis.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

On the Universality of the Energy Response Function in the Long-Range Spin Glass Model with Sparse, Modular Couplings

We consider energy relaxation of the long-range spin glass model with sparse couplings, the so-called dilute Sherrington-Kirkpatrick (SK) model, starting from a random initial state. We consider the effect that modularity of the coupling matrix has on this relaxation dynamics. In the absence of finite size effects, the relaxation dynamics appears independent of modularity. For finite sizes, a more modular system reaches a less favorable energy at long times. For small sizes, a more modular system also has a less favorable energy at short times. For large sizes, modularity appears to lead to slightly more favorable energies at intermediate times. We discuss these results in the context of evolutionary theory, where horizontal gene transfer, absent in the Glauber equilibration dynamics of the SK model studied here, endows modular organisms with larger response functions at short times.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Physical Model of the Immune Response of Bacteria Against Bacteriophage Through the Adaptive CRISPR-Cas Immune System

Bacteria and archaea have evolved an adaptive, heritable immune system that recognizes and protects against viruses or plasmids. This system, known as the CRISPR-Cas system, allows the host to recognize and incorporate short foreign DNA or RNA sequences, called `spacers' into its CRISPR system. Spacers in the CRISPR system provide a record of the history of bacteria and phage coevolution. We use a physical model to study the dynamics of this coevolution as it evolves stochastically over time. We focus on the impact of mutation and recombination on bacteria and phage evolution and evasion. We discuss the effect of different spacer deletion mechanisms on the coevolutionary dynamics. We make predictions about bacteria and phage population growth, spacer diversity within the CRISPR locus, and spacer protection against the phage population.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Prediction of heart rate response to conclusion of spontaneous breathing trial by fluctuation dissipation theory

The non-equilibrium fluctuation dissipation theorem is applied to predict how critically ill patients respond to treatment, based upon data currently collected by standard hospital monitoring devices. This framework is demonstrated on a common procedure in critical care: the spontaneous breathing trial. It is shown that the responses of groups of similar patients to the spontaneous breathing trial can be predicted by the non-equilibrium fluctuation dissipation approach. This mathematical framework, when fully formed and applied to other clinical interventions, may serve as part of the basis for personalized critical care.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Evolutionary Processes in Finite Populations

We consider the evolution of large but finite populations on arbitrary fitness landscapes. We describe the evolutionary process by a Markov, Moran process. We show that to $\mathcal O(1/N)$, the time-averaged fitness is lower for the finite population than it is for the infinite population. We also show that fluctuations in the number of individuals for a given genotype can be proportional to a power of the inverse of the mutation rate. Finally, we show that the probability for the system to take a given path through the fitness landscape can be non-monotonic in system size.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A Multi-Scale Model for Correlation in B Cell VDJ Usage of Zebrafish

The zebrafish (\emph{Danio rerio}) is one of the model animals for study of immunology because the dynamics in the adaptive immune system of zebrafish are similar to that in higher animals. In this work, we built a multi-scale model to simulate the dynamics of B cells in the primary and secondary immune responses of zebrafish. We use this model to explain the reported correlation between VDJ usage of B cell repertoires in individual zebrafish. We use a delay ordinary differential equation (ODE) system to model the immune responses in the 6-month lifespan of a zebrafish. This mean field theory gives the number of high affinity B cells as a function of time during an infection. The sequences of those B cells are then taken from a distribution calculated by a "microscopic" random energy model. This generalized $NK$ model shows that mature B cells specific to one antigen largely possess a single VDJ recombination. The model allows first-principles calculation of the probability, $p$, that two zebrafish responding to the same antigen will select the same VDJ recombination. This probability $p$ increases with the B cell population size and the B cell selection intensity. The probability $p$ decreases with the B cell hypermutation rate. The multi-scale model predicts correlations in the immune system of the zebrafish that are highly similar to that from experiment.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A Novel Sequence-Based Antigenic Distance Measure for H1N1, with Application to Vaccine Effectiveness and the Selection of Vaccine Strains

H1N1 influenza causes substantial seasonal illness and was the subtype of the 2009 influenza pandemic. Precise measures of antigenic distance between the vaccine and circulating virus strains help researchers design influenza vaccines with high vaccine effectiveness. We here introduce a sequence-based method to predict vaccine effectiveness in humans. Historical epidemiological data show that this sequence-based method is as predictive of vaccine effectiveness as hemagglutination inhibition (HI) assay data from ferret animal model studies. Interestingly, the expected vaccine effectiveness is greater against H1N1 than H3N2, suggesting a stronger immune response against H1N1 than H3N2. The evolution rate of hemagglutinin in H1N1 is also shown to be greater than that in H3N2, presumably due to greater immune selection pressure.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Low-dimensional clustering detects incipient dominant influenza strain clusters

Influenza has been circulating in the human population and has caused three pandemics in the last century (1918 H1N1, 1957 H2N2, 1968 H3N2). The 2009 A(H1N1) was classified by the World Health Organization (WHO) as the fourth pandemic. Influenza has a high evolution rate, which makes vaccine design challenging. We here consider an approach for early detection of new dominant strains. By clustering the 2009 A(H1N1) sequence data, we found two main clusters. We then define a metric to detect the emergence of dominant strains. We show on historical H3N2 data that this method is able to identify a cluster around an incipient dominant strain before it becomes dominant. For example, for H3N2 as of March 30, 2009, the method detects the cluster for the new A/British Columbia/RV1222/2009 strain. This strain detection tool would appear to be useful for annual influenza vaccine selection.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Multicanonical molecular dynamics by variable-temperature thermostats and variable-pressure barostats

Sampling from flat energy or density distributions has proven useful in equilibrating complex systems with large energy barriers. Several thermostats and barostats are presented to sample these flat distributions by molecular dynamics. These methods use a variable temperature or pressure that is updated on the fly in the thermodynamic controller. These methods are illustrated on a Lennard-Jones system and a structure-based model of proteins.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Predicting Fixation Tendencies of the H3N2 Influenza Virus by Free Energy Calculation

Influenza virus evolves to escape from immune system antibodies that bind to it. We used free energy calculations with Einstein crystals as reference states to calculate the difference of antibody binding free energy ($ΔΔG$) induced by amino acid substitution at each position in epitope B of the H3N2 influenza hemagglutinin, the key target for antibody. A substitution with positive $ΔΔG$ value decreases the antibody binding constant. On average an uncharged to charged amino acid substitution generates the highest $ΔΔG$ values. Also on average, substitutions between small amino acids generate $ΔΔG$ values near to zero. The 21 sites in epitope B have varying expected free energy differences for a random substitution. Historical amino acid substitutions in epitope B for the A/Aichi/2/1968 strain of influenza A show that most fixed and temporarily circulating substitutions generate positive $ΔΔG$ values. We propose that the observed pattern of H3N2 virus evolution is affected by the free energy landscape, the mapping from the free energy landscape to virus fitness landscape, and random genetic drift of the virus. Monte Carlo simulations of virus evolution are presented to support this view.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Quantifying Selection and Diversity in Viruses by Entropy Methods, with Application to the Hemagglutinin of H3N2 Influenza

Many viruses evolve rapidly. For example, hemagglutinin of the H3N2 influenza A virus evolves to escape antibody binding. This evolution of the H3N2 virus means that people who have previously been exposed to an influenza strain may be infected by a newly emerged virus. In this paper, we use Shannon entropy and relative entropy to measure the diversity and selection pressure by antibody in each amino acid site of H3 hemagglutinin between the 1992-1993 season and the 2009-2010 season. Shannon entropy and relative entropy are two independent state variables that we use to characterize H3N2 evolution. The entropy method estimates future H3N2 evolution and migration using currently available H3 hemagglutinin sequences. First, we show that the rate of evolution increases with the virus diversity in the current season. The Shannon entropy of the sequence in the current season predicts relative entropy between sequences in the current season and those in the next season. Second, a global migration pattern of H3N2 is assembled by comparing the relative entropy flows of sequences sampled in China, Japan, the USA, and Europe. We verify this entropy method by describing two aspects of historical H3N2 evolution. First, we identify 54 amino acid sites in hemagglutinin that have evolved in the past to evade the immune system. Second, the entropy method shows that epitopes A and B in the top of hemagglutinin evolve most vigorously to escape antibody binding. Our work provides a novel entropy-based method to predict and quantify future H3N2 evolution and to describe the evolutionary history of H3N2.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The Emergence of Modularity in Biological Systems

In this review, we discuss modularity and hierarchy in biological systems. We review examples from protein structure, genetics, and biological networks of modular partitioning of the geometry of biological space. We review theories to explain modular organization of biology, with a focus on explaining how biology may spontaneously organize to a structured form. That is, we seek to explain how biology nucleated from among the many possibilities in chemistry. The emergence of modular organization of biological structure will be described as a symmetry-breaking phase transition, with modularity as the order parameter. Experimental support for this description will be reviewed. Examples will be presented from pathogen structure, metabolic networks, gene networks, and protein-protein interaction networks. Additional examples will be presented from ecological food networks, developmental pathways, physiology, and social networks.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Heterogeneous diversity of spacers within CRISPR

Clustered regularly interspaced short palindromic repeats (CRISPR) in bacterial and archaeal DNA have recently been shown to be a new type of anti-viral immune system in these organisms. We here study the diversity of spacers in CRISPR under selective pressure. We propose a population dynamics model that explains the biological observation that the leader-proximal end of CRISPR is more diversified and the leader-distal end of CRISPR is more conserved. This result is shown to be in agreement with recent experiments. Our results show thatthe CRISPR spacer structure is influenced by and provides a record of the viral challenges that bacteria face.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Quasispecies theory for finite populations

We present stochastic, finite-population formulations of the Crow-Kimura and Eigen models of quasispecies theory, for fitness functions that depend in an arbitrary way on the number of mutations from the wild type. We include back mutations in our description. We show that the fluctuation of the population numbers about the average values are exceedingly large in these physical models of evolution. We further show that horizontal gene transfer reduces by orders of magnitude the fluctuations in the population numbers and reduces the accumulation of deleterious mutations in the finite population due to Muller's ratchet. Indeed the population sizes needed to converge to the infinite population limit are often larger than those found in nature for smooth fitness functions in the absence of horizontal gene transfer. These analytical results are derived for the steady-state by means of a field-theoretic representation. Numerical results are presented that indicate horizontal gene transfer speeds up the dynamics of evolution as well.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Structure and Response in the World Trade Network

We examine how the structure of the world trade network has been shaped by globalization and recessions over the last 40 years. We show that by treating the world trade network as an evolving system, theory predicts the trade network is more sensitive to evolutionary shocks and recovers more slowly from them now than it did 40 years ago, due to structural changes in the world trade network induced by globalization. We also show that recession-induced change to the world trade network leads to an \emph{increased} hierarchical structure of the global trade network for a few years after the recession.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Theoretical Aspects of Immunity

The immune system recognizes a myriad of invading pathogens and their toxic products. It does so with a finite repertoire of antibodies and T cell receptors. We here describe theories that quantify the immune system dynamics. We describe how the immune system recognizes antigens by searching the large space of receptor molecules. We consider in some detail the theories that quantify the immune response to influenza and dengue fever. We review theoretical descriptions of the complementary evolution of pathogens that occurs in response to immune system pressure. Methods including bioinformatics, molecular simulation, random energy models, and quantum field theory contribute to a theoretical understanding of aspects of immunity.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Solution of the Crow-Kimura and Eigen models for alphabets of arbitrary size by Schwinger spin coherent states

To represent the evolution of nucleic acid and protein sequence, we express the parallel and Eigen models for molecular evolution in terms of a functional integral representation with an $h$-letter alphabet, lifting the two-state, purine/pyrimidine assumption often made in quasi-species theory. For arbitrary $h$ and a general mutation scheme, we obtain the solution of this model in terms of a maximum principle. Euler's theorem for homogeneous functions is used to derive this `thermodynamic' formulation of evolution. The general result for the parallel model reduces to known results for the purine/pyrimidine $h=2$ alphabet and the nucleic acid $h=4$ alphabet for the Kimura 3 ST mutation scheme. Examples are presented for the $h=4$ and $h=20$ cases. We derive the maximum principle for the Eigen model for general $h$. The general result for the Eigen model reduces to a known result for $h=2$. Examples are presented for the nucleic acid $h=4$ and the amino acid $h=20$ alphabet. An error catastrophe phase transition occurs in these models, and the order of the phase transition changes from second to first order for smooth fitness functions when the alphabet size is increased beyond two letters to the generic case. As examples, we analyze the general analytic solution for sharp peak, linear, quadratic, and quartic fitness functions.

0 citations0 reviewsNo code linkedOpen access
preprint1998arXiv

Numerical Observation of Disorder-Induced Anomalous Kinetics in the A + A -> 0 Reaction

We address via numerical simulation the two-dimensional bimolecular annihilation reaction $A + A \to \emptyset$ in the presence of quenched, random impurities. Renormalization group calculations have suggested that this reaction displays anomalous kinetics at long times, $c_{A}(t) \sim at^{δ-1}$, for certain types of topological or charged reactants and impurities. Both the exponent and the prefactor depend on the strength of disorder. The decay exponents determined from our simulations agree well with the values predicted by theory. The observed renormalization of the prefactor also agrees well with the values predicted by theory.

0 citations0 reviewsNo code linkedOpen access
preprint1997arXiv

A Configurational Bias Monte Carlo Method for Linear and Cyclic Peptides

In this manuscript, we describe a new configurational bias Monte Carlo technique for the simulation of peptides. We focus on the biologically relevant cases of linear and cyclic peptides. Our approach leads to an efficient, Boltzmann-weighted sampling of the torsional degrees of freedom in these biological molecules, a feat not possible with previous Monte Carlo and molecular dynamics methods.

0 citations0 reviewsNo code linkedOpen access