Coupling biomechanics to a cellular level model: An approach to patient-specific image driven multi-scale and multi-physics tumor simulation

Modeling of tumor growth has been performed according to various approaches addressing different biocomplexity levels and spatiotemporal scales. Mathematical treatments range from partial differential equation based diffusion models to rule-based cellular level simulators, aiming at both improving our quantitative understanding of the underlying biological processes and, in the mid- and long term, constructing reliable multi-scale predictive platforms to support patient-individualized treatment planning and optimization. The aim of this paper is to establish a multi-scale and multi-physics approach to tumor modeling taking into account both the cellular and the macroscopic mechanical level. Therefore, an already developed biomodel of clinical tumor growth and response to treatment is self-consistently coupled with a biomechanical model. Results are presented for the free growth case of the imageable component of an initially point-like glioblastoma multiforme tumor. The composite model leads to significant tumor shape corrections that are achieved through the utilization of environmental pressure information and the application of biomechanical principles. Using the ratio of smallest to largest moment of inertia of the tumor material to quantify the effect of our coupled approach, we have found a tumor shape correction of 20\% by coupling biomechanics to the cellular simulator as compared to a cellular simulation without preferred growth directions. We conclude that the integration of the two models provides additional morphological insight into realistic tumor growth behavior. Therefore, it might be used for the development of an advanced oncosimulator focusing on tumor types for which morphology plays an important role in surgical and/or radio-therapeutic treatment planning.

Theoretical considerations to the verification of dynamic multileaf collimated fields with a SLIC-type portal image detector

The verification possibilities of dynamically collimated treatment beams with a scanning liquid ionization chamber electronic portal image device (SLIC-EPID) are investigated. The ion concentration in the liquid of a SLIC-EPID and therefore the read-out signal is determined by two parameters of a differential equation describing the creation and recombination of the ions. Due to the form of this equation, the portal image detector describes a nonlinear dynamic system with memory. In this work, the parameters of the differential equation were experimentally determined for the particular chamber in use and for an incident open 6 MV photon beam. The mathematical description of the ion concentration was then used to predict portal images of intensity-modulated photon beams produced by a dynamic delivery technique, the sliding window approach. Due to the nature of the differential equation, a mathematical condition for 'reliable leaf motion verification' in the sliding window technique can be formulated. It is shown that the time constants for both formation and decay of the equilibrium concentration in the chamber is in the order of seconds. In order to guarantee reliable leaf motion verification, these time constants impose a constraint on the rapidity of the image-read out for a given maximum leaf speed. For a leaf speed of 2 cm s(-1), a minimum image acquisition frequency of about 2 Hz is required. Current SLIC-EPID systems are usually too slow since they need about a second to acquire a portal image. However, if the condition is fulfilled, the memory property of the system can be used to reconstruct the leaf motion. It is shown that a simple edge detecting algorithm can be employed to determine the leaf positions. The method is also very robust against image noise.

Comparative theoretical analysis of continuous wave laser cutting of metals at 1 and 10 um wavelength

We present a derivation and, based on it, an extension of a model originally proposed by V.G. Niziev to describe continuous wave laser cutting of metals. Starting from a local energy balance and by incorporating heat removal through heat conduction to the bulk material, we find a differential equation for the cutting profile. This equation is solved numerically and yields, besides the cutting profiles, the maximum cutting speed, the absorptivity profiles, and other relevant quantities. Our main goal is to demonstrate the model’s capability to explain some of the experimentally observed differences between laser cutting at around 1 and 10 μm wavelengths. To compare our numerical results to experimental observations, we perform simulations for exactly the same material and laser beam parameters as those used in a recent comparative experimental study. Generally, we find good agreement between theoretical and experimental results and show that the main differences between laser cutting with 1- and 10-μm beams arise from the different absorptivity profiles and absorbed intensities. Especially the latter suggests that the energy transfer, and thus the laser cutting process, is more efficient in the case of laser cutting with 1-μm beams.

Optimized CUDA-Based PDE Solver for Reaction Diffusion Systems on Arbitrary Surfaces

Partial differential equation (PDE) solvers are commonly employed to study and characterize the parameter space for reaction-diffusion (RD) systems while investigating biological pattern formation. Increasingly, biologists wish to perform such studies with arbitrary surfaces representing ‘real’ 3D geometries for better insights. In this paper, we present a highly optimized CUDA-based solver for RD equations on triangulated meshes in 3D. We demonstrate our solver using a chemotactic model that can be used to study snakeskin pigmentation, for example. We employ a finite element based approach to perform explicit Euler time integrations. We compare our approach to a naive GPU implementation and provide an in-depth performance analysis, demonstrating the significant speedup afforded by our optimizations. The optimization strategies that we exploit could be generalized to other mesh based processing applications with PDE simulations.

Differential effects of dexamethasone on the chondrogenesis of mesenchymal stromal cells: influence of microenvironment, tissue origin and growth factor

Mesenchymal stromal cells (MSCs), which reside within various tissues, are utilized in the engineering of cartilage tissue. Dexamethasone (DEX)--a synthetic glucocorticoid--is almost invariably applied to potentiate the growth-factor-induced chondrogenesis of MSCs in vitro, albeit that this effect has been experimentally demonstrated only for transforming-growth-factor-beta (TGF-β)-stimulated bone-marrow-derived MSCs. Clinically, systemic glucocorticoid therapy is associated with untoward side effects (e.g., bone loss and increased susceptibility to infection). Hence, the use of these agents should be avoided or limited. We hypothesize that the influence of DEX on the chondrogenesis of MSCs depends upon their tissue origin and microenvironment [absence or presence of an extracellular matrix (ECM)], as well as upon the nature of the growth factor. We investigated its effects upon the TGF-β1- and bone-morphogenetic-protein 2 (BMP-2)-induced chondrogenesis of MSCs as a function of tissue source (bone marrow vs. synovium) and microenvironment [cell aggregates (no ECM) vs. explants (presence of a natural ECM)]. In aggregates of bone-marrow-derived MSCs, DEX enhanced TGF-β1-induced chondrogenesis by an up-regulation of cartilaginous genes, but had little influence on the BMP-2-induced response. In aggregates of synovial MSCs, DEX exerted no remarkable effect on either TGF-β1- or BMP-2-induced chondrogenesis. In synovial explants, DEX inhibited BMP-2-induced chondrogenesis almost completely, but had little impact on the TGF-β1-induced response. Our data reveal that steroids are not indispensable for the chondrogenesis of MSCs in vitro. Their influence is context dependent (tissue source of the MSCs, their microenvironment and the nature of the growth-factor). This finding has important implications for MSC based approaches to cartilage repair.

Social cognition as a mediator variable between neurocognition and functional outcome in schizophrenia: empirical review and new results by structural equation modeling

Cognitive impairments are currently regarded as important determinants of functional domains and are promising treatment goals in schizophrenia. Nevertheless, the exact nature of the interdependent relationship between neurocognition and social cognition as well as the relative contribution of each of these factors to adequate functioning remains unclear. The purpose of this article is to systematically review the findings and methodology of studies that have investigated social cognition as a mediator variable between neurocognitive performance and functional outcome in schizophrenia. Moreover, we carried out a study to evaluate this mediation hypothesis by the means of structural equation modeling in a large sample of 148 schizophrenia patients. The review comprised 15 studies. All but one study provided evidence for the mediating role of social cognition both in cross-sectional and in longitudinal designs. Other variables like motivation and social competence additionally mediated the relationship between social cognition and functional outcome. The mean effect size of the indirect effect was 0.20. However, social cognitive domains were differentially effective mediators. On average, 25% of the variance in functional outcome could be explained in the mediation model. The results of our own statistical analysis are in line with these conclusions: Social cognition mediated a significant indirect relationship between neurocognition and functional outcome. These results suggest that research should focus on differential mediation pathways. Future studies should also consider the interaction with other prognostic factors, additional mediators, and moderators in order to increase the predictive power and to target those factors relevant for optimizing therapy effects.

Impact of salt on cardiac differential gene expression and coronary lesion in normotensive mineralocorticoid-treated mice

We previously reported that excess of deoxycorticosterone-acetate (DOCA)/salt-induced cardiac hypertrophy in the absence of hypertension in one-renin gene mice. This model allows us to study molecular mechanisms of high-salt intake in the development of cardiovascular remodeling, independently of blood pressure in a high mineralocorticoid state. In this study, we compared the effect of 5-wk low- and high-salt intake on cardiovascular remodeling and cardiac differential gene expression in mice receiving the same amount of DOCA. Differential gene and protein expression was measured by high-density cDNA microarray assays, real-time PCR and Western blot analysis in DOCA-high salt (HS) vs. DOCA-low salt (LS) mice. DOCA-HS mice developed cardiac hypertrophy, coronary perivascular fibrosis, and left ventricular dysfunction. Differential gene and protein expression demonstrated that high-salt intake upregulated a subset of genes encoding for proteins involved in inflammation and extracellular matrix remodeling (e.g., Col3a1, Col1a2, Hmox1, and Lcn2). A major subset of downregulated genes encoded for transcription factors, including myeloid differentiation primary response (MyD) genes. Our data provide some evidence that vascular remodeling, fibrosis, and inflammation are important consequences of a high-salt intake in DOCA mice. Our study suggests that among the different pathogenic factors of cardiac and vascular remodeling, such as hypertension and mineralocorticoid excess and sodium intake, the latter is critical for the development of the profibrotic and proinflammatory phenotype observed in the heart of normotensive DOCA-treated mice.

Differential gene and protein expression in abluminal sprouting and intraluminal splitting forms of angiogenesis

In adult skeletal muscle, abluminal sprouting or longitudinal splitting of capillaries can be initiated separately by muscle overload and elevated microcirculation shear stress respectively. In the present study, gene and protein expression patterns associated with the different forms of angiogenesis were examined using a targeted gene array (Superarray), validated by quantitative RT (reverse transcription)-PCR and immunoblots. Sprouting angiogenesis induced large changes in expression levels in genes associated with extracellular matrix remodelling, such as MMP-2 (matrix metalloproteinase-2), TIMP (tissue inhibitor of metalloproteinases), SPARC (secreted protein, acidic and rich in cysteine) and thrombospondin. Changes in neuropilin, midkine and restin levels, which may underpin changes in endothelial morphology, were seen during splitting angiogenesis. Up-regulation of VEGF (vascular endothelial growth factor), Flk-1, angiopoietin-2 and PECAM-1 (platelet/endothelial cell adhesion molecule-1) was seen in both forms of angiogenesis, representing a common angiogenic response of endothelial cells. In conclusion, the present study demonstrates that general angiogenic signals from growth factors can be influenced by the local microenvironment resulting in differing forms of capillary growth to produce a co-ordinated expansion of the vascular bed.

Differential regulation of metzincins in experimental chronic renal allograft rejection: potential markers and novel therapeutic targets

Chronic renal allograft rejection is characterized by alterations in the extracellular matrix compartment and in the proliferation of various cell types. These features are controlled, in part by the metzincin superfamily of metallo-endopeptidases, including matrix metalloproteinases (MMPs), a disintegrin and metalloproteinase (ADAM) and meprin. Therefore, we investigated the regulation of metzincins in the established Fisher to Lewis rat kidney transplant model. Studies were performed using frozen homogenates and paraffin sections of rat kidneys at day 0 (healthy controls) and during periods of chronic rejection at day +60 and day +100 following transplantation. The messenger RNA (mRNA) expression was examined by Affymetrix Rat Expression Array 230A GeneChip and by real-time Taqman polymerase chain reaction analyses. Protein expression was studied by zymography, Western blot analyses, and immunohistology. mRNA levels of MMPs (MMP-2/-11/-12/-14), of their inhibitors (tissue inhibitors of metalloproteinase (TIMP)-1/-2), ADAM-17 and transforming growth factor (TGF)-beta1 significantly increased during chronic renal allograft rejection. MMP-2 activity and immunohistological staining were augmented accordingly. The most important mRNA elevation was observed in the case of MMP-12. As expected, Western blot analyses also demonstrated increased production of MMP-12, MMP-14, and TIMP-2 (in the latter two cases as individual proteins and as complexes). In contrast, mRNA levels of MMP-9/-24 and meprin alpha/beta had decreased. Accordingly, MMP-9 protein levels and meprin alpha/beta synthesis and activity were downregulated significantly. Members of metzincin families (MMP, ADAM, and meprin) and of TIMPs are differentially regulated in chronic renal allograft rejection. Thus, an altered pattern of metzincins may represent novel diagnostic markers and possibly may provide novel targets for future therapeutic interventions.

Mechanical signals regulating extracellular matrix gene expression in fibroblasts

Mechanical forces are essential for connective tissue homeostasis. The extracellular matrix (ECM) plays a key role in the transmission of forces generated by the organism (e.g. muscle contraction) and externally applied (e.g. gravity). The expression of specific ECM proteins such as collagens and tenascin-C, as well as of matrix metalloproteinases, involved in their turnover, is influenced by mechanical stimuli. The precise mechanisms by which mechanical strains are translated into chemical signals and lead to differential gene expression are however not fully understood. Cell-matrix adhesion sites are good candidates for hosting a "mechanosensory switch", as they transmit forces from the ECM to the cytoskeleton and vice versa by physically linking the cytoskeleton to the ECM. Integrins, transmembrane proteins located to these adhesion sites, have been shown to trigger a set of internal signaling cascades after mechanical stimulation. We have shown that the expression level of tenascin-C directly correlates with externally applied mechanical stress, as well as with RhoA/RhoA-dependent kinase-mediated cytoskeletal tension. Presumably other genes are regulated in a similar manner. The changes in ECM composition and mechanical properties derived from mechanical stress are relevant in medical intervention after ligament and tendon injury.

Differential effect of p47phox and gp91phox deficiency on the course of Pneumococcal Meningitis.

Bacterial meningitis is a severe inflammatory disease of the central nervous system and is characterized by massive infiltration of granulocytes into the cerebrospinal fluid (CSF). To assess the role of NADPH oxidase-derived reactive oxygen species (ROS) in pneumococcal meningitis, mice deficient in either the gp91 subunit (essential for functioning of the phagocyte enzyme) or the p47 subunit (essential for functioning of homologous enzymes in nonphagocytic cells) were intracisternally infected with live Streptococcus pneumoniae, and defined disease parameters were measured during the acute stage of infection. While none of the parameters measured (including CSF bacterial titers) were significantly different in gp91(-/-) and wild-type mice, the infection in p47(-/-) mice was associated with significantly increased inflammation of the subarachnoid and ventricular space, disruption of the blood-brain barrier, and the presence of interleukin-1 beta, tumor necrosis factor alpha, and matrix metalloproteinase 9 in the cortex. These changes were associated with approximately 10-fold-higher CSF bacterial titers in p47(-/-) mice than in wild-type mice (P < 0.001). In contrast to infection with live bacteria, the inflammatory response, including CSF leukocytosis, was significantly attenuated in p47(-/-) mice (but not gp91(-/-) mice) challenged with a fixed number of heat-inactivated pneumococci. Impairment of the host defense appeared to be responsible for the higher bacterial titers in p47(-/-) mice. Therefore, these results indicate that ROS generated by a gp91-independent NADPH oxidase(s) are important for establishing an adequate inflammatory response to pneumococcal CSF infection.

Dichotomous Hamiltonians with unbounded entries and solutions of Riccati equations

An operator Riccati equation from systems theory is considered in the case that all entries of the associated Hamiltonian are unbounded. Using a certain dichotomy property of the Hamiltonian and its symmetry with respect to two different indefinite inner products, we prove the existence of nonnegative and nonpositive solutions of the Riccati equation. Moreover, conditions for the boundedness and uniqueness of these solutions are established.

Differential inflammatory response of dental pulp explants and fibroblasts to saliva

Abstract AIM: To investigate the inflammatory response of dental pulp fibroblasts and the respective explants to whole saliva. METHODOLOGY: Explants from human and porcine dental pulp tissue and isolated dental pulp fibroblasts were used to investigate the inflammatory response to sterile saliva. Cytokine and chemokine expression was assessed by RT-PCR. Western blot analysis and pharmacologic inhibitors were used to determine the involvement of signalling pathways. RESULTS: Dental pulp explants of human and porcine origin exposed to human saliva exhibited no major changes of IL-6 and IL-8 mRNA expression (P > 0.05). In contrast, isolated porcine and human dental pulp fibroblasts, when stimulated with human saliva, exhibited a vastly increased expression of IL-6 and IL-8 mRNA (P < 0.05). In pulp fibroblasts, saliva also increased the expression of other cytokines and chemokines via activation of NFkappaB, ERK and p38 signalling. Notably, a significantly reduced inflammatory response was elicited when pulp fibroblasts were transiently exposed to saliva. CONCLUSIONS: Saliva has a potential impact on inflammation of dental pulp fibroblasts in vitro but not when cells are embedded in the intrinsic extracellular matrix of the explant tissue.