Predictive oncology: a review of multidisciplinary, multiscale in silico modeling linking phenotype, morphology and growth.

Empirical evidence and theoretical studies suggest that the phenotype, i.e., cellular- and molecular-scale dynamics, including proliferation rate and adhesiveness due to microenvironmental factors and gene expression that govern tumor growth and invasiveness, also determine gross tumor-scale morphology. It has been difficult to quantify the relative effect of these links on disease progression and prognosis using conventional clinical and experimental methods and observables. As a result, successful individualized treatment of highly malignant and invasive cancers, such as glioblastoma, via surgical resection and chemotherapy cannot be offered and outcomes are generally poor. What is needed is a deterministic, quantifiable method to enable understanding of the connections between phenotype and tumor morphology. Here, we critically assess advantages and disadvantages of recent computational modeling efforts (e.g., continuum, discrete, and cellular automata models) that have pursued this understanding. Based on this assessment, we review a multiscale, i.e., from the molecular to the gross tumor scale, mathematical and computational "first-principle" approach based on mass conservation and other physical laws, such as employed in reaction-diffusion systems. Model variables describe known characteristics of tumor behavior, and parameters and functional relationships across scales are informed from in vitro, in vivo and ex vivo biology. We review the feasibility of this methodology that, once coupled to tumor imaging and tumor biopsy or cell culture data, should enable prediction of tumor growth and therapy outcome through quantification of the relation between the underlying dynamics and morphological characteristics. In particular, morphologic stability analysis of this mathematical model reveals that tumor cell patterning at the tumor-host interface is regulated by cell proliferation, adhesion and other phenotypic characteristics: histopathology information of tumor boundary can be inputted to the mathematical model and used as a phenotype-diagnostic tool to predict collective and individual tumor cell invasion of surrounding tissue. This approach further provides a means to deterministically test effects of novel and hypothetical therapy strategies on tumor behavior.

A LAPLACE TRANSFORM PAIR MODEL TO DETERMINE BREMSSTRAHLUNG SPECTRA FROM ATTENUATION DATA (RADIATION, X-RAY, SPECTRAL MODELING)

Every x-ray attenuation curve inherently contains all the information necessary to extract the complete energy spectrum of a beam. To date, attempts to obtain accurate spectral information from attenuation data have been inadequate.^ This investigation presents a mathematical pair model, grounded in physical reality by the Laplace Transformation, to describe the attenuation of a photon beam and the corresponding bremsstrahlung spectral distribution. In addition the Laplace model has been mathematically extended to include characteristic radiation in a physically meaningful way. A method to determine the fraction of characteristic radiation in any diagnostic x-ray beam was introduced for use with the extended model.^ This work has examined the reconstructive capability of the Laplace pair model for a photon beam range of from 50 kVp to 25 MV, using both theoretical and experimental methods.^ In the diagnostic region, excellent agreement between a wide variety of experimental spectra and those reconstructed with the Laplace model was obtained when the atomic composition of the attenuators was accurately known. The model successfully reproduced a 2 MV spectrum but demonstrated difficulty in accurately reconstructing orthovoltage and 6 MV spectra. The 25 MV spectrum was successfully reconstructed although poor agreement with the spectrum obtained by Levy was found.^ The analysis of errors, performed with diagnostic energy data, demonstrated the relative insensitivity of the model to typical experimental errors and confirmed that the model can be successfully used to theoretically derive accurate spectral information from experimental attenuation data. ^

Elevated albumin in retinas of monkeys with experimental glaucoma.

PURPOSE: To establish the identity of a prominent protein, approximately 70 kDa, that is markedly increased in the retina of monkeys with experimental glaucoma compared with the fellow control retina, the relationship to glaucoma severity, and its localization in the retina. METHODS: Retinal extracts were subjected to 2-D gel electrophoresis to identify differentially expressed proteins. Purified peptides from the abundant 70 kDa protein were analyzed and identified by liquid chromatography/mass spectrometry/mass spectrometry (LC/MS/MS) separation, and collision-induced dissociation sequencing. Protein identity was performed on MASCOT (Matrix Science, Boston, MA) and confirmed by Western blot. The relationship between the increase in this protein and glaucoma severity was investigated by regression analyses. Protein localization in retina was evaluated by immunohistochemistry with confocal imaging. RESULTS: The abundant protein was identified as Macaca mulatta serum albumin precursor (67 kDa) from eight non-overlapping proteolytic fragments, and the identity was confirmed by Western blot. The average increase in retinal albumin content was 2.3 fold (P = 0.015). In glaucoma eyes, albumin was localized to some neurons of the inner nuclear layer, in the inner plexiform layer, and along the vitreal surface, but it was only found in blood vessels in control retinas. CONCLUSIONS: Albumin is the abundant protein found in the glaucomatous monkey retinas. The increased albumin is primarily localized to the inner retina where oxidative damage associated with experimental glaucoma is known to be prominent. Since albumin is a major antioxidant, the increase of albumin in the retinas of eyes with experimental glaucoma may serve to protect the retina against oxidative damage.

{\it In vivo\/} magnetic resonance studies of experimental liver disease: Carbon tetrachloride hepatotoxicity and alcohol-induced fatty liver in rat

Magnetic resonance imaging (MRI) and magnetic resonance spectroscopy (MRS) were used to non-invasively determine if cirrhosis induced by carbon tetrachloride (CCl$\sb4$) and phospholipase-D (PLD) could be distinguished from fatty infiltration in rat. MRS localization and water suppression methods were developed, implemented and evaluated in terms of their application to in vivo proton NMR studies of experimental liver disease. MRS studies were also performed to quantitate fatty infiltration resulting from carbon tetrachloride (CCl$\sb4$) or alcohol (ethanol) administration and the MRS results were confirmed using biochemical total lipid analysis and histology. $\rm T\sb1$ weighted MR images acquired weekly, 48 hours post administration, demonstrated only a slight increase in overall liver intensity with CCl$\sb4$ or alcohol administration, which is consistent with previously reported results. The MR images were able to detect nodules resulting from CCl$\sb4$+PLD induced cirrhosis as hypointense regions, also consistent with previous reports. Localized in vivo water and lipid proton $\rm T\sb1$ relaxation time measurements were performed and demonstrated no statistically significant trends for either agent. In vivo proton spectra were also acquired using stimulated echo techniques to quantitatively follow the changes in liver lipid content. The changes in liver lipid content observed using MRS were verified by total lipid analysis using the Folch technique and histology. The in vivo $\rm T\sb1$ and lipid quantification data str inconsistent with the previous hypothesis that the changes in $\rm T\sb1$ weighted images were the result of increased "free" water content and, therefore, increased water $\rm T\sb1$ relaxation times. These data indicate that the long term changes are more likely the result of changes in lipid content. The data are also shown to agree with the accepted hypothesis that the time course and mechanism of fatty infiltration are different for CCl$\sb4$ and alcohol. The hypothesis that the lipids resulting from either protocol are from the same lipid fraction(s), presumably triglycerides, is also supported. And lastly, on the basis of MR images and quantitative MRS lipid information, it was shown that cirrhosis could be distinguished from fatty infiltration. ^

The dynamics of calmodulin signaling revealed using optical methods

The Ca2+-binding protein calmodulin (CaM) is a key transducer of Ca2+ oscillations by virtue of its ability to bind Ca 2+ selectively and then interact specifically with a large number of downstream enzymes and proteins. It remains unclear whether Ca2+ -dependent signaling alone can activate the full range of Ca 2+/CaM regulated processes or whether other regulatory schemes in the cell exist that allow specific targeting of CaM to subsets of Ca 2+/CaM binding sites or regions of the cell. Here we investigate the possibility that alterations of the availability of CaM may serve as a potential cellular mechanism for regulating the activation of CaM-dependent targets. By utilizing sensitive optical techniques with high spatial and temporal resolution, we examine the intracellular dynamics of CaM signaling at a resolution previously unattainable. After optimizing and characterizing both the optical methods and fluorescently labeled probes for intracellular measurements, the diffusion of CaM in the cytoplasm of HEK293 cells was analyzed. It was discovered that the diffusion characteristics of CaM are similar to that of a comparably sized inert molecule. Independent manipulation of experimental parameters, including increases in total concentrations of CaM and intracellular Ca2+ levels, did not change the diffusion of CaM in the cytoplasm. However, changes in diffusion were seen when the concentration of Ca2+/CaM-binding targets was increased in conjunction with elevated Ca2+. This indicates that CaM is not normally limiting for the activation of Ca 2+/CaM-dependent enzymes in HEK293 cells but reveals that the ratio of CaM to CaM-dependent targets is a potential mechanism for changing CaM availability. Next we considered whether cellular compartmentalization may act to regulate concentrations of available Ca2+/CaM in hippocampal neurons. We discovered changes in diffusion parameters of CaM under elevated Ca2+ conditions in the soma, neurite and nucleus which suggest that either the composition of cytoplasm is different in these compartments and/or they are composed of unique families of CaM-binding proteins. Finally, we return to the HEK293 cell and for the first time directly show the intracellular binding of CaM and CaMKII, an important target for CaM critical for neuronal function and plasticity. Furthermore, we analyzed the complex binding stoichiometry of this molecular interaction in the basal, activated and autophosphorylated states of CaMKII and determined the impact of this binding on CaM availability in the cell. Overall these results demonstrate that regulation of CaM availability is a viable cellular mechanism for regulating the output of CaM-dependent processes and that this process is tuned to the specific functional needs of a particular cell type and subcellular compartment. ^

Fluorophotometry is a sensitive method/technique to evaluate epithelial permeability in murine experimental dry eye

Purpose. Fluorophotometry is a well validated method for assessing corneal permeability in human subjects. However, with the growing importance of basic science animal research in ophthalmology, fluorophotometry’s use in animals must be further evaluated. The purpose of this study was to evaluate corneal epithelial permeability following desiccating stress using the modified Fluorotron Master™. ^ Methods. Corneal permeability was evaluated prior to and after subjecting 6-8 week old C57BL/6 mice to experimental dry eye (EDE) for 2 and 5 days (n=9/time point). Untreated mice served as controls. Ten microliters of 0.001% sodium fluorescein (NaF) were instilled topically into each mouse’s left eye to create an eye bath, and left to permeate for 3 minutes. The eye bath was followed by a generous wash with Buffered Saline Solution (BSS) and alignment with the Fluorotron Master™. Seven corneal scans using the Fluorotron Master were performed during 15 minutes (1 st post-wash scans), followed by a second wash using BSS and another set of five corneal scans (2nd post-wash scans) during the next 15 minutes. Corneal permeability was calculated using data calculated with the FM™ Mouse software. ^ Results. When comparing the difference between the Post wash #1 scans within the group and the Post wash #2 scans within the group using a repeated measurement design, there was a statistical difference in the corneal fluorescein permeability of the Post-wash #1 scans after 5 days (1160.21±108.26 vs. 1000.47±75.56 ng/mL, P<0.016 for UT-5 day comparison 8 [0.008]), but not after only 2 days of EDE compared to Untreated mice (1115.64±118.94 vs. 1000.47±75.56 ng/mL, P>0.016 for UT-2 day comparison [0.050]). There was no statistical difference between the 2 day and 5 day Post wash #1 scans (P=.299). The Post-wash #2 scans demonstrated that EDE caused a significant NaF retention at both 2 and 5 days of EDE compared to baseline, untreated controls (1017.92±116.25, 1015.40±120.68 vs. 528.22±127.85 ng/mL, P<0.05 [0.0001 for both]). There was no statistical difference between the 2 day and 5 day Post wash #2 scans (P=.503). The comparison between the Untreated post wash #1 with untreated post wash #2 scans using a Paired T-test showed a significant difference between the two sets of scans (P=0.000). There is also a significant difference between the 2 day comparison and the 5 day comparison (P values = 0.010 and 0.002, respectively). ^ Conclusion. Desiccating stress increases permeability of the corneal epithelium to NaF, and increases NaF retention in the corneal stroma. The Fluorotron Master is a useful and sensitive tool to evaluate corneal permeability in murine dry eye, and will be a useful tool to evaluate the effectiveness of dry eye treatments in animal-model drug trials.^

Statistical Methods for Differential Expressions of Genes Detected in Multiple-Condition Experiment of Microarray

Most studies of differential gene-expressions have been conducted between two given conditions. The two-condition experimental (TCE) approach is simple in that all genes detected display a common differential expression pattern responsive to a common two-condition difference. Therefore, the genes that are differentially expressed under the other conditions other than the given two conditions are undetectable with the TCE approach. In order to address the problem, we propose a new approach called multiple-condition experiment (MCE) without replication and develop corresponding statistical methods including inference of pairs of conditions for genes, new t-statistics, and a generalized multiple-testing method for any multiple-testing procedure via a control parameter C. We applied these statistical methods to analyze our real MCE data from breast cancer cell lines and found that 85 percent of gene-expression variations were caused by genotypic effects and genotype-ANAX1 overexpression interactions, which agrees well with our expected results. We also applied our methods to the adenoma dataset of Notterman et al. and identified 93 differentially expressed genes that could not be found in TCE. The MCE approach is a conceptual breakthrough in many aspects: (a) many conditions of interests can be conducted simultaneously; (b) study of association between differential expressions of genes and conditions becomes easy; (c) it can provide more precise information for molecular classification and diagnosis of tumors; (d) it can save lot of experimental resources and time for investigators.^