Multi-Scale Modeling for Image Analysis of Brain Tumor Studies

Image-based modeling of tumor growth combines methods from cancer simulation and medical imaging. In this context, we present a novel approach to adapt a healthy brain atlas to MR images of tumor patients. In order to establish correspondence between a healthy atlas and a pathologic patient image, tumor growth modeling in combination with registration algorithms is employed. In a first step, the tumor is grown in the atlas based on a new multi-scale, multi-physics model including growth simulation from the cellular level up to the biomechanical level, accounting for cell proliferation and tissue deformations. Large-scale deformations are handled with an Eulerian approach for finite element computations, which can operate directly on the image voxel mesh. Subsequently, dense correspondence between the modified atlas and patient image is established using nonrigid registration. The method offers opportunities in atlasbased segmentation of tumor-bearing brain images as well as for improved patient-specific simulation and prognosis of tumor progression.

Pseudodominant Inheritance of Goitrous Congenital Hypothyroidism Caused by TPO Mutations: Molecular and in Silico Studies

Context and Objective: Most cases of goitrous congenital hypothyroidism (CH) from thyroid dyshormonogenesis 1) follow a recessive mode of inheritance and 2) are due to mutations in the thyroid peroxidase gene (TPO). We report the genetic mechanism underlying the apparently dominant inheritance of goitrous CH in a nonconsanguineous family of French Canadian origin. Design, Setting, and Participants: Two brothers identified by newborn TSH screening had severe hypothyroidism and a goiter with increased (99m)Tc uptake. The mother was euthyroid, but the father and two paternal uncles had also been diagnosed with goitrous CH. After having excluded PAX8 gene mutations, we hypothesized that the underlying defect could be TPO mutations. Results: Both compound heterozygous siblings had inherited a mutant TPO allele carried by their mother (c.1496delC; p.Pro499Argfs2X), and from their father, one brother had inherited a missense mutation (c.1978C-->G; p.Gln660Glu) and the other an insertion (c.1955insT; p.Phe653Valfs15X). The thyroid gland of one uncle who is a compound heterozygote for TPO mutations (p.Phe653Valfs15X/p.Gln660Glu) was removed because of concurrent multiple endocrine neoplasia type 2A. Immunohistochemistry revealed normal TPO staining, implying that Gln660Glu TPO is expressed properly. Modeling of this mutant in silico suggests that its three-dimensional structure is conserved, whereas the electrostatic binding energy between the Gln660Glu TPO and its heme group becomes repulsive. Conclusion: We report a pedigree presenting with pseudodominant goitrous CH due to segregation of three different TPO mutations. Although goitrous CH generally follows a recessive mode of inheritance, the high frequency of TPO mutations carriers may lead to pseudodominant inheritance.

Modeling of human P450 oxidoreductase structure by in silico mutagenesis and MD simulation

P450 oxidoreductase (POR) is the obligate electron donor for microsomal cytochrome P450s and mutations in POR cause several metabolic disorders. We have modeled the structure of human P450 oxidoreductase by in silico amino acid replacements in the rat POR crystal structure. The rat POR has 94% homology with human POR and 38 amino acids were replaced to make its sequence identical to human POR. Several rounds of molecular dynamic simulations refined the model and removed structural clashes from side chain alterations of replaced amino acids. This approach has the advantage of keeping the cofactor contacts and structural features of the core enzyme intact which could not be achieved by homology based approaches. The final model from our approach was of high quality and compared well with experimentally determined structures of other PORs. This model will be used for analyzing the structural implications of mutations and polymorphisms in human POR.

Determination of the molecular subtypes of diffuse large B-cell lymphomas using immunohistochemistry: a case series from the Inselspital, Bern, and a critical appraisal of this determination in Switzerland

The two major subtypes of diffuse large B-cell lymphoma (DLBCL) (germinal centre B-cell - like (GCB-DLBCL) and activated B-cell - like (ABC-DLBCL)) are defined by means of gene expression profiling (GEP). Patients with GCB-DLBCL survive longer with the current standard regimen R-CHOP than patients with ABC-DLBCL. As GEP is not part of the current routine diagnostic work-up, efforts have been made to find a substitute than involves immunohistochemistry (IHC). Various algorithms achieved this with 80-90% accuracy. However, conflicting results on the appropriateness of IHC have been reported. Because it is likely that the molecular subtypes will play a role in future clinical practice, we assessed the determination of the molecular DLBCL subtypes by means of IHC at our University Hospital, and some aspects of this determination elsewhere in Switzerland. The most frequently used Hans algorithm includes three antibodies (against CD10, bcl-6 and MUM1). From records of the routine diagnostic work-up, we identified 51 of 172 (29.7%) newly diagnosed and treated DLBCL cases from 2005 until 2010 with an assigned DLBCL subtype. DLBCL subtype information was expanded by means of tissue microarray analysis. The outcome for patients with the GCB subtype was significantly better compared with those with the non-GC subtype, independent of the age-adjusted International Prognostic Index. We found a lack of standardisation in the subtype determination by means of IHC in Switzerland and significant problems of reproducibility. We conclude that the Hans algorithm performs well in our hands and that awareness of this important matter is increasing. However, outside clinical trials, vigorous efforts to standardise IHC determination are needed as DLBCL subtype-specific therapies emerge.

Adaptive Algorithms for Personalized Diabetes Treatment

Dynamic systems, especially in real-life applications, are often determined by inter-/intra-variability, uncertainties and time-varying components. Physiological systems are probably the most representative example in which population variability, vital signal measurement noise and uncertain dynamics render their explicit representation and optimization a rather difficult task. Systems characterized by such challenges often require the use of adaptive algorithmic solutions able to perform an iterative structural and/or parametrical update process towards optimized behavior. Adaptive optimization presents the advantages of (i) individualization through learning of basic system characteristics, (ii) ability to follow time-varying dynamics and (iii) low computational cost. In this chapter, the use of online adaptive algorithms is investigated in two basic research areas related to diabetes management: (i) real-time glucose regulation and (ii) real-time prediction of hypo-/hyperglycemia. The applicability of these methods is illustrated through the design and development of an adaptive glucose control algorithm based on reinforcement learning and optimal control and an adaptive, personalized early-warning system for the recognition and alarm generation against hypo- and hyperglycemic events.

Simulating Cortical Development as a Self Constructing Process: A Novel Multi-Scale Approach Combining Molecular and Physical Aspects

Current models of embryological development focus on intracellular processes such as gene expression and protein networks, rather than on the complex relationship between subcellular processes and the collective cellular organization these processes support. We have explored this collective behavior in the context of neocortical development, by modeling the expansion of a small number of progenitor cells into a laminated cortex with layer and cell type specific projections. The developmental process is steered by a formal language analogous to genomic instructions, and takes place in a physically realistic three-dimensional environment. A common genome inserted into individual cells control their individual behaviors, and thereby gives rise to collective developmental sequences in a biologically plausible manner. The simulation begins with a single progenitor cell containing the artificial genome. This progenitor then gives rise through a lineage of offspring to distinct populations of neuronal precursors that migrate to form the cortical laminae. The precursors differentiate by extending dendrites and axons, which reproduce the experimentally determined branching patterns of a number of different neuronal cell types observed in the cat visual cortex. This result is the first comprehensive demonstration of the principles of self-construction whereby the cortical architecture develops. In addition, our model makes several testable predictions concerning cell migration and branching mechanisms.

Resolving diffusion in clay minerals at different time scales: Combination of experimental and modeling approaches

Since no single experimental or modeling technique provides data that allow a description of transport processes in clays and clay minerals at all relevant scales, several complementary approaches have to be combined to understand and explain the interplay between transport relevant phenomena. In this paper molecular dynamics simulations (MD) were used to investigate the mobility of water in the interlayer of montmorillonite (Mt), and to estimate the influence of mineral surfaces and interlayer ions on the water diffusion. Random Walk (RW) simulations based on a simplified representation of pore space in Mt were used to estimate and understand the effect of the arrangement of Mt particles on the meso- to macroscopic diffusivity of water. These theoretical calculations were complemented with quasielastic neutron scattering (QENS) measurements of aqueous diffusion in Mt with two pseudo-layers of water performed at four significantly different energy resolutions (i.e. observation times). The size of the interlayer and the size of Mt particles are two characteristic dimensions which determine the time dependent behavior of water diffusion in Mt. MD simulations show that at very short time scales water dynamics has the characteristic features of an oscillatory motion in the cage formed by neighbors in the first coordination shell. At longer time scales, the interaction of water with the surface determines the water dynamics, and the effect of confinement on the overall water mobility within the interlayer becomes evident. At time scales corresponding to an average water displacement equivalent to the average size of Mt particles, the effects of tortuosity are observed in the meso- to macroscopic pore scale simulations. Consistent with the picture obtained in the simulations, the QENS data can be described using a (local) 3D diffusion at short observation times, whereas at sufficiently long observation times a 2D diffusive motion is clearly observed. The effects of tortuosity measured in macroscopic tracer diffusion experiments are in qualitative agreement with RW simulations. By using experimental data to calibrate molecular and mesoscopic theoretical models, a consistent description of water mobility in clay minerals from the molecular to the macroscopic scale can be achieved. In turn, simulations help in choosing optimal conditions for the experimental measurements and the data interpretation. (C) 2014 Elsevier B.V. All rights reserved.

Modeling immune functions of the mouse blood-cerebrospinal fluid barrier in vitro: primary rather than immortalized mouse choroid plexus epithelial cells are suited to study immune cell migration across this brain barrier.

BACKGROUND The blood-cerebrospinal fluid barrier (BCSFB) established by the choroid plexus (CP) epithelium has been recognized as a potential entry site of immune cells into the central nervous system during immunosurveillance and neuroinflammation. The location of the choroid plexus impedes in vivo analysis of immune cell trafficking across the BCSFB. Thus, research on cellular and molecular mechanisms of immune cell migration across the BCSFB is largely limited to in vitro models. In addition to forming contact-inhibited epithelial monolayers that express adhesion molecules, the optimal in vitro model must establish a tight permeability barrier as this influences immune cell diapedesis. METHODS We compared cell line models of the mouse BCSFB derived from the Immortomouse(®) and the ECPC4 line to primary mouse choroid plexus epithelial cell (pmCPEC) cultures for their ability to establish differentiated and tight in vitro models of the BCSFB. RESULTS We found that inducible cell line models established from the Immortomouse(®) or the ECPC4 tumor cell line did not express characteristic epithelial proteins such as cytokeratin and E-cadherin and failed to reproducibly establish contact-inhibited epithelial monolayers that formed a tight permeability barrier. In contrast, cultures of highly-purified pmCPECs expressed cytokeratin and displayed mature BCSFB characteristic junctional complexes as visualized by the junctional localization of E-cadherin, β-catenin and claudins-1, -2, -3 and -11. pmCPECs formed a tight barrier with low permeability and high electrical resistance. When grown in inverted filter cultures, pmCPECs were suitable to study T cell migration from the basolateral to the apical side of the BCSFB, thus correctly modelling in vivo migration of immune cells from the blood to the CSF. CONCLUSIONS Our study excludes inducible and tumor cell line mouse models as suitable to study immune functions of the BCSFB in vitro. Rather, we introduce here an in vitro inverted filter model of the primary mouse BCSFB suited to study the cellular and molecular mechanisms mediating immune cell migration across the BCSFB during immunosurveillance and neuroinflammation.