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.

A structurally and compositionally well-defined and spectrally tunable artificial light-harvesting system has been constructed in which multiple organic dyes attached to a three-arm-DNA nanostructure serve as an antenna conjugated to a photosynthetic reaction center isolated from Rhodobacter sphaeroides 2.4.1. The light energy absorbed by the dye molecules is transferred to the reaction center, where charge separation takes place. The average number of DNA three-arm junctions per reaction center was tuned from 0.75 to 2.35. This DNA-templated multichromophore system serves as a modular light-harvesting antenna that is capable of being optimized for its spectral properties, energy transfer efficiency, and photostability, allowing one to adjust both the size and spectrum of the resulting structures. This may serve as a useful test bed for developing nanostructured photonic systems.

Time-resolved fluorescence spectroscopy was used to explore the pathway and kinetics of energy transfer in photosynthetic membrane vesicles (chromatophores) isolated from Rhodobacter (Rba.) sphaeroides cells harvested 2, 4, 6 or 24 hours after a transition from growth in high to low level illumination. As previously observed, this light intensity transition initiates the remodeling of the photosynthetic apparatus and an increase in the number of light harvesting 2 (LH2) complexes relative to light harvesting 1 (LH1) and reaction center (RC) complexes. It has generally been thought that the increase in LH2 complexes served the purpose of increasing the overall energy transmission to the RC. However, fluorescence lifetime measurements and analysis in terms of energy transfer within LH2 and between LH2 and LH1 indicate that, during the remodeling time period measured, only a portion of the additional LH2 generated are well connected to LH1 and the reaction center. The majority of the additional LH2 fluorescence decays with a lifetime comparable to that of free, unconnected LH2 complexes. The presence of large LH2-only domains has been observed by atomic force microscopy in Rba. sphaeroides chromatophores (Bahatyrova et al., Nature, 2004, 430, 1058), providing structural support for the existence of pools of partially connected LH2 complexes. These LH2-only domains represent the light-responsive antenna complement formed after a switch in growth conditions from high to low illumination, while the remaining LH2 complexes occupy membrane regions containing mixtures of LH2 and LH1–RC core complexes. The current study utilized a multi-parameter approach to explore the fluorescence spectroscopic properties related to the remodeling process, shedding light on the structure-function relationship of the photosynthetic assembles. Possible reasons for the accumulation of these largely disconnected LH2-only pools are discussed.

There are an increasing variety of applications in which peptides are both synthesized and used attached to solid surfaces. This has created a need for high throughput sequence analysis directly on surfaces. However, common sequencing approaches that can be adapted to surface bound peptides lack the throughput often needed in library-based applications. Here we describe a simple approach for sequence analysis directly on solid surfaces that is both high speed and high throughput, utilizing equipment available in most protein analysis facilities. In this approach, surface bound peptides, selectively labeled at their N-termini with a positive charge-bearing group, are subjected to controlled degradation in ammonia gas, resulting in a set of fragments differing by a single amino acid that remain spatially confined on the surface they were bound to. These fragments can then be analyzed by MALDI mass spectrometry, and the peptide sequences read directly from the resulting spectra.

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.

Chronic manganese (Mn) exposure is associated with neuromotor and neurocognitive deficits, but the exact mechanism of Mn neurotoxicity is still unclear. With the advent of magnetic resonance imaging (MRI), in-vivo analysis of brain structures has become possible. Among different sub-cortical structures, the basal ganglia (BG) has been investigated as a putative anatomical biomarker in MR-based studies of Mn toxicity. However, previous investigations have yielded inconsistent results in terms of regional MR signal intensity changes. These discrepancies may be due to the subtlety of brain alterations caused by Mn toxicity, coupled to analysis techniques that lack the requisite detection power. Here, based on brain MRI, we apply a 3D surface-based morphometry method on 3 bilateral basal ganglia structures in school-age children chronically exposed to Mn through drinking water to investigate the effect of Mn exposure on brain anatomy. Our method successfully pinpointed significant enlargement of many areas of the basal ganglia structures, preferentially affecting the putamen. Moreover, these areas showed significant correlations with fine motor performance, indicating a possible link between altered basal ganglia neurodevelopment and declined motor performance in high Mn exposed children.

Cells routinely compartmentalize enzymes for enhanced efficiency of their metabolic pathways. Here we report a general approach to construct DNA nanocaged enzymes for enhancing catalytic activity and stability. Nanocaged enzymes are realized by self-assembly into DNA nanocages with well-controlled stoichiometry and architecture that enabled a systematic study of the impact of both encapsulation and proximal polyanionic surfaces on a set of common metabolic enzymes. Activity assays at both bulk and single-molecule levels demonstrate increased substrate turnover numbers for DNA nanocage-encapsulated enzymes. Unexpectedly, we observe a significant inverse correlation between the size of a protein and its activity enhancement. This effect is consistent with a model wherein distal polyanionic surfaces of the nanocage enhance the stability of active enzyme conformations through the action of a strongly bound hydration layer. We further show that DNA nanocages protect encapsulated enzymes against proteases, demonstrating their practical utility in functional biomaterials and biotechnology.

The apolipoprotein E (APOE) e4 allele is the most prevalent genetic risk factor for Alzheimer's disease (AD). Hippocampal volumes are generally smaller in AD patients carrying the e4 allele compared to e4 noncarriers. Here we examined the effect of APOE e4 on hippocampal morphometry in a large imaging database—the Alzheimer's Disease Neuroimaging Initiative (ADNI). We automatically segmented and constructed hippocampal surfaces from the baseline MR images of 725 subjects with known APOE genotype information including 167 with AD, 354 with mild cognitive impairment (MCI), and 204 normal controls. High-order correspondences between hippocampal surfaces were enforced across subjects with a novel inverse consistent surface fluid registration method. Multivariate statistics consisting of multivariate tensor-based morphometry (mTBM) and radial distance were computed for surface deformation analysis. Using Hotelling's T2 test, we found significant morphological deformation in APOE e4 carriers relative to noncarriers in the entire cohort as well as in the nondemented (pooled MCI and control) subjects, affecting the left hippocampus more than the right, and this effect was more pronounced in e4 homozygotes than heterozygotes. Our findings are consistent with previous studies that showed e4 carriers exhibit accelerated hippocampal atrophy; we extend these findings to a novel measure of hippocampal morphometry. Hippocampal morphometry has significant potential as an imaging biomarker of early stage AD.

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.