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.

Background:
Data assimilation refers to methods for updating the state vector (initial condition) of a complex spatiotemporal model (such as a numerical weather model) by combining new observations with one or more prior forecasts. We consider the potential feasibility of this approach for making short-term (60-day) forecasts of the growth and spread of a malignant brain cancer (glioblastoma multiforme) in individual patient cases, where the observations are synthetic magnetic resonance images of a hypothetical tumor.

Results:
We apply a modern state estimation algorithm (the Local Ensemble Transform Kalman Filter), previously developed for numerical weather prediction, to two different mathematical models of glioblastoma, taking into account likely errors in model parameters and measurement uncertainties in magnetic resonance imaging. The filter can accurately shadow the growth of a representative synthetic tumor for 360 days (six 60-day forecast/update cycles) in the presence of a moderate degree of systematic model error and measurement noise.

Conclusions:
The mathematical methodology described here may prove useful for other modeling efforts in biology and oncology. An accurate forecast system for glioblastoma may prove useful in clinical settings for treatment planning and patient counseling.

The electronic structure of eight zinc-centered porphyrin macrocyclic molecules are investigated using density functional theory for ground-state properties, time-dependent density functional theory (TDDFT) for excited states, and Franck-Condon (FC) analysis for further characterization of the UV-vis spectrum. Symmetry breaking was utilized to find the lowest energy of the excited states for many states in the spectra. To confirm the theoretical modeling, the spectroscopic result from zinc phthalocyanine (ZnPc) is used to compare to the TDDFT and FC result. After confirmation of the modeling, five more planar molecules are investigated: zinc tetrabenzoporphyrin (ZnTBP), zinc tetrabenzomonoazaporphyrin (ZnTBMAP), zinc tetrabenzocisdiazaporphyrin (ZnTBcisDAP), zinc tetrabenzotransdiazaporphyrin (ZnTBtransDAP), and zinc tetrabenzotriazaporphyrin (ZnTBTrAP). The two latter molecules are then compared to their phenylated sister molecules: zinc monophenyltetrabenzotriazaporphyrin (ZnMPTBTrAP) and zinc diphenyltetrabenzotransdiazaporphyrin (ZnDPTBtransDAP). The spectroscopic results from the synthesis of ZnMPTBTrAP and ZnDPTBtransDAP are then compared to their theoretical models and non-phenylated pairs. While the Franck-Condon results were not as illuminating for every B-band, the Q-band results were successful in all eight molecules, with a considerable amount of spectral analysis in the range of interest between 300 and 750 nm. The π-π^* transitions are evident in the results for all of the Q bands, while satellite vibrations are also visible in the spectra. In particular, this investigation finds that, while ZnPc has a D_4h symmetry at ground state, a C_4v symmetry is predicted in the excited-state Q band region. The theoretical results for ZnPc found an excitation energy at the Q-band 0-0 transition of 1.88 eV in vacuum, which is in remarkable agreement with published gas-phase spectroscopy, as well as our own results of ZnPc in solution with Tetrahydrofuran that are provided in this paper.

We report a study of the magnetic domain structure of nanocrystalline thin films of nickel-zinc ferrite. The ferrite films were synthesized using aqueous spin-spray coating at low temperature (∼90 °C) and showed high complex permeability in the GHz range. Electron microscopy and microanalysis revealed that the films consisted of columnar grains with uniform chemical composition. Off-axis electron holography combined with magnetic force microscopy indicated a multi-grain domain structure with in-plane magnetization. The correlation between the magnetic domain morphology and crystal structure is briefly discussed.

We develop a framework to uncover and analyse dynamical anomalies from massive, nonlinear and non-stationary time series data. The framework consists of three steps: preprocessing of massive datasets to eliminate erroneous data segments, application of the empirical mode decomposition and Hilbert transform paradigm to obtain the fundamental components embedded in the time series at distinct time scales, and statistical/scaling analysis of the components. As a case study, we apply our framework to detecting and characterizing high-frequency oscillations (HFOs) from a big database of rat electroencephalogram recordings. We find a striking phenomenon: HFOs exhibit on–off intermittency that can be quantified by algebraic scaling laws. Our framework can be generalized to big data-related problems in other fields such as large-scale sensor data and seismic data analysis.

The ongoing Zika virus (ZIKV) epidemic in the Americas poses a major global public health emergency. While ZIKV is transmitted from human to human by bites of Aedes mosquitoes, recent evidence indicates that ZIKV can also be transmitted via sexual contact with cases of sexually transmitted ZIKV reported in Argentina, Canada, Chile, France, Italy, New Zealand, Peru, Portugal, and the USA. Yet, the role of sexual transmission on the spread and control of ZIKV infection is not well-understood. We introduce a mathematical model to investigate the impact of mosquito-borne and sexual transmission on the spread and control of ZIKV and calibrate the model to ZIKV epidemic data from Brazil, Colombia, and El Salvador. Parameter estimates yielded a basic reproduction number R_0 = 2.055 (95% CI: 0.523–6.300), in which the percentage contribution of sexual transmission is 3.044% (95% CI: 0.123–45.73). Our sensitivity analyses indicate that R_0 is most sensitive to the biting rate and mortality rate of mosquitoes while sexual transmission increases the risk of infection and epidemic size and prolongs the outbreak. Prevention and control efforts against ZIKV should target both the mosquito-borne and sexual transmission routes.

Successful identification of directed dynamical influence in complex systems is relevant to significant problems of current interest. Traditional methods based on Granger causality and transfer entropy have issues such as difficulty with nonlinearity and large data requirement. Recently a framework based on nonlinear dynamical analysis was proposed to overcome these difficulties. We find, surprisingly, that noise can counterintuitively enhance the detectability of directed dynamical influence. In fact, intentionally injecting a proper amount of asymmetric noise into the available time series has the unexpected benefit of dramatically increasing confidence in ascertaining the directed dynamical influence in the underlying system. This result is established based on both real data and model time series from nonlinear ecosystems. We develop a physical understanding of the beneficial role of noise in enhancing detection of directed dynamical influence.

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.

Background: Androgens bind to the androgen receptor (AR) in prostate cells and are essential survival factors for healthy prostate epithelium. Most untreated prostate cancers retain some dependence upon the AR and respond, at least transiently, to androgen ablation therapy. However, the relationship between endogenous androgen levels and cancer etiology is unclear. High levels of androgens have traditionally been viewed as driving abnormal proliferation leading to cancer, but it has also been suggested that low levels of androgen could induce selective pressure for abnormal cells. We formulate a mathematical model of androgen regulated prostate growth to study the effects of abnormal androgen levels on selection for pre-malignant phenotypes in early prostate cancer development.

Results: We find that cell turnover rate increases with decreasing androgen levels, which may increase the rate of mutation and malignant evolution. We model the evolution of a heterogeneous prostate cell population using a continuous state-transition model. Using this model we study selection for AR expression under different androgen levels and find that low androgen environments, caused either by low serum testosterone or by reduced 5α-reductase activity, select more strongly for elevated AR expression than do normal environments. High androgen actually slightly reduces selective pressure for AR upregulation. Moreover, our results suggest that an aberrant androgen environment may delay progression to a malignant phenotype, but result in a more dangerous cancer should one arise.

Conclusions: The model represents a useful initial framework for understanding the role of androgens in prostate cancer etiology, and it suggests that low androgen levels can increase selection for phenotypes resistant to hormonal therapy that may also be more aggressive. Moreover, clinical treatment with 5α-reductase inhibitors such as finasteride may increase the incidence of therapy resistant cancers.

The development of non-volatile logic through direct coupling of spontaneous ferroelectric polarization with semiconductor charge carriers is nontrivial, with many issues, including epitaxial ferroelectric growth, demonstration of ferroelectric switching and measurable semiconductor modulation. Here we report a true ferroelectric field effect—carrier density modulation in an underlying Ge(001) substrate by switching of the ferroelectric polarization in epitaxial c-axis-oriented BaTiO₃ grown by molecular beam epitaxy. Using the density functional theory, we demonstrate that switching of BaTiO₃ polarization results in a large electric potential change in Ge. Aberration-corrected electron microscopy confirms BaTiO₃ tetragonality and the absence of any low-permittivity interlayer at the interface with Ge. The non-volatile, switchable nature of the single-domain out-of-plane ferroelectric polarization of BaTiO₃ is confirmed using piezoelectric force microscopy. The effect of the polarization switching on the conductivity of the underlying Ge is measured using microwave impedance microscopy, clearly demonstrating a ferroelectric field effect.