ASU Scholarship Showcase

Five immunocompetent C57BL/6-cBrd/cBrd/Cr (albino C57BL/6) mice were injected with GL261-luc2 cells, a cell line sharing characteristics of human glioblastoma multiforme (GBM). The mice were imaged using magnetic resonance (MR) at five separate time points to characterize growth and development of the tumor. After 25 days, the final tumor volumes of the mice varied from 12 mm³ to 62 mm³, even though mice were inoculated from the same tumor cell line under carefully controlled conditions. We generated hypotheses to explore large variances in final tumor size and tested them with our simple reaction-diffusion model in both a 3-dimensional (3D) finite difference method and a 2-dimensional (2D) level set method. The parameters obtained from a best-fit procedure, designed to yield simulated tumors as close as possible to the observed ones, vary by an order of magnitude between the three mice analyzed in detail. These differences may reflect morphological and biological variability in tumor growth, as well as errors in the mathematical model, perhaps from an oversimplification of the tumor dynamics or nonidentifiability of parameters. Our results generate parameters that match other experimental in vitro and in vivo measurements. Additionally, we calculate wave speed, which matches with other rat and human measurements.

Extreme events, a type of collective behavior in complex networked dynamical systems, often can have catastrophic consequences. To develop effective strategies to control extreme events is of fundamental importance and practical interest. Utilizing transportation dynamics on complex networks as a prototypical setting, we find that making the network “mobile” can effectively suppress extreme events. A striking, resonance-like phenomenon is uncovered, where an optimal degree of mobility exists for which the probability of extreme events is minimized. We derive an analytic theory to understand the mechanism of control at a detailed and quantitative level, and validate the theory numerically. Implications of our finding to current areas such as cybersecurity are discussed.

The relation between flux and fluctuation is fundamental to complex physical systems that support and transport flows. A recently obtained law predicts monotonous enhancement of fluctuation as the average flux is increased, which in principle is valid but only for large systems. For realistic complex systems of small sizes, this law breaks down when both the average flux and fluctuation become large. Here we demonstrate the failure of this law in small systems using real data and model complex networked systems, derive analytically a modified flux-fluctuation law, and validate it through computations of a large number of complex networked systems. Our law is more general in that its predictions agree with numerics and it reduces naturally to the previous law in the limit of large system size, leading to new insights into the flow dynamics in small-size complex systems with significant implications for the statistical and scaling behaviors of small systems, a topic of great recent interest.

Gompertz’s empirical equation remains the most popular one in describing cancer cell population growth in a wide spectrum of bio-medical situations due to its good fit to data and simplicity. Many efforts were documented in the literature aimed at understanding the mechanisms that may support Gompertz’s elegant model equation. One of the most convincing efforts was carried out by Gyllenberg and Webb. They divide the cancer cell population into the proliferative cells and the quiescent cells. In their two dimensional model, the dead cells are assumed to be removed from the tumor instantly. In this paper, we modify their model by keeping track of the dead cells remaining in the tumor. We perform mathematical and computational studies on this three dimensional model and compare the model dynamics to that of the model of Gyllenberg and Webb. Our mathematical findings suggest that if an avascular tumor grows according to our three-compartment model, then as the death rate of quiescent cells decreases to zero, the percentage of proliferative cells also approaches to zero. Moreover, a slow dying quiescent population will increase the size of the tumor. On the other hand, while the tumor size does not depend on the dead cell removal rate, its early and intermediate growth stages are very sensitive to it.

Evolutionary dynamical models for cyclic competitions of three species (e.g., rock, paper, and scissors, or RPS) provide a paradigm, at the microscopic level of individual interactions, to address many issues in coexistence and biodiversity. Real ecosystems often involve competitions among more than three species. By extending the RPS game model to five (rock-paper-scissors-lizard-Spock, or RPSLS) mobile species, we uncover a fundamental type of mesoscopic interactions among subgroups of species. In particular, competitions at the microscopic level lead to the emergence of various local groups in different regions of the space, each involving three species. It is the interactions among the groups that fundamentally determine how many species can coexist. In fact, as the mobility is increased from zero, two transitions can occur: one from a five- to a three-species coexistence state and another from the latter to a uniform, single-species state. We develop a mean-field theory to show that, in order to understand the first transition, group interactions at the mesoscopic scale must be taken into account. Our findings suggest, more broadly, the importance of mesoscopic interactions in coexistence of great many species.

Dynamical systems based on the minority game (MG) have been a paradigm for gaining significant insights into a variety of social and biological behaviors. Recently, a grouping phenomenon has been unveiled in MG systems of multiple resources (strategies) in which the strategies spontaneously break into an even number of groups, each exhibiting an identical oscillation pattern in the attendance of game players. Here we report our finding of spontaneous breakup of resources into three groups, each exhibiting period-three oscillations. An analysis is developed to understand the emergence of the striking phenomenon of triple grouping and period-three oscillations. In the presence of random disturbances, the triple-group/period-three state becomes transient, and we obtain explicit formula for the average transient lifetime using two methods of approximation. Our finding indicates that, period-three oscillation, regarded as one of the most fundamental behaviors in smooth nonlinear dynamical systems, can also occur in much more complex, evolutionary-game dynamical systems. Our result also provides a plausible insight for the occurrence of triple grouping observed, for example, in the U.S. housing market.

Understanding the dynamics of human movements is key to issues of significant current interest such as behavioral prediction, recommendation, and control of epidemic spreading. We collect and analyze big data sets of human movements in both cyberspace (through browsing of websites) and physical space (through mobile towers) and find a superlinear scaling relation between the mean frequency of visit〈f〉and its fluctuation σ : σ ∼〈f⟩^β with β ≈ 1.2. The probability distribution of the visiting frequency is found to be a stretched exponential function. We develop a model incorporating two essential ingredients, preferential return and exploration, and show that these are necessary for generating the scaling relation extracted from real data. A striking finding is that human movements in cyberspace and physical space are strongly correlated, indicating a distinctive behavioral identifying characteristic and implying that the behaviors in one space can be used to predict those in the other.

We sought to evaluate the reproducibility of a liquid chromatography-tandem mass spectrometry (LC-MS/MS)-based approach to measure the stable-isotope enrichment of in vivo-labeled muscle ATP synthase β subunit (β-F1-ATPase), a protein most directly involved in ATP production, and whose abundance is reduced under a variety of circumstances. Muscle was obtained from a rat infused with stable-isotope-labeled leucine. The muscle was homogenized, β-F1-ATPase immunoprecipitated, and the protein was resolved using 1D-SDS PAGE. Following trypsin digestion of the isolated protein, the resultant peptide mixtures were subjected to analysis by HPLC-ESI-MS/MS, which resulted in the detection of multiple β-F1-ATPase peptides. There were three β-F1-ATPase unique peptides with a leucine residue in the amino acid sequence, and which were detected with high intensity relative to other peptides and assigned with >95% probability to β-F1-ATPase. These peptides were specifically targeted for fragmentation to access their stable-isotope enrichment based on MS/MS peak areas calculated from extracted ion chromatographs for selected labeled and unlabeled fragment ions. Results showed best linearity (R[superscript 2] = 0.99) in the detection of MS/MS peak areas for both labeled and unlabeled fragment ions, over a wide range of amounts of injected protein, specifically for the β-F1-ATPase[subscript 134-143] peptide. Measured stable-isotope enrichment was highly reproducible for the β-F1-ATPase[subscript 134-143] peptide (CV = 2.9%). Further, using mixtures of synthetic labeled and unlabeled peptides we determined that there is an excellent linear relationship (R[superscript 2] = 0.99) between measured and predicted enrichment for percent enrichments ranging between 0.009% and 8.185% for the β-F1-ATPase[subscript 134-143] peptide. The described approach provides a reliable approach to measure the stable-isotope enrichment of in-vivo-labeled muscle β-F1-ATPase based on the determination of the enrichment of the β-F1-ATPase[subscript 134-143] peptide.

Background: Healthy individuals on the lower end of the insulin sensitivity spectrum also have a reduced gene expression response to exercise for specific genes. The goal of this study was to determine the relationship between insulin sensitivity and exercise-induced gene expression in an unbiased, global manner.

Methods and Findings: Euglycemic clamps were used to measure insulin sensitivity and muscle biopsies were done at rest and 30 minutes after a single acute exercise bout in 14 healthy participants. Changes in mRNA expression were assessed using microarrays, and miRNA analysis was performed in a subset of 6 of the participants using sequencing techniques. Following exercise, 215 mRNAs were changed at the probe level (Bonferroni-corrected P<0.00000115). Pathway and Gene Ontology analysis showed enrichment in MAP kinase signaling, transcriptional regulation and DNA binding. Changes in several transcription factor mRNAs were correlated with insulin sensitivity, including MYC, r=0.71; SNF1LK, r=0.69; and ATF3, r= 0.61 (5 corrected for false discovery rate). Enrichment in the 5’-UTRs of exercise-responsive genes suggested regulation by common transcription factors, especially EGR1. miRNA species of interest that changed after exercise included miR-378, which is located in an intron of the PPARGC1B gene.

Conclusions: These results indicate that transcription factor gene expression responses to exercise depend highly on insulin sensitivity in healthy people. The overall pattern suggests a coordinated cycle by which exercise and insulin sensitivity regulate gene expression in muscle.

Although insulin resistance in skeletal muscle is well-characterized, the role of circulating whole blood in the metabolic syndrome phenotype is not well understood. We set out to test the hypothesis that genes involved in inflammation, insulin signaling and mitochondrial function would be altered in expression in the whole blood of individuals with metabolic syndrome. We further wanted to examine whether similar relationships that we have found previously in skeletal muscle exist in peripheral whole blood cells. All subjects (n=184) were Latino descent from the Arizona Insulin Resistance registry. Subjects were classified based on the metabolic syndrome phenotype according to the National Cholesterol Education Program’s Adult Treatment Panel III. Of the 184 Latino subjects in the study, 74 were classified with the metabolic syndrome and 110 were without. Whole blood gene expression profiling was performed using the Agilent 4x44K Whole Human Genome Microarray. Whole blood microarray analysis identified 1,432 probes that were altered in expression ≥1.2 fold and P<0.05 after Benjamini-Hochberg in the metabolic syndrome subjects. KEGG pathway analysis revealed significant enrichment for pathways including ribosome, oxidative phosphorylation and MAPK signaling (all Benjamini-Hochberg P<0.05). Whole blood mRNA expression changes observed in the microarray data were confirmed by quantitative RT-PCR. Transcription factor binding motif enrichment analysis revealed E2F1, ELK1, NF-kappaB, STAT1 and STAT3 significantly enriched after Bonferroni correction (all P<0.05). The results of the present study demonstrate that whole blood is a useful tissue for studying the metabolic syndrome and its underlying insulin resistance although the relationship between blood and skeletal muscle differs.